Skip to main content
Download PDF

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

Genetics Selection Evolution volume 49, Article number: 19 (2017) Cite this article

Abstract

Background

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.

Results

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.

Conclusions

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.

Background

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 [4–6]. 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
Full size table

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.

Methods

Ethics statement

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

Animals

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 (http://www.illumina.com). All the Awassi and some of the Afec-Assaf samples were genotyped as part of the International Sheep Genomic Consortium initiative (http://www.sheephapmap.org). 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 (https://www.ncbi.nlm.nih.gov) 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.

Bioinformatics

Candidate regions were examined using the sheep genome browsers Build v3.1 [26–28]. 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
Full size table
Fig. 1
figure 1

Manhattan plot of GWAS comparing Awassi and Afec-Assaf sheep

Full size image

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 [10–12, 18, 34–39].

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 [49–51] 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 [55–57], 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 [60–64] (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.

Conclusions

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].

References

  1. Gootwine E. Mini review: breeding Awassi and Assaf sheep for diverse management conditions. Trop Anim Health Prod. 2011;43:1289–96.

    Article  PubMed  Google 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.

  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.

    Article  Google 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.

    Article  CAS  PubMed  PubMed Central  Google 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.

    Article  CAS  PubMed  Google 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.

    Article  CAS  PubMed  Google 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.

    Article  PubMed  Google Scholar 

  8. Gootwine E, Pollott GE. Factors affecting milk production in Improved Awassi dairy ewes. Anim Sci. 2000;71:607–15.

    Article  Google 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.

    Article  CAS  PubMed  Google 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.

    Article  CAS  PubMed  PubMed Central  Google 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.

    Article  PubMed  PubMed Central  Google 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.

    Article  PubMed  Google 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.

    Article  PubMed  PubMed Central  Google 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Astle W, Balding DJ. Population structure and criptic relatedness in genetic association studies. Stat Sci. 2009;24:451–71.

    Article  Google 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.

    Article  CAS  PubMed  Google 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.

    Article  PubMed  PubMed Central  Google 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.

    Article  CAS  PubMed  Google 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Bonfield JK, Whitwham A. Gap5–editing the billion fragment sequence assembly. Bioinformatics. 2010;26:1699–703.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Li H, Durbin R. Fast and accurate long-read alignment with Burrows–Wheeler transform. Bioinformatics. 2010;26:589–95.

    Article  PubMed  PubMed Central  Google 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Randhawa IA, Khatkar MS, Thomson PC, Raadsma HW. A meta-assembly of selection signatures in cattle. PLoS One. 2016;11:e0153013.

    Article  PubMed  PubMed Central  Google 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.

    Article  CAS  PubMed  PubMed Central  Google 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. http://www.livestockgenomics.csiro.au/sheep/oar3.1.php.

  27. http://genome.ucsc.edu. Accessed 14 June 2016.

  28. http://asia.ensembl.org. Accessed 14 June 2016.

  29. http://www.ncbi.nlm.nih.gov/. Accessed 14 June 2016.

  30. http://www.genecards.org/. 14 June 2016.

  31. http://sphinx.vet.unimelb.edu.au/QTLdb/. 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.

    Article  CAS  PubMed  Google 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.

    Article  PubMed  PubMed Central  Google 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.

    Article  PubMed  PubMed Central  Google 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.

    Article  CAS  PubMed  PubMed Central  Google 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.

    Article  CAS  PubMed  PubMed Central  Google 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.

    Article  PubMed  PubMed Central  Google 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.

    Article  PubMed  PubMed Central  Google 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.

    Article  CAS  PubMed  Google 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.

    Article  PubMed  PubMed Central  Google Scholar 

  41. Shimshony A. Pneumonia in East-Friesian lambs and crosses in Israel. Rev Sci Tech Off Int Epiz. 1983;2:467–71.

    Article  Google Scholar 

  42. Duroudier NP, Tulah AS, Sayers I. Leukotriene pathway genetics and pharmacogenetics in allergy. Allergy. 2009;64:823–39.

    Article  CAS  PubMed  Google 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.

    Article  CAS  PubMed  Google 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.

  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.

    Article  CAS  PubMed  PubMed Central  Google 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.

  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.

    Article  CAS  PubMed  PubMed Central  Google 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.

    Article  PubMed  PubMed Central  Google 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.

    Article  CAS  PubMed  Google 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.

    Article  CAS  PubMed  Google Scholar 

  51. Bakhtiarizadeh MR, Moradi-Shahrbabak M, Ebrahimie E. Underlying functional genomics of fat deposition in adipose tissue. Gene. 2013;521:122–8.

    Article  CAS  PubMed  Google 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.

    Article  CAS  PubMed  Google 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.

    Article  CAS  PubMed  Google 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.

    Article  CAS  Google 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.

    Article  PubMed  PubMed Central  Google 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.

    Article  CAS  PubMed  Google Scholar 

  57. Sturm RA. Molecular genetics of human pigmentation diversity. Hum Mol Genet. 2009;18:R9–17.

    Article  CAS  PubMed  Google 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.

    PubMed  PubMed Central  Google Scholar 

  59. Color genes. 2014. http://www.espcr.org/micemut. 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.

    Article  CAS  PubMed  Google 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.

    Article  Google 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.

    Article  CAS  PubMed  PubMed Central  Google 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.

    Article  CAS  PubMed  PubMed Central  Google 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.

    Article  CAS  PubMed  PubMed Central  Google 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.

    Article  Google 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.

    Article  PubMed  PubMed Central  Google 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.

    Article  CAS  PubMed  Google 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.

    Article  CAS  PubMed  PubMed Central  Google 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.

    Article  PubMed  PubMed Central  Google 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.

    Article  PubMed  PubMed Central  Google 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.

    Article  PubMed  PubMed Central  Google 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.

  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.

    Article  CAS  PubMed  PubMed Central  Google 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.

    Article  PubMed  PubMed Central  Google 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.

    Article  CAS  PubMed  PubMed Central  Google 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.

    Article  CAS  Google 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.

    Article  PubMed  PubMed Central  Google 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.

    Article  CAS  Google 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.

  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.

    Article  Google 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.

    Article  PubMed  Google 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.

    Article  CAS  PubMed  PubMed Central  Google 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.

    Article  CAS  PubMed  Google 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.

    Article  PubMed  PubMed Central  Google 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.

    Article  CAS  PubMed  PubMed Central  Google 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.

    Article  PubMed  PubMed Central  Google 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.

    Article  CAS  PubMed  PubMed Central  Google 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.

    Article  PubMed  PubMed Central  Google 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.

    Article  PubMed  PubMed Central  Google Scholar 

  91. Mucha S, Bunger L, Conington J. Genome-wide association study of footrot in Texel sheep. Genet Sel Evol. 2015;47:35.

    Article  PubMed  PubMed Central  Google 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.

    Article  CAS  Google 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.

    Article  PubMed  PubMed Central  Google 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.

    Article  PubMed  Google 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.

    Article  PubMed  PubMed Central  Google 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.

    Article  CAS  PubMed  PubMed Central  Google 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  99. Schleinitz D, Bottcher Y, Bluher M, Kovacs P. The genetics of fat distribution. Diabetologia. 2014;57:1276–86.

    Article  CAS  PubMed  Google Scholar 

  100. Cristancho AG, Lazar MA. Forming functional fat: a growing understanding of adipocyte differentiation. Nat Rev Mol Cell Biol. 2011;12:722–34.

    Article  CAS  PubMed  Google Scholar 

Download references

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.

Acknowledgements

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.

Author information

Authors and Affiliations

  1. Institute of Animal Science, ARO, The Volcani Center, PO Box 15159, 7528809, Rishon LeZion, Israel

    Eyal Seroussi, Alexander Rosov, Andrey Shirak, Alon Lam & Elisha Gootwine

Authors
  1. Eyal Seroussi
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Alexander Rosov
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Andrey Shirak
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Alon Lam
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Elisha Gootwine
    View author publications

    You can also search for this author in PubMed Google Scholar

Corresponding author

Correspondence to Elisha Gootwine.

Additional files

12711_2017_296_MOESM1_ESM.docx

Additional file 1: Table S1. Genome-wide association studies in sheep reviewed in the current study [10–14, 17, 18, 32–39, 47, 48, 54, 65, 67–95].

Additional file 2: Table S2. Genomic relationship within and between local Awassi flocks.

12711_2017_296_MOESM3_ESM.docx

Additional file 3: Table S3. Variable genomic regions (VGR) that differ between the Awassi and the Afec-Assaf according to PLINK and EMMAX analyses. SNPs were ranked according to their P-values, and the selected regions included the top 0.1% markers and their neighboring (1.5 Mb upstream and 1.5 Mb downstream) markers ranked in the top 1% of the markers. Two or more genomic regions were combined if they overlapped. Region position was determined using sheep genome assembly version 3.1. Regions identified by both analyses are marked in gray.

12711_2017_296_MOESM4_ESM.docx

Additional file 4: Figure S1. Different horn phenotypes in Awassi (A to D) and Afec-Assaf (E to H) sheep. (A) Polled, female. (B) Knobs, female. (C) Horns, female. (D) Horns, male. (E) Polled, female. (F) Knobs, male. (G) Knobs, male. (H) Scurs, male.

12711_2017_296_MOESM5_ESM.docx

Additional file 5: Table S4. Association between presence of a 1.8-kb insertion at the 3′-UTR of RXFP2 and horn phenotype in Awassi and Afec-Assaf sheep. A 428-bp PCR product indicates the presence of the insertion [19].

12711_2017_296_MOESM6_ESM.docx

Additional file 6: Table S5. Candidate genes for sheep coat colors and patterns [10, 11, 16, 34, 35, 38, 45, 52–54, 96, 97].

12711_2017_296_MOESM7_ESM.docx

Additional file 7. Candidate gene sequences. This file contains sequences of candidate genes in a FASTA-like format. Each sequence header consists of the gene symbol followed by its reference GenBank accession number and the breed name. Polymorphic nucleotide positions are indicated in blue and red font colors.

12711_2017_296_MOESM8_ESM.docx

Additional file 8: Figure S2. Variability in fat-tail phenotypes in Assaf sheep. Awassi-like fat tail (A). Fat tails with different amounts of fat deposition and different lengths (B to H).

12711_2017_296_MOESM9_ESM.docx

Additional file 9: Table S6. Candidate genes and SNP markers associated with carrying fat tail in sheep [38, 39, 47, 49–51].

12711_2017_296_MOESM10_ESM.docx

Additional file 10: Table S7. Genes known to affect coat pigmentation in humans, cattle and other mammals [55–58].

12711_2017_296_MOESM11_ESM.docx

Additional file 11: Table S8. Human genes associated with fat distribution [61, 63, 64, 98–100].

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), 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 (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Seroussi, E., Rosov, A., Shirak, A. et al. Unveiling genomic regions that underlie differences between Afec-Assaf sheep and its parental Awassi breed. Genet Sel Evol 49, 19 (2017). https://doi.org/10.1186/s12711-017-0296-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12711-017-0296-3

Keywords

Genetics Selection Evolution

ISSN: 1297-9686