Distinct functional enrichment of transcriptional signatures in pigs with high and low IFN-gamma responses after vaccination with a porcine reproductive and respiratory syndrome virus (PRRSV)
Veterinary Research volume 47, Article number: 104 (2016)
Little is known about the host factor in the response to PRRSV vaccination. For this purpose, piglets were immunized with a commercial PRRSV-live vaccine and classified as high responders (HR) or low responders (LR) as regards to the frequencies of virus-specific IFN-γ-secreting cells. Six weeks post vaccination, PBMCs isolated from three individuals with the most extreme responses in each HR and LR groups and 3 unvaccinated controls, were either stimulated with phytohaemagglutinin, challenged with the vaccine or mock treated for 24 h, prior conducting transcriptional studies, gene ontology and pathway analyses. The LR group had very low neutralizing antibody levels and showed a higher number of down-regulated transcripts compared with the HR group (FDR < 0.2, P < 0.001). Down-regulated genes encoded chemoattractants, proinflammatory cytokines and the interferon-inducible GBP family, and showed enrichment in wounding (FDR < 3.6E-13), inflammation (FDR < 8E-12), defence (FDR < 8.7E-09) and immunity (FDR < 7.6E-08), suggesting immune response impairment. In the HR group, down-regulated genes were involved in protein transport (FDR < 4.77E-03), locomotory behavior (FDR < 5.47E-3), regulation of protein localization (FDR < 1.02E-02), and regulation of TNF superfamily member 15 and miR181. In contrast, the HR group presented up-regulated transcripts associated with wounding (FDR < 4.95). Moreover, IFN-γ was predicted to be an inhibited upstream regulator since IFN-γ pathways were associated with higher number of down-regulated genes in the LR (n = 40) than the HR (n = 10). Divergent responses to PRRSV-vaccination may be the result of the genetic background of the host.
Porcine reproductive and respiratory syndrome (PRRS) is arguably the most endemic infectious disease challenge for the pig industry and causes very significant economic losses worldwide [1, 2]. The causative agent is an enveloped, single-stranded 15-kb positive-sense RNA virus belonging to the Arteriviridae family in the order Nidovirales  known as PRRS virus (PRRSV). The extensive genetic diversity characteristic of PRRSV encodes considerable antigenic diversity which presents a highly variable challenge to the host’s immune system [4–12]. Moreover PRRSV is known to modify the host’s innate immune response, particularly by inhibiting type-I interferons (IFNs), affecting cytokine production or by down regulating toll-like receptor (TLRs) expression [13–19].
The most widely used strategy for reducing the incidence and limiting the impact of PRRSV-infection damages is vaccination. Several commercial vaccines (including live attenuated and inactivated) are available on the market but their efficacies are considered to be at best partial as they cannot provide full protection against the wide diversity of PRRSV strains circulating in the field [20, 21]. PRRSV vaccines are routinely applied to breeding herds to prevent reproductive problems but the beneficial effect of vaccination to piglets is still controversial. After immunization of a group of pigs with a given vaccine, a variety of immune response are often generated among pigs and only some of them are being protected against infection. Thus, the genetics of the host (pigs) and of the pathogen (PRRSV) appear to play an important role on the efficacy of vaccination in the control of PRRS.
The interaction of PRRSV with the immune system is believed to be of critical importance for defining immunological and clinical outcomes of the infection and, along with other factors, is associated with an inefficient development of the adaptive immunity. The prototypical adaptive response to PRRSV is characterized by a weak and delayed production of neutralizing antibodies and by relatively low levels of IFN-γ producing T-cells [22, 23]. Moreover, there is increasing evidence that the genetic variation of both host and virus also have impacts on the response to infection [24–29].
Therefore, in order to understand immune responses to PRRSV infection and/or vaccination it is necessary to examine both sources of variation: host (pig) and pathogen (PRRSV). Analysis of the host response using transcriptional profiling of pigs with different immune status against PRRSV may help to understand why vaccines offer only partial protection. The objective of this study was to identify transcriptional signals and molecular pathways associated with responses to a PRRSV vaccine . Systems biology approaches have been used successfully to identify early gene signatures that predict immune responses in humans vaccinated with the live attenuated vaccine YF-17D against yellow fever virus and inactivated seasonal influenza vaccine against influenza virus with accuracy [31, 32]. In order to explore this approach for PRRSV vaccines pigs were vaccinated with a modified live (ML) PRRSV vaccine and later classified into two groups, namely “high responders” (HR) and “low responders” (LR) on the basis of the frequencies of virus-specific interferon-γ secreting cells (IFN-γ SC). Subsequently, the humoral response against the virus and the transcriptional profiles of peripheral blood mononuclear cells (PBMC) of both groups were examined in order to identify genes, molecular functions and biological pathways associated with the vaccine response.
Materials and methods
Animals and study design
Twenty-two healthy four-week-old Landrace × Large White PRRSV-negative pigs (as determined by serology and RT-PCR) were vaccinated intramuscularly with a genotype-1 commercial ML vaccine (Pyrsvac-183) according to the manufacturer’s instructions. Pigs were followed up for 42 days after vaccination. In parallel, 10 control pigs were subjected to a mock vaccination with sterile saline solution (SS). Pigs were bled for immunological determinations at day 0 (pre-vaccination) and then once a week post-vaccination (PV) for 6 weeks: days 7, 14, 21, 28, 35 and 42 PV (Figure 1).
Peripheral blood mononuclear cell isolation
Heparinized blood samples were used for separation of PBMC by density-gradient centrifugation with Histopaque 1.077 (Sigma). PBMC were subsequently cultured using RPMI medium supplemented with 10% foetal calf serum (FCS) (Invitrogen, Madrid, Spain), 1 mM non-essential amino acids (Invitrogen), 1 mM sodium pyruvate (Invitrogen), 5 Mm 2-mercaptoethanol (Sigma), 50 000 IU penicillin (Invitrogen), 50 mg streptomycin (Invitrogen) and 50 mg gentamicin (Sigma) (complete RPMI). Trypan blue was used to assess cell viability.
ELISPOT assays were performed as described previously [12, 33] using commercial mAbs (P2G10 and biotin P2C11; BD Biosciences Pharmingen). Peripheral blood mononuclear cells (PBMC) were stimulated for 24 h (5 × 105 PBMC/well, in triplicate) with the ML vaccine strain at a multiplicity of infection (m.o.i.) of 0.1 as a recall antigen, phytohaemagglutinin (PHA, 10 μg/mL) or culture medium. To calculate the frequencies of virus specific IFN-γ-SC, average counts of spots in unstimulated wells were subtracted from average counts obtained in antigen stimulated wells. Results were expressed as the number of PRRSV-specific IFN-γ-SC per 5 × 105 PBMC.
Three pigs with the highest virus-specific IFN-γ-SC average responses throughout the timecourse of the study, designated as HR (High responders), and 3 pigs with the lowest average levels, or LR (Low responders), were selected for further analysis (Figure 1). Average HR and LR frequencies at 21 days PV were 26 ± 3/5 × 105 cells versus 12 ± 5.6/5 × 105 cells and at 42 days, 54.3 ± 8.1/5 × 105 cells versus 22.3 ± 4.5/5 × 105 cells (Figure 2). Three healthy unvaccinated pigs (SS) were randomly taken as a control group.
A commercially available ELISA (Idexx PRRS X3 Ab Test) was used for measuring non neutralizing anti-PRRSV antibodies in blood samples taken at the different time-points of the study. Results were reported as a ratio of optical densities between the results of a given sample and the positive control (S/P ratio) included in the kit (cut-off: S/P ≥ 0.4). All samples were tested the same day and in the same batch of plates.
In parallel, homologous (against the vaccine virus) neutralizing antibodies (NA) titers were determined in a viral neutralization test (VNT). The VNT assays were done in MARC-145 cells using log2 serial dilutions (from 2−1 to 2−8) as described before running samples in duplicate [34, 35].
In vitro treatment of peripheral blood mononuclear cells
At 6 weeks PV (Figure 1), PBMC isolated from the selected pigs of each group (HR and LR) and from unvaccinated pigs were seeded at 5 × 105 cells/well in 96 well-plates and were cultured in triplicates with complete RPMI alone (negative controls), complete RPMI with phytohaemagglutinin (PHA, 10 μg/mL, positive controls) or the PRRSV vaccine strain (m.o.i. 0.1, specific stimulus) for 24 h. After that period, cells were harvested, centrifuged and cell pellets were resuspended in 1 mL of Trizol (Life Technologies) and frozen at −80 °C until processing for RNA extraction, microarray and downstream analysis.
Total RNA from PBMC of selected group pigs at 6 weeks PV was extracted using Trizol (Invitrogen, Paisley UK) according to the manufacturer’s instructions and was later purified using the Qiagen RNeasy minikit (Qiagen, Crawley, UK). RNA was eluted from the spin column in 30 µL of RNase-free water and stored as aliquots at −80 °C. The quantity and quality of RNA were assessed using a Nanodrop spectrophotometer (NanoDrop Technologies Inc, Wilmington, DE, USA) and Agilent 2100 bioanalyser (Agilent Technologies, Palo Alto, CA, USA) respectively.
Microarray platform and data analysis
The extracted RNA was then analysed using the Affymetrix Snowball GeneChip® . This GeneChip comprises 23 937 probe sets that interrogate approximately 23 256 transcripts from 20 201 Sus scrofa genes. Sense-strand cDNA was generated from total RNA (500 ng) subjected to two rounds of amplification (Ambion® WT Expression Kit). The resulting cDNA was used for biotin labelling and fragmentation according to the Affymetrix GeneChip® WT Terminal Labelling and Hybridization protocol (Affymetrix UK, High Wycombe). Biotin-labelled fragments of cDNA (5.5 μg) were hybridized to Affymetrix Snowball arrays using the Affymetrix HybWashStain kit and following the manufacturer’s recommendations. After hybridization, the arrays were washed and stained using the Affymetrix Fluidics Station 450 and then scanned in an Affymetrix 7G scanner. Image generation and the resulting CEL files for analysis were produced in AGCC—Affymetrix GeneChip Command Console Software. Initial QCs were performed in Expression Console. All microarray data used in the analyses herein are available from the Array Express repository . The Affymetrix.CEL files were imported into the Partek Genomics Suite software package version 6.13.0213 (Partek, St. Louis, USA) for data analysis. Transcriptional responses were normalised to those from unvaccinated pigs prior to running an ANOVA analysis of the data (Additional file 1). Up-regulated and down-regulated differentially expressed transcripts in HR and LR were selected for further consideration if the false discovery rate (FDR) was ≤0.2.
For network analysis, the normalised array data were uploaded to the software Biolayout Express3D as described previously [36, 38]. Expression variations across treated groups were generated using the gplots package in R.
Quantitative real time RT-PCR validation (qRT-PCR)
The differential expression of several selected genes, as identified from the microarray data, was verified at various time points using qRT-PCR. Reverse transcription was performed as described previously [26, 27]. Briefly, one microgram of total RNA was reverse transcribed using a TaqMan kit (Applied Biosystems, Foster City, CA, USA). For qRT-PCR, Platinum SYBR Green PCR SuperMix UDG was used, as described above. The qRT-PCR was performed with a Stratagene MX3000P (Stratagene, La Jolla, CA, USA). In Table 1 is showed the information about primer sequences of the selected genes. Samples were tested in triplicate, GAPDH served as the housekeeping gene and results were calculated as described previously [26, 27].
Gene ontology and pathway analysis
Gene ontology and pathway analysis were carried out using DAVID bioinformatics resources  and Ingenuity pathway analysis , respectively. In particular for the identification of upstream regulators using IPA, expression datasets for the HR and the LR pigs were analysed independently with Fisher’s exact test as described by the software. The activation z-score was used to infer the state of activation of upstream regulators based on a comparison with a model that assigns random regulations. The P value overlap, which indicates possible upstream regulators, represents the significance of the overlap between the dataset genes identified and the known targets of transcriptional regulators. Differences were considered significant with P-value <0.05 and activation Z-score ≥2.
Comparison of frequencies of IFN-SC, S/P ratios or titres of NA between groups was done using StatsDirect3 (v3.0.97) and the Kruskal–Wallis test. The coefficient of correlation (R2) between microarray and qRT-PCR techniques was produced by plotting the relative level of gene expression as assessed by microarray against the transcript levels generated by qRT-PCR.
Characterization of the humoral response in HR and LR pigs
When the S/P ratios, corresponding to non-neutralizing antibodies (NA), were examined, values were always higher (p < 0.05) in the LR vaccinated group compared to the HR vaccinated group from day 14 PV onwards (Figure 3A). In contrast, HR had an earlier and stronger NA response to PRRSV vaccine at 28, 35 and 42 days PV in comparison to the LR group (Figure 3B). Regarding the viremia, CT means for HR were 28.3 ± 1.8 (day 14); 27.8 ± 1.7 (day 21); 29.9 ± 3.9 (day 28); and 34.0 ± 5.0 (day 35). For LR, CT means were 30.9 ± 3.4 (day 14); 28.4 ± 5.7 (day 21); 32.9 ± 1.8 (day 28); and 34.7 ± 1.6 (day 35). All animals were negative for PRRSV vaccine virus in blood at day 42, with the exception of one animal of the LR group that showed a CT = 35 that corresponds to 1.2–1.3 log10 genomic equivalences per mL.
Host transcript regulation during PRRSV vaccine challenge
From the total number of evaluated genes, 29 and 127 transcripts were identified as differentially expressed in the mock stimulated PBMC cultures of HR and LR vaccinated pigs, respectively (ANOVA FDR <0.2, P < 0.001). The proportion of genes which showed reduced levels of expression either in the presence of PHA treatment or in mock-stimulated cultures was larger in the LR vaccinated group compared with the HR (Table 2). In contrast, PBMCs from HR and LR vaccinated pigs and treated in vitro with the vaccine did not harbor statistically different regulated transcripts when compared to the controls (Table 2). To investigate further this observation, we examined differential expression using network analysis with Biolayout Express3D. Nine probesets mapping to 7 unique gene symbols showed a differential expression of >2 Log fold change between the PBMCs of LR and HR vaccinated pigs and unvaccinated control groups after vaccine stimulation (Figure 4, red squared section). Most of them showed increased expression in LR and HR groups compared to controls. Among these genes we found three chemokines, CXCL9, CXCL10 and CXCL11, which activate monocytes and T-lymphocytes, two Cytochromes, CYP3A46 and CYP3A29 of which the latter has been reported to be regulated by interferon , and OAS1 which is an interferon-induced antiviral enzyme .
Validation of differentially regulated transcripts using quantitative real time PCR (qRT-PCR)
Nine genes for which changes in transcript levels had been observed were selected for validation by qRT-PCR analysis. The selection criteria were: (i), differential regulation in the both groups (fold change ≥2.0, FDR ≤0.2) and ii), involvement in immune functions such as inflammation, antiviral responses and cell growth and recruitment molecules. Thus we selected the following genes: GBP2, IL1A, IL8, CCL4, SAA1, IGF1, CCL2, IL1B and LTBP1. The results showed positive correlation coefficients between the microarray and qRT-PCR techniques in the HR (R2 = 0.76) and LR (R2 = 0.86) groups (Additional file 2).
Gene ontology analyses
Transcripts differentially regulated in the HR and LR vaccinated groups (Additional file 1) were analysed with respect to their biological processes, molecular functions and cellular (Additional file 3). Interestingly, the LR vaccinated group presented a higher number of down-regulated genes (6–25 genes) per biological processes category than the HR vaccinated (3–6 genes) but in a more restricted range of categories (9 categories vs 18 categories in the HR) (Figure 5). In the LR group mostly down-regulated transcripts showed significant fold functional enrichment (FFE) in 9 biological process involved in inflammation and immune responses as follow: response to wounding (8.6 FFE, FDR < 3.6E-13), inflammatory response (11.2 FFE, FDR < 8E-12), defence response (6.5 FFE, FDR < 8.7E-09) and immune response (5.4 FFE, FDR < 7.6E-08) (Additional file 3). Our results suggest that two broad molecular functions involved in chemokine receptor binding/activity and cytokine activity were down-regulated in LR and included genes encoding a wide range of chemoattractants (CCL2, CCL17, CCL22 and CCL23) and the proinflammatory cytokines IL1A and ILIB and receptor antagonist IL1RN. Among the most down-regulated genes involved in immune response included several members of interferon-inducible guanylate binding protein (GBP) family GBP2 (13.1 fold), GBP4 (4.6 fold) and GBP6 (9.9 fold). Finally, whereas no significant function dominated the category of up-regulated genes for the LR group, the greatest fold changes were observed for insulin-like growth factor 1 (IGF1, 13.0 fold).
In the HR vaccinated group, although the number of genes that were down-regulated per category was lower than the LR, the number of biological processes involved was twice as great as in the LR vaccinated pigs (Figure 5). The biological functional fold enrichments were associated with the regulation of protein transport, locomotory behaviour and regulation of protein localization which were enriched 42.4 FFE (FDR < 4.77E-03), 21.2 FFE (FDR < 5.47E-3) and 35 FFE (FDR < 1.02E-02), respectively. It is noteworthy that down-regulated genes in this group harboured the same trend than in the LR vaccinated group. For example SAA1, CCL3L1 and CCL4 were down-regulated group (Additional file 3). Also transcripts encoding for the tumor necrosis factor (TNF) superfamily member 15 (TNFSF15) (−2.2 fold, P < 7.9E-05) and miR181 (−2.4 fold, P < 2.1E-04) appeared specifically down-regulated in a statistical significant manner. In contrast, up-regulated transcripts in the HR vaccinated group were barely dominated by a biological function that was the response to wounding (50.2 FFE, FDR < 4.95) involving purinergic P2Y G-protein coupled receptor (P2RY12), vanin1 (VNN1) and IGF1 (Additional file 3).
Upstream regulators using Ingenuity pathway analysis (IPA)
To identify the cascade of upstream transcriptional regulators that could explain the observed gene expression data for the HR and LR vaccinated groups we performed in silico analysis using IPA. Table 3 shows top 10 upstream regulators, each predicted to act upon different sets of target molecules, likely mediating different biological effects. HR and LR datasets shared a similar set of inhibited upstream regulators (P < 1E-06) with the exception of IL1A and STAT3 which were specific to the LR vaccinated group (Table 3; Additional file 4). It is noteworthy that IFNG was identified as an inhibited upstream regulator with a greater activation z-score in the LR vaccinated group (−4.4, P value < 3.36E-20) than in the HR group (−2.4, P value < 1.26E-06). The results of the in silico analysis of the relationship between IFNG and target genes show that 10 and 40 predicted relationships were identified with IFNG and target genes in the HR and LR vaccinated group, respectively (Figure 6). Thus, IFNG was predicted to be an inhibited upstream regulator to explain the trend of regulation in HR and LR vaccinated groups.
It is increasingly recognized that the control of PRRS will require a global approach that will consider relevant factors of the pathogen, the host and the environment. Although valuable tools, current PRRSV vaccines are not fully efficacious and some vaccinated animals may still become infected . It is thought that this partial lack of efficacy may be attributable to the genetic and antigenic diversity of the virus. However, several studies have highlighted that there is high inter-individual variability in the development of the immune response and protection after vaccination with a PRRSV vaccine [12, 44, 45] which may account for some lack of efficacy. Therefore one of the challenges of PRRS control is the development of tools to predict the efficacy of vaccines with greater accuracy.
To date, measurement of humoral immunity (neutralizing antibodies) and cell-mediated immunity (IFN-γ responses levels) have been used for monitoring PRRSV infection [20, 21, 46]. We hypothesized that the generation and the maintenance of protective immunity to PRRSV vaccine are a complex process that entails the regulation of transcriptional pathways leading to the production of neutralizing antibody and the induction of cell-mediated immunity. Thus, in the present study we have focused primarily in the PRRSV-specific IFN-γ responses after vaccination as a mean of discriminating between animals. The three animals with the most extreme results in the ELISPOT were selected. Our results show that in a group of cross-bred pigs of the same origin, responses to PRRSV vaccination may range from low IFN-γ responses and low or nil NA, as in the Low Responders (LR) group, to higher IFN-γ responses and moderate titers of NA, in the High Responders (HR) group. These differences indicate that inter-individual variability in markers of PRRSV vaccine-induced immunity, including neutralizing antibody levels and INF-γ responses, may be regulated, at least in part by host genetic factors . This observation is indeed reminiscent with current trends suggesting that host genetic variation in the outcomes of PRRSV infections occurs in growing pigs as well as during reproduction . Specifically, a QTL on chromosome SSC4 was found to be associated with viraemia levels and the growth rate in experimentally infected pigs . Further, the SNPs defining this QTL have been shown to be associated with resistance and performance in several independent populations [48–50]. Collectively, these studies suggest that immune response to PRRSV is likely the results of multigenic influences and possibly not of a single dominant gene/allele and advocate the need for further genetic investigation of the host response to PRRSV vaccine.
To further delineate which genes and pathways are associated with HR and LR vaccinated pigs we have conducted a genome–wide transcriptional study of PBMCs 6 weeks post vaccination. The selection of the time-points was based on the likelihood of having the minimal influence with residual vaccine virus in blood  and to choose a moment in the development of the immune response when cell-mediated immunity was more or less in a steady state [11, 22].
The number of examined animals per group can represent a critical factor to reach statistical power for comparisons for such studies. Indeed, the high inter-individual variability in gene expression of control animals noticed using Biolayout Express3D software could explain why statistical significance was not achieved among vaccinated and unvaccinated groups after vaccine treatment (Table 2; Figure 4). However, the results of gene expression after vaccine stimulation showed an enrichment of immune-related genes in the LR and HR vaccinated groups as compared to the unvaccinated controls. Most of these genes, CXCL9, CXCL10 and CXCL11 were related with monocyte and T-lymphocyte activation and with interferon-modulation antiviral responses (CYP3A46, CYP3A29 and OAS1) [41, 42]. For this reason we considered that these changes in gene expression were potentially biologically relevant despite not being selected using the Anova test with Partek software. On the other hand and despite the limited number of pigs tested, statistical significance was reached for the other group comparisons (mock and PHA stimulated) suggesting that future transcriptional study of vaccine efficacy could benefit from larger number of animals per experimental group.
Gene regulation in the LR vaccinated group suggested that in those animals there was a functional impairment in at least some components of the inflammatory response such as CCL2 or IL1A but also that transcript levels or transcription of genes encoding some chemokines that act as strong chemotactic compounds for T cells, NK cell, dendritic cells and monocytes/macrophages are affected (CCL17, CCL22, CCL23). Interestingly down-regulated transcripts encoding guanylate-binding proteins GBP2, 4 and 6 were identified. Members of the GBP gene family are associated with antiviral activity in humans and as discussed previously by Boddicker et al.  they have been highlighted as potential candidate genes for the QTL identified on SSC4 which explains 15% of the variance in genomic estimated breeding values for viral load and 11% for weight gain 42 days post infection . Recently, an intronic SNP in GBP5 has been identify as a strong candidate causal mutation for the SSC4 QTL that controls variation in host response to PRRSV . As regards the HR pigs, either the lower proportion of inhibited up stream regulators of the IFN-γ cascade and the up-regulation of transcripts encoding for molecules involved in anti-apoptosis mechanisms (P2RY12, VNN1 and IGF1), suggests increased chemotactic activity, viral-induced-apoptosis resistance and a better induction of cellular immune response.
Curiously, HR vaccinated animals presented a down-regulation of a member of the tumor necrosis factor superfamily, TNFSF15, and miR181 transcripts. TNFSF15 is produced by inflammatory cytokine-stimulated endothelial cells and by TLR ligands-stimulated antigen presenting cells, and it is involved in apoptosis, cell proliferation and polarization to Th1 and Th17 responses . Recently, increased levels in sera and genetic polymorphisms of TNFSF15 have been observed in inflammatory diseases that are caused by altered immunological reactions such as Crohn’s disease or ulcerative colitis . As regards to miR181, over-expression of this molecule has been recently identified in minimally PRRSV-permissive cells or tissues, and this molecule is able to inhibit PRRSV replication by both mechanisms: binding to a conserved region in the downstream of open reading frame 4  and downregulating the PRRSV receptor CD163 in blood monocytes and porcine alveolar macrophages . The fact that TNFSF15 and miR181 were found down-regulated in HR vaccinated pigs is noteworthy and warrant further studies. It is possible that these factors may be associated with either protective mechanism induced against harmful effects of an exacerbated inflammatory response and/or a higher cell permissibility to the vaccine virus and thus, better priming and immune stimulation.
It is likely that at 6 weeks PV the transcriptional response we have measured may have been mostly associated with the resolution phase rather than the primary response. However, the finding that the transcriptional signature in LR pigs was associated with a down-regulation of the IFN-γ response was interesting for several reasons: this specific transcriptional signature was consistent with the low IFN-γ level found in the LR or the HR pigs maybe resolving some of the vaccine response. On the other hand, the results of the experiment are compatible with either a faster or a stronger response in the HR group. In each case, implications would be a little different. If the differences were related just to the intensity, the outcome of vaccination in LR could maybe be ameliorated by use of adjuvants; if the differences relied in the velocity of the immune response this could be more much difficult to solve. However, in practical terms, in situations where the risk of entering in contact with the virus are high, the point is to achieve the fastest and strongest possible response. As shown here, at 42 days PV a group of animals was far from having that type of immunity and this seemed to be related to some elements in the genetics of the pig.
In summary, this works showed that variable responses to PRRSV vaccination may be the result of the genetic characteristics of the host. If this was to be proven in future challenge studies, genetic selection of pigs could be used to increase resistance to the infection through vaccination. Taken together this work advocate the need for a greater understanding of the host genomics/genetics of PRRSV vaccine response.
Neumann EJ, Kliebenstein JB, Johnson CD, Mabry JW, Bush EJ, Seitzinger AH, Green AL, Zimmerman JJ (2005) Assessment of the economic impact of porcine reproductive and respiratory syndrome on swine production in the United States. J Am Vet Med Assoc 227:385–392
Nieuwenhuis N, Duinhof TF, van Nes A (2012) Economic analysis of outbreaks of porcine reproductive and respiratory syndrome virus in nine sow herds. Vet Rec 170:225
Cavanagh D (1997) Nidovirales: a new order comprising coronaviridae and arteriviridae. Arch Virol 142:629–633
Meng XJ, Paul PS, Halbur PG, Lum MA (1995) Phylogenetic analyses of the putative M (ORF 6) and N (ORF 7) genes of porcine reproductive and respiratory syndrome virus (PRRSV): implication for the existence of two genotypes of PRRSV in the USA and Europe. Arch Virol 140:745–755
Murtaugh MP, Faaberg KS, Laber J, Elam M, Kapur V (1998) Genetic variation in the PRRS virus. Adv Exp Med Biol 440:787–794
Kim WI, Lee DS, Johnson W, Roof M, Cha SH, Yoon KJ (2007) Effect of genotypic and biotypic differences among PRRS viruses on the serologic assessment of pigs for virus infection. Vet Microbiol 123:1–14
Stadejek T, Oleksiewicz MB, Scherbakov AV, Timina AM, Krabbe JS, Chabros K, Potapchuk D (2008) Definition of subtypes in the European genotype of porcine reproductive and respiratory syndrome virus: nucleocapsid characteristics and geographical distribution in Europe. Arch Virol 153:1479–1488
Darwich L, Díaz I, Mateu E (2010) Certainties, doubts and hypotheses in porcine reproductive and respiratory syndrome virus immunobiology. Virus Res 154:123–132
Darwich L, Gimeno M, Sibila M, Diaz I, de la Torre E, Dotti S, Kuzemtseva L, Martin M, Pujols J, Mateu E (2011) Genetic and immunobiological diversities of porcine reproductive and respiratory syndrome genotype I strains. Vet Microbiol 150:49–62
Gimeno M, Darwich L, Diaz I, de la Torre E, Pujols J, Martín M, Inumaru S, Cano E, Domingo M, Montoya M, Mateu E (2011) Cytokine profiles and phenotype regulation of antigen presenting cells by genotype-I porcine reproductive and respiratory syndrome virus isolates. Vet Res 42:9
Díaz I, Darwich L, Pappaterra G, Pujols J, Mateu E (2006) Different European-type vaccines against porcine reproductive and respiratory syndrome virus have different immunological properties and confer different protection to pigs. Virology 351:249–259
Díaz I, Gimeno M, Darwich L, Navarro N, Kuzemtseva L, Lopez S, Galindo I, Segales J, Martin M, Pujols J, Mateu E (2012) Characterization of homologous and heterologous adaptive immune responses in porcine reproductive and respiratory syndrome virus infection. Vet Res 43:30
Albina E, Piriou L, Hutet E, Cariolet R, L’Hospitalier R (1998) Immune responses in pigs infected with porcine reproductive and respiratory syndrome virus (PRRSV). Vet Immunol Immunopathol 61:49–66
Buddaert W, Van Reeth K, Pensaert M (1998) In vivo and in vitro interferon (IFN) studies with the porcine reproductive and respiratory syndrome virus (PRRSV). Adv Exp Med Biol 440:461–467
Van den Heuvel MM, Tensen CP, van As JH, Van den Berg TK, Fluitsma DM, Dijkstra CD, Döpp EA, Droste A, Van Gaalen FA, Sorg C, Högger P, Beelen RH (1999) Regulation of CD163 on human macrophages: cross-linking of CD163 induces signaling and activation. J Leukoc Biol 66:858–866
Beura LK, Sarkar SN, Kwon B, Subramaniam S, Jones C, Pattnaik AK, Osorio FA (2010) Porcine reproductive and respiratory syndrome virus nonstructural protein 1beta modulates host innate immune response by antagonizing IRF3 activation. J Virol 84:1574–1584
Subramaniam S, Kwon B, Beura LK, Kuszynski CA, Pattnaik AK, Osorio FA (2010) Porcine reproductive and respiratory syndrome virus non-structural protein 1 suppresses tumor necrosis factor-alpha promoter activation by inhibiting NF-κB and Sp1. J Virol 406:270–279
Chaung HC, Chen CW, Hsieh BL, Chung WB (2010) Toll-like receptor expressions in porcine alveolar macrophages and dendritic cells in responding to poly IC stimulation and porcine reproductive and respiratory syndrome virus (PRRSV) infection. Comp Immunol Microbiol Infect Dis 33:197–213
Kuzemtseva L, de la Torre E, Martín G, Soldevila F, Ait-Ali T, Mateu E, Darwich L (2014) Regulation of toll-like receptors 3, 7 and 9 in porcine alveolar macrophages by different genotype 1 strains of porcine reproductive and respiratory syndrome virus. Vet Immunol Immunopathol 158:189–198
Mateu E, Diaz I (2008) The challenge of PRRS immunology. Vet J 177:345–351
Martínez-Lobo FJ, Díez-Fuertes F, Simarro I, Castro JM, Prieto C (2011) Porcine reproductive and respiratory syndrome virus isolates differ in their susceptibility to neutralization. Vaccine 29:6928–6940
Meier WA, Galeota J, Osorio FA, Husmann RJ, Schnitzlein WM, Zuckermann FA (2003) Gradual development of the interferon-gamma response of swine to porcine reproductive and respiratory syndrome virus infection or vaccination. Virology 309:18–31
Díaz I, Darwich L, Pappaterra G, Pujols J, Mateu E (2005) Immune responses of pigs after experimental infection with a European strain of Porcine reproductive and respiratory syndrome virus. J Gen Virol 86:1943–1951
Petry DB, Holl JW, Weber JS, Doster AR, Osorio FA, Johnson RK (2005) Biological responses to porcine respiratory and reproductive syndrome virus in pigs of two genetic populations. J Anim Sci 83:1494–1502
Lewis CR, Ait-Ali T, Clapperton M, Archibald AL, Bishop S (2007) Genetic perspectives on host responses to porcine reproductive and respiratory syndrome (PRRS). Viral Immunol 20:343–358
Ait-Ali T, Wilson AD, Westcott DG, Clapperton M, Waterfall M, Mellencamp MA, Drew TW, Bishop SC, Archibald AL (2007) Innate immune responses to replication of porcine reproductive and respiratory syndrome virus in isolated swine alveolar macrophages. Viral Immunol 20:105–118
Ait-Ali T, Wilson AD, Carré W, Westcott DG, Frossard JP, Mellencamp MA, Mouzaki D, Matika O, Waddington D, Drew TW, Bishop SC, Archibald AL (2011) Host inhibits replication of European porcine reproductive and respiratory syndrome virus in macrophages by altering differential regulation of type-I interferon transcriptional response. Immunogenetics 63:437–448
Boddicker N, Waide EH, Rowland RR, Lunney JK, Garrick DJ, Reecy JM, Dekkers JC (2012) Evidence for a major QTL associated with host response to porcine reproductive and respiratory syndrome virus challenge. J Anim Sci 90:1733–1746
Badaoui B, Tuggle CK, Hu Z, Reecy JM, Ait-Ali T, Anselmo A, Botti S (2013) Pig immune response to general stimulus and to porcine reproductive and respiratory syndrome virus infection: a meta-analysis approach. BMC Genomics 14:220
O’Connor D, Pollard AJ (2013) Characterizing vaccine responses using host genomic and transcriptomic analysis. Clin Infect Dis 57:860–869
Querec TD, Akondy RS, Lee EK, Cao W, Nakaya HI, Teuwen D, Pirani A, Gernert K, Deng J, Marzolf B, Kennedy K, Wu H, Bennouna S, Oluoch H, Miller J, Vencio RZ, Mulligan M, Aderem A, Ahmed R, Pulendran B (2009) Systems biology approach predicts immunogenicity of the yellow fever vaccine in humans. Nat Immunol 10:116–125
Nakaya HI, Wrammert J, Lee EK, Racioppi L, Marie-Kunze S, Haining WN, Means AR, Kasturi SP, Khan N, Li GM, McCausland M, Kanchan V, Kokko KE, Li S, Elbein R, Mehta AK, Aderem A, Subbarao K, Ahmed R, Pulendran B (2011) Systems biology of vaccination for seasonal influenza in humans. Nat Immunol 12:786–795
Zuckermann FA, Husmann RJ, Schwartz R, Brandt J, Mateu De Antonio E, Martin S (1998) Interleukin-12 enhances the virus-specific interferon gamma response of pigs to an inactivated pseudorabies virus vaccine. Vet Immunol Immunopathol 63:57–67
Yoon IJ, Joo HS, Goyal SM, Molitor TW (1994) A modified serum neutralization test for the detection of antibody to porcine reproductive and respiratory syndrome virus in swine sera. J Vet Diagn Invest 6:289–292
Jusa ER, Inaba Y, Kouno M, Hirose O, Shibata I, Kubota M, Yasuhara H (1996) Slow-reacting and complement-requiring neutralizing antibody in swine infected with porcine reproductive and respiratory syndrome (PRRS) virus. J Vet Med Sci 58:749–753
Freeman TC, Ivens A, Baillie JK, Beraldi D, Barnett MW, Dorward D, Downing A, Fairbairn L, Kapetanovic R, Raza S, Tomoiu A, Alberio R, Wu C, Su AI, Summer KM, Tuggle CK, Archibald AL, Hume DA (2012) A gene expression atlas of the domestic pig. BMC Biol 10:90
ArrayExpress – functional genomics data. http://www.ebi.ac.uk/arrayexpress. Accession number E-MTAB-3237. Accessed 14 Aug 2015
Theocharidis A, van Dongen S, Enright AJ, Freeman TC (2009) Network visualization and analysis of gene expression data using BioLayout express (3D). Nat Protoc 4:1535–1550
DAVID Bioinformatics Resources 6.7. http://david.abcc.ncifcrf.gov. Accessed 14 Aug 2015
Ingenuity. www.ingenuity.com. Accessed 14 Aug 2015
Li X, Jin X, Zhou X, Wang X, Shi D, Xiao Y, Bi D (2014) Pregnane X receptor is required for IFN-α-mediated CYP3A29 expression in pigs. Biochem Biophys Res Commun 445:469–474
Marques J, Anwar J, Eskildsen-Larsen S, Rebouillat D, Paludan SR, Sen G, Williams BR, Hartmann R (2008) The p59 oligoadenylate synthetase-like protein possesses antiviral activity that requires the C-terminal ubiquitin-like domain. J Gen Virol 89:2767–2772
Renukaradhya GJ, Meng XJ, Calvert JG, Roof M, Lager KM (2015) Live porcine reproductive and respiratory syndrome virus vaccines: current status and future direction. Vaccine 33:4069–4080
Mair KH, Koinig H, Gerner W, Höhne A, Bretthauer J, Kroll JJ, Roof MB, Saalmüller A, Stadler K, Libanova R (2015) Carbopol improves the early cellular immune responses induced by the modified-life vaccine Ingelvac PRRS® MLV. Vet Microbiol 176:352–357
Pileri E, Gibert E, Soldevila F, García-Saenz A, Pujols J, Diaz I, Darwich L, Casal J, Martín M, Mateu E (2015) Vaccination with a genotype 1 modified live vaccine against porcine reproductive and respiratory syndrome virus significantly reduces viremia, viral shedding and transmission of the virus in a quasi-natural experimental model. Vet Microbiol 175:7–16
Murtaugh MP, Genzow M (2011) Immunological solutions for treatment and prevention of porcine reproductive and respiratory syndrome (PRRS). Vaccine 29:8192–8204
Haralambieva IH, Ovsyannikova IG, Pankratz VS, Kennedy RB, Jacobson RM, Poland GA (2013) The genetic basis for interindividual immune response variation to measles vaccine: new understanding and new vaccine approaches. Expert Rev Vaccines 12:57–70
Boddicker NJ, Garrick DJ, Rowland RR, Lunney JK, Reecy JM, Dekkers JC (2014) Validation and further characterization of a major quantitative trait locus associated with host response to experimental infection with porcine reproductive and respiratory syndrome virus. Anim Genet 45:48–58
Koltes JE, Fritz-Waters E, Eisley CJ, Choi I, Bao H, Kommadath A, Serão NV, Boddicker NJ, Abrams SM, Schroyen M, Loyd H, Tuggle CK, Plastow GS, Guan L, Stothard P, Lunney JK, Liu P, Carpenter S, Rowland RR, Dekkers JC, Reecy JM (2015) Identification of a putative quantitative trait nucleotide in guanylate binding protein 5 for host response to PRRS virus infection. BMC Genomics 16:412
Gol S, Estany J, Fraile LJ, Pena RN (2015) Expression profiling of the GBP1 gene as a candidate gene for porcine reproductive and respiratory syndrome resistance. Anim Genet 46:599–606
Martínez-Lobo FJ, de Lome LC, Díez-Fuertes F, Segalés J, García-Artiga C, Simarro I, Castro JM, Prieto C (2013) Safety of Porcine Reproductive and Respiratory Syndrome Modified Live Virus (MLV) vaccine strains in a young pig infection model. Vet Res 44:115
Hitomi Y, Kawashima M, Aiba Y, Nishida N, Matsuhashi M, Okazaki H, Nakamura M, Tokunaga K (2015) Human primary biliary cirrhosis-susceptible allele of rs4979462 enhances TNFSF15 expression by binding NF-1. Hum Genet 134:737–747
Guo XK, Zhang Q, Gao L, Li N, Chen XX, Feng WH (2013) Increasing expression of microRNA 181 inhibits porcine reproductive and respiratory syndrome virus replication and has implications for controlling virus infection. J Virol 87:1159–1171
Gao L, Guo XK, Wang L, Zhang Q, Li N, Chen XX, Wang Y, Feng WH (2013) MicroRNA 181 suppresses porcine reproductive and respiratory syndrome virus (PRRSV) infection by targeting PRRSV receptor CD163. J Virol 87:8808–8812
Elass-Rochard E, Rombouts Y, Coddeville B, Maes E, Blervaque R, Hot D, Kremer L, Guérardel Y (2012) Structural determination and Toll-like receptor 2-dependent proinflammatory activity of dimycolyl-diarabino-glycerol from Mycobacterium marinum. J Biol Chem 287:3432–3444
The authors declare that they have no competing interests.
ID, FS, EC, ADW, YL, LD contributed to the laboratory analysis. BG analyzed the data using BioLayout Express3D analysis. EM contributed to the study design, discussion of results and manuscript preparation. LD, ALA and TAA contributed to the study design, data analysis, discussion and manuscript preparation. All authors read and approved the final manuscript.
This work was supported by funds from the European Project PoRRSCON-EU (FP7-KBBE-2009-3-245141), the COST-Action (FA-0902) and the Spanish Ministry Science and Innovation (AGL2008-05708-C02-02). LD has been supported by travel grants from the “Agència de Gestió d’Ajuts Universitaris i de Recerca” of the “Generalitat de Catalunya” (2011-BE-100125) and the “Ministerio de Educación, Cultura y Deporte” (Programa Estatal del Talento y su Empleabilidad en I+D+I, Subprograma Estatal de Movilidad, del Plan Estatal de Investigación Científica y Técnica y de Innovación 2013–2016: CAS14/00139). ADW and TAA were supported by BBSRC Institute Strategic Programme grants. We are indebted for the help provided by ARK-Genomics (now part of Edinburgh Genomics) for the microarray hybridizations and support for the data analysis using Partek Software.
Additional file 1. List of regulated transcripts with P value < 0.001, FDR < 0.2 and fold change above 2 and below -2 in PBMCs of HR and LR pigs are displayed. The lists of transcripts were generated using Partek software and only non-redundant records are shown HR and LR vaccinated pigs and HR and LR vaccinated pigs in presence of PHA. To enhance the comparison between the HR and LR datasets P value and fold change are shown across multiple experimental groups.
Additional file 2. Validation of the microarray data using real time qPCR. A subset of 9 regulated transcripts was reassessed using real-time qPCR. Data were used to run a regression analysis between the microarray and the real time qPCR experiment and generate the correlation coefficients R2.
Additional file 3. Gene ontology analysis of the data using DAVID bioinformatics resources . Analysis was run to generate the most likely gene ontology term for the biological processes, cellular components and molecular functions for each group of regulated transcripts in HR and LR vaccinated pigs based on P-value, number of gene component and fold enrichment. Only data with FDR <5% are shown.
Additional file 4. Upstream regulators analysis using Ingenuity Pathway Analysis for the HR and LR vaccinated pigs. As described in material and methods section the expression dataset for HR and LR pigs were analysed independently with Fisher’s exact test. To create this dataset, activation z-score| ≥2 for HR and LR pigs were generated and ranked from the highest to the lowest. Data were considered significant with P value < 0.05 and ¦activation Z-score| ≥2. The activation z-score was used to infer the state of activation of upstream regulators based on a comparison with a model that assigns random regulations. Target molecules are shown for each upstream regulator.
About this article
Cite this article
Ait-Ali, T., Díaz, I., Soldevila, F. et al. Distinct functional enrichment of transcriptional signatures in pigs with high and low IFN-gamma responses after vaccination with a porcine reproductive and respiratory syndrome virus (PRRSV). Vet Res 47, 104 (2016). https://doi.org/10.1186/s13567-016-0392-3
- False Discovery Rate
- High Responder
- Guanylate Binding Protein
- Biolayout Express3D
- High Responder Group