Genomic inbreeding and selection signatures in the local dairy breed Icelandic Cattle.

Icelandic Cattle is the only dairy cattle breed native to Iceland. It currently numbers ca. 26 000 breeding females. We used 50k genotypes of over 8000 Icelandic Cattle to estimate genomic and pedigree-based inbreeding and to detect selection signatures using the integrated haplotype score. We used 47 Icelandic bulls genotyped with a 770k SNP chip to estimate LD decay for comparison with other Nordic dairy cattle breeds. We detected ROHs on the autosomes and computed ROH-based autosomal inbreeding coefficients. Average inbreeding coefficients according to pedigree and ROHs were 0.0621 and 0.101. Effective population sizes for the years 2009-2017 according to pedigree, ROHs, genomic relationship matrix, excess of homozygosity and individual increase in inbreeding were 81, 65, 60, 58 and 92 respectively. We identified three regions and six candidate genes that showed a signature of selection according to the integrated haplotype score (P < 0.05) on chromosomes 1, 16 and 23. The LD structure of Icelandic Cattle is shaped by a long period of isolation and a small founder population. The estimate of LD at distances closer than 0.3 Mb is lower in Icelandic Cattle than in Danish Jersey, but is higher than in Danish Holstein and Red Nordic Dairy Cattle breeds. Our findings show that inbreeding rates in Icelandic Cattle currently are sustainable according to FAO guidelines, and our results do not indicate severe historical inbreeding.

Genome-wide association analysis of milk yield traits in Nordic Red Cattle using imputed whole genome sequence variants.

BACKGROUND: The Nordic Red Cattle consisting of three different populations from Finland, Sweden and Denmark are under a joint breeding value estimation system. The long history of recording of production and health traits offers a great opportunity to study production traits and identify causal variants behind them. In this study, we used whole genome sequence level data from 4280 progeny tested Nordic Red Cattle bulls to scan the genome for loci affecting milk, fat and protein yields. RESULTS: Using a genome-wise significance threshold, regions on Bos taurus chromosomes 5, 14, 23, 25 and 26 were associated with fat yield. Regions on chromosomes 5, 14, 16, 19, 20 and 25 were associated with milk yield and chromosomes 5, 14 and 25 had regions associated with protein yield. Significantly associated variations were found in 227 genes for fat yield, 72 genes for milk yield and 30 genes for protein yield. Ingenuity Pathway Analysis was used to identify networks connecting these genes displaying significant hits. When compared to previously mapped genomic regions associated with fertility, significantly associated variations were found in 5 genes common for fat yield and fertility, thus linking these two traits via biological networks. CONCLUSION: This is the first time when whole genome sequence data is utilized to study genomic regions affecting milk production in the Nordic Red Cattle population. Sequence level data offers the possibility to study quantitative traits in detail but still cannot unambiguously reveal which of the associated variations is causative. Linkage disequilibrium creates difficulties to pinpoint the causative genes and variations. One solution to overcome these difficulties is the identification of the functional gene networks and pathways to reveal important interacting genes as candidates for the observed effects. This information on target genomic regions may be exploited to improve genomic prediction.

Fine mapping of a calving QTL on Bos taurus autosome 18 in Holstein cattle.

Decreased calving performance not only directly impacts the economic efficiency of dairy cattle farming but also influences public concern for animal welfare. Previous studies have revealed a QTL on Bos taurus autosome (BTA) 18 that has a large effect on calving traits in Holstein cattle. In this study, fine mapping of this QTL was performed using imputed high-density SNP chip (HD) genotypes followed by imputed next-generation sequencing (NGS) variants. BTA18 was scanned for seven direct calving traits in 6113 bulls with imputed HD genotypes. SNP rs136283363 (BTA18: 57 548 213) was consistently the most significantly associated SNP across all seven traits [e.g. p-value = 2.04 × 10(-59) for birth index (BI)]. To finely map the QTL region and to explore pleiotropic effects, we studied NGS variants within the targeted region (BTA18: 57 321 450-57 625 355) for associations with direct calving traits and with three conformation traits. Significant variants were prioritized, and their biological relevance to the traits was interpreted. Considering their functional relationships with direct calving traits, SIGLEC12, CD33 and CEACAM18 were proposed as candidate genes. In addition, pleiotropic effects of this QTL region on direct calving traits and conformation traits were observed. However, the extent of linkage disequilibrium combined with the lack of complete annotation and potential errors in the Bos taurus genome assembly hampered our efforts to pinpoint the causal mutation.

Quantitative trait loci markers derived from whole genome sequence data increases the reliability of genomic prediction.

This study investigated the effect on the reliability of genomic prediction when a small number of significant variants from single marker analysis based on whole genome sequence data were added to the regular 54k single nucleotide polymorphism (SNP) array data. The extra markers were selected with the aim of augmenting the custom low-density Illumina BovineLD SNP chip (San Diego, CA) used in the Nordic countries. The single-marker analysis was done breed-wise on all 16 index traits included in the breeding goals for Nordic Holstein, Danish Jersey, and Nordic Red cattle plus the total merit index itself. Depending on the trait's economic weight, 15, 10, or 5 quantitative trait loci (QTL) were selected per trait per breed and 3 to 5 markers were selected to tag each QTL. After removing duplicate markers (same marker selected for more than one trait or breed) and filtering for high pairwise linkage disequilibrium and assaying performance on the array, a total of 1,623 QTL markers were selected for inclusion on the custom chip. Genomic prediction analyses were performed for Nordic and French Holstein and Nordic Red animals using either a genomic BLUP or a Bayesian variable selection model. When using the genomic BLUP model including the QTL markers in the analysis, reliability was increased by up to 4 percentage points for production traits in Nordic Holstein animals, up to 3 percentage points for Nordic Reds, and up to 5 percentage points for French Holstein. Smaller gains of up to 1 percentage point was observed for mastitis, but only a 0.5 percentage point increase was seen for fertility. When using a Bayesian model accuracies were generally higher with only 54k data compared with the genomic BLUP approach, but increases in reliability were relatively smaller when QTL markers were included. Results from this study indicate that the reliability of genomic prediction can be increased by including markers significant in genome-wide association studies on whole genome sequence data alongside the 54k SNP set.

Short communication: Improving accuracy of Jersey genomic evaluations in the United States and Denmark by sharing reference population bulls.

The effect on prediction accuracy for Jersey genomic evaluations of Danish and US bulls from using a larger reference population was assessed. Each country contributed genotypes from 1,157 Jersey bulls to the reference population of the other. Data were separated into reference (US only, Danish only, and combined US-Danish) and validation (US only and Danish only) populations. Depending on trait (milk, fat, and protein yields and component percentages; productive life; somatic cell score; daughter pregnancy rate; 14 conformation traits; and net merit), the US reference population included 2,720 to 4,772 bulls and cows with traditional evaluations as of August 2009; the Danish reference population included 635 to 996 bulls. The US validation population included 442 to 712 bulls that gained a traditional evaluation between August 2009 and December 2013; the Danish validation population included 105 to 196 bulls with multitrait across-country evaluations on the US scale by December 2013. Genomic predicted transmitting abilities (GPTA) were calculated on the US scale using a selection index that combined direct genomic predictions with either traditional predicted transmitting ability for the reference population or traditional parent averages (PA) for the validation population and a traditional evaluation based only on genotyped animals. Reliability for GPTA was estimated from the reference population and August 2009 traditional PA and PA reliability. For prediction of December 2013 deregressed daughter deviations on the US scale, mean August 2009 GPTA reliability for Danish validation bulls was 0.10 higher when based on the combined US-Danish reference population than when the reference population included only Danish bulls; for US validation bulls, mean reliability increased by 0.02 when Danish bulls were added to the US reference population. Exchanging genotype data to increase the size of the reference population is an efficient approach to increasing the accuracy of genomic prediction when the reference population is small.

Identifying QTL and genetic correlations between fur quality traits in mink (Neovison vison).

Mapping of QTL affecting fur quality traits (guard hair length, guard hair thickness, density of wool, surface of the fur and quality) and skin length was performed in a three-generation mink population (F2 design). In the parental generation, Nordic Brown mink were crossed reciprocally with American Black short nap mink. In all, 1082 mink encompassing three generations were used for the analyses. The mink were genotyped for 104 microsatellites covering all 14 autosomes. The QTL analyses were performed by least-square regression implemented in gridqtl software. Genetic and phenotypic correlations and heritabilities were estimated using the average information-restricted maximum-likelihood method. Evidence was found for QTL affecting fur quality traits on nine autosomes. QTL were detected for guard hair thickness on chromosomes 1, 2, 3, 6 and 13; for guard hair length on chromosomes 2, 3 and 6; for wool density on chromosomes 6 and 13; for surface on chromosomes 7, 12 and 13; for quality on chromosomes 6, 7, 11 and 13; and for skin length on chromosomes 7 and 9. Proximity of locations of QTL for guard hair length, guard hair thickness and for wool density and quality suggests that some of the traits are in part under the influence of the same genes. Traits under the influence of QTL at close or identical positions also were traits that were strongly genotypically correlated. Based on the results of correlation analyses, the most important single traits influencing the quality were found to be density of wool, guard hair thickness and appearance of the surface.

Adding cows to the reference population makes a small dairy population competitive.

Small dairy breeds are challenged by low reliabilities of genomic prediction. Therefore, we evaluated the effect of including cows in the reference population for small dairy cattle populations with a limited number of sires in the reference population. Using detailed simulations, 2 types of scenarios for maintaining and updating the reference population over a period of 15yr were investigated: a turbo scheme exclusively using genotyped young bulls and a hybrid scheme with mixed use of genotyped young bulls and progeny-tested bulls. Two types of modifications were investigated: (1) number of progeny-tested bulls per year was tested at 6 levels: 15, 40, 60, 100, 250, and 500; and (2) each year, 2,000 first-lactation cows were randomly selected from the cow population for genotyping or, alternatively, an additional 2,000 first-lactation cows were randomly selected and typed in the first 2yr. The effects were evaluated in the 2 main breeding schemes. The breeding schemes were chosen to mimic options for the Danish Jersey cattle population. Evaluation criteria were annual monetary genetic gain, rate of inbreeding, reliability of genomic predictions, and variance of response. Inclusion of cows in the reference population increased monetary genetic gain and decreased the rate of inbreeding. The increase in genetic gain was larger for the turbo schemes with shorter generation intervals. The variance of response was generally higher in turbo schemes than in schemes using progeny-tested bulls. However, the risk was reduced by adding cows to the reference population. The annual genetic gain and the reliability of genomic predictions were slightly higher with more cows in the reference population. Inclusion of cows in the reference population is a rapid way to increase reliabilities of genomic predictions and hence increase genetic gain in a small population. An economic evaluation shows that genotyping of cows is a profitable investment.

Genomic selection strategies in a small dairy cattle population evaluated for genetic gain and profit.

The objective of this study was to evaluate a genomic breeding scheme in a small dairy cattle population that was intermediate in terms of using both young bulls (YB) and progeny-tested bulls (PB). This scheme was compared with a conventional progeny testing program without use of genomic information and, as the extreme case, a juvenile scheme with genomic information, where all bulls were used before progeny information was available. The population structure, cost, and breeding plan parameters were chosen to reflect the Danish Jersey cattle population, being representative for a small dairy cattle population. The population consisted of 68,000 registered cows. Annually, 1,500 bull dams were screened to produce the 500 genotyped bull calves from which 60 YB were selected to be progeny tested. Two unfavorably correlated traits were included in the breeding goal, a production trait (h(2)=0.30) and a functional trait (h(2)=0.04). An increase in reliability of 5 percentage points for each trait was used in the default genomic scenario. A deterministic approach was used to model the different breeding programs, where the primary evaluation criterion was annual monetary genetic gain (AMGG). Discounted profit was used as an indicator of the economic outcome. We investigated the effect of varying the following parameters: (1) increase in reliability due to genomic information, (2) number of genotyped bull calves, (3) proportion of bull dam sires that are young bulls, and (4) proportion of cow sires that are young bulls. The genomic breeding scheme was both genetically and economically superior to the conventional breeding scheme, even in a small dairy cattle population where genomic information causes a relatively low increase in reliability of breeding values. Assuming low reliabilities of genomic predictions, the optimal breeding scheme according to AMGG was characterized by mixed use of YB and PB as bull sires. Exclusive use of YB for production cows increased AMGG up to 3 percentage points. The results from this study supported our hypothesis that strong interaction effects exist. The strongest interaction effects were obtained between increased reliabilities of genomic estimated breeding values and more intensive use of YB. The juvenile scheme was genetically inferior when the increase in reliability was low (5 percentage points), but became genetically superior at higher reliabilities of genomic estimated breeding values. The juvenile scheme was always superior according to discounted profit because of the shorter generation interval and minimizing costs for housing and feeding waiting bulls.

Genomic predictions based on a joint reference population for the Nordic Red cattle breeds.

The main aim of this study was to compare accuracies of imputation and genomic predictions based on single and joint reference populations for Norwegian Red (NRF) and a composite breed (DFS) consisting of Danish Red, Finnish Ayrshire, and Swedish Red. The single nucleotide polymorphism (SNP) data for NRF consisted of 2 data sets: one including 25,000 markers (NRF25K) and the other including 50,000 markers (NRF50K). The NRF25K data set had 2,572 bulls, and the NRF50K data set had 1,128 bulls. Four hundred forty-two bulls were genotyped in both data sets (double-genotyped bulls). The DFS data set (DSF50K) included 50,000 markers of 13,472 individuals, of which around 4,700 were progeny-tested bulls. The NRF25K data set was imputed to 50,000 density using the software Beagle. The average error rate for the imputation of NRF25K decreased slightly from 0.023 to 0.021, and the correlation between observed and imputed genotypes changed from 0.935 to 0.936 when comparing the NRF50K reference and the NRF50K-DFS50K joint reference imputations. A genomic BLUP (GBLUP) model and a Bayesian 4-component mixture model were used to predict genomic breeding values for the NRF and DFS bulls based on the single and joint NRF and DFS reference populations. In the multiple population predictions, accuracies of genomic breeding values increased for the 3 production traits (milk, fat, and protein yields) for both NRF and DFS. Accuracies increased by 6 and 1.3 percentage points, on average, for the NRF and DFS bulls, respectively, using the GBLUP model, and by 9.3 and 1.3 percentage points, on average, using the Bayesian 4-component mixture model. However, accuracies for health or reproduction traits did not increase from the multiple population predictions. Among the 3 DFS populations, Swedish Red gained most in accuracies from the multiple population predictions, presumably because Swedish Red has a closer genetic relationship with NRF than Danish Red and Finnish Ayrshire. The Bayesian 4-component mixture model performed better than the GBLUP model for most production traits for both NRF and DFS, whereas no advantage was found for health or reproduction traits. In general, combining NRF and DFS reference populations was useful in genomic predictions for both the NRF and DFS bulls.

Genome-wide association study using high-density single nucleotide polymorphism arrays and whole-genome sequences for clinical mastitis traits in dairy cattle.

Mastitis is a mammary disease that frequently affects dairy cattle. Despite considerable research on the development of effective prevention and treatment strategies, mastitis continues to be a significant issue in bovine veterinary medicine. To identify major genes that affect mastitis in dairy cattle, 6 chromosomal regions on Bos taurus autosome (BTA) 6, 13, 16, 19, and 20 were selected from a genome scan for 9 mastitis phenotypes using imputed high-density single nucleotide polymorphism arrays. Association analyses using sequence-level variants for the 6 targeted regions were carried out to map causal variants using whole-genome sequence data from 3 breeds. The quantitative trait loci (QTL) discovery population comprised 4,992 progeny-tested Holstein bulls, and QTL were confirmed in 4,442 Nordic Red and 1,126 Jersey cattle. The targeted regions were imputed to the sequence level. The highest association signal for clinical mastitis was observed on BTA 6 at 88.97 Mb in Holstein cattle and was confirmed in Nordic Red cattle. The peak association region on BTA 6 contained 2 genes: vitamin D-binding protein precursor (GC) and neuropeptide FF receptor 2 (NPFFR2), which, based on known biological functions, are good candidates for affecting mastitis. However, strong linkage disequilibrium in this region prevented conclusive determination of the causal gene. A different QTL on BTA 6 located at 88.32 Mb in Holstein cattle affected mastitis. In addition, QTL on BTA 13 and 19 were confirmed to segregate in Nordic Red cattle and QTL on BTA 16 and 20 were confirmed in Jersey cattle. Although several candidate genes were identified in these targeted regions, it was not possible to identify a gene or polymorphism as the causal factor for any of these regions.

MBL1 genotypes in wild boar populations from Sweden, Austria, the Czech Republic, and Japan.

The single nucleotide polymorphism (SNP) G949T in the mannose-binding lectin ( MBL ) 1 gene has been associated with low MBL-A concentration in serum and detected at different frequencies in various European pig populations. However, the origin of this SNP is not known. Part of the MBL1 gene was sequenced in 12 wild boar/Large White crossbred pigs from the second backcross (BC 2 ) generation in a family material originating from two wild boar x Large White intercrosses. Also, MBL-A serum concentration was measured in the entire BC 2 generation (n = 45). Furthermore, the genotypes of 68 wild boars from Sweden, Austria, the Czech Republic, and Japan were determined in regard to five previously described SNPs in MBL1 . The T allele of G949T was present among the BC 2 animals. MBL-A serum concentration in the BC 2 animals showed a bimodal distribution, with one-third of the animals at levels between 0.7 and 1.6 µg mL(-1) and the remaining pigs at levels around 13 µg mL(-1) . There was a co-variation between the presence of the T allele and low MBL-A concentration in serum. The genotyping of the wild boars revealed differences between populations. The T allele of G949T was not detected in the Austrian and Japanese samples and is thus unlikely to be an original feature of wild boars. In contrast, it was present at high frequency (0.35) among the Swedish wild boars, probably representing a founder effect. Five MBL1 haplotypes were resolved. Only two of these were present among the Japanese wild boars compared to four in each of the European populations. This difference may reflect differences in selection pressure and population history.

Confirmation and fine-mapping of clinical mastitis and somatic cell score QTL in Nordic Holstein cattle.

A genome-wide association study of 2098 progeny-tested Nordic Holstein bulls genotyped for 36 387 SNPs on 29 autosomes was conducted to confirm and fine-map quantitative trait loci (QTL) for mastitis traits identified earlier using linkage analysis with sparse microsatellite markers in the same population. We used linear mixed model analysis where a polygenic genetic effect was fitted as a random effect and single SNPs were successively included as fixed effects in the model. We detected 143 SNP-by-trait significant associations (P < 0.0001) on 20 chromosomes affecting mastitis-related traits. Among them, 21 SNP-by-trait combinations exceeded the genome-wide significant threshold. For 12 chromosomes, both the present association study and the previous linkage study detected QTL, and of these, six were in the same chromosomal locations. Strong associations of SNPs with mastitis traits were observed on bovine autosomes 6, 13, 14 and 20. Possible candidate genes for these QTL were identified. Identification of SNPs in linkage disequilibrium with QTL will enable marker-based selection for mastitis resistance. The candidate genes identified should be further studied to detect candidate polymorphisms underlying these QTL.

Effect of marker-data editing on the accuracy of genomic prediction.

Genomic selection is a method to predict breeding values using genome-wide single-nucleotide polymorphism (SNP) markers. High-quality marker data are necessary for genomic selection. The aim of this study was to investigate the effect of marker-editing criteria on the accuracy of genomic predictions in the Nordic Holstein and Jersey populations. Data included 4429 Holstein and 1071 Jersey bulls. In total, 48,222 SNP for Holstein and 44,305 SNP for Jersey were polymorphic. The SNP data were edited based on (i) minor allele frequencies (MAF) with thresholds of no limit, 0.001, 0.01, 0.02, 0.05 and 0.10, (ii) deviations from Hardy-Weinberg proportions (HWP) with thresholds of no limit, chi-squared p-values of 0.001, 0.02, 0.05 and 0.10, and (iii) GenCall (GC) scores with thresholds of 0.15, 0.55, 0.60, 0.65 and 0.70. The marker data sets edited with different criteria were used for genomic prediction of protein yield, fertility and mastitis using a Bayesian variable selection and a GBLUP model. De-regressed EBV were used as response variables. The result showed little difference between prediction accuracies based on marker data sets edited with MAF and deviation from HWP. However, accuracy decreased with more stringent thresholds of GC score. According to the results of this study, it would be appropriate to edit data with restriction of MAF being between 0.01 and 0.02, a p-value of deviation from HWP being 0.05, and keeping all individual SNP genotypes having a GC score over 0.15.

Comparison of genomic predictions using medium-density (∼54,000) and high-density (∼777,000) single nucleotide polymorphism marker panels in Nordic Holstein and Red Dairy Cattle populations.

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.

Influence of model specifications on the reliabilities of genomic prediction in a Swedish-Finnish red breed cattle population.

Using a combined multi-breed reference population, this study explored the influence of model specification and the effect of including a polygenic effect on the reliability of genomic breeding values (DGV and GEBV). The combined reference population consisted of 2986 Swedish Red Breed (SRB) and Finnish Ayrshire (FAY) dairy cattle. Bayesian methodology (common prior and mixture models with different prior distribution settings for the marker effects) as well as a best linear unbiased prediction with a genomic relationship matrix [genomic best linear unbiased predictor (GBLUP)] was used in the prediction of DGV. Mixture models including a polygenic effect were used to predict GEBV. In total, five traits with low, high and medium heritability were analysed. For the models using a mixture prior distribution, reliabilities of DGV tended to decrease with an increasing proportion of markers with small effects. The influence of the inclusion of a polygenic effect on the reliability of DGV varied across traits and model specifications. Average correlation between DGV with the Mendelian sampling term, across traits, was highest (R(2) = 0.25) for the GBLUP model and decreased with increasing proportion of markers with large effects. Reliabilities increased when DGV and parent average information were combined in an index. The GBLUP model with the largest gain across traits in the reliability of the index achieved the highest DGV mean reliability. However, the polygenic models showed to be less biased and more consistent in the estimation of DGV regardless of the model specifications compared with the mixture models without the polygenic effect.

Genome-wide association study for calving traits in Danish and Swedish Holstein cattle.

A total of 22 quantitative trait loci (QTL) were detected on 19 chromosomes for direct and maternal calving traits in cattle using a genome-wide association study. Calving performance is affected by the genotypes of both the calf (direct effect) and dam (maternal effect). To identify the QTL contributing these effects to calving characteristics, we performed a genome-wide association study using a mixed-model analysis in Danish and Swedish Holstein cattle. The analysis incorporated 2,062 progeny-tested bulls, and 36,387 single nucleotide polymorphism markers on 29 bovine autosomes were analyzed for association with 14 calving traits. Strong evidence for the presence of QTL that affect calving traits was observed on chromosomes 4, 6, 12, 18, 20, and 25. The QTL intervals were generally smaller than those described in earlier linkage studies. The identification of calving trait-associated single nucleotide polymorphisms and mapping of the corresponding QTL in small chromosomal regions will facilitate the search for candidate calving performance genes and polymorphisms.

Reliabilities of genomic prediction using combined reference data of the Nordic Red dairy cattle populations.

This study investigated the possibility of increasing the reliability of direct genomic values (DGV) by combining reference populations. The data were from 3,735 bulls from Danish, Swedish, and Finnish Red dairy cattle populations. Single nucleotide polymorphism markers were fitted as random variables in a Bayesian model, using published estimated breeding values as response variables. In total, 17 index traits were analyzed. Reliabilities were estimated using a 5-fold cross validation, and calculated as the within-year squared correlation between estimated breeding values and DGV. Marker effects were estimated using reference populations from individual countries, as well as using a combined reference population from all 3 countries. Single-country reference populations gave mean reliabilities across 17 traits of 0.19 to 0.23, whereas the combined reference gave mean reliabilities of 0.26 for all populations. Using marker effects from 1 population to predict the other 2 gave a loss in mean reliability of 0.14 to 0.21 when predicting Swedish or Finnish animals with Danish marker effects, or vice versa. Using Swedish or Finnish marker effects to predict each other only showed a loss in mean reliability of 0.03 to 0.05. A combined Swedish-Finnish reference population led to an average reliability as high as that from the 3-country reference population, but somewhat different for individual traits. The results from this study show that it is possible to increase the reliability of DGV by combining reference populations from related populations.

Effect of imputing markers from a low-density chip on the reliability of genomic breeding values in Holstein populations.

The purpose of this study was to investigate the imputation error and loss of reliability of direct genomic values (DGV) or genomically enhanced breeding values (GEBV) when using genotypes imputed from a 3,000-marker single nucleotide polymorphism (SNP) panel to a 50,000-marker SNP panel. Data consisted of genotypes of 15,966 European Holstein bulls from the combined EuroGenomics reference population. Genotypes with the low-density chip were created by erasing markers from 50,000-marker data. The studies were performed in the Nordic countries (Denmark, Finland, and Sweden) using a BLUP model for prediction of DGV and in France using a genomic marker-assisted selection approach for prediction of GEBV. Imputation in both studies was done using a combination of the DAGPHASE 1.1 and Beagle 2.1.3 software. Traits considered were protein yield, fertility, somatic cell count, and udder depth. Imputation of missing markers and prediction of breeding values were performed using 2 different reference populations in each country: either a national reference population or a combined EuroGenomics reference population. Validation for accuracy of imputation and genomic prediction was done based on national test data. Mean imputation error rates when using national reference animals was 5.5 and 3.9% in the Nordic countries and France, respectively, whereas imputation based on the EuroGenomics reference data set gave mean error rates of 4.0 and 2.1%, respectively. Prediction of GEBV based on genotypes imputed with a national reference data set gave an absolute loss of 0.05 in mean reliability of GEBV in the French study, whereas a loss of 0.03 was obtained for reliability of DGV in the Nordic study. When genotypes were imputed using the EuroGenomics reference, a loss of 0.02 in mean reliability of GEBV was detected in the French study, and a loss of 0.06 was observed for the mean reliability of DGV in the Nordic study. Consequently, the reliability of DGV using the imputed SNP data was 0.38 based on national reference data, and 0.48 based on EuroGenomics reference data in the Nordic validation, and the reliability of GEBV using the imputed SNP data was 0.41 based on national reference data, and 0.44 based on EuroGenomics reference data in the French validation.

Consequences of elimination of the Rendement Napole allele from Danish Hampshire.

An elimination programme was carried out to remove the dominant Rendement Napole mutation (RN(-) ) from Danish Hampshire pigs. We reasoned that during and after the elimination of the RN(-) allele, genetic gain of production traits decreased while rate of inbreeding in the population increased compared to the period prior to elimination. The hypothesis was tested by estimating the genetic gain in seven production traits and measuring the rate of inbreeding in the population prior to and during the elimination period. Genetic gain was reduced for quantitative traits daily gain(30-100 kg) and feed conversion ratio, while gain for ultimate-pH, lean meat percentage and slaughter loss were increased slightly. There were no changes in genetic gain for daily gain(birth-30 kg) and conformation. RN polymorphism affected several of the quantitative traits. The RN(-) mutation had a dominant effect on the traits daily gain(birth-30 kg) , daily gain(30-100 kg) , slaughter loss, lean meat percentage and ultimate-pH. It exhibited overdominance for feed conversion ratio and additive effect for conformation. Rate of inbreeding decreased during the elimination of RN(-) . Our findings indicate that the consequences of the elimination programme were not as serious as were feared and that a carefully designed preselection strategy may avoid unacceptable loss of genetic gain and excessive loss of genetic variation.

Genome-wide association mapping for female fertility traits in Danish and Swedish Holstein cattle.

A genome-wide association study was conducted using a mixed model analysis for QTL for fertility traits in Danish and Swedish Holstein cattle. The analysis incorporated 2,531 progeny tested bulls, and a total of 36,387 SNP markers on 29 bovine autosomes were used. Eleven fertility traits were analyzed for SNP association. Furthermore, mixed model analysis was used for association analyses where a polygenic effect was fitted as a random effect, and genotypes at single SNPs were successively included as a fixed effect in the model. The Bonferroni correction for multiple testing was applied to adjust the significance threshold. Seventy-four SNP-trait combinations showed chromosome-wide significance, and five of these were significant genome-wide. Twenty-four QTL regions on 14 chromosomes were detected. Strong evidence for the presence of QTL that affect fertility traits were observed on chromosomes 3, 5, 10, 13, 19, 20, and 24. The QTL intervals were generally smaller than those described in earlier linkage studies. The identification of fertility trait-associated SNPs and mapping of the corresponding QTL in small chromosomal regions reported here will facilitate searches for candidate genes and candidate polymorphisms.