The availability of SNP chips for massive genotyping has proven to be useful to genetically characterize populations of domestic cattle and to assess their degree of divergence. In this study, the Illumina BovineHD BeadChip genotyping array was used to describe the genetic variability and divergence among 7 important autochthonous Spanish beef cattle breeds. The within-breed genetic diversity, measured as the marker expected heterozygosity, was around 0.30, similar to other European cattle breeds. The analysis of molecular variance revealed that 94.22% of the total variance was explained by differences within individuals whereas only 4.46% was the result of differences among populations. The degree of genetic differentiation was small to moderate as the pairwise fixation index of genetic differentiation among breeds (F) estimates ranged from 0.026 to 0.068 and the Nei's D genetic distances ranged from 0.009 to 0.016. A neighbor joining (N-J) phylogenetic tree showed 2 main groups of breeds: Pirenaica, Bruna dels Pirineus, and Rubia Gallega on the one hand and Avileña-Negra Ibérica, Morucha, and Retinta on the other. In turn, Asturiana de los Valles occupied an independent and intermediate position. A principal component analysis (PCA) applied to a distance matrix based on marker identity by state, in which the first 2 axes explained up to 17.3% of the variance, showed a grouping of animals that was similar to the one observed in the N-J tree. Finally, a cluster analysis for ancestries allowed assigning all the individuals to the breed they belong to, although it revealed some degree of admixture among breeds. Our results indicate large within-breed diversity and a low degree of divergence among the autochthonous Spanish beef cattle breeds studied. Both N-J and PCA groupings fit quite well to the ancestral trunks from which the Spanish beef cattle breeds were supposed to derive.

BackgroundThe single-step covariance matrix H combines the pedigree-based relationship matrix with the more accurate information on realized relatedness of genotyped individuals represented by the genomic relationship matrix . In particular, to improve convergence behavior of iterative approaches and to reduce inflation, two weights and have been introduced in the definition of , which blend the inverse of a part of with the inverse of . Since the definition of this blending is based on the equation describing , its impact on the structure of is not obvious. In a joint discussion, we considered the question of the shape of for non-trivial and .ResultsHere, we present the general matrix as a function of these parameters and discuss its structure and properties. Moreover, we screen for optimal values of and with respect to predictive ability, inflation and iterations up to convergence on a well investigated, publicly available wheat data set.ConclusionOur results may help the reader to develop a better understanding for the effects of changes of and on the covariance model. In particular, we give theoretical arguments that as a general tendency, inflation will be reduced by increasing or by decreasing .

Background Procedures for the detection of signatures of selection can be classified according to the source of information they use to reject the null hypothesis of absence of selection. Three main groups of tests can be identified that are based on: (1) the analysis of the site frequency spectrum, (2) the study of the extension of the linkage disequilibrium across the length of the haplotypes that surround the polymorphism, and (3) the differentiation among populations. The aim of this study was to compare the performance of a subset of these procedures by using a dataset on seven Spanish autochthonous beef cattle populations.ResultsAnalysis of the correlations between the logarithms of the statistics that were obtained by 11 tests for detecting signatures of selection at each single nucleotide polymorphism confirmed that they can be clustered into the three main groups mentioned above. A factor analysis summarized the results of the 11 tests into three canonical axes that were each associated with one of the three groups. Moreover, the signatures of selection identified with the first and second groups of tests were shared across populations, whereas those with the third group were more breed-specific. Nevertheless, an enrichment analysis identified the metabolic pathways that were associated with each group; they coincided with canonical axes and were related to immune response, muscle development, protein biosynthesis, skin and pigmentation, glucose metabolism, fat metabolism, embryogenesis and morphology, heart and uterine metabolism, regulation of the hypothalamic–pituitary–thyroid axis, hormonal, cellular cycle, cell signaling and extracellular receptors.ConclusionsWe show that the results of the procedures used to identify signals of selection differed substantially between the three groups of tests. However, they can be classified using a factor analysis. Moreover, each canonical factor that coincided with a group of tests identified different signals of selection, which could be attributed to processes of selection that occurred at different evolutionary times. Nevertheless, the metabolic pathways that were associated with each group of tests were similar, which suggests that the selection events that occurred during the evolutionary history of the populations probably affected the same group of traits.Electronic supplementary materialThe online version of this article (doi:10.1186/s12711-016-0258-1) contains supplementary material, which is available to authorized users.