PubMed@Masakazu SAGA:23580252 JSONTXT

{"target":"https://pubannotation.org/docs/sourcedb/PubMed@Masakazu%20SAGA/sourceid/23580252","sourcedb":"PubMed@Masakazu SAGA","sourceid":"23580252","text":"Shotgun Proteomic Analysis on the Diapause and Non-Diapause Eggs of Domesticated Silkworm Bombyx mori \r\n\r\nTo clarify the molecular mechanisms of silkworm diapause, it is necessary to investigate the molecular basis at protein level. Here, the spectra of peptides digested from silkworm diapause and non-diapause eggs were obtained from liquid chromatography tandem mass spectrometry (LC-MS/MS) and were analyzed by bioinformatics methods. A total of 501 and 562 proteins were identified from the diapause and non-diapause eggs respectively, of which 309 proteins were shared commonly. Among these common-expressed proteins, three main storage proteins (vitellogenin precursor, egg-specific protein and low molecular lipoprotein 30 K precursor), nine heat shock proteins (HSP19.9, 20.1, 20.4, 20.8, 21.4, 23.7, 70, 90-kDa heat shock protein and heat shock cognate protein), 37 metabolic enzymes, 22 ribosomal proteins were identified. There were 192 and 253 unique proteins identified in the diapause and non-diapause eggs respectively, of which 24 and 48 had functional annotations, these unique proteins indicated that the metabolism, translation of the mRNA and synthesis of proteins were potentially more highly represented in the non-dipause eggs than that in the diapause eggs. The relative mRNA levels of four identified proteins in the two kinds of eggs were also compared using quantitative reverse transcription PCR (qRT-PCR) and showed some inconsistencies with protein expression. GO signatures of 486 out of the 502 and 545 out of the 562 proteins identified in the diapause and non-diapause eggs respectively were available. In addition, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis showed the Metabolism, Translation and Transcription pathway were potentially more active in the non-dipause eggs at this stage.\r\n\r\nIntroduction\r\n\r\nDiapause is a special physiological state of arrested development by many insects to avoid unfavorable environments such as low temperature, drought or food shortage and to synchronize their life cycles to these changes. Diapause is widespread among insects and can occur at any stage of the life cycle, i.e., adult, pupa, larva or egg, and each species enters diapause at fixed stages. Diapause is endogenously controlled, and this dormancy typically begins well before conditions become too harsh to support normal growth and development.\r\n\r\nIn the silkworm Bombyx mori, the development of diapause-destined embryos is arrested during the G2 cell cycle stage immediately after formation of the cephalic lobe and telson and sequential segmentation of the mesoderm. In Bombyx mori, embryonic diapause is determined by a diapause hormone that is secreted by the suboesophageal ganglion (SG) of the mother moth during the pupal period to act on her developing ovaries and is responsible for induction of embryonic diapause of the silkworm, Bombyx mori . The termination of diapause requires a low temperature of 5°C for 2–3 months, various stimuli can also artificially prevent or terminate diapause, such as HCl, physical stress, corona discharge, or higher oxygen pressure. Once diapause terminates, the embryos are competent to resume development at 25°C and cells enter the M phase; Cell division in the embryos then resumes.\r\n\r\nHowever, the molecular mechanisms involved in generating, maintaining, and breaking diapause have not yet been fully elucidated. Currently, only a few studies have elaborated the molecular mechanisms that regulated diapause. Most of these studies have mentioned the molecular regulation of diapause in the larvae, pupae and adults stage. However, the molecular regulation of embryonic diapause is still unclear, perhaps because of the difficulty of extracting the limited amounts of RNA from a developmental stage that has relatively little tissue and large amounts of protein and lipid.\r\n\r\nThe gene expression profiles linked with diapause have been researched in many insects, such as northern malt fly species, flesh fly, Chinese oak silkworm, Allonemobius socius, Drosophila melanogaster, Colorado potato beetle. Nevertheless, gene expression profile alone is not sufficient to explain gene functions, and mRNA levels incompletely correlate with the protein levels because of the alternative splicing and dynamics of gene translation. Therefore, with the gradual completion of genome sequencing of insect model organisms, proteomics has become the focal point in entomological research including the diapause mechanism. For a long time, two-dimensional gel electrophoresis (2-DE) combined with mass spectrometry (MS) has been frequently used in insect proteomic research. Another relevant approach in proteomics for large-scale characterization of proteome profiles is shotgun proteomics, which is based on in-gel or gel free digestion of protein mixtures followed by liquid chromatography tandem mass spectrometry (LC-MS/MS) separation, MS detection and database searching, and provides an extremely sensitive and high-throughput approach to determine the proteome components in a complex biological sample. This approach has been widely used in many organisms, such as Bombyx mori , Homo sapiens and Orientia tsutsugamushi .\r\n\r\nFollowing the rapid development of proteomics and bioinformatics approaches, the credibility of results is greatly promoted. In this study, we utilized the shotgun LC–MS/MS approach combined with bioinformatics analysis to illuminate the differences among protein identification profiles of the diapause and non-diapause eggs of the silkworm and to find valuable clues regarding the molecular regulation of diapause at the early stage of embryogenesis in the silkworm.\r\n\r\nMaterials and Methods\r\n\r\nSilkworm Rearing and Sample Collection\r\n\r\nBivoltine silkworm strain 932 was protected at 20°C with half on a short-day photoperiod (8-h-light:16-h-dark for the generation of non-diapause eggs) and half on a long-day photoperiod (16-h-light:8-h-dark for the generation of diapause eggs) separately during the hatching period of silkworm eggs. Larvae were reared on fresh mulberry leaves under an environment of 24 to 25°C and 85% relative humidity. Pupae were kept at 25°C. After female moths emerged, they copulated with males (usually within 5 h) and then laid the eggs at 25°C. Diapause and non-diapause eggs which laid by 30 female moths seperately within 1 h were collected to obtain synchronously developing egg batches. All the diapause and non-diapause eggs were collected 24 h after oviposition and divided into 0.1 g each sample (corresponding to about 200 eggs, three copies). All samples were stored at −80°C until use.\r\n\r\nSample Preparation and SDS-PAGE Separation\r\n\r\nSamples of diapause eggs (D) and non-diapause eggs (ND) were mechanically homogenized on ice for 10 min in 10 µL lyses buffer (comprising 8 M urea, 2 M thiourea, 2% CHAPS, 20 mM Tris–HCl, 30 mM DTT) per 1 mg tissue. The samples were then sonicated in an ice-bath for five circles and each circle contained a 30 s sonication with a 30 s interval. The samples were then centrifuged at 12,000×g at 4°C for 15 min. The supernatants were then collected and were centrifuged again and the resultant supernatants were stored at −20°C for further study. The concentrations of protein samples were determined by the Bradford methods.The samples were boiled for 2 min and centrifuged at 12,000×g for 10 min before they were subjected to SDS-PAGE, using a 5% stacking gel and a 12.5% resolving gel. For each sample, a total of 300 µg protein was separated using SDS-PAGE on three lanes with 100 µg each lane. The gels were stained with Coomassie Brilliant Blue (CBB) R250 (Sigma, USA) after electrophoresis.\r\n\r\nIn-gel Digestion\r\n\r\nEach gel lane was manually cut into 8 bands according to the deepness of Coomassie staining (Fig. 1). The gel bands were sliced into 1×1 mm pieces and subjected to in-gel tryptic digestion which was essentially carried out as described by Shevchenko et al.. Briefly, the gel pieces were rinsed thrice using Milli-Q water (Millipore) and destained twice with 25 mM NH4HCO3 in 50% acetonitrile (ACN, Amersham) at 37°C until the color depigmented completely. The dried gels were incubated with 50 mM Tris[2-carboxyethyl]phosphine (TCEP, Sigma) in 25 mM NH4HCO3 at 56°C for 1 h to reductively cleave the disulfide bonds of proteins and then the resulting sulfhydryl functional groups were alkylated by 100 mM iodoacetamide (IAA, Amersham) in 25 mM NH4HCO3 at room temperature in the dark for 0.5 h. Gel pieces were washed twice with 25 mM ammonium bicarbonate in 50% acetonitrile solution, dehydrated twice with 100% acetonitrile, and dried in a vacuum centrifuge. Subsequently, the proteins were digested with 20 ng/µL modified trypsin (Sigma) for 20 h at 37°C.The resulting peptides were extracted twice from the gel pieces with 5% trifluoroacetic acid (TFA, Fluka, Milwaukee, WI, USA) in 50% ACN solution. The pooled extracts were evaporated in a vacuum centrifuge (Labconco, Kansas, MO), and resuspended with 0.1% methanoic acid (Sigma). After trypsin digested, 8 slices were combined to 4 groups for each sample (slice 1 and 2, slice 3 and 4, slice 5 and 6, slice 7 and 8) prior to the LC separation and MS detection.\r\n\r\nSDS-PAGE patterns of the diapause (D) and non-diapause (ND) eggs protein samples respectively.\r\n\r\nThe samples were separated by 12.5% resolving gel in triplicate. The three batches for both ND and D were pooled into two sets of eight prior to in-gel digest. Then eight slices were combined to four groups for each sample (comprising slices 1 and 2, slices 3 and 4, slices 5 and 6, slices 7 and 8 for both ND and D) prior to the LC separation and MS detection.\r\n\r\nShotgun LC-MS/MS Analysis\r\n\r\nAll digested peptide mixtures were separated by reverse phase (RP) HPLC followed by tandem MS analysis. RP-HPLC was performed on a surveyor LC system (Thermo Finnigan, San Jose, CA). Samples were loaded into a trap column (Zorbax 300SB-C18 peptide traps, 300 µm × 65 mm, Agilent Technologies, Wilmington, DE) first at a 3 µL/min flow rate for peptides enrichment and desalting. After flow-splitting down to about 1.5 µL/min, peptides were transferred to the analytical column (RP-C18, 150 µm × 150 mm, Column Technology, Inc., Fremont, CA) for separation with a 195 min linear gradient from 95% buffer A (0.1% methanoic acid in water) to 50% buffer B (84% ACN, 0.1% methanoic acid in water) at a flow rate of 250 nL/min. The analytical column was regenerated for 20 min with buffer A before loading the next sample. A Finnigan LTQ linear ion trap mass spectrometer equipped with a nanospray source was used for the MS/MS experiment in the positive ion mode. The temperature of the ion transfer capillary was set at 170°C. The spray voltage was 3.0 kV and normalized collision energy was set at 35.0% for MS/MS. The MS analysis was performed with one full MS scan (m/z 400–1800 with a resolution R = 60,000 at m/z 400) followed by 10 MS/MS scans on the 10 most intense ions from the MS spectrum with the dynamic exclusion settings: repeat count 2, repeat duration 30 s, exclusion duration 90 s. Data were acquired in data dependent mode using Xcalibur software (version 2.0.7, Thermo Fisher Scientific).\r\n\r\nDatabase Search\r\n\r\nDatabase search was carried out against the in-house database, with the protein sequences downloaded from NCBInr protein database (http://www.ncbi.nlm.nih.gov/) including the domesticated silkworm (B.mori,2224 entries), the wild silkworm (B. mandarina,13 entries) and the fruit fly (D. melanogaster, 27777 entries), along with predicted B. mori protein sequence data (14,623 entries). _The four raw datasets of ND and D samples were searched against the in-house database on a local server using Turbo SEQUEST software (Bioworks version 3.2, Thermo Finnigan). The mass tolerances of precursor ion and fragmentation ion were set to 1.5 Da and 1.0 Da, respectively. The trypsin-cleavage was at both ends of the protein and two missing cleavage sites were allowed. Only b and y fragment ions were taken into account. Static (carbamidomethyl) modification on cysteine, and variable modifications (oxidation) on methionine were set for all searches. The results for each dataset were stored as.out format files. All the.out files were filtered by Buildsummary software. The protein identification criteria that we used were based on Delta CN (≥0.1) and Xcorr (one charge ≥1.9, two charges ≥2.2, three charges ≥3.75).\r\n\r\nQuantitative RT-PCR\r\n\r\nTotal RNA was extracted from diapause and non-diapause eggs samples using TRIzol reagent (Invitrogen, Carlsbad, CA, USA) according to the manufacturer's instructions and then reverse-transcribed into cDNA with random primer and M-MLV reverse transcriptase (Promega, Madison, WI, USA). The qRT-PCR was performed in 25 µL reactions with 100 ng reverse transcription product, 200 nM each of the forward and reverse primer, and the SYBR® Premix Ex Taq™ (Takara, Tokyo, Japan). The cDNA was amplified in a Rotor-Gene 3000 real-time PCR system (Corbett Research, Sydney, Australia) according to the following program: initial denaturation at 95°C for 10 min, and 40 cycles for amplification at 95°C for 30 s, 56°C for 30 s and 72°C for 30 s followed by an additional steps for melting curve with the increase of temperature from 72 to 95°C at 0.5°C/6 s and 30°C for 30 s. The expression levels were calculated using Rotor-Gene software (version 6.0.19) on the basis of the difference of Ct value by normalizing with the reference gene (BmactinA3, accession No.X04507). The diapause-unique gene BGIBMGA001595 (Aliphatic nitrilase, BmnitA)and the non-diapause-unique genes BGIBMGA002594 (adenylate kinase 2, BmADK2), BGIBMGA011412 (Isocitrate dehydrogenase, BmIDH), BGIBMGA002521 (gamma-glutamyltransferase, BmGGT) were chosen for the qRT-PCR. The qRT-PCR was performed using the following primers: BmnitA forward, TCGGGAAACATCGCAAGAAC and reverse, GCCGTATCTGGTCGCAAATAC, BmADK2 forward, ATCACGCTCAAACGGTTCCTT and reverse, AGACGTCATCAGCGGCTTTC, BmIDH forward, CGTTGCGACCAGACATAAGGA and reverse, TTCACGGATTCGTGTTCCAA, BmGGT forward, GACAGCCTCAAACCCAATCAG and reverse, GCGGCCATAAAGCCATCTC.\r\n\r\nBioinformatics Analysis\r\n\r\nProtein sequences were searched against InterPro member databases using the InterProScan software to identify protein signatures. The compiled RAW outputs were subjected to gene ontology (GO) category analysis using the Web Gene Ontology Annotation Plot (WEGO). The two groups of datasets were simultaneously subjected to online analysis (http://wego.genomics.org.cn/cgi-bin/wego/index.pl) to conveniently compare them in one graph. The P-values were calculated using Pearson Chi-Square test where available. The proteins identified were classified into cellular component, molecular function and biological process. The pathways related to the identified proteins were predicted by searching against KEGG reference pathway database (http://www.genome.jp/kegg/tool/search_pathway.html) with the available protein sequences.\r\n\r\nResults and Discussion\r\n\r\nProteome Profiles of the D Compared with ND Samples\r\n\r\nThe shotgun proteomics strategy, based on proteolytic digestion of complex protein mixtures, peptides LC separation and tandem MS sequencing, has been widely utilized. However, database searching remains the bottleneck for many shotgun proteomics experiments. So in our research, the 4 raw data of each sample sets from LC-MS/MS were subjected to in-house database search using SEQUEST. A total of 501 and 562 proteins were identified from the D and ND respectively (Fig. 2), with 8091 peptides including 1560 unique peptides and 8125 peptides including 1826 unique peptides respectively. The number of common-expressed proteins among the two samples was 309, which was 61.68% and 54.98% of the proteins in the D and ND respectively. Moreover, there were 135 common-expressed proteins with functional annotations (Table S1), the others were predicted proteins from the silkworm genome database. There were 192 and 253 unique proteins identified in the D and ND respectively, of which 24 and 48 had functional annotations (Table 1 and 2).\r\n\r\nVenn diagram shows the numbers of identified proteins from diapause and non-diapause eggs of the silkworms.\r\n\r\nEach number with no overlap of circles shows the number of proteins uniquely observed in that sample, while overlapping circle shows the numbers of identified proteins common to two of the analyzes.\r\n\r\n\nTable 1 the unique expressed proteins of D with functional annotation.\r\n\r\nGene ID/GI number\tTheor.pI/Mw(kDa)\tNo. of peptides\tNo. of unique peptides\tCover percent (%)\tProtein description\tUnique peptide sequences\t \tBGIBMGA013139 gi|153791691|\t7.36/58913.37\t1\t1\t3.28\tAMP-activated proteinkinase\tK.IADFGLSNM*MMDGEFLR.T\t \tBGIBMGA012151 gi|114053021|\t6.62/42674.16\t1\t1\t4.52\tARP1 actin-related protein1-like protein A\tR.AQYVLPDGTQLEIGPAR.F\t \tBGIBMGA007039 gi|112983070|\t6.06/12248.78\t2\t1\t11.43\tBCP inhibitor\tK.FADYTPEEQQSR.L\t \tBGIBMGA006850 gi|112983188|\t6/55968.01\t2\t2\t2.22\tcalcineurin A\tK.NAEHFSHNSVR.G\t \tgi|160358387|\t8.82/59918.81\t1\t1\t2.13\tcytochrome P450 CYP6AE9\tR.EVMEDYVLLDK.I\t \tgi|160358391|\t8.76/59701.51\t1\t1\t2.52\tcytochrome P450 CYP6AE7\tK.TMGYLFFESTTRR.G\t \tBGIBMGA010671 gi|114053117|\t5.38/25082.82\t2\t2\t11.47\teukaryotic translation initiation factor 3 subunit k\tK.NNLANLLGGIDDVTLK.H R.ISGFHDSIR.K\t \tBGIBMGA000619 gi|169234934|\t4.77/23880.86\t2\t2\t11.57\tFK506-binding protein\tK.LTIPASLGYGER.G K.LVEEIFQHEDKDK.N\t \tBGIBMGA004040 gi|112983048|\t8.51/14345.72\t2\t1\t8.87\thypothetical protein LOC692620\tK.LPDLWEELAIK.E\t \tBGIBMGA005301 gi|148298789|\t4.25/28434.35\t2\t2\t13.97\thypothetical protein LOC778506\tK.IDDSVKPTEVAAATEEK.K K.KDDIAPEDSDIAKPETVPEVK.T\t \tgi|114051870|\t6.76/28608.9\t1\t1\t8.33\tlutathione transferase o1\t-.MTYFHSVNAGVIPPPALTDK.L\t \tBGIBMGA014404 gi|112983100|\t4.96/45937.83\t1\t1\t2.38\tmasquerade-like serine proteinase homolog\tR.AGEWDTQNTK.E\t \tBGIBMGA00145 gi|112982906|\t6.53/44709.45\t1\t1\t3.71\tMAP kinse-ERK kinase\tR.IPESILGTITSAVLK.G\t \tBGIBMGA00052 gi|112983290|\t5.8/87764.03\t1\t1\t1.68\tneutral endopeptidase 24.11\tR.VIDLDQASLGLSR.E\t \tgi|182509212|\t8.8/51474.1\t1\t1\t2.24\tolfactory receptor-like receptor\tR.LLLNLLLANR.R\t \tBGIBMGA008164 gi|112983896|\t6.33/50013.91\t2\t2\t6.42\tputative paralytic peptide-binding protein\tK.VGDQQLFLIENR.E R.IPVVSVDDSDTSSYIR.D\t \tBGIBMGA004657 gi|114052160|\t6.44/27143\t1\t1\t4.47\tproteasome subunit alpha type 6-A\tR.GTDAAVVAAQR.K\t \tBGIBMGA000867 gi|112982669|\t10.46/28069.61\t1\t1\t4.21\tribosomal protein S2\tK.TYAYLTPDLWR.D\t \tBGIBMGA010178 gi|112982956|\t11.37/35506.38\t1\t1\t3.59\tsplicing factor arginine/serine-rich 6\tR.VYVGGLPFGVR.E\t \tgi|169234928|\t5.06/369030.66\t1\t1\t0.40\ttitin1\tK.PEEANTLRALVDK.I\t \tBGIBMGA000624 gi|148298772|\t5.24/380464.73\t1\t1\t0.36\ttitin2\tK.PEGPVIM*REISR.E\t \tBGIBMGA011615 gi|112983266|\t8.9/176184.13\t2\t1\t1.03\ttopoisomerase II\tK.FKM*EKLEDDIVLLMSR.R\t \tBGIBMGA001584 gi|160333857|\t4.81/29389.15\t1\t1\t5.16\ttropomyosin isoform 3\tK.LSEASQAADESER.I\t \tgi|112984130|\t6.96/33312.27\t2\t1\t4.93\tyellow6\tR.QNLIPYNFIDDVIR.T\t \t\r\n\r\n\nTable 2 the unique expressed proteins of ND with functional annotation.\r\n\r\nGene ID/GI number\tTheor.pI/Mw(kDa)\tNo. of peptides\tNo. of unique peptides\tCover percent (%)\tProtein description\tUnique peptide sequences\t \tgi|114051770|\t5.99/43014.71\t4\t3\t11.17\t26S proteasome non-ATPase regulatory subunit 13\tK.FAVIDVSDFLTK.K R.ALAGAGAVLTAAR.V R.GDNLIQLYNNFLTTFESK.I\t \tBGIBMGA001490 gi|112983564|\t4.75/44677.53\t3\t3\t8.21\t45 kDa immunophilin FKBP45\tK.ALSGGVQIEDLK.L K.TGDSIAFLTNGK.C K.VVMVYYEGR.L\t \tBGIBMGA012298gi|114053253|\t7.12/52872.82\t2\t1\t3.31\t6-phosphogluconate dehydrogenase\tK.NPQLTNLLLDPFFSER.I\t \tBGIBMGA008749 gi|114052224|\t8.37/19958.11\t1\t1\t6.32\tadenylate cyclase\tR.FYADSAATLVR.Y\t \tBGIBMGA011922 gi|114052488|\t6.47/40033.17\t1\t1\t3.99\talcohol dehydrogenase\tR.IIGVDINPDKFEVAK.K\t \tBGIBMGA007656 gi|148298845|\t11.42/20174.53\t1\t1\t7.95\targinine/serine-rich splicing factor 7\tR.NPPGFAFVEFEDPR.D\t \tBGIBMGA008670 gi|114052262|\t8.31/27604.56\t3\t2\t13.58\tATP synthase\tK.KENVLLQLEAAYR.E K.LAAWLDKEVEATENEWNEGR.N\t \tBGIBMGA000622 gi|168823429|\t5.93/544514.58\t1\t1\t0.23\tBm kettin\tK.TDIIDESVTIK.P\t \tBGIBMGA011218 gi|114052306|\t4.77/84822.51\t4\t3\t5.95\tCarboxylesterase\tK.TETPEQAPSVK.E K.VLETSNNLFGLVIEK.E R.TANIPILVGTTSLEYACER.K\t \tBGIBMGA006916 gi|134254470|\t8.54/61666\t1\t1\t5.36\tcytochrome P450 18a1\tK.LPPGPWGPPVVGYLPFLGVRHKTFLQLAR.N\t \tBGIBMGA003391 gi|114052917|\t9.94/12748.96\t1\t1\t8.40\tDnaJ domain-containing protein\tK.FDSETLANSK.Y\t \tBGIBMGA014398 gi|114052402|\t5.49/49150.93\t1\t1\t2.06\teukaryotic translation termination factor 1\tK.SQEGSQFVR.G\t \tBGIBMGA007889 gi|114051800|\t5.71/36908.18\t2\t1\t8.51\teukaryotic translation initiation factor 3 subunit 2 beta\tK.IHSVKEHTHQITDM*QLSRDGTMFVTASK.D\t \tBGIBMGA002526 gi|112983824|\t5.34/74917.21\t2\t1\t1.54\tecdysteroid-inducible angiotensin-converting enzyme-related gene product\tK.NSIIEKPTDR.E\t \tBGIBMGA008780 gi|112982960|\t6.75/26033.48\t3\t2\t6.11\tferritin\tK.ALASLYLK.R R.IFFIHR.E\t \tBGIBMGA011438 gi|112983348|\t8.72/22373.68\t1\t1\t6.03\tglutathione peroxidase\tR.HGPNTDPLDLVK.S\t \tBGIBMGA005781 gi|169646838|\t5.15/25392.59\t1\t1\t8.07\theat shock protein 25.4\tK.IQNLPWDVNSEGSWVYEK.D\t \tBGIBMGA008736 gi|112982982|\t5.12/44820.95\t3\t3\t7.80\themolin\tK.EAPAEVLFR.E K.SDFGVASTR.A R.TLATQGEDVTIPCK.A\t \tBGIBMGA004112 gi|114051598|\t5.34/69314.08\t1\t1\t2.81\tleukotriene A4 hydrolase\tR.YLEDLIGGPEVFDDFLR.S\t \tgi|114052571|\t6.43/23644.23\t2\t1\t11.82\tlysophospholipase\tR.HTASLIFLHGLGDTGHGWASTIAGIR.G\t \tBGIBMGA006419 gi|153792270|\t6.12/67258.39\t3\t2\t5.39\tmalate dehydrogenase\tK.YCTFNDDIQGTAAVAVAGLLASLR.L R.FAQDHAPVR.T\t \tBGIBMGA004838 gi|114052779|\t10.08/48824.97\t1\t1\t3.29\tmitochondrial ribosomal protein S5\tK.LWKSVTSVSNAGAK.K\t \tBGIBMGA002937 gi|182509194|\t6.51/124059.24\t1\t1\t1.82\tnitric-oxide synthase like protein\tK.QFVSCTVKANKDLGDASAER.S\t \tgi|112984438|\t5.55/203642.64\t7\t6\t4.16\tovarian serine protease\tK.GPYEQIYK.V K.ISLLHLK.S K.KPLFAVSNTEGFLHR.N K.LTIEQIEQEANHICHYLGFSSAR.T K.VYVEGNYR.C R.TLYTDFDETPLFLR.E\t \tBGIBMGA008602 gi|112984330|\t6.26/262354.41\t1\t1\t0.64\tp260\tR.QPYCGLYGLVKKLSK.Q\t \tBGIBMGA001292 gi|114052687|\t5.44/26024.4\t1\t1\t5.75\tprogesterone receptor membrane component 2\t-.MDSPPEVENVQAK.H\t \tBGIBMGA003028 gi|114052034|\t8.29/27871.9\t1\t1\t5.62\tproteasome beta-subunit\tK.LTVEDPVTLEYITR.Y\t \tBGIBMGA013898 gi|114050993|\t4.98/26874.64\t2\t1\t5.76\tproteasome zeta subunit\tK.AIGSGSEGAQQSLK.E\t \tBGIBMGA005994 gi|114052713|\t5.92/37494.59\t1\t1\t3.94\tproteasome 26S non-ATPase subunit 7\tR.DIKDTTVGSLSQR.I\t \tBGIBMGA007169 gi|112983398|\t6.33/21828.97\t1\t1\t5.73\tras oncoprotein\tK.LVVVGAGGVGK.S\t \tBGIBMGA001800 gi|112982800|\t11.66/48110.93\t1\t1\t2.57\tribosomal protein L4\tK.IPELPLVVADK.V\t \tBGIBMGA002811 gi|160333861|\t9.98/25492.78\t1\t1\t13.24\tribosomal protein L10\tR.ATVDDFPLCVHLVSDEYEQLSSEALEAGR.I\t \tBGIBMGA010867 gi|112984078|\t10.27/29623.84\t6\t4\t11.41\tribosomal protein S4\tK.LGGVYAPR.P K.LRECLPLVIFLR.N K.YALTGNEVLK.I R.ECLPLVIFLR.N\t \tgi|112982661|\t10.97/29058.07\t1\t1\t6.72\tribosomal protein S6\tR.KGAQEIPGLTDGNVPRR.L\t \tgi|114052054|\t5.9/19293.48\t1\t1\t7.69\trotamase Pin1\tR.TKEEALDILQEYR.R\t \tBGIBMGA004340 gi|114052803|\t8.6/11739.45\t1\t1\t15.32\tsalivary cysteine-rich peptide\tK.SCSEPAAYGGGTGSSSK.H\t \tBGIBMGA007079 gi|114052783|\t7.56/51220.85\t1\t1\t4.09\tserine hydroxymethyltransferase\tK.LLNSNLWEADPELFDIIVK.E\t \tBGIBMGA010212 gi|114051043|\t5.28/51672.02\t1\t1\t3.06\tserine protease inhibitor serpin\tR.EVRIKFSTIIDSLK.K\t \tBGIBMGA006043 gi|114052122|\t7.74/33484.08\t1\t1\t4.12\tstathmin\tR.LTLEQQTAEVYK.A\t \tBGIBMGA011499 gi|112984508|\t9.03/50600.85\t2\t1\t2.78\ttranslation initiation factor 2 gamma subunit\tR.TTQSNLHQQDLSK.L\t \tBGIBMGA010283 gi|148298831|\t7/88104.06\t1\t1\t1.02\ttuftelin interacting protein 11\tK.M*FPEDILK.E\t \tBGIBMGA005794 gi|151301088|\t7.17/22241.48\t1\t1\t5.10\tubiquinone biosynthesis protein COQ7-like protein\tR.TLM*QDPNVDK.E\t \tgi|114052396|\t6.83/20898.87\t3\t2\t10.38\tubiquitin-conjugating enzyme E2M\tK.DLNELNLPK.T K.EAADELQSNR.R\t \tBGIBMGA005628 gi|114050831|\t7.52/52953.59\t1\t1\t4.86\turidine 5'-monophosphate synthase\tK.KDDTCLIIEDVITSGSSILETVK.D\t \tBGIBMGA013190 gi|114053067|\t5.95/20271.21\t2\t1\t5.59\tvacuolar protein sorting 29\tK.TLASDVHVVR.G\t \tBGIBMGA013964 gi|114050729|\t8.46/44007.51\t1\t1\t4.15\tvacuolar ATPase subunit C\tR.YGLPVNFQAVVM*VPAR.K\t \tBGIBMGA010247 gi|114052088|\t8.98/26119.27\t1\t1\t5.75\tvacuolar ATP synthase subunit E\tR.LELIAQQLLPEIR.N\t \tBGIBMGA004036 gi|112983588|\t5.83/65983.72\t1\t1\t1.83\tVlg protein\tK.VAVAYGGTAVR.H\t \t\r\n\r\nTheoretical Two-dimensional Distribution of the Identified Proteins\r\n\r\nThe theoretical isoelectric point (pI) and molecular weight (Mw) of the identified proteins were calculated using the Compute pI/Mw tool (http://cn.expasy.org/tools/pi_tool.html) according to the predicted amino acid sequences. It showed that 87.30% of the total proteins were distributed in the range of pI 4–7 and 8–10 (Fig. 3a). About 60.87% of the proteins were distributing in the range of 15–60 kDa (Fig. 3b). In the two kinds of eggs, less than 7.62% of proteins showed pI 7–8. Furthermore, 25 and 29 proteins with higher pI (more than 10), usually difficult to be separated by 2-DE, were also identified from D and ND respectively. It also revealed that the protein distributions were nearly identical in the D and ND samples.\r\n\r\nTheoretical pI and Mw distribution of the identified proteins.\r\n\r\n(a) Distribution of pI. (b) Distribution of Mw. The theoretical isoelectric point (pI) and molecular weight (Mw) of the proteins were calculated using the Compute pI/Mw tool (http://cn.expasy.org/tools/pi_tool.html) according to the predicted amino acid sequences.\r\n\r\nProfiling of Common-expressed Proteins in the D and ND Samples\r\n\r\nOn the list of common-expressed proteins of D and ND with annotations, many functional proteins were detected (Table S1). Among these common-expressed proteins we identified, vitellogenin precursor (gi:112983746), egg-specific protein (gi:187281695) and low molecular lipoprotein 30 K precursor (gi:156119320, gi:162461355 and gi:112984502) which comprise the three main storage proteins in the silkworm eggs. The relative concentration of a protein identified by MS in a biological sample is directly related to the number of identified peptides, neglecting possible effects such as the enzymatic digestion constraint, the detection mass range of the mass spectrometer and differential post-translational modifications. Therefore the number of identified unique peptides assembled into a protein may reflect the protein’s relative abundance. So according to the amount of identified unique peptides, vitellogenin precursor, egg-specific protein, low molecular lipoprotein 30 K precursor, and some predicted proteins such as BGIBMGA013342-PA, BGIBMGA004585-PA and BGIBMGA004399-PA are highly represented in the samples. In these predicted proteins, the BGIBMGA013342-PA and BGIBMGA004585-PA were involved in lipid transport (GO: 0006869) with the molecular Function of lipid transporter activity (GO: 0005319), while the BGIBMGA004399-PA was identified to be a low molecular weight lipoprotein which located at the extracellular region (GO: 0005576).These highly-represented proteins were important at the early stage of the embryogenesis.\r\n\r\nMoreover, many heat shock proteins were identified, such as heat shock protein HSP19.9, HSP20.1, HSP20.4, HSP20.8, HSP23.7, HSP70, HSP90 and heat shock cognate protein. These proteins, HSP19.9, HSP20.1, HSP20.4, HSP20.8, HSP23.7 belong to a family of small heat shock protein (sHSP) which mainly function as molecular chaperones to protect proteins from being denatured during extreme conditions, especially under high temperature stress. The sHSP family is functionally more diverse than other HSPs. Moreover sHSPs can also develop a protection function under other stress conditions, such as cold, drought, oxidation, hypertonic stress, UV, and heavy metals, and even under high population densities of organisms.The HSPs including sHSPs also play an important role in normal development. However, insect orthologous sHSPs may not be associated with response to environmental stresses and may be involved in basic metabolic processes. Moreover, silkworm sHSPs may play an important role in the development of the germocyte and have functions in immune defense mechanisms. HSP70, heat shock cognate protein,HSP90 and sHSPs are associated with diapause in a number of species.But in this study, some HSPs were highly expressed both in diapause and non-diapause eggs, including HSP90, HSP20.8, HSP20.4 and HSP70. A possible explanation is that these HSPs may play an important function in the initial embryos development regardless of diapause or non-diapause.\r\n\r\nFurthermore, among the common-expressed proteins with functional annotations, 37 metabolic enzymes were identified (11.97% of the common-expressed proteins with functional annotations). The metabolic enzymes include ADP/ATP translocase, fibroinase, glucose-6-phosphate isomerase, glutamate dehydrogenase, isocitrate dehydrogenase, salivary secreted ribonuclease, transketolase etc. The identification of these essential metabolic enzymes indicated that during embryogenesis stage metabolism was active in the eggs. Many ribosomal proteins (L7, L7A, L9, L10A, L11, L13, L13A, L15, L17, L18, L19, L23A, L38, S3, S5, S7, S8, S9, S15A, S27, S28, and S30) were found in D and ND samples. This abundance of ribosomal proteins on the first day after oviposition might be closely related to protein synthesis during embryogenesis.\r\n\r\nIn addition, 62 proteins that may also play important roles in the egg development were identified commonly, such as 14-3-3 epsilon protein, antennal binding protein, calreticulin, cellular retinoic acid binding protein, chemosensory protein 11, exuperantia, perilipin, profiling, transferring and etc. 14-3-3 epsilon protein (Bm14-3-3ε) is a member of the 14-3-3 protein family. The 14-3-3 proteins along with partner proteins(for example CDK11, PFTK1, GSK3beta, Chk1) are involved in the regulation of several crucial cellular processes including metabolism, signal transduction, cell development, differentiation, apoptosis, transcription, stress responses and malignant transformation. In this study, we identified Bm14-3-3ε in the D and ND samples. This result indicated that Bm14-3-3ε may play a role in regulating embryonic development of silkworm. Another important protein was translationally controlled tumor protein (TCTP), a highly conserved protein upregulated in various tumours. Hsu et al reported that reducing Drosophila TCTP (dTCTP) levels will reduces its cell size, cell number and organ size. Further more, calreticulin located in the endoplasmic reticulum, has been implicated in many diverse functions, including: regulation of intracellular Ca2+ homeostasis, chaperone activity, steroid-mediated gene regulation, and cell adhesion. The highest level of mRNA expression of calreticulin was exhibited in the fat body of Bombyx mori . Our result showed that calreticulin was identified during the initial embryogenesis both in D and ND samples, which means that calreticulin is essential during embryonic development.\r\n\r\nProfiling of D-Unique Proteins\r\n\r\nThe characteristic proteins with functional annotations unique to the D samples are shown in Table 1. Among the D-unique proteins, BCP inhibitor (BCP1, Bombyx cysteine proteinase inhibitor) is inhibitory towards the processing of the enzymatically inactive proform of BCP (pro-BCP) to the activated mature BCP but has no effects on trypsin and pepsin, and is highly expressed in the metamorphosis stage of B.mori . In our study, BCPI which was exclusively identified in the D samples may be involved in regulating the development of diapause-destined embryos to arrest during the G2 cell cycle stage.\r\n\r\nCalcineurin A and FK506-binding protein were other proteins which were uniquely identified in the D sample. Calcineurin A with four EF-hand type calcium-binding structures is localized in the cytoplasm of the pheromone-producing cells and participates in the intracellular signal transduction of PBAN (Pheromone biosynthesis activating neuropeptide) in B. mori . The FK506-binding protein (FKBP) belongs to a ubiquitous family of proteins which participates in a variety of pathways, including protein folding, down-regulation of T-cell activation and inhibition of cell-cycle progression. Both Calcineurin A and FK506-binding protein identified in the diapause eggs are related with the Ca2+ release, suggesting that in the diapause eggs the regulation of the Ca2+ might be more active than in the non-diapause eggs. Other proteins identified in the diapause eggs were cytochrome P450 CYP6AE9, cytochrome P450 CYP6AE7, titin1, titin2, topoisomerase II etc (Table 1).\r\n\r\nIn the diapause eggs, the cephalic lobe and telson and sequential segmentation of the mesoderm are formed before the embryos enters into diapause. In another words, the development of the diapause eggs are on one hand moving toward embryonic development, while on the other hand preparing to enter diapause after oviposition. Thus, these D-unique proteins may be adapted to this physiological need.\r\n\r\nProfiling of ND-Unique Proteins\r\n\r\nCompared with the unique proteins of the D samples, the characteristic proteins unique to the ND samples are shown in Table 2. There were more metabolic enzymes identified in the ND samples than in the D samples, including 6-phosphogluconate dehydrogenase, adenylate cyclase, alcohol dehydrogenase, ATP synthase(subunit B), carboxylesterase, glutathione peroxidase, malate dehydrogenase, ovarian serine protease, serine hydroxylmethyltransferase etc (Table 2). In the non-diapause eggs, the development of the embryos was consistent and without interruption, so the metabolism of the non-diapause eggs was potentially more active than the diapause eggs.\r\n\r\nIn addition, we also identified eukaryotic translation initiation factor 3 subunit 2 beta, translation initiation factor 2 gamma subunit, eukaryotic translation termination factor 1 and more ribosomal proteins (L4, L10, S4, S6) in the non-diapause eggs. These components are important in regulating translation and protein synthesis, indicating that translation of the mRNA and synthesis of proteins was potentially more active in the non-dipause eggs than that in the diapause eggs. Proteasome beta-subunit, proteasome zeta subunit and proteasome 26S non-ATPase subunit 7 were also identified only in the non-dipause eggs, which play important roles in degrading cytosolic and nuclear proteins previously labeled with ubiquitin molecules. Other proteins identified in the non-diapause eggs were ecdysteroid-inducible angiotensin-converting enzyme-related gene product whose transcription was directly induced by 20-hydroxyecdysone (20E), ubiquitin-conjugating enzyme E2M, stathmin and a ras oncoprotein, all of which are associated with embryonic development.\r\n\r\nQuantitative RT-PCR Analysis\r\n\r\nThe dipause-unique gene BmnitA and the non-diapause-unique genes BmADK2, BmIDH and BmGGT were selected to detect their relative mRNA levels by qRT-PCR analysis. Relative gene expression levels were showed in Fig. 4. In the diapause eggs, BmnitA, BmADK2 and BmGGT mRNA relative expression levels were significantly increased compared with those in non-diapause eggs according to the development stage, while BmIDH mRNA relative expression levels were remained constant compared with those in non-diapause eggs according to the development stage, which indicated that BmIDH is a possible non-diapause-unique gene. Interestingly at the initial embryogenesis, the mRNA relative expressions of these four genes were all significantly higher in the non-diapause eggs than in the diapause eggs. In the non-diapause eggs, BmIDH mRNA relative expression was significantly higher more than 21-fold than that in diapause eggs, BmnitA mRNA relative expression was significantly higher more than 21-fold than that in diapause eggs, BmADK2 mRNA relative expression was significantly higher more than 1795-fold than that in diapause eggs, BmGGT mRNA relative expression was significantly higher more than 57-fold than that in diapause eggs. However, the mRNA levels were not always consistent with their protein expressions. In the mRNA level, the expression of these four genes can be detected whether in diapause or non-dipause eggs.\r\n\r\nGene expression level differences of the identified proteins in D and ND eggs of 1, 3, 5, 7 days after oviposition.\r\n\r\nThe mRNA levels of the target genes were normalized by subtracting Ct value of the reference gene, BmactinA3 (No.X04507). D means diapause and N means non-diapause.The (2-ΔΔCT values represents relative gene expression level. Values are means ± SD (standard deviations).\r\n\r\nGene Ontology Analysis of the Functional Categories\r\n\r\nGene ontology (GO) is now widely used to describe the function of genes and gene products in a standardized format (http://www.geneontology.org). To understand the functions of the proteins we identified, the protein sequences were queried against the InterPro databases and the resultant proteins were functionally categorized based on universal GO annotation terms using the online GO tool WEGO. GO signatures of 486 out of the 501 and 545 out of the 562 proteins identified in the diapause and non-diapause eggs respectively were available. They were classified into Cellular Component, Molecular Function and Biological Process according to the GO hierarchy using WEGO (Fig. 5).\r\n\r\nGene ontology categories for the identified proteins of D and ND.\r\n\r\nThe identified proteins were classified into cellular component, molecular function and biological process by WEGO according to the GO terms. The number of genes is the number of times the GO term is used to annotate genes in the cluster. The left-hand shows its proportion in total genes of related genes with GO terms.\r\n\r\nIn the Cellular Component category, proteins mapping to cell, cell part, macromolecular complex and organelle related GO terms were the most abundant, mapping to membrane-enclosed lumen were the fewest. In the subcategory of cell part, 114 and 131 proteins of D and ND respectively were ascribed to intracellular. In the subcategory of organelle, 30 and 35 proteins were separately assigned to membrane-bounded organelle, and 39 and 48 proteins were separately assigned to non-membrane-bounded organelle.\r\n\r\nAccording to the Molecular Function category, most proteins were addressed to binding (194/231 in the D and ND samples respectively) and catalytic activity (177/215 in the D and ND samples respectively), especially the nucleotide, nucleoside, chromatin, ion and nucleic acid binding and hydrolase, oxidoreductase, transferase activities. The groups with much fewer terms (only one protein annotated) include small protein activating enzyme activity, chromatin binding, odorant binding, metal cluster binding, enzyme activator activity, phosphatase regulator activity, transcription repressor activity, transcription initiation factor activity proteins. Moreover, one protein with deaminase activity, two proteins with cyclase activity and one protein with thioredoxin-disulfide reductase activity were annotated functional proteins only in the non-dipause eggs.\r\n\r\nConsidering the Biological Process category, most of proteins were involved in the metabolic and cellular process. Much fewer proteins were involved in the reproduction, reproductive, biological adhesion and developmental process. In the metabolic process subcategories, most proteins were related to primary metabolic, cellular metabolic and macromolecule metabolic process. In the cellular process subcategories, most proteins were associated with the cellular metabolic process followed by the cell communication process. Besides, only one protein was functionally annotated to each of secondary metabolic process, cellular component disassembly, translational initiation, membrane docking and secretion by cell in the diapause eggs.\r\n\r\nGO analysis on the identified proteins presented an overall view on the functional categories of the diapause and non-diapause egg proteomes. In addition, proteins with GO annotation in the diapause eggs were not significantly different (p\u003c0.05) to those in the non-diapause eggs, indicating that the expressed proteins during the early stages of embryonic development have many functional similarities between diapause and non-diapause eggs.\r\n\r\nKEGG Pathway Analysis\r\n\r\nKEGG is a large resource contains information for various model organisms about Molecular Interactions, Reaction Networks, Cellular Processes and Human Diseases. In the present study, a total of 689 of 754 proteins were subjected to query against the KEGG reference pathway database and 192 non-redundant pathways were indicated (Table S2) that most of them were related to Metabolism and Organism Systems. Other pathways such as the Genetic Information Processing, Environmental Information Processing, Cellular Processes and Human Diseases were also detected. 42 and 29 KEGG pathways were involve only in the diapause eggs unique proteins and non-diapause eggs unique proteins, respectively (Table 3, 4).\r\n\r\n\nTable 3 KEGG pathways of diapause unique proteins.\r\n\r\nAnnotation\tTotal Genes\t \tMetabolism\t\t \tLipid Metabolism; Steroid hormone biosynthesis\t3\t \tLipid Metabolism; Linoleic acid metabolism\t3\t \tXenobiotics Biodegradation and Metabolism; Caprolactam degradation\t20\t \tLipid Metabolism; Biosynthesis of unsaturated fatty acids\t14\t \tEnvironmental Information Processing\t\t \tSignal Transduction; ErbB signaling pathway\t20\t \tSignal Transduction; Phosphatidylinositol signaling system\t20\t \tSignal Transduction; mTOR signaling pathway\t20\t \tSignal Transduction; VEGF signaling pathway\t40\t \tCellular Processes\t\t \tCell Growth and Death; Cell cycle - Caulobacter\t20\t \tTransport and Catabolism; Regulation of autophagy\t20\t \tCell Communication; Focal adhesion\t21\t \tCell Motility; Regulation of actin cytoskeleton\t21\t \tGenetic Information Processing\t\t \tFolding, Sorting and Degradation; Sulfur relay system\t20\t \tOrganismal Systems\t\t \tDevelopment; Dorso-ventral axis formation\t20\t \tDevelopment; Axon guidance\t21\t \tDevelopment; Osteoclast differentiation\t40\t \tImmune System; Complement and coagulation cascades\t9\t \tImmune System; Toll-like receptor signaling pathway\t20\t \tEnvironmental Adaptation; Plant-pathogen interaction\t20\t \tImmune System; Natural killer cell mediated cytotoxicity\t40\t \tImmune System; T cell receptor signaling pathway\t40\t \tImmune System; B cell receptor signaling pathway\t40\t \tImmune System; Fc epsilon RI signaling pathway\t20\t \tImmune System; Leukocyte transendothelial migration\t1\t \tEndocrine System; Adipocytokine signaling pathway\t39\t \tHuman Diseases\t\t \tPathogenic Escherichia coli infection\t1\t \tInfectious Diseases; Shigellosis\t1\t \tInfectious Diseases; Hepatitis C\t19\t \tInfectious Diseases; Measles\t19\t \tCancers; Pathways in cancer\t20\t \tCancers; Colorectal cancer\t20\t \tCancers; Renal cell carcinoma\t20\t \tCancers; Pancreatic cancer\t20\t \tCancers; Endometrial cancer\t20\t \tCancers; Glioma\t20\t \tCancers; Prostate cancer\t20\t \tCancers; Thyroid cancer\t20\t \tCancers; Melanoma\t20\t \tCancers; Bladder cancer\t20\t \tCancers; Chronic myeloid leukemia\t20\t \tCancers; Acute myeloid leukemia\t20\t \tCancers; Non-small cell lung cancer\t20\t \t\r\n\r\nKEGG pathways of non-diapause unique proteins.\r\n\r\nAnnotation\tTotal Genes\t \tMetabolism\t\t \tLipid Metabolism; Fatty acid biosynthesis\t18\t \tMetabolism of Cofactors and Vitamins; Ubiquinone and other terpenoid-quinone biosynthesis\t19\t \tAmino Acid Metabolism; Phenylalanine metabolism\t6\t \tMetabolism of Other Amino Acids; Selenocompound metabolism\t40\t \tMetabolism of Other Amino Acids; Cyanoamino acid metabolism\t35\t \tLipid Metabolism; Arachidonic acid metabolism\t51\t \tCarbohydrate Metabolism; Glyoxylate and dicarboxylate metabolism\t20\t \tMetabolism of Terpenoids and Polyketides; Terpenoid backbone biosynthesis\t20\t \tBiosynthesis of Other Secondary Metabolites; Indole alkaloid biosynthesis\t6\t \tBiosynthesis of Other Secondary Metabolites; Isoquinoline alkaloid biosynthesis\t6\t \tBiosynthesis of Other Secondary Metabolites; Betalain biosynthesis\t6\t \tMetabolism of Terpenoids and Polyketides; Insect hormone biosynthesis\t37\t \tGenetic Information Processing\t\t \tTranslation; RNA transport\t20\t \tTranslation; mRNA surveillance pathway\t20\t \tTranscription; Spliceosome\t20\t \tEnvironmental Information Processing\t\t \tSignaling Molecules and Interaction; Cytokine-cytokine receptor interaction\t2\t \tSignal Transduction; Hedgehog signaling pathway\t20\t \tCellular Processes\t\t \tCell Growth and Death; Meiosis - yeast\t20\t \tTransport and Catabolism; Endocytosis\t19\t \tOrganismal Systems\t\t \tImmune System; RIG-I-like receptor signaling pathway\t20\t \tSensory System; Olfactory transduction\t20\t \tSensory System; Taste transduction\t20\t \tSensory System; Phototransduction - fly\t20\t \tExcretory System; Aldosterone-regulated sodium reabsorption\t19\t \tExcretory System; Endocrine and other factor-regulated calcium reabsorption\t58\t \tDigestive System; Carbohydrate digestion and absorption\t19\t \tHuman Diseases\t\t \tInfectious Diseases; Bacterial invasion of epithelial cells\t19\t \tInfectious Diseases; Toxoplasmosis\t1\t \tInfectious Diseases; Amoebiasis\t20\t \t\r\n\r\nThe metabolism pathways in non-diapause eggs were more highly represented than which in diapause eggs. The amino acid metabolism, carbohydrate metabolism and biosynthesis of other secondary metabolites was more highly represented in the non-diapause eggs, meanwhile the lipid metabolism pathway such as steroid hormone biosynthesis, linoleic acid metabolism and biosynthesis of unsaturated fatty acids was more highly represented in the diapause eggs. The translation and transcription pathway such as RNA transport, mRNA surveillance pathway and spliceosome were only detected in the non-diapause eggs. The sensory system (olfactory transduction, taste transduction and phototransduction), excretory system and digestive system were also detected in the non-diapause eggs, which indicated that the embryonic development was in progress.\r\n\r\nErbB signaling pathway, phosphatidylinositol signaling system, mTOR signaling pathway and VEGF signaling pathway which belonged to signal transduction pathway were more active in the diapause eggs. ErbB signaling regulates diverse cellular functions, such as proliferation, migration, differentiation and survival/death, and participates in various developmental processes during both invertebrate and vertebrate early embryogenesis. Interestingly, development and immune system pathways were also more represented in the diapause eggs, maybe the trigger of immune system pathway was the self-protection of diapause eggs during the long diapause stage.\r\n\r\nConclusions\r\n\r\nThis study provides a catalogue and an initial analysis of the proteomes of the diapause and non-diapause eggs during embryogenesis of the silkworm by shotgun proteomic analysis. Unique proteins identified in the two kinds of eggs and common proteins shared by each other were identified and analyzed. The relative mRNA levels of four identified proteins in the two kinds of eggs were also compared using qRT-PCR and showed some inconsistencies with protein expression. GO analysis of these proteins also provided us with a global view of their functions. In addition, KEGG pathway analysis showed the Metabolism, Translation and Transcription pathway were potentially more highly represented in the non-dipause eggs at this stage. These results will also help further research on finding the diapause-associated proteins during egg development of silkworm. However, the shotgun proteomic analysis has some shortcomings, for example, it is too complex for protein assembly and it mainly depends on bioinformatic methods. With the development of genomics and bioinformatics, the shotgun LC-MS/MS will be a promising strategy in proteomics research.\r\n\r\nSupporting Information\r\n\r\n\r\n","tracks":[]}