Genetic parameters and genomic breeding values for digital dermatitis in Holstein Friesian dairy cattle: host susceptibility, infectivity and the basic reproduction ratio

International audienceBackground : For infectious diseases, the probability that an animal gets infected depends on its own susceptibility, and on the number of infectious herd mates and their infectivity. Together with the duration of the infectious period, susceptibility and infectivity determine the basic reproduction ratio of the disease ( R0 ). R0 is the average number of secondary cases caused by a typical infectious individual in an otherwise uninfected population. An infectious disease dies out when R0 is lower than 1. Thus, breeding strategies that aim at reducing disease prevalence should focus on reducing R0 , preferably to a value lower than 1. In animal breeding, however, R0 has received little attention. Here, we estimate the additive genetic variance in host susceptibility, host infectivity, and R0 for the endemic claw disease digital dermatitis (DD) in Holstein Friesian dairy cattle, and estimate genomic breeding values (GEBV ) for these traits. We recorded DD disease status of both hind claws of 1513 cows from 12 Dutch dairy farms, every 2 weeks, 11 times. The genotype data consisted of 75,904 single nucleotide polymorphisms (SNPs) for 1401 of the cows. We modelled the probability that a cow got infected between recordings, and compared four generalized linear mixed models. All models included a genetic effect for susceptibility; Models 2 and 4 also included a genetic effect for infectivity, while Models 1 and 2 included a farm*period interaction. We corrected for variation in exposure to infectious herd mates via an offset.Results : GEBV for R0 from the model that included genetic effects for susceptibility only had an accuracy of ~ 0.39 based on cross-validation between farms, which is very high given the limited amount of data and the complexity of the trait. Models with a genetic effect for infectivity showed a larger bias, but also a slightly higher accuracy of GEBV. Additive genetic standard deviation for R0 was large, i.e. ~ 1.17, while the mean R0 was 2.36.Conclusions : GEBV for R0 showed substantial variation. The mean R0 was only about one genetic standard deviation greater than 1. These results suggest that lowering DD prevalence by selective breeding is promising

A model to estimate effects of SNPs on host susceptibility and infectivity for an endemic infectious disease

International audienceAbstractBackgroundInfectious diseases in farm animals affect animal health, decrease animal welfare and can affect human health. Selection and breeding of host individuals with desirable traits regarding infectious diseases can help to fight disease transmission, which is affected by two types of (genetic) traits: host susceptibility and host infectivity. Quantitative genetic studies on infectious diseases generally connect an individual’s disease status to its own genotype, and therefore capture genetic effects on susceptibility only. However, they usually ignore variation in exposure to infectious herd mates, which may limit the accuracy of estimates of genetic effects on susceptibility. Moreover, genetic effects on infectivity will exist as well. Thus, to design optimal breeding strategies, it is essential that genetic effects on infectivity are quantified. Given the potential importance of genetic effects on infectivity, we set out to develop a model to estimate the effect of single nucleotide polymorphisms (SNPs) on both host susceptibility and host infectivity. To evaluate the quality of the resulting SNP effect estimates, we simulated an endemic disease in 10 groups of 100 individuals, and recorded time-series data on individual disease status. We quantified bias and precision of the estimates for different sizes of SNP effects, and identified the optimum recording interval when the number of records is limited.ResultsWe present a generalized linear mixed model to estimate the effect of SNPs on both host susceptibility and host infectivity. SNP effects were on average slightly underestimated, i.e. estimates were conservative. Estimates were less precise for infectivity than for susceptibility. Given our sample size, the power to estimate SNP effects for susceptibility was 100% for differences between genotypes of a factor 1.56 or more, and was higher than 60% for infectivity for differences between genotypes of a factor 4 or more. When disease status was recorded 11 times on each animal, the optimal recording interval was 25 to 50% of the average infectious period.ConclusionsOur model was able to estimate genetic effects on susceptibility and infectivity. In future genome-wide association studies, it may serve as a starting point to identify genes that affect disease transmission and disease prevalence

Quantifying transmission of Mycobacterium avium subsp. paratuberculosis among group-housed dairy calves

International audienceAbstractJohne’s disease (JD) is a chronic enteritis caused by Mycobacterium avium subsp. paratuberculosis (MAP), with control primarily aimed at preventing new infections among calves. The aim of the current study was to quantify calf-to-calf transmission of MAP among penmates in an experimental trial. Newborn Holstein bull calves (n = 32) were allocated into pens of 4, with 2 inoculated (IN) calves and 2 calves that were contact exposed (CE). Calves were group-housed for 3 months, with frequent collection of fecal and blood samples and tissue collection after euthanasia. The basic reproduction ratio (R0) was estimated using a final size (FS) model with a susceptible-infected model, based on INF-γ ELISA and tissue culture followed by qPCR. In addition, the transmission rate parameter (β) for new shedding events was estimated using a general linearized method (GLM) model with a susceptible-infected-susceptible model based on culture, followed by qPCR, of fecal samples collected during group housing. The R0 was derived for IN and CE calves separately, due to a difference in susceptibility, as well as differences in duration of shedding events. Based on the FS model, interferon-γ results from blood samples resulted in a R0IG of 0.90 (0.24, 2.59) and tissue culture resulted in a R0T of 1.36 (0.45, 3.94). Based on the GLM model, the R0 for CE calves to begin shedding (R0CE) was 3.24 (1.14, 7.41). We concluded that transmission of MAP infection between penmates occurred and that transmission among calves may be an important cause of persistent MAP infection on dairy farms that is currently uncontrolled for in current JD control programs

Habitat loss exacerbates pathogen spread : An Agent-based model of avian influenza infection in migratory waterfowl

Habitat availability determines the distribution of migratory waterfowl along their flyway, which further influences the transmission and spatial spread of avian influenza viruses (ATVs). The extensive habitat loss in the East Asian-Australasian Flyway (EAAF) may have potentially altered the virus spread and transmission, but those consequences are rarely studied. We constructed 6 fall migration networks that differed in their level of habitat loss, wherein an increase in habitat loss resulted in smaller networks with fewer sites and links. We integrated an agent-based model and a susceptible-infected-recovered model to simulate waterfowl migration and AIV transmission. We found that extensive habitat loss in the EAAF can 1) relocate the outbreaks northwards, responding to the distribution changes of wintering waterfowl geese, 2) increase the outbreak risk in remaining sites due to larger goose congregations, and 3) facilitate AIV transmission in the migratory population. In addition, our modeling output was in line with the predictions from the concept of "migratory escape", i.e., the migration allows the geese to "escape" from the location where infection risk is high, affecting the pattern of infection prevalence in the waterfowl population. Our modeling shed light on the potential consequences of habitat loss in spreading and transmitting AIV at the flyway scale and suggested the driving mechanisms behind these effects, indicating the importance of conservation in changing spatial and temporal patterns of AIV outbreaks.Peer reviewe

Contribution of cats and dogs to SARS-CoV-2 transmission in households

INTRODUCTION: SARS-CoV-2 is known to jump across species. The occurrence of transmission in households between humans and companion animals has been shown, but the contribution of companion animals to the overall transmission within a household is unknown. The basic reproduction number ( R 0) is an important indicator to quantify transmission. For a pathogen with multiple host species, such as SARS-CoV-2, the basic reproduction number needs to be calculated from the partial reproduction numbers for each combination of host species. METHOD: In this study, the basic and partial reproduction numbers for SARS-CoV-2 were estimated by reanalyzing a survey of Dutch households with dogs and cats and minimally one SARS-CoV-2-infected human. RESULTS: For households with cats, a clear correlation between the number of cats and the basic reproduction number (Spearman's correlation: p 0.40, p-value: 1.4 × 10 -5) was identified, while for dogs, the correlation was smaller and not significant (Spearman's correlation: p 0.12, p-value: 0.21). Partial reproduction numbers from cats or dogs to humans were 0.3 (0.0-2.0) and 0.3 (0.0-3.5) and from humans to cats or dogs were 0.6 (0.4-0.8) and 0.6 (0.4-0.9). DISCUSSION: Thus, the estimations of within-household transmission indicated the likelihood of transmission from these companion animals to humans and vice versa, but the observational nature of this study limited the ability to establish conclusive evidence. This study's findings support the advice provided during the pandemic to COVID-19 patients to maintain distance from companion animals as a precautionary measure and given the possibility of transmission, although there is an overall relatively limited impact on the pandemic when compared to human-to-human transmission

Quantifying Rift Valley fever virus transmission efficiency in a lamb-mosquito-lamb model

Rift Valley fever virus (RVFV) is a (re)emerging mosquito-borne pathogen impacting human and animal health. How RVFV spreads through a population depends on population-level and individual-level interactions between vector, host and pathogen. Here, we estimated the probability for RVFV to transmit to naive animals by experimentally exposing lambs to a bite of an infectious mosquito, and assessed if and how RVFV infection subsequently developed in the exposed animal. Aedes aegypti mosquitoes, previously infected via feeding on a viremic lamb, were used to expose naive lambs to the virus. Aedes aegypti colony mosquitoes were used as they are easy to maintain and readily feed in captivity. Other mosquito spp. could be examined with similar methodology. Lambs were exposed to either 1-3 (low exposure) or 7-9 (high exposure) infectious mosquitoes. All lambs in the high exposure group became viremic and showed characteristic signs of Rift Valley fever within 2-4 days post exposure. In contrast, 3 out of 12 lambs in the low exposure group developed viremia and disease, with similar peak-levels of viremia as the high exposure group but with some heterogeneity in the onset of viremia. These results suggest that the likelihood for successful infection of a ruminant host is affected by the number of infectious mosquitoes biting, but also highlights that a single bite of an infectious mosquito can result in disease. The per bite mosquito-to-host transmission efficiency was estimated at 28% (95% confidence interval: 15 - 47%). We subsequently combined this transmission efficiency with estimates for life traits of Aedes aegypti or related mosquitoes into a Ross-McDonald mathematical model to illustrate scenarios under which major RVFV outbreaks could occur in naïve populations (i.e., R0 >1). The model revealed that relatively high vector-to-host ratios as well as mosquitoes feeding preferably on competent hosts are required for R0 to exceed 1. Altogether, this study highlights the importance of experiments that mimic natural exposure to RVFV. The experiments facilitate a better understanding of the natural progression of disease and a direct way to obtain epidemiological parameters for mathematical models

Inferring bovine tuberculosis transmission between cattle and badgers via the environment and risk mapping

Bovine tuberculosis (bTB), caused by Mycobacterium bovis, is one of the most challenging and persistent health issues in many countries worldwide. In several countries, bTB control is complicated due to the presence of wildlife reservoirs of infection, i.e. European badger (Meles meles) in Ireland and the UK, which can transmit infection to cattle. However, a quantitative understanding of the role of cattle and badgers in bTB transmission is elusive, especially where there is spatial variation in relative density between badgers and cattle. Moreover, as these two species have infrequent direct contact, environmental transmission is likely to play a role, but the quantitative importance of the environment has not been assessed. Therefore, the objective of this study is to better understand bTB transmission between cattle and badgers via the environment in a spatially explicit context and to identify high-risk areas. We developed an environmental transmission model that incorporates both within-herd/territory transmission and between-species transmission, with the latter facilitated by badger territories overlapping with herd areas. Model parameters such as transmission rate parameters and the decay rate parameter of M. bovis were estimated by maximum likelihood estimation using infection data from badgers and cattle collected during a 4-year badger vaccination trial. Our estimation showed that the environment can play an important role in the transmission of bTB, with a half-life of M. bovis in the environment of around 177 days. Based on the estimated transmission rate parameters, we calculate the basic reproduction ratio (R) within a herd, which reveals how relative badger density dictates transmission. In addition, we simulated transmission in each small local area to generate a first between-herd R map that identifies high-risk areas.Department of Agriculture, Food and the Marin

Efficient Direct and Limited Environmental Transmission of SARS-CoV-2 Lineage B.1.22 in Domestic Cats

The susceptibility of domestic cats to infection with SARS-CoV-2 has been demonstrated by several experimental studies and field observations. We performed an extensive study to further characterize the transmission of SARS-CoV-2 between cats, through both direct and indirect contact. To that end, we estimated the transmission rate parameter and the decay parameter for infectivity in the environment. Using four groups of pair-transmission experiment, all donor (inoculated) cats became infected, shed virus, and seroconverted, while three out of four direct contact cats got infected, shed virus, and two of those seroconverted. One out of eight cats exposed to a SARS-CoV-2-contaminated environment became infected but did not seroconvert. Statistical analysis of the transmission data gives a reproduction number R0 of 2.18 (95% CI = 0.92 to 4.08), a transmission rate parameter b of 0.23 day21 (95% CI = 0.06 to 0.54), and a virus decay rate parameter m of 2.73 day21 (95% CI = 0.77 to 15.82). These data indicate that transmission between cats is efficient and can be sustained (R0 . 1), however, the infectiousness of a contaminated environment decays rapidly (mean duration of infectiousness 1/2.73 days). Despite this, infections of cats via exposure to a SARS-CoV-2-contaminated environment cannot be discounted if cats are exposed shortly after contamination. IMPORTANCE This article provides additional insight into the risk of infection that could arise from cats infected with SARS-CoV-2 by using epidemiological models to determine transmission parameters. Considering that transmission parameters are not always provided in the literature describing transmission experiments in animals, we demonstrate that mathematical analysis of experimental data is crucial to estimate the likelihood of transmission. This article is also relevant to animal health professionals and authorities involved in risk assessments for zoonotic spill-overs of SARS-CoV-2. Last but not least, the mathematical models to calculate transmission parameters are applicable to analyze the experimental transmission of other pathogens between animals

Mass Mortality Caused by Highly Pathogenic Influenza A(H5N1) Virus in Sandwich Terns, the Netherlands, 2022

We collected data on mass mortality in Sandwich terns (Thalasseus sandvicensis) during the 2022 breeding season in the Netherlands. Mortality was associated with at least 2 variants of highly pathogenic avian influenza A(H5N1) virus clade 2.3.4.4b. We report on carcass removal efforts relative to survival in colonies. Mitigation strategies urgently require structured research

Update on the correlation of the highest energy cosmic rays with nearby extragalactic matter

Data collected by the Pierre Auger Observatory through 31 August 2007 showed evidence for anisotropy in the arrival directions of cosmic rays above the Greisen-Zatsepin-Kuz'min energy threshold, \nobreak{$6\times 10^{19}$eV}. The anisotropy was measured by the fraction of arrival directions that are less than $3.1^\circ$ from the position of an active galactic nucleus within 75 Mpc (using the V\'eron-Cetty and V\'eron $12^{\rm th}$ catalog). An updated measurement of this fraction is reported here using the arrival directions of cosmic rays recorded above the same energy threshold through 31 December 2009. The number of arrival directions has increased from 27 to 69, allowing a more precise measurement. The correlating fraction is $(38^{+7}_{-6})$, compared with $21$ expected for isotropic cosmic rays. This is down from the early estimate of $(69^{+11}_{-13})$. The enlarged set of arrival directions is examined also in relation to other populations of nearby extragalactic objects: galaxies in the 2 Microns All Sky Survey and active galactic nuclei detected in hard X-rays by the Swift Burst Alert Telescope. A celestial region around the position of the radiogalaxy Cen A has the largest excess of arrival directions relative to isotropic expectations. The 2-point autocorrelation function is shown for the enlarged set of arrival directions and compared to the isotropic expectation.Comment: Accepted for publication in Astroparticle Physics on 31 August 201