Cereal weevils' antimicrobial peptides: at the crosstalk between development, endosymbiosis and immune response

Interactions between animals and microbes are ubiquitous in nature and strongly impact animal physiology. These interactions are shaped by the host immune system, which responds to infections and contributes to tailor the associations with beneficial microorganisms. In many insects, beneficial symbiotic associations not only include gut commensals, but also intracellular bacteria, or endosymbionts. Endosymbionts are housed within specialized host cells, the bacteriocytes, and are transmitted vertically across host generations. Host-endosymbiont co-evolution shapes the endosymbiont genome and host immune system, which not only fights against microbial intruders, but also ensures the preservation of endosymbionts and the control of their load and location. The cereal weevil Sitophilus spp. is a remarkable model in which to study the evolutionary adaptation of the immune system to endosymbiosis owing to its binary association with a unique, relatively recently acquired nutritional endosymbiont, Sodalis pierantonius. This Gram-negative bacterium has not experienced the genome size shrinkage observed in long-term endosymbioses and has retained immunogenicity. We focus here on the sixteen antimicrobial peptides (AMPs) identified in the Sitophilus oryzae genome and their expression patterns in different tissues, along host development or upon immune challenges, to address their potential functions in the defensive response and endosymbiosis homeostasis along the insect life cycle. This article is part of the theme issue 'Sculpting the microbiome: how host factors determine and respond to microbial colonization'.

Genetic legacies of mega-landslides: Cycles of isolation and contact across flank collapses in an oceanic island

Catastrophic flank collapses are recognized as important drivers of insular biodiversity dynamics, through the disruption of species ranges and subsequent allopatric divergence. However, little empirical data supports this conjecture, with their evolutionary consequences remaining poorly understood. Using genome-wide data within a population genomics and phylogenomics framework, we evaluate how mega-landslides have impacted evolutionary and demographic history within a species complex of weevils (Curculionidae) within the Canary Island of Tenerife. We reveal a complex genomic landscape, within which individuals of single ancestry were sampled in areas characterized by long-term geological stability, relative to the timing of flank collapses. In contrast, individuals of admixed ancestry were almost exclusively sampled within the boundaries of flank collapses. Estimated divergence times among ancestral populations aligned with the timings of mega-landslide events. Our results provide first evidence for a cyclical dynamic of range fragmentation and secondary contact across flank collapse landscapes, with support for a model where this dynamic is mediated by Quaternary climate oscillations. The context within which we reveal climate and topography to interact cyclically through time to shape the geographic structure of genetic variation, together with related recent work, highlights the importance of topoclimatic phenomena as an agent of diversification within insular invertebrates.

Comparative transcriptomic analysis of chemoreceptors in two sympatric scarab beetles, Hylamorpha elegans and Brachysternus prasinus

Chemoreception through odorant receptors (ORs), ionotropic receptors (IRs) and gustatory receptors (GRs) represents the functions of key proteins in the chemical ecology of insects. Recent studies have identified chemoreceptors in coleopterans, facilitating the evolutionary analysis of not only ORs but also IRs and GRs. Thus, Cerambycidae, Tenebrionidae and Curculionidae have received increased attention. However, knowledge of the chemoreceptors from Scarabaeidae is still limited, particularly for those that are sympatric. Considering the roles of chemoreceptors, this analysis could shed light on evolutionary processes in the context of sympatry. Therefore, the aim of this study was to identify and compare the repertoires of ORs, GRs and IRs between two sympatric scarab beetles, Hylamorpha elegans and Brachysternus prasinus. Here, construction of the antennal transcriptomes of both scarab beetle species and analyses of their phylogeny, molecular evolution and relative expression were performed. Thus, 119 new candidate chemoreceptors were identified for the first time, including 17 transcripts for B. prasinus (1 GR, 3 IRs and 13 ORs) and 102 for H. elegans (22 GRs, 14 IRs and 66 ORs). Orthologs between the two scarab beetle species were found, revealing specific expansions as well as absence in some clades. Purifying selection appears to have occurred on H. elegans and B. prasinus ORs. Further efforts will be focused on target identification to characterize kairomone and/or pheromone receptors.

Exploring the functional profiles of odorant binding proteins crucial for sensing key odorants in the new leaves of coconut palms in Rhynchophorus ferrugineus

The red palm weevil (RPW), Rhynchophorus ferrugineus (Curculionidae: Coleoptera) is a highly destructive global pest of coconut trees, with a preference for laying its eggs on new leaves. Females can identify where to lay eggs by using their sense of smell to detect specific odorants found in new leaves. In this study, we focused on the two odorants commonly found in new leaves by GC-MS: trans, trans-2,4-nonadienal and trans-2-nonenal. Our behavioral assays demonstrated a significant attraction of females to both of these odorants, with their electrophysiological responses being dose-dependent. Furthermore, we examined the expression patterns induced by these odorants in eleven RferOBP genes. Among them, RferOBP3 and RferOBP1768 exhibited the most significant and simultaneous upregulation. To further understand the role of these two genes, we conducted experiments with females injected with OBP-dsRNA. This resulted in a significant decrease in the expression of RferOBP3 and RferOBP1768, as well as impaired the perception of the two odorants. A fluorescence competitive binding assay also showed that both RferOBPs strongly bound to the odorants. Additionally, sequence analysis revealed that these two RferOBPs belong to the Minus-C family and possess four conserved cysteines. Molecular docking simulations showed strong interactions between these two RferOBPs and the odorant molecules. Overall, our findings highlight the crucial role of RferOBP3 and RferOBP1768 in the olfactory perception of the key odorants in coconut palm new leaves. This knowledge significantly improves our understanding of how RPW females locate sites for oviposition and lays the foundation for future research on the development of environmentally friendly pest attractants.

Comparative transcriptomics reveals the conservation and divergence of reproductive genes across three sympatric Tomicus bark beetles

Three tree-killing bark beetles belonging to the genus Tomicus, Tomicus yunnanensis, Tomicus brevipilosus and Tomicus minor (Coleoptera; Curculionidae, Scolytinae), are serious wood-borers with larvae feeding on the phloem tissues of Pinus yunnanensis. The three Tomicus beetles, in some cases, coexist in a same habitat, providing a best system for exploring the conservation and divergence of reproductive genes. Here, we applied comparative transcriptomics and molecular biology approaches to characterize reproductive-related genes in three sympatric Tomicus species. Illumina sequencing of female and male reproductive systems and residual bodies generated a large number of clean reads, representing 185,920,232 sequences in T. yunnanensis, 169,153,404 in T. brevipilosus and 178,493,176 in T. minor that were assembled into 32,802, 56,912 and 33,670 unigenes, respectively. The majority of the genes had detectable expression in reproductive tissues (FPKM >1), particularly those genes in T. brevipilosus accounting for 76.61 % of the total genes. From the transcriptomes, totally 838 genes encoding 463 detoxification enzymes, 339 chemosensory membrane proteins and 36 ionotropic glutamate receptors (iGluRs) were identified, including 622 reproductive tissue-expressed genes. Of these, members of carboxylesterases (COEs), ionotropic receptors (IRs), sensory neuron membrane proteins (SNMPs) and iGluRs were highly conserved in gene numbers and sequence identities across three Tomicus species. Further, expression profiling analyses revealed a number of genes expressed in reproductive tissues and the diverse expression characteristics in these beetles. The results provide evidence for the conservation and differences of reproductive genes among three sympatric closely related beetles, helping understand their different reproductive strategies and the maximization of the reproductive success.

The genome of the invasive and broadly polyphagous Diaprepes root weevil, Diaprepes abbreviatus (Coleoptera), reveals an arsenal of putative polysaccharide-degrading enzymes

The Diaprepes root weevil (DRW), Diaprepes abbreviatus, is a broadly polyphagous invasive pest of agriculture in the southern United States and the Caribbean. Its genome was sequenced, assembled, and annotated to study genomic correlates of specialized plant-feeding and invasiveness and to facilitate the development of new methods for DRW control. The 1.69 Gb D. abbreviatus genome assembly was distributed across 653 contigs, with an N50 of 7.8 Mb and the largest contig of 62 Mb. Most of the genome was comprised of repetitive sequences, with 66.17% in transposable elements, 5.75% in macrosatellites, and 2.06% in microsatellites. Most expected orthologous genes were present and fully assembled, with 99.5% of BUSCO genes present and 1.5% duplicated. One hundred and nine contigs (27.19 Mb) were identified as putative fragments of the X and Y sex chromosomes, and homology assessment with other beetle X chromosomes indicated a possible sex chromosome turnover event. Genome annotation identified 18,412 genes, including 43 putative horizontally transferred (HT) loci. Notably, 258 genes were identified from gene families known to encode plant cell wall degrading enzymes and invertases, including carbohydrate esterases, polysaccharide lyases, and glycoside hydrolases (GH). GH genes were unusually numerous, with 239 putative genes representing 19 GH families. Interestingly, several other beetle species with large numbers of GH genes are (like D. abbreviatus) successful invasive pests of agriculture or forestry.

Near-chromosomal-level genome of the red palm weevil (Rhynchophorus ferrugineus), a potential resource for genome-based pest control

The red palm weevil (RPW) is a highly destructive pest that mainly affects palms, particularly date palms (Phoenix dactylifera), in the Arabian Gulf region. In this study, we present a near-chromosomal-level genome assembly of the RPW using a combination of PacBio HiFi and Dovetail Omini-C reads. The final genome assembly is around 779 Mb in size, with an N50 of ~43 Mb, consistent with our previous flow cytometry estimates. The completeness of the genome was confirmed through BUSCO analysis, which indicates the presence of 99.5% of BUSCO single copy orthologous genes. The genome annotation identified a total of 29,666 protein-coding, 1,091 tRNA and 543 rRNA genes. Overall, the proposed genome assembly is significantly superior to existing assemblies in terms of contiguity, integrity, and genome completeness.

Bacillus thuringiensis Bt_UNVM-84, a Novel Strain Showing Insecticidal Activity against Anthonomus grandis Boheman (Coleoptera: Curculionidae)

Bacillus thuringiensis is a Gram-positive bacterium known for its insecticidal proteins effective against various insect pests. However, limited strains and proteins target coleopteran pests like Anthonomous grandis Boheman, causing substantial economic losses in the cotton industry. This study focuses on characterizing a Bacillus sp. strain, isolated from Oncativo (Argentina), which exhibits ovoid to amorphous parasporal crystals and was designated Bt_UNVM-84. Its genome encodes genes for the production of two pairs of binary Vpb1/Vpa2 proteins and three Cry-like proteins showing similarity with different Cry8 proteins. Interestingly, this gene content was found to be conserved in a previously characterized Argentine isolate of B. thuringiensis designated INTA Fr7-4. SDS-PAGE analysis revealed a major band of 130 kDa that is proteolytically processed to an approximately 66-kDa protein fragment by trypsin. Bioassays performed with spore-crystal mixtures demonstrated an interesting insecticidal activity against the cotton boll weevil A. grandis neonate larvae, resulting in 91% mortality. Strain Bt_UNVM-84 is, therefore, an interesting candidate for the efficient biological control of this species, causing significant economic losses in the cotton industry in the Americas.

Interaction of insecticidal proteins from Pseudomonas spp. and Bacillus thuringiensis for boll weevil management

Cotton crop yields are largely affected by infestations of Anthonomus grandis, which is its main pest. Although Bacillus thuringiensis (Bt) derived proteins can limit insect pest infestations, the diverse use of control methods becomes a viable alternative in order to prolong the use of technology in the field. One of the alternative methods to Bt technology has been the utilization of certain Pseudomonas species highly efficient in controlling coleopteran insects have been used to produce highly toxic insecticidal proteins. This study aimed to evaluate the toxicity of IPD072Aa and PIP-47Aa proteins, isolated from Pseudomonas spp., in interaction with Cry1Ia10, Cry3Aa, and Cry8B proteins isolated from B. thuringiensis, to control A. grandis in cotton crops. The genes IPD072Aa and PIP-47Aa were synthesized and cloned into a pET-SUMO expression vector. Moreover, Cry1Ia10, Cry3Aa, and Cry8B proteins were obtained by inducing recombinant E. coli clones, which were previously acquired by our research group from the Laboratory of Bacteria Genetics and Applied Biotechnology (LGBBA). These proteins were visualized in SDS-PAGE, quantified, and incorporated into an artificial diet to estimate their lethal concentrations (LC) through individual or combined bioassays. The results of individual toxicity revealed that IPD072Aa, PIP-47Aa, Cry1Ia10, Cry3Aa, and Cry8B were efficient in controlling A. grandis, with the latter being the most toxic. Regarding interaction assays, a high synergistic interaction was observed between Cry1Ia10 and Cry3Aa. All interactions involving Cry3Aa and PIP-47Aa, when combined with other proteins, showed a clear synergistic effect. Our findings highlighted that the tested proteins in combination, for the most part, increase toxicity against A. grandis neonate larvae, suggesting possible constructions for pyramiding cotton plants to the manage and the control boll weevils.

Molecular and Functional Characterization of Pheromone Binding Protein 2 from Cyrtotrachelus buqueti (Coleoptera: Curculionidae)

Pheromone-binding proteins (PBPs) play important roles in binding and transporting sex pheromones. However, the PBP genes identified in coleopteran insects and their information sensing mechanism are largely unknown. Cyrtotrachelus buqueti (Coleoptera: Curculionidae) is a major insect pest of bamboo plantations. In this study, a novel PBP gene, CbuqPBP2, from C. buqueti was functionally characterized. CbuqPBP2 was more abundantly expressed in the antennae of both sexes than other body parts, and its expression level was significantly male-biased. Fluorescence competitive binding assays showed that CbuqPBP2 exhibited the strongest binding affinity to dibutyl phthalate (Ki = 6.32 μM), followed by styrene (Ki = 11.37 μM), among twelve C. buqueti volatiles. CbuqPBP2, on the other hand, showed high binding affinity to linalool (Ki = 10.55), the main volatile of host plant Neosinocalamus affinis. Furthermore, molecular docking also demonstrated the strong binding ability of CbuqPBP2 to dibutyl phthalate, styrene, and linalool, with binding energy values of -5.7, -6.6, and -6.0 kcal/mol, respectively, and hydrophobic interactions were the prevailing forces. The knockdown of CbuqPBP2 expression via RNA interference significantly reduced the electroantennography (EAG) responses of male adults to dibutyl phthalate and styrene. In conclusion, these results will be conducive to understanding the olfactory mechanisms of C. buqueti and promoting the development of novel strategies for controlling this insect pest.

Dryocoetiops krivetsae sp. n. (Coleoptera: Curculionidae: Scolytinae: Dryocoetini), the northernmost species of the genus: conflicts between molecular and morphological data in the tribe Dryocoetini

A new species of bark beetle, Dryocoetiops krivetsae Kerchev, Mandelshtam, Bykov et Ilinsky, sp. n., from the southern part of Primorsky Krai (Russian Far East), is described. This is the northernmost discovery of the Oriental genus Dryocoetiops Schedl, 1957, and the first record of the genus in the Russian fauna. Phylogenetic analysis of nuclear and mitochondrial genes i) confirms an independent lineage of Dryocoetiops krivetsae and ii) indicates numerous conflicts of genera relationships in the tribe Dryocoetini.

Transcriptome profiling and in silico docking analysis of phosphine resistance in rice weevil, Sitophilus oryzae (Coleoptera: Curculionidae)

The rice weevil, Sitophilus oryzae (Linnaeus, Coleoptera: Curculionidae), is a serious cosmopolitan pest that affects grain in storage and has developed high levels of resistance toward phosphine. In this study, RNA-seq data was used to study the phosphine resistance mechanisms in S. oryzae. Resistant and susceptible populations of S. oryzae were identified based on phosphine bioassays conducted in 32 populations collected across Tamil Nadu, India. Differential expression of mitochondrial (COX1, COX2, COX3, ND2, ND3, ATP6, and ATP8) and detoxification genes (Cyps, Gsts, and Cbe) were observed in the resistant and susceptible populations of S. oryzae. The previously characterized phosphine resistant gene, dld (dihydrolipoamide dehydrogenase) linked to the rph2 locus, was found to be up-regulated in resistant S. oryzae population (ISO-TNAU-RT) treated with phosphine. Also, the genes involved in Tricarboxylic acid (TCA) cycle were significantly down-regulated. In addition, a significant up-regulation in the expression of the antioxidant enzymes superoxide dismutase (2.5×) and catalase (2.1×) in ISO-TNAU-RT populations was recorded. Furthermore, a distinct amino acid substitution, Lysine > Glutamic acid (K141E) was identified in resistant phenotypes. In silico docking studies of both resistant and susceptible DLD protein with phosphine molecule revealed that the amino acid residues involved in the interaction were different. This suggested that the amino acid substitution might lead to structural modifications which reduces the affinity of the target (phosphine). This study provides insight on the various genes, pathways, and functional mechanisms having a significant role in phosphine resistance in S. oryzae.

Phylogenomics illuminates the phylogeny of flower weevils (Curculioninae) and reveals ten independent origins of brood-site pollination mutualism in true weevils

Weevils are an unusually species-rich group of phytophagous insects for which there is increasing evidence of frequent involvement in brood-site pollination. This study examines phylogenetic patterns in the emergence of brood-site pollination mutualism among one of the most speciose beetle groups, the flower weevils (subfamily Curculioninae). We analysed a novel phylogenomic dataset consisting of 214 nuclear loci for 202 weevil species, with a sampling that mainly includes flower weevils as well as representatives of all major lineages of true weevils (Curculionidae). Our phylogenomic analyses establish a uniquely comprehensive phylogenetic framework for Curculioninae and provide new insights into the relationships among lineages of true weevils. Based on this phylogeny, statistical reconstruction of ancestral character states revealed at least 10 independent origins of brood-site pollination in higher weevils through transitions from ancestral associations with reproductive structures in the larval stage. Broadly, our results illuminate the unexpected frequency with which true weevils-typically specialized phytophages and hence antagonists of plants-have evolved mutualistic interactions of ecological significance that are key to both weevil and plant evolutionary fitness and thus a component of their deeply intertwined macroevolutionary success.

Museomics unveil systematics, diversity and evolution of Australian cycad-pollinating weevils

Weevils have been shown to play significant roles in the obligate pollination of Australian cycads. In this study, we apply museomics to produce a first molecular phylogeny estimate of the Australian cycad weevils, allowing an assessment of their monophyly, placement and relationships. Divergence dating suggests that the Australian cycad weevils originated from the Late Oligocene to the Middle Miocene and that the main radiation of the cycad-pollinating groups occurred from the Middle to the Late Miocene, which is congruent with the diversification of the Australian cycads, thus refuting any notion of an ancient ciophilous system in Australia. Taxonomic studies reveal the existence of 19 Australian cycad weevil species and that their associations with their hosts are mostly non-species-specific. Co-speciation analysis shows no extensive co-speciation events having occurred in the ciophilous system of Australian cycads. The distribution pattern suggests that geographical factors, rather than diversifying coevolution, constitute the overriding process shaping the Australian cycad weevil diversity. The synchronous radiation of cycads and weevil pollinators is suggested to be a result of the post-Oligocene diversification common in Australian organisms.

Identification of Cry toxin receptor genes homologs in a de novo transcriptome of Premnotrypes vorax (Coleoptera: Curculionidae)

The white potato worm Premnotrypes vorax (Hustache) (Coleoptera: Curculionidae) is one of the most destructive insect pests of potato crops in South America. Like many coleopteran insects, P. vorax shows low susceptibility to Cry insecticidal proteins produced by the bacterium Bacillus thuringiensis (Bt). However, the presence of Cry toxin receptors in the midgut of this this insect has never been studied. The main Cry-binding proteins described in other insect species are cadherin (CAD), aminopeptidase N (APN), alkaline phosphatase (ALP) and ATP-binding cassette (ABC) transporters. In this study, we analyzed and validated a de novo assembled transcriptome of Illumina sequencing data to identify and to characterize homologs of Cry toxin receptors. We identified the protein sequences in P. vorax that show high identity with their orthologous sequences of the Cry toxin binding proteins in other coleopteran larvae such as APN, ALP, CAD and ABC transporter. This study provides preliminary identification of putative receptor genes of Cry proteins that would be useful for future studies involving biocontrol of this important potato crop pest.

Phylogenomics of weevils revisited: data curation and modelling compositional heterogeneity

Weevils represent one of the most prolific radiations of beetles and the most diverse group of herbivores on land. The phylogeny of weevils (Curculionoidea) has received extensive attention, and a largely satisfactory framework for their interfamilial relationships has been established. However, a recent phylogenomic study of Curculionoidea based on anchored hybrid enrichment (AHE) data yielded an abnormal placement for the family Belidae (strongly supported as sister to Nemonychidae + Anthribidae). Here we reanalyse the genome-scale AHE data for Curculionoidea using various models of molecular evolution and data filtering methods to mitigate anticipated systematic errors and reduce compositional heterogeneity. When analysed with the infinite mixture model CAT-GTR or using appropriately filtered datasets, Belidae are always recovered as sister to the clade (Attelabidae, (Caridae, (Brentidae, Curculionidae))), which is congruent with studies based on morphology and other sources of molecular data. Although the relationships of the 'higher Curculionidae' remain challenging to resolve, we provide a consistent and robust backbone phylogeny of weevils. Our extensive analyses emphasize the significance of data curation and modelling across-site compositional heterogeneity in phylogenomic studies.

Transcriptome analysis of Euwallacea interjectus reveals differentially expressed unigenes related to developmental stages and egg laying

Euwallacea interjectus (Curculionidae: Scolytinae) is an ambrosia beetle species in its early stages of research. Therefore, studying the related molecular mechanism associated with the development and egg laid is essential. Transcriptome sequencing was used in this study to compare the gene expression of the beetles at different developmental stages and female adults before and after oviposition. A total of 40,047 annotated unigenes were obtained. There were 4225 differentially expressed unigenes (DEUs) from larva to prepupa stage, 3651 DEUs between prepupa and pupa, 1675 DEUs generated from pupa to adult, and 4762 DEUs between females before and after oviposition. The most significant pathway differences between different development stages and before and after oviposition were selected through functional annotation of DEUs between different stages. Among them, there were many pathways related to protein metabolism including: neuroactive ligand-receptor interaction, endoplasmic reticulum and RNA transport. This study provides valuable information on the molecular regulation mechanism of development and the egg laid of E. interjectus.

Acceptance and utilization efficiency of a purple durum wheat genotype by Sitophilus granarius (L.)

The granary weevil (Sitophilus granarius L.) is a major primary pest of stored cereals throughout the world. Among the major classes of plant secondary metabolites, flavonoids can affect insect feeding behaviour and their growth rate. In this study, the susceptibility of an anthocyanin-rich purple durum wheat genotype (T1303) to the granary weevil was evaluated in comparison with two yellow durum (Ofanto) and bread (Mec) wheat varieties. The feeding response and food utilisation efficiency by adult insects was also investigated by calculating nutritional indices in whole flour disk bioassays. Different levels of susceptibility to granary weevil emerged among genotypes tested. The mean food consumption by an insect, F1 progeny, and female parental offspring calculated for the T1303 genotype were significantly lower than those of yellow kernel wheat varieties. Moreover, T1303 genotype induced deterrence in the adult insects as demonstrated by the positive values of the food deterrence index. Besides, relative grow rate and efficiency conversion of ingested food indices were negative for T1303 and positive for both yellow wheat varieties indicating respectively a decrease and an increase of insect body weight during the bioassays. Finally, a higher mortality rate was recorded for insects fed on T1303 flour disks compared to disks obtained from yellow wheat varieties. These results provide evidence for the antifeedant and toxic effects of anthocyanins present in the T1303 pericarp against the granary weevil. Overall, this study contributes new insights into the mechanisms of host acceptance and food utilization by S. granarius and would be useful to identify antifeedant flavonoids as well as to develop varietal resistance-based strategies against this pest.

A new species of Bracon (Braconidae: Braconinae) from central Mexico, probable parasitoid of a weevil that feeds on roots of Argemone ochroleuca Sweet (Papaveraceae)

A new species of the braconine genus Bracon (subgenus Bracon), B. hidalguensis sp. nov., is described from the locality of Tasquillo in the state of Hidalgo, central Mexico. The new species was reared from roots of Argemone ochroleuca Sweet (Papaveraceae), where specimens of the weevil species Conotrachelus leucophaeus (Champion) (Curculionidae) were also obtained and thus probably it represents its host. The new Bracon species was characterised molecularly with DNA barcoding (COI) and a fragment of the variable D23 region of the nuclear ribosomal 28S gene.

An accurate, efficient, and economical identification technology for black twig borer based on species-specific cytochrome C oxidase subunit I PCR assay

Xylosandrus compactus (Eichhoff) (Coleoptera: Curculionidae, Scolytinae) is a worldwide invasive species that causes huge economic loss and environmental damage in many countries. Traditional morphological characteristics make it hard to identify scolytines due to their tiny size. Besides, the intercepted insect samples are incomplete, and the limitation of insect (larvae and pupae) morphology makes morphological identification more difficult. The majority of the damage is caused by adults and fungi that serve as nutrition for their larvae. They destroy plant trunks, branches, and twigs, affecting plant transport tissues in both weak and healthy plants. An accurate, efficient, and economical molecular identification technique for X. compactus not restricted by professional taxonomic knowledge is necessary. In the present study, a molecular identification tool based on the mitochondrial DNA gene, cytochrome C oxidase subunit I (COI) was developed. A species-specific COI (SS-COI) PCR assay was designed to identify X. compactus regardless of the developmental stage. Twelve scolytines commonly found in eastern China, namely Xylosandrus compactus, X. crassiusculus, X. discolor, X. germanus, X. borealis, X. amputates, X. eupatorii, X. mancus, Xyleborinus saxesenii, Euwallacea interjectus, E. fornicatus, and Acanthotomicus suncei, were included in the study. Additionally, specimens of X. compactus from 17 different areas in China, as well as a specimen collected from the United Stated, were also analyzed. Results demonstrated the accuracy and high efficiency of the assay, regardless of the developmental stage or the type of specimen. These features provide a good application prospect for fundamental departments and can be used to prevent the harmful consequences of the spread of X. compactus.

Chromosome-level genome assemblies of two parasitoid biocontrol wasps reveal the parthenogenesis mechanism and an associated novel virus

BACKGROUND

Biocontrol is a key technology for the control of pest species. Microctonus parasitoid wasps (Hymenoptera: Braconidae) have been released in Aotearoa New Zealand as biocontrol agents, targeting three different pest weevil species. Despite their value as biocontrol agents, no genome assemblies are currently available for these Microctonus wasps, limiting investigations into key biological differences between the different species and strains.

METHODS AND FINDINGS

Here we present high-quality genomes for Microctonus hyperodae and Microctonus aethiopoides, assembled with short read sequencing and Hi-C scaffolding. These assemblies have total lengths of 106.7 Mb for M. hyperodae and 129.2 Mb for M. aethiopoides, with scaffold N50 values of 9 Mb and 23 Mb respectively. With these assemblies we investigated differences in reproductive mechanisms, and association with viruses between Microctonus wasps. Meiosis-specific genes are conserved in asexual Microctonus, with in-situ hybridisation validating expression of one of these genes in the ovaries of asexual Microctonus aethiopoides. This implies asexual reproduction in these Microctonus wasps involves meiosis, with the potential for sexual reproduction maintained. Investigation of viral gene content revealed candidate genes that may be involved in virus-like particle production in M. aethiopoides, as well as a novel virus infecting M. hyperodae, for which a complete genome was assembled.

CONCLUSION AND SIGNIFICANCE

These are the first published genomes for Microctonus wasps which have been deployed as biocontrol agents, in Aotearoa New Zealand. These assemblies will be valuable resources for continued investigation and monitoring of these biocontrol systems. Understanding the biology underpinning Microctonus biocontrol is crucial if we are to maintain its efficacy, or in the case of M. hyperodae to understand what may have influenced the significant decline of biocontrol efficacy. The potential for sexual reproduction in asexual Microctonus is significant given that empirical modelling suggests this asexual reproduction is likely to have contributed to biocontrol decline. Furthermore the identification of a novel virus in M. hyperodae highlights a previously unknown aspect of this biocontrol system, which may contribute to premature mortality of the host pest. These findings have potential to be exploited in future in attempt to increase the effectiveness of M. hyperodae biocontrol.

A new species of Scolytus Geoffroy (Coleoptera, Curculionidae, Scolytinae) from Yunnan, China

Scolytus unicornis, a new species of Scolytus Geoffroy from Yunnan, China, is described and illustrated. Three DNA barcoding sequences (COI, 28S, CAD) of this species are provided. The new species is distinguished from other Asian Scolytus species by the longitudinal wrinkles on the frons only in the area below the eyes, a large median spine situated in the middle of the ventrite 2 base, and female frons with a slightly raised blunt tubercle above the epistoma.

Cloning and expression of the mitochondrial cytochrome c oxidase subunit II gene in Sitophilus zeamais and interaction mechanism with allyl isothiocyanate

In the United States, allyl isothiocyanate (AITC) has been registered as an insecticide, bactericide, and nematicide. And it has been confirmed that AITC has significant insecticidal activities against four stored product pests including Sitophilus zeamais Mostchulky (Coleoptera: Curculionidae). This study aimed to verify the mechanism of action of AITC on cytochrome c oxidase core subunits II in S. zeamais. Enzyme - catalyzed reactions and Fourier transform infrared spectrometer (FTIR) analysis revealed that the expressed COX II proteins could competitively bind and inhibit the activity of COX II. Furthermore, molecular docking results showed that a sulfur atom of AITC could form a 2.9 Å hydrogen bond with Ile-30, having a binding energy of -2.46 kcal/mol.

De novo transcriptome assemblies of five major European oilseed rape insect pests

OBJECTIVE

Insect pests can cause severe losses in oilseed rape yields across Europe. Genomic and transcriptomic information is very limited for these insects. The aim of our study was to provide transcriptomic resources on several oilseed rape herbivores that will support research into their biology and help develop new methods of sustainable pest management.

DATA

Transcriptomes for larval stages of five major European pest species were de novo assembled by Trinity assembler. Total number of transcripts ranged from 112,247 for Ceutorhynchus pallidactylus to 225,110 for Ceutorhyncus napi. Intermediate numbers of 140,588, 140,998 and 144,504, were found for Psylliodes chrysocephala, Dasineura brassicae, and Brassicogethes aeneus, respectively. Bench-marking universal single-copy orthologues analyses for each dataset indicated high degree of completeness for all five species. The transcriptomes extend the list of genomic data on insect larvae that constitute major pests of oilseed rape. The data provide information on larval physiology and form a basis to develop highly specific RNA interference-based plant protection.

Identification and sex-specific expression of chemosensory genes in the antennal transcriptomes of Pachyrhinus yasumatsui (Coleoptera: Curculionidae)

Pachyrhinus yasumatsui Kono et Morimoto is a major pest of Chinese jujube, which is widespread in northern China and causes severe economic losses in the jujube industry. Chemosensory genes play crucial roles in insect behaviors. Currently, little is known about chemosensory genes in P. yasumatsui. In the present study, antennal transcriptomes of female and male adult P. yasumatsui were annotated. In total, 113 genes involved in chemosensory functions were identified, including 41 odorant receptors, 28 odorant-binding proteins, 16 ionotropic receptors, 15 chemosensory proteins, 9 gustatory receptors, and 4 sensory neuron membrane proteins. Subsequently, the phylogenetic analyses of these olfactory-related proteins in P. yasumatsui were conducted using multiple sequence alignment. Furthermore, sex-specific expression levels of 113 genes were analyzed based on fragments per kilobase of transcript per million mapped reads (FPKM). Then, the quantitative real-time PCR (RT-qPCR) was used to quantify gene expression profiles of 28 P. yasumatsui OBPs (PyasOBPs) and 15 CSPs (PyasCSPs). The results revealed that 20 PyasOBPs and 13 PyasCSPs exhibited significantly higher expression in the antennae than in the bodies, suggesting that they might have functions in olfaction. Moreover, some OBPs and CSPs (PyasOBP6, PyasOBP7, PyasOBP16, PyasOBP21, and PyasCSP4) exhibited female-biased expression, indicating that they might take part in several female-specific behaviors. This study will promote the understanding of olfactory mechanism in P. yasumatsui, and our findings lay the groundwork for developing environmentally friendly pest management measures.

Insight into weevil biology from a reference quality genome of the boll weevil, Anthonomus grandis grandis Boheman (Coleoptera: Curculionidae)

The boll weevil, Anthonomus grandis grandis Boheman, is one of the most historically impactful insects due to its near destruction of the US cotton industry in the early 20th century. Contemporary efforts to manage this insect primarily use pheromone baited traps for detection and organophosphate insecticides for control, but this strategy is not sustainable due to financial and environmental costs. We present a high-quality boll weevil genome assembly, consisting of 306 scaffolds with approximately 24,000 annotated genes, as a first step in the identification of gene targets for novel pest control. Gene content and transposable element distribution are similar to those found in other Curculionidae genomes; however, this is the most contiguous and only assembly reported to date for a member in the species-rich genus Anthonomus. Transcriptome profiles across larval, pupal, and adult life stages led to identification of several genes and gene families that could present targets for novel control strategies.

Odorant-Binding and Chemosensory Proteins in Anthonomus eugenii (Coleoptera: Curculionidae) and Their Tissue Expression

The pepper weevil Anthonomus eugenii is one of the most damaging pests to the pepper crop. To offer alternative management strategies to insecticides, several studies have identified the semiochemicals that are involved in the pepper weevil's aggregation and mating behavior; however, there is no information on its perireceptor molecular mechanism, to date. In this study, bioinformatics tools were used to functionally annotate and characterize the A. eugenii head transcriptome and their probable coding proteins. We identified twenty-two transcripts belonging to families related to chemosensory processes, seventeen corresponding to odorant-binding proteins (OBP), and six to chemosensory proteins (CSP). All results matched with closely related Coleoptera: Curculionidae homologous proteins. Likewise, twelve OBP and three CSP transcripts were experimentally characterized by RT-PCR in different female and male tissues. The results by sex and tissue display the different expression patterns of the AeugOBPs and AeugCSPs; some are present in both sexes and all tissues, while others show expressions with higher specificity, which suggests diverse physiological functions in addition to chemo-detection. This study provides information to support the understanding of odor perception in the pepper weevil.

Complete Mitochondrial Genome of Scolytoplatypodini Species (Coleoptera: Curculionidae: Scolytinae) and Phylogenetic Implications

The complete mitochondrial genomes (mitogenomes) of beetles in the tribe Scolytoplatypodini (genus Scolytoplatypus) were sequenced and annotated. These included Scolytoplatypus raja (15,324 bp), Scolytoplatypus sinensis (15,394 bp), Scolytoplatypus skyliuae (15,167 bp), and Scolytoplatypus wugongshanensis (15,267 bp). The four mitogenomes contained 37 typical genes, including 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), and 2 ribosomal RNA genes (rRNAs). The gene orientation and arrangement of the four mitogenomes were similar to other Coleoptera mitogenomes. PCGs mostly started with ATN and terminated with TAA. The Ka/Ks ratio of 13 PCGs in the four species revealed that cox1 had the slowest evolutionary rate and atp8 and nad6 had a higher evolutionary rate. All tRNAs had typical cloverleaf secondary structures, but trnS1 lacked dihydrouridine arm. Partial tRNAs lost the discriminator nucleotide. The trnY did not possess the discriminator nucleotide and also lost three bases, showing a special amino-acyl arm. Bayesian inference (BI) and maximum likelihood (ML) methods were conducted for phylogenetic analyses using 13 PCGs. Scolytoplatypodini was clustered with Hylurgini and Hylastini, and the monophyly of Scolytoplatypodini was supported. The four newly sequenced mitogenomes increase understanding of the evolutionary relationships of Scolytoplatypodini and other Scolytinae species.

Specific Enriched Acinetobacter in Camellia Weevil Gut Facilitate the Degradation of Tea Saponin: Inferred from Bacterial Genomic and Transcriptomic Analyses

Beneficial gut bacteria can enhance herbivorous arthropod adaptation to plant secondary compounds (PSMs), and specialist herbivores provide excellent examples of this. Tea saponin (TS) of Camellia oleifera is triterpenoids toxic to seed-feeding weevil pest, Curculio chinensis (CW). Previous studies disclosed that Acinetobacter, which was specific enriched in the CW's gut, was involved in helping CW evade TS toxicity of C. oleifera. However, it is still not clear whether Acinetobacter is associated with other anti-insect compounds, and the molecular mechanism of Acinetobacter degradation of TS has not been clarified. To address these questions, we explored the relationship between host plant toxin content and Acinetobacter of CW gut bacteria. Results demonstrated that TS content significantly affected the CW gut microbiome structure and enriched bacteria functional for TS degradation. We further isolated Acinetobacter strain and conducted its genome and transcriptome analyses for bacterial characterization and investigation on its role in TS degradation. Biological tests were carried out to verify the ability of the functional bacterium within CW larvae to detoxify TS. Our results showed that TS-degrading bacteria strain (Acinetobacter sp. AS23) genome contains 47 genes relating to triterpenoids degradation. The AS23 strain improved the survival rate of CW larvae, and the steroid degradation pathway could be the key one for AS23 to degrade TS. This study provides the direct evidence that gut bacteria mediate adaptation of herbivorous insects to phytochemical resistance. IMPORTANCE Microorganism is directly exposed to the plant toxin environment and play a crucial third party in herbivores gut. Although previous studies have proved the existence of gut bacteria that help CWs degrade TS, the specific core flora and its function have not been explored. In this study, we investigated the correlation between the larva gut microbiome and plant secondary metabolites. Acinetobacter genus was the target flora related to TS degradation. There were many terpenoids genes in Acinetobacter sp. AS23 genome. Results of transcriptome analysis and biological tests suggested that steroid degradation pathway be the key pathway of AS23 to degrade TS. This study not only provides direct evidence that gut microbes mediate the rapid adaptation of herbivorous insects to phytochemical resistance, but also provides a theoretical basis for further research on the molecular mechanism of intestinal bacteria cooperating with pests to adapt to plant toxins.

Building resiliency in conifer forests: Interior spruce crosses among weevil resistant and susceptible parents produce hybrids appropriate for multi-trait selection

Tree planting programs now need to consider climate change increasingly, therefore, the resistance to pests plays an essential role in enabling tree adaptation to new ranges through tree population movement. The weevil Pissodes strobi (Peck) is a major pest of spruces and substantially reduces lumber quality. We revisited a large Interior spruce provenance/progeny trial (2,964 genotypes, 42 families) of varying susceptibility, established in British Columbia. We employed multivariate mixed linear models to estimate covariances between, and genetic control of, juvenile height growth and resistance traits. We performed linear regressions and ordinal logistic regressions to test for impact of parental origin on growth and susceptibility to the pest, respectively. A significant environmental component affected the correlations between resistance and height, with outcomes dependent on families. Parents sourced from above 950 m a.s.l. elevation negatively influenced host resistance to attacks, probably due to higher P. engelmannii proportion. For the genetic contribution of parents sourced from above 1,200 m a.s.l., however, we found less attack severity, probably due to a marked mismatch in phenologies. This clearly highlights that interspecific hybrid status might be a good predictor for weevil attacks and delineates the boundaries of successful spruce population movement. Families resulting from crossing susceptible parents generally showed fast-growing trees were the most affected by weevil attacks. Such results indicate that interspecific 'hybrids' with a higher P. glauca ancestry might be genetically better equipped with an optimized resource allocation between defence and growth and might provide the solution for concurrent improvement in resistance against weevil attacks, whilst maintaining tree productivity.

Natural allelic variation confers high resistance to sweet potato weevils in sweet potato

Sweet potato (Ipomoea batatas L.) is a major root crop worldwide. Sweet potato weevils (SPWs) pose one of the most significant challenges to sweet potato production in tropical and subtropical regions, causing deleterious economic and environmental effects. Characterizing the mechanisms underlying natural resistance to SPWs is therefore crucial; however, the genetic basis of host SPW resistance (SPWR) remains unclear. Here we obtained two sweet potato germplasm with high SPWR and, by map-based cloning, revealed two major SPW-resistant genes-SPWR1 and SPWR2-that are important regulators of natural defence against SPWs. The SPW-induced WRKY transcriptional factor SPWR1 directly activates the expression of SPWR2, and SPWR2, the conserved dehydroquinate synthase, promotes the accumulation of quinate derivative metabolites that confer SPWR in sweet potato. Generally, our results provide new insights into the molecular mechanism underlying sweet potato-SPW interactions and will aid future efforts to achieve eco-friendly SPW management.

Tandemly duplicated genes encoding polygalacturonase inhibitors are associated with bruchid (Callosobruchus chinensis) resistance in moth bean (Vigna aconitifolia)

Bruchids are stored-grain insect pests responsible for serious seed loss in legume crops. A previous study using an F₂ population (F2OA) derived from a cross between wild moth-bean (Vigna aconitifolia [Jacq.] Maréchal) accession TN67 (resistant) and cultivated moth-bean accession ICPMO056 (susceptible) revealed that resistance to the azuki bean weevil (Callosobruchus chinensis L.) in TN67 was regulated by a single gene located in the major quantitative trait locus-qVacBrc2.1. In this study, qVacBrc2.1 was finely mapped and candidate genes in this locus were identified using F2OA and another large F₂ population (F2NB) derived from the cross mentioned previously. In contrast to the previous study, segregation analysis in the F2NB population revealed that resistance against this pest was controlled by two genes. Furthermore, the addition of novel markers to qVacBrc2.1 and reanalysis of the QTL in the F2OA population demonstrated that qVacBrc2.1 constituted two linked QTLs-qVacBrc2.1-A and qVacBrc2.1-B. The presence of qVacBrc2.1-B was verified using the population F2NB. Comparative genomics using three Vigna spp. strongly suggested the presence of two tandemly duplicated genes, VacPGIP1 and VacPGIP2, which encoded polygalacturonase inhibitors (polygalacturonase-inhibiting proteins) as the candidates for conferring resistance, but only VacPGIP1 could be successfully cloned and sequenced. The alignment of VacPGIP1 coding sequences of TN67 and ICPMO056 revealed eight single nucleotide polymorphisms, three of which altered the amino-acid sequence of the predicted domains of polygalacturonase inhibitors in ICPMO056. Overall, these findings indicate that VacPGIP1 and VacPGIP2 regulated C. chinensis resistance in TN67.

Discovery and functional characterization of novel cotton promoters with potential application to pest control

pGhERF105 and pGhNc-HARBI1 promoters are highly responsive to CBW infestation and exhibit strong activity in vegetative and reproductive tissues, increasing their potential application in GM crop plants for pest control. The main challenge to cotton (Gossypium hirsutum) crop productivity is the constant attack of several pests, including the cotton boll weevil (CBW, Anthonomus grandis), which uses cotton floral buds for feeding and egg-laying. The endophytic nature of the early developmental stages of CBW makes conventional pesticide-based control poorly efficient. Most biotechnological assets used for pest control are based on Bacillus thurigiensis insecticidal Cry toxins or the silencing of insect-pest essential genes using RNA-interference technology. However, suitable plant promoter sequences are required to efficiently drive insecticidal molecules to the target plant tissue. This study selected the Ethylene Responsive Factor 105 (GhERF105) and Harbinger transposase-derived nuclease (GhNc-HARBI1) genes based on available transcriptome-wide data from cotton plants infested by CBW larvae. The GhERF105 and GhNc-HARBI1 genes showed induction kinetics from 2 to 96 h under CBW's infestation in cotton floral buds, uncovering the potential application of their promoters. Therefore, the promoter regions (1,500 base pairs) were assessed and characterized using Arabidopsis thaliana transgenic plants. The pGhERF105 and pGhNc-HARBI1 promoters showed strong activity in plant vegetative (leaves and roots) and reproductive (flowers and fruits) tissues, encompassing higher GUS transcriptional activity than the viral-constitutive Cauliflower Mosaic Virus 35S promoter (pCaMV35S). Notably, pGhERF105 and pGhNc-HARBI1 promoters demonstrated more efficiency in driving reporter genes in flowers than other previously characterized cotton flower-specific promoters. Overall, the present study provides a new set of cotton promoters suitable for biotechnological application in cotton plants for pest resistance.

Identification of miRNAs in Response to Sweet Potato Weevil (Cylas formicarius) Infection by sRNA Sequencing

The sweet potato weevil (Cylas formicarius) is an important pest in the growing and storage of sweet potatoes. It is a common pest in the sweet potato production areas of southern China, causing serious harm to the development of the sweet potato industry. For the existing cultivars in China and abroad, there is no sweet potato variety with complete resistance to the sweet potato weevil. Thus, understanding the regulation mechanisms of sweet potato weevil resistance is the prerequisite for cultivating sweet potato varieties that are resistant to the sweet potato weevil. However, very little progress has been made in this field. In this study, we inoculated adult sweet potato weevils into sweet potato tubers. The infected sweet potato tubers were collected at 0, 24, 48, and 72 h. Then, a miRNA library was constructed for Eshu 6 and Guang 87 sweet potato tubers infected for different lengths of time. A total of 407 known miRNAs and 298 novel miRNAs were identified. A total of 174 differentially expressed miRNAs were screened out from the known miRNAs, and 247 differentially expressed miRNAs were screened out from the new miRNAs. Moreover, the targets of the differentially expressed miRNAs were predicted and their network was further investigated through GO analysis and KEGG analysis using our previous transcriptome data. More importantly, we screened 15 miRNAs and their target genes for qRT-PCR verification to confirm the reliability of the high-throughput sequencing data, which indicated that these miRNAs were detected and most of the expression results were consistent with the sequencing results. These results provide theoretical and data-based resources for the identification of miRNAs in response to sweet potato weevil infection and an analysis of the molecular regulatory mechanisms involved in insect resistance.

Analysis of Genetic Diversity of Banana Weevils (Cosmopolites sordidus) (Coleoptera: Curculionidae) Using Transcriptome-Derived Simple Sequence Repeat Markers

The banana weevil, Cosmopolites sordidus (Germar) (Coleoptera: Curculionidae) is an economically important insect pest of bananas. It causes up to 100% yield losses and substantial lifespan reduction in bananas. Advances in genomics, proteomics, and sequencing technologies have provided powerful pathways to genotyping disastrous pests such as C. sordidus. However, such technologies are often not available to the majority of rural subtropical African banana growers and pest control managers. This study was therefore motivated by the need to create cheap and easily accessible C. sordidus genotyping methods that could be deployed by banana pest control managers to the benefit of C. sordidus control programs in the tropics where such advanced technologies are not readily accessible. We used an in-house C. sordidus transcriptome from the an-ongoing study from which we mined an array of simple sequence repeat (SSR) markers. Of these, six highly polymorphic transcriptome-derived SSR markers were used to successfully genotype within and among banana weevil population genetic diversity of 12 C. sordidus populations collected from four banana-growing agro-ecological zones (AEZs) in Uganda. The developed transcriptome-derived SSR markers can be used by researchers in population genetics for characterization of the C. sordidus and identification of new genes that are linked to traits of particular interest. The significant genetic diversity revealed in C. sordidus provides pertinent information for integrated pest management strategies.

Population Structure and Genetic Diversity of the Pepper Weevil (Coleoptera: Curculionidae) Using the COI Barcoding Region

The pepper weevil Anthonomus eugenii Cano (Coleoptera: Curculionidae) is a pest of economic importance for Capsicum species pepper in North America that attacks the reproductive structures of the plant. The insect is distributed across Mexico, the United States, and the Caribbean, and is occasionally found during the pepper growing season in southern Ontario, Canada. Continuous spread of the insect to new areas is partially the result of global pepper trade. Here, we describe the genetic diversity of the pepper weevil using the mitochondrial COI barcoding region across most of its geographic range. In this study, 44 (H1-H44) highly similar haplotypes were identified, the greatest number of haplotypes and haplotype diversity were observed among specimens from its native Mexico, followed by specimens from the United States. Unlike Mexico, a low haplotype diversity was found among specimens from Canada, the Dominican Republic, Italy, and the Netherlands. Out of these 44 haplotypes, 29 are reported for the first time. Haplotype diversity in the Canadian population suggests either multiple and continuous introductions of the pepper weevil into this area or a single introduction of genetically diverse individuals. We discuss the importance of such population genetic data in tailoring pepper weevil management programs, using Canada as an example.

Mutation in the Cadherin Gene Is a Key Factor for Pink Bollworm Resistance to Bt Cotton in China

Transgenic crops producing Bacillus thuringiensis (Bt) toxins are widely planted for insect control, but their efficacy may decrease as insects evolve resistance. Understanding the genetic basis of insect resistance is essential for developing an integrated strategy of resistance management. To understand the genetic basis of resistance in pink bollworm (Pectinophora gossypiella) to Bt cotton in the Yangtze River Valley of China, we conducted an F2 screening for alleles associated with resistance to the Bt (Cry1Ac) protein for the first time. A total of 145 valid single-paired lines were screened, among which seven lines were found to carry resistance alleles. All field parents in those seven lines carried recessive resistance alleles at the cadherin locus, including three known alleles, r1, r13 and r15, and two novel alleles, r19 and r20. The overall frequency of resistance alleles in 145 lines was 0.0241 (95% CI: 0.0106-0.0512). These results demonstrated that resistance was rare and that recessive mutation in the cadherin gene was the primary mechanism of pink bollworm resistance to Bt cotton in the Yangtze River Valley of China, which will provide a scientific basis for implementing targeted resistance management statics of pink bollworm in this region.

Metabolomics and transcriptomics of pheromone biosynthesis in an aggressive forest pest Ips typographus

Eurasian spruce bark beetle, Ips typographus, is a destructive pest in spruce forests. The ability of I. typographus to colonise host trees depends on its massive aggregation behaviour mediated by aggregation pheromones, consisting of 2-methyl-3-buten-2-ol and cis-verbenol. Other biologically active compounds such as ipsdienol and verbenone have also been detected in the beetle. Biosynthesis of 2-methyl-3-buten-2-ol and ipsdienol de novo from mevalonate and that of cis-verbenol from α-pinene sequestrated from the host have been reported in preliminary studies. However, knowledge on the molecular mechanisms underlying pheromone biosynthesis in this pest is currently limited. In this study, we performed metabolomic and differential gene expression (DGE) analysis for the pheromone-producing life stages of I. typographus. The highest amounts of 2-methyl-3-buten-2-ol (238 ng/gut) and cis-verbenol (23 ng/gut) were found in the fed male gut (colonisation stage) and the immature male gut (early stage), respectively. We also determined the amount of verbenyl oleate (the possible storage form of cis-verbenol), a monoterpenyl fatty acid ester, to be approximately 1604 ng/mg in the immature stage in the beetle body. DGE analysis revealed possible candidate genes involved in the biosynthesis of the quantified pheromones and related compounds. A novel hemiterpene-synthesising candidate isoprenyl-di-phosphate synthase Ityp09271 gene proposed for 2-methyl-3-buten-2-ol synthesis was found to be highly expressed only in the fed male beetle gut. Putative cytochrome P450 genes involved in cis/trans-verbenol synthesis and an esterase gene Ityp11977, which could regulate verbenyl oleate synthesis, were identified in the immature male gut. Our findings from the molecular analysis of pheromone-producing gene families are the first such results reported for I. typographus. With further characterisation of the identified genes, we can develop novel strategies to disrupt the aggregation behaviour of I. typographus and thereby prevent vegetation loss.

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.

Silencing of vitellogenin gene contributes to the promise of controlling red palm weevil, Rhynchophorus ferrugineus (Olivier)

Red palm weevil [Rhynchophorus ferrugineus (Olivier)], is native to South Asia and expanding its distribution range globally. Recent invasions of red palm weevil around the world, including Saudi Arabia, has become a global constraint for the production of palm species. Although, several control measures have been tested, none of them seemed successful against this invasive species. Therefore, we focused on silencing the reproduction control gene vitellogenin (Vg) based on RNA interference (RNAi) strategy for its possible application in the management of R. ferrugineus. The Vg is a major yolk protein precursor critical for oogenesis. To do this, fat body transcriptome of R. ferrugineus female adults was sequenced, which provided partial Vg gene transcript (FPKM 5731.60). A complete RfVg gene transcript of 5504 bp encoding 1787 amino acids was then sequenced using RCAE-PCR strategy and characterized. Phylogenetic analysis suggested that RfVg has closer ancestry to the coleopteran insects. The RfVg-based RNAi significantly suppressed the expressions of Vg gene. The 15, 20 and 25 days post-injection periods suppressed Vg expressions by 95, 96.6 and 99%, respectively. The suppressed Vg expressions resulted in the dramatic failure of Vg protein expression, which caused atrophied ovaries or no oogenesis and ultimately eggs were not hatched. These results suggest that knockdown of Vg gene involved in R. ferrugineus reproduction is a promising target for RNAi-based management of R. ferrugineus.

Function of four mitochondrial genes in fumigation lethal mechanisms of Allyl Isothiocyanate against Sitophilus zeamais adults

Allyl isothiocyanate (AITC) is a promising alternative to chemical fumigants, and mitochondrial dysfunction has been proposed to play a crucial role in its lethal mechanisms; however, the specific lethal mechanisms of AITC remain unknown. Four mitochondrial electron transport chain genes, nd5, nd6, cox1, and cox5, were selected from adult Sitophilus zeamais and processed with RNA interference experiments. Then, the biochemical and biophysical effects were compared between double-stranded RNA (dsRNA)-mediated insects and wild-type insects after AITC fumigation at the concentration of LC₅₀ values. The bioactivity of AITC against dsnd6-mediated insects increased, while the bioactivity against dcox1-mediated insects decreased. Compared with the wild-type insects, the increase of reactive oxygen species (ROS) levels by AITC in mitochondria from dsnd6-mediated insects increased by 18.95%, while that of dscox1-mediated insects decreased by 27.45%. The effects of AITC on mRNA expression levels of detoxifying enzymes including CAT (down-regulation effect) and CuZnSOD (overexpression effect) partly recovered in the dsnd5-mediated insects, while a greater effect was observed for dscox1-mediated insects. Molecular docking results indicated that ASN511 at the cox1 subunit was the binding site of AITC by one hydrogen bond, with a bond distance of 2.1 Å. These findings provide insight for further applications of AITC and could provide a novel strategy to investigate lethal mechanisms of insecticides.

Field-evolved resistance to beta-cyfluthrin in the boll weevil: Detection and characterization

BACKGROUND

Insecticide resistance in arthropods is an inherited trait that has become a major cause of insect pest control failure. Monitoring the level of susceptibility and characterization of the type of resistance of key pest species aims to determine the risk of resistance selection in time to take action to mitigate control failures. Seven populations of the boll weevil, Anthonomus grandis grandis, collected from cotton fields in the Semiarid and Cerrado areas of Brazil, were screened for their resistance to malathion and beta-cyfluthrin, insecticides widely recommended for control of boll weevil and other pests.

RESULTS

The levels of adult mortality were variable for beta-cyfluthrin (0-82%) but invariant (100%) for malathion. Bioassays of concentration-mortality were used to determine lethal concentrations (LCs) for each insecticide. The LC-values corroborate the lack of resistance to field rates of malathion but high levels of resistance to beta-cyfluthrin from 62.7- to 439.7-fold. Weevils resistant to beta-cyfluthrin were found through genome sequencing to possess a kdr mutation through the L1014F substitution in the voltage gated-sodium channel gene.

CONCLUSIONS

This study found boll weevil resistance to beta-cyfluthrin to be not mediated by carboxylesterases, but with cross-resistance to DDT and carbaryl, and kdr mutation as the major mechanism of the resistance in our samples. Caution is recommended in further use of beta-cyfluthrin against boll weevil due to potential resistance. Monitoring studies using other boll weevil populations are recommended to determine the geographic pattern and extent of pyrethroid resistance. © 2021 Society of Chemical Industry.

A highly-contiguous genome assembly of the Eurasian spruce bark beetle, Ips typographus, provides insight into a major forest pest

Conifer-feeding bark beetles are important herbivores and decomposers in forest ecosystems. These species complete their life cycle in nutritionally poor substrates and some can kill enormous numbers of trees during population outbreaks. The Eurasian spruce bark beetle (Ips typographus) can destroy >100 million m3 of spruce in a single year. We report a 236.8 Mb I. typographus genome assembly using PacBio long-read sequencing. The final phased assembly has a contig N50 of 6.65 Mb in 272 contigs and is predicted to contain 23,923 protein-coding genes. We reveal expanded gene families associated with plant cell wall degradation, including pectinases, aspartyl proteases, and glycosyl hydrolases. This genome sequence from the genus Ips provides timely resources to address questions about the evolutionary biology of the true weevils (Curculionidae), one of the most species-rich animal families. In forests of today, increasingly stressed by global warming, this draft genome may assist in developing pest control strategies to mitigate outbreaks.

Epicuticular Wax Rice Mutants Show Reduced Resistance to Rice Water Weevil (Coleoptera: Curculionidae) and Fall Armyworm (Lepidoptera: Noctuidae)

Plant structural traits can act as barriers for herbivore attachment, feeding, and oviposition. In particular, epicuticular waxes (EWs) on the aerial surfaces of many land plants offer protection from biotic and abiotic stresses. In rice (Oryza sativa L.), mutations that reduce EWs have been previously reported. However, whether such mutations affect rice water weevil (Lissorhoptrus oryzophilus Kuschel) and fall armyworm (Spodoptera frugiperda Smith) performance has not been investigated yet. These pests cause significant economic problems in important rice-producing areas of the United States. The aim of our study was to characterize the EWs of EW mutants and wild-type rice plants by gas chromatography-mass spectrometry and compare the resistance of mutant and wild-type plants against rice water weevil and fall armyworm. We hypothesized that mutants with reduced EWs would have weaker resistance to pests than wild-type plants. Three mutant lines (6-1A, 7-17A, and 11-39A) and their wild-type parent (cv. 'Sabine') were used to test this hypothesis. Levels of EWs were significantly lower in mutant lines than in the wild-type, and qualitative differences in EW composition were also observed. Reduction in EWs significantly affected performance of insects in experiments conducted under greenhouse conditions. Experiments with rice water weevils were conducted in arenas in which females were given a choice of the mutants and the wild-type for oviposition. Number of first instars emerging from the three EW mutants (an indication of oviposition preference) was higher on the three EW mutants than on wild-type plants with normal wax levels. Similarly, in no-choice experiments using whole plants or detached leaves, weight gains of armyworms on leaves were higher on the mutant lines than on the wild-type. These results indicate that EW traits are involved in rice resistance to weevils and armyworms. Understanding the plant traits that contribute to resistance to rice pests will be helpful for the development of resistant varieties for reducing pest insect damage.

Function of mevalonate pathway genes in the synthesis of frontalin in Chinese white pine beetle, Dendroctonus armandi (curculionidae: Scolytinae)

The Chinese white pine beetle (Dendroctonus armandi Tsai and Li) is a significant pest of pine forests in the Qinling and Bashan Mountains of China. Adult males commonly produce frontalin using precursors synthesized through the mevalonate pathway, which is regulated by juvenile hormone III (JHIII). In this study, the expression levels of mevalonate pathway genes were quantified after phloem feeding and topical application of the JHIII solution. The frontalin was quantified by gas chromatography-mass spectrometry. Both the phloem feeding and JHIII treatments produced an evident upregulation in the male gut, mainly in 3-hydroxy-3-methylglutaryl-CoA synthase (HMGS) and 3-hydroxy-3-methylglutaryl-CoA reductase (HMGR). Moreover, HMGS, HMGR, isopentenyl diphosphate isomerase, and geranyl diphosphate synthase/farnesyl diphosphate synthase were upregulated in fed and JHIII-stimulated males of D. armandi under both conditions (solitary and paired). The expression levels were higher in paired compared to solitary males. Males had higher expression levels compared with females. Correspondingly, the phloem-feeding males produced more frontalin than JHIII-treated males, and the production of frontalin was higher in paired males than in solitary males. The knockdown of mevalonate pathway genes using RNAi in vivo effectively reduced the messenger RNA level of these genes and inhibited the production of frontalin. Among them, the silencing of HMGR or HMGS genes reduced the synthesis of frontalin most significantly.

The Easter Egg Weevil (Pachyrhynchus) genome reveals syntenic patterns in Coleoptera across 200 million years of evolution

Patterns of genomic architecture across insects remain largely undocumented or decoupled from a broader phylogenetic context. For instance, it is unknown whether translocation rates differ between insect orders. We address broad scale patterns of genome architecture across Insecta by examining synteny in a phylogenetic framework from open-source insect genomes. To accomplish this, we add a chromosome level genome to a crucial lineage, Coleoptera. Our assembly of the Pachyrhynchus sulphureomaculatus genome is the first chromosome scale genome for the hyperdiverse Phytophaga lineage and currently the largest insect genome assembled to this scale. The genome is significantly larger than those of other weevils, and this increase in size is caused by repetitive elements. Our results also indicate that, among beetles, there are instances of long-lasting (>200 Ma) localization of genes to a particular chromosome with few translocation events. While some chromosomes have a paucity of translocations, intra-chromosomal synteny was almost absent, with gene order thoroughly shuffled along a chromosome. This large amount of reshuffling within chromosomes with few inter-chromosomal events contrasts with patterns seen in mammals in which the chromosomes tend to exchange larger blocks of material more readily. To place our findings in an evolutionary context, we compared syntenic patterns across Insecta in a phylogenetic framework. For the first time, we find that synteny decays at an exponential rate relative to phylogenetic distance. Additionally, there are significant differences in decay rates between insect orders, this pattern was not driven by Lepidoptera alone which has a substantially different rate.

Host-specific gene expression as a tool for introduction success in Naupactus parthenogenetic weevils

Food resource access can mediate establishment success in invasive species, and generalist herbivorous insects are thought to rely on mechanisms of transcriptional plasticity to respond to dietary variation. While asexually reproducing invasives typically have low genetic variation, the twofold reproductive capacity of asexual organisms is a marked advantage for colonization. We studied host-related transcriptional acclimation in parthenogenetic, invasive, and polyphagous weevils: Naupactus cervinus and N. leucoloma. We analyzed patterns of gene expression in three gene categories that can mediate weevil-host plant interactions through identification of suitable host plants, short-term acclimation to host plant defenses, and long-term adaptation to host plant defenses and their pathogens. This approach employed comparative transcriptomic methods to investigate differentially expressed host detection, detoxification, immune defense genes, and pathway-level gene set enrichment. Our results show that weevil gene expression responses can be host plant-specific, and that elements of that response can be maintained in the offspring. Some host plant groups, such as legumes, appear to be more taxing as they elicit a complex gene expression response which is both strong in intensity and specific in identity. However, the weevil response to taxing host plants shares many differentially expressed genes with other stressful situations, such as host plant cultivation conditions and transition to novel host, suggesting that there is an evolutionarily favorable shared gene expression regime for responding to different types of stressful situations. Modulating gene expression in the absence of other avenues for phenotypic adaptation may be an important mechanism of successful colonization for these introduced insects.

Phenotypic and molecular analyses in rice weevil, Sitophilus oryzae (Linneaus) (Coleoptera: Curculionidae): identification of a super kdr mutation, T929I, conferring resistance to deltamethrin

BACKGROUND

The rice weevil, Sitophilus oryzae (L.) (Coleoptera: Curculionidae) is a cosmopolitan pest of stored cereal grains and other commodities globally. Infestations caused by S. oryzae makes grains unsuitable for consumption, processing, and export. Deltamethrin, a synthetic pyrethroid insecticide, is widely used in major grain storages in India as a prophylactic treatment to control this pest. However, recurrent use of this insecticide had led to genetic resistance in S. oryzae, questioning its ongoing use at the current recommended concentration.

RESULTS

Dose response analysis of resistant (Delta-R) and susceptible (Lab-S) strains of S. oryzae collected from grain storages across southern India, revealed that Delta-R was 134-fold more resistant than the Lab-S at median lethal concentration (LC₅₀ ). A concentration of 180 ppm over 48 h effectively discriminated 16 resistant field populations from Lab-S with per cent resistance ranging from 8.72% to 75.86%. Exposing all the resistant populations to 1000 ppm over 48 h identified 12 populations with strongly resistant individuals and confirmed the existence of two distinct resistance phenotypes, 'weak' and 'strong' in S. oryzae. Furthermore, sequence analysis of the voltage-gated sodium channel (vgsc) gene in Delta-R identified a single target site mutation, T929I conferring resistance in S. oryzae. CAPS (Cleaved Amplified Polymorphic Sequence) marker analysis of this allele confirmed that frequency of resistance is high (up to 0.96) supporting the results of phenotypic analysis.

CONCLUSION

Both phenotype and molecular marker analyses clearly demonstrated that deltamethrin at 180 and 1000 ppm can be used to discriminate weakly and strongly resistant populations in S. oryzae, respectively. Resistance diagnostics based on the mutation, T929I, supports our phenotypic data and indicates that resistance to deltamethrin in S. oryzae is prevalent in southern parts of India, stressing the need to identify a synergist or suitable alternatives. © 2021 Society of Chemical Industry.

De novo transcriptome analysis and identification of reproduction control genes from the red palm weevil Rhynchophorus ferrugineus

Recent attacks by the red palm weevil, Rhynchophorus ferrugineus (Olivier), have become a severe problem for palm species. In present work, fat body transcriptome of adult female red palm weevil was analyzed, focusing on the identification of reproduction control genes. Transcriptome study was completed by means of next-generation sequencing (NGS) using Illumina Hiseq 2000 sequencing system. A total of 105,938,182 raw reads, 102,645,544 clean reads, and 9,238,098,960 clean nucleotides with a guanine-cytosine content of 40.31%, were produced. The processed transcriptome data resulted in 43,789 unique transcripts (with mean lengths of 1,172 bp). It was found that 20% of total unique transcripts shared up to 80%-100% sequence identity with homologous species, mainly the mountain pine beetle Dendroctonus ponderosae (59.9%) and red flour beetle Tribolium castaneum (26.9%). Nearly 25 annotated genes were predicted to be involved in red palm weevil reproduction, including five vitellogenin (Vg) transcripts. Among the five Vg gene transcripts, one was highly expressed compared with the other four (FPKM values of 1.963, 1.471, 1.028, and 1.017, respectively), and the five Vg gene transcripts were designated as RfVg, RfVg-equivalent1, RfVg-equivalent2, RfVg-equivalent3, and RfVg-equivalent4, respectively. The high expression level of RfVg verified by RT-polymerase chain reaction analysis suggested that RfVg is the primary functional Vg gene in red palm weevil. A high similarity of RfVg with other Coleopterans was also reflected in a phylogenetic tree, where RfVg was placed within the clade of the order Coleoptera. Awareness of the major genes that play critical roles in reproduction and proliferation of red palm weevil is valuable to understand their reproduction mechanism at a molecular level. In addition, for future molecular studies, the NGS dataset obtained will be useful and will promote the exploration of biotech-based control strategies against red palm weevil, a primary pest of palm trees.