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.

The aim of this study was to estimate genetic parameters for several female fertility criteria and to choose the most suitable selection index in Spanish Florida and Payoya goat breeds. In this study, we analyzed as fertility traits, the age at first kidding (AgFiKid), and the interval between the first and second kiddings (Int12Kid), between the second, third, and remaining kiddings (Int3toKid), and between all kiddings (IntAllKid) in 51,123 and 22,049 Florida and Payoya females, respectively. Genetic parameters were estimated by fitting animal models using restricted maximum likelihood (REML) methodology. We proposed six selection indices to compare the genetic responses for all traits included, based on a new selection index theory. The heritability and repeatability estimates of the traits were low, as expected. The genetic correlations among fertility traits covered a wide range of values from 0.07 (AgFiKid-Int12Kid) to 0.71 (Int3toKid-IntAllKid) in Florida and from −0.02 (AgFiKid-Int12Kid) to 0.82 (Int3toKid-IntAllKid) in Payoya. Overall, the results of this study indicate that IntAllKid gives the highest genetic responses in both breeds but is expressed late in a female’s life. However, AgFiKid and Int12Kid could be recommended as early selection criteria for female fertility in both breeds.

A total of 1,031,143 records of daily dairy control test of Spanish Florida goats were used for this study. The database was edited, and only the records of the first three lactations were kept. The final database contained 340,654 daily-test somatic cell counts from 27,749 daughters of 941 males and 16,243 goats. The evolution of this count in the last 14 years was analyzed following French and American international associations’ criteria for the risk of mastitis in goats, and confirmed the slight increase in SCS in the last years and the importance of this problem (50% of dairy control tests show a risk of suffering mastitis). For the genetic analysis, the SCS records were log-transformed to normalize this variable. Two strategies were used for the genetic analysis: a univariate animal model for the SCS assuming that SCS does not vary throughout the parities, and a multi-character animal model, where SCS is not considered as the same character in the different parities. The heritabilities (h²) were higher in the multiple traits models, showings an upward trend from the first to the third parity (h² between 0.245 to 0.365). The genetic correlations of the same trait, as well as between breeding values (GVs) between different parities, were different from unity. The breeding values (EBVs) obtained for both models were subjected to a PCA: the first eigenvector (${\lambda }_1)$ explained most of the variations (between 74% to 90%), while the second ${\lambda }_2$accounted for between 9% to 20% of the variance, which shows that the selection will be proportionally favorable but not equivalent in all parities and that there are some variations in the type of response.

We studied the effect of the Temperature Humidity Index (THI) (i.e., the average of temperature and relative humidity registered at meteorological stations) closest to the farms taken during the test day (TD), for total daily protein and fat yields (fpy) of the three main Spanish dairy goats. The data were from Florida (11,244 animals and 126,825 TD), Malagueña (12,215 animals and 141,856 TD) and Murciano Granadina (5162 animals and 62,834 TD) breeding programs and were studied by different linear models to estimate the nature of the fpy response throughout the THI and the weeks of lactation (Days in Milk, DIM) trajectories. The results showed an antagonism between THI and DIM, with a marked depression in the fpy level in animals kept in the hot zone of the THI values (THI > 25) compared with those in the cold zone (THI ≤ 16), with a negative impact equivalent to production of 13 to 30 days. We used a Reaction Norm model (RN), including THI and DIM as fixed covariates and a Test Day Model (TDM), to estimate the genetic (co)variance components. The heritability and genetic correlations estimated with RN and TDM showed a decreased pattern along the scale of THI and DIM, with slight differences between breeds, meaning that there was significant genetic variability in the animal’s ability to react to different levels of THI, which is not constant throughout the DIM, showing the existence of genotype-environment interaction. The breeding values (BV) of all animals for each level of THI and DIM were subject to a principal component analysis, and the results showed that 89 to 98% of the variance between the BV was explained by the two first eigenvalues. The standardized BV were weighted with the corresponding eigenvector coefficients to construct an index that showed, in a single indicator, the most complete expression of the existing genetic variability in the animals’ ability to produce fpy along the trajectories of THI and DIM. This new option will make it easier to select animals which are more productive, and with better adaptability to heat stress, as well as enabling us to exploit genetic variations in the form of the response to heat stress to be adapted to different production systems.

The aim of this study was to compare selection indices for important traits in intensive Spanish goat breeds in four economic scenarios, using the Florida as most representative breed of this production system in Spain. For this analysis, we considered the following traits: milk yield (MY), fat plus protein yields (FPY), casein yield (CY), somatic cell score (SCS), reproductive efficiency (RE), litter size (LS), mammary system (MS), final score (FS), body capacity index (BCI), and length of productive life (LPL). We estimated the genetic parameters and EBVs of most of these traits with REML methodology, while LPL was modeled through survival analysis. Four scenarios were proposed, depending on the overall objective for improvement: (1) milk production, (2) milk production and cheese extract, (3) cheese extract, and (4) milk production, cheese extract and sale of animals. Then, within each scenario, three different types of indices were designed using the different primary and secondary objectives/criteria considered suitable to improve the overall objective. The results indicated that selecting only for primary traits yielded the highest genetic response for all the scenarios. Including secondary traits led to positive correlated responses in those traits, but a decrease in the responses in the primary criteria.