Genetics of the phenotypic evolution in sheep: a molecular look at diversity-driving genes
Genetics Selection Evolution volume 54, Article number: 61 (2022)
After domestication, the evolution of phenotypically-varied sheep breeds has generated rich biodiversity. This wide phenotypic variation arises as a result of hidden genomic changes that range from a single nucleotide to several thousands of nucleotides. Thus, it is of interest and significance to reveal and understand the genomic changes underlying the phenotypic variation of sheep breeds in order to drive selection towards economically important traits.
Various traits contribute to the emergence of variation in sheep phenotypic characteristics, including coat color, horns, tail, wool, ears, udder, vertebrae, among others. The genes that determine most of these phenotypic traits have been investigated, which has generated knowledge regarding the genetic determinism of several agriculturally-relevant traits in sheep. In this review, we discuss the genomic knowledge that has emerged in the past few decades regarding the phenotypic traits in sheep, and our ultimate aim is to encourage its practical application in sheep breeding. In addition, in order to expand the current understanding of the sheep genome, we shed light on research gaps that require further investigation.
Although significant research efforts have been conducted in the past few decades, several aspects of the sheep genome remain unexplored. For the full utilization of the current knowledge of the sheep genome, a wide practical application is still required in order to boost sheep productive performance and contribute to the generation of improved sheep breeds. The accumulated knowledge on the sheep genome will help advance and strengthen sheep breeding programs to face future challenges in the sector, such as climate change, global human population growth, and the increasing demand for products of animal origin.
One of the most intriguing questions in the biological sciences is how a genome can result in wonderfully varied phenotypes? Throughout and following sheep domestication, various genotypes have been generated. Certain genomic patterns have been selected naturally or artificially to give rise to various phenotypes [1,2,3], which have culminated in the generation of distinct currently known sheep breeds (Fig. 1). Among the main phenotypic traits that vary widely in sheep, are included tails, horns, and coat colors. For instance, large-fat tails in sheep appeared as an adaptive response to harsh environmental conditions and food scarcity . However, the thin-tail phenotype has currently become more desirable, considering that the fat-tail phenotype carries several negative implications, which likely influence animal locomotion, mating, food conversion, breeding costs, among others. Sheep horns are commonly used for protection against predators and as sexual weaponry. In modern breeding practices, the hornless (or polled) phenotype has become increasingly desirable as a way to avoid the need for dehorning , which is usually performed to protect animals and producers from accidental injuries. In addition, coat color serves in the adaptation to environmental conditions and avoiding predators. Currently, selection for specific sheep coat color patterns is performed based on the production system, e.g., white color is desirable in several wool-producing sheep breeds . Combined with other phenotypic traits, coat color is considered a determinant of sheep breeds. Several phenotypic traits are significantly associated with production and animal welfare parameters. Hence, understanding the genetic background of these phenotypic traits is of great importance for a rapid and direct selection and improvement of sheep breeds. The sequencing of the sheep genome and the elucidation of its content and variations have offered an unprecedented opportunity for a better understanding of the underlying genetic traits controlling sheep phenotype [7,8,9,10].
Recent advances in sheep genomics have helped in identifying genes that are implicated in determining major phenotypic traits in sheep. For instance, it has been reported that the PDGFD, BMP2, HOXB13, and TBXT genes regulate sheep tail phenotype [11,12,13,14], RXFP2 and HOXD1 regulate horn and multi-horn traits [5, 15], whereas ASIP, MC1R, and TYRP1 regulate coat color patterns [16,17,18]. However, the causative genomic variations within or around some of these genes are currently unknown. Given the importance of the genes controlling sheep phenotypic traits that are linked with adaption, production, and welfare traits, in this review, we collectively discuss the genomic history of these traits and describe genes and variants that are currently known for these traits and highlight traits which need further investigation to elucidate their genetic determinism. We discuss potential genes and their variants that control main phenotypic traits in sheep, including tail, horns, coat color, bone color, fat color, wool length, fleece type, ear size, number of nipples, number of thoracic and lumbar vertebrae, body height, and stature. We anticipate that exploring genes that control significant phenotypic traits in sheep will foment selective breeding programs and genetic-based strategies for the rapid improvement of productivity and welfare of sheep populations.
The long-standing conundrum of sheep tails: a genomic perspective
Sheep tail phenotypic patterns vary among different sheep breeds in terms of length, directionality, and fat accumulation level  (Fig. 2). The genetic basis of sheep tail phenotypic variations has been elucidated based on genomic studies. Moradi et al.  first applied a genomic scan approach to identify genomic loci that underlie the phenotypic differences between fat- and thin-tailed Iranian sheep breeds, and revealed seven genomic regions that overlap with quantitative trait loci (QTL) shown to affect fat and carcass yield. In addition, three genomic regions located on Ovis aries chromosome 5 (OAR5), OAR7, and OARX were suggested to be associated with the sheep tail phenotype . A subsequent genome-wide analysis by Wei et al.  using ten representatives of Chinese sheep breeds revealed the association of PPP1CC and PDGFD with significant differences between fat- and thin-tailed animals. Specifically, PDGFD was positively selected in fat-tailed sheep breeds. In another study, BMP2 and VRTN were identified as the genes that are most probably involved in the determination of sheep tail phenotype .
Subsequently, Yuan et al.  proposed that the HOXA11, BMP2, PPP1CC, SP3, SP9, WDR92, PROKR1, and ETAA1 genes potentially play an important role in determining the fat-tailed phenotype. An integrated analysis of genome-wide association study (GWAS) and X-chromosome-wide association analysis (XWAS) in small-tailed and large-tailed Han sheep revealed a set of genes (namely, CREB1, STEAP4, CTBP1, and NRIP1) that could be potentially associated with lipid metabolism or adipocyte development . A candidate genomic region, i.e., OARX 88–89 Mb, was also identified, which was shown to contain significant single nucleotide polymorphisms (SNPs) that could be associated with tail phenotype . In addition, by analyzing genotypic data of 11 Ethiopian indigenous sheep breeds and two populations of thin-tailed sheep from Sudan, a set of potential genes associated with sheep tail phenotype was proposed . This includes putative candidate regions that spanned genes affecting growth traits and fat deposition (NPR2, HINT2, SPAG8, and INSR), tail formation (ALX4, HOXB13, and BMP4), embryonic development of tendons, bones, and cartilages (EYA2 and SULF2), lipogenesis and intracellular long-chain fatty acids transport (FABP3) . In another study, genomic data of three Iranian sheep breeds (i.e., Baluchi, Lori-Bakhtiari, and Zel) were also analyzed, which revealed that the RPS6KA3, HOXB13, HOXB9, ACAN, and CHRDL1 genes as potentially influencing the skeletal system and tail traits when comparing fat-tailed Baluchi sheep and thin-tailed Zel sheep; conversely, NPR2 and EN1 were suggested as potentially influencing energy metabolism, skeletal system, and tail traits when comparing fat-tailed Lori-Bakhtiari sheep and thin-tailed Zel sheep .
Pan et al.  analyzed the whole-genome sequences of 89 sheep and identified two significant signals on OAR13 and OAR15. The sweep signal on OAR13 was consistent with previous studies [11, 12]. In particular, Moioli et al.  suggested that the BMP2 gene could be potentially causative for the sheep fat tail formation, whereas Wei et al.  suggested the PPP1CC gene, which is a paralogue of the original PPP1CC gene on OAR13 (NCBI gene identifier: LOC101117953). The other signal on OAR15 contained the PDGFD gene, which was previously highlighted by Wei et al. . The newly characterized LOC101117953 has been suggested to originate from a retro-copy (r-PPP1CC) of the PPP1CC gene via clear intron losses . Since r-PPP1CC is a promoter-less non-functional pseudo-gene, it was assumed to be less likely the causative gene. However, the locus can potentially serve as a regulatory element, which changes the expressional patterns of the distant BMP2 gene . Mastrangelo et al.  compared the genomic data of two fat- and thin-tailed sheep breeds, and again revealed previously reported genes, i.e., VRTN and BMP2, with a potential role in determining variations in the number of vertebrae and tail formation, respectively, and PDGFRA and PDGFD, which seem to be involved in preadipocyte differentiation. Moreover, Mastrangelo et al.  also reported a set of genes, namely, ANAPC1, MSRB3, CXXC5, PSD2, SLIT2, EPHA5, CHP1, OIP5, PCDH9, CDS2, BMP2, and OAS2, which could be associated with fat deposition in sheep breeds from Africa and Eurasia.
Recently, a number of large-scale genomic studies have been performed to elucidate the phenotypic differences in sheep tail formation. Li et al.  investigated the genetic background of the configuration of sheep tails by performing separate pairwise-population selection tests and analyzing genotype patterns located in the promoter region; in this study, the PDGFD gene was consistently selected by multiple comparisons and its expression varied between fat- and thin-tailed sheep breeds. Dong et al.  highlighted a 6.8-kb region within the first intron of the PDGFD gene that contained the most differentiated SNPs between fat- and thin-tailed sheep breeds and that likely harbors regulatory mutation(s) determining the fat-tail configuration. Furthermore, it has been demonstrated that expression of PDGFD was high on days 60 and 70 of embryonic development in the tail adipose tissue, which suggests a potential role of PDGFD in the initiation and/or commitment of preadipocytes; in contrast, expression of PDGFD decreased between days 80 and 90 of embryonic development when fat started to accumulate in the tail, thus suggesting a negative association with fat maturation .
Baazaoui et al.  employed a combination of runs of homozygosity (ROH) hotpots analysis and an FST-outlier approach to identify genomic regions under selection in Tunisian sheep breeds. An outstanding selection signal harbored several relevant genes linked with lipid metabolic processes (i.e., CDS2, PROKR1, and BMP2). The study also emphasized the potential involvement of BMP2 in the sheep tail formation . Genes linked with the skeletal tail phenotype, including HOXB13, NR5A1, NR6A1, KDM6B, and NR3C1, have also been proposed to play a major role in the adaptation of African sheep . In addition, allele-specific expression (ASE) in thin-tailed Texel sheep and fat-rumped Kazakh sheep was investigated . In this study, 2405 ASE genes were described from adipose tissue; among these, PDGFD and LYPLA1 were highlighted as possibly linked with lipid-related biological processes. In addition, seven ASE SNPs were found in the PDGFD gene, and a distinct genotype pattern located in the upstream region of PDGFD was detected by selective sweep analysis between mouflon, fat-tailed, and thin-tailed sheep. Moreover, potential exon splicing events within the PDGFD gene in thin-tailed sheep have been described . More recently, Li et al.  emphasized the significance of both PDGFD and BMP2 signals in sheep tail formation, and described a new variant based on a 169-bp insertion close to the 5′-UTR region of HOXB13, which is likely to be associated with the long tail phenotype. This variant was further confirmed as a 167-bp segment with a linked missense mutation (c.23C>G) in Merinolandschaf sheep . It is likely that multiple genes may contribute to the configuration of sheep tails since tail phenotype and related fat deposition levels vary greatly between different sheep breeds from different geographical regions.
Although the level of fat deposition mainly determines sheep tail phenotype, the length and number of caudal vertebrae have also been considered as important factors that influence tail morphology (Fig. 3). Zhi et al.  analyzed genomic differences between two lines of Chinese Hulunbuir sheep of short tail and normal tail, which revealed the TBXT gene (also known as T; Brachyury) as a candidate for determining the number of caudal vertebrae in sheep. Subsequently, Han et al.  confirmed the previously reported missense mutation (c.334G>T) with a completely linked synonymous variant (c.333G>C) to be associated with the characteristic short caudal vertebrae in fat-rumped sheep breeds. Interestingly, the functional role of the TBXT gene in forming the sheep tail length was validated using the CRISPR/Cas9 system . Collectively, the abovementioned findings provide significant information obtained in research studies from the past few decades to understand the genetic causes of phenotypic variations in the configuration of sheep tails. These genomic studies identified sets of genes that could potentially be associated with phenotypic variations in sheep tails (see Additional file 1: Fig. S1a), and narrowed down the candidate genes and genomic regions involved in sheep tail formation (e.g., BMP2, LOC101117953, PDGFD, HOXB13, and TBXT). Further genomic investigations of sheep tails are shown in Additional file 2: Table S1. Moreover, sheep tail formation has been investigated using transcriptomic [mRNAs (see Additional file 2: Table S2 and Additional file 1: Fig. S1b), miRNAs (see Additional file 2: Table S3), and lncRNAs (see Additional file 2: Table S4)] and proteomic (see Additional file 2: Table S5) approaches, as well as direct gene expression and potential SNP associations (see Additional file 2: Table S6). These studies further identified potential genes, regulatory factors, and molecular pathways that might contribute to the formation and phenotype variation of sheep tails.
Development of sheep horns: the remainders of their wild past
The RXFP2 gene and its role in sheep horn development
Horns (cranial appendages or headgear)  are typically used by male animals in mating conflicts for dominance . Sheep can be horned, polled (hornless), or scurred (i.e., vestigial and deformed horns) (Fig. 4). Hornedness and polledness phenotypes have varying patterns among different sheep breeds. In certain cases, horns are found in both sexes, but their size is larger in males; males can be horned but females polled; and both sexes can be polled . In other cases, horned and polled male and female representatives can be found within the same breed. Early attempts to understand the genetics of sheep horn development date from over a century ago [37,38,39,40,41]. Since then, several efforts have been made to understand the mode of inheritance of this unique phenotypic trait in sheep [42,43,44,45,46,47,48]. With the emergence of genomic-based approaches, the underlying genetic basis of this phenotypic trait in sheep has been progressively revealed. In 1996, Montgomery et al.  refuted the existence of a genetic locus determining horn traits in sheep on OAR1 as seen in cattle. However, a horn-determining locus was mapped on OAR10 in sheep that originate from Merino and Romney crosses. Subsequently, using a free-living population of Soay sheep (Fig. 5) in St. Kilda, Beraldi et al.  also mapped the horn locus on OAR10 at a location similar to that previously described in modern domestic sheep. In addition, Pickering et al.  localized the horn locus within a 14-Mb region using six genetic markers. Subsequently, by using fine-mapping with additional genetic markers, the horn locus was localized within a 50-kb region. This study also reported a large (approximately 3 kb) retrotransposed insertion in the non-coding 3′-UTR region of a then unnamed candidate gene, which was found exclusively in polled individuals .
Johnston et al.  further narrowed the horn locus down to a ~ 7.4-cM interval on OAR10 in a Soay sheep population. Subsequently, Johnston et al.  performed a GWAS using ~ 36,000 SNPs and determined that the RXFP2 gene was the main candidate for determining horn phenotype in sheep. Consistently, when comparing the genomes of global sheep breeds using the OvineSNP50K BeadChip, the strongest SNP under selection was determined to be located immediately adjacent to RXFP2, which suggested that RXFP2 can be potentially targeted by breeding for horn removal, being likely one of the oldest morphological changes occurring during sheep domestication . Johnston et al.  later described two RXFP2 alleles in Soay sheep, namely, Ho+ and HoP, with the Ho+ allele conferring larger, normal horns and the HoP allele conferring smaller horns. Ho+Ho+, Ho+HoP, and HoPHoP males were horned, and approximately half of the HoPHoP males were scurred. In contrast, the Ho+Ho+ genotype led to horned, Ho+HoP to scurred, and HoPHoP to polled females. In addition, it has been proposed that the Ho+ allele that confers larger horns is associated with higher reproductive success, whereas the smaller horn HoP allele is associated with increased survival. Collectively, these pieces of evidence suggest a trade-off between reproduction and survival for this allele in male Soay sheep. Several subsequent genomic studies highlighted the role of RXFP2 in determining horn phenotype when comparing different regional sheep populations. Pan et al.  described RXFP2 as a consistently strong marker of positive selection based on a GWAS in two semi-feral sheep populations from the Qinghai-Tibetan Plateau. In addition, the expression of RXFP2 in sheep has been studied in 13 tissues, and expression was exclusively in horn-related tissues . Moreover, an introgression of genomic regions into the RXFP2 locus in Tibetan sheep from Argali sheep has been suggested, in which two missense variants (c.1999C>T and c.1957C>T) were also detected . Recently, the same region was suggested to be introgressed from Iranian mouflon into Tibetan sheep . Nevertheless, in this region, three major highly divergent haplogroups were found in haplotype patterns across sheep breeds with different horn phenotypes .
Wiedemar et al.  described a 1833-bp genomic insertion (namely, an EEF1A1-like insertion) in the flanking 3′-region of the sheep RXFP2 gene that is potentially associated with polledness. In this study, this region was genotyped in seven Swiss sheep breeds. It was shown that horned individuals were homozygous for the absence of the insertion, whereas polled individuals were homozygous for its presence. It has been speculated that when the EEF1A1-like insertion is present, EEF1A1 transcripts bind to the 3′-UTR of the RXFP2 mRNA, leading to post-transcriptional downregulation [56, 57]. In Bündner Oberländer sheep with both polled and horned phenotypes, horned individuals were homozygous for the absence of the insertion; however, two polled individuals were found to be heterozygous for the insertion. Later, Lühken et al.  reported that the actual length of the EEF1A1-like insertion was 1780 bp. By genotyping this insertion in a large number of sheep breeds and individuals with different horn statuses, it was concluded that this insertion cannot be considered as the exclusive determinant of polledness in sheep, and thus could not be considered as a universal marker for determining the genetic causality of horns in sheep. The insertion was shown to segregate with horn status (such as polled or horned) only to a certain extent and within certain sheep breeds. In other breeds with variable horn statuses in both sexes or with a sex-dependent horn status, the segregation of this potential insertion is weaker. It has been suggested that this insertion is among the potential various causative factors for polledness in sheep . He et al.  performed PCR amplification of the DNA of 49 Chinese sheep with three horn statuses (polled, two-horned, and four-horned), and found the EEF1A1-like insertion in 68.4% (13/19) of the polled individuals and 43.3% (13/30) of the horned individuals. These findings suggest a yet unverified association between the EEF1A1-like insertion and the polledness phenotype in sheep.
Large-scale investigations have also been performed to elucidate the origins of mammalian horns [61, 62]. Wang et al.  identified 624 horn-related genes that are highly co-expressed in bone, skin, testis, and brain tissues. In addition, 201 genes were commonly found in both horn and antler tissues, and these genes were enriched in bone and skin development, and neurogenesis pathways. This study provided important insights into the origin of horns and further proposed a number of genes potentially associated with horn development regulation. Among these genes, SOX10, SNAI1, and TFAP2A were highly expressed in fetal sheep horn buds but not in adjacent skin tissues. In addition, this study also suggested a single evolutionary origin for the ruminant headgear, as well as proposed that the horn-stimulating gene RXFP2 is under convergent pseudogenization in secondarily headgearless lineages . Furthermore, using transcriptomic data, seven highly expressed keratin genes, namely, KRT1, KRT2, KRT3, KRT5, KRT10, KRT14, and KRT84, were found to be associated with the formation of keratinous sheath . Type II a-keratin proteins expressed by the above-mentioned genes (except KRT10 and KRT14) have been suggested to play an essential role in the development of the horn keratinous sheath . Further investigations on the sheep horn phenotype are shown in Additional file 2: Table S7. Collectively, these findings highlight the contribution of a gene set to the determination of sheep horn status; in particular, the prominent role of RXFP2 with the presence/absence of horns in sheep has been suggested. Although a significant association between RXFP2 and horn development has been indicated, the exact variants and the underlying mechanisms by which RXFP2 contributes to phenotypic differences in horn status in sheep have not yet been elucidated.
Beyond the horned or polled dichotomy: polyceraty in sheep
Polyceraty is related to the condition where the animal carries multiple horns (Fig. 6), and which is considered rare in sheep . Among native sheep populations, only few breeds carry multiple horns (Fig. 7), and polycerate sheep usually carry four or six horns. Ren et al.  performed a GWAS comparing 24 two-horned and 22 four-horned Chinese Sishui Fur sheep, which enabled the identification of a potential genomic region (132.0–133.1 Mb) in OAR2 that contained the top ten SNPs associated with the polycerate phenotype. In humans and mice, this genomic region includes the HOXD gene cluster, which comprises nine genes, including HOXD1, HOXD3, HOXD4, and HOXD8-HOXD13, as well as other functional genes, such as EVX2 and KIAA171, which was shown to have a potential association with the formation of limbs and genital buds . Subsequently, Kijas et al.  performed a GWAS mapping 125 two- and four-horned individuals from Jacob and Navajo-Churro sheep breeds to investigate the genetic cause of four-hornedness as well as examine whether it was also related with eyelid abnormality at various degrees of severity. This study confirmed that OAR2 harbors a strong marker linked with polyceraty (131.9–132.6 Mb). The MTX2 gene and the HOXD gene cluster (especially HOXD1) were the closest to the genetic marker that was most strongly associated with polyceraty. In Jacob and Navajo-Churro sheep breeds, eyelid abnormality was strictly observed in polled individuals and in polycerate individuals, which suggests an absence of eyelid abnormality in two-horned animals and a potential developmental link between eyelid abnormality and polyceraty ; thus, this study proposed a pleiotropic effect for the development of the horn trait in sheep. Greyvenstein et al.  conducted a GWAS on 43 Damara sheep (19 four-horned, 7 three-horned, 16 two-horned, and one polled), and confirmed that the proposed genomic region (128–135 Mb) spanning nine embryonic development HOXD genes is causal for the polyceraty phenotype.
He et al.  performed a GWAS on 34 two-horned and 32 four-horned individuals from three Chinese breeds (namely, Altay, Mongolian, and Sishui Fur). The potential genomic region (132.6–132.7 Mb) on OAR2 was similar to that previously identified in polycerate sheep breeds, i.e., Sishui Fur, Jacob, Navajo-Churro, and Damara, which suggests a common genomic cause for the polyceraty phenotype among different sheep breeds. In a subsequent study, the same research group performed another GWAS using 72 polycerate and 24 two-horned Tibetan sheep, to map the genomic regions of typical (four horns) and atypical polycerate (three, five, and six horns) phenotypes . The same candidate genomic region (131.990–133.525 Mb) was identified, thus suggesting that this region affects both typical and atypical polycerate phenotypes. In a more recent study, Allais-Bonnet et al.  re-visited the Illumina OvineHD Beadchip genotyping data of 111 cases and 87 control sheep previously reported [65, 66], and identified the polycerate locus between the positions 132,717,593 and 133,151,166 bp on OAR2. Subsequently, by comparing whole-genome sequences of 11 polycerate specimens and 1179 control sheep from various populations, a four-nucleotide deletion variant (g.132,832,249_132,832,252del) was revealed at position +4 to +7 bp after the first exon of the HOXD1 gene. Moreover, this potential variant further showed a perfect association when 236 animals from various populations, containing polycerate specimens, were investigated . This variant was proposed to influence the splicing of HOXD1 precursor RNAs, thus generating the multi-horn phenotype in sheep.
Sheep coat color: a major breed determinant
The ASIP gene
Sheep coat color has been investigated since a long time [68,69,70,71,72,73] and its variations have sparked the curiosity of a wide spectrum of researchers and specialists to understand its biological and genetic basis [6, 74, 75]. Coat colors in sheep vary widely among and within breeds (see Additional file 1: Fig. S2). In Australian Merino sheep, a uniform white fleece is a desirable trait, and the small proportion of colored-fleece lambs born annually in these flocks leads to significant financial loss for sheep breeders and wool processors. To reveal genes linked to pigmentation in the sheep genome, Parsons et al.  used comparative mapping information to investigate the ovine homologous position of the Agouti gene. The study showed that the locus responsible for recessive pigmentation in Australian Merino sheep is positioned within OAR13 in a region homologous to murine chromosome 2 and human chromosome 20 regions that contain the Agouti gene locus. In a subsequent report , the same authors used the available sequencing data of human and mouse homologs to isolate the coding region of the ovine Agouti gene, which revealed that the isolated ovine sequence has an authentic agouti homolog at the nucleotide and amino acid levels to other previously reported mammalian species. In addition, by developing a linkage map for loci responsible for certain phenotypic traits in Soay sheep, the ‘coat pattern’ (i.e., intensified coat color in the ‘self’ individuals vs. the contrast in color between the belly and the rest of the body in ‘wild’ individuals) mapped to OAR13 and close to the candidate Agouti locus .
Smit et al.  amplified a 142-bp segment in exon 2 of the Agouti gene (ASIP) using 361 samples obtained from sheep. Bands of 142 bp were present in all animals and bands of 137 bp (agouti-like sequence) were shown in 173 of the investigated sheep. The sequences of the 142-bp amplicons were identical to the previously published ovine agouti sequence  and a 5-bp deletion (D5; positions 98–102 from the start codon in the normal agouti transcript) was identified in the 137-bp amplicon. Subsequently, D5 was surveyed in ten sheep breeds, and the presence of D5 was observed in white-faced sheep breeds but not in black-faced sheep breeds . Royo et al.  reanalyzed the potential of D5 by sequencing the ASIP exon 2 and intron 2 using samples obtained from 188 individuals of the rare Xalda sheep breed, and the deletion (c.100_104delAGGAA) was verified in the surveyed individuals. All Xalda individuals carrying the deletion in a homozygous mode were phenotypically black; however, other black-coated Xalda individuals (91%) were not homozygous, which suggests that the homozygous deletion could not be the only determinant of black coat color in sheep. In addition, mRNA levels of ASIP was evaluated in 24 (14 black-coated and 10 white-coated) Xalda individuals, which revealed that ASIP expression in black-coated individuals was significantly lower than in white-coated individuals, and that D5 did not influence ASIP expression .
The role of ASIP in determining coat color pattern was primarily proposed by classical genetics. The dominant white or tan (AWt) ASIP allele is responsible for the yellow/red phenotype (caused by phaeomelanin), whereas the recessive allele non-agouti (Aa) is responsible for the black/brown phenotype (caused by eumelanin). In addition to other patterns such as badgerface (Ab) where the dorsal region is pale and the ventral region is darker. To reveal the genomic basis of these patterns, Norris et al.  performed sequence analysis of genomic DNA, bacterial artificial chromosome clones, and RT-PCR products from a group of sheep breeds with different coat colors. Using these molecular tools, the study proposed a model for the evolution of the ovine ASIP locus. This includes (i) a 190-kb tandem duplication encompassing the ASIP and AHCY genes and the ITCH promoter region, which was identified as the genetic determinant of the white coat color of dominant white/tan (AWt) agouti sheep; (ii) a single copy of ASIP with a silenced ASIP promoter was identified as the genetic cause of the recessive black sheep; and (iii) a single copy of wild-type ASIP (A+) with a functional promoter was identified as the genetic cause of the ancient Barbary sheep coat color (tan body and pale belly) . Notably, a 9-bp deletion (D9) was detected in this study in exon 2 of ASIP. Collectively, the findings of this study showed the complexity of the mechanism underlying ASIP regulation and provided novel insights into its role in determining sheep coat color . In addition, the copy number variation (CNV) in the Agouti locus was also used as a model for a segregation analysis to estimate allelic copy number .
Li et al.  performed a GWAS to identify the genetic determinants of the white vs. non-white coat color variation in a Finnsheep population. Among the 35 significant SNPs identified, 25 were located in or adjacent to five previously reported genes (ASIP, TYRP1, KIT, MC1R, and MITF) known to be involved in sheep coat color development. In addition, the s66432.1 SNP in the ASIP gene showed the strongest association signal with the development of coat color patterns between white and non-white individuals. In particular, the white coat color pattern was associated with the A allele, whereas the non-white coat color pattern was associated with the G allele. Additional studies on the sheep ASIP gene are reported in Additional file 2: Table S8. Interestingly, with the emergence of gene-editing endonucleases, the ASIP gene was targeted using the CRISPR/Cas9-based knockout approach [82, 83]. Among the generated Merino sheep, three individuals harboring a 4-bp deletion in the ASIP gene showed (i) a badger-face with black body coat color (n = 2) and (ii) brown coat color with light ventral pigmentation (n = 1), whereas (iii) two individuals with a 20-bp deletion showed a black-white spotted coat color  (see Additional file 1: Fig. S3). Collectively, these findings led to the conclusion that the CRISPR/Cas9-mediated ASIP disruption together with D9/D5 mutations and ASIP duplications contributed to the emergence of coat color in ASIP-targeted white Merino sheep. Thus, the manipulation of the ASIP gene using gene-editing tools enabled the determination of its functional role in sheep coat pigmentation.
The MC1R gene
It was expected that various loci or genes would epistatically interact to generate a wide variety of sheep coat colors. In addition to the Agouti locus that encodes ASIP, the Extension locus encoding the MC1R gene was shown to play a role in the formation of coat color in sheep . It plays the opposite role to ASIP, i.e., its high functional expression is responsible for the dark coat color phenotype. Våge et al.  analyzed the ovine MC1R gene in Norwegian Dala sheep, which comprise individuals with completely white or black coat colors. The study identified the complete co-segregation of c.218 T>A (p.73Met>Lys) and c.361G>A (p.121Asp > Asn) mutations with dominant black coat color (the dominant ED allele indicates the black coat pattern and the wild-type E+ indicates the white coat pattern). When these two mutations are induced in the mice MC1 receptor gene, the recapitulated p.73Met>Lys mutation alone was able to constitutively activate the receptor, whereas the p.121Asp >Asn mutation lacked this functionality, being thus proposedly required for high-affinity ligand binding . In addition, Våge et al.  confirmed the co-segregation of these two mutations (the dominant ED allele) with the black coat color pattern in Damara, Black Merino, and Black Corriedale. Moreover, Royo et al.  detected the presence of the dominant black ED allele in Karakul, Black Merino, and Black Castellana sheep.
In addition, a GWAS was performed using the OvineSNP50K BeadChip to determine the major pigmentation gene using 42 samples obtained from Manchega and Rasa Aragonesa sheep, in which the solid black and solid white coat colors segregate . In this study, the SNP s26449 was the most strongly associated with coat pigmentation and the closest to MC1R . In Brazilian Morada Nova sheep, a GWAS also revealed the occurrence of MC1R in the highly significant SNP-containing region . In this study, white-coated and red-coated animals were homozygous for the E+ allele, whereas black-coated animals carried at least one copy of the ED allele . In a study performed to identify pigmentation-specific genes associated with skin photosensitization, MC1R was also the most significant genomic region under selection, and a number of other genes were also suggested to be potentially associated with this trait, including KIT, PDGFRA, IRF4, SOX10, PICK1, and EDN3 . Moreover, a GWAS using the OvineSNP600K Beadchip was performed on 75 Chinese Tan sheep individuals (29 black-headed and 46 white-coated; see Additional file 1: Fig. S4), which revealed two significant SNPs (rs409651063 and rs408511664) in MC1R that were associated with coat color development . Further information from previous studies on the sheep MC1R gene is shown in Additional file 2: Table S9. Collectively, these results highlight the role of MC1R and its variants in the development of coat color pattern in sheep.
The TYRP1 gene
Considering the complexity of the coat color trait, it has been suggested that more than one gene could determine the phenotypic variation of coat colors. In Soay sheep, coat color density can be classified into dark and light patterns (in a ratio of approximately 3:1). Beraldi et al.  developed a complete linkage map to identify the loci responsible for determining several phenotypic traits, including coat color, which revealed TYRP1 as a strong candidate gene influencing sheep coat color. Subsequently, Gratten et al.  analyzed TYRP1 in Soay sheep, and showed that a nonsynonymous substitution (c.869G>T) in exon 4 was highly associated with coat color pattern. For this SNP, light-colored sheep were homozygous for T, whereas dark-colored sheep were either homozygous for G or heterozygous. In addition, no significant difference was found when analyzing TYRP1 expression in dark- and light-coated Soay sheep, which hindered the assessment of the potential role of this SNP in determining sheep coat color . Furthermore, an association between the TYRP1 c.869G>T mutation and body weight in Soay sheep has been reported . It has been shown that dark-coated sheep with GG and GT genotypes were heavier at birth than light-coated sheep with the TT genotype. Moreover, dark-coated sheep with the GG genotype showed a fitness disadvantage compared to dark-coated sheep with the GT genotype. In this case, it was also difficult to predict how TYPR1 influences the body weight or fitness . It has been proposed that other genetic factors in the vicinity of TYRP1 are likely responsible for microevolutionary constraints that influence these traits. It has also been suggested that the changing ratio between dark- and light-coated Soay sheep is a result of climate change [92, 93]; however, such a claim is not well evidenced .
In addition, Raadsma et al.  performed an extended QTL study on 13 skin and fiber pigmentation traits in Awassi × Merino × Merino backcross sheep. Significant QTL signals were identified on OAR2 where TYRP1 is located. Notably, the nonsynonymous TYRP1 SNP proposed by Gratten et al.  in Soay sheep was not confirmed in the examined flock. However, this study proposed an influence of the nonsynonymous SNP (c.2240C>G) in exon 2 of TYRP1 on coat pigmentation . In addition, two mutations, g.80608128G>T and g.80611700C>T within the TYRP1 gene have been reportedly associated with brown coat color development in American sheep populations . Furthermore, a GWAS was performed using Valais Red sheep with the aim to compare brown and black individuals, which showed a strong signal in the TYRP1 locus, and three variants were shown to be associated with the recessively inherited brown coat color, namely (i) the previously reported missense variant (c.869G>T), (ii) a frame-shift variant (c.86_87delGA), and (iii) a nonsense variant (c.1066C>T) . Taken together, these studies provide evidence for the potential involvement of TYRP1 in the development of coat pigmentation in sheep. Additional information on sheep coat color variation is shown in Additional file 2: Table S10, which highlights the potential contribution of a set of genes to the formation of coat color in sheep.
Black bone color or hyperpigmentation
Black bone color or hyperpigmentation has been observed in few sheep breeds. In this condition, the pigment is found not only in the coat but also in the skin, muscles, inner organs, and bones. In the indigenous Chinese Nanping black-boned sheep, two synonymous mutations (c.192G>C and c.462C>T) in exon 1 of the TYR gene, as well as two synonymous mutations (c.12A>G and c.144G>C) in the MC1R gene were detected . In addition, four combinations of a synonymous (c.90C>T) mutation and a missense (c.203C>T) mutation were found in the TYRP1 gene in both black-boned and normal sheep . Two other SNPs, a synonymous (c.102C>T) and a missense (c.215T>C) were also indicated in exon 2 of TYRP1 in black-boned sheep . Patterns of these SNPs showed elevated plasma tyrosinase activity [99,100,101]. In general, black-boned sheep show higher levels of melanin and related pigments in tissues and organs compared to the other examined breeds with white bone color . Another gene, EDN3, known to cause hyperpigmentation in black-boned chickens (the silky fowl), was also analyzed in sheep . However, the results of this study showed non-significant differences in genotype and allele frequencies of the identified SNPs and copy number of EDN3 gene among black-boned and non-black-boned sheep. Interestingly, a different bone pigmentation, e.g., red pigmentation, has also been described in goats [103,104,105]. Such pigmentation characteristics might influence carcass quality and standards and, subsequently, consumer preference.
Yellow fat color is an inherited recessive trait caused by the accumulation of carotenoids in adipose tissues [106,107,108,109,110]. This trait is known to influence carcass quality; thus it is of great interest to understand its genetic background. Våge et al.  analyzed two genes (BCMO1 and BCO2) known to be important for carotenoid degradation in Norwegian sheep, and highlighted the potential of a nonsense mutation (c.196C>T) to be associated with the emergence of the yellow fat phenotype. However, individuals who lacked this mutation had yellow fat, which suggested that another mechanism might be involved in the emergence of this trait. To further explore this, a more recent investigation was conducted in Norwegian spælsau sheep with aberrations in the BCO2 mRNA, which revealed the presence of a 7.9-kb insertion of an endogenous Jaagsiekte Sheep retrovirus (enJSRV) sequence into the first intron of BCO2 . This insertion led to the inactivation of BCO2 due to the generation of a truncated product of 58 amino acids (29 from BCO2 and 29 from enJSRV), whereas the functional BCO2 product is known to contain 576 amino acids . Interestingly, the ovine BCO2 gene was previously targeted using the CRISPR/Cas9 gene-editing system in a study aimed at generating sheep with multiplex editing in three genes (MSTN, ASIP, and BCO2) [82, 113]. The BCO2 gene was targeted to functionally validate its biological role, and was used as a gene with an obvious phenotype to validate the efficiency of CRISPR/Cas9 in sheep. The study revealed that the biallelic modification of BCO2 resulted in yellow fat compared to the monoallelic and wild-type individuals, confirming its functional role in fat color formation in sheep .
Genetic aspects of other phenotypic traits
Several genes are known to influence sheep wool quality and quantity, e.g., those included in the KAP families [114,115,116]. In the following sections, we will focus on genes related to the wool phenotype rather than genes that influence quantitatively wool aspects. In addition to the KAP families, a number of genes were shown to greatly influence fiber traits, e.g., FGF5. The inactivation of the FGF5 gene was first proven in mice to cause abnormally long hair . Natural variants in the sheep FGF5 gene causing the long hair phenotype have not yet been described . However, the potential of FGF5 has been investigated using gene modification tools. In order to induce wool growth by antagonizing the FGF5 function, transgenic Chinese Merino sheep with ectopic expression of short alternative FGF5 transcripts (FGF5s) were generated by zygotic injection of recombinant lentivirus; wool length was significantly longer in the transgenic group, and greasy wool weight was increased in yearling . In addition, direct knockout of the FGF5 was also achieved using the CRISPR/Cas9 system, and FGF5-disrupted sheep showed the desired long wool phenotype [120,121,122,123]. Interestingly, similar findings were also found in goats [124,125,126]. Using the FGF5-disrupted sheep model, a crosstalk between androgen and the Wnt/β-catenin signaling pathway in influencing wool increase and active hair-follicle density was revealed . This interplay has been linked to the activation of the c-MYC and KRT genes associated with the inner root sheath. Furthermore, the Shh pathway being located downstream of the Wnt/β-catenin signaling pathway, it was also shown to be involved in wool length . Moreover, a recent study on putative SNPs within the FGF5 full-length gene in fine-wool sheep revealed potential molecular markers associated with wool length and greasy wool weight . Based on these results, FGF5 has been highlighted as a potential gene that could lead to improving yield in wool-producing sheep breeds, thus promoting the development of the wool industry.
Chinese Tan sheep have curly fleece during the first months following birth (Fig. 8g), and this phenotype gradually disappears with aging. To understand this phenotypic change, a transcriptomic study involving sheep at two different developmental stages was performed, which enabled the identification of 600 differentially-expressed genes (DEG) . The top 20 highly expressed transcripts belonged to genes included in KAP families as well as other genes, namely, KAP3.2, KAP2.3, KRT25, KAP1.1, KRT5, KRT71, KRT14, KRT2.11, KRT83, KRT31, KRT34, KRT33, KRT27, KRT33B, TCHH, RPS8, PRR9, CST6, ASIP, and CO3. Furthermore, isoform analysis revealed that the MT3 isoform was upregulated in lamb skin compared to adult individuals, suggesting its potential involvement in the development of curly fleece phenotype . In addition, using an integrated miRNA-mRNA analysis of samples obtained from sheep at these two physiological stages, 49 differentially-expressed (DE) miRNAs were identified . In addition, 36 potential miRNA-mRNA target pairs were identified (including 25 DE miRNAs and 10 DEG). Among DE miRNAs, it was confirmed that miR-432 could inhibit KRT83 expression . Moreover, using suppression subtractive hybridization analysis, 67 DEG were identified in the samples of sheep at two developmental stages, and KRT71 was significantly highly expressed in individuals within the curly fleece group compared to the non-curling group . The MFZ1 transcription factor has also been suggested to regulate the expression of KRT71 . It has also been suggested that KRT83 expression was higher in Chinese Tan lambs with curly fleece compared to adults, and that the expression of the CAP1 transcription factor was low . When CAP1 was overexpressed, KRT83 expression was repressed, which suggests that CAP1 has a role in the regulation of KRT83 expression .
In a genome-wide detection of CNV on 48 Chinese Tan sheep, a CNV region overlapping the homeobox transcription factor DLX3 was suggested as a candidate CNV for the curly fleece phenotype . SNPs (c.118T>C, c.228T>C, c.688A>G, and c.1,038_1,039 insC) within the 3′-UTR of DLX3 were also previously described as associated with wool crimp in sheep . In Chinese Merino sheep, a set of genes (namely, CDC6, RARA, IGFBP4, and GHR) was identified as candidates for crimp number, in particular IGFBP4 . In addition, the levels of PAPPA2 expression in dermal papilla cells (DPC) of straight-wool Hu lambs was significantly higher than in DPC of curved wool Hu lambs . It has been speculated that PAPPA2 affects the proliferation of DPC by regulating IGFBP5 expression, which influences coat phenotypic patterns . In Hu sheep, miR-143 has been speculated to influence wool curvature by targeting CUX1 thereby inhibiting the proliferation of DPC and KRT71, which ultimately regulates the growth and development of hair follicles . The higher expression of CUX1 and KRT71 was proposed to be related to wool curvature . Furthermore, an analysis of DNA methylation levels in Chinese Tan sheep at three growth stages (1-, 24-, and 48-month old) revealed 51 differentially methylated genes (DMG) associated with age growth and curly fleece formation . Among these DMG, KRT71 and CD44 were highly methylated in one-month-old lambs, whereas ROR2 and ZDHHC13 were highly methylated in 48-month-old sheep. Interestingly, a similar phenotype of curly fibers was also found in goats [138,139,140,141]. Further information on the curly fleece phenotype in sheep is provided in Additional file 2: Table S11. Several genes and small regulatory RNAs were proposed to be involved in the formation of the curly fleece phenotype. Understanding the genetic basis of wool-related traits in sheep is of significant importance for the genetic customization of wool characteristics and quality for economic purposes.
Fine vs. coarse fleece
Wool fiber diameter varies widely among sheep breeds. Fine wool fiber is characterized by a small diameter, whereas coarse wool fiber has a larger diameter compared to fine fiber. Fine wool is employed in manufacturing soft clothes, whereas coarse wool is employed in carpet weaving. As it is significantly determinant for wool characteristics, early research efforts have concentrated on understanding the genetic determinants leading to differences in wool fiber diameter [142,143,144,145]. In a GWAS using 765 Chinese Merino sheep, nine significant SNPs were suggested to be associated with wool fiber diameter; among these, four SNPs were located within previously known or predicted genes, including TSPEAR, PIK3R4, KRTCAP3, and YWHAZ . These genes were found to be directly or indirectly associated with skin development, and their role in wool follicle development and fiber diameter has been proposed. To determine the wool fiber diameter coefficient of variation, three significant SNPs were detected, among which one was found to be located within the KIF16B gene . In addition, when 18 markers under selection in fine-wool Australian Merino and coarse-wool Spanish Churra sheep breeds were analyzed, divergent variants potentially involved in phenotypic variation of wool type, growth, and meat production/quality traits were identified . However, in this study, the investigated two breeds had varied phenotypic traits (wool, color, carcass, milk production, parasite resistance, among others), thus the identified markers under selection are likely not exclusively related to wool traits . In addition, in another study, BNC1 has been suggested as significantly associated with the mean fiber diameter in Chinese Merino sheep .
To further understand wool fiber diameter at the transcriptome level, a study was conducted using Gansu Alpine Merino sheep, which have wool of two fiber diameters (super-fine and fine), were analyzed, and 40 DEG were identified . Also, when analyzing natural antisense transcripts, three of 40 DEG (AQP8, Bos_d2, and SPRR) had significant natural antisense transcripts and were all significantly downregulated in the investigated super-fine wool group, thus indicating the correlation between gene transcripts and natural antisense transcripts in gene expression regulation . Additional information on wool diameter is provided in Additional file 2: Table S12. Collectively, the findings discussed above point towards the existence of sets of potential genes that regulate wool fiber diameter, but further studies are required to elucidate their underlying genetic mechanism.
Hairy vs. woolly fleece
Several modern domestic sheep breeds have been selected for short and woolly fleece coat phenotype, which differs substantially from the ancestral long and hairy fleece coat phenotype. To identify the molecular basis of this evolutionary phenotypic variation, Demars et al.  analyzed the Romane sheep breed in which fleece pattern segregates, to determine the causal mutation that controls fleece variation. An insertion of an antisense EIF2S2 retrogene (asEIF2S2; ~ 1.5 kb) into the 3′-UTR of IRF2BP2 was identified as leading to the production of an abnormal IRF2BP2/asEIF2S2 transcript. This chimeric messenger targets the genuine sense of EIF2S2 RNA and alters the expression of both EIF2S2 and IRF2BP2 mRNA by a retroposition mechanism-led RNA-RNA hybridization. It was shown that woolly sheep carried an insertion (IRF2BP2asEIF2S2 allele), whereas hairy sheep carried the ancestral IRF2BP2wt allele . This study provided the keys to understand the unique mechanism that underlies a causal variant leading to phenotype variation related to fiber folliculogenesis in sheep . Lv et al.  identified a novel mutation in OAR25 (g.7068586T>C) in the 3′-UTR of IRF2BP2 as a causal variant influencing fleece fiber diameter. In addition, the expression of a bindable miRNA (oar-miR-20a-3p) was detected exclusively in the skin tissues of coarse-wool sheep carrying the ancestral T allele. Furthermore, a strong selection marker in the secondary wool follicle-associated VEGFA gene was detected, which was shown to be co-activated by IRF2BP2. It has been proposed that oar-miR-20a-3p binds to the wild-type 3′-UTR of IRF2BP2, regulating its expression and subsequently inhibiting VEGFA expression in the skin tissues of coarse-wool sheep .
Interestingly, when genotyping the hairy/woolly locus in a ~ 1600-year-old naturally mummified sheep found in a salt mine (Chehrābād) in northwestern Iran, the ancestral “hairy” allele (likely IRF2BP2wt/wt or IRF2BP2wt/asEIF2S2) was identified . In addition, using scanning electron microscopic observations of the mummified hair fibers showed characteristics that are specific to sheep hair fibers of mouflon and medium-wool breeds. This suggested that, during that ancient period, the herd was maintained for meat or milk production rather than wool extraction, and ovicaprids were destined to be used as food source for workers in that location . This study highlighted the potentiality of using a locus for the phenotypic characterization of ancestral breeds using ancient DNA. Fleece pattern in sheep influences adaptation to environmental conditions, thus it is expected that the discovered variants related to different fleece patterns will be used in selection- or molecular genetics-based strategies to customize sheep breeds in order to better adapt to future environmental changes and sustainable sheep production.
Sheep ear morphology varies widely in terms of length (large, medium, small or “microtia”, and earless or “anotia”) and erectness (pricked or “erect” and dropped or “floppy”), and is considered a sheep breed determinant (Fig. 8b–e). In sheep breeds with varying ear length, such as Awassi, sheep farmers prefer to breed individuals with the normal large ear phenotype owing to their better physical constitution. In addition, it is indicated that earless Awassi sheep tend to show a more flighty behavior, which likely results from impaired hearing ability . Ear phenotype has also been associated with environmental adaptation, which highlights the importance of investigating the molecular basis of the ear phenotype in sheep. A SNP within the MSRB3 gene was to strongly segregate within the Chinese Duolang sheep with large ear phenotype and not in other breeds with normal or smaller ears, thus suggesting a role for MSRB3 in determining the ear phenotype in sheep . In a GWAS performed to investigate microtia in Jordanian Awassi sheep (in which 12 normal-eared and 8 earless individuals were included), a SNP adjacent to GATA6 was identified . GATA6 was proposed based on its location and function, which is related to developmental processes associated with outer ear shaping (e.g., chondrogenesis, the process formation of cartilage) .
To further investigate genes that influence the ear phenotype in Duolang sheep, a GWAS was carried out using 115 adult individuals with different ear sizes . The most significant SNP was located within the DCC gene which is known to be associated with ear development; two other significant SNPs (PTPRD and SOX5) showed a physical association with ear development. In this study, no association with GATA6 was found . When genotyping 40 individuals (20 with microtia and 20 normal) of the Valle del Belice sheep breed, a strong marker within intron 3 of CLRN1 was identified . CLRN1 is associated with Usher syndrome type 3A, which is a progressive hearing and vision loss, and was thus proposed as a novel candidate gene for determining microtia in sheep . In another GWAS using Chinese Altay sheep (26 microtia individuals and 29 normal-ear size individuals) combined with direct genotyping of a larger number of individuals, a 76-bp duplication in the evolutionarily conserved region (ECR) of HMX1 was determined . This gene is known to be associated with a congenital deformity of the outer ear in species such as cattle and rats, and was proposed to determine microtia in sheep . Furthermore, by analyzing previously studied genes in 84 Awassi sheep (16 earless, 41 short-eared, and 27 normal-eared), the ECR 76-bp duplication of HMX1 and a missense SNP (g.34498242A>G) in exon 7 of GATA6 were shown to be significantly associated with microtia . In addition, using a pruned dataset of 515 sheep from 17 breeds with different ear phenotypes, the most relevant marker for determining sheep ear size was found within the previously described MSRB3 gene . The strongest markers for ear erectness (dropped vs. pricked ears) were also located within the 3’-region of MSRB3 . Potential introgressive haplotypes of MSRB3 were found to be associated with ear morphology . Xu et al.  reported potential novel genes associated with ear size, which included PTCH1 and EMX2. Information on genetic variants known to determine the ear phenotype can be used in genetic screening to improve the selection against undesirable ear phenotypes (e.g., microtia) in certain sheep breeds.
Interestingly, other regions (including e.g., the SUPT3H, STXBP5L, DENND1A, and VPS13B genes) that have been introgressed from wild Ovis species into domestic sheep were suggested to be potentially associated with facial shape traits (e.g., nose bridge breadth, nose morphology, chin dimples, and forehead protrusion) .
Number of nipples
The supernumerary nipple/teat phenotype (also known as polythelia; see Additional file 1: Fig. S5a–c) was observed in sheep and documented over a century ago [160,161,162,163,164,165,166]. Ewes with supernumerary nipples usually carry three, four, five, or six nipples, and extra nipples are extremely small and undeveloped. In the past, extra-nippled ewes were believed to be more fertile, and selective breeding for supernumerary nipples could have been developed to ultimately lead them from their rudimentary condition to fully functional milk-producing nipples . However, it has been suggested that this trait has no practical value in sheep production , and now it is considered a morphological defect that might impede milking efficiency by increasing milking time and the risk of injuries , as well as negative implications for overall udder health . A GWAS was performed using 144 ewes of the Chinese Wadi breed (75 cases with four teats and 69 control cases with two teats), which revealed 63 significant SNPs and a set of functional genes (BBX, SEMA3D, IRF2BP2, SETBP1, NAV3, CD47, LRP1B, CERS5, CASK, and CADM2); a significant SNP near the BBX and CD47 genes was determined to be highly associated with the supernumerary nipple phenotype . In addition, Gene Ontology (GO) enrichment analysis revealed a set of functional genes that are potentially involved in determining nipple number and development, including CASK, a gene that was characterized as a risk factor associated with breast cancer development .
In addition to the five genes associated with breast cancer (LRP1B, GRM3, MACROD2, SETBP1, and GPC3) , seven novel genes (PHGDH, KDM3A, GLIS3, FSHR, CSNV2, CSN1S1, and ROBO2) previously reported in mice were associated with mammary and nipple development in sheep . Furthermore, by performing GWAS using CNV data, the associated functional gene, GPC5, was found to be associated with nipple number . In Hu sheep, LHFP, DPYSL2, and TDP-43 genes were also found to be related to nipple number . The identification of genes and variants that control nipple number trait will likely enable the molecular characterization of sheep breeds for the desired nipple phenotype.
Number of thoracic and lumbar vertebrae
Vertebrae are the interlocking bones that constitute the spinal cord. Sheep, as a tetrapod, have five types of vertebrae, namely cervical, thoracic, lumbar, sacral, and coccygeal or caudal vertebrae (see Additional file 1: Fig. S6). In addition to studies performed to identify genes associated with caudal vertebrae in sheep [14, 33], it has been suggested that the number of thoracolumbar vertebrae is associated with carcass length and weight [171, 172]. Thus, it is of interest to understand the genetic basis of vertebrae-related traits. Six hundred Chinese Kazakh sheep (400 with a larger number of thoracic vertebrae and 200 with an average number of thoracic vertebrae) were investigated by DNA sequencing. Eleven putatively selected genomic regions related to the indicated phenotype were revealed, including the vertebral development-related VRTN gene and the HOXA gene family . Further investigation showed that VRTN is the major promising gene involved in the determination of the number of thoracic vertebrae in sheep. In addition, during fetal development, the expression level of VRTN was shown to be significantly higher in sheep with more thoracic vertebrae compared to sheep in the control group .
A recent study identified genes associated with the number of ribs in 36 Hu sheep with 14 ribs (supernumerary condition) and 36 Hu sheep with 13 ribs (normal condition), and a combined effect of CPOX (g.178730623T>G; Oar_rambouillet_v1.0), KCNH1 (g.75716237C>G), and CPQ (g.88323841A>G) was proposed to be associated with rib number . Furthermore, using whole-genome resequencing of Chinese Duolang sheep with a different number of lumber vertebrate (6 vs. 7) revealed a significant region under selection (49.68–49.74 Mb; Oar_v4.0) that is contained within the SFRP4 gene . In addition, two SNPs (rs600370085C>T and rs415133338A>G) within this gene were significantly associated with multi-lumbar vertebrae . Additional file 2: Tables S13 and S14 contains further information on SNPs associated with the number of thoracic and lumbar vertebrae in sheep. Increasing the understanding of genes and variants controlling this phenotype in sheep will enable the design of selective breeding programs for increasing carcass length and weight.
Body height and stature
Dwarfism, whether proportionate or disproportionate, is included among the conditions that influence sheep phenotype [176, 177]. Ancon sheep is one of the breeds known to be affected by dwarfism, as a result of a putative dwarf mutation, leading to individuals with elongated bodies and very short legs . This phenotype was subjected to selection by breeders since such affected sheep are unable to jump fences. Ancon sheep were raised for a brief period of time, but went extinct as their phenotype was no longer desired [179, 180]. Another physical condition including chondrodysplasia, characterized by dwarfism and angular deformities of forelimbs, was reported in Texel as a result of a 1-bp deletion of T (g.25513delT) at position 107 in exon 3 of the SLC13A1 gene . Another condition of dysplasia was also reported in the Brazilian Cabugi sheep but its potential genetic basis has not been elucidated; in this context, dwarf lambs are characterized by a small size, short legs, domed heads, and other abnormalities . The spider lamb syndrome that is associated with severe skeletal deformities (e.g., abnormal spines and long splayed legs) is characterized by a loss-of-function nonsynonymous T>A transversion (p.700Val>Glu) in the FGFR3 gene . Subsequently, it has been proposed that the decreased expression of FGFR3 increased the proliferation and growth of plate chondrocytes by increasing telomerase reverse transcriptase (TERT) expression levels in vitro, suggesting a potential mechanism for skeletal overgrowth in vivo . In another report, it has been proposed that heterozygosity of the inactivated FGFR3 gene was used to improve skeletal growth . Furthermore, several genomic analyses identified a selective sweep on OAR6 harboring the LCORL locus that is implicated in determining height and stature in sheep [186,187,188,189,190,191,192,193].
Generic notes on other phenotypic traits
Other phenotypic traits in sheep are not yet well understood at the genetic level. Among these are included belly bareness (referring to the naturally bare area of skin on the belly [194, 195]), breech bareness (referring to the area of naturally bare skin around the perineum [196,197,198]), skin wrinkles (the folds or creases in the skin [199,200,201]), albinism (the congenital absence of any pigmentation [202, 203]), dewlap (a fold of loose skin which hangs from the neck or the throat as found in Persian sheep; Fig. 8a), manes (the growth of long hair on the neck as found in Barbados Blackbelly sheep; Fig. 8h), wattles (skin appendages found in the cervical region [204, 205] as found in Madras Red sheep; Fig. 8i), udder conformation , and other wool-related traits [207, 208].
Other significant traits for which the genes involved are well known have a great influence on sheep phenotype. These include, e.g., MSTN that causes increased body weight by generating the double-muscle phenotype [209,210,211], callipyge that causes the muscular rump phenotype [212,213,214], and SOCS2 that is associated with body weight and size [215, 216]. Since these traits are production traits rather than phenotypic traits, we have not presented them in detail in this review. Collectively, it can be stated that all the above-mentioned, well-validated genes and their variants function in an orchestrated manner to give rise to the phenotypic variation among sheep breeds.
In the present review, several genes known to be or potentially involved in determining sheep phenotype variation have been discussed. Understanding the genetic basis of phenotypic trait variation in sheep is of great importance, since these are breed determinants and highly linked with production and welfare traits. In spite of the number of genes discussed here, the genetic link and expression patterns of several of these genes are not yet well understood. Further studies are required to identify novel genomic loci and genes for traits for which the genetic determinism has not yet been described. The identification of genes, their causative variants, and their expression characteristics will help to improve sheep breeding through genomic-based selection and modern molecular genetic tools.
Availability of data and materials
Chessa B, Pereira F, Arnaud F, Amorim A, Goyache F, Mainland I, et al. Revealing the history of sheep domestication using retrovirus integrations. Science. 2009;324:532–6.
Alberto FJ, Boyer F, Orozco-terWengel P, Streeter I, Servin B, de Villemereuil P, et al. Convergent genomic signatures of domestication in sheep and goats. Nat Commun. 2018;9:813.
Deng J, Xie X-L, Wang D-F, Zhao C, Lv F-H, Li X, et al. Paternal origins and migratory episodes of domestic sheep. Curr Biol. 2020;30:4085–95.
Moradi MH, Nejati-Javaremi A, Moradi-Shahrbabak M, Dodds KG, McEwan JC. Genomic scan of selective sweeps in thin and fat tail sheep breeds for identifying of candidate regions associated with fat deposition. BMC Genet. 2012;13:10.
Johnston S, McEwan J, Pickering N, Kijas J, Beraldi D, Pilkington J, et al. Genome-wide association mapping identifies the genetic basis of discrete and quantitative variation in sexual weaponry in a wild sheep population. Mol Ecol. 2011;20:2555–66.
Ryder ML. Fleece colour in sheep and its inheritance. Anim Breed Abstr. 1980;48:305–24.
Dalrymple BP, Kirkness EF, Nefedov M, McWilliam S, Ratnakumar A, Barris W, et al. Using comparative genomics to reorder the human genome sequence into a virtual sheep genome. Genome Biol. 2007;8:R152.
Kijas JW, Lenstra JA, Hayes B, Boitard S, Porto Neto LR, San Cristobal M, et al. Genome-wide analysis of the world’s sheep breeds reveals high levels of historic mixture and strong recent selection. PLoS Biol. 2012;10:e1001258.
Jiang Y, Xie M, Chen W, Talbot R, Maddox JF, Faraut T, et al. The sheep genome illuminates biology of the rumen and lipid metabolism. Science. 2014;344:1168–73.
Davenport KM, Bickhart DM, Worley K, Murali SC, Salavati M, Clark EL, et al. An improved ovine reference genome assembly to facilitate in-depth functional annotation of the sheep genome. Gigascience. 2022;11:giab096.
Wei C, Wang H, Liu G, Wu M, Cao J, Liu Z, et al. Genome-wide analysis reveals population structure and selection in Chinese indigenous sheep breeds. BMC Genomics. 2015;16:194.
Moioli B, Pilla F, Ciani E. Signatures of selection identify loci associated with fat tail in sheep. J Anim Sci. 2015;93:4660–9.
Li R, Gong M, Zhang X, Wang F, Liu Z, Zhang L, et al. The first sheep graph-based pan-genome reveals the spectrum of structural variations and their effects on tail phenotypes. BioRxiv. 2022. https://doi.org/10.1101/2021.12.22.472709.
Zhi D, Da L, Liu M, Cheng C, Zhang Y, Wang X, et al. Whole genome sequencing of Hulunbuir short-tailed sheep for identifying candidate genes related to the short-tail phenotype. G3 (Bethesda). 2018;8:377–83.
Allais-Bonnet A, Hintermann A, Deloche M-C, Cornette R, Bardou P, Naval-Sanchez M, et al. Analysis of polycerate mutants reveals the evolutionary co-option of HOXD1 for horn patterning in bovidae. Mol Biol Evol. 2021;38:2260–72.
Parsons YM, Fleet MR, Cooper DW. The Agouti gene: a positional candidate for recessive self-colour pigmentation in Australian Merino sheep. Aust J Agric Res. 1999;50:1099–103.
Våge DI, Klungland H, Lu D, Cone RD. Molecular and pharmacological characterization of dominant black coat color in sheep. Mamm Genome. 1999;10:39–43.
Beraldi D, McRae AF, Gratten J, Slate J, Visscher PM, Pemberton JM. Development of a linkage map and mapping of phenotypic polymorphisms in a free-living population of Soay sheep (Ovis aries). Genetics. 2006;173:1521–37.
Kalds P, Luo Q, Sun K, Zhou S, Chen Y, Wang X. Trends towards revealing the genomic architecture of sheep tail patterning: Promising genes and investigatory pathways. Anim Genet. 2021;52:799–812.
Yuan Z, Liu E, Liu Z, Kijas JW, Zhu C, Hu S, et al. Selection signature analysis reveals genes associated with tail type in Chinese indigenous sheep. Anim Genet. 2017;48:55–66.
Xu S-S, Ren X, Yang G-L, Xie X-L, Zhao Y-X, Zhang M, et al. Genome-wide association analysis identifies the genetic basis of fat deposition in the tails of sheep (Ovis aries). Anim Genet. 2017;48:560–9.
Ahbara A, Bahbahani H, Almathen F, Al Abri M, Agoub MO, Abeba A, et al. Genome-wide variation, candidate regions and genes associated with fat deposition and tail morphology in Ethiopian indigenous sheep. Front Genet. 2019;9:699.
Manzari Z, Mehrabani-Yeganeh H, Nejati-Javaremi A, Moradi MH, Gholizadeh M. Detecting selection signatures in three Iranian sheep breeds. Anim Genet. 2019;50:298–302.
Pan Z, Li S, Liu Q, Wang Z, Zhou Z, Di R, et al. Rapid evolution of a retro-transposable hotspot of ovine genome underlies the alteration of BMP2 expression and development of fat tails. BMC Genomics. 2019;20:261.
Mastrangelo S, Moioli B, Ahbara A, Latairish S, Portolano B, Pilla F, et al. Genome-wide scan of fat-tail sheep identifies signals of selection for fat deposition and adaptation. Anim Prod Sci. 2018;59:835–48.
Mastrangelo S, Bahbahani H, Moioli B, Ahbara A, Al Abri M, Almathen F, et al. Novel and known signals of selection for fat deposition in domestic sheep breeds from Africa and Eurasia. PLoS One. 2019;14:e0209632.
Li X, Yang J, Shen M, Xie X-L, Liu G-J, Xu Y-X, et al. Whole-genome resequencing of wild and domestic sheep identifies genes associated with morphological and agronomic traits. Nat Commun. 2020;11:2815.
Dong K, Yang M, Han J, Ma Q, Han J, Song Z, et al. Genomic analysis of worldwide sheep breeds reveals PDGFD as a major target of fat-tail selection in sheep. BMC Genomics. 2020;21:800.
Baazaoui I, Bedhiaf-Romdhani S, Mastrangelo S, Ciani E. Genome-wide analyses reveal population structure and identify candidate genes associated with tail fatness in local sheep from a semi-arid area. Animal. 2021;15:100193.
Ahbara AM, Robert C, Kebede A, Abebe A, Latairish S, Agoub MO, et al. Natural adaptation and human selection of northeast African sheep genomes. Genomics. 2022;114:110448.
Wang F, Shao J, He S, Guo Y, Pan X, Wang Y, et al. Allele-specific expression and splicing provides insight into the phenotypic differences between thin- and fat-tailed sheep breeds. J Genet Genomics. 2022;49:583–6.
Lagler DK, Hannemann E, Eck K, Klawatsch J, Seichter D, Russ I, et al. Fine-mapping and identification of candidate causal genes for tail length in the Merinolandschaf breed. BioRxiv. 2022. https://doi.org/10.1101/2022.02.27.481613.
Han J, Yang M, Guo T, Niu C, Liu J, Yue Y, et al. Two linked TBXT (brachyury) gene polymorphisms are associated with the tailless phenotype in fat-rumped sheep. Anim Genet. 2019;50:772–7.
Li X, He S-G, Li W-R, Luo L-Y, Yan Z, Mo D-X, et al. Genomic analyses of Pamir argali, Tibetan sheep, and their hybrids provide insights into chromosome evolution, phenotypic variation, and germplasm innovation. Genome Res. 2022. https://doi.org/10.1101/gr.276769.122.
Zhang Y, Huang W, Hayashi C, Gatesy J, McKittrick J. Microstructure and mechanical properties of different keratinous horns. J R Soc Interface. 2018;15:20180093.
Castle WE. Genetics of horns in sheep. J Hered. 1940;31:486–7.
Wood TB. Note on the inheritance of horns and face colour in sheep. J Agric Sci. 1905;1:364–5.
Wood TB. The inheritance of horns and face colour in sheep. J Agric Sci. 1909;3:145–54.
Arkell TR, Davenport CB. Horns in sheep as a typical sex-limited character. Science. 1912;35:375–7.
Arkell TR, Davenport CB. The nature of the inheritance of horns in sheep. Science. 1912;35:927–927.
Castle WE. Are horns in sheep a sex-limited character. Science. 1912;35:574–5.
Dolling C. The inheritance of horns in sheep. Proc Aust Soc Anim Prod. 1959;1:161–4.
Dolling C. Hornedness and polledness in sheep. I. The inheritance of polledness in the Merino. Aust J Agric Res. 1960;11:427–38.
Dolling C. Hornedness and polledness in sheep. II. The inheritance of horns in Merino ewes. Aust J Agric Res. 1960;11:618–27.
Dolling C. Hornedness na polledness in sheep. III. The inheritance of horns in Dorset Horn ewes. Aust J Agric Res. 1960;11:845–50.
Dolling C. Hornedness and polledness in sheep. IV. Triple alleles affecting horn growth in the Merino. Aust J Agric Res. 1961;12:353.
Dolling C, Carter N. Hornedness and polledness in sheep. V. Some effects of the horn alleles on characters in the Merino. Aust J Agric Res. 1961;12:483–97.
Dolling C. Hornedness and polledness in sheep. VI. The inheritance of polledness in the Border Leicester. Aust J Agric Res. 1968;19:649–55.
Montgomery GW, Henry HM, Dodds KG, Beattie AE, Wuliji T, Crawford AM. Mapping the Horns (Ho) locus in sheep: a further locus controlling horn development in domestic animals. J Hered. 1996;87:358–63.
Pickering NK, Johnson PL, Auvray B, Dodds KG, McEwan JC. Mapping the horns locus in sheep. Proc Assoc Advmt Anim Breed Genet. 2009;18:88–91.
Johnston SE, Beraldi D, McRae AF, Pemberton JM, Slate J. Horn type and horn length genes map to the same chromosomal region in Soay sheep. Heredity (Edinb). 2010;104:196–205.
Johnston SE, Gratten J, Berenos C, Pilkington JG, Clutton-Brock TH, Pemberton JM, et al. Life history trade-offs at a single locus maintain sexually selected genetic variation. Nature. 2013;502:93–5.
Pan Z, Li S, Liu Q, Wang Z, Zhou Z, Di R, et al. Whole-genome sequences of 89 Chinese sheep suggest role of RXFP2 in the development of unique horn phenotype as response to semi-feralization. Gigascience. 2018;7:giy019.
Hu X-J, Yang J, Xie X-L, Lv F-H, Cao Y-H, Li W-R, et al. The genome landscape of Tibetan sheep reveals adaptive introgression from Argali and the history of early human settlements on the Qinghai-Tibetan plateau. Mol Biol Evol. 2019;36:283–303.
Cheng H, Zhang Z, Wen J, Lenstra JA, Heller R, Cai Y, et al. Long divergent haplotypes introgressed from wild sheep are associated with distinct morphological and adaptive characteristics in domestic sheep. BioRxiv. 2022. https://doi.org/10.1101/2022.05.17.492311.
Wiedemar N, Drögemüller C. A 1.8-kb insertion in the 3′-UTR of RXFP2 is associated with polledness in sheep. Anim Genet. 2015;46:457–61.
Derks MFL, Steensma M. Review: balancing selection for deleterious alleles in livestock. Front Genet. 2021;12:2526.
Lühken G, Krebs S, Rothammer S, Küpper J, Mioč B, Russ I, et al. The 1.78-kb insertion in the 3′-untranslated region of RXFP2 does not segregate with horn status in sheep breeds with variable horn status. Genet Sel Evol. 2016;48:78.
Simon R, Drögemüller C, Lühken G. The complex and diverse genetic architecture of the absence of horns (polledness) in domestic ruminants, including goats and sheep. Genes (Basel). 2022;13:832.
He X, Zhou Z, Pu Y, Chen X, Ma Y, Jiang L. Mapping the four-horned locus and testing the polled locus in three Chinese sheep breeds. Anim Genet. 2016;47:623–7.
Wang Y, Zhang C, Wang N, Li Z, Heller R, Liu R, et al. Genetic basis of ruminant headgear and rapid antler regeneration. Science. 2019;364:eaav6335.
Chen L, Qiu Q, Jiang Y, Wang K, Lin Z, Li Z, et al. Large-scale ruminant genome sequencing provides insights into their evolution and distinct traits. Science. 2019;364:eaav6202.
Dýrmundsson ÓR. Four-hornedness; a rare peculiarity still found in Icelandic sheep. Icelandic Sheep Breeders North Am Newsl. 2005;9:6–8.
Ren X, Yang G-L, Peng W-F, Zhao Y-X, Zhang M, Chen Z-H, et al. A genome-wide association study identifies a genomic region for the polycerate phenotype in sheep (Ovis aries). Sci Rep. 2016;6:21111.
Kijas JW, Hadfield T, Naval Sanchez M, Cockett N. Genome-wide association reveals the locus responsible for four-horned ruminant. Anim Genet. 2016;47:258–62.
Greyvenstein OFC, Reich CM, van Marle-Koster E, Riley DG, Hayes BJ. Polyceraty (multi-horns) in Damara sheep maps to ovine chromosome 2. Anim Genet. 2016;47:263–6.
He X, Song S, Chen X, Song T, Lobsang T, Guan W, et al. Genome-wide association analysis reveals the common genetic locus for both the typical and atypical polycerate phenotype in Tibetan sheep. Anim Genet. 2018;49:142–3.
Roberts JAF. Colour inheritance in sheep. I. Black colour and badger-face pattern in welh mountain sheep. J Genet. 1924;14:367–74.
Roberts JAF. Colour inheritance in sheep. II. The piebald pattern of the piebald breed. J Genet. 1926;17:77–83.
Roberts JAF. Colour inheritance in sheep. III. Face and leg colour. J Genet. 1928;19:261–8.
Roberts JAF. Colour inheritance in sheep. IV. White colour, recessive black colour, recessive brown colour, badger-face pattern, and reversed badger-face pattern. J Genet. 1930;22:165–80.
Roberts JAF, White RG. Colour inheritance in sheep. V. Dominant black. J Genet. 1930;22:181–90.
Roberts JAF. Colour inheritance in sheep. VI. The genetic constitution of the wild Mouflon. J Genet. 1931;25:1–15.
Lundie RS. The genetics of colour in fat-tailed sheep: a review. Trop Anim Health Prod. 2011;43:1245–65.
Koseniuk A, Ropka-Molik K, Rubiś D, Smołucha G. Genetic background of coat colour in sheep. Arch Anim Breed. 2018;61:173–8.
Parsons Y, Fleet M, Cooper D. Isolation of the ovine agouti coding sequence. Pigment cell Res. 1999;12:394–7.
Smit MA, Shay TL, Beever JE, Notter DR, Cockett NE. Identification of an agouti-like locus in sheep. Anim Genet. 2002;33:383–5.
Royo LJ, Alvarez I, Arranz JJ, Fernández I, Rodríguez A, Pérez-Pardal L, et al. Differences in the expression of the ASIP gene are involved in the recessive black coat colour pattern in sheep: evidence from the rare Xalda sheep breed. Anim Genet. 2008;39:290–3.
Norris BJ, Whan VA. A gene duplication affecting expression of the ovine ASIP gene is responsible for white and black sheep. Genome Res. 2008;18:1282–93.
Henshall JM, Whan VA, Norris BJ. Reconstructing CNV genotypes using segregation analysis: combining pedigree information with CNV assay. Genet Sel Evol. 2010;42:34.
Li M-H, Tiirikka T, Kantanen J. A genome-wide scan study identifies a single nucleotide substitution in ASIP associated with white versus non-white coat-colour variation in sheep (Ovis aries). Heredity (Edinb). 2014;112:122–31.
Wang X, Niu Y, Zhou J, Yu H, Kou Q, Lei A, et al. Multiplex gene editing via CRISPR/Cas9 exhibits desirable muscle hypertrophy without detectable off-target effects in sheep. Sci Rep. 2016;6:32271.
Zhang X, Li W, Liu C, Peng X, Lin J, He S, et al. Alteration of sheep coat color pattern by disruption of ASIP gene via CRISPR Cas9. Sci Rep. 2017;7:8149.
Klungland H, Vage DI. Pigmentary switches in domestic animal species. Ann N Y Acad Sci. 2003;994:331–8.
Våge DI, Fleet MR, Ponz R, Olsen RT, Monteagudo LV, Tejedor MT, et al. Mapping and characterization of the dominant black colour locus in sheep. Pigment Cell Res. 2003;16:693–7.
Kijas JW, Serrano M, McCulloch R, Li Y, Salces Ortiz J, Calvo JH, et al. Genomewide association for a dominant pigmentation gene in sheep. J Anim Breed Genet. 2013;130:468–75.
Muniz MMM, Caetano AR, McManus C, Cavalcanti LCG, Façanha DAE, Leite JHGM, et al. Application of genomic data to assist a community-based breeding program: a preliminary study of coat color genetics in Morada Nova sheep. Livest Sci. 2016;190:89–93.
Baazaoui I, McEwan J, Anderson R, Brauning R, McCulloch A, Van Stijn T, et al. GBS data identify pigmentation-specific genes of potential role in skin-photosensitization in two Tunisian sheep breeds. Animals (Basel). 2019;10:5.
Gebreselassie G, Liang B, Berihulay H, Islam R, Abied A, Jiang L, et al. Genomic mapping identifies two genetic variants in the MC1R gene for coat colour variation in Chinese Tan sheep. PLoS One. 2020;15:e0235426.
Gratten J, Beraldi D, Lowder BV, McRae AF, Visscher PM, Pemberton JM, et al. Compelling evidence that a single nucleotide substitution in TYRP1 is responsible for coat-colour polymorphism in a free-living population of Soay sheep. Proc Biol Sci. 2007;274:619–26.
Gratten J, Wilson AJ, McRae AF, Beraldi D, Visscher PM, Pemberton JM, et al. A localized negative genetic correlation constrains microevolution of coat color in wild sheep. Science. 2008;319:318–20.
Maloney SK, Fuller A, Mitchell D. Climate change: is the dark Soay sheep endangered? Biol Lett. 2009;5:826–9.
Maloney SK, Fuller A, Mitchell D. A warming climate remains a plausible hypothesis for the decrease in dark Soay sheep. Biol Lett. 2010;6:680–1.
Gratten J, Wilson AJ, McRae AF, Beraldi D, Visscher PM, Pemberton JM, et al. No evidence for warming climate theory of coat colour change in Soay sheep: a comment on Maloney et al. Biol Lett. 2010;6:678–9.
Raadsma HW, Jonas E, Fleet MR, Fullard K, Gongora J, Cavanagh CR, et al. QTL and association analysis for skin and fibre pigmentation in sheep provides evidence of a major causative mutation and epistatic effects. Anim Genet. 2013;44:547–59.
Posbergh C, Huson H. Making Moorit: Mutations in TYRP1 are responsible for brown coat color in different United States sheep breeds. In: Proceedings of the 11th World Congress on Genetics Applied to Livestock Production: 11–16 February 2018; Auckland; 2018.
Paris JM, Letko A, Häfliger IM, Ammann P, Flury C, Drögemüller C. Identification of two TYRP1 loss-of-function alleles in Valais Red sheep. Anim Genet. 2019;50:778–82.
Deng WD, Yang SL, Huo YQ, Gou X, Shi XW, Mao HM. Physiological and genetic characteristics of black-boned sheep (Ovis aries). Anim Genet. 2006;37:586–8.
Deng WD, Xi DM, Gou X, Yang SL, Shi XW, Mao HM. Pigmentation in Black-boned sheep (Ovis aries): association with polymorphism of the Tyrosinase gene. Mol Biol Rep. 2008;35:379–85.
Deng WD, Shu W, Yang SL, Shi XW, Mao HM. Pigmentation in Black-boned sheep (Ovis aries): association with polymorphism of the MC1R gene. Mol Biol Rep. 2009;36:431–6.
Li G, Xiong H, Xi D, Memon S, Wang L, Liu X, et al. An examination of melanogenic traits and TYRP1 polymorphism in Nanping and Romney Marsh sheep breeds. Arch Anim Breed. 2018;61:131–41.
Darwish HYA, Zhang Y, Cui K, Yang Z, Han D, Dong X, et al. Molecular cloning and characterization of the endothelin 3 gene in black bone sheep. J Anim Sci Biotechnol. 2018;9:57.
Wu C, Li X, Han T, Li P, Wang J, Liu G, et al. Dietary pseudopurpurin improves bone geometry architecture and metabolism in red-bone Guishan goats. PLoS One. 2012;7:e37469.
Yue Y, Yang B, Li Y, Zhang W, Zhang H, Wang J, et al. Chinese indigenous goat breeds. In: Sustainable goat production in adverse environments. Vol II. Cham: Springer International Publishing; 2017. p. 41–54.
He Y, Hong Q, Zhou D, Wang S, Yang B, Yuan Y, et al. Genome-wide selective detection of Mile red-bone goat using next-generation sequencing technology. Ecol Evol. 2021;11:14805–12.
Castle WE. Yellow fat in sheep. J Hered. 1934;25:246.
Hill F. Xanthophyll pigmentation in sheep fat. Nature. 1962;194:865–6.
Crane B, Clare NT. Nature of carotenoid pigments in yellow fat of sheep. New Zeal J Agric Res. 1975;18:273–5.
Kirton AH, Crane B, Paterson DJ, Clare NT. Yellow fat in lambs caused by carotenoid pigmentation. New Zeal J Agric Res. 1975;18:267–72.
Baker RL, Steine T, Våben AW, Breines D. The inheritance and incidence of yellow fat in Norwegian sheep. Acta Agric Scand. 1985;35:389–97.
Våge DI, Boman IA. A nonsense mutation in the beta-carotene oxygenase 2 (BCO2) gene is tightly associated with accumulation of carotenoids in adipose tissue in sheep (Ovis aries). BMC Genet. 2010;11:10.
Kent M, Moser M, Boman IA, Lindtveit K, Árnyasi M, Sundsaasen KK, et al. Insertion of an endogenous Jaagsiekte sheep retrovirus element into the BCO2 - gene abolishes its function and leads to yellow discoloration of adipose tissue in Norwegian Spælsau (Ovis aries). BMC Genomics. 2021;22:492.
Niu Y, Jin M, Li Y, Li P, Zhou J, Wang X, et al. Biallelic β-carotene oxygenase 2 knockout results in yellow fat in sheep via CRISPR/Cas9. Anim Genet. 2017;48:242–4.
Gong H, Zhou H, Forrest R, Li S, Wang J, Dyer J, et al. Wool keratin-associated protein genes in sheep—a review. Genes (Basel). 2016;7:24.
Yu Z, Plowman JE, Maclean P, Wildermoth JE, Brauning R, McEwan JC, et al. Ovine keratome: identification, localisation and genomic organisation of keratin and keratin-associated proteins. Anim Genet. 2018;49:361–70.
Zhou H, Gong H, Wang J, Luo Y, Li S, Tao J, et al. The complexity of the ovine and caprine keratin-associated protein genes. Int J Mol Sci. 2021;22:12838.
Hébert JM, Rosenquist T, Götz J, Martin GR. FGF5 as a regulator of the hair growth cycle: evidence from targeted and spontaneous mutations. Cell. 1994;78:1017–25.
Zhang L, He S, Liu M, Liu G, Yuan Z, Liu C, et al. Molecular cloning, characterization, and expression of sheep FGF5 gene. Gene. 2015;555:95–100.
Li W, He S, Liu C, Zhang X, Wang L, Lin J, et al. Ectopic expression of FGF5s induces wool growth in Chinese merino sheep. Gene. 2017;627:477–83.
Hu R, Fan ZY, Wang BY, Deng SL, Zhang XS, Zhang JL, et al. Rapid communication: generation of FGF5 knockout sheep via the CRISPR/Cas9 system. J Anim Sci. 2017;95:2019–24.
Li W, Liu C, Zhang X, Chen L, Peng X, He S, et al. CRISPR/Cas9-mediated loss of FGF5 function increases wool staple length in sheep. FEBS J. 2017;284:2764–73.
Zhang R, Wu H, Lian Z. Bioinformatics analysis of evolutionary characteristics and biochemical structure of FGF5 Gene in sheep. Gene. 2019;702:123–32.
Zhang R, Li Y, Jia K, Xu X, Li Y, Zhao Y, et al. Crosstalk between androgen and Wnt/β-catenin leads to changes of wool density in FGF5-knockout sheep. Cell Death Dis. 2020;11:407.
Wang X, Yu H, Lei A, Zhou J, Zeng W, Zhu H, et al. Generation of gene-modified goats targeting MSTN and FGF5 via zygote injection of CRISPR/Cas9 system. Sci Rep. 2015;5:13878.
Wang X, Cai B, Zhou J, Zhu H, Niu Y, Ma B, et al. Disruption of FGF5 in cashmere goats using CRISPR/Cas9 results in more secondary hair follicles and longer fibers. PLoS One. 2016;11:e0164640.
Li G, Zhou S, Li C, Cai B, Yu H, Ma B, et al. Base pair editing in goat: nonsense codon introgression into FGF5 results in longer hair. FEBS J. 2019;286:4675–92.
Zhao H, Hu R, Li F, Yue X. Five SNPs within the FGF5 gene significantly affect both wool traits and growth performance in fine-wool sheep (Ovis aries). Front Genet. 2021;12:732097.
Kang X, Liu G, Liu Y, Xu Q, Zhang M, Fang M. Transcriptome profile at different physiological stages reveals potential mode for curly fleece in Chinese Tan sheep. PLoS One. 2013;8: e71763.
Liu Y, Zhang J, Xu Q, Kang X, Wang K, Wu K, et al. Integrated miRNA-mRNA analysis reveals regulatory pathways underlying the curly fleece trait in Chinese tan sheep. BMC Genomics. 2018;19:360.
Kang X, Liu Y, Zhang J, Xu Q, Liu C, Fang M. Characteristics and expression profile of KRT71 screened by suppression subtractive hybridization cDNA library in curly fleece Chinese Tan sheep. DNA Cell Biol. 2017;36:552–64.
Liu Y, Kang X, Yang W, Xie M, Zhang J, Fang M. Differential expression of KRT83 regulated by the transcript factor CAP1 in Chinese Tan sheep. Gene. 2017;614:15–20.
Ma Q, Liu X, Pan J, Ma L, Ma Y, He X, et al. Genome-wide detection of copy number variation in Chinese indigenous sheep using an ovine high-density 600 K SNP array. Sci Rep. 2017;7:912.
Rong E, Zhang Z, Qiao S, Yang H, Yan X, Li H, et al. Functional characterization of a single nucleotide polymorphism in the 3’ untranslated region of sheep DLX3 gene. PLoS One. 2015;10:e0137135.
Zhao B, Luo H, Huang X, Wei C, Di J, Tian Y, et al. Integration of a single-step genome-wide association study with a multi-tissue transcriptome analysis provides novel insights into the genetic basis of wool and weight traits in sheep. Genet Sel Evol. 2021;53:56.
Wu T, Wang S, Jin Q, Lv X, Sun W. PAPPA2 promote the proliferation of dermal papilla cells in Hu sheep (Ovis aries) by regulating IGFBP5. Genes (Basel). 2021;12:1490.
Hu T, Huang S, Lv X, Wang S, Getachew T, Mwacharo JM, et al. MiR-143 targeting CUX1 to regulate proliferation of dermal papilla cells in Hu sheep. Genes (Basel). 2021;12:2017.
Liu Y, Xu Q, Kang X, Wang K, Wang J, Feng D, et al. Dynamic changes of genomic methylation profiles at different growth stages in Chinese Tan sheep. J Anim Sci Biotechnol. 2021;12:118.
Li Y, Zhou G, Zhang R, Guo J, Li C, Martin G, et al. Comparative proteomic analyses using iTRAQ-labeling provides insights into fiber diversity in sheep and goats. J Proteomics. 2018;172:82–8.
Xiao P, Zhong T, Liu Z, Ding Y, Guan W, He X, et al. Integrated analysis of methylome and transcriptome changes reveals the underlying regulatory signatures driving curly wool transformation in Chinese Zhongwei goats. Front Genet. 2020;10:1263.
Ding Y, Xue X, Liu Z, Ye Y, Xiao P, Pu Y, et al. Expression profiling and functional characterization of miR-26a and miR-130a in regulating Zhongwei goat hair development via the TGF-β/SMAD pathway. Int J Mol Sci. 2020;21:5076.
Liu Y, Ding Y, Liu Z, Chen Q, Li X, Xue X, et al. Integration analysis of transcriptome and proteome reveal the mechanisms of goat wool bending. Front Cell Dev Biol. 2022;10:836913.
Parsons YM, Cooper DW, Piper LR. Evidence of linkage between high-glycine-tyrosine keratin gene loci and wool fibre diameter in a Merino half-sib family. Anim Genet. 1994;25:105–8.
Allain D, Lantier I, Elsen JM, François D, Brunel JC, Weisbecker JL, et al. A design aiming at detecting QTL controlling wool traits and other traits in the INRA401 sheep line. In: Proceedings of the 6th World Congress on Genetics Applied to Livestock Production: 11–16 January 1998; Armidale; 1998.
Purvis I, Franklin I. Major genes and QTL influencing wool production and quality: a review. Genet Sel Evol. 2005;37:S97–107.
Purvis IW, Jeffery N. Genetics of fibre production in sheep and goats. Small Rumin Res. 2007;70:42–7.
Wang Z, Zhang H, Yang H, Wang S, Rong E, Pei W, et al. Genome-wide association study for wool production traits in a Chinese Merino sheep population. PLoS One. 2014;9:e107101.
Gutiérrez-Gil B, Esteban-Blanco C, Wiener P, Chitneedi PK, Suarez-Vega A, Arranz J-J. High-resolution analysis of selection sweeps identified between fine-wool Merino and coarse-wool Churra sheep breeds. Genet Sel Evol. 2017;49:81.
Yue Y, Guo T, Liu J, Guo J, Yuan C, Feng R, et al. Exploring differentially expressed genes and natural antisense transcripts in sheep (Ovis aries) skin with different wool fiber diameters by digital gene expression profiling. PLoS One. 2015;10:e0129249.
Demars J, Cano M, Drouilhet L, Plisson-Petit F, Bardou P, Fabre S, et al. Genome-wide identification of the mutation underlying fleece variation and discriminating ancestral hairy species from modern woolly sheep. Mol Biol Evol. 2017;34:1722–9.
Caspermeyer J. So sheep may safely graze: honing in on culprit behind fleece variation in domesticated sheep. Mol Biol Evol. 2017;34:1821–1821.
Lv F-H, Cao Y-H, Liu G-J, Luo L-Y, Lu R, Liu M-J, et al. Whole-genome resequencing of worldwide wild and domestic sheep elucidates genetic diversity, introgression and agronomically important loci. Mol Biol Evol. 2021;39:msab353.
Rossi C, Ruß-Popa G, Mattiangeli V, McDaid F, Hare AJ, Davoudi H, et al. Exceptional ancient DNA preservation and fibre remains of a Sasanian saltmine sheep mummy in Chehrābād. Iran Biol Lett. 2021;17:20210222.
Jawasreh K, Boettcher PJ, Stella A. Genome-wide association scan suggests basis for microtia in Awassi sheep. Anim Genet. 2016;47:504–6.
Gao L, Xu S-S, Yang J-Q, Shen M, Li M-H. Genome-wide association study reveals novel genes for the ear size in sheep (Ovis aries). Anim Genet. 2018;49:345–8.
Mastrangelo S, Sottile G, Sutera AM, Di Gerlando R, Tolone M, Moscarelli A, et al. Genome-wide association study reveals the locus responsible for microtia in Valle del Belice sheep breed. Anim Genet. 2018;49:636–40.
He S, Zhang Z, Sun Y, Ren T, Li W, Zhou X, et al. Genome-wide association study shows that microtia in Altay sheep is caused by a 76 bp duplication of HMX1. Anim Genet. 2020;51:132–6.
Jawasreh KIZ, Al-Omari HD-A. Duplication in ECR near HMX1 and a SNP in GATA6 genes regulate microtia in Awassi sheep. Genes (Basel). 2020;11:597.
Paris JM, Letko A, Häfliger IM, Ammann P, Drögemüller C. Ear type in sheep is associated with the MSRB3 locus. Anim Genet. 2020;51:968–72.
Xu S-S, Gao L, Shen M, Lyu F. Whole-genome selective scans detect genes associated with important phenotypic traits in sheep (Ovis aries). Front Genet. 2021;12:738879.
Bell AG. Bell on the development by selection of supernumerary mammæ in sheep. Science. 1899;9:637–9.
Bell AG. The multi-nippled sheep of Beinn Bhreagh. Science. 1904;19:767–8.
Bell AG. Sheep-breeding experiments on Beinn Bhreagh. Science. 1912;36:378–84.
Bell AG. Saving the six-nippled breed: Mr. Bell’s last contribution to science, with an introduction by Mrs. Bell. J Hered. 1923;14:99–111.
Castle WE. The genetics of multi-nippled sheep—an analysis of the sheep-breeding experiments of Dr. and Mrs Alexander Graham Bell at Beinn Bhreagh, NS. J Hered. 1924;15:75–85.
Phillips RW, Schott RG, Spencer DA. The genetics, physiology, and economic importance of the multinipple trait in sheep. Washington: United States Department of Agriculture. Tech Bull. 1945;909:1–16.
Phillips RW, Schott RG, Spencer DA. The multinipple trait in sheep. J Hered. 1946;37:19–26.
Martin P, Palhière I, Tosser-Klopp G, Rupp R. Heritability and genome-wide association mapping for supernumerary teats in French Alpine and Saanen dairy goats. J Dairy Sci. 2016;99:8891–900.
Hardwick LJA, Phythian CJ, Fowden AL, Hughes K. Size of supernumerary teats in sheep correlates with complexity of the anatomy and microenvironment. J Anat. 2020;236:954–62.
Peng W-F, Xu S-S, Ren X, Lv F-H, Xie X-L, Zhao Y-X, et al. A genome-wide association study reveals candidate genes for the supernumerary nipple phenotype in sheep (Ovis aries). Anim Genet. 2017;48:570–9.
Zhao Y, Pu Y, Liang B, Bai T, Liu Y, Jiang L, et al. A study using single-locus and multi-locus genome-wide association study to identify genes associated with teat number in Hu sheep. Anim Genet. 2022;53:203–11.
Li C, Zhang X, Cao Y, Wei J, You S, Jiang Y, et al. Multi-vertebrae variation potentially contribute to carcass length and weight of Kazakh sheep. Small Rumin Res. 2017;150:8–10.
Zhang Z, Sun Y, Du W, He S, Liu M, Tian C. Effects of vertebral number variations on carcass traits and genotyping of Vertnin candidate gene in Kazakh sheep. Asian Australas J Anim Sci. 2017;30:1234–8.
Li C, Li M, Li X, Ni W, Xu Y, Yao R, et al. Whole-genome resequencing reveals loci associated with thoracic vertebrae number in sheep. Front Genet. 2019;10:674.
Zhang D, Zhang X, Li F, Liu T, Hu Z, Gao N, et al. Whole-genome resequencing identified candidate genes associated with the number of ribs in Hu sheep. Genomics. 2021;113:2077–84.
Li C, Liu K, Dai J, Li X, Liu X, Ni W, et al. Whole-genome resequencing to investigate the determinants of the multi-lumbar vertebrae trait in sheep. Gene. 2021;809:146020.
Thompson KG, Piripi SA, Dittmer KE. Inherited abnormalities of skeletal development in sheep. Vet J. 2008;177:324–33.
Boegheim IJM, Leegwater PAJ, van Lith HA, Back W. Current insights into the molecular genetic basis of dwarfism in livestock. Vet J. 2017;224:64–75.
Humphreys D. On a new variety in the breeds of sheep. Philos Trans R Soc Lond. 1813;103:88–95.
Landauer W, Chang TK. The Ancon or Otter sheep: history and genetics. J Hered. 1949;40:105–12.
Gidney L. Earliest archaeological evidence of the Ancon mutation in sheep from Leicester. UK Int J Osteoarchaeol. 2007;17:318–21.
Zhao X, Onteru SK, Piripi S, Thompson KG, Blair HT, Garrick DJ, et al. In a shake of a lamb’s tail: using genomics to unravel a cause of chondrodysplasia in Texel sheep. Anim Genet. 2012;43:9–18.
Dantas F, Medeiros G, Figueiredo A, Thompson K, Riet-Correa F. Skeletal dysplasia with craniofacial deformity and disproportionate dwarfism in hair sheep of northeastern Brazil. J Comp Pathol. 2014;150:245–52.
Beever JE, Smit MA, Meyers SN, Hadfield TS, Bottema C, Albretsen J, et al. A single-base change in the tyrosine kinase II domain of ovine FGFR3 causes hereditary chondrodysplasia in sheep. Anim Genet. 2006;37:66–71.
Smith LB, Belanger JM, Oberbauer AM. Fibroblast growth factor receptor 3 effects on proliferation and telomerase activity in sheep growth plate chondrocytes. J Anim Sci Biotechnol. 2012;3:39.
Smith LB, Dally MR, Sainz RD, Rodrigue KL, Oberbauer AM. Enhanced skeletal growth of sheep heterozygous for an inactivated fibroblast growth factor receptor 3. J Anim Sci. 2006;84:2942–9.
Al-Mamun HA, Kwan P, Clark SA, Ferdosi MH, Tellam R, Gondro C. Genome-wide association study of body weight in Australian Merino sheep reveals an orthologous region on OAR6 to human and bovine genomic regions affecting height and weight. Genet Sel Evol. 2015;47:66.
Matika O, Riggio V, Anselme-Moizan M, Law AS, Pong-Wong R, Archibald AL, et al. Genome-wide association reveals QTL for growth, bone and in vivo carcass traits as assessed by computed tomography in Scottish Blackface lambs. Genet Sel Evol. 2016;48:11.
Bolormaa S, Hayes BJ, van der Werf JHJ, Pethick D, Goddard ME, Daetwyler HD. Detailed phenotyping identifies genes with pleiotropic effects on body composition. BMC Genomics. 2016;17:224.
Manunza A, Cardoso TF, Noce A, Martínez A, Pons A, Bermejo LA, et al. Population structure of eleven Spanish ovine breeds and detection of selective sweeps with BayeScan and hapFLK. Sci Rep. 2016;6:27296.
Rochus CM, Tortereau F, Plisson-Petit F, Restoux G, Moreno-Romieux C, Tosser-Klopp G, et al. Revealing the selection history of adaptive loci using genome-wide scans for selection: an example from domestic sheep. BMC Genomics. 2018;19:71.
Naval-Sanchez M, Nguyen Q, McWilliam S, Porto-Neto LR, Tellam R, Vuocolo T, et al. Sheep genome functional annotation reveals proximal regulatory elements contributed to the evolution of modern breeds. Nat Commun. 2018;9:859.
Ruiz-Larrañaga O, Langa J, Rendo F, Manzano C, Iriondo M, Estonba A. Genomic selection signatures in sheep from the Western Pyrenees. Genet Sel Evol. 2018;50:9.
Signer-Hasler H, Burren A, Ammann P, Drögemüller C, Flury C. Runs of homozygosity and signatures of selection: a comparison among eight local Swiss sheep breeds. Anim Genet. 2019;50:512–25.
Craven A, Scobie D, Ashby M, Nixon A. Wool follicle characteristics in sheep with differing ventral wool length phenotypes. Proc New Zeal Soc Anim Prod. 2008;68:49–52.
Scobie DR, O’Connell D, Morris CA, Hickey SM. Inherited bareness of the belly reduces the time taken to shear ewes and hoggets. Anim Prod Sci. 2011;51:176–82.
Scobie DR, O’Connell D, Morris CA, Hickey SM. A preliminary genetic analysis of breech and tail traits with the aim of improving the welfare of sheep. Aust J Agric Res. 2007;58:161–7.
Scobie DR, O’Connell D, Morris CA, Hickey SM, Scobie DR, O’Connell D, et al. Dag score is negatively correlated with breech bareness score of sheep. Aust J Exp Agric. 2008;48:999–1003.
Scobie DR, O’Connell D. Breech bareness reduces flystrike in New Zealand crossbred sheep. Anim Prod Sci. 2010;50:599–602.
Scobie DR, Young SR, O’Connell D, Gurteen S. Skin wrinkles affect wool characteristics and the time taken to harvest wool from Merino and Halfbred sheep. New Zeal J Agric Res. 2005;48:177–85.
Scobie DR, Young SR, O’Connell D, Eythorsdottir E. Skin wrinkles of the sire adversely affect Merino and halfbred pelt characteristics and other production traits. Aust J Exp Agric. 2006;45:1551–7.
Scholtz AJ, Cloete SWP, van Wyk JB, Misztal I, du Toit E, van der Linde TC de K, et al. Genetic (co)variances between wrinkle score and absence of breech strike in mulesed and unmulesed Merino sheep, using a threshold model. Anim Prod Sci. 2010;50:210–8.
Adalsteinsson S. Albinism in Icelandic sheep. J Hered. 1977;68:347–9.
Rowett MA, Fleet MR. Albinism in a Suffolk sheep. J Hered. 1993;84:67–9.
Wassin B. Linkage studies in sheep. J Hered. 1931;22:9–13.
Lancaster MJ, Medwell WD. Neck wattles in lambs. Aust Vet J. 1991;68:75–6.
Legarra A, Ugarte E. Genetic parameters of udder traits, somatic cell score, and milk yield in Latxa sheep. J Dairy Sci. 2005;88:2238–45.
Ponzoni RW, Rogan IM, James PJ. Genetic improvement of apparel and carpet wool production. In: Proceedings of the Fourth World Congress on Genetics Applied to Livestock Production: 23–27 July 1990; Edinburgh; 1990.
Bolormaa S, Swan AA, Brown DJ, Hatcher S, Moghaddar N, van der Werf JH, et al. Multiple-trait QTL mapping and genomic prediction for wool traits in sheep. Genet Sel Evol. 2017;49:62.
Clop A, Marcq F, Takeda H, Pirottin D, Tordoir X, Bibé B, et al. A mutation creating a potential illegitimate microRNA target site in the myostatin gene affects muscularity in sheep. Nat Genet. 2006;38:813–8.
Boman IA, Klemetsdal G, Blichfeldt T, Nafstad O, Våge DI. A frameshift mutation in the coding region of the myostatin gene (MSTN) affects carcass conformation and fatness in Norwegian White Sheep (Ovis aries). Anim Genet. 2009;40:418–22.
Boman IA, Klemetsdal G, Nafstad O, Blichfeldt T, Våge DI. Impact of two myostatin (MSTN) mutations on weight gain and lamb carcass classification in Norwegian White Sheep (Ovis aries). Genet Sel Evol. 2010;42:4.
Cockett NE, Jackson SP, Shay TL, Nielsen D, Moore SS, Steele MR, et al. Chromosomal localization of the callipyge gene in sheep (Ovis aries) using bovine DNA markers. Proc Natl Acad Sci USA. 1994;91:3019–23.
Cockett NE, Jackson SP, Shay TL, Farnir F, Berghmans S, Snowder GD, et al. Polar overdominance at the ovine callipyge locus. Science. 1996;273:236–8.
Freking BA, Murphy SK, Wylie AA, Rhodes SJ, Keele JW, Leymaster KA, et al. Identification of the single base change causing the callipyge muscle hypertrophy phenotype, the only known example of polar overdominance in mammals. Genome Res. 2002;12:1496–506.
Rupp R, Senin P, Sarry J, Allain C, Tasca C, Ligat L, et al. A point mutation in suppressor of cytokine signalling 2 (Socs2) increases the susceptibility to inflammation of the mammary gland while associated with higher body weight and size and higher milk production in a sheep model. PLoS Genet. 2015;11:e1005629.
Zhou S, Cai B, He C, Wang Y, Ding Q, Liu J, et al. Programmable base editing of the sheep genome revealed no genome-wide off-target mutations. Front Genet. 2019;10:215.
China National Commission of Animal Genetic Resources. Animal genetic resources in China: sheep and goats. Beijing: China Agriculture Press; 2011.
Zhu C, Fan H, Yuan Z, Hu S, Ma X, Xuan J, et al. Genome-wide detection of CNVs in Chinese indigenous sheep with different types of tails using ovine high-density 600K SNP arrays. Sci Rep. 2016;6:27822.
Zhang T, Gao H, Sahana G, Zan Y, Fan H, Liu J, et al. Genome-wide association studies revealed candidate genes for tail fat deposition and body size in the Hulun Buir sheep. J Anim Breed Genet. 2019;136:362–70.
Zhao F, Deng T, Shi L, Wang W, Zhang Q, Du L, et al. Genomic scan for selection signature reveals fat deposition in Chinese indigenous sheep with extreme tail types. Animals (Basel). 2020;10:773.
Zhu C, Li N, Cheng H, Ma Y. Genome wide association study for the identification of genes associated with tail fat deposition in Chinese sheep breeds. Biol Open. 2021. https://doi.org/10.1242/bio.054932.
Tao L, Liu YF, Zhang H, Li HZ, Zhao FP, Wang FY, et al. Genome-wide association study and inbreeding depression on body size traits in Qira black sheep (Ovis aries). Anim Genet. 2021;52:560–4.
Luo R, Zhang X, Wang L, Zhang L, Li G, Zheng Z. GLIS1, a potential candidate gene affect fat deposition in sheep tail. Mol Biol Rep. 2021;48:4925–31.
Taghizadeh S, Gholizadeh M, Rahimi-Mianji G, Moradi MH, Costilla R, Moore S, et al. Genome-wide identification of copy number variation and association with fat deposition in thin and fat-tailed sheep breeds. Sci Rep. 2022;12:8834.
Moradi MH, Nejati-Javaremi A, Moradi-Shahrbabak M, Dodds KG, Brauning R, McEwan JC. Hitchhiking mapping of candidate regions associated with fat deposition in Iranian thin and fat tail sheep breeds suggests new insights into molecular aspects of fat tail selection. Animals (Basel). 2022;12:1423.
Luo R, Dai X, Zhang L, Li G, Zheng Z. Genome-wide DNA methylation patterns of muscle and tail-fat in DairyMeade sheep and Mongolian sheep. Animals (Basel). 2022;12:1399.
Wang X, Zhou G, Xu X, Geng R, Zhou J, Yang Y, et al. Transcriptome profile analysis of adipose tissues from fat and short-tailed sheep. Gene. 2014;549:252–7.
Miao X, Luo Q, Qin X, Guo Y, Zhao H. Genome-wide mRNA-seq profiling reveals predominant down-regulation of lipid metabolic processes in adipose tissues of Small Tail Han than Dorset sheep. Biochem Biophys Res Commun. 2015;467:413–20.
Kang D, Zhou G, Zhou S, Zeng J, Wang X, Jiang Y, et al. Comparative transcriptome analysis reveals potentially novel roles of Homeobox genes in adipose deposition in fat-tailed sheep. Sci Rep. 2017;7:14491.
Li B, Qiao L, An L, Wang W, Liu J, Ren Y, et al. Transcriptome analysis of adipose tissues from two fat-tailed sheep breeds reveals key genes involved in fat deposition. BMC Genomics. 2018;19:338.
Ma L, Li Z, Cai Y, Xu H, Yang R, Lan X. Genetic variants in fat- and short-tailed sheep from high-throughput RNA-sequencing data. Anim Genet. 2018;49:483–7.
Bakhtiarizadeh MR, Salehi A, Alamouti AA, Abdollahi-Arpanahi R, Salami SA. Deep transcriptome analysis using RNA-Seq suggests novel insights into molecular aspects of fat-tail metabolism in sheep. Sci Rep. 2019;9:9203.
Fan H, Hou Y, Sahana G, Gao H, Zhu C, Du L, et al. A transcriptomic study of the tail fat deposition in two types of Hulun Buir sheep according to tail size and sex. Anim Multidiscip. 2019;9:655.
Yuan Z, Xiang R, Li W, Li F, Yue X. Transcriptomic analyses revealed common tailed and perirenal adipose differentially expressed genes in four Chinese indigenous sheep breeds. Livest Sci. 2019;230:103832.
Bakhtiarizadeh MR, Alamouti AA. RNA-Seq based genetic variant discovery provides new insights into controlling fat deposition in the tail of sheep. Sci Rep. 2020;10:13525.
Fatemeh Hosseini S, Bakhtiarizadeh MR, Salehi A. Meta-analysis of RNA-Seq datasets highlights novel genes/pathways involved in fat deposition in fat-tail of sheep. Res Sq. 2022. https://doi.org/10.21203/rs.3.rs-1770948/v1.
Farhadi S, ShodjaGhias J, Hasanpur K, Mohammadi SA, Ebrahimie E. Molecular mechanisms of fat deposition: IL-6 is a hub gene in fat lipolysis, comparing thin-tailed with fat-tailed sheep breeds. Arch Anim Breed. 2021;64:53–68.
Wang X, Fang C, He H, Cao H, Liu L, Jiang L, et al. Identification of key genes in sheep fat tail evolution based on RNA-seq. Gene. 2021;781:145492.
Zhang W, Xu M, Wang J, Wang S, Wang X, Yang J, et al. Comparative transcriptome analysis of key genes and pathways activated in response to fat deposition in two sheep breeds with distinct tail phenotype. Front Genet. 2021;12:639030.
Han J, Ma S, Liang B, Bai T, Zhao Y, Ma Y, et al. Transcriptome profiling of developing ovine fat tail tissue reveals an important role for MTFP1 in regulation of adipogenesis. Front Cell Dev Biol. 2022;10:839731.
Jin M, Fei X, Li T, Lu Z, Chu M, Di R, et al. Transcriptome study digs out BMP2 involved in adipogenesis in sheep tails. BMC Genomics. 2022;23:457.
Miao X, Luo Q, Qin X, Guo Y. Genome-wide analysis of microRNAs identifies the lipid metabolism pathway to be a defining factor in adipose tissue from different sheep. Sci Rep. 2015;5:18470.
Zhou G, Wang X, Yuan C, Kang D, Xu X, Zhou J, et al. Integrating miRNA and mRNA expression profiling uncovers miRNAs underlying fat deposition in sheep. Biomed Res Int. 2017;2017:1857580.
Pan Y, Jing J, Zhao J, Jia X, Qiao L, et al. MicroRNA expression patterns in tail fat of different breeds of sheep. Livest Sci. 2018;207:7–14.
Yang G, Zhang H, Zhang S, Li Z, Gao F, Wang G, et al. Identification of the genetic basis for the large-tailed phenotypic trait in Han sheep through integrated mRNA and miRNA analysis of tail fat tissue samples. Res Sq. 2020. https://doi.org/10.21203/rs.3.rs-107294/v1.
Fei X, Jin M, Wang Y, Li T, Lu Z, Yuan Z, et al. Transcriptome reveals key microRNAs involved in fat deposition between different tail sheep breeds. PLoS One. 2022;17:e0264804.
Jin M, Fei X, Li T, Lu Z, Chu M, Di R, et al. Oar-miR-432 regulates fat differentiation and promotes the expression of BMP2 in ovine preadipocytes. Front Genet. 2022;13: 844747.
Wang Q, Cao H, Su X, Liu W. Identification of key miRNAs regulating fat metabolism based on RNA-seq from fat-tailed sheep and F2 of wild Argali. Gene. 2022;834:146660.
Meale SJ, Romao JM, He ML, Chaves AV, McAllister TA, Guan LL. Effect of diet on microRNA expression in ovine subcutaneous and visceral adipose tissues. J Anim Sci. 2014;92:3328–37.
Deng K, Ren C, Fan Y, Liu Z, Zhang G, Zhang Y, et al. miR-27a is an important adipogenesis regulator associated with differential lipid accumulation between intramuscular and subcutaneous adipose tissues of sheep. Domest Anim Endocrinol. 2020;71:106393.
Li B, Huang X, Yang C, Ge T, Zhao L, Zhang X, et al. miR-27a regulates sheep adipocyte differentiation by targeting CPT1B gene. Animals (Basel). 2021;12:28.
Jin X, Hao Z, Zhao M, Shen J, Ke N, Song Y, et al. microRNA-148a regulates the proliferation and differentiation of ovine preadipocytes by targeting PTEN. Animals (Basel). 2021;11:820.
Jin X, Wang J, Hu J, Liu X, Li S, Lu Y, et al. MicroRNA-200b regulates the proliferation and differentiation of ovine preadipocytes by targeting p27 and KLF9. Animals (Basel). 2021;11:2417.
Han F, Zhou L, Zhao L, Wang L, Liu L, Li H, et al. Identification of miRNA in sheep intramuscular fat and the role of miR-193a-5p in proliferation and differentiation of 3T3-L1. Front Genet. 2021;12:633295.
Ma L, Zhang M, Jin Y, Erdenee S, Hu L, Chen H, et al. Comparative transcriptome profiling of mRNA and lncRNA related to tail adipose tissues of sheep. Front Genet. 2018;9:365.
Bakhtiarizadeh MR, Salami SA. Identification and expression analysis of long noncoding RNAs in fat-tail of sheep breeds. G3 (Bethesda). 2019;9:1263–76.
He X, Wu R, Yun Y, Qin X, Chen L, Han Y, et al. Transcriptome analysis of messenger RNA and long noncoding RNA related to different developmental stages of tail adipose tissues of sunite sheep. Food Sci Nutr. 2021;9:5722–34.
Yuan Z, Ge L, Sun J, Zhang W, Wang S, Cao X, et al. Integrative analysis of Iso-Seq and RNA-seq data reveals transcriptome complexity and differentially expressed transcripts in sheep tail fat. PeerJ. 2021;9:e12454.
Su X-H, He H-Y, Fang C, Liu L-L, Liu W-J. Transcriptome profiling of LncRNAs in sheep tail fat deposition. Anim Biotechnol. 2021;4:1–11.
Wang J, Xu M, Wang X, Yang J, Gao L, Zhang Y, et al. Comparative proteome analysis reveals lipid metabolism-related protein networks in response to rump fat mobilization. Int J Mol Sci. 2018;19:2556.
Han J, Guo T, Yue Y, Lu Z, Liu J, Yuan C, et al. Quantitative proteomic analysis identified differentially expressed proteins with tail/rump fat deposition in Chinese thin- and fat-tailed lambs. PLoS One. 2021;16:e0246279.
Zhu C, Cheng H, Li N, Liu T, Ma Y. Isobaric tags for relative and absolute quantification-based proteomics reveals candidate proteins of fat deposition in Chinese indigenous sheep with morphologically different tails. Front Genet. 2021;12:710449.
Bakhtiarizadeh MR, Moradi-Shahrbabak M, Ebrahimie E. Underlying functional genomics of fat deposition in adipose tissue. Gene. 2013;521:122–8.
Xu X, Wei X, Yang Y, Niu W, Kou Q, Wang X, et al. PPARγ, FAS, HSL mRNA and protein expression during Tan sheep fat-tail development. Electron J Biotechnol. 2015;18:122–7.
Bahnamiri HZ, Zali A, Ganjkhanlou M, Sadeghi M, Shahrbabak HM. Regulation of lipid metabolism in adipose depots of fat-tailed and thin-tailed lambs during negative and positive energy balances. Gene. 2018;641:203–11.
ZakariapourBahnamiri H, Ganjkhanlou M, Zali A, Sadeghi M, Moradi Shahrbabak H, NehzatiPaghaleh GA. Expression of genes related to liver fatty acid metabolism in fat-tailed and thin-tailed lambs during negative and positive energy balances. J Anim Physiol Anim Nutr (Berl). 2019;103:427–35.
Zeng J, Zhou S-W, Zhao J, Jin M-H, Kang D-J, Yang Y-X, et al. Role of OXCT1 in ovine adipose and preadipocyte differentiation. Biochem Biophys Res Commun. 2019;512:779–85.
Li B, Qiao L, Yan X, Shi T, Ren D, Zhao Y, et al. mRNA expression of genes related to fat deposition during in vitro ovine adipogenesis. Can J Anim Sci. 2019;99:764–71.
Li Q, Lu Z, Jin M, Fei X, Quan K, Liu Y, et al. Verification and analysis of sheep tail type-associated PDGF-D gene polymorphisms. Animals (Basel). 2020;10:89.
Lu Z, Liu J, Han J, Yang B. Association between BMP2 functional polymorphisms and sheep tail type. Animals (Basel). 2020;10:739.
Zeng J, Zhou S, Yang Y, Du J, Kang D, Wang X, et al. Effect of dietary nutrition on tail fat deposition and evaluation of tail-related genes in fat-tailed sheep. Electron J Biotechnol. 2020;46:30–7.
Liang C, Qiao L, Han Y, Liu J, Zhang J, Liu W. Regulatory roles of SREBF1 and SREBF2 in lipid metabolism and deposition in two Chinese representative fat-tailed sheep breeds. Animals (Basel). 2020;10:1317.
Wang YL, Li X, Cao Y, Xiao C, Liu Y, Jin HG, et al. Effect of the ACAA1 gene on preadipocyte differentiation in sheep. Front Genet. 2021;12:649140.
Zhao L, Li F, Liu T, Yuan L, Zhang X, Zhang D, et al. Ovine ELOVL5 and FASN genes polymorphisms and their correlations with sheep tail fat deposition. Gene. 2022;807:145954.
Zhao Y, Zhang D, Zhang X, Li F, Xu D, Zhao L, et al. Expression features of the ovine FTO gene and association between FTO polymorphism and tail fat deposition related-traits in Hu sheep. Gene. 2022;826:146451.
Cui P, Wang W, Zhang D, Li C, Huang Y, Ma Z, et al. Identification of TRAPPC9 and BAIAP2 gene polymorphisms and their association with fat deposition-related traits in Hu sheep. Front Vet Sci. 2022;9:928375.
Kumar B, Francis SM, Suttie JM, Thompson MP. Expression of obese mRNA in genetically lean and fat selection lines of sheep. Comp Biochem Physiol B Biochem Mol Biol. 1998;120:543–8.
Bonnet M, Leroux C, Faulconnier Y, Hocquette JF, Bocquier F, Martin P, et al. Lipoprotein lipase activity and mRNA are up-regulated by refeeding in adipose tissue and cardiac muscle of sheep. J Nutr. 2000;130:749–56.
Lemor A, Mielenz M, Altmann M, von Borell E, Sauerwein H. mRNA abundance of adiponectin and its receptors, leptin and visfatin and of G-protein coupled receptor 41 in five different fat depots from sheep. J Anim Physiol Anim Nutr (Berl). 2010;94:e96-101.
Zhang J, Jing J-J, Jia X-L, Qiao L-Y, Liu J-H, Liang C, et al. mRNA expression of ovine angiopoietin-like protein 4 gene in adipose tissues. Asian-Australas J Anim Sci. 2016;29:615–23.
Dervishi E, Serrano C, Joy M, Serrano M, Rodellar C, Calvo JH. The effect of feeding system in the expression of genes related with fat metabolism in semitendinous muscle in sheep. Meat Sci. 2011;89:91–7.
Yuan Y-N, Liu W-Z, Liu J-H, Qiao L-Y, Wu J-L. Cloning and ontogenetic expression of the uncoupling protein 1 gene UCP1 in sheep. J Appl Genet. 2012;53:203–12.
Bai JL, Xu HW, Zang RX, He HJ, Cai Y, Cao X, et al. Cloning of the heart fatty acid-binding protein (H-FABP) gene and its tissue-specific expression profile in the Lanzhou fat-tailed sheep, Ovis aries. Small Rumin Res. 2013;112:114–22.
Jiao X-L, Jing J-J, Qiao L-Y, Liu J-H, Li L-A, Zhang J, et al. Ontogenetic expression of Lpin2 and Lpin3 genes and their associations with traits in two breeds of Chinese tat-tailed sheep. Asian-Australas J Anim Sci. 2016;29:333–42.
Deng K, Ren C, Liu Z, Gao X, Fan Y, Zhang G, et al. Characterization of RUNX1T1, an adipogenesis regulator in ovine preadipocyte differentiation. Int J Mol Sci. 2018;19:1300.
Zhang D, Bai J, Ma Z, Ma X, Cao X, Li F. Regulatory roles of adiponectin receptor 1 and 2 in sheep preadipocytes during adipocyte differentiation. Ital J Anim Sci. 2019;18:704–12.
Cao Y, Wang S, Liu S, Wang Y, Jin H, Ma H, et al. Effects of long-chain fatty acyl-CoA synthetase 1 on diglyceride synthesis and arachidonic acid metabolism in sheep adipocytes. Int J Mol Sci. 2020;21:2044.
Jing X, Zhou J, Degen A, Wang W, Guo Y, Kang J, et al. Comparison between Tibetan and Small-tailed Han sheep in adipocyte phenotype, lipid metabolism and energy homoeostasis regulation of adipose tissues when consuming diets of different energy levels. Br J Nutr. 2020;124:668–80.
Ma M, Wang C, Ao Y, He N, Hao F, Liang H, et al. HOXC10 promotes proliferation and attenuates lipid accumulation of sheep bone marrow mesenchymal stem cells. Mol Cell Probes. 2020;49:101491.
Mohammadi H, Shahrebabak MM, Sadeghi M. Association between single nucleotide polymorphism in the ovine DGAT1 gene and carcass traits in two Iranian sheep breeds. Anim Biotechnol. 2013;24:159–67.
Mohammadi F, Farhadi A, Rahimi-Mianji G, Gholizadeh M. Molecular genetics and bioinformatics analysis of EDG1 and AKIRIN2 genes in Iranian fat-tailed and nonfat-tailed sheep breeds. Small Rumin Res. 2016;144:263–8.
Wang S, Liu S, Yuan T, Sun X. Genetic effects of FTO gene insertion/deletion (InDel) on fat-tail measurements and growth traits in Tong sheep. Anim Biotechnol. 2021;32:229–39.
Al-Thuwaini TM, Al-Shuhaib MBS, Lepretre F, Mahdi ZA. Co-inherited novel SNPs of the LIPE gene associated with increased carcass dressing and decreased fat-tail weight in Awassi breed. Trop Anim Health Prod. 2020;52:3631–8.
Wang J, Zhang X, Wang X, Li F, Zhang D, Li X, et al. Polymorphism and expression of the HMGA1 gene and association with tail fat deposition in Hu sheep. Anim Biotechnol. 2021; In press.
Dominik S, Henshall JM, Hayes BJ. A single nucleotide polymorphism on chromosome 10 is highly predictive for the polled phenotype in Australian Merino sheep. Anim Genet. 2012;43:468–70.
Wang X, Zhou G, Li Q, Zhao D, Chen Y. Discovery of SNPs in RXFP2 related to horn types in sheep. Small Rumin Res. 2014;116:133–6.
Duijvesteijn N, Weerasinghe S, van der Werf JHJ. Production accuracies for polled and horned merino sheep using different genetic models. Proc Assoc Adv Anim Breed Genet. 2017;22:349–52.
Duijvesteijn N, Bolormaa S, Daetwyler HD, van der Werf JHJ. Genomic prediction of the polled and horned phenotypes in Merino sheep. Genet Sel Evol. 2018;50:28.
He X, Chen X, Pu Y, Guan W, Song S, Zhao Q, et al. iTRAQ-based quantitative proteomic analysis reveals key pathways responsible for scurs in sheep (Ovis aries). J Integr Agric. 2018;17:1843–51.
Gratten J, Pilkington JG, Brown EA, Clutton-Brock TH, Pemberton JM, Slate J. Selection and microevolution of coat pattern are cryptic in a wild population of sheep. Mol Ecol. 2012;21:2977–90.
Gratten J, Pilkington JG, Brown EA, Beraldi D, Pemberton JM, Slate J. The genetic basis of recessive self-colour pattern in a wild sheep population. Heredity (Edinb). 2010;104:206–14.
Fontanesi L, Dall’Olio S, Beretti F, Portolano B, Russo V. Coat colours in the Massese sheep breed are associated with mutations in the agouti signalling protein (ASIP) and melanocortin 1 receptor (MC1R) genes. Animal. 2011;5:8–17.
Fontanesi L, Rustempašić A, Brka M, Russo V. Analysis of polymorphisms in the agouti signalling protein (ASIP) and melanocortin 1 receptor (MC1R) genes and association with coat colours in two Pramenka sheep types. Small Rumin Res. 2012;105:89–96.
Hepp D, Gonçalves GL, Moreira GRP, de Freitas TRO. Epistatic interaction of the melanocortin 1 receptor and agouti signaling protein genes modulates wool color in the Brazilian Creole sheep. J Hered. 2016;107:544–52.
Rochus CM, Westberg Sunesson K, Jonas E, Mikko S, Johansson AM. Mutations in ASIP and MC1R: dominant black and recessive black alleles segregate in native Swedish sheep populations. Anim Genet. 2019;50:712–7.
Fontanesi L, Beretti F, Riggio V, Dall’Olio S, Calascibetta D, Russo V, et al. Sequence characterization of the melanocortin 1 receptor (MC1R) gene in sheep with different coat colours and identification of the putative e allele at the ovine Extension locus. Small Rumin Res. 2010;91:200–7.
Yang G-L, Fu D-L, Lang X, Wang Y-T, Cheng S-R, Fang S-L, et al. Mutations in MC1R gene determine black coat color phenotype in Chinese sheep. Sci World J. 2013;2013:675382.
Mahmoud AH, Mashaly AM, Rady AM, Al-Anazi KM, Saleh AA. Allelic variation of melanocortin-1 receptor locus in Saudi indigenous sheep exhibiting different color coats. Asian-Australas J Anim Sci. 2016;30:154–9.
Stamatis C, Giannoulis T, Galliopoulou E, Billinis C, Mamuris Z. Genetic analysis of melanocortin 1 receptor gene in endangered Greek sheep breeds. Small Rumin Res. 2017;157:54–6.
García-Gámez E, Reverter A, Whan V, McWilliam SM, Arranz JJ, International Sheep Genomics Consortium ISG, et al. Using regulatory and epistatic networks to extend the findings of a genome scan: Identifying the gene drivers of pigmentation in merino sheep. PLoS One. 2011;6:e21158.
Mastrangelo S, Criscione A, Sottile G, Portolano B, Marletta D, Bordonaro S. Genome-wide analysis identifies potentially causative genes explaininga the phenotypic variability in Pinzirita sheep. Anim Genet. 2019;50:189–90.
Posbergh CJ, Staiger EA, Huson HJ. A stop-gain mutation within MLPH is responsible for the lilac dilution observed in Jacob sheep. Genes (Basel). 2020;11:618.
Zsolnai A, Egerszegi I, Rózsa L, Anton I. Genetic status of lowland-type Racka sheep colour variants. Animal. 2021;15:100080.
Naya H, Urioste JI, Chang Y-M, Rodrigues-Motta M, Kremer R, Gianola D. A comparison between Poisson and zero-inflated Poisson regression models with an application to number of black spots in Corriedale sheep. Genet Sel Evol. 2008;40:379–94.
Peñagaricano F, Urioste JI, Naya H, De Los CG, Gianola D. Assessment of Poisson, probit and linear models for genetic analysis of presence and number of black spots in Corriedale sheep. J Anim Breed Genet. 2011;128:105–13.
Peñagaricano F, Zorrilla P, Naya H, Robello C, Urioste JI. Gene expression analysis identifies new candidate genes associated with the development of black skin spots in Corriedale sheep. J Appl Genet. 2012;53:99–106.
Fan R, Xie J, Bai J, Wang H, Tian X, Bai R, et al. Skin transcriptome profiles associated with coat color in sheep. BMC Genomics. 2013;14:389.
Yao L, Bao A, Hong W, Hou C, Zhang Z, Liang X, et al. Transcriptome profiling analysis reveals key genes of different coat color in sheep skin. PeerJ. 2019;7:e8077.
Shi X, Wu J, Lang X, Wang C, Bai Y, Riley DG, et al. Comparative transcriptome and histological analyses provide insights into the skin pigmentation in Minxian black fur sheep (Ovis aries). PeerJ. 2021;9:e11122.
Han J, Yang M, Guo T, Yue Y, Liu J, Niu C, et al. Molecular characterization of two candidate genes associated with coat color in Tibetan sheep (Ovis arise). J Integr Agric Elsevier. 2015;14:1390–7.
Wang H, Xue L, Li Y, Zhao B, Chen T, Liu Y, et al. Distribution and expression of SLC45A2 in the skin of sheep with different coat colors. Folia Histochem Cytobiol. 2016;54:143–50.