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

Rémi La Polla, Marie-Claire Testard, Océane Garcia, Abdelghafar Goumaidi, Catherine Legras-Lachuer, Blandine de Saint-Vis

► To cite this version:

Rémi La Polla, Marie-Claire Testard, Océane Garcia, Abdelghafar Goumaidi, Catherine Legras-Lachuer, et al.. Involvement of the Wnt pathway in BVDV cytopathogenic strain replication in primary bovine cells. Virology Journal, 2022, 19 (1), pp.134. 10.1186/s12985-022-01863-6. hal-03788382

HAL Id: hal-03788382

https://hal.inrae.fr/hal-03788382

Submitted on 26 Sep 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

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

Rémi La Polla<fnref n="1" />,², Marie-Claire Testard<fnref n="1" />, Océane Garcia<fnref n="1" />, Abdelghafar Goumaidi³, Catherine Legras-Lachuer²† and Blandine de Saint-Vis<fnref n="1" />* ID

**Abstract**

**Background:** 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.

**Methods:** 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.

**Results:** 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.

**Conclusions:** 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.

**Keywords:** BVDV-1, Cytopathogenic, Transcriptomics, Wnt pathway, Viral replication

**Introduction**

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

<footnote>†Catherine Legras-Lachuer and Blandine D. E. Saint-Vis equally contributed to this work.</footnote>

<footnote>*Correspondence: blandine.de_saint-vis@boehringer-ingelheim.com</footnote>

<footnote n="1">¹ Site Lyon porte des Alpes, Boehringer Ingelheim Animal Health, 813 cours du 3eme Millenaire, 69800 Saint Priest, France
Full list of author information is available at the end of the article</footnote>
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–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–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/Ca²⁺ 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 permisiveness 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.

## Methods

### 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 cm² flask at 37 °C with 5% CO₂.

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 × 10⁵ 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.

## Transcriptomic study

Primary cells were cultured for 3 days at 37 °C in 24 well plates and at a starting density of 3 × 10⁵ 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 Analysis (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 × 10⁵ 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.

<table><caption>**Table 1** Primers and probes sequences</caption><thead><tr><td>Target (Accession number)</td><td>Foward (5'-&gt; 3')</td><td>Reverse (5'-&gt; 3')</td><td>Probe (5'-&gt; 3')</td></tr></thead><tbody><tr><td>BVDV NADL (E2 region)</td><td>CCAAAGATGCACGAGAGAAACC</td><td>CTCTTGCTTTCGCCCATCAA</td><td>TGCATACAAGAGCCTTGCCGACCAGTGT</td></tr><tr><td>TBP (NM0010745742)</td><td>GCGTTTTGCTGCTGTAATCA</td><td>GGAACTTCACATCACAGCTC</td><td>CGCACCACTGCACTGATATTCAGTTCT</td></tr><tr><td>TLR2 (KT601038)</td><td>GTGCTTTGTGGACAGCATGG</td><td>TTGGACAGGTCAAGGCTTTTCA</td><td>CAGGCTTCTTCTCTGTCTTGTGACCC</td></tr><tr><td>TLR3 (AJ812026)</td><td>AAAGGGACTTTGAGGCAGGT</td><td>GATTTAAACATYCCTCTTCG</td><td>ATGCAGTTCAGCAAGCTATTGAACAA</td></tr><tr><td>NFKB (NM_001076409)</td><td>CTGGAAGCACGAATGACAGA</td><td>AACATGAGCCGTACCACGCT</td><td>CTATAATCCCGGACTTTT GGT GC A</td></tr><tr><td>IL6 (EU276071)</td><td>CAAAATCTCTGCAATGAGAAAGG</td><td>GCTTCAGGATCTGGATCAGTG</td><td>GATGCTTCCAATCTGGGTTCAATC</td></tr><tr><td>IL8 (AF232704)</td><td>ATGACTTCCAAGCTGGCTGTT</td><td>ACAACCTTCTGCACCCACTTT</td><td>CATTCCACACCTTTCCACCCCAAATTTATC</td></tr><tr><td>Mxl (AF047692)</td><td>TTCCACCTGAAGAAGGGCTA</td><td>GCCATACTTCTGTAACTCCTCTGTT</td><td>ATCCCCCTGCCTGGCAGAAAGACT</td></tr><tr><td>ISG15 (MI 74,366)</td><td>TGCAGAACTGC ATCTC CATC</td><td>CATGAACACGGTGCACCCCT</td><td>AGACCAGTTCTGGCTGTCTTTTGAAG</td></tr></tbody></table>
# Results

## 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 × 10⁵ cells were infected at an MOI of 1 by the NADL strain and then, observed at 8, 24, 48, 72 hpi, by microcopy. 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.

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.

## 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* (10⁴.⁷ in BPC infected vs. 10².⁶ 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* (10⁵ copies/µL for BPC vs 10³ copies/µL for MDBK) and *IL8* (10⁴.⁷ vs 10².⁶ copies/µL) genes that were lower in infected MDBK compared to infected BPC.

## 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

<figure><img src="image_3.png" /><figcaption>**Fig. 1** Microscopic observation of cytopathic effect on MDBK cells and BPC. Cells were infected at an MOI of 1 and observed with microscope in white light at 48 hpi</figcaption></figure>
<figure><img src="image_4.png" /><figcaption>**Fig. 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)</figcaption></figure>

type of study, we observe a homogeneous distribution of upregulated and downregulated genes for both time points (Fig. 4).

To complete this global analysis, a focus using Ingenuity Pathway Analysis 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).

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).
<figure><img src="image_5.png" /><figcaption>**Fig. 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)</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>**Fig. 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</figcaption></figure>
<table><caption>**Table 2** Transcripts fold change of genes involved in Wnt signaling pathway</caption><thead><tr><th>Wnt pathway genes</th><th>NADL 10 hpi vs NADL 0 hpi</th><th>p-value</th><th>NADL 30 hpi vs NADL 0 hpi</th><th>p-value</th><th>Wnt pathway genes</th><th>NADL 10 hpi vs NADL 0 hpi</th><th>p-value</th><th>NADL 30 hpi vs NADL 0 hpi</th><th>p-value</th></tr></thead><tbody><tr><td>ACVR2A</td><td>-2,65</td><td>2,64E-14</td><td>-2,67</td><td>5,26E-48</td><td>PIN1</td><td></td><td></td><td>1,88</td><td>4,22E-18</td></tr><tr><td>AKT3</td><td>1,58</td><td>1,96E-12</td><td></td><td></td><td>PPARD</td><td>-2,99</td><td>8,25E-05</td><td></td><td></td></tr><tr><td>APC</td><td>-2,82</td><td>1,08E-03</td><td></td><td></td><td>PPP2CB</td><td></td><td></td><td>-2,14</td><td>3,89E-37</td></tr><tr><td>APPL1</td><td>-1,61</td><td>1,69E-34</td><td>-1,78</td><td>2,57E-49</td><td>PPP2R3A</td><td></td><td></td><td>-1,85</td><td>5,67E-18</td></tr><tr><td>BMPR2</td><td>-1,59</td><td>1,56E-10</td><td>-1,94</td><td>4,89E-10</td><td>PPP2R1B</td><td>1,54</td><td>2,21E-38</td><td></td><td></td></tr><tr><td>CCND1</td><td>2,74</td><td>2,75E-77</td><td></td><td></td><td>PPP2R5B</td><td></td><td></td><td>2,00</td><td>4,78E-32</td></tr><tr><td>CDH1</td><td>-6,84</td><td>5,86E-69</td><td>-4,21</td><td>3,50E-42</td><td>PPP2R5D</td><td></td><td></td><td>-1,53</td><td>2,75E-29</td></tr><tr><td>CDH2</td><td></td><td></td><td>-2,77</td><td>6,28E-17</td><td>RARA</td><td></td><td></td><td>-5,35</td><td>3,69E-41</td></tr><tr><td>CREBBP</td><td>-4,94</td><td>4,12E-07</td><td></td><td></td><td>RARG</td><td>-1,99</td><td>3,19E-06</td><td>-2,61</td><td>9,86E-05</td></tr><tr><td>CSNK2A1</td><td></td><td></td><td>-1,54</td><td>7,54E-38</td><td>RPS27A</td><td></td><td></td><td>2,96</td><td>5,38E-09</td></tr><tr><td>CSNK1G2</td><td>1,51</td><td>4,98E-22</td><td></td><td></td><td>SFRP1</td><td>2,19</td><td>2,72E-19</td><td>5,25</td><td>8,56E-87</td></tr><tr><td>DKK3</td><td></td><td></td><td>2,81</td><td>7,50E-05</td><td>SFRP2</td><td>-2,26</td><td>8,62E-43</td><td>-10,10</td><td>0</td></tr><tr><td>DKKL1</td><td>1,61</td><td>3,99E-06</td><td>3,30</td><td>8,23E-35</td><td>SFRP4</td><td></td><td></td><td>2,32</td><td>1,41E-103</td></tr><tr><td>DVL3</td><td>-3,07</td><td>1,32E-22</td><td>-3,49</td><td>5,77E-28</td><td>SOX4</td><td></td><td></td><td>1,68</td><td>1,15E-17</td></tr><tr><td>FZD2</td><td>-2,89</td><td>1,46E-82</td><td>-5,48</td><td>3,86E-173</td><td>TCF3</td><td></td><td></td><td>-1,92</td><td>8,41E-15</td></tr><tr><td>FZD3</td><td></td><td></td><td>2,02</td><td>5,08E-15</td><td>TCF4</td><td></td><td></td><td>-2,16</td><td>3,62E-15</td></tr><tr><td>FZD4</td><td>-2,11</td><td>4,86E-26</td><td>-3,34</td><td>1,64E-61</td><td>TGFB1</td><td></td><td></td><td>-1,58</td><td>5,90E-03</td></tr><tr><td>FZD7</td><td>-3,06</td><td>3,24E-179</td><td>-3,31</td><td>2,84E-201</td><td>TGFB2</td><td>-1,69</td><td>1,98E-07</td><td></td><td></td></tr><tr><td>FZD8</td><td></td><td></td><td>-3,03</td><td>1,30E-18</td><td>TGFB3</td><td>-1,57</td><td>6,87E-12</td><td>-1,85</td><td>3,02E-21</td></tr><tr><td>GJA1</td><td>1,63</td><td>3,24E-34</td><td></td><td></td><td>TGFBR2</td><td></td><td></td><td>2,11</td><td>1,66E-02</td></tr><tr><td>GNAQ</td><td>-1,69</td><td>2,08E-30</td><td>-2,49</td><td>7,50E-84</td><td>TGFBR3</td><td>-2,17</td><td>9,53E-24</td><td></td><td></td></tr><tr><td>GSK3A</td><td></td><td></td><td>-1,54</td><td>8,17E-08</td><td>TLE3</td><td></td><td></td><td>-7,38</td><td>3,71E-35</td></tr><tr><td>GSK3B</td><td></td><td></td><td>-2,03</td><td>4,51E-26</td><td>UBA52</td><td></td><td></td><td>1,52</td><td>3,87E-27</td></tr><tr><td>KREMEN1</td><td>-1,70</td><td>1,51E-05</td><td>-3,56</td><td>2,17E-25</td><td>Ubb</td><td>1,63</td><td>6,58E-145</td><td>2,84</td><td>0</td></tr><tr><td>LRP1</td><td>-5,59</td><td>2,82E-05</td><td>-13,25</td><td>5,63E-11</td><td>UBC</td><td>-2,93</td><td>3,38E-13</td><td>-2,42</td><td>1,81E-09</td></tr><tr><td>LRP6</td><td></td><td></td><td>-1,88</td><td>1,71E-05</td><td>WNT2</td><td>9,54</td><td>2,02E-131</td><td>5,39</td><td>2,92E-71</td></tr><tr><td>MARK2</td><td></td><td></td><td>-2,13</td><td>1,83E-10</td><td>WNT2B</td><td>-2,56</td><td>7,03E-69</td><td>-3,40</td><td>1,90E-112</td></tr><tr><td>MDM2</td><td>2,65</td><td>6,90E-04</td><td>8,78</td><td>3,10E-16</td><td>WNT5A</td><td>15,30</td><td>0</td><td>9,74</td><td>0</td></tr><tr><td>MYC</td><td>1,81</td><td>1,44E-19</td><td></td><td></td><td>WNT5B</td><td>2,46</td><td>1,01E-120</td><td>1,56</td><td>4,74E-29</td></tr><tr><td>NLK</td><td>1,98</td><td>7,53E-09</td><td>1,86</td><td>1,05E-07</td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>

Wnt signaling pathway is studied after cell infection by NADL at 10 hpi or 30 hpi versus 0 hpi. The genes upregulated are colored in red and the genes downregulated are colored in green. A threshold of 100 for reads counts and ≥ 1.5 for fold change is set.

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 microcopy at 48 hpi to evaluate the presence of ECP. Results showed a fivefold diminution of virus titers when the iCRT14 inhibitor is added 2 h post infection (10<sup>6.4</sup> FAID50/mL vs. 10<sup>5.8</sup> FAID50/mL for the control) and a tenfold virus reduction with the FZM1 inhibitor (10<sup>4.3</sup> FAID<sub>50</sub>/mL compared to 10<sup>5.46</sup> 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).

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.

## Discussion

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 missiveness and viral production of BPC are equivalent to those of MDBK cells, on the over hand, the antiviral
<figure><img src="image_7.png" /></figure>

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
<figure><img src="image_8.png" /><figcaption>**Fig. 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</figcaption></figure>

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.

## Conclusion

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.
<figure><img src="image_9.png" /><figcaption>**Fig. 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</figcaption></figure>

**Abbreviations**

BPC: Bovine lung primary cells; BVDV: Bovine viral diarrhea virus; CPE: Cytopathic effects; MDBK: Madin-Darby bovine kidney cells; MOI: Multiplicity of infection; PI: Persistently infected; hpi: Hours post-infection.

**Acknowledgements**

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

**Author contributions**

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.

**Funding**

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

**Availability of data and material**

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

**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.

**Author details**

¹Site Lyon porte des Alpes, Boehringer Ingelheim Animal Health, 813 cours du 3eme Millenaire, 69800 Saint Priest, France. ²Laboratoire d'Écologie Microbienne - UMR 5557, Université Claude Bernard Lyon 1, 43 Boulevard du 11 Novembre 1918, 69622 Villeurbanne Cedex, France. ³Faculté de Médecine et de Pharmacie, Viroscan3D, 8 Avenue Rockefeller, 69373 Lyon, France.
Received: 23 February 2022 Accepted: 2 August 2022

Published online: 19 August 2022

## References

1. Tautz N, Tews BA, Meyers G. The molecular biology of pestiviruses, In: Advances in virus research. Elsevier; 2015. pp. 47–160. https://doi.org/10.1016/bs.aivir.2015.03.002

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. https://doi.org/10.1177/104063879700900310.

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. https://doi.org/10.4314/ovj.v8i4.23.

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. https://doi.org/10.1002/rmv.1892.

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

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. https://doi.org/10.1128/jvi.77.22.12276-12284.2003.

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. https://doi.org/10.1371/journal.ppat.1002959.

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. https://doi.org/10.1002/jor.22196.

9. Weiser K, Barton M, Gershooony D, Dasgupta R, Cardozo T. HIV's Nef interacts with β-catenin of the Wnt signaling pathway in HEK293 cells. PLoS ONE. 2013;8: e77865. https://doi.org/10.1371/journal.pone.0077865.

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 for hepatitis C virus-induced liver pathogenesis. J Hepatol. 2009;51:853–64. https://doi.org/10.1126/jhep.2009.06.026.

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. https://doi.org/10.1128/JVI.02388-16.

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. https://doi.org/10.1128/JVI.00543-16.

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. https://doi.org/10.1371/journal.ppat.1003416.

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. https://doi.org/10.1371/journal.pone.0191010.

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. https://doi.org/10.1016/s0092-8674(00)81664-5.

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. https://doi.org/10.1126/science.281.5382.1509.

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. https://doi.org/10.1073/pnas.96.10.5522.

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. https://doi.org/10.1038/sj.onc.1208303.

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.

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. https://doi.org/10.1038/sj.onc.1207892.

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. https://doi.org/10.3389/fimmu.2019.02521.

22. Villalba M, Fredericksen F, Otth C, Olavarria V. Transcriptomic analysis of responses to cytopathic bovine viral diarrhea virus-1 (BVDV-1) infection in MDBK cells. Mol Immunol. 2016;71:192–202. https://doi.org/10.1016/j.molimm.2016.01.009.

23. Hilbe M, Girao V, Bachofen C, Schweizer M, Ziinszky 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. https://doi.org/10.1177/0300985812447826.

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. https://doi.org/10.1016/j.virol.2021.12.005 (Epub ahead of print).

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

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. https://doi.org/10.1007/s00705-011-0943-x.

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. https://doi.org/10.1128/JVI.79.8.5006-5016. 2005.

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. https://doi.org/10.1128/JVI.02971-15.

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. https://doi.org/10.1016/j.virol.2016.10.014.

## Publisher's Note

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

Ready to submit your research? Choose BMC and benefit from:

* fast, convenient online submission

* thorough peer review by experienced researchers in your field

* rapid publication on acceptance

* support for research data, including large and complex data types

* gold Open Access which fosters wider collaboration and increased citations

* maximum visibility for your research: over 100M website views per year

At BMC, research is always in progress.

Learn more biomedcentral.com/submissions

<figure><img src="image_10.png" /></figure>