Abstract
Flaviviruses are positive-stranded RNA viruses that are a public health problem because of their widespread distribution and their ability to cause a variety of diseases in humans. West Nile virus is a mosquito-borne member of this genus and is the etiologic agent of West Nile encephalitis. Clinical manifestations of West Nile virus infection are diverse, and their pathogenic mechanisms depend on complex virus-cell interactions. In the present work, we used proteomics technology to analyze early Vero cell response to West Nile infection. The differential proteomes were resolved 24 h postinfection using two-dimensional DIGE followed by mass spectrometry identification. Quantitative analysis (at least 2-fold quantitative alteration, p < 0.05) revealed 127 differentially expressed proteins with 68 up-regulated proteins and 59 down-regulated proteins of which 93 were successfully identified. The implication for mammalian cellular responses to this neurotropic flavivirus infection was analyzed and made possible more comprehensive characterization of the virus-host interactions involved in pathogenesis. The present study thus provides large scale protein-related information that should be useful for understanding how the host metabolism is modified by West Nile infection and for identifying new potential targets for antiviral therapy.
West Nile virus (WNV)1 is a mosquito-borne flavivirus belonging to the Japanese encephalitis virus (JEV) serocomplex. The virus is maintained in nature in enzootic cycles in which it is transmitted between ornithophilic mosquitoes and avian hosts. In mammals, including humans, WNV is an encephalitic flavivirus and can cause natural infections of the central nervous system (CNS) with a neuropathogenesis involving neuroinvasiveness (ability to enter the CNS) and neurovirulence (replication within the CNS) (1). To date, no pharmacological treatment exists for WNV, and a vaccine is only available for horses.
First isolated in 1937, WNV has become endemic in Africa, the Middle East, and parts of Asia and Europe (2, 3). Phylogenetics analysis groups WNV strains into two distinct lineages. Viruses in lineage 2 are found only in Africa, whereas viruses in lineage 1 are present both in Africa and in other areas, particularly Asia and Europe. Since 1999, WNV from lineage 1 (NY99) has reached North America where, in 2002, it caused the largest arboviral meningoencephalitis outbreak ever recorded in this area (4).
It is known that flavivirus replication can cause extensive rearrangement of host cell cytoskeletal and membrane compartments leading to a “cytopathic effect” in various cell cultures of human, primate, rodent, and insect origin (5). Recent studies have revealed specific effects of viruses on cellular processes. It has been demonstrated that flaviviruses can induce cell death directly through viral replication and the production of proapoptotic proteins (6–11), but the mechanism of pathogenesis has not been elucidated.
Although neurons are regarded as the major target of WNV in vivo (2), WNV infection has been shown to induce apoptosis in different cell lines in a similar manner in vitro (12, 13). This includes a wide range of different cell types with, in particular, the African green monkey kidney continuous cell line (Vero) recommended by the World Health Organization Collaborating Center for systematic research and isolation of arboviruses as well as a substrate to develop live attenuated and inactivated vaccines. Acute infection of Vero cells by WNV produces a lytic infection with a characteristic rounding cytopathic effect and the production of a large number of infectious particles in the culture fluid within 3 days postinfection (14). Although this permissive mammalian cell system is widely used for flavivirus isolation, propagation, and titration, to date no studies have focused on identifying Vero cellular proteins whose expression has been altered by WNV infection. We considered that Vero cells could be a good model for in vitro identification of cell protein alterations with possible implication in certain pathogenic mechanisms.
In the present work, fluorescent 2D DIGE technology combined with MS analysis was used to examine the consequences of Vero cell infection by WNV. To evaluate early mammalian cell response after infection and to avoid the effect of cell death and protein degradation, the culture conditions (e.g. infectious dose and incubation time) were optimized. A total of 93 differentially expressed protein spots were identified (over ±2-fold, p < 0.05) and confirmed by fluorescent Western blot analysis. The implication for cellular responses to this flavivirus infection as well as the potential roles of certain altered identified proteins are discussed to characterize the pathophysiologic processes. This study can also provide useful clues for antiviral research.
EXPERIMENTAL PROCEDURES
Reagents
N-Hydroxysuccinimide esters Cy2, Cy3, and Cy5; urea; glycerol; mineral oil; Immobiline DryStrip gels (18 cm, pH 4–7 and 6–11); and IPG buffer solutions (pH 4–7 and 6–11) were purchased from GE Healthcare. Acrylamide, DTT, Tris, glycine, SDS, and SYPRO Ruby gel stain solution were purchased from Bio-Rad. Dimethylformamide, CHAPS, l-lysine, ammonium persulfate, iodoacetamide, agarose, bromphenol blue, and TFA were purchased from Aldrich. Thiourea, TEMED, acetone, ACN, and ethanol were purchased from Fluka (Buchs, Switzerland). Trypsin (sequencing grade) was purchased from Promega (Madison, WI). All buffers were prepared with Milli-Q water (Millipore, Bedford, MA).
Virus Production and Titration
All work with infectious virus was carried out in a biosafety level 3 laboratory. WNV strain New York 99 was provided by X. de Lamballerie (Unité des virus émergents, Faculté de médecine de la Timone, Marseille, France). The batches of viral inoculum used in this study were prepared by two passages in C6/36 cells (from Aedes albopictus; ATCC clone CRL 1660) as described previously (15). Finally cell culture supernatants were collected and frozen at −80 °C in 50 mm Hepes. Viral titers were estimated by plaque assay on Vero cells (ATTC clone CCL 81) according to the dilution method reported previously (16) and a maximum likelihood-based estimating method (17).
Infection of Vero Cells
Briefly Vero cells (ATTC clone CCL 81) were cultured in commercial Medium 199 (BioWhittaker) supplemented with 5% fetal calf serum. One day after seeding, when cultures had reached confluency, cells were infected at a multiplicity of infection (m.o.i.) of 1 for 1.5 h at 37 °C in a 5% CO2 atmosphere. Finally fresh Medium 199 was added to all wells, and the cultures were incubated at 37 °C. Twenty-four hours after virus infection, the supernatants and cells were harvested and analyzed. Mock Vero cells served as the controls and were cultured in parallel to WNV-infected cells.
Viral RNA Extraction and Real Time RT-PCR Quantification
Cell culture supernatants were collected, and viral RNA was extracted using the High Pure Viral RNA kit (Roche Diagnostics) following the manufacturer's protocol. RT-PCRs were performed on a Lightcycler thermocycler (Roche Diagnostics) using the Superscript III Platinum One-Step qRT-PCR kit (Invitrogen) as described previously (18). Specific WNV TaqMan® probes and primers used were synthesized (Eurogentec S.A., Herstal, France) in accordance with published detection assays (19).
Indirect Immunofluorescence Assay
Vero cells were grown in an eight-microchamber Lab-Tek slide (Nunc, Roskilde, Denmark) for 24 h. Cells were then either infected at an m.o.i. of 1 with WNV or left untreated. Twenty-four hours after virus infection, supernatants were collected for virus titration, and cells were fixed with 3% paraformaldehyde in PBS. Viral antigens in infected cells were detected by indirect immunofluorescence assay using WNV-specific hyperimmune mouse ascitic fluids (1:100) and Alexa Fluor 488-conjugated goat anti-mouse IgG (1:100) (Molecular Probes) (20). A DAPI counterstaining was performed to visualize cell nuclei, and the percentage of positive cells was determined using a fluorescence microscope (Zeiss Axioskop 40).
Detection of Apoptosis
Vero cells were grown in an eight-microchamber Lab-Tek slide (Nunc) and infected as described above. Twenty-four and 36 h after infection, WNV-infected cells and mock controls were fixed with 3% paraformaldehyde in PBS. After a DAPI counterstaining, apoptosis was detected by means of TUNEL analysis with fluorescein-dUTP (In Situ Cell Death Detection kit, Roche Diagnostics). Infected and mock cells were then counted under a fluorescence microscope (Zeiss Axioskop 40).
Protein Sample Preparation
Twenty-four hours after WNV infection, infected and mock cells were lysed using standard cell lysis buffer containing 8 m urea, 2 m thiourea, 4% (w/v) CHAPS, and 30 mm Tris, adjusted to pH 8.5. The cells were disrupted by ultrasonication (Vibracell 72412, Bioblock Scientific, Illkirch, France) five times for 60 s on ice at maximum amplitude. The resulting homogenate was centrifuged for 15 min at 16,000 × g at 4 °C, and the supernatant was further precipitated with 100% acetone. Protein concentrations for each sample were determined in duplicate using the Lowry method (DC Protein Assay kit, Bio-Rad) according to the manufacturer's instructions and confirmed by SDS-PAGE. Samples were resuspended in standard cell lysis buffer to obtain a protein concentration adjusted to 2.5 μg/μl.
CyDye Labeling
Proteins in each sample were minimally labeled with CyDye according to the manufacturer's recommended protocols. In short, 50 μg of each protein sample were labeled with 400 pmol of either Cy3, Cy5, or Cy2 (Cy2 was used to label the internal standard) freshly dissolved in anhydrous dimethylformamide and incubated on ice for 30 min in the dark. The reaction was quenched with 1 μl of free lysine (10 nm) by incubation for 10 min on ice. Cy3-, Cy5-, and Cy2-labeled samples were then pooled, and an equal volume of 2× sample buffer was added (8 m urea, 2 m thiourea, 4% (w/v) CHAPS, 10 mm DTT, and 1.2% (v/v) IPG buffer corresponding to the IPG strips used).
Isoelectric Focusing
DeStreak buffer containing 1.2% (v/v) IPG buffer (pH 4–7 or 6–11) was used for overnight rehydration of 18-cm IPG strips with different pH gradient ranges (4–7 or 6–11). The samples were applied at the acidic end of the IPG strips using a cup-loading technique. IEF was carried out at 20 °C for a total of 65 kV-h on an Ettan IPGphor II electrophoresis unit (GE Healthcare).
Second Dimension Electrophoresis
Prior to separation in the second dimension, the IPG strips were reduced and alkylated in equilibration buffer containing 50 mm Tris-HCl, pH 8.6 buffer, 6 m urea, 2% SDS, and 30% glycerol supplemented with 1% (w/v) DTT or 2.5% (w/v) iodoacetamide instead of DTT for 15 min. Equilibrated IPG strips were then transferred onto 10% uniform polyacrylamide gels cast in low fluorescence glass plates using an Ettan-DALT Six system (GE Healthcare). Strips were overlaid with 0.5% agarose in 1× running buffer containing bromphenol blue, and electrophoresis was run at 10 watts/gel until the dye front migrated out of the gel at 20 °C.
Image Analysis
After electrophoresis, the gels with CyDye-labeled proteins were scanned three times with a TyphoonTM Trio image scanner (GE Healthcare) each time at different excitation wavelengths (Cy3, 580 BP 30/green (532 nm); Cy5, 670 BP 30/red (633 nm); Cy2, 520 BP 40/blue (488 nm)). Prescans were performed to adjust the photomultiplier tube voltage to obtain images with a maximum intensity of 60,000–80,000 U. Images were cropped with ImageQuantTM software (GE Healthcare) and further analyzed using the software package DeCyder v6.5 (GE Healthcare). Intragel spot detection and quantification were performed using the differential in-gel analysis mode, whereas images from different gels were matched using the biological variance analysis mode. Matching between gels was performed using the in-gel standard from each image pair. The paired t test was used for statistical analysis of the data. A false discovery rate (FDR) correction was applied to eliminate false positives. Protein spots that were expressed differentially between two experimental groups (|ratio| ≥ 2, p ≤ 0.05) were marked. After 2D DIGE imaging and image analysis, the protein spot pattern was visualized by SYPRO Ruby staining (Bio-Rad) according to the manufacturer's protocol.
In-gel Digestion
Based on the DeCyder v6.5 analysis, spots of interest were excised using a manual spot picker from SYPRO Ruby stain gels. Excised plugs were digested overnight at 37 °C by sequencing grade trypsin (12.5 μg/ml; Promega) in 50 mm NH4HCO3. The resulting extracted peptides were then stored at −20 °C before their analysis by MS.
Mass Spectrometry Analysis
The samples were analyzed by nanoscale capillary LC-MS/MS. Purification and analysis were performed on a C18 capillary column using a CapLC system (Waters, Milford, MA) coupled to a hybrid quadrupole orthogonal acceleration time-of-flight tandem mass spectrometer (Q-TOF Ultima, Waters). Chromatographic separations were conducted on a reversed-phased capillary column (AtlantisTM dC18, 3 μm, 75 μm × 150 mm; Nano EaseTM, Waters) with a 180–200 nl min−1 flow. The gradient profile consisted of a linear gradient from 95% A (H2O, 0.1% HCOOH) to 60% B (80% acetonitrile, 0.1% HCOOH) in 60 min followed by a linear gradient to 95% B in 10 min. Mass data acquisitions were piloted by MassLynx 4.0 software using automatic switching between MS and MS/MS modes. The internal parameters of Q-TOF were set as follows. The electrospray capillary voltage was set to 3.2 kV, the cone voltage was set to 30V, and the source temperature was set to 80 °C. The MS survey scan was m/z 400–1300 with a scan time of 1 s and an interscan time of 0.1 s. When the intensity of a peak rose above a threshold of 15 counts, tandem mass spectra were acquired. Normalized collision energies for peptide fragmentation were set using the charge state recognition files for +2 and +3 peptide ions. The scan range for MS/MS acquisition was from m/z 50 to 1500 with a scan time of 1 s and an interscan time of 0.1 s. Fragmentation was performed using argon as the collision gas and with the collision energy profile optimized for various mass ranges and charges of precursor ions. Mass data collected during a nanoscale capillary LC-MS/MS analysis were processed using ProteinLynx Global Server 2.2 software (Waters) with the following parameters: no background substraction, smooth 3/2 Savitzky Golay, and no deisotoping to generate peak lists in the Micromass pkl format. pkl files were then fed into a local search engine, Mascot Daemon v2.2.2 (Matrix Science, London, UK). The data were searched against the primates (295,078 sequences) and viruses (435,150 sequences) National Center for Biotechnology Information non-redundant (NCBInr) protein database (December 2, 2007). Search parameters allowed for one missed tryptic cleavage site, the carbamidomethylation of cysteine, and the possible oxidation of methionine; precursor and product ion mass error tolerance was <0.2 Da. All identified proteins have a Mascot score greater than 69, corresponding to a statistically significant (p < 0.05) confident identification. Moreover among the positive matches, only protein identifications based on at least two different non-overlapping peptide sequences of more than 6 amino acids and with a mass tolerance <0.05 Da were accepted. These additional validation criteria are a good compromise to limit the number of false positive matches without missing real proteins of interest.
SDS-PAGE, Blotting, and Analysis Procedures
Immunoblotting with fluorescence-based methods including the ECL Plex-CyDye immunoblot detection system (GE Healthcare) was used to detect both the total protein expression profile and the specific immunoreactive proteins. The same protein samples used for 2D DIGE were minimally labeled with Cy3 cyanine dye as described above (see “CyDye Labeling”). Labeled samples were separated by 10% SDS-PAGE in a Mini-PROTEAN Cell (Bio-Rad) according to Laemmli (21). Gels were transferred to a nitrocellulose membrane (0.45 μm; GE Healthcare) using a semidry blotting system at 200 mA for 30 min (22). Blots were saturated with 5% nonfat dried milk in PBS containing 0.1% (v/v) Tween 20 for 1 h. Western blot analyses were carried out with anti-tubulin (1:500, mouse IgG isotype FITC conjugate) (clone DM1A, Sigma), anti-vimentin (1:500, mouse IgG isotype) (clone V9, Sigma), anti-enolase (1:1000, rabbit antibody) (Sigma), and anti-transglutaminase type 2 (1:2000, goat antibody) (Sigma) diluted in PBS containing 0.2% (v/v) Tween 20 with 5% nonfat dried milk. After 2 h of incubation at room temperature, anti-vimentin and anti-enolase were revealed with ECL Plex goat anti-mouse IgG Cy2- and goat anti-rabbit IgG Cy5-conjugated secondary antibody (GE Healthcare) (1:2000), respectively. Anti-transglutaminase was revealed with mouse anti-goat IgG FITC conjugate (1:2000) (Jackson ImmunoResearch Laboratories, West Grove, PA). All manipulations were protected from light. The gels, just after the electrophoresis and immunoblots, were scanned using a Typhoon Trio image scanner as described above (see “Image Analysis”). Immunoreactive bands were analyzed using Quantity One software (Bio-Rad). To evaluate the expression level of the different proteins, immunoreactive band intensities were normalized to the intensities of a global protein pattern labeled with Cy3. Band intensities were also corrected for adjacent background. Standard molecular weight markers were loaded in each gel (Bio-Rad).
RESULTS
Virus Infection Conditions
Optimal culture conditions (e.g. infectious dose and incubation time) were tested to determine the early host responses of Vero cells following WNV infection at the proteome level. Chu and Ng (12) have previously studied the effect of WNV infectious dose on the process of cell death. The authors demonstrated that a low infectious dose of virus (m.o.i. ≤ 1) induced a delay in the kinetics of virus replication and in apoptotic cell death in Vero cells. To obtain a high percentage of infected cells with minimal cell death in our assays, a m.o.i. of 1 and 24 h postinfection were the conditions chosen for subsequent proteomics analysis. Under these conditions, more than 90% of cells were infected by the virus, and cell death was less than 2%, comparable to the mock control cells. DNA fragmentation by TUNEL assay was controlled with mock- and WNV-infected Vero cells cultivated for 36 h with a m.o.i. of 1. In these conditions, up to 10% of infected cells were labeled with fluorescein-dUTP. A representative result of three independent experiments is shown in Fig. 1. Moreover real time RT-PCR quantification of WNV RNA in the cell culture supernatants indicated a productive virus infection (data not shown).
Evaluation of virus infection and cell viability. Mock- (a, c, and e) or WNV-infected (b, d, and f) Vero cells at 24 (a, b, c, and d) or 36 h (e and f) postinfection were observed under a fluorescence microscope. a and b, viral antigens (green) were detected by immunofluorescence assay using WNV-specific hyperimmune mouse ascitic fluids (1:100) and Alexa Fluor 488-conjugated goat anti-mouse IgG (1:100). c, d, e, and f, apoptosis (green) was revealed by the incorporation of fluorescein-dUTP at the free 3′-hydroxyl ends of fragmented DNA by TUNEL staining assay. A DAPI counterstaining was performed to visualize cell nuclei (blue).
Analysis of Differentially Expressed Proteins after WNV Infections
To identify the cell proteins involved in response to WNV infection, 2D DIGE experiments were performed. Four independent cultures of mock- and WNV-infected Vero cells cultivated for 24 h with a m.o.i. of 1 were included in this analysis. After 2D DIGE, the Cy2, Cy3, and Cy5 channels of each gel were individually imaged, and the images were analyzed using DeCyder 6.5 software.
For proteins separated in the pH 4–7 range, 1950 protein spots were matched. Of them, 79 spots were significantly modified between the mock- and WNV-infected Vero cells (33 spots were up-regulated, and 46 spots were down-regulated; |ratioWNV/mock| ≥ 2, p ≤ 0.05 after FDR correction; spots shown in Fig. 2). For proteins separated in the pH 6–11 range, 1674 protein spots were matched. Of them, 48 spots were significantly modified between the mock- and WNV-infected Vero cells (35 spots were up-regulated, and 13 spots were down-regulated; |ratioWNV/mock| ≥ 2, p ≤ 0.05 after FDR correction; spots shown in Fig. 3). A total of 127 protein spots were found to have been modified following Vero cell infection by WNV under our conditions.
2D DIGE analysis of mock- or WNV-infected Vero cells at 24 h postinfection (m.o.i. = 1). Representative data from a 2D DIGE experiment using a 10% homogenous SDS-polyacrylamide gel with the pH range from 4 to 7 is shown. Proteins from mock-infected Vero cells were conjugated to Cy5, and proteins from WNV-infected Vero cells were conjugated to Cy3. As determined by DeCyder software, protein spots that were differentially expressed between two experimental conditions (|ratio| ≥ 2, p ≤ 0.05 after an FDR correction) were marked with master numbers. Normal and italic numbers correspond to protein spots significantly up- and down-regulated following Vero cell infection by WNV, respectively. Proteins identified by mass spectrometer are listed in Table I.
2D-DIGE analysis of mock- or WNV-infected Vero cells at 24 h postinfection (m.o.i. = 1). A, representative data from a 2D DIGE experiment using a 10% homogenous SDS-polyacrylamide gel with the pH range from 6 to 11 is shown. Proteins from mock-infected Vero cells were conjugated to Cy3, and proteins from WNV-infected Vero cells were conjugated to Cy5. As determined by DeCyder software, protein spots that were differentially expressed between two experimental conditions (|ratio| ≥ 2, p ≤ 0.05 after an FDR correction) were marked with master numbers. Normal and italic numbers correspond, respectively, to protein spots up- and down-expressed following Vero cell infection by WNV. Numbers correspond to proteins identified by mass spectrometer in Table II. B, enlarged area of the gel; spots differentially expressed between mock- and WNV-infected cells are noted.
Identification of Differentially Expressed Proteins after WNV Infection
One hundred and twenty-seven differentially detected protein spots were excised manually from gels, subjected to in-gel digestion, and analyzed by LC-MS/MS. The resulting fragment ion spectra were searched against primate and viral protein databases (NCBInr). Tables I and II list proteins that were identified by tandem mass spectrometry from gels with pH ranges of pH 4–7 and pH 6–11, respectively. Ninety-three protein spots (73%) were identified with a high degree of confidence (Tables I and II and supplemental Tables 1 and 2). As expected, 15 protein spots with a significant average ratio -fold change were identified as WNV proteins such as polyprotein precursor or nonstructural protein 1 (NS1), NS3, and NS5 (Tables I and II).
The proteins were identified by mass spectrometry following in-gel trypsin digestion. The spot numbers correspond to the same numbers as indicated on Fig. 2. The identities of the spots, their NCBI accession numbers, and the theoretical molecular masses and pI values as well as the number of peptide sequences, the corresponding percent sequence coverage, and the Mascot score are listed for MS/MS analysis. (Protein scores greater than 69 are considered as significant (p < 0.05).) Paired average volume ratio and p values (ttest) between Vero cells infected by WNV versus uninfected cells (WNV/mock) were quantified using DeCyder software. n.i., no identification; TNF, tumor necrosis factor; CTBP, C-terminal-binding protein; CBP, cAMP-response element-binding protein.
The proteins were identified by mass spectrometry following in-gel trypsin digestion. The spot numbers correspond to the same numbers as indicated in Fig. 3. The identities of the spots, their NCBI accession numbers, and the theoretical molecular masses and pI values as well as the number of peptide sequences, the corresponding percent sequence coverage, and the Mascot score are listed for MS/MS analysis. (Protein scores greater than 69 are considered as significant (p < 0.05).) Paired average volume ratio and p values (t test) between Vero cells infected by WNV versus uninfected cells (WNV/mock) were quantified using DeCyder software n.i., no identification.
In addition to these viral proteins, 54 distinct cellular proteins were identified according to their NCBI number (Tables I and II). The main proteins with profile modification were localized in the cytoplasm (60.7%) and nucleus (19.6%) (Fig. 4 A). These proteins were grouped into functional categories, such as transcription and translation, cytoskeletal and cell mobility, chaperones, energy pathways, apoptosis, lipid metabolism, ubiquitin-proteasome pathway, regulation and metabolism of nucleotide, binding, and transport (Tables I and II and Fig. 4B). Several proteins (e.g. vimentin and transglutaminase 2) were detected in more than one spot, suggesting different isoforms and/or post-translational modifications. Tumor rejection antigen (gp96) protein (Table I) was detected in seven spots with an up-representation for three of them (that is 436, 437, and 442) and a down-representation for four of them (that is 1187, 1189, 1192, and 1195). These results indicated that following WNV infection of Vero cells a group of host protein isoforms (e.g. 436, 437, and 442) was privileged to the detriment of the second group of other isoforms from the same host protein. Thirty-four protein spots were not identified probably because of insufficient amounts of protein or of low MS spectra qualities.
Pie chart representing up- and down-regulated proteins identified by mass spectrometry following Vero cell infection by WNV. Proteins have been classified according to their subcellular location (A) and biological function (B). The percentages of proteins associated with each category are indicated in parentheses.
Protein Validation by Western Blotting Analysis
To verify the 2D DIGE results, Western blot analyses were performed. Several proteins identified by MS were selected on the basis of their representativity and the availability of the corresponding antibodies. For all Western blots, each protein sample was labeled with Cy3 cyanine dye revealing the variations in sample loading that were taken into account for the calculation of the average band volume ratio. The same samples used for 2D DIGE analysis were subjected to one-dimensional SDS-PAGE and blotted with the antibodies against tubulin, vimentin, transglutaminase, and enolase.
Anti-tubulin FITC-conjugated antibody confirmed that this protein was up-regulated following WNV infection of Vero cells (Table I and Fig. 5 ). For vimentin, the antibody detected two molecular weight isoforms (Fig. 5A). This is consistent with the results of 2D DIGE analysis that showed two groups of protein spots (1482 and 1483 and 1574, 1575, and 1577) corresponding to vimentin with different molecular weights. Both analyses indicated, for the two groups of protein spots, a down-representation of vimentin following WNV infection (Fig. 5B). Anti-transglutaminase antibodies detected two bands by Western blot that were up-regulated following WNV infection (Fig. 5, A and B). Although proteomics analysis identified two protein spots (741 and 742) corresponding to transglutaminase that were also up-expressed following WNV infection (Fig. 2 and Table I), these two protein spots were transglutaminase isoforms with the same molecular weight. These two protein spots could thus correspond to only one of the two bands detected on the Western blot. It could effectively be hypothesized that the expression level of the highest band and its -fold change detected on the Western blot were not sufficient for the inclusion of this isoform in the list of statistically significant protein spots differently regulated following WNV infection using proteomics analysis. Using anti-enolase antibodies, we confirmed the down-expression of this protein (Fig. 5). Thus, Western blot experiments confirmed the results obtained by 2D DIGE analysis and protein identification by mass spectrometry.
Western blot validations of differentially expressed proteins. A, the same protein samples used for DIGE analysis were minimally labeled with Cy3 cyanine dye. At the top, a representative protein profile of the four biological replicates from lysates of mock- and WNV-infected Vero cells separated by 10% SDS-PAGE is shown. Immunoblotting with fluorescence-based methods was used to detect an overlaid fluorescent scan of the general protein patterns (Cy3 dye (green)) and the specific immunoreactive proteins (FITC/Cy2 dye (yellow) or Cy5 dye (blue)). The arrows indicate the specific protein bands revealed. Standard molecular masses are indicated on the left. MW, molecular weight; kDa, kilodalton. B, the graphs correspond to the mean ± S.D. of protein quantity measured by densitometry of antigenic bands. Densitometry analyses were done using Quantity One software (Bio-Rad), and data were normalized to levels of global protein pattern intensity. Results are mean ± S.D. of the four biological replicates, and -fold change from WNV-/mock-infected is shown at the top of each graph. A.U., arbitrary units.
DISCUSSION
A clear understanding of WNV pathogenesis is lacking, so the discovery of new virus-host interactions is necessary to identify novel strategies to treat or prevent this viral infection. Despite the importance of mosquito-borne virus diseases, few studies have examined host metabolism alterations following flavivirus infection. In this way, microarray studies have analyzed changes at the transcription level in response to dengue virus (DENV) or WNV infections (23–28). To date, only two studies have analyzed the interplay between flavivirus and host cells using proteomics analysis (31, 32), and more investigations are needed to understand the mechanism of virus replication and pathogenesis. However, studies combining transcriptomics and proteomics analysis indicated that significant divergence could be observed between protein levels and mRNA abundance (29, 30). In this way, our study provides global information on protein alterations in mammalian cells infected by WNV, an emergent pathogenic flavivirus. Of the 127 spots differentially expressed in Vero cells after WNV infection according to 2D DIGE technology, 93 protein spots were identified. The identified cellular proteins function in diverse biological processes, and some functional groups are formed through cluster and pathway analysis. The potential roles of some of these altered proteins belonging to major clusters in response to WNV infection are discussed below in relation with pathogenesis and early host antiviral response.
Alteration of Cytoskeleton Networks
The altered microtubule-associated protein identified in the present work, i.e. tubulin α, was up-regulated, whereas the intermediate filament protein vimentin was down-regulated in WNV-infected Vero cells. Like many viruses (33), flaviviruses have to manipulate and utilize the host cytoskeleton to promote viral infection, but this mechanism is still unclear.
Microtubules and microtubule-associated proteins are known to play an important role in the intracellular trafficking of viral components as wells as virions in the infected host cell (34). The actin filaments were shown to be involved in mediating the internalization of WNV particles into mammalian cells by clathrin-dependent endocytosis (35). Moreover it has been demonstrated that actin filaments also play a key role in the release of West Nile (Sarafend) virions (36). Microtubules are cytoskeleton polar structures constituting α- and β-tubulin subunits that perform general functions including organelle movement and cargo transport in all kinds of cells. A previous study has shown that JEV infection induced ultrastructural changes in Vero cells with microtubule rearrangement and redistribution (37). This phenomenon could facilitate the transport of viral proteins from the rough endoplasmic reticulum (ER) and Golgi apparatus to the convoluted membrane, which may serve as a reservoir for viral proteins during JEV multiplication. Furthermore as tubulin exhibits the functions of a chaperone (38), overexpression of its cellular component may also help viral proteins to maintain their conformation in the convoluted membrane.
Vimentin is a major component of the type III intermediate filaments found in many cell lines (39). This 57-kDa protein serves to maintain cell shape and to anchor and position the nucleus within the cell. Major changes in the distribution of vimentin are observed when cells move or divide (40, 41), but this protein is also redistributed in cells expressing misfolded proteins and during virus infection (42, 43). The vimentin rearrangements with the formation of a cage structure often observed could represent a protective response by the cell during virus infection (44). In fact, the vimentin cage may have a cytoprotective role by preventing the diffusion of viral components into the cytoplasm. Several studies have shown that human immunodeficiency virus type 1 protease cleaves the intermediate filament vimentin and induces the collapse of vimentin in infected cells (45, 46). Further elucidation is in progress to determine whether the flaviviral NS2B-NS3 protease can cleave vimentin, resulting in strongly decreased levels and the collapse of the vimentin network. On the other hand, vimentin is also implicated in the regulation of cell death via its caspase-mediated cleavage (47, 48), and the observed down-regulation of the protein could also be related to the induction of apoptosis in WNV-infected Vero cells (49).
Stress Response Proteins
The quantities of two specific cellular proteins related to oxidative stress were found to have been modified in WNV-infected Vero cells. The 70-kDa heat shock protein (Hsp70) was highly expressed, whereas Hsp90 was down-regulated. Eukaryotic Hsp70s are highly abundant cytosolic and nuclear molecular chaperones that play essential roles in various aspects of protein homeostasis (50). Moreover the antiviral activity of some drugs has been associated with the induction of the specific Hsp70 (51, 52). Additionally it has been shown that Hsp70 was able to prevent the cytotoxic effects of the WNV capsid protein, suggesting a protective cell function for this molecular chaperone against WNV infection (53). However, the specific overexpression of Hsp70 by flavivirus infection has not been described until now. Thus, the present study provides the first demonstration of the up-regulation of Hsp70 after WNV infection, highlighting the possible role of this specific chaperone as a protector of the cell from apoptosis (54).
Hsp90 is a molecular chaperone that is highly expressed in most cell types in various organisms (55). It displays ATP-dependent folding capacity (56) and appears to have a specific set of client proteins in vivo unlike its promiscuous cousin chaperone Hsp70 (57, 58). Hsp90 involvement in viral replication has been reported for many viruses, and it has been demonstrated that Hsp90 inhibition blocks viral replication (59). Several viruses require viral and host molecular chaperones for entry, replication, and assembly as well as other steps in viral production (60, 61). Recently a role for Hsp90 in the control of hepatitis C, flock house, and influenza virus polymerase function has been shown (62–65), and it has been proposed that Hsp90 is a major host factor that is of central importance for viral replication for a wide spectrum of RNA viruses (66). The chaperone Hsp90 has also been identified as an essential factor in the folding and maturation of picornavirus capsid proteins (67). The importance of Hsp90 for the replication of multiple viruses opens up an interesting possibility for developing new antiviral therapies that have not yielded drug-resistant viruses (68). Although the down-regulation of Hsp90 observed in WNV-infected Vero cells was not expected, this result may still be an innate cellular response reflecting an early cellular survival mechanism against WNV infection.
RNA Processing Machinery and Host Translation Pathway
Among the differentially expressed host proteins, many of them are known to participate in viral replication and translation (Tables I and II). Positive-strand RNA viruses must recruit normal components of host cellular RNA processing or translation machineries for viral RNA synthesis and protein synthesis (69).
Human hnRNP C1/C2 are members of the heterogeneous nuclear ribonucleoprotein family, which consists of 20 major hnRNP proteins (designated hnRNP A1 through U) with molecular sizes of ∼36–120 kDa (70). These proteins are involved in mRNA biogenesis and contain important conserved motifs essential for RNA binding, protein-protein interaction, and nuclear localization. Several ribonucleoproteins, hnRNP Q, A1, A2/B, H, and YB-1 have been found to bind specifically to DENV 3′-untranslated region, suggesting that these molecules may play a biologically significant role in the DENV life cycle (71). Moreover YB-1 is directly involved in repression of DENV translation and can participate in early innate immune responses during DENV infection. A specific association has also been demonstrated between hnRNP C1/C2 and dengue virus NS1 proteins that may be favorable for virus survival in host cells (72). In the present work, hnRNP L and hnRNP U were identified for the first time as being significantly up-regulated in Vero WNV-infected cells, suggesting that these two specific host cellular proteins play a key role in WNV replication. Interestingly we also observed a depletion of the poly(rC)-binding protein (αCP-2) after WNV infection. This protein has multiple functions in post-transcriptional control of host and viral gene expression. With regard to the hepatitis A virus, it has been demonstrated that αCP-2 is cleaved by the viral protease 3C, which could down-regulate viral protein synthesis, giving way to viral RNA synthesis (73). Because of several differences between picornavirus and flavivirus replication strategies (genome cyclization and cap-dependent translation), the role of αCP-2 in the flavivirus regulation life cycle remains to be elucidated. In this study, another important regulating protein, nucleolin, was found to be more abundant in WNV-infected cells. Nucleolin is one of the most abundant nonribosomal proteins in the nucleolus, and it shuttles between nucleoli, nucleoplasm, cytoplasm, and the cell surface (74). Previous works have reported that nucleolin plays a role in virus replication with direct interaction with the hepatitis C viral RNA-dependent RNA polymerase or overexpression in herpes simplex virus-infected cells (75, 76). Based on these data, an up-regulation of nucleolin in WNV-infected Vero cells suggests that the host protein plays an important part in flavivirus biology, making possible future drug target development.
Translation factors have been well documented as playing crucial roles in viral RNA and protein synthesis (77). Moreover it is known that flaviviruses do not shut off host protein synthesis and must compete with the cellular translation machinery for limiting factors (78). In this way, the eukaryotic translation elongation factors 1A (eEF1A) have been identified as up-regulated proteins or important host cell factors for viral RNA replication during DENV or WNV infection (32, 79). In the present study, an increase in the eEF2 level was observed in the WNV-infected Vero cells suggesting that this protein may also play an important role in the translation of WNV viral RNA. Additionally several tRNA synthetases (glutamyl-prolyl/glutaminyl/tryptophanyl-tRNA synthetases) were shown to be overexpressed during WNV infection, suggesting that their overall up-regulation is required for virus propagation in Vero cells. Because flaviviruses prevent the shutoff of host cell translation, the decrease in the eukaryotic translation initiation factor 3 (eIF3) was not expected in WNV-infected Vero cells. eIF3 is one of the 11 or more initiation factors that are involved in the first stage of protein synthesis in eukaryotes. This initiation factor plays an important role by bridging the 43 S preinitiation complex to mRNA via the cap-binding complex eIF4F (80). In coronavirus-infected cells, no obvious inhibition of host protein synthesis was observed at least at early stages of the virus infection cycle. However, eIF3f has been identified as a regulator of the translation of virus-induced genes at late stages of the virus infection cycle (81). Although proteases of many picornaviruses and caliciviruses are known to cleave several translation factors (82, 83), this phenomenon had not yet been reported during flavivirus infection. Thus, the eIF3 down-expression observed in WNV-infected Vero cells may be the reflection of translation regulation mechanisms exploited by flavivirus as a means of controlling viral pathogenesis that have not been described yet.
Apoptotic Pathways
Apoptosis has opposing effects on viral pathogenesis by either preventing viral dissemination or promoting viral spread by release of the progeny virus. However, in the case of neurotropic viruses such WNV, neuron cell death has serious consequences for the host, and elucidation of virus-triggered apoptotic pathways may be critical for the understanding of the pathogenicity of flavivirus (49).
The programmed cell death 8 protein or apoptosis-inducing factor (PDCD8/AIF) is a flavoprotein localized in the mitochondrial intermembrane space, and its translocation to the nucleus induces chromosome condensation and fragmentation (84). Up-regulation of this protein in WNV-infected Vero cells reflects the induction of virus-mediated apoptosis.
Flavivirus replication is associated with extensive proliferation of intracellular membranes, which are, in part, derived from the ER. Translation of the flavivirus polyproteins is associated with the ER membranes, and the ER is also the site of viral encapsidation and envelopment. Therefore, it is likely that flavivirus infection triggers not only perturbation of ER homeostasis but also the unfolded protein response as evidenced by substantial induction of the expression of several chaperones (11). During ER stress, tumor rejection antigen, also known as gp96, is one of the important molecular chaperones that has been involved in unfolded protein response and apoptosis (85, 86). In our study, we observed a decrease in protein levels of gp96 in WNV-infected Vero cells. As already described for JEV infection (87), this phenomenon accelerates ER stress-induced apoptosis and may represent a host defense mechanism to limit viral replication.
Conclusion
Several Vero cell proteins that were modulated early in WNV infection were identified in the present study. To our knowledge, this work made it possible to identify the largest number of cellular proteins regulated differently in response to a flavivirus infection. Most of these proteins were involved in transcription/translation processes, alteration of the cytoskeleton networks, stress cellular response, and apoptotic pathways. Although the roles of some of these alterations could be related to cellular antiviral response or pathogenic mechanisms, functional significance of other altered proteins remains unclear and needs further investigation. However, the present study provides large scale protein-related information that should be useful for understanding the modulation of host metabolism in West Nile infection and for identifying new potential targets for antiviral therapy.
Acknowledgments
We thank Matthieu Pophillat and Samuel Granjeaud for technical assistance on database searching. This manuscript was reviewed and corrected by K. Snaith (Médicis Traduction), member of the French Society of Translators.
Footnotes
-
↵* This work was supported by the French Armed Forces Medical Service and the Délégation Générale pour l'Armement (Vi-Ve-Pro-Guyane Grants 06CO006 and 05CO008-05).
-
↵
The on-line version of this article (available at http://www.mcponline.org) contains supplemental
-
↵The abbrevations used are:
- WNV
- West Nile Virus
- 2D
- two-dimensional
- CNS
- central nervous system
- DENV
- dengue virus
- ER
- endoplasmic reticulum
- FDR
- false discovery rate
- JEV
- Japanese encephalitis virus
- m.o.i.
- multiplicity of infection
- TEMED
- N,N,N′,N′-tetramethylethylenediamine
- DAPI
- 4′ 6-diamidino-2-phenylindole
- TUNEL
- terminal deoxynucleotidyltransferase dUTP nick-end labeling
- BP
- band pass
- NCBInr
- National Center for Biotechnology Information non-redundant
- NS
- nonstructural protein
- Hsp
- heat shock protein
- hnRNP
- heterogeneous nuclear ribonucleoprotein
- eEF
- eukaryotic translation elongation factor
- eIF
- eukaryotic translation initiation factor
- Received December 10, 2008.
- Accepted April 21, 2009.
- © 2009 The American Society for Biochemistry and Molecular Biology
REFERENCES
- 1.↵
- 2.↵
- 3.↵
- 4.↵
- 5.↵
- 6.↵
- 7.↵
- 8.↵
- 9.↵
- 10.↵
- 11.↵
- 12.↵
- 13.↵
- 14.↵
- 15.↵
- 16.↵
- 17.↵
- 18.↵
- 19.↵
- 20.↵
- 21.↵
- 22.↵
- 23.↵
- 24.↵
- 25.↵
- 26.↵
- 27.↵
- 28.↵
- 29.↵
- 30.↵
- 31.↵
- 32.↵
- 33.↵
- 34.↵
- 35.↵
- 36.↵
- 37.↵
- 38.↵
- 39.↵
- 40.↵
- 41.↵
- 42.↵
- 43.↵
- 44.↵
- 45.↵
- 46.↵
- 47.↵
- 48.↵
- 49.↵
- 50.↵
- 51.↵
- 52.↵
- 53.↵
- 54.↵
- 55.↵
- 56.↵
- 57.↵
- 58.↵
- 59.↵
- 60.↵
- 61.↵
- 62.↵
- 63.↵
- 64.↵
- 65.↵
- 66.↵
- 67.↵
- 68.↵
- 69.↵
- 70.↵
- 71.↵
- 72.↵
- 73.↵
- 74.↵
- 75.↵
- 76.↵
- 77.↵
- 78.↵
- 79.↵
- 80.↵
- 81.↵
- 82.↵
- 83.↵
- 84.↵
- 85.↵
- 86.↵
- 87.↵