Genome-wide assessment of diversity and differentiation between original and modern Brown cattle populations.

Identifying genomic regions involved in the differences between breeds can provide information on genes that are under the influence of both artificial and natural selection. The aim of this study was to assess the genetic diversity and differentiation among four different Brown cattle populations (two original vs. two modern populations) and to characterize the distribution of runs of homozygosity (ROH) islands using the Illumina Bovine SNP50 BeadChip genotyping data. After quality control, 34 735 SNPs and 106 animals were retained for the analyses. Larger heterogeneity was highlighted for the original populations. Patterns of genetic differentiation, multidimensional scaling, and the neighboring joining tree distinguished the modern from the original populations. The FST -outlier identified several genes putatively involved in the genetic differentiation between the two groups, such as stature and growth, behavior, and adaptability to local environments. The ROH islands within both the original and the modern populations overlapped with QTL associated with relevant traits. In modern Brown (Brown Swiss and Italian Brown), ROH islands harbored candidate genes associated with milk production traits, in evident agreement with the artificial selection conducted to improve this trait in these populations. In original Brown (Original Braunvieh and Braunvieh), we identified candidate genes related with fat deposition, confirming that breeding strategies for the original Brown populations aimed to produce dual-purpose animals. Our study highlighted the presence of several genomic regions that vary between Brown populations, in line with their different breeding histories.

Comparative selection signature analyses identify genomic footprints in Reggiana cattle, the traditional breed of the Parmigiano-Reggiano cheese production system.

Reggiana is an autochthonous cattle breed reared mainly in the province of Reggio Emilia, located in the North of Italy. Reggiana cattle (originally a triple-purpose population largely diffused in the North of Italy) are characterised by a typical solid red coat colour. About 2500 cows of this breed are currently registered to its herd book. Reggiana is now considered a dual-purpose breed even if it is almost completely dedicated to the production of a mono-breed branded Protected Designation of Origin Parmigiano-Reggiano cheese, which is the main driver of the sustainable conservation of this local genetic resource. In this study, we provided the first overview of genomic footprints that characterise Reggiana and define the diversity of this local cattle breed. A total of 168 Reggiana sires (all bulls born over 35 years for which semen was available) and other 3321 sires from 3 cosmopolitan breeds (Brown, Holstein and Simmental) were genotyped with the Illumina BovineSNP50 panel. ADMIXTURE analysis suggested that Reggiana breed might have been influenced, at least in part, by the other three breeds included in this study. Selection signatures in the Reggiana genome were identified using three statistical approaches based on allele frequency differences among populations or on properties of haplotypes segregating in the populations (fixation index (FST); integrated haplotype score; cross-population extended haplotype homozygosity). We identified several regions under peculiar selection in the Reggiana breed, particularly on bovine chromosome (BTA) 6 in the KIT gene region, that is known to be involved in coat colour pattern distribution, and within the region of the LAP3, NCAPG and LCORL genes, that are associated with stature, conformation and carcass traits. Another already known region that includes the PLAG1 gene (BTA14), associated with conformation traits, showed a selection signature in the Reggiana cattle. On BTA18, a signal of selection included the MC1R gene that causes the red coat colour in cattle. Other selection sweeps were in regions, with high density of quantitative trait loci for milk production traits (on BTA20) and in several other large regions that might have contributed to shape and define the Reggiana genome (on BTA17 and BTA29). All these results, overall, indicate that the Reggiana genome might still contain several signs of its multipurpose and non-specialised utilisation, as already described for other local cattle populations, in addition to footprints derived by its ancestral origin and by its adaptation to the specialised Parmigiano-Reggiano cheese production system.

Genome-wide analyses reveal the regions involved in the phenotypic diversity in Sicilian pigs.

Nero Siciliano (Sicilian Black, SB) is a local pig breed generally of uniform black color. In addition to this officially recognized breed, there are animals showing morphological characteristics resembling the SB but with gray hair (Sicilian Grey, SG). The SG, compared with the SB, also shows a more compact structure with greater transverse diameters, higher average daily gains and lower thickness of the back fat. In this study, using the Illumina PorcineSNP60 BeadChip, we run genome-wide analyses to identify regions that may explain the phenotypic differences between SB (n = 21) and SG (n = 27) individuals. Combining the results of the two case-control approaches (GWAS and FST ), we identified two significant regions, one on SSC5 (95 401 083 bp) and one on SSC15 (55 051 435 bp), which contains several candidate genes related to growth traits in pig. The results of the Bayesian population differentiation approach identified a marker near the MGAT4C, a gene associated with average daily gain in pigs. Finally, scanning the genome for runs of homozygosity islands, we found that the two groups have different runs of homozygosity islands, with several candidate genes involved in coat color (in SG) or related to different pig performance traits (in SB). In summary, the two analyzed groups differed for several phenotypic traits, and genes involved in these traits (growth, meat traits and coat color) were detected. This study provided another contribution to the identification of genomic regions involved in phenotypic variability in local pig populations.

Genome-wide association studies for milk production traits in Valle del Belice sheep using repeated measures.

Genome-wide association studies (GWASes) have become a powerful tool for identifying genomic regions associated with important traits in livestock. Milk production traits in dairy sheep are measured at different time points during their life span. Using phenotypic data generated from longitudinal traits could improve the power of association studies but until now have received less attention in GWASes as a methodology and has not been implemented. The aim of this study was to carry out a GWAS for milk production traits in Valle del Belice sheep using repeated measures. After quality control, 469 ewes and 37 228 SNPs were retained for the analysis, and phenotypic data included 5586 test-day records for five milk production traits (milk yield, MY; fat yield and percentage, FY and F%; protein yield and percentage, PY and P%). Nine SNPs located within or close to known genes were found to be associated with milk production traits. In particular, rs398340969, associated with both milk yield and protein yield, is located within the DCPS gene. In addition, rs425417915 and rs417079368, both associated with both fat percentage and protein percentage, are located within the TTC7B gene and at 0.37 Mb within the SUCNR1 gene respectively. In summary, the use of repeated records was beneficial for mapping genomic regions affecting milk production traits in the Valle del Belice sheep.

Genome-wide association study reveals the locus responsible for microtia in Valle del Belice sheep breed.

Microtia is a congenital deformity of the outer ear with phenotypes varying from a small auricle to total absence (anotia). The genetic basis is still poorly understood, and very few studies have been performed in sheep. Valle del Belice sheep is a breed showing microtia. The aim of this study was to identify the potential genomic regions involved in microtia in sheep. A total of 40 individuals, 20 with microtia and 20 normal, were genotyped with the Illumina OvineSNP50 BeadChip. The comparison among the results from a genome-wide association study, Fisher's exact test and FST analysis revealed a single strong association signal: rs419889303 on chromosome 1, located within intron 3 of the CLRN1 gene. Our study suggests for the first time that this novel candidate gene is responsible for microtia in sheep. Additional analysis based on the sequencing would help confirm our findings and allow for the proposal of a precise genetic basis for microtia in sheep.

Genome-wide identification of runs of homozygosity islands and associated genes in local dairy cattle breeds.

Runs of homozygosity (ROH) are widely used as predictors of whole-genome inbreeding levels in cattle. They identify regions that have an unfavorable effect on a phenotype when homozygous, but also identify the genes associated with traits of economic interest present in these regions. Here, the distribution of ROH islands and enriched genes within these regions in four dairy cattle breeds were investigated. Cinisara (71), Modicana (72), Reggiana (168) and Italian Holstein (96) individuals were genotyped using the 50K v2 Illumina BeadChip. The genomic regions most commonly associated with ROHs were identified by selecting the top 1% of the single nucleotide polymorphisms (SNPs) most commonly observed in the ROH of each breed. In total, 11 genomic regions were identified in Cinisara and Italian Holstein, and eight in Modicana and Reggiana, indicating an increased ROH frequency level. Generally, ROH islands differed between breeds. The most homozygous region (>45% of individuals with ROH) was found in Modicana on chromosome 6 within a quantitative trail locus affecting milk fat and protein concentrations. We identified between 126 and 347 genes within ROH islands, which are involved in multiple signaling and signal transduction pathways in a wide variety of biological processes. The gene ontology enrichment provided information on possible molecular functions, biological processes and cellular components under selection related to milk production, reproduction, immune response and resistance/susceptibility to infection and diseases. Thus, scanning the genome for ROH could be an alternative strategy to detect genomic regions and genes related to important economic traits.

Runs of homozygosity reveal genome-wide autozygosity in Italian sheep breeds.

The availability of dense single nucleotide polymorphism (SNP) assays allows for the determination of autozygous segments based on runs of consecutive homozygous genotypes (ROH). The aim of the present study was to investigate the occurrence and distribution of ROH in 21 Italian sheep breeds using medium-density SNP genotypes in order to characterize autozygosity and identify genomic regions that frequently appeared in ROH within individuals, namely ROH islands. After filtering, the final number of animals and SNPs retained for analyses were 502 and 46 277 respectively. A total of 12 302 ROH were identified. The mean number of ROH per breed ranged from 10.58 (Comisana) to 44.54 (Valle del Belice). The average length of ROH across breeds was 4.55 Mb and ranged from 3.85 Mb (Biellese) to 5.51 Mb (Leccese). Valle del Belice showed the highest value of inbreeding on the basis of ROH (FROH  = 0.099), whereas Comisana showed the lowest (FROH  = 0.016), and high standard deviation values revealed high variability in autozygosity levels within each breed. Differences also existed in the length of ROH. Analysis of the distribution of ROH according to their size showed that, for all breeds, the majority of the detected ROH were <10 Mb in length, with a few long ROH >25 Mb. The levels of ROH that we estimated here reflect the inbreeding history of the investigated sheep breeds. These results also highlight that ancient and recent inbreeding have had an impact on the genome of the Italian sheep breeds and suggest that several animals have experienced recent autozygosity events. Comisana and Bergamasca appeared as the less consanguineous breeds, whereas Barbaresca, Leccese and Valle del Belice showed ROH patterns typically produced by recent inbreeding. Moreover, within the genomic regions most commonly associated with ROH, several candidate genes were detected.

Penalized classification for optimal statistical selection of markers from high-throughput genotyping: application in sheep breeds.

The identification of individuals' breed of origin has several practical applications in livestock and is useful in different biological contexts such as conservation genetics, breeding and authentication of animal products. In this paper, penalized multinomial regression was applied to identify the minimum number of single nucleotide polymorphisms (SNPs) from high-throughput genotyping data for individual assignment to dairy sheep breeds reared in Sicily. The combined use of penalized multinomial regression and stability selection reduced the number of SNPs required to 48. A final validation step on an independent population was carried out obtaining 100% correctly classified individuals. The results using independent analysis, such as admixture, F st, principal component analysis and random forest, confirmed the ability of these methods in selecting distinctive markers. The identified SNPs may constitute a starting point for the development of a SNP based identification test as a tool for breed assignment and traceability of animal products.

Preselection statistics and Random Forest classification identify population informative single nucleotide polymorphisms in cosmopolitan and autochthonous cattle breeds.

Commercial single nucleotide polymorphism (SNP) arrays have been recently developed for several species and can be used to identify informative markers to differentiate breeds or populations for several downstream applications. To identify the most discriminating genetic markers among thousands of genotyped SNPs, a few statistical approaches have been proposed. In this work, we compared several methods of SNPs preselection (Delta, F st and principal component analyses (PCA)) in addition to Random Forest classifications to analyse SNP data from six dairy cattle breeds, including cosmopolitan (Holstein, Brown and Simmental) and autochthonous Italian breeds raised in two different regions and subjected to limited or no breeding programmes (Cinisara, Modicana, raised only in Sicily and Reggiana, raised only in Emilia Romagna). From these classifications, two panels of 96 and 48 SNPs that contain the most discriminant SNPs were created for each preselection method. These panels were evaluated in terms of the ability to discriminate as a whole and breed-by-breed, as well as linkage disequilibrium within each panel. The obtained results showed that for the 48-SNP panel, the error rate increased mainly for autochthonous breeds, probably as a consequence of their admixed origin lower selection pressure and by ascertaining bias in the construction of the SNP chip. The 96-SNP panels were generally more able to discriminate all breeds. The panel derived by PCA-chrom (obtained by a preselection chromosome by chromosome) could identify informative SNPs that were particularly useful for the assignment of minor breeds that reached the lowest value of Out Of Bag error even in the Cinisara, whose value was quite high in all other panels. Moreover, this panel contained also the lowest number of SNPs in linkage disequilibrium. Several selected SNPs are located nearby genes affecting breed-specific phenotypic traits (coat colour and stature) or associated with production traits. In general, our results demonstrated the usefulness of Random Forest in combination to other reduction techniques to identify population informative SNPs.

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 breed's 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.

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.

Study of polymorphisms in the promoter region of ovine ß-lactoglobulin gene and phylogenetic analysis among the Valle del Belice breed and other sheep breeds considered as ancestors.

The aim of this work was to sequence the promoter region of ß-lactoglobulin (BLG) gene in four sheep breeds, in order to identify polymorphisms, infer and analyze haplotypes, and phylogenetic relationship among the Valle del Belice breed and the other three breeds considered as ancestors. Sequencing analysis and alignment of the obtained sequences showed the presence of 36 single nucleotide polymorphisms (SNPs) and one deletion. A total of 22 haplotypes found in "best" reconstruction were inferred considering the 37 polymorphic sites identified. Haplotypes were used for the reconstruction of a phylogenetic tree using the Neighbor-Joining algorithm. The number of polymorphisms identified showed high variability within breeds. Analysis of genetic diversity indexes showed that the Sarda breed presented the lowest nucleotide diversity, whereas the Comisana breed presented the highest one. Comparing the nucleotide diversity among breeds, the highest value was obtained between Valle del Belice and Pinzirita breeds, whereas the lowest one was between Valle del Belice and Sarda breeds. Considering that polymorphisms in the promoter region of BLG gene could have a functional role associated with milk composition, the lowest value of nucleotide diversity between Valle del Belice and Sarda breeds may be related to a higher similarity of milk composition of these two breeds compared to the others. Further analyses will be conducted in order to evaluate the possible correlation between the genetic diversity indexes and the BLG content in milk of our breeds.

Polymorphisms of beta-lactoglobulin promoter region in three Sicilian goat breeds.

Several beta-lactoglobulin (BLG) polymorphisms have been described within the proximal promoter region and coding region of the caprine gene, although no genetic variants affecting the protein amino acid composition and/or expression level have been characterized so far. Binding sites for several transcription factors (TFs) are present in the BLG promoter region. The aims of this work were to sequence the full-length promoter region of three Sicilian goat breeds in order to identify polymorphisms, analyze the identified haplotypes, search for differences between breeds for the presence of polymorphisms in this gene region, search for putative TFs binding sites, and check if polymorphisms lay within the identified TFs binding sites. The promoter region of BLG gene in Sicilian goat breeds showed high level of polymorphism due to the presence of 36 single nucleotide polymorphisms (SNPs). Association between polymorphic sites was computed within the whole sample analyzed and 18 haplotypes were inferred. Binding sites for three milk protein binding factors (MPBFs) and four nuclear factor-I (NF-I) were found within BLG promoter region based on the ovine sequence. The identification of some SNPs within TFs binding sites allowed hypothesizing the loss of TFs. Further studies are in progress to evaluate the effect of these mutations on binding affinity of TFs, the functional interaction of the TFs with the goat BLG promoter, and the relationship of the polymorphisms with BLG gene expression and milk production and composition.

Hierarchical structure of the Sicilian goats revealed by Bayesian analyses of microsatellite information.

Genetic structure and relationship amongst the main goat populations in Sicily (Girgentana, Derivata di Siria, Maltese and Messinese) were analysed using information from 19 microsatellite markers genotyped on 173 individuals. A posterior Bayesian approach implemented in the program STRUCTURE revealed a hierarchical structure with two clusters at the first level (Girgentana vs. Messinese, Derivata di Siria and Maltese), explaining 4.8% of variation (amovaФ(ST) estimate). Seven clusters nested within these first two clusters (further differentiations of Girgentana, Derivata di Siria and Maltese), explaining 8.5% of variation (amovaФ(SC) estimate). The analyses and methods applied in this study indicate their power to detect subtle population structure.

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.

A first comparative map of copy number variations in the sheep genome.

We carried out a cross species cattle-sheep array comparative genome hybridization experiment to identify copy number variations (CNVs) in the sheep genome analysing ewes of Italian dairy or dual-purpose breeds (Bagnolese, Comisana, Laticauda, Massese, Sarda, and Valle del Belice) using a tiling oligonucleotide array with ~385,000 probes designed on the bovine genome. We identified 135 CNV regions (CNVRs; 24 reported in more than one animal) covering ~10.5 Mb of the virtual sheep genome referred to the bovine genome (0.398%) with a mean and a median equal to 77.6 and 55.9 kb, respectively. A comparative analysis between the identified sheep CNVRs and those reported in cattle and goat genomes indicated that overlaps between sheep and both other species CNVRs are highly significant (P<0.0001), suggesting that several chromosome regions might contain recurrent interspecies CNVRs. Many sheep CNVRs include genes with important biological functions. Further studies are needed to evaluate their functional relevance.

Copy number variation and missense mutations of the agouti signaling protein (ASIP) gene in goat breeds with different coat colors.

In goats, classical genetic studies reported a large number of alleles at the Agouti locus with effects on coat color and pattern distribution. From these early studies, the dominant A(Wt) (white/tan) allele was suggested to cause the white color of the Saanen breed. Here, we sequenced the coding region of the goat ASIP gene in 6 goat breeds (Girgentana, Maltese, Derivata di Siria, Murciano-Granadina, Camosciata delle Alpi, and Saanen), with different coat colors and patterns. Five single nucleotide polymorphisms (SNPs) were identified, 3 of which caused missense mutations in conserved positions of the cysteine-rich carboxy-terminal domain of the protein (p.Ala96Gly, p.Cys126Gly, and p.Val128Gly). Allele and genotype frequencies suggested that these mutations are not associated or not completely associated with coat color in the investigated goat breeds. Moreover, genotyping and sequencing results, deviation from Hardy-Weinberg equilibrium, as well as allele copy number evaluation from semiquantitative fluorescent multiplex PCR, indicated the presence of copy number variation (CNV) in all investigated breeds. To confirm the presence of CNV and evaluate its extension, we applied a bovine-goat cross-species array comparative genome hybridization (aCGH) experiment using a custom tiling array based on bovine chromosome 13. aCGH results obtained for 8 goat DNA samples confirmed the presence of CNV affecting a region of less that 100 kb including the ASIP and AHCY genes. In Girgentana and Saanen breeds, this CNV might cause the A(Wt) allele, as already suggested for a similar structural mutation in sheep affecting the ASIP and AHCY genes, providing evidence for a recurrent interspecies CNV. However, other mechanisms may also be involved in determining coat color in these 2 breeds.

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

Short communication: casein haplotype variability in sicilian dairy goat breeds.

In the Mediterranean region, goat milk production is an important economic activity. In the present study, 4 casein genes were genotyped in 5 Sicilian goat breeds to 1) identify casein haplotypes present in the Argentata dell'Etna, Girgentana, Messinese, Derivata di Siria, and Maltese goat breeds; and 2) describe the structure of the Sicilian goat breeds based on casein haplotypes and allele frequencies. In a sample of 540 dairy goats, 67 different haplotypes with frequency >or=0.01 and 27 with frequency >or=0.03 were observed. The most common CSN1S1-CSN2-CSN1S2-CSN3 haplotype for Derivata di Siria and Maltese was FCFB (0.17 and 0.22, respectively), whereas for Argentata dell'Etna, Girgentana and Messinese was ACAB (0.06, 0.23, and 0.10, respectively). According to the haplotype reconstruction, Argentata dell'Etna, Girgentana, and Messinese breeds presented the most favorable haplotype for cheese production, because the casein concentration in milk of these breeds might be greater than that in Derivata di Siria and Maltese breeds. Based on a cluster analysis, the breeds formed 2 main groups: Derivata di Siria, and Maltese in one group, and Argentata dell'Etna and Messinese in the other; the Girgentana breed was between these groups but closer to the latter.