Merino and Merino-derived sheep breeds have been widely distributed across the world, both as purebred and admixed populations. They represent an economically and historically important genetic resource which over time has been used as the basis for the development of new breeds. In order to examine the genetic influence of Merino in the context of a global collection of domestic sheep breeds, we analyzed genotype data that were obtained with the OvineSNP50 BeadChip (Illumina) for 671 individuals from 37 populations, including a subset of breeds from the Sheep HapMap dataset. Based on a multi-dimensional scaling analysis, we highlighted four main clusters in this dataset, which corresponded to wild sheep, mouflon, primitive North European breeds and modern sheep (including Merino), respectively. The neighbor-network analysis further differentiated North-European and Mediterranean domestic breeds, with subclusters of Merino and Merino-derived breeds, other Spanish breeds and other Italian breeds. Model-based clustering, migration analysis and haplotype sharing indicated that genetic exchange occurred between archaic populations and also that a more recent Merino-mediated gene flow to several Merino-derived populations around the world took place. The close relationship between Spanish Merino and other Spanish breeds was consistent with an Iberian origin for the Merino breed, with possible earlier contributions from other Mediterranean stocks. The Merino populations from Australia, New Zealand and China were clearly separated from their European ancestors. We observed a genetic substructuring in the Spanish Merino population, which reflects recent herd management practices. Our data suggest that intensive gene flow, founder effects and geographic isolation are the main factors that determined the genetic makeup of current Merino and Merino-derived breeds. To explain how the current Merino and Merino-derived breeds were obtained, we propose a scenario that includes several consecutive migrations of sheep populations that may serve as working hypotheses for subsequent studies.

Among the European countries, Italy counts the largest number of local goat breeds. Thanks to the recent availability of a medium-density SNP (single nucleotide polymorphism) chip for goat, the genetic diversity of Italian goat populations was characterized by genotyping samples from 14 Italian goat breeds that originate from different geographical areas with more than 50 000 SNPs evenly distributed on the genome. Analysis of the genotyping data revealed high levels of genetic polymorphism and an underlying North–south geographic pattern of genetic diversity that was highlighted by both the first dimension of the multi-dimensional scaling plot and the Neighbour network reconstruction. We observed a moderate and weak population structure in Northern and Central-Southern breeds, respectively, with pairwise FST values between breeds ranging from 0.013 to 0.164 and 7.49 % of the total variance assigned to the between-breed level. Only 2.11 % of the variance explained the clustering of breeds into geographical groups (Northern, Central and Southern Italy and Islands). Our results indicate that the present-day genetic diversity of Italian goat populations was shaped by the combined effects of drift, presence or lack of gene flow and, to some extent, by the consequences of traditional management systems and recent demographic history. Our findings may constitute the starting point for the development of marker-assisted approaches, to better address future breeding and management policies in a species that is particularly relevant for the medium- and long-term sustainability of marginal regions.

Bovine milk is important for human nutrition, but its fat content is often criticized as a risk factor in cardiovascular disease. Selective breeding programs could be used to alter the fatty acid (FA) composition of bovine milk to improve the healthiness of dairy products for human consumption. Here, we performed a genome-wide association study (GWAS) on bovine milk to identify genomic regions or specific genes associated with FA profile and to investigate genetic differences between the Italian Simmental (IS) and Italian Holstein (IH) breeds. To achieve this, we first characterized milk samples from 416 IS cows and 436 IH cows for their fat profile by gas chromatography. Subjects were genotyped with single nucleotide polymorphism array and a single-marker regression model for GWAS was performed. Our findings confirm previously reported quantitative trait loci strongly associated with bovine milk fat composition. More specifically, our GWAS results revealed significant signals on chromosomes Bos taurus autosome 19 and 26 for milk FA. Further analysis using a gene-centric approach and pathway meta-analysis identified not only some well-known genes underlying quantitative trait loci for milk FA components, such as FASN, SCD, and DGAT1, but also other significant candidate genes, including some with functional roles in pathways related to "Lipid metabolism." Highlighted genes related to FA profile include ECI2, PCYT2, DCXR, G6PC3, PYCR1, and ALG12 in IS, and CYP17A1, ACO2, PI4K2A, GOT1, GPT, NT5C2, PDE6G, POLR3H, and COX15 in IH. Overall, the breed-specific association outcomes reflect differences in the genetic backgrounds of the IS and IH breeds and their selective breeding histories.