Next Article in Journal
Global Potential Distribution of Invasive Species Pseudococcus viburni (Hemiptera: Pseudococcidae) under Climate Change
Previous Article in Journal
RNAi-Mediated Silencing of Laccase 2 in Culex pipiens Pupae via Dehydration and Soaking Results in Multiple Defects in Cuticular Development
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Insects
Volume 15
Issue 3
10.3390/insects15030194
insects-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Presence of Multiple Genetic Mutations Related to Insecticide Resistance in Chinese Field Samples of Two Phthorimaea Pest Species

by
Jiang Zhu
1,2,
Ruipeng Chen
2,
Juan Liu
1,
Weichao Lin
1,
Jiaxin Liang
1,
Ralf Nauen
3,
Suhua Li
1,* and
Yulin Gao
2,*
1
Shenzhen Branch, Guangdong Laboratory of Lingnan Modern Agriculture, Key Laboratory of Synthetic Biology, Ministry of Agriculture and Rural Affairs, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen 518000, China
2
State Key Laboratory for Biology of Plant Diseases and Insect Pests, Institute of Plant Protection, Chinese Academy of Agricultural Sciences, Beijing 100193, China
3
Crop Science Division, R&D, Pest Control, Bayer AG, Alfred-Nobel-Str. 50, 40789 Monheim, Germany
*
Authors to whom correspondence should be addressed.
Insects 2024, 15(3), 194; https://doi.org/10.3390/insects15030194
Submission received: 2 February 2024 / Revised: 4 March 2024 / Accepted: 11 March 2024 / Published: 14 March 2024
(This article belongs to the Section Insect Molecular Biology and Genomics)
Browse Figures
Review Reports Versions Notes

Abstract

:

Simple Summary

Insecticide resistance stands as one of the most notable evolutionary phenomena for researchers. Two closely related species pests, Phthorimaea operculella and Phthorimaea absoluta, both feeding on potato crops, have developed distinct resistance mechanisms. In this study, we compared the presence of target-site mutations in P. operculella and P. absoluta in three common insecticide targets, Ace1 (acetylcholinesterase), VGSC (voltage-gated sodium channel), and RDL (GABA-gated chloride channel). Unexpected genetic divergence in target-site mutations was observed between the two species. P. absoluta had three Ace1 mutations (A201S, L231V, and F290V), four VGSC mutations (M918T, L925M, T928I, and L1014F), and one RDL mutation (A301S). On the other hand, P. operculella had Ace1 mutations (F158Y, A201S, and L231V) and only one VGSC mutation (L1014F) at lower frequencies, with no mutations detected in RDL. These findings deepen our understanding of evolutionary dynamics in pest species and offer potential strategies for more precise and sustainable pest control in potatoes.

Abstract

Potatoes hold the distinction of being the largest non-cereal food crop globally. The application of insecticides has been the most common technology for pest control. The repeated use of synthetic insecticides of the same chemical class and frequent applications have resulted in the emergence of insecticide resistance. Two closely related pests that feed on potato crops are the potato tuber moth, Phthorimaea operculella, and the tomato leafminer, Phthorimaea absoluta (syn. Tuta absoluta). Previous studies indicated the existence of insecticide resistance to various classes of insecticides including organophosphates, carbamates, and pyrethroids in field populations of P. operculella and P. absoluta. However, the exact mechanisms of insecticide resistance in P. operculella and to a lesser extent P. absoluta remain still poorly understood. Detecting resistance genotypes is crucial for the prediction and management of insecticide resistance. In this study, we identified multiple genetic mutations related to insecticide resistance in two species of Phthorimaea. An unexpected genetic divergence on target-site mutations was observed between P. operculella and P. absoluta. Three mutations (A201S, L231V, and F290V) in Ace1 (acetylcholinesterase), four mutations (M918T, L925M, T928I, and L1014F) in VGSC (voltage-gated sodium channel), and one mutation (A301S) in RDL (GABA-gated chloride channel) have been detected with varying frequencies in Chinese P. absoluta field populations. In contrast, P. operculella field populations showed three mutations (F158Y, A201S, and L231V) in Ace1, one mutation (L1014F) in VGSC at a lower frequency, and no mutation in RDL. These findings suggest that pyrethroids, organophosphates, and carbamates are likely to be ineffective in controlling P. absoluta, but not P. operculella. These findings contributed to a deeper understanding of the presence of target-site mutations conferring resistance to commonly used (and cheap) classes of insecticides in two closely related potato pests. It is recommended to consider the resistance status of both pests for the implementation of resistance management strategies in potatoes.

1. Introduction

Potatoes (Solanum tuberosum L.) rank as the fourth most important food crop globally, next to rice, wheat, and maize [1]. It is a critical crop and widely cultivated in terms of its rich nutritive value. China, the biggest producer of potatoes since 1993 [2], is playing an increasing role in the global potato market. FAO data show that it continues to increase in its importance with an annual production of potatoes in China up to approximately 10 million tons in 2022 (http://www.fao.org/statistics/en/, accessed on 13 March 2023). Most of the potato production takes place in the northern and southwest regions of China, notably in Yunnan province.
The potato tuber moth, Phthorimaea operculella (Zeller) (Lepidoptera: Gelechiidae), is considered one of the most devastating insect pests on potatoes, feeding on roots, stems, foliage, and, more importantly, potato tubers. Severe infestations result in substantial yield and quality losses, even up to 100% losses to potato plants in fields as well as in storage [3]. P. operculella originated from the tropical mountainous region of South America and then distributed over a wide range of tropical and subtropical countries in North America, Europe, Africa, and Asia [4,5]. It was first introduced to China in 1937 and found on tobacco in Yulin, Guangxi province [6]. To our knowledge, it is widely distributed in more than twenty provinces in China, including Yunnan, Guizhou, and Sichuan [7,8].
The tomato leafminer, Phthorimaea absoluta (Meyrick) (Lepidoptera: Gelechiidae), is also known as Tuta absoluta, another important invasive pest within the Gelechiidae family, causing destructive harm to Solanaceous crops including potatoes. The larvae primarily feed on the mesophyll of all aerial parts of the potato plant, leading to leaf damage and plant withering, and resulting in extreme yield losses and economic damage [9]. It originated from Peru and quickly spread across Europe, Africa, and Asia in the past decades [10,11,12,13], and is now found in more than 110 countries and regions [14]. Since its first detection in Xinjiang province in 2017 [15], P. absoluta has subsequently spread throughout China and has been detected in more than ten provinces, especially in Yunnan province [16,17].
P.operculella and P. absoluta are two closely related species that feed on potato plants, and it is difficult to distinguish them by morphological characteristics. Their larvae feed on the tubers and/or foliage of potato crops, which causes cumulative damage. Both pests can occur on potato crops at the same time, rendering control methods to keep them under economic damage thresholds important. The cultivation of dominant potato varieties [18], the use of sex pheromones [19], and deeper seed planting are recommended to prevent pest infestation. However, the application of chemical insecticides has been the mainstay of controlling P. operculella and P. absoluta in diverse agricultural production systems [8,14]. To date, only cypermethrin and lufenuron have been registered to control P. operculella in China [20], with no options available to control P. absoluta. Several applied problems, such as poor chemical selection and substandard application practices, have exacerbated control failures and contribute to the evolution of a variety of resistance mechanisms in response to insecticide selection. Frequent applications of insecticides of the same chemical class or mode of action potentially select for insecticide resistance [21].
Current knowledge supports that insecticide resistance is primarily driven by metabolic enzymes (i.e., cytochrome P450-monooxygenases (P450s)) and target-site mutations affecting insecticide binding such as mutations in Ace1 (acetylcholinesterase), VGSC (voltage-gated sodium channels), and RDL (GABA-gated chloride channel) [22,23,24], conferring resistance to organophosphates, carbamates, pyrethroids, and fiproles, respectively. Metabolic and target-site resistance may be observed simultaneously or isolated in a population under continuous selection pressure. Insecticide resistance mediated using P450s is highly prevalent and has been extensively documented in association with metabolic detoxification under insecticide selection [24,25,26]. Simultaneously, the studies on target-site mutations conferring insecticide resistance have seen increasing importance in recent decades [27,28,29,30,31,32,33,34] and are considered of major concern. It has been reported that mutations in Ace1 (G247S [28], A302S [29], F402V [30], F439W/Y [31]) lead to resistance against organophosphates (OPs) and carbamates (CBs). Mutations in VGSC (V410L [32], M918I + L1014F [33]) have been shown to confer resistance to pyrethroids (PYs). Additionally, mutations in RDL (A301S [34]) have been linked to resistance against dieldrin and fipronil.
Previous studies revealed the existence of insecticide resistance to organophosphates, carbamates, pyrethroids, and diamides in P. operculella and P. absoluta field populations [6,11,35,36,37,38,39,40,41]. However, there are much fewer reports on the molecular basis of insecticide resistance in P. operculella than in P. absoluta. As these two closely related species both cause damage to potato crops, this study detected their multiple genetic mutations related to insecticide resistance by PCR and investigated the distribution and frequencies of mutations in the southwest and northwest regions of China. We attempted to explore the resistance mechanism in field-collected populations of P. operculella and P. absoluta with a specific focus on target-site mutations conferring resistance to commonly used chemical classes of insecticides, including OPs, CBs, PYs, and fiproles. These data provide a basis for developing resistance management strategies in potatoes, covering both pests under applied conditions.

2. Methods

2.1. Insects

P. operculella and P. absoluta larvae were collected at six different locations in China during April and August 2023 (Table S3). Field-collected insects were stored in 100% ethanol, brought to the laboratory, and kept at −20 °C until use. The species determination was based on PCR and sequencing according to Kwon et al. [42].

2.2. Genomic DNA Extraction

Thirty individuals from each sampling site were selected randomly. The genomic DNA (gDNA) of individual larvae was prepared according to the method of Rinkevich et al. [43]. Briefly, the larva was placed in a tube with 500 µL of lysis buffer containing 100 mM Tris-HCl (pH = 8.0), 50 mM NaCl, 10 mM EDTA, 1% (w/v) SDS, and 0.1 mg/mL proteinase K, homogenized using a plastic pestle, and then incubated at 56 °C for 1 h. The mixture was set in an ice bath for 10 min after the addition of 50 µL 8 M potassium acetate, and then centrifuged at 12,000× g for 10 min. The supernatant was then transferred to a new tube and the gDNA was precipitated by adding ice-cold ethanol (2:1) at room temperature for 10 min. The pellet was washed with 75% ethanol and resuspended in ddH2O. The prepared gDNA samples were stored at −20 °C.

2.3. PCR Amplification of Target Gene Fragments and Sanger Sequencing

Fragments of the PoVGSC (~2500 bp), PoAce1 (~3500 bp), PoRdl (~650 bp), PaVGSC (~1600 bp), PaAce1 (~3500 bp), and PaRdl (~1000 bp) genes were amplified by PCR using primer pairs shown in Table S1. The PCR reactions were performed with 2× Phanta Max Master Mix (Vazyme Biotech Co., Ltd., Nanjing, China) in a T-100 thermal cycler (Bio-Rad, Hercules, CA, USA) as follows: 35 cycles of 5 s at 98 °C, 10 s at 55 °C, and 30 s at 72 °C. The purified PCR products from heterozygotes were ligated with the pUC19 vector and transformed into competent cells of the Escherichia coli DH5α strain (Shenzhen Kangti Life Technology Co., Ltd., Shenzhen, China). The PCR products from homozygotes and three clones from heterozygotes were Sanger-sequenced by Sangon Biotech (Guangzhou, China).

2.4. DNA Sequence Analysis and Phylogenetic Analysis

DNA sequencing data were checked and curated manually. All confirmed sequences from filed populations, combined with data from the NCBI database, were aligned with MUSCLE [44]. The haplotypes were constructed using DnaSPv5 software from the alignment data [45]. The genealogical relations among haplotypes were analyzed using Networkv10 [46]. The phylogenetic tree was generated using MEGA X [47] with the maximum likelihood method and the Tamura 3-parameter model with 1000 bootstrap replicates based on the DNA alignment.

2.5. Molecular Clock Estimates

In order to examine the relationships among Gelechiidae, the mitochondrial cytochrome c oxidase subunit I (COI) gene sequences were sourced from the GenBank database (Table S2). The divergence time estimation was performed using the RelTime method following a standard protocol as implemented in MEGA X under the GTR + G model with 1000 bootstrap replicates [47]. Drosophila melanogaster was set as an outgroup. There were a total of 1531 positions in the final dataset. The calibration was set to 101.9–125.2 MYA (Phthorimaea operculella versus Helicoverpa armigera), 79.8–127.3 MYA (Phthorimaea operculella versus Plutella xylostella), and 38–60.3 MYA (Helicoverpa armigera versus Spodoptera frugiperda) following the TimeTree database (http://timetree.org/, accessed on 12 February 2024). The other parameters were defined by default.

2.6. Syntenic Relationship of the Insecticide Target Gene Loci

To detect the collinearity of target genes between P. operculella and P. absoluta, the local comparative genomics analysis was conducted with TBtool software using the Genome Region Compare Suite plugin with the following blastp parameters: maxHsp500 and minLen100 [48]. The complete genome of P. operculella (GenBank assembly number: GCA_024500475.1) and P. absoluta (GenBank assembly number: GCA_027580185.1), along with their annotation information, were downloaded from the NCBI dataset (https://www.ncbi.nlm.nih.gov/genome/, accessed on 2 February 2024).

3. Results

3.1. Distribution and Frequency of Insecticide Resistance-Conferring Mutations

Three target genes expressing commonly known insecticide targets were analyzed for mutations conferring insecticide resistance. Based on previous studies, the mutations G119S, A201S, G328A, and F331W (Torpedo californica numbering) in Ace1 were particularly analyzed in amplified sequences containing Ace1 codons 119–331 of P. operculella (PoAce1) and P. absoluta (PaAce1). An amplified fragment of around 700 bp was analyzed from each species and, in total, four mutation sites were detected, i.e., F158Y, A201S, L231V, and F290V, including two formerly validated Ace1 target-site mutations (A201S and F290V). The frequencies of the detected mutant alleles are detailed in Table 1. Interestingly, the F158Y allele was exclusively detected in Ace1 of P. operculella, whereas the two known resistance mutations (A201S and F290V) were predominantly detected in P. absoluta. The commonly known Ace1 mutations G119S and F331W were not detected in field populations of either P. operculella or P. absoluta.
The VGSC gene of P. operculella and P. absoluta contains 31 and 34 exons, respectively, and encodes the target for pyrethroid insecticides. Here, we analyzed an amplified region (~2000 bp) containing codons 918 and 1014 of VGSC for the presence of mutations formerly described to confer pyrethroid resistance. Four common mutations (M918T, L925M, T929I, and L1014F, Musca domestica numbering) were observed in P. absoluta field populations at variable frequencies (Table 2). Among these, the M918T and T929I resistance alleles were shown to be widely distributed in sampled larvae of P. absoluta but completely absent in field populations of P. operculella (Table 2). The only known VGSC mutation detected in P. operculella was L1014F at quite low frequency, whereas the frequency of L1014F in P. absoluta was 100%.
Lastly, an approximately 750 bp fragment of the RDL sequence encompassing the A301S mutation (D. melanogaster numbering) was examined. This mutation was not found in P. operculella, but in P. absoluta. The frequency of the resistance allele A301S was 27%, 25.5%, and 34.5% in larvae collected at GD, PL, and XS sites, respectively (Table 3).

3.2. Haplotype Identification and Phylogeny

Based on the sequence data from field populations (30 individuals for each population) and the NCBI database, 6 PoAce1 (Figure 1) and 4 PaAce1 (Figure 2) haplotypes were identified. Fifteen synonymous and two non-synonymous mutations were identified in P. operculella (Figure 1A). The sequence of PoAce1-H1 originated from a laboratory population without exposure to any insecticides for decades. Blast searches revealed that haplotype PoAce1-H2 shared an identical sequence with the sequence (access number: MG265690.1) from Syria and all three field populations (Figure 1B,C). The PoAce1-H2 and PoAce1-H5 were found with no mutation in combination with either F158Y or A201S (Figure 1B,C). Similarly, two synonymous and three non-synonymous mutations were observed in P. absoluta (Figure 2A). The PaAce1-H1 was a common haplotype in China, the USA (access number: JARADA010000055.1, JAQAHW010000022.1, SNMR01040031.1), and Brazil (access number: KU985167.1) (Figure 2B). It was assumed to be an ancestral haplotype based on network results (Figure 2C).
For a better understanding, the haplotype identification and evolutionary origin of codons 918 and 1014 in VGSC were analyzed separately. A fragment of 569 bp containing VGSC codon 918 from 90 individuals in P. operculella was grouped into one haplotype. Nine haplotypes were identified with a shorter length (55 bp) of PaVGSC. Three non-synonymous mutations (M918T, L925M, and T929I) and one synonymous mutation in the encoding region were found (Figure 3A). The haplotypes PoVGSC-H1, PoVGSC-H2, and PoVGSC-H4 were observed in all field populations, also partially shared with the USA (access number: JARADA010000037.1, JAQAHW010000009.1, SNMR01042398.1) and Iran (access number: KY767010.1) haplotype (Figure 3B). Based on the phylogenetic tree, the M918T mutation has a different evolutionary origin (Figure 3B). As to codon 1014, a fragment of 340 bp containing VGSC codon 1014 from 90 individuals in P. operculella was grouped into three haplotypes. Three SNPs (single-nucleotide polymorphisms) were observed in the intron region, and only one non-synonymous mutation (L1014F) was identified. Similarly, the length of 414 bp containing the VGSC codon 1014 from 90 individuals in P. absoluta was grouped into one haplotype.
Only one haplotype of the RDL gene was identified from 90 individuals in P. operculella, with a length of 567 bp. In P. absoluta, one non-synonymous mutation in the encoding region and two polymorphic sites in the intron were observed (Figure 4A). From a total of 126 sequences obtained from 90 individuals, 5 haplotypes were identified. Of them, PaRDL-H2 was detected in all populations, including the USA (ARADA010000003.1, JAQAHW010000001.1, SNMR01042068.1), while PaRDL-H4 was uniquely distributed in GD city (Figure 4B,C). The phylogenetic tree and network results indicated that haplotypes with the A301S mutation were clustered into different branches (Figure 4B,C), suggesting their multiple evolutionary origins.

3.3. The Divergence Time among Gelechiidae

Based on the phylogenetic timetree, the divergence time of Gelechiidae was estimated due to the absence of fossil records. It has illuminated the temporal aspects of the evolutionary history of P. operculella and P. absoluta, around 46.56 million years ago (MYA) (Figure 5). Remarkably, P. absoluta exhibited an early divergence of approximately 41.27 MYA from P. operculella (36.51 MYA).

3.4. The Collinearity of Specific Genomic Regions in P. operculella and P. absoluta

Three target genes, VGSC (PoScaffold2 versus PaScaffold9), RDL (PoScaffold3 versus PaScaffold1), and Ace1 (PoScaffold113 versus PaScaffold22), were located in different genomic regions. It appeared to be the orthologous genomic regions that span 280~410 kb between P. operculella and P. absoluta with high collinearity. Interestingly, some gene duplications, insertions or deletions, and rearrangements were observed between P. operculella and P. absoluta (Figure 6).

4. Discussion

Research on insecticide resistance has proliferated since the 1960s [49], and investigators worldwide are striving to unravel the mechanisms driving the evolution of insecticide resistance. The most rapid and convenient way to detect the resistance levels in field populations is to genotype field-collected samples for known target-site mutations against commonly used classes of chemistry, thus allowing for the mapping of the distribution and allele frequencies of resistant haplotypes. P. operculella and P. absoluta are two major moth pests that feed on potato plants and particularly P. absoluta has been shown to evolve target-site resistance against several classes of insecticides [37,40,41], whereas such knowledge for P. operculella is largely lacking. Phylogenetic reconstructions suggested that P. absoluta and P. operculella diverged approximately 46.56 MYA (Figure 5). The orthologous genomic regions, including the target gene, between P. operculella and P. absoluta showed high collinearity (Figure 6). However, the target-site mutations exhibited unexpected genetic divergence under insecticide selection pressure.
Mutations in the AChE leading to insensitivity towards OP and CB insecticides are known and largely contribute to OP and CB field resistance in numerous pest insects [50]. Insensitive AChE is conferred by one or more point mutations in the Ace1 gene. Such mutations include, for example, G119S [28], A201S [51], A280T [31], G328A [31], and F331C/Y/W [31,52], to name those most reported. Several of these mutations have been functionally validated in recombinantly expressed wildtype and mutant AChE variants [53]. These results showed a strong correlation between frequencies of the mutation and phenotypic levels of resistance. Among them, A201S has been reported in field-collected P. absoluta populations from Brazil and Europe [41] and is known to confer OP and CB resistance. In our survey, the A201S mutation was also observed in field-collected Chinese P. absoluta samples at high frequency without susceptible homozygous haplotypes (Table 1). Although we did not have bioassay results, the extremely high conservation of the alanine (A201) in all other insect species and previous reports of the same mutation in other OP-resistant insects strongly suggest that this mutation also confers resistance in P. absoluta. Additionally, the F290V mutation was also found in Yunnan and Xinjiang provinces. Several studies functionally validated that these mutations confer insecticide resistance, thus suggesting that the application of OPs and CBs is likely to be ineffective in controlling Chinese populations of P. absoluta. Compared with P. absoluta, only A201S was observed in P. operculella field populations, and even absent in the XW colony, among these common mutations about Ace1 related to insecticide resistance. Du et al. found that P. operculella collected from Yunnan Province developed significant resistance to dichlorvos, methamidophos, and ethyl-parathion in 2006 [6]. In view of the correlation between mutations in Ace1 and OP resistance, the A201S mutation may be one of the reasons for mediating OP resistance in P. operculella in Yunnan Province. Intriguingly, the F158Y mutation was only observed in P. operculella, whereas there is no previous report of the F158Y mutation. Whether the natural F158Y mutation is related to insecticide resistance remains unclear and needs to be clarified in the future.
Pyrethroids, a major class of neurotoxic insecticides, have been extensively used to control various agricultural pests. Among these, cypermethrin is most widely used to control P. operculella in China [20]. Resistance to PYs often arises through knockdown resistance (kdr) resulting from point mutations in the VGSC that reduce the sensitivity of the insect nervous system to these compounds [54]. Three target resistance-conferring VGSC mutations, including kdr (L1014F), kdr-his (L1014H), and super-kdr (M918T + L1014F), have been identified in multiple species, such as Diptera [55], Lepidoptera [56], and Siphonaptera [57]. Based on our survey, the kdr and super-kdr mutations have been commonly detected in P. absoluta field populations. The construction of the haplotype network and phylogeny analysis indicated that the substitution of methionine (M918) to threonine (T918) has different origins (Figure 3). These indications suggest that the P. absoluta in China would exhibit resistance to PYs. Additionally, the T929I mutation was observed with relatively high frequencies in P. absoluta field populations. The T929I mutation was first reported in the pyrethroid-resistant diamondback moth, P. xylostella, and is presumed to play an important role in pyrethroid resistance [58]. T929I in combination with L1014F was found to make the sodium channel highly insensitive to a range of type I and type II pyrethroid insecticides and DDT [59]. The combination of M918T, T929I, and L1014F mutations was observed in P. absoluta field populations from Yunnan and Xinjiang provinces and reported in populations from Europe and Brazil [40]. Generally, the combined mutations often result in higher resistance levels, affecting a broad spectrum of pyrethroids. Our findings suggested that P. absoluta has historically undergone strong selection with a range of pyrethroid insecticides, resulting in significant control problems under applied conditions in the field. In P. operculella, the kdr mutation was at low frequency in field populations, with other mutations associated with PY resistance barely present. Long-term use of chemical insecticides over decades has led to ineffective pest control, resulting in widespread distribution throughout China. It is inferred that P450s may play an important role in insecticide resistance in P. operculella field populations in China.
GABA receptors, targeted by cyclodienes and fiproles, are ligand-gated homopentameric chloride channels consisting of RDL subunits encoded by the RDL gene. One common target-site mutation is A301S, conferring high levels of cyclodiene resistance and lower levels of fiprole resistance, e.g., in Nilaparvata lugens [34]. Here, we detected the A301S mutation at moderate frequency in P. absoluta field populations, but not in P. operculella field populations. The A301S mutation in three P. absoluta field populations was narrowly distributed at a low frequency, and its contribution to insecticide failures in tomato leafminer control remains unclear.
Overall, we detected at varying frequencies multiple target-site mutations in P. operculella and P. absoluta field populations in China likely to compromise the efficacy of commonly used chemical classes of insecticides such as organophosphates and pyrethroids. Nevertheless, P. operculella showed a much lower frequency of target-site mutations than P. absoluta, suggesting that the latter is more difficult to control with less costly or less expensive insecticides such as pyrethroids. However, we have not investigated the potential contribution of upregulated metabolic detoxification mechanisms which can modulate and contribute substantially to target-site mediated resistance to pyrethroids [60]. The current control strategy in most areas sampled in China involves spraying insecticides to protect the potato crops, but some of the recommended chemical classes such as pyrethroids and OPs are likely to fail under field conditions, particularly against P. absoluta. Therefore, resistance monitoring, such as that shown here by genotyping samples for the presence of target-site mutations, is crucial for the implementation of rational strategies of insecticide application and resistance management.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/insects15030194/s1, Table S1: Primers informations; Table S2: Sequence informations; Table S3: Informations on field strain of P. operculella and P. absoluta used in the present study.

Author Contributions

Conceptualization, Y.G., S.L. and J.Z.; methodology, J.Z., Y.G. and S.L.; software, J.Z., J.L. (Jiaxin Liang) and R.C.; validation, J.Z., J.L. (Juan Liu) and W.L.; formal analysis, Y.G., S.L., R.N. and J.Z.; investigation, J.Z. and J.L. (Jiaxin Liang); resources, J.Z., R.C., J.L. (Juan Liu) and W.L.; data curation, J.Z.; writing—original draft preparation, J.Z.; writing—review & editing, Y.G., S.L. and R.N.; visualization, J.Z.; supervision, Y.G. and S.L.; project administration, Y.G.; funding acquisition, Y.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Guangdong Major Project of Basic and Applied Basic Research (2021B0301030004), the National Natural Science Foundation of China (grants 32372535), and the Science and Technology Innovation Project of the Chinese Academy of Agricultural Sciences.

Data Availability Statement

The data supporting the conclusions of this article are provided within the article. The original datasets analyzed in this study are available from the corresponding author upon request. The nucleotide sequences reported in this study have been deposited in the GenBank database listed in Table S2.

Acknowledgments

We would like to thank Yanli Yang (Yunnan Agricultural University) for providing the locations and sample collection.

Conflicts of Interest

The author (Ralf Nauen) was employed by Bayer, a manufacturer of crop protection products. The remaining 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.

References

  1. Si, H.J.; Zhang, N.; Tang, X.; Yang, J.W.; Wen, Y.K.; Wang, L.; Zhou, X.Y. Transgenic Research in Tuber and Root Crops: A Review. In Genetic Engineering of Horticultural Crops; Rout, G.R., Peter, K.V., Eds.; Academic Press: Cambridge, MA, USA, 2018; pp. 225–248. [Google Scholar] [CrossRef]
  2. Wang, Q.B.; Zhang, W. China’s potato industry and potential impacts on the global market. Am. J. Potato Res. 2004, 81, 101–109. [Google Scholar] [CrossRef]
  3. Amiri, A.N.; Bakhsh, A. An effective pest management approach in potato to combat insect pests and herbicide. 3 Biotech. 2019, 9, 16. [Google Scholar] [CrossRef] [PubMed]
  4. Rondon, S.I. The potato tuberworm: A literature review of its biology, ecology, and control. Am. J. Potato Res. 2010, 87, 149–166. [Google Scholar] [CrossRef]
  5. Rondon, S.I.; Gao, Y.L. Potato tuberworm. In Insect Pests of Potato: Global Perspectives on Biology and Management; Alyokhin, A., Rondon, S.I., Gao, Y.L., Eds.; Academic Press: Oxford, UK, 2022; pp. 149–160. Available online: https://www.sciencedirect.com/book/9780123868954/insect-pests-of-potato (accessed on 2 February 2024).
  6. Du, L.T.; Li, Z.Y.; Zhou, L.M.; Li, C.M.; Tang, X. Comparsion of five chemicals in three groups for their effects controlling larva of potato moth. China Potato J. 2006, 20, 92–93. [Google Scholar] [CrossRef]
  7. Gao, Y.L. Green control techniques for potato tuberworm (Phthorimaea operculella). Sci. Agric. Sin. 2021, 54, 533–535. [Google Scholar] [CrossRef]
  8. Yan, J.J.; Guo, W.C.; Li, G.Q.; Pan, H.P.; Chen, B.; Zhou, W.W.; Li, S.H.; Gao, Y.L. Current status and prospects of the management of important insect pests on potato in China. Plant Protect. 2023, 49, 190–195. [Google Scholar] [CrossRef]
  9. Yang, J.; Wang, F.G.; Gao, S.H.; Ma, H.M. Identification and integrated prevention and control of tomato leaf miner moth, Tuta absoluta. Appl. Tech. Rural Area 2021, 3, 88–89. [Google Scholar]
  10. Desneux, N.; Wajnberg, E.; Wyckhuys, K.A.G.; Burgio, G.; Arpaia, S.; Narváez-Vasquez, C.A.; Gonzalez-Cabrera, J.; Ruescas, D.C.; Tabone, E.; Frandon, J.; et al. Biological invasion of European tomato crops by Tuta absoluta: Ecology, geographic expansion and prospects for biological control. J. Pest. Sci. 2010, 83, 197–215. [Google Scholar] [CrossRef]
  11. Desneux, N.; Luna, M.G.; Guillemaud, T.; Urbaneja, A. The invasive South American tomato pinworm, Tuta absoluta, continues to spread in Afro-Eurasia and beyond: The new threat to tomato world production. J. Pest Sci. 2011, 84, 403–408. [Google Scholar] [CrossRef]
  12. Brévault, T.; Bernadas, G.; Sylla, S.; Diatte, M.; Diarra, K. Tuta absoluta Meyrick (Lepidoptera: Gelechiidae): A new threat to tomato production in sub-Saharan Africa. Afr. Entomol. 2014, 22, 441–444. [Google Scholar] [CrossRef]
  13. Zhang, G.F.; Ma, D.Y.; Wang, Y.S.; Gao, Y.H.; Liu, W.X.; Zhang, R.; Fu, W.J.; Xian, X.Q.; Wang, J.; Kuang, M.; et al. First report of the South American tomato leafminer, Tuta absoluta (Meyrick), in China. J. Integr. Agr. 2020, 19, 1912–1917. [Google Scholar] [CrossRef]
  14. Desneux, N.; Han, P.; Mansour, R.; Arnó, J.; Brévault, T.; Campos, M.R.; Chailleux, A.; Guedes, R.N.C.; Karimi, J.; Konan, K.A.J.; et al. Integrated pest management of Tuta absoluta: Practical implementations across different world regions. J. Pest Sci. 2022, 95, 17–39. [Google Scholar] [CrossRef]
  15. Zhang, G.F.; Ma, D.Y.; Liu, W.X.; Wang, Y.S.; Fu, W.J.; Wang, J.; Gao, Y.H.; Wan, F.H. The arrival of Tuta absoluta (Meyrick) (Lepidoptera: Gelechiidae), in China. J. Biosaf. 2019, 28, 200–203. [Google Scholar] [CrossRef]
  16. Lu, Y.Y. Invasion alert: Rapid spread of the tomato leafminer, Tuta absoluta (Lepidoptera: Gelechiidae), an invasive pest in China. J. Environ. Entomol. 2021, 43, 526–528. [Google Scholar] [CrossRef]
  17. Zhang, G.F.; Zhang, Y.B.; Xian, X.Q.; Liu, W.X.; Li, P.; Liu, W.C.; Liu, H.; Fen, X.D.; Lv, Z.C.; Wang, Y.S.; et al. Damage of an important and newly invaded agricultural pest, Phthorimaea absoluta, and its prevention and management measures. Plant Protect. 2022, 48, 51–58. [Google Scholar]
  18. Zhang, M.D.; Yan, J.J.; Ali, A.; Gao, Y.L. Potato plant variety affects the performance and oviposition preference of Phthorimaea operculella Zeller (Lepidoptera: Gelechiidae). Pest Manag. Sci. 2022, 78, 3912–3919. [Google Scholar] [CrossRef] [PubMed]
  19. Yan, J.J.; Zhang, M.D.; Ali, A.; Du, X.; Mei, X.D.; Gao, Y.L. Optimization and field evaluation of sex-pheromone of potato tuber moth, Phthorimaea operculella Zeller (Lepidoptera: Gelechiidae). Pest Manag. Sci. 2022, 78, 3903–3911. [Google Scholar] [CrossRef]
  20. Yan, J.J.; Zhang, M.D.; Gao, Y.L. Biology, ecology and integrated management of potato tuber moth, Phthorimaea operculella (Lepidoptera: Gelechiidae). Acta Entomol. Sin. 2019, 62, 1469–1482. [Google Scholar] [CrossRef]
  21. Raymond, M.; Berticat, C.; Weill, M.; Pasteur, N.; Chevillon, C. Insecticide resistance in the mosquito culex pipiens: What have we learned about adaptation? Genetica 2001, 8, 287–296. [Google Scholar] [CrossRef]
  22. Feyereisen, R.; Dermauw, W.; Van-Leeuwen, T. Genotype to phenotype, the molecular and physiological dimensions of resistance in arthropods. Pestic Biochem. Physiol. 2015, 121, 61–77. [Google Scholar] [CrossRef]
  23. Mugenzi, L.M.J.; Menze, B.D.; Tchouakui, M.; Wondji, M.J.; Irving, H.; Tchoupo, M.; Hearn, J.; Weedall, G.D.; Riveron, J.M.; Wondji, C.S. Cis-regulatory CYP6P9b P450 variants associated with loss of insecticide-treated bed net efficacy against Anopheles funestus. Nat. Commun. 2019, 10, 4652. [Google Scholar] [CrossRef] [PubMed]
  24. Nauen, R.; Bass, C.; Feyereisen, R.; Vontas, J. The role of cytochrome P450s in insect toxicology and resistance. Annu. Rev. Entomol. 2022, 67, 105–124. [Google Scholar] [CrossRef] [PubMed]
  25. Dermauw, W.; Van-Leeuwen, T.; Feyereisen, R. Diversity and evolution of the P450 family in arthropods. Insect Biochem. Mol. Biol. 2020, 127, 103490. [Google Scholar] [CrossRef] [PubMed]
  26. Feyereisen, R. Insect CYP genes and P450 enzymes. In Insect Molecular Biology and Biochemistry; Gilbert, L., Ed.; Elsevier: Amsterdam, The Netherlands, 2012; pp. 236–316. [Google Scholar] [CrossRef]
  27. 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] [PubMed]
  28. Engdahl, C.; Knutsson, S.; Fredriksson, S.A.; Linusson, A.; Bucht, G.; Ekström, F. Acetylcholinesterases from the disease vectors Aedes aegypti and Anopheles gambiae: Functional characterization and comparisons with vertebrate orthologues. PLoS ONE 2015, 10, e0138598. [Google Scholar] [CrossRef] [PubMed]
  29. Andrews, M.C.; Callaghan, A.; Field, L.M.; Williamson, M.S.; Moores, G.D. Identification of mutations conferring insecticide-insensitive AChE in the cotton-melon aphid, Aphis gossypii Glover. Insect Mol. Biol. 2004, 13, 555–561. [Google Scholar] [CrossRef]
  30. Chang, C.; Cheng, X.; Huang, X.Y.; Dai, S.M. Amino acid substitutions of acetylcholinesterase associated with carbofuran resistance in Chilo suppressalis. Pest Manag. Sci. 2014, 70, 1930–1935. [Google Scholar] [CrossRef]
  31. Ilias, A.; Vontas, J.; Tsagkarakou, A. Global distribution and origin of target site insecticide resistance mutations in Tetranychus urticae. Insect Biochem. Mol. Biol. 2014, 48, 17–28. [Google Scholar] [CrossRef]
  32. Saavedra-Rodriguez, K.; Campbell, C.L.; Lozano, S.; Penilla-Navarro, P.; Lopez-Solis, A.; Solis-Santoyo, F. Rodriguez, A.D.; Perera, R.; Black Lv, W.C. Permethrin resistance in Aedes aegypti: Genomic variants that confer knockdown resistance, recovery, and death. PLoS Genet. 2021, 17, e1009606. [Google Scholar] [CrossRef]
  33. Zhao, Y.H.; Feng, X.Y.; Li, M.; Qiu, X.H. The double-mutation (M918I + L1014F) kdr allele is fixed in Cimex hemipterus populations in Guangxi, China. Bull Entomol. Res. 2020, 110, 506–511. [Google Scholar] [CrossRef]
  34. Garrood, W.T.; Zimmer, C.T.; Gutbrod, O.; Lüke, B.; Williamson, M.S.; Bass, C.; Nauen, R.; Davies, T.G.E. Influence of the RDL A301S mutation in the brown planthopper Nilaparvata lugens on the activity of phenylpyrazole insecticides. Pestic Biochem. Physiol. 2017, 142, 1–8. [Google Scholar] [CrossRef]
  35. Doğramacı, M.; Tingey, W.M. Comparison of insecticide resistance in a North American field population and a laboratory colony of potato tuberworm (Lepidoptera: Gelechiidae). J. Pest Sci. 2008, 81, 17–22. [Google Scholar] [CrossRef]
  36. El-Kady, H. Insecticide resistance in potato tuber moth Phthorimaea operculella Zeller in Egypt. J. Am. Sci. 2011, 7, 263–266. Available online: http://www.americanscience.org (accessed on 2 February 2024).
  37. Roditakis, E.; Steinbach, D.; Moritz, G.; Vasakis, E.; Stavrakaki, M.; Ilias, A.; García-Vidal, L.; Aguirre, M.R.M.; Bielza, P.; Morou, E.; et al. Ryanodine receptor point mutations confer diamide insecticide resistance in tomato leafminer, Tuta absoluta (Lepidoptera: Gelechiidae). Insect Biochem. Mol. Biol. 2017, 80, 11–20. [Google Scholar] [CrossRef]
  38. Campos, M.R.; Silva, T.B.; Silva, W.M.; Silva, J.E.; Siqueira, H.A. Spinosyn resistance in the tomato borer Tuta absoluta (Meyrick) (Lepidoptera: Gelechiidae). J. Pest Sci. 2015, 88, 405–412. [Google Scholar] [CrossRef]
  39. Silva, J.E.; Ribeiro, L.M.D.S.; Vinasco, N.; Guedes, R.N.C.; Siqueira, H.Á.A. Field-evolved resistance to chlorantraniliprole in the tomato pinworm Tuta absoluta: Inheritance, cross-resistance profile, and metabolism. J. Pest Sci. 2019, 92, 1421–1431. [Google Scholar] [CrossRef]
  40. Haddi, K.; Berger, M.; Bielza, P.; Cifuentes, D.; Field, L.M.; Gorman, K.; Rapisarda, C.; Williamson, M.S.; Bass, C. Identification of mutations associated with pyrethroid resistance in the voltage-gated sodium channel of the tomato leaf miner (Tuta absoluta). Insect Biochem. Mol. Biol. 2012, 42, 506–513. [Google Scholar] [CrossRef] [PubMed]
  41. Haddi, K.; Berger, M.; Bielza, P.; Rapisarda, C.; Williamson, M.S.; Moores, G.; Bass, C. Mutation in the ace-1 gene of the tomato leaf miner (Tuta absoluta) associated with organophosphates resistance. J. Appl. Entomol. 2017, 141, 612–619. [Google Scholar] [CrossRef]
  42. Kwon, M.; Kim, J.; Maharjan, R.; Choi, J.Y.; Kim, G.H. Change in the distribution of the potato tuber moth, Phthorimaea operculella (Zeller) (Lepidoptera: Gelechiidae), in Korea. J. Asia-Pac. Entomol. 2017, 20, 1249–1253. [Google Scholar] [CrossRef]
  43. Rinkevich, F.D.; Zhang, L.; Hamm, R.L.; Brady, S.G.; Lazzaro, B.P.; Scott, J.G. Frequencies of the pyrethroid resistance alleles of Vssc1 and CYP6D1 in house flies from the eastern United States. Insect Mol. Biol. 2006, 15, 157–167. [Google Scholar] [CrossRef] [PubMed]
  44. Edgar, R.C. MUSCLE: A multiple sequence alignment method with reduced time and space complexity. BMC Bioinform. 2004, 5, 113. [Google Scholar] [CrossRef] [PubMed]
  45. Librado, P.; Rozas, J. DnaSP v5: A software for comprehensive analysis of DNA polymorphism data. Bioinformatics 2009, 25, 1451–1452. [Google Scholar] [CrossRef]
  46. Paradis, E. Analysis of haplotype networks: The randomized minimum spanning tree method. Methods Ecol. Evol. 2018, 9, 1308–1317. [Google Scholar] [CrossRef]
  47. Kumar, S.; Stecher, G.; Li, M.; Knyaz, C.; Tamura, K. MEGA X: Molecular Evolutionary Genetics Analysis across Computing Platforms. Mol. Biol. Evol. 2018, 35, 1547–1549. [Google Scholar] [CrossRef]
  48. Chen, C.; Chen, H.; Zhang, Y.; Thomas, H.R.; Frank, M.H.; He, Y.; Xia, R. TBtools: An integrative toolkit developed for interactive analyses of big biological data. Mol. Plant. 2020, 13, 1194–1202. [Google Scholar] [CrossRef]
  49. Bielza, P. Insecticide resistance in natural enemies. In Advances in Insect Control and Resistance Management; Horowitz, A.R., Ishaaya, I., Eds.; Springer: Berlin/Heidelberg, Germany, 2016; pp. 313–329. [Google Scholar] [CrossRef]
  50. Lee, S.H.; Kim, Y.H.; Kwon, D.H.; Cha, D.J.; Kim, J.H. Mutation and duplication of arthropod acetylcholinesterase: Implications for pesticide resistance and tolerance. Pestic Biochem. Physiol. 2015, 120, 118–124. [Google Scholar] [CrossRef] [PubMed]
  51. Zhen, C.; Miao, L.; Liang, P.; Gao, X.W. Survey of organophosphate resistance and an Ala216Ser substitution of acetylcholinesterase-1 gene associated with chlorpyrifos resistance in Apolygus lucorum (Meyer-Dür) collected from the transgenic Bt cotton fields in China. Pestic Biochem. Physiol. 2016, 132, 29–37. [Google Scholar] [CrossRef] [PubMed]
  52. Khajehali, J.; Van-Leeuwen, T.; Grispou, M.; Morou, E.; Alout, H.; Weill, M.; Tirry, L.; Vontas, J.; Tsagkarakou, A. Acetylcholinesterase point mutations in European strains of Tetranychus urticae (Acari: Tetranychidae) resistant to organophosphates. Pest Manag. Sci. 2010, 66, 220–228. [Google Scholar] [CrossRef] [PubMed]
  53. Kim, S.; Yoon, K.A.; Lee, S.H. Short communication: Screening of proof-of-concept mutations of honey bee acetylcholinesterase 2 conferring resistance to organophosphorus and carbamate insecticides. Comp. Biochem. Physiol. C Toxicol. Pharmacol. 2023, 264, 109524. [Google Scholar] [CrossRef]
  54. Dong, K.; Du, Y.; Rinkevich, F.; Nomura, Y.; Xu, P.; Wang, L.X.; Sliver, K.; Zhorov, B. Molecular biology of insect sodium channels and pyrethroid resistance. Insect Biochem. Mol. Biol. 2014, 50, 1–17. [Google Scholar] [CrossRef] [PubMed]
  55. Lynd, A.; Oruni, A.; Van’t, H.A.E.; Morgan, J.C.; Naego, L.B.; Pipini, D.; O’Kines, K.A.; Bobanga, T.L.; Donnelly, M.J.; Weetman, D. Insecticide resistance in Anopheles gambiae from the northern Democratic Republic of Congo, with extreme knockdown resistance (kdr) mutation frequencies revealed by a new diagnostic assay. Malar. J. 2018, 17, 412. [Google Scholar] [CrossRef]
  56. Tian, X.R.; Liu, J.D.; Guo, Z.M.; Hu, B.; Kibe, M.D.; Wang, S.Y.; Wei, Q.; Su, J.Y. The characteristics of voltage-gated sodium channel and the association with lambda cyhalothrin resistance in Spodoptera exigua. J. Asia-Pac. Entomol. 2018, 21, 1020–1027. [Google Scholar] [CrossRef]
  57. Ghavami, M.B.; Haghi, F.P.; Alibabaei, Z.; Enayati, A.A.; Vatandoost, H. First report of target site insensitivity to pyrethroids in human flea, Pulex irritans (Siphonaptera: Pulicidae). Pestic Biochem. Physiol. 2018, 146, 97–105. [Google Scholar] [CrossRef] [PubMed]
  58. Schuler, T.H.; Martinez-Torres, D.; Thompson, A.J.; Denholm, I.; Devonshire, A.L.; Duce, I.R.; Williamson, M.S. Toxicological, electrophysiological, and molecular characterisation of knockdown resistance to pyrethroid insecticides in the diamondback moth, Plutella xylostella (L.). Pestic Biochem Phys. 1998, 59, 169–182. [Google Scholar] [CrossRef]
  59. Usherwood, P.N.; Davies, T.G.; Mellor, I.R.; O’Reilly, A.O.; Peng, F.; Vais, H.; Khambay, B.P.S.; Field, L.M.; Williamson, M.S. Mutations in DIIS5 and the DIIS4-S5 linker of Drosophila melanogaster sodium channel define binding domains for pyrethroids and DDT. FEBS Lett. 2007, 581, 5485–5492. [Google Scholar] [CrossRef]
  60. Samantsidis, G.R.; Panteleri, R.; Denecke, S.; Kounadi, S.; Christou, I.; Nauen, R.; Douris, V.; Vontas, J. ‘What I cannot create, I do not understand’: Functionally validated synergism of metabolic and target site insecticide resistance. Proc. Biol. Sci. 2020, 287, 20200838. [Google Scholar] [CrossRef]
Insects 15 00194 g001
Figure 1. The evolution of the Ace1 haplotypes in different populations of P. operculella. (A) The nucleotide alignment of six haplotypes of Ace1, H1–H6, is presented. Mutations are labeled in red dotted boxes. (B) The phylogenetic tree is inferred using the maximum likelihood method. Bootstrap values are conducted for 1000 replicates and only scores above 50 are shown. The abbreviations of locations are LL (Luliang), XW (Xuanwei), ML (Malong), LS (laboratory population) from China, and CD (Damascus) from Syria. (C) The network presents the genealogy of Ace1 haplotypes, and the short solid lines indicate the possible mutational steps between haplotypes.
Figure 1. The evolution of the Ace1 haplotypes in different populations of P. operculella. (A) The nucleotide alignment of six haplotypes of Ace1, H1–H6, is presented. Mutations are labeled in red dotted boxes. (B) The phylogenetic tree is inferred using the maximum likelihood method. Bootstrap values are conducted for 1000 replicates and only scores above 50 are shown. The abbreviations of locations are LL (Luliang), XW (Xuanwei), ML (Malong), LS (laboratory population) from China, and CD (Damascus) from Syria. (C) The network presents the genealogy of Ace1 haplotypes, and the short solid lines indicate the possible mutational steps between haplotypes.
Insects 15 00194 g001
Insects 15 00194 g002
Figure 2. The evolution of the Ace1 haplotypes in different populations of P. absoluta. (A) The nucleotide alignment of four haplotypes of Ace1, H1–H4, is presented. Mutations are labeled in red dotted boxes. (B) The phylogenetic tree is inferred using the maximum likelihood method and presented. Bootstrap values are conducted for 1000 replicates and only scores above 50 are shown. The abbreviations of locations are PL (Panlong), GD (Guandu), and XS (Xinshi) from China, FL (Florida) and CA (California) from the USA, and MG (Minas Gerais) from Brazil. (C) The network presents the genealogy of Ace1 haplotypes, and the short solid lines indicate the possible mutational steps between haplotypes.
Figure 2. The evolution of the Ace1 haplotypes in different populations of P. absoluta. (A) The nucleotide alignment of four haplotypes of Ace1, H1–H4, is presented. Mutations are labeled in red dotted boxes. (B) The phylogenetic tree is inferred using the maximum likelihood method and presented. Bootstrap values are conducted for 1000 replicates and only scores above 50 are shown. The abbreviations of locations are PL (Panlong), GD (Guandu), and XS (Xinshi) from China, FL (Florida) and CA (California) from the USA, and MG (Minas Gerais) from Brazil. (C) The network presents the genealogy of Ace1 haplotypes, and the short solid lines indicate the possible mutational steps between haplotypes.
Insects 15 00194 g002
Insects 15 00194 g003
Figure 3. The evolution of the VGSC haplotypes in different populations of P. absoluta. (A) The nucleotide alignment of nine haplotypes of VGSC containing codon 918, H1–H9, is presented. Mutations are labeled in red dotted boxes. (B) The phylogenetic tree is inferred using the maximum likelihood method. Bootstrap values are conducted for 1000 replicates and only scores above 50 are shown. The abbreviations of locations are PL (Panlong), GD (Guandu), and XS (Xinshi) from China, FL (Florida) and CA (California) from the USA, and IRN (Iran) from Iran. (C) The network presents the genealogy of VGSC haplotypes, and the short solid lines indicate the possible mutational steps between haplotypes.
Figure 3. The evolution of the VGSC haplotypes in different populations of P. absoluta. (A) The nucleotide alignment of nine haplotypes of VGSC containing codon 918, H1–H9, is presented. Mutations are labeled in red dotted boxes. (B) The phylogenetic tree is inferred using the maximum likelihood method. Bootstrap values are conducted for 1000 replicates and only scores above 50 are shown. The abbreviations of locations are PL (Panlong), GD (Guandu), and XS (Xinshi) from China, FL (Florida) and CA (California) from the USA, and IRN (Iran) from Iran. (C) The network presents the genealogy of VGSC haplotypes, and the short solid lines indicate the possible mutational steps between haplotypes.
Insects 15 00194 g003
Insects 15 00194 g004
Figure 4. The evolution of the RDL haplotypes in different populations of P. absoluta. (A) The nucleotide alignment of five haplotypes of RDL, H1–H5, is presented. Mutations are labeled in red dotted boxes. The intron region is marked in a blue vertical line, and the GT-AG is presented in italics. (B) The phylogenetic tree is inferred using the maximum likelihood method. Bootstrap values are conducted for 1000 replicates and only scores above 50 are shown. The abbreviations of locations are PL (Panlong), GD (Guandu), and XS (Xinshi) from China, and FL (Florida) and CA (California) from the USA. (C) The network presents the genealogy of RDL haplotypes, and the short solid lines indicate the possible mutational steps between haplotypes.
Figure 4. The evolution of the RDL haplotypes in different populations of P. absoluta. (A) The nucleotide alignment of five haplotypes of RDL, H1–H5, is presented. Mutations are labeled in red dotted boxes. The intron region is marked in a blue vertical line, and the GT-AG is presented in italics. (B) The phylogenetic tree is inferred using the maximum likelihood method. Bootstrap values are conducted for 1000 replicates and only scores above 50 are shown. The abbreviations of locations are PL (Panlong), GD (Guandu), and XS (Xinshi) from China, and FL (Florida) and CA (California) from the USA. (C) The network presents the genealogy of RDL haplotypes, and the short solid lines indicate the possible mutational steps between haplotypes.
Insects 15 00194 g004
Insects 15 00194 g005
Figure 5. Timetree analysis for Gelechiidae members using the RelTime method. Drosophila melanogaster was used as an outgroup. Helicoverpa armigera, Spodoptera frugiperda, and Plutella xylostella were set as calibration for Gelechiidae members lacking fossil histories. The divergence time of P. operculella and P. absoluta was marked in red circle. The red font represents the two target species of this study, The gray font represents the use of Drosophila mela-nogaster as an outgroup.
Figure 5. Timetree analysis for Gelechiidae members using the RelTime method. Drosophila melanogaster was used as an outgroup. Helicoverpa armigera, Spodoptera frugiperda, and Plutella xylostella were set as calibration for Gelechiidae members lacking fossil histories. The divergence time of P. operculella and P. absoluta was marked in red circle. The red font represents the two target species of this study, The gray font represents the use of Drosophila mela-nogaster as an outgroup.
Insects 15 00194 g005
Insects 15 00194 g006
Figure 6. The collinearity map depicting synteny of target gene loci in P. operculella and P. absoluta. The complete genomes of P. operculella and P. absoluta are available under the GenBank assembly numbers GCA_024500475.1 and GCA_027580185.1, respectively. The genomic region, 280~410 kb, analyzed included Ace1, VGSC, and RDL. The obtained synteny between P. operculella (marked in marine blue) and P. absoluta (marked in tiffany blue) is shown in (A), (B), and (C), respectively. The arrow indicates the direction of protein coding. The synteny of target genes (Ace1, VGSC, and RDL) is marked in a dashed line.
Figure 6. The collinearity map depicting synteny of target gene loci in P. operculella and P. absoluta. The complete genomes of P. operculella and P. absoluta are available under the GenBank assembly numbers GCA_024500475.1 and GCA_027580185.1, respectively. The genomic region, 280~410 kb, analyzed included Ace1, VGSC, and RDL. The obtained synteny between P. operculella (marked in marine blue) and P. absoluta (marked in tiffany blue) is shown in (A), (B), and (C), respectively. The arrow indicates the direction of protein coding. The synteny of target genes (Ace1, VGSC, and RDL) is marked in a dashed line.
Insects 15 00194 g006
Table 1. Frequencies of Ace1 genotypes in different field populations of P. operculella and P. absoluta.
Table 1. Frequencies of Ace1 genotypes in different field populations of P. operculella and P. absoluta.
Species Populations n F158Y A201S L231V F290V
FF/YYAA/SSLL/VVFF/VV
P. operculellaLL3083107871300010010000
ML30867790730010010000
XW3090100100000010010000
P. absolutaPL3010000077237723040600
GD3010000067336337073270
XS301000006337973070300
Note: The abbreviations of locations are LL (Luliang), ML (Malong), XW (Xuanwei), PL (Panlong), GD (Guandu) and XS (Xinshi) of China. The abbreviation of amino acids are F (phenylalanine), Y (tyrosine), A (alanine), S (serine), L (leucine) and V (valine).
Table 2. Frequencies of VGSC genotypes in different field populations of P. operculella and P. absoluta.
Table 2. Frequencies of VGSC genotypes in different field populations of P. operculella and P. absoluta.
Species Populations n M918T L925M T929I L1014F
MM/TTLL/MMTT/IILL/FF
P. operculellaLL3010000100001000084133
ML3010000100001000087103
XW3010000100001000077203
P. absolutaPL30304723937023473000100
GD30534079010010474300100
XS30503713937013404700100
Note: The abbreviations of locations are LL (Luliang), ML (Malong), XW (Xuanwei), PL (Panlong), GD (Guandu) and XS (Xinshi) of China. The abbreviations of amino acids are M (methionine), T (threonine), L (leucine), F (phenylalanine) and I (isoleucine).
Table 3. Frequencies of RDL genotypes in different field populations of P. operculella and P. absoluta.
Table 3. Frequencies of RDL genotypes in different field populations of P. operculella and P. absoluta.
Species Populations n A301S
A/AA/SS/S
P. operculellaLL3010000
ML3010000
XW3010000
P. absolutaPL30632017
GD30661717
XS30542323
Note: The abbreviations of locations are LL (Luliang), ML (Malong), XW (Xuanwei), PL (Panlong), GD (Guandu) and XS (Xinshi) of China. The abbreviations of amino acids are A (alanine) and S (serine).
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhu, J.; Chen, R.; Liu, J.; Lin, W.; Liang, J.; Nauen, R.; Li, S.; Gao, Y. Presence of Multiple Genetic Mutations Related to Insecticide Resistance in Chinese Field Samples of Two Phthorimaea Pest Species. Insects 2024, 15, 194. https://doi.org/10.3390/insects15030194

AMA Style

Zhu J, Chen R, Liu J, Lin W, Liang J, Nauen R, Li S, Gao Y. Presence of Multiple Genetic Mutations Related to Insecticide Resistance in Chinese Field Samples of Two Phthorimaea Pest Species. Insects. 2024; 15(3):194. https://doi.org/10.3390/insects15030194

Chicago/Turabian Style

Zhu, Jiang, Ruipeng Chen, Juan Liu, Weichao Lin, Jiaxin Liang, Ralf Nauen, Suhua Li, and Yulin Gao. 2024. "Presence of Multiple Genetic Mutations Related to Insecticide Resistance in Chinese Field Samples of Two Phthorimaea Pest Species" Insects 15, no. 3: 194. https://doi.org/10.3390/insects15030194

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