Skip to main content

Identification of serum proteomic biomarkers for early porcine reproductive and respiratory syndrome (PRRS) infection

Abstract

Background

Porcine reproductive and respiratory syndrome (PRRS) is one of the most significant swine diseases worldwide. Despite its relevance, serum biomarkers associated with early-onset viral infection, when clinical signs are not detectable and the disease is characterized by a weak anti-viral response and persistent infection, have not yet been identified. Surface-enhanced laser desorption ionization time of flight mass spectrometry (SELDI-TOF MS) is a reproducible, accurate, and simple method for the identification of biomarker proteins related to disease in serum. This work describes the SELDI-TOF MS analyses of sera of 60 PRRSV-positive and 60 PRRSV-negative, as measured by PCR, asymptomatic Large White piglets at weaning. Sera with comparable and low content of hemoglobin (< 4.52 μg/mL) were fractionated in 6 different fractions by anion-exchange chromatography and protein profiles in the mass range 1–200 kDa were obtained with the CM10, IMAC30, and H50 surfaces.

Results

A total of 200 significant peaks (p < 0.05) were identified in the initial discovery phase of the study and 47 of them were confirmed in the validation phase. The majority of peaks (42) were up-regulated in PRRSV-positive piglets, while 5 were down-regulated. A panel of 14 discriminatory peaks identified in fraction 1 (pH = 9), on the surface CM10, and acquired at low focus mass provided a serum protein profile diagnostic pattern that enabled to discriminate between PRRSV-positive and -negative piglets with a sensitivity and specificity of 77% and 73%, respectively.

Conclusions

SELDI-TOF MS profiling of sera from PRRSV-positive and PRRSV-negative asymptomatic piglets provided a proteomic signature with large scale diagnostic potential for early identification of PRRSV infection in weaning piglets. Furthermore, SELDI-TOF protein markers represent a refined phenotype of PRRSV infection that might be useful for whole genome association studies.

Background

Porcine reproductive and respiratory syndrome (PRRS) is one of the most important infectious swine diseases throughout the world [13] and is still having, more than two decades after its emergence, major impacts on pig health and welfare (reviewed by [4]). The responsible agent is an enveloped, ca. 15 kb long positive-stranded RNA virus (PRRSV) that belongs to the Arteriviridae family [5] and that can cause late-term abortions in sows and respiratory symptoms and mortality in young or growing pigs. Once this virus has entered a herd it tends to remain present and active indefinitely causing severe economic losses and marketing problems due to high direct medication costs and considerable animal health costs needed to control secondary pathogens [6, 7].

Pigs of all ages are susceptible to this highly infectious virus, which has been shown to be present in most pigs for the first 105 days post infection [8]. However clinical manifestations vary with physiological status and age [9], as the virus uses several immune evasion ways to complicate the ability of the host to respond to the infection process [4, 10, 11]. Weaning piglets, in particular, are likely to be exposed to the infection. Although PRRSV viraemia is often asymptomatic in these piglets, their productive performance is significantly decreased. Indeed, despite being sero-negative, persistently infected piglets still harbor PRRSV and have been shown to be a source of virus for susceptible animals [12].

SELDI-TOF MS analysis allows the comparison of protein profiles obtained from a large number of diverse biological samples by combining two principles, chromatography by retention on chip surface on the basis of defined properties (e.g. charge, surface hydrophobicity, or biospecific interaction with ligands) and mass spectrometry. It is thus distinct from common non-selective techniques, such as two-dimensional polyacrilamide gel electrophoresis (2D-PAGE) and matrix-assisted laser desorption ionisation (MALDI) MS. SELDI-TOF MS has been widely used for diagnostic biomarker discovery and validation across studies in blood serum/plasma, particularly in cancer research (reviewed by [13]), but also to characterize and identify biomarkers associated with viral and other infectious diseases [1419]. The protein signatures identified by SELDI-TOF MS analysis have thus many potential applications in animal health, including early diagnosis of diseases, prediction of disease states, as well as monitoring of disease progression, recovery, and response to vaccination. Few reports have been published for livestock applications [1922].

Current needs in veterinary medicine and animal husbandry include the identification of tools that allow the early warning of diseases, especially during the incubation periods and before the onset of clinical signs. Therefore, the objective of this study was to identify by SELDI-TOF MS a proteomic profile able to differentiate PPRSV-positive from -negative weaning piglets raised in commercial farms and without clinical symptoms of the disease. We optimized the experimental conditions previously described [20] and validated 47 statistically significant discriminatory biomarkers. Among these, a combination of 14 biomarkers identified in F1 on CM10 at low focus mass permitted to correctly assign the piglets to the PPRSV-positive or PRRSV-negative groups with sensitivity and specificity of 77% and 73%, respectively.

Results

To enable identification of medium-low abundant proteins, only samples with a total content of hemoglobin lower than 4.52 μg/mL were included in the study. Total hemoglobin absorbance and the resulting hemoglobin content were calculated for all the piglet sera in both discovery (n = 50) and validation (n = 70) phases of the study [ Additional file 1: Table S1 and Additional file 2: Table S2, respectively].

Fractioning of the sera resulted in six different pH fractions; F1 = pH9, F2 = pH7, F3 = pH5, F4 = pH4, F5 = pH3, and F6 = organic solvent. The fractions F1, F4, and F6 were analyzed on the three surfaces CM10, IMAC30, and H50 at both low and high focus masses. Fractions F2 and F3 were excluded from further analyses because preliminary data with 3 serum samples showed that they still contained elevated quantities of abundant proteins (such as albumin), as well as the quality of the spectra and the number of signals detected were very low. Fraction F5 was excluded because no signals were detected.

The fractions F1, F4, and F6 on the surfaces CM10, IMAC30, and H50 showed generally good signal intensities and low coefficient of variation (CV) values (< 30%) in both the discovery and validation phases. Exceptions were fraction F1 on IMAC30 (analyzed at high focus mass) and H50 (both low and high focus masses), as well as fraction F4 on H50 (low focus mass), which were therefore excluded from further analyses.

Discovery phase

A total of 50 pig sera, 25 from PRRSV-positive and 25 from PRRSV-negative piglets were analyzed during the discovery phase of the study [ Additional file 1: Table S1].

We found a total of 785 protein peaks in the sera of all samples (Table 1). The most represented pH fraction was F6 (n = 381), followed by F4 (n = 223), and F1 (n = 181). On surface CM10 we identified 317 peaks, on IMAC30 302 peaks, and on H50 166 peaks. Furthermore, a much higher number of peaks (n = 512) was found on low mass range (1–20 kDa) compared to the high (n = 273; 20–200 kDa).

Table 1 Protein peaks identified by SELDI-TOF MS in the discovery phase of the study

Of the total 785 peaks, 200 were statistically significant (p < 0.05) and permitted to discriminate between PRRSV-positive and PRRSV-negative piglets. Discriminatory peaks were found in F1 (n = 80), F4 (n = 49), and F6 (n = 71) on the surfaces CM10 (n = 107), IMAC50 (n = 58), and H50 (n = 35), as well with low (n = 110) and high (n = 90) focus masses (Table 1).

The highest sensitivity (80%) and specificity (76%) were obtained with the 22 discriminatory peaks of F1 on CM10 at low focus mass. Higher sensitivities were found with the 18 peaks of F4 on CM10 at low focus mass (87%), the 7 peaks of F6 on CM10 at low focus mass (85%), and the 12 peaks of F6 on CM10 at high focus mass (87%), however the specificities of these peaks were lower (64%, 66%, and 66%, respectively).

Validation phase

The validation phase was performed on 35 new PRRSV-positive and 35 new PRRSV-negative piglets using the same experimental conditions applied in the discovery phase [ Additional file 2: Table S2]. Of the total 200 peaks that were significant in the discovery phase, 47 were confirmed in the validation phase (Table 2).

Table 2 Discriminatory protein peaks identified in the discovery phase and confirmed in the validation phase

In particular, 28 peaks were confirmed on CM10, 19 on IMAC30, whereas none of the peaks could be validated on the surface H50. In the 3 fractions with different pH tested, F1 contained 28 peaks, F4 3 peaks, and F6 16 peaks. A higher number of peaks (n = 36) corresponded to small peptides (acquired at low focus mass 1–20 kDa), compared to big peptides (n = 11) that were acquired at high focus mass (20–200 kDa).

The vast majority (42) of the peaks were up-regulated in PRRSV-positive piglets compared to the negative, while only 5 peaks (F1 on CM10: 5,468 and 5,536 Da; F6 on CM10: 14,843 Da; and F6 on IMAC30: 27,806 and 27,606 Da) were down-regulated (Table 2). In line with the results of the discovery phase, the combination of peaks with the highest sensitivities (77% and 64.5%) and specificities (73% and 69.7%) were found on CM10 at low focus mass with the 14 discriminatory peaks of F1 and the 6 discriminatory peaks of F6, respectively (Table 2). The correctly and incorrectly assigned piglets using these peaks are graphically illustrated in the heat map of Figure 1; part 1A shows the 14 peaks of F1 and part 1B the 6 peaks identified in F6.

Figure 1
figure 1

Heat map showing cluster analysis of the PRRSV-positive and PRRSV-negative piglets tested with the 2 combinations of discriminatory peaks that showed the highest sensitivity and specificity values. The x-axis of the heat maps shows the piglets analyzed in the validation phase (blue: PRRSV-positive; red: PRRSV-negative), while the y-axis displays the molecular weights in Dalton of the 14 significant discriminatory peaks identified in F1 (A) and the 6 peaks in F6 (B) both on the surface CM10 at low focus mass. The maps contain peak fold changes Z-score normalized over all piglets. They are color coded, with red corresponding to up-regulation and green to down-regulation in PRRSV-positive piglets. As expected, piglets from the two different groups clustered together, although some incorrectly assigned piglets could be observed (as confirmed by the calculated sensitivities and specificities values, see text).

Principal component analysis (PCA) was performed on the profiles of the 47 discriminatory peaks identified during the discovery and confirmed during the validation phase to identify and quantify independent sources of variation observed in the data. PCA analysis showed that 58.2% (PCA1), 17.9% (PCA2), and 12.9% (PCA3) of the total variability within the data was accounted for the X, Y, and Z axes, respectively. These axes were used to plot the data (Figure 2) and they provide an overview of the variation between the individual samples and show how samples grouped. Figure 2A showed three-dimensionally that the PCA peak profiles of piglets positive to PRRSV differed from piglets negative to PRRSV and revealed a good separation among the profiles of the two different groups, especially considering the high heterogeneity of the samples included in the study, as reported in the MM section and in [ Additional file 1: Table S1 and Additional file 2: Table S2]. Furthermore, with the exception of few outliers, PCA1 combined with PCA2 also separated well the two piglet populations (Figure 2B).

Figure 2
figure 2

Principal component analysis (PCA) showing the effects of the 47 significant discriminatory peaks on piglets positive or negative to PRRSV infection. The figure shows a projection of the measured peak intensities profiles onto the plane spanned by the three principal components (PCAs) that are the axes along which the data vary the most, for the 35 PRRSV-positive (blue) and the 35 PRRSV-negative (red) piglets of the validation study. PCA1, PCA2, and PCA3 accounted for 58.2%, 17.9%, and 12.9% of the variability in the data, respectively. PCA analysis illustrates a 3-dimentional plot comparison of PCA1, PCA2 and PCA3 in the three axes (A), as well as 2-dimentional score plot comparisons between PCA1 and PCA2 (B).

Comparison with relevant protein peaks and immunity genes related to PRRSV infection in other studies

To provide an overview of the current literature and to try to correlate the discriminatory peaks identified in this study with relevant proteins, we summarized in Table 3 the molecular weights of several peaks that have been shown to be related to PRRSV infection.

Table 3 Comparison between relevant PPRSV-related and pig proteins identified in other studies and the discriminatory peaks found in this study

First of all, we summarized the available information on the PRRS viral proteins. The PRRSV genome is ca. 15 kb in size and consists of the 5' untranslated region (UTR), at least nine open reading frames (ORFs), and the 3' UTR followed by a polyadenylation tail. The expected and experimentally identified MWs for each viral protein from different studies are reported in Table 3, along with the MW of the closest discriminatory peak identified in the current study.

Interestingly, the MW of the viral proteins ORF2b, ORF4, and ORF7 were very similar (difference of MW ≤0.3 kDa) to up-regulated discriminatory peaks identified here (Table 3).

As next, we compared proteins related to PRRSV infection that were identified in additional studies (Table 3); interestingly, all the 9 peaks found by [28], and in particular the only up-regulated in PRRSV infected (corresponding to the Alpha 1 S (a1S)-subunit of porcine Haptoglobin), showed minimal MW differences (≤0.3 kDa) with up-regulated peaks identified in this study (Table 3).

Additional discriminatory peaks found in the current study were very similar (MW differences ≤0.3 kDa) to those identified in other PRRS-related proteomic studies (Table 3). They corresponded to the following proteins: Glyceraldehyde-3-phosphate dehydrogenase, Proteasome activator hPA28 subunit beta, S100 calcium binding protein A10, Galectin 1, and Gastric-associated differentially expressed protein YA61P [26]; Heat shock 27 kDa protein 1, Superoxide dismutase 2, Myoglobin, and Vacuolar protein sorting 29 [29]; Heat shock protein 27 kDa and Nucleoside diphosphate kinase A [30]; Heat shock 27 kDa protein 1, Galectin 1, and Ubiquitin [31].

Discussion

In the present work, we show that proteomic fingerprint profiling is useful in researches on PRRS immuno-pathogenesis and might also be a robust, large scale diagnostic tool for the assessment of the proportion of PRRSV-positive weaning piglets without clinical symptoms in a herd. Indeed, we confirmed that the high-throughput capacity of the SELDI-TOF MS technology allows the screening for disease biomarkers of hundred of samples in a relative short-time period and with minimal sample preparation (as previously also reported by [32]).

Our results indicate that from the 200 significant peaks found in the discovery phase, a total of 47 could be confirmed in the validation phase. These values are comparable with another study where similar experimental conditions were applied to ovine sera [19].

Our findings also show that the combination of 14 discriminatory peaks in F1 on CM10 at low focus mass provided the highest sensitivity of 77% and specificity of 73% to correctly assign the piglets to the PPRSV-positive or PRRSV-negative groups. These percentages are in line with recent studies in humans using the same technology [33, 34]. Also the PCA results showed a good separation of the piglets in the two groups under examination. This was reached even though the tested piglets had large variability and heterogeneity, as they were collected from several farms located in different regions, and underwent high environmental pressures, typical of the field conditions. This is mainly due to the careful choice of the serum samples, where we tried to minimize the environmental differences by using same experimental parameters (e.g. sample collection procedures, storage, handling) and by including a similar number of pigs from the same breed (Large White) and with very similar sex ratios and ages (at weaning).

In a preliminary work [20] we had successfully transferred the experimental conditions used in profiling experiments of human sera to pig sera. However, in that work, none of the potential biomarkers identified in the discovery phase could be validated in the subsequent validation phase, because of high samples heterogeneity and high content of serum (e.g. albumin) and contaminant proteins (e.g. hemoglobin), having a negative effects on the detection of significant biomarkers, particularly those corresponding to the medium-low abundant proteins. It has been reported that low abundant proteins constitute about 1% of the entire human serum proteome, with the remaining 99% being comprised of only 22 proteins [35]. As it was therefore necessary to reduce the level of abundant proteins, in this follow up study, particular relevance was given to the content of the contaminant protein hemoglobin. Only non-hemolytic samples with similar, low contents of hemoglobin were included in the study. Additionally, to further increase the likelihood to identify statistically significant discriminatory biomarkers, we introduced a fractioning step based on anion-exchange chromatography. In a similar study performed with MALDI-TOF [28], where serum samples were analyzed in the first weeks (2–16) of PRRSV infection (also verified by PCR), a significantly lower number of peaks were identified compared to the present work. While protein peaks with M/Z values of 4.165, 4.460, 5.560, 8.330, 8.825, 12.250/12.550, and 14.010 kDa were found in 94 serum samples from 59 pigs, only one peak (9.244 kDa), corresponding to the alpha 1 S (a1S)-subunit of porcine Haptoglobin (Hp), was differentially up-regulated in PRRSV infected pigs. Interestingly, all these peaks were very similar (MW difference ≤0.3 kDa) with discriminatory peaks identified here (details in Table 3). Furthermore, two peaks identified in this study (23.162 and 14.843 kDa) were similar to peaks identified elsewhere (corresponding to Heat shock 27 kDa protein 1 [2931] and Galectin 1 [26, 31], respectively). In accordance with [31], the identified peak corresponding to Heat shock 27 kDa protein 1 was up-regulated, while the peak corresponding to Galectin 1 was down-regulated. Thus, these proteins seem to be very interesting and suitable candidates for future investigations.

The preponderance of the significant biomarkers had a molecular mass lower than 20 kDa, confirming that small peptides are a rich source of relevant biomarkers in SELDI-TOF MS analyses as previously reported in human [36] and ovine [19] sera. This may also partly be caused by the fact that the low molecular weight region (LMW) of the serum proteome, called peptidome, is an assortment of small intact proteins and proteolytic fragments of larger proteins, including several classes of physiologically important proteins like peptide hormones and components of both the innate and adaptive immune systems (i.e. cytokines and chemokines) [35, 37]. This is particularly interesting as the patho-physiological state of the body’s tissue is predominantly reflected in the LMW and low abundance region of the serum proteome, and specific protein fragments of the serum peptidome have been shown to contain a rich source of disease-specific diagnostic information and they have been correlated with disease stages in several studies (reviewed by [37]).

In agreement with other studies [29, 31], we found that the majority of the discriminatory biomarkers were up-regulated in PRRSV-positive piglets. This seems to suggest that the corresponding proteins might be of viral origin or related to the innate or adaptive immune responses (e.g. cytokines, chemokines, acute phase proteins, toll like receptors). In fact, several peaks showed high similarities (MW differences ≤0.3 kDa) with previous works, in particular regarding viral proteins (Table 3). The assignment of the discriminatory peak to a specific protein will require additional work, because the SELDI-TOF technology can only detect masses/peaks of proteins that are differentially expressed between samples but can not directly identify the proteins. This represents one of the major drawbacks of this technology compared to other methods. However, an advantage of the SELDI-TOF MS in this regard is that the results of this technique might lead to the identification of new proteins that were previously not correlated to the disease, and this might hopefully lead to the identification of new biomarkers representing the field situation. The interpretation of these results and the continuation of this project will benefit from the very imminent termination and publication of the sequence of the swine genome [38], which will definitely contribute to a more precise annotation and a better identification of genes and proteins and thus will greatly facilitate genome wide mapping association studies.

Conclusions

Although a combination of peaks identified with different experimental conditions (e.g. using different fractions and different surfaces) might have provided higher discriminatory power, here we developed a PRRSV diagnostic test based on peaks identified with the same experimental conditions (e.g. fraction, surface, and focus mass), which can be reproduced at high-throughput at reasonable costs. These results provide a set of proteomic biomarkers and related, optimized experimental conditions for high-throughput profiling of pig populations by SELDI-TOF MS for whole genome association studies, where identification of proteins underlying the phenotype can be made a posteriori. SELDI-TOF MS might therefore represent a complementary test or a possible alternative to classical (PCR) and more recent diagnostic methods (e.g. antibody detection in saliva) for profiling large flocks of pigs at reasonable costs, using blood samples that are routinely collected for general veterinary inspections. As well, these SELDI-TOF MS based tests could complement and provide a broader reference for emerging diagnostic methods and have potential applications for the detection of relevant proteins having highly heritable traits (e.g. acute phase proteins).

Methods

Piglets

A total of 120 serum samples of Large White piglets were selected from a well defined and characterized repository database, presently containing more than 20,000 swine samples from 18 different farms of the Lombardy region, Italy. Selection of the piglets aimed to minimize environmental factors and experimental conditions that might influence the results [39]. Hence, all piglets were from the same breed (Large White), had similar ages (weaning: 45–50 days), and their sera showed a low and comparable amount of hemoglobin (calculated as shown below).

In the discovery phase of the study, a total of 50 pig sera, 25 from PRRSV-positive and 25 from PRRSV-negative piglets, as determined by PCR (see below), were analyzed [ Additional file 1: Table S1]. The validation phase was performed with the same experimental conditions as the discovery phase. A total of 35 new PRRSV-positive and 35 new PRRSV-negative piglets were examined [ Additional file 2: Table S2]. The actual duration of infection for each individual PRRSV-positive piglet was unknown, as sera were collected and analyzed once for each piglet (at weaning: 45–50 days of age). None of the piglets was treated, as they did not show any symptom of the disease.

To ensure large variability and heterogeneity of the samples and minimize environmental differences, we included in the PRRSV-positive and -negative groups similar numbers of piglets with the same sex that originated from several farms located in different regions. In fact, PRRSV-positive piglets originated from 6 farms of the Lodi region (n = 8) and 7 farms of the Mantua region (n = 52), while PRRSV-negative piglets were collected in 5 farms around Lodi (n = 19) and 9 farms around Mantua (n = 41). Sex ratios males/females (44/76) were very similar in PRRSV-positive (21 vs. 39) and -negative (23 vs. 37) piglets, respectively.

Veterinary inspections of the overall clinical status of the piglets at the day of serum collection did not evidence any clinical symptoms of PRRS, including respiratory distress or sneezing.

Serum samples

All the serum samples were collected, stored, and handled in the same way. They were obtained for each piglet by storing two mL of whole blood without anticoagulants at room temperature (RT) for 4 h followed by centrifugation at 3,500 rpm for 4 min. As suggested in a previous work [20], an abundant quantity of hemoglobin in the serum can hide early diagnostic biomarkers of PRRSV by competing with the other serum components for the binding site of the chromatographic surfaces. To avoid the consequent signal suppression of the medium-low abundant proteins, only non-hemolytic samples were included in the present study.

A total of 200 clear, transparent sera without red pigmentation (low hemoglobin content) were first selected by visual screening from the total sera available in the database. Hemoglobin content of each serum sample was then determined according to [40] with minor modifications. A calibration curve was generated using five standard solutions (concentrations: 1.8, 3.6, 5.4, 7.2, and 9 μg/ml) of porcine hemoglobin diluted in 400 μL commercially available porcine serum (Sigma Aldrich, St Louis, MO, USA). Triplicate samples were incubated for 5·min at RT, then absorbance (E) was measured at 380, 415, and 440 nm. Absorbance at 380 and 440 nm was used to discern background absorbance flanking the absorbance peak (415·nm) of oxygenated hemoglobin. Absorbance due to hemoglobin was calculated as: E415–[(E380 + E440)/2]. Hemoglobin absorbance values of the samples were converted to μg/mL of hemoglobin by means of the calibration curve. Of the 200 initial samples, a total of 120 samples having an absorbance ≤ 0.085 (corresponding to a hemoglobin content below 4.52 μg/mL) were included in the study; 50 in the discovery and 70 in the validation phases, respectively.

Viral RNA extraction from the sera was performed following standard Roche procedures (High Pure Viral RNA Kit, Roche Diagnostics GmbH, Germany). Presence or absence of PRRSV was determined by multiplex PCR of conserved regions of viral ORF7 using primers and conditions previously described [41, 42]. The test also enabled to discriminate European and American genotypes and could detect all the different viral strains present in the Lombardy region at the time of sample collection.

Serum fractionation

All the detailed steps of the SELDI-TOF MS process performed here are schematically represented [see Additional file 3: Figure S1]. The protocol follows the manufacturer’s instructions with minor modifications (Bio-Rad Laboratories, ProteinChip® Serum Fractionation Kit manual).

Briefly, serum samples were pre-fractionated with U9 buffer (9 M urea, 2% 3-[(3-Cholamidopropyl)dimethylammonio]-1-propanesulfonate (CHAPS), 50 mM Tris–HCl, pH = 9) to favor dissociation of protein complexes [ Additional file 3: Figure S1A].

Sera were fractionated using a ProteinChip Q strong anion-exchange resin filtration plate (Bio-Rad Laboratories, Hercules, CA). The filtration plate was re-hydrated and equilibrated with rehydration buffer (50 mM Tris–HCl, pH = 9) and the resin washed with rehydration buffer and U1 solution (1 M urea, 0.2% CHAPS, 50 mM Tris–HCl, pH = 9) [ Additional file 3: Figure S1B]. Serum samples were then mixed with U1 solution and added to the equilibrated filtration plate. Successive elutions with different buffers with decreasing pH and a final organic solvent (= different fractions) were collected by centrifugation. The buffers used included pH = 9 (50 mM Tris–HCl, 0.1% n-octyl β-D-glucopyranoside (OGP)), pH = 7 (50 mM 4-(2-Hydroxyethyl)piperazine-1-ethanesulfonic acid (HEPES), 0.1% OGP), pH = 5 (100 mM Na acetate, 0.1% OGP), pH = 4 (100 mM Na acetate, 0.1% OGP), pH = 3 (100 mM Na citrate, 0.1% OGP), and organic solvent (33.3% isopropanol, 16.7% acetonitrile, 0.1% trifluoroacetic acid) [ Additional file 3: Figure S1C].

ProteinChip arrays

The six pH fractions obtained (F1 = pH9, F2 = pH7, F3 = pH5, F4 = pH4, F5 = pH3, and F6 = organic solvent) were profiled on weak cation-exchange (CM10), immobilized metal affinity capture-copper (IMAC30-CU), and reverse-phase (H50) ProteinChip® arrays. The arrays were initially placed in a Bioprocessor (C50-30011, Bio-Rad Laboratories) and then treated according to their surface [ Additional file 3: Figure S1D]. Each sample fraction was then bound/spotted randomly to the different ProteinChip® arrays using array-specific binding buffers [ Additional file 3: Figure S1E]. A 50% saturated sinapinic acid (SPA) matrix solution was finally added to each spot on the ProteinChip array prior to the final analysis [ Additional file 3: Figure S1F].

SELDI-TOF MS analysis

ProteinChip arrays were read using a Ciphergen Protein-Chip Reader PCS4000 model and data were analyzed with Ciphergen Express Software (Ciphergen Biosystems).

Profiles were collected in the range 1–200 kDa at the two different ion focus mass 10 kDa (“low focus mass”) and 50 kDa (“high focus mass”). The instrument was calibrated for dataset collection using all-in-one peptide standard (Bio-Rad Laboratories) in the 1–20 kDa range for 10 kDa low ion focus mass and all-in-one protein standard in the 20–200 kDa range for 50 kDa high ion focus mass [ Additional file 3: Figure S1G].

Ciphergen Express software analysis

Spectra were normalized by total ion current, starting and ending at the M/Z of the collection ranges (1–20 or 20–200 kDa) after baseline subtraction and noise calculation. Outlier spectra were removed. The spectra were aligned to a reference spectrum with the normalization factor nearest 1.0. The spectra were aligned only if the percentage coefficient of variation was reduced after the alignment. Peaks from the different spectra were aligned using the cluster wizard function of the Ciphergen Express 3.0.6 software. The peak detection was automated within the M/Z range of analysis. Peaks were detected on the first pass when the signal-to-noise (S/N) ratio was 7 and the peak was 5 times the valley depth. Peaks below threshold were deleted and all first-pass peaks were preserved. Clusters were created within 0.15% of M/Z for each peak detected in the first pass. The clusters were completed by adding peaks with S/N ratio of 2 and two times the valley depth. P-values and ROC/AUC (Receiver Operating Characteristic/ Area Under Curve) values were calculated by using the P-value wizard.

A 2-tailed t-test was used for statistical analysis of differences in peak intensity between groups. P-values below 0.05 were considered statistically significant. Principal component analysis (PCA) and agglomerative hierarchical clustering algorithm were applied to investigate the pattern among the different statistically significant peaks.

PCA is a multivariate data analysis that transforms without a loss of essential information a number of correlated variables into a smaller number of uncorrelated variables called principal components (PCs), which can explain sufficiently the data structure. PCA transformation allows studying many variables simultaneously, showing how similar samples are correlated and grouped together. The data structure is visualized directly in a graphical way by projection of objects onto the space defined by the selected PCAs (for details see [43]).

Finally, to evaluate the influence of external variables (e.g. sample processing and acquisition) on the system under study and to calculate the dispersion of the acquired data, the coefficient of variation (CV), which is the normalized measure of dispersion of a probability distribution and shows the% dispersion of the data in rapport to the media (intensity variation), was also calculated. Six serum samples commercially available were prepared and analyzed in parallel with the pig samples of both, discovery and validation phases. The CV was calculated for all fractions and surfaces by choosing 6 peaks evenly distributed along the entire range.

References

  1. Davies G, Genini S, Bishop SC, Giuffra E: An assessment of opportunities to dissect host genetic variation in resistance to infectious diseases in livestock. Animal 2009, 3: 415–436. 10.1017/S1751731108003522

    Article  CAS  PubMed  Google Scholar 

  2. An TQ, Tian ZJ, Xiao Y, Li R, Peng JM, Wei TC, Zhang Y, Zhou YJ, Tong GZ: Origin of highly pathogenic porcine reproductive and respiratory syndrome virus, China. Emerg Infect Dis 2010, 16: 365–367. 10.3201/eid1602.090005

    Article  PubMed Central  PubMed  Google Scholar 

  3. Zhou L, Yang H: Porcine reproductive and respiratory syndrome in China. Virus Res 2010, 154: 31–37. 10.1016/j.virusres.2010.07.016

    Article  CAS  PubMed  Google Scholar 

  4. Lunney JK, Benfield DA, Rowland RR: Porcine reproductive and respiratory syndrome virus: an update on an emerging and re-emerging viral disease of swine. Virus Res 2010, 154: 1–6. 10.1016/j.virusres.2010.10.009

    Article  CAS  PubMed  Google Scholar 

  5. Snijder EJ, Meulenberg JJ: The molecular biology of arteriviruses. J Gen Virol 1998, 79: 961–979.

    CAS  PubMed  Google Scholar 

  6. Cho JG, Dee SA: Porcine reproductive and respiratory syndrome virus. Theriogenology 2006, 66: 655–662. 10.1016/j.theriogenology.2006.04.024

    Article  PubMed  Google Scholar 

  7. Neumann EJ, Kliebenstein JB, Johnson CD, Mabry JW, Bush EJ, Seitzinger AH, Green AL, Zimmerman JJ: Assessment of the economic impact of porcine reproductive and respiratory syndrome on swine production in the United States. J Am Vet Med Assoc 2005, 227: 385–392. 10.2460/javma.2005.227.385

    Article  PubMed  Google Scholar 

  8. Horter DC, Pogranichniy RM, Chang CC, Evans RB, Yoon KJ, Zimmerman JJ: Characterization of the carrier state in porcine reproductive and respiratory syndrome virus infection. Vet Microbiol 2002, 86: 213–228. 10.1016/S0378-1135(02)00013-5

    Article  PubMed  Google Scholar 

  9. Klinge KL, Vaughn EM, Roof MB, Bautista EM, Murtaugh MP: Age-dependent resistance to Porcine reproductive and respiratory syndrome virus replication in swine. Virol J 2009, 6: 177. 10.1186/1743-422X-6-177

    Article  PubMed Central  PubMed  Google Scholar 

  10. Kimman TG, Cornelissen LA, Moormann RJ, Rebel JM, Stockhofe-Zurwieden N: Challenges for porcine reproductive and respiratory syndrome virus (PRRSV) vaccinology. Vaccine 2009, 27: 3704–3718. 10.1016/j.vaccine.2009.04.022

    Article  CAS  PubMed  Google Scholar 

  11. Genini S, Delputte PL, Malinverni R, Cecere M, Stella A, Nauwynck HJ, Giuffra E: Genome-wide transcriptional response of primary alveolar macrophages following infection with porcine reproductive and respiratory syndrome virus. J Gen Virol 2008, 89: 2550–2564. 10.1099/vir.0.2008/003244-0

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  12. Wills RW, Doster AR, Galeota JA, Sur JH, Osorio FA: Duration of infection and proportion of pigs persistently infected with porcine reproductive and respiratory syndrome virus. J Clin Microbiol 2003, 41: 58–62. 10.1128/JCM.41.1.58-62.2003

    Article  PubMed Central  PubMed  Google Scholar 

  13. Somasundaram K, Nijaguna MB, Kumar DM: Serum proteomics of glioma: methods and applications. Expert Rev Mol Diagn 2009, 9: 695–707. 10.1586/erm.09.52

    Article  CAS  PubMed  Google Scholar 

  14. Hodgetts A, Levin M, Kroll JS, Langford PR: Biomarker discovery in infectious diseases using SELDI. Future Microbiol 2007, 2: 35–49. 10.2217/17460913.2.1.35

    Article  CAS  PubMed  Google Scholar 

  15. Devitt EJ, Power KA, Lawless MW, Browne JA, Gaora PO, Gallagher WM, Crowe J: Early proteomic analysis may allow noninvasive identification of hepatitis C response to treatment with pegylated interferon alpha-2b and ribavirin. Eur J Gastroenterol Hepatol 2011, 23: 177–183. 10.1097/MEG.0b013e3283424e3e

    Article  CAS  PubMed  Google Scholar 

  16. Toro-Nieves DM, Rodriguez Y, Plaud M, Ciborowski P, Duan F, Perez Laspiur J, Wojna V, Melendez LM: Proteomic analyses of monocyte-derived macrophages infected with human immunodeficiency virus type 1 primary isolates from Hispanic women with and without cognitive impairment. J Neurovirol 2009, 15: 36–50. 10.1080/13550280802385505

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  17. Pang RT, Poon TC, Chan KC, Lee NL, Chiu RW, Tong YK, Wong RM, Chim SS, Ngai SM, Sung JJ, Lo YM: Serum proteomic fingerprints of adult patients with severe acute respiratory syndrome. Clin Chem 2006, 52: 421–429. 10.1373/clinchem.2005.061689

    Article  CAS  PubMed  Google Scholar 

  18. Fenollar F, Goncalves A, Esterni B, Azza S, Habib G, Borg JP, Raoult D: A serum protein signature with high diagnostic value in bacterial endocarditis: results from a study based on surface-enhanced laser desorption/ionization time-of-flight mass spectrometry. J Infect Dis 2006, 194: 1356–1366. 10.1086/508429

    Article  CAS  PubMed  Google Scholar 

  19. Rioux MC, Carmona C, Acosta D, Ward B, Ndao M, Gibbs BF, Bennett HP, Spithill TW: Discovery and validation of serum biomarkers expressed over the first twelve weeks of Fasciola hepatica infection in sheep. Int J Parasitol 2008, 38: 123–136. 10.1016/j.ijpara.2007.07.017

    Article  CAS  PubMed  Google Scholar 

  20. Genini S, Cantu M, Botti S, Malinverni R, Costa A, Marras D, Giuffra E: Diagnostic markers for diseases: SELDI-TOF profiling of pig sera for PRRS. Dev Biol 2008, 132: 399–403.

    CAS  Google Scholar 

  21. Zhong L, Taylor DL, Whittington RJ: Proteomic profiling of ovine serum by SELDI-TOF MS: optimisation, reproducibility and feasibility of biomarker discovery using routinely collected samples. Comp Immunol Microbiol Infect Dis 2010, 33: 47–63. 10.1016/j.cimid.2008.07.009

    Article  PubMed  Google Scholar 

  22. Leak LV, Liotta LA, Krutzsch H, Jones M, Fusaro VA, Ross SJ, Zhao Y, Petricoin EF: Proteomic analysis of lymph. Proteomics 2004, 4: 753–765. 10.1002/pmic.200300573

    Article  CAS  PubMed  Google Scholar 

  23. Nelsen CJ, Murtaugh MP, Faaberg KS: Porcine reproductive and respiratory syndrome virus comparison: divergent evolution on two continents. J Virol 1999, 73: 270–280.

    PubMed Central  CAS  PubMed  Google Scholar 

  24. Zimmerman J, Karriker L, Ramirez A, Schwartz K, Stevenson G: Diseases of Swine. New York: John Wiley & Sons Inc; 2012.

    Google Scholar 

  25. Meulenberg JJ, Petersen-den Besten A, De Kluyver EP, Moormann RJ, Schaaper WM, Wensvoort G: Characterization of proteins encoded by ORFs 2 to 7 of Lelystad virus. Virology 1995, 206: 155–163. 10.1016/S0042-6822(95)80030-1

    Article  CAS  PubMed  Google Scholar 

  26. Zhang C, Xue C, Li Y, Kong Q, Ren X, Li X, Shu D, Bi Y, Cao Y: Profiling of cellular proteins in porcine reproductive and respiratory syndrome virus virions by proteomics analysis. Virol J 2010, 7: 242. 10.1186/1743-422X-7-242

    Article  PubMed Central  PubMed  Google Scholar 

  27. Wu WH, Fang Y, Farwell R, Steffen-Bien M, Rowland RR, Christopher-Hennings J, Nelson EA: A 10-kDa structural protein of porcine reproductive and respiratory syndrome virus encoded by ORF2b. Virology 2001, 287: 183–191. 10.1006/viro.2001.1034

    Article  CAS  PubMed  Google Scholar 

  28. Gnanandarajah JS, Dvorak CM, Johnson CR, Murtaugh MP: Presence of free haptoglobin alpha 1 S-subunit in acute porcine reproductive and respiratory syndrome virus infection. J Gen Virol 2008, 89: 2746–2753. 10.1099/vir.0.2008/003954-0

    Article  CAS  PubMed  Google Scholar 

  29. Zhang H, Guo X, Ge X, Chen Y, Sun Q, Yang H: Changes in the cellular proteins of pulmonary alveolar macrophage infected with porcine reproductive and respiratory syndrome virus by proteomics analysis. J Proteome Res 2009, 8: 3091–3097. 10.1021/pr900002f

    Article  CAS  PubMed  Google Scholar 

  30. Sagong M, Lee C: Differential cellular protein expression in continuous porcine alveolar macrophages regulated by the porcine reproductive and respiratory syndrome virus nucleocapsid protein. Virus Res 2010, 151: 88–96. 10.1016/j.virusres.2010.04.003

    Article  CAS  PubMed  Google Scholar 

  31. Ding Z, Li ZJ, Zhang XD, Li YG, Liu CJ, Zhang YP, Li Y: Proteomic alteration of Marc-145 cells and PAMs after infection by porcine reproductive and respiratory syndrome virus. Vet Immunol Immunopathol 2012, 145: 206–213. 10.1016/j.vetimm.2011.11.005

    Article  CAS  PubMed  Google Scholar 

  32. Issaq HJ, Veenstra TD, Conrads TP, Felschow D: The SELDI-TOF MS approach to proteomics: protein profiling and biomarker identification. Biochem Biophys Res Commun 2002, 292: 587–592. 10.1006/bbrc.2002.6678

    Article  CAS  PubMed  Google Scholar 

  33. Monari E, Casali C, Cuoghi A, Nesci J, Bellei E, Bergamini S, Fantoni LI, Natali P, Morandi U, Tomasi A: Enriched sera protein profiling for detection of non-small cell lung cancer biomarkers. Proteome Sci 2011, 9: 55. 10.1186/1477-5956-9-55

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  34. Hu Q, Huang Y, Wang Z, Tao H, Liu J, Yan L, Wang K: Application of Surface-Enhanced Laser Desorption/Ionization Time-of-Flight Mass Spectrometry Coupled With an Artificial Neural Network Model for the Diagnosis of Hepatocellular Carcinoma. Hepatogastroenterology 2012, 59: 1902–1906.

    Article  CAS  PubMed  Google Scholar 

  35. Tirumalai RS, Chan KC, Prieto DA, Issaq HJ, Conrads TP, Veenstra TD: Characterization of the low molecular weight human serum proteome. Mol Cell Proteomics 2003, 2: 1096–1103. 10.1074/mcp.M300031-MCP200

    Article  CAS  PubMed  Google Scholar 

  36. Hortin GL: The MALDI-TOF mass spectrometric view of the plasma proteome and peptidome. Clin Chem 2006, 52: 1223–1237. 10.1373/clinchem.2006.069252

    Article  CAS  PubMed  Google Scholar 

  37. Petricoin EF, Belluco C, Araujo RP, Liotta LA: The blood peptidome: a higher dimension of information content for cancer biomarker discovery. Nat Rev Cancer 2006, 6: 961–967. 10.1038/nrc2011

    Article  CAS  PubMed  Google Scholar 

  38. Archibald AL, Bolund L, Churcher C, Fredholm M, Groenen MA, Harlizius B, Lee KT, Milan D, Rogers J, Rothschild MF, Uenishi H, Wang J, Schook LB: Swine Genome Sequencing Consortium: Pig genome sequence–analysis and publication strategy. BMC Genomics 2010, 11: 438. 10.1186/1471-2164-11-438

    Article  PubMed Central  PubMed  Google Scholar 

  39. Biffani S, Botti S, Bishop SC, Stella A, Giuffra E: Using SNP array data to test for host genetic and breed effects on Porcine Reproductive and Respiratory Syndrome Viremia. BMC Proc 2011,5(Suppl 4):S28. 10.1186/1753-6561-5-S4-S28

    Article  PubMed Central  PubMed  Google Scholar 

  40. Rider CV, Gorr TA, Olmstead AW, Wasilak BA, LeBlanc GA: Stress signaling: coregulation of hemoglobin and male sex determination through a terpenoid signaling pathway in a crustacean. J Exp Biol 2005, 208: 15–23. 10.1242/jeb.01343

    Article  CAS  PubMed  Google Scholar 

  41. Suarez P, Zardoya R, Prieto C, Solana A, Tabares E, Bautista JM, Castro JM: Direct detection of the porcine reproductive and respiratory syndrome (PRRS) virus by reverse polymerase chain reaction (RT-PCR). Arch Virol 1994, 135: 89–99. 10.1007/BF01309767

    Article  CAS  PubMed  Google Scholar 

  42. Kono Y, Kanno T, Shimizu M, Yamada S, Ohashi S, Nakamine M, Shirai J: Nested PCR for detection and typing of porcine reproductive and respiratory syndrome (PRRS) virus in pigs. J Vet Med Sci 1996, 58: 941–946. 10.1292/jvms.58.10_941

    Article  CAS  PubMed  Google Scholar 

  43. Beebe KR, Pell RJ, Seasholtz MB: Chemometrics: a Practical Guide. New York: John Wiley & Sons Inc; 1998.

    Google Scholar 

Download references

Acknowledgements

The authors are very thankful to Maria Cecere, Roberto Grande, Roberto Malinverni, Silvia Rossini, and Elisa Filippi for technical assistance in the laboratory and for helpful comments. They also are very thankful to Joan K. Lunney for critical reading and revision of the manuscript. This project was supported by a grant of the Italian Ministry of Research (MIUR project, art.10 D.M. 593/00).

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Sem Genini.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

SG established the experimental settings, prepared the samples, helped conceiving the study design, and wrote the manuscript. TP helped with the analysis of the data and the manuscript writing. ACO conducted the SELDI-TOF MS experiment and performed statistical analyses. SB arranged for sample collection and helped to draft the manuscript. MVL developed the primers to detect PRRSV and performed the PCR analyses to assess presence/absence of the virus. ACA was responsible for data storage and database maintenance. EG coordinated the overall project and participated in the design of the study and helped to draft the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material

12953_2012_370_MOESM1_ESM.doc

Additional file 1: Table S1. Pigs tested with SELDI-TOF MS during the discovery phase of the study. List of the 25 positive and 25 negative pigs to PRRS (PCR-tested) analyzed with SELDI-TOF MS during the discovery phase of the study. The pig ID is reported with the total absorbance and the total amount of hemoglobin present in the sample, the status regarding the PRRS virus, as well as the sex and the number and location of the farm (MA = Mantua region, LO = Lodi region). (DOC 78 KB)

12953_2012_370_MOESM2_ESM.doc

Additional file 2: Table S2. Pigs tested with SELDI-TOF MS during the validation phase of the study. List of the 35 positive and 35 negative pigs to PRRS (PCR-tested) analyzed with SELDI-TOF MS during the validation phase of the study. The pig ID is reported with the total absorbance and the total amount of hemoglobin present in the sample, the status regarding the PRRS virus, as well as the sex and the number and location of the farm (MA = Mantua region, LO = Lodi region). (DOC 87 KB)

12953_2012_370_MOESM3_ESM.doc

Additional file 3: Figure S1. Detailed protocol of the SELDI-TOF MS analysis. Schematic illustration of the protocol used for SELDI-TOF MS, which follows the manufacturer’s instruction manual with minor modifications (Bio-Rad Laboratories, ProteinChip® Serum Fractionation Kit manual). The protocol is divided in 7 main steps: A) Pre-fractionation of the sera; B) Rehydration and equilibration of the Protein Chip Q strong anion-exchange resin filtration plate; C) Fractionation of the sera; D) Preparation of ProteinChip® Arrays; E) Binding of the serum fractions to the arrays; F) Preparation and application of the matrix; and G) SELDI-TOF MS analysis. (DOC 270 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Genini, S., Paternoster, T., Costa, A. et al. Identification of serum proteomic biomarkers for early porcine reproductive and respiratory syndrome (PRRS) infection. Proteome Sci 10, 48 (2012). https://doi.org/10.1186/1477-5956-10-48

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1477-5956-10-48

Keywords