Skip to main content

Involvement of the Wnt pathway in BVDV cytopathogenic strain replication in primary bovine cells



Bovine viral diarrhea virus 1 (BVDV-1) of the pestivirus genus is an economically crippling virus in the cattle industry; this positive RNA virus causes mucosal disease resulting in reproductive losses and other disease syndromes. The pathogenesis mechanism of the disease caused by BVDV infection is not well understood; for a better understanding of in vivo host BVDV-1 interactions, we conducted a transcriptomic study of infected cells at different times post-infection.


We compared the permissiveness and cellular response of a BVDV-1 cytopathogenic strain on Madin-Darby Bovine Kidney cells (MDBK) and bovine lung primary cells, a model closer to in vivo infection. Then a RNAseq analysis was realized on the infected bovine lung primary cells, at 10 hpi and 30 hpi (hours post-infection), to identify transcriptomic signatures.


RNAseq analysis on BVDV-1 infected bovine primary cells showed 2,759 and 5,376 differentially expressed genes at respectively 10 hpi and 30 hpi with an absolute Fold Change  ≥ 2. Among the different pathways deregulated, data analysis revealed a deregulation of Wnt signaling pathway, a conserved process that play a critical role in embryogenesis, cellular proliferation, and differentiation as well as in viral responses against viruses such as Influenza or Hepatitis C. We demonstrated here that the deregulation of the Wnt/βcatenin signaling pathway plays a role in viral replication of BVDV cp strain. Interestingly, we showed that the inhibition of this Wnt pathway using two inhibitors, FZM1 and iCRT14, induced a delay in onset of the establishment of a cytopathic effect of primary cells.


Thereby, this study highlighted a role of the Wnt signaling pathway in the BVDV-1 viral replication in bovine cells, suggesting an interesting option to explore as a new therapeutic target.


Bovine viral diarrhea virus (BVDV), belonging to the Pestivirus genus of the Flaviviridae family, is causative organism of bovine viral diarrhea and mucosal disease, a disease that causes reproductive losses and a range of other syndromes in cattle like congenital malformation and immunosuppression [1]. The virus is contracted from contact between infected cattle through saliva excretion or other contacts with fluids and, following infection, BVDV can infect almost every organ system such as gastrointestinal or respiratory tracts but appears to be dependent on BVDV species [2, 3]. BVDV is a positive-sense RNA virus of 12 kb with a 5’ and a 3’ untranslated region framing one open reading frame which will be translated into only one polyprotein then 12 proteins [1]. Two biotypes of BVDV are recognized: cytopathogenic (cp) or non-cytopathogenic (ncp), a criterion defined by the strain’s capacity to induce a cytopathic effect (CPE) in cell culture. In cattle, both biotypes are involved in mucosal disease, a disease emerging in an animal Persistently Infected (PI) by a ncp strain which is superinfected by a cp strain [3].

As per most viruses, BVDV uses and deregulates the host components and cellular pathways at different steps of the life cycle, such as antiviral pathways, inflammatory response and apoptosis. Several studies report also a deregulation of the Wnt canonical pathway, a conserved process playing a critical role in embryogenesis, cellular proliferation and differentiation,

This Wnt pathway has been recently reported as playing a role in host responses to infection including viral infection [4] both for DNA viruses such as Herpes Simplex virus [5], Epstein Bar virus [6], Cytomegalovirus [7], Hepatitis virus type B (HBV) and RNA viruses with Human Immunodeficient virus [8, 9], Hepatitis virus type C [10], Dengue virus [11], Rift Valley Fever virus [12], Sendai virus [13] or Influenza [14].

Wnt signaling pathway is triggered by WNT ligands binding to Frizzled (FZD) 7-transmembrane domain receptor, low-density lipoprotein-relation protein (LRP5, LRP6) receptors on the cell surface as well as receptor tyrosine kinases ROR and RYK. Depending on this receptor context, WNT ligands then activate distinct intracellular pathways, which can be grouped into β-catenin dependent and β-catenin independent signaling events. Briefly, in the β-catenin dependent pathway, binding of WNT to FDZ receptor and LRP co-receptor promotes the destruction complex stabilization of cytoplasmic β-catenin that then is translocated to nucleus and enables to active transcription factors of the TCF/LEF family. The Wnt response genes are regulated, such as the MYC and PPAR transcription factors involved in cell proliferation and lipid metabolism respectively, or the Cyclin D gene involved in cell cycle [15,16,17]. Among these genes, the receptors FZD and LRP, are generally downregulated while antagonist DKK are overexpressed to establish a feedback control of the Wnt pathway [18,19,20].

In the independent β-catenin signaling pathways, Wnt/Jnk signaling in conjunction with ROR or RYK co-receptors activates via RAC1 and RHOA the transcription factors of the JUN family, whereas Wnt/Ca2+ signaling activates the transcription factors of the NF-AT family [21].

A transcriptomic study of in vitro BVDV-1 infected cells showed an activation of the Wnt/β-catenin signaling pathway at a transcriptional level [22]. However, this study was performed on MDBK cells, an epithelial cell line derived from bovine kidney and did not investigate the effect of the Wnt pathway inhibition in viral replication.

In the present study, to be closer to an in vivo infection, we have chosen to work with primary pulmonary cells from bovine lung, an in vivo target cell type of BVDV [2, 23]. First, a preliminary study was carried out on MDBK and primary cells to evaluate the permissiveness and the cellular response after BVDV-1 NADL infection. Then, we performed by RNAseq, at 10 hpi and 30 hpi (hours post-infection), a transcriptomic analysis of lung bovine cells infected with the cytopathogenic BVDV-1 NADL. A complementary study has showed the activation of the Wnt pathway after infection of primary cells with cytopathogenic strain compared to non-cytopathogenic strain [24]. Here, data analysis confirmed a significant activation of the Wnt signaling pathway and demonstrated an implication of the β-catenin dependent pathway. To investigate the role of this Wnt pathway by BVDV-1 we used two inhibitors, iCRT14 and FZM1, and revealed both a reduction of viral production and a delay in apparition of CPE, suggesting a key role of the Wnt/β-catenin signaling pathway in the viral replication and pathogenicity of the BVDV-1 virus in mucosal disease.


Cells and viruses

Madin-Darby Bovine Kidney cells (MDBK) were cultured in a rich glucose medium (BI medium) supplemented with 5% of fetal bovine serum (FBS) in a 150 cm2 flask at 37 °C with 5% CO2.

Several bovine lungs were collected from healthy animals (Centre de Recherche de Saint Vulbas, BI France). Several pieces were placed in a GentleMACS tubes (Miltenyi Biotec) for grinding. The lysates were then filtered on a 40 µm filter and centrifugated at 700 g during 10 min. The cell pellets were resuspended in a rich glucose medium (BI medium) supplemented with 10% of FBS. The obtained bovine lung primary cells (BPC) were pooled from 4 different lungs then cultured and conserved in nitrogen liquid with 10% DMSO.

The strain BVDV-1 NADL cp (AJ133738) were amplified internally on bovine cells. Viral titers in fluorescent antibodies infectious dose 50% (FAID50/mL) were determined on BPC or MDBK after 3 days with a pool of anti-Pestivirus monoclonal antibodies (BI collection) and Goat anti-Mouse IgG Secondary Antibody, Alexa Fluor 488 (Invitrogen). The CPE were observed on BPC or MDBK with ZOE Fluorescent Cell Imager (Biorad).

Kinetics protocol

BPC and MDBK were cultured for 3 days at 37 °C with 5% CO2 in 24 well plates at a starting density of 3 × 105 cells/well. After 3 days, cells were inoculated with BVDV-1 NADL cp at a multiplicity of infection (MOI) of 1 or non-infected. The infection was done in triplicates.

To study cellular responses, cell cultures were harvested for RNA extraction at 0, 6-, 15-, 24-, and 30-h post-infection. At each time, cell lysates were harvested with RLT buffer (Qiagen), and total RNA was extracted using RNeasy Plus Mini Kit (Qiagen) and stored at − 70 °C.

Droplet digital PCR

Collected cells were screened for the detection of the E2 region of BVDV-1 and 7 genes involved in antiviral response. The PCR primers used in this study were designed with UGENE software (Table 1) [25]. Droplet digital PCR was performed on 5µL of RNA extract with One-step RT-ddPCR Advanced Kit for Probes, AutoDG and QX200 system (Biorad). Reaction mixtures were incubated for 60 min at 50 °C, 10 min at 95 °C, followed by 40 cycles of 30 s at 95 °C, 1 min at 60 °C, and finally 10 min at 98 °C. For each sample, housekeeping gene TBP (Tata Binding Protein) was used to normalize.

Table 1 Primers and probes sequences

Transcriptomic study

Primary cells were cultured for 3 days at 37 °C in 24 well plates and at a starting density of 3 × 105 cells/well in triplicate. The medium used is a rich glucose medium (BI) supplemented with 5% of FBS without anti-Pestivirus antibodies. Then, the cells were inoculated at a MOI of 1 with BVDV-1 NADL and harvested at 0–10–30 h post-infection. Total RNA was extracted from cells with RNeasy Plus Mini kit (Qiagen), quantified using the Quantifluor RNA system (Promega) and qualified using a Bioanalyzer 2100 (Agilent). The library was prepared using the NextFlex Rapid Directional RNAseq V2.0 (Perkin Elmer), quantified and qualified using the Quantus Quantification kit (Promega) and the Fragment Analyzer (AATI). Sequencing was performed on a NextSeq 500 Flow Cell High Output SR75 instrument (Illumina) with 12 samples per flow cell (corresponding to 46 million reads/sample).

Differential expression analysis

Trimming step was performed with CutAdapt software and trimmed reads were mapped using kallisto v0.44.0 against Bovine transcriptome (ARS-UCD1.2 from Ensembl). A second analysis was performed with a reads subset against the bovine genome using Tophat v2.1.1, PicardTools v2.1.0 and FeatureCount v1.5.0 software. Then, normalization was performed with DESeq2 R packages v1.26.0 and differential expression was analyzed by Ingenuity Pathway Analyis (Qiagen). Threshold set to a minimum of 100 reads with an absolute FC ≥ 1.5 and a p-value ≤ 0.05.

Wnt pathway inhibition

BPC were plated in 24 wells for 3 days at 37 °C at a starting density of 3 × 105 cells/well in triplicate. Cells were inoculated at a MOI of 1 with BVDV-1 NADL and treated with 25 µM of iCRT14 (Sigma-SML0203) or 10 µM FZM1 (Merck-534358) at 1-, 2- or 3-h post-infection. The supernatants were harvested at 24 h post-infection and titred, in FAID50/mL, on MDBK with a pool of anti-Pestivirus monoclonal antibodies (BI collection) and Goat anti-Mouse IgG Secondary Antibody, Alexa Fluor 488 (Invitrogen).

The CPE were observed on BPC under microscope in white light.


Comparison of BVDV-1 NADL susceptibility and viral replication of BPC cells in comparison to MDBK cells

In order to compare the infection efficiency of BPC cells to BVDV-1 NADL, in comparison to the MDBK cells, 3 × 105 cells were infected at an MOI of 1 by the NADL strain and then, observed at 8, 24, 48, 72 hpi, by microscopy. Results (Fig. 1) showed presence of CPE at similar levels for both infected cell types after 48 h only, confirming the permissiveness of BPC cells for the NADL strain. At 72 h post infection the CPE is total for the two cell types.

Fig. 1
figure 1

Microscopic observation of cytopathic effect on MDBK cells and BPC. Cells were infected at a MOI of 1 and observed with microscope in white light at 48 hpi

To evaluate viral replication, RNAs were extracted from both MDBK and BPC lysates, at different times post infection (0, 6, 15, 24, 30) and quantified by droplet PCR using specific primer of E2 region (See Method paragraph). Results (Fig. 2) showed a similar production of viral RNA in both MDBK and BPC.

Fig. 2
figure 2

BVDV replication in bovine primary cells and MDBK. Viral RNA production follow by droplet digital PCR. The kinetics was done in triplicate. In dark blue corresponds viral RNA in BPC, in light blue viral RNA in MDBK and in green non infected BPC and MDBK cells. All quantifications were normalized with TBP (n = 3)

Comparison of cellular antiviral responses induced by NADL between MDBK and BPC cells

To assess antiviral cellular responses of BPC cells to NADL and compared them with BPC cells, we chose a panel of 7 genes involved in antiviral response TLR2, TLR3, NFKB, IL6, IL8, ISG15 and MX1, and monitored by droplet digital PCR, their expression during time course of infection 0, 6, 15, 24, 30 h post-infection.

Results (Fig. 3) showed an overexpression of IL-6, IL-8, ISG15, MX1 and TLR2 genes, in infected MDBK and BPC cells, but an overexpression only in BPC cells, of the NFKB and TLR3 (104.7 in BPC infected vs. 102.6 copies/µL in uninfected cells) genes suggesting a better antiviral response of the target BPC cells to NADL infection. This is supported by the levels of expression of IL6 (105 copies/µL for BPC vs 103 copies/µL for MDBK) and IL8 (104.7 vs 102.6 copies/µL) genes that were lower in infected MDBK compared to infected BPC.

Fig. 3
figure 3

Cellular gene expression of BPC and MDBK followed by droplet digital PCR. The kinetics was done in triplicate. In dark blue corresponds BPC infected, in light blue MDBK infected and in green non infected BPC and MDBK cells. All quantifications were normalized with TBP (n = 3)

Analysis of the Wnt signaling pathway genes responses after BVDV-1 NADL infection

In order to better understand the host–pathogen interactions, the cellular gene expression induced by BVDV-1 was analyzed by transcriptomic analysis, using RNASeq. For this, total cellular RNA was extracted from infected BVDV-1 NADL (MOI 1) or non-infected cells, at 0, 10, and 30 hpi and then sequenced using Illumina NextSeq 500 technology.

Data analysis revealed 2759 differentially expressed genes at 10 hpi and 5376 at 30 hpi with an absolute Fold Change (FC) ≥ 2 (p-value ≤ 0.05). As expected for this type of study, we observe a homogeneous distribution of upregulated and downregulated genes for both time points (Fig. 4).

Fig. 4
figure 4

Homogeneous distribution of genes deregulation in infected BPC. Genes significantly downregulated are represented by green spots and genes significantly upregulated are represented by red spots

To complete this global analysis, a focus using Ingenuity Pathway Analyis was performed on the Wnt signaling pathway identified previously by PANTHER database and activated after infection with the cytopathogenic biotype [24]. In response to viral infection, WNT ligand are involved in the positive or negative regulation of type I interferon. For example, WNT2B and WNT9B were identified as negative regulators of Sendai virus-inducted interferon-β; WNT5A as positive, and WNT7B, WNT9B as negative regulators of Influenza A replication; WNT2 and WNT3 as positive regulators of Dengue virus but their functions in the context of immune responses are incompletely understood and role of this Wnt pathway in BVDV infection has been yet never described.

We found a down expression of the receptors belonging to the Wnt pathway including FZD 2, 4, 7, 8 genes (FC − 5.48, − 3.34, − 3.31 and − 3.03, respectively at 30 hpi), a down expression of LRP1 and 6 (− 13.25 and − 1.88 at 30 hpi), involved in the feedback control, and an overexpression of DKK 3 and DKKL1 genes (2.81 and 3.30 at 30 hpi), antagonist of the WNT ligands, confirming the involvement of the Wnt pathway (Table 2) (Fig. 5).

Table 2 Transcripts fold change of genes involved in Wnt signaling pathway
Fig. 5
figure 5

Schematic representation of the Wnt canonical pathway. The Wnt canonical pathway is activated after WNT ligand binding to receptors complex Frizzled-LRP5/6. This interaction allows the formation of a complex, the “signalosome”, that blocked β-catenin degradation leading to its accumulation in cytosol and translocation in nucleus. The β-catenin forms a transcription complex with the transcription factor TCF and LEF to activate genes expression. On the OFF state, β-catenin is continually degraded by the destruction complex which promotes β-catenin phosphorylation leading to its degradation by the proteasome. Some interesting genes identified in Table 2 are reported with arrows (red color for upregulated genes and green for downregulated genes)

Moreover, we found a dysregulation of gene expression involved the dependent β-catenin signaling including an overexpression of MYC (1.81) and CCND1 (2.74), Wnt2 (9.54/5.39) and Wnt5A (15.30/9.74) ligands, a downexpression of GSK3A (− 1.54) and GSK3B (− 2.03) involved in the ß-Catenin degradation (Table 2).

Analysis of Wnt signaling pathway inhibition on BVDV replication

In order to understand the role of this Wnt/ß-catenin signaling pathway in BVDV-1 infection, an inhibition of this pathway was performed in infected cells, using two inhibitors of this pathway: the iCRT14 inhibitor that inhibits the of β-catenin-TCF interaction and the FZM1 inhibitor that promotes the phosphorylation of the β-catenin and its degradation. BPC cells were infected at a MOI of 1 with BVDV-1 NADL cp strain and 25 µM of iCRT14 and 10 µM of FZM1 was added at 1-, 2- or 3-h post-infection. Cell viability in the different conditions was evaluated using a Vi-Cell cell counter and viability analyzer and no impact was observed (data not shown).

Then, supernatants were collected from BPC infected cells, at 24 h post infection and then titrated on MDBK cells (Fig. 6). Moreover, BPC infected cells treated 1-h post-infection with inhibitors were observed by microscopy at 48 hpi to evaluate the presence of ECP. Results showed a fivefold diminution of virus titers when the iCRT14 inhibitor is add 2 h post infection (106.4 FAID50/mL vs. 105.8 FAID50/mL for the control) and a tenfold virus reduction with the FZM1 inhibitor (104.3 FAID50/mL compared to 105.46 FAID50/mL for the control) (Fig. 6). For microscopic observations we observed a complete absence of CPE with both FZM1 and iCRT14 inhibitors (Fig. 7).

Fig. 6
figure 6

BVDV NADL viral production on BPC treated with iCRT14. BPC were plated in 24 wells for 3 days then infected with BVDV NADL cp strain at a MOI of 1. FZM1 or iCRT14 was added in culture at 1-, 2- or 3-hours post-infection and cells were harvested at 24 hpi. The viral protein is labeled with a pool of anti-Pestivirus antibodies then revealed with Alexa488 (green). Titration was performed on MDBK in duplicate for each replicate (n=3). Student test p-value * <0.05 and ** <0.01

Fig. 7
figure 7

Microscopic observation of cytopathic effect on BPC. Cells were infected at a MOI of 1 and treated 1 hpi with inhibitor FZM1 or iCRT14 at 10 µM and 25 µM respectively. Cells were observed under microscope in white light at 48 hpi. Titration was performed on MDBK in duplicate for each replicate (n = 3). Student test p-value * < 0.05

In conclusion, altogether of our results demonstrated for the first time, the implication of the Wnt/ß-Catenin pathway in the viral replication of BVDV-1 NADL cp strain.


This present study has the objective to improve our knowledge about the cell response after infection by a cytopathogenic strain of BVDV, in primary bovine pulmonary cells, the in vivo target cells of the virus. Although previous work reports a transcriptomic analysis of cells infected with BVDV1 [2, 23] and under lighted the deregulation of different gene pathways including the Wnt pathway, this study was carried out on MDBK cells, an epithelial cell line derived from bovine kidney, and author did not investigate the role of the Wnt pathway in BVDV replication.

In a first time, we confirmed the interest of working on primary cells. Effectively, if we showed that the permissiveness and viral production of BPC are equivalent to those of MDBK cells, on the over hand, the antiviral response to BVDV is more important in BPC. Indeed, a targeted analysis of gene expression by droplet digital PCR, of 7 genes involved in antiviral response allowed us to show an overexpression of IL-6, IL-8, ISG15, MX1 and TLR2 genes in both MDBK cells and BPC but higher, in BPC cells (up to 100-fold for the IL6 and IL8 genes) and an overexpression of the NFKB and TLR3 genes only in BPC cells.

A previous study has demonstrated the activation of the Wnt pathway depending on the biotype of BVDV [24]. In this study RNAseq analysis confirmed, in BVDV-1 NADL infected BPC cells, an activation of the Wnt/β-catenin pathway with an overexpression of genes coding for WNT ligands (Wnt2 and Wnt5A) and Wnt response genes (MYC, CCND1, DKK3, DKKL1), down regulation of FZD and LRP receptor genes (Fig. 7). The involvement of the Wnt pathway was confirmed using iCRT14 and FZM1 inhibitors, that induced a reduction of viral production and a delay in apparition of CPE, suggesting a key role of the Wnt/β-catenin signaling pathway in the viral replication and production. WNT network in response to infections has driven significant interest. WNT ligands (and other ligands for WNT receptors) are known to contribute to the host control in peculiar, in polarization of the antigen-presenting cells (APCs), orchestrating phagocytosis, and inflammatory cytokine responses; they are also involved in viral infection by both DNA and RNA viruses, including as example Hepatitis C [10], HIV [8, 9], Influenza [14], and Dengue RNA viruses [11].

For HCV, a flavivirus closer to Pestivirus, Liu et al. have shown that the HCV Core protein activates the β-catenin/Tcf-4-dependent pathway and increases active β-catenin expression and upregulates gene expression of canonical WNT ligands and Frizzled receptors [26]. It has been also reported that the HCV NS5A protein activates β-catenin signaling cascades (Street) and that the interaction of NS5 with β-catenin leads to the stabilization of the cellular protein in cytoplasm [10, 27].

Concerning bovine viruses, it has been reported that bovine herpesvirus (HSV)1 stimulates the Wnt/b-catenin pathway [28] and that the iCRT14 inhibitor reduces productive infection [29].

The correlation between the apoptosis process and mucosal disease lesions due to cp strain superinfection has been established previously in BVDV [23]. Mucosal disease is triggered following the superinfection of a PI animal by a cp strain and is always fatal in cattle. So, it could be of interest to study the effect of the Wnt pathway inhibition on an in vivo model presenting mucosal disease. Understanding the functional nature and temporal regulation of WNT responses in the host response to infection to BVDV, is essential for identifying therapeutic opportunities.


The pathogenesis mechanism of mucosal disease caused by BVDV infection is not well understood. Herein, through the identification of the involvement of the Wnt/β catenin pathway, we identified two specific inhibitors that would allow to better understand the pathogeny of the disease, and that could be evaluated as potential therapeutics target to minimize the BVDV load in animal.

Availability of data and material

The datasets analyzed during the current study are available from the corresponding author on reasonable request.



Bovine lung primary cells


Bovine viral diarrhea virus


Cytopathic effects


Madin-Darby bovine kidney cells


Multiplicity of infection


Persistently infected


Hours post-infection


  1. Tautz N, Tews BA, Meyers G. The molecular biology of pestiviruses, In: Advances in virus research. Elsevier; 2015. pp. 47–160.

  2. Spagnuolo-Weaver M, Allan GM, Kennedy S, Foster JC, Adair BM. Distribution of cytopathic and noncytopathic bovine viral diarrhea virus antigens in tissues of calves following acute experimental infection. J VET Diagn Invest. 1997;9:287–97.

    CAS  Article  PubMed  Google Scholar 

  3. Spetter MJ, Uriarte ELL, Altamiranda EAG, Leunda MR, Pereyra SB, Verna AE, Odeón AC. Dual natural infection with bovine viral diarrhea virus -1 and -2 in a stillborn calf: tissue distribution and molecular characterization. Open Vet J. 2018;8:493.

    Article  PubMed  PubMed Central  Google Scholar 

  4. van Zuylen WJ, Rawlinson WD, Ford CE. The Wnt pathway: a key network in cell signalling dysregulated by viruses: Wnt signalling and viral infection. Rev Med Virol. 2016;26:340–55.

    Article  PubMed  Google Scholar 

  5. Zhu L, Jones C. The canonical Wnt/β-catenin signaling pathway stimulates herpes simplex virus 1 productive infection. Virus Res. 2018;256:29–37.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  6. Morrison JA, Klingelhutz AJ, Raab-Traub N. Epstein-Barr virus latent membrane protein 2A activates beta-catenin signaling in epithelial cells. J Virol. 2003;77:12276–84.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  7. Angelova M, Zwezdaryk K, Ferris M, Shan B, Morris CA, Sullivan DE. Human cytomegalovirus infection dysregulates the canonical Wnt/β-catenin signaling pathway. PLoS Pathog. 2012;8: e1002959.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Butler JS, Dunning EC, Murray DW, Doran PP, O’Byrne JM. HIV-1 protein induced modulation of primary human osteoblast differentiation and function via a Wnt/β-catenin-dependent mechanism. J Orthop Res. 2013;31:218–26.

    CAS  Article  PubMed  Google Scholar 

  9. Weiser K, Barton M, Gershoony D, Dasgupta R, Cardozo T. HIV’s Nef interacts with β-catenin of the Wnt signaling pathway in HEK293 cells. PLoS ONE. 2013;8: e77865.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  10. Park C-Y, Choi S-H, Kang S-M, Kang J-I, Ahn B-Y, Kim H, Jung G, Choi K-Y, Hwang SB. Nonstructural 5A protein activates β-catenin signaling cascades: Implication of hepatitis C virus-induced liver pathogenesis. J Hepatol. 2009;51:853–64.

    CAS  Article  PubMed  Google Scholar 

  11. Smith JL, Jeng S, McWeeney SK, Hirsch AJ. A MicroRNA screen identifies the wnt signaling pathway as a regulator of the interferon response during flavivirus infection. J Virol. 2017;91:e02388-e2416.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  12. Harmon B, Bird SW, Schudel BR, Hatch AV, Rasley A, Negrete OA. A genome-wide RNA interference screen identifies a role for Wnt/β-catenin signaling during rift valley fever virus infection. J Virol. 2016;90:7084–97.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  13. Baril M, Es-Saad S, Chatel-Chaix L, Fink K, Pham T, Raymond V-A, Audette K, Guenier A-S, Duchaine J, Servant M, Bilodeau M, Cohen E, Grandvaux N, Lamarre D. Genome-wide RNAi screen reveals a new role of a WNT/CTNNB1 signaling pathway as negative regulator of virus-induced innate immune responses. PLoS Pathog. 2013;9: e1003416.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  14. More S, Yang X, Zhu Z, Bamunuarachchi G, Guo Y, Huang C, Bailey K, Metcalf JP, Liu L. Regulation of influenza virus replication by Wnt/β-catenin signaling. PLoS ONE. 2018;13: e0191010.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  15. He TC, Chan TA, Vogelstein B, Kinzler KW. PPARdelta is an APC-regulated target of nonsteroidal anti-inflammatory drugs. Cell. 1999;99:335–45.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  16. He TC, Sparks AB, Rago C, Hermeking H, Zawel L, da Costa LT, Morin PJ, Vogelstein B, Kinzler KW. Identification of c-MYC as a target of the APC pathway. Science. 1998;281:1509–12.

    CAS  Article  PubMed  Google Scholar 

  17. Shtutman M, Zhurinsky J, Simcha I, Albanese C, D’Amico M, Pestell R, Ben-Ze’ev A. The cyclin D1 gene is a target of the beta-catenin/LEF-1 pathway. Proc Natl Acad Sci USA. 1999;96:5522–7.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  18. González-Sancho JM, Aguilera O, García JM, Pendás-Franco N, Peña C, Cal S, García de Herreros A, Bonilla F, Muñoz A. The Wnt antagonist DICKKOPF-1 gene is a downstream target of beta-catenin/TCF and is downregulated in human colon cancer. Oncogene. 2005;24:1098–103.

    CAS  Article  PubMed  Google Scholar 

  19. Müller HA, Samanta R, Wieschaus E. Wingless signaling in the Drosophila embryo: zygotic requirements and the role of the frizzled genes. Development. 1999;126:577–86.

    Article  Google Scholar 

  20. Niida A, Hiroko T, Kasai M, Furukawa Y, Nakamura Y, Suzuki Y, Sugano S, Akiyama T. DKK1, a negative regulator of Wnt signaling, is a target of the beta-catenin/TCF pathway. Oncogene. 2004;23:8520–6.

    CAS  Article  PubMed  Google Scholar 

  21. Ljungberg JK, Kling JC, Tran TT, Blumenthal A. Functions of the WNT signaling network in shaping host responses to infection. Front Immunol. 2019;10:2521.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  22. Villalba M, Fredericksen F, Otth C, Olavarría V. Transcriptomic analysis of responses to cytopathic bovine viral diarrhea virus-1 (BVDV-1) infection in MDBK cells. Mol Immunol. 2016;71:192–202.

    CAS  Article  PubMed  Google Scholar 

  23. Hilbe M, Girao V, Bachofen C, Schweizer M, Zlinszky K, Ehrensperger F. Apoptosis in Bovine viral diarrhea virus (BVDV)-induced mucosal disease lesions: a histological, immunohistological, and virological investigation. Vet Pathol. 2013;50:46–55.

    CAS  Article  PubMed  Google Scholar 

  24. La Polla R, Testard MC, Goumaidi A, Chapot E, Legras-Lachuer C, de Saint-Vis B. Identification of differentially expressed gene pathways between cytopathogenic and non-cytopathogenic BVDV-1 strains by analysis of the transcriptome of infected primary bovine cells. Virology. 2021;18(567):34–46. (Epub ahead of print).

    CAS  Article  Google Scholar 

  25. Okonechnikov K, Golosova O, Fursov M, UGENE team. Unipro UGENE: a unified bioinformatics toolkit. Bioinformatics. 2012;28:1166–1167.

  26. Liu J, Wang Z, Tang J, Tang R, Shan X, Zhang W, Chen Q, Zhou F, Chen K, Huang A, Tang N. Hepatitis C virus core protein activates Wnt/β-catenin signaling through multiple regulation of upstream molecules in the SMMC-7721 cell line. Arch Virol. 2011;156:1013–23.

    CAS  Article  PubMed  Google Scholar 

  27. Street A, Macdonald A, McCormick C, Harris M. Hepatitis C virus NS5A-mediated activation of phosphoinositide 3-kinase results in stabilization of cellular β-catenin and stimulation of β-catenin-responsive transcription. J Virol. 2005;79:5006.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  28. Liu Y, Hancock M, Workman A, Doster A, Jones C. β-Catenin, a transcription factor activated by canonical Wnt signaling, is expressed in sensory neurons of calves latently infected with bovine herpesvirus 1. J Virol. 2016;90:3148–59.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  29. Zhu L, Thunuguntla P, Liu Y, Hancock M, Jones C. The β-catenin signaling pathway stimulates bovine herpesvirus 1 productive infection. Virology. 2017;500:91–5.

    CAS  Article  PubMed  Google Scholar 

Download references


The authors thank Philippe PICHON and CRSV technicians for the lung sampling and Stéphanie BLANCHARD, Philippe CHEROT, Sylvie MBELO-NGONKANGA for the technical support.


This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors.

Author information

Authors and Affiliations



RLP was involved in the conceptualization, methodology, writing—original draft, formal analysis, and methodology. BSV contributed to writing—review and editing, funding acquisition, methodology, and supervision. CLL contributed to resources and writing—review and editing. MCT helped in the investigation, data collection and analysis. AG was involved in the data collection and curation and formal analysis. OG was involved in analysis and methodology. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Blandine de Saint-Vis.

Ethics declarations

Ethical approval

The animal sampling conducted for this study was reviewed by internal Ethics Committee and conformed to the relevant regulatory standards defined by European and French Laws.

Consent for publication & competing interest

All authors have approved the manuscript for submission and have no conflicts of interest to disclose.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

La Polla, R., Testard, MC., Garcia, O. et al. Involvement of the Wnt pathway in BVDV cytopathogenic strain replication in primary bovine cells. Virol J 19, 134 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • BVDV-1
  • Cytopathogenic
  • Transcriptomics
  • Wnt pathway
  • Viral replication