Genomic adaptation to polyphagy and insecticides in a major East Asian noctuid pest

Chromosome-level genome of black cutworm provides novel insights into polyphagy and seasonal migration in insects

BACKGROUND

The black cutworm, Agrotis ipsilon, is a serious global underground pest. Its distinct phenotypic traits, especially its polyphagy and ability to migrate long distances, contribute to its widening distribution and increasing difficulty of control. However, knowledge about these traits is still limited.

RESULTS

We generated a high-quality chromosome-level assembly of A. ipsilon using PacBio and Hi-C technology with a contig N50 length of ~ 6.7 Mb. Comparative genomic and transcriptomic analyses showed that detoxification-associated gene families were highly expanded and induced after insects fed on specific host plants. Knockout of genes that encoded two induced ABC transporters using CRISPR/Cas9 significantly reduced larval growth rate, consistent with their contribution to host adaptation. A comparative transcriptomic analysis between tethered-flight moths and migrating moths showed expression changes in the circadian rhythm gene AiCry2 involved in sensing photoperiod variations and may receipt magnetic fields accompanied by MagR and in genes that regulate the juvenile hormone pathway and energy metabolism, all involved in migration processes.

CONCLUSIONS

This study provides valuable genomic resources for elucidating the mechanisms involved in moth migration and developing innovative control strategies.

De novo assembly of a chromosome-level reference genome of the ornamental butterfly Sericinus montelus based on nanopore sequencing and Hi-C analysis

Sericinus montelus (Lepidoptera, Papilionidae, Parnassiinae) is a high-value ornamental swallowtail butterfly species widely distributed in Northern and Central China, Japan, Korea, and Russia. The larval stage of this species feeds exclusively on Aristolochia plants. The Aristolochia species is well known for its high levels of aristolochic acids (AAs), which have been found to be carcinogenic for numerous animals. The swallowtail butterfly is among the few that can feed on these toxic host plants. However, the genetic adaptation of S. montelus to confer new abilities for AA tolerance has not yet been well explored, largely due to the limited genomic resources of this species. This study aimed to present a chromosome-level reference genome for S. montelus using the Oxford Nanopore long-read sequencing, Illumina short-read sequencing, and Hi-C technology. The final assembly was composed of 581.44 Mb with an expected genome size of 619.27 Mb. Further, 99.98% of the bases could be anchored onto 30 chromosomes. The N50 of contigs and scaffolds was 5.74 and 19.12 Mb, respectively. Approximately 48.86% of the assembled genome was suggested to be repeat elements, and 13,720 protein-coding genes were predicted in the current assembly. The phylogenetic analysis indicated that S. montelus diverged from the common ancestor of swallowtails about 58.57-80.46 million years ago. Compared with related species, S. montelus showed a significant expansion of P450 gene family members, and positive selections on eloa, heatr1, and aph1a resulted in the AA tolerance for S. montelus larva. The de novo assembly of a high-quality reference genome for S. montelus provided a fundamental genomic tool for future research on evolution, genome genetics, and toxicology of the swallowtail butterflies.

Genomic analyses provide insights into the genome evolution and environmental adaptation of the tobacco moth Ephestia elutella

Ephestia elutella is a major pest responsible for significant damage to stored tobacco over many years. Here, we conduct a comparative genomic analysis on this pest, aiming to explore the genetic bases of environmental adaptation of this species. We find gene families associated with nutrient metabolism, detoxification, antioxidant defense and gustatory receptors are expanded in the E. elutella genome. Detailed phylogenetic analysis of P450 genes further reveals obvious duplications in the CYP3 clan in E. elutella compared to the closely related species, the Indianmeal moth Plodia interpunctella. We also identify 229 rapidly evolving genes and 207 positively selected genes in E. elutella, respectively, and highlight two positively selected heat shock protein 40 (Hsp40) genes. In addition, we find a number of species-specific genes related to diverse biological processes, such as mitochondria biology and development. These findings advance our understanding of the mechanisms underlying processes of environmental adaptation on E. elutella and will enable the development of novel pest management strategies.

Diversity, duplication, and genomic organization of homeobox genes in Lepidoptera

Homeobox genes encode transcription factors with essential roles in patterning and cell fate in developing animal embryos. Many homeobox genes, including Hox and NK genes, are arranged in gene clusters, a feature likely related to transcriptional control. Sparse taxon sampling and fragmentary genome assemblies mean that little is known about the dynamics of homeobox gene evolution across Lepidoptera or about how changes in homeobox gene number and organization relate to diversity in this large order of insects. Here we analyze an extensive data set of high-quality genomes to characterize the number and organization of all homeobox genes in 123 species of Lepidoptera from 23 taxonomic families. We find most Lepidoptera have around 100 homeobox loci, including an unusual Hox gene cluster in which the lab gene is repositioned and the ro gene is next to pb A topologically associating domain spans much of the gene cluster, suggesting deep regulatory conservation of the Hox cluster arrangement in this insect order. Most Lepidoptera have four Shx genes, divergent zen-derived loci, but these loci underwent dramatic duplication in several lineages, with some moths having over 165 homeobox loci in the Hox gene cluster; this expansion is associated with local LINE element density. In contrast, the NK gene cluster content is more stable, although there are differences in organization compared with other insects, as well as major rearrangements within butterflies. Our analysis represents the first description of homeobox gene content across the order Lepidoptera, exemplifying the potential of newly generated genome assemblies for understanding genome and gene family evolution.

The oriental armyworm genome yields insights into the long-distance migration of noctuid moths

The oriental armyworm, Mythimna separata, is known for its long-distance seasonal migration and environment-dependent phase polymorphisms. Here, we present a chromosome-level genome reference and integrate multi-omics, functional genetics, and behavioral assays to explore the genetic bases of the hallmark traits of M. separata migration. Gene family comparisons show expansion of gustatory receptor genes in this cereal crop pest. Functional investigation of magnetoreception-related genes and associated flight behaviors suggest that M. separata may use the geomagnetic field to guide orientation in its nocturnal flight. Comparative transcriptome characterizes a suite of genes that may confer the observed plasticity between phases, including genes involved in protein processing, hormone regulation, and dopamine metabolism. We further report molecular signatures that underlie the dynamic regulation of a migratory syndrome coordinating reproduction and flight. Our study yields insights into environment-dependent developmental plasticity in moths and advances our understanding of long-distance migration in nocturnal insect pests.

A General Signal Pathway to Regulate Multiple Detoxification Genes Drives the Evolution of Helicoverpa armigera Adaptation to Xenobiotics

The study of insect adaptation to the defensive metabolites of host plants and various kinds of insecticides in order to acquire resistance is a hot topic in the pest-control field, but the mechanism is still unclear. In our study, we found that a general signal pathway exists in H. armigera which can regulate multiple P450s, GSTs and UGTs genes to help insects decrease their susceptibility to xenobiotics. Knockdown of HaNrf2 and HaAhR expression could significantly increase the toxicity of xenobiotics to H. armigera, and simultaneously decrease the gene expression of P450s, GSTs and UGTs which are related to the xenobiotic metabolism and synthesis of insect hormone pathways. Then, we used EMSA and dual luciferase assay to verify that a crosstalk exists between AhR and Nrf2 to regulate multiple P450s, GSTs and UGTs genes to mediate H. armigera susceptibility to plant allelochemicals and insecticides. The detoxification genes' expression network which can be regulated by Nrf2 and AhR is still unknown, and there were also no reports about the crosstalk between AhR and Nrf2 that exist in insects and can regulate multiple detoxification genes' expression. Our results provide a new general signaling pathway to reveal the adaptive mechanism of insects to xenobiotics and provides further insight into designing effective pest-management strategies to avoid the overuse of insecticides.

Reference Genome Sequences of the Oriental Armyworm, Mythimna separata (Lepidoptera: Noctuidae)

Lepidopteran insects are an important group of animals, including those used as biochemical and physiological model species in the insect and silk industries as well as others that are major agricultural pests. Therefore, the genome sequences of several lepidopteran insects have been reported. The oriental armyworm, Mythimna separata, is an agricultural pest commonly used to study insect immune reactions and interactions with parasitoid wasps as hosts. To improve our understanding of these research topics, reference genome sequences were constructed in the present study. Using long-read and short-read sequence data, de novo assembly and polishing were performed and haplotigs were purged. Subsequently, gene predictions and functional annotations were performed. To search for orthologs of the Toll and Immune Deficiency (IMD) pathways and for C-type lectins, annotation data analysis, BLASTp, and Hummer scans were performed. The M. separata genome is 682 Mbp; its contig N50 was 2.7 Mbp, with 21,970 genes and 24,452 coding sites predicted. All orthologs of the core components of the Toll and IMD pathways and 105 C-type lectins were identified. These results suggest that the genome data were of sufficient quality for use as reference genome data and could contribute to promoting M. separata and lepidopteran research at the molecular and genome levels.

Systematic identification and functional analysis of long noncoding RNAs involved in indoxacarb resistance in Spodoptera litura

Long noncoding RNAs (lncRNAs) are noncoding transcripts that are more than 200 nucleotides long. They play essential roles in regulating a variety of biological processes in many species, including insects, and some lncRNAs have been found to be associated with insecticide resistance. However, the characteristics and biological functions of lncRNAs involved in indoxacarb resistance are unknown in Spodoptera litura. We performed RNA sequencing in the SS, InRS, and FInRS of S. litura and identified 11 978 lncRNAs, including 3 136 intergenic lncRNAs, 7 393 intronic lncRNAs, and 1 449 anti-sense lncRNAs. Compared with the SS, 51 lncRNAs were upregulated and 134 lncRNAs were downregulated in the two resistant strains, and 908 differentially expressed mRNAs were predicted as the target genes of the 185 differentially expressed lncRNAs. Further analysis showed that 112 of differentially expressed lncRNAs may be associated with indoxacarb resistance by regulating the expression of 14 P450s, seven CCEs, one GST, six UGTs, five ABC transporters, and 24 cuticle protein genes, and 79 of differentially expressed lncRNAs may regulate the expression of 14 detoxification genes and 19 cuticle protein genes to participate in indoxacarb resistance by sponging 10 microRNAs. Interestingly, 47 of differentially expressed lncRNAs may mediate indoxacarb resistance through both lncRNA-mRNA and lncRNA-miRNA-mRNA regulatory pathways. Furthermore, quantitative PCR, RNA interference, and indoxacarb bioassay analyses indicated that overexpressed LNC_004867 and LNC_006576 were involved in indoxacarb resistance. This study provides comprehensive information for lncRNAs of S. litura, and presents evidence that lncRNAs have key roles in conferring insecticide resistance in S. litura.

Gut transcriptome analysis of P450 genes and cytochrome P450 reductase in three moth species feeding on gymnosperms or angiosperms

Cytochrome P450 enzymes (P450s, CYPs) are a superfamily of heme–thiolate proteins involved in the metabolism of endogenous and exogenous substances in insects. In this study, the identification of putative P450 proteins was done and the elimination of the repeated sequences resulted in 57 proteins from Gastropacha populifolia, 63 proteins from Dendrolimus punctatus, and 53 proteins from Dendrolimus tabulaeformis. The putative P450 proteins were aligned together with seven other insect species based on five conserved domains. A total of ten co-orthologous groups were identified. Interestingly, one co-orthologous gene, CYP4g15 in CYP4 clan, was identified and its 3D structure analysis showed that the highly conserved sites of the predicted motifs were close to the active sites of P450. Furthermore, this study revealed that insect CYP4g15 and two bacteria cytochrome P450 were monophyletic. This suggests that insects CYP4g15 are not only functionally conserved but also an ancient gene originating from different bacteria species.

Host-plant adaptation as a driver of incipient speciation in the fall armyworm (Spodoptera frugiperda)

Background

Divergent selection on host-plants is one of the main evolutionary forces driving ecological speciation in phytophagous insects. The ecological speciation might be challenging in the presence of gene flow and assortative mating because the direction of divergence is not necessarily the same between ecological selection (through host-plant adaptation) and assortative mating. The fall armyworm (FAW), a major lepidopteran pest species, is composed of two sympatric strains, corn and rice strains, named after two of their preferred host-plants. These two strains have been hypothesized to undergo incipient speciation, based on (i) several lines of evidence encompassing both pre- and post-zygotic reproductive isolation, and (ii) the presence of a substantial level of genetic differentiation. Even though the status of these two strains has been established a long time ago, it is still yet to be found whether these two strains indeed exhibit a marked level of genetic differentiation from a large number of genomic loci. Here, we analyzed whole genome sequences from 56 FAW individuals either collected from pasture grasses (a part of the favored host range of the rice strain) or corn to assess the role of host-plant adaptation in incipient speciation.

Results

Principal component analysis of whole genome data shows that the pattern of divergence in the fall armyworm is predominantly explained by the genetic differentiation associated with host-plants. The level of genetic differentiation between corn and rice strains is particularly marked in the Z chromosome. We identified one autosomal locus and two Z chromosome loci targeted by selective sweeps specific to rice strain and corn strain, respectively. The autosomal locus has both increased DXY and FST while the Z chromosome loci had decreased DXY and increased FST.

Conclusion

These results show that the FAW population structure is dominated by the genetic differentiation between corn and rice strains. This differentiation involves divergent selection targeting at least three loci, which include a locus potentially causing reproductive isolation. Taken together, these results suggest the evolutionary scenario that host-plant speciation is a driver of incipient speciation in the fall armyworm.

Effects of phytocompound Precocene 1 on the expression and functionality of the P450 gene in λ-cyhalothrin-resistant Spodoptera litura (Fab.)

Spodoptera litura (Fabricius) is an agriculturally significant polyphagous insect pest that has evolved a high level of resistance to conventional insecticides. A dietary assay was used in this work to assess the resilience of field populations of S. litura to λ-cyhalothrin. Analysis of the function and expression of the cytochrome P450 gene was used to test the sensitivity of S. litura larvae to sub-lethal concentrations of the insecticidal plant chemical Precocene 1, both by itself and in combination with λ-cyhalothrin. The activity of esterase enzymes (α and β) was found to decrease 48 h post treatment with Precocene 1. The activity of GST enzyme and cytochrome P450 increased with Precocene 1 treatment post 48 h, however. Expression studies revealed the modulation by Precocene 1 of cytochrome P450 genes, CYP4M16, CYP4M15, CYP4S8V4, CYP4G31, and CYP4L10. While CYP4M16 expression was stimulated the most by the synergistic Precocene 1 + λ-cyhalothrin treatment, expression of CYP4G31 was the most down-regulated by Precocene 1 exposure. Hence, it is evident that λ-cyhalothrin-resistant pest populations are still sensitive to Precocene 1 at a sublethal concentration that is nevertheless capable of hindering their development. Precocene 1 can therefore be considered a potent candidate for the effective management of insecticide-resilient S. litura.

Pathogenicity of Beauveria bassiana PfBb and Immune Responses of a Non-Target Host, Spodoptera frugiperda (Lepidoptera: Noctuidae)

Exploring the pathogenicity of a new fungus strain to non-target host pests can provide essential information on a large scale for potential application in pest control. In this study, we tested the pathogenicity of Beauveria bassiana PfBb on the important agricultural pest Spodoptera frugiperda (Lepidoptera: Noctuidae) by determining the relative activities of protective enzymes and detoxifying enzymes in different larval instars. Our results show that the B. bassiana PfBb strain could infect all six larval instars of S. frugiperda, and its virulence to S. frugiperda larvae gradually increased with an increase in spore concentration. Seven days after inoculation, the LC₅₀ of B. bassiana PfBb was 7.7 × 10⁵, 5.5 × 10⁶, 2.2 × 10⁷, 3.1 × 10⁸, 9.6 × 10⁸, and 2.5 × 10¹¹ spores/mL for first to sixth instars of S. frugiperda, respectively, and the LC₅₀ and LC₉₀ of B. bassiana PfBb for each S. frugiperda instar decreased with infection time, indicating a significant dose effect. Furthermore, the virulence of B. bassiana PfBb to S. frugiperda larvae gradually decreased with an increase in larval instar. The activities of protective enzymes (i.e., catalase, peroxidase, and superoxide dismutase) and detoxifying enzymes (i.e., glutathione S-transferases, carboxylesterase, and cytochrome P450) in S. frugiperda larvae of the first three instars infected with B. bassiana PfBb changed significantly with infection time, but such variations were not obvious in the fifth and sixth instars. Additionally, after being infected with B. bassiana PfBb, the activities of protective enzymes and detoxification enzymes in S. frugiperda larvae usually lasted from 12 to 48 h, which was significantly longer than the control. These results indicate that the pathogenicity of B. bassiana PfBb on the non-target host S. frugiperda was significant but depended on the instar stage. Therefore, the findings of this study suggest that B. bassiana PfBb can be used as a bio-insecticide to control young larvae of S. frugiperda in an integrated pest management program.

Functional analysis of a bitter gustatory receptor highly expressed in the larval maxillary galea of Helicoverpa armigera

Many plant secondary substances are feeding deterrents for insects and play a key role in the selection of host plants. The taste sensilla of phytophagous insects contain gustatory sensory neurons sensitive to deterrents but the molecular basis of deterrent chemoreception remains unknown. We investigated the function of Gr180, the most highly expressed bitter gustatory receptor in the maxillary galea of Helicoverpa armigera larvae. Functional analyses using the Xenopus oocyte expression system and two-electrode voltage clamp revealed that the oocytes expressing Gr180 responded to coumarin. Tip recording results showed that the medial sensilla styloconica of the maxilla of fifth instar larvae exhibited electrophysiological responses to coumarin. Two-choice feeding bioassays confirmed that coumarin inhibited larval feeding. A homozygous mutant strain of H. armigera with truncated Gr180 proteins (Gr180-/-) was established using the CRISPR-Cas9 system. The responses of the medial sensilla styloconica in Gr180-/- to coumarin were almost abolished, and the responses to sinigrin and strychnine were also significantly decreased. Knockout of Gr180 alleviated the feeding deterrent effects of coumarin, sinigrin, and strychnine. Thus, we conclude that Gr180 is a receptor responding to coumarin,and also participates in sensing sinigrin and strychnine. These results enhance our understanding of the gustatory sensing mechanisms of phytophagous insects to deterrents.

Two Odorant-Binding Proteins Involved in the Recognition of Sex Pheromones in Spodoptera litura Larvae

Usually, the recognition of sex pheromone signals is restricted to adult moths. Here, our behavioral assay showed that fourth-instar Spodoptera litura larvae are attracted to cabbage laced with minor sex pheromones Z9,E12-tetradecadienyl acetate (Z9,E12-14:Ac) or Z9-tetradecenyl acetate (Z9-14:Ac). Seven odorant-binding proteins (OBPs) were upregulated after exposure to Z9,E12-14:Ac, and one OBP was upregulated after exposure to Z9-14:Ac. Fluorescence competitive binding assays showed that GOBP2 and OBP7 bound to sex pheromones. RNAi treatment significantly downregulated GOBP2 and OBP7 mRNA expression by 70.37 and 63.27%, respectively. The siOBP-treated larvae were not attracted to Z9,E12-14:Ac or Z9-14:Ac, and the corresponding preference indices were significantly lower than those in siGFP-treated larvae. Therefore, we concluded that GOBP2 and OBP7 are involved in the attraction of S. litura larvae to food containing Z9,E12-14:Ac and Z9-14:Ac. These results provide an important basis for exploring the olfactory mechanisms underlying sex pheromone attraction in moth larvae.

Identification and expression profiles of gustatory receptor genes in Bactrocera minax larvae (Diptera: Tephritidae): Role of BminGR59f in larval growth

Insects employ various types of gustatory receptors (GRs) to identify nutrient-rich food and avoid toxic substances. The larval gustatory system is the critical checkpoint for food acceptance or rejection. As a specialist herbivore, the larvae of Bactrocera minax feed only on unripe citrus fruits. However, how larvae use GRs to check and adapt to the secondary metabolites in unripe citrus fruits remains unknown. In this study, we first performed developmental expression profiles showing that most BminGRs genes were highly expressed in 1st and 2nd instar larvae and that tissue-specific expression indicated high expression of most BminGRs genes in the mouthparts of 2nd instar larvae. Furthermore, we found that silencing BminGR59f by RNA interference (RNAi) affected the growth of 2nd instar B. minax larvae. Hesperidin and naringin were screened as ligands of BminGR59f via RNAi and cell calcium imaging, and the combination of these two flavones increased the body weight of larvae. In summary, we identified a novel gustatory perception pattern in B. minax for detecting hesperidin and naringin, which boosted the growth of B. minax larvae. These results shed light on how specialist herbivores detect and adapt to host metabolites in adverse environments depending on larval GRs.

Activation of the ROS/CncC Signaling Pathway Regulates Cytochrome P450 CYP4BQ1 Responsible for (+)-α-Pinene Tolerance in Dendroctonus armandi

Bark beetles mainly rely on detoxification enzymes to resist the host tree's defense against oleoresin terpenes. Cytochrome P450 enzymes (CYPs) play an important role in the detoxification of plant allelochemicals and pesticides in insect. One P450 gene (DaCYP4BQ1) is associated with the response of (+)-α-pinene in Dendroctonus armandi. However, the regulatory mechanism of this P450 gene response to (+)-α-pinene is still unknown. In this study, spatiotemporal expression profiling indicated that CYP4BQ1 was highly expressed in adult and larval stages of D. armandi, and it was predominantly expressed in fat body, midgut, and Malpighian tubules of adults. Moreover, the expression of CYP4BQ1 significantly increased after exposure to (+)-α-pinene, and depletion of it decreased the tolerance of adults to (+)-α-pinene. In addition, (+)-α-pinene treatment induced the expression of the transcription factors cap 'n' collar isoform C (CncC) and its binding factor muscle aponeurosis fibromatosis (Maf), elevated the level of hydrogen peroxide (H₂O₂), and increased the activities of antioxidant enzymes. Silencing CncC suppressed CYP4BQ1 expression and enhanced the susceptibility of beetles to (+)-α-pinene. Similarly, application of the reactive oxygen species (ROS) scavenger N-acetylcysteine reduced the production and accumulation of H₂O₂, suppressed the expression of CncC, Maf, and CYP4BQ1 and led to decreased tolerance of adults to (+)-α-pinene. In contrast, ingestion of the CncC agonist curcumin elevated CYP4BQ1 expression and enhanced (+)-α-pinene tolerance. The results demonstrate that, in D. armandi, (+)-α-pinene induces CYP4BQ1 via activation of the ROS/CncC signaling pathway.

Identification and characterization of candidate detoxification genes in Pharsalia antennata Gahan (Coleoptera: Cerambycidae)

The wood-boring beetles, including the majority of Cerambycidae, have developed the ability to metabolize a variety of toxic compounds derived from host plants and the surrounding environment. However, detoxification mechanisms underlying the evolutionary adaptation of a cerambycid beetle Pharsalia antennata to hosts and habitats are largely unexplored. Here, we characterized three key gene families in relation to detoxification (cytochrome P450 monooxygenases: P450s, carboxylesterases: COEs and glutathione-S-transferases: GSTs), by combinations of transcriptomics, gene identification, phylogenetics and expression profiles. Illumina sequencing generated 668,701,566 filtered reads in 12 tissues of P. antennata, summing to 100.28 gigabases data. From the transcriptome, 215 genes encoding 106 P450s, 77 COEs and 32 GSTs were identified, of which 107 relatives were differentially expressed genes. Of the identified 215 genes, a number of relatives showed the orthology to those in Anoplophora glabripennis, revealing 1:1 relationships in 94 phylogenetic clades. In the trees, P. antennata detoxification genes mainly clustered into one or two subfamilies, including 64 P450s in the CYP3 clan, 33 COEs in clade A, and 20 GSTs in Delta and Epsilon subclasses. Combining transcriptomic data and PCR approaches, the numbers of detoxification genes expressed in abdomens, antennae and legs were 188, 148 and 141, respectively. Notably, some genes exhibited significantly sex-biased levels in antennae or legs of both sexes. The findings provide valuable reference resources for further exploring xenobiotics metabolism and odorant detection in P. antennata.

CRISPR/Cas9-Mediated Knockout Reveals the Involvement of CYP304F1 in β-Cypermethrin and Chlorpyrifos Resistance in Spodoptera litura

Functions of insect CYP2 clan P450s in insecticide resistance are relatively less reported. In Spodoptera litura, a gene from the CYP2 clan (CYP304F1) was validated to be up-regulated significantly in a pyrethroid- and organophosphate-resistant population (QJ) than a susceptible population by RNA-Seq and qRT-PCR. Spatial-temporal expression indicated the high expression of CYP304F1 in the fourth, fifth, and sixth instar larvae and the metabolism-related tissue fat body and malpighian tubules. CYP304F1 was knocked out by CRISPR/Cas9, and a homozygous population (QJ-CYP304F1) with a G-base deletion at exon 2 was obtained after selection. Bioassay results showed that the LD₅₀ values to β-cypermethrin and chlorpyrifos in the QJ-CYP304F1 population decreased significantly, and the resistance ratio was both 1.81-fold in the QJ population compared with that in the QJ-CYP304F1 population. The toxicity of fenvalerate, cyhalothrin, or phoxim showed no significant change. These results suggested that CYP304F1 is involved in β-cypermethrin and chlorpyrifos resistance in S. litura.

Detoxification Gene Families at the Genome-Wide Level of Rhus Gall Aphid Schlechtendalia chinensis

The Rhus gall aphid Schlechtendalia chinensis uses the species Rhus chinensis as its primary host plant, on which galls are produced. The galls have medicinal properties and can be used in various situations due to their high tannin content. Detoxification enzymes play significant roles in the insect lifecycle. In this study, we focused on five detoxification gene families, i.e., glutathione-S-transferase (GST), ABC transporter (ABC), Carboxylesterase (CCE), cyto-chrome P450 (CYP), and UDP-glycosyltransferase (UDP), and manually annotated 144 detoxification genes of S. chinensis using genome-wide techniques. The detoxification genes appeared mostly on chromosome 1, where a total of two pair genes were identified to show tandem duplications. There were 38 gene pairs between genomes of S. chinensis and Acyrthosiphon pisum in the detoxification gene families by collinear comparison. Ka/Ks ratios showed that detoxification genes of S. chinensis were mainly affected by purification selection during evolution. The gene expression numbers of P450s and ABCs by transcriptome sequencing data were greater, while gene expression of CCEs was the highest, suggesting they might be important in the detoxification process. Our study has firstly identified the genes of the different detoxification gene families in the S. chinensis genome, and then analyzed their general features and expression, demonstrating the importance of the detoxification genes in the aphid and providing new information for further research.

Haplotype-resolved genome assembly of Bletilla striata (Thunb.) Reichb.f. to elucidate medicinal value

Bletilla striata, commonly known as baiji, is a species used in traditional Chinese medicine; it is highly regarded for its medicinal applications and therefore has high economic value. Here, we report a high-quality haplotype-resolved genome of B. striata, haplotype A (2.37 Gb, with a scaffold N50 of 146.39 Mb and a contig N50 of 1.65 Mb) and haplotype B (2.43 Gb, with a scaffold N50 of 150.22 Mb and a contig N50 of 1.66 Mb), assembled from high-fidelity (HiFi) reads and chromosome conformation capture (Hi-C) reads. We find evidence that B. striata has undergone two whole-genome duplication (WGD) events: an ancient WGD event shared by most monocots and a recent WGD event unique to all orchids. We also reconstructed the ancestral orchid karyotype (AOK) of 18 ancient chromosomes and the evolutionary trajectories of 16 modern B. striata chromosomes. Comparative genomic analysis suggests that the expanded gene families of B. striata might play important roles in secondary metabolite biosynthesis and environmental adaptation. By combining genomic and transcriptomic data, we identified the 10 core members from nine gene families that were probably involved in B. striata polysaccharide (BSP) biosynthesis. Based on virus-induced gene silencing (VIGS) and yeast two-hybrid experiments, we present an MYB transcription factor (TF), BsMYB2, that can regulate BSP biosynthesis by directly interacting with eight key BSP-related genes: sacA1, HK1, scrK1, scrK2, GPI1, manA1, GMPP1 and UGP2_1. Our study will enhance the understanding of orchid evolution and accelerate the molecular-assisted breeding of B. striata for improving traits of medicinal value.

Genome-Wide Analysis of Odorant and Gustatory Receptors in Six Papilio Butterflies (Lepidoptera: Papilionidae)

The chemical interactions of insects and host plants are shaping the evolution of chemosensory receptor gene families. However, the correlation between host range and chemoreceptor gene repertoire sizes is still elusive in Papilionidae. Here, we addressed the issue of whether host plant diversities are correlated with the expansions of odorant (ORs) or gustatory (GRs) receptors in six Papilio butterflies. By combining genomics, transcriptomics and bioinformatics approaches, 381 ORs and 328 GRs were annotated in the genomes of a generalist P. glaucus and five specialists, P. xuthus, P. polytes, P. memnon, P. machaon and P. dardanus. Orthologous ORs or GRs in Papilio had highly conserved gene structure. Five Papilio specialists exhibited a similar frequency of intron lengths for ORs or GRs, but which was different from those in the generalist. Phylogenetic analysis revealed 60 orthologous OR groups, 45 of which shared one-to-one relationships. Such a single gene in each butterfly also occurred in 26 GR groups. Intriguingly, bitter GRs had fewer introns than other GRs and clustered into a large clade. Focusing on the two chemoreceptor gene families in P. xuthus, most PxutORs (52/58) were expressed in antennae and 31 genes in reproductive tissues. Eleven out of 28 foretarsus-expressed PxutGRs were female-biased genes, as strong candidates for sensing oviposition stimulants. These results indicate that the host range may not shape the large-scale expansions of ORs and GRs in Papilio butterflies and identify important molecular targets involved in olfaction, oviposition or reproduction in P. xuthus.

CRISPR/Cas9-mediated Serine protease 2 disruption induces male sterility in Spodoptera litura

Male fertility is essential for reproduction and population growth in animals. Many factors affect male fertility, such as courtship behavior, sperm quantity, and sperm motility, among others. Seminal Fluid Proteins (SFPs) are vital components of seminal fluid in the male ejaculate, which affect male fertility, sperm activation, and female ovulation. However, the knowledge of SFPs is insufficient; the function of many SFPs remains unknown, and most described functions were mainly characterized in Drosophila or other laboratory models. Here, we focus on the Serine protease 2 (Ser2) gene in the lepidopteran pest Spodoptera litura. The Ser2 gene was specifically expressed in male adults. Disruption of the Ser2 gene mediated by CRISPR/Cas9 induced male sterility but females remained fertile. PCR-based detection of the next-generation mutants showed that male sterility was stably inherited. The qRT-PCR analysis of SlSer2 mutants showed that motor protein family genes and structural protein family genes were down-regulated, while protein modification family genes were up-regulated, suggesting that SlSer2 may be involved in sperm movement and activity. These results demonstrate that Ser2 is an important component of SFPs in seminal fluid and was identified for a useful sterile gene for pest control that may lead to new control strategies for lepidopteran insect pests such as S. litura.

Genome-wide identification and characterization of basic helix-loop-helix transcription factors in Spodoptera litura upon pathogen infection

Basic helix-loop-helix (bHLH) transcription factors play an important role in a wide range of metabolic and developmental processes in eukaryotes, and bHLH proteins also participate in immune responses, especially in plants. However, their roles in insects upon entomopathogen infection are unknown. In this study, 54 bHLH genes in 41 families were identified in a polyphagous pest, Spodoptera litura, including a new bHLH gene in group B, which is specifically present in Lepidoptera and was thus named Lep. The conserved amino acids in the bHLH domain, structural architecture, and chromosomal distribution of bHLH genes in S. litura were analyzed. The bHLH genes in Plutella xylostella and Apis mellifera were also updated, and genome-wide comparison and phylogenetic analysis of bHLH members in 5 holometabolous insects were performed. The expression profiles of S. litura bHLH (SlbHLH) genes in 3 tissues at different developmental stages and their responses to S. litura nucleopolyhedrovirus (SpltNPV), Nomuraea rileyi (Nr), and Bacillus thuringiensis (Bt) infection were investigated. More SlbHLHs in group B were expressed and differentially expressed during pathogen infections, and SlbHLHs tended to be downregulated in the midgut of S. litura larvae after B. thuringiensis treatment. Our study provides an overview of bHLH family members in S. litura and their responses to different pathogens used for pest biocontrol. These findings on bHLH members may contribute to uncovering the mechanism of host-pathogen interaction.

Three cytochrome P450 CYP4 family genes regulated by the CncC signaling pathway mediate phytochemical susceptibility in the red flour beetle, Tribolium castaneum

BACKGROUND

Insect cytochrome P450 monooxygenases (P450s) play a crucial role in phytochemical metabolism and tolerance. Three P450 genes (TcCYP4Q3, TcCYP4Q5, and TcCYP4Q7) are associated with the response of eugenol in Tribolium castaneum. However, the responding mechanisms of these P450 genes to eugenol remain unknown.

RESULTS

Here, spatiotemporal expression profiling revealed that TcCYP4Q3 and TcCYP4Q5 were most highly expressed in late adult, while TcCYP4Q7 was predominantly expressed in late larva; and all of these three P450 genes were mainly expressed in the fat body of larvae. Furthermore, the expressions of these three P450 genes were significantly up-regulated after exposure to eugenol, and depletion of them enhanced the susceptibility of beetles to eugenol. Interestingly, RNA interference (RNAi) against the CncC gene, a transcription factor of CncC signaling pathway associated with regulation of insect P450s in response to phytochemicals, reduced the transcripts of these three P450 genes following exposure to eugenol. Investigation of CncC signaling pathway showed that this pathway could be activated by eugenol.

CONCLUSION

Altogether, the results indicate that these three P450 genes are regulated by CncC signaling pathway to participate in the susceptibility of Tribolium castaneum to phytochemicals. These findings will aid implications for the development of novel therapeutics to control pest. © 2022 Society of Chemical Industry.

Spodoptera littoralis genome mining brings insights on the dynamic of expansion of gustatory receptors in polyphagous noctuidae

The bitter taste, triggered via gustatory receptors, serves as an important natural defense against the ingestion of poisonous foods in animals, and the increased host breadth is usually linked to an increase in the number of gustatory receptor genes. This has been especially observed in polyphagous insect species, such as noctuid species from the Spodoptera genus. However, the dynamic and physical mechanisms leading to these gene expansions and the evolutionary pressures behind them remain elusive. Among major drivers of genome dynamics are the transposable elements but, surprisingly, their potential role in insect gustatory receptor expansion has not been considered yet. In this work, we hypothesized that transposable elements and possibly positive selection would be involved in the highly dynamic evolution of gustatory receptor in Spodoptera spp. We first sequenced de novo the full 465 Mb genome of S. littoralis, and manually annotated the main chemosensory genes, including a large repertoire of 373 gustatory receptor genes (including 19 pseudogenes). We also improved the completeness of S. frugiperda and S. litura gustatory receptor gene repertoires. Then, we annotated transposable elements and revealed that a particular category of class I retrotransposons, the SINE transposons, was significantly enriched in the vicinity of gustatory receptor gene clusters, suggesting a transposon-mediated mechanism for the formation of these clusters. Selection pressure analyses indicated that positive selection within the gustatory receptor gene family is cryptic, only 7 receptors being identified as positively selected. Altogether, our data provide a new good quality Spodoptera genome, pinpoint interesting gustatory receptor candidates for further functional studies and bring valuable genomic information on the mechanisms of gustatory receptor expansions in polyphagous insect species.

Chromosomal-level reference genome of the moth Heortia vitessoides (Lepidoptera: Crambidae), a major pest of agarwood-producing trees

The moth Heortia vitessoides Moore (Lepidoptera: Crambidae) is a major pest of ecologically, commercially and culturally important agarwood-producing trees in the genus Aquilaria. In particular, H. vitessoides is one of the most destructive defoliating pests of the incense tree Aquilaria sinesis, which produces a valuable fragrant wood used as incense and in traditional Chinese medicine [33]. Nevertheless, a genomic resource for H. vitessoides is lacking. Here, we present a chromosomal-level assembly for H. vitessoides, consisting of a 517 megabase (Mb) genome assembly with high physical contiguity (scaffold N50 of 18.2 Mb) and high completeness (97.9% complete BUSCO score). To aid gene annotation, 8 messenger RNA transcriptomes from different developmental stages were generated, and a total of 16,421 gene models were predicted. Expansion of gene families involved in xenobiotic metabolism and development were detected, including duplications of cytosolic sulfotransferase (SULT) genes shared among lepidopterans. In addition, small RNA sequencing of 5 developmental stages of H. vitessoides facilitated the identification of 85 lepidopteran conserved microRNAs, 94 lineage-specific microRNAs, as well as several microRNA clusters. A large proportion of the H. vitessoides genome consists of repeats, with a 29.12% total genomic contribution from transposable elements, of which long interspersed nuclear elements (LINEs) are the dominant component (17.41%). A sharp decrease in the genome-wide percentage of LINEs with lower levels of genetic distance to family consensus sequences suggests that LINE activity has peaked in H. vitessoides. In contrast, opposing patterns suggest a substantial recent increase in DNA and LTR element activity. Together with annotations of essential sesquiterpenoid hormonal pathways, neuropeptides, microRNAs and transposable elements, the high-quality genomic and transcriptomic resources we provide for the economically important moth H. vitessoides provide a platform for the development of genomic approaches to pest management, and contribute to addressing fundamental research questions in Lepidoptera.

Genome-Wide Analysis of Odorant-Binding Proteins and Chemosensory Proteins in the Bean bug Riptortus pedestris

Insects have sensitive olfactory systems to interact with environment and respond to the change in host plant conditions. Key genes in the system can be potential targets for developing new and efficient pest behaviour control methods. Riptortus pedestris is an important soybean pest in East Asia and has caused serious damage to the soybean plants in Huang-Huai-Hai region of China. However, the current treatment of pests is dominated by chemical insecticides and lacks efficient sustainable prevention and control technologies. In this study, we identified 49 putative odorant-binding proteins (OBPs) (43 were new genes) and 25 chemosensory proteins (CSPs) (17 were new genes) in R. pedestris genome. These OBP and CSP genes are clustered in highly conserved groups from other hemipteran species in phylogenetic trees. Most RpedOBPs displayed antennal-biased expression. Among the 49 RpedOBPs, 33 were significantly highly expressed in the antennae, including three male-biased and nine female-biased. While many RpedCSPs were detected both in the antennae and in non-antennal tissues, only 11 RpedCSPs displayed antennal-biased expression, in which four RpedCSPs were male-biased and five RpedCSPs were female-biased. Some OBP and CSP genes showed sex-biased expression profiles. Our results not only provide a foundation for future exploration of the functions of RpedOBPs and RpedCSPs but also aid in developing environmentally friendly insecticides in the future.

Genome-Wide Identification and Stage-Specific Expression Profile Analysis Reveal the Function of Ribosomal Proteins for Oogenesis of Spodoptera litura

Ribosomal proteins (Rps) are indispensable in ribosome biogenesis and protein synthesis, which tightly correlate with cell growth and proliferation in different physiological processes across species. Up to now, genes coding for Rps have been identified and studied in many species, however, their information still remains elusive in many insect species, especially in Spodoptera litura. In this study, 81 Rp genes were identified from S. litura genome and were mapped to their positions on the chromosomes. In addition, their physical and chemical properties, gene structure, phylogenetic relationships, targeted microRNAs were also analyzed. Gene ontology analysis disclosed that Rp genes were closely associated with processes related to ribosome biosynthesis, proteins translation processing, molecular binding activities. The quantitative real-time PCR (qRT-PCR) revealed expression profiles of Rp genes varied in different stages of oogenesis, and found that most Rp genes accumulated in previtellogenesis stage. This study described the comprehensive genome-wide analysis of Rp gene family in agricultural pests, which provided foundation for further characterizing the roles of Rps in oogenesis of insects, and some Rp genes may further serve as targets for innovative pest control.

Pesticide resistance in arthropods: Ecology matters too

Pesticide resistance development is an example of rapid contemporary evolution that poses immense challenges for agriculture. It typically evolves due to the strong directional selection that pesticide treatments exert on herbivorous arthropods. However, recent research suggests that some species are more prone to evolve pesticide resistance than others due to their evolutionary history and standing genetic variation. Generalist species might develop pesticide resistance especially rapidly due to pre‐adaptation to handle a wide array of plant allelochemicals. Moreover, research has shown that adaptation to novel host plants could lead to increased pesticide resistance. Exploring such cross‐resistance between host plant range evolution and pesticide resistance development from an ecological perspective is needed to understand its causes and consequences better. Much research has, however, been devoted to the molecular mechanisms underlying pesticide resistance while both the ecological contexts that could facilitate resistance evolution and the ecological consequences of cross‐resistance have been under‐studied. Here, we take an eco‐evolutionary approach and discuss circumstances that may facilitate cross‐resistance in arthropods and the consequences cross‐resistance may have for plant–arthropod interactions in both target and non‐target species and species interactions. Furthermore, we suggest future research avenues and practical implications of an increased ecological understanding of pesticide resistance evolution.

Bracovirus Sneaks Into Apoptotic Bodies Transmitting Immunosuppressive Signaling Driven by Integration-Mediated eIF5A Hypusination

A typical characteristics of polydnavirus (PDV) infection is a persistent immunosuppression, governed by the viral integration and expression of virulence genes. Recently, activation of caspase-3 by Microplitis bicoloratus bracovirus (MbBV) to cleave Innexins, gap junction proteins, has been highlighted, further promoting apoptotic cell disassembly and apoptotic body (AB) formation. However, whether ABs play a role in immune suppression remains to be determined. Herein, we show that ABs transmitted immunosuppressive signaling, causing recipient cells to undergo apoptosis and dismigration. Furthermore, the insertion of viral–host integrated motif sites damaged the host genome, stimulating eIF5A nucleocytoplasmic transport and activating the eIF5A-hypusination translation pathway. This pathway specifically translates apoptosis-related host proteins, such as P53, CypA, CypD, and CypJ, to drive cellular apoptosis owing to broken dsDNA. Furthermore, translated viral proteins, such Vank86, 92, and 101, known to complex with transcription factor Dip3, positively regulated DHYS and DOHH transcription maintaining the activation of the eIF5A-hypusination. Mechanistically, MbBV-mediated extracellular vesicles contained inserted viral fragments that re-integrated into recipients, potentially via the homologous recombinant repair system. Meanwhile, this stimulation regulated activated caspase-3 levels via PI3K/AKT 308 and 473 dephosphorylation to promote apoptosis of granulocyte-like recipients Sf9 cell; maintaining PI3K/AKT 473 phosphorylation and 308 dephosphorylation inhibited caspase-3 activation leading to dismigration of plasmatocyte-like recipient High Five cells. Together, our results suggest that integration-mediated eIF5A hypusination drives extracellular vesicles for continuous immunosuppression.

Lacking of sex-lethal gene lowers the fertility of male reproduction in Spodoptera litura (Lepidoptera)

Sex-lethal (Sxl) encodes an RNA-binding protein that acts as the switch of sex determination in Drosophila and influences the genitalia formation and gonadal development. However, its sex-determination roles are not conserved in all insects and its role in the gonadal development of Lepidoptera is not well documented. In this study, three splicing variants of Sxl mRNA were identified in Spodoptera litura and they highly expressed in gonads, particularly in the testis. The mRNA levels of SlSxl exhibited higher expression in the spermatid than the testis sheaths, and gradually increased with the spermiogenesis. Sex-lethal protein (SlSXL) is mainly distributed in the cytoplasm of spermatocytes and the head of spermatid. Knockout of SlSxl resulted in fewer eupyrene sperm bundles and apyrene sperm bundles in the testes of moth and a large number of undeveloped spermatocysts retained in the moth of mutant testis, and leading to the reduction of oviposition and hatch rate in the offsprings after mating with female. These results suggest that SlSxl is a critical player in the spermiogenesis of S. litura.

Genome-wide identification and comparative analysis of Cry toxin receptor families in 7 insect species with a focus on Spodoptera litura

Cadherin, aminopeptidase N (APN) and alkaline phosphatase (ALP) have been characterized as Cry receptors. In this study, comparative genomic analysis of the 3 receptor families was performed in 7 insects. ALPs and APNs are divided into three and eight clades in phylogenetic trees, respectively. ALPs in clade 3 and APNs in clade 1 contain multiple paralogs within each species and most paralogs are located closely in chromosomes. Drosophila melanogaster has expanded APNs in clade 5 and were lowly expressed in midgut. Cadherins are divided into 16 clades; they may diverge before holometabolous insect speciation except for BtR and Cad89D-like clades. Eight insects from different orders containing BtR orthologs are sensitive to Cry1A or Cry3A, while five species without BtR are insensitive to both toxins. Most APNs in clade 1, several ALPs in clade 3, BtR and Cad89D-like genes were highly or moderately expressed in larval midgut of Spodoptera litura and the other six species, and several members in these clades have been identified as Cry receptors. Expressions of putative S. litura Cry receptors in the midgut after exposing to Bt toxins were also analyzed.

Transcriptomics and metagenomics of common cutworm (Spodoptera litura) and fall armyworm (Spodoptera frugiperda) demonstrate differences in detoxification and development

BACKGROUND

Spodoptera litura is an important polyphagous pest that causes significant damage to the agricultural sector. We performed RNA-seq of 15 S. litura individuals from larval (fifth and sixth instar larvae), chrysalis, and adult developmental stages. We also compared the S. litura transcriptome data with Spodoptera frugiperda across the same developmental stages, which was sequenced in our previous study.

RESULTS

A total of 101,885 differentially expressed transcripts (DETs) were identified in S. litura. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses indicated that S. litura may undergo active xenobiotic and detoxifying metabolism during its larval and adult stages, which may explain difficulties with current population control measures. We also found that DETs of single-copy orthologous genes between S. litura and S. frugiperda were involved in basic metabolism and development. However, energy and metabolic processes genes had a higher expression in S. litura, whereas nervous and olfactory function genes had a higher expression in S. frugiperda. Metagenomics analysis in larval S. litura and S. frugiperda revealed that microbiota participate in the detoxification and metabolism processes, but the relative abundance of detoxification-related microbiota was more abundant in S. frugiperda. Transcriptome results also confirmed the detoxification-related pathway of S. frugiperda was more abundant than in S. litura.

CONCLUSIONS

Significant changes at transcriptional level were identified during the different development stages of S. litura. Importantly, we also identified detoxification associated genes and gut microbiota between S. litura and S. frugiperda at different developmental stages, which will be valuable in revealing possible mechanisms of detoxification and development in these two lepidopterans.

Genomic features of the polyphagous cotton leafworm Spodoptera littoralis

BACKGROUND

The cotton leafworm, Spodoptera littoralis, is a highly polyphagous pest of many cultivated plants and crops in Africa and Europe. The genome of this pest will help us to further understand the molecular mechanisms of polyphagy.

RESULTS

Herein, the high-quality genome of S. littoralis was obtained by Pacific Bioscience (PacBio) sequencing. The assembled genome size of S. littoralis is 436.55 Mb with a scaffold N50 of 6.09 Mb, consisting of 17,207 annotated protein-coding genes. Phylogenetic analysis shows that S. littoralis and its sibling species S. litura diverged about 5.44 million years ago. Expanded gene families were mainly involved in metabolic detoxification and tolerance to toxic xenobiotics based on GO (Gene Ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway analysis. Comparative genomics analysis showed that gene families involved in detoxification and chemosensation were significantly expanded in S. littoralis, representing genetic characteristics related to polyphagy and an extensive host range.

CONCLUSIONS

We assembled and annotated the reference genome of S. littoralis, and revealed that this pest has the genetic features of strong detoxification capacity, consistent with it being a significant risk to a wide range of host crops. These data resources will provide support for risk assessment and early warning monitoring of major polyphagous agricultural pests.

Contribution of multiple overexpressed carboxylesterase genes to indoxacarb resistance in Spodoptera litura

BACKGROUND

As an important family of detoxification enzymes, carboxylesterases (CarEs) have important roles in the development of insecticide resistance in almost all agricultural pests. Previous studies have suggested that enhancement of CarE activity is an important mechanism mediating indoxacarb resistance in Spodoptera litura, and several CarE genes have been found to be overexpressed in indoxacarb-resistant strains. However, the functions of these CarE genes in indoxacarb resistance needs to be further investigated.

RESULTS

The synergist triphenyl phosphate effectively reduced the resistance of S. litura to indoxacarb, suggesting an involvement of CarEs in indoxacarb resistance. Among seven identified S. litura CarE genes (hereafter SlituCOE), six were overexpressed in two indoxacarb-resistant strains, but there were no significant differences in gene copy number. Knockdown of SlituCOE009 and SlituCOE050 enhanced indoxacarb sensitivity in both susceptible and resistant strains, whereas knockdown of SlituCOE090, SlituCOE093 and SlituCOE074 enhanced indoxacarb sensitivity in only the resistant strain. Knockdown of the sixth gene, SlituCOE073, did not have any effect. Furthermore, simultaneous knockdown of the five SlituCOE genes had a greater effect on increasing indoxacarb sensitivity than silencing them individually. By contrast, overexpression of the five SlituCOE genes individually in Drosophila melanogaster significantly decreased the toxicity of indoxacarb to transgenic fruit flies. Furthermore, modeling and docking analysis indicated that the catalytic pockets of SlituCOE009 and SlituCOE074 were ideally shaped for indoxacarb and N-decarbomethoxylated metabolite (DCJW), but the binding affinity for DCJW was stronger than for indoxacarb.

CONCLUSION

This study reveals that multiple overexpressed CarE genes are involved in indoxacarb resistance in S. litura.

The Genomics and Population Genomics of the Light Brown Apple Moth, Epiphyas postvittana, an Invasive Tortricid Pest of Horticulture

Simple Summary

In this study, we produced a genomic resource for the light brown apple moth, Epiphyas postvittana, to understand the biological basis of adaptation to a high number of hosts (polyphagy) and the invasive nature of this and other lepidopteran pests. The light brown apple moth is an invasive pest of horticultural plants, with over 500 recorded plant hosts. With origins in Australia, the pest has subsequently spread to New Zealand, Hawaii, California and Europe, causing significant economic losses for fruit producers. Comparative genomic analyses with other lepidopteran genomes indicate that a high proportion of the genome is made up of repetitive sequences, with the majority of the known elements being DNA transposable elements and retrotransposons. Twenty gene families show significant expansions, including some likely to have a role in its pest status. Finally, population genomics, investigated by a RAD-tag approach, indicated likely patterns of invasion and admixture, with Californian moths most probably being derived from Australia.

Abstract

The light brown apple moth, Epiphyas postvittana is an invasive, polyphagous pest of horticultural systems around the world. With origins in Australia, the pest has subsequently spread to New Zealand, Hawaii, California and Europe, where it has been found on over 500 plants, including many horticultural crops. We have produced a genomic resource, to understand the biological basis of the polyphagous and invasive nature of this and other lepidopteran pests. The assembled genome sequence encompassed 598 Mb and has an N50 of 301.17 kb, with a BUSCO completion rate of 97.9%. Epiphyas postvittana has 34% of its assembled genome represented as repetitive sequences, with the majority of the known elements made up of longer DNA transposable elements (14.07 Mb) and retrotransposons (LINE 17.83 Mb). Of the 31,389 predicted genes, 28,714 (91.5%) were assigned to 11,438 orthogroups across the Lepidoptera, of which 945 were specific to E. postvittana. Twenty gene families showed significant expansions in E. postvittana, including some likely to have a role in its pest status, such as cytochrome p450s, glutathione-S-transferases and UDP-glucuronosyltransferases. Finally, using a RAD-tag approach, we investigated the population genomics of this pest, looking at its likely patterns of invasion.

The phylogenetic and evolutionary analyses of detoxification gene families in Aphidinae species

Detoxification enzymes play significant roles in the interactions between insects and host plants, wherein detoxification-related genes make great contributions. As herbivorous pests, aphids reproduce rapidly due to parthenogenesis. They are good biological materials for studying the mechanisms that allow insect adaptation to host plants. Insect detoxification gene families are associated with insect adaptation to host plants. The Aphidinae is the largest subfamily in the Aphididae with at least 2483 species in 256 genera in 2 tribes: the Macrosiphini (with 3/4 of the species) and the Aphidini. Most aphid pests on crops and ornamental plants are Aphidinae. Members of the Aphidinae occur in nearly every region of the world. The body shape and colour vary significantly. To research the role that detoxification gene families played in the process of aphid adaptation to host evolution, we analyzed the phylogeny and evolution of these detoxification gene families in Aphidinae. In general, the P450/GST/CCE gene families contract, whereas the ABC/UGT families are conserved in Aphidinae species compared to these families in other herbivorous insects. Genus-specific expansions of P450 CYP4, and GST Delta have occurred in the genus Acyrthosiphon. In addition, the evolutionary rates of five detoxification gene families in the evolution process of Aphidinae are different. The comparison of five detoxification gene families among nine Aphidinae species and the estimated relative evolutionary rates provided herein support an understanding of the interaction between and the co-evolution of Aphidinae and plants.

Population genomic signatures of the oriental fruit moth related to the Pleistocene climates

The Quaternary climatic oscillations are expected to have had strong impacts on the evolution of species. Although legacies of the Quaternary climates on population processes have been widely identified in diverse groups of species, adaptive genetic changes shaped during the Quaternary have been harder to decipher. Here, we assembled a chromosome-level genome of the oriental fruit moth and compared genomic variation among refugial and colonized populations of this species that diverged in the Pleistocene. High genomic diversity was maintained in refugial populations. Demographic analysis showed that the effective population size of refugial populations declined during the penultimate glacial maximum (PGM) but remained stable during the last glacial maximum (LGM), indicating a strong impact of the PGM rather than the LGM on this pest species. Genome scans identified one chromosomal inversion and a mutation of the circadian gene Clk on the neo-Z chromosome potentially related to the endemicity of a refugial population. In the colonized populations, genes in pathways of energy metabolism and wing development showed signatures of selection. These different genomic signatures of refugial and colonized populations point to multiple impacts of Quaternary climates on adaptation in an extant species.

The oriental fruit moth is a pest species native to East Asia with refugial and colonized populations throughout the region. Here, a chromosome-level assembly for the species is reported and used to identify genomic signatures related to Quaternary climate change.

Advances in organophosphorus pesticides pollution: Current status and challenges in ecotoxicological, sustainable agriculture, and degradation strategies

Organophosphorus pesticides (OPPs) are one of the most widely used types of pesticide that play an important role in the production process due to their effects on preventing pathogen infection and increasing yield. However, in the early development and application of OPPs, their toxicological effects and the issue of environmental pollution were not considered. With the long-term overuse of OPPs, their hazards to the ecological environment (including soil and water) and animal health have attracted increasing attention. Therefore, this review first clarified the classification, characteristics, applications of various OPPs, and the government's restriction requirements on various OPPs. Second, the toxicological effects and metabolic mechanisms of OPPs and their metabolites were introduced in organisms. Finally, the existing methods of degrading OPPs were summarized, and the challenges and further addressing strategy of OPPs in the sustainable development of agriculture, the environment, and ecology were prospected. However, methods to solve the environmental and ecological problems caused by OPPs from the three aspects of use source, use process, and degradation methods were proposed, which provided a theoretical basis for addressing the stability of the ecological environment and improving the structure of the pesticide industry in the future.

Identification and dynamic expression profiling of circadian clock genes in Spodoptera litura provide new insights into the regulation of sex pheromone communication

Spodoptera litura is an important pest that causes significant economic damage to numerous crops worldwide. Sex pheromones (SPs) mediate sexual communication in S. litura and show a characteristic degree of rhythmic activity, occurring mainly during the scotophase; however, the specific regulatory mechanisms remain unclear. Here, we employed a genome-wide analysis to identify eight candidate circadian clock genes in S. litura. Sequence characteristics and expression patterns were analyzed. Our results demonstrated that some circadian clock genes might regulate the biosynthesis and perception of SPs by regulating the rhythmic expression of SP biosynthesis-related genes and SP perception-related genes. Interestingly, all potential genes exhibited peak expression in the scotophase, consistent with the SP could mediate courtship and mating behavior in S. litura. Our findings are helpful in elucidating the molecular mechanism by which circadian clock genes regulate sexual communication in S. litura.

Genomic insight into the scale specialization of the biological control agent Novius pumilus (Weise, 1892)

BACKGROUND

Members of the genus Novius Mulsant, 1846 (= Rodolia Mulsant, 1850) (Coleoptera, Coccinellidae), play important roles in the biological control of cotton cushion scale pests, especially those belonging to Icerya. Since the best-known species, the vedalia beetle Novius cardinalis (Mulsant, 1850) was introduced into California from Australia, more than a century of successful use in classical biological control, some species of Novius have begun to exhibit some field adaptations to novel but related prey species. Despite their economic importance, relatively little is known about the underlying genetic adaptations associated with their feeding habits. Knowledge of the genome sequence of Novius is a major step towards further understanding its biology and potential applications in pest control.

RESULTS

We report the first high-quality genome sequence for Novius pumilus (Weise, 1892), a representative specialist of Novius. Computational Analysis of gene Family Evolution (CAFE) analysis showed that several orthogroups encoding chemosensors, digestive, and immunity-related enzymes were significantly expanded (P < 0.05) in N. pumilus compared to the published genomes of other four ladybirds. Furthermore, some of these orthogroups were under significant positive selection pressure (P < 0.05). Notably, transcriptome profiling demonstrated that many genes among the significantly expanded and positively selected orthogroups, as well as genes related to detoxification were differentially expressed, when N. pumilus feeding on the nature prey Icerya compared with the no feeding set. We speculate that these genes are vital in the Icerya adaptation of Novius species.

CONCLUSIONS

We report the first Novius genome thus far. In addition, we provide comprehensive transcriptomic resources for N. pumilus. The results from this study may be helpful for understanding the association of the evolution of genes related to chemosensing, digestion, detoxification and immunity with the prey adaptation of insect predators. This will provide a reference for future research and utilization of Novius in biological control programs. Moreover, understanding the possible molecular mechanisms of prey adaptation also inform mass rearing of N. pumilus and other Novius, which may benefit pest control.

Change of nutrients, microorganisms, and physical properties of exposed extensive green roof substrate

Green roofs (GR) can be used as a nature-based solution to tackle eco-environmental problems caused by climate change and rapid urbanization. The substrate in the GRs is the growing medium for vegetation, and its properties directly affect the ecosystem services of GRs. To investigate the characteristic changes of an exposed substrate after the removal of vegetation, a one-year field experiment was conducted. Substrate properties were comprehensively compared for areas in GR that were planted with Sedum lineare and those with bare substrate. Results show that vegetation cover not only prevented substrate loss by 5.14% (p < 0.05) but also protected the chemical, microbial, and physical properties of the substrate. Moreover, the structure of the substrate changed, as evidenced by a significant increase in fine sand (p < 0.05). The results highlight that attention should be paid to maintaining vegetation cover during GR management. In addition, extensive GRs may not be suitable for fallowing. Once a GR has been established, it needs regular maintenance. Otherwise, the ecological and economic benefits of the GR may be reduced. The findings of the present study can be used to determine the life-cycle costs. Further research should focus on differences in the substrate loss rates, runoff, and temperatures of the substrates under exposure and vegetation cover. The microbial changes after revegetation should also be studied to clarify the role of vegetation in GR ecosystems. The present study provides a reference for improving GR management and ensuring their sustainability.

Genome resequencing reveals rapid, repeated evolution in the Colorado potato beetle

Insecticide resistance and rapid pest evolution threatens food security and the development of sustainable agricultural practices, yet the evolutionary mechanisms that allow pests to rapidly adapt to control tactics remains unclear. Here we examine how a global super-pest, the Colorado potato beetle (CPB), Leptinotarsa decemlineata, rapidly evolves resistance to insecticides. Using whole genome resequencing and transcriptomic data focused on its ancestral and pest range in North America, we assess evidence for three, non-mutually exclusive models of rapid evolution: pervasive selection on novel mutations, rapid regulatory evolution, and repeated selection on standing genetic variation. Population genomic analysis demonstrates that CPB is geographically structured, even among recently established pest populations. Pest populations exhibit similar levels of nucleotide diversity, relative to non-pest populations, and show evidence of recent expansion. Genome scans provide clear signatures of repeated adaptation across CPB populations, with especially strong evidence of selection on insecticide resistance genes in different populations. Analyses of gene expression show that constitutive upregulation of candidate insecticide resistance genes drives distinctive population patterns. CPB evolves insecticide resistance repeatedly across agricultural regions, leveraging similar genetic pathways but different genes, demonstrating a polygenic trait architecture for insecticide resistance that can evolve from standing genetic variation. Despite expectations, we do not find support for strong selection on novel mutations, or rapid evolution from selection on regulatory genes. These results suggest that integrated pest management practices must mitigate the evolution of polygenic resistance phenotypes among local pest populations, in order to maintain the efficacy and sustainability of novel control techniques.

Bombyx mori-derived aldo-keto reductase AKR2E8 detoxifies aldehydes present in mulberry leaves

Lepidopterans are agricultural pests. Since the silkworm is a model for lepidopterans, analysis of the enzymes produced by silkworms is of great interest for developing methods of pest control. The aldo-keto reductase (AKR) superfamily catalyzes the reduction of aldehydes by converting a carbonyl group to an alcohol group. Here, we characterized a new AKR present in the silkworm Bombyx mori, which has been designated as AKR2E8. Amino acid sequence and phylogenetic analyses showed that AKR2E8 is similar to human AKR1B1 and AKR1B10. Three amino acid residues in the active site were identical among AKR2E8, AKR1B1, and AKR1B10. Recombinant AKR2E8 overexpressed in Escherichia coli used nicotinamide adenine dinucleotide phosphate as a coenzyme to reduce the aldehydes present in mulberry (Morus alba) leaves. AKR2E8 was found to reduce benzaldehyde, hexanal, heptanal, nonanal, trans-2-nonenal, and citral. No nicotinamide adenine dinucleotide-dependent activity was detected. Akr2e8 mRNA was detected in the testes, ovaries, and fat body; the highest expression was found in the midgut. The substrate specificity and highest observed expression of AKR2E8 in the midgut suggests that AKR2E8 may play a major role in aldehyde detoxification in silkworms. The findings of this study may assist in the development of pest control methods for controlling the population of lepidopterans, such as silkworms, that damage crops.

Expanding the Menu: Are Polyphagy and Gene Family Expansions Linked across Lepidoptera?

Evolutionary expansions and contractions of gene families are often correlated with key innovations and/or ecological characteristics. In butterflies and moths (Lepidoptera), expansions of gene families involved in detoxification of plant specialized metabolites are hypothesized to facilitate a polyphagous feeding style. However, analyses supporting this hypothesis are mostly based on a limited number of lepidopteran species. We applied a phylogenomics approach, using 37 lepidopteran genomes, to analyze if gene family evolution (gene gain and loss) is associated with the evolution of polyphagy. Specifically, we compared gene counts and evolutionary gene gain and loss rates of gene families involved in adaptations with plant feeding. We correlated gene evolution to host plant family range (phylogenetic diversity) and specialized metabolite content of plant families (functional metabolite diversity). We found a higher rate for gene loss than gene gain in Lepidoptera, a potential consequence of genomic rearrangements and deletions after (potentially small-scale) duplication events. Gene family expansions and contractions varied across lepidopteran families, and were associated to host plant use and specialization levels. Within the family Noctuidae, a higher expansion rate for gene families involved in detoxification can be related to the large number of polyphagous species. However, gene family expansions are observed in both polyphagous and monophagous lepidopteran species and thus seem to be species-specific in the taxa sampled. Nevertheless, a significant positive correlation of gene counts of the carboxyl- and choline esterase and glutathione-S-transferase detoxification gene families with the level of polyphagy was identified across Lepidoptera.

Spodoptera frugiperda Sf9 cells as a model system to investigate the role of detoxification gene expression in response to xenobiotics

Spodoptera frugiperda (fall armyworm) is a highly destructive invasive pest that feeds on numerous crops including maize and rice. It has developed sophisticated mechanisms to detoxify xenobiotics such as secondary plant metabolites as well as manmade insecticides. The aim of the study was to explore the detoxification response to plant secondary metabolites and insecticides employing a S. frugiperda Sf9 cell model exposed to indole 3-carbinol (I3C) and methoprene. The cell Inhibitory Concentration 50 (IC₅₀) for these molecules was determined and IC₁₀, IC₂₀ and IC₃₀ doses were used to monitor the induction profiles of detoxification genes. Cytochrome P450 monooxygenases (P450s) of the CYP9A subfamily were the most inducible genes of the seven examined. Our results also showed the induction of the transcription factor Cap'n'collar isoform C (CncC). Transient transformation of Sf9 cells overexpressing CncC and its partner muscle aponeurosis fibromatosis (Maf) induces overexpression of CYP4M14, CYP4M15, CYP321A9 and GSTE1 while CYP9As were not induced. Next, we determined the capacity of recombinantly expressed CYP9A30, CYP9A31 and CYP9A32 to interact with methoprene and I3C. Fluorescence-based biochemical assays revealed an interaction of methoprene with functionally expressed CYP9A30, CYP9A31 and CYP9A32 whereas almost no interaction was detected for I3C, suggesting the ability of CYP9As to metabolize methoprene. Our results showed that Sf9 cells could be a useful model to decipher detoxification pathways of S. frugiperda.

Parallel Evolution of C-Type Lectin Domain Gene Family Sizes in Insect-Vectored Nematodes

The dispersal stage of pathogens is crucial for the successful spread and infection of their hosts. Some plant-parasitic nematodes (PPNs) have evolved specialized dispersal stages to reach healthy hosts by being carried out by insect vectors. Because gene gain and loss is a major factor contributing to the evolution of novel characteristics, it is essential to clarify the gene family characteristics among nematodes with different dispersal modes to disentangle the evolution of insect-mediated dispersal. Here, the size of the C-type lectin (CTL) family genes of insect-vectored nematodes was found to be drastically reduced compared with those of self-dispersing nematodes, whereas the diversity of their functional domains was significantly higher. The gene family sizes of vector-dispersed nematodes were only a twentieth of the size of that of a self-dispersing (i.e., without a biotic vector) nematode model Caenorhabditis elegans, and these genes were inactive during the dispersal stage. Phylogenetic analysis showed that some CTL genes of vector-borne PPNs shared higher homology to the animal parasitic nematodes compared with other PPNs. Moreover, homology modeling predicted that the CTLs of insect-vectored nematodes bear remarkable structural similarity to the lectin genes of their vector's immune system. Because CTL genes are important sugar-binding proteins for the innate immune response of C. elegans, the loss of some CTL genes of vector-transmitted PPNs might be responsible for their parallel adaptations to a mutualistic relationship with their vector. These results expand our understanding of the evolutionary benefits of vector-mediated transmission for the nematode and vector-nematode co-evolution.

Silencing of cytochrome P450 gene CYP321A1 effects tannin detoxification and metabolism in Spodoptera litura

Cytochrome P450 monooxygenase (P450 or CYP) plays an important role in the metabolism of insecticides and plant allelochemicals by insects. CYP321B1, a novel Spodoptera litura P450 gene, was identified and characterized. CYP321B1 contains a 1488 bp open reading frame (ORF) that encodes a 495 amino acid protein. In fourth instar larvae, the highest CYP321B1 expression levels were found in the midgut and fat body. In the tannin feeding test, tannin can significantly induce the expression of CYP321B1 in the midgut and fat body of 4th instar larvae. To verify the function of CYP321B1, RNA interference and metabolome analysis were performed. The results showed that silencing CYP321B1 significantly reduced the rate of weight gain under tannin induction. Metabolome analysis showed silencing affected 47 different metabolites, mainly involved in secondary metabolite biosynthesis and amino acid metabolism, including amino acids, lipid fatty acids, organic acids and their derivatives. Henoxyacetic acid and cysteamine are the most highly regulated metabolites, respectively. These findings demonstrate that CYP321B1 plays an important role in tannin detoxification and metabolism. Functional knowledge about metabolite detoxification genes in this major herbivorous insect pest can provide new insights into this biological process and provide new targets for agricultural pest control.