Unique features of a global human ectoparasite identified through sequencing of the bed bug genome

Under the radar: Transcriptomic responses of bed bugs to an entomopathogen, environmental bacteria, and a human pathogen

Bed bugs (Hemiptera: Cimicidae) are widely distributed, obligately blood-feeding insects, but they have never been linked to pathogen transmission in humans. Most other hematophagous insects that frequently bite humans transmit pathogens, and it is unclear why bed bugs do not. One hypothesis is that bed bugs have evolved a highly robust immune system because their mating system, traumatic insemination, exposes females to consistent wounding and bacterial infections. Although this has been proposed, very little is known about the bed bug immune system and how bed bugs respond to microbial challenges introduced by wounding. Similarly, there is little known about how the bed bug immune system responds to human pathogens. Understanding the bed bug immune system could give insight to why bed bugs appear not to transmit disease and under what circumstances they could, while also facilitating biological control efforts involving microbes. To investigate the transcriptomic response of bed bugs to immune challenges, we exposed female bed bugs to three bacterial challenges. 1.) Pseudomonas fluorescens, an entomopathogen known to have harmful effects to bed bugs, 2.) bacteria cultured from a bed bug enclosure (99.9 % Bacillus spp.), likely encountered during traumatic insemination, and 3.) Borrelia duttoni, a human vector-borne pathogen that causes relapsing fever. We compared the transcriptomes of infected bed bugs with uninfected matched controls in a pairwise fashion, focusing on immune-related genes. We found many known antimicrobial effector genes upregulated in response to P. fluorescens and traumatic insemination-associated bacteria, but interestingly, not in response to B. duttoni. In the differentially expressed genes that were shared between experiments, we found significant overlap in the P. fluorescens treatment and the traumatic insemination bacteria treatment, and between the P. fluorescens and B. duttoni treatments, but not between the traumatic insemination bacteria treatment and the B. duttoni treatment. Finally, we identify previously overlooked candidates for future studies of immune function in bed bugs, including a peroxidase-like gene, many putative cuticle-associated genes, a laccase-like gene, and a mucin-like gene. By taking a comprehensive transcriptomic approach, our study is an important step in understanding how bed bugs respond to diverse immune challenges.

Semiochemicals modulating bed bug behaviour

The common bed bug, Cimex lectularius (Cimicidae, Hemiptera), is obligate haematophagous and can carry pathogens but is not known to transmit diseases to humans in natural settings. Bed bugs are activated and guided by semiochemicals and heat by potential hosts, as well as by an aggregation pheromone when searching for food and home. No sexual pheromone has been detected in bed bugs, which mate through traumatic insemination. First, after mounting, males are able to distinguish between large nymphs, other males, and females. To avoid unwanted traumatic inseminations, nymphs and adults emit an alarm pheromone. The olfactory system of bed bugs has relatively few odorant and ionotropic receptors, which likely reflects the simple environment that bed bugs live in, and use to search for hosts and conspecifics.

Species-specific duplicated FMRFaR-like gene A62 regulates spontaneous locomotion in Apolygus lucorum

BACKGROUND

Apolygus lucorum, a major cotton pest, has undergone a significant expansion of the FMRFaR gene within the GPCR superfamily, resulting in two classes of GPCR, namely FMRFaR (A54-55) and newly duplicated FMRFaR-like (A56-62). Notably, FMRFaR-like genes, particularly A62, show enhanced expression in the legs and wings of adults, indicating their potential role in locomotion. Employing A62 as a representative of FMRFaR-like, our study investigates the influence of FMRFa, FMRFaR, and FMRFaR-like on locomotion and development of A. lucorum.

RESULTS

FMRFaR and FMRFa exhibit comparable temporal and tissue expression patterns, whereas the FMRFaR-like genes within A. lucorum exhibit completely distinct evolutionary and expression patterns compared to classical FMRFaR. RNA interference (RNAi) experiments revealed that suppressing FMRFa expression results in complete lethality in A. lucorum, but neither FMRFaR nor A62 exhibit the same effect after RNAi. Suppressing the expression of FMRFa only decreases the expression of the A54 gene simultaneously, suggesting that A54 may function as a classical FMRFaR activated by FMRFa. RNAi of A62 leads to wing malformation and a significant reduction in spontaneous movement behavior in A. lucorum. Further transcriptomic analysis revealed that A62 affects the A. lucorum's movement behavior through energy metabolism pathways and motor protein pathways.

CONCLUSION

Our study unveils the unique and complex roles of FMRFa and its receptor in A. lucorum. These findings provide valuable insights into potential targets for pest control strategies aimed at managing A. lucorum populations in cotton fields. © 2024 Society of Chemical Industry.

Deployment and transcriptional evaluation of nitisinone, an FDA-approved drug, to control bed bugs

Bed bugs are blood-feeders that rapidly proliferate into large indoor infestations. Their bites can cause allergies, secondary infections and psychological stress, among other problems. Although several tactics for their management have been used, bed bugs continue to spread worldwide wherever humans reside. This is mainly due to human-mediated transport and their high resistance to several classes of insecticides. New treatment options with novel modes of action are required for their control. In this study, we evaluated the use of nitisinone (NTBC), an FDA-approved drug, for bed bug control in an insecticide-susceptible (HH) and an insecticide-resistant (CIN) population. Although NTBC was lethal to both populations when administered orally or applied topically in very low doses, we observed a slight but significant resistance in the CIN population. Transcriptomic analysis in both populations indicated that NTBC treatment elicited a broad suppression of genes associated with RNA post-transcriptional modifications, translation, endomembrane system, protein post-translational modifications and protein folding. The CIN population exhibited higher ATP production and xenobiotic detoxification. Feeding studies on a mouse model highlight that NTBC could be used as a control method of bed bugs by host treatment. The results demonstrate that NTBC can be used as a new active ingredient for bed bug control by topical or oral treatment and shed light on the molecular mechanisms of suppressed tyrosine metabolism following NTBC treatment.

A novel prolixicin identified in common bed bugs with activity against both bacteria and parasites

The hematophagous common bed bug, Cimex lectularius, is not known to transmit human pathogens outside laboratory settings, having evolved various immune defense mechanisms including the expression of antimicrobial peptides (AMPs). We unveil three novel prolixicin AMPs in bed bugs, exhibiting strong homology to the prolixicin of kissing bugs, Rhodnius prolixus, and to diptericin/attacin AMPs. We demonstrate for the first time sex-specific and immune mode-specific upregulation of these prolixicins in immune organs, the midgut and rest of body, following injection and ingestion of Gr+ (Bacillus subtilis) and Gr- (Escherichia coli) bacteria. Synthetic CL-prolixicin2 significantly inhibited growth of E. coli strains and killed or impeded Trypanosoma cruzi, the Chagas disease agent. Our findings suggest that prolixicins are regulated by both IMD and Toll immune pathways, supporting cross-talk and blurred functional differentiation between major immune pathways. The efficacy of CL-prolixicin2 against T. cruzi underscores the potential of AMPs in Chagas disease management.

Genome of tropical bed bug Cimex hemipterus (Cimicidae, Hemiptera) reveals tetraspanin expanded in bed bug ancestor

Cimex species are ectoparasites that exclusively feed on warm-blooded animals such as birds and mammals. Three cimicid species are known to be persistent pests for humans, including the tropical bed bug Cimex hemipterus, common bed bug Cimex lectularius, and Eastern bat bug Leptocimex boueti. To date, genomic information is restricted to the common bed bug C. lectularius, which limits understanding their biology and to provide controls of bed bug infestations. Here, a chromosomal-level genome assembly of C. hemipterus (495 Mb [megabase pairs]) contained on 16 pseudochromosomes (scaffold N50 = 34 Mb), together with 9 messenger RNA and small RNA transcriptomes were obtained. In comparison between hemipteran genomes, we found that the tetraspanin superfamily was expanded in the Cimex ancestor. This study provides the first genome assembly for the tropical bed bug C. hemipterus, and offers an unprecedented opportunity to address questions relating to bed bug infestations, as well as genomic evolution to hemipterans more widely.

Dynamics of cuticle-associated transcript profiles during moulting of the bed bug Cimexlectularius

The bed bug Cimex lectularius is a worldwide human pest. The sequenced genome allows molecular analyses of all aspects of bed bug biology. The present work was conducted to contribute to bed bug cuticle biology. As in other insect species, the C. lectularius cuticle consists of the three horizontal layers procuticle, epicuticle and envelope. To analyse the genes needed for the establishment of the stratified cuticle, we studied the expression pattern of 42 key cuticle-related genes at the transition of the penultimate nymphal stage to adult animals when a new cuticle is formed. Based on gene expression dynamics, in simplified model, we distinguish two key events during cuticle renewal in C. lectularius. First, upon blood feeding, modulation of ecdysone signalling culminates in the transcriptional activation of the transcription factor Clec-Ftz-F1 that possibly controls the expression of 32 of the 42 genes tested. Second, timed expression of Clec-Ftz-F1 seems to depend also on the insulin signalling pathway as RNA interference against transcripts of the insulin receptor delays Clec-Ftz-F1 expression and stage transition. An important observation of our transcript survey is that genes needed for the construction of the three cuticle layers are largely expressed simultaneously. Based on these data, we hypothesise a considerable synchronous mechanism of layer formation rather than a strictly sequential one. Together, this work provides a basis for functional analyses of cuticle formation in C. lectularius.

The puncture mechanics: an example from the bed bug Cimex lectularius showing traumatic insemination using the paramere

Cimicidae are well-known for traumatic insemination, and males pierce females with their parameres and transfer sperm through them. The shape of parameres is relatively stable in the family, but in some genera, the paramere is elongated, appearing less resistant against lateral deflection. To understand the mechanical limitations of the paramere, we studied its penetration mechanics of the common bed bug, Cimex lectularius. We examined the post-abdominal morphology, paramere geometry and material properties and conducted breaking stress experiments on the paramere under wet and dry conditions. Mechanical property gradients are present with the paramere tip as the stiffest region and the base as the most flexible one. These mechanical properties relate to the presence of Ca, Zn and Si. The basal wing-shaped structure is flexible, enabling it to interlock with the anal region during mating. The paramere is slightly twisted; the tip region is circular in cross-section, and the geometry of the rest is rather complex. In the mechanical tests, wet parameres mainly buckled, while dried parameres broke off. The level of structural failures depended on directions from which the compression forces were applied. Structural, material and mechanical strengthening mechanisms preventing the paramere from mechanical failure are discussed.

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.

Chromosome-level genome of the poultry shaft louse Menopon gallinae provides insight into the host-switching and adaptive evolution of parasitic lice

BACKGROUND

Lice (Psocodea: Phthiraptera) are one important group of parasites that infects birds and mammals. It is believed that the ancestor of parasitic lice originated on the ancient avian host, and ancient mammals acquired these parasites via host-switching from birds. Here we present the first chromosome-level genome of Menopon gallinae in Amblycera (earliest diverging lineage of parasitic lice). We explore the transition of louse host-switching from birds to mammals at the genomic level by identifying numerous idiosyncratic genomic variations.

RESULTS

The assembled genome is 155 Mb in length, with a contig N50 of 27.42 Mb. Hi-C scaffolding assigned 97% of the bases to 5 chromosomes. The genome of M. gallinae retains a basal insect repertoire of 11,950 protein-coding genes. By comparing the genomes of lice to those of multiple representative insects in other orders, we discovered that gene families of digestion, detoxification, and immunity-related are generally conserved between bird lice and mammal lice, while mammal lice have undergone a significant reduction in genes related to chemosensory systems and temperature. This suggests that mammal lice have lost some of these genes through the adaption to environment and temperatures after host-switching. Furthermore, 7 genes related to hematophagy were positively selected in mammal lice, suggesting their involvement in the hematophagous behavior.

CONCLUSIONS

Our high-quality genome of M. gallinae provides a valuable resource for comparative genomic research in Phthiraptera and facilitates further studies on adaptive evolution of host-switching within parasitic lice.

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.

Characterization of the common bed bug's eggshell and egg glue proteins

An insect egg is one of the most vulnerable stages of insect life, and the evolutionary success of a species depends on the eggshell protecting the embryo and the egg glue securing the attachment. The common bed bug (Cimex lectularius), notorious for its painful and itchy bites, infests human dwellings to feed on blood. They are easier to find these days as they adapt to develop resistance against commonly used insecticides. In this study, we identify and characterize the eggshell protein and the probable egg glue protein (i.e. keratin associated protein 5-10 like protein) of the bed bug by using mass spectrometry and bioinformatics analysis. Furthermore, by using transcription profiling and in vivo RNA interference, we show evidences that the keratin associated protein 5-10 like protein functions as the glue protein. Finally, structural characterizations on the two proteins are performed using recombinant proteins. Amino acid sequences of various insect eggshell and egg glue proteins support their independent evolution among different insect groups. Hence, inhibiting the function of these proteins related to the earliest stage of life can achieve species-specific population control. In this respect, our results would be a starting point in developing new ways to control bed bug population.

Transcriptome analysis reveals salivary gland-specific neuropeptide signaling genes in the predatory stink bug, Picromerus lewisi

Predatory stink bugs derive from phytophagous stink bugs and evolved enhanced predation skills. Neuropeptides are a diverse class of ancient signaling molecules that regulate physiological processes and behavior in animals, including stink bugs. Neuropeptide evolution might be important for the development of predation because neuropeptides can be converted to venoms that impact prey. However, information on neuropeptide signaling genes in predatory stink bugs is lacking. In the present study, neuropeptide signaling genes of Picromerus lewisi, an important predatory stink bug and an effective biological agent, were comprehensively identified by transcriptome analysis, with a total of 59 neuropeptide precursor genes and 58 potential neuropeptide receptor genes found. In addition, several neuropeptides and their receptors enriched in salivary glands of P. lewisi were identified. The present study and subsequent functional research contribute to an in-depth understanding of the biology and behavior of the predatory bugs and can provide basic information for the development of better pest management strategies, possibly including neuropeptide receptors as insecticide targets and salivary gland derived venom toxins as novel killing moleculars.

Parallel Losses of Blue Opsin Correlate with Compensatory Neofunctionalization of UV-Opsin Gene Duplicates in Aphids and Planthoppers

Expanding on previous efforts to survey the visual opsin repertoires of the Hemiptera, this study confirms that homologs of the UV- and LW-opsin subfamilies are conserved in all Hemiptera, while the B-opsin subfamily is missing from the Heteroptera and subgroups of the Sternorrhyncha and Auchenorrhyncha, i.e., aphids (Aphidoidea) and planthoppers (Fulgoroidea), respectively. Unlike in the Heteroptera, which are characterized by multiple independent expansions of the LW-opsin subfamily, the lack of B-opsin correlates with the presence of tandem-duplicated UV-opsins in aphids and planthoppers. Available data on organismal wavelength sensitivities and retinal gene expression patterns lead to the conclusion that, in both groups, one UV-opsin paralog shifted from ancestral UV peak sensitivity to derived blue sensitivity, likely compensating for the lost B-opsin. Two parallel bona fide tuning site substitutions compare to 18 non-corresponding amino acid replacements in the blue-shifted UV-opsin paralogs of aphids and planthoppers. Most notably, while the aphid blue-shifted UV-opsin clade is characterized by a replacement substitution at one of the best-documented UV/blue tuning sites (Rhodopsin site 90), the planthopper blue-shifted UV-opsin paralogs retained the ancestral lysine at this position. Combined, the new findings identify aphid and planthopper UV-opsins as a new valuable data sample for studying adaptive opsin evolution.

Sleep: An Essential and Understudied Process in the Biology of Blood-Feeding Arthropods

Understanding the biology of blood-feeding arthropods is critical to managing them as vectors of etiological agents. Circadian rhythms act in the regulation of behavioral and physiological aspects such as blood feeding, immunity, and reproduction. However, the impact of sleep on these processes has been largely ignored in blood-feeding arthropods, but recent studies in mosquitoes show that sleep-like states directly impact host landing and blood feeding. Our focus in this review is on discussing the relationship between sleep and circadian rhythms in blood-feeding arthropods along with how unique aspects such as blood gluttony and dormancy can impact sleep-like states. We highlight that sleep-like states are likely to have profound impacts on vector-host interactions but will vary between lineages even though few direct studies have been conducted. A myriad of factors, such as artificial light, could directly impact the time and levels of sleep in blood-feeding arthropods and their roles as vectors. Lastly, we discuss underlying factors that make sleep studies in blood-feeding arthropods difficult and how these can be bypassed. As sleep is a critical factor in the fitness of animal systems, a lack of focus on sleep in blood-feeding arthropods represents a significant oversight in understanding their behavior and its role in pathogen transmission.

Lipid metabolism dysfunction following symbiont elimination is linked to altered Kennedy pathway homeostasis

Lipid metabolism is critical for insect reproduction, especially for species that invest heavily in the early developmental stages of their offspring. The role of symbiotic bacteria during this process is understudied but likely essential. We examined the role of lipid metabolism during the interaction between the viviparous tsetse fly (Glossina morsitans morsitans) and its obligate endosymbiotic bacteria (Wigglesworthia glossinidia) during tsetse pregnancy. We observed increased CTP:phosphocholine cytidylyltransferase (cct1) expression during pregnancy, which is critical for phosphatidylcholine biosynthesis in the Kennedy pathway. Experimental removal of Wigglesworthia impaired lipid metabolism via disruption of the Kennedy pathway, yielding obese mothers whose developing progeny starve. Functional validation via experimental cct1 suppression revealed a phenotype similar to females lacking obligate Wigglesworthia symbionts. These results indicate that, in Glossina, symbiont-derived factors, likely B vitamins, are critical for the proper function of both lipid biosynthesis and lipolysis to maintain tsetse fly fecundity.

Metabolic resistance to deltamethrin is mediated by P450 and esterases in common bed bugs Cimex lectularius L. (Heteroptera: Cimicidae)

The infestations of Cimex lectularius L. (Heteroptera: Cimicidae) registered in the last decades have been influenced by several human activities, including international tourism and commerce. Moreover, the development of insecticide resistance and careless pest control strategies contributed to the dispersal of bed bugs. Given the complexity of the topic, distinguishing physiological and molecular mechanisms involved in resistance can help design proper control tools and limit the resistance spread. Here we determined the susceptibility to deltamethrin and imidacloprid in bed bugs collected in Italy. Also, we assessed the role of esterases and P450 monooxygenases by direct enzymatic activity measurement and inhibition by synergism bioassays. Our results showed that the field-collected colony exhibited high resistance ratios to imidacloprid and deltamethrin (757 and >60,000 times, respectively) compared to the susceptible colony. Moreover, resistant bed bugs showed increased activity of esterases and P450 monooxygenases. The synergistic effect of piperonyl butoxide (PBO) suggests the significant contribution of both enzymatic groups as detoxification pathways implicated in pyrethroid-resistant bed bugs. Further investigations are needed to unravel the biochemical and molecular basis involved in the resistant phenotype for developing novel strategies for pest control.

The antennal transcriptome of Triatoma infestans reveals substantial expression changes triggered by a blood meal

BACKGROUND

Triatoma infestans is the main vector of Chagas disease in the Americas, currently transmitting it in Argentina, Paraguay, and Bolivia. Many T. infestans populations present insecticide resistance, reducing the efficiency of control campaigns. Alternative vector control methods are needed, and molecular targets mediating fundamental physiological processes can be a promising option to manipulate kissing bug behavior. Therefore, it is necessary to characterize the main sensory targets, as well as to determine whether they are modulated by physiological factors. In order to identify gene candidates potentially mediating host cue detection, the antennal transcripts of T. infestans fifth instar larvae were sequenced and assembled. Besides, we evaluated whether a blood meal had an effect on transcriptional profiles, as responsiveness to host-emitted sensory cues depends on bug starvation.

RESULTS

The sensory-related gene families of T. infestans were annotated (127 odorant receptors, 38 ionotropic receptors, 11 gustatory receptors, 41 odorant binding proteins, and 25 chemosensory proteins, among others) and compared to those of several other hemipterans, including four triatomine species. Several triatomine-specific lineages representing sensory adaptations developed through the evolution of these blood-feeding heteropterans were identified. As well, we report here various conserved sensory gene orthogroups shared by heteropterans. The absence of the thermosensor pyrexia, of pickpocket receptor subfamilies IV and VII, together with clearly expanded takeout repertoires, are revealed features of the molecular bases of heteropteran antennal physiology. Finally, out of 2,122 genes whose antennal expression was significantly altered by the ingestion of a blood meal, a set of 41 T. infestans sensory-related genes (9 up-regulated; 32 down-regulated) was detected.

CONCLUSIONS

We propose that the set of genes presenting nutritionally-triggered modulation on their expression represent candidates to mediate triatomine host-seeking behavior. Besides, the triatomine-specific gene lineages found represent molecular adaptations to their risky natural history that involves stealing blood from an enormously diverse set of vertebrates. Heteropteran gene orthogroups identified may represent unknown features of the sensory specificities of this largest group of hemipteroids. Our work is the first molecular characterization of the peripheral modulation of sensory processes in a non-dipteran vector of human disease.

Comparative analysis of detoxification-related gene superfamilies across five hemipteran species

Background

Hemiptera is one of the most speciose orders of insects, and the most speciose considering Hemimetabola. Through their evolutive history, hemipterans with different feeding habits have adapted to deal with different chemical challenges. Three major gene families are involved in xenobiotic detoxification in insects: the cytochromes P450 (CYPs), carboxyl/cholinesterases (CCEs), and glutathione transferases (GSTs). Here we perform a comparative analysis on the complement of these gene superfamilies across five hemipteran species; four heteropterans (the pentatomid plant feeders Nezara viridula and Halyomorpha halys; the hematophagous Cimex lectularius, Cimicidae, and Rhodnius prolixus, Reduviidae), and one Auchenorrhyncha plant feeder (Nilaparvata lugens).

Results

Our results point to an expansion of several enzyme families associated with xenobiotic detoxification in heteropterans with respect to other species and the existence of a dynamic evolution pattern including CYP3 clan, hormone and pheromone processing class in the CCE superfamily, and sigma class in GST superfamily. Other detoxification-related families are reduced in the hemipteran species analyzed here: reduction or even absence of epsilon class and reduced delta class in GST superfamily; absence of mitochondrial CYP12 family; absence of CYP9 family in CYP3 clan; and reduction or even absence of some dietary/detoxification groups of CCEs. Interestingly, the most polyphagous species analyzed here (H. halys) is also the one that presents the largest repertoire of detoxification enzymes. Gene cluster analysis suggests that this could be due to gene duplication events.

Conclusions

The evolutionary analysis performed here reveals characteristics that are both common and particular for heteropterans. The composition and organization of detoxification-related gene families could shed light on evolutionary forces that shaped their divergence. These families are important for both the detoxification of diet products and for conferring tolerance or resistance to synthetic insecticides. Furthermore, we present the first comprehensive analysis of detoxification gene superfamilies in N. viridula, an understudied species in spite of its economic relevance as a crop pest. The information obtained is of interest for basic insect science as well as for the control of harmful species and the management of insecticide resistance.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08974-y.

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.

Characterization of New Defensin Antimicrobial Peptides and Their Expression in Bed Bugs in Response to Bacterial Ingestion and Injection

Common bed bugs, Cimex lectularius, can carry, but do not transmit, pathogens to the vertebrate hosts on which they feed. Some components of the innate immune system of bed bugs, such as antimicrobial peptides (AMPs), eliminate the pathogens. Here, we determined the molecular characteristics, structural properties, and phylogenetic relatedness of two new defensins (CL-defensin1 (XP_024085718.1), CL-defensin2 (XP_014240919.1)), and two new defensin isoforms (CL-defensin3a (XP_014240918.1), CL-defensin3b (XP_024083729.1)). The complete amino acid sequences of CL-defensin1, CL-defensin2, CL-defensin3a, and CL-defensin3b are strongly conserved, with only minor differences in their signal and pro-peptide regions. We used a combination of comparative transcriptomics and real-time quantitative PCR to evaluate the expression of these defensins in the midguts and the rest of the body of insects that had been injected with bacteria or had ingested blood containing the Gram-positive (Gr+) bacterium Bacillus subtilis and the Gram-negative (Gr–) bacterium Escherichia coli. We demonstrate, for the first time, sex-specific and immunization mode-specific upregulation of bed bug defensins in response to injection or ingestion of Gr+ or Gr– bacteria. Understanding the components, such as these defensins, of the bed bugs’ innate immune systems in response to pathogens may help unravel why bed bugs do not transmit pathogens to vertebrates.

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

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

Comparative Efficacy of a Fungal Entomopathogen with a Broad Host Range against Two Human-Associated Pests

The ability of a fungal entomopathogen to infect an insect depends on a variety of factors, including strain, host, and environmental conditions. Similarly, an insect’s ability to prevent fungal infection is dependent on its biology, environment, and evolutionary history. Synanthropic pests have adapted to thrive in the indoor environment, yet they arose from divergent evolutionary lineages and occupy different feeding guilds. The hematophagous bed bug (Cimex lectularius) and omnivorous German cockroach (Blattella germanica) are highly successful indoors, but have evolved different physiological and behavioral adaptations to cope with the human-built environment, some of which also reduce the efficacy of fungal biopesticides. In order to gain greater insight into the host barriers that prevent or constrain fungal infection in bed bugs and German cockroaches, we tested different doses of Beauveria bassiana GHA through surface contact, topical application, feeding, and injection. Bed bugs were generally more susceptible to infection by B. bassiana with the mode of delivery having a significant impact on infectivity. The German cockroach was highly resilient to infection, requiring high doses of fungal conidia (>8.8 × 104) delivered by injection into the hemocoel to cause mortality. Mortality occurred much faster in both insect species after exposure to surfaces dusted with dry conidia than surfaces treated with conidia suspended in water or oil. These findings highlight the importance of developing innovative delivery techniques to enhance fungal entomopathogens against bed bugs and cockroaches.

The Unfolded-Protein Response Triggers the Arthropod Immune Deficiency Pathway

The insect immune deficiency (IMD) pathway is a defense mechanism that senses and responds to Gram-negative bacteria. Ticks lack genes encoding upstream components that initiate the IMD pathway. Despite this deficiency, core signaling molecules are present and functionally restrict tick-borne pathogens. The molecular events preceding activation remain undefined. Here, we show that the unfolded-protein response (UPR) initiates the IMD network. The endoplasmic reticulum (ER) stress receptor IRE1α is phosphorylated in response to tick-borne bacteria but does not splice the mRNA encoding XBP1. Instead, through protein modeling and reciprocal pulldowns, we show that Ixodes IRE1α complexes with TRAF2. Disrupting IRE1α-TRAF2 signaling blocks IMD pathway activation and diminishes the production of reactive oxygen species. Through in vitro, in vivo, and ex vivo techniques, we demonstrate that the UPR-IMD pathway circuitry limits the Lyme disease-causing spirochete Borrelia burgdorferi and the rickettsial agents Anaplasma phagocytophilum and A. marginale (anaplasmosis). Altogether, our study uncovers a novel linkage between the UPR and the IMD pathway in arthropods. IMPORTANCE The ability of an arthropod to harbor and transmit pathogens is termed "vector competency." Many factors influence vector competency, including how arthropod immune processes respond to the microbe. Divergences in innate immunity between arthropods are increasingly being reported. For instance, although ticks lack genes encoding key upstream molecules of the immune deficiency (IMD) pathway, it is still functional and restricts causative agents of Lyme disease (Borrelia burgdorferi) and anaplasmosis (Anaplasma phagocytophilum). How the IMD pathway is activated in ticks without classically defined pathway initiators is not known. Here, we found that a cellular stress response network, the unfolded-protein response (UPR), functions upstream to induce the IMD pathway and restrict transmissible pathogens. Collectively, this explains how the IMD pathway can be activated in the absence of canonical pathway initiators. Given that the UPR is highly conserved, UPR-initiated immunity may be a fundamental principle impacting vector competency across arthropods.

Allatostatin A Signalling: Progress and New Challenges From a Paradigmatic Pleiotropic Invertebrate Neuropeptide Family

Neuropeptides have gained broad attraction in insect neuroscience and physiology, as new genetic tools are increasingly uncovering their wide-ranging pleiotropic functions with high cellular resolution. Allatostatin A (AstA) peptides constitute one of the best studied insect neuropeptide families. In insects and other panarthropods, AstA peptides qualify as brain-gut peptides and have regained attention with the discovery of their role in regulating feeding, growth, activity/sleep and learning. AstA receptor homologs are found throughout the protostomia and group with vertebrate somatostatin/galanin/kisspeptin receptors. In this review, we summarise the current knowledge on the evolution and the pleiotropic and cell-specific non-allatostatic functions of AstA. We speculate about the core functions of AstA signalling, and derive open questions and challengesfor future research on AstA and invertebrate neuropeptides in general.

The Adipokinetic Peptides of Hemiptera: Structure, Function, and Evolutionary Trends

The Hemiptera comprise the most species-rich order of the hemimetabolous insects. Members of a number of superfamilies, most notably especially the more basal ones such as white flies, psyllids and aphids, belong to the most destructive agricultural insects known worldwide. At the other end of the phylogenetic tree are hemipterans that are notorious medical pests (e.g. kissing bugs). Most of the hemipteran species are good flyers, and lipid oxidation plays a pivotal role to power the contraction of flight muscles and, in aquatic water bugs, also deliver the ATP for the extensive swimming action of the leg muscles. Mobilization of stored lipids (mostly triacylglycerols in the fat body) to circulating diacylglycerols in the hemolymph is regulated by a set of small neuropeptides, the adipokinetic hormones (AKHs). We searched the literature and publicly available databases of transcriptomes and genomes to present here AKH sequences from 191 hemipteran species. Only few of these peptides were sequenced via Edman degradation or mass spectrometry, and even fewer were characterized with molecular biology methods; thus, the majority of the AKHs we have identified by bioinformatics are merely predicted sequences at this stage. Nonetheless, a total of 42 AKH primary sequences are assigned to Hemiptera. About 50% of these structures occur also in other insect orders, while the remaining 50% are currently unique for Hemiptera. We find 9 novel AKHs not shown to be synthesized before in any insect. Most of the hemipteran AKHs are octapeptides (28) but there is an impressive number of decapeptides (12) compared to other speciose orders such as Diptera and Lepidoptera. We attempt to construct a hypothetical molecular peptide evolution of hemipteran AKHs and find quite a bit of overlapping with current phylogenetic ideas of the Hemiptera. Lastly, we discuss the possibility to use the sequence of the aphid AKH as lead peptide for the research into a peptide mimetic fulfilling criteria of a green insecticide.

Analysis of Transcriptome Difference between Blood-Fed and Starved Tropical Bed Bug, Cimex hemipterus (F.) (Hemiptera: Cimicidae)

Simple Summary

Bed bugs are well known for their extreme resilience to starvation. The molecular mechanisms behind this ability, however, are little known. Thus, the whole transcriptomes of blood-fed and starved bed bugs from the species Cimex hemipterus (tropical bed bugs) were sequenced and compared. The transcriptome of tropical bed bugs was initially annotated. Following differentially expressed genes (DEGs) analysis, regulated transcripts were mostly identified in biological processes during blood-feeding and starvation. The results provide an overview of the functional genes proportion of this species and a deeper understanding of the bed bug’s molecular mechanism of resistance to blood feeding and starvation.

Abstract

The reference transcriptome for Cimex hemipterus (tropical bed bug) was assembled de novo in this study, and differential expression analysis was conducted between blood-fed and starved tropical bed bug. A total of 24,609 transcripts were assembled, with around 79% of them being annotated against the Eukaryotic Orthologous Groups (KOG) database. The transcriptomic comparison revealed several differentially expressed genes between blood-fed and starved bed bugs, with 38 of them being identifiable. There were 20 and 18 genes significantly upregulated in blood-fed and starved bed bugs, respectively. Differentially expressed genes (DEGs) were revealed to be associated with regulation, metabolism, transport, motility, immune, and stress response; endocytosis; and signal transduction. The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis showed an enrichment of genes encoding steroid biosynthesis, glycosaminoglycan biosynthesis, butanoate metabolism, and autophagy in both blood-fed and starved bed bugs. However, in blood-fed bed bugs, genes involved in histidine metabolism, caffeine metabolism, ubiquinone/terpenoid-quinone biosynthesis, and sulfur relay system were enriched. On the other hand, starvation activates genes related to nicotinate and nicotinamide metabolism, fatty acid elongation, terpenoid backbone biosynthesis, metabolism of xenobiotics by cytochrome P450, riboflavin metabolism, apoptosis, and protein export. The present study is the first to report a de novo transcriptomic analysis in C. hemipterus and demonstrated differential responses of bed bugs in facing blood-feeding and starvation.

Molecular (co)evolution of hymenopteran social parasites and their hosts

Social parasitism describes a fascinating way of life in which species exploit the altruistic behaviour of closely related, social species. Social parasites have repeatedly evolved in the social Hymenoptera, including ants, bees, and wasps. The common ancestry and shared (social) environment with their hosts facilitates the study of molecular adaptations to the parasitic lifestyle. Moreover, when social parasites are widespread and virulent, they exert strong selection pressure on their hosts, leading to the evolution of defense mechanisms and triggering a coevolutionary arms race. Recent advances in sequencing technology now make it possible to study the molecular basis of this coevolutionary process. In addition to describing the latest developments, we highlight open research questions that could be tackled with genomic, transcriptomic, or epigenetic data.

De novo transcriptome sequencing of the northern fowl mite, Ornithonyssus sylviarum, shed light on parasitiform poultry mites evolution and its chemoreceptor repertoires

The northern fowl mite (NFM), Ornithonyssus sylviarum, and the poultry red mite (PRM), Dermanyssus gallinae, are the most serious pests of poultry, both of which have an expanding global prevalence. Research on NFM has been constrained by a lack of genomic and transcriptomic data. Here, we report and analyze the first global transcriptome data across all mite live stages and sexes. A total of 28,999 unigenes were assembled, of which 19,750 (68.10%) were annotated using seven functional databases. The biological function of these unigenes was classified using the GO, KOG, and KEGG databases. To gain insight into the chemosensory receptor-based system of parasitiform mites, we furthermore assessed the gene repertoire of gustatory receptors (GRs) and ionotropic receptors (IRs), both of which encode putative ligand-gated ion channel proteins. While these receptors are well characterized in insect model species, our understanding of chemosensory detection in mites and ticks is in its infancy. To address this paucity of data, we identified 9 IR/iGluRs and 2 GRs genes by analyzing transcriptome data in the NFM, while 9 GRs and 41 IR/iGluRs genes were annotated in the PRM genome. Taken together, the transcriptomic and genomic annotation of these two species provide a valuable reference for studies of parasitiform mites and also help to understand how chemosensory gene family expansion/contraction events may have been reshaped by an obligate parasitic lifestyle compared with their free-living closest relatives. Future studies should include additional species to validate this observation and functional characterization of the identified proteins as a step forward in identifying tools for controlling these poultry pests.

The closer the better: Sensory tools and host-association in blood-sucking insects

Many hematophagous insects acquire medical and veterinary relevance because they transmit disease causing pathogens to humans. Hematophagy is only fulfilled once a blood feeder successfully locates a vertebrate host by means of fine sensory systems. In nature, blood-sucking insects can exploit environments with differential association with their hosts. Given the relevance of the sensory systems during host searching, we review the current state of knowledge of the sensory machinery of four blood-sucking insects: human lice, bed bugs, kissing bugs and mosquitoes. Each one is representative of highly anthropophilic behaviours and a different degree of association with human hosts. We compare the number, arrangement and functional type of cuticular sensory structures dispersed on the main sensory organs. We also compare the genetic machinery potentially involved in the detection of host stimuli. Finally, we discuss the sensory diversity of the insects studied here.

Characterization of heat exposure-associated escape behaviors and HSP gene expression in bed bugs (Cimex lectularius L.)

BACKGROUND

Heat can be effective for bed bug elimination. However, in some cases bed bugs survive heat treatments. The objectives of this study were to determine the behavioral responses of bed bugs to rising harborage temperatures (23.0-49.0 °C) and identify which heat shock protein (HSP) genes are expressed after heat exposure. First, a custom-made copper arena and harborage were used to determine the escape behaviors of six bed bug populations. Next, HSP gene expression responses of select populations were determined after heat exposure using real time quantitative polymerase chain reaction (RT-qPCR).

RESULTS

Analysis of the 25 min behavioral experiment data found that harborage top temperatures associated with 25%, 50% and 75% probabilities of bed bugs to flee the harborage did not differ significantly between populations. Also, the percentage of insects that escaped from heated areas and survived (4.0-12.0%) was not different between populations. However, when specific temperatures at which successful escapes occurred were statistically compared, the Poultry House population was found to flee the harborage at statistically higher temperatures (43.6 ± 0.5 °C) than others (40.5 ± 0.6-42.0 ± 0.7 °C). The RT-qPCR experiments revealed that the HSP70.1, HSP70.3, and Putative Small HSP genes were significantly up-regulated 15 min, 2, and 4 h post-heat exposure and decreased back to baseline levels by 24 h.

CONCLUSIONS

This study shows that when harborage top temperatures approach 40.0-43.0 °C, bed bugs will disperse in search for cooler areas. This work implicates the HSP70.1, HSP70.3, and Putative Small HSP genes in heat induced stress recovery of bed bugs. © 2021 Society of Chemical Industry.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Next-Generation Sequencing and Its Impacts on Entomological Research in Ecology and Evolution

The advent of NGS-based methods has been profoundly transforming entomological research. Through continual development and improvement of different methods and sequencing platforms, NGS has promoted mass elucidation of partial or whole genetic materials associated with beneficial insects, pests (of agriculture, forestry and animal, and human health), and species of conservation concern, helping to unravel ecological and evolutionary mechanisms and characterizing survival, trophic interactions, and dispersal. It is shifting the scale of biodiversity and environmental analyses from individuals and biodiversity indicator species to the large-scale study of communities and ecosystems using bulk samples of species or a mixed "soup" of environmental DNA. As the NGS-based methods have become more affordable, complexity demystified, and specificity and sensitivity proven, their use in entomological research has spread widely. This article presents several examples on how NGS-based methods have been used in entomology to provide incentives to apply them when appropriate and to open our minds to the expected advances in entomology that are yet to come.

Chromosome-level genome assembly of the mirid predator Cyrtorhinus lividipennis Reuter (Hemiptera: Miridae), an important natural enemy in the rice ecosystem

Though the genomes of many rice herbivorous pests have recently been well characterized, little is known about the genome of their natural enemies. Here, by using the Illumina and PacBio platforms, we sequenced and assembled the whole genome of the mirid species Cyrtorhinus lividipennis Reuter (Hemiptera: Miridae), which is an economically and ecologically important natural enemy in the rice ecosystem acting as a dominant predator for planthoppers and leafhoppers in the field. Through Hi-C scaffolding, 1615 scaffolds with a total size of 338.08 Mb were successfully anchored onto 13 chromosomes. The assembled genome size was 345.75 Mb with a final scaffold N50 of 27.58 Mb. Approximately 107.51 Mb of sequences accounting for 31.10% of the genome were identified as repeat elements, and 14,644 protein-coding genes were annotated. Phylogenetic analysis showed that C. lividipennis clustered with other Hemipteran species and diverged from Apolygus lucorum about 66.7 million years ago. Gene families related to detoxification, environmental adaptation and digestion were analysed comparatively with other Hemipteran species, but no significant expansion or contraction was found in C. lividipennis. We also observed male meiosis in C. lividipennis, which showed a typical post-reduction of sex chromosomes and a karyotype of 2n = 22 + XY. As the first natural-enemy genome in the rice ecosystem, the genomic resource of C. lividipennis not only expands our understanding of the multitrophic interactions (host plant-prey-predator), but also provides a genomic basis for better understanding this dominant predator and therefore promotes sustainable rice pest management and food grain production.

The Olfactory Chemosensation of Hematophagous Hemipteran Insects

As one of the most abundant insect orders on earth, most Hemipteran insects are phytophagous, with the few hematophagous exceptions falling into two families: Cimicidae, such as bed bugs, and Reduviidae, such as kissing bugs. Many of these blood-feeding hemipteran insects are known to be realistic or potential disease vectors, presenting both physical and psychological risks for public health. Considerable researches into the interactions between hemipteran insects such as kissing bugs and bed bugs and their human hosts have revealed important information that deepens our understanding of their chemical ecology and olfactory physiology. Sensory mechanisms in the peripheral olfactory system of both insects have now been characterized, with a particular emphasis on their olfactory sensory neurons and odorant receptors. This review summarizes the findings of recent studies of both kissing bugs (including Rhodnius prolixus and Triatoma infestans) and bed bugs (Cimex lectularius), focusing on their chemical ecology and peripheral olfactory systems. Potential chemosensation-based applications for the management of these Hemipteran insect vectors are also discussed.

The Behavioral Response to Heat in the Common Bed Bug, Cimex lectularius (Hemiptera: Cimicidae)

The bed bug, Cimex lectularius L., is a common ectoparasite found to live among its vertebrate hosts. Antennal segments in bugs are critical for sensing multiple cues in the environment for survival. To determine whether the thermo receptors of bed bugs are located on their antennae; innovative bioassays were created to observe the choice between heated and unheated stimuli and to characterize the response of bugs to a heat source. Additionally, the effect of complete antenectomized segments on heat detection were evaluated. Heat, carbon dioxide, and moisture are cues that are found to activate bed bug behavior; a temperature at 38°C was used to assess the direction/degree at which the insect reacts to the change in distance from said stimulus. Using a lightweight spherical ball suspended by air through a vacuum tube, bed bugs and other insects are able to move in 360° while on a stationary point. Noldus EthoVision XT was used to capture video images and to track the bed bugs during 5-min bioassays. A bioassay was created using four Petri dish arenas to observe bed bug attraction to heat based on antennae segments at 40°C. The purpose of this study was to evaluate the effects of heat on complete antenectomized segments of the antennae. The results in this experiment suggest that bed bugs detect and are attracted to heat modulated by nutritional status. Learning the involvement of antennae segments in heat detection will help identify the location and role of thermoreceptors for bed bug host interaction.

Annotation of segmentation pathway genes in the Asian citrus psyllid, Diaphorina citri

Insects have a segmented body plan that is established during embryogenesis when the anterior-posterior (A-P) axis is divided into repeated units by a cascade of gene expression. The cascade is initiated by protein gradients created by translation of maternally provided mRNAs, localized at the anterior and posterior poles of the embryo. Combinations of these proteins activate specific gap genes to divide the embryo into distinct regions along the anterior-posterior axis. Gap genes then activate pair-rule genes, which are usually expressed in parts of every other segment. The pair-rule genes, in turn, activate expression of segment polarity genes in a portion of each segment. The segmentation genes are generally conserved among insects, although there is considerable variation in how they are deployed. We annotated 25 segmentation gene homologs in the Asian citrus psyllid, Diaphorina citri. Most of the genes expected to be present in D. citri based on their phylogenetic distribution in other insects were identified and annotated. Two exceptions were eagle and invected, which are present in at least some hemipterans, but were not found in D. citri. Many of the segmentation pathway genes are likely to be essential for D. citri development, and thus they may be useful targets for gene-based pest control methods.

Utilizing a chromosomal-length genome assembly to annotate the Wnt signaling pathway in the Asian citrus psyllid, Diaphorina citri

The Asian citrus psyllid, Diaphorina citri, is an insect vector that transmits Candidatus Liberibacter asiaticus, the causal agent of the Huanglongbing (HLB), or citrus greening disease. This disease has devastated Florida's citrus industry, and threatens California's industry as well as other citrus producing regions around the world. To find novel solutions to the disease, a better understanding of the vector is needed. The D. citri genome has been used to identify and characterize genes involved in Wnt signaling pathways. Wnt signaling is utilized for many important biological processes in metazoans, such as patterning and tissue generation. Curation based on RNA sequencing data and sequence homology confirms 24 Wnt signaling genes within the D. citri genome, including homologs for beta-catenin, Frizzled receptors, and seven Wnt-ligands. Through phylogenetic analysis, we classify D. citri Wnt ligands as Wg/Wnt1, Wnt5, Wnt6, Wnt7, Wnt10, Wnt11, and WntA. The D. citri version 3.0 genome with chromosomal length scaffolds reveals a conserved Wnt1-Wnt6-Wnt10 gene cluster with a gene configuration like that in Drosophila melanogaster. These findings provide greater insight into the evolutionary history of D. citri and Wnt signaling in this important hemipteran vector. Manual annotation was essential for identifying high quality gene models. These gene models can be used to develop molecular systems, such as CRISPR and RNAi, which target and control psyllid populations to manage the spread of HLB. Manual annotation of Wnt signaling pathways was done as part of a collaborative community annotation project.

Annotation of yellow genes in Diaphorina citri, the vector for Huanglongbing disease

Huanglongbing (HLB), also known as citrus greening disease, is caused by the bacterium Candidatus Liberibacter asiaticus (CLas). It is a serious threat to global citrus production. This bacterium is transmitted by the Asian citrus psyllid, Diaphorina citri (Hemiptera). There are no effective in planta treatments for CLas. Therefore, one strategy is to manage the psyllid population. Manual annotation of the D. citri genome can identify and characterize gene families that could be novel targets for psyllid control. The yellow gene family is an excellent target because yellow genes, which have roles in melanization, are linked to development and immunity. Combined analysis of the genome with RNA-seq datasets, sequence homology, and phylogenetic trees were used to identify and annotate nine yellow genes in the D. citri genome. Manual curation of genes in D. citri provided in-depth analysis of the yellow family among hemipteran insects and provides new targets for molecular control of this psyllid pest. Manual annotation was done as part of a collaborative Citrus Greening community annotation project.

Genome and transcriptome analysis of the mealybug Maconellicoccus hirsutus: Correlation with its unique phenotypes

Mealybugs are aggressive pests with world-wide distribution and are suitable for the study of different phenomena like genomic imprinting and epigenetics. Genomic approaches facilitate these studies in absence of robust genetics in this system. We sequenced, de novo assembled, annotated Maconellicoccus hirsutus genome. We carried out comparative genomics it with four mealybug and eight other insect species, to identify expanded, specific and contracted gene classes that relate to pesticide and desiccation resistance. We identified horizontally transferred genes adding to the mutualism between the mealybug and its endosymbionts. Male and female transcriptome analysis indicates differential expression of metabolic pathway genes correlating with their physiology and the genes for sexual dimorphism. The significantly lower expression of endosymbiont genes in males relates to the depletion of endosymbionts in males during development.

Relaxed selection underlies genome erosion in socially parasitic ant species

Inquiline ants are highly specialized and obligate social parasites that infiltrate and exploit colonies of closely related species. They have evolved many times convergently, are often evolutionarily young lineages, and are almost invariably rare. Focusing on the leaf-cutting ant genus Acromyrmex, we compared genomes of three inquiline social parasites with their free-living, closely-related hosts. The social parasite genomes show distinct signatures of erosion compared to the host lineages, as a consequence of relaxed selective constraints on traits associated with cooperative ant colony life and of inquilines having very small effective population sizes. We find parallel gene losses, particularly in olfactory receptors, consistent with inquiline species having highly reduced social behavioral repertoires. Many of the genomic changes that we uncover resemble those observed in the genomes of obligate non-social parasites and intracellular endosymbionts that branched off into highly specialized, host-dependent niches.

Smelling in the dark: Phylogenomic insights into the chemosensory system of a subterranean beetle

The chemosensory system has experienced relevant changes in subterranean animals, facilitating the perception of specific chemical signals critical to survival in their particular environment. However, the genomic basis of chemoreception in cave-dwelling fauna has been largely unexplored. We generated de novo transcriptomes for antennae and body samples of the troglobitic beetle Speonomus longicornis (whose characters suggest an extreme adaptation to a deep subterranean environment) in order to investigate the evolutionary origin and diversification of the chemosensory gene repertoire across coleopterans through a phylogenomic approach. Our results suggested a diminished diversity of odourant and gustatory gene repertoires compared to polyphagous beetles that inhabit surface habitats. Moreover, S. longicornis showed a large diversity of odourant-binding proteins, suggesting an important role of these proteins in capturing airborne chemical cues. We identified a gene duplication of the ionotropic coreceptor IR25a, a highly conserved single-copy gene in protostomes involved in thermal and humidity sensing. In addition, no homologous genes to sugar receptors or the ionotropic receptor IR41a were detected. Our findings suggest that the chemosensory gene repertoire of this cave beetle may result from adaptation to the highly specific ecological niche it occupies, and that gene duplication and loss may have played an important role in the evolution of gene families involved in chemoreception. Altogether, our results shed light on the genomic basis of chemoreception in a cave-dwelling invertebrate and pave the road towards understanding the genomic underpinnings of adaptation to the subterranean lifestyle at a deeper level.

Jekyll or Hyde? The genome (and more) of Nesidiocoris tenuis, a zoophytophagous predatory bug that is both a biological control agent and a pest

Nesidiocoris tenuis (Reuter) is an efficient predatory biological control agent used throughout the Mediterranean Basin in tomato crops but regarded as a pest in northern European countries. From the family Miridae, it is an economically important insect yet very little is known in terms of genetic information and no genomic or transcriptomic studies have been published. Here, we use a linked-read sequencing strategy on a single female N. tenuis. From this, we assembled the 355 Mbp genome and delivered an ab initio, homology-based and evidence-based annotation. Along the way, the bacterial "contamination" was removed from the assembly. In addition, bacterial lateral gene transfer (LGT) candidates were detected in the N. tenuis genome. The complete gene set is composed of 24 688 genes; the associated proteins were compared to other hemipterans (Cimex lectularis, Halyomorpha halys and Acyrthosiphon pisum). We visualized the genome using various cytogenetic techniques, such as karyotyping, CGH and GISH, indicating a karyotype of 2n = 32. Additional analyses include the localization of 18S rDNA and unique satellite probes as well as pooled sequencing to assess nucleotide diversity and neutrality of the commercial population. This is one of the first mirid genomes to be released and the first of a mirid biological control agent.

Effects of NF-kB Signaling Inhibitors on Bed Bug Resistance to Orally Provisioned Entomopathogenic Bacteria

Bed bugs are globally important pests and there is an ongoing need for the development and improvement of bed bug control tools. Though promising against other insect pests, the exploration of biological methods for bed bug control is limited. Previously, we identified several species of bacteria that have entomopathogenic effects against bed bugs when ingested. We also described the conservation of several antibacterial responses in bed bugs, including the expression of immune effector genes regulated by NF-kB transcription factors through the Toll and immune deficiency (IMD) signaling pathways. Accordingly, we predicted that chemical inhibition of NF-kB signaling could reduce bed bug resistance to orally provisioned entomopathogenic bacteria, potentially improving their effectiveness as biological control agents. In the present study, we administered four small molecule inhibitors of NF-kB signaling (BMS345541, IKK16, IMD0354, Takinib) to bed bugs by feeding them in a blood meal. We then quantified basal mortality and mortality in response to oral infection with two different entomopathogenic bacteria (Pseudomonas entomophila and Bacillus thuringiensis israelensis). None of the NF-kB signaling inhibitors tested increased mortality above control levels when administered alone, suggesting a lack of direct toxicity. However, one inhibitor (IKK16) significantly enhanced the rate of mortality from oral infection with P. entomophila. Enhanced mortality was independent of direct effects of IKK16 on P. entomophila growth in vitro but was associated with higher bacterial loads in vivo (i.e., reduced resistance). Together, these results provide new insight into the regulation of the bed bug immune system and suggest that administration of entomopathogens in combination with inhibition of immune signaling pathways to reduce infection resistance may be effective for biological control of bed bugs.

The genome of the stable fly, Stomoxys calcitrans, reveals potential mechanisms underlying reproduction, host interactions, and novel targets for pest control

Background

The stable fly, Stomoxys calcitrans, is a major blood-feeding pest of livestock that has near worldwide distribution, causing an annual cost of over $2 billion for control and product loss in the USA alone. Control of these flies has been limited to increased sanitary management practices and insecticide application for suppressing larval stages. Few genetic and molecular resources are available to help in developing novel methods for controlling stable flies.

Results

This study examines stable fly biology by utilizing a combination of high-quality genome sequencing and RNA-Seq analyses targeting multiple developmental stages and tissues. In conjunction, 1600 genes were manually curated to characterize genetic features related to stable fly reproduction, vector host interactions, host-microbe dynamics, and putative targets for control. Most notable was characterization of genes associated with reproduction and identification of expanded gene families with functional associations to vision, chemosensation, immunity, and metabolic detoxification pathways.

Conclusions

The combined sequencing, assembly, and curation of the male stable fly genome followed by RNA-Seq and downstream analyses provide insights necessary to understand the biology of this important pest. These resources and new data will provide the groundwork for expanding the tools available to control stable fly infestations. The close relationship of Stomoxys to other blood-feeding (horn flies and Glossina) and non-blood-feeding flies (house flies, medflies, Drosophila) will facilitate understanding of the evolutionary processes associated with development of blood feeding among the Cyclorrhapha.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-00975-9.

Filling in the gaps: A reevaluation of the Lygus hesperus peptidome using an expanded de novo assembled transcriptome and molecular cloning

Peptides are the largest and most diverse class of molecules modulating physiology and behavior. Previously, we predicted a peptidome for the western tarnished plant bug, Lygus hesperus, using transcriptomic data produced from whole individuals. A potential limitation of that analysis was the masking of underrepresented genes, in particular tissue-specific transcripts. Here, we reassessed the L. hesperus peptidome using a more comprehensive dataset comprised of the previous transcriptomic data as well as tissue-specific reads produced from heads and accessory glands. This augmented assembly significantly improves coverage depth providing confirmatory transcripts for essentially all of the previously identified families and new transcripts encoding a number of new peptide precursors corresponding to 14 peptide families. Several families not targeted in our initial study were identified in the expanded assembly, including agatoxin-like peptide, CNMamide, neuropeptide-like precursor 1, and periviscerokinin. To increase confidence in the in silico data, open reading frames of a subset of the newly identified transcripts were amplified using RT-PCR and sequence validated. Further PCR-based profiling of the putative L. hesperus agatoxin-like peptide precursor revealed evidence of alternative splicing with near ubiquitous expression across L. hesperus development, suggesting the peptide serves functional roles beyond that of a toxin. The peptides predicted here, in combination with those identified in our earlier study, expand the L. hesperus peptidome to 42 family members and provide an improved platform for initiating molecular and physiological investigations into peptidergic functionality in this non-model agricultural pest.

Identification and Expression Analysis of G Protein-Coupled Receptors in the Miridae Insect Apolygus lucorum

G protein-coupled receptors (GPCRs) are the largest and most versatile family of transmembrane receptors in the cell and they play a vital role in the regulation of multiple physiological processes. The family Miridae (Hemiptera: Heteroptera) is one of the most diverse families of insects. Until now, information on GPCRs has been lacking in Miridae. Apolygus lucorum, a representative species of the Miridae, is an omnivorous pest that occurs worldwide and is notorious for causing serious damage to various crops and substantial economic losses. By searching the genome, 133 GPCRs were identified in A. lucorum. Compared with other model insects, we have observed GPCR genes to be remarkably expanded in A. lucorum, especially focusing on biogenic amine receptors and neuropeptide receptors. Among these, there is a novel large clade duplicated from known FMRFamide receptors (FMRFaRs). Moreover, the temporal and spatial expression profiles of the 133 genes across developmental stages were determined by transcriptome analysis. Most GPCR genes showed a low expression level in the whole organism of A. lucorum. However, there were a few highly expressed GPCR genes. The highly expressed LW opsins in the head probably relate to nocturning of A. lucorum, and the expression of Cirl at different times and in different tissues indicated it may be involved in growth and development of A. lucorum. We also found C2 leucine-rich repeat-containing GPCRs (LGRs) were mainly distributed in Hemiptera and Phthiraptera among insects. Our study was the first investigation on GPCRs in A. lucorum and it provided a molecular target for the regulation and control of Miridae pests.

Rhodnius prolixus uses the peptidoglycan recognition receptor rpPGRP-LC/LA to detect Gram-negative bacteria and activate the IMD pathway

Insects rely on an innate immune system to recognize and eliminate pathogens. Key components of this system are highly conserved across all invertebrates. To detect pathogens, insects use Pattern recognition receptors (PRRs) that bind to signature motifs on the surface of pathogens called Pathogen Associated Molecular Patterns (PAMPs). In general, insects use peptidoglycan recognition proteins (PGRPs) in the Immune Deficiency (IMD) pathway to detect Gram-negative bacteria, and other PGRPs and Gram-negative binding proteins (GNBPs) in the Toll pathway to detect Gram-positive bacteria and fungi, although there is crosstalk and cooperation between these and other pathways. Once pathogens are recognized, these pathways activate the production of potent antimicrobial peptides (AMPs). Most PRRs in insects have been reported from genome sequencing initiatives but few have been characterized functionally. The initial studies on insect PRRs were done using established dipteran model organisms such as Drosophila melanogaster, but there are differences in the numbers and functional role of PRRs in different insects. Here we describe the genomic repertoire of PGRPs in Rhodnius prolixus, a hemimetabolous hemipteran vector of the parasite Trypanosoma cruzi that causes Chagas disease in humans. Using a de novo transcriptome from the fat body of immune activated insects, we found 5 genes encoding PGRPs. Phylogenetic analysis groups R. prolixus PGRPs with D. melanogaster PGRP-LA, which is involved in the IMD pathway in the respiratory tract. A single R. prolixus PGRP gene encodes isoforms that contain an intracellular region or motif (cryptic RIP Homotypic Interaction Motif-cRHIM) that is involved in the IMD signaling pathway in D. melanogaster. We characterized and silenced this gene using RNAi and show that the PGRPs that contain cRHIMs are involved in the recognition of Gram-negative bacteria, and activation of the IMD pathway in the fat body of R. prolixus, similar to the PGRP-LC of D. melanogaster. This is the first functional characterization of a PGRP containing a cRHIM motif that serves to activate the IMD pathway in a hemimetabolous insect.

Genome Mining and Expression Analysis of Carboxylesterase and Glutathione S-Transferase Genes Involved in Insecticide Resistance in Eggplant Shoot and Fruit Borer, Leucinodes orbonalis (Lepidoptera: Crambidae)

The shoot and fruit borer, Leucinodes orbonalis (Lepidoptera: Crambidae) is the major cause of low productivity in eggplant and insecticides being the mainstay of management of L. orbonalis. However, field control failures are widespread due to the evolution of insecticide resistance. Taking advantage of the whole genome sequence information, the present study investigated the level of insecticide resistance and the expression pattern of individual carboxylesterase (CE) and glutathione S-transferases (GSTs) genes in various field collected populations of L. orbonalis. Dose-mortality bioassays revealed a very high level of resistance development against fenvalerate (48.2–160-fold), phosalone (94-534.6-fold), emamectin benzoate (7.2–55-fold), thiodicarb (9.64–22.7-fold), flubendiamide (187.4–303.0-fold), and chlorantraniliprole (1.6–8.6-fold) in field populations as compared to laboratory-reared susceptible iso-female colony (Lo-S). Over-production of detoxification enzymes viz., CE and GST were evident upon enzyme assays. Mining of the draft genome of L. orbonalis yielded large number of genes potentially belonging to the CE and GST gene families with known history of insecticide resistance in other insects. Subsequent RT-qPCR studies on relative contribution of individual genes revealed over-expression of numerous GSTs and few CEs in field populations, indicating their possible involvement of metabolic enzymes in insecticide resistance. The genomic information will facilitate the development of novel resistance management strategies against this pest.