Complete blood count data and leukocyte expression of cytokine genes and cytokine receptor genes associated with bovine respiratory disease in calves

Differential gene expression in peripheral leukocytes of pre-weaned Holstein heifer calves with respiratory disease

Bovine respiratory disease (BRD) is a leading cause of calf morbidity and mortality, and prevalence remains high despite current management practices. Differential gene expression (DGE) provides detailed insight into individual immune responses and can illuminate enriched pathways and biomarkers that contribute to disease susceptibility and outcomes. The aims of this study were to investigate differences in peripheral leukocyte gene expression in Holstein preweaned heifer calves 1) with and without BRD, and 2) across weeks of age. Calves were enrolled for this short-term longitudinal study on two commercial dairies in Washington State. Calves were assessed every two weeks throughout the pre-weaning period using clinical respiratory scoring (CRS) and thoracic ultrasonography (TUS), and blood samples were collected. Calves were selected that were either healthy (n = 10) or had BRD diagnosed by CRS (n = 7), TUS (n = 6), or both (n = 6) in weeks 5 or 7 of life). Three consecutive time point samples were analyzed for each BRD calf consisting of PRE, ONSET, and POST samples. Nineteen genes of interest were selected based on previous gene expression studies in cattle: ALOX15, BPI, CATHL6, CXCL8, DHX58, GZMB, HPGD, IFNG, IL17D, IL1R2, ISG15, LCN2, LIF, MX1, OAS2, PGLYRP1, S100A8, SELP, and TNF. Comparisons were made between age and disease time point matched BRD and healthy calves as well as between calf weeks of age. No DGE was observed between diseased and healthy calves; however, DGE was observed between calf weeks of age regardless of disease state. Developmental differences in leukocyte gene expression, phenotype, and functionality make pre-weaned calves immunologically distinct from mature cattle, and early life shifts in calf leukocyte populations likely contribute to the age-related gene expression differences we observed. Age overshadows disease impacts to influence gene expression in young calves, and immune development progresses upon a common trajectory regardless of disease during the preweaning period.

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 × 10⁷/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.

Sexual Dimorphic Innate Immune Response to a Viral-Bacterial Respiratory Disease Challenge in Beef Calves

The potential for sexually dimorphic innate immune responses to respiratory disease was evaluated, where eight steers and seven heifers (280 ± 4 kg) were subjected to a viral−bacterial respiratory disease challenge utilizing bovine herpesvirus-1 (BoHV-1; intranasal; 1 × 108 PFU/nostril) and Mannheimia haemolytica (MH; intratracheal; 1.3 × 107 CFU/head) administered 72 h later. Body temperature was lesser in heifers than steers (p < 0.01). There was a sex × time interaction (p = 0.05) for white blood cells where heifers had reduced concentrations compared with steers at −72 and 0 h but greater concentrations from 36 to 60 h post-MH. Concentrations of neutrophils were lesser in heifers compared to steers from 0 to 4 h, and from 8 to 12 h (p = 0.03). Lymphocytes were greater in heifers compared to steers at 12 h and from 36 to 60 h post-MH (p < 0.01). The neutrophil−lymphocyte ratio was lesser in heifers compared to steers from 2 to 24 h and at 48 h post-MH (p < 0.01). Monocytes were greater in heifers compared to steers from 24 to 60 h post-MH (p < 0.01), while eosinophils were greater in heifers compared to steers at 48 and 60 h (p < 0.01). Serum IL-4 was lesser in heifers compared to steers at 0 h and from 2 to 72 h post-MH challenge (p = 0.02). Non-esterified fatty acid concentrations were lesser (p < 0.01) in heifers compared to steers from 2 to 4 h post-MH challenge. Urea nitrogen concentrations were greater (p < 0.01) in heifers than steers at 36 h post-MH challenge. Data from this study reveal distinct differences in the acute phase response following a respiratory disease challenge where steers produced an early response, while the response in heifers appeared to be delayed.

Hematological and gene co-expression network analyses of high-risk beef cattle defines immunological mechanisms and biological complexes involved in bovine respiratory disease and weight gain

Bovine respiratory disease (BRD), the leading disease complex in beef cattle production systems, remains highly elusive regarding diagnostics and disease prediction. Previous research has employed cellular and molecular techniques to describe hematological and gene expression variation that coincides with BRD development. Here, we utilized weighted gene co-expression network analysis (WGCNA) to leverage total gene expression patterns from cattle at arrival and generate hematological and clinical trait associations to describe mechanisms that may predict BRD development. Gene expression counts of previously published RNA-Seq data from 23 cattle (2017; n = 11 Healthy, n = 12 BRD) were used to construct gene co-expression modules and correlation patterns with complete blood count (CBC) and clinical datasets. Modules were further evaluated for cross-populational preservation of expression with RNA-Seq data from 24 cattle in an independent population (2019; n = 12 Healthy, n = 12 BRD). Genes within well-preserved modules were subject to functional enrichment analysis for significant Gene Ontology terms and pathways. Genes which possessed high module membership and association with BRD development, regardless of module preservation (“hub genes”), were utilized for protein-protein physical interaction network and clustering analyses. Five well-preserved modules of co-expressed genes were identified. One module (“steelblue”), involved in alpha-beta T-cell complexes and Th2-type immunity, possessed significant correlation with increased erythrocytes, platelets, and BRD development. One module (“purple”), involved in mitochondrial metabolism and rRNA maturation, possessed significant correlation with increased eosinophils, fecal egg count per gram, and weight gain over time. Fifty-two interacting hub genes, stratified into 11 clusters, may possess transient function involved in BRD development not previously described in literature. This study identifies co-expressed genes and coordinated mechanisms associated with BRD, which necessitates further investigation in BRD-prediction research.

The effects of administering different metaphylactic antimicrobials on growth performance and health outcomes of high-risk, newly received feedlot steers

Bovine respiratory disease (BRD) is the primary animal health concern facing feedlot producers. Many antimicrobial mitigation strategies are available, but few studies have compared feedlot performance during both the receiving and finishing periods following application of different antimicrobials used as metaphylaxis at arrival. The objective of this study was to compare antimicrobial metaphylaxis methods on clinical health and growth performance across both the receiving and finishing periods. A total of 238 multiple-sourced steers in two source blocks were used in a generalized complete block design. The four treatments included: 1) a negative control, 5 mL of sterile saline injected subcutaneously (CON); 2) subcutaneous administration of florfenicol at 40 mg/kg of BW (NUF); 3) subcutaneous administration of ceftiofur in the posterior aspect of the ear at 6.6 mg/kg of BW (EXC); and 4) subcutaneous administration of tulathromycin at 2.5 mg/kg of BW (DRA). The morbidity rate for the first treatment of BRD was decreased for the DRA and EXC treatments compared to CON and NUF (P &lt; 0.01). Additionally, average daily gain (ADG), dry matter intake (DMI), and gain-to-feed (G:F) were greater (P ≤ 0.02) in the DRA treatment during the receiving period compared to all other treatments. The ADG was also greater (P &lt; 0.05) for EXC than the CON treatment throughout the finishing period. Nonetheless, other growth performance variables did not differ among metaphylactic treatments during the finishing period (P ≥ 0.14). Likewise, no differences in carcass characteristics or liver abscess score were observed (P ≥ 0.18). All complete blood count (CBC) variables were affected by day (P ≤ 0.01) except mean corpuscular hemoglobin concentration (P = 0.29). Treatment × time interactions were observed for platelet count, white blood cell (WBC) count, monocyte count and percentage, and lymphocyte percentage (P ≤ 0.03). However, there were no observed hematological variables that differed among treatment (P ≥ 0.10). The results indicate that some commercially available antimicrobials labeled for metaphylactic use are more efficacious than others in decreasing morbidity rate.

Review of Associated Health Benefits of Algal Supplementation in Cattle with Reference to Bovine Respiratory Disease Complex in Feedlot Systems

Within the Australian beef industry bovine respiratory disease is considered one of the most common disease and costs the industry an average net loss of $1647.53 Australian dollars per animal death to bovine respiratory disease complex (BRD). This is due to the disease overwhelming the animal’s immune system during a period where they experience multiple stressors that consequently increase the animal’s susceptivity to disease. Research into the bioactive compounds commonly found in marine algae is rapidly increasing due to its positive health benefits and potential immune modulating properties. Algal supplementation within previous studies has resulted in improved reproduction potential, growth performance, increases antioxidant activity and decreased proinflammatory cytokine concentrations. Additional research is required to further understand the aetiology of BRD and complete analysis of the bioavailability of these bioactive compounds within marine algae to fully explore the potential of marine algae supplementation.

Integrated Network Analysis to Identify Key Modules and Potential Hub Genes Involved in Bovine Respiratory Disease: A Systems Biology Approach

Background: Bovine respiratory disease (BRD) is the most common disease in the beef and dairy cattle industry. BRD is a multifactorial disease resulting from the interaction between environmental stressors and infectious agents. However, the molecular mechanisms underlying BRD are not fully understood yet. Therefore, this study aimed to use a systems biology approach to systematically evaluate this disorder to better understand the molecular mechanisms responsible for BRD. Methods: Previously published RNA-seq data from whole blood of 18 healthy and 25 BRD samples were downloaded from the Gene Expression Omnibus (GEO) and then analyzed. Next, two distinct methods of weighted gene coexpression network analysis (WGCNA), i.e., module-trait relationships (MTRs) and module preservation (MP) analysis were used to identify significant highly correlated modules with clinical traits of BRD and non-preserved modules between healthy and BRD samples, respectively. After identifying respective modules by the two mentioned methods of WGCNA, functional enrichment analysis was performed to extract the modules that are biologically related to BRD. Gene coexpression networks based on the hub genes from the candidate modules were then integrated with protein-protein interaction (PPI) networks to identify hub-hub genes and potential transcription factors (TFs). Results: Four significant highly correlated modules with clinical traits of BRD as well as 29 non-preserved modules were identified by MTRs and MP methods, respectively. Among them, two significant highly correlated modules (identified by MTRs) and six nonpreserved modules (identified by MP) were biologically associated with immune response, pulmonary inflammation, and pathogenesis of BRD. After aggregation of gene coexpression networks based on the hub genes with PPI networks, a total of 307 hub-hub genes were identified in the eight candidate modules. Interestingly, most of these hub-hub genes were reported to play an important role in the immune response and BRD pathogenesis. Among the eight candidate modules, the turquoise (identified by MTRs) and purple (identified by MP) modules were highly biologically enriched in BRD. Moreover, STAT1, STAT2, STAT3, IRF7, and IRF9 TFs were suggested to play an important role in the immune system during BRD by regulating the coexpressed genes of these modules. Additionally, a gene set containing several hub-hub genes was identified in the eight candidate modules, such as TLR2, TLR4, IL10, SOCS3, GZMB, ANXA1, ANXA5, PTEN, SGK1, IFI6, ISG15, MX1, MX2, OAS2, IFIH1, DDX58, DHX58, RSAD2, IFI44, IFI44L, EIF2AK2, ISG20, IFIT5, IFITM3, OAS1Y, HERC5, and PRF1, which are potentially critical during infection with agents of bovine respiratory disease complex (BRDC). Conclusion: This study not only helps us to better understand the molecular mechanisms responsible for BRD but also suggested eight candidate modules along with several promising hub-hub genes as diagnosis biomarkers and therapeutic targets for BRD.

Low serum vitamin D concentrations in Spring-born dairy calves are associated with elevated peripheral leukocytes

A role for vitamin D in the immune system is emerging from human research but data in the bovine is limited. In the current study, 48 Holstein–Friesian calves were randomly assigned to one of 4 groups designed to expose calves to divergent vitamin D levels for a 7 month period and to determine its effects on circulating immunity in young calves. Concentrations of circulating 25-hydroxyvitamin D (25OHD) was measured in serum using a commercial ELISA with validated bovine standards. Results showed that mean circulating concentrations of 25OHD at birth was 7.64 ± 3.21 ng/ml indicating vitamin D deficiency. Neither the injection of Vit D₃ at birth nor the elevated levels in milk replacer yield discernible changes to pre-weaning circulating concentration of 25OHD. No calf reached the recommended level of vitamin D immune sufficiencyof 30 ng/ml of 25OHD until at least 3 months of age (T4). Increasing dietary Vit D₃ via ration in the post-weaning period significantly elevated 25OHD concentrations in serum in VitD-In calves. Maximal levels of circulating 25OHD were achieved in VitD-Out calves, reaching 60.86 ± 7.32 ng/ml at 5 months of age (T7). Greatest divergence in haematology profile was observed between Ctl-In vs VitD-In groups with Ctl-In calves showing an elevated count of neutrophils, eosinophils, and basophils associated with reduced 25OHD concentrations. Neither IL-8 expression nor ROS production in serum were significantly different between calves with high and low 25OHD, indicating that other vitamin D-dependent mechanisms may contribute to the divergent circulating cellular profiles observed. This novel data on the vitamin D status of neonatal calves identifies a significant window of vitamin D insufficiency which is associated with significant differences in circulating immune cell profiles. Vitamin D insufficiency may therefore exacerbate pre-weaning disease susceptibility, and further work in now warranted.

Effect of treatment of pneumonia and otitis media with tildipirosin or florfenicol + flunixin meglumine on health and upper respiratory tract microbiota of preweaned Holstein dairy heifers

The objective of this randomized clinical study was to compare the effect of 2 antimicrobial interventions, tildipirosin or florfenicol + flunixin meglumine, used for treatment of pneumonia and extralabel treatment for otitis on health parameters and upper respiratory tract (URT) microbiota of preweaned Holstein calves. Housed preweaned Holstein heifers diagnosed with either otitis or pneumonia were assigned into 1 of 2 treatment groups, receiving a single subcutaneous injection of either 4 mg/kg of tildipirosin (TLD; n = 444) or 40 mg/kg of florfenicol combined with 2.2 mg/kg of a nonsteroidal anti-inflammatory, flunixin meglumine (FLF; n = 442). Calves were enrolled and treated on the day of diagnosis of the first case of pneumonia or otitis. If a calf had a recurrent case, the opposite drug was administered, respecting an interval of 5 d between drug injections. Blood samples for leukocyte counts were collected at 0, 2, 4, and 6 d after treatment, and rectal temperature was measured daily during the 5 d after treatment. Ear scores were observed from calves with otitis. Additionally, swabs of the URT were collected from a subset of 20 calves in each treatment group at d 0, 3, 6, 9, and 11 following enrollment for analysis of URT microbiota through next-generation sequencing of the 16S rRNA gene and quantitative PCR. Swabs were also collected from a comparative group of 20 healthy calves that did not receive any drug. No differences were observed between groups for recurrence risk of either pneumonia (TLD = 32.4%; FLF = 29.7%) or otitis (TLD = 72.7%; FLF = 73.6%). Similarly, no differences were observed for the total number of treatments for pneumonia (TLD = 1.45; FLF = 1.42) or otitis (TLD = 2.96; FLF = 3.07). On the other hand, both drugs reduced rectal temperature, ear scores, and leukocyte counts, with FLF calves having a greater reduction in rectal temperature within 4 d after treatment. Both TLD and FLF reduced the total bacterial load when compared with healthy untreated calves, but no differences were observed between treatment groups. Furthermore, compared with the untreated group, treated calves had lower mean relative abundances (MRA) of the genera Mannheimia, Moraxella, and Pasteurella within 11, 9, and 3 d after treatment, respectively; however, no significant differences were observed between TLD and FLF. On the other hand, MRA of Mycoplasma was not decreased by both treatments compared to untreated animals, and a higher MRA was observed in the TLD group during 11 d after treatment in comparison to FLF and untreated calves. Based on this data, we concluded that both drugs used in the study were effective in reducing rectal temperature, ear scores, leukocyte counts, and MRA of the genera Mannheimia, Pasteurella, and Moraxella in the URT, and calves treated with FLF had a greater reduction in rectal temperature.

Whole-Blood Transcriptome Analysis of Feedlot Cattle With and Without Bovine Respiratory Disease

Bovine respiratory disease (BRD) is one of the main factors leading to morbidity and mortality in feedlot operations in North America. A complex of viral and bacterial pathogens can individually or collectively establish BRD in cattle, and to date, most disease characterization studies using transcriptomic techniques examine bronchoalveolar and transtracheal fluids, lymph node, and lung tissue as well as nasopharyngeal swabs, with limited studies investigating the whole-blood transcriptome. Here, we aimed to identify differentially expressed (DE) genes involved in the host immune response to BRD using whole blood and RNA sequencing. Samples were collected from heifers (average arrival weight = 215.0 ± 5.3 kg) with (n = 25) and without (n = 18) BRD at a commercial feedlot in Western Canada. RNAseq analysis showed a distinct whole-blood transcriptome profile between BRD and non-BRD heifers. Further examination of the DE genes revealed that those involved in the host inflammatory response and infectious disease pathways were enriched in the BRD animals, while gene networks associated with metabolism and cell growth and maintenance were downregulated. Overall, the transcriptome profile derived from whole blood provided evidence that a distinct antimicrobial peptide-driven host immune response was occurring in the animals with BRD. The blood transcriptome of the BRD animals shows similarities to the transcriptome profiles obtained from lung and bronchial lymph nodes in other studies. This suggests that the blood transcriptome is a potential diagnostic tool for the identification of biomarkers of BRD infection and can be measured in live animals and used to further understand infection and disease in cattle. It may also provide a useful tool to increase the understanding of the genes involved in establishing BRD in beef cattle and be used to investigate potential therapeutic applications.

Growth performance and hematological changes of weaned beef calves diagnosed with respiratory disease using respiratory scoring and thoracic ultrasonography

This study investigated 1) the effect of clinical bovine respiratory disease (BRD) and associated lung consolidations on growth performance and hematological profiles of recently weaned beef calves and 2) the relationship between clinical respiratory signs and lung consolidation detected by thoracic ultrasonography (TUS). One hundred and fifty-three weaned beef calves (209 days old [SD: 35.8] and 306 kg [SD: 26.3], at arrival) purchased and transported from auction markets were accommodated indoors in concrete slatted floor pens. Calves were weighed weekly from arrival until day 28 and on day 65 post-arrival. Assessment of BRD and blood sample collection for hematological profiles were performed on scheduled days (at arrival, on days 7, 14, and 28) and on other days upon BRD diagnosis. Animals were assessed for BRD using a total clinical respiratory score (CRS) of five clinical signs (rectal temperature, ear position, cough, nasal secretion, and eye secretion with each ranging from normal [0] to abnormal [3]) and TUS scores (normal [0] to lung consolidation ≥ 1 cm2 [2]). Based on CRS, 35% of calves were CRS+ (CRS ≥ 5) and 65% were CRS- (CRS < 5). Although no lung consolidations (TUS-) were detected at arrival, 34% of calves developed lung consolidation (≥1 cm2) (TUS+) during the first 28 d post-arrival. Only fever (>39.6 °C) and nasal discharge were weakly associated (r = 0.19, P <0.05) with lung consolidation. On the day of BRD detection, neutrophil number and neutrophil:lymphocyte ratio were 58% and 73% greater, respectively, in BRD calves with lung consolidation compared with healthy calves. From day 0 to 65, calf average daily gain (ADG) did not differ (P >0.05) between CRS+ and CRS- calves but was 0.09 kg/d lower (P < 0.05) for TUS+ compared with TUS- calves. Calves classified as BRD (CRS + TUS ≥ 5) with lung consolidation had lower (P < 0.05) ADG from arrival until day 28 than healthy calves and BRD calves without lung consolidation (0.11 ± 0.10 vs. 0.53 ± 0.07 vs. 0.57 ± 0.10 kg/d, respectively); however, no differences in ADG were observed from day 0 to 65. Conventional methods to diagnose BRD failed to detect calves with lung lesions. TUS is a useful tool to detect lung lesions and its implementation in combination with CRS should provide a more accurate and early diagnosis of BRD, which is fundamental to successful treatment, animal welfare, and growth performance.

Influence of a Polyherbal Mixture in Dairy Calves: Growth Performance and Gene Expression

A polyherbal feed mixture containing (Achyrantes aspera, Trachyspermum ammi, Citrullus colocynthis, Andrographis paniculata, and Azadirachta indica) was evaluated in growing calves through blood chemistry, blood biometry, and gene expression during the pre-ruminant to weaning period. Forty Holstein calves (initial BW 45.6 ± 3.2 kg; 22.8 ± 0.9 days post birth) from a dairy farm were randomly assigned to the following treatments: 0, 3, 4, and 5 g/d of a polyherbal mixture, dosed in colloid gels with gelatin. Calves were housed in individual outdoor boxes with ad libitum access to a 21.5% CP calf starter and water and fed individually with a mixture of milk and a non-medicated milk replacer (22% CP). Blood samples were collected on day 59 for blood chemistry, blood biometry, and gene expression analysis in leukocyte through microarray assays. Immunoglobulins were quantified by enzyme-linked immunosorbent assay. The animals treated with the polyherbal mixture showed a quadratic effect on final body weight, daily weight gain, final hip height, and final thoracic girth. The best performance results were obtained with a treatment dose of 4 g/d. The serum IgG increased linearly with the treatment doses. Gene set enrichment analysis of upregulated genes revealed that the three biological processes with higher fold change were tight junction, mucin type O-Glycan biosynthesis, and intestinal immune network for IgA production. Also, these upregulated genes influenced arachidonic acid metabolism, and pantothenate and CoA biosynthesis. Gene ontology enrichment analysis indicated that the pathways enriched were PELP1 estrogen receptor interacting protein pathways, nuclear receptors in lipid metabolism and toxicity, tight junction, ECM-receptor interaction, thyroid hormone signaling pathways, vascular smooth muscle contraction, ribosome function, glutamatergic synapse pathway, focal adhesion, Hippo, calcium, and MAPK signaling pathways.

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 (10^3.5 TCID₅₀/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.