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 (h2) were higher in the multiple traits models, showings an upward trend from the first to the third parity (h2 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 (λ1) explained most of the variations (between 74% to 90%), while the second λ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.

Genetic parameters of litter size (LS) and number of kids born alive (NBA) were estimated across parities in Florida and Payoya goat breeds using a threshold random regression model (RRM). We analysed the reproductive records for the Florida and Payoya breeds separately, and a total of 130.849 and 67.478 reproductive records from the first to the seventh parities from Florida and Payoya, respectively, were included in the analysis. Random regressions on Legendre polynomials of standardised parity were included for permanent environmental and additive genetic effects. We based our estimation of all the covariance components on Bayesian inference under categorical distribution and considered heterogeneous residual variances in the model. The estimates of heritabilities ranged from null to 0.37 for LS and from 0.02 to 0.20 for NBA, while the repeatability estimates were between 0.008-0.60 and 0.03-0.44 for LS and NBA, respectively. Phenotypic correlations within-trait along parities ranged widely from À0.31 to 0.97, while the genetic correlations within-trait through parities showed a wider range and magnitude of values, from À0.98 to 0.99 for LS and from À0.72 to 0.99 for NBA. The estimates of genetic correlations between LS and NBA in the different parities were positive in the first two parities and negative in the later ones. Due to the noticeable changes in variance and the incomplete genetic correlations, the use of a random regression model to analyse LS and NBA is recommended in goat breeds. HIGHLIGHTSA bivariate threshold random regression model has been used for the first time to estimate genetic parameters of litter size and number of kids born alive across parities in Spanish goat breeds. The evolution of genetic parameters along parities is different between Florida and Payoya that have different production systems, being more intensive in the case of the Florida breed. Estimates of variance were not constant across parities for both traits and breeds. Genetic correlations within-trait across parities were different from unity. Thus, LS and NBA should be treated as different traits in successive parities. Given that the estimates of variance and covariances are not constant across parities for both traits and breeds, the use of RRM is highly recommended.

This study aimed to estimate genetic parameters of reproductive efficiency over a wide age range of females using a random regression model (RRM) in Spanish goat breeds. A total of 138,139 and 64,638 reproductive records from the first to the sixth parities from Florida and Payoya, respectively, were included in the analysis. Random regressions on Legendre polynomials of standardised age were included for permanent environmental and additive genetic effects. Estimation of the covariance components was based on the Bayesian inference using the GIBBS3F90 software. Differences among genetic variance components for reproductive efficiency were observed over the animal's lifetime. The estimates of heritabilities were moderate, ranging from 0.21 to 0.32 for Florida and from 0.25 to 0.35 for Payoya, while the fractions of phenotypic variance explained by the permanent environmental effect were high, varying between 0.45-0.68 and 0.58-0.71 for Florida and Payoya, respectively. The correlations for permanent environmental effect over age ranged from 0.37 to 0.99, while the genetic correlations between the different ages varied from 0.36 to 0.98 for Florida and from 0.80 to 0.99 for Payoya. The results from this study support the validity of using an RRM to genetically analyse reproductive efficiency in Spanish dairy goats following the changes in variances and the genetic correlations different from unity over the animal's lifetime. Moreover, reproductive efficiency is a highly heritable trait that is expressed early in a female's life, and it could be used as a precocious selection criterion to improve female fertility in Spanish dairy goats. HIGHLIGHTSReproductive efficiency (RE) has been proposed as a trait to improve female fertility in dairy goats. RE was genetically analysed using an RRM over a range of ages in females from the Florida and Payoya breeds raised under different production systems. Variance components for RE were not constant over age, while the genetic correlations between RE in the different ages were different from unity. The use of RRM for RE is therefore justified in both goat breeds.