A model species for agricultural pest genomics: the genome of the Colorado potato beetle, Leptinotarsa decemlineata (Coleoptera: Chrysomelidae)

Role of CYP9E2 and a long non-coding RNA gene in resistance to a spinosad insecticide in the Colorado potato beetle, Leptinotarsa decemlineata

Spinosads are insecticides used to control insect pests, especially in organic farming where limited tools for pest management exist. However, resistance has developed to spinosads in economically important pests, including Colorado potato beetle (CPB), Leptinotarsa decemlineata. In this study, we used bioassays to determine spinosad sensitivity of two field populations of CPB, one from an organic farm exposed exclusively to spinosad and one from a conventional farm exposed to a variety of insecticides, and a reference insecticide naïve population. We found the field populations exhibited significant levels of resistance compared with the sensitive population. Then, we compared transcriptome profiles between the two field populations to identify genes associated primarily with spinosad resistance and found a cytochrome P450, CYP9E2, and a long non-coding RNA gene, lncRNA-2, were upregulated in the exclusively spinosad-exposed population. Knock-down of these two genes simultaneously in beetles of the spinosad-exposed population using RNA interference (RNAi) resulted in a significant increase in mortality when gene knock-down was followed by spinosad exposure, whereas single knock-downs of each gene produced smaller effects. In addition, knock-down of the lncRNA-2 gene individually resulted in significant reduction in CYP9E2 transcripts. Finally, in silico analysis using an RNA-RNA interaction tool revealed that CYP9E2 mRNA contains multiple binding sites for the lncRNA-2 transcript. Our results imply that CYP9E2 and lncRNA-2 jointly contribute to spinosad resistance in CPB, and lncRNA-2 is involved in regulation of CYP9E2 expression. These results provide evidence that metabolic resistance, driven by overexpression of CYP and lncRNA genes, contributes to spinosad resistance in CPB.

Protein complexes from edible mushrooms as a sustainable potato protection against coleopteran pests

Protein complexes from edible oyster mushrooms (Pleurotus sp.) composed of pleurotolysin A2 (PlyA2) and pleurotolysin B (PlyB) exert toxicity in feeding tests against Colorado potato beetle (CPB) larvae, acting through the interaction with insect-specific membrane sphingolipid. Here we present a new strategy for crop protection, based on in planta production of PlyA2/PlyB protein complexes, and we exemplify this strategy in construction of transgenic potato plants of cv Désirée. The transgenics in which PlyA2 was directed to the vacuole and PlyB to the endoplasmic reticulum are effectively protected from infestation by CPB larvae without impacting plant performance. These transgenic plants showed a pronounced effect on larval feeding rate, the larvae feeding on transgenic plants being on average five to six folds lighter than larvae feeding on controls. Further, only a fraction (11%-37%) of the larvae that fed on transgenic potato plants completed their life cycle and developed into adult beetles. Moreover, gene expression analysis of CPB larvae exposed to PlyA2/PlyB complexes revealed the response indicative of a general stress status of larvae and no evidence of possibility of developing resistance due to the functional inactivation of PlyA2/PlyB sphingolipid receptors.

Highly dynamic evolution of the chemosensory system driven by gene gain and loss across subterranean beetles

Chemical cues in subterranean habitats differ highly from those on the surface due to the contrasting environmental conditions, such as absolute darkness, high humidity or food scarcity. Subterranean animals underwent changes to their sensory systems to facilitate the perception of essential stimuli for underground lifestyles. Despite representing unique systems to understand biological adaptation, the genomic basis of chemosensation across cave-dwelling species remains unexplored from a macroevolutionary perspective. Here, we explore the evolution of chemoreception in three beetle tribes that underwent at least six independent transitions to the underground, through a phylogenomics spyglass. Our findings suggest that the chemosensory gene repertoire varies dramatically between species. Overall, no parallel changes in the net rate of evolution of chemosensory gene families were detected prior, during, or after the habitat shift among subterranean lineages. Contrarily, we found evidence of lineage-specific changes within surface and subterranean lineages. However, our results reveal key duplications and losses shared between some of the lineages transitioning to the underground, including the loss of sugar receptors and gene duplications of the highly conserved ionotropic receptors IR25a and IR8a, involved in thermal and humidity sensing among other olfactory roles in insects. These duplications were detected both in independent subterranean lineages and their surface relatives, suggesting parallel evolution of these genes across lineages giving rise to cave-dwelling species. Overall, our results shed light on the genomic basis of chemoreception in subterranean beetles and contribute to our understanding of the genomic underpinnings of adaptation to the subterranean lifestyle at a macroevolutionary scale.

Crucial roles of specialized chitinases in elytral and hindwing cuticles construction in Leptinotarsa decemlineata

BACKGROUND

The Colorado potato beetle (CPB), Leptinotarsa decemlineata, is a major potato (Solanum tuberosum) pest, infesting over 16 million km2 and causing substantial economic losses. The insect cuticle forms an apical extracellular matrix (ECM) envelope covering exposed organs to direct morphogenesis and confer structural protection. While select chitinase (Cht) genes have proven essential for larval development, their potential activities directing ECM remodeling underlying adult wing maturation remain undefined.

RESULTS

We investigated the expression patterns and performed an oral RNA interference (RNAi) screen targeting 19 LdChts in late-instar L. decemlineata larvae. Subsequently, we assessed their effects on adult eclosion and wing characteristics. Knockdown of LdCht5, LdCht7, LdCht10, LdIDGF2, and LdIDGF4, as well as others from Group IV (LdCht15, LdCht12, LdCht17, and LdCht13) and Groups VII-X (LdCht2, LdCht11, LdCht1, and LdCht3), resulting in shrunken, misshapen elytra with reduced areal density, as well as transverse wrinkling and impaired wing-tip folding in hindwings. Scanning electron micrographs revealed eroded elytral ridges alongside thinned, ruptured hindwing veins, indicative of mechanical fragility post-LdCht suppression. Spectroscopic analysis uncovered biomolecular alterations underlying the elytral anomalies, including decreases in peaks representing chitin, proteins, and lipids. This loss of essential ECM components provides evidence for the fragility, wrinkling, and shrinkage observed in the RNAi groups.

CONCLUSION

Our findings elucidate the crucial role of chitinases in the turnover of chitinous cuticles on beetle wings, offering insights into RNAi-based control strategies against this invasive pest. © 2024 Society of Chemical Industry.

A near-complete genome assembly of Monochamus alternatus a major vector beetle of pinewood nematode

The Japanese sawyer beetle, Monochamus alternatus, is not only one of the most important wood boring pest itself, but also a major vector of the invasive pinewood nematode (PWN), which is the causal agent of the devastative pine wilt disease (PWD) and threats the global pine forest. Here, we present a near-complete genome of M. alternatus at the chromosome level. The assembled genome was 792.05 Mb with contig N50 length of 55.99 Mb, which is the largest N50 size among the sequenced Coleoptera insects currently. 99.57% of sequence was anchored onto ten pseudochromosomes (one X-chromosome and nine autosomes), and the final genome harbored only 13 gaps. BUSCO evaluation revealed the presence of 99.0% of complete core genes. Thus, our genome assembly represented the highest-contiguity genome assembly as well as high completeness in insects so far. We identified 20,471 protein-coding genes, of which 20,070 (98.04%) were functionally annotated. The genome assembly of M. alternatus provides a valuable resource for exploring the evolution of the symbiosis between PWN and the vector insects.

De novo assembly and annotation of Popillia japonica's genome with initial clues to its potential as an invasive pest

BACKGROUND

The spread of Popillia japonica in non-native areas (USA, Canada, the Azores islands, Italy and Switzerland) poses a significant threat to agriculture and horticulture, as well as to endemic floral biodiversity, entailing that appropriate control measures must be taken to reduce its density and limit its further spread. In this context, the availability of a high quality genomic sequence for the species is liable to foster basic research on the ecology and evolution of the species, as well as on possible biotechnologically-oriented and genetically-informed control measures.

RESULTS

The genomic sequence presented and described here is an improvement with respect to the available draft sequence in terms of completeness and contiguity, and includes structural and functional annotations. A comparative analysis of gene families of interest, related to the species ecology and potential for polyphagy and adaptability, revealed a contraction of gustatory receptor genes and a paralogous expansion of some subgroups/subfamilies of odorant receptors, ionotropic receptors and cytochrome P450s.

CONCLUSIONS

The new genomic sequence as well as the comparative analyses data may provide a clue to explain the staggering invasive potential of the species and may serve to identify targets for potential biotechnological applications aimed at its control.

Genome evolution and divergence in cis-regulatory architecture is associated with condition-responsive development in horned dung beetles

Phenotypic plasticity is thought to be an important driver of diversification and adaptation to environmental variation, yet the genomic mechanisms mediating plastic trait development and evolution remain poorly understood. The Scarabaeinae, or true dung beetles, are a species-rich clade of insects recognized for their highly diversified nutrition-responsive development including that of cephalic horns-evolutionarily novel, secondary sexual weapons that exhibit remarkable intra- and interspecific variation. Here, we investigate the evolutionary basis for horns as well as other key dung beetle traits via comparative genomic and developmental assays. We begin by presenting chromosome-level genome assemblies of three dung beetle species in the tribe Onthophagini (> 2500 extant species) including Onthophagus taurus, O. sagittarius, and Digitonthophagus gazella. Comparing these assemblies to those of seven other species across the order Coleoptera identifies evolutionary changes in coding sequence associated with metabolic regulation of plasticity and metamorphosis. We then contrast chromatin accessibility in developing head horn tissues of high- and low-nutrition O. taurus males and females and identify distinct cis-regulatory architectures underlying nutrition- compared to sex-responsive development, including a large proportion of recently evolved regulatory elements sensitive to horn morph determination. Binding motifs of known and new candidate transcription factors are enriched in these nutrition-responsive open chromatin regions. Our work highlights the importance of chromatin state regulation in mediating the development and evolution of plastic traits, demonstrates gene networks are highly evolvable transducers of environmental and genetic signals, and provides new reference-quality genomes for three species that will bolster future developmental, ecological, and evolutionary studies of this insect group.

Identification of chitinase genes and roles in the larval-pupal transition of Leptinotarsa decemlineata

BACKGROUND

Insect chitinases play crucial roles in degrading chitin in the extracellular matrix, affecting insect development and molting. However, our understanding of the specific functions of various chitinases in Leptinotarsa decemlineata is limited, hindering the deployment of novel gene-targeting technologies as pest management strategies.

RESULTS

We identified and characterized 19 full-length complementary DNA (cDNA) sequences of chitinase genes (LdChts) in Leptinotarsa decemlineata. Despite having varying domain architectures, all these chitinases contained at least one chitinase catalytic domain. Phylogenetic analysis classified the chitinase proteins into ten distinct clusters (groups I-X). Expression profiles showed the highest expression in chitin-rich tissues or during specific developmental stages from the larva-to-pupa transition. Gene-specific RNA interference (RNAi) experiments provided valuable insight into chitinase gene function. Silencing of group II LdCht10 prevented larval-larval molting, larval-prepupal, and prepupal-pupal processes. Moreover, our study revealed that LdCht5, LdCht2, LdCht11, LdCht1, and LdCht3 from groups I and VII-X were specifically essential for the transition from prepupal to pupal stage, whereas LdIDGF2 from group V was necessary for the larval-prepupal metamorphic process. The chitinase gene LdCht7 from group III and LdIDGF4 from group V were involved in both the larva-to-prepupa and the prepupa-to-pupa shift. Additionally, our findings also shed light on the exclusive expression of nine chitinase genes within group IV in the digestive system, suggesting their potential role in regulating larval body weight and larva-to-pupa transition.

CONCLUSION

Our results provide a comprehensive understanding of the functional specialization of chitinase genes during the molting process of various stages and identify potential targets for RNAi-based management of Leptinotarsa decemlineata. © 2023 Society of Chemical Industry.

Investigation of genetic diversity and polyandry of Leptinotarsa decemlineata using X-linked microsatellite markers

A panel of X-linked microsatellite markers was newly designed using the data from a previous sequencing project available in NCBI and used for a study of the Colorado potato beetle (CPB, Leptinotarsa decemlineata) X-haplotype variability. The analysis of scaffolds 49 and 61 (newly identified as fragments of CPB chromosome X) found ten high-quality markers, which were arranged in two PCR multiplexes and evaluated in both 420 CPB adults, collected from 14 localities of Czechia and Slovakia, and 866 larvae from five single-female families from two more Czech localities. Length polymorphisms found in 6 loci have predicted 192 potential X-haplotypes, however, only 36 combinations were detected in the adult males (N = 189), and seven additional ones in the larvae. The X-haplotypes were also generally unevenly distributed; five of the most frequent haplotypes were detected in 55% of males, 19 repeating up to ten-times in 38.7% of males and the remained 12 occurred uniquely in 6.3% of males. Bulk analysis of X-haplotypes dissimilarity indicated seven haplotype groups diversified by mutations and recombinations. Two haplotypes showed a distinctive regional distribution, which indicates an east-west disruption of CPB migration probably caused by different environments of localities in the South Bohemia region and Vysocina region. On the contrary, the results indicate a south-north migration corridor alongside the Vltava River. In the single-female families, from 6 to 13 distinct paternal haplotypes were detected, which proved and quantified a frequented polyandry in CPB.

Roles of insect odorant binding proteins in communication and xenobiotic adaptation

Odorant binding proteins (OBPs) are small water-soluble proteins mainly associated with olfaction, facilitating the transport of odorant molecules to their relevant receptors in the sensillum lymph. While traditionally considered essential for olfaction, recent research has revealed that OBPs are engaged in a diverse range of physiological functions in modulating chemical communication and defense. Over the past 10 years, emerging evidence suggests that OBPs play vital roles in purifying the perireceptor space from unwanted xenobiotics including plant volatiles and pesticides, potentially facilitating xenobiotic adaptation, such as host location, adaptation, and pesticide resistance. This multifunctionality can be attributed, in part, to their structural variability and effectiveness in transporting, sequestering, and concealing numerous hydrophobic molecules. Here, we firstly overviewed the classification and structural properties of OBPs in diverse insect orders. Subsequently, we discussed the myriad of functional roles of insect OBPs in communication and their adaptation to xenobiotics. By synthesizing the current knowledge in this field, our review paper contributes to a comprehensive understanding of the significance of insect OBPs in chemical ecology, xenobiotic adaptation, paving the way for future research in this fascinating area of study.

Lack of fitness costs in dsRNA-resistant Leptinotarsa decemlineata ([Coleoptera]: [Chrysomelidae])

The Colorado potato beetle, Leptinotarsa decemlineata (Say) ([Coleoptera]: [Chrysomelidae]), is the most important defoliator of solanaceous plants worldwide. This insect displays a notorious ability in adapting to biological and synthetic insecticides, although in some cases this adaptation carries relevant fitness costs. Insecticidal gene silencing by RNA interference is a novel mode of action pesticide against L. decemlineata that is activated by ingestion of a double stranded RNA (dsRNA) targeting a vital L. decemlineata gene. We previously reported laboratory selection of a > 11,000-fold resistant strain of L. decemlineata to a dsRNA delivered topically to potato leaves. In this work, we tested the existence of fitness costs in this dsRNA-resistant colony by comparing biological parameters to the parental strain and an additional susceptible reference strain. Biological parameters included length of egg incubation period, number of eggs per clutch, egg viability, larval viability, length of larval and pupal periods, adult emergence, number of eggs laid per day, sex ratio, and adult longevity. Comparisons between the 3 beetle strains detected no fitness costs associated with resistance to dsRNA. This information is important to guide effective insect resistance management plans for dsRNA insecticides against L. decemlineata applied topically to potato leaves.

Evolution of chemosensory and detoxification gene families across herbivorous Drosophilidae

Herbivorous insects are exceptionally diverse, accounting for a quarter of all known eukaryotic species, but the genomic basis of adaptations that enabled this dietary transition remains poorly understood. Many studies have suggested that expansions and contractions of chemosensory and detoxification gene families-genes directly mediating interactions with plant chemical defenses-underlie successful plant colonization. However, this hypothesis has been challenging to test because the origins of herbivory in many insect lineages are ancient (>150 million years ago (mya)), obscuring genomic evolutionary patterns. Here, we characterized chemosensory and detoxification gene family evolution across Scaptomyza, a genus nested within Drosophila that includes a recently derived (<15 mya) herbivore lineage of mustard (Brassicales) specialists and carnation (Caryophyllaceae) specialists, and several nonherbivorous species. Comparative genomic analyses revealed that herbivorous Scaptomyza has among the smallest chemosensory and detoxification gene repertoires across 12 drosophilid species surveyed. Rates of gene turnover averaged across the herbivore clade were significantly higher than background rates in over half of the surveyed gene families. However, gene turnover was more limited along the ancestral herbivore branch, with only gustatory receptors and odorant-binding proteins experiencing strong losses. The genes most significantly impacted by gene loss, duplication, or changes in selective constraint were those involved in detecting compounds associated with feeding on living plants (bitter or electrophilic phytotoxins) or their ancestral diet (fermenting plant volatiles). These results provide insight into the molecular and evolutionary mechanisms of plant-feeding adaptations and highlight gene candidates that have also been linked to other dietary transitions in Drosophila.

Global patterns of genomic and phenotypic variation in the invasive harlequin ladybird

BACKGROUND

The harlequin ladybird Harmonia axyridis (Coleoptera: Coccinellidae), native to Asia, has been introduced to other major continents where it has caused serious negative impacts on local biodiversity. Though notable advances to understand its invasion success have been made during the past decade, especially with then newer molecular tools, the conclusions reached remain to be confirmed with more advanced genomic analyses and especially using more samples from larger geographical regions across the native range. Furthermore, although H. axyridis is one of the best studied invasive insect species with respect to life history traits (often comparing invasive and native populations), the traits responsible for its colonization success in non-native areas warrant more research.

RESULTS

Our analyses of genome-wide nuclear population structure indicated that an eastern Chinese population could be the source of all non-native populations and revealed several putatively adaptive candidate genomic loci involved in body color variation, visual perception, and hemolymph synthesis. Our estimates of evolutionary history indicate (1) asymmetric migration with varying population sizes across its native and non-native range, (2) a recent admixture between eastern Chinese and American populations in Europe, (3) signatures of a large progressive, historical bottleneck in the common ancestors of both populations and smaller effective sizes of the non-native population, and (4) the southwest origin and subsequent dispersal routes within its native range in China. In addition, we found that while two mitochondrial haplotypes-Hap1 and Hap2 were dominant in the native range, Hap1 was the only dominant haplotype in the non-native range. Our laboratory observations in both China and USA found statistical yet slight differences between Hap1 and Hap2 in some of life history traits.

CONCLUSIONS

Our study on H. axyridis provides new insights into its invasion processes into other major continents from its native Asian range, reconstructs a geographic range evolution across its native region China, and tentatively suggests that its invasiveness may differ between mitochondrial haplotypes.

Tracking Adaptive Pathways of Invasive Insects: Novel Insight from Genomics

Despite the huge human and economic costs of invasive insects, which are the main group of invasive species, their environmental impacts through various mechanisms remain inadequately explained in databases and much of the invasion biology literature. High-throughput sequencing technology, especially whole-genome sequencing, has been used as a powerful method to study the mechanisms through which insects achieve invasion. In this study, we reviewed whole-genome sequencing-based advances in revealing several important invasion mechanisms of invasive insects, including (1) the rapid genetic variation and evolution of invasive populations, (2) invasion history and dispersal paths, (3) rapid adaptation to different host plant ranges, (4) strong environmental adaptation, (5) the development of insecticide resistance, and (6) the synergistic damage caused by invasive insects and endosymbiotic bacteria. We also discussed prevention and control technologies based on whole-genome sequencing and their prospects.

Chromosome-level genome assembly of Dastarcus helophoroides provides insights into CYP450 genes expression upon insecticide exposure

BACKGROUND

Dastarcus helophoroides is an important natural enemy of cerambycids, and is wildly used in biological control of pests. Nevertheless, the absence of complete genomic information limits the investigation of the underlying molecular mechanisms. Here, a chromosome-level of Dastarcus helophoroides genome is assembled using a combination strategy of Illumina, PacBio, 10x™ Genomics, and Hi-C.

RESULTS

The final assembly is 609.09 Mb with contig N50, scaffold N50 and GC content of 5.46 Mb, 42.56 Mb and 31.50%, respectively, and 95.25% of the contigs anchor into 13 chromosomes. In total 14 890 protein-coding genes and 65.37% repeat sequences are predicted in the assembly genome. The phylogenetic analysis of single-copy gene families shared among 20 insect species indicates that Dastarcus helophoroides is placed as the sister species to clade (Nitidulidae+Curculionoidea+Chrysomeloidea) + Tenebrionoidea, and diverges from the related species ~242.9 Mya. In total 36 expanded gene families are identified in Dastarcus helophoroides genome, and are functionally related to drug metabolism and metabolism of xenobiotics by cytochrome P450. Some members of CYP4 Clade and CYP6 Clade are up-regulated in Dastarcus helophoroides adults upon insecticide exposure, of which expressions of DhCYP4Q, DhCYP6A14X1 and DhCYP4C1 are significantly up-regulated. The silencing of the three genes leads to adults more sensitive to insecticide and increased knocked-down rate, which may indicate their critical roles in stress resistance and detoxication.

CONCLUSION

Our study systematically integrated the chromosome-level genome, transcriptome and gene expression of Dastarcus helophoroides, which will provide valuable resources for understanding mechanisms of pesticide metabolism, growth and development, and utilization of the natural enemy in integrated control. © 2022 Society of Chemical Industry.

The effects of short-term glyphosate-based herbicide exposure on insect gene expression profiles

Glyphosate-based herbicides (GBHs) are the most frequently used herbicides worldwide. The use of GBHs is intended to tackle weeds, but GBHs have been shown to affect the life-history traits and antioxidant defense system of invertebrates found in agroecosystems. Thus far, the effects of GBHs on detoxification pathways among invertebrates have not been sufficiently investigated. We performed two different experiments-1) the direct pure glyphosate and GBH treatment, and 2) the indirect GBH experiment via food-to examine the possible effects of environmentally relevant GBH levels on the survival of the Colorado potato beetle (Leptinotarsa decemlineata) and the expression profiles of their detoxification genes. As candidate genes, we selected four cytochrome P450 (CYP), three glutathione-S-transferase (GST), and two acetylcholinesterase (AChE) genes that are known to be related to metabolic or target-site resistances in insects. We showed that environmentally relevant levels of pure glyphosate and GBH increased the probability for higher mortality in the Colorado potato beetle larvae in the direct experiment, but not in the indirect experiment. The GBHs or glyphosate did not affect the expression profiles of the studied CYP, GST, or AChE genes; however, we found a large family-level variation in expression profiles in both the direct and indirect treatment experiments. These results suggest that the genes selected for this study may not be the ones expressed in response to glyphosate or GBHs. It is also possible that the relatively short exposure time did not affect gene expression profiles, or the response may have already occurred at a shorter exposure time. Our results show that glyphosate products may affect the survival of the herbivorous insect already at lower levels, depending on their sensitivity to pesticides.

Identification and expression analyses of the olfactory-related genes in different tissues' transcriptome of a predacious soldier beetle, Podabrus annulatus (Coleoptera, Cantharidae)

We sequenced and analyzed the transcriptomes from different tissues of the soldier beetle, Podabrus annulatus (Coleoptera: Cantharidae), and obtained 75.74 Gb clean reads which were assembled into 95,274 unigenes. Among these transcripts, 25,484 unigenes of highly quality were annotated. Based on annotation and tBLASTn results, we identified a total of 101 candidate olfactory-related genes for the first time, including 11 putative odorant-binding proteins (OBPs), 6 chemosensory proteins (CSP), 50 olfactory receptors (ORs), 25 gustatory receptors (GRs), 6 ionotropic receptors (IRs), and 3 sensory neuron membrane proteins (SNMPs). BLASTX best-hit results indicated that these chemosensory genes were most identical to their respective orthologs from Photinus pyralis. Phylogenetic analyses also revealed that the ORs, GRs, and IRs of Podabrus annulatus are closely related to those of Photinus pyralis. The fragment per kilobase per million mapped fragments (FPKM) values showed that the PannOBP2, PannOBP3, and PannOBP10 were predominantly expressed in the antennae, PannOBP1 in the abdomen-thorax, while others were not identified to be tissue-specific. These olfactory-related differentially expressed genes (DEGs) demonstrated different roles in the olfactory system of Podabrus annulatus. This study establishes the groundwork for future research into the molecular mechanism of olfactory recognition in Podabrus annulatus.

Evolution of chemosensory and detoxification gene families across herbivorous Drosophilidae

Herbivorous insects are exceptionally diverse, accounting for a quarter of all known eukaryotic species, but the genetic basis of adaptations that enabled this dietary transition remains poorly understood. Many studies have suggested that expansions and contractions of chemosensory and detoxification gene families - genes directly mediating interactions with plant chemical defenses - underlie successful plant colonization. However, this hypothesis has been challenging to test because the origins of herbivory in many lineages are ancient (>150 million years ago [mya]), obscuring genomic evolutionary patterns. Here, we characterized chemosensory and detoxification gene family evolution across Scaptomyza, a genus nested within Drosophila that includes a recently derived (<15 mya) herbivore lineage of mustard (Brassicales) specialists and carnation (Caryophyllaceae) specialists, and several non-herbivorous species. Comparative genomic analyses revealed that herbivorous Scaptomyza have among the smallest chemosensory and detoxification gene repertoires across 12 drosophilid species surveyed. Rates of gene turnover averaged across the herbivore clade were significantly higher than background rates in over half of the surveyed gene families. However, gene turnover was more limited along the ancestral herbivore branch, with only gustatory receptors and odorant binding proteins experiencing strong losses. The genes most significantly impacted by gene loss, duplication, or changes in selective constraint were those involved in detecting compounds associated with feeding on plants (bitter or electrophilic phytotoxins) or their ancestral diet (yeast and fruit volatiles). These results provide insight into the molecular and evolutionary mechanisms of plant-feeding adaptations and highlight strong gene candidates that have also been linked to other dietary transitions in Drosophila .

Spinosad-associated modulation of select cytochrome P450s and glutathione S-transferases in the Colorado potato beetle, Leptinotarsa decemlineata

The Colorado potato beetle (Leptinotarsa decemlineata (Say)) is an insect pest that threatens potato crops. Multiple options exist to limit the impact of this pest even though insecticides remain a primary option for its control. Insecticide resistance has been reported in Colorado potato beetles and a better understanding of the molecular players underlying such process is of utmost importance to optimize the tools used to mitigate the impact of this insect. Resistance against the insecticide spinosad has been reported in this insect and this work thus aims at exploring the expression of targets previously associated with insecticide response in Colorado potato beetles exposed to this compound. Amplification and quantification of transcripts coding for cytochrome P450s and glutathione S-transferases were conducted via qRT-PCR in insects treated with varying doses of spinosad and for different time duration. This approach notably revealed differential expression of CYP6a23 and CYP12a5 in insects exposed to low doses of spinosad for 4 h as well as modulation of CYP6a13, CYP6d4, GST, GST1, and GST1-Like in insects treated with high doses of spinosad for the same duration. RNAi-based targeting of CYP4g15 and CYP6a23 was associated with marked reduction of transcript expression 7 days following dsRNA injection and reduction of the former had a marked impact on insect viability. In general, results presented here provide novel information regarding the expression of transcripts relevant to spinosad response in Colorado potato beetles and reveal a novel target to consider in the development of RNAi-based strategies aimed at this potato pest.

Role of Insect and Mammal Glutathione Transferases in Chemoperception

Glutathione transferases (GSTs) are ubiquitous key enzymes with different activities as transferases or isomerases. As key detoxifying enzymes, GSTs are expressed in the chemosensory organs. They fulfill an essential protective role because the chemosensory organs are located in the main entry paths of exogenous compounds within the body. In addition to this protective function, they modulate the perception process by metabolizing exogenous molecules, including tastants and odorants. Chemosensory detection involves the interaction of chemosensory molecules with receptors. GST contributes to signal termination by metabolizing these molecules. By reducing the concentration of chemosensory molecules before receptor binding, GST modulates receptor activation and, therefore, the perception of these molecules. The balance of chemoperception by GSTs has been shown in insects as well as in mammals, although their chemosensory systems are not evolutionarily connected. This review will provide knowledge supporting the involvement of GSTs in chemoperception, describing their localization in these systems as well as their enzymatic capacity toward odorants, sapid molecules, and pheromones in insects and mammals. Their different roles in chemosensory organs will be discussed in light of the evolutionary advantage of the coupling of the detoxification system and chemosensory system through GSTs.

The redlegged earth mite draft genome provides new insights into pesticide resistance evolution and demography in its invasive Australian range

Genomic data provide valuable insights into pest management issues such as resistance evolution, historical patterns of pest invasions and ongoing population dynamics. We assembled the first reference genome for the redlegged earth mite, Halotydeus destructor (Tucker, 1925), to investigate adaptation to pesticide pressures and demography in its invasive Australian range using whole-genome pool-seq data from regionally distributed populations. Our reference genome comprises 132 autosomal contigs, with a total length of 48.90 Mb. We observed a large complex of ace genes, which has presumably evolved from a long history of organophosphate selection in H. destructor and may contribute towards organophosphate resistance through copy number variation, target-site mutations and structural variants. In the putative ancestral H. destructor ace gene, we identified three target-site mutations (G119S, A201S and F331Y) segregating in organophosphate-resistant populations. Additionally, we identified two new para sodium channel gene mutations (L925I and F1020Y) that may contribute to pyrethroid resistance. Regional structuring observed in population genomic analyses indicates that gene flow in H. destructor does not homogenize populations across large geographic distances. However, our demographic analyses were equivocal on the magnitude of gene flow; the short invasion history of H. destructor makes it difficult to distinguish scenarios of complete isolation vs. ongoing migration. Nonetheless, we identified clear signatures of reduced genetic diversity and smaller inferred effective population sizes in eastern vs. western populations, which is consistent with the stepping-stone invasion pathway of this pest in Australia. These new insights will inform development of diagnostic genetic markers of resistance, further investigation into the multifaceted organophosphate resistance mechanism and predictive modelling of resistance evolution and spread.

The Metagenomic Analysis of Viral Diversity in Colorado Potato Beetle Public NGS Data

The Colorado potato beetle (CPB) is one of the most serious insect pests due to its high ecological plasticity and ability to rapidly develop resistance to insecticides. The use of biological insecticides based on viruses is a promising approach to control insect pests, but the information on viruses which infect leaf feeding beetles is scarce. We performed a metagenomic analysis of 297 CPB genomic and transcriptomic samples from the public National Center for Biotechnology Information Sequence Read Archive (NCBI SRA) database. The reads that were not aligned to the reference genome were assembled with metaSPAdes, and 13314 selected contigs were analyzed with BLAST tools. The contigs and non-aligned reads were also analyzed with Kraken2 software. A total of 3137 virus-positive contigs were attributed to different viruses belonging to 6 types, 17 orders, and 32 families, matching over 97 viral species. The annotated sequences can be divided into several groups: those that are homologous to genetic sequences of insect viruses (Adintoviridae, Ascoviridae, Baculoviridae, Dicistroviridae, Chuviridae, Hytrosaviridae, Iflaviridae, Iridoviridae, Nimaviridae, Nudiviridae, Phasmaviridae, Picornaviridae, Polydnaviriformidae, Xinmoviridae etc.), plant viruses (Betaflexiviridae, Bromoviridae, Kitaviridae, Potyviridae), and endogenous retroviral elements (Retroviridae, Metaviridae). Additionally, the full-length genomes and near-full length genome sequences of several viruses were assembled. We also found sequences belonging to Bracoviriform viruses and, for the first time, experimentally validated the presence of bracoviral genetic fragments in the CPB genome. Our work represents the first attempt to discover the viral genetic material in CPB samples, and we hope that further studies will help to identify new viruses to extend the arsenal of biopesticides against CPB.

Chromosome-level genome assembly of the Colorado potato beetle, Leptinotarsa decemlineata

The Colorado potato beetle (Leptinotarsa decemlineata) is one of the most notorious insect pests of potatoes globally. Here, we generated a high-quality chromosome-level genome assembly of L. decemlineata using a combination of the PacBio HiFi sequencing and Hi-C scaffolding technologies. The genome assembly (-1,008 Mb) is anchored to 18 chromosomes (17 + XO), with a scaffold N50 of 58.32 Mb. It contains 676 Mb repeat sequences and 29,606 protein-coding genes. The chromosome-level genome assembly of L. decemlineata provides in-depth knowledge and will be a helpful resource for the beetle and invasive biology research communities.

A draft Diabrotica virgifera virgifera genome: insights into control and host plant adaption by a major maize pest insect

BACKGROUND

Adaptations by arthropod pests to host plant defenses of crops determine their impacts on agricultural production. The larval host range of western corn rootworm, Diabrotica virgifera virgifera (Coleoptera: Chrysomelidae), is restricted to maize and a few grasses. Resistance of D. v. virgifera to crop rotation practices and multiple insecticides contributes to its status as the most damaging pest of cultivated maize in North America and Europe. The extent to which adaptations by this pest contributes to host plant specialization remains unknown.

RESULTS

A 2.42 Gb draft D. v. virgifera genome, Dvir_v2.0, was assembled from short shotgun reads and scaffolded using long-insert mate-pair, transcriptome and linked read data. K-mer analysis predicted a repeat content of ≥ 61.5%. Ortholog assignments for Dvir_2.0 RefSeq models predict a greater number of species-specific gene duplications, including expansions in ATP binding cassette transporter and chemosensory gene families, than in other Coleoptera. A majority of annotated D. v. virgifera cytochrome P450s belong to CYP4, 6, and 9 clades. A total of 5,404 transcripts were differentially-expressed between D. v. virgifera larvae fed maize roots compared to alternative host (Miscanthus), a marginal host (Panicum virgatum), a poor host (Sorghum bicolor) and starvation treatments; Among differentially-expressed transcripts, 1,908 were shared across treatments and the least number were between Miscanthus compared to maize. Differentially-expressed transcripts were enriched for putative spliceosome, proteosome, and intracellular transport functions. General stress pathway functions were unique and enriched among up-regulated transcripts in marginal host, poor host, and starvation responses compared to responses on primary (maize) and alternate hosts.

CONCLUSIONS

Manual annotation of D. v. virgifera Dvir_2.0 RefSeq models predicted expansion of paralogs with gene families putatively involved in insecticide resistance and chemosensory perception. Our study also suggests that adaptations of D. v. virgifera larvae to feeding on an alternate host plant invoke fewer transcriptional changes compared to marginal or poor hosts. The shared up-regulation of stress response pathways between marginal host and poor host, and starvation treatments may reflect nutrient deprivation. This study provides insight into transcriptomic responses of larval feeding on different host plants and resources for genomic research on this economically significant pest of maize.

RNAi of Mannosidase-Ia in the Colorado potato beetle and changes in the midgut and peritrophic membrane

BACKGROUND

In addition to its role in the digestive system, the peritrophic membrane (PM) provides a physical barrier protecting the intestine from abrasion and against pathogens. Because of its sensitivity to RNA interference (RNAi), the notorious pest insect, the Colorado potato beetle (CPB, Leptinotarsa decemlineata), has become a model insect for functional studies. Previously, RNAi-mediated silencing of Mannosidase-Ia (ManIa), a key enzyme in the transition from high-mannose glycan moieties to paucimannose N-glycans, was shown to disrupt the transition from larva to pupa and the metamorphosis into adult beetles. While these effects at the organismal level were interesting in a pest control context, the effects at the organ or tissue level and also immune effects have not been investigated yet. To fill this knowledge gap, we performed an analysis of the midgut and PM in ManIa-silenced insects.

RESULTS

As marked phenotype, the ManIa^RNAi insects, the PM pore size was found to be decreased when compared to the control GFP^RNAi insects. These smaller pores are related to the observation of thinner microvilli (Mv) on the epithelial cells of the midgut of ManIa^RNAi insects. A midgut and PM proteome study and reverse transcription quantitative polymerase chain reaction (RT-qPCR) analysis with a selection of marker genes was performed to characterize the midgut cells and understand their response to the silencing of ManIa. In agreement with the loss of ManIa activity, an accumulation of high-mannose N-glycans was observed in the ManIa-silenced insects. As a pathogen-associated molecular pattern (PAMP), the presence of these glycan structures could trigger the activation of the immune pathways.

CONCLUSION

The observed decrease in PM pore size could be a response to prevent potential pathogens to access the midgut epithelium. This hypothesis is supported by the strong increase in transcription levels of the anti-fungal peptide drosomycin-like in ManIa^RNAi insects, although further research is required to elucidate this possibility. The potential immune response in the midgut and the smaller pore size in the PM shed a light on the function of the PM as a physical barrier and provide evidence for the relation between the Mv and PM. © 2022 Society of Chemical Industry.

Bee-safe peptidomimetic acaricides achieved by comparative genomics

The devastating Varroa mite (Varroa destructor Anderson and Trueman) is an obligatory ectoparasite of the honey bee, contributing to significant colony losses in North America and throughout the world. The limited number of conventional acaricides to reduce Varroa mites and prevent disease in honey bee colonies is challenged with wide-spread resistance and low target-site selectivity. Here, we propose a biorational approach using comparative genomics for the development of honey bee-safe and selective acaricides targeting the Varroa mite-specific neuropeptidergic system regulated by proctolin, which is lacking in the honey bee. Proctolin is a highly conserved pentapeptide RYLPT (Arg-Tyr-Leu-Pro-Thr) known to act through a G protein-coupled receptor to elicit myotropic activity in arthropod species. A total of 33 different peptidomimetic and peptide variants were tested on the Varroa mite proctolin receptor. Ligand docking model and mutagenesis studies revealed the importance of the core aromatic residue Tyr2 in the proctolin ligand. Peptidomimetics were observed to have significant oral toxicity leading to the paralysis and death of Varroa mites, while there were no negative effects observed for honey bees. We have demonstrated that a taxon-specific physiological target identified by advanced genomics information offers an opportunity to develop Varroa mite-selective acaricides, hence, expedited translational processes.

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.

The N-glycosylation-related genes as potential targets for RNAi-mediated pest control of the Colorado potato beetle (Leptinotarsa decemlineata)

BACKGROUND

N-glycosylation is one of the most common and important post-translational modifications in the eukaryotic cell. The study of protein N-glycosylation in several model insects confirmed the importance of this process in insect development, immunity, survival and fertility. The Colorado potato beetle (Leptinotarsa decemlineata) (CPB) is a common pest of Solanaceae crops. With the infamous title of champion of insecticide resistance, novel pest control strategies for this insect are needed. Luckily this pest insect is reported as very sensitive for the post-genomic technology of RNA interference (RNAi).

RESULTS

In this project, we investigated the importance of N-glycosylation in the survival and development of CPB using RNAi-mediated gene silencing of N-glycosylation-related genes (NGRGs) during the different transition steps from the larva, through the pupa to the adult stage. High mortality was observed in the larval stage with the silencing of early NGRGs, as STT3a, DAD1 and GCS1. With dsRNA against middle NGRGs, abnormal phenotypes at the ecdysis process and adult formation were observed, while the silencing of late NGRGs did not cause mortality.

CONCLUSION

The lethal phenotypes observed on silencing of the genes involved in the early processing steps of the N-glycosylation pathway suggest these genes are good candidates for RNAi-mediated control of CPB. Next to the gene-specific mechanism of RNAi for biosafety and possible implementation in integrated pest management, we believe these early NGRGs provide a possible alternative to the well-known target genes Snf7 and vacuolar ATPases that are now used in the first commercial RNAi-based products and thus they may be useful in the context of proactive resistance management. © 2021 Society of Chemical Industry.

Population structure and pattern of geographic differentiation of Colorado potato beetle, Leptinotarsa decemlineata (Coleoptera: Chrysomelidae) in Turkey

BACKGROUND

The Colorado potato beetle (CPB) is the most harmful pest of potato in potato cultivation regions globally. Although it is an economically important agricultural pest, the population structure and colonization route of this species in Turkey are uncertain. We used microsatellite and mitochondrial DNA (mtDNA) markers to obtain information about the population source, structure and bio-invasion route of CPB populations in Turkey.

RESULTS

The common single mtDNA haplotype in European CPB populations was obtained in all Turkish CPB populations based on mtDNA data analysis. However, microsatellites revealed a low level of genetic variation in CPB populations. The results of microsatellite analysis [factorial correspondence analysis (FCA), Bayesian analysis of genetic population structure (BAPS), unweighted pair group method with arithmetic mean (UPGMA) dendrogram, F-statistics and Nei's distances] indicated three groups for invasive CPB: Thrace-Marmara and Aegean; Black Sea, Central Anatolia and Mediterranean; Northeastern Anatolia. Region-specific alleles have been identified in regions, where commercial potato cultivation and insecticide use are intensive.

CONCLUSION

The detection of a single fixed European haplotype in all Turkish populations has proved that CPB in Turkey originated from Europe as a result of a founder event occurred in European populations. Low genetic variation was due to the short time period since the spread of CPB from America to Europe. The highest number of private alleles were found in the top commercial potato cultivation region-Central Anatolia from where the CPB populations spread to other parts of Turkey. © 2021 Society of Chemical Industry.

Unveiling gene expression regulation of the Bacillus thuringiensis Cry3Aa toxin receptor ADAM10 by the potato dietary miR171c in Colorado potato beetle

BACKGROUND

The Colorado potato beetle (CPB) is a worldwide devastating pest of potato plants and other Solanaceae characterized by its remarkable ability to evolve resistance to insecticides. Bacillus thuringiensis (Bt) Cry3Aa toxin represents an environmentally safe alternative for CPB control but larvae susceptibility to this toxin has been reported to vary depending on the host plant on which larvae feed. To gain more insight into how nutrition mediates Bt tolerance through effects on gene expression, here we explored the post-transcriptional regulation by microRNAs (miRNAs) of the CPB-ADAM10 gene encoding the Cry3Aa toxin functional receptor ADAM10.

RESULTS

The lower CPB-ADAM10 gene expression in CPB larvae fed on potato plants cv. Vivaldi than those fed on potato cv. Monalisa or tomato plants was inversely related to Cry3Aa toxicity. By high-throughput sequencing we identified seven CPB miRNAs and one potato miRNA predicted to base pair with the CPB-ADAM10 messenger RNA. No differential expression of the endogenous lde-miR1175-5p was found in larvae feeding on any of the two potato plant varieties. However, statistically significant increased amounts of potato stu-miR171c-5p were detected in CPB larvae fed on potato cv. Vivaldi compared to larvae fed on potato cv. Monalisa.

CONCLUSION

Our results support a role for dietary miRNAs in Bt toxicity by regulating the CPB-ADAM10 gene encoding the Cry3Aa toxin receptor ADAM10 in CPB larvae and opening up the possibility of exploiting plant natural variation in miRNAs to provide more sustainable potato crop protection against CPB. © 2021 The Authors. Pest Management Science published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

RNA interference targeting Ras GTPase gene Ran causes larval and adult lethality in Leptinotarsa decemlineata

BACKGROUND

RNA interference (RNAi) is a breakthrough technology in pest control. It is highly efficient to Coleopteran pests such as the Colorado potato beetle Leptinotarsa decemlineata, a serious pest defoliator mainly attacking potatoes worldwide. The first step for effective pest control by RNAi is the development of effective and reliable target genes.

RESULTS

Our results revealed that continuous ingestion of dsLdRan for 3 days successfully silenced the target gene, inhibited larval growth and killed 100% L. decemlineata larvae. When the bioassay began at the second-, third/fourth-instar larval stages, the larval lethality mainly occurred at the fourth larval instar and prepupal stages, respectively. Importantly, consumption of dsLdRan for 3 days by the newly-emerged males and females effectively knocked down the target transcript, reduced fresh weights and caused 100% of lethality within a week. The LdRan females possessed underdeveloped ovaries.

CONCLUSION

Considering that the larvae, adults and eggs are simultaneously sited on the potato plants, bacterially-expressed dsLdRan is a potential RNAi-based strategy for managing L. decemlineata in the potato field. © 2022 Society of Chemical Industry.

Chromosome-level genome assembly and population genomic analyses provide insights into adaptive evolution of the red turpentine beetle, Dendroctonus valens

BACKGROUND

Biological invasions are responsible for substantial environmental and economic losses. The red turpentine beetle (RTB), Dendroctonus valens LeConte, is an important invasive bark beetle from North America that has caused substantial tree mortality in China. The lack of a high-quality reference genome seriously limits deciphering the extent to which genetic adaptions resulted in a secondary pest becoming so destructive in its invaded area.

RESULTS

Here, we present a 322.41 Mb chromosome-scale reference genome of RTB, of which 98% of assembled sequences are anchored onto fourteen linkage groups including the X chromosome with a N50 size of 24.36 Mb, which is significantly greater than other Coleoptera species. Repetitive sequences make up 45.22% of the genome, which is higher than four other Coleoptera species, i.e., Mountain pine beetle Dendroctonus ponderosae, red flour beetle Tribolium castaneum, blister beetle Hycleus cichorii, and Colorado potato beetle Leptinotarsa decemlineata. We identify rapidly expanded gene families and positively selected genes in RTB, which may be responsible for its rapid environmental adaptation. Population genetic structure of RTB was revealed by genome resequencing of geographic populations in native and invaded regions, suggesting substantial divergence of the North American population and illustrates the possible invasion and spread route in China. Selective sweep analysis highlighted the enhanced ability of Chinese populations in environmental adaptation.

CONCLUSIONS

Overall, our high-quality reference genome represents an important resource for genomics study of invasive bark beetles, which will facilitate the functional study and decipher mechanism underlying invasion success of RTB by integrating the Pinus tabuliformis genome.

Insights Into Chemosensory Proteins From Non-Model Insects: Advances and Perspectives in the Context of Pest Management

Nowadays, insect chemosensation represents a key aspect of integrated pest management in the Anthropocene epoch. Olfaction-related proteins have been the focus of studies due to their function in vital processes, such ashost finding and reproduction behavior. Hence, most research has been based on the study of model insects, namely Drosophila melanogaster, Bombyx mori or Tribolium castaneum. Over the passage of time and the advance of new molecular techniques, insects considered non-models have been studied, contributing greatly to the knowledge of insect olfactory systems and enhanced pest control methods. In this review, a reference point for non-model insects is proposed and the concept of model and non-model insects is discussed. Likewise, it summarizes and discusses the progress and contribution in the olfaction field of both model and non-model insects considered pests in agriculture.

Recombinant Cathepsin L of Tribolium castaneum and Its Potential in the Hydrolysis of Immunogenic Gliadin Peptides

Wheat gliadins contain a large amount of glutamine- and proline-rich peptides which are not hydrolyzed by human digestive peptidases and can cause autoimmune celiac disease and other forms of gluten intolerance in predisposed people. Peptidases that efficiently cleave such immunogenic peptides can be used in enzyme therapy. The stored product insect pest Tribolium castaneum efficiently hydrolyzes gliadins. The main digestive peptidase of T. castaneum is cathepsin L, which is from the papain C1 family with post-glutamine cleavage activity. We describe the isolation and characterization of T. castaneum recombinant procathepsin L (rpTcCathL1, NP_001164001), which was expressed in Pichia pastoris cells. The activation of the proenzyme was conducted by autocatalytic processing. The effects of pH and proenzyme concentration in the reaction mixture on the processing were studied. The mature enzyme retained high activity in the pH range from 5.0 to 9.0 and displayed high pH-stability from 4.0 to 8.0 at 20 °C. The enzyme was characterized according to electrophoretic mobility under native conditions, activity and stability at various pH values, a sensitivity to various inhibitors, and substrate specificity, and its hydrolytic effect on 8-, 10-, 26-, and 33-mer immunogenic gliadins peptides was demonstrated. Our results show that rTcCathL1 is an effective peptidase that can be used to develop a drug for the enzyme therapy of various types of gluten intolerance.

Genome size evolution in the beetle genus Diabrotica

Diabrocite corn rootworms are one of the most economically significant pests of maize in the United States and Europe and an emerging model for insect-plant interactions. Genome sizes of several species in the genus Diabrotica were estimated using flow cytometry along with that of Acalymma vittatum as an outgroup. Genome sizes ranged between 1.56 and 1.64 gigabase pairs and between 2.26 and 2.59 Gb, respectively, for the Diabrotica subgroups fucata and virgifera; the Acalymma vittatum genome size was around 1.65 Gb. This result indicated that a substantial increase in genome size occurred in the ancestor of the virgifera group. Further analysis of the fucata group and the virgifera group genome sequencing reads indicated that the genome size difference between the Diabrotica subgroups could be attributed to a higher content of transposable elements, mostly miniature inverted-transposable elements and gypsy-like long terminal repeat retroelements.

Recruitment of lysosomal cathepsins B, L and D as digestive enzymes in Coleoptera

The recruitment of the lysosomal cathepsins B (CAB), L (CAL) and D (CAD) as luminal digestive enzymes was investigated in 3 species of beetles. Gene expression was determined by RNA-seq in different regions of the midgut and in the carcasses from the transcriptomes of Dermestes maculatus, Tenebrio molitor and Zabrotes subfasciatus. These data together with phylogenetic analyses, allowed us to identify the sequences of the gene coding for digestive and lysosomal CABs, CADs and CALs in T. molitor and Z. subfasciatus and observe the absence of digestive cathepsins in D. maculatus. Comparisons of structures based on the overall similarity of modelled structures were performed and subsite residues in the lysosomal and digestive CALs were identified by molecular docking. The data showed that S2 subsites are very variable, probably as an adaption to a luminal digestive role. The survey of sequences of the gene coding for cathepsins in the genomes of 13 beetle species from different phylogenetic groups showed that expansion of CAL and CAB genes occurred only in the Cucujiformia clade. Several digestive CABs have a reduced occluding loop, probably to act as digestive enzymes. Pollen-feeding was proposed to be the selective pressure to recruit cathepsins as digestive enzymes in Cucujiformia beetles.

Persistent Parental RNAi in the Beetle Tribolium castaneum Involves Maternal Transmission of Long Double-Stranded RNA

Parental RNA interference (pRNAi) is a powerful and widely used method for gene-specific knockdown. Yet in insects its efficacy varies between species, and how the systemic response is transmitted from mother to offspring remains elusive. Using the beetle Tribolium castaneum, an RT-qPCR strategy to distinguish the presence of double-stranded RNA (dsRNA) from endogenous mRNA is reported. It is found that injected dsRNA is directly transmitted into the egg and persists throughout embryogenesis. Despite this depletion of dsRNA from the mother, it is shown that strong pRNAi can persist for months before waning at strain-specific rates. In seeking the receptor proteins for cellular uptake of long dsRNA into the egg, a phylogenomics profiling approach of candidate proteins is also presented. A visualization strategy based on taxonomically hierarchical assessment of orthology clustering data to rapidly assess gene age and copy number changes, refined by sequence-based evidence, is demonstrated. Repeated losses of SID-1-like channel proteins in the arthropods, including wholesale loss in the Heteroptera (true bugs), which are nonetheless highly sensitive to pRNAi, are thereby documented. Overall, practical considerations for insect pRNAi against a backdrop of outstanding questions on the molecular mechanism of dsRNA transmission for long-term, systemic knockdown are elucidated.

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.

Diversity and Molecular Evolution of Odorant Receptor in Hemipteran Insects

Simple Summary

Insects’ behavior and ecology are closely related to their chemosensory systems, during which odorant receptors (ORs) play an essential role in host recognition. Although OR gene evolution has been studied in many insect orders, a comprehensive evolutionary analysis and expression of OR gene gain and loss events among diverse hemipteran species are still needed. In this study, we identified and analyzed the OR genes from hemipteran species systematically. The number of OR genes discovered in each species ranged from less than ten to hundreds. Gene gain and loss events of OR have occurred in several species in the seven major clades classified through phylogenetic analysis. Then, we discovered the amino acid differences between species to understand the molecular evolution of OR in the order Hemiptera through positive selection. This study lays a foundation for subsequent investigations into the molecular mechanisms of Hemiptera olfactory receptors involved in host recognition.

Abstract

Olfaction is a critical physiologic process for insects to interact with the environment, especially plant-emitted volatiles, during which odorant receptors (ORs) play an essential role in host recognition. Although OR gene evolution has been studied in many insect orders, a comprehensive evolutionary analysis and expression of OR gene gain and loss events among diverse hemipteran species are still required. In this study, we identified and analyzed 887 OR genes from 11 hemipteran species. The number of OR genes discovered in each species ranged from less than ten to hundreds. Phylogenetic analysis revealed that all identified Hemiptera OR genes were classified into seven major clades. Gene gain and loss events of OR have occurred in several species. Then, by positive selection, we discovered the amino acid differences between species to understand the molecular evolution of OR in the order Hemiptera. Additionally, we discussed how evolutionary analysis can aid the study of insect–plant communication. This study lays a foundation for subsequent investigations into the molecular mechanisms of Hemiptera olfactory receptors involved in host recognition.

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.

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.

A scaffold-level genome assembly of a minute pirate bug, Orius laevigatus (Hemiptera: Anthocoridae), and a comparative analysis of insecticide resistance-related gene families with hemipteran crop pests

BACKGROUND

Orius laevigatus, a minute pirate bug, is a highly effective beneficial predator of crop pests including aphids, spider mites and thrips in integrated pest management (IPM) programmes. No genomic information is currently available for O. laevigatus, as is the case for the majority of beneficial predators which feed on crop pests. In contrast, genomic information for crop pests is far more readily available. The lack of publicly available genomes for beneficial predators to date has limited our ability to perform comparative analyses of genes encoding potential insecticide resistance mechanisms between crop pests and their predators. These mechanisms include several gene/protein families including cytochrome P450s (P450s), ATP binding cassette transporters (ABCs), glutathione S-transferases (GSTs), UDP-glucosyltransferases (UGTs) and carboxyl/cholinesterases (CCEs).

METHODS AND FINDINGS

In this study, a high-quality scaffold level de novo genome assembly for O. laevigatus has been generated using a hybrid approach with PacBio long-read and Illumina short-read data. The final assembly achieved a scaffold N50 of 125,649 bp and a total genome size of 150.98 Mb. The genome assembly achieved a level of completeness of 93.6% using a set of 1658 core insect genes present as full-length genes. Genome annotation identified 15,102 protein-coding genes - 87% of which were assigned a putative function. Comparative analyses revealed gene expansions of sigma class GSTs and CYP3 P450s. Conversely the UGT gene family showed limited expansion. Differences were seen in the distributions of resistance-associated gene families at the subfamily level between O. laevigatus and some of its targeted crop pests. A target site mutation in ryanodine receptors (I4790M, PxRyR) which has strong links to diamide resistance in crop pests and had previously only been identified in lepidopteran species was found to also be present in hemipteran species, including O. laevigatus.

CONCLUSION AND SIGNIFICANCE

This assembly is the first published genome for the Anthocoridae family and will serve as a useful resource for further research into target-site selectivity issues and potential resistance mechanisms in beneficial predators. Furthermore, the expansion of gene families often linked to insecticide resistance may be an indicator of the capacity of this predator to detoxify selective insecticides. These findings could be exploited by targeted pesticide screens and functional studies to increase effectiveness of IPM strategies, which aim to increase crop yields by sustainably, environmentally-friendly and effectively control pests without impacting beneficial predator populations.

Primary digestive cathepsins L of Tribolium castaneum larvae: Proteomic identification, properties, comparison with human lysosomal cathepsin L

We previously described the most highly expressed enzymes from the gut of the red flour beetle, Tribolium castaneum, as cathepsins L. In the present study, two C1 family-specific cysteine cathepsin L enzymes from the larval midgut were isolated and identified using MALDI-TOF MS analysis. The isolated T. castaneum cathepsins were characterized according to their specificity against chromogenic and fluorogenic peptide substrates, and the most efficiently hydrolyzed substrate was Z-FR-pNA with Arg in the P1 subsite. The specificity of insect digestive cathepsins was compared with human lysosomal cathepsin L, the well-studied peptidase of the C1 family cathepsins. T. castaneum digestive cathepsins efficiently hydrolyzed substrates with small and uncharged amino acid residues at P1 (Ala, Gln) more than human cathepsin L. In particular, these insect digestive cathepsins cleaved with higher efficiency the analogs of immunogenic peptides of gliadins, which contribute to autoimmune celiac disease in susceptible people, and thus insect enzymes may be useful in enzymatic treatments for this disease. A bioinformatic study supported by the proteomic analysis of the primary structures of the isolated cathepsins was used to compare tertiary models. The phylogenetic analysis of coleopteran and human cathepsins from the L subfamily indicated that insect digestive cathepsins grouped separately from lysosomal cathepsins.

Gut transcriptome of two bark beetle species stimulated with the same kairomones reveals molecular differences in detoxification pathways

Dendroctonus bark beetles are the most destructive agents in coniferous forests. These beetles come into contact with the toxic compounds of their host's chemical defenses throughout their life cycle, some of which are also used by the insects as kairomones to select their host trees during the colonization process. However, little is known about the molecular mechanisms by which the insects counteract the toxicity of these compounds. Here, two sibling species of bark beetles, D. valens and D. rhizophagus, were stimulated with vapors of a blend of their main kairomones (α-pinene, β-pinene and 3-carene), in order to compare the transcriptional response of their gut. A total of 48 180 unigenes were identified in D. valens and 43 704 in D. rhizophagus, in response to kairomones blend. The analysis of differential gene expression showed a transcriptional response in D. valens (739 unigenes, 0.58–10.36 Log2FC) related to digestive process and in D. rhizophagus (322 unigenes 0.87–13.08 Log2FC) related to xenobiotics metabolism. The expression profiles of detoxification genes mainly evidenced the up-regulation of COEs and GSTs in D. valens, and the up-regulation of P450s in D. rhizophagus. Results suggest that terpenes metabolism comes accompanied by an integral hormetic response, result of compensatory mechanisms, including the activation of other metabolic pathways, to ensure the supply of energy and the survival of organisms which is specific for each species, according to its life history and ecological strategy.

The transposable element-rich genome of the cereal pest Sitophilus oryzae

BACKGROUND

The rice weevil Sitophilus oryzae is one of the most important agricultural pests, causing extensive damage to cereal in fields and to stored grains. S. oryzae has an intracellular symbiotic relationship (endosymbiosis) with the Gram-negative bacterium Sodalis pierantonius and is a valuable model to decipher host-symbiont molecular interactions.

RESULTS

We sequenced the Sitophilus oryzae genome using a combination of short and long reads to produce the best assembly for a Curculionidae species to date. We show that S. oryzae has undergone successive bursts of transposable element (TE) amplification, representing 72% of the genome. In addition, we show that many TE families are transcriptionally active, and changes in their expression are associated with insect endosymbiotic state. S. oryzae has undergone a high gene expansion rate, when compared to other beetles. Reconstruction of host-symbiont metabolic networks revealed that, despite its recent association with cereal weevils (30 kyear), S. pierantonius relies on the host for several amino acids and nucleotides to survive and to produce vitamins and essential amino acids required for insect development and cuticle biosynthesis.

CONCLUSIONS

Here we present the genome of an agricultural pest beetle, which may act as a foundation for pest control. In addition, S. oryzae may be a useful model for endosymbiosis, and studying TE evolution and regulation, along with the impact of TEs on eukaryotic genomes.

Structural and Functional Characterization of One Unclassified Glutathione S-Transferase in Xenobiotic Adaptation of Leptinotarsa decemlineata

Arthropod Glutathione S-transferases (GSTs) constitute a large family of multifunctional enzymes that are mainly associated with xenobiotic or stress adaptation. GST-mediated xenobiotic adaptation takes place through direct metabolism or sequestration of xenobiotics, and/or indirectly by providing protection against oxidative stress induced by xenobiotic exposure. To date, the roles of GSTs in xenobiotic adaptation in the Colorado potato beetle (CPB), a notorious agricultural pest of plants within Solanaceae, have not been well studied. Here, we functionally expressed and characterized an unclassified-class GST, LdGSTu1. The three-dimensional structure of the LdGSTu1 was solved with a resolution up to 1.8 Å by X-ray crystallography. The signature motif VSDGPPSL was identified in the “G-site”, and it contains the catalytically active residue Ser14. Recombinant LdGSTu1 was used to determine enzyme activity and kinetic parameters using 1-chloro-2, 4-dinitrobenzene (CDNB), GSH, p-nitrophenyl acetate (PNA) as substrates. The enzyme kinetic parameters and enzyme-substrate interaction studies demonstrated that LdGSTu1 could catalyze the conjugation of GSH to both CDNB and PNA, with a higher turnover number for CDNB than PNA. The LdGSTu1 enzyme inhibition assays demonstrated that the enzymatic conjugation of GSH to CDNB was inhibited by multiple pesticides, suggesting a potential function of LdGSTu1 in xenobiotic adaptation.

Sequence variation and regulatory variation in acetylcholinesterase genes contribute to insecticide resistance in different populations of Leptinotarsa decemlineata

Although insect herbivores are known to evolve resistance to insecticides through multiple genetic mechanisms, resistance in individual species has been assumed to follow the same mechanism. While both mutations in the target site insensitivity and increased amplification are known to contribute to insecticide resistance, little is known about the degree to which geographic populations of the same species differ at the target site in a response to insecticides. We tested structural (e.g., mutation profiles) and regulatory (e.g., the gene expression of Ldace1 and Ldace2, AChE activity) differences between two populations (Vermont, USA and Belchow, Poland) of the Colorado potato beetle, Leptinotarsa decemlineata in their resistance to two commonly used groups of insecticides, organophosphates, and carbamates. We established that Vermont beetles were more resistant to azinphos-methyl and carbaryl insecticides than Belchow beetles, despite a similar frequency of resistance-associated alleles (i.e., S291G) in the Ldace2 gene. However, the Vermont population had two additional amino acid replacements (G192S and F402Y) in the Ldace1 gene, which were absent in the Belchow population. Moreover, the Vermont population showed higher expression of Ldace1 and was less sensitive to AChE inhibition by azinphos-methyl oxon than the Belchow population. Therefore, the two populations have evolved different genetic mechanisms to adapt to organophosphate and carbamate insecticides.