Search

Additional file 2: Figure S1. Manhattan plots of lesion susceptibility for sole haemorrhage (SH), sole ulcers (SU), and white line disease (WL). The standardized SNP effects were calculated by the estimated SNP effects divided by their empirical standard deviation. The data provided represent the Manhattan plots of the standardized SNP effects for lesion susceptibility of sole haemorrhage (SH), sole ulcers (SU), and white line disease (WL) from the genome-wide association analysis.

Additional file 1: Table S1. The top ten markers associated with lesion severity for sole haemorrhage (SH), sole ulcer (SU), and white line disease (WL) from single-marker GWA analyses, including chromosome (BTA), position (Pos), minor allele frequency (MAF), and P-value. The data provided represent the position, allele frequency, and the P value of the top ten markers associated with lesion severity for sole haemorrhage (SH), sole ulcer (SU), and white line disease (WL) from the genome-wide association analysis. Table S2. The top ten markers with the largest standardized SNP effects (Std_SNP_effect) on lesion susceptibility for sole haemorrhage (SH), sole ulcer (SU), and white line disease (WL) from single-marker GWA analyses, including chromosome (BTA), position (Pos), and minor allele frequency. The data provided represent the position, allele frequency, and the effect size of the top ten markers associated with lesion susceptibility for sole haemorrhage (SH), sole ulcer (SU), and white line disease (WL) from the genome-wide association analysis.

This study aimed to determine the association between the lameness advantage genetic index and four outcomes: sole haemorrhage (SH), sole ulcers (SU), white line lesions (WL), and lameness during mobility scoring.We enrolled 2352 Holstein cows from four predominantly housed dairy herds in the UK. Cows were mobility scored and foot lesions recorded at four time points from before calving to late lactation. Cows were genotyped and genetic indexes were assigned to each cow following national genetic evaluations. Lameness records and genetic indexes were matched for 2107 cows. Four separate multivariable logistic regression models, which included farm and parity as covariables, were used to quantify the association between the lameness advantage index and whether animals were affected by SH, SU, WL, or lameness.The odds ratios (95% confidence intervals) for one-point increases in the lameness advantage index were 0.79 (0.72-0.86), 0.68 (0.59-0.78), 0.94 (0.84-1.04), and 0.82 (0.74-0.91) for SH, SU, WL, and lameness, respectively. The same trends were present when the sire's lameness advantage index was evaluated in place of the animal's own, although the strength of this association was generally weaker.The lameness advantage index is associated with SH, SU, and lameness, therefore selection on the lameness advantage index could be considered in herds aiming to reduce lameness. Where genomic testing of heifers is not conducted, sire lameness advantage index may still be effective to reduce SH and SU incidence.

MHCY is a candidate region for influencing immune responses in chickens. MHCY contains multiple specialized, polymorphic MHC class I loci along with loci belonging to 4 additional gene families. In this study, MHCY haplotypes were tested for association with cecal colonization after Campylobacter jejuni infection of a backcross [(Line 61 × Line N) × Line N] population derived from 2 White Leghorn research lines, Line 61 and Line N, that were previously shown to exhibit heritable differences in colonization. Samples were obtained for 51 birds challenged with 108 CFU Campylobacter jejuni at 3 wk of age. Viable C. jejuni in the ceca were enumerated 5 d postinfection and counts were log-transformed for analysis. Birds were assigned to either low or high colonization groups based on the individual count being below or above the mean bacterial count for all birds. The mean bacterial count of the low infection group differed significantly from the high infection group. Sex and MHCB haplotype had similar distributions within the 2 groups. Overall, 7 MHCY haplotypes were found to be segregating. Two were significantly associated with C. jejuni colonization. MHCY Y18 was associated with low colonization (P = 3.00 × 10-5); whereas MHCY Y11a was associated with high colonization (P = 0.008). The MHCY haplotype impacted the mean bacterial count among all birds with MHCY Y18 having the lowest bacterial count compared with MHCY Y11a and all other MHCY (Y5, Y7, Y8, Y11b, and Y11c) haplotypes. These findings support further investigation of the contribution of chicken MHCY in resistance to Campylobacter colonization.

MHCY is a candidate region for influencing immune responses in chickens. MHCY contains multiple specialized, polymorphic MHC class I loci along with loci belonging to 4 additional gene families. In this study, MHCY haplotypes were tested for association with cecal colonization after Campylobacter jejuni infection of a backcross [(Line 61 × Line N) × Line N] population derived from 2 White Leghorn research lines, Line 61 and Line N, that were previously shown to exhibit heritable differences in colonization. Samples were obtained for 51 birds challenged with 108 CFU Campylobacter jejuni at 3 wk of age. Viable C. jejuni in the ceca were enumerated 5 d postinfection and counts were log-transformed for analysis. Birds were assigned to either low or high colonization groups based on the individual count being below or above the mean bacterial count for all birds. The mean bacterial count of the low infection group differed significantly from the high infection group. Sex and MHCB haplotype had similar distributions within the 2 groups. Overall, 7 MHCY haplotypes were found to be segregating. Two were significantly associated with C. jejuni colonization. MHCY Y18 was associated with low colonization (P = 3.00 × 10−5); whereas MHCY Y11a was associated with high colonization (P = 0.008). The MHCY haplotype impacted the mean bacterial count among all birds with MHCY Y18 having the lowest bacterial count compared with MHCY Y11a and all other MHCY (Y5, Y7, Y8, Y11b, and Y11c) haplotypes. These findings support further investigation of the contribution of chicken MHCY in resistance to Campylobacter colonization.

The objective of this study was to investigate the association between (sub)clinical mastitis (CM) in the first 30 days in milk (DIM) and the presence of sole ulcers (SU) later in lactation.Holstein cows and heifers were examined for presence of sole haemorrhage and SU before calving, in the first 14 days postcalving and in early lactation (after 30 DIM). CM episodes and somatic cell counts (SCC) measurements were obtained from farm records. Multivariable logistic regression was used for data analysis.Odds of SU in early lactation were 2.44 times greater (95% confidence interval [CI] 0.97-5.54) in cows that had CM in the first 30 DIM compared to cows that did not have CM in the first 30 DIM. When cows that had SU precalving or at the calving check were excluded from the dataset, an association of CM in the first 30 DIM with later presence of SU was no longer statistically significant but the same numeric trend still existed (odds ratio [OR] 2.25, 95% CI 0.81-5.34). The odds of SU in early lactation were 1.70 times greater in cows that had high SCC compared to cows that did not have high SCC in the first 100 DIM (95% CI 1.13-2.55).An association was found between CM in the first 30 DIM and presence of SU in early lactation (after 30 DIM). Elucidating the mechanism behind this relationship could improve our understanding of the aetiopathogenesis of both diseases and lead to new preventive strategies.

The study of host-pathogen interactions has been accelerating rapidly in both medical and veterinary contexts due to methodological advances and increasing infectious disease threats. This Research Topic attracted a wide breadth of submissions spanning analyses of host-pathogen interactions using in vitro and in vivo models of human disease; through to a range of studies in livestock species, comprising cattle, pigs, chickens, and farmed fish. An impressive breadth of technical approaches was also in evidence, from GWAS, SNP analysis, transcriptomics, RNAi, proteomics, and reverse genetics; often deployed in combination in a single study. Overall, there was as a strong bias toward “host-centric” articles (n = 12) compared with “pathogen-centric” (n = 6) studies, perhaps reflecting the worldwide drive to improve disease resistance in livestock.

Background Campylobacter jejuni is the leading cause of bacterial gastroenteritis in humans and the handling or consumption of contaminated poultry meat is a key source of infection. Selective breeding of poultry that exhibit elevated resistance to Campylobacter is an attractive control strategy. Here we studied the global transcriptional response of inbred chicken lines that differ in resistance to C. jejuni colonisation at a key site of bacterial persistence. Results Three-week-old chickens of line 61 and N were inoculated orally with C. jejuni strain M1 and caecal contents and tonsils were sampled at 1 and 5 days post-infection. Caecal colonisation was significantly lower in line 61 compared to line N at 1 day post-infection, but not 5 days post-infection. RNA-Seq analysis of caecal tonsils of both lines revealed a limited response to C. jejuni infection compared to age-matched uninfected controls. In line N at days 1 and 5 post-infection, just 8 and 3 differentially expressed genes (DEGs) were detected (fold-change > 2 and false-discovery rate of 1, a broader response to C. jejuni was observed, with 69 DEGs relating to immune regulation, cell signalling and metabolism at 1 day post-infection. However, by day 5 post-infection, no DEGs were detected. By far, the greatest number of DEGs were between uninfected birds of the two lines implying that differential resistance to C. jejuni is intrinsic. Of these genes, several Major Histocompatibility Complex class I-related genes (MHCIA1, MHCBL2 and MHCIY) and antimicrobial peptides (MUC2, AvBD10 and GZMA) were expressed to a greater extent in line N. Two genes within quantitative trait loci associated with C. jejuni colonisation were also more highly expressed in line N (ASIC4 and BZFP2). Quantitative reverse-transcriptase PCR analysis of a subset of transcripts confirmed the RNA-Seq results. Conclusions Our data indicate a limited transcriptional response in the caecal tonsils of inbred chickens to intestinal colonisation by Campylobacter but identify a large number of differentially transcribed genes between lines 61 and N that may underlie variation in heritable resistance to C. jejuni.

Canine Inflammatory bowel disease (IBD) is a chronic multifactorial disease, resulting from complex interactions between the intestinal immune system, microbiota and environmental factors in genetically predisposed dogs. Previously, we identified several single nucleotide polymorphisms (SNP) and regions on chromosomes (Chr) 7, 9, 11 and 13 associated with IBD in German shepherd dogs (GSD) using GWAS and FST association analyses. Here, building on our previous results, we performed a targeted next-generation sequencing (NGS) of a two Mb region on Chr 9 and 11 that included 14 of the newly identified candidate genes, in order to identify potential functional SNPs that could explain these association signals. Furthermore, correlations between genotype and treatment response were estimated. Results revealed several SNPs in the genes for canine EEF1A1, MDH2, IL3, IL4, IL13 and PDLIM, which, based on the known function of their corresponding proteins, further our insight into the pathogenesis of IBD in dogs. In addition, several pathways involved in innate and adaptive immunity and inflammatory responses (i.e. T helper cell differentiation, Th1 and Th2 activation pathway, communication between innate and adaptive immune cells and differential regulation of cytokine production in intestinal epithelial cells by IL-17A and IL-17F), were constructed involving the gene products in the candidate regions for IBD susceptibility. Interestingly, some of the identified SNPs were present in only one outcome group, suggesting that different genetic factors are involved in the pathogenesis of IBD in different treatment response groups. This also highlights potential genetic markers to predict the response in dogs treated for IBD.

Application: Exertional rhabdomyolysis is a common disease of horses (prevalence 7–10% in racehorses [5], and presents both a significant welfare concern and huge financial cost to equine sport. Prophylaxis and disease management are challenging, and aetiopathogenesis remains unclear

Application: The equine muscle disorder exertional rhabdomyolysis (ER) is an animal welfare concern that has substantial financial implications for equine industries. Breeders may inadvertently be selecting for the disease (Isgren et al., 2010), and identifying genetic markers could inform future breeding programmes. Recognition of specific clinical presentations across different breeds might reveal separate causes enabling more refined treatments and prophylaxis.

Our aims were to (1) determine how interdigital skin temperature (IST), measured using infrared thermography, was associated with different stages of digital dermatitis (DD) lesions and (2) develop and validate models that can use IST measurements to identify cows with an active DD lesion. Between March 2019 and March 2020, infrared thermographic images of hind feet were taken from 2,334 Holstein cows across 4 farms. We recorded the maximum temperature reading from infrared thermographic images of the interdigital skin between the heel bulbs on the hind feet. Pregnant animals were enrolled approximately 1 to 2 mo precalving, reassessed 1 wk after calving, and again at approximately 50 to 100 d postpartum. At these time points, IST and the clinical stage of DD (M-stage scoring system: M1–M4.1) were recorded in addition to other data such as the ambient environmental temperature, height, body condition score, parity, and the presence of other foot lesions. A mixed effect linear regression model with IST as the dependent variable was fitted. Interdigital skin temperature was associated with DD lesions; compared to healthy feet, IST was highest in feet with M2 lesions, followed by M1 and M4.1 lesions. Subsequently, the capacity of IST measurements to detect the presence or absence of an active DD lesion (M1, M2, or M4.1) was explored by fitting logistic regression models, which were tested using 10-fold validation. A mixed effect logistic regression model with the presence of active DD as the dependent variable was fitted first. The average area under the curve for this model was 0.80 when its ability to detect presence of active DD was tested on 10% of the data that were not used for the model's training; an average sensitivity of 0.77 and an average specificity of 0.67 was achieved. This model was then restricted so that only explanatory variables that could be practically recorded in a nonresearch, external setting were included. Validation of this model demonstrated an average area under the curve of 0.78, a sensitivity of 0.88, and a specificity of 0.66 for 1 of the time points (precalving). Lower sensitivity and specificity were achieved for the other 2 time points. Our study adds further evidence to the relationship between DD and foot skin temperature using a large data set with multiple measurements per animal. Additionally, we highlight the potential for infrared thermography to be used for routine on-farm diagnosis of active DD lesions.

Funder: Biotechnology, Campylobacter is the leading cause of bacterial foodborne gastroenteritis worldwide. Handling or consumption of contaminated poultry meat is a key risk factor for human campylobacteriosis. One potential control strategy is to select poultry with increased resistance to Campylobacter. We associated high-density genome-wide genotypes (600K single nucleotide polymorphisms) of 3000 commercial broilers with Campylobacter load in their caeca. Trait heritability was modest but significant (h2 = 0.11 ± 0.03). Results confirmed quantitative trait loci (QTL) on chromosomes 14 and 16 previously identified in inbred chicken lines, and detected two additional QTLs on chromosomes 19 and 26. RNA-Seq analysis of broilers at the extremes of colonisation phenotype identified differentially transcribed genes within the QTL on chromosome 16 and proximal to the major histocompatibility complex (MHC) locus. We identified strong cis-QTLs located within MHC suggesting the presence of cis-acting variation in MHC class I and II and BG genes. Pathway and network analyses implicated cooperative functional pathways and networks in colonisation, including those related to antigen presentation, innate and adaptive immune responses, calcium, and renin-angiotensin signalling. While co-selection for enhanced resistance and other breeding goals is feasible, the frequency of resistance-associated alleles was high in the population studied and non-genetic factors significantly influenced Campylobacter colonisation.

Telomere length is predictive of adult health and survival across vertebrate species. However, we currently do not know whether such associations result from among-individual differences in telomere length determined genetically or by early-life environmental conditions, or from differences in the rate of telomere attrition over the course of life that might be affected by environmental conditions. Here, we measured relative leukocyte telomere length (RLTL) multiple times across the entire lifespan of dairy cattle in a research population that is closely monitored for health and milk production and where individuals are predominantly culled in response to health issues. Animals varied in their change in RLTL between subsequent measurements and RLTL shortened more during early life and following hotter summers which are known to cause heat stress in dairy cows. The average amount of telomere attrition calculated over multiple repeat samples of individuals predicted a shorter productive lifespan, suggesting a link between telomere loss and health. TL attrition was a better predictor of when an animal was culled than their average TL or the previously for this population reported significant TL at the age of 1 year. Our present results support the hypothesis that TL is a flexible trait that is affected by environmental factors and that telomere attrition is linked to animal health and survival traits. Change in telomere length may represent a useful biomarker in animal welfare studies.

Telomere length measured in blood cells is predictive of subsequent adult health and survival across a range of vertebrate species. However, we currently do not know whether such associations result from among-individual differences in telomere length determined genetically or by environmental factors early in life, or from differences in the rate of telomere attrition over the course of life. Here, we measured relative leukocyte telomere length (RLTL) multiple times across the entire lifespan of dairy cattle in a research population that is closely monitored for health and milk production and where individuals are only culled in response to health issues and less due to poor milk production than on purely commercial farms. Our results clearly show that the average amount of telomere attrition over an individual’s life, not their average or early life telomere length predicted when an individual was culled. Within-individual telomere length attrition could reflect environmental or physiological insults which may accumulate to predict individual health-span. We also show that animals with more telomere attrition in their first year of life were culled at a younger age, indicating that early life stressors may have a prolonged effect on adult life.