Preprint
Article

This version is not peer-reviewed.

Porcine Circovirus Modulates Swine Influenza Virus Replication in Pig Tracheal Epithelial Cells and Porcine Alveolar Macrophages

A peer-reviewed article of this preprint also exists.

Submitted:

22 April 2023

Posted:

23 April 2023

You are already at the latest version

Abstract
The pathogenesis of porcine circovirus type 2b (PCV2b) and swine influenza A virus (SwIV) during co-infection in swine respiratory cells is poorly understood. To elucidate the impact of PCV2b/SwIV co-infection, newborn porcine tracheal epithelial cells (NPTr) and immortalized porcine alveolar macrophages (iPAM 3D4/21) were co-infected with PCV2b and SwIV (H1N1 or H3N2 genotype). Viral replication, cell viability and cytokine mRNA expression were determined and compared between single-infected and co-infected cells. Finally, 3’mRNA sequencing was performed to identify the modulation of gene expression and cellular pathways in co-infected cells. It was found that PCV2b significantly decreased and improved SwIV replication, in co-infected NPTr and iPAM 3D4/21 cells respectively, compared to single infected cells. Interestingly, PCV2b/SwIV co-infection synergistically up-regulated IFN expression in NPTr cells whereas in iPAM 3D4/21 cells, PCV2b impaired the SwIV IFN induced response, both correlating with SwIV replication modulation. RNA-sequencing analyses revealed that the modulation of gene expression and enriched cellular pathways during PCV2b/SwIV H1N1 co-infection is regulated in a cell type-dependent-manner. This study revealed different outcomes of PCV2b/SwIV co-infection in porcine epithelial cells and macrophages and provides new insights on porcine viral co-infections pathogenesis.
Keywords: 
;  ;  ;  ;  ;  ;  

1. Introduction

The porcine respiratory disease complex (PRDC) is cause of enormous economic loss to the swine industry worldwide [1]. Porcine reproductive and respiratory syndrome virus (PRRSV), porcine circovirus type 2 (PCV2) and swine influenza A virus (SwIV) are the main viruses associated to PRDC and are involved into costly outbreaks on swine production sites [2].
PCV2 is a non-enveloped single-stranded circular DNA virus that belongs to the Circoviridae viral family, genus Circovirus [3]. PCV2 is the etiological agent of several diseases and syndromes such as post-weaning multisystemic wasting syndrome (PMWS), porcine dermatitis and nephropathy syndrome (PDNS), among others which are collectively named porcine circovirus-associated disease (PCVAD) [4]. The prevalence of PCV2 infection is high in swine farms all over the world [5]. However, most of PCV2 infections are reported to be subclinical [5]. It is noteworthy to mention that the virus can alter the innate immune response and cause immunosuppression, which favors co- infection and/or the secondary infection with other pathogens [6]. Eight genotypes of PCV2 have been described so far (PCV2a-PCV2h) [7,8,9]. The emergence of PCV2b genotype in Canada was associated with a significant death rate increase during the post-weaning multisystemic wasting syndrome (PMWS) outbreak in late 2004 [10] and was reported as the predominant genotype identified during the same period in samples submitted to the Diagnostic service of the Faculty of Veterinary Medicine (Faculté de médecine vétérinaire) of the University of Montreal [10,11].
Swine influenza A viruses (SwIV are enveloped single stranded negative sense RNA viruses belonging to the Orthomyxoviridae viral family, genus Alphainfluenzavirus [3]. The subtypes of influenza A virus (IAV) are determined by the hemagglutinin and neuraminidase proteins which are embedded into the envelope of the virion. Eighteen HA subtypes (H1–H18) and eleven NA subtypes (N1–N11) have been reported so far [12,13]. SwIV subtypes H1N1, H1N2, H3N2, and A(H1N1)pdm09 have been the most frequently influenza A virus subtypes reported in swine worldwide [14]. In Canada, SwIV subtypes H1N1, H1N2 and H3N2 have been circulating in swine herds since at least 1980’s [15,16,17,18,19,20,21]. In 2009, A(H1N1)pdm09 was reported in swine in canadian provinces Manitoba, Alberta, Saskatchewan and Quebec [22].
Infections caused by SwIV are characterized by high morbidity and very low mortality [23]. It is known that IAV induces high levels of cytokine release very often associated to viral loads and severity of the disease [24] .
PCV2 co-infections with SwIV or PRRSV have been extensively studied in vivo [25,26,27,28,29,30,31,32,33] and in vitro [34,35,36,37]. Even though PCV2 and SwIV co-infections are also prevalent in pigs [1] , pathogenesis studies concerning this co-infection are scarce. Interestingly, an in vivo epidemiologic assessment showed that PCV2-positive pigs were more likely to be infected with SwIV than PCV2-negative pigs [38]. Other authors have reported that PCV2/SwIV H1N1 co-infection did not affect PCV2 replication in experimentally infected pigs, whereas PCV2 infection did increase SwIV-related clinical disease in the infected animals [39].
PCV targets mostly cells of monocyte/macrophage lineage in vivo but can infect other cell types [40,41]. SwIV preferentially targets epithelial cells of the respiratory tract, but alveolar macrophages can also be infected [42,43,44]. Epithelial cells detect viral pathogens via pattern recognition receptors (PRRs) and release anti/proinflammatory cytokines to recruit and activate innate immune cells which finally trigger innate and adaptive immune responses [45,46]. On the other hand, the importance of macrophages for host defense against invading pathogens is well known [47]. They are professional phagocytic cells that clear infectious particles and apoptotic cells and participate in the adaptive immune response by T cells activation via antigen presentation [47]. It is well known that concomitant infections modulate viral pathogenesis and affects host immune cells responses [1]. However, there is a lack of information regarding the effect of PCV2b/SwIV co-infections at the cellular and molecular level. The objective of the present study was to evaluate the effects of PCV2b/SwIV co-infections on virus targeted host cells such as swine respiratory epithelial cells and porcine alveolar macrophages. The results of the present study revealed that PCV2b modulates SwIV replication during co-infection, meanwhile affecting SwIV mRNA cytokine expression and cellular genes modulation in infected cells, in a cell-type dependent manner.

2. Materials and Methods

2.1. Cells

Newborn pig tracheal epithelial cell line (NPTr) and immortalized porcine alveolar macrophage (iPAM 3D4/21) cell line were used for all PCV2b and SwIV single and co-infections. Madin-Darby Canine Kidney (MDCK) cell line was used for SwIV titration and propagation. The NPTr cell line was kindly provided by Dr. M. Ferrari (Instituto Zooprofilattico Sperimental, Brescia, Italy) [48]. The NPTr and the MDCK (ATCC CCL-34) cells were cultured in Eagle's Minimum Essential Medium (EMEM) (Wisent Bioproducts, QC, Canada) supplemented with 10% fetal bovine serum (FBS) (Wisent Bioproducts, QC, Canada), 1 mM sodium pyruvate, 10 I.U./mL of penicillin, 10 μg/mL of streptomycin and 250 g/L amphotericin B solution (Wisent Bioproducts, QC, Canada) [48]. The iPAM 3D4/21 cell line (ATCC CRL-2843) was maintained in RPMI 1640 medium (Invitrogen Corporation, GibcoBRL, ON, Canada) with 10% FBS and adjusted to contain 2 mM L-glutamine (Invitrogen), 10 mM HEPES (Invitrogen), 1 mM sodium pyruvate, 1 mM nonessential amino acids (Invitrogen) and 0.1 mg/ml streptomycin/100U penicillin solution (Invitrogen). All cells were cultivated in a humidified incubator in 5% CO2 atmosphere at 37°C.

2.2. Viruses

The PCV2b strain (FMV-06-0732) was isolated from a clinical case in Quebec in 2006 [10]. It was serially propagated in NPTr cells then purified and concentrated following ultracentrifugation on a 30% sucrose cushion using the SW28 Beckman Coulter rotor (Beckman Coulter Canada Inc., Mississauga, ON, Canada) at 25000 rpm for 4 hours.
The SwIV H1N1 (A/swine/St-Hyacinthe/148/1990) was isolated from pigs with respiratory symptoms during a 1990/91outbreak of respiratory disease in Quebec [16]. The SwIV H3N2 (A/swine/Quebec/1708732/2015(H3N2) was isolated from a pig with acute respiratory disease. The sample was submitted by veterinary practitioners to Diagnostic Veterinary Virology Laboratory (DVVL) of the University of Montreal in 2015 as part of routine diagnostic testing. Viral titers of PCV2b and SwIV stocks were determined in NPTr cells and MDCK, respectively, by the Spearman-Kärber method [49,50]. Viral titers were expressed in tissue culture infectious dose 50% per mL (TCID50/mL).

2.3. Immunofluorescence assay

Single infection and co-infection were confirmed using an immunofluorescence assay (IFA) as previously described [51]. Briefly, the cell medium was removed and infected cells were fixed with a mixture of acetone-methanol (50/50, v/v) and incubated for a 20 minute period. Thereafter, fixed infected cells were washed three times with a PBS solution and then permeabilized with a solution containing 0.1% triton X-100 in PBS and incubated for 10 minutes. After incubation with a blocking solution (1% bovine serum in PBS-Tween for 20 minutes), cells were incubated with a 1/200 dilution of the polyclonal PCV2 porcine antiserum [51] at 37°C for 90 minutes and/or with 1/200 diluted monoclonal mouse anti-NP SwIV primary antibody (Bio-Rad). Then, the cells were washed three times with a solution containing 1% bovine serum in PBS-T and incubated with a 1/75 dilution of a goat anti-swine rhodamine conjugated secondary antibody (Jackson ImmunoResearch, PA, USA) and/or a 1/200 dilution of goat anti-mouse FITC conjugated secondary antibody (Invitrogen, CA, USA), at 37°C for 60 minutes. After three washing steps, the cells were visualized using a Leica DMI 4000 inverted widefield fluorescence microscope (Leica Microsystems Inc, Richmond Hill, Canada). Pictures were acquired with a DFC 490 digital camera (Leica Microsystems Inc.) and images were analyzed using Leica Application Suite Software, version 2.4.0 (Leica Microsystems Inc.) and ImageJ software (Laboratory for Optical and Computational Instrumentation, LOCI, University of Wisconsin) [52].

2.4. PCV2b/SwIV co-infection

To obtain PCV2b/SwIV co-infected cells, NPTr or iPAM 3D4/21 cells were firstly infected with PCV2b at a MOI of 0.05. Afterwards, PCV2b-infected cells were passaged at least three times in the presence of the virus to achieve a stable persistent infection. PCV2b-infected cells were thereafter dispensed into 24-well plates (1 x 105 cells/well) and then co-infected with SwIV H1N1 or H3N2 subtypes at a MOI of 1 for 1 hour in the presence of trypsine (1ug/mL). Cells were then washed three times with PBS. Fresh medium with 2 % FBS was added to the wells and plates were incubated in a humidified incubator with 5% CO2 at 37°C. After 4, 12, 24, 48 and 72 hours SwIV post-infection, cells were subjected to 2-3 rounds of freeze-thaw cycles to release virus particles and cell debris which were then removed by centrifugation at 8000 rpm at 4oC for 15 minutes. The supernatant was kept at -80°C until virus titer determination. SwIV titers were determined by the Spearman- Kärber method [49,50] in MDCK cells and expressed as tissue culture infectious dose 50% per mL (TCID50/mL). PCV2b quantification was performed by qPCR assay as described by Gagnon et al. (2008), using the primers and probe presented in Supplemental Table S1. PCV2b concentrations were expressed as DNA copy amount per mL.

2.5. Cell viability assay

Cell viability assay was performed at 24 hours PCV2b/SwIV post-infection with the Celltiter 96 Aqueous One Solution Cell Proliferation Assay kit (Promega, Madison, WI, USA), according to manufacturer’s instruction. Briefly, 20 µL of the reagent were added to each well of the microplate, the cells were then incubated for 2 hours at 37 °C and the absorbance was measured at 490 nm (Biotek® Synergy HT plaque reader, Vermont, USA ). The infected cells viability percentages were calculated using the non-infected cells as control. Experiments were done in triplicate and repeated at least two times.

2.6. Cytokine mRNAs expression in PCV2b and/or SwIV infected cells

The modulation of mRNAs expression of the following interleukins (ILs), IL6, IL8, IL10 and interferons (IFNs), (IFN)-α, IFN-β and IFN- γ in co-infected cells versus single-infected cells was assessed by RT-qPCR. Total cellular RNA was extracted from infected and mock-infected NPTr or iPAM 3D4/21 cells and purified with the RNeasy Mini Kit (Qiagen), according to the manufacturer’s instructions. Total RNA concentration was measured with a Qubit fluorometer (Thermo Fisher, Walthman, MA, USA). One μg of total cellular RNA was reverse transcribed using M-MLV reverse transcriptase (Invitrogen, Burlington, ON, CA), according to the manufacturer’s protocol. After the reverse transcription reaction, the cDNA was used in qPCR reactions with PowerTrack SYBR Green Master Mix kit (Thermo Fisher Scientific) on a QuantStudio 3 Real-Time PCR System (Thermo Fisher Scientific). Quantification of differences between groups was done using the 2−ΔΔCt method. β2-microglobulin (B2M), β-actin (ACTB) and peptidylprolyl isomerase A (PPIA) were used as normalizing genes to compensate for potential differences in cDNA amounts. Mock-infected cells were used as a calibrator reference in the analysis. The primers used for the specific amplification of the targeted cDNA are described in Supplemental Table S1.

2.6. 3’ mRNA-seq Library Preparation and Sequencing

Total RNA was purified for RNA-seq analysis from infected cells as described in the previous section. RNA concentration and quality were assessed using Agilent 2100 Bioanalyzer apparatus with the RNA 6000 Nano Kit (Agilent Technologies, CA, USA). Only RNA samples with acceptable RNA integrity number (RIN>7) were submitted to library preparation. The cDNA libraries were constructed using Lexogen’s QuantSeqTM 3’ mRNA-Seq Kit (Lexogen GmbH, Vienna, Austria) according to the manufacturer’s recommendations. A total of 12 libraries; including 3 libraries from PCV2b-infected cells, 3 from SwIV H1N1-infected cells, 3 from PCV2b/SwIV H1N1-infected cells and 3 from mock infected cells, were synthesized. Concentration and quality of the purified libraries were assessed using the Qubit fluorometer (Thermo Fisher, Walthman, MA, USA) and an Agilent high-sensitivity DNA kit with a Bioanalyzer (Agilent, CA, USA), respectively. All libraries were sequenced with a MiSeq high throughput apparatus (Illumina, CA, USA) using v3 cartridge (150-cycles) according to manufacturer’s instructions.

2.7. Bioinformatic analysis

Raw sequencing data was imported in FASTQ format into CLC Genomics Workbench (version 22.0.1, Qiagen, CA, USA). Reads were trimmed for quality and adaptors using the CLC Genomics Workbench software and trimmed reads were mapped to Sus scrofa rRNA 12S, 16S and 18S. After removing the rRNA mapped reads the remaining unmapped reads were mapped to the reference genome (Sscrofa11.1). The mapped reads were then used for differential expression calculations. Differentially expressed genes (DEGs) with a False Discovery Rate (FDR) threshold < 0.05 and a fold change >1.5 were selected and used for volcano plots and venn diagrams visualization in CLC Genomics workbench. Pathway enrichment analysis was performed with the Database for Annotation, Visualization and Integrated Discovery (DAVID, 2021) [53,54]. Additionally, the list of DEGs from the co-infected cells were analyzed with ClueGo application (version 2.5.5) available in Cytoscape (version 3.9.1) to identify protein-protein interaction networks and associated cellular pathways [41,42,43,44]. All sequences were deposited in GEO repository with de accession number GSE229215.

2.8. Statistical analyses

All statistical analyses were performed using GraphPad Prism software (GraphPad Prism 7.0.0). The different statistical tests being used within each analysis, are described in figure legends. Differences were considered significant at p-value <0.05.

3. Results

3.1. Modulation of swine influenza A virus replication in co-infected PCV2b/SwIV H1N1 cells

The immunofluorescence assay (IFA) performed on co-infected NPTr and iPAM 3D4/21 cells confirmed the simultaneous presence of PCV2b and SwIV H1N1 antigens in the cells (Figure 1). Similar results were observed in PCV2b/SwIV H3N2 co-infected cells (data not shown). PCV2b internalized with similar efficiency in NPTr and iPAM 3D4/21 cells as the percentage of PCV2b positive cells by IFA was similar for both types of cells. However, the number of SwIV H1N1 positive cells was significantly higher in NPTr infected cells than in iPAM 3D4/21 infected cells (p < 0.001).
The infectious titer of SwIV H1N1/H3N2 subtypes and the PCV2b viral genome quantity in co-infected cells was determined and compared to single infected cells at different time post-infection (Figure 2). The SwIV H1N1 and H3N2 infectious titer at 24 hours post-infection (hpi) was significantly higher (p < 0.0001) in single infected NPTr cells compared to single infected iPAM 3D4/21 cells. These results were expected considering that SwIV preferentially infects epithelial cells of the respiratory tract [42]. However, PCV2b qPCR quantification results were similar in both NPTr and iPAM 3D4/21 infected cells (Figure 3), confirming the IFA results (Figure 1).
During co-infection, it was found that PCV2b modulated the replication kinetics of both subtypes of SwIV, H1N1 and H3N2, in the infected cells. In fact, the SwIV infectious titer decreased significantly after 24 hpi in PCV2b/SwIV co-infected NPTr cells compared to cells infected only with SwIV (Figure 2A and 2C). Inversely, in iPAM 3D4/21 infected cells the infectious titers of SwIV H1N1 and H3N2 were significantly higher in co-infected cells versus single-infected cells (Figure 2B and 2D). However, the presence of SwIV did not seems to influence the PCV2b replication during co-infection. (Figure 3).

3.2. Cell viability assay

To evaluate if the modulation of SwIV replication during co-infection could affect cell viability compared to single infection, a cell viability assay was performed at different times post-infection in single infected and co-infected NPTr and iPAM 3D4/21 cells. As expected, the cell viability decreased from 24 hpi to 72 hpi in co-infected as well as in single infected cells compared to mock infected cells (Figure 4). Moreover, the dual infection of PCV2b and SwIV in both NPTr and iPAM 3D4/21 cells significantly decreased the cell viability over time compared to SwIV single infected cells. Nevertheless, no significant differences were found between PCV2b co-infected and single infected cells.

3.3. Modulation of cytokines mRNAs expression in co-infected cells

The levels of mRNA expression of pro-inflammatory (IL-6 and IL-8) and anti-inflammatory (IL-10) cytokines, as well as of type I IFN (IFN-α/β) and type II IFN (IFN- γ) were measured at 24 hours post-infection in single-infected and co-infected cells. The results obtained with NPTr infected cells revealed that IL-6, IL-8, and IL-10 mRNA expressions were higher in PCV2b/SwIV H1N1 and PCV2b/SwIV H3N2 co-infected cells compared to SwIV H1N1 and SwIV H3N2 single infected cells, respectively (Figure 5A). However, no difference was found for mRNA expression of interleukins tested between PCV2b/SwIV H1N1 and PCV2b infected cells. In PCV2b/SwIV H3N2 NPTr infected cells, IL-6 and IL-10, but not IL-8, mRNA expression was reduced compared to PCV2b single infected cells (Figure 5A).
In the case of IFN type I response, a clear synergistic effect was observed in co-infected NPTr cells as IFN-α and IFN-β mRNAs expression was significantly higher in co-infected cells compared to in PCV2b or SwIV (H1N1 or H3N2) single-infected cells (Figure 5C). Similar results were obtained for IFN-γ in PCV2b/SwIV H1N1 co-infected cells compared to single infected cells. However, although significant differences were observed in the mRNA expression of IFN- γ between PCV2b/SwIV H3N2 and SwIV H3N2 infections, no difference was found when compared to PCV2b single infected cells (Figure 5C).
As shown in Figure 5B, in iPAM 3D4/21 cells, the co-infection with PCV2b and SwIV H1N1 did not significantly impact the IL-6 mRNA expression. However, in co-infected PCV2b/SwIV H3N2 cells, IL-6 expression was reduced compared to SwIV H3N2 infected cells, but not when comparing to PCV2b single infected cells. IL-10 expression decreased in PCV2b/SwIV H1N1 co-infected iPAM 3D4/21 cells compared to PCV2b or SwIV H1N1 infected cells. The dual infection of PCV2b and SwIV H3N2 increased the expression of IL-10 regarding to SwIV H3N2 alone, however, the transcriptional level of this cytokine was significantly reduced in these infections compared to PCV2b single infected cells. Regarding the IFN response in co-infected iPAM 3D4/21 cells (Figure 5D), no difference was observed on the mRNA expression of the IFNs tested compared to PCV2b infected cells. Nevertheless, IFN-β and IFN- γ mRNA expressions were significantly higher in SwIV (H1N1 or H3N2) single infected cells compared to PCV2b/SwIV co-infected cells. Also, IFN-α mRNA expression was higher in SwIV H3N2 infection than in co-infected cells. It is interesting to note that the cell type influenced the modulation of IFN mRNA expression. For example, the mRNA expression levels of IFN-α, IFN-β and IFN- γ were up-regulated in NPTr cells following PCV2b infection (Figure 5C) compared to mock infected cells, but in PCV2b infected iPAM 3D4/21, there was a trend toward down-regulation of those mRNAs(Figure 5D). This difference in IFNs transcriptional response was also observed between co-infected NPTr and iPAM 3D4/21 cells (Figure 5C and Figure 5D). Overall, these results at least revealed that dual PCV2b and SwIV infection modulates cytokine transcriptional responses in a cell-type dependent-manner.

3.4. Differentially Expressed Genes (DEGs) and pathway enrichment analysis in PCV2b, SwIV H1N1 and PCV2b/SwIV H1N1 infected cells

RNA-seq analysis was performed to characterize and compare the transcriptomic response in NPTr and iPAM 3D4/21 co-infected and single infected cells. The transcriptomic analyses were performed at 24 hpi to ensure having a significant amount of suitable viable cells following viral infection. SwIV H1N1 was chosen for the present experiment because during SwIV H1N1/PCV2b co-infection, the modulation of SwIV replication was greater at 24 hpi compared to SwIV H3N2/PCV2b co-infection (Figure 2).
First, 3′ mRNA sequencing was used to identify the DEGs in infected cells compared to mock-infected cells. In the NPTr cells, 365, 98 and 627 DEGs were identified in PCV2b, SwIV H1N1 and PCV2b/SwIV H1N1 infected cells, respectively, when comparing to mock-infected cells. (Figure 6A). The number of up-regulated and down-regulated genes in each infection are listed in Supplemental Table S2 and S3. In iPAM 3D4/21 cells, 316, 64 and 164 DEGs were identified in PCV2b, SwIV H1N1 and PCV2b/SwIV H1N1 infected cells, respectively, compared to mock-infected cells. (Figure 6B). The number of up-regulated and down-regulated genes in each infection are shown in Supplemental Table S4 and S5.
The identified DEGs were used to perform pathway enrichment analysis to determine the cellular pathways impacted during viral infection in single and co-infected cells. After the enrichment analysis, the top 10 over-represented pathways having an FDR < 0.05 were retained. Only pathways having a fold enrichment > 2 are illustrated in infected NPTr cells (Figure 6). In the case of iPAM 3D4/21 cells, the number of enriched pathways having the same FDR ≥ 2 was extremely high. That is why only pathways having a fold enrichment > 4 are depicted in iPAM 3D4/21 cells.
The pathway enrichment analysis in NPTr cells revealed the following cellular pathways: ECM-receptor interaction, protein processing in endoplasmic reticulum, focal adhesion and regulation of actin cytoskeleton were among the most enriched pathways in PCV2b infected cells as well as in co-infected cells (Figure 7A, E7C). Other over-represented pathways in NPTr co-infected cells included mRNA surveillance, HIF-1 signaling pathway, spliceosome, phagosome, NOD-like receptor signaling pathway and influenza A pathway (Figure 7E). In SwIV H1N1 NPTr infected cells, most enriched pathways were related to the protein synthesis machinery and included pathways like: translation, ribosome, metabolism of proteins, among others involved in similar functions (Figure 7C). Also, Nonsense Mediated Decay (NMD) pathways were identified in SwIV H1N1 infected NPTr cells. The pathway enrichment analysis performed on iPAM 3D4/21 cells showed that in single infected as well as in co-infected cells, NF-Kappa B and TNF signaling pathways were the most enriched pathways in the set of DEGs analyzed. However, in co-infected cells, additional enriched cellular pathways were found, these included apoptosis, influenza A, HIF-1 signaling pathway and MAPK signaling pathways (Figure 7B, D7D, 7E and 7F).
An additional analysis was performed in Cytoscape’s ClueGo application with the list of DEGs from the co-infected cells to identify protein-protein interaction networks and to confirm the cellular pathways modulated during PCV2b/SwIV H1N1 co-infection. This analysis confirmed the enrichment of cellular pathways such as mRNA surveillance, spliceosome, protein processing in endoplasmic reticulum, phagosome and regulation of actin cytoskeleton in NPTr co-infected cells (Supplemental Figure S1). In co-infected iPAM 3D4/21 cells, HIF-1 signaling pathway, TNF signaling, regulation of innate immune response, IL-17 signaling, adaptive immune response based on somatic recombination of immune receptors built from immunoglobulin superfamily domains, apoptosis and influenza A pathways were identified (Supplemental Figure S2). The Supplemental Figures S2 and S3 illustrate several DEGs that are involved in more than one cellular pathway in co-infected NPTr and iPAM 3D4/21 cells, respectively.

4. Discussion

In the present study, PCV2b replication was not modulated in the presence of both SwIV strains (H1N1 or H3N2) and in both infected cell lines being tested (Figure 3). These results are in accordance with previous studies. For example, Wei et al, 2010 concluded that infection of Caesarean-derived colostrum-deprived pigs with SIV H1N1, one week after a previous infection with PCV2, did not influence PCV2 replication in dually infected pigs, while PCV2 infection increases SwIV-related clinical disease [39]. Unfortunately, this previous study did not provide any information regarding the effect of PCV2/SwIV co-infection on SwIV replication. Interestingly, our results demonstrated that PCV2b modulated SwIV replication in a cell-type dependent-manner (Figure 2). In fact, PCV2b decreased the replication of SwIV in NPTr cells whereas the SwIV titer was enhanced in iPAM 3D4/21 co-infected cells. Productive replication of SwIV in macrophages can alter antiviral macrophage functions like phagocytosis, resulting in enhanced disease severity and impaired bacterial clearance during secondary bacterial infection [47,55].
PCV2b has been previously reported to induce up-regulation of cytokines such as IL-6, IL-8 and IL-10 in vitro [56,57,58]. In the present study, IL-6 and IL-10 mRNA expressions were increased by PCV2b infection in NPTr and iPAM 3D4/21 cells in accordance with previous studies reporting up-regulation of IL-6 in epithelial cells [58] and increased levels of IL-6 and IL-10 in porcine alveolar macrophages (PAMs) [59,60]. In the present study, IL-8 mRNA was upregulated in PCV2b NPTr infected cells, as expected, whereas it was not in PCV2b iPAM 3D4/21 infected cells (Figure 5). These results do not correlate with previous reports showing a significant up-regulation of IL-8 mRNA expression and/or protein production in PCV2-inoculated swine alveolar macrophages [61,62]. The contradictory results obtained in the present study maybe due to the time points selected to measure the expression of this cytokine in macrophages after PCV2b infection or the use of an immortalized PAM cell line (3D4/21) instead of primary porcine alveolar macrophages. It is noteworthy to mention that cytokine mRNA expression was determined in persistently PCV2b-infected cells, to mimic the superinfection events that occurs in PCV2b previously infected pigs. Therefore, a persistently PCV2b infection in iPAM 3D4/21 cells could have a differential effect on mRNA expression compared to an acute PCV2b cell infection.
Induction of pro-inflammatory cytokines during swine influenza virusSwIV infection has been correlated with viral replication and clinical signs [63]. In the current study, IL-6, IL-8, and IL-10 mRNA expressions were higher in PCV2b/SwIV (H1N1 or H3N2) co-infected NPTr cells compared to SwIV H1N1 and SwIV H3N2 single infected cells. However, no difference was observed regarding the tested interleukin mRNA expression levels between PCV2b/SwIV H1N1 and PCV2b infected cells. Nevertheless, in PCV2b/SwIV H3N2 infected cells, IL-6, and even more IL-10, mRNA expression levels were significantly reduced compared to PCV2b single infected NPTr cells (Figure 5A). These results reveal that the modulation of interleukin expression by PCV2b/SwIV co-infection on NPTr cells might be influenced by the SwIV genotype being investigated.
Regarding the mRNA expression of antiviral cytokines such as IFNs, their expression level modulations in co-infected cells was found to be cell-type dependent. In NPTr cells, the IFN type I (α/β) mRNA expression levels were synergistically up-regulated in co-infected cells compared to single infected cells, whereas in iPAM 3D4/21 cells, a trend towards a reduced IFN type I and II mRNA expression was observed in PCV2b single infected and co-infected cells (Figure 5). Several researchers have previously reported contradictory data in regards to the modulation of IFN response following infection with PCV2b. Wang et al. (2022) reported that PCV2 infection interferes with the activation of type I IFNs signaling pathway and inhibits the IFN-induced ISGs expression in vivo and in vitro [64]. Others have reported the inhibition of IFN-α in porcine peripheral blood mononuclear cells (PBMCs) [65,66] and the inhibition of type I IFN (IFN-α /IFN-β) in PK-15 cells [67,68] following PCV2 infection or in the presence of PCV2 viral protein or DNA. In addition, according to Gao et al. (2014), PCV2 significantly inhibited pseudorabies virus (PRV)-induced IFN-γ mRNA expression in swine PBMC in vitro [69]. However, other authors have reported that porcine circovirus type 2 induces type I interferon production in porcine alveolar macrophages [70] and IFN-β in PK-15 infected cells [71]. Interestingly, it is known that PCV2 genome contains CpG motifs with both IFN-α inhibitory and stimulatory properties [72]. Kekarainen et al. (2008b) suggested that PCV2 viral elements can distinctly regulate cytokine production according to the cell population [73]. Overall, it is easy to conclude that IFN response modulation following PCV2 infection is complex and could vary depending on the experimental models being studied (in the present study: tracheal epithelial cells versus macrophages). It is known that type I and III IFNs are rapidly induced after influenza virus sensing by pattern recognition receptors (Toll-like receptors, RIG-I-like receptors, NOD-like receptors) on respiratory epithelial cells. The released IFNs bind to their cognate receptors and activates the JAK-STAT signaling pathway. This triggers the transcription of hundreds of IFN-stimulatory genes (ISGs) that induce an “antiviral state” in infected and nearby cells to ultimately restrict viral replication and propagation [74,75,76,77,78]. In a study performed by Wu et al. (2022), they overexpressed microRNA let-7 in A549 cells which resulted in an increase of interferon type I mRNA expression and consequently inhibited influenza virus infection [79]. Moreover, Fong et al. (2022) demonstrated that IFN-γ inhibits influenza virus replication in respiratory epithelial cells by reducing viral binding onto the cells [80]. The results obtained during the present study suggest that the increased expression of IFN mRNA in PCV2b/SwIV infected NPTr cells could lead to a reduction of SwIV viral replication compared to single SwIV infected cells. It is also well known that influenza virus possesses several strategies to counteract IFN response [81,82]. That would explain why the increased IFN mRNA expression in co-infected NPTr cells could be responsible for a partial reduction, and not a complete shutdown, of SwIV replication (Figure 2). Interestingly, Czerkies et al. (2022) have reported interference between respiratory syncytial virus (RSV) and influenza A virus H1N1 in human alveolar epithelial cells (A549). The authors of the report found that previous infection of cells with RSV does not prevent a subsequent influenza virus infection, whereas RSV indeed protects bystander cells against influenza virus infection by triggering secretion of type I and type III IFNs [83]. A similar effect in PCV2b/SwIV co-infected NPTr cells could be involved and explain the PCV2b impact on the apparent inhibition of SwIV replication. However, this PCV2b bystander phenomenon effect on SwIV replication is highly plausible but still needs to be demonstrated. In iPAM 3D4/21 co-infected cells, PCV2b completely impaired the modulation of the transcriptional expression of IFN-β and IFN-γ by SwIV. Gao et al. (2014) have reported a similar PCV2 effect on IFN expression during dual infection with pseudorabies virus (PRV) in PBMC [69]. In this previous study, PCV2 significantly inhibited the ability of inactivated PRV to induce IFN-γ expression. The authors suggested that PCV2 could affect the cell immune response to PRV [69]. In the context of the present study, it is believed that the impaired IFN response observed in PCV2b persistently infected iPAM 3D4/21 cells could favor SwIV infection and replication.
Differential transcriptome analysis of PCV2b/SwIV H1N1 and single infected porcine cells revealed distinct signatures on host gene expression that may have impacted SwIV H1N1 replication during the co-infection. Several of the top enriched pathways identified in co-infected cells in this study are known to be involved in the host antiviral response and influenza virus pathogenesis and propagation [84] such as NOD-like receptor signaling pathway in co-infected NPTr cells and apoptosis and MAPK signaling pathways in co-infected iPAM 3D4/21 cells (Figure 7). NOD-like receptors (NLR) is a specific family of cytosolic pattern-recognition receptors (PRRs) that contains more than 20 members in mammals and plays a pivotal role in the recognition of intracellular ligands [85]. A group of NLRs is associated with the formation of the inflammasome, a cytosolic multiprotein complex that, once activated, induce the production of IL-1β and IL-18 and the recruitment of immune cells to the site of infection [86,87]. The NLRP3 inflammasome has an important role in the immune response against influenza A virus infection. Its role in both protective and detrimental immune responses during influenza A virus infection has been studied and extensively reviewed [87,88,89,90]. Interestingly, NLRP3 was found to be up regulated in PCV2b/SwIV H1N1 infected NPTr cells (FC 6.16, FDR =0.04) but not in SwIV H1N1 and PCV2b single infected cells. Moreover, DDX3X, a host protein which has been implicated in the coordination of host defense against influenza virus by activating the NLRP3 inflammasome and type I interferon response [91,92] was upregulated in co-infected NPTr cells (FC 3, FDR =5.84e-4) and in PCV2b infected cells (FC 2.42, FDR = 0.01), but not in SwIV H1N1-infected cells.
Apoptosis and MAPK pathways can be modulated during the replication cycle of influenza virus to favor viral infection. Wurzer et al. (2003) have found that caspase 3 activation during the onset of apoptosis is important for efficient influenza virus propagation [93]. The MAPK pathway, on the other hand, has been shown to positively regulate multiple steps of influenza virus replication such as RNA synthesis, vRNP export and release (budding) of virus, clathrin-independent endocytosis and viral internalization, nuclear export and protein synthesis, (see the review of Kumar et al, 2018 [94] for information and references regarding the role of MAPK signaling in the replication of several virus, including influenza virus).
HIF-1 signaling pathway was identified among the most modulated pathways in the co-infected NPTr and iPAM 3D4/21 cells but not in PCV2b or SwIV H1N1 single infected cells (Figure 7). HIF-1-related pathways have been reported to be involved in viral infections and the innate and the adaptive immune responses of the host against viruses [95]. Previous reports on HIF-1 role into influenza virus pathogenesis are contradictory. It was reported that the H1N1 virus activates HIF-1 pathway and then induces glycolysis to promote viral replication whereas knockdown of this pathway significantly reduced H1N1 replication in A549 cells [96]. However, during another study it was demonstrated that HIF-1α knockdown in A549 cells promoted influenza A virus replication by promoting autophagy in cells [97].

5. Conclusions

In conclusion, it was found that PCV2b decreases SwIV replication in porcine tracheal epithelial cells while enhancing SwIV replication in alveolar porcine macrophages, probably through modulation of IFN expression. In addition, PCV2b and SwIV were found to synergistically enhance IFNs mRNA expression in NPTr infected cells whereas in iPAM 3D4/21 cells, PCV2b impaired the capacity of SwIV to promote IFNs mRNA expression. Cellular genes and pathways were found to be differentially modulated in PCV2b/SwIV H1N1 co-infected cells compared to PCV2b and SwIV H1N1 single-infected cells. However, further studies are needed to establish the role of identified DEGs and cellular pathways into the pathogenesis of PCV2b/SwIV co-infection.

Supplementary Materials

The following supporting information can be downloaded at the website of this paper posted on Preprints.org. Figure S1: Protein-protein interaction networks and associated cellular pathways on PCV2b/swH1N1 co-infected NPTr cells; Figure S2: Protein-protein interaction networks and associated cellular pathways on PCV2b/swH1N1 co-infected iPAM 3D4/21 cells; Table S1: List of primers and probe used for PCV2b quantification and primers used in RT-qPCR for cytokines mRNA quantification in NPTr and iPAM 3D4/21 cells; Table S2: Up-regulated genes in NPTr cells; Table S3: Down-regulated genes in NPTr cells; Table S4: Up-regulated genes in iPAM 3D4/21 cells; Table S5: Down-regulated genes in iPAM 3D4/21 cells.

Author Contributions

Yaima Burgher: Development of all the protocols, acquisition, analysis and interpretation of all the results, writing original draft and revision of the manuscript. Dr Chantale Provost: Project and protocols design, optimization of the library preparation and sequencing, revision of the manuscript. Marie-Jeanne Pesant: Optimization of the library preparation and sequencing, revision of the manuscript. Dr. Carl A. Gagnon: Design of the project and the experiments, funding acquisition, interpretation of the results, revision of the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This project was funded by a Natural Sciences and Engineering Research Council of Canada (NSERC) discovery grant (number RGPIN-2017-05240). C.A. Gagnon is financially supported by a Canadian Swine Research and Development Cluster (CSRDC) research project, Les éleveurs de porcs du Québec and by Probiotech International (number 1781). Y. Burgher was a recipient of a graduate student scholarship from the Fonds de recherche du Québec – Nature et technologie (FRQNT). Y. Burgher and M.J. Pesant were also recipients of graduate student scholarships from the Swine and poultry infectious diseases research center (Centre de recherche en infectiologie porcine et avicole, CRIPA-FRQ) and of different graduate student merit scholarships for Graduate and Postdoctoral Studies from the Université de Montréal.

Acknowledgments

The authors would like to acknowledge the support of current and former post-graduate students of Dr Gagnon’s Laboratories, Faculté de médecine vétérinaire, Université de Montréal.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Opriessnig, T.; Gimenez-Lirola, L.G.; Halbur, P.G. Polymicrobial respiratory disease in pigs. Anim Health Res Rev 2011, 12, 133–148. [Google Scholar] [CrossRef] [PubMed]
  2. Brockmeier, S.L.; Halbur, P.G.; Thacker, E.L. Porcine respiratory disease complex. Wiley Online Library, 2002: p. 231-258.
  3. Lefkowitz, E.J.; Dempsey, D.M.; Hendrickson, R.C.; Orton, R.J.; Siddell, S.G.; Smith, D.B.J.N.a.r. Virus taxonomy: the database of the International Committee on Taxonomy of Viruses (ICTV). 2018, 46, D708–D717.
  4. Chae, C. A review of porcine circovirus 2-associated syndromes and diseases. Vet J 2005, 169, 326–336. [Google Scholar] [CrossRef] [PubMed]
  5. Segales, J. Porcine circovirus type 2 (PCV2) infections: clinical signs, pathology and laboratory diagnosis. Virus Res 2012, 164, 10–19. [Google Scholar] [CrossRef] [PubMed]
  6. Shi, R.; Hou, L.; Liu, J. Host immune response to infection with porcine circoviruses. Animal Diseases 2021, 1, 1–10. [Google Scholar] [CrossRef]
  7. Franzo, G.; Segales, J. Porcine circovirus 2 (PCV-2) genotype update and proposal of a new genotyping methodology. PLoS One 2018, 13, e0208585. [Google Scholar]
  8. Wang, Y.; Noll, L.; Lu, N.; Porter, E.; Stoy, C.; Zheng, W.; Liu, X.; Peddireddi, L.; Niederwerder, M.; Bai, J. Genetic diversity and prevalence of porcine circovirus type 3 (PCV3) and type 2 (PCV2) in the Midwest of the USA during 2016-2018. Transbound Emerg Dis 2020, 67, 1284–1294. [Google Scholar] [CrossRef]
  9. Link, E.K.; Eddicks, M.; Nan, L.; Ritzmann, M.; Sutter, G.; Fux, R. Discriminating the eight genotypes of the porcine circovirus type 2 with TaqMan-based real-time PCR. Virol J 2021, 18, 70. [Google Scholar] [CrossRef]
  10. Gagnon, C.A.; Tremblay, D.; Tijssen, P.; Venne, M.H.; Houde, A.; Elahi, S.M. The emergence of porcine circovirus 2b genotype (PCV-2b) in swine in Canada. Can Vet J 2007, 48, 811–819. [Google Scholar]
  11. Gagnon, C.A.; del Castillo, J.R.; Music, N.; Fontaine, G.; Harel, J.; Tremblay, D. Development and use of a multiplex real-time quantitative polymerase chain reaction assay for detection and differentiation of Porcine circovirus-2 genotypes 2a and 2b in an epidemiological survey. J Vet Diagn Invest 2008, 20, 545–558. [Google Scholar] [CrossRef]
  12. Louten, J. Essential human virology. 2022: Academic Press.
  13. Zhuang, Q.; Wang, S.; Liu, S.; Hou, G.; Li, J.; Jiang, W.; Wang, K.; Peng, C.; Liu, D.; Guo, A. Diversity and distribution of type A influenza viruses: an updated panorama analysis based on protein sequences. Virology Journal 2019, 16, 1–38. [Google Scholar] [CrossRef]
  14. Chauhan, R.P.; Gordon, M.L. A systematic review analyzing the prevalence and circulation of influenza viruses in swine population worldwide. Pathogens 2020, 9, 355. [Google Scholar]
  15. Sanford, S.; Josephson, G.; Key, D. An epizootic of swine influenza in Ontario. The Canadian Veterinary Journal 1983, 24, 167. [Google Scholar]
  16. Bikour, M.H.; Frost, E.H.; Deslandes, S.; Talbot, B.; Elazhary, Y. Persistence of a 1930 swine influenza A (H1N1) virus in Quebec. Journal of general virology 1995, 76, 2539–2547. [Google Scholar] [CrossRef] [PubMed]
  17. Bikour, M.H.; Cornaglia, E.; Weber, J.M.; Elazhary, Y. Antigenic characterization of an H3N2 swine influenza virus isolated from pigs with proliferative and necrotizing pneumonia in Quebec. Canadian Journal of Veterinary Research 1994, 58, 287. [Google Scholar]
  18. Carman, S.; Stansfield, C.; Weber, J.; Bildfell, R.; Van Dreumel, T. H3N2 influenza A virus recovered from a neonatal pig in Ontario--1997. The Canadian Veterinary Journal 1999, 40, 889. [Google Scholar] [PubMed]
  19. Grgić, H.; Costa, M.; Friendship, R.M.; Carman, S.; Nagy, É.; Wideman, G.; Weese, S.; Poljak, Z. Molecular characterization of H3N2 influenza A viruses isolated from Ontario swine in 2011 and 2012. Virology journal, 2014. 11: p. 1-14.
  20. Grgić, H.; Costa, M.; Friendship, R.M.; Carman, S.; Nagy, É.; Poljak, Z. Genetic characterization of H1N1 and H1N2 influenza A viruses circulating in Ontario pigs in 2012. PLoS One 2015, 10, e0127840. [Google Scholar] [CrossRef] [PubMed]
  21. Nelson, M.I.; Culhane, M.R.; Trovão, N.S.; Patnayak, D.P.; Halpin, R.A.; Lin, X.; Shilts, M.H.; Das, S.R.; Detmer, S.E. The emergence and evolution of influenza A (H1α) viruses in swine in Canada and the United States. The Journal of general virology 2017, 98, 2663. [Google Scholar] [CrossRef]
  22. Nfon, C.K.; Berhane, Y.; Hisanaga, T.; Zhang, S.; Handel, K.; Kehler, H.; Labrecque, O.; Lewis, N.S.; Vincent, A.L.; Copps, J. Characterization of H1N1 swine influenza viruses circulating in Canadian pigs in 2009. Journal of virology 2011, 85, 8667–8679. [Google Scholar] [CrossRef]
  23. Ma, W. Swine influenza virus: Current status and challenge. Virus research, 2020. 288: p. 198118.
  24. La Gruta, N.L.; Kedzierska, K.; Stambas, J.; Doherty, P.C. A question of self-preservation: immunopathology in influenza virus infection. Immunology and cell biology 2007, 85, 85–92. [Google Scholar] [CrossRef]
  25. Sangpratum, N.; Kedkovid, R.; Woonwong, Y.; Sirisereewan, C.; Kesdangsakonwut, S.; Thanawongnuwech, R.; Arunorat, J. Dual infection of a Thai isolate HP-PRRSV and the pdmH1N1 2009 SIV in weanling pigs. The Thai Journal of Veterinary Medicine 2019, 49, 71–79. [Google Scholar] [CrossRef]
  26. Kitikoon, P.; Vincent, A.L.; Jones, K.R.; Nilubol, D.; Yu, S.; Janke, B.H.; Thacker, B.J.; Thacker, E.L. Vaccine efficacy and immune response to swine influenza virus challenge in pigs infected with porcine reproductive and respiratory syndrome virus at the time of SIV vaccination. Vet Microbiol, 2009. 139(3-4): p. 235-44.
  27. Van Reeth, K.; Nauwynck, H.; Pensaert, M. Clinical effects of experimental dual infections with porcine reproductive and respiratory syndrome virus followed by swine influenza virus in conventional and colostrum-deprived pigs. J Vet Med B Infect Dis Vet Public Health 2001, 48, 283–292. [Google Scholar] [CrossRef]
  28. Van Reeth, K.; Nauwynck, H.; Pensaert, M. Dual infections of feeder pigs with porcine reproductive and respiratory syndrome virus followed by porcine respiratory coronavirus or swine influenza virus: a clinical and virological study. Vet Microbiol, 1996. 48(3-4): p. 325-35.
  29. Czyzewska-Dors, E.; Pomorska-Mol, M.; Dors, A.; Pluta, A.; Podgorska, K.; Kwit, K.; Stasiak, E.; Lukomska, A. Proinflammatory Cytokine Changes in Bronchoalveolar Lavage Fluid Cells Isolated from Pigs Infected Solely with Porcine Reproductive and Respiratory Syndrome Virus or Co-infected with Swine Influenza Virus. J Vet Res 2019, 63, 489–495. [Google Scholar] [CrossRef]
  30. Niederwerder, M.C.; Jaing, C.J.; Thissen, J.B.; Cino-Ozuna, A.G.; McLoughlin, K.S.; Rowland, R.R. Microbiome associations in pigs with the best and worst clinical outcomes following co-infection with porcine reproductive and respiratory syndrome virus (PRRSV) and porcine circovirus type 2 (PCV2). Vet Microbiol, 2016. 188: p. 1-11.
  31. Fan, P.; Wei, Y.; Guo, L.; Wu, H.; Huang, L.; Liu, J.; Liu, C. Synergistic effects of sequential infection with highly pathogenic porcine reproductive and respiratory syndrome virus and porcine circovirus type 2. Virol J, 2013. 10: p. 265.
  32. Opriessnig, T.; Gauger, P.C.; Faaberg, K.S.; Shen, H.; Beach, N.M.; Meng, X.J.; Wang, C.; Halbur, P.G. Effect of porcine circovirus type 2a or 2b on infection kinetics and pathogenicity of two genetically divergent strains of porcine reproductive and respiratory syndrome virus in the conventional pig model. Vet Microbiol, 2012. 158(1-2): p. 69-81.
  33. Sinha, A.; Shen, H.G.; Schalk, S.; Beach, N.M.; Huang, Y.W.; Meng, X.J.; Halbur, P.G.; Opriessnig, T. Porcine reproductive and respiratory syndrome virus (PRRSV) influences infection dynamics of porcine circovirus type 2 (PCV2) subtypes PCV2a and PCV2b by prolonging PCV2 viremia and shedding. Vet Microbiol, 2011. 152(3-4): p. 235-46.
  34. Provost, C.; Hamonic, G.; Gagnon, C.A.; Meurens, F. Dual infections of CD163 expressing NPTr epithelial cells with influenza A virus and PRRSV. Veterinary microbiology, 2017. 207: p. 143-148.
  35. Dobrescu, I.; Levast, B.; Lai, K.; Delgado-Ortega, M.; Walker, S.; Banman, S.; Townsend, H.; Simon, G.; Zhou, Y.; Gerdts, V. In vitro and ex vivo analyses of co-infections with swine influenza and porcine reproductive and respiratory syndrome viruses. Veterinary microbiology, 2014. 169(1-2): p. 18-32.
  36. Richmond, O.; Cecere, T.E.; Erdogan, E.; Meng, X.J.; Pineyro, P.; Subramaniam, S.; Todd, S.M.; LeRoith, T. The PD-L1/CD86 ratio is increased in dendritic cells co-infected with porcine circovirus type 2 and porcine reproductive and respiratory syndrome virus, and the PD-L1/PD-1 axis is associated with anergy, apoptosis, and the induction of regulatory T-cells in porcine lymphocytes. Vet Microbiol, 2015. 180(3-4): p. 223-9.
  37. Cecere, T.E.; Meng, X.J.; Pelzer, K.; Todd, S.M.; Beach, N.M.; Ni, Y.Y.; Leroith, T. Co-infection of porcine dendritic cells with porcine circovirus type 2a (PCV2a) and genotype II porcine reproductive and respiratory syndrome virus (PRRSV) induces CD4(+)CD25(+)FoxP3(+) T cells in vitro. Vet Microbiol, 2012. 160(1-2): p. 233-9.
  38. Dorr, P.M.; Baker, R.B.; Almond, G.W.; Wayne, S.R.; Gebreyes, W.A. Epidemiologic assessment of porcine circovirus type 2 coinfection with other pathogens in swine. J Am Vet Med Assoc 2007, 230, 244–250. [Google Scholar] [CrossRef] [PubMed]
  39. Wei, H.; Lenz, S.D.; Van Alstine, W.G.; Stevenson, G.W.; Langohr, I.M.; Pogranichniy, R.M. Infection of cesarean-derived colostrum-deprived pigs with porcine circovirus type 2 and Swine influenza virus. Comp Med 2010, 60, 45–50. [Google Scholar] [PubMed]
  40. Misinzo, G.; Meerts, P.; Bublot, M.; Mast, J.; Weingartl, H.M.; Nauwynck, H.J. Binding and entry characteristics of porcine circovirus 2 in cells of the porcine monocytic line 3D4/31. J Gen Virol, 2005. 86(Pt 7): p. 2057-2068.
  41. Sanchez, R.E., Jr. Meerts, P.; Nauwynck, H.J.; Pensaert, M.B.; Change of porcine circovirus 2 target cells in pigs during development from fetal to early postnatal life. Vet Microbiol, 2003. 95(1-2): p. 15-25.
  42. Wu, N.-H.; Yang, W.; Beineke, A.; Dijkman, R.; Matrosovich, M.; Baumgärtner, W.; Thiel, V.; Valentin-Weigand, P.; Meng, F.; Herrler, G. The differentiated airway epithelium infected by influenza viruses maintains the barrier function despite a dramatic loss of ciliated cells. Scientific reports 2016, 6, 1–12. [Google Scholar] [CrossRef] [PubMed]
  43. Nicol, M.Q.; Dutia, B.M. The role of macrophages in influenza A virus infection. Future Virology 2014, 9, 847–862. [Google Scholar]
  44. Yu, W.C.; Chan, R.W.; Wang, J.; Travanty, E.A.; Nicholls, J.M.; Peiris, J.M.; Mason, R.J.; Chan, M.C. Viral replication and innate host responses in primary human alveolar epithelial cells and alveolar macrophages infected with influenza H5N1 and H1N1 viruses. Journal of virology 2011, 85, 6844–6855. [Google Scholar] [CrossRef] [PubMed]
  45. Günther, J. and Seyfert, H.-M. The first line of defence: insights into mechanisms and relevance of phagocytosis in epithelial cells. in Seminars in immunopathology. 2018. Springer .
  46. Vareille, M.; Kieninger, E.; Edwards, M.R.; Regamey, N. The airway epithelium: soldier in the fight against respiratory viruses. Clinical microbiology reviews 2011, 24, 210–229. [Google Scholar] [CrossRef] [PubMed]
  47. Cline, T.D.; Beck, D.; Bianchini, E. Influenza virus replication in macrophages: balancing protection and pathogenesis. The Journal of general virology 2017, 98, 2401. [Google Scholar] [CrossRef]
  48. Ferrari, M.; Scalvini, A.; Losio, M.N.; Corradi, A.; Soncini, M.; Bignotti, E.; Milanesi, E.; Ajmone-Marsan, P.; Barlati, S.; Bellotti, D.; et al. Establishment and characterization of two new pig cell lines for use in virological diagnostic laboratories. J Virol Methods 2003, 107, 205–212. [Google Scholar] [CrossRef]
  49. Kärber, G. Tabellen zur näherungsweisen Bestimmung von Titern. Naunyn Schmiedebergs Arch. Exp. Pathol. Pharmakol. 1931. 162: p. 480.
  50. Spearman, C. Review of The Method of ‘Right and Wrong Cases’(‘Constant Stimuli’) without Gauss’s Formula. Psychol Bull 1909, 6, 27–28. [Google Scholar]
  51. Racine, S.; Kheyar, A.; Gagnon, C.A.; Charbonneau, B.; Dea, S. Eucaryotic expression of the nucleocapsid protein gene of porcine circovirus type 2 and use of the protein in an indirect immunofluorescence assay for serological diagnosis of postweaning multisystemic wasting syndrome in pigs. Clin Diagn Lab Immunol 2004, 11, 736–741. [Google Scholar] [CrossRef]
  52. Schneider, C.A.; Rasband, W.S.; Eliceiri, K.W. NIH Image to ImageJ: 25 years of image analysis. Nature methods 2012, 9, 671–675. [Google Scholar] [CrossRef] [PubMed]
  53. Huang, D.; Sherman, B.; Lempicki, R. Systematic and integrative analysis of large gene lists using DAVID Bioinformatics Resources. Nature Protoc. 2009, 4, 44–57. [Google Scholar] [CrossRef] [PubMed]
  54. Huang, D.W.; Sherman, B.T.; Lempicki, R.A. Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic acids research 2009, 37, 1–13. [Google Scholar] [CrossRef] [PubMed]
  55. Jakab, G.J. Immune impairment of alveolar macrophage phagocytosis during influenza virus pneumonia. American Review of Respiratory Disease 1982, 126, 778–782. [Google Scholar] [PubMed]
  56. Du, Q.; Huang, Y.; Wang, T.; Zhang, X.; Chen, Y.; Cui, B.; Li, D.; Zhao, X.; Zhang, W.; Chang, L.; et al. Porcine circovirus type 2 activates PI3K/Akt and p38 MAPK pathways to promote interleukin-10 production in macrophages via Cap interaction of gC1qR. Oncotarget 2016, 7, 17492–17507. [Google Scholar] [CrossRef]
  57. Liu, S.; Li, Q.; Qiao, J.; Wang, J.; Cui, D.; Gu, K.; Zhou, S.; Li, H. Endothelial IL-8 induced by porcine circovirus type 2 affects dendritic cell maturation and antigen-presenting function. Virol J 2019, 16, 154. [Google Scholar] [CrossRef]
  58. Gu, C.; Gao, X.; Guo, D.; Wang, J.; Wu, Q.; Nepovimova, E.; Wu, W.; Kuca, K. Combined Effect of Deoxynivalenol (DON) and Porcine Circovirus Type 2 (Pcv2) on Inflammatory Cytokine mRNA Expression. Toxins 2021, 13. [Google Scholar] [CrossRef]
  59. Yang, S.; Liu, B.; Yin, S.; Shang, Y.; Zhang, X.; Khan, M.U.Z.; Liu, X.; Cai, J. Porcine Circovirus Type 2 Induces Single Immunoglobulin Interleukin-1 Related Receptor (SIGIRR) Downregulation to Promote Interleukin-1beta Upregulation in Porcine Alveolar Macrophage. Viruses 2019, 11. [Google Scholar] [CrossRef]
  60. Wu, X.; Wang, X.; Shi, T.; Luo, L.; Qiao, D.; Wang, Z.; Han, C.; Du, Q.; Tong, D.; Huang, Y. Porcine Circovirus Type 2 Rep Enhances IL-10 Production in Macrophages via Activation of p38-MAPK Pathway. Viruses 2019, 11. [Google Scholar] [CrossRef] [PubMed]
  61. Chang, H.W.; Jeng, C.R.; Lin, T.L.; Liu, J.J.; Chiou, M.T.; Tsai, Y.C.; Chia, M.Y.; Jan, T.R.; Pang, V.F. Immunopathological effects of porcine circovirus type 2 (PCV2) on swine alveolar macrophages by in vitro inoculation. Vet Immunol Immunopathol, 2006. 110(3-4): p. 207-19.
  62. Li, W.; Liu, S.; Wang, Y.; Deng, F.; Yan, W.; Yang, K.; Chen, H.; He, Q.; Charreyre, C.; Audoneet, J.C. Transcription analysis of the porcine alveolar macrophage response to porcine circovirus type 2. BMC Genomics, 2013. 14: p. 353.
  63. Van Reeth, K.; Van Gucht, S.; Pensaert, M. Correlations between lung proinflammatory cytokine levels, virus replication, and disease after swine influenza virus challenge of vaccination-immune pigs. Viral immunology 2002, 15, 583–594. [Google Scholar] [CrossRef] [PubMed]
  64. Wang, Z.; Chen, J.; Zhang, Q.-g.; Huang, K.; Ma, D.; Du, Q.; Tong, D.; Huang, Y. Porcine circovirus type 2 infection inhibits the activation of type I interferon signaling via capsid protein and host gC1qR. Veterinary Microbiology, 2022. 266: p. 109354.
  65. Hasslung, F.C.; Berg, M.; Allan, G.M.; Meehan, B.M.; McNeilly, F.; Fossum, C. Identification of a sequence from the genome of porcine circovirus type 2 with an inhibitory effect on IFN-alpha production by porcine PBMCs. J Gen Virol 2003, 84, 2937–2945. [Google Scholar] [CrossRef] [PubMed]
  66. Vincent, I.E.; Balmelli, C.; Meehan, B.; Allan, G.; Summerfield, A.; McCullough, K.C. Silencing of natural interferon producing cell activation by porcine circovirus type 2 DNA. Immunology 2007, 120, 47–56. [Google Scholar] [CrossRef] [PubMed]
  67. Choi, C.Y.; Choi, Y.C.; Park, I.B.; Lee, C.H.; Kang, S.J.; Chun, T. The ORF5 protein of porcine circovirus type 2 enhances viral replication by dampening type I interferon expression in porcine epithelial cells. Vet Microbiol, 2018, 226, 50-58.
  68. Li, J.; Lu, M.; Huang, B.; Lv, Y. Porcine circovirus type 2 inhibits interferon-β expression by targeting Karyopherin alpha-3 in PK-15 cells. Virology 2018, 520, 75–82. [Google Scholar] [CrossRef] [PubMed]
  69. Gao, F.; Xie, J.L.; Jia, C.W.; Ren, H.Y.; Zhou, S.H. Effects of porcine circovirus type 2 and pseudorabies vaccine co-inoculation on regulatory cytokine mRNA expression in pig peripheral blood mononuclear cells. Genet Mol Res 2014, 13, 1540–1547. [Google Scholar] [CrossRef]
  70. Chen, M.; Han, J.; Zhang, Y.; Duan, D.; Zhang, S. Porcine circovirus type 2 induces type I interferon production via MyD88–IKKα–IRFs signaling rather than NF-κB in porcine alveolar macrophages in vitro. Res. Vet. Sci. 2016, 104, 188–194. [Google Scholar] [CrossRef]
  71. Huang, B.; Zhang, L.; Lu, M.; Li, J.; Lv, Y. PCV2 infection activates the cGAS/STING signaling pathway to promote IFN-beta production and viral replication in PK-15 cells. Vet Microbiol 2018, 227, 34–40. [Google Scholar] [CrossRef]
  72. Wikstrom, F.H.; Meehan, B.M.; Berg, M.; Timmusk, S.; Elving, J.; Fuxler, L.; Magnusson, M.; Allan, G.M.; McNeilly, F.; Fossum, C. Structure-dependent modulation of alpha interferon production by porcine circovirus 2 oligodeoxyribonucleotide and CpG DNAs in porcine peripheral blood mononuclear cells. J. Virol. 2007, 81, 4919–4927. [Google Scholar] [CrossRef]
  73. Kekarainen, T.; Montoya, M.; Dominguez, J.; Mateu, E.; Segales, J. Porcine circovirus type 2 (PCV2) viral components immunomodulate recall antigen responses. Vet Immunol Immunopathol, 2008. 124(1-2), 41-9.
  74. Bakre, A.A.; Jones, L.P.; Murray, J.; Reneer, Z.B.; Meliopoulos, V.A.; Cherry, S.; Schultz-Cherry, S.; Tripp, R.A. Innate Antiviral Cytokine Response to Swine Influenza Virus by Swine Respiratory Epithelial Cells. J. Virol. 2021, 95, e00692-21. [Google Scholar] [CrossRef]
  75. Crosse, K.M.; Monson, E.A.; Beard, M.R.; Helbig, K.J. Interferon-stimulated genes as enhancers of antiviral innate immune signaling. J. Innate Immun. 2018, 10, 85–93. [Google Scholar] [CrossRef] [PubMed]
  76. De Veer, M.J.; Holko, M.; Frevel, M.; Walker, E.; Der, S.; Paranjape, J.M.; Silverman, R.H.; Williams, B.R. Functional classification of interferon-stimulated genes identified using microarrays. Journal of leukocyte biology 2001, 69, 912–920. [Google Scholar] [CrossRef] [PubMed]
  77. Schoggins, J.W.; Wilson, S.J.; Panis, M.; Murphy, M.Y.; Jones, C.T.; Bieniasz, P.; Rice, C.M. A diverse range of gene products are effectors of the type I interferon antiviral response. Nature 2011, 472, 481–485. [Google Scholar] [CrossRef] [PubMed]
  78. Schoggins, J.W.; Rice, C.M. Interferon-stimulated genes and their antiviral effector functions. Current opinion in virology 2011, 1, 519–525. [Google Scholar]
  79. Wu, W.; Wang, C.; Xia, C.; Liu, S.; Mei, Q. MicroRNA let-7 Suppresses Influenza A Virus Infection by Targeting RPS16 and Enhancing Type I Interferon Response. Front. Cell. Infect. Microbiol. 2022, 937. [Google Scholar] [CrossRef]
  80. Fong, C.H.-Y.; Lu, L.; Chen, L.-L.; Yeung, M.-L.; Zhang, A.J.; Zhao, H.; Yuen, K.-Y.; To, K.K.-W. Interferon-gamma inhibits influenza A virus cellular attachment by reducing sialic acid cluster size. Iscience 2022, 25, 104037. [Google Scholar] [CrossRef]
  81. Garcia-Sastre, A. Induction and evasion of type I interferon responses by influenza viruses. Virus research 2011, 162, 12–18. [Google Scholar] [CrossRef]
  82. Pauli, E.-K.; Schmolke, M.; Wolff, T.; Viemann, D.; Roth, J.; Bode, J.G.; Ludwig, S. Influenza A virus inhibits type I IFN signaling via NF-κB-dependent induction of SOCS-3 expression. PLoS pathogens 2008, 4, e1000196. [Google Scholar] [CrossRef]
  83. Czerkies, M.; Kochańczyk, M.; Korwek, Z.; Prus, W.; Lipniacki, T. RSV protects bystander cells against IAV infection by triggering secretion of type I and type III interferons. bioRxiv 2022, 11, 463877. [Google Scholar]
  84. Gaur, P.; Munjal, A.; Lal, S.K. Influenza virus and cell signaling pathways. Med. Sci. Monit. : Int. Med. J. Exp. Clin. Res. 2011, 17, RA148. [Google Scholar] [CrossRef]
  85. Kanehisa, M. and Goto, S.; KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Research 2000, 28, 27–30. [Google Scholar] [CrossRef]
  86. Coutermarsh-Ott, S.; Eden, K.; Allen, I.C. Beyond the inflammasome: regulatory NOD-like receptor modulation of the host immune response following virus exposure. J. Gen. Virol. 2016, 97 Pt 4, 825. [Google Scholar] [CrossRef] [PubMed]
  87. Pandey, K.P.; Zhou, Y. Influenza A Virus Infection Activates NLRP3 Inflammasome through Trans-Golgi Network Dispersion. Viruses 2022, 14, 88. [Google Scholar] [PubMed]
  88. Tate, M.D.; Mansell, A. An update on the NLRP3 inflammasome and influenza: the road to redemption or perdition? Current opinion in immunology 2018, 54, 80–85. [Google Scholar]
  89. Malik, G.; Zhou, Y. Innate immune sensing of influenza A virus. Viruses 2020, 12, 755. [Google Scholar] [PubMed]
  90. Tate, M.D.; Ong, J.D.; Dowling, J.K.; McAuley, J.L.; Robertson, A.B.; Latz, E.; Drummond, G.R.; Cooper, M.A.; Hertzog, P.J.; Mansell, A. Reassessing the role of the NLRP3 inflammasome during pathogenic influenza A virus infection via temporal inhibition. Scientific reports 2016, 6, 1–8. [Google Scholar] [CrossRef] [PubMed]
  91. Kesavardhana, S.; Samir, P.; Zheng, M.; Malireddi, R.S.; Karki, R.; Sharma, B.R.; Place, D.E.; Briard, B.; Vogel, P.; Kanneganti, T.-D. DDX3X coordinates host defense against influenza virus by activating the NLRP3 inflammasome and type I interferon response. Journal of Biological Chemistry, 2021. 296.
  92. Thulasi Raman, S.N.; Liu, G.; Pyo, H.M.; Cui, Y.C.; Xu, F.; Ayalew, L.E.; Tikoo, S.K.; Zhou, Y. DDX3 interacts with influenza A virus NS1 and NP proteins and exerts antiviral function through regulation of stress granule formation. Journal of Virology 2016, 90, 3661–3675. [Google Scholar] [CrossRef] [PubMed]
  93. Wurzer, W.J.; Planz, O.; Ehrhardt, C.; Giner, M.; Silberzahn, T.; Pleschka, S.; Ludwig, S. Caspase 3 activation is essential for efficient influenza virus propagation. The EMBO journal 2003, 22, 2717–2728. [Google Scholar] [CrossRef]
  94. Kumar, R.; Khandelwal, N.; Thachamvally, R.; Tripathi, B.N.; Barua, S.; Kashyap, S.K.; Maherchandani, S.; Kumar, N. Role of MAPK/MNK1 signaling in virus replication. Virus research, 2018. 253, 48-61.
  95. Reyes, A.; Corrales, N.; Gálvez, N.M.; Bueno, S.M.; Kalergis, A.M.; González, P.A. Contribution of hypoxia inducible factor-1 during viral infections. Virulence 2020, 11, 1482–1500. [Google Scholar] [CrossRef]
  96. Ren, L.; Zhang, W.; Zhang, J.; Zhang, J.; Zhang, H.; Zhu, Y.; Meng, X.; Yi, Z.; Wang, R. Influenza A virus (H1N1) infection induces glycolysis to facilitate viral replication. Virologica Sinica 2021, 36, 1532–1542. [Google Scholar] [CrossRef]
  97. Zhao, C.; Chen, J.; Cheng, L.; Xu, K.; Yang, Y.; Su, X. Deficiency of HIF-1α enhances influenza A virus replication by promoting autophagy in alveolar type II epithelial cells. Emerging Microbes Infections 2020, 9, 691–706. [Google Scholar] [CrossRef]
Figure 1. PCV2b and SwIV H1N1 co-localization in both NPTr and iPAM 3D4/21 co-infected cells. IFA was performed at 24 hours post-infection to detect PCV2b (red) and SwIV H1N1 (green) in co-infected cells. The cells were infected with 0.05 MOI of PCV2b and passaged 2-3 times. Then, PCV2b-infected cells were infected with 1 MOI of SwIV H1N1. Nuclear staining with DAPI is shown in blue. Co-localization pictures were realized with ImageJ. The red arrows point to single cells expressing both PCV2b and SwIV H1N1 antigens. Statistical analyses were done using a two-way ANOVA followed by Sidak's Multiple Comparison Test (GraphPad Prism software, version 7.00).
Figure 1. PCV2b and SwIV H1N1 co-localization in both NPTr and iPAM 3D4/21 co-infected cells. IFA was performed at 24 hours post-infection to detect PCV2b (red) and SwIV H1N1 (green) in co-infected cells. The cells were infected with 0.05 MOI of PCV2b and passaged 2-3 times. Then, PCV2b-infected cells were infected with 1 MOI of SwIV H1N1. Nuclear staining with DAPI is shown in blue. Co-localization pictures were realized with ImageJ. The red arrows point to single cells expressing both PCV2b and SwIV H1N1 antigens. Statistical analyses were done using a two-way ANOVA followed by Sidak's Multiple Comparison Test (GraphPad Prism software, version 7.00).
Preprints 71630 g001
Figure 2. PCV2b co-infection effects on SwIV replication in NPTr and iPAM 3D4/21 infected cells. The SwIV titer was determined in MDCK cells by the Spearman- Kärber method and expressed in tissue culture infectious dose 50% per mL (TCID50/mL). The PCV2b/SwIV H1N1 (A) and PCV2b/SwIV H3N2 (C) co-infection experiment in NPTr cells as well as the PCV2b/SwIV H1N1 (B) and PCV2b/SwIV H3N2 (D) co-infection experiment in iPAM 3D4/21 cells were repeated 3 times. Data are presented with standard deviation (SD) values. Statistical analyses were done using two-way repeated measures ANOVA (GraphPad Prism software, version 7.00). *p < 0.05, **p < 0.01, *** p < 0.001, **p < 0.01.
Figure 2. PCV2b co-infection effects on SwIV replication in NPTr and iPAM 3D4/21 infected cells. The SwIV titer was determined in MDCK cells by the Spearman- Kärber method and expressed in tissue culture infectious dose 50% per mL (TCID50/mL). The PCV2b/SwIV H1N1 (A) and PCV2b/SwIV H3N2 (C) co-infection experiment in NPTr cells as well as the PCV2b/SwIV H1N1 (B) and PCV2b/SwIV H3N2 (D) co-infection experiment in iPAM 3D4/21 cells were repeated 3 times. Data are presented with standard deviation (SD) values. Statistical analyses were done using two-way repeated measures ANOVA (GraphPad Prism software, version 7.00). *p < 0.05, **p < 0.01, *** p < 0.001, **p < 0.01.
Preprints 71630 g002
Figure 3. SwIV co-infection effects on PCV2b replication in NPTr and iPAM 3D4/21 infected cells. PCV2b quantification in co-infected NPTr (A) and iPAM 3D4/21 (B) cells was performed by qPCR assay and expressed as PCV2b genome copies per mL of sample. The experiments were repeated 3 times. Data are presented with standard deviation (SD) values. Statistical analyses were done using two-way repeated measures ANOVA (GraphPad Prism software, version 7.0.0).
Figure 3. SwIV co-infection effects on PCV2b replication in NPTr and iPAM 3D4/21 infected cells. PCV2b quantification in co-infected NPTr (A) and iPAM 3D4/21 (B) cells was performed by qPCR assay and expressed as PCV2b genome copies per mL of sample. The experiments were repeated 3 times. Data are presented with standard deviation (SD) values. Statistical analyses were done using two-way repeated measures ANOVA (GraphPad Prism software, version 7.0.0).
Preprints 71630 g003
Figure 4. Cell viability was determined in single-infected and co-infected PCV2b/SwIV cells at 24 hpi (A and B), 48 hpi (C and D) and 72 hpi (E and F) in NPTr (A, C and E) and iPAM 3D4/21 cells (B, D and F). The experiments were repeated three times. The data represents percentage of cell viability in infected cells with respect to mock-infected cells and are presented with standard deviation (SD) values. Statistical analyses were done using one way ANOVA followed by Tukey's Multiple Comparison Test (GraphPad Prism software, version 7.0.0). Different superscripts indicate a statistically significant difference (p < 0.05) between groups.
Figure 4. Cell viability was determined in single-infected and co-infected PCV2b/SwIV cells at 24 hpi (A and B), 48 hpi (C and D) and 72 hpi (E and F) in NPTr (A, C and E) and iPAM 3D4/21 cells (B, D and F). The experiments were repeated three times. The data represents percentage of cell viability in infected cells with respect to mock-infected cells and are presented with standard deviation (SD) values. Statistical analyses were done using one way ANOVA followed by Tukey's Multiple Comparison Test (GraphPad Prism software, version 7.0.0). Different superscripts indicate a statistically significant difference (p < 0.05) between groups.
Preprints 71630 g004
Figure 5. PCV2b/SwIV co-infection effects on the modulation of cytokines mRNA expression. The mRNA expressions of IL-6, IL-8 and IL-10 in NPTr cells (A) and iPAM 3D4/21 (B) as well as of IFN-α, IFN-β and IFN-γ in NPTr cells (C) and iPAM 3D4/21 (D) were determined by RT-qPCR assays. The 2−ΔΔCt method was used to calculate the fold change of cytokine mRNA expression in infected cells with respect to mock-infected cells at 24 hours post-infection. The experiments were repeated at least three times. All data are presented with standard deviation (SD) values. Statistical analyses were done using an ordinary two-way ANOVA followed by Tukey's Multiple Comparison Test (GraphPad Prism software, version 7.0.0). Different superscripts indicate a statistically significant difference (p < 0.01) between groups within the same tested cytokine mRNA.
Figure 5. PCV2b/SwIV co-infection effects on the modulation of cytokines mRNA expression. The mRNA expressions of IL-6, IL-8 and IL-10 in NPTr cells (A) and iPAM 3D4/21 (B) as well as of IFN-α, IFN-β and IFN-γ in NPTr cells (C) and iPAM 3D4/21 (D) were determined by RT-qPCR assays. The 2−ΔΔCt method was used to calculate the fold change of cytokine mRNA expression in infected cells with respect to mock-infected cells at 24 hours post-infection. The experiments were repeated at least three times. All data are presented with standard deviation (SD) values. Statistical analyses were done using an ordinary two-way ANOVA followed by Tukey's Multiple Comparison Test (GraphPad Prism software, version 7.0.0). Different superscripts indicate a statistically significant difference (p < 0.01) between groups within the same tested cytokine mRNA.
Preprints 71630 g005
Figure 6. Venn diagram of differentially expressed genes (DEGs) in PCV2b/SwIV co-infected cells. Mock-infected cells were used as control to identify DEGs in NPTr (A) and iPAM 3D4/21 (B) single infected and co-infected cells. The numbers in overlapping areas represent the number of DEGs shared among the different virus infection experimental groups. A False Discovery Rate (FDR) <0.05 and a fold change cut-off of 1.5 was used to identify the DEGs.
Figure 6. Venn diagram of differentially expressed genes (DEGs) in PCV2b/SwIV co-infected cells. Mock-infected cells were used as control to identify DEGs in NPTr (A) and iPAM 3D4/21 (B) single infected and co-infected cells. The numbers in overlapping areas represent the number of DEGs shared among the different virus infection experimental groups. A False Discovery Rate (FDR) <0.05 and a fold change cut-off of 1.5 was used to identify the DEGs.
Preprints 71630 g006
Figure 7. KEGG pathway enrichment analysis of the DEGs identified in single infected and co-infected cells. Enriched KEGG pathways after PCV2b, SwIV H1N1 and PCV2b/SwIV H1N1 infection in NPTr cells (A, C and E, respectively) and in iPAM 3D4/21 cells (B, D and F, respectively) are illustrated. The FDR value < 0.05 was considered for statistically significant enriched pathways.
Figure 7. KEGG pathway enrichment analysis of the DEGs identified in single infected and co-infected cells. Enriched KEGG pathways after PCV2b, SwIV H1N1 and PCV2b/SwIV H1N1 infection in NPTr cells (A, C and E, respectively) and in iPAM 3D4/21 cells (B, D and F, respectively) are illustrated. The FDR value < 0.05 was considered for statistically significant enriched pathways.
Preprints 71630 g007
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

Subscribe

Disclaimer

Terms of Use

Privacy Policy

Privacy Settings

© 2025 MDPI (Basel, Switzerland) unless otherwise stated