Proteome Profile of Swine Testicular Cells Infected with Porcine Transmissible Gastroenteritis Coronavirus

The interactions occurring between a virus and a host cell during a viral infection are complex. The purpose of this paper was to analyze altered cellular protein levels in porcine transmissible gastroenteritis coronavirus (TGEV)-infected swine testicular (ST) cells in order to determine potential virus-host interactions. A proteomic approach using isobaric tags for relative and absolute quantitation (iTRAQ)-coupled two-dimensional liquid chromatography-tandem mass spectrometry identification was conducted on the TGEV-infected ST cells. The results showed that the 4-plex iTRAQ-based quantitative approach identified 4,112 proteins, 146 of which showed significant changes in expression 48 h after infection. At 64 h post infection, 219 of these proteins showed significant change, further indicating that a larger number of proteomic changes appear to occur during the later stages of infection. Gene ontology analysis of the altered proteins showed enrichment in multiple biological processes, including cell adhesion, response to stress, generation of precursor metabolites and energy, cell motility, protein complex assembly, growth, developmental maturation, immune system process, extracellular matrix organization, locomotion, cell-cell signaling, neurological system process, and cell junction organization. Changes in the expression levels of transforming growth factor beta 1 (TGF-β1), caspase-8, and heat shock protein 90 alpha (HSP90α) were also verified by western blot analysis. To our knowledge, this study is the first time the response profile of ST host cells following TGEV infection has been analyzed using iTRAQ technology, and our description of the late proteomic changes that are occurring after the time of vigorous viral production are novel. Therefore, this study provides a solid foundation for further investigation, and will likely help us to better understand the mechanisms of TGEV infection and pathogenesis.


Introduction
Porcine transmissible gastroenteritis coronavirus (TGEV) is an animal coronavirus that causes severe gastroenteritis in young TGEV-seronegative pigs. Various breeds of pigs, regardless of age, are susceptible to TGEV; however, the mortality rate for piglets under 2 weeks of age is the highest, reaching almost 100%. Diseased pigs often present with vomiting, dehydration, and severe diarrhea. Further, the disease is known to affect pigs in many countries throughout the world and an outbreak can cause enormous losses in the pig industry [1,2]. The pathogen, TGEV, which belongs to the Alphacoronavirus genus of the Coronavirinae subfamily within the family Coronaviridae, is an enveloped, nonsegmented, single-stranded positive-sense RNA virus [3,4]. The envelop, core, and nucleocapsid of the TGEV virion contain four major structural proteins: the nucleocapsid (N) protein, the membrane (M) glycoprotein, the small envelope (E) protein, and the spike (S) protein [5]. The tropism and pathogenicity of the virus are influenced by the S protein, which has four major antigenic sites, A, B, C, and D, with site A being the major inducer of antibody neutralization [3,5]. The M protein, which plays a central role in virus assembly by interacting with viral ribonucleoprotein (RNP) and S glycoproteins [6], is embedded within the virus membrane and interacts with the nucleocapsid, forming the core of TGEV virion. In addition, the N-terminal domain of the M protein is essential for interferon alpha (IFN-a) induction [7], which is involved in the host's innate immune response. The E protein, a transmembrane protein that acts as a minor structural component in TGEV and affects virus morphogenesis, is essential for virion assembly and release [8].
TGEV RNA, along with the N protein, is infectious and invades the organism through the digestive and respiratory tracts, resulting in infection of the small intestinal enterocytes, villous atrophy, and severe watery diarrhea. These changes in intestinal health are known to be important during the pathogenesis of TGEV infection [9]. Furthermore, corresponding to these pathologic changes observed in vivo, TGEV can also propagate and cause cytopathic effects (CPEs) in multiple types of cultured cells, such as swine testicular (ST) cells, PK-15 cells, and villous enterocytes. Notably, ST cells are more susceptible to TGEV, and higher levels of virus replication have been observed in this cell line [10,11].
The full RNA genome of TGEV is approximately 28.5 kb in length and has a 59-cap structure and a poly(A) tail at the 39 end. The 9 open reading frame (ORF) genes included in the TGEV genome are arranged in the following order 59-la-lb-S-3a-3b-E-M-N-7-39. The first gene at the 59 end consists of two large ORFs, ORF la and ORF lb, which constitute the replicase gene, known for its RNA-dependent RNA-polymerase and helicase activities, as well as other enzymes, such as endoribonuclease, 39-59exoribonuclease, 29-O-ribose methyltransferase, ribose ADP 1'' phosphatase, etc. [12]. ORF2, ORF4, ORF5, and ORF6 encode the S, E, M, and N proteins, respectively, while ORF3a, ORF3b, and ORF7 encode non-structural proteins [13]. Some investigators have suggested that ORF3 may be related to viral virulence and pathogenesis [12], while ORF7 may interact with host cell proteins and play a role in TGEV replication [14]. In fact, a recent study indicates that plasmid-transcribed small hairpin (sh) RNAs targeting the ORF7 gene of TGEV is capable of inhibiting virus replication and expression of the viral target gene in ST cells in vitro [15]. Although we have some knowledge concerning the translation and function of these viral proteins, the interactions that occur between these proteins and host cell proteins are not fully understood.
Importantly, recent advances in proteomic technology have allowed for more in depth investigation of virus-host interactions, and different techniques have been successfully applied to identify altered proteins in infected host cells and tissues. For example, Sun et al. [16] have identified 35 differentially expressed proteins in PK-15 cells infected with classical swine fever virus (CSFV) using two-dimensional polyacrylamide gel electrophoresis (2D PAGE) followed by matrix-assisted laser desorption-ionization time-offlight tandem mass spectrometry (MALDI-TOF-MS/MS). In addition, two-dimensional fluorescence difference gel electrophoresis (2D-DIGE) and MS/MS proteomic approaches have been applied to characterize protein changes occurring in host cells in response to porcine circovirus type 2 (PCV2) infection [17]. The same methods have also been studied for many other pathogenic animal viruses, including porcine reproductive and respiratory syndrome virus (PRRSV) [18], coronavirus infectious bronchitis virus (IBV) [19], severe acute respiratory syndrome-associated coronavirus (SARS-CoV) [20], and TGEV [21]. However, these conventional approaches based on 2D gel electrophoresis are not suitable for detecting low abundance, hydrophobic, or very acidic/ basic proteins. On the other hand, the isobaric tags for relative and absolute quantitation (iTRAQ) technique, in association with liquid chromatograph (LC), is a more advanced method for proteomic research, and is capable of detecting a much larger number of proteins, even those with low abundance, in addition to identifying and quantifying the proteins simultaneously [22]. To this end, Lu et al. [23] previously used the iTRAQ method to identify 160 significantly altered proteins in pulmonary alveolar macrophages (PAMs) infected with PRRSV. Similarly, this method has been used to investigate influenza virus infection in primary human macrophages [24], human immunodeficiency virus 1 (HIV-1) infection in CD4 + T cells [25], and Epstein-Barr virus (EBV) infection in nasopharyngeal carcinoma cell line [26].
Here, we report the first differential proteomic analysis of TGEV-infected and uninfected ST cells using iTRAQ labeling followed by 2D-LC-MS and bioinformatic analyses. The proteomic data obtained in this study will help to enhance our understanding of the host response to TGEV infection, but also provide new insights on the mechanisms of disease onset.

Materials and Methods
Cell culture and viral replication ST cells were obtained from the American Type Culture Collection (ATCC). The cells were cultured in high-glucose Dulbecco's modified Eagle's medium (DMEM; GIBCO, UK) containing 1% L-glutamine and 10% fetal bovine serum (FBS) (Hyclone, Logan, UT) at 37uC in 5% CO 2 . Culture medium was replaced two to three times per week. The TGEV TH-98 strain was isolated from a suburb of Harbin, Heilongjiang province, China. The virus was propagated in ST cells and preserved at 2 70uC in our laboratory.

TGEV infection
The monolayer of confluent ST cells was dispersed with 0.25% trypsin and 0.02% ethylenediaminetetraacetic acid (EDTA) and seeded in 6-cm cell culture flasks. After a 24 h incubation period, the culture medium was removed and the ST cells were washed with phosphate buffered saline (PBS, pH 7.4). The cells were then infected with the TGEV TH-98 strain at a 50% tissue culture infectious dose (TCID 50

Reverse transcription polymerase chain reaction (RT-PCR) and real time quantitative PCR (qRT-PCR)
To determine the extent of TGEV infection, conventional RT-PCR and qRT-PCR assays were performed to detect the viral N gene. Monolayers of ST cells were infected with TGEV as described above. Cells were collected from 24 to 80 hpi at 8 h intervals, and the total RNA of the infected cells was extracted using Trizol (Invitrogen). RNA samples were reverse-transcribed using PrimeScript RT reagent Kit (Takara Bio, Dalian, China), according to the manufacturer's instructions. The RT reaction was incubated at 37uC for 15 min followed by 85uC for 5 s. A mixture of oligo dT primers and random 6 mers was used in the RT step. The cDNA was stored at 220uC until further use.
PCR was performed for the TGEV N gene in a 25 ml reaction mixture containing 1 ml of the cDNA, 0.5 ml of each forward (F) and reverse (R) primer, 12.5 ml of Premix Taq (Takara Bio, Dalian, China), and 10.5 ml DEPC water, starting with a 5 min denaturation at 95 C followed by 32 cycles of 30 s denaturation at 95 C, 30 s annealing at 56 C, and 40 s extension at 72 C. A final extension step was carried out at 72 C for 10 min. RT-PCR products were resolved on a 15 g/L agarose gel. The following PCR primers were used in this study: TGEV N (F, 59-GAGC-AGTGCCAAGCATTACCC-39 and R, 59-GACTTCTAT CT-GGTCGCCATCTTC-39) and b-actin (F, 59-GCAAGGACCTC-TACGCCAA-39 and R, 59-CTGGAAGGTGGACAGCGAG-39).
The mRNA expression level of the TGEV N gene was quantified using a SYBR Green assay on a Bio-Rad iQ5 real time PCR detection system as described previously [27]. We used the same primers listed above for qRT-PCR. Reactions were carried out in 50 ml volumes containing 0.5 ml of 20 6 SYBR Green I, 2 ml of cDNA template, 1 ml of each F and R primer, 25 ml of 2 6 PCR buffer, and 20.5 ml DEPC water. The cycling conditions were 94uC for 4 min, followed by 35 cycles of 94uC for 20 s, 60uC for 30 s, 72uC for 30 s, and then a final extension of 10 min at 72uC. The relative gene expression was determined with the 2 (2DDCt) method [28], and the tests were performed in triplicate. Following ST cell infection, cells were collected at 48 and 64 hpi  by centrifugation at 3,000 rpm for 5 min at 4uC, washed twice with PBS, and 1 mL of iTRAQ lysis solution (8 M urea, 1% (w/v) dithiothreitol (DTT)) containing protease inhibitor was added. Then, the cells were put in an ice bath and broken up by sonication. The solution was then mixed for 30 min at 4uC. The soluble protein fraction was harvested by centrifugation at 40,000 6 g for 30 min at 4uC and the debris was discarded. The protein concentration was determined with the Bradford protein assay (2-D Quant Kit, Bestbio, China). A 100 mg aliquot of protein from each sample was reduced, alkylated, and trypsin-digested as described in the iTRAQ protocol (AB Sciex, American), followed by labeling with the 4-plex iTRAQ Reagents Multiplex Kit according to the manufacturer's instructions (AB Sciex, American). Two virus-free samples at 48 h and 64 h were labeled with iTRAQ tags 114 and 115, while two TGEV-infected samples at 48 h and 64 h were labeled with tags 116 and 117. The labeled digests were then pooled, dried using a vacuum freeze drier (Christ RVC 2225, Germany), and preserved at 220uC for later use.

2D LC-MS/MS analysis
The combined peptide mixtures were separated by reversed phase high-performance liquid chromatography (HPLC) (Ekspert ultraLC 100, AB Sciex, USA) on a Durashell-C18 reverse phase column (4.6 mm 6 250 mm, 5 mm 100 Å , Agela). The mobile phases used were composed of 20 mM ammonium formate (pH 10) in water (labeled mobile phase A) and 20 mM ammonium formate (pH 10) in acetonitrile(ACN) (mobile phase B). The flow rate was 0.8 mL/min, and the elutant was collected into 48 centrifuge tubes at each minute after the first 5 min. Each aliquot was then dried by vacuum freezing.
The peptides were then analyzed with a nanoflow reversedphase liquid chromatography-tandem mass spectrometry (nano-RPLC-MS/MS) system (TripleTOF 5600, AB Sciex, USA). The above 48 tubes were merged into 10 components dissolved in 2% ACN and 0.1% formic acid (FA), then centrifuged at 12,000 6 g for 10 min. The supernatant (8 ml) was used for loading at a rate of 2 ml/min, with a separation rate of 0.3 ml/min. The mobile phase A used in this analysis was composed of 2% ACN and 0.2% FA, while mobile phase B was composed of 98% ACN and 0.1% FA. The following MS parameters were utilized: source gas parameters (ion spray voltage: 2.

Data analysis and bioinformatics
Protein identification and quantification were performed with the ProteinPilot software (version 4.0, AB Sciex) using the Paragon algorithm. Each MS/MS spectrum was searched against a database of Sus scrofa protein sequences (NCBI nr, released in March 2011, downloaded from ftp://ftp.ncbi.nih.gov/genomes/ Sus_scrofa/protein/). The following search parameters were used: iTRAQ 4-plex (peptide labeled), cysteine alkylation with methyl methanethiosulfonate(MMTS), trypsin digestion, biological modifications allowed, a thorough search, a detected protein threshold of 95% confidence (unused Protscore $1.3), and a critical false discovery rate (FDR) of 1%. The peptide and protein selection criteria for relative quantitation were performed as described previously, whereby only peptides unique for a given protein were considered [29]. In addition, proteins with an iTRAQ ratio higher than 20 or lower than 0.05 as well as proteins in reverse database were removed [30].
To assign enriched Gene Ontology (GO) terms to the identified proteins, the differentially expressed proteins identified from iTRAQ experiments and all of the 4,112 measured proteins were classified based on their GO annotations using QuickGO (http:// www.ebi.ac.uk/QuickGO/), with UniProt ID (http://www. uniprot.org/?tab=mapping) as the data source. GO enrichment analysis of the differentially regulated proteins was evaluated using all of the 4,112 quantified proteins as background with hypergeometric distribution [31]. Categories belonging to biological processes, molecular functions, and cellular components that were identified at a confidence level of 95% were included in the analysis. The protein-protein interaction network for a select group of proteins was analyzed using the STRING 9.1 database (http://string-db.org/). Network analysis was set at medium confidence (STRING score .0.4).

Western blot analysis
Following ST cell infection with TGEV, the culture medium was removed after incubating for 48 h and 64 h; then, the cells were washed with cold PBS and collected after centrifugation at 3,000 rpm for 10 min. Cells were then lysed in RIPA lysis buffer with protease inhibitors (Applygen Technologies Inc., China). Cellular debris was removed by centrifugation at 12,000 6 g for 5 min at 4uC, and the protein concentration was measured by Coomassie blue G250 staining. An equal amount (20 mg) of cell lysate from each sample was separated using 10% SDS-PAGE and then transferred to polyvinyl difluoride (PVDF) membranes (Millipore, Bedford, USA). The PVDF membranes were then blocked with 5% (w/v) de-fatted milk powder dissolved in tris buffered saline and tween 20 (TBST) buffer (150 mM NaCl, 50 mM Tris, 0.05% Tween 20) for 1 h at 37uC. After blocking, membranes were incubated with anti-glyceraldehyde 3-phosphate dehydrogenase (GAPDH) mouse monoclonal antibody (1:3000; Western Biotechnology, China), anti-heat shock protein 90 alpha (Hsp90a/HSP90AA1) antibody (1:300; Abcam, Cambridge, UK), anti-caspase 8 antibody (1:300; Abcam, Cambridge, UK), or antitransforming growth factor b 1 (TGF-b1/TGFB1) antibody (1:300; Abcam, Cambridge, UK) overnight at 4uC, followed by HRP-conjugated secondary antibody (1:5000; Western Biotechnology, China) for 1.5 h at 37uC. The membranes were then washed four times in TBST buffer for 5 min each time. Protein band detection was performed using ECL reagents (Applygen Technologies Inc., China), and the band intensities were analyzed using Labworks 4.6 software.

Confirmation of TGEV infection in ST Cells
After introducing TGEV into the ST cells, we observed the induction of typical CPEs, including cell rounding, swelling, granular degeneration of the cytoplasm, cell detachment, and severely diseased cell morphology, from 40 to 64 h after inoculation (

Dynamic changes in viral gene expression in infected cells
To further identify the extent of TGEV infection, the mRNA expression levels of viral genes in infected cells were determined using qRT-PCR. Comparative threshold (Ct) cycle values in three independent experiments were calculated and the results indicated that the average Ct value for the TGEV N gene ranged from 25.2 to 27.5. Correspondingly, the average Ct value observed for the bactin control gene ranged from 19.6 to 21.0. The relative expression of TGEV N mRNA was calculated using the 2 (-DDCT) method [28], and the change in expression at each time point is indicated in Figure 2B. These data show that, following infection, the viral mRNA levels increased gradually over time, and reached a peak at 48 hpi. Following this time point, the viral mRNA levels appear to decrease.

Protein identification by MS
In the infected ST cells, a total of 29,214 peptides and 4,364 proteins were detected (Table S1); however, only 4,112 proteins were quantified reliably (Table S2). Notably, the abnormal proteins, such as the proteins with iTRAQ ratio higher than 20 or lower than 0.05, which are not quantifiable [30], were removed and only proteins with reasonable ratios across all channels were investigated further. Figure 3A depicts the scatter plots for the log 10 116/114 and log 10 117/115 ratios in the iTRAQ experiment. Linear regression analysis showed that correlation (R 2 ) was 0.58, with a p-value less than 0.05. These results suggest that the alterations in protein abundance due to virus infection were nearlinear dependency between the two time points. In order to  identify the proteins that were significantly different at each time point (infected/uninfected) or between the different time points, we analyzed the distribution of ratios for the identified proteins as shown in the Figure 3B. For the distribution range of the differentially expressed proteins identified at 48 hpi, shown in Figure 3C, a ratio higher than 3.35 or lower than 21.35 was defined as a statistically significant difference in protein expression. At 64 hpi, a ratio higher than 4.55 or lower than 22.15 was defined as a statistically significant difference in protein expression. According to analyses, the differentially expressed proteins identified were considered to show a significant upward or downward trend if their expression ratios were greater than 4.0 or less than 0.25 compared to the control group.
Using the criterion listed above, the expression of 146 proteins was significantly changed at 48 hpi (95 upregulated and 51 downregulated), while 219 proteins were significantly changed at 64 hpi (172 upregulated and 47 downregulated). Further, 72 proteins were identified to be significantly different between the two time points (54 upregulated and 18 downregulated), resulting in a total of 316 unique proteins being significantly altered during TGEV infection, including 162 predicted proteins (Table S3 and  Table 1 (excluding the predicted proteins)). Because the current pig genome database is poorly annotated compared to the human genome database, there were numerous proteins that were unassigned or uncharacterized, resulting in a large number of predicted proteins in our analysis. However, our ability to detect the unannotated proteins by MS demonstrates that they do existence in this species, and additional research concerning their function is warranted.

GO enrichment analysis
Biological process-based enrichment analysis of the differentially expressed proteins revealed that six common GO terms were significantly enriched in this set of proteins (p,0.05). Thus, it appears that in TGEV-infected ST cells at 48 and 64 hpi there are expression changes in proteins that are related to cell adhesion, neurological system processes, extracellular matrix organization, locomotion, cell junction organization, and cell-cell signaling. Moreover, at the later time point, 64 hpi, our GO term analysis also indicated that a significant number of the differentially expressed proteins were related to cellular stress (p = 8.18E-4), generation of precursor metabolites and energy (p = 2.74E-3), cell motility (p = 6.71E-3), protein complex assembly (p = 4.69E-2), growth (p = 3.87E-2), developmental maturation (p = 1.53E-2), and immune system processes (p = 4.67E-2) ( Table 2).
To further investigate the localization pattern of these differentially expressed genes, a cellular component-based enrichment analysis was performed. At 48 hpi, we observed the significant enrichments in extracellular region (p = 1.29E-4), proteinaceous extracellular matrix (p = 1.62E-4), and extracellular  space (p = 1.52E-2) (Table S4). In addition, 37 differentially expressed proteins were also significantly enriched (p = 8.65E-3) in mitochondrion at 64 hpi (Table S5). The final step of our GO enrichment analysis consisted of investigating the mechanistic role these genes play in the cell. To do so, we performed a molecular function-based enrichment analysis. This analysis showed that two GO terms, unfolded protein binding (p = 2.67E-2) and transmembrane transporter activity (p = 3.55E-2), were significantly enriched at 64 hpi (Table  S5). Further GO analysis of the differentially expressed proteins between the two time points indicated that there were no significant enriched terms.

Protein-protein interaction analysis
In order to understand the interactions between TGEV and host cell proteins, we further analyzed the differentially expressed proteins by searching the STRING 9.1 database (http://string-db. org/) for protein-protein interactions (Figure 4). In this STRING analysis, the interactions (edges) of the submitted proteins (nodes) were scored according to known and predicted protein-protein interactions. We created three protein network maps: one for proteins changed significantly at 48 hpi (30 nodes and 15 edges; Figure 4A), one for proteins changed significantly at 64 hpi (66 nodes and 70 edges; Figure 4B), and one for the proteins that were significantly changed when the viral infection was prolonged from 48 to 64 h (24 nodes and 9 edges; Figure 4C). Notably, the protein network constructed for the 64 hpi time point is clearly much more extensive than the two other networks, and these proteinprotein interactions suggest the existence of reported functional linkages. GO enrichment analysis for the STRING protein network at 64 hpi showed that several biological processes were significantly affected (p,0.05 based on the FDR correction) in this network, including the regulation of viral genome replication, the innate immune response, negative regulation of viral genome replication, positive and negative regulation of viral processes, and ATP biosynthetic processes (Table 3). However, at 48 hpi, the most enriched biological process was related to cell recognition during phagocytosis(p = 8.02E-1). In Figure 4C, we have shown that the majority proteins in these protein networks, such as radical S-adenosyl methionine domain containing protein 2 (RSAD2), Mx dynamin-like GTPase 1 (Mx1), 29-59-oligoadenylate synthetase 1 (OAS1), Mx dynamin-like GTPase 2 (Mx2), are involved in the innate immune response. These data suggest that some entirely different host proteins, interactions, or processes, including the immune response, were perturbed at these times during TGEV infection. edges. Network analysis was set at medium confidence (STRING score = 0.4). Seven different colored lines were used to represent the types of evidence for the association: green, neighborhood evidence; red, gene fusion; blue, co-occurrence; black, co-expression; purple, experimental; light blue, database; yellow, text mining. doi:10.1371/journal.pone.0110647.g004 Western blot confirmation of altered expression for three of the differentially expressed proteins To further confirm the proteomic data for three of the proteins, western blot analysis was performed to investigate the changes in the expression of HSP90a, caspase 8, and TGF-b1. The proteins were selected based on three criteria: 1) the expression of the protein was increased or decreased during TGEV infection according to our proteomics data; 2) the protein is known to be relevant during viral infection; and 3) each protein analyzed needs to be involved in a special biological process as determined by our GO enrichment analysis [32]. HSP90a, caspase 8, and TGF-b1 all filled these criteria and their protein expression was analyzed via western blot analysis of the cell lysate. As shown in Figure 5, the expression of HSP90a was significantly downregulated in TGEVinfected cells at 64 hpi, while the expression of caspase-8 was upregulated from 48 to 64 hpi in these cells. The expression of TGF-b1 was also significantly induced in TGEV-infected cells following infection. Thus, these results confirm the altered expression observed in the proteomic data for these three representative proteins during TGEV infection.

Discussion
The interactions between a virus and a host cell during a viral infection are complex, involving numerous genes and signaling pathways. ST cells are known to be sensitive to TGEV, resulting in increased viral multiplication and CPEs [15]. In order to better understand the interactions between the host proteome and TGEV, we adopted an iTRAQ quantitative proteomic approach to investigate the altered cellular proteins of the ST cells during TGEV infection in vitro. Compared with the 2-DE and 2D-DIGE methods often used, the 2D-LC-MS/MS method utilized here provides more quantitative and qualitative information about the proteins, and can also detect membrane proteins, hydrophobic proteins, higher molecular weight proteins, and low-abundance proteins, which are often missed by other methods. iTRAQ also has more advantages compared to isotope-coded affinity tags (ICAT) and stable isotope labeling by amino acids in cell culture (SILAC) methods, which both allow multiple labeling and quantitation of four to eight samples simultaneously with high sensitivity [22,33,34]. Further, the iTRAQ technique has been widely used for quantitative proteomics, including protein expression analysis and biomarker identification [23][24][25][26]35].
Prior to proteomic analysis, we determined which time points to investigate following infection by observing the morphological changes and analyzing viral gene expression dynamics in the TGEV infected cells. The results indicated that TGEV induced significant CPEs from 40 to 64 hpi in infected cells compared to the mock infected cells. At 40 hpi, less than 50% of the infected cells were morphologically altered, while at 48 hpi more than 80% infected cells showed rounding and granular degeneration. Further, the mRNA level of the viral N gene in ST cells continuously increased in the infected cells until 48 h, at which time we observed the highest viral replication level. At 64 hpi, the morphological effects observed were much more pronounced, characterized by even more cellular rounding and detachment. However, the mRNA levels of the viral N gene decreased rapidly from 48 to 64 h, a phenomenon we believe may be attributed to the host's immune response or a decrease in infected cell viability as the TGEV infection progressed. Based on our qRT-PCR and CPE analyses, we choose to more deeply investigate the proteomic changes occurring in the TGEV-infected ST cells at 48 hpi and 64 hpi using a 4-plex iTRAQ analysis.
In our analysis, we observed a statistically significant change in the expression of 316 proteins during TGEV infection in vitro. This number includes protein changes that were unique for a specific time point as well as those shared at these different time conditions. For example, the expression level of HSP90a expression was unchanged at 48 hpi, but decreased at 64 hpi, Note: The significance of the GO biological process is derived from the network in Figure 4B and was determined using the FDR correction (p,0.05). doi:10.1371/journal.pone.0110647.t003 making this change unique for the latter time point. On the other hand, TGF-b1 was observed to increase at both of the time points, and was thus labeled a shared protein change. Moreover, the 316 altered proteins also includes proteins that changed from 48 hpi to 64 hpi, rather than one of these time points compared to noninfected cells. For example, mitochondrial aldehyde dehydrogenase 2 (ALDH2) and MHC class I antigen (PD1) were not changed at 48 or 64 hpi compared to the control group, but increased at 64 hpi compared with 48 hpi. We also observed a larger proteomic shift at 64 hpi compared to the 48 hpi time point in the infected ST cells. Further, some proteins previously reported to play a role in virus-induced host cell death, such as caspase-8, caspase-3, caspase-9, and porcine aminopeptidase-N (pAPN) [36][37][38], were also identified using this iTRAQ technique. These caspase proteins are known to be involved in TGEV-induced cell apoptosis processes, while pAPN is the cell receptor for TGEV. Our results indicate that TGEV infection caused significant upregulation of caspase-8 expression at two time points (approximately 7-fold at 48 hpi and 16-fold at 64 hpi) in the virus-infected ST cells, and this change was verified by western blotting analysis. However, the expression of caspase-3, caspase-9, and pAPN was not significantly altered, indicating that the pathways involving these genes are not altered or that other proteins are compensating for their lack of change. In this regard, we identified an additional 15 proteins involved in cell death pathways that had significantly altered expression levels (p = 4.46E-2) (Table S6), including melanoma differentiation associated protein-5 (MDA5), monocyte chemoattractant protein 1 (CCL2), thioredoxin-dependent peroxide reductase, mitochondrial (PRDX3), peroxiredoxin-2 (PRDX2), predicted protein CYR61 (CYR61), keratin, type II cytoskeletal 8 (KRT8), predicted bcl-2-like protein 13 (BCL2L13), predicted integrin alpha-5 isoform 1 (ITGA5), TGF-b1, amyloid beta A4 protein (APP), clusterin (CLU), C-C motif chemokine 5 (CCL5), heat shock 70 kDa protein 1B (HSPA1B), alpha-crystallin B chain (CRYAB), voltage-dependent anion-selective channel protein 1 (VDAC1), all of which, with the exception of PRDX2 and BCL2L13 were upregulated at one or two time points. Regulation of cell death is known to be important for replication and pathogenesis in various coronaviruses [39], and we believe that further research on these proteins will lead to a better understanding of cell death regulation during TGEV infection.
In order to determine what other processes, in addition to cell death, were affected by TGEV infection, we performed a GO enrichment analysis for the different temporal conditions. This analysis indicated that six biological processes were significantly affected at 48 and 64 hpi, and the differentially expressed proteins involved in these processes were almost the same. The large overlap between the two time points suggests that some of the same sets of host proteins or processes were disturbed at these times. However, it is also likely that some processes were affected solely at one time point or the other. At 48 hpi, serine/threonineprotein phosphatase PP1-beta-catalytic subunit (PPP1CB), scavenger receptor class B member 1 (SCARB1), transforming growth factor-beta-induced protein ig-h3 (TGFBI), and predicted inositol 1,4,5-trisphosphate receptor type 3 (ITPR3) were uniquely altered, likely indicating changes in cell adhesion and/or cell-cell signaling processes. At 64 hpi, on the other hand, calreticulin (CALR), predicted tumor-associated calcium signal transducer 2-like (TACSTD2), vascular cell adhesion molecule, galectin-3 (LGALS3), glutamate dehydrogenase 1 (GLUD1), and C-X-C motif chemokine 16 (CXCL16) were uniquely changed, also indicating changes in cell adhesion and/or cell-cell signaling as well as extracellular matrix organization and locomotion. We believe that these uniquely altered proteins reflect changes in specific/specialized processes at each time point that are tightly linked to the temporal changes observed in the host cell morphology and gene/protein expression after TGEV infection.
The most significantly enriched GO category related to the differentially expressed proteins was stress, which included 12 differentially expressed proteins at 48 hpi and 27 different proteins at 64 hpi. The increased number of proteins association with this GO term at 48 hpi likely highlights the initial upregulation of the cellular stress response, while the higher number at 64 hpi indicates that the stress response to TGEV infection is likely more fully induced at this later stage. HSPs, also known as stress proteins, are often involved in the cellular response to stress, influencing changes in the state or activity of the cell or organism. HSP90, which has two isoforms (HSP90a and HSP90b), is one of the most abundant molecular chaperones that is induced in response to cellular stress, and it functions to stabilize proteins involved in cell growth and anti-apoptotic signaling [40]. The expression of HSP90a has been reported to play an important role in the replication of some viruses, such as Ebola virus (EBOV) [41], hepatitis C virus (HCV) [42], influenza virus [43], and Japanese encephalitis virus [44]. On the other hand, the reduction of HSP90b has been reported to decrease the correct assembly of human enterovirus 71 viral particles [40]. In this study, HSP90a and heat shock 90kD protein 1, beta (HSPCB/HSP90b) were significantly downregulated at 64 hpi in the TGEV-infected ST cells, but were unchanged at 48 hpi, indicating that they may play a similar role in TGEV infection. Interestingly, a member of the HSP70 protein family, heat shock 70 kDa protein 1B (HSPA1B), as well as mitochondrial 60 kDa heat shock protein (HSP60) were both upregulated in infected ST cells at 48 and/or 64 hpi. HSP60 is a mitochondrial chaperonin protein involved in protein folding and a number of extracellular immunomodulatory activities. Elevated expression of HSP60 is associated with a number of inflammatory disorders [45]. HSP70 plays an important role in multiple processes within cells, including protein translation, folding, intracellular trafficking, and degradation. A previous study has revealed that HSP70 is involved in all steps of the viral life cycle, including replication, and is highly specific in regards to viral response, differing from one cell to another for any given virus type [46]. For example, silencing HSP70 expression has been associated with an increase in viral protein levels, while an increase in HSP70 has been suspected to be the initial cellular response to protect against viral infection in rotavirus-infected cells [47]. Further, a recent study showed that HSP70 is an essential host factor for the replication of PRRSV as the silence of HSP70 significantly reduced PRRSV replication [48]. Our results provide new experimental evidence relating the expression of HSP90, HSP70, and HSP60 to TGEV infection, and we speculate that these proteins play a potential role in TGEV replication. Additional work is required to investigate the detailed role of these proteins during TGEV infection.
Furthermore, another significantly enriched GO process we observed that 11 significantly altered proteins was immune system processes. Most of these proteins were significantly upregulated at 64 hpi in response to the viral infection, while some were first upregulated at 48 hpi, including CCL5 and TGF-b1. Chemokines, such as CCL2, CCL5, and CXCL16, whose main function is macrophage recruitment and activation, are potentially involved in host-mediated immunopathology. A recent study showed that coronavirus infection of transgenic mice expressing CCL2 led to a dysregulated immune response without effective virus clearance and enhanced death [49]. In additional, TGEV-infection can induce the expression of proinflammatory genes, including CCL2, CCL5, and probable ATP-dependent RNA helicase DDX58 (DDX58/RIG-1), in cell culture and in vivo in the absence of viral protein 7 [50]. In this study, we observed an upregulation of CCL2, CCL5, CXCL16, TGF-b1, and DDX58 expression. TGF-b1 is a multifunctional cytokine, secreted from various cells, and, in immunology, it regulates cellular proliferation, differentiation, and other cellular functions for a variety of cell types, especially regulatory T cells [51]. Some research has indicated that SARS-CoV papain-like protease (PLpro) increases TGF-b1 mRNA expression and protein production in human promonocytes [52]. Further, Gomez-Laguna et al. [53] inferred that the upregulation of the TGF-b may impair the host immune response during PRRSV infection by limiting the overproduction of proinflammatory cytokines necessary to decrease PRRSV replication. In response to viral infection, DDX58 plays important roles in the recognition of RNA viruses in various cells, and has been identified as a candidate for a cytoplasmic viral dsRNA receptor [54]. Further, upregulation of this gene activates cells to produce type I interferons, which may increase the antiviral status of cells to protect against viral infection. In this regard, we found that interferon-inducible antiviral proteins, RSAD2, OAS1, were also upregulated in the period of late infection, suggesting that many of the proteins identified in this study are associated with inflammation, IFN activation, and the innate immune response. Increased expression of these proteins may help the virus enter the cell as well as potentially enhance TGEV replication or the host response against the virus, during the late stages of infection.
In conclusion, we used the iTRAQ method to identify 316 significantly altered proteins in TGEV-infected ST cells. A larger number of protein expression changes occurred at 64 hpi compared to 48 hpi, indicating a larger shift in the proteome in the later stages of infection. GO analysis of these differentially expressed proteins indicated that a number of diverse biological processes are affected. In addition, many of the significant immune response related changes in protein expression we discovered are novel and, to our knowledge, have not been detected in previous proteome study. Results from this study complement the previous proteomics data obtained concerning the host response to a viral infection, and further facilitates a better understanding of the pathogenic mechanisms of TGEV infection and molecular responses of host cells to this virus.