Analysis of the genomic landscape of inbreeding in two divergent groups of Spanish Florida goats.

The analysis of the genomic landscape of inbreeding using runs of homozygosity (ROH) patterns is becoming an interesting tool to partially understand phenotypic differences among individuals. In this study, we analysed genome-wide ROH patterns in two groups of Florida goats. We first determined the inbreeding levels of each individual by calculating ROH-based inbreeding coefficients (FROH ). Then, the individuals were divided into two groups based on FROH : high inbreeding (HI, FROH >0.1) and low inbreeding (LI, FROH <0.03). Finally, we performed an extensive in-depth analysis of ROH distribution in each group separately. We found a higher abundance of short ROH in LI, whereas long ROH was more frequent in HI. Furthermore, ROH abundance was not evenly distributed among chromosomes within groups, with some chromosomes showing larger numbers of ROH, like CHI6, CHI7 and CHI27. A different landscape was observed in recent inbreeding (ROH >8 Mb), with significant increases in CHI6, CHI11 and CHI28. Determination of genomic regions with significantly increased ROH (ROH islands-ROHi) showed 13 ROHi related to whole inbreeding and five ROHi associated with recent inbreeding analysis. Within these genomic regions, 123 and 101 genes were identified in HI and LI, respectively, including 10 and seven candidate genes previously related to production, fertility and heat resistance in goats and livestock species.

Genetic Parameters of Somatic Cell Score in Florida Goats Using Single and Multiple Traits Models.

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.

Use of Principal Component Analysis to Combine Genetic Merit for Heat Stress and for Fat and Protein Yield in Spanish Autochthonous Dairy Goat Breeds.

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.

Selection Criteria for Improving Fertility in Spanish Goat Breeds: Estimation of Genetic Parameters and Designing Selection Indices for Optimal Genetic Responses.

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.