Next Article in Journal
Development of a Highly Efficient Hybrid Peptide That Increases Immunomodulatory Activity Via the TLR4-Mediated Nuclear Factor-κB Signaling Pathway
Next Article in Special Issue
Flavescence Dorée Phytoplasma Has Multiple ftsH Genes that Are Differentially Expressed in Plants and Insects
Previous Article in Journal
Mechanisms of Metastatic Tumor Dormancy and Implications for Cancer Therapy
Previous Article in Special Issue
Dual RNA Sequencing of Vitis vinifera during Lasiodiplodia theobromae Infection Unveils Host–Pathogen Interactions
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 20
Issue 24
10.3390/ijms20246160
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Quantitative Proteomic Analysis Uncovers the Mediation of Endoplasmic Reticulum Stress-Induced Autophagy in DHAV-1-Infected DEF Cells

by
Jingjing Lan
1,2,
Ruihua Zhang
1,2,
Honglei Yu
1,2,
Jingyu Wang
1,2,
Wenxiang Xue
1,2,
Junhao Chen
1,3,
Shaoli Lin
4,
Yu Wang
5,
Zhijing Xie
1,2 and
Shijin Jiang
1,2,*
1
College of Veterinary Medicine, Shandong Agricultural University, Taian 271000, China
2
Shandong Provincial Key Laboratory of Animal Biotechnology and Disease Control and Prevention, Taian 271000, China
3
College of Public Health and Management, Weifang Medical University, Weifang 261042, China
4
Molecular Virology Laboratory, VA-MD College of Veterinary Medicine and Maryland Pathogen Research Institute, University of Maryland, College Park, MD 20742, USA
5
Department of Basic Medical Sciences, Taishan Medical College, Taian 271000, China
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2019, 20(24), 6160; https://doi.org/10.3390/ijms20246160
Submission received: 16 October 2019 / Revised: 3 December 2019 / Accepted: 4 December 2019 / Published: 6 December 2019
(This article belongs to the Special Issue Host–Pathogen Interaction)
Browse Figures
Versions Notes

Abstract

:
Autophagy is a tightly regulated catabolic process and is activated in cells in response to stress signals. Despite extensive study, the interplay between duck hepatitis A virus type 1 (DHAV-1) and the autophagy of host cells is not clear. In this study, we applied proteomics analysis to investigate the interaction mechanism between DHAV-1 and duck embryo fibroblast (DEF) cells. In total, 507 differentially expressed proteins (DEPs) were identified, with 171 upregulated proteins and 336 downregulated proteins. The protein expression level of heat shock proteins (Hsps) and their response to stimulus proteins and zinc finger proteins (ZFPs) were significantly increased while the same aspects of ribosome proteins declined. Bioinformatics analysis indicated that DEPs were mainly involved in the “response to stimulus”, the “defense response to virus”, and the “phagosome pathway”. Furthermore, Western blot results showed that the conversion of microtubule-associated protein 1 light chain 3-I (LC3-I) to the lipidation form of LC3-II increased, and the conversion rate decreased when DEF cells were processed with 4-phenylbutyrate (4-PBA). These findings indicated that DHAV-1 infection could cause endoplasmic reticulum (ER) stress-induced autophagy in DEF cells, and that ER stress was an important regulatory factor in the activation of autophagy. Our data provide a new clue regarding the host cell response to DHAV-1 and identify proteins involved in the DHAV-1 infection process or the ER stress-induced autophagy process.

1. Introduction

Duck hepatitis A virus (DHAV) is a member of the genus Avihepatovirus in the Picornaviridae family. DHAV is divided into three serotypes, namely, the world-wide traditional serotype called DHAV-1 [1,2], a serotype isolated in Taiwan called DHAV-2 [3], and a serotype isolated in South Korea and China called DHAV-3 [4,5]; no antigenic relationships have been found among them [3]. As a fatal rapidly spreading disease, DHAV-1 infection is characterized by liver petechiae and hepatitis in young ducklings and egg drop in laying duck flocks [6,7,8,9]. In order to control DHAV-1 infection, researchers have made great efforts to study the interactive mechanism between DHAV-1 and host cells [10,11,12].
Autophagy is a traditional mechanism that degrades protein and waste in cells. It has been reported previously that the endoplasmic reticulum (ER), as an integral and elaborate organelle for folding and modifying secretory proteins, can induce autophagy if it is damaged [13]. Viral infection can lead to disorder of the intracellular environment, such as the accumulation of misfolded protein or unfolded protein in ER or Ca2+ balance, so cells launch the ER stress and autophagy reaction against the infection. ER autophagy is a selective autophagy process with a key role in regulating the unfolded protein response (UPR), which is responsible for maintaining cell homeostasis [14]. Researchers have affirmed that ER stress and autophagy participate in numerous cell processes during virus infection, such as cell death, the immune response, and viral replication [15,16,17]. Moreover, Toll-like receptor and type I interferon production are triggered by autophagosome fusion with the lysosomal pathway [18]. Therefore, it is not surprising that viruses have evolved some evasion mechanisms to achieve infection. Recent studies have shown that some viruses can inhibit or evade autophagy, whereas some viruses can not only induce autophagy but even take advantage of it to promote virus replication [19,20,21]. Although DHAV-1 has been reported to induce apoptosis and the immune response [22,23], it is still necessary to find more sufficient evidence to elucidate the phenomenon in ER stress-induced autophagy, especially because there are currently no reports on this process.
Proteomic methods are a highly specific, effective, and universal technique, which do not require multi-step sample preparation [24]. Compared to RNA-seq, proteomic techniques can accurately reflect the abundances of downstream proteins, thus strategies focusing on protein quantification or/and post-translational modification have been widely applied in this area [25,26]. To date, several studies have focused on the interaction of viruses and host cells based on the proteomic method [27,28,29,30]. In this study, we focused on proteome changes of host proteins that were possibly involved in ER stress-induced autophagy in duck embryo fibroblast (DEF) cells, which are a natural primary target for DHAV-1. The quantification results, followed by gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, showed that differentially expressed proteins (DEPs) were mainly involved in cellular processes, cell stimulation, the immune response, lysosomes, phagosomes, and others. Transmission electron microscopy (TEM) analysis confirmed that double-membraned autophagy-like vesicles were generated in DHAV-1-infected DEF cells. Proteomics results and Western blot results indicated that DHAV-1 was involved in the process of ER stress-induced autophagy. Overall, these findings enhance the understanding of the pathogenic mechanism at the cellular level during DHAV-1 infection.

2. Results

2.1. DHAV-1 Infection

After DHAV-1 infection, the cell morphologies were observed. The results showed that DEF cells shrunk at 48 h post-infection (hpi) and shedded at 60 hpi (Figure 1A). An indirect immunofluorescence assay (IFA) showed that DHAV-1 successfully infected DEF cells, observed by green fluorescence, while the non-infected DEF cells showed no fluorescence (Figure 1B). The qRT-PCR results showed that the optimal amount of DHAV-1 infection was a multiplicity of infection (MOI) of 2 (Figure 1C). The messenger RNA (mRNA) copies of DHAV-1 reached a maximum at 48 hpi, decreased at 60 hpi, and reached the lowest level at 72 hpi (Figure 1D). Therefore, DEF cells at 48 hpi were chosen to further investigate the changes in the host proteins during DHAV-1 infection.

2.2. Tandem Mass Tag (TMT)-Based Proteomic Analysis

Based on TMT labeling-based proteomics technology, we aimed to investigate the global pattern of DEF cells either infected with DHAV-1 or not. Our experimental workflow is shown in Figure 2. In the identification, the pairwise Pearson’s correlation coefficients were greater than 0.7, suggesting sufficient reproducibility of the experiment (Figure 3A). From tandem mass spectrometry (LC-MS/MS) analysis, a total of 273,692 (47,151 matched) spectra were obtained. Of these spectra, 29,975 identified peptides (28,990 unique peptides) and 5250 identified proteins (4573 quantified proteins) were detected (Table 1, Table S1), and the average mass error of the peptides was <10 ppm, suggesting a high mass accuracy of the MS data (Figure 3B). In accordance with the properties of tryptic peptides, the length of the peptides was mainly distributed between 8 and 20 amino acid residues (Figure 3C), confirming that the sample preparation met the standard for further analysis. DEPs with significant change (p-value < 0.05) were selected, where the cut-off point was >1.3-fold change (p-value < 0.05) for upregulated proteins and <0.77-fold change (p-value < 0.05) for downregulated proteins. A total of 507 DEPs were identified, including 171 upregulated proteins and 336 downregulated proteins (Table S2).

2.3. Functional Classification and Subcellular Location Analysis of DEPs

To characterize the functions and subcellular locations of DEPs between the DHAV-1-infected cells (D group) and the non-infected cells (N group), GO functional classification was performed and subcellular locations were analyzed (Figure 4, Tables S3 and S4). DEPs were divided into three different groups according to the cell component, molecular function, and biological process. The cell component results showed that the DEPs were mainly distributed in cells (25%), organelles (24%), the membrane (16%), and the extracellular region (13%) (Figure 4A). Regarding molecular function, the largest two GO categories were binding and catalytic activity, accounting for 47% (330 DEPs) and 28% (193 DEPs), respectively (Figure 4B). Biological process was the category we were most concerned with; DEPs were involved in cellular processes, single-organism processes, and biological regulation, accounting for 15%, 15%, and 12%, respectively (Figure 4C). These results demonstrated that the DEPs identified in the study participated in multiple cellular components, diverse molecular functions, and varied biological processes. Subcellular location predictions showed that the DEPs were mainly distributed in the cytoplasm (35%) and the nucleus (23%) (Figure 4D), highlighting that the two organelles play significant roles in DHAV-1 infection.

2.4. Enrichment Analysis of DEPs in GO, KEGG, and Protein Domain

To gain greater insight into the preferred functional enrichment of DEPs in response to DHAV-1 infection, three types of enrichment-based analyses were performed, namely, GO, the KEGG pathway, and protein domain analysis (Figure 5). We tested the dataset for enrichment in three GO categories, i.e., cellular components, molecular function, and biological processes (Figure 5A, Tables S5 and S6). Within the cellular component category, the proteins involved in the extracellular space were upregulated while intrinsic components of the membrane-related proteins were downregulated. Molecular function analysis showed that DEPs involved in oxygen binding and signal transduction were enriched in up- and downregulated protein clusters, respectively. Regarding biological processes, both up- and downregulated proteins were distributed in various metabolic processes. Notably, proteins related to the immune response, heat shock proteins (Hsps), zinc finger proteins (ZFPs), and proteins involved in the response to biotic stimuli were significantly enriched in upregulated proteins, whereas ribosomal proteins were decreased (Table 2). We speculate that DHAV-1 infection causes the accumulation of folded proteins and then activates ER stress in DEF cells.
The KEGG pathway analysis revealed that lysosomes and phagosomes were significantly enriched, except in the immune-related pathway (Figure 5B, Tables S7 and Table S8). Combined with the results showing sequestosome 1 (SQSTM1/p62) and death-associated protein kinase (DAPK) proteins were significantly enhanced, we speculate that DHAV-1 mediates the autophagy process. Domain enrichment analysis showed that the “poly (ADP-ribose) polymerase (PARP), catalytic domain”, and “Globin-like” were enriched in upregulated proteins, whereas the “immunoglobulin-like domain” and “G protein alpha subunit, helical insertion” were extremely abundant in downregulated proteins (Figure 5C, Tables S9 and S10).

2.5. Protein Interaction Network Analysis between DEPs and Ribosome Proteins

A network analysis of protein–protein interactions was performed to investigate the important interactions between DEPs (Table S11). The results revealed tight communications among ribosome proteins (green node), immune-related proteins (U3I1B1 (IFIT5, interferon-induced protein with tetratricopeptide repeats 5), U3J8S3 (MX1, myxovirus resistance gene 1), U3J642 (IFIH1, interferon-induced helicase), and U3J6Y2 (ISG15, interferon-stimulated genes 15)) and response to stimulus proteins (U3IED8 (RSAD2, Radical S-adenosyl methionine domain containing 2), U3IRZ8 (PARP9), U3J8K7 (Notch2), and U3IJI0 (ADAM17, ADAM metallopeptidase domain 17)), and highly connected clusters in ribosome proteins (Figure 6).

2.6. Validation of TMT/MS Data Using Alternative Methods

To verify the TMT/MS data, we performed parallel reaction monitoring (PRM) analysis for eight selected DEPs (four upregulated proteins and four downregulated proteins). The relative abundance data are presented as a linear box-and whiskers plot and the expression of all identified DEPs was identical to the proteomic data (Figure 7A). The value of the D/N ratio between PRM and TMT/MS showed no noticeable differences in DEPs, except the ISG15 protein, which may be a reflection of the different detection methods. To further verify the TMT/MS data, a Western blot assay was performed and the results showed that integrin beta 2 (ITGB2), lymphocyte cytosolic protein 1 (LCP1), myosin light chain 1 (MYL1), and fibulin 2 (FBLN2) were downregulated in DHAV-1-infected cells, whereas tripartite motif containing 25 (TRIM25), tenascin C (TNC), ISG15, and Mov10 RISC complex RNA helicase (MOV10) were upregulated (Figure 7B). These findings were consistent with our TMT/MS data. Taken together, both the PRM and the Western blot results supported our TMT/MS data.

2.7. Induction of ER Stress-Induced Autophagy by DHAV-1 Infection

To further verify the speculation whether DHAV-1 infection mediates ER stress-induced autophagy, we first observed the formation of autophagosomes in DHAV-1-infected cells with rapamycin (RAPA)-treated cells and mock cells used as controls. Ultrastructural results showed that double-membraned autophagy-like vesicles formed in DHAV-1-infected cells and RAPA-treated cells but not in mock cells (Figure 8A). Notably, both the proteomic analysis and Western blot revealed that the expression levels of DAPK2 and sequestosome 1 (p62/SQSTM1) protein increased with DHAV-1 infection while their expression levels decreased with 4-phenylbutyrate (4-PBA) treatment (Figure 8B). We also found that the transformation from microtubule-associated protein 1 light chain 3-I (LC3I) to light chain 3-II (LC3II) was obviously increased in DHAV-1-infected cells (Figure 8C). These results revealed that autophagy could be initiated by DHAV-1, which might limit the degradation process.
With Hsps upregulated and ribosome proteins downregulated in our proteomic data, the expression level of the ER stress chaperone, glucose-regulated protein 78 (GRP78), was detected. Western blot results showed that the GRP78 protein was highly expressed in DHAV-1 infection, suggesting that DHAV-1 initiated the ER stress response (Figure 8D). To further verify the relationship between ER stress and autophagy that might exist in DHAV-1-infected cells, we examined the expression change in the transformation of LC3II/LC3I when cells were treated with 4-phenylbutyrate (4-PBA). The results showed that the transformation rate of LC3II/LC3I significantly increased compared with the mock cells, but it appeared to decrease in the 4-PBA-treated cells (Figure 8E).

3. Discussion

There is increasing application in proteomic analysis to identify the global changes of host proteins that are affected by viruses, stress, or abiotic factors [31,32,33]. Proteomic methods in virology can uncover antiviral responses, virus replication, and even the viral mechanisms of escape from host immune responses at the protein level. DEF cells have been used extensively to explore the molecular biological properties of DHAV-1 [34]. In this study, a quantitative proteomic analysis of DHAV-1-infected DEF cells at the whole cell protein level was performed using the TMT-labeling proteomic technique. We found 507 DEPs, among which 171 were significantly higher in abundance and 336 were in lower abundance in the D group. By referring to previous virus–host interaction analyses regarding the different virulence of DHAV-1 [22], immune-related genes, such as retinoic acid-inducible gene I (RIG-1) and melanoma differentiation associated gene 5 (MDA5), were significantly expressed in both studies. Common regulatory biological processes were also observed. For example, both studies revealed that host proteins were involved in “lysosome abundance”. The identification of these regulated proteins and processes in both studies confirmed the changes in lysosome abundance during DHAV-1 infection. Moreover, DEPs involved in ER stress and other metabolic processes were also found in our study. Proteomics analysis has previously showed changes and interactions in mitochondrial, ER, and ribosomal proteins of SH-SY5Y cells after oxygen and glucose deprivation [35]. In the present study, DEPs involved in ER stress and other metabolic processes were also identified. Further studies regarding these processes should be performed to increase our understanding of the pathogenic mechanism of DHAV-1.
It is essential in host defense that the cell maintains homeostatic processes regarding cell response to virus infection. By referring to GO enrichment, we found many biological processes and protein families that were significantly regulated and participated in cellular metabolism during DHAV-1 infection (Figure 5). It has been reported previously that the UPR regulator and ER chaperone, GRP78/BiP, was required for stress-induced autophagy [36]. A major proteomic response in the ER was identified after DHAV-1 infection. Activating transcription factor 3 (ATF3), which is involved in the external stimulus process, is integral to the double-stranded RNA-activated protein kinase (PKR)-like endoplasmic reticulum kinase (PERK)/eukaryotic initiation factor-2-α (eIF2α) signaling branch of the UPR [37]. In the present study, a major proteomic response in the ER was identified and ATF3 was upregulated in DHAV-1-infected DEF cells. Hsps, proteins that function in response to stimuli, and the ribosome protein family have been reported to be involved in cell stress, according to previous studies [38,39,40,41]. These proteins were identified as DEPs in the DHAV-1-infected DEF cells (Table 2), therefore they might be responsible for cell stress adaptation caused by DHAV-1 infection.
The acute decline of ribosome proteins is an implication that the host cell is deprived of a physiological equilibrium. The upregulation of GRP78 in DHAV-1-infected DEF cells (Figure 8) further supported the theory that ER stress was activated by DHAV-1. Ribosomal proteins, such as RPL31, RPL11, RPL38, RPL35, and RPL35A, were downregulated and formed a high interaction network (Figure 6), which may play roles in protein translation and ER docking, as per previous results [42]. Importantly, these ribosome proteins have relationships with immune-related proteins (red node in Figure 6), suggesting that the maintenance of cellular homeostasis requires complex co-regulatory interaction networks.
Autophagy, a key process in housekeeping and balance of the cell, can be initiated to maintain the synthesis, degradation, and subsequent recycling of cellular components when ER homeostasis cannot be rescued by the unfolded protein response (UPR) caused by a virus [43]. DAPK2 has been found to be a novel regulator of autophagy under stress and steady-state conditions by suppressing the mammalian target of rapamycin (mTOR) activity [44]. In the present study, both the proteomic data and Western blot assay results showed that the level of autophagy increased and DAPK2 was highly expressed in DHAV-1-infected cells. Previous studies showed that ER stress was an essential step for autophagy initiation, and pretreatment with 4-PBA decreased autophagy, illustrating that there was a close relationship between autophagy activation and ER stress [45,46]. However, whether the initiation of autophagy induced by DHAV-1 is related to ER stress is still unclear. In DHAV-1-infected cells, we observed the formation of autophagosomes, an increase in expression from LC3I to LC3II, and a decrease in expression levels of LC3 II protein when treated with 4-PBA (Figure 8). These results indicated that the DHAV-1 infection triggered the ER stress and autophagy processes, and ER stress was a key promoter in the induction of autophagy in the DHAV-1-infected DEF cells.
It was reported that the inhibition of autophagic flux may provide membranous surfaces for the assembly of viral RNA replication complexes and allow the extracellular release of viruses without cell lysis [47]. In theory, the formation of autophagosomes in DHAV-1-infected DEF cells could lead to a decrease in their degradation. However, in the present study, the less abundant proteins were significantly enriched in the lysosome and phagosome pathways, and p62 protein was more abundant, signifying that the fusion of autophagosomes with lysosomes might be limited in the DHAV-1 infection process. These results indicated that DHAV-1 may have adapted to evade autophagy to promote its infection and replication by preventing the degradation of autolysosomes, which was consistent with results regarding the coxsackievirus B3 infection [48].

4. Materials and Methods

4.1. Cells, Viruses, and Antibodies

DEF cells were obtained from duck embryos that were specifically pathogen-free for 10 days, as described previously [49]. DEF cells were cultured into Dulbecco’s modified eagle medium (DMEM) (Gibco, Carlsbad, CA, USA) with 10% fetal bovine serum (FBS, Gibco, USA) and antibiotics (100 IU/mL penicillin and 100 IU/mL streptomycin) at 37 °C under an atmosphere of 5% CO2/95% air.
DHAV-1 LY0801 strain (accession no. FJ436047) is a virulent strain, which was isolated in 2008 from an outbreak of severe duck virus hepatitis (DVH) in the Shandong province of China and kept in the veterinary molecular etiology laboratory of Shandong Agricultural University [2].
Anti-DHAV-1 monoclonal antibody (mAb) 4F8, whose epitope is “75GEILT80” in the VP1 protein of DHAV-1, was stored in our laboratory [50]. Fluorescein isothiocyanate (FITC)-labeled goat anti-mouse antibody was purchased from KPL (Gaithersburg, MD, USA). Other antibodies used in the Western blot assay in this study were purchased from Abcam (Cambridge, MA, USA).

4.2. Indirect Immunofluorescence Assay

DEF cells in 12-well plates were infected with DHAV-1. At 48 hpi, cells were washed with phosphate-buffered saline (PBS) and fixed with a mixture of acetone and formaldehyde (1:1, v/v) for 15 min at room temperature. Then, the cells were incubated with anti-DHAV-1 mAb 4F8 (dilution of 1:1000 with PBS) for 1 h at 37 °C. After being washed five times with PBS, FITC-conjugated goat anti-mouse IgG (Sigma, St. Louis, MO, USA) at a 1:200 dilution was added and incubated for 1 h at 37 °C in the dark. Finally, the plates were washed five times with PBS and observed using fluorescence microscopy (Leica AF6000, Leica, Wetzlar, Germany).

4.3. Quantitative of DHAV-1

Non-infected and DHAV-1-infected DEF cells were collected every 12 h from 12 to 72 hpi, and intracellular mRNA was extracted using Trizol reagent (Promega, Madison, WI, USA). The extracted RNA was then used for qRT-PCR quantification using the ABI 7500 (Applied Biosystems, Carlsbad, CA, USA) real-time PCR system according to the methods constructed by our lab [51]. Each experiment was replicated three times. The cycle threshold (Ct) value of each sample was used for data analysis.

4.4. Sample Preparation and TMT Labeling

DEF cells were infected with 2 MOI DHAV-1 for 2 h at 37 °C, then cells were washed with PBS and replaced with fresh DMEM with 2% FBS. At 48 hpi, cells of the D group or the N group were collected and centrifuged at 2000 rpm for 10 min to obtain the cell pellet. Each group was processed with three independent biological replications. These samples were sonicated three times on ice by using a high intensity ultrasonic processor (Scientz) in lysis buffer (8 M urea, 1% Protease Inhibitor Cocktail). After centrifugation at 12,000 g for 10 min, the remaining debris was removed, and the supernatant was collected. The concentration of proteins was determined with a bicinchoninic acid (BCA) kit (Beyotime Biotechnology, Shanghai, China), according to the manufacturer’s instructions. The extracted proteins were reduced with 5 mM dithiothreitol (DTT) at 56 °C for 30 min and alkylated with 11 mM iodoacetamide (IAA) in the darkness for 15 min at room temperature. After being diluted with 100 mM triethylammonium bicarbonate (TEAB), the proteins were digested with trypsin (Promega) at a ratio of 1:50 (trypsin/protein w/w). After trypsin digestion at 37 °C overnight, the peptide was desalted using a Strata X C18 SPE column (Phenomenex, Torrance, CA, USA) and vacuum-dried. The peptide was reconstituted in 0.5 M TEAB (Sigma, St. Louis, MI, USA) and processed according to the manufacturer’s protocol for the TMT kit (Thermo Fisher Scientific, Waltham, MA, USA). The peptide mixtures were then incubated for 2 h at room temperature, pooled, and desalted, and then dried using vacuum centrifugation.

4.5. HPLC Fractionation and LC-MS/MS Analysis

The labeled peptides were fractionated using high pH reverse-phase high pressure liquid chromatography (HPLC) with an Agilent 300Extend C18 column (5 μm particles, 4.6 mm ID, 250 mm length, Agilent, Santa Clara, USA). Briefly, peptides were first separated with a gradient of 8% to 32% acetonitrile (pH 9.0) over 60 min into 60 fractions, and the peptides were combined into 18 fractions and dried via vacuum centrifugation. Peptides were dissolved in 0.1% formic acid (solvent A, FA) and directly loaded onto a home-made reversed-phase analytical column (15-cm length, 75 μm i.d.). The gradient was comprised of an increase from 6% to 23% solvent B (0.1% FA in 98% acetonitrile) over 26 min, 23% to 35% in 8 min, and climbing to 80% in 4 min. Finally, a holding phase at 80% for 3 min was performed, all at a constant flow rate of 400 nL/min on an EASY-nLC 1000 ultraperformance liquid chromatography (UPLC) system (Thermo Fisher Scientific, Waltham, MA, USA).
The peptides were subjected to national system of innovation (NSI) sourcing, followed by MS/MS in Q ExactiveTM Plus (Thermo Fisher Scientific, Waltham, MA, USA), which was coupled online to the UPLC. Peptides were then selected for MS/MS using a normalized collisional energy (NCE) setting of 28 and the fragments were detected in the Orbitrap at a resolution of 17,500. A data-dependent procedure was performed, which involved alternating between one MS scan followed by 20 MS/MS scans with 15.0 s of dynamic exclusion. The automatic gain control (AGC) was set at 5E4 and the fixed first mass was set at 100 m/z. Three biological replicates were performed for each group.

4.6. Database Processing

The resulting MS/MS data were processed using the Maxquant search engine (v1.5.2.8). Tandem mass spectra were searched against the UniProt Anas platyrhynchos protein database (http://www.uniprot.org, 16377 sequences) and concatenated with a reverse decoy database. This target-decoy search strategy was employed to calculate the false-discovery rates (FDRs) [52]. The setting parameters of the proteome analysis were as follows: (1) The trypsin/P, trypsin; (2) 2 missing cleavages were allowed; (3) the minimum length of the peptide was 7; (4) the maximum number of modifications per peptide was 5; (5) the mass tolerance for precursor ions of the first search and the main search were 20 and 5 ppm, respectively; (6) the mass tolerance for fragment ions was 0.02 Da; (7) the fixed modification was Cys carbamidomethyl; and (8) the variable modification was Met oxidation. The FDR was adjusted to <1% and the minimum score for the peptides was set to >40.

4.7. Bioinformatics Analysis

GO annotation analysis was derived from the UniProt-GOA database (http://www.ebi.ac.uk/GOA/) [53]. The KEGG database was used to annotate the protein pathway and a two-tailed Fisher’s exact test was used to test the enrichment of the differentially expressed protein against all identified proteins [54]. The functional descriptions of the identified protein domains were annotated using InterProScan (a sequence analysis application) based on the protein sequence alignment method, and the InterPro domain database was used [55]. GO, the KEGG pathway, and protein domain enrichment were all performed with a corrected P-value of <0.05. This filtered P-value matrix was transformed via the function x = −log10 (P-value). Finally, these x values were z-transformed for each functional category and the z scores were clustered by one-way hierarchical clustering (Euclidean distance, average linkage clustering) in Genesis. Cluster membership was visualized with a heat map using the “heatmap.2” function from the “gplots” R-package. All DEP sequences were searched against the STRING database version 10.5 for protein–protein interactions (PPIs). All interactions that had a confidence score ≥0.7 (high confidence) in the STRING database were retrieved for the analysis. The interaction network obtained from STRING was visualized in Cytoscape software (version 3.0.1) [56].

4.8. PRM Analysis

The peptide samples were prepared according to the whole cell proteome analysis methodology described above. The proteins selected for PRM validation were based on the proteome results. The peptide sample was initially run in a data-dependent acquisition (DDA) mode using Q ExactiveTM Plus (Thermo Fisher Scientific, Waltham, MA, USA), which was coupled to the UPLC online with an identical gradient to subsequent PRM analysis. The electrospray voltage applied was 2.0 kV. A full mass spectrum was detected using the Orbitrap at a resolution of 35,000 (the m/z scan range was 350 to 1000), followed by 20 MS/MS scans using the Orbitrap at a resolution of 17,500 (the automatic gain control (AGC) target was 3E6 and the maximum injection time was 20 ms) in a data-independent procedure. The isolation window for MS/MS was set at 2.0 m/z. The NCE was 27% with higher energy collisional dissociation (HCD). Three biological replicates were performed.

4.9. Transmission Electron Microscopy

DEF cells were infected with DHAV-1 (MOI = 2) and maintained in 2% DMEM. At 48 hpi, cells were processed for EM according to standard procedure. Simply, cells were fixed with fixing solution for 2 h at 4 °C. The cell pellet was embedded with 1% agarose and washed three times with 0.1 M PBS (pH 7.4). Cells were post-fixed with 1% OsO4 in 0.1 M PBS (pH 7.4) at room temperature for 2 h. After rinsing in 0.1 M PBS, samples were processed via a series of dehydration steps ranging from 50% to 100% ethanol, with a final 100% acetone incubation. Samples were infiltrated with 1:1 acetone:EMBed 812 for 2 to 4 h, then with 2:1 acetone:EMBed812 overnight, and finally with pure EMBed 812 for 5 to 8 h. The resin was polymerized with embedded cells at 60 °C for 48 h. Sections measuring 80 nm were cut with ultramicrotome and stained with uranyl acetate in pure ethanol for 15 min, then rinsed with distilled water. Finally, images were observed with TEM.

4.10. Western Blotting

DEF cells were infected with DHAV-1 (2 MOI) as described above, and cell proteins were extracted at the indicated time. DEF cells were treated with 2.5 μg/mL tunicamycin (Tm), 6 μM RAPA, or 4ug/mL 4-PBA for 24 h. Extracted protein mixtures were separated using sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) gel and transferred to poly-vinylidene difluoride (PVDF) membranes using the Trans-Blot Turbo Transfer System (BIO-RAD Laboratories, Berkeley, CA, USA), which was then blocked with 5% nonfat milk. The membranes were incubated with a specific primary antibody diluted in primary antibody dilution buffer (Beyotime Biotechnology, Shanghai, China) overnight at 4 °C. After washing with tris-buffered saline with 0.1% Tween 20 (TBST), the membranes were incubated with specific secondary antibodies for an additional 2 h at room temperature. Finally, the immune complexes were visualized with an enhanced chemiluminescence (ECL) detection kit (Beyotime Biotechnology, Shanghai, China), and protein expression levels were analyzed using Image J software (Rawak Software, Inc., Munich, Germany).

5. Conclusions

In conclusion, using proteomic technology, a total of 507 DEPs were identified in DHAV-1-infected DEF cells. The proteomic data, TEM, and Western blot results indicated that DHAV-1 infection triggered ER stress and the autophagy process, and that ER stress was a key promoter in the induction of autophagy in DHAV-1-induced DEF cells. To the best of our knowledge, this was the first report regarding the utilization of proteomic technology to investigate DHAV-1–host interactions, and these findings could pave a new way toward uncovering the regulatory mechanisms behind DHAV-1 infection.

Supplementary Materials

Supplementary materials can be found at https://www.mdpi.com/1422-0067/20/24/6160/s1. Table S1. List of proteins identified and quantified by the TMT analysis in the experiment. Table S2. List of differentially expressed proteins. Table S3. GO annotation. Table S4. Subcellular location. Table S5. Enrichment analysis of DEPs based on GO annotation in upregulated proteins. Table S6. Enrichment analysis of DEPs based on GO annotation in downregulated proteins. Table S7. Enrichment analysis of the DEPs according to KEGG pathway (upregulated). Table S8. Enrichment analysis of the DEPs according to KEGG pathway (downregulated). Table S9. Enrichment of the DEPs according to domain (upregulated). Table S10. Enrichment of the DEPs according to domain (downregulated). Table S11. Interaction networks of DEPs.

Author Contributions

Conceptualization, J.L. and S.J.; data curation, R.Z. and J.W.; formal analysis, J.L., R.Z., H.Y., J.C., and S.L.; funding acquisition, S.J.; investigation, J.L., R.Z., H.Y., and W.X.; methodology, J.L., W.X., J.C., and Y.W.; project administration, S.J.; resources, H.Y. and J.W.; software, J.W., W.X., and Y.W.; supervision, J.C. and S.J.; writing—original draft, J.L.; writing—review and editing, R.Z., S.L., Z.X., and S.J.

Funding

This study was supported by grants from the National Natural Science Foundation of China (31772754); Shandong Modern Agricultural Technology & Industry System, China (SDAIT-11-15) and Funds of Shandong “Double Tops” Program, China (SYL2017YSTD11).

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

4-PBA4-phenylbutyrate
ADAM17ADAM metallopeptidase domain 17
AGCAutomatic gain control
AGCAutomatic gain control
ATF3Activating transcription factor 3
BCABicinchoninic acid
DAPK2Death-associated protein kinase 2
DDAData-dependent acquisition
DEFDuck embryo fibroblast
DEPsDifferentially expressed proteins
DHAV-1Duck hepatitis A virus type 1
DMEMDulbecco’s modified eagle medium
DTTDithiothreitol
DVHDuck virus hepatitis
ECLEnhanced chemiluminescence
EREndoplasmic reticulum
FAFormic acid
FBLN2Fibulin 2
FBSFetal bovine serum
FDRFalse discovery rates
FITCFluorescein isothiocyanate
GOGene Ontology
GRP78Glucose-regulated protein 78
HCDHigher energy collisional dissociation
HpiHours post-infection
HPCHigh pressure liquid chromatography
HspsHeat shock proteins
IAAIodoacetamide
IFIH1Interferon-induced helicase
IFIT5Interferon-induced protein with tetratricopeptide repeats 5
ISG15Interferon-stimulated gene 15
ITGB2Integrin beta 2
KEGGKyoto Encyclopedia of Genes and Genomes
LC3-IMicrotubule-associated protein 1 light chain 3-I
LC3-IIMicrotubule-associated protein 1 light chain 3-II
LCP1Lymphocyte cytosolic protein 1
mAbMonoclonal antibody
MDA5Melanoma differentiation associated gene 5
MOIMultiplicity of infection
MOV10Mov10 RISC complex RNA helicase
MS/MSTandem mass spectrometry
Mx1Myxovirus resistance gene 1
MYL1Myosin light chain 1
NCENormalized collisional energy
NSINational system of innovation
p62/SQSTM1Sequestosome 1
PARP9Poly(ADP-ribose) polymerase 9
PBSPhosphate-buffered saline
PERKDouble stranded RNA activated protein kinase (PKR)-like endoplasmic reticulum kinase
PPIProtein-protein interactions
PRMParallel reaction monitoring
PVDFPoly-vinylidene difluoride
RAPARapamycin
RIG-1Retinoic acid-inducible gene 1
RSAD2Radical S-adenosyl methionine domain containing 2
SDS-PAGESodium dodecyl sulfate-polyacrylamide gel electrophoresis
TBSTTris-buffered saline with 0.1% Tween 20
TEABTriethylammonium bicarbonate
TEMTransmission electron microscopy
TMTTandem Mass Tag
TNCTenascin C
TRIM25Tripartite motif containing 25
UPLCUltraperformance liquid chromatography
UPRUnfolded protein response
ZFPZinc finger proteins

References

  1. Kim, M.C.; Kwon, Y.K.; Joh, S.J.; Lindberg, A.M.; Kwon, J.H.; Kim, J.H.; Kim, S.J. Molecular analysis of duck hepatitis virus type 1 reveals a novel lineage close to the genus Parechovirus in the family Picornaviridae. J. Gen. Virol. 2006, 87 Pt 11, 3307–3316. [Google Scholar] [CrossRef]
  2. Gao, J.; Chen, J.; Si, X.; Xie, Z.; Zhu, Y.; Zhang, X.; Wang, S.; Jiang, S. Genetic variation of the VP1 gene of the virulent duck hepatitis A virus type 1 (DHAV-1) isolates in Shandong province of China. Virol. Sin. 2012, 27, 248–253. [Google Scholar] [CrossRef] [PubMed]
  3. Tseng, C.H.; Tsai, H.J. Molecular characterization of a new serotype of duck hepatitis virus. Virus Res. 2007, 126, 19–31. [Google Scholar] [CrossRef] [PubMed]
  4. Kim, M.C.; Kwon, Y.K.; Joh, S.J.; Kim, S.J.; Tolf, C.; Kim, J.H.; Sung, H.W.; Lindberg, A.M.; Kwon, J.H. Recent Korean isolates of duck hepatitis virus reveal the presence of a new geno- and serotype when compared to duck hepatitis virus type 1 type strains. Arch. Virol. 2007, 152, 2059–2072. [Google Scholar] [CrossRef] [PubMed]
  5. Xu, Q.; Zhang, R.; Chen, L.; Yang, L.; Li, J.; Dou, P.; Wang, H.; Xie, Z.; Wang, Y.; Jiang, S. Complete genome sequence of a duck hepatitis a virus type 3 identified in eastern China. J. Virol. 2012, 86, 13848. [Google Scholar] [CrossRef] [Green Version]
  6. Yugo, D.M.; Hauck, R.; Shivaprasad, H.L.; Meng, X.J. Hepatitis Virus Infections in Poultry. Avian Dis. 2016, 60, 576–588. [Google Scholar] [CrossRef]
  7. Song, C.; Yu, S.; Duan, Y.; Hu, Y.; Qiu, X.; Tan, L.; Sun, Y.; Wang, M.; Cheng, A.; Ding, C. Effect of age on the pathogenesis of DHV-1 in Pekin ducks and on the innate immune responses of ducks to infection. Arch. Virol. 2014, 159, 905–914. [Google Scholar] [CrossRef]
  8. Zhang, R.; Chen, J.; Zhang, J.; Yang, Y.; Li, P.; Lan, J.; Xie, Z.; Jiang, S. Novel duck hepatitis A virus type 1 isolates from adult ducks showing egg drop syndrome. Vet. Microbiol. 2018, 221, 33–37. [Google Scholar] [CrossRef]
  9. Levine, P.; Fabricant, J. A hitherto-undescribed virus disease of ducks in North America. Cornell Vet. 1950, 40, 71–86. [Google Scholar]
  10. Ou, X.; Mao, S.; Cao, J.; Cheng, A.; Wang, M.; Zhu, D.; Chen, S.; Jia, R.; Liu, M.; Sun, K.; et al. Comparative analysis of virus-host interactions caused by a virulent and an attenuated duck hepatitis A virus genotype 1. PLoS ONE 2017, 12, e0178993. [Google Scholar] [CrossRef]
  11. Mao, S.; Wang, M.; Ou, X.; Sun, D.; Cheng, A.; Zhu, D.; Chen, S.; Jia, R.; Liu, M.; Sun, K.; et al. Virologic and Immunologic Characteristics in Mature Ducks with Acute Duck Hepatitis A Virus 1 Infection. Front. Immunol. 2017, 8, 1574. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Song, C.; Liao, Y.; Gao, W.; Yu, S.; Sun, Y.; Qiu, X.; Tan, L.; Cheng, A.; Wang, M.; Ma, Z.; et al. Virulent and attenuated strains of duck hepatitis A virus elicit discordant innate immune responses in vivo. J. Gen. Virol. 2014, 95 Pt 12, 2716–2726. [Google Scholar] [CrossRef]
  13. Lipatova, Z.; Segev, N. A Role for Macro-ER-Phagy in ER Quality Control. PLoS Genet. 2015, 11, e1005390. [Google Scholar] [CrossRef] [PubMed]
  14. Bernales, S.; Schuck, S.; Walter, P. ER-phagy: Selective autophagy of the endoplasmic reticulum. Autophagy 2007, 3, 285–287. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Jheng, J.R.; Ho, J.Y.; Horng, J.T. ER stress, autophagy, and RNA viruses. Front. Microbiol. 2014, 5, 388. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Deretic, V.; Levine, B. Autophagy, immunity, and microbial adaptations. Cell Host Microbe 2009, 5, 527–549. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Levine, B.; Deretic, V. Unveiling the roles of autophagy in innate and adaptive immunity. Nat. Rev. Immunol. 2007, 7, 767–777. [Google Scholar] [CrossRef]
  18. Deretic, V. Multiple regulatory and effector roles of autophagy in immunity. Curr. Opin. Immunol. 2009, 21, 53–62. [Google Scholar] [CrossRef] [Green Version]
  19. Ait-Goughoulte, M.; Kanda, T.; Meyer, K.; Ryerse, J.S.; Ray, R.B.; Ray, R. Hepatitis C virus genotype 1a growth and induction of autophagy. J. Virol. 2008, 82, 2241–2249. [Google Scholar] [CrossRef] [Green Version]
  20. Orvedahl, A.; Levine, B. Viral evasion of autophagy. Autophagy 2008, 4, 280–285. [Google Scholar] [CrossRef] [Green Version]
  21. Sir, D.; Chen, W.L.; Choi, J.; Wakita, T.; Yen, T.S.; Ou, J.H. Induction of incomplete autophagic response by hepatitis C virus via the unfolded protein response. Hepatology 2008, 48, 1054–1061. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Xie, J.; Zeng, Q.; Wang, M.; Ou, X.; Ma, Y.; Cheng, A.; Zhao, X.X.; Liu, M.; Zhu, D.; Chen, S.; et al. Transcriptomic Characterization of a Chicken Embryo Model Infected with Duck Hepatitis A Virus Type 1. Front. Immunol. 2018, 9, 1845. [Google Scholar] [CrossRef] [PubMed]
  23. Sheng, X.D.; Zhang, W.P.; Zhang, Q.R.; Gu, C.Q.; Hu, X.Y.; Cheng, G.F. Apoptosis induction in duck tissues during duck hepatitis A virus type 1 infection. Poult. Sci. 2014, 93, 527–534. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Zhang, K.; Zheng, S.; Yang, J.S.; Chen, Y.; Cheng, Z. Comprehensive profiling of protein lysine acetylation in Escherichia coli. J. Proteome Res. 2013, 12, 844–851. [Google Scholar] [CrossRef] [PubMed]
  25. Brewis, I.A.; Brennan, P. Proteomics technologies for the global identification and quantification of proteins. Adv. Protein Chem. Struct. Biol. 2010, 80, 1–44. [Google Scholar] [PubMed]
  26. Wu, W.W.; Wang, G.; Baek, S.J.; Shen, R.F. Comparative study of three proteomic quantitative methods, DIGE, cICAT, and iTRAQ, using 2D gel- or LC-MALDI TOF/TOF. J. Proteome Res. 2006, 5, 651–658. [Google Scholar] [CrossRef]
  27. Liu, F.; Zheng, H.; Tong, W.; Li, G.X.; Tian, Q.; Liang, C.; Li, L.W.; Zheng, X.C.; Tong, G.Z. Identification and Analysis of Novel Viral and Host Dysregulated MicroRNAs in Variant Pseudorabies Virus-Infected PK15 Cells. PLoS ONE 2016, 11, e0151546. [Google Scholar] [CrossRef]
  28. Das, P.P.; Lin, Q.; Wong, S.M. Comparative proteomics of Tobacco mosaic virus-infected Nicotiana tabacum plants identified major host proteins involved in photosystems and plant defence. J. Proteom. 2019, 194, 191. [Google Scholar] [CrossRef]
  29. Yang, S.; Pei, Y.; Zhao, A. iTRAQ-based Proteomic Analysis of Porcine Kidney Epithelial PK15 cells Infected with Pseudorabies virus. Sci. Rep. 2017, 7, 45922. [Google Scholar] [CrossRef] [Green Version]
  30. Wang, Z.; Zhao, K.; Pan, Y.; Wang, J.; Song, X.; Ge, W.; Yuan, M.; Lei, T.; Wang, L.; Zhang, L.; et al. Genomic, expressional, protein-protein interactional analysis of Trihelix transcription factor genes in Setaria italia and inference of their evolutionary trajectory. BMC Genom. 2018, 19, 665. [Google Scholar] [CrossRef] [Green Version]
  31. Di Carli, M.; Benvenuto, E.; Donini, M. Recent insights into plant-virus interactions through proteomic analysis. J. Proteome Res. 2012, 11, 4765–4780. [Google Scholar] [CrossRef] [PubMed]
  32. Zandalinas, S.I.; Mittler, R.; Balfagon, D.; Arbona, V.; Gomez-Cadenas, A. Plant adaptations to the combination of drought and high temperatures. Physiol. Plant. 2018, 162, 2–12. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Perez-Clemente, R.M.; Vives, V.; Zandalinas, S.I.; Lopez-Climent, M.F.; Munoz, V.; Gomez-Cadenas, A. Biotechnological approaches to study plant responses to stress. BioMed Res. Int. 2013, 2013, 654120. [Google Scholar] [CrossRef] [PubMed]
  34. Cao, J.; Ou, X.; Zhu, D.; Ma, G.; Cheng, A.; Wang, M.; Chen, S.; Jia, R.; Liu, M.; Sun, K.; et al. The 2A2 protein of Duck hepatitis A virus type 1 induces apoptosis in primary cell culture. Virus Genes 2016, 52, 780–788. [Google Scholar] [CrossRef] [PubMed]
  35. Herrmann, A.G.; Deighton, R.F.; Le Bihan, T.; McCulloch, M.C.; Searcy, J.L.; Kerr, L.E.; McCulloch, J. Adaptive changes in the neuronal proteome: Mitochondrial energy production, endoplasmic reticulum stress, and ribosomal dysfunction in the cellular response to metabolic stress. J. Cerebr. Blood Flow Metab. 2013, 33, 673–683. [Google Scholar] [CrossRef] [Green Version]
  36. Li, J.; Ni, M.; Lee, B.; Barron, E.; Hinton, D.R.; Lee, A.S. The unfolded protein response regulator GRP78/BiP is required for endoplasmic reticulum integrity and stress-induced autophagy in mammalian cells. Cell Death Differ. 2008, 15, 1460–1471. [Google Scholar] [CrossRef]
  37. Jiang, H.Y.; Wek, S.A.; McGrath, B.C.; Lu, D.; Hai, T.; Harding, H.P.; Wang, X.; Ron, D.; Cavener, D.R.; Wek, R.C. Activating transcription factor 3 is integral to the eukaryotic initiation factor 2 kinase stress response. Mol. Cell. Biol. 2004, 24, 1365–1377. [Google Scholar] [CrossRef] [Green Version]
  38. Haslbeck, M.; Vierling, E. A first line of stress defense: Small heat shock proteins and their function in protein homeostasis. J. Mol. Biol. 2015, 427, 1537–1548. [Google Scholar] [CrossRef] [Green Version]
  39. Arrigo, A.P. The cellular “networking” of mammalian Hsp27 and its functions in the control of protein folding, redox state and apoptosis. Adv. Exp. Med. Biol. 2007, 594, 14–26. [Google Scholar]
  40. Boulon, S.; Westman, B.J.; Hutten, S.; Boisvert, F.M.; Lamond, A.I. The nucleolus under stress. Mol. Cell 2010, 40, 216–227. [Google Scholar] [CrossRef]
  41. Bai, P. Biology of Poly (ADP-Ribose) Polymerases: The Factotums of Cell Maintenance. Mol. Cell 2015, 58, 947–958. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Li, C.; Guo, Z.; Zhao, R.; Sun, W.; Xie, M. Proteomic Analysis of Liver Proteins in a Rat Model of Chronic Restraint Stress-Induced Depression. BioMed Res. Int. 2017, 2017, 7508316. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Teng, R.J.; Du, J.; Welak, S.; Guan, T.; Eis, A.; Shi, Y.; Konduri, G.G. Cross talk between NADPH oxidase and autophagy in pulmonary artery endothelial cells with intrauterine persistent pulmonary hypertension. AM. J. Physiol.-Lung Cell. Mol. Physiol. 2012, 302, L651–L663. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Ber, Y.; Shiloh, R.; Gilad, Y.; Degani, N.; Bialik, S.; Kimchi, A. DAPK2 is a novel regulator of mTORC1 activity and autophagy. Cell Death Differ. 2015, 22, 465–475. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Kouroku, Y.; Fujita, E.; Tanida, I.; Ueno, T.; Isoai, A.; Kumagai, H.; Ogawa, S.; Kaufman, R.J.; Kominami, E.; Momoi, T. ER stress (PERK/eIF2alpha phosphorylation) mediates the polyglutamine-induced LC3 conversion, an essential step for autophagy formation. Cell Death Differ. 2007, 14, 230–239. [Google Scholar] [CrossRef] [Green Version]
  46. Qin, L.; Wang, Z.; Tao, L.; Wang, Y. ER stress negatively regulates AKT/TSC/mTOR pathway to enhance autophagy. Autophagy 2010, 6, 239–247. [Google Scholar] [CrossRef] [Green Version]
  47. Taylor, M.P.; Kirkegaard, K. Potential subversion of autophagosomal pathway by picornaviruses. Autophagy 2008, 4, 286–289. [Google Scholar] [CrossRef] [Green Version]
  48. Kemball, C.C.; Alirezaei, M.; Flynn, C.T.; Wood, M.R.; Harkins, S.; Kiosses, W.B.; Whitton, J.L. Coxsackievirus infection induces autophagy-like vesicles and megaphagosomes in pancreatic acinar cells in vivo. J. Virol. 2010, 84, 12110–12124. [Google Scholar] [CrossRef] [Green Version]
  49. Jacolot, S.; Ferec, C.; Mura, C. Iron responses in hepatic, intestinal and macrophage/monocyte cell lines under different culture conditions. Blood Cell. Mol. Dis. 2008, 41, 100–108. [Google Scholar] [CrossRef]
  50. Zhang, R.; Zhou, G.; Xin, Y.; Chen, J.; Lin, S.; Tian, Y.; Xie, Z.; Jiang, S. Identification of a conserved neutralizing linear B-cell epitope in the VP1 proteins of duck hepatitis A virus type 1 and 3. Vet. Microbiol. 2015, 180, 196–204. [Google Scholar] [CrossRef]
  51. Lin, S.L.; Cong, R.C.; Zhang, R.H.; Chen, J.H.; Xia, L.L.; Xie, Z.J.; Wang, Y.; Zhu, Y.L.; Jiang, S.J. Circulation and in vivo distribution of duck hepatitis A virus types 1 and 3 in infected ducklings. Arch. Virol. 2016, 161, 405–416. [Google Scholar] [CrossRef] [PubMed]
  52. Tyanova, S.; Temu, T.; Cox, J. The MaxQuant computational platform for mass spectrometry-based shotgun proteomics. Nat. Protoc. 2016, 11, 2301–2319. [Google Scholar] [CrossRef] [PubMed]
  53. Barrell, D.; Dimmer, E.; Huntley, R.P.; Binns, D.; O’Donovan, C.; Apweiler, R. The GOA database in 2009—An integrated Gene Ontology Annotation resource. Nucleic Acids Res. 2009, 37, D396–D403. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Kanehisa, M.; Furumichi, M.; Tanabe, M.; Sato, Y.; Morishima, K. KEGG: New perspectives on genomes, pathways, diseases and drugs. Nucleic Acids Res. 2017, 45, D353–D361. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Jones, P.; Binns, D.; Chang, H.Y.; Fraser, M.; Li, W.; McAnulla, C.; McWilliam, H.; Maslen, J.; Mitchell, A.; Nuka, G.; et al. InterProScan 5: Genome-scale protein function classification. Bioinformatics 2014, 30, 1236–1240. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Shannon, P.; Markiel, A.; Ozier, O.; Baliga, N.S.; Wang, J.T.; Ramage, D.; Amin, N.; Schwikowski, B.; Ideker, T. Cytoscape: A software environment for integrated models of biomolecular interaction networks. Genome Res. 2003, 13, 2498–2504. [Google Scholar] [CrossRef] [PubMed]
Ijms 20 06160 g001 550
Figure 1. Characteristics of duck hepatitis A virus 1 (DHAV-1) infection on duck embryo fibroblast (DEF) cells. (A) Cell morphology of DEF cells with or without DHAV-1 infection were observed at 24, 48, 60, and 72 hpi. (B) DHAV-1 successfully infected DEF cells. (C) The mRNA copy numbers of DHAV-1 on DEF cells with different multiplicity of infection (MOI) at 12, 24, 36, 48, and 60 hpi. (D) The mRNA copy numbers of DHAV-1 on DEF cells with a MOI of 2 at 12, 24, 36, 48, 60, and 72 hpi.
Figure 1. Characteristics of duck hepatitis A virus 1 (DHAV-1) infection on duck embryo fibroblast (DEF) cells. (A) Cell morphology of DEF cells with or without DHAV-1 infection were observed at 24, 48, 60, and 72 hpi. (B) DHAV-1 successfully infected DEF cells. (C) The mRNA copy numbers of DHAV-1 on DEF cells with different multiplicity of infection (MOI) at 12, 24, 36, 48, and 60 hpi. (D) The mRNA copy numbers of DHAV-1 on DEF cells with a MOI of 2 at 12, 24, 36, 48, 60, and 72 hpi.
Ijms 20 06160 g001
Ijms 20 06160 g002 550
Figure 2. A schematic diagram of the experimental design flow. Proteins were extracted from DHAV-1-infected DEF cells and non-infected DEF cells. N: Non-infected DEF cells; D: DHAV-1-infected DEF cells. Experiments were performed as biological triplicates.
Figure 2. A schematic diagram of the experimental design flow. Proteins were extracted from DHAV-1-infected DEF cells and non-infected DEF cells. N: Non-infected DEF cells; D: DHAV-1-infected DEF cells. Experiments were performed as biological triplicates.
Ijms 20 06160 g002
Ijms 20 06160 g003 550
Figure 3. Quality control validation of mass spectrometry (MS) data. (A) Mass data of all identified peptides. It is a negative correlation or positive correlation when the Pearson coefficient is closer to -1 or 1, respectively. The green indicates the Pearson coefficient closer to -1, while the red indicates the Pearson coefficient closer to 1. (B) Average peptide mass error. (C) Length distribution of all identified peptides.
Figure 3. Quality control validation of mass spectrometry (MS) data. (A) Mass data of all identified peptides. It is a negative correlation or positive correlation when the Pearson coefficient is closer to -1 or 1, respectively. The green indicates the Pearson coefficient closer to -1, while the red indicates the Pearson coefficient closer to 1. (B) Average peptide mass error. (C) Length distribution of all identified peptides.
Ijms 20 06160 g003
Ijms 20 06160 g004a 550Ijms 20 06160 g004b 550
Figure 4. Gene ontology (GO) functional classification and subcellular functional annotation of differentially expressed proteins (DEPs). (A) GO annotation in terms of cellular components. (B) GO annotation in terms of molecular function. (C) GO annotation in terms of biological processes. (D) Subcellular locations of DEPs. GO: Gene ontology.
Figure 4. Gene ontology (GO) functional classification and subcellular functional annotation of differentially expressed proteins (DEPs). (A) GO annotation in terms of cellular components. (B) GO annotation in terms of molecular function. (C) GO annotation in terms of biological processes. (D) Subcellular locations of DEPs. GO: Gene ontology.
Ijms 20 06160 g004aIjms 20 06160 g004b
Ijms 20 06160 g005a 550Ijms 20 06160 g005b 550
Figure 5. Enrichment clustering analyses of DEPs. (A) DEPs were classified by GO enrichment based on the following three categories: biological processes, cellular components, and molecular function. (B) DEPs were annotated based on the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway database. (C) DEPs were annotated based on the protein domain database.
Figure 5. Enrichment clustering analyses of DEPs. (A) DEPs were classified by GO enrichment based on the following three categories: biological processes, cellular components, and molecular function. (B) DEPs were annotated based on the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway database. (C) DEPs were annotated based on the protein domain database.
Ijms 20 06160 g005aIjms 20 06160 g005b
Ijms 20 06160 g006 550
Figure 6. Interaction network between DEPs and ribosome proteins when comparing the D group to the N group. Red nodes represent upregulated DEPs while green nodes represent downregulated DEPs.
Figure 6. Interaction network between DEPs and ribosome proteins when comparing the D group to the N group. Red nodes represent upregulated DEPs while green nodes represent downregulated DEPs.
Ijms 20 06160 g006
Ijms 20 06160 g007 550
Figure 7. Verification of the DEPs by parallel reaction monitoring (PRM) and Western blot. (A) DEPs, including integrin beta 2 (ITGB2), lymphocyte cytosolic protein 1(LCP1), myosin light chain 1 (MYL1), fibulin 2 (FBLN2), tripartite motif containing 25 (TRIM25), tenascin C (TNC), interferon-stimulated genes 15 (ISG15), and Mov10 RISC complex RNA helicase (MOV10), in DHAV-1-infected and non-infected DEF cells were verified using PRM analysis. (B) Western blot analysis was used to reconfirm the expression level of the above-mentioned proteins. The expression trends of the DEPs that were identified in Western blot analysis matched those observed in the proteomics data.
Figure 7. Verification of the DEPs by parallel reaction monitoring (PRM) and Western blot. (A) DEPs, including integrin beta 2 (ITGB2), lymphocyte cytosolic protein 1(LCP1), myosin light chain 1 (MYL1), fibulin 2 (FBLN2), tripartite motif containing 25 (TRIM25), tenascin C (TNC), interferon-stimulated genes 15 (ISG15), and Mov10 RISC complex RNA helicase (MOV10), in DHAV-1-infected and non-infected DEF cells were verified using PRM analysis. (B) Western blot analysis was used to reconfirm the expression level of the above-mentioned proteins. The expression trends of the DEPs that were identified in Western blot analysis matched those observed in the proteomics data.
Ijms 20 06160 g007
Ijms 20 06160 g008 550
Figure 8. Mediation of ER stress-induced autophagy by DHAV-1. (A) The formation of double-membraned autophagy-like vesicles in DEF cells with mock infection, DHAV-1 infection, and rapamycin (RAPA) treatment was observed. (B) Expression levels of DAPK2 and p62/SQSTM1 were identified using Western blot. (C) Western blot detected the transformation from microtubule-associated protein 1 light chain 3-I (LC3I) to light chain 3-II (LC3II) caused by DHAV-1 infection in DEF cells. (D) Expression level analysis of the ER stress marker protein glucose-regulated protein 78 (GRP78) in DHAV-1-infected DEF cells. (E) Inhibiting ER stress could decrease the transformation from LC3I to LC3II caused by DHAV-1 infection.
Figure 8. Mediation of ER stress-induced autophagy by DHAV-1. (A) The formation of double-membraned autophagy-like vesicles in DEF cells with mock infection, DHAV-1 infection, and rapamycin (RAPA) treatment was observed. (B) Expression levels of DAPK2 and p62/SQSTM1 were identified using Western blot. (C) Western blot detected the transformation from microtubule-associated protein 1 light chain 3-I (LC3I) to light chain 3-II (LC3II) caused by DHAV-1 infection in DEF cells. (D) Expression level analysis of the ER stress marker protein glucose-regulated protein 78 (GRP78) in DHAV-1-infected DEF cells. (E) Inhibiting ER stress could decrease the transformation from LC3I to LC3II caused by DHAV-1 infection.
Ijms 20 06160 g008
Table
Table 1. Summary of MS/MS spectra database search analysis.
Table 1. Summary of MS/MS spectra database search analysis.
Total SpectraMatched SpectraPeptidesUnique PeptidesIdentified ProteinsQuantifiable Proteins
273,69247,151 (17.2%)29,97528,99052504573
Table
Table 2. A partial result of the differentially expressed proteins of the DHAV-1-infected group (D) compared to the non-infected group (N).
Table 2. A partial result of the differentially expressed proteins of the DHAV-1-infected group (D) compared to the non-infected group (N).
Protein IDDescriptionD/N RatioP-Value
Heat Shock Proteins
R0K012Activator of 90 kDa heat shock protein ATPase-like protein 21.350.0086
U3IGK2DnaJ heat shock protein family (Hsp40) member C211.3550.00424
U3IPI8Heat shock protein family B (small) member 81.7730.0000227
U3IZQ3Heat shock protein family B (small) member 72.0340.000119
Ribosome Related Proteins
U3IZF8Ribosomal protein L110.6910.000158
U3IPM7Ribosomal protein L380.7320.000925
U3J6F1Ribosomal protein L35a0.7510.000819
R0L8K360S ribosomal protein L350.6430.000863
U3JA05Ribosomal protein L310.740.0000179
U3IU4460S ribosomal protein L360.5160.000141
U3IA42Mitochondrial ribosomal protein L270.7650.0119
U3IK48Ribosomal protein S50.6580.000638
U3J834Ribosomal protein S190.7570.000105
Zinc Finger Proteins
U3IIQ5Zinc finger E-box binding homeobox 11.3450.0414
U3IUP4NFX1-type zinc finger-containing protein 12.6750.0000234
Response to Stimulus
U3IWA4Apolipoprotein B1.3880.00256
U3ITV8Interferon induced with helicase C domain 11.9860.0000433
U3I515eukaryotic translation initiation factor 2 alpha kinase 23.4170.00000138
U3IQI3F-box protein 182.4680.0000779
S4SM19ATP-dependent RNA helicase2.8060.0000825
U3J4N1SAM and HD domain containing deoxynucleoside triphosphate triphosphohydrolase 11.8280.0000226
U3IQ81Adenosine deaminase, RNA specific1.670.0000588
U3J485ATP binding cassette subfamily A member 11.4140.001
U3IRD3Ankyrin repeat domain 11.5670.00122
U3I3I3Beta-2-microglobulin2.7180.0000783
U3IED8Radical S-adenosyl methionine domain containing 24.4540.0000988
R0LWV1Activating transcription factor 31.7370.000138
Immune-Related Proteins
U3IXS5Tripartite motif containing 351.50.0000634
U3I1L8Tripartite motif containing 591.380.000201
U3J0V8Signal transducer and activator of transcription 12.9850.0000419
U3IB66Tripartite motif containing 253.4690.0000000192
U3I8M2Interferon induced protein 352.7160.000245
U3ITV8Interferon induced with helicase C domain 11.9860.0000433
U3I515eukaryotic translation initiation factor 2 alpha kinase 23.4170.00000138
U3IRN3S100 calcium binding protein A122.6870.0000754
U3IGN0N-myc and STAT interactor1.4080.000302
U3J6Y2ISG15 ubiquitin-like modifier9.8580.00000348
U3I1B1IFN-induced protein with tetratricopeptide repeats 58.5060.00104
Autophagy-Associated Proteins
U3I5T9V-type proton ATPase subunit0.4930.0000955
U3IIS7Integrin beta0.5830.00000173
U3IN31N-acetylglucosamine-6-sulfatase0.7320.000137
U3IYU5Cathepsin S0.630.000337
U3I5T9V-type proton ATPase subunit D0.4930.0000955
U3I832V-type proton ATPase subunit H0.7690.00176
U3IST9V-type proton ATPase subunit C0.5410.00000192
U3IMC2Sequestosome 1/SQSTM11.3360.000459
U3I4S4Death associated protein kinase 21.3740.000564

Share and Cite

MDPI and ACS Style

Lan, J.; Zhang, R.; Yu, H.; Wang, J.; Xue, W.; Chen, J.; Lin, S.; Wang, Y.; Xie, Z.; Jiang, S. Quantitative Proteomic Analysis Uncovers the Mediation of Endoplasmic Reticulum Stress-Induced Autophagy in DHAV-1-Infected DEF Cells. Int. J. Mol. Sci. 2019, 20, 6160. https://doi.org/10.3390/ijms20246160

AMA Style

Lan J, Zhang R, Yu H, Wang J, Xue W, Chen J, Lin S, Wang Y, Xie Z, Jiang S. Quantitative Proteomic Analysis Uncovers the Mediation of Endoplasmic Reticulum Stress-Induced Autophagy in DHAV-1-Infected DEF Cells. International Journal of Molecular Sciences. 2019; 20(24):6160. https://doi.org/10.3390/ijms20246160

Chicago/Turabian Style

Lan, Jingjing, Ruihua Zhang, Honglei Yu, Jingyu Wang, Wenxiang Xue, Junhao Chen, Shaoli Lin, Yu Wang, Zhijing Xie, and Shijin Jiang. 2019. "Quantitative Proteomic Analysis Uncovers the Mediation of Endoplasmic Reticulum Stress-Induced Autophagy in DHAV-1-Infected DEF Cells" International Journal of Molecular Sciences 20, no. 24: 6160. https://doi.org/10.3390/ijms20246160

APA Style

Lan, J., Zhang, R., Yu, H., Wang, J., Xue, W., Chen, J., Lin, S., Wang, Y., Xie, Z., & Jiang, S. (2019). Quantitative Proteomic Analysis Uncovers the Mediation of Endoplasmic Reticulum Stress-Induced Autophagy in DHAV-1-Infected DEF Cells. International Journal of Molecular Sciences, 20(24), 6160. https://doi.org/10.3390/ijms20246160

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop