Differential Proteome Analysis of Chikungunya Virus Infection on Host Cells

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.


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 . 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.

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 postinfection (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.

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 TM 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.

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).

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).

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 21.42 and 1.72 CDK1 and PDHA1 respectively ( Figures 5C and 5D), which was comparable to the observed 21.77 and 1.96 fold-changes in 2-DGE analysis.

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 23.1 to 23.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 manipula- tion 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, ubiquitinproteasome 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 posttranscriptional 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]. CHIKVinduced 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.
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  The mean % spot volume (n = 5) was used for the analysis of fold difference between mock control and CHIKV-infected protein spots. SD represents standard deviation of five biological replicates. 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 cyclinactivated 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 downregulated, 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 cyclindependent 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.

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 FITCconjugated 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 280 uC. 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.

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.5610 5 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 uC. Thereafter, the cells were washed with PBS and incubated in FITC-conjugated secondary antibody (1:1000 dilution) for 30 min at 37 uC. 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 uC. 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 uC. 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,0006g 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).

Differential gel analysis
Gels were scanned with ImageScanner TM III (GE Healthcare, Uppsala, Sweden) and analyzed using ImageMaster TM 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 uC 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 ml of 10 ng/ml trypsin at 37 uC. 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 ml elution solution (50% ACN/0.1% FA) and mixed with saturated a-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 ml aliquots in duplicates. Mass spectra for each peptide were obtained on a MALDI-TOF/TOF (ABI 4800 Plus, Applied Biosystems TM , 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 TM , 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 mg) 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 uC 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 TM 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 Gene-Quant TM 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 mg of high quality RNA was converted to cDNA using High Capacity RNA-to-cDNA Kit (Applied Biosystems TM , 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 SYBRHGreen Master Mix (Applied Biosystems TM , Foster City, CA). Real-time qPCR was performed using StepOnePlus TM Real-Time PCR System (Applied Biosystems TM , 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. Figure S1 The proteome maps of differentially expressed whole cell proteins in mock control and CHIKV-infected WRL-68 cells.