Figures
Abstract
Background
Chikungunya virus (CHIKV) is an emerging mosquito-borne alphavirus that has caused multiple unprecedented and re-emerging outbreaks in both tropical and temperate countries. Despite ongoing research efforts, the underlying factors involved in facilitating CHIKV replication during early infection remains ill-characterized. The present study serves to identify host proteins modulated in response to early CHIKV infection using a proteomics approach.
Methodology and Principal Findings
The whole cell proteome profiles of CHIKV-infected and mock control WRL-68 cells were compared and analyzed using two-dimensional gel electrophoresis (2-DGE). Fifty-three spots were found to be differentially modulated and 50 were successfully identified by MALDI-TOF/TOF. Eight were significantly up-regulated and 42 were down-regulated. The mRNA expressions of 15 genes were also found to correlate with the corresponding protein expression. STRING network analysis identified several biological processes to be affected, including mRNA processing, translation, energy production and cellular metabolism, ubiquitin-proteasome pathway (UPP) and cell cycle regulation.
Conclusion/Significance
This study constitutes a first attempt to investigate alteration of the host cellular proteome during early CHIKV infection. Our proteomics data showed that during early infection, CHIKV affected the expression of proteins that are involved in mRNA processing, host metabolic machinery, UPP, and cyclin-dependent kinase 1 (CDK1) regulation (in favour of virus survival, replication and transmission). While results from this study complement the proteomics results obtained from previous late host response studies, functional characterization of these proteins is warranted to reinforce our understanding of their roles during early CHIKV infection in humans.
Citation: Thio CL-P, Yusof R, Abdul-Rahman PSA, Karsani SA (2013) Differential Proteome Analysis of Chikungunya Virus Infection on Host Cells. PLoS ONE 8(4): e61444. https://doi.org/10.1371/journal.pone.0061444
Editor: Paul Digard, University of Edinburgh, United Kingdom
Received: January 5, 2013; Accepted: March 10, 2013; Published: April 10, 2013
Copyright: © 2013 Thio et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: This work was supported by the Postgraduate Researh Fund (Grant no.: PV088/2012A) and the University Malaya Research Grant (Grant no.: RG013-09AFR). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Chikungunya (CHIK) is a long-neglected disease that only recently began to garner attention from the scientific community following devastating outbreaks that struck India and the Indian Ocean Islands from 2004 to 2007. This disease causes substantial morbidity and an estimated death rate of 1∶1,000 [1]. Despite being perceived as a tropical disease, recent CHIK cases and sporadic outbreaks were documented in temperate regions, suggesting that this infectious disease is no longer geographically restricted to tropical countries [2]. In Malaysia, three separate outbreaks have been reported over the past 15 years [3], [4], [5].
The causative agent for CHIK infection is the chikungunya virus (CHIKV), an alphavirus belonging to the family Togaviridae [6]. CHIKV is transmitted by the mosquito Aedes aegypti and Aedes albopictus. CHIKV can be genotypically classified into the East Central South African, West African and Asian genotypes [7]. Upon infection, CHIKV causes an acute illness characterized by the classical triad of symptoms of fever, rash and debilitating arthralgia which can persist for years. However, cases from recent outbreaks saw an increasing occurrence of atypical clinical manifestations such as neurological and cardiovascular complications [8]. As there is currently no effective vaccine or antiviral regimen to combat this disease, treatment is solely palliative. All things considered, it is not surprising that CHIK is now regarded as a potential health problem in need of a solution.
Recent research efforts have focused on understanding the viral tropism and mechanisms associated with the pathogenesis of CHIK infection. In vitro studies using a panel of mammalian cell lines showed rapid induction of cytopathic effects and cell death via apoptosis in most adherent cell lines with the exception of blood-derived cell lines [9]. Autophagic process and apoptosis were also recently shown to facilitate CHIKV dissemination [10], [11]. At the molecular level, proteomics studies on CHIKV interaction with vector and mammalian host proteins have unravelled new clues in elucidating the mechanisms involved in viral replication and transmission from vector to host as well as disease progression in host cells [12], [13], [14]. Despite the extensive research, much remains to be discovered to fully comprehend the pathogenesis of CHIKV.
Contrary to the aforementioned proteomics research which investigated the late host response to CHIKV infection [13], our present study aims to identify proteins altered during early infection in the host cells by means of 2-dimensional gel electrophoresis (2-DGE). The global proteome profile of CHIKV-infected WRL-68 cells was compared with uninfected mock control cells to single out differentially expressed spots for mass spectrometric (MS) identification with subsequent Western blot validation, as well as transcript expression analysis. Results showed widespread alteration of proteins involved in several biological processes known to play essential roles in virus replication. While this study provides new insights into CHIKV pathogenesis, functional characterization of these proteins will be required to better understand their roles during early infection.
Results
Cytopathogenicity of CHIKV
The cytopathic nature of CHIKV infection in mammalian cell lines, which was reported in several studies [9], [15], [16], was observed in WRL-68 cells infected with the virus at varying MOI (MOI of 0.5, 1.0, 5.0 and 10.0) and time-points (24 and 48 h). This isolate was found to induce cytopathic effects (CPE), characterized by cell shrinkage and detachment, within 48 h of infection, as depicted in Figure 1A. CPE induction was also determined to be MOI-dependent, as cells infected at higher MOI (MOI of 5.0 and 10.0) showed more profound CPE than that of cells infected at low MOI (MOI of 0.5 and 1.0), at 48 h post-infection (p.i.). On the contrary, no significant changes in morphology were observed at 24 h p.i. at the MOI of 0.5, 1.0 and 5.0, while mild CPE was observed at the MOI of 10.0. Mock control cells were cultured in parallel and served as negative control.
(A) Morphological examination of WRL-68 cells infected at the MOI of 0.5, 1.0, 5.0 and 10.0 at 24 and 48 h incubation revealed a MOI and time-dependent induction of CPE by CHIKV. All images were captured at 100X magnification. (B) Flow cytometric quantification of percentage of cell death by AV/PI double staining of cells. Error bars indicate standard deviation of three biological replicates. (C) Flow cytometric quantification of percentage of infection by immunostaining of cells with anti-CHIK E2 mAB 3E4 (1∶100 dilution). Error bars indicate standard deviation of three biological replicates. (D) Confirmation of infection via indirect immunofluorescence assay at the optimized MOI of 5.0 at 24 h p.i. Mock cells served as negative control. All images were captured at100X magnification.
Optimization of the infection conditions for early infection study
As the aim of this study was to investigate alterations in the host cellular proteome during early CHIKV infection (i.e., the stages preceding cell death), the infection conditions (MOI and incubation time-point) were meticulously optimized to maximize infection while maintaining cell death at a minimum level. Relative quantification of percentages of infection and cell death of WRL-68 cells infected at various MOI (MOI of 0.5, 1.0, 5.0 and 10.0) for 24 and 48 h was determined by flow cytometric analysis.
The results showed that WRL-68 cells infected at the MOI of 5.0 for 24 h recorded significantly high percentage of infection at 74.77% (Figure 1B). Percentage of cell death (25.90%), albeit higher than mock control cells (14.33%), showed no significant differences when compared with cells infected at lower MOI (MOI of 0.1 and 0.5) at 24 h p.i. (Figure 1C). Furthermore, prolonging the incubation period significantly increased the percentage of cell death to more than 50%, irrespective of the MOI used. Immunostaining with anti-CHIK E2 mAb 3E4 revealed intense cytoplasmic staining in infected cells at the selected conditions, confirming infection, whereas no staining was apparent with the mock control cells (Figure 1D). Taken together, the MOI of 5.0 and 24 h incubation time-point were determined to be the optimal conditions for early CHIKV infection study.
2-DGE profiles of CHIKV infected WRL-68 cells
Comparative proteomics analysis between mock control and CHIKV-infected WRL-68 whole cell proteome was carried out using 2-DGE. Five biological replicates (n = 5) were analysed for each group. A typical gel profile for WRL-68 whole cell proteome is shown in Figure 2 (The representative proteome maps for mock control and CHIKV-infected WRL-68 cells are shown in Supplementary Figure S1). Image analysis using the ImageMaster™ 2D Platinum v7.0 software detected more than 1300 spots in each gel. Comparison of the normalized percentage spot volume between both groups revealed 53 differentially expressed spots (Fold-change>1.3, p<0.05). Of these, 44 demonstrated reduced spot intensity whereas nine exhibited increased spot intensity. All 53 protein spots were manually excised for subsequent tryptic digestion and tandem MS identification.
Forty µg of protein sample were focused on 13 cm, pH 3–10 linear IPG drystrips, followed by second dimension SDS-PAGE separation on 12.5% polyacrylamide gel which was silver stained. Five biological replicates (n = 5) for each group (Mock control and CHIKV-infected) were analyzed using ImageMaster™ 2D Platinum v7.0 software. Fifty-three spots were determined to be differentially expressed (Fold-change >1.3, p<0.05). The position of each spot is indicated by circles on the proteome map. The uppercase ‘U’ and ‘D’ denote up-regulated and down-regulated spots, respectively.
Mass spectrometric identification of differentially expressed proteins
Of the 53 protein spots subjected to MALDI-TOF/TOF identification, 50 were successfully identified, corresponding to 45 proteins (Table 1). Unique peptides identified for each protein are listed in Supplementary Table S1. Three protein spots were not identified most likely due to low abundance, resulting in low confidence score. More than one spot was identified for four proteins; guanine nucleotide-binding protein subunit beta-2-like 1 (GNB2L1), Rab GDP dissociation inhibitor beta (GDI2), eukaryotic elongation factor-2 (EEF2) and triosephosphate isomerase (TPI1)). These spots are most likely different isoforms of the protein. Functional classification based on existing information from Swiss-Prot/TrEMBL database identified proteins involved in metabolism (42.22%) and transcription/translation (17.78%) to be mainly affected by CHIKV infection (Figure 3A), whereas classification based on sub-cellular localization showed that most altered proteins to be of cytoplasmic (56.90%) and nuclear origin (17.24%) (Figure 3B).
(A) Functional categorization and (B) Sub-cellular localization of differentially modulated proteins were determined based on Swiss-Prot/TrEMBL database search.
Protein network analysis
STRING network analysis of protein-protein interactions was performed to identify functionally linked proteins and determine the potential biological processes affected [17]. The network is presented under confidence view, whereby stronger associations are represented by thicker lines or edges and vice versa, whereas proteins are represented as nodes. Twenty additional interacting proteins were added to provide a more comprehensive view of the interactions. The protein names and gene symbols used in this network are listed in Supplementary Table S2. All gene symbols were derived from the HUGO Gene Nomenclature Committee (HGNC) (http://www.genenames.org). Figure 4 shows the interaction between 45 identified proteins and the additional interactors. Thirty seven proteins were found to be linked either directly or indirectly through one or more interacting proteins, suggesting the existence of reported functional linkages. Eight biological processes were determined to be significantly involved (p<0.05 based on false discovery rate (FDR) correction) in this network, including energy production, cell cycle regulation, gene expression, mRNA metabolism, protein metabolism and modification, DNA replication and ubiquitin-protein ligase activity (Table 2).
Interaction map was generated using default settings (Medium confidence of 0.4 and 7 criteria for linkage: neighbourhood, gene fusion, co-occurrence, co-expression, experimental evidences, existing databases and text mining). Twenty additional interplay proteins were also added to each network. The protein names and gene symbols used in this network are listed in Supplementary Table S2.
Immunoblot validation of proteomics data
Two proteins, CDK1 and PDHA1, representing the down- and up-regulated groups respectively, were randomly selected for Western blot validation. GAPDH was used as the loading control for PDHA1 as both PDHA1 and ACTB have similar molecular mass of ∼43 kDa, and thus, cannot be stained together on the same blot. Immunoblots confirmed their down- and up-regulation, as shown in Figures 5A and 5B. Densitometric analysis revealed fold differences of −1.42 and 1.72 CDK1 and PDHA1 respectively (Figures 5C and 5D), which was comparable to the observed −1.77 and 1.96 fold-changes in 2-DGE analysis.
Confirmation of the expression profiles for CDK1 (A) and PDHA1 (B) was performed via immunoblot analysis. Densitometric analysis of the mean relative intensity (n = 3) for each target protein showed down-regulation of CDK1 by 1.42 fold (C) and up-regulation of PDHA1 by 1.72 fold (D). The intensity for CDK1 and PDHA1 was normalized against ACTB and GAPDH, respectively. Error bars indicate standard deviation of three biological replicates. * indicates significant difference in expression (p<0.05).
Transcript expression analysis of selected altered proteins
The transcript expression of 36 selected proteins was evaluated using real-time qPCR (The gene names and primer sequences are listed in Supplementary Table S3). All primers had amplification efficiencies within the acceptable range of 90 to 110% (Slope values between −3.1 to −3.6). It is a known fact that mRNA expression do not always correlate with protein expression [18]. In our study however, the direction of mRNA and protein expression changes of 15 proteins including CDK1 and PDHA1 were the same (Table 3). On the other hand, the transcript expression of four other proteins; adenine phosphoribosyltransferase (APRT), electron transport flavoprotein subunit alpha (ETFA), actin-related protein 2/3 complex subunit 2 (ARPC2) and cyclophilin A (PPIA), showed the opposite direction of expression change despite being statistically significant. Meanwhile, the mRNA expression levels of 17 other proteins showed no statistically significant differences.
Discussion
It is well-established that CHIKV induces rapid and profound CPE in human host cells which culminate in cell death via apoptosis. The events preceding the inevitable cell demise, however, remain ill-characterized. A previous proteomic study on new-born mice focused on investigating the dynamic overview of altered protein expression during late stages of CHIKV infection, whereby alterations of stress, inflammation, urea cycle, energy metabolism and apoptotic-related proteins were implicated in the observed disease pathogenesis [13]. In this study, we shifted the focus to examining global changes of the host cell proteome during early CHIKV infection, with aims of identifying key proteins that are potentially involved in facilitating CHIKV replication. It has been reported that during early infection, viral replication and dissemination occurs rapidly through manipulation of the host cell machinery owing to the simplicity of the viral makeup [19]. By collating data from proteomics and bioinformatics analyses, we inferred the potential manipulation or subversion of various important cellular processes including mRNA and protein metabolism, energy production, ubiquitin-proteasome pathway (UPP) and cell cycle regulation by CHIKV.
Alteration of proteins involved in mRNA processing and translation machinery
Virus hijacking of the host mRNA processing and translational machinery is an essential process for virus replication. Viruses with positive sense RNA in particular, have been shown to recruit components of the host protein biosynthesis machineries for viral RNA and protein synthesis [20]. In the current study, we identified several deregulated proteins involved in mRNA processing and translation, including heterogeneous nuclear ribonucleoproteins C1/C2 (hnRNP C1/C2), poly(rC)-binding protein 1 (hnRNP E1), elongation factor- 2 (EEF-2), translation initiation factor EIF-2B subunit alpha (eIF2B1) and eukaryotic translation initiation factor 3 subunit H (eIF3H).
Heterogeneous ribonucleoproteins (hnRNPs) are complexes of RNA and proteins involved in an array of cellular functions such as transcription, pre-mRNA processing and cytoplasmic mRNA translation and turnover [21]. In our study, hnRNP C1/C2 was found to be up-regulated by 3.10 fold while hnRNP E1 was down-regulated by 1.42 fold. Transcript level of hnRNP E1 mRNA showed similar down-regulation while the mRNA expression of hnRNP C1/C2 was not significantly altered, suggesting that post-transcriptional and post-translational modification may play a role in modulating the expression of the latter protein. In a previous study, hnRNP C1/C2 was shown to promote dengue virus survival in host cells [22] while hnRNP E1 inhibits vesicular somatitis virus replication [23]. Ergo, the up-regulation of hnRNP C1/C2 in the present study may signify its recruitment by CHIKV whereas hnRNP E1 may possibly exert negative effects towards CHIKV propagation which is counteracted by its inhibition.
Translation factors are known to play crucial roles in viral RNA and protein synthesis and different viruses exert different mechanisms to modulate host translational proteins to their benefit, as shown in several studies [24], [25], [26]. Alphaviruses have been shown to induce global shutoff of protein synthesis by inhibiting or modifying host translational factors [27]. CHIKV-induced host translational shutoff was recently shown to occur, through an unidentified protein kinase R (PKR)-independent mechanism [28]. In this study, down-regulation of proteins involved in initiation of translation (eIF2B1 and eIF3SH) and elongation of the newly synthesised polypeptide chain (EEF-2) was observed, although at the transcript level, only EEF-2 and eIF2B1 genes were down-regulated. The exact roles of these proteins in host translational shutoff, however, cannot be ascertained at this point. Nonetheless, down-regulation of these proteins may inhibit the host translational machinery to a certain extent, possibly contributing to the observed down-regulation of most altered proteins in this study.
Differential expression of proteins involved in cellular energy production and metabolism
Of the 19 regulated proteins identified to be involved in cellular metabolism, 18 were down-regulated. Only PDHA1, a subunit of the pyruvate dehydrogenase complex involved in transforming pyruvate to acetyl-CoA in the tricarboxylic acid (TCA) cycle [29], was up-regulated by 1.96 fold. Up-regulation of this protein was further confirmed by immunoblot (Figures 5B and 5D). Transcript expression study on 14 selected genes revealed that 8 genes; PHDA1, alpha-enolase (ENO1), isocitrate dehydrogenase (IDH1), isopentyl-diphosphate Delta-isomerase 1 (IDI1), adenylosuccinate synthetase isozyme 2 (ADSS), ribose-phosphate pyrophosphokinase 1 (PRPS1), S-methyl-5-thioadenosine phosphorylase (MTAP) and phosphoserine aminotransferase (PSAT1), had expression changes of the same directionality as the protein expression (Table 1).
Based on the proteomics analysis, energy production in WRL-68 cells was expected to be significantly affected through reduced expression of glycolytic enzymes including ENO1, TPI1 and phosphoglycerate mutase 1 (PGAM1), as well as down-regulation of IDH1 which catalyzes the oxidative decarboxylation of isocitrate to alpha-ketoglutarate in the TCA cycle [30]. Four proteins associated with the adenine salvage pathway, namely PRPS1, ADSS, MTAP and adenine phosphoribosyltransferase (APRT), were also down-regulated. Similar dysregulation was observed with IDI1 and hydroxymethylglutaryl-CoA synthase (HMGCS1), two key enzymes involved in the biosynthesis of cholesterol, coenzyme Q and isoprenylated proteins through the mevalonate pathway [31].
Effects on proteins involved in the UPP
UPP is an essential intracellular system for protein degradation, with multiple cellular functions including cell cycle regulation, apoptosis, DNA repair, signal transduction and transcriptional regulation [32]. Many viruses have been reported to evolve different strategies to utilize the UPP for various purposes, including avoidance of host immune surveillance, viral maturation, viral progeny release and transcriptional regulation [33], [34], [35]. Our proteomics data showed down-regulation of two UPP associated proteins; ubiquitin-conjugating enzyme E2 N (UBE2N) and proteasome subunit alpha type-6 (PSMA6). At the transcript level however, only UBE2N showed the same direction of expression change as the protein expression. UBE2N is a ubiquitin-carrier enzyme which carries and binds ubiquitin to the ubiquitin-ligase enzyme for subsequent ubiquitination of targeted proteins. PSMA6 is the subunit of the 20S proteasome subcomplex which forms the multicatalytic 26S proteasome that degrades polyubiquitinated proteins into smaller peptides [32].
Down-regulation of proteins involved in cell cycle regulation
Cyclin-dependent kinases (CDKs) are a family of cyclin-activated serine/threonine kinases involved in various cellular processes including regulation of cell cycle (CDK1, 2, 3, 4, 6 and 7), neuronal functions (CDK5) and transcription (CDK7, 8 and 9) [36]. While CDKs are commonly associated with nuclear replication of DNA and RNA viruses, several studies have expanded the role of CDKs to cytoplasmic replication of RNA viruses as well [37], [38]. In this study, CDK1 was found to be down-regulated, both at the protein and gene expression level. CDK1 is activated by cyclin B and functions in allowing entry into mitosis from the G2 phase [39]. Inhibition of this protein would cause cell cycle arrest at G2 phase. Meanwhile, SET protein is a phosphoprotein found to regulate the cell cycle by inhibiting cyclin B-CDK1 activity [40]. In our study, SET protein was found to be up-regulated at both the protein and transcript levels, which favours the inhibition of cyclin B-CDK1 activity.
In conclusion, our proteomics data suggested that during early infection, CHIKV affects the expression of proteins involved in mRNA processing, host metabolic machinery, UPP, and cyclin-dependent kinase 1 (CDK1) regulation (in favour of virus survival, replication and transmission). While results from this study complement the proteomics results obtained from previous late host response studies, functional characterization of these proteins is warranted to reinforce our understanding of their roles during early CHIKV infection in humans.
Materials and Methods
Cell lines
WRL-68 human hepatic cells, a HeLa derivative cell line that is highly susceptible to CHIKV infection (ATCC Cat No. CL-48), Vero cells (ATCC Cat. No. CCL-81), and C6/36 Aedes albopictus cells (ATCC Cat. No. CRL-1660) were used in this study. WRL-68 and Vero cells were cultured in DMEM medium (GIBCO, Grand Island, NY) supplemented with 10% heat-inactivated fetal bovine serum (FBS) (GIBCO, Grand Island, NY) at 37 °C. C6/36 cells were grown in L-15 medium (Sigma Aldrich, St Louis, MO) supplemented with 10% tryptose phosphate broth (TPB) (Sigma Aldrich, St Louis, MO) and 10% FBS at 28 °C.
Antibodies
The antibodies used for indirect immunofluorescence assay (IIFA) and immunostaining by flow cytometry were anti-CHIK E2 monoclonal antibody (mAb) 3E4 (a kind gift from Dr. Philippe Desprès from the Pasteur Institute of France) and FITC-conjugated goat anti-mouse IgG secondary antibody (Novus Biologicals, Littleton, CO). The primary antibodies used for Western blot validation were mouse mAb to beta-actin (ACTB), cyclin-dependent kinase 1 (CDK1), glyceraldehyde 3-phosphate dehydrogenase (GAPDH) or pyruvate dehydrogenase (PDHA1). Horseradish peroxidise (HRP)-conjugated goat anti-mouse IgG was used as the secondary antibody. All antibodies used for validation were purchased from Santa Cruz Biotechnology, Santa Cruz, CA.
Virus stock propagation and titration
CHIK/06/08 clinical isolate of the ECSA genotype was propagated twice in C6/36 cell line and virus stock was harvested from the culture supernatant and stored at −80 °C. Mock control cells were cultured in parallel but without virus infection and processed in the same manner. Virus titer was determined by standard plaque assay procedure on Vero cells. Titers were expressed as plaque-forming units (PFU)/ml.
Infection of WRL-68 cells with CHIKV
WRL-68 cells were infected with CHIKV at the MOI of 0.5, 1.0, 5.0 and 10.0 for 2 h at 37 °C. Mock control cells were incubated in parallel with culture supernatant of mock control C6/36 cells. Viral inoculum was subsequently removed and the cells were further incubated in DMEM maintenance medium containing 2% FBS for 24 and 48 h. The optimal MOI and time-point for early infection study were selected based on flow cytometric quantitative analysis of percentage of cell infection and cell death [24].
IIFA
Prior to flow cytometric quantification, CHIKV infection in WRL-68 cells was confirmed by IIFA, as previously described [41] with modifications. WRL-68 cells were seeded overnight at a density of 1.5×105 cells/well in a 24-well culture dish, and subsequently infected at various MOI. Mock control cells were cultured in parallel. After 24 and 48 h incubation, the cells were fixed with 3.7% formaldehyde in phosphate buffered saline (PBS) for 20 min, washed with PBS and permeabilized with 0.15 M glycine for 10 min. Permeabilized cells were washed extensively and further incubated with anti-CHIK E2 mAb 3E4 (1∶100 dilution) for 30 min at 37 °C. Thereafter, the cells were washed with PBS and incubated in FITC-conjugated secondary antibody (1∶1000 dilution) for 30 min at 37 °C. The cells were observed under an inverted microscope (Nikon Eclipse Ti-5, Japan) and fluorescent pictures were acquired using NIS-Elements imaging software (Nikon, Japan).
Flow cytometric quantification of percentage CHIKV infection and cell death
Quantification of percentage infection was carried out as previously described [24] with modifications. Mock control and CHIKV-infected cells were harvested at appropriate time-points and fixed with 3.7% formaldehyde for 30 min. The cells were washed with staining buffer (0.1% (w/v) sodium azide in 1% FBS, pH 7.5), and incubated with anti-CHIK E2 mAb 3E4 (1∶100 dilution) for 90 min at 37 °C. Thereafter, the cells were washed and further incubated in FITC-conjugated goat anti-mouse IgG secondary antibody (1∶1000 dilution) for 60 min at 37 °C. After extensive washing, the cells were resuspended in PBS and analyzed with BD FACSCanto II flow cytometer (BD Biosciences, San Jose, CA) using FACSDiva v6.1 software.
Percentage cell death was determined using FITC Annexin V Apoptosis Detection Kit I (BD Biosciences, San Jose, CA) according to the manufacturer's protocol. Annexin V/propidium iodide stained cells were analyzed by flow cytometry.
Protein sample processing
Whole cell proteome were extracted on ice with lysis buffer (7 M Urea, 2 M Thiourea, 4% CHAPS, 2% IPG Buffer, 40 mM DTT). Cellular debris was pelleted at 17,000×g and protein supernatant was cleaned using 2-D Clean-Up Kit (GE Heathcare, Uppsala, Sweden) as described by the manufacturer. Protein estimation was performed using Bradford Protein assay (Bio-Rad Laboratories, Richmond, CA).
2-DGE
Forty µg (for analytical gel) and 160 µg (for preparative gel) of protein was mixed with rehydration buffer (7 M Urea, 2 M Thiourea, 2% CHAPS, 0.5% IPG Buffer, 1% Bromophenol blue) to a final volume of 250 µl and left overnight to rehydrate into 13 cm pH 3–10 linear immobilized pH gradient DryStrips (GE Healthcare, Uppsala, Sweden). First dimension isoelectric focusing was performed at 20 °C according to the following protocol: (i) 500 Vh, 500 V (Step-and-hold), (ii) 1,000 Vh, 1,000 V (Gradient), (iii) 16,000 Vh, 8,000 V (Gradient) and (iv) 12,000 Vh, 8,000 V (Step-and-hold). The strips were subsequently equilibrated with equilibration buffer (6 M Urea, 75 mM Tris-HCl pH 8.8, 29.3% Glycerol, 2% SDS, 0.002% Bromophenol Blue) containing 1% DTT for 15 min, followed by equilibration with equilibration buffer containing 2.5% iodoacetamide for another 15 min. Proteins were resolved on 12.5% SDS-PAGE homogenous gels at 50 V for 30 min, and 500 V for 2 h. Gels were silver stained according to a modified, MS-compatible silver staining protocol [42].
Differential gel analysis
Gels were scanned with ImageScanner™ III (GE Healthcare, Uppsala, Sweden) and analyzed using ImageMaster™ 2D Platinum v7.0 software (Amersham Biosciences, Sweden). Ten gels were used for analysis (five biological replicates per group). The volume of each spot was normalized against the total volume of all spots in the gel, and the normalized values were expressed as percentage spot volume. Spots having a fold-change of at least 1.3 and p<0.05 (as determined by one-way ANOVA and Student's t-test) were excised from multiple preparative gels for in-gel digestion.
In-gel tryptic digestion
In-gel digestion was performed using Trypsin Gold (Promega, Madison, WI) as previously described [43], [44]. Briefly, excised spots were destained with destaining solution (15 mM potassium ferricyanide/50 mM sodium thiosulphate), followed by reduction with 10 mM DTT/100 mM ammonium bicarbonate for 30 min at 60 °C and alkylation with 55 mM iodoacetamide/100 mM ammonium bicarbonate for 20 min in the dark. The gel plugs were washed trice with 50% acetonitrile (ACN)/100 mM ammonium bicarbonate, 20 min each wash, and dehydrated with 100% ACN for 20 min. The gel plugs were subsequently dried using a vacuum centrifuge (HetoVac VR-1 vacuum concentrator, Birkercd, Denmark), and digested overnight in 25 µl of 10 ng/µl trypsin at 37 °C. Tryptic peptides were then extracted twice, first with 50% ACN for 15 min, followed by 100% ACN for another 15 min. The extracted solutions were pooled together into a clean tube and dried using a vacuum centrifuge.
MALDI-TOF/TOF analysis
Dried peptides were reconstituted in 0.1% formic acid (FA) and desalted using ZipTip C18 (Millipore, Billerica, USA), according to the manufacturer's protocol. Following ZipTip cleanup, the peptides were eluted out in 2 µl elution solution (50% ACN/0.1% FA) and mixed with saturated α-cyano-4-hydroxycinnamic acid (CHCA) matrix prepared in 50% ACN/0.1% trifluoroacetic acid (TFA), at a 1∶1 ratio. Peptides were spotted on stainless-steel sample target plate in 0.7 µl aliquots in duplicates. Mass spectra for each peptide were obtained on a MALDI-TOF/TOF (ABI 4800 Plus, Applied Biosystems™, Foster City, CA) mass spectrometer using a previously established setting [43]. The spectra were analyzed with the Global Protein Server (GPS) explorer 3.6 software (Applied Biosystems™, Foster City, CA), which uses an internal MASCOT program (Matrix Science, London, UK) to match the MS and MS/MS data against existing database information. The data obtained were searched against human databases downloaded from the Swiss-Prot/TrEMBL homepage (http://www.expasy.ch/sprot).
Bioinformatics
Categorization of functional and sub-cellular distribution of proteins was performed based on Swiss-Prot/TrEMBL database search. Protein-protein interactions were predicted using Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) database v9.0 (http://www.string-db.org/). The Swiss-Prot identifier for the genes (eg. ENOA_HUMAN for alpha-enolase), in ‘Protein mode’, was used to search against the STRING database. Network analysis was set at medium stringency (STRING score = 0.4). Proteins were linked based on seven criteria; neighbourhood, gene fusion, co-occurrence, co-expression, experimental evidences, existing databases and textmining.
Western blot
Samples of CHIKV-infected and mock control cells from three independent biological replicates (not used for 2-DGE analysis) were lysed with RIPA buffer (25 mM Tris pH 7.6, 150 mM NaCl, 1.0% Triton-X, 1.0% Sodium deoxycholate, 0.1% SDS) and quantified using BCA Protein Assay Kit (Pierce, Rockford, IL). Denatured proteins (20 µg) from each sample were loaded into each lane and resolved on 12.5% polyacrylamide gels at a constant voltage of 100 V. The resolved proteins were electroblotted onto PVDF-membranes at a constant current of 80 mA for 1 h 30 min. Non-specific bindings were blocked overnight at 4 °C with 5% w/v non-fat powdered milk in Tris-buffered saline with Tween-20 (TBST) solution (50 mM Tris pH 7.4, 150 mM NaCl, 0.05% Tween-20). After extensive washing, the membranes were incubated with either mouse mAb to ACTB, CDK1, GAPDH or PDHA1 (1∶500 dilution) for 1 h 30 min at room temperature. Subsequently, the membranes were incubated with HRP-conjugated goat anti-mouse IgG (1∶2,500 dilution) for 1 h at room temperature. Target proteins were detected with TMB Stabilized Substrate for HRP (Promega, Madison, WI). The blots were scanned using ImageScanner™ III in reflective mode and densitometric quantification was performed using ImageJ v1.45 freeware (http://rsbweb.nih.gov/ij). The mean relative density for each target band was normalized against ACTB or GAPDH.
Real-time quantitative PCR (qPCR)
Total RNA of CHIKV-infected and mock control cells from three biological replicates was extracted using Qiagen RNeasy Mini Kit (Qiagen, Valencia, CA) as described by the manufacturer. Purity of extracted RNA was determined by measuring the A260/A280 and A230/A260 absorbance ratio using GeneQuant™ 1300 spectrophotometer (GE Healthcare, Uppsala, Sweden). RNA integrity was confirmed by visualization of distinct 18S and 28S ribosomal RNA bands resolved on 1% agarose gel electrophoresis. One µg of high quality RNA was converted to cDNA using High Capacity RNA-to-cDNA Kit (Applied Biosystems™, Foster City, CA), following the manufacturer's protocol. Primers specific for the gene of interest were designed with Primer3 Input v4.0 (http://frodo.wi.mit.edu/primer3/) and primer efficiency test was performed for each primer pair to confirm specificity towards the gene of interest. RNA sample (10 ng) was mixed with the respective primer pair and Fast SYBR®Green Master Mix (Applied Biosystems™, Foster City, CA). Real-time qPCR was performed using StepOnePlus™ Real-Time PCR System (Applied Biosystems™, Foster City, CA). The expression level of each target gene was normalized against ACTB. Statistical significance of altered gene expression was determined using Student's t-test, where the significance was defined at p<0.05.
Supporting Information
Figure S1.
The proteome maps of differentially expressed whole cell proteins in mock control and CHIKV-infected WRL-68 cells.
https://doi.org/10.1371/journal.pone.0061444.s001
(DOCX)
Table S1.
List of peptide sequences identified by MALDI-TOF/TOF MS.
https://doi.org/10.1371/journal.pone.0061444.s002
(DOCX)
Table S2.
Protein names and abbreviations used in STRING network analysis.
https://doi.org/10.1371/journal.pone.0061444.s003
(DOCX)
Table S3.
List of primer sequences used in real-time qPCR analysis.
https://doi.org/10.1371/journal.pone.0061444.s004
(DOCX)
Acknowledgments
We would like to thank Dr. Philippe Després from the Pasteur Institute of France for his generous contribution of the anti-CHIK E2 mAb 3E4. This research had been facilitated by access to the University of Malaya Centre for Proteomics Research (UMCPR) and Medical Biotechnology Laboratory, Faculty of Medicine, University of Malaya.
Author Contributions
Conceived and designed the experiments: SAK PSA. Performed the experiments: CLPT. Analyzed the data: CLPT SAK PSA. Contributed reagents/materials/analysis tools: SAK PSA RY. Wrote the paper: CLPT SAK PSA.
References
- 1. Kariuki Njenga M, Nderitu L, Ledermann JP, Ndirangu A, Logue CH, et al. (2008) Tracking epidemic Chikungunya virus into the Indian Ocean from East Africa. J Gen Virol 89: 2754–2760.
- 2. Nkoghe D, Kassa RF, Caron M, Grard G, Mombo I, et al. (2012) Clinical forms of chikungunya in Gabon, 2010. PLoS Negl Trop Dis 6: e1517.
- 3. Lam SK, Chua KB, Hooi PS, Rahimah MA, Kumari S, et al. (2001) Chikungunya infection-an emerging disease in Malaysia. Southeast Asian Journal of Tropical Medicine and Public Health 32: 447–451.
- 4. Chua KB (2010) Epidemiology of chikungunya in Malaysia: 2006-2009. Medical Journal of Malaysia 65: 277–282.
- 5. Noridah O, Paranthaman V, Nayar SK, Masliza M, Ranjit K, et al. (2007) Outbreak of chikungunya due to virus of Central/East African genotype in Malaysia. Medical Journal of Malaysia 62: 323–328.
- 6. Schuffenecker I, Iteman I, Michault A, Murri S, Frangeul L, et al. (2006) Genome microevolution of chikungunya viruses causing the Indian Ocean outbreak. PLoS Med 3: e263.
- 7. Volk SM, Chen R, Tsetsarkin KA, Adams AP, Garcia TI, et al. (2010) Genome-scale phylogenetic analyses of chikungunya virus reveal independent emergences of recent epidemics and various evolutionary rates. J Virol 84: 6497–6504.
- 8. Rajapakse S, Rodrigo C, Rajapakse A (2010) Atypical manifestations of chikungunya infection. Trans R Soc Trop Med Hyg 104: 89–96.
- 9. Sourisseau M, Schilte C, Casartelli N, Trouillet C, Guivel-Benhassine F, et al. (2007) Characterization of reemerging chikungunya virus. PLoS Pathog 3: e89.
- 10. Krejbich-Trotot P, Gay B, Li-Pat-Yuen G, Hoarau JJ, Jaffar-Bandjee MC, et al. (2011) Chikungunya triggers an autophagic process which promotes viral replication. Virol J 8: 432.
- 11. Krejbich-Trotot P, Denizot M, Hoarau JJ, Jaffar-Bandjee MC, Das T, et al. (2011) Chikungunya virus mobilizes the apoptotic machinery to invade host cell defenses. FASEB J 25: 314–325.
- 12.
Tchankouo-Nguetcheu S, Khun H, Pincet L, Roux P, Bahut M, et al.. (2010) Differential protein modulation in midguts of Aedes aegypti infected with chikungunya and dengue 2 viruses. PLoS One 5.
- 13. Dhanwani R, Khan M, Alam SI, Rao PV, Parida M (2011) Differential proteome analysis of Chikungunya virus-infected new-born mice tissues reveal implication of stress, inflammatory and apoptotic pathways in disease pathogenesis. Proteomics 11: 1936–1951.
- 14. Abere B, Wikan N, Ubol S, Auewarakul P, Paemanee A, et al. (2012) Proteomic analysis of chikungunya virus infected microgial cells. PLoS One 7: e34800.
- 15. Ozden S, Huerre M, Riviere JP, Coffey LL, Afonso PV, et al. (2007) Human muscle satellite cells as targets of Chikungunya virus infection. PLoS One 2: e527.
- 16. Salvador B, Zhou Y, Michault A, Muench MO, Simmons G (2009) Characterization of Chikungunya pseudotyped viruses: Identification of refractory cell lines and demonstration of cellular tropism differences mediated by mutations in E1 glycoprotein. Virology 393: 33–41.
- 17. von Mering C, Huynen M, Jaeggi D, Schmidt S, Bork P, et al. (2003) STRING: a database of predicted functional associations between proteins. Nucleic Acids Res 31: 258–261.
- 18. Guo Y, Xiao P, Lei S, Deng F, Xiao GG, et al. (2008) How is mRNA expression predictive for protein expression? A correlation study on human circulating monocytes. Acta Biochimica et Biophysica Sinica 40: 11.
- 19.
Krasner R (2009) The Microbial Challenge: Science, Disease, and Public Health, Second Edition. Massachusetts: Jones and Bartlett Publishers. 476 p.
- 20. Li HP, Huang P, Park S, Lai MM (1999) Polypyrimidine tract-binding protein binds to the leader RNA of mouse hepatitis virus and serves as a regulator of viral transcription. Journal of Virology 73: 772–777.
- 21. Krecic AM, Swanson MS (1999) hnRNP complexes: composition, structure, and function. Current Opinion in Cell Biology 11: 363–371.
- 22. Noisakran S, Sengsai S, Thongboonkerd V, Kanlaya R, Sinchaikul S, et al. (2008) Identification of human hnRNP C1/C2 as a dengue virus NS1-interacting protein. Biochemical and Biophysical Research Communications 372: 67–72.
- 23. Dinh PX, Beura LK, Panda D, Das A, Pattnaik AK (2011) Antagonistic effects of cellular poly(C) binding proteins on vesicular stomatitis virus gene expression. Journal of Virology 85: 9459–9471.
- 24. Pattanakitsakul SN, Rungrojcharoenkit K, Kanlaya R, Sinchaikul S, Noisakran S, et al. (2007) Proteomic analysis of host responses in HepG2 cells during dengue virus infection. J Proteome Res 6: 4592–4600.
- 25. Pastorino B, Boucomont-Chapeaublanc E, Peyrefitte CN, Belghazi M, Fusai T, et al. (2009) Identification of cellular proteome modifications in response to West Nile virus infection. Mol Cell Proteomics 8: 1623–1637.
- 26. Kushner DB, Lindenbach BD, Grdzelishvili VZ, Noueiry AO, Paul SM, et al. (2003) Systematic, genome-wide identification of host genes affecting replication of a positive-strand RNA virus. Proceedings of the National Academy of Sciences of the United States of America 100: 15764–15769.
- 27. McInerney GM, Kedersha NL, Kaufman RJ, Anderson P, Liljestrom P (2005) Importance of eIF2alpha phosphorylation and stress granule assembly in alphavirus translation regulation. Molecular Biology of the Cell 16: 3753–3763.
- 28. White LK, Sali T, Alvarado D, Gatti E, Pierre P, et al. (2011) Chikungunya virus induces IPS-1-dependent innate immune activation and protein kinase R-independent translational shutoff. Journal of Virology 85: 606–620.
- 29. Dahl HH, Hunt SM, Hutchison WM, Brown GK (1987) The human pyruvate dehydrogenase complex. Isolation of cDNA clones for the E1 alpha subunit, sequence analysis, and characterization of the mRNA. Journal of Biological Chemistry 262: 7398–7403.
- 30. Christensen BC, Smith AA, Zheng S, Koestler DC, Houseman EA, et al. (2011) DNA methylation, isocitrate dehydrogenase mutation, and survival in glioma. Journal of the National Cancer Institute 103: 143–153.
- 31. Bentinger M, Tekle M, Dallner G (2010) Coenzyme Q-biosynthesis and functions. Biochem Biophys Res Commun 396: 74–79.
- 32. Glickman MH, Ciechanover A (2002) The ubiquitin-proteasome proteolytic pathway: destruction for the sake of construction. Physiological Reviews 82: 373–428.
- 33. Gao G, Luo H (2006) The ubiquitin-proteasome pathway in viral infections. Canadian Journal of Physiology and Pharmacology 84: 5–14.
- 34. Leong WF, Chow VT (2006) Transcriptomic and proteomic analyses of rhabdomyosarcoma cells reveal differential cellular gene expression in response to enterovirus 71 infection. Cell Microbiol 8: 565–580.
- 35. Kanlaya R, Pattanakitsakul SN, Sinchaikul S, Chen ST, Thongboonkerd V (2010) The ubiquitin-proteasome pathway is important for dengue virus infection in primary human endothelial cells. J Proteome Res 9: 4960–4971.
- 36. Schang LM (2004) Effects of pharmacological cyclin-dependent kinase inhibitors on viral transcription and replication. Biochimica et Biophysica Acta 1697: 197–209.
- 37. Poggioli GJ, DeBiasi RL, Bickel R, Jotte R, Spalding A, et al. (2002) Reovirus-induced alterations in gene expression related to cell cycle regulation. Journal of Virology 76: 2585–2594.
- 38. Ray RB, Steele R, Meyer K, Ray R (1997) Transcriptional repression of p53 promoter by hepatitis C virus core protein. J Biol Chem 272: 10983–10986.
- 39. Lindqvist A, Rodriguez-Bravo V, Medema RH (2009) The decision to enter mitosis: feedback and redundancy in the mitotic entry network. J Cell Biol 185: 193–202.
- 40. Canela N, Rodriguez-Vilarrupla A, Estanyol JM, Diaz C, Pujol MJ, et al. (2003) The SET protein regulates G2/M transition by modulating cyclin B-cyclin-dependent kinase 1 activity. Journal of Biological Chemistry 278: 1158–1164.
- 41. Brehin AC, Rubrecht L, Navarro-Sanchez ME, Marechal V, Frenkiel MP, et al. (2008) Production and characterization of mouse monoclonal antibodies reactive to Chikungunya envelope E2 glycoprotein. Virology 371: 185–195.
- 42. Yan JX, Wait R, Berkelman T, Harry RA, Westbrook JA, et al. (2000) A modified silver staining protocol for visualization of proteins compatible with matrix-assisted laser desorption/ionization and electrospray ionization-mass spectrometry. Electrophoresis 21: 3666–3672.
- 43. Dahlan HM, Karsani SA, Rahman MA, Hamid NA, Top AG, et al. (2012) Proteomic analysis reveals that treatment with tocotrienols reverses the effect of H(2)O(2) exposure on peroxiredoxin expression in human lymphocytes from young and old individuals. J Nutr Biochem 23: 741–751.
- 44. Tan EC, Karsani SA, Foo GT, Wong SM, Rahman NA, et al. (2012) Proteomic analysis of cell suspension cultures of Boesenbergia rotunda induced by phenylalanine: identification of proteins involved in flavonoid and phenylpropanoid biosynthesis pathways. Plant Cell Tissue and Organ Culture 111: 219–229.