Microarray analysis reveals distinct signaling pathways transcriptionally activated by infection with bovine viral diarrhea virus in different cell types

Tribbles pseudokinase 3 promotes enterovirus A71 infection via dual mechanisms

Enterovirus A71 (EV-A71) is the main pathogen causing hand, foot and mouth disease (HFMD) in children and occasionally associated with neurological diseases such as aseptic meningitis, brainstem encephalitis (BE) and acute flaccid paralysis. We report here that cellular pseudokinase tribbles 3 (TRIB3) facilitates the infection of EV-A71 via dual mechanisms. In one hand, TRIB3 maintains the metabolic stability of scavenger receptor class B member 2 (SCARB2), the bona fide receptor of EV-A71, to enhance the infectious entry and spreading of the virus. On the other hand, TRIB3 facilitates the replication of EV-A71 RNA in a SCARB2-independent manner. The critical role of TRIB3 in EV-A71 infection and pathogenesis was further demonstrated in vivo in mice. In comparison to wild-type C57BL/6 mice, EV-A71 infection in TRIB3 knockdown mice (Trib3+/-) resulted in significantly lower viral loads in muscular tissues and reduced lethality and severity of clinical scores and tissue pathology. In addition, TRIB3 also promoted the replication of coxsackievirus B3 (CVB3) and coxsackievirus A16 (CVA16) in vitro. In conclusion, our results suggest that TRIB3 is one of key host cellular proteins required for the infection and pathogenesis of EV-A71 and some other human enteroviruses and may thus be a potential therapeutic target for combating the infection of those viruses.

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.

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.

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

Background

Bovine viral diarrhea virus (BVDV) is an economically important viral pathogen of domestic and wild ruminants. Apart from cattle, small ruminants (goats and sheep) are also the susceptible hosts for BVDV. BVDV infection could interfere both of the innate and adaptive immunity of the host, while the genes and mechanisms responsible for these effects have not yet been fully understood. Peripheral blood mononuclear cells (PBMCs) play a pivotal role in the immune responses to viral infection, and these cells were the target of BVDV infection. In the present study, the transcriptome of goat peripheral blood mononuclear cells (PBMCs) infected with BVDV-2 was explored by using RNA-Seq technology.

Results

Goat PBMCs were successfully infected by BVDV-2, as determined by RT-PCR and quantitative real-time RT-PCR (qRT-PCR). RNA-Seq analysis results at 12 h post-infection (hpi) revealed 499 differentially expressed genes (DEGs, fold-change ≥ ± 2, p < 0.05) between infected and mock-infected PBMCs. Of these genes, 97 were up-regulated and the remaining 352 genes were down-regulated. The identified DEGs were found to be significantly enriched for locomotion/ localization, immune response, inflammatory response, defense response, regulation of cytokine production, etc., under GO enrichment analysis. Cytokine-cytokine receptor interaction, TNF signaling pathway, chemokine signaling pathway, etc., were found to be significantly enriched in KEGG pathway database. Protein-protein interaction (PPI) network analysis indicated most of the DEGs related to innate or adaptive immune responses, inflammatory response, and cytokine/chemokine-mediated signaling pathway. TNF, IL-6, IL-10, IL-12B, GM-CSF, ICAM1, EDN1, CCL5, CCL20, CXCL10, CCL2, MAPK11, MAPK13, CSF1R and LRRK1 were located in the core of the network and highly connected with other DGEs.

Conclusions

BVDV-2 infection of goat PBMCs causes the transcription changes of a series of DEGs related to host immune responses, including inflammation, defense response, cell locomotion, cytokine/chemokine-mediated signaling, etc. The results will be useful for exploring and further understanding the host responses to BVDV-2 infection in goats.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5830-y) contains supplementary material, which is available to authorized users.

Transcriptomic analysis of responses to cytopathic bovine viral diarrhea virus-1 (BVDV-1) infection in MDBK cells

The bovine viral diarrhea virus (BVDV) is responsible for significant economic losses in the dairy and cattle industry; however, little is known about the protective and pathological responses of hosts to infection. The present study determined the principal molecular markers implicated in viral infection through meta-transcriptomic analysis using MDBK cells infected for two hours with a field isolate of BVDV-1. While several immune regulator genes were induced, genes involved in cell signaling, metabolic processes, development, and integrity were down-regulated, suggesting an isolation of infected cells from cell-to-cell interactions and responses to external signals. Analysis through RT-qPCR confirmed the expression of more than one hundred markers. Interestingly, there was a significant up-regulation of two negative NF-κB regulators, IER3 and TNFAIP3, indicating a possible blocking of this signaling pathway mediated by BVDV-1 infection. Additionally, several genes involved in the metabolism of reactive oxygen species were down-regulated, suggesting increased oxidative stress. Notably, a number of genes involved in cellular growth and development were also regulated during infection, including MTHFD1L, TGIF1, and Brachyury. Moreover, there was an increased expression of the genes β-catenin, caprin-2, GSK3β, and MMP-7, all of which are crucial to the Wnt signaling pathway that is implicated in the embryonic development of a variety of organisms. This meta-transcriptomic analysis provides the first data towards understanding the infection mechanisms of cytopathic BVDV-1 and the putative molecular relationship between viral and host components.

Global transcriptomic profiling of bovine endometrial immune response in vitro. II. Effect of bovine viral diarrhea virus on the endometrial response to lipopolysaccharide

Infection with noncytopathic bovine viral diarrhea virus (ncpBVDV) is associated with uterine disease and infertility. This study investigated the influence of ncpBVDV on immune functions of the bovine endometrium by testing the response to bacterial lipopolysaccharide (LPS). Primary cultures of mixed epithelial and stromal cells were divided into four treatment groups (control [CONT], BVDV, CONT+LPS, and BVDV+LPS) and infected with ncpBVDV for 4 days followed by treatment with LPS for 6 h. Whole-transcriptomic gene expression was measured followed by Ingenuity Pathway Analysis. Differential expression of 184 genes was found between CONT and BVDV treatments, showing interplay between induction and inhibition of responses. Up-regulation of TLR3, complement, and chemotactic and TRIM factors by ncpBVDV all suggested an ongoing immune response to viral infection. Down-regulation of inflammatory cytokines, chemokines, CXCR4, and serine proteinase inhibitors suggested mechanisms by which ncpBVDV may simultaneously counter the host response. Comparison between BVDV+LPS and CONT+LPS treatments showed 218 differentially expressed genes. Canonical pathway analysis identified the key importance of interferon signaling. Top down-regulated genes were RSAD2, ISG15, BST2, MX2, OAS1, USP18, IFIT3, IFI27, SAMD9, IFIT1, and DDX58, whereas TRIM56, C3, and OLFML1 were most up-regulated. Many of these genes are also regulated by IFNT during maternal recognition of pregnancy. Many innate immune genes that typically respond to LPS were inhibited by ncpBVDV, including those involved in pathogen recognition, inflammation, interferon response, chemokines, tissue remodeling, cell migration, and cell death/survival. Infection with ncpBVDV can thus compromise immune function and pregnancy recognition, thereby potentially predisposing infected cows to postpartum bacterial endometritis and reduced fertility.

Cytopathic BVDV-1 strain induces immune marker production in bovine cells through the NF-κB signaling pathway

The bovine viral diarrhea virus (BVDV-1) is a pathogen responsible for high economic losses in the cattle industry worldwide. This virus has the capacity to modulate the immune system of several higher vertebrates, but there is little information available on the cell infection mechanism. To further investigate the effects of BVDV-1 on the activation of the immune response, the Madin-Darby bovine kidney cell line was infected with the cytopathic CH001 field isolate of BVDV-1, and the proinflammatory and antiviral cytokine expression profiles were analyzed. The results showed that BVDV-1 was able to induce the production of BCL3, IL-1β, IL-8, IL-15, IL-18, Mx-1, IRF-1, and IRF-7 in a way similar to polyinosinic-polycytidylic acid. Interestingly, all BVDV-1 activities were blocked by pharmacological inhibitors of the NF-κB signaling pathway. These results, together with in silico analyses showing the presence of several regulatory consensus target motifs, suggest that BVDV-1 regulates gene expression in bovines through the activation of several key transcription factors. Collectively, these data identified BVDV-1 as a viral regulator of immune marker expression, even from early infection. Additionally, this is the first report to find BVDV-1 modulating the activation of cytokine production and transcriptions factors mainly through the NF-κB pathway in vertebrates.

The Molecular Biology of Pestiviruses

Pestiviruses are among the economically most important pathogens of livestock. The biology of these viruses is characterized by unique and interesting features that are both crucial for their success as pathogens and challenging from a scientific point of view. Elucidation of these features at the molecular level has made striking progress during recent years. The analyses revealed that major aspects of pestivirus biology show significant similarity to the biology of human hepatitis C virus (HCV). The detailed molecular analyses conducted for pestiviruses and HCV supported and complemented each other during the last three decades resulting in elucidation of the functions of viral proteins and RNA elements in replication and virus-host interaction. For pestiviruses, the analyses also helped to shed light on the molecular basis of persistent infection, a special strategy these viruses have evolved to be maintained within their host population. The results of these investigations are summarized in this chapter.

Expression of type I interferon-induced antiviral state and pro-apoptosis markers during experimental infection with low or high virulence bovine viral diarrhea virus in beef calves

The objective of this study was to compare the mRNA expression of host genes involved in type-I interferon-induced antiviral state (IFN-α, IFN-β, Mx-1, PKR, OAS-1 and ISG-15), and apoptosis (caspase-3, -8, and -9), after experimental infection of beef calves with low or high virulence noncytopathic (ncp) bovine viral diarrhea virus (BVDV) strains. Thirty BVDV-naïve, clinically normal calves were randomly assigned to three groups. Calves were intranasally inoculated with low (LV; n=10, strain SD-1) or high (HV; n=10, strain 1373) virulence ncp BVDV or BVDV-free cell culture medium (Control, n=10). Quantitative RT-PCR was used to determine the target gene expression in tracheo-bronchial lymph nodes and spleen 5 days after infection. Interferon-α and -β mRNA levels were up-regulated in tracheo-bronchial lymph nodes (P<0.05) in the HV group, but not in the LV group, compared with the control group. There was an up-regulation of type I interferon-induced genes in spleen and tracheo-bronchial lymph nodes of HV and LV groups, compared with the control group (P<0.01). mRNA levels of OAS-1 and ISG-15 were significantly higher in LV than HV calves (P<0.05). A significant up-regulation of caspase-8 and -9 was observed in tracheo-bronchial lymph nodes in the LV group (P=0.01), but not in the HV group. In conclusion, experimental infection with either high or low virulence BVDV strains induced a significant expression of the type I interferon-induced genes in beef calves. There was a differential expression of some interferon-induced genes (OAS-1 and ISG-15) and pro-apoptosis markers based on BVDV virulence and genotype.

A novel biclustering approach with iterative optimization to analyze gene expression data

OBJECTIVE

With the dramatic increase in microarray data, biclustering has become a promising tool for gene expression analysis. Biclustering has been proven to be superior over clustering in identifying multifunctional genes and searching for co-expressed genes under a few specific conditions; that is, a subgroup of all conditions. Biclustering based on a genetic algorithm (GA) has shown better performance than greedy algorithms, but the overlap state for biclusters must be treated more systematically.

RESULTS

We developed a new biclustering algorithm (binary-iterative genetic algorithm [BIGA]), based on an iterative GA, by introducing a novel, ternary-digit chromosome encoding function. BIGA searches for a set of biclusters by iterative binary divisions that allow the overlap state to be explicitly considered. In addition, the average of the Pearson's correlation coefficient was employed to measure the relationship of genes within a bicluster, instead of the mean square residual, the popular classical index. As compared to the six existing algorithms, BIGA found highly correlated biclusters, with large gene coverage and reasonable gene overlap. The gene ontology (GO) enrichment showed that most of the biclusters are significant, with at least one GO term over represented.

CONCLUSION

BIGA is a powerful tool to analyze large amounts of gene expression data, and will facilitate the elucidation of the underlying functional mechanisms in living organisms.

Validation and application of normalization factors for gene expression studies in rubella virus-infected cell lines with quantitative real-time PCR

Reference genes are generally employed in real-time quantitative PCR (RT-qPCR) experiments to normalize variability between different samples. The aim of this study was to identify and validate appropriate reference genes as internal controls for RT-qPCR experiments in rubella virus (RV)-infected Vero and MCF-7 cell lines using SYBR green fluorescence. The software programs geNorm and NormFinder and the DeltaDeltaC(t) calculation were used to determine the expression stability and thus reliability of nine suitable reference genes. HPRT1 and HUEL, and HUEL and TBP were identified to be most suitable for RT-qPCR analysis of RV-infected Vero and MCF-7 cells, respectively. These genes were used as normalizers for transcriptional activity of selected cellular genes. The results confirm previously published microarray and Northern blot data, particularly on the transcriptional activity of the cyclin-dependent kinase inhibitor p21 and the nuclear body protein SP100. Furthermore, the mRNA level of the mitochondrial protein p32 is increased in RV-infected cells. The effect on cellular gene transcription by RV-infection seems to be cell line-specific, but genes of central importance for viral life cycle appear to be altered to a similar degree. This study does not only provide an accurate and flexible tool for the quantitative analysis of gene expression patterns in RV-infected cell lines. It also indicates, that the suitability of a reference gene as normalizer of RT-qPCR data and the host-cell response to RV-infection are strictly cell-line specific.

Bovine viral diarrhea virus non-structural protein 5A interacts with NIK- and IKKbeta-binding protein

Bovine viral diarrhea virus (BVDV) is a positive-sense, single-stranded RNA virus that causes an economically important livestock disease worldwide. Previous studies have suggested that non-structural protein 5A (NS5A) from hepatitis C virus (HCV) and BVDV plays a similar role during virus infection. Extensive reports are available on HCV NS5A and its interactions with the host cellular proteins; however, the role of NS5A during BVDV infection remains largely unclear. To identify the cellular proteins that interact with the N terminus of NS5A and could be involved in its function, we conducted a yeast two-hybrid screening. As a result, we identified a cellular protein termed bovine NIK- and IKKbeta-binding protein (NIBP), which is involved in protein trafficking and nuclear factor kappa B (NF-kappaB) signalling in cells. The interaction of NS5A with NIBP was confirmed both in vitro and in vivo. Complementing our glutathione S-transferase pull-down and immunoprecipitation data are the confocal immunofluorescence results, which indicate that NS5A colocalized with NIBP on the endoplasmic reticulum in the cytoplasm of BVDV-infected cells. Moreover, the minimal residues of NIBP that interact with NS5A were mapped as aa 597-623. In addition, overexpression of NS5A inhibited NF-kappaB activation in HEK293 and LB9.K cells as determined by luciferase reporter-gene assay. We further showed that inhibition of endogenous NIBP by small interfering RNA molecules enhanced virus replication, indicating the importance of NIBP implications in BVDV pathogenesis. Being the first reported interaction between NIBP and a viral protein, this finding suggests a novel mechanism whereby viruses may subvert host-cell machinery for mediating trafficking as well as NF-kappaB signalling.

Activation of extracellular signal-regulated kinase in MDBK cells infected with bovine viral diarrhea virus

Our efforts to identify the cellular signaling cascades triggered by bovine viral diarrhea virus (BVDV) infection in MDBK cells revealed marked activation of extracellular signal-regulated kinase 1/2 (ERK). Enhanced phosphorylation of ERK was detected following infection with cytopathogenic (cp) BVDV, but not with noncytopathogenic BVDV. It appears that cp BVDV-induced ERK phosphorylation is caused by oxidative stress, because ERK phosphorylation was inducible by treatment with hydrogen peroxide or serum deprivation and was attenuated by addition of antioxidants. These results indicate that BVDV infection influences the ERK signaling pathway via oxidative stress, depending on the biotype.