RNA Sequencing (RNA-Seq) Based Transcriptome Analysis in Immune Response of Holstein Cattle to Killed Vaccine against Bovine Viral Diarrhea Virus Type I

Veterinary systems biology for bridging the phenotype-genotype gap via computational modeling for disease epidemiology and animal welfare

Veterinary systems biology is an innovative approach that integrates biological data at the molecular and cellular levels, allowing for a more extensive understanding of the interactions and functions of complex biological systems in livestock and veterinary science. It has tremendous potential to integrate multi-omics data with the support of vetinformatics resources for bridging the phenotype-genotype gap via computational modeling. To understand the dynamic behaviors of complex systems, computational models are frequently used. It facilitates a comprehensive understanding of how a host system defends itself against a pathogen attack or operates when the pathogen compromises the host's immune system. In this context, various approaches, such as systems immunology, network pharmacology, vaccinology and immunoinformatics, can be employed to effectively investigate vaccines and drugs. By utilizing this approach, we can ensure the health of livestock. This is beneficial not only for animal welfare but also for human health and environmental well-being. Therefore, the current review offers a detailed summary of systems biology advancements utilized in veterinary sciences, demonstrating the potential of the holistic approach in disease epidemiology, animal welfare and productivity.

RNA-Seq Analysis of Peripheral Whole Blood from Dairy Bulls with High and Low Antibody-Mediated Immune Responses-A Preliminary Study

Enhancing the immune response through breeding is regarded as an effective strategy for improving animal health, as dairy cattle identified as high immune responders are reported to have a decreased prevalence of economically significant diseases. The identification of differentially expressed genes (DEGs) associated with immune responses might be an effective tool for breeding healthy dairy cattle. In this study, antibody-mediated immune responses (AMIRs) were induced by the immunization of hen egg white lysozyme (HEWL) in six Chinese Holstein dairy bulls divided into high- and low-AMIR groups based on their HEWL antibody level. Then, RNA-seq was applied to explore the transcriptome of peripheral whole blood between the two comparison groups. As a result, several major upregulated and downregulated genes were identified and attributed to the regulation of locomotion, tissue development, immune response, and detoxification. In addition, the result of the KEGG pathway analysis revealed that most DEGs were enriched in pathways related to disease, inflammation, and immune response, including antigen processing and presentation, Staphylococcus aureus infection, intestinal immune network for IgA production, cytokine-cytokine receptor interaction, and complement and coagulation cascades. Moreover, six genes (BOLA-DQA5, C5, CXCL2, HBA, LTF, and COL1A1) were validated using RT-qPCR, which may provide information for genomic selection in breeding programs. These results broaden the knowledge of the immune response mechanism in dairy bulls, which has strong implications for breeding cattle with an enhanced AMIR.

Raw potato starch diet supplement in weaned pigs could reduce Salmonella Typhimurium infection by altering microbiome composition and improving immune status

Backgorund

Salmonella enterica serovar Typhimurium (ST) is one of the causative agents of gastroenteritis in pigs. Pigs fed a diet supplemented with raw potato starch (RPS) have improved gut health by the alteration of the microbiota composition and production of short-chain fatty acids (SCFAs). This study aimed to evaluate the effects of RPS supplementation in reducing infection severity and fecal shedding in ST-infected pigs.

Methods

The weaned experimental pigs were divided into two groups: CON (n = 6) fed a corn/soybean-based diet and TRT (n = 6) supplemented with 5% RPS. After 21 d, the pigs were inoculated with ST, and their body weight, clinical signs, and fecal shedding of ST were monitored for 14 d. At 14 d post-inoculation (dpi), the jejunum, cecum, ileum, and colon tissues were collected from euthanized pigs, and histopathological lesions and cytokine gene expression were compared. Additionally, blood samples at 2 dpi were analyzed for gene ontology enrichment. Moreover, the gutmicrobiome was analyzed using 16S rRNA metagenomic sequencing, and the SCFA concentration was measured using gas chromatography.

Results

The average daily weight gain was significantly higher in TRT than in CON during the ST infection period; however, histopathological lesion scores were significantly lower in TRT than in CON. The relative abundance of nine genera of butyrate- and acetate-producing bacteria significantly increased in TRT compared with that of only two acetate-producing bacteria in CON. Among the genes involved in the immune response, IL-18 expression level was significantly lower in the jejunum and colon in TRT than in CON. Furthermore, Reg3γ expression was significantly different in the cecum and colon of both groups.

Conclusion

The diet supplemented with RPS in weaned pigs could result in predominance of butyrate- and acetate-producing bacteria, reducing the severity of ST infection by improving the immune status.

Comparative transcriptome analysis reveals the immune response of turbot (Scophthalmus maximus) induced by inactivated bivalent vaccine

Vibrio species are important pathogens that affect a wide range of farmed fish. Vaccination is regarded as the most efficacious strategy for fighting bacterial infections. However, the underlying mechanisms remain to be elucidated. In the present study, a comparative transcriptome analysis was performed on the spleens from turbot (Scophthalmus maximus) induced by an inactivated bivalent vaccine (Vibrio anguillarum and Vibrio harveyi, IVVah1) at 4 week and 1 day post further challenge. Strong immune responses were induced by the bivalent vaccine, besides differentially expressed genes (DEGs) associated with adaptive immunity, more innate immunity-related DEGs were detected. At the late stage of vaccination, immune-related molecules associated with pattern recognition receptors, inflammatory factors, complement and coagulation cascade-related components, and antigen processing and presentation were significantly regulated, and some of them were even further up-regulated after the bacterial challenge, indicating the cooperation of multiple immune processes during the vaccine immunization process. In addition to the terms or pathways associated with the immune response, enrichment analysis revealed multiple significantly enriched terms/pathways associated with the response to stimulus/stress, homeostasis, metabolism, and biosynthesis, suggesting that a defensive status was established by the bivalent vaccine. This study furnishes new insights into the internal mechanism of immunity upon a combined vaccine administrating in turbot and lays a foundation for developing highly immunogenic vaccines in teleost.

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.

Non-structural proteins of bovine viral diarrhea virus

Bovine viral diarrhea virus (BVDV) belongs to the family Flaviviridae genus pestivirus. The viral genome is a single-stranded, positive-sense RNA that encodes four structural proteins (i.e., C, Erns, E1, and E2) and eight non-structural proteins (NSPs) (i.e., Npro, p7, NS2, NS3, NS4A, NS4B, NS5A, and NS5B). Cattle infected with BVDV exhibit a number of different clinical signs including diarrhea, abortion, and other reproductive disorders which have a serious impact on the cattle industry worldwide. Research on BVDV mainly focuses on its structural protein, however, progress in understanding the functions of the NSPs of BVDV has also been made in recent decades. The knowledge gained on the BVDV non-structural proteins is helpful to more fully understand the viral replication process and the molecular mechanism of viral persistent infection. This review focuses on the functions of BVDV NSPs and provides references for the identification of BVDV, the diagnosis and prevention of Bovine viral diarrhea mucosal disease (BVD-MD), and the development of vaccines.

Th17 cell differentiation induced by cytopathogenic biotype BVDV-2 in bovine PBLCs

Background

Bovine viral diarrhea virus (BVDV) is a major pathogen that causes bovine viral diarrhea/mucosal disease (BVD-MD), which has become a global infectious disease due to its wide spread and the lack of effective treatment. The process of BVDV infection is complex. Once infected, host immune cells are activated and modulated. As a major immune cell, peripheral blood lymphocyte cells (PBLCs) are the primary target of BVDV. In order to further understand the mechanism of BVDV- host interaction, the expression profiles of host lymphocytes mRNAs associated with BVDV infection were investigated by transcriptomic sequencing analysis.

Results

The transcriptomic sequencing analysis was performed on bovine PBLCs infected with CP BVDV-2 GS2018 after 12 h of infection. Gene expression profiling demonstrated that 1052 genes were differentially expressed in GS2018 infected PBLCs compared with the control group. Of these genes, 485 genes were up-regulated and 567 were down-regulated. The 19 differential expressed genes (DEGs) were selected for validation using quantitative real-time PCR and the results were consistent with the results of RNA-Seq. Gene ontology enrichment and KEGG pathway analysis showed that 1052 DEGs were significantly enriched in 16 pathways, including cytokine-cytokine receptor interaction, IL17, PI3K-Akt, MAPK and TNF signaling pathway. PPI network analysis showed that IL17A, IFN-γ and TNF-α interacted with various proteins and may play crucial roles in BVDV-2 infection. Of note, we confirmed that GS2018 induced Th17 cell differentiation in PBLCs and persistently increased the expression levels of IL17A. In turn, the replication of GS2018 was inhibited by IL17A.

Conclusion

In this study, the transcription changes of DEGs related to host immune responses in bovine PBLCs were caused by CP BVDV-2 infection. In particular, the effector molecules IL17A of Th17 cells were significantly up-regulated, which inhibited viral replication. These results will contribute to exploration and further understanding of the host immune response mechanism and interaction between host and BVDV-2.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08194-w.

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.

Post-Translational Protein Deimination Signatures in Plasma and Plasma EVs of Reindeer (Rangifer tarandus)

Simple Summary

Reindeer are an important wild and domesticated species of the Arctic, Northern Europe, Siberia and North America. As reindeer have developed various strategies to adapt to extreme environments, this makes them an interesting species for studies into diversity of immune and metabolic functions in the animal kingdom. Importantly, while reindeer carry natural infections caused by viruses (including coronaviruses), bacteria and parasites, they can also act as carriers for transmitting such diseases to other animals and humans, so called zoonosis. Reindeer are also affected by chronic wasting disease, a neuronal disease caused by prions, similar to scrapie in sheep, mad cows disease in cattle and Creutzfeldt-Jakob disease in humans. The current study assessed a specific protein modification called deimination/citrullination, which can change how proteins function and allow them to take on different roles in health and disease processes. Profiling of deiminated proteins in reindeer showed that many important pathways for immune defenses, prion diseases and metabolism are enriched in deiminated proteins, both in plasma, as well as in plasma extracellular vesicles. This study provides a platform for the development of novel biomarkers to assess wild life health status and factors relating to zoonotic disease.

Abstract

The reindeer (caribou) Rangifer tarandus is a Cervidae in the order Artiodactyla. Reindeer are sedentary and migratory populations with circumpolar distribution in the Arctic, Northern Europe, Siberia and North America. Reindeer are an important wild and domesticated species, and have developed various adaptive strategies to extreme environments. Importantly, deer have also been identified to be putative zoonotic carriers, including for parasites, prions and coronavirus. Therefore, novel insights into immune-related markers are of considerable interest. Peptidylarginine deiminases (PADs) are a phylogenetically conserved enzyme family which causes post-translational protein deimination by converting arginine into citrulline in target proteins. This affects protein function in health and disease. Extracellular vesicles (EVs) participate in cellular communication, in physiological and pathological processes, via transfer of cargo material, and their release is partly regulated by PADs. This study assessed deiminated protein and EV profile signatures in plasma from sixteen healthy wild female reindeer, collected in Iceland during screening for parasites and chronic wasting disease. Reindeer plasma EV profiles showed a poly-dispersed distribution from 30 to 400 nm and were positive for phylogenetically conserved EV-specific markers. Deiminated proteins were isolated from whole plasma and plasma EVs, identified by proteomic analysis and protein interaction networks assessed by KEGG and GO analysis. This revealed a large number of deimination-enriched pathways for immunity and metabolism, with some differences between whole plasma and EVs. While shared KEGG pathways for whole plasma and plasma EVs included complement and coagulation pathways, KEGG pathways specific for EVs were for protein digestion and absorption, platelet activation, amoebiasis, the AGE–RAGE signaling pathway in diabetic complications, ECM receptor interaction, the relaxin signaling pathway and the estrogen signaling pathway. KEGG pathways specific for whole plasma were pertussis, ferroptosis, SLE, thyroid hormone synthesis, phagosome, Staphylococcus aureus infection, vitamin digestion and absorption, and prion disease. Further differences were also found between molecular function and biological processes GO pathways when comparing functional STRING networks for deiminated proteins in EVs, compared with deiminated proteins in whole plasma. This study highlights deiminated proteins and EVs as candidate biomarkers for reindeer health and may provide information on regulation of immune pathways in physiological and pathological processes, including neurodegenerative (prion) disease and zoonosis.

Gene signatures of SARS-CoV/SARS-CoV-2-infected ferret lungs in short- and long-term models

Coronaviruses (CoVs) consist of six strains, and the severe acute respiratory syndrome coronavirus (SARS-CoV), newly found coronavirus (SARS-CoV-2) has rapidly spread leading to a global outbreak. The ferret (Mustela putorius furo) serves as a useful animal model for studying SARS-CoV/SARS-CoV-2 infection and developing therapeutic strategies. A holistic approach for distinguishing differences in gene signatures during disease progression is lacking. The present study discovered gene expression profiles of short-term (3 days) and long-term (14 days) ferret models after SARS-CoV/SARS-CoV-2 infection using a bioinformatics approach. Through Gene Ontology (GO) and MetaCore analyses, we found that the development of stemness signaling was related to short-term SARS-CoV/SARS-CoV-2 infection. In contrast, pathways involving extracellular matrix and immune responses were associated with long-term SARS-CoV/SARS-CoV-2 infection. Some highly expressed genes in both short- and long-term models played a crucial role in the progression of SARS-CoV/SARS-CoV-2 infection, including DPP4, BMP2, NFIA, AXIN2, DAAM1, ZNF608, ME1, MGLL, LGR4, ABHD6, and ACADM. Meanwhile, we revealed that metabolic, glucocorticoid, and reactive oxygen species-associated networks were enriched in both short- and long-term infection models. The present study showed alterations in gene expressions from short-term to long-term SARS-CoV/SARS-CoV-2 infection. The current result provides an explanation of the pathophysiology for post-infectious sequelae and potential targets for treatment.

The Effect of Bovine Viral Diarrhea Virus (BVDV) Strains and the Corresponding Infected-Macrophages' Supernatant on Macrophage Inflammatory Function and Lymphocyte Apoptosis

Bovine viral diarrhea virus (BVDV) is an important viral disease of cattle that causes immune dysfunction. Macrophages are the key cells for the initiation of the innate immunity and play an important role in viral pathogenesis. In this in vitro study, we studied the effect of the supernatant of BVDV-infected macrophage on immune dysfunction. We infected bovine monocyte-derived macrophages (MDM) with high or low virulence strains of BVDV. The supernatant recovered from BVDV-infected MDM was used to examine the functional activity and surface marker expression of normal macrophages as well as lymphocyte apoptosis. Supernatants from the highly virulent 1373-infected MDM reduced phagocytosis, bactericidal activity and downregulated MHC II and CD14 expression of macrophages. Supernatants from 1373-infected MDM induced apoptosis in MDBK cells, lymphocytes or BL-3 cells. By protein electrophoresis, several protein bands were unique for high-virulence, 1373-infected MDM supernatant. There was no significant difference in the apoptosis-related cytokine mRNA (IL-1beta, IL-6 and TNF-a) of infected MDM. These data suggest that BVDV has an indirect negative effect on macrophage functions that is strain-specific. Further studies are required to determine the identity and mechanism of action of these virulence factors present in the supernatant of the infected macrophages.

Post-Translational Protein Deimination Signatures in Serum and Serum-Extracellular Vesicles of Bos taurus Reveal Immune, Anti-Pathogenic, Anti-Viral, Metabolic and Cancer-Related Pathways for Deimination

The bovine immune system is known for its unusual traits relating to immunoglobulin and antiviral responses. Peptidylarginine deiminases (PADs) are phylogenetically conserved enzymes that cause post-translational deimination, contributing to protein moonlighting in health and disease. PADs also regulate extracellular vesicle (EV) release, forming a critical part of cellular communication. As PAD-mediated mechanisms in bovine immunology and physiology remain to be investigated, this study profiled deimination signatures in serum and serum-EVs in Bos taurus. Bos EVs were poly-dispersed in a 70-500 nm size range and showed differences in deiminated protein cargo, compared with whole sera. Key immune, metabolic and gene regulatory proteins were identified to be post-translationally deiminated with some overlapping hits in sera and EVs (e.g., immunoglobulins), while some were unique to either serum or serum-EVs (e.g., histones). Protein-protein interaction network analysis of deiminated proteins revealed KEGG pathways common for serum and serum-EVs, including complement and coagulation cascades, viral infection (enveloped viruses), viral myocarditis, bacterial and parasitic infections, autoimmune disease, immunodeficiency intestinal IgA production, B-cell receptor signalling, natural killer cell mediated cytotoxicity, platelet activation and hematopoiesis, alongside metabolic pathways including ferroptosis, vitamin digestion and absorption, cholesterol metabolism and mineral absorption. KEGG pathways specific to EVs related to HIF-1 signalling, oestrogen signalling and biosynthesis of amino acids. KEGG pathways specific for serum only, related to Epstein-Barr virus infection, transcription mis-regulation in cancer, bladder cancer, Rap1 signalling pathway, calcium signalling pathway and ECM-receptor interaction. This indicates differences in physiological and pathological pathways for deiminated proteins in serum-EVs, compared with serum. Our findings may shed light on pathways underlying a number of pathological and anti-pathogenic (viral, bacterial, parasitic) pathways, with putative translatable value to human pathologies, zoonotic diseases and development of therapies for infections, including anti-viral therapies.