Genome-wide association study for milking speed in Fleckvieh cattle.

Milking speed is an important trait influencing udder health of dairy cows as well as labor efficiency. Yet, it has received little attention in genomic association studies. The main objective of this study was to determine regions and genes on the genome with a potential effect on milking speed in Fleckvieh (dual purpose Simmental) cattle. Genome-wide association studies were conducted using de-regressed breeding values of bulls as phenotypes. Six SNP on 4 autosomes were significantly associated with milking speed for additive effects. Significant regions on BTA4 and BTA19 correspond with findings for other dairy cattle breeds. Based on the observation of Fleckvieh breed managers, variation of milking speed in batches of daughters of some bulls is much higher than in daughter groups of other bulls. This difference in within family variation may be caused by transmission of alternative alleles of bulls being heterozygous for a gene affecting milking speed. To check on this, we considered standard deviation of yield deviations in milking speed of half-sib daughters as a new trait and performed GWAS for dominance effects. One signal on BTA5 passed the genome wide Bonferroni threshold that corresponded to the significant signal from standard GWAS on de-regressed breeding values. The key conclusion of this study is that several strong genomic signals were found for milking speed in Fleckvieh cattle and that the strongest of them are supported by similar findings in Brown Swiss and Holstein Friesian cattle. Milking speed is a complex trait whose sub-processes have not yet been elucidated in detail. Hence, it remains a challenge to link the associated regions on the genome with causal genes and their functions.

How pedigree errors affect genetic evaluations and validation statistics.

Pedigrees used in genetic evaluations contain errors. Because of such errors, assumptions regarding the relatedness among individuals in genetic evaluation models are wrong. Consequences of that have been investigated in earlier studies focusing on models that did not account for genomic information yet. The objective of this work was to investigate the effects of pedigree errors on the results from genetic evaluations using the single-step model, and the effect of such effects on results from validation studies with forward prediction. We used a real pedigree (n = 361,980) and real genotypes (n = 25,950) of Fleckvieh cattle, sampled in a way to provide a good consistency between pedigree and genomic relationships. Given the real pedigree and genotypes, true breeding values (TBV) were simulated to have a covariance structure equal to the matrix H assumed in a single-step model. Based on TBV, phenotypes were simulated with a heritability of 0.25. Genetic evaluations were conducted with a conventional animal model (i.e., without genomic information) and a single-step animal model under scenarios using either the correct pedigree or a pedigree containing 5%, 10%, or 20% of wrong records. Wrong records were simulated by randomly assigning wrong sires to nongenotyped females. The increasing rates of pedigree errors led to decreasing correlations between TBV and EBV and lower standard deviations of predictions. Less variation was observed because pedigree errors operate actually as a random exchange of daughters among bulls, making them look more similar to each other than they actually are. This occurs of course only when animals have progeny. Therefore, this decreased variation was more pronounced for progeny tested bulls than for young selection candidates. In a forward prediction validation scenario, the stronger decrease in variation when animals get progeny caused an apparent inflation of early predictions. This phenomenon may contribute to the usually observed problem of inflation of early predictions observed in validation studies.

Genotype-by-environment interactions at the trait level and total merit index level for milk production and functional traits in Brown Swiss cattle.

The production environments of the German-Austrian Brown Swiss population show a wide range due to differences in topography, landscapes, local climates, and different farm management systems. Extensive production systems such as organic farming have become increasingly popular in recent decades because of interest in sustainability and consumer preferences. Compared with conventional farmers, organic farmers put more weight on fitness traits. Besides the official total merit index (TMI), a selection index applying relative economic weights (REWs) suitable for organic production systems is provided for Brown Swiss cattle in Germany. The aim of the study was to investigate genotype-by-environment interactions (GxE) for milk production traits and functional traits (including longevity, fertility traits, and calving traits) in a sample of the German-Austrian Brown Swiss population housed in Baden-Wuerttemberg (southern Germany) by applying bivariate and random regression sire models. For bivariate analyses, the production environment was binary classified by farm management system (organic and conventional) and altitude of farm location (above or below 800 m above sea level (ASL)). Milk energy yields (MEY) obtained from herd effects were used as continuously scaled environmental descriptor in the reaction norm approach. The TMIs for sires were calculated based on breeding values estimated with different models and environment-specific REWs to determine possible GxE at TMI levels and rerankings of sires. In bivariate analyses, genetic correlations at the trait level were high and ranged from rg = 0.99 (calving to first insemination, cystic ovaries, and maternal stillbirth rate) to rg = 0.79 (first insemination to conception for altitude). Except for the latter, no severe GxE were found at the trait level using the bivariate models. Fat yield was the only trait showing minor GxE in the reaction norm model approach. Investigating the environmental sensitivity at the TMI level revealed rank correlations between the different environment-specific TMIs that were close to unity, implying no severe reranking effects. The results show no need to account for different environments in Brown Swiss cattle breeding programs.

Short communication: Calculating analytical reliabilities for single-step predictions.

It has been shown that single-step genomic BLUP (ssGBLUP) can be reformulated, resulting in an equivalent SNP model that includes the explicit imputation of gene contents of all ungenotyped animals in the pedigree. This reformulation reveals the underlying mechanism enabling ungenotyped animals to contribute information to genotyped animals via estimates of marker effects and consequently to the reliability of genomic predictions, a key feature generally associated with the single-step approach. Irrespective of which BLUP formulation is used for genomic prediction, with increasing numbers of genotyped animals, the marker-oriented model is recommended when calculating the reliabilities of genomic predictions. This approach has the advantage of a manageable and stable size of the model matrix that needs to be inverted to calculate analytical prediction error variances of marker effects, an advantage that also holds for prediction with the single-step model. However, when including imputed genotypes in the design matrix of marker effects, an additional imputation residual term has to be considered to account for the prediction error of imputation. We summarize some of the theoretical aspects associated with the calculation of analytical reliabilities of single-step predictions. Derivations are based on the equivalent reformulation of ssGBLUP as a marker-oriented model and the calculation of prediction error variances of marker effects. We propose 2 approximations that allow for a substantial reduction of the complexity of the matrix operations involved, while retaining most of the relevant information required for reliability calculations. We additionally provide a general framework for an implementation of single-step reliability approximation using standard animal model reliabilities as a starting point. Finally, we demonstrate the effectiveness of the proposed approach using a small example extracted from data of the routine evaluation on dual-purpose Fleckvieh (Simmental) cattle.

Technical note: Methods for interim prediction of single-step breeding values for young animals.

Single-step genomic evaluations have the advantage of simultaneously combining all pedigree, phenotypic, and genotypic information available. However, systems with a large number of genotyped animals have some computational challenges. In many genomic breeding programs, genomic predictions of young animals should become available for selection decisions in the shortest time possible, which requires either a very effective estimation or an approximation with negligible loss in accuracy. We investigated different procedures for predicting breeding values of young genotyped animals without setting up the full single-step system augmented for the additional genotypes. Methods were based on transmitting the information from single-step breeding values of genotyped animals that took part in the previous full run to young animals, either through genomic relationships or through a marker-based model. The different procedures were tested on real data from the April 2017 run of the German-Austrian official genomic evaluation for Fleckvieh. The data set included 62,559 genotyped animals and was used to run single-step evaluations for 23 conformation traits. A further data set comprising 1,768 young animals was used for interim prediction and we called it the validation set. The reference values for validation were the predicted breeding values of the young animals from a full single-step run containing the genotypes of all 64,327 animals. Correlations between the approximated predictions and those from the full single-step run also containing genotypes from young animals averaged 0.9932 for the best method (from 0.990 to 0.995 across traits). In conclusion, prediction of single-step breeding values for young animals can be well approximated using systems of size equal to the number of markers.

Considering dominance in reduced single-step genomic evaluations.

Single-step models including dominance can be an enormous computational task and can even be prohibitive for practical application. In this study, we try to answer the question whether a reduced single-step model is able to estimate breeding values of bulls and breeding values, dominance deviations and total genetic values of cows with acceptable quality. Genetic values and phenotypes were simulated (500 repetitions) for a small Fleckvieh pedigree consisting of 371 bulls (180 thereof genotyped) and 553 cows (40 thereof genotyped). This pedigree was virtually extended for 2,407 non-genotyped daughters. Genetic values were estimated with the single-step model and with different reduced single-step models. Including more relatives of genotyped cows in the reduced single-step model resulted in a better agreement of results with the single-step model. Accuracies of genetic values were largest with single-step and smallest with reduced single-step when only the cows genotyped were modelled. The results indicate that a reduced single-step model is suitable to estimate breeding values of bulls and breeding values, dominance deviations and total genetic values of cows with acceptable quality.

Short communication: The role of genotypes from animals without phenotypes in single-step genomic evaluations.

In a 2-step genomic system, genotypes of animals without phenotypes do not influence genomic prediction of other animals, but that might not be the case in single-step systems. We investigated the effects of including genotypes from culled bulls on the reliability of genomic predictions from single-step evaluations. Four scenarios with a constant amount of phenotypic information and increasing numbers of genotypes from culled bulls were simulated and compared with respect to prediction reliability. With increasing numbers of genotyped culled bulls, there was a corresponding increase in prediction reliability. For instance, in our simulation scenario the reliability for selection candidates was twice as large when all culled bulls from the last 4 generations were included in the analysis. Single-step evaluations imply the imputation of all nongenotyped animals in the pedigree. We showed that this imputation was increasingly more accurate as increasingly more genotypic information from the culled bulls was taken into account. This resulted in higher prediction reliabilities. The extent of the benefit from including genotypes from culled bulls might be more relevant for small populations with low levels of reliabilities.

Short communication: The effect of genotyping cows to improve the reliability of genomic predictions for selection candidates.

In this study we investigate the potential of enlarging the reference population for genomic prediction in dairy cattle by routinely genotyping a random sample of the first-crop daughters of every AI bull in the breeding program. We analyzed small nuclear pedigrees, each consisting of a genotyped selection candidate and 3 generations of genotyped male ancestors. Genotypes were taken from the genomic routine evaluation of Fleckvieh cattle in Germany and Austria. The phenotypic information of a daughter of any one male in each of these pedigrees was either considered to be part of the daughter yield deviation of the corresponding sire, or was assumed to be an individually observed genotyped daughter of this sire. Daughter genotypes in this case were simulated from phased haplotypes of their sires and random maternal gametes drawn from a haplotype library. We measured the gain from genotyping daughters as the increase in model-based theoretical reliability of the genomic prediction for a putative selection candidate. We expressed the improvements as a marginal increase, corresponding to an increase in reliability at a reliability baseline level of zero, to simplify comparisons. Results were encouraging with 2 to 40% of marginal reliability increase for selection candidates depending on the assumed heritability of the trait and the number of daughters modeled to be genotyped in the design.

How imputation errors bias genomic predictions.

The objective of this study was to investigate in detail the biasing effects of imputation errors on genomic predictions. Direct genomic values (DGV) of 3,494 Brown Swiss selection candidates for 37 production and conformation traits were predicted using either their observed 50K genotypes or their 50K genotypes imputed from a mimicked 6K chip. Changes in DGV caused by imputation errors were shown to be systematic. The DGV of top animals were, on average, underestimated and that of bottom animals were, on average, overestimated when imputed genotypes were used instead of observed genotypes. This pattern might be explained by the fact that imputation algorithms will usually suggest the most frequent haplotype from the sample whenever a haplotype cannot be determined unambiguously. That was empirically shown to cause an advantage for the bottom animals and a disadvantage for the top animals.

Short communication: Importance of introgression for milk traits in the German Vorderwald and Hinterwald cattle.

The subject of the present study was to analyze the influence of genetic introgression on milk yield performance of the German local Vorderwald and Hinterwald cattle breeds. Deviations of milk yield, fat yield, and protein yield of cows as well as pedigree information were analyzed. A sire model was used to estimate genetic trend and effects of the migrant breeds. Migrant contributions to Vorderwald cattle were high and have been rising even in the recent past. The effects of these breeds on milk yield performance were positive. Montbéliarde cattle not only had the largest effect on milk production of Vorderwald cattle but also the highest genetic contribution to this breed. Genetic introgression with Montbéliarde continued until recently. This suggests that introgression of high-yielding breeds is still a preferred method for genetic improvement of local breeds, even though it diminishes their value for conservation. Hence, the current population management has too little focus on the preservation of genetic uniqueness. In comparison, migrant breed contributions to the Hinterwald cattle, a breed with a unique phenotype and an own niche, were moderate and almost constant over the time. For the Hinterwald cattle, no significant effect of migrant breeds could be detected, which suggests that population management has different priorities in different endangered breeds. We conclude that not only the registration of animals from local breeds but also the breeding programs themselves should be supported and need to be controlled.

On the limited increase in validation reliability using high-density genotypes in genomic best linear unbiased prediction: observations from Fleckvieh cattle.

This study investigated reliability of genomic predictions using medium-density (40,089; 50K) or high-density (HD; 388,951) marker sets. We developed an approximate method to test differences in validation reliability for significance. Model-based reliability and the effect of HD genotypes on inflation of predictions were analyzed additionally. Genomic breeding values were predicted for at least 1,321 validation bulls based on phenotypes and genotypes of at least 5,324 calibration bulls by means of a linear model in milk, fat, and protein yield; somatic cell score; milkability; muscling; udder, feet, and legs score as well as stature. In total, 1,485 bulls were actually HD genotyped and HD genotypes of the other animals were imputed from 50K genotypes using FImpute software. Validation reliability was measured as the coefficient of determination of the weighted regression of daughter yield deviations on predicted breeding values divided by the reliability of daughter yield deviations and inflation was evaluated by the slope of this regression. Model-based reliability was calculated from the model. Distributions for validation reliability of 50K markers were derived by repeated sampling of 50,000-marker samples from HD to test differences in validation reliability statistically. Additionally, the benefit of HD genotypes in validation reliability was tested by repeated sampling of validation groups and calculation of the difference in validation reliability between HD and 50K genotypes for the sampled groups of bulls. The mean benefit in validation reliability of HD genotypes was 0.015 compared with real 50K genotypes and 0.028 compared with 50K samples from HD affected by imputation error and was significant for all traits. The model-based reliability was, on average, 0.036 lower and the regression coefficient was 0.036 closer to the expected value with HD genotypes. The observed gain in validation reliability with HD genotypes was similar to expectations based on the number of markers and the effective number of segregating chromosome segments. Sampling error in the marker-based relationship coefficients causing overestimation of the model-based reliability was smaller with HD genotypes. Inflation of the genomic predictions was reduced with HD genotypes, accordingly. Similar effects on model-based reliability and inflation, but not on the validation reliability, were obtained by shrinkage estimation of the realized relationship matrix from 50K genotypes.

Revisiting the quantitative trait loci for milk production traits on BTA6.

A parallel association study was performed in two independent cattle populations based on 41 validated, targeted single nucleotide polymorphisms (SNPs) and four microsatellite markers to re-evaluate the multiple quantitative trait loci (QTL) architecture for milk performance on bovine chromosome 6 (BTA6). Two distinct QTL located in the vicinity of the middle region of BTA6, but differing unambiguously regarding their effects on milk composition and yield traits were validated in the German Holstein population. A highly significant association of the protein variant ABCG2 p.Tyr581Ser with milk composition traits reconfirmed the causative molecular relevance of the ABCG2 gene in QTL region 1, whereas in QTL region 2, significant and tentative associations between gene variants RW070 and RW023 (located in the promoter region and exon 9 of the PPARGC1A gene for milk yield traits) were detected. For the German Fleckvieh population, only RW023 showed a tentative association with milk yield traits, whereas those loci with significant effects in German Holsteins (ABCG2 p.Tyr581Ser, RW070) showed fixed alleles. Even though our new data highlight two variants in the PPARGC1A gene (RW023, RW070) in QTL region 2, based on the results of our study, currently no unequivocal conclusion about the causal background of this QTL affecting milk yield traits can be drawn. Notably, the German Holstein and Fleckvieh populations, known for their divergent degree of dairy type, differ substantially in their allele frequencies for the growth-associated NCAPG p.Ile442Met locus.

Short communication: Genomic selection using a multi-breed, across-country reference population.

Three breeds (Fleckvieh, Holstein, and Jersey) were included in a reference population, separately and together, to assess the accuracy of prediction of genomic breeding values in single-breed validation populations. The accuracy of genomic selection was defined as the correlation between estimated breeding values, calculated using phenotypic data, and genomic breeding values. The Holstein and Jersey populations were from Australia, whereas the Fleckvieh population (dual-purpose Simmental) was from Austria and Germany. Both a BLUP with a multi-breed genomic relationship matrix (GBLUP) and a Bayesian method (BayesA) were used to derive the prediction equations. The hypothesis tested was that having a multi-breed reference population increased the accuracy of genomic selection. Minimal advantage existed of either GBLUP or BayesA multi-breed genomic evaluations over single-breed evaluations. However, when the goal was to predict genomic breeding values for a breed with no individuals in the reference population, using 2 other breeds in the reference was generally better than only 1 breed.

Confirmation and refinement of a QTL on BTA5 affecting milk production traits in the Fleckvieh dual purpose cattle breed.

We analysed a QTL affecting milk yield (MY), milk protein yield (PY) and milk fat yield (FY) in the dual purpose cattle breed Fleckvieh on BTA5. Twenty-six microsatellite markers covering 135 cM were selected to analyse nine half-sib families containing 605 sons in a granddaughter design. We thereby assigned two new markers to the public linkage map using the CRI-MAP program. Phenotypic records were daughter yield deviations (DYD) originating from the routinely performed genetic evaluations of breeding animals. To determine the position of the QTL, three different approaches were applied: interval mapping (IM), linkage analysis by variance component analysis (LAVC), and combined linkage disequilibrium (LD) and linkage (LDL) analysis. All three methods mapped the QTL in the same marker interval (BM2830-ETH152) with the greatest test-statistic value at 118, 119.33 and 119.33 cM respectively. The positive QTL allele simultaneously increases DYD in the first lactation by 272 kg milk, 7.1 kg milk protein and 7.0 kg milk fat. Although the mapping accuracy and the significance of a QTL effect increased from IM over LAVC to LDL, the confidence interval was large (13, 20 and 24 cM for FY, MY and PY respectively) for the positional cloning of the causal gene. The estimated averages of pair wise marker LD with a distance <5 cM were low (0.107) and reflect the large effective population size of the Fleckvieh subpopulation analysed. This low level of LD suggests a need for increase in marker density in following fine mapping steps.

Genetic parameters for milkability from the first three lactations in Fleckvieh cows.

Test-day records for average flow rate (AFR) from the routine dairy recording from Bavarian Fleckvieh cows were analysed. Two data sets with observations on approximately 20 000 cows each were sampled from the total data set. For the estimation of variance parameters, a two-step approach was applied. In a first step multiple-trait restricted maximum likelihood (REML) analyses were carried out. For each of the first three lactations, six time periods with up to 33 days were defined. An algorithm for iterative summing of expanded part matrices was applied in order to combine the estimates. In a second step covariance functions (CF) for additive-genetic variances and non-genetic animal variances were derived using second-order Legendre polynomials plus an exponential term. Estimates of test-day heritability for AFR ranged from 0.21 to 0.40, and were largest in lactation 1. For lactations 1 and 3, heritabilities decreased considerably towards the end of lactation. Genetic correlation estimates within lactation decreased as the distance between days in milk (DIM) increased. Genetic correlations between corresponding DIM in the three lactations were generally large, ranging from 0.80 to 0.99. The largest estimates were found between DIM from lactations 2 and 3. Results from this study suggest that including AFR data from second and third lactations in genetic evaluation systems could the improve accuracy of genetic selection.

Strategies for estimating genetic parameters in marker-assisted best linear unbiased predictor models in dairy cattle.

An appropriate strategy to estimate variance components and breeding values in genetic models with quantitative trait loci (QTL) was developed for a dairy cattle breeding scheme by utilizing simulated data. Reliable estimates for variance components in QTL models are a prerequisite in fine-mapping experiments and for marker-assisted genetic evaluations. In cattle populations, only a small fraction of the population is genotyped at genetic markers, and only these animals are included in marker-assisted genetic evaluation models. Phenotypic information in these models are precorrected phenotypes [daughter yield deviations (DYD) for bulls, yield deviations (YD) for cows] estimated by standard animal models from the entire population. Because DYD and YD may represent different amounts of information, the problem of weighting these 2 types of information appropriately arises. To detect the best combination of phenotypes and weighting factors, a stochastic simulation for a trait representing milk yield was used. The results show that DYD models are generally optimal for estimating QTL variance components, but properties of estimates depend strongly on weighting factors. An example for the benefit in selection of using YD is shown for the selection among paternal half-sibs inheriting alternative QTL alleles. Even if QTL effects are small, marker-assisted best unbiased linear prediction can improve the selection among half-sibs, because the Mendelian sampling variance within family can be exploited, especially in DYD-YD models. Marker-assisted genetic evaluation models should also include YD for cows to ensure that marker-assisted selection improves selection even for moderate QTL effects (> or =10%). A useful strategy for practical implementation is to estimate variance components in DYD models and breeding values in DYD-YD models.

Multiplicative random regression model for heterogeneous variance adjustment in genetic evaluation for milk yield in Simmental.

A multiplicative random regression (M-RRM) test-day (TD) model was used to analyse daily milk yields from all available parities of German and Austrian Simmental dairy cattle. The method to account for heterogeneous variance (HV) was based on the multiplicative mixed model approach of Meuwissen. The variance model for the heterogeneity parameters included a fixed region x year x month x parity effect and a random herd x test-month effect with a within-herd first-order autocorrelation between test-months. Acceleration of variance model solutions after each multiplicative model cycle enabled fast convergence of adjustment factors and reduced total computing time significantly. Maximum Likelihood estimation of within-strata residual variances was enhanced by inclusion of approximated information on loss in degrees of freedom due to estimation of location parameters. This improved heterogeneity estimates for very small herds. The multiplicative model was compared with a model that assumed homogeneous variance. Re-estimated genetic variances, based on Mendelian sampling deviations, were homogeneous for the M-RRM TD model but heterogeneous for the homogeneous random regression TD model. Accounting for HV had large effect on cow ranking but moderate effect on bull ranking.

Analysis of genotype by environment interaction for milk yield traits in first lactation of Simmental cattle.

The breeding goal for Simmental cattle is derived for intensively managed dairy farms. Its suitability for extensive farms was addressed by analysing possible genotype by environment interaction (G x E) between the management levels for first lactation milk yield traits. A first analysis was performed with the data collected from 300 000 purebred daughters of 278 second crop bulls born in Bavaria in 1993 and 1994. The farms were classified by herd-year-effect, using the sum of fat and protein yields into two levels of management, either with 33 or 10% quantiles, corresponding to approximately 100 000 cows and 30 000 cows, respectively. The comparison was based on 'daughter yield' deviations (DYD). Correlations between DYD of extensive and intensive environments were 0.90, 0.91 and 0.87 for milk, fat and protein yield (kg) for 33% quantiles, respectively. Corresponding correlations for 10% quantiles were 0.85, 0.83 and 0.77. Despite high correlations, 50 out of 149 sires showed significant differences between DYD in different environments. Bulls with higher DYD for milk yield on intensive farms were superior in all environments. For the second analysis extensive and intensive farms in northern and southern Bavaria were chosen at random. Approximately 20 000 cows in each management class were used for the estimation of genetic parameters. In both regions phenotypic and additive-genetic variances were higher in the intensively managed herds. Likewise heritabilities were higher for fat and protein yield, but not for milk where higher heritabilities were observed in 33% quantiles. Genetic correlations between extensive and intensive environments were 0.97 and above (33% quantiles). Ten per cent quantiles led to lower genetic correlations (0.90-0.95). Although no serious re-ranking effects of sires were evident, the scale effect and the differences in genetic parameters should be taken into consideration in practical breeding.