Transcriptome analysis reveals differential immune related genes expression in bovine viral diarrhea virus-2 infected goat peripheral blood mononuclear cells (PBMCs)

Screening for immune biomarkers associated with infection or protection against Ehrlichia ruminantium by RNA-sequencing analysis

Heartwater is one of the most economically important tick-borne fatal diseases of livestock. The disease is caused by the bacteria Ehrlichia ruminantium transmitted by Amblyomma ticks. Although there is evidence that interferon-gamma controls E. ruminantium growth and that cellular immune responses are protective, an effective recombinant vaccine for this disease is lacking. Analyses of markers associated with infection as well as protection will lead to a better understanding of the E. ruminantium immune response and corresponding pathways induced in sheep peripheral blood mononuclear cells (PBMC) will assist in development of such a vaccine. In this study, Biomarkers of infection (BMI) were identified as uniquely expressed genes during primary infection and biomarkers of protection (BMP) associated with immune to heartwater were identified post challenge. Sheep were experimentally infected and challenged with E. ruminantium infected ticks. The immune phenotypic and transcriptome profile of their PBMC were compared to their own naïve PBMC collected before infection. The study revealed 305 differentially expressed genes (DEGs) as BMI, of these 17 were upregulated at all three time-points investigated. These DEGs, form part of the bacterial invasion of epithelial cells Kyoto Encyclopaedia of Genes and Genomes (KEGG) pathway, and others detected from day 1 post infection and are considered predictive markers for early heartwater infection in ruminants. Similarly, a total of 332 DEGs were identified as BMP, of these 100 were upregulated and 75 were downregulated at all three time-points investigated. However, at D1PC most DEGs were downregulated (n = 1312) that correlated with a reduction in the % CD4 and CD8 T cells detected with flow cytometry. KEGG pathway analyses showed complete down regulation of T cell specific pathways possibly due to homing of immune cells to the site of infection after acquired immunity developed. At D4PC, expression levels of most of these downregulated genes increased and by D6PC they were upregulated. This indicates that the sampling time-point for biomarker analyses is important when results for acquired immune responses are inferred. This data identified DEGs that could be considered as biomarkers of protective immunity that can be used for identification of vaccine antigens and provides a strong foundation to further development of heartwater recombinant vaccines.

Non-cytopathic bovine viral diarrhea virus (BVDV) inhibits innate immune responses via induction of mitophagy

Bovine viral diarrhea virus (BVDV) belongs to the genus Pestivirus within the family Flaviviridae. Mitophagy plays important roles in virus-host interactions. Here, we provide evidence that non-cytopathic (NCP) BVDV shifts the balance of mitochondrial dynamics toward fission and induces mitophagy to inhibit innate immune responses. Mechanistically, NCP BVDV triggers the translocation of dynamin-related protein (Drp1) to mitochondria and stimulates its phosphorylation at Ser616, leading to mitochondrial fission. In parallel, NCP BVDV-induced complete mitophagy via Parkin-dependent pathway contributes to eliminating damaged mitochondria to inhibit MAVS- and mtDNA-cGAS-mediated innate immunity responses, mtROS-mediated inflammatory responses and apoptosis initiation. Importantly, we demonstrate that the LIR motif of ERNS is essential for mitophagy induction. In conclusion, this study is the first to show that NCP BVDV-induced mitophagy plays a central role in promoting cell survival and inhibiting innate immune responses in vitro.

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.

Construction and Analysis of ceRNA Networks Reveal the Key Genes Associated with Bovine Herpesvirus Type 1 Infection

Background

Virus infection can cause the changes of lncRNA expression levels to regulate the interaction between virus and host, but the relationship between BHV-1 infection and lncRNA has not been reported.

Methods

In this study, in order to reveal the molecular mechanism of RNA in BoHV-1 infection, the Madin-Darby bovine kidney (MDBK) cells were infected with BoHV-1, transcriptome sequencing were performed by next-generation sequencing at 18 h or 24 h or 33 h of viral infection and then based on the competitive endogenous RNA (ceRNA) theory, lncRNA-miRNA-mRNA networks were constructed using these high-throughput sequencing data. The network analysis, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were performed for functional annotation and exploration of ncRNA ceRNAs in BoHV-1 infection.

Results

The results showed that 48 lncRNAs, 123 mRNAs and 20 miRNAs as differentially expressed genes, and the mitogen activated protein kinase (MAPK) pathway and calcium signaling pathway were significantly enriched in the ceRNA network. Some differentially expressed lncRNA genes were randomly selected for verification by RT-qPCR, and the results showed that their expression trend was consistent with the results of transcriptome sequencing data.

Conclusion

This study revealed that BoHV-1 infection can affect the expression of RNAs in MDBK cells and the regulation of ceRNA network to carry out corresponding biological functions in the host, but further experimental studies are still necessary to prove the hub genes function in ceRNA network and the molecular mechanism in BoHV-1 infection.

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.

Transcriptomic and Metabolomic Changes Reveal the Immunomodulatory Function of Casein Phosphopeptide-Selenium Chelate in Beagle Dogs

Casein phosphopeptide-selenium chelate (CPP-Se) is an organic compound produced by the chelation of casein phosphopeptide with selenium. This compound showed the ability to modulate canine immune response in our previous study; but its effect on the peripheral blood transcriptome and serum metabolome was unknown. This study aims to reveal the potential mechanism behind the immunomodulatory function of CPP-Se. We have identified 341 differentially expressed genes (DEGs) in CPP-Se groups as compared to the control group which comprised 110 up-regulated and 231 down-regulated genes. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis found that DEGs were mainly involved in immune-related signaling pathways. Moreover, the immune-related DEGs and hub genes were identified. Similarly, metabolomics identified 53 differentially expressed metabolites (DEMs) in the CPP-Se group, of which 17 were up-regulated and 36 were down-regulated. The pathways mainly enriched by DEMs were primary bile acid biosynthesis, tryptophan metabolism, and other amino acids metabolic pathways. Combined analysis of transcriptomic and metabolomic data showed that the DEGs and DEMs were commonly enriched in fatty acid biosynthesis, pyrimidine metabolism, glutathione metabolism, and glycerolipid metabolic pathways. Taken together, our findings provided a theoretical basis for further understanding of the immunomodulatory function of CPP-Se as well as a scientific reference for the future use of CPP-Se in pet foods as a dietary supplement to modulate the immunity.

Accumulation of Fat Not Responsible for Femoral Head Necrosis, Revealed by Single-Cell RNA Sequencing: A Preliminary Study

The etiology of osteonecrosis of the femoral head (ONFH) is not yet fully understood. However, ONFH is a common disease with high morbidity, and approximately one-third of cases are caused by glucocorticoids. We performed single-cell RNA sequencing of bone marrow to explore the effect of glucocorticoid on ONFH. Bone marrow samples of the proximal femur were extracted from four participants during total hip arthroplasty, including two participants diagnosed with ONFH for systemic lupus erythematosus (SLE) treated with glucocorticoids (the case group) and two participants with femoral neck fracture (the control group). Unbiased transcriptome-wide single-cell RNA sequencing analysis and computational analyses were performed. Seventeen molecularly defined cell types were identified in the studied samples, including significantly dysregulated neutrophils and B cells in the case group. Additionally, fatty acid synthesis and aerobic oxidation were repressed, while fatty acid beta-oxidation was enhanced. Our results also preliminarily clarified the roles of the inflammatory response, substance metabolism, vascular injury, angiogenesis, cell proliferation, apoptosis, and dysregulated coagulation and fibrinolysis in glucocorticoid-induced ONFH. Notably, we list the pathways that were markedly altered in glucocorticoid-induced ONFH with SLE compared with femoral head fracture, as well as their common genes, which are potential early therapeutic targets. Our results provide new insights into the mechanism of glucocorticoid-induced ONFH and present potential clues for effective and functional manipulation of human glucocorticoid-induced ONFH, which could improve patient outcomes.

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.

Goat CCL5 promotes cell viability and inflammatory factors production in lung fibroblasts and macrophages

Inflammatory chemokine CCL5 can mediate the occurrence of inflammatory reactions and participate in various disease processes. (Ch)CCL5 gene of Jintang black goat (Capra hircus, C. hircus) was cloned. The CDS (coding sequences) was 276 bp in length and encoded 91 amino acids. The 26.5 kDa recombinant protein was expressed by Escherichia coli system and purified by Ni-Agarose. The viabilities of primary goat lung fibroblasts could be enhanced after treating with ChCCL5 protein (12.5, 25, 50 μg/mL) (P < 0.05). The expression levels of interleukin-1beta (IL-1β), interleukin 6 (IL-6), tumor necrosis factor (TNF-α), C-C motif chemokine ligand 2 (CCL2) and heat-shock proteins (Hsp70) genes were upregulated after treating with ChCCL5 protein (12.5, 25, 50 μg/mL). Besides, the viabilities and phagocytic abilities of primary mouse peritoneal macrophages could be enhanced after treating with ChCCL5 protein (12.5, 25, 50 μg/mL) (P < 0.05). The expression levels of IL-1β, IL-6, toll-like receptor 4 (TLR4), inducible nitric oxide synthase (iNOs) and TNF-α genes were upregulated after treating with ChCCL5 protein (12.5, 25, 50 μg/mL) (P < 0.05). These results indicated that goat CCL5 might play a role in the inflammatory response by regulating the inflammatory cytokines produced by lung fibroblasts and macrophages.

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.

A Comparative Transcriptomic Analysis Reveals That HSP90AB1 Is Involved in the Immune and Inflammatory Responses to Porcine Deltacoronavirus Infection

PDCoV is an emerging enteropathogenic coronavirus that mainly causes acute diarrhea in piglets, seriously affecting pig breeding industries worldwide. To date, the molecular mechanisms of PDCoV-induced immune and inflammatory responses or host responses in LLC-PK cells in vitro are not well understood. HSP90 plays important roles in various viral infections. In this study, HSP90AB1 knockout cells (HSP90AB1^KO) were constructed and a comparative transcriptomic analysis between PDCoV-infected HSP90AB1^WT and HSP90AB1^KO cells was conducted using RNA sequencing to explore the effect of HSP90AB1 on PDCoV infection. A total of 1295 and 3746 differentially expressed genes (DEGs) were identified in PDCoV-infected HSP90AB1^WT and HSP90AB1^KO cells, respectively. Moreover, most of the significantly enriched pathways were related to immune and inflammatory response-associated pathways upon PDCoV infection. The DEGs enriched in NF-κB pathways were specifically detected in HSP90AB1^WT cells, and NF-κB inhibitors JSH-23, SC75741 and QNZ treatment reduced PDCoV infection. Further research revealed most cytokines associated with immune and inflammatory responses were upregulated during PDCoV infection. Knockout of HSP90AB1 altered the upregulated levels of some cytokines. Taken together, our findings provide new insights into the host response to PDCoV infection from the transcriptome perspective, which will contribute to illustrating the molecular basis of the interaction between PDCoV and HSP90AB1.

Age-dependent immune response in pigs against foot-and-mouth disease virus in vitro

Foot-and-mouth disease, one of the most contagious diseases in cloven-hoofed animals, causes significant economic losses. The pathogenesis of foot-and-mouth disease virus (FMDV) infection is known to differ with age of the animals. In this study, we aimed to reveal the difference in immunological response in the initial stage of FMDV infection between piglets and adult pigs. Peripheral blood mononuclear cells (PBMCs) were isolated from 3 piglets (8 weeks old) and 3 pigs (35 weeks old) that were not vaccinated against FMDV. O-type FMDV (2 × 10² median tissue culture infectious dose) was inoculated into porcine PBMCs and the cells were incubated at 37.0°C under 5% CO₂ for various time periods (0, 1, 3, 6, 12, 24, and 48 h). The total RNA was obtained from the FMDV-inoculated PBMCs after each time point, and the virus titer was investigated in these RNA samples. Furthermore, dynamics of mRNA expression of the six tested cytokines (interferon [IFN]-α, IFN-γ, ‎interleukin [IL]-6, IL-8, IL-10, and tumor necrosis factor [TNF]-α) in FMDV-inoculated porcine PBMCs were evaluated by time-series analysis to determine the differences, if any, based on the age of the pigs. The PBMCs of piglets contained the highest quantity of FMDV mRNA at 6 hours post-inoculation (hpi), and the PBMCs of pigs had the highest quantity of FMDV mRNA at 3 hpi. The mean cycle threshold-value in the PBMCs steadily decreased after the peak time point in the piglets and pigs (6 and 3 hpi, respectively). The dynamics of mRNA expression of all cytokines except TNF-α showed age-dependent differences in FMDV-inoculated PBMCs. The mRNA expression of most cytokines was more pronounced in the piglets than in the pigs, implying that the immune response against FMDV showed an age-dependent difference in pigs. In conclusion, within 48 hpi, the 8-week-old piglets responded more rapidly and were more sensitive to FMDV infection than the 35-week-old pigs, which could be associated with the difference in the pathogenesis of FMDV infection among the pigs. These results provide valuable insights into the mechanisms underlying the age-dependent differences in immune response in pigs against FMDV infection.

Identification of differentially expressed gene pathways between cytopathogenic and non-cytopathogenic BVDV-1 strains by analysis of the transcriptome of infected primary bovine cells

The bovine viral diarrhea virus 1 (BVDV-1), belonging to the Pestivirus genus, is characterized by the presence of two biotypes, cytopathogenic (cp) or non-cytopathogenic (ncp). For a better understanding of the host pathogen interactions, we set out to identify transcriptomic signatures of bovine lung primary cells (BPCs) infected with a cp or a ncp strain. For this, we used both a targeted approach by reverse transcription droplet digital PCR and whole genome approach using RNAseq. Data analysis showed 3571 differentially expressed transcripts over time (Fold Change >2) and revealed that the most deregulated pathways for cp strain are signaling pathways involved in responses to viral infection such as inflammatory response or apoptosis pathways. Interestingly, our data analysis revealed a deregulation of Wnt signaling pathway, a pathway described in embryogenesis, that was specifically seen with the BVDV-1 cp but not the ncp suggesting a role of this pathway in viral replication.

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.

Transcriptome and Proteomic Analysis Reveals Up-Regulation of Innate Immunity-Related Genes Expression in Caprine Herpesvirus 1 Infected Madin Darby Bovine Kidney Cells

Caprine herpesvirus 1 (CpHV-1) is a member of the alpha subfamily of herpesviruses, which is responsible for genital lesions and latent infections in goat populations worldwide. In this study, for the first time, the transcriptome and proteomics of CpHV-1 infected Madin Darby bovine kidney (MDBK) cells were explored using RNA-Sequencing (RNA-Seq) and isobaric tags for relative and absolute quantitation-liquid chromatography tandem mass spectrometry (iTRAQ-LC-MS/MS) technology, respectively. RNA-Seq analysis revealed 81 up-regulated and 19 down-regulated differentially expressed genes (DEGs) between infected and mock-infected MDBK cells. Bioinformatics analysis revealed that most of these DEGs were mainly involved in the innate immune response, especially the interferon stimulated genes (ISGs). Gene Ontology (GO) enrichment analysis results indicated that the identified DEGs were significantly mainly enriched for response to virus, defense response to virus, response to biotic stimulus and regulation of innate immune response. Viral carcinogenesis, the RIG-I-like receptor signaling pathway, the cytosolic DNA-sensing pathway and pathways associated with several viral infections were found to be significantly enriched in the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway database. Eleven selected DEGs (Mx1, RSAD2, IFIT1, IFIT2, IFIT5, IFIH1, IFITM3, IRF7, IRF9, OAS1X and OAS1Y) associated with immune responses were selected, and they exhibited a concordant direction both in RNA-Seq and quantitative real-time RT-PCR analysis. Proteomic analysis also showed significant up-regulation of innate immunity-related proteins. GO analysis showed that the differentially expressed proteins were mostly enriched in defense response and response to virus, and the pathways associated with viral infection were enriched under KEGG analysis. Protein-protein interaction network analysis indicated most of the DEGs related to innate immune responses, as DDX58(RIG-I), IFIH1(MDA5), IRF7, Mx1, RSAD2, OAS1 and IFIT1, were located in the core of the network and highly connected with other DGEs. Our findings support the notion that CpHV-1 infection induced the transcription and protein expression alterations of a series of genes related to host innate immune response, which helps to elucidate the resistance of host cells to viral infection and to clarify the pathogenesis of CpHV-1.

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

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

Common and Differential Dynamics of the Function of Peripheral Blood Mononuclear Cells between Holstein and Jersey Cows in Heat-Stress Environment

Simple Summary

Seasonal change, particularly changing to hot and humid season, has a negative effect on dairy cows in various ways, including productivity, reproduction, metabolism, and immunity. In high-temperature and humid weather, dairy cows are vulnerable to diseases by weakened immune system. However, the cause of this has not been fully described. Therefore, this study aims to understand changes of specific gene expression and immune pathways based on transcriptome analysis from peripheral blood mononuclear cells of Holstein and Jersey dairy cows between normal and heat-stress environmental conditions. We observed that the two breeds of dairy cow have common and different immune shifts according to the changes of temperature and humidity condition. Overall, the findings of this study improve the understanding of the underlying mechanisms by which seasonal changes affect dairy cow immunity.

Abstract

Heat stress has been reported to affect the immunity of dairy cows. However, the mechanisms through which this occurs are not fully understood. Two breeds of dairy cow, Holstein and Jersey, have distinct characteristics, including productivity, heat resistance, and disease in high-temperature environments. The objective of this study is to understand the dynamics of the immune response of two breeds of dairy cow to environmental change. Ribonucleic acid sequencing (RNA-seq) results were analyzed to characterize the gene expression change of peripheral blood mononuclear cells (PBMCs) in Holstein and Jersey cows between moderate temperature-humidity index (THI) and high THI environmental conditions. Many of the differentially expressed genes (DEGs) identified are associated with critical immunological functions, particularly phagocytosis, chemokines, and cytokine response. Among the DEGs, CXCL3 and IL1A were the top down-regulated genes in both breeds of dairy cow, and many DEGs were related to antimicrobial immunity. Functional analysis revealed that cytokine and chemokine response-associated pathways in both Holstein and Jersey PBMCs were the most important pathways affected by the THI environmental condition. However, there were also breed-specific genes and pathways that altered according to THI environmental condition. Collectively, there were both common and breed-specific altered genes and pathways in Holstein and Jersey cows. The findings of this study expand our understanding of the dynamics of immunity in different breeds of dairy cow between moderate THI and high THI environmental conditions.

The diverse roles of RIP kinases in host-pathogen interactions

Receptor Interacting Protein Kinases (RIPKs) are cellular signaling molecules that are critical for homeostatic signaling in both communicable and non-communicable disease processes. In particular, RIPK1, RIPK2, RIPK3 and RIPK7 have emerged as key mediators of intracellular signal transduction including inflammation, autophagy and programmed cell death, and are thus essential for the early control of many diverse pathogenic organisms. In this review, we discuss the role of each RIPK in host responses to bacterial and viral pathogens, with a focus on studies that have used pathogen infection models rather than artificial stimulation with purified pathogen associated molecular patterns. We also discuss the intricate mechanisms of host evasion by pathogens that specifically target RIPKs for inactivation, and finally, we will touch on the controversial issue of drug development for kinase inhibitors to treat chronic inflammatory and neurological disorders, and the implications this may have on the outcome of pathogen infections.

Identification and Validation of Stage-Associated PBMC Biomarkers in Breast Cancer Using MS-Based Proteomics

Background: It is well-described that the transcriptome of peripheral blood mononuclear cells (PBMCs) can be altered in the context of many malignancies to allow them avoid the effective immune response, which leads to cancer invasiveness. Here, we used an MS-based strategy to discover biomarkers in the PBMCs of breast cancer (BC) patients and validated them at different stages of BC. Methods: PBMCs were isolated from the breast cancer patients and were cultured alone or co-cultured with breast cancer cell lines. The role of PBMC in the invasion property of breast cancer cells was explored. NF-kB activity was also measured in the co-cultured breast cancer cells. Identification of protein profiles in the secretome and proteome of the co-cultured PBMCs was performed using SWATH mass spectrometry. Pathway enrichment and gene ontology analyses were carried out to look for the molecular pathways correlated with the protein expression profile of PBMCs in the breast cancer patients. Quantitative real-time polymerase chain reaction (qPCR) was performed to validate the candidate genes in the PBMC fraction of the breast cancer patients at the primary and metastatic stages. In silico survival analysis was performed to assess the potential clinical biomarkers in these PBMC subtypes. Results: PBMCs could significantly increase the invasion property of the BC cells concomitant with a decrease in E-cadherin and an increase in both Vimentin and N-cadherin expression. The NF-kB activity in the BC cells significantly increased following co-culturing implying the role of PBMCs in EMT induction. Enrichment analysis showed that the differentially expressed proteins in PBMCs are mainly associated with IL-17, PI3K-Akt, and HIF-1 signaling pathway, in which a set of seven proteins including TMSB4X, HSPA4, S100A9, SRSF6, THBS1, CUL4A, and CANX were frequently expressed. Finally, in silico analysis confirmed that a gene set consisting of S100A9, SRSF6, THBS1, CUL4A, and CANX were found to provide an insight for the identification of metastasis in breast cancer patients. Conclusion: In conclusion, our study revealed that the protein expression profile in PBMCs is a reflection of the proteins expressed in the BC tissue itself; however, the abundance level is different due to the stage of cancer.

Viperin protein inhibits the replication of caprine parainfluenza virus type 3 (CPIV 3) by interaction with viral N protein

Caprine parainfluenza virus type3 (CPIV3) is a newly identified member of Paramyxoviridae family. CPIV3 is highly prevalence in China and showed pathogenicity to goats; in addition, CPIV3 infection causes severe clinical disease under stress and/or co-infection conditions. Viperin is one of the hundreds of interferon-stimulated genes (ISGs), and possesses a wide range of antiviral activities. The aim of this study was to systemically explore the anti-CPIV3 activity of ruminants' Viperin. CPIV3 infection up-regulated Viperin transcription but not protein expression in MDBK cells. Bovine and caprine Viperin genes (bVi and gVi) were amplified and analyzed by BLAST and multiple alignment. The obtained bVi/gVi amino acid sequences showed 99.5%-100% identity with previously submitted sequences and has variants at N-terminal domain (1-70aa) between each other. The pcDNA3.1 plasmids containing bVi and gVi genes were constructed to over-express the target proteins. CPIV3 was inoculated in MDBK cells over-expressing bVi/gVi and viral load was detected by qRT-PCR, virus titration and Western blot. Both of the bVi and gVi significantly inhibited CPIV3 genome copy numbers and viral titers at 24 and 48 hpi (P < 0.01); and viral N protein expression was also decreased, comparing with those of mock transfected group. The last 50aa C-terminal region was crucial for its anti-CPIV3 activity. In addition, the over-expression of bVi/gVi did not influence CPIV3 binding, entry and release in the cells. These results indicated the anti-CPIV3 activity occurred in viral RNA/protein synthesis progress of the viral replication cycle. The Viperin also showed similar inhibitory effect on different CPIV3 strains. The potential interaction of Viperin with viral proteins (N, P, C and V) was determined by confocal laser scanning microscopy and Co-IP assay. Co-localization of Viperin with N, P or C, but not V, was observed; while only N protein direct interacted with Viperin in Co-IP test, no matter using viral protein expressing plasmids transfected or CPIV3 infected cell samples. In conclusion, the bVi and gVi Viperin effectively inhibited CPIV3 replication potentially via the interaction of Viperin with viral N protein. The present results gave more information about antiviral activity of ruminants Viperin and provided foundation for further studies of the interaction of Viperin with CPIV3 and other related viruses.

Development and evaluation of a real-time PCR panel for the detection of 20 immune markers in cattle and sheep

The establishment of a panel of immune markers is of paramount importance to understand the different transcription patterns of infectious diseases in livestock. The array of commercially available immunological assays for cattle and sheep is currently limited, due to the lack of antibodies for these species. Even though SYBR Green based real time quantitative PCR (qPCR) is the most commonly used method to study cytokine transcription in ruminants, a lack of standardization impairs its implementation in the study of different ruminant diseases. In order to obtain reliable qPCR results, several variables need to be considered: choice of reference genes for optimal normalization, variation of annealing temperature among primer sets, and assay specificity and sensitivity. In this study, we developed and validated a panel of immune markers in bovine and ovine samples using SYBR Green based qPCR in a cost-effective way with multiple primer sets optimised to amplify at a common thermal cycling temperature. Twenty primer sets were designed to quantify immune markers (IL-1b, IL-2, IL-4, IL-5, IL-6, IL-10, IL-12, IL-13, IL-15, IL-18, IL-23, TNF-α, IFN-γ, IFN-α, Ki-67, NFkB-65, TLR-3, TLR-4, TLR-8 and Rig-1) in ovine and bovine templates. For optimal normalization and selection of suitable reference genes, primer sets that measure the transcription of five reference genes were also included in the panel. The amplification efficiency, linearity and specificity was validated for all target genes. Optimal amplification conditions were achieved in both ovine and bovine samples for all gene targets, with the exception of Ki67. Relative quantification studies were performed on ovine and bovine mRNA obtained from sheep peripheral blood mononuclear cells (PBMCs) stimulated with three different treatments (PMA/Ionomycin, Concanavalin A (Con A) and pokeweed mitogen (PWM)). Pokeweed and ConA efficiently induced gene transcription of most of the targeted genes, while PMA/Ionomycin showed a weaker induction. Finally, we further assessed usability of our panel by running it on bovine monocyte derived dendritic cells (MoDCs) stimulated with different vaccines. Results confirmed the induction of a specific pro-inflammatory gene transcription pattern by rabies vaccine, which resembles the one occurring during viral infection. Altogether, we validated the efficiency and usability of an extended real-time PCR panel that gives the possibility to rapidly measure a broad spectrum of ovine and bovine immune markers by using a single set of reagents and protocol thus representing a valid and cost-effective tool for research purposes.

Identification of protein inhibitor of activated STAT 4, a novel host interacting partner that involved in bovine viral diarrhea virus growth

Background

Bovine viral diarrhea virus (BVDV) belongs to the Flaviviridae family and the pestivius virus group. BVDV is responsible for significant economic loss in cattle industry worldwide because of reducing reproductive performance, increasing incidence of other diseases and mortality among young stock. The core (C) protein of the Flaviviridae family member is involved in host antiviral immune response through activation of related signaling pathways that affect the viral replication. However, the influence of C protein-interaction partners in BVDV infections is poorly defined.

Methods

To explore C-protein-interacting partners, yeast two-hybrid was used to screen the interaction protein of C protein using bovine peripheral blood mononuclear cell (PBMC) cDNA library. The co-immunoprecipitation and confocal assays were manipulated to determine the interaction between potential partners and C protein. Knockdown and overexpression of the partner were used to examine whether the C-protein-interacting partner plays a role in BVDV proliferation and virulence. Meanwhile, qRT-PCR and western blot assays were used to investigate the effect of C protein and C-protein-interacting partner on the immune response of host cells.

Results

We identified protein inhibitor of activated STAT 4 (PIAS4) as a novel interacting partner of the BVDV C protein. Co-immunoprecipitation and confocal assays demonstrated a strong interaction between C protein and PIAS4. Silencing of PIAS4 with small interfering RNA suppressed C protein expression and BVDV growth, while overexpression of PISA4 increased C protein expression and BVDV growth. The overexpression of PIAS4 increased the cell apoptosis. Meanwhile, the expressions of STAT4, SOCS3, IFITM, IFN-α were negatively regulated by the expression of PIAS4. The expression of C protein suppressed the antiviral proteins expression, and the inhibition effect was enhanced by interaction of PIAS4 and C protein. These results highlighted the beneficial properties of cellular PIAS4 for BVDV protein expression and growth.

Conclusions

This study provides reliable clues for understanding the roles of PIAS4 in the regulation of BVDV growth.

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.