RNA-Seq based transcriptome analysis during bovine viral diarrhoea virus (BVDV) infection

Preclinical Efficacy of Cap-Dependent and Independent mRNA Vaccines against Bovine Viral Diarrhea Virus-1

Bovine viral diarrhea virus (BVDV) is an RNA virus associated with severe economic losses in animal production. Effective vaccination and viral surveillance are urgent for the prevention and control of BVDV infection. However, the application of traditional modified live vaccines and inactivated vaccines is faced with tremendous challenges. In the present study, we describe the preclinical efficacy of two BVDV mRNA vaccines tested in mice and guinea pigs, followed by a field trial in goats, where they were compared to a commercial vaccine (formaldehyde inactivated). The two mRNAs were engineered to express the envelope protein E2 of BVDV-1, the most prevalent subtype across the world, through a 5' cap-dependent or independent fashion. Better titers of neutralizing antibodies against BVDV-1 were achieved using the capped RNA in the sera of mice and guinea pigs, with maximum values reaching 9.4 and 13.7 (by -log2), respectively, on the 35th day post-vaccination. At the same time point, the antibody levels in goats were 9.1 and 10.2 for the capped and capless RNAs, respectively, and there were no significant differences compared to the commercial vaccine. The animals remained healthy throughout the experiment, as reflected by their normal leukogram profiles. Collectively, our findings demonstrate that mRNA vaccines have good safety and immunogenicity, and we laid a strong foundation for the further exploitation of efficient and safe BVDV vaccines.

Construction of an IFNAR1 knockout MDBK cell line using CRISPR/Cas9 and its effect on bovine virus replication

The type I interferon (IFN) pathway is important for eukaryotic cells to resist viral infection, as well as an impediment to efficient virus replication. Therefore, this study aims to create an IFNAR1 knockout (KO) Madin-Darby bovine kidney (MDBK) cell line using CRISPR/Cas9 and investigate its application and potential mechanism in increasing viral replication of bovines. The IFNAR1 KO cells showed increased titers of bovine viral diarrhea virus (BVDV) (1.5 log10), with bovine enterovirus and bovine parainfluenza virus type 3 (0.5-0.8 log10). RNA-seq revealed reduced expression of the genes related IFN-I pathways including IFNAR1, STAT3, IRF9, and SOCS3 in IFNAR1 KO cells compared with WT cells. In WT cells, 306 differentially expressed genes (DEGs) were identified between BVDV-infected and -uninfected cells. Of these, 128 up- and 178 down-regulated genes were mainly associated with growth cycle and biosynthesis, respectively. In IFNAR1 KO cells, 286 DEGs were identified, with 82 up-regulated genes were associated with signaling pathways, and 204 down-regulated genes. Further, 92 DEGs were overlapped between WT and IFNAR1 KO cells including ESM1, IL13RA2, and SLC25A34. Unique DEGs in WT cells were related to inflammation and immune regulation, whereas those unique in IFNAR1 KO cells involved in cell cycle regulation through pathways such as MAPK. Knocking down SLC25A34 and IL13RA2 in IFNAR1 KO cells increased BVDV replication by 0.3 log10 and 0.4 log10, respectively. Additionally, we constructed an IFNAR1/IFNAR2 double-knockout MDBK cell line, which further increased BVDV viral titers compared with IFNAR1 KO cells (0.6 log10). Overall, the IFNAR1 KO MDBK cell line can support better replication of bovine viruses and therefore provides a valuable tool for bovine virus research on viral pathogenesis and host innate immune response.

PI3K/AKT mediated De novo fatty acid synthesis regulates RIG-1/MDA-5-dependent type I IFN responses in BVDV-infected CD8+T cells

Bovine viral diarrhea virus (BVDV) has caused massive economic losses in the cattle business worldwide. Fatty acid synthase (FASN), a key enzyme of the fatty acid synthesis (FAS) pathway, has been shown to support virus replication. To investigate the role of fatty acids (FAs) in BVDV infection, we infected CD8+T lymphocytes obtained from healthy cattle with BVDV in vitro. During early cytopathic (CP) and noncytopathic (NCP) BVDV infection in CD8+ T cells, there is an increase in de novo lipid biosynthesis, resulting in elevated levels of free fatty acids (FFAs) and triglycerides (TG). BVDV infection promotes de novo lipid biosynthesis in a dose-dependent manner. Treatment with the FASN inhibitor C75 significantly reduces the phosphorylation of PI3K and AKT in BVDV-infected CD8+ T cells, while inhibition of PI3K with LY294002 decreases FASN expression. Both CP and NCP BVDV strains promote de novo fatty acid synthesis by activating the PI3K/AKT pathway. Further investigation shows that pharmacological inhibitors targeting FASN and PI3K concurrently reduce FFAs, TG levels, and ATP production, effectively inhibiting BVDV replication. Conversely, the in vitro supplementation of oleic acid (OA) to replace fatty acids successfully restored BVDV replication, underscoring the impact of abnormal de novo fatty acid metabolism on BVDV replication. Intriguingly, during BVDV infection of CD8+T cells, the use of FASN inhibitors prompted the production of IFN-α and IFN-β, as well as the expression of interferon-stimulated genes (ISGs). Moreover, FASN inhibitors induce TBK-1 phosphorylation through the activation of RIG-1 and MDA-5, subsequently activating IRF-3 and ultimately enhancing the IFN-1 response. In conclusion, our study demonstrates that BVDV infection activates the PI3K/AKT pathway to boost de novo fatty acid synthesis, and inhibition of FASN suppresses BVDV replication by activating the RIG-1/MDA-5-dependent IFN response.

Efficacy of H2O2 inactivated bovine virus diarrhoea virus (BVDV) type 1 vaccine in mice

BACKGROUND

Bovine viral diarrhea (BVD) is an acute febrile infectious disease caused by the bovine viral diarrhea virus (BVDV), which has brought huge economic losses to the world's cattle industry. At present, commercial inactivated BVDV vaccines may cause some adverse reactions during use. This study aims to develop a safer and more efficient inactivated BVDV vaccine.

METHODS

Here, we described the generation and preclinical efficacy of a hydrogen peroxide (H2O2) inactivated BVDV type 1 vaccine in mice, and administered it separately with commercial vaccine (formaldehyde inactivated) in mice to study its efficacy.

RESULTS

The BVDV type 1 IgG, IFN- γ, IL-4 and neutralizing antibody in the serum of the H2O2 inactivated vaccine group can be maintained in mice for 70 days. The IgG level reached its maximum value of 0.67 on the 42nd day, significantly higher than the commercial formaldehyde inactivated BVDV type 1 vaccine. IFN- γ and IL-4 reached their maximum values on the 28th day after immunization, at 123.16 pg/ml and 143.80 pg/ml, respectively, slightly higher than commercial vaccines, but the effect was not significant. At the same time the BVDV-1 neutralizing antibody titer reached a maximum of 12 Nu on the 42nd day post vaccination.

CONCLUSIONS

The H2O2 inactivated BVDV vaccine has good safety and immunogenicity, which provides a potential solution for the further development of an efficient and safe BVDV vaccine.

Immune evasion strategies of bovine viral diarrhea virus

Bovine viral diarrhea virus (BVDV) is a significant pathogen that causes great economic losses in the global livestock industry. During the long-term interactions between BVDV and its hosts, the virus has evolved multiple strategies to evade the host's innate immunity and adaptive immunity, thereby promoting viral survival and replication. This review focuses on the most recent research on immune evasion strategies employed by BVDV, including evading type I IFN signaling pathway, evading host adaptive immunity, mediating NF-κB signaling pathway, mediating cell apoptosis and inducing autophagy. Unraveling BVDV's immune evasion strategies will enhance our understanding of the pathogenesis of BVDV and contribute to the development of more effective therapies for the prevention, control and eradication of BVDV.

The combination of vaccines and adjuvants to prevent the occurrence of high incidence of infectious diseases in bovine

As the global population grows, the demand for beef and dairy products is also increasing. The cattle industry is facing tremendous pressures and challenges. The expanding cattle industry has led to an increased risk of disease in cattle. These diseases not only cause economic losses but also pose threats to public health and safety. Hence, ensuring the health of cattle is crucial. Vaccination is one of the most economical and effective methods of preventing bovine infectious diseases. However, there are fewer comprehensive reviews of bovine vaccines available. In addition, the variable nature of bovine infectious diseases will result in weakened or even ineffective immune protection from existing vaccines. This shows that it is crucial to improve overall awareness of bovine vaccines. Adjuvants, which are crucial constituents of vaccines, have a significant role in enhancing vaccine response. This review aims to present the latest advances in bovine vaccines mainly including types of bovine vaccines, current status of development of commonly used vaccines, and vaccine adjuvants. In addition, this review highlights the main challenges and outstanding problems of bovine vaccines and adjuvants in the field of research and applications. This review provides a theoretical and practical basis for the eradication of global bovine infectious diseases.

Integrative Transcriptomics and Proteomics Analysis Reveals Immune Response Process in Bovine Viral Diarrhea Virus-1-Infected Peripheral Blood Mononuclear Cells

Bovine viral diarrhea virus (BVDV) causes bovine viral diarrhea-mucosal disease, inflicting substantial economic losses upon the global cattle industry. Peripheral blood mononuclear cells (PBMCs) are the central hub for immune responses during host-virus infection and have been recognized as crucial targets for BVDV infection. In order to elucidate the dynamics of host-BVDV-1 interaction, this study harnessed RNA-seq and iTRAQ methods to acquire an extensive dataset of transcriptomics and proteomics data from samples of BVDV-1-infected PBMCs at the 12-h post-infection mark. When compared to mock-infected PBMCs, we identified 344 differentially expressed genes (DEGs: a total of 234 genes with downregulated expression and 110 genes with upregulated expression) and 446 differentially expressed proteins (DEPs: a total of 224 proteins with downregulated expression and 222 proteins with upregulated expression). Selected DEGs and DEPs were validated through quantitative reverse transcriptase-polymerase chain reaction and parallel reaction monitoring. Gene ontology annotation and KEGG enrichment analysis underscored the significant enrichment of DEGs and DEPs in various immunity-related signaling pathways, including antigen processing and presentation, complement and coagulation cascades, cytokine-cytokine receptor interaction, and the NOD-like receptor signaling pathway, among others. Further analysis unveiled that those DEGs and DEPs with downregulated expression were predominantly associated with pathways such as complement and coagulation cascades, the interleukin-17 signaling pathway, cytokine-cytokine receptor interaction, the PI3K-Akt signaling pathway, the tumor necrosis factor signaling pathway, and the NOD-like receptor signaling pathway. Conversely, upregulated DEGs and DEPs were chiefly linked to metabolic pathways, oxidative phosphorylation, complement and coagulation cascades, and the RIG-I-like receptor signaling pathway. These altered genes and proteins shed light on the intense host-virus conflict within the immune realm. Our transcriptomics and proteomics data constitute a significant foundation for delving further into the interaction mechanism between BVDV and its host.

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.

Transcriptome analysis reveals the regulatory effects of Bacillus amyloliquefaciens and Bacillus pumilus on immune and digestive related genes in the spleen of weanling black goats

The aim of the present study was to evaluate the effects of Bacillus amyloliquefaciens fsznc-06 and Bacillus pumilus fsznc-09 on the expressions of spleen genes in weanling Jintang black goats. Bacillus amyloliquefaciens fsznc-06 (BA-treated group) and Bacillus pumilus fsznc-09 (BP-treated group) were directly fed to goats, and the spleens were harvested for transcriptome analysis. The KEGG pathway analysis showed that the differentially expressed genes (DEGs) in BA-treated vs CON group were mainly involved in digestive system and immune system, while those in BP-treated vs CON group were mainly involved in immune system, and those in BA-treated vs BP-treated group were mainly involved in digestive system. In conclusion, Bacillus amyloliquefaciens fsznc-06 might promote the expressions of genes related to immune system and digestive system, reduce the expressions of disease genes related to digestive system and might promote mutual accommodation of some immune genes in weanling black goat. Bacillus pumilus fsznc-09 might promote the expressions of genes related to immune system and mutual accommodation of some immune genes in weanling black goat. Bacillus amyloliquefaciens fsznc-06 has advantages over Bacillus pumilus fsznc-09 in promoting the expressions of genes related to digestive system and mutual accommodation of some immune genes.

Upregulation of host genes during disease progression in bovine leukemia virus infection is independent of overexpression of viral transcriptional regulators in vitro

Bovine leukemia virus (BLV) is a member of the genus Deltaretrovirus within the family Retroviridae that infects bovine B cells, causing persistent lymphocytosis and enzootic bovine leukosis (EBL) in a small fraction of infected cattle. As changes in the transcriptome of infected cells are important for BLV disease progression, comprehensive analysis of gene expression in different disease states is required. In this study, we performed an RNA-seq analysis using samples from non-EBL cattle with and without BLV infection. Subsequently, a transcriptome analysis was conducted in combination with previously obtained RNA-seq data from EBL cattle. We found several differentially expressed genes (DEGs) between the three groups. After screening and confirmation of target DEGs using real-time reverse transcription polymerase chain reaction, we found that 12 target genes were significantly upregulated in EBL cattle compared to BLV-infected cattle without lymphoma. In addition, the expression levels of B4GALT6, ZBTB32, EPB4L1, RUNX1T1, HLTF, MKI67, and TOP2A were significantly and positively correlated with the proviral load in BLV-infected cattle. Overexpression experiments revealed that these changes were independent of BLV tax or BLV AS1-S expression in vitro. Our study provides additional information on host gene expression during BLV infection and EBL development, which may be helpful for understanding the complexity of transcriptome profiles during disease progression.

Bacteroides vulgatus attenuates experimental mice colitis through modulating gut microbiota and immune responses

Introduction

Bacteroides vulgatus is one of the predominant Bacteroides species in the human gut and exerts a series of beneficial effects. The aim of this study was to investigate the protective role of B. vulgatus Bv46 in a dextran sodium sulfate (DSS) induced colitis mouse model.

Methods

Female C57BL/6J mice were given 3% DSS in drinking water to induce colitis and simultaneously treated with B. vulgatus Bv46 by gavage for 7 days. Daily weight and disease activity index (DAI) of mice were recorded, and the colon length and histological changes were evaluated. The effects of B. vulgatus Bv46 on gut microbiota composition, fecal short chain fatty acids (SCFAs) concentration, transcriptome of colon, colonic cytokine level and cytokine secretion of RAW 264_·7 macrophage cell line activated by the lipopolysaccharide (LPS) were assessed.

Results and Discussion

B. vulgatus Bv46 significantly attenuated symptoms of DSS-induced colitis in mice, including reduced DAI, prevented colon shortening, and alleviated colon histopathological damage. B. vulgatus Bv46 modified the gut microbiota community of colitis mice and observably increased the abundance of Parabacteroides, Bacteroides, Anaerotignum and Alistipes at the genus level. In addition, B. vulgatus Bv46 treatment decreased the expression of colonic TNF-α, IL-1β and IL-6 in DSS-induced mouse colitis in vivo, reduced the secretion of TNF-α, IL-1β and IL-6 in macrophages stimulated by LPS in vitro, and downregulated the expression of Ccl19, Cd19, Cd22, Cd40 and Cxcr5 genes in mice colon, which mainly participate in the regulation of B cell responses. Furthermore, oral administration of B. vulgatus Bv46 notably increased the contents of fecal SCFAs, especially butyric acid and propionic acid, which may contribute to the anti-inflammatory effect of B. vulgatus Bv46. Supplementation with B. vulgatus Bv46 serves as a promising strategy for the prevention of colitis.

Polypharmacology-based approach for screening TCM against coinfection of Mycoplasma gallisepticum and Escherichia coli

Natural products and their unique polypharmacology offer significant advantages for finding novel therapeutics particularly for the treatment of complex diseases. Meanwhile, Traditional Chinese Medicine exerts overall clinical benefits through a multi-component and multi-target approach. In this study, we used the previously established co-infection model of Mycoplasma gallisepticum and Escherichia coli as a representative of complex diseases. A new combination consisting of 6 herbs were obtained by using network pharmacology combined with transcriptomic analysis to reverse screen TCMs from the Chinese medicine database, containing Isatdis Radix, Forsythia Fructus, Ginkgo Folium, Mori Cortex, Licorice, and Radix Salviae. The results of therapeutic trials showed that the Chinese herbal compounds screened by the target network played a good therapeutic effect in the case of co-infection. In summary, these data suggested a new method to validate target combinations of natural products that can be used to optimize their multiple structure-activity relationships to obtain drug-like natural product derivatives.

Comparative transcriptome analysis of MDBK cells reveals that BoIFN-γ augmented host immune responses to bovine herpesvirus 1 infection

Bovine herpesvirus 1 (BoHV-1) is an alphaherpesvirus that causes infectious bovine rhinotracheitis and infectious pustular vulvovaginitis in cattle. Ιnterferon-gamma (IFN-γ) is a pleiotropic cytokine with antiviral activity that modulates the innate and adaptive immune responses. In this study, we prepared high-purity bovine interferon gamma (BoIFN-γ) dimer protein using prokaryotic expression system and affinity chromatography. We subsequently investigated the effect of BoIFN-γ on BoHV-1 infection in Madin-Darby bovine kidney (MDBK) cells. The results showed that BoIFN-γ pre-treament not only decreased the production of BoHV-1 but also reduced the cytopathic effect of the virus. Differential gene expression profiles of BoHV-1 infected MDBK cells were then analysed through high-throughput RNA sequencing. The data showed that BoIFN-γ pre-treatment reduced lipid metabolism disorder and DNA damage caused by BoHV-1 infection. Furthermore, BoIFN-γ treatment upregulated the transcription of interferon regulatory transcription factors (IRF1 and GBP5) and interferon-stimulated genes (ISGs) of MDBK cells. Additionally, BoIFN-γ promotes expression of cellular protein involved in complement activation and coagulation cascades response as well as antigen processing and presentation process, while BoHV-1 infection dramatically downregulates transcription of these immune components including C3, C1r, C1s, PLAT, ITGB2, PROCR, BoLA, CD74, B2M, PA28, BoLA-DRA, and TAPBP. Collectively, our findings revealed that BoIFN-γ pre-treatment can improve host resistance to BoHV-1 infection and regulate transcription or expression of host protein associated with cellular metabolism and innate immune response. This provides insights into the development of prophylactic agents for prevention and control of BoHV-1 infection.

Transcriptomic Analysis of MDBK Cells Infected with Cytopathic and Non-Cytopathic Strains of Bovine Viral Diarrhea Virus (BVDV)

Bovine viral diarrhea virus (BVDV) belongs to the Flaviviridae family and the Pestivirus genus. Infection with BVDV causes a disease with a wide spectrum of clinical symptoms, most often mild, although infections with this virus constitute a serious economic problem all over the world. The virus is characterized by a high genetic variability, while the accumulation of single mutations leads to the formation of its new variants. The aim of this study was to better understand the complicated pathogenesis of this disease at the molecular level via the analysis of the transcriptome of cells infected with this virus. The bovine kidney cell line (MDBK), the cytopathic (cp) reference strain, and two non-cytopathic (ncp) BVD virus field strains were used in transcriptomic studies. The cell transcriptome was tested 24 and 72 h after infection. The results of the microarray analysis revealed changes in the expression levels of numerous genes. Genes with changed expression as a result of infection with the cp strain caused changes in the expression levels of a large number of genes and enriched a number of pathways. Genes with increased expression levels were enriched among other pathways involved in the cell cycle, while genes with reduced expression levels enriched pathways mostly related to metabolism. Genes with increased expression levels as a result of infection with ncp strains enriched a much smaller number of pathways, among them, pathways related to signaling activity 24 h post-infection and serine biosynthetic pathways both 24 and 72 h post-infection. Pathways enriched by genes with reduced expression levels were related to the innate immune response (72 h post-infection) or metabolism (24 and 72 h post-infection). The results of microarray studies can help us to better understand the host’s response to BVDV infection.

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.

Integrative Transcriptomics and Proteomics Analysis Provide a Deep Insight Into Bovine Viral Diarrhea Virus-Host Interactions During BVDV Infection

Bovine viral diarrhea virus (BVDV) is the causative agent of bovine viral diarrhea-mucosal disease (BVD-MD), an important viral disease in cattle that is responsible for extensive economic losses to the cattle industry worldwide. Currently, several underlying mechanisms involved in viral replication, pathogenesis, and evading host innate immunity of BVDV remain to be elucidated, particularly during the early stage of virus infection. To further explore the mechanisms of BVDV-host interactions, the transcriptomics and proteomics profiles of BVDV-infected MDBK cells were sequenced using RNA-seq and iTRAQ techniques, respectively, and followed by an integrative analysis. Compared with mock-infected MDBK cells, a total of 665 differentially expressed genes (DEGs) (391 down-regulated, 274 up-regulated) and 725 differentially expressed proteins (DEPs) (461 down-regulated, 264 up-regulated) were identified. Among these, several DEGs and DEPs were further verified using quantitative RT-PCR and western blot. Following gene ontology (GO) annotation and KEGG enrichment analysis, we determined that these DEGs and DEPs were significantly enriched in multiple important cellular signaling pathways including NOD-like receptor, Toll-like receptor, TNF, NF-κB, MAPK, cAMP, lysosome, protein processing in endoplasmic reticulum, lipid metabolism, and apoptosis signaling pathways. Significantly, the down-regulated DEGs and DEPs were predominantly associated with apoptosis-regulated elements, inflammatory factors, and antiviral elements that were involved in innate immunity, thus, indicating that BVDV could inhibit apoptosis and the expression of host antiviral genes to facilitate viral replication. Meanwhile, up-regulated DEGs and DEPs were primarily involved in metabolism and autophagy signaling pathways, indicating that BVDV could utilize the host metabolic resources and cell autophagy to promote replication. However, the potential mechanisms BVDV-host interactions required further experimental validation. Our data provide an overview of changes in transcriptomics and proteomics profiles of BVDV-infected MDBK cells, thus, providing an important basis for further exploring the mechanisms of BVDV-host interactions.

Comparative transcriptome analysis reveals potential anti-viral immune pathways of turbot (Scophthalmus maximus) subverted by megalocytivirus RBIV-C1 for immune evasion

Successful viral infection and multiplication chiefly rely on virus subversion mechanisms against host anti-viral immune responses. In this study, in order to reveal the anti-viral immune-related pathways suppressed by megalocytivirus infection, transcriptome analysis was performed on the head-kidney of turbot (Scophthalmus maximus) infected with lethal dose of RBIV-C1 at 3, 6 and 9 days post challenge (dpc). The results showed that, compared to unchallenged groups, 190, 1220, and 3963 DEGs were detected in RBIV-C1 infected groups at 3, 6 and 9 dpc, respectively, of which, DEGs of complement components and pattern recognition proteins were up-regulated at 3 dpc and down-regulated at 6 and 9 dpc, DEGs of cytokines were up-regulated at 6 dpc and down-regulated at 9 dpc. Expression trend analysis revealed that DEGs of profiles 9 and 13 featured decreased expression patterns and were significantly enriched into 10 immune-related pathways, i.e., complement and coagulation cascades, cytokine-cytokine receptor interaction, chemokine signaling pathway, B/T cell receptor signaling pathway, antigen processing and presentation, and so on. Further co-expression network analysis (WGCNA) revealed positive correlated innate immune related pathways at 3 and 6 dpc, and negative correlated innate and adaptive immune related pathways at 9 dpc. This study revealed a set of anti-viral immune genes/pathways that would also be potential targets subverted by RBIV-C1 for immune evasion, which can serve as a valuable resource for future studies on the molecular mechanisms of anti-viral immune defense of turbot and immune escape of megalocytivirus.

Development of a One-Step Multiplex Real-Time PCR Assay for the Detection of Viral Pathogens Associated With the Bovine Respiratory Disease Complex

Bovine respiratory disease complex (BRDC) occurs widely in cattle farms. The main viral pathogens include bovine viral diarrhea virus (BVDV), Bovine herpesvirus 1 (BoHV-1), bovine parainfluenza virus type 3 (BPIV3), and bovine respiratory syncytial virus (BRSV), and the newly emerged influenza D virus (IDV). In this study, we have developed a one-step multiplex real-time Polymerase Chain Reaction (PCR) capable of simultaneously detecting these five viral pathogens causing BRDC. The established assay could specifically detect targeted viruses without cross-reaction with others. The detection limit was ~10 copies/reaction for single real-time PCR and 100 copies/ reaction for multiplex real-time PCR assay. A total of 213 nasal samples from cattle with signs of respiratory tract disease were then collected for performance evaluation of the established platform, proving that the method has good specificity and sensitivity. The surveillance data suggested that BVDV and BoHV-1 infections are the dominant cause of BRDC in the herd, whereas the detection rate of IDV, BIPV3, and BRSV is relatively lower. In summary, the established assay provides technical support for rapid clinical detection of BRDC associated viral pathogens to guide the formulation of BRDC prevention and control measures.

Genome-wide transcriptome analysis of porcine epidemic diarrhea virus virulent or avirulent strain-infected porcine small intestinal epithelial cells

Porcine epidemic diarrhea virus (PEDV) is the main cause of diarrhea, vomiting, and mortality in pigs, which results in devastating economic loss to the pig industry around the globe. In recent years, the advent of RNA-sequencing technologies has led to delineate host responses at late stages of PEDV infection; however, the comparative analysis of host responses to early-stage infection of virulent and avirulent PEDV strains is currently unknown. Here, using the BGI DNBSEQ RNA-sequencing, we performed global gene expression profiles of pig intestinal epithelial cells infected with virulent (GDS01) or avirulent (HX) PEDV strains for 3, 6, and 12 ​h. It was observed that over half of all significantly dysregulated genes in both infection groups exhibited a down-regulated expression pattern. Functional enrichment analyses indicated that the differentially expressed genes (DEGs) in the GDS01 group were predominantly related to autophagy and apoptosis, whereas the genes showing the differential expression in the HX group were strongly enriched in immune responses/inflammation. Among the DEGs, the functional association of TLR3 and IFIT2 genes with the HX and GDS01 strains replication was experimentally validated by TLR3 inhibition and IFIT2 overexpression systems in cultured cells. TLR3 expression was found to inhibit HX strain, but not GDS01 strain, replication by enhancing the IFIT2 expression in infected cells. In conclusion, our study highlights similarities and differences in gene expression patterns and cellular processes/pathways altered at the early-stage infection of PEDV virulent and avirulent strains. These findings may provide a foundation for establishing novel therapies to control PEDV infection.

•

Differentially expressed genes (DEGs) in the GDS01 group were predominantly related to autophagy and apoptosis.

•

DEGs in the HX group were strongly enriched in immune responses/inflammation.

•

TLR3 inhibits HX strain replication, but not GDS01 strain, by enhancing the IFIT2 expression.

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.

Screening of host genes regulated by ID1 and ID3 proteins during foot-and-mouth disease virus infection

Foot-and-mouth disease virus (FMDV) is an important pathogen that harms cloven-hoofed animals and has caused serious losses to livestock production since its discovery. Furthermore, inhibitor of DNA binding (ID) proteins have been thoroughly studied in tumorigenesis, differentiation and metastasis, but its role in viral infection is rarely known. In this study, three gene knockout cell lines ID1 KO, ID3 KO, ID1/3 KO were obtained based on BHK-21 cells. We found that ID1 and ID3 genes single or double knockout promote the replication of FMDV. Moreover, compared with negative control cells during virus infection, there were 551 up-regulated genes and 1222 down-regulated genes in the ID1 KO cell line; 916 up-regulated genes and 1845 down-regulated genes in the ID3 KO cell line; 810 up-regulated genes and 1566 down-regulated genes in ID1/3 KO cell line. Further genes expression patterns verification results also showed a good correlation between the data of RT-qRCR and RNA-seq. These findings provide a basis for studying the relevant mechanisms between host genes and ID genes during FMDV infection.

Network-based analysis on genetic variants reveals the immunological mechanism underlying Alzheimer's disease

Alzheimer's Disease (AD) is a neurodegenerative disorder characterized by the impairment of cognitive function and loss of memory. Previous studies indicate an essential role of immune response in AD, but the detailed mechanisms remain unclear. In this study, we obtained 1664 credible risk variants (CRVs) based on the most significant SNP detected by International Genomics of Alzheimer's Project, from which 99 genes (CRVs-related genes) were identified. Function analysis revealed that these genes were mainly involved in immune response and amyloid-β and its precursor metabolisms, indicating a potential role of immune response in regulating neurobiological processes in the etiology of neurodegenerative disease. Pathway crosstalk analysis revealed the complicated connections between immune-related pathways. Further, we found that the CRVs-related genes showed temporal-specific expression in the thalamus in adolescence developmental period. Cell type-specific expression analysis found that CRVs-related genes might be specifically expressed in brain cells such as astrocytes and oligodendrocytes. Protein-protein interaction network analysis identified the highly interconnected 'hub' genes, all of which were susceptible loci of AD. These results indicated that the CRVs may exert a potential influence in AD by regulating immune response, thalamus development, astrocytes activities, and amyloid-β binding. Our results provided hints for further experimental verification of AD pathophysiology.

Comparison of the Ability of High and Low Virulence Strains of Non-cytopathic Bovine Viral Diarrhea Virus-1 to Modulate Expression of Interferon Tau Stimulated Genes in Bovine Endometrium

Bovine Viral Diarrhea virus (BVDV) is a pestivirus with a single-stranded, positive sense RNA genome. It is endemic in many cattle populations, causing major economic losses in part due to reduced fertility. BVDV exhibits great genetic diversity and is classified as type 1 or 2 (BVDV-1, BVDV-2) with either non-cytopathogenic (ncp) or cytopathogenic (cp) biotypes. Differing strains of ncpBVDV differ in virulence, affecting clinical outcome. BVDV replicates in the reproductive tract, affecting host immunity and embryo survival. This study used an in vitro model of primary bovine endometrial cell cultures to compare the effects of two BVDV ncp type 1a strains of differing virulence (termed HO and KY) on endometrial transcription of candidate interferon stimulated genes (ISG) using qPCR. Half the cultures were stimulated with interferon tau (IFNT, the conceptus produced pregnancy recognition factor) in the presence or absence of viral infection. Cultures were replicated on cells from 10 BVDV-free cows. IFNT treatment stimulated transcription of 10 candidate ISGs, whereas both ncpBVDV-1 strains alone inhibited transcription of 8/10 ISGs. In combined BVDV-1+IFNT cultures, the stimulatory effect of IFNT on expression of GBP4, ISG15, HERC5, RSAD2, IFIH1, IFIT3, and MX1 was significantly inhibited by HO, but only ISG15, RSAD2, IFI27, and IFIT3 were decreased by KY. Inhibition by HO was generally greater. The IFNT-induced expression of TRIM56 was, however, increased by HO. These data show that HO, the more virulent ncpBVDV-1 strain, has a greater capacity to inhibit key antiviral pathways. These differences need confirmation at the protein level but may influence immune tolerance of the host. They could also reduce fertility by increasing uterine susceptibility to bacterial infection and disrupting IFNT-mediated pregnancy recognition.

Comprehensive Analyses of circRNA Expression Profiles and Function Prediction in Chicken Cecums After Eimeria tenella Infection

Coccidiosis is an important intestinal parasitic disease that causes great economic losses to the global poultry production industry. Circular RNAs (circRNAs) are long non-coding RNAs that play important roles in various infectious diseases and inflammatory responses. However, the expression profiles and functions of circRNAs during Eimeria tenella (E. tenella) infection remain unclear. In this study, high-throughput sequencing was carried out to detect circRNAs in chicken cecal tissues from the control (JC), resistant (JR), and susceptible (JS) groups on day 4.5 postinfection (pi), respectively. A total of 104 circRNAs were differentially expressed, including 47 circRNAs between the JS and JC groups, 38 between the JR and JS groups, and 19 between the JR and JC groups. Functional analyses indicated that these differentially expressed circRNAs were involved in pathways related to E. tenella infection; the adaptive immune response was enriched in the JS vs JC group, the NF-kappa B signaling and natural killer cell-mediated cytotoxicity pathways were enriched in the JS vs JC and JR vs JC groups, while the B cell receptor signaling pathway was enriched in only the JR vs JC group. Moreover, the coexpression network of differentially expressed circRNAs and mRNAs suggested that circRNA2202 and circRNA0759 associated with DTX1 in the JS vs JC group, circRNA4338 associated with VPREB3 and CXCL13L3 in the JR vs JC group, and circRNA2612 associated with IL8L1 and F2RL2 in the JR vs JS group were involved in the immune response upon E. tenella infection. In conclusion, our results provide valuable information on the circRNAs involved in the progression of chicken E. tenella infection and advance our understanding of the circRNA regulatory mechanisms of host resistance and susceptibility to E. tenella infection in chickens.

Immunogene expression analysis in betanodavirus infected-Senegalese sole using an OpenArray® platform

The transcriptomic response of Senegalese sole (Solea senegalensis) triggered by two betanodaviruses with different virulence to that fish species has been assessed using an OpenArray® platform based on TaqMan™ quantitative PCR. The transcription of 112 genes per sample has been evaluated at two sampling times in two organs (head kidney and eye/brain-pooled samples). Those genes were involved in several roles or pathways, such as viral recognition, regulation of type I (IFN-1)-dependent immune responses, JAK-STAT cascade, interferon stimulated genes, protein ubiquitination, virus responsive genes, complement system, inflammatory response, other immune system effectors, regulation of T-cell proliferation, and proteolysis and apoptosis. The highly virulent isolate, wSs160.3, a wild type reassortant containing a RGNNV-type RNA1 and a SJNNV-type RNA2 segments, induced the expression of a higher number of genes in both tested organs than the moderately virulent strain, a recombinant harbouring mutations in the protruding domain of the capsid protein. The number of differentially expressed genes was higher 2 days after the infection with the wild type isolate than at 3 days post-inoculation. The wild type isolate also elicited an exacerbated interferon 1 response, which, instead of protecting sole against the infection, increases the disease severity by the induction of apoptosis and inflammation-derived immunopathology, although inflammation seems to be modulated by the complement system. Furthermore, results derived from this study suggest a potential important role for some genes with high expression after infection with the highly virulent virus, such as rtp3, sacs and isg15. On the other hand, the infection with the mutant does not induce immune response, probably due to an altered recognition by the host, which is supported by a different viral recognition pathway, involving myd88 and tbkbp1.

Dynamic Transcriptome Reveals the Mechanism of Liver Injury Caused by DHAV-3 Infection in Pekin Duck

Duck hepatitis A virus 3 (DHAV-3) is a wild endemic virus, which seriously endangers the duck industry in China. The present study aims to elucidate the mechanism of duck resistance to DHAV-3 infection. Both resistant and susceptible ducks were challenged with DHAV-3 in this experiment. The histopathological features and serum biochemical indices (ALT and AST) were analyzed to estimate liver injury status at 6, 12, 15, and 24 h post-infection (hpi). The dynamic transcriptomes of liver were analyzed to explain the molecular regulation mechanism in ducks against DHAV-3. The result showed that the liver injury in susceptible ducks was more serious than that in the resistant ducks throughout the four time points. A total of 2,127 differentially expressed genes (DEGs) were identified by comparing the transcriptome of the two populations. The expression levels of genes involved in innate immune response increased rapidly in susceptible ducks from 12 hpi. Similarly, the expression of genes involved in cytokine regulation also increased at the same time points, while the expression levels of these genes in resistant ducks remained similar between the various time points. KEGG enrichment analysis of the DEGs revealed that the genes involved in cytokine regulation and apoptosis were highly expressed in susceptible ducks than that in resistant ducks, suggesting that excessive cytokine storm and apoptosis may partially explain the mechanism of liver injury caused by DHAV-3 infection. Besides, we found that the FUT9 gene may contribute to resistance towards DHAV-3 in resistant ducklings. These findings will provide insight into duck resistance and susceptibility to DHAV-3 infection in the early phases, facilitate the development of a strategy for DHAV-3 prevention and treatment, and enhance genetic resistance via genetic selection in animal breeding.

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.

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

Simple Summary

Due to the undeniable detrimental impact of bovine viral diarrhea virus (BVDV) on cattle worldwide, various preventive approaches are carried out to control the spread of this disease. Among the established preventive approaches, vaccination remains the most widely used cost-effective method of control. Hence, a deeper study into the host immune response to vaccines will further refine the efficacy of these vaccines; the identification of differentially expressed genes (DEGs) related to immune response might bring a long-lasting solution. Thus far, studies showing the genes related to the immune response of cattle to vaccines are still limited. Therefore, this study identified DEGs in animals with high and low sample to positive (S/P) ratio based on the BVDV antibody level, using RNA sequencing (RNA-seq) transcriptome analysis, and functional enrichment analysis in gene ontology (GO) annotations and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway. Results revealed that several upregulated and downregulated genes were significantly annotated to antigen processing and presentation (MHC class I), immune response, and interferon-gamma production, indicating the immune response of the animals related to possible shaping of their adaptive immunity against the BVDV type I. Moreover, significant enrichment to various KEGG pathways related to the development of adaptive immunity was observed.

Abstract

Immune response of 107 vaccinated Holstein cattle was initially obtained prior to the ELISA test. Five cattle with high and low bovine viral diarrhea virus (BVDV) type I antibody were identified as the final experimental animals. Blood samples from these animals were then utilized to determine significant differentially expressed genes (DEGs) using the RNA-seq transcriptome analysis and enrichment analysis. Our analysis identified 261 DEGs in cattle identified as experimental animals. Functional enrichment analysis in gene ontology (GO) annotations and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways revealed the DEGs potentially induced by the inactivated BVDV type I vaccine, and might be responsible for the host immune responses. Our findings suggested that inactivated vaccine induced upregulation of genes involved in different GO annotations, including antigen processing and presentation of peptide antigen (via MHC class I), immune response, and positive regulation of interferon-gamma production. The observed downregulation of other genes involved in immune response might be due to inhibition of toll-like receptors (TLRs) by the upregulation of the Bcl-3 gene. Meanwhile, the result of KEGG pathways revealed that the majority of DEGs were upregulated and enriched to different pathways, including cytokine-cytokine receptor interaction, platelet activation, extracellular matrix (ECM) receptor interaction, hematopoietic cell lineage, and ATP-binding cassette (ABC) transporters. These significant pathways supported our initial findings and are known to play a vital role in shaping adaptive immunity against BVDV type 1. In addition, type 1 diabetes mellitus pathways tended to be significantly enriched. Thus, further studies are needed to investigate the prevalence of type 1 diabetes mellitus in cattle vaccinated with inactivated and live BVDV vaccine.