In dairy cattle, the widespread use of artificial insemination has resulted in increased selection intensity, which has led to spectacular increase in productivity. However, cow fertility has concomitantly severely declined. It is generally assumed that this reduction is primarily due to the negative energy balance of high-producing cows at the peak of lactation. We herein describe the fine-mapping of a major fertility QTL in Nordic Red cattle, and identify a 660-kb deletion encompassing four genes as the causative variant. We show that the deletion is a recessive embryonically lethal mutation. This probably results from the loss of RNASEH2B, which is known to cause embryonic death in mice. Despite its dramatic effect on fertility, 13%, 23% and 32% of the animals carry the deletion in Danish, Swedish and Finnish Red Cattle, respectively. To explain this, we searched for favorable effects on other traits and found that the deletion has strong positive effects on milk yield. This study demonstrates that embryonic lethal mutations account for a non-negligible fraction of the decline in fertility of domestic cattle, and that associated positive effects on milk yield may account for part of the negative genetic correlation. Our study adds to the evidence that structural variants contribute to animal phenotypic variation, and that balancing selection might be more common in livestock species than previously appreciated. We report the identification of a large deletion encompassing four genes and the demonstration of its negative effect on fertility in Nordic Red dairy cattle. We show that this deletion is recessively lethal (homozygous embryos die) and therefore, when carrier cows are mated to carrier bulls, there is a high risk of embryonic mortality. As a result, chances of insemination failure are higher for such matings. Surprisingly, despite its negative effect, the deletion is frequent in Nordic Red cattle. We show that this high frequency may be a consequence of the fact that the deletion is associated with increased milk production and therefore selected for. Due to increased levels of inbreeding resulting from the widespread use of artificial insemination, such recessive lethal alleles may account for a non-negligible fraction of the reduction in fertility observed in cattle.

This study investigated genomic prediction using medium-density (∼54,000; 54K) and high-density marker panels (∼777,000; 777K), based on data from Nordic Holstein and Red Dairy Cattle (RDC). The Holstein data comprised 4,539 progeny-tested bulls, and the RDC data 4,403 progeny-tested bulls. The data were divided into reference data and test data using October 1, 2001, as a cut-off date (birth date of the bulls). This resulted in about 25% genotyped bulls in the Holstein test data and 20% in the RDC test data. For each breed, 3 data sets of markers were used to predict breeding values: (1) 54K data set with missing genotypes, (2) 54K data set where missing genotypes were imputed, and (3) imputed high-density (HD) marker data set created by imputing the 54K data to the HD data based on 557 bulls genotyped using a 777K single nucleotide polymorphism chip in Holstein, and 706 bulls in RDC. Based on the 3 marker data sets, direct genomic breeding values (DGV) for protein, fertility, and udder health were predicted using a genomic BLUP model (GBLUP) and a Bayesian mixture model with 2 normal distributions. Reliability of DGV was measured as squared correlations between deregressed proofs (DRP) and DGV corrected for reliability of DRP. Unbiasedness was assessed by regression of DRP on DGV, based on the bulls in the test data sets. Averaged over the 3 traits, reliability of DGV based on the HD markers was 0.5% higher than that based on the 54K data in Holstein, and 1.0% higher than that in RDC. In addition, the HD markers led to an improvement of unbiasedness of DGV. The Bayesian mixture model led to 0.5% higher reliability than the GBLUP model in Holstein, but not in RDC. Imputing missing genotypes in the 54K marker data did not improve genomic predictions for most of the traits.

This study investigated the accuracy of direct genomic breeding values (DGV) using a genomic BLUP model, genomic enhanced breeding values (GEBV) using a one-step blending approach, and GEBV using a selection index blending approach for 15 traits of Nordic Red Cattle. The data comprised 6,631 bulls of which 4,408 bulls were genotyped using Illumina Bovine SNP50 BeadChip (Illumina, San Diego, CA). To validate reliability of genomic predictions, about 20% of the youngest genotyped bulls were taken as test data set. Deregressed proofs (DRP) were used as response variables for genomic predictions. Reliabilities of genomic predictions in the validation analyses were measured as squared correlations between DRP and genomic predictions corrected for reliability of DRP, based on the bulls in the test data sets. A set of weighting (scaling) factors was used to construct the combined relationship matrix among genotyped and nongenotyped bulls for one-step blending, and to scale DGV and its expected reliability in the selection index blending. Weighting (scaling) factors had a small influence on reliabilities of GEBV, but a large influence on the variation of GEBV. Based on the validation analyses, averaged over the 15 traits, the reliability of DGV for bulls without daughter records was 11.0 percentage points higher than the reliability of conventional pedigree index. Further gain of 0.9 percentage points was achieved by combining information from conventional pedigree index using the selection index blending, and gain of 1.3 percentage points was achieved by combining information of genotyped and nongenotyped bulls simultaneously applying the one-step blending. These results indicate that genomic selection can greatly improve the accuracy of preselection for young bulls in Nordic Red population, and the one-step blending approach is a good alternative to predict GEBV in practical genetic evaluation program.

The objectives of this study were to evaluate the feasibility of use of the test-day (TD) single-step genomic BLUP (ssGBLUP) using phenotypic records of Nordic Red Dairy cows. The critical point in ssGBLUP is how genomically derived relationships (G) are integrated with population-based pedigree relationships (A) into a combined relationship matrix (H). Therefore, we also tested how different weights for genomic and pedigree relationships affect ssGBLUP, validation reliability, and validation regression coefficients. Deregressed proofs for 305-d milk, protein, and fat yields were used for a posteriori validation. The results showed that the use of phenotypic TD records in ssGBLUP is feasible. Moreover, the TD ssGBLUP model gave considerably higher validation reliabilities and validation regression coefficients than the TD model without genomic information. No significant differences were found in validation reliability between the different TD ssGBLUP models according to bootstrap confidence intervals. However, the degree of inflation in genomic enhanced breeding values is affected by the method used in construction of the H matrix. The results showed that ssGBLUP provides a good alternative to the currently used multi-step approach but there is a great need to find the best option to combine pedigree and genomic information in the genomic matrix.

The availability of whole genome sequencing (WGS) data enables the discovery of causative single nucleotide polymorphisms (SNPs) or SNPs in high linkage disequilibrium with causative SNPs. This study investigated effects of integrating SNPs selected from imputed WGS data into the data of 54K chip on genomic prediction in Danish Jersey. The WGS SNPs, mainly including peaks of quantitative trait loci, structure variants, regulatory regions of genes, and SNPs within genes with strong effects predicted with variant effect predictor, were selected in previous analyses for dairy breeds in Denmark–Finland–Sweden (DFS) and France (FRA). Animals genotyped with 54K chip, standard LD chip, and customized LD chip which covered selected WGS SNPs and SNPs in the standard LD chip, were imputed to 54K together with DFS and FRA SNPs. Genomic best linear unbiased prediction (GBLUP) and Bayesian four-distribution mixture models considering 54K and selected WGS SNPs as one (a one-component model) or two separate genetic components (a two-component model) were used to predict breeding values. For milk production traits and mastitis, both DFS (0.025) and FRA (0.029) sets of additional WGS SNPs improved reliabilities, and inclusions of all selected WGS SNPs generally achieved highest improvements of reliabilities (0.034). A Bayesian four-distribution model yielded higher reliabilities than a GBLUP model for milk and protein, but extra gains in reliabilities from using selected WGS SNPs were smaller for a Bayesian four-distribution model than a GBLUP model. Generally, no significant difference was observed between one-component and two-component models, except for using GBLUP models for milk.

The aim of this study was to explore the possibilities of using body condition score (BCS) or dairy character (DC) as indicators of mastitis and diseases other than mastitis in first-parity Danish Holsteins. The dataset included 28,948 observations on conformation scores and 365,136 disease observations. The analysis was performed using a multitrait linear sire model. Heritability estimates for BCS and DC were moderate (0.25 and 0.22), and heritability estimates for mastitis and diseases other than mastitis were low (0.038 and 0.022). Between BCS and diseases other than mastitis, the genetic correlation was -0.22, whereas the genetic correlation was -0.16 between BCS and mastitis. The genetic correlation between DC and diseases other than mastitis was 0.43, and between DC and mastitis it was 0.27. The genetic correlation between BCS and DC was -0.61. Residual correlations were close to 0, except between BCS and DC (-0.37). Including DC as an indicator of diseases other than mastitis will increase the accuracy of the predicted breeding value for diseases, especially when the progeny group is small. Using BCS as an additional indicator of diseases did not increase the accuracy. Breeding for less DC will increase resistance to diseases.