- Research article
- Open Access
Immune-related genes of the larval Holotrichia parallela in response to entomopathogenic nematodes Heterorhabditis beicherriana LF
BMC Genomics volume 22, Article number: 192 (2021)
Entomopathogenic nematodes (EPNs) emerge as compatible alternatives to conventional insecticides in controlling Holotrichia parallela larvae (Coleoptera: Scarabaeidae). However, the immune responses of H. parallela against EPNs infection remain unclear.
In present research, RNA-Seq was firstly performed. A total of 89,427 and 85,741 unigenes were achieved from the midgut of H. parallela larvae treated with Heterorhabditis beicherriana LF for 24 and 72 h, respectively; 2545 and 3156 unigenes were differentially regulated, respectively. Among those differentially expressed genes (DEGs), 74 were identified potentially related to the immune response. Notably, some immune-related genes, such as peptidoglycan recognition protein SC1 (PGRP-SC1), pro-phenoloxidase activating enzyme-I (PPAE-I) and glutathione s-transferase (GST), were induced at both treatment points. Bioinformatics analysis showed that PGRP-SC1, PPAE-I and GST were all involved in anti-parasitic immune process. Quantitative real-time PCR (qRT-PCR) results showed that the three immune-related genes were expressed in all developmental stages; PGRP-SC1 and PPAE-I had higher expressions in midgut and fat body, respectively, while GST exhibited high expression in both of them. Moreover, in vivo silencing of them resulted in increased susceptibility of H. parallela larvae to H. beicherriana LF.
These results suggest that H. parallela PGRP-SC1, PPAE-I and GST are involved in the immune responses to resist H. beicherriana LF infection. This study provides the first comprehensive transcriptome resource of H. parallela exposure to nematode challenge that will help to support further comparative studies on host-EPN interactions.
The dark black chafer (Holotrichia parallela Motschulsky) is a polyphagous pest that impacts many major crops, pastures and herbs in East Asia, particularly in China [1,2,3]. The larvae, white grubs, live in soil and prefer to feed on plant roots, causing an average economic loss of more than 15% per year ; in serious cases, the losses may exceed 50% . In recent years, the implementation of agricultural measures, such as no-tillage and shallow tillage systems, straw return, has created unique conditions for the survival and quantity of white grubs [6, 7]. The application of chemical insecticides has been widely used as an effective measure in white grub management but resulted in soil and groundwater pollution [8, 9].
Entomopathogenic nematodes (EPNs), including Steinernematidae and Heterorhabditidae, are regarded as potential biological control agents to control a wide range of insect pests, especially those that occur in soil . Similar to predators, the infective juveniles (IJs) of EPNs possess chemoreceptors that can actively search for susceptible hosts living in the soil [11, 12]. Once IJs have captured susceptible hosts, they will enter the host through their natural body openings such as stoma, valve and anus, eventually reaching the host’s haemocoel and gut region, and then the IJs release the symbiotic bacteria that secrete proteases, protoxins and other insecticidal substances complex to establish infection . More importantly, EPNs are highly virulent, killing the pest quickly, meanwhile they are safe to vertebrates and plants . In 1929, EPNs were first discovered as the parasites of Japanese beetle (Popillia japonica) in New Jersey  and produced a high mortality to white grubs . To date, approximately seven EPN species, including Steinernema scarabaei, Steinernema glaseri, Heterorhabditis bacteriophora H06, Heterorhabditis spp., Steinernema longicaudum X-7, Heterorhabditis indica LN2 and Heterorhabditis beicherriana, have shown potential control abilities for P. japonica, masked chafers (Cyclocephala borealis Arrow), Holotrichia oblita and H. parallela [7, 17,18,19,20,21]. However, the pathogenicity of nematodes is unstable in practical applications, and the control effect is unsatisfactory [22, 23]. The deadliest cause was that white grubs possess a strong immune system , which may effectively resist nematodes to establish infection . Specifically, white grubs can recognize EPNs through the sensitive immune system and can encapsulate or melanize invading nematodes [26,27,28].
Generally, the immune system of insects will be activated after infection by parasites or pathogens. First, pattern recognition proteins (PRPs) primarily recognize and tag invaders [29,30,31,32]. Then, PRPs rely on the cascade reaction of serine proteases (SPs) to continuously expand signals and eventually transmit the signals to the nucleus . Finally, antimicrobial peptides (AMPs) were used to perform phagocytosis, nodulation and encapsulation depending on the invader [34, 35], and other immune-related genes, including GSTs, heat shock protein (HSPs) and superoxide dismutase (SOD), were also involved in the anti-parasitic immune response .
Recently, mRNA sequencing (RNA-Seq or transcriptome sequencing) has been widely applied to identify immune-related genes and study the molecular basis of host-bacterial and/or host-parasite interactions. Moreover, this method also provides comprehensive insight into the immune gene repertoire of non-model insects, and this technology has become a reliable method to identify and analyze differentially expressed genes (DEGs) as targets for pest control [36, 37]. For instance, transcriptome sequencing results of Drosophila response to H. bacteriophora H222 infection indicated that most of the strongly induced genes were implicated in immune responses ; Many immune-related genes (such as PRPs, immune-related signal transduction proteins, AMPs and cellular response proteins) were induced in the white-spotted flower chafer (Protaetia brevitarsis seulensis Kolbe) with Escherichia coli and Saccharomyces cerevisiae challenge ; Heliothis virescens response to H. bacteriophora infection results indicated that insect immune response genes were induced upon nematode invasion, but the majority of these genes were suppressed after the release of symbiotic bacteria by the nematode . However, the molecular immune mechanism of H. parallela against EPNs infection has not been elucidated.
In this study, we primarily generated the first comprehensive transcriptome resource regarding H. parallela larvae in response to the nematode infection and screened the DEGs associated with immune defence. Then, we predicted the roles of three key immune-related genes in anti-parasite immune processes based on bioinformatics analysis. After that, we performed qRT-PCR to investigate the spatiotemporal expression patterns of these key immune-related genes. Finally, we silenced these immune-related genes to test their effects on the susceptibility of H. parallela larvae to H. beicherriana LF.
Transcriptome overview of Holotrichia parallela midgut
Through sequencing, a total of 89,427 and 85,741 unigenes were generated from the midgut samples of H. parallela larvae treated with the nematode H. beicherriana LF for 24 and 72 h, respectively. The specific statistics of transcriptomes are shown in Table 1.
After treating the larvae of H. parallela with nematodes for 24 h, 2545 genes were significantly differentially expressed compared to the sterile water treatment, among which 925 were upregulated and 1620 were downregulated. At 72 h postexposing, 3156 genes were significantly differentially expressed, among which 891 were upregulated and 2265 were downregulated (Fig. 1a). In addition, we found that 159 genes were upregulated and 333 genes were downregulated at both time- poins; 214 genes were downregulated at 24 h and upregulated at 72 h; 217 genes were upregulated at 24 h and downregulated at 72 h (Fig. 1b).
Of all the unigenes, approximately 14.90% were annotated to seven databases, including NT, NR, UniProt database, RNAMMER, eggNOG, KEGG and GO. Using the GO classification (Fig. 1c-d), these DEGs participated in at least 58 biological activities and were characterized into three groups: cellular component, biological process and molecular function. For the cellular component, the high percentage of genes were concentrated in the cell part category (24 h, up: 1242 genes, 17.46%; down: 3423 genes, 16.03%. 72 h, up: 1078 genes, 27.05%; down: 7244 genes, 18.51%). The biological process group showed a significant percentage of genes assigned to metabolic process category (24 h, up: 1148 genes, 16.14%; down: 2559 genes, 11.98%. 72 h, up: 1022 genes, 25.65%; down: 5490 genes, 14.03%). For molecular function, the most represented ontology was catalytic activity (24 h, up: 1240 genes, 17.44%; down: 2077 genes, 9.73%. 72 h, up: 996 genes, 24.99%; down: 4187 genes, 10.7%).
Immune defensive pathways
Considering that the immune response between H. parallela and H. beicherriana LF interaction is the key factor in the successful control of nematode parasitism, the transcripts related to immune defence responses upon nematode infection were selected. Through a complete search of the immune-related terms from the differential gene expression data, we identified a total of 74 DEGs that might be related to the immune defence mechanism in H. parallela. In this instance, differentially expressed transcripts were distributed among different immune response pathways, including recognition, activation of signalling pathways and production of effector molecules (Table 2). The results showed that most of the genes showed no significant expression change after treatment with nematodes for 24 h, and many genes showed depressed expression after treatment with nematodes for 72 h relative to the control. However, some transcripts were induced for both treatment points, such as PGRP-SC1, PPAE-I and GST, indicating that those genes may be involved in the immune responses to effectively resist nematode infection. The annotations of these genes in transcriptome data are summarized in Supplementary Table S1. These sequence data have been submitted to the GenBank databases under accession number DN15104, DN15190 and DN16733, respectively.
Bioinformatics analysis of PGRP-SC1, PPAE-I and GST
Based on the sequence analysis, PGRP-SC1 lacked a signal peptide but contained a transmembrane segment, which indicated that it may function on the cell membrane to recognize and bind pathogens. Moreover, PGRP-SC1 contains a highly conserved homologous PGRP domain consisting of approximately 130 amino acid residues at the shed terminal (Fig. 2a). It has also been noticed that two of five key residues responsible for zinc binding and amidase catalytic activity (His18, Tyr47, His123, Lys129 (Thr in D. melanogaster PGRP-LB/SC1/SC2) and Cys131 in T7 lysozyme  were substituted from His18 to Val191 and from Cys131 to Ser303 in H. parallela PGRP-SC1 (Fig. 2b). The results of Pro-CHECK showed that 3D model of PGRP-SC1, PPAE-I and GST were all reasonably constructed (Supplementary Fig. S1). Among them, the 3D structure of PGRP-SC1 indicated that three peripheral α-helices and five β-strands constitute the active domain center (Fig. 2c), and binding pocket was found in the protein surface (Fig. 2d). The phylogenetic analysis showed that H. parallela PGRP-SC1 (Hp SC1) clustered with O. taurus PGRP 2 (Ot 2) (Fig. 2e).
Cleavage activation of pro-phenoloxidase mediated by clip domain SPs is critical to melanization, which could originate from melanin to promote wound healing, epidermal sclerosis, free radicals, and the aggregation and encapsulation of pathogens . Sequence analysis showed that the open reading frame of H. parallela PPAE-I (a typical SP) contains 1095 nucleotides and 365 amino acids. We found that PPAE-I possessed a catalytic triangle composed of His155, Asp221 and Ser316. PROSITE analysis showed that PPAE-I was composed of a signal peptide (22 amino acid residues), a clip domain (51 amino acid residues) and a catalytic domain (255 amino acid residues). PPAE-I contains 12 extremely conserved cysteine residues and 6 in each of the two domains, which can form intramolecular disulphide bonds (Fig. 3a). In addition, the cutting site of PPAE-I was DEEK^ILGG, where it was recognized and activated by the upstream protease (Fig. 3b). The three highly conserved amino acids Asp310, Ser335 and G1y337 made up the active pocket of the PPAE-I substrate, which determines the specificity of PPAE-I substrate recognition (Fig. 3d). A calcium ion was anchored by Glu175 and Asp183, which may stabilize the overall structure of the PPAE-I domain (Fig. 3c). Further phylogenetic tree analysis showed that H. parallela PPAE-I (Hp PPAE-I) clustered close to Holotrichia diomphalia PPAF-I (Hd PPAF-I) with bootstrap values of 100 (Fig. 3e).
GSTs are the superfamily of multifunctional detoxification isoenzymes involved in the regulation of redox homeostasis and play crucial roles in innate immunity . Sequence analysis showed that H. parallela GST genes possessed N/C-terminal domains, G/H-binding sites and dimer interfaces (polypeptide binding sites) (Fig. 4a-b). The GST structure exhibited the overall folding of sigma-class. Each monomer of H. parallela GST included nine α-helices and four β-strands (Fig. 4c). One molecule of glutathione (GSH) bound to each H. parallela GST monomer and the binding site was located in the deep cleft between the two domains (Fig. 4d). Phylogenetic tree analysis indicated that most species were well differentiated with high bootstrap values, and H. parallela GST (Hp GST) clustered close to A. pisum GST (Api GST) (Fig. 4e).
Spatiotemporal analysis of PGRP-SC1, PPAE-I and GST
The expression patterns of PGRP-SC1, PPAE-I and GST genes were further determined in different developmental stages and tissues, respectively. The PGRP-SC1 was expressed throughout all developmental stages, exhibiting a low level at third instar larvae (Fig. 5a). PPAE-I exhibited a high level at first instar larvae (Fig. 5b). GST exhibited a high level at third instar larvae (Fig. 5c). For different tissues, PGRP-SC1 and PPAE-I had higher expressions in midgut and fat body, respectively, while GST exhibited high expression in both of them (Fig. 5d-f).
To validate the transcriptome results, qRT-PCR was performed. The qRT-PCR data shown the PGRP-SC1, PPAE-I and GST gene had a higher expression in the larvae of H. parallela treated with nematodes for 24 h with 4.66-, 3.37- and 4.40-fold, respectively, compared to the larvae treated with sterile water. After the larvae treated with nematodes for 72 h, the fold changes were 2.92-, 1.06- and 3.00-fold, respectively (Fig. 6). The results obtained through qRT-PCR were in agreement with the transcriptome data on the fold changes of three candidate immune genes (R2 = 0.847, P = 0.194 for 24 h and R2 = 0.909, P = 0.256 for 72 h at 95% confidence interval). The consistency of qRT-PCR results with transcriptome data confirmed the reliability and accuracy of sequencing results.
The expression of PGRP-SC1, PPAE-I and GST after RNAi
RNAi was used to study the functions and their effects on nematode infection resistance of PGRP-SC1, PPAE-I and GST. qRT-PCR confirmed that we successfully knocked down the genes following injection of dsRNAs into the second instar larvae of H. parallela (Fig. 7). The expression difference were not significant between control groups (dsGFP injected and water-injected) of those genes in different treatments. However, the expression of PGRP-SC1 and GST in the treatments (dsPGRP-SC1 and dsGST injected larvae) was both significantly reduced in comparison to the water injected larvae at 48 h. The relative expression levels of these two genes were depressed by 92.18% (PGRP-SC1) and 77.00% (GST), respectively. The expression of PPAE-I in dsPPAE-I injected larvae was significantly depressed in comparison to the water injected larvae at 24, 48 and 72 h. The relative expression levels of PPAE-I were depressed by 95.83, 99.10 and 66.33%, respectively.
Susceptibility to Heterorhabditis beicherriana LF after RNAi
The susceptibility of H. parallela second instar larvae was investigated after we silenced PGRP-SC1, PPAE-I and GST (Fig. 8). After injection of dsPGRP-SC1, dsPPAE-I or dsGST for 48 h, the larvae were exposed to nematodes (100 and 200 IJs/grub)-incorporated diets for 24 h, and the larvae exhibited a higher mortality than the control larvae injected with dsGFP and water. While dsGFP did not show significant differences to water control, the mortality of H. parallela larvae injected with dsPGRP-SC1 or dsPPAE-I was significantly higher than water injected control by increased 41.66% or 33.33% when exposed to nematodes at 100 IJs/grub and 24.08% or 25.00% at 200 IJs/grub. The mortality of H. parallela larvae was significantly different between the injection of dsGST and water injected control by increased 16.67% when exposed to nematodes at 100 IJs/grub, whereas the larvae only exhibited a slightly higher mortality by increased 12.97% at 200 IJs/grub.
As one of the most important destructive belowground herbivores in China, the larvae of H. parallela have caused considerable damage to crop production . To date, many studies have proven that EPNs are compatible alternatives to conventional insecticides for the larvae of H. parallela because they not only can actively search for their hosts in soil but also possess the property of fast-acting than other biological agents [12, 43]. Meanwhile, combating to the host immune response are regarded as crucial to successful parasitism . Thus, the immune and anti-immune mechanisms between H. parallela and EPNs are urgently to clarify. In this study, we carried out RNA-Seq to screen the DEGs involved in the immune response of H. parallela against EPNs infection and functionally validated three crucial immune-related genes using in vivo RNAi technology.
Through comparative analysis, 925 transcripts were induced upon nematode challenge for 24 h in our grubs, whereas only 314 genes were induced in H. virescens larvae after exposure to H. bacteriophora for 23 h . However, homologue of thioester-containing complement protein 3, basement membrane component (glutactin) and recognition protein (GNBP-like 3) which were not identified in the 74 DEGs implicating in immune responses of H. parallela has been proven to be involved in the immune responses of Drosophila larvae against H. bacteriophora (H222 strain) infection . Therefore, we think the immune responses between white grubs and Drosophila against nematodes infection were different. In addition, some sequences respond to symbiotic bacterium were detected after treatment with nematodes for 72 h (Supplementary Table S2). Previous studies have shown that 10% of H. bacteriophora can escape the encapsulation of P. japonica larvae after entering for 24 h  and the escaped nematodes will release symbiotic bacteria 30 min later . Thus, the symbiotic bacteria of H. beicherriana LF was predicted to participate in immune responses after H. parallela larvae treated with the nematode for 72 h.
Identification of parasite/pathogen invasion in the host is the initial step in the initiation of immune responses, and this process relies on biosensor proteins called pattern recognition receptors (PRRs), which can detect and bind to certain pathogen-associated molecular patterns (PAMPs) on the surface of invading microbes or foreign bodies [46, 47]. PGRPs are important members of insect PRRs and play an important role in innate immunity [48, 49]. Previous studies have shown that PGRP-LB1, −LB2, −LC2, −LC3 and -LD in Monochamus alternatus are likely to encode membrane proteins ; PGRP-SC1 of H. parallela contains a similar structure (without a signal peptide but containing a transmembrane segment), indicating that it may function on the cell membrane to recognize and trigger immune responses. Moreover, 2 of 5 key residues responsible for zinc binding and amidase catalytic activity were substituted, indicating that H. parallela PGRP-SC1 has no amidase activity and cannot cleave the bacterial peptidoglycan (PGN) . However, this type of PGRP may function as a sensor for ligand-dependent signalling to activate the pro-phenoloxidase (PPO) cascade . For example, overexpression of PGRP-LE in D. melanogaster led to the activation of the PPO cascade . Therefore, we speculate that H. parallela PGRP-SC1 may participate in the PPO cascade reaction as a signal recognition receptor. Further phylogenetic analysis of PGRP showed that H. parallela PGRP-SC1 and B. mori PGRP-L1 were in a cluster. It is already known that PGRP-L1 in B. mori could induce the expression of antigen processing machinery (APM) genes and may be players in the immune deficiency (IMD) pathway . Therefore, we speculate that H. parallela PGRP-SC1 may also be involved in the IMD pathway. However, the detailed functionality of PGRP-SC1 in H. parallela merits further studies with more biochemical methods.
Once the insect PRRs recognize the pathogen, the second step of innate immune-signal modulation is immediately initiated. At this time, the cascade reaction of SPs mainly relies on continuously expanding signals and eventually transmitting the signals to the nucleus . Our research indicated that H. parallela PPAE-I was a typical serine protease with a clip domain and a catalytic domain, and further phylogenetic tree analysis showed that H. parallela PPAE-I clustered close to H. diomphalia PPAF-I. Previous findings have demonstrated that clip domain SPs are actively involved in the prophenoloxidase-activation cascade [56, 57]. Moreover, H. diomphalia PPAF-I could activate H. diomphalia PPO in the presence of PPAF-II . In summary, H. parallela PPAE-I may function in inducing pro-phenoloxidase activity and mediating the melanization response.
GST is a widespread multigene family of detoxification enzymes and acts as an effector molecule after induction by immunization in innate immunity . GST has a variety of functions, including detoxification of endogenous and/or xenobiotic compounds and protection of cells from oxidative damage, whereas other functions may include contributing to metabolic and signalling pathways [60, 61]. A number of studies have shown that GST is generally induced by immunization and acts as an effector to regulate the metabolism of exogenous toxic substances. For instance, T. castaneum GST was activated to defend against lipopolysaccharide (LPS)-challenge , and P. brevitarsis GST was upregulated in response to E. coli challenge . In addition, sequence analysis showed that H. parallela GST possessed the characteristics of the GST family, and phylogenetic analysis indicated that H. parallela GST and D. melanogaster/N. vespilloides GST were in a cluster, which has been proven to respond to microbial challenge [63, 64]. In general terms, H. parallela GST may play a partial immune protection role in H. parallela against nematode infection.
Then, the roles of these three genes were determined by dsRNA injection-mediated RNAi and pathogenicity assays. The results showed that silencing of PGRP-SC1 and PPAE-I resulted in significant susceptibility of H. parallela to nematodes, and the mortality of H. parallela larvae was increased to different degrees compared to that of the control (water-injected). In addition, silencing of GST also resulted in increased susceptibility of H. parallela to nematodes. These results suggest that PGRP-SC1, PPAE-I, and GST are involved in immune responses to resist nematode infection.
Overall, our study strongly support that PGRP-SC1, PPAE-I and GST are genes related to nematode resistance in H. parallela, and the detailed functionalities of these genes deserve further investigation. Meanwhile, these results may help to elucidate the nematode resistance mechanism and provide valuable insights in designing appropriate resistance management strategies for this famous agricultural pest. In addition, more studies are needed to elucidate the immune mechanisms between H. parallela and symbiotic bacteria of H. beicherriana LF or other kinds of EPNs.
Nematode and insect culture
Entomopathogenic nematode H. beicherriana LF were cultured in vivo on larvae of the greater wax moth (Galleria mellonella) at 25 °C and 80% relative humidity (RH). Released infective juveniles were collected and stored at 7 ± 1 °C for no more than 15 days before use . The larvae of H. parallela were fed germinated wheat seeds, and actively feeding grubs (approximately 13 days old and at the beginning of second instar) were kept individually in the holes of six-well bioassay trays (Corning Inc., NY, USA) filled with 15 g and 18% soil moisture (w/w) of high temperature sterilized sandy loam soil (76% sand, 16% silt, 6% clay, 0.5% organic matter; pH 6.7) for use.
Nematode infection to Holotrichia parallela
One hundred IJs were added into each hole of the six-well bioassay trays described above. A total of 60 H. parallela larvae were treated with nematodes, and control group was treated with sterile water. The six-well bioassay trays were then transferred into a dark incubator under the following conditions: 25 °C and 80% RH. The midgut samples (n = 3 for each treatment) of H. parallela used for RNA-Seq were collected in RNase-free tubes at 24 h and 72 h after treatment with nematodes or sterile water following the method below.
Disintegration of Holotrichia parallela midgut
H. parallela larvae treated with nematodes or sterile water were transferred on ice for 1 min and then placed on a precooled petri dish (d = 15 cm). After adjusting the focal length of SMZ-168 Stereo Zoom microscope (Sinolab Inc., Beijing, China) until the insect body could be observed clearly, the head and tail of the larvae were cut off with scissors and the whole midgut was carefully pulled out using tweezers. The malpighian tubule and fat around the midgut were removed with tweezers, and the peritrophic membrane was pulled to remove all food particles in the midgut. After three rinses with RNase-free water, three samples for each treatment were collected in a 1.5 ml RNase-free tube and stored at − 80 °C for RNA isolation. The experiment was performed with two replicates of each treatment for RNA-Seq.
RNA isolation and cDNA library preparation
Total RNA was extracted using TRIzol Reagent (Invitrogen, Carlsbad, CA, USA) following the manufacturer’s protocol followed by ethanol/isopropanol precipitation and suspended in a final volume of 30 μl with RNase-free water. After extraction, the total RNA quality of different treatments were measured with NanoDrop 2000 spectrophotometer (Thermo Scientific) and in a 1.5% agarose gel. RNA integrity and concentration were subsequently analyzed using 2100 RNA Nano 6000 Assay Kit (Agilent Technologies, CA, USA).
Following the Illumina protocols, 2 μg RNA of different treatment was used for cDNA library construction, and each treatment was replicated two times. Following the use of oligo (dT) magnetic beads, fragmentation of mRNA, generation of the first strand cDNA with a random hexamer primer and synthesis of second strands, cDNA end repair and adenylation at the 3′ end, and adapter ligation and cDNA fragment enrichment were performed. Then, the cDNA products were purified and amplified by PCR to create the final cDNA library.
RNA-Seq and de novo assembly
Paired-end sequencing was conducted on an Illumina HiSeq 2000 (Illumina Inc., USA) at Anoroad Genomics Co., Ltd. (Beijing, China). After sequencing, the raw reads were filtered to generate clean reads. First, the reads with adapters were removed, and then the reads with > 5% unknown bases (N) were discarded. Finally, the low-quality reads (with > 50% of nucleotides for which the Phred Quality Score Q was less than or equal to 19) were filtered out.
The clean reads of each treatment were loaded to Trinity (Trinity Release v2.4.0)  for de novo assembly under the paired-end mode with default parameters. Based on the filtered clean data, the full-length transcript sequence was assembled with Trinity, and each gene was taken based on the transcript sequence. The longest transcript sequence was recorded as ‘Unigene’.
RPKMs (Reads Per Kilobase Million Mapped Reads)  were used in DESeq2  to compare the differences of gene expression between the nematodes and sterile water treat. Genes with |log2 Fold Change| ≥ 1 and q < 0.05 were selected as significantly DEGs. Then, DEGs were further annotated by Gene Ontology (GO) function (http://www.geneontology.org/).
Identification of immune-related genes in Holotrichia parallela
To identify the immune-related genes involved in the process of H. parallela larvae defence against nematode infection, we query with the terms that are reported as the immune relative genes from the differential gene expression data that obtained above. The terms mainly include ‘immune’, ‘recognition protein’, ‘haemocyte’, ‘glutactin’, ‘antimicrobial’, ‘lysozyme’, ‘lipopolysaccharide’, ‘lectins’, ‘peptidoglycans’ etc. All candidate immune genes were manually verified with the BLASTx program at the National Center for Biotechnology Information (NCBI: https://www.ncbi.nlm.nih.gov/). Thereafter, we focused on immune genes that were significantly induced at both time points, as these potential candidate immune-related genes may be involved in the immune responses to effectively resist nematode infection.
Bioinformatics analysis of three potential candidate immune-related genes
The deduced amino acid sequences of three potential candidate immune-related genes (PGRP-SC1, PPAE-I and GST) were obtained using the Translate tool provided by the Swiss Institute Bioinformatics (https://web.expasy.org/translate/). The signal peptide and transmembrane domain were analyzed by SignalP 5.0 (http://www.cbs.dtu.dk/services/SignalP) and TMHMM server v. 2.0 (http://www.cbs.dtu.dk/services/TMHMM/). In addition, domain analysis of the retrieved protein sequences were executed by NCBI Conserve-Domain Tool (https://www.ncbi.nlm.nih.gov/cdd) and PROSITE (http://au.expasy.org/prosite/). The 3D structures were obtained using templates of Bumblebee PGRP-SA (5xz4.1.A, 34.76% identity), Holotrichia diomphalia PPAF-I (2olg.1.A, 66.67% identity) and Nilaparvata lugens GST S2 (5h5l.1.A, 40.30% identity), respectively, employing online Swiss-model software. Pro-CHECK (http://servicesn.mbi.ucla.edu/PROCHECK/) was used to evaluate the rationality of the 3D structure of the constructed model. The generated PDB files were then visualized with Pymol Molecular Graphics System.
Phylogenetic analysis of PGRP-SC1 from H. parallela and 53 PGRPs from 12 other insect species; PPAE-I from H. parallela and 26 SPs from 12 other insect species; GST from H. parallela and 26 GSTs from 12 other insect species, which were used to construct a phylogenetic tree using MEGA5.0. Phylogenetic analysis were conducted by the neighbour-joining method with p-distance under the default parameters. Bootstrap values were obtained by the bootstrap method using 1000 repetitions. The species name, gene name, abbreviation and GenBank accession involved in phylogenetic analysis are listed in Supplementary Table S3.
Spatiotemporal analysis of PGRP-SC1, PPAE-I and GST
To determine the expression profiles of these genes in different development stages and tissues, we collected eggs (5 days old, ~ 60 mg), first instar larvae (10 days old, ~ 60 mg), second instar larvae (18 days old), third instar larvae (28 days old) and various tissues (hemolymph, malpighian tubule, midgut and fat body) of the first instar larvae for RNA isolation. Total RNA samples of each treatment were individually extracted using TRIzol Reagent (TIANGEN, Biotech (Beijing) Co., Ltd., China) and each treatment repeat three times. First-strand cDNA was synthesized from 1 μg of total RNA following the instructions of TransScript One-Step gDNA Removal and cDNA Synthesis SuperMix (TransGen Biotech (Beijing) Co., Ltd., China). The cDNA products were diluted 10-fold for use as template in qRT-PCR.
qRT-PCR was performed with Tip Green qPCR SuperMix (TransGen, China) on an ABI Prism 7500 Fast Detection System using a GO Taq®qPCR kit (Invitrogen, USA) according to the manufacturer’s instructions. Gene-specific primers designed for target genes are listed in Supplementary Table S4. qRT-PCR reactions were performed in a 20 μl mixture containing 2 μl of cDNA, 10 μl of Tip Green qPCR SuperMix, 0.4 μl of each primer (10 μM), and 7.2 μl of nuclease-free water. The optimized qRT-PCR programme conditions were 94 °C for 30 s followed by 40 cycles of 94 °C for 5 s and 60 °C for 34 s. After the cycling protocol, melting curves were obtained by increasing the temperature from 60 to 95 °C (0.2 °C sec− 1) to denature the double-stranded DNA. H. parallela DAPDH was used as an internal standard to normalize the expression level. Each qRT-PCR experiment was performed with three biological replicates and four technical replicates. The relative expression level of target genes was calculated with the 2-ΔΔCt method .
To confirm the transcriptome results, these immune-related genes (i.e., PGRP-SC1, PPAE-I and GST) of H. parallela responding to nematode infection were selected for qRT-PCR analysis. Fresh total RNA was obtained as described in RNA-Seq and qRT-PCR was performed as described above.
RNA interference to Holotrichia parallela larvae
dsRNA synthesis: RNAi was used to study the functions of PGRP-SC1, PPAE-I, and GST involved in the immune defence system of H. parallela against H. beicherriana LF infection. Specific primers (Supplementary Table S4) with T7 promoter sequences were designed to synthesize dsRNA of PGRP-SC1, PPAE-I, GST and GFP using the T7 RiboMAX™ Express RNAi kit (Promega, WI, USA) according to the manufacturer’s instructions. The size of the dsRNA products was confirmed by electrophoresis on 1.2% agarose gel, and the final concentration of dsRNA was adjusted to 1 μg/μl.
In vivo RNAi: For in vivo RNAi, 1.5 μg of each dsRNA, dsGFP and the same volume of RNase-free water were injected into the haemocoel of H. parallela second instar larvae (18–20 days old) through the abdomen (WPI microinjection system, USA). After injection, the larvae were kept individually in the holes of six-well bioassay trays and fed germinated wheat seeds. The six-well bioassay trays were transferred into a dark incubator under the following conditions: 25 °C and 80% RH. The midgut samples (n = 3 for each treatment) from challenged or control groups (dsGFP-injected and water-injected) was collected after treatment for 24, 48 and 72 h, respectively, and the experiment was performed with three biological replicates. The qRT-PCR process was executed as described above.
Susceptibility of Holotrichia parallela larvae to Heterorhabditis beicherriana LF after RNAi
Susceptibility of H. parallela larvae to H. beicherriana LF after RNAi of PGRP-SC1, PPAE-I and GST was conducted by determining the pathogenicity change of H. beicherriana LF on H. parallela larvae. After injection following our previous procedures, the larvae (18–20 days old) were kept individually in the holes of six-well bioassay trays, fed germinated wheat seeds and then transferred into a dark incubator as described above. Forty eight hours later, two concentrations (100 and 200 IJs/grub) of H. beicherriana LF were added to each hole of six-well bioassay trays, respectively. Each injection treatment consisted of 36 individuals (grub) and the experiment was performed three times. The grub mortality was monitored after 24 h, and the infected grubs were marked as inactive and red.
Data statistics were performed using GraphPad Prism software version 8.0. All proportions were transformed by arcsine square root transformation in Microsoft Excel 5.0 before analysis . The significance levels of spatiotemporal analysis were analyzed using one-way analysis of variance (ANOVA) followed by Tukey significant difference test (P < 0.05). The significance levels of in vivo RNAi and white grubs mortality were analyzed using multiple t-test at the significance levels set at *P < 0.05, **P < 0.01, ***P < 0.001, NS, not significant.
Availability of data and materials
These sequence data of H. parallela PGRP-SC1, PPAE-I and GST have been submitted to the GenBank databases under accession number DN15104, DN15190 and DN16733, respectively.
Analysis of variance
Antigen processing machinery
Cellular component organization or biogenesis
Differentially expressed genes
Heat shock protein
Immune deficiency pathway
National Center for Biotechnology Information
Pathogen-associated molecular patterns
Peptidoglycan recognition protein
Peptidoglycan recognition protein SC1
Pro-phenoloxidase activating enzyme-I
Presynaptic process involved in chemical synaptic transmission
Pattern recognition proteins
Pattern recognition receptors
Quantitative real-time PCR
Reads Per Kilobase Million Mapped Reads
Toepfer S, Li HM, Pak SG, Son KM, Ryang YS, Kang SI, Han RC, Holmes K. Soil insect pests of cold temperate zones of East Asia, including DPR Korea: a review. J Pest Sci. 2014;87(4):567–95. https://doi.org/10.1007/s10340-013-0540-8.
Cao YZ, Wu YQ. Soil dwelling insect pests. In: Guo YY, editor. Crop diseases and insects in China. Beijing: China Agriculture Press; 2015. p. 1541–3.
Pan Y, Zhao SW, Tang XL, Wang S, Wang X, Zhang XX, Zhou JJ, Xi JH. Transcriptome analysis of maize reveals potential key genes involved in the response to belowground herbivore Holotrichia parallela larvae feeding. Genome. 2020;63(1):1–12. https://doi.org/10.1139/gen-2019-0043.
Luo ZX, Li KB, Cao YZ, Yin J, Zhang J, Zhang JH, Shang GQ. Investigations on soil-inhabiting pests in peanut fields in Henan. Plant Prot. 2009;35(2):104–8.
Liu QZ, Li JX, Xu XJ, Sun CM, Kang YG, Zhou HY, Hu DX, Ma J, Li SX. The preliminary study on grub control with Rhabditis (Oscheius) spp in peanut fields. Acta Agric Boreali-Sin. 2007;22(S2):250–3.
Pei GY, Ma SF, Liu J, Liu BC. Effects of different cultivation patterns on grubs occurrence and yield in soybean fields. Soybean Bull. 2010;4:19–20.
Li ET, Cao YZ, Zhang S, Li KB, Li JY, Wang QL, Liu CQ, Yin J. Bioefficacy of the combined application of entomopathogenic nematode Heterorhabditis beicherriana strain LF and Bacillus thuringiensis strain HBF-18 against Holotrichia oblita (Coleoptera: Scarabaeidae) larvae. Acta Entomol Sin. 2019;62(5):605–14.
Liu SS, Li KB, Yin J, Cao YZ. Review of the researches on biological control of grubs. Chin J Biol Control. 2008;24:168–73.
Li ET. Overview of the effectiveness of entomopathogenic nematodes. In: Chen WQ, editor. Green plant protection and rural revitalization. Beijing: China Agricultural Sciences and Technology Press; 2018. p. 184–91.
Grewal PS, Ehlers R-U, Shapiro-Ilan DI. Nematodes as biocontrol agents. In: Lichtfouse E, editor. Sociology, organic farming, climate change and soil science. Oxfordshire: CABI Press; 2005. p. 347–78.
Toepfer S, Kurtz B, Kuhlmann U. Influence of soil on the efficacy of entomopathogenic nematodes in reducing Diabrotica virgifera in maize. J Pest Sci. 2010;83(3):257–64. https://doi.org/10.1007/s10340-010-0293-6.
Bal HK, Grewal PS. Lateral dispersal and foraging behavior of entomopathogenic nematodes in the absence and presence of mobile and non-mobile hosts. PLoS One. 2015;10(6):e0129887. https://doi.org/10.1371/journal.pone.0129887.
Toubarro D, Lucena-Robles M, Nascimento G, Costa G, Montiel R, Coelho AV, Simōes N. An apoptosis inducing serine protease secreted by the entomopathogenic nematode Steinernema carpocapsae. Int J Parasitol. 2009;39(12):1319–30. https://doi.org/10.1016/j.ijpara.2009.04.013.
Kaya HK. Entomopathogenic nematodes in biological control. In: Gaugler R, Kaya HK, editors. Soil ecology. Boca Raton: CRC Press; 1990. p. 93–115.
Glaser RW, Fox H. A nematode parasite of the Japanese beetle (Popillia japonica Newm.). Science. 1930;71(1827):16–7. https://doi.org/10.1126/science.71.1827.16-b.
Glaser RW, Farrell CC. Field experiments with the Japanese beetle and its nematode parasite. J New York Entomol. 1935;43(3):345–71.
Li SC, Ding WD, Han FS. Application of entomophagous nematode Taishan no. 1 against peanut grubs. J Plant Prot. 1993;20:55–60.
Koppenhöfer AM, Brown IM, Gaugler R, Grewal PS, Kaya HK, Klein MG. Synergism of entomopathogenic nematodes and imidacloprid against white grubs: greenhouse and field evaluation. Biol Control. 2000;19(3):245–51. https://doi.org/10.1006/bcon.2000.0863.
Koppenhöfer AM, Fuzy EM. Steinernema scarabaei for the control of white grubs. Biol Control. 2003;28(1):47–59. https://doi.org/10.1016/S1049-9644(03)00048-3.
Guo WX, Yan X, Zhao GY, Han RC. Efficacy of entomopathogenic Steinernema and Heterorhabditis nematodes against white grubs (Coleoptera: Scarabaeidae) in peanut fields. J Econ Entomol. 2013;106(3):1112–7. https://doi.org/10.1603/EC12477.
Guo WX, Yan X. Zhao GY, Chen JH, Han RC. Efficacy of entomopathogenic Steinernema and Heterorhabditis nematodes against Holotrichia oblita. J Pest Sci 2015; 88(2): 359–368, DOI: https://doi.org/10.1007/s10340-014-0626-y.
Shapiro-Ilan DI, Gouge DH, Koppenhöfer AM. Entomopathogenic nematology. In: Gaugler R, editor. Factor affecting commercial success: Case studies in cotton turf and citrus. Oxfordshire: CABI Press; 2002. p. 333–56.
Shapiro-Ilan DI, Gough DH, Piggott SJ, Patterson JF. Application technology and environmental considerations for use of entomopathogenic nematodes in biological control. Biol Control. 2006;38(1):124–33. https://doi.org/10.1016/j.biocontrol.2005.09.005.
Bang K, Sj H, Lee J, Cho S. Identification of immunity-related genes in the larvae of Protaetia brevitarsis seulensis (Coleoptera: Cetoniidae) by a next-generation sequencing-based transcriptome analysis. J Insect Sci. 2015;15(1):142.
Castillo JC, Reynolds SE, Eleftherianos I. Insect immune responses to nematode parasites. Trends Parasitol. 2011;27(12):537–47. https://doi.org/10.1016/j.pt.2011.09.001.
Wang Y, Gaugler R, Cui L. Variations in immune response of Popillia Japonica and Acheta domesticus to Heterorhabditis bacteriophora and Steinernema species. J Nematol. 1994;26(1):11–8.
Li X-Y, Cowles RS, Cowles EA, Gaugler R, Cox-Foster DL. Relationship between the successful infection by entomopathogenic nematodes and the host immune response. Int J Parasitol. 2007;37(3–4):365–74. https://doi.org/10.1016/j.ijpara.2006.08.009.
An R, Voss M, Jagdale GB, Grewal PS. Differences in immune defense evasion of selected inbred lines of Heterorhabditis bacteriophora in two white grub species. Insects. 2012;3(2):378–89. https://doi.org/10.3390/insects3020378.
Janeway CA, Medzhitov R. Innate immune recognition. Annu Rev Immunol. 2002;20(1):197–216. https://doi.org/10.1146/annurev.immunol.20.083001.084359.
Steiner H. Peptidoglycan recognition proteins: on and off switches for innate immunity. Immunol Rev. 2004;198(1):83–96. https://doi.org/10.1111/j.0105-2896.2004.0120.x.
Montaño AM, Tsujino F, Takahata N, Satta Y. Evolutionary origin of peptidoglycan recognition proteins in vertebrate innate immune system. BMC Evol Biol. 2011;11(1):79. https://doi.org/10.1186/1471-2148-11-79.
Dziarski R, Gupta D. How innate immunity proteins kill bacteria and why they are not prone to resistance. Curr Genet. 2018;64(1):125–9. https://doi.org/10.1007/s00294-017-0737-0.
Cerenius L, Söderhäll K. The prophenoloxidase activating system in vertebrates. Immunol Rev. 2004;198(1):116–26. https://doi.org/10.1111/j.0105-2896.2004.00116.x.
Bang K, Park S, Yoo JY, Cho S. Characterization and expression of attacin, an antibacterial protein-encoding gene, from the beet armyworm, Spodoptera exigua (Hübner) (Insecta: Lepidoptera: Noctuidae). Mol Biol Rep. 2012;39(5):5151–9. https://doi.org/10.1007/s11033-011-1311-3.
Kwon HJ, Bang K, Lee M, Cho S. Molecular cloning and characterization of a lysozyme cDNA from the mole cricket Gryllotalpa orientalis (Orthoptera: Gryllotalpidae). Mol Biol Rep. 2014;41(9):5745–54. https://doi.org/10.1007/s11033-014-3446-5.
Liu Y, Shen DX, Zhou F, Wang GR, An CJ. Identification of immunity-related genes in Ostrinia furnacalis against entomopathogenic fungi by RNA-seq analysis. PLoS One. 2014;9(1):e86436. https://doi.org/10.1371/journal.pone.0086436.
Vogel H, Badapanda C, Knorr E, Vilcinskas A. RNA-sequencing analysis reveals abundant developmental stage-specific and immunity-related genes in the pollen beetle Meligethes aeneus. Insect Mol Biol. 2014;23(1):98–112. https://doi.org/10.1111/imb.12067.
Arefin B, Kucerova L, Dobes P, Markus R, Strnad H, Wang Z, Hyrsl P, Zurovec M, Theopold U. Genome-wide transcriptional analysis of Drosophila larvae infected by entomopathogenic nematodes shows involvement of complement, recognition and extracellular matrix proteins. J Innate Immun. 2014;6(2):192–204. https://doi.org/10.1159/000353734.
An R, Suri KS, Jurat-Fuentes JL, Grewal PS. Dynamics of transcriptomic response to infection by the nematode Heterorhabditis bacteriophora and its bacterial symbiont Photorhabdus temperata in Heliothis virescens larvae. Insect Mol Biol. 2017;26(5):584–600. https://doi.org/10.1111/imb.12321.
Dziarski R, Gupta D. The peptidoglycan recognition proteins (PGRPs). Genome Biol. 2006;7(8):232. https://doi.org/10.1186/gb-2006-7-8-232.
Zou Z, Shin SW, Alvarez KS, Kokoza V, Raikhel AS. Distinct melanization pathways in the mosquito Aedes aegypti. Immun. 2010;32(1):41–53. https://doi.org/10.1016/j.immuni.2009.11.011.
Wang CY, Zhao JM, Mu CK, Wang Q, Wu HF, Wang CL. cDNA cloning and mRNA expression of four glutathione s-transferase (GST) genes from Mytilus galloprovincialis. Fish Shellfish Immun. 2013;34(20):697e703.
Wu SH, Roger R, Loke TK, Curt AL, Douglas GP. Interaction between entomopathogenic nematodes and entomopathogenic fungi applied to third instar southern masked chafer white grubs, Cyclocephala lurida (Coleoptera: Scarabaeidae), under laboratory and greenhouse conditions. Biol Control. 2014;76:65–73. https://doi.org/10.1016/j.biocontrol.2014.05.002.
Wang Y, Campbell JF, Gaugler R. Infection of entomopathogenic nematodes Steinernema glaseri and Heterorhabditis bacteriophora against Popillia japonica (Coleoptera, Scarabaeidae) larvae. J Invertebr Pathol. 1995;66(2):178–84. https://doi.org/10.1006/jipa.1995.1081.
Bowen D, Rocheleau TA, Blackburn M, Andreev O, Golubeva E, Bhartia R, Ffrench-Constant RH. Insecticidal toxins from the bacterium Photorhabdus luminescens. Science. 1998;280(5372):2129–32. https://doi.org/10.1126/science.280.5372.2129.
Stokes BA, Yadav S, Shokal U, Smith LC, Eleftherianos I. Bacterial and fungal pattern recognition receptors in homologous innate signaling pathways of insects and mammals. Front Microbiol. 2015;6:1–12.
Meng Q, Zhang JH, Zhang H, Zhou GL, Ni RY, Zhao YN, Qin QL, Zou Z. Comparative analysis of C-type lectin domain proteins in the ghost moth, Thitarodes xiaojinensis (Lepidoptera: Hepialidae). Insect Sci. 2019;26(3):453–65. https://doi.org/10.1111/1744-7917.12564.
Sumathipala N, Jiang HB. Involvement of Manduca sexta peptidoglycan recognition protein-1 in the recognition of bacteria and activation of prophenoloxidase system. Insect Biochem Molec. 2010;40(6):487–95. https://doi.org/10.1016/j.ibmb.2010.04.008.
Xia XF, You MS, Rao XJ, Yu XQ. Insect C-type lectins in innate immunity. Dev Comp Immunol. 2018;83:70–9. https://doi.org/10.1016/j.dci.2017.11.020.
Zhou J, Zhao LL, Yu HY, Zhang W, Ahmad F, Hu SN, Zou Z, Sun JH. Comparative analysis of the Monochamus alternatus immune system. Insect Sci. 2018;25(4):581–603. https://doi.org/10.1111/1744-7917.12453.
Zaidman-Rémy A, Hervé M, Poidevin M, Pili-Floury S, Kim MS, Blanot D, Oh BH, Ueda R, Mengin-Lecreulx D, Lemaitre B. The Drosophila amidase PGRP-LB modulates the immune response to bacterial infection. Immunity. 2006;24(4):463–73. https://doi.org/10.1016/j.immuni.2006.02.012.
Kurata S. Peptidoglycan recognition proteins in Drosophila immunity. Dev Comp Immunol. 2014;42(1):36–41. https://doi.org/10.1016/j.dci.2013.06.006.
Schmidt RL, Trejo TR, Plummer TB, Platt JL, Tang AH. Infection-induced proteolysis of PGRP-LC controls the IMD activation and melanization cascades in Drosophila. FASEB J. 2008;22(3):918–29. https://doi.org/10.1096/fj.06-7907com.
Zhan MY, Yang PJ, Rao XJ. Molecular cloning and analysis of PGRP-L1 and IMD from silkworm Bombyx mori. Comp Biochem Phys B. 2018;215:19–30. https://doi.org/10.1016/j.cbpb.2017.10.002.
Chu Y, Liu Y, Shen DX, Hong F, Wang GR, An CJ. Serine proteases SP1 and SP13 mediate the melanization response of Asian corn borer, Ostrinia furnacalis, against entomopathogenic fungus Beauveria bassiana. J Invertebr Pathol. 2015;128:64–72. https://doi.org/10.1016/j.jip.2015.02.010.
Kambris Z, Brun S, Jang IH, Nam HJ, Romeo Y, Takahashi K. Drosophila immunity: a large-scale in vivo RNAi screen identifies five serine proteases required for toll activation. Curr Biol. 2006;16(8):808–13. https://doi.org/10.1016/j.cub.2006.03.020.
Roh KB, Kim CH, Lee H, Kwon HM, Park JW, Ryu JH, Kurokawa K, Ha NC, Lee WJ, Lemaitre B, Söderhäll K, Lee BL. Proteolytic cascade for the activation of the insect toll pathway induced by the fungal cell wall component. J Biol Chem. 2009;284(29):19474–81. https://doi.org/10.1074/jbc.M109.007419.
Lee SY, Kwon TH, Hyun JH, Choi JS, Kawabata S, Iwanaga S, Lee BL. In vitro activation of prophenoloxidase by two kinds of pro-phenol oxidase activating factors isolated from hemolymph of coleopteran, Holotrichia diomphalia larvae. Eur J Biochem. 1998;254(1):50–7. https://doi.org/10.1046/j.1432-1327.1998.2540050.x.
Liu CY, Song YY, Ren HN, Sun GG, Liu RD, Jiang P, Long SR, Zhang X, Wang ZQ, Cui J. Cloning and expression of a Trichinella spiralis putative glutathione s-transferase and its elicited protective immunity against challenge infections. Parasite Vector. 2017;10(1):448. https://doi.org/10.1186/s13071-017-2384-1.
Feng QL, Davey KG, Pang ASD, Primavera M, Ladd TR, Zheng SC, Sohi SS, Retnakaran A, Palli SR. Glutathione s-transferase from the spruce bud worm, Choristoneura fumiferana: identification, characterization, localization, cDNA cloning, and expression. Insect Biochem Molec. 1999;29(9):779–93. https://doi.org/10.1016/S0965-1748(99)00048-X.
Ketterman AJ, Saisawang C, Wongsantichon J. Insect glutathione transferases. Drug Metab Rev. 2011;43(2):253–65. https://doi.org/10.3109/03602532.2011.552911.
Altincicek B, Elashry A, Guz N, Grundle FA, Vilcinskas A, Dehne HW. Next generation sequencing based transcriptome analysis of septic-injury responsive genes in the beetle Tribolium castaneum. PLoS One. 2013;8:e52004.
Lemaitre B, Hoffmann J. The host defense of Drosophila melanogaster. Annu Rev Immunol. 2007;25(1):697–743. https://doi.org/10.1146/annurev.immunol.25.022106.141615.
Vogel H, Badapanda C, Knorr E, Vilcinskas A. Identifification of immunity-related genes in the burying beetle Nicrophorus vespilloides by suppression subtractive hybridization. Insect Mol Biol. 2011;20(6):787–800. https://doi.org/10.1111/j.1365-2583.2011.01109.x.
Malan AP, Konetze R, Moore SD. Isolation and identification of entomopathogenic nematodes from citrus orchards in South Africa and their biocontrol potential against false codling moth. J Invertebr Pathol. 2011;108(2):115–25. https://doi.org/10.1016/j.jip.2011.07.006.
Grabherr MG, Haas BJ, Yassour M, Levin J, Thompson DA, Amit I, et al. Full-length transcriptome assembly from RNA-seq data without a reference genome. Nat Biotechnol. 2011;29(7):644–52. https://doi.org/10.1038/nbt.1883.
Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B. Mapping and quantifying mammalian transcriptomes by RNA-seq. Nat Methods. 2008;5(7):621–8. https://doi.org/10.1038/nmeth.1226.
Anders S, Huber W. Differential expression analysis for sequence count data. Genome Biol. 2010;11(10):R106. https://doi.org/10.1186/gb-2010-11-10-r106.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(−Delta Delta C(T)) method. Methods. 2001;25(4):402–8. https://doi.org/10.1006/meth.2001.1262.
Abbott WS. A method of computing the effectiveness of an insecticide. J Econ Entomol. 1925;18(2):265–7. https://doi.org/10.1093/jee/18.2.265a.
We thank the Annoroad Gene Technology (Beijing) Co. Ltd. (Beijing, China) for the transcriptome sequencing and analysis.
This work was supported by the National Key R&D Program of China [2017YFE0130400, 2018YFD0201000], and National Natural Science Foundation of China .
Ethics approval and consent to participate
Animals used in this research were important agricultural pests, therefore approval from any ethics committee or institutional review board was not necessary. No ethics approval was required for the collection or experimentation of the animal used in this study.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1: Table S1.
The annotations of PGRP-SC1, PPAE-I and GST. Table S2:. The annotations of symbiotic bacteria in the transcriptomes of the Holotrichia parallela larvae postexposure to H. beicherriana LF for 72 h. Table S3. Summary of all species name, gene name, abbreviation and GenBank accession involved in phylogenetic analysis. Table S4. Oligonucleotides used for qRT-PCR and RNAi. The T7 promoter sequence is bolded in RNAi. Fig. S1. Rationality evaluation results of PGRP-SC1 (a), PPAE-I (b) and GST (c) 3D model. The red regions (A, B, L) represents the residues in most favoured regions; Bright yellow regions (a, b, l, p) represents the residues in additional allowed regions; Dark yellow regions (~a, ~b, ~l, ~p) represents the residues in generously allowed regions.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Li, E., Qin, J., Feng, H. et al. Immune-related genes of the larval Holotrichia parallela in response to entomopathogenic nematodes Heterorhabditis beicherriana LF. BMC Genomics 22, 192 (2021). https://doi.org/10.1186/s12864-021-07506-4
- Entomopathogenic nematode
- Holotrichia parallela
- Immune response