Genome of the small hive beetle (Aethina tumida, Coleoptera: Nitidulidae), a worldwide parasite of social bee colonies, provides insights into detoxification and herbivory

Comparative analyses of the banded alder borer (Rosalia funebris) and Asian longhorned beetle (Anoplophora glabripennis) genomes reveal significant differences in genome architecture and gene content among these and other Cerambycidae

Rosalia funebris (RFUNE; Cerambycidae), the banded alder borer, is a longhorn beetle whose larvae feed on the wood of various economically and ecologically significant trees in western North America. Adults are short-lived and not known to consume plant material substantially. We sequenced, assembled, and annotated the RFUNE genome using HiFi and RNASeq data. We documented genome architecture and gene content, focusing on genes putatively involved in plant feeding (phytophagy). Comparisons were made to the well-studied genome of the Asian longhorned beetle (AGLAB; Anoplophora glabripennis) and other Cerambycidae. The 814 Mb RFUNE genome assembly was distributed across 42 contigs, with an N50 of 30.18 Mb. Repetitive sequences comprised 60.27% of the genome, and 99.0% of expected single-copy orthologous genes were fully assembled. We identified 12,657 genes, fewer than in the four other species studied, and 46.4% fewer than for Aromia moschata (same subfamily as RFUNE). Of the 7,258 orthogroups shared between RFUNE and AGLAB, 1,461 had more copies in AGLAB and 1,023 had more copies in RFUNE. We identified 240 genes in RFUNE that putatively arose via horizontal transfer events. The RFUNE genome encoded substantially fewer putative plant cell wall degrading enzymes than AGLAB, which may relate to the longer-lived plant-feeding adults of the latter species. The RFUNE genome provides new insights into cerambycid genome architecture and gene content and provides a new vantage point from which to study the evolution and genomic basis of phytophagy in beetles.

Matrix metalloproteinase 2 contributes to adult eclosion and immune response in the small hive beetle, Aethina tumida

During the pupal-adult eclosion process of holometabolous insects, the old cuticle is shed and replaced by a completely different new cuticle that requires tanning and expansion, along with extensive extracellular matrix (ECM) remodeling. In vertebrates, matrix metalloproteinases (MMPs), a class of zinc-dependent endopeptidases, play key roles in regulating the ECM that surrounds cells. However, little is known about these extracellular proteinases available in insects. The small hive beetle (SHB), Aethina tumida, is a widespread invasive parasite of honey bees. In this study, 6 MMP homologs were identified in the SHB genome. RNA interference experiments showed that all 6 AtMmps are not required for the larval-pupal transition, only AtMmp2 was essential for pupal-adult eclosion in SHB. Knockdown of AtMmp2 resulted in eclosion defects and wing expansion failure, as well as mortality within 3 d of adult eclosion. Transcriptomic analysis revealed that knockdown of AtMmp2 significantly increased expression of the Toll and Imd pathways, chitin metabolism, and cross-linking (such as the pro-phenoloxidase activating cascade pathway and the tyrosine-mediated cuticle sclerotization and pigmentation pathway). These data revealed evolutionarily conserved functions of Mmp2 in controlling adult eclosion and wing expansion, also provided a preliminary exploration of the novel function of regulating Toll and Imd pathways, as well as new insights into how MMPs regulate insect development and defense barriers.

Biotic and abiotic stresses on honeybee health

Honeybees are the most critical pollinators providing key ecosystem services that underpin crop production and sustainable agriculture. Amidst a backdrop of rapid global change, this eusocial insect encounters a succession of stressors during nesting, foraging, and pollination. Ectoparasitic mites, together with vectored viruses, have been recognized as central biotic threats to honeybee health, while the spread of invasive giant hornets and small hive beetles also increasingly threatens colonies worldwide. Cocktails of agrochemicals, including acaricides used for mite treatment, and other pollutants of the environment have been widely documented to affect bee health in various ways. Additionally, expanding urbanization, climate change, and agricultural intensification often result in the destruction or fragmentation of flower-rich bee habitats. The anthropogenic pressures exerted by beekeeping management practices affect the natural selection and evolution of honeybees, and colony translocations facilitate alien species invasion and disease transmission. In this review, the multiple biotic and abiotic threats and their interactions that potentially undermine bee colony health are discussed, while taking into consideration the sensitivity, large foraging area, dense network among related nestmates, and social behaviors of honeybees.

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

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

Thiamethoxam soil contaminations reduce fertility of soil-dwelling beetles, Aethina tumida

There in increasing evidence for recent global insect declines. This is of major concern as insects play a critical role in ecosystem functionality and human food security. Even though environmental pollutants are known to reduce insect fertility, their potential effects on insect fitness remain poorly understood - especially for soil-dwelling species. Here, we show that fertility of soil-dwelling beetles, Aethina tumida, is reduced, on average, by half due to field-realistic neonicotinoid soil contaminations. In the laboratory, pupating beetles were exposed via soil to concentrations of the neonicotinoid thiamethoxam that reflect global pollution of agricultural and natural habitats. Emerged adult phenotypes and reproduction were measured, and even the lowest concentration reported from natural habitats reduced subsequent reproduction by 50%. The data are most likely a conservative estimate as the beetles were only exposed during pupation. Since the tested concentrations reflect ubiquitous soil pollution, the data reveal a plausible mechanism for ongoing insect declines. An immediate reduction in environmental pollutants is urgently required if our aim is to mitigate the prevailing loss of species biodiversity.

Tracking Adaptive Pathways of Invasive Insects: Novel Insight from Genomics

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

Multiple mating by both sexes in an invasive insect species, Aethina tumida (Coleoptera: Nitidulidae)

Multiple mating by both sexes is common among sexually reproducing animals. Small hive beetles (SHB), Aethina tumida, are parasites of bee nests endemic to sub-Saharan Africa and have become a widespread invasive species. Despite the considerable economic damages they can cause, their basic biology remains poorly understood. Here we show that male and female small hive beetles can mate multiple times, suggesting that costs for mating are low in this species. In an invasive A. tumida population in the United States, a combination of laboratory experiments for males and paternity analysis with eight polymorphic DNA microsatellite markers for field-caught females were used to estimate the number of mating by both sexes. The data show that females and males can mate multiple times-females mated with up to eight males, whereas males mated with at least seven females. The results also showed that A. tumida displayed a skewed paternity, although this was not consistent among the tested females. Thus, first or last male advantage seem to be unlikely in A. tumida. Our observations that individuals of both sexes of A. tumida can mate multiple times opens new research avenues for examining drivers of multiple mating and determining the role it may play in promoting biological invasions.

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

BACKGROUND

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

RESULTS

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

CONCLUSION

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

Role of Insect and Mammal Glutathione Transferases in Chemoperception

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

Identification and characterization of CYP307A1 as a molecular target for controlling the small hive beetle, Aethina tumida

BACKGROUND

The molting hormone 20-hydroxyecdysone (20E) plays a key role in insect development, metamorphosis, and reproduction. Previous studies have shown that ecdysteroid metabolism is regulated by a series of CYP genes in most of the insect species. However, the roles of these CYP genes in a Coleopteran beetle, Aethina tumida (small hive beetle, SHB) have not yet been explored.

RESULTS

In the current study, we identified seven CYP genes (six Halloween genes and one AtCYP18A1 gene) related to 20E metabolism. Reverse transcription quantitative polymerase chain reaction (RT-qPCR) showed that AtCYP307A1 and AtCYP307B1 were primarily expressed in the embryonic stage and in the cephalothorax of larvae. RNA interference (RNAi) screening revealed that suppression of AtCYP307A1 expression caused a lethal phenotype during the larval-pupal metamorphosis. Furthermore, Hematoxylin and Eosin staining of the integument showed that the RNAi of AtCYP307A1 inhibited the apolysis and degradation of the old cuticle. In addition, silencing of AtCYP307A1 resulted in significant down-regulation of 20E titers and the expression levels of 20E signaling pathway genes. Finally, the AtCYP307A1 RNAi phenotype was rescued by topical application of 20E.

CONCLUSION

Our studies suggest that AtCYP307A1 involved in 20E synthesis is indispensable during the larval-pupal metamorphosis of beetles, which could serve as a putative insecticide target for pest control. © 2022 Society of Chemical Industry.

Identification of sex chromosomes and primary sex ratio in the small hive beetle, a worldwide parasite of honey bees

BACKGROUND

The small hive beetle (SHB), Aethina tumida, has emerged as a worldwide threat to honey bees in the past two decades. These beetles harvest nest resources, feed on larval bees, and ultimately spoil nest resources with gelatinous slime together with the fungal symbiont Kodamaea ohmeri.

RESULTS

Here, we present the first chromosome-level genome assembly for the SHB. With a 99.1% representation of conserved (BUSCO) arthropod genes, this resource enables the study of chemosensory, digestive, and detoxification traits critical for SHB success and possible control. We use this annotated assembly to characterize features of SHB sex chromosomes and a female-skewed primary sex ratio. We also found chromosome fusion and a lower recombination rate in sex chromosomes than in autosomes.

CONCLUSIONS

Genome-enabled insights will clarify the traits that allowed this beetle to exploit hive resources successfully and will be critical for determining the causes of observed sex ratio asymmetries.

Genome-wide and expression-profiling analyses of the cytochrome P450 genes in Tenebrionidea

Cytochrome P450 monooxygenases (CYPs) are present in almost all areas of the tree of life. As one of the largest and most diverse superfamilies of multifunctional enzymes, they play important roles in the metabolism of xenobiotics and biosynthesis of endogenous compounds, shaping the success of insects. In this study, the CYPome (an omics term for all the CYP genes in a genome) diversification was examined in the four Tenebrionidea species through genome-wide analysis. A total of 483 CYP genes were identified, of which 103, 157, 122, and 101 were respectively deciphered from the genomes of Tebebrio molitor, Asbolus verucosus, Hycleus cichorii and Hycleus phaleratus. These CYPs were classified into four major clans (mitochondrial, CYP2, CYP3, and CYP4), and clans CYP3 and CYP4 are most diverse. Phylogenetic analysis showed that most CYPs of these Tenebrionidea beetles from each clan had a very close 1:1 orthology to each other, suggesting that they originate closely and have evolutionally conserved function. Expression analysis at different developmental stages and in various tissues showed the life stage-, gut-, salivary gland-, fat body-, Malpighian tubule-, antennae-, ovary- and testis-specific expression patterns of T. molitor CYP genes, implying their various potential roles in development, detoxification, immune response, digestion, olfaction, and reproduction. Our studies provide a platform to understand the evolution of Tenebrionidea CYP gene superfamily, and a basis for further functional investigation of the T. molitor CYPs involved in various biological processes.

Comparative genomic analysis of carboxylesterase genes in Tenebrio molitor and other four tenebrionids

Carboxylesterases (COEs) have various functions in wide taxons of organisms. In insects, COEs are important enzymes involved in the hydrolysis of a variety of ester-containing xenobiotics, neural signal transmission, pheromone degradation, and reproductive development. Understanding the diversity of COEs is basic to illustrate their functions. In this study, we identified 53, 105, 37, and 39 COEs from the genomes of Tenebrio molitor, Asbolus verucosus, Hycleus cichorii, and H. phaleratus in the superfamily of Tenebrionidea, respectively. Phylogenetic analysis showed that 234 COEs from these four species and those reported in Tribolium castaneum (63) could be divided into 12 clades and three major classes. The α-esterases significantly expanded in T. molitor, A. verucosus, and T. castaneum compared to dipteran and hymenopteran insects. In T. molitor, most COEs showed tissue and stage-specific but not a sex-biased expression. Our results provide insights into the diversity and evolutionary characteristics of COEs in tenebrionids, and lay a foundation for the functional characterization of COEs in the yellow mealworm.

Comparative genomic analysis of ABC transporter genes in Tenebrio molitor and four other tenebrionid beetles (Coleoptera: Tenebrionidea)

ATP-binding cassette (ABC) transporters, one of the largest transmembrane protein families, transport a diverse number of substate across membranes. Details of their diverse physiological functions have not been established. Here, we identified 87 ABC transporter genes in the genomes of Tenebrio molitor along with those from Asbolus verrucosus (104), Hycleus cichorii (65), and Hycleus phaleratus (80). Combining these genes (336 in total) with genes reported in Tribolium castaneum (73), we analyzed the phylogeny of ABC transporter genes in all five Tenebrionids. They are assigned into eight subfamilies (ABCA-H). In comparison to other species, the ABCC subfamily in this group of beetles appears expanded. The expression profiles of the T. molitor genes at different life stages and in various tissues were also investigated using transcriptomic analysis. Most of them display developmental specific expression patterns, suggesting to us their possible roles in development. Most of them are highly expressed in detoxification-related tissues including gut and Malpighian tubule, from which we infer their roles in insecticide resistance. We detected specific or abundant expressions of many ABC transporter genes in various tissues such as salivary gland, ovary, testis, and antenna. This new information helps generate new hypotheses on their biological significance within tissues.

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

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

A chromosome-level genome assembly and intestinal transcriptome of Trypoxylus dichotomus (Coleoptera: Scarabaeidae) to understand its lignocellulose digestion ability

Lignocellulose, as the key structural component of plant biomass, is a recalcitrant structure, difficult to degrade. The traditional management of plant waste, including landfill and incineration, usually causes serious environmental pollution and health problems. Interestingly, the xylophagous beetle, Trypoxylus dichotomus, can decompose lignocellulosic biomass. However, the genomics around the digestion mechanism of this beetle remain to be elucidated. Here, we assembled the genome of T. dichotomus, showing that the draft genome size of T. dichotomus is 636.27 Mb, with 95.37% scaffolds anchored onto 10 chromosomes. Phylogenetic results indicated that a divergent evolution between the ancestors of T. dichotomus and the closely related scarabaeid species Onthophagus taurus occurred in the early Cretaceous (120 million years ago). Through gene family evolution analysis, we found 67 rapidly evolving gene families, within which there were 2 digestive gene families (encoding Trypsin and Enoyl-(Acyl carrier protein) reductase) that have experienced significant expansion, indicating that they may contribute to the high degradation efficiency of lignocellulose in T. dichotomus. Additionally, events of chromosome breakage and rearrangement were observed by synteny analysis during the evolution of T. dichotomus due to chromosomes 6 and 8 of T. dichotomus being intersected with chromosomes 2 and 10 of Tribolium castaneum, respectively. Furthermore, the comparative transcriptome analyses of larval guts showed that the digestion-related genes were more commonly expressed in the midgut or mushroom residue group than the hindgut or sawdust group. This study reports the well-assembled and annotated genome of T. dichotomus, providing genomic and transcriptomic bases for further understanding the functional and evolutionary mechanisms of lignocellulose digestion in T. dichotomus.

Phylogenetic Analysis of Small Hive Beetles From Native to Introduced Populations

The small hive beetle (SHB), a social parasite of beehives, is native to sub-Saharan Africa and has spread to America, Europe, and Australia. Recently, these beetles invaded China, causing widespread colony collapses in the honeybee, Apis cerana. In this study, single nucleotide polymorphisms (SNPs) were identified in the beetle genome from its native range (Africa), a region that was invaded by SHBs nearly 30 years ago (America), and more recent invasions (Asia). The beetles in the United States formed the earliest branch and show signs of two decades of gene flow and local adaptation to differentiate this population from the native ones. The beetles in China were deep branched and showed the highest fixation index when compared to the US populations. The number of SNPs in overexpressed genes was significantly higher than the transcriptome. Gene-expression profiles presented here distinguish the characters between adult and larvae SHBs.

Molecular Detection and Differentiation of Arthropod, Fungal, Protozoan, Bacterial and Viral Pathogens of Honeybees

The honeybee Apis mellifera is highly appreciated worldwide because of its products, but also as it is a pollinator of crops and wild plants. The beehive is vulnerable to infections due to arthropods, fungi, protozoa, bacteria and/or viruses that manage to by-pass the individual and social immune mechanisms of bees. Due to the close proximity of bees in the beehive and their foraging habits, infections easily spread within and between beehives. Moreover, international trade of bees has caused the global spread of infections, several of which result in significant losses for apiculture. Only in a few cases can infections be diagnosed with the naked eye, by direct observation of the pathogen in the case of some arthropods, or by pathogen-associated distinctive traits. Development of molecular methods based on the amplification and analysis of one or more genes or genomic segments has brought significant progress to the study of bee pathogens, allowing for: (i) the precise and sensitive identification of the infectious agent; (ii) the analysis of co-infections; (iii) the description of novel species; (iv) associations between geno- and pheno-types and (v) population structure studies. Sequencing of bee pathogen genomes has allowed for the identification of new molecular targets and the development of specific genotypification strategies.

The Genome of Rhyzopertha dominica (Fab.) (Coleoptera: Bostrichidae): Adaptation for Success

The lesser grain borer, Rhyzopertha dominica (F.) (Coleoptera: Bostrichidae), is a major global pest of cereal grains. Infestations are difficult to control as larvae feed inside grain kernels, and many populations are resistant to both contact insecticides and fumigants. We sequenced the genome of R. dominica to identify genes responsible for important biological functions and develop more targeted and efficacious management strategies. The genome was assembled from long read sequencing and long-range scaffolding technologies. The genome assembly is 479.1 Mb, close to the predicted genome size of 480.4 Mb by flow cytometry. This assembly is among the most contiguous beetle assemblies published to date, with 139 scaffolds, an N₅₀ of 53.6 Mb, and L₅₀ of 4, indicating chromosome-scale scaffolds. Predicted genes from biologically relevant groups were manually annotated using transcriptome data from adults and different larval tissues to guide annotation. The expansion of carbohydrase and serine peptidase genes suggest that they combine to enable efficient digestion of cereal proteins. A reduction in the copy number of several detoxification gene families relative to other coleopterans may reflect the low selective pressure on these genes in an insect that spends most of its life feeding internally. Chemoreceptor genes contain elevated numbers of pseudogenes for odorant receptors that also may be related to the recent ontogenetic shift of R. dominica to a diet consisting primarily of stored grains. Analysis of repetitive sequences will further define the evolution of bostrichid beetles compared to other species. The data overall contribute significantly to coleopteran genetic research.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Comparative genomics suggests local adaptations in the invasive small hive beetle

Invasive species are a major driver of ecological and environmental changes that affect human health, food security, and natural biodiversity. The success and impact of biological invasions depend on adaptations to novel abiotic and biotic selective pressures. However, the molecular mechanisms underlying adaptations in invasive parasitic species are inadequately understood. Small hive beetles, Aethina tumida, are parasites of bee nests. Originally endemic to sub-Saharan Africa, they are now found nearly globally. Here, we investigated the molecular bases of the adaptations to novel environments underlying their invasion routes. Genomes of historic and recent adults A. tumida from both the endemic and introduced ranges were compared. Analysis of gene-environment association identified 3049 candidate loci located in 874 genes. Functional annotation showed a significant bias toward genes linked to growth and reproduction. One of the genes from the apoptosis pathway encodes an "ecdysone-related protein," which is a crucial regulator in controlling body size in response to environmental cues for holometabolous insects during cell death and renewal. Genes whose proteins regulate organ size, ovary activation, and oviposition were also detected. Functions of these enriched pathways parallel behavioral differences between introduced and native A. tumida populations, which may reflect patterns of local adaptation. The results considerably improve our understanding of the underlying mechanisms and ecological factors driving adaptations of invasive species. Deep functional investigation of these identified loci will help clarify the mechanisms of local adaptation in A. tumida.

Genome Assembly of the Cold-Tolerant Leaf Beetle Gonioctena quinquepunctata, an Important Resource for Studying Its Evolution and Reproductive Barriers between Species

Coleoptera is the most species-rich insect order, yet is currently underrepresented in genomic databases. An assembly was generated for ca. 1.7 Gb genome of the leaf beetle Gonioctena quinquepunctata by first assembling long-sequence reads (Oxford Nanopore; ± 27-fold coverage) and subsequently polishing the resulting assembly with short sequence reads (Illumina; ± 85-fold coverage). The unusually large size (most Coleoptera species are associated with a reported size below 1 Gb) was at least partially attributed to the presence of a large fraction of repeated elements (73.8%). The final assembly was characterized by an N50 length of 432 kb and a BUSCO score of 95.5%. The heterozygosity rate was ± 0.6%. Automated genome annotation informed by RNA-Seq resulted in 40,568 predicted proteins, which is much larger than the typical range 17,000–23,000 predicted for other Coleoptera. However, no evidence of a genome duplication was detected. This new reference genome will contribute to our understanding of genetic variation in the Coleoptera. Among others, it will also allow exploring reproductive barriers between species, investigating introgression in the nuclear genome, and identifying genes involved in resistance to extreme climate conditions.

A high-quality carabid genome assembly provides insights into beetle genome evolution and cold adaptation

The hyperdiverse order Coleoptera comprises a staggering ~25% of known species on Earth. Despite recent breakthroughs in next generation sequencing, there remains a limited representation of beetle diversity in assembled genomes. Most notably, the ground beetle family Carabidae, comprising more than 40,000 described species, has not been studied in a comparative genomics framework using whole genome data. Here we generate a high-quality genome assembly for Nebria riversi, to examine sources of novelty in the genome evolution of beetles, as well as genetic changes associated with specialization to high-elevation alpine habitats. In particular, this genome resource provides a foundation for expanding comparative molecular research into mechanisms of insect cold adaptation. Comparison to other beetles shows a strong signature of genome compaction, with N. riversi possessing a relatively small genome (~147 Mb) compared to other beetles, with associated reductions in repeat element content and intron length. Small genome size is not, however, associated with fewer protein-coding genes, and an analysis of gene family diversity shows significant expansions of genes associated with cellular membranes and membrane transport, as well as protein phosphorylation and muscle filament structure. Finally, our genomic analyses show that these high-elevation beetles have endosymbiotic Spiroplasma, with several metabolic pathways (e.g., propanoate biosynthesis) that might complement N. riversi, although its role as a beneficial symbiont or as a reproductive parasite remains equivocal.

Evolution: Symbiotic Microbes Mediate Host Range of Herbivorous Beetles

A new study characterizes the evolution of metabolic interactions between herbivorous tortoise beetles and their obligate bacterial symbionts. Expansion of symbiont metabolic range is shown to alter beetle metabolism of plant tissues and is implicated in broadening the diversity of plants eaten.

The draft genome of the specialist flea beetle Altica viridicyanea (Coleoptera: Chrysomelidae)

BACKGROUND

Altica (Coleoptera: Chrysomelidae) is a highly diverse and taxonomically challenging flea beetle genus that has been used to address questions related to host plant specialization, reproductive isolation, and ecological speciation. To further evolutionary studies in this interesting group, here we present a draft genome of a representative specialist, Altica viridicyanea, the first Alticinae genome reported thus far.

RESULTS

The genome is 864.8 Mb and consists of 4490 scaffolds with a N50 size of 557 kb, which covered 98.6% complete and 0.4% partial insect Benchmarking Universal Single-Copy Orthologs. Repetitive sequences accounted for 62.9% of the assembly, and a total of 17,730 protein-coding gene models and 2462 non-coding RNA models were predicted. To provide insight into host plant specialization of this monophagous species, we examined the key gene families involved in chemosensation, detoxification of plant secondary chemistry, and plant cell wall-degradation.

CONCLUSIONS

The genome assembled in this work provides an important resource for further studies on host plant adaptation and functionally affiliated genes. Moreover, this work also opens the way for comparative genomics studies among closely related Altica species, which may provide insight into the molecular evolutionary processes that occur during ecological speciation.

Horizontally acquired antibacterial genes associated with adaptive radiation of ladybird beetles

BACKGROUND

Horizontal gene transfer (HGT) has been documented in many herbivorous insects, conferring the ability to digest plant material and promoting their remarkable ecological diversification. Previous reports suggest HGT of antibacterial enzymes may have contributed to the insect immune response and limit bacterial growth. Carnivorous insects also display many evolutionary successful lineages, but in contrast to the plant feeders, the potential role of HGTs has been less well-studied.

RESULTS

Using genomic and transcriptomic data from 38 species of ladybird beetles, we identified a set of bacterial cell wall hydrolase (cwh) genes acquired by this group of beetles. Infection with Bacillus subtilis led to upregulated expression of these ladybird cwh genes, and their recombinantly produced proteins limited bacterial proliferation. Moreover, RNAi-mediated cwh knockdown led to downregulation of other antibacterial genes, indicating a role in antibacterial immune defense. cwh genes are rare in eukaryotes, but have been maintained in all tested Coccinellinae species, suggesting that this putative immune-related HGT event played a role in the evolution of this speciose subfamily of predominant predatory ladybirds.

CONCLUSION

Our work demonstrates that, in a manner analogous to HGT-facilitated plant feeding, enhanced immunity through HGT might have played a key role in the prey adaptation and niche expansion that promoted the diversification of carnivorous beetle lineages. We believe that this represents the first example of immune-related HGT in carnivorous insects with an association with a subsequent successful species radiation.

Genome-wide identification of the aspartic protease gene family and their response under powdery mildew stress in wheat

Aspartic proteases (APs) are one of the four main protease super families. In plants, they are involved in many biological processes, such as biotic and abiotic stress resistance, protein processing and degradation, senescence, and programmed cell death. By performing a database (TGACv1) search and domain prediction, we identified 263 wheat AP (TaAP) proteins and observed 38 TaAP genes exhibiting alternative splicing. Moreover, by constructing a phylogenetic tree, we found that the TaAP proteins can be divided into three families and have a certain close evolutionary relationship to Arabidopsis thaliana and rice AP proteins. Transcriptome analysis showed that 29 genes in the TaAP family were up-regulated after being induced by powdery mildew. The expression of TaAP224 showed the most significant difference in transcriptome and qRT-PCR analyses. Subsequently, the promoters of these 29 genes were analysed, and we found that they contained multiple disease resistance and hormone elements, such as WRKY71OS, a common disease resistance element that is also involved in the GA signalling pathway and inhibits starch hydrolysis. The comprehensive annotation and expression profiling performed in this study increased our understanding of the TaAP family genes in wheat growth and development, and the results can be used as a basis for further study of candidate TaAP genes involved in powdery mildew resistance mechanisms.

Multi-faceted analysis provides little evidence for recurrent whole-genome duplications during hexapod evolution

BACKGROUND

Gene duplication events play an important role in the evolution and adaptation of organisms. Duplicated genes can arise through different mechanisms, including whole-genome duplications (WGDs). Recently, WGD was suggested to be an important driver of evolution, also in hexapod animals.

RESULTS

Here, we analyzed 20 high-quality hexapod genomes using whole-paranome distributions of estimated synonymous distances (KS), patterns of within-genome co-linearity, and phylogenomic gene tree-species tree reconciliation methods. We observe an abundance of gene duplicates in the majority of these hexapod genomes, yet we find little evidence for WGD. The majority of gene duplicates seem to have originated through small-scale gene duplication processes. We did detect segmental duplications in six genomes, but these lacked the within-genome co-linearity signature typically associated with WGD, and the age of these duplications did not coincide with particular peaks in KS distributions. Furthermore, statistical gene tree-species tree reconciliation failed to support all but one of the previously hypothesized WGDs.

CONCLUSIONS

Our analyses therefore provide very limited evidence for WGD having played a significant role in the evolution of hexapods and suggest that alternative mechanisms drive gene duplication events in this group of animals. For instance, we propose that, along with small-scale gene duplication events, episodes of increased transposable element activity could have been an important source for gene duplicates in hexapods.

Identification and characterization of detoxification genes in two cerambycid beetles, Rhaphuma horsfieldi and Xylotrechus quadripes (Coleoptera: Cerambycidae: Clytini)

The longhorned beetles, Rhaphuma horsfieldi and Xylotrechus quadripes, are two polyphagous insects with larvae feeding on different host plants. In this study, we identified and characterized three gene superfamilies of cytochrome P450s (CYPs), carboxylesterases (COEs) and glutathione-S-transferases (GSTs) involved in the detoxification of endobiotics (e.g., hormones and steroids) and xenobiotics (e.g., insecticides, sex pheromones and plant allelochemicals) through a combination approach of bioinformatics, phylogenetics, expression profiles and genomics. Transcriptome analyses led to the identification of 281 transcripts encoding 135 P450s, 108 COEs and 38 GSTs from the two beetles, coupled with comparative studies of detoxification genes among coleopteran species, suggesting a correlation between host range and the sizes of P450 or COE gene repertoires. The P450s of two beetles were phylogenetically classified into four clades, representing the majority of genes in the CYP3 clan. The COEs from R. horsfieldi and X. quadripes were separately grouped into 11 and 10 clades, and the GST superfamily was assigned into six clades. Expression profiles revealed that the detoxification genes were broadly expressed in various tissues as an implication of functional diversities. Ultimately and more importantly, five alternative splicing events in the Epsilon GSTs, including RhorGSTe7.1/GSTe7.2 and XquaGSTe3.1/GST3.2, were acquired in Coleoptera, in which these genes and their orthologs shared highly conserved gene structure. Our current study has complemented the resources for the detoxification genes in the family Cerambycidae, and allows for functional experiments to identify candidate molecular targets involved in pest resistance to insecticides like organophosphates, organochlorines and pyrethroids.

In silico identification and assessment of insecticide target sites in the genome of the small hive beetle, Aethina tumida

BACKGROUND

The small hive beetle, Aethina tumida, is a rapidly emerging global pest of honey bee colonies. Small hive beetle infestation can be extremely destructive, which may cause honey bees to abscond and render colony infrastructure unusable. Due to the impacts small hive beetles have on honey bees, a wide variety of physical, cultural, and chemical control measures have been implemented to manage small hive beetle infestations. The use of insecticides to control small hive beetle populations is an emerging management tactic. Currently, very little genomic information exists on insecticide target sites in the small hive beetle. Therefore, the objective of this study is to utilize focused in silico comparative genomics approaches to identify and assess the potential insecticide sensitivity of the major insecticide target sites in the small hive beetle genome.

RESULTS

No previously described resistance mutations were identified in any orthologs of insecticide target sites. Alternative exon use and A-to-I RNA editing were absent in AtumSC1. The ryanodine receptor in small hive beetle (Atum_Ryr) was highly conserved and no previously described resistance mutations were identified. A total of 12 nAChR subunits were identified with similar alternative exon use in other insects. Alternative exon use and critical structural features of the GABA-gated chloride channel subunits (Atum_RDL, Atum_GRD, and Atum_LCCH3) were conserved. Five splice variants were found for the glutamate-gated chloride channel subunit. Exon 3c of Atum_GluCl may be a beetle-specific alternative exon. The co-occurrence of exons 9a and 9b in the pH-sensitive chloride channel (Atum_pHCl) is a unique combination that introduces sites of post-translational modification. The repertoire and alternative exon use for histamine-gated chloride channels (Atum-HisCl), octopamine (Atum_OctR) and tyramine receptors (Atum_TAR) were conserved.

CONCLUSIONS

The recently published small hive beetle genome likely serves as a reference for insecticide-susceptible versions of insecticide target sites. These comparative in silico studies are the first step in discovering targets that can be exploited for small hive beetle-specific control as well as tracking changes in the frequency of resistance alleles as part of a resistance monitoring program. Comparative toxicity alongside honey bees is required to verify these in silico predictions.

Draft Genome of a Blister Beetle Mylabris aulica

Mylabris aulica is a widely distributed blister beetle of the Meloidae family. It has the ability to synthesize a potent defensive secretion that includes cantharidin, a toxic compound used to treat many major illnesses. However, owing to the lack of genetic studies on cantharidin biosynthesis in M. aulica, the commercial use of this species is less extensive than that of other blister beetle species in China. This study reports a draft assembly and possible genes and pathways related to cantharidin biosynthesis for the M. aulica blister beetle using nanopore sequencing data. The draft genome assembly size was 288.5 Mb with a 467.8 Kb N50, and a repeat content of 50.62%. An integrated gene finding pipeline performed for assembly obtained 16,500 protein coding genes. Benchmarking universal single-copy orthologs assessment showed that this gene set included 94.4% complete Insecta universal single-copy orthologs. Over 99% of these genes were assigned functional annotations in the gene ontology, Kyoto Encyclopedia of Genes and Genomes, or Genbank non-redundant databases. Comparative genomic analysis showed that the completeness and continuity of our assembly was better than those of Hycleus cichorii and Hycleus phaleratus blister beetle genomes. The analysis of homologous orthologous genes and inference from evolutionary history imply that the Mylabris and Hycleus genera are genetically close, have a similar genetic background, and have differentiated within one million years. This M. aulica genome assembly provides a valuable resource for future blister beetle studies and will contribute to cantharidin biosynthesis.

The evolution and genomic basis of beetle diversity

The order Coleoptera (beetles) is arguably the most speciose group of animals, but the evolutionary history of beetles, including the impacts of plant feeding (herbivory) on beetle diversification, remain poorly understood. We inferred the phylogeny of beetles using 4,818 genes for 146 species, estimated timing and rates of beetle diversification using 89 genes for 521 species representing all major lineages and traced the evolution of beetle genes enabling symbiont-independent digestion of lignocellulose using 154 genomes or transcriptomes. Phylogenomic analyses of these uniquely comprehensive datasets resolved previously controversial beetle relationships, dated the origin of Coleoptera to the Carboniferous, and supported the codiversification of beetles and angiosperms. Moreover, plant cell wall-degrading enzymes (PCWDEs) obtained from bacteria and fungi via horizontal gene transfers may have been key to the Mesozoic diversification of herbivorous beetles-remarkably, both major independent origins of specialized herbivory in beetles coincide with the first appearances of an arsenal of PCWDEs encoded in their genomes. Furthermore, corresponding (Jurassic) diversification rate increases suggest that these novel genes triggered adaptive radiations that resulted in nearly half of all living beetle species. We propose that PCWDEs enabled efficient digestion of plant tissues, including lignocellulose in cell walls, facilitating the evolution of uniquely specialized plant-feeding habits, such as leaf mining and stem and wood boring. Beetle diversity thus appears to have resulted from multiple factors, including low extinction rates over a long evolutionary history, codiversification with angiosperms, and adaptive radiations of specialized herbivorous beetles following convergent horizontal transfers of microbial genes encoding PCWDEs.

Structural and functional insights into the Diabrotica virgifera virgifera ATP-binding cassette transporter gene family

BACKGROUND

The western corn rootworm, Diabrotica virgifera virgifera, is a pervasive pest of maize in North America and Europe, which has adapted to current pest management strategies. In advance of an assembled and annotated D. v. virgifera genome, we developed transcriptomic resources to use in identifying candidate genes likely to be involved in the evolution of resistance, starting with members of the ATP-binding cassette (ABC) transporter family.

RESULTS

In this study, 65 putative D. v. virgifera ABC (DvvABC) transporters were identified within a combined transcriptome assembly generated from embryonic, larval, adult male, and adult female RNA-sequence libraries. Phylogenetic analysis placed the deduced amino-acid sequences of the DvvABC transporters into eight subfamilies (A to H). To supplement our sequence data with functional analysis, we identified orthologs of Tribolium castaneum ABC genes which had previously been shown to exhibit overt RNA interference (RNAi) phenotypes. We identified eight such D. v. virgifera genes, and found that they were functionally similar to their T. castaneum counterparts. Interestingly, depletion of DvvABCB_39715 and DvvABCG_3712 transcripts in adult females produced detrimental reproductive and developmental phenotypes, demonstrating the potential of these genes as targets for RNAi-mediated insect control tactics.

CONCLUSIONS

By combining sequence data from four libraries covering three distinct life stages, we have produced a relatively comprehensive de novo transcriptome assembly for D. v. virgifera. Moreover, we have identified 65 members of the ABC transporter family and provided the first insights into the developmental and physiological roles of ABC transporters in this pest species.

Identification of Aethina tumida Kir Channels as Putative Targets of the Bee Venom Peptide Tertiapin Using Structure-Based Virtual Screening Methods

Venoms are comprised of diverse mixtures of proteins, peptides, and small molecules. Identifying individual venom components and their target(s) with mechanism of action is now attainable to understand comprehensively the effectiveness of venom cocktails and how they collectively function in the defense and predation of an organism. Here, structure-based computational methods were used with bioinformatics tools to screen and identify potential biological targets of tertiapin (TPN), a venom peptide from Apis mellifera (European honey bee). The small hive beetle (Aethina tumida (A. tumida)) is a natural predator of the honey bee colony and was found to possess multiple inwardly rectifying K⁺ (Kir) channel subunit genes from a genomic BLAST search analysis. Structure-based virtual screening of homology modelled A. tumida Kir (atKir) channels found TPN to interact with a docking profile and interface “footprint” equivalent to known TPN-sensitive mammalian Kir channels. The results support the hypothesis that atKir channels, and perhaps other insect Kir channels, are natural biological targets of TPN that help defend the bee colony from infestations by blocking K⁺ transport via atKir channels. From these in silico findings, this hypothesis can now be subsequently tested in vitro by validating atKir channel block as well as in vivo TPN toxicity towards A. tumida. This study highlights the utility and potential benefits of screening in virtual space for venom peptide interactions and their biological targets, which otherwise would not be feasible.

Impact of Biotic and Abiotic Stressors on Managed and Feral Bees

Large-scale declines in bee abundance and species richness over the last decade have sounded an alarm, given the crucial pollination services that bees provide. Population dips have specifically been noted for both managed and feral bee species. The simultaneous increased cultivation of bee-dependent agricultural crops has given rise to additional concern. As a result, there has been a surge in scientific research investigating the potential stressors impacting bees. A group of environmental and anthropogenic stressors negatively impacting bees has been isolated. Habitat destruction has diminished the availability of bee floral resources and nest habitats, while massive monoculture plantings have limited bee access to a variety of pollens and nectars. The rapid spread and increased resistance buildup of various bee parasites, pathogens, and pests to current control methods are implicated in deteriorating bee health. Similarly, many pesticides that are widely applied on agricultural crops and within beehives are toxic to bees. The global distribution of honey bee colonies (including queens with attendant bees) and bumble bee colonies from crop to crop for pollination events has been linked with increased pathogen stress and increased competition with native bee species for limited resources. Climatic alterations have disrupted synchronous bee emergence with flower blooming and reduced the availability of diverse floral resources, leading to bee physiological adaptations. Interactions amongst multiple stressors have created colossal maladies hitting bees at one time, and in some cases delivering additive impacts. Initiatives including the development of wild flower plantings and assessment of pesticide toxicity to bees have been undertaken in efforts to ameliorate current bee declines. In this review, recent findings regarding the impact of these stressors on bees and strategies for mitigating them are discussed.

Shared and unique microbes between Small hive beetles (Aethina tumida) and their honey bee hosts

The small hive beetle (SHB) is an opportunistic parasite that feeds on bee larvae, honey, and pollen. While SHBs can also feed on fruit and other plant products, like its plant‐feeding relatives, SHBs prefer to feed on hive resources and only reproduce inside bee colonies. As parasites, SHBs are inevitably exposed to bee‐associated microbes, either directly from the bees or from the hive environment. These microbes have unknown impacts on beetles, nor is it known how extensively beetles transfer microbes among their bee hosts. To identify sets of beetle microbes and the transmission of microbes from bees to beetles, a metagenomic analysis was performed. We identified sets of herbivore‐associated bacteria, as well as typical bee symbiotic bacteria for pollen digestion, in SHB larvae and adults. Deformed wing virus was highly abundant in beetles, which colonize SHBs as suggested by a controlled feeding trial. Our data suggest SHBs are vectors for pathogen transmission among bees and between colonies. The dispersal of host pathogens by social parasites via floral resources and the hive environment increases the threats of these parasites to honey bees.

Draft Genome Sequence of the Yeast Kodamaea ohmeri, a Symbiont of the Small Hive Beetle

Kodamaea ohmeri is a symbiont of the small hive beetle (SHB), which is a scavenger of honey bee colonies. The SHB causes absconding of the economically important honey bee (Apis mellifera) and deposits K. ohmeri in the honeycomb. We describe long-read sequencing and further analyses of the K. ohmeri genome.

Genome of the small hive beetle (Aethina tumida, Coleoptera: Nitidulidae), a worldwide parasite of social bee colonies, provides insights into detoxification and herbivory

Background

The small hive beetle (Aethina tumida; ATUMI) is an invasive parasite of bee colonies. ATUMI feeds on both fruits and bee nest products, facilitating its spread and increasing its impact on honey bees and other pollinators. We have sequenced and annotated the ATUMI genome, providing the first genomic resources for this species and for the Nitidulidae, a beetle family that is closely related to the extraordinarily species-rich clade of beetles known as the Phytophaga. ATUMI thus provides a contrasting view as a neighbor for one of the most successful known animal groups.

Results

We present a robust genome assembly and a gene set possessing 97.5% of the core proteins known from the holometabolous insects. The ATUMI genome encodes fewer enzymes for plant digestion than the genomes of wood-feeding beetles but nonetheless shows signs of broad metabolic plasticity. Gustatory receptors are few in number compared to other beetles, especially receptors with known sensitivity (in other beetles) to bitter substances. In contrast, several gene families implicated in detoxification of insecticides and adaptation to diverse dietary resources show increased copy numbers. The presence and diversity of homologs involved in detoxification differ substantially from the bee hosts of ATUMI.

Conclusions

Our results provide new insights into the genomic basis for local adaption and invasiveness in ATUMI and a blueprint for control strategies that target this pest without harming their honey bee hosts. A minimal set of gustatory receptors is consistent with the observation that, once a host colony is invaded, food resources are predictable. Unique detoxification pathways and pathway members can help identify which treatments might control this species even in the presence of honey bees, which are notoriously sensitive to pesticides.