Prediction of genetic correlations and international breeding values for missing traits

Prediction of genetic merit for missing traits is possible by combining available indicator traits. Indicator traits were combined using genetic correlations obtained from multiple regression equations of estimated genetic correlations among available indicator traits on variables explaining production circumstances and trait definitions. This prediction of missing traits was closer to actual breeding values than breeding values for any of the indicator traits. This was verified by evaluating clinical mastitis in each of the Nordic countries as a missing trait. The derived methodology was used to predict breeding values for clinical mastitis in the United States for local and international bulls with an average reliability of 43%.     

Mating allocations in Holstein combining genomic information and linear programming optimization at the herd level

In this study, we explored mating allocation in Holstein using genomic information for 24,333 Holstein females born in Denmark, Finland, and Sweden. We used 2 data sets of bulls: the top 50 genotyped bulls and the top 25 polled genotyped bulls on the Nordic total merit scale. We used linear programming to optimize economic scores within each herd, considering genetic level, genetic relationship, semen cost, the economic impact of genetic defects, polledness, and β-casein. We found that it was possible to reduce genetic relationships and eliminate expression of genetic defects with minimal effect on the genetic level in total merit index. Compared with maximizing only Nordic total merit index, the relative frequency of polled offspring increased from 13.5 to 22.5%, and that of offspring homozygous for β-casein (A2A2) from 66.7 to 75.0% in one generation, without any substantial negative impact on other comparison criteria. Using only semen from polled bulls, which might become necessary if dehorning is banned, considerably reduced the genetic level. We also found that animals carrying the polled allele were less likely to be homozygous for β-casein (A2A2) and more likely to be carriers of the genetic defect HH1. Hence, adding economic value to a monogenic trait in the economic score used for mating allocation sometimes negatively affected another monogenetic trait. We recommend that the comparison criteria used in this study be monitored in a modern genomic mating program.     

Derivation of sustainable breeding goals for dairy cattle using selection index theory

The objective was to present 2 methods for the derivation of nonmarket values for functional traits in dairy cattle using deterministic simulation and selection index theory. A nonmarket value can be a value representing animal welfare and societal influences for animal production, which can be added to market economic values in the breeding goal to define sustainable breeding goals. The first method was restricted indices. A consequence of adding a nonmarket value to a market economic value for a given functional trait is less selection emphasis on milk yield. In the second method, the loss in selection response in milk resulting from greater emphasis on functional traits was quantified. The 2 methods were demonstrated using a breeding goal for dairy cattle with 4 traits (milk yield, mastitis resistance, conception rate, and stillbirth). Nonmarket values derived separately using restricted indices were 0.4 and 2.6 times the value of market economic values for mastitis resistance and conception rate, respectively. Nonmarket values for mastitis resistance and conception rate were both lower when derived simultaneously than when derived separately. This was due to the positive genetic correlation between mastitis resistance and conception rate, and because both traits are negatively correlated with milk yield. Using the second method and accepting a 5% loss in selection response for milk yield, nonmarket values for mastitis, conception rate, and stillbirth were 0.3, 1.4, and 2.9 times the market economic values. It was concluded that the 2 methods could be used to derive nonmarket values for functional traits in dairy cattle.     

Simulating consequences of choosing a breeding goal for organic dairy production

In Denmark, Finland, and Sweden, the Nordic Total Merit index is used as the breeding selection tool for both organic and conventional dairy farmers based on common economic models for conventional dairy farming. Organic farming is based on the principles of organic agriculture (POA) defined by the International Federation of Organic Agriculture Movements. These principles are not set up with an economic point of view, and therefore it may be questionable to use a breeding goal (BG) for organic dairy production based on economic models. In addition to economics and the principles of organic agriculture, it is important to look at farmers' preferences for improving BG traits when setting up a BG for organic farming. The aim of this research was to set up, simulate, and compare long-term effects of different BG for organic and conventional dairy production systems based on economic models, farmers' preferences, and POA, with particular emphasis on disease resistance or on roughage consumption and feed efficiency. The BG based on economic models and on farmers' preferences were taken from previous studies. The other BG were desired gains indices, set up by means of a questionnaire about relatedness between the POA and BG traits. Each BG was simulated in the stochastic simulation program ADAM. The BG based on POA, with particular emphasis on disease resistance or on roughage consumption and feed efficiency, caused favorable genetic gain in all 12 traits included in this study compared with 6 traits for the other BG. The BG based on POA, with particular emphasis on disease resistance or on roughage consumption and feed efficiency, were very different from BG for organic and conventional production based on economic models and farmers' preferences in both simulated genetic change and correlations between BG. The BG that was created based on the principles of organic agriculture could be used as a specific index for organic dairy farming in Denmark, but this index was economically not very sustainable. Hence, an intermediate breeding goal could be developed by breeding companies to address both economics and the principles of organic agriculture.     

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.     

Selection for profit in cattle: I. Economic weights for purebred dairy cattle in the Czech Republic

A bioeconomic model for dairy cattle production was used to estimate economic values of 18 traits for dairy sires in purebred Holstein and Czech Fleckvieh populations. Economic values were defined as partial derivatives of the profit function with respect to each trait in a closed production system with dairy cow herds and integrated fattening of bulls. All revenues and costs associated with cows calving in the herds within one year and with their progeny were discounted at 5% per annum back to the date of calving. Calculations were carried out for the situation in the Czech Republic in 2005 (scenario 1: market quotas for milk yield and fat percentage) and for the expected situation in 2015 (scenario 2: free market). The relative economic importance of each trait was expressed as a ratio of the standardized economic value of that trait (its marginal economic value multiplied by its genetic standard deviation) to the standardized economic value of 305-d milk yield, with average fat and protein percentages. In addition to milk yield, somatic cell score was the second most important trait, achieving 32% to 43% of the value for milk yield in both scenarios. The relative importance of milk components differed notably between scenarios. The relative importance was approximately zero for protein and from −14 to −23% for fat percentage in scenario 1, but changed to 38% for protein and 27 to 31% for fat percentage in scenario 2. In both scenarios and for both breeds, the relative economic values for somatic cell score and length of productive life of cows were similar to those for fat and protein percentages in scenario 2. The smallest relative economic values (less than 4% of the relative importance of milk yield) were for birth weight, conception rate of heifers, and carcass traits. In conclusion, relative emphasis among traits in the breeding objective for Czech dairy cattle should be reassessed according to the expected situation after shifting to a free market economy in 2015.     

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² = 0.30) and a functional trait (h² = 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.     

Development of a breeding program for improving the milk yield performance of Butana cattle under smallholder production conditions using a stochastic simulation approach

Butana is one of the local dairy cattle breeds of Sudan commonly kept by smallholder producers. This breed has been strongly promoted to advance the dairy production sector in the country. The main problem, however, is the lack of a systematic breeding program that involves smallholder producers. The aim of the current study was to identify the most promising design for a breeding program to improve the milk yield performance of Butana cattle under smallholder production conditions. In total, 3 breeding scenarios, including (1) the use of farm bulls, (2) the use of village bulls, and (3) the rotational use of village bulls within village groups, were simulated using a stochastic simulation approach. For each breeding scenario, 3 selection methods for bulls were considered, namely random mating, phenotypic selection, and selection based on estimated breeding value (EBV). The results showed that no genetic gain was realized with random mating in all breeding scenarios. In the farm bull breeding scenario, annual genetic gain (standard deviation units) ranged from 0.01 to 0.19 (phenotypic selection) and from 0.01 to 0.39 (selection based on EBV). In the village bull breeding scenarios, the annual genetic gain ranged from 0.01 to 0.21 (phenotypic selection) and 0.01 to 0.45 (selection based on EBV). The lowest genetic gain was realized for the rotational use of village bulls among villages within groups. Through the rotational use of village bulls, however, a higher genetic variance was maintained than in the farm and village bull breeding scenarios. We concluded that a village bull breeding program with selection based on EBV of young bulls was the most promising breeding design for achieving the breeding goal. Further studies are needed to assess the organizational feasibility of such a breeding program to ensure the participation of smallholder producers and its sustainability.     

Economic values of production and functional traits,including residual feed intake,in Finnish milk production

Improving the feed efficiency of dairy cattle has a substantial effect on the economic efficiency and on the reduction of harmful environmental effects of dairy production through lower feeding costs and emissions from dairy farming. To assess the economic importance of feed efficiency in the breeding goal for dairy cattle, the economic values for the current breeding goal traits and the additional feed efficiency traits for Finnish Ayrshire cattle under production circumstances in 2011 were determined. The derivation of economic values was based on a bioeconomic model in which the profit of the production system was calculated, using the generated steady state herd structure. Considering beef production from dairy farms, 2 marketing strategies for surplus calves were investigated: (A) surplus calves were sold at a young age and (B) surplus calves were fattened on dairy farms. Both marketing strategies were unprofitable when subsidies were not included in the revenues. When subsidies were taken into account, a positive profitability was observed in both marketing strategies. The marginal economic values for residual feed intake (RFI) of breeding heifers and cows were −25.5 and −55.8 €/kg of dry matter per day per cow and year, respectively. The marginal economic value for RFI of animals in fattening was −29.5 €/kg of dry matter per day per cow and year. To compare the economic importance among traits, the standardized economic weight of each trait was calculated as the product of the marginal economic value and the genetic standard deviation; the standardized economic weight expressed as a percentage of the sum of all standardized economic weights was called relative economic weight. When not accounting for subsidies, the highest relative economic weight was found for 305-d milk yield (34% in strategy A and 29% in strategy B), which was followed by protein percentage (13% in strategy A and 11% in strategy B). The third most important traits were calving interval (9%) and mature weight of cows (11%) in strategy A and B, respectively. The sums of the relative economic weights over categories for RFI were 6 and 7% in strategy A and B, respectively. Under production conditions in 2011, the relative economic weights for the studied feed efficiency traits were low. However, it is possible that the relative importance of feed efficiency traits in the breeding goal will increase in the future due to increasing requirements to mitigate the environmental impact of milk production.     

Invited review: Genomic selection in dairy cattle: Progress and challenges

A new technology called genomic selection is revolutionizing dairy cattle breeding. Genomic selection refers to selection decisions based on genomic breeding values (GEBV). The GEBV are calculated as the sum of the effects of dense genetic markers, or haplotypes of these markers, across the entire genome, thereby potentially capturing all the quantitative trait loci (QTL) that contribute to variation in a trait. The QTL effects, inferred from either haplotypes or individual single nucleotide polymorphism markers, are first estimated in a large reference population with phenotypic information. In subsequent generations, only marker information is required to calculate GEBV. The reliability of GEBV predicted in this way has already been evaluated in experiments in the United States, New Zealand, Australia, and the Netherlands. These experiments used reference populations of between 650 and 4,500 progeny-tested Holstein-Friesian bulls, genotyped for approximately 50,000 genome-wide markers. Reliabilities of GEBV for young bulls without progeny test results in the reference population were between 20 and 67%. The reliability achieved depended on the heritability of the trait evaluated, the number of bulls in the reference population, the statistical method used to estimate the single nucleotide polymorphism effects in the reference population, and the method used to calculate the reliability. A common finding in 3 countries (United States, New Zealand, and Australia) was that a straightforward BLUP method for estimating the marker effects gave reliabilities of GEBV almost as high as more complex methods. The BLUP method is attractive because the only prior information required is the additive genetic variance of the trait. All countries included a polygenic effect (parent average breeding value) in their GEBV calculation. This inclusion is recommended to capture any genetic variance not associated with the markers, and to put some selection pressure on low-frequency QTL that may not be captured by the markers. The reliabilities of GEBV achieved were significantly greater than the reliability of parental average breeding values, the current criteria for selection of bull calves to enter progeny test teams. The increase in reliability is sufficiently high that at least 2 dairy breeding companies are already marketing bull teams for commercial use based on their GEBV only, at 2 yr of age. This strategy should at least double the rate of genetic gain in the dairy industry. Many challenges with genomic selection and its implementation remain, including increasing the accuracy of GEBV, integrating genomic information into national and international genetic evaluations, and managing long-term genetic gain.     

Weighting factors of sire daughter information in international genetic evaluations

International genetic evaluations of dairy bulls are currently based on national genetic evaluation results. Total number of daughters in a country is used to weight national information, but may not optimally reflect the precision of a sire's daughter contribution to international genetic evaluations. This study investigates the impact of alternative weighting factors on international evaluation results. A conventional progeny test scheme was simulated for two dairy cattle populations, with semen exchange at a fixed rate after each generation. True breeding values for both populations were generated as bivariate normal deviates. Each cow had three lactation records in one country only. After 10 generations of selection, all records were used in national breeding value prediction. National breeding values of bulls were used as input to international evaluations. Seven different weighting factors were evaluated: 1) total number of daughters; 2) total number of lactations; 3) as (one) also adjusted for finite contemporary group size; 4) as (three) also adjusted for distribution of daughters over contemporary groups; 5) effective daughter contribution considering finite contemporary group size and correlation between repeated records; 6) as (five) also considering the reliability of the daughter dam evaluation; and 7) as (five) also considering the reliability of the daughter female ancestors' evaluations. Using the last two weighting factors yielded empirically unbiased estimates of sire variance. Using total number of daughters overestimated genetic variance by up to 7%. In general, international breeding values were marginally affected by choice of weighting factor. The effect was larger when different national evaluation models had been applied in the two countries. International reliabilities for the last two weighting factors were close to expectation, whereas using total number of daughters resulted in 1 to 4% negative bias. In practice, different countries apply a wide range of national evaluation models, and genetic ties may be weak between some populations, thereby increasing the potential effect of weighting factors on international comparisons. The weighting factor developed in this study, which considers contemporary group structure, correlation between repeated records, and reliability of dams of daughters, should replace total number of daughters in international genetic evaluations of dairy sires.     

Retrospective analysis of the accuracy of conversion equations and multiple-trait, across-country evaluations of Holstein bulls used internationally

In 1995, the multiple-trait across country genetic evaluation procedure replaced regression-based conversion equations as the preferred method for international genetic comparisons of dairy bulls. In the present study, February 1999 estimated breeding values of 632 foreign Holstein bulls that were used in Canada, Germany, Italy, The Netherlands, Sweden, and the US were compared with January 1995 predictions from home country data only. January 1995 predicted breeding values for each importing country were calculated using three methods: the multiple-trait, across-country evaluation procedure; conversion equations based on the multiple-trait, across-country evaluations; and conversion equations based on the Wilmink method. Mean correlations between 1999 estimated breeding values in the importing countries and 1995 predictions from international data were from 0.76 to 0.81 for all methods. The multiple-trait, across-country evaluation procedure is expected to lead to selection of different bulls, because bulls were allowed to be ranked differently in each country, but no significant increase in accuracy of selection was observed. The lack of improvement in accuracy of prediction was most likely due to limitations in data structure. International genetic comparisons are largely driven by data from a relatively small number of evaluated bulls with exported semen. Data from siblings and more distant relatives provide only weak, indirect genetic links between countries, and inclusion of such data seems to provide a minimal improvement in accuracy. Limitations in data structure might be alleviated by methods that define environments by climate or management factors rather than country borders.     

Derivation and calculation of approximate reliabilities and daughter yield-deviations of a random regression test-day model for genetic evaluation of dairy cattle

Test-day milk, fat, protein yield, and somatic cell score (SCS) were analyzed separately using data from the first 3 lactations and a random regression model. Data used in the model were from Austria, Germany, and Luxembourg and from Holstein, Red, and Jersey dairy cattle. For reliability approximation, a multiple-trait effective daughter contribution (MTEDC) method was developed under general multiple trait models, including random regression test-day models, by extending the single-trait daughter equivalents concept. The MTEDC was applied to the very large dairy population, with about 15.5 million animals. The calculation of reliabilities required less computer memory than the corresponding iteration program and a significantly lower computing time equivalent to 24 rounds of iteration. A formula for daughter-yield deviations was derived for bulls under multiple-trait models. Reliability associated with daughter-yield deviations was approximated using the MTEDC method. Both the daughter-yield deviation formula and associated reliability method were verified in a simulation study using the random regression test-day model. Correlations of lactation daughter-yield deviations with estimated breeding values calculated from a routine genetic evaluation were 0.996 for all bulls and 0.95 for young bulls having only daughters with short lactations.     

Short communication: Multivariate outlier detection for routine Nordic dairy cattle genetic evaluation in the Nordic Holstein and Red population

It is of practical importance to ensure the data quality from a milk-recording system before use for genetic evaluation. A procedure was developed for detection of multivariate outliers based on an approximation for Mahalanobis distance and was implemented in the Nordic Holstein and Red population. The general target of this procedure is based on the Nordic Cattle Genetic Evaluation yield model, which is a 9-trait model for milk, protein, and fat in the first 3 lactations. The procedure is based on the phenotypic correlation structure as a function of days in milk (DIM) and on computation of trait means and standard deviations within a production year, lactation, and DIM. For each record in the data, a Mahalanobis distance value was computed based on the trait mean and the covariance matrix for the actual production year, lactation, and DIM. A set of cutoff values, ranging from 10 to 100 with steps of 10, for discarding multivariate outliers was investigated. Prediction accuracy was calculated as the Pearson correlations between estimated breeding values predicted by full data set and estimated breeding values predicted by reduced data set for cows without records in the reduced data set and with 1 or more records deleted due to the editing rules on Mahalanobis distance. The results showed that, averaged over all scenarios, gains of 0.005 to 0.048 on prediction accuracy have been obtained by deleting the multivariate outliers. The improvements were more profound for progeny of young bulls compared with progeny of proven bulls. It is easy to implement this multivariate outlier-detection procedure in the routine genetic evaluation for different dairy cattle breeds; however, an optimal cutoff value for Mahalanobis distance needs to be defined to achieve an acceptable compromise between genetic evaluation accuracy and data deletion.     

Emphasis on resilience in dairy cattle breeding: Possibilities and consequences

This study aimed to investigate dairy cattle breeding goals with more emphasis on resilience. We simulated the consequences of increasing weight on resilience indicators and an assumed true resilience trait (TR). Two environments with different breeding goals were simulated to represent the variability of production systems across Europe. Ten different scenarios were stochastically simulated in a so-called pseudogenomic simulation approach. We showed that many modern dairy cattle breeding goals most likely have negative genetic gain for TR and promising resilience indicators such as the log-transformed, daily deviation from the lactation curve (LnVAR). In addition, there were many ways of improving TR by increasing the breeding goal weight of different resilience indicators. The results showed that adding breeding goal weight to resilience indicators, such as body condition score and LnVAR, could reverse the negative trend observed for resilience indicators. Loss in the aggregate genotype calculated with only current breeding goal traits was 12 to 76%. This loss was mainly due to a reduction in genetic gain in milk production. We observed higher genetic gain in beef production, fertility, and udder health when breeding for more resilience, but from an economical point of view, this was not high enough to compensate for the reduction in genetic gain in milk production. The highest genetic gain in TR was obtained when adding the highest breeding goal weight to LnVAR or TR, both with 0.29 genetic standard deviation units. The indicators we used, body condition score and LnVAR, can be measured on a large scale today with relatively cheap methods, which is crucial if we want to improve these traits through breeding. Economic values for resilience have to be estimated to find the most optimal breeding goal for a more resilient dairy cow in the future.     

Derivation of economic values for German dairy breeds by means of a bio-economic model—with special emphasis on functional traits

To assess the economic importance of breeding traits, economic values (EV) were derived for 3 German dairy cattle breeds: German Holstein (HOL), Angler (ANG), and Red and White Dual-Purpose (RDN). For that purpose, the stochastic bio-economic model SimHerd (SimHerd A/S, Viborg, Denmark) was used, which simulates the expected monetary gain in dairy herds. The EV was calculated as the alteration in average net return of the herd responding to a marginal change in the trait of interest. When deriving EV using SimHerd, economic consequences resulting from changes in the age structure of a dairy herd (i.e., structural herd effects) are considered. However, this requires the simulation of relationships between traits in the bio-economic model. To avoid double counting, the EV of a trait was corrected for effects from alterations in correlated traits using multiple regression analysis. The EV were derived for 23 traits in terms of production, conformation and workability, dairy health, calf survival, and reproduction performance. Furthermore, the relative economic importance of the breeding traits was calculated. Relative emphasis on production was between 39.9 and 44.4% in the breeds studied. Total costs per case of ketosis and metritis ranged from €167 to €196 and €173 to €182, respectively. Highest marginal EV of direct health traits were found for mastitis (€257 to €271 per case) and lameness (€270 to €310 per case). Consequently, relative emphasis on direct health traits was between 15.7 and 17.9%. The EV of reproduction performance showed largest differences among the cattle breeds. Overall relative emphasis on reproduction was 10.5% in HOL, 10.8% in ANG, and 6.5% in RDN. The relative economic importance of cow mortality ranged from 15.5 to 16.0% across the breeds. Collectively, the study showed the high economic importance of functional traits in the cattle breeds studied.     

Economic evaluation of the breeding goal for Norwegian Red dairy cattle

In this paper, a translog profit function was applied to estimate the economic values of the traits included in the breeding goal for Norwegian Red dairy cattle. The following 10 traits are included in the breeding goal: milk, meat, mastitis resistance, fertility, calving difficulties, stillbirths, other diseases, udder, temperament, and legs. An empirical implementation that locally approximates the unknown true profit function was suggested and estimated, taking farm heterogeneity into account. The model was applied to a panel data set of 3,259 Norwegian dairy farms over the period 1999 to 2003. Panel data, also called longitudinal or cross-sectional time-series data, are multiple cases (cows, farms, countries, etc.) observed over 2 or more time periods. The data set consisted of farm-level data, including production and economic data from the farm and the estimated breeding values for each cow's sire. The estimated economic values make it possible to test whether genetic selection has been profitable for the farmer, and the extent to which the currently used economic values were optimal during the period 1999 to 2003. Although the translog profit function is quite flexible, it is rather complex, and a simplified version of the model, a Cobb-Douglas profit function, was also estimated. However, the hypothesis that this simpler function adequately describes the data compared with the full translog model was rejected. Further, the hypothesis that the estimated breeding values are profit neutral was rejected (i.e., the hypothesis that there are no interactions between input and output prices on one hand and estimated breeding values on the other). These results indicated that selection not only leads to a parallel shift in profits, but also to changes in input use. Seven of the 10 traits had a significant effect on the farmers’ profit. The 3 traits that were not significant were calving difficulty, stillbirth, and other diseases. The results showed that the breeding program for Norwegian Red cattle has been fairly successful in improving farmers’ profits. However, a slight modification of the breeding goal, such as a reduction in the weights for stillbirths and other diseases and an increase in the weights for meat and temperament, would increase farm profits.     

Progeny testing and selection intensity for Holstein bulls in different countries

International Bull Evaluation Service (Interbull) Holstein evaluations from February 1995 through February 2003 were used to determine characteristics of progeny testing for Holstein bulls in Australia, Canada, Denmark, France, Germany, Italy, New Zealand, Sweden, The Netherlands, and the United States. The decision to graduate a bull from progeny test (PT) was assumed to have been made based on the second Interbull evaluation, and graduation was defined as the addition of 200 daughters in the period 2.5 to 4.5 yr later. Mean bull age at PT decision varied across countries by 12 mo. Mean numbers of herds and daughters ranged from 39 to 111 and 54 to 144, respectively. Countries with higher requirements for official evaluations generally had more herds and daughters but older bulls at PT decision. Mean estimated breeding values for yield traits of sires of tested bulls were most similar across countries for fat, differing by only 6.4 kg. The four countries highest for sire protein differed only by 1 kg; however, the range was 12 kg. Percentages of bulls graduated ranged from 4.4 to 14.7 across countries. Selection intensities (standardized selection differentials) tended to be about 1.0 for yield traits. Selection intensities for somatic cell score were generally unfavorable, reflecting selection for negatively correlated yield traits. Reflecting variation in national breeding goals, selection intensities for stature were positive for most countries and highly negative for New Zealand. Selection intensity for fore udder was generally the lowest among the traits examined. All but one country showed positive selection for udder support. These statistics permit comparison of the components of PT programs across country, illustrating possible opportunities for improvement.     

Assessment of the value of international genetic evaluations for yield in predicting domestic breeding values for foreign Holstein bulls

International genetic evaluations are a valuable source of information for decisions about the importation of (the semen of) foreign bulls. This study analyzed data from 6 countries (Australia, Canada, Italy, France, the Netherlands, and the United States) and compared international evaluations for production traits of foreign bulls (i.e., when no national daughter information was available) to their national breeding values in August 2009, which were based only on domestic daughters’ data. A total of 821 bulls with highly reliable estimated breeding values (EBV) for milk, fat, and protein yield were analyzed. No evidence of systematic over- or underestimation was found in most of the countries analyzed. Observed correlations between national and international evaluations were close to 0.9 and, for most countries, generally close to their expected values (calculated from national and international EBV reliabilities). In Italy, however, higher differences between observed and expected correlations and significant mean differences between EBV for more than one trait were observed in bulls progeny-tested in the United States and in other European countries (with differences up to 33.1% of the genetic standard deviation). These results were probably induced by a relatively recent change in the model for national evaluation. The findings in this study reflect a conservative estimate of the real value of international evaluations, as changes in methodologies in either the national or the international evaluations decreased the ability of past international evaluations to predict current national evaluations. Nevertheless, our results indicate that international evaluations based on foreign information for Holstein bulls were reasonably accurate predictors of the future national breeding values based only upon domestic daughters.     

Association analysis for young stock survival index with imputed whole-genome sequence variants in Nordic Holstein cattle

Identification of the genetic variants associated with calf survival in dairy cattle will aid in the elimination of harmful mutations from the cattle population and the reduction of calf and young stock mortality rates. We used de-regressed estimated breeding values for the young stock survival (YSS) index as response variables in a genome-wide association study with imputed whole-genome sequence variants. A total of 4,610 bulls with estimated breeding values were genotyped with the Illumina BovineSNP50 (Illumina, San Diego, CA) single nucleotide polymorphism (SNP) genotyping array. Genotypes were imputed to whole-genome sequence variants. After quality control, 15,419,550 SNP on 29 Bos taurus autosomes (BTA) were used for association analysis. A modified mixed-model association analysis was used for a genome scan, followed by a linear mixed-model analysis for selected genetic variants. We identified 498 SNP on BTA5 and BTA18 that were associated with the YSS index in Nordic Holstein. The SNP rs440345507 (Chr5:94721790) on BTA5 was the putative causal mutation affecting YSS. Two haplotype-based models were used to identify haplotypes with the largest detrimental effects on YSS index. For each association signal, 1 haplotype region with harmful effects and the lead associated SNP were identified. Detected haplotypes on BTA5 and BTA18 explained 1.16 and 1.20%, respectively, of genetic variance for the YSS index. We examined whether YSS quantitative trait loci (QTL) on BTA5 and BTA18 were associated with stillbirth. YSS QTL on BTA18 overlapped a QTL region for stillbirth, but most likely 2 different causal variants were responsible for these 2 QTL. Four component traits of the YSS index, defined by sex and age, were analyzed separately by the modified mixed-model approach. The same genomic regions were associated with both bull and heifer calf mortality. Several genes (EPS8, LOC100138951, and KLK family genes) contained a lead associated SNP or were included in haplotypes with large detrimental effects on YSS in Nordic Holstein cattle.