Skip to main content

Proteomics analysis reveals heat shock proteins involved in caprine parainfluenza virus type 3 infection



Caprine parainfluenza virus type 3 (CPIV3) is major pathogen of goat herds causing serious respiratory tract disease and economic losses to the goat industry in China. We analyzed the differential proteomics of CPIV3-infected Madin-Darby bovine kidney (MDBK) cells using quantitative iTRAQ coupled LC-MS/MS. In addition, four DEPs were validated by qRT-PCR and western blot analysis.


Quantitative proteomics analysis revealed 163 differentially expressed proteins (DEPs) between CPIV3-infected and mock-infected groups (p-value < 0.05 and fold change > 1.2), among which 91 were down-regulated and 72 were up-regulated. Gene ontology (GO) analysis showed that these DEPs were involved in molecular functions, cellular components and biological processes. Biological functions in which the DEPs were involved in included diseases, genetic information processing, metabolism, environmental information processing, cellular processes, and organismal systems. STRING analysis revealed that four heat shock proteins (HSPs) included HSPA5, HSPA1B, HSP90B1 and HSPA6 may be associated with proliferation of CPIV3 in MDBK cells. qRT-PCR and western blot analysis showed that the selected HSPs were identical to the quantitative proteomics data.


To our knowledge, this is the first report of the proteomic changes in MDBK cells after CPIV3 infection.


In August 2013, an outbreak of severe goat respiratory disease occurred throughout the major goat herd regions of eastern China. The causative agent was identified as a novel strain of parainfluenza virus type 3 (PIV3) and was designated as caprine parainfluenza virus type 3 (CPIV3) strain JS2013 [1]. The infected goats exhibited high fever, coughing, nasal discharge and dyspnea. Necropsy of the infected goats showed mild to moderate gross lesions in the lungs, and increased amounts of secretion in the tracheas and bronchia were also observed. Genome sequence alignment and phylogenetic analysis revealed that the genome of CPIV3 strain JS2013 showed only 73.3–75.5% identity with BPIV3 and HPIV3 strains [2]. Based on phylogenetic analysis, this pathogen was designated as CPIV3, a member of the PIV3 group belonging to the Respirovirus genus within the Paramyxiviridae family. Moreover, we further demonstrated that CPIV3 strain JS2013 can be transferred horizontally between adjacent pens [3]. Recently, a seroprevalence study using 2919 serum samples in China reported a CPIV3 prevalence of 39.9% in goats [4]. Another study reported that 35% of nasal swabs and serum samples from clinically diseased goats were positive for CPIV3 by quantitative RT-PCR (qRT-PCR) [5]. It is noteworthy that the spread of CPIV3 has caused heavy economic losses in China [6].

To understand the pathogenesis of viral infection, research on virus-host interaction is critical. Virus infection can dramatically affect host cell morphology, transcription and translation patterns, the cytoskeleton, the cell cycle and innate immune responses of the host, the apoptosis pathway, and may also cause inflammation and alter stress responses [7]. Many functional and morphological changes in host cells are associated with significant changes in the patterns of expression of host cells. Therefore, information on proteome changes in the host following CPIV3 infection may be crucial to understand the host response to viral pathogenesis. In recent years, comparative proteomic analysis has emerged as a valuable tool for the establishment of the global host protein profiles in response to virus infection [8]. This technique has been widely used to investigate proteome changes in cow, yak, buffalo, goat and camel milk [9], and peste des petits ruminants virus (PPRV)-infected Vero cells [10], based on the isobaric tags for relative and absolute quantification (iTRAQ) method. In addition, this technique has also been widely employed to examine the mechanisms of viral infection through comparative investigation of the proteome changes, for example, in the case of Crimean-Congo hemorrhagic fever virus (CCHFV) [11] and bovine respiratory syncytial viruses (BRSV) [12].

However, to the best of our knowledge, no previous study has analyzed the proteomic changes in CPIV3-infected MDBK cells. Proteomic techniques are effective tools to characterize protein expression profiles, and have been widely used to investigate disease-associated proteins [13, 14]. Among current proteomics methods, quantitative high-throughput proteomics approaches are useful for the analysis of infection-associated proteins [15, 16]. In our current study, we used a quantitative proteomics approach based on an iTRAQ tandem mass spectrometry (MS/MS) technique to identify differentially expressed proteins (DEPs) between CPIV3-infected and mock-infected MDBK cells. The functions of the DEPs were analyzed to determine whether they might be associated with CPIV3 infection [17]. Our findings provide valuable insight into the changes in cellular processes that occur during CPIV3 infection.


CPIV3 propagation in MDBK cells

The kinetics of CPIV3 propagation in MDBK cells were observed by monitoring the CPE at 24, 48 and 72 h post infection (hpi) (Fig. 1a), a minimal CPE was visible at 24 hpi, whereas an obvious CPE was observed at 48 hpi, and at 72hpi, almost all cells were disrupted. The TCID50 showed that the viral titer reached 103.5 TCID50/ml at 24 hpi, peaked at 107.0 TCID50/ml at 72 hpi and then declined (Fig. 1b). To ensure a higher proportion of infected cells and to avoid an excessive CPE, we selected 24 hpi as the time point under our infection conditions for further proteomic analysis.

Fig. 1

Confirmation of CPIV3 infection in MDBK cells. a A CPE was observed in MDBK cells at 24, 48 and 72 h after CPIV3 infection (MOI = 1), with mock-infected cells included as a control. b One-step growth curve of CPIV3 strain JS2013 in MDBK cells

Identification and annotation of proteins

We detected 8153 proteins and quantified 4109 proteins, including 28,815 peptides (Additional file 1: Figure S1). Detected proteins were annotated according to the GO database in the following categories: cellular components (CC), biological processes (BP), and molecular functions (MF) (Additional file 2: Figure S2). The top 20 pathways containing the largest number of proteins among the 8153 proteins were annotated according to KEGG (Additional file 3: Figure S3). Based on the KOG, 830 of the proteins were annotated as being involved in information storage and processing, 1545 were annotated as cellular processes and signaling, 581 were annotated as metabolism, and 699 were annotated as poorly characterized (Additional file 4: Figure S4 and Data Sheet 5). Furthermore, the cutoff criteria considered for the DEPs were set with an adjusted p-value of < 0.05 and a ratio of > 1.2-fold difference. Among the DEPs, 163 proteins from the two sets of biological replicates overlapped and were subsequently adjusted for multiple testing according to the stringent method of Benjamini and Hochberg [18]. Of these, 72 proteins were up-regulated and 91 proteins were down-regulated based on our criteria for the identification of DEPs in the MDBK-infected and mock-infected groups using the iTRAQ-MS/MS approach. Protein ratios were presented as CPIV3-infected/mock-infected. An average V/C ratio > 1 represented up-regulated proteins and an average V/C ratio < 1 represented down-regulated proteins. A list of DEPs information is shown in Table 1. DEPs between the two groups are shown as heat map and scatterplot (Additional file 5: Figure S6 and S7). Finally, the DEPs displaying the greatest increase and decrease in expression in the CPIV3-infected MDBK cells were FAM81B protein (1:0.118) and the DEP displaying the greatest decrease in expression in the CPIV3-infected MDBK cells was carboxypeptidase (1:1.206).

Table 1 Statistically significant DEPs identified by iTRAQ analysis of MDBK cells infected with CPIV3

GO analysis of the DEPs

The molecular functional classes and subcellular locations of the 163 DEPs were analyzed using UniProt and the GO database. The 163 DEPs were annotated into the categories: cellular component, biological process, or molecular function, and the distribution of up-regulated and down-regulated proteins among these GO annotations are shown in Additional file 6: Figure S8. GO enrichment annotation comparisons were performed to elucidate the characteristics of the altered proteins in MDBK cells induced by CPIV3 infection, to determine any associations with virulence and pathogenicity. In terms of biological process annotation, DEPs were mainly involved in cell aggregation, cellular processes, cellular component organization or biogenesis, locomotion, metabolic processes, multicellular organismal processes and reproductive processes; in terms cellular component annotation, DEPs were mainly involved in the cell part, extracellular region part, membrane part, organelle part, protein-containing complex and supramolecular complex; in terms of molecular function annotation, DEPs were mainly involved in binding, catalytic activity, molecular carrier activity and transporter activity (Fig. 2).

Fig. 2

GO enriched histogram of DEPs. Each column in the figure is a GO terms, the abscissa text indicates the name and classification of GO, and the height of the column indicates the enrichment rate. The color indicates the significance of the enrichment (p-value). The darker the color, the more significant the enrichment of the GO term (*P < 0.05; **P < 0.01; ***P < 0.001)

KEGG (Kyoto encyclopedia of genes and genomes) pathway analysis of the DEPs

The KEGG pathway is a collection of pathway maps that represent molecular interactions and reaction networks in cell line. The 93 DEPs identified were annotated, and mapped to a total of six KEGG pathway categories, which included metabolism, disease, genetic information processing, cellular processes, environmental information processing, and organismal systems pathway categories (Additional file 7 Data Sheet 9). The enrichment annotation protein pathway information is shown in Fig. 3. The results showed that most of the abundant KEGG terms were involved in biological processes such as the p53 signaling pathway, microRNAs in cancer, alanine, aspartate and glutamate metabolism, nitrogen metabolism, the estrogen signaling pathway, mineral absorption and thyroid hormone synthesis. Functional classification by KEGG showed that the upregulated and downregulated proteins could be divided among six distinct functional sets: environmental information processing, cellular processes, metabolism, genetic information processing, organismal systems and human diseases (Fig. 4).

Fig. 3

KEGG enrichment annotation of the DEPs. Each column in the figure is a pathway. The abscissa text indicates the name and classification of the pathway, and the height of the column indicates the enrichment rate. The color indicates the significance of the enrichment (p-value). The darker the color, the more significant the enrichment of the pathway (*P < 0.05; **P < 0.01; ***P < 0.001)

Fig. 4

Functional characterization of DEPs. a Cellular processes, metabolism and organismal systems. b Environmental information processing, genetic information processing and human diseases. More information is available in Additional file 5: Figure S7

STRING analysis of the relationships between DEPs

With the goal of exploring the potential protein network connections for the differentially regulated proteins in detail, the STRING tool was used. The differentially regulated proteins were mainly mapped to four functional networks (Fig. 5). A specific network had at least four “focus” proteins (HSPA5, HSPA1B, HSP90B1 and HSPA6). The networks of interest corresponded to: cell-to-cell signaling, hereditary disorder, cell death and survival, cardiovascular disease, cellular developmental, RNA post-transcriptional modification, cellular growth and proliferation.

Fig. 5

Specific network analysis of proteins significantly altered in CPIV3-infected cells. The network of DEPs with STRING analysis. Each node represents a protein in the graph, each line represents the interaction between proteins, and the wider the line, the closer the relationship

Confirmation of proteomic data by qRT-PCR

Alterations in the expression of a protein may be owing to a change in its mRNA levels. To confirm the results of the proteomic analysis by mRNA expression, transcriptional alterations in four selected proteins were measured by qRT-PCR. The qRT-PCR analysis showed that no difference in the ratio of these mRNAs between the CPIV3 infected group and the mock infected group were consistent with those obtained using quantitative proteomics analysis (Fig. 6). The mRNA expression of HSPA5, HSP90B1, HSPA1B and HSPA6 were increased in CPIV3-infected MDBK cells. Therefore, the trends in the mRNA expression were consistent with those in their corresponding proteins.

Fig. 6

qRT-PCR analysis of mRNA expression in the CPIV3-infected and mock-infected groups. The cells were collected at 24 hpi for qRT-PCR to analyze the relative mRNA expression of the HSPA5, HSP90B1, HSPA1B and HSPA6 genes. The GAPDH gene was included as a control housekeeping gene for the normalization of samples . Error bars represent standard deviations

Western blot analysis of HSPA1B

We analyzed the expression levels of HSPA1B (up-regulated) in CPIV3-infected MDBK cells (Fig. 7) by western blot at 24 h and 48 h. Figure 7 shows that HSPA1B was up-regulated in CPIV3-infected MDBK cells at 24 h and 48 h. The results were consistent with those obtained using the iTRAQ labeled LC-MS/MS system.

Fig. 7

Analysis of HSPA1B expression levels in CPIV3-infected and control cells by western blot analysis at 24 h and 48 h. Protein samples were separated by SDS-PAGE. Western blot analysis was performed using antibodies to the HSPA1B protein.β-actin protein was detected as an internal control


Proteomic techniques have become significant methodologies for determining cellular protein interactions and host cellular pathophysiological processes following virus infection [19, 20]. As a general rule, no important host cell membrane rearrangement or cytoskeleton collapse is observed following virus infection but the point at which a high virus yield is obtained is considered as the best time for proteomic analysis [21, 22]. Taking this substantial evidence into consideration, cell samples at 24 hpi were chosen for further proteomic analysis. Based on our study, the expression levels of 163 DEPs were found to be significantly altered in CPIV3-infected cells. The results of GO, KEGG pathway and STRING analysis predicted that these DEPs pertaining to different types of functional categories and signal pathways. Western blot and qRT-PCR were also applied to validate some differential proteins at the mRNA and protein levels. To date, no analysis has been reported of the differential proteomes of MDBK cells infected with CPIV3. Our data may provide an overview of the proteins altered in expression during the host response to CPIV3 infection and may provide insight in the process of CPIV3 pathogenesis.

Studies have shown that HSPs may play an important role in virus host cell interactions during in vivo and in vitro infection [23, 24]. Inhibitors of HSP90 can inhibit herpes simplex virus type 1 (HSV-1) infection [25]. Bovine viral diarrhea virus (BVDV) structural proteins comprise the C nucleocapsid protein and three envelope glycoproteins, Erns, E1 and E2 [26]. A previous study found that HSP110 enhanced the presentation of E2 to CD4 T cells in vitro to improve the immunogenicity of an E2 vaccine in cattle [27]. Previous work demonstrated that HSP70 is actively released into the extracellular milieu and acts as a cytokine and peptide adjuvant, thereby promoting both the innate and adaptive immune responses [28]. In our analysis, four proteins (HSPA5, HSPA1B, HSP90B1 and HSPA6) were identified following CPIV3 infection. HSP90B1 is proposed to be associated with poor survival from hepatocellular carcinoma (HCC), whereas high levels of HSPA5 and HSPA6 may be associated with earlier recurrence of HCC [29]. HSPA1B, also known as heat shock protein 72, is a member of the HSP70 family. HSP70 expression levels rapidly increased in response to cellular stresses such as heat shock, or in response to certain viral infections [30,31,32,33].

In the current study, HSP70 was rarely detected in the mock-infected group, whereas it was notably present in the CPIV3 group. CPIV3 infection resulted in the up-regulated secretion of exosomes and packaging of the viral proteins into exosomes, and these results suggested that CPIV3 infection may enhance HSP70-mediated exosome release (unpublished data). In addition, HSP70 is actively released into the extracellular milieu, thereby promoting innate and adaptive immune responses [34]. In this study, HSPA5, HSPA1B, HSP90B1 and HSPA6 were up-regulated at 24 hpi to various degrees following CPIV3 -infection of MDBK cells. Different expression levels of HSPA1B were detected by western blot analysis at 24 hpi and 48 hpi after CPIV3 -infection of MDBK cells. This may indicate that HSPA1B affects the proliferation of CPIV3 in MDBK cells. HSPA1B is an endogenous ligand for toll-like receptor TLR4, thereby stimulating innate immunity [35], and HSPA1B regulates the NF-κB pathway via TLR2 and TLR4 in fibroblasts. However, fibroblasts and macrophages interact with each other to mediate the immune response. Activation of the NF-κB pathway then results to in enhanced secretion of pro-inflammatory cytokines (TNF-α, IL-6 and IL-1β) and neutrophil chemoattractant MIP-2 and Cxcl1 from macrophages [36]. This evidence indicates that HSPA1B may be associated with the proliferation of CPIV3 in MDBK cells through an ability to interact with key components of the NF-κB pathway, moreover, those involved in innate immunity, but the detailed mechanism remains unknown. However, the detailed functions of these pathways and proteins changes in CPIV3 infection therefore requires further verification.


The proteomic changes in CPIV3-infected MDBK cells were analyzed using iTRAQ combined with LC-MS/MS. To the best of our knowledge, this is the first time proteomics has been used to explore the virus–host protein interaction network in CPIV3-infected MDBK cells. The results revealed 163 DEPs, among which 72 were up-regulated and 91 were down-regulated. In addition, four DEPs were validated by qRT-PCR and HSPA1B was validated by western blot analysis. These results were consistent with those of label-free LC-MS analysis. Our analyses of the DEPs were descriptive, and further functional investigations are required to elucidate the pathogenic mechanisms and cellular responses to CPIV3 infection.


Cell culture and virus infection

CPIV3 strain JS2013 isolated in Jiangsu Province was used for virus infection. MDBK cells were cultured in Dulbecco’s modified Eagle’s medium (DMEM; Sigma, CA, USA) supplemented with 10% fetal bovine serum (FBS; HyClone, UT, USA), at 37 °C in an atmosphere of 5% CO2 [2]. When the cells grow to 70–80% confluence, they were inoculated with CPIV3 at a multiplicity of infection (MOI) of 1. After 1 h of adsorption, infected cells were maintained in fresh medium containing 2% FBS. Uninfected cells were used as a control. The CPIV3- or mock-infected cells were collected at 24 hpi. Viral propagation was confirmed by the observation of a cytopathic effect (CPE).

Protein sample preparation and labeling with iTRAQ reagent

The CPIV3- and mock-infected cell samples were washed three times with cold phosphate-buffered saline (PBS) and then treated with lysis buffer containing 8 M urea, 4% CHAPS, 2 M thiourea, and 30 mM Tris-HCl on ice for 30 min until the cell line were completely lysed. The supernatant was collected by centrifugation at 12000×g for 30 min at 4 °C after ultrasonication treatment for 2 min. The protein concentration in the supernatants was quantified using the Bradford protein assay. After reduction and cysteine-blocking as described in the iTRAQ protocol (AB Sciex, Concord, ON, USA), solutions containing 100 μg protein were digested overnight at 37 °C with sequence grade modified trypsin (Promega, Madison, WI, USA) and then labeled with different iTRAQ tags. The labeled samples were then mixed and dried with a rotary vacuum concentrator.

LC-MS analysis

Ten microliters (μl) of each fraction were analyzed by Q Exactive (Thermo, USA) mass spectrometer coupled to a Proxeon Biosystem Easy-nLC 1200 (Thermo Fisher Scientific, Waltham, MA, USA) in the LC-MS experiments. The peptide mixture (5 g) was loaded onto a C18 column (75 μm × 25 cm, Thermo,USA) packed with RP-C18 (5 m) resin in buffer A (2% ACN with 0.1% formic acid), and eluted with a linear gradient of buffer B (80% ACN with 0.1% formic acid) at a flow rate of 300 nl/min for 120 min using IntelliFlow technology. The equate underwent electrospray ionization for LC-MS analysis. The MS/MS instrument was run in the peptide recognition mode, and the spectra were acquired using a data-dependent top-20 method based on the selection of the most abundant precursor ions from the survey scan (350–1300 m/z) for HCD fragmentation. Determination of the target value was based on the predictive automatic gain control, and the dynamic exclusion duration was 18 s. Survey scans were acquired at a resolution of 70,000 at m/z 200, and the resolution for the HCD spectra was set to 17,500 at m/z 200. The 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%. Thermo Xcalibur 4.0 (Thermo, USA) was used to collect MS analysis data via DDA mode.

Data analysis

The MS data were analyzed using Proteome Discoverer™ software 2.1. When the library was searched, the raw file was submitted to the Proteome Discoverer server searched against the Uniprot Bos taurus database (197,939 total sequences, downloaded April 26, 2018). The following parameters were used for protein identification: a precursor mass tolerance of 20 ppm; a fragment mass tolerance of 0.05 Da; trypsin digestion; max. Missed cleavage sites of 2; the variable dynamic modifications included oxidation (M), iTRAQ8plex (Y) and acetyl (protein N-terminus), and the fixed static modifications included carbamidomethyl (C), iTRAQ8plex (K) and iTRAQ8plex (N-term). The cutoff for the global false discovery rate (FDR) for peptide and protein identification was set to 0.01. The value of the quantitative ratio for each protein relative to the internal reference was calculated, and averaged to obtain the quantitative ratio (V/C) of the proteins identified in the treatment groups [37]. Proteins with a fold change > 1.2 and a p-value < 0.05 were considered to shows significantly different expression. Auto bias-correction was executed to decrease the artificial error. Statistical analysis was performed using Excel 2007 software. The DEPs were annotated using gene ontology (GO) and KEGG database. The Cluster of Orthologous Groups of proteins (COG or KOG) were retrieved, and mapped to pathways in the KEGG database [38]. In addition, DEPs were analyzed using STRING for predicting functional association networks of proteins.

CPIV3 yield quantification

MDBK cells were seeded in 96-well plates and incubated for 24 h. Then, CPIV3 samples were 10-fold serially diluted and added to each well in quadruplicate. MDBK cells exhibit CPE were scored positive for viral growth and the TCID50 was calculated by the Reed–Muench method [39].

mRNA quantitation by qRT-PCR

Total cellular RNA was extracted from the CPIV3-infected and mock-infected MDBK cells using Transzol UP reagent (Transgen Co. Ltd., Beijing, China) according to the manufacturer’s protocol. Specific primers for amplifying various genes were as follows: for GAPDH mRNA analysis, 5′-GATTGTCAGCAATGCCTCCT-3′ (forward) and 5′-GGTCATA AGTCCCTCCACGA-3′ (reverse) were used; for HSPA5 mRNA analysis, 5′-GTGCCCACCA AGAAGTCTCA-3′ (forward) and 5′-CTTTCGTCAGGGGTCGTTC A-3′ (reverse) were used; for HSP90B1 mRNA analysis, 5′-TCAAGGGTGTTGTGGACTCG-3′ (forward) and 5′-GCT GAAGTGTCTCACGGG AA-3′ (reverse) were used; for HSPA1B mRNA analysis, 5′-AGTC GGACATGAAGCACTGG-3′ (forward) and 5′-TCACCTGCACCTTAGGCTTG-3′ (reverse) were used; and for HSPA6 mRNA analysis, 5′-AGGACAGGCGCAAAGTACAA-3′ (forward) and 5′-TGCTCCAGCTCCCTCTTTTG-3′ (reverse) were used. GAPDH was employed as an internal reference gene. The first-strand cDNA was synthesized via PrimeScript™ RT Master Mix (TaKaRa, Dalian, China). Then qRT-PCR was performed using the SYBR Premix Ex Taq™ II Kit (TaKaRa) on an ABI Step One thermocycler (Applied Biosystems, CA, USA). The relative expression level of each mRNA was calculated by the 2-ΔΔct method. Three independent biological replicates were performed for each gene.

Western blot analysis

To further verify the variation in the DEPs identified by the proteomic approaches, HSPA1B was selected for western blot analysis. The CPIV3- and mock-infected cells were collected at 24 and 48 hpi. Equivalent amounts of cell lysate from each sample were collected. After measuring the protein concentrations, equivalent amounts of cellular proteins were separated by SDS-PAGE and transferred onto nitrocellulose PVDF membranes (Millipore, USA). The membranes were incubated overnight at 4 °C with primary rabbit polyclonal antibodies of anti-HSPA1B (Biyotime, Shanghai, China). Then the membranes were further incubated for 1 h with horseradish peroxidase-conjugated goat anti-rabbit secondary antibody (BIOSS, Beijing, China). The protein bands were detected using the ECL Detection Kit (Vazyme, Nanjing, China). β-actin protein was used as an internal control.



Cluster of orthologous groups of proteins


Cytopathic effect


Caprine parainfluenza virus type 3


Differentially expressed proteins


Dulbecco’s modified Eagle’s medium


False discovery rate


Gene ontology


Heat shock proteins


Isobaric tags for relative and absolute quantification


Kyoto encyclopedia of genes and genomes


Madin-Darby bovine kidney


Multiplicity of infection


  1. 1.

    Li W, Mao L, Cheng S, Wang Q, Huang J, Deng J, Wang Z, Zhang W, Yang L, Hao F, et al. A novel parainfluenza virus type 3 (PIV3) identified from goat herds with respiratory diseases in eastern China. Vet Microbiol. 2014;174(1–2):100–6.

    Article  Google Scholar 

  2. 2.

    Yang L, Li W, Mao L, Hao F, Wang Z, Zhang W, Deng J, Jiang J. Analysis on the complete genome of a novel caprine parainfluenza virus 3. Infect Genet Evol. 2016;38:29–34.

    Article  Google Scholar 

  3. 3.

    Li W, Hao F, Mao L, Wang Z, Zhou T, Deng J, Li J, Zhang W, Yang L, Lv Y, et al. Pathogenicity and horizontal transmission studies of caprine parainfluenza virus type 3 JS2013 strain in goats. Virus Res. 2016;223:80–7.

    CAS  Article  Google Scholar 

  4. 4.

    Mao L, Li W, Zhou T, Yang L, Hao F, Li J, Zhang W, Luo X, Jiang J. Development of a blocking ELISA for caprine parainfluenza virus type 3. J Virol Methods. 2017;250:59–65.

    CAS  Article  Google Scholar 

  5. 5.

    Li J, Li W, Mao L, Hao F, Yang L, Zhang W, Jiang J. Rapid detection of novel caprine parainfluenza virus type 3 (CPIV3) using a TaqMan-based RT-qPCR. J Virol Methods. 2016;236:126–31.

    CAS  Article  Google Scholar 

  6. 6.

    Mao L, Yang L, Li W, Liang P, Zhang S, Li J, Sun M, Zhang W, Wang L, Zhong C, et al. Epidemiological investigation and phylogenetic analysis of caprine parainfluenza virus type 3 in sheep of China. Transbound Emerg Dis. 2019.

  7. 7.

    Zheng J, Sugrue RJ, Tang K. Mass spectrometry based proteomic studies on viruses and hosts--a review. Anal Chim Acta. 2011;702(2):149–59.

    CAS  Article  Google Scholar 

  8. 8.

    Han K, Zhao D, Liu Y, Liu Q, Huang X, Yang J, An F, Li Y. Quantitative proteomic analysis of duck ovarian follicles infected with duck Tembusu virus by label-free LC-MS. Front Microbiol. 2016;7:463.

    PubMed  PubMed Central  Google Scholar 

  9. 9.

    Yang Y, Bu D, Zhao X, Sun P, Wang J, Zhou L. Proteomic analysis of cow, yak, buffalo, goat and camel milk whey proteins: quantitative differential expression patterns. J Proteome Res. 2013;12(4):1660–7.

    CAS  Article  Google Scholar 

  10. 10.

    Pandey A, Sahu AR, Wani SA, Saxena S, Kanchan S, Sah V, Rajak KK, Khanduri A, Sahoo AP, Tiwari AK, et al. Modulation of host miRNAs transcriptome in lung and spleen of Peste des Petits ruminants virus infected sheep and goats. Front Microbiol. 2017;8:1146.

    Article  Google Scholar 

  11. 11.

    Fernandez de Mera IG, Chaligiannis I, Hernandez-Jarguin A, Villar M, Mateos-Hernandez L, Papa A, Sotiraki S, Ruiz-Fons F, Cabezas-Cruz A, Gortazar C, et al. Combination of RT-PCR and proteomics for the identification of Crimean-Congo hemorrhagic fever virus in ticks. Heliyon. 2017;3(7):e00353.

    Article  Google Scholar 

  12. 12.

    Hagglund S, Blodorn K, Naslund K, Vargmar K, Lind SB, Mi J, Arainga M, Riffault S, Taylor G, Pringle J, et al. Proteome analysis of bronchoalveolar lavage from calves infected with bovine respiratory syncytial virus-insights in pathogenesis and perspectives for new treatments. PLoS One. 2017;12(10):e0186594.

    Article  Google Scholar 

  13. 13.

    Sun D, Zhang H, Guo D, Sun A, Wang H. Shotgun proteomic analysis of plasma from dairy cattle suffering from footrot: characterization of potential disease-associated factors. PLoS One. 2013;8(2):e55973.

    CAS  Article  Google Scholar 

  14. 14.

    He Y, Li W, Liao G, Xie J. Mycobacterium tuberculosis-specific phagosome proteome and underlying signaling pathways. J Proteome Res. 2012;11(5):2635–43.

    CAS  Article  Google Scholar 

  15. 15.

    Zeng S, Zhang H, Ding Z, Luo R, An K, Liu L, Bi J, Chen H, Xiao S, Fang L. Proteome analysis of porcine epidemic diarrhea virus (PEDV)-infected Vero cells. Proteomics. 2015;15(11):1819–28.

    CAS  Article  Google Scholar 

  16. 16.

    Linde ME, Colquhoun DR, Ubaida Mohien C, Kole T, Aquino V, Cotter R, Edwards N, Hildreth JE, Graham DR. The conserved set of host proteins incorporated into HIV-1 virions suggests a common egress pathway in multiple cell types. J Proteome Res. 2013;12(5):2045–54.

    CAS  Article  Google Scholar 

  17. 17.

    Sun D, Shi H, Guo D, Chen J, Shi D, Zhu Q, Zhang X, Feng L. Analysis of protein expression changes of the Vero E6 cells infected with classic PEDV strain CV777 by using quantitative proteomic technique. J Virol Methods. 2015;218:27–39.

    CAS  Article  Google Scholar 

  18. 18.

    Cho YE, Singh TS, Lee HC, Moon PG, Lee JE, Lee MH, Choi EC, Chen YJ, Kim SH, Baek MC. In-depth identification of pathways related to cisplatin-induced hepatotoxicity through an integrative method based on an informatics-assisted label-free protein quantitation and microarray gene expression approach. Mol Cell Proteomics. 2012;11(1):M111 010884.

    Article  Google Scholar 

  19. 19.

    Zhang X, Zhou J, Wu Y, Zheng X, Ma G, Wang Z, Jin Y, He J, Yan Y. Differential proteome analysis of host cells infected with porcine circovirus type 2. J Proteome Res. 2009;8(11):5111–9.

    CAS  Article  Google Scholar 

  20. 20.

    Maxwell KL, Frappier L. Viral proteomics. Microbiol Mol Biol Rev. 2007;71(2):398–411.

    CAS  Article  Google Scholar 

  21. 21.

    An K, Fang L, Luo R, Wang D, Xie L, Yang J, Chen H, Xiao S. Quantitative proteomic analysis reveals that transmissible gastroenteritis virus activates the JAK-STAT1 signaling pathway. J Proteome Res. 2014;13(12):5376–90.

    CAS  Article  Google Scholar 

  22. 22.

    Zhang LK, Chai F, Li HY, Xiao G, Guo L. Identification of host proteins involved in Japanese encephalitis virus infection by quantitative proteomics analysis. J Proteome Res. 2013;12(6):2666–78.

    CAS  Article  Google Scholar 

  23. 23.

    Braga ACS, Carneiro BM, Batista MN, Akinaga MM, Bittar C, Rahal P. Heat shock proteins HSPB8 and DNAJC5B have HCV antiviral activity. PLoS One. 2017;12(11):e0188467.

    Article  Google Scholar 

  24. 24.

    Rathore AP, Haystead T, Das PK, Merits A, Ng ML, Vasudevan SG. Chikungunya virus nsP3 & nsP4 interacts with HSP-90 to promote virus replication: HSP-90 inhibitors reduce CHIKV infection and inflammation in vivo. Antivir Res. 2014;103:7–16.

    CAS  Article  Google Scholar 

  25. 25.

    Zhong M, Zheng K, Chen M, Xiang Y, Jin F, Ma K, Qiu X, Wang Q, Peng T, Kitazato K, et al. Heat-shock protein 90 promotes nuclear transport of herpes simplex virus 1 capsid protein by interacting with acetylated tubulin. PLoS One. 2014;9(6):e99425.

    Article  Google Scholar 

  26. 26.

    Patton JT, Chizhikov V, Taraporewala Z, Chen D. Virus replication. Methods Mol Med. 2000;34:33–66.

    CAS  PubMed  Google Scholar 

  27. 27.

    McLaughlin K, Carr VB, Iqbal M, Seago J, Lefevre EA, Robinson L, Prentice H, Charleston B. Hsp110-mediated enhancement of CD4+ T cell responses to the envelope glycoprotein of members of the family Flaviviridae in vitro does not occur in vivo. Clin Vaccine Immunol. 2011;18(2):311–7.

    CAS  Article  Google Scholar 

  28. 28.

    Hunter-Lavin C, Davies EL, Bacelar MM, Marshall MJ, Andrew SM, Williams JH. Hsp70 release from peripheral blood mononuclear cells. Biochem Biophys Res Commun. 2004;324(2):511–7.

    CAS  Article  Google Scholar 

  29. 29.

    Yang Z, Zhuang L, Szatmary P, Wen L, Sun H, Lu Y, Xu Q, Chen X. Upregulation of heat shock proteins (HSPA12A, HSP90B1, HSPA4, HSPA5 and HSPA6) in tumour tissues is associated with poor outcomes from HBV-related early-stage hepatocellular carcinoma. Int J Med Sci. 2015;12(3):256–63.

    CAS  Article  Google Scholar 

  30. 30.

    Cheung RK, Dosch HM. The growth transformation of human B cells involves superinduction of hsp70 and hsp90. Virology. 1993;193(2):700–8.

    CAS  Article  Google Scholar 

  31. 31.

    Lefeuvre A, Contamin H, Decelle T, Fournier C, Lang J, Deubel V, Marianneau P. Host-cell interaction of attenuated and wild-type strains of yellow fever virus can be differentiated at early stages of hepatocyte infection. Microbes Infect. 2006;8(6):1530–8.

    CAS  Article  Google Scholar 

  32. 32.

    Liao WJ, Fan PS, Fu M, Fan XL, Liu YF. Increased expression of 70 kD heat shock protein in cultured primary human keratinocytes induced by human papillomavirus 16 E6/E7 gene. Chin Med J. 2005;118(24):2058–62.

    CAS  PubMed  Google Scholar 

  33. 33.

    Mayer MP. Recruitment of Hsp70 chaperones: a crucial part of viral survival strategies. Rev Physiol Biochem Pharmacol. 2005;153:1–46.

    CAS  Article  Google Scholar 

  34. 34.

    Mansilla MJ, Costa C, Eixarch H, Tepavcevic V, Castillo M, Martin R, Lubetzki C, Aigrot MS, Montalban X, Espejo C. Hsp70 regulates immune response in experimental autoimmune encephalomyelitis. PLoS One. 2014;9(8):e105737.

    Article  Google Scholar 

  35. 35.

    Rusai K, Banki NF, Prokai A, Podracka L, Szebeni B, Tulassay T, Reusz GS, Sallay P, Kormendy R, Szabo AJ, et al. Heat shock protein polymorphism predisposes to urinary tract malformations and renal transplantation in children. Transplant Proc. 2010;42(6):2309–11.

    CAS  Article  Google Scholar 

  36. 36.

    Zulaziz N, Azhim A, Himeno N, Tanaka M, Satoh Y, Kinoshita M, Miyazaki H, Saitoh D, Shinomiya N, Morimoto Y. Photodynamic therapy mediates innate immune responses via fibroblast-macrophage interactions. Hum Cell. 2015;28(4):159–66.

    CAS  Article  Google Scholar 

  37. 37.

    Unwin RD, Griffiths JR, Whetton AD. Simultaneous analysis of relative protein expression levels across multiple samples using iTRAQ isobaric tags with 2D nano LC-MS/MS. Nat Protoc. 2010;5(9):1574–82.

    CAS  Article  Google Scholar 

  38. 38.

    Kanehisa M, Goto S, Sato Y, Furumichi M, Tanabe M. KEGG for integration and interpretation of large-scale molecular data sets. Nucleic Acids Res. 2012;40(Database issue):D109–14.

    CAS  Article  Google Scholar 

  39. 39.

    Reed LJ, Muench H. A simple method of estimation of fifty percent endpoint. Am J Epidemiol. 1938;27:493–7.

    Article  Google Scholar 

Download references


We thank Kate Fox, DPhil, from Liwen Bianji, Edanz Group China (, for editing the English text of a draft of this manuscript.


This work was supported by the National Natural Science Foundation of China (31702272, 31802196), Natural Science Foundation of Jiangsu Province, China (BK20170595), Natural Science Foundation of Shandong Province, China (ZR2016CP08), and the National Key R&D Program of China (2016YFD0500908, 2018YFD0502100). The funders had no role in study design, in the collection, analysis and interpretation of data, in the writing the manuscript, or in the decision to submit the article for publication.

Availability of data and materials

The datasets contained in this study are available from the corresponding author upon request.

Author information




CZ and JL took part in all the experiments and wrote the manuscript. LM, ML, XZ and WL helped to designed the whole project and draft the manuscript. MS, FX, LY, WZ and ZL conducted cell culture and sample processing for sequencing. XJ conducted data analysis. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Jizong Li.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

Figure S1. Information on the detected proteins in CPIV3-infected MDBK cells (JPG 357 kb)

Additional file 2:

Figure S2. Detected proteins were annotated in the GO database (JPG 1406 kb)

Additional file 3:

Figure S3. The top 20 pathways annotated by KEGG (JPG 738 kb)

Additional file 4:

Figure S4 and Data Sheet 5. Proteins were annotated based on the KOG (ZIP 669 kb)

Additional file 5:

Figure S6 and S7 Heat map and scatterplot (ZIP 1116 kb)

Additional file 6:

Figure S8. GO annotations for the up-regulated and down-regulated proteins (JPG 2482 kb)

Additional file 7:

Data Sheet 9. The 93 DEPs were annotated into six KEGG pathway categories (XLS 32 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Zhong, C., Li, J., Mao, L. et al. Proteomics analysis reveals heat shock proteins involved in caprine parainfluenza virus type 3 infection. BMC Vet Res 15, 151 (2019).

Download citation


  • Caprine parainfluenza virus type 3
  • Madin-Darby bovine kidney cells
  • Proteomic analysis
  • iTRAQ
  • LC-MS/MS