Abstract

The arenaviruses include a number of important pathogens including Lassa virus and Junin virus. Presently, the only treatment is supportive care and the antiviral Ribavirin. In the event of an epidemic, patient triage may be required to more effectively manage resources; the development of prognostic biomarker signatures, correlating with disease severity, would allow rational triage. Using a pair of arenaviruses, which cause mild or severe disease, we analyzed extracts from infected cells using SELDI mass spectrometry to characterize potential biomarker profiles. EDGE analysis was used to analyze longitudinal expression differences. Extracts from infected guinea pigs revealed protein peaks which could discriminate between mild or severe infection and between times post-infection. Tandem mass-spectrometry identified several peaks, including the transcriptional regulator prothymosin- . Further investigation revealed differences in secretion of this peptide. These data show proof of concept that proteomic profiling of host markers could be used as prognostic markers of infectious disease.

1. Introduction

Viral hemorrhagic fevers (VHFs) represent an important class of emerging diseases and are classified as category A priority pathogens by the National Institute for Allergy and Infectious Diseases. The VHFs Lassa fever and the Argentine, Bolivian, Venezuelan, and Brazilian hemorrhagic fevers are caused by members of the Arenaviridae [1]. Despite their status as emerging infectious diseases, little is known about the molecular basis of the pathology of these diseases, and this has hindered the design of antiviral therapeutics. Arenavirus pathogenesis is believed to involve the dysregulation of cytokines; for example, infection with the virulent Lassa virus (LASV) prevents macrophage activation and TNF- production, whereas infection with the avirulent Mopeia virus induces proinflammatory cytokine production [2], a similar effect has been shown with attenuated and virulent isolates of Pichindé virus (PICV) [3]; proteomic and kinomic level analyses have also shown differential cell signaling events induced by attenuated and virulent arenavirus infection [46].

Currently, the only drug used in the treatment of these infections is the broad-spectrum antiviral, Ribavirin, which requires administration within the first week of infection to be most efficacious [7]. As only a small proportion of patients infected with LASV develop a hemorrhagic fever, indicators of prognosis will be critical in the efficient management of healthcare resources during epidemics. Currently, the best prognostic marker for severe disease is viral load, with high levels of virus correlating with increased likelihood of severe disease [8]. However, quicker methods of prognosis determination would be of significant use in a clinical setting.

We have previously used surface-enhanced laser desorption ionization time-of-flight (SELDI-TOF) mass spectrometry to characterize proteins bound to a thioaptamer (oligonucleotides with thiophosphate backbone substitutions) on a ProteinChip array [9]. SELDI is a high-throughput profiling platform that combines retentive chromatography and mass spectrometry in one platform REFS. The process involves selective extraction and retention of proteins on chromatographic chip surfaces arranged in an array of spots and their subsequent analysis by mass spectrometry.

In this study, we have used SELDI-TOF to compare cell extracts for representative protein “signatures” following infection. PICV is a genetically characterized model for Lassa fever and recent construction infectious clones may allow the role of individual sequence changes in determining pathogenesis to be elucidated [10, 11]. Using two variants of PICV, P2, which causes a mild, self-limiting infection, and P18, which causes a lethal hemorrhagic disease [12], we can identify signatures unique to mild and severe infection, and which reflect common responses.

In this report, we describe the host-cell protein “signatures” following mock, P2, and P18 PICV infection over time. We show that, consistent with previous data from several groups, the virulent variant more closely resembles mock-infection, suggesting a suppression of cellular responses to infection [2, 5, 13]. In addition, we identify one of these protein peaks by tandem mass spectrometry as prothymosin- , which may indicate a potential therapeutic target for the treatment of hemorrhagic arenavirus infection.

2. Materials and Methods

2.1. Cell Culture and Viruses

Murine monocyte-like P388D1 cells were maintained in RPMI supplemented with 5% fetal bovine serum and 5 mM glutamine. Cells were cultured in the absence of antibiotics to ensure no underlying contamination. Cells were infected in triplicate with PEG purified P2 or P18 PICV at a multiplicity of infection of 5 or mock infected with virus purification medium. Cells were harvested at 8 timepoints and cytoplasmic and nuclear extracts prepared. For animal experiments, male Hartley guinea pigs were infected intraperitoneally with 1000 plaque forming units of the P2 or P18 variant of PICV or mock-infected with phosphate buffered saline ( per group/timepoint) and sacrificed at 1 and 6 days post infection. Animal experiments were performed following institutional animal care and use approved guidelines and protocols.

2.2. Cytoplasmic and Nuclear Fractionation of Cells

Cells were fractionated into cytoplasmic and nuclear fractions as previously described [14] with the addition of a nuclear purification step using Optiprep (Axis-Shield, Oslo, Norway) gradients. Briefly, lysates were underlayered with 10 mL 30% Optiprep and 5 mL 35% Optiprep and centrifuged at C for 30 minutes at 4300 g. The interface was removed and placed in a fresh tube which was filled with sucrose buffer I (described in Dyer & Herzog, 1995) plus 1.5 mM CaCl2. Following centrifugation at C for 15 minutes at 1900 g, the pellet was resuspended in sucrose buffer I and the centrifugation repeated. Nuclear lysis was completed following the referenced protocol.

2.3. Fractionation of Cytoplasmic Portion of P388D1 Cells

5.4 mL buffer containing 10 mM Tris and 9 M urea (pH 9) was added to 2.7 mL combined cytoplasmic portion of mock-, P2-, and P18-infected macrophages, and the mixture was incubated at C for 20 minutes. It then was dialyzed against 50 mM Tris (pH 8), after which the dialyzed solution was added to 30  L packed Q HyperD F beads (Pall corporation, PA), and incubated for 30 minutes at C. The beads were collected by centrifugation (3000 g, 1 minutes) and the supernatant was saved as flow-through. The proteins bound to the beads were eluted by lowering the pH of the buffer in five fractions (90  L of each eluent in 50 mM of each pH buffer): Q1 (pH 6.5), Q2 (pH 5.5), Q3 (pH 4), Q4 (pH 3) and Q5 (pH  M NaCl). 10  L of each fraction was added to 90  L 50 mM Tris (pH 8) buffer and was captured on Q10 anion exchange ProteinChip array surfaces (Bio-Rad Laboratories, Hercules, CA) with a bioprocessor Biomek 3000 (Beckman Coulter, Fullerton, CA).

2.4. SELDI-TOF Mass Spectrometry and Spectral Analysis

The mass spectra of the proteins captured on Q10 chips were recorded on the PCS 4000 ProteinChip array reader (Bio-Rad laboratories, Hercules, CA). Experiments were run in triplicate at the 8 time points mentioned above, leading to a total of 72 samples each for the cytoplasmic and nuclear fractions. SELDI spectra were recorded on the PCS4000 ProteinChip System at a laser intensity setting of 2600. Spectra were analyzed with the Ciphergen Express (CE) software. Analysis was limited to the molecular weight (MW) range 2 kDa–20 kDa. Peak intensities were scaled using CE.

2.5. Statistical Analysis

Log-transformed peak intensity values were subjected to a 2-way ANOVA test with Matlab v. R2007a (The Mathworks, Natick, MA). -values were calculated for three null hypotheses: (a) that all samples from the Mock, P2, and P18 infections came from the same population, (b) that all samples from the eight different time points came from the same population, and (c) that the effects due to the two factors, infection and time, are additive.

Spotfire DecisionSite 9.0 (TIBCO Software Inc., Palo Alto, CA) was used for hierarchical clustering (HC) to group the biomarkers detected by the timecourse analysis described above. SELDI peak intensities were first Z-score normalized and HC was then run using the unweighted average method, with Euclidean distance as the similarity measure and average value as the ordering function. Following HC, Principal Components Analysis (PCA) was performed transforming the data into a small number of dimensions to visualize the clustering of the biomarkers.

We also analyzed data using “Extraction of Differential Genomic Expression” (EDGE) as this technique is designed for the analysis of timecourse data [15, 16]. The method builds on the false-discovery rate (FDR) developed by Benjamini and Hochberg to control the expected proportion of falsely rejected null hypotheses [17]. For the multiple significance testing performed in the analysis of-omics data, the -values provide a more meaningful and practical measure than the -value [18]. Appropriate cutoffs are imposed on the -value to control the desired proportion of falsely rejected hypotheses.

2.6. Protein Identification and Verification

The cytoplasmic samples were treated with 9 M urea to denature the proteins, fractionated by Q HyperD F beads, and eluted into flow-through (FT), and five pH factions at pH 6.5, 5.5, 4.0, 3.0, and pH  M NaCl. The fractions were captured on Q10 chips. The chip was loaded on the high-resolution QSTAR XL spectrophotometer and data collection was performed. The fragment masses were used as probes to search the mouse protein database for protein identification.

2.7. Biomarker Identification Using SELDI Quadrupole Time of Flight Tandem Mass Spectrometry

Tandem mass spectrometric peptide sequencing was accomplished using QSTAR XL (Applied Biosystems Inc., Foster City, CA) quadrupole time of flight instrument equipped with Ciphergen ProteinChip Interface PCI 1000. The instrument was calibrated externally using the known MS/MS spectrum of ACTH peptide (1–24) at 2465.2 m/z, where four fragment ions and the parent ion were used as calibration points. All mass spectra were acquired in positive ion mode with collision energy that follows the rule of approximate 50 eV/kDa parent peptide mass. Raw data were analyzed using the instrument’s software. Fragmentation information from tandem MS spectrum was first submitted to the sequence query program in Matrix Science (http://www.matrixscience.com/) for possible identities with a precision tolerance of 0.4 Da for MS/MS fragments and 70 ppm for the parent peptide.

2.8. Verification of the Identified Biomarkers

60 uL of packed protein G agarose beads (Santa Cruz Biotechnology, CA) were incubated with 30 uL mouse monoclonal antibody against prothymosin alpha (Alexis Biochemicals, CA) for 1 hour at 4C, followed by three times wash of 1XPBS, 50 uL combined cytoplasmic portion of mock-, P2-, and P18-infected macrophages was then added to the beads, and the solution was incubated at 4C for 2 hours, followed by three times wash of 1XPBS + 0.1% Triton X-100. After a brief wash with water, proteins bound to the antibody against prothymosin alpha were eluted with 60 uL of elution buffer containing 0.3% TFA and 50% acetonitrile. The eluant was analyzed by PCS 4000 ProteinChip array reader (Bio-Rad laboratories, Hercules, CA).

2.9. Enzyme-Linked Immunosorbant Assay

Supernatants from infected triplicate cultures were added to 96-well plates in triplicate wells (Nunc) and incubated overnight at C. Wells were then blocked with 5% w/v bovine serum albumin for 1 hour. Primary antibody (antiprothymosin- 1, Santa Cruz, Santa Cruz, CA, 1/2000) was added to wells in blocking solution and incubated for 2 hours. Wells were then washed four times in phosphate buffered saline + 0.05% tween-20 for 5 minutes prior to incubation with secondary antibody (anti-rabbit Ig HRP conjugate, 1/2000) in blocking solution for 1 hour. Wells were washed as before and levels of [pro] thymosin are indirectly assayed by the addition of the colorimetric substrate TMB (Sigma, St Louis, MO).

3. Results

3.1. SELDI Mass Spectrometry of Guinea Pig Peritoneal Cells

Peritoneal cell extracts from infected guinea pigs were analyzed by surface enhanced laser desorption ionization (SELDI) mass spectrometry (Figure 1). Peaks were identified that could discriminate between type of infection (e.g., Figure 1, box A) and stage of infection (e.g., Figure 1, box B). Consistent with our previous findings, infection with the virulent P18 variant of the virus resulted in a spectrum that more closely resembled mock infection than infection with the attenuated P2 virus.

We next infected triplicate cultures of murine P388D1 macrophage-like cells for a timecourse of infection from 15 minutes to 16 hours. We used a cell line for this study in order to increase the sample quantity for downstream assays and database searching for identification by mass spectrometry. We have previously used this cell line to investigate PICV-host interactions and are familiar with how cellular responses in this cell line model those seen in primary guinea pig macrophages. Triplicate nuclear and cytoplasmic samples were prepared following mock-infection or infection with P2 or P18 PICV at 15 m, 30 m 1, 2, 4, 8, 12, and 16 hours post infection. Samples were analyzed by SELDI mass spectrometry. Statistical analysis by ANOVA and EDGE identified 27 and 35 significantly different peaks, respectively.

Hierarchical cluster analysis revealed similarities in peak patterns between P2 and P18 infection at late times (12 and 16 hours) post infection (data not shown), but that P18 was most similar to mock infection at early times post infection, consistent with our previous results. We next selected several significant peaks for identification by mass spectrometry. We identified several peaks as two SH3-domain binding proteins, stathmin, prothymosin- , thioredoxin, and various ribosomal proteins. Prothymosin- was identified as a significantly differentially expressed protein peak with and EDGE value of 0.0164. We have previously identified prothymosin- as a node in a signaling network constructed following kinomics analysis of PICV-infected cells [5] and so was selected for additional study.

3.2. Verification of the 3787 Peak as Prothymosin- from On-Chip Digestion and MS/MS

The pattern of prothymosin- peak expression in a representative sample of cytoplasmic extracts is shown in Figure 2. Lower peak intensity following infection with P18 PICV can be seen at 1, 2, and 4 hours post-infection. However, prothymosin- can shuttle between the cytoplasm and the nucleus where it functions as a histone H1 binding protein and is associated with cell growth [1921]. Prothymosin- and the peptide thymosin- 1, produced by proteolytic processing of prothymosin- and an immune-modulator, can also be secreted [2224]; for these reasons, the changes in peak intensities may not reveal the functional changes and effects in prothymosin- activity and regulation.

In order to verify the marker with molecular mass at 3787 Da, the cytoplasmic portion of the combined mock-, P2-, and P18-infected macrophages was fractionated on Q HyperD F beads with six fractions including the flow-through; each fraction was further captured on Q10 Protein Chip array surfaces. The marker protein was highly enriched in Q3 fraction on Q10 array surfaces (data not shown). The molecular mass of the marker is only 3787 Da; since the QSTAR XL instrument can analyze parent ion up to 5400 m/z, the peak could be identified as prothymosin- without any further purification (Figure 3). The MS/MS spectrum is shown in Figure 4(a). We confirmed the identification of prothymosin- using immunoprecipitation. Protein G agarose beads coupled to a monoclonal antibody against prothymosin alpha were used to immunoprecipitate prothymosin- from extracts; the same protein peak at 3787 Da was identified (Figure 4(b)).

3.3. Verification of Protein Identification and (Pro)Thymosin- Secretion

Prothymosin- contains the secreted peptide thymosin- 1 [2527]. We harvested supernatants from triplicate infected cultures over a 16-hour timecourse of infection and assayed the levels of (pro)thymosin- in the supernatants (Figure 4(a)). As the thymosin peptide is contained within prothymosin- , we are unable to distinguish between these peptides. Infection with the attenuated P2 variant of PICV induced consistently higher levels of (pro)thymosin- in supernatants than mock infection or infection with the virulent P18 variant, ranging a 2–2.5-fold increase over baseline (mock-infected) levels (Figure 5(a)). Levels of secreted (pro)thymosin- were significantly higher ( ) from 15 minutes to 8 hours postinfection. Consistent with many of the significant peaks in the SELDI analysis and many of our previous observations, the response induced by P18 infection more closely resembled mock-infection than infection with P2 PICV [35], although P18 did begin to induce higher levels of secreted (pro)thymosin- than mock-infection at late times postinfection.

We infected triplicate cultures of P388D1 cells with P2 or P18 PICV or mock infected with PEG purification medium and prepared whole cell extracts at 24 hours post infection. Stathmin expression was assayed by immunoblot. Stathmin expression was significantly reduced in both P2 and P18 infected cells compared to mock, with P18 infection showing reduced stathmin expression compared to P2 across triplicate experiments (Figure 5(b)). Stathmin expression is controlled by a number of transcription factors, including a signaling module based around NF-κB, AP-1, c-Myc, and p53, which we have previously shown to be differentially regulated in PICV infection [4, 5], and which may suggest the molecular basis for differential stathmin expression.

4. Discussion

Many virus infections present with nonspecific flu-like symptoms complicating diagnosis; additionally, very few markers exist for determining the prognosis of infection. In particular cases, such as in the response to an epidemic or an intentional release, the ability to quickly triage patients on the basis of likely outcome could play a significant role in the management of what might be limited healthcare resources. The development of biomarkers that correlate with disease prognosis would be a valuable clinical tool. Additionally, certain biomarkers may reveal an underlying mechanism of pathogenesis. These can feed back into basic research to identify potential therapies.

Biological phenomena, such as the progression of viral infection studied here, are temporal in nature, and the resulting observations need to be analyzed with statistical techniques that can reveal temporal differential expression. EDGE analysis of our data uncovered 35 differentially expressing proteins as compared to 27 by 2-way ANOVA. Thus our results here indicate for proteomic data, as Storey et al. (2005) first demonstrated with genomic data, that dynamic methods such as EDGE are more effective at discovering candidate biomarkers that may be hidden to static methods. It is worth noting, however, that EDGE is a univariate method and it examines the data protein by protein (peak by peak, in the case of the present method) to evaluate differential expression. These studies demonstrate that SELDI mass spectrometry can be used to determine differences between attenuated and virulent virus infection, and that the development of characteristic diagnostic and prognostic “biomarker signatures” may be of use in managing these infections. In addition, we have shown that this approach can be used as a discovery tool to determine potential mechanisms of pathogenesis and identify targets for therapeutic intervention at the level of modulation of the host immune response.

In this model system, we have used two variants of a PICV which result in a mild infection, from which the animals recover, or a severe hemorrhagic disease. Further studies will use only the virulent virus variant and investigate the differential responses between animals that succumb to infection and those which recover. This study represents an important proof of principle that mass spectrometry methods can be used to differentiate between infections that result in severe disease and those that are self-limiting, suggesting that diagnostic and prognostic assays can be developed which function by detecting host responses to infection, rather than virus material directly. This may allow the development of assays which can identify the agent of infection a priori and which may be less prone to reduction of assay sensitivity by viral mutation.

We also identified differential secretion of a host immunomodulatory peptide between infection with attenuated and virulent virus variants although we could not discriminate between prothymosin- and thymosin- 1 in this study. As infection with the attenuated virus leads to a mild-infection and viral clearance, it is likely that the correct immune signaling events are activated, leading to the development of protective immunity, events which may be suppressed in virulent infection. The use of immunomodulation and targeting of the host-responses to, and critical cellular-factors for, infection may be a viable therapeutic option in the treatment of infectious disease and may allow existing drugs to be repositioned [2832]. Decreased thymosin- 1 secretion in virulent infection may contribute to the lack of development of the protective immune response, potentially through failure of proper induction of NF-κB, MAP kinase, and MyD88 pathways [33, 34]. Thymosin- 1 is already being tested as a therapeutic against HIV/AIDS, sepsis, hepatitis B and C viruses, and several cancers [24, 3540]. Its characterized pharmacodynamics and kinetics may facilitate its development as a therapeutic for emerging infectious diseases.

Prothymosin- can also function as a transcriptional regulator. Interestingly, prothymosin- has been shown to interact with cAMP response element-binding (CREB) binding protein and can stimulate NF-κB and AP-1 driven transcription [41]. We have previously identified these transcription factors as being centrally involved in mediating host-responses to PICV infection [35]. NF-κB and AP-1 have also been implicated in mediating host-responses to other Lassa fever models and hemorrhagic fever infections such as Ebola and Rift Valley fever [4245]. These pathways have also suggested further targets for immunomodulatory therapies [31]. Modulation of CREB binding protein function is the mechanism used by the A238L protein of the hemorrhagic DNA virus African swine fever virus, to modulate pro-inflammatory gene transcription. The protein translocates to the nucleus and alters CREB binding protein and p300 function, preventing recruitment to enhanceosomes and inhibiting expression of genes including TNF- and COX-2 [4649].

We have used SELDI mass-spectrometry to compare protein “signatures” following infection of guinea pigs with attenuated and virulent virus variants. This study demonstrates proof of principle that proteomic methods can distinguish between attenuated and virulent forms of infection in vivo. We have shown that, consistent with previous studies, infection with the virulent virus results in responses more similar to mock-infection, suggesting an inhibition of host responses. We identified a peptide-prothymosin- -which contains the soluble immune mediator thymosin- 1 as being differentially expressed between infections and confirmed differential secretion. These studies shed further light on the mechanisms of hemorrhagic arenavirus pathogenesis, identify a potential strategy for therapy, and illustrate how proteomic approaches for biomarker discovery can generate additional hypotheses for further, targeted investigation.

Disclosure Statement

H. Wang and L. O. Lomas disclose employment by Vermillion Inc. during the period of performance of this study.

Acknowledgments

This research was supported by the Defense Advanced Research Projects Agency (P42296LS0000) and the National Institutes of Health (UO1 A1054827). The authors are grateful to Barry Elsom for technical assistance.