- Open Access
Quantitative phosphoproteome on the silkworm (Bombyx mori) cells infected with baculovirus
Virology Journalvolume 14, Article number: 117 (2017)
Bombyx mori has become an important model organism for many fundamental studies. Bombyx mori nucleopolyhedrovirus (BmNPV) is a significant pathogen to Bombyx mori, yet also an efficient vector for recombinant protein production. A previous study indicated that acetylation plays many vital roles in several cellular processes of Bombyx mori while global phosphorylation pattern upon BmNPV infection remains elusive.
Employing tandem mass tag (TMT) labeling and phosphorylation affinity enrichment followed by high-resolution LC-MS/MS analysis and intensive bioinformatics analysis, the quantitative phosphoproteome in Bombyx mori cells infected by BmNPV at 24 hpi with an MOI of 10 was extensively examined.
Totally, 6480 phosphorylation sites in 2112 protein groups were identified, among which 4764 sites in 1717 proteins were quantified. Among the quantified proteins, 81 up-regulated and 25 down-regulated sites were identified with significant criteria (the quantitative ratio above 1.3 was considered as up-regulation and below 0.77 was considered as down-regulation) and with significant p-value (p < 0.05). Some proteins of BmNPV were also hyperphosphorylated during infection, such as P6.9, 39 K, LEF-6, Ac58-like protein, Ac82-like protein and BRO-D.
The phosphorylated proteins were primary involved in several specific functions, out of which, we focused on the binding activity, protein synthesis, viral replication and apoptosis through kinase activity.
Bombyx mori, the domesticated silkworm, is an economically significant insect in silk production . Farmers in many developing countries such as China, India, Brazil, Vietnam, and Thailand are evolving sericulture as a foremost income source. China, the biggest cocoon producer (almost 80%) worldwide, generated a sericulture-based income of 22.4 billion ¥ (Yuan) (approximately equal to USD 3.24 billion by March 2017) by producing 6.61 × 108 kg of cocoons in 2011 . Moreover, this lepidopteran insect has become a model organism for several fundamental studies in biochemistry, molecular genetics, and genomics [3, 4]. In-depth studies of silkworm were extensively carried out since its genome was completely sequenced in 2004 . H3 N-terminus, an N-terminal tail of histone H3 in holocentric chromosomes of the silkworm, is reported to be extensively acetylated and methylated, which is almost similar to that of mammals . Furthermore, the proteome of B. mori was detected to harbor a total of 342 acetylated proteins with 667 Kac (Lysine acetylation) sites . However, some pathogenic agents such as viruses, fungi, and bacteria pose serious biological challenges for this sericulture causing almost 20% losses of the potential cocoon production each year. Bombyx mori nucleopolyhedrovirus (BmNPV), a member of baculoviridae family, is the pivotal viral pathogen of B. mori that could elicit almost 80% economic loses in sericultural industry . On the other hand, BmNPV has been used as an effective, feasible, and safe technology for the production of recombinant proteins in insect or insect cultured cells called BEVS (baculovirus expression vector system) . This system produces a large amount of the desired proteins by infecting insect cells with a virus encoding the desired transgene under the robust baculovirus polyhedrin promoter .
Post-translational modification (PTM) is one of three levels of protein activity regulation that could efficiently drive adaptive cellular responses by adding or removing functional groups from protein residues . In addition to acetylation and methylation as the variety of PTM, reversible phosphorylation is an important regulatory mechanism that occurs in both prokaryotic and eukaryotic organisms . Reversible phosphorylation results in a conformational change affecting some regulatory processes such as biological thermodynamics of energy-requiring reactions, enzyme inhibition and activation, protein interaction via domain recognition, and protein degradation [12,13,14].
Major efforts have been invested to study the phosphorylation in several species including insects. It has been reported that BRO proteins of BmNPV undergo phosphorylation when the virus infects the host cells and reaches the maximum of phosphorylation level between 14 and 20 h p.i. which was proposed that it has effects on DNA and RNA binding activity regulations . Although phosphorylation has been examined for its crucial role in the regulation of metabolism, the investigation on the protein for its phosphorylation in silkworm is yet lacking. Moreover, BmNPV is not just the most influential factor on B. mori but also useful vector on BEVS system. Hence, the study of the effect of infection on the global phosphoproteome pattern in B. mori is important to comprehend and discover advanced inventions in both fields. Employing TMT (tandem mass tag) labeling and phosphorylation affinity enrichment followed by high-resolution LC-MS/MS analysis, we analyzed the quantitative phosphoproteome in the pair of experimental conditions (control and infected by BmNPV) to elucidate the comprehensive profiling of phosphorylation on the proteome of B. mori after the infection of BmNPV. Altogether, 6480 phosphorylation sites in 2112 protein groups were identified, among which 4764 sites in 1717 proteins were quantified. Some proteins of BmNPV also underwent post-translational modification such as p6.9 and 39 K proteins which interestingly have 2 and 4 up-regulation sites, respectively. Thus, our results will provide a fundamental resource for the prospective study of the phosphorylation in B. mori.
Cell culture and viral infection
B. mori cells, originating from the ovary of silkworm, are maintained in our laboratory. Cells were grown at 27 °C in Sf-900 medium (Thermo Fisher Scientific, USA) complemented with 10% fetal bovine serum. The BmN cells were cultured at a density of 1 × 106 cells in 25 cm2 flask. In the present study, two categories of cells were examined: control cells (normal cell) and those that were infected with the preserved virus, B. mori nucleopolyhedrovirus (BmNPV), in 24 h using 10 MOI (multiplicity of infection). The cell culture and viral infection experiments were done as biological triplicates. After infection, the cells were collected and centrifuged at 3000 rpm. The control cells were treated similarly. Both pellets were washed twice with PBS after the supernatant was discarded.
The samples were mildly sonicated three times on ice employing a high intensity ultrasonic processor (Scientz, Ningbo Scientz Biotechnology Co., Ltd) in lysis buffer (8 M urea, 2 mM EDTA, 10 mM dithiothreitol [DTT], 2% phosphatase inhibitor cocktail V). Then, the samples were centrifuged at 20000 g at 4 °Cfor 10 min to eliminate the remaining debris. The proteins were precipitated using cold 15% trichloroacetic acid (TCA) for 2 h at −20 °C. The supernatant was discarded after the proteins were pelleted by centrifugation at 4 °Cfor10 min. The remaining precipitate was washed three times with cold acetone. The protein extract was re-suspended in buffer (8 M urea, 100 mM NH4HCO3, pH 8.0). The 2-D Quant kit (GE health care) was used to estimate the protein concentration according to the manufacturer’s instructions.
The protein solution was treated with 10 mM DTT for 1 h at 37 °C to obtain a reduced form followed by alkylation using 20 mM iodoacetamide (IAA) for 45 min at room temperature in darkness. Subsequently, for trypsin digestion, the protein was diluted with 100 mM triethylammonium bicarbonate (TEAB) to urea concentration less than 2 M. For the first digestion, the protein was added to trypsin at 1:50 trypsin-to-protein mass ratio overnight and for the second digestion at 1: 100 trypsin-to-protein mass ratio for 4 h.
Following the trypsin digestion, the peptides then were desalted by Strata-X C18 SPE column (Phenomenex) and vacuum-dried. After reconstitution in 0.5 M TEAB, the peptides were labeled using the tandem mass tag (TMT) according to the manufacturer’s instructions (Thermo Scientific) of the 6-plex TMT kit. Concisely, one unit of TMT reagent (defined as the amount of reagent required to label 1 mg of protein) was liquefied and reconstituted in ACN. The peptide mixtures were incubated for 2 h at room temperature, pooled, desalted, and dried by vacuum centrifugation.
High pH reverse-phase HPLC with Agilent 300 Extend C18 column (5 μm particles, 4.6 mm ID, 250 mm length) was used to fractionate the peptide samples. The peptides were separated with a gradient of 2–60% acetonitrile in 10 mM ammonium bicarbonate pH 10 over a period of 80 min into 80 fractions. The peptides then were combined into 8 fractions and dried by vacuum centrifugation.
First, the peptide mixtures were incubated with IMAC microspheres suspension under vibration. The centrifugation precipitated the IMAC microspheres enriched with phosphopeptides. The IMAC microspheres were then washed with 50% CAN/6% TFA and 30% CAN/0.1% TFA, sequentially, in order to remove the nonspecifically absorbed peptides. The elution buffer containing 10% NH4OH was added to elute the enriched phosphopeptides from the IMAC microspheres with vibration. The supernatant containing phosphatides were collected and lyophilized for LC-MC/MS analysis.
The peptide samples were thawed in solvent A (0.1% FA in 2% CAN) and directly loaded onto a reversed-phase pre-column (Acclaim PepMap 100, Thermo Scientific). A reversed-phase analytical column (Acclaim PepMap RSLC C18 column, 50 μm × 15 mm, 2 μm, 100 Å, Thermo Fisher Scientific) was used to separate the peptides with a linear gradient of 4–22% solvent B (0.1% FA in 98% ACN) for 50 min, 22–35% solvent B for 12 min, 35–85% solvent B for 4 min, and holding at 85% for the last 4 min at a constant flow rate of 300 nL/min on an EASY-nLC 1000 UPLC system. The peptides analysis was conducted using Q Exactive™ Plus hybrid Quadrupole-Orbitrap mass spectrometer (ThermoFisher Scientific). The peptides were subjected to NSI (nanospray ion) source followed by tandem mass spectrometry (MS/MS) in Q Exactive™ Plus coupled online to the UPLC. The intact peptides were detected in the orbitrap at a resolution of 70,000. The peptides were selected for MS/MS using the NCE setting as 28 and ion fragments were detected in the Orbitrap at a resolution of 17,500. The top 20 precursor ions above a threshold ion count of 5.0E3 were applied with a data-dependent procedure which alternated between 1 MS scan followed by 20 MS/MS scans in the MS survey scan with 15.0 s dynamic exclusion. The applied electrospray voltage was 2.0 kV. Overfilling of the Orbitrap was prevented by employing an automatic gain control (AGC), and 5E4 ions were accumulated for the generation of MS/MS spectra. The fixed first mass was set as 100 m/z, and the MS scan range was 350–1800 m/z.
MS/MS data analysis
The MS/MS data were analyzed using MaxQuant with integrated Andromeda search engine (v.188.8.131.52). The tandem mass spectra was searched against B. mori database concatenated with reverse decoy database. The peptides were cleaved by trypsin/P enzyme which allowed up to 2 missing cleavages, 5 modifications per peptide, and 5 charges. The mass error was set to 10 ppm for precursor ions and 0.02 Da for the fragment ions. The oxidation of methionine, the phosphorylation on serine, threonine, or tyrosine, and acetylation on the N-terminal of the protein were determined as variable modifications while carbamidomethylation on cysteine was determined as fixed modification. The false discovery rate (FDR) thresholds for protein, peptide, and modification site were specified at 1%. The minimum peptide length was set at 7 and TMT-6plex was selected for quantification. All the other parameters in MaxQuant were configured to default values. The site localization probability was set as >0.5.
Gene Ontology (GO) annotation proteome was derived from the UniProt-GOA database [http://www.ebi.ac.uk/GOA/]. The identified protein ID was first converted to UniProt ID and then mapped to GO IDs. The InterProScan software was employed to annotated protein’s GO functional database based on the sequence alignment method in case some identified proteins were not annotated by UniProt-GOA database. The proteins were then classified by Gene Ontology annotation based on three categories: biological processes, cellular components, and molecular functions. The KEGG (Kyoto encyclopedia of genes and genomes) pathway annotation was performed using KEGG online service tools KAAS according to the annotated protein’s KEGG database description. The domain functional description annotation was conducted by InterProScan based on the protein sequence alignment method and the InterPro domain database.
A functional enrichment analysis was performed based on the three types of annotation (GO, KEEG pathway, and domain annotation). For each category, a two-tailed Fisher’s exact test was employed to test the enrichment or depletion of the differentially expressed protein against all the resulting protein clusters. Standard FDR control methods were used to carry out the correction of multiple hypothesis testing. Any category with a p-value <0.05 in any of the clusters was considered to be significant.
For predicting the subcellular localization, we used WLF PSORT, an updated version of PSORT/PSORT II for the prediction of eukaryotic sequences. Soft motif-x was used to carry out the analysis of sequences model comprised of the amino acids in specific positions of modifying 21-mers (6 amino acids upstream and downstream of the site) in all protein sequences.
Quantification of phosphoproteome
Using TMT labeling and phosphorylation affinity enrichment followed by high-resolution LC-MS/MS analysis, we performed the very first quantitative phosphoproteome analysis aiming to compare the global pattern of phosphoproteome of B. mori after infection of BmNPV (Fig. 1a). Focusing on those phosphoproteins that were either up- or down-regulated after BmNPV infection, we identified 6480 phosphorylation sites in 2112 protein groups, among which 4764 sites in 1717 proteins were quantified. The fold-change cut-off was set when proteins with the quantitative ratios above 1.3 or below 0.77 were deemed to be significant. Among the quantified proteins, 81 hyperphosphorylated sites in 66 proteins and 25 hypophosphorylated sites in 24 proteins were identified after BmNPV infection. The distribution of peptides length from the prepared sample was between 8 and 20 (Fig. 1b), which was complied with the standard and in agreement with the property of tryptic peptides. The mass error of the peptide weight was almost zero, and most of them were less than 0.02 Da (Fig. 1c), indicating that the accuracy of MS data was as required.
Among the regulated phosphoproteins, some proteins possessed a remarkable elevation ratio, for example, origin recognition complex subunit 2 with the highest phosphorylation ratio 5.69 among the hyperphosphorylated proteins and uncharacterized protein LOC101739044 with the lowest phosphorylation ratio 0.37 among the hypophosphorylated proteins. Moreover, other proteins also occurred with more than one regulated site, i.e., thymosin isoform XI with 4up-regulated sites and upstream activation factor subunit spp27 with 2down-regulated sites (Table 1). There are also some notable proteins of BmNPV which were modified during infection. For instance, P6.9 has two up-regulation sites and 39 K protein has 4 up-regulation sites, one of which had a very high phosphorylation ratio of 16.68 fold. Other modulated BmNPV proteins also have more than one hyperphosphorylated site such as LEF-6, Ac58-like protein, Ac82-like protein, and BRO-D in most of which a relative high phosphorylation ratio was observed (Table 2). Those proteins with the recognizable characteristic either with the remarkable phosphorylation ratio or more than one regulated sites are strongly suggested to have a significant role in the alteration of several biological functions in silkworm cells after BmNPV infection.
The GO annotation classified either the up-regulated or down-regulated proteins into three categories: biological processes, cellular components, and molecular functions. The up-regulated proteins in the biological process category were mainly identified at metabolic (9 proteins) and cellular processes (8 proteins); the down-regulated proteins were also mostly identified at cellular (7 proteins) and metabolic processes (6 proteins). The pie chart of the modulated proteins with respect to the biological processes, either the up- or down-regulated protein (Fig. 2a), showed that proteins distributed in metabolic and cellular processes have the same percentage (26%). Regarding cellular components, the up-regulated proteins were highly identified at cell (7 proteins) while the down-regulated proteins were apparently distributed evenly in the three components: membrane, organelle, and cell (2 proteins in each). Altogether, the distribution of up- or down-regulated proteins in the term of cellular component was mostly identified in the cell at 37% (Fig. 2b). Moreover, the most abundant modulated proteins in terms of molecular function were distributed in the binding activity: 32 up-regulated and 12down-regulated proteins, which amounted to 65% (Fig. 2c), the highest among the three main categories. These results strongly suggested that the binding activity in silkworm cell after BmNPV infection may be comprehensively regulated by phosphorylation. In addition, phosphorylation also seemed to exhibit an essential regulatory role in the cellular and metabolic processes.
The functional enrichment was carried out in three types of analysis, i.e., gene ontology, domain, and KEGG pathway analysis. The GO enrichment analysis showed four modulated molecular functions listed in Fig. 3a, which indicated that most of them were involved in the binding activity. The four domains were reported to be up-regulated by domain enrichment analysis (Fig. 3b) while only bmor04214 apoptosis-fly pathway in the silkworm cells was enriched in the KEGG pathway analysis (Fig. 3c). The subcellular localization classification revealed that phosphorylated proteins in the treated B. mori cells compared with the control group were primarily located in the nucleus: either up-regulated (40 proteins) or down-regulated proteins (13 proteins). In addition, the cytosolic phosphorylated proteins were also highly identified in both up-regulated (13 proteins) and down-regulated (6 proteins) phosphoproteome (Table 3). The motif analysis was carried out to identify the over represented residues which are amino acids that often found both upstream and downstream of a modification site. Using soft motif-x resulting in four distinct phosphorylation motifs (Fig. 4), the pattern of amino acids distribution was observed from position −6 to +6 around the phosphorylation sites. We found that the hydrophobic amino acid, proline, was drastically over represented at position +1 following the phosphorylated serine at each of the six patterns. Moreover, proline was also enriched at +2 position in motif 3 indicating that the serine-rich region and hydrophobic amino acid (proline) region of the protein were preferred for phosphorylation. In addition, some basic amino acid was also over represented, i.e., arginine in −3 position at motif 1 and lysine in +3 position at motif 2 implying that polar, hydrophobic, and positively charged amino acids may be functionally critical for phosphorylation.
The comprehensive phosphoproteome quantification has been extensively studied in several species attributable to its pivotal role leading to the altered gene expression through conformational changes in some functional proteins. B. mori and BmNPV were a topic of intense focus among several groups due to their correlation in the interaction of insect host and baculovirus pathogen with respect to the mechanism of virus invasion of the host  and the corresponding host immune system responses . The present study analyzed the impact of infection of BmNPV to the global pattern of proteomic phosphorylation in silkworm cells. However, to our knowledge, this is the first study determining the phosphoproteome of B. mori after BmNPV infection resulting in a significant amount of candidate target phosphorylated proteins for deeper insight into the mechanism of phosphorylation to control the alterations on the host cells during viral infection.
Overall, the results showed that BmNPV infection greatly influenced the phosphorylation of diverse proteins related to several functions (Table 1). The notable alterations on the phosphorylation of proteins were involved in binding activity, protein synthesis, viral replication, and apoptosis through kinase activity and are further discussed below.
Phosphorylation may have a pivotal function on the regulation of binding activity
As shown by the GO annotation, the most significant altered phosphoproteins (65%) were found in the binding activity for molecular function, either up-regulated or down-regulated (Fig. 2c). Among those modulated proteins, 34% were involved in nucleic acid binding; 10 proteins were up-regulated and 5 proteins were down-regulated. Although the phosphorylation ratio was little, more than 1/3 of the total number of proteins in the binding activity were involved in nucleic acid binding activity, indicating that phosphorylation may have a significant role in either the DNA or RNA binding activity in the silkworm cells after BmNPV infection. For instance, DNA topoisomerase II, one of the up-regulated proteins with 1.43 ratio phosphorylation, is a well-known enzyme in the process of transcription, replication, chromatin remodeling, and segregation . This enzyme has been reported to have the ability to simplify the DNA topology by relaxing supercoiled DNA and unlinking knotted or catenated DNA which is in accordance with its physiological role in DNA replication and chromosome segregation [19, 20]. Moreover, the enzyme is also a predominant component of the nuclear scaffold and matrix [18, 21]. The matrix associated with the DNA regions was also reported to consist of topoisomerase II binding sites . According to a search of the BmNPV genome sequence, it does not encode an ortholog of topoisomerase II. Furthermore, a recent report has proposed a plausible hypothesis that the intertwined baculovirus DNA was decatenated by topoisomerase 2 during viral replication . Hence we hypothesize that the virus employs this enzyme for its genome replication by activating the enzyme through up-regulation. Study about topoisomerase for baculovirus is still lack. Further studies are necessitated to substantiate our hypothesis.
Another remarkable protein is the subunitssp27, a component of the upstream activation factor (UAF) complex, which interestingly has two phosphorylation sites in the position 84th and 88th of amino acid sequence with the elevation ratio 0.57- and 0.73-fold, respectively. Although in B. mori this protein is not yet examined, in yeast Saccharomyces cerevisiae, the upstream activation factor (UAF) is well-known as a multi-protein complex that tightly binds to the upstream element of the rDNA promoter and also strongly stimulates the transcription of RNA polymerase I (Pol I) . Subunit spp27 of this protein is a homolog from Schizosaccharomyces pombe which has UAF-like activity  such as the ability to interact with rDNA upstream sequence binding protein; Acr1, an accumulation of condensing at rDNA 1 which has been reported to be weakly similar to RRN9, a component of Saccharomyces cereviseae UAF complex . Subunit spp27 also has been reported to have good alignment with UAF30p, another subunit of UAF complex of S. cereviseae . Furthermore, UAF also plays a critical role in the silencing of transcription through RNA polymerase II (Pol II) by stabilizing Pol I state . The down-regulation of this protein may cause the inhibition of the silencing, which leads to the activation of transcription by Pol II. This activation of Pol II might correspond with the BmNPV-based manipulation of the host gene by miRNA that is canonically transcribed by Pol II. By encoding a miRNA, the BmNPV represses the expression of the silkworm GTP-binding nuclear protein Ran, resulting in the decrease of the host small-RNA population . The down-regulation of UAF causing all those possible consequents which are eventually leading to the alleviation of the amount of host small-RNA might be required for the enhancement of BmNPV load in the infected cells.
Another intriguing result is that thymosin isoform XI (tiXI) is an up-regulated protein with four sites exhibiting high phosphorylation ratio. Those four sites were located at position 106th, 30th, 31st, and 107th on the amino acid sequence and the phosphorylation elevation ratios compared with the control group were 1.99, 2.0, 2.9, and 1.82, respectively. Although tiXI is not widely focused, the family of this protein, thymosin beta 4 (Tβ4), has been extensively studied in several species and is known to have a pivotal role in the organization of cytoskeleton by binding or sequestering to actin monomers . This protein regulates the actin cytoskeleton dynamics through retaining a large pool of actin monomers which interact with globular actin to produce a 1:1 complex for controlling the filamentous actin assembly . Tβ4 regulates the equilibrium between the globular and filamentous actin, which is vital for rapid rearrangement of the cytoskeleton . The hyperphosphorylation of this protein after BmNPV infection might activate the cellular Tβ4 protein which is responsible to the filamentous actin to globular actin ratio alteration . This regulation appears close relation with the infection process of BmNPV. The actin and microtubule have significant roles in every virus life cycle. Viruses have ability to reorganize the cytoskeleton and restructure the host transport equipment to fulfill their needs , for instance, H5N1 could manipulate and exploit diverse host cytoskeletal protein to promote their infection . AcMNPV has been reported to harness actin polymerization based motility in two infection phases. In the initial phase right after cell entry, they utilize motility for cytoplasm exploration and translocation to the nucleus through nuclear periphery for initiation the viral gene. In the next phase, after early gene expression, AcMNPV required motility to accumulate the distinct set of nucleocapsid at the tips of actin-rich surface spikes . The up-regulation occurred on four sites of this tiXI protein with a high phosphorylation ratio after BmNPV infection strongly suggested that phosphorylation might have a predominant influence on cytoskeleton change through actin binding which is important for life cycle of every virus. Further studies are highly recommended to authenticate our hypothesis.
Phosphorylation may represent an obligatory regulation in protein synthesis
Similar to other viruses, baculoviruses employ the host chaperones to assist the rapid synthesis of a large quantity of their viral proteins . Baculoviruses completely rely on the host translational equipment such as ribosomes, tRNA, amino acid metabolism and transport, chaperones for protein folding and other translation factors. The presence of AT-rich regions and unstructured 5′ untranslated (UTR) regions of the virus mRNA allow the baculovirus to compete with host mRNA for translation by its higher binding activity to host translation factor . A previous study revealed some translational factors in silkworm cells that were commonly up-regulated after BmNPV infection such as eukaryotic initiation factor (eIF3–6, eIF1A, aIF3-2b) and an elongation factor (EF1d)  corresponding to the present study, resulting in the hyperphosphorylation of EF1d with fold elevation ratio 1.69. Many studies have reported that the EF1d phosphorylation could alter translational efficiency [37,38,39]. Furthermore, in mammalian cell infected by herpes viruses, EF1d also phosphorylated by conserved protein kinases encoded by these viruses suggesting that it plays an important role of efficient replication of these viruses . The up-regulation of BmNPV EF1d suggested that this virus might employ EF1d and enhance its translational efficiency for viral replication through phosphorylation.
Another notably up-regulated phosphoprotein is ATP-dependent RNA helicase DDX10. RNA helicase is well-known as an essential factor in most of the RNA metabolism processes including ribosome biogenesis, pre-mRNA splicing and translation initiation [41, 42]. DDX10 is categorized in the group of DEAD-box family, the largest group of helicases . Although this protein from silkworm is yet less studied, the other proteins in the family from different species have been elucidated for the correlation with the viral infection. DDX1 is the common helicase from DEAD-box family, which has been reported as a cellular co-factor for viral replication by several viruses such as coronavirus , JC virus , hepatitis C virus , and also HIV-1 for replication and nuclear export . With a high fold elevation ratio (2.31), the hyperphosphorylation of ATP-dependent RNA helicase DDX10 protein is designated to play a major role in protein synthesis in silkworm cells after BmNPV infection. Thus, further studies are imperative to better understand the role of hyperphosphorylation in the ATP-dependent RNA helicase DDX10 for protein synthesis of B. mori cells after BmNPV infection.
Phosphorylation may be responsible for viral replication during infection
During infection, while the host DNA replication is detained, viral replication is rapidly conducted harnessing the host replication machinery and concentrating in a virogenic stroma, which is a separate compartment within the nucleus [48, 49]. DNA ligase 1-like with the fold elevation ratio 1.44, is one of the up-regulated proteins that is associated with the replication of DNA. DNA ligase has been highlighted as a crucial factor in the linkage of Okazaki fragment on the lagging strand . Interestingly, the Viral DNA ligase of vaccinia virus has been reported to harness cellular Topoisomerase 2 to site the viral replication and assembly . However, BmNPV apparently does not encode DNA ligase, only some baculoviruses do, such as LdMNPV (Lymantria dispar multicapsid nuclear polyhedrosis virus), according to the BLAST result, its DNA ligase is most optimally aligned with the B. mori DNA ligase 3 isoform X1 after other lepidopteran suggesting that the DNA ligase of LdMNPV was probably acquired from host. The fact that BmNPV genome does not have the ortholog of DNA ligase suggested that the up-regulation of cellular DNA ligase may be important either for viral replication or assembly.
The origin recognition complex subunit 2 (orc2), with a high phosphorylation ratio (5.69), is the most hyperphosphorylated protein for DNA replication and the highest up-regulated protein among all the different classes of proteins. Although the mechanism in detail remains elusive, the viral replication of BmNPV has been reported to affect the expression of BmORCs . In the case of human, ORC is known as important protein which is able to initiate DNA replication by facilitating the establishment of the pre-replication complex (pre-RC) at the origin of DNA replication . ORC is also known to be associated with an EBV (Epstein-Barr virus) replication origin for the propagation of its genome by utilizing antibodies against the three different subunits of human ORC in order to precipitate the cross-linked chromatin . The up-regulation of B. mori orc2 protein after BmNPV infection was probably important for the activation of this protein. Together with DNA ligase, we speculated that they have an essential function for the viral genome replication of BmNPV. Further research is highly recommended to ascertain this hypothesis.
Phosphorylation may regulate anti-apoptotic cells through the kinases activity after BmNPV infection
Apoptosis is a genetically programmed cell death as an emergency response to a variety of stimuli such as radiation damage, aberrant growth or viral infection . Some viruses such as Oropouche virus (OROV) induce apoptosis in the cells right after their infection that causes degeneration in cells shortly after infection . However, some viruses evolve strategies to evade early apoptosis in order to avoid the elimination of progeny virus spread and the limitation of virus production . For instance, an inhibitor of apoptosis proteins (IAPs), a family of proteins which were originally discovered in baculovirus has been reported to be involved in suppressing the host cell death in response to viral infection [57, 58]. Protein kinases C (PKC) family members are associated with an extensive range of cellular responses including cell permeability, migration, hypertrophy, proliferation, secretion, and also apoptosis. In particular, PKC cascades are apparently essential modulators of the apoptotic response . A recent study showed a high up-regulated phosphorylation of PKC in the silkworm cells after BmNPV infection. The hyperphosphorylated PKC occurred in the epsilon isoform (PKCε). Among several members of PKC protein family, some isoforms induceapoptosis, and the others inhibit it. PKCε possesses the ability to enhance tumorigenesis that can associate with the suppression of apoptosis . PKCε has also been reported to trigger cellular survival by initiating the NF-κB pathway , which leads to the activation of gene transcription encoding IAPs . With PKCε high phosphorylation ratio (3.53), we strongly suggested that the hyperphosphorylation of this protein has a great implication on the survival of host cells during BmNPV infection.
The phosphorylation of BmNPV protein
All of the quantified proteins of BmNPV are up-regulated; some of them have a notable phosphorylation ratio and are even at more than one site. 39 K protein, also known as Bmorf 27, is one of the prominent hyperphosphorylated proteins of BmNPV which has 4 phosphorylated sites and one of them has the highest phosphorylation ratio among all of BmNPV phosphorylated protein and is even higher than it of B. mori phosphorylated protein. It has been well known that 39 K promoters have important utility for insect cell engineering such as BEV system by providing the highest level of tightly regulated protein production . It also has been reported that BmNPV 39 K plays an important role in viral late gene transcription. Moreover, the deletion of this gene resulted in very poor budded virus production, low polyhedral production, oral infectivity attenuation, and deprived virogenic stroma formation . The homolog of 39 k, AcMNPV orf 36 also called as pp31 has been reported to be phosphorylated as well . Although the function of pp31 AcMNPV phosphorylation remains elusive, the great phosphorylation ratio of 39 K phosphorylation strongly suggested that phosphorylation may have important role for this protein in promoting protein production. Further study is necessary for supporting this hypothesis.
BmNPV protamine-like protein named P6.9 underwent phosphorylation during infection; it has two phosphorylated sites with phosphorylation ratio of 2.73 and 2.07, respectively. Although in BmNPV this protein has not been well assessed but in Autographa californica multiple nucleopolyhedrovirus (AcMNPV), P6.9 has been reported to be transiently phosphorylated in infected cells prior to nucleocapsid assembly . This protein binds and condenses baculoviral DNA for packing into capsids . By the time nucleocapsids enter nuclei through nuclear pores, P6.9 is phosphorylated which could repel the negatively charged DNA, leading to the viral genome release associated with capsid and kinase [68, 69]. A recent study indicated that during infection there are 13 Ser/Thr phosphorylation sites on P6.9 AcNMPV which interestingly 7 of those are reported to be dependent to PK1 (protein kinase 1). Moreover, replacement of those 7 Ser/Thr phosphorylated residues with Ala in P6.9 AcNMPV significantly diminishes the transcription of the very late viral genes and viral infectivity . Further research is necessary to find whether phosphorylation of P6.9 BmNPV has an important role on DNA binding activity and if PK1 also has influence to the phosphorylation of this protein.
Another remarkable protein of BmNPV is LEF-6 with 5 phosphorylated sites; the phosphorylation ratio in every site was relatively high. In AcMNPV lef-6 has been reported to show ability in accelerating the viral lifecycle and increasing virus yields . The study on this protein phosphorylation is still lacking. However, considering its 5 hyperphosphorylated sites, we strongly suggest that phosphorylation on this protein may play significant role in the viral yielding and lifecycle.
The current study comprehensively identified the global pattern of phosphorylation of B. mori after BmNPV infection, harnessing a highly sensitive proteomic method. In total, 4764 phosphorylation sites in 1717 proteins were quantified in silkworm cells, which were extensively involved in the regulation of several aspects of metabolism and molecular function during BmNPV infection. These results also suggest a large number of candidate target phosphorylated proteins for in-depth studies on the mechanism of phosphorylation modulating the alterations on the host cells during viral infection. Furthermore, we classified and discussed most of the differentially and remarkably phosphorylated proteins into several functional groups, i.e., binding activity, protein synthesis, viral replication, and apoptosis through kinase activity. There are also remarkable proteins of BmNPV which were modified during infectionsuch as P6.9, 39 K, LEF-6, Ac58-like protein, Ac82-like protein, and BRO-D among which most of those sites have relative high phosphorylation ratio and have more than one hyperphosphorylated site. Nevertheless, there will be undeniable cross-talk between these groups despite the phosphorylated proteins are involved in multiple molecular functions and cellular processes.
Autographa californica multiple nucleopolyhedrovirus
Automatic gain control
Baculovirus expression vector system
B. mori nucleopolyhedrovirus
DNA binding protein
Eukaryotic initiation factor
The false discovery rate
Inhibitor of apoptosis proteins
Kyoto encyclopedia of genes and genome
Lymantria dispar multicapsid nuclear polyhedrosis virus
Multiplicity of infection
The origin recognition complex subunit 2
Protein kinase C epsilon
- Pol I:
RNA polymerase I
- Pol II:
RNA polymerase II
Trichloro acetic acid
Triethyl ammonium bicarbonate
Thymosin isoform XI
Tandem mass tag
Upstream activation factor
Xia Q, Zhou Z, Lu C, Cheng D, Dai F, Li B, et al. A draft sequence for the genome of the domesticated silkworm (Bombyx mori). Science. 2004;306:1937–40.
Jiang L, Xia Q. The progress and future of enhanching antiviral capacity by transgenic technology in the silkworm Bombyx mori. Insect Biochem Mol Biol Rev. 2014;48:1–7.
Sekimizu N, Paudel A, Hamamoto H. Animal welfare and use of silkworm as a model animal. Drug Discov Ther. 2012;6(4):226–9.
Mita K. Genome of a lepidopteran model insect, the silkworm Bombyx mori. Seikagaku. 2009;81(5):353–60.
Mon H, Izumi M, Mitsunobu H, Tatsuke T, Liyama K, Jikuya H, et al. Post-translational modification of the N-terminal tail of histone H3 in holocentric chromosomes of Bombyx mori. Insect Biochem Mol Biol. 2011;41:902–8.
Nie Z, Zhu H, Zhou Y, Wu C, Liu Y, Sheng Q, et al. Comprehensive profiling of lysine acetylation suggest the widespread function is regulated by protein acetylation in the silkworm, Bombyx mori. Proteomics. 2015;15:3253–66.
Jiang L, Zhao P, Cheng T, Sun Q, Peng Z, Dang Y, et al. A transgenic animal with antiviral properties that might inhibit multiple stages of infection. Antivir Res. 2013;98:171–3.
Lee TI, Young RA. Transcription of eukaryotic protein-coding genes. Annu Rev Genet. 2000;34:77–137.
Lee DF, Chen CC, Hsu TA, Juang JL. A baculovirus super infection system: efficient vehicle for gene transfer into drosophila S2 cells. J Virol. 2000;74:11873–80.
Deribe YL, Pawson T, Dikic I. Post-translational modifications in signal integration. Nat Struct Mol Biol. 2010;17:666–72.
Chang C, Stewart RC. The two-component system. Regulation of diverse signaling pathways in prokaryotes and eukaryotes. Plant Physiol. 1998;117(3):723–31.
van Weeren PC, de Bruyn KM, de Vries-Smits AM, van Lint J, Burgering BM. Essential role for protein kinase B (PKB) in insulin-induced glycogen synthase kinase 3 inactivation. Characterization of dominant-negative mutant of PKB. J Biol Chem. 1998;273:13150–6.
Cole PA, Shen K, Qiao Y, Wang D. Protein tyrosine kinases Src and Csk: a tail's tale. Curr Opin Chem Biol. 2003;7:580–5.
Babior BM. NADPH oxidase: an update. Blood. 1999;93(5):1464–76.
Zemskov EA, Kang W, Maeda S. Evidence for nucleic acid binding ability and nucleosome association of Bombyx mori nucleopolyhedrovirus BRO proteins. J Virol. 2000;74(15):6784–9.
Rahman MM, Gopinathan KP. Systemic and in vitro infection process of Bombyx mori nucleopolyhedrovirus. Virus Res. 2004;101:109–18.
Xue J, Qiao N, Zhang W, Cheng RL, Zhang XQ, Bao YY, et al. Dynamic interaction between Bombyx mori nucleopolyhedrovirus and its host cells revealed by transcriptome analysis. J Virol. 2012;86:7345–59.
Earnshaw WC, Heck MM. Localization of topoisomerase II in mitotic chromosomes. J Cell Biol. 1985;100:1716–25.
Liu LF, Liu CC, Alberts BM. Type II DNA topoisomerases: enzymes that can unknot a topologically knotted DNA molecule via a reversible double-strand break. Cell. 1980;19(3):697–707.
Ullsperger CJ, Vologodskii AV, Cozzarelli NR. Unlinking of DNA topoisomerase during DNA replication. In: Eckstein F, David M, Lilley J, editors. Nucleic acid and molecular biology. Berlin: Springer; 1995. p. 115.
Berrios M, Osheroff N, Fisher PA. In situ localization of DNA topoisomerase II, a major polypeptide component of the drosophila nuclear matrix fraction. Proc Natl Acad Sci U S A. 1985;82:4142–6.
Wu Y, Wu Y, Wu Y, Tang H, Wu H, Zhang G, et al. Screening of candidate proteins interacting with IE-2 of BmNPV. Mol Biol Rep. 2013;40:5797–804.
Volkman LE. Baculoviruses and nucleosome management. Virology. 2015;476:257–63.
Vu L, Siddiqi I, Lee BS, Josaitis CA, Nomura M. RNA polymerase switch in transcription of yeast rDNA: role of transcription factor UAF (upstream activation factor) in silencing rDNA transcription by RNA polymerase II. Proc Natl Acad Sci U S A. 1999;96:4390–5.
Liu M, Guo A, Boukhgalter B, Van Den Heuvel K, Tripp M. Et el. Characterization of the fission yeast ribosomal DNA binding factor: components share homology with upstream activating factor and with SWI/SNF subunits. Nucleic Acids Res. 2002;30(24):5347–59.
Nakazawa N, Nakamura T, Kokubu A, Ebe M, Nagao K, Yanagida M. Dissection of the essential steps for condensin accumulation at kinetochores and rDNAs during fission yeast mitosis. J Cell Biol. 2008;180(6):1115–31.
Singh CP, Singh J, Nagaraju J. A baculovirus-encoded microRNA (miRNA) suppresses its host miRNA biogenesis by regulating the exportin-5 cofactor ran. J Virol. 2012;86:7867–79.
Sanger JM, Golla R, Safer D, Choi JK, Yu KR, Sanger JW, et al. Increasing intracellular concentrations of thymosin beta 4 in PtK2 cells: effects on stress fibers, cytokinesis, and cell spreading. Cell Motil Cytoskeleton. 1995;31:307–22.
Safer D, Elzinga M, Nachmias VT. Thymosin beta 4 and Fx, an actin-sequestering peptide, are indistinguishable. J Biol Chem. 1991;266(7):4029–32.
Grant DS, Kinsella JL, Kibbey MC, LaFlamme S, Burbelo PD, Goldstein AL, et al. Matrigel induces thymosin beta 4 gene in differentiating endothelial cells. J Cell Sci. 1995;108(12):3685–94.
Kim S, Kwon J. Actin cytoskeletal rearrangement and dysfunction due to activation of the receptor for advanced glycation end products is inhibited by thymosin beta 4. J Physiol. 2015;593(8):1873–86.
Radtke K, Döhner K, Sodeik B. Viral interactions with the cytoskeleton: a hitchhiker's guide to the cell. Cell Microbiol. 2006;8(3):387–400.
Li Y, Ming F, Huang H, Guo K, Chen H, Jin M, et al. Proteome response of chicken embryo fibroblast cells to recombinant H5N1 avian influenza viruses with different neuraminidase stalk lengths. Sci Rep. 2017;12(7):40698.
Ohkawa T, Volkman LE, Welch MD. Actin-based motility drives baculovirus transit to the nucleus and cell surface. J Cell Biol. 2010;190(2):187–95.
Champoux JJ. DNA topoisomerases: structure, function, and mechanism. Annu Rev Biochem. 2001;70:369–413.
Scheper GC, Vries RG, Broere M, Usmany M, Voorma HO, Vlak JM, et al. Translational properties of the untranslated regions of the p10 messenger RNA of Autographa californica multicapsid nucleopolyhedrovirus. J Gen Virol. 1997;78:687–96.
Chang YW, Traugh JA. Phosphorylation of elongation factor 1 and ribosomal protein S6 by multipotential S6 kinase and insulin stimulation of translational elongation. J Biol Chem. 1997;272:28252–7.
Minella O, Cormier P, Morales J, Poulhe R, Belle R, Mulner-Lorillon O. cdc2 kinase sets a memory phosphorylation signal on elongation factor EF-1 delta during meiotic cell division, which perdures in early development. Cell Mol Biol (Noisy-Le-Grand). 1994;40:521–5.
Venema RC, Peters HI, Traugh JA. Phosphorylation of elongation factor 1 (EF-1) and valyl-tRNA synthetase by protein kinase C and stimulation of EF-1 activity. J Biol Chem. 1991;266:12574–80.
Kawaguchi Y, Kato K, Tanaka M, Kanamori M, Nishiyama Y, Yamanashi Y. Conserved protein Kinases encoded by Herpesviruses and cellular protein Kinase cdc2 target the same Phosphorylation site in eukaryotic elongation factor 1δ. J Virol. 2003;77(4):2359–68.
Jankowsky E. RNA Helicases at work: binding and rearranging. Trends Biochem Sci. 2010;36:19–29.
Fairman-Williams ME, Guenther UP, Jankowsky E. SF1 and SF2 helicases: family matters. Curr Opin Struct Biol. 2010;20:313–24.
Gorbalenya AE, Koonin EV. Helicases: amino acid comparisons and structure-function relationships. Curr Opin Struct Biol. 1993;3:419–29.
Xu L, Khadijah S, Fang S, Wang L, Tay FP, Liu DX. The cellular RNA helicase DDX1 interacts with coronavirus nonstructural protein 14 and enhances viral replication. J Virol. 2010;84:8571–83.
Sunden Y, Semba S, Suzuki T, Okada Y, Orba Y, Nagashima K, et al. Identification of DDX1 as a JC virus transcriptional control region-binding protein. Microbiol Immunol. 2007;51:327–37.
Pan TT, Fang CY, Yi ZG, Yang PY, Yuan ZG. Subproteomic analysis of the cellular proteins associated with the 3′ untranslated region of the hepatitis C virus genome in human liver cells. Biochem Biophys Res Commun. 2006;347:683–91.
Edgcomb SP, Carmel AB, Naji S, Ambrus-Aikelin G, Reyes JR, Saphire AC, et al. DDX1 is an RNA-dependent ATPase involved in HIV-1 rev function and virus replication. J Mol Biol. 2012;415:61–74.
Braunagel SC, Parr R, Belyavskyi M, Summers MD. Autographa californica Nucleopolyhedrovirus infection results in Sf9 cell cycle arrest at G2/M phase. Virology. 1998;244:195–211.
Vanarsdall AL, Okano K, Rohrmann GF. Characterization of the replication of a baculovirus mutant lacking the DNA polymerase gene. Virology. 2005;331:175–80.
Nguyen Q, Nielsen LK, Reid S. Genome scale transcriptomics of baculovirus-insect interactions. Viruses. 2013;5:2721–47.
Lin YCJ, Li J, Irwin CR, Jenkins H, DeLange L, Evans DH. Vaccinia virus DNA ligase recruits cellular topoisomerase II to sites of viral replication and assembly. J Virol. 2008;82(12):5922–32.
Yang HP, Luo SJ, Li YN, Zhang YZ, Zhang ZF. Identification and characterization of the DNA replication origin recognition complex gene family in the silkworm Bombyx mori. Biosci Rep. 2011;31(5):353–61.
Shen Z, Prasanth SG. Orc2 protects ORCA from ubiquitin-mediated degradation. Cell Cycle. 2012;11:3578–89.
Chaudhuri B, Xu H, Todorov I, Dutta A, Yates JL. Human DNA replication initiation factors, ORC and MCM, associate with oriP of Epstein-Barr virus. Proc Natl Acad Sci U S A. 2001;98:10085–9.
Teodoro JG, Branton PE. Regulation of apoptosis by viral gene products. J Virol. 1997;71(3):1739–46.
Acrani GO, Gomes R, Proença-Módena JL, da Silva AF, Carminati PO, Silva ML, et al. Apoptosis induced by Oropouche virus infection in HeLa cells is dependent on virus protein expression. Virus Res. 2010;149:56–63.
Crook NE, Clem RJ, Miller LK. An apoptosis-inhibiting baculovirus gene with a zinc finger-like motif. J Virol. 1993;67(4):2168–74.
RJ BMJC, Miller LK. An apoptosis-inhibiting gene from a nuclear polyhedrosis virus encoding a polypeptide with Cys/his sequence motifs. J Virol. 1994;68(4):2521–8.
Utz P, Anderson P. Life and death decisions: regulation of apoptosis by proteolysis of signaling molecules. Cell Death Differ. 2000;7:589–602.
Wheeler DL, Martin KE, Ness KJ, Li Y, Dreckschmidt NE, Wartman M, et al. Protein kinase C epsilon is an endogenous photosensitizer that enhances ultraviolet radiation-induced cutaneous damage and development of squamous cell carcinomas1. Cancer Res. 2004;64:7756–65.
Catley MC, Cambridge LM, Nasuhara Y, Ito K, Chivers JE, Baton A, et al. Inhibitors of protein kinase C (PKC) prevent activated transcription: role of events downstream of NF-KappaB DNA binding. J Biol Chem. 2004;279:18457–66.
Danial NN, Korsmeyer SJ. Cell death: critical control points. Cell. 2004;116:205–16.
Lin CH, Jarvis DL. Utility of temporally distinct baculovirus promoters for constitutive and baculovirus-inducible transgene expression in transformed insect cells. J Biotechnol. 2013;165(1):11–7.
Gomi S, Zhou CE, Yih W, Majima K, Maeda S. Deletion analysis of four of eighteen late gene expression factor gene homologues of the baculovirus, BmNPV. Virology. 1997;230(1):35–47.
Broussard DR, Guarino LA, Jarvis DL. Dynamic phosphorylation of Autographa californica nuclear polyhedrosis virus pp31. J Virol. 1996;70(10):6767–74.
Oppenheimer DI, Volkman LE. Proteolysis of p6.9 induced by cytochalasin D in Autographa californica M nuclear polyhedrosis virus-infected cells. Virology. 1995;207(1):1–11.
Tweeten KA, Bulla LA, Consigli RA. Characterization of an extremely basic protein derived from granulosis virus nucleocapsids. J Virol. 1980;33(2):866–76.
Funk CJ, Consigli RA. Phosphate cycling on the basic protein of Plodia interpunctella granulosis virus. Virology. 1993;193(1):396–402.
Wilson ME, Consigli RA. Functions of a protein kinase activity associated with purified capsids of the granulosis virus infecting Plodia interpunctella. Virology. 1985;143(2):526–35.
Li A, Zhao H, Lai Q, Huang Z, Yuan M, Yang K. Posttranslational modifications of Baculovirus Protamine-like protein P6.9 and the significance of its hyperphosphorylation for viral very late Gene Hyperexpression. Sandri-Goldin RM ed. J Virol. 2015;89(15):7646–59.
Lin G, Blissard GW. Analysis of an Autographa californica multicapsid nucleopolyhedrovirus LEF-6-null virus: LEF-6 is not essential for viral replication but appears to accelerate late gene transcription. J Virol. 2002;76(11):5503–14.
We would like to thank Dr. Chen Bu for his excellent technical assistance in statistical analysis of the data.
This work was supported by Science foundation of Zhejiang province (No. LY17C170006) and the National High-tech R&D program (863 Program) (No. 2011AA100603).
Availability of data and materials
All data and information are available without restriction and included in the manuscript.
WY conceived and designed the experiments. JS and CZ performed the experiments. JS, SX, DH, MW and YQ participated in data analysis. WY and JS drafted the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interest.
Consent for publication
Ethics approval and consent to participate
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.