In dairy cattle, resistance, tolerance and resilience refer to the adaptation ability to a broad range of environmental conditions, implying stable performances (e.g. production level, fertility status) independent from disease or infection pressure. All three mechanisms resistance, tolerance and resilience contribute to overall robustness, implying the evaluation of phenotyping and breeding strategies for improved robustness in dairy cattle populations. Classically, breeding approaches on improved robustness rely on simple production traits, in combination with detailed environmental descriptors and enhanced statistical modelling to infer possible genotype by environment interactions. In this regard, innovative environmental descriptors were heat stress indicators, and statistical modelling focussed on random regression or reaction norm methodology. A robust animal has high breeding values over a broad spectra of environmental levels. During the last years, direct health traits were included into selection indices, implying advances in genetic evaluations for traits being linked to resistance or tolerance against infectious and non-infectious diseases. Up to now, genetic evaluation for health traits is primarily based on subjectively measured producer-recorded data, with disease trait heritabilities in a low-to-moderate range. Thus, it is imperative to identify objectively measurable phenotypes as suitable biomarkers. New technologies (e.g. mid-infrared spectrometry) offer possibilities to determine potential biomarkers via laboratory analyses. Novel biomarkers include measurable physiological traits (e.g. serum metabolites, hormone levels) as indicators for a current infection, or the host’s reaction to environmental stressors. The rumen microbiome composition is proposed as a biomarker to detect interactions between host genotype and environmental effects. The understanding of host genetic variation in disease resistance and individual expression of robustness encourages analyses on the underlying immune response (IR) system. Recent advances have been made in order to infer the genetic background of IR traits and cows immunological competence in relation to functional and production traits. Thus, a last aspect of this review addresses the genetic background and current state of genetic control for resistance to economically relevant infectious and non-infectious dairy cattle diseases by considering immune-related factors.

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.

To narrow the gap of missing knowledge on Rickettsia spp. and Anaplasma phagocytophilum infections in ticks in northwestern Germany and, at the same time, to provide first prevalence data on these pathogens in the city of Hamburg, a total of 1,400 questing Ixodes ricinus ticks were collected at ten different public green areas from April until October 2011. Ticks were examined using probe-based quantitative real-time PCR. A percentage of 3.6% (51/1,400) ticks were tested positive for A. phagocytophilum infections divided into 2.1% (3/141) adults [1.7% (1/60) females and 2.5% (2/81) males] and 3.8% (48/1,259) nymphs. The percentage of infected ticks per sampling site varied statistically significantly from 0.7% (1/140) to 12.1% (17/140), whereas between sampling months, no statistically significant differences were observed (2.0-6.5%, 4-13/140). The overall Rickettsia spp. infection rate was 52.5% (735/1,400). In adult ticks, Rickettsia spp. infection rate was 56% (79/141) divided into 61.7% (37/60) infected females and 51.9% (42/81) infected males. Nymphs showed an infection rate of 52.1% (656/1,259). In contrast to A. phagocytophilum infections, no statistically significant differences in Rickettsia spp. infection rates among sampling sites (44.3-63.6%, 62-89/140) were observed, whereas seasonal variations were obvious: the percentage of Rickettsia-positive ticks was significantly lower in April (36.5%, 73/200) and May (29.5%, 59/200) compared to the summer and fall months (55.0-64.5%, 110-129/200). Rickettsia species differentiation via real-time pyrosequencing revealed Rickettsia helvetica as the only occurring species. Co-infections with both Rickettsia spp. and A. phagocytophilum were detected in 2.0% (28/1,400) of the ticks. The present study revealed that in the city of Hamburg, the tick infection rate with A. phagocytophilum is comparable with other German data, whereas the Rickettsia spp. infection rate of 52.5% is by far the highest prevalence detected in Germany so far. As the city of Hamburg has 1.8 million inhabitants and attracts millions of tourists every year, the potential health risk should not be underestimated.

To obtain initial data on Borrelia burgdorferi sensu lato (Spirochaetales: Spirochaetaceae) in Ixodes ricinus (Ixodida: Ixodidae) ticks in Hamburg, Germany, 1400 questing ticks were collected by flagging at 10 different public recreation areas in 2011 and analysed using probe-based quantitative real-time polymerase chain reaction. The overall rate of infection with B. burgdorferi s.l. was 34.1%; 30.0% of adults were infected (36.7% of females and 26.0% of males), as were 34.5% of nymphs. Significant differences in tick infection rates were observed between the spring and summer/autumn months, as well as among sampling locations. Borrelia genospecies identification by reverse line blotting was successful in 43.6% of positive tick samples. The most frequent genospecies was Borrelia garinii/Borrelia bavariensis, followed by Borrelia afzelii, Borrelia valaisiana, B. burgdorferi sensu stricto, Borrelia spielmanii, Borrelia bissettii and Borrelia lusitaniae. Based on previously published data, co-infection of Borrelia and Rickettsiales spp. was determined in 25.8% of ticks. Overall, 22.9% of ticks were co-infected with Rickettsia spp. (Rickettsiales: Rickettsiaceae), 1.7% with Anaplasma phagocytophilum (Rickettsiales: Anaplasmataceae), and 1.2% with both pathogens. Study results show a high prevalence of Borrelia-positive ticks in recreation areas in the northern German city of Hamburg and the potential health risk to humans in these areas should not be underestimated.

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...

German Black Pied cattle (DSN) is an endangered population of about 2,550 dual-purpose cattle in Germany. Having a milk yield of about 2,500 kg less than the predominant dairy breed Holstein, the preservation of DSN is supported by the German government and the EU. The identification of the genomic loci affecting milk production in DSN can provide a basis for selection decisions for genetic improvement of DSN in order to increase market chances through the improvement of milk yield. A genome-wide association analysis of 30 milk traits was conducted in different lactation periods and numbers. Association using multiple linear regression models in R was performed on 1,490 DSN cattle genotyped with BovineSNP50 SNP-chip. 41 significant and 20 suggestive SNPs affecting milk production traits in DSN were identified, as well as 15 additional SNPs for protein content which are less reliable due to high inflation. The most significant effects on milk yield in DSN were detected on chromosomes 1, 6, and 20. The region on chromosome 6 was located nearby the casein gene cluster and the corresponding haplotype overlapped the CSN3 gene (casein kappa). Associations for fat and protein yield and content were also detected. High correlation between traits of the same lactation period or number led to some SNPs being significant for multiple investigated traits. Half of all identified SNPs have been reported in other studies, previously. 15 SNPs were associated with the same traits in other breeds. The other associated SNPs have been reported previously for traits such as exterior, health, meat and carcass, production, and reproduction traits. No association could be detected between DGAT1 and other known milk genes with milk production traits despite the close relationship between DSN and Holstein. The results of this study confirmed that many SNPs identified in other breeds as associated with milk traits also affect milk traits in dual-purpose DSN cattle and can be used for further genetic analysis to identify genes and causal variants that affect milk production in DSN cattle.

No abstract

SummaryThe aim of this study was to detect selection signatures considering cows from the German Holstein (GH) and the local dual‐purpose black and white (DSN) population, as well as from generated sub‐populations. The 4654 GH and 261 DSN cows were genotyped with the BovineSNP50 Genotyping BeadChip. The geographical herd location was used as an environmental descriptor to create the East‐DSN and West‐DSN sub‐populations. In addition, two further sub‐populations of GH cows were generated, using the extreme values for solutions of residual effects of cows for the claw disorder dermatitis digitalis. These groups represented the most susceptible and most resistant cows. We used cross‐population extended haplotype homozygosity methodology (XP‐EHH) to identify the most recent selection signatures. Furthermore, we calculated Wright’s fixation index (FST). Chromosomal segments for the top 0.1 percentile of negative or positive XP‐EHH scores were studied in detail. For gene annotations, we used the Ensembl database and we considered a window of 250 kbp downstream and upstream of each core SNP corresponding to peaks of XP‐EHH. In addition, functional interactions among potential candidate genes were inferred via gene network analyses. The most outstanding XP‐EHH score was on chromosome 12 (at 77.34 Mb) for DSN and on chromosome 20 (at 36.29–38.42 Mb) for GH. Selection signature locations harbored QTL for several economically important milk and meat quality traits, reflecting the different breeding goals for GH and DSN. The average FST value between GH and DSN was quite low (0.068), indicating shared founders. For group stratifications according to cow health, several identified potential candidate genes influence disease resistance, especially to dermatitis digitalis.