- Research article
- Open Access
Expressed sequence tags from larval gut of the European corn borer (Ostrinia nubilalis): Exploring candidate genes potentially involved in Bacillus thuringiensis toxicity and resistance
BMC Genomics volume 10, Article number: 286 (2009)
Lepidoptera represents more than 160,000 insect species which include some of the most devastating pests of crops, forests, and stored products. However, the genomic information on lepidopteran insects is very limited. Only a few studies have focused on developing expressed sequence tag (EST) libraries from the guts of lepidopteran larvae. Knowledge of the genes that are expressed in the insect gut are crucial for understanding basic physiology of food digestion, their interactions with Bacillus thuringiensis (Bt) toxins, and for discovering new targets for novel toxins for use in pest management. This study analyzed the ESTs generated from the larval gut of the European corn borer (ECB, Ostrinia nubilalis), one of the most destructive pests of corn in North America and the western world. Our goals were to establish an ECB larval gut-specific EST database as a genomic resource for future research and to explore candidate genes potentially involved in insect-Bt interactions and Bt resistance in ECB.
We constructed two cDNA libraries from the guts of the fifth-instar larvae of ECB and sequenced a total of 15,000 ESTs from these libraries. A total of 12,519 ESTs (83.4%) appeared to be high quality with an average length of 656 bp. These ESTs represented 2,895 unique sequences, including 1,738 singletons and 1,157 contigs. Among the unique sequences, 62.7% encoded putative proteins that shared significant sequence similarities (E-value ≤ 10-3)with the sequences available in GenBank. Our EST analysis revealed 52 candidate genes that potentially have roles in Bt toxicity and resistance. These genes encode 18 trypsin-like proteases, 18 chymotrypsin-like proteases, 13 aminopeptidases, 2 alkaline phosphatases and 1 cadherin-like protein. Comparisons of expression profiles of 41 selected candidate genes between Cry1Ab-susceptible and resistant strains of ECB by RT-PCR showed apparently decreased expressions in 2 trypsin-like and 2 chymotrypsin-like protease genes, and 1 aminopeptidase genes in the resistant strain as compared with the susceptible strain. In contrast, the expression of 3 trypsin- like and 3 chymotrypsin-like protease genes, 2 aminopeptidase genes, and 2 alkaline phosphatase genes were increased in the resistant strain. Such differential expressions of the candidate genes may suggest their involvement in Cry1Ab resistance. Indeed, certain trypsin-like and chymotrypsin-like proteases have previously been found to activate or degrade Bt protoxins and toxins, whereas several aminopeptidases, cadherin-like proteins and alkaline phosphatases have been demonstrated to serve as Bt receptor proteins in other insect species.
We developed a relatively large EST database consisting of 12,519 high-quality sequences from a total of 15,000 cDNAs from the larval gut of ECB. To our knowledge, this database represents the largest gut-specific EST database from a lepidopteran pest. Our work provides a foundation for future research to develop an ECB gut-specific DNA microarray which can be used to analyze the global changes of gene expression in response to Bt protoxins/toxins and the genetic difference(s) between Bt- resistant and susceptible strains. Furthermore, we identified 52 candidate genes that may potentially be involved in Bt toxicity and resistance. Differential expressions of 15 out of the 41 selected candidate genes examined by RT-PCR, including 5 genes with apparently decreased expression and 10 with increased expression in Cry1Ab-resistant strain, may help us conclusively identify the candidate genes involved in Bt resistance and provide us with new insights into the mechanism of Cry1Ab resistance in ECB.
The genomic information on insects has increased tremendously during last several years. Whole genomes have been sequenced for several insect species, including the fruit fly (Drosophila melanogaster) , African malaria mosquito (Anopheles gambiae) , yellow fever mosquito (Aedes aegypti) , honey bee (Apis mellifera) , silkworm (Bombyx m ori) [5, 6], red flour beetle (Tribolium castaneum) , and 11 other Drosophila species [8, 9]. Genome sequencing of other insect species, including pea aphid (Acyrthosiphon pisum), northern house mosquito (Culex pipiens), three species of parasitoid wasp (Nasonia sp.), Hessian fly (Mayetiola destructor), blood sucking bug (Rhodnius prolixus), and body louse (Pediculus humanus), are currently in progress [10–12]. The red flour beetle is the only agricultural insect pest whose whole genome sequence has become available to date.
Lepidoptera, the second most biodiverse group of insect species after Coleoptera, represents more than 160,000 species including many of the most devastating pests of crops, forests and stored products . The silkworm was the first lepidopteran insect to have its complete genome sequenced . However, genomic information for other lepidopterans, particularly agricultural pest species is limited but urgently needed due to their economic importance and biodiversity. Sequencing of the expressed sequence tags (ESTs) has been recognized as an economical approach to identify a large number of expressed genes that can be used in gene expression and other genomic studies [14–16]. Indeed, ESTs have been generated from several lepidopteran insects including the silkworm , spruce budworm (Choristoneura fumiferana) , cotton bollworm (Helicoverpa armigera) , diamondback moth (Plutella xylostella) , tobacco hawkmoth (Manduca sexta) [21, 22], and fall armyworm (Spodoptera frugiperda) [10, 23].
It has been long recognized that the insect gut is an important target for developing new strategies for insect pest management. Until now, however, only a few studies have focused on the development of gut-specific EST libraries of lepidopterans as a tool to identify candidate genes involved in the toxicity of insecticides and the development of insecticide resistance. Gut-specific EST libraries were reported for light brown apple moth (Epiphyas postvittana) (6,416 ESTs) , bertha armyworm (Mamestra configurata) (30 serine protease-related sequences) , and European corn borer (ECB, Ostrinia nubilalis) (1,745 ESTs) .
ECB is one of the most destructive pests of corn and can cause as much as $1 billion of economic loss annually in the United States alone [27, 28]. ECB also represents a complex of stalk borers, such as the southwestern corn borer (Diatraea grandiosella) and the sugarcane borer (Diatraea saccharalis). These stalk borers share similar ecosystem and create similar damage to corn plants. Although ECB has been successfully managed using transgenic Bt corn hybrids (plants that express insecticidal toxins of Bacillus thuringiensis or Bt), there are increasing concerns about the potential development of Bt resistance in ECB because of the widespread use of Bt corn [28, 29]. Indeed, several ECB colonies have developed resistance to Bt toxins under laboratory selection conditions [30, 31].
The main target for Bt toxins is the insect midgut, where Bt protoxins are activated by gut proteases to produce activated Bt toxins. The activated toxins then bind to specific receptor(s) to confer toxicity . This means that insect resistance to Bt toxins could be conferred by protease-mediated and receptor-mediated mechanisms [33–37]. Because Bt toxins and insect gut interactions are determined by many gene products in the insect gut, including many proteins/enzymes involved in Bt protoxin activation, toxin binding to receptors and toxin degradation, any change in these systems has the potential to affect a particular Bt's specificity and efficacy, and could lead to Bt resistance in insects.
Our goals are to develop a gut-specific EST database from ECB larvae and explore candidate genes that are potentially involved in insect-Bt interactions and Bt resistance. In this paper, we report the analysis and annotations of 15,000 ESTs derived from the gut of ECB larvae. We discuss the putative identities of the ESTs, their potential biological and molecular functions, and present comparative analyses of our ESTs with sequences from other insects. This work provides the opportunity for developing an ECB gut-specific microarray that can be used to study insect-Bt interactions and genetic basis of Bt resistance in ECB. Furthermore, we revealed 52 candidate genes that could be involved in Bt toxicity and resistance. Among the 41 selected candidate genes examined by RT-PCR, we found 5 genes with apparently decreased expressions and 10 with increased expressions in Cry1Ab-resistant strain of ECB as compared with the susceptible strain of ECB. Differential expressions of these genes in a Cry1Ab-resistant strain may suggest possible involvement of these genes in Cry1Ab resistance, and therefore provides us with new insights into the mechanism of Cry1Ab resistance in ECB. This study may serve as a model for studying Bt resistance mechanisms and for developing bio-pesticides for all closely related corn stalk borers.
Results and discussion
Development and analysis of the ECB gut ESTs
We first used pPCR-XL-TOPO plasmid vector to prepare a cDNA library using total RNA purified from the whole guts of fifth-instar larvae of ECB. After we sequenced a total of 1,152 cDNA clones, we found that the cDNA inserts in the vector were not sufficient long (average length: 441 bp). Therefore, we used lambda Uni-ZAP RX vector to prepare a second cDNA library using mRNA purified from the guts of fifth-instar larvae of ECB. This library provided us with much longer cDNA inserts (average length: 674 bp). Because of this significantly improved quality of the ESTs generated from the lambda library, we used the lambda library for our further sequencing of ESTs. Among the 15,000 random cDNA clones sequenced, only <8% were from the plasmid library whereas >92% were from the lambda library (Table 1).
Our analysis of the 15,000 sequences resulted in 13,066 readable sequences (i.e., 87.1% success rate). These sequences were first trimmed for removal of vector sequences and then were subjected to filtration to exclude the sequences of <100 bp. Further analysis, using RepeatMasker and Organelle Masker programs , removed an additional 547 sequences. Thus, the total number of high quality sequences obtained was 12,519 (83.4%) with an average length of 656 bp (Table 1). These high quality sequences have been deposited in the EST database (dbEST) with GenBank accession numbers from GH987145 to GH999663 at the National Center for Biotechnology Information (NCBI). Redundancy and assembly analyses of the high quality sequences using Sequencher software (Gene Codes Corp., Ann Arbor, MI, USA) resulted in 2,895 unique ESTs, including 1,157 contiguous sequences (contigs) that consist of 2 or more sequences, and 1,738 singletons that represent single sequences. The majority of the contigs were assembled from 10 or fewer ESTs (Figure 1A). On average, however, each contig was assembled from 10.1 sequences due to a few highly redundant ESTs. Putative identities of the unique sequences were determined by searching the non-redundant database in GenBank using BLASTx. Among the 2,895 unique sequences, 1,816 (62.7%) showed significant matches at E-values of ≤ 10-3, whereas the remaining 1,077 (37.3%) did not exhibit meaningful matches (Figure 1B).
The abundance of transcripts for a particular gene of an organism can be estimated from the corresponding EST abundance in a cDNA library . The most abundant ESTs in our cDNA libraries were those encoding trypsin-like proteases and chymotrypsin-like proteases (Table 2). As this cDNA library was constructed from the gut of ECB, the high number of transcripts from the digestive enzymes was expected. The most abundant contig was #0038 which consisted of 525 ESTs, and it included 4% of the total sequences. This contig shared maximum sequence similarity with the trypsin-like protease, T25 precursor, characterized previously in ECB . Other abundant transcripts were contigs #0026 and #0062. Contig #0026 included 197 ESTs and encoded a putative chymotrypsin-like protease. Contig #0062 included 137 ESTs and encoded a putative trypsin-like serine protease. The highly expressed genes in ECB larval gut may have important implications for their growth and development. However, care must be taken in making general conclusions about the redundancy of EST's because some artifacts could also be involved .
Identification of the ORF and putative secretary proteins
The 2,895 contigs and singletons were subjected to the ORF predictor software to identify the open reading frame (ORF) of the sequences. This was done to identify the novel gene candidates, which have clear coding capacity. Among 2,895 unique ESTs, 1,119 (38.7%) had ORFs of at least 450 bp. Among 1,119 ORFs, 994 putative protein sequences (88.8%) shared sequence similarity (E-value ≤ 10-3) with known proteins in the non-redundant (NR) protein database in GenBank, whereas 125 (11.2%) did not share significant similarity with any known protein in the same database (Figure 2A). Thus, at least 11.2% of the protein-coding genes in the gut of ECB are potentially new genes. The remaining 1,553 contigs and singletons (53.6%) had an ORF of <450 bp and 223 (7.7%) did not have an ORF. Among the ESTs with ORFs of <450 bp, 452 (29.2%) had matches in the NR protein database, whereas 1,011 (70.8%) did not have matches. Many sequences did not have ORF of ≥ 450 bp because the sequences were too short (approximately 650 sequences were less than 450 bp). The lack of the ORFs in other sequences can be due to frame shift errors, 5' truncation of cDNA clones and the ESTs that were not derived from mRNA .
To identify the secretory proteins, putative protein sequences were examined to identify potential secretion signal peptide using SignalP software . A total of 439 (15.2%) putative proteins were predicted to contain signal peptides (Figure 2B). Among the putative secretory proteins, 298 sequences (67.9%) had matches with known proteins in the NR protein database, whereas 141 putative secretory proteins (32.1%) were unique, sharing no significant sequence similarity with any known protein. This information is valuable since secretory proteins are important components of biological processes in the gut [44, 45].
Comparative analyses of ECB gut ESTs
The development of EST databases has been recognized as a rapid method of sampling an organism's transcriptome and is complementary to a whole genome-sequencing project . Indeed, a large number of ESTs have been generated from other model organisms. The 2,895 contigs and singletons obtained from the larval gut of ECB were compared with the sequences from other organisms. The first hits (highest score) of the sequences in the NR database were taken into account to determine the most similar organism. The largest number of first hit sequences (390; 13.5%) came up with B. mori (Figure 3). This can be explained by the fact that the genome of B. mori has been sequenced and partially annotated, and that both ECB and B. mori are lepidopterans. The second largest number of first hit sequences (290; 10.0%) was with T. castaneum, followed by Ae. aegypti (109; 3.8%), Culex pipiens (91; 3.1%), and A. gambiae (81; 3.8%). Only 2.5% of the sequences (72) were found to be most similar to predicted protein sequences from O. nubilalis. This is simply due to the very small number of sequences currently available in NCBI database from ECB.
In order to compare our ECB gut ESTs with the 1,745 ECB ESTs that are already available in NCBI database, we performed BLASTN searches. Among our 2,895 contigs and singletons, 1,279 (44.2%) had significant matches at a cutoff E-value of ≤ 10-3 whereas 1,616 (55.8%) did not show any significant matches in NCBI database using BLASTN search. We compared our ECB ESTs with the ECB ESTs available in NCBI dbEST database. We found 475 sequences (16.4%) that had significant matches with E-values less than E-150 (Figure 4A). Within this category, 88 ESTs (3.0%) had matches with E-values less than E-150, 23 (0.8%) had E-values between E-150 and E-100, 131 (4.5%) had E-values of E-100 and E-50, 152 (5.2%) had E-values of E-50 and E-20, and 81 (2.7%) had E-values between E-20 and E-5 (Figure 4B). A total of 2,420 ESTs (83.6%) had no hits with currently available midgut ESTs in NCBI database.
Because B. mori genome has not been fully annotated, we have also compared our ESTs with all available B. mori ESTs using BLASTN. Among the 2,895 contigs and singletons, 579 (20.0%) had hits with B. mori sequences at E-value < 10-3 (Figure 4A). The remaining 2,316 ESTs (80.0%) did not show a significant match with the B. mori sequences. Among the 579 unique ESTs which had hits in the database, 43 (7.4%) had matches with E-value less than E-150, 64 (11.1%) had E-values between E-150 and E-100, 156 (26.9%) had E-values between E-100 and E-50, 135 (23.3%) had E-values between E-50 and E-20, and 181 sequences (31.3%) had E-values between E-20 and E-5 (Figure 4B).
Blast2GO software was used to obtain the gene ontology (GO) terms for the unique sequences by comparing them through the Gene Ontology Consortium . Among the 2,895 contigs and singletons, 1,815 showed blast hits at E-value ≤ 10-3 and 1,119 ESTs of the 1,815 were mapped. A total of 120 mapped ESTs showed both the GO terms and Enzyme Commission (EC) numbers. Figure 5 shows the EST functional categories, where the ECB unique ESTs were assigned to putative biological processes, molecular functions, and cellular components. Within the biological process category, 24.0% belong to cellular processes, followed by 17.0% metabolic processes, 11.0% developmental processes, 11.0% multi-cellular processes, and 8.0% each for biological regulation and localization. In the molecular function category, the maximum GO terms (40.0%) are included in catalytic activity, followed by binding (31.0%), transporter activity (10.0%), and 5.0% each for enzyme regulation activity and structural molecular activity (9.0%). In cellular components category, cell part, cell, and organelle had 27.0%, 24.0%, and 18.0% of the GO terms, respectively. They were followed by organelle part (13.0%), macromolecular complex (11.0%), envelope (4.0%), and membrane-enclosed lumen (3.0%).
Identification of ESTs potentially relevant to the Bt toxicity and resistance
The mode of Bt action in insects includes the ingestion of Bt protoxins, solubilization of Bt protoxins in insect gut, proteolytic activation of protoxins, binding of toxins to Bt receptors, membrane integration, pore formation, cell lysis, and insect death . According to this mode of action, a target insect could potentially develop resistance to Bt protoxins or toxins via one or more changes in the Bt-receptor interaction pathway. Indeed, the two most commonly identified Bt resistance mechanisms are protease-mediated and receptor-mediated resistance . Our analysis of ESTs derived from the larval gut of ECB revealed a number of genes that are potentially involved in Bt toxicity and resistance (Table 3). Specifically, we identified 18 ESTs putatively encoding trypsin-like proteases and 18 ESTs putatively encoding chymotrypsin-like proteases with E-value ranges from 2e-26 to 3e-137 and E-value 3e-27 to 3e-149, respectively. Changes in the proteolytic activity of digestive enzymes can alter the toxicity of Bt protoxins or toxins through effects on crystal solubilization and/or activation of protoxins, as well as degradation of activated toxin [33, 50–56]. A previous study from our lab has shown that Bt resistance in a Dipel-resistant strain of ECB was primarily associated with reduced trypsin-like protease activity [35, 40]. These trypsin-like proteases were also revealed in our EST analysis. Thus, our analysis of the ESTs generated from the guts of ECB larvae revealed many more candidate genes that deserve further analysis for their roles in Bt toxicity and resistance in ECB.
Our EST analysis also revealed 13 ESTs putatively encoding aminopeptidases (E-value 1e-64 to 1e-116), 1 encoding a cadherin-like protein (E-value 1e-35), and 2 encoding alkaline phosphatases (E-value 1e-115 to 1e-131). Aminopeptidase N, cadherin-like proteins, and alkaline phosphatases have been found to serve as Bt toxin binding receptors in other insect species [57–59]. To verify the function of aminopeptidase N as a receptor for Bt Cry1Ac toxin in Spodoptera litura, RNAi technology was used to reduce the expression of aminopeptidase N. This resulted in a significant reduction in the susceptibility of the insect to Cry1Ac toxin . Gahan et al.  showed that in a resistant strain (YHD2) of Heliothis virescens, there was a disruption of a cadherin-superfamily gene by a retrotransposon-mediated insertion that resulted in high levels of resistance to the Bt toxin Cry1Ac. Fernandez et al.  also reported that a GPI (glycosylphosphatidyl-inositol)-anchored ALP (alkaline phosphatase) was an important receptor molecule involved in Cry11Aa interactions with midgut cells and toxicity to Ae. aegypti larvae. These studies demonstrate that aminopeptidases, cadherin-like proteins, and alkaline phosphatases can serve as Bt toxin receptors involved in Bt toxicity and resistance. Thus, identification of these candidate Bt receptor genes in this study will allow us to further examine whether receptor-mediated resistance is involved in Bt resistance in ECB.
Comparison of expression profiles between Cry1Ab-susceptible and resistant strains of ECB
We performed RT-PCR to compare the expression patterns of the candidate genes relevant to Bt toxicity and resistance between Cry1Ab-susceptible and resistant strains of ECB. Among 41 selected genes from the 52 candidate genes, which included 15 that putatively code for trypsin-like serine proteases, 13 for chymotrypsin-like serine proteases, 10 for aminopeptidases, 2 for alkaline phosphatases, and 1 for cadherin-like protein, we found apparently decreased expressions in 2 trypsin-like and 2 chymotrypsin-like protease genes, and 1 aminopeptidase genes in the resistant strain as compared with the susceptible strain (Figure 6). Among these genes, 2 trypsin-like protease genes (contig  and ECB-30-C08) were virtually absent in the resistant strain. In contrast, we found apparently increased expressions in 3 trypsin-like and 3 chymotrypsin-like protease genes, 2 aminopeptidase genes, and 2 alkaline phosphatase genes in the resistant strain. The most noticeable increases were found in 1 trypsin-like protease (contig ), 3 chymotrypsin-like protease (ECB-V-25_E02, contig , and ECB-23_F02), 1 alkaline phosphatase (contig ), and 1 aminopeptidase (ECB-D12) genes.
Although RT-PCR is not quantitative, reproducible results of such differential expression patterns for these candidate genes in the Cry1Ab-susceptible and resistant strains of ECB may imply their potential roles in conferring or contributing to Cry1Ab resistance as well as genetic differences between the susceptible and resistant strains of ECB. Indeed, certain trypsin-like and chymotrypsin-like proteases have previously been found to activate or degrade Bt protoxins and toxins, whereas several aminopeptidases, cadherin-like proteins and alkaline phosphatases have been demonstrated to serve as Bt receptor proteins in other insect species. Thus, our results may help conclusively identify the candidate genes involved in Cry1Ab resistance and provide us with new insights into the mechanism of Cry1Ab resistance in ECB. Nevertheless, further research will be needed to confirm their involvements and to elucidate their roles in Cry1Ab resistance in ECB.
Our study resulted in a gut-specific EST database containing 12,519 high-quality ESTs from a total of 15,000 ESTs sequenced in an agriculturally important lepidopteran pest. To our knowledge, this database represents the largest gut-specific EST database from a lepidopteran pest. Our analysis using ORF predictor software showed that approximately 11.2% of the protein coding genes in our database may be specific to ECB as these sequences have an ORF of at least 450 bp but did not have significant matches with known sequences in NCBI database. We have also identified 52 candidate genes that are relevant to Bt toxicity and resistance. These genes encode trypsin-like proteases, chymotrypsin-like proteases, aminopeptidases, cadherin-like protein, and alkaline phosphatases. Furthermore, we showed differential expressions of 15 out of the 41 representative candidate genes that were examined by RT-PCR, including 5 genes with apparently decreased expressions and 10 with increased expressions in Cry1Ab-resistant strain. These results may help us further narrow down the candidate genes possibly involved in Cry1Ab resistance, and provide us with new insights into the mechanism of Bt resistance in general in ECB.
We are in the process of developing a microarray using our unique ESTs together with the ECB gut-specific sequences which are already available in the GenBank. The microarray technology will help us analyze the global change of gene expression in response to Bt protoxins/toxins. It will also allow us to analyze any genetic differences between Bt resistant and -susceptible strains of ECB. Our genomic information on ECB could also serve as a valuable resource for identifying critical/vulnerable genes from the gut of ECB that would make useful physiological targets for new toxins that could be developed for use in pest management.
Insects rearing and dissection
The KS-SC Bt-susceptible ECB colony was used for generating EST libraries. This colony originated from the egg masses collected from the cornfields near St. John, Kansas, in 1995. The colony has been reared since then on artificial diets in the laboratory at Kansas State University according to Huang et al. . The resistant ECB strain originated from a field collection of 126 diapausing larvae obtained from non-Bt hybrids in Kandiyohi Co., MN in 2001. The resistant strain was initiated from 14 larvae that survived exposure to a diagnostic Cry1Ab concentration used to identify potential changes in susceptibility to Cry1Ab [64, 65]. To minimize inbreeding or founder effects, the resistant insects were backcrossed twice with the susceptible strain which originated from the same collection. Because the resistance was incompletely recessive and involved multiple factors , the F1 progeny were randomly mated to obtain recombination of resistance factors in the F2 progeny to allow selection of resistant genotypes. The insects were then subjected to selection at a Cry1Ab concentration corresponding to two- to three-fold the LC50 for the F1 progeny (150 ng/cm2) . This selection event was designed to eliminate all the susceptible homozygotes and most of the heterozygotes. The resistant survivors from this selection event were then subjected to a second cycle of backcrossing, random mating, and selection. After six generations, the Cry1Ab concentration used in selections was gradually increased to achieve 750 ng/cm2 at generation F10, a concentration that kills virtually all F1 progeny. At generation F17, the resistance to Cry1Ab in the re-selected strain was in excess of 800-fold. The guts were dissected from fifth-instar larvae in DEPC (diethylpyrocarbonate)-treated distilled water and were stored in TRI reagent™ (Molecular Research, Inc., Cincinnati, OH) at -80°C until used.
cDNA library construction and sequencing
Total RNA was isolated from the whole guts of ECB larvae using TRI reagent™. The plasmid library was constructed using Creator SMART™ cDNA library construction kit from Clontech (Palo Alto, CA) following the manufacturer's protocols with one modification; instead of using the original phage vector, PCR fragments were cloned directly into a pPCR-XL-TOPO plasmid using a TOPO TA cloning kit (Invitrogen, Carlsbad, CA). The λ-library was constructed using ZAP-cDNA synthesis kit and ZAP-cDNA Gigapack III gold cloning kit (Stratagene, La Jolla, CA) according to the manufacturer's protocols. Briefly, double stranded cDNA was synthesized from poly(A) RNA, size-fractionated through a Sepharose CL-2B gel filtration column, and ligated into λ Uni-ZAP XR vector. The ligated DNA was packaged with the Gigapack III gold packaging extract and the library was plated on LB/agar plates. Recombinant plasmid within the lambda Uni-ZAP XR vector was in vivo excised using the ExAssist helper phage and recircularized to generate subclones in the pBluescript SK phagemid vector. To sequence the clones, M13R and M13F primers were used for 5' and 3' sequencing, respectively. Plasmid DNA was isolated using Qiagen Bio Robot 3000 and sequenced using an ABI 3700 DNA analyzer.
EST analyses and annotations
The DNA sequences were preprocessed by using the online software EGassembler . Specifically, sequence cleaning process was employed to trim the vector and adaptor sequences from the ESTs. RepeatMasker process was used to mask the interspersed repeats and low complexity regions of the sequences by using Drosophila Repbase repeat library. The sequences were further masked by using vector masking against NCBI's vector library and organelle masking against mitochondrial library. The preprocessed ESTs were then assembled by using Sequencher software (Gene Codes Corp., Ann Arbor, MI). The ORF regions of the assembled ESTs were identified by using the ORF predictor software  and secretory proteins were identified by looking for signal peptide sequence using SignalP software . Gene ontology (GO) annotation was derived using Blast2GO software http://www.blast2go.de/.
Comparative analysis of ESTs
The ECB unique ESTs were comparatively analyzed for their sequence similarities against other organisms. The organism associated with the EST showing the highest BLAST score in GenBank databases was selected. The ECB gut ESTs were also compared with sequences from the silkworm and ECB that are currently available in the database by using BLASTN with a cutoff E-value of 10-3.
Expression profiling by RT-PCR
Forty-one out of the 52 candidate genes were selected for comparing their apparent gene expression profiles between the Cry1Ab-susceptible and resistant strains of ECB by using RT-PCR. These genes were selected solely based on their representations among different gene groups from our EST analysis. After total RNA was isolated from four midguts dissected from one-day-old fifth-instar larvae of each strain (Cry1Ab-susceptible and resistant strains) of ECB by using TRI reagentTM (Sigma, St. Louis, MO), it was treated with TURBO™ DNase (Ambion, Austin, TX)to remove any genomic DNA contaminations. Three micrograms of total RNA was used for synthesis of first strand cDNA using SuperScript® III First-Strand Synthesis System (Invitrogen, Carlsbad, CA). cDNA prepared from total RNA was used as a template for RT-PCR. A minimum of two biological replications was used for all the PCR primer pairs. For all trypsin-like (except for ECB-30_C08) and chymotrypsin-like serine protease, alkaline phosphatase, and RPS3 genes, 25 PCR cycles were used whereas for aminopeptidase and cadherin-like protein, 27 PCR cycles were used. For one trypsin-like serine protease gene (ECB-30_C08), however, 33 PCR cycles were used as the expression of this gene using fewer cycles was not visible on agarose gels. Each PCR was performed for above mentioned number of cycles, each consisting of 94°C for 30s, 55°C for 60s, and 72°C for 60s. The sequences of forward and reverse PCR primers, and expected size of PCR product for each of 41 candidate genes are provided in Additional file 1.
Adams MD, Celniker SE, Holt RA, Evans CA, Gocayne JD, Amanatides PG, Scherer SE, Li PW, Hoskins RA, Galle RF, et al: The genome sequence of Drosophila melanogaster. Science. 2000, 287: 2185-2195. 10.1126/science.287.5461.2185.
Holt RA, Subramanian GM, Halpern A, Sutton GG, Charlab R, Nusskern DR, Wincker P, Clark AG, Ribeiro JM, Wides R, et al: The genome sequence of the malaria mosquito Anopheles gambiae. Science. 2002, 298 (5591): 129-149. 10.1126/science.1076181.
Nene V, Wortman JR, Lawson D, Haas B, Kodira C, Tu ZJ, Loftus B, Xi Z, Megy K, Grabherr M, et al: Genome sequence of Aedes aegypti, a major arbovirus vector. Science. 2007, 316: 1718-1723. 10.1126/science.1138878.
The Honeybee Genome Sequencing Consortium: Insights into social insects from the genome of the honeybee Apis mellifera. Nature. 2006, 443: 931-949. 10.1038/nature05260.
Mita K, Kasahara M, Sasaki S, Nagayasu Y, Yamada T, Kanamori H, Namiki N, Kitagawa M, Yamashita H, Yasukochi Y, et al: The genome sequence of silkworm, Bombyx mori. DNA Res. 2004, 11: 27-35. 10.1093/dnares/11.1.27.
Xia Q, Zhou Z, Lu C, Cheng D, Dai F, Li B, Zhao P, Zha X, Cheng T, Chai C, et al: A draft sequence for the genome of the domesticated silkworm (Bombyx mori). Science. 2004, 306: 1937-1940. 10.1126/science.1102210.
Tribolium Genome Sequencing Consortium: The genome of the model beetle and pest Tribolium castaneum. Nature. 2008, 452: 949-955. 10.1038/nature06784.
Crosby MA, Goodman JL, Strelets VB, Zhang P, Gelbart WM, the FlyBase Consortium: FlyBase: genomes by the dozen. Nucleic Acids Res. 2007, 35: D486-D491. 10.1093/nar/gkl827.
Lin MF, Carlson JW, Crosby MA, Matthews BB, Yu C, Park S, Wan KH, Schroeder AJ, Gramates LS, St Pierre SE, et al: Revisiting the protein-coding gene catalog of Drosophila melanogaster using 12 fly genomes. Genome Res. 2007, 17: 1823-1836. 10.1101/gr.6679507.
Deng Y, Dong Y, Thodima V, Clem RJ, Passarelli AL: Analysis and functional annotation of expressed sequence tags from the fall armyworm Spodoptera frugiperda. BMC Genomics. 2006, 7: 264-10.1186/1471-2164-7-264.
Grimmelikhuijzen CJ, Cazzamali G, Williamson M, Hauser F: The promise of insect genomics. Pest Manag Sci. 2007, 63: 413-416. 10.1002/ps.1352.
Sattelle DB, Jones AK, Buckingham SD: Insect genomes: challenges and opportunities for neuroscience. Invert Neurosci. 2007, 7: 133-136. 10.1007/s10158-007-0054-2.
Pierce NE: Predatory and parasitic Lepidoptera: carnivores living on plants. J Lepido Soc. 1995, 49: 412-453.
Gerhold D, Caskey CT: It's the genes! EST access to human genome content. BioEssays. 1996, 18: 973-981. 10.1002/bies.950181207.
Dimopoulos G, Casavant TL, Chang S, Scheetz T, Roberts C, Donohue M, Schultz J, Benes V, Bork P, Ansorge W: Anopheles gambiae pilot gene discovery project: Identification of mosquito innate immunity genes from expressed sequence tags generated from immune-competent cell lines. Proc Natl Acad Sci. 2000, 97: 6619-6624. 10.1073/pnas.97.12.6619.
Porcel BM, Tran AN, Tammi M, Nyarady Z, Rydaker M, Urmenyi TP, Rondinelli E, Pettersson U, Andersson B, Aslund L: Gene survey of the pathogenic protozoan Trypanosoma cruzi. Genome Res. 2000, 10: 1103-1107. 10.1101/gr.10.8.1103.
Mita K, Morimyo M, Okano K, Koike Y, Nohata J, Kawasaki H, Kadono-Okuda K, Yamamoto K, Suzuki MG, Shimada T, Goldsmith MR, Maeda S: The construction of an EST database for Bombyx mori and its application. Proc Natl Acad Soc USA. 2003, 100: 14121-14126. 10.1073/pnas.2234984100.
Li L, Krell PJ, Arif BM, Feng Q, Doucet D: Integration and Analysis of an EST database from the insect Choristoneura fumiferana. 2003, [http://www.pestgenomics.org/database.htm]
Dong D-J, He H-J, Chai LQ, Jiang XJ, Wang JX, Zhao XF: Identification of genes differentially expressed during larval molting and metamorphosis of Helicoverpa armigera. BMC Dev Biol. 2007, 7: 73-10.1186/1471-213X-7-73.
Eum J, Kang S, Han S: Annotated expressed sequence tags for studies on the expression pattern of the immunized Plutella xylostella. [http://www.uniprot.org/uniprot/Q60FS0]
Robertson H, Martos MR, Sears CR, Todres EZ, Walden KKO, Nardi JB: Diversity of odourant binding proteins revealed by an expressed sequence tag project on male Manduca sexta moth antennae. Insect Mol Biol. 1999, 8: 501-518. 10.1046/j.1365-2583.1999.00146.x.
Zou Z, Najar F, Wang Y, Roe B, Jiang H: Pyrosequence analysis of expressed sequence tags for Manduca sexta hemolymph proteins involved in immune responses. Insect Biochem Molec Biol. 2008, 38: 677-682. 10.1016/j.ibmb.2008.03.009.
Nègre V, Hôtelier T, Volkoff AN, Gimenez S, Cousserans F, Mita K, Sabau X, Rocher J, Lopez-Ferber M, D'Alençon E, Audant P, Sabourault C, Bidegainberry V, Hilliou F, Fournier P: SPODOBASE: an EST database for the lepidopteran crop pest Spodoptera. BMC Bioinformatics. 2006, 7: 322-10.1186/1471-2105-7-322.
Simpson R, Newcomb RD, Gatehouse HS, Crowhurst RN, Chagné D, Gatehouse LN, Markwick NP, Beuning LL, Murray C, Marshall SD, et al: Expressed sequence tags from the midgut of Epiphyas postvittana (Walker). (Lepidoptera: Tortricidae). Insect Mol Biol. 2007, 16: 675-690.
Hegedus D, Baldwin D, O'Grady M, Braun L, Gleddie S, Sharpe A, Lydiate D, Erlandson M: Midgut proteases from Mamestra configurata (Lepidoptera: Noctuidae) larvae: characterization, cDNA cloning, and expressed sequence tag analysis. Arch Insect Biochem Physiol. 2003, 53: 30-47. 10.1002/arch.10084.
Coates BS, Sumerford DV, Hellmich RL, Lewis LC: Mining an Ostrinia nubilalis midgut expressed sequence tag (EST) library for candidate genes and single nucleotide polymorphisms (SNPs). Insect Molec Biol. 2008, 17: 607-620. 10.1111/j.1365-2583.2008.00833.x.
Ostlie KR, Hutchison WD, Hellmich RL: Bt corn and European corn borer: long-term success through resistance management. 1997, NCR Publication 602. University of Minnesota, St. Paul, MN
Gould F: Sustainability of transgenic insecticidal cultivars: integrating pest genetics and ecology. Annu Rev Entomol. 1998, 43: 701-726. 10.1146/annurev.ento.43.1.701.
Wolfenbarger LL, Phifer PR: Biotechnology and ecology – The ecological risks and benefits of genetically engineered plants. Science. 2000, 290: 208-2093. 10.1126/science.290.5499.2088.
Huang F, Higgins RA, Buschman LL: Baseline susceptibility and changes in susceptibility to Bacillus thuringiensis subsp. kurstaki under selection pressure in European corn borer (Lepidoptera: Pyralidae). J Econ Entomol. 1997, 90: 1137-43.
Bolin PC, Hutchison WD, Andow DA: Long-term selection for resistance to Bacillus thuringiensis Cry1Ac endotoxin in a Minnesota population of European corn borer (Lepidoptera: Crambidae). J Econ Entomol. 1999, 92: 1021-1030.
Gill SS, Cowles EA, Pietrantonio PV: The mode of action of Bacillus thuringiensis δ-endotoxin. Annu Rev Entomol. 1992, 37: 615-636. 10.1146/annurev.en.37.010192.003151.
Oppert B, Kramer KJ, Beeman RW, Johnson DE, McGaughey WH: Proteinase-mediated insect resistance to Bacillus thuringiensis toxins. J Biol Chem. 1997, 272: 23473-23476. 10.1074/jbc.272.38.23473.
Huang F, Zhu KY, Buschman LL, Higgins RA, Oppert B: Comparison of midgut proteinases in Bacillus thuringiensis susceptible and -resistant European corn borer, Ostrinia nubilalis (Lepidoptera: Pyralidae). Pestic Biochem Physiol. 1999, 65: 132-139. 10.1006/pest.1999.2438.
Li H, Oppert B, Higgins RA, Huang F, Zhu KY, Buschman LL: Comparative analysis of proteinase activities of Bacillus thuringiensis-resistant and -susceptible Ostrinia nubilalis (Lepidoptera: Crambidae). Insect Biochem Mol Biol. 2004, 34: 753-762.
Lee MK, Rajamohan F, Gould F, Dean DH: Resistance to Bacillus thuringiensis CryIA ±-endotoxins in a laboratory-selected Heliothis virescens strain is related to receptor alteration. Appl Environ Microbiol. 1995, 61: 3836-42.
Herrero S, Oppert B, Ferre J: Different mechanisms of resistance to Bacillus thuringiensis toxins in the Indianmeal moth. Appl Environ Microbiol. 2001, 67: 1085-89. 10.1128/AEM.67.3.1085-1089.2001.
Masoudi-Nejad A, Tonomura K, Kawashima S, Moriya Y, Suzuki M, Itoh M, Kanehisa M, Endo T, Goto S: EGassembler: online bioinformatics service for large-scale processing, clustering and assembling ESTs and genomic DNA fragments. Nucleic Acids Res. 2006, 34: W459-462. 10.1093/nar/gkl066.
Papanicolaou A, Joron M, McMillan WO, Blaxter ML, Jiggins CD: Genomic tools and cDNA derived markers for butterflies. Mol Ecol. 2005, 14: 2883-2897. 10.1111/j.1365-294X.2005.02609.x.
Li H, Oppert B, Higgins RA, Huang F, Buschman LL, Gao JR, Zhu KY: Characterizations of cDNAs encoding three trypsin-like proteinases and mRNA quantitative analyses in Bt-resistant and -susceptible Ostrinia nubilalis (Lepidoptera: Crambidae). Insect Biochem Mol Biol. 2005, 35: 845-860.
Kikuchi T, Aikawa T, Kosaka H, Pritchard L, Ogura N, Jones JT: Expressed sequence tag (EST) analysis of the pine wood nematode Bursaphelenchus xylophilus and B. mucronatus. Mol Biochem Parasitol. 2007, 155: 9-17. 10.1016/j.molbiopara.2007.05.002.
Whitfield CW, Band MR, Bonaldo MF, Kumar CG, Liu L, Pardinas JR, Robertson HM, Soares MB, Robinson GE: Annotated expressed sequence tags and cDNA microarrays for studies of brain and behavior in the honey bee. Genome Res. 2002, 12: 555-566. 10.1101/gr.5302.
Bendtsen JD, Nielsen H, von Heijne G, Brunak S: Improved prediction of signal peptides: SignalP 3.0. J Mol Biol. 2004, 340: 783-795. 10.1016/j.jmb.2004.05.028.
O'Donnell RA, Blackman MJ: The role of malaria merozoite proteases in red blood cell invasion. Curr Opinion Microbiol. 2005, 8: 422-427. 10.1016/j.mib.2005.06.018.
Adams MD, Soares MB, Kerlavage AR, Fields C, Venter JC: Rapid cDNA sequencing (expressed sequence tags) from a directionally cloned human infant brain cDNA library. Nat Genet. 1993, 4: 373-380. 10.1038/ng0893-373.
Tosini F, Trasarti E, Pozio E: Apicomplexa genes involved in the host cell invasion: the Cpa135 protein family. Parassitologia. 2006, 48: 105-107.
Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, et al: Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 2000, 25: 25-29. 10.1038/75556.
Bravo A, Gill SS, Soberón M: Mode of action of Bacillus thuringiensis Cry and Cyt toxins and their potential for insect control.
Ferré J, Van RJ: Biochemistry and genetics of insect resistance to Bacillus thuringiensis. Annu Rev Entomol. 2002, 47: 501-533. 10.1146/annurev.ento.47.091201.145234.
Milne R, Kaplan H: Purification and characterization of a trypsin-like digestive enzyme from spruce budworm (Choristoneura fumiferana) responsible for the activation of d-endotoxin from Bacillus thuringiensis. Insect Biochem Mol Biol. 1993, 23: 663-673. 10.1016/0965-1748(93)90040-Y.
Oppert B, Kramer KJ, Johnson DE, MacIntosh SC, McGaughey WH: Altered protoxin activation by midgut enzymes from a Bacillus thuringiensis resistant strain of Plodia interpunctella. Biochem Biophys Res Commun. 1994, 198: 940-947. 10.1006/bbrc.1994.1134.
Oppert B, Kramer KJ, Johnson DE, Upton SJ, McGaughey WH: Luminal proteinases from Plodia interpunctella and the hydrolysis of Bacillus thuringiensis Cry1Ac protoxin. Insect Biochem Mol Biol. 1996, 26: 571-583. 10.1016/S0965-1748(96)00013-6.
Martínez-Ramírez AC, Real MD: Proteolytic processing of Bacillus thuringiensis CryIIIA toxin and specific binding to brush border membrane vesicles of Leptinotarsa decemlineata (Colorado potato beetle). Pestic Biochem Physiol. 1996, 54: 115-122. 10.1006/pest.1996.0015.
Keller M, Sneh B, Strizhov N, Prudovsky E, Regev A, Koncz C, Schell J, Zilberstein A: Digestion of δ-endotoxin by gut proteases may explain reduced sensitivity of advanced instar larvae of Spodoptera littoralis to CryIC. Insect Biochem Mol Biol. 1996, 26: 365-373. 10.1016/0965-1748(95)00102-6.
Forcada C, Alcacer E, Garcera MD, Martinez R: Differences in the midgut proteolytic activity of two Heliothis virescens strains, one susceptible and one resistant to Bacillus thuringiensis. Arch Insect Biochem Physiol. 1996, 31: 257-272. 10.1002/(SICI)1520-6327(1996)31:3<257::AID-ARCH2>3.0.CO;2-V.
Forcada C, Alcacer E, Garcera MD, Tato A, Martinez R: Resistance to Bacillus thuringiensis Cry1Ac toxin in three strains of Heliothis virescens: proteolytic and SEM study of the larval midgut. Arch Insect Biochem Physiol. 1999, 42: 51-63. 10.1002/(SICI)1520-6327(199909)42:1<51::AID-ARCH6>3.0.CO;2-6.
Herrero S, Gechev T, Bakker PL, Moar WJ, Maagd RA: Bacillus thuringiensis Cry1Ca-resistant Spodoptera exigua lacks expression of one of four aminopeptidase N genes. BMC Genomics. 2005, 6: 96-10.1186/1471-2164-6-96.
Hara H, Atsumi S, Yaoi K, Nakanishi K, Higurashi S, Miura N, Tabunoki H, Sato R: A cadherin-like protein functions as a receptor for Bacillus thuringiensis Cry1Aa and Cry1Ac toxins on midgut epithelial cells of Bombyx mori larvae. FEBS Lett. 2003, 538: 29-34. 10.1016/S0014-5793(03)00117-0.
Jurat-Fuentes JL, Adang MJ: Characterization of a Cry1Acreceptor alkaline phosphatase in susceptible and resistant Heliothis virescens larvae. Eur J Biochem. 2004, 271: 3127-3135. 10.1111/j.1432-1033.2004.04238.x.
Rajagopal R, Sivakumar S, Agrawal N, Malhotra P, Bhatnagar RK: Silencing of midgut aminopeptidase N of Spodoptera litura by double-stranded RNA establishes its role as Bacillus thuringiensis toxin receptor. J Biol Chem. 2002, 277: 46849-46851. 10.1074/jbc.C200523200.
Gahan LJ, Gould F, Heckel DG: Identification of a gene associated with Bt resistance in Heliothis virescens. Science. 2001, 293: 857-860. 10.1126/science.1060949.
Fernandez LE, Aimanova KG, Gill SS, Bravo A, Soberón M: A GPI-anchored alkaline phosphatase is a functional midgut receptor of Cry11Aa toxin in Aedes aegypti larvae. Biochem J. 2006, 394: 77-84. 10.1042/BJ20051517.
Huang F, Higgins RA, Buschman LL: Baseline susceptibility and changes in susceptibility to Bacillus thuringiensis subsp. kurstaki under selection pressure in European corn borer, Ostrinia nubilalis Hübner (Lepidoptera: Pyralidae). J Econ Entomol. 1997, 90: 1137-1143.
Siegfried BD, Spencer T, Crespo ALB, Storer NP, Head GP, Owens ED, Guyer D: Ten years of monitoring for Bt resistance in the European corn borer: What we know, what we don't know and what we can do better. Amer Entomol. 2007, 53: 208-214.
Crespo ALB, Spencer T, Alves AP, Hellmich RL, Blankenship EE, Magalhaes LC, Siegfried BD: On-plant survival and inheritance of resistance to Cry1Ab toxin from Bacillus thuringiensis in a field-derived strain of European corn borer, Ostrinia nubilalis. Pest Manag Sci. 2009,
Marçon PCRG, Young LJ, Steffey KL, Siegfried BD: Baseline susceptibility of European corn borer (Lepidoptera: Crambidae) to Bacillus thuringiensis toxins. J Econ Entomo. 1999, 92: 279-285.
Min XJ, Butler G, Storms R, Tsang A: OrfPredictor: predicting protein-coding regions in EST-derived sequences. Nucleic Acids Res. 2005, 33: W677-680. 10.1093/nar/gki394.
Conesa A, Gotz S, Garcia-Gomez JM, Terol J, Talon M, Robles M: Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005, 21: 3674-3676. 10.1093/bioinformatics/bti610.
The authors thank Yoonseong Park for his helpful comments on an earlier draft of this manuscript, Lisa Tan for maintaining the European corn borer colonies, Xiang Liu for his technical assistance. Mention of trade names or commercial products in this publication is solely for the purpose of providing specific information and does not imply recommendation or endorsement by U.S. Department of Agriculture or Kansas State University. This study was supported in part by the Kansas Agricultural Experiment Station and the Arthropod Genomics Center funded by K-State Targeted Excellence program at Kansas State University. This paper is contribution No 09-105-J from the Kansas Agricultural Experiment Station. The Ostrinia nubilalis voucher specimens (voucher No. 079) are located in the Kansas State University Museum of Entomological and Prairie Arthropod Research, Manhattan, Kansas.
CK conducted the major part of this study including experimental design, construction of the cDNA libraries, EST analysis, RT-PCR analysis, and manuscript preparation. YCZ participated in experimental design, EST sequencing and preliminary analysis of EST data. MSC assisted in the development of the project, the establishment of the collaboration in EST sequencing, and manuscript preparation. LLB participated in experimental design, maintenance of the insect culture, and manuscript preparation. RAH participated in the development of the project and experimental design. JY assisted in EST sequencing and analysis. BDS and ALBC contributed materials and participated in data analysis and manuscript preparation. SM participated in experimental design and manuscript preparation. KYZ coordinated the project and participated in experimental design, EST analysis, and manuscript preparation. All authors read and approved the final manuscript.