Transcriptomic Analysis of Circulating Leukocytes Obtained during the Recovery from Clinical Mastitis Caused by Escherichia coli in Holstein Dairy Cows

Simple Summary Escherichia coli is a bacterium which infects cow udders causing clinical mastitis, a potentially severe disease with welfare and economic consequences. During an infection, white blood cells (leukocytes) enter the udder to provide immune defence and assist tissue repair. We sequenced RNA derived from circulating leukocytes to investigate which genes are up- or down-regulated in dairy cows with naturally occurring cases of clinical mastitis in comparison with healthy control cows from the same farm. We also looked for genetic variations between infected and healthy cows. Blood samples were taken either EARLY (around 10 days) or LATE (after 4 weeks) during the recovery phase after diagnosis. Many genes (1090) with immune and inflammatory functions were up-regulated during the EARLY phase. By the LATE phase only 29 genes were up-regulated including six haemoglobin subunits, possibly important for the production of new red blood corpuscles. Twelve genetic variations which were associated with an increased or decreased expression of some important immune genes were identified between the infected and control cows. These results show that the initial inflammatory response to E. coli continued for at least 10 days despite the cows having received prompt veterinary treatment, but they had largely recovered within 4 weeks. Genetic differences between cows may predispose some animals to infection. Abstract The risk and severity of clinical infection with Escherichia coli as a causative pathogen for bovine mastitis is influenced by the hosts’ phenotypic and genotypic variables. We used RNA-Seq analysis of circulating leukocytes to investigate global transcriptomic profiles and genetic variants from Holstein cows with naturally occurring cases of clinical mastitis, diagnosed using clinical symptoms and milk microbiology. Healthy lactation-matched cows served as controls (CONT, n = 6). Blood samples were collected at two time periods during the recovery phase post diagnosis: EARLY (10.3 ± 1.8 days, n = 6) and LATE (46.7 ± 11 days, n = 3). Differentially expressed genes (DEGs) between the groups were identified using CLC Genomics Workbench V21 and subjected to enrichment analysis. Variant calling was performed following GATKv3.8 best practice. The comparison of E. coli(+) EARLY and CONT cows found the up-regulation of 1090 DEGs, mainly with immune and inflammatory functions. The key signalling pathways involved NOD-like and interleukin-1 receptors and chemokines. Many up-regulated DEGs encoded antimicrobial peptides including cathelicidins, beta-defensins, S100 calcium binding proteins, haptoglobin and lactoferrin. Inflammation had largely resolved in the E. coli(+) LATE group, with only 29 up-regulated DEGs. Both EARLY and LATE cows had up-regulated DEGs encoding ATP binding cassette (ABC) transporters and haemoglobin subunits were also up-regulated in LATE cows. Twelve candidate genetic variants were identified in DEGs between the infected and CONT cows. Three were in contiguous genes WIPI1, ARSG and SLC16A6 on BTA19. Two others (RAC2 and ARHGAP26) encode a Rho-family GTPase and Rho GTPase-activating protein 26. These results show that the initial inflammatory response to E. coli continued for at least 10 days despite prompt treatment and provide preliminary evidence for genetic differences between cows that may predispose them to infection.