# Transcriptomic analysis of Spodoptera frugiperda Sf9 cells resistant to Bacillus thuringiensis Cry1Ca toxin reveals that extracellular Ca²⁺, Mg²⁺ and production of cAMP are involved in toxicity

Claude Castella, David Pauron, Frederique Hilliou, van Tran Trang, Nathalie Zucchini-Pascal, Armel Gallet, Pierre Barbero

<figure><img src="image_1.png" /></figure>

To cite this version:

Claude Castella, David Pauron, Frederique Hilliou, van Tran Trang, Nathalie Zucchini-Pascal, et al. Transcriptomic analysis of Spodoptera frugiperda Sf9 cells resistant to Bacillus thuringiensis Cry1Ca toxin reveals that extracellular Ca²⁺, Mg²⁺ and production of cAMP are involved in toxicity. Biology Open, 2019, 8 (4), pp.bio037085. 10.1242/bio.037085. hal-02322389

HAL Id: hal-02322389

https://hal.science/hal-02322389

Submitted on 26 May 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
RESEARCH ARTICLE

# Transcriptomic analysis of *Spodoptera frugiperda* Sf9 cells resistant to *Bacillus thuringiensis* Cry1Ca toxin reveals that extracellular Ca²⁺, Mg²⁺ and production of cAMP are involved in toxicity

Claude Castella, David Pauron, Frédérique Hilliou, Van Tran Trang, Nathalie Zucchini-Pascal, Armel Gallet and Pierre Barbero*

## ABSTRACT

*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 Mg²⁺-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 Ca²⁺ and/or Mg²⁺ 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.

**KEY WORDS:** *Bacillus thuringiensis*, Cry1Ca, Sf9 cells, Resistance, Toxins, Transcriptomic

## INTRODUCTION

*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 Mg²⁺-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 (cAMPᵢ). 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

Université Côte d'Azur, INRA, CNRS, ISA, France.

*Author for correspondence (barbero@unice.fr)

<figure><img src="image_3.png" /></figure>

P.B., 0000-0002-7362-3968

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

Received 10 July 2018; Accepted 22 March 2019
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.

## RESULTS

### 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

$LD_{50}$ 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 ($LD_{50}$) 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

<figure><img src="image_4.png" /><figcaption>**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 (V). Main graph shows cell viability measurements of selected resistant Sf9 cell lines $Sf9-LD_{50}$ (O) and $Sf9-LD_{80}$ (▼) 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 ($LD_{50}$) 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.</figcaption></figure>
program to develop cell strains resistant to this toxin. Doses calculated for Sf9wt cells and corresponding to the LD₅₀ and LD₈₀ were used to select cells over 40 generations. At the end of the selection program, resistant cell lines obtained using LD₅₀ or LD₈₀ doses were named Sf9-LD₅₀ and Sf9-LD₈₀ respectively. Cell viability was tested at the 40th generation and LD₅₀ measured for Sf9-LD₅₀ (white circles in Fig. 1A) and Sf9-LD₈₀ 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 LD₅₀Sf9-resistant/LD₅₀Sf9wt ratio: RF=4.9±0.8 for Sf9-LD₅₀ cells and RF=58.9±2.8 for Sf9-LD₈₀ 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-LD₅₀ cells were affected and no change in Sf9-LD₈₀ cell density and morphology was observed.

### 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 2 few 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-LD₅₀ and Sf9-LD₈₀ cells is probably not caused by loss of a plasma membrane receptor.

### 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-LD₅₀ to 8% in Sf9-LD₈₀) and to transcription processes (from 6% in Sf9-LD₅₀ to 14% in Sf9-LD₈₀) 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-LD₅₀ and in Sf9-LD₈₀ 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

<figure><img src="image_5.png" /><figcaption>**Fig. 2.** Cry1Ca binding to the plasma membrane of Sf9 and S2 cell lines, Sf9wt cells or resistant cell lines Sf9-LD₅₀ or Sf9-LD₈₀ were incubated 40 min in 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.</figcaption></figure>
<figure><img src="image_6.png" /><figcaption>A</figcaption></figure>

**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-LD₅₀ and Sf9-LD₈₀ 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-LD₅₀ only (82), in Sf9-LD₈₀ only (125) and common in both resistant cell lines (49).

<figure><img src="image_7.png" /><figcaption>B</figcaption></figure>

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-LD₈₀ only; the selected cell lines with the highest calculated RF.

### 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 Mg²⁺ and Ca²⁺, whereas EGTA preferentially chelates Ca²⁺. No detectable mortality was observed when Sf9wt or Sf9-LD₈₀ 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 LD₅₀s 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 LD₅₀s 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-LD₈₀ 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-LD₈₀ 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 Mg²⁺ and Ca²⁺ in Cry1Ca action. This protective effect was not likely due to loss of interaction of Cry1Ca with plasma membrane binding regions.

### Cry1Ca toxin induces Cyclic AMP production

Previous studies proposed that Cry protein toxicity involves activation of an Mg²⁺-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
<table><caption>Table 1. List of genes with the same expression profile in both resistant cell lines</caption><thead><tr><th>Accession No.</th><th>Description</th><th>Expression</th><th>Oligo ID</th></tr></thead><tbody><tr><td>KJ622099.1</td><td>Agrobiotis segetum enoyl-CoA hydratase</td><td>+</td><td>adhoc-25862</td></tr><tr><td>NP_001124379</td><td>Bombix mori diuretic hormone 31</td><td>+</td><td>adhoc-40459</td></tr><tr><td>NP_001091832</td><td>Bombyx mori eukaryotic initiation factor 4E-1</td><td>+</td><td>adhoc-27125</td></tr><tr><td><strong>XM_021336943</strong></td><td><strong>Helicoverpa armigera organic cation transporteur protein-like</strong></td><td><strong>+</strong></td><td><strong>adhoc-27137</strong></td></tr><tr><td><strong>NG_053148.1</strong></td><td><strong>Homo sapiens calmodulin binding transcription activator 1</strong></td><td><strong>+</strong></td><td><strong>adhoc-38874</strong></td></tr><tr><td>XM_004922449.3</td><td>Predicted: Bombyx mori zinc finger protein 608</td><td>+</td><td>adhoc-34656</td></tr><tr><td>XM_002763546.4</td><td>Predicted: Callithrix jacchus UDP glucuronosyltransferase family 2</td><td>+</td><td>adhoc-38447</td></tr><tr><td>XM_021328952.1</td><td>Predicted: Helicoverpa armigera eukaryotic translation initiation factor 3</td><td>+</td><td>adhoc-40626</td></tr><tr><td>XM_021340684.1</td><td>Predicted: Helicoverpa armigera receptor expression-enhancing protein 5</td><td>+</td><td>adhoc-38756</td></tr><tr><td>XM_021342409.1</td><td>Predicted: Helicoverpa armigera RNA-binding protein pno1</td><td>+</td><td>adhoc-35838</td></tr><tr><td>XM_021329746.1</td><td>Predicted: Helicoverpa armigera silencing protein 10</td><td>+</td><td>adhoc-40628</td></tr><tr><td>XM_021325195.1</td><td>Predicted: Helicoverpa armigera uncharacterized LOC110369688</td><td>+</td><td>adhoc-33957</td></tr><tr><td>XM_021335084.1</td><td>Predicted: Helicoverpa armigera uncharacterized LOC110376574</td><td>+</td><td>adhoc-44429</td></tr><tr><td>XM_014576615.1</td><td>Predicted: Pelodiscus sinensis cartilage intermediate layer protein</td><td>+</td><td>adhoc-34293</td></tr><tr><td>XM_011570118.1</td><td>Predicted: Plutella xylostella coatomer subunit delta</td><td>+</td><td>adhoc-34845</td></tr><tr><td>XM_011560908.1</td><td>Predicted: Plutella xylostella protein BCCIP homolog</td><td>+</td><td>adhoc-40531</td></tr><tr><td>XM_011159698.1</td><td>Predicted: Solonopsis invicta venom serine carboxypeptidase</td><td>+</td><td>adhoc-40864</td></tr><tr><td>XP_967064</td><td>Predicted: Tribolium castaneum ran-binding protein homolog</td><td>+</td><td>adhoc-36521</td></tr><tr><td><strong>KC188795.1</strong></td><td><strong>Spodoptera exigua 14-3-3 protein epsilon</strong></td><td><strong>+</strong></td><td><strong>adhoc-36480</strong></td></tr><tr><td><strong>EU259815.1</strong></td><td><strong>Spodoptera exigua activated C kinase 1 receptor (RACK1)</strong></td><td><strong>+</strong></td><td><strong>adhoc-36385</strong></td></tr><tr><td>JN164650.1</td><td>Spodoptera frugiperda cap'n'collar isoform C</td><td>+</td><td>adhoc-39822</td></tr><tr><td><strong>DQ295917.1</strong></td><td><strong>Spodoptera frugiperda G protein-coupled receptor</strong></td><td><strong>+</strong></td><td><strong>adhoc-40926</strong></td></tr><tr><td>FO681374.1</td><td>Spodoptera frugiperda sequence from BAC clone 96D18</td><td>+</td><td>adhoc-41072</td></tr><tr><td>JQ979047.1</td><td>Spodoptora littoralis cytochrome P450 CYP305B1</td><td>+</td><td>adhoc-39913</td></tr><tr><td>XM_013328524.1</td><td>Predicted Amyelois transitella epidermal retinol dehydrogenase 2- like</td><td>—</td><td>adhoc-38974</td></tr><tr><td>XM_013330116.1</td><td>Predicted Amyelois transitella translation initiation factor 3 subunit M</td><td>—</td><td>adhoc-34731</td></tr><tr><td><strong>XM_021334894.1</strong></td><td><strong>Predicted: Helicoverpa armigera calcineurin B homologous protein 1</strong></td><td><strong>—</strong></td><td><strong>adhoc-34328</strong></td></tr><tr><td>XM_021342567.1</td><td>Predicted: Helicoverpa armigera eukaryotic translation initiation factor 5A</td><td>—</td><td>adhoc-39550</td></tr><tr><td>XM_021326538.1</td><td>Predicted: Helicoverpa armigera retinoblastoma-like protein 1</td><td>—</td><td>adhoc-27605</td></tr><tr><td>XM_021325970.1</td><td>Predicted: Helicoverpa armigera transcription factor E2F7-like</td><td>—</td><td>adhoc-38939</td></tr><tr><td>XM_021338618.1</td><td>Predicted: Helicoverpa armigera ubiquinone monooxygenase COQ6</td><td>—</td><td>adhoc-38013</td></tr><tr><td>XM_021337093.1</td><td>Predicted: Helicoverpa armigera uncharacterized protein</td><td>—</td><td>adhoc-44377</td></tr><tr><td>XM_021333277.1</td><td>Predicted: Helicoverpa armigera importin subunit alpha-7</td><td>—</td><td>adhoc-27546</td></tr><tr><td>KT218668.1</td><td>Spodoptera frugiperda ribosomal protein L23A</td><td>—</td><td>adhoc-34322</td></tr><tr><td>AY072287.1</td><td>Spodoptera frugiperda ribosomal protein L3</td><td>—</td><td>adhoc-41009</td></tr><tr><td>AY646107.1</td><td>Trichoplusia ni putative cytoplasmic actin variant 2</td><td>—</td><td>adhoc-44378</td></tr></tbody></table>

36 identified genes out of 49 genes common to Sf9-LD₅₀ and Sf9-LD₈₀ cells and regulated the same way in both resistant selected cell lines are listed. Bold lines highlight genes potentially involved in cellular pathways requiring or modifying cation intracellular concentrations. Upregulated, +; downregulated, -. Oligo IDs correspond to the identity of the oligonucleotide sequence spotted onto the S. frugiperda microarray.

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 Ca²⁺.

### Sf9-LD₈₀ cells are impaired in cAMP production

Since Cry1Ca action seemed to require increase of intracellular concentration of cAMP ([cAMP]ᵢ) we wondered whether the resistant cell line Sf9-LD₈₀ was able to produce this second messenger. Sf9wt sensitive cells or Sf9-LD₈₀ 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-LD₈₀ 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-LD₈₀ cells could

overcome Cry1Ca toxicity by blocking cAMP-dependent signalling normally triggered by the toxin in sensitive cells.

### DISCUSSION

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 LD₅₀ 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 LD₅₀ or LD₈₀ Cry1Ca doses and we obtained cell lines with an increasing resistance factor, named Sf9-LD₅₀ and Sf9-LD₈₀.

While the so-called 'mode I 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 immunofluorescence 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-LD₈₀
<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /><figcaption>**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.</figcaption></figure>

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-LD₅₀ or Sf9-LD₈₀ 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

<figure><img src="image_10.png" /><figcaption>**Fig. 5.** Cry1Ca plasma membrane binding is not affected by the presence of chelators.

Sf9wt cells (left row) or Sf9-LD₈₀ 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.</figcaption></figure>

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 ([Ca²⁺]ᵢ). Both Ca²⁺ 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 Ca²⁺ binding sites. Increased intracellular Ca²⁺ 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 [Ca²⁺]ᵢ 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
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, Mg²⁺ could have a more central role than Ca²⁺ 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 Mg²⁺ 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 Mg²⁺-independent AC whose regulation requires Ca²⁺ (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]ᵢ in favour of a peculiar mode of action of this toxin (Fig. 6). Moreover, we observed (Fig. 7) that the Sf9-LD₈₀ resistant cell line was unable to elevate [cAMP]ᵢ in response to Cry1Ca treatment. Thus, the [cAMP]ᵢ 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 Ca²⁺ as another important molecule involved in the Cry1Ca mode of action, we might think that through toxin pore formation, Ca²⁺ entry in the cell would stimulate Ca²⁺-dependent events interacting with the PKA pathway. The protein kinase C (PKC) pathway could be such Ca²⁺-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 Ca²⁺ 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 Ca²⁺ 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-LD₈₀ cells, which were the most resistant to Cry1Ca toxicity, were affected in cAMP production. Furthermore, our results have shown the importance of Ca²⁺ and Mg²⁺ 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 Ca²⁺ 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.

## MATERIALS AND METHODS

### Insect cell lines culture and selection

Sf9 cell lines (ATCC® Number: CRL-1711TM; 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 UI/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 (LD₅₀), or at the lethal dose giving 80% dead cells (LD₈₀). 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 LD₅₀ or LD₈₀ toxin doses. Recovered resistant cells after the selection program were referred to as Sf9-LD₅₀ and Sf9-LD₈₀. 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)phHT320 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 Na₂CO₃/NaHCO₃, 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×10⁵ 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% (LD₅₀) or 80% (LD₈₀) dead cells were directly deduced from the curves. Resistance ratio (RR) due to treatments were obtained by the formula: $LD₅₀$ (Sf9+treatment) / $LD₅₀$ (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 T<sub>m</sub>: 46°C). The experimental design of each comparison (Sf9-LD₅₀ versus Sf9p40, Sf9-LD₈₀ 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-LD₅₀, specifically in Sf9-LD₈₀ 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 (–).

### Competing interests

The authors declare no competing or financial interests.

### 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.

### Supplementary information

Supplementary information available online at http://bio.biologists.org/lookup/doi/10.1242/bio.037085.supplemental

### References

* 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. doi:10.1002/ps.1798

* 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. doi:10.1128/AEM.68.9.4583-4592.2002

* 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. doi:10.1002/arch.10128

* 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. doi:10.1242/jcs.02440

* 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. doi:10.1042/BJ20101010

* 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. doi:10.1016/0005-2744(77)90155-3

* 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. doi:10.1016/j.ibmb.2016.07.008

* 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. doi:10.1074/mcp.M200069-MCP200
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. doi:10.1186/s12864-015-2240-7

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. doi:10.1007/s00122-002-0942-0

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. doi:10.3390/toxins7020516

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. doi:10.1007/s13205-018-1223-1

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. doi:10.1038/nrm1365

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. doi:10.1603/EC14190

Finkler, A., Ashery-Padan, R. and Fromm, H. (2007). CAMTAs: Calmodulin-binding transcription activators from plants to human. *FEBS Lett.* **581**, 3893-3899. doi:10.1016/j.febslet.2007.07.051

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. doi:10.1126/science.1060949

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. doi:10.1186/gb-2004-5-10-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. doi:10.1371/journal.pone.0165632.eCollection 2016

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. doi:10.1007/s002320001060

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. doi:10.1016/j.jip.2007.03.008

Hess, A. and Iyer, H. (2007). Fisher's combined p-value for detecting differentially expressed genes using Affymetrix expression arrays. *BMC Genomics* **8**, 96. doi:10.1186/1471-2164-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. doi:10.1074/jbc.M701314200

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. doi:10.1016/j.jip.2006.01.010

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. doi:10.1371/journal.pone.0017606

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. doi:10.1042/bj2530235

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. doi:10.1006/jipa.1997.4723

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. doi:10.1007/s00232-008-9106-8

Le Brigand, K. and Barbry, P. (2007). Mediante: a web-based microarray data manager. *Bioinformatics* **23**, 1304-1306. doi:10.1093/bioinformatics/btm106

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. doi:10.1016/j.ibmb.2006.05.009

Li, H., Rao, A. and Hogan, P. G. (2011). Interaction of calcineurin with substrates and targeting proteins. *Trends Cell Biol.* **21**, 91-103. doi:10.1016/j.tcb.2010.09.011

Litvin, T. N., Kamenetsky, M., Zarifyan, A., Buck, J. and Levin, L. R. (2003). Kinetic properties of "soluble" adenyllyl cyclase. *J. Biol. Chem.* **278**, 15922-15926. doi:10.1074/jbc.M212475200

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. doi:10.1603/0022-0493-93.1.1

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. doi:10.1128/AEM.67.7.3216-3219.2001

Luo, K., Lu, Y.-J. and Adang, M. J. (1996). A 106 kDa form of aminopeptidase is a receptor for Bacillus thuringiensis Cry1C δ-endotoxin in the brush border membrane of Manduca sexta. *Insect Biochem. Mol. Biol.* **26**, 783-791. doi:10.1016/S0965-1748(96)00027-6

Masson, L., Lu, Y.-J., Mazza, A., Brousseau, R. and Adang, M. J. (1995). The Cry1A(c) receptor purified from manduca sexta displays multiple specificities. *J. Biol. Chem.* **270**, 20309-20315. doi:10.1074/jbc.270.35.20309

Mochly-Rosen, D. (1995). Localization of protein kinases by anchoring proteins: a theme in signal transduction. *Science* **268**, 247-251. doi:10.1126/science.7716516

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. doi:10.1073/pnas.0831036100

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. doi:10.1186/s12864-015-1751-6

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. doi:10.1016/S0965-1748(01)00089-3

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. doi:10.3390/toxins6123296

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. doi:10.1007/s00232-001-0071-8

Pieroni, J. P., Harry, A., Chen, J., Jacobowitz, O., Magnusson, R. P. and Iyengar, R. (1995). Distinct characteristics of the basal activities of adenyllyl Cyclases 2 and z6. *J. Biol. Chem.* **270**, 21368-21373. doi:10.1074/jbc.270.36.21368

Pigott, C. R. and Ellar, D. J. (2007). Role of receptors in bacillus thuringiensis crystal toxin activity. *Microbiol. Mol. Biol. Rev.* **71**, 255-281. doi:10.1128/MMBR.00034-06

Portugal, L., Munoz-Garay, C., Martinez 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. doi:10.1016/j.ibmb.2016.11.004

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. doi:10.1023/A:1006946027619

Raymond, B., Johnston, P. R., Nielsen-LeRoux, C., Lereclus, D. and Crickmore, N. (2010). Bacillus thuringiensis: an impotent pathogen? *Trends Microbiol.* **18**, 189-194. doi:10.1016/j.tim.2010.02.006

Robinson-White, A. and Stratakis, C. A. (2002). Protein kinase a signaling. *Ann. N. Y. Acad. Sci.* **968**, 256-270. doi:10.1111/j.1749-6632.2002.tb04340.x

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. doi:10.1016/0005-2736(91)90237-3

Smyth, G. K. (2004). Linear models and empirical bayes methods for assessing differential expression in microarray experiments. *Stat. Appl. Genet. Mol. Biol.* **3**, Article3. doi:10.2202/1544-6115.1027

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. doi:10.1016/S0378-1097(00)00394-3

Sprang, S. R. (1997). G proteins, effectors and GAPs: structure and mechanism. *Curr. Opin. Struct. Biol.* **7**, 849-856. doi:10.1016/S0959-440X(97)80157-1

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. doi:10.1603/EC10040

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. doi:10.1002/ps.4740

Tesmer, J. J., Sunahara, R. K., Johnson, R. A., Gosselin, G., Gilman, A. G. and Sprang, S. R. (1999). Two-metal-ion catalysis in adenyllyl cyclase. *Science* **285**, 756-760. doi:10.1126/science.285.5428.756

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. doi:10.1016/j.jip.2012.05.001

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. doi:10.1016/S0005-2736(97)00184-3

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. doi:10.1038/sj.bjp.0703452

Wettenhall, J. M. and Smyth, G. K. (2004). limmaGUI: a graphical user interface for linear modeling of microarray data. *Bioinformatics* **20**, 3705-3706. doi:10.1093/bioinformatics/bth449

Willoughby, D. and Cooper, D. M. (2007). Organization and Ca 2+ regulation of adenyllyl cyclases in cAMP microdomains. *Physiol. Rev.* **87**, 965-1010. doi:10.1152/physrev.00049.2006

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. doi:10.1111/j.1472-765X.2005.01712.x

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. doi:10.1021/bi001321q

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. doi:10.1021/jf025863l

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. doi:10.1074/jbc.274.21.14909

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. doi:10.1016/j.bbrc.2007.08.002

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 adenyllyl cyclase/PKA signaling pathway is induced by the Cry1Ab toxin of Bacillus thuringiensis. *Proc. Natl. Acad. Sci. USA* **103**, 9897-9902. doi:10.1073/pnas.0604017103

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. doi:10.1074/jbc.M110057200