Skip to main content

ORIGINAL RESEARCH article

Front. Physiol., 29 October 2024
Sec. Invertebrate Physiology
This article is part of the Research Topic Insect Physiology Aspects of Environmentally Friendly Strategies for Crop Pests and Invertebrate Vectors Control, Volume II View all 8 articles

Cry1Ac toxin binding in the velvetbean caterpillar Anticarsia gemmatalis: study of midgut aminopeptidases N

  • 1Departamento de Bioquímica, Instituto de Química, Universidade Federal do Rio de Janeiro, Rio de Janeiro, Brazil
  • 2Instituto de Biofísica Carlos Chagas Filho, Universidade Federal do Rio de Janeiro, Rio de Janeiro, Brazil
  • 3Instituto Nacional de Entomologia Molecular, Rio de Janeiro, Brazil

The velvetbean caterpillar Anticarsia gemmatalis is one of the main soybean defoliators in Brazil. Currently, the main biopesticide used to control insect pests worldwide is the bacteria Bacillus thuringiensis (Bt), which produces entomopathogenic Crystal toxins (Cry) that act in the midgut of susceptible insects, leading them to death. The mode of action of Cry toxins in the midgut involves binding to specific receptors present on the brush border of epithelial cells such as aminopeptidase N (APN), alkaline phosphatase (ALP), cadherin, and others. Mutations in these receptors, among other factors, may be involved in the development of resistance; identification of functional Cry receptors in the midgut of A. gemmatalis is crucial to develop effective strategies to overcome this possible scenario. This study’s goal is to characterize APNs of A. gemmatalis and identify a receptor for Cry1Ac in the midgut. The interaction of Bt spores with the midgut epithelium was observed in situ by immunohistochemistry and total aminopeptidase activity was estimated in brush border membrane vesicle (BBMV) samples, presenting higher activity in challenged individuals than in control ones. Ten APN sequences were found in a A. gemmatalis’ transcriptome and subjected to different in silico analysis, such as phylogenetic tree, multiple sequence alignment and identification of signal peptide, activity domains and GPI-anchor signal. BBMV proteins from 5th instar larvae were submitted to a ligand blotting using activated Cry1Ac toxin and a commercial anti-Cry polyclonal antibody; corresponding bands of proteins that showed binding to Cry toxin were excised from the SDS-PAGE gel and subjected to mass spectrometry analysis, which resulted in the identification of seven of those APNs. Quantitative PCR was realized to compare expression levels between individuals subjected to sublethal infection with Bt spores and control ones, presenting up- and downregulations upon Bt infection. From these results, we can infer that aminopeptidases N in A. gemmatalis could be involved in the mode of action of Cry toxins in its larval stage.

1 Introduction

Bacillus thuringiensis (Berliner, 1915; Bacillales: Bacillaceae) (Bt) is a gram-positive, spore-forming bacteria known for its entomopathogenic effect against susceptible insects (Bravo et al., 2007; Sanchis, 2011). During its sporulation phase of growth, it produces insecticidal proteins as crystal inclusions (Cry toxins) that are pathogenic to several insect models, including known insect pests from the orders Lepidoptera, Coleoptera, and Hymenoptera (Bravo et al., 2007). For this reason, these proteins have been used worldwide in insecticidal sprays and, in later years, incorporated in genetically engineered crops (Horikoshi et al., 2021b; Pozebon et al., 2020; Bernardi et al., 2012).

Soybean (Glycine max (L.) Merrill, 1917) is one of the major agricultural commodities worldwide, with great importance in the global market, with projections that indicate that the brazilian share of global trade could increase to over 60% until 2033 (dos Reis et al., 2020; Valdes et al., 2023). Brazil is currently the main producer and exporter of soybean (CONAB, 2023; FAO, 2024), with 154.6 million tons of grains produced in the 2022/23 harvest (CONAB, 2023). The velvetbean caterpillar Anticarsia gemmatalis (Hübner, 1818; Lepidoptera: Erebidae) is one of the soybean’s main defoliators in the Americas (Hoffman-Campo et al., 2000; Bernardi et al., 2012), causing great damage to the production of this grain. While chemical pesticides have typically been used in A. gemmatalis control, Bt-based biopesticides have been thoroughly used in the management of this species’ populations in the field (Knaak and Fiuza, 2005; Fernandes et al., 2021). In recent years, genetically modified soybean cultivars expressing insecticidal Cry proteins from B. thuringiensis (Bt crops) are quickly becoming a key tool in the management of pests, after its commercial availability (Pozebon et al., 2020; Bernardi et al., 2012; Chattopadhyay and Banerjee, 2018). Bt soybean MON 87701 × MON 89788 (Intacta RR2 PRO®) provides protection against A. gemmatalis in a high-dose manner, presenting high levels of Cry1Ac expression throughout the planting season conferring complete neonate mortality and effectively managing their populations in the field (Tabashnik et al., 2023; Horikoshi et al., 2021a; Bernardi et al., 2012).

Once ingested, Bt crystals are solubilized in the insect midgut alkaline environment, originating Cry protoxins, which are further activated by midgut proteases. Upon activation, Cry toxins bind to specific receptors in the midgut brush border membranes of epithelial cells and undergo conformational changes that cause the formation of oligomers, which are then inserted into the cell membrane. Cry toxin oligomerization at the plasma membrane forms pores, which leads to osmotic imbalance and cell lysis, which ultimately leads to insect death (Bravo et al., 2005; Adang et al., 2014; Palma et al., 2014; Melo et al., 2014; Liu et al., 2021). Four classes of proteins seem to be the major gut receptors for Cry toxins: aminopeptidase N (APN), alkaline phosphatase (ALP), cadherin (CAD), and ATP binding cassette subfamily C member 2 (ABCC2) transporters (Endo, 2022; Soberón et al., 2018; Bravo et al., 2004; Bravo et al., 2007; Gómez et al., 2002; Gómez et al., 2006; Pacheco et al., 2009). After activation, Cry toxins bind to APNs, abundant in lipid rafts in the membrane, which promotes the localization and concentration of activated toxins in these regions (Xu et al., 2014; Liu et al., 2021). The abundance of APN and its lower affinity for Cry toxins compared to receptors like CAD may allow it to act as a “toxin sink” that concentrates Cry proteins at the midgut membrane surface. APN can then pass the toxin to other receptors while also facilitating insertion of the oligomeric pore complex into the membrane. Accordingly, a decrease in total aminopeptidase activity was observed in strains resistant to Cry toxins (Yang et al., 2010). The APN (EC.3.4.11.2) consists of a class of metalloproteases that act in the midgut brush border of insect larvae cleaving N-terminal amino acids from peptides during digestion (Adang, 2004). Knight et al. (1994) first identified an APN as a Cry toxin-binding protein and putative receptor in Manduca sexta and since then several other works have identified these proteins as Cry receptors in different insect species (Adang, 2004; Pigott and Ellar, 2007; Liu et al., 2021). APNs present features such as aminopeptidase motif “GAMENWG,” Zn++-binding motif “HEXXHX18E,” a signal peptide in the N-terminal end and GPI-anchor peptide in the C-terminal end, which facilitates their attachment to the brush border, along with several O- and N-glycosylation sites (Pigott and Ellar, 2007; Adang, 2004). Based on their amino acid sequence similarity, they are classified into 13 clusters in Lepidoptera, with several isoforms expressed in the midgut (Guo et al., 2020; Crava et al., 2010; Hughes, 2014; Lin et al., 2014). The identification of lepidopteran APNs has allowed studies of the role of these proteins in the mode of action of Cry toxins, demonstrating their involvement in the pathogenesis of Bt toxins, mainly as functional receptors for these toxins (Gill and Ellar, 2002; Rajagopal et al., 2003; Wang et al., 2005).

Despite the importance of the velvetbean caterpillar (A. gemmatalis) as a major soybean pest, limited research has been conducted on the functional receptors for Cry toxins in this species. Previous studies have demonstrated the binding of various Cry toxins to A. gemmatalis’ midgut brush border membrane vesicles (BBMV) through competition-binding assays (Bel et al., 2017; Bel et al., 2019) and to the midgut epithelial tissue using biotinylated toxins (Fiuza et al., 2013). Additionally, a membrane-associated alkaline phosphatase (ALP) characterized in the midgut of A. gemmatalis showed interaction with Cry1Ac toxin in vitro through enzyme-linked immunosorbent assay (da Silva et al., 2019), suggesting its potential role as a Cry toxin ligand. Even though APNs are commonly described and studied as potential Cry receptors in other lepidopteran species, no works described these proteins in A. gemmatalis. In this study, we characterized A. gemmatalis’ aminopeptidases (AgAPNs) and identified those that positively bound to Cry toxin in vitro; analyses of AgAPN gene expression after Bt spore feeding bioassays demonstrated that some of these proteins could be involved in Cry’s mode of action, due to changes in expression in this condition. These findings will contribute to future works, aiding the understanding of potential APN role in Cry toxin mode of action in A. gemmatalis.

2 Materials and methods

2.1 Insects

A colony of A. gemmatalis was established in the lab using eggs obtained from EMBRAPA SOJA, Londrina, PR, Brazil. Larvae were reared on an artificial diet previously described by Hoffmann-Campo et al. (1985) and maintained under 25°C ± 3°C, 70% ± 10% humidity, and 14:10 h (light/dark) photoperiod.

2.2 Immunohistochemistry

For the immunohistochemistry assay, 0.5 mg/mL of Bt spores was mixed in the artificial diet fed to 4th-5th instar larvae; the control group was fed in an artificial diet mixed with distilled water. Samples were obtained after 12 h of Bt exposure and controls without exposure to the spores were also monitored under the same conditions. Challenged and control larvae were dissected and the gut of each was washed in phosphate-buffered saline [PBS; 137 mM NaCl, 2.7 mM KCl, 10 mM phosphate (pH 7.4)] to remove any unbound material before sample fixation. Tissues were fixed in 4% formaldehyde, 0.1% glutaraldehyde in 0.1 M sodium cacodylate buffer (pH 7.2) (Gomes et al., 2013) and stored in 4°C until use, for no longer than a week. After fixation, tissues were washed in 0.1 M sodium cacodylate buffer and incubated in blocking buffer (2% bovine serum albumin, 0.3% Triton X-100, PBS) for 2 h, followed by incubation in 1:250 commercial polyclonal primary anti-B. thuringiensis Cry1Ab Toxin antibody (Abcam Inc.; catalog number #ab51586) for 2 h in blocking buffer. After being washed in the blocking buffer, the tissues were incubated in 1:500 anti-rabbit Alexa-488-conjugated secondary antibody for 2 h, washed and incubated in 0.1 μg/mL DAPI. Whole gut samples were observed in a Zeiss 910 LSM confocal microscope.

2.3 In silico analysis of APN sequences

Raw RNA-Seq data of A. gemmatalis was obtained from a previously published study by Pezenti et al. (2021) (SRA accession number: PRJNA387150). The transcriptome was reassembled following the pipeline described in the original publication. APN amino acid sequences were identified through eggNOG-mapper (Cantalapiedra et al., 2021) and aligned using the Clustal W server (https://www.ebi.ac.uk/jdispatcher/msa/clustalo). For the phylogenetic tree, recovered APN sequences from nine lepidopteran species (M. sexta, Heliothis virescens, Spodoptera frugiperda, Spodoptera litura, Helicoverpa armigera, Trichoplusia ni, Bombyx mori, Plutella xylostella, and Ostrinia furnacalis), along with five sequences of Homo sapiens’ APNs to form an outgroup (Supplementary Table S1). Alignment was conducted with Clustal W on MEGA 11 (Tamura et al., 2021) and the phylogenetic analysis was performed in RAxML (Random Axelerated Maximum Likelihood) (Stamatakis, 2014) software using the maximum likelihood (ML) method with a bootstrapping procedure with 1,000 replicates. The tree was visualized by FigTree (http://tree.bio.ed.ac.uk/software/figtree/). The presence and location of signal peptide cleavage sites in the APN amino acid sequences were determined using the SignalP 5.0 Server (https://services.healthtech.dtu.dk/services/SignalP-5.0/). The presence of glycosylphosphatidylinositol (GPI) anchors was predicted using the PredGPI software (https://busca.biocomp.unibo.it/predgpi/). Potential N-linked and O-linked glycosylation sites were analyzed using the NetNGlyc1.0 (https://services.healthtech.dtu.dk/services/NetNGlyc-1.0/) and NetOGlyc4.0 (https://services.healthtech.dtu.dk/services/NetOGlyc-4.0/) programs, respectively. Schematic representations of the APN sequences, including the locations of signal peptide cleavage sites, GPI anchors, and glycosylation sites, were generated using the Illustrator for Biological Sequences (IBS, 2022) software (https://ibs.renlab.org/#/home).

2.4 Preparation of midgut brush border membrane vesicles (BBMVs)

Midguts from 5th instar larvae of A. gemmatalis were longitudinally dissected and washed in saline solution (NaCl 0.15 M). The samples were stored in ice-cold SET buffer (0.15 M sucrose, 17 mM tris[hydroxymethyl]aminomethane (Tris), 5 mM ethylene glycol-bis (β-aminoethyl ether)-N,N,N′,N′-tetraacetic acid (EGTA); pH 7.5) and kept at −20°C until use, for no longer than a week. BBMVs were prepared following the MgCl2 differential precipitation method described in Wolfersberger et al. (1987) with some modifications. Briefly, the dissected midgut-epithelium samples were homogenized in SET buffer with a glass potter, and an equal volume of 24 mM MgCl2 was added to the homogenate. Samples were incubated on ice for 15 min and then centrifuged at 3,000 g for 15 min at 4°C. The pellet was discarded, and the supernatant was centrifuged at 36,603 g for 30 min at 4°C. The pellet was resuspended in SET buffer and the centrifugation cycle was repeated. The final pellet was resuspended in 4-(2-hydroxyethyl) piperazine-1-ethane-sulfonic acid (HEPES) buffer (50 mM; pH 7.2). This resuspended pellet (P3) containing the proteins from BBMVs was quantified for total protein content using Pierce™ 660 nm Protein Assay (Thermo Scientific); bovine serum albumin (BSA; Sigma Chemical Company) was used as standard. The BBMV protein profile was observed by (12%) sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) (Laemmli, 1970).

2.5 Total leucine aminopeptidase activity in midgut BBMV samples

For total leucine aminopeptidase activity assay, a sublethal concentration of Bt spores (0.1323 mg/mL) was mixed in the artificial diet fed to 4th-5th instar larvae; the control group was fed in the artificial diet mixed with distilled water. This sublethal concentration was estimated in a bioassay that determined lethal concentrations of Bt spores to A. gemmatalis colony kept in the lab (data not shown). After 48 h of exposure to Bt, midguts were dissected and pooled (4 midguts per sample), with three biological replicates performed for the assay. Samples were then subjected to the preparation of BBMVs protocol and used in the total aminopeptidase activity assay. Total protein concentration was determined using the Pierce™ 660 nm Protein Assay Kit (Thermo Scientific) according to the manufacturer’s instructions, with bovine serum albumin (BSA) as the standard. Samples were diluted to a final concentration of 0.05 mg/mL in 0.1 M Tris-HCl buffer (pH 8.6) for the measurement of total aminopeptidase activity. Total aminopeptidase activity was assessed using the chromogenic substrate leucine-p-nitroanilide (LpNA) (Sigma, St. Louis, MO), as previously described by Valaitis et al. (1999) and Erlanger et al. (1961). Briefly, 5 µg of each diluted sample was mixed with 0.5 mM LpNA in 0.1 M Tris-HCl buffer (pH 8.6) in a microplate well. The enzymatic reaction was monitored by measuring the increase in optical absorbance at 405 nm using a SpectraMax® M2e microplate reader (Molecular Devices). Absorbance readings were taken every 30 s for a total of 15 min.

The mean velocity (Mean V) was calculated as an increase of absorbance per min in the linear portion of initial velocity of the enzymatic reaction using the GraphPad Prism 8.0. One unit of the enzyme activity was defined as the amount of enzyme that hydrolyzed 1 µmol of substrate to chromogenic product per min. Specific aminopeptidase activity calculations were performed based on the extinction coefficient (9.9 mM−1 cm−1) for p-nitroaniline (Hafkenscheid, 1984). APN activity was estimated from 3 reactions (replicates) for each condition (control and challenged) and sample type. Specific APN activities were presented as means and standard errors of the mean (SEM). Two-way ANOVA was used to determine differences in APN activity between the two conditions (control and challenged) for each sample.

2.6 Cry toxin purification and activation

B. thuringiensis subesp. kurstaki LFB-FIOCRUZ 475 (Bt) was kindly provided by Leon Rabinovitch (CCGB/IOC/Fiocruz/MS). This strain was used to produce Cry toxins, following Coleção de Culturas do Gênero Bacillus e Gêneros Correlatos (CCGB, Fundação Oswaldo Cruz) instructions. Bt spores were grown during 48 h in Nutrient Broth medium (Himedia) at 200 rpm at 37°C. This pre-inoculum was transferred to a new volume of nutrient broth and was shaken in the same growth conditions for 5 days until complete sporulation, as recommended instruction. After that, the culture was centrifuged at 7,500 g for 30 min at 4°C and the pellet was collected as spores/crystals-enriched fraction. Crystalline inclusions were solubilized by incubating the pellet in a 50 mM sodium carbonate buffer (pH 9.6) containing 0.1% 2-mercaptoethanol (SC buffer) under agitation for 2 h at room temperature. The mixture was centrifuged at 16,000 g for 30 min at 4°C. The supernatant was labeled as whole Cry1Ac toxin (WCT; about 130 kDa) and was either stored or further activated. WCT activation was performed by overnight incubation at 4°C with bovine pancreas trypsin (2 mg/mL). Following incubation, activated WCT was concentrated in an Amicon 30,000 MWCO spin filter for the removal of reminiscent trypsin. The molecular weight and the toxin integrity were checked by 10% SDS-PAGE (data not shown). This sample was named Cry1Ac toxin and used in the toxin overlay assay described below.

2.7 Toxin overlay assay in midgut BBMV proteins

The midgut BBMV proteins (50 µg) were separated by 12% SDS-PAGE (Laemmli, 1970), transferred to a nitrocellulose membrane, and Ponceau S was used to confirm protein transfer (Towbin et al., 1979). Following, the membrane was incubated in Blocking Buffer [5% (w/v) skimmed milk powder in Tris buffered-saline (TBS, pH 7.4) containing 0.1% Tween-20] overnight. For ligand blotting analysis (Toxin Overlay Assay - TOA), the membrane was incubated in 5 mL of TOA Blocking Buffer [3% (w/v) BSA in Tris buffered-saline (TBS, pH 7.4) containing 0.1% Tween-20] and containing WCT (Cry1Ac toxin) for 7 h at room temperature to interact with respective BBMV proteins. BSA was used as a negative control. The unbound toxins were removed by washing the blots with 5 mL of TOA Blocking Buffer, three times. Following, this blot was individually incubated with a polyclonal anti-B. thuringiensis Cry1Ab Toxin antibody (Abcam Inc.) (1:5,000 dilution in TOA Blocking Buffer) overnight; the section containing the activated Cry1Ac toxin was cut off from the membrane and incubated only with blocking buffer as a negative control. The membrane was then incubated with an anti-rabbit alkaline phosphatase-conjugated antibody (Sigma Aldrich). Finally, the blot was developed using NBT/BCIP (0.5 mM nitro blue tetrazolium, 0.57 mM 5-bromo-4-chloro-3-indolyl-phosphate) substrate in alkaline phosphatase buffer (10 mM Tris [pH 9,6], 100 mM NaCl, 5 mM MgCl2).

2.8 Sample preparation for LC-MS/MS analysis

Protein samples for mass spectrometry analysis were obtained from a 12% SDS-PAGE gel. Two bands corresponding to those that demonstrated binding to Cry1Ac in the ligand blot were excised from the gel. Coomassie R stain was removed from the excised bands. Reduction was performed by incubating the samples with 10 mM dithiothreitol (DTT) in 50 mM ammonium bicarbonate for 1 h at 30°C. Following reduction, alkylation was carried out by adding 40 mM iodoacetamide and incubating the samples for 30 min in the dark. After alkylation, the samples were washed with Milli-Q® water and dehydrated using 90% acetonitrile. The samples were then vacuum-dried to remove all traces of acetonitrile. Protein digestion was performed using commercial trypsin (Promega) at a 1:50 (w/w) enzyme-to-protein ratio. The samples were incubated with trypsin overnight at 37°C. Following digestion, 0.1% trifluoroacetic acid (TFA) was added to the samples, and they were incubated at room temperature for 1 h.

Tryptic peptides were purified using manual reversed-phase chromatography with Poros 50 R2 resin (PerSeptive Biosystems) and subsequently vacuum-dried. The peptides were then solubilized in 0.1% formic acid for injection into the mass spectrometer. Samples were analyzed using a nanoLC-MS/MS system consisting of an EASY II-nano LC (Proxeon Biosystem) coupled to a Q-Exactive Plus (Thermo Fisher Scientific) mass spectrometer. A total of 1 µg of peptides was loaded onto a manually packed pre-column (100 μm × 2 cm) containing C-18 ReproSil 5 μm resin (Dr. Maisch) and then separated on a 20 cm analytical column packed with Reprosil-pur C18-AQ 3 μm resin (Dr. Maisch). Chromatographic separation was performed using a linear gradient of solution B (95% ACN, 0.1% TFA) from 5% to 20 s, followed by an increase to 40% over 8 min, then to 95% over 4 min, and finally maintained at 95% for 8 min. The flow rate was set to 250 nL/min. Mass spectra % over 40 min were acquired in positive mode using a Top 10 data-dependent acquisition (DDA) method. MS1 scans were acquired in the Orbitrap analyzer with a mass range of 350–1,800 m/z, a resolution of 60,000 (at m/z 400), a minimum signal threshold of 10,000, and an isolation window of 2.0. The 10 most intense ions were selected for fragmentation by collision-induced dissociation (CID) with a normalized collision energy of 30 and a dynamic exclusion time of 30 s.

2.9 LC-MS/MS data analysis for protein identification

The spectra obtained from the LC-MS/MS analyses were processed using Proteome Discoverer 2.1 Software (Thermo Scientific) with the Sequest HT search engine against an A. gemmatalis protein database. This database was generated through a de novo assembly of the A. gemmatalis transcriptome based on raw data from a published study by Pezenti et al. (2021) (SRA accession number: PRJNA387150). For the search, the following parameters were used: precursor tolerance of 10 ppm, fragment tolerance of 0.1 Da, tryptic cleavage specificity, two maximum missed cleavage sites allowed, fixed modification of carbamidomethyl (Cys) and oxidation (Met). Peptides with high confidence were selected, and only identifications with q values equal to or less than 0.01 FDR were considered.

2.10 Gene expression analysis

Samples for analysis of gene expression were obtained in a feeding bioassay, in which a sublethal concentration of Bt spores (0.1323 mg/mL) was mixed in the artificial diet fed to 4th-5th instar larvae; the control group was fed in the artificial diet mixed with distilled water. After 24, 48 and 72 h of exposure to Bt, midguts were dissected and pooled (4 midguts per sample), with five biological replicates performed for the assay. Specimens were dissected and midgut epithelia were cleaned and used for total RNA extraction with TRI Reagent® (Sigma-Aldrich). Total RNA was quantified using NanoDrop One/Onec Microvolume UV-Vis Spectrophotometer (ThermoFisher Scientific), following the manufacturer’s recommendations. Total RNA (500 ng) was reverse transcribed using GoScript™ Reverse Transcriptase kit (Promega) following the manufacturer’s recommendations. For cDNA amplification, primers were designed following the re-annotation of a previously published A. gemmatalis transcriptome (Pezenti et al., 2021) (Table 1).

Table 1
www.frontiersin.org

Table 1. Aminopeptidase N contigs from A. gemmatalis identified in the transcriptome (Pezenti et al., 2021).

Primers were validated by RT-PCR using GoTaq Hot Start Colorless Master Mix (Promega) and the Veriti® 96-Well Thermal Cycler (Applied Biosystems). The cycling conditions were as follows: 94°C for 2 min, followed by 40 cycles of 94°C for 30 s, 52°C for 30 s, 72°C for 1 min. This stage was followed by 72°C for 10 min. Resulting amplification was visualized in 2% agarose gels using GelRed® Nucleic Acid Gel Stain (Biotium) following manufacturer’s recommendations.

Quantitative real-time PCR (qRT-PCR) was performed using qPCRBIO SyGreen Mix (PCR Biosystems) and the ViiA™ 7 Real-Time PCR System (Applied Biosystems). The qRT-PCR cycling conditions were as follows: 50°C for 2 min, 95°C for 10 min, followed by 40 cycles of 95°C for 30 s, 60°C for 30 s, and 72°C for 30 s. Relative expression levels were analyzed using the 2−ΔΔCT method, using the expression level of the housekeeping gene elongation factor 1-alpha (ELF-1α) as a reference gene normalizer.

2.11 Statistical analysis

The results were submitted to one- and two-Way ANOVA (followed by Sidak’s multiple comparisons test), using the GraphPad Prism 8 software. Differences were considered significant at p < 0.05.

3 Results

3.1 Cry1Ac binds to epithelial cells in the midgut of Anticarsia gemmatalis

To verify the immobilization of Cry toxins in the A. gemmatalis gut, we conducted an experiment where 5th instar A. gemmatalis larvae were fed with a sublethal Bt spore dose (0.5 mg/mL) for 12 h. After incubation with anti-B. thuringiensis Cry1Ab polyclonal antibody, Bt toxin was identified as a punctate pattern along the whole midgut epithelia (Figures 1B, C), suggesting a possible interaction with the apical portion of epithelial cells. Fluorescence was specific to the midgut and not observed in negative controls (Figure 1A). These findings suggest that Cry toxins interact with A. gemmatalis’ brush border membrane proteins, where APNs are known to be localized.

Figure 1
www.frontiersin.org

Figure 1. Immunofluorescence detection of Cry toxins in A. gemmatalis midgut tissues. Midgut samples were collected from larvae subjected to a feeding bioassay. The challenged group was fed an artificial diet containing Bacillus thuringiensis (Bt) spores at a concentration of 0.5 mg/mL, while the control group received only the artificial diet without spores. After 12 h of exposure, larval midguts were dissected, fixed, and incubated with an anti-Bacillus thuringiensis Cry1Ab Toxin antibody (1:250 dilution). Detection was performed using a confocal microscope following incubation with an Alexa Fluor 488-conjugated anti-rabbit secondary antibody. Nuclei were counterstained with DAPI for visualization. (A) Midgut from the control group showing no detectable signal from Cry antibodies at the brush border, included as a reference. (B, C) Detection of Cry toxins (green dots) at the apical brush border of the midgut epithelium using the Cry antibody. Images are representative of ten midguts.

3.2 In silico analysis of aminopeptidases N from Anticarsia gemmatalis

To identify putative APN expressed by A. gemmatalis, we analyzed a previously published A. gemmatalis whole-body RNA-Seq dataset (Pezenti et al., 2021). The analysis of the de novo assembly of the transcriptome data resulted in the identification of 10 APN contigs, which are listed in Table 1. Of those, eight were full length sequences, hence being used in all posterior analysis; aminopeptidase N13 and aminopeptidase N-like isoform X2 were left out of most analyses since the sequences were only partial, impairing the results.

Alignment of deduced translated amino acid sequences of AgAPNs through multiple sequence alignment and the search for the presence of characteristic gluzincin aminopeptidase activity motif “GAMENWG” and the consensus zinc-binding motif “HEXXHX18H” showed the presence of both motifs in nearly all sequences (Figure 2A).

Figure 2
www.frontiersin.org

Figure 2. In silico identification and analysis of AgAPNs. AgAPN sequences were identified from a de novo assembled A. gemmatalis transcriptome and aligned against APN homologs. (A) Multiple sequence alignment of AgAPN sequences, partial. Red boxes indicate aminopeptidase motif “GAMENWG” and Zn++ binding motif “HEXXHX18E,” respectively. Alignment was conducted in Clustal Omega server, and the results were visualized as colored classification (percentage identity) using Jalview software (version 2.11.3.3). (B) Schematic representation of A. gemmatalis’ APNs. Green boxes correspond to the signal peptide; yellow boxes correspond to the aminopeptidase motif; red boxes correspond to the Zn++ binding motif; blue circles correspond to the GPI-anchor signal site.

The predicted N-terminal signal peptide cleavage site location was present in most of the sequences, varying between them; no signal peptide was predicted for AgAPN11 (Figure 2B). Further, a putative C-terminal GPI-anchor signal was predicted for AgAPN3, AgAPN4, AgAPN8, AgAPN6 and AgAPN2. No GPI-anchor site was predicted for AgAPN5, AgAPN10 and AgAPN11. Several potential O-glycosylation sites and one N-glycosylation site were identified in all sequences.

For the phylogenetic tree, APN proteins from different Lepidoptera species, belonging to the thirteen APN classes (Guo et al., 2020), were downloaded from NCBI and UniProt databases (Supplementary Table S1) to construct the tree using the maximum-likelihood method (Figure 3). This type of analysis has been broadly utilized to demonstrate the clustering of APN genes into different clusters, with high evolutionary conservation in each class (Guo et al., 2020; Wang et al., 2023). Following the genome-wide unified nomenclature and classification of APN genes in lepidopoteran insects (Guo et al., 2020), we designated APN genes in A. gemmatalis as AgAPN2, AgAPN3, AgAPN4, AgAPN5, AgAPN6, AgAPN8, AgAPN10 and AgAPN11.

Figure 3
www.frontiersin.org

Figure 3. Phylogenetic tree of AgAPN sequences. Analysis conducted in RAxML software, using the maximum likelihood method with 1,000 bootstraps. Five Homo sapiens’ APNs were used as an outgroup and APNs from nine lepidopteran species (Bombyx mori, Helicoverpa armigera, Heliothis virescens, Manduca sexta, Ostrinia furnacalis, Plutella xylostella, Spodoptera frugiperda, Spodoptera litura and Trichoplusia ni) were used to conduct the analysis (Supplementary Table S1). Eigth APNs of A. gemmatalis were indicated by bold.

3.3 Total aminopeptidase activity in the midgut

To evaluate aminopeptidase microvilli association and Bt modulation of aminopeptidase activity, we obtained fractionated samples from the midgut BBMVs of challenged and control individuals. Accordingly, we observed a Bt-upregulation of the aminopeptidase activity associated with the microvilli fraction (Figure 4).

Figure 4
www.frontiersin.org

Figure 4. Total aminopeptidase N activity of BBMV preparation samples of A. gemmatalis larvae. 5th instar larvae were subjected to a sublethal concentration of Bt spores (0.1323 mg/mL) for 48 h, when midguts were dissected and subjected to the preparation of midgut brush border membrane vesicles protocol. Control groups were kept unchallenged. A higher aminopeptidase activity was expected in the BBMV fraction, where aminopeptidases N were expected to be detected. S1, supernatant 1, of the mild centrifugation; S2, supernatant 2, of the first ultracentrifugation; S3, supernatant 3, of the second ultracentrifugation; P3, pellet 3, of the second ultracentrifugation, ressuspended in HEPES. Bars represent the means and standard deviation of 3 biological replicates (pool of 4 midguts). Asterisks indicate significant difference (P < 0.0001, two-way ANOVA).

3.4 Identification of Cry1Ac-binding proteins by LC-MS/MS

To identify Cry binding proteins in A. gemmatalis, we conducted a toxin overlay assay using Cry1Ac activated toxin on SDS-PAGE separated midgut BBMV proteins. The results presented in Figure 4 clearly show that Cry1Ac toxin interacted with membrane proteins, specifically with two bands of approximately 110-kDa and 65-kDa (Figure 5).

Figure 5
www.frontiersin.org

Figure 5. Toxin overlay assay of Cry1Ac toxin and BBMV proteins of A. gemmatalis. Midguts from 5th instar larvae were dissected (N = 200) and subjected to the preparation of midgut brush border membrane vesicles protocol. Proteins were separated by 12% SDS-PAGE and transferred to a nitrocellulose membrane; after blocking, the membrane was incubated with TOA Blocking Buffer containing activated Cry1Ac toxin, following incubation with Cry antibody. (A) Nitrocellulose membrane stained with Ponceau S; (B) membrane after blotting. Arrows indicate bands in P3 that present a positive result to the interaction with anti-Cry antibodies; Bt+ positively interacted with the antibody, as a positive control. P, protein standard; P3, pellet 3, of the second ultracentrifugation; BSA, bovine serum albumin (negative control); Bt+, activated Cry toxin (positive control); Bt−, activated Cry toxin not incubated with Cry antibody (negative control).

A new electrophoresis was performed from the same BBMV sample and bands in the same positions where there was binding to Cry1Ac in the TOA were excised and analyzed together by mass spectrometry. Results identified 1,152 proteins (Supplementary Table S1); the samples for this assay contained all proteins present in the two excised bands. Among them, seven AgAPNs were found (Table 2), suggesting that they could act as receptors for this toxin in A. gemmatalis’ midgut. The high number of unique peptides increases reliability in these results. It was possible to identify an alkaline phosphatase, which can also act as a receptor in the midgut, with a similar role in the progression of the infection.

Table 2
www.frontiersin.org

Table 2. A. gemmatalis’ midgut proteins that demonstrated binding to Cry1Ac identified by mass spectrometry.

3.5 APN expression analysis after exposure to Bt spores

To evaluate differences in expression of APNs between control individuals and challenged with Bt spores individuals, we performed RT-qPCR analysis of all AgAPNs described in this work upon 24, 48 and 72 h of exposure (Figure 6). Regarding the different conditions of exposure, AgAPN10, AgAPN4, AgAPN5, AgAPN6 and AgAPN11 did not present any significant difference in expression between the challenged and control groups in any of the time points. Genes of AgAPN2 and AgAPN3 of the challenged group presented an upregulation in 24 h of exposure in comparison with the control group, with a decrease in expression in the following time points and no significant difference from the control group. AgAPN8, on the other hand, presented downregulation of expression from 48 to 72 h in comparison with the control group. AgAPN5 and AgAPN6 did not present any significant difference of expression between treatments in each time point but did present a tendency downregulation in the challenged group from 48 to 72 h in relation to the initial period of exposure.

Figure 6
www.frontiersin.org

Figure 6. Expression analysis of AgAPNs in A. gemmatalis larvae midguts. Samples were obtained from midguts of 5th instar larvae subjected to a sublethal concentration of Bt spores in the challenged group, replaced by water in the control group. Midguts were dissected after 24, 48 and 72 h of exposure, following RNA extraction. Each graph shows the relative expression levels of AgAPNs in the midgut of A. gemmatalis’ larvae in each condition. (A) AgAPN2, (B) AgAPN3, (C) AgAPN4, (D) AgAPN5, (E) AgAPN6, (F) AgAPN8, (G) AgAPN10, (H) AgAPN11. Relative expression calculated through the 2−ΔΔCt method using ELF-1α as a reference gene. Gray bars correspond to control individuals and black bars correspond to challenged individuals. The analysis was conducted with five biological replicates, each consisting in a pool of four midguts, and two technical replicates. Asterisks indicate significant difference (*p < 0.02; **p < 0.005; ***p < 0.0005, one-way ANOVA).

4 Discussion

B. thuringiensis (Bt) is the most widely used bacteria for control of insect pests due to its production of parasporal entomopathogenic crystals during sporulation, known as Cry toxins (Liu et al., 2021). These are toxic by ingestion, acting in the midgut of susceptible insects where, after solubilization and activation, the Cry proteins oligomerize and form pores, which are then inserted in the epithelial membrane. This disrupts the cellular membrane, causing leakage of cellular contents and osmotic imbalance and allowing gut bacteria to invade the cells and the hemolymph, leading to septicemia; this ultimately triggers the infected insect to stop eating, alongside with the infection events that culminates in insect death (Liu et al., 2021). This infection triggers host responses, modulating different pathways and molecules, to counteract these effects (Pinos et al., 2021).

The binding of Cry1Ac toxin to brush border membrane vesicles (BBMV) from A. gemmatalis midguts was demonstrated by the fractionation of total aminopeptidase activity, revealing increased activity in BBMV fractions of larvae fed a diet containing Bt spores. This rise in aminopeptidase activity is consistent with previous reports that describe elevated APN activity following Cry toxin exposure (Ingle et al., 2001). These findings were further supported by immunohistochemistry, which confirmed the interaction of Cry1Ac with the midgut epithelium, as previously shown in other lepidopterans (Chauhan et al., 2021; Valaitis, 2011; Chen et al., 2005; Aimanova, Zhuang, and Gill, 2006).

Despite its importance as a major soybean defoliator, few studies addressed the molecular basis of Cry toxins’ mode of action on the midgut of A. gemmatalis. Bel et al. (2017), Bel et al. (2019) and Fiuza et al. (2013) identified Cry toxin binding to the midgut brush border tissue but did not identify a functional receptor for these toxins. Da Silva et al. (2019) characterized a membrane-associated alkaline phosphatase that binds to Cry1Ac toxin, but further studies would be required to describe it as a functional receptor. APNs have been extensively described as Cry toxin receptors in many insect pest species; these proteins consist in a group of glycoproteins attached to the midgut epithelial membrane through a GPI-anchor (Sato, 2002; Pardo-Lopez et al., 2013; Flores-Escobar et al., 2013; Herrero et al., 2016; Liu et al., 2021).

Ten APNs were identified from the transcriptome data of A. gemmatalis (Pezenti et al., 2021) with eight full-length sequences that were subsequently analyzed. A phylogenetic tree containing APNs from seven other lepidopteran species, belonging to the eight described classes of APNs in insects (Crava et al., 2010; Hughes, 2014; Lin et al., 2014), showed that the eight AgAPNs were distributed into different classes and named accordingly (AgAPN2, AgAPN3, AgAPN4, AgAPN5, AgAPN6, AgAPN8, AgAPN10 and AgAPN11), following the genome-wide unified nomenclature and classification of APN genes in lepidopteran insects (Guo et al., 2020). In silico analysis of the sequences revealed the presence of characteristic motifs “GAMENWG” and “HEXXHX18E” in most of the sequences; AgAPN10, AgAPN5 and AgAPN6 do not present the entire “GAMENWG” motif and AgAPN10 also does not present the entire “HEXXHX18E” motif. C-terminal GPI-anchor signal was identified for AgAPN2, AgAPN3, AgAPN4, AgAPN6 and AgAPN8 and N-terminal signal peptide cleavage site was identified in all sequences except for AgAPN11. It was interesting to note that AgAPN11 lacked both C-terminal GPI-anchor signal and N-terminal signal peptide cleavage site; a similar occurrence has been described for AjAPN9 in Achaea janata (Chauhan et al., 2021).

APN was initially shown as a Cry1Ac binding protein through TOA ligand blot analysis (Knight et al., 1994; Sangadala et al., 1994). Our present ligand-binding study revealed interaction of A. gemmatalis midgut epithelial cell membrane proteins with the activated Cry1Ac toxin obtained from B. thuringiensis sor. kurstaki, with a prominent interaction primarily seen with a ∼100 kDa protein, as well as a minor interaction with a ∼60 kDa protein. Corresponding bands from an SDS-PAGE were then excised and the identities of the proteins were confirmed by mass spectrometry analysis, in which several AgAPN sequences were identified; the same experiment was conducted for M. sexta, in which mass spectrometry of a band that positively bound to Cry2Ab toxin identified the protein as an APN (Onofre et al., 2017), and for Athethis lepigone, another lepidopteran species (Wang et al., 2017). Of the seven AgAPNs identified, four figured among the 100 highest PSMs (peptide spectrum match), which corresponds to a score of matches between experimental MS/MS spectra to the theoretical spectra predicted for the tryptic peptides, in which proteins with peptides best matched with the experimental spectrum are considered the most likely candidates (Wu et al., 2006). Therefore, supporting the presence of those proteins in the samples analyzed. An alkaline phosphatase (ALP), present in the transcriptome data, was also identified in the mass spectrometry analysis; this protein is frequently grouped with APNs, mainly because it is also anchored to the membrane by an GPI-anchor, besides displaying a similar role during the toxin’s mode of action (Ziga-Navarrete et al., 2021; Liu et al., 2021). As mentioned before, our group published a work characterizing an ALP from A. gemmatalis and describing an in vitro binding to Cry1Ac (da Silva et al., 2019); our data corroborate this finding.

In resistant strains, expression of genes encoding Cry receptors is commonly downregulated, whereas non-receptor paralogs are upregulated, believed to compensate for the absence of the down-regulated ones (Guo et al., 2022). The cabbage looper, Trichloplusia ni, presented down-regulation of TnAPN1 in Cry-resistant strains (Tiewsiri and Wang, 2011); tolerant strains of the castor semilooper larvae (A. janata) presented reduced expression of AjAPN2 while AjAPN4 was up-regulated (Chauhan et al., 2021); resistant strains of Spodoptera exigua reared in-lab lacked expression of a SeAPN1, whereas susceptible insects positively expressed this gene (Herrero et al., 2005). For Chilo suppressalis, knockdown of CsAPN6 and CsAPN8 reduced the larvae’s sensitivity to transgenic rice expressing Cry1 toxins (Sun et al., 2020), and previous studies had demonstrated the involvement of APNs in the toxin’s mode of action (Qiu et al., 2017; Wang et al., 2017). In our analysis we could identify an upregulation of AgAPN2 and AgAPN3 genes at 24 h of exposure to Bt spores, followed by a decrease in the following time points, as well as AgAPN5 and AgAPN6. Expression levels of AgAPN8 also demonstrated a downregulation, with great significance values. All these results could indicate a possible role of some AgAPNs in the mode of action of the Cry1Ac toxin in A. gemmatalis, but further studies are required for their characterization as putative receptors in the midgut, such as silencing through RNA interference and CRISPR/Cas9.

In conclusion, the identification of some AgAPNs as Cry-binding proteins in the ligand blot, as well as the expression analysis of all isoforms, could give us a direction in the sense of which of those proteins could be in fact involved in this process in A. gemmatalis midgut. Overall, this study gives a perspective of Cry1Ac possible receptors in A. gemmatalis, enhancing our understanding of this mode of action in this great agricultural pest.

Data availability statement

The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE (Perez-Riverol et al., 2022) partner repository with the dataset identifier PXD055041.

Ethics statement

Ethical approval was not required for the study involving animals in accordance with the local legislation and institutional requirements because this manuscript presents research on animals that do not require ethical approval for their study.

Author contributions

ML: Conceptualization, Investigation, Methodology, Supervision, Validation, Visualization, Writing–original draft, Writing–review and editing. IP: Investigation, Visualization, Writing–review and editing. LR: Investigation, Methodology, Writing–review and editing. AMd: Investigation, Writing–review and editing. AMe: Investigation, Writing–review and editing. YS: Investigation, Writing–review and editing. MM: Investigation, Writing–review and editing. MJ: Resources, Writing–review and editing. FN: Resources, Writing–review and editing. CA: Funding acquisition, Resources, Writing–review and editing. GD: Resources, Writing–review and editing. FG: Funding acquisition, Methodology, Resources, Visualization, Writing–review and editing. DO: Conceptualization, Funding acquisition, Methodology, Project administration, Resources, Supervision, Validation, Visualization, Writing–original draft, Writing–review and editing.

Funding

The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This work was supported by Carlos Chagas Filho Foundation for Supporting Research in the State of Rio de Janeiro (FAPERJ grants numbers E-26/210.173/2023 and E-26/210.996/2019); Coordination for the Improvement of Higher Education Personnel (CAPES - Finance Code 001).

Acknowledgments

We are grateful to Professor Rafael Dias Mesquita (Federal University of Rio de Janeiro) for the aid with the phylogenetic analysis, providing servers to run our experiment. We are also grateful to Professors Ana Claudia do Amaral Melo, Anderson de Sá Pinheiro, Bianca Cruz Neves, Marcia Regina Soares da Silva, Mônica Ferreira Moreira Carvalho Cardoso and Rodrigo Volcan Almeida (Federal University of Rio de Janeiro) for kindly providing laboratory supplies and facilities. To João Henrique de Oliveira Rangel, Igor Oliveira de Almeida and Danielle Aline da Silva Ribeiro for the help in maintaining the insect colony.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2024.1484489/full#supplementary-material

References

Adang M. J. (2004). Insect aminopeptidase N. Handb. Proteolytic Enzym, 296–299. doi:10.1016/b978-0-12-079611-3.50079-3

CrossRef Full Text | Google Scholar

Adang M. J., Crickmore N., Jurat-Fuentes J. L. (2014). Diversity of Bacillus thuringiensis crystal toxins and mechanism of action. Adv. Insect Physiology, 39–87. doi:10.1016/b978-0-12-800197-4.00002-6

CrossRef Full Text | Google Scholar

Aimanova K. G., Zhuang M., Gill S. S. (2006). Expression of cry1ac cadherin receptors in insect midgut and cell lines. J. Invertebr. Pathology 92 (3), 178–187. doi:10.1016/j.jip.2006.04.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Bel Y., Sheets J. J., Tan S. Y., Narva K. E., Escriche B. (2017). Toxicity and binding studies of Bacillus thuringiensis cry1ac, cry1f, cry1c, and cry2a proteins in the soybean pests Anticarsia gemmatalis and Chrysodeixis (Pseudoplusia) includens. Appl. Environ. Microbiol. 83 (11), 003266–e417. doi:10.1128/aem.00326-17

PubMed Abstract | CrossRef Full Text | Google Scholar

Bel Y., Zack M., Narva K., Escriche B. (2019). Specific binding of Bacillus thuringiensis cry1ea toxin, and cry1ac and Cry1Fa competition analyses in Anticarsia gemmatalis and Chrysodeixis includens. Sci. Rep. 9 (1), 18201. doi:10.1038/s41598-019-54850-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Bernardi O., Malvestiti G. S., Dourado P. M., Oliveira W. S., Martinelli S., Berger G. U., et al. (2012). Assessment of the high-dose concept and level of control provided by Mon 87701 × Mon 89788 soybean against Anticarsia gemmatalis and Pseudoplusia includens (Lepidoptera: noctuidae) in Brazil. Pest Manag. Sci. 68 (7), 1083–1091. doi:10.1002/ps.3271

PubMed Abstract | CrossRef Full Text | Google Scholar

Bravo A., Gill S. S., Soberón M. (2007). Mode of action of Bacillus thuringiensis cry and cyt toxins and their potential for insect control. Toxicon 49 (4), 423–435. doi:10.1016/j.toxicon.2006.11.022

PubMed Abstract | CrossRef Full Text | Google Scholar

Bravo A., Gómez I., Conde J., Muñoz-Garay C., Sánchez J., Miranda R., et al. (2004). Oligomerization triggers binding of a Bacillus thuringiensis Cry1Ab pore-forming toxin to aminopeptidase N receptor leading to insertion into membrane microdomains. Biochimica Biophysica Acta (BBA) - Biomembr. 1667 (1), 38–46. doi:10.1016/j.bbamem.2004.08.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Bravo A., Soberón M., Gill S. S. (2005). Bacillus thuringiensis: mechanisms and use. Compr. Mol. Insect Sci., 175–205. doi:10.1016/b0-44-451924-6/00081-8

CrossRef Full Text | Google Scholar

Cantalapiedra C. P., Hernández-Plaza A., Letunic I., Bork P., Huerta-Cepas J. (2021). EggNOG-mapper v2: functional annotation, orthology assignments, and domain prediction at the metagenomic scale. Mol. Biol. Evol. 38 (12), 5825–5829. doi:10.1093/molbev/msab293

PubMed Abstract | CrossRef Full Text | Google Scholar

Chattopadhyay P., Banerjee G. (2018). Recent advancement on chemical arsenal of BT toxin and its application in pest management system in Agricultural Field. 3 Biotech. 8 (4), 201. doi:10.1007/s13205-018-1223-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Chauhan V. K., Dhania N. K., Lokya V., Bhuvanachandra B., Padmasree K., Dutta-Gupta A. (2021). Midgut aminopeptidase N expression profile in castor semilooper (Achaea janata) during sublethal cry toxin exposure. J. Biosci. 46 (1), 29. doi:10.1007/s12038-021-00148-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen G., Wang Y., Liu Y., Chen F., Han L. (2020). Differences in midgut transcriptomes between resistant and susceptible strains of Chilo suppressalis to cry1c toxin. BMC Genomics 21 (1), 634. doi:10.1186/s12864-020-07051-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen J., Brown M. R., Hua G., Adang M. J. (2005). Comparison of the localization of Bacillus thuringiensis Cry1A δ-endotoxins and their binding proteins in larval midgut of tobacco hornworm, Manduca sexta. Cell Tissue Res. 321 (1), 123–129. doi:10.1007/s00441-005-1124-6

PubMed Abstract | CrossRef Full Text | Google Scholar

CONAB (2023). Acompanhamento da Safra Brasileira de Grãos.

Google Scholar

Crava C. M., Bel Y., Lee S. F., Manachini B., Heckel D. G., Escriche B. (2010). Study of the aminopeptidase N gene family in the Lepidopterans Ostrinia nubilalis (Hübner) and Bombyx mori (L.): sequences, mapping and expression. Insect Biochem. Mol. Biol. 40 (7), 506–515. doi:10.1016/j.ibmb.2010.04.010

PubMed Abstract | CrossRef Full Text | Google Scholar

da Silva G., Costa Ramos L. F., Dos Santos Seckler H., Mendonça Gomes F., Reis Cortines J., Ramos I., et al. (2019). Biochemical characterization of digestive membrane-associated alkaline phosphatase from the velvet bean caterpillar Anticarsia gemmatalis’, archives of insect biochemistry and physiology. Arch. Insect Biochem. Physiol. 102 (1), e21591. doi:10.1002/arch.21591

PubMed Abstract | CrossRef Full Text | Google Scholar

dos Reis L. C., Santos e Silva C. M., Bezerra B. G., Spyrides M. H. C. (2020). Caracterização da variabilidade da precipitação no MATOPIBA, região produtora de soja. Rev. Bras. Geogr. Física 13 (4), 1425–1441. doi:10.26848/rbgf.v13.4.p1425-1441

CrossRef Full Text | Google Scholar

Endo H. (2022). Molecular and kinetic models for pore formation of Bacillus thuringiensis cry toxin. Toxins 14 (7), 433. doi:10.3390/toxins14070433

PubMed Abstract | CrossRef Full Text | Google Scholar

Erlanger B. F., Kokowsky N., Cohen W. (1961). The preparation and properties of two new chromogenic substrates of Trypsin. Archives Biochem. Biophysics 95 (2), 271–278. doi:10.1016/0003-9861(61)90145-x

PubMed Abstract | CrossRef Full Text | Google Scholar

FAO (2024). Faostat. Available at: https://www.fao.org/faostat/ (Accessed August 09, 2024).

Google Scholar

Fernandes F. O., de Souza T. D., Sanches A. C., Dias N. P., Desiderio J. A., Polanczyk R. A. (2021). Sub-lethal effects of a Bt-based bioinsecticide on the biological conditioning of Anticarsia gemmatalis. Ecotoxicol. Lond. Engl. 30 (10), 2071–2082. doi:10.1007/s10646-021-02476-5

CrossRef Full Text | Google Scholar

Fiuza L. M., Knaak N., da Silva R. F. P., Henriques J. A. P. (2013). Receptors and lethal effect of Bacillus thuringiensis insecticidal crystal proteins to the Anticarsia gemmatalis (lepidoptera, Noctuidae). ISRN Microbiol. 2013, 940284–940287. doi:10.1155/2013/940284

PubMed Abstract | CrossRef Full Text | Google Scholar

Flores-Escobar B., Rodríguez-Magadan H., Bravo A., Soberón M., Gómez I. (2013). Differential role of Manduca sexta aminopeptidase-N and alkaline phosphatase in the mode of action of cry1aa, Cry1Ab, and cry1ac toxins from Bacillus thuringiensis. Appl. Environ. Microbiol. 79 (15), 4543–4550. doi:10.1128/aem.01062-13

PubMed Abstract | CrossRef Full Text | Google Scholar

Gill M., Ellar D. (2002). Transgenic drosophila reveals a functional in vivo receptor for the Bacillus thuringiensis toxin cry1ac1. Insect Mol. Biol. 11 (6), 619–625. doi:10.1046/j.1365-2583.2002.00373.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Gomes F. M., Carvalho D. B., Machado E. A., Miranda K. (2013). Ultrastructural and functional analysis of secretory goblet cells in the midgut of the Lepidopteran Anticarsia gemmatalis. Cell Tissue Res. 352 (2), 313–326. doi:10.1007/s00441-013-1563-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Gómez I., Arenas I., Benitez I., Miranda-Ríos J., Becerril B., Grande R., et al. (2006). Specific epitopes of domains II and III of Bacillus thuringiensis Cry1Ab toxin involved in the sequential interaction with cadherin and aminopeptidase-N receptors in Manduca sexta. J. Biol. Chem. 281 (45), 34032–34039. doi:10.1074/jbc.m604721200

PubMed Abstract | CrossRef Full Text | Google Scholar

Gómez I., Sánchez J., Miranda R., Bravo A., Soberón M. (2002). Cadherin-like receptor binding facilitates proteolytic cleavage of helix α-1 in domain I and oligomer pre-pore formation of Bacillus thuringiensis cry1ab toxin. FEBS Lett. 513 (2–3), 242–246. doi:10.1016/s0014-5793(02)02321-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Guo Z., Guo L., Qin J., Ye F., Sun D., Wu Q., et al. (2022). A single transcription factor facilitates an insect host combating Bacillus thuringiensis infection while maintaining fitness. Nat. Commun. 13 (1), 6024. doi:10.1038/s41467-022-33706-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Guo Z., Kang S., Sun D., Gong L., Zhou J., Qin J., et al. (2020). MAPK-dependent hormonal signaling plasticity contributes to overcoming Bacillus thuringiensis toxin action in an insect host. Nat. Commun. 11 (1), 3003. doi:10.1038/s41467-020-16608-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Hafkenscheid J. C. M. (1984). “Aminopeptidases and amino acid arylamidases in methods of enzymatic analysis,” in Enzymes 3; peptidases, proteinases and their inhibitors. Editor J. Bergemayer (Weinheim: Verlag Chemie), 2–34.

Google Scholar

Herrero S., Bel Y., Hernández-Martínez P., Ferré J. (2016). Susceptibility, mechanisms of response and resistance to Bacillus thuringiensis toxins in Spodoptera spp. Curr. Opin. Insect Sci. 15, 89–96. doi:10.1016/j.cois.2016.04.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Herrero S., Gechev T., Bakker P. L., Moar W. J., de Maagd R. A. (2005). Bacillus thuringiensis Cry1ca-resistant Spodoptera exigua lacks expression of one of four aminopeptidase N genes. BMC Genomics 6 (1), 96. doi:10.1186/1471-2164-6-96

PubMed Abstract | CrossRef Full Text | Google Scholar

Hoffman-Campo C. B., Moscardi F., Correa-Ferreira B. S., Oliveira L. J., Sosa-Gomez D. R., Panizzi A. R., et al. (2000). Pragas da soja no Brasil e seu manejo integrado. Londrina: EMBRAPA-CNPSo, 70.

Google Scholar

Hoffmann-Campo C. B., de Oliveira E. B., Moscardi F. (1985). Criação massal da lagarta da soja (Anticarsia gemmatalis). Available at: https://www.infoteca.cnptia.embrapa.br/infoteca/handle/doc/445420.

Google Scholar

Horikoshi R. J., Bernardi O., Godoy D. N., Semeão A. A., Willse A., Corazza G. O., et al. (2021a). Resistance status of Lepidopteran soybean pests following large-scale use of Mon 87701 × Mon 89788 soybean in Brazil. Sci. Rep. 11 (1), 21323. doi:10.1038/s41598-021-00770-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Horikoshi R. J., Dourado P. M., Berger G. U., de S Fernandes D., Omoto C., Willse A., et al. (2021b). Large-scale assessment of Lepidopteran soybean pests and efficacy of cry1ac soybean in Brazil. Sci. Rep. 11 (1), 15956. doi:10.1038/s41598-021-95483-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Hughes A. L. (2014). Evolutionary diversification of aminopeptidase N in Lepidoptera by conserved clade-specific amino acid residues. Mol. Phylogenetics Evol. 76, 127–133. doi:10.1016/j.ympev.2014.03.014

PubMed Abstract | CrossRef Full Text | Google Scholar

IBS (2022). 2.0: an upgraded illustrator for the visualization of biological sequences. Nucleic Acids Res. doi:10.1093/nar/gkac373

CrossRef Full Text | Google Scholar

Ingle S. S., Trivedi N., Prasad R., Kuruvilla J., Rao K. K., Chhatpar H. S. (2001). Aminopeptidase-N from the Helicoverpa armigera (Hubner) brush border membrane vesicles as a receptor of Bacillus thuringiensis cry1ac δ-endotoxin. Curr. Microbiol. 43 (4), 255–259. doi:10.1007/s002840010297

PubMed Abstract | CrossRef Full Text | Google Scholar

Knaak N., Fiuza L. M. (2005). Histopathology of Anticarsia gemmatalis Hübner (Lepidoptera; Noctuidae) treated with Nucleopolyhedrovirus and Bacillus thuringiensis serovar kurstaki. Braz. J. Microbiol. 36 (2), 196–200. doi:10.1590/s1517-83822005000200017

CrossRef Full Text | Google Scholar

Knight P. J., Crickmore N., Ellar D. J. (1994). The receptor for Bacillus thuringiensis cryla(c) delta-endotoxin in the brush border membrane of the Lepidopteran Manduca sexta is aminopeptidase N. Mol. Microbiol. 11 (3), 429–436. doi:10.1111/j.1365-2958.1994.tb00324.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Laemmli U. K. (1970). Cleavage of structural proteins during the assembly of the head of bacteriophage T4. Nature 227 (5259), 680–685. doi:10.1038/227680a0

PubMed Abstract | CrossRef Full Text | Google Scholar

Lei Y., Zhu X., Xie W., Wu Q., Wang S., Guo Z., et al. (2014). Midgut transcriptome response to a Cry toxin in the diamondback moth, Plutella xylostella (Lepidoptera: plutellidae). Gene 533 (1), 180–187. doi:10.1016/j.gene.2013.09.091

PubMed Abstract | CrossRef Full Text | Google Scholar

Lin P., Cheng T., Jin S., Jiang L., Wang C., Xia Q. (2014). Structural, evolutionary and functional analysis of APN genes in the lepidoptera Bombyx mori. Gene 535 (2), 303–311. doi:10.1016/j.gene.2013.11.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu L., Li Z., Luo X., Zhang X., Chou S. H., Wang J., et al. (2021). Which is stronger? A continuing battle between cry toxins and insects. Front. Microbiol. 12, 665101. doi:10.3389/fmicb.2021.665101

PubMed Abstract | CrossRef Full Text | Google Scholar

Melo A. L., Soccol V. T., Soccol C. R. (2014). Bacillus thuringiensis: mechanism of action, resistance, and new applications: a Review. Crit. Rev. Biotechnol. 36 (2), 317–326. doi:10.3109/07388551.2014.960793

PubMed Abstract | CrossRef Full Text | Google Scholar

Onofre J., Gaytán M. O., Peña-Cardeña A., García-Gomez B. I., Pacheco S., Gómez I., et al. (2017). Identification of aminopeptidase-N2 as a Cry2Ab binding protein in Manduca sexta. Peptides 98, 93–98. doi:10.1016/j.peptides.2017.01.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Pacheco S., Gómez I., Gill S. S., Bravo A., Soberón M. (2009). Enhancement of insecticidal activity of Bacillus thuringiensis cry1a toxins by fragments of a toxin-binding cadherin correlates with oligomer formation. Peptides 30 (3), 583–588. doi:10.1016/j.peptides.2008.08.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Palma L., Muñoz D., Berry C., Murillo J., Caballero P. (2014). Bacillus thuringiensis toxins: an overview of their biocidal activity. Toxins 6 (12), 3296–3325. doi:10.3390/toxins6123296

PubMed Abstract | CrossRef Full Text | Google Scholar

Pardo-López L., Soberón M., Bravo A. (2013). Bacillus thuringiensis insecticidal three-domain cry toxins: mode of action, insect resistance and consequences for crop protection. FEMS Microbiol. Rev. 37 (1), 3–22. doi:10.1111/j.1574-6976.2012.00341.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Perez-Riverol Y., Bai J., Bandla C., García-Seisdedos D., Hewapathirana S., Kamatchinathan S., et al. (2022). The PRIDE database resources in 2022: a Hub for mass spectrometry-based proteomics evidences. Nucleic Acids Res. 50 (D1), D543–D552. doi:10.1093/nar/gkab1038

PubMed Abstract | CrossRef Full Text | Google Scholar

Pezenti L. F., Sosa-Gómez D. R., de Souza R. F., Vilas-Boas L. A., Gonçalves K. B., da Silva C. R. M., et al. (2021). Transcriptional profiling analysis of susceptible and resistant strains of Anticarsia gemmatalis and their response to Bacillus thuringiensis. Genomics 113 (4), 2264–2275. doi:10.1016/j.ygeno.2021.05.012

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Pinos D., Andrés-Garrido A., Ferré J., Hernández-Martínez P. (2021). Response mechanisms of invertebrates to Bacillus thuringiensis and its pesticidal proteins. Microbiol. Mol. Biol. Rev. 85 (1), e00007–e00020. doi:10.1128/mmbr.00007-20

PubMed Abstract | CrossRef Full Text | Google Scholar

Pozebon H., Marques R. P., Padilha G., O Neal M., Valmorbida I., Bevilaqua J. G., et al. (2020). Arthropod invasions versus soybean production in Brazil: a Review. J. Econ. Entomology 113 (4), 1591–1608. doi:10.1093/jee/toaa108

PubMed Abstract | CrossRef Full Text | Google Scholar

Qiu L., Cui S., Liu L., Zhang B., Ma W., Wang X., et al. (2017). Aminopeptidase N1 is involved in Bacillus thuringiensis cry1ac toxicity in the beet armyworm, Spodoptera exigua. Sci. Rep. 7 (1), 45007. doi:10.1038/srep45007

PubMed Abstract | CrossRef Full Text | Google Scholar

Rajagopal R., Agrawal N., Selvapandiyan A., Sivakumar S., Ahmad S., Bhatnagar R. K. (2003). Recombinantly expressed isoenzymic aminopeptidases from Helicoverpa armigera (American cotton bollworm) midgut display differential interaction with closely related Bacillus thuringiensis insecticidal proteins. Biochem. J. 370 (3), 971–978. doi:10.1042/bj20021741

PubMed Abstract | CrossRef Full Text | Google Scholar

Sanchis V. (2011). From microbial sprays to insect-resistant transgenic plants: history of the biospesticide Bacillus thuringiensis. A Review. Agron. Sustain. Dev. 31 (1), 217–231. doi:10.1051/agro/2010027

CrossRef Full Text | Google Scholar

Sangadala S., Walters F. S., English L. H., Adang M. J. (1994). A mixture of Manduca sexta aminopeptidase and phosphatase enhances Bacillus thuringiensis insecticidal cryia(c) toxin binding and 86RB(+)-K+ efflux in vitro. J. Biol. Chem. 269 (13), 10088–10092. doi:10.1016/s0021-9258(17)36993-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Sato R. (2002). Aminopeptidase N as a receptor for Bacillus thuringiensis cry toxins. Adv. Microb. Control Insect Pests, 1–13. doi:10.1007/978-1-4757-4437-8_1

CrossRef Full Text | Google Scholar

Soberón M., Monnerat R., Bravo A. (2018). Mode of action of cry toxins from Bacillus thuringiensis and resistance mechanisms. Toxinology, 15–27. doi:10.1007/978-94-007-6449-1_28

CrossRef Full Text | Google Scholar

Stamatakis A. (2014). RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinforma. Oxf. Engl. 30 (9), 1312–1313. doi:10.1093/bioinformatics/btu033

PubMed Abstract | CrossRef Full Text | Google Scholar

Sun Y., Yang P., Jin H., Liu H., Zhou H., Qiu L., et al. (2020). Knockdown of the aminopeptidase N genes decreases susceptibility of Chilo suppressalis larvae to cry1ab/cry1ac and cry1ca. Pesticide Biochem. Physiology 162, 36–42. doi:10.1016/j.pestbp.2019.08.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Tabashnik B. E., Fabrick J. A., Carrière Y. (2023). Global patterns of insect resistance to transgenic BT crops: the first 25 years. J. Econ. Entomology 116 (2), 297–309. doi:10.1093/jee/toac183

CrossRef Full Text | Google Scholar

Tamura K., Stecher G., Kumar S. (2021). MEGA11: molecular evolutionary genetics analysis version 11. Mol. Biol. Evol. 38 (7), 3022–3027. doi:10.1093/molbev/msab120

PubMed Abstract | CrossRef Full Text | Google Scholar

Tiewsiri K., Wang P. (2011). Differential alteration of two aminopeptidases n associated with resistance to Bacillus thuringiensis toxin cry1ac in cabbage looper. Proc. Natl. Acad. Sci. 108 (34), 14037–14042. doi:10.1073/pnas.1102555108

PubMed Abstract | CrossRef Full Text | Google Scholar

Towbin H., Staehelin T., Gordon J. (1979). Electrophoretic transfer of proteins from polyacrylamide gels to nitrocellulose sheets: procedure and some applications. Proc. Natl. Acad. Sci. Sep 76 (9), 4350–4354. doi:10.1073/pnas.76.9.4350

PubMed Abstract | CrossRef Full Text | Google Scholar

Valaitis A. P. (2011). Localization of Bacillus thuringiensis cry1a toxin-binding molecules in gypsy moth larval gut sections using fluorescence microscopy. J. Invertebr. Pathology 108, 69–75. doi:10.1016/j.jip.2011.07.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Valaitis A. P., Augustin S., Clancy K. M. (1999). Purification and characterization of the western spruce budworm larval midgut proteinases and comparison of gut activities of laboratory-reared and field-collected insects. Insect Biochem. Mol. Biol. 29 (5), 405–415. doi:10.1016/s0965-1748(99)00017-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Valdes C., Gillespie J., Dohlman E. (2023). Soybean production, marketing costs, and export competitiveness in Brazil and the United States. Economic Information Bulletin: United States Department of Agriculture, Economic Research Service

Google Scholar

Visweshwar R., Sharma H. C., Akbar S. M. D., Sreeramulu K. (2015). Elimination of gut microbes with antibiotics confers resistance to Bacillus thuringiensis toxin proteins in Helicoverpa armigera (Hubner). Appl. Biochem. Biotechnol. 177 (8), 1621–1637. doi:10.1007/s12010-015-1841-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang C., Deng Z., Yuan J., Xu K., Sha L., Guan X., et al. (2023). Removal of an aminopeptidase N from midgut brush border does not affect susceptibility of Spodoptera litura (Lepidoptera: Noctuidae) larvae to four insecticidal proteins of Bacillus thuringiensis (Bacillales: Bacillaceae). J. Econ. Entomol. 116 (1), 223–232. doi:10.1093/jee/toac184

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang L., Gu S. H., Nangong Z. Y., Song P., Wang Q. Y. (2017). Aminopeptidase N5 (APN5) as a putative functional receptor of cry1ac toxin in the larvae of Athetis lepigone. Curr. Microbiol. 74 (4), 455–459. doi:10.1007/s00284-017-1215-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang P., Zhang X., Zhang J. (2005). Molecular characterization of four midgut aminopeptidase N isozymes from the Cabbage Looper, Trichoplusia ni. Insect Biochem. Mol. Biol. 35 (6), 611–620. doi:10.1016/j.ibmb.2005.02.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang X. Y., Du L. X., Liu C. X., Gong L., Han L. Z., Peng Y. F. (2017). RNAi in the Striped Stem Borer, Chilo suppressalis, establishes a functional role for aminopeptidase N in cry1ab intoxication. J. Invertebr. Pathology 143, 1–10. doi:10.1016/j.jip.2016.11.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Wolfersberger M., Luethy P., Maurer A., Parenti P., Sacchi F., Giordana B., et al. (1987). Preparation and partial characterization of amino acid transporting brush border membrane vesicles from the larval midgut of the cabbage butterfly (Pieris brassicae). A, Comp. Physiol. 86 (2), 301–308. doi:10.1016/0300-9629(87)90334-3

CrossRef Full Text | Google Scholar

Wu F.-X., Gagné P., Droit A., Poirier G. G. (2006). RT-PSM, a real-time program for peptide-spectrum matching with statistical significance. Rapid Commun. mass Spectrom. RCM 20 (8), 1199–1208. doi:10.1002/rcm.2435

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu C., Wang B. C., Yu Z., Sun M. (2014). Structural insights into Bacillus thuringiensis Cry, Cyt and parasporin toxins. Toxins 6 (9), 2732–2770. doi:10.3390/toxins6092732

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang Y., Zhu Y. C., Ottea J., Husseneder C., Leonard B. R., Abel C., et al. (2010). Molecular characterization and RNA interference of three midgut aminopeptidase N isozymes from Bacillus thuringiensis-susceptible and -resistant strains of sugarcane borer, Diatraea saccharalis. Insect Biochem. Mol. Biol. 40 (8), 592–603. doi:10.1016/j.ibmb.2010.05.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Ziga-Navarrete F., Bravo A., Sobern M., Gmez I. (2021). Role of GPI-anchored membrane receptors in the mode of action of Bacillus thuringiensis cry toxins. Integr. Pest Manag. Pest Control - Curr. Future Tactics. doi:10.5772/30718

CrossRef Full Text | Google Scholar

Keywords: aminopeptidase N, Cry toxin, receptor, Bacillus thuringiensis, Anticarsia gemmatalis

Citation: Lanzaro MD, Padilha I, Ramos LFC, Mendez APG, Menezes A, Silva YM, Martins MR, Junqueira M, Nogueira FCS, AnoBom CD, Dias GM, Gomes FM and Oliveira DMP (2024) Cry1Ac toxin binding in the velvetbean caterpillar Anticarsia gemmatalis: study of midgut aminopeptidases N. Front. Physiol. 15:1484489. doi: 10.3389/fphys.2024.1484489

Received: 21 August 2024; Accepted: 09 October 2024;
Published: 29 October 2024.

Edited by:

Guenter Artur Schaub, Ruhr University Bochum, Germany

Reviewed by:

Zhaojiang Guo, Chinese Academy of Agricultural Sciences, China
Carlos Peres Silva, Federal University of Santa Catarina, Brazil

Copyright © 2024 Lanzaro, Padilha, Ramos, Mendez, Menezes, Silva, Martins, Junqueira, Nogueira, AnoBom, Dias, Gomes and Oliveira. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: D. M. P. Oliveira, ZGFuaW9saXZlaXJhQGlxLnVmcmouYnI=

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.