Skip to content


Genetics Selection Evolution

Open Access

Unveiling genomic regions that underlie differences between Afec-Assaf sheep and its parental Awassi breed

  • Eyal Seroussi1,
  • Alexander Rosov1,
  • Andrey Shirak1,
  • Alon Lam1 and
  • Elisha Gootwine1Email author
Genetics Selection Evolution201749:19

Received: 20 June 2016

Accepted: 6 February 2017

Published: 10 February 2017



Sheep production in Israel has improved by crossing the fat-tailed local Awassi breed with the East Friesian and later, with the Booroola Merino breed, which led to the formation of the highly prolific Afec-Assaf strain. This strain differs from its parental Awassi breed in morphological traits such as tail and horn size, coat pigmentation and wool characteristics, as well as in production, reproductive and health traits. To identify major genes associated with the formation of the Afec-Assaf strain, we genotyped 41 Awassi and 141 Afec-Assaf sheep using the Illumina Ovine SNP50 BeadChip array, and analyzed the results with PLINK and EMMAX software. The detected variable genomic regions that differed between Awassi and Afec-Assaf sheep (variable genomic regions; VGR) were compared to selection signatures that were reported in 48 published genome-wide association studies in sheep. Because the Afec-Assaf strain, but not the Awassi breed, carries the Booroola mutation, association analysis of BMPR1B used as the test gene was performed to evaluate the ability of this study to identify a VGR that includes such a major gene.


Of the 20 detected VGR, 12 were novel to this study. A ~7-Mb VGR was identified on Ovies aries chromosome OAR6 where the Booroola mutation is located. Similar to other studies, the most significant VGR was detected on OAR10, in a region that contains candidate genes affecting horn type (RXFP2), climate adaptation (ALOX5AP), fiber diameter (KATNAl1), coat pigmentation (FRY) and genes associated with fat distribution. The VGR on OAR2 included BNC2, which is also involved in controlling coat pigmentation in sheep. Six other VGR contained genes that were shown to be involved in coat pigmentation by analyzing their mammalian orthologues. Genes associated with fat distribution in humans, including GRB14 and COBLL1, were located in additional VGR. Sequencing DNA from Awassi and Afec-Assaf individuals revealed non-synonymous mutations in some of these candidate genes.


Our results highlight VGR that differentiate the Awassi breed from the Afec-Assaf strain, some of which may include genes that confer an advantage to Afec-Assaf and Assaf over Awassi sheep with respect to intensive sheep production under Mediterranean conditions.


Sheep BreedMerino SheepSignificant Single Nucleotide PolymorphismCoat PigmentationVariable Genomic Region


In Israel, within-breed selection, crossbreeding and gene introgression have contributed to the transition of the sheep industry from traditional extensive production using the native fat-tailed local Awassi breed to highly intensive production with the Assaf and Afec-Assaf sheep [1]. Selection within the local Awassi sheep that began in the 1930s led to the formation of the improved Awassi dairy strain [2]. The Assaf dairy breed was formed about 30 years later by crossing the improved Awassi with the European East Friesian (EF) breed [3]. In 1986, breeding of the highly prolific Afec-Assaf strain was initiated by introgressing the Booroola mutation (B) at the BMPR1B gene into the Assaf breed by crossing Assaf ewes with Booroola Merino rams [46]. The high prolificacy of the Afec-Assaf strain is due to the presence of the B allele, which is inherited in an almost completely dominant mode [7], since BB and B+ ewes have similar prolificacy. However, since the perinatal lamb survival rate is higher in B+ than in BB ewes, B+ is the recommended genotype for females in commercial flocks [7]. The B allele segregates in the Afec-Assaf population, thus genotyping lambs for the Booroola mutation and selecting only the B+ ewe lambs as replacements is a common practice in Afec-Assaf flocks to maintain high prolificacy in the long term. Afec-Assaf sheep differ widely from the Awassi sheep in a number of phenotypic characteristics, such as the presence of horns, tail shape, coat color pigmentation and wool quality, as well as seasonality, prolificacy, and milk production [8, 9] (Table 1).
Table 1

Phenotypic, production and reproductive traits of local Awassi, improved Awassi, Assaf and Afec-Assaf sheep


Local Awassi

Improved Awassi



Body size





Tail shape

Fat tail

Fat tail

Partial fat tail

Partial fat tail

Horn shape

Most rams and ewes are horned or carry scurs

Most males and females are polled, but some carry short knobs or scurs

Wool type





Coat pigmentation

Most sheep have brown face and white fleece

Most sheep have white face and white fleece




Moderately seasonal

Moderately seasonal

Adaptation to Middle-Eastern climate





Milk production





Prolificacy (lambs born/ewe lambing)





aSmall: ram body weight ~70 kg; large: ram body weight ~110 kg

In sheep, genome-wide association studies (GWAS) have contributed to better understanding of their domestication [10, 11], as well as to the identification of specific genomic regions that are associated with differences in various traits between breeds (see Additional file 1: Table S1). GWAS have made it possible to identify candidate genes that underlie various phenotypic differences, for instance in Texel sheep RXFP2 for horn type [12] and disease-resistance genes such as PITX3 and DMP1 which are involved in microphthalmia [13] and in Corriedale sheep for rickets [14].

Cryptic relatedness, which occurs when there is unknown kinship within the sample, and population stratification due to random genetic drift in the sample’s subpopulations are two major confounding effects that cause spurious associations in GWAS analyses [15]. Hence, incorporating a known genetic marker that is located within a major gene that segregates within the tested population into the GWAS analysis may assist in verifying the power of the study for detecting selection signatures. This has been the case in several sheep GWAS, in which the MC1R gene that affects coat pigmentation [16] was used as the ‘test’ gene, and polymorphism at the GDF8 gene associated with muscle hypertrophy [17] was used to establish criteria to detect selection signatures.

By using the ovine single nucleotide polymorphism (SNP) 50 BeadChip array (Illumina Inc., San Diego, CA) and BMPR1B as a test gene, the aims of our study were to: (1) compare the Awassi and Afec-Assaf genomes, searching for variable genomic regions (VGR) that differ between the two breeds; and (2) link these VGR to genes and selection signatures that were previously described in sheep GWAS.


Ethics statement

Experimental protocols were approved (Approval No. IL 415-12) by the Volcani Center Institutional Animal Care and Use Committee.


Local Awassi sheep are raised by Bedouin growers in small unconnected flocks that are maintained under extensive traditional management with no records. The Awassi cohort consisted of 24 local Awassi rams that included all of the rams from five flocks, and 17 improved Awassi ewes from the Ein-Harod flock. Since the degree of relationship between animals from the same flock is unknown, the effective number of local Awassi individuals might be smaller than sampled. However, this was not the case, since a genomic relationship analysis (see Additional file 2: Table S2) showed no strong kinship between individuals within or among flocks. The Afec-Assaf cohort consisted of 22 males and 119 females from the experimental flock of the Volcani Center at Bet Dagan and two commercial farms. Genotyping indicated that 35, 95 and 11 animals of the Afec-Assaf cohort were homozygous BB, heterozygous B+ or non-carrier ++ for the Booroola mutation, respectively. As expected [7], the average prolificacy of the BB and B+ females (based on 3 to 9 parity records per animal) was high, i.e. 3.1 and 3.0 lambs born/lambing, respectively. No lambing records were available for non-carrier ++ females.

Genotyping using the ovine 50 K beadchip

DNA was extracted from blood samples using the DNeasy Blood & Tissue Kit (Qiagen, Valencia, CA) according to the manufacturer’s instructions. DNA samples were genotyped by the Ovine SNP50 BeadChip assay using standard procedures ( All the Awassi and some of the Afec-Assaf samples were genotyped as part of the International Sheep Genomic Consortium initiative ( Other Afec-Assaf samples were genotyped at the Centre for Reproduction and Genomics (CRG), Invermay, New Zealand.

Of the 59,454 SNPs on the ovine SNP50 BeadChip, 46,563 SNPs that were randomly distributed across the genome passed the quality-control filter criteria (call rate higher than 99%, genotyping frequency higher than 95%, minor allele frequency lower than 0.05 and Hardy–Weinberg equilibrium P > 0.001), and were used for further analyses.

Genotyping at the RXFP2 gene

Phenotypic variation in horn appearance includes the presence of normal horns (in females they are smaller than in males), deformed horns (scurs), short knobs at the site of horn growth, and a polled-non-horned growth phenotype, which may include a concave depression in the skull bone at the horn site [18]. DNA from Awassi and Afec-Assaf sheep with different horn phenotypes was extracted from blood samples or buccal swabs using standard DNA-extraction protocols. Primers, PCR conditions and separation of PCR products on an agarose gel to genotype animals that carry a 1833-bp genomic insertion located in the 3′-UTR of RXFP2, were performed according to Wiedemar and Drogemuller [19]. Association between the genotype at the RXFP2 gene and horn appearance was tested by Chi square analysis.

DNA sequencing

To obtain preliminary information on variation of the coding regions of genes that are included in VGR, DNA of one Afec-Assaf individual and one Awassi individual, selected randomly, was sequenced using the Illumina 100-bp paired-end technology (HiSeq 2000, Illumina) with one or two lanes, respectively. The sequences were deposited in the European Nucleotide Archive (ENA) under accession number PRJEB12018. Achieved coverage was 11- and 26-fold for the Afec-Assaf and Awassi individuals, respectively.

Identification of candidate polymorphisms

Exon sequences of the candidate genes were downloaded from NCBI ( and genomic reads were mapped to these templates using GAP5 software [20]. BWA options for this mapping were set to bam bwasw -t 8 -T 30 [21]. Genetic markers were identified and called by comparing the consensus contig sequences of the Awassi and Assaf individuals with the GAP5 contig comparator and the heterozygosity search option (value = 10).

Scoring amino acid substitutions

Amino acid substitutions were evaluated using PROVEAN analysis v1.1 [22] which scores the effect of a variation in protein sequence on the function of that protein.

Statistical analysis

Various methods have been applied to verify selection signatures in cattle [23] as well as in sheep. To minimize the risk of detecting spurious associations in our GWAS analysis, we applied two tools: PLINK, which is suitable for calculating either F st values or probability values for genetic association test and EMMAX that calculates these probability values using a relationship matrix that accounts for relatedness within the sample. Association between SNPs and breed (Awassi or Afec-Assaf) were evaluated using PLINK v1.90b3.38 64-bit [24]. Probability (Fisher’s exact test of significance) and F st values for differences in allele frequency were estimated with the following options on the command line: –sheep –hwe 0.001 –geno 0.1 –maf 0.05 –fst –assoc fisher. F st estimates were obtained only for the autosomes and were computed using the method introduced by Weir and Cockerham. The Efficient Mixed-Model Association eXpedited (EMMAX) program [25] was used to identify genomic regions on autosomes and on the X chromosome. For the EMMAX analysis, the following options were used on the command lines: emmax-kin -v -h -s -d 10; emmax -v -d 10, for the creation of the identity by state (IBS) relationship matrix and for the association test, respectively.

Genomic inflation factors were equal to 9.35 and 1.34 for the unadjusted PLINK Fisher and EMMAX analyses, respectively. Using each of the programs, SNPs were ranked according to their association with one of the breeds by chance. The top ~0.1% (n = 50) significant SNPs and their adjacent neighboring SNPs, which ranked in the top 0.5% and were found within 1.5 Mb upstream or downstream of the top marker, were used to define VGR. Neighboring genomic regions were combined into one region if they overlapped. The Bonferroni-corrected significance threshold of 1.07 × 10−6 (0.05/46,563) was used to correct for multiple comparisons and the threshold for significance was set at P < 0.05. The position of regions was determined according to version 3.1 of the sheep genome assembly.


Candidate regions were examined using the sheep genome browsers Build v3.1 [2628]. We used the following databases: NCBI [29], GeneCards® [30], and SheepQTL [31], as well as literature searches, to obtain gene annotations and functions.

Results and discussion

PLINK analysis showed a correlation of 0.98 between the values of F st and the minus logarithm of the probability according to Fisher’s exact test (data not shown). Thus, since both parameters yielded similar estimates for allele frequency differences, we chose to consider only the F st values in the further analyses. PLINK and EMMAX-based GWAS revealed 36 and 32 VGR, respectively (see Additional file 3: Table S3). F st values for the most significant SNPs in each region in the PLINK analysis ranged from 0.46 to 0.80, and the corresponding P values in the EMMAX analysis ranged from <10−5 to <10−15. The 20 VGR that were detected by both PLINK and EMMAX analyses were further analyzed (Table 2 and Fig. 1).
Table 2

Variable genomic regions detected in Awassi and Afec-Assaf sheep according to both PLINK and EMMAX analyses

Variable genomic region number

OAR number

Region position (Mb)

PLINK - Peak SNP position (Mb)

EMMAX -Peak SNP position (Mb)


EMMAX P value




















































































































Fig. 1

Manhattan plot of GWAS comparing Awassi and Afec-Assaf sheep

Variable genomic regions on OAR6 (between 24.97 and 32.04 Mb)

As expected, a VGR (#10) that spanned the BMPR1B gene (between positions 29.36 and 29.45 Mb) was identified on OAR6 (Table 2), with the most significant SNP, s24937.1, being located about 3.0 Mb upstream of the gene. None of the four SNPs of the Illumina Ovine SNP50 BeadChip that are located within non-coding regions of BMPR1B namely: OAR6_33210385.1, OAR6_33229928.1, OAR6_33291928.1 or s21306, reached significance (P > 0.05). DNA sequencing showed that the Afec-Assaf individual was heterozygous for the Booroola mutation, while the Awassi individual did not carry the mutation.

A similar level of precision in the detection of a point mutation was obtained in a GWAS applied in Texel sheep to map the causative mutation for microphthalmia, which revealed a region of about 2.6 Mb [13]. In the case of the autosomal recessive lethal disorder brachygnathia, cardiomegaly and renal hypoplasia syndrome in polled Merino sheep, SNPs that are associated with this syndrome are spread across 1.1 Mb [32]. However, the precision that was achieved in a GWAS for detecting the causative mutation for dominant pigmentation in the MC1R gene was higher [16], with only one significant SNP that was 25 kb away from MC1R, being identified for the Manchega sheep population, in which this trait segregates.

The VGR on OAR6 that includes the BMPR1B gene also emerged in a GWAS that searched for genomic differences between dairy and non-dairy sheep breeds [17], and in another study that focused on quantitative trait loci (QTL) for carcass traits [33]. Since the Booroola mutation was introgressed into the Assaf breed from the small-size non-dairy Booroola Merino breed, the genomic difference that was found between the Awassi and Assaf breeds at around 29 Mb on OAR6 might also be related to the difference in allelic frequency of genes that control these traits. One relevant gene was ABCG2, which has a major effect on milk production, and is mapped to OAR6 at 36.5 Mb, about 3.0 Mb downstream to the 3′ end of that VGR [17].

Variable genomic regions on OAR10 that span 30 Mb

A major VGR (#15) was found on OAR10 (between positions 27.65 and 31.64 Mb, Table 2). Significant selection signatures in this chromosomal region have been detected in several GWAS in which sheep breeds were compared for horn phenotype, climate adaptation, fiber diameter, coat pigmentation or fat distribution [1012, 18, 3439].

Horn phenotype

Awassi rams are mostly horned, as well as some of the ewes, whereas most Assaf sheep are polled, similar to the situation observed in the EF parental breed. Indeed, the RXFP2 gene (on OAR10 between 29.45 and 29.50 Mb), which is involved in horn inheritance [12], was identified in the middle of this OAR10 VGR (Table 2). To further corroborate this, we genotyped 61 horned and non-horned Awassi and Afec-Assaf sheep (see Additional file 4: Figure S1), to identify those that carried a 1833-bp genomic insertion located in the 3′-UTR of RXFP2. It has been shown that this insertion is linked to horn inheritance in Swiss sheep breeds [19]. Our results (see Additional file 5: Table S4) show that in Awassi and Afec-Assaf sheep, there is also an association between the RXFP2 polymorphism in the 3′-UTR and the hornness-polledness phenotype (P < 0.0001), with absence of the insertion being associated with the horned phenotype. However, the range of phenotypes observed in the group of heterozygous animals indicated that other genetic modifications may also underlie the horned phenotype. A recent study of 38 sheep breeds [40] reached similar conclusions i.e. that the 1.78-kb insertion in the 3′-UTR of RXFP2 cannot be considered as the only cause of polledness in sheep.

Climate adaptation

The EF and Awassi breeds—the parental breeds of the Afec-Assaf—differ in their adaptability to the hot and humid Mediterranean climate. Susceptibility to pneumonia in the summer is a major death threat in Israel for pure EF but not Awassi lambs [41]. A GWAS that investigated environmental adaptive selection in sheep [36] suggested ALOX5AP (on OAR10 between 30.36 and 30.38 Mb) as a candidate gene for climate adaptation. Located in the distal part of VGR #15 on OAR10, this gene encodes a protein that is required for the synthesis of leukotrienes, which are biologically active lipid mediators involved in various types of inflammatory responses. In humans, polymorphism in ALOX5AP has been associated with lung function [42, 43].

Fiber diameter

In a GWAS that included Chinese sheep breeds, differences in the VGR on OAR10 were associated with differences in the coefficient of variation of fiber diameter [37]. Indeed, the Awassi and Assaf breeds differ in their variation in fiber diameter, with the Awassi’s fibers being more uniform [44]. The KATNAL1 gene (on OAR10 between 30.77 and 30.78 Mb) was the closest gene mapped to the most significant SNP (OAR10: 30,700,856).

Coat pigmentation

As indicated in Table 1, the coat color of most Afec-Assaf sheep is white, while most Awassi sheep have reddish-brown pigmentation on their head, neck and legs. The mode of inheritance of coat pigmentation pattern in sheep is complex, as already shown for the piebald phenotype in Merino sheep [45]. Of the 11 genes that have been suggested as candidates that affect sheep coat pigmentation (see Additional file 6: Table S5), FRY (on OAR10 between 28.98 and 29.19 Mb) was located in the middle of VGR #15 on OAR10. FRY is a key candidate gene for the piebald phenotype in Merino sheep [45] and for coat color differences between Rambouillet (white coat color) and Suffolk (black head and legs) sheep [35], and is associated with the black spot phenotype in Valley-type Tibetan sheep [38]. Thus, polymorphism associated with the FRY gene may contribute to the differences in coat color pigmentation distribution between the Awassi and Afec-Assaf sheep.

Sequence analysis

Comparing the coding sequences of the three candidate genes, namely: ALOX5AP, KATNAL1 and FRY, of Awassi and Assaf individuals revealed two synonymous mutations, no sequence variation and three conservative amino-acid substitutions: p.V383I, p.V1207I and p.A1722S, respectively (see Additional file 7).

Fat distribution

While the Awassi sheep have a relatively large fat tail, which usually hangs down to the hock, most Assaf and Afec-Assaf sheep carry a reduced-size fat tail which varies in size (see Additional file 8: Figure S2). Archeological evidence has documented the presence of fat-tailed sheep in the Fertile Crescent as early as 4000 years BC. In this region, sheep were the major source of animal fat, which was consumed as a delicacy and used for industrial purposes. Selection for fat tail was driven, among other reasons, by the convenience of harvesting the fat tissue from the carcasses. A study of body and carcass conformation of Awassi and Merino–Awassi crossbred sheep [46] demonstrated that the reduction in fat-tail size in Awassi crosses is associated with a change in body fat distribution: carrying a fat tail was associated with a decrease in the relative amount of body fat (omental, mesenteric, kidney, channel and scrotal fat), as well as a decrease in the relative amount of subcutaneous and intramuscular fat.

GWAS comparing fat- and thin-tailed sheep [38, 39, 47, 48] or studies on differential gene expression in such breeds [4951] revealed several chromosomal regions that are associated with the presence of a fat tail (see Additional file 9: Table S6). Only two of these regions, i.e. regions numbered 14 and 15, both on OAR10, overlapped with VGR that distinguish Awassi from Afec-Assaf sheep. These two regions were identified in a study that compared Cyprus fat-tailed sheep, a breed related to the Awassi breed, and Laticauda sheep, to several Italian thin-tailed breeds [39].

Taken together, the VGR on OAR10 that spans a region between 27.65 and 31.64 Mb and detected in the current study, has been associated in various sheep GWAS with several phenotypic traits that differ between Awassi and Afec-Assaf sheep.

Additional variable genomic regions associated with differences in coat pigmentation

Two of the coat pigmentation-related genes known in sheep: BNC2 and TYRP1 (see Additional file 6: Table S5) were localized, respectively, within or close to VGR #6 that spans a region between 83.11 and 84.53 Mb on OAR2 (Table 2). BNC2 (on OAR2 between 84.29 and 84.69 Mb) was identified as a candidate gene for coat pigmentation in Italian breeds [11] with the most differentiated breed being the Comisana, which resembles the Awassi breed with a white body and brick-red face. The other breeds that contrasted with the Comisana were the Altamurana and Leccese white-coat breeds, and the Sardinian ancestral black sheep. Notably, the TYRP1 gene (on OAR2 between 80.60 and 80.62 Mb), which is associated with skin and fiber pigmentation differences between Awassi and Merino sheep [52], and with pigmentation polymorphism in Soay sheep [53] and in Finnsheep [54], was mapped closely to VGR #6 (on OAR2 between 83.11 and 84.53). Comparing sequences of Awassi and Afec-Assaf BNC2- and TYRP1-coding regions (see Additional file 7) yielded one neutral substitution (p.W70R) and one conservative replacement, respectively.

To further search for possible candidate genes associated with coat pigmentation differences between Awassi and Afec-Assaf sheep, we compiled a list of 172 genes known to affect coat pigmentation in humans [5557], cattle [58] and other mammals [59] (see Additional file 10: Table S7). Nine of these genes namely: NTRK1, RPL24, SEMA4A, KRT1, KRT4, KRT75, BMPR1B, MAB21L2, and EDARADD were located within, or close to one of the VGR that were identified on six ovine chromosomes.

Taken together, our results indicate that two of the coat pigmentation-related genes known in sheep: FRY and BNC2, and probably TYRP1, as well as some other genes have possible roles in the differences in coat pigmentation between Awassi and Afec-Assaf sheep.

Additional variable genomic regions associated with differences in fat distribution

Fat distribution between subcutaneous and visceral adipose tissue in humans is genetically controlled [60]. After annotating 107 genes associated with fat distribution in humans [6064] (see Additional file 11: Table S8), we found that the DCST2 gene was annotated close to VGR #3 on OAR1 (between 103.40 and 103.45 Mb), and that two genes, GRB14 and COBLL1, were both located close to VGR #7 on OAR2 (between 145.65 and 148.41 Mb). Only a few synonymous and conservative substitutions were identified in DCST2 and GRB14, respectively, when the corresponding coding sequences of Awassi and Afec-Assaf sheep were compared (see Additional file 7). However, in the case of COBLL1, four non-synonymous substitutions were found, one of which, p.I448S (ref XP_012011424.1), was deleterious according to the PROVEN analysis [22]. If COBLL1 is indeed involved in the control of differences in fat deposition in sheep, fat-tailed sheep could serve as a natural large animal model for studying the possible adverse metabolic and physiological outcomes of body fat distribution.

Variable genomic regions not associated with a specific trait

Our analysis revealed additional VGR that were previously identified in other GWAS. However, these were not associated with a specific production, reproductive or health trait. VGR #13 on OAR8 (between 62.62 and 62.83 Mb) was detected in a GWAS that included multiple worldwide breeds [10], while VGR #16 on OAR10 (between 33.75 and 36.24 Mb) was detected in a GWAS that comprised five American breeds [35]. VGR #19 on OAR25 (between 6.65 and 8.66 Mb) was characterized as an “ancestral signature of selection” that differentiates European from Mediterranean breeds early in domestication [11]. It is worth noting that the EF and Awassi breeds, which are the parental breeds of the Afec-Assaf strain, belong to these two groups, respectively. This VGR on OAR25 was also identified in several studies in which different breeds were compared [10, 23, 38]. Finally, VGR #20 on OAR26 (36.58 Mb) was identified in a GWAS that analyzed the highly prolific Afec-Assaf ewes with different perinatal lamb viability rates [65].

Twelve of the 20 VGR identified in the current study, namely VGR #1, 2, 3, 4, 5, 7, 8, 9, 11, 12, 17 and 18 (Table 2) were not described in any of the 48 GWAS sheep studies listed in Table S1 (see Additional file 1: Table S1). These VGR, which may have arisen due to either selection or random genomic drift, may contain chromosomal regions that contribute to the genetic differences between the Awassi breed and its derived Afec-Assaf strain.


Using the Illumina Ovine SNP50 BeadChip array, we revealed 20 VGR that are associated with differences in morphological, production and reproductive traits between the Afec-Assaf strain and the Awassi breed, a hardy Middle Eastern fat-tailed breed. Focusing on BMPR1B as a test gene allowed us to evaluate the ability of our study to identify VGR that harbor major genes. We compared our results to those of 48 GWAS in sheep, and GWAS in other domestic animals and in humans. Thus, we were able to highlight candidate genes that are involved in differences between Awassi and Afec-Assaf sheep regarding morphological traits, such as horn appearance, presence of a fat tail and coat pigmentation. In most cases, a comparison between the exon sequences of the candidate genes from Awassi and Afec-Assaf individuals revealed synonymous or neutral variations.

It is of interest that 60% of the VGR identified in the current study were not described in previous sheep GWAS that included a variety of sheep breeds. Bearing in mind that VGR may also arise from random drift, selection for genes that are included in these VGR may have further improved the Assaf breed and contributed to establish it as the major dairy breed in Israel and Spain, where it outperforms the local dairy breeds [66].


Authors’ contributions

EG and ES designed the project. AR was responsible for breeding and data collection for Afec-Assaf sheep. AL, AR and EG were responsible for collecting the DNA samples. ES was responsible for all bioinformatic analyses. AS carried out the genotyping. EG drafted the manuscript, and ES and AS contributed to the manuscript in its final version. All authors read and approved the final manuscript.


Contribution from the Agricultural Research Organization, Institute of Animal Science, Bet-Dagan, Israel, is acknowledged.

Competing interests

The authors declare that they have no competing interests.

Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Authors’ Affiliations

Institute of Animal Science, ARO, The Volcani Center, Rishon LeZion, Israel


  1. Gootwine E. Mini review: breeding Awassi and Assaf sheep for diverse management conditions. Trop Anim Health Prod. 2011;43:1289–96.View ArticlePubMedGoogle Scholar
  2. Epstein H. The Awassi sheep with special reference to the improved dairy type. FAO Animal production and health paper 57. Rome: Food and Agriculture Organization of the United Nations; 1985.Google Scholar
  3. Gootwine E, Goot H. Lamb and milk production of Awassi and East-Friesian sheep and their crosses under Mediterranean environment. Small Rumin Res. 1966;20:255–60.View ArticleGoogle Scholar
  4. Mulsant P, Lecerf F, Fabre S, Schibler L, Monget P, Lanneluc I, et al. Mutation in bone morphogenetic protein receptor-IB is associated with increased ovulation rate in Booroola Merino ewes. Proc Natl Acad Sci USA. 2001;98:5104–9.View ArticlePubMedPubMed CentralGoogle Scholar
  5. Souza CJ, MacDougall C, Campbell BK, McNeilly AS, Baird DT. The Booroola (FecB) phenotype is associated with a mutation in the bone morphogenetic receptor type 1 B (BMPR1B) gene. J Endocrinol. 2001;169:R1–6.View ArticlePubMedGoogle Scholar
  6. Wilson T, Wu XY, Juengel JL, Ross IK, Lumsden JM, Lord EA, et al. Highly prolific Booroola sheep have a mutation in the intracellular kinase domain of bone morphogenetic protein IB receptor (ALK-6) that is expressed in both oocytes and granulosa cells. Biol Reprod. 2001;64:1225–35.View ArticlePubMedGoogle Scholar
  7. Gootwine E, Reicher S, Rosov A. Prolificacy and lamb survival at birth in Awassi and Assaf sheep carrying the FecB (Booroola) mutation. Anim Reprod Sci. 2008;108:402–11.View ArticlePubMedGoogle Scholar
  8. Gootwine E, Pollott GE. Factors affecting milk production in Improved Awassi dairy ewes. Anim Sci. 2000;71:607–15.View ArticleGoogle Scholar
  9. Pollott GE, Gootwine E. Reproductive performance and milk production of Assaf sheep in an intensive management system. J Dairy Sci. 2004;87:3690–703.View ArticlePubMedGoogle Scholar
  10. 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.View ArticlePubMedPubMed CentralGoogle Scholar
  11. Fariello MI, Servin B, Tosser-Klopp G, Rupp R, Moreno C, San Cristobal M, et al. Selection signatures in worldwide sheep populations. PLoS One. 2014;9:e103813.View ArticlePubMedPubMed CentralGoogle Scholar
  12. Johnston SE, McEwan JC, Pickering NK, Kijas JW, Beraldi D, Pilkington JG, 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.View ArticlePubMedGoogle Scholar
  13. Becker D, Tetens J, Brunner A, Burstel D, Ganter M, Kijas J, et al. Microphthalmia in Texel sheep is associated with a missense mutation in the paired-like homeodomain 3 (PITX3) gene. PLoS One. 2010;5:e8689.View ArticlePubMedPubMed CentralGoogle Scholar
  14. Zhao X, Dittmer KE, Blair HT, Thompson KG, Rothschild MF, Garrick DJ. A novel nonsense mutation in the DMP1 gene identified by a genome-wide association study is responsible for inherited rickets in Corriedale sheep. PLoS One. 2011;6:e21739.View ArticlePubMedPubMed CentralGoogle Scholar
  15. Astle W, Balding DJ. Population structure and criptic relatedness in genetic association studies. Stat Sci. 2009;24:451–71.View ArticleGoogle Scholar
  16. 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.View ArticlePubMedGoogle Scholar
  17. Gutierrez-Gil B, Arranz JJ, Pong-Wong R, Garcia-Gamez E, Kijas J, Wiener P. Application of selection mapping to identify genomic regions associated with dairy production in sheep. PLoS One. 2014;9:e94623.View ArticlePubMedPubMed CentralGoogle Scholar
  18. 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.View ArticlePubMedGoogle Scholar
  19. Wiedemar N, Drogemuller C. A 1.8-kb insertion in the 3′-UTR of RXFP2 is associated with polledness in sheep. Anim Genet. 2015;46:457–61.View ArticlePubMedPubMed CentralGoogle Scholar
  20. Bonfield JK, Whitwham A. Gap5–editing the billion fragment sequence assembly. Bioinformatics. 2010;26:1699–703.View ArticlePubMedPubMed CentralGoogle Scholar
  21. Li H, Durbin R. Fast and accurate long-read alignment with Burrows–Wheeler transform. Bioinformatics. 2010;26:589–95.View ArticlePubMedPubMed CentralGoogle Scholar
  22. Choi Y, Sims GE, Murphy S, Miller JR, Chan AP. Predicting the functional effect of amino acid substitutions and indels. PLoS One. 2012;7:e46688.View ArticlePubMedPubMed CentralGoogle Scholar
  23. Randhawa IA, Khatkar MS, Thomson PC, Raadsma HW. A meta-assembly of selection signatures in cattle. PLoS One. 2016;11:e0153013.View ArticlePubMedPubMed CentralGoogle Scholar
  24. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007;81:559–75.View ArticlePubMedPubMed CentralGoogle Scholar
  25. Kang M, Sul JH, Service SK, Zaitlen NA, Kong SY, Freimer NB, et al. Variance component model to account for sample structure in genome-wide association studies. Nat Genet. 2010;42:348–54.View ArticlePubMedPubMed CentralGoogle Scholar
  27. Accessed 14 June 2016.
  28. Accessed 14 June 2016.
  29. Accessed 14 June 2016.
  30. 14 June 2016.
  31. 14 June 2016.
  32. Shariflou MR, Wade CM, Kijas J, McCulloch R, Windsor PA, Tammen I, et al. Brachygnathia, cardiomegaly and renal hypoplasia syndrome (BCRHS) in Merino sheep maps to a 1.1-megabase region on ovine chromosome OAR2. Anim Genet. 2013;44:231–3.View ArticlePubMedGoogle Scholar
  33. 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.View ArticlePubMedPubMed CentralGoogle Scholar
  34. Randhawa IA, Khatkar MS, Thomson PC, Raadsma HW. Composite selection signals can localize the trait specific genomic regions in multi-breed populations of cattle and sheep. BMC Genet. 2014;15:34.View ArticlePubMedPubMed CentralGoogle Scholar
  35. Zhang L, Mousel MR, Wu X, Michal JJ, Zhou X, Ding B, et al. Genome-wide genetic diversity and differentially selected regions among Suffolk, Rambouillet, Columbia, Polypay, and Targhee sheep. PLoS One. 2013;8:e65942.View ArticlePubMedPubMed CentralGoogle Scholar
  36. Lv FH, Agha S, Kantanen J, Colli L, Stucki S, Kijas JW, et al. Adaptations to climate-mediated selective pressures in sheep. Mol Biol Evol. 2014;31:3324–43.View ArticlePubMedPubMed CentralGoogle Scholar
  37. 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.View ArticlePubMedPubMed CentralGoogle Scholar
  38. 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.View ArticlePubMedPubMed CentralGoogle Scholar
  39. Moioli B, Pilla F, Ciani E. Signatures of selection identify loci associated with fat tail in sheep. J Anim Sci. 2015;93:4660–9.View ArticlePubMedGoogle Scholar
  40. 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.View ArticlePubMedPubMed CentralGoogle Scholar
  41. Shimshony A. Pneumonia in East-Friesian lambs and crosses in Israel. Rev Sci Tech Off Int Epiz. 1983;2:467–71.View ArticleGoogle Scholar
  42. Duroudier NP, Tulah AS, Sayers I. Leukotriene pathway genetics and pharmacogenetics in allergy. Allergy. 2009;64:823–39.View ArticlePubMedGoogle Scholar
  43. Ro M, Kim S, Pyun JA, Shin C, Cho NH, Lee JY, et al. Association between arachidonate 5-lipoxygenase-activating protein (ALOX5AP) and lung function in a Korean population. Scand J Immunol. 2012;76:151–7.View ArticlePubMedGoogle Scholar
  44. Goot H, Gootwine E. Wool traits of Awassi and Assaf and their respective Booroola Merino F1 crosses. In Major genes for reproduction in sheep. Elsen JM, Bodin L, Thimonier J, editors. Proceedings of the International Workshop Les Colloques de l’INRA: 16–18 July 1990; Paris. 1991. p. 341–7.Google Scholar
  45. Garcia-Gamez E, Reverter A, Whan V, McWilliam SM, Arranz JJ, Kijas J. 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.View ArticlePubMedPubMed CentralGoogle Scholar
  46. Goot H, Bor A, Hasdai A, Zenou A, Gootwine E. Body and carcass composition of Awassi, Assaf, Booroola X Awassi and Booroola X Assaf ram lambs. In Major genes for reproduction in sheep. Elsen JM, Bodin L, Thimonier J, editors. Proceedings of the International Workshop Les Colloques de l’INRA: 16–18 July 1990; Paris. 1991. p. 335–40.Google Scholar
  47. 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.View ArticlePubMedPubMed CentralGoogle Scholar
  48. Wang H, Zhang L, Cao J, Wu M, Ma X, Liu Z, et al. Genome-wide specific selection in three domestic sheep breeds. PLoS One. 2015;10:e0128688.View ArticlePubMedPubMed CentralGoogle Scholar
  49. 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.View ArticlePubMedGoogle Scholar
  50. Aali M, Moradi-Shahrbabak M, Moradi-Shahrbabak H, Sadeghi M. Detecting novel SNPs and breed-specific haplotypes at calpastatin gene in Iranian fat- and thin-tailed sheep breeds and their effects on protein structure. Gene. 2014;537:132–9.View ArticlePubMedGoogle Scholar
  51. Bakhtiarizadeh MR, Moradi-Shahrbabak M, Ebrahimie E. Underlying functional genomics of fat deposition in adipose tissue. Gene. 2013;521:122–8.View ArticlePubMedGoogle Scholar
  52. 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.View ArticlePubMedGoogle Scholar
  53. 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.View ArticlePubMedGoogle Scholar
  54. Li MH, 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.View ArticleGoogle Scholar
  55. Han J, Kraft P, Nan H, Guo Q, Chen C, Qureshi A, et al. A genome-wide association study identifies novel alleles associated with hair color and skin pigmentation. PLoS Genet. 2008;4:e1000074.View ArticlePubMedPubMed CentralGoogle Scholar
  56. Sulem P, Gudbjartsson DF, Stacey SN, Helgason A, Rafnar T, Magnusson KP, et al. Genetic determinants of hair, eye and skin pigmentation in Europeans. Nat Genet. 2007;39:1443–52.View ArticlePubMedGoogle Scholar
  57. Sturm RA. Molecular genetics of human pigmentation diversity. Hum Mol Genet. 2009;18:R9–17.View ArticlePubMedGoogle Scholar
  58. Gutierrez-Gil B, Arranz JJ, Wiener P. An interpretive review of selective sweep studies in Bos taurus cattle populations: identification of unique and shared selection signals across breeds. Front Genet. 2015;6:167.PubMedPubMed CentralGoogle Scholar
  59. Color genes. 2014. Accessed 14 June 2016.
  60. Schleinitz D, Böttcher Y, Blüher M, Kovacs P. The genetics of fat distribution. Diabetologia. 2014;57:1276–86.View ArticlePubMedGoogle Scholar
  61. Wagner R, Machicao F, Fritsche A, Stefan N, Haring HU, Staiger H. The genetic influence on body fat distribution. Drug Discov Today Dis Mech. 2013;10:e5–13.View ArticleGoogle Scholar
  62. Shungin D, Winkler TW, Croteau-Chonka DC, Ferreira T, Locke AE, Mägi R, et al. New genetic loci link adipose and insulin biology to body fat distribution. Nature. 2015;518:187–96.View ArticlePubMedPubMed CentralGoogle Scholar
  63. Heid IM, Jackson AU, Randall JC, Winkler TW, Qi L, Steinthorsdottir V, et al. Meta-analysis identifies 13 new loci associated with waist-hip ratio and reveals sexual dimorphism in the genetic basis of fat distribution. Nat Genet. 2010;42:949–60.View ArticlePubMedPubMed CentralGoogle Scholar
  64. Lu Y, Day FR, Gustafsson S, Buchkovich ML, Na J, Bataille V, et al. New loci for body fat percentage reveal link between adiposity and cardiometabolic disease risk. Nat Commun. 2016;7:10495.View ArticlePubMedPubMed CentralGoogle Scholar
  65. Lam A, Rosov A, Shirak A, Seroussi E, Gootwine E. Genetic variation in Afec-Assaf ewes differing in their lamb survival rate at birth. J Anim Sci. 2013;91(E-Suppl 2):295.Google Scholar
  66. Legaz E, Álvarez I, Royo LJ, Fernández I, Gutiérrez JP, Goyache F. Genetic relationships between Spanish Assaf (Assaf.E) and Spanish native dairy sheep breeds. Small Rumin Res. 2008;80:39–44.View ArticleGoogle Scholar
  67. Fariello MI, Boitard S, Naya H, SanCristobal M, Servin B. Detecting signatures of selection through haplotype differentiation among hierarchically structured populations. Genetics. 2013;193:929–41.View ArticlePubMedPubMed CentralGoogle Scholar
  68. Ciani E, Crepaldi P, Nicoloso L, Lasagna E, Sarti FM, Moioli B, et al. Genome-wide analysis of Italian sheep diversity reveals a strong geographic pattern and cryptic relationships between breeds. Anim Genet. 2014;45:256–66.View ArticlePubMedGoogle Scholar
  69. Zhu C, Fan H, Yuan Z, Hu S, Zhang L, Wei C, et al. Detection of selection signatures on the X chromosome in three sheep breeds. Int J Mol Sci. 2015;16:20360–74.View ArticlePubMedPubMed CentralGoogle Scholar
  70. Moioli B, Scata MC, Steri R, Napolitano F, Catillo G. Signatures of selection identify loci associated with milk yield in sheep. BMC Genet. 2013;14:76.View ArticlePubMedPubMed CentralGoogle Scholar
  71. Garcia-Gámez E, Gutiérrez-Gill B, Sahana G, Sánchez JP, Bayón Y, Arranz JJ. GWA Analysis for milk production traits in dairy sheep and genetic support for a QTN influencing milk protein percentage in the LALBA gene. PLoS One. 2012;7:e47782.View ArticlePubMedPubMed CentralGoogle Scholar
  72. 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.View ArticlePubMedPubMed CentralGoogle Scholar
  73. Rowe S, McEwan JC, Hickey SM, Anderson RA, Hyndman D, Young E, et al. Genomic selection as a tool to decrease greenhouse gas emission from dual purpose New Zealand sheep. In Proceedings of the 10th World Congress of Genetics Applied to Livestock Production: 17–22 August 2014; Vancouver. 2014.Google Scholar
  74. Zhang L, Liu J, Zhao F, Ren H, Xu L, Lu J, et al. Genome-wide association studies for growth and meat production traits in sheep. PLoS One. 2013;8:e66569.View ArticlePubMedPubMed CentralGoogle Scholar
  75. McRae KM, McEwan JC, Dodds KG, Gemmell NJ. Signatures of selection in sheep bred for resistance or susceptibility to gastrointestinal nematodes. BMC Genomics. 2014;15:637.View ArticlePubMedPubMed CentralGoogle Scholar
  76. Riggio V, Pong-Wong R, Sallé G, Usai MG, Casu S, Moreno CR, et al. A joint analysis to identify loci underlying variation in nematode resistance in three European sheep populations. J Anim Breed Genet. 2014;131:426–36.View ArticlePubMedPubMed CentralGoogle Scholar
  77. Riggio V, Matika O, Pong-Wong R, Stear MJ, Bishop SC. Genome-wide association and regional heritability mapping to identify loci underlying variation in nematode resistance and body weight in Scottish Blackface lambs. Hered (Edinb). 2013;110:420–9.View ArticleGoogle Scholar
  78. 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.View ArticlePubMedPubMed CentralGoogle Scholar
  79. Kemper KE, Daetwyler HD, Visscher PM, Goddard ME. Comparing linkage and association analyses in sheep points to a better way of doing GWAS. Genet Res (Camb). 2012;94:191–203.View ArticleGoogle Scholar
  80. Elbeltagy AR, Kim ES, Aboul-Naga AM, Rischkowsky B, Rothschild MF. Heat stress in desert sheep and goats: signatures of natural selection and initiation of GWAS. In Proceedings of the International Plant and Animal Genome XXII Conference: 10–15 January 2014; San Diego. 2014.Google Scholar
  81. Gholizadeh M, Rahimi-Mianji G, Nejati-Javaremi A, Jan de Koning D, Jonas E. Genomewide association study to detect QTL for twinning rate in Baluchi sheep. J Genet. 2012;93:489–93.View ArticleGoogle Scholar
  82. Gholizadeh M, Rahimi-Mianji G, Nejati-Javaremi A. Genomewide association study of body weight traits in Baluchi sheep. J Genet. 2015;94:143–6.View ArticlePubMedGoogle Scholar
  83. Ren X, Yang GL, Peng WF, Zhao YX, Zhang M, Chen ZH, et al. A genome-wide association study identifies a genomic region for the polycerate phenotype in sheep (Ovis aries). Sci Rep. 2016;6:21111.View ArticlePubMedPubMed CentralGoogle Scholar
  84. Kijas JW, Ortiz JS, McCulloch R, James A, Brice B, Swain B, et al. Genetic diversity and investigation of polledness in divergent goat populations using 52 088 SNPs. Anim Genet. 2013;44:325–35.View ArticlePubMedGoogle Scholar
  85. Garcia-Gámez E, Reverter A, Whan V, McWilliam SM, Arranz JJ, International Sheep Genomics Consortium, 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.View ArticlePubMedPubMed CentralGoogle Scholar
  86. Demars J, Fabre S, Sarry J, Rossetti R, Gilbert H, Persani L, et al. Genome-wide association studies identify two novel BMP15 mutations responsible for an atypical hyperprolificacy phenotype in sheep. PLoS Genet. 2013;9:e1003482.View ArticlePubMedPubMed CentralGoogle Scholar
  87. Hazard D, Moreno C, Foulquié D, Delval E, Francois D, Bouix J, et al. Identification of QTLs for behavioral reactivity to social separation and humans in sheep using the OvineSNP50 BeadChip. BMC Genomics. 2014;15:778.View ArticlePubMedPubMed CentralGoogle Scholar
  88. White SN, Mousel MR, Herrmann-Hoesing LM, Reynolds JO, Leymaster KA, Neibergs HL, et al. Genome-wide association identifies multiple genomic regions associated with susceptibility to and control of ovine lentivirus. PLoS One. 2012;7:e47829.View ArticlePubMedPubMed CentralGoogle Scholar
  89. Mömke S, Kerkmann A, Wöhlke A, Ostmeier M, Hewicker-Trautwein M, Ganter M, et al. A frameshift mutation within LAMC2 is responsible for Herlitz type junctional epidermolysis bullosa (HJEB) in black headed mutton sheep. PLoS One. 2011;6:e18943.View ArticlePubMedPubMed CentralGoogle Scholar
  90. Pickering NK, Auvray B, Dodds KG, McEwan JC. Genomic prediction and genome-wide association study for dagginess and host internal parasite resistance in New Zealand sheep. BMC Genomics. 2015;16:958.View ArticlePubMedPubMed CentralGoogle Scholar
  91. Mucha S, Bunger L, Conington J. Genome-wide association study of footrot in Texel sheep. Genet Sel Evol. 2015;47:35.View ArticlePubMedPubMed CentralGoogle Scholar
  92. Moioli B, D’ Andrea S, De Grossi L, Sezzi E, De Sanctis B, Catillo G, et al. Genomic scan for identifying candidate genes for paratuberculosis resistance in sheep. Anim Prod Sci. 2016;56:1046–55.View ArticleGoogle Scholar
  93. Suarez-Vega A, Gutierrez-Gil B, Benavides J, Perez V, Tosser-Klopp G, Klopp C, et al. Combining GWAS and RNA-Seq approaches for detection of the causal mutation for hereditary junctional epidermolysis bullosa in sheep. PLoS One. 2015;10:e0126416.View ArticlePubMedPubMed CentralGoogle Scholar
  94. 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(Suppl 1):9–18.View ArticlePubMedGoogle Scholar
  95. Atlija M, Arranz JJ, Martinez-Valladares M, Gutierrez-Gil B. Detection and replication of QTL underlying resistance to gastrointestinal nematodes in adult sheep using the ovine 50 K SNP array. Genet Sel Evol. 2016;48:4.View ArticlePubMedPubMed CentralGoogle Scholar
  96. 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.View ArticlePubMedPubMed CentralGoogle Scholar
  97. Yang GL, Fu DL, Lang X, Wang YT, Cheng SR, Fang SL, et al. Mutations in MC1R gene determine black coat color phenotype in Chinese sheep. Sci World J. 2013;2013:675382.Google Scholar
  98. Shungin D, Winkler TW, Croteau-Chonka DC, Ferreira T, Locke AE, Magi R, et al. New genetic loci link adipose and insulin biology to body fat distribution. Nature. 2015;518:187–96.View ArticlePubMedPubMed CentralGoogle Scholar
  99. Schleinitz D, Bottcher Y, Bluher M, Kovacs P. The genetics of fat distribution. Diabetologia. 2014;57:1276–86.View ArticlePubMedGoogle Scholar
  100. Cristancho AG, Lazar MA. Forming functional fat: a growing understanding of adipocyte differentiation. Nat Rev Mol Cell Biol. 2011;12:722–34.View ArticlePubMedGoogle Scholar


© The Author(s) 2017