Search

Disease-related milk losses directly affect dairy herds' profitability and the production efficiency of the dairy industry. Therefore, this study aimed to quantify phenotypic variability in milk fluctuation periods related to diseases and to explore milk fluctuation traits as indicators of disease resilience. By combining high-frequency daily milk yield data with disease records of cows that were treated and recovered from the disease, we estimated milk variability trends within a fixed period around the treatment day of each record for 5 diseases: udder health, reproductive disorders, metabolic disorders, digestive disorders, and hoof health. The average milk yield decreased rapidly from 6 to 8 d before the treatment day for all diseases, with the largest milk reduction observed on the treatment day. Additionally, we assessed the significance of milk fluctuation periods highly related to diseases by defining milk fluctuations as a period of at least 10 consecutive days in which milk yield fell below 90% of the expected milk production values at least once. We defined the development and recovery phases of milk fluctuations using 3,847 milk fluctuation periods related to disease incidences, and estimated genetic parameters of milk fluctuation traits, including milk losses, duration of the fluctuation, variation rate in daily milk yield, and standard deviation of milk deviations for each phase and their genetic correlation with several important traits. In general, the disease-related milk fluctuation periods lasted 21.19 ± 10.36 d with a milk loss of 115.54 ± 92.49 kg per lactation. Compared with the development phase, the recovery phase lasted an average of 3.3 d longer, in which cows produced 11.04 kg less milk and exhibited a slower variation rate in daily milk yield of 0.35 kg/d. There were notable differences in milk fluctuation traits depending on the disease, and greater milk losses were observed when multiple diseases occurred simultaneously. All milk fluctuation traits evaluated were heritable with heritability estimates ranging from 0.01 to 0.10, and moderate to high genetic correlations with milk yield (0.34 to 0.64), milk loss throughout the lactation (0.22 to 0.97), and resilience indicator (0.39 to 0.95). These results indicate that cows with lower milk losses and higher resilience tend to have more stable milk fluctuations, which supports the potential for breeding for more disease-resilient cows based on milk fluctuation traits. Overall, this study confirms the high effect of diseases on milk yield variability and provides insightful information about their relationship with relevant traits in Holstein cattle. Furthermore, this study shows the potential of using high-frequency automatic monitoring of milk yield to assist on breeding practices and health management in dairy cows. [ABSTRACT FROM AUTHOR]

The elite bull plays an extremely important role in the genetic progression of the dairy cow population. The previous results indicated the potential positive relationship of large scrotal circumference (SC) with improved semen volume, concentration, and motility. In order to improve bull's semen quantity and quality by selection, it is necessary to estimate the genetic parameters of semen traits and their correlations with other conformation traits such as SC that could be used for an indirect selection. In this study, the genetic parameters of seven semen traits (n = 66,260) and nine conformation traits (n = 3,642) of Holstein bulls (n = 453) were estimated by using the bivariate repeatability animal model with the average information-restricted maximum likelihood (AI-REML) approach. The results showed that the estimated heritabilities of semen traits ranged from 0.06 (total number of motile sperm, TNMS) to 0.37 (percentage of abnormal sperm, PAS) and conformation traits ranged from 0.23 (pin width, PW) to 0.69 (hip height, HH). The highest genetic correlations were found between semen volume per ejaculation (SVPE), semen concentration per ejaculation (SCPE), total number of sperm (TNS), and TNMS traits that were 0.97, 0.98, 1.00, and 0.99, respectively. Phenotypic correlations between SC and SVPE, SCPE, TNS, and TNMS were 0.35, 0.35, 0.48, and 0.42, respectively. In summary, the moderate or high heritability of semen traits indicates that genetic improvement of semen quality by selection is feasible, where SC could be a useful trait for indirect selection or as correlated information to improve semen quantity and production in the practical bull breeding programs. [ABSTRACT FROM AUTHOR]

Reproductive traits of dairy cattle are bound to the actual efficiency of farm operation, which therefore show great economic importance. Among them, some traits were deemed to be simultaneously affected by service sire and mating cow. Service sires are proved to play an important role in reproduction process of cows. However, limited study explored the genetic effect of service sire (GESS), let alone the genomic prediction of this effect. In the present study, 2244 genotyped bulls together with phenotypic records were used to predict the GESS on conception rate, 56‐day non‐return rate, calving ease, stillbirth and gestation length. The feasibilities of multi‐step genomic best linear unbiased predictor (msGBLUP) and single‐step genomic best linear unbiased predictor (ssGBLUP) were investigated under different scenarios, that is, different marker densities and validation population. The predictive accuracies and unbiasedness for GESS ranged from 0.159 to 0.647 and from 0.202 to 2.018, respectively, when validated on young bulls, while the accuracies and unbiasedness ranged from 0.409 to 0.802 and 0.333 to 1.146 when validated on random split data sets. It is feasible to predict GESS on reproductive traits by using a linear mixed model and genomic data, and high‐density marker panel had limited contribution to the prediction. This research investigated the potential factors that influence the genomic prediction of GESS on reproductive traits and indicated the possibility of genomic selection on GESS, both in ideal and practical circumstances. [ABSTRACT FROM AUTHOR]

This study investigated the reliability of genomic prediction (GP) using breed origin of alleles (BOA) approach in the Nordic Red (RDC) population, which has an admixed population structure. The RDC population consists of animals with varying degrees of genetic materials from the Danish Red (RDM), Swedish Red (SRB), Finnish Ayrshire (FAY), and Holstein (HOL) because bulls have been used across the breeds. The BOA approach was tested using 39,550 RDC animals in the reference population and 11,786 in the validation population. Deregressed proofs (DRP) of milk, fat and protein were used as response variable for GP. Direct genomic breeding values (DGV) for animals in the validation population were calculated with (BOA model) or without (joint model) considering breed origin of alleles. The joint model assumed homogeneous marker effects and a single set of marker effects were estimated, whereas BOA model assumed heterogeneous marker effects, and different sets of marker effects were estimated across the breeds. For the BOA approach, we tested scenarios assuming both correlated (BOA_cor) and uncorrelated (BOA_uncor) marker effects between the breeds. Additionally, we investigated GP using a standard Illumina 50K chip and including SNP selected from imputed whole-genome sequencing (50K+WGS). We also studied the effect of estimating (co)variances for genome regions of different sizes to exploit the information of the genome regions contributing to the (co)variance between the breeds. Region sizes were set as 1 SNP, a group of 30 or 100 adjacent SNP, or the whole genome. Reliability of DGV was measured as squared correlations between DGV and DRP divided by the reliability of DRP. Across the 3 traits, in general, RS30 and RS100 SNP yielded the highest reliabilities. Including WGS SNP improved reliabilities in almost all scenarios (0.297 on average for 50K and 0.307 on average for 50K+WGS). The BOA_uncor (0.233 on average) was inferior to the joint model (0.339 on average), but the reliabilities obtained using BOA_cor (0.334 on average) in most cases were not significantly different from those obtained using the joint model. The results indicate that both including additional whole-genome sequencing SNP and dividing the genome into fixed regions improve GP in the RDC. The BOA models have the potential to increase the reliability of GP, but the benefit is limited in populations with a high exchange of genetic material for a long time, as is the case for RDC. [ABSTRACT FROM AUTHOR]

Background: In dairy cattle, genomic selection has been implemented successfully for purebred populations, but, to date, genomic estimated breeding values (GEBV) for crossbred cows are rarely available, although they are valuable for rotational crossbreeding schemes that are promoted as efficient strategies. An attractive approach to provide GEBV for crossbreds is to use estimated marker effects from the genetic evaluation of purebreds. The effects of each marker allele in crossbreds can depend on the breed of origin of the allele (BOA), thus applying marker effects based on BOA could result in more accurate GEBV than applying only proportional contribution of the purebreds. Application of BOA models in rotational crossbreeding requires methods for detecting BOA, but the existing methods have not been developed for rotational crossbreeding. Therefore, the aims of this study were to develop and test methods for detecting BOA in a rotational crossbreeding system, and to investigate methods for calculating GEBV for crossbred cows using estimated marker effects from purebreds. Results: For detecting BOA in crossbred cows from rotational crossbreeding for which pedigree is recorded, we developed the AllOr method based on the comparison of haplotypes in overlapping windows. To calculate the GEBV of crossbred cows, two models were compared: a BOA model where marker effects estimated from purebreds are combined based on the detected BOA; and a breed proportion model where marker effects are combined based on estimated breed proportions. The methods were tested on simulated data that mimic the first four generations of rotational crossbreeding between Holstein, Jersey and Red Dairy Cattle. The AllOr method detected BOA correctly for 99.6% of the marker alleles across the four crossbred generations. The reliability of GEBV was higher with the BOA model than with the breed proportion model for the four generations of crossbreeding, with the largest difference observed in the first generation. Conclusions: In rotational crossbreeding for which pedigree is recorded, BOA can be accurately detected using the AllOr method. Combining marker effects estimated from purebreds to predict the breeding value of crossbreds based on BOA is a promising approach to provide GEBV for crossbred dairy cows. [ABSTRACT FROM AUTHOR]

This study investigated effects of integrating single-nucleotide polymorphisms (SNPs) selected based on previous genome-wide association studies (GWASs), from imputed whole-genome sequencing (WGS) data, in the conventional 54K chip on genomic prediction reliability of young stock survival (YSS) traits in dairy cattle. The WGS SNPs included two groups of SNP sets that were selected based on GWAS in the Danish Holstein for YSS index (YSS_SNPs, n = 98) and SNPs chosen as peaks of quantitative trait loci for the traits of Nordic total merit index in Denmark–Finland–Sweden dairy cattle populations (DFS_SNPs, n = 1,541). Additionally, the study also investigated the possibility of improving genomic prediction reliability for survival traits by modeling the SNPs within recessive lethal haplotypes (LET_SNP, n = 130) detected from the 54K chip in the Nordic Holstein. De-regressed proofs (DRPs) were obtained from 6,558 Danish Holstein bulls genotyped with either 54K chip or customized LD chip that includes SNPs in the standard LD chip and some of the selected WGS SNPs. The chip data were subsequently imputed to 54K SNP together with the selected WGS SNPs. Genomic best linear unbiased prediction (GBLUP) models were implemented to predict breeding values through either pooling the 54K and selected WGS SNPs together as one genetic component (a one-component model) or considering 54K SNPs and selected WGS SNPs as two separate genetic components (a two-component model). Across all the traits, inclusion of each of the selected WGS SNP sets led to negligible improvements in prediction accuracies (0.17 percentage points on average) compared to prediction using only 54K. Similarly, marginal improvement in prediction reliability was obtained when all the selected WGS SNPs were included (0.22 percentage points). No further improvement in prediction reliability was observed when considering random regression on genotype code of recessive lethal alleles in the model including both groups of the WGS SNPs. Additionally, there was no difference in prediction reliability from integrating the selected WGS SNP sets through the two-component model compared to the one-component GBLUP. [ABSTRACT FROM AUTHOR]

Icelandic Cattle is a local dairy cattle breed in Iceland. With about 26,000 breeding females, it is by far the largest among the indigenous Nordic cattle breeds. The objective of this study was to investigate the feasibility of genomic selection in Icelandic Cattle. Pedigree-based best linear unbiased prediction (PBLUP) and single-step genomic best linear unbiased prediction (ssGBLUP) were compared. Accuracy, bias, and dispersion of estimated breeding values (EBV) for milk yield (MY), fat yield (FY), protein yield (PY), and somatic cell score (SCS) were estimated in a cross validation-based design. Accuracy (⁠ r ^ ⁠) was estimated by the correlation between EBV and corrected phenotype in a validation set. The accuracy (⁠ r ^ ⁠) of predictions using ssGBLUP increased by 13, 23, 19, and 20 percentage points for MY, FY, PY, and SCS for genotyped animals, compared with PBLUP. The accuracy of nongenotyped animals was not improved for MY and PY, but increased by 0.9 and 3.5 percentage points for FY and SCS. We used the linear regression (LR) method to quantify relative improvements in accuracy, bias (⁠ Δ ^ ⁠), and dispersion (⁠ b ^ ⁠) of EBV. Using the LR method, the relative improvements in accuracy of validation from PBLUP to ssGBLUP were 43%, 60%, 50%, and 48% for genotyped animals for MY, FY, PY, and SCS. Single-step GBLUP EBV were less underestimated (⁠ Δ ^ ⁠), and less overdispersed (⁠ b ^ ⁠) than PBLUP EBV for FY and PY. Pedigree-based BLUP EBV were close to unbiased for MY and SCS. Single-step GBLUP underestimated MY EBV but overestimated SCS EBV. Based on the average accuracy of 0.45 for ssGBLUP EBV obtained in this study, selection intensities according to the breeding scheme of Icelandic Cattle, and assuming a generation interval of 2.0 yr for sires of bulls, sires of dams and dams of bulls, genetic gain in Icelandic Cattle could be increased by about 50% relative to the current breeding scheme. [ABSTRACT FROM AUTHOR]

Background: In dairy cattle populations in which crossbreeding has been used, animals show some level of diversity in their origins. In rotational crossbreeding, for instance, crossbred dams are mated with purebred sires from different pure breeds, and the genetic composition of crossbred animals is an admixture of the breeds included in the rotation. How to use the data of such individuals in genomic evaluations is still an open question. In this study, we aimed at providing methodologies for the use of data from crossbred individuals with an admixed genetic background together with data from multiple pure breeds, for the purpose of genomic evaluations for both purebred and crossbred animals. A three-breed rotational crossbreeding system was mimicked using simulations based on animals genotyped with the 50 K single nucleotide polymorphism (SNP) chip. Results: For purebred populations, within-breed genomic predictions generally led to higher accuracies than those from multi-breed predictions using combined data of pure breeds. Adding admixed population's (MIX) data to the combined pure breed data considering MIX as a different breed led to higher accuracies. When prediction models were able to account for breed origin of alleles, accuracies were generally higher than those from combining all available data, depending on the correlation of quantitative trait loci (QTL) effects between the breeds. Accuracies varied when using SNP effects from any of the pure breeds to predict the breeding values of MIX. Using those breed-specific SNP effects that were estimated separately in each pure breed, while accounting for breed origin of alleles for the selection candidates of MIX, generally improved the accuracies. Models that are able to accommodate MIX data with the breed origin of alleles approach generally led to higher accuracies than models without breed origin of alleles, depending on the correlation of QTL effects between the breeds. Conclusions: Combining all available data, pure breeds' and admixed population's data, in a multi-breed reference population is beneficial for the estimation of breeding values for pure breeds with a small reference population. For MIX, such an approach can lead to higher accuracies than considering breed origin of alleles for the selection candidates, and using breed-specific SNP effects estimated separately in each pure breed. Including MIX data in the reference population of multiple breeds by considering the breed origin of alleles, accuracies can be further improved. Our findings are relevant for breeding programs in which crossbreeding is systematically applied, and also for populations that involve different subpopulations and between which exchange of genetic material is routine practice. [ABSTRACT FROM AUTHOR]

Due to the complexity of longevity trait in dairy cattle, two groups of trait definitions are widely used to measure longevity, either covering the full lifespan or representing only a part of it to achieve an early selection. Usually, only one group of longevity definition is used in breeding program for one population, and genetic studies on the comparisons of two groups of trait definitions are scarce. Based on the data of eight traits well representing the both groups of trait definitions, the current study investigated genetic parameters and genetic architectures of longevity in Holsteins. Heritabilities and correlations of eight longevity traits were estimated using single-trait and multi-trait animal models, with the data from 103,479 cows. Among the cows with phenotypes, 2,630 cows were genotyped with the 150K-SNP panel. A single-trait fixed and random Circuitous Probability Unification model was performed to detect candidate genes for eight longevity traits. Generally, all eight longevity traits had low heritabilities, ranging from 0.038 for total productive life and herd life to 0.090 for days from the first calving to the end of first lactation or culling. High genetic correlations were observed among the traits within the same definition group: from 0.946 to 0.997 for three traits reflecting full lifespan and from 0.666 to 0.997 for five traits reflecting partial productive life. Genetic correlations between two groups of traits ranged from 0.648 to 0.963, and increased gradually with the extension of lactations number regarding the partial productive life traits. A total of 55 SNPs located on 25 chromosomes were found genome-wide significantly associated with longevity, in which 12 SNPs were associated with more than one trait, even across traits of different definition groups. This is the first study to investigate the genetic architecture of longevity representing both full and the partial lifespan simultaneously, which will assist the selection of an appropriate trait definition for genetic improvement of longevity. Because of high genetic correlations with the full lifespan traits and higher heritability, the partial productive life trait measured as the days from the first calving to the end of the third lactation or culling could be a good alternative for early selection on longevity. The candidate genes identified by this study, such as RPRM, GRIA3, GTF2H5, CA5A, CACNA2D1, FGF10, and DNAJA3, could be used to pinpoint causative mutations for longevity and further benefit the genomic improvement of longevity in dairy cattle. [ABSTRACT FROM AUTHOR]

Background: A better understanding of the genetic architecture underlying complex traits (e.g., the distribution of causal variants and their effects) may aid in the genomic prediction. Here, we hypothesized that the genomic variants of complex traits might be enriched in a subset of genomic regions defined by genes grouped on the basis of "Gene Ontology" (GO), and that incorporating this independent biological information into genomic prediction models might improve their predictive ability. Results: Four complex traits (i.e., milk, fat and protein yields, and mastitis) together with imputed sequence variants in Holstein (HOL) and Jersey (JER) cattle were analysed. We first carried out a post-GWAS analysis in a HOL training population to assess the degree of enrichment of the association signals in the gene regions defined by each GO term. We then extended the genomic best linear unbiased prediction model (GBLUP) to a genomic feature BLUP (GFBLUP) model, including an additional genomic effect quantifying the joint effect of a group of variants located in a genomic feature. The GBLUP model using a single random effect assumes that all genomic variants contribute to the genomic relationship equally, whereas GFBLUP attributes different weights to the individual genomic relationships in the prediction equation based on the estimated genomic parameters. Our results demonstrate that the immune-relevant GO terms were more associated with mastitis than milk production, and several biologically meaningful GO terms improved the prediction accuracy with GFBLUP for the four traits, as compared with GBLUP. The improvement of the genomic prediction between breeds (the average increase across the four traits was 0.161) was more apparent than that it was within the HOL (the average increase across the four traits was 0.020). Conclusions: Our genomic feature modelling approaches provide a framework to simultaneously explore the genetic architecture and genomic prediction of complex traits by taking advantage of independent biological knowledge. [ABSTRACT FROM AUTHOR]

For genomic prediction of crossbred animals, models that account for the breed origin of alleles (BOA) in marker genotypes can allow the effects of marker alleles to differ depending on their ancestral breed. Previous studies have shown that genomic estimated breeding values for crossbred cows can be calculated using the marker effects that are estimated in the contributing pure breeds and combined based on estimated BOA in the genotypes of the crossbred cows. In the presented study, we further exploit the BOA information for improving the prediction of genomic breeding values of crossbred dairy cows. We investigated 2 types of BOA-derived breed proportions: global breed proportions, defined as the proportion of marker alleles assigned to each breed across the whole genome; and local breed proportions (LBP), defined as the proportions of alleles on chromosome segments which were assigned to each breed. Further, we investigated 2 BOA-derived measures of heterozygosity for the prediction of total genetic value. First, global breed heterozygosity, defined as the proportion of marker loci that have alleles originating in 2 different breeds over the whole genome. Second, local breed heterozygosity (LBH), defined as proportions of marker loci on chromosome segments that had alleles originating in 2 different breeds. We estimated variance related to LBP and LBH on the remaining variation after accounting for prediction with solutions from the genomic evaluations of the pure breeds and validated alternative models for production traits in 5,214 Danish crossbred dairy cows. The estimated LBP variances were 0.9, 1.2, and 1.0% of phenotypic variance for milk, fat, and protein yield, respectively. We observed no clear LBH effect. Cross-validation showed that models with LBP effects had a numerically small but statistically significantly higher predictive ability than models only including global breed proportions. We observed similar improvement in accuracy by the model having an across crossbred residual additive genetic effect, accounting for the additive genetic variation that was not accounted for by the solutions from purebred. For genomic predictions of crossbred animals, estimated BOA can give useful information on breed proportions, both globally in the genome and locally in genome regions, and on breed heterozygosity. [ABSTRACT FROM AUTHOR]

Genomic predictions have been applied for dairy cattle for more than a decade with great success, but genomic estimated breeding values (GEBV) are not widely available for crossbred dairy cows. The large reference populations already in place for genomic evaluations of many pure breeds makes it interesting to use the accurate solutions, in particular the estimated marker effects, from these evaluations for calculation of GEBV for crossbred heifers and cows. Effects of marker alleles in crossbred animals can depend on breed origin of the alleles (BOA). Therefore, our aim was to investigate if reliable GEBV for crossbred dairy cows can be obtained by combining estimated marker effects from purebred evaluations based on BOA. We used data on 5,467 Danish crossbred dairy cows with contributions from Holstein, Jersey, and Red Dairy Cattle breeds. We assessed BOA assignment on their genotypes and found that we could assign 99.3% of the alleles to a definite breed of origin. We compared GEBV for 2 traits, protein yield and interval between first and last insemination of cows, with 2 models that both combine estimated marker effects from the genomic evaluations of the pure breeds: a breed of origin model that accounts for BOA and a breed proportion model that only accounts for genomic breed proportions in the crossbred animals. We accounted for the difference in level between the purebred evaluations by including intercepts in the models based on phenotypic averages. The predictive ability for protein yield was significantly higher from the breed of origin model, 0.45 compared with 0.43 from the breed proportion model. Furthermore, for the breed proportion model, the GEBVs had level bias, which made comparison across groups with different breed composition skewed. We therefore concluded that reliable genomic predictions for crossbred dairy cows can be obtained by combining estimated marker effects from the genomic evaluations of purebreds using a model that accounts for BOA. [ABSTRACT FROM AUTHOR]