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

Pest Status, Bio-Ecology, and Area-Wide Management of Mirids in East Asia

Mirids (Hemiptera: Heteroptera: Miridae) feed upon a wide variety of cultivated and wild plants and can be economically important crop pests. They have traditionally been perceived as innocuous herbivores in East Asia; however, population levels of various mirid species have dramatically increased over the past decades. High-profile pests such as Apolygus spp., Adelphocoris spp., and Lygus spp. are now widely distributed across the region, and their infestation pressure is associated with climate, agroecological conditions, and farming practices. This review outlines how an in-depth understanding of pest biology, a systems-level characterization of pest ecology, and a comprehensive evaluation of integrated pest management tactics have enabled sustainable management of mirids across crop boundaries and harvest cycles. This work underscores how more holistic, integrative research approaches can accelerate the implementation of area-wide management of generalist pests, effectively prevent pest population build-up and yield impact, and shrink the environmental footprint of agriculture. In addition to highlighting the merits of interdisciplinary systems approaches, we discuss prospects and challenges for the sustainable management of polyphagous mirid pests in landscape matrices.

The chromosome-scale reference genome of mirid bugs (Adelphocoris suturalis) genome provides insights into omnivory, insecticide resistance, and survival adaptation

BACKGROUND

Adelphocoris suturalis (Hemiptera: Miridae) is a notorious agricultural pest, which causes serious economic losses to a diverse range of agricultural crops around the world. The poor understanding of its genomic characteristics has seriously hindered the establishment of sustainable and environment-friendly agricultural pest management through biotechnology and biological insecticides.

RESULTS

Here, we report a chromosome-level assembled genome of A. suturalis by integrating Illumina short reads, PacBio, 10x Chromium, and Hi-C mapping technologies. The resulting 1.29 Gb assembly contains twelve chromosomal pseudomolecules with an N50 of 1.4 and 120.6 Mb for the contigs and scaffolds, respectively, and carries 20,010 protein-coding genes. The considerable size of the A. suturalis genome is predominantly attributed to a high amount of retrotransposons, especially long interspersed nuclear elements (LINEs). Transcriptomic and phylogenetic analyses suggest that A. suturalis-specific candidate effectors, and expansion and expression of gene families associated with omnivory, insecticide resistance and reproductive characteristics, such as digestion, detoxification, chemosensory receptors and long-distance migration likely contribute to its strong environmental adaptability and ability to damage crops. Additionally, 19 highly credible effector candidates were identified and transiently overexpressed in Nicotiana benthamiana for functional assays and potential targeting for insect resistance genetic engineering.

CONCLUSIONS

The high-quality genome of A. suturalis provides an important genomic landscape for further investigations into the mechanisms of omnivory, insecticide resistance and survival adaptation, and for the development of integrated management strategies.

The role of tetradecane in the identification of host plants by the mirid bugs Apolygus lucorum and Adelphocoris suturalis and potential application in pest management

The mirid bugs Apolygus lucorum and Adelphocoris suturalis are considered serious pests of many crops in China, the host plant recognition of these pests remains unclear. The current study investigated the vital odor cues of two mirid bugs and evaluated the role of olfactory recognition in host recognition. The GC-EAD response of mirid bugs to volatiles of their host plant Phaseolus vulgaris was tested. Tetradecane, 2-propyl-1-pentanol, and dodecanal elicited strong EAG responses by mirid bugs and were tested with field experiments. The results indicated tetradecane was significantly more attractive than other attractants, yielding 30.33 ± 2.19 mirid bugs trapped during 7 days. The selected response rates to tetradecane were above 60%, which was most attractive to female A. lucorum at 1.5 mg/ml. Among seven tetradecane derivatives, tetradecane and tetradecanoic acid were the most potent attractants to A. lucorum and A. suturalis. Tetradecane was present in the volatiles of 10 common hosts, and their difference in relative content was significant. The presence of tetradecane seemed relevant to the olfactory response intensity of two mirid bugs towards the different host plants. The artificial supplement of tetradecane increased the attractive effect of host plants. These results suggested that tetradecane plays a vital role in the olfactory selection by two mirid bugs, and it can be made into field baits as a novel ecological strategy to manage these pests with widely reported pesticide resistance. However, results suggested host recognition is not entirely dependent on odor cues. We demonstrated that A. suturalis and A. lucorum adults have similar olfactory recognition mechanisms to their hosts in long-distance host selection. While, the differences in host plant selection between the two pests should occur in close range due to differences in gustatory or tactile sensory organs of A. lucorum and A. suturalis.

Functional Role of AsAP in the Reproduction of Adelphocoris suturalis (Hemiptera: Miridae)

Adelphocoris suturalis Jakovlev (Hemiptera: Miridae) is an omnivorous agricultural pest that has severe economic impacts on a diverse range of agricultural crops. Although the targeted disruption of reproductive development among insects has been proposed as a novel control strategy for pest species, the current understanding of the physiology and molecular mechanisms of A. suturalis reproduction is very limited. In this study, we isolated a putative A. suturalisaspartic protease (AsAP) gene that is highly expressed in the fat body and ovaries of sexually mature females. The double-stranded RNA (dsRNA)-mediated knockdown of AsAP suppressed ovarian development and negatively impacted female fertility, which suggested that it plays an essential role in A. suturalis reproduction. The results of this study could help to expand our understanding of A. suturalis reproductive development and have the potential to facilitate the development of effective strategies for the better control of this pest species.

Potential Acetylcholinesterase Inhibitor Acting on the Pesticide Resistant and Susceptible Cotton Pests

Gossypium spp., produces economically important cotton fiber, and its yield is highly affected due to pest attacks. Insecticidal target site mutation is one of the reasons behind insecticide resistance to a wide range of pesticides. Acetylcholinesterase (AChE) protein sequences from major pests of cotton were analyzed to assess various physicochemical properties, presence of motifs, and understand evolutionary relationship. The impact of three mutant AChE1, A. lucorum A216S, B. tabaci F392W, and A. gossypii A302S, on the strucutral stability was assessed, and F392W_AChE1 was selected based on 100 ns molecular dynamics simulation. Virtual screening of the zinc database and high-throughput virtual screening, standard precision, and extra precision docking resulted in the identification of six compounds. The six identified compounds and six known commercial pesticdes were docked with three mutant and three wild type AChE1, and one (C1) was selected based on Tice criteria. The conformational and interaction stability of the AChE1-C1 and F392W_AChE1-C1 complexes were monitored at 100 ns Gromacs simulation and were found to be thermodynamically favorable. Therefore, C1 may have the potential to bind to the resistant and susceptible strains of cotton pest, and the resistance developed by insects could be arrested. Furthermore, synthesis and field study of C1 will lead us to a better understanding of the efficacy of the identified compound.

Gut Bacterial Diversity in Different Life Cycle Stages of Adelphocoris suturalis (Hemiptera: Miridae)

Bacteria and insects have a mutually beneficial symbiotic relationship. Bacteria participate in several physiological processes such as reproduction, metabolism, and detoxification of the host. Adelphocoris suturalis is considered a pest by the agricultural industry and is now a major pest in cotton, posing a serious threat to agricultural production. As with many insects, various microbes live inside A. suturalis. However, the microbial composition and diversity of its life cycle have not been well-studied. To identify the species and community structure of symbiotic bacteria in A. suturalis, we used the HiSeq platform to perform high-throughput sequencing of the V3-V4 region in the 16S rRNA of symbiotic bacteria found in A. suturalis throughout its life stages. Our results demonstrated that younger nymphs (1st and 2nd instar nymphs) have higher species richness. Proteobacteria (87.06%) and Firmicutes (9.43%) were the dominant phyla of A. suturalis. At the genus level, Erwinia (28.98%), Staphylococcus (5.69%), and Acinetobacter (4.54%) were the dominant bacteria. We found that the relative abundance of Erwinia was very stable during the whole developmental stage. On the contrary, the relative abundance of Staphylococcus, Acinetobacter, Pseudomonas, and Corynebacterium showed significant dynamic changes at different developmental stages. Functional prediction of symbiotic bacteria mainly focuses on metabolic pathways. Our findings document symbiotic bacteria across the life cycle of A. suturalis, as well as differences in both the composition and richness in nymph and adult symbiotic bacteria. Our analysis of the bacteria in A. suturalis provides important information for the development of novel biological control strategies.

Characterization and functional analysis of two acetylcholinesterase genes in Bradysia odoriphaga Yang et Zhang (Diptera: Sciaridae)

Two acetylcholinesterase genes (Boace1 and Boace2) were cloned from Bradysia odoriphaga, a devastating soil pest that mainly damages Chinese chives. The Boace1 encodes BoAChE1 protein consisting of 696 amino acid residues, while Boace2 encodes BoAChE2 containing 638 amino acids. Phylogenetic analysis showed that Boace1 and Boace2 are appeared to be distinct clusters. The gene expression patterns at different development stages and various body parts tissues were examined, and their biological functions were characterized by RNA interference and analog docking prediction. The results showed that both Boace genes were expressed in all developmental stages and examined tissues. The transcript level of Boace2 was significantly higher than Boace1 in all tested samples, and Boace1 was found most abundant in the head while Boace2 was highly expressed in the fat body of B. odoriphaga. The silencing of Boace1 and Boace2 significantly decreased the AChE activity of 36.6% and 14.8% separately, and increased the susceptibility of B. odoriphaga to phoxim, with 60.8% and 44.7% mortality. Besides, overexpression and gene duplication of Boace1 were found in two field resistant populations, and two major mutations, A319S and G400V, were detected in Boace1. Moreover, the docking results revealed that BoAChE1 had a higher affinity towards organophosphorus than BoAChE2. It is concluded that Boace2 is the most abundant ace type in B. odoriphaga, while both Boace play vital roles. Boace1 might play a major neurological function and more likely be the prime target for insecticides, while Boace2 might play some important unidentified roles.

Genetic analysis and molecular detection of resistance to chlorpyrifos mediated by the A216S substitution in acetylcholinesterase-1 in the plant bug Apolygus lucorum

The green plant bug Apolygus lucorum is a major pest of Bacillus thuringiensis cotton in China. Previously, we reported that chlorpyrifos resistance in a laboratory-selected strain of A. lucorum (BZ-R) is associated with the homozygosis of an allele in the ace-1 gene encoding an alanine to serine substitution at position 216 of acetylcholinesterase-1. Here we describe the results of crosses between the resistant BZ-R strain (41-fold to chlorpyrifos) and the unselected susceptible BZ-S strain homozygous for the wild type alanine allele at position 216. Resistance to chlorpyrifos was inherited as a semi-dominant trait mainly controlled by a single autosomal gene and co-segregates strongly but not completely with the serine substitution in ace-1. Synergism bioassays and enzyme assays showed that minor contributions to resistance are also made by enhanced cytochrome P450 and carboxylesterase activities. A survey of 25 field populations from five Chinese provinces showed strong positive correlations between 50% lethal concentration against chlorpyrifos and S216 allele and genotype frequencies, although the most tolerant populations still only show 40%-50% S216 allele frequencies. The results above provide important information for designing effective resistance monitoring and management strategies for A. lucorum in China.

Selection and Validation of Suitable Reference Genes for RT-qPCR Analysis in Apolygus lucorum (Hemiptera: Miridae)

Apolygus lucorum (Meyer-Dür) is a destructive pest to >280 plants. Major economic significance and pesticide resistance issues have created a need for integrated pest management (e.g., RNAi, entomopathogen-based bioinsecticides) for A. lucorum. To better develop these control strategies, large-scale genetic studies involving gene-expression analysis are required and reverse transcription-quantitative polymerase chain reaction (RT-qPCR) is the most commonly used method. However, there have been no reports on appropriate reference genes in A. lucorum. Here, we evaluated nine widely utilized reference genes including EF1γ, RPL32, RPL27, SDH, TBP, ACT, ACT2, GAPDH, and βTUB for their expression stabilities in A. lucorum under five different conditions i.e., life stage, tissue, sex, dsRNA injection, and entomopathogen infection. Based on the gene stability ranking calculated by RefFinder, which integrates four algorithms (geNorm, delta Ct method, NormFinder, and BestKeeper), we recommend RPL27 and RPL32 as the most appropriate reference genes for molecular studies in different life stages and tissues; GAPDH and EF1γ for different sexes and entomopathogen infection studies; and RPL27 and EF1γ for RNAi studies. The results of this study will help improve the accuracy and reliability for normalizing the RT-qPCR data for further molecular analysis in A. lucorum.

Genome Analysis of Cytochrome in Dinotefuran-Treated Apolygus lucorum (Meyer-Dür)

In this study, two CYP genes, CYP395G1 and CYP4EY1, were analyzed in Apolygus lucorum (Hemiptera: Miridae). The expression pattern in different developmental stages of both CYP395G1 and CYP4EY1 revealed that first instar nymphs possessed the highest gene-transcript levels. After 12 h of dinotefuran treatments, the expression levels of CYP395G1 increased by 1.92-fold, while the CYP4EY1 expression decreased. It was observed that the sensitivities of dinotefuran to laboratory strains from F₀ to F₉ and the mRNA expression levels of CYP395G1 in the F₉ dinotefuran selected strain were higher when compared to the control strain. The results suggest that CYP395G1 was possibly a candidate P450 that was involved in dinotefuran detoxification. The dinotefuran resistant strain of A. lucorum was difficult to increase, due to its special molecular structure of dinotefuran compared with other neonicotinoids.

Transcript-Level Analysis of Detoxification Gene Mutation-Mediated Chlorpyrifos Resistance in Laodelphax striatellus (Hemiptera: Delphacidae)

Enhanced detoxification and target mutations that weaken insecticide binding ability are major mechanisms of insecticide resistance. Among these, over-expression or site mutations of carboxylesterase (CarE), cytochrome P450s (CYP450), and glutathione-S-transferase (GST) were the main form responsible for insecticide detoxification; however, transcript-level analysis of the relationship of detoxification gene mutations with chlorpyrifos (an organophosphorus insecticide) resistance is scarce thus far. In this study, multiple sites exhibiting polymorphisms within three detoxification genes were firstly examined via sequencing among different chlorpyrifos-resistant and susceptible individuals of Laodelphax striatellus. For example, the mutation frequencies of A374V in LsCarE16 were 83, 33, and 3%, S277A in LsCarE24 were 88, 28, and 3%, E36K in LsCYP426A1 were 100, 65, and 0% for chlorpyrifos-resistant, resistant decay, and susceptible individuals, respectively. Analysis also found expression levels of GSTd1, GSTt1, GSTs2, CYP4DE1U1, and CYP425B1 are coordinated with chlorpyrifos resistance levels; moreover, we found the deficiencies of 43S and 44A as well as two point mutations of E60D and Q61H at N-terminal region of the OP potential target acetylcholinesterase (AChE) in high resistant but not in low-chlorpyrifos resistant individuals. The results above all demonstrated the dynamic evolutionary process of insecticide resistance and revealed some resistance factors that only played roles at certain resistance level; high insecticide resistance in this example is the result of synergistic impact from multiple resistance factors.

Transcriptome analysis of three cotton pests reveals features of gene expressions in the mesophyll feeder Apolygus lucorum

The green mirid bug Apolygus lucorum is an agricultural pest that is known to cause damage to more than 150 plant species. Here, we report the transcriptomes of A. lucorum at three different developmental stages (the second and fifth instar nymphs and adults). A total of 98,236 unigenes with an average length of 1,335 nt was obtained, of which 50,640 were annotated, including those encoding digestive enzymes and cytochrome P450s. Comparisons with cotton bollworm and cotton aphid transcriptomes revealed distinct features of A. lucorum as a mesophyll feeder. The gene expression dynamics varied during development from young nymphs to adults. The high-quality transcriptome data and the gene expression dynamics reported here provide valuable data for a more comprehensive understanding of the physiology and development of mirid bugs, and for mining targets for their control.