Next Article in Journal
Sterile Insect Technique: Successful Suppression of an Aedes aegypti Field Population in Cuba
Next Article in Special Issue
Resistance in the Genus Spodoptera: Key Insect Detoxification Genes
Previous Article in Journal
Utilizing Red Spotted Apollo Butterfly Transcriptome to Identify Antimicrobial Peptide Candidates against Porphyromonas gingivalis
Previous Article in Special Issue
The Essential and Enigmatic Role of ABC Transporters in Bt Resistance of Noctuids and Other Insect Pests of Agriculture
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Geographic Monitoring of Insecticide Resistance Mutations in Native and Invasive Populations of the Fall Armyworm

1
DGIMI, University of Montpellier, INRAE, F-34095 Montpellier, France
2
CIRAD, UPR AIDA, F-34398 Montpellier, France
3
PHIM Plant Health Institute, University of Montpellier, IRD, CIRAD, INRAE, Institut Agro, F-34398 Montpellier, France
4
AIDA, University of Montpellier, CIRAD, F-34398 Montpellier, France
5
Black Mountain Laboratories, CSIRO, Canberra, ACT 2601, Australia
*
Authors to whom correspondence should be addressed.
Insects 2021, 12(5), 468; https://doi.org/10.3390/insects12050468
Submission received: 31 March 2021 / Revised: 11 May 2021 / Accepted: 12 May 2021 / Published: 18 May 2021

Abstract

:

Simple Summary

The moth fall armyworm (Spodoptera frugiperda) is a major agricultural pest insect damaging a wide range of crops, especially corn. Field evolved resistance against Bacillus thuringiensis (Bt) toxins and synthetic insecticides has been repeatedly reported. While the fall armyworm is native to the Americas, its biological invasion was first reported from West Africa in 2016. Since then, this pest has been detected across sub-Saharan and North Africa, Asia, and Oceania. Here, we examine the geographical distribution of mutations causing resistance against Bt or synthetic insecticides to test if the invasion was accompanied by the spread of resistance mutations using 177 individuals collected from 12 geographic populations including North and South America, West and East Africa, India, and China. We observed that Bt resistance mutations generated in Puerto Rico or Brazil were found only from their native populations, while invasive populations had higher copy numbers of cytochrome P450 genes and higher proportions of resistance mutations at AChE, which are known to cause resistance against synthetic insecticides. This result explains the susceptibility to Bt insecticides and the resistance against synthetic insecticides in invasive Chinese populations. This information will be helpful in investigating the cause and consequence associated with insecticide resistance.

Abstract

Field evolved resistance to insecticides is one of the main challenges in pest control. The fall armyworm (FAW) is a lepidopteran pest species causing severe crop losses, especially corn. While native to the Americas, the presence of FAW was confirmed in West Africa in 2016. Since then, the FAW has been detected in over 70 countries covering sub-Saharan Africa, the Middle East, North Africa, South Asia, Southeast Asia, and Oceania. In this study, we tested whether this invasion was accompanied by the spread of resistance mutations from native to invasive areas. We observed that mutations causing Bt resistance at ABCC2 genes were observed only in native populations where the mutations were initially reported. Invasive populations were found to have higher gene numbers of cytochrome P450 genes than native populations and a higher proportion of multiple resistance mutations at acetylcholinesterase genes, supporting strong selective pressure for resistance against synthetic insecticides. This result explains the susceptibility to Bt insecticides and resistance to various synthetic insecticides in Chinese populations. These results highlight the necessity of regular and standardized monitoring of insecticide resistance in invasive populations using both genomic approaches and bioassay experiments.

1. Introduction

Insecticide resistance is one of the main challenges for the control of insect pests. Commonly used insecticides to control pest insects can be classified into two main types. The first group is Bt (Bacillus thuringiensis) toxins, which are generally produced by transgenic crops. The global areas of planted Bt crops is positively correlated with the number of Bt resistance pest species, implying that field-evolved Bt resistance is prevalent [1]. Field-evolved Bt resistance is a particularly serious issue in corn and cotton because the majority of those plants are Bt crops in the USA and because only GM-cotton expressing Bt toxins are grown in Australia. Field-evolved Bt resistance has been observed in several major insect Lepidoptera pest species including maize stalk borer (Busseola fusca), western corn rootworm (Diabrotica virgifera virgifera), cotton bollworm (Helicoverpa zea), Old World cotton bollworm (Helicoverpa armigera), Helicoverpa punctigera, fall armyworm (Spodoptera frugiperda), and pink bollworm (Pectinophora gossypiella) [1]. The other group is synthetic insecticides such as organophosphates, pyrethroids, neonicotinoids, organochlorides, carbamates, ryanoids, and spinosyns. By 2021, more than 600 arthropod species have been included in the arthropod pesticide resistance database (https://www.pesticideresistance.org/) (accessed on 21 January 2021).
Several genetic mechanisms of Bt resistance have been reported (reviewed in [2]). If a mutation in a midgut receptor causes reduced physical binding with Bt proteins, the toxicity of Bt can be decreased. These receptors include cadherin, aminopeptidase-N, alkaline phosphatases, and ATP-binding cassette transporters. Alternatively, altered processing of Bt prototoxin, increased immune status, sequestration of Bt protein, and accelerated recovery of epithelial integrity can also increase Bt resistance. The resistance against synthetic insecticides can be caused by detoxification genes [3] such as cytochrome P450 gene, esterase, and glutathione S transferase. In addition, mutations in calcium channels, acetylcholinesterase, and nicotinic acetylcholine receptors may also cause resistance [4,5,6,7].
Fall armyworm (FAW, Spodoptera frugiperda, J.E.Smith) (Lepidoptera: Noctuidae: Noctuinae) is one of the most damaging pest insects, partly due to the extreme polyphagy by feeding on at least 353 species in 76 plant families, which include several economically important crops such as corn, rice, sorghum, sugarcane, cotton, and soybean [8] Genomic analyses have demonstrated a rapid expansion of detoxification genes [9,10,11], potentially associated with this extreme polyphagy by overcoming plant defense toxins from diverse plants.
Field-evolved resistance to Bt insecticide in the FAW has been reported from Puerto Rico [12,13], Brazil [14,15], Argentina [16], and the USA [17]. Five mutations causing Bt resistance have been reported. Interestingly, all these mutations were observed in the ATP Binding Cassette Subfamily C Member 2 (ABCC2) gene, which encodes ATP-binding cassette transporter proteins, implying parallel evolution. In a population from Puerto Rico, 2bp GC insertion at ABCC2 causes a frameshift mutation and a premature stop codon, which results in improper binding between ABCC2 and Bt Cry1F and Cry1A proteins [18,19]. Boaventura et al. identified three causal mutations of Cry1F resistance at ABCC2 including GY deletion, P799K/R substitution, and G1088D substitution from a Brazilian population [20]. Guan et al. also identified a 12bp insertion, which causes a frameshift mutation and a premature stop codon, from another Brazilian population [21]. As well as the Bt resistance, field-evolved resistance to synthetic insecticides has been widely reported. For example, a population in Puerto Rico was reported to have resistance against a wide range of synthetic insecticides including flubendiamide, chlorantraniliprole, methomyl, thiodicarb, permethrin, chlorpyriphos, zeta-cypermethrin, deltamethrin, triflumuron, and spinetoram [22]. We observed that P450 genes cause resistance against deltamethrin [23]. Boaventura et al. reported that I4734M and G4946E substitutions at the ryanodine receptor (RyR) cause resistance against diamides [24]. Carvalho et al. reported that resistance against organophosphate and carbamates is caused by A314S, G340A, and F402V mutations at acetylcholinesterase (AChE) and that pyrethroid resistance is caused by T929I, L932F, and L1014F mutations at the voltage-gated sodium channel (VGSC) gene [4].
FAW is native to North and South America, and its presence was first confirmed in West Africa in 2016 [25]. Since then, FAW has been detected across almost the entire Sub-Saharan Africa, in northern Africa/the Near East (Egypt), Middle East Asia, South Asia, South East Asia, East Asia, Australia, and more recently, the Canary Islands located off the coast between Western Africa and Morocco, and New Caledonia in the South Pacific (https://www.cabi.org/ISC/fallarmyworm) (accessed on 21 January 2021). The invasion caused a serious reduction in corn production, between 21–53% in Africa [26]. While native populations are composed of two sympatric strains, corn (sfC) and rice strains (sfR), named after their preferred host-plants [27,28,29], the invasive populations are observed as hybrids between corn and rice strains [30,31,32], but are predominantly observed in corn fields. Genomic analyses support the multiple introductions of FAW from multiple native populations including Florida and Brazil to multiple invasive areas including China and Africa [31,33]. However, the invasive populations exhibited much more homogeneous DNA sequences than native populations, implying that these invasive populations had experienced extensive admixture since their establishment across the Old World [31].
The invasion process opens up the possibility that resistance mutations have been spread as well in the invasive area. Recently, Boaventura et al. observed that resistance mutations at AChE were frequently observed from two invasive populations in Kenya and Indonesia and two populations in Brazil and Puerto Rico, while GY deletion of ABCC2 genes was observed only in the Brazilian population [34]. Boaventura et al. also reported the presence of the VGSC L1014F allele with low frequency in an Indonesian population that has, to-date, yet to be reported in other invasive populations [34]. Zhang et al. reported that the 2 bp GC insertion was not observed from Chinese populations [32], while Guan et al. reported that the 12 bp insertion in exon 15 of the ABCC2 gene was not detected in FAW populations from Uganda, Malawi, and from six Chinese provinces [21]. However, we still lack a comprehensive understanding of the geographic distribution of these resistance mutations because global profiling of the causal resistance mutations has not yet been performed. We performed whole-genome resequencing from 177 individuals from both invasive and native populations sampled from 12 geographic populations in North America, South America, Africa, India, and China to investigate adaptive evolution associated with the invasive success in our previous study [31]. In this study, we undertook reanalysis of the resequencing data of the FAW populations to test the presence of causal resistance mutations against insecticides in the FAWs.

2. Methods

2.1. Resequencing Data

Resequencing data generated in our previous study were reused [31,33]. More specifically, this dataset included 29 sfC and 49 sfR individuals from invasive populations (Benin, Malawi, Uganda, China, and India) and 70 sfC and 29 sfR individuals from native populations (Mississippi, Florida, Puerto Rico, Brazil, French Guinea, Guadeloupe, and Mexico) (Figure 1) (please see Table S1 for more information). This dataset includes four Bt-resistant and six susceptible Brazilian individuals. Paired-end whole genome resequencing was performed with approximately 20X coverage with 300 bp insert size, and 150 bp read length using Illumina (San Diego, CA, USA) HiSeq 2500, HiSeq 4000, and Novaseq 6000. Variant calling was performed using the GATK-4.0.11.0 package [35] for single nucleotide variant (SNV) and CNVCaller [36] for copy number variation (CNV). The numbers of variations were 27,117,672 and 22,916 for SNV and CNV, respectively.

2.2. Gene Annotation and Statistical Analysis

ABCC2, RyR, AChE, and VGSC were annotated by aligning protein sequences obtained from NCBI and reference genomes using exonerate 2.2.0 [38] with the protein2genome model. The accession numbers of these genes were QGS83596, XP_022819835, MK226188, and KC435025, for ABCC2, AChE, RyR, and VGSC, respectively. Genetic variations within these genes were identified from the resequencing data [31] using IGV tools [39]. P450 genes were identified from the gene annotation files generated in our previous study (OGS7.0) [31]. The identification of the P450 clan was performed by blasting against OGS2.2 [9], which includes manually annotated P450 genes, with 98% cutoff using ncbi-blast-2.10.0+ [40]. The Fisher’s exact test was performed using the R package to test a statistical significance of increased gene duplication of P450 genes. Weir and Cockerham’s FST [41] between Bt susceptible and resistant individuals was calculated at the ABCC2 gene using vcftools v0.1.13 [42]. Increased FST at ABCC2 was tested by calculating the proportion of FST at randomly chosen loci with the same length of ABCC2 from whole genome sequences with 100,000 replications.

3. Results

3.1. Bt Insecticide Resistance–ABCC2 Gene

We tested the presence of five reported mutations at the ABCC2 gene causing the resistance against Bt insecticide. These mutations include (i) 2 bp GC insertion leading to a premature stop codon due to the frameshift (originally identified from a population in Puerto Rico) [18,19]; (ii) GY deletion (Brazil) [20]; (iii) 12 bp insertion causing premature stop codon due to frameshift (Brazil) [21]; (iv) P799K/R substitution (Brazil) [20]; and (v) G1088D substitution (Brazil) [20] (Figure 2). We identified the 2 bp GC insertion from 11 individuals from Puerto Rico in the resequencing data (Table 1). As this mutation was originally reported from a population in Puerto Rico, this result supports that this mutation remained in this original population in Puerto Rico and did not spread to the Old World by invasion.
GY deletion was identified from two individuals in the Brazilian population, and P799K/R was observed from the same individuals. This result implies that this mutation also remained in the original population. Intriguingly, these mutations exist at the homozygous state in these two Bt susceptible individuals, suggesting a possibility that these mutations alone are not sufficient for Bt resistance.
We identified two individuals with 12 bp insertion only from Brazil, also supporting that this mutation remains in the original population and confirms the analysis of Guan et al. [21]. G1088D was not observed from the resequencing data.
As unidentified mutations causing Bt resistance might exist, we investigated sequences with complete genetic differentiation between resistant and susceptible Brazilian individuals. In total, 49 SNPs had complete genetic differentiation (i.e., FST = 1), and all these positions were found to be intronic. This result shows that other causal mutations are not likely to exist in the tested individuals.
We tested the possibility that natural selection on the ABCC2 gene is responsible for the field-evolved Bt resistance in the Brazilian population. In this case, ABCC2 is expected to have increased genetic differentiation between resistant and susceptible groups. FST calculated between these two groups was 0.2922 at ABCC2 gene. In total, 7.7% of 100,000 randomly chosen loci with the same length of the ABCC2 gene had FST higher than 0.2922 (equivalent to p-value = 0.077) (Figure 3), indicating marginally significant statistical support for increased genetic differentiation. This result implies that natural selection might have increased the Bt resistance of the Brazilian population.

3.2. Synthetic Insecticide Resistance

P450 genes increase the level of synthetic insecticide resistance in FAW, and P450 gene duplication has been suggested as a possible cause underpinning resistance [23]. In total, 99 P450 genes were identified from the reference genome (Table S2). Copy number variation compared with the reference genome was observed from 47 P450 genes, which include4, 12, and 22 genes belonging to clan2, clan3, and clan4, respectively. In all clades, invasive populations have higher numbers of duplicated P450 genes per individual than native populations for all clades (Figure 4). On the other hand, native populations have higher numbers of deleted P450 genes per individuals than invasive populations for all clades. Consequently, invasive populations have overrepresented numbers of duplicated P450 genes compared with deleted ones than the native populations (p-value = 6.018 × 10−21, 0.0004553, 0.04414 for clan2, clan3, and clan4, respectively; two-tailed Fisher’s exact test). This result supports the hypothesis that invasive populations have higher P450 gene numbers than native populations through CNVs.
The resequencing data show all three reported resistance mutations at the AChE (Figure 5), as shown by Boaventura et al. [34]. F209V mutations were observed from populations in Uganda, Malawi, Brazil, and China, as reported from Guan et al. [21] as well as India, Benin, Puerto Rico, Mexico, Mississippi, and Florida. Invasive populations have a higher proportion of individuals with resistance F209V mutations (80.26%) than native populations (72.37%) (p = 0.001097, two-tailed Fisher’s exact test). A201S mutation was found from populations in Uganda, Malawi, and Brazil, as shown by Guan et al., but also from Benin, Mississippi, and Florida. Invasive populations have a higher proportion of individuals with A201S resistance mutations (20.0%) than native populations (7.14%) (p = 0.01966). Interestingly, two codons were found for the susceptible ‘A’ in the population from Benin (GCG and GCA). In Benin, the GCG codon was observed from 28 individuals (as well as other populations). The GCA codon was observed from three individuals in Benin. The most parsimonious explanation is that the GCA codon was newly generated from GCG by a point mutation at the third codon position. The G227A mutation was found in Brazil, as per Guan et al. [21] as well as in Puerto Rico and Florida. The proportion of individuals with this resistance mutation was 17.65% in native populations. Resistance mutations at RyR (I4734M and G4946E substitutions) or the VGSC (T929I, L932F, and L1014F substitutions) were not observed in the analyzed populations of FAW.

4. Discussion

The global spread of FAW through invasion may be accompanied by the spread of resistance mutations to Bt or synthetic insecticides to invasive areas. This is known to occur in other species. The population of the cotton bollworm Helicoverpa armigera is thought to have invaded South America with a pyrethroid resistance allele, which was then shared with a sister species, H. zea [43,44]. Therefore, monitoring the geographical distribution of resistance mutations in the context of invasion can provide useful information for pest FAW control and insecticide resistance management. We identified known resistance mutations from the whole genome sequences of 177 individuals collected from 12 geographic populations in North America, South America, Africa, India, and China. We showed that resistance Bt mutations at ABCC2 remained at their origin in Puerto Rico and Brazil, meaning that these mutations do not spread to invasive areas. However, invasive populations have increased copy numbers of P450 genes compared with native populations through copy number variations. We also observed that invasive populations had increased proportions of resistance mutations at AChE. These results can be best explained by selective pressure causing increased resistance against synthetic insecticides in invasive populations.
Bioassay experiments demonstrated that Chinese populations are susceptible to Bt Cry1Fa toxins [31]. The analyzed Chinese samples in this study did not have mutations causing Cry1Fa resistance, explaining the susceptibility to this Bt toxin. Interestingly, a South African population showed moderate resistance against Bt Cry1A [32], while GC insertion causing Cry1A resistance [18,19] was not identified from invasive populations in the resequencing data. Therefore, the South African population could potentially have a resistance mutation that is yet to be characterized in native populations, and we cannot exclude the possibility of on-going gene flow of unidentified resistance mutations either from native to invasive populations or from other invasive populations that represent independent introduction from South African FAWs. More studies are necessary to identify causal resistance mutation in the South African populations and to identify the origin of this resistance mutation.
If a resistance Bt mutation exists in invasive populations with low frequency, individuals with the resistance mutations might not be included in the analyzed resequencing data. If this is true, existing resistance mutations in invasive populations might be undetected. The probability of not finding a resistance mutation is (1 − p)2N, where p is the allele frequency of resistance mutations and N is the number of diploid organisms. If the allele frequency of resistance Bt mutation is 0.01, representing a rare allele, the probability of missing this mutation in our resequencing data is only (1 − 0.01)2 × 99 = 13.67%. Therefore, we do not believe that a resistance Bt mutation is likely to exist. Even in the case that this possibility of 13.67% is true, the resistance mutation is likely to be eliminated in a population by genetic drift due to the low allele frequency (0.01).
A whole-genome analysis is necessary to contrast the abundance of P450 genes between invasive and native populations, because P450 genes exist in multi-copies scattered across the genome [9,10,11]. Our whole genome analysis at a population level shows an increased copy number of P450 genes in the invasive FAW populations. Bioassay experiments demonstrated that Chinese populations have resistance to pyrethroid insecticides [32,45]. As P450 genes are known to cause resistance against deltamethrin [23], a type of pyrethroid in FAWs, the increased P450 gene numbers might be responsible for this resistance. Chinese populations also have resistance against organophosphate [32,45] as well as oxadiazine and diamide [32]. As P450 proteins are probably able to detoxify other types of insecticides than pyrethroid, the increased P450 genes may be responsible for the resistance against a wider range of insecticides.
In addition, the increased proportion of resistance mutations at AChE can cause the increased resistance against synthetic insecticides in invasive populations. Multiple resistance mutations with the increased proportion imply strong selective pressure for resistance against synthetic insecticide. The vast majority of beneficial mutations are quickly removed in a population due to stochastic fluctuation of allele frequency by the genetic bottleneck, and the survival probability of beneficial mutations is just two times that of the selective coefficient. Therefore, the probability of observing a single resistance mutation at AChE is proportional to the extent of increased fitness due to the resistance. The probability of observing multiple resistance mutations is a function of the multiplication of each selective coefficient, ranging between 0 to 1, at each resistance mutation. Therefore, multiple resistance mutations are difficult to observe without very high selective coefficients, and we argue that invasive populations are under strong selective pressure for the resistance against synthetic insecticides.
Glutathione S transferases or UDP glycosyltransferase might also cause resistance against synthetic insecticides. However, the direct relationship between these two gene families and resistance is yet to be tested in FAW. Genomic analyses and bioassay experiments will generate a more comprehensive list of causal resistance mutations against synthetic insecticides.
The differential usage of synthetic insecticides may explain the observations that invasive populations have increased allele frequency of resistance mutations at AchE genes and increased number of P450 genes by CNVs. The average annual insecticide use is particularly high in East Asia including China [46]. We speculate that initially introduced FAW insects in East Asia might have survived through adaptive evolution by the increased P450 gene number or the increased allele frequency of resistance mutations at AchE genes, while these FAWs remained undetected due to the low population sizes before explosive population expansion in West Africa, where genetic admixture occurred between populations with different invasive origins [31]. More studies are necessary to test the role of adaptive evolution associated with synthetic insecticide resistance in invasive success.
The FST result supports natural selection for Bt resistance in the Brazilian population. If selective pressure with the same directionality exists in invasive populations, newly generated Bt resistance mutations can be selected in the near future. Theoretically, species with large population sizes have higher adaptive potential than species with small population sizes because of a higher population-scaled rate of mutation [47,48]. Since the estimated population size of FAWs is as high as 3.75 million [49], the population-scaled mutation rate can be sufficient to generate a genetic variation for drastically rapid Darwinian positive selection. We observed two synonymous codons at a resistance position in AChE from a population Benin, implying that mutational influx could not be a limiting factor of adaptive evolution through positive selection. Therefore, invasive populations might experience adaptive evolution, causing resistance against Bt mutations in the near future.

5. Conclusions

In this study, we showed that Bt resistance mutations did not spread during this period to the Old World by the surveyed invasive FAW populations and that invasive populations have increased copy numbers of P450 genes and increased proportions of resistance mutations at AchE compared with native populations. These results imply that invasive populations might be under selective pressure for the resistance against synthetic insecticides.
The samples used in this study were obtained between 2017 to 2018. Since the first reported FAWs in Western Africa in 2016 [31], the genetic information in this study represents snapshots at the early phase of invasions. As Bt crops and synthetic insecticides are commonly used in the invasive area, the level of insecticide resistance can be changed together with underlying genotypes. This change is likely because our result is in line with selective pressure both for Bt and, in particular, synthetic insecticides. Therefore, monitoring changes in field susceptibility to various insecticides will be helpful in FAW control, together with the analysis of changes in genomics sequences to identify known or new causal resistance mutations.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/insects12050468/s1 (Table S1: Sample collection site and year, Table S2: The list of analyzed P450 genes in this study).

Author Contributions

Conceptualization, S.Y., T.W., and K.N.; Formal analysis, S.Y. and K.N.; Resources, N.N., P.J.S., T.B., W.T.T., K.G., an E.d.; Writing—original draft preparation, S.Y. and K.N.; Writing—review and editing, N.N., P.J.S., T.B., W.T.T., K.G., E.d., and T.W.; Supervision, K.N., T.B., and E.d.; Funding acquisition, K.N. All authors have read and agreed to the published version of the manuscript.

Funding

This work (ID 1702-018, given to K.N.) was publicly funded through ANR (the French National Research Agency) under the “Investissements d’avenir” program with the reference ANR-10-LABX-001-01 Labex Agro and coordinated by Agropolis Fondation under the frame of I-SITE MUSE (ANR-16-IDEX-0006). In addition, a grant from the department of Santé des Plantes et Environnement at Institut national de la recherche agronomique was awarded to K.N. (adaptivesv). This work was also financially supported by EUPHRESCO (FAW-spedcom, given to Anne-Nathalie Volkoff) and by CSIRO Health & Biosecurity (given to W.T.T., T.W., and K.G.). S.Y. was supported by a CIRAD-INRAE PhD fellowship.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

We will provide all files generated in this study upon request. These files include vcf files, gene annotation files, and computer programming scripts.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Tabashnik, B.E.; Brévault, T.; Carrière, Y. Insect Resistance to Bt Crops: Lessons from the First Billion Acres. Nat. Biotechnol. 2013, 31, 510–521. [Google Scholar] [CrossRef]
  2. Jurat-Fuentes, J.L.; Heckel, D.G.; Ferré, J. Mechanisms of Resistance to Insecticidal Proteins from Bacillus thuringiensis. Annu. Rev. Entomol. 2021, 66, 121–140. [Google Scholar] [CrossRef]
  3. Hemingway, J.; Field, L.; Vontas, J. An Overview of Insecticide Resistance. Science 2002, 298, 96–97. [Google Scholar] [CrossRef]
  4. Carvalho, R.A.; Omoto, C.; Field, L.M.; Williamson, M.S.; Bass, C. Investigating the Molecular Mechanisms of Organophosphate and Pyrethroid Resistance in the Fall Armyworm Spodoptera frugiperda. PLoS ONE 2013, 8, e62268. [Google Scholar] [CrossRef] [Green Version]
  5. Shan, T.; Chen, C.; Ding, Q.; Chen, X.; Zhang, H.; Chen, A.; Shi, X.; Gao, X. Molecular Characterization and Expression Profiles of Nicotinic Acetylcholine Receptors in Bradysia odoriphaga. Pestic. Biochem. Physiol. 2020, 165, 104563. [Google Scholar] [CrossRef] [PubMed]
  6. Kim, A.-Y.; Kwon, D.H.; Jeong, I.H.; Koh, Y.H. An Investigation of the Molecular and Biochemical Basis Underlying Chlorantraniliprole-Resistant Drosophila Strains and Their Cross-Resistance to Other Insecticides. Arch. Insect Biochem. Physiol. 2018, 99, e21514. [Google Scholar] [CrossRef]
  7. Lin, L.; Hao, Z.; Cao, P.; Yuchi, Z. Homology Modeling and Docking Study of Diamondback Moth Ryanodine Receptor Reveals the Mechanisms for Channel Activation, Insecticide Binding and Resistance. Pest Manag. Sci. 2020, 76, 1291–1303. [Google Scholar] [CrossRef]
  8. Montezano, D.G.; Specht, A.; Sosa-Gómez, D.R.; Roque-Specht, V.F.; Sousa-Silva, J.C.; de Paula-Moraes, S.V.; Peterson, J.A.; Hunt, T.E. Host Plants of Spodoptera frugiperda (Lepidoptera: Noctuidae) in the Americas. Afr. Entomol. 2018, 26, 286–300. [Google Scholar] [CrossRef] [Green Version]
  9. Gouin, A.; Bretaudeau, A.; Nam, K.; Gimenez, S.; Aury, J.-M.; Duvic, B.; Hilliou, F.; Durand, N.; Montagné, N.; Darboux, I.; et al. Two Genomes of Highly Polyphagous Lepidopteran Pests (Spodoptera frugiperda, Noctuidae) with Different Host-Plant Ranges. Sci. Rep. 2017, 7, 11816. [Google Scholar] [CrossRef]
  10. Gui, F.; Lan, T.; Zhao, Y.; Guo, W.; Dong, Y.; Fang, D.; Liu, H.; Li, H.; Wang, H.; Hao, R. Genomic and Transcriptomic Analysis Unveils Population Evolution and Development of Pesticide Resistance in Fall Armyworm Spodoptera frugiperda. Protein Cell 2020. [Google Scholar] [CrossRef]
  11. Xiao, H.; Ye, X.; Xu, H.; Mei, Y.; Yang, Y.; Chen, X.; Yang, Y.; Liu, T.; Yu, Y.; Yang, W. The Genetic Adaptations of Fall Armyworm Spodoptera frugiperda Facilitated Its Rapid Global Dispersal and Invasion. Mol. Ecol. Resour. 2020, 20, 1050–1068. [Google Scholar] [CrossRef] [PubMed]
  12. Jakka, S.R.K.; Gong, L.; Hasler, J.; Banerjee, R.; Sheets, J.J.; Narva, K.; Blanco, C.A.; Jurat-Fuentes, J.L. Field-Evolved Mode 1 Resistance of the Fall Armyworm to Transgenic Cry1Fa-Expressing Corn Associated with Reduced Cry1Fa Toxin Binding and Midgut Alkaline Phosphatase Expression. Appl. Environ. Microbiol. 2016, 82, 1023–1034. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Storer, N.P.; Babcock, J.M.; Schlenz, M.; Meade, T.; Thompson, G.D.; Bing, J.W.; Huckaba, R.M. Discovery and Characterization of Field Resistance to Bt Maize: Spodoptera frugiperda (Lepidoptera: Noctuidae) in Puerto Rico. J. Econ. Entomol. 2010, 103, 1031–1038. [Google Scholar] [CrossRef]
  14. Monnerat, R.; Martins, E.; Macedo, C.; Queiroz, P.; Praça, L.; Soares, C.M.; Moreira, H.; Grisi, I.; Silva, J.; Soberon, M.; et al. Evidence of Field-Evolved Resistance of Spodoptera frugiperda to Bt Corn Expressing Cry1F in Brazil That Is Still Sensitive to Modified Bt Toxins. PLoS ONE 2015, 10, e0119544. [Google Scholar] [CrossRef] [Green Version]
  15. Omoto, C.; Bernardi, O.; Salmeron, E.; Sorgatto, R.J.; Dourado, P.M.; Crivellari, A.; Carvalho, R.A.; Willse, A.; Martinelli, S.; Head, G.P. Field-Evolved Resistance to Cry1Ab Maize by Spodoptera frugiperda in Brazil. Pest Manag. Sci. 2016, 72, 1727–1736. [Google Scholar] [CrossRef] [PubMed]
  16. Chandrasena, D.I.; Signorini, A.M.; Abratti, G.; Storer, N.P.; Olaciregui, M.L.; Alves, A.P.; Pilcher, C.D. Characterization of Field-Evolved Resistance to Bacillus Thuringiensis-Derived Cry1F δ-Endotoxin in Spodoptera frugiperda Populations from Argentina. Pest Manag. Sci. 2018, 74, 746–754. [Google Scholar] [CrossRef] [Green Version]
  17. Schlum, K.A.; Lamour, K.; de Bortoli, C.P.; Banerjee, R.; Meagher, R.; Pereira, E.; Murua, M.G.; Sword, G.A.; Tessnow, A.E.; Viteri Dillon, D.; et al. Whole Genome Comparisons Reveal Panmixia among Fall Armyworm (Spodoptera frugiperda) from Diverse Locations. BMC Genom. 2021, 22, 179. [Google Scholar] [CrossRef]
  18. Banerjee, R.; Hasler, J.; Meagher, R.; Nagoshi, R.; Hietala, L.; Huang, F.; Narva, K.; Jurat-Fuentes, J.L. Mechanism and DNA-Based Detection of Field-Evolved Resistance to Transgenic Bt Corn in Fall Armyworm (Spodoptera frugiperda). Sci. Rep. 2017, 7, 10877. [Google Scholar] [CrossRef] [Green Version]
  19. Flagel, L.; Lee, Y.W.; Wanjugi, H.; Swarup, S.; Brown, A.; Wang, J.; Kraft, E.; Greenplate, J.; Simmons, J.; Adams, N.; et al. Mutational Disruption of the ABCC2 Gene in Fall Armyworm, Spodoptera frugiperda, Confers Resistance to the Cry1Fa and Cry1A.105 Insecticidal Proteins. Sci. Rep. 2018, 8, 7255. [Google Scholar] [CrossRef] [Green Version]
  20. Boaventura, D.; Ulrich, J.; Lueke, B.; Bolzan, A.; Okuma, D.; Gutbrod, O.; Geibel, S.; Zeng, Q.; Dourado, P.M.; Martinelli, S.; et al. Molecular Characterization of Cry1F Resistance in Fall Armyworm, Spodoptera frugiperda from Brazil. Insect Biochem. Mol. Biol. 2020, 116, 103280. [Google Scholar] [CrossRef]
  21. Guan, F.; Zhang, J.; Shen, H.; Wang, X.; Padovan, A.; Walsh, T.K.; Tay, W.T.; Gordon, K.H.J.; James, W.; Czepak, C.; et al. Whole-Genome Sequencing to Detect Mutations Associated with Resistance to Insecticides and Bt Proteins in Spodoptera frugiperda. Insect Sci. 2020. [Google Scholar] [CrossRef] [PubMed]
  22. Gutiérrez-Moreno, R.; Mota-Sanchez, D.; Blanco, C.A.; Whalon, M.E.; Terán-Santofimio, H.; Rodriguez-Maciel, J.C.; DiFonzo, C. Field-Evolved Resistance of the Fall Armyworm (Lepidoptera: Noctuidae) to Synthetic Insecticides in Puerto Rico and Mexico. J. Econ. Entomol. 2019, 112, 792–802. [Google Scholar] [CrossRef] [PubMed]
  23. Gimenez, S.; Abdelgaffar, H.; Goff, G.L.; Hilliou, F.; Blanco, C.A.; Hänniger, S.; Bretaudeau, A.; Legeai, F.; Nègre, N.; Jurat-Fuentes, J.L.; et al. Adaptation by Copy Number Variation Increases Insecticide Resistance in the Fall Armyworm. Commun. Biol. 2020, 3, 664. [Google Scholar] [CrossRef] [PubMed]
  24. Boaventura, D.; Bolzan, A.; Padovez, F.E.; Okuma, D.M.; Omoto, C.; Nauen, R. Detection of a Ryanodine Receptor Target-Site Mutation in Diamide Insecticide Resistant Fall Armyworm, Spodoptera frugiperda. Pest Manag. Sci. 2020, 76, 47–54. [Google Scholar] [CrossRef]
  25. Goergen, G.; Kumar, P.L.; Sankung, S.B.; Togola, A.; Tamò, M. 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 2016, 11, e0165632. [Google Scholar] [CrossRef] [Green Version]
  26. Day, R.; Abrahams, P.; Bateman, M.; Beale, T.; Clottey, V.; Cock, M.; Colmenarez, Y.; Corniani, N.; Early, R.; Godwin, J. Fall Armyworm: Impacts and Implications for Africa. Outlooks Pest Manag. 2017, 28, 196–201. [Google Scholar] [CrossRef] [Green Version]
  27. Pashley, D.P. Host-Associated Genetic Differentiation in Fall Armyworm (Lepidoptera: Noctuidae): A Sibling Species Complex? Ann. Entomol. Soc. Am. 1986, 79, 898–904. [Google Scholar] [CrossRef]
  28. Pashley, D.P. Host-associated differentiation in armyworms (Lepidoptera: Noctuidae): An allozymic and mitochondrial DNA perspective. In Electrophoretic Studies on Agricultural Pests; Loxdale, H.D., Hollander, J.D., Eds.; Clarendon Press: Oxford, UK, 1989; Volume 39, pp. 103–144. [Google Scholar]
  29. Dumas, P.; Legeai, F.; Lemaitre, C.; Scaon, E.; Orsucci, M.; Labadie, K.; Gimenez, S.; Clamens, A.-L.; Henri, H.; Vavre, F.; et al. Spodoptera frugiperda (Lepidoptera: Noctuidae) Host-Plant Variants: Two Host Strains or Two Distinct Species? Genetica 2015, 143, 305–316. [Google Scholar] [CrossRef] [Green Version]
  30. Nagoshi, R.N.; Goergen, G.; Plessis, H.D.; van den Berg, J.; Meagher, R. Genetic Comparisons of Fall Armyworm Populations from 11 Countries Spanning Sub-Saharan Africa Provide Insights into Strain Composition and Migratory Behaviors. Sci. Rep. 2019, 9, 8311. [Google Scholar] [CrossRef] [Green Version]
  31. Yainna, S.; Tay, W.T.; Fiteni, E.; Legeai, F.; Clamens, A.-L.; Gimenez, S.; Frayssinet, M.; Asokan, R.; Kalleshwaraswamy, C.M.; Deshmukh, S.; et al. Genomic Balancing Selection Is Key to the Invasive Success of the Fall Armyworm. bioRxiv 2020. [Google Scholar] [CrossRef]
  32. Zhang, L.; Liu, B.; Zheng, W.; Liu, C.; Zhang, D.; Zhao, S.; Li, Z.; Xu, P.; Wilson, K.; Withers, A.; et al. Genetic Structure and Insecticide Resistance Characteristics of Fall Armyworm Populations Invading China. Mol. Ecol. Resour. 2020, 20, 1682–1696. [Google Scholar] [CrossRef] [PubMed]
  33. Tay, W.T.; Rane, R.; Padovan, A.; Walsh, T.; Elfekih, S.; Downes, S.; Nam, K.; d’Alençon, E.; Zhang, J.; Wu, Y.; et al. Whole Genome Sequencing of Global Spodoptera frugiperda Populations: Evidence for Complex, Multiple Introductions across the Old World. bioRxiv 2020. [Google Scholar] [CrossRef]
  34. Boaventura, D.; Martin, M.; Pozzebon, A.; Mota-Sanchez, D.; Nauen, R. Monitoring of Target-Site Mutations Conferring Insecticide Resistance in Spodoptera frugiperda. Insects 2020, 11, 545. [Google Scholar] [CrossRef]
  35. McKenna, A.; Hanna, M.; Banks, E.; Sivachenko, A.; Cibulskis, K.; Kernytsky, A.; Garimella, K.; Altshuler, D.; Gabriel, S.; Daly, M.; et al. The Genome Analysis Toolkit: A MapReduce Framework for Analyzing next-Generation DNA Sequencing Data. Genome Res. 2010, 20, 1297–1303. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Wang, X.; Zheng, Z.; Cai, Y.; Chen, T.; Li, C.; Fu, W.; Jiang, Y. CNVcaller: Highly Efficient and Widely Applicable Software for Detecting Copy Number Variations in Large Populations. GigaScience 2017, 6. [Google Scholar] [CrossRef]
  37. Create Your Own Custom Map. Available online: https://mapchart.net/index.html (accessed on 17 May 2021).
  38. Slater, G.S.C.; Birney, E. Automated Generation of Heuristics for Biological Sequence Comparison. BMC Bioinform. 2005, 6, 31. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Robinson, J.T.; Thorvaldsdóttir, H.; Winckler, W.; Guttman, M.; Lander, E.S.; Getz, G.; Mesirov, J.P. Integrative Genomics Viewer. Nat. Biotechnol. 2011, 29, 24–26. [Google Scholar] [CrossRef] [Green Version]
  40. Camacho, C.; Coulouris, G.; Avagyan, V.; Ma, N.; Papadopoulos, J.; Bealer, K.; Madden, T.L. BLAST+: Architecture and Applications. BMC Bioinform. 2009, 10, 421. [Google Scholar] [CrossRef] [Green Version]
  41. Weir, B.S.; Cockerham, C.C. Estimating F-Statistics for the Analysis of Population Structure. Evolution 1984, 38, 1358–1370. [Google Scholar]
  42. Danecek, P.; Auton, A.; Abecasis, G.; Albers, C.A.; Banks, E.; DePristo, M.A.; Handsaker, R.E.; Lunter, G.; Marth, G.T.; Sherry, S.T.; et al. The Variant Call Format and VCFtools. Bioinformatics 2011, 27, 2156–2158. [Google Scholar] [CrossRef]
  43. Walsh, T.K.; Joussen, N.; Tian, K.; McGaughran, A.; Anderson, C.J.; Qiu, X.; Ahn, S.-J.; Bird, L.; Pavlidi, N.; Vontas, J.; et al. Multiple Recombination Events between Two Cytochrome P450 Loci Contribute to Global Pyrethroid Resistance in Helicoverpa armigera. PLoS ONE 2018, 13, e0197760. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Valencia-Montoya, W.A.; Elfekih, S.; North, H.L.; Meier, J.I.; Warren, I.A.; Tay, W.T.; Gordon, K.H.; Specht, A.; Paula-Moraes, S.V.; Rane, R. Adaptive Introgression across Semipermeable Species Boundaries between Local Helicoverpa zea and Invasive Helicoverpa armigera Moths. Mol. Biol. Evol. 2020, 37, 2568–2583. [Google Scholar] [CrossRef] [PubMed]
  45. Zhang, D.; Xiao, Y.; Xu, P.; Yang, X.; Wu, Q.; Wu, K. Insecticide Resistance Monitoring for the Invasive Populations of Fall Armyworm, Spodoptera frugiperda in China. J. Integr. Agric. 2021, 20, 783–791. [Google Scholar] [CrossRef]
  46. Zhang, W. Global Pesticide Use: Profile, Trend, Cost/Benefit and More. Proc. Int. Acad. Ecol. Environ. Sci. 2018, 8, 1–27. [Google Scholar]
  47. Lanfear, R.; Kokko, H.; Eyre-Walker, A. Population Size and the Rate of Evolution. Trends Ecol. Evol. 2014, 29, 33–41. [Google Scholar] [CrossRef] [PubMed]
  48. Nam, K.; Munch, K.; Mailund, T.; Nater, A.; Greminger, M.P.; Krützen, M.; Marquès-Bonet, T.; Schierup, M.H. Evidence That the Rate of Strong Selective Sweeps Increases with Population Size in the Great Apes. Proc. Natl. Acad. Sci. USA 2017, 114, 1613–1618. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Nam, K.; Nhim, S.; Robin, S.; Bretaudeau, A.; Nègre, N.; d’Alençon, E. Positive Selection Alone Is Sufficient for Whole Genome Differentiation at the Early Stage of Speciation Process in the Fall Armyworm. BMC Evol. Biol. 2020, 20, 152. [Google Scholar] [CrossRef]
Figure 1. Countries where the analyzed individuals were sampled. The green color indicates native populations. Purple (2016), orange (2017), and yellow colors (2018) indicate the reported years of detection in countries from where individuals were collected. The map was generated using mapchart [37].
Figure 1. Countries where the analyzed individuals were sampled. The green color indicates native populations. Purple (2016), orange (2017), and yellow colors (2018) indicate the reported years of detection in countries from where individuals were collected. The map was generated using mapchart [37].
Insects 12 00468 g001
Figure 2. Full amino acid sequence of the ABCC2 gene with information on all mutations studied. The green and orange arrows indicate the identified and unidentified mutations from the resequencing data, respectively.
Figure 2. Full amino acid sequence of the ABCC2 gene with information on all mutations studied. The green and orange arrows indicate the identified and unidentified mutations from the resequencing data, respectively.
Insects 12 00468 g002
Figure 3. FST calculated between CC and rCC (red vertical line) and from random grouping among CC and rCC.
Figure 3. FST calculated between CC and rCC (red vertical line) and from random grouping among CC and rCC.
Insects 12 00468 g003
Figure 4. Average numbers of deleted or duplicated P450 genes per individual in invasive and native populations for clan2, clan3, and clan4.
Figure 4. Average numbers of deleted or duplicated P450 genes per individual in invasive and native populations for clan2, clan3, and clan4.
Insects 12 00468 g004
Figure 5. The proportion of individuals with susceptibility and resistance at AChE for (A) F209V, (B) A201S, and (C) G227A in invasive and native populations.
Figure 5. The proportion of individuals with susceptibility and resistance at AChE for (A) F209V, (B) A201S, and (C) G227A in invasive and native populations.
Insects 12 00468 g005
Table 1. Individuals with observed resistant mutations (* indicates individuals heterozygous for the mutations). PR: Puerto Rico, CC: susceptible Brazilian individuals, rCC: resistant Brazilian individuals.
Table 1. Individuals with observed resistant mutations (* indicates individuals heterozygous for the mutations). PR: Puerto Rico, CC: susceptible Brazilian individuals, rCC: resistant Brazilian individuals.
MutationsIndividuals with Resistance Mutations
2 bp insertion and a frameshift mutation* PR1, PR5, PR12, PR14, PR16, * PR18, * PR19, * PR27, PR30, * PR31, * PR33
GY deletionCC44, CC69
P799K/RCC44, CC69
12 bp insertion and a frameshift mutation* rCC25, * rCC5
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Yainna, S.; Nègre, N.; Silvie, P.J.; Brévault, T.; Tay, W.T.; Gordon, K.; dAlençon, E.; Walsh, T.; Nam, K. Geographic Monitoring of Insecticide Resistance Mutations in Native and Invasive Populations of the Fall Armyworm. Insects 2021, 12, 468. https://doi.org/10.3390/insects12050468

AMA Style

Yainna S, Nègre N, Silvie PJ, Brévault T, Tay WT, Gordon K, dAlençon E, Walsh T, Nam K. Geographic Monitoring of Insecticide Resistance Mutations in Native and Invasive Populations of the Fall Armyworm. Insects. 2021; 12(5):468. https://doi.org/10.3390/insects12050468

Chicago/Turabian Style

Yainna, Sudeeptha, Nicolas Nègre, Pierre J. Silvie, Thierry Brévault, Wee Tek Tay, Karl Gordon, Emmanuelle dAlençon, Thomas Walsh, and Kiwoong Nam. 2021. "Geographic Monitoring of Insecticide Resistance Mutations in Native and Invasive Populations of the Fall Armyworm" Insects 12, no. 5: 468. https://doi.org/10.3390/insects12050468

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop