Investigation of phenotypic, genetic and genomic background of Milk spectra in Sarda dairy sheep.

Aim of this study was to analyse the genetic background of milk Fourier transform infrared (FTIR) spectra in dairy sheep. Individual milk FTIR spectra, with 1060 wavenumbers each, were available for 793 adult Sarda breed ewes genotyped at 45,813 SNP. The absorbance values of each wavenumber was analysed using a linear mixed model that included dim class, parity and lambing month as fixed effects and flock-test date and animal as random effects. The model was applied to estimate variance components and heritability and to perform a genome-wide association study for each wavenumber. Average h2 of wavenumbers absorbance was 0.13 ± 0.08, with the largest values observed in the regions associated with the characteristic bonds of carbonylic and methylenic groups of milk fat (h2 = 0.57 at 1724-1728 cm-1; and h2 = 0.34 at 2811-2834 cm-1, respectively). The absorbance values of wavenumbers were moderately correlated with the estimated heritabilities. After the Bonferroni correction, a total of nine markers were found to be significantly associated with 32 different wavenumbers. Of particular interest was the SNP s63269.1, mapped on chromosome 2, that was found to be associated with 27 wavenumbers. Genes previously found to be related to traits of interest (e.g. disease resistance, milk yield and quality, cheese firmness) are located close to the significant SNP. As expected, the heritability estimated for the absorbance of each wavenumbers seems to be associated with the related milk components.

Relationship between inbreeding and milk production traits in two Italian dairy sheep breeds.

The effects of inbreeding in livestock species breeds have been well documented and they have a negative impact on profitability. The objective of this study was to evaluate the levels of inbreeding in Sarda (SAR, n = 785) and Valle del Belice (VdB, n = 473) dairy sheep breeds and their impact on milk production traits. Two inbreeding coefficients (F) were estimated: using pedigree (FPED ), or runs of homozygosity (ROH; FROH ) at different minimum ROH lengths and different ROH classes. After the quality control, 38,779 single nucleotide polymorphisms remained for further analyses. A mixed-linear model was used to evaluate the impact of inbreeding coefficients on production traits within each breed. VdB showed higher inbreeding coefficients compared to SAR, with both breeds showing lower estimates as the minimum ROH length increased. Significant inbreeding depression was found only for milk yield, with a loss of around 7 g/day (for SAR) and 9 g/day (VdB) for a 1% increase of FROH . The present study confirms how the use of genomic information can be used to manage intra-breed diversity and to calculate the effects of inbreeding on phenotypic traits.

Analysis of ddRAD-seq data provides new insights into the genomic structure and patterns of diversity in Italian donkey populations.

With more than 150 recognized breeds, donkeys assume relevant economic importance, especially in developing countries. Even if the estimated number of heads worldwide is 53M, this species received less attention than other livestock species. Italy has traditionally been considered one of the cradles of European donkey breeding, and despite a considerable loss of biodiversity, today still counts nine autochthonous populations. A total of 220 animals belonging to nine different populations were genotyped using the double-digest restriction site associated DNA (ddRAD) sequencing to investigate the pattern of diversity using a multi-technique approach. A total of 418,602,730 reads were generated and successfully demultiplexed to obtain a medium-density SNP genotypes panel with about 27K markers. The diversity indices showed moderate levels of variability. The genetic distances and relationships, largely agree with the breeding history of the donkey populations under investigation. The results highlighted the separation of populations based on their genetic origin or geographical proximity between breeding areas, showed low to moderate levels of admixture, and indicated a clear genetic difference in some cases. For some breeds, the results also validate the success of proper management conservation plans. Identified runs of homozygosity islands, mapped within genomic regions related to immune response and local adaptation, are consistent with the characteristics of the species known for its rusticity and adaptability. This study is the first exhaustive genome-wide analysis of the diversity of Italian donkey populations. The results emphasized the high informativeness of genome-wide markers retrieved through the ddRAD approach. The findings take on great significance in designing and implementing conservation strategies. Standardized genotype arrays for donkey species would make it possible to combine worldwide datasets to provide further insights into the evolution of the genomic structure and origin of this important genetic resource.

Donkeys assume relevant economic importance in several countries worldwide. However, the genetic structure of these populations is less investigated compared to other species. The aim of this study was to investigate the genetic background of nine different Italian donkey populations. A total of 220 animals were genotyped with about 27K markers extracted by the double-digest restriction site associated DNA sequencing. The consistency of the results across different approaches agreed with the demographic history, the origin, and previous results on the nine donkey populations, suggesting that our conclusions are robust. Moreover, the results of the present study highlighted low to moderate levels of admixture and, for some breeds, confirmed the success of proper management conservation plans.