Characterization of cephalosporin and fluoroquinolone resistant Enterobacterales from Irish farm waste by whole genome sequencing.

Background: The Enterobacterales are a group of Gram-negative bacteria frequently exhibiting extended antimicrobial resistance (AMR) and involved in the transmission of resistance genes to other bacterial species present in the same environment. Due to their impact on human health and the paucity of new antibiotics, the World Health Organization (WHO) categorized carbapenem resistant and ESBL-producing as critical. Enterobacterales are ubiquitous and the role of the environment in the transmission of AMR organisms or antimicrobial resistance genes (ARGs) must be examined in tackling AMR in both humans and animals under the one health approach. Animal manure is recognized as an important source of AMR bacteria entering the environment, in which resistant genes can accumulate. Methods: To gain a better understanding of the dissemination of third generation cephalosporin and fluoroquinolone resistance genes between isolates in the environment, we applied whole genome sequencing (WGS) to Enterobacterales (79 E. coli, 1 Enterobacter cloacae, 1 Klebsiella pneumoniae, and 1 Citrobacter gillenii) isolated from farm effluents in Ireland before (n = 72) and after (n = 10) treatment by integrated constructed wetlands (ICWs). DNA was extracted using the MagNA Pure 96 system (Roche Diagnostics, Rotkreuz, Switzerland) followed by WGS on a MiSeq platform (Illumina, Eindhoven, Netherlands) using v3 chemistry as 300-cycle paired-end runs. AMR genes and point mutations were identified and compared to the phenotypic results for better understanding of the mechanisms of resistance and resistance transmission. Results: A wide variety of cephalosporin and fluoroquinolone resistance genes (mobile genetic elements (MGEs) and chromosomal mutations) were identified among isolates that mostly explained the phenotypic AMR patterns. A total of 31 plasmid replicon types were identified among the 82 isolates, with a subset of them (n = 24), identified in E. coli isolates. Five plasmid replicons were confined to the Enterobacter cloacae isolate and two were confined to the Klebsiella pneumoniae isolate. Virulence genes associated with functions including stress, survival, regulation, iron uptake secretion systems, invasion, adherence and toxin production were identified. Conclusion: Our study showed that antimicrobial resistant organisms (AROs) can persist even following wastewater treatment and could transmit AMR of clinical relevance to the environment and ultimately pose a risk to human or animal health.

Whole blood transcriptome analysis in dairy calves experimentally challenged with bovine herpesvirus 1 (BoHV-1) and comparison to a bovine respiratory syncytial virus (BRSV) challenge.

Bovine herpesvirus 1 (BoHV-1), is associated with several clinical syndromes in cattle, among which bovine respiratory disease (BRD) is of particular significance. Despite the importance of the disease, there is a lack of information on the molecular response to infection via experimental challenge with BoHV-1. The objective of this study was to investigate the whole-blood transcriptome of dairy calves experimentally challenged with BoHV-1. A secondary objective was to compare the gene expression results between two separate BRD pathogens using data from a similar challenge study with BRSV. Holstein-Friesian calves (mean age (SD) = 149.2 (23.8) days; mean weight (SD) = 174.6 (21.3) kg) were either administered BoHV-1 inoculate (1 × 107/mL × 8.5 mL) (n = 12) or were mock challenged with sterile phosphate buffered saline (n = 6). Clinical signs were recorded daily from day (d) -1 to d 6 (post-challenge), and whole blood was collected in Tempus RNA tubes on d six post-challenge for RNA-sequencing. There were 488 differentially expressed (DE) genes (p < 0.05, False Discovery rate (FDR) < 0.10, fold change ≥2) between the two treatments. Enriched KEGG pathways (p < 0.05, FDR <0.05); included Influenza A, Cytokine-cytokine receptor interaction and NOD-like receptor signalling. Significant gene ontology terms (p < 0.05, FDR <0.05) included defence response to virus and inflammatory response. Genes that are highly DE in key pathways are potential therapeutic targets for the treatment of BoHV-1 infection. A comparison to data from a similar study with BRSV identified both similarities and differences in the immune response to differing BRD pathogens.

Genomic diversity, virulence and source of Campylobacter jejuni contamination in Irish poultry slaughterhouses by whole genome sequencing.

AIMS: The aim was to exploit whole genome sequencing (WGS) to assess genomic diversity, identify virulence genes and deduce the proportion of Campylobacter colonized broilers that directly contaminate their carcasses. METHODS AND RESULTS: Campylobacter jejuni isolates (107) from caeca and carcass neck skin samples (50 pairs from the same batch plus 7 individual caeca) sampled at three poultry slaughterhouses over a one-year period were selected for sequencing (MiSeq; Illumina). FastQ files were submitted to BioNumerics for analysis using the wgMLST scheme for allele calling. Campylobacter cgMLST and hierarchical clustering was performed by applying the single linkage algorithm. Sequence types (STs) were determined in silico from the WGS data and isolates were assigned into clonal complexes (CCs) using the Campylobacter PubMLST.org database. Virulence genes were determined by downloading core sequences from the virulence factor database (VFDB) and the National Center for Biotechnology Information (NCBI). A high degree of diversity was observed with 23 different STs identified. ST257 and CC-21 were the most common STs and CCs, respectively. cgMLST analysis suggested that 56% of carcass contamination was a direct result of contamination from caeca from the same batch. Virulence genes known to play a role in human C. jejuni infection were identified such as the wlaN gene and the genes associated with lipooligosaccharide synthesis, which were identified in 30% of isolates. CONCLUSIONS: Caecal colonization was the more plausible occurring source of C. jejuni contamination of broiler carcasses, compared with cross-contamination from another batch or the environment. The high rate of genetic diversity observed amongst caecal isolates is consistent with a wide variety of Campylobacter strains circulating in poultry flocks in Ireland. SIGNIFICANCE AND IMPACT OF STUDY: The results will further inform broiler processors and regulators about the influence and importance of on-farm colonization versus slaughterhouse cross-contamination and the relationship between C. jejuni in caeca and carcasses during processing.

Elucidation of the Host Bronchial Lymph Node miRNA Transcriptome Response to Bovine Respiratory Syncytial Virus.

Bovine respiratory disease (BRD) causes substantial morbidity and mortality, affecting cattle of all ages. One of the main causes of BRD is an initial inflammatory response to bovine respiratory syncytial virus (BRSV). MicroRNAs are novel and emerging non-coding small RNAs that regulate many biological processes and are implicated in various inflammatory diseases. The objective of the present study was to elucidate the changes in the bovine bronchial lymph node miRNA transcriptome in response to BRSV following an experimental viral challenge. Holstein-Friesian calves were either administered a challenge dose of BRSV (103.5 TCID50/ml × 15 ml) (n = 12) or were mock inoculated with sterile phosphate buffered saline (n = 6). Daily scoring of clinical signs was performed and calves were euthanized at day 7 post-challenge. Bronchial lymph nodes were collected for subsequent RNA extraction and sequencing (75 bp). Read counts for known miRNAs were generated using the miRDeep2 package using the UMD3.1 reference genome and the bovine mature miRNA sequences from the miRBase database (release 22). EdgeR was used for differential expression analysis and Targetscan was used to identify target genes for the differentially expressed (DE) miRNAs. Target genes were examined for enriched pathways and gene ontologies using Ingenuity Pathway Analysis (Qiagen). Multi-dimensional scaling (MDS) based on miRNA gene expression changes, revealed a clearly defined separation between the BRSV challenged and control calves, although the clinical manifestation of disease was only mild. One hundred and nineteen DE miRNAs (P < 0.05, FDR < 0.1, fold change > 1.5) were detected between the BRSV challenged and control calves. The DE miRNAs were predicted to target 465 genes which were previously found to be DE in bronchial lymph node tissue, between these BRSV challenged and control calves. Of the DE predicted target genes, 455 had fold changes that were inverse to the corresponding DE miRNAs. There were eight enriched pathways among the DE predicted target genes with inverse fold changes to their corresponding DE miRNA including: granulocyte and agranulocyte adhesion and diapedesis, interferon signalling and role of pathogen recognition receptors in recognition of bacteria and viruses. Functions predicted to be increased included: T cell response, apoptosis of leukocytes, immune response of cells and stimulation of cells. Pathogen recognition and proliferation of cytotoxic T cells are vital for the recognition of the virus and its subsequent elimination.

Messenger RNA biomarkers of Bovine Respiratory Syncytial Virus infection in the whole blood of dairy calves.

Bovine Respiratory Syncytial Virus (BRSV) is a primary viral cause of Bovine Respiratory Disease (BRD) in young calves, which is responsible for substantial morbidity and mortality. Infection with BRSV induces global gene expression changes in respiratory tissues. If these changes are observed in tissues which are more accessible in live animals, such as whole blood, they may be used as biomarkers for diagnosis of the disease. Therefore, the objective of the current study was to elucidate the whole blood transcriptomic response of dairy calves to an experimental challenge with BRSV. Calves (Holstein-Friesian) were either administered BRSV inoculate (103.5 TCID50/ml × 15 ml) (n = 12) or sterile phosphate buffered saline (n = 6). Clinical signs were scored daily and whole blood was collected in Tempus RNA tubes immediately prior to euthanasia, at day 7 post-challenge. RNA was extracted from blood and sequenced (150 bp paired-end). The sequence reads were aligned to the bovine reference genome (UMD3.1) and EdgeR was subsequently employed for differential gene expression analysis. Multidimensional scaling showed that samples from BRSV challenged and control calves segregated based on whole blood gene expression changes, despite the BRSV challenged calves only displaying mild clinical symptoms of the disease. There were 281 differentially expressed (DE) genes (p < 0.05, FDR < 0.1, fold change > 2) between the BRSV challenged and control calves. The top enriched KEGG pathways and gene ontology terms were associated with viral infection and included "Influenza A", "defense response to virus", "regulation of viral life cycle" and "innate immune response". Highly DE genes involved in these pathways may be beneficial for the diagnosis of subclinical BRD from blood samples.

ATAC-Seq identifies regions of open chromatin in the bronchial lymph nodes of dairy calves experimentally challenged with bovine respiratory syncytial virus.

BACKGROUND: Bovine Respiratory Syncytial Virus (BRSV) is a cause of Bovine Respiratory Disease (BRD). DNA-based biomarkers contributing to BRD resistance are potentially present in non-protein-coding regulatory regions of the genome, which can be determined using ATAC-Seq. The objectives of this study were to: (i) identify regions of open chromatin in DNA extracted from bronchial lymph nodes (BLN) of healthy dairy calves experimentally challenged with BRSV and compare them with those from non-challenged healthy control calves, (ii) elucidate the chromatin regions that were differentially or uniquely open in the BRSV challenged relative to control calves, and (iii) compare the genes found in regions proximal to the differentially open regions to the genes previously found to be differentially expressed in the BLN in response to BRSV and to previously identified BRD susceptibility loci. This was achieved by challenging clinically healthy Holstein-Friesian calves (mean age 143 ± 14 days) with either BRSV inoculum (n = 12) or with sterile phosphate buffered saline (PBS) (n = 6) and preparing and sequencing ATAC-Seq libraries from fresh BLN tissues. RESULTS: Using Diffbind, 9,144 and 5,096 differentially accessible regions (P < 0.05, FDR < 0.05) were identified between BRSV challenged and control calves employing DeSeq2 and EdgeR, respectively. Additionally, 8,791 chromatin regions were found to be uniquely open in BRSV challenged calves. Seventy-six and 150 of the genes that were previously found to be differentially expressed using RNA-Seq, were located within 2 kb downstream of the differentially accessible regions, and of the regions uniquely open in BRSV challenged calves, respectively. Pathway analyses within ClusterProfiler indicated that these genes were involved in immune responses to infection and participated in the Th1 and Th2 pathways, pathogen recognition and the anti-viral response. There were 237 differentially accessible regions positioned within 40 previously identified BRD susceptibility loci. CONCLUSIONS: The identified open chromatin regions are likely to be involved in the regulatory response of gene transcription induced by infection with BRSV. Consequently, they may contain variants which impact resistance to BRD that could be used in breeding programmes to select healthier, more robust cattle.

Electrochemical assay of sorbitol dehydrogenase at PEDOT modified electrodes - a new milk biomarker for confirmation of pregnancy in dairy cattle.

A robust electrochemical assay for sorbitol dehydrogenase (SORDH) activity in milk was developed using voltammetry and chronocoulometry at bare and polymer modified transducers. The motivation for the work was to evaluate the potential of SORDH as an early biomarker of bovine pregnancy using milk as sample matrix. SORDH is an enzyme involved in carbohydrate metabolism converting sorbitol, the sugar alcohol form of glucose, into fructose, with NAD+ as a cofactor being simultaneously reduced to NADH. The assay was optimised via direct NADH oxidation on glassy carbon and screen printed carbon electrodes followed by electropolymerisation of 3,4-ethylenedioxythiophene (EDOT) monomer to form an NADH responsive PEDOT surface which operated well in undiluted milk samples. Assay conditions such as incubation time and temperature were optimised resulting in a 3 min assay at 37 °C in the presence of 10 mM NAD+ and 20 mM sorbitol co-substrates, enabling NADH electro oxidation (linear range 0.25-5 mM, sensitivity 9.17 µC cm-2 mM-1 in undiluted milk). SORDH determination followed over the range 0.31-10 U mL-1 in milk samples with sensitivity 5.45 µC cm-2 U-1 mL with LOD 0.0787 U mL-1. The assay was applied to milk sample testing acquired as part of an approved animal study involving control and breeding cycles of dairy cows with focus on analysis at day 19 post artificial insemination. Significant differences between control and pregnant SORDH levels in whole milk animal samples were found (average values 2.57 and 4.07 ng mL-1 respectively), as verified using a commercial SORDH ELISA optical assay. Finally, progesterone monitoring over days 16-21 of the oestrous cycle employed an optical ELISA assay and confirmed maintenance of progesterone levels from day 19 onwards.

Bovine hepatic miRNAome profiling and differential miRNA expression analyses between beef steers with divergent feed efficiency phenotypes.

MicroRNAs (miRNAs) are small RNA molecules involved in regulation of multiple biological processes through modulating expression of their target genes. Here we employed RNAseq to profile liver tissue miRNAome of 60 steers from Angus, Charolais, and Kinsella Composite (KC) populations. Of these animals, 36 animals (n = 12 for each breed) were utilized to identify differentially expressed (DE) miRNAs between animals with high (n = 6) or low (n = 6) phenotypic values of residual feed intake (RFI), a common measurement of feed efficiency. At a threshold of fold-change > 1.5 and P-value < 0.05, we detected 12 (7 up- and 5 downregulated in low-RFI animals), 18 (12 up- and 6 downregulated), and 13 (8 up- and 5 downregulated) DE miRNAs for Angus, Charolais, and KC steers, respectively. Most of the DE miRNAs were breed specific, with bta-miR-449a and bta-miR-AB-2 being differentially expressed in all three breeds. The predicted target genes of the identified DE miRNA are mainly involved in cell cycle, cell death and survival, cell signaling, cellular growth and proliferation, protein trafficking, cell morphology, cell-to-cell signaling and interaction, cellular development, molecular transport, post-translational modification, as well as nutrient metabolism (lipids, carbohydrates, protein and amino acid). Our results provide insights into the bovine hepatic miRNAome and their potential roles in molecular regulation of RFI in beef cattle.

Genome wide association study of passive immunity and disease traits in beef-suckler and dairy calves on Irish farms.

Calves with lower concentrations of immunoglobulin G (IgG) in their blood, have a greater risk of developing diseases. There is a lack of knowledge on genetic markers known to be associated with immunological variability or disease resistance. Therefore, the objective of this study was to identify SNP markers associated with passive immunity measures (serum IgG, serum protein, albumin, globulin and total protein concentrations, total solids Brix percentage, zinc sulphate turbidity units) and disease (pneumonia, diarrhoea, crude illness) traits in Irish commercial beef-suckler and dairy calves through genome wide association studies (GWAS). Genotyping was performed on DNA samples from beef-suckler (n = 698) and dairy (n = 1178) calves, using the IDBv3 chip. Heritability of passive immunity associated traits (range 0.02-0.22) and the disease traits (range 0.03-0.20) were low-to-moderate. Twenty-five and fifteen SNPs approached genome wide significance (P < 5 × 10-5) for the passive immunity and the disease traits, respectively. One SNP "ARS-BFGL-BAC-27914" reached Bonferroni genome wide significance (P < 1.15 × 10-6) for an association with serum IgG concentration in beef calves. Further work will evaluate these SNPs in larger cattle populations and assess their contribution to genomic selection breeding strategies, aimed towards producing more disease resistant livestock.

Experimental challenge with bovine respiratory syncytial virus in dairy calves: bronchial lymph node transcriptome response.

Bovine Respiratory Disease (BRD) is the leading cause of mortality in calves. The objective of this study was to examine the response of the host's bronchial lymph node transcriptome to Bovine Respiratory Syncytial Virus (BRSV) in a controlled viral challenge. Holstein-Friesian calves were either inoculated with virus (103.5 TCID50/ml × 15 ml) (n = 12) or mock challenged with phosphate buffered saline (n = 6). Clinical signs were scored daily and blood was collected for haematology counts, until euthanasia at day 7 post-challenge. RNA was extracted and sequenced (75 bp paired-end) from bronchial lymph nodes. Sequence reads were aligned to the UMD3.1 bovine reference genome and differential gene expression analysis was performed using EdgeR. There was a clear separation between BRSV challenged and control calves based on gene expression changes, despite an observed mild clinical manifestation of the disease. Therefore, measuring host gene expression levels may be beneficial for the diagnosis of subclinical BRD. There were 934 differentially expressed genes (DEG) (p < 0.05, FDR <0.1, fold change >2) between the BRSV challenged and control calves. Over-represented gene ontology terms, pathways and molecular functions, among the DEG, were associated with immune responses. The top enriched pathways included interferon signaling, granzyme B signaling and pathogen pattern recognition receptors, which are responsible for the cytotoxic responses necessary to eliminate the virus.

Risk factors associated with exposure to bovine respiratory disease pathogens during the peri-weaning period in dairy bull calves.

BACKGROUND: Bovine respiratory disease (BRD) remains among the leading causes of death of cattle internationally. The objective of this study was to identify risk factors associated with exposure to BRD pathogens during the peri-weaning period (day (d)-14 to d 14 relative to weaning at 0) in dairy bull calves using serological responses to these pathogens as surrogate markers of exposure. Clinically normal Holstein-Friesian and Jersey breed bull calves (n = 72) were group housed in 4 pens using a factorial design with calves of different breeds and planes of nutrition in each pen. Intrinsic, management and clinical data were collected during the pre-weaning (d - 56 to d - 14) period. Calves were gradually weaned over 14 days (d - 14 to d 0). Serological analysis for antibodies against key BRD pathogens (BRSV, BPI3V, BHV-1, BHV-4, BCoV, BVDV and H. somni) was undertaken at d - 14 and d 14. Linear regression models (for BVDV, BPI3V, BHV-1, BHV-4, BCoV and H. somni) and a single mixed effect random variable model (for BRSV) were used to identify risk factors for changes in antibody levels to these pathogens. RESULTS: BRSV was the only pathogen which demonstrated clustering by pen. Jersey calves experienced significantly lower changes in BVDV S/P than Holstein-Friesian calves. Animals with a high maximum respiratory score (≥8) recorded significant increases in H. somni S/P during the peri-weaning period when compared to those with respiratory scores of ≤3. Haptoglobin levels of between 1.32 and 1.60 mg/ml at d - 14 were significantly associated with decreases in BHV-1 S/N during the peri-weaning period. Higher BVDV S/P ratios at d - 14 were significantly correlated with increased changes in serological responses to BHV-4 over the peri-weaning period. CONCLUSIONS: Haptoglobin may have potential as a predictor of exposure to BHV-1. BRSV would appear to play a more significant role at the 'group' rather than 'individual animal' level. The significant associations between the pre-weaning levels of antibodies to certain BRD pathogens and changes in the levels of antibodies to the various pathogens during the peri-weaning period may reflect a cohort of possibly genetically linked 'better responders' among the study population.

Illumina MiSeq 16S amplicon sequence analysis of bovine respiratory disease associated bacteria in lung and mediastinal lymph node tissue.

BACKGROUND: Bovine respiratory disease (BRD) is caused by growth of single or multiple species of pathogenic bacteria in lung tissue following stress and/or viral infection. Next generation sequencing of 16S ribosomal RNA gene PCR amplicons (NGS 16S amplicon analysis) is a powerful culture-independent open reference method that has recently been used to increase understanding of BRD-associated bacteria in the upper respiratory tract of BRD cattle. However, it has not yet been used to examine the microbiome of the bovine lower respiratory tract. The objective of this study was to use NGS 16S amplicon analysis to identify bacteria in post-mortem lung and lymph node tissue samples harvested from fatal BRD cases and clinically healthy animals. Cranial lobe and corresponding mediastinal lymph node post-mortem tissue samples were collected from calves diagnosed as BRD cases by veterinary laboratory pathologists and from clinically healthy calves. NGS 16S amplicon libraries, targeting the V3-V4 region of the bacterial 16S rRNA gene were prepared and sequenced on an Illumina MiSeq. Quantitative insights into microbial ecology (QIIME) was used to determine operational taxonomic units (OTUs) which corresponded to the 16S rRNA gene sequences. RESULTS: Leptotrichiaceae, Mycoplasma, Pasteurellaceae, and Fusobacterium were the most abundant OTUs identified in the lungs and lymph nodes of the calves which died from BRD. Leptotrichiaceae, Fusobacterium, Mycoplasma, Trueperella and Bacteroides had greater relative abundances in post-mortem lung samples collected from fatal cases of BRD in dairy calves, compared with clinically healthy calves without lung lesions. Leptotrichiaceae, Mycoplasma and Pasteurellaceae showed higher relative abundances in post-mortem lymph node samples collected from fatal cases of BRD in dairy calves, compared with clinically healthy calves without lung lesions. Two Leptotrichiaceae sequence contigs were subsequently assembled from bacterial DNA-enriched shotgun sequences. CONCLUSIONS: The microbiomes of the cranial lung lobe and mediastinal lymph node from calves which died from BRD and from clinically healthy H-F calves have been characterised. Contigs corresponding to the abundant Leptotrichiaceae OTU were sequenced and found not to be identical to any known bacterial genus. This suggests that we have identified a novel bacterial species associated with BRD.

Early recognition of bovine respiratory disease in calves using automated continuous monitoring of cough sounds.

Bovine respiratory disease (BRD) complex in calves impairs health and welfare and causes severe economic losses for the Stockperson. Early recognition of BRD should lead to earlier veterinary (antibiotic/anti-inflammatory) treatment interventions thereby reducing the severity of the disease and associated costs. Coughing is one of the clinical manifestations of BRD. It is believed that by automatically and continuously monitoring the sounds within calf houses, and analysing the coughing frequency, early recognition of BRD in calves is possible. Therefore, the objective of the present study was to develop an automated calf cough monitor and examine its potential as an early warning system for BRD in artificially reared dairy calves. The coughing sounds of 62 calves were continuously recorded by a microphone over a three-month period. A sound analysis algorithm was developed to distinguish calf coughs from other sounds (e.g. mechanical sounds). During the sound recording period the health of the calves was assessed and scored periodically per week by a trained human observer. Calves presenting with BRD received antibiotic and/or anti-inflammatory treatment and the dates of treatment were recorded. This treatment date reference served as a comparison for the investigation of whether an increase in coughing frequency could be related to calves developing BRD. The calf cough detection algorithm achieved 50.3% sensitivity, 99.2% specificity and 87.5% precision. Four out of five periods, where coughing frequency was observed to be increased, coincided with the development of BRD in more than one calf. This period of increased coughing frequency was always observed before the calves were treated. Therefore, the calf cough monitor has the potential to identify early onset of BRD in calves.