Insect genomes: progress and challenges

The impact of microplastics on tissue-specific gene expression in the tropical house cricket, Gryllodes sigillatus

Microplastics are ubiquitous in our environment, resulting in animal exposure and consumption via food, water, and air. Animals that consume microplastics may suffer from physiological effects like immunotoxicity or mitochondrial dysfunction, but how specific tissues may differentially respond to plastic consumption is poorly understood, particularly in terrestrial insects. Here, we measured transcriptomic responses of tissues (midgut, hindgut, fat body and ovaries) to microplastic consumption in a generalist ground-dwelling insect, the tropical house cricket, Gryllodes sigillatus. Using this approach, we provide insights on how microplastics may impact specific organ systems. We generated a de novo transcriptome, a useful resource for further studies on this emerging model insect, that we then used to infer differential gene expression due to microplastic consumption in individual organs. Ingestion of microplastics elicited unique changes in gene expression depending on the tissue of focus, with notable differentially-expressed genes related to survival and stress pathways as well as those related to metabolism, immunity, and cancer.

Insect Cytochrome P450 Database: An Integrated Resource of Genetic Diversity, Evolution and Function

Insects, the most numerous and diverse group of animal species on Earth, have important interactions with humans through providing resources, transmitting diseases and damaging agricultural cultivars. Cytochrome P450 monooxygenases (P450s) are one of the most important protein families in insects implicated in the endogenous metabolism and detoxification of xenobiotics, including allelochemicals, insecticides and environmental pollutants. To better understand the evolution and function of insect P450s and support the development and application of insecticides for pest control, an integrated bioinformatics platform is highly desirable. Here, we present the Insect Cytochrome P450 database (ICPD, http://www.insectp450.net/), which contains 66,477 P450s collected from public databases and predicted from the genomes of 682 insect species using a standardised bioinformatics protocol. Phylogenetic relationships between P450 genes are constructed for each species. The structures of all P450 proteins in the database are predicted using ESMFold, then visualised using WeView. Web services, such as BLAST, homogeneous modelling and molecular docking, are provided for determining the catalytic activities of P450 proteins. The ICPD will facilitate systematic investigations of the evolution and functions of the complete insect P450 complement, and represents a powerful tool for guiding insecticide design and application.

Chromosome-level genome assembly of Dendrothrips minowai and genomic analysis highlights distinct adaptations to high polyphenols in tea plants

BACKGROUND

Dendrothrips minowai Priesner, a significant pest in tea-producing regions of Asia, particularly China, damages tea plants (Camellia sinensis) by feeding on their tender leaves rich in polyphenols. Seven assembled genomes from the order Thysanoptera are currently available.

RESULTS

This study presents the first chromosome-level genome assembly of D. minowai generated by PacBio Revio, Oxford Nanopore Technologies, MGI, and Hi-C technology. The assembled genome measures 350.11 Mb with 1269 contigs with a contig N50 of 536.34 Kb and a scaffold N50 of 16.86 Mb, organized across 19 chromosomes. A total of 16 730 protein-coding genes were identified, with 92.28% functionally annotated. The phylogenetic analysis reveals that D. minowai diverged approximately 103.2 million years ago, preceding all reported genomes of Thripidae species. Comparative genomic analysis identified 12 expanded and 172 contracted gene families of D. minowai, with expanded gene families linked to host plant metabolite processing and detoxification enzymes. Additionally, oligophagous thrips, D. minowai and Stenchaetothrips biformis, possess fewer chemosensory genes (gustatory receptors, odorant receptors, ionotropic receptors, chemosensory proteins, and odorant binding proteins) and detoxification genes (P450s, carboxyl/cholinesterases, UDP-glycosyltransferases) than polyphagous species (Frankliniella occidentalis and Thrips palmi). Interestingly, D. minowai exhibits an expansion in ABC transporter families, especially ABCG and ABCC, which is likely essential for detoxifying the high polyphenol content in tea plants.

CONCLUSION

This study provides another genome sequence for oligophagous thrips species, which enriches the genomic data for further studies on the evolution, host adaptation, and novel control strategies of thrips. © 2025 Society of Chemical Industry.

Chromosome-level genome assembly of the Adonis ladybird Hippodamia variegata

The Adonis ladybird (Hippodamia variegata), an important predator in agricultural ecosystems, plays a crucial role in biological control and is a significant model for evolutionary and genomic studies within Coccinellidae. Despite its ecological importance, the lack of a reference genome for H. variegata has limited in-depth investigations into its biology and potential as a biocontrol agent. Here, we present a high-quality, chromosome-level genome assembly of H. variegata. The final assembly spans 493.01 Mb, with a scaffold N50 of 28.19 Mb and a GC content of 36.41%. Hi-C sequencing data enabled the anchoring of 343.20 Mb to 10 chromosomes. Repetitive elements accounted for 258.56 Mb (52.44%) of the genome, with long interspersed nuclear elements (LINEs) being the most prevalent. We identified 37,348 protein-coding genes, of which 78.55% were functionally annotated in public protein databases. This high-quality genome assembly will serve as a valuable resource for furthering our understanding of Adonis ladybird's evolutionary biology, enhancing its utility in pest management, and supporting future research on ladybird genomics.

Single specimen genome assembly of Culicoides stellifer shows evidence of a non-retroviral endogenous viral element

BACKGROUND

Advancing our knowledge of vector species genomes is a key step in our battle against the spread of diseases. Biting midges of the genus Culicoides are vectors of arboviruses that significantly affect livestock worldwide. Culicoides stellifer is a suspected vector with a wide range distribution in North America, for which cryptic diversity has been described.

RESULTS

With just one specimen of C. stellifer, we assembled and annotated the nuclear and mitochondrial genome using the ultra-low input DNA PacBio protocol. The genome assembly is 119 Mb in length with a contig N50 value of 479.3 kb, contains 11% repeat sequences and 18,895 annotated protein-coding genes. To further elucidate the role of this species as a vector, we provide genomic evidence of a non-retroviral endogenous viral element integrated into the genome that corresponds to rhabdovirus nucleocapsid proteins, the same family as the vesicular stomatitis virus.

CONCLUSIONS

This genomic information will pave the way for future investigations into this species's putative vector role. We also demonstrate the practicability of completing genomic studies in small dipterans using single specimens preserved in ethanol as well as introduce a workflow for data analysis that considers the challenges of insect genome assembly.

High-quality reference genome of predatory mite Neoseiulus californicus McGregor (Acari: Phytoseiidae) provides insights into its biological traits and potential RNAi off-target effects

BACKGROUND

Neoseiulus californicus is a predatory mite that can control various spider mites and other small arthropods. Despite its acknowledged effectiveness in the natural enemy market, a crucial knowledge gap exists in understanding the genomic features related to its predatory traits and adaptation. With the increasing emphasis on modern pest management strategies and dynamic environmental changes in plant production trends, constructing a reliable genomic resource for N. californicus becomes imperative.

RESULTS

In this study, we provided a high-quality genome assembly and annotation of N. californicus, with a size of 188.43 Mb and 12 946 predicted genes. We identified genomic features and traits related to its detoxification, stress response, sensory system, mobility, secretory toxins, digestive enzymes and horizontally transferred genes (HGTs) by comparative genomics. We also predicted the potential off-target effects of double-stranded RNA (dsRNA) that could be used to control spider mites, based on public data and the N. californicus genome.

CONCLUSION

Our study provided a valuable genomic resource for a commercialized predatory mite, offering useful insights for the design of integrated pest management (IPM) strategies in the new era. Further studies are needed to explore the functional roles of key gene families in predatory mites, as well as their interactions with their prey and IPM strategies. © 2024 Society of Chemical Industry.

Properties of "Stable" Mosquito Cytochrome P450 Enzymes

The use of insecticides is widespread in the control of debilitating mosquito-borne diseases. P450 enzymes (CYPs) play essential roles in mosquito physiological function but also in the enzymatic detoxification of xenobiotics. Broadly speaking, CYPs can be classified as "stable", meaning those that have no or very few paralogs, and "labile", constituting gene families with many paralogous members. The evolutionary dichotomy between "stable" and "labile" P450 genes is fuzzy and there is not a clear phylogenetic demarcation between P450s involved in detoxification and P450s involved in essential metabolic processes. In this study, bioinformatic methods were used to explore differences in the sequences of "stable" and "labile" P450s that may facilitate their functional classification. Genomic and sequence data of Anopheles gambiae (Agam), Aedes aegypti (Aaeg), and Culex quinquefasciatus (Cqui) CYPs were obtained from public databases. The results of this study show that "stable" CYPs are encoded by longer genes, have longer introns and more exons, and contain a higher proportion of hydrophobic amino acids than "labile" CYPs. Compared to "labile" CYPs, a significantly higher proportion of "stable" CYPs are associated with biosynthetic and developmental processes.

CRISPR/Cas9-mediated knockout of Tektin 4-like gene (TEKT4L) causes male sterility of Cydia pomonella

The sterile insect technique (SIT) is a well-established and environmentally benign method for population control. Identifying genes that regulate insect fertility while preserving growth and development is crucial for implementing a novel SIT-based pest management approach utilizing CRISPR/Cas9 to target these genes for genetic manipulation. Tektin (TEKT), an essential alpha-helical protein pivotal in sperm formation due to its role in cilia and flagella assembly, has garnered attention. In this study, we identified 7 TEKT genes in the testis of Cydia pomonella, a globally invasive fruit pest. Notably, Tektin4-like (TEKT4L) displayed the highest expression level in male adult especially the testes, suggesting its significance in reproductive processes. By utilizing CRISPR/Cas9 technology to knockout TEKT4L, male sterility was induced, showcasing dominant inherited. When wild-type (WT) females mated with TEKT4L-/- males, eggs laying proceeded normally, but the hatching rate was dramatically reduced, with only 15.49% progressing to the eyespot stage and 68.86% failing to develop normally. The reproductive fitness of TEKT4L-/- males was robust enough to facilitate the transmission of genetic modifications efficiently within the C.pomonella population, yielding a small number of viable offspring. Subsequent cage trials demonstrated the effectiveness of this population in suppressing laboratory populations of C.pomonella, achieving notable results with a relatively low release ratio (TEKT4L-/-♂: WT♂: WT♀ = 5:1:5). Consequently, the targeted disruption of the TEKT4L gene holds promise as a fundamental element in a novel pest control strategy against C. pomonella.

The Genomics Revolution Drives a New Era in Entomology

Thanks to the fast development of sequencing techniques and bioinformatics tools, sequencing the genome of an insect species for specific research purposes has become an increasingly popular practice. Insect genomes not only provide sets of gene sequences but also represent a change in focus from reductionism to systemic biology in the field of entomology. Using insect genomes, researchers are able to identify and study the functions of all members of a gene family, pathway, or gene network associated with a trait of interest. Comparative genomics studies provide new insights into insect evolution, addressing long-lasting controversies in taxonomy. It is also now feasible to uncover the genetic basis of important traits by identifying variants using genome resequencing data of individual insects, followed by genome-wide association analysis. Here, we review the current progress in insect genome sequencing projects and the application of insect genomes in uncovering the phylogenetic relationships between insects and unraveling the mechanisms of important life-history traits. We also summarize the challenges in genome data sharing and possible solutions. Finally, we provide guidance for fully and deeply mining insect genome data.

Template-based modeling of insect odorant receptors outperforms AlphaFold3 for ligand binding predictions

Insects rely on odorant receptors (ORs) to detect and respond to volatile environmental cues, so the ORs are attracting increasing interest as potential targets for pest control. However, experimental analysis of their structures and functions faces significant challenges. Computational methods such as template-based modeling (TBM) and AlphaFold3 (AF3) could facilitate the structural characterisation of ORs. This study first showed that both models accurately predicted the structural fold of MhOR5, a jumping bristletail OR with known experimental 3D structures, although accuracy was higher in the extracellular region of the protein and binding mode of their cognate ligands with TBM. The two approaches were then compared for their ability to predict the empirical binding evidence available for OR-odorant complexes in two economically important fruit fly species, Bactrocera dorsalis and B. minax. Post-simulation analyses including binding affinities, complex and ligand stability and receptor-ligand interactions (RLIs) revealed that TBM performed better than AF3 in discriminating between binder and non-binder complexes. TBM's superior performance is attributed to hydrophobicity-based helix-wise multiple sequence alignment (MSA) between available insect OR templates and the ORs for which the binding data were generated. This MSA identified conserved residues and motifs which could be used as anchor points for refining the alignments.

Identifying the sex chromosome and sex determination genes in the cattle tick, Rhipicephalus (Boophilus) microplus

Rhipicephalus (Boophilus) microplus is globally one of the most economically important ectoparasites of cattle costing the agriculture industry billions of dollars annually. Resistance to chemical control measures has prompted the development of novel methods of control. Recent advancements in genetic control measures for human and other animal vectors have utilized sex determination research to manipulate sex ratios, which have shown promising results in mosquitoes namely Aedes aegypti and Anopheles stephensi. Here, we use R. (B.) microplus as a model to provide foundational research to allow similar avenues of investigation in ticks using R. (B.) microplus as a model. Karyotypes for R. (B.) microplus show an XX:XO sex determining system with the largest chromosome being the sex chromosome. Using flow cytometric methods, the size of the sex chromosome was estimated at 526.91 Mb. All measures to identify the sex chromosome within the cattle tick genome assembly associated sex chromosomal characteristics to two chromosomes. This discrepancy between the assembly and karyotypes of the tick led to generating a new genome assembly with a single adult male specimen. The two chromosomes in question aligned with a single scaffold within the new genome that had a length of 513.29 Mb and was the first time the sex chromosome was identified in an Ixodid genome assembly.

Transcriptomic and Gene Expression Analysis of Chemosensory Genes from White Grubs of Hylamorpha elegans (Coleoptera: Scarabaeidae), a Subterranean Pest in South America

Olfaction and gustation processes play key roles in the life cycle of insects, such as finding and accepting food sources, oviposition sites, and mates, among other fundamental aspects of insect development. In this context, chemosensory genes found in sensory organs (e.g., antennae and maxillary palps) are crucial for understanding insect behaviour, particularly the phytophagous behaviour of insect pests that attack economically important crops. An example is the scarab beetle Hylamorpha elegans, which feeds on the roots of several crops important for livestock in its larval stage. In this study, chemosensory gene candidates of H. elegans white grubs identified through the head transcriptome and phylogenetic and tissue-biased gene expression (antennae, head without antennae, and legs) have been reported. Overall, 47 chemosensory genes were identified (2 ORs, 1 GR, 11 IRs, 9 CSPs, and 24 OBPs). Gene expression analysis revealed the predominant presence of IRs in the legs, whereas ORs and the GR were present in the heads and/or antennae. Particularly, HeleOBP9 and HeleCSP2 were significantly expressed in the head but not in the antennae or legs; these and other genes are discussed as potential targets in the context of H. elegans management.

SENSE-PPI reconstructs interactomes within, across, and between species at the genome scale

Ab initio computational reconstructions of protein-protein interaction (PPI) networks will provide invaluable insights into cellular systems, enabling the discovery of novel molecular interactions and elucidating biological mechanisms within and between organisms. Leveraging the latest generation protein language models and recurrent neural networks, we present SENSE-PPI, a sequence-based deep learning model that efficiently reconstructs ab initio PPIs, distinguishing partners among tens of thousands of proteins and identifying specific interactions within functionally similar proteins. SENSE-PPI demonstrates high accuracy, limited training requirements, and versatility in cross-species predictions, even with non-model organisms and human-virus interactions. Its performance decreases for phylogenetically more distant model and non-model organisms, but signal alteration is very slow. In this regard, it demonstrates the important role of parameters in protein language models. SENSE-PPI is very fast and can test 10,000 proteins against themselves in a matter of hours, enabling the reconstruction of genome-wide proteomes.

Analysis of maternal effect genes from maternal mRNA in eggs of Sogatella furcifera

To understand how many kinds of mRNAs in female adults can be transferred into the eggs and the molecular basis of embryonic axis specification in Sogatella furcifera, we performed de novo transcriptome sequencing of six cDNA libraries of female adults and unfertilized eggs. Total 60,306 unigenes were obtained, with an average length of 1114.51 bp and N50 length of 2112 bp. Total 2900 differentially expressed genes with 496 upregulated and 2404 downregulated were found in unfertilized egg compared to female adult. Most of mRNAs in female adult could be passed into its eggs. Total 65 maternal effect genes were identified, including many homologous genes involved in embryonic axis specialization of D. melanogaster. This study provide transcriptome resources to elucidate the functions of maternal effect genes and the molecular basis of embryonic axis specification in S. furcifera in the future.

Genome-Wide Identification and Expression Pattern of Sugar Transporter Genes in the Brown Planthopper, Nilaparvata lugens (Stål)

Sugar transporters play important roles in controlling carbohydrate transport and are responsible for mediating the movement of sugars into cells in numerous organisms. In insects, sugar transporters not only play a role in sugar transport but may also act as receptors for virus entry and the accumulation of plant defense compounds. The brown planthopper, Nilaparvata lugens, inflicts damage on rice plants by feeding on their phloem sap, which is rich in sugars. In the present study, we identified 34 sugar transporters in N. lugens, which were classified into three subfamilies based on phylogenetic analysis. The motif numbers varied from seven to eleven, and motifs 2, 3, and 4 were identified in the functional domains of all 34 NlST proteins. Chromosome 1 was found to possess the highest number of NlST genes, harboring 15. The gut, salivary glands, fat body, and ovary were the different tissues enriched with NlST gene expression. The expression levels of NlST2, 3, 4, 7, 20, 27, 28, and 31 were higher in the gut than in the other tissues. When expressed in a Saccharomyces cerevisiae hexose transporter deletion mutant (strain EBY.VW4000), only ApST4 (previously characterized) and NlST4, 28, and 31 were found to transport glucose and fructose, resulting in functional rescue of the yeast mutant. These results provide valuable data for further studies on sugar transporters in N. lugens and lay a foundation for finding potential targets to control N. lugens.

Chromosome-Level Genome Assembly of Apoderus dimidiatus Voss (Coleoptera: Attelabidae): Insights into Evolution and Behavior

Attelabidae insects have attracted much attention due to their unique leaf rolling behavior before oviposition. However, the lack of genomic data makes it difficult to understand the molecular mechanism behind their behavior and their evolutionary relationship with other species. To address this gap, we utilized Illumina and Nanopore sequencing platforms along with Hi-C technology to establish a highly accurate whole genome of A. dimidiatus at the chromosome level. The resulting genome size was determined to be 619.26 Mb, with a contig N50 of 50.89 Mb and GC content of 33.89%. Moreover, a total of 12,572 genes were identified, with 82.59% being functionally annotated, and 64.78% designated as repeat sequences. Our subsequent phylogenetic tree analysis revealed that Attelabidae's divergence from Curculionidae occurred approximately 161.52 million years ago. Furthermore, the genome of A. dimidiatus contained 334 expanded gene families and 1718 contracted gene families. In addition, using Phylogenetic Analysis by Maximum Likelihood (PAML), we identified 106 rapidly evolved genes exhibiting significant signals and 540 positively selected genes. Our research endeavors to serve as an invaluable genomic data resource for the study of Attelabidae, offering fresh perspectives for the exploration of its leaf rolling behavior.

Protamines and the sperm nuclear basic proteins Pandora's Box of insects

Insects are the largest group of animals when it comes to the number and diversity of species. Yet, with the exception of Drosophila, no information is currently available on the primary structure of their sperm nuclear basic proteins (SNBPs). This paper represents the first attempt in this regard and provides information about six species of Neoptera: Poecillimon thessalicus, Graptosaltria nigrofuscata, Apis mellifera, Nasonia vitripennis, Parachauliodes continentalis, and Tribolium castaneum. The SNBPs of these species were characterized by acetic acid urea gel electrophoresis (AU-PAGE) and high-performance liquid chromatography fractionated. Protein sequencing was obtained using a combination of mass spectrometry sequencing, Edman N-terminal degradation sequencing and genome mining. While the SNBPs of several of these species exhibit a canonical arginine-rich protamine nature, a few of them exhibit a protamine-like composition. They appear to be the products of extensive cleavage processing from a precursor protein which are sometimes further processed by other post-translational modifications that are likely involved in the chromatin transitions observed during spermiogenesis in these organisms.

Transcriptome-informed identification and characterization of Planococcus citri cis- and trans-isoprenyl diphosphate synthase genes

Insect physiology and reproduction depend on several terpenoid compounds, whose biosynthesis is mainly unknown. One enigmatic group of insect monoterpenoids are mealybug sex pheromones, presumably resulting from the irregular coupling activity of unidentified isoprenyl diphosphate synthases (IDSs). Here, we performed a comprehensive search for IDS coding sequences of the pest mealybug Planococcus citri. We queried the available genomic and newly generated short- and long-read P. citri transcriptomic data and identified 18 putative IDS genes, whose phylogenetic analysis indicates several gene family expansion events. In vitro testing confirmed regular short-chain coupling activity with five gene products. With the candidate with highest IDS activity, we also detected low amounts of irregular coupling products, and determined amino acid residues important for chain-length preference and irregular coupling activity. This work therefore provides an important foundation for deciphering terpenoid biosynthesis in mealybugs, including the sex pheromone biosynthesis in P. citri.

Activation of immune pathways in common bed bugs, Cimex lectularius, in response to bacterial immune challenges - a transcriptomics analysis

The common bed bug, Cimex lectularius, is an urban pest of global health significance, severely affecting the physical and mental health of humans. In contrast to most other blood-feeding arthropods, bed bugs are not major vectors of pathogens, but the underlying mechanisms for this phenomenon are largely unexplored. Here, we present the first transcriptomics study of bed bugs in response to immune challenges. To study transcriptional variations in bed bugs following ingestion of bacteria, we extracted and processed mRNA from body tissues of adult male bed bugs after ingestion of sterile blood or blood containing the Gram-positive (Gr+) bacterium Bacillus subtilis or the Gram-negative (Gr-) bacterium Escherichia coli. We analyzed mRNA from the bed bugs' midgut (the primary tissue involved in blood ingestion) and from the rest of their bodies (RoB; body minus head and midgut tissues). We show that the midgut exhibits a stronger immune response to ingestion of bacteria than the RoB, as indicated by the expression of genes encoding antimicrobial peptides (AMPs). Both the Toll and Imd signaling pathways, associated with immune responses, were highly activated by the ingestion of bacteria. Bacterial infection in bed bugs further provides evidence for metabolic reconfiguration and resource allocation in the bed bugs' midgut and RoB to promote production of AMPs. Our data suggest that infection with particular pathogens in bed bugs may be associated with altered metabolic pathways within the midgut and RoB that favors immune responses. We further show that multiple established cellular immune responses are preserved and are activated by the presence of specific pathogens. Our study provides a greater understanding of nuances in the immune responses of bed bugs towards pathogens that ultimately might contribute to novel bed bug control tactics.

Antennal Transcriptome Evaluation and Analysis for Odorant-Binding Proteins, Chemosensory Proteins, and Suitable Reference Genes in the Leaf Beetle Pest Diorhabda rybakowi Weise (Coleoptera: Chrysomelidae)

Diorhabda rybakowi Weise is one of the dominant pests feeding on Nitraria spp., a pioneer plant used for windbreaking and sand fixation purposes, and poses a threat to local livestock and ecosystems. To clarify the key olfactory genes of D. rybakowi and provide a theoretical basis for attractant and repellent development, the optimal reference genes under two different conditions (tissue and sex) were identified, and the bioinformatics and characterization of the tissue expression profiles of two categories of soluble olfactory proteins (OBPs and CSPs) were investigated. The results showed that the best reference genes were RPL13a and RPS18 for comparison among tissues, and RPL19 and RPS18 for comparison between sexes. Strong expressions of DrybOBP3, DrybOBP6, DrybOBP7, DrybOBP10, DrybOBP11, DrybCSP2, and DrybCSP5 were found in antennae, the most important olfactory organ for D. rybakowi. These findings not only provide a basis for further in-depth research on the olfactory molecular mechanisms of host-specialized pests but also provide a theoretical basis for the future development of new chemical attractants or repellents using volatiles to control D. rybakowi.

GC Content Across Insect Genomes: Phylogenetic Patterns, Causes and Consequences

The proportions of A:T and G:C nucleotide pairs are often unequal and can vary greatly between animal species and along chromosomes. The causes and consequences of this variation are incompletely understood. The recent release of high-quality genome sequences from the Darwin Tree of Life and other large-scale genome projects provides an opportunity for GC heterogeneity to be compared across a large number of insect species. Here we analyse GC content along chromosomes, and within protein-coding genes and codons, of 150 insect species from four holometabolous orders: Coleoptera, Diptera, Hymenoptera, and Lepidoptera. We find that protein-coding sequences have higher GC content than the genome average, and that Lepidoptera generally have higher GC content than the other three insect orders examined. GC content is higher in small chromosomes in most Lepidoptera species, but this pattern is less consistent in other orders. GC content also increases towards subtelomeric regions within protein-coding genes in Diptera, Coleoptera and Lepidoptera. Two species of Diptera, Bombylius major and B. discolor, have very atypical genomes with ubiquitous increase in AT content, especially at third codon positions. Despite dramatic AT-biased codon usage, we find no evidence that this has driven divergent protein evolution. We argue that the GC landscape of Lepidoptera, Diptera and Coleoptera genomes is influenced by GC-biased gene conversion, strongest in Lepidoptera, with some outlier taxa affected drastically by counteracting processes.

Unsupervised AI reveals insect species-specific genome signatures

Insects are a highly diverse phylogeny and possess a wide variety of traits, including the presence or absence of wings and metamorphosis. These diverse traits are of great interest for studying genome evolution, and numerous comparative genomic studies have examined a wide phylogenetic range of insects. Here, we analyzed 22 insects belonging to a wide phylogenetic range (Endopterygota, Paraneoptera, Polyneoptera, Palaeoptera, and other insects) by using a batch-learning self-organizing map (BLSOM) for oligonucleotide compositions in their genomic fragments (100-kb or 1-Mb sequences), which is an unsupervised machine learning algorithm that can extract species-specific characteristics of the oligonucleotide compositions (genome signatures). The genome signature is of particular interest in terms of the mechanisms and biological significance that have caused the species-specific difference, and can be used as a powerful search needle to explore the various roles of genome sequences other than protein coding, and can be used to unveil mysteries hidden in the genome sequence. Since BLSOM is an unsupervised clustering method, the clustering of sequences was performed based on the oligonucleotide composition alone, without providing information about the species from which each fragment sequence was derived. Therefore, not only the interspecies separation, but also the intraspecies separation can be achieved. Here, we have revealed the specific genomic regions with oligonucleotide compositions distinct from the usual sequences of each insect genome, e.g., Mb-level structures found for a grasshopper Schistocerca americana. One aim of this study was to compare the genome characteristics of insects with those of vertebrates, especially humans, which are phylogenetically distant from insects. Recently, humans seem to be the "model organism" for which a large amount of information has been accumulated using a variety of cutting-edge and high-throughput technologies. Therefore, it is reasonable to use the abundant information from humans to study insect lineages. The specific regions of Mb length with distinct oligonucleotide compositions have also been previously observed in the human genome. These regions were enriched by transcription factor binding motifs (TFBSs) and hypothesized to be involved in the three-dimensional arrangement of chromosomal DNA in interphase nuclei. The present study characterized the species-specific oligonucleotide compositions (i.e., genome signatures) in insect genomes and identified specific genomic regions with distinct oligonucleotide compositions.

Direct sequencing of insect symbionts via nanopore adaptive sampling

Insect symbionts can alter their host phenotype and their effects can range from beneficial to pathogenic. Moreover, many insects exhibit co-infections, making their study more challenging. Less than 1% of insect species have high-quality referenced genomes available and fewer still also have their symbionts sequenced. Two methods are commonly used to sequence symbionts: whole-genome sequencing to concomitantly capture the host and bacterial genomes, or isolation of the symbiont's genome before sequencing. These methods are limited when dealing with rare or poorly characterized symbionts. Long-read technology is an important tool to generate high-quality genomes as they can overcome high levels of heterozygosity, repeat content, and transposable elements that confound short-read methods. Oxford Nanopore (ONT) adaptive sampling allows a sequencing instrument to select or reject sequences in real time. We describe a method based on ONT adaptive sampling (subtractive) approach that readily permitted the sequencing of the complete genomes of mitochondria, Buchnera and its plasmids (pLeu, pTrp), and Wolbachia genomes in two aphid species, Aphis glycines and Pentalonia nigronervosa. Adaptive sampling is able to retrieve organelles such as mitochondria and symbionts that have high representation in their hosts such as Buchnera and Wolbachia, but is less successful at retrieving symbionts in low concentrations.

Analysis of Genes Associated with Feeding Preference and Detoxification in Various Developmental Stages of Aglais urticae

Herbivorous insects and host plants have developed a close and complex relationship over a long period of co-evolution. Some plants provide nutrients for insects, but plants' secondary metabolites also influence their growth and development. Urtica cannabina roots and leaves are poisonous, yet Aglais urticae larvae feed on them, so we aimed to clarify the mechanism enabling this interaction. At present, studies on the detoxification mechanism of the A. urticae are rare. In our study, first, we used the A. urticae larval odor selection behavior bioassay and choice feeding preference assay to analyze the feeding preferences of A. urticae on its host plant, U. cannabina. Next, we used transcriptome sequencing to obtain the unigenes annotated and classified by various databases, such as KEGG and GO. In this study, we found that U. cannabina could attract A. urticae larvae to feed via scent, and the feeding preference assay confirmed that larvae preferred U. cannabina leaves over three other plants: Cirsium japonicum, Cannabis sativa, and Arctium lappa. The activity of detoxifying enzymes GST and CarE changed in larvae that had consumed U. cannabina. Furthermore, through transcriptomic sequencing analysis, 77,624 unigenes were assembled from raw reads. The numbers of differentially expressed genes were calculated using pairwise comparisons of all life stages; the expression of detoxification enzyme genes was substantially higher in larvae than in the pupal and adult stages. Finally, we identified and summarized 34 genes associated with detoxification enzymes, such as UDP-glucose 4-epimerase gene, 5 Glutathione S-transferase genes, 4 Carboxylesterase genes, 4 Cytochrome P450 genes, 10 ATP-binding cassette genes, 4 Superoxide dismutase, and Peroxidase. Moreover, we identified 28 genes associated with the development of A. urticae. The qRT-PCR results were nearly consistent with the transcriptomic data, showing an increased expression level of four genes in larvae. Taken together, this study examines the correlation between A. urticae and host plants U. cannabina, uncovering a pronounced preference for A. urticae larvae toward host plants. Consistent with RNA-seq, we investigated the mechanism of A. urticae's interaction with host plants and identified detoxification-related genes. The present study provides theoretical support for studying insect adaptation mechanisms and biological control.

A chromosome-level genome assembly of the Rhus gall aphid Schlechtendalia chinensis provides insight into the endogenization of Parvovirus-like DNA sequences

The Rhus gall aphid, Schlechtendalia chinensis, feeds on its primary host plant Rhus chinensis to induce galls, which have economic importance in medicines and the food industry. Rhus gall aphids have a unique life cycle and are economically beneficial but there is huge gap in genomic information about this group of aphids. Schlechtendalia chinensis induces rich-tannin galls on its host plant and is emerging as a model organism for both commercial applications and applied research in the context of gall production by insects. Here, we generated a high-quality chromosome-level assembly for the S. chinensis genome, enabling the comparison between S. chinensis and non-galling aphids. The final genome assembly is 344.59 Mb with 91.71% of the assembled sequences anchored into 13 chromosomes. We predicted 15,013 genes, of which 14,582 (97.13%) coding genes were annotated, and 99% of the predicted genes were anchored to the 13 chromosomes. This assembly reveals the endogenization of parvovirus-related DNA sequences (PRDs) in the S. chinensis genome, which could play a role in environmental adaptations. We demonstrated the characterization and classification of cytochrome P450s in the genome assembly, which are functionally crucial for sap-feeding insects and have roles in detoxification and insecticide resistance. This genome assembly also revealed the whole genome duplication events in S. chinensis, which can be considered in comparative evolutionary analysis. Our work represents a reference genome for gall-forming aphids that could be used for comparative genomic studies between galling and non-galling aphids and provides the first insight into the endogenization of PRDs in the genome of galling aphids. It also provides novel genetic information for future research on gall-formation and insect-plant interactions.

The Development and Expansion of in vivo Germline Editing Technologies in Arthropods: Receptor-Mediated Ovary Transduction of Cargo (ReMOT Control) and Beyond

In the past 20 years, sequencing technologies have led to easy access to genomic data from nonmodel organisms in all biological realms. Insect genetic manipulation, however, continues to be a challenge due to various factors, including technical and cost-related issues. Traditional techniques such as microinjection of gene-editing vectors into early stage embryos have been used for arthropod transgenesis and the discovery of Clustered regularly interspaced short palindromic repeats and CRISPR-associated protein (CRISPR-Cas) technologies allowed for targeted mutagenesis and the creation of knockouts or knock-ins in arthropods. Receptor-Mediated Ovary Transduction of Cargo (ReMOT Control) acts as an alternative to embryonic microinjections, which require expensive equipment and extensive hands-on training. ReMOT Control's main advantage is its ease of use coupled with the ability to hypothetically target any vitellogenic species, as injections are administered to the egg-laying adult rather than embryos. After its initial application in the mosquito Aedes aegypti, ReMOT Control has successfully produced mutants not only for mosquitoes but for multiple arthropod species from diverse orders, such as ticks, mites, wasps, beetles, and true bugs, and is being extended to crustaceans, demonstrating the versatility of the technique. In this review, we discuss the current state of ReMOT Control from its proof-of-concept to the advances and challenges in the application across species after 5 years since its development, including novel extensions of the technique such as direct parental (DIPA)-CRISPR.

Advances in genome sequencing reveal changes in gene content that contribute to arthropod macroevolution

Current sequencing technology allows for the relatively affordable generation of highly contiguous genomes. Technological advances have made it possible for researchers to investigate the consequences of diverse sorts of genomic variants, such as gene gain and loss. With the extraordinary number of high-quality genomes now available, we take stock of how these genomic variants impact phenotypic evolution. We take care to point out that the identification of genomic variants of interest is only the first step in understanding their impact. Painstaking lab or fieldwork is still required to establish causal relationships between genomic variants and phenotypic evolution. We focus mostly on arthropod research, as this phylum has an impressive degree of phenotypic diversity and is also the subject of much evolutionary genetics research. This article is intended to both highlight recent advances in the field and also to be a primer for learning about evolutionary genetics and genomics.

The loci of insect phenotypic evolution

Insects are important elements of terrestrial ecosystems because they pollinate plants, destroy crops, transmit diseases to livestock and humans, and are important components of food chains. Here, I used Gephebase, a manually curated database of genetic variants associated with natural and domesticated trait variation, to explore current knowledge about the genes and the mutations known to contribute to natural phenotypic variation in insects. Analysis of over 600 mutations reveals that data are concentrated toward certain species and traits and that experimental approaches have changed over time. The distribution of coding and cis-regulatory changes varies with traits, experimental approaches, and identified gene loci. Recent studies highlight the important role of standing variation, repeated mutations in hotspot genes, recombination, inversions, and introgression.

The highly improved genome of Ixodes scapularis with X and Y pseudochromosomes

Ixodes scapularis, the black-legged tick, is the principal vector of the Lyme disease spirochete, Borrelia burgdorferi, and is responsible for most of the ∼470,000 estimated Lyme disease cases annually in the USA. Ixodes scapularis can transmit six additional pathogens of human health significance. Because of its medical importance, I. scapularis was the first tick genome to be sequenced and annotated. However, the first assembly, I. scapularis Wikel (IscaW), was highly fragmented because of the technical challenges posed by the long, repetitive genome sequences characteristic of arthropod genomes and the lack of long-read sequencing techniques. Although I. scapularis has emerged as a model for tick research because of the availability of new tools such as embryo injection and CRISPR-Cas9-mediated gene editing yet the lack of chromosome-scale scaffolds has slowed progress in tick biology and the development of tools for their control. Here we combine diverse technologies to produce the I. scapularis Gulia-Nuss (IscGN) genome assembly and gene set. We used DNA from eggs and male and female adult ticks and took advantage of Hi-C, PacBio HiFi sequencing, and Illumina short-read sequencing technologies to produce a chromosome-level assembly. In this work, we present the predicted pseudochromosomes consisting of 13 autosomes and the sex pseudochromosomes: X and Y, and a markedly improved genome annotation compared with the existing assemblies and annotations.

In silico SNP prediction of selected protein orthologues in insect models for Alzheimer's, Parkinson's, and Huntington's diseases

Alzheimer's, Parkinson's, and Huntington's are the most common neurodegenerative diseases that are incurable and affect the elderly population. Discovery of effective treatments for these diseases is often difficult, expensive, and serendipitous. Previous comparative studies on different model organisms have revealed that most animals share similar cellular and molecular characteristics. The meta-SNP tool includes four different integrated tools (SIFT, PANTHER, SNAP, and PhD-SNP) was used to identify non synonymous single nucleotide polymorphism (nsSNPs). Prediction of nsSNPs was conducted on three representative proteins for Alzheimer's, Parkinson's, and Huntington's diseases; APPl in Drosophila melanogaster, LRRK1 in Aedes aegypti, and VCPl in Tribolium castaneum. With the possibility of using insect models to investigate neurodegenerative diseases. We conclude from the protein comparative analysis between different insect models and nsSNP analyses that D. melanogaster is the best model for Alzheimer's representing five nsSNPs of the 21 suggested mutations in the APPl protein. Aedes aegypti is the best model for Parkinson's representing three nsSNPs in the LRRK1 protein. Tribolium castaneum is the best model for Huntington's disease representing 13 SNPs of 37 suggested mutations in the VCPl protein. This study aimed to improve human neural health by identifying the best insect to model Alzheimer's, Parkinson's, and Huntington's.

Chromosome study of the Hymenoptera (Insecta): from cytogenetics to cytogenomics

A brief overview of the current stage of the chromosome study of the insect order Hymenoptera is given. It is demonstrated that, in addition to routine staining and other traditional techniques of chromosome research, karyotypes of an increasing number of hymenopterans are being studied using molecular methods, e.g., staining with base-specific fluorochromes and fluorescence in situ hybridization (FISH), including microdissection and chromosome painting. Due to the advent of whole genome sequencing and other molecular techniques, together with the "big data" approach to the chromosomal data, the current stage of the chromosome research on Hymenoptera represents a transition from Hymenoptera cytogenetics to cytogenomics.

Chromosome-level genome assembly of Bactrocera correcta provides insights into its adaptation and invasion mechanisms

Bactrocera correcta is an invasive polyphagous pest with significant ecological and economic implications. Understanding its genetic characteristics and the molecular mechanisms that drive its rapid adaptation to new environments requires genomic information. In this study, we successfully assembled the chromosome-level genome of B. correcta using PacBio long-read sequencing, Illumina sequencing, and chromatin conformation capture (Hi-C) methods. The final genome assembly spans a total length of 702.65 Mb. We managed to anchor approximately 86.88% of the assembled contigs into 6 linkage groups, ranging from 17.97 Mb to 166.49 Mb. Additionally, our analysis predicted a total of 21,015 genes, with repetitive sequences accounting for 58.22% of the genome. We further identified retroelements and DNA transposons as the major contributors to the larger size of the B. correcta genome, constituting 36.06% and 30.92% of the repetitive sequences, respectively. Our divergence time estimation placed B. correcta's split from other Bactrocera species at around 5.99-16.71 million years ago. Through gene family analyses, we discovered significant expansions in sensing-related gene families (IR, GR), heat shock proteins (HSP60), and resistance-related gene families (ABC) in B. correcta compared to its closest relatives. Transcriptomic analysis revealed substantial upregulation of HSP genes, especially those from the HSP20 subfamily, in response to high temperatures. The availability of this reference genome serves as a foundation for the identification of precise target genes in B. correcta, facilitating molecular prevention and control strategies.

The Relationship between Cognition and Brain Size or Neuron Number

The comparative approach is a powerful way to explore the relationship between brain structure and cognitive function. Thus far, the field has been dominated by the assumption that a bigger brain somehow means better cognition. Correlations between differences in brain size or neuron number between species and differences in specific cognitive abilities exist, but these correlations are very noisy. Extreme differences exist between clades in the relationship between either brain size or neuron number and specific cognitive abilities. This means that correlations become weaker, not stronger, as the taxonomic diversity of sampled groups increases. Cognition is the outcome of neural networks. Here we propose that considering plausible neural network models will advance our understanding of the complex relationships between neuron number and different aspects of cognition. Computational modelling of networks suggests that adding pathways, or layers, or changing patterns of connectivity in a network can all have different specific consequences for cognition. Consequently, models of computational architecture can help us hypothesise how and why differences in neuron number might be related to differences in cognition. As methods in connectomics continue to improve and more structural information on animal brains becomes available, we are learning more about natural network structures in brains, and we can develop more biologically plausible models of cognitive architecture. Natural animal diversity then becomes a powerful resource to both test the assumptions of these models and explore hypotheses for how neural network structure and network size might delimit cognitive function.

The Chromosome-Level Genome Assembly of Bean Blossom Thrips (Megalurothrips usitatus) Reveals an Expansion of Protein Digestion-Related Genes in Adaption to High-Protein Host Plants

Megalurothrips usitatus (Bagnall) is a destructive pest of legumes, such as cowpea. The biology, population dynamics and control strategies of this pest have been well studied. However, the lack of a high-quality reference genome for M. usitatus has hindered the understanding of key biological questions, such as the mechanism of adaptation to feed preferentially on high-protein host plants and the resistance to proteinase inhibitors (PIs). In this study, we generated a high-resolution chromosome-level reference genome assembly (247.82 Mb, 16 chromosomes) of M. usitatus by combining Oxford Nanopore Technologies (ONT) and Hi-C sequencing. The genome assembly showed higher proportions of GC and repeat content compared to other Thripinae species. Genome annotation revealed 18,624 protein-coding genes, including 4613 paralogs that were preferentially located in TE-rich regions. GO and KEGG enrichment analyses of the paralogs revealed significant enrichment in digestion-related genes. Genome-wide identification uncovered 506 putative digestion-related enzymes; of those, proteases, especially their subgroup serine proteases (SPs), are significantly enriched in paralogs. We hypothesized that the diversity and expansion of the digestion-related genes, especially SPs, could be driven by mobile elements (TEs), which promote the adaptive evolution of M. usitatus to high-protein host plants with high serine protease inhibitors (SPIs). The current study provides a valuable genomic resource for understanding the genetic variation among different pest species adapting to different plant hosts.

Functional roles of two novel P450 genes in the adaptability of Conogethes punctiferalis to three commonly used pesticides

Introduction: Insect cytochrome P450 (CYP450) genes play important roles in the detoxification and metabolism of xenobiotics, such as plant allelochemicals, mycotoxins and pesticides. The polyphagous Conogethes punctiferalis is a serious economic pest of fruit trees and agricultural crops, and it shows high adaptability to different living environments. Methods: The two novel P450 genes CYP6CV1 and CYP6AB51 were identified and characterized. Quantitative real-time PCR (qRT-PCR) technology was used to study the expression patterns of the two target genes in different larval developmental stages and tissues of C. punctiferalis. Furthermore, RNA interference (RNAi) technology was used to study the potential functions of the two P450 genes by treating RNAi-silenced larvae with three commonly used pesticides. Results: The CYP6CV1 and CYP6AB51 genes were expressed throughout various C. punctiferalis larval stages and in different tissues. Their expression levels increased along with larval development, and expression levels of the two target genes in the midgut were significantly higher than in other tissues. The toxicity bioassay results showed that the LC₅₀ values of chlorantraniliprole, emamectin benzoate and lambda-cyhalothrin on C. punctiferalis larvae were 0.2028 μg/g, 0.0683 μg/g and 0.6110 mg/L, respectively. After treating with different concentrations of chlorantraniliprole, emamectin benzoate and lambda-cyhalothrin (LC₁₀, LC₃₀, LC₅₀), independently, the relative expressions of the two genes CYP6CV1 and CYP6AB51 were significantly induced. After the dsRNA injection, the expression profiles of the two CYP genes were reduced 72.91% and 70.94%, respectively, and the mortality rates of the larvae significantly increased when treated with the three insecticides independently at LC₁₀ values. Discussion: In the summary, after interfering with the CYP6CV1 and CYP6AB51 in C. punctiferalis, respectively, the sensitivity of C. punctiferalis to chlorantraniliprole, emamectin benzoate and lambda-cyhalothrin was significantly increased, indicating that the two CYP6 genes were responsible for the adaptability of C. punctiferalis to the three chemical insecticides in C. punctiferalis. The results from this study demonstrated that CYP6CV1 and CYP6AB51 in C. punctiferalis play crucial roles in the detoxification of chlorantraniliprole, emamectin benzoate and lambda-cyhalothrin.

Fully sequencing the cassava full-length cDNA library reveals unannotated transcript structures and alternative splicing events in regions with a high density of single nucleotide variations, insertions-deletions, and heterozygous sequences

The primary transcript structure provides critical insights into protein diversity, transcriptional modification, and functions. Cassava transcript structures are highly diverse because of alternative splicing (AS) events and high heterozygosity. To precisely determine and characterize transcript structures, fully sequencing cloned transcripts is the most reliable method. However, cassava annotations were mainly determined according to fragmentation-based sequencing analyses (e.g., EST and short-read RNA-seq). In this study, we sequenced the cassava full-length cDNA library, which included rare transcripts. We obtained 8,628 non-redundant fully sequenced transcripts and detected 615 unannotated AS events and 421 unannotated loci. The different protein sequences resulting from the unannotated AS events tended to have diverse functional domains, implying that unannotated AS contributes to the truncation of functional domains. The unannotated loci tended to be derived from orphan genes, implying that the loci may be associated with cassava-specific traits. Unexpectedly, individual cassava transcripts were more likely to have multiple AS events than Arabidopsis transcripts, suggestive of the regulated interactions between cassava splicing-related complexes. We also observed that the unannotated loci and/or AS events were commonly in regions with abundant single nucleotide variations, insertions-deletions, and heterozygous sequences. These findings reflect the utility of completely sequenced FLcDNA clones for overcoming cassava-specific annotation-related problems to elucidate transcript structures. Our work provides researchers with transcript structural details that are useful for annotating highly diverse and unique transcripts and alternative splicing events.

Putting hornets on the genomic map

Hornets are the largest of the social wasps, and are important regulators of insect populations in their native ranges. Hornets are also very successful as invasive species, with often devastating economic, ecological and societal effects. Understanding why these wasps are such successful invaders is critical to managing future introductions and minimising impact on native biodiversity. Critical to the management toolkit is a comprehensive genomic resource for these insects. Here we provide the annotated genomes for two hornets, Vespa crabro and Vespa velutina. We compare their genomes with those of other social Hymenoptera, including the northern giant hornet Vespa mandarinia. The three hornet genomes show evidence of selection pressure on genes associated with reproduction, which might facilitate the transition into invasive ranges. Vespa crabro has experienced positive selection on the highest number of genes, including those putatively associated with molecular binding and olfactory systems. Caste-specific brain transcriptomic analysis also revealed 133 differentially expressed genes, some of which are associated with olfactory functions. This report provides a spring-board for advancing our understanding of the evolution and ecology of hornets, and opens up opportunities for using molecular methods in the future management of both native and invasive populations of these over-looked insects.

The genome of the soybean gall midge (Resseliella maxima)

The cecidomyiid fly, soybean gall midge, Resseliella maxima Gagné, is a recently discovered insect that feeds on soybean plants in the Midwestern United States. R. maxima larvae feed on soybean stems that may induce plant death and can cause considerable yield losses, making it an important agricultural pest. From three pools of 50 adults each, we used long-read nanopore sequencing to assemble a R. maxima reference genome. The final genome assembly is 206 Mb with 64.88× coverage, consisting of 1,009 contigs with an N50 size of 714 kb. The assembly is high quality with a Benchmarking Universal Single-Copy Ortholog (BUSCO) score of 87.8%. Genome-wide GC level is 31.60%, and DNA methylation was measured at 1.07%. The R. maxima genome is comprised of 21.73% repetitive DNA, which is in line with other cecidomyiids. Protein prediction annotated 14,798 coding genes with 89.9% protein BUSCO score. Mitogenome analysis indicated that R. maxima assembly is a single circular contig of 15,301 bp and shares highest identity to the mitogenome of the Asian rice gall midge, Orseolia oryzae Wood-Mason. The R. maxima genome has one of the highest completeness levels for a cecidomyiid and will provide a resource for research focused on the biology, genetics, and evolution of cecidomyiids, as well as plant-insect interactions in this important agricultural pest.

Genome and Genetic Engineering of the House Cricket (Acheta domesticus): A Resource for Sustainable Agriculture

Background: The house cricket, Acheta domesticus, is one of the most farmed insects worldwide and the foundation of an emerging industry using insects as a sustainable food source. Edible insects present a promising alternative for protein production amid a plethora of reports on climate change and biodiversity loss largely driven by agriculture. As with other crops, genetic resources are needed to improve crickets for food and other applications. Methods: We present the first high quality annotated genome assembly of A. domesticus from long read data and scaffolded to chromosome level, providing information needed for genetic manipulation. Results: Gene groups related to immunity were annotated and will be useful for improving value to insect farmers. Metagenome scaffolds in the A. domesticus assembly, including Invertebrate Iridescent Virus 6 (IIV6), were submitted as host-associated sequences. We demonstrate both CRISPR/Cas9-mediated knock-in and knock-out of A. domesticus and discuss implications for the food, pharmaceutical, and other industries. RNAi was demonstrated to disrupt the function of the vermilion eye-color gene producing a useful white-eye biomarker phenotype. Conclusions: We are utilizing these data to develop technologies for downstream commercial applications, including more nutritious and disease-resistant crickets, as well as lines producing valuable bioproducts, such as vaccines and antibiotics.

DNA sequencing in the classroom: complete genome sequence of two earwig (Dermaptera; Insecta) species

BACKGROUND

Despite representing the largest fraction of animal life, the number of insect species whose genome has been sequenced is barely in the hundreds. The order Dermaptera (the earwigs) suffers from a lack of genomic information despite its unique position as one of the basally derived insect groups and its importance in agroecosystems. As part of a national educational and outreach program in genomics, a plan was formulated to engage the participation of high school students in a genome sequencing project. Students from twelve schools across Chile were instructed to capture earwig specimens in their geographical area, to identify them and to provide material for genome sequencing to be carried out by themselves in their schools.

RESULTS

The school students collected specimens from two cosmopolitan earwig species: Euborellia annulipes (Fam. Anisolabididae) and Forficula auricularia (Fam. Forficulidae). Genomic DNA was extracted and, with the help of scientific teams that traveled to the schools, was sequenced using nanopore sequencers. The sequence data obtained for both species was assembled and annotated. We obtained genome sizes of 1.18 Gb (F. auricularia) and 0.94 Gb (E. annulipes) with the number of predicted protein coding genes being 31,800 and 40,000, respectively. Our analysis showed that we were able to capture a high percentage (≥ 93%) of conserved proteins indicating genomes that are useful for comparative and functional analysis. We were also able to characterize structural elements such as repetitive sequences and non-coding RNA genes. Finally, functional categories of genes that are overrepresented in each species suggest important differences in the process underlying the formation of germ cells, and modes of reproduction between them, features that are one of the distinguishing biological properties that characterize these two distant families of Dermaptera.

CONCLUSIONS

This work represents an unprecedented instance where the scientific and lay community have come together to collaborate in a genome sequencing project. The versatility and accessibility of nanopore sequencers was key to the success of the initiative. We were able to obtain full genome sequences of two important and widely distributed species of insects which had not been analyzed at this level previously. The data made available by the project should illuminate future studies on the Dermaptera.

The Genome of the Soybean Gall Midge ( Resseliella maxima )

The cecidomyiid fly, soybean gall midge, Resseliella maxima Gagné, is a recently discovered insect that feeds on soybean plants in the Midwest US. Resseliella maxima larvae feed on soybean stems which may induce plant death and can cause considerable yield losses, making it an important agricultural pest. From three pools of 50 adults each, we used long-read nanopore sequencing to assemble a R. maxima reference genome. The final genome assembly is 206 Mb with 64.88X coverage, consisting of 1009 contigs with an N50 size of 714 kb. The assembly is high quality with a BUSCO score of 87.8%. Genome-wide GC level is 31.60% and DNA methylation was measured at 1.07%. The R. maxima genome is comprised of 21.73% repetitive DNA, which is in line with other cecidomyiids. Protein prediction annotated 14,798 coding genes with 89.9% protein BUSCO score. Mitogenome analysis indicated that R. maxima assembly is a single circular contig of 15,301 bp and shares highest identity to the mitogenome of the Asian rice gall midge, Orseolia oryzae (Wood-Mason). The R. maxima genome has one of the highest completeness levels for a cecidomyiid and will provide a resource for research focused on the biology, genetics, and evolution of cecidomyiids, as well as plant-insect interactions in this important agricultural pest.

An Efficient Feature Selection Algorithm for Gene Families Using NMF and ReliefF

Gene families, which are parts of a genome's information storage hierarchy, play a significant role in the development and diversity of multicellular organisms. Several studies have focused on the characteristics of gene families, such as function, homology, or phenotype. However, statistical and correlation analyses on the distribution of gene family members in the genome have yet to be conducted. Here, a novel framework incorporating gene family analysis and genome selection based on NMF-ReliefF is reported. Specifically, the proposed method starts by obtaining gene families from the TreeFam database and determining the number of gene families within the feature matrix. Then, NMF-ReliefF is used to select features from the gene feature matrix, which is a new feature selection algorithm that overcomes the inefficiencies of traditional methods. Finally, a support vector machine is utilized to classify the acquired features. The results show that the framework achieved an accuracy of 89.1% and an AUC of 0.919 on the insect genome test set. We also employed four microarray gene data sets to evaluate the performance of the NMF-ReliefF algorithm. The outcomes show that the proposed method may strike a delicate balance between robustness and discrimination. Additionally, the proposed method's categorization is superior to state-of-the-art feature selection approaches.

Migratory behaviour is positively associated with genetic diversity in butterflies

Migration is typically associated with risk and uncertainty at the population level, but little is known about its cost-benefit trade-offs at the species level. Migratory insects in particular often exhibit strong demographic fluctuations due to local bottlenecks and outbreaks. Here, we use genomic data to investigate levels of heterozygosity and long-term population size dynamics in migratory insects, as an alternative to classical local and short-term approaches such as regional field monitoring. We analyse whole-genome sequences from 97 Lepidoptera species and show that individuals of migratory species have significantly higher levels of genome-wide heterozygosity, a proxy for effective population size, than do nonmigratory species. Also, we contribute whole-genome data for one of the most emblematic insect migratory species, the painted lady butterfly (Vanessa cardui), sampled across its worldwide distributional range. This species exhibits one of the highest levels of genomic heterozygosity described in Lepidoptera (2.95 ± 0.15%). Coalescent modelling (PSMC) shows historical demographic stability in V. cardui, and high effective population size estimates of 2-20 million individuals 10,000 years ago. The study reveals that the high risks associated with migration and local environmental fluctuations do not seem to decrease overall genetic diversity and demographic stability in migratory Lepidoptera. We propose a "compensatory" demographic model for migratory r-strategist organisms in which local bottlenecks are counterbalanced by reproductive success elsewhere within their typically large distributional ranges. Our findings highlight that the boundaries of populations are substantially different for sedentary and migratory insects, and that, in the latter, local and even regional field monitoring results may not reflect whole population dynamics. Genomic diversity patterns may elucidate key aspects of an insect's migratory nature and population dynamics at large spatiotemporal scales.

Rapid adaptation in a fast-changing world: Emerging insights from insect genomics

Many researchers have questioned the ability of biota to adapt to rapid anthropogenic environmental shifts. Here, we synthesize emerging genomic evidence for rapid insect evolution in response to human pressure. These new data reveal diverse genomic mechanisms (single locus, polygenic, structural shifts; introgression) underpinning rapid adaptive responses to a variety of anthropogenic selective pressures. While the effects of some human impacts (e.g. pollution; pesticides) have been previously documented, here we highlight startling new evidence for rapid evolutionary responses to additional anthropogenic processes such as deforestation. These recent findings indicate that diverse insect assemblages can indeed respond dynamically to major anthropogenic evolutionary challenges. Our synthesis also emphasizes the critical roles of genomic architecture, standing variation and gene flow in maintaining future adaptive potential. Broadly, it is clear that genomic approaches are essential for predicting, monitoring and responding to ongoing anthropogenic biodiversity shifts in a fast-changing world.

Genome Assembly and Comparative Analysis of the Egg Parasitoid Wasp Trichogramma dendrolimi Shed Light on the Composition and Evolution of Olfactory Receptors and Venoms

Trichogramma dendrolimi is one of the most successfully industrialized Trichogramma species used to control agricultural and forestry pests in China. However, the molecular mechanisms underlying its host recognition and parasitism remain largely unknown, partially due to the limited genome information of this parasitoid wasp. Here, we present a high-quality de novo assembly of T. dendrolimi through a combination of Illumina and PacBio sequencing technologies. The final assembly had a length of 215.2 Mb and contains 316 scaffolds with a scaffold N50 size of 1.41 Mb. Repetitive sequences with a length of 63.4 Mb and 12,785 protein-coding genes were identified. Significantly expanded gene families were identified to be involved in the development and regulatory processes, while remarkably contracted gene families were involved in the transport processes in T. dendrolimi. The olfactory and venom-associated genes were identified in T. dendrolimi and 24 other hymenopteran species, using uniform methods combining BLAST and HMM profiling. The identified venom genes of T. dendrolimi were enriched in antioxidant activity, tricarboxylic acid cycle, response to oxidative stress and cell redox homeostasis. Our study provides an important resource for comparative genomics and functional studies to interpret the molecular mechanisms underlying host recognition and parasitism of Trichogramma species.

Linkage mapping and genome annotation give novel insights into gene family expansions and regional recombination rate variation in the painted lady (Vanessa cardui) butterfly

Characterization of gene family expansions and crossing over is crucial for understanding how organisms adapt to the environment. Here, we develop a high-density linkage map and detailed genome annotation of the painted lady butterfly (Vanessa cardui) - a non-diapausing, highly polyphagous species famous for its long-distance migratory behavior and almost cosmopolitan distribution. Our results reveal a complex interplay between regional recombination rate variation, gene duplications and transposable element activity shaping the genome structure of the painted lady. We identify several lineage specific gene family expansions. Their functions are mainly associated with protein and fat metabolism, detoxification, and defense against infection - critical processes for the painted lady's unique life-history. Furthermore, the detailed recombination maps allow us to characterize the regional recombination landscape, data that reveal a strong effect of chromosome size on the recombination rate, a limited impact of GC-biased gene conversion and a positive association between recombination and short interspersed elements.

Molecular advances to study the function, evolution and spectral tuning of arthropod visual opsins

Visual opsins of vertebrates and invertebrates diversified independently and converged to detect ultraviolet to long wavelengths (LW) of green or red light. In both groups, colour vision largely derives from opsin number, expression patterns and changes in amino acids interacting with the chromophore. Functional insights regarding invertebrate opsin evolution have lagged behind those for vertebrates because of the disparity in genomic resources and the lack of robust in vitro systems to characterize spectral sensitivities. Here, we review bioinformatic approaches to identify and model functional variation in opsins as well as recently developed assays to measure spectral phenotypes. In particular, we discuss how transgenic lines, cAMP-spectroscopy and sensitive heterologous expression platforms are starting to decouple genotype-phenotype relationships of LW opsins to complement the classical physiological-behavioural-phylogenetic toolbox of invertebrate visual sensory studies. We illustrate the use of one heterologous method by characterizing novel LW Gq opsins from 10 species, including diurnal and nocturnal Lepidoptera, a terrestrial dragonfly and an aquatic crustacean, expressing them in HEK293T cells, and showing that their maximum absorbance spectra (λmax) range from 518 to 611 nm. We discuss the advantages of molecular approaches for arthropods with complications such as restricted availability, lateral filters, specialized photochemistry and/or electrophysiological constraints. This article is part of the theme issue 'Understanding colour vision: molecular, physiological, neuronal and behavioural studies in arthropods'.

A chromosome-level genome assembly provides insights into the environmental adaptability and outbreaks of Chlorops oryzae

Chlorops oryzae is a pest of rice that has caused severe damage to crops in major rice-growing areas in recent years. We generated a 447.60 Mb high-quality chromosome-level genome with contig and scaffold N50 values of 1.17 Mb and 117.57 Mb, respectively. Hi-C analysis anchored 93.22% scaffolds to 4 chromosomes. The relatively high expression level of Heat Shock Proteins (HSPs) and antioxidant genes in response to thermal stress suggests these genes may play a role in the environmental adaptability of C. oryzae. The identification of multiple pathways that regulate reproductive development (juvenile hormone, 20-hydroxyecdsone, and insulin signaling pathways) provides evidence that these pathways also play an important role in vitellogenesis and thus insect population maintenance. These findings identify possible reasons for the increased frequency of outbreaks of C. oryzae in recent years. Our chromosome-level genome assembly may provide a basis for further genetic studies of C. oryzae, and promote the development of novel, sustainable strategies to control this pest.

A chromosome-level genome assembly for the rice pest, Chlorops oryzae, pinpoints molecular pathways that might contribute toward increased outbreaks for this important crop pest.

Conferring the Midas Touch on Integrative Taxonomy: A Nanogold-Oligonucleotide Conjugate-Based Quick Species Identification Tool

Nanogold or functionalized gold nanoparticles (GNPs) have myriad applications in medical sciences. GNPs are widely used in the area of nanodiagnostics and nanotherapeutics. Applications of GNPs in taxonomic studies have not been studied vis-à-vis its extensive medical applications. GNPs have great potential in the area of integrative taxonomy. We have realized that GNPs can be used to visually detect animal species based on molecular signatures. In this regard, we have synthesized gold nanoparticles (<20 nm) and have developed a method based on interactions between thiolated DNA oligonucleotides and small-sized GNPs, interactions between DNA oligonucleotides and target DNA molecules, and self-aggregating properties of small-sized GNPs under high salt concentrations leading to a visible change in colour. Exploiting these intermolecular and interparticle interactions under aqueous conditions, in the present work, we have demonstrated the application of our procedure by using a DNA oligonucleotide probe designed against a portion of the mitochondrial genome of the codling moth Cydia pomonella. This method is accurate, quick, and easy to use once devised and can be used as an additional tool along with DNA barcoding. This tool can be used for distinguishing cryptic species, identification of morphovariants of known species, diet analysis, and identification of pest species in quarantine facilities without any need of performing repetitive DNA sequencing. We suggest that designing and selecting a highly specific DNA probe is crucial in increasing the specificity of the procedure. Present work may be considered as an effort to introduce nanotechnology as a new discipline to the extensive field of integrative taxonomy with which disciplines like palaeontology, embryology, anatomy, ethology, ecology, biochemistry, and molecular biology are already associated for a long time.