Differential alteration of two aminopeptidases N associated with resistance to Bacillus thuringiensis toxin Cry1Ac in cabbage looper
Kasorn Tiewsiri and Ping Wang1
Department of Ent
omology, Cornell University, New York State Agricultural Experiment Station, Geneva, NY 14456 Edited* by May R. Berenbaum, University of Illinois at Urbana–Champaign, Urbana, IL, and approved July 22, 2011 (received for review February 19, 2011)
The soil bacterium Bacillus thuringiensis (Bt) is the most successfully used biopesticide in agriculture, and its insecticidal protein genes are the primary transgenes used for insect control in transgenic crops. However, evolution of insect resistance to Bt toxins threatens the long-term future of Bt applications. To date, cases of resistance to Bt toxins have been reported in agricultural situations in six insect species, but the molecular basis for these cases of resistance remains unclear. Here we report that the resistance to the Bt toxin Cry1Ac in the cabbage looper, Trichoplusia ni, evolved in greenhouses, is associated with differential alteration of two midgut aminopeptidases N, APN1 and APN6, conferred by a trans-regulatory mechanism. Biochemical, proteomic, and molecular analyses showed that in the Cry1Ac-resistant T. ni, APN1 was signi?cantly down-regulated, whereas APN6 was signi?cantly up-regulated. The Cry1Ac resistance was correlated with downregulation of APN1 but not with the up-regulation of APN6. The concurrent up-regulation of APN6 and down-regulation of APN1 might play a role in compensating for the loss of APN1 to minimize the ?tness costs of the resistance. Along with identifying reduced expression of APN1 as the molecular basis of Bt resistance selected in an agricultural setting, our ?ndings demonstrate the importance of APN1 to the mode of action of Bt toxin Cry1Ac.
he soil bacterium Bacillus thuringiensis (Bt) has long been used as a biopesticide for insect control in agriculture and public health (1, 2). Bt genes coding for insecticidal toxins are the primary transgenes in current transgenic crops (Bt crops) conferring insect resistance (3). With the increasing scale and prolonged planting of Bt crops, evolution of insect resistance to Bt toxins in agricultural systems has become the foremost threat to the long-term future of Bt biopesticides and Bt crops. To date, cases of insect resistance to Bt toxins in open ?elds or greenhouses have been reported in six species (2). Cry toxins are the major insecticidal proteins in Bt. Once ingested by insects, the Cry protoxins are activated by the insect digestive proteases. The active toxins penetrate through the insect midgut peritrophic membrane and reach the midgut brush border membrane, where they interact with speci?c binding sites, leading to cell lysis after a multistep cascade that is incompletely understood. The Bt pathogenesis pathway is complex, and thus mechanisms of Bt resistance can be diverse. To date, numerous insect strains resistant to Bt toxins have been established by selection with Bt toxins under laboratory conditions (4), and various resistance mechanisms have been reported in laboratory-selected Bt-resistant insects, including alterations of midgut digestive proteases, decreased peritrophic membrane permeability, heightened immune response, enhanced esterase production, and reduced Cry toxin binding (1, 2, 5). More recently, mutation of an ABC transporter in Heliothis virescens was found to correlate with resistance to Cry1A toxins (6). However, it has become clear that resistance mechanisms in laboratory-selected Bt-resistant insects do not necessarily represent Bt resistance mechanisms evolved in ?eld populations (7). Molecular mechanisms of Bt resistance evolved in agriculture have not yet been reported.
Alteration of Midgut Brush Border Membrane Vesicle Proteins in Cry1Ac-Resistant Larvae Identi?ed by SDS/PAGE Separation. Proteins
from the midgut brush border membrane vesicle (BBMV) proteins of the susceptible Cornell strain and the resistant GLENCry1Ac-BCS strain exhibited highly similar protein pro?les on SDS/PAGE analysis, except for the absence of a protein band at 110 kDa in BBMVs from resistant larvae (Fig. 1). In this band, APN1 was the primary protein with 8 and 46 minor proteins, respectively, determined by nano liquid chromatography–tandem mass spectroscopy (nano LC-MS/MS) analyses with two different mass spectrometry systems, the Synapt HDMS (Waters) and the LTQ Orbitrap Velos (Thermo-Fisher Scienti?c). The relative abundance of APN1 in this 110-kDa band was estimated as 77 and 42 mol %, respectively, with the data obtained from the two mass spectrometry systems using exponentially modi?ed protein abundance index (emPAI) values of
Author contributions: K.T. and P.W. designed research; K.T. performed research; K.T. and P.W. analyzed data; and K.T. and P.W. wrote the paper. The authors declare no con?ict of interest. *This Direct Submission article had a prearranged editor. Data deposition: The DNA sequences reported in this paper have been deposited in the GenBank database (accession nos. JF303656–JF303663).
To whom correspondence should be addressed. E-mail: email@example.com or pingwang@ cornell.edu.
PNAS | August 23, 2011 | vol. 108 | no. 34 | 14037–14042
“Mode 1” type resistance is the most common type of resistance to Bt observed in insects. It is characterized by a high level of resistance (>500-fold) to at least one Cry1A toxin but not to Cry1C toxins, recessive inheritance, and reduced binding of one or more Cry1A toxins to the midgut brush border membrane (8). Mode 1 type resistance has been associated with mutations of the midgut cadherin gene in three Lepidoptera: H. virescens, Pectinophora gossypiella, and Helicoverpa armigera (9–11). The midgut cadherin has been suggested to serve as the ?rst Cry toxin receptor in the sequential interactions of Cry toxins with midgut brush border membrane proteins (12, 13). Consequently, mutations in the cadherin gene may result in Bt resistance in insects (9–11). The ?eld- and greenhouse-evolved resistance to Cry1Ac in Plutella xylostella and T. ni is of the typical mode 1 type. However, the cadherin gene in P. xylostella or T. ni is not involved in the resistance mechanism (14, 15), indicating that the mode 1 type resistance selected in agricultural systems involves a different but yet to be known molecular genetic mechanism. The Cry1Ac resistance that has evolved in T. ni populations in commercial greenhouses is monogenic and recessive in inheritance and is conferred by loss of toxin-binding sites in the larval midgut brush border membrane (16, 17). In this study, the greenhouse-derived Cry1Ac-resistant T. ni (16) was used as a unique system to identify the alteration of Cry1Ac-binding proteins and its association with the resistance to Cry1Ac. This study reports the identi?cation of the molecular basis of insect resistance to Bt toxins evolved in an agricultural system.
MW (kDa) 250 150 100 75 50 37 APN1 (77%, 42%)
(APN6) 0.03% 0.2% 1.1% 1.3% 2.1% 2.7%
could bind to multiple proteins on the membrane blot, and that the patterns for the samples from the resistant and susceptible larvae were similar, except for the lack of a positive band at 110 kDa from the resistant larvae (Fig. 2A). Identi?cation of the proteins in this 110-kDa band by nano LC-MS/MS con?rmed that APN1 was the primary protein (77 mol %). The positive band at 200 kDa was not visible on the SDS/PAGE gel stained with Coomassie blue, but analysis of the gel slice corresponding to the position of the positive band on the membrane blot showed that the band contained cadherin. The primary protein in the >200-kDa positive band was identi?ed as polycalin, which accounted for 26 mol % of the proteins in the band by emPAIbased estimation. Western blot analysis of the midgut BBMV proteins con?rmed that the 110-kDa Cry1Ac-binding protein in the midgut BBMVs was APN1, and that this 110-kDa APN1 was lacking in the resistant larvae (Fig. 2B).
Transcript Levels of APN1 and APN6 Genes Are Differentially Altered in Cry1Ac-Resistant Larvae. Quantitative RT-PCR (qRT-PCR)
Fig. 1. SDS/PAGE analysis of midgut BBMV proteins from Cry1Ac-resistant and -susceptible larvae. (A) BBMV proteins from the resistant larvae lacked a 110-kDa protein band. APN1 was identi?ed to be the primary protein in the band, accounting for 77 and 42 mol %, respectively, by two nano LC-MS/ MS systems. (B) APN6 was identi?ed in 6 of the 8 protein bands from the resistant larvae with relative abundance ranging from 0.03 to 2.7 mol % and in 1 of the 17 protein bands from the susceptible larvae with a relative abundance of 0.03 mol %. Arrows indicate the protein bands analyzed by nano LC-MS/MS for protein identi?cation.
analysis revealed a signi?cantly lower transcript level of apn1 in the resistant strain, with a ratio of 0.026 to the transcript level of apn1 in the susceptible strain. In contrast, the transcript level of apn6 in the resistant strain was signi?cantly up-regulated, with a ratio of 39 to that in the susceptible strain. The transcript levels of apn2, apn3, apn4, and apn5 did not differ signi?cantly between the susceptible and resistant larvae (Fig. 3).
APN1 and APN6 Genes Show No Genetic Linkage with Cry1Ac Resistance. For linkage analysis of apn1 and apn6 with Cry1Ac
the identi?ed proteins [protein content (mol %) = emPAI/ ∑(emPAI) × 100] (18). Further identi?cation of proteins from 17 protein bands ranging from 33 to >250 kDa from the susceptible strain and 8 matching bands ranging from 88 to >250 kDa from the resistant strain showed that APN6 was widely present in the protein bands analyzed from the resistant strain, in abundances ranging from 0.03 to 2.7 mol %. However, among the 17 protein bands analyzed from the susceptible strain, only one band minimally visible on the gel at 150 kDa was found to contain APN6 at low abundance (0.03 mol %) (Fig. 1B).
Alteration of Midgut BBMV Proteins in Cry1Ac-Resistant Larvae Identi?ed by Liquid-Based Quantitative Proteomic Analysis. Com-
parative proteomic analysis of midgut BBMV proteins from the susceptible and resistant larvae by isobaric tagging for relative and absolute quanti?cation (iTRAQ) identi?ed 1,464 and 1,440 proteins at a 95% con?dence interval in two independently prepared sample sets, including 908 proteins identi?ed in both sample sets. Of these 908 proteins, those containing at least one unique peptide (496 proteins) were used for quantitative analysis to ensure quantitative data reliability. Based on the variations observed between the two sample sets [Δlog2 ratio = │log2 ratio(sample set 1) ? log2 ratio(sample set 2)│], the iTRAQ quantitative ratio [│log2 (115 + 117)/(114 + 116)│] cutoff point was determined to be 1.5 (i.e., 95% of the proteins identi?ed fell within a deviation of 1.5). With a signi?cance cutoff threshold of 1.5, only two proteins, APN1 and APN6, were identi?ed to be signi?cantly different in quantity between the susceptible and the resistant strains in both sample sets. The log2ratio of APN1 in the BBMV proteins between the resistant and susceptible strains was found to be ?3.1 and ?3.3 in the two sample sets. Therefore, the ratio of APN1 between the resistant and susceptible strains was 0.11 (quantitative ratio = 2[log2ratio(sample set 1) + log2ratio(sample set 2)]/2). In contrast, the quantity of APN6 was higher in the BBMV proteins of the resistant strain. The log2ratio of APN6 between the resistant and susceptible strains was 2.6 in both sample sets; therefore, the ratio of APN6 between the resistant and susceptible strains was 6.0.
Midgut BBMV Proteins from Cry1Ac-Resistant Larvae Lack the 110kDa Cry1Ac- Binding Protein APN1. A Cry1Ac toxin overlay binding
resistance, a single-pair cross between a male from the GLENCry1Ac-BCS strain (apn1Gapn1G apn6Gapn6G) and a female from the Benzon strain (apn1Bapn1B apn6Bapn6B) was done to produce F1 progeny (apn1Gapn1B apn6Gapn6B), after which a female F1 progeny was backcrossed with a male from the GLENCry1Ac-BCS strain to produce backcross progenies. Genotyping of 12 Cry1Ac-selected larvae and 12 Cry1Ac-nonselected larvae from the backcross family showed that in both the Cry1Acselected and -nonselected groups, six individuals had the genotype apn1Gapn1B apn6Gapn6B, and the other six individuals had the genotype apn1Gapn1G apn6Gapn6G. The ratios of apn1Gapn1G to apn1Gapn1B and of apn6Gapn6G to apn6Gapn6B were 1:1, as would be expected with perfect random assortment (P > 0.10 by the χ2 test), in both the Cry1Ac-selected and -nonselected groups, demonstrating that apn1 and apn6 did not cosegregate with the resistance trait. This indicates that there is no genetic linkage between apn1 and apn6 with Cry1Ac resistance. Neither genotype apn1Gapn1G apn6Gapn6B nor apn1Gapn1B apn6Gapn6G was present in the 24 individuals analyzed, indicating that apn1 and
MW kDa 200 Polycalin Cadherin
MW kDa 200
analysis of the midgut BBMV proteins showed that Cry1Ac
14038 | www.pnas.org/cgi/doi/10.1073/pnas.1102555108
Fig. 2. Cry1Ac toxin overlay binding analysis and identi?cation of toxinbinding proteins from midgut BBMV proteins. (A) Three Cry1Ac-binding proteins from the midgut BBMV proteins were identi?ed as polycalin, cadherin, and APN1 by nano LC-MS/MS. The 110-kDa Cry1Ac-binding protein APN1 was lacking in the resistant strain. (B) Western blot analysis with antibodies speci?c to T. ni APN1.
Tiewsiri and Wang
Susceptible strain Resistant strain
0.01 apn1 apn2 apn3 apn4 apn5 apn6
Fig. 3. Relative transcript levels of APN genes in the Cry1Ac-resistant and -susceptible larvae by qRT-PCR analysis. Error bars indicate SEMs from the analysis of ?ve individuals. Asterisks indicate that transcript levels in the resistant larvae were signi?cantly different from the susceptible control larvae by the t test (P < 0.01; n = 5).
and APN1 was detected by Western blot analysis. The ratio between the numbers of individuals of the two groups was 8:9 in both families, which was statistically similar to the predicted random assortment ratio 1:1 (P > 0.10 by the χ2 test). In contrast, all of the larvae in the Cry1Ac-selected groups from both backcross families (22 larvae in all) had a signi?cantly reduced apn1 transcript level with a ratio <0.1 to the susceptible larvae (Fig. 4), demonstrating the linkage (cosegregation) between the signi?cantly reduced apn1 transcript level (apn1 transcript level <0.1) and Cry1Ac resistance (P < 0.001, χ2 test). APN6 gene transcript levels in both Cry1Ac-selected and -nonselected larvae were similarly up-regulated, with a ratio to that of the susceptible strain ranging from 15 to >75, with no association with resistance (Fig. 4). Although the unselected backcross individuals in both backcross families exhibited two distinct groups with differing apn1 transcript levels and differing presence or absence of APN1 (Fig. 4), these two groups of larvae had similar apn6 transcript levels (P > 0.10, t test of the means of apn6 transcript levels in the two groups). Therefore, there was no correlation between the levels of apn1 and apn6 transcripts in the larvae. Discussion In this study, comparative biochemical and proteomic analyses have shown that the Cry1Ac-resistant T. ni strain differs from its near-isogenic susceptible strain by two proteins in the midgut brush border membrane, APN1 and APN6. In the Cry1Ac-resistant strain, APN1 expression was signi?cantly reduced to a ratio of 0.11 at the protein level and a ratio of 0.026 at the transcript level to the susceptible strain. In contrast, APN6 was signi?cantly increased, to a ratio of 6.0 at the protein level and a ratio of 39 at the transcript level. Furthermore, the 110-kDa APN1 protein was lacking in the resistant larvae (Figs. 1 and 2). The differential changes of APN1 and APN6 proteins in the resistant strain were not conferred by mutations of the APN1 and APN6 genes, but instead were regulated at the transcription level by an as-yet unidenti?ed trans-regulatory mechanism (Fig. 3). The down-regulation of APN1 expression was linked to the resistance to Cry1Ac in T. ni (Fig. 4).
apn6 are in the same linkage group, as is known in other Lepidoptera (19).
Association Between Resistance and Transcription of APN Genes. Two reciprocal F2 backcross families were generated between the GLEN-Cry1Ac-BCS strain and its near-isogenic Cornell strain and treated with or without Cry1Ac selection. qRT-PCR analysis of apn1 transcript levels in larval midgut showed that in both backcross families, the non–Cry1Ac-selected individuals clearly exhibited two distinct groups in apn1 transcript level (Fig. 4). One group demonstrated a signi?cantly reduced apn1 transcript level to a ratio of no more than 0.1 to the susceptible strain and no detectable 110-kDa APN1 by Western blot analysis. The other group had an apn1 transcript level ratio of >0.4 to that of the susceptible strain, which was similar to that in F1 individuals,
Relative transcript level
Relative apn1 transcript level
1.2 1 0.8 0.6 0.4 0.2 0
Relative apn1 transcript level
1.6 1.4 1.2 1 0.8 0.6 0.4 0.2 0
backcross untreated (rs:rr)
backcross treated (rr) Relative apn6 transcript level
backcross untreated (rs:rr)
backcross treated (rr)
Relative apn6 transcript level
backcross untreated (rs:rr)
backcross treated (rr)
backcross untreated (rs:rr)
backcross treated (rr)
Fig. 4. Transcript levels of APN1 and APN6 genes in F1, Cry1Ac-selected and -nonselected larvae from backcross family a (A) and family b (B) analyzed by qRTPCR. The transcript levels are relative to the transcript levels in the susceptible Cornell strain. The expression of APN1 in nonselected larvae was also examined by Western blot analysis (Upper).
Tiewsiri and Wang
PNAS | August 23, 2011 | vol. 108 | no. 34 | 14039
Recent studies have indicated that Cry1A toxins interact with at least two types of receptors on the midgut brush border membrane (12, 13, 20, 21). Activated Cry1A toxins bind to the ?rst receptor (the midgut cadherin) with high af?nity, and the interaction with cadherin facilitates oligomerization of the toxins via a proteolytic process. The Cry1A oligomers have a high binding af?nity to the secondary receptor (APN or alkaline phosphatase) and thus bind to the secondary receptor, eventually leading to insertion of the oligomers into the midgut cell membrane, with resulting cell lysis. It also has been proposed that binding of Cry1A toxins to the cadherin may activate a cellular signaling pathway leading to cell death without the involvement of APN (21). Nevertheless, in the sequential toxin-binding events, cadherin has been suggested to serve as the ?rst receptor for Cry1A toxins. Loss-of-function (i.e., binding af?nity to Cry1A) mutations of the cadherin gene are expected to confer mode 1 type resistance as observed in H. virescens, P. gossypiella, and H. armigera (9–11); however, the mode 1 type resistance in T. ni is associated with a loss of APN1. In previous studies, we found that speci?c binding of Cry1Ab or Cry1Ac to the midgut BBMVs from the Cry1Ac-resistant T. ni larvae was undetectable, and that cadherin was not involved in the resistance (15, 16). Therefore, results from our studies on Cry1Ac resistance in T. ni indicate that cadherin alone without APN1 is not suf?cient to serve as an ef?cient binding site for Cry1Ab or Cry1Ac toxins in the midgut of T. ni, and that the toxicity of Cry1Ab or Cry1Ac in T. ni requires crucial involvement of APN1. APN1 has been well documented as a Cry1A-binding protein associated with toxicity of Cry1A toxins in lepidopterans (20). The Cry1Ac overlay binding analysis in this study con?rmed Cry1Ac binding to APN1 of T. ni (Fig. 2). More recently, APN6 has been identi?ed from lepidopteran larvae, and af?nity puri?cation of H. armigera midgut proteins using Cry1Ac as a ligand indicated that APN6 was not among the Cry1Ac-binding proteins (22). Similarly, results of the present study indicate that APN6 in the midgut BBMVs from the Cry1Ac-resistant T. ni is signi?cantly more abundant than that from the Cry1Ac-susceptible strain (Fig. 1), but no speci?c binding of Cry1Ac or Cry1Ab to the midgut BBMVs from resistant larvae could be detected (16), nor additional midgut BBMV protein bands with Cry1Ac binding af?nity were detected from resistant larvae by the Cry1Ac overlay binding assay (Fig. 2). Therefore, APN6 does not appear to be a Cry1A-binding protein. In a recent report, Pacheco et al. (23) proposed a “ping-pong” binding mechanism based on observations on the binding of Cry1Ab mutants to Manduca sexta cadherin and APN. In this ping-pong binding model, the Cry toxin ?rst binds to the abundant but low-af?nity APN, facilitating concentration of the toxin at the midgut brush border membrane for the consequent binding to cadherin. The Cry toxin binds to the APN with a high af?nity after oligomerization on interaction with cadherin (23). Our results indicate that binding of Cry1Ab or Cry1Ac to the midgut brush border membrane in T. ni requires APN1, because no speci?c binding sites for Cry1Ab or Cry1Ac could be detected in the midgut BBMVs from the Cry1Ac-resistant T. ni (16). It has been well documented that Cry1Ab or Cry1Ac can bind to isolated cadherin protein without APN (12, 24, 25). Cry1Ac also can bind to the cadherin from both the Cry1Ac-susceptible and the Cry1Ac-resistant T. ni larvae, as demonstrated by the toxin overlay binding assay (Fig. 2), in disagreement with the previous observation that Cry1Ac does not bind to the BBMVs from the resistant T. ni (16). Apparently, the binding of toxin to cadherin in the midgut brush border membrane differs from that of the separated proteins in vitro. Since the ?rst report of ?eld-evolved insect resistance to Bt in 1990 (26), populations of two insect species, P. xylostella, and T. ni, have developed resistance to Bt sprays in open-?eld or commercial greenhouses, and four other species, Busseola fusca,
14040 | www.pnas.org/cgi/doi/10.1073/pnas.1102555108
Spodoptera frugiperda, Helicoverpa zea, and P. gossypiella, have developed resistance to Bt crops (2, 3). However, the molecular mechanisms of the ?eld- and greenhouse-evolved resistance have not yet been identi?ed. Here we report the identi?cation of a molecular mechanism of Bt resistance evolved in agricultural systems. With alteration of APNs as the molecular basis for Cry1Ac resistance, the greenhouse-evolved T. ni not only is highly resistant to Cry1Ac on an arti?cial diet, but also is highly resistant to Cry1Ac plants (16). Interestingly, different alleles and differing transcript levels of APN genes have been reported in laboratory-selected Cry1A toxin-resistant strains of several lepidopterans (27, 28). However, genetic linkages between the different APN alleles or differential transcript levels and the resistance in these insects have not been established to con?rm the correlation with resistance. The heterozygous individuals and the homozygous Cry1Acresistant individuals clearly showed signi?cantly different levels of apn1 expression (unselected individuals in Fig. 4). However, the expression level of apn6 in these same individuals did not demonstrate a correlated pattern (Fig. 4). Therefore, the upregulated apn6 expression is unlikely to be a response to the reduction of apn1 expression, but may be coregulated with apn1 expression. The factor or factors regulating the differential transcription of apn1 and apn6 remain to be identi?ed. The possibility that two different but genetically linked pathways regulate the expression of apn1 and apn6 cannot be excluded. Although up-regulation of apn6 transcription is not directly linked to Cry1Ac resistance, the up-regulation of apn6 expression could play a role in minimizing possible ?tness costs associated with the resistance by compensating for the loss of APN1 with APN6. Indeed, midgut aminopeptidase activities were similar in the resistant and susceptible T. ni (16). The genetic basis for the alteration of Cry1A-binding sites in mode 1 type resistance is reportedly associated with target site (cadherin) gene mutations (9–11). Our present results show that Bt toxin-binding site alteration can be conferred by a trans-regulatory mechanism as well. Therefore, the molecular bases for mode 1 type high-level Bt resistance are not limited to binding site mutations, and the genetic bases for Bt resistance are more diverse than previously thought. Materials and Methods
Insects. A highly inbred laboratory strain of T. ni (Cornell strain) was used as the susceptible strain (17). A Cry1Ac-resistant strain of T. ni near-isogenic to the susceptible Cornell strain, GLEN-Cry1Ac-BCS, was derived from a Bt-resistant population collected from commercial greenhouses (16). The GLENCry1Ac-BCS strain had been backcrossed with the Cornell strain eight times at the time that study was carried out. In addition, a susceptible T. ni strain from Benzon Research (Benzon strain) was used for genetic linkage analysis of APN genes with resistance to Cry1Ac in T. ni. Preparation of Bt Toxin Cry1Ac. Cry1Ac protoxin was prepared from the Bt kurstaki strain HD-73, as described by Kain et al. (17). The Cry1Ac protoxin was activated by incubation with N-p-tosyl-L-phenylalanine chloromethyl ketone-treated trypsin at a ratio of 1: 20 (wt/wt) in 50 mM Na2CO3 (pH 10.0) at 37 °C for 1–16 h, and completion of digestion was examined by SDS/PAGE analysis. Activated Cry1Ac toxin was puri?ed by anion-exchange chromatography on a UNO Q column (Bio-Rad) with a linear gradient of 0–1 M NaCl in 20 mM Na2CO3 (pH 10.0). Midgut BBMV Preparation. Mid-?fth instar larvae were dissected in cold dissection buffer [17 mM Tris-HCl (pH 7.5), 5 mM EGTA, 300 mM mannitol, 1 mM PMSF] to isolate the midgut epithelium. Midgut BBMV proteins were prepared as described by Wolfersberger et al. (29). The protein concentrations of the BBMV protein preparations were determined using the Bio-Rad Protein Assay Kit II. The activities of brush border membrane marker enzymes alkaline phosphatase and aminopeptidase in the BBMV protein preparations and their initial midgut tissue homogenates were determined (30) to con?rm the enrichment of brush border membranes in the BBMV protein
Tiewsiri and Wang
preparations. The enrichment of the two enzyme activities was typically 5- to 6-fold and 7- to 10-fold, respectively. iTRAQ Labeling and Quantitative Proteomic Analysis of Midgut BBMV Proteins. BBMV proteins were solubilized in 0.5 M Hepes (pH 7.4) with 5 mM EGTA, 0.3 M mannitol, and 1% SDS and then reduced with 5 mM Tris-(2-carboxyethyl)phosphine at 37 °C for 1 h. Subsequently, the thiol groups were blocked with 8 mM methyl methanethiosulfonate at room temperature for 10 min. The protein samples were then digested with sequencing-grade modi?ed trypsin at 37 °C for 16 h, and the resulting tryptic peptides were labeled using iTRAQ reagent (Applied Biosystems). The sample from the Cornell strain was labeled separately with reporter ion reagents 114 and 116, and that from the GLEN-Cry1Ac-BCS strain was labeled with reporter ion reagents 115 and 117 as technical repeats. The labeled samples were combined and fractionated by OFFGEL isoeletric focusing electrophoresis using the Agilent 3100 OFFGEL Fractionator (Agilent) with Immobiline DryStrip (pH 3–10; 24 cm) (GE Healthcare). The fractions collected were pooled into 10 ?nal fractions and analyzed by nano LC-MS/MS analysis after desalting by solid-phase extraction using the Sep-Pak C18 cartridge (Waters). Nano LC-MS/MS analysis was performed using the LTQ Orbitrap Velos mass spectrometer equipped with nano ion source with high-energy collisional dissociation at Cornell University’s Proteomics and Mass Spectrometry Core Facility. Samples were injected onto a PepMap C18 trap column (5 μm; 300 μm × 5 mm) (Dionex) for online desalting and then separated on a PepMap C18 RP nano column (3 μm; 75 μm × 15 cm) (Dionex). The eluted peptides were detected in the LTQ Orbitrap Velos through a nano ion source with a 10-μm analyte emitter (New Objective). Data were acquired with Xcalibur 2.1 software (Thermo Fisher Scienti?c). The MS/MS raw spectra from iTRAQ experiments were processed with Proteome Discoverer 1.1 (Thermo Fisher Scienti?c), and a subsequent database search was performed using Mascot Deamon version 2.2.04 (Matrix Science) with a T. ni protein sequence database containing 15,536 sequence entries generated by combining 12,457 sequences (including 12,294 ESTs) downloaded from GenBank (http://www.ncbi.nlm.nih.gov/genbank/) on November 20, 2009, and 3,079 sequences (including 2,992 ESTs) generated in our laboratory. For the iTRAQ quantitative analysis and protein identi?cation, peptide mass tolerance and fragment mass tolerance values were 10 ppm and 30 mDa, respectively. The signi?cance threshold was set at a 95% con?dence interval, and only those peptides passing this ?lter were used for protein identi?cation. Furthermore, proteins identi?ed as containing at least two peptides with a P value < 0.001 determined by Mascot probability analysis were analyzed further. Intensities of the reporter ions (114, 115, 116, and 117) from iTRAQ tags on fragmentation were used for quanti?cation, and the relative protein ratios were normalized at the median ratio for the fourplex in each set of experiments. Internal errors for the iTRAQ analysis were examined with the two technical replicates for each BBMV protein sample included in the iTRAQ runs to ensure the quality of the analysis. The internal error was de?ned such that for 95% of the identi?ed proteins, the difference in protein ratio between the two technical replicates [│log2(117/115)│ and │log2(116/114)│] was below this value (31). For the two sample sets analyzed in this study, the internal errors of the iTRAQ analyses were 0.25 and 0.27 (log2 ratio), respectively. Similarly, biological variation was determined with the data from two independently prepared sets of midgut BBMV samples [Δlog2 ratio = │log2 ratio (sample set 1) ? log2 ratio(sample set 2)│]. The signi?cance cutoff point was set such that the variations of 95% of the identi?ed proteins between the two sample sets was below this value (31). SDS/PAGE Analysis of BBMV Proteins and Protein Identi?cation by Nano LC-MS/ MS. BBMV proteins (50 μg) were separated by 10% SDS/PAGE and then stained with Coomassie brilliant blue R250. Selected protein bands were excised and subjected to in-gel digestion with trypsin for nano LC-MS/MS analysis with the LTQ Orbitrap Velos or the Synapt HDMS MS system at Cornell University’s Proteomics and Mass Spectrometry Core Facility. MS/ MS data analysis for protein identi?cation was performed as described above.
Toxin Overlay Binding Assay and Western Blot Analysis. Toxin overlay binding assays of Cry1Ac to T. ni larval midgut BBMV proteins were conducted as described by Bravo et al. (32). The binding of Cry1Ac to the BBMV proteins on the membrane blot was detected by probing with a Cry1Ac-speci?c rabbit antibody, followed by the secondary anti-rabbit IgG antibody conjugated with alkaline phosphatase, and then visualized with a colorimetric reaction with nitroblue tetrazolium/bromochloroindolyl phosphate. Duplicate BBMV protein samples from the susceptible and resistant strains were included in the SDS/PAGE, and the gel was stained with Coomassie brilliant blue R250. The protein bands matching the positive bands from the toxin overlay binding assay were excised for protein identi?cation by nano LC-MS/ MS analysis as described above. For Western blot analysis of APN1 in T. ni larvae, proteins from midgut tissue or midgut BBMV proteins were separated by SDS/PAGE and transferred onto Immobilon-P membrane (Millipore). The membrane was treated in 5% nonfat milk and 0.5% Tween-20 and then incubated with a polyclonal antibody speci?c to T. ni APN1 at 4 °C overnight, followed by incubation with the alkaline phosphatase–conjugated secondary antibodies after several washes of the membrane with PBS. Positive antibody reactions were visualized by a colorimetric reaction as described above. qRT-PCR. Total RNA was isolated from individual midguts of ?fth instar larvae using the Qiagen RNeasy Mini Kit coupled with an on-column digestion procedure with DNase, following the manufacturer’s instructions, and then used for cDNA synthesis with the Promega ImProm-II reverse-transcription system. The cDNA preparations from individual larvae were used for qRTPCR analysis. Real-time PCR samples were prepared in iQ SYBR Green Supermix (Bio-Rad), and reactions were carried out using the iQ5 RT-PCR detection system (Bio-Rad). qRT-PCR included an initial hot start at 95 °C for 3 min, followed by 40 cycles of ampli?cation at 95 °C for 15 s, 58 °C for 30 s, and 72 °C for 30 s. Relative transcript levels of six T. ni APN genes in the test larvae compared with the control Cornell strain were calculated using the ΔΔCT method, with an actin gene transcript as an internal control for sample normalization. Genetic Linkage Analyses. For linkage analysis of APN1 and APN6 genes with Cry1Ac resistance, a backcross family between the GLEN-Cry1Ac-BCS strain and the Benzon strain was prepared by a single-pair cross of a female F1 progeny with a male from the GLEN-Cry1Ac-BCS strain. Individuals from the backcross family were reared on an arti?cial diet (non–Cry1Ac-selected) or treated with 500 μg/mL of Cry1Ac for 7 d as described by Wang et al. (16) to select Cry1Ac-resistant individuals (Cry1Ac-selected). Larval genomic DNA was prepared from individual larvae using a rapid DNA isolation method (33). A 665-bp genomic DNA fragment for an APN1 gene and a 406-bp genomic DNA fragment for an APN6 gene were ampli?ed from each larva by PCR and subsequently sequenced to determine the allele types. For linkage analysis of APN1 and APN6 gene transcript levels with Cry1Ac resistance, a single-pair cross was prepared between a male from GLENCry1Ac-BCS and a female from its near-isogenic Cornell strain to generate F1 progeny. An F1 female was backcrossed with a male from the GLEN-Cry1AcBCS strain, and an F1 male was backcrossed with a female from the GLENCry1Ac-BCS strain to generate two backcross families, known as backcross family a and backcross family b, respectively. The progenies from each backcross family were divided into two groups with and without selection with 500 μg/mL of Cry1Ac toxin as described above. Fifth instar larvae from the Cry1Ac-selected and -nonselected groups were dissected to isolate the midgut tissue, and the midgut tissue was processed for Western blot analysis of APN1 protein and qRT-PCR analysis for APN gene transcription as described above. ACKNOWLEDGMENTS. We thank Wendy Kain for excellent technical assistance, Sheng Zhang for critical advice and help with the proteomic analysis, and Xin Zhang for the antibodies used in this study. This study was supported by Agriculture and Food Research Initiative Competitive Grant 2008-35302-18806 from the US Department of Agriculture’s National Institute of Food and Agriculture.
1. Schnepf E, et al. (1998) Bacillus thuringiensis and its pesticidal crystal proteins. Microbiol Mol Biol Rev 62:775–806. 2. Bravo A, Likitvivatanavong S, Gill SS, Soberón M (2011) Bacillus thuringiensis: A story of a successful bioinsecticide. Insect Biochem Mol Biol 41:423–431. 3. Tabashnik BE, Van Rensburg JB, Carrière Y (2009) Field-evolved insect resistance to Bt crops: De?nition, theory, and data. J Econ Entomol 102:2011–2025. 4. Ferré J, Van Rie J (2002) Biochemistry and genetics of insect resistance to Bacillus thuringiensis. Annu Rev Entomol 47:501–533.
5. Heckel DG, et al. (2007) The diversity of Bt resistance genes in species of Lepidoptera. J Invertebr Pathol 95:192–197. 6. Gahan LJ, Pauchet Y, Vogel H, Heckel DG (2010) An ABC transporter mutation is correlated with insect resistance to Bacillus thuringiensis Cry1Ac toxin. PLoS Genet 6:e1001248. 7. Tabashnik BE, et al. (2003) Insect resistance to transgenic Bt crops: Lessons from the laboratory and ?eld. J Econ Entomol 96:1031–1038. 8. Tabashnik BE, et al. (1998) Insect resistance to Bacillus thuringiensis: Uniform or diverse? Philos Trans R Soc Lond B Biol Sci 353:1751–1756.
Tiewsiri and Wang
PNAS | August 23, 2011 | vol. 108 | no. 34 | 14041
9. Gahan LJ, Gould F, Heckel DG (2001) Identi?cation of a gene associated with Bt resistance in Heliothis virescens. Science 293:857–860. 10. Morin S, et al. (2003) Three cadherin alleles associated with resistance to Bacillus thuringiensis in pink bollworm. Proc Natl Acad Sci USA 100:5004–5009. 11. Xu X, Yu L, Wu Y (2005) Disruption of a cadherin gene associated with resistance to Cry1Ac δ-endotoxin of Bacillus thuringiensis in Helicoverpa armigera. Appl Environ Microbiol 71:948–954. 12. Bravo A, et al. (2004) Oligomerization triggers binding of a Bacillus thuringiensis Cry1Ab pore-forming toxin to aminopeptidase N receptor leading to insertion into membrane microdomains. Biochim Biophys Acta 1667:38–46. 13. Jurat-Fuentes JL, Adang MJ (2006) Cry toxin mode of action in susceptible and resistant Heliothis virescens larvae. J Invertebr Pathol 92:166–171. 14. Baxter SW, et al. (2005) Novel genetic basis of ?eld-evolved resistance to Bt toxins in Plutella xylostella. Insect Mol Biol 14:327–334. 15. Zhang X (2007) Sequence variation in cadherin alleles from the cabbage looper, Trichoplusia ni. MSc thesis (Cornell University, Ithaca, NY). 16. Wang P, et al. (2007) Mechanism of resistance to Bacillus thuringiensis toxin Cry1Ac in a greenhouse population of the cabbage looper, Trichoplusia ni. Appl Environ Microbiol 73:1199–1207. 17. Kain WC, et al. (2004) Inheritance of resistance to Bacillus thuringiensis Cry1Ac toxin in a greenhouse-derived strain of cabbage looper (Lepidoptera: Noctuidae). J Econ Entomol 97:2073–2078. 18. Ishihama Y, et al. (2005) Exponentially modi?ed protein abundance index (emPAI) for estimation of absolute protein amount in proteomics by the number of sequenced peptides per protein. Mol Cell Proteomics 4:1265–1272. 19. Baxter SW, Zhao JZ, Shelton AM, Vogel H, Heckel DG (2008) Genetic mapping of Bttoxin binding proteins in a Cry1A toxin-resistant strain of diamondback moth Plutella xylostella. Insect Biochem Mol Biol 38:125–135. 20. Pigott CR, Ellar DJ (2007) Role of receptors in Bacillus thuringiensis crystal toxin activity. Microbiol Mol Biol Rev 71:255–281. 21. Zhang X, Candas M, Griko NB, Taussig R, Bulla LA, Jr. (2006) A mechanism of cell death involving an adenylyl cyclase/PKA signaling pathway is induced by the Cry1Ab toxin of Bacillus thuringiensis. Proc Natl Acad Sci USA 103:9897–9902. 22. Angelucci C, et al. (2008) Diversity of aminopeptidases, derived from four lepidopteran gene duplications, and polycalins expressed in the midgut of Helicoverpa
armigera: Identi?cation of proteins binding the δ-endotoxin, Cry1Ac of Bacillus thuringiensis. Insect Biochem Mol Biol 38:685–696. Pacheco S, et al. (2009) Domain II loop 3 of Bacillus thuringiensis Cry1Ab toxin is involved in a “ping-pong” binding mechanism with Manduca sexta aminopeptidase-N and cadherin receptors. J Biol Chem 284:32750–32757. Hua G, Jurat-Fuentes JL, Adang MJ (2004) Bt-R1a extracellular cadherin repeat 12 mediates Bacillus thuringiensis Cry1Ab binding and cytotoxicity. J Biol Chem 279: 28051–28056. Flannagan RD, et al. (2005) Identi?cation, cloning and expression of a Cry1Ab cadherin receptor from European corn borer, Ostrinia nubilalis (Hubner) (Lepidoptera: Crambidae). Insect Biochem Mol Biol 35:33–40. Tabashnik BE, Cushing NL, Finson N, Johnson MW (1990) Field development of resistance to Bacillus thuringiensis in diamondback moth (Lepidoptera: Plutellidae). J Econ Entomol 83:1671–1676. Yang Y, et al. (2010) Molecular characterization and RNA interference of three midgut aminopeptidase N isozymes from Bacillus thuringiensis-susceptible and -resistant strains of sugarcane borer, Diatraea saccharalis. Insect Biochem Mol Biol 40:592–603. Zhang S, et al. (2009) Mutation of an aminopeptidase N gene is associated with Helicoverpa armigera resistance to Bacillus thuringiensis Cry1Ac toxin. Insect Biochem Mol Biol 39:421–429. Wolfersberger M, et al. (1987) Preparation and partial characterization of amino acid–transporting brush border membrane vesicles from the larval midgut of the cabbage butter?y (Pieris brassicae). Comp Biochem Physiol 86A:301–308. Jurat-Fuentes JL, Adang MJ (2004) Characterization of a Cry1Ac receptor alkaline phosphatase in susceptible and resistant Heliothis virescens larvae. Eur J Biochem 271: 3127–3135. Redding AM, Mukhopadhyay A, Joyner DC, Hazen TC, Keasling JD (2006) Study of nitrate stress in Desulfovibrio vulgaris Hildenborough using iTRAQ proteomics. Brief Funct Genomics Proteomics 5:133–143. Bravo A, Miranda R, Gómez I, Soberón M (2002) Pore formation activity of Cry1Ab toxin from Bacillus thuringiensis in an improved membrane vesicle preparation from Manduca sexta midgut cell microvilli. Biochim Biophys Acta 1562:63–69. Kikkert JR, et al. (2006) Detection of Contarinia nasturtii (Diptera: Cecidomyiidae) in New York, a new pest of cruciferous plants in the United States. J Econ Entomol 99: 1310–1315.
14042 | www.pnas.org/cgi/doi/10.1073/pnas.1102555108
Tiewsiri and Wang