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

Integrated Analysis of Transcriptome Profiles and lncRNA-miRNA-mRNA Competing Endogenous RNA Regulatory Network to Identify Biological Functional Effects of Genes and Pathways Associated with Johne's Disease in Dairy Cattle

Paratuberculosis or Johne's disease (JD), a chronic granulomatous gastroenteritis caused by Mycobacterium avium subsp. paratuberculosis (MAP), causes huge economic losses and reduces animal welfare in dairy cattle herds worldwide. At present, molecular mechanisms and biological functions involved in immune responses to MAP infection of dairy cattle are not clearly understood. Our purpose was to integrate transcriptomic profiles and competing endogenous RNA (ceRNA) network analyses to identify key messenger RNAs (mRNAs) and regulatory RNAs involved in molecular regulation of peripheral blood mononuclear cells (PBMCs) for MAP infection in dairy cattle. In total, 28 lncRNAs, 42 miRNAs, and 370 mRNAs were identified by integrating gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses. In this regard, we identified 21 hub genes (CCL20, CCL5, CD40, CSF2, CXCL8, EIF2AK2, FOS, IL10, IL17A, IL1A, IL1B, IRF1, MX2, NFKB1, NFKBIA, PTGS2, SOCS3, TLR4, TNF, TNFAIP3, and VCAM1) involved in MAP infection. Furthermore, eight candidate subnets with eight lncRNAs, 29 miRNAs, and 237 mRNAs were detected through clustering analyses, whereas GO enrichment analysis of identified RNAs revealed 510, 22, and 11 significantly enriched GO terms related to MAP infection in biological process, molecular function, and cellular component categories, respectively. The main metabolic-signaling pathways related to MAP infection that were enriched included the immune system process, defense response, response to cytokine, leukocyte migration, regulation of T cell activation, defense response to bacterium, NOD-like receptor, B cell receptor, TNF, NF-kappa B, IL-17, and T cell receptor signaling pathways. Contributions of transcriptome profiles from MAP-positive and MAP-negative sample groups plus a ceRNA regulatory network underlying phenotypic differences in the intensity of pathogenicity of JD provided novel insights into molecular mechanisms associated with immune system responses to MAP infection in dairy cattle.

Strategies for Bovine Respiratory Disease (BRD) Diagnosis and Prognosis: A Comprehensive Overview

Despite significant advances in vaccination strategies and antibiotic therapy, bovine respiratory disease (BRD) continues to be the leading disease affecting the global cattle industry. The etiology of BRD is complex, often involving multiple microbial agents, which lead to intricate interactions between the host immune system and pathogens during various beef production stages. These interactions present environmental, social, and geographical challenges. Accurate diagnosis is essential for effective disease management. Nevertheless, correct identification of BRD cases remains a daunting challenge for animal health technicians in feedlots. In response to current regulations, there is a growing interest in refining clinical diagnoses of BRD to curb the overuse of antimicrobials. This shift marks a pivotal first step toward establishing a structured diagnostic framework for this disease. This review article provides an update on recent developments and future perspectives in clinical diagnostics and prognostic techniques for BRD, assessing their benefits and limitations. The methods discussed include the evaluation of clinical signs and animal behavior, biomarker analysis, molecular diagnostics, ultrasound imaging, and prognostic modeling. While some techniques show promise as standalone diagnostics, it is likely that a multifaceted approach-leveraging a combination of these methods-will yield the most accurate diagnosis of BRD.

Modeling Human Respiratory Syncytial Virus (RSV) Infection: Recent Contributions and Future Directions Using the Calf Model of Bovine RSV Disease

The human orthopneumovirus (human respiratory syncytial virus [RSV]) is a leading cause of respiratory disease in children worldwide and a significant cause of infant mortality in low- and middle-income countries. The natural immune response to the virus has a preponderant role in disease progression, with a rapid neutrophil infiltration and dysbalanced T cell response in the lungs associated with severe disease in infants. The development of preventive interventions against human RSV has been difficult partly due to the need to use animal models that only partially recapitulate the immune response as well as the disease progression seen in human infants. In this brief review, we discuss the contributions of the calf model of RSV infection to understanding immunity to RSV and in developing vaccine and drug candidates, focusing on recent research areas. We propose that the bovine model of RSV infection is a valuable alternative for assessing the translational potential of interventions aimed at the human population.

Integrated Comparative Transcriptome and circRNA-lncRNA-miRNA-mRNA ceRNA Regulatory Network Analyses Identify Molecular Mechanisms Associated with Intramuscular Fat Content in Beef Cattle

Intramuscular fat content (IMF), one of the most important carcass traits in beef cattle, is controlled by complex regulatory factors. At present, molecular mechanisms involved in regulating IMF and fat metabolism in beef cattle are not well understood. Our objective was to integrate comparative transcriptomic and competing endogenous RNA (ceRNA) network analyses to identify candidate messenger RNAs (mRNAs) and regulatory RNAs involved in molecular regulation of longissimus dorsi muscle (LDM) tissue for IMF and fat metabolism of 5 beef cattle breeds (Angus, Chinese Simmental, Luxi, Nanyang, and Shandong Black). In total, 34 circRNAs, 57 lncRNAs, 15 miRNAs, and 374 mRNAs were identified by integrating gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses. Furthermore, 7 key subnets with 16 circRNAs, 43 lncRNAs, 7 miRNAs, and 237 mRNAs were detected through clustering analyses, whereas GO enrichment analysis of identified RNAs revealed 48, 13, and 28 significantly enriched GO terms related to IMF in biological process, molecular function, and cellular component categories, respectively. The main metabolic-signaling pathways associated with IMF and fat metabolism that were enriched included metabolic, calcium, cGMP-PKG, thyroid hormone, and oxytocin signaling pathways. Moreover, MCU, CYB5R1, and BAG3 genes were common among the 10 comparative groups defined as important candidate marker genes for fat metabolism in beef cattle. Contributions of transcriptome profiles from various beef breeds and a competing endogenous RNA (ceRNA) regulatory network underlying phenotypic differences in IMF provided novel insights into molecular mechanisms associated with meat quality.

Influence of the At-Arrival Host Transcriptome on Bovine Respiratory Disease Incidence during Backgrounding

Bovine respiratory disease (BRD) remains the leading disease within the U.S. beef cattle industry. Marketing decisions made prior to backgrounding may shift BRD incidence into a different phase of production, and the importance of host gene expression on BRD incidence as it relates to marketing strategy is poorly understood. Our objective was to compare the influence of marketing on host transcriptomes measured on arrival at a backgrounding facility on the subsequent probability of being treated for BRD during a 45-day backgrounding phase. This study, through RNA-Seq analysis of blood samples collected on arrival, evaluated gene expression differences between cattle which experienced a commercial auction setting (AUCTION) versus cattle directly shipped to backgrounding from the cow–calf phase (DIRECT); further analyses were conducted to determine differentially expressed genes (DEGs) between cattle which remained clinically healthy during backgrounding (HEALTHY) versus those that required treatment for clinical BRD within 45 days of arrival (BRD). A profound difference in DEGs (n = 2961) was identified between AUCTION cattle compared to DIRECT cattle, regardless of BRD development; these DEGs encoded for proteins involved in antiviral defense (increased in AUCTION), cell growth regulation (decreased in AUCTION), and inflammatory mediation (decreased in AUCTION). Nine and four DEGs were identified between BRD and HEALTHY cohorts in the AUCTION and DIRECT groups, respectively; DEGs between disease cohorts in the AUCTION group encoded for proteins involved in collagen synthesis and platelet aggregation (increased in HEALTHY). Our work demonstrates the clear influence marketing has on host expression and identified genes and mechanisms which may predict BRD risk.

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.

Impact of preweaning vaccination on host gene expression and antibody titers in healthy beef calves

The impact of preweaning vaccination for bovine respiratory viruses on cattle health and subsequent bovine respiratory disease morbidity has been widely studied yet questions remain regarding the impact of these vaccines on host response and gene expression. Six randomly selected calves were vaccinated twice preweaning (T1 and T3) with a modified live vaccine for respiratory pathogens and 6 randomly selected calves were left unvaccinated. Whole blood samples were taken at first vaccination (T1), seven days later (T2), at revaccination and castration (T3), and at weaning (T4), and utilized for RNA isolation and sequencing. Serum from T3 and T4 was analyzed for antibodies to BRSV, BVDV1a, and BHV1. Sequenced RNA for all 48 samples was bioinformatically processed with a HISAT2/StringTie pipeline, utilizing reference guided assembly with the ARS-UCD1.2 bovine genome. Differentially expressed genes were identified through analyzing the impact of time across all calves, influence of vaccination across treatment groups at each timepoint, and the interaction of time and vaccination. Calves, regardless of vaccine administration, demonstrated an increase in gene expression over time related to specialized proresolving mediator production, lipid metabolism, and stimulation of immunoregulatory T-cells. Vaccination was associated with gene expression related to natural killer cell activity and helper T-cell differentiation, enriching for an upregulation in Th17-related gene expression, and downregulated genes involved in complement system activity and coagulation mechanisms. Type-1 interferon production was unaffected by the influence of vaccination nor time. To our knowledge, this is the first study to evaluate mechanisms of vaccination and development in healthy calves through RNA sequencing analysis.

Assessment of Rapid MinION Nanopore DNA Virus Meta-Genomics Using Calves Experimentally Infected with Bovine Herpes Virus-1

Bovine respiratory disease (BRD), which is the leading cause of morbidity and mortality in cattle, is caused by numerous known and unknown viruses and is responsible for the widespread use of broad-spectrum antibiotics despite the use of polymicrobial BRD vaccines. Viral metagenomics sequencing on the portable, inexpensive Oxford Nanopore Technologies MinION sequencer and sequence analysis with its associated user-friendly point-and-click Epi2ME cloud-based pathogen identification software has the potential for point-of-care/same-day/sample-to-result metagenomic sequence diagnostics of known and unknown BRD pathogens to inform a rapid response and vaccine design. We assessed this potential using in vitro viral cell cultures and nasal swabs taken from calves that were experimentally challenged with a single known BRD-associated DNA virus, namely, bovine herpes virus 1. Extensive optimisation of the standard Oxford Nanopore library preparation protocols, particularly a reduction in the PCR bias of library amplification, was required before BoHV-1 could be identified as the main virus in the in vitro cell cultures and nasal swab samples within approximately 7 h from sample to result. In addition, we observed incorrect assignment of the bovine sequence to bacterial and viral taxa due to the presence of poor-quality bacterial and viral genome assemblies in the RefSeq database used by the EpiME Fastq WIMP pathogen identification software.

Use of nCounter mRNA profiling to identify at-arrival gene expression patterns for predicting bovine respiratory disease in beef cattle

Background

Transcriptomics has identified at-arrival differentially expressed genes associated with bovine respiratory disease (BRD) development; however, their use as prediction molecules necessitates further evaluation. Therefore, we aimed to selectively analyze and corroborate at-arrival mRNA expression from multiple independent populations of beef cattle. In a nested case-control study, we evaluated the expression of 56 mRNA molecules from at-arrival blood samples of 234 cattle across seven populations via NanoString nCounter gene expression profiling. Analysis of mRNA was performed with nSolver Advanced Analysis software (p < 0.05), comparing cattle groups based on the diagnosis of clinical BRD within 28 days of facility arrival (n = 115 Healthy; n = 119 BRD); BRD was further stratified for severity based on frequency of treatment and/or mortality (Treated_1, n = 89; Treated_2+, n = 30). Gene expression homogeneity of variance, receiver operator characteristic (ROC) curve, and decision tree analyses were performed between severity cohorts.

Results

Increased expression of mRNAs involved in specialized pro-resolving mediator synthesis (ALOX15, HPGD), leukocyte differentiation (LOC100297044, GCSAML, KLF17), and antimicrobial peptide production (CATHL3, GZMB, LTF) were identified in Healthy cattle. BRD cattle possessed increased expression of CFB, and mRNA related to granulocytic processes (DSG1, LRG1, MCF2L) and type-I interferon activity (HERC6, IFI6, ISG15, MX1). Healthy and Treated_1 cattle were similar in terms of gene expression, while Treated_2+ cattle were the most distinct. ROC cutoffs were used to generate an at-arrival treatment decision tree, which classified 90% of Treated_2+ individuals.

Conclusions

Increased expression of complement factor B, pro-inflammatory, and type I interferon-associated mRNA hallmark the at-arrival expression patterns of cattle that develop severe clinical BRD. Here, we corroborate at-arrival mRNA markers identified in previous transcriptome studies and generate a prediction model to be evaluated in future studies. Further research is necessary to evaluate these expression patterns in a prospective manner.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12917-022-03178-8.

Alterations in arthropod and neuronal exosomes reduce virus transmission and replication in recipient cells

Aim

Targeting the modes of pathogen shedding/transmission via exosomes or extracellular vesicles has been envisioned as the best approach to control vector-borne diseases. This study is focused on altering exosomes stability to affect the pathogen transmission from infected to naïve recipient cells.

Methods

In this study, neuronal or arthropod exosomes were treated at different temperatures or with different salts or pH conditions to analyze their ability and efficiency in the transmission of tick-borne Langat virus (LGTV) from infected to naïve recipient cells.

Results

Quantitative real-time PCR (qRT-PCR) and immunoblotting analyses revealed that treatment of neuronal or tick exosomes at warmer temperatures of 37 °C or 23 °C, respectively, or with sulfate salts such as Magnesium or Ammonium sulfates or with highly alkaline pH of 9 or 11.5, dramatically reduced transmission of LGTV via infectious exosomes (human or tick cells-derived) to human neuronal (SH-SY5Y) cells or skin keratinocytes (HaCaT cells), respectively.

Conclusion

Overall, this study suggests that exosome-mediated viral transmission of vector-borne pathogens to the vertebrate host or the viral dissemination and replication within or between the mammalian host can be reduced by altering the ability of exosomes with basic changes in temperatures, salts or pH conditions.

A Review of UK-Registered and Candidate Vaccines for Bovine Respiratory Disease

Vaccination is widely regarded as a cornerstone in animal or herd health and infectious disease management. Nineteen vaccines against the major pathogens implicated in bovine respiratory disease are registered for use in the UK by the Veterinary Medicines Directorate (VMD). However, despite annual prophylactic vaccination, bovine respiratory disease is still conservatively estimated to cost the UK economy approximately £80 million per annum. This review examines the vaccine types available, discusses the surrounding literature and scientific rationale of the limitations and assesses the potential of novel vaccine technologies.

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.