Genome sequence of the pea aphid Acyrthosiphon pisum

Rapid Targeted Assembly of the Proteome Reveals Evolutionary Variation of GC Content in Avian Lice

Nucleotide base composition plays an influential role in the molecular mechanisms involved in gene function, phenotype, and amino acid composition. GC content (proportion of guanine and cytosine in DNA sequences) shows a high level of variation within and among species. Many studies measure GC content in a small number of genes, which may not be representative of genome-wide GC variation. One challenge when assembling extensive genomic data sets for these studies is the significant amount of resources (monetary and computational) associated with data processing, and many bioinformatic tools have not been optimized for resource efficiency. Using a high-performance computing (HPC) cluster, we manipulated resources provided to the targeted gene assembly program, automated target restricted assembly method (aTRAM), to determine an optimum way to run the program to maximize resource use. Using our optimum assembly approach, we assembled and measured GC content of all of the protein-coding genes of a diverse group of parasitic feather lice. Of the 499 426 genes assembled across 57 species, feather lice were GC-poor (mean GC = 42.96%) with a significant amount of variation within and between species (GC range = 19.57%-73.33%). We found a significant correlation between GC content and standard deviation per taxon for overall GC and GC3, which could indicate selection for G and C nucleotides in some species. Phylogenetic signal of GC content was detected in both GC and GC3. This research provides a large-scale investigation of GC content in parasitic lice laying the foundation for understanding the basis of variation in base composition across species.

Genome assembly and transcriptomic analysis to elucidate the ability of Nasonovia ribisnigri to break host plant resistance

Aphid genomic resources enable the study of complex life history traits and provide information on vector biology, host adaption and speciation. The currant-lettuce aphid (Nasonovia ribisnigri (Hemiptera: Aphididae) (Mosley)) is a cosmopolitan pest of outdoor lettuce (Lactuca sativa (Asterales: Asteraceae) (Linnaeus)). Until recently, the use of resistant cultivars was an effective method for managing N. ribisnigri. A resistant cultivar containing a single gene (Nr-locus), introduced in the 1980s, conferred complete resistance to feeding. Overreliance of this Nr-locus in lettuce resulted in N. ribisnigri's ability to break resistance mechanism, with first reports during 2003. Our work attempts to understand which candidate gene(s) are associated with this resistance-breaking mechanism. We present two de novo draft assembles for N. ribisnigri genomes, corresponding to both avirulent (Nr-locus susceptible) and virulent (Nr-locus resistant) biotypes. Changes in gene expression of the two N. ribisnigri biotypes were investigated using transcriptomic analyses of RNA-sequencing (RNA-seq) data to understand the potential mechanisms of resistance to the Nr-locus in lettuce. The draft genome assemblies were 94.2% and 91.4% complete for the avirulent and virulent biotypes, respectively. Out of the 18,872 differentially expressed genes, a single gene/locus was identified in N. ribisnigri that was shared between two resistant-breaking biotypes. This locus was further explored and validated in Real-Time Quantitative Reverse Transcription PCR (qRT-PCR) experiments and has predicted localisations in both the cytoplasm and nucleus. This is the first study to provide evidence that a single gene/locus is likely responsible for the ability of N. ribisnigri to overcome the Nr-locus resistance in the lettuce host.

Pathogen-microbiome interactions and the virulence of an entomopathogenic fungus

Bacteria shape interactions between hosts and fungal pathogens. In some cases, bacteria associated with fungi are essential for pathogen virulence. In other systems, host-associated microbiomes confer resistance against fungal pathogens. We studied an aphid-specific entomopathogenic fungus called Pandora neoaphidis in the context of both host and pathogen microbiomes. Aphids host several species of heritable bacteria, some of which confer resistance against Pandora. We first found that spores that emerged from aphids that harbored protective bacteria were less virulent against subsequent hosts and did not grow on plate media. We then used 16S amplicon sequencing to study the bacterial microbiome of fungal mycelia and spores during plate culturing and host infection. We found that the bacterial community is remarkably stable in culture despite dramatic changes in pathogen virulence. Last, we used an experimentally transformed symbiont of aphids to show that Pandora can acquire host-associated bacteria during infection. Our results uncover new roles for bacteria in the dynamics of aphid-pathogen interactions and illustrate the importance of the broader microbiological context in studies of fungal pathogenesis.

IMPORTANCE

Entomopathogenic fungi play important roles in the population dynamics of many insect species. Understanding the factors shaping entomopathogen virulence is critical for agricultural management and for the use of fungi in pest biocontrol. We show that heritable bacteria in aphids, which confer protection to their hosts against fungal entomopathogens, influence virulence against subsequent hosts. Aphids reproduce asexually and are typically surrounded by genetically identical offspring, and thus these effects likely shape the dynamics of fungal disease in aphid populations. Furthermore, fungal entomopathogens are known to rapidly lose virulence in lab culture, complicating their laboratory use. We show that this phenomenon is not driven by changes in the associated bacterial microbiome. These results contribute to our broader understanding of the aphid model system and shed light on the biology of the Entomophthorales-an important but understudied group of fungi.

Immune Gene Repertoire of Soft Scale Insects (Hemiptera: Coccidae)

Insects possess an effective immune system, which has been extensively characterized in several model species, revealing a plethora of conserved genes involved in recognition, signaling, and responses to pathogens and parasites. However, some taxonomic groups, characterized by peculiar trophic niches, such as plant-sap feeders, which are often important pests of crops and forestry ecosystems, have been largely overlooked regarding their immune gene repertoire. Here we annotated the immune genes of soft scale insects (Hemiptera: Coccidae) for which omics data are publicly available. By using immune genes of aphids and Drosophila to query the genome of Ericerus pela, as well as the transcriptomes of Ceroplastes cirripediformis and Coccus sp., we highlight the lack of peptidoglycan recognition proteins, galectins, thaumatins, and antimicrobial peptides in Coccidae. This work contributes to expanding our knowledge about the evolutionary trajectories of immune genes and offers a list of promising candidates for developing new control strategies based on the suppression of pests' immunity through RNAi technologies.

Transcriptomic and Metatranscriptomic Analyses Provide New Insights into the Response of the Pea Aphid Acyrthosiphon pisum (Hemiptera: Aphididae) to Acetamiprid

Acetamiprid is a broad-spectrum neonicotinoid insecticide used in agriculture to control aphids. While recent studies have documented resistance to acetamiprid in several aphid species, the underlying mechanisms are still not fully understood. In this study, we analyzed the transcriptome and metatranscriptome of a laboratory strain of the pea aphid, Acyrthosiphon pisum (Harris, 1776), with reduced susceptibility to acetamiprid after nine generations of exposure to identify candidate genes and the microbiome involved in the adaptation process. Sequencing of the transcriptome of both selected (RS) and non-selected (SS) strains allowed the identification of 14,858 genes and 4938 new transcripts. Most of the differentially expressed genes were associated with catalytic activities and metabolic pathways involving carbon and fatty acids. Specifically, alcohol-forming fatty acyl-CoA reductase (FAR) and acyl-CoA synthetase (ACSF2), both involved in the synthesis of epidermal wax layer components, were significantly upregulated in RS, suggesting that adaptation to acetamiprid involves the synthesis of a thicker protective layer. Metatranscriptomic analyses revealed subtle shifts in the microbiome of RS. These results contribute to a deeper understanding of acetamiprid adaptation by the pea aphid and provide new insights for aphid control strategies.

Toxicological, biochemical, and in silico investigations of three trehalase inhibitors for new ways to control aphids

Insect trehalases have been identified as promising new targets for pest control. These key enzymes are involved in trehalose hydrolysis and plays an important role in insect growth and development. In this contribution, plant and microbial compounds, namely validamycin A, amygdalin, and phloridzin, were evaluated for their effect, through trehalase inhibition, on Acyrthosiphon pisum aphid. The latter is part of the Aphididae family, main pests as phytovirus vectors and being very harmful for crops. Validamycin A was confirmed as an excellent trehalase inhibitor with an half maximal inhibitory concentration and inhibitor constant of 2.2 × 10-7 and 5 × 10-8 M, respectively, with a mortality rate of ~80% on a A. pisum population. Unlike validamycin A, the insect lethal efficacy of amygdalin and phloridzin did not correspond to their trehalase inhibition, probably due to their hydrolysis by insect β-glucosidases. Our docking studies showed that none of the three compounds can bind to the trehalase active site, unlike their hydrolyzed counterparts, that is, validoxylamine A, phloretin, and prunasin. Validoxylamine A would be by far the best trehalase binder, followed by phloretin and prunasin.

Ingested soil bacteria breach gut epithelia and prime systemic immunity in an insect

Insects lack acquired immunity and were thought to have no immune memory, but recent studies reported a phenomenon called immune priming, wherein sublethal dose of pathogens or nonpathogenic microbes stimulates immunity and prevents subsequential pathogen infection. Although the evidence for insect immune priming is accumulating, the underlying mechanisms are still unclear. The bean bug Riptortus pedestris acquires its gut microbiota from ambient soil and spatially structures them into a multispecies and variable community in the anterior midgut and a specific, monospecies Caballeronia symbiont population in the posterior region. We demonstrate that a particular Burkholderia strain colonizing the anterior midgut stimulates systemic immunity by penetrating gut epithelia and migrating into the hemolymph. The activated immunity, consisting of a humoral and a cellular response, had no negative effect on the host fitness, but on the contrary protected the insect from subsequent infection by pathogenic bacteria. Interruption of contact between the Burkholderia strain and epithelia of the gut weakened the host immunity back to preinfection levels and made the insects more vulnerable to microbial infection, demonstrating that persistent acquisition of environmental bacteria is important to maintain an efficient immunity.

Localization of nitric oxide-producing hemocytes in Aedes and Culex mosquitoes infected with bacteria

Mosquitoes are significant vectors of various pathogens. Unlike vertebrates, insects rely solely on innate immunity. Hemocytes play a crucial role in the cellular part of the innate immune system. The gaseous radical nitric oxide (NO) produced by hemocytes acts against pathogens and also functions as a versatile transmitter in both the immune and nervous systems, utilizing cyclic guanosine monophosphate (cGMP) as a second messenger. This study conducted a parallel comparison of NO synthase (NOS) expression and NO production in hemocytes during Escherichia coli K12 infection in four vector species: Aedes aegypti, Aedes albopictus, Culex pipiens molestus, and Culex pipiens quinquefasciatus. Increased NOS expression by NADPH diaphorase (NADPHd) staining and NO production by immunofluorescence against the by-product L-citrulline were observed in infected mosquito hemocytes distributed throughout the abdomens. NADPHd activity and citrulline labeling were particularly found in periostial hemocytes near the heart, but also on the ventral nerve chord (VNC). Pericardial cells of Ae. aegypti and Cx. p. molestus showed increased citrulline immunofluorescence, suggesting their involvement in the immune response. Oenocytes displayed strong NADPHd and citrulline labeling independent of infection status. This comparative study, consistent with findings in other species, suggests a widespread phenomenon of NO's role in hemocyte responses during E. coli infection. Found differences within and between genera highlight the importance of species-specific investigations.

Endosymbiont Tremblaya phenacola influences the reproduction of cotton mealybugs by regulating the mechanistic target of rapamycin pathway

The intricate evolutionary dynamics of endosymbiotic relationships result in unique characteristics among the genomes of symbionts, which profoundly influence host insect phenotypes. Here, we investigated an endosymbiotic system in Phenacoccus solenopsis, a notorious pest of the subfamily Phenacoccinae. The endosymbiont, "Candidatus Tremblaya phenacola" (T. phenacola PSOL), persisted throughout the complete life cycle of female hosts and was more active during oviposition, whereas there was a significant decline in abundance after pupation in males. Genome sequencing yielded an endosymbiont genome of 221.1 kb in size, comprising seven contigs and originating from a chimeric arrangement between betaproteobacteria and gammaproteobacteria. A comprehensive analysis of amino acid metabolic pathways demonstrated complementarity between the host and endosymbiont metabolism. Elimination of T. phenacola PSOL through antibiotic treatment significantly decreased P. solenopsis fecundity. Weighted gene coexpression network analysis demonstrated a correlation between genes associated with essential amino acid synthesis and those associated with host meiosis and oocyte maturation. Moreover, altering endosymbiont abundance activated the host mechanistic target of rapamycin pathway, suggesting that changes in the amino acid abundance affected the host reproductive capabilities via this signal pathway. Taken together, these findings demonstrate a mechanism by which the endosymbiont T. phenacola PSOL contributed to high fecundity in P. solenopsis and provide new insights into nutritional compensation and coevolution of the endosymbiotic system.

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

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

Phylloxera and Aphids Show Distinct Features of Genome Evolution Despite Similar Reproductive Modes

Genomes of aphids (family Aphididae) show several unusual evolutionary patterns. In particular, within the XO sex determination system of aphids, the X chromosome exhibits a lower rate of interchromosomal rearrangements, fewer highly expressed genes, and faster evolution at nonsynonymous sites compared with the autosomes. In contrast, other hemipteran lineages have similar rates of interchromosomal rearrangement for autosomes and X chromosomes. One possible explanation for these differences is the aphid's life cycle of cyclical parthenogenesis, where multiple asexual generations alternate with 1 sexual generation. If true, we should see similar features in the genomes of Phylloxeridae, an outgroup of aphids which also undergoes cyclical parthenogenesis. To investigate this, we generated a chromosome-level assembly for the grape phylloxera, an agriculturally important species of Phylloxeridae, and identified its single X chromosome. We then performed synteny analysis using the phylloxerid genome and 30 high-quality genomes of aphids and other hemipteran species. Unexpectedly, we found that the phylloxera does not share aphids' patterns of chromosome evolution. By estimating interchromosomal rearrangement rates on an absolute time scale, we found that rates are elevated for aphid autosomes compared with their X chromosomes, but this pattern does not extend to the phylloxera branch. Potentially, the conservation of X chromosome gene content is due to selection on XO males that appear in the sexual generation. We also examined gene duplication patterns across Hemiptera and uncovered horizontal gene transfer events contributing to phylloxera evolution.

Untangling an insect's virome from its endogenous viral elements

BACKGROUND

Insects are an important reservoir of viral biodiversity, but the vast majority of viruses associated with insects have not been discovered. Recent studies have employed high-throughput RNA sequencing, which has led to rapid advances in our understanding of insect viral diversity. However, insect genomes frequently contain transcribed endogenous viral elements (EVEs) with significant homology to exogenous viruses, complicating the use of RNAseq for viral discovery.

METHODS

In this study, we used a multi-pronged sequencing approach to study the virome of an important agricultural pest and prolific vector of plant pathogens, the potato aphid Macrosiphum euphorbiae. We first used rRNA-depleted RNAseq to characterize the microbes found in individual insects. We then used PCR screening to measure the frequency of two heritable viruses in a local aphid population. Lastly, we generated a quality draft genome assembly for M. euphorbiae using Illumina-corrected Nanopore sequencing to identify transcriptionally active EVEs in the host genome.

RESULTS

We found reads from two insect-specific viruses (a Flavivirus and an Ambidensovirus) in our RNAseq data, as well as a parasitoid virus (Bracovirus), a plant pathogenic virus (Tombusvirus), and two phages (Acinetobacter and APSE). However, our genome assembly showed that part of the 'virome' of this insect can be attributed to EVEs in the host genome.

CONCLUSION

Our work shows that EVEs have led to the misidentification of aphid viruses from RNAseq data, and we argue that this is a widespread challenge for the study of viral diversity in insects.

Phylloxera and aphids show distinct features of genome evolution despite similar reproductive modes

Genomes of aphids (family Aphididae) show several unusual evolutionary patterns. In particular, within the XO sex determination system of aphids, the X chromosome exhibits a lower rate of interchromosomal rearrangements, fewer highly expressed genes, and faster evolution at nonsynonymous sites compared to the autosomes. In contrast, other hemipteran lineages have similar rates of interchromosomal rearrangement for autosomes and X chromosomes. One possible explanation for these differences is the aphid's life cycle of cyclical parthenogenesis, where multiple asexual generations alternate with one sexual generation. If true, we should see similar features in the genomes of Phylloxeridae, an outgroup of aphids which also undergoes cyclical parthenogenesis. To investigate this, we generated a chromosome-level assembly for the grape phylloxera, an agriculturally important species of Phylloxeridae, and identified its single X chromosome. We then performed synteny analysis using the phylloxerid genome and 30 high-quality genomes of aphids and other hemipteran species. Unexpectedly, we found that the phylloxera does not share aphids' patterns of chromosome evolution. By estimating interchromosomal rearrangement rates on an absolute time scale, we found that rates are elevated for aphid autosomes compared to their X chromosomes, but this pattern does not extend to the phylloxera branch. Potentially, the conservation of X chromosome gene content is due to selection on XO males that appear in the sexual generation. We also examined gene duplication patterns across Hemiptera and uncovered horizontal gene transfer events contributing to phylloxera evolution.

Gene identification and RNAi-silencing of p62/SQSTM1 in the vector Rhodnius prolixus reveals a high degree of sequence conservation but no apparent deficiency-related phenotypes in vitellogenic females

Autophagy and the ubiquitin-proteasome system (UPS) are important cellular mechanisms that coordinate protein degradation essential for proteostasis. P62/SQSTM1 is a receptor cargo protein able to deliver ubiquitinated targets to the proteasome proteolytic complex and/or to the autophagosome. In the insect vector of Chagas disease, Rhodnius prolixus, previous works have shown that the knockdown of different autophagy-related genes (ATGs) and ubiquitin-conjugating enzymes resulted in abnormal oogenesis phenotypes and embryo lethality. Here, we investigate the role of the autophagy/UPS adaptor protein p62 during the oogenesis and reproduction of this vector. We found that R. prolixus presents one isoform of p62 encoded by a non-annotated gene. The predicted protein presents the domain architecture anticipated for p62: PB1 (N-term), ZZ-finger, and UBA (C-term) domains, and phylogenetic analysis showed that this pattern is highly conserved within insects. Using parental RNAi, we found that although p62 is expressed in the ovary, midgut, and fat body of adult females, systemic silencing of this gene did not result in any apparent phenotypes under in-house conditions. The insects' overall levels of blood meal digestion, lifespan, yolk protein production, oviposition, and embryo viability were not altered when compared to controls. Because it is known that autophagy and UPS can undergo compensatory mechanisms, we asked whether the silencing of p62 was triggering adaptative changes in the expression of genes of the autophagy, UPS, and the unfolded protein response (UPR) and found that only ATG1 was slightly up regulated in the ovaries of silenced females. In addition, experiments to further investigate the role of p62 in insects previously silenced for the E1-conjugating enzyme (a condition known to trigger the upregulation of p62), also did not result in any apparent phenotypes in vitellogenic females.

Hybridisation has shaped a recent radiation of grass-feeding aphids

BACKGROUND

Aphids are common crop pests. These insects reproduce by facultative parthenogenesis involving several rounds of clonal reproduction interspersed with an occasional sexual cycle. Furthermore, clonal aphids give birth to live young that are already pregnant. These qualities enable rapid population growth and have facilitated the colonisation of crops globally. In several cases, so-called "super clones" have come to dominate agricultural systems. However, the extent to which the sexual stage of the aphid life cycle has shaped global pest populations has remained unclear, as have the origins of successful lineages. Here, we used chromosome-scale genome assemblies to disentangle the evolution of two global pests of cereals-the English (Sitobion avenae) and Indian (Sitobion miscanthi) grain aphids.

RESULTS

Genome-wide divergence between S. avenae and S. miscanthi is low. Moreover, comparison of haplotype-resolved assemblies revealed that the S. miscanthi isolate used for genome sequencing is likely a hybrid, with one of its diploid genome copies closely related to S. avenae (~ 0.5% divergence) and the other substantially more divergent (> 1%). Population genomics analyses of UK and China grain aphids showed that S. avenae and S. miscanthi are part of a cryptic species complex with many highly differentiated lineages that predate the origins of agriculture. The complex consists of hybrid lineages that display a tangled history of hybridisation and genetic introgression.

CONCLUSIONS

Our analyses reveal that hybridisation has substantially contributed to grain aphid diversity, and hence, to the evolutionary potential of this important pest species. Furthermore, we propose that aphids are particularly well placed to exploit hybridisation events via the rapid propagation of live-born "frozen hybrids" via asexual reproduction, increasing the likelihood of hybrid lineage formation.

Knockdown of carnitine palmitoyltransferase I (CPT1) reduces fat body lipid mobilization and resistance to starvation in the insect vector Rhodnius prolixus

The energy stored in fatty acids is essential for several critical activities of insects, such as embryogenesis, oviposition, and flight. Rhodnius prolixus is an obligatory hematophagous hemipteran and vector of Chagas disease, and it feeds infrequently on very large blood meals. As digestion slowly occurs, lipids are synthesized and accumulate in the fat body, mainly as triacylglycerol, in lipid droplets. Between feeding bouts, proper mobilization and oxidation of stored lipids are crucial for survival, and released fatty acids are oxidized by mitochondrial β-oxidation. Carnitine palmitoyl transferase I (CPT1) is the enzyme that catalyzes the first reaction of the carnitine shuttle, where the activated fatty acid, acyl-CoA, is converted to acyl-carnitine to be transported into the mitochondria. Here, we investigated the role of CPT1 in lipid metabolism and in resistance to starvation in Rhodnius prolixus. The expression of the CPT1 gene (RhoprCpt1) was determined in the organs of adult females on the fourth day after a blood meal, and the flight muscle showed higher expression levels than the ovary, fat body, and anterior and posterior midgut. RhoprCpt1 expression in the fat body dramatically decreased after feeding, and started to increase again 10 days later, but no changes were observed in the flight muscle. β-oxidation rates were determined in flight muscle and fat body homogenates with the use of ³H-palmitate, and in unfed females, they were higher in the flight muscle. In the fat body, lipid oxidation activity did not show any variation before or at different days after feeding, and was not affected by the presence of etomoxir or malonyl-CoA. We used RNAi and generated RhoprCPT1-deficient insects, which surprisingly did not show a decrease in measured ³H-palmitate oxidation rates. However, the RNAi-knockdown females presented increased amounts of triacylglycerol and larger lipid droplets in the fat body, but not in the flight muscle. When subjected to starvation, these insects had a shorter lifespan. These results indicated that the inhibition of RhoprCpt1 expression compromised lipid mobilization and affected resistance to starvation.

Pest status, molecular evolution, and epigenetic factors derived from the genome assembly of Frankliniella fusca, a thysanopteran phytovirus vector

BACKGROUND

The tobacco thrips (Frankliniella fusca Hinds; family Thripidae; order Thysanoptera) is an important pest that can transmit viruses such as the tomato spotted wilt orthotospovirus to numerous economically important agricultural row crops and vegetables. The structural and functional genomics within the order Thysanoptera has only begun to be explored. Within the > 7000 known thysanopteran species, the melon thrips (Thrips palmi Karny) and the western flower thrips (Frankliniella occidentalis Pergrande) are the only two thysanopteran species with assembled genomes.

RESULTS

A genome of F. fusca was assembled by long-read sequencing of DNA from an inbred line. The final assembly size was 370 Mb with a single copy ortholog completeness of ~ 99% with respect to Insecta. The annotated genome of F. fusca was compared with the genome of its congener, F. occidentalis. Results revealed many instances of lineage-specific differences in gene content. Analyses of sequence divergence between the two Frankliniella species' genomes revealed substitution patterns consistent with positive selection in ~ 5% of the protein-coding genes with 1:1 orthologs. Further, gene content related to its pest status, such as xenobiotic detoxification and response to an ambisense-tripartite RNA virus (orthotospovirus) infection was compared with F. occidentalis. Several F. fusca genes related to virus infection possessed signatures of positive selection. Estimation of CpG depletion, a mutational consequence of DNA methylation, revealed that F. fusca genes that were downregulated and alternatively spliced in response to virus infection were preferentially targeted by DNA methylation. As in many other insects, DNA methylation was enriched in exons in Frankliniella, but gene copies with homology to DNA methyltransferase 3 were numerous and fragmented. This phenomenon seems to be relatively unique to thrips among other insect groups.

CONCLUSIONS

The F. fusca genome assembly provides an important resource for comparative genomic analyses of thysanopterans. This genomic foundation allows for insights into molecular evolution, gene regulation, and loci important to agricultural pest status.

Silencing of the 20S proteasomal subunit-α6 triggers full oogenesis arrest and increased mRNA levels of the selective autophagy adaptor protein p62/SQSTM1 in the ovary of the vector Rhodnius prolixus

The high reproductive rates of insects contribute significantly to their ability to act as vectors of a variety of vector-borne diseases. Therefore, it is strategically critical to find molecular targets with biotechnological potential through the functional study of genes essential for insect reproduction. The ubiquitin-proteasome system is a vital degradative pathway that contributes to the maintenance of regular eukaryotic cell proteostasis. This mechanism involves the action of enzymes to covalently link ubiquitin to proteins that are meant to be delivered to the 26S proteasome and broken down. The 26S proteasome is a large protease complex (including the 20S and 19S subcomplexes) that binds, deubiquitylates, unfolds, and degrades its substrates. Here, we used bioinformatics to identify the genes that encode the seven α and β subunits of the 20S proteasome in the genome of R. prolixus and learned that those transcripts are accumulated into mature oocytes. To access proteasome function during oogenesis, we conducted RNAi functional tests employing one of the 20S proteasome subunits (Prosα6) as a tool to suppress 20S proteasomal activity. We found that Prosα6 silencing resulted in no changes in TAG buildup in the fat body and unaffected availability of yolk proteins in the hemolymph of vitellogenic females. Despite this, the silencing of Prosα6 culminated in the impairment of oocyte maturation at the early stages of oogenesis. Overall, we discovered that proteasome activity is especially important for the signals that initiate oogenesis in R. prolixus and discuss in what manner further investigations on the regulation of proteasome assembly and activity might contribute to the unraveling of oogenesis molecular mechanisms and oocyte maturation in this vector.

Pigment-dispersing factor is present in circadian clock neurons of pea aphids and may mediate photoperiodic signalling to insulin-producing cells

The neuropeptide pigment-dispersing factor (PDF) plays a pivotal role in the circadian clock of most Ecdysozoa and is additionally involved in the timing of seasonal responses of several photoperiodic species. The pea aphid, Acyrthosiphon pisum, is a paradigmatic photoperiodic species with an annual life cycle tightly coupled to the seasonal changes in day length. Nevertheless, PDF could not be identified in A. pisum so far. In the present study, we identified a PDF-coding gene that has undergone significant changes in the otherwise highly conserved insect C-terminal amino acid sequence. A newly generated aphid-specific PDF antibody stained four neurons in each hemisphere of the aphid brain that co-express the clock protein Period and have projections to the pars lateralis that are highly plastic and change their appearance in a daily and seasonal manner, resembling those of the fruit fly PDF neurons. Most intriguingly, the PDF terminals overlap with dendrites of the insulin-like peptide (ILP) positive neurosecretory cells in the pars intercerebralis and with putative terminals of Cryptochrome (CRY) positive clock neurons. Since ILP has been previously shown to be crucial for seasonal adaptations and CRY might serve as a circadian photoreceptor vital for measuring day length, our results suggest that PDF plays a critical role in aphid seasonal timing.

Chromosome-level genome assembly of the spotted alfalfa aphid Therioaphis trifolii

The spotted alfalfa aphid (SAA, Therioaphis trifolii) (Hemiptera: Aphididae) is a destructive pest of cultivated alfalfa (Medicago sativa L.) that leads to large financial losses in the livestock industry around the world. Here, we present a chromosome-scale genome assembly of T. trifolii, the first genome assembly for the aphid subfamily Calaphidinae. Using PacBio long-read sequencing, Illumina sequencing, and Hi-C scaffolding techniques, a 541.26 Mb genome was generated, with 90.01% of the assembly anchored into eight scaffolds, and the contig and scaffold N50 are 2.54 Mb and 44.77 Mb, respectively. BUSCO assessment showed a completeness score of 96.6%. A total of 13,684 protein-coding genes were predicted. The high-quality genome assembly of T. trifolii not only provides a genomic resource for the more complete analysis of aphid evolution, but also provides insights into the ecological adaptation and insecticide resistance of T. trifolii.

The transcriptome of Icerya aegyptiaca (Hemiptera: Monophlebidae) and comparison with neococcoids reveal genetic clues of evolution in the scale insects

BACKGROUND

Scale insects are worldwide sap-sucking parasites, which can be distinguished into neococcoids and non-neococcoids. Neococcoids are monophyletic with a peculiar reproductive system, paternal genome elimination (PGE). Different with neococcoids, Iceryini, a tribe in non-neococcoids including several damaging pests, has abdominal spiracles, compound eyes in males, relatively abundant wax, unique hermaphrodite system, and specific symbionts. However, the current studies on the gene resources and genomic mechanism of scale insects are mainly limited in the neococcoids, and lacked of comparison in an evolution frame.

RESULT

We sequenced and de novo assembled a transcriptome of Icerya aegyptiaca (Douglas), a worldwide pest of Iceryini, and used it as representative of non-neococcoids to compare with the genomes or transcriptomes of other six species from different families of neococcoids. We found that the genes under positive selection or negative selection intensification (simplified as "selected genes" below) in I. aegyptiaca included those related to neurogenesis and development, especially eye development. Some genes related to fatty acid biosynthesis were unique in its transcriptome with relatively high expression and not detected in neococcoids. These results may indicate a potential link to the unique structures and abundant wax of I. aegyptiaca compared with neococcoids. Meanwhile, genes related to DNA repair, mitosis, spindle, cytokinesis and oogenesis, were included in the selected genes in I. aegyptiaca, which is possibly associated with cell division and germ cell formation of the hermaphrodite system. Chromatin-related process were enriched from selected genes in neococcoids, along with some mitosis-related genes also detected, which may be related to their unique PGE system. Moreover, in neococcoid species, male-biased genes tend to undergo negative selection relaxation under the PGE system. We also found that the candidate horizontally transferred genes (HTGs) in the scale insects mainly derived from bacteria and fungi. bioD and bioB, the two biotin-synthesizing HTGs were exclusively found in the scale insects and neococcoids, respectively, which possibly show potential demand changes in the symbiotic relationships.

CONCLUSION

Our study reports the first I. aegyptiaca transcriptome and provides preliminary insights for the genetic change of structures, reproductive systems and symbiont relationships at an evolutionary aspect. This will provide a basis for further research and control of scale insects.

The mechanism of Cry41-related toxin against Myzus persicae based on its interaction with Buchnera-derived ATP-dependent 6-phosphofructokinase

BACKGROUND

Myzus persicae (Hemiptera: Aphididae) is one of the most notorious pests of many crops worldwide. Most Cry toxins produced by Bacillus thuringiensis show very low toxicity to M. persicae; however, a study showed that Cry41-related toxin had moderate toxic activity against M. persicae. In our previous work, potential Cry41-related toxin-binding proteins in M. persicae were identified, including cathepsin B, calcium-transporting ATPase, and Buchnera-derived ATP-dependent 6-phosphofructokinase (PFKA). Buchnera is an endosymbiont present in almost all aphids and it provides necessary nutrients for aphid growth. This study investigated the role of Buchnera-derived PFKA in Cry41-related toxicity against M. persicae.

RESULTS

In this study, recombinant PFKA was expressed and purified, and in vitro assays revealed that PFKA bound to Cry41-related toxin, and Cry41-related toxin at 25 μg ml^-1 significantly inhibited the activity of PFKA. In addition, when M. persicae was treated with 30 μg ml^-1 of Cry41-related toxin for 24 h, the expression of dnak, a single-copy gene in Buchnera, was significantly decreased, indicating a decrease in the number of Buchnera.

CONCLUSION

Our results suggest that Cry41-related toxin interacts with Buchnera-derived PFKA to inhibit its enzymatic activity and likely impair cell viability, resulting in a decrease in the number of Buchnera, and finally leading to M. persicae death. These findings open up new perspectives in our understanding of the mode of action of Cry toxins and are useful in helping improve Cry toxicity for aphid control. © 2023 Society of Chemical Industry.

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.

Comparative transcriptomics of aphid species that diverged > 22 MYA reveals genes that are important for the maintenance of their symbiosis

Most plant-sap feeding insects have obligate relationships with maternally transmitted bacteria. Aphids require their nutritional endosymbiont, Buchnera aphidicola, for the production of essential amino acids. Such endosymbionts are harbored inside of specialized insect cells called bacteriocytes. Here, we use comparative transcriptomics of bacteriocytes between two recently diverged aphid species, Myzus persicae and Acyrthosiphon pisum, to identify key genes that are important for the maintenance of their nutritional mutualism. The majority of genes with conserved expression profiles in M. persicae and A. pisum are for orthologs previously identified in A. pisum to be important for the symbiosis. However, asparaginase which produces aspartate from asparagine was significantly up-regulated only in A. pisum bacteriocytes, potentially because Buchnera of M. persicae encodes its own asparaginase enzyme unlike Buchnera of A. pisum resulting in Buchnera of A. pisum to be dependent on its aphid host for aspartate. One-to-one orthologs that explained the most amount of variation for bacteriocyte specific mRNA expression for both species includes a collaborative gene for methionine biosynthesis, multiple transporters, a horizontally transmitted gene, and secreted proteins. Finally, we highlight species-specific gene clusters which may contribute to host adaptations and/or accommodations in gene regulation to changes in the symbiont or the symbiosis.

Horizontally transferred genes as RNA interference targets for aphid and whitefly control

RNA interference (RNAi)-based technologies are starting to be commercialized as a new approach for agricultural pest control. Horizontally transferred genes (HTGs), which have been transferred into insect genomes from viruses, bacteria, fungi or plants, are attractive targets for RNAi-mediated pest control. HTGs are often unique to a specific insect family or even genus, making it unlikely that RNAi constructs targeting such genes will have negative effects on ladybugs, lacewings and other beneficial predatory insect species. In this study, we sequenced the genome of a red, tobacco-adapted isolate of Myzus persicae (green peach aphid) and bioinformatically identified 30 HTGs. We then used plant-mediated virus-induced gene silencing (VIGS) to show that several HTGs of bacterial and plant origin are important for aphid growth and/or survival. Silencing the expression of fungal-origin HTGs did not affect aphid survivorship but decreased aphid reproduction. Importantly, although there was uptake of plant-expressed RNA by Coccinella septempunctata (seven-spotted ladybugs) via the aphids that they consumed, we did not observe negative effects on ladybugs from aphid-targeted VIGS constructs. To demonstrate that this approach is more broadly applicable, we also targeted five Bemisia tabaci (whitefly) HTGs using VIGS and demonstrated that knockdown of some of these genes affected whitefly survival. As functional HTGs have been identified in the genomes of numerous pest species, we propose that these HTGs should be explored further as efficient and safe targets for control of insect pests using plant-mediated RNA interference.

Buchnera breaks the specialization of the cotton-specialized aphid (Aphis gossypii) by providing nutrition through zucchini

The cotton aphid, Aphis gossypii Glover, is a species of polyphagous aphid with many biotypes, and its host transfer has always been the focus of research on the control of cotton aphid. An important factor affecting aphid specialization is the nutritional association with microbial symbionts that provide the host with nutrients lacking in the diet. We analyzed the microbial composition and biodiversity of reared on zucchini for 10 generations (T1-T10) and cotton as a control (CK), by high-throughput Illumina sequencing of 16S ribosomal RNA genes. The findings showed that the change in plant hosts decreased the richness and variety of microbial species. Regardless of whether the plant host is altered or not, Proteobacteria and Firmicutes are the predominate phyla in cotton-specialized aphid. Additionally, cotton-specialized aphids that live in zucchini had considerably lower relative abundances of non-dominant phyla (Bacteroidetes) than cotton hosts. At the genus level the dominant communities were Buchnera, Acinetobacter, and Arsenophonus. The relative abundance of Buchnera was significantly higher in aphids reared on zucchini than those on cotton, whereas the opposite was observed for Acinetobacter, as well as for some non-dominant communities (Stenotrophomonas, Pseudomons, Flavobacterium, Novosphingobium). Collectively, this study clarifies the dynamic changes of symbiotic bacteria in cotton-specialized aphids reared on zucchini for multiple generations. Among them, Buchnera is crucial for the cotton-specialized aphid to get nutrients during the transfer of the host and has a favorable impact on the colonization of cotton-specialized aphid populations on zucchini hosts. It not only enriches our understanding of the relationship between the bacterial microbiota of aphids and their adaptability to new hosts, zucchini, but also expands the current body of research on the mechanisms underlying the host shifting ability of cotton-specialized aphids.

What Are the Functional Roles of Piwi Proteins and piRNAs in Insects?

Research on Piwi proteins and piRNAs in insects has focused on three experimental models: oogenesis and spermatogenesis in Drosophila melanogaster, the antiviral response in Aedes mosquitoes and the molecular analysis of primary and secondary piRNA biogenesis in Bombyx mori-derived BmN4 cells. Significant unique and complementary information has been acquired and has led to a greater appreciation of the complexity of piRNA biogenesis and Piwi protein function. Studies performed in other insect species are emerging and promise to add to the current state of the art on the roles of piRNAs and Piwi proteins. Although the primary role of the piRNA pathway is genome defense against transposons, particularly in the germline, recent findings also indicate an expansion of its functions. In this review, an extensive overview is presented of the knowledge of the piRNA pathway that so far has accumulated in insects. Following a presentation of the three major models, data from other insects were also discussed. Finally, the mechanisms for the expansion of the function of the piRNA pathway from transposon control to gene regulation were considered.

16S rRNA Gene Sequencing of Six Psyllid Species of the Family Carsidaridae Identified Various Bacteria Including Symbiopectobacterium

Psyllids (Hemiptera: Sternorrhyncha: Psylloidea) are plant sap-sucking insects that are closely associated with various microbes. To obtain a more detailed understanding of the ecological and evolutionary behaviors of microbes in Psylloidea, the bacterial populations of six psyllid species, belonging to the family Carsidaridae, were analyzed using high-throughput amplicon sequencing of the 16S rRNA gene. The majority of the secondary symbionts identified in the present study were gammaproteobacteria, particularly those of the order Enterobacterales, including Arsenophonus and Sodalis, which are lineages found in a wide variety of insect hosts. Additionally, Symbiopectobacterium, another Enterobacterales lineage, which has recently been recognized and increasingly shown to be vertically transmitted and mutualistic in various invertebrates, was identified for the first time in Psylloidea. This lineage is closely related to Pectobacterium spp., which are plant pathogens, but forms a distinct clade exhibiting no pathogenicity to plants. Non-Enterobacterales gammaproteobacteria found in the present study were Acinetobacter, Pseudomonas (both Pseudomonadales), Delftia, Comamonas (both Burkholderiales), and Xanthomonas (Xanthomonadales), a putative plant pathogen. Regarding alphaproteobacteria, three Wolbachia (Rickettsiales) lineages belonging to supergroup B, the major group in insect lineages, were detected in four psyllid species. In addition, a Wolbachia lineage of supergroup O, a minor group recently found for the first time in Psylloidea, was detected in one psyllid species. These results suggest the pervasive transfer of bacterial symbionts among animals and plants, providing deeper insights into the evolution of the interactions among these organisms.

Aphid facultative symbionts confer no protection against the fungal entomopathogen Batkoa apiculata

Fungi in the family Entomophthoraceae are prevalent pathogens of aphids. Facultative symbiotic bacteria harbored by aphids, including Spiroplasma sp. and Regiella insecticola, have been shown to make their hosts more resistant to infection with the fungal pathogen Pandora neoaphidis. How far this protection extends against other species of fungi in the family Entomophthoraceae is unknown. Here we isolated a strain of the fungal pathogen Batkoa apiculata infecting a natural population of pea aphids (Acyrthosiphon pisum) and confirmed its identity by sequencing the 28S rRNA gene. We then infected a panel of aphids each harboring a different species or strain of endosymbiotic bacteria to test whether aphid symbionts protect against B. apiculata. We found no evidence of symbiont-mediated protection against this pathogen, and our data suggest that some symbionts make aphids more susceptible to infection. This finding is relevant to our understanding of this important model of host-microbe interactions, and we discuss our results in the context of aphid-microbe ecological and evolutionary dynamics.

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.

Effects of abiotic stresses on the expression of chitinase-like genes in Acyrthosiphon pisum

Insect chitinases play a crucial part to digest chitin in the exoskeleton during the molting process. However, research on insect chitinase related to the environmental stress response is very limited. This study was the first conducted to expression analysis of chitinase- related genes in A. pisum under abiotic stresses. Here, we identified five chitinase-like proteins (ApIDGF, ApCht3, ApCht7, ApCht10 and ApENGase), and clustered them into five groups (group II, III, V, Ⅹ, and ENGase). Developmental expression analysis revealed that the five A. pisum chitinase-related genes were expressed at whole developmental stages with different relative expression patterns. When aphids were exposed to various abiotic stresses including temperature, insecticide and the stress 20-hydroxyecdysone (20E), all five chitinase genes were differentially expressed in A. pisum. The results showed that insecticide such as imidacloprid down-regulated the expression of these five Cht-related genes. Analysis of temperature stress of A. pisum chitinase suggested that ApCht7 expression was high at 10°C, which demonstrates its important role in pea aphids under low temperature. Conversely, ApCht10 was more active under high temperature stress, as it was significantly up-regulated at 30°C. Besides, 20E enhanced ApCht3 and ApCht10 expression in A. pisum, but reduced ApCht7 expression. These findings provide basic information and insights for the study of the role of these genes under abiotic stress, which advances our knowledge in the management of pea aphids under multiple stresses.

Bacteriocyte plasticity in pea aphids facing amino acid stress or starvation during development

An important contributing factor to the evolutionary success of insects is nutritional association with microbial symbionts, which provide the host insects with nutrients lacking in their unbalanced diets. These symbionts are often compartmentalized in specialized cells of the host, the bacteriocytes. Even though bacteriocytes were first described more than a century ago, few studies have explored their dynamics throughout the insect life cycle and in response to environmental stressors. Here, we use the Buchnera aphidicola/pea aphid symbiotic system to study how bacteriocytes are regulated in response to nutritional stress throughout aphid development. Using artificial diets, we analyzed the effects of depletion or excess of phenylalanine or leucine, two amino acids essential for aphid growth and whose biosynthetic pathways are shared between the host and the symbiont. Bacteriocytes responded dynamically to those treatments, while other tissues showed no obvious morphological change. Amino acid depletion resulted in an increase in bacteriocyte numbers, with the extent of the increase depending on the amino acid, while excess either caused a decrease (for leucine) or an increase (for phenylalanine). Only a limited impact on survival and fecundity was observed, suggesting that the adjustment in bacteriocyte (and symbiont) numbers is sufficient to withstand these nutritional challenges. We also studied the impact of more extreme conditions by exposing aphids to a 24 h starvation period at the beginning of nymphal development. This led to a dramatic drop in aphid survival and fecundity and a significant developmental delay. Again, bacteriocytes responded dynamically, with a considerable decrease in number and size, correlated with a decrease in the number of symbionts, which were prematurely degraded by the lysosomal system. This study shows how bacteriocyte dynamics is integrated in the physiology of insects and highlights the high plasticity of these cells.

Coevolution of Metabolic Pathways in Blattodea and Their Blattabacterium Endosymbionts, and Comparisons with Other Insect-Bacteria Symbioses

Many insects harbor bacterial endosymbionts that supply essential nutrients and enable their hosts to thrive on a nutritionally unbalanced diet. Comparisons of the genomes of endosymbionts and their insect hosts have revealed multiple cases of mutually-dependent metabolic pathways that require enzymes encoded in 2 genomes. Complementation of metabolic reactions at the pathway level has been described for hosts feeding on unbalanced diets, such as plant sap. However, the level of collaboration between symbionts and hosts that feed on more variable diets is largely unknown. In this study, we investigated amino acid and vitamin/cofactor biosynthetic pathways in Blattodea, which comprises cockroaches and termites, and their obligate endosymbiont Blattabacterium cuenoti (hereafter Blattabacterium). In contrast to other obligate symbiotic systems, we found no clear evidence of "collaborative pathways" for amino acid biosynthesis in the genomes of these taxa, with the exception of collaborative arginine biosynthesis in 2 taxa, Cryptocercus punctulatus and Mastotermes darwiniensis. Nevertheless, we found that several gaps specific to Blattabacterium in the folate biosynthetic pathway are likely to be complemented by their host. Comparisons with other insects revealed that, with the exception of the arginine biosynthetic pathway, collaborative pathways for essential amino acids are only observed in phloem-sap feeders. These results suggest that the host diet is an important driving factor of metabolic pathway evolution in obligate symbiotic systems. IMPORTANCE The long-term coevolution between insects and their obligate endosymbionts is accompanied by increasing levels of genome integration, sometimes to the point that metabolic pathways require enzymes encoded in two genomes, which we refer to as "collaborative pathways". To date, collaborative pathways have only been reported from sap-feeding insects. Here, we examined metabolic interactions between cockroaches, a group of detritivorous insects, and their obligate endosymbiont, Blattabacterium, and only found evidence of collaborative pathways for arginine biosynthesis. The rarity of collaborative pathways in cockroaches and Blattabacterium contrasts with their prevalence in insect hosts feeding on phloem-sap. Our results suggest that host diet is a factor affecting metabolic integration in obligate symbiotic systems.

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.

Chromosomal-level assembly of Bactericera cockerelli reveals rampant gene family expansions impacting genome structure, function and insect-microbe-plant-interactions

Lineage specific expansions and gene duplications are some of the most important sources of evolutionary novelty in eukaryotes. Although not as prevalent in eukaryotes compared to bacteria, horizontal gene transfer events can also result in key adaptations for insects, especially for those involved in insect-microbe interactions. In this study we assemble the first chromosomal assembly of the psyllid Bactericera cockerelli and reveal that the B. cockerelli genome has experienced significantly more gene expansion events compared to other Hemipteran representatives with fully sequenced genomes. We also reveal that B. cockerelli's genome is the largest psyllid genome (567 Mb) sequenced to date and is ~15% larger than the other two psyllid species genomes sequenced (Pachypsylla venusta and Diaphorina citri). Structurally, B. cockerelli appears to have an additional chromosome compared to the distantly related psyllid species P. venusta due to a previous chromosomal fission or fusion event. The increase in genome size and dynamic nature of the B. cockerelli genome may largely be contributed to the widespread expansion of type I and type II repeat elements that are rampant across all of B. cockerelli's. chromosomes. These repeat elements are distributed near equally in both euchromatic and heterochromatic regions. Furthermore, significant gene family expansions and gene duplications were uncovered for genes that are expected to be important in its adaptation to insect-plant and microbe interactions, which include transcription factors, proteases, odorant receptors, and horizontally transferred genes that are involved in the nutritional symbioses with their long-term nutritional endosymbiont Carsonella.

Germline specification and axis determination in viviparous and oviparous pea aphids: conserved and divergent features

Aphids are hemimetabolous insects that undergo incomplete metamorphosis without pupation. The annual life cycle of most aphids includes both an asexual (viviparous) and a sexual (oviparous) phase. Sexual reproduction only occurs once per year and is followed by many generations of asexual reproduction, during which aphids propagate exponentially with telescopic development. Here, we discuss the potential links between viviparous embryogenesis and derived developmental features in the pea aphid Acyrthosiphon pisum, particularly focusing on germline specification and axis determination, both of which are key events of early development in insects. We also discuss potential evolutionary paths through which both viviparous and oviparous females might have come to utilize maternal germ plasm to drive germline specification. This developmental strategy, as defined by germline markers, has not been reported in other hemimetabolous insects. In viviparous females, furthermore, we discuss whether molecules that in other insects characterize germ plasm, like Vasa, also participate in posterior determination and how the anterior localization of the hunchback orthologue Ap-hb establishes the anterior-posterior axis. We propose that the linked chain of developing oocytes and embryos within each ovariole and the special morphology of early embryos might have driven the formation of evolutionary novelties in germline specification and axis determination in the viviparous aphids. Moreover, based upon the finding that the endosymbiont Buchnera aphidicola is closely associated with germ cells throughout embryogenesis, we propose presumptive roles for B. aphidicola in aphid development, discussing how it might regulate germline migration in both reproductive modes of pea aphids. In summary, we expect that this review will shed light on viviparous as well as oviparous development in aphids.

Downregulation of NADPH-cytochrome P450 reductase via RNA interference increases the susceptibility of Acyrthosiphon pisum to desiccation and insecticides

Nicotinamide adenine dinucleotide phosphate (NADPH)-cytochrome P450 reductase (CPR) is involved in the metabolism of endogenous and exogenous substances, and detoxification of insecticides. RNA interference (RNAi) of CPR in certain insects causes developmental defects and enhanced susceptibility to insecticides. However, the CPR of Acyrthosiphon pisum has not been characterized, and its function is still not understood. In this study, we investigated the biochemical functions of A. pisum CPR (ApCPR). ApCPR was found to be transcribed in all developmental stages and was abundant in the embryo stage, and in the gut, head, and abdominal cuticle. After optimizing the dose and silencing duration of RNAi for downregulating ApCPR, we found that ApCPR suppression resulted in a significant decrease in the production of cuticular and internal hydrocarbon contents, and of cuticular waxy coatings. Deficiency in cuticular hydrocarbons (CHCs) decreased the survival rate of A. pisum under desiccation stress and increased its susceptibility to contact insecticides. Moreover, desiccation stress induced a significant increase in ApCPR mRNA levels. We further confirmed that ApCPR participates in CHC production. These results indicate that ApCPR modulates CHC production, desiccation tolerance, and insecticide susceptibility in A. pisum, and presents a novel target for pest control.

Ploidy dynamics in aphid host cells harboring bacterial symbionts

Aphids have evolved bacteriocytes or symbiotic host cells that harbor the obligate mutualistic bacterium Buchnera aphidicola. Because of the large cell size (approximately 100 μm in diameter) of bacteriocytes and their pivotal role in nutritional symbiosis, researchers have considered that these cells are highly polyploid and assumed that bacteriocyte polyploidy may be essential for the symbiotic relationship between the aphid and the bacterium. However, little is known about the ploidy levels and dynamics of aphid bacteriocytes. Here, we quantitatively analyzed the ploidy levels in the bacteriocytes of the pea-aphid Acyrthosiphon pisum. Image-based fluorometry revealed the hyper polyploidy of the bacteriocytes ranging from 16- to 256-ploidy throughout the lifecycle. Bacteriocytes of adult parthenogenetic viviparous females were ranged between 64 and 128C DNA levels, while those of sexual morphs (oviparous females and males) were comprised of 64C, and 32–64C cells, respectively. During post-embryonic development of viviparous females, the ploidy level of bacteriocytes increased substantially, from 16 to 32C at birth to 128–256C in actively reproducing adults. These results suggest that the ploidy levels are dynamically regulated among phenotypes and during development. Our comprehensive and quantitative data provides a foundation for future studies to understand the functional roles and biological significance of the polyploidy of insect bacteriocytes.

Genome of the parasitoid wasp Cotesia chilonis sheds light on amino acid resource exploitation

Background

A fundamental feature of parasitism is the nutritional exploitation of host organisms by their parasites. Parasitoid wasps lay eggs on arthropod hosts, exploiting them for nutrition to support larval development by using diverse effectors aimed at regulating host metabolism. However, the genetic components and molecular mechanisms at the basis of such exploitation, especially the utilization of host amino acid resources, remain largely unknown. To address this question, here, we present a chromosome-level genome assembly of the parasitoid wasp Cotesia chilonis and reconstruct its amino acid biosynthetic pathway.

Results

Analyses of the amino acid synthetic pathway indicate that C. chilonis lost the ability to synthesize ten amino acids, which was confirmed by feeding experiments with amino acid-depleted media. Of the ten pathways, nine are known to have been lost in the common ancestor of animals. We find that the ability to synthesize arginine was also lost in C. chilonis because of the absence of two key genes in the arginine synthesis pathway. Further analyses of the genomes of 72 arthropods species show that the loss of arginine synthesis is common in arthropods. Metabolomic analyses by UPLC-MS/MS reveal that the temporal concentrations of arginine, serine, tyrosine, and alanine are significantly higher in host (Chilo suppressalis) hemolymph at 3 days after parasitism, whereas the temporal levels of 5-hydroxylysine, glutamic acid, methionine, and lysine are significantly lower. We sequence the transcriptomes of a parasitized host and non-parasitized control. Differential gene expression analyses using these transcriptomes indicate that parasitoid wasps inhibit amino acid utilization and activate protein degradation in the host, likely resulting in the increase of amino acid content in host hemolymph.

Conclusions

We sequenced the genome of a parasitoid wasp, C. chilonis, and revealed the features of trait loss in amino acid biosynthesis. Our work provides new insights into amino acid exploitation by parasitoid wasps, and this knowledge can specifically be used to design parasitoid artificial diets that potentially benefit mass rearing of parasitoids for pest control.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01313-3.

Chromosome‐level genome assembly for the horned‐gall aphid provides insights into interactions between gall‐making insect and its host plant

The aphid Schlechtendalia chinensis is an economically important insect that can induce horned galls, which are valuable for the medicinal and chemical industries. Up to now, more than twenty aphid genomes have been reported. Most of the sequenced genomes are derived from free‐living aphids. Here, we generated a high‐quality genome assembly from a galling aphid. The final genome assembly is 271.52 Mb, representing one of the smallest sequenced genomes of aphids. The genome assembly is based on contig and scaffold N50 values of the genome sequence are 3.77 Mb and 20.41 Mb, respectively. Nine‐seven percent of the assembled sequences was anchored onto 13 chromosomes. Based on BUSCO analysis, the assembly involved 96.9% of conserved arthropod and 98.5% of the conserved Hemiptera single‐copy orthologous genes. A total of 14,089 protein‐coding genes were predicted. Phylogenetic analysis revealed that S. chinensis diverged from the common ancestor of Eriosoma lanigerum approximately 57 million years ago (MYA). In addition, 35 genes encoding salivary gland proteins showed differentially when S. chinensis forms a gall, suggesting they have potential roles in gall formation and plant defense suppression. Taken together, this high‐quality S. chinensis genome assembly and annotation provide a solid genetic foundation for future research to reveal the mechanism of gall formation and to explore the interaction between aphids and their host plants.

Characterization of transposable elements within the Bemisia tabaci species complex

Background

Whiteflies are agricultural pests that cause negative impacts globally to crop yields resulting at times in severe economic losses and food insecurity. The Bemisia tabaci whitefly species complex is the most damaging in terms of its broad crop host range and its ability to serve as vector for over 400 plant viruses. Genomes of whiteflies belonging to this species complex have provided valuable genomic data; however, transposable elements (TEs) within these genomes remain unexplored. This study provides the first accurate characterization of TE content within the B. tabaci species complex.

Results

This study identified that an average of 40.61% of the genomes of three whitefly species (MEAM1, MEDQ, and SSA-ECA) consists of TEs. The majority of the TEs identified were DNA transposons (22.85% average) while SINEs (0.14% average) were the least represented. This study also compared the TE content of the three whitefly genomes with three other hemipteran genomes and found significantly more DNA transposons and less LINEs in the whitefly genomes. A total of 63 TE superfamilies were identified to be present across the three whitefly species (39 DNA transposons, six LTR, 16 LINE, and two SINE). The sequences of the identified TEs were clustered which generated 5766 TE clusters. A total of 2707 clusters were identified as uniquely found within the whitefly genomes while none of the generated clusters were from both whitefly and non-whitefly TE sequences.

This study is the first to characterize TEs found within different B. tabaci species and has created a standardized annotation workflow that could be used to analyze future whitefly genomes.

Conclusion

This study is the first to characterize the landscape of TEs within the B. tabaci whitefly species complex. The characterization of these elements within the three whitefly genomes shows that TEs occupy significant portions of B. tabaci genomes, with DNA transposons representing the vast majority. This study also identified TE superfamilies and clusters of TE sequences of potential interest, providing essential information, and a framework for future TE studies within this species complex.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13100-022-00270-6.

Annotation of Hox cluster and Hox cofactor genes in the Asian citrus psyllid, Diaphorina citri, reveals novel features

Hox genes and their cofactors are essential developmental genes specifying regional identity in animals. Hox genes have a conserved arrangement in clusters in the same order in which they specify identity along the anterior-posterior axis. A few insect species have breaks in the cluster, but these are exceptions. We annotated the 10 Hox genes of the Asian citrus psyllid Diaphorina citri, and found a split in its Hox cluster between the Deformed and Sex combs reduced genes - the first time a break at this position has been observed in an insect Hox cluster. We also annotated D. citri orthologs of the Hox cofactor genes homothorax, PKNOX and extradenticle and found an additional copy of extradenticle in D. citri that appears to be a retrogene. Expression data and sequence conservation suggest that the extradenticle retrogene may have retained the original extradenticle function and allowed divergence of the parental extradenticle gene.

Cholecystokinin/sulfakinin peptide signaling: conserved roles at the intersection between feeding, mating and aggression

Neuropeptides are the most diverse messenger molecules in metazoans and are involved in regulation of daily physiology and a wide array of behaviors. Some neuropeptides and their cognate receptors are structurally and functionally well conserved over evolution in bilaterian animals. Among these are peptides related to gastrin and cholecystokinin (CCK). In mammals, CCK is produced by intestinal endocrine cells and brain neurons, and regulates gall bladder contractions, pancreatic enzyme secretion, gut functions, satiety and food intake. Additionally, CCK plays important roles in neuromodulation in several brain circuits that regulate reward, anxiety, aggression and sexual behavior. In invertebrates, CCK-type peptides (sulfakinins, SKs) are, with a few exceptions, produced by brain neurons only. Common among invertebrates is that SKs mediate satiety and regulate food ingestion by a variety of mechanisms. Also regulation of secretion of digestive enzymes has been reported. Studies of the genetically tractable fly Drosophila have advanced our understanding of SK signaling mechanisms in regulation of satiety and feeding, but also in gustatory sensitivity, locomotor activity, aggression and reproductive behavior. A set of eight SK-expressing brain neurons plays important roles in regulation of these competing behaviors. In males, they integrate internal state and external stimuli to diminish sex drive and increase aggression. The same neurons also diminish sugar gustation, induce satiety and reduce feeding. Although several functional roles of CCK/SK signaling appear conserved between Drosophila and mammals, available data suggest that the underlying mechanisms differ.

Genome sequence of the English grain aphid, Sitobion avenae and its endosymbiont Buchnera aphidicola

The English grain aphid, Sitobion avenae, is a major agricultural pest of wheat, barley and oats, and one of the principal vectors of barley yellow dwarf virus leading to significant reductions in grain yield, annually. Emerging resistance to and increasing regulation of insecticides has resulted in limited options for their control. Using PacBio HiFi data, we have produced a high-quality draft assembly of the S. avenae genome; generating a primary assembly with a total assembly size of 475.7 Mb, and an alternate assembly with a total assembly size of 430.8 Mb. Our primary assembly was highly contiguous with only 326 contigs and a contig N50 of 15.95 Mb. Assembly completeness was estimated at 97.7% using BUSCO analysis and 31,007 and 29,037 protein-coding genes were predicted from the primary and alternate assemblies, respectively. This assembly, which is to our knowledge the first for an insecticide resistant clonal lineage of English grain aphid, will provide novel insight into the molecular and mechanistic determinants of resistance and will facilitate future research into mechanisms of viral transmission and aphid behavior.

Annotation of glycolysis, gluconeogenesis, and trehaloneogenesis pathways provide insight into carbohydrate metabolism in the Asian citrus psyllid

Citrus greening disease is caused by the pathogen Candidatus Liberibacter asiaticus and transmitted by the Asian citrus psyllid, Diaphorina citri. No curative treatment or significant prevention mechanism exists for this disease, which causes economic losses from reduced citrus production. A high-quality genome of D. citri is being manually annotated to provide accurate gene models to identify novel control targets and increase understanding of this pest. Here, we annotated 25 D. citri genes involved in glycolysis and gluconeogenesis, and seven in trehaloneogenesis. Comparative analysis showed that glycolysis genes in D. citri are highly conserved but copy numbers vary. Analysis of expression levels revealed upregulation of several enzymes in the glycolysis pathway in the thorax, consistent with the primary use of glucose by thoracic flight muscles. Manually annotating these core metabolic pathways provides accurate genomic foundation for developing gene-targeting therapeutics to control D. citri.

Genomic identification, annotation, and comparative analysis of Vacuolar-type ATP synthase subunits in Diaphorina citri

The hemipteran insect Diaphorina citri, or Asian citrus psyllid, is a vector for Candidatus Liberibacter asiaticus (CLas), the bacterium causing citrus greening disease, or Huanglongbing (HLB). Millions of citrus trees have been destroyed, and every grove in Florida, USA, has been directly affected by this disease. In eukaryotes, vacuolar-type ATP synthase (V-ATPase) is an abundant heterodimeric enzyme that serves the cell with essential compartment acidification through the active processes that transport protons across the membrane. Fifteen putative V-ATPase genes in the D. citri genome were manually curated. Comparative genomic analysis revealed that D. citri V-ATPase subunits share domains and motifs with other insects, including the V-ATPase-A superfamily domain. Phylogenetic analysis separates D. citri V-ATPase subunits into expected clades with orthologous sequences. Annotation of the D. citri genome is a critical step towards developing directed pest management strategies to reduce the spread of HLB throughout the citrus industry.

TmIKKε Is Required to Confer Protection Against Gram-Negative Bacteria, E. coli by the Regulation of Antimicrobial Peptide Production in the Tenebrio molitor Fat Body

The inhibitor of nuclear factor-kappa B (NF-κB) kinase (IKK) is the core regulator of the NF-κB pathway against pathogenic invasion in vertebrates or invertebrates. IKKβ, -ε and -γ have pivotal roles in the Toll and immune deficiency (IMD) pathways. In this study, a homolog of IKKε (TmIKKε) was identified from Tenebrio molitor RNA sequence database and functionally characterized for its role in regulating immune signaling pathways in insects. The TmIKKε gene is characterized by two exons and one intron comprising an open reading frame (ORF) of 2,196 bp that putatively encodes a polypeptide of 731 amino acid residues. TmIKKε contains a serine/threonine protein kinases catalytic domain. Phylogenetic analysis established the close homology of TmIKKε to Tribolium castaneum IKKε (TcIKKε) and its proximity with other IKK-related kinases. The expression of TmIKKε mRNA was elevated in the gut, integument, and hemocytes of the last-instar larva and the fat body, Malpighian tubules, and testis of 5-day-old adults. TmIKKε expression was significantly induced by Escherichia coli, Staphylococcus aureus, and Candida albicans challenge in whole larvae and tissues, such as hemocytes, gut, and fat body. The knockdown of the TmIKKε messenger RNA (mRNA) expression significantly reduced the survival of the larvae against microbial challenges. Further, we investigated the induction patterns of 14 T. molitor antimicrobial peptides (AMPs) genes in TmIKKε gene-silencing model after microbial challenges. While in hemocytes, the transcriptional regulation of most AMPs was negatively regulated in the gut and fat body tissue of T. molitor, AMPs, such as TmTenecin 1, TmTenecin 4, TmDefensin, TmColeoptericin A, TmColeoptericin B, TmAttacin 1a, and TmAttacin 2, were positively regulated in TmIKKε-silenced individuals after microbial challenge. Collectively, the results implicate TmIKKε as an important factor in antimicrobial innate immune responses in T. molitor.

Dissecting cricket genomes for the advancement of entomology and entomophagy

Significant advances in biophysical methods such as next-generation sequencing technologies have now opened the way to conduct evolutionary and applied research based on the genomic information of greatly diverse insects. Crickets belonging to Orthoptera (Insecta: Polyneoptera), one of the most flourishing groups of insects, have contributed to the development of multiple scientific fields including developmental biology and neuroscience and have been attractive targets in evolutionary ecology for their diverse ecological niches. In addition, crickets have recently gained recognition as food and feed. However, the genomic information underlying their biological basis and application research toward breeding is currently underrepresented. In this review, we summarize the progress of genomics of crickets. First, we outline the phylogenetic position of crickets in insects and then introduce recent studies on cricket genomics and transcriptomics in a variety of fields. Furthermore, we present findings from our analysis of polyneopteran genomes, with a particular focus on their large genome sizes, chromosome number, and repetitive sequences. Finally, how the cricket genome can be beneficial to the food industry is discussed. This review is expected to enhance greater recognition of how important the cricket genomes are to the multiple biological fields and how basic research based on cricket genome information can contribute to tackling global food security.

A scaffold-level genome assembly of a minute pirate bug, Orius laevigatus (Hemiptera: Anthocoridae), and a comparative analysis of insecticide resistance-related gene families with hemipteran crop pests

BACKGROUND

Orius laevigatus, a minute pirate bug, is a highly effective beneficial predator of crop pests including aphids, spider mites and thrips in integrated pest management (IPM) programmes. No genomic information is currently available for O. laevigatus, as is the case for the majority of beneficial predators which feed on crop pests. In contrast, genomic information for crop pests is far more readily available. The lack of publicly available genomes for beneficial predators to date has limited our ability to perform comparative analyses of genes encoding potential insecticide resistance mechanisms between crop pests and their predators. These mechanisms include several gene/protein families including cytochrome P450s (P450s), ATP binding cassette transporters (ABCs), glutathione S-transferases (GSTs), UDP-glucosyltransferases (UGTs) and carboxyl/cholinesterases (CCEs).

METHODS AND FINDINGS

In this study, a high-quality scaffold level de novo genome assembly for O. laevigatus has been generated using a hybrid approach with PacBio long-read and Illumina short-read data. The final assembly achieved a scaffold N50 of 125,649 bp and a total genome size of 150.98 Mb. The genome assembly achieved a level of completeness of 93.6% using a set of 1658 core insect genes present as full-length genes. Genome annotation identified 15,102 protein-coding genes - 87% of which were assigned a putative function. Comparative analyses revealed gene expansions of sigma class GSTs and CYP3 P450s. Conversely the UGT gene family showed limited expansion. Differences were seen in the distributions of resistance-associated gene families at the subfamily level between O. laevigatus and some of its targeted crop pests. A target site mutation in ryanodine receptors (I4790M, PxRyR) which has strong links to diamide resistance in crop pests and had previously only been identified in lepidopteran species was found to also be present in hemipteran species, including O. laevigatus.

CONCLUSION AND SIGNIFICANCE

This assembly is the first published genome for the Anthocoridae family and will serve as a useful resource for further research into target-site selectivity issues and potential resistance mechanisms in beneficial predators. Furthermore, the expansion of gene families often linked to insecticide resistance may be an indicator of the capacity of this predator to detoxify selective insecticides. These findings could be exploited by targeted pesticide screens and functional studies to increase effectiveness of IPM strategies, which aim to increase crop yields by sustainably, environmentally-friendly and effectively control pests without impacting beneficial predator populations.

Microbiome analyses of 12 psyllid species of the family Psyllidae identified various bacteria including Fukatsuia and Serratia symbiotica, known as secondary symbionts of aphids

Background

Psyllids (Hemiptera: Psylloidea) comprise a group of plant sap-sucking insects that includes important agricultural pests. They have close associations not only with plant pathogens, but also with various microbes, including obligate mutualists and facultative symbionts. Recent studies are revealing that interactions among such bacterial populations are important for psyllid biology and host plant pathology. In the present study, to obtain further insight into the ecological and evolutionary behaviors of bacteria in Psylloidea, we analyzed the microbiomes of 12 psyllid species belonging to the family Psyllidae (11 from Psyllinae and one from Macrocorsinae), using high-throughput amplicon sequencing of the 16S rRNA gene.

Results

The analysis showed that all 12 psyllids have the primary symbiont, Candidatus Carsonella ruddii (Gammaproteobacteria: Oceanospirillales), and at least one secondary symbiont. The majority of the secondary symbionts were gammaproteobacteria, especially those of the family Enterobacteriaceae (order: Enterobacteriales). Among them, symbionts belonging to “endosymbionts3”, which is a genus-level monophyletic group assigned by the SILVA rRNA database, were the most prevalent and were found in 9 of 11 Psyllinae species. Ca. Fukatsuia symbiotica and Serratia symbiotica, which were recognized only as secondary symbionts of aphids, were also identified. In addition to other Enterobacteriaceae bacteria, including Arsenophonus, Sodalis, and “endosymbionts2”, which is another genus-level clade, Pseudomonas (Pseudomonadales: Pseudomonadaceae) and Diplorickettsia (Diplorickettsiales: Diplorickettsiaceae) were identified. Regarding Alphaproteobacteria, the potential plant pathogen Ca. Liberibacter europaeus (Rhizobiales: Rhizobiaceae) was detected for the first time in Anomoneura mori (Psyllinae), a mulberry pest. Wolbachia (Rickettsiales: Anaplasmataceae) and Rickettsia (Rickettsiales: Rickettsiaceae), plausible host reproduction manipulators that are potential tools to control pest insects, were also detected.

Conclusions

The present study identified various bacterial symbionts including previously unexpected lineages in psyllids, suggesting considerable interspecific transfer of arthropod symbionts. The findings provide deeper insights into the evolution of interactions among insects, bacteria, and plants, which may be exploited to facilitate the control of pest psyllids in the future.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12866-021-02429-2.