Label‐Free Quantitative Proteomic Analysis of Differentially Expressed Membrane Proteins of Pulmonary Alveolar Macrophages Infected with Highly Pathogenic Porcine Reproductive and Respiratory Syndrome Virus and Its Attenuated Strain

Abstract Significant differences exist between the highly pathogenic (HP) porcine reproductive and respiratory syndrome virus (PRRSV) and its attenuated pathogenic (AP) strain in the ability to infect host cells. The mechanisms by which different virulent strains invade host cells remain relatively unknown. In this study, pulmonary alveolar macrophages (PAMs) are infected with HP‐PRRSV (HuN4) and AP‐PRRSV (HuN4‐F112) for 24 h, then harvested and subjected to label‐free quantitative MS. A total of 2849 proteins are identified, including 95 that are differentially expressed. Among them, 26 proteins are located on the membrane. The most differentially expressed proteins are involved in response to stimulus, metabolic process, and immune system process, which mainly have the function of binding and catalytic activity. Cluster of differentiation CD163, vimentin (VIM), and nmII as well as detected proteins are assessed together by string analysis, which elucidated a potentially different infection mechanism. According to the function annotations, PRRSV with different virulence may mainly differ in immunology, inflammation, immune evasion as well as cell apoptosis. This is the first attempt to explore the differential characteristics between HP‐PRRSV and its attenuated PRRSV infected PAMs focusing on membrane proteins which will be of great help to further understand the different infective mechanisms of HP‐PRRSV and AP‐PRRSV.


Introduction
Porcine reproductive and respiratory syndrome virus (PRRSV) has been a leading economically significant viral pathogen of swine worldwide for almost 28 years. [1][2][3] PRRSV, equine arteritis virus, simian hemorrhagic fever virus, and lactate dehydrogenaseelevating virus are members of the family Arteriviridae. [4,5] PRRSV is a positive-sense, single-stranded, RNA virus with a full-length genome of 15 kb that has a 5 cap and a 3 poly (A) tail. [4,[6][7][8] PRRSV was first reported in the United States in the late 1980s. [9] In 2006, several large-scale, severe outbreaks of highly pathogenic (HP) atypical PRRSV (HP-PRRSV) were reported in China and neighboring Asian countries. [2,3,10] Reported HP-PRRSV morbidity and mortality rates were much higher than previous pandemic PRRSV strains and associated with more severe clinical presentations and higher rectal temperature (>41°C). [11] The emergence of HP-PRRSV has caused great economic loss to the swine industry in China and made preventing and controlling PRRSV outbreaks difficult. Therefore, elucidating the causes of the greater virulence of PRRSV and the differences between the HP and attenuated pathogenic (AP) strains has become even more important. To this end, several trials have been conducted to identify virulence factors; these studies have resulted in some successes. [12,13] However, changes in virulence and pathogenic mechanisms are difficult to discern. Other than virulence in vivo, many distinctions in the biological aspects of HP-PRRSV and AP-PRRSV have been noted, such as in viral binding and entry into pulmonary alveolar macrophages (PAMs). PRRSV exhibits highly restricted cell tropism both in vivo and in vitro. [14] The virus can be detected only in well-differentiated macrophages of lungs, lymph nodes, Peyer's patches, spleen, tonsils, and thymus. PAMs are the main target cells of PRRSV. [15] PRRSV can also replicate in vitro in the African www.advancedsciencenews.com www.proteomics-journal.com

Significance of the study
Membrane proteins (MPs) of PAMs infected by highly proteomic (HP)-and attenuated proteomic (AP)-PRRSV have been elucidated by LC-MS/MS for label-free quantitative proteomics. Ninety-five differentially expressed proteins were identified and characterized. The most significant difference in the biological process between PAMs infected with HP-and AP-PRRSV is the metabolic process. Most different molecular functions were classified as binding and catalytic activities. Cellular component categories showed that 26 differentially expressed proteins were confirmed as MPs based on the annotation of UniProt database, such as RAP2A, VCL (Vinculin), IFITM3 function in cell-cell junctions, ERK signaling pathway, G protein signaling pathways, biotic stimulus, and so on. Among them, VCL is a kind of F-actin-binding protein which is involved in cell-matrix adhesion and cell-cell adhesion in humans. It was demonstrated that over expression of VCL could inhibit the replication of both HP-PRRSV and the attenuated PRRSV in the mRNA level. There were obvious differences in the inhibiting ability for HP-PRRSV and its attenuated strain. This is the first attempt to explore the differential characteristics between HP-PRRSV and its attenuated PRRSV-infected PAMs focusing on MPs which will be of great help to further understand the different infective mechanisms of HP-PRRSV and AP-PRRSV.
green monkey kidney cell line MA-104 and its derivatives, MARC-145 and CL-2621, which are considered permissive cell lines for PRRSV. [16,17] Reportedly, PRRSV targets cellular membrane proteins (MPs) and enters target cells through receptormediated endocytosis during viral infection. [18] Studies have investigated possible mechanisms employed by PRRSV to infect PAMs. [19] Elucidating the characteristics of viruses and interactions between viruses and host cells is increasingly important. However, exploring individual proteins within the many proteins in cells is difficult. Nonetheless, gradual advancements have come through use of proteomic techniques. Of these, MS-based quantitative proteomic techniques offer the advantage of better accuracy and sensitivity, and have been widely used to analyze host cell responses to viral infection. Among them, liquid chromatographytandem MS (LC-MS/MS) for label-free quantitative proteomics (LFQP) is an important mass spectrometric tool to detect and quantify large amounts of proteins. [20] Compared with quantitative proteomics using stable isotope labeling such as stable isotope labeling by amino acids in cell culture and isobaric tags for relative and absolute quantitation, LFQP detects greater amounts of proteins and important signaling pathways and networks. The aim of this study was to determine potentially different infection mechanisms used by the HP strain vHuN4 [10] and its derivative attenuated strain vHuN4-F112 [21,22] using LFQP.

Ethics Statement
The animal study protocols were approved by the Animal Care and Use Committee of Shanghai Veterinary Research Institute, Chinese Academy of Agricultural Sciences.

Virus Inoculation
After PAMs were washed with PBS three times, dead and nonadherent cells were removed when confluency exceeded 95%. Three dishes were inoculated with vHuN4 and another three with vHuN4-F112 at multiplicity of infection 1. An additional three dishes were inoculated with DMEM (Gibco-BRL) as a blank control. All dishes were incubated at 37°C in an atmosphere of 5% CO 2 , as described previously. [13] After incubation for 1 h, inocula were discarded and PAMs were washed with PBS three times. Cell monolayers in all dishes were overlaid with RPMI-1640 medium containing 2% fetal bovine serum and incubated at 37°C in a 5% CO 2 atmosphere for 24 h.

Extraction of MPs of PAMs
PAMs were digested with 2 mL 0.25% trypsinethylenediaminetetraacetic acid solution (Gibco/BRL), collected by gently pipetting, centrifuged at 1000 × g for 5 min and lysed using the ProteoExtract Transmembrane Protein Extraction Kits (NOVAGEN, EMD Biosciences, Inc., Madison, WI, USA), [24,25] according to the manufacturer's instructions. Cells were resuspended in Extraction Buffer 1 and protease inhibitor cocktail, incubated for 10 min at 4°C with gentle agitation, and centrifuged at 1000 × g for 5 min at 4°C. After removing supernatants, pellets were resuspended in 0.2 mL Extraction Buffer 2A and protease inhibition cocktail, incubated for 45 min at room temperature with gentle agitation, and centrifuged at www.advancedsciencenews.com www.proteomics-journal.com 16 000 × g for 15 min at 4°C. Supernatants were precipitated with 1 mL acetone and centrifuged at 12 000 × g for 10 min. After evaporating to dryness, 150 μL SDT buffer (4% sodium dodecyl sulfate, 100 mm Tris/HCl at pH 7.6, 0.1 m dithiothreitol) was added and mixtures were heated in boiling water for 5 min. After centrifugation, supernatants were collected and quantified with a BCA Protein Assay Kit (Bio-Rad, USA). Fig. 1

Protein Digestion
Digestion of protein (250 μg for each sample) was performed according to the FASP (Filter-Aided Sample Preparation) procedure. Briefly, the detergent, DTT and other low-molecular-weight components were removed using 200 μL UA buffer (8 m Urea, 150 mm Tris-HCl pH 8.0) by repeated ultrafiltration (Microcon units, 10 kD) facilitated by centrifugation. Then 100 μL 0.05 m iodoacetamide in UA buffer was added to block reduced cysteine residues and the samples were incubated for 20 min in darkness. The filter was washed with 100 μL UA buffer three times and then 100 μL 25 mm NH 4 HCO 3 twice. Finally, the protein suspension was digested with 3 μg trypsin (Promega) in 40 μL 25 mm NH 4 HCO 3 overnight at 37°C, and the resulting peptides were collected as a filtrate. The peptide content was estimated by UV light spectral density at 280 nm using an extinctions coefficient of 1.1 of 0.1% (g L −1 ) solution that was calculated on the basis of the frequency of tryptophan and tyrosine in vertebrate proteins. [26]

LC-MS/MS Analysis
The peptide of each sample was desalted on C18 Cartridges (Empore SPE Cartridges C18 (standard density), bed id 7 mm, volume 3 mL, Sigma), then concentrated by vacuum centrifugation and reconstituted in 40 μL of 0.1% (v/v) trifluoroacetic acid. MS experiments were performed on a Q Exactive mass spectrometer that was coupled to Easy nLC (Proxeon Biosystems, now Thermo Fisher Scientific). Five microgram peptide was loaded onto a C18-reversed phase column (Thermo Scientific Easy Column, 10 cm long, 75 μm inner diameter, 3 μm resin) in buffer A (2% acetonitrile and 0.1% Formic acid) and separated with a linear gradient of buffer B (80% acetonitrile and 0.1% Formic acid) at a flow rate of 250 nL min −1 controlled by IntelliFlow technology over 120 min. MS data was acquired using a datadependent top ten method dynamically choosing the most abundant precursor ions from the survey scan (300-1800 m/z) for HCD fragmentation. Determination of the target value is based on predictive automatic gain control. Dynamic exclusion duration was 25 s. Survey scans were acquired at a resolution of 70 000 at m/z 200 and resolution for HCD spectra was set to 17 500 at m/z 200. Normalized collision energy was 30 eV and the underfill ratio, which specifies the minimum percentage of the target value likely to be reached at maximum fill time, was defined as 0.1%. The instrument was run with peptide recognition mode enabled. MS experiments were performed triply for each sample. [27]

Sequence Database Searching and Data Analysis
The MS data were analyzed using MaxQuant software version 1.3.0.5. MS data were searched against the uniprot Sus scrofa sequence database (including 34 253 sequences downloaded on 12/27/2014). An initial search was set at a precursor mass window of 6 ppm. The search followed an enzymatic cleavage rule of Trypsin/P and allowed maximal two missed cleavage sites and a mass tolerance of 20 ppm for fragment ions. Carbamidomethylation of cysteines was defined as fixed modification, while protein N-terminal acetylation and methionine oxidation were defined as variable modifications for database searching. The cutoff of global false discovery rate for peptide and protein identification was set to 0.01. Label-free quantification was carried out in MaxQuant as previously described. Protein abundance was calculated on the basis of the normalized spectral protein intensity (LFQ intensity).
All statistical analyses were performed using unpaired t-tests. A p-value <0.05 and ratio >2 or <0.5 were considered to indicate significant differences. Gene Ontology (GO) annotation and functional classification of identified proteins was with Blast2GO ver. V2.6.2 with the current public database b2g_aug 12 (www.Blast2go.com). Identified proteins were classified www.advancedsciencenews.com www.proteomics-journal.com using Blast2go steps under default parameters: blast, mapping, and annotation. Protein-protein interaction networks were analyzed using (String software string-db.org/) . Confidence view was assigned a score of 0.4, indicating medium confidence.

Expression Vector Construction and Transfection
The porcine vinculin (VCL) were amplified from the cDNA obtained from PAM cells. Restriction enzyme sites were incorporated into the primer sequences to facilitate molecular cloning. PCR products were cloned into the pCAGGS vector to produce the porcine VCL expression vector. For transfection, cells were seeded in 6-well plates (Corning) and transfected at 70-80% confluency with respective constructed plasmids DNA by using Lipofectamine 3000 (Life Technologies), according to the manufacturer's instructions. The HP-HuN4 or HuN4-F112 was infected in MARC-145 cells at 36 h post transfection. Empty vector transfection samples served as controls in the experiment.

RNA Extraction, cDNA Synthesis, and Real-time RT-PCR
At 48 h post infection, total RNAs of PRRSV-infected cells were extracted using the RNeasy mini kit (Qiagen), the viral RNAs in supernatants were isolated using QIAamp Viral RNA Mini kit (Qiagen), according to the instruction manual. All the isolated RNAs were used as the template for synthesis of firststrand cDNA by RT-PCR using RT primed by oligo (dT) 18 primer using the PrimeScript RT Master Mix (Perfect Real Time, TaKaRa), according to the manufacturer's instructions. Then the cDNA templates were quantified using PRRSV-specific real-time RT-qPCR. [28] 3. Results

Label-Free LC-MS/MS for Quantitative Analysis of MPs of PAMs Infected with PRRSVs of Different Virulence
A total of 2849 proteins were detected by LFQP and are displayed in a heatmap (Figure 2). Statistical significance was determined using unpaired t-tests. For all tests, a p-value of <0.05 and ratio of >2 or <0.5 was considered to indicate a significant difference. Glyceraldehyde 3-phosphate dehydrogenase was the internal normalization control and the ratio of glyceraldehyde 3-phosphate dehydrogenase between HP and AP-PRRSV was 0.95 ± 0.55. The "only one exists" group indicated that proteins were detected only in one group but not another group due to low expression level ( Table 1). Correlation analysis indicated good repeatability of the technology (Figure 3).
A total of 95 differentially expressed proteins were identified ( Table 2). A control group was used to exclude false-positive interference. Data from the control group was also used to elucidate functions of target proteins. Differentially expressed proteins between control and HP-PRRSV-infected cells (Con/HP) and the www.advancedsciencenews.com www.proteomics-journal.com p <0.01 and ratio >2 or <0.5 indicated quantitative difference between two groups. The "only one exits" group indicated that proteins were detected three times in one, but not in the other group. HP-PRRSV group means proteins identified in the HP-PRRSV infected PAMs, while AP-PRRSV was proteins detected in the attenuated pathogenic PRRSV infected PAMs. Control was displayed as the mock.
control and AP-PRRSV-infected cells (Con/AP) are in Supporting Information.

Subcellular and Functional Characterization and Bioinformatics of Differently Expressed Proteins between the HP-PRRSV and AP-PRRSV Groups
To extend the molecular characterization of quantitative differences and only-one-exists groups, UniProt and GO databases were used to characterize information about biological processes (BPs), molecular functions (MF), and cellular components (CC).
BPs of H (vH/vA >2, including 41 proteins) and A (vH/vA < 0.5, including 54 proteins) groups are in Figures 4A and B BP. In group H, GO annotations were primarily distributed in response to stimulus (70.7%), metabolic process (68.3%), and immune system process (43.9%). Ratios were 63.0% response to stimulus, 90.7% metabolic process, and 33.3% immune system process in group A. The most significant difference in BP between PAMs infected with HP-PRRSV and AP-PRRSV was seen for metabolic process, which may be the major reason for the large differences among animals challenged with different virulence of PRRSV. Molecular function categories of H group and A group were shown in Figure 4A,B MF. Most different molecular functions were classified as binding (87.8 and 90.7%) and catalytic activity (31.7 and 44.4%). Binding of H group included enzyme binding (31.7%), nucleic acid binding (31.7%) and protein complex binding (26.8%), while group A mainly involved nucleic acid binding (31.5%), nucleoside phosphate binding (27.8%) and nucleotide binding (27.8%) from the analysis of GO distribution by level 4. Enzyme code distribution suggested there were five transferases (12.2%), one hydrolase (2.4%), one lyase (2.4%), and two ligases (4.9%) detected in group H. while four oxidoreductases (1.9%), five transferases (9.3%), five hydrolases (9.3%), three lyases (5.6%), and three isomerases (5.6%) in group A.
CC categories were illustrated in Figure 4A B CC. Ninety-five detected differentially expressed proteins were annotated and categorized to CCs of macromolecular complex (58.9%), membrane (57.9%), membrane-enclosed lumen (57.9%), and extracellular region (44.2%). Because of technological problems, we were unable to conclude that all the detected proteins were indeed MPs. However, 26 were confirmed based on the annotation of UniProt database and there also existed proteins partially anchored on the membrane or binding with the MPs.

Validation of Differentially Expressed MPs by Western Blots
To verify the differentially expressed proteins via LC-MS/MS for LFQP, Western blots were conducted for two proteins partially located on membranes. Expression of HSP70 and the KDEL receptor from cell lysates of vHuN4-infected and vHuN4-F112-infected PAMs, and DMEM-inoculated PAMs were tested with antibodies to the proteins. LFQP showed that the ratios between vHuN4infected and vHuN4-F112-infected PAMs reached 0.67 and 0.51, respectively. Western blots confirmed LFQP results ( Figure 5).

VCL Protein Inhibits Virus Replication
To examine whether the differentially expressed proteins detected affects virus infection, the VCL transient overexpression vector was transfected into MARC-145 cells, followed by the HP HuN4 or its attenuated strain infection. The PRRSV-specific RT-qPCR results showed that VCL protein could inhibit both viruses, especially for HP-HuN4 strain. Compared with the empty vector

Discussion
In this study, LFQP of MPs of HP-PRRSV-, AP-PRRSV-infected PAMs and the control was performed. Important information about target proteins related to virus infection was obtained. The HP-PRRSV strain vHuN4 and its derivative, the serially cellpassaged attenuated strain vHuN4-F112, [29] revealed different infection mechanisms in PAMs. A total of 2849 proteins were identified among the control, AP-PRRSV-infected, and HP-PRRSV-infected PAMs. Of these, 2400 were detected in all three groups (Figure 7). We focused on 95 differentially expressed proteins of AP-PRRSV-and HP-PRRSV-infected PAMs; among these, 43 were detected in only one group.

Identification and Potential Function of MPs in Immunology and Inflammation
PRRSV has been a threat to the global pig economy for several years because of its persistent infection, immune escape, and high mortality from inflammation and high fever. [30] The attenuated PRRSV vHuN4-F112 vaccine strain attenuated from HP-PRRSV vHuN4 by serial passages, which is now used in China. [31]  We analyzed MPs to identify factors associated with immunological effects and determine differences between HP-PRRSV and AP-PRRSV. MPs classified based on GO analysis are in Table 3.
A heatmap based on the UniProt database was constructed to comprehend the functions and BPs of differently expressed proteins ( Figure 8B), which revealed clustering and abundance of the 26 detected proteins that existed confidentially on the membrane based on UniProt database annotation. We first focused on proteins associated with immunology and inflammation with higher abundance in AP-PRRSV-infected PAMs. PTPN2 (vH/vA = 0.419), a member of the protein tyrosine phosphatase family, functions as signaling molecule that regulates cellular processes related to the Jak-STAT, IL-3, IL-5, and granulocyte-macrophage colony-stimulating factor (GM-CSF) signaling pathways; dephosphorylation of nonreceptor kinases including JAK1, JAK3, STAT3, and STAT6; and negative regulation of IL-2-, IL-4, IL-6, and IFN-mediated signaling. PTPN2 also functions in the response to inflammation via NF-κB. [32][33][34] SIPA1 (H-/A+) is a mitogen-induced GT-Pase activating protein for Ras-related regulatory proteins. It was related to G-protein signaling and blood-brain barrier and immune cell transmigration: VCAM01/CD106 (cluster of differentiation) signaling pathways. [35,36] C5AR1 (H-/A+) is the receptor for the chemotactic and inflammatory peptide anaphylatoxin C5a, which stimulates chemotaxis, granule enzyme release, intracellular calcium release, and superoxide anion production, participates in the innate and adaptive immune responses to www.advancedsciencenews.com www.proteomics-journal.com Table 3. Statistics analysis of proteins that existed definitely on the membrane. the lectin-induced complement pathway. [37][38][39][40][41][42] Polyribonucleotide 5 -hydroxyl-kinase Clp1 (CLP1) (vH/vA = 0.136) mainly acts as a kinase that binds ATP hosting 5 -hydroxyl-kinase activity, and functions in mRNA cleavage and in siRNA loading onto the RNAinduced silencing complex involved in RNA interference and its destruction. The kinase hClp1 phosphorylates and licenses synthetic siRNAs to assemble into an RNA-induced silencing complex for cleavage of target RNA. [43] Expression of these molecules was higher in the AP-PRRSV-infected group and had a similar level to the control group. Therefore, we concluded that they might be related to the immune escape of HP-PRRSV and that the attenuated vaccine strain would not have side effects on these immune factors expression in the host cell, and these proteins above would contribute to the resistance of PRRSV. Some proteins related to immunology and inflammation had higher expression in HP-PRRSV-infected cells. Raftlin (vH/vA = 2.554) protein is pivotal for maintenance of lipid rafts and may be involved in regulation of B-cell antigen receptor-mediated signaling. Raftlin promotes binding of double-stranded RNA, activations of B cell receptors and toll-like receptor 3 signaling pathways, and is involved in IL-17 production to release proinflammatory cytokines. Hence, the higher abundance of raftlin in the HP-PRRSV group compared to the AP-PRRSV and control groups may explain the more severe inflammation triggered by HP-PRRSV. [44,45] IARS (H+/A-), Isoleucyl-tRNA synthetase is a target of autoantibodies in autoimmune diseases. [46,47] VCL (H+/A-), Vinculin is a cytoskeletal protein associated with cellcell and cell-matrix junctions, and is related to the IL-3, IL-5, and GM-CSF signaling pathways. [48] RAP2A (H+/A-) is a member of the Ras oncogene family, small GTP-binding protein, of which active form interacts with several effectors [49][50][51] related to the ERK and G protein signaling pathways. [52,53] It was reported that PRRSV could induces prostaglandin E2 production through cyclooxygenase 1 and is related to ERK signaling. [54] The differential expression of RAP2A may provide information for future research on the interaction between PRRSV and the ERK signaling pathway. The abundance of RAP2A was lower in AP-PRRSV-infected PAMs than in HP-PRRSV-infected PAMs and control cells, which had similar levels. Thus, we hypothesized that RAP2A was associated with immunization with the attenuated vaccine. IFN-induced transmembrane protein 3, IFITM3 (vH/vA = 2.861), responds to biotic stimulus and had the highest expression levels among these proteins in the HP-PRRSV group, as compared to the proteins identified in the AP-PRRSV group but it was not detected in the control group. In humans, this protein negatively regulates the entry of multiple viruses, including influenza A virus, SARS coronavirus, Marburg virus (MARV), Ebola virus, Dengue virus, West Nile virus (WNV), human immunodeficiency virus (HIV) type 1, and vesicular stomatitis virus (VSV), into host cells. [55,56] It could inhibit hemagglutinin protein -mediated entry of influenza virus, GP1, two-mediated viral entry of Marburg virus and Ebola virus, S protein-mediated viral entry of SARS coronavirus, and G protein-mediated viral entry of VSV, thereby playing a critical role in the structural stability and function of vacuolar ATPase (v-ATPase). Establishing physical contact with the v-ATPase of endosomes is critical for proper clathrin localization and is required for v-ATPase to lower the pH in phagocytic endosomes, thus establishing an antiviral state. [57] High expression of both HP-PRRSV and AP-PRRSV suggests that IFITM3 may help inhibit PRRSV entry and promote resistance to HP-PRRSV infection. Peptidyl-prolyl cis-trans isomerase B (PPIB, vH/vA = 0.406) has peptidyl-prolyl cis-trans isomerase activity and is involved in protein folding and protein peptidyl-prolyl isomerization, which can accelerate protein folding. [58] PPIB is positively regulated by viral genome replication, is involved with the viral processes of hepatitis C virus (HCV), and interacts with and stimulates the RNA-binding activity of HCV NS5B. PPIB is critical for efficient replication of the HCV genome. [59] Compared with the control group, PPIB was downregulated in the HP-PRRSV and AP-PRRSV groups. The abundance of PPIB was lower in the HP-PRRSV group than the AP-PRRSV group, suggesting an association with PRRSV replication.

Potential Different Approaches Used by Both PRRSVs
String network analysis was used to elucidate interactions among differentially expressed proteins. The essential factors and receptors involved in the entry of PRRSV are reportedly CD163 (Q2VL90), nmHC II-A (MYH9, F1SKJ1), sialoadhesin (SN, A7LCJ3), CD151 (F1RYZ1), and vimentin (VIM) (P02543). [60][61][62][63][64][65] We assessed the involvement of these MPs together with proteins detected by string analysis on viral entry. Receptor proteins involved in PRRSV entry are in Figure 9  To better understand the characteristics of pathway-like regions, the GenomRNAi database was analyzed to annotate detected proteins. [66] Regions containing proteins with higher abundance in HP-PRRSV-infected PAMs are green and in AP-PRRSV-infected PAMs are red ( Figure 9A). A search of the GenomRNAi database determined that linked proteins IK (IK cytokine), WD repeat domain 12, dyskerin pseudouridine synthase 1, and adenylosuccinate lyase (ADSL) in the green region of Figure 9A shared similar annotations and functioned to decrease expression of NF-κB or IL-8, [67] which are related to inflammation. [44,45] Hence, further analysis of these linked proteins may help explain differential mechanisms between HP-PRRSV and AP-PRRSV infection. Other proteins in the green-colored region of 9A are reported to influence virus infection. Nucleolar and coiled-body phosphoprotein 1 increases Sindbis virus infection [68] and ADSL increases human papilloma virus 16-GFP infection. [69] The lower abundance of nucleolar and coiled-body phosphoprotein and ADSL in the HP-PRRSV-infected group may be related to the host immune response repression of PRRSV and may be used by HP-PRRSV during infection to self-upregulate. The red region of Figure 9A contains some linked proteins with some Figure 9. Protein-protein interaction network determined by String software showing interactions among differentially expressed proteins between two virulent groups with CD163, VIM, and MYH9 added. Red, protein abundance higher in HP-PRRSV than AP-PRRSV; green, protein abundance in HP-PRRSV lower than in AP-PRRSV. Line colors represent type of evidence for association: green, neighborhood; red, fusion; purple, experimental; light blue, database; black, expression; blue, co-occurrence; yellow, text mining. Gene abbreviations are shown. vH, protein abundance in HP-PRRSV; vA protein abundance in AP-PRRSV. vH/vA = 0, detected only in AP-PRRSV; vH/vA = null, detected only in HP-PRRSV.
information. Uridine monophosphate synthetase decreases HIV-1 infection, [70] while programmed cell death 5 decreases HCV replication, [71] both downregulate NF-κB expression. [67] The lower abundance of these proteins in the HP-PRRSV group suggested that HP-PRRSV escaped inhibition through an unknown mechanism. Myosin-9 (MYH9) appears to function in cytokinesis, cell shape, and specialized functions such as secretion and capping. [72] MYH9 is an important factor for PRRSV infection and interacts with GP5 of PRRSV, [65] although the underlying mechanisms remain unclear because of limited research. As shown by our results, MYH9 was related to VCL in AP-PRRSV and ACTL6 in HP-PRRSV. VCL is an actin filament (F-actin)binding protein involved in cell-matrix adhesion and cell-cell adhesion in humans, regulation of cell-surface E-cadherin expression, and potentiation of mechanosensing, and may be important in cell morphology and locomotion by promoting binding with actin, alpha-catenin, cadherin, dystroglycan, and ubiquitin protein ligase. [73] In HIV research, transient overexpression of VCL reduced the susceptibility of human cells to infection with HIV-1 and negatively affected paxillin phosphorylation and limited retroviral infection. [74] Just like HIV-1, in our study, it was demonstrated that over expression of VCL could inhibit the replication of both HP-PRRSV and the attenuated PRRSV in the mRNA level. There were obvious differences in the inhibiting ability for HP-PRRSV and its attenuated strain. ACTL6A, which is an actin-like protein 6A, is involved in transcriptional www.advancedsciencenews.com www.proteomics-journal.com activation and repression of select genes by chromatin remodeling (alteration of DNA-nucleosome topology) and mainly functions in chromatin binding and transcription coactivator activities. [75] We found that MYH9 may regulate AP-PRRSV and HP-PRRSV infection via different pathways. Through Super Pathways annotation (www.genecards.org), VCL and ACTL6A were identified as participants in the IL-3, IL-5, GM-CSF, and TNF-α/NF-κB signaling pathways, where they may help with differential mechanisms of PRRSV infection with different virulence.
Using the GenomRNAi database, some linked proteins were found to be related to viruses or inflammation ( Figure 9B). In the green region, GNAS complex locus (GNAS) increased human papilloma virus 16-GFP, [69] and serine/arginine repetitive matrix 2 decreased HCV infection, influenza A replication and viral numbers, and IL-8 expression, [76,77] Small nuclear ribonucleoprotein U1 subunit 70 decreased influenza A replication and viral numbers. [77] IRF-3 decreased infection by HCV, West Nile virus, and Dengue virus. [76,78] In the red region, HSF1 decreased HCV replication. [76] protein kinase CAMP-activated catalytic subunit alpha decreased VSV infection. [79] Signal transducer and activator of transcription 6 decreased IL-8 expression. RAP2A, MINK1, and GNAS are regulatory proteins in the RAS pathway. RAP2A and MINK1 were detected only in the HP-PRRSV group, whereas GNAS was detected only in the AP-PRRSV group, in accordance with a previous report that stimulation of RAS increases the replication ability of HCV by reducing IFN-JAK-STAT pathway activity. [80] We proposed that PRRSV was also related to the RAS pathway, and differed between HP-PRRSV and AP-PRRSV.
CREB-binding protein (CREBBP) is composed of doublestranded RNA-activated transcription factor with IRF-3, and double-stranded RNA-activated transcription factor is activated in many virus-infected cells to promote apoptosis. [81,82] CREBBP may interact with human herpes virus 8 vIRF-1, which could inhibit the binding of CREBBP to IRF-3. [83] Our results showed that IRF-3 abundance was tenfold lower in the HP-PRRSV than the AP-PRRSV group, which we proposed was a protective mechanism of PRRSV to escape from host immunity and ensure survival after viral infection. The ability of HP-PRRSV HuN4 to induce cell apoptosis is stronger than classical PRRSV (CH-1a) in immune organs and lungs of piglets. [84] We concluded that HP-PRRSV might interact with CREBBP and decrease IRF-3 expression to escape from host immunity and cause severe damage to the cell, highlighting a difference from AP-PRRSV.
Notable differences exist during the infection of AP-PRRSV and HP-PRRSV. HP-PRRSV could inhibit host immune function and evade the immune response via unknown mechanism. [85,86] Label-free MS was performed using AP-PRRSV-infected and HP-PRRSV-infected PAMs. This is the first attempt to explore the differential characteristics between HP-PRRSV and its attenuated PRRSV infected PAMs focusing on membrane proteins. By analyzing detected proteins in HP-infected, AP-PRRSV-infected, and control group, proteins related to the immune response or virus replication were identified that may elucidate unique pathways used by different virulent PRRSVs for cell entry, virus replication, and immune escape mechanisms. Researches on these detected proteins will help with the elucidation of the identity, the expression abundance, and significance of them, future study will be focused on functions of these key membrane proteins to deepen our understanding of differential mechanisms between HP-PRRSV and AP-PRRSV infection.