B. Portolano

Coat colours in the Massese sheep breed are associated with mutations in the agouti signalling protein (ASIP) and melanocortin 1 receptor (MC1R) genes

Massese is an Italian dairy sheep breed characterized by animals with black skin and horns and black or apparent grey hairs. Owing to the presence of these two coat colour types, this breed can be considered an interesting model to evaluate the effects of coat colour gene polymorphisms on this phenotypic trait. Two main loci have been already shown to affect coat colour in sheep: Agouti and Extension coding for the agouti signalling protein (ASIP) and melanocortin 1 receptor (MC1R) genes, respectively. The Agouti locus is affected by a large duplication including the ASIP gene that may determine the Agouti white and tan allele (A(Wt)). Other disrupting or partially inactivating mutations have been identified in exon 2 (a deletion of 5 bp, D(5); and a deletion of 9 bp, D(9)) and in exon 4 (g.5172T>A, p.C126S) of the ASIP gene. Three missense mutations in the sheep MC1R gene cause the dominant black E(D) allele (p.M73K and p.D121N) and the putative recessive e allele (p.R67C). Here, we analysed these ASIP and MC1R mutations in 161 Massese sheep collected from four flocks. The presence of one duplicated copy allele including the ASIP gene was associated with grey coat colour (P = 9.4E-30). Almost all animals with a duplicated copy allele (37 out of 41) showed uniform apparent grey hair and almost all animals without a duplicated allele (117 out of 120) were completely black. Different forms of duplicated alleles were identified in Massese sheep including, in almost all cases, copies with exon 2 disrupting or partially inactivating mutations making these alleles different from the A(Wt) allele. A few exceptions were observed in the association between ASIP polymorphisms and coat colour: three grey sheep did not carry any duplicated copy allele and four black animals carried a duplicated copy allele. Of the latter four sheep, two carried the E(D) allele of the MC1R gene that may be the cause of their black coat colour. The coat colour of all other black animals may be determined by non-functional ASIP alleles (non-agouti alleles, A(a)) and in a few cases by the E(D) Extension allele. At least three frequent ASIP haplotypes ([D(5):g.5172T], [N:g.5172A] and [D(5):g.5172A]) were detected (organized into six different diplotypes). In conclusion, the results indicated that coat colours in the Massese sheep breed are mainly derived by combining ASIP and MC1R mutations.

Genomic inbreeding estimation in small populations: evaluation of runs of homozygosity in three local dairy cattle breeds

In the local breeds with small population size, one of the most important problems is the increase of inbreeding coefficient (F). High levels of inbreeding lead to reduced genetic diversity and inbreeding depression. The availability of high-density single nucleotide polymorphism (SNP) arrays has facilitated the quantification of F by genomic markers in farm animals. Runs of homozygosity (ROH) are contiguous lengths of homozygous genotypes and represent an estimate of the degree of autozygosity at genome-wide level. The current study aims to quantify the genomic F derived from ROH (F ROH) in three local dairy cattle breeds. F ROH values were compared with F estimated from the genomic relationship matrix (F GRM), based on the difference between observed v. expected number of homozygous genotypes (F HOM) and the genomic homozygosity of individual i (F MOL i ). The molecular coancestry coefficient (f MOL ij ) between individuals i and j was also estimated. Individuals of Cinisara (71), Modicana (72) and Reggiana (168) were genotyped with the 50K v2 Illumina BeadChip. Genotypes from 96 animals of Italian Holstein cattle breed were also included in the analysis. We used a definition of ROH as tracts of homozygous genotypes that were >4 Mb. Among breeds, 3661 ROH were identified. Modicana showed the highest mean number of ROH per individual and the highest value of F ROH, whereas Reggiana showed the lowest ones. Differences among breeds existed for the ROH lengths. The individuals of Italian Holstein showed high number of short ROH segments, related to ancient consanguinity. Similar results showed the Reggiana with some extreme animals with segments covering 400 Mb and more of genome. Modicana and Cinisara showed similar results between them with the total length of ROH characterized by the presence of large segments. High correlation was found between F HOM and F ROH ranged from 0.83 in Reggiana to 0.95 in Cinisara and Modicana. The correlations among F ROH and other estimated F coefficients were generally lower ranged from 0.45 (F MOL i -F ROH) in Cinisara to 0.17 (F GRM-F ROH) in Modicana. On the basis of our results, recent inbreeding was observed in local breeds, considering that 16 Mb segments are expected to present inbreeding up to three generations ago. Our results showed the necessity of implementing conservation programs to control the rise of inbreeding and coancestry in the three Italian local dairy cattle breeds.

The genome-wide structure of two economically important indigenous Sicilian cattle breeds.

Genomic technologies, such as high-throughput genotyping based on SNP arrays, provided background information concerning genome structure in domestic animals. The aim of this work was to investigate the genetic structure, the genome-wide estimates of inbreeding, coancestry, effective population size (Ne), and the patterns of linkage disequilibrium (LD) in 2 economically important Sicilian local cattle breeds, Cinisara (CIN) and Modicana (MOD), using the Illumina Bovine SNP50K v2 BeadChip. To understand the genetic relationship and to place both Sicilian breeds in a global context, genotypes from 134 other domesticated bovid breeds were used. Principal component analysis showed that the Sicilian cattle breeds were closer to individuals of Bos taurus taurus from Eurasia and formed nonoverlapping clusters with other breeds. Between the Sicilian cattle breeds, MOD was the most differentiated, whereas the animals belonging to the CIN breed showed a lower value of assignment, the presence of substructure, and genetic links with the MOD breed. The average molecular inbreeding and coancestry coefficients were moderately high, and the current estimates of Ne were low in both breeds. These values indicated a low genetic variability. Considering levels of LD between adjacent markers, the average r(2) in the MOD breed was comparable to those reported for others cattle breeds, whereas CIN showed a lower value. Therefore, these results support the need of more dense SNP arrays for a high-power association mapping and genomic selection efficiency, particularly for the CIN cattle breed. Controlling molecular inbreeding and coancestry would restrict inbreeding depression, the probability of losing beneficial rare alleles, and therefore the risk of extinction. The results generated from this study have important implications for the development of conservation and/or selection breeding programs in these 2 local cattle breeds.

Genome-wide analysis in endangered populations: A case study in Barbaresca sheep

Analysis of genomic data is becoming increasingly common in the livestock industry and the findings have been an invaluable resource for effective management of breeding programs in small and endangered populations. In this paper, with the goal of highlighting the potential of genomic analysis for small and endangered populations, genome-wide levels of linkage disequilibrium, measured as the squared correlation coefficient of allele frequencies at a pair of loci, effective population size, runs of homozygosity (ROH) and genetic diversity parameters, were estimated in Barbaresca sheep using Illumina OvineSNP50K array data. Moreover, the breeds genetic structure and its relationship with other breeds were investigated. Levels of pairwise linkage disequilibrium decreased with increasing distance between single nucleotide polymorphisms. An average correlation coefficient <0.25 was found for markers located up to 50 kb apart. Therefore, these results support the need to use denser single nucleotide polymorphism panels for high power association mapping and genomic selection efficiency in future breeding programs. The estimate of past effective population size ranged from 747 animals 250 generations ago to 28 animals five generations ago, whereas the contemporary effective population size was 25 animals. A total of 637 ROH were identified, most of which were short (67%) and ranged from 1 to 10 Mb. The genetic analyses revealed that the Barbaresca breed tended to display lower variability than other Sicilian breeds. Recent inbreeding was evident, according to the ROH analysis. All the investigated parameters showed a comparatively narrow genetic base and indicated an endangered status for Barbaresca. Multidimensional scaling, model-based clustering, measurement of population differentiation, neighbor networks and haplotype sharing distinguished Barbaresca from other breeds, showed a low level of admixture with the other breeds considered in this study, and indicated clear genetic differences compared with other breeds. Attention should be given to the conservation of Barbaresca due to its critical conservation status. In this context, genomic information may have a crucial role in management of small and endangered populations.

Comparison between different selection criteria in the genetic evaluation of Valle del Belice sheep

Lactation length in dairy sheep affects milk yield like other genetic and environmental factors. The length of the production period is affected by management decisions such as culling, mating and particularly ranking of animals with different parity and lambing in different months or seasons. Moreover the low heritability of lactation length (Barillet and Boichard, 1987; Dahlin et al., 1998) does not allow its use as a selection criterion. For this reason to achieve a good reliability in phenotypic and genetic evaluation of dairy species, production variability caused by systematic environmental effects must be removed. This is of particular interest for dairy sheep and goats reared in Sicily, where the typical production system is based on pasture, and related food availability is strongly affected by seasonal and annual climatic variations, which results in considerable variations in daily yields........

Effect of somatic cell count level on functional longevity in Valle del Belice dairy sheep assessed using survival analysis.

The objectives of this study were to evaluate the effect of somatic cell count (SCC) on functional longevity and to estimate the heritability of functional longevity using survival analysis in Valle del Belice dairy sheep. A total of 4,880 lactations of 2,190 ewes from 11 flocks were used. In this study, SCC was considered as an indication of subclinical mastitis. In case of clinical cases, identified by the technicians at milking time, test-day weights and milk samples of those ewes were not considered. Somatic cells were analyzed as counts, without any transformation, and were grouped in 3 classes based on the observed SCC maximum (mxSCC). The mxSCC classes, expressed as 10(3) cells/mL, were classified as 1 if mxSCC <or= 500, 2 if 500 < mxSCC < 1,000, and 3 if mxSCC >or= 1,000. An increase in SCC was associated with an increased hazard of being culled. Ewes in the highest class of SCC on a test-day had a 20% higher hazard of being culled than those in the lowest class. Therefore, SCC played a role in culling decisions of Valle del Belice dairy sheep farmers. The heritability estimate for functional longevity was 7% on the logarithmic scale and 11% on the real scale, indicating that selection for this trait is possible in sheep. The flock-year-season effect explained 19% of the variation on the logarithmic scale and 27% of the variation on the real scale.

Effect of weather conditions on somatic cell score in Sicilian Valle del Belice ewes

Abstract Mastitis susceptibility of Valle del Belice ewes from the south of Sicily to temperature, humidity, precipitation, solar radiation, sun hours, air pressure, wind-speed and wind-direction measured by a public weather station was investigated using 65,848 test-day somatic cell score (SCS) records of 5,237 ewes. All weather parameters had an effect on SCS in a regression approach. Extreme values of maximum and minimum temperature-humidity indices resulted in increased SCS. Higher precipitation, solar radiation and sun hours resulted in increased SCS, whereas higher air pressure and wind speed resulted in reduced SCS.

Effects of different simplified milk recording methods on genetic evaluation with Test-Day animal model

Abstract The aims of the present study were to compare estimated breeding values (EBV) for milk yield using different testing schemes with a test-day animal model and to evaluate the effect of different testing schemes on the ranking of top sheep. Alternative recording schemes that use less information than that currently obtained with a monthly test-day schedule were employed to estimate breeding values. A random regression animal mixed model that used a spline function of days in milk was fitted. EBVs obtained with alternative recording schemes showed different degrees of Spearman correlation with EBVs obtained using the monthly recording scheme. These correlations ranged from 0.77 to 0.92. A reduction in accuracy and intensity of selection could be anticipated if these alternative schemes are used; more research in this area is needed to reduce the costs of test-day recording.

Variance component estimation for SCS in different local breeds and species

The Italian Bracco is one of the oldest pointing dog breed, used for hunting ever since the Renaissance time; paintings of the 14th century show hunting sceneries with dogs similar to the present day Bracco. The breed has been officially registered by ENCI (the Italian cynological club) in 1949, when the definitive standard was established. In this work, we report the first results of a study aimed at measuring morphological traits in this breed, as a part of a more comprehensive study whose objective is to identify the characters that have the potential of being genetically improved. R. Ciampolini1, F. Cecchi1, A. Bramante1, F. Casetti2, S. Presciuttini1

Genome-wide scan of fat-tail sheep identifies signals of selection for fat deposition and adaptation

Fat tail in sheep represents a valuable energy reserve for facing future climate changes. The identification of genes with a role in the fat-tail phenotype may contribute to understanding the physiology of fat deposition and the mechanisms of adaptation. Genotypic data obtained with the OvineSNP50K array in 13 thin-tail sheep breeds from Italy were used to identify selection signatures of fat tail through pairwise thin- versus fat-tail sheep breed comparisons, with the following fat-tail breeds of the Mediterranean area: two unique Italian fat-tail breeds (Barbaresca and Laticauda), a Barbary sheep breed from Libya, Ossimi breed from Egypt, Cyprus Fat-Tail and Chios from the Greek islands Cyprus and Chios, respectively. Fst and χ2 values obtained for >40 000 polymorphic markers allowed confirmation of 12 fat-tail associations that were previously reported in Chinese and Iranian breeds. Two of these signals – on OAR 7 and OAR 13 – are in the proximity of two genes – VRTN and BMP2 – with a role in the variation of vertebral number and in fat-tail formation respectively. Two identified signals on OAR 6 and OAR 15 encompass two genes, PDGFRA and PDGFD, involved in the differentiation of preadipocytes. Further signals detected herein were reported in Chinese sheep as signatures of adaptation to desert areas. For several of the detected associations, the known role in either fat deposition or adaptation, thus contributing to revealing the molecular basis underlying mechanisms of energy storage and climate adaptation.