Single-step genomic BLUP with many metafounders

Single-step genomic BLUP (ssGBLUP) model for routine genomic prediction of breeding values is developed intensively for many dairy cattle populations. Compatibility between the genomic (G) and the pedigree (A) relationship matrices remains an important challenge required in ssGBLUP. The compatibility relates to the amount of missing pedigree information. There are two prevailing approaches to account for the incomplete pedigree information: unknown parent groups (UPG) and metafounders (MF). unknown parent groups have been used routinely in pedigree-based evaluations to account for the differences in genetic level between groups of animals with missing parents. The MF approach is an extension of the UPG approach. The MF approach defines MF which are related pseudo-individuals. The MF approach needs a Γ matrix of the size number of MF to describe relationships between MF. The UPG and MF can be the same. However, the challenge in the MF approach is the estimation of Γ having many MF, typically needed in dairy cattle. In our study, we present an approach to fit the same amount of MF as UPG in ssGBLUP with Woodbury matrix identity (ssGTBLUP). We used 305-day milk, protein, and fat yield data from the DFS (Denmark, Finland, Sweden) Red Dairy cattle population. The pedigree had more than 6 million animals of which 207,475 were genotyped. We constructed the preliminary gamma matrix (Γpre) with 29 MF which was expanded to 148 MF by a covariance function (Γ148). The quality of the extrapolation of the Γpre matrix was studied by comparing average off-diagonal elements between breed groups. On average relationships among MF in Γ148 were 1.8% higher than in Γpre. The use of Γ148 increased the correlation between the G and A matrices by 0.13 and 0.11 for the diagonal and off-diagonal elements, respectively. [G]EBV were predicted using the ssGTBLUP and Pedigree-BLUP models with the MF and UPG. The prediction reliabilities were slightly higher for the ssGTBLUP model using MF than UPG. The ssGBLUP MF model showed less overprediction compared to other models.

Genomic Estimated Breeding Value of Milk Performance and Fertility Traits in the Russian Black-and-White Cattle Population

A breakthrough in cattle breeding was achieved with the incorporation of animal genomic data into breeding programs. The introduction of genomic selection has a major impact on traditional genetic assessment systems and animal genetic improvement programs. Since 2010, genomic selection has been officially introduced in the evaluation of the breeding and genetic potential of cattle in Europe, the U.S., Canada, and many other developed countries. The purpose of this study is to develop a system for a genomic evaluation of the breeding value of the domestic livestock of Black-and-White and Russian Holstein cattle based on 3 milk performance traits: daily milk yield (kg), daily milk fat (%), and daily milk protein content (%) and 6 fertility traits: age at first calving (AFC), calving interval (CI), calving to first insemination interval (CFI), interval between first and last insemination (IFL), days open (DO), and number of services (NS). We built a unified database of breeding animals from 523 breeding farms in the Russian Federation. The database included pedigree information on 2,551,529 cows and 69,131 bulls of the Russian Holstein and Black-and-White cattle breeds, as well as information on the milk performance of 1,597,426 cows with 4,771,366 completed lactations. The date of birth of the animals included in the database was between 1975 and 2017. Genotyping was performed in 672 animals using a BovineSNP50 v3 DNA Analysis BeadChip microarray (Illumina, USA). The genomic estimated breeding value (GEBV) was evaluated only for 644 animals (427 bulls and 217 cows) using the single-step genomic best linear unbiased prediction – animal model (ssGBLUP-AM). The mean genetic potential was +0.88 and +1.03 kg for the daily milk yield, -0.002% for the milk fat content, and –0.003 and 0.001% for the milk protein content in the cows and bulls, respectively. There was negative genetic progress in the fertility traits in the studied population between 1975 and 2017. The reliability of the estimated breeding value (EBV) for genotyped bulls ranged from 89 to 93% for the milk performance traits and 85 to 90% for the fertility traits, whereas the reliability of the genomic estimated breeding value (GEBV) varied 54 to 64% for the milk traits and 23 to 60% for the fertility traits. This result shows that it is possible to use the genomic estimated breeding value with rather high reliability to evaluate the domestic livestock of Russian Holstein and Black-and-White cattle breeds for fertility and milk performance traits. This system of genomic evaluation may help bring domestic breeding in line with modern competitive practices and estimate the breeding value of cattle at birth based on information on the animal’s genome.

Autoregressive single-step model for genomic evaluation of longitudinal reproductive traits in portuguese holstein cattle

We investigated the applicability of ssGBLUP methodology under the autoregressive model (H-AR) for genomic evaluation of longitudinal reproductive traits in Portuguese Holstein cattle. The genotype data of 1,230 bulls and 1,645 cows were considered in our study. The reproductive traits evaluated were interval from calving to first service (ICF), calving interval (CI) and daughter pregnancy rate (DPR) measured during the first four parities. Reliability and rank correlation were used to compare the H-AR with the traditional pedigree-based autoregressive models (A-AR). In addition, a validation study was performed considering different scenarios. Higher genomic estimated breeding values (GEBV) reliabilities were obtained for genotyped bulls when evaluated under the H-AR model, with emphasis on bulls with less than 9 daughters. For this group, the averages of GEBV reliabilities corresponded to 0.62, 0.69 and 0.62 for ICF, CI and DPR, respectively, while the averages obtained by the A-AR model were 0.27, 0.15 and 0.16. The validation study was favourable to H-AR. The best results were observed in the scenario where genotyped cows were combined with contributing bulls (genotyped bulls with daughter or relationship information in the population). Overall, the results suggest that ssGBLUP methodology under the autoregressive model is a feasible and applicable approach to be used in genomic analyses of longitudinal reproductive traits in Portuguese Holstein cattle.

GWAS and gene networks for milk-related traits from test-day multiple lactations in Portuguese Holstein cattle

This study focused on the identification of QTL regions, candidate genes, and network related genes based on the first 3 lactations (LAC3) of milk, fat, and protein yields, and somatic cell score (SCS) in Portuguese Holstein cattle. Additionally, the results were compared with those from only first lactation (LAC1) data. The analyses were performed using the weighted single-step GWAS under an autoregressive test-day (TD) multiple lactations model. A total of 11,434,294 and 4,725,673 TD records from LAC3 and LAC1, respectively, including 38,323 autosomal SNPs and 1338 genotyped animals were used in GWAS analyses. A total of 51 (milk), 5 (fat), 24 (protein), and 4 (SCS) genes were associated to previously annotated relevant QTL regions for LAC3. The CACNA2D1 at BTA4 explained the highest proportion of genetic variance respectively for milk, fat, and protein yields. For SCS, the TRNAG-CCC at BTA14, MAPK10, and PTPN3 genes, both at BTA6 were considered important candidate genes. The accessed network refined the importance of the reported genes. CACNA2D1 regulates calcium density and activation/inactivation kinetics of calcium transport in the mammary gland; whereas TRNAG-CCC, MAPK10, and PTPN3 are directly involved with inflammatory processes widely derived from mastitis. In conclusion, potential candidate genes (TRNAG-CCC, MAPK10, and PTPN3) associated with somatic cell were highlighted, which further validation studies are needed to clarify its mechanism action in response to mastitis. Moreover, most of the candidate genes identified were present in both (LAC3 and LAC1) for milk, fat and protein yields, except for SCS, in which no candidate genes were shared between LAC3 and LAC1. The larger phenotypic information provided by LAC3 dataset was more effective to identify relevant genes, providing a better understanding of the genetic architecture of these traits over all lactations simultaneously.

Autoregressive repeatability model for genetic evaluation of longitudinal reproductive traits in dairy cattle

We investigated the efficiency of the autoregressive repeatability model (AR) for genetic evaluation of longitudinal reproductive traits in Portuguese Holstein cattle and compared the results with those from the conventional repeatability model (REP). The data set comprised records taken during the first four calving orders, corresponding to a total of 416, 766, 872 and 766 thousand records for interval between calving to first service, days open, calving interval and daughter pregnancy rate, respectively. Both models included fixed (month and age classes associated to each calving order) and random (herd-year-season, animal and permanent environmental) effects. For AR model, a first-order autoregressive (co)variance structure was fitted for the herd-year-season and permanent environmental effects. The AR outperformed the REP model, with lower Akaike Information Criteria, lower Mean Square Error and Akaike Weights close to unity. Rank correlations between estimated breeding values (EBV) with AR and REP models ranged from 0.95 to 0.97 for all studied reproductive traits, when the total bulls were considered. When considering only the top-100 selected bulls, the rank correlation ranged from 0.72 to 0.88. These results indicate that the re-ranking observed at the top level will provide more opportunities for selecting the best bulls. The EBV reliabilities provided by AR model was larger for all traits, but the magnitudes of the annual genetic progress were similar between two models. Overall, the proposed AR model was suitable for genetic evaluations of longitudinal reproductive traits in dairy cattle, outperforming the REP model.