Transcriptome immune analysis of the invasive beetle Octodonta nipae (Maulik) (Coleoptera: Chrysomelidae) parasitized by Tetrastichus brontispae Ferrière (Hymenoptera: Eulophidae)

The Symbiotic Bacteria-Xenorhabdus nematophila All and Photorhabdus luminescens H06 Strongly Affected the Phenoloxidase Activation of Nipa Palm Hispid, Octodonta nipae (Coleoptera: Chrysomelidae) Larvae

Symbiotic bacteria form a mutualistic relationship with nematodes and are pathogenic to many insect pests. They kill insects using various strategies to evade or suppress their humoral and cellular immunity. Here we evaluate the toxic effects of these bacteria and their secondary metabolites on the survival and phenoloxidase (PO) activation of Octodonta nipae larvae using biochemical and molecular methods. The results show P. luminescens H06 and X. nematophila All treatments caused significant reductions in the number of O. nipae larvae in a dose-dependent manner. Secondly, the O. nipae immune system recognizes symbiotic bacteria at early and late stages of infection via the induction of C-type lectin. Live symbiotic bacteria significantly inhibit PO activity in O. nipae whereas heat-treated bacteria strongly increase PO activity. Additionally, expression levels of four O. nipae proPhenoloxidase genes following treatment with P. luminescens H06 and X. nematophila All were compared. We found that the expression levels of all proPhenoloxidase genes were significantly down-regulated at all-time points. Similarly, treatments of O. nipae larvae with metabolites benzylideneacetone and oxindole significantly down-regulated the expression of the PPO gene and inhibited PO activity. However, the addition of arachidonic acid to metabolite-treated larvae restored the expression level of the PPO gene and increased PO activity. Our results provide new insight into the roles of symbiotic bacteria in countering the insect phenoloxidase activation system.

The Entomopathogenic Nematodes H. bacteriophora and S. carpocapsae Inhibit the Activation of proPO System of the Nipa Palm Hispid Octodonta nipae (Coleoptera: Chrysomelidae)

Entomopathogenic nematodes are biocontrol agents of invasive insect pests in soil and cryptic habitats. Nipa palm hispid, Octodonta nipae, is a pest of palm trees in Sothern China. To address its increasing damage, environmentally friendly control methods are required. This study aimed to test efficacy of Heterorhabditis bacteriophora and Steinernema carpocapsae on O. nipae and investigated the influence of secondary metabolites, nematodes, and their isolated cuticles on the activation of O. nipae’s prophenoloxidase system using qPCR analysis. Our data revealed that O. nipae were less susceptible to H. bacteriophora than S. carpocapsae and penetrations of infective juveniles were higher with S. carpocapsae treatment than H. bacteriophora. Moreover, expression levels of the serine protease P56, prophenoloxidase activation factor 1, PPO and serine protease inhibitor 28 upon S. carpocapsae and H. bacteriophora infections were generally downregulated at all times. However, upon heating, the cuticles lost their inhibitory effects and resulted in upregulation of the PPO gene. Similarly, the addition of arachidonic acid reversed the process and resulted in the upregulation of the PPO gene compared to the control. Further work is needed to identify toxic substances secreted by these EPNs to evade O. nipae’s immune system.

Variation in Parasitoid Virulence of Tetrastichus brontispae during the Targeting of Two Host Beetles

In host-parasitoid interactions, antagonistic relationship drives parasitoids to vary in virulence in facing different hosts, which makes these systems excellent models for stress-induced evolutionary studies. Venom compositions varied between two strains of Tetrastichus brontispae, Tb-Bl and Tb-On. Tb-Bl targets Brontispa longissima pupae as hosts, and Tb-On is a sub-population of Tb-Bl, which has been experimentally adapted to a new host, Octodonta nipae. Aiming to examine variation in parasitoid virulence of the two strains toward two hosts, we used reciprocal injection experiments to compare effect of venom/ovarian fluids from the two strains on cytotoxicity, inhibition of immunity and fat body lysis of the two hosts. We found that Tb-Onvenom was more virulent towards plasmatocyte spreading, granulocyte function and phenoloxidase activity than Tb-Blvenom. Tb-Blovary was able to suppress encapsulation and phagocytosis in both hosts; however, Tb-Onovary inhibition targeted only B. longissima. Our data suggest that the venom undergoes rapid evolution when facing different hosts, and that the wasp has good evolutionary plasticity.

Isolation, Identification, and Bioinformatic Analysis of Antibacterial Proteins and Peptides from Immunized Hemolymph of Red Palm Weevil Rhynchophorus ferrugineus

Red palm weevil (Rhynchophorus ferrugineus Olivier, 1791, Coleoptera: Curculionidae) is a destructive pest of palms, rapidly extending its native geographical range and causing large economic losses worldwide. The present work describes isolation, identification, and bioinformatic analysis of antibacterial proteins and peptides from the immunized hemolymph of this beetle. In total, 17 different bactericidal or bacteriostatic compounds were isolated via a series of high-pressure liquid chromatography steps, and their partial amino acid sequences were determined by N-terminal sequencing or by mass spectrometry. The bioinformatic analysis of the results facilitated identification and description of corresponding nucleotide coding sequences for each peptide and protein, based on the recently published R. ferrugineus transcriptome database. The identified compounds are represented by several well-known bactericidal factors: two peptides similar to defensins, one cecropin-A1-like peptide, and one attacin-B-like protein. Interestingly, we have also identified some unexpected compounds comprising five isoforms of pheromone-binding proteins as well as seven isoforms of odorant-binding proteins. The particular role of these factors in insect response to bacterial infection needs further investigation.

Cloning and Immunosuppressive Properties of an Acyl-Activating Enzyme from the Venom Apparatus of Tetrastichus brontispae (Hymenoptera: Eulophidae)

Venom injected into the host plays vital roles in facilitating successful parasitization and development for parasitoid wasps, especially those devoid of polydnavirus, and the abundant venom proteins appear to be most likely involved in parasitization success. Previously, we found the four most abundant venom proteins, including 4-coumarate:CoA ligase-like 4 (4CL4-like), in the Tetrastichus brontispae (Hymenoptera: Eulophidae) venom apparatus. In this study, we cloned, expressed T. brontispae 4CL4-like (Tb4CL4-like) in Escherichia coli, and investigated its immunosuppressive properties. The deduced amino acid sequence for Tb4CL4-like shares high identity at conserved amino acids associated with the acyl-activating enzyme (AAE) consensus motif but shows only <40% identity with the members in the AAE superfamily. mRNA abundance analysis indicated that Tb4CL4-like was transcribed mainly in the venom apparatus. Recombinant Tb4CL4-like inhibited Octodonta nipae (Coleoptera: Chrysomelidae) pupal cellular encapsulation and spreading by targeting the hemocyte cytoskeleton and reduced the hemocyte-mediated phagocytosis of E. coli in vivo. Moreover, Tb4CL4-like exhibited greater affinity to palmitic acid and linolenic acid based on the molecular docking assay and is hypothesized to be involved in fatty acid metabolism. In conclusion, our results suggest that Tb4CL4-like may be an immunity-related AAE protein that is involved in the regulation of host immunity through fatty acid metabolism-derived signaling pathways.

Exploring the Role of Relish on Antimicrobial Peptide Expressions (AMPs) Upon Nematode-Bacteria Complex Challenge in the Nipa Palm Hispid Beetle, Octodonta nipae Maulik (Coleoptera: Chrysomelidae)

The humoral immune responses of the nipa palm hispid beetle Octodonta nipae involves the inducible expression of the genes coding for antimicrobial peptides (AMPs) which are mediated by immune deficiency signaling pathways. In insects, the nuclear factor-κB (NF-κB) transcription factor, Relish, has been shown to regulate AMP gene expressions upon microbial infections. Here, we dissect the expression patterns of some AMPs in O. nipae during infections by entomopathogenic nematodes (EPNs) and their symbionts, before and after Relish knock down. Our results indicate that, prior to gene silencing, the AMPs attacin C1, attacin C2, and defensin 2B were especially expressed to great extents in the insects challenged with the nematodes Steinernema carpocapsae and Heterorhabditis bacteriophora as well as with their respective symbionts Xenorhabdus nematophila and Photorhabdus luminescens. The study also established the partial sequence of OnRelish/NF-κB p110 subunit in O. nipae, with an open reading frame coding for a protein with 102 amino acid residues. A typical Death domain-containing protein was detected (as seen in Drosophila) at the C-terminus of the protein. Phylogenetic analysis revealed that in O. nipae, Relish is clustered with registered Relish/NF-κB p110 proteins from other species of insect especially Leptinotarsa decemlineata from the same order Coleoptera. Injection of OnRelish dsRNA remarkably brought down the expression of OnRelish and also reduced the magnitude of transcription of attacin C1 and defensin 2B upon S. carpocapsae and H. bacteriophora and their symbionts infections. Altogether, our data unveil the expression pattern of OnRelish as well as that of some AMP genes it influences during immune responses of O. nipae against EPNs and their symbionts.

Beetles as Model Organisms in Physiological, Biomedical and Environmental Studies - A Review

Model organisms are often used in biological, medical and environmental research. Among insects, Drosophila melanogaster, Galleria mellonella, Apis mellifera, Bombyx mori, Periplaneta americana, and Locusta migratoria are often used. However, new model organisms still appear. In recent years, an increasing number of insect species has been suggested as model organisms in life sciences research due to their worldwide distribution and environmental significance, the possibility of extrapolating research studies to vertebrates and the relatively low cost of rearing. Beetles are the largest insect order, with their representative - Tribolium castaneum - being the first species with a completely sequenced genome, and seem to be emerging as new potential candidates for model organisms in various studies. Apart from T. castaneum, additional species representing various Coleoptera families, such as Nicrophorus vespilloides, Leptinotarsa decemlineata, Coccinella septempunctata, Poecilus cupreus, Tenebrio molitor and many others, have been used. They are increasingly often included in two major research aspects: biomedical and environmental studies. Biomedical studies focus mainly on unraveling mechanisms of basic life processes, such as feeding, neurotransmission or activity of the immune system, as well as on elucidating the mechanism of different diseases (neurodegenerative, cardiovascular, metabolic, or immunological) using beetles as models. Furthermore, pharmacological bioassays for testing novel biologically active substances in beetles have also been developed. It should be emphasized that beetles are a source of compounds with potential antimicrobial and anticancer activity. Environmental-based studies focus mainly on the development and testing of new potential pesticides of both chemical and natural origin. Additionally, beetles are used as food or for their valuable supplements. Different beetle families are also used as bioindicators. Another important research area using beetles as models is behavioral ecology studies, for instance, parental care. In this paper, we review the current knowledge regarding beetles as model organisms and their practical application in various fields of life science.

The transcriptome analysis of Protaetia brevitarsis Lewis larvae

Larvae of the pest Protaetia brevitarsis are used to treat infections in traditional Chinese medicine. However, genomic information about this non-model species is currently lacking. To better understand the fundamental biology of this non-model species, its transcriptome was obtained using next generation sequencing and then analyzed. A total of 7.62 Gb of clean reads were obtained, which were assembled into 169,087 transcripts corresponding to 142,000 annotated unigenes. These unigenes were functionally classified according to Gene Ontology (GO), euKaryotic Ortholog Groups of proteins (KOG), and Kyoto Encyclopedia of Genes and Genomes (KEGG) annotations. A total of 41,921 unigenes were assigned to 56 GO terms, 21,454 unigenes were divided among 26 KOG categories, and 16,368 unigenes were assigned to 32 KEGG pathways. In addition, 19,144 simple sequence repeats (SSRs) were identified. Furthermore, several kinds of natural antimicrobial peptides and proteins, 4 histones with potential antimicrobial activity, and 41 potential antimicrobial peptide sequences were identified. These data are the first reported whole transcriptome sequence of P. brevitarsis larvae, which represents a valuable genomic resource for studying this species, thus promoting the utilization of its medical potential.

The First Complete Genome Sequence of a Novel Tetrastichus brontispae RNA Virus-1 (TbRV-1)

The complete sequence of a novel RNA virus isolated from Tetrastichus brontispae (TbRV-1) was determined to be 12,239 nucleotides in length with five non-overlapping, linearly arranged coding sequences (CDS), potentially encoding nucleoproteins, hypothetical proteins, matrix proteins, glycoproteins, and RNA-dependent RNA polymerases. Sequence analysis indicated that the RNA-dependent RNA polymerase of TbRV-1 shares a 65% nucleotide and 67% amino acid sequence identity with Hubei dimarhabdovirus 2, suggesting that TbRV-1 is a member of the dimarhabdovirus supergroup. This corresponded to the result of the phylogenetic analysis. The affiliation of TbRV-1 with members of the family Rhabdoviridae was further validated by similar transcription termination motifs (GGAACUUUUUUU) to the Drosophila sigmavirus. The prevalence of TbRV-1 in all tissues suggested that the virus was constitutive of, and not specific to, any wasp tissue. To our knowledge, this is the first report on the complete genome sequence of a dimarhabdovirus in parasitoids.

The Endosymbiotic Wolbachia and Host COI Gene Enables to Distinguish Between Two Invasive Palm Pests; Coconut Leaf Beetle, Brontispa longissima and Hispid Leaf Beetle, Octodonta nipae

To elucidate taxonomic eminence of identical pest species is essential for many ecological and conservation studies. Without proficient skills, accurate molecular identification and characterization are laborious and time-consuming. The coconut leaf beetle, Brontispa longissima (Gestro) (Coleoptera: Chrysomelidae), is biologically and morphologically identical to hispid leaf beetle, Octodonta nipae (Maulik) (Coleoptera: Chrysomelidae), and is known as the most harming nuisances of palm cultivation worldwide. The present examination was to establish Wolbachia genotyping analysis along with host cytochrome oxidase subunit I (COI) gene for accurate identification between these individuals of the same family (Chrysomelidae). Here, we have cloned and sequenced a gene coding Wolbachia surface protein (wsp) and COI gene regions amplified from both species by polymerase chain reaction. The nucleotide sequences were directly determined (≈600 bp for wsp and ≈804 bp for COI) and aligned using the multiple alignment algorithms in the ESPript3 package and the MEGA5 program. Comparative sequence analysis indicated that the representative of wsp and COI sequences from these two beetles were highly variable. To ensure this bacterial variation, multilocus sequence typing (MLST) of bacterial genes was conducted, and the results vindicated the same trend of variations. Furthermore, the phylogenetic analysis also indicates that B. longissima and O. nipae being the two different species harbors two distinct Wolbachia Hertig and Burt (Rickettsiales: Anaplamataceae) supergroups B and A, respectively. The present outcomes quickly discriminate between these two species. Considering its simplicity and cost-effectiveness, it can be used as a diagnostic tool for discriminating such invasive species particularly B. longissima and O. nipae which has overlapping morphologic characters.

Development of Microsatellite Markers for the Nipa Palm Hispid Beetle, Octodonta nipae (Maulik)

The nipa palm hispid beetle, Octodonta nipae (Maulik) (Coleoptera: Chrysomelidae), is an important invasive pest on palm plants in southern China. Based on existing transcriptome data, polymorphism simple sequence repeat (SSR) loci were identified. In total, 1274 SSR loci were identified from 49919 unigenes. The majority of them contained mononucleotide, dinucleotide, and trinucleotide motifs (43.56%, 26.14%, and 28.18%), in which A/T (41.21%) and AT/TA (15.86%) were the most abundant motifs. 104 pairs of the SSR primers produced amplification bands of expected sizes in O. nipae, 80 pairs of SSR primers were tested randomly for polymorphism, 9 loci of them were validated to be polymorphic markers, and the number of alleles ranged from 2 to 3, with an average of 2.56 per locus. The population of Zhangzhou and Fuzhou was analyzed by the 9 loci (On1-On9). These SSR transcriptome data can provide invaluable resource for SSR development, population genetics research, invasion and expansion mechanism, paternity testing, and other research on O. nipae and its related species.

Descriptions of immature stages of Octodonta nipae (Maulik) (Coleoptera, Chrysomelidae, Cassidinae, Cryptonychini)

Octodonta nipae (Maulik, 1921), a hispid that damages several species of palm trees, was introduced accidently into China in 2001. The egg, larva, prepupa and pupa of O. nipae are illustrated and described in detail and compared with another invasive species, Brontispa longissima (Gestro, 1885); the scanning electron micrographs of the head capsule, antenna, maxilla, labium and lateral scoli are provided, as well as photos of body of all larval instars and pupa. It is the second description of immature stages in the genus Octodonta Chapuis.

A novel bacterial symbiont association in the hispid beetle, Octodonta nipae (Coleoptera: Chrysomelidae), their dynamics and phylogeny

The hispid leaf beetle, Octodonta nipae (Maulik), (Coleoptera: Chrysomelidae), is a devastating pest of palm cultivation worldwide. Endosymbiotic bacteria in the genus Wolbachia are arguably one of the most abundant bacterial group associated with arthropods. Owing to its critical effects on host reproduction, Wolbachia has garnered much attention as a prospective future tool for insect pest management. However, their association, infection dynamics, and functionality remain unknown in this insect pest. Here, we diagnosis for the first time, the infection prevalence, and occurrence of Wolbachia in O. nipae. Experimental evidence by the exploration of wsp gene vindicate that O. nipae is naturally infected with bacterial symbiont of genus Wolbachia, showing a complete maternal inheritance with shared a common Wolbachia strain (wNip). Moreover, MLST (gatB, fbpA, coxA, ftsZ, and hcpA) analysis enabled the detections of new sequence type (ST-484), suggesting a particular genotypic association of O. nipae and Wolbachia. Subsequently, quantitative real-time PCR (qPCR) assay demonstrated variable infection density across different life stages (eggs, larvae, pupae and adult male and female), body parts (head, thorax, abdomen), and tissues (ovaries, testes, and guts). Infection density was higher in egg and female adult stage, as well as abdomen and reproductive tissues as compared to other samples. Interestingly, Wolbachia harbored dominantly in a female than the male adult, while, no significant differences were observed between male and female body parts and tissues. Phylogeny of Wolbachia infection associated with O. nipae rectified from all tested life stages were unique and fall within the same monophyletic supergroup-A of Wolbachia clades. The infection density of symbiont is among the valuable tool to understand their biological influence on hosts, and this latest discovery would facilitate the future investigations to understand the host-symbiont complications and its prospective role as a microbiological agent to reduce pest populations.

Identification of three prophenoloxidase-activating factors (PPAFs) from an invasive beetle Octodonta nipae Maulik (Coleoptera: Chrysomelidae) and their roles in the prophenoloxidase activation

A typical characteristic of the insect innate immune system is the activation of the serine protease cascade in the hemolymph. As being the terminal component of the extracellular serine protease cascade in the prophenoloxidase (proPO) activating system, proPO-activating factors (PPAFs) activated by the upstream cascade may generate active phenoloxidase, which then induces downstream melanization. In the present study, we reported three PPAFs from the nipa palm hispid beetle Octodonta nipae (Maulik) (designated as OnPPAF1, OnPPAF2, OnPPAF3). All three OnPPAFs contained a single clip domain at the amino-terminus followed by a trypsin-like serine protease domain at the carboxyl-terminus, except the Ser in the active sites of OnPPAF2 and OnPPAF3 was substituted with Gly. Transcript expression analysis revealed that all OnPPAFs were highly expressed in hemolymph, whereas OnPPAF2 showed an extremely low mRNA abundance compared with that of OnPPAF1 and OnPPAF3, and that the abundance of all three OnPPAFs was dramatically increased upon bacterial challenge. Knockdown of OnPPAF1 or OnPPAF3 resulted in a reduction of hemolymph phenoloxidase activity and an inhibition of hemolymph melanization, whereas the knockdown of OnPPAF2 did not affect the proPO cascade. Our work thus implies that the three OnPPAFs may have different functions and regulation during immune responses in O. nipae.

Parasitism Performance of Tetrastichus brontispae Ferriere over the Coconut Hispine Beetle, Brontispa longissima (Gestro)

In this study, the effect of host density, host, and parasitoid ages in choice and no-choice tests on the parasitism performance of Tetrastichus brontispae Ferriere, one of the major parasitoid of Brontispa longissima (Gestro), was investigated in the laboratory. The results revealed that an increased host density resulted in no increased parasitism of B. longissima by T. brontispae; the optimal host density was three host pupae per parasitoid when considering the costs for mass rearing. Moreover, parasitoid age was quite crucial for effective parasitism and affected the emergence rate. Although 2-h to 4-day-old parasitoids successfully parasitized the host pupae, younger parasitoids (within 2-day-old) presented higher parasitism capacity than older parasitoids. More importantly, both choice and no-choice tests confirmed that all host stages tested from 2-h to 4-day-old were suitable for T. brontispae parasitization, although 2-h to 2-day-old hosts were preferred. We also demonstrated that sex ratio, emergence rate, and egg to adult developmental time were not influenced by host density, parasitoid, and host age in both choice and no-choice tests. Our data will allow for more accurate prediction and interpretation on the parasitization by T. brontispae, supporting mass-production initiatives and mass release in programs of B. longissima.

Deep Sequencing-Based Transcriptome Analysis Reveals the Regulatory Mechanism of Bemisia tabaci (Hemiptera: Aleyrodidae) Nymph Parasitized by Encarsia sophia (Hymenoptera: Aphelinidae)

The whitefly Bemisia tabaci is a genetically diverse complex with multiple cryptic species, and some are the most destructive invasive pests of many ornamentals and crops worldwide. Encarsia sophia is an autoparasitoid wasp that demonstrated high efficiency as bio-control agent of whiteflies. However, the immune mechanism of B. tabaci parasitization by E. sophia is unknown. In order to investigate immune response of B. tabaci to E. Sophia parasitization, the transcriptome of E. sophia parasitized B. tabaci nymph was sequenced by Illumina sequencing. De novo assembly generated 393,063 unigenes with average length of 616 bp, in which 46,406 unigenes (15.8% of all unigenes) were successfully mapped. Parasitization by E. sophia had significant effects on the transcriptome profile of B. tabaci nymph. A total of 1482 genes were significantly differentially expressed, of which 852 genes were up-regulated and 630 genes were down-regulated. These genes were mainly involved in immune response, development, metabolism and host signaling pathways. At least 52 genes were found to be involved in the host immune response, 33 genes were involved in the development process, and 29 genes were involved in host metabolism. Taken together, the assembled and annotated transcriptome sequences provided a valuable genomic resource for further understanding the molecular mechanism of immune response of B. tabaci parasitization by E. sophia.

THE ENDOPARASITOID Pteromalus puparum INFLUENCES HOST GENE EXPRESSION WITHIN FIRST HOUR OF PARASITIZATION

The small cabbage butterfly, Pieris rapae, is an important pest of cruciferous corps, and Pteromalus puparum is a predominant pupal endoparasitoid wasp of this butterfly. For successful development of parasitoid offspring, female parasitoids usually introduce one or several kinds of maternal factors into the hemocoels during oviposition to suppress host immunity. To investigate the early changes in host immune-related genes following parasitization, we analyzed transcriptomes of parasitized and unparasitized, control, host pupae. Approximately 17.7 and 19.3 million paired-end reads were generated from nonparasitized and parasitized host pupae, and assembled de novo into 45,639 transcripts and 27,659 nonredundant unigenes. The average unigene length was 790 bp. A total 18,377 of 27,659 unigenes were annotated and we identified 557 differentially expressed unigenes in host pupae at 1 h after parasitization, of which 21 were immune-related. Parasitization led to downregulation of most pattern recognition receptors and upregulation of all serine protease inhibitors. The transcirptomic profile of P. rapae is considerably affected by parasitization. This study provides valuable sources for future investigations of the molecular interaction between P. puparum and its host P. rapae.

A Rapid Diagnostic Technique to Discriminate between Two Pests of Palms, Brontispa longissima and Octodonta nipae (Coleoptera: Chrysomelidae), for Quarantine Applications

Octodonta nipae (Maulik) is morphologically and biologically similar to Brontispa longissima (Gestro), one of the most damaging pests of coconut. The two species share several palm hosts and produce similar symptoms. They are easily confused when they invade a new area, and without professional taxonomic expertise accurate identification is arduous and time-consuming. Thus, a method of rapidly distinguishing these two invasive insects rapidly is critical for quarantine. Based on the first internal transcribed spacer and the mitochondrial cytochrome oxidase I gene, two pairs of special primers and corresponding polymerase chain reaction processes have been developed to enhance a single objective band only from the O. nipae DNA template. This will quickly discriminate between these two species. The present results provide a rapid method of distinguishing O. nipae from B. longissima to help avoid misidentification, and furthermore, to facilitate rapid and appropriate quarantine decisions and effective treatments to the pest.