# Proteomic Analysis of Larval Midgut from the Silkworm (Bombyx mori)

Authors: Sai Zhang; Yunmin Xu; Qiang Fu; Ling Jia; Zhonghuai Xiang; Ningjia He
Journal: Comparative and Functional Genomics
 (2011)

---

## Abstract

The midgut is the major organ for food digestion, nutrient absorption and also a barrier for foreign substance. The 5th-instar larval stage of silkworm is very important for larval growth, development, and silk production. In the present study, we used 2-DE and matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF-MS) to analyze the midgut proteins from the 5th-instar larvae as well as the midgut proteins under starvation condition. A total of 96 proteins were identified in this study; and among them, 69 proteins were observed in midgut for the first time. We also found that the silkworm larval midgut responded to starvation by producing a 10 kDa heat shock protein and a diapause hormone precursor.


---

## Body

## 1. Introduction

Silkworm,Bombyx mori, is a domesticated insect for silk production. It is a lab-reared animal, and its relatively large size allows silkworm to be an important model in molecular genetics and in structural and functional genomics [1]. In addition, the silkworm also is a phytophagous insect and therefore a representative of lepidopteran pest insects.The silkworm, as a holometabolous insect, has four remarkably different developmental stages: egg, larva, pupa, and moth. The larva is the only stage in which the animals feed and excrete in the whole life cycle. This stage, especially the 5th instar, is very important because the larvae have enough nutrient for growth, development, and silk production. The larval midgut is formed of an epithelial cell monolayer composed of columnar, goblet, and stem cells [2]. The columnar cells are mainly responsible for food digestion and nutrition absorption. Digestion is usually controlled by digestive enzymes and is dependent on their localization in gut [3]. The apical membrane of columnar cells is characterized by a well-developed brush border [4]. Goblet cells participate in the ionic regulation in gut, and the regenerative cells are responsible for renewing epithelial cells during metamorphosis [5].Midgut is also a barrier to the foreign substances during food digestion. It has been found that some proteins such as lipase [6] and SP-2 [7] in midgut have antiviral activity against Bombyx mori nuclear polyhedrosis virus (BmNPV). Moreover, midgut has early been recognized as one of the important targets for insect control. One successful example is the transgenic crops that produce B.  thuringiensis crystal δ-endotoxins. These toxins bind to their receptors and then form a prepore oligomeric structure through which cell content leaking leads to the death of insect [8]. An alternative method for more specific control of insect is to silence the expression of genes using RNA interference (RNAi) through the midgut [9].In order to understand the molecular mechanisms for nutrition digestion and midgut-derived defense, Yao et al. [10] and Kajiwara et al. [11] separately reported the proteomic analysis of the silkworm larval midgut. The former analysis was only performed on the 5th-instar day-3 silkworm larvae, and the later work focused on the midgut proteins from 5th-instar day-2 larvae. It is well known that biochemical metabolism is dramatically changed during the development of larval midgut in the 5th instar. In our current study, we are interested in ascertaining which proteins participate in midgut development in the penultimate instar larvae. Furthermore, when starvation is acted on animals, the animals will take violent responses to starvation by triggering emergent change in protein metabolism [12, 13]. These also can be revealed by proteomic analysis.In the present study, we used 2-DE combined with MALDI-TOF-MS to analyze the midgut proteins of the 5th-instar larvae under normal and starvation conditions. As a result, 96 proteins were identified in this study and among them 69 proteins were observed for the first time. The results also indicated that the starved silkworm larval produced a 10 kDa heat shock protein and a diapause hormone precursor in response to this disadvantage condition.


## 2. Materials and Methods

### 2.1. Insect and Protein Preparation

The silkworm strainp50 (DaZao) used in this study was provided by the Institute of Sericulture and System Biology at the Southwest University of China. The larvae were reared on the fresh mulberry leaves at 25°C to the fifth instar day 3. Then, the larvae were randomly divided into two groups. The larvae in group I still were reared on fresh mulberry leaves, while those in group II were starved until wandering stage (about 240 h). Larvae randomly sampled from each group were dissected at 24-hour intervals, and the midguts were washed with 0.75% NaCl solution and stored at −20°C. Larval midguts were ground to powder in liquid nitrogen with mortar and pestle. The proteins from silkworm midguts were extracted with 1 mL lysis buffer (8 M urea, 4% (m/v) 3[(3-cholamidopropyl) dimethylammonio]-1-propane sulphonate (CHAPS, Sangon, shanghai, China), and 4% (m/v) dithiothreitol (DDT, BBI, Canada)) at 4°C for 40 min. After centrifuging twice at 12,000 g for 15 min, the supernatants were collected for proteomics analysis. The concentration of proteins was measured by Bradford kit as described in manufacturer’s instructions (TIANGEN, Beijing, China).


### 2.2. Two-Dimensional Gel Electrophoresis

The midgut proteins were separated by two-dimensional electrophoresis according to the manufacturer’s instructions (GE Healthcare, USA). Briefly, proteins were dissolved in rehydration solution (8 M urea, 2% CHAPS, 0.28% DDT, and 1% IPG buffer pH 3–10). The protein solution was loaded onto an immobilized PH gradients strip (IPG strip, 18 cm, line, pH 3–10, GE Healthcare, Sweden). Isoelectric focusing (IEF) was performed as follows: 50 V for 12 h, 100 V for 1 h, 200 V for 1 h, 500 V for 30 min, 1000 V for 30 min, 3000 V for 30 min, and 8000 V for a total of 70000 Vhr (voltage multiplied by hours). The current did not exceed 50 mA per strip. After IEF, the strips were equilibrated in equilibration buffer I (50 mM Tris-HCl, pH 8.8, 6 mol/L urea, 30% glycerol, 2% SDS, and 1% DTT) for 15 min and then in buffer II (replaced 1% DTT with 2.5% iodacetamide) for another 15 min. The equilibrated gel strip was loaded onto a 12.5% SDS-polyacrylamide gel. Electrophoresis was carried out in an Ettan DALTsix vertical electrophoresis system (Amersham Biosciences, Sweden) at 10 mA per gel for approximately 1 h, and then switched to 40 mA per gel until the bromophenol blue reached the bottom of the gel. Protein was visualized by silver staining.


### 2.3. Gel Image Analysis

The resultant 2D gels were scanned by Image Scanner III (GE Healthcare, USA) with LabScan 6.0 software (GE Healthcare, USA) at an optical resolution of 300 dpi. The scanned gels were analyzed using the ImageMaster 2D platinum 6.0 (GE Healthcare, USA) software supplied by the manufacturer, and the parameters for spots detecting were set to smooth: 2–5, minarea: 30–50, and saliency: 5–15.


### 2.4. In-Gel Digestion and MALDI-TOF-MS Analysis

The protein spots excised from the gels were cut into small pieces (about 1 mm3) and transferred into 0.5 mL Eppendorf tubes. In-gel digestion was carried out according to what was previously described in [14]. In brief, 100 μL of destaining solution (15 mmol/L K3[Fe(CN)6], 50 mmol/L Na2S2O3) was added to each tube for 10 min destaining. Then, the gel pieces were washed with ddH2O for 15 min for three times. The gel pieces were dehydrated with 100 μL ACN until the gel pieces turned opaque white. The ACN was discarded after centrifugation. The gel pieces were dried at room temperature for more than 40 min. Fifty microgram trypsin (10 mg/mL in 25 mM NH4HCO3, PH 8.0, Sigma, USA) was added to each tube and the incubation was carried out at 37°C for 16–20 h. The tryptic peptides were extracted twice with 30 μL 50% ACN/5% trifluoroacetic acid (TFA, Sigma, USA). The supernatants were pooled and dried in a Centrivap Cold Trap (LABCONCO, USA). Prior to mass spectrum analysis, the peptides were redissolved in 3 μL of 50% ACN/0.5% TFA and equally mixed with saturated α-cyano-4-hydroxycinamic acid (CHCA, Sigma, USA). The mass spectra of the peptides were recorded on a Voyager DE PRO MALDI-TOF-MS (Applied Biosystems, USA). The database used in this analysis contained two parts of protein sequence: 6642 silkworm protein sequences downloaded from NCBI (http://www.ncbi.nlm.nih.gov/) with keyword of Bombyx and 14623 sequences generated by the annotation of silkworm genomic sequences (http://silkworm.swu.edu.cn/silkdb/doc/download.html). The analysis of mass spectra and the process of protein identification were carried out according to a previous study [14].


### 2.5. RNA Extraction and Real-Time PCR

Total RNA was extracted from the midgut tissues in two groups using Trizol reagent (Invitrogen, USA). The RNA samples were dissolved in RNase-free water, and the concentrations were measured using spectrophotometer (Amersham Biosciences, Sweden). According to the protocol, 3μg of total RNAs were used to synthesize the first-strand cDNA in 25 μL of reverse transcription PCR system with AMV Reverse Transcriptase (Promega, USA). For real-time PCR, the primers designed for 10 kDa small heat shock protein gene were 5′-GTTCCTCTTCTGGACCGTG-3′ (forward) and 5′-GACCGACCGCTACTACTTCT-3′ (reverse), and for diapause hormone precursor gene were 5′-GGAAATGTACCAACCTGACCC-3′  (forward) and 5′-AACGCTTCTGGCAACCCTAT-3′ (reverse). Real-time PCR was performed in triplicate for each gene of interest in 20 μL reactions with reverse transcription product, 2 × PCR buffer, 50 × ROX Reference Dye, 100 nM of each of the forward and reverse primer, and RNase free ddH2O. The PCR reaction was performed on an ABI Prism 7000 Sequence Detection System (Applied Biosystems) using the following program: initial denaturation at 95°C for 10 sec and 40 cycles of 95°C for 5 sec 60°C for 31 sec. β-actin3 was used as inner control.


## 2.1. Insect and Protein Preparation

The silkworm strainp50 (DaZao) used in this study was provided by the Institute of Sericulture and System Biology at the Southwest University of China. The larvae were reared on the fresh mulberry leaves at 25°C to the fifth instar day 3. Then, the larvae were randomly divided into two groups. The larvae in group I still were reared on fresh mulberry leaves, while those in group II were starved until wandering stage (about 240 h). Larvae randomly sampled from each group were dissected at 24-hour intervals, and the midguts were washed with 0.75% NaCl solution and stored at −20°C. Larval midguts were ground to powder in liquid nitrogen with mortar and pestle. The proteins from silkworm midguts were extracted with 1 mL lysis buffer (8 M urea, 4% (m/v) 3[(3-cholamidopropyl) dimethylammonio]-1-propane sulphonate (CHAPS, Sangon, shanghai, China), and 4% (m/v) dithiothreitol (DDT, BBI, Canada)) at 4°C for 40 min. After centrifuging twice at 12,000 g for 15 min, the supernatants were collected for proteomics analysis. The concentration of proteins was measured by Bradford kit as described in manufacturer’s instructions (TIANGEN, Beijing, China).


## 2.2. Two-Dimensional Gel Electrophoresis

The midgut proteins were separated by two-dimensional electrophoresis according to the manufacturer’s instructions (GE Healthcare, USA). Briefly, proteins were dissolved in rehydration solution (8 M urea, 2% CHAPS, 0.28% DDT, and 1% IPG buffer pH 3–10). The protein solution was loaded onto an immobilized PH gradients strip (IPG strip, 18 cm, line, pH 3–10, GE Healthcare, Sweden). Isoelectric focusing (IEF) was performed as follows: 50 V for 12 h, 100 V for 1 h, 200 V for 1 h, 500 V for 30 min, 1000 V for 30 min, 3000 V for 30 min, and 8000 V for a total of 70000 Vhr (voltage multiplied by hours). The current did not exceed 50 mA per strip. After IEF, the strips were equilibrated in equilibration buffer I (50 mM Tris-HCl, pH 8.8, 6 mol/L urea, 30% glycerol, 2% SDS, and 1% DTT) for 15 min and then in buffer II (replaced 1% DTT with 2.5% iodacetamide) for another 15 min. The equilibrated gel strip was loaded onto a 12.5% SDS-polyacrylamide gel. Electrophoresis was carried out in an Ettan DALTsix vertical electrophoresis system (Amersham Biosciences, Sweden) at 10 mA per gel for approximately 1 h, and then switched to 40 mA per gel until the bromophenol blue reached the bottom of the gel. Protein was visualized by silver staining.


## 2.3. Gel Image Analysis

The resultant 2D gels were scanned by Image Scanner III (GE Healthcare, USA) with LabScan 6.0 software (GE Healthcare, USA) at an optical resolution of 300 dpi. The scanned gels were analyzed using the ImageMaster 2D platinum 6.0 (GE Healthcare, USA) software supplied by the manufacturer, and the parameters for spots detecting were set to smooth: 2–5, minarea: 30–50, and saliency: 5–15.


## 2.4. In-Gel Digestion and MALDI-TOF-MS Analysis

The protein spots excised from the gels were cut into small pieces (about 1 mm3) and transferred into 0.5 mL Eppendorf tubes. In-gel digestion was carried out according to what was previously described in [14]. In brief, 100 μL of destaining solution (15 mmol/L K3[Fe(CN)6], 50 mmol/L Na2S2O3) was added to each tube for 10 min destaining. Then, the gel pieces were washed with ddH2O for 15 min for three times. The gel pieces were dehydrated with 100 μL ACN until the gel pieces turned opaque white. The ACN was discarded after centrifugation. The gel pieces were dried at room temperature for more than 40 min. Fifty microgram trypsin (10 mg/mL in 25 mM NH4HCO3, PH 8.0, Sigma, USA) was added to each tube and the incubation was carried out at 37°C for 16–20 h. The tryptic peptides were extracted twice with 30 μL 50% ACN/5% trifluoroacetic acid (TFA, Sigma, USA). The supernatants were pooled and dried in a Centrivap Cold Trap (LABCONCO, USA). Prior to mass spectrum analysis, the peptides were redissolved in 3 μL of 50% ACN/0.5% TFA and equally mixed with saturated α-cyano-4-hydroxycinamic acid (CHCA, Sigma, USA). The mass spectra of the peptides were recorded on a Voyager DE PRO MALDI-TOF-MS (Applied Biosystems, USA). The database used in this analysis contained two parts of protein sequence: 6642 silkworm protein sequences downloaded from NCBI (http://www.ncbi.nlm.nih.gov/) with keyword of Bombyx and 14623 sequences generated by the annotation of silkworm genomic sequences (http://silkworm.swu.edu.cn/silkdb/doc/download.html). The analysis of mass spectra and the process of protein identification were carried out according to a previous study [14].


## 2.5. RNA Extraction and Real-Time PCR

Total RNA was extracted from the midgut tissues in two groups using Trizol reagent (Invitrogen, USA). The RNA samples were dissolved in RNase-free water, and the concentrations were measured using spectrophotometer (Amersham Biosciences, Sweden). According to the protocol, 3μg of total RNAs were used to synthesize the first-strand cDNA in 25 μL of reverse transcription PCR system with AMV Reverse Transcriptase (Promega, USA). For real-time PCR, the primers designed for 10 kDa small heat shock protein gene were 5′-GTTCCTCTTCTGGACCGTG-3′ (forward) and 5′-GACCGACCGCTACTACTTCT-3′ (reverse), and for diapause hormone precursor gene were 5′-GGAAATGTACCAACCTGACCC-3′  (forward) and 5′-AACGCTTCTGGCAACCCTAT-3′ (reverse). Real-time PCR was performed in triplicate for each gene of interest in 20 μL reactions with reverse transcription product, 2 × PCR buffer, 50 × ROX Reference Dye, 100 nM of each of the forward and reverse primer, and RNase free ddH2O. The PCR reaction was performed on an ABI Prism 7000 Sequence Detection System (Applied Biosystems) using the following program: initial denaturation at 95°C for 10 sec and 40 cycles of 95°C for 5 sec 60°C for 31 sec. β-actin3 was used as inner control.


## 3. Results

### 3.1.  2-DE Profiles of Silkworm Larval Midgut Proteins

The silkworm midgut proteins were extracted from the fifth-instar day-4 (V4) to day-8 (V8) feeding larvae and from the fifth-instar day-4 (V4N) to day-13 (V13N) nonfeeding larvae. A total of 15 extractions were separated by 2D technology, as shown in Supplemental Data 1 and 2 in Supplemental Material available at doi:10.1155/2011/876064. The protein spots were then detected using ImageMaster 2D platinum 6.0 software. On the gels V4 to V8, 439, 461, 498, 606, and 620 protein spots were detected, respectively. And on the gels V4N to V13N, 420, 455, 500, 610, 630, 635, 643, 640, 660, and 684 protein spots were detected, respectively. The number of protein spots increased gradually with the development of larvae in two groups. As shown in Figure1, the protein spots on the V7 gel were marked by number and further analyzed by a Voyager DE PRO MALDI-TOF-MS. After being searched by GPMAW software, proteins with more than five matched peptides and the coverage of peptides larger than 20% were accepted. The results were summarized in Table 1. A total of 96 proteins were identified. Among them, protein spots no. 95 and 96 were only found in nonfeeding midgut. When compared to the data obtained from the previous studies [10, 11, 15, 16], 69 proteins, accounting for 71% of all proteins, were found at the first time in our present study. Based on the molecular function and GO annotation (see below), these 96 proteins were classified into 12 categories including: (1) eight are for nutrition storage, (2) three for cell growth, (3) sixteen for protein metabolism, (4) three for carbohydrate metabolism, (5) nine for lipid metabolism, (6) eighteen for energy metabolism, (7) three are immune-related proteins, (8) two are antioxidation proteins, (9) eleven are muscle-related proteins, (10) five are for hormone metabolism, (11) sixteen for other functions, and (12) two are unannotated proteins. The identified proteins have been submitted to http://silkworm.swu.edu.cn/cgi-bin/wego/index.pl for GO (gene ontology) annotation. As shown in Figure 2, 96 proteins were classified into three groups including cellular component, molecular function containing binding, catalytic, transporter, electron carrier, and antioxidant, and biological processes involved in metabolism, cellular process, and biological regulation.Table 1
List of proteins identified from midgut ofBombyx mori.
Protein classificationSpot no.Protein descriptionAccession numberMW (kDa)/ theory PICoverage(%)V4/V4N(F/N)V5/V5N(F/N)V6/VN (F/N)V7/V7N(F/N)V8/V8N(F/N)V9N(N)V10N(N)V11N  (N)V12N(N)V13N (N)Nutrition storage830 kDa protein [Bombyx mori]ACO3575029.83/5.9023.8+/++/++/++/++/++++++3730 K lipoprotein precursor [Bombyx mori]CAA3043330.26/7.0832.6-/-+/-+/+-/--/-−−−−−38Low-molecular-mass 30 kDa lipoprotein 21G1Precursor [Bombyx mori]CAA3853132.98/6.1926.5-/--/--/-+/++/++++++57Low-molecular-mass 30 kDa lipoprotein 19G1 [Bombyx mori]CAA3853332.41/9.5865.5-/--/--/-+/-+/+−−−−−60Low-molecular-mass 30 kDa lipoprotein 21G1 Precursor [Bombyx mori]Q0080130.35/6.3850.2-/--/--/-+/++/++++++82Low molecular 30 kDa lipoprotein PBMHPC-19 precursor [Bombyx mori]NP_00109519729.41/7.2033.9-/--/--/-+/++++/++++++85Low-molecular-mass 30 kDa lipoprotein 19G1 [Bombyx mori]CAA3853329.56/7.9139.1+/--/++/++/++/++++++90Mature 30 K lipoprotein [Bombyx mori]CAA3043428.57/6.3744.0-/--/--/-+/++/++++−+Cell growth45Imaginal disk growth factor [Bombyx mori]BAF7362348.39/7.7928.3-/--/--/-++/+++/++++++++++++50Actin-depolymerizing factor1 [Bombyx mori]NP_00109327817.23/6.1941.2+/++/++/++/+-/++++++75Transferrin [Bombyx mori]NP_00103701477.20/6.8427.5-/-+/++/++/++/++++++Protein metabolism360 kDa heat shock protein, mitochondrial [Drosophila melanogaster]NP_51111561.23/5.3723.6+/++/++/++/++/++++++7Protein phosphatase 1 catalytic subunit [Bombyx mori]NP_00104048032.28/5.8040.5+/++/++/++/++/++++++18Acireductone dioxygenase [Aedes aegypti]XP_001649899.21.91/4.8540.8+/-+/-+/++/++/++++++28Glyoxylate reductase/hydroxypyruvate reductase [Culex quinquefasciatus]XP_00186566037.34/7.9432.4-/-+/-+/++/++/++++++40Protein disulfide isomerase [Aedes aegypti]ABF1837846.34/5.1133.4-/--/-+/++/++/++++++43, 44DnaJ-like protein [Bombyx mori]NP_00104018559.96/7.3229.9-/--/--/-+/++/++++++46Dihydrolipoamide dehydrogenase [Bombyx mori]AAM9325552.94/8.9027.8-/--/--/-+/++/++++++47Aspartate aminotransferase [Bombyx mori]ABF5121440.32/5.8732.3-/--/--/-+/-+/++++++52Protein disulfide-isomerase-like protein ERp57 [Bombyx mori]NP_00103699755.48/5.1737.1-/-+/++/+++++/+++/+++++++58Ribosomal protein P0 [Bombyx mori]NP_00103712334.24/5.6141.8-/--/--/-+/++/++++++63Ubiquitin carboxyl-terminal hydrolase [Drosophila melanogaster]NP_47694025.52/4.8729.6+/++/++/++/+-/++++++70Ribosomal protein P2 [Bombyx mori]NP_00103721321.20/5.7938.3+/++/++/++/++/++++++74Heat shock protein 70 [Bombyx mori]NP_00103683773.21/4.9826.0-/-+/++/+++/+++/++++++++++87Peptidylprolyl isomerase B [Bombyx mori]NP_00104047922.39/8.9142.9-/--/--/-+/++/++++++9510 kDa heat shock protein, mitochondrial [Lepeophtheirus salmonis]ACO1179611.18/7.6849.5-/--/--/--/+-/+++++−Carbohydrate metabolism6Enolase [Bombyx mori]NP_00109183147.19/5.5541.6+/++/++/++/++/++++++13Triosephosphate isomerase [Bombyx mori]AAU3418526.94/5.6050.8+/++/++/++/++/++++++31Glyceraldehyde-3-phosphate dehydrogenase [Bombyx mori]NP_00103738635.42/8.5529.2-/-+/-+/++/++/++++++Lipid metabolism30Enoyl-CoA hydratase precursor 1 [Bombyx mori]ABD3610732.13/8.3345.9-/-+/-+/++/++/++++++35Fatty acid binding protein [Bombyx mori]NP_00103757415.01/5.0824.2+/++/+++/+++/+++/+++++++++36Apolipophorin-III precursor [Bombyx mori]AAB0285220.82/9.5343.0+/++/++/++/+-/++++++48Acyl-coenzyme A dehydrogenase [Bombyx mori]NP_00103767247.53/7.1622.7-/--/--/-+/++/++++++54Acyl-coenzyme A dehydrogenase [Aedes aegypti]EAT4276635.81/7.1844.5-/--/--/-+/-+/++++++563-hydroxy acyl-CoA dehydrogenase [Bombyx mori]NP_00104041433.80/8.8728.4-/--/--/-+/++/++++−+65Lysophospholipase [Bombyx mori]NP_00104025525.97/6.9436.2-/--/--/-+/++/++++++79Acetoacetyl-CoA thiolase [Bombyx mori]NP_00109329640.00/7.5642.6-/--/--/-+/+-/++++++91Phosphatidylethanolamine binding protein isoform 2 [Bombyx mori]NP_00109326721.92/5.9447.7-/--/--/-+/++/++++−+Energy metabolism1H+-transporting ATP synthase beta subunit isoform 2 [Bombyx mori]NP_00104170554.89/5.2037.8+/++/++/++/++/++++++5H+-transporting ATP synthase beta subunit isoform 1 [Bombyx mori]NP_00104045055.04/5.1333.9+/++/++/++/++/++++++15Abnormal wing disc-like protein [Bombyx mori]NP_00109328417.37/7.3048.1++/++++/+++++/++++/+++/+++++++++++++++20Vacuolar ATP synthase subunit F [Bombyx mori]NP_00104044813.92/5.9250.0+/-+/-+/++/++/++++++22Cytochrome c oxidase subunit Va [Bombyx mori]NP_00110674217.24/5.5931.8+/++/++/++/++/++++−+23H+-transporting ATP synthase delta subunit [Bombyx mori]NP_00109309117.16/5.6649.1+/++/++/++/++/++++++33H+-transporting ATP synthase subunit d [Bombyx mori]NP_00109327920.20/5.4939.7+/++/++/+-/--/-−−−−−51Vacuolar ATPase B subunit [Bombyx mori]ACE7827155.19/5.2038.1-/--/-+/++++/+++/+++++++53Ubiquinol-cytochrome c reductase [Bombyx mori]ABF6022529.43/8.5223.1-/--/--/-+/++/++++++55Vacuolar ATP synthase catalytic subunit A [Bombyx mori]NP_00109182968.60/5.1338.7-/-+/++/++++/+++/++++++++++++64ATP synthase subunit alpha [Bombyx mori]BAA0372528.66/8.4341.9-/--/--/-+/++/++++++66ATP synthase subunit alpha [Bombyx mori]BAA0372528.66/8.4341.9-/--/--/-++/+++/++++++++714-Hydroxybutyrate CoA-transferase, putative [Aedes aegypti]EAT4492143.62/6.6128.1-/--/--/-+++/+++/++++++++++++76Arginine kinase [Bombyx mori]ACI0104840.33/5.8436.9+/++/++/++/++/++++++77Ubiquinol-cytochrome c reductase core protein II [Bombyx mori]NP_00110622546.18/9.5522.4-/--/--/-+/++/++++++86ATP synthase subunit alpha [Bombyx mori]Q0740528.66/8.4341.9-/--/--/-+/++/+−−−++88ATP synthase subunit alpha [Bombyx mori]Q0740528.66/8.4341.9+/++/++/++/++/++++++92NADH dehydrogenase (ubiquinone) Fe-S protein 8 [Bombyx mori]NP_00104031625.49/6.1431.2-/--/--/-+/++/++++−+Immune related4Immune-related Hdd13 [Hyphantria cunea]AAD0928128.38/6.8660.8+/++/++/++/++/++++++16Cyclophilin A [Bombyx mori]NP_00103730118.10/7.9773.9+/++/++/++/++/++++++17Cyclophilin  [Bombyx mori]ABA5417318.12/7.9762.4+/++/++/++/--/-−−−−−Antioxidation34Glutathione S-transferase 10 [Bombyx mori]NP_00110846624.47/8.1126.7-/-+/--/--/--/-−−−−−89Thiol peroxiredoxin [Bombyx mori]AAR1542022.08/6.1035.9+/++/++/++/+-/++++++Muscle related2Actin A3 [Bombyx mori]CAA2819242.23/5.3831.9+/++/++++/++++++/++++/+++++++++++++9Actin, muscle-type A1 [Bombyx mori]Gi|113216|42.27/5.1929.3+/++/++/++/++/++++++14Actin, cytoplasmic A4 [Bombyx mori]NP_00111972730.05/5.8024.7+/++/++/++/++/++++++21Myosin 1 light chain [Lonomia oblique]AAV9141122.77/4.3133.0+/++/++/++/++/++++++39Actin, cytoplasmic A4 [Bombyx mori]NP_00111972730.05/5.8025.5-/--/--/-+/--/++++++59Beta-tubulin [Bombyx mori]NP_00103688725.18/4.5130.8+/++/++/++/++/++++++61Myosin heavy chain [Bombyx mori]ACQ7282596.12/4.9820.2-/--/--/-+/-+/-−−+++62Myosin light chain 2 [Bombyx mori]ACF5973522.04/4.5232.8-/--/--/-+/--/++++++67Muscular protein 20 [Bombyx mori]NP_00104047620.30/8.8145.7-/--/--/-+/-+/-−−−++68Transgelin [Bombyx mori]NP_00104037220.91/8.4085.0-/--/--/++/++/++++++93Myosin light chain 2 [Bombyx mori]NP_00109181322.04/4.5230.3-/--/--/-+/++/++++++Hormone metabolism24, 25Juvenile hormone diol kinase [Bombyx mori]NP_00103708020.79/4.4239.5+/++/++/++/++/++++++83Putative farnesoic acid O-methyl transferase [Bombyx mori]NP_00109181526.83/6.7335.1-/--/--/-+/++/++++++94Putative farnesoic acid O-methyl transferase [Bombyx mori]NP_00109181524.26/5.0328.7+/++/++/++/++/++++++96Diapause hormone precursor [Bombyx mori]BAA0375522.39/9.4460.9-/--/--/+-/--/+−−−++Other proteins1027 kDa glycoprotein, p27K [Bombyx mori]Q8T11325.58/4.9525.1+/++/++/++/++/++++++11Cellular retinoic acid binding protein [Bombyx mori]NP_00103736414.97/5.5356.8+/++/++/++/++/++++++12Stardust, isoform B [Drosophila melanogaster]NP_57246332.15/5.8126.3+/++/++/++/++/++++++19NADPH oxidase [Culex quinquefasciatus]XP_00185905930.17/8.3535.1++/+++/+++/++-/--/-−−−−−29Similar to CG31790-PA [Tribolium castaneum]XP_97196050.42/6.0122.0-/-+/-+/++/--/++++++32Serpin-2 [Bombyx mori]NP_00103702141.86/4.7221.1-/-+/++/++/+-/++++++41GE24263 [Drosophila yakuba]EDW9753047.04/5.9529.3-/--/--/-+/++/++++++42PREDICTED: hypothetical protein [Nasonia vitripennis].XP_00160257351.08/5.6828.5-/--/-+/++/++/++++++49AGAP011050-PA [Anopheles gambiae]EAA4551134.16/6.0939.7-/--/--/-+/++/++++++69GI18783 [Drosophila mojavensis]EDW0993918.98/4.7439.9+/++/++/++/++/++++++72Mitochondrial aldehyde dehydrogenase [Bombyx mori]NP_00104047556.24/7.5529.6-/--/--/-+/++/++++++73Soluble guanylyl cyclase beta-1 subunit [Manduca sexta]AAC6126436.77/5.5520.1-/--/--/-+/-+/-−−−−+78Putative acetyl transferase [Bombyx mori]BAH9656142.01/7.8431.0-/--/--/-+/++/++++++80PREDICTED: similar to CG10122-PA [Tribolium castaneum]XP_97091843.23/5.0320.2-/--/-+/++/++/+−−+−+81Ef1alpha-like factor isoform 1 [Bombyx mori]NP_00104011831.48/7.8446.4-/--/-+/++/++/++++++84Voltage-dependent anion-selective channel [Drosophila melanogaster]NP_00103389930.13/7.4535.5-/--/--/-+/++/++++++Unannotated26Unknown [Bombyx mori]BGIBMGA00411612.80/4.4943.1+/++/++/+-/--/-−−−−−27Unknown [Bombyx mori]BGIBMGA01336053.34/5.8033.1-/-+/-+/++/++/++++++Total9637/3447/4052/5485/8775/848080827885“+” and “−” indicate protein with and without expression, respectively. The number of “+” means the x-fold up expression compared with the “+”. “Vx” means fifth-instar day-x larvae, “F” means feeding, and “N” means nonfeeding.Figure 1
The 2-DE profile of midgut proteins at the fifth-instar day 7 larvae. A total of 300μg protein from midgut was applied onto a pH 3–10 Immobiline dryStrip (18 cm) for the first-dimensional electrophoresis and then separated by 12.5% SDS-PAGE. The 2-DE gels were stained by silver. Proteins identified in this study are indicated, and the numbers corresponded to those are in Table 1.Figure 2
GO categories of the identified proteins by Wego software.


### 3.2. Changes in Protein Pattern in the Development of the Midguts of Normal-Feeding Silkworm Larvae

As shown in Figure1 and Supplemental data 1, the composition of the silkworm midgut proteins was different from V4 to V8. The number of protein spots, especially those ranging in size from 30 kDa to 90 kDa at higher pI region, gradually increased with the growth of larvae. Table 1 showed that the abundances of 64 midgut proteins were different from V4 to V8, mainly including those involved in nutrition storage, protein/lipid/energy metabolism, and muscle-related proteins. Most prominent among them were thirteen proteins with more than two-fold differential expression. Two spots corresponding to a juvenile hormone diol kinase (JHDK) and a fatty-acid-binding protein were present from V4 to V8, and displayed the highest expression level in V6 (Figures 3(a) and 3(b)). An abnormal wing disc-like protein had a peak in V7 (Figure 3(c)). From Figures 3(d) to 3(i), one can see that ten proteins among which two were low-molecular-mass 30 kDa lipoproteins 19G1 and PBMHPC-19 precursor, three were ATP synthase subunits, and five were, respectively, protein disulfide-isomerase-like protein (PDI) ERp57, enolase, arginine kinase, Enoyl-CoA hydratase precursor 1, and actin were expressed at relatively higher levels in V8. From Table 1, we also found that 6 proteins (spot no. 17, 19, 26, 33, 34, and 37), a 30 K lipoprotein precursor, an H+-transporting ATP synthase subunit, a cyclophilin, a glutathione S-transferase10, a NADPH oxidase, and an unknown protein, were barely detectable before metamorphosis in feeding larvae.Figure 3
The partial enlargement of midgut proteins changed with development in normal feeding silkworm. Red arrows pointed out the target proteins and their names were displayed at right. Vx means fifth-instar day-x larvae.


### 3.3. Midgut Proteins in Response to Starvation

After being fed with mulberry leaves for three days, a part of these fed silkworm larvae were separated for starvation treatment. The starved larvae survived for 10 days. Most of them were dead during starvation. At the day 13 (V13N), a few larvae struggled to pupate. The midguts were carefully collected at 24-hour intervals from the starved larvae until day 13 (V13N). We subsequently investigated these 10 protein samples from V4N to V13N by 2D electrophoresis coupled with MALDI-TOF MS. Supplemental Data 2 showed that the protein profiles of V4N and V5N were quite similar to each other, and the number of the protein spots on both gels was decreased as starvation time increased, compared to those in V4 and V5. However, the spot number was dramaticly increased at V6N compared to that in V5N. On the V6N to V12N gels, the numbers of protein spots were also increased and varied from 610 to 660. The protein distribution pattern of V6N to V12N resembled those of V5 to V8, suggesting that the starved larvae might keep their midgut intact and functional for a period of time. It is also important to note that the spot number was increased at V13N, up to 684, presumably coincident with the striking alteration of starved larvae at that stage.The proteins identified in starved larval midgut were summarized in Table1. As shown in Table 1, over 50% of identified proteins in nonfeeding silkworm larvae had a similar variation trend to that in normal feeding ones. On the contrary, 17 proteins delay their presentation compared with those in feeding group. Moreover, nine proteins (spot no. 29, 32, 36, 39, 50, 62, 63, 79, and 89) still present in nonfeeding larvae were undetectable before metamorphosis in feeding larvae; as a representative data, the expression pattern of actin-depolymerizing factor 1 was shown in Figure 4(g).Figure 4
The partial enlargement of the midgut proteins in nonfeeding silkworms. Red arrows pointed out the target proteins and their names were displayed at right. Vx means fifth-instar day-x larvae.In Figures4(a) and 4(b), a 10 kDa heat shock protein and a diapause hormone precursor were only detected in the nonfeeding groups. Real-time PCR was carried out to reveal changing expression patterns of two genes encoding the 10 kDa heat shock protein and diapauses hormone precursor. The data showed that this small heat shock protein gene was expressed in the normal feeding and nonfeeding larvae, but its peak signal can be detected after spending nine days in starvation (Supplementary Data 3A). As shown in supplementary data 3B, diapauses hormone precursor gene was expressed only in starved larvae at V5N and V13N.Besides these two proteins, the glutathione S-transferase 10 was identified at V5 but not found from V4N to V13N (Figure4(c)). The expressions of a PDI-like protein ERp57 and a vacuolar ATPase B subunit were upregulated at V6N and thereafter downregulated and returned to the same level as feeding silkworm (Figure 4(d)). The amounts of H+-transporting ATP synthase beta subunit isoform 2 and imaginal disc growth factor were at first lower in nonfeeding larvae but finally increased to the same level as in feeding larvae (Figures 4(d) and 4(e)). In Figure 4(f), we found that the amount of an abnormal wing disc-like protein decreased with development progress in feeding larvae, but after starvation, it kept a stable level (Figure 4(g)).


## 3.1.  2-DE Profiles of Silkworm Larval Midgut Proteins

The silkworm midgut proteins were extracted from the fifth-instar day-4 (V4) to day-8 (V8) feeding larvae and from the fifth-instar day-4 (V4N) to day-13 (V13N) nonfeeding larvae. A total of 15 extractions were separated by 2D technology, as shown in Supplemental Data 1 and 2 in Supplemental Material available at doi:10.1155/2011/876064. The protein spots were then detected using ImageMaster 2D platinum 6.0 software. On the gels V4 to V8, 439, 461, 498, 606, and 620 protein spots were detected, respectively. And on the gels V4N to V13N, 420, 455, 500, 610, 630, 635, 643, 640, 660, and 684 protein spots were detected, respectively. The number of protein spots increased gradually with the development of larvae in two groups. As shown in Figure1, the protein spots on the V7 gel were marked by number and further analyzed by a Voyager DE PRO MALDI-TOF-MS. After being searched by GPMAW software, proteins with more than five matched peptides and the coverage of peptides larger than 20% were accepted. The results were summarized in Table 1. A total of 96 proteins were identified. Among them, protein spots no. 95 and 96 were only found in nonfeeding midgut. When compared to the data obtained from the previous studies [10, 11, 15, 16], 69 proteins, accounting for 71% of all proteins, were found at the first time in our present study. Based on the molecular function and GO annotation (see below), these 96 proteins were classified into 12 categories including: (1) eight are for nutrition storage, (2) three for cell growth, (3) sixteen for protein metabolism, (4) three for carbohydrate metabolism, (5) nine for lipid metabolism, (6) eighteen for energy metabolism, (7) three are immune-related proteins, (8) two are antioxidation proteins, (9) eleven are muscle-related proteins, (10) five are for hormone metabolism, (11) sixteen for other functions, and (12) two are unannotated proteins. The identified proteins have been submitted to http://silkworm.swu.edu.cn/cgi-bin/wego/index.pl for GO (gene ontology) annotation. As shown in Figure 2, 96 proteins were classified into three groups including cellular component, molecular function containing binding, catalytic, transporter, electron carrier, and antioxidant, and biological processes involved in metabolism, cellular process, and biological regulation.Table 1
List of proteins identified from midgut ofBombyx mori.
Protein classificationSpot no.Protein descriptionAccession numberMW (kDa)/ theory PICoverage(%)V4/V4N(F/N)V5/V5N(F/N)V6/VN (F/N)V7/V7N(F/N)V8/V8N(F/N)V9N(N)V10N(N)V11N  (N)V12N(N)V13N (N)Nutrition storage830 kDa protein [Bombyx mori]ACO3575029.83/5.9023.8+/++/++/++/++/++++++3730 K lipoprotein precursor [Bombyx mori]CAA3043330.26/7.0832.6-/-+/-+/+-/--/-−−−−−38Low-molecular-mass 30 kDa lipoprotein 21G1Precursor [Bombyx mori]CAA3853132.98/6.1926.5-/--/--/-+/++/++++++57Low-molecular-mass 30 kDa lipoprotein 19G1 [Bombyx mori]CAA3853332.41/9.5865.5-/--/--/-+/-+/+−−−−−60Low-molecular-mass 30 kDa lipoprotein 21G1 Precursor [Bombyx mori]Q0080130.35/6.3850.2-/--/--/-+/++/++++++82Low molecular 30 kDa lipoprotein PBMHPC-19 precursor [Bombyx mori]NP_00109519729.41/7.2033.9-/--/--/-+/++++/++++++85Low-molecular-mass 30 kDa lipoprotein 19G1 [Bombyx mori]CAA3853329.56/7.9139.1+/--/++/++/++/++++++90Mature 30 K lipoprotein [Bombyx mori]CAA3043428.57/6.3744.0-/--/--/-+/++/++++−+Cell growth45Imaginal disk growth factor [Bombyx mori]BAF7362348.39/7.7928.3-/--/--/-++/+++/++++++++++++50Actin-depolymerizing factor1 [Bombyx mori]NP_00109327817.23/6.1941.2+/++/++/++/+-/++++++75Transferrin [Bombyx mori]NP_00103701477.20/6.8427.5-/-+/++/++/++/++++++Protein metabolism360 kDa heat shock protein, mitochondrial [Drosophila melanogaster]NP_51111561.23/5.3723.6+/++/++/++/++/++++++7Protein phosphatase 1 catalytic subunit [Bombyx mori]NP_00104048032.28/5.8040.5+/++/++/++/++/++++++18Acireductone dioxygenase [Aedes aegypti]XP_001649899.21.91/4.8540.8+/-+/-+/++/++/++++++28Glyoxylate reductase/hydroxypyruvate reductase [Culex quinquefasciatus]XP_00186566037.34/7.9432.4-/-+/-+/++/++/++++++40Protein disulfide isomerase [Aedes aegypti]ABF1837846.34/5.1133.4-/--/-+/++/++/++++++43, 44DnaJ-like protein [Bombyx mori]NP_00104018559.96/7.3229.9-/--/--/-+/++/++++++46Dihydrolipoamide dehydrogenase [Bombyx mori]AAM9325552.94/8.9027.8-/--/--/-+/++/++++++47Aspartate aminotransferase [Bombyx mori]ABF5121440.32/5.8732.3-/--/--/-+/-+/++++++52Protein disulfide-isomerase-like protein ERp57 [Bombyx mori]NP_00103699755.48/5.1737.1-/-+/++/+++++/+++/+++++++58Ribosomal protein P0 [Bombyx mori]NP_00103712334.24/5.6141.8-/--/--/-+/++/++++++63Ubiquitin carboxyl-terminal hydrolase [Drosophila melanogaster]NP_47694025.52/4.8729.6+/++/++/++/+-/++++++70Ribosomal protein P2 [Bombyx mori]NP_00103721321.20/5.7938.3+/++/++/++/++/++++++74Heat shock protein 70 [Bombyx mori]NP_00103683773.21/4.9826.0-/-+/++/+++/+++/++++++++++87Peptidylprolyl isomerase B [Bombyx mori]NP_00104047922.39/8.9142.9-/--/--/-+/++/++++++9510 kDa heat shock protein, mitochondrial [Lepeophtheirus salmonis]ACO1179611.18/7.6849.5-/--/--/--/+-/+++++−Carbohydrate metabolism6Enolase [Bombyx mori]NP_00109183147.19/5.5541.6+/++/++/++/++/++++++13Triosephosphate isomerase [Bombyx mori]AAU3418526.94/5.6050.8+/++/++/++/++/++++++31Glyceraldehyde-3-phosphate dehydrogenase [Bombyx mori]NP_00103738635.42/8.5529.2-/-+/-+/++/++/++++++Lipid metabolism30Enoyl-CoA hydratase precursor 1 [Bombyx mori]ABD3610732.13/8.3345.9-/-+/-+/++/++/++++++35Fatty acid binding protein [Bombyx mori]NP_00103757415.01/5.0824.2+/++/+++/+++/+++/+++++++++36Apolipophorin-III precursor [Bombyx mori]AAB0285220.82/9.5343.0+/++/++/++/+-/++++++48Acyl-coenzyme A dehydrogenase [Bombyx mori]NP_00103767247.53/7.1622.7-/--/--/-+/++/++++++54Acyl-coenzyme A dehydrogenase [Aedes aegypti]EAT4276635.81/7.1844.5-/--/--/-+/-+/++++++563-hydroxy acyl-CoA dehydrogenase [Bombyx mori]NP_00104041433.80/8.8728.4-/--/--/-+/++/++++−+65Lysophospholipase [Bombyx mori]NP_00104025525.97/6.9436.2-/--/--/-+/++/++++++79Acetoacetyl-CoA thiolase [Bombyx mori]NP_00109329640.00/7.5642.6-/--/--/-+/+-/++++++91Phosphatidylethanolamine binding protein isoform 2 [Bombyx mori]NP_00109326721.92/5.9447.7-/--/--/-+/++/++++−+Energy metabolism1H+-transporting ATP synthase beta subunit isoform 2 [Bombyx mori]NP_00104170554.89/5.2037.8+/++/++/++/++/++++++5H+-transporting ATP synthase beta subunit isoform 1 [Bombyx mori]NP_00104045055.04/5.1333.9+/++/++/++/++/++++++15Abnormal wing disc-like protein [Bombyx mori]NP_00109328417.37/7.3048.1++/++++/+++++/++++/+++/+++++++++++++++20Vacuolar ATP synthase subunit F [Bombyx mori]NP_00104044813.92/5.9250.0+/-+/-+/++/++/++++++22Cytochrome c oxidase subunit Va [Bombyx mori]NP_00110674217.24/5.5931.8+/++/++/++/++/++++−+23H+-transporting ATP synthase delta subunit [Bombyx mori]NP_00109309117.16/5.6649.1+/++/++/++/++/++++++33H+-transporting ATP synthase subunit d [Bombyx mori]NP_00109327920.20/5.4939.7+/++/++/+-/--/-−−−−−51Vacuolar ATPase B subunit [Bombyx mori]ACE7827155.19/5.2038.1-/--/-+/++++/+++/+++++++53Ubiquinol-cytochrome c reductase [Bombyx mori]ABF6022529.43/8.5223.1-/--/--/-+/++/++++++55Vacuolar ATP synthase catalytic subunit A [Bombyx mori]NP_00109182968.60/5.1338.7-/-+/++/++++/+++/++++++++++++64ATP synthase subunit alpha [Bombyx mori]BAA0372528.66/8.4341.9-/--/--/-+/++/++++++66ATP synthase subunit alpha [Bombyx mori]BAA0372528.66/8.4341.9-/--/--/-++/+++/++++++++714-Hydroxybutyrate CoA-transferase, putative [Aedes aegypti]EAT4492143.62/6.6128.1-/--/--/-+++/+++/++++++++++++76Arginine kinase [Bombyx mori]ACI0104840.33/5.8436.9+/++/++/++/++/++++++77Ubiquinol-cytochrome c reductase core protein II [Bombyx mori]NP_00110622546.18/9.5522.4-/--/--/-+/++/++++++86ATP synthase subunit alpha [Bombyx mori]Q0740528.66/8.4341.9-/--/--/-+/++/+−−−++88ATP synthase subunit alpha [Bombyx mori]Q0740528.66/8.4341.9+/++/++/++/++/++++++92NADH dehydrogenase (ubiquinone) Fe-S protein 8 [Bombyx mori]NP_00104031625.49/6.1431.2-/--/--/-+/++/++++−+Immune related4Immune-related Hdd13 [Hyphantria cunea]AAD0928128.38/6.8660.8+/++/++/++/++/++++++16Cyclophilin A [Bombyx mori]NP_00103730118.10/7.9773.9+/++/++/++/++/++++++17Cyclophilin  [Bombyx mori]ABA5417318.12/7.9762.4+/++/++/++/--/-−−−−−Antioxidation34Glutathione S-transferase 10 [Bombyx mori]NP_00110846624.47/8.1126.7-/-+/--/--/--/-−−−−−89Thiol peroxiredoxin [Bombyx mori]AAR1542022.08/6.1035.9+/++/++/++/+-/++++++Muscle related2Actin A3 [Bombyx mori]CAA2819242.23/5.3831.9+/++/++++/++++++/++++/+++++++++++++9Actin, muscle-type A1 [Bombyx mori]Gi|113216|42.27/5.1929.3+/++/++/++/++/++++++14Actin, cytoplasmic A4 [Bombyx mori]NP_00111972730.05/5.8024.7+/++/++/++/++/++++++21Myosin 1 light chain [Lonomia oblique]AAV9141122.77/4.3133.0+/++/++/++/++/++++++39Actin, cytoplasmic A4 [Bombyx mori]NP_00111972730.05/5.8025.5-/--/--/-+/--/++++++59Beta-tubulin [Bombyx mori]NP_00103688725.18/4.5130.8+/++/++/++/++/++++++61Myosin heavy chain [Bombyx mori]ACQ7282596.12/4.9820.2-/--/--/-+/-+/-−−+++62Myosin light chain 2 [Bombyx mori]ACF5973522.04/4.5232.8-/--/--/-+/--/++++++67Muscular protein 20 [Bombyx mori]NP_00104047620.30/8.8145.7-/--/--/-+/-+/-−−−++68Transgelin [Bombyx mori]NP_00104037220.91/8.4085.0-/--/--/++/++/++++++93Myosin light chain 2 [Bombyx mori]NP_00109181322.04/4.5230.3-/--/--/-+/++/++++++Hormone metabolism24, 25Juvenile hormone diol kinase [Bombyx mori]NP_00103708020.79/4.4239.5+/++/++/++/++/++++++83Putative farnesoic acid O-methyl transferase [Bombyx mori]NP_00109181526.83/6.7335.1-/--/--/-+/++/++++++94Putative farnesoic acid O-methyl transferase [Bombyx mori]NP_00109181524.26/5.0328.7+/++/++/++/++/++++++96Diapause hormone precursor [Bombyx mori]BAA0375522.39/9.4460.9-/--/--/+-/--/+−−−++Other proteins1027 kDa glycoprotein, p27K [Bombyx mori]Q8T11325.58/4.9525.1+/++/++/++/++/++++++11Cellular retinoic acid binding protein [Bombyx mori]NP_00103736414.97/5.5356.8+/++/++/++/++/++++++12Stardust, isoform B [Drosophila melanogaster]NP_57246332.15/5.8126.3+/++/++/++/++/++++++19NADPH oxidase [Culex quinquefasciatus]XP_00185905930.17/8.3535.1++/+++/+++/++-/--/-−−−−−29Similar to CG31790-PA [Tribolium castaneum]XP_97196050.42/6.0122.0-/-+/-+/++/--/++++++32Serpin-2 [Bombyx mori]NP_00103702141.86/4.7221.1-/-+/++/++/+-/++++++41GE24263 [Drosophila yakuba]EDW9753047.04/5.9529.3-/--/--/-+/++/++++++42PREDICTED: hypothetical protein [Nasonia vitripennis].XP_00160257351.08/5.6828.5-/--/-+/++/++/++++++49AGAP011050-PA [Anopheles gambiae]EAA4551134.16/6.0939.7-/--/--/-+/++/++++++69GI18783 [Drosophila mojavensis]EDW0993918.98/4.7439.9+/++/++/++/++/++++++72Mitochondrial aldehyde dehydrogenase [Bombyx mori]NP_00104047556.24/7.5529.6-/--/--/-+/++/++++++73Soluble guanylyl cyclase beta-1 subunit [Manduca sexta]AAC6126436.77/5.5520.1-/--/--/-+/-+/-−−−−+78Putative acetyl transferase [Bombyx mori]BAH9656142.01/7.8431.0-/--/--/-+/++/++++++80PREDICTED: similar to CG10122-PA [Tribolium castaneum]XP_97091843.23/5.0320.2-/--/-+/++/++/+−−+−+81Ef1alpha-like factor isoform 1 [Bombyx mori]NP_00104011831.48/7.8446.4-/--/-+/++/++/++++++84Voltage-dependent anion-selective channel [Drosophila melanogaster]NP_00103389930.13/7.4535.5-/--/--/-+/++/++++++Unannotated26Unknown [Bombyx mori]BGIBMGA00411612.80/4.4943.1+/++/++/+-/--/-−−−−−27Unknown [Bombyx mori]BGIBMGA01336053.34/5.8033.1-/-+/-+/++/++/++++++Total9637/3447/4052/5485/8775/848080827885“+” and “−” indicate protein with and without expression, respectively. The number of “+” means the x-fold up expression compared with the “+”. “Vx” means fifth-instar day-x larvae, “F” means feeding, and “N” means nonfeeding.Figure 1
The 2-DE profile of midgut proteins at the fifth-instar day 7 larvae. A total of 300μg protein from midgut was applied onto a pH 3–10 Immobiline dryStrip (18 cm) for the first-dimensional electrophoresis and then separated by 12.5% SDS-PAGE. The 2-DE gels were stained by silver. Proteins identified in this study are indicated, and the numbers corresponded to those are in Table 1.Figure 2
GO categories of the identified proteins by Wego software.


## 3.2. Changes in Protein Pattern in the Development of the Midguts of Normal-Feeding Silkworm Larvae

As shown in Figure1 and Supplemental data 1, the composition of the silkworm midgut proteins was different from V4 to V8. The number of protein spots, especially those ranging in size from 30 kDa to 90 kDa at higher pI region, gradually increased with the growth of larvae. Table 1 showed that the abundances of 64 midgut proteins were different from V4 to V8, mainly including those involved in nutrition storage, protein/lipid/energy metabolism, and muscle-related proteins. Most prominent among them were thirteen proteins with more than two-fold differential expression. Two spots corresponding to a juvenile hormone diol kinase (JHDK) and a fatty-acid-binding protein were present from V4 to V8, and displayed the highest expression level in V6 (Figures 3(a) and 3(b)). An abnormal wing disc-like protein had a peak in V7 (Figure 3(c)). From Figures 3(d) to 3(i), one can see that ten proteins among which two were low-molecular-mass 30 kDa lipoproteins 19G1 and PBMHPC-19 precursor, three were ATP synthase subunits, and five were, respectively, protein disulfide-isomerase-like protein (PDI) ERp57, enolase, arginine kinase, Enoyl-CoA hydratase precursor 1, and actin were expressed at relatively higher levels in V8. From Table 1, we also found that 6 proteins (spot no. 17, 19, 26, 33, 34, and 37), a 30 K lipoprotein precursor, an H+-transporting ATP synthase subunit, a cyclophilin, a glutathione S-transferase10, a NADPH oxidase, and an unknown protein, were barely detectable before metamorphosis in feeding larvae.Figure 3
The partial enlargement of midgut proteins changed with development in normal feeding silkworm. Red arrows pointed out the target proteins and their names were displayed at right. Vx means fifth-instar day-x larvae.


## 3.3. Midgut Proteins in Response to Starvation

After being fed with mulberry leaves for three days, a part of these fed silkworm larvae were separated for starvation treatment. The starved larvae survived for 10 days. Most of them were dead during starvation. At the day 13 (V13N), a few larvae struggled to pupate. The midguts were carefully collected at 24-hour intervals from the starved larvae until day 13 (V13N). We subsequently investigated these 10 protein samples from V4N to V13N by 2D electrophoresis coupled with MALDI-TOF MS. Supplemental Data 2 showed that the protein profiles of V4N and V5N were quite similar to each other, and the number of the protein spots on both gels was decreased as starvation time increased, compared to those in V4 and V5. However, the spot number was dramaticly increased at V6N compared to that in V5N. On the V6N to V12N gels, the numbers of protein spots were also increased and varied from 610 to 660. The protein distribution pattern of V6N to V12N resembled those of V5 to V8, suggesting that the starved larvae might keep their midgut intact and functional for a period of time. It is also important to note that the spot number was increased at V13N, up to 684, presumably coincident with the striking alteration of starved larvae at that stage.The proteins identified in starved larval midgut were summarized in Table1. As shown in Table 1, over 50% of identified proteins in nonfeeding silkworm larvae had a similar variation trend to that in normal feeding ones. On the contrary, 17 proteins delay their presentation compared with those in feeding group. Moreover, nine proteins (spot no. 29, 32, 36, 39, 50, 62, 63, 79, and 89) still present in nonfeeding larvae were undetectable before metamorphosis in feeding larvae; as a representative data, the expression pattern of actin-depolymerizing factor 1 was shown in Figure 4(g).Figure 4
The partial enlargement of the midgut proteins in nonfeeding silkworms. Red arrows pointed out the target proteins and their names were displayed at right. Vx means fifth-instar day-x larvae.In Figures4(a) and 4(b), a 10 kDa heat shock protein and a diapause hormone precursor were only detected in the nonfeeding groups. Real-time PCR was carried out to reveal changing expression patterns of two genes encoding the 10 kDa heat shock protein and diapauses hormone precursor. The data showed that this small heat shock protein gene was expressed in the normal feeding and nonfeeding larvae, but its peak signal can be detected after spending nine days in starvation (Supplementary Data 3A). As shown in supplementary data 3B, diapauses hormone precursor gene was expressed only in starved larvae at V5N and V13N.Besides these two proteins, the glutathione S-transferase 10 was identified at V5 but not found from V4N to V13N (Figure4(c)). The expressions of a PDI-like protein ERp57 and a vacuolar ATPase B subunit were upregulated at V6N and thereafter downregulated and returned to the same level as feeding silkworm (Figure 4(d)). The amounts of H+-transporting ATP synthase beta subunit isoform 2 and imaginal disc growth factor were at first lower in nonfeeding larvae but finally increased to the same level as in feeding larvae (Figures 4(d) and 4(e)). In Figure 4(f), we found that the amount of an abnormal wing disc-like protein decreased with development progress in feeding larvae, but after starvation, it kept a stable level (Figure 4(g)).


## 4. Discussion

The domesticated silkworm,Bombyx mori, is a herbivorous insect and a model organism for Lepidoptera. Chinese and Japanese scientists made great efforts to accomplish the sequencing project of the whole silkworm genome [17, 18], which offers Bombyx researchers an opportunity to identify peptides using proteomic method. Midgut, the largest digestive organ in the silkworm body, is important for the proteomic research [10, 11]. Besides cuticle, silkworm midgut is another large interface between intra and extra. It is an entry for pathogens, toxins, and pesticides. Hence, any knowledge of protein components of silkworm midgut would be essential for developing new lepidopteran pest control strategies.In this study, we used two-dimensional electrophoresis combined with MALDI-TOF-MS to analyse the midgut proteins in the normal fifth-instar silkworm larvae as well as in the starved larvae. The numbers of proteins in two groups were gradually increased from day 3 to the last larval day, suggesting that more and more midgut proteins may participate in the growth and metamorphosis of this tissue. Of the 96 proteins identified in this proteomics analysis, 69 proteins had not been observed prior to this study [10, 11, 15, 16], and the majority of midgut proteins detected showed developmental differences. For example, JHDK is responsible for the degradation of JH. Li et al. found that the mRNA of JHDK was expressed throughout 4th- and 5th-instar larvae at a constant level and distributed in foregut and midgut [19]. While JHDK protein reaches its peak at day 6 of 5th-instar larvae, it seems likely that it might be responsible for the decreased titer of JH at that stage. The fatty-acid-binding protein is involved in fatty acid transport. According to the microarray data, the gene for fatty-acid-binding protein was expressed in the 5th-instar silkworm larvae with a higher level at days 4-5, and then decreased to a constant level. The highest abundance of this protein was observed at day 6 on our 2-DE profiles, suggesting the active metabolism of lipids at that time. The abnormal wing disc-like protein detected was dramatically decreased at the day 8. It has the nucleoside diphosphate kinase activity and is involved in epithelial integrity [20]. Our observation indicated that it was regulated by hormone. Our analysis also led us to find a transferrin. Transferrin is a multifunction protein working in iron binding, transportation, cell growth, differentiation, and protection cells by inhibiting apoptosis [21].Different insects have different responses to starvation. In some insects likeManduca sexta, the subsequent development of their larvae would be delayed when starved [22], but Onthophagus  taurus larvae respond to food deprivation by shortening the length of the instar, becoming premature, pupating, and eclosing early to be small adults [23]. Bombyx mori belongs to the former type. When the 5th-instar silkworm larvae are starved from day 3, the duration of larval stage would be prolonged. It was found that the threshold weight is needed before metamorphosis in Psacothea  hilaris [24]. This may be a reason why most of starved silkworm larvae were unable to pupate.Our data showed that nine proteins undetectable in feeding larvae before metamorphosis were still detected in nonfeeding silkworm. Of these proteins, actin-depolymerizing factor 1 plays roles in regulating F-actin organization, in cell and organ expansion [25]. More detailed analysis will uncover its function in silkworm midgut. Proteomic analysis also allowed us to identify two proteins only existing in the midgut of the starved silkworm larvae. One is a small heat shock protein and the other is a diapause hormone precursor. Small heat shock proteins are known to act as molecular chaperones to refold polypeptides trapped in protein aggregates and to degrade the misfolded proteins [26]. It has been reported that small heat shock protein has additional functions such as  the regulation of programmed cell death and the promotion or inhibition of apoptosis to maintain homeostasis [27–31]. This small heat shock protein was undetectable in normal feeding larvae but can be detected in starvation ones, suggesting that it may play some roles in starved larvae for survival. Diapause hormone precursor which is relative to the synthesis of diapause hormone is another protein identified in the starved larvae. Insects can undergo diapause under extreme conditions. Gene encoding diapause hormone precursor is expressed by the induction of starvation, suggesting that, when confronts starvation, larva would stop its development against starvation for survival.Interestingly, the starved larvae produced more JHDK than the normal feeding ones. The more JHDK, the lower JH titer. At the end of penultimate instar of larva, metamorphosis is characterized by a sharp decrease in JH. However, most of the starved larvae did not undergo metamorphosis, indicating that low level of JH titer associated with high level of ecdysteroid titer in nutrition state is essential for initiating larval-to-pupal metamorphosis.Larval midgut is a main place where digestion action usually occurs due to digestive enzymes. In this study, we found that, in starvation situation, the digestive enzymes in larvae midgut were not disappeared. These digested enzymes such as glycosidase, amylase, and trehalase were synthesized, secreted, and trapped into the glycocalyx, as first proposed by Santos et al. [32].Since the chemical pesticides cause serious environmental problems, new strategies for pest control need to be developed. So far, there have been two major methods for pests control via insect midgut [33]. One is to use Cry toxins from B.  thuringiensis as biocontrol agents to lyse midgut epithelial. Another is to use RNAi technology. Selected RNAi targets in the midgut include mRNAs coding for vacuolar ATPase, ribosomal protein S4, actin, and α-tubulin [33]. Our present study provided more targets derived from midgut for the pest management.



---

*Source: 876064-2011-05-18.xml*