Bacillus thuringiensis (Bt) produces pore forming toxins that have been used for pest control in agriculture for many years. However, their molecular and cellular mode of action is still unclear. While a first model – referred to as the pore forming model – is the most widely accepted scenario, a second model proposed that toxins could trigger an Mg2+-dependent intracellular signalling pathway leading to cell death. Although Cry1Ca has been shown to form ionic pores in the plasma membrane leading to cell swelling and death, we investigated the existence of other cellular or molecular events involved in Cry1Ca toxicity. The Sf9 insect cell line, derived from Spodoptera frugiperda, is highly and specifically sensitive to Cry1Ca. Through a selection program we developed various levels of laboratory-evolved Cry1Ca-resistant Sf9 cell lines. Using a specific S. frugiperda microarray we performed a comparative transcriptomic analysis between sensitive and resistant cells and revealed genes differentially expressed in resistant cells and related to cation-dependent signalling pathways. Ion chelators protected sensitive cells from Cry1Ca toxicity suggesting the necessity of both Ca2+ and/or Mg2+ for toxin action. Selected cells were highly resistant to Cry1Ca while toxin binding onto their plasma membrane was not affected. This suggested a resistance mechanism different from the classical ‘loss of toxin binding’. We observed a correlation between Cry1Ca cytotoxicity and the increase of intracellular cAMP levels. Indeed, Sf9 sensitive cells produced high levels of cAMP upon toxin stimulation, while Sf9 resistant cells were unable to increase their intracellular cAMP. Together, these results provide new information about the mechanism of Cry1Ca toxicity and clues to potential resistance factors yet to discover.

Bacillus thuringiensis (Bt) is a Gram-positive bacterium that produces proteins with a wide variety of insecticidal properties. These microbial insecticides have been used for decades as pest control agents and they represent an alternative to chemical pesticides in a modern agriculture that strives to be more respectful to the environment and to human health. Moreover, observations of insect resistance to classical chemical pesticides favoured the development and use of the insecticidal weapons produced by Bt (Chattopadhyay and Banerjee, 2018).

The major insecticidal weapons of Bt are two multigenic families of toxins, cry and cyt (Crickmore et al., 1998). Cry proteins are produced as protoxins in crystal inclusions during Bt sporulation phase. They belong to the pore forming toxins (PFT) class of bacterial toxins (Palma et al., 2014). After spore and crystal ingestion they are delivered to the insect intestinal tract where their activation occurs allowing binding to midgut epithelial cells that results in cell lysis and death of the target insect (Raymond et al., 2010). Two different modes of action on intestinal cells have been proposed and particularly well documented for Cry1A toxins. The first and well-established model, referred to as the pore-forming model, requires the sequential binding to two specific receptors localized at the plasma membrane of insect intestinal cells: a cadherin receptor protein (CADR) and a glycosyl-phosphatidylinositol (GPI) membrane-anchored aminopeptidase N (APN). This sequential binding allows pre-pore complex formation and membrane insertion where they act as functional cationic-specific pores causing osmolytic lysis of targeted cells (Jiménez-Juárez et al., 2007; Soberón et al., 2000; Zhuang et al., 2002). The second model of Cry action, completely independent of pore formation, is referred to as the signal transduction model. Zhang and colleagues showed that an Mg2+-dependent signalling pathway is essential to Cry1A-induced cell death. This model also starts with the binding of Cry1A to the primary receptor CADR triggering the recruitment and activation of a heterotrimeric G protein, activation of an adenylyl cyclase (AC), and elevation of intracellular cyclic AMP (cAMPi). This second messenger then activates a protein kinase A (PKA) whose activity is shown to be important for toxin-induced cell death (Zhang et al., 2005, 2006). If CADR and APN were the first proteins identified as Cry receptors in insects, numerous other molecules that specifically bind Cry toxins, such as alkaline phosphatase or ABC transporter have been reported (Heckel, 2012; Pigott and Ellar, 2007). The existence of these many potential receptors makes it more difficult to demonstrate a single mode of action of Cry toxins.

Despite all the studies published on Cry1A toxins, numerous events are still missing in the scenario of toxin action leading to insect cell death (Vachon et al., 2012). Cry1C has been described as a pore forming toxin able to oligomerize and form ionic channels after membrane insertion (Laflamme et al., 2008; Peyronnet et al., 2001). Previous studies using histological sections or purified plasma membranes of insect epithelial midgut cells revealed specific Cry1C receptors with low or no competition with Cry1A toxins (Agrawal et al., 2002; Alcantara et al., 2004; Kwa et al., 1998). Cry1C and Cry1A toxins specifically bind to distinct isoforms of APN present in the brush border membrane of Manduca sexta insect (Luo et al., 1996; Masson et al., 1995). Moreover, Liu and colleagues have shown that resistance of the diamondback moth to Cry1C was not the result of reduced binding of this toxin to insect midgut membranes, i.e. a resistance mechanism different from that observed for Cry1A-resistant insect (Liu et al., 2000). Finally, Cry1C has been shown to be effective against Cry1A-resistant insects or to act synergistically with Cry1A on target insects (Abdullah et al., 2009; Xue et al., 2005). Pyramiding of the cry1A and cry1C genes in transgenic plants has been shown to be an effective strategy for delaying the evolution of insect (Tang et al., 2017; Cao et al., 2002). However, field population of diamondback moth repeatedly selected in laboratory for their resistance to Cry1C presented various levels of cross-resistance with members of the Cry1A toxin sub-class (Liu et al., 2001). Taken together, all these data are revealing the complexity of Cry toxin action and are suggesting that Cry1A and Cry1C toxins may have different mechanisms of toxicity. Knowledge of the molecular events involved in Cry1C mode of action would allow a better use of this toxin as a potential alternative to Cry1A and the understanding of resistance mechanisms developed by insects.

The fall armyworm, Spodoptera frugiperda (J. E. Smith), is one of the most serious pests ravaging economically important crops across the American continent and was recently spotted in Africa (Goergen et al., 2016; Nagoshi and Meagher, 2008). Transgenic crops expressing Bt toxins have been widely and efficiently used against S. frugiperda. However, numerous cases of laboratory- or field-evolved resistance to Cry toxin have been reported (Farias et al., 2014; Storer et al., 2010). In order to mimic strategies developed by S. frugiperda to reduce its susceptibility to toxins we engineered a resistant cellular model using Sf9 cells. These cells are derived from S. frugiperda and have been intensively used as an insect model to study Cry1C toxicity. Among all the Lepidoptera cell lines studied by Kwa and colleagues, Sf9 cells appeared highly and specifically sensitive to Cry1C with more than 98% of the cells affected (swollen cells) after 3 h of toxin treatment (Kwa et al., 1998). However, a cell cycle-dependent resistance to Cry1C was observed in Sf9 cells. Cells arrested in G2-M phase were insensitive to Cry1C and gradually regained their sensitivity when the arresting agent was removed (Avisar et al., 2005).

In the present study we selected two levels of resistance to Cry1Ca in Sf9 cells. Genes differentially expressed between resistant selected cell lines and the wild-type sensitive one were studied. Our analysis led us to focus on the molecular mechanism involving cation-dependent cellular pathway that could be involved in Cry1Ca toxicity.

Development and characterization of Sf9 cells resistant to Cry1Ca toxin

Sf9 cells are known to be susceptible to the Cry1C Bt toxin (Kwa et al., 1998). We first verified if the Sf9 cell line we were using was still highly and specifically sensitive to the Cry1Ca toxin. Viability tests were performed onto Sf9 wild-type (Sf9wt) cells in the presence of increasing concentrations of different Bt toxins (Fig. 1A). Solubility problems prevented us from calculating LD50 of Cry1Ab on Sf9 cells. Nevertheless, as observed in the inset in Fig. 1 and as reported by others (Bretschneider et al., 2016) Cry1Ab (white triangles) was not efficient in killing cells, while Cry1Ca (black circles) efficiently killed the cells with a lethal dose, giving 50% dead cells (LD50) evaluated at 0.9±0.15 mg/l. In order to study the mode of action of Cry1Ca and to investigate resistance mechanisms observed in S. frugiperda, we started a selection program to develop cell strains resistant to this toxin. Doses calculated for Sf9wt cells and corresponding to the LD50 and LD80 were used to select cells over 40 generations. At the end of the selection program, resistant cell lines obtained using LD50 or LD80 doses were named Sf9-LD50 and Sf9-LD80 respectively. Cell viability was tested at the 40th generation and LD50 measured for Sf9-LD50 (white circles in Fig. 1A) and Sf9-LD80 resistant cells (black triangles in Fig. 1A) were 4.4±1.6 mg/l and 53±13 mg/l, respectively. We calculated resistance factors (RF) using LD50Sf9-resistant/LD50Sf9wt ratio: RF=4.9±0.8 for Sf9-LD50 cells and RF=58.9±2.8 for Sf9-LD80 cells. These resistant cells presented no significant time growth difference as compared to Sf9wt cells used for the selection program (data not shown). Cells were incubated for 18 h with 0.9 mg/l of Cry1Ca and observed using phase contrast microscopy (Fig. 1B). Both untreated Sf9 resistant cells presented no significant morphological difference compared to Sf9 sensitive cells without toxin in the media. When Cry1Ca was added in the media, Sf9wt cells showed morphological changes with vacuole formation and lysis. While a large number of Sf9wt cells were killed, only a few Sf9-LD50 cells were affected and no change in Sf9-LD80 cell density and morphology was observed.

Fig. 1.

Sensitivity of Sf9wt cells to Cry toxins and development of populations resistant to Cry1Ca. (A) Cells were incubated in complete media containing various concentrations of Bt toxins. Cell viability was evaluated by SYBR Gold staining. Inset shows the sensitivity of Sf9wt cells to Cry1Ca toxin (●) as compared to Cry1Ab toxin (). Main graph shows cell viability measurements of selected resistant Sf9 cell lines Sf9-LD50 (○) and Sf9-LD80 (▾) compared to unselected Sf9wt (●) in the presence of Cry1C toxin. Data are presented as the mean±s.d. of three independent experiments performed in duplicate. Cry1Ca lethal concentration giving 50% dead cells (LD50) is deduced from the curves. (B) Sensitive or resistant Sf9 cell lines were incubated for 18 h in toxin-free media (-Tox) or containing 0.9 mg/l of Cry1Ca toxin (+Tox) and observed using phase contrast microscopy. Shown results are those from one experiment representative of three independent experiments. Black arrows point to examples of affected cells. Scale bar: 20 µm.

Fig. 1.

Sensitivity of Sf9wt cells to Cry toxins and development of populations resistant to Cry1Ca. (A) Cells were incubated in complete media containing various concentrations of Bt toxins. Cell viability was evaluated by SYBR Gold staining. Inset shows the sensitivity of Sf9wt cells to Cry1Ca toxin (●) as compared to Cry1Ab toxin (). Main graph shows cell viability measurements of selected resistant Sf9 cell lines Sf9-LD50 (○) and Sf9-LD80 (▾) compared to unselected Sf9wt (●) in the presence of Cry1C toxin. Data are presented as the mean±s.d. of three independent experiments performed in duplicate. Cry1Ca lethal concentration giving 50% dead cells (LD50) is deduced from the curves. (B) Sensitive or resistant Sf9 cell lines were incubated for 18 h in toxin-free media (-Tox) or containing 0.9 mg/l of Cry1Ca toxin (+Tox) and observed using phase contrast microscopy. Shown results are those from one experiment representative of three independent experiments. Black arrows point to examples of affected cells. Scale bar: 20 µm.

Cry1Ca still binds to the plasma membrane of resistant Sf9 cell lines

Among all known resistance mechanisms to Bacillus entomopathogenic toxins, so far loss of plasma membrane receptor binding appears to be one of the major mechanisms developed by insects (Heckel et al., 2007). In order to verify if our resistant cell lines established such a resistant mechanism, sensitive or resistant cells were incubated for 40 min with complete media containing 3 mg/l of Cry1Ca and immunostaining of the toxin was performed (Fig. 2). We found no detectable difference in Cry1Ca membrane binding between sensitive and resistant cells. When compared to the control where no toxin was present in the media (left row), sensitive and resistant cell lines showed comparable numbers of cells stained by antibodies to Cry1Ca. Most of the plasma membranes were uniformly labelled with 2few punctuates eventually caused by toxin clustering and/or oligomerization. The same experiment was performed with Drosophila Schneider 2 (S2) cells known to be insensitive to Cry1Ca toxin (bottom row). We found no Cry1Ca signal associated with the plasma membrane of S2 cells supporting the idea that toxin membrane binding is necessary for toxin action. Our data suggested that the observed resistance in Sf9-LD50 and Sf9-LD80 cells is probably not caused by loss of a plasma membrane receptor.

Fig. 2.

Cry1Ca binding to the plasma membrane of Sf9 and S2 cell lines,Sf9wt cells or resistant cell lines Sf9-LD50 or Sf9-LD80 were incubated 40 minin complete media free of toxin (−Cry1C) or containing 3 mg/l of Cry1Ca toxin (+Cry1C). Immunostaining was performed using primary antibodies to Cry1Ca and secondary antibodies coupled to Alexa Fluor 488. Toxin localization was visualized using an inverted epifluorescence microscope. As a control, immunostaining was carried out onto cells incubated without toxin (−Cry1C) and S2 cells (S2+Cry1C) insensitive to Cry1Ca toxin. The numbers indicate the percentage of cells immunostained. Shown results are those from one experiment representative of two independent experiments. Scale bar: 20 µm.

Fig. 2.

Cry1Ca binding to the plasma membrane of Sf9 and S2 cell lines,Sf9wt cells or resistant cell lines Sf9-LD50 or Sf9-LD80 were incubated 40 minin complete media free of toxin (−Cry1C) or containing 3 mg/l of Cry1Ca toxin (+Cry1C). Immunostaining was performed using primary antibodies to Cry1Ca and secondary antibodies coupled to Alexa Fluor 488. Toxin localization was visualized using an inverted epifluorescence microscope. As a control, immunostaining was carried out onto cells incubated without toxin (−Cry1C) and S2 cells (S2+Cry1C) insensitive to Cry1Ca toxin. The numbers indicate the percentage of cells immunostained. Shown results are those from one experiment representative of two independent experiments. Scale bar: 20 µm.

Comparative transcriptomic analysis of Sf9 resistant versus sensitive cell lines

Comparative transcriptomic analysis was performed onto total mRNA isolated from sensitive Sf9 cells and resistant selected cells. Expression levels were analysed on a specific S. frugiperda microarray representing 9721 genes. To evaluate the potential effects of cell aging and culture condition on gene expression, we compared the expression patterns of the 8th (Sf9p8) and 40th (Sf9p40) generations of wild-type cells. Only three showed genes-coding for a Rho GTPase, a fatty acid-binding protein 2 (FABP2) and an unidentified protein among the 9721 tested genes that were significantly differently transcribed (Table S1, sheet 1). The variations observed for the Rho GTPase and the FABP2 are in accordance with previous studies showing that mRNA and protein levels of these two molecules are regulated throughout the cell cycle and cell aging, and are involved in the regulation of key components of the cell cycle (Coleman et al., 2004; Pu et al., 1999; Yoon et al., 2007). The low gene expression variation observed between Sf9p8 and Sf9p40 validated our global approach to identify genes with expression profiles specifically modified during the establishment of Cry1Ca resistance in Sf9 cells. Classification of the differentially expressed genes using Gene Ontology terms revealed numerous biological processes modified in both resistant selected cells: transcription, translation, transport, binding, enzymatic activity, secretion, signalling and structural (Fig. 3A). Among all these biological processes, we observed that as RF increases, an enrichment of terms related to the cell signalling pathway (from 4% in Sf9-LD50 to 8% in Sf9-LD80) and to transcription processes (from 6% in Sf9-LD50 to 14% in Sf9-LD80) are seen. Moreover, as shown in the Venn diagram proposed in Fig. 3B, while 82 and 125 genes presented variation of expression specifically in Sf9-LD50 and in Sf9-LD80 cells respectively, only 49 genes appeared to be common to both selected cell lines. Thirty-six out of these 49 common genes differentially up- or downregulated in the same way in both resistant cell lines and were identified (Table 1). While numerous cellular processes were modified, we thought that these common genes could play a central role in the mechanism of resistance to Cry1Ca or in its mode of action. Among these 36 identified common genes, 17% appeared to be possibly involved in cellular pathways requiring or modifying cation intracellular concentrations. Moreover, Omic analysis revealed that cation homeostasis and/or cation-dependent cell signalling processes were modified in cells or insects resistant to Bt toxins (Candas et al., 2003; Canton et al., 2015; Nanoth Vellichirammal et al., 2015). All these observations led us to investigate the importance of cations in Cry1Ca toxicity using Sf9-LD80 only; the selected cell lines with the highest calculated RF.

Fig. 3.

Transcriptomic analysis of resistant Sf9 cell lines. (A) Pie charts of the biological functions of detected genes over or under-expressed in Sf9 resistant cell lines. Sf9-LD50 and Sf9-LD80 are compared to Sf9wt sensitive cells. Genes identified on the Spodo Database are classified into seven categories using Gene Ontology terms: transport, binding, enzymatic activity, secretion, signalling, structural, transcription and translation. (B) Venn diagram showing number of genes differentially expressed in Sf9-LD50 only (82), in Sf9-LD80 only (125) and common in both resistant cell lines (49).

Fig. 3.

Transcriptomic analysis of resistant Sf9 cell lines. (A) Pie charts of the biological functions of detected genes over or under-expressed in Sf9 resistant cell lines. Sf9-LD50 and Sf9-LD80 are compared to Sf9wt sensitive cells. Genes identified on the Spodo Database are classified into seven categories using Gene Ontology terms: transport, binding, enzymatic activity, secretion, signalling, structural, transcription and translation. (B) Venn diagram showing number of genes differentially expressed in Sf9-LD50 only (82), in Sf9-LD80 only (125) and common in both resistant cell lines (49).

Table 1.

List of genes with the same expression profile in both resistant cell lines

List of genes with the same expression profile in both resistant cell lines
List of genes with the same expression profile in both resistant cell lines

Protective effect of divalent ions chelators

Divalent ions have been described as common key players in the main proposed modes of action of Bt toxins (Zhang et al., 2006). To investigate their contribution to Cry1Ca toxicity, we carried out Sf9 cells intoxication during 18 h with different concentrations of toxin in the presence of 5 mM EDTA or 5 mM EGTA and we evaluated cell viability. EDTA chelates Mg2+ and Ca2+, whereas EGTA preferentially chelates Ca2+. No detectable mortality was observed when Sf9wt or Sf9-LD80 cells were incubated with chelators only (data not shown). As shown in Fig. 4A, both EDTA and EGTA displayed a protective effect regarding toxin action when compared to Sf9 cells intoxicated without ion chelators in the media. The LD50s measured for Sf9wt (black circles), Sf9wt with EDTA (black triangles) and Sf9wt with EGTA (white circles) were 2.1±0.4, 35.8±6 and 18±5 mg/l of Cry1Ca, respectively. Calculation of LD50s ratios revealed an 8.57±0.76-fold decrease in Cry1Ca toxicity when EGTA was present in the media and a 17.05±0.43-fold decrease for EDTA. Observation under a light microscope revealed a slight effect on the morphology of Sf9wt protected cells (Fig. 4B) and Sf9-LD80 resistant cells (data not shown). Most of the Sf9wt cells protected by chelator presence (+Tox +EDTA or +Tox +EGTA) were still intact compared to cells with Cry1Ca only (+Tox), but they appeared less flattened and round shaped than untreated Sf9 (−Tox). These minor morphological changes were not observed when both cell lines were incubated with media containing chelators only (data not shown) suggesting that these changes were not the result of cell adhesion alteration caused by extracellular ions clustering. In order to verify whether the observed protective effect of both EDTA and EGTA was not the result of Cry1Ca binding loss onto Sf9 plasma membrane, we performed immunostaining using Cry1Ca primary antibodies in presence of 5 mM chelators (Fig. 5). Control cells (CT) were incubated with EDTA only. The percentage of membrane stained cells was calculated in the different conditions. While 0% of total cells were membrane-labelled in the control, we did not observe significant difference in toxin binding onto Sf9wt membranes in presence or absence of EDTA or EGTA. Once again, Sf9-LD80 cells presented comparable results to Sf9wt cells with about one-third of total cells positive for toxin binding. All these results revealed a significant protective effect of cationic chelators suggesting an important contribution of both Mg2+ and Ca2+ in Cry1Ca action. This protective effect was not likely due to loss of interaction of Cry1Ca with plasma membrane binding regions.

Fig. 4.

Sf9wt cells are protected from Cry1Ca toxicity by ion chelators. (A) Sf9wt cells were incubated for 18 h in complete media containing various concentrations of Cry1Ca toxin alone, Cry1Ca toxin +5 mM EGTA or Cry1Ca toxin +5 mM EDTA. Cell viability was evaluated by SYBR Gold staining. Data are presented as the mean±s.d. of three independent experiments done in duplicate. (B) Sf9wt cells were observed using phase contrast microscopy 18 h after incubation with 5 mg/l of Cry1Ca toxin (+Tox), Cry1Ca toxin +5 mM EDTA (+Tox +EDTA) or Cry1Ca toxin +5 mM EGTA (+Tox +EGTA). Cell morphology was compared to Sf9wt cells in complete media only (−Tox). Shown results are those from one experiment representative of three independent experiments. Scale bar: 20 µm.

Fig. 4.

Sf9wt cells are protected from Cry1Ca toxicity by ion chelators. (A) Sf9wt cells were incubated for 18 h in complete media containing various concentrations of Cry1Ca toxin alone, Cry1Ca toxin +5 mM EGTA or Cry1Ca toxin +5 mM EDTA. Cell viability was evaluated by SYBR Gold staining. Data are presented as the mean±s.d. of three independent experiments done in duplicate. (B) Sf9wt cells were observed using phase contrast microscopy 18 h after incubation with 5 mg/l of Cry1Ca toxin (+Tox), Cry1Ca toxin +5 mM EDTA (+Tox +EDTA) or Cry1Ca toxin +5 mM EGTA (+Tox +EGTA). Cell morphology was compared to Sf9wt cells in complete media only (−Tox). Shown results are those from one experiment representative of three independent experiments. Scale bar: 20 µm.

Fig. 5.

Cry1Ca plasma membrane binding is not affected by the presence of chelators. Sf9wt cells (left row) or Sf9-LD80 cells (right row) were incubated for 40 min in complete media containing 3 mg/l of Cry1Ca toxin alone (+Tox), Cry1Ca toxin +5 mM EDTA (+Tox +EDTA) or Cry1Ca toxin +5 mM EGTA (+Tox +EGTA). Immunostaining was performed using primary antibodies to Cry1Ca and secondary antibodies coupled to Alexa Fluor 488. Toxin localization was visualized using an inverted epifluorescence microscope. As a control (CT), immunostaining was carried out onto cells incubated without toxin. The numbers indicate the percentage of cells immunostained. Shown results are those from one experiment representative of two independent experiments. Scale bar: 20 µm.

Fig. 5.

Cry1Ca plasma membrane binding is not affected by the presence of chelators. Sf9wt cells (left row) or Sf9-LD80 cells (right row) were incubated for 40 min in complete media containing 3 mg/l of Cry1Ca toxin alone (+Tox), Cry1Ca toxin +5 mM EDTA (+Tox +EDTA) or Cry1Ca toxin +5 mM EGTA (+Tox +EGTA). Immunostaining was performed using primary antibodies to Cry1Ca and secondary antibodies coupled to Alexa Fluor 488. Toxin localization was visualized using an inverted epifluorescence microscope. As a control (CT), immunostaining was carried out onto cells incubated without toxin. The numbers indicate the percentage of cells immunostained. Shown results are those from one experiment representative of two independent experiments. Scale bar: 20 µm.

Cry1Ca toxin induces Cyclic AMP production

Previous studies proposed that Cry protein toxicity involves activation of an Mg2+-dependent signalling pathway leading to the production of the second messenger cAMP (Zhang et al., 2006). Since (i) divalent ions chelators appeared to have a protective effect upon Cry1Ca toxicity and (ii) numerous crosstalk have been described between cAMP pathway and calcium signal (Hofer, 2012), we incubated Sf9 cells with the toxin in media with or without chelators and we evaluated cAMP production (Fig. 6). Activation with forskolin, a direct adenylate cyclase (AC) activator, was used as a positive control for cAMP production. As expected, Sf9wt cells treated with 100 µM forskolin presented a 35-fold increase of cAMP level when compared to untreated cells. A significantly stronger effect was observed in Sf9wt cells treated with Cry1Ca, with a 96-fold increase. When stimulations were achieved in presence of EDTA in the media, we observed a specific inhibitory effect on cAMP production affecting forskolin treatment while cAMP production triggered by Cry1Ca was only slightly affected. Indeed, the fold increase in cAMP using forskolin or Cry1Ca were respectively reduced to 18-fold and 92-fold. Finally, when EGTA was present in the media of stimulated cells we observed a specific inhibitory effect onto cAMP production due to Cry1Ca toxin with a cAMP fold increase significantly reduced from 96-fold to 35-fold. We observed no EGTA inhibitory effect on cAMP levels obtained after forskolin treatment. Taken together these results showed that Cry1Ca was able to trigger the production of cAMP in Sf9 cells through a mechanism requiring extracellular Ca2+.

Fig. 6.

cAMP is produced in Sf9wt cells after toxin treatment. Sf9wt cells were incubated 20 min in complete media only (−) or media supplemented with (+) 100 µM forskolin (Fors.), 1 mg/l Cry1Ca toxin (Tox.), 5 mM EDTA or 5 mM EGTA. After incubation, cytosolic HCl extract were performed and cAMP amount measured using an immunoassay kit. cAMP fold increase was evaluated compared to non-stimulated cells (see Materials and Methods). Data are presented as the mean±s.e. of three independent experiments done in duplicate. Student's t-tests were used for comparisons among groups. ***P<0.001; n.s., not significant compared to the condition with forskolin only. “P<0.05, “““P<0.001, compared to the condition with toxin only.

Fig. 6.

cAMP is produced in Sf9wt cells after toxin treatment. Sf9wt cells were incubated 20 min in complete media only (−) or media supplemented with (+) 100 µM forskolin (Fors.), 1 mg/l Cry1Ca toxin (Tox.), 5 mM EDTA or 5 mM EGTA. After incubation, cytosolic HCl extract were performed and cAMP amount measured using an immunoassay kit. cAMP fold increase was evaluated compared to non-stimulated cells (see Materials and Methods). Data are presented as the mean±s.e. of three independent experiments done in duplicate. Student's t-tests were used for comparisons among groups. ***P<0.001; n.s., not significant compared to the condition with forskolin only. “P<0.05, “““P<0.001, compared to the condition with toxin only.

Sf9-LD80 cells are impaired in cAMP production

Since Cry1Ca action seemed to require increase of intracellular concentration of cAMP ([cAMP]i) we wondered whether the resistant cell line Sf9-LD80 was able to produce this second messenger. Sf9wt sensitive cells or Sf9-LD80 resistant cells were stimulated with forskolin or Cry1Ca and the increase in cAMP was measured (Fig. 7). In contrast to what we observed for Sf9wt cells, Sf9-LD80 had an altered response to stimulation by both forskolin and the toxin. While 100 µM forskolin increased the cAMP amount eightfold (less than one quarter of what we observed for Sf9wt cells), no cAMP increase was detected after stimulation with Cry1Ca. This experiment suggests that Sf9-LD80 cells could overcome Cry1Ca toxicity by blocking cAMP-dependent signalling normally triggered by the toxin in sensitive cells.

Fig. 7.

Sf9-DL80 cells are affected in cAMP production. Sf9wt or Sf9-LD80 cells were incubated for 20 min in complete media only (CT) or media supplemented with 100 µM forskolin (+Fors.) or 1 mg/l Cry1Ca toxin (+Tox.). After incubation, cytosolic HCl extracts were performed and cAMP amount measured using an immunoassay kit. cAMP fold increase was evaluated compared to non-stimulated cells (see Materials and Methods). Data are presented as the mean±s.e. of three independent experiments done in duplicate. Student's t-tests were used for comparisons among groups. *P<0.05, compared to Sf9wt cells stimulated with forskolin only. “P<0.05, compared to Sf9-LD80 stimulated with toxin only.

Fig. 7.

Sf9-DL80 cells are affected in cAMP production. Sf9wt or Sf9-LD80 cells were incubated for 20 min in complete media only (CT) or media supplemented with 100 µM forskolin (+Fors.) or 1 mg/l Cry1Ca toxin (+Tox.). After incubation, cytosolic HCl extracts were performed and cAMP amount measured using an immunoassay kit. cAMP fold increase was evaluated compared to non-stimulated cells (see Materials and Methods). Data are presented as the mean±s.e. of three independent experiments done in duplicate. Student's t-tests were used for comparisons among groups. *P<0.05, compared to Sf9wt cells stimulated with forskolin only. “P<0.05, compared to Sf9-LD80 stimulated with toxin only.

In this report, we describe Sf9 cells resistant to the Cry1Ca toxin obtained after a selection program. We tested Sf9wt cells for their sensitivity to another Cry toxin (Cry1Ab) at toxin concentrations up to 100 mg/l in the same time conditions as with Cry1Ca and confirmed that they were almost ineffective on these cells, suggesting a peculiar mode of action for Cry1Ca. The LD50 measured with these cells that are specifically sensitive to Cry1Ca were in the same concentration range to what was observed in previous studies (Kwa et al., 1998). Therefore, we started our selection program using LD50 or LD80 Cry1Ca doses and we obtained cell lines with an increasing resistance factor, named Sf9-LD50 and Sf9-LD80.

While the so-called ‘mode 1 resistance’, the most frequently mechanism of Cry toxins resistance, results in toxin receptor binding defects (Gahan et al., 2001; Jurat-Fuentes et al., 2011; Morin et al., 2003) we observed no detectable difference in Cry1Ca binding to the plasma membrane of Sf9 sensitive or resistant cells (Fig. 2). It is important to note that this approach using immuno-fluorescence microscopy could not rule out an undetectable small diminution of receptor affinity or loss of a secondary receptor, which could be sufficient to perturb toxin toxicity onto Sf9-LD80 cells. However, our binding results are strongly correlated with our transcriptomic analysis. Indeed, we did not detect any significant variation of expression of previously described and widely accepted plasma membrane toxin receptor as CADR or APN. All these results strongly suggest that Cry1Ca resistance observed in the Sf9-LD50 or Sf9-LD80 cells is probably not caused by loss or alteration of plasma membrane binding.

Use of differential transcriptomic analysis to understand resistance mechanisms developed by insects to counteract toxicity is a challenging task. It generates a large amount of data and shows that cellular responses are likely diverse and affect many cellular processes as observed in Fig. 3A and in other studies (Nanoth Vellichirammal et al., 2015; Canton et al., 2015). Therefore, we decided to focus on genes differentially expressed in the same way in both selected resistant cells, as they may be key players in resistance processes (Table 1). Hence, we highlighted genes involved in cation-dependent signalling pathways or cation homeostasis. For example, calmodulin-binding transcription activator (CAMTA) is activated after a rise in intracellular calcium concentration ([Ca2+]i). Both Ca2+ and calmodulin (Cam) bind to CAMTA, which translocates into the nucleus and activates its target genes (Finkler et al., 2007). Like CAMTA, the predicted Calcineurin B homologue was differentially expressed in both selected resistant cell lines. Calcineurin is a calcium-activated protein phosphatase that participates in diverse biological functions through interaction with various partners (Li et al., 2011). It is composed of two subunits, a catalytic subunit (CnA) and a regulatory subunit (CnB) containing Ca2+ binding sites. Increased intracellular Ca2+ allows binding of the cations to the regulatory sites, which initiates a series of conformational changes and the activation of calcineurin (Yang and Klee, 2000). Thus, in both examples a rise in [Ca2+]i appears as a crucial event in the activation of our potential key players in Cry1Ca resistance.

Earlier studies showed the involvement of divalent cations in Cry activity. Cry1Ca toxicity increases with the presence of calcium in a dose-dependent manner and the toxin induces a rapid elevation of intracellular calcium in treated cells (Monette et al., 1997). We followed intoxication of Sf9 sensitive cells in the presence of EDTA or EGTA and observed a protective effect for both ion chelators (Fig. 4). Similar to the observations for Cry1Ab binding onto S5 cells (Zhang et al., 2005), neither EDTA nor EGTA affected the interaction of Cry1Ca with the plasma membrane of Sf9wt or Sf9-LD80 cell lines (Fig. 5). The protection obtained with the Ca2+ chelator EGTA was not the result of a loss of toxin binding capacity to Sf9wt cells and could be the result of the failure of the cells to increase their [Ca2+]i in response to the toxin. If pore formation or the signal transduction pathway are the two main models of toxin action leading to cell death, a third model was proposed in which formation of toxin oligomers followed by pore formation due to oligomer insertion in cell membrane could participate in establishment of a signalling pathway (Jurat-Fuentes and Adang, 2006). Pores composed of Cry1Ca toxin do not induce only potassium-selective channels but increase the permeability of the plasma membrane of Sf9 cells to numerous solutes (Guihard et al., 2000; Villalon et al., 1998). The toxin has been shown to trigger an intracellular calcium surge in Sf9 cells and is still toxic when Sf9 cells are incubated in media containing barium, a divalent cation used as potassium channel blocker (Monette et al., 1997; Schwartz et al., 1991). So, the protective effect of calcium chelators could result from blocking Ca2+ entry in the cell through poorly selective channels obtained by toxin membrane insertion and the incapacity of this calcium surge to contribute to a Ca2+-dependent signal pathway, leading to cell death. When we used the Mg2+ chelator EDTA, we observed a stronger decrease of Cry1Ca toxicity compared to protection obtained with EGTA (Fig. 4). Like Zhang and colleagues, who observed that Cry1Ab still binds to the plasma membrane of High five cells in the presence of EDTA (Zhang et al., 2005), we saw no binding loss of Cry1Ca to Sf9wt sensitive or Sf9-LD80 resistant cells when incubated with 5 mM EDTA (Fig. 5). So, as described for Cry1Ab, depletion of Mg2+ using EDTA could prevent Cry1Ca-induced cellular response and stop the establishment of a Mg2+-dependent intracellular pathway critical to promote cell death (Zhang et al., 2006). Since Mg2+ is essential for adenylyl cyclase (AC) to catalyse the formation of cAMP (Tesmer et al., 1999; Pieroni et al., 1995) and since AC activity has been shown by Zhang et al. (2005) to participate in Cry1Ab toxicity, we monitored cAMP production after Sf9 cell intoxication with Cry1Ca (Fig. 6). Treatment of Sf9 sensitive cells with forskolin, used as positive control, produced a 38-fold increase of intracellular [cAMP] while Cry1Ca presented an even stronger stimulation of cAMP production. [cAMP]i increase after toxin treatment was not affected by the presence of EDTA, whereas the cAMP increase was reduced two fold in Sf9wt cells stimulated with forskolin plus EDTA as compared to forskolin alone. In contrast, EGTA affected only Cry1Ca-induced cAMP production. These results suggest the importance of Ca2+ ions in toxin action via formation of cAMP. Since numerous AC activities have be shown to be regulated by [Ca2+]i (Willoughby and Cooper, 2007), Cry1Ca, after formation of pores in the Sf9 membrane, could activate AC through entry in cells of extracellular Ca2+. Näsman's work showed that applying 1 mM EGTA to Sf9 cell media blocked the raise in [Ca2+]i and drastically reduced the increase of cAMP levels after stimulation with octopamine (Näsman et al., 2002). Thus, the increase in cAMP levels observed in our cells treated with Cry1Ca could be due to the activation of calcium-dependent AC, with EGTA protecting cells by preventing the increase in calcium and AC activation.

If EGTA was the most efficient to reduce cAMP elevation in response to Cry1Ca stimulation, EDTA which preferentially binds Mg2+ appeared the most protective against Cry1Ca toxicity. This effect of EDTA could be explained, on one hand, by the necessity of Mg2+ for catalytic activity of ACs. Indeed, Mg2+ which binds to the active domain of ACs is essential for their activity and therefore for the synthesis of cAMP (Pieroni et al., 1995; Tesmer et al., 1999). On the other hand, Mg2+ could have a more central role than Ca2+ in participating in the activation of numerous proposed other partner proteins involved in toxin-induced signalling pathway, beginning with heterotrimeric G protein activation (Sprang, 1997). However, while in the absence of Mg2+ we observed an effect on cAMP production by forskolin stimulation, EDTA did not seem to prevent the increase of cAMP after the action of the toxin. Cry1Ca might be triggering the production of cAMP by activation of an Mg2+-independent AC whose regulation requires Ca2+ (Braun et al., 1977; Litvin et al., 2003).

It has been proposed that the mechanism of action of the toxin can be independent of cAMP production (Knowles and Farndale, 1988; Portugal et al., 2017). Portugal and colleagues clearly demonstrated that both Cry1Ab and Cry1Ac toxins triggered cell death in CF1 cells via pore formation activity. While they saw neither the activation of PKA nor an increase in intracellular cAMP concentration during Cry1Ab or Cry1Ac intoxication, in our hands Cry1Ca triggered an increase of the [cAMP]i in favour of a peculiar mode of action of this toxin (Fig. 6). Moreover, we observed (Fig. 7) that the Sf9-LD80 resistant cell line was unable to elevate [cAMP]i in response to Cry1Ca treatment. Thus, the [cAMP]i increase seems to be an essential element to the toxicity induced by Cry1Ca but appeared insufficient since forskolin, although it may increase cAMP levels, presents no toxicity onto Sf9 cells. While Zhang and colleagues observed a similar lack of forskolin toxicity, Cry1Ab cytotoxic effects were potentialised when sensitive cells were pre-treated with forskolin (Zhang et al., 2006). In light of these observations, they proposed that if cAMP is a key player in toxin effect through activation of PKA, specific PKA-dependent effector(s) mediating downstream cell death activity might be stimulated by Cry1Ab but not by forskolin. All these observations raise the question of how specificity is maintained in this cAMP/PKA system, activated by both forskolin and toxins (Cry1Ab and Cry1Ca) but where only Bt toxins are capable of generating cellular events leading to cell death. Since we identified Ca2+ as another important molecule involved in the Cry1Ca mode of action, we might think that through toxin pore formation, Ca2+ entry in the cell would stimulates Ca2+-dependent events interacting with the PKA pathway. The protein kinase C (PKC) pathway could be such Ca2+-dependent event which, together with the PKA activation, would allow the toxin to perform its cytotoxic effect. PKC is a large superfamily of protein kinases comprising numerous members activated by elevation of intracellular Ca2+ concentration (Webb et al., 2000). Coordinated activity and cross-talk between the PKA and PKC pathways have been shown in multiple cellular processes (Robinson-White and Stratakis, 2002). It might be possible that such a scenario was encountered in our intoxicated cells where forskolin was not sufficient, and Ca2+ was necessary, to promote cell death using Cry1Ca. It is worth noting that since the intracellular responses to Cry1 toxins studied using cell lines in culture remain controversial, this scenario will have to be validated in larvae of target insects. Finally, we observed in both selected resistant cell lines the upregulation of a gene described as the homologous of S. exigua activated C kinase receptor (RACK1). This PKC-activated binding protein (Mochly-Rosen, 1995) has been shown as a specific binding partner of a cAMP-phosphodiesterase (PDE) and was proposed as a potential point of cross-talk between the cAMP and PKC signalling pathways (Bird et al., 2010; Yarwood et al., 1999). RACK1 and its partners could, after intoxication by Cry1Ca, be the node of convergence of the two signalling pathways mentioned above.

In summary, Sf9-LD80 cells, which were the most resistant to Cry1Ca toxicity, were affected in cAMP production. Furthermore, our results have shown the importance of Ca2+ and Mg2+ ions in the toxicity of Cry1Ca. These ions may participate in different ways to the establishment of a cAMP-dependent signalling pathway following the binding of the toxin to the membrane of Sf9 cells. Both the signalling model and the pore forming model may contribute to Cry1Ca toxicity as was observed in the third model of toxin action proposed by Jurat-Fuentes and colleagues for Cry1A toxin, or more recently for the pore-forming α-Toxin from Clostridium septicum (Chakravorty et al., 2015; Jurat-Fuentes and Adang 2006). At low concentrations, Cry1Ca insertion in the plasma membrane could create pores, allowing Ca2+ entry in the cell and resulting in the activation of a cAMP-dependent cell death pathway. Future studies focusing on the identification and characterization of proteins involved in cAMP turnover, such as PDEs and their regulators, may identify new resistance factors, and contribute to a better understanding of the molecular mode of action of Cry1Ca and insect resistance mechanisms.

Insect cell lines culture and selection

Sf9 cell lines (ATCC® Number: CRL-1711™; Invitrogen) from ovarian tissue of pupal stage lepidopteran S.frugiperda were routinely maintained in complete medium (Grace's Insect Medium, Cambrex) containing 10% foetal bovine serum (FBS), penicillin (50 Ul/ml) and streptomycin (50 μg/ml), at 28°C, and spread when 80% confluence was reached. For resistant strain selection, cells were plated onto six-well plates until they reached 80% confluence and treated for 24 h at 28°C, either with Cry1Ca toxin at the lethal dose giving 50% dead cells (LD50), or at the lethal dose giving 80% dead cells (LD80). These lethal doses were calculated for Sf9wt when treated with various concentrations of Cry1Ca toxin (see below). Treated cells were further rinsed with Dulbecco's phosphate buffer saline (D-PBS) and allowed to recover for 2 days. After this recovery time, cells were harvested and spread in a new six-well plate for another round of selection. Forty rounds were done using the LD50 or LD80 toxin doses. Recovered resistant cells after the selection program were referred to as Sf9-LD50 and Sf9-LD80. Potential effects of aging were evaluated on cells grown for 40 rounds without toxin (Sf9p40) as compared to eight rounds (Sf9p8). S2 cells (R690-07, Invitrogen) derived from primary culture of late stage Drosophila melanogaster embryos were cultured at 24°C on Schneider medium (Invitrogen).

Cry1 toxin production and purification

Cry1Ab or Cry1Ca toxins were independently produced from Bt strain 407-(minus)pHT320 transformed with plasmid expressing the full sequence of the cry1Ca gene (Sanchis et al., 1999) or expressing the full sequence of the cry1Ab gene (gift from M. Gohar, Micalis Institute, INRA, France). Cells were grown for 65 h at 30°C in HCT media containing 0.3% glucose and 10 mg/l erythromycin, spun down by centrifugation at 6500 g for 15 min at 4°C and were passed through a Potter homogenizer. After five cycles of centrifugation/homogenisation, Cry inclusion bodies were loaded onto a discontinuous sucrose gradient (weight/volume: 79/72/67%) which was then centrifuged at 23,000 rpm for 19 h at 4°C in a SW27 rotor (Beckman). Cry inclusion bodies were collected at the interface, washed with distilled water and solubilized 1 h at 37°C in 50 mM Na2CO3/NaHCO3, 10 mM DTT pH 10.2. Cry concentration was determined using the Bradford method (Protein Assay Reagent, BioRad) and toxin activation was further performed via a 16 h incubation at 37°C under gentle agitation in 2% trypsin (weight/weight). Activated Cry were purified by size-exclusion chromatography on Äkta-FPLC (GE Healthcare).

Viability test

1.5×105 cells were plated per well onto a 96-well microassay plate (BD Falcon) and left for 24 h at 28°C. Then cells were incubated for 18 h in complete media containing increasing concentration of activated Cry1Ab or Cry1Ca toxins (from 10−1 to 10−5 g/l). Before addition of the toxins, cells were pre-incubated 30 min in fresh media or, when divalent ions effects were measured, in fresh media containing 5 mM divalent ion chelating agents; ethylene-diamine-tetra-acetic-acid (EDTA) or ethylene-glycol-bis(2-aminoethylether)-N,N,N,N′-tetra-acetic-acid (EGTA) (Sigma-Aldrich). After toxin incubation, cells were rapidly incubated with SYBR GoldTM dye (Molecular Probes) at a 1:10,000 dilution and observed using an inverted epifluorescence equipped microscope (Eclipse TE2000-U, Nikon) coupled to a CCD camera (Orca ER, Hamamatsu). Two pictures per toxin concentration were used to count both stained and unstained cells (300–400 total cells). Cell viability was calculated by the formula: (stained cells/total cells)×100. Data were then treated on the SigmaPlot 10.0 software (Systat Software) and typical dose-response sigmoid curves were obtained with non-linear estimation. Cry1Ca lethal doses giving 50% (LD50) or 80% (LD80) dead cells were directly deduced from the curves. Resistance ratio (RR) due to treatments were obtained by the formula: LD50 (Sf9+treatment)/LD50 (Sf9 in free media).

Microarray experimental design

Oligonucleotides were designed from 79148 ESTs sequences recovered from eight S. frugiperda tissues. The assembly analysis programme CAP3 identified 10092 contigs and singleton from these ESTs. The S. frugiperda microarray consisted of 9773 60-mer oligonucleotides (synthesized by Sigma-Aldrich) designed to match unique contigs or singletons and to suit hybridization conditions (average GC content: 46%; average Tm: 46°C). The experimental design of each comparison (Sf9-LD50 versus Sf9p40, Sf9-LD80 versus Sf9p40, Sf9p8 versus Sf9p40) was made using Cy3- and Cy5-labelled cDNA. Each comparison consisted of six microarrays, three biological replicates hybridized with dye swap (fully balanced dye swap design) and duplicate spots. The microarray were hybridized with cDNA prepared as described earlier (Le Goff et al., 2006) and scanned using GenePixPro scanner (Axon,version 3.01). Experimental data and associated microarray designs have been deposited in the NCBI Gene Expression Omnibus (accession number: GSE20507, http://www.ncbi.nlm.nih.gov/geo/) using Mediante database for data transfer (Le Brigand and Barbry, 2007).

Microarray data analysis

Data analysis was performed using the Bioconductor suite of statistical packages and limmaGUI (Gentleman et al., 2004; Wettenhall and Smyth, 2004). Specific expression intensity was obtained by subtracting the background intensity from the foreground intensity for each non-flag spot (all flagged spots were eliminated). The expression data were normalised by within-array normalization with the ‘loess method’ and between-array normalization using the ‘quantile method’ (Yang et al., 2003). The linear model for the series of arrays and empirical Bayes method was then applied for assessing differential expression (Smyth, 2004). The false discovery rate of the P-value for multiple tests was controlled using the Benjamini-Hochberg method. Differentially expressed genes were selected if the absolute value of log2-fold-change was greater than 1 and adjusted P-value below 0.01 and if the average intensity was greater than twice the average background. In order to provide an overall measure of evidence of differential expression, we used the Fisher's method for combining adjusted P-values from independent tests of significance of duplicate spots (Hess and Iyer, 2007). Lists of genes differentially up- or downregulated specifically in Sf9-LD50, specifically in Sf9-LD80 or common to both cell lines are available in the Supplementary information (Table S1).

Immunocytochemistry and live imaging

Cells were incubated for 40 min with complete media containing 3 mg/l of activated Cry1Ca toxin. Then cells were washed with D-PBS and fixed for 20 min in 4% paraformaldehyde/PBS solution at room temperature. After two rinses in D-PBS cells were quickly rinsed in 100% methanol and incubated in PBS containing 100 mM glycine. Cells were further incubated 2 h at room temperature with gentle agitation in PBS 10% normal goat serum. Immunostainings were performed with primary antibodies to Cry1Ca (rat serum obtained from an immunization program, Eurogentec, Belgium) at the 1:2000 dilution and secondary antibodies goat anti-rat IgG, Alexa Fluor 488 conjugated (Molecular Probes). Immunostaining without primary antibodies was performed as control. Cells were then rinsed with PBS and mounted under Mowiol® 4-88 based media (Polysciences). Images were obtained using an inverted epifluorescence equipped microscope (Eclipse TE2000-U, Nikon) coupled to a CCD camera (Orca ER, Hamamatsu). Four images per experiment were used to count both stained and unstained cells (150–200 total cells) and the percentage of stained cells was calculated by the formula: (stained cells/total cells)×100. Cell morphology of living Sf9 cells was examined directly in culture wells under light microscopy after 18 h incubation in complete media containing 5 mg/l of Cry1Ca with or without 5 mM of chelator agents.

cAMP measurement

Cells were plated onto a 24-well microassay plate (BD Falcon) until reaching 80% confluency, washed with D-PBS and incubated for 20 min at 28°C with complete media only (control, CT) or containing Cry1Ca toxin (1 mg/l) or forskolin (100 µM). Cells were further incubated 10 min in HCl 0.1 M at 4°C, scrapped, vigorously shaken and centrifuged at 600×g for 10 min at 4°C. Supernatants were collected and proteins were quantified by the Bradford method assay. Amounts of cAMP/mg of protein were measured in HCl extracts using the Direct cAMP enzyme immunoassay kit according to the manufacturer instruction (Sigma-Aldrich). The cAMP fold increases were calculated by the ratio of cAMP amount in stimulated conditions divided by the cAMP amount in non stimulated cells (−).

Author contributions

Conceptualization: C.C., D.P., F.H., P.B.; Methodology: C.C., D.P., F.H., V.T.T., N.Z.-P., P.B.; Software: C.C., F.H., V.T.T.; Validation: C.C., V.T.T., P.B.; Formal analysis: C.C., V.T.T., P.B.; Investigation: C.C., P.B.; Resources: D.P., A.G.; Data curation: P.B.; Writing - original draft: D.P., A.G., P.B.; Supervision: P.B.

Funding

This research received no specific grant from any funding agency in the public, commercial, or not-for-profit sectors.

Abdullah
,
M. A.
,
Moussa
,
S.
,
Taylor
,
M. D.
and
Adang
,
M. J.
(
2009
).
Manduca sexta (Lepidoptera: Sphingidae) cadherin fragments function as synergists for Cry1A and Cry1C Bacillus thuringiensis toxins against noctuid moths Helicoverpa zea, Agrotis ipsilon and Spodoptera exigua
.
Pest Manag. Sci.
65
,
1097
-
1103
.
Agrawal
,
N.
,
Malhotra
,
P.
and
Bhatnagar
,
R. K.
(
2002
).
Interaction of gene-cloned and insect cell-expressed aminopeptidase N of spodoptera litura with insecticidal crystal protein Cry1C
.
Appl. Environ. Microbiol.
68
,
4583
-
4592
.
Alcantara
,
E. P.
,
Aguda
,
R. M.
,
Curtiss
,
A.
,
Dean
,
D. H.
and
Cohen
,
M. B.
(
2004
).
Bacillus thuringiensis?-endotoxin binding to brush border membrane vesicles of rice stem borers
.
Arch. Insect Biochem. Physiol.
55
,
169
-
177
.
Avisar
,
D.
,
Segal
,
M.
,
Sneh
,
B.
and
Zilberstein
,
A.
(
2005
).
Cell-cycle-dependent resistance to Bacillus thuringiensis Cry1C toxin in Sf9 cells
.
J. Cell Sci.
118
,
3163
-
3171
.
Bird
,
R. J.
,
Baillie
,
G. S.
and
Yarwood
,
S. J.
(
2010
).
Interaction with receptor for activated C-kinase 1 (RACK1) sensitizes the phosphodiesterase PDE4D5 towards hydrolysis of cAMP and activation by protein kinase C
.
Biochem. J.
432
,
207
-
216
.
Braun
,
T.
,
Frank
,
H.
,
Dods
,
R.
and
Sepsenwol
,
S.
(
1977
).
Mn2+-sensitive, soluble adenylate cyclase in rat testis Differentiation from other testicular nucleotide cyclases
.
Biochim. Biophys. Acta
481
,
227
-
235
.
Bretschneider
,
A.
,
Heckel
,
D. G.
,
Pauchet
,
Y.
(
2016
).
Three toxins, two receptors, one mechanism: Mode of action of Cry1A toxins from Bacillus thuringiensis in Heliothis virescens
.
Insect Biochem Mol Biol.
76
,
109
-
117
.
Candas
,
M.
,
Loseva
,
O.
,
Oppert
,
B.
,
Kosaraju
,
P.
and
Bulla
,
L. A.
 Jr.
(
2003
).
Insect resistance to bacillus thuringiensis
.
Mol. Cell. Proteomics
2
,
19
-
28
.
Canton
,
P. E.
,
Cancino-Rodezno
,
A.
,
Gill
,
S. S.
,
Soberón
,
M.
and
Bravo
,
A.
(
2015
).
Transcriptional cellular responses in midgut tissue of Aedes aegypti larvae following intoxication with Cry11Aa toxin from Bacillus thuringiensis
.
BMC Genomics
16
,
1042
.
Cao
,
J.
,
Zhao
,
J.-Z.
,
Tang
,
J. D.
,
Shelton
,
A. M.
and
Earle
,
E. D.
(
2002
).
Broccoli plants with pyramided Cry1Ac and Cry1C Bt genes control diamondback moths resistant to Cry1A and Cry1C proteins
.
Theor. Appl. Genet.
105
,
258
-
264
.
Chakravorty
,
A.
,
Awad
,
M. M.
,
Cheung
,
J. K.
,
Hiscox
,
T. J.
,
Lyras
,
D.
and
Rood
,
J. I.
(
2015
).
The pore-forming α-toxin from clostridium septicum activates the MAPK pathway in a Ras-c-Raf-dependent and independent manner
.
Toxins
7
,
516
-
534
.
Chattopadhyay
,
P.
and
Banerjee
,
G.
(
2018
).
Recent advancement on chemical arsenal of Bt toxin and its application in pest management system in agricultural field
.
3 Biotech
,
8
,
201
.
Coleman
,
M. L.
,
Marshall
,
C. J.
and
Olson
,
M. F.
(
2004
).
RAS and RHO GTPases in G1-phase cell-cycle regulation
.
Nat. Rev. Mol. Cell Biol.
5
,
355
-
366
.
Crickmore
,
N.
,
Zeigler
,
D. R.
,
Feitelson
,
J.
,
Schnepf
,
E.
,
Van Rie
,
J.
,
Lereclus
,
D.
,
Baum
,
J.
,
Dean
,
D. H.
(
1998
).
Revision of the nomenclature for the Bacillus thuringiensis pesticidal crystal proteins
.
Microbiol. Mol. Biol. Rev.
62
,
807
-
813
.
Farias
,
J. R.
,
Horikoshi
,
R. J.
,
Santos
,
A. C.
and
Omoto
,
C.
(
2014
).
Geographical and temporal variability in susceptibility to Cry1F toxin from bacillus thuringiensis in Spodoptera frugiperda (Lepidoptera: Noctuidae) populations in Brazil
.
J. Econ. Entomol.
107
,
2182
.
Finkler
,
A.
,
Ashery-Padan
,
R.
and
Fromm
,
H.
(
2007
).
CAMTAs: Calmodulin-binding transcription activators from plants to human
.
FEBS Lett.
581
,
3893
-
3898
.
Gahan
,
L. J.
,
Gould
,
F.
and
Heckel
,
D. G.
(
2001
).
Identification of a gene associated with Bt Resistance in Heliothis virescens
.
Science
293
,
857
-
860
.
Gentleman
,
R. C.
,
Carey
,
V. J.
,
Bates
,
D. M.
,
Bolstad
,
B.
,
Dettling
,
M.
,
Dudoit
,
S.
,
Ellis
,
B.
,
Gautier
,
L.
,
Ge
,
Y.
,
Gentry
,
J.
, et al. 
(
2004
).
Bioconductor: open software development for computational biology and bioinformatics
.
Genome Biol.
5
,
R80
.
Goergen
,
G.
,
Kumar
,
P. L.
,
Sankung
,
S. B.
,
Togola
,
A.
and
Tamo
,
M.
(
2016
).
First Report of Outbreaks of the Fall Armyworm Spodoptera frugiperda (J E Smith) (Lepidoptera, Noctuidae), a New Alien Invasive Pest in West and Central Africa
.
PLoS ONE
11
,
e0165632
.
Guihard
,
G.
,
Vachon
,
V.
,
Laprade
,
R.
and
Schwartz
,
J. L.
(
2000
).
Kinetic properties of the channels formed by the bacillus thuringiensis insecticidal crystal protein Cry1C in the plasma membrane of Sf9 cells
.
J. Membr. Biol.
175
,
115
-
122
.
Heckel
,
D. G.
(
2012
).
Learning the ABCs of Bt: ABC transporters and insect resistance to Bacillus thuringiensis provide clues to a crucial step in toxin mode of action
.
Pesticide Biochemistry and Physiology Special Issue: Molecular Approaches to Pest Control, Toxicology and Resistance. Pestic Biochem Phys.
104
,
103
-
110
.
Heckel
,
D. G.
,
Gahan
,
L. J.
,
Baxter
,
S. W.
,
Zhao
,
J. Z.
,
Shelton
,
A. M.
,
Gould
,
F.
and
Tabashnik
,
B. E.
(
2007
).
The diversity of Bt resistance genes in species of Lepidoptera
.
J. Invertebr. Pathol.
95
,
192
-
197
.
Hess
,
A.
and
Iyer
,
H.
(
2007
).
Fisher's combined p-value for detecting differentially expressed genes using Affymetrix expression arrays
.
BMC Genomics
8
,
96
.
Hofer
,
A. M.
(
2012
).
Interactions between calcium and cAMP signaling
.
Curr
19
,
5768
-
5773
.
Jiménez-Juárez
,
N.
,
Muñoz-Garay
,
C.
,
Gómez
,
I.
,
Saab-Rincon
,
G.
,
Damian-Almazo
,
J. Y.
,
Gill
,
S. S.
,
Soberón
,
M.
and
Bravo
,
A.
(
2007
).
Bacillus thuringiensis Cry1Ab mutants affecting oligomer formation are non-toxic to manduca sexta larvae
.
J. Biol. Chem.
282
,
21222
-
21229
.
Jurat-Fuentes
,
J. L.
and
Adang
,
M. J.
(
2006
).
Cry toxin mode of action in susceptible and resistant Heliothis virescens larvae
.
J. Invertebr. Pathol.
92
,
166
-
171
.
Jurat-Fuentes
,
J. L.
,
Karumbaiah
,
L.
,
Jakka
,
S. R.
,
Ning
,
C.
,
Liu
,
C.
,
Wu
,
K.
,
Jackson
,
J.
,
Gould
,
F.
,
Blanco
,
C.
,
Portilla
,
M.
, et al. 
(
2011
).
Reduced levels of membrane-bound alkaline phosphatase are common to lepidopteran strains resistant to Cry toxins from bacillus thuringiensis
.
PLoS ONE
6
,
e17606
.
Knowles
,
B. H.
and
Farndale
,
R. W.
(
1988
).
Activation of insect cell adenylate cyclase by Bacillus thuringiensis δ-endotoxins and melittin. Toxicity is independent of cyclic AMP
.
Biochem. J.
253
,
235
-
241
.
Kwa
,
M. S.
,
de Maagd
,
R. A.
,
Stiekema
,
W. J.
,
Vlak
,
J. M.
and
Bosch
,
D.
(
1998
).
Toxicity and binding properties of the bacillus thuringiensis delta-endotoxin Cry1C to cultured insect cells
.
J. Invertebr. Pathol.
71
,
121
-
127
.
Laflamme
,
E.
,
Badia
,
A.
,
Lafleur
,
M.
,
Schwartz
,
J. L.
and
Laprade
,
R.
(
2008
).
Atomic force microscopy imaging of bacillus thuringiensis Cry1 toxins interacting with insect midgut apical membranes
.
J. Membr. Biol.
222
,
127
-
139
.
Le Brigand
,
K.
and
Barbry
,
P.
(
2007
).
Mediante: a web-based microarray data manager
.
Bioinformatics
23
,
1304
-
1306
.
Le Goff
,
G.
,
Hilliou
,
F.
,
Siegfried
,
B. D.
,
Boundy
,
S.
,
Wajnberg
,
E.
,
Sofer
,
L.
,
Audant
,
P.
,
ffrench-Constant
,
R. H.
and
Feyereisen
,
R.
(
2006
).
Xenobiotic response in Drosophila melanogaster: sex dependence of P450 and GST gene induction
.
Insect Biochem. Mol. Biol.
36
,
674
-
682
.
Li
,
H.
,
Rao
,
A.
and
Hogan
,
P. G.
(
2011
).
Interaction of calcineurin with substrates and targeting proteins
.
Trends Cell Biol.
21
,
91
-
103
.
Litvin
,
T. N.
,
Kamenetsky
,
M.
,
Zarifyan
,
A.
,
Buck
,
J.
and
Levin
,
L. R.
(
2003
).
Kinetic properties of “soluble” adenylyl cyclase
.
J. Biol. Chem.
278
,
15922
-
15926
.
Liu
,
Y.-B.
,
Tabashnik
,
B. E.
,
Masson
,
L.
,
Escriche
,
B.
and
Ferré
,
J.
(
2000
).
Binding and toxicity of bacillus thuringiensis protein Cry1C to susceptible and resistant diamondback Moth (Lepidoptera: Plutellidae)
.
J. Econ. Entomol.
93
,
1
-
6
.
Liu
,
Y.-B.
,
Tabashnik
,
B. E.
,
Meyer
,
S. K.
and
Crickmore
,
N.
(
2001
).
Cross-resistance and stability of resistance to bacillus thuringiensis Toxin Cry1C in diamondback moth
.
Appl. Environ. Microbiol.
67
,
3216
-
3219
.
Luo
,
K.
,
Lu
,
Y.-J.
and
Adang
,
M. J.
(
1996
).
A 106 kDa form of aminopeptidase is a receptor for Bacillus thuringiensis CryIC δ-endotoxin in the brush border membrane of Manduca sexta
.
Insect Biochem. Mol. Biol.
26
,
783
-
791
.
Masson
,
L.
,
Lu
,
Y.-J.
,
Mazza
,
A.
,
Brousseau
,
R.
and
Adang
,
M. J.
(
1995
).
The CryIA(c) receptor purified from manduca sexta displays multiple specificities
.
J. Biol. Chem.
270
,
20309
-
20315
.
Mochly-Rosen
,
D.
(
1995
).
Localization of protein kinases by anchoring proteins: a theme in signal transduction
.
Science
268
,
247
-
251
.
Monette
,
R.
,
Potvin
,
L.
,
Baines
,
D.
,
Laprade
,
R.
and
Schwartz
,
J. L.
(
1997
).
Interaction between Calcium Ions and Bacillus thuringiensis Toxin Activity against Sf9 Cells (Spodoptera frugiperda, Lepidoptera)
.
Appl. Environ. Microbiol.
63
,
440
-
447
.
Morin
,
S.
,
Biggs
,
R. W.
,
Sisterson
,
M. S.
,
Shriver
,
L.
,
Ellers-Kirk
,
C.
,
Higginson
,
D.
,
Holley
,
D.
,
Gahan
,
L. J.
,
Heckel
,
D. G.
,
Carriere
,
Y.
, et al. 
(
2003
).
Three cadherin alleles associated with resistance to Bacillus thuringiensis in pink bollworm
.
Proc. Natl. Acad. Sci. USA
100
,
5004
-
5009
.
Nagoshi
,
R. N.
and
Meagher
,
R. L.
(
2008
).
Review of Fall Armyworm (Lepidoptera: Noctuidae) Genetic Complexity and Migration
.
Fla. Entomol.
91
,
546
-
554
.
Nanoth Vellichirammal
,
N.
,
Wang
,
H.
,
Eyun
,
S. I.
,
Moriyama
,
E. N.
,
Coates
,
B. S.
,
Miller
,
N. J.
and
Siegfried
,
B. D.
(
2015
).
Transcriptional analysis of susceptible and resistant European corn borer strains and their response to Cry1F protoxin
.
BMC Genomics
16
,
558
.
Näsman
,
J.
,
Kukkonen
,
J. P.
and
Åkerman
,
K. E.
(
2002
).
Dual signalling by different octopamine receptors converges on adenylate cyclase in Sf9 cells
.
Insect Biochem. Mol. Biol.
32
,
28
-
93
.
Palma
,
L.
,
Muñoz
,
D.
,
Berry
,
C.
,
Murillo
,
J.
and
Caballero
,
P.
(
2014
).
Bacillus thuringiensis toxins: an overview of their biocidal activity
.
Toxins
6
,
3296
-
3325
.
Peyronnet
,
O.
,
Vachon
,
V.
,
Schwartz
,
J.-L.
and
Laprade
,
R.
(
2001
).
Ion channels induced in planar lipid bilayers by the bacillus thuringiensis toxin Cry1Aa in the presence of gypsy moth (Lymantria dispar) brush border membrane
.
J. Membr. Biol.
184
,
45
-
54
.
Pieroni
,
J. P.
,
Harry
,
A.
,
Chen
,
J.
,
Jacobowitz
,
O.
,
Magnusson
,
R. P.
and
Iyengar
,
R.
(
1995
).
Distinct characteristics of the basal activities of adenylyl Cyclases 2 and z6
.
J. Biol. Chem.
270
,
21368
-
21373
.
Pigott
,
C. R.
and
Ellar
,
D. J.
(
2007
).
Role of receptors in bacillus thuringiensis crystal toxin activity
.
Microbiol. Mol. Biol. Rev.
71
,
255
-
281
.
Portugal
,
L.
,
Munoz-Garay
,
C.
,
Martínez de Castro
,
D. L.
,
Soberón
,
M.
and
Bravo
,
A.
(
2017
).
Toxicity of Cry1A toxins from Bacillus thuringiensis to CF1 cells does not involve activation of adenylate cyclase/PKA signaling pathway
.
Insect Biochem. Mol. Biol.
80
,
21
-
31
.
Pu
,
L.
,
Igbavboa
,
U.
,
Wood
,
W. G.
,
Roths
,
J. B.
,
Kier
,
A. B.
,
Spener
,
F.
and
Schroeder
,
F.
(
1999
).
Expression of fatty acid binding proteins is altered in aged mouse brain
.
Mol. Cell. Biochem.
198
,
69
-
78
.
Raymond
,
B.
,
Johnston
,
P. R.
,
Nielsen-LeRoux
,
C.
,
Lereclus
,
D.
and
Crickmore
,
N.
(
2010
).
Bacillus thuringiensis: an impotent pathogen?
Trends Microbiol.
18
,
189
-
194
.
Robinson-White
,
A.
and
Stratakis
,
C. A.
(
2002
).
Protein kinase a signaling
.
Ann. N. Y. Acad. Sci.
968
,
256
-
270
.
Sanchis
,
V.
,
Gohar
,
M.
,
Chaufaux
,
J.
,
Arantes
,
O.
,
Meier
,
A.
,
Agaisse
,
H.
,
Cayley
,
J.
and
Lereclus
,
D.
(
1999
).
Development and field performance of a broad-spectrum nonviable asporogenic recombinant strain of Bacillus thuringiensis with greater potency and UV resistance
.
Appl. Environ. Microbiol.
65
,
4032
-
4039
.
Schwartz
,
J. L.
,
Garneau
,
L.
,
Masson
,
L.
and
Brousseau
,
R.
(
1991
).
Early response of cultured lepidopteran cells to exposure to δ-endotoxin from Bacillus thuringiensis: Involvement of calcium and anionic channels
.
Biochim. Biophys. Acta
1065
,
250
-
260
.
Smyth
,
G. K.
(
2004
).
Linear models and empirical bayes methods for assessing differential expression in microarray experiments
.
Stat. Appl. Genet. Mol. Biol.
3
,
Article3
.
Soberón
,
M.
,
Perez
,
R. V.
,
Nunez-Valdez
,
M. E.
,
Lorence
,
A.
,
Gomez
,
I.
,
Sanchez
,
J.
and
Bravo
,
A.
(
2000
).
Evidence for intermolecular interaction as a necessary step for pore-formation activity and toxicity of Bacillus thuringiensis Cry1Ab toxin
.
FEMS Microbiol. Lett.
191
,
221
-
225
.
Sprang
,
S. R.
(
1997
).
G proteins, effectors and GAPs: structure and mechanism
.
Curr. Opin. Struct. Biol.
7
,
849
-
856
.
Storer
,
N. P.
,
Babcock
,
J. M.
,
Schlenz
,
M.
,
Meade
,
T.
,
Thompson
,
G. D.
,
Bing
,
J. W.
and
Huckaba
,
R. M.
(
2010
).
Discovery and characterization of field resistance to Bt Maize: Spodoptera frugiperda (Lepidoptera: Noctuidae) in puerto rico
.
J. Econ. Entomol.
103
,
1031
-
1038
.
Tang
,
H.
,
Chen
,
G.
,
Chen
,
F.
,
Han
,
L.
and
Peng
,
Y.
(
2017
).
Development and relative fitness of Cry1C resistance in Chilo suppressalis
.
Pest Manag. Sci.
74
,
590
-
597
.
Tesmer
,
J. J.
,
Sunahara
,
R. K.
,
Johnson
,
R. A.
,
Gosselin
,
G.
,
Gilman
,
A. G.
and
Sprang
,
S. R.
(
1999
).
Two-metal-ion catalysis in adenylyl cyclase
.
Science
285
,
756
-
760
.
Vachon
,
V.
,
Laprade
,
R.
,
Schwartz
,
J. L.
,
Fortea
,
E.
,
Lemieux
,
V.
,
Potvin
,
L.
,
Chikwana
,
V.
,
Griffin
,
S.
,
Hey
,
T.
,
McCaskill
,
D.
, et al. 
(
2012
).
Current models of the mode of action of Bacillus thuringiensis insecticidal crystal proteins: a critical review
.
J. Invertebr. Pathol.
111
,
1
-
12
.
Villalon
,
M.
,
Vachon
,
V.
,
Brousseau
,
R.
,
Schwartz
,
J. L.
and
Laprade
,
R.
(
1998
).
Video imaging analysis of the plasma membrane permeabilizing effects of Bacillus thuringiensis insecticidal toxins in Sf9 cells
.
Biochim. Biophys. Acta
1368
,
27
-
34
.
Webb
,
B. L.
,
Hirst
,
S. J.
and
Giembycz
,
M. A.
(
2000
).
Protein kinase C isoenzymes: a review of their structure, regulation and role in regulating airways smooth muscle tone and mitogenesis
.
Br. J. Pharmacol.
130
,
1433
-
1452
.
Wettenhall
,
J. M.
and
Smyth
,
G. K.
(
2004
).
limmaGUI: a graphical user interface for linear modeling of microarray data
.
Bioinformatics
20
,
3705
-
3706
.
Willoughby
,
D.
and
Cooper
,
D. M.
(
2007
).
Organization and Ca 2+ regulation of adenylyl cyclases in cAMP microdomains
.
Physiol. Rev.
87
,
965
-
1010
.
Xue
,
J.-L.
,
Cai
,
Q.-X.
,
Zheng
,
D.-S.
and
Yuan
,
Z.-M.
(
2005
).
The synergistic activity between Cry1Aa and Cry1c from Bacillus thuringiensis against Spodoptera exigua and Helicoverpa armigera
.
Lett. Appl. Microbiol.
40
,
460
-
465
.
Yang
,
S.-A.
and
Klee
,
C. B.
(
2000
).
Low affinity Ca 2+ -binding sites of calcineurin B mediate conformational changes in calcineurin A
.
Biochemistry
39
,
16147
-
16154
.
Yang
,
C.-Y.
,
Pang
,
J.-C.
,
Kao
,
S.-S.
and
Tsen
,
H.-Y.
(
2003
).
Enterotoxigenicity and cytotoxicity of bacillus thuringiensis strains and development of a process for Cry1Ac production
.
J. Agric. Food Chem.
51
,
100
-
105
.
Yarwood
,
S. J.
,
Steele
,
M. R.
,
Scotland
,
G.
,
Houslay
,
M. D.
and
Bolger
,
G. B.
(
1999
).
The RACK1 signaling scaffold protein selectively interacts with the cAMP-specific phosphodiesterase PDE4D5 isoform
.
J. Biol. Chem.
274
,
14909
-
14917
.
Yoon
,
Y. S.
,
Choo
,
J. H.
,
Yoo
,
T.
,
Kang
,
K.
and
Chung
,
J. H.
(
2007
).
RhoB is epigenetically regulated in an age- and tissue-specific manner
.
Biochem. Biophys. Res. Commun.
362
,
164
-
169
.
Zhang
,
X.
,
Candas
,
M.
,
Griko
,
N. B.
,
Rose-Young
,
L.
and
Bulla
,
L. A.
(
2005
).
Cytotoxicity of Bacillus thuringiensis Cry1Ab toxin depends on specific binding of the toxin to the cadherin receptor BT-R1 expressed in insect cells
.
Cell Death Differ.
12
,
1407
-
1416
.
Zhang
,
X.
,
Candas
,
M.
,
Griko
,
N. B.
,
Taussig
,
R.
and
Bulla
,
L. A.
 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
.
Zhuang
,
M.
,
Oltean
,
D. I.
,
Gómez
,
I.
,
Pullikuth
,
A. K.
,
Soberón
,
M.
,
Bravo
,
A.
and
Gill
,
S. S.
(
2002
).
Heliothis virescens and manduca sexta lipid rafts are involved in Cry1A toxin binding to the midgut epithelium and subsequent pore formation
.
J. Biol. Chem.
277
,
13863
-
13872
.

Competing interests

The authors declare no competing or financial interests.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution and reproduction in any medium provided that the original work is properly attributed.

Supplementary information