Comparative proteomic analysis reveals different responses in porcine lymph nodes to virulent and attenuated homologous African swine fever virus strains
Veterinary Research volume 49, Article number: 90 (2018)
African swine fever (ASF) is a pathology of pigs against which there is no treatment or vaccine. Understanding the equilibrium between innate and adaptive protective responses and immune pathology might contribute to the development of strategies against ASFV. Here we compare, using a proteomic approach, the course of the in vivo infection caused by two homologous strains: the virulent E75 and the attenuated E75CV1. Our results show a progressive loss of proteins by day 7 post-infection (pi) with E75, reflecting tissue destruction. Many signal pathways were affected by both infections but in different ways and extensions. Cytoskeletal remodelling and clathrin-endocytosis were affected by both isolates, while a greater number of proteins involved on inflammatory and immunological pathways were altered by E75CV1. 14-3-3 mediated signalling, related to immunity and apoptosis, was inhibited by both isolates. The implication of the Rho GTPases by E75CV1 throughout infection is also evident. Early events reflected the lack of E75 recognition by the immune system, an evasion strategy acquired by the virulent strains, and significant changes at 7 days post-infection (dpi), coinciding with the peak of infection and the time of death. The protein signature at day 31 pi with E75CV1 seems to reflect events observed at 1 dpi, including the upregulation of proteosomal subunits and molecules described as autoantigens (vimentin, HSPB1, enolase and lymphocyte cytosolic protein 1), which allow the speculation that auto-antibodies could contribute to chronic ASFV infections. Therefore, the use of proteomics could help understand ASFV pathogenesis and immune protection, opening new avenues for future research.
African swine fever virus (ASFV) is the causal agent of a haemorrhagic and often-lethal porcine disease, African swine fever (ASF), which causes affected countries important economic losses. There is no vaccine available against the disease, albeit promising developments for future implementation are being currently developed . ASF may range from an acute, highly lethal infection to subclinical chronic forms, depending on a complex contribution of viral and host factors . The pig immune response to ASFV has been widely studied [3, 4], showing that the virus has effective mechanisms of evading pig defensive systems, thus contributing to the immune pathology observed during acute ASF, and to virus persistence in its hosts . Studies about virus-cells interaction have contributed significantly to unravel the mechanisms involved in pig response [3, 6,7,8,9,10]. In this regard, it has been shown that the ASFV genome encodes a large number of genes that have been identified as playing a role in host immune evasion including: interferon (IFN) inhibition by several multigene family members , the NF-κB and NFAT inhibitor A238L or the apoptosis inhibitor A179L, among others. All these aspects have recently been reviewed . In addition, it is known that ASFV controls host transcription and cellular machinery of protein synthesis  thus contributing to pathology. The high complexity of ASFV, together with its tropism for the immune system of the host, complicates the understanding of ASF pathogenesis. Both ASFV-specific antibodies  and CD8+ T-cells [3, 14], have been postulated as key players in the protection against ASFV. The dual role of the immune system during ASFV infection becomes evident again during chronic ASF-infections, characterized by mild clinical signs associated, on occasions, with immunopathological processes, such as immune complexes and swine IgM, IgG and C1q depositions . Swine macrophages are the main target for ASFV, and depending on multiple factors, including the virulence of the ASFV strain, the immune system can play a dual role during both, ASFV-protection but also contributing to the virus pathogenesis. As mentioned above, virulent ASFV strains can evade the early recognition of the immune system, strategies that facilitate their replication and in vivo dissemination that, in the last stages of acute infection, provoke tissue destruction, leukopenia and total dysregulation of the immune system, reflected with the typical storm of cytokines thus contributing to ASF acute pathogenesis . Conversely, infections with attenuated ASFV strains normally render subclinical infections that are rapidly recognized by the innate immune system and cleared from the body, yielding pigs capable of resisting homologous lethal challenge .
Here we extend these studies by presenting a comparative proteomic analysis using gastrohepatic lymph nodes (GLN) from pigs infected in vivo with either the attenuated E75CV1 strain or with E75, its parental virulent ASFV strain, at different days post-infection (dpi); a model system previously used to unmask some of the mechanisms involved in homologous protection against ASFV . Since ASFV-infection modifies the expression patterns of host cell proteins, the application of proteomic approaches might help to clarify the intrinsic mechanism involved in ASFV-host interactions, as has been demonstrated [10, 17]. These methodologies have contributed to better understandings of how some viruses reprogram cell resources for their own benefit, or avoid host defensive mechanisms to survive in their host , but a limited number of proteomic studies have been carried out with swine pathogens. Thus, proteomic strategies have been used to study the differences observed after in vitro infection of Vero cells with virulent or attenuated strains of porcine epidemic diarrhoea virus , or PK-15 cells infected with classical swine fever virus . Additionally, our group has previously described the protein profile of lymph nodes of piglets inoculated with Porcine circovirus type 2 (PCV2) . Two proteomic studies have been published so far for ASFV infection: the first one focussed on proteins differentially expressed after ASFV-macrophage infection, although they were not specifically identified  and the second one, using Vero cells infected with the 608VR13 ASFV isolate, identifying proteins involved in apoptosis or in the transcription modulation .
In this work, we have conducted the first 2-DE proteomic approach and bioinformatic data analysis using lymph node tissue from pigs infected with two homologous ASFV strains with distinct virulence: the tissue culture adapted E75CV1, an attenuated virus, and its parental virulent E75 virus. In order to reach a deeper understanding of pig-ASFV interactions and swine immune responses, samples were obtained at different times post-infection. This large- scale proteomic study should provide a complete view of the major and important functions and pathways altered during the infection process, and hopefully, it could be used in the future to develop anti-ASFV strategies.
Materials and methods
Experimental design for the in vivo experiment
The in vivo experiment has been previously described and the animals and samples used in this analysis are the same as those used in experiment 2 of our previous study . Briefly, a group of 24 pigs was infected with 104 HAU50 of E75CV1, a second group of 18 pigs was infected using the same dose and route of the virulent E75 and finally a third group of 12 pigs remained uninfected (control group). Pigs from the three groups were sacrificed at days 1, 3 and 7 pi (6 or 3 per group and day) and 6 pigs from the attenuated E75CV1 and 3 from the control group were also sacrificed at day 31 pi. Lymph node samples of all animals (54 pigs) were analysed by 2-DE and mass spectrometry (see below).
Two-dimensional gel electrophoresis (2-DE) and image analysis
GLN were subjected to mechanical dissociation in sample buffer (7 M urea, 2 M thiourea, 4% w/v CHAPS, 1% w/v DTT, 0.8% ampholytes, 0.2 M PMSF) by scraping and gentle squeezing. Each supernatant was recovered and, after clean up precipitation, lymph node protein extracts for each condition (control, E75-infected or E75CV1-infected) and day (1, 3, 7 and 31 dpi) were pooled and analyzed by 2-DE following previously published methods . Briefly, protein extracts were diluted in rehydration buffer (7 M urea, 2 M thiourea, 4% w/v CHAPS, 20 mM DTT, 0.5% Triton X-100, 0.5% ampholytes) and immobilized pH gradient strips of 17 cm (5–8 linear pH gradient, Bio-Rad) were rehydrated with 300 μL (500 μg) of each protein solution and focused in a PROTEAN IEF Cell (Bio-Rad) using the following parameters: (1) active rehydration at 50 V for 12 h; (2) at 250 V for 15 min without pause after rehydratation; (3) rapid ramp until reaching 10 000 V/h and (4) until 80 000 V/h with slow ramped voltage. Previously, protein profile was analyzed by 2-DE on 3–10 pH gradient strips, finding that the resolved spots were concentrated in the 5–8 pH range. Second dimension was performed on 12% SDS–polyacrylamide gels using PROTEAN PlusTM Dodeca (Bio-Rad). Four replicated for each of the conditions (virulent, attenuated or control) were analyzed simultaneously for each dpi.
Gels were stained with SYPRO Ruby protein gel stain (Bio-Rad). Gel images were digitized with the FX Pro Plus Multi imager system (Bio-Rad) and analyzed with the PD Quest version 7.3 software (Bio-Rad). Spots detected by the program were matched between each gel in each group. Normalized protein spot volume (area multiplied by stain intensity) was calculated for each spot in the control sample and compared to its counterpart (sample from ASFV-infected animals). Intensity data were used to calculate differences in protein expression between groups (controls vs. E75 or E75CV1 infected samples), for each dpi, using a Student’s t test (p < 0.05) (two-tailed with unequal variances) after checking normality by the Wilks–Shapiro test. Differentially expressed spots were selected for protein identification by mass spectrometry.
Protein identification by mass spectrometry
Spots were automatically excised in a ProPic station (Genomic Solutions) and digested with modified porcine trypsin (sequencing grade; Promega), by using a ProGest digestion station (Genomic Solutions), as previously described . Briefly, gel pieces were destained with ammonium bicarbonate/acetonitrile (ACN), and subsequently subjected to dehydration/rehydration cycles and dried. Gel pieces were digested with trypsin, peptides extracted with trichloroacetic acid and desalted and concentrated by using μC-18 ZipTip columns (Millipore) in a ProMS station (Genomic Solutions). After that were loaded onto a 4800 MALDI-TOF/TOF Analyzer (Applied Biosystems) in automatic mode with the following setting: for the Mass Spectrometry (MS) data, m/z range 800–4000 with an accelerating voltage of 20 kV, peak density of maximum 50 peaks per 200 Da, minimal S/N ratio of 10 and maximum peak at 65. Spectra were internally calibrated with peptides from trypsin autolysis (M + H+= 842.509, M + H+= 2211.104). For the MS/MS data, fragment selection criteria were a minimum signal/noise (S/N) ratio of 5, a maximum number of peaks set at 65 and peak density of maximum 50 peaks per 200 Da. For each precursor selected for MS/MS analysis, fragment mass values in the range from 60 to 10 Da below precursor mass were used to peptide identification. Protein identification was assigned by peptide mass fingerprinting and confirmed by MS/MS analysis of at least three peptides in each sample. Mascot 1.9 search engine (Matrixscience) was used for protein identification running on GPS software (Applied Biosystems) against the NCBI mammalian database (updated monthly). Only those proteins with a significant protein score (p < 0.05) according to Mascot were taken into account in subsequent analysis.
Functional analysis of the proteins
Differentially regulated proteins were analyzed through the use of Ingenuity Pathway Analysis software (IPA, Ingenuity® Systems). This system identifies the functions and canonical pathways that are most significant to the data set, supported by canonical information stored in the Ingenuity Knowledge Base. Fisher’s exact test was used to calculate a p value determining the probability that the association between the proteins in the dataset and the functions and canonical pathways is explained by chance alone. These p-values are calculated based on the number of proteins that participate in a given function or pathway relative to the total number of occurrences of proteins in functions or pathway annotations stored in the Ingenuity Pathways Knowledge Base . The whole Ingenuity Knowledge Base was considered as reference set.
Two-dimensional electrophoresis reveals differential protein expression kinetics after in vivoinfection with either virulent E75 or attenuated E75CV1 ASFV-isolates
To widen and deepen our understanding of the mechanisms underlying ASFV-pathogenesis and immune response to ASFV, we characterized the variation in GLN protein profiles of pigs infected with E75 (high virulence) and E75CV (low virulence) at different times post-infection using a 2-DE proteomics approach. Our previous study  allowed to verify homogeneity between animals and that no outlier there was among the pigs studied. In these conditions, the technical variation can be the dominate source of variation . To overcome the inherent experimental variations of this technique, we performed simultaneously four replicate (four gels) of each day and condition and more than 85% of the spots were common between gels. After image analysis and visual confirmation of gels, around 800 spots were detected. Paired analyses between infected and control groups detected 80, 94, 62 and 39 differently changed spots with E75CV1 attenuated strain at 1, 3, 7 and 31 dpi respectively, corresponding to 42, 72, 53 and 33 different proteins. With E75 virulent isolate 57, 83 and 34 spots result differentially expressed at 1, 3 years, 7 dpi, corresponding to 33, 60 and 24 different proteins. In some cases, multiple spots were unambiguously identified as the same protein, it could arrive due to post-translational modification (crucial in the control of numerous regulatory pathways, degradation of proteins, biochemical alterations or pathogenesis), different isoforms derived from different genes of a multigen family, proteolytic damage or chemical modification of protein during sample preparation.
Differences in protein profiles were not only observed between virulent and attenuated ASFV strains, but also were shown between the different times post-infection tested. The most dramatic changes found affected lymph nodes from pigs infected with E75 at day 7pi, coinciding with the late phase of the E75 lethal virus infection, probably reflects the massive tissue destruction observed at this time post-infection, with only a few intact cells being present, mostly corresponding to infiltrates of ASFV-infected macrophages . GLNs from day 7 E75-infected pigs showed most proteins downregulated while infection with the attenuated E75CV1 virus upregulates a larger number of proteins than the virulent E75 strain at days 1, 3 and 7 pi, respectively (Table 1). As expected, there are a considerable number of spots in common that change after infection with the two isolates (36, 62 and 23 at 1, 3 and 7 dpi respectively), while the rest became altered only during the infection with one of the isolates (Additional files 1, 2, 3, 4 show the changes of expression of all proteins altered after infection and additional MS/MS information is available in Ref. ). The proteins with the largest change in fold change are shown in Tables 2, 3, 4, 5. Similarly, the expression of some proteins became apparent at specific times post-infection, while others were affected throughout the infection. As good examples, albumin and cytoskeletal proteins such as actin-like proteins were commonly affected by both ASFV strains, while many other proteins were differentially affected, such as heat shock protein B1 (HSPB1), involved in stabilizing actin filaments after stress, that results up-regulated exclusively at 1 dpi with E75CV1 . Of special interest were the proteins differentially regulated after the infection with the E75CV1 attenuated strain, reflecting changes in the immune system. Interestingly, some of these proteins became upregulated as soon as day 1 pi with E75CV1 (Additional file 1), including proteasome activator complex (PSME 1 and 2). Two proteins were inversely regulated after infection with the attenuated E75CV1 strain: HSPA1B and vimentin. While HPSA1B showed a downregulation, vimentin was upregulated at day 1 pi. Other proteins specifically upregulated at day 1 pi with E75CV1 were HSPB1, enolase or lymphocyte cytosolic protein 1 (LCP-1) (Additional file 1). As expected, the early events described at day 1 pi with E75CV1, correlate with the profile of upregulated proteins specifically found by day 31 pi and related with the immune response (Additional file 4), including immunoglobulin component fractions, such as IGKC and IGHG3, or components of the proteasome such as PMSA1 and PMSA6. The up-regulation of retinol-binding protein 4 (responsible for retinol transport), galactose mutarotase or calreticulin are also of interest.
Functional analysis of differentially expressed proteins after in vivo infection with either virulent E75 or attenuated E75CV1 ASFV-isolates
Bioinformatic tools were employed to biologically interpret the data set of protein obtained from 2-DE analysis, with the aim of gaining an insight into biological functions and pathways associated with the proteome response of porcine GLN to ASFV, as well as to discover differences in processes that occur after infection with each one of the ASFV-strains. Therefore, we have analyzed our data set using IPA, focused on cell functions and on canonical pathways, so-called because they contain well-established knowledge about specific relationships between groups of proteins. Differentially expressed proteins were involved in different aspects of the host–pathogen interaction (Additional files 5, 6, 7, 8, 9, 10, 11). With respect to functions, inflammatory and immunological disease were functions altered by both virus, but a greater number of proteins involved in these functions were differentially expressed after E75CV1 than E75 infection, at 1 dpi. Thus, we have found the up-regulation of proteins such as: HSPA5 and HSPB1, LCP1, vimentin, C reactive protein or hemopexin, a protein that positively regulates the interferon-gamma-mediated signaling pathway . Accordingly, our previous qPCR results  showed an over-expression of interferon-gamma after infection with E75CV1 but not due to E75 infection, at 1 dpi. Conversely, GLNs from E75-infected pigs showed, at day 1 pi, the specific downregulation of SERPINA3, an acute phase protein that is induced during inflammation  or vitamin D-binding protein (implicated in macrophage activation and inflammation) , most probably contributing to a delayed inflammatory response to the virulent isolate.
Figure 1 also included other functions (derived from the Ingenuity Pathways Analysis) significantly altered at early times after infection with both virus strains, e.g. those associated to free radical scavenging, cell death and cell-to-cell signaling and interaction, while cellular assembly and organization or cellular functions showed fewer modifications. Nevertheless, by day 7 pi with E75CV1, coinciding with the recovery of the infected pigs , these latter functions became more relevant. Conversely, GLN from E75-infected pigs by day 7 pi, showed protein profiles enriched in functions related to tissue destruction when compared to E75CV1; these were mainly associated with connective tissue disorders, skeletal and muscular disorders or with organismal injury and abnormalities. Some of these abnormalities became evident as early as at day 3 pi, coinciding with the ASFV replication in the GLN .
The association of the proteins significantly affected by ASFV infection with canonical pathways is shown in Table 6. As expected, the remodelling of epithelial adherent junction pathways was enriched for both viruses. In this canonical pathway, several members of the actin family participate, shown in here as differentially regulated. Also, clathrin-mediated endocytosis signaling pathways, involving molecules such as actin, heat shock proteins, serpin or transferrin, were regulated at all times post-infection with the exception of day 1 pi with E75CV1. One of the most remarkable and novel observations of our work affects the 14-3-3 mediated signaling, inhibited by both isolates at 1 and 3 dpi. It is also worth highlighting the differential regulation of Rho GTPases in GLNs from E75CV1 infected pigs between days 1 and 7, coinciding with the induction of innate immunity and the resolution of the subclinical ASFV infection . Regarding pathways regulated by only one virus, or in some cases at only times post-infection, is worthy to mention the G2/M DNA damage checkpoint regulation, involving the different 14.3.3 isoforms which we have been found downregulated after infection.
Proteome approaches are being increasingly used in many different systems to investigate host–microbe interactions and it has provided important information about the protein profile of cells infected with ASFV [10, 17]. The results presented here could complement the aforementioned ones since they have been performed with infected pig tissue. In this context, in vivo models could help to reflect the multiple events undergone by the host upon pathogen infection. Here, we applied a proteomics approach, based on 2-DE coupled to an in vivo experiment, to obtain new insights about the differential responses triggered by attenuated and virulent homologous ASFV strains. Lymphnodes are complex tissues composed by different cell types and, despite ASFV main targets are the macrophages, it is known that the effect of virus infection is not only restricted to the infected cell but also to the cascade of events it provokes the massive apoptosis induced in surrounding non-infected lymphocytes, this being a key event for ASF pathogenesis, responsible for lymphoid organ impairment in acute ASFV infection . Besides, the lesions found in ours previous manuscript  after necropsy together with a comparative transcription profile of key immune mediators that were differentially modulated throughout the infections point gastrohepatic lymph node as the organ of choice to improve our knowledge about of the differential porcine response against virulent or attenuated ASFV isolated, to further compare globally the immunopathogenesis of both virus strains.
Also, we used IPA in a hypothesis generative manner aiming to unmask the most relevant functions and pathways altered throughout the infection in GLN, one of its main target organs. Overall, the analysis shows that, together with pathways involved in host-immune responses, a significant number of other host functions are modified, demonstrating the complex host-virus interactions that occur in vivo. Despite the much lower ASFV load observed after infection with the attenuated ASFV strain , major protein changes were detected early after infection with E75CV1, affecting not only the infected macrophages, but also surrounding cell-types in an indirect manner. Therefore, E75 and E75CV1 not only differed in the kinetics and in the clinical symptoms they provoke, but also in the pathways activated throughout the infections. The fewer number of proteins with virulent isolated at 1 dpi might help explain the failure of the innate immune system to detect and control the first rounds of E75-replication, thus allowing its rapid spread. Conversely, efficient regulation of the innate immune system became evident for E75CV1 as early as at 1 day pi. As example, the implication of the Rho GTPases signalling pathway by E75CV1 throughout infection, (Table 6), perhaps confirming the relevance of this pathway in the innate immune system. Although the best-known function of this protein family is regulate and coordinates of actin and microtubule cytoskeletal dynamics and adhesion , Rho GTPases regulate numerous basic cell functions including regulation of the signaling pathways and cellular responses that enable to phagocytes perform their innate immune functions to respond invading pathogens. Thus, they are key regulators of cell migration (through various cell-surface receptors as TLRs), reactive oxygen species (ROS) production by NADPH oxidase, phagocytosis and degranulations well as a essential, and perhaps unique, roles in the motile responses of leukocytes . Another studies also indicate that Rho GTPases provide alternative pathways to regulate NF-kB transcriptional activity in cells of the innate immune system .
Several other proteins became differentially regulated by day 1 pi, as vimentin, enolase, HSPB1 or lymphocyte cytosolic protein 1 (LCP-1) (Additional file 1) from E75CV1-pigs. Vimentin has been described as an important molecule during ASFV morphogenesis by changing its localization to viral factories after in vitro infection . Thus, the upregulation of vimentin, at day 1 pi, might be an indirect effect provoked by E75CV1 infection on surrounding non-infected cells. Coincidently, vimentin and other proteins specifically upregulated at day 1 pi with E75CV1, such as HSPB1, enolase or LCP-1 (Additional file 1), share both their potential role in immune defense and their description as autoantigens in autoimmune disorders . In particular, the role of HSPs as immunogenic molecules able to activate T cells has been long known  and several reports have shown that pretreatment with HSPs protects from autoimmune disease . Also, T-cell response to α-enolase could be involved in the pathogenesis of autoimmune diseases . We are currently trying to unmask the potential presence of auto-antibodies in pigs infected with E75CV1 by using sera from day 31 pi and the upregulated autoantigens found at day 1 pi. These markers might be of utility to diagnose chronic infections and/or to better understand ASFV pathogenesis during chronic infections.
Of particular interest could be the complex regulation found for several members of the heterogeneous nuclear ribonucleoprotein (hnRNP) complex, RNA-binding proteins involved in downstream gene regulation and G2/M DNA damage checkpoint regulation, apoptosis and immune regulation . hnRNPH1, a protein previously shown as being involved in virus replication , was upregulated at day 1 pi with both viruses, overexpression that was maintained for the attenuated E75CV1 virus at day 3 and 7 pi. It has been reported that upregulation of hnRNPH causes a decrease in HIV virion production . Interestingly, it is worth noting the involvement of hnRNPH expression with TNF-α production  and NF-κB activation, both proteins being involved in inflammation and immune response. In Hepatitis C virus infection, it has been proposed that hnRNP might limit the amount of viral RNA genomes available for incorporation into virus . These changes in expression could be significant and add information to the only described interaction between the ASFV p30 antigen and hnRNPK in Vero infected cells . The authors suggest that the interaction hnRNPK-p30 could contribute to the host cell shut-off and represent a possible additional mechanism by which ASFV down-regulates host cell mRNA translation.
Together with pathways involved in immune responses, it is worth highlighting, the pathways involved in cytoskeletal and epithelial adherent junction remodeling (Table 6), reflecting the use that ASFV makes of cytoskeleton from virus entry to virus morphogenesis and cell egress [9, 44]. This result perfectly correlates with the function that this pathway plays during ASFV entry in pig cells , involving molecules such as actin, heat shock proteins, serpin or transferrin, which have been found altered in our study. Similarly, the clathrin-mediated endocytosis pathways here highlighted, have already been described as essential for ASFV entry in susceptible cells , giving consistency and validity to our results. Clathrin-mediated endocytosis, a strategy used by many viruses for cell entry , is commonly activated during both attenuated and virulent ASFV infection (Table 6) from day 1 pi with E75 and from day 3 pi with the E75CV1 strain, coinciding with the differential kinetics of ASFV in vivo replication observed for both strains . The differential expression of several other proteins also involved in cytoskeleton formation it is worthy to be discussed. Thereby, the Rho GTPase family is also involved in the regulation of microtubules during dynein-mediated capsid transport of herpes virus associated Kaposi’s sarcoma , a pathway also required for ASFV entry, morphogenesis and exit from the infected cells [31, 44]. The inhibition of RhoGTPases by RhoGDI by attenuated isolate might have a negative effect in all these processes thus impairing the in vivo transmission of E75CV1. Conversely, the down-expression of RhoGDI at 7 dpi, might contribute to the successful systemic dissemination of E75. Also, the results observed with RP/EB microtubule-associated protein (which negatively regulates microtubule formation) point in this direction. Thus, interestingly, in our study this protein is down-regulated at 1 dpi with both virus (which would facilitate the formation of the microtubules and therefore traffic virus at the onset of infection) but over-expressed to 3 dpi only with attenuated isolate (inhibiting the formation of the microtubules and thus the transport of the attenuated isolate).
A chapter apart deserves discussing the sub-expression found at early times post-infection (day 1 and 3 pi) with both viruses, of several of the 7 isoforms of 14-3-3 protein, a novel finding for ASFV. 14-3-3 interactome studies have demonstrated that 14-3-3 proteins participate in many events associated with infection in other viruses , mainly activated by dsRNA. The downregulation observed for both E75CV1 and E75 in vivo might reflect the need of ASFV to evade the innate immune responses triggered by 14-3-3, including the activation of TNF and NFκB signalling or any other antimicrobial responses triggered by activating the TLR-14-3-3 pathways . While 14-3-3 has been strongly associated with the innate immunity activated in response to dsRNA , its negative effect on virus morphogenesis , and ASFV-exit from the cell, , have been also described. On the other hand, an important function of 14-3-3 proteins is to inhibit apoptosis  and downregulation of 14-3-3 might also have a direct effect on the G2/M DNA damage checkpoint regulation , a pathway that prevents cells with damaged DNA entering the M phase of cell division before repairing. So, a defective G2/M checkpoint leads DNA damaged cells to apoptosis . Curiously enough, E75CV1 downregulates the G2/M cell cycle control checkpoint both at day 1 and 3 pi, while E75 does it at day 3 pi (Table 6) as described for other viruses . The concomitant downregulation of the G2/M DNA damage checkpoint and 14-3-3 observed might contribute to differentially activate the apoptosis of the infected cells, perhaps contributing to the deficient in vivo dissemination of E75CV1 . All together, these results seem to reflect a very complex regulation of apoptosis during ASFV infection, as has been previously postulated . These results complement previous work demonstrating the effect of DNA damage and apoptosis in ASFV-in vitro replication , confirming the role that 14-3-3 play during apoptosis inhibition , as has been demonstrated for other viruses . Studies to confirm the relevance of 14-3-3 in this inhibition might help to design novel antiviral strategies.
The last part of the discussion will be dedicated to the results obtained with samples harvested at day 31 after the infection with E75CV1, that might be very useful to understand the intrinsic mechanisms involved in protection against E75-virulent challenge. As expected, the proteomic analysis performed with these samples seem to confirm the key relevance that both antibodies  and CD8 T-cells [15, 60], play in protection against ASFV. The immunoglobulin (Ig) isoforms detected could play a dual role participating either in ASFV-antibody mediated protection , or also in the formation of immune complexes found in chronically ASFV infected pigs . PMSA1 and PMS6 are involved in swine leukocyte antigen class I (SLAI) presentation and CD8-T cell induction but also have endoribonuclease activity, playing important defensive roles in response to external stimuli . The up-regulation of other proteins with implication in the antigenic presentation, such as retinol-binding protein 4 (responsible for retinol transport) and galactose mutarotase (GALM) are also of interest. Retinol is a potent regulator of B-cell receptor function and B-cell activation  and regulates GALM gene expression, which could have an important role in cell adhesion and antigen presentation . Other components of the B cell response were upregulated, including members of Ig family, proteins implicated in antigenic processing and SLAI antigen presentation (as proteasome subunits), calreticulin or serpin. Calreticulin (CRT) is a Ca2+-binding protein involved in more than 40 functions, including the unfolded protein response (UPR) and antigen presentation. Interestingly, CRT has been found upregulated in vitro as part of endoplasmic reticulum (ER) stress and the UPR provoked by ASFV in infected cell . UPR is a antiviral mechanism, against which many viruses develop multiple strategies  and on the other hand, CRT transport antigens to the ER, facilitating antigen presentation in association with the major histocompatibility complex class I (MHC-I) molecules to elicit peptide-specific CD8+ T cell responses . In accordance with the latter, several other proteins involved in antigen presentation and T-cell activation have been found up-regulated at day 31 pi with E75CV1, including SERPINB9 and proteasomal subunits. Serpin B anti-proteases have been defined as regulators of the immune response up-regulated during several virus infections, including HIV  and Epstein–Barr virus infection . SERPINB9 has been shown to be involved in protection of antigen presenting cells, enhancing T cell activation and immune response , including in IFN-γ production and antiviral cytopathic responses  and survival of CD8+ memory T cells . The proteomic data obtained with samples from E75CV1-recoverd pigs (31 dpi) perfectly fits with the increasing evidence that Th1 and specific CD8 T-cells play in protection , including our previous data using this same ASFV-infection model . Altogether, our results could indicate that the over-expression of CRT, SERPINB9 and subunits from the proteasome reflect the relevance that SLAI presentation and CD8 T-cell activation play during ASFV infection, opening new avenues to fight the disease. Interestingly, these three components have already been used as genetic adjuvants to improve the specific immunity against several pathogens . We are currently extending our studies to the field of ASFV vaccinology. Thus, today we know that targeting antigens to the proteasome improves the protection against ASFV even in the absence of antibodies [14, 60]. The knowledge gained here opens new avenues to improve these strategies in the near future.
In conclusion, the data presented here are the first to compare kinetics of protein expression profiles from pigs infected by homologous virulent or attenuated strains of ASFV, through a proteomics approach coupled with a large-scale in vivo infection, in order to allow the generation of advances in our understanding about the pig immune response to virus and pathogenesis of ASF over time. Our results confirm a differential interaction with the immune system for both viruses. Thus, GLN from E75CV1 infected pigs showed the largest number of differentially upregulated proteins as early as 1 dpi, many of them involved in the activation of different innate immune pathways, including autoantigens. In addition to a lower replication efficiency at early time post-infection by attenuated isolated, the induction of specific antibody and T-cell responses at 31 dpi, were observed, once E75CV1 has been cleared. We believe that the increased information yielded by this global approach could improve our knowledge about the major point underlying host–pathogen interactions and might be the important for the development of an efficacious ASF vaccine.
Rock DL (2017) Challenges for African swine fever vaccine development-”… perhaps the end of the beginning.”. Vet Microbiol 206:52–58
Portugal R, Coelho J, Höper D, Little NS, Smithson C, Upton C, Martins C, Leitão A, Keil GM (2015) Related strains of African swine fever virus with different virulence: genome comparison and analysis. J Gen Virol 96:408–419
Takamatsu HH, Denyer MS, Lacasta A, Stirling C, Argilaguet J, Netherton CL, Oura C, Martins C, Rodriguez F (2013) Cellular immunity in ASFV responses. Virus Res 173:110–121
Escribano JM, Galindo I, Alonso C (2013) Antibody-mediated neutralization of African swine fever virus: myths and facts. Virus Res 173:101–109
Reis AL, Netherton C, Dixon LK (2017) Unraveling the armor of a killer: evasion of host defenses by African swine fever virus. J Virol 91:e02338
Zhang F, Hopwood P, Abrams CC, Downing A, Murray F, Talbot R, Archibald A, Lowden S, Dixon LK (2006) Macrophage transcriptional responses following in vitro infection with a highly virulent African swine fever virus isolate. J Virol 80:10514–10521
Ramiro-Ibañez F, Ortega A, Ruiz-Gonzalvo F, Escribano JM, Alonso C (1997) Modulation of immune cell population and activation markers in the pathogenesis of African swine fever virus infection. Virus Res 47:31–40
Gil S, Sepúlveda N, Albina E, Leitão A, Martins C (2008) The low-virulent African swine fever virus (ASFV/NH/P68) induces enhanced expression and production of relevant regulatory cytokines (IFNalpha, TNFalpha and IL12p40) on porcine macrophages in comparison to the highly virulent ASFV/L60. Arch Virol 153:1845–1854
Alonso C, Galindo I, Cuesta-Geijo MA, Cabezas M, Hernaez B, Muñoz-Moreno R (2013) African swine fever virus-cell interactions: from virus entry to cell survival. Virus Res 173:42–57
Alfonso P, Rivera J, Hernaez B, Alonso C, Escribano J (2004) Identification of cellular proteins modified in response to African swine fever virus infection by proteomics. J Proteomics 4:2037–2046
Golding JP, Goatley L, Goodbourn S, Dixon LK, Taylor G, Netherton CL (2016) Sensitivity of African swine fever virus to type I interferon is linked to genes within multigene families 360 and 505. Virology 493:154–161
Galindo I, Alonso C (2017) African swine fever virus: a review. Viruses 9:E103
Sánchez EG, Quintas A, Nogal M, Castelló A, Revilla Y (2013) African swine fever virus controls the host transcription and cellular machinery of protein synthesis. Virus Res 173:58–75
Lacasta A, Ballester M, Monteagudo PL, Rodriguez JM, Salas ML, Accensi F, Pina-Pedrero S, Bensaid A, Argilaguet J, Lopez-Soria S, Hutet E, Le Potier MF, Rodriguez F (2014) Expression library immunization can confer protection against lethal challenge with African swine fever virus. J Virol 88:13322–13332
Fernandez A, Perez J, de Martin J, Carrasco L, Dominguez J, Sierra MA (1992) Localization of African swine fever viral antigen, swine IgM, IgG and C1q in lung and liver tissues of experimentally infected pigs. J Comp Pathol 107:81–90
Lacasta A, Monteagudo P, Jiménez-Marin A, Accensi F, Ballester M, Argilaguet J, Galindo-Cardiel I, Segalés J, Salas ML, Dominguez J, Moreno A, Garrido JJ, Rodríguez F (2015) Live attenuated African swine fever viruses as ideal tools to dissect the mechanisms involved in viral pathogenesis and immune protection. Vet Res 46:135
Rodriguez J, Salas M, Santaren J (2001) African swine fever virus induced polypeptides in porcine alveolar macrophages and in Vero cells: two-dimensional gel analysis. Proteomics 1:1447–1456
Rodrigo G, Daros JA, Elena SF (2017) Virus-host interactome: putting the accent on how it changes. J Proteomics 156:1–4
Guo X, Hu H, Chen F, Li Z, Ye S, Cheng S, Zhang M, He Q (2016) iTRAQ-based comparative proteomic analysis of Vero cells infected with virulent and CV777 vaccine strain-like strains of porcine epidemic diarrhea virus. J Proteomics 130:65–75
Sun J, Jiang Y, Shi Z, Yan Y, Guo H, He F, Tu C (2008) Proteomic alteration of PK-15 cells after infection by classical swine fever virus. J Proteome Res 7:5263–5269
Ramírez-Boo M, Nuñez E, Jorge I, Navarro P, Fernandes L, Segales J, Garrido JJ, Vázquez J, Moreno A (2011) Quantitative proteomics by 2-DE, 16O/18O labelling and linear ion trap mass spectrometry analysis of lymph nodes from piglets inoculated by porcine circovirus type 2. Proteomics 11:3452–3469
Collado-Romero M, Prado-Martins R, Arce C, Moreno A, Lucena C, Carvajal A, Garrido JJ (2012) An in vivo proteomic study of the interaction between Salmonella Typhimurium and porcine ileum mucosa. J Proteomics 75:2015–2026
Pallen C, Friry-Santini C, Herouet-Guicheney C, Capt A (2014) Technical variability of 2D gel electrophoresis—application to soybean allergens. Toxic Rep 1:734–742
Central Service for Research Support. http://www.uco.es/scai/geles_2017.html
Wettstein G, Bellaye PS, Micheau O, Bonniaud P (2012) Small heat shock proteins and the cytoskeleton: an essential interplay for cell integrity? Int J Biochem Cell Biol 44:1680–1686
Xu M, Tan C, Hu J, Alwahsh SM, Yan J, Hu J, Dai Z, Wang Z, Zhou J, Fan J, Huang X (2014) Expression of hemopexin in acute rejection of rat liver allograft identified by serum proteomic analysis. Shock 42:65–74
Janciauskiene S, Wright HT (1998) Inflammation, antichymotrypsin, and lipid metabolism: autogenic etiology of Alzheimer’s disease. BioEssays 20:1039–1046
Pacini S, Punzi T, Morucci G, Gulisano M, Ruggiero M (2012) Effects of vitamin D-binding protein-derived macrophage-activating factor on human breast cancer cells. Anticancer Res 32:45–52
Bidaud-Meynard A, Binamé F, Lagrée V, Moreau V (2017) Regulation of Rho GTPase activity at the leading edge of migrating cells by p190RhoGAP. Small GTPases 13:1–12
Ramiro-Ibañez F, Ortega A, Brun A, Escribano JM, Alonso C (1996) Apoptosis: a mechanism of cell killing and lymphoid organ impairment during acute African swine fever virus infection. J Gen Virol 77:2209–2219
Quetglas J, Hernaez B, Galindo L, Muñoz R, Cuesta M, Alonso C (2012) Small Rho GTPases and cholesterol biosynthetic pathway intermediates in African swine fever virus infection. J Virol 86:1758–1767
Bokoch GM (2005) Regulation of innate immunity by Rho GTPases. Trends Cell Biol 15:163–171
Stefanovic S, Windsor M, Nagata KI, Inagaki M, Wileman T (2005) Vimentin rearrangement during African swine fever virus infection involves retrograde transport along microtubules and phosphorylation of vimentin by calcium calmodulin kinase II. J Virol 79:11766–11775
Moy L, Levine J (2014) Autoimmune hepatitis: a classic autoimmune liver disease. Curr Probl Pediatr Adolesc Health Care 44:341–346
Zugel U, Kaufmann SH (1999) Role of heat shock proteins in protection from and pathogenesis of infectious diseases. Clin Microbiol 12:19–39
Fijak M, Iosub R, Schneider E, Linder M, Respondek K, Klug J, Meinhardt A (2005) Identification of immunodominant autoantigens in rat autoimmune orchitis. J Pathol 207:127–138
Nahm DH, Lee KL, Shin JY, Ye JM, Kang Y, Park HS (2006) Identification of alpha-enolase as an autoantigen associated with severe asthma. J Allergy Clin Immunol 118:376–381
Yabas M, Elliott H, Hoyne GF (2015) The role of alternative splicing in the control of immune homeostasis and cellular differentiation. Int J Mol Sci 17:E3
Lee JW, Liao PC, Young KC, Chang CL, Chen SS, Chang TT, Lai MD, Wang SW (2011) Identification of hnRNPH1, NF45, and C14orf166 as novel host interacting partners of the mature hepatitis C virus core protein. J Proteome Res 10:4522–4534
Jablonski JA, Caputi M (2009) Role of cellular RNA processing factors in human immunodeficiency virus type 1 mRNA metabolism, replication, and infectivity. J Virol 83:981–992
Mishra KP, Shweta Diwaker S, Ganju L (2012) Dengue virus infection induces up-regulation of hn RNP-H and PDIA3 for its multiplication in the host cell. Virus Res 163:573–579
Poenisch M, Metz P, Blankenburg H, Ruggieri A, Lee JY, Rupp D, Rebhan I, Diederich K, Kaderali L, Domingues F, Albrecht M, Lohmann V, Erfle H, Bartenschlager R (2015) Identification of HNRNPK as regulator of hepatitis C virus particle production. PLoS Pathog 11:e1004573
Hernaez B, Escribano JM, Alonso C (2008) African swine fever virus protein p30 interaction with heterogeneous nuclear ribonucleoprotein K (hnRNP-K) during infection. FEBS Lett 582:3275–3280
Netherton C, Wileman T (2013) African swine fever virus organelle rearrangements. Virus Res 173:76–86
Andrés G (2017) African swine fever virus gets undressed: new insights on the entry pathway. J Virol 91:e01906–e01916
Daecke JO, Fackler T, Dittmar MT, Krausslich HG (2005) Involvement of clathrin-mediated endocytosis in human immunodeficiency virus type 1 entry. J Virol 79:1581–1594
Naranatt P, Krishnan H, Smith M, Chandran B (2005) Kaposis sarcoma associated herpesvirus modulated microtubule dynamics via RhoA-GTP-diaphanous 2 signaling and utilizes the dynein motors to deliver its DNA to the nucleus. J Virol 79:1197–1206
Johnson C, Tinti M, Wood NT, Campbell DG, Toth R, Dubois F, Geraghty K, Wong BH, Brown LJ, Tyler J, Gernez A, Chen S, Synowsky S, MacKintosh C (2011) Visualization and biochemical analyses of the emerging mammalian 14-3-3-phosphoproteome. Mol Cell Proteomics 10(M110):005751
Ohman T, Söderholm S, Hintsanen P, Välimäki E, Lietzén N, MacKintosh C, Aittokallio T, Matikainen S, Nyman T (2014) Phosphoproteomics combined with quantitative 14-3-3-affinity capture identifies SIRT1 and RAI as novel regulators of cytosolic double-stranded RNA recognition pathway. Mol Cell Proteomics 13:2604–2617
Ohman T, Lietzen N, Valimaki E, Melchiorsen J, Matikainen S, Nyman T (2010) Cytosolic RNA recognition pathway activates 14-3-3 protein mediated signaling and caspase dependent disruption of cytokeratin network in human keratinocytes. J Proteome Res 9:1549–1564
Pei Z, Harrison M, Schmitt A (2011) Parainfluenza virus 5 m protein interaction with host protein 14-3-3 negatively affects virus particle formation. J Virol 85:2050–2059
Salas ML, Andrés G (2013) African swine fever virus morphogenesis. Virus Res 173:29–41
Xing H, Zhang S, Weinheimer C, Kovacs A (2000) 14-3-3 proteins block apoptosis and differentially regulate MAPK cascades. EMBO J 19:349–358
Samuel T, Weber H, Rauch P, Verdoodt B, Eppel JT, McShea A, Hermeking H, Funk JO (2001) The G2/M regulator 14-3-3sigma prevents apoptosis through sequestration of Bax. J Biol Chem 276:45201–45206
Cuddihy AR, O’Connell MJ (2003) Cell-cycle responses to DNA damage in G2. Int Rev Cytol 222:99–140
Mendes M, Pérez-Hernandez D, Vázquez J, Coelho A, Cunhaa C (2013) Proteomic changes in HEK-293 cells induced by hepatitis delta virus replication. J Proteomics 89:24–38
Tulman E, Delhon G, Ku B, Rock D (2009) African swine fever virus. Curr Top Microbiol Immunol 328:43–87
Granja AG, Nogal ML, Hurtado C, Salas J, Salas ML, Carrascosa AL, Revilla Y (2004) Modulation of p53 cellular function and cell death by African swine fever virus. J Virol 78:7165–7174
Chan YK, Gack MU (2016) A phosphomimetic-based mechanism of dengue virus to antagonize innate immunity. Nat Immunol 17:523–530
Argilaguet JM, Pérez-Martín E, Nofrarías M, Gallardo C, Accensi F, Lacasta A, Mora M, Ballester M, Galindo-Cardiel I, López-Soria S, Escribano JM, Reche PA, Rodríguez F (2012) DNA vaccination partially protects against African swine fever virus lethal challenge in the absence of antibodies. PLoS One 7:e40942
Kulichkova VA, Tsimokha AS, Fedorova OA, Moiseeva TN, Bottril A, Lezina L, Gauze LN, Konstantinova IM, Mittenberg AG, Barlev NA (2010) 26S proteasome exhibits endoribonuclease activity controlled by extra-cellular stimuli. Cell Cycle 9:840–849
Chen Q, Ross AC (2005) Vitamin A and immune function: retinoic acid modulates population dynamics in antigen receptor and CD38-stimulated splenic B cells. Proc Natl Acad Sci U S A 102:14142–14149
Pai T, Chen Q, Zhang Y, Zolfaghari R, Ross AC (2007) Galacto mutarotase and other galactose-related genes are rapidly induced by retinoic acid in human myeloid cells. Biochemistry 46:15198–15207
Galindo I, Hernaez B, Muñoz-Moreno R, Cuesta-Geijo MA, Dalmau-Mena I, Alonso C (2012) The ATF6 branch of unfolded protein response and apoptosis are activated to promote African swine fever virus infection. Cell Death Dis 3:e341
Fung T, Huang M, Liu DX (2014) Coronavirus-induced ER stress response and its involvement in regulation of coronavirus–host interactions. Virus Res 194:110–123
Gomez-Gutierrez J, Elpek KG, Montes de Oca-Luna R, Shirwan H, Zhou S, McMasters K (2007) Vaccination with an adenoviral vector expressing calreticulin-human papillomavirus 16 E7 fusion protein eradicates E7 expressing established tumors in mice. Cancer Immunol Immunother 56:997–1007
Aboud L, Ball TB, Tjernlund A, Burgener A (2014) The role of serpin and cystatin antiproteases in mucosal innate immunity and their defense against HIV. Am J Reprod Immunol 71:12–23
Classen CF, Bird I, Debatin KM (2006) Modulation of the granzyme B inhibitor proteinase inhibitor 9 (PI-9) by activation of lymphocytes and monocytes in vitro and by Epstein–Barr virus and bacterial infection. Clin Exp Immunol 143:534–542
Barrie MB, Stout HW, Abougergi MS, Miller BC, Thiele DL (2004) Antiviral cytokines induce hepatic expression of the granzyme B inhibitors, proteinase inhibitor 9 and serine proteinase inhibitor 6. J Immunol 172:6453–6459
Phillips T, Opferman JT, Shah R, Liu N, Froelich CJ, Ashton-Rickardt PG (2004) A role for the granzyme B inhibitor serine protease inhibitor 6 in CD8+ memory cell homeostasis. J Immunol 173:3801–3809
The authors declare that they have no competing interests.
JHU carried out the proteomic analysis. AJM performed the bioinformatic analysis. AL, SP and PLM performed the in vivo experiment, collected and processed the tissue samples for proteomic analysis. JHU, AJM and AL participated in analyses and interpretation of results, drafted and edited the manuscript. AM, FR, and JJG were in charge project design, discussion of the results, as well as writing the manuscript. All authors read and approved the final manuscript.
We thank Erena Ruiz Mora and María Jesús Navas for skillful technical assistance. We thank Kevin Dalton for editing the manuscript. Rights to use IPA bioinformatic tools was granted by the Andalusian Platform of Bioinformatics (University of Málaga, Spain). The proteomics analysis was done at the University of Córdoba Proteomics Core Facility, who are members of the ProteoRed network.
This work has been funded by the Spanish Government (Projects Reference Number AGL201022229 and AGL2017-87415-R). Anna Lacasta and Júber Herrera were financially supported by an FPU fellowship and an FPI fellowship, respectively, both from the Spanish Government. CERCA Programme from the Generalitat de Catalunya is also acknowledged. IRTA is supported by CERCA Programme/Generalitat de Catalunya.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Differentially regulated proteins in response to ASFV infection (1 dpi). Changes of expression of all proteins altered after this time after infection.
Differentially regulated proteins in response to ASFV infection (3 dpi). Changes of expression of all proteins altered after this time after infection.
Differentially regulated proteins in response to ASFV infection (7 dpi). Changes of expression of all proteins altered after this time after infection.
Differentially regulated proteins in response to ASFV infection (31 dpi). Changes of expression of all proteins altered after this time after infection.
Canonical pathways and functions significantly regulated by attenuated ASFV in porcine lymph node at 1 dpi. The analysis is derived from the Ingenuity Pathways Analysis.
Canonical pathways and functions significantly regulated by virulent ASFV in porcine lymph node at 1 dpi. The analysis is derived from the Ingenuity Pathways Analysis.
Canonical pathways and functions significantly regulated by attenuated ASFV in porcine lymph node at 3 dpi. The analysis is derived from the Ingenuity Pathways Analysis.
Canonical pathways and functions significantly regulated by virulent ASFV in porcine lymph node at 3 dpi. The analysis is derived from the Ingenuity Pathways Analysis.
Canonical pathways and functions significantly regulated by attenuates ASFV in porcine lymph node at 7 dpi. The analysis is derived from the Ingenuity Pathways Analysis.
Canonical pathways and functions significantly regulated by virulent ASFV in porcine lymph node at 7 dpi. The analysis is derived from the Ingenuity Pathways Analysis.
Canonical pathways and functions significantly regulated by attenuated ASFV in porcine lymph node at 31 dpi. The analysis is derived from the Ingenuity Pathways Analysis.
About this article
Cite this article
Herrera-Uribe, J., Jiménez-Marín, Á., Lacasta, A. et al. Comparative proteomic analysis reveals different responses in porcine lymph nodes to virulent and attenuated homologous African swine fever virus strains. Vet Res 49, 90 (2018). https://doi.org/10.1186/s13567-018-0585-z