Next Article in Journal
Anticancer Effects of Mitoquinone via Cell Cycle Arrest and Apoptosis in Canine Mammary Gland Tumor Cells
Next Article in Special Issue
Challenges of Robust RNAi-Mediated Gene Silencing in Aedes Mosquitoes
Previous Article in Journal
The Development of CDC25A-Derived Phosphoseryl Peptides That Bind 14-3-3ε with High Affinities
Previous Article in Special Issue
Physiological Measurements and Transcriptomics Reveal the Fitness Costs of Monochamus saltuarius to Bursaphelenchus xylophilus
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

1
Department of Agricultural Sciences, University of Naples Federico II, 80126 Naples, Italy
2
BAT Center—Interuniversity Center for Studies on Bioinspired Agro-Environmental Technology, University of Naples Federico II, 80126 Naples, Italy
3
Research Centre for Olive, Fruit and Citrus Crops, Council for Agricultural Research and Economics, 81100 Caserta, Italy
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2024, 25(9), 4922; https://doi.org/10.3390/ijms25094922
Submission received: 30 March 2024 / Revised: 25 April 2024 / Accepted: 29 April 2024 / Published: 30 April 2024

Abstract

:
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.

1. Introduction

Due to their constant exposure to a diverse range of pathogens, parasites, and environmental stressors, insects exhibit a multi-faceted immune system [1,2]. In addition to behavioral adaptations and mechanical barriers, such as the exoskeleton, insects have evolved effective innate immunity that acts through both cellular and humoral responses directed against invaders [1,3]. Activation of these responses occurs through the recognition of pathogen-associated molecular patterns (PAMPs) by receptors, known as pattern recognition receptors (PRRs), located on hemocytes (immune cells) and epithelial cells from barrier sites throughout the insect’s body [4,5]. Cellular immune responses include phagocytosis, nodulation, encapsulation, and melanization events mediated by hemocytes [6,7,8]. Humoral responses orchestrated by signaling pathways such as Imd, Toll, Jak/Stat, and JNK lead to the synthesis of various defense enzymes, complement-like proteins, and antimicrobial peptides (AMPs) in response to infection [9,10,11]. Immune responses vary depending on the size and PAMPs of the intruder. Encapsulation and melanization function as defensive strategies against larger intruders, such as the eggs laid by endophagous parasitoids [12], whereas AMPs play a central role in mitigating the impact of pathogenic microorganisms [13,14].
Although the immune system of insects has been extensively studied, the majority of the studies on this subject are focused on holometabolous species, including the fruit fly Drosophila melanogaster Meigen (Diptera: Drosophilidae), honey bees, mosquitoes, beetles, and moths [15,16,17,18,19]. One of the first heterometabolous insects for which the immune system was characterized at both the genomic and functional levels was the pea aphid Acyrthosiphon pisum (Harris) (Hemiptera: Sternorrhyncha; Aphididae) [20,21]. Several studies indicate that this agricultural pest exhibits a reduced immune system [22], which is enhanced by its secondary bacterial symbionts [23]. The pea aphid is widely accepted as a model insect for studying a wide range of biological and physiological processes [24], including those deriving from their interaction with parasites and symbionts [25]. Furthermore, its genetic information [26] has the potential to provide valuable insights into other hemipteran pests that are less extensively studied. Indeed, recent omics studies pointed out that immune system reduction is not restricted to aphids, but is a common feature of several hemipteran species, such as Diaphorina citri (Hemiptera: Liviidae) [27], Plautia stali (Hemiptera: Pentatomidae) [28], and Rodnius prolixus (Hemiptera: Reduviidae) [29].
Alongside closely related hemipterans such as aphids (Aphidomorpha) and whiteflies (Aleyrodomorpha) [30], scale insects (Coccomorpha) are sap-sucking and obligate plant parasites [31,32]. The feeding behavior of scale insects delays plant growth and, in severe infestations, can lead to the death of the entire plant [31]. Indirect damage derives from the production of honeydew, which results in the growth of saprophytic fungi, thereby reducing the rate of plant photosynthesis and causing decline [31]. In addition, some species may act as vectors of pathogenic viruses [33]. This taxonomic group shows remarkable diversity in both external and internal morphology, as well as in reproduction and symbiotic systems, making them a fascinating subject for scientific study [31]. Scale insects, whose name derives from the commonly produced protective covering (“scale”), exhibit peculiar adaptations. They have sexual dimorphism, characterized by ephemeral alate males lacking functional mouthparts and stationary non-winged adult females, which produce a variety of protective waxy secretions [31].
The majority of these insects, owing to their imbalanced diet rich in carbon, engage in obligate symbiotic relationships with different species of bacteria or fungi [34,35,36]. Some species exclusively harbor a single obligate symbiont responsible for the synthesis of essential nutrients [37,38], while others form additional associations with facultative symbiotic organisms [39,40]. Symbiotic microorganisms are housed in specialized cells known as bacteriocytes (or mycetocytes), fat body cells, or the midgut epithelium or are dispersed in the hemolymph [34,41].
Among scale insects, the family Coccidae, commonly known as “soft scale insects”, includes approximately 1180 species worldwide [42]. Some of these are important pests of crops [43] or forest plants [44]. These pests have a diverse complex of natural enemies that can control their populations, including commercially available predatory insects and naturally occurring parasitoids [45]. The most important group of their antagonists comprises chalcidoid wasps (Hymenoptera: Chalcidoidea) [46,47], which mainly belong to the Encyrtidae, Aphelinidae, and Eulophidae families [46]. Ladybirds (Coleoptera: Coccinellidae), particularly Cryptolaemus montrouzieri and Chilocorus sp., are well-known predators of Coccidae and have proven effective in various biological control programs [48,49,50]. Furthermore, a number of entomopathogenic fungi can infect and exert a detrimental effect on scale insects. Besides the best known entomopathogens, such as species of Akanthomyces/Lecanicillium, which are commonly reported to haunt populations of Coccidae [51,52], new species have been characterized in recent years for their pathogenicity to these pests [53,54,55,56]. Moreover, species in certain fungal genera that are frequently reported to establish an endophytic association with plants, such as Fusarium and Cladosporium, could play an ecological role in the containment of scale insects [51,57,58]. However, the most intriguing relationship concerns species of Ophiocordyceps. Also referred to with the anamorphic name Hirsutella, these fungi are commonly regarded as specialized entomopathogens [59]. Notwithstanding, evidence from several independent studies has demonstrated that they develop an intimate symbiotic relationship with soft scales involving their transovarial transmission between generations [35,60,61], which deserves to be examined more in depth in view of possible applications of its disruption in pest control.
As with most sternorrhynchan species, the immunity of soft scale insects has not been thoroughly investigated and remains to be clarified. To date, only a few studies have been conducted on the immune response of Coccidae, mainly aimed at determining the encapsulation response against eggs released by encyrtid parasitoid wasps [62,63,64]. The limited understanding of soft scales’ immunity hinders our comprehension of their interactions with symbionts and pathogens, as well as the potential development of new control strategies based, for example, on the suppression of the immune response through RNA interference, as recently proposed [65,66].
In the present work, we aimed to identify and annotate immune genes in soft scale insects by searching the recently published genome of the Chinese white wax scale insect (Ericerus pela) (Chavannes), as well as the transcriptomes of Coccus sp. and Ceroplastes cirripediformis Comstock (Hemiptera: Coccidae).

2. Results and Discussion

2.1. Overview of Immune Genes’ Annotation

We focused our annotation efforts on a subset of genes involved in the three phases of the insect immune response: recognition, signaling, and response. All annotations are based on the recently completed sequencing of E. pela (colony RIRI-1) [67]. By using protein sequences from A. pisum and D. melanogaster, known to be involved in insect immunity, as queries in BLAST searches, we successfully identified 66 potential immune genes in E. pela. Specifically, we identified 12 genes related to recognition, 35 involved in signaling, and 19 associated with the response to pathogenic microorganisms (Figure 1).

2.2. Annotation of Recognition Genes

Our analysis pointed out the occurrence of 12 genes in Coccidae with significant matches with A. pisum and Drosophila genes involved in recognition (Table 1). As occurs in aphids, Coccidae species lack peptidoglycan recognition proteins (PGRPs), class C scavengers, and Nimrod and eater receptors.

2.2.1. Peptidoglycan Receptor Proteins

Peptidoglycans are essential cell wall components of almost all bacteria, which are recognized by the immune system through pathogen recognition receptors (PRRs). In insects, several families of pattern recognition molecules that detect peptidoglycans have been identified, and the role of peptidoglycan receptor proteins (PGRPs) in host defense is relatively well-characterized in Drosophila [68]. PGRP-based recognition activates both the Toll and IMD/JNK pathways, leading to proPO activation or the synthesis of antimicrobial peptides [69]. Most insect species investigated possess several PGRP genes that differ both structurally and functionally. For example, Drosophila has 13 PGRP genes encoding 19 proteins, while Anopheles gambiae has 7 PGRP genes encoding 9 proteins [68]. However, like the pea aphid [22], Coccidae appear to have no PGRPs (Table 1).

2.2.2. Gram-Negative Binding Proteins

Gram-negative binding proteins’ (GNBPs) architectures consist of a carbohydrate-binding module (CBM) at the N-terminus and a glucanase-like domain (Glu) in the C-terminus [70]. The CBM interacts with microbial polysaccharides, while the Glu domain interacts with downstream proteases, thereby initiating immune pathways [71]. GNBPs recognize both bacterial and fungal pathogens, resulting in the activation of immune signaling pathways in insects [72]. Specifically, in Drosophila, GNBP1 and peptidoglycan-recognition protein-SA (PGRP-SA) collaboratively activate the Toll pathway in response to Gram-positive bacterial infections [73], whereas GNBP3 is essential for Toll pathway activation in response to fungal infections [74].
A study based on gene knockdown revealed that the two GNBPs predicted in the genome database of pea aphids [22] are involved in the antibacterial response in the pea aphid, likely acting as PRRs in the prophenoloxidase pathway [75]. Our analysis identified a single gene encoded in the transcriptome and the genome of C. cirripediformis and E. pela, respectively, while two different genes seem to occur in the transcriptome of Coccus sp. (Figure 2).

2.2.3. Lectins

Lectins, a diverse group of sugar-binding proteins, are integral to the immune response of several insect species. They are known for their broad spectrum of pathogen binding and involvement in various immune processes such as opsonization, melanization, antibacterial peptide synthesis, encapsulation, and direct killing of bacteria [76]. Drosophila c-type lectins (CTLs) have been implicated in facilitating the encapsulation of parasitoid invaders by marking surfaces for hemocyte recruitment [77]. Interestingly, as in A. pisum, in Coccidae, no homologs of D. melanogaster DL1 (AAF53793.1) have been found (Table 1).
Our phylogenetic reconstruction (Figure 3) revealed that the two CTLs identified in Coccidae are more closely related to DL2 (NP_001014489.1) than to DL3 (NP_001014490.1). This suggests that the phylogeny of CTLs is characterized by species-specific contraction and expansion events, influenced by factors such as environmental pressure, pathogen interactions, and microbiota [78].
Galectins, another widely distributed group of lectins [79], are upregulated in mosquitoes in response to both bacterial and malaria parasite infection [80]. Insect galectins are thought to be involved in pathogen recognition, agglutination, and phagocytosis [79,81]. Genome-wide analyses have revealed variation in galectin transcripts across insect species, with 5 in D. melanogaster, 8 in A. gambiae, 12 in Aedes aegypti [82], 4 in many Lepidoptera species [83], and 1 in aphids. In contrast, Coccidae lack galectin putative homologs (Table 1).

2.2.4. Thioester-Containing Proteins

Thioester-containing proteins (TEPs) are a family of proteins structurally related to vertebrate complement proteins, including an intramolecular β-cysteinyl-γ-glutamyl thioester bond [84,85]. As observed in complement proteins, some TEPs are involved in the opsonization of microbes and pathogens, ‘marking’ them for phagocytosis, melanization, and the formation of lytic complexes [86,87]. Due to their involvement in microbe recognition, TEPs can be classified as PRRs [85].
As in aphids, Coccidae omics data showed the presence of one TEP gene encoding two isoforms, except in Coccus sp., where only one isoform was found (Figure 4). The most closely related proteins are the two isoforms encoded by the only TEP ortholog in the A. pisum genome. Indeed, in contrast to what is reported in [22], referring to an old annotation, only one TEP ortholog was identified by our analysis in the A. pisum genome using four Drosophila homologs (TepI, TepII, TepIII, and TepIV) as the query (Table 1).

2.2.5. Class C Scavenger and Nimrod Receptors

As observed in pea aphids, both class C scavenger and Nimrod receptors are absent in Coccidae (Table 1). Class C scavenger receptors, which have been identified only in Drosophila, exhibit a broad affinity toward both Gram-positive and Gram-negative bacteria [88].
The Nimrod family of proteins is characterized by the presence of epidermal growth factor (EGF)-like domains, also called ‘NIM repeats’ [89]. Several members of the Nimrod superfamily appear to function as receptors in phagocytosis and bacterial binding [90,91]. NimC1 and Eater, two EGF-like repeat Nimrod surface receptors specifically expressed in hemocytes, synergistically contribute to bacterial phagocytosis [92] and are both absent in Coccidae.

2.3. Annotation of Signaling Pathways

Our analysis revealed the occurrence of 35 genes in Coccidae with significant matches with genes of Drosophila and A. pisum involved in signaling. Coccidae lack MyD88, TNF-receptor-associated factor 3, and cactus genes, belonging to the Toll pathway, as well as several members of IMD signaling pathway, which are present in Drosophila and A. pisum genomes (Table 2).

2.3.1. The Toll Signaling Pathway

The Toll pathway in Drosophila functions in both development and innate immunity. Deletion of its component genes increases the susceptibility to various pathogens, including Gram-positive bacteria, fungal pathogens, some Gram-negative bacteria, and viruses [93]. Moreover, upregulation of Toll pathway components occurs in response to parasitoid wasp invasion [94]. The Toll pathway appears to be intact in Coccidae, except the MyD88 adaptor and the inhibitor molecule cactus (a homolog of IkB) (Table 2), which are instead present in A. pisum. We found convincing matches for genes encoding the extracellular cytokine spätzle, the transmembrane receptor Toll (Figure 4), the tube adaptor, the kinase pelle, cactin, pellino, Traf, and the transactivator dorsal (Table 2). Coccidae seem to have multiple spätzles, putative homologs of Drosophila spätzles 1, 2, 3, 4, and 6 (Table 2), for which a phylogenetic reconstruction was not possible due to high divergence between the different spätzle subfamilies [95].
Coccidae also have multiple genes encoding Toll receptors (Figure 5), which function as transmembrane receptors in both mammals and insects. While nine single-copy Toll genes have been identified in D. melanogaster (Toll1 to Toll9), it seems that Coccidae, like other insects, lack some of these genes, but have multiple isoforms of others. Notably, no Toll6 and Toll2/7 homologs have been found in the C. cirripediformis transcriptome, while Coccus sp. lacks Toll10 homologs (Figure 5).
In other organisms, some Toll subfamilies are involved in immune function, while others function in developmental processes [96]. However, an accurate homology-based approach including different species is essential for understanding Toll functions in Coccidae.

2.3.2. The JAK/STAT Signaling Pathway

In Drosophila, the JAK/STAT pathway, similar to the Toll pathway, plays roles in both development and immunity. Despite being the least understood of the core insect immune pathways, it appears to induce hemocyte overproliferation and antiviral responses [97]. Additionally, changes in gene expression observed after parasitoid wasp invasion of Drosophila larvae indicate the involvement of the JAK/STAT pathway in the response to parasitoids [98].
As in A. pisum, Coccidae have homologs of all core JAK/STAT genes, including genes encoding the cytokine receptor domeless (Figure 6), JAK tyrosine kinase (also known as Hopscotch), and the STAT92E transcription factor (Table 2). However, no homologs were found for upd (unpaired), considered a key ligand in Drosophila JAK/STAT induction. This ligand is also missing in other insects (e.g., A. mellifera) [99]. The presence of the core JAK/STAT pathway members (Table 2) suggests that JAK/STAT remains functional in Coccidae and is triggered by a currently unrecognized ligand.

2.3.3. IMD and JNK Signaling Pathways

The IMD pathway is critical for fighting Gram-negative bacteria in Drosophila [93], and IMD pathway member knockouts influence the susceptibility to some Gram-positive bacteria and fungi as well [100]. As in A. pisum, Coccidae appear to be missing many crucial components of the IMD signaling pathway, such as IMD, dFADD, kenny, and Relish (Rel) (Table 2).
Pea aphids lack genes associated with the IMD pathway but possess orthologs for most components of the JNK pathway. In Drosophila, the JNK pathway is involved in various developmental processes, along with wound healing, and has been suggested to regulate antimicrobial peptide gene expression and cellular immune responses [93]. The genes involved include hep, basket, and JRA. Searches for homologs to the Drosophila kayak (kay) gene found no hits in Coccidae. Considering that this gene is also involved in controlling viral infections [101], understanding its role in the Coccidae family, whose members are known for their capacity to transmit plant viruses [33], may be of importance in the management of vector-borne plant pathogens.
Given that in Drosophila, the IMD pathway activates components of the JNK pathway [93], the presence of JNK but absence of the IMD signaling pathway in Coccidae suggests that an alternative pathway for JNK activation, independent of IMD, involving the inducer Eiger, may occur [22].

2.4. Annotation of Response Genes

Our analysis revealed the occurrence of 19 genes in Coccidae with significant matches with genes of Drosophila and A. pisum involved in the immune response. As occurs in aphids, Coccidae species lack antimicrobial peptides; moreover, they do not possess thaumathins, which are present in the A. pisum genome (Table 3).

2.4.1. Antimicrobial Peptides

Antimicrobial peptides (AMPs) play a key role in the immune response of many organisms, including insects. They are the most widely studied humoral effectors and can be produced constitutively or following induction through particular signaling pathways. In the Drosophila genome, there are currently seven well-characterized families of inducible AMPs, including 21 AMP/AMP-like genes, which play an important role both in counteracting the onset of infections and in maintaining homeostasis with symbiotic microorganisms [102]. The main site of production of these molecules is represented by fat body cells; however, hemocytes and cells of the cuticular epithelium, intestinal epithelium, and reproductive tract are also involved in the synthesis of AMPs. Most AMPs are 15–50 amino acids long and have an amphipathic structure; they are able to alter the permeability of the microbial membrane, generating an alteration of the osmotic balance, with consequent lysis [10]. Currently, several structural families of AMPs from insects are known (defensins, cecropins, drosocins, attacins, diptericins, metchnikowins, and melittins), some of which are peculiar to a particular taxonomic group. For example, drosomycins have only been identified in Drosophila, while gloverins have been specifically described in Lepidoptera; cecropins, and attacins, meanwhile, are known in several insect species, while defensins are widely distributed throughout the animal kingdom [14].
As observed in pea aphids [22], Coccidae are missing many of the antimicrobial peptides common to other insects. Extensive searches for genes encoding AMPs previously identified in Hemiptera species (thanatin, pentatomicin, lugensin, cicadin, cryptonin, pyrrhocoricin, oncocin, hemiptericin) also revealed no hits (Table S1). Even the six thaumatin homologs in the A. pisum genome, which show overall sequence and predicted structural similarities to plant thaumatins [22], are absent in Coccidae (Table 3). Recently, it was suggested that the impact of selection on the innate immune system can act on AMPs, indicating that some AMPs can be deleterious molecules in the absence of microbial challenges, due to their costly production and/or their toxic effects [102]. We speculate that such AMP loss derives from a reduced pathogen pressure in Coccidae and other hemipterans, which feed on the generally microbe-free plant phloem [103].

2.4.2. Lysozyme and Peptidoglycan Degradation

Lysozymes, enzymes responsible for breaking down bacterial cell walls by targeting the polysaccharide component of peptidoglycan, are classified into two classes in insects: the c-type, with muramidase activity, and the i-type, which possess both muramidase and isopeptidase activities [104,105]. As observed in A. pisum [22], Coccidae lack genes for several lysozymes (LysB, LysD, LysE, and LysP), which are highly expressed in the gut of Drosophila and are involved in regulating the microbial composition and in degrading peptidoglycan from dietary bacteria [106].
Only three genes encoding i-type Lys were identified in the genome of E. pela (Table 3) and in the transcriptome of Coccus sp., while there were only two in the transcriptome of C. cirripediformis. The i-type Lys sequences of Coccidae formed three monophyletic groups: two contain an invertebrate (I)-type lysozyme domain profile (LYSOZYME_I) and EF-hand calcium-binding domain, and one contains the LYSOZYME_I domain alone (Figure 7). Notably, C. cirripediformis lacks the lysozyme with the LYSOZYME_I domain alone.
Although numerous insect genes encoding both c-type and i-type lysozymes have been identified through genome and transcriptome analyses, i-type lysozymes have been poorly investigated from a functional perspective [107]. An i-type lysozyme of the beetle Harmonia axyridis was recombinantly expressed in the yeast Pichia pastoris, but the purified protein showed no muramidase and no isopeptidase activity [108]. Transcription and immunofluorescence analysis revealed that this i-type lysozyme is produced in the fat body cells but is not inducible by immune challenge. These findings suggest that i-type lysozymes in insects may have acquired novel and as yet undetermined functions during evolution [108].
One of the defining characteristics of the Hemiptera biology is their mutualistic symbiosis with microorganisms. Symbiotic bacteria and fungi play crucial roles in tasks such as nutrition and defense, often residing within specialized cells (bacteriocytes) or dispersed throughout the hemolymph [35,109].
In the bacteriocytes of A. pisum, two genes thought to be involved in the degradation and recycling of peptidoglycan, LD carboxypeptidase (ldcA) and rare lipoprotein A (rlpA), are expressed at high levels [110]. These genes have been acquired from bacteria of the genus Wolbachia or Rickettsia through horizontal gene transfer [111,112,113].
Our analysis did not yield any significant matches for rlpA, but for ldcA, we found a single hit in the genome of E. pela. However, our phylogenetic reconstruction revealed that the protein is encoded by ldcA clusters for Rickettsia and is poorly related with sequences found in aphid genomes (Figure 8). This suggests that the identified protein is likely encoded by the genome of the bacterial symbiont Rickettsia sp., known to colonize the gut of E. pela [114]. The presence of this sequence in the assembled genome of E. pela is likely due to contamination of the sample with DNA from the gut symbiont. The lack of ldcA acquisition by the E. pela genome is in line with previous studies on symbiotic microorganisms of Coccidae, which, rather than being intracellular bacteria as in aphids, are mostly fungi localized in hemolymph, fat body cells, and ovarioles [36].
Moreover, we assessed the presence of homologs of bacteriocyte-specific cysteine-rich (BCR) peptides, which are small disulfide bond-rich proteins expressed exclusively in aphid bacteriocytes putatively involved in endosymbionts’ control and belonging to a structural class of defensins [115,116]. However, our analysis did not yield any significant matches, supporting the hypothesis that Coccidae lack bacteriocytes and their specific mechanisms of symbiosis mediation.

2.4.3. Chitinases

Chitinases, which are glycosyl hydrolases, are enzymes responsible for breaking down chitin, a polymer of N-acetyl-d-glucosamine and the second most abundant biopolymer worldwide. Chitin serves as a structural component in various biological matrices, such as arthropod exoskeletons and fungal cell walls [117]. In arthropods, chitinases fulfill dual roles: aiding in molting processes and serving as defense mechanisms against parasites like fungi and nematodes [118,119]. These enzymes target chitin by hydrolyzing the 1,4-β-linkages between its constituent glucosamine units. Chitinases and lysozymes belong to the same superfamily of hydrolases, exhibiting similar catalytic activities. In fact, certain chitinases possess lysozyme activity, and vice versa [120].
Chitinase-like proteins of Coccidae are included in the majority of the groups classified on a phylogenetic basis [118], except for groups IV and VII where no Coccidae sequences have been identified. In group II, the only member is from C. cirripediformis, while group I includes only sequences from E. pela and Coccus sp. (Figure 9). Further studies are required to determine the biochemical properties and enzymatic activities of these chitinase-like proteins in Coccidae.

2.4.4. Prophenoloxidase

Phenoloxidase-mediated melanin formation is a characteristic feature accompanying wound clotting, phagocytosis, and encapsulation of pathogens and parasites [122]. In aphids, the inactive enzyme prophenoloxidase (ProPO) is activated by serine proteases to produce phenoloxidase [123,124]. ProPO-activating proteinases, such as phenoloxidase activating factor 2 (paf2) in A. pisum, are upregulated in response to parasitization by parasitic wasps [125]. Coccidae appear to possess more than one ProPO homolog, as observed in A. pisum.
Furthermore, our analysis revealed the presence of two paralogs belonging to the paf2 family in the considered species of Coccidae. The sequences from the three Coccidae species form two monophyletic groups, with the most closely related sequences being their orthologs in the A. pisum genome (Figure 10).

2.4.5. Nitric Oxide Synthase

Nitric oxide is a highly unstable free radical gas that has been shown to be toxic to both parasites and pathogens [126]. Production of nitric oxide is mediated by the enzyme nitric oxide synthase (Nos) and triggers the activation of the Toll/IMD signal pathway [127]. In insects, Nos is upregulated after both parasite [128] and bacterial infection [129]. Like pea aphids, E. pela have one Nos homolog (Table 3).

3. Materials and Methods

Immune gene candidates from A. pisum, identified by [22] and reported in Table 1, Table 2 and Table 3, were used to query the E. pela genome (GenBank: GCA_011428145.1). Most searches utilized the tblastn search to search for hits against the assembled E. pela genome, considering as positive only the hits with an e-value less than 1 × 10−5. For genes absent in the A. pisum genome, genes from D. melanogaster were used to query the E. pela genome.
To identify immune-related genes in Coccidae, we aligned the A. pisum major immune protein sequences on the E. pela genome, using exonerate [130], and we extracted the translated CDS, using getorf [131]. In brief, the protein sequences of A. pisum immune genes were used as inputs in exonerate using the protein2genome model, which allows introns in the alignment, but also allows frameshifts, and exon phase changes when a codon is split by an intron [130]. The resulting sequences were translated using getorf, which finds and outputs the sequences of open reading frames (ORFs) in nucleotide sequences [131].
To reconstruct the immune gene phylogeny of Coccidae, the identified immune-related protein sequences of E. pela were used as a query in tblastn searches against the assembled transcriptome of C. cirripediformis (TSA project accession: GCWZ01) and Coccus sp. (TSA project accession: GCWW01). Putative homologous sequences in other insect species (Apis mellifera, Tribolium castaneum, and Bombyx mori) were identified by sequence similarity searches through BlastP, using D. melanogaster proteins (Table 1, Table 2 and Table 3) as the query versus the non-redundant NCBI database (nr NCBI db). However, using the ldcA protein sequence as a query resulted in no hits in the above-mentioned species. Therefore, a blastP search against the whole nr NCBI db was also performed to reconstruct the ldcA phylogeny. One best hit per query was selected and all the protein sequences were aligned using Muscle 3.8 [132], with default settings.
Alignments were automatically trimmed using Gblocks version 0.91b [133] to avoid comparisons of non-conserved regions present only in a subset of the taxa. The best-fit model of amino acid substitution and phylogenetic reconstruction was generated using RAxML 8.2.12 [134]. The maximum likelihood tree was run for 1000 bootstrap replicates and the tree figure was plotted using FigTree v1.4.3. Protein sequences were analyzed with ScanProsite (https://prosite.expasy.org/scanprosite/, accessed on 5 March 2024) in order to identify active sites and conserved patterns [135].

4. Conclusions

Our annotation of Coccidae immune genes sheds light on the poorly explored repertoire of defense molecules and mechanisms used by a group of insects with a peculiar morphology and habits. The immune gene loss pattern recalls what is observed in aphids, with some exceptions (i.e., genes that are absent in Coccidae, but present in aphids). These are (1) galectins, involved in recognition; (2) some members of the Toll and JNK pathways; and (3) thaumatins, which are antimicrobial effectors also identified in Coleoptera and even in plants. Understanding if this erosion of the immune repertoire is the result of a reduced pressure due to the fungal endosymbiont presence or sterile lifestyles, such as plant-sap feeding, requires further investigations.
Our approach consisted of running searches using the immune gene repertoire of aphids as the query against Coccidae genomics and transcriptomics assemblies. Indeed, aphids are one of the closest relatives of Coccidae, which have a deeply characterized immune gene repertoire. This approach has the advantage of being robust, although somewhat conservative, because it does not take into account alternative pathways and uncharacterized genes that may have evolved in Coccidae (and scale insects, in general), also as an adaptation to the observed loss of immune genes.
However, this work represents the first overview of the immune gene diversity of Coccidae, which we hope will inspire future studies aimed at functionally characterizing the identified genes. Indeed, understanding the role of Coccidae immune genes is a key step for developing new strategies of pest management based on the suppression of the immune response, to enhance the killing activity of entomopathogens.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/ijms25094922/s1.

Author Contributions

Conceptualization, A.B.; formal analysis, E.R. and I.D.L.; investigation, A.B.; data curation, A.B. and I.D.L.; writing—original draft preparation, A.B. and E.R.; writing—review and editing, A.B. and R.N.; visualization, E.R.; supervision, R.N.; funding acquisition, R.N. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All fasta sequences used and alignments are available from Zenodo at https://doi.org/10.5281/zenodo.10864136.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Siva-Jothy, M.T.; Moret, Y.; Rolff, J. Insect immunity: An evolutionary ecology perspective. In Advances in Insect Physiology; Simpson, S.J., Ed.; Academic Press: Cambridge, MA, USA, 2005; Volume 32, pp. 1–48. [Google Scholar]
  2. Ali Mohammadie Kojour, M.; Han, Y.S.; Jo, Y.H. An overview of insect innate immunity. Entomol. Res. 2020, 50, 282–291. [Google Scholar] [CrossRef]
  3. Beckage, N.E. Insect Immunology; Academic Press: Cambridge, MA, USA, 2011. [Google Scholar]
  4. Wang, X.; Zhang, Y.; Zhang, R.; Zhang, J. The diversity of pattern recognition receptors (PRRs) Involved with insect defense against pathogens. Curr. Opin. Insect Sci. 2019, 33, 105–110. [Google Scholar] [CrossRef] [PubMed]
  5. Zhao, L.; Niu, J.; Feng, D.; Wang, X.; Zhang, R. Immune functions of pattern recognition receptors in Lepidoptera. Front. Immunol. 2023, 14, 1203061. [Google Scholar] [CrossRef] [PubMed]
  6. Becchimanzi, A.; Di Lelio, I.; Pennacchio, F.; Caccia, S. Analysis of cellular immune responses in lepidopteran larvae. In Immunity in Insects; Sandrelli, F., Tettamanti, G., Eds.; Springer Protocols Handbooks; Springer: New York, NY, USA, 2020; pp. 97–111. [Google Scholar]
  7. Eleftherianos, I.; Heryanto, C.; Bassal, T.; Zhang, W.; Tettamanti, G.; Mohamed, A. Haemocyte-mediated immunity in insects: Cells, processes and associated components in the fight against pathogens and parasites. Immunology 2021, 164, 401–432. [Google Scholar] [CrossRef] [PubMed]
  8. Strand, M.R. The insect cellular immune response. Insect Sci. 2008, 15, 1–14. [Google Scholar] [CrossRef]
  9. Medzhitov, R.; Janeway, C.A., Jr. Innate immunity: Impact on the adaptive immune response. Curr. Opin. Immunol. 1997, 9, 4–9. [Google Scholar] [CrossRef]
  10. Wu, Q.; Patočka, J.; Kuča, K. Insect antimicrobial peptides, a mini review. Toxins 2018, 10, 461. [Google Scholar] [CrossRef]
  11. Zhang, W.; Tettamanti, G.; Bassal, T.; Heryanto, C.; Eleftherianos, I.; Mohamed, A. Regulators and signalling in insect antimicrobial innate immunity: Functional molecules and cellular pathways. Cell Signal. 2021, 83, 110003. [Google Scholar] [CrossRef] [PubMed]
  12. Iacovone, A.; Ris, N.; Poirié, M.; Gatti, J.-L. Time-course analysis of Drosophila suzukii interaction with endoparasitoid wasps evidences a delayed encapsulation response compared to D. melanogaster. PLoS ONE 2018, 13, e0201573. [Google Scholar] [CrossRef]
  13. Manniello, M.D.; Moretta, A.; Salvia, R.; Scieuzo, C.; Lucchetti, D.; Vogel, H.; Sgambato, A.; Falabella, P. Insect antimicrobial peptides: Potential weapons to counteract the antibiotic resistance. Cell. Mol. Life Sci. 2021, 78, 4259–4282. [Google Scholar] [CrossRef]
  14. Stączek, S.; Cytryńska, M.; Zdybicka-Barabas, A. Unraveling the role of antimicrobial peptides in insects. Int. J. Mol. Sci. 2023, 24, 5753. [Google Scholar] [CrossRef] [PubMed]
  15. Broderick, N.A.; Welchman, D.P.; Lemaitre, B. Recognition and response to microbial infection in Drosophila. In Insect Infection and Immunity: Evolution, Ecology, and Mechanisms; Rolff, J., Reynolds, S., Eds.; Oxford University Press: Oxford, UK, 2009. [Google Scholar]
  16. Jiang, H.; Vilcinskas, A.; Kanost, M.R. Immunity in lepidopteran insects. In Invertebrate Immunity; Springer: Berlin/Heidelberg, Germany, 2010; pp. 181–204. [Google Scholar]
  17. King, J.G. Developmental and comparative perspectives on mosquito immunity. Dev. Comp. Immunol. 2020, 103, 103458. [Google Scholar] [CrossRef] [PubMed]
  18. Liu, Q.-X.; Su, Z.-P.; Liu, H.-H.; Lu, S.-P.; Zhao, Y.; Ma, B.; Hou, Y.-M.; Shi, Z.-H. Current understanding and perspectives on the potential mechanisms of immune priming in beetles. Dev. Comp. Immunol. 2022, 127, 104305. [Google Scholar] [CrossRef] [PubMed]
  19. Morfin, N.; Anguiano-Baez, R.; Guzman-Novoa, E. Honey bee (Apis mellifera) immunity. Vet. Clin. North Am. Food Anim. Pract. 2021, 37, 521–533. [Google Scholar] [CrossRef] [PubMed]
  20. Laughton, A.M.; Garcia, J.R.; Altincicek, B.; Strand, M.R.; Gerardo, N.M. Characterisation of immune responses in the pea aphid, Acyrthosiphon pisum. J. Insect Physiol. 2011, 57, 830–839. [Google Scholar] [CrossRef] [PubMed]
  21. Schmitz, A.; Anselme, C.; Ravallec, M.; Rebuf, C.; Simon, J.-C.; Gatti, J.-L.; Poirié, M. The cellular immune response of the pea aphid to foreign intrusion and symbiotic challenge. PLoS ONE 2012, 7, e42114. [Google Scholar] [CrossRef] [PubMed]
  22. Gerardo, N.M.; Altincicek, B.; Anselme, C.; Atamian, H.; Barribeau, S.M.; de Vos, M.; Duncan, E.J.; Evans, J.D.; Gabaldón, T.; Ghanim, M.; et al. Immunity and other defenses in pea aphids, Acyrthosiphon pisum. Genome Biol. 2010, 11, R21. [Google Scholar] [CrossRef] [PubMed]
  23. Renoz, F.; Noël, C.; Errachid, A.; Foray, V.; Hance, T. Infection dynamic of symbiotic bacteria in the pea aphid Acyrthosiphon pisum gut and host immune response at the early steps in the infection process. PLoS ONE 2015, 10, e0122099. [Google Scholar] [CrossRef] [PubMed]
  24. Brisson, J.A.; Stern, D.L. The pea aphid, Acyrthosiphon pisum: An emerging genomic model system for ecological, developmental and evolutionary studies. Bioessays 2006, 28, 747–755. [Google Scholar] [CrossRef]
  25. Russo, E.; Di Lelio, I.; Shi, M.; Becchimanzi, A.; Pennacchio, F. Aphidius ervi venom regulates Buchnera contribution to host nutritional suitability. J. Insect Physiol. 2023, 147, 104506. [Google Scholar] [CrossRef]
  26. The International Aphid Genomics Consortium. Genome sequence of the pea aphid Acyrthosiphon pisum. PLoS Biol. 2010, 8, e1000313. [Google Scholar]
  27. Arp, A.P.; Hunter, W.B.; Pelz-Stelinski, K.S. Annotation of the Asian citrus psyllid genome reveals a reduced innate immune system. Front. Physiol. 2016, 7, 205526. [Google Scholar] [CrossRef] [PubMed]
  28. Nishide, Y.; Kageyama, D.; Yokoi, K.; Jouraku, A.; Tanaka, H.; Futahashi, R.; Fukatsu, T. Functional crosstalk across IMD and toll pathways: Insight into the evolution of incomplete immune cascades. Proc. R. Soc. B Biol. Sci. 2019, 286, 20182207. [Google Scholar] [CrossRef]
  29. Salcedo-Porras, N.; Guarneri, A.; Oliveira, P.L.; Lowenberger, C. Rhodnius prolixus: Identification of missing components of the IMD immune signaling pathway and functional characterization of its role in eliminating bacteria. PLoS ONE 2019, 14, e0214794. [Google Scholar] [CrossRef]
  30. Hodgson, C. A Review of neococcid scale insects (Hemiptera: Sternorrhyncha: Coccomorpha) based on the morphology of the adult males. Zootaxa 2020, 4765, 1–264. [Google Scholar] [CrossRef]
  31. Gullan, P.J.; Kosztarab, M. Adaptations in scale insects. Annu. Rev. Entomol. 1997, 42, 23–50. [Google Scholar] [CrossRef] [PubMed]
  32. Kono, M.; Koga, R.; Shimada, M.; Fukatsu, T. Infection dynamics of coexisting beta- and gammaproteobacteria in the nested endosymbiotic system of mealybugs. Appl. Environ. Microbiol. 2008, 74, 4175–4184. [Google Scholar] [CrossRef]
  33. Perilla-Henao, L.M.; Casteel, C.L. Vector-borne bacterial plant pathogens: Interactions with hemipteran insects and plants. Front. Plant Sci. 2016, 7, 1163. [Google Scholar] [CrossRef] [PubMed]
  34. Rosenblueth, M.; Martínez-Romero, J.; Ramírez-Puebla, S.T.; Vera-Ponce de León, A.; Rosas-Pérez, T.; Bustamante-Brito, R.; Rincón-Rosales, R.; Martínez-Romero, E.; Rosenblueth, M.; Martínez-Romero, J.; et al. Endosymbiotic microorganisms of scale insects. TIP 2018, 21, 53–69. [Google Scholar]
  35. Szklarzewicz, T.; Michalik, K.; Grzywacz, B.; Kalandyk-Kołodziejczyk, M.; Michalik, A. Fungal associates of soft scale insects (Coccomorpha: Coccidae). Cells 2021, 10, 1922. [Google Scholar] [CrossRef]
  36. Szklarzewicz, T.; Michalik, A.; Michalik, K. The diversity of symbiotic systems in scale insects. In Symbiosis: Cellular, Molecular, Medical and Evolutionary Aspects; Kloc, M., Ed.; Results and Problems in Cell Differentiation; Springer International Publishing: Cham, Switzerland, 2020; pp. 469–495. [Google Scholar]
  37. Baumann, P. Biology of bacteriocyte-associated endosymbionts of plant sap-sucking insects. Annu. Rev. Microbiol. 2005, 59, 155–189. [Google Scholar] [CrossRef] [PubMed]
  38. Vera-Ponce de León, A.; Ormeño-Orrillo, E.; Ramírez-Puebla, S.T.; Rosenblueth, M.; Degli Esposti, M.; Martínez-Romero, J.; Martínez-Romero, E. Candidatus Dactylopiibacterium carminicum, a nitrogen-fixing symbiont of Dactylopius cochineal insects (Hemiptera: Coccoidea: Dactylopiidae). Genome Biol. Evol. 2017, 9, 2237–2250. [Google Scholar] [CrossRef] [PubMed]
  39. Husnik, F.; McCutcheon, J.P. Repeated replacement of an intrabacterial symbiont in the tripartite nested mealybug symbiosis. Proc. Natl. Acad. Sci. USA 2016, 113, E5416–E5424. [Google Scholar] [CrossRef] [PubMed]
  40. Rosas-Pérez, T.; Rosenblueth, M.; Rincón-Rosales, R.; Mora, J.; Martínez-Romero, E. Genome sequence of “Candidatus Walczuchella monophlebidarum” the flavobacterial endosymbiont of Llaveia axin axin (Hemiptera: Coccoidea: Monophlebidae). Genome Biol. Evol. 2014, 6, 714–726. [Google Scholar] [CrossRef]
  41. Matsuura, Y.; Koga, R.; Nikoh, N.; Meng, X.-Y.; Hanada, S.; Fukatsu, T. Huge symbiotic organs in giant scale insects of the genus Drosicha (Coccoidea: Monophlebidae) harbor flavobacterial and enterobacterial endosymbionts. Zool. Sci. 2009, 26, 448–456. [Google Scholar] [CrossRef] [PubMed]
  42. García Morales, M.; Denno, B.D.; Miller, D.R.; Miller, G.L.; Ben-Dov, Y.; Hardy, N.B. ScaleNet: A literature-based model of scale insect biology and systematics. Database 2016, 2016, bav118. [Google Scholar] [PubMed]
  43. Mansour, R.; Grissa-Lebdi, K.; Suma, P.; Mazzeo, G.; Russo, A. Key scale insects (Hemiptera: Coccoidea) of high economic importance in a Mediterranean area: Host plants, bio-ecological characteristics, natural enemies and pest management strategies-a review. Plant Prot. Sci. 2017, 53, 1–14. [Google Scholar] [CrossRef]
  44. Nicoletti, R.; De Masi, L.; Migliozzi, A.; Calandrelli, M.M. Analysis of dieback in a coastal pinewood in Campania, Southern Italy, through high-resolution remote sensing. Plants 2024, 13, 182. [Google Scholar] [CrossRef]
  45. Ramos, A.S.d.J.C.; Lemos, R.N.S.d.; Costa, V.A.; Peronti, A.L.B.G.; Silva, E.A.d.; Mondego, J.M.; Moreira, A.A. Hymenopteran parasitoids associated with scale insects (Hemiptera: Coccoidea) in tropical fruit trees in the Eastern Amazon, Brazil. Fla. Entomol. 2018, 101, 273–278. [Google Scholar] [CrossRef]
  46. Amouroux, P.; Crochard, D.; Correa, M.; Groussier, G.; Kreiter, P.; Roman, C.; Guerrieri, E.; Garonna, A.; Malausa, T.; Zaviezo, T. Natural enemies of armored scales (Hemiptera: Diaspididae) and soft scales (Hemiptera: Coccidae) in Chile: Molecular and morphological identification. PLoS ONE 2019, 14, e0205475. [Google Scholar] [CrossRef]
  47. De León, J.H.; Neumann, G.; Follett, P.A.; Hollingsworth, R.G. Molecular markers discriminate closely related species Encarsia diaspidicola and Encarsia berlesei (Hymenoptera: Aphelinidae): Biocontrol candidate agents for white peach scale in Hawaii. J. Econ. Entomol. 2010, 103, 908–916. [Google Scholar] [CrossRef]
  48. Kundoo, A.A.; Khan, A.A. Coccinellids as biological control agents of soft bodied insects: A review. J. Entomol. Zool. Stud. 2017, 5, 1362–1373. [Google Scholar]
  49. Li, W.; Chen, B.; Toulakhom, C.; Wang, X. Two new species of Chilocorus Leach, 1815 from Laos (Coleoptera Coccinellidae Chilocorini). Biodivers. Data J. 2021, 9, e72966. [Google Scholar] [CrossRef] [PubMed]
  50. Mani, M. Hundred and sixty years of Australian lady bird beetle Crypotolaemus montrouzieri Mulsant—A global view. Biocontrol Sci. Technol. 2018, 28, 938–952. [Google Scholar] [CrossRef]
  51. Liu, W.; Xie, Y.; Dong, J.; Xue, J.; Zhang, Y.; Lu, Y.; Wu, J. Pathogenicity of three entomopathogenic fungi to Matsucoccus matsumurae. PLoS ONE 2014, 9, e103350. [Google Scholar] [CrossRef] [PubMed]
  52. Nicoletti, R.; Becchimanzi, A. Endophytism of Lecanicillium and Akanthomyces. Agriculture 2020, 10, 205. [Google Scholar] [CrossRef]
  53. Dao, H.T.; Beattie, G.A.C.; Rossman, A.Y.; Burgess, L.W.; Holford, P. Four putative entomopathogenic fungi of armoured scale insects on citrus in Australia. Mycol. Prog. 2016, 15, 47. [Google Scholar] [CrossRef]
  54. Khonsanit, A.; Noisripoom, W.; Mongkolsamrit, S.; Phosrithong, N.; Luangsa-ard, J.J. Five new species of Moelleriella infecting scale insects (Coccidae) in Thailand. Mycol. Prog. 2021, 20, 847–867. [Google Scholar] [CrossRef]
  55. Urbina, H.; Ahmed, M.Z. Characterization of the entomopathogenic fungal species Conoideocrella luteorostrata on the scale insect pest Fiorinia externa infesting the Christmas tree Abies fraseri in the USA. Fla. Entomol. 2022, 105, 10–21. [Google Scholar] [CrossRef]
  56. Xu, X.-L.; Zeng, Q.; Lv, Y.-C.; Jeewon, R.; Maharachchikumbura, S.S.; Wanasinghe, D.N.; Hyde, K.D.; Xiao, Q.-G.; Liu, Y.-G.; Yang, C.-L. Insight into the systematics of novel entomopathogenic fungi associated with armored scale insect, Kuwanaspis howardi (Hemiptera: Diaspididae) in China. J. Fungi 2021, 7, 628. [Google Scholar] [CrossRef]
  57. De Lima, I.J.; Carneiro Leão, M.P.; Da Silva Santos, A.C.; Da Costa, A.F.; Tiago, P.V. Production of conidia by entomopathogenic isolates of Fusarium caatingaense on different vegetable substrates. Biocontrol Sci. Technol. 2021, 31, 206–218. [Google Scholar] [CrossRef]
  58. Nicoletti, R.; Russo, E.; Becchimanzi, A. Cladosporium—Insect relationships. J. Fungi 2024, 10, 78. [Google Scholar] [CrossRef]
  59. Qu, J.; Zou, X.; Cao, W.; Xu, Z.; Liang, Z. Two new species of Hirsutella (Ophiocordycipitaceae, Sordariomycetes) that are parasitic on lepidopteran insects from China. MycoKeys 2021, 82, 81. [Google Scholar] [CrossRef] [PubMed]
  60. Deng, J.; Yu, Y.; Wang, X.; Liu, Q.; Huang, X. The ubiquity and development-related abundance dynamics of Ophiocordyceps fungi in soft scale insects. Microorganisms 2021, 9, 404. [Google Scholar] [CrossRef] [PubMed]
  61. Gomez-Polo, P.; Ballinger, M.J.; Lalzar, M.; Malik, A.; Ben-Dov, Y.; Mozes-Daube, N.; Perlman, S.J.; Iasur-Kruh, L.; Chiel, E. An exceptional family: Ophiocordyceps-allied fungus dominates the microbiome of soft scale insects (Hemiptera: Sternorrhyncha: Coccidae). Mol. Ecol. 2017, 26, 5855–5868. [Google Scholar] [CrossRef] [PubMed]
  62. Blumberg, D. Encapsulation of parasitoid eggs in soft scales (Homoptera: Coccidae)*. Ecol. Entomol. 1977, 2, 185–192. [Google Scholar] [CrossRef]
  63. Kapranas, A.; Federici, B.A.; Luck, R.F.; Johnson, J. Cellular immune response of brown soft scale Coccus hesperidum L. (Hemiptera: Coccidae) to eggs of Metaphycus luteolus Timberlake (Hymenoptera: Encyrtidae). Biol. Control 2009, 48, 1–5. [Google Scholar] [CrossRef]
  64. Russo, J.; Allo, M.-R.; Nenon, J.-P.; Brehélin, M. The hemocytes of the mealybugs Phenacoccus manihoti and Planococcus citri (Insecta: Homoptera) and their role in capsule formation. Can. J. Zool. 1994, 72, 252–258. [Google Scholar] [CrossRef]
  65. Caccia, S.; Astarita, F.; Barra, E.; Di Lelio, I.; Varricchio, P.; Pennacchio, F. Enhancement of Bacillus thuringiensis toxicity by feeding Spodoptera littoralis larvae with bacteria expressing immune suppressive dsRNA. J. Pest Sci. 2020, 93, 303–314. [Google Scholar] [CrossRef]
  66. Caccia, S.; Di Lelio, I.; La Storia, A.; Marinelli, A.; Varricchio, P.; Franzetti, E.; Banyuls, N.; Tettamanti, G.; Casartelli, M.; Giordana, B.; et al. Midgut microbiota and host immunocompetence underlie Bacillus thuringiensis killing mechanism. Proc. Natl. Acad. Sci. USA 2016, 113, 9486–9491. [Google Scholar] [CrossRef]
  67. Yang, P.; Yu, S.; Hao, J.; Liu, W.; Zhao, Z.; Zhu, Z.; Sun, T.; Wang, X.; Song, Q. Genome sequence of the Chinese white wax scale insect Ericerus pela: The first draft genome for the Coccidae family of scale insects. GigaScience 2019, 8, giz113. [Google Scholar] [CrossRef] [PubMed]
  68. Kurata, S. Peptidoglycan recognition proteins in Drosophila immunity. Dev. Comp. Immunol. 2014, 42, 36–41. [Google Scholar] [CrossRef] [PubMed]
  69. Cerenius, L.; Söderhäll, K. Immune properties of invertebrate phenoloxidases. Dev. Comp. Immunol. 2021, 122, 104098. [Google Scholar] [CrossRef] [PubMed]
  70. Rao, X.-J.; Zhan, M.-Y.; Pan, Y.-M.; Liu, S.; Yang, P.-J.; Yang, L.-L.; Yu, X.-Q. Immune functions of insect βGRPs and their potential application. Dev. Comp. Immunol. 2018, 83, 80–88. [Google Scholar] [CrossRef] [PubMed]
  71. Takahashi, D.; Garcia, B.L.; Kanost, M.R. Initiating Protease with modular domains interacts with β-glucan recognition protein to trigger innate immune response in insects. Proc. Natl. Acad. Sci. USA 2015, 112, 13856–13861. [Google Scholar] [CrossRef] [PubMed]
  72. Warr, E.; Das, S.; Dong, Y.; Dimopoulos, G. The Gram-negative bacteria-binding protein gene family: Its role in the innate immune system of Anopheles gambiae and in anti-Plasmodium defence. Insect Mol. Biol. 2008, 17, 39–51. [Google Scholar] [CrossRef] [PubMed]
  73. Pili-Floury, S.; Leulier, F.; Takahashi, K.; Saigo, K.; Samain, E.; Ueda, R.; Lemaitre, B. In vivo RNA interference analysis reveals an unexpected role for GNBP1 in the defense against Gram-positive bacterial infection in Drosophila adults *. J. Biol. Chem. 2004, 279, 12848–12853. [Google Scholar] [CrossRef] [PubMed]
  74. Gottar, M.; Gobert, V.; Matskevich, A.A.; Reichhart, J.-M.; Wang, C.; Butt, T.M.; Belvin, M.; Hoffmann, J.A.; Ferrandon, D. Dual detection of fungal infections in Drosophila via recognition of glucans and sensing of virulence factors. Cell 2006, 127, 1425–1437. [Google Scholar] [CrossRef]
  75. Ji, J.; Zhou, L.; Xu, Z.; Ma, L.; Lu, Z. Two atypical Gram-negative bacteria-binding proteins are involved in the antibacterial response in the pea aphid (Acyrthosiphon pisum). Insect Mol. Biol. 2021, 30, 427–435. [Google Scholar] [CrossRef]
  76. Palmer, W.J.; Jiggins, F.M. Comparative genomics reveals the origins and diversity of arthropod immune systems. Mol. Biol. Evol. 2015, 32, 2111–2129. [Google Scholar] [CrossRef]
  77. Ao, J.; Ling, E.; Yu, X.-Q. Drosophila C-type lectins enhance cellular encapsulation. Mol. Immunol. 2007, 44, 2541–2548. [Google Scholar] [CrossRef] [PubMed]
  78. Xia, X.; You, M.; Rao, X.-J.; Yu, X.-Q. Insect C-type lectins in innate immunity. Dev. Comp. Immunol. 2018, 83, 70–79. [Google Scholar] [CrossRef] [PubMed]
  79. Pace, K.E.; Baum, L.G. Insect galectins: Roles in immunity and development. Glycoconj. J. 2002, 19, 607–614. [Google Scholar] [CrossRef] [PubMed]
  80. Dimopoulos, G.; Richman, A.; della Torre, A.; Kafatos, F.C.; Louis, C. Identification and characterization of differentially expressed cDNAs of the vector mosquito, Anopheles gambiae. Proc. Natl. Acad. Sci. USA 1996, 93, 13066–13071. [Google Scholar] [CrossRef]
  81. Rao, X.-J.; Wu, P.; Shahzad, T.; Liu, S.; Chen, L.; Yang, Y.-F.; Shi, Q.; Yu, X.-Q. Characterization of a dual-CRD galectin in the silkworm Bombyx Mori. Dev. Comp. Immunol. 2016, 60, 149–159. [Google Scholar] [CrossRef] [PubMed]
  82. Waterhouse, R.M.; Kriventseva, E.V.; Meister, S.; Xi, Z.; Alvarez, K.S.; Bartholomay, L.C.; Barillas-Mury, C.; Bian, G.; Blandin, S.; Christensen, B.M.; et al. Evolutionary dynamics of immune-related genes and pathways in disease-vector mosquitoes. Science 2007, 316, 1738–1743. [Google Scholar] [CrossRef] [PubMed]
  83. Lin, Z.; Wang, J.-L.; Cheng, Y.; Wang, J.-X.; Zou, Z. Pattern recognition receptors from lepidopteran insects and their biological functions. Dev. Comp. Immunol. 2020, 108, 103688. [Google Scholar] [CrossRef] [PubMed]
  84. Blandin, S.A.; Marois, E.; Levashina, E.A. Antimalarial responses in Anopheles gambiae: From a complement-like protein to a complement-like pathway. Cell Host Microbe 2008, 3, 364–374. [Google Scholar] [CrossRef] [PubMed]
  85. Theopold, U.; Schmid, M. Thioester-containing proteins: At the crossroads of immune effector mechanisms. Virulence 2017, 8, 1468–1470. [Google Scholar] [CrossRef] [PubMed]
  86. Lagueux, M.; Perrodou, E.; Levashina, E.A.; Capovilla, M.; Hoffmann, J.A. Constitutive expression of a complement-like protein in toll and JAK gain-of-function mutants of Drosophila. Proc. Natl. Acad. Sci. USA 2000, 97, 11427–11432. [Google Scholar] [CrossRef]
  87. Levashina, E.A.; Moita, L.F.; Blandin, S.; Vriend, G.; Lagueux, M.; Kafatos, F.C. Conserved role of a complement-like protein in phagocytosis revealed by dsRNA knockout in cultured cells of the mosquito, Anopheles gambiae. Cell 2001, 104, 709–718. [Google Scholar] [CrossRef] [PubMed]
  88. Rämet, M.; Pearson, A.; Manfruelli, P.; Li, X.; Koziel, H.; Göbel, V.; Chung, E.; Krieger, M.; Ezekowitz, R.A. Drosophila scavenger receptor CI is a pattern recognition receptor for bacteria. Immunity 2001, 15, 1027–1038. [Google Scholar] [CrossRef] [PubMed]
  89. Somogyi, K.; Sipos, B.; Pénzes, Z.; Kurucz, E.; Zsámboki, J.; Hultmark, D.; Andó, I. Evolution of genes and repeats in the nimrod superfamily. Mol. Biol. Evol. 2008, 25, 2337–2347. [Google Scholar] [CrossRef] [PubMed]
  90. Ju, J.S.; Cho, M.H.; Brade, L.; Kim, J.H.; Park, J.W.; Ha, N.-C.; Söderhäll, I.; Söderhäll, K.; Brade, H.; Lee, B.L. A novel 40-kDa protein containing six repeats of an epidermal growth factor-like domain functions as a pattern recognition protein for lipopolysaccharide. J. Immunol. 2006, 177, 1838–1845. [Google Scholar] [CrossRef] [PubMed]
  91. Kurucz, E.; Márkus, R.; Zsámboki, J.; Folkl-Medzihradszky, K.; Darula, Z.; Vilmos, P.; Udvardy, A.; Krausz, I.; Lukacsovich, T.; Gateff, E.; et al. Nimrod, a putative phagocytosis receptor with EGF repeats in Drosophila plasmatocytes. Curr. Biol. 2007, 17, 649–654. [Google Scholar] [CrossRef] [PubMed]
  92. Melcarne, C.; Ramond, E.; Dudzic, J.; Bretscher, A.J.; Kurucz, É.; Andó, I.; Lemaitre, B. Two nimrod receptors, NimC1 and Eater, synergistically contribute to bacterial phagocytosis in Drosophila melanogaster. FEBS J. 2019, 286, 2670–2691. [Google Scholar] [CrossRef] [PubMed]
  93. Boutros, M.; Agaisse, H.; Perrimon, N. Sequential activation of signaling pathways during innate immune responses in Drosophila. Dev. Cell 2002, 3, 711–722. [Google Scholar] [CrossRef] [PubMed]
  94. Yang, L.; Qiu, L.-M.; Fang, Q.; Stanley, D.W.; Ye, G.-Y. Cellular and humoral immune interactions between Drosophila and its parasitoids. Insect Sci. 2021, 28, 1208–1227. [Google Scholar] [CrossRef] [PubMed]
  95. Lima, L.F.; Torres, A.Q.; Jardim, R.; Mesquita, R.D.; Schama, R. Evolution of Toll, Spatzle and MyD88 in insects: The problem of the Diptera bias. BMC Genom. 2021, 22, 562. [Google Scholar] [CrossRef]
  96. Leulier, F.; Lemaitre, B. Toll-like receptors—Taking an evolutionary approach. Nat. Rev. Genet. 2008, 9, 165–178. [Google Scholar] [CrossRef]
  97. Agaisse, H.; Perrimon, N. The roles of JAK/STAT signaling in Drosophila immune responses. Immunol. Rev. 2004, 198, 72–82. [Google Scholar] [CrossRef] [PubMed]
  98. Schlenke, T.A.; Morales, J.; Govind, S.; Clark, A.G. Contrasting infection strategies in generalist and specialist wasp parasitoids of Drosophila melanogaster. PLoS Pathog. 2007, 3, e158. [Google Scholar] [CrossRef] [PubMed]
  99. Evans, J.D.; Aronstein, K.; Chen, Y.P.; Hetru, C.; Imler, J.-L.; Jiang, H.; Kanost, M.; Thompson, G.J.; Zou, Z.; Hultmark, D. Immune pathways and defence mechanisms in honey bees Apis mellifera. Insect Mol. Biol. 2006, 15, 645–656. [Google Scholar] [CrossRef] [PubMed]
  100. Dionne, M.S.; Schneider, D.S. Models of infectious diseases in the fruit fly Drosophila melanogaster. Dis. Models Mech. 2008, 1, 43–49. [Google Scholar] [CrossRef] [PubMed]
  101. Chowdhury, A.; Modahl, C.M.; Tan, S.T.; Xiang, B.W.W.; Missé, D.; Vial, T.; Kini, R.M.; Pompon, J.F. JNK pathway restricts DENV2, ZIKV and CHIKV infection by activating complement and apoptosis in mosquito salivary glands. PLoS Pathog. 2020, 16, e1008754. [Google Scholar] [CrossRef]
  102. Hanson, M.A.; Lemaitre, B. New insights on Drosophila antimicrobial peptide function in host defense and beyond. Curr. Opin. Immunol. 2020, 62, 22–30. [Google Scholar] [CrossRef] [PubMed]
  103. Jing, X.; Wong, A.C.-N.; Chaston, J.M.; Colvin, J.; McKenzie, C.L.; Douglas, A.E. The bacterial communities in plant phloem-sap-feeding insects. Mol. Ecol. 2014, 23, 1433–1444. [Google Scholar] [CrossRef] [PubMed]
  104. Bachali, S.; Jager, M.; Hassanin, A.; Schoentgen, F.; Jollès, P.; Fiala-Medioni, A.; Deutsch, J.S. Phylogenetic analysis of invertebrate lysozymes and the evolution of lysozyme function. J. Mol. Evol. 2002, 54, 652–664. [Google Scholar] [CrossRef] [PubMed]
  105. Van Herreweghe, J.M.; Michiels, C.W. Invertebrate lysozymes: Diversity and distribution, molecular mechanism and in vivo function. J. Biosci. 2012, 37, 327–348. [Google Scholar] [CrossRef]
  106. Marra, A.; Hanson, M.A.; Kondo, S.; Erkosar, B.; Lemaitre, B. Drosophila antimicrobial peptides and lysozymes regulate gut microbiota composition and abundance. mBio 2021, 12, 10–1128. [Google Scholar] [CrossRef]
  107. Paskewitz, S.M.; Li, B.; Kajla, M.K. Cloning and molecular characterization of two invertebrate-type lysozymes from Anopheles gambiae. Insect Mol. Biol. 2008, 17, 217–225. [Google Scholar] [CrossRef] [PubMed]
  108. Beckert, A.; Wiesner, J.; Schmidtberg, H.; Lehmann, R.; Baumann, A.; Vogel, H.; Vilcinskas, A. Expression and characterization of a recombinant I-type lysozyme from the harlequin ladybird beetle Harmonia axyridis. Insect Mol. Biol. 2016, 25, 202–215. [Google Scholar] [CrossRef] [PubMed]
  109. Moran, N.A.; McCutcheon, J.P.; Nakabachi, A. Genomics and evolution of heritable bacterial symbionts. Annu. Rev. Genet. 2008, 42, 165–190. [Google Scholar] [CrossRef] [PubMed]
  110. Nikoh, N.; Nakabachi, A. Aphids acquired symbiotic genes via lateral gene transfer. BMC Biol. 2009, 7, 12. [Google Scholar] [CrossRef]
  111. Templin, M.F.; Ursinus, A.; Höltje, J. A defect in cell wall recycling triggers autolysis during the stationary growth phase of Escherichia coli. EMBO J. 1999, 18, 4108–4117. [Google Scholar] [CrossRef] [PubMed]
  112. Jorgenson, M.A.; Chen, Y.; Yahashiri, A.; Popham, D.L.; Weiss, D.S. The bacterial septal ring protein RlpA is a lytic transglycosylase that contributes to rod shape and daughter cell separation in Pseudomonas aeruginosa. Mol. Microbiol. 2014, 93, 113–128. [Google Scholar] [CrossRef] [PubMed]
  113. Nakabachi, A. Horizontal gene transfers in insects. Curr. Opin. Insect Sci. 2015, 7, 24–29. [Google Scholar] [CrossRef]
  114. Sun, T.; Wang, X.-Q.; Zhao, Z.-L.; Yu, S.-H.; Yang, P.; Chen, X.-M. A lethal fungus infects the Chinese white wax scale insect and causes dramatic changes in the host microbiota. Sci. Rep. 2018, 8, 5324. [Google Scholar] [CrossRef] [PubMed]
  115. Shigenobu, S.; Stern, D.L. Aphids evolved novel secreted proteins for symbiosis with bacterial endosymbiont. Proc. R. Soc. B Biol. Sci. 2013, 280, 20121952. [Google Scholar] [CrossRef]
  116. Loth, K.; Parisot, N.; Paquet, F.; Terrasson, H.; Sivignon, C.; Rahioui, I.; Ribeiro Lopes, M.; Gaget, K.; Duport, G.; Delmas, A.F.; et al. Aphid BCR4 structure and activity uncover a new defensin peptide superfamily. Int. J. Mol. Sci. 2022, 23, 12480. [Google Scholar] [CrossRef]
  117. Shahidi, F.; Abuzaytoun, R. Chitin, chitosan, and co-products: Chemistry, production, applications, and health effects. In Advances in Food and Nutrition Research; Academic Press: Cambridge, MA, USA, 2005; Volume 49, pp. 93–135. [Google Scholar]
  118. Arakane, Y.; Muthukrishnan, S. Insect chitinase and chitinase-like proteins. Cell. Mol. Life Sci. 2010, 67, 201–216. [Google Scholar] [CrossRef] [PubMed]
  119. Ding, J.-L.; Hou, J.; Feng, M.-G.; Ying, S.-H. Transcriptomic analyses reveal comprehensive responses of insect hemocytes to mycopathogen Beauveria bassiana, and fungal virulence-related cell wall protein assists pathogen to evade host cellular defense. Virulence 2020, 11, 1352–1365. [Google Scholar] [CrossRef] [PubMed]
  120. Jollès, P. Lysozymes—Model Enzymes in Biochemistry and Biology; Birkhäuser: Basel, Switzerland, 1996. [Google Scholar]
  121. Sustar, A.E.; Strand, L.G.; Zimmerman, S.G.; Berg, C.A. Imaginal disk growth factors are Drosophila chitinase-like proteins with roles in morphogenesis and CO2 response. Genetics 2023, 223, iyac185. [Google Scholar] [CrossRef]
  122. Nappi, A.J.; Christensen, B.M. Melanogenesis and associated cytotoxic reactions: Applications to insect innate immunity. Insect Biochem. Mol. Biol. 2005, 35, 443–459. [Google Scholar] [CrossRef]
  123. Xu, L.; Ma, L.; Wang, W.; Li, L.; Lu, Z. Phenoloxidases are required for the pea aphid’s defence against bacterial and fungal infection. Insect Mol. Biol. 2019, 28, 176–186. [Google Scholar] [CrossRef] [PubMed]
  124. Ma, L.; Chen, F.; Wang, W.; Xu, L.; Lu, Z.-Q. Identification of two clip domain serine proteases involved in the pea aphid’s defense against bacterial and fungal infection. Insect Sci. 2020, 27, 735–744. [Google Scholar] [CrossRef]
  125. McLean, A.H.C.; Parker, B.J. Variation in intrinsic resistance of pea aphids to parasitoid wasps: A transcriptomic basis. PLoS ONE 2020, 15, e0242159. [Google Scholar] [CrossRef]
  126. Rivero, A. Nitric oxide: An antiparasitic molecule of invertebrates. Trends Parasitol. 2006, 22, 219–225. [Google Scholar] [CrossRef]
  127. Sadekuzzaman, M.; Kim, Y. Nitric oxide mediates antimicrobial peptide gene expression by activating eicosanoid signaling. PLoS ONE 2018, 13, e0193282. [Google Scholar] [CrossRef]
  128. Nappi, A.J.; Vass, E.; Frey, F.; Carton, Y. Nitric oxide involvement in Drosophila immunity. Nitric Oxide 2000, 4, 423–430. [Google Scholar] [CrossRef]
  129. Ma, L.; Yan, X.; Zhou, L.; Wang, W.; Chen, K.; Hao, C.; Lu, Z.; Qie, X. Nitric oxide synthase is required for the pea aphid’s defence against bacterial infection. Insect Mol. Biol. 2023, 32, 187–199. [Google Scholar] [CrossRef] [PubMed]
  130. Slater, G.S.C.; Birney, E. Automated generation of heuristics for biological sequence comparison. BMC Bioinform. 2005, 6, 31. [Google Scholar] [CrossRef] [PubMed]
  131. Rice, P.; Longden, I.; Bleasby, A. EMBOSS: The European molecular biology open software suite. Trends Genet. 2000, 16, 276–277. [Google Scholar] [CrossRef] [PubMed]
  132. Edgar, R.C. MUSCLE: Multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32, 1792–1797. [Google Scholar] [CrossRef] [PubMed]
  133. Castresana, J. Selection of conserved blocks from multiple alignments for their use in phylogenetic analysis. Mol. Biol. Evol. 2000, 17, 540–552. [Google Scholar] [CrossRef] [PubMed]
  134. Stamatakis, A. RAxML-VI-HPC: Maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics 2006, 22, 2688–2690. [Google Scholar] [CrossRef]
  135. de Castro, E.; Sigrist, C.J.A.; Gattiker, A.; Bulliard, V.; Langendijk-Genevaux, P.S.; Gasteiger, E.; Bairoch, A.; Hulo, N. ScanProsite: Detection of PROSITE signature matches and prorule-associated functional and structural residues in proteins. Nucleic Acids Res. 2006, 34, W362–W365. [Google Scholar] [CrossRef]
Figure 1. Overview of gene families involved in recognition, signaling, and response pathways in Ericerus pela, Acyrthosiphon pisum, and Drosophila melanogaster. Numbers in blocks indicate the different genes identified for each protein family. Abbreviations: peptidoglycan recognition (PGRP), Gram-negative binding (GNBP), thioester-containing (TEP), scavenger receptor class C (SRC), antimicrobial peptide (AMP), pro-phenoloxidase (PPO), phenoloxidase (PO), nitric oxide (NO).
Figure 1. Overview of gene families involved in recognition, signaling, and response pathways in Ericerus pela, Acyrthosiphon pisum, and Drosophila melanogaster. Numbers in blocks indicate the different genes identified for each protein family. Abbreviations: peptidoglycan recognition (PGRP), Gram-negative binding (GNBP), thioester-containing (TEP), scavenger receptor class C (SRC), antimicrobial peptide (AMP), pro-phenoloxidase (PPO), phenoloxidase (PO), nitric oxide (NO).
Ijms 25 04922 g001
Figure 2. Phylogenetic tree based on maximum likelihood analysis of GNBP1 putative homologs identified in E. pela genome (in red) and in the transcriptomes of C. cirripediformis (in blue) and Coccus sp. (in purple). GNBP1 sequences of Coccidae family form a monophyletic group, closely related to GNBP1 of A. pisum. Sequences of E. pela are two isoforms encoded by the same gene, while Coccus sp. shows two different genes encoding homologs of GNBP1. One of these proteins, GCWW01035969.1, is more closely related to the GNBP1 identified in C. cirripediformis. The longest branch of the unrooted tree is used as the outgroup. Bootstrap support values are indicated at each node. The scale bar indicates the number of amino acid substitutions per site.
Figure 2. Phylogenetic tree based on maximum likelihood analysis of GNBP1 putative homologs identified in E. pela genome (in red) and in the transcriptomes of C. cirripediformis (in blue) and Coccus sp. (in purple). GNBP1 sequences of Coccidae family form a monophyletic group, closely related to GNBP1 of A. pisum. Sequences of E. pela are two isoforms encoded by the same gene, while Coccus sp. shows two different genes encoding homologs of GNBP1. One of these proteins, GCWW01035969.1, is more closely related to the GNBP1 identified in C. cirripediformis. The longest branch of the unrooted tree is used as the outgroup. Bootstrap support values are indicated at each node. The scale bar indicates the number of amino acid substitutions per site.
Ijms 25 04922 g002
Figure 3. Phylogenetic tree based on maximum likelihood analysis of CTL putative homologs identified in E. pela genome (in red) and in the transcriptomes of C. cirripediformis (in blue) and Coccus sp. (in purple). CTL sequences of Coccidae family form two monophyletic groups, one with c-type domain signature (C_TYPE_LECTIN_1) and profile (C_TYPE_LECTIN_2), and another with c-type domain profile alone (C_TYPE_LECTIN_2), marked in yellow and green, respectively. Sequences of E. pela in the yellow group are two isoforms encoded by the same gene. The longest branch of the unrooted tree (DL3) is used as the outgroup. Bootstrap support values are indicated at each node. The scale bar indicates the number of amino acid substitutions per site.
Figure 3. Phylogenetic tree based on maximum likelihood analysis of CTL putative homologs identified in E. pela genome (in red) and in the transcriptomes of C. cirripediformis (in blue) and Coccus sp. (in purple). CTL sequences of Coccidae family form two monophyletic groups, one with c-type domain signature (C_TYPE_LECTIN_1) and profile (C_TYPE_LECTIN_2), and another with c-type domain profile alone (C_TYPE_LECTIN_2), marked in yellow and green, respectively. Sequences of E. pela in the yellow group are two isoforms encoded by the same gene. The longest branch of the unrooted tree (DL3) is used as the outgroup. Bootstrap support values are indicated at each node. The scale bar indicates the number of amino acid substitutions per site.
Ijms 25 04922 g003
Figure 4. Phylogenetic tree based on maximum likelihood analysis of thioester-containing protein (TEP) putative homologs identified in genome of E. pela (in red) and in the transcriptomes of C. cirripediformis (in blue) and Coccus sp. (in purple). Sequences of C. cirripediformis are two isoforms encoded by the same gene. The longest branch of the unrooted tree is used as the outgroup. Bootstrap support values are indicated at each node. The scale bar indicates the number of amino acid substitutions per site.
Figure 4. Phylogenetic tree based on maximum likelihood analysis of thioester-containing protein (TEP) putative homologs identified in genome of E. pela (in red) and in the transcriptomes of C. cirripediformis (in blue) and Coccus sp. (in purple). Sequences of C. cirripediformis are two isoforms encoded by the same gene. The longest branch of the unrooted tree is used as the outgroup. Bootstrap support values are indicated at each node. The scale bar indicates the number of amino acid substitutions per site.
Ijms 25 04922 g004
Figure 5. Phylogenetic tree based on maximum likelihood analysis of Toll-1, Toll-6, Toll-7, Toll-10, and Tollo putative homologs identified in E. pela genome (in red) and in the transcriptomes of C. cirripediformis (in blue) and Coccus sp. (in purple). E. pela sequences in the Toll-1 group are two isoforms encoded by the same gene. The A. pisum sequence XP_003248960.1, which is annotated as the Toll-6 receptor, is here annotated as Toll-10 following the classification of [95]. The longest branch of the unrooted tree is used as the outgroup. Bootstrap support values are indicated at each node. The scale bar indicates the number of amino acid substitutions per site.
Figure 5. Phylogenetic tree based on maximum likelihood analysis of Toll-1, Toll-6, Toll-7, Toll-10, and Tollo putative homologs identified in E. pela genome (in red) and in the transcriptomes of C. cirripediformis (in blue) and Coccus sp. (in purple). E. pela sequences in the Toll-1 group are two isoforms encoded by the same gene. The A. pisum sequence XP_003248960.1, which is annotated as the Toll-6 receptor, is here annotated as Toll-10 following the classification of [95]. The longest branch of the unrooted tree is used as the outgroup. Bootstrap support values are indicated at each node. The scale bar indicates the number of amino acid substitutions per site.
Ijms 25 04922 g005
Figure 6. Phylogenetic tree based on maximum likelihood analysis of interleukin JAK/STAT receptor (domeless) putative homologs identified in E. pela genome (in red) and in the transcriptomes of C. cirripediformis (in blue) and Coccus sp. (in purple). Sequences of E. pela are two isoforms encoded by the same gene. The closely related dscam2 (Down syndrome cell adhesion) sequences are used as the outgroup. Bootstrap support values are indicated at each node. The scale bar indicates the number of amino acid substitutions per site.
Figure 6. Phylogenetic tree based on maximum likelihood analysis of interleukin JAK/STAT receptor (domeless) putative homologs identified in E. pela genome (in red) and in the transcriptomes of C. cirripediformis (in blue) and Coccus sp. (in purple). Sequences of E. pela are two isoforms encoded by the same gene. The closely related dscam2 (Down syndrome cell adhesion) sequences are used as the outgroup. Bootstrap support values are indicated at each node. The scale bar indicates the number of amino acid substitutions per site.
Ijms 25 04922 g006
Figure 7. Phylogenetic tree based on maximum likelihood analysis of lysozyme putative homologs identified in E. pela genome (in red) and in the transcriptomes of C. cirripediformis (in blue) and Coccus sp. (in purple). Lysozyme sequences of Coccidae family form three monophyletic groups, one with the invertebrate (I)-type lysozyme domain profile (LYSOZYME_I) alone, and another with LYSOZYME_I and EF-hand calcium-binding domain (EF_HAND_1), marked in light and dark blue, respectively. The longest branch of the unrooted tree is used as the outgroup. Bootstrap support values are indicated at each node. The scale bar indicates the number of amino acid substitutions per site.
Figure 7. Phylogenetic tree based on maximum likelihood analysis of lysozyme putative homologs identified in E. pela genome (in red) and in the transcriptomes of C. cirripediformis (in blue) and Coccus sp. (in purple). Lysozyme sequences of Coccidae family form three monophyletic groups, one with the invertebrate (I)-type lysozyme domain profile (LYSOZYME_I) alone, and another with LYSOZYME_I and EF-hand calcium-binding domain (EF_HAND_1), marked in light and dark blue, respectively. The longest branch of the unrooted tree is used as the outgroup. Bootstrap support values are indicated at each node. The scale bar indicates the number of amino acid substitutions per site.
Ijms 25 04922 g007
Figure 8. Phylogenetic tree based on maximum likelihood analysis of ldca (LD carboxypeptidase) putative homologs identified in E. pela genome (in red). The sequence of E. pela groups with LD carboxypeptidases of Rickettsia sp., which is a bacterial symbiont of several insect species, including E. pela. Sequences identified in Wolbachia species are used as the outgroup. Bootstrap support values are indicated at each node. The scale bar indicates the number of amino acid substitutions per site.
Figure 8. Phylogenetic tree based on maximum likelihood analysis of ldca (LD carboxypeptidase) putative homologs identified in E. pela genome (in red). The sequence of E. pela groups with LD carboxypeptidases of Rickettsia sp., which is a bacterial symbiont of several insect species, including E. pela. Sequences identified in Wolbachia species are used as the outgroup. Bootstrap support values are indicated at each node. The scale bar indicates the number of amino acid substitutions per site.
Ijms 25 04922 g008
Figure 9. Phylogenetic tree based on maximum likelihood analysis of chitinase putative homologs identified in the E. pela genome (in red) and in the transcriptomes of C. cirripediformis (in blue) and Coccus sp. (in purple). Chitinase sequences of the Coccidae family are included in all the groups identified in [118], except for groups IV and VII. Chitinase groups are marked by a vertical colored bar, along with their domain architecture and length, except for sequence GCWZ01013619.1_1 of C. cirripediformis, which is the only member identified in Coccidae belonging to group II. Domain architecture identified by ScanProsite consists of the glycosyl hydrolases family 18 (GH18) domain (GH18_2) and active site (GH18_1) and chitin-binding type-2 domain (CHIT_BIND_II). Imaginal disk growth factor (group V) chitinase-like proteins, involved in morphogenesis and CO2 response rather than immunity [121], are used as the outgroup. Bootstrap support values are indicated at each node. The scale bar indicates the number of amino acid substitutions per site.
Figure 9. Phylogenetic tree based on maximum likelihood analysis of chitinase putative homologs identified in the E. pela genome (in red) and in the transcriptomes of C. cirripediformis (in blue) and Coccus sp. (in purple). Chitinase sequences of the Coccidae family are included in all the groups identified in [118], except for groups IV and VII. Chitinase groups are marked by a vertical colored bar, along with their domain architecture and length, except for sequence GCWZ01013619.1_1 of C. cirripediformis, which is the only member identified in Coccidae belonging to group II. Domain architecture identified by ScanProsite consists of the glycosyl hydrolases family 18 (GH18) domain (GH18_2) and active site (GH18_1) and chitin-binding type-2 domain (CHIT_BIND_II). Imaginal disk growth factor (group V) chitinase-like proteins, involved in morphogenesis and CO2 response rather than immunity [121], are used as the outgroup. Bootstrap support values are indicated at each node. The scale bar indicates the number of amino acid substitutions per site.
Ijms 25 04922 g009
Figure 10. Phylogenetic tree based on maximum likelihood analysis of phenoloxidase activating factor 2 (paf2) putative homologs identified in the E. pela genome (in red) and in the transcriptomes of C. cirripediformis (in blue) and Coccus sp. (in purple). Each Coccidae species has two paralogs belonging to the paf2 family. The sequences of the three Coccidae species form two monophyletic groups, whose most closely related sequences are their orthologs in the A. pisum genome. The longest branch of the unrooted tree is used as the outgroup. Bootstrap support values are indicated at each node. The scale bar indicates the number of amino acid substitutions per site.
Figure 10. Phylogenetic tree based on maximum likelihood analysis of phenoloxidase activating factor 2 (paf2) putative homologs identified in the E. pela genome (in red) and in the transcriptomes of C. cirripediformis (in blue) and Coccus sp. (in purple). Each Coccidae species has two paralogs belonging to the paf2 family. The sequences of the three Coccidae species form two monophyletic groups, whose most closely related sequences are their orthologs in the A. pisum genome. The longest branch of the unrooted tree is used as the outgroup. Bootstrap support values are indicated at each node. The scale bar indicates the number of amino acid substitutions per site.
Ijms 25 04922 g010
Table 1. Immune genes of Ericerus pela involved in recognition. Genes not found in E. pela are colored in red.
Table 1. Immune genes of Ericerus pela involved in recognition. Genes not found in E. pela are colored in red.
Role in Insect ImmunityGene SymbolGene Name 1D. melanogaster
NCBI Protein ID
A. pisum
NCBI Protein ID
Best Matches (E-Value) 2
bacterial recognitionPGRP-LCpeptidoglycan recognition proteinAAF50302.3not foundnot found
activation of PPO cascade and autophagyPGRP-LEpeptidoglycan recognition proteinNP_573078.1not foundnot found
bacterial recognitionPGRP-SApeptidoglycan recognition proteinAAF48056.1not foundnot found
bacterial recognitionPGRP-SDpeptidoglycan recognition proteinCAD89193.1not foundnot found
bacterial recognitionPGRP-LBpeptidoglycan recognition proteinNP_650079.1not foundnot found
bacterial recognitionPGRP-SC1apeptidoglycan recognition proteinCAD89161.1not foundnot found
bacterial recognitionPGRP-SC2peptidoglycan recognition proteinCAD89187.1not foundnot found
pgn degradation and antibacterial activityPGRP-SB1peptidoglycan recognition proteinCAD89136.1not foundnot found
blocking of imd pathwayPGRP-LFpeptidoglycan recognition proteinNP_648299.3not foundnot found
activation of imd pathwayPGRP-LApeptidoglycan recognition proteinAAF50304.2not foundnot found
bacterial and fungal pattern recognitionGNBP1Gram-negative binding protein 1Q9NHB0.2XP_001944473.2QBOQ01000878.1 (3 × 10−8)
QBOQ01000589.1 (2 × 10−6)
bacterial and fungal pattern recognitionGNBP2Gram-negative binding protein 2ACU30172.1XP_001944473.2QBOQ01000878.1 (3 × 10−8)
QBOQ01000589.1 (2 × 10−6)
bacterial and fungal pattern recognitionGNBP3Gram-negative binding protein 3CAJ18910.1XP_029342159.1QBOQ01000878.1 (1 × 10−18)
bacterial recognition, induction of PPO cascadeDL1c-type lectin 1AAF53793.1not foundnot found
bacterial recognition, induction of PPO cascadeDL2c-type lectin 2NP_001014489.1XP_016663197.1
XP_001950803.2
XP_001945032.2
QBOQ01000461.1 (1 × 10−52) QBOQ01000461.1 (3 × 10−53)
QBOQ01000466.1 (4 × 10−24)
bacterial recognition, induction of PPO cascadeDL3c-type lectin 3NP_001014490.1XP_016663197.1
XP_001950803.2
QBOQ01000461.1 (1 × 10−52) QBOQ01000461.1 (3 × 10−53)
several roles have been hypothesizedgalectingalectin 9ADZ99399.1XP_001943769.2not found
mark pathogens for phagocytosisTep1thioester containing protein 1CAB87807.1XP_029348718.1QBOQ01000202.1 (2 × 10−38)
mark pathogens for phagocytosisTep2thioester containing protein 2CAB87808.1XP_029348718.1QBOQ01000202.1 (2 × 10−38)
mark pathogens for phagocytosisTep3thioester containing protein 3AAL39195.1XP_029348714.1QBOQ01000202.1 (2 × 10−38)
mark pathogens for phagocytosisTep4thioester containing protein 4NP_523603.2XP_029348718.1QBOQ01000202.1 (2 × 10−38)
bacterial and fungal recognitionpespeste, scavenger receptor class bAHN54246.1XP_029341846.1QBOQ01001218.1 (2 × 10−33)
bacterial and fungal recognitioncrqcroquemortAAF51494.1XP_001944867.2QBOQ01000024.1 (1 × 10−35)
bacterial and fungal recognitiondrprdraperNP_477450.1XP_001942552.2QBOQ01001915.1 (2 × 10−29)
bind to lipoproteins and bacteriasr-CIscavenger receptor class c, type iAAW79470.1not foundnot found
bind to lipoproteins and bacteriasr-CIIscavenger receptor class c, type iiAAF58551.1not foundnot found
bind to lipoproteins and bacteriasr-CIIIscavenger receptor class c, type iiiAAF37564.1not foundnot found
bind to lipoproteins and bacteriasr-CIVscavenger receptor class c, type ivAAF51092.1not foundnot found
receptor in phagocytosis and microbial bindingeatereaterAAF56664.5not foundnot found
receptor in phagocytosis and microbial bindingnim-C1Nimrod c1AAF53364.2not foundnot found
1 Alternative names are separated by commas. 2 E-values refer to A. pisum proteins, if present, or D. melanogaster if the gene was not found in A. pisum.
Table 2. Immune genes of Ericerus pela involved in signaling. Genes not found in E. pela are colored in red.
Table 2. Immune genes of Ericerus pela involved in signaling. Genes not found in E. pela are colored in red.
Role in Insect
Immunity
Gene
Symbol
Gene Name 1D. melanogaster
NCBI Protein ID
A. pisum
NCBI Protein ID
Best Matches (E-Value) 2
Toll pathwayspz1-1spätzle 1BNP_733188.1NP_001153589QBOQ01001621.1 (7 × 10−13)
Toll pathwayspz1-2spätzle 1BiiNP_001138116.1NP_001153590QBOQ01001132.1 (9 × 10−8)
Toll pathwayspz2spätzle 2, neurotrophin 1NP_001261417.1XP_001948459.1QBOQ01001810.1 (4 × 10−20)
Toll pathwayspz3spätzle 3NP_609160.2XP_029341989.1QBOQ01001601.1 (2 × 10−21)
QBOQ01001203.1 (2 × 10−14)
QBOQ01000537.1 (9 × 10−13)
Toll pathwayspz4spätzle 4NP_609504.2NP_001153592QBOQ01000537.1 (2 × 10−28)
QBOQ01001203.1 (4 × 10−28)
QBOQ01001601.1 (2 × 10−14)
Toll pathwaySpz5spätzle 5NP_647753.1XP_001947495.2QBOQ01001431.1 (2 × 10−17)
Toll pathwayspz6spätzle 6NP_611961.1XP_001944046QBOQ01001423.1 (6 × 10−58)
Toll pathwayToll-1protein TollNP_524518.1XP_008182102.1QBOQ01001036.1 (3 × 10−57)
QBOQ01001364.1 (4 × 10−54)
QBOQ01000985.1 (8 × 10−48)
Toll pathwayToll-1protein TollNP_524518.1XP_001942733.2QBOQ01001036.1 (6 × 10−54)
QBOQ01001364.1 (7 × 10−49)
QBOQ01000985.1 (2 × 10−38)
Toll pathway18w18 wheeler, Toll-2NP_476814.1XP_001946943.2QBOQ01000059.1 (0.0)
QBOQ01000985.1 (0.0)
QBOQ01000048.1 (0.0)
Toll pathwayToll-6Toll-6NP_001246766.1XP_001947324.1QBOQ01000985.1 (0.0)
QBOQ01000059.1 (0.0)
QBOQ01000048.1 (0.0)
Toll pathwayToll-6Toll-6NP_001246766.1XP_003248960.1QBOQ01000048.1 (0.0)
QBOQ01000985.1 (0.0)
QBOQ01000059.1 (0.0)
Toll pathwayToll-7Toll-7NP_523797.1XP_001946943_2QBOQ01000059.1 (0.0)
QBOQ01000985.1 (0.0)
QBOQ01000048.1 (0.0)
Toll pathwayTolloTollo, Toll-8NP_524757.1XP_001948566.1QBOQ01000985.1 (0.0)
QBOQ01000059.1 (0.0)
QBOQ01000048.1 (0.0)
Toll pathwaytubtube, interleukin-1 receptor-associated kinase 4NP_001189164.1BAH72505.1QBOQ01000327.1 (8 × 10−15)
Toll pathwayMyd88myeloid differentiation primary response geneAAF58953.1XP_001948320.2not found
Toll pathwaypllpelleAAF56686.1XP_029346632.1QBOQ01000327.1 (2 × 10−33)
QBOQ01002061.1 (3 × 10−12)
QBOQ01001518.1 (3 × 10−11)
Toll pathwaycactcactusAAN10936.1NP_001156668.1not found
Toll pathwaycactincactinNP_523422.4XP_001952287.2QBOQ01001452.1 (1 × 10−88)
Toll pathwayPlipellinoNP_524466.1XP_001946282.3QBOQ01001351.1 (4 × 10−29)
Toll pathwayTraf1, Traf4TNF-receptor-associated factor 1AAD34346.1XP_001948355.1QBOQ01000448.1 (3 × 10−67)
Toll pathwayTraf2, Traf6TNF-receptor-associated factor 2AAF46338.1XP_029347356.1QBOQ01001366.1 (1 × 10−17)
Toll pathwayTraf3, Traf-likeTNF-receptor-associated factor 3NP_727976.1not foundnot found
Toll pathwaydldorsalAAF53611.1XP_001949498.2QBOQ01000587.1 (1 × 10−68)
Toll pathwayDifdorsal-related immunity factor, embryonic polarity proteinNP_523589.2XP_001949498.2QBOQ01000587.1 (1 × 10−68)
Jak/stat pathwaydomedomeless 1, interleukine JAK/STAT receptorCAD12503.1XP_029341085.1QBOQ01000913.1 (1 × 10−114)
Jak/stat pathwaydome2domeless 2Not foundXP_029341036.1QBOQ01000913.1 (6 × 10−111)
Jak/stat pathwayhops, jakhopscotch, Janus kinaseNP_511119.2XP_008188128.1QBOQ01001628.1 (2 × 10−29)
QBOQ01002061.1 (4 × 10−28)
QBOQ01000952.1 (6 × 10−24)
Jak/stat pathwayStat92Esignal-transducer and activator of transcription, marelleAAX33462.1XP_008188159.1QBOQ01001541.1 (5 × 10−44)
QBOQ01000405.1 (1 × 10−42)
Jak/stat pathwayupd1unpaired 1NP_525095.2not foundnot found
Jak/stat pathwayupd2unpaired 2NP_001356882.1not foundnot found
Jak/stat pathwayupd3unpaired 3NP_001097014.1not foundnot found
Imd pathwayimdimmune deficiencyNP_573394.1not foundnot found
Imd pathwaydFadddFaddNP_651006.1not foundnot found
Imd pathwayDredddeath related ced-3, caspase-1NP_477249.3XP_029344969.1QBOQ01001252.1 (4 × 10−61)
Imd pathwayRelRelishNP_477094.1not foundnot found
Imd pathwayTab2TAK1-associated binding protein 2NP_611408.2XP_003244590.1QBOQ01001392.1 (7 × 10−8)
Imd pathwayTak1TGF-β activated kinase 1AAF50895.1XP_029347425.1QBOQ01001920.1 (3 × 10−35)
QBOQ01001518.1 (6 × 10−29)
QBOQ01000779.1 (1 × 10−17)
Imd pathwaykeykennyNP_523856.2not foundnot found
Imd pathwayDiap2death-associated inhibitor of apoptosis 2NP_477127.1XP_016661891.1QBOQ01002166.1 (3 × 10−20)
QBOQ01001565.1 (2 × 10−12)
QBOQ01001600.1 (4 × 10−12)
Imd pathwayird5immune response deficiency 5, IK-β, IKKB, I-kappaB kinase betaNP_524751.3XP_001946184.1QBOQ01000860.1 (0.0)
Jnk pathwayhephemipterousNP_727661.1XP_008180171.1QBOQ01001476.1 (1 × 10−126)
Jnk pathwaybskbasketP92208.1XP_001945460.2QBOQ01001118.1 (6 × 10−45)
QBOQ01001795.1 (4 × 10−40)
QBOQ01001783.1 (4 × 10−26)
Jnk pathwayJraJun-related antigenAAF58845.1XP_001947556.1QBOQ01001648.1 (8 × 10−16)
Jnk pathwaykaykayakNP_001027579.1XP_016663984.1not found
Jnk pathwayegrEigerAAF58848.2XP_008178962.1QBOQ01000141.1 (7 × 10−5)
1 Alternative names are separated by commas. 2 E-values refer to A. pisum proteins, if present, or D. melanogaster if the gene was not found in A. pisum.
Table 3. Immune genes of Ericerus pela involved in response. Genes not found in E. pela are colored in red.
Table 3. Immune genes of Ericerus pela involved in response. Genes not found in E. pela are colored in red.
Role in Insect
Immunity
Gene
Symbol
Gene Name 1D. melanogaster
NCBI Protein ID
A. pisum
NCBI Protein ID
Best Matches (E-Value) 2
antimicrobial peptideAttattacinNP_523745.1not foundnot found
antimicrobial peptideCeccecropinC0HKQ7.1not foundnot found
antimicrobial peptideDefdefensinANY27112.1not foundnot found
antimicrobial peptideDrodrosocinXP_016946682.1not foundnot found
antimicrobial peptideMtkmetchnikowinAAO72489.1not foundnot found
antimicrobial peptide andropinP21663.1not foundnot found
antimicrobial peptide diptericinQER92349.1not foundnot found
antimicrobial peptideDrsdrosomycinANY27466.1not foundnot found
antimicrobial peptide holotricinXP_051861657.1not foundnot found
antimicrobial peptide bomaninA1ZB62.1not foundnot found
antimicrobialLOC100164856thaumatin-like proteinnot foundXP_001942718.2not found
antimicrobialLOC100160062thaumatin-like protein 1bnot foundXP_001942572.1not found
antimicrobialLOC100570639thaumatin-like protein 1not foundXP_003248856.4not found
antimicrobialLOC100162111uncharacterized LOC100162111, thaumatin familynot foundNP_001155516not found
antimicrobialLOC100168942TLP-PA-domain proteinnot foundNP_001156304.1not found
antimicrobialLOC100169496pathogenesis-related protein 5-likenot foundNP_001313585.1not found
microbial degradationLysXlysozyme X, i-typeCAL85493.1not foundnot found
microbial degradationLysBlysozyme B, i-typeNP_001261245.1not foundnot found
microbial degradationLysPlysozyme, i-typeNP_476828.1not foundnot found
microbial degradationLysClysozymeCAA80228not foundnot found
microbial degradationLysDlysozymeNP_476823.1not foundnot found
microbial degradationLysElysozymeNP_476827.2not foundnot found
microbial degradationLysSlysozymeNP_476829.1not foundnot found
microbial degradationlysozyme, i-typeLOC100167742, lysozymeACD99447.1XP_001949318.2QBOQ01000327.1 (6 × 10−19)
QBOQ01001156.1 (5 × 10−18)
QBOQ01000040.1 (7 × 10−18)
microbial degradationlysozyme, i-typeLOC100168424, destabilaseNP_611164.3NP_001156290.1QBOQ01002128.1 (5 × 10−17)
QBOQ01000327.1 (6 × 10−8)
QBOQ01000040.1 (9 × 10−8)
microbial degradationlysozyme, i-typeLOC100160909, destabilaseNP_611163.2NP_001155465.1QBOQ01000040.1 (3 × 10−11)
QBOQ01002128.1 (8 × 10−8)
fungal degradationCht2chitinase-like protein 2, mucinNP_001261282.1XP_016663378.1QBOQ01000205.1 (1 × 10−40)
QBOQ01000535.1 (1 × 10−19)
QBOQ01001282.1 (3 × 10−19)
fungal degradationCht4chitinase-like protein 4, flocculation proteinNP_524962.2XP_029343203.1QBOQ01001282.1 (3 × 10−40)
QBOQ01000205.1 (3 × 10−20)
QBOQ01000535.1 (5 × 10−19)
fungal degradationCht5chitinase-like protein 5, endochitinaseNP_650314.1XP_008181779.1QBOQ01000410.1 (3 × 10−47)
QBOQ01000205.1 (3 × 10−16)
QBOQ01001282.1 (5 × 10−21)
fungal degradationCht6chitinase-like protein 6, flocculation proteinNP_001245602.1XP_029343203.1QBOQ01001282.1 (3 × 10−40)
QBOQ01000205.1 (3 × 10−20)
QBOQ01000535.1 (5 × 10−19)
fungal degradationCht7chitinase-like protein 7, chitinase 10NP_647768.3XP_001950380.1QBOQ01000535.1 (9 × 10−87)
QBOQ01000062.1 (2 × 10−61)
QBOQ01001282.1 (3 × 10−24)
fungal degradationCht7chitinase 3-like, LOC100169240NP_647768.3XP_008182858.1QBOQ01002026.1 (6 × 10−19)
QBOQ01000535.1 (8 × 10−14)
QBOQ01000205.1 (4 × 10−8)
fungal degradationCht6LOC100162732NP_001245599.1XP_001945470.2QBOQ01001292.1 (2 × 10−7)
fungal degradationidgf6idgfNP_001286499.1NP_001162142.1QBOQ01000713.1 (3 × 10−55)
prophenoloxidase responsePPO1prophenoloxidase 1NP_476812.1XP_001949307.1QBOQ01000687.1 (5 × 10−45)
QBOQ01000496.1 (8 × 10−33)
QBOQ01002014.1 (3 × 10−27)
prophenoloxidase responsePPO2prophenoloxidase 2NP_610443.1XP_001951137.1QBOQ01000687.1 (2 × 10−43)
QBOQ01000496.1 (2 × 10−37)
QBOQ01002014.1 (5 × 10−27)
phenoloxidase activationPAF2, PPAF2phenoloxidase-activating factor 2AAO24923.1XP_003244500.1QBOQ01000870.1 (3 × 10−13)
QBOQ01001093.1 (2 × 10−12)
QBOQ01001364.1 (2 × 10−9)
phenoloxidase activationPAF2, PPAF2phenoloxidase-activating factor 2AAO24923.1XP_001952301.1QBOQ01000870.1 (2 × 10−24)
QBOQ01001093.1 (2 × 10−28)
QBOQ01002068.1 (1 × 10−15)
phenoloxidase activationSPserine protease-like precursorNP_001097766.1NP_001155379.1QBOQ01000877.1 (2 × 10−24)
QBOQ01000234.1 (2 × 10−23)
QBOQ01001364.1 (6 × 10−8)
cell aggregationHmct, hemolectinhemocytinNP_001261809.1XP_001952865.2QBOQ01002094.1 (3 × 10−22)
production of nitric oxide, a toxic gasNosnitric oxide synthaseNP_001027243.2XP_029343919.1QBOQ01001094.1 (4 × 10−29)
peptidoglycan degradationldcaputative LD carboxypeptidasenot foundXP_029341985.1QBOQ01000175.1 (2 × 10−43)
1 Alternative names are separated by commas. 2 E-values refer to A. pisum proteins, if present, or D. melanogaster if the gene was not found in A. pisum.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Becchimanzi, A.; Nicoletti, R.; Di Lelio, I.; Russo, E. Immune Gene Repertoire of Soft Scale Insects (Hemiptera: Coccidae). Int. J. Mol. Sci. 2024, 25, 4922. https://doi.org/10.3390/ijms25094922

AMA Style

Becchimanzi A, Nicoletti R, Di Lelio I, Russo E. Immune Gene Repertoire of Soft Scale Insects (Hemiptera: Coccidae). International Journal of Molecular Sciences. 2024; 25(9):4922. https://doi.org/10.3390/ijms25094922

Chicago/Turabian Style

Becchimanzi, Andrea, Rosario Nicoletti, Ilaria Di Lelio, and Elia Russo. 2024. "Immune Gene Repertoire of Soft Scale Insects (Hemiptera: Coccidae)" International Journal of Molecular Sciences 25, no. 9: 4922. https://doi.org/10.3390/ijms25094922

APA Style

Becchimanzi, A., Nicoletti, R., Di Lelio, I., & Russo, E. (2024). Immune Gene Repertoire of Soft Scale Insects (Hemiptera: Coccidae). International Journal of Molecular Sciences, 25(9), 4922. https://doi.org/10.3390/ijms25094922

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop