Genetic parameters for both a liver damage phenotype caused by and antibody response to phenotype in dairy and beef cattle

Genetic and genomic analysis of Belgian Blue's susceptibility for psoroptic mange

BACKGROUND

Psoroptic mange, caused by Psoroptes ovis mites, is affecting Belgian Blue cattle's welfare and production potential. The Belgian Blue cattle-known for its high degree of muscling, low feed conversion ratio and high beef quality-is highly susceptible for this disease.

RESULTS

In this study, we phenotyped 1975 Belgian Blue cattle from more than 100 different groups on commercial beef farms for their psoroptic mange susceptibility. Substantial individual differences were observed within these management groups, with lesion extent differences up to ± 15%. Animal models showed that estimated heritabilities were low for lesion extent and severe lesion extent (0.07 and 0.09, respectively) and 0.12 for the number of mites. A genome wide association study for mange susceptibility revealed signals on BTA6, BTA11, BTA15 and BTA24. In these regions, candidate genes GBA3, RAG2, and TRAF6 were identified.

CONCLUSIONS

Despite the challenges in phenotyping for psoroptic mange due to the timing of screening, the continuous evolution of lesions and different management conditions, we successfully conducted a study on the genetic susceptibility to psoroptic mange in Belgian Blue cattle. Our results clearly indicate that psoroptic mange is under polygenic control and the underlying candidate genes should be studied more thoroughly. This is the first study providing candidate genes for this complex disease. These results are already valuable for Belgian Blue breeding, however, further research is needed to unravel the architecture of this disease and to identify causal mutations.

Phenotypic and genetic associations between feeding behavior and carcass merit in crossbred growing cattle

In growing cattle, the phenotypic and genetic relationships between feeding behavior and both live animal ultrasound measures and subsequent carcass merit are generally poorly characterized. The objective of the current study was to quantify the phenotypic and genetic associations between a plethora of feeding behavior traits with both pre-slaughter ultrasound traits and post-slaughter carcass credentials in crossbred Bos taurus cattle. Carcass data were available on 3,146 young bulls, steers, and heifers, of which 2,795 and 2,445 also had records for pre-slaughter ultrasound muscle depth and intramuscular fat percentage, respectively; a total of 1,548 steers and heifers had information on all of the feeding behavior, ultrasound, and carcass traits. Young bulls were fed concentrates, while steers and heifers were fed a total mixed ration. Feeding behavior traits were defined based on individual feed events or meal events (i.e., individual feed events grouped into meals). Animal linear mixed models were used to estimate (co)variance components. Phenotypic correlations between feeding behavior and both ultrasound and carcass traits were generally weak and not different from zero, although there were phenotypic correlations of 0.40, 0.26, and 0.37 between carcass weight and feeding rate, energy intake per feed event, and energy intake per meal, respectively. Genetically, cattle that had heavier carcass weights, better carcass conformation, or a higher dressing percentage fed for a shorter time per day (genetic correlations [±SE] of -0.46 ± 0.12, -0.39 ± 0.11, and -0.50 ± 0.10, respectively). Genetic correlations of 0.43 ± 0.12 and 0.68 ± 0.13 were estimated between dressing difference (i.e., differential between live weight pre-slaughter and carcass weight) and energy intake per feed event and energy intake per meal, respectively. Neither intramuscular fat percentage measured on live animals nor carcass fat score (i.e., a measure of subcutaneous fat cover of the carcass) was genetically associated with any of the feeding behavior traits. The genetic associations between some feeding behavior traits and both ultrasound and carcass traits herein suggested that indirect responses in the feeding behavior of growing cattle from selection for improved carcass merit would be expected. Such changes in feeding behavior patterns in cattle may be reduced by measuring and including feeding behavior in a multiple-trait selection index alongside carcass traits.

Genetic variability in the feeding behavior of crossbred growing cattle and associations with performance and feed efficiency

The objectives of the present study were to estimate genetic parameters for several feeding behavior traits in growing cattle, as well as the genetic associations among and between feeding behavior and both performance and feed efficiency traits. An additional objective was to investigate the use of feeding behavior traits as predictors of genetic merit for feed intake. Feed intake and live-weight data on 6,088 growing cattle were used of which 4,672 had ultrasound data and 1,548 had feeding behavior data. Feeding behavior traits were defined based on individual feed events or meal events (where individual feed events were grouped into meals). Univariate and bivariate animal linear mixed models were used to estimate (co)variance components. Heritability estimates (± SE) for the feeding behavior traits ranged from 0.19 ± 0.08 for meals per day to 0.61 ± 0.10 for feeding time per day. The coefficient of genetic variation per trait varied from 5% for meals per day to 22% for the duration of each feed event. Genetically heavier cattle, those with a higher daily energy intake (MEI), or those that grew faster had a faster feeding rate, as well as a greater energy intake per feed event and per meal. Better daily feed efficiency (i.e., lower residual energy intake) was genetically associated with both a shorter feeding time per day and shorter meal time per day. In a validation population of 321 steers and heifers, the ability of estimated breeding values (EBV) for MEI to predict (adjusted) phenotypic MEI was demonstrated; EBVs for MEI were estimated using multi-trait models with different sets of predictor traits such as liveweight and/or feeding behaviors. The correlation (± SE) between phenotypic MEI and EBV for MEI marginally improved (P < 0.001) from 0.64 ± 0.03 to 0.68 ± 0.03 when feeding behavior phenotypes from the validation population were included in a genetic evaluation that already included phenotypic mid-test metabolic live-weight from the validation population. This is one of the largest studies demonstrating that significant exploitable genetic variation exists in the feeding behavior of young crossbred growing cattle; such feeding behavior traits are also genetically correlated with several performance and feed efficiency metrics. Nonetheless, there was only a marginal benefit to the inclusion of time-related feeding behavior phenotypes in a genetic evaluation for MEI to improve the precision of the EBVs for this trait.

Cross-sectional analyses of a national database to determine if superior genetic merit translates to superior dairy cow performance

Various studies have validated that genetic divergence in dairy cattle translates to phenotypic differences; nonetheless, many studies that consider the breeding goal, or associated traits, have generally been small scale, often undertaken in controlled environments, and they lack consideration for the entire suite of traits included in the breeding goal. Therefore, the objective of the present study was to fill this void, and in doing so, provide producers with confidence that the estimated breeding values (EBV) included in the breeding goal do (or otherwise) translate to desired changes in performance among commercial cattle; an additional outcome of such an approach is the identification of potential areas for improvements. Performance data on 536,923 Irish dairy cows (and their progeny) from 13,399 commercial spring-calving herds were used. Association analyses between the cow's EBV of each trait included in the Irish total merit index for dairy cows (which was derived before her own performance data accumulated) and her subsequent performance were undertaken using linear mixed models; milk production, fertility, calving, maintenance (i.e., liveweight), beef, health, and management traits were all considered in the analyses. Results confirm that excelling in EBV for individual traits, as well as on the total merit index, generally delivers superior phenotypic performance; examples of the improved performance for genetically elite animals include a greater yield and concentration of both milk fat and milk protein, despite a lower milk volume, superior reproductive performance, better survival, improved udder and hoof health, lighter cows, and fewer calving complications; all these gains were achieved with minimal to no effect on the beef merit of the dairy cow's progeny. The associated phenotypic change in each performance trait per unit change in its respective EBV was largely in line with the direction and magnitude of expectation, the exception being for calving interval. Per unit change in calving interval EBV, the direction of phenotypic response was as anticipated but the magnitude of the response was only half of what was expected. Despite the deviation from expectation between the calving interval EBV and its associated phenotype, a superior total merit index or a superior fertility EBV was indeed associated with an improvement in all detailed fertility performance phenotypes investigated. Results substantiate that breeding is a sustainable strategy of improving phenotypic performance in commercial dairy cattle and, by extension, profit.

Stakeholder perceptions of non-regulatory bovine health issues in Ireland: past and future perspectives

Background

In recent years, there have been multiple (political, environmental, cultural) drivers of change in Irish agriculture, including the establishment of Animal Health Ireland (AHI) in 2009, to provide leadership of non-regulatory livestock health issues (diseases and conditions of livestock that are endemic in Ireland but which are not currently subject to international legislation). In this study, we describe the opinion of stakeholders (farmers, veterinary practitioners and agricultural industry professional service providers), elicited by means of a survey, on their perceptions of changes in selected non-regulatory bovine health issues over the last 10 years and priority issues relevant to non-regulatory bovine health to be tackled over the next 10 years.

Results

A total of 673 individuals participated in the online questionnaire. For the majority of the non-regulatory bovine health issues, most participants felt there had been improvements over the last 10 years. However, professional service providers were generally more conservative in their response to improvements on-farm compared to farmers. Several issues, particularly BVD and udder health/milk quality, were viewed more positively by all relevant respondents. There was reasonable agreement between responses from different respondent types and sectors regarding the top three priorities relevant to non-regulatory bovine animal health for the next 10 years in Ireland, which included antimicrobial resistance (highlighting measures to reduce both on-farm usage and resistance), anthelmintic resistance, greenhouse emissions and calf welfare.

Conclusions

The results are encouraging, demonstrating a perception of improvement in a number of non-regulatory bovine health issues in Ireland over the last ten years. With respect to the next 10 years, stakeholders prioritised antimicrobial and anthelmintic resistance, greenhouse gas emissions and calf welfare, which aligns closely with broader societal concerns. This information is useful to AHI, particularly with respect to future priorities. However, these concerns are broad in scope and will require further considerations, including collaborations, between AHI and partnering organisations. Given that there were differences between farmers and professional service providers in responses, it is useful to consider how the aims and the benefits of future AHI programmes are framed and communicated to all stakeholders.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13620-020-00178-8.

Validation of a beef cattle maternal breeding objective based on a cross-sectional analysis of a large national cattle database

Despite the importance of validating any technology prior to recommendation for use, few studies exist in the scientific literature which have demonstrated the superior performance of high-ranking animals in a given total merit index; this is especially true for maternal cattle selection indexes. The objective of the present study was to demonstrate the impact of the Irish total merit maternal-based index and provide the benefits of using the Irish total merit maternal-based beef index as part of a breeding policy. The validation exercise was undertaken using 269,407 records (which included the cow's own records and her progeny records) from 92,300 females differing in a total merit index for maternal value; a comparison was also made with the Irish terminal index. Association analyses were undertaken within the framework of linear and threshold mixed models; the traits analyzed were fertility (e.g., calving interval), slaughter (e.g., harvest weight), live weight (e.g., weaning weight), and producer-recorded traits (e.g., docility). All traits were analyzed with the maternal index and terminal index fitted as covariate(s) separately. Depending on the independent variable analyzed, the other fixed effects included: parity of cow, heterosis and recombination loss of cow and/or progeny, gender of progeny, and the estimated breeding value of the sire; contemporary group was included as a random effect. The results demonstrate the effectiveness of using total merit indexes to improve performance in a whole range of different traits, despite the often antagonistic genetic correlations among traits that underpin the index. Cows excelling on the maternal index had less calving difficulty, superior fertility performance, lighter carcasses, and live weight, as well as being more easily managed. Additionally, progeny of higher maternal index cows were lighter at birth and more docile albeit with a small impact on slaughter traits. In contrast, higher terminal index cows had more calving difficulty, compromised fertility and had heavier carcasses themselves as well as their progeny. While the differences in phenotypic performance between groups on maternal index was, in most instances, relatively small, the benefits are: (1) expected to be greater when more genetically extreme groups of animals are evaluated and (2) expected to accumulate over time given the cumulative and permanent properties of breeding schemes.

Factors associated with the weight of individual primal cuts and their inter-relationship in cattle

Input parameters for decision support tools are comprised of, amongst others, knowledge of the associated factors and the extent of those associations with the animal-level feature of interest. The objective of the present study was to quantify the association between animal-level factors with primal cut yields in cattle and to understand the extent of the variability in primal cut yields independent carcass weight. The data used consisted of the weight of 14 primal carcass cuts (as well as carcass weight, conformation, and fat score) on up to 54,250 young cattle slaughtered between the years 2013 and 2017. Linear mixed models, with contemporary group of herd-sex-season of slaughter as a random effect, were used to quantify the associations between a range of model fixed effects with each primal cut separately. Fixed effects in the model were dam parity, heterosis coefficient, recombination loss, a covariate per breed representing the proportion of Angus, Belgian Blue, Charolais, Jersey, Hereford, Limousin, Simmental, and Holstein-Friesian and a three-way interaction between whether the animal was born in a dairy or beef herd, sex, and age at slaughter, with or without carcass weight as a covariate in the mixed model. The raw correlations among all cuts were all positive varying from 0.33 (between the bavette and the striploin) to 0.93 (between the topside and knuckle). The partial correlation among cuts, following adjustment for differences in carcass weight, varied from -0.36 to 0.74. Age at slaughter, sex, dam parity, and breed were all associated (P < 0.05) with the primal cut weight. Knowledge of the relationship between the individual primal cuts, and the solutions from the models developed in the study, could prove useful inputs for decision support systems to increase performance.

Genome-wide association study of endo-parasite phenotypes using imputed whole-genome sequence data in dairy and beef cattle

BACKGROUND

Quantitative genetic studies suggest the existence of variation at the genome level that affects the ability of cattle to resist to parasitic diseases. The objective of the current study was to identify regions of the bovine genome that are associated with resistance to endo-parasites.

METHODS

Individual cattle records were available for Fasciola hepatica-damaged liver from 18 abattoirs. Deregressed estimated breeding values (EBV) for F. hepatica-damaged liver were generated for genotyped animals with a record for F. hepatica-damaged liver and for genotyped sires with a least one progeny record for F. hepatica-damaged liver; 3702 animals were available. In addition, individual cow records for antibody response to F. hepatica on 6388 genotyped dairy cows, antibody response to Ostertagia ostertagi on 8334 genotyped dairy cows and antibody response to Neospora caninum on 4597 genotyped dairy cows were adjusted for non-genetic effects. Genotypes were imputed to whole-sequence; after edits, 14,190,141 single nucleotide polymorphisms (SNPs) and 16,603,644 SNPs were available for cattle with deregressed EBV for F. hepatica-damaged liver and cows with an antibody response to a parasitic disease, respectively. Association analyses were undertaken using linear regression on one SNP at a time, in which a genomic relationship matrix accounted for the relationships between animals.

RESULTS

Genomic regions for F. hepatica-damaged liver were located on Bos taurus autosomes (BTA) 1, 8, 11, 16, 17 and 18; each region included at least one SNP with a p value lower than 10^-6. Five SNPs were identified as significant (q value < 0.05) for antibody response to N. caninum and were located on BTA21 or 25. For antibody response to F. hepatica and O. ostertagi, six and nine quantitative trait loci (QTL) regions that included at least one SNP with a p value lower than 10^-6 were identified, respectively. Gene set enrichment analysis revealed a significant association between functional annotations related to the olfactory system and QTL that were suggestively associated with endo-parasite phenotypes.

CONCLUSIONS

A number of novel genomic regions were suggestively associated with endo-parasite phenotypes across the bovine genome and two genomic regions on BTA21 and 25 were associated with antibody response to N. caninum.

Genetic parameters for variability in the birth of persistently infected cattle following likely in utero exposure to bovine viral diarrhea virus

Genetic selection is an inexpensive and complementary strategy to traditional methods of improving animal health and welfare. Nonetheless, endeavors to incorporate animal health and welfare traits in international breeding programs have been hampered by the availability of informative phenotypes. The recent eradication program for bovine viral diarrhea (BVD) in the Republic of Ireland has provided an opportunity to quantify the potential benefits that genetic selection could offer BVD eradication programs elsewhere, as well as inform possible eradication programs for other diseases in the Republic of Ireland. Using a dataset of 188,085 Irish calves, the estimated direct and maternal heritability estimates for the birth of persistently infected calves following likely in utero exposure to BVD virus ranged from not different from zero (linear model) to 0.29 (SE = 0.075; threshold model) and from essentially zero (linear model) to 0.04 (SE = 0.033; threshold model), respectively. The corresponding genetic SD for the direct and maternal effect of the binary trait (0, 1) ranged from 0.005 (linear model) to 0.56 (threshold model) units and ranged from 0.00008 (linear model) to 0.20 (threshold model) units, respectively. The coefficient of direct genetic variation based on the linear model was 2.56% indicating considerable genetic variation could be exploited. Based on results from the linear model in the present study, there is the potential to reduce the incidence of persistent infection in cattle by on average 0.11 percentage units per year which is cumulative and permanent. Therefore, genetic selection can contribute to reducing the incidence of persistent infection in cattle. Moreover, where populations are free from persistent infection, inclusion of the estimated genetic merit for BVD in national breeding indexes could contribute to a preservation of a BVD-free status. Results from the present study can be used to inform breeding programs of the potential genetic gains achievable. Moreover, the approaches used in the present study can be applied to other diseases when data become available.

Genome-wide associations and functional gene analyses for endoparasite resistance in an endangered population of native German Black Pied cattle

Background

Gastrointestinal nematodes (GIN), liver flukes (Fasciola hepatica) and bovine lungworms (Dictyocaulus viviparus) are the most important parasitic agents in pastured dairy cattle. Endoparasite infections are associated with reduced milk production and detrimental impacts on female fertility, contributing to economic losses in affected farms. In quantitative-genetic studies, the heritabilities for GIN and F. hepatica were moderate, encouraging studies on genomic scales. Genome-wide association studies (GWAS) based on dense single nucleotide polymorphism (SNP) marker panels allow exploration of the underlying genomic architecture of complex disease traits. The current GWAS combined the identification of potential candidate genes with pathway analyses to obtain deeper insights into bovine immune response and the mechanisms of resistance against endoparasite infections.

Results

A 2-step approach was applied to infer genome-wide associations in an endangered dual-purpose cattle subpopulation [Deutsches Schwarzbuntes Niederungsrind (DSN)] with a limited number of phenotypic records. First, endoparasite traits from a population of 1166 Black and White dairy cows [including Holstein Friesian (HF) and DSN] naturally infected with GIN, F. hepatica and D. viviparus were precorrected for fixed effects using linear mixed models. Afterwards, the precorrected phenotypes were the dependent traits (rFEC-GIN, rFEC-FH, and rFLC-DV) in GWAS based on 423,654 SNPs from 148 DSN cows. We identified 44 SNPs above the genome-wide significance threshold (p_Bonf = 4.47 × 10^− 7), and 145 associations surpassed the chromosome-wide significance threshold (range: 7.47 × 10^− 6 on BTA 1 to 2.18 × 10^− 5 on BTA 28). The associated SNPs identified were annotated to 23 candidate genes. The DAVID analysis inferred four pathways as being related to immune response mechanisms or involved in host-parasite interactions. SNP effect correlations considering specific chromosome segments indicate that breeding for resistance to GIN or F. hepatica as measured by fecal egg counts is genetically associated with a higher risk for udder infections.

Conclusions

We detected a large number of loci with small to moderate effects for endoparasite resistance. The potential candidate genes regulating resistance identified were pathogen-specific. Genetic antagonistic associations between disease resistance and productivity were specific for specific chromosome segments. The 2-step approach was a valid methodological approach to infer genetic mechanisms in an endangered breed with a limited number of phenotypic records.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5659-4) contains supplementary material, which is available to authorized users.

Candidate genes associated with the heritable humoral response to Mycobacterium avium ssp. paratuberculosis in dairy cows have factors in common with gastrointestinal diseases in humans

Infection of cattle with bovine paratuberculosis (i.e., Johne's disease) is caused by Mycobacterium avium ssp. paratuberculosis (MAP) and results in a chronic incurable gastroenteritis. This disease, which has economic ramifications for the cattle industry, is increasing in detected prevalence globally; subclinically infected animals can silently shed the bacterium into the environment for years, exposing contemporaries and hampering disease-control programs. The objective of the present study was to first quantify the genetic parameters for humoral response to MAP in dairy cattle. This was followed by a genome-based association analysis and subsequent downstream bioinformatic analyses from imputed whole genome sequence SNP data. After edits, ELISA test records were available on 136,767 cows; analyses were also undertaken on a subset of 33,818 of these animals from herds with at least 5 MAP ELISA-positive cows, with at least 1 of those positive cows being homebred. Variance components were estimated using univariate animal and sire linear mixed models. The heritability calculated from the animal model for humoral response to MAP using alternative phenotype definitions varied from 0.02 (standard error = 0.003) to 0.05 (standard error = 0.008). The genome-based associations were undertaken within a mixed model framework using weighted deregressed estimated breeding values as a dependent variable on 1,883 phenotyped animals that were ≥87.5% Holstein-Friesian. Putative susceptibility quantitative trait loci (QTL) were identified on Bos taurus autosome 1, 3, 5, 6, 8, 9, 10, 11, 13, 14, 18, 21, 23, 25, 26, 27, and 29; mapping the most significant SNP to genes within and overlapping these QTL revealed that the most significant associations were with the 10 functional candidate genes KALRN, ZBTB20, LPP, SLA2, FI3A1, LRCH3, DNAJC6, ZDHHC14, SNX1, and HAS2. Pathway analysis failed to reveal significantly enriched biological pathways, when both bovine-specific pathway data and human ortholog data were taken into account. The existence of genetic variation for MAP susceptibility in a large data set of dairy cows signifies the potential of breeding programs for reducing MAP susceptibility. Furthermore, the identification of susceptible QTL facilitates greater biological understanding of bovine paratuberculosis and potential therapeutic targets for future investigation. The novel molecular similarities identified between bovine paratuberculosis and human inflammatory bowel disease suggest potential for human therapeutic interventions to be translated to veterinary medicine and vice versa.

Little genetic variability in resilience among cattle exists for a range of performance traits across herds in Ireland differing in Fasciola hepatica prevalence

It is anticipated that in the future, livestock will be exposed to a greater risk of infection from parasitic diseases. Therefore, future breeding strategies for livestock, which are generally long-term strategies for change, should target animals adaptable to environments with a high parasitic load. Covariance components were estimated in the present study for a selection of dairy and beef performance traits over herd-years differing in Fasciola hepatica load using random regression sire models. Herd-year prevalence of F. hepatica was determined by using F. hepatica-damaged liver phenotypes which were recorded in abattoirs nationally. The data analyzed consisted up to 83,821 lactation records from dairy cows for a range of milk production and fertility traits, as well as 105,054 young animals with carcass-related information obtained at slaughter. Reaction norms for individual sires were derived from the random regression coefficients. The heritability and additive genetic standard deviations for all traits analyzed remained relatively constant as herd-year F. hepatica prevalence gradient increased up to a prevalence level of 0.7; although there was a large increase in heritability and additive genetic standard deviation for milk and fertility traits in the observed F. hepatica prevalence levels >0.7, only 5% of the data existed in herd-year prevalence levels >0.7. Very little rescaling, therefore, exists across differing herd-year F. hepatica prevalence levels. Within-trait genetic correlations among the performance traits across different herd-year F. hepatica prevalence levels were less than unity for all traits. Nevertheless, within-trait genetic correlations for milk production and carcass traits were all >0.8 for F. hepatica prevalence levels between 0.2 and 0.8. The lowest estimate of within-trait genetic correlations for the different fertility traits ranged from -0.03 (SE = 1.09) in age of first calving to 0.54 (SE = 0.22) for calving to first service interval. Therefore, there was reranking of sires for fertility traits across different F. hepatica prevalence levels. In conclusion, there was little or no genetic variability in sensitivity to F. hepatica prevalence levels among cattle for milk production and carcass traits. But, some genetic variability in sensitivity among dairy cows did exist for fertility traits measured across herds differing in F. hepatica prevalence.

Genetic correlations between endo-parasite phenotypes and economically important traits in dairy and beef cattle

Parasitic diseases have economic consequences in cattle production systems. Although breeding for parasite resistance can complement current control practices to reduce the prevalence globally, there is little knowledge of the implications of such a strategy on other performance traits. Records on individual animal antibody responses to Fasciola hepatica, Ostertagia ostertagi, and Neospora caninum were available from cows in 68 dairy herds (study herds); national abattoir data on F. hepatica-damaged livers were also available from dairy and beef cattle. After data edits, 9,271 dairy cows remained in the study herd dataset, whereas 19,542 dairy cows and 68,048 young dairy and beef animals had a record for the presence or absence of F. hepatica-damaged liver in the national dataset. Milk, reproductive, and carcass phenotypes were also available for a proportion of these animals as well as their contemporaries. Linear mixed models were used to estimate variance components of antibody responses to the three parasites; covariance components were estimated between the parasite phenotypes and economically important traits. Heritability of antibody responses to the different parasites, when treated as a continuous trait, ranged from 0.07 (O. ostertagi) to 0.13 (F. hepatica), whereas the coefficient of genetic variation ranged from 4% (O. ostertagi) to 20% (F. hepatica). The antibody response to N. caninum was genetically correlated with the antibody response to both F. hepatica (-0.29) and O. ostertagi (-0.67); a moderately positive genetic correlation existed between the antibody response to F. hepatica and O. ostertagi (0.66). Genetic correlations between the parasite phenotypes and the milk production traits were all close to zero (-0.14 to 0.10), as were the genetic correlations between F. hepatica-damaged livers and the carcass traits of carcass weight, conformation, and fat score evaluated in cows and young animals (0.00 to 0.16). The genetic correlation between F. hepatica-damaged livers in cows and milk somatic cell score was 0.32 (SE = 0.20). Antibody responses to F. hepatica and O. ostertagi had favorable genetic correlations with fertility traits, but conversely, antibody response to N. caninum and F. hepatica-damaged livers were unfavorably genetically correlated with fertility. This study provides the necessary information to undertake national multitrait genetic evaluations for parasite phenotypes.

Genetic line comparisons and genetic parameters for endoparasite infections and test-day milk production traits

Keeping dairy cows in grassland systems relies on detailed analyses of genetic resistance against endoparasite infections, including between- and within-breed genetic evaluations. The objectives of this study were (1) to compare different Black and White dairy cattle selection lines for endoparasite infections and (2) the estimation of genetic (co)variance components for endoparasite and test-day milk production traits within the Black and White cattle population. A total of 2,006 fecal samples were taken during 2 farm visits in summer and autumn 2015 from 1,166 cows kept in 17 small- and medium-scale organic and conventional German grassland farms. Fecal egg counts were determined for gastrointestinal nematodes (FEC-GIN) and flukes (FEC-FLU), and fecal larvae counts for the bovine lungworm Dictyocaulus viviparus (FLC-DV). The lowest values for gastrointestinal nematode infections were identified for genetic lines adopted to pasture-based production systems, especially selection lines from New Zealand. Heritabilities were low for FEC-GIN (0.05-0.06 ± 0.04) and FLC-DV (0.05 ± 0.04), but moderate for FEC-FLU (0.33 ± 0.06). Almost identical heritabilities were estimated for different endoparasite trait transformations (log-transformation, square root). The genetic correlation between FEC-GIN and FLC-DV was 1.00 ± 0.60, slightly negative between FEC-GIN and FEC-FLU (-0.10 ± 0.27), and close to zero between FLC-DV and FEC-FLU (0.03 ± 0.30). Random regression test-day models on a continuous time scale [days in milk (DIM)] were applied to estimate genetic relationships between endoparasite and longitudinal test-day production traits. Genetic correlations were negative between FEC-GIN and milk yield (MY) until DIM 85, and between FEC-FLU and MY until DIM 215. Genetic correlations between FLC-DV and MY were negative throughout lactation, indicating improved disease resistance for high-productivity cows. Genetic relationships between FEC-GIN and FEC-FLU with milk protein content were negative for all DIM. Apart from the very early and very late lactation stage, genetic correlations between FEC-GIN and milk fat content were negative, whereas they were positive for FEC-FLU. Genetic correlations between FEC-GIN and somatic cell score were positive, indicating similar genetic mechanisms for susceptibility to udder and endoparasite infections. The moderate heritabilities for FEC-FLU suggest inclusion of FEC-FLU into overall organic dairy cattle breeding goals to achieve long-term selection response for disease resistance.