Bayesian inference of genetic parameters for test-day milk yield, milk quality traits, and somatic cell score in Burlina cows

The aim of the study was to infer (co)variance components for daily milk yield, fat and protein contents, and somatic cell score (SCS) in Burlina cattle (a local breed in northeast Italy). Data consisted of 13576 monthly test-day records of 666 cows (parities 1 to 8) collected in 10 herds between 1999 and 2009. Repeatability animal models were implemented using Bayesian methods. Flat priors were assumed for systematic effects of herd test date, days in milk, and parity, as well as for permanent environmental, genetic, and residual effects. On average, Burlina cows produced 17.0 kg of milk per day, with 3.66 and 3.33% of fat and protein, respectively, and 358000 cells per mL of milk. Marginal posterior medians (highest posterior density of 95%) of heritability were 0.18 (0.09–0.28), 0.28 (0.21–0.36), 0.35 (0.25–0.49), and 0.05 (0.01–0.11) for milk yield, fat content, protein content, and SCS, respectively. Marginal posterior medians of genetic correlations between the traits were low and a 95% Bayesian confidence region included zero, with the exception of the genetic correlation between fat and protein contents. Despite the low number of animals in the population, results suggest that genetic variance for production and quality traits exists in Burlina cattle.     

Milk yield persistency in Brazilian Gyr cattle based on a random regression model

With the objective of evaluating measures of milk yield persistency, 27,000 test-day milk yield records from 3362 first lactations of Brazilian Gyr cows that calved between 1990 and 2007 were analyzed with a random regression model. Random, additive genetic and permanent environmental effects were modeled using Legendre polynomials of order 4 and 5, respectively. Residual variance was modeled using five classes. The average lactation curve was modeled using a fourth-order Legendre polynomial. Heritability estimates for measures of persistency ranged from 0.10 to 0.25. Genetic correlations between measures of persistency and 305-day milk yield (Y305) ranged from -0.52 to 0.03. At high selection intensities for persistency measures and Y305, few animals were selected in common. As the selection intensity for the two traits decreased, a higher percentage of animals were selected in common. The average predicted breeding values for Y305 according to year of birth of the cows had a substantial annual genetic gain. In contrast, no improvement in the average persistency breeding value was observed. We conclude that selection for total milk yield during lactation does not identify bulls or cows that are genetically superior in terms of milk yield persistency. A measure of persistency represented by the sum of deviations of estimated breeding value for days 31 to 280 in relation to estimated breeding value for day 30 should be preferred in genetic evaluations of this trait in the Gyr breed, since this measure showed a medium heritability and a genetic correlation with 305-day milk yield close to zero. In addition, this measure is more adequate at the time of peak lactation, which occurs between days 25 and 30 after calving in this breed.     

Relation of Milk Production Loss to Milk Somatic Cell Count

Milk production loss was studied in relation to increased somatic cell count (SCC). Available data were weekly test-day milk yields and SCC (in 1,000 cells/ml), and mastitis incidences. In total, 18,131 records from 274 cows were used. Production loss was determined for test-day kg milk, kg protein, and kg energy-corrected milk. Least-squares analysis of variance was used to estimate the direct effect of Log10(SCC) on production. The recorded measures of production were first corrected for fixed effects, with adjustment factors estimated from a healthy data-set. The average daily milk yield was 19.7 kg/day in first lactation and 22.0 in later lactations. The geometric mean of SCC was 63.1 in first lactation and 107.2 in later lactations. The incidence of clinical mastitis treated by a veterinarian was 19.8% of the lactations-at-risk. Linear relationships were found between the production parameters and Log10(SCC). Quadratic and cubic effects were evaluated, but were found to contribute little to the overall fit of the models. The individual milk yield loss was 1.29 kg/day for each unit increase in Log10(SCC) for cows in first lactation. Milk yield decreased by 2.04 kg/day per unit Log10(SCC) for older cows. Corresponding values for protein yield were 0.042 and 0.067 kg/day for first and later lactations, respectively.     

Random regression models using different functions to model test-day milk yield of Brazilian Holstein cows

We analyzed 152,145 test-day records from 7317 first lactations of Holstein cows recorded from 1995 to 2003. Our objective was to model variations in test-day milk yield during the first lactation of Holstein cows by random regression model (RRM), using various functions in order to obtain adequate and parsimonious models for the estimation of genetic parameters. Test-day milk yields were grouped into weekly classes of days in milk, ranging from 1 to 44 weeks. The contemporary groups were defined as herd-test-day. The analyses were performed using a single-trait RRM, including the direct additive, permanent environmental and residual random effects. In addition, contemporary group and linear and quadratic effects of the age of cow at calving were included as fixed effects. The mean trend of milk yield was modeled with a fourth-order orthogonal Legendre polynomial. The additive genetic and permanent environmental covariance functions were estimated by random regression on two parametric functions, Ali and Schaeffer and Wilmink, and on B-spline functions of days in milk. The covariance components and the genetic parameters were estimated by the restricted maximum likelihood method. Results from RRM parametric and B-spline functions were compared to RRM on Legendre polynomials and with a multi-trait analysis, using the same data set. Heritability estimates presented similar trends during mid-lactation (13 to 31 weeks) and between week 37 and the end of lactation, for all RRM. Heritabilities obtained by multi-trait analysis were of a lower magnitude than those estimated by RRM. The RRMs with a higher number of parameters were more useful to describe the genetic variation of test-day milk yield throughout the lactation. RRM using B-spline and Legendre polynomials as base functions appears to be the most adequate to describe the covariance structure of the data.     

Association of the DGAT1 gene polymorphism in stud bulls with milk productivity in cows

The substitution of lysine for alanine (K232A) in the acyl-CoA-diacylglycerol acyltransferase, which is encoded by the DGAT1 gene, was tested for the significance for breeding evaluation of stud bulls of the holsteinized Black-and-White breed. The breeding value (deviation of milk productivity in daughters compared with cows of the same age) was estimated by the DYDC (daughter yield deviation to contemporary) method with modification. The frequency of allele 232K in the bulls examined was 0.28, lower than in Holstein bulls (0.4-0.6). The greatest effect of the A232K substitution was observed for the percent fat content in milk (1.4 sigma) and milk yield (0.76 sigma), and a lower effect was established for the milk protein yield (0.47 sigma) and percent protein content in milk (0.44 sigma). In the case of milk fat yield, the effect was nonsignificant. A method was proposed for converting the data on fat yield in order to obtain significant results in this case as well. The effect of the A232K substitution was estimated at 154 kg for milk yield, 2.8 kg for protein yield, 0.079% for percent fat content, and 0.015% for percent protein content. The results are discussed in the context of multiple pleiotropic effects of the K232A substitution in the DGAT1 gene. It is proposed that the K232A substitution of the DGAT1 gene be used as a golden standard in comparisons of the effect on milk productivity for the total gene set. This approach will allow a meta-analysis of the gene effects in spite of the differences in dairy cattle breeds and methods used to analyze their breeding value. In view of more than 30-year experience of using sperm of Holstein stud bulls, including those bred in North America, it was noted that the effect of the A232K substitution on milk productivity traits agreed well with the data reported for the North American commercial population of Holstein cattle.     

Genetic variation in milk urea nitrogen concentration of dairy cattle and its implications for reducing urinary nitrogen excretion

Nitrogen (N) leached into groundwater from urine patches of cattle grazing in situ is an environmental problem in pasture-based dairy industries. One potential mitigation is to breed cattle for lower urinary nitrogen (UN) excretion. Urinary nitrogen is difficult to measure, while milk urea nitrogen concentration (MUN) is relatively easy to measure. For animals fed diets of differing N content in confinement, MUN is moderately heritable and is positively related to UN. However, there is little information on the heritability of MUN, and its relationship with other traits such as milk yield and composition, for animals grazing fresh pasture. Milk urea nitrogen concentration data together with milk yield, fat, protein and lactose composition and somatic cell count was collected from 133 624 Holstein-Friesian (HF), Jersey (J) and HF×J (XBd) cows fed predominantly pasture over three full lactations and one part lactation. Mean MUN was 14.0; and 14.4, 13.2 and 13.9 mg/dl for HF, J and XBd cows, respectively. Estimates of heritability of MUN were 0.22 using a repeatability model that fitted year-of-lactation by month-of-lactation by cow-age with days-in-milk within month-of-lactation and cow-age, and 0.28 using a test-day model analysis with Gibbs sampling methods. Sire breeding values (BVs) ranged from −2.8 to +3.2 indicating that MUN could be changed by selection. The genetic correlation between MUN and percent true protein in milk was −0.22; −0.29 for J cows and −0.16 for HF cows. Should the relationship between MUN and UN observed in dietary manipulation studies hold similarly when MUN is manipulated by genetic selection, UN excretion could be reduced by 6.6 kg/cow per year in one generation of selection using sires with low MUN BVs. Although J cows had lower MUN than HF, total herd UN excretion may be similar for the same fixed feed supply because more J cows are required to utilise the available feed. The close relationship between blood plasma urea N concentration and MUN may enable early selection of bulls to breed progeny that excrete less UN.     

Effect of leptin gene polymorphism on the breeding value of milk production traits in Iranian Holstein

New molecular techniques focused on genome analysis, open new possibilities for more accurate evaluation of economiclly important traits in farm animals. Milk production traits are typical quantitative characteristics controlled by a number of genes. Mutations in their sequences may alter animal performance as well as their breeding values. In this study, we investigated the effect of Kpn2I restriction fragment length polymorphisms in the leptin gene, on bull breeding values for milk yield, fat, and protein yield, and their percentage. In order to test for an association between the leptin single-nucleotide polymorphism in exon 2 and milk productivity, we genotyped 134 Iranian Holstein bulls. Breeding values for milk-related traits (milk yield, fat, and protein yield and percentage) were estimated using the BLUP based on an animal model. The effect of the genotypes of Kpn2I polymorphism on the breeding values for milk-related traits was examined using least square methods. The T allele frequency was 0.425. Genotypes were distributed according to the Hardy-Weinberg equilibrium. Bulls with TT genotype had higher milk, fat and protein yield compared with TC and CC bulls (P < 0.05). Bulls with CC genotype had higher protein percentage compared with TT and TC bulls (P < 0.05). The association between leptin polymorphism with milk production traits suggests that this marker may be useful for selection based on molecular information.     

Bovine renal lipofuscinosis: Prevalence, genetics and impact on milk production and weight at slaughter in Danish cattle

Background

Bovine renal lipofuscinosis (BRL) is an incidental finding in cattle at slaughter. Condemnation of the kidneys as unfit for human consumption was until recently considered the only implication of BRL. Recent studies have indicated a negative influence on the health of affected animals. The present study investigated the prevalence, genetics and effect of BRL on milk yield and weight at slaughter.

Methods

BRL status of slaughter cattle was recorded at four abattoirs during a 2-year-period. Data regarding breed, age, genetic descent, milk yield and weight at slaughter were extracted from the Danish Cattle Database. The prevalence of BRL was estimated stratified by breed and age-group. Furthermore, total milk yield, milk yield in last full lactation and weight at slaughter were compared for BRL-affected and non-affected Danish Holsteins and Danish Red cattle.

Results

433,759 bovines were slaughtered and 787 of these had BRL. BRL was mainly diagnosed in Danish Red, Danish Holstein and crossbreds. The age of BRL affected animals varied from 11 months to 13 years, but BRL was rarely diagnosed in cattle less than 2 years of age. The total lifelong energy corrected milk (ECM) yields were 3,136 and 4,083 kg higher for BRL affected Danish Red and Danish Holsteins, respectively. However, the median life span of affected animals was 4.9 months longer, and age-corrected total milk yield was 1,284 kg lower for BRL affected Danish Red cows. These cows produced 318 kg ECM less in their last full lactation. Weight at slaughter was not affected by BRL status. The cases occurred in patterns consistent with autosomal recessive inheritance and several family clusters of BRL were found. Analysis of segregation ratios demonstrated the expected ratio for Danish Red cattle, but not for Danish Holsteins.

Conclusion

The study confirmed that BRL is a common finding in Danish Holsteins and Danish Red cattle at slaughter. The disorder is associated with increased total milk yield due to a longer production life. However, a reduced milk yield was detected in the end of the production life in Danish Red. The study supports that BRL is inherited autosomal recessively in the Danish Red breed and Danish Holsteins, but with incomplete penetrance of the genotype in Danish Holsteins.     

Inclusion of cow records in genomic evaluations and impact on bias due to preferential treatment

Background

Today, genomic evaluations are an essential feature of dairy cattle breeding. Initially, genomic evaluation targeted young bulls but recently, a rapidly increasing number of females (both heifers and cows) are being genotyped. A rising issue is whether and how own performance of genotyped cows should be included in genomic evaluations. The purpose of this study was to assess the impact of including yield deviations, i.e. own performance of cows, in genomic evaluations.

Methods

Two different genomic evaluations were performed: one including only reliable daughter yield deviations of proven bulls based on their non-genotyped daughters, and one including both daughter yield deviations for males and own yield deviations for genotyped females. Milk yield, the trait most prone to preferential treatment, and somatic cell count, for which such a bias is very unlikely, were studied. Data consisted of two groups of animals from the three main dairy breeds in France: 11 884 elite females genotyped by breeding companies and 7032 cows genotyped for a research project (and considered as randomly selected from the commercial population).

Results

For several measures that could be related to preferential treatment bias, the elite group presented a different pattern of estimated breeding values for milk yield compared to the other combinations of trait and group: for instance, for milk yield, the average difference between estimated breeding values with or without own yield deviations was significantly different from 0 for this group. Correlations between estimated breeding values with or without yield deviations were lower for elite females than for randomly selected cows for milk yield but were very similar for somatic cell count.

Conclusions

This study demonstrated that including own milk performance of elite females leads to biased (over-estimated) genomic evaluations. Thus, milk production records of elite cows require specific treatment in genomic evaluation.     

Expression profile of genes associated with mastitis in dairy cattle

In order to characterize the expression of genes associated with immune response mechanisms to mastitis, we quantified the relative expression of the IL-2, IL-4, IL-6, IL-8, IL-10, IFN-γ and TNF- α genes in milk cells of healthy cows and cows with clinical mastitis. Total RNA was extracted from milk cells of six Black and White Holstein (BW) cows and six Gyr cows, including three animals with and three without mastitis per breed. Gene expression was analyzed by real-time PCR. IL-10 gene expression was higher in the group of BW and Gyr cows with mastitis compared to animals free of infection from both breeds (p < 0.05). It was also higher in BW Holstein animals with clinical mastitis (p < 0.001), but it was not significant when Gyr cows with and without mastitis were compared (0.05 < p < 0.10). Among healthy cows, BW Holstein animals tended to present a higher expression of all genes studied, with a significant difference for the IL-2 and IFN- γ genes (p < 0.001). For animals with mastitis no significant difference in gene expression was observed between the two breeds. These findings suggest that animals with mastitis develop a preferentially cell-mediated immune response. Further studies including larger samples are necessary to better characterize the gene expression profile in cows with mastitis.     

Evaluation of reference lactation length in Chios dairy sheep

Definition and establishment of a fixed reference lactation length could provide useful tools for on-farm comparison of ewes and flock management as well as genetic evaluations for the breeding programme. The objectives of this study were to (i) evaluate different reference lactation lengths for the Chios dairy sheep and (ii) define the most suitable reference length for the breed. A total of 260 042 test-day milk records from 24 474 ewes in 130 flocks collected between 2003 and 2014 were used; 15 different lactation lengths were evaluated ranging from 120 to 260 days, defined at 10-day intervals as reference for the Chios sheep. The evaluation criteria included: (a) heritability and repeatability of milk yield in each reference lactation, (b) genetic correlation of reference lactation milk yield with actual lactation milk yield and yield at first test-day record and (c) correlated response in reference lactation milk yield from selection based on first test-day milk yield. The latter emulates genetic gains achieved for milk yield based on early lactation selection. Heritability and repeatability estimates of reference lactation milk yield and genetic correlation with actual lactation yield favoured long reference lactations (180 to 230 days). On the contrary, correlation with first test-day record milk yield was higher for short lactations (120 to 170 days). Moreover, selection on first test-day record milk yield would lead to a correlated response in reference yield in 220 days equal to 85% of the highest estimate achieved in the maximum reference length of 260 days (190 days when only considering first lactation milk yield). Based on the results of the present study, an overall reference lactation length for the Chios breed of 220 days post-lambing and a first lactation reference length of 190 days post-lambing are recommended.     

Suitability of cross-bred cows for organic farms based on cross-breeding effects on production and functional traits

Data from 113 Dutch organic farms were analysed to determine the effect of cross-breeding on production and functional traits. In total, data on 33 788 lactations between January 2003 and February 2009 from 15 015 cows were available. Holstein–Friesian pure-bred cows produced most kg of milk in 305 days, but with the lowest percentages of fat and protein of all pure-bred cows in the data set. Cross-breeding Holstein dairy cows with other breeds (Brown Swiss, Dutch Friesian, Groningen White Headed, Jersey, Meuse Rhine Yssel, Montbéliarde or Fleckvieh) decreased milk production, but improved fertility and udder health in most cross-bred animals. In most breeds, heterosis had a significant effect (P < 0.05) on milk (kg in 305 days), fat and protein-corrected milk production (kg in 305 days) and calving interval (CI) in the favourable direction (i.e. more milk, shorter CI), but unfavourably for somatic cell count (higher cell count). Recombination was unfavourable for the milk production traits, but favourable for the functional traits (fertility and udder health). Farm characteristics, like soil type or housing system, affected the regression coefficients on breed components significantly. The effect of the Holstein breed on milk yield was twice as large in cubicle housing as in other housing systems. Jerseys had a negative effect on fertility only on farms on sandy soils. Hence, breed effects differ across farming systems in the organic farming and farmers can use such information to dovetail their farming system with the type of cow they use.     

Mediation analysis to estimate direct and indirect milk losses associated with bacterial load in bovine subclinical mammary infections

Milk losses associated with mastitis can be attributed to either effects of pathogens per se (i.e. direct losses) or to effects of the immune response triggered by the presence of mammary pathogens (i.e. indirect losses). Test-day milk somatic cell counts (SCC) and number of bacterial colony forming units (CFU) found in milk samples are putative measures of the level of immune response and of the bacterial load, respectively. Mediation models, in which one independent variable affects a second variable which, in turn, affects a third one, are conceivable models to estimate direct and indirect losses. Here, we evaluated the feasibility of a mediation model in which test-day SCC and milk were regressed toward bacterial CFU measured at three selected sampling dates, 1 week apart. We applied this method on cows free of clinical signs and with records on up to 3 test-days before and after the date of the first bacteriological samples. Most bacteriological cultures were negative (52.38%), others contained either staphylococci (23.08%), streptococci (9.16%), mixed bacteria (8.79%) or were contaminated (6.59%). Only losses mediated by an increase in SCC were significantly different from null. In cows with three consecutive bacteriological positive results, we estimated a decreased milk yield of 0.28 kg per day for each unit increase in log₂-transformed CFU that elicited one unit increase in log₂-transformed SCC. In cows with one or two bacteriological positive results, indirect milk loss was not significantly different from null although test-day milk decreased by 0.74 kg per day for each unit increase of log₂-transformed SCC. These results highlight the importance of milk losses that are mediated by an increase in SCC during mammary infection and the feasibility of decomposing total milk loss into its direct and indirect components.     

Strategic test-day recording regimes to estimate lactation yield in tropical dairy animals

Background

In developing dairy sectors, genetic improvement programs have limited resources and recording of herds is minimal. This study evaluated different methods to estimate lactation yield and sampling schedules with fewer test-day records per lactation to determine recording regimes that (1) estimate lactation yield with a minimal impact on the accuracy of selection and (2) optimise the available resources.

Methods

Using Sahiwal cattle as a tropical dairy breed example, weekly milk records from 464 cows were used in a simulation study to generate different shaped lactation curves. The daily milk yields from these simulated lactation curves were subset to equally spaced (weekly, monthly and quarterly) and unequally spaced (with four, five or six records per lactation) test-day intervals. Lactation yield estimates were calculated from these subsets using two methods: the test-interval method and Wood’s (Nature 216:164-165, 1967) lactation curve model. Using the resulting lactation yields, breeding values were predicted and comparisons were made between the sampling regimes and estimation methods.

Results

The results show that, based on the mean square error of prediction, use of Wood’s lactation curve model to estimate total yield was more accurate than use of the test-interval method. However, the differences in the ranking of animals were small, i.e. a 1 to 5% difference in accuracy. Comparisons between the different test-day sampling regimes showed that, with the same number of records per lactation (for example, quarterly and four test-days), strategically timed test-days can result in more accurate estimates of lactation yield than test-days at equal intervals.

Conclusions

An important outcome of these results is that combining Wood’s model for lactation yield estimation and as few as four, five or six strategically placed test-day records can produce estimates of lactation yield that are comparable with estimates based on monthly test-day records using the test-interval method. Furthermore, calculations show that although using fewer test-days results in a decrease in the accuracy of selection, it does provide an opportunity to progeny-test more sires. Thus, using strategically timed test-days and Wood’s model to estimate lactation yield, can lead to a more efficient use of the allocated resources.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-014-0078-0) contains supplementary material, which is available to authorized users.     

Effects of stage of pregnancy on variance components,daily milk yields and 305-day milk yield in Holstein cows,as estimated by using a test-day model

Pregnancy and calving are elements indispensable for dairy production, but the daily milk yield of cows decline as pregnancy progresses, especially during the late stages. Therefore, the effect of stage of pregnancy on daily milk yield must be clarified to accurately estimate the breeding values and lifetime productivity of cows. To improve the genetic evaluation model for daily milk yield and determine the effect of the timing of pregnancy on productivity, we used a test-day model to assess the effects of stage of pregnancy on variance component estimates, daily milk yields and 305-day milk yield during the first three lactations of Holstein cows. Data were 10 646 333 test-day records for the first lactation; 8 222 661 records for the second; and 5 513 039 records for the third. The data were analyzed within each lactation by using three single-trait random regression animal models: one model that did not account for the stage of pregnancy effect and two models that did. The effect of stage of pregnancy on test-day milk yield was included in the model by applying a regression on days pregnant or fitting a separate lactation curve for each days open (days from calving to pregnancy) class (eight levels). Stage of pregnancy did not affect the heritability estimates of daily milk yield, although the additive genetic and permanent environmental variances in late lactation were decreased by accounting for the stage of pregnancy effect. The effects of days pregnant on daily milk yield during late lactation were larger in the second and third lactations than in the first lactation. The rates of reduction of the 305-day milk yield of cows that conceived fewer than 90 days after the second or third calving were significantly (P<0.05) greater than that after the first calving. Therefore, we conclude that differences between the negative effects of early pregnancy in the first, compared with later, lactations should be included when determining the optimal number of days open to maximize lifetime productivity in dairy cows.     

Quantitative trait loci linked to PRNP gene controlling health and production traits in INRA 401 sheep

In this study, the potential association of PrP genotypes with health and productive traits was investigated. Data were recorded on animals of the INRA 401 breed from the Bourges-La Sapinière INRA experimental farm. The population consisted of 30 rams and 852 ewes, which produced 1310 lambs. The animals were categorized into three PrP genotype classes: ARR homozygous, ARR heterozygous, and animals without any ARR allele. Two analyses differing in the approach considered were carried out. Firstly, the potential association of the PrP genotype with disease (Salmonella resistance) and production (wool and carcass) traits was studied. The data used included 1042, 1043 and 1013 genotyped animals for the Salmonella resistance, wool and carcass traits, respectively. The different traits were analyzed using an animal model, where the PrP genotype effect was included as a fixed effect. Association analyses do not indicate any evidence of an effect of PrP genotypes on traits studied in this breed. Secondly, a quantitative trait loci (QTL) detection approach using the PRNP gene as a marker was applied on ovine chromosome 13. Interval mapping was used. Evidence for one QTL affecting mean fiber diameter was found at 25 cM from the PRNP gene. However, a linkage between PRNP and this QTL does not imply unfavorable linkage disequilibrium for PRNP selection purposes.     

The impact of a leptin gene SNP on beef calf weaning weights

Prior research indicates that a SNP at position 305 of exon 2 in the leptin gene affects milk production in dairy cows. Dairy cows with at least one copy of the T allele have been shown to have higher milk production than CC cows. If that effect carries over to beef breeds, it is reasonable to expect that CT and TT beef cows will wean heavier calves than CC beef cows. We tested this hypothesis for a herd of mixed breed cows using anova. Results indicated that both crossbred CT and TT beef cows wean significantly heavier beef calves than CC crossbred beef cows. A lack of observations generally hinders detection of significance in other breeds. However, two other comparisons were found to be significant. The results suggest further investigation into the link between leptin genotype and calf weaning weights. Aside from interest to animal scientists, these results have the potential to alter mating and replacement selection decisions by cow-calf producers, given the importance of weaning weights on profitability.     

Paratuberculosis: decrease in milk production of German Holstein dairy cows shedding Mycobacterium avium ssp. paratuberculosis depends on within-herd prevalence

Paratuberculosis impairs productivity of infected dairy cows because of reduced milk production and fertility and enhanced risk of culling. The magnitude of the milk yield depression in individual cows is influenced by factors such as parity, the stage of the disease and the choice of test used. The objectives of this case–control study were to substantiate the influence of the different levels of the within-herd prevalence (WHP) on individual milk yield of fecal culture (FC)-positive cows (FC+) compared with FC-negative herd-mates (FC−), and to estimate the magnitude of the deviation of the milk yield, milk components and somatic cell count (SCC) in an FC-based study. Of a total of 31 420 cows from 26 Thuringian dairy herds tested for paratuberculosis by FC, a subset of 1382 FC+ and 3245 FC− with milk recording data were selected as cases and controls, respectively. The FC− cows were matched for the same number and stage of lactation (±10 days in milk) as one FC+ from the same herd. Within a mixed model analysis using the fixed effects of Mycobacterium avium ssp. paratuberculosis (MAP) status, lactation number, days in milk, prevalence class of farm and the random effect of farm on milk yield per day (kg), the amount of fat and protein (mg/dl) and lactose (mg/dl) as well as the SCC (1000/ml) were measured. On the basis of least square means, FC+ cows had a lower test-day milk yield (27.7±0.6 kg) compared with FC− (29.0±0.6 kg), as well as a lower milk protein content and a slightly diminished lactose concentration. FC status was not associated with milk fat percentage or milk SCC. In FC+ cows, reduction in milk yield increased with increasing WHP. An interaction of FC status and farm was found for the test-day milk yield, and milk protein percentage, respectively. We conclude that the reduction in milk yield of FC+ cows compared with FC− herd-mates is significantly influenced by farm effects and depends on WHP class. Owners of MAP-positive dairy herds may benefit from the reduction in WHP not only by reducing number of infected individuals but also by diminishing the individual losses in milk production per infected cow, and therefore should establish control measures.     

The effect of missing marker genotypes on the accuracy of gene-assisted breeding value estimation: a comparison of methods

In livestock populations, missing genotypes on a large proportion of the animals is a major problem when implementing gene-assisted breeding value estimation for genes with known effect. The objective of this study was to compare different methods to deal with missing genotypes on accuracy of gene-assisted breeding value estimation for identified bi-allelic genes using Monte Carlo simulation. A nested full-sib half-sib structure was simulated with a mixed inheritance model with one bi-allelic quantitative trait loci (QTL) and a polygenic effect due to infinite number of polygenes. The effect of the QTL was included in gene-assisted BLUP either by random regression on predicted gene content, i.e. the number of positive alleles, or including haplotype effects in the model with an inverse IBD matrix to account for identity-by-descent relationships between haplotypes using linkage analysis information (IBD-LA). The inverse IBD matrix was constructed using segregation indicator probabilities obtained from multiple marker iterative peeling. Gene contents for unknown genotypes were predicted using either multiple marker iterative peeling or mixed model methodology. For both methods, gene-assisted breeding value estimation increased accuracies of total estimated breeding value (EBV) with 0% to 22% for genotyped animals in comparison to conventional breeding value estimation. For animals that were not genotyped, the increase in accuracy was much lower (0% to 5%), but still substantial when the heritability was 0.1 and when the QTL explained at least 15% of the genetic variance. Regression on predicted gene content yielded higher accuracies than IBD-LA. Allele substitution effects were, however, overestimated, especially when only sires and males in the last generation were genotyped. For juveniles without phenotypic records and traits measured only on females, the superiority of regression on gene content over IBD-LA was larger than when all animals had phenotypes. Missing gene contents were predicted with higher accuracy using multiple-marker iterative peeling than with using mixed model methodology, but the difference in accuracy of total EBV was negligible and mixed model methodology was computationally much faster than multiple iterative peeling. For large livestock populations it can be concluded that gene-assisted breeding value estimation can be practically best performed by regression on gene contents, using mixed model methodology to predict missing marker genotypes, combining phenotypic information of genotyped and ungenotyped animals in one evaluation. This technique would be, in principle, also feasible for genomic selection. It is expected that genomic selection for ungenotyped animals using predicted single nucleotide polymorphism gene contents might be beneficial especially for low heritable traits.