Genetic evaluation of mastitis liability and recovery through longitudinal analysis of transition probabilities

Review: Livestock disease resilience: from individual to herd level

Infectious diseases are a major threat to the sustainable production of high-producing animals. Control efforts, such as vaccination or breeding approaches often target improvements to individual resilience to infections, i.e., they strengthen an animal's ability to cope with infection, rather than preventing infection per se. There is increasing evidence for the contribution of non-clinical carriers (animals that become infected and are infectious but do not develop clinical signs) to the overall health and production of livestock populations for a wide range of infectious diseases. Therefore, we strongly advocate a shift of focus from increasing the disease resilience of individual animals to herd disease resilience as the appropriate target for sustainable disease control in livestock. Herd disease resilience not only captures the direct effects of vaccination or host genetics on the health and production performance of individuals but also the indirect effects on the environmental pathogen load that herd members are exposed to. For diseases primarily caused by infectious pathogens shed by herd members, these indirect effects on herd resilience are mediated both by individual susceptibility to infection and by characteristics (magnitude of infectiousness, duration of infectious period) that influence pathogen shedding from infected individuals. We review what is currently known about how vaccination and selective breeding affect herd disease resilience and its underlying components, and outline the changes required for improvement. To this purpose, we also seek to clarify and harmonise the terminology used in the different animal science disciplines to facilitate future collaborative approaches to infectious disease control in livestock.

Estimating individuals' genetic and non-genetic effects underlying infectious disease transmission from temporal epidemic data

Individuals differ widely in their contribution to the spread of infection within and across populations. Three key epidemiological host traits affect infectious disease spread: susceptibility (propensity to acquire infection), infectivity (propensity to transmit infection to others) and recoverability (propensity to recover quickly). Interventions aiming to reduce disease spread may target improvement in any one of these traits, but the necessary statistical methods for obtaining risk estimates are lacking. In this paper we introduce a novel software tool called SIRE (standing for "Susceptibility, Infectivity and Recoverability Estimation"), which allows for the first time simultaneous estimation of the genetic effect of a single nucleotide polymorphism (SNP), as well as non-genetic influences on these three unobservable host traits. SIRE implements a flexible Bayesian algorithm which accommodates a wide range of disease surveillance data comprising any combination of recorded individual infection and/or recovery times, or disease diagnostic test results. Different genetic and non-genetic regulations and data scenarios (representing realistic recording schemes) were simulated to validate SIRE and to assess their impact on the precision, accuracy and bias of parameter estimates. This analysis revealed that with few exceptions, SIRE provides unbiased, accurate parameter estimates associated with all three host traits. For most scenarios, SNP effects associated with recoverability can be estimated with highest precision, followed by susceptibility. For infectivity, many epidemics with few individuals give substantially more statistical power to identify SNP effects than the reverse. Importantly, precise estimates of SNP and other effects could be obtained even in the case of incomplete, censored and relatively infrequent measurements of individuals' infection or survival status, albeit requiring more individuals to yield equivalent precision. SIRE represents a new tool for analysing a wide range of experimental and field disease data with the aim of discovering and validating SNPs and other factors controlling infectious disease transmission.

A genetic study of new udder health indicator traits with data from automatic milking systems

The current study aimed to investigate new udder health traits based on data from automatic milking systems (AMS) for use in routine genetic evaluations. Data were from 77 commercial herds; out of these, 24 had equipment for measuring online cell count (OCC), whereas all had data on electrical conductivity (EC). A total of 4,714 Norwegian Red dairy cows and 2,363,928 milkings were included in the genetic analyses. Electrical conductivity was available on quarter level for each milking, whereas OCC was measured per milking. The AMS traits analyzed were log-transformed online cell count (lnOCC), maximum conductivity (ECmax), mean conductivity (ECmean), elevated mastitis risk (EMR), and log-transformed EMR (lnEMR). In addition, lactation mean somatic cell score (LSCS) was collected from the Norwegian dairy herd recording system. Elevated mastitis risk expresses the probability of a cow having mastitis and was calculated from smoothed lnOCC values according to individual trend and level of the OCC curve. The udder health traits from AMS were analyzed as repeated milkings from 30 to 320 DIM, and LSCS as repeated parities. In addition, both ECmax and lnOCC were analyzed as multiple traits by splitting the lactation into 5 periods. (Co)variance components were estimated from bivariate mixed linear animal models, and investigated traits showed genetic variation. Estimated heritabilities of ECmean, ECmax, and lnEMR were 0.35, 0.23, and 0.12, respectively, whereas EMR and lnOCC both showed heritabilities of 0.09. Heritability varied between periods of lactation, from 0.04 to 0.13 for lnOCC and from 0.12 to 0.27 for ECmax, although standard errors of certain periods were large. Genetic correlations among the AMS traits ranged from 0 to 0.99. The genetic correlations between EC-based traits and OCC-based traits in AMS were 0. Genetic correlations with LSCS were favorable, ranging from 0.37 to 0.80 (±0.11-0.22). The strongest correlation (0.80 ± 0.13) was found between LSCS and lnEMR. Results question the value of ECmax and ECmean as indicators of udder health in genetic evaluations and suggest OCC to be more valuable in this manner. This study demonstrates a potential of using AMS data as additional information on udder health for genetic evaluations, although further investigation is recommended before these traits can be implemented.

Genome-Wide Association Study for Susceptibility to and Recoverability From Mastitis in Danish Holstein Cows

Because mastitis is very frequent and unavoidable, adding recovery information into the analysis for genetic evaluation of mastitis is of great interest from economical and animal welfare point of view. Here we have performed genome-wide association studies (GWAS) to identify associated single nucleotide polymorphisms (SNPs) and investigate the genetic background not only for susceptibility to – but also for recoverability from mastitis. Somatic cell count records from 993 Danish Holstein cows genotyped for a total of 39378 autosomal SNP markers were used for the association analysis. Single SNP regression analysis was performed using the statistical software package DMU. Substitution effect of each SNP was tested with a t-test and a genome-wide significance level of P-value < 10^-4 was used to declare significant SNP-trait association. A number of significant SNP variants were identified for both traits. Many of the SNP variants associated either with susceptibility to – or recoverability from mastitis were located in or very near to genes that have been reported for their role in the immune system. Genes involved in lymphocyte developments (e.g., MAST3 and STAB2) and genes involved in macrophage recruitment and regulation of inflammations (PDGFD and PTX3) were suggested as possible causal genes for susceptibility to – and recoverability from mastitis, respectively. However, this is the first GWAS study for recoverability from mastitis and our results need to be validated. The findings in the current study are, therefore, a starting point for further investigations in identifying causal genetic variants or chromosomal regions for both susceptibility to – and recoverability from mastitis.

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.

Integration of epidemiology into the genetic analysis of mastitis in Swedish Holstein

Heritability of mastitis (and diseases in general) tends to be low. One possible cause is that no clear distinction can be made between resistant and nonresistant animals, because healthy animals include animals that have not been exposed to pathogens and resistant animals. To account for this, we quantified the prevalence of clinical mastitis (CM) and subclinical mastitis (SCM) in 2,069 Swedish Holstein herds as a measure of exposure. Herd prevalence averaged 26.5% for SCM and 6.4% for CM; 61% of the first lactations of 177,309 cows were classified as having at least one case of SCM and 10% as having CM. In a reaction norm approach, heritability of (S)CM was quantified as a function of herd prevalence of (S)CM. The best-fitting model was a second-order polynomial of first-lactation cow SCM as a function of herd prevalence SCM, and a first-order (linear) polynomial of first-lactation cow CM as a function of CM herd prevalence. Heritability for SCM ranged from 0.069 to 0.105 and for CM from 0.016 to 0.032. For both, we found no clear effect of herd prevalence on their heritability. Genetic correlations within traits across herd prevalences were all greater than 0.92. Whether relationships among prevalence, exposure, disease, and genetics were as expected is a matter of discussion, but reaction norm analyses may be a valuable tool for epidemiological genetics.