How does a beef × dairy calving affect the dairy cow's following lactation?

For beef semen usage on dairy cows, much of the research has focused on the performance of the crossbred calves, yet little focus has been given to the subsequent performance of the cow herself. This study aimed to evaluate the performance of dairy cows for milk yield, fertility, and survival traits after giving birth to beef × dairy crossbred calves and compare this with the performance after giving birth to purebred dairy calves. Further, we aimed to study if the effect of a difficult calving was the same regardless of whether the calf was purebred dairy or beef × dairy crossbred. Phenotypic records from 587,288 calving events from 1997 to 2020 were collected from the Swedish milk recording system from cows of the dairy breeds Swedish Red (SR) and Swedish Holstein. The sire beef breeds studied were Aberdeen Angus, Hereford (combined in category LHT), Charolais, Limousin, and Simmental (category HVY). Sixteen traits were defined and grouped in 3 categories: cumulative and 305-d milk, fat, and protein yield, daily milk yield, and 75-d milk yield as yield traits; calving to first insemination interval, calving to last insemination interval, first to last insemination interval, calving interval, and number of inseminations as fertility traits; and survival to 75 d or to next calving and lactation length as measures of survival. The data were analyzed for all traits for first and second parities separately using mixed linear models, with a focus on the estimates of cow breed by service sire breed combinations. All traits in parity 2 were adjusted for previous 305-d milk yield based on the expectation that low-yielding cows would more likely to be inseminated with beef semen. Overall, milk yield was lower after beef × dairy calvings compared with the purebred dairy calvings. The largest effects were found on cumulative yields and in second parity, with lower effects for yields early in lactation and yields in first parity. The largest decrease was 13 to 14 kg (0.12 phenotypic SD) for cumulative fat yield when breeding beef breed sires with purebred SR dams. For fertility traits, for most breed combinations, the effects were not large enough to be significant. Conversely, all beef × dairy crossbred combinations showed significantly lower results for survival to the next lactation, and mostly also for lactation length. There was some indication that dairy cows with beef × dairy calvings in parity 2 that were the result of maximum 2 inseminations in parity 1, had lower survival than corresponding calvings resulting from more than 2 inseminations. This could indicate that the former cows were marked for culling already when inseminated. There was generally an unfavorable effect of a difficult calving on all traits, however, there were almost no significant interactions between calving performance and dam by sire breed combination, and these interactions were never significant in first parity.

Multibreed genomic prediction using summary statistics and a breed-origin-of-alleles approach.

Because of an increasing interest in crossbreeding between dairy breeds in dairy cattle herds, farmers are requesting breeding values for crossbred animals. However, genomically enhanced breeding values are difficult to predict in crossbred populations because the genetic make-up of crossbred individuals is unlikely to follow the same pattern as for purebreds. Furthermore, sharing genotype and phenotype information between breed populations are not always possible, which means that genetic merit (GM) for crossbred animals may be predicted without the information needed from some pure breeds, resulting in low prediction accuracy. This simulation study investigated the consequences of using summary statistics from single-breed genomic predictions for some or all pure breeds in two- and three-breed rotational crosses, rather than their raw data. A genomic prediction model taking into account the breed-origin of alleles (BOA) was considered. Because of a high genomic correlation between the breeds simulated (0.62-0.87), the prediction accuracies using the BOA approach were similar to a joint model, assuming homogeneous SNP effects for these breeds. Having a reference population with summary statistics available from all pure breeds and full phenotype and genotype information from crossbreds yielded almost as high prediction accuracies (0.720-0.768) as having a reference population with full information from all pure breeds and crossbreds (0.753-0.789). Lacking information from the pure breeds yielded much lower prediction accuracies (0.590-0.676). Furthermore, including crossbred animals in a combined reference population also benefitted prediction accuracies in the purebred animals, especially for the smallest breed population.

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.

Cow characteristics associated with the variation in number of contacts between dairy cows.

In modern freestall barns where large groups of cows are housed together, the behavior displayed by herd mates can influence the welfare and production of other individuals. Therefore, understanding social interactions in groups of dairy cows is important to enhance herd management and optimize the outcomes of both animal health and welfare in the future. Many factors can affect the number of social contacts in a group. This study aimed to identify which characteristics of a cow are associated with the number of contacts it has with other group members in 2 different functional areas (feeding and resting area) to increase our understanding of the social behavior of dairy cows. Inside 2 herds housed in freestall barns with around 200 lactating cows each, cow positions were recorded with an ultra-wideband real-time location system collecting all cows' positions every second over 2 wk. Using the positioning data of the cows, we quantified the number of contacts between them, assuming that cows spending time in proximity to one another (within a distance of 2.5 m for at least 10 min per day) were interacting socially. We documented in which barn areas these interactions occurred and used linear mixed models to investigate if lactation stage, parity, breed, pregnancy status, estrus, udder health, and claw health affect the number of contacts. We found variation in the number of contacts a cow had between individuals in both functional areas. Cows in later lactation had more contacts in the feeding area than cows in early lactation. Furthermore, in one herd, higher parity cows had fewer contacts in the feeding area than first parity cows, and in the other herd, cows in third parity or higher had more contacts in the resting area. This study indicates that cow characteristics such as parity and days in milk are associated with the number of contacts a cow has daily to its herd mates and provides useful information for further research on social interactions of dairy cows.

Mating allocations in Nordic Red Dairy Cattle using genomic information.

In this study, we compared mating allocations in Nordic Red Dairy Cattle using genomic information. We used linear programming to optimize different economic scores within each herd, considering genetic level, semen cost, the economic impact of recessive genetic defects, and genetic relationships. We selected 9,841 genotyped females born in Denmark, Finland, or Sweden in 2019 for mating allocations. We used 2 different pedigree relationship coefficients, the first tracing the pedigree 3 generations back from the parents of the potential mating and the second based on all available pedigree information. We used 3 different genomic relationship coefficients, 1 SNP-by-SNP genomic relationship and 2 based on shared genomic segments. We found high correlations (≥0.83) between the pedigree and genomic relationship measures. The mating results showed that it was possible to reduce the different genetic relationships between parents with minimal effect on genetic level. Including the cost of known recessive genetic defects eliminated expression of genetic defects. It was possible to reduce genomic relationships between parents with pedigree measures, but it was best done with genomic measures. Linear programming maximized the economic score for all herds studied within seconds, which means that it is suitable for implementation in mating software to be used by advisors and farmers.

Genetic consequences of terminal crossbreeding, genomic test, sexed semen, and beef semen in dairy herds.

The development of breeding tools, such as genomic selection and sexed semen, has progressed rapidly in dairy cattle breeding during the past decades. In combination with beef semen, these tools are adopted increasingly at herd level. Dairy crossbreeding is emerging, but the economic and genetic consequences of combining it with the other breeding tools are relatively unknown. We investigated 5 different sexed semen schemes where 0, 50, and 90% of the heifers; 50% of the heifers + 25% of the first-parity cows; and 90% of the heifers + 45% of the first-parity cows were bred to sexed semen. The 5 schemes were combined in scenarios managing pure-breeding or terminal crossbreeding, including genomic testing of all newborn heifers or no testing, and keeping Swedish Red or Swedish Holstein as an initial breed. Thus, 40 scenarios were simulated, combining 2 stochastic simulation models: SimHerd Crossbred (operational returns) and ADAM (genetic returns). The sum of operational and genetic returns equaled the total economic return. Beef semen was used in all scenarios to limit the surplus of replacement heifers. Terminal crossbreeding implied having a nucleus of purebred females, where some were inseminated with semen of the opposite breed. The F1 crossbred females were inseminated with beef semen. The reproductive performance played a role in improving the benefit of any of the tools. The most considerable total economic returns were achieved when all 4 breeding tools were combined. For Swedish Holstein, the highest total economic return compared with a pure-breeding scenario, without sexed semen and genomic test, was achieved when 90% sexed semen was used in heifers and 45% sexed semen was used for first-parity cows combined with genomic test and crossbreeding (+€58, 33% crossbreds in the herd). The highest total economic return for Swedish Red compared with a pure-breeding scenario, without sexed semen and genomic test, was achieved when 90% sexed semen was used in heifers combined with genomic test and crossbreeding (+€94, 46% crossbreds in the herd). Terminal crossbreeding resulted in lower genetic returns across the herd compared with the corresponding pure-breeding scenarios but was compensated by a higher operational return.

Association of genomically enhanced and parent average breeding values with cow performance in Nordic dairy cattle.

This study compared the abilities of virgin heifer genomically enhanced breeding values (GEBV) and parent average breeding values (PA) to predict future cow performance. To increase confidence in genomic technology among farmers, a clear demonstration of the relationship between genomic predictions and future phenotypes is needed. We analyzed 12 different traits in first parity, including production, conformation, fertility, and other functional traits. Phenotype data were obtained from national milk recording schemes and breeding values from the Nordic Cattle Genetic Evaluation. Direct genomic breeding values were calculated using genomic BLUP and combined with traditional breeding values, using bivariate blending. The data covered 14,862 Red Dairy Cattle, 17,145 Holstein, and 7,330 Jersey genotyped virgin heifers born between 2013 and 2015 in Denmark, Finland, and Sweden. Phenotypes adjusted for systematic environmental effects were used as measures of cow performance. Two methods were used to compared virgin heifer GEBV and PA regarding their ability to predict future cow performance: (1) correlations between breeding values and adjusted phenotypes, (2) ranking cows into 4 quartiles for their virgin heifer GEBV or PA, and calculating actual cow performance for each quartile. We showed that virgin heifer GEBV predicted cow performance significantly better than PA for the vast majority of analyzed traits. The correlations with adjusted phenotypes were 38 to 136% higher for GEBV than for PA in Red Dairy Cattle, 42 to 194% higher for GEBV in Holstein, and 11 to 78% higher for GEBV in Jersey. The relative change between GEBV bottom and top quartiles compared with that between PA bottom and top quartiles ranged from 9 to 261% for RDC, 42 to 138% for Holstein, and 4 to 90% for Jersey. Hence, farmers in Denmark, Finland, and Sweden can have confidence in using genomic technology on their herds.

Economic consequences of dairy crossbreeding in conventional and organic herds in Sweden.

This study simulated the consequences of crossbreeding between Swedish Holstein and Swedish Red on herd dynamics and herd profitability under Swedish conditions. Two base herds were simulated using a stochastic herd simulation model, SimHerd Crossbred. The herds reflected average Swedish conventional and organic herds having purebred Swedish Holstein. For each base herd, 3 breeding strategies were simulated: pure-breeding, 2-breed terminal crossbreeding, and 2-breed rotational crossbreeding. The terminal crossbreeding strategy implied having a nucleus of Swedish Holstein and a proportion of F1 Swedish Red × Swedish Holstein crossbred cows within the same herd. The crossbreds in this herd did not produce replacement heifers but exclusively beef × dairy cross calves. Beef semen was also used in the pure-breeding (10-20% in cows) and the rotational crossbreeding (40% in cows) strategies to retain a limited surplus of replacement heifers. To ensure an adequate number of crossbreds in the terminal crossbreeding strategy, X-sorted sexed semen was used for insemination in all the purebred heifers. The outcome was 67% purebred and 31% F1 crossbreds in the herd. In addition, 31% heterosis was expressed compared with 67% heterosis expressed using a 2-breed rotational crossbreeding strategy. Compared with the pure-breeding strategy, crossbreeding increased the annual contribution margin per cow by €20 to €59, with the rotational crossbreeding strategy creating the largest profitability. The increased profitability was mainly due to improved functional traits, especially fertility. For the conventional production system, the replacement rate was 39.3% in the pure-breeding strategy and decreased to 35.8 and 30.1% in the terminal and rotational crossbreeding strategy, respectively. Similar changes happened in the organic production system. Additionally, the crossbreeding strategies earned €22 to €42 more annually per cow from selling live calves for slaughter due to the extended use of beef semen. Milk production was similar between pure-breeding and terminal crossbreeding, and only decreased 1 to 2% in rotational crossbreeding. These results show that crossbreeding between Swedish Holstein and Swedish Red can be profitable in both conventional and organic Swedish herds using the strategies we have simulated. However, some aspects remain to be investigated, such as the economically optimal breeding strategy, genetic improvement, and transition strategies.

Genetic parameters for noncoagulating milk, milk coagulation properties, and detailed milk composition in Swedish Red Dairy Cattle.

The rennet-induced coagulation ability of milk is important in cheese production. For Swedish Red Dairy Cattle (RDC), this ability is reduced because of a high prevalence of noncoagulating (NC) milk. In this study, we simultaneously combined genetic parameters for NC milk, milk coagulation properties, milk composition, physical traits, and milk protein composition. Our aim was to estimate heritability and genetic and phenotypic correlations for NC milk and 24 traits (milk coagulation properties, milk composition, physical traits, and milk protein composition). Phenotypes and ∼7,000 SNP genotypes were available for all 600 Swedish RDC. The genotypes were imputed from ∼7,000 SNP to 50,000 SNP. Variance components and genetic parameters were estimated with an animal model. In Swedish RDC, a moderate heritability estimate of 0.28 was found for NC milk. For the other 24 traits, heritability estimates ranged from 0.12 to 0.77 (standard errors from 0.08 to 0.18). A total of 300 phenotypic and genetic correlations were estimated. For phenotypic and genetic correlations, 172 and 95 were significant, respectively. In general, most traits showing significant genetic correlations also showed significant phenotypic correlations. In this study, phenotypic and genetic correlations with NC milk suggest that many correlations between traits exist, making it difficult to predict the real consequences on the composition of milk, if selective breeding is applied on NC milk. We speculate that some of these consequences may lead to changes in the composition of milk, most likely affecting its physical and organoleptic properties. However, our results suggest that κ-casein could be used as an indicator trait to predict the occurrence of NC milk at the herd level.

Effects of milk proteins and posttranslational modifications on noncoagulating milk from Swedish Red dairy cattle.

Milk that does not coagulate after rennet addition, also called noncoagulating (NC) milk, is unwanted in cheese production due to prolonged processing time. Amounts of whey and casein proteins, genetic variants, as well as posttranslational modifications (PTM) of proteins are all contributing factors in rennet-induced coagulation of milk. In this study, we conducted a wide-ranging investigation of milk proteins in milk samples from 616 Swedish Red dairy cattle using liquid chromatography-high resolution mass spectrometry. Relative concentration of proteins, genetic variants, and PTM were compared between NC milk and coagulating milk. The PTM investigated were phosphorylation of caseins and glycosylation of κ-casein. Several genetic variants and PTM were found, including rare phosphorylation variants of the αS-caseins. Genetic variants were found to effect the expressed amount of different proteins. Further, the effect of protein amounts and PTM on a binary NC milk trait was modeled using a generalized linear model. The model showed that NC milk significantly correlated with higher relative concentrations of α-lactalbumin and ß-casein and lower relative concentrations of ß-lactoglobulin and κ-casein. Regarding PTM of caseins, an effect on NC milk from a lower relative concentration of αS1-casein with 8 phosphate groups were found, even though an effect from total relative concentration of αS1-casein was not found. This study has provided insights into protein variants and PTM important for NC milk to improve this undesirable property.

Genetic heterogeneity of feed intake, energy-corrected milk, and body weight across lactation in primiparous Holstein, Nordic Red, and Jersey cows.

In this study, we aimed to estimate and compare the genetic parameters of dry matter intake (DMI), energy-corrected milk (ECM), and body weight (BW) as 3 feed efficiency-related traits across lactation in 3 dairy cattle breeds (Holstein, Nordic Red, and Jersey). The analyses were based on weekly records of DMI, ECM, and BW per cow across lactation for 842 primiparous Holstein cows, 746 primiparous Nordic Red cows, and 378 primiparous Jersey cows. A random regression model was applied to estimate variance components and genetic parameters for DMI, ECM, and BW in each lactation week within each breed. Phenotypic means of DMI, ECM, and BW observations across lactation showed to be in very similar patterns between breeds, whereas breed differences lay in the average level of DMI, ECM, and BW. Generally, for all studied breeds, the heritability for DMI ranged from 0.2 to 0.4 across lactation and was in a range similar to the heritability for ECM. The heritability for BW ranged from 0.4 to 0.6 across lactation, higher than the heritability for DMI or ECM. Among the studied breeds, the heritability estimates for DMI shared a very similar range between breeds, whereas the heritability estimates for ECM tended to be different between breeds. For BW, the heritability estimates also tended to follow a similar range between breeds. Among the studied traits, the genetic variance and heritability for DMI varied across lactation, and the genetic correlations between DMI at different lactation stages were less than unity, indicating a genetic heterogeneity of feed intake across lactation in dairy cattle. In contrast, BW was the most genetically consistent trait across lactation, where BW among all lactation weeks was highly correlated. Genetic correlations between DMI, ECM, and BW changed across lactation, especially in early lactation. Energy-corrected milk had a low genetic correlation with both DMI and BW at the beginning of lactation, whereas ECM was highly correlated with DMI in mid and late lactation. Based on our results, genetic heterogeneity of DMI, ECM, and BW across lactation generally was observed in all studied dairy breeds, especially for DMI, which should be carefully considered for the recording strategy of these traits. The genetic correlations between DMI, ECM, and BW changed across lactation and followed similar patterns between breeds.

Validation of consistency of Mendelian sampling variance.

Experiences from international sire evaluation indicate that the multiple-trait across-country evaluation method is sensitive to changes in genetic variance over time. Top bulls from birth year classes with inflated genetic variance will benefit, hampering reliable ranking of bulls. However, none of the methods available today enable countries to validate their national evaluation models for heterogeneity of genetic variance. We describe a new validation method to fill this gap comprising the following steps: estimating within-year genetic variances using Mendelian sampling and its prediction error variance, fitting a weighted linear regression between the estimates and the years under study, identifying possible outliers, and defining a 95% empirical confidence interval for a possible trend in the estimates. We tested the specificity and sensitivity of the proposed validation method with simulated data using a real data structure. Moderate (M) and small (S) size populations were simulated under 3 scenarios: a control with homogeneous variance and 2 scenarios with yearly increases in phenotypic variance of 2 and 10%, respectively. Results showed that the new method was able to estimate genetic variance accurately enough to detect bias in genetic variance. Under the control scenario, the trend in genetic variance was practically zero in setting M. Testing cows with an average birth year class size of more than 43,000 in setting M showed that tolerance values are needed for both the trend and the outlier tests to detect only cases with a practical effect in larger data sets. Regardless of the magnitude (yearly increases in phenotypic variance of 2 or 10%) of the generated trend, it deviated statistically significantly from zero in all data replicates for both cows and bulls in setting M. In setting S with a mean of 27 bulls in a year class, the sampling error and thus the probability of a false-positive result clearly increased. Still, overall estimated genetic variance was close to the parametric value. Only rather strong trends in genetic variance deviated statistically significantly from zero in setting S. Results also showed that the new method was sensitive to the quality of the approximated reliabilities of breeding values used in calculating the prediction error variance. Thus, we recommend that only animals with a reliability of Mendelian sampling higher than 0.1 be included in the test and that low heritability traits be analyzed using bull data sets only.

Detection of evaluation bias caused by genomic preselection.

The aim of this simulation study was to investigate whether it is possible to detect the effect of genomic preselection on Mendelian sampling (MS) means or variances obtained by the MS validation test. Genomic preselection of bull calves is 1 additional potential source of bias in international evaluations unless adequately accounted for in national evaluations. Selection creates no bias in traditional breeding value evaluation if the data of all animals are included. However, this is not the case with genomic preselection, as it excludes culled bulls. Genomic breeding values become biased if calculated using a multistep procedure instead of, for example, a single-step method. Currently, about 60% of the countries participating in international bull evaluations have already adopted genomic selection in their breeding schemes. The data sent for multiple across-country evaluation can, therefore, be very heterogeneous, and a proper validation method is needed to ensure a fair comparison of the bulls included in international genetic evaluations. To study the effect of genomic preselection, we generated a total of 50 replicates under control and genomic preselection schemes using the structures of the real data and pedigree from a medium-size cow population. A genetic trend of 15% of the genetic standard deviation was created for both schemes. In carrying out the analyses, we used 2 different heritabilities: 0.25 and 0.10. From the start of genomic preselection, all bulls were genomically preselected. Their MS deviations were inflated with a value corresponding to selection of the best 10% of genomically tested bull calves. For cows, the MS deviations were unaltered. The results revealed a clear underestimation of bulls' breeding values (BV) after genomic preselection started, as well as a notable deviation from zero both in true and estimated MS means. The software developed recently for the MS validation test already produces yearly MS means, and they can be used to devise an appropriate test. Mean squared true MS of genomically preselected bulls was clearly inflated. After correcting for the simulated preselection bias, the true genetic variance was smaller than the parametric value used to simulate BV, and also below the variance based on the estimated BV. Based on this study, the lower the trait's heritability, the stronger the bias in estimated BV and MS means and variances. Daughters of genomically preselected bulls had higher true and estimated BV compared with the control scheme and only slightly elevated MS means, but no effect on genetic variances was observed.

Bivariate threshold models for genetic evaluation of susceptibility to and ability to recover from mastitis in Danish Holstein cows.

Mastitis in dairy cows is an unavoidable problem and genetic variation in recovery from mastitis, in addition to susceptibility, is therefore of interest. Genetic parameters for susceptibility to and recovery from mastitis were estimated for Danish Holstein-Friesian cows using data from automatic milking systems equipped with online somatic cell count measuring units. The somatic cell count measurements were converted to elevated mastitis risk, a continuous variable [on a (0-1) scale] indicating the risk of mastitis. Risk values >0.6 were assumed to indicate that a cow had mastitis. For each cow and lactation, the sequence of health states (mastitic or healthy) was converted to a weekly transition: 0 if the cow stayed within the same state and 1 if the cow changed state. The result was 2 series of transitions: one for healthy to diseased (HD, to model mastitis susceptibility) and the other for diseased to healthy (DH, to model recovery ability). The 2 series of transitions were analyzed with bivariate threshold models, including several systematic effects and a function of time. The model included effects of herd, parity, herd-test-week, permanent environment (to account for the repetitive nature of transition records from a cow) plus two time-varying effects (lactation stage and time within episode). In early lactation, there was an increased risk of getting mastitis but the risk remained stable afterwards. Mean recovery rate was 45% per lactation. Heritabilities were 0.07 [posterior mean of standard deviations (PSD) = 0.03] for HD and 0.08 (PSD = 0.03) for DH. The genetic correlation between HD and DH has a posterior mean of -0.83 (PSD = 0.13). Although susceptibility and recovery from mastitis are strongly negatively correlated, recovery can be considered as a new trait for selection.

Neglect of lactation stage leads to naive assessment of residual feed intake in dairy cattle.

Residual feed intake (RFI) is a candidate trait for feed efficiency in dairy cattle. We investigated the influence of lactation stage on the effect of energy sinks in defining RFI and the genetic parameters for RFI across lactation stages for primiparous dairy cattle. Our analysis included 747 primiparous Holstein cows, each with recordings on dry matter intake (DMI), milk yield, milk composition, and body weight (BW) over 44 lactation weeks. For each individual cow, energy-corrected milk (ECM), metabolic BW (MBW), and change in BW (ΔBW) were calculated in each week of lactation and were taken as energy sinks when defining RFI. Two RFI models were considered in the analyses; RFI model [1] was a 1-step RFI model with constant partial regression coefficients of DMI on energy sinks (ECM, MBW, and ΔBW) over lactation. In RFI model [2], data from 44 lactation weeks were divided into 11 consecutive lactation periods of 4 wk in length. The RFI model [2] was identical to model [1] except that period-specific partial regressions of DMI on ECM, MBW, and ΔBW in each lactation period were allowed across lactation. We estimated genetic parameters for RFI across lactation by both models using a random regression method. Using RFI model [2], we estimated the period-specific effects of ECM, MBW, and ΔBW on DMI in all lactation periods. Based on results from RFI model [2], the partial regression coefficients of DMI on ECM, MBW, and ΔBW differed across lactation in RFI. Constant partial regression coefficients of DMI on energy sinks over lactation was not always sufficient to account for the effects across lactation and tended to give roughly average information from all period-specific effects. Heritability for RFI over 44 lactation weeks ranged from 0.10 to 0.29 in model [1] and from 0.10 to 0.23 in model [2]. Genetic variance and heritability estimates for RFI from model [2] tended to be slightly lower and more stable across lactation than those from model [1]. In both models, RFI was genetically different over lactation, especially between early and later lactation stages. Genetic correlation estimates for RFI between early and later lactation tended to be higher when using model [2] compared with model [1]. In conclusion, partial regression coefficients of DMI on energy sinks differed across lactation when modeling RFI. Neglect of lactation stage when defining RFI could affect the assessment of RFI and the estimation of genetic parameters for RFI across lactation.

Genetic trends for fertility, udder health and protein yield in Swedish red cattle estimated with different models.

The aim of this study was to estimate and compare genetic trends in Swedish Red cattle using a full multiple-trait (MT) model and trait-group-wise models for female fertility, udder health and protein yield. Field data for maiden heifers from 1989 and cows with a first and second lactation between 1990 and 2007 were included. (Co)variance components were estimated prior to prediction of breeding values. The estimated genetic trends were clearly favourable for protein yield and udder conformation, and in most cases neutral to favourable for clinical mastitis and calving to first insemination. In maiden heifers, the trends were neutral for number of inseminations per service period. Unfavourable genetic trends were estimated for number of inseminations in the first two lactations, but the trends seemed less unfavourable from evaluations within trait groups compared with when using the full MT model. Excluding maiden heifer data affected genetic trends less than using trait-group-wise analyses instead of a full MT model. Unfavourable genetic trends in functional traits may be missed unless the traits are evaluated in a MT model including traits under strong selection.

Merging pedigree databases to describe and compare mating practices and gene flow between pedigree dogs in France, Sweden and the UK.

Merging pedigree databases across countries may improve the ability of kennel organizations to monitor genetic variability and health-related issues of pedigree dogs. We used data provided by the Société Centrale Canine (France), Svenska Kennelklubben (Sweden) and the Kennel Club (UK) to study the feasibility of merging pedigree databases across countries and describe breeding practices and international gene flow within the following four breeds: Bullmastiff (BMA), English setter (ESE), Bernese mountain dog (BMD) and Labrador retriever (LBR). After merging the databases, genealogical parameters and founder contributions were calculated according to the birth period, breed and registration country of the dogs. Throughout the investigated period, mating between close relatives, measured as the proportion of inbred individuals (considering only two generations of pedigree), decreased or remained stable, with the exception of LBR in France. Gene flow between countries became more frequent, and the origins of populations within countries became more diverse over time. In conclusion, the potential to reduce inbreeding within purebred dog populations through exchanging breeding animals across countries was confirmed by an improved effective population size when merging populations from different countries.

Genetic parameters for dry matter intake in primiparous Holstein, Nordic Red, and Jersey cows in the first half of lactation.

Dry matter intake (DMI) is a key component of feed efficiency in dairy cattle. In this study, we estimated genetic parameters of DMI over the first 24 lactation weeks in 3 dairy cattle breeds: Holstein, Nordic Red, and Jersey. In total, 1,656 primiparous cows (717 Holstein, 663 Nordic Red, and 276 Jersey) from Denmark, Finland, and Sweden were studied. For each breed, variance components, heritability, and repeatability for weekly DMI were estimated in 6 consecutive periods of the first 24 lactation weeks based on a repeatability animal model. Genetic correlations for DMI between different lactation periods were estimated using bivariate models. Based on our results, Holstein and Nordic Red cows had similar DMI at the beginning of lactation, but later in lactation Holstein cows had a slightly higher DMI than Nordic Red cows. In comparison, Jersey cows had a significantly lower DMI than the other 2 breeds within the first 24 lactation weeks. Heritability estimates for DMI ranged from 0.20 to 0.40 in Holsteins, 0.25 to 0.41 in Nordic Red, and 0.17 to 0.42 in Jerseys within the first 24 lactation weeks. Genetic and phenotypic variances for DMI varied along lactation within each breed and tended to be higher in the middle of lactation than at the beginning of the lactation. High genetic correlations were noted for DMI in lactation wk 5 to 24 in all 3 breeds, whereas DMI at early lactation (lactation wk 1 to 4) tended to be genetically different from DMI in the middle of lactation. The 3 breeds in this study might differ in their genetic variances for DMI, but the differences were not statistically significant in most of the studied periods. Breed differences for the genetic variance tended to be more obvious than for heritability. The potential breed differences in genetic variation for DMI should be considered in a future study using feed intake information from multiple breeds.

Estimates of genetic variance and variance of predicted genetic merits using pedigree or genomic relationship matrices in six Brown Swiss cattle populations for different traits.

The amount of variance captured in genetic estimations may depend on whether a pedigree-based or genomic relationship matrix is used. The purpose of this study was to investigate the genetic variance as well as the variance of predicted genetic merits (PGM) using pedigree-based or genomic relationship matrices in Brown Swiss cattle. We examined a range of traits in six populations amounting to 173 population-trait combinations. A main aim was to determine how using different relationship matrices affect variance estimation. We calculated ratios between different types of estimates and analysed the impact of trait heritability and population size. The genetic variances estimated by REML using a genomic relationship matrix were always smaller than the variances that were similarly estimated using a pedigree-based relationship matrix. The variances from the genomic relationship matrix became closer to estimates from a pedigree relationship matrix as heritability and population size increased. In contrast, variances of predicted genetic merits obtained using a genomic relationship matrix were mostly larger than variances of genetic merit predicted using pedigree-based relationship matrix. The ratio of the genomic to pedigree-based PGM variances decreased as heritability and population size rose. The increased variance among predicted genetic merits is important for animal breeding because this is one of the factors influencing genetic progress.

Breakeven prices for recording of indicator traits to reduce the environmental impact of milk production.

A breeding scheme using genomic selection and an indicator trait for environmental impact (EI) was studied to find the most effective recording strategy in terms of annual monetary genetic gain and breakeven price for the recording of indicator traits. The breakeven price shows the investment space for developing a recording system for an indicator trait. The breeding goal consisted of three traits ­ milk production, functional trait and environmental impact ­ with economic values of €83, €82 and €-83, respectively. The first scenario included only breeding goal traits and no indicator traits (NoIT). The other scenarios included all three breeding goal traits and one indicator trait (IT) for EI. The indicator traits were recorded on a large scale (stayability after first lactation and stature), medium scale (live weight and greenhouse gases (GHG) measured in the breath of the cow during milking) or small scale (residual feed intake and total enteric methane measured in a respiration chamber). In the scenario with stayability, the genetic gain in EI was over 11% higher than it was in NoIT. The breakeven price of recording stayability was €8 per record. Stayability is easy to record in the national milk recording system, and its use as an indicator trait for EI would not generate any additional recording costs. Therefore, stayability would be a good indicator trait to use to mitigate EI. The highest genetic gain in EI (23% higher compared to NoIT) was achieved when the GHG measured in the breath of the cow was used as indicator trait. The breakeven price for this indicator trait was €29 per record in the reference population. Ideally the recording of a specific indicator trait for EI would take place when: (i) the genetic correlation between the IT and EI is high; and (ii) the number of phenotypic records for the indicator trait is high enough to achieve a moderately high reliability of direct genomic values.