The main goal of the current work was to infer the demographic history of seven Spanish goat breeds (Malagueña, Murciano-Granadina, Florida, Palmera, Mallorquina, Bermeya and Blanca de Rasquera) based on genome-wide diversity data generated with the Illumina Goat SNP50 BeadChip (population size, N = 176). Five additional populations from Europe (Saanen and Carpathian) and Africa (Tunisian, Djallonké and Sahel) were also included in this analysis (N = 80) for comparative purposes. Our results show that the genetic background of Spanish goats traces back mainly to European breeds although signs of North African admixture were detected in two Andalusian breeds (Malagueña and Murciano-Granadina). In general, observed and expected heterozygosities were quite similar across the seven Spanish goat breeds under analysis irrespective of their population size and conservation status. For the Mallorquina and Blanca de Rasquera breeds, which have suffered strong population declines during the past decades, we observed increased frequencies of large-sized (ROH), a finding that is consistent with recent inbreeding. In contrast, a substantial part of the genome of the Palmera goat breed comprised short ROH, which suggests a strong and ancient founder effect. Admixture with African goats, genetic drift and inbreeding have had different effects across the seven Spanish goat breeds analysed in the current work. This has generated distinct patterns of genome-wide diversity that provide new clues about the demographic history of these populations.
Background: Inbreeding is caused by mating between related individuals and is associated with reduced fitness and performance (inbreeding depression). Several studies have detected heterogeneity in inbreeding depression among founder individuals. Recently, a procedure was developed to predict hidden inbreeding depression load that is associated with founders using the Mendelian sampling of non-founders. The objectives of this study were to: (1) analyse the population structure and general inbreeding, and (2) test this recent approach for predicting hidden inbreeding depression load for four morphological traits and two morphology defects in the Pura Raza Española (PRE) horse breed. Results: The regression coefficients that were calculated between trait performances and inbreeding coefficients demonstrated the existence of inbreeding depression. In total, 58,772,533 partial inbreeding coefficients ( Fij ) were estimated for the whole PRE population (328,706 horses). We selected the descendants of horses with a Fij ≥ 6.25% that contributed to at least four offspring and for which morphological traits were measured for the subsequent analysis of inbreeding depression load (639 horses). A pedigree was generated with the last five generations (5026 animals) used as the reference population (average inbreeding coefficient of 8.39% and average relatedness coefficient of 10.76%). Heritability estimates ranged from 0.08 (cresty neck) to 0.80 (height at withers), whereas inbreeding depression load ratios ranged from 0.01 (knock knee) to 0.40 (length of shoulder), for an inbreeding coefficient of 10%. Most of the correlations between additive and inbreeding depression load genetic values and correlations between inbreeding depression load genetic values for the different traits were positive or near 0. Conclusions: Although the average inbreeding depression loads presented negative values, a certain percentage of the animals showed neutral or even positive values. Thus, high levels of inbreeding do not always lead to a decrease in mean phenotypic value or an increase in morphological defects. Hence, individual inbreeding depression loads could be used as a tool to select the most appropriate breeding animals. The possibility of selecting horses that have a high genetic value and are more resistant to the deleterious effects of inbreeding should help improve selection outcomes.
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. Analysis 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. We 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.
Fertility is one of the key factors in the economic and productive success of the equine industry. Despite this, studies on the genetic causes affecting reproductive performance are scarce, especially in mares, where the genetic architecture of the reproductive traits is extremely complex. Today, with the increasing availability of new genomic methodologies for this species, we are presented with an interesting opportunity to understand the genetic basis of equine reproductive disorders. These include, among others, novel techniques for detecting chromosomal abnormalities, whose association with infertility in horses was established over 50 years ago; new sequencing technologies permitting an accurate detection of point mutations influencing fertility, as well as the study of inbreeding and molecular homozygosity, which has been widely suggested as one of the main causes of low reproductive performance in horses. Finally, over the last few years, reproductive performance has also been associated with copy number variants and candidate genes detected by genome-wide association studies on fertility traits. However, such studies are still scarce, probably because they depend on the existence of large and accurate phenotypic datasets of reproductive and/or fertility traits, which are still difficult to obtain in equines.
The aim of this work was to analyze genetic parameters such as the inbreeding coefficient (F), relatedness coefficient (AR) and partial inbreeding coefficient (Fij) of the whole PRE population, and the ancestors which account for 50% of the total genetic variability of the current population, from genealogical information. The average F of the whole PRE population (328,706 animals) has decreased from 8.45% to 7.51% in the least 20 years. The Fij was estimated for the whole PRE population, resulting in a database of 58,772,533 records containing one record for each Fij that each animal receives from a certain common ancestor (CA). A total of 10,244 CAs contributed to the Fij with an average of 5370 descendants, with each descendant having an average of 170 CAs. Over the generations, the number of CAs has increased, while the proportion of Fij by each one has decreased. In addition, the contributions of the more influential ancestors have changed. The increased census, the limited use of artificial insemination and our increased knowledge about inbreeding depression and the animals’ breeding values allow breeders to select horses more for their functionality and conformation than for their pedigree reputation, which is the cause of all these changes.