Haplotype-resolved genome assembly enables gene discovery in the red palm weevil Rhynchophorus ferrugineus

Identification of Multifunctional Putative Bioactive Peptides in the Insect Model Red Palm Weevil (Rhynchophorus ferrugineus)

Innate immunity, the body's initial defense against bacteria, fungi, and viruses, heavily depends on antimicrobial peptides (AMPs), which are small molecules produced by all living organisms. Insects, with their vast biodiversity, are one of the most abundant and innovative sources of AMPs. In this study, AMPs from the red palm weevil (RPW) Rhynchophorus ferrugineus (Coleoptera: Curculionidae), a known invasive pest of palm species, were examined. The AMPs were identified in the transcriptomes from different body parts of male and female adults, under different experimental conditions, including specimens collected from the field and those reared in the laboratory. The RPW transcriptomes were examined to predict antimicrobial activity, and all sequences putatively encoding AMPs were analyzed using several machine learning algorithms available in the CAMPR3 database. Additionally, anticancer, antiviral, and antifungal activity of the peptides were predicted using iACP, AVPpred, and Antifp server tools, respectively. Physicochemical parameters were assessed using the Antimicrobial Peptide Database Calculator and Predictor. From these analyses, 198 putatively active peptides were identified, which can be tested in future studies to validate the in silico predictions. Genome-wide analysis revealed that several AMPs have predominantly emerged through gene duplication. Noticeably, we detect a newly originated defensin allele from an ancestral defensin via the deletion of two amino acids following gene duplication in RPW, which may confer an enhanced resilience to microbial infection. Our study shed light on AMP gene families and shows that high duplication and deletion rates are essential to achieve a diversity of antimicrobial mechanisms; hence, we propose the RPW AMPs as a model for exploring gene duplication and functional variations against microbial infection.

Silencing sensory neuron membrane protein RferSNMPu1 impairs pheromone detection in the invasive Asian Palm Weevil

The red palm weevil (RPW), Rhynchophorus ferrugineus (Olivier), also known as the Asian palm weevil, is an invasive pest that causes widespread damage to palm trees around the globe. As pheromone communication is crucial for their mass attack and survival on palm trees, the olfactory concept of pest control strategies has been widely explored recently. We aim to understand the molecular basis of olfaction in RPW by studying one of the key olfactory proteins in insect pheromone communication, sensory neuron membrane proteins (SNMPs). SNMPs belong to the CD36 (cluster of differentiation 36) family that perform two distinct olfactory roles in insects, either in pheromone (odorant) transfer to the odorant receptors (SNMP1) or in the pheromone clearing process (SNMP2). In this study, we performed antennal transcriptomic screening and identified six SNMPs, mapping them on the R. ferrugineus genome, and confirmed four distinct SNMPs. Both SNMP1 proteins in RPW, viz., RferSNMPu1 and RferSNMPu2, were mapped onto the same scaffold in different loci in the RPW genome. To further understand the function of these proteins, we first classified them using phylogenetic analysis and checked their tissue-specific expression patterns. Further, we measured the relative transcript abundance of SNMPs in laboratory-reared, field-collected adults and pheromone-exposure experiments, ultimately identifying RferSNMPu1 as a potential candidate for functional analysis. We mapped RferSNMPu1 expression in the antennae and found that expression patterns were similar in both sexes. We used RNAi-based gene silencing to knockdown RferSNMPu1 and tested the changes in the RPW responses to aggregation pheromone compounds, 4-methyl-5-nonanol (ferrugineol) and 4-methyl-5-nonanone (ferrugineone), and a kairomone, ethyl acetate using electroantennogram (EAG) recordings. We found a significant reduction in the EAG recordings in the RferSNMPu1 knockdown strain of adult RPWs, confirming its potential role in pheromone detection. The structural modelling revealed the key domains in the RferSNMPu1 structure, which could likely be involved in pheromone detection based on the identified ectodomain tunnels. Our studies on RferSNMPu1 with a putative role in pheromone detection provide valuable insight into understanding the olfaction in R. ferrugineus as well as in other Curculionids, as SNMPs are under-explored in terms of its functional role in insect olfaction. Most importantly, RferSNMPu1 can be used as a potential target for the olfactory communication disruption in the R. ferrugineus control strategies.

Deorphanizing an odorant receptor tuned to palm tree volatile esters in the Asian palm weevil sheds light on the mechanisms of palm tree selection

The Asian palm weevil, Rhynchophorus ferrugineus, is a tremendously important agricultural pest primarily adapted to palm trees and causes severe destruction, threatening sustainable palm cultivation worldwide. The host plant selection of this weevil is mainly attributed to the functional specialization of odorant receptors (ORs) that detect palm-derived volatiles. Yet, ligands are known for only two ORs of R. ferrugineus, and we still lack information on the mechanisms of palm tree detection. This study identified a highly expressed antennal R. ferrugineus OR, RferOR2, thanks to newly generated transcriptomic data. The phylogenetic analysis revealed that RferOR2 belongs to the major coleopteran OR group 2A and is closely related to a sister clade containing an R. ferrugineus OR (RferOR41) tuned to the non-host plant volatile and antagonist, α-pinene. Functional characterization of RferOR2 via heterologous expression in Drosophila olfactory neurons revealed that this receptor is tuned to several ecologically relevant palm-emitted odors, most notably ethyl and methyl ester compounds, but not to any of the pheromone compounds tested, including the R. ferrugineus aggregation pheromone. We did not evidence any differential expression of RferOR2 in the antennae of both sexes, suggesting males and females detect these compounds equally. Next, we used the newly identified RferOR2 ligands to demonstrate that including synthetic palm ester volatiles as single compounds and in combinations in pheromone-based mass trapping has a synergistic attractiveness effect to R. ferrugineus aggregation pheromone, resulting in significantly increased weevil catches. Our study identified a key OR from a palm weevil species tuned to several ecologically relevant palm volatiles and represents a significant step forward in understanding the chemosensory mechanisms of host detection in palm weevils. Our study also defines RferOR2 as an essential model for exploring the molecular basis of host detection in other palm weevil species. Finally, our work showed that insect OR deorphanization could aid in identifying novel behaviorally active volatiles that can interfere with weevil host-searching behavior in sustainable pest management applications.

Genome-wide identification, molecular evolution and gene expression of P450 gene family in Cyrtotrachelus buqueti

BACKGROUND

Insect Cytochrome P450 monooxygenase (CYPs or P450s) plays an important role in detoxifying insecticides, causing insect populations to develop resistance. However, the molecular functions of P450 gene family in Cyrtotrachelus buqueti genome are still lacking.

RESULTS

In this study, 71 CbuP450 genes have been identified. The amino acids length of CbuP450 proteins was between 183 aa ~ 1041 aa. They are proteins with transmembrane domains. The main component of their secondary structure is α-helix and random coils. Phylogenetic analysis showed that C. buqueti and Rhynchophorus ferrugineus were the most closely related. This gene family has 29 high-frequency codons, which tend to use A/T bases and A/T ending codons. Gene expression analysis showed that CbuP450_23 in the female adult may play an important role on high temperature resistance, and CbuP450_17 in the larval may play an important role on low temperature tolerance. CbuP450_10, CbuP450_17, CbuP450_23, CbuP450_10, CbuP450_16, CbuP450_20, CbuP450_23 and CbuP450_ 29 may be related to the regulation of bamboo fiber degradation genes in C. buqueti. Protein interaction analysis indicates that most CbuP450 proteins are mainly divided into three aspects: encoding the biosynthesis of ecdysteroids, participating in the decomposition of synthetic insecticides, metabolizing insect hormones, and participating in the detoxification of compounds.

CONCLUSIONS

We systematically analyzed the gene and protein characteristics, gene expression, and protein interactions of CbuP450 gene family, revealing the key genes involved in the stress response of CbuP450 gene family in the resistance of C. buqueti to high or low temperature stress, and identified the key CbuP450 proteins involved in important life activity metabolism. These results provided a reference for further research on the function of P450 gene family in C. buqueti.

Taxonomy, Biology, Symbionts, Omics, and Management of Rhynchophorus Palm Weevils (Coleoptera: Curculionidae: Dryophthorinae)

Palm weevils, Rhynchophorus spp., are destructive pests of native, ornamental, and agricultural palm species. Of the 10 recognized species, two of the most injurious species, Rhynchophorus ferrugineus and Rhynchophorus palmarum, both of which have spread beyond their native range, are the best studied. Due to its greater global spread and damage to edible date industries in the Middle East, R. ferrugineus has received more research interest. Integrated pest management programs utilize traps baited with aggregation pheromone, removal of infested palms, and insecticides. However, weevil control is costly, development of resistance to insecticides is problematic, and program efficacy can be impaired because early detection of infestations is difficult. The genome of R. ferrugineus has been sequenced, and omics research is providing insight into pheromone communication and changes in volatile and metabolism profiles of weevil-infested palms. We outline how such developments could lead to new control strategies and early detection tools.

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.

Metagenomic binning of PacBio HiFi data prior to assembly reveals a complete genome of Cosmopolites sordidus (Germar) (Coleopterea: Curculionidae, Dryophthorinae) the most damaging arthropod pest of bananas and plantains

PacBio HiFi sequencing was employed in combination with metagenomic binning to produce a high-quality reference genome of Cosmopolites sordidus. We compared k-mer and alignment reference based pre-binning and post-binning approaches to remove contamination. We were also interested to know if the post-binning approach had interspersed bacterial contamination within intragenic regions of Arthropoda binned contigs. Our analyses identified 3,433 genes that were composed with reads identified as of putative bacterial origins. The pre-binning approach yielded a C. sordidus genome of 1.07 Gb genome composed of 3,089 contigs with 98.6% and 97.1% complete and single copy genome and protein BUSCO scores respectively. In this article we demonstrate that in this case the pre-binning approach does not sacrifice assembly quality for more stringent metagenomic filtering. We also determine post-binning allows for increased intragenic contamination increased with increasing coverage, but the frequency of gene contamination increased with lower coverage. Future work should focus on developing reference free pre-binning approaches for HiFi reads produced from eukaryotic based metagenomic samples.

Nuclear and Mitochondrial Genome Assemblies of the Beetle, Zygogramma bicolorata, a Globally Important Biocontrol Agent of Invasive Weed Parthenium hysterophorus

Implementing a genetic-based approach to achieve the full potential of classical biocontrol programs has been advocated for decades. The availability of genome-level information brings the opportunity to scrutinize biocontrol traits for their efficacy and evolvability. However, implementation of this advocacy remains limited to few instances. Biocontrol of a globally noxious weed, Parthenium hysterophorus, by the leaf-feeding beetle, Zygogramma bicolorata, has been in place for more than four decades now, with varying levels of success. As the first step in providing genetic-based improvement to this biocontrol program, we describe the nuclear and mitochondrial assemblies of Z. bicolorata. We assembled the genome from the long-read sequence data, error corrected with high-throughput short reads and checked for contaminants and sequence duplication to produce a 936 Mb nuclear genome. With 96.5% Benchmarking Universal Single-Copy Orthologs completeness and the long terminal repeat assembly index 12.91, we present a reference-quality assembly that appeared to be repeat rich at 62.7% genome-wide and consists of 29,437 protein-coding regions. We detected signature of nuclear insertion of mitochondrial fragments in 80 nuclear positions comprising 13 kb out of 17.9 kb mitochondria genome sequence. This genome, along with its annotations, provides a valuable resource to gain further insights into the biocontrol traits of Z. bicolorata for improving the control of the invasive weed P. hysterophorus.

Genome assembly of the Korean intertidal mud-creeper Batillaria attramentaria

Batillaridae is a common gastropod family that occurs abundantly in the shallow coastal zone of the intertidal mudflats of the northwest Pacific Ocean, Australasia, and North America. In this family, Batillaria attramentaria is known for its biological invasion and colonization in estuarine and intertidal zones. It can endure and adapt the harsh intertidal conditions such as frequent temperature alteration, salinity, and air exposure. Therefore, we sequenced and assembled this Korean batillariid genome to get insight into its intertidal adaptive features. Approximately 53 Gb of DNA sequences were generated, and 863 scaffolds were assembled into a draft genome of 0.715 Gb with 97.1% BUSCO completeness value. A total of 40,596 genes were predicted. We estimated that B. attramentaria and Conus consors diverged about 230 million years ago (MYA) based on the phylogenetic analysis of closely related gastropod species. This genome study sets the footstep for genomics studies among native and introduced Batillaria populations and the Batillaridae family members.

Omics in the Red Palm Weevil Rhynchophorus ferrugineus (Olivier) (Coleoptera: Curculionidae): A Bridge to the Pest

The red palm weevil (RPW), Rhynchophorus ferrugineus (Coleoptera: Curculionidae), is the most devastating pest of palm trees worldwide. Mitigation of the economic and biodiversity impact it causes is an international priority that could be greatly aided by a better understanding of its biology and genetics. Despite its relevance, the biology of the RPW remains poorly understood, and research on management strategies often focuses on outdated empirical methods that produce sub-optimal results. With the development of omics approaches in genetic research, new avenues for pest control are becoming increasingly feasible. For example, genetic engineering approaches become available once a species's target genes are well characterized in terms of their sequence, but also population variability, epistatic interactions, and more. In the last few years alone, there have been major advances in omics studies of the RPW. Multiple draft genomes are currently available, along with short and long-read transcriptomes, and metagenomes, which have facilitated the identification of genes of interest to the RPW scientific community. This review describes omics approaches previously applied to RPW research, highlights findings that could be impactful for pest management, and emphasizes future opportunities and challenges in this area of research.

Local assembly of long reads enables phylogenomics of transposable elements in a polyploid cell line

Animal cell lines often undergo extreme genome restructuring events, including polyploidy and segmental aneuploidy that can impede de novo whole-genome assembly (WGA). In some species like Drosophila, cell lines also exhibit massive proliferation of transposable elements (TEs). To better understand the role of transposition during animal cell culture, we sequenced the genome of the tetraploid Drosophila S2R+ cell line using long-read and linked-read technologies. WGAs for S2R+ were highly fragmented and generated variable estimates of TE content across sequencing and assembly technologies. We therefore developed a novel WGA-independent bioinformatics method called TELR that identifies, locally assembles, and estimates allele frequency of TEs from long-read sequence data (https://github.com/bergmanlab/telr). Application of TELR to a ∼130x PacBio dataset for S2R+ revealed many haplotype-specific TE insertions that arose by transposition after initial cell line establishment and subsequent tetraploidization. Local assemblies from TELR also allowed phylogenetic analysis of paralogous TEs, which revealed that proliferation of TE families in vitro can be driven by single or multiple source lineages. Our work provides a model for the analysis of TEs in complex heterozygous or polyploid genomes that are recalcitrant to WGA and yields new insights into the mechanisms of genome evolution in animal cell culture.

Chromosome-scale assembly of the yellow mealworm genome

Background: The yellow mealworm beetle, Tenebrio molitor, is a promising alternative protein source for animal and human nutrition and its farming involves relatively low environmental costs. For these reasons, its industrial scale production started this century. However, to optimize and breed sustainable new T. molitor lines, the access to its genome remains essential. Methods: By combining Oxford Nanopore and Illumina Hi-C data, we constructed a high-quality chromosome-scale assembly of T. molitor. Then, we combined RNA-seq data and available coleoptera proteomes for gene prediction with GMOVE. Results: We produced a high-quality genome with a N50 = 21.9Mb with a completeness of 99.5% and predicted 21,435 genes with a median size of 1,780 bp. Gene orthology between T. molitor and Tribolium castaneum showed a highly conserved synteny between the two coleoptera and paralogs search revealed an expansion of histones in the T. molitor genome. Conclusions: The present genome will greatly help fundamental and applied research such as genetic breeding and will contribute to the sustainable production of the yellow mealworm.

Chromosome-scale assembly of the yellow mealworm genome

Background: The yellow mealworm beetle, Tenebrio molitor, is a promising alternative protein source for animal and human nutrition and its farming involves relatively low environmental costs. For these reasons, its industrial scale production started this century. However, to optimize and breed sustainable new T. molitor lines, the access to its genome remains essential. Methods: By combining Oxford Nanopore and Illumina Hi-C data, we constructed a high-quality chromosome-scale assembly of T. molitor. Then, we combined RNA-seq data and available coleoptera proteomes for gene prediction with GMOVE. Results: We produced a high-quality genome with a N50 = 21.9Mb with a completeness of 99.5% and predicted 21,435 genes with a median size of 1,780 bp. Gene orthology between T. molitor and Tribolium castaneum showed a highly conserved synteny between the two coleoptera and paralogs search revealed an expansion of histones in the T. molitor genome. Conclusions: The present genome will greatly help fundamental and applied research such as genetic breeding and will contribute to the sustainable production of the yellow mealworm.

Microsatellite Variation in the Most Devastating Beetle Pests (Coleoptera: Curculionidae) of Agricultural and Forest Crops

Weevils, classified in the family Curculionidae (true weevils), constitute a group of phytophagous insects of which many species are considered significant pests of crops. Within this family, the red palm weevil (RPW), Rhynchophorus ferrugineus, has an integral role in destroying crops and has invaded all countries of the Middle East and many in North Africa, Southern Europe, Southeast Asia, Oceania, and the Caribbean Islands. Simple sequence repeats (SSRs), also termed microsatellites, have become the DNA marker technology most applied to study population structure, evolution, and genetic diversity. Although these markers have been widely examined in many mammalian and plant species, and draft genome assemblies are available for many species of true weevils, very little is yet known about SSRs in weevil genomes. Here we carried out a comparative analysis examining and comparing the relative abundance, relative density, and GC content of SSRs in previously sequenced draft genomes of nine true weevils, with an emphasis on R. ferrugineus. We also used Illumina paired-end sequencing to generate draft sequence for adult female RPW and characterized it in terms of perfect SSRs with 1–6 bp nucleotide motifs. Among weevil genomes, mono- to trinucleotide SSRs were the most frequent, and mono-, di-, and hexanucleotide SSRs exhibited the highest GC content. In these draft genomes, SSR number and genome size were significantly correlated. This work will aid our understanding of the genome architecture and evolution of Curculionidae weevils and facilitate exploring SSR molecular marker development in these species.

Satellitome of the Red Palm Weevil, Rhynchophorus ferrugineus (Coleoptera: Curculionidae), the Most Diverse Among Insects

The red palm weevil, Rhynchophorus ferrugineus, is the most harmful species among those pests affecting palm trees. Its impact causes important economic losses around the World. Nevertheless, the genetic information of Rh. ferrugineus is very scarce. Last year, the first genome assembly was published including a rough description of its repeatome. However, no information has been added about one of the main components of repeated DNA, the satellite DNA. Herein, we presented the characterization of the satellitome of this important species that includes 112 satellite DNA families, the largest number in an insect genome. These satellite DNA families made up around 25% of the genome while the most abundant family, RferSat01-169, alone represented 20.4%. Chromosomal location of most abundant satellite DNA families performed by fluorescence in situ hybridization showed that all of them are dispersed in the euchromatin on all chromosomes but some of them are also specifically accumulated either on the pericentromeric heterochromatic regions of all chromosomes or on specific chromosomes. Finally, the transcription of satellitome families was analyzed through Rh. ferrugineus development. It was found that 55 out of 112 satellite DNA families showed transcription, some families seemed to be transcribed across all stages while a few appeared to be stage-specific, indicating a possible role of those satellite DNA sequences in the development of this species.

Population genetic structure of the carrot weevil (Listronotus oregonensis) in North America

Population genetic studies of insect pests enhance our ability to anticipate problems in agroecosystems, such as pest outbreaks, insecticide resistance, or expansions of the host range. This study focuses on geographic distance and host plant selection as potential determinants of genetic differentiation of the carrot weevil Listronotus oregonensis, a major pest of several apiaceous crops in North America. To undertake genetic studies on this species, we assembled the first complete genome sequence for L. oregonensis. Then, we used both haplotype discrimination with mitochondrial DNA (mtDNA) and a genotyping-by-sequencing (GBS) approach to characterize the genetic population structure. A total of 220 individuals were sampled from 17 localities in the provinces of Québec, Ontario, Nova Scotia (Canada), and the state of Ohio (USA). Our results showed significant genetic differences between distant populations across North America, indicating that geographic distance represents an important factor of differentiation for the carrot weevil. Furthermore, the GBS analysis revealed more different clusters than COI analysis between Québec and Nova Scotia populations, suggesting a recent differentiation in the latter province. In contrast, we found no clear evidence of population structure associated with the four cultivated apiaceous plants tested (carrot, parsley, celery, and celeriac) using populations from Québec. This first characterization of the genetic structure of the carrot weevil contributes to a better understanding of the gene flow of the species and helps to adapt local pest management measures to better control this agricultural pest.

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.