Skip to main content

Proteomics analysis of differentially expressed proteins in chicken trachea and kidney after infection with the highly virulent and attenuated coronavirus infectious bronchitis virus in vivo

Abstract

Background

Infectious bronchitis virus (IBV) is first to be discovered coronavirus which is probably endemic in all regions with intensive impact on poultry production. In this study, we used two-dimensional gel electrophoresis (2-DE) and two-dimensional fluorescence difference gel electrophoresis (2-DIGE), coupled with matrix-assisted laser desorption/ionization time-of-flight tandem mass spectrometry (MALDI-TOF/TOF-MS), to explore the global proteome profiles of trachea and kidney tissues from chicken at different stages infected in vivo with the highly virulent ck/CH/LDL/97I P5 strain of infectious bronchitis virus (IBV) and the embryo-passaged, attenuated ck/CH/LDL/97I P115 strain.

Results

Fifty-eight differentially expressed proteins were identified. Results demonstrated that some proteins which had functions in cytoskeleton organization, anti-oxidative stress, and stress response, showed different change patterns in abundance from chicken infected with the highly virulent ck/CH/LDL/97I P5 strain and those given the embryo-passaged, attenuated P115 stain. In addition, the dynamic transcriptional alterations of 12 selected proteins were analyzed by the real-time RT-PCR, and western blot analysis confirmed the change in abundance of heat shock proteins (HSP) beta-1, annexin A2, and annexin A5.

Conclusions

The proteomic alterations described here may suggest that these changes to protein expression correlate with IBV virus' virulence in chicken, hence provides valuable insights into the interactions of IBV with its host and may also assist with investigations of the pathogenesis of IBV and other coronavirus infections.

Background

Coronaviruses (CoVs) are enveloped single-stranded positive sense RNA viruses that belong to the family Coronaviridae in the order Nidovirales. They are able to infect humans as well as other animals, including cows, pigs, mice, and chickens, they generally cause respiratory infection, gastrointestinal, and neurological disorders of varying severity. Infectious bronchitis virus (IBV) was the first coronavirus to be discovered, and is classed among the Gamma coronaviruses on the basis of antigenic and genetic relatedness [1]. It is a major poultry pathogen and is probably endemic in all chicken-raising regions; it has a severe impact on poultry production, causing heavy economic losses. All strains of IBV are capable of infecting a large range of epithelial surfaces of chickens, such as those of the trachea, kidney, oviduct and proventriculus [2].

Coronavirus infection has dramatic effects on host cell morphology, transcription and translation patterns, the cell cycle, cytoskeleton, suppression of interferon, and apoptosis pathways. Coronavirus infection may also cause inflammation, alter the immune and stress responses, and modify the coagulation pathways [3]. Such profound functional and morphological changes in host cells are associated with significant changes in the patterns of expression of host cell genes. Several studies have described changes in host gene expression associated with coronavirus infection, as documented by microarray technologies [4–9]. But ultimately, protein expression and post-translational modification (PTM) determine virus replication. Furthermore, transcriptome analyses only provide a snapshot of gene expression patterns; they also suffer from several limitations, including inconsistencies with the levels of expression of the corresponding proteins as well as lacking the ability to provide information on PTM. Approaches that use proteomics are promising because they can circumvent some of the issues associated with transcriptomics approaches [10]. More recently, comparative proteomics analysis has emerged as a valuable tool for the establishment of the global host protein profile in response to virus infection. It has been used to study enveloped RNA viruses such as influenza virus, respiratory syncytial virus (RSV), parainfluenza virus (PIV), human metapneumovirus (hMPV), SARS-CoV, and mouse hepatitis virus (MHV) [11–18]. It provides invaluable information on the cellular signaling pathways involved in either the cellular response to viral infections, or the viral manipulation of cellular machinery to ensure their own survival. For IBV, to the best of our knowledge, only some recent studies have investigated the changes in the expression of cellular proteins during IBV infection in ex vivo or in ovo [19–21]. However, the in vivo infection model could yield more biologically relevant insights into pathogenesis.

In this study, we used two-dimensional gel electrophoresis (2-DE) and two-dimensional fluorescence difference gel electrophoresis (2-DIGE), coupled with matrix-assisted laser desorption/ionization time-of-flight tandem mass spectrometry (MALDI-TOF/TOF-MS), to explore global changed proteome profiles of trachea and kidney tissues from chicken at different stages infected in vivo with the highly virulent ck/CH/LDL/97I P5 strain of IBV and an embryo-passaged strain of attenuated virulence, ck/CH/LDL/97I P115. In total, 58 differentially expressed proteins were identified and classified into several functional categories, including cytoskeleton organization, anti-oxidative stress, the stress response, acute phase response, and energy metabolism. In addition, the dynamic transcriptional alterations of 12 selected proteins were analyzed by the real-time RT-PCR method. Simultaneously, western blot analysis confirmed the change in abundance of the heat shock proteins (HSP) beta-1, annexin A2, and annexin A5. The potential roles of some of these identified proteins are discussed in order to characterize their potential functional roles during IBV infection in vivo. These results provide valuable insights into the interactions of IBV with its host, and may also be useful in investigations of the pathogenesis of IBV and other coronaviruses.

Results

IBV antibody detection and observed clinical signs

All chickens exhibited respiratory clinical signs at about 4-14 dpi with the IBV ck/CH/LDL/97I P5 strain. The clinical signs included tracheal rales, watery eyes, nasal mucus, and sneezing. The clinical signs shown by the inoculated birds tended to disappear gradually after 14 dpi. Gross lesions of the chickens killed in the P5-infected group were confined mainly to the kidneys. The kidney parenchyma of the dead birds was pale, swollen and mottled; the tubules and urethras were distended with uric acid crystals [22]. For chickens in the control group and those inoculated with the IBV ck/CH/LDL/97I P115 strain, no respiratory clinical signs and no gross lesions were observed during the experimental period.

As summarized in Additional file 1, no chickens inoculated with the ck/CH/LDL/97I P5 or the ck/CH/LDL/97I P115 strain of IBV showed seroconversion at 4 dpi. Antibodies appeared at 7 dpi and all of the chickens showed seroconversion after 14 dpi with both IBV strains. The chickens in the non-inoculated control group showed a negative serum antibody response.

Analysis of viral load in the trachea and kidney of IBV-infected chickens

Successful IBV infection was also verified using real-time RT-PCR. Results are presented in Figure 1. Virus was not detected from the trachea and kidney of chickens from control group. In the P5-infected group, virus was detected from trachea and kidney at 4, 7, 14, and 21 dpi, the peak of viral copy number was reached at 4 dpi, after which time viral load fell. In trachea of P115-infected group, the peak of viral copy number was at 4 dpi, then fell until 14 dpi, and a little fluctuated at 21 dpi. In kidney of P115-infected group, the peak of viral copy number was also at 4 dpi, then fell at 7 dpi, and only a little fluctuated at 14 and 21 dpi. Furthermore, both in trachea and kidney, the P115-infected group had the lower viral genome copies than P5-infected group.

Figure 1
figure 1

The viral load in each sample was quantified using real-time RT-PCR. The average viral copy number (VCN) per g tissue of each group was calculated. Error bars indicate standard error of the mean, and dpi represent days post-inoculation.

Analysis of differentially changed proteins in abundance by 2-DE and 2-DIGE

The global protein changes in the chicken trachea and kidney tissues at different stages after infection with the IBV ck/CH/LDL/97I P5 and ck/CH/LDL/97I P115 strain were investigated. For the tracheal protein samples, Figure 2 shows representative images of the tracheal samples: 1366 ± 39, 1536 ± 126, and 1600 ± 167 protein spots were detected in gels from the control group, P5-infected group, and P115-infected group at 4 dpi; 1355 ± 300, 1518 ± 175, and 1078 ± 122 protein spots were detected in gels from the control group, P5-infected group, and P115-infected group at 7 dpi; 1293 ± 91, 1365 ± 126, and 1220 ± 56 protein spots were detected in gels from the control group, P5-infected group, and P115-infected group at 14 dpi; 1204 ± 91, 1236 ± 42, and 1111 ± 50 protein spots were detected in gels from the control group, P5-infected group, and P115-infected group at 21 dpi. For the kidney protein samples, Figure 3 shows representative images of the kidney samples at 4, 7, 14, and 21 dpi: 2315 ± 87, 2482 ± 189, 2607 ± 238, and 2593 ± 192 protein spots were detected, respectively.

Figure 2
figure 2

Representative 2-DE maps of tracheal tissues from the IBV ck/CH/LDL/97I P 5 -infected group, P 115 -infected group and control group. (A) 4 dpi; (B) 7 dpi; (C) 14 dpi; (D) 21 dpi. Protein samples were separated on 13 cm pH 3-10 linear IPG strips, followed by SDS-PAGE, and stained with Coomassie Blue R-350. The images were analyzed using Image Master 2D Platinum 6.0 software. The differentially expressed protein spots identified were marked with circles and labeled with the respective Match ID listed in Tables 1 and 2.

Figure 3
figure 3

Representative 2-DIGE maps of kidney tissues from the IBV ck/CH/LDL/97I P 5 -infected group, P 115 -infected group and control group. (A) 4 dpi; (B) 7 dpi; (C) 14 dpi; (D) 21 dpi. Protein samples labeled with fluor were separated on 24 cm pH 3-10 linear IPG strips, followed by SDS-PAGE, gels were scanned on a Typhoon 9400 scanner, and image analysis was performed with Ettanâ„¢ DeCyder Software version v6.5. The protein spots identified were marked with circles and labeled with the respective Match ID listed in Tables 3 and 4.

The number of differentially expressed protein spots in the P5 and P115 infected groups is summarized in Additional file 2. In the tracheal total proteome at 4 dpi, the host response to infection with both P5 and P115 strains appears to involve predominantly increase of proteins abundance. By 7 dpi, the pattern is different, with infection with the P115 strain resulting in a predominant increase of proteins abundance, and with P5 exhibiting a majority of decrease events. By 14 and 21 dpi, both viruses produce a broadly similar response, with the vast majority of changes in protein abundance involving increase. For the kidney total proteome, at 4 dpi, P5 infection induced an approximately equivalent number of proteins with increased and decreased abundance, whereas P115 infection induced mainly increase of proteins abundance. By 7, 14, and 21 dpi, both P5 and P115 infection induced mainly increase of host proteins abundance.

Identification of differentially changed proteins in abundance by MS

All interesting changed protein spots in trachea and kidney tissues were excised, and identified by MALDI-TOF-TOF MS. From trachea and kidney tissues, 24 and 34 proteins were successfully identified, respectively. Detailed information on the identified proteins is provided in Tables 1, 2, 3 and 4, Additional file 3, and Additional file 4. Figures 2 and 3 show representative gels indicating protein spots identified in tracheal and kidney tissues respectively. As shown in Tables and Figures, several proteins were identified in multiple spots with the same molecular weight but different isoelectric points, such as three spots contained lamin-A (spot # 4317, 4310, 4321 in Figure 2).

Table 1 Similar abundance changed proteins in tracheal from chicken infected with IBV ck/CH/LDL/97I P5 and ck/CH/LDL/97I P115
Table 2 Differential abundance changed proteins in tracheal from chicken infected with IBV ck/CH/LDL/97I P5 and ck/CH/LDL/97I P115
Table 3 Similar abundance changed proteins in kidney from chicken infected with IBV ck/CH/LDL/97I P5 and ck/CH/LDL/97I P115
Table 4 Differential abundance changed proteins in kidney from chicken infected with IBV ck/CH/LDL/97I P5 and ck/CH/LDL/97I P115

According to annotations from the UniProt Knowledgebase (UniProtKB) and the Gene Ontology databases, most of the identified proteins were involved in cytoskeleton organization, anti-oxidative stress, the stress response, acute phase response, energy metabolism, macromolecular biosynthesis, signal transduction and ion transport (summarized in Additional file 5). Among them, the abundance of annexin A2, annexin A5, pyruvate kinase (PKM2), alpha-enolase, mitochondrial phosphoenolpyruvate carboxykinase (PCK2), triosephosphate isomerase, heat shock protein beta-1, manganese-containing superoxide dismutase (MnSOD), vimentin, lamin-A, cartilage matrix protein, alpha-tropomyosin, nucleoside diphosphate kinase, sulfotransferase, and low molecular weight phosphotyrosine proteinphosphatase, were induced to be differentially patterns changed in chickens infected with the P5 and P115 strains, suggested that infection with the P5 and P115 strains produces different host response. Furthermore, this difference in the pattern of change was induced predominantly in the early stages of the infection cycle, which suggests that critical events early in infection are likely to be of key importance in determining the fate of the host. In addition, fibrinogen β and γ chains showed a similar pattern of change, with increase in the P5- and P115-infected groups.

Gene ontology annotations of differentially changed proteins in abundance

In order to generate an overview of the subcellular location and biological processes of this abundance changed proteins in trachea and kidney, P5-infected and P115-infected group, and different dpi, categorization of these proteins was performed on the basis of Gene Ontology (GO) annotations. As shown in Figure 4A, cellular component ontology revealed that the majority of the identified proteins were associated with intracellular (GO:0005622) and cytoplasm (GO:0005737) both in trachea and kidney, P5-infected and P115-infected group, and different dpi. As shown in Figure 4B, biological process ontology revealed that the majority of the identified proteins were associated with metabolic processes (GO:0008152), response to stimulus (GO:0050896), and regulation of biological process (GO:0050789) both in trachea and kidney, P5-infected and P115-infected group, and different dpi.

Figure 4
figure 4

Gene Ontology annotation analysis of differentially expressed proteins according to their cellular component (A) and biological process (B). This classification was produced on the basis of an analysis using the GOSlimViewer tool at the Agbase database http://www.agbase.msstate.edu/. The results were compared between trachea and kidney, P5-infected and P115-infected group, and different day's post-inoculation.

Analysis of mRNA level by real-time RT-PCR

The alterations in the mRNA level of 12 selected proteins in the trachea or kidney from the control, P5-infected, and P115-infected groups were analyzed at 4, 7, 14, and 21 dpi. For genes from tracheal tissue, as shown in Figure 5A, at 4 dpi, the mRNA level was up-regulated in the P5-infected group, when compared with the P115-infected group. By 7 and 14 dpi, the majority of the mRNA exhibited similar down-regulation in both the P5-infected group and the P115-infected group. By 21 dpi, their change patterns were different: the mRNA level of these genes was up-regulated in the P5-infected group, but for the P115-infected group, they remained low at this time point. For genes from kidney tissue, as shown in Figure 5B, at 4 dpi, infection with the P115 strain resulted in a stronger up-regulation in the mRNA level than infection with the P5 strain. However, by 7 dpi, infection with P5 induced a stronger down-regulation than P115 infection. At 14 dpi, both virus strains resulted in a similar response, with the majority of mRNAs showing a down-regulation. In contrast, at 21 dpi, the majority of mRNAs were still down-regulated in the P115-infected group, but P5 infection induced up-regulation of the mRNA level of several genes.

Figure 5
figure 5

Transcript analysis of the 12 proteins differentially expressed in trachea (A) and kidney (B) of chickens infected with the IBV ck/CH/LDL/97I P 5 and ck/CH/LDL/97I P 115 strains by real-time RT-PCR. Total RNA extracts were prepared from the trachea and kidney of chickens in the control, P5-infected, and P115-infected group's at all four time points. Data represent means of three biological replicates per group. Error bars indicate standard error. Samples were normalized with the expression of the 18S ribosomal RNA gene. For symbols indicating different genes, refer to Tables 1, 2, 3, and 4.

In comparison with the results obtained using proteomics (summarized in Additional file 6 and Additional file 7), the trends in the mRNA levels of these genes were not completely consistent with the change patterns of their corresponding proteins in 2-DE or 2-DIGE gels. The disparity between the levels of mRNA and their corresponding proteins may occur because posttranscriptional and posttranslational modifications, as well as differential mRNA and protein degradation rates, may also contribute to these discrepancies. Several papers have addressed this question performing parallel proteomic/gene expression studies [23–26].

Protein validation by western blot analysis

To confirm the dynamic alterations of protein abundance during infection with IBV ck/CH/LDL/97I P5 and ck/CH/LDL/97I P115, three proteins including annexin A2 (ANXA2), annexin A5 (ANXA5), and heat shock protein beta-1 (HSPB1) were selected for western blot analysis. As shown in Figure 6A and 6B, the abundance of HSPB1 was increased in trachea tissues from the P5-infected group at 4, 7, and 21 dpi, relative to the control group. It was also increased in the P115-infected group, relative to the control group. The abundance of Annexin A2 was increased in tracheal tissues from the P115-infected group at 4 and 7 dpi, but decreased in tracheal tissues from the P5-infected group at 7 dpi, relative to the control group. In addition, the abundance of annexin A5 was increased in P5-infected kidney tissues at 7 dpi, relative to the control group. These data were in agreement with the results obtained from the 2-DE and 2-DIGE analysis.

Figure 6
figure 6

(A) Western blot analysis of heat shock protein beta-1 (HSPB1) and annexin A2 (ANXA2) in trachea tissue, annexin A5 (ANXA5) in kidney tissue of chickens infected with the IBV ck/CH/LDL/97I P 5 and ck/CH/LDL/97I P 115 strain. (B) Protein band density was analyzed with the BandScan software. Beta-Actin was used as the internal control. Mean values ± SE were calculated from three independent samples. AU, arbitrary units.

Discussion

Viruses are obligatory intracellular pathogens that rely on the host cell for essentially all steps of their life cycle. Although coronaviruses use host proteins as part of their replication strategies, it has also become clear that the immune, metabolic, stress, cell cycling and other pathways are activated by infection. Determination, using genomics and proteomics, of the extent to which virus-host interaction is coronavirus-specific and organ-specific, will be of importance [14, 27]. Our previous study on the global profile of host protein alterations in response to IBV infection was focused on an in ovo infection model system [21]. In this study, we applied a comparative proteomics technical platform for the first time to explore the abundance changed protein in trachea and kidney tissues from chicken infected in vivo with the highly virulent ck/CH/LDL/97I P5 strain of IBV and the embryo-passaged, attenuated ck/CH/LDL/97I P115 strain of IBV. Some identified proteins are likely to be important in the host response to virus infection, including cytoskeletal proteins, stress response proteins, and anti-oxidative proteins. Interestingly, the abundance of these proteins showed different change patterns with IBV strains differ in virulence, suggesting that some of these differences might be responsible for virulence, and consistent with our previous study which demonstrated the differences in pathology and virulence for these two different viruses [22, 28]. These results provide an overview of the proteome profile of the host in response to different virulent IBV infection in vivo.

The cytoskeletal network is a cellular scaffold system whose functions include maintenance of cellular shape, enabling of cellular migration, cell division, intracellular transport, signaling, and membrane organization. Some host cytoskeletal proteins have been reported as differentially altered by virus infection in quantitative proteomic studies [16–18]. Recently, study has revealed that several different cellular proteins involved in cell morphology and the cytoskeleton changed in abundance in IBV infected cells [19, 20]. IBV infection resulted in a number of changes to the nucleolus both in terms of gross morphology and protein content [29]. Our results presented a tendency that the abundance of several cytoskeletal proteins showed increased with the degree of virulence getting strong, which is possibly due to the collapse and dispersal of the cytoskeleton in the IBV-infected cells, as demonstrated by other viruses [30]. Vimentin is a major component of type III intermediate filaments that has been reported to be redistributed in cells around sites of virus replication and assembly during virus infection. For instance, infection with African swine fever virus (ASFV) can lead to the rearrangement of vimentin into a cage surrounding a virus factory, which may prevent movement of viral components into the cytoplasm and concentrate late structural proteins at sites of virus assembly [31]. Specifically, vimentin was observed to be increased in abundance in the cytoplasmic proteome of IBV-infected cells [19, 20]. In our study, the abundance of vimentin was increased more strongly in the highly virulent P5 strain infected group, compared with the attenuated P115 strain infected group at 21 dpi. Given that vimentin can confer rigidity to domains of the cytoplasm, the vimentin cage may provide a physical scaffold to facilitate the construction of the virus factory. The stronger induction of vimentin could be a result of the large replication capacity of the highly virulent IBV strain. Tropomyosin belongs to the family of actin-binding proteins that serve important functions in microfilament stabilization, regulation of microfilament branching, actin polymerization, and intracellular transport. The abundance of alpha-tropomyosin was found increased in IBV-infected cell [19, 20]. In our study, the abundance of alpha-tropomyosin was increased much more strongly with highly virulent IBV ck/CH/LDL/97I P5 infection than with the attenuated P115 strain, it is likely that it help to orchestrate virus assembly, release and efficient cell-to-cell spread, also may due to the different virulence between IBV ck/CH/LDL/97I P5 and P115 strain. In addition, in our previous in ovo infection model [21], the abundance of some cytoskeletal including alpha-tropomyosin was decreased in the IBV-infected embryonic tissues, this difference need to be further investigated.

In both our previous [21] and current study, the abundance of annexin A2 and annexin A5 were found to be changed upon IBV infection. Annexins are a family of conserved proteins characterized by their ability to bind and order charged phospholipids in membranes, often in response to elevated intracellular calcium. These family members are involved specifically in a diverse range of cellular functions both inside the cell and extracellularly [32]. Annexin A2 (ANXA2) can associate with actin filaments and mediates membrane trafficking and membrane-cytoskeletal interactions. It has been identified as an important host factor for several viruses and at different stages of their life cycle. ANXA2 is incorporated into cytomegalovirus and influenza virus particles, promotes the entry of virus, and plays a role in Human immunodeficiency virus -1 (HIV-1) assembly, Bluetongue virus (BTV) release, and hepatitis C virus (HCV) replication [33–40]. Annexin A5 (ANXA5) is involved in various intra- and extracellular processes including signal transduction, anti-inflammatory processes, membrane trafficking, and ion channel activity [41], it also acts to regulate blood coagulation, binding to and shielding exposed phospholipids and masking their pro-thrombotic properties [42]. The expression of ANXA2 and ANXA5 has been found to be altered during other IBV infection system [19, 20]. In this study, the abundance of ANXA2 showed weaker increase at early stages of infection with the highly virulent IBV ck/CH/LDL/97I P5 strain than with the attenuated P115 strain. For ANXA5, the level of expression induced at 7 dpi with the highly virulent P5 strain of IBV was greater than that induced by the attenuated P115 strain. These trends were also validated by western blot analysis. Whether this difference between P5 and P115 infectioned groups is related with the virulence of IBV or various responses of the host cell deserves further study.

MnSOD is the primary antioxidant enzyme in the mitochondria that catalyzes the conversion of superoxide molecules to hydrogen peroxide and molecular oxygen and therefore forms one of the cell's major defense mechanisms against oxidative stress [43]. Studies have reported that viral infections cause oxidative stress, which is associated with the activation of phagocytes and an increase in the release of reactive oxygen species (ROS) that play a positive modulatory role in immune activation, the inflammatory response, eradication of viral infections and immunity-induced cellular injury [44]. The expression of MnSOD was altered during infection with some viruses, such as porcine reproductive and respiratory syndrome virus, influenza A virus, and SARS-CoV [45–47]. The abundance of the anti-oxidative protein MnSOD was increased much significantly in highly virulent IBV ck/CH/LDL/97I P5 infected group than the attenuated P115 strain infected group. The much amplitude increased production of MnSOD during highly virulent strain infection can affect the host cell pro-/anti-oxidant balance, which probably results in more significant immune activation and a stress response induced by reactive oxygen species. It is important in the eradication of viral infections and immune-induced cellular injury.

When comparing the protein profiles of groups infected with different strains, we found that the abundance of HSPB1, a well-known heat-shock protein, in highly virulent IBV ck/CH/LDL/97I P5 strain infected group was significantly higher than the embryo-passaged, attenuated P115 strain at the early stage of infection. Western blot and real-time RT-PCR further confirmed this alteration. HSPB1 has different cytoprotective roles, including acting as a molecular chaperone, maintaining the normal function of cells through interaction with and stabilization of the cytoskeleton, regulation of translational initiation, modulation of inflammation, inhibition of apoptosis, stimulation of innate and adaptive immune responses, and responding to a wide variety of stressful stimuli [48–50]. It is probably involved in all phases of the viral life cycle, including cell entry, virion disassembly, viral genome transcription, replication and morphogenesis. The abundance of HSPB1 was observed to be increased in IBV-infected Vero cells using stable isotope labeling with amino acids in cell culture (SILAC) [20]. IBV infection induces cell cycle arrest at both S and G2/M phases and caspase-dependent apoptosis at late stages of the viral infection cycle [51–53]. Results presented in our study were probably due to the highly virulent virus triggering a markedly more robust inflammation and stress response. It also suggests that many important, and probably different events in IBV pathogenesis and immunology, such as the stress response, inflammation, and apoptosis, occur early in infection, and that these events may contribute to development of an appropriate immune response and the outcome of viral infection. In our previous study [21], the abundance of HSPB1 was found to be decreased after IBV infection in ovo. This difference in different infection model need to be further examined.

In this study, infection with both the highly virulent strain and the embryo-passaged, attenuated strain of IBV was accompanied by elevation of proteins related to energy metabolism. Viruses are obligate parasites that are completely dependent on their host's cellular metabolism for reproduction. Viral infection has been shown to modulate the levels of numerous host metabolic components in pathways such as glycolysis, the tricarboxylic acid (TCA) cycle, pentose phosphate pathway, and macromolecular biosynthesis in order to redirect valuable resources to their own mass production. Studies have demonstrated that the rate of glycolysis in cells infected by Rous sarcoma virus, feline leukemia virus, and poliomyelitis virus was increased by as much as 370% [54]. Glycolytic enzymes are involved not only in carbohydrate metabolism and increased ATP production, also take part in the induction of anti-oxidative stress in host cells and contribute to transcription of RNA virus genomes [55]. The abundance of some proteins involved in glycolysis, such as aldose reductase, pyruvate kinase, alpha-enolase, and triosephosphate isomerase was found increased in IBV infected cell [19, 20]. Our study also showed that the abundance of some glycolytic enzymes, including phosphoenolpyruvate carboxykinase (PCK1), mitochondrial phosphoenolpyruvate carboxykinase (PCK2), pyruvate kinase, and alpha-enolase, were increased significantly in the early stage of infection with both the highly virulent and the attenuated IBV strain. These findings support the view that IBV infection probably triggers activation of host energy metabolic components via systemic or global mechanisms, to keep up with the energy demands of its own replication.

The post-translational modification (PTM) plays critical roles in cellular regulation. It has been reported that PTM and expression of highly related gene sequences can induce multiple protein spots in animals [56]. An advantage of 2-DE is the monitoring of multiple forms of a protein species, which offers the opportunity to investigate the effects of the virus infection on protein modification. Our study revealed that several proteins were expressed by multiple spots in our study. Although we do not know at present how these multiple spots were generated, some of them were more likely due to post-translational modifications (PTMs), or highly related gene sequences coding isoforms [57].

Conclusion

In summary, we investigated the proteome profiles of tracheal and kidney tissues from chicken infected with highly virulent IBV ck/CH/LDL/97I P5 and embryo-passaged, attenuated IBV ck/CH/LDL/97I P115. Some proteins involved in cytoskeleton organization, stress response, and anti-oxidative stress, showed different change in abundance with IBV strains differing in virulence. While the roles of some identified alterations could be related to host antiviral response or pathogenic mechanisms, functional significance of other altered proteins remains unclear and needs further investigation. However, our findings provide proteome-related information on a large scale that should be useful in increasing our understanding of the pathogenic and immune mechanisms of IBV infection. In addition, they will also provide reference for similar research into other coronaviruses.

Materials and methods

Experimental animals and infecting virus strains

One-day-old White Leghorn specific-pathogen-free (SPF) chickens were obtained from the Laboratory Animal Center, Harbin Veterinary Research Institute, the Chinese Academy of Agricultural Sciences, China. The chickens were maintained in isolators with negative pressure, and food and water were provided ad libitum. This study was approved by the Animal Welfare Committee of Heilongjiang Province, China.

The highly virulent ck/CH/LDL/97I P5 strain of IBV and attenuated IBV strain ck/CH/LDL/97I P115 were obtained by passaging the IBV ck/CH/LDL/97I strain, which was isolated and identified by our laboratory as described previously [22]. The IBV strains were propagated once in 9- to 11-day-old embryonated SPF chicken eggs for preparation of seed stock. The presence of viral particles in the allantoic fluid of inoculated eggs was confirmed with a negative contrast electronic microscope (JEM-1200, EX) as described previously [22].

Experimental design

Sixty-six one-day-old SPF White Leghorn chickens were divided randomly into three groups: the P5-infected group, P115-infected group and control group. Twenty-two chickens were allocated to each group and housed in different isolators. Chickens in the P5-infected and P115-infected groups were inoculated with the IBV ck/CH/LDL/97I P5 strain and ck/CH/LDL/97I P115 strain, respectively, by oculo-nasal application at 11 days of age with a dose of log106 EID50 per chicken. The chickens in the control group were mock-inoculated in parallel with sterile allantoic fluid.

Blood samples from 10 chickens in each group were collected at 4, 7, 14, and 21 days post inoculation. The sera were collected for enzyme-linked immunosorbent assay (ELISA) testing. Three chickens were selected randomly from each group and killed humanely at 4, 7, 14, and 21 days post inoculation, respectively. Trachea and kidney tissues were separated rapidly and washed with ice-cold PBS buffer, snap-frozen in liquid nitrogen, and kept subsequently at -80°C for 2-DE or 2-DIGE, real-time RT-PCR, and western blot analysis.

Serum antibody detection

Serum samples were assayed using a commercial IBV antibody test kit (IDEXX Corporation, Westbrook, Maine, USA) according to the manufacturer's instructions. Each sample was tested in triplicate. Serum-to-positive ratios (S/P ratios) were calculated as described previously [2, 58]. Individual serum titers were calculated from these S/P ratios, evaluated as positive or negative, and expressed as OD650 nm values according to the manufacturer's instructions.

Quantitative analysis of IBV in trachea and kidney by real-time RT-PCR

Viral load of IBV was analyzed using R. M. Jones described methods [59]. Tissue samples were ground and homogenized, 100 mg tissue homogenates were suspended in 500 μl phosphate-buffered saline (PBS) containing 100 μg penicillin and 100 μg streptomycin/ml, the suspension were freezed thawing three times, then centrifuged at 13,000 × g at 4°C for 5 minutes. RNA was extracted from 200 μl tissue supernatant using TRIzol Reagent (Invitrogen) following the manufacturer's instructions. The real-time RT-PCR assay used the following primer and probe sequences: IBVF forward primer CTA TCG CCA GGG AAA TGT C, IBVR reverse primer GCG TCC TAG TGC TGT ACC C, IBV TaqMan® probe FAM-CCTGGAAACGAACGGTAGACCCT-TAMRA [59]. One-step real-time RT-PCR reactions were performed using One Step PrimeScript® RT-PCR kit (TaKaRa Biotech Co. Ltd., Dalian) on LightCycler® 480 real-time PCR system (Roche) according to the following steps: reverse transcription at 42°C for 10 min, denaturation at 95°C for 10 s. and 40 cycles with 95°C for 5 s, 55°C for 20 s, 72°C for 10 s, followed by a 40°C for 10 s cooling step. All of the samples were tested in triplicate in each reaction. The data were analyzed using LightCycler® 480 Software Version 1.5.

Preparation of protein samples

The frozen trachea or kidney tissues were placed in liquid nitrogen and ground thoroughly to a very fine powder. Samples of 100 mg tissue powder were dissolved in 500 μl lysing solution containing 7 M urea, 2 M thiourea, 4% CHAPS, 40 mM DTT, 2% IPG buffer pH 3-10, 1% Nuclease Mix and 1% Protease Inhibitor Mix (GE Healthcare), then incubated for 2 h at 4°C with vortexing once every 15 min, and centrifuged at 15,000 × g for 1 h at 4°C. The supernatant was collected and purified with the PlusOne 2D Clean-up kit (GE Healthcare). The concentration of each protein sample was determined with the PlusOne 2D Quant Kit (GE Healthcare). Protein samples were aliquoted and stored at -80°C for 2-DE or 2-DIGE analysis.

2-DE analysis of trachea protein samples

Thirty-six samples of tracheal protein from the three groups (P5-infected, P115-infected and control) at 4, 7, 14, and 21 days post-inoculation (dpi) were analyzed by 2-DE using the methods previously described [21]. First, 400 μg protein samples were added to rehydration solution (7 M urea, 2 M thiourea, 40 mM DTT, 2% CHAPS, 0.5% pH 3-10 or pH 4-7 IPG buffer, and 0.002% bromophenol blue) to make the final volume up to 250 μl, following Isoelectric focusing (IEF), the IPG strips were equilibrated and the second dimension separation was performed on 12.5% SDS-polyacrylamide gels on the SE600 Ruby system (GE Healthcare). Then, the gels were stained with PlusOne Coomassie Blue R-350 (GE Healthcare) and scanned with an ImageScanner III (GE Healthcare). Quantification analyses were performed with Image Master 2D Platinum software v6.0 (GE Healthcare). For image analysis, three independent gels from the P5-infected and P115-infected groups were compared with those from the corresponding control group at 4, 7, 14, and 21 dpi respectively. The normalized volume values (vol %) of matched protein spots were subjected to Student's t test using the SPSS statistical software package version 16.0. The criterion used to define differential expression of spots was that the ratio of the vol % in the P5-infected group or the P115-infected group vs. the control group was more than 1.5 (p < 0.05) or less than 0.67 (p < 0.05). Differentially changed protein spots were excised manually from the gels and subjected to MS analysis.

2-DIGE analysis of kidney protein samples

The kidney protein samples from each group at 4, 7, 14, and 21 dpi were used for 2-DIGE analysis respectively. The pH of the protein samples was adjusted to 8.5 and the protein concentration was adjusted to 5 μg/μl. Equal amounts of protein from each sample were pooled together as the internal standard. The proteins were minimally labeled according to the manufacturer's instructions (CyDye DIGE fluor minimal labeling kit, GE Healthcare). The Cy2 was used to label the pooled internal standard, and Cy3 and Cy5 were used randomly to label samples from the control group, P5-infected group and P115-infected group. To minimize system error and inherent biological variation, the sample multiplexing was also randomized (see Additional file 8) to produce unbiased results. Following the labeling reaction, 50 μg of each Cy2, Cy3, and Cy5 labeled sample was mixed, the pooled sample of each gel was adjusted to 450 μl with rehydration buffer (7 M urea, 2 M thiourea, 2%CHAPS, 2%v/v IPG buffer,130 mM DTT), and loaded subsequently onto 24 cm, linear pH 3-10 IPG strips (GE Healthcare). Isoelectric focusing (IEF) was performed on an Ettan™ IPGphor 3 (GE Healthcare) using the following protocol: 30 V for 12 h, 200 V for 2 h, 500 V for 2 h, 1000 V for 2 h, 8000 V for 3 h, and 8000 V for 65000 Vh. Then, the IPG strips were equilibrated and the second dimension separation was conducted on 12.5% SDS-polyacrylamide gels using an Ettan™ DALT six system (GE Healthcare). The CyDye-labeled gels were scanned using a Typhoon 9400 (GE Healthcare). Image analysis was performed using Ettan™ DeCyder Software version v6.5 (GE Healthcare) as described in the user's manual. The statistical analysis of changes in protein abundance in different gels was performed on the basis of the spot volumes. Protein spots with statistically significant results for the Student's t-test (p < 0.05) and Average Ratio more than1.5 or less than 0.67 were considered to be differentially changed protein spots.

The preparative gels were made with 1200 μg of unlabeled proteins (400 μg from each group). These mixed proteins were loaded and separated under the same conditions as described above. The preparative gels were stained with Coomassie Blue R-350. Each spot of interest, defined on the basis of CyDye images, was matched with a Coomassie Blue R-350 image, then excised manually from the gel and subjected to MS analysis.

MALDI-TOF/TOF MS and database search

The gel samples were placed in a tube and washed twice with 500 μl and 250 μl ddH2O for 15 min. For trypsin digestion, the gel samples were washed twice with 50 mM w/v NH4HCO3 and covered with 10 mg/ml Porcine Trypsin solution (Promega, Madison, WI, USA) in 50 mM w/v NH4HCO3. After incubation overnight at 37°C, the supernatant was removed into a second tube and 40 μl 50 mM w/v NH4HCO3 was added. Gel samples were washed with 40 μl of 50 mM w/v NH4HCO3, the supernatant was collected, and both collected supernatants were combined. The gel was washed with 70% v/v ACN and dried in a Speed Vac (Vacuum Concentrator, Bachhofer). The peptide mixtures were desalted using ZipTip C-18 RP tips (Millipore, Billerica, MA, USA) which were wetted with 100% ACN and equilibrated with 0.1% TFA. Peptide samples, which were redissolved in 10 ml 0.5% TFA, were eluted with 50% ACN/0.1% TFA and dried in a Speed Vac (Vacuum Concentrator).

The purified peptides were spotted on a MALDI plate and covered with 0.7 μl of 2 mg/ml 3, 5-Dimethoxy-4-hydroxycinnamic acid matrix (Sigma) with 10 mM NH4H2PO4 in 60% ACN. All samples were analyzed by MALDI-TOF/TOF MS with a 4700 Proteomics Analyzer (Applied Biosystems, Foster City, CA). Monoisotopic peak masses were acquired in a mass range of 800-4000 Da, with a signal-to-noise ratio (S/N) of 200. Five of the most intense ion signals, excluding common trypsin autolysis peaks and matrix ion signals, were selected as precursors for MS/MS acquisition. The peptide mass fingerprint (PMF) combined MS/MS data were submitted to MASCOT version 3.0 (Matrix Science) for identification according to the NCBInr database (release 16/01/2010, 10343571 sequences, 3528215794 residues). The search parameters were set as follows: Gallus, trypsin cleavage (one missed cleavage allowed), carbamidomethylation of cysteine as fixed modification, oxidation of methionine as variable modification, peptide mass tolerance set at 100 ppm, fragment tolerance set at 0.8 Da. The criterion for successful identification of a protein was the protein score confidence interval (C.I. %) ≥ 95%.

Gene ontology (GO) annotation of differentially expressed proteins

Spot identities were submitted to GORetriever http://www.agbase.msstate.edu/ to obtain the GO annotations. If no annotation was returned, GOanna was used to retrieve GO annotations assigned on the basis of sequence similarities. The resulting annotations were summarized according to the GOA and whole proteome GOSlim set using GOSlimViewer [60].

Analysis of mRNA levels by real-time RT-PCR

Specific primers were designed according to the corresponding gene sequences of the MS-identified proteins using Beacon Designer software 7.5 (Primer Biosoft International). Information on the primers is listed in Additional file 9. Total RNA was extracted using TRIzol Reagent (Invitrogen) according to the manufacturer's instructions. Two micrograms of total RNA was reverse transcribed with 200 U M-MLV Reverse Transcriptase (Invitrogen) and 500 ng Oligo(dT)18 as the first strand primer in 20 μl reaction solution. Real-time RT-PCR was carried out using SYBR® Premix Ex Taq™ II kit (Takara) on the LightCycler® 480 real-time PCR system (Roche) according to the following steps: 30 s at 95°C, 40 cycles of denaturation at 95°C for 10 s, and annealing and extension at 55°C for 45 s. Each sample was amplified in triplicate. Quantitative analysis of the data was performed using the 2-ΔΔCt method [61]; samples from the control group at 4 dpi were used as the calibrator (relative expression = 1), and the 18S ribosomal RNA gene was used as an internal reference gene.

Western blot analysis

Samples of tracheal and kidney proteins from the P5-infected group, P115-infected group, and control group at all four time points were prepared, and the protein concentration was determined as described above. Equivalent amounts of total protein were subjected to 12% SDS-PAGE and then transferred to a nitrocellulose membrane. After blocking for 1 h at 37°C, the membranes were incubated separately with mouse monoclonal antibody to annexin A5 (sc-32321, Santa Cruz Biotechnology, USA), mouse monoclonal antibody to HSPB1 (sc-51956, Santa Cruz Biotechnology, USA), and goat polyclonal antibody to annexin A2 (sc-1924, Santa Cruz Biotechnology, USA) overnight at 4°C. The membranes were incubated subsequently with IRDye700DX conjugated anti-mouse secondary antibody (610-130-121, Rockland, Gilbertsville, PA) for 1 h at 37°C, and scanned finally on a LI-COR infrared imaging system using their Odyssey software (Li-Cor Bioscience, Lincoln, NE). β-actin (sc-47778, Santa Cruz Biotechnology, USA) was used as a reference protein to check equal loading. Triplicate experiments were performed for each sample. Densitometric analysis of protein bands was done by using BandScan software (Glyko).

Abbreviations

DPI:

Days post infection

ACN:

Acetonitrile

CHAPS:

3-[(3-cholamidopropyl) dimethyl-ammonio] -1-propanesulfonate

DTT:

Dithiothreitol

IBV:

Infectious bronchitis coronavirus

IEF:

Isoelectric focusing

IPG:

Immobilized pH gradient

MALDI-TOF-TOF/MS:

Matrix-assisted laser desorption/ionization time-of-flight tandem mass spectrometry

PMF:

Peptide mass fingerprinting

SDS-PAGE:

Sodium dodecyl sulfate polyacrylamide gel electrophoresis

SPF:

Specific pathogen free

TFA:

Trifluoroacetic acid

2-DE:

Two-dimensional gel electrophoresis

2-DIGE:

Two-dimensional fluorescence difference gel electrophoresis.

References

  1. Carstens E: Report from the 40th meeting of the Executive Committee of the International Committee of Taxonomy of Viruses. 2009, 1571–1574.

    Google Scholar 

  2. Han Z, Sun C, Yan B, Zhang X, Wang Y, Li C, Zhang Q, Ma Y, Shao Y, Liu Q, et al.: A 15-year analysis of molecular epidemiology of avian infectious bronchitis coronavirus in China. Infect Genet Evol 2011, 11: 190–200. 10.1016/j.meegid.2010.09.002

    Article  Google Scholar 

  3. Enjuanes L, Almazan F, Sola I, Zuniga S: Biochemical aspects of coronavirus replication and virus-host interaction. Annu Rev Microbiol 2006, 60: 211–230. 10.1146/annurev.micro.60.080805.142157

    Article  CAS  Google Scholar 

  4. Tang BS, Chan KH, Cheng VC, Woo PC, Lau SK, Lam CC, Chan TL, Wu AK, Hung IF, Leung SY, Yuen KY: Comparative host gene transcription by microarray analysis early after infection of the Huh7 cell line by severe acute respiratory syndrome coronavirus and human coronavirus 229E. J Virol 2005, 79: 6180–6193. 10.1128/JVI.79.10.6180-6193.2005

    Article  CAS  Google Scholar 

  5. Dar A, Munir S, Vishwanathan S, Manuja A, Griebel P, Tikoo S, Townsend H, Potter A, Kapur V, Babiuk LA: Transcriptional analysis of avian embryonic tissues following infection with avian infectious bronchitis virus. Virus Res 2005, 110: 41–55. 10.1016/j.virusres.2005.01.006

    Article  CAS  Google Scholar 

  6. Wang X, Rosa AJ, Oliverira HN, Rosa GJ, Guo X, Travnicek M, Girshick T: Transcriptome of local innate and adaptive immunity during early phase of infectious bronchitis viral infection. Viral Immunol 2006, 19: 768–774. 10.1089/vim.2006.19.768

    Article  CAS  Google Scholar 

  7. Ng LF, Hibberd ML, Ooi EE, Tang KF, Neo SY, Tan J, Murthy KR, Vega VB, Chia JM, Liu ET, Ren EC: A human in vitro model system for investigating genome-wide host responses to SARS coronavirus infection. BMC Infect Dis 2004, 4: 34. 10.1186/1471-2334-4-34

    Article  Google Scholar 

  8. Leong WF, Tan HC, Ooi EE, Koh DR, Chow VT: Microarray and real-time RT-PCR analyses of differential human gene expression patterns induced by severe acute respiratory syndrome (SARS) coronavirus infection of Vero cells. Microb Infect/Institut Pasteur 2005, 7: 248–259. 10.1016/j.micinf.2004.11.004

    Article  CAS  Google Scholar 

  9. Raaben M, Groot Koerkamp MJ, Rottier PJ, de Haan CA: Mouse hepatitis coronavirus replication induces host translational shutoff and mRNA decay, with concomitant formation of stress granules and processing bodies. Cellular Microbiol 2007, 9: 2218–2229. 10.1111/j.1462-5822.2007.00951.x

    Article  CAS  Google Scholar 

  10. Haq K, Brisbin JT, Thanthrige-Don N, Heidari M, Sharif S: Transcriptome and proteome profiling of host responses to Marek's disease virus in chickens. Vet Immunol Immunopathol 2010, 138: 292–302. 10.1016/j.vetimm.2010.10.007

    Article  CAS  Google Scholar 

  11. Jiang XS, Tang LY, Dai J, Zhou H, Li SJ, Xia QC, Wu JR, Zeng R: Quantitative analysis of severe acute respiratory syndrome (SARS)-associated coronavirus-infected cells using proteomic approaches: implications for cellular responses to virus infection. Mol Cell Proteomics 2005, 4: 902–913. 10.1074/mcp.M400112-MCP200

    Article  CAS  Google Scholar 

  12. Chen JH, Chang YW, Yao CW, Chiueh TS, Huang SC, Chien KY, Chen A, Chang FY, Wong CH, Chen YJ: Plasma proteome of severe acute respiratory syndrome analyzed by two-dimensional gel electrophoresis and mass spectrometry. Proc Natl Acad Sci USA 2004, 101: 17039–17044. 10.1073/pnas.0407992101

    Article  CAS  Google Scholar 

  13. Zhang L, Zhang ZP, Zhang XE, Lin FS, Ge F: Quantitative proteomics analysis reveals BAG3 as a potential target to suppress severe acute respiratory syndrome coronavirus replication. J Virol 2010, 84: 6050–6059. 10.1128/JVI.00213-10

    Article  CAS  Google Scholar 

  14. Vogels MW, van Balkom BW, Kaloyanova DV, Batenburg JJ, Heck AJ, Helms JB, Rottier PJ, de Haan CA: Identification of host factors involved in coronavirus replication by quantitative proteomics analysis. Proteomics 2011, 11: 64–80. 10.1002/pmic.201000309

    Article  CAS  Google Scholar 

  15. Emmott E, Wise H, Loucaides EM, Matthews DA, Digard P, Hiscox JA: Quantitative proteomics using SILAC coupled to LC-MS/MS reveals changes in the nucleolar proteome in influenza A virus-infected cells. J Proteome Res 2010, 9: 5335–5345. 10.1021/pr100593g

    Article  CAS  Google Scholar 

  16. Munday DC, Emmott E, Surtees R, Lardeau CH, Wu W, Duprex WP, Dove BK, Barr JN, Hiscox JA: Quantitative proteomic analysis of A549 cells infected with human respiratory syncytial virus. Mol Cell Proteomics 2010, 9: 2438–2459. 10.1074/mcp.M110.001859

    Article  CAS  Google Scholar 

  17. Coombs KM, Berard A, Xu W, Krokhin O, Meng X, Cortens JP, Kobasa D, Wilkins J, Brown EG: Quantitative proteomic analyses of influenza virus-infected cultured human lung cells. J Virol 2010, 84: 10888–10906. 10.1128/JVI.00431-10

    Article  CAS  Google Scholar 

  18. van Diepen A, Brand HK, Sama I, Lambooy LH, van den Heuvel LP, van der Well L, Huynen M, Osterhaus AD, Andeweg AC, Hermans PW: Quantitative proteome profiling of respiratory virus-infected lung epithelial cells. J Proteomics 2010, 73: 1680–1693. 10.1016/j.jprot.2010.04.008

    Article  CAS  Google Scholar 

  19. Emmott E, Smith C, Emmett SR, Dove BK, Hiscox JA: Elucidation of the avian nucleolar proteome by quantitative proteomics using SILAC and changes in cells infected with the coronavirus infectious bronchitis virus. Proteomics 2010, 10: 3558–3562. 10.1002/pmic.201000139

    Article  CAS  Google Scholar 

  20. Emmott E, Rodgers MA, Macdonald A, McCrory S, Ajuh P, Hiscox JA: Quantitative proteomics using stable isotope labeling with amino acids in cell culture reveals changes in the cytoplasmic, nuclear, and nucleolar proteomes in Vero cells infected with the coronavirus infectious bronchitis virus. Mol Cell Proteomics 2010, 9: 1920–1936. 10.1074/mcp.M900345-MCP200

    Article  CAS  Google Scholar 

  21. Cao Z, Han Z, Shao Y, Geng H, Kong X, Liu S: Proteomic analysis of chicken embryonic trachea and kidney tissues after infection in ovo by avian infectious bronchitis coronavirus. Proteome Science 2011, 9: 11. 10.1186/1477-5956-9-11

    Article  CAS  Google Scholar 

  22. Liu S, Zhang X, Wang Y, Li C, Liu Q, Han Z, Zhang Q, Kong X, Tong G: Evaluation of the protection conferred by commercial vaccines and attenuated heterologous isolates in China against the CK/CH/LDL/97I strain of infectious bronchitis coronavirus. Vet J 2009, 179: 130–136. 10.1016/j.tvjl.2007.08.018

    Article  CAS  Google Scholar 

  23. Conrads KA, Yi M, Simpson KA, Lucas DA, Camalier CE, Yu LR, Veenstra TD, Stephens RM, Conrads TP, Beck GR Jr: A combined proteome and microarray investigation of inorganic phosphate-induced pre-osteoblast cells. Mol Cell Proteomics 2005, 4: 1284–1296. 10.1074/mcp.M500082-MCP200

    Article  CAS  Google Scholar 

  24. Kuo CC, Kuo CW, Liang CM, Liang SM: A transcriptomic and proteomic analysis of the effect of CpG-ODN on human THP-1 monocytic leukemia cells. Proteomics 2005, 5: 894–906. 10.1002/pmic.200401144

    Article  CAS  Google Scholar 

  25. Scheurer SB, Rybak JN, Rosli C, Neri D, Elia G: Modulation of gene expression by hypoxia in human umbilical cord vein endothelial cells: A transcriptomic and proteomic study. Proteomics 2004, 4: 1737–1760. 10.1002/pmic.200300689

    Article  CAS  Google Scholar 

  26. Martin SA, Mohanty BP, Cash P, Houlihan DF, Secombes CJ: Proteome analysis of the Atlantic salmon (Salmo salar) cell line SHK-1 following recombinant IFN-gamma stimulation. Proteomics 2007, 7: 2275–2286. 10.1002/pmic.200700020

    Article  CAS  Google Scholar 

  27. Perlman S, Netland J: Coronaviruses post-SARS: update on replication and pathogenesis. Nat Rev Microbiol 2009, 7: 439–450. 10.1038/nrmicro2147

    Article  CAS  Google Scholar 

  28. Liu S, Han Z, Chen J, Liu X, Shao Y, Kong X, Tong G, Rong J: S1 gene sequence heterogeneity of a pathogenic infectious bronchitis virus strain and its embryo-passaged, attenuated derivatives. Avian Pathol 2007, 36: 231–234. 10.1080/03079450701338730

    Article  CAS  Google Scholar 

  29. Dove BK, You JH, Reed ML, Emmett SR, Brooks G, Hiscox JA: Changes in nucleolar morphology and proteins during infection with the coronavirus infectious bronchitis virus. Cell Microbiol 2006, 8: 1147–1157. 10.1111/j.1462-5822.2006.00698.x

    Article  CAS  Google Scholar 

  30. Radtke K, Dohner K, Sodeik B: Viral interactions with the cytoskeleton: a hitchhiker's guide to the cell. Cell Microbiol 2006, 8: 387–400. 10.1111/j.1462-5822.2005.00679.x

    Article  CAS  Google Scholar 

  31. Stefanovic S, Windsor M, Nagata KI, Inagaki M, Wileman T: Vimentin rearrangement during African swine fever virus infection involves retrograde transport along microtubules and phosphorylation of vimentin by calcium calmodulin kinase II. J Virol 2005, 79: 11766–11775. 10.1128/JVI.79.18.11766-11775.2005

    Article  CAS  Google Scholar 

  32. Hayes MJ, Longbottom RE, Evans MA, Moss SE: Annexinopathies. Subcell Biochem 2007, 45: 1–28. 10.1007/978-1-4020-6191-2_1

    Article  CAS  Google Scholar 

  33. Wright JF, Kurosky A, Pryzdial EL, Wasi S: Host cellular annexin II is associated with cytomegalovirus particles isolated from cultured human fibroblasts. J Virol 1995, 69: 4784–4791.

    CAS  Google Scholar 

  34. LeBouder F, Morello E, Rimmelzwaan GF, Bosse F, Pechoux C, Delmas B, Riteau B: Annexin II incorporated into influenza virus particles supports virus replication by converting plasminogen into plasmin. J Virol 2008, 82: 6820–6828. 10.1128/JVI.00246-08

    Article  CAS  Google Scholar 

  35. Gonzalez-Reyes S, Garcia-Manso A, del Barrio G, Dalton KP, Gonzalez-Molleda L, Arrojo-Fernandez J, Nicieza I, Parra F: Role of annexin A2 in cellular entry of rabbit vesivirus. J Gen Virol 2009, 90: 2724–2730. 10.1099/vir.0.013276-0

    Article  CAS  Google Scholar 

  36. Raynor CM, Wright JF, Waisman DM, Pryzdial EL: Annexin II enhances cytomegalovirus binding and fusion to phospholipid membranes. Biochemistry 1999, 38: 5089–5095. 10.1021/bi982095b

    Article  CAS  Google Scholar 

  37. Harrist AV, Ryzhova EV, Harvey T, Gonzalez-Scarano F: Anx2 interacts with HIV-1 Gag at phosphatidylinositol (4,5) bisphosphate-containing lipid rafts and increases viral production in 293T cells. PLoS One 2009, 4: e5020. 10.1371/journal.pone.0005020

    Article  Google Scholar 

  38. Ryzhova EV, Vos RM, Albright AV, Harrist AV, Harvey T, Gonzalez-Scarano F: Annexin 2: a novel human immunodeficiency virus type 1 Gag binding protein involved in replication in monocyte-derived macrophages. J Virol 2006, 80: 2694–2704. 10.1128/JVI.80.6.2694-2704.2006

    Article  CAS  Google Scholar 

  39. Beaton AR, Rodriguez J, Reddy YK, Roy P: The membrane trafficking protein calpactin forms a complex with bluetongue virus protein NS3 and mediates virus release. Proc Natl Acad Sci USA 2002, 99: 13154–13159. 10.1073/pnas.192432299

    Article  CAS  Google Scholar 

  40. Backes P, Quinkert D, Reiss S, Binder M, Zayas M, Rescher U, Gerke V, Bartenschlager R, Lohmann V: Role of annexin A2 in the production of infectious hepatitis C virus particles. J Virol 2010, 84: 5775–5789. 10.1128/JVI.02343-09

    Article  CAS  Google Scholar 

  41. Ravassa S, Bennaghmouch A, Kenis H, Lindhout T, Hackeng T, Narula J, Hofstra L, Reutelingsperger C: Annexin A5 down-regulates surface expression of tissue factor: a novel mechanism of regulating the membrane receptor repertoir. J Biol Chem 2005, 280: 6028–6035.

    Article  CAS  Google Scholar 

  42. Rand JH, Wu XX, Giesen P: A possible solution to the paradox of the "lupus anticoagulant": antiphospholipid antibodies accelerate thrombin generation by inhibiting annexin-V. Thromb Haemost 1999, 82: 1376–1377.

    CAS  Google Scholar 

  43. Yesilkaya H, Kadioglu A, Gingles N, Alexander JE, Mitchell TJ, Andrew PW: Role of manganese-containing superoxide dismutase in oxidative stress and virulence of Streptococcus pneumoniae. Infect Immun 2000, 68: 2819–2826. 10.1128/IAI.68.5.2819-2826.2000

    Article  CAS  Google Scholar 

  44. Suresh DR, Annam V, Pratibha K, Prasad BV: Total antioxidant capacity-a novel early bio-chemical marker of oxidative stress in HIV infected individuals. J Biomed Sci 2009, 16: 61. 10.1186/1423-0127-16-61

    Article  CAS  Google Scholar 

  45. Zhang H, Guo X, Ge X, Chen Y, Sun Q, Yang H: Changes in the cellular proteins of pulmonary alveolar macrophage infected with porcine reproductive and respiratory syndrome virus by proteomics analysis. J Proteome Res 2009, 8: 3091–3097. 10.1021/pr900002f

    Article  CAS  Google Scholar 

  46. Lai CC, Jou MJ, Huang SY, Li SW, Wan L, Tsai FJ, Lin CW: Proteomic analysis of up-regulated proteins in human promonocyte cells expressing severe acute respiratory syndrome coronavirus 3C-like protease. Proteomics 2007, 7: 1446–1460. 10.1002/pmic.200600459

    Article  CAS  Google Scholar 

  47. Vester D, Rapp E, Kluge S, Genzel Y, Reichl U: Virus-host cell interactions in vaccine production cell lines infected with different human influenza A virus variants: a proteomic approach. J Proteomics 2010, 73: 1656–1669. 10.1016/j.jprot.2010.04.006

    Article  CAS  Google Scholar 

  48. Samali A, Robertson JD, Peterson E, Manero F, van Zeijl L, Paul C, Cotgreave IA, Arrigo AP, Orrenius S: Hsp27 protects mitochondria of thermotolerant cells against apoptotic stimuli. Cell Stress Chaperones 2001, 6: 49–58. 10.1379/1466-1268(2001)006<0049:HPMOTC>2.0.CO;2

    Article  CAS  Google Scholar 

  49. Concannon CG, Gorman AM, Samali A: On the role of Hsp27 in regulating apoptosis. Apoptosis 2003, 8: 61–70. 10.1023/A:1021601103096

    Article  CAS  Google Scholar 

  50. Ferns G, Shams S, Shafi S: Heat shock protein 27: its potential role in vascular disease. Int J Exp Pathol 2006, 87: 253–274. 10.1111/j.1365-2613.2006.00484.x

    Article  CAS  Google Scholar 

  51. Liu C, Xu HY, Liu DX: Induction of Caspase-Dependent Apoptosis in Cultured Cells by the Avian Coronavirus Infectious Bronchitis Virus. J Virol 2001, 75: 6402–6409. 10.1128/JVI.75.14.6402-6409.2001

    Article  CAS  Google Scholar 

  52. Li F: Cell cycle arrest and apoptosis induced by the coronavirus infectious bronchitis virus in the absence of p53. Virology 2007, 365: 435–445. 10.1016/j.virol.2007.04.015

    Article  CAS  Google Scholar 

  53. Dove B, Brooks G, Bicknell K, Wurm T, Hiscox JA: Cell Cycle Perturbations Induced by Infection with the Coronavirus Infectious Bronchitis Virus and Their Effect on Virus Replication. J Virol 2006, 80: 4147–4156. 10.1128/JVI.80.8.4147-4156.2006

    Article  CAS  Google Scholar 

  54. Maynard ND, Gutschow MV, Birch EW, Covert MW: The virus as metabolic engineer. Biotechnol J 2010, 5: 686–694. 10.1002/biot.201000080

    Article  CAS  Google Scholar 

  55. Kim JW, Dang CV: Multifaceted roles of glycolytic enzymes. Trends Biochem Sci 2005, 30: 142–150. 10.1016/j.tibs.2005.01.005

    Article  CAS  Google Scholar 

  56. Tabuse Y, Nabetani T, Tsugita A: Proteomic analysis of protein expression profiles during Caenorhabditis elegans development using two-dimensional difference gel electrophoresis. Proteomics 2005, 5: 2876–2891. 10.1002/pmic.200401154

    Article  CAS  Google Scholar 

  57. Gorg A, Weiss W, Dunn MJ: Current two-dimensional electrophoresis technology for proteomics. Proteomics 2004, 4: 3665–3685. 10.1002/pmic.200401031

    Article  Google Scholar 

  58. Liu S, Zhang X, Gong L, Yan B, Li C, Han Z, Shao Y, Li H, Kong X: Altered pathogenicity, immunogenicity, tissue tropism and 3'-7 kb region sequence of an avian infectious bronchitis coronavirus strain after serial passage in embryos. Vaccine 2009, 27: 4630–4640. 10.1016/j.vaccine.2009.05.072

    Article  CAS  Google Scholar 

  59. Jones RM, Ellis RJ, Cox WJ, Errington J, Fuller C, Irvine RM, Wakeley PR: Development and Validation of RT-PCR Tests for the Detection and S1 Genotyping of Infectious Bronchitis Virus and Other Closely Related Gammacoronaviruses Within Clinical Samples. Transbound Emerg Dis 2011, 58: 411–420. 10.1111/j.1865-1682.2011.01222.x

    Article  CAS  Google Scholar 

  60. McCarthy FM, Bridges SM, Wang N, Magee GB, Williams WP, Luthe DS, Burgess SC: AgBase: a unified resource for functional analysis in agriculture. Nucleic Acids Res 2007, 35: D599-D603. 10.1093/nar/gkl936

    Article  CAS  Google Scholar 

  61. Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods 2001, 25: 402–408. 10.1006/meth.2001.1262

    Article  CAS  Google Scholar 

Download references

Acknowledgements

This work was funded by a grant from the earmarked fund for Modern Agro-industry Technology Research System.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Shengwang Liu.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

SL designed the study. SL and ZC drafted the manuscript. ZC, ZH, YS, and XL carried out virus infection, serum antibody detection, and quantitative analysis of IBV. ZC, JS and DY carried out the 2-DE experiments, image analysis, excised the protein spots, data analysis and interpretation, confirmed the differential expression by real-time RT-PCR and Western blotting analysis. SL wrote the manuscript. XK revised the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material

12953_2011_337_MOESM1_ESM.DOC

Additional file 1: Table S1 Serological results post inoculation with IBV ck/CH/LDL/97I P 5 and ck/CH/LDL/97I P 115 . (DOC 30 KB)

12953_2011_337_MOESM2_ESM.PPT

Additional file 2: Figure S1 Summary of changes in protein levels over time following infection with IBV ck/CH/LDL/97I P 5 and P 115 . The y axis shows the number of differentially expressed protein spots; individual spots can be found in Tables 1, 2, 3, and 4. (PPT 176 KB)

12953_2011_337_MOESM3_ESM.XLS

Additional file 3: This includes the PMF spectrum and confirmed MALDI-TOF-TOF spectrum of differentially expressed protein spots in IBV-infected chick tracheal tissues. (XLS 10 MB)

12953_2011_337_MOESM4_ESM.XLS

Additional file 4: This includes the PMF spectrum and confirmed MALDI-TOF-TOF spectrum of differentially expressed protein spots in IBV-infected chick kidney tissues. (XLS 13 MB)

12953_2011_337_MOESM5_ESM.DOC

Additional file 5: Table S2 Biological function of differentially expressed proteins reported in viral infection. (DOC 74 KB)

12953_2011_337_MOESM6_ESM.DOC

Additional file 6: Table S3 Comparison of the fold changes for protein abundance observed by 2-DE gel analysis and mRNA level obtained by real-time RT-PCR in trachea tissues. (DOC 52 KB)

12953_2011_337_MOESM7_ESM.DOC

Additional file 7: Table S4 Comparison of the fold changes for protein abundance observed by 2-DIGE gel analysis and mRNA expression obtained by real-time RT-PCR in kidney tissues. (DOC 48 KB)

Additional file 8: Table S5 Experiment design of the different fluorescent dye labeling. (DOC 51 KB)

Additional file 9: Table S6 The primers of Real-time RT-PCR. (DOC 42 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Cao, Z., Han, Z., Shao, Y. et al. Proteomics analysis of differentially expressed proteins in chicken trachea and kidney after infection with the highly virulent and attenuated coronavirus infectious bronchitis virus in vivo . Proteome Sci 10, 24 (2012). https://doi.org/10.1186/1477-5956-10-24

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1477-5956-10-24

Keywords