Increased inbreeding is an inevitable consequence of selection in livestock populations. The analysis of high-density single nucleotide polymorphisms (SNPs) facilitates the identification of long and uninterrupted runs of homozygosity (ROH) that can be used to identify chromosomal regions that are identical by descent. In this work, the distribution of ROH of different lengths in five Italian cattle breeds is described. A total of 4095 bulls from five cattle breeds (2093 Italian Holstein, 749 Italian Brown, 364 Piedmontese, 410 Marchigiana and 479 Italian Simmental) were genotyped at 54K SNP loci. ROH were identified and used to estimate molecular inbreeding coefficients (FROH ), which were compared with inbreeding coefficients estimated from pedigree information (FPED ) and using the genomic relationship matrix (FGRM ). The average number of ROH per animal ranged from 54 ± 7.2 in Piedmontese to 94.6 ± 11.6 in Italian Brown. The highest number of short ROH (related to ancient consanguinity) was found in Piedmontese, followed by Simmental. The Italian Brown and Holstein had a higher proportion of longer ROH distributed across the whole genome, revealing recent inbreeding. The FPED were moderately correlated with FROH > 1 Mb (0.662, 0.700 and 0.669 in Italian Brown, Italian Holstein and Italian Simmental respectively) but poorly correlated with FGRM (0.134, 0.128 and 0.448 for Italian Brown, Italian Holstein and Italian Simmental respectively). The inclusion of ROH > 8 Mb in the inbreeding calculation improved the correlation of FROH with FPED and FGRM . ROH are a direct measure of autozygosity at the DNA level and can overcome approximations and errors resulting from incomplete pedigree data. In populations with high linkage disequilibrium (LD) and recent inbreeding (e.g. Italian Holstein and Italian Brown), a medium-density marker panel, such as the one used here, may provide a good estimate of inbreeding. However, in populations with low LD and ancient inbreeding, marker density would have to be increased to identify short ROH that are identical by descent more precisely.

Milk fatty acid composition is a parameter of great interest for evaluation of nutritional quality of milk. Stearoyl-CoA desaturase (SCD) is a key enzyme in mammary lipid metabolism because it is able to add a double bond in the cis delta9-position in a large spectrum of medium- and long-chain fatty acids. A polymorphism with 2 alleles (A and V) in the fifth exon of the SCD gene has been reported. The effect of SCD genotype on individual milk fatty acid composition and on cis-9 unsaturated/saturated fatty acid ratios of 297 Holstein Italian Friesian cows was investigated in this paper. The SCD genotypes were determined by using a single strand conformation polymorphism method. Relative frequencies of SCD genotypes were 27, 60, and 13% for AA, AV, and VV, respectively. Milk of AA cows had a greater content of cis-9 C18:1 and total monounsaturated fatty acids and a higher C14:1/C14 ratio than did milk of VV cows. The relative contribution of SCD genotype to variation of monounsaturated fatty acids, cis-9 C18:1, and cis-9 C14:1 was 5, 4, and 7.7%, respectively. No significant differences were detected between SCD genotypes in the milk content of cis-9, trans-11 C18:2. Results of the present work provide some indication of an association between SCD locus and the fatty acid profile in the examined sample of Italian Holsteins, thus suggesting a possible role of this gene in the genetic variation of milk nutritional properties.

Italy counts several sheep breeds, arisen over centuries as a consequence of ancient and recent genetic and demographic events. To finely reconstruct genetic structure and relationships between Italian sheep, 496 subjects from 19 breeds were typed at 50K single nucleotide polymorphism loci. A subset of foreign breeds from the Sheep HapMap dataset was also included in the analyses. Genetic distances (as visualized either in a network or in a multidimensional scaling analysis of identical by state distances) closely reflected geographic proximity between breeds, with a clear north-south gradient, likely because of high levels of past gene flow and admixture all along the peninsula. Sardinian breeds diverged more from other breeds, a probable consequence of the combined effect of ancient sporadic introgression of feral mouflon and long-lasting genetic isolation from continental sheep populations. The study allowed the detection of previously undocumented episodes of recent introgression (Delle Langhe into the endangered Altamurana breed) as well as signatures of known, or claimed, historical introgression (Merino into Sopravissana and Gentile di Puglia; Bergamasca into Fabrianese, Appenninica and, to a lesser extent, Leccese). Arguments that would question, from a genomic point of view, the current breed classification of Bergamasca and Biellese into two separate breeds are presented. Finally, a role for traditional transhumance practices in shaping the genetic makeup of Alpine sheep breeds is proposed. The study represents the first exhaustive analysis of Italian sheep diversity in an European context, and it bridges the gap in the previous HapMap panel between Western Mediterranean and Swiss breeds.

We investigated the potential of using multivariate factor analysis to extract metabolic information from data on the quantity and quality of milk produced under different management systems. We collected data from individual milk samples taken from 1,158 Brown Swiss cows farmed in 85 traditional or modern herds in Trento Province (Italy). Factor analysis was carried out on 47 individual fatty acids, milk yield, and 5 compositional milk traits (fat, protein, casein, and lactose contents, somatic cell score). According to a previous study on multivariate factor analysis, a variable was considered to be associated with a specific factor if the absolute value of its correlation with the factor was ≥0.60. The extracted factors were representative of the following 12 groups of fatty acids or functions: de novo fatty acids, branched fatty acid-milk yield, biohydrogenation, long-chain fatty acids, desaturation, short-chain fatty acids, milk protein and fat contents, odd fatty acids, conjugated linoleic acids, linoleic acid, udder health, and vaccelenic acid. Only 5 fatty acids showed small correlations with these groups. Factor analysis suggested the existence of differences in the metabolic pathways for de novo short- and medium-chain fatty acids and Δ-desaturase products. An ANOVA of factor scores highlighted significant effects of the dairy farming system (traditional or modern), season, herd/date, parity, and days in milk. Factor behavior across levels of fixed factors was consistent with current knowledge. For example, compared with cows farmed in modern herds, those in traditional herds had higher scores for branched fatty acids, which were inversely associated with milk yield; primiparous cows had lower scores than older cows for de novo fatty acids, probably due to a larger contribution of lipids mobilized from body depots on milk fat yield. The statistical approach allowed us to reduce a large number of variables to a few latent factors with biological meaning and able to represent groups of fatty acids with a common origin and function. Multivariate factor analysis would therefore be a valuable tool for studying the influence of different production environments and individual animal factors on milk fatty acid composition, and for developing nutritional strategies able to manipulate the milk fatty acid profile according to consumer demand.

This review gives an overview of the mathematical modelling of lactation curves in dairy cattle. Over the last ninety years, the development of this field of study has followed the main requirements of the dairy cattle industry. Non-linear parametric functions have represented the preferred tools for modelling average curves of homogeneous groups of animals, with the main aim of predicting yields for management purposes. The increased availability of records per individual lactations and the genetic evaluation based on test day records has shifted the interest of modellers towards more flexible and general linear functions, as polynomials or splines. Thus the main interest of modelling is no longer the reconstruction of the general pattern of the phenomenon but the fitting of individual deviations from an average curve. Other specific approaches based on the modelling of the correlation structure of test day records within lactation, such as mixed linear models or principal component analysis, have been used to test the statistical significance of fixed effects in dairy experiments or to create new variables expressing main lactation curve traits. The adequacy of a model is not an absolute requisite, because it has to be assessed according to the specific purpose it is used for. Occurrence of extended lactations and of new productive and functional traits to be described and the increase of records coming from automatic milking systems likely will represent some of the future challenges for the mathematical modelling of the lactation curve in dairy cattle

The aim of this study was to elucidate the structure of relationships between milk yield, composition, and coagulation properties of Brown Swiss cattle. Multivariate factor analysis was used to derive new synthetic variables that can be used for selection purposes. For this reason, genetic parameters of these new variables were estimated. Individual records on milk yield, fat and protein percentages, casein content, lactose percentage, somatic cell count, titratable acidity, and pH were taken on 1,200 Italian Brown Swiss cows located in 38 herds. Factor analysis was able to extract 4 latent variables with an associated communality equal to 70% of the total original variance. The 4 latent factors were interpreted as indicators of milk composition, coagulation, acidity, and mammary gland health, respectively. Factor scores calculated for each animal exhibited coherent patterns along the lactation and across different parities. Estimation of genetic parameters of factor scores carried out with a multiple-trait Bayesian hierarchical model showed moderate to low heritabilities (raging from 0.10 to 0.23) and genetic correlations (from -0.15 to 0.46). Results of the present study support the hypothesis of a simpler structure that controls, at least in part, the covariance of milk composition and coagulation properties. Moreover, extracted variables may be useful for both breeding and management purposes, being able to represent, with a single value for each animal, complex traits such as milk coagulation properties or health status of the mammary gland.

BackgroundIn the last 50 years, the diversity of cattle breeds has experienced a severe contraction. However, in spite of the growing diffusion of cosmopolite specialized breeds, several local cattle breeds are still farmed in Italy. Genetic characterization of breeds represents an essential step to guide decisions in the management of farm animal genetic resources. The aim of this work was to provide a high-resolution representation of the genome-wide diversity and population structure of Italian local cattle breeds using a medium-density single nucleotide polymorphism (SNP) array.ResultsAfter quality control filtering, the dataset included 31,013 SNPs for 800 samples from 32 breeds. Our results on the genetic diversity of these breeds agree largely with their recorded history. We observed a low level of genetic diversity, which together with the small size of the effective populations, confirmed that several breeds are threatened with extinction. According to the analysis of runs of homozygosity, evidence of recent inbreeding was strong in some local breeds, such as Garfagnina, Mucca Pisana and Pontremolese. Patterns of genetic differentiation, shared ancestry, admixture events, and the phylogenetic tree, all suggest the presence of gene flow, in particular among breeds that originate from the same geographical area, such as the Sicilian breeds. In spite of the complex admixture events that most Italian cattle breeds have experienced, they have preserved distinctive characteristics and can be clearly discriminated, which is probably due to differences in genetic origin, environment, genetic isolation and inbreeding.ConclusionsThis study is the first exhaustive genome-wide analysis of the diversity of Italian cattle breeds. The results are of significant importance because they will help design and implement conservation strategies. Indeed, efforts to maintain genetic diversity in these breeds are needed. Improvement of systems to record and monitor inbreeding in these breeds may contribute to their in situ conservation and, in view of this, the availability of genomic data is a fundamental resource.Electronic supplementary materialThe online version of this article (10.1186/s12711-018-0406-x) contains supplementary material, which is available to authorized users.

Milk composition and its technological properties are traits of interest for the dairy sheep industry because almost all milk produced is processed into cheese. However, several variables define milk technological properties and a complex correlation pattern exists among them. In the present work, we measured milk composition, coagulation properties, and individual cheese yields in a sample of 991 Sarda breed ewes in 47 flocks. The work aimed to study the correlation pattern among measured variables and to obtain new synthetic indicators of milk composition and cheese-making properties. Multivariate factor analysis was carried out on individual measures of milk coagulation parameters; cheese yield; fat, protein, and lactose percentages; somatic cell score; casein percentage; NaCl content; pH; and freezing point. Four factors that were able to explain about 76% of the original variance were extracted. They were clearly interpretable: the first was associated with composition and cheese yield, the second with udder health status, the third with coagulation, and the fourth with curd characteristics. Factor scores were then analyzed by using a mixed linear model that included the fixed effect of parity, lambing month, and lactation stage, and the random effect of flock-test date. The patterns of factor scores along lactation stages were coherent with their technical meaning. A relevant effect of flock-test date was detected, especially on the 2 factors related to milk coagulation properties. Results of the present study suggest the existence of a simpler latent structure that regulates relationships between variables defining milk composition and coagulation properties in sheep. Heritability estimates for the 4 extracted factors were from low to moderate, suggesting potential use of these new variables as breeding goals.