- Research article
- Open Access
High throughput profiling of the cotton bollworm Helicoverpa armigera immunotranscriptome during the fungal and bacterial infections
- Guang-Hua Xiong†1, 4,
- Long-Sheng Xing†1, 4,
- Zhe Lin1,
- Tusar T Saha2,
- Chengshu Wang3,
- Haobo Jiang5 and
- Zhen Zou1Email author
© Xiong et al. 2015
- Received: 19 December 2014
- Accepted: 31 March 2015
- Published: 18 April 2015
Innate immunity is essential in defending against invading pathogens in invertebrates. The cotton bollworm, Helicoverpa armigera (Hübner) is one of the most destructive lepidopteran pests, which causes enormous economic losses in agricultural production worldwide. The components of the immune system are largely unknown in this insect. The application of entomopathogens is considered as an alternative to the chemical insecticides for its control. However, few studies have focused on the molecular mechanisms of host-pathogen interactions between pest insects and their pathogens. Here, we investigated the immunotranscriptome of H. armigera larvae and examined gene expression changes after pathogen infections. This study provided insights into the potential immunity-related genes and pathways in H. armigera larvae.
Here, we adopted a high throughput RNA-seq approach to determine the immunotranscriptome of H. armigera larvae injected with buffer, fungal pathogen Beauveria bassiana, or Gram-negative bacterium Enterobacter cloacae. Based on sequence similarity to those homologs known to participate in immune responses in other insects, we identified immunity-related genes encoding pattern recognition receptors, signal modulators, immune effectors, and nearly all members of the Toll, IMD and JAK/STAT pathways. The RNA-seq data indicated that some immunity-related genes were activated in fungus- and bacterium-challenged fat body while others were suppressed in B. bassiana challenged hemocytes, including the putative IMD and JAK-STAT pathway members. Bacterial infection elevated the expression of recognition and modulator genes in the fat body and signal pathway genes in hemocytes. Although fat body and hemocytes both are important organs involved in the immune response, our transcriptome analysis revealed that more immunity-related genes were induced in the fat body than that hemocytes. Furthermore, quantitative real-time PCR analysis confirmed that, consistent with the RNA-seq data, the transcript abundances of putative PGRP-SA1, Serpin1, Toll-14, and Spz2 genes were elevated in fat body upon B. bassiana infection, while the mRNA levels of defensin, moricin1, and gloverin1 were up-regulated in hemocytes.
In this study, a global survey of the host defense against fungal and bacterial infection was performed on the non-model lepidopteran pest species. The comprehensive sequence resource and expression profiles of the immunity-related genes in H. armigera are acquired. This study provided valuable information for future functional investigations as well as development of specific and effective agents to control this pest.
- Helicoverpa armigera
- Fungal infection
- Bacterial challenge
- Innate immunity
- Fat body
The cotton bollworm Helicoverpa armigera, is a destructive and highly polyphagous insect pest in Asia, Europe, Africa, and Australia. It causes serious damage to cotton, sorghum, corn and several other crops. Further aggravating the situation, this pest is spreading geographically by invading newer territory. Recently, their presence outside the Asian continent has been confirmed with reports of establishment of this species in Brazil . Primary methods of controlling this devastating agricultural pest include traditional spraying of pesticides and cultivation of transgenic pest resistant crop varieties developed in the recent past. However, apart from causing environmental pollution these methods also lead to resistance development in pests, thus necessitating immediate development of novel biological control methods. Entomopathogenic fungi like Beauveria bassiana has been widely considered as a potent eco-friendly bio-control agent . The potential of B. bassiana as an effective suppressant of H. armigera and other lepidopteran pests is under consideration . Before B. bassiana kills the pest, it needs to invade and overcome the host immune responses that have not been well understood in H. armigera.
The invasion of pathogens and parasites into an insect is defended by innate immune system, a strong universal defense mechanism shared by both vertebrates and invertebrates . This system consists of initial physical barriers to pathogen penetration followed by humoral and cellular response against subsequent microbial infection. While the cellular response involves phagocytosis, encapsulation, and nodule formation by the circulating hemocytes, the humoral response refers to the process of melanization and the production of immune effector molecules, which are mainly produced in the fat body. Anti-bacterial immunity depends on two principal signaling pathways, Toll and IMD, which are conserved across various insect species indicating their central importance throughout arthropod evolution [5,6]. These pathways are activated by the binding of pattern recognition receptors (PRRs) such as peptidoglycan recognition protein (PGRP), β-1,3-glucan recognition protein (βGRP), C-type lectin (CTL) to pathogen associated molecular pattern (PAMP) on the surface of invading microorganisms . Specifically, Gram-negative bacteria trigger the IMD pathway leading to an acute response, while Gram-positive bacteria and fungi activate the Toll pathway eliciting a more sustained and complex reaction . Both of these signaling pathways result in the induction of various NF-кB transcription factors, which then drive a comprehensive synthesis of immune effectors including antimicrobial peptides (AMPs) . Like immune signaling pathways, melanization which results from the activation of prophenoloxidase (PPO), is also considered as an universal immune response against bacteria in insects . It involves melanin synthesis, sequestering and killing the invading microorganisms or parasites .
After invading the host, B. bassiana triggers the cellular and humoral immune mechanism of the host resulting in melanized nodules and activation of Toll pathway [9,10]. The Toll signaling pathway is mediated by an extracellular protease cascade composed of a series of clip-domain serine protease (cSP), which helps to amplify the initial recognition signal . Several immune molecules like GNBP3, Persephone, and Drosomycin are involved in the fungus activated Toll pathway . Besides the Toll pathway, melanization is also implicated to play a role in antifungal immunity. B. bassiana infection of mosquito Aedes aegypti enhances the cleavage of hemolymph PPO and its conversion to phenoloxidase (PO), the rate limiting enzyme in melanogenesis . Similarly, induced melanization in another mosquito Anopheles gambiae retarded the growth and dissemination of B. bassiana . However, the thorough understanding of interaction between melanization and fungal challenge is largely unclear.
B. bassiana parasitizes a wide range of hosts including H. armigera and other lepidopteran pests. Most studies on insect immunity were performed in model organisms, such as Drosophila melanogaster , Manduca sexta , and Tenebrio molitor . Systematic analysis of immunity-related genes has been conducted in several other insects including A. gambiae, A. aegypti, Bombyx mori, Apis mellifera and Tribolium castaneum, whose genome sequences are available. It is only recently that H. armigera immunity started to attract some researchers’ attentions . Fourteen genes involved in the H. armigera interaction with Bacillus thuringiensis and Autographa californica multiple nucleocapsid nucleopolyhedrovirus were obtained . The expression and characterization of CTLs and PGRPs in response to pathogen were also reported [15-17]. However, H. armigera genome is still unavailable and its immune response to B. bassiana remains largely unexplored, restricting its further development and adoption as a biological agent in the control of lepidopteran pests.
More interestingly, systematic survey on H. armigera is of great importance for understanding the origin and evolution of immune systems, which is a challenging task since genetic backgrounds of non-model insects are unclear, let alone expression profiling. In this work, RNA-seq was employed to perform immunotranscriptome analysis of H. armigera larval hemocytes and fat body in response to the challenges of entomopathogenic fungi B. bassiana and Gram-negative bacteria Enterobactor cloacae. In an attempt to gain insights into the H. armigera immune response against pathogens, we characterized 233 such genes and divided them into pattern recognition, signal transduction, execution, and cellular responses. Hierarchical clustering analyses of differentially expressed gene (DEG) cohorts indicated that the fungi and bacteria elicited distinct expression profiles of genes in hemocytes and fat body. These results provide an overview of the tissue-specific expression profiles in response to infections and a platform for further exploring the molecular basis of host antimicrobial response.
Immunization, sequencing assembly, and gene identification
In D. melanogaster and mosquitoes, immunity-related genes regulated by IMD pathway showed an acute expression pattern involved in the anti-bacterial responses, whereas Toll pathway displayed sustained activation of immunity-related genes against fungal infection [5,18]. To gain detailed information about the H. armigera transcriptome, the hemocytes and fat body samples were collected from E. cloacae infected larvae at 6 h and B. bassiana infected larvae at 48 h, respectively. Six cDNA libraries were constructed from the RNA samples extracted from control fat body (FB_Mock), fungus-induced fat body (FB_Bb), bacterium-induced fat body (FB_Ec), control hemocytes (HC_Mock), fungus-induced hemocytes (HC_Bb), and bacterium-induced hemocytes (HC_Ec) (in Additional file 1: Figure S2). These libraries were sequenced on an Illumina Hiseq 2000 system. After removal of adaptor sequences and low-quality reads (Q < 20), these libraries yielded 72.6, 79.2, 65.6, 80.1, 46.7, and 80.4 million high-quality reads with 7.54, 7.70, 5.90, 7.53, 8.25, and 7.45 Gb data, respectively (in Additional file 2: Table S1). All the reads were de novo assembled into 57,588, 61,452, 55,394, 60,873, 42,073, and 58,262 contigs with average lengths of 1,102, 1,139, 1,034, 1,128, 995, and 1,188 nt, respectively [19,20]. Pooled clean reads from the six libraries and three others (E. cloacae, B. bassiana, and PBS injected 2nd instar larvae), totally yielded 150,606 contigs (average size: 1,124 nt) (in Additional file 2: Table S1).
To predict functions, all the obtained unigenes were annotated according to BLASTX searches against the non-redundant sequence database resulting 60,458 (43.2%) of them displaying homology to the known proteins (E-value < 1e-5) (in Additional file 1: Figure S3) . For species distribution, nearly 17,000 (11.2%) annotated unigenes were homologous to B. mori, followed by Danaus plexippus (6.6%). Only fewer transcripts were matched to the ones from H. armigera (0.6%), probably due to the unavailability of its genome sequences (in Additional file 1: Figure S3). Additionally, 713 transcripts with high similarity to non-coding RNAs by homology search were obtained (in Additional file 3: Table S2) . After removing redundant sequences and those shorter than 150 nt, we retained the unigenes with significant BLAST hits and generated a non-redundant dataset containing 37,694 sequences (average length: 1,711 nt) (in Additional file 2: Table S1). Length distribution of the H. armigera transcripts indicates that this final dataset has a high proportion of long transcripts (>2000 nt) and increases the chances of full length gene predictions. Based on the unigene information (e.g. gene ID, length, putative function) (in Additional file 4: Table S3), the immunity-related genes were identified and filtered out for further analysis. The detailed workflow from immune challenge to immunity-related gene identification is provided in a flowchart (in Additional file 1: Figure S2).
Identification and functional classification of differentially expressed genes in response to B. bassiana and E. cloacae infections
To gain insights into the tissue-specific transcriptional changes in H. armigera larvae infected by B. bassiana and E. cloacae, we performed pairwise comparisons between libraries to identify the differentially expressed transcripts (DETs) (in Additional file 1: Figure S4A) . Relative to the control, transcripts with greater than 2-fold change and p value less than 0.001 were considered as differentially expressed. The complete list of DETs with FPKM (fragments per kilobase of transcript per million) values (in Additional file 5: Table S4) revealed that more unigenes exhibited remarkable changes in mRNA levels in fat body (530 up- and 415 down-regulated) than in hemocytes (235 up- and 320 down-regulated) 6 h after E. cloacae injection (in Additional file 1: Figure S4B) indicating a tissue-specific affect. More pronounced changes occurred to the larvae 48 h after B. bassiana injection. Compared to the control group, there were 2,014 (1,469 up- and 545 down-regulated) transcripts that were significantly changed in the fat body elicited by B. bassiana, while more transcripts (496 up- and 2081 down-regulated) suppressed in hemocytes were observed.
In order to compare gene expression levels in hemocytes and fat body in response to the pathogen invasion, we performed hierarchical clustering of DETs and identified five discrete clusters showing expression trends relevant to the infection. The 487 and 2,015 genes in clusters 1 and 2 had the highest and the lowest RNA levels in hemocytes from B. bassiana infected larvae, respectively, among all the experimental conditions (Figure 1C and in Additional file 5: Table S4). Clusters 3 (1,280 genes) and 4 (370 genes) showed the most significant mRNA level increase in fat body after the fungal and bacterial infection, respectively. Furthermore, cluster 5 represented gene cohorts significantly suppressed in hemocytes and fat body after the infections. Consistent with the information revealed by hierarchical clustering, a Venn diagram analysis indicated that 1,251 transcripts were exclusively regulated in the fungus challenged fat body (1,109 up- and 142 down-regulated), while 1,857 other transcripts were specifically regulated in the fungus challenged hemocytes (290 up- and 1,567 down-regulated) (Figure 1D). In comparison, fewer gene transcripts (222 + 67 + 101 + 15) were exclusively regulated in the two tissues after the bacterial infection. Considering that the gene repertoire of fat body had small overlap with that of hemocytes (Figure 1B), we speculate that tissue type instead of the pathogen is more important in regulating the gene transcription under the same experimental condition. Besides, the data showed consistent up-regulation of 12 genes in both fat body and hemocytes after the fungal and bacterial infection and most of them were immunity-related.
Taken together, our results clearly indicated that global gene expression profiles varied among different tissues and treatments (Figures 1 and 2). While tissue is the major determinant, other factors also have minor but significant impacts on the transcriptome. The distinct GO and KEGG distribution of up- and down-regulated gene cohorts showed similar trends. Detailed expression profiling of the immunity-related genes was still required for further investigation of tissue-specific immune responses against B. bassiana and E. cloacae infection.
For in-depth analysis of the interactions between H. armigera and pathogens, we annotated the insect defense genes and examined their phylogenetic relationships. By sequencing the tissues that were enriched in the immune-related transcripts at sufficient depth, followed by the BLAST search , we were able to detect 233 immunity-related genes (Figure 2D, in Additional file 6: Table S5). This number was significantly higher when compared with other insects from holometabolous orders, such as A. mellifera, but similar to that in B. mori. This result indicated that H. armigera, well adapted to various plant hosts and environment, has a sizable repertoire of genes related to the cellular and humoral immune responses against the wounding and infection.
Comparative analysis of genes for immune signal recognition
When pathogens infect insects, multiple PRRs may recognize conserved determinants (e.g. lipopolysaccharide, peptidoglycan, lipoteichoic acid, β-1,3-glucan) on the pathogen surface to trigger host defense reactions . While some PRRs are constitutively expressed as surveillance molecules at certain levels in hemolymph, others are synthesized in response to the entry of microorganisms. Upon binding to their target molecules, PRRs undergo a conformational change required for recruiting plasma factors and hemocytes that eliminate the invading pathogens , finally leading to PPO activation and synthesis of immune effectors (e.g. AMPs) via immune signaling pathways . The H. armigera immunotranscriptome has at least 41 PRR transcripts, including 9 PGRPs, 5 βGRPs, 24 CTLs, and 3 galectins (in Additional file 6: Table S5, and in Additional file 7: Table S6).
Some thioester-containing proteins (TEPs), belonging to C3/α-2 macroglobulin superfamily, harbor an intrachain β-cysteinyl-γ-glutamyl thioester bond, which is used for self- or non-self-identification . α2-macroglobulin is a broad-spectrum protease inhibitor using a trapping mechanism through its thioester linkage. During the immune responses, complement factors C3, C4, and C5 bind to self/non-self surfaces via the same bond. In mosquitos, TEP family presents extensive gene duplication, with 15 genes in A. gambiae, 8 in A. aegypti, and 10 in C. pipiens . In A. gambiae, TEP1 played the key role in blocking Plasmodium oocyst proliferation in the midgut . Helicoverpa, however, has a relatively small TEP family with only four members which are homologous to B. mori TEPs (Additional file 1: Figure S6, in Additional file 6: Table S5). Since H. armigera TEP4 contains the signature motif GCGEQ, it is probably involved in immune surveillance against pathogens or parasites.
Based on sequence homology search, we have also identified other H. armigera PRRs encoded by 10 scavenger receptors class B (ScR-Bs), 1 hemocytin, 1 hemolin, 1 DSCAM, 1 Draper, 1 RhoL, 1 Rap1, 1 enabled, 1 Fascin, 1 Scar, 1 Dizzy, 1 TM9SF4, 1 Integrin, 1 Noduler, and 4 Eater, which are involved in phagocytosis, encapsulation, and various other cellular immune responses (in Additional file 6: Table S5) [46,47]. The phylogenetic analysis of the ScR-Bs (in Additional file 1: Figure S7) demonstrated that H. armigera and B. mori ScR-Bs were most likely to form the orthologous pairs. The H. armigera ScR-B9 might function in phagocytosis of the apoptotic bodies, since it is orthologous to D. melanogaster Croquemort (CG4280). The presence of a large number of PRR genes in the immunotranscriptome of H. armigera, indicated their involvement in the immune response to the infection of B. bassiana and E. cloacae.
Genes involved in melanization and extracellular signal modulation
Upon PRR binding, the invading microbes may lead to the PPO activation. As key enzymes in the melanization reaction, the copper-containing phenoloxidases are activated from its inactive zymogens by specific proteolytic cleavage . Active phenoloxidases catalyze the conversion of monophenols to o-diphenols and further oxidation of o-diphenols to quinones, which then polymerize to form melanin that entraps and kills microbial pathogens and parasites [8,48]. Except for mosquitoes, most insects had one to three PPOs . Similar to other lepidopteran species, two PPO transcripts were identified in the H. armigera transcriptome. The H. armigera PPO1 and PPO2 share 80.3% and 79% of identity in amino acid sequence to the B. mori PPO1 and PPO2 (in Additional file 6: Table S5), respectively. The phylogenetic analysis establishes the orthologous groups of HaPPO1-BmPPO1-MsPPO1 and HaPPO2-BmPPO2-MsPPO2 (in Additional file 1: Figure S8).
Serine proteases (SPs) and their non-catalytic homologs (SPHs) constitute one of the largest protein family in insects that are functionally involved in various physiological processes such as digestion, development, hemolymph clotting, and immune responses . Hemolymph SP cascades, comprised of multiple clip-domain SPs (cSPs), are essential for proteolytic activation of PPO and the processing of Toll pathway ligand, spätzle, which will be discussed in detail in the next section. In total, we have identified 46 SPs and 7 SPHs in H. armigera transcriptome (in Additional file 8: Table S7). Among them, 41 have complete sequences containing entire open reading frames (ORFs). Apart from these cSPs, five other SPs (SP41-44, 47) were shown to contain the structural modules important for protein-protein interactions, such as LDLA, Sushi, and SRCR. Detailed analysis of SPs revealed the orthology among H. armigera SP42, M. sexta HP14 and D. melanogaster modular SP. The respective SPs from M. sexta and D. melanogaster apparently are the first key enzymes triggering PPO activation and Spätzle processing [51,52]. Thus we postulated a similarly important role for H. armigera SP42, which needs further experimental verification.
The PPO activation is delicately regulated by serpins. Serpins are a superfamily of proteins, which have around 350–400 amino acid residues and the conserved tertiary structure comprising nine α-helices and three β-sheets . Most of them inhibit SP-mediated processes by forming covalent complexes with their cognate SPs. The P1 residue at an exposed reactive site loop (RSL) near the carboxyl terminus determines its inhibitory selectivity. Sequence similarity search of H. armigera unigenes yields 22 serpin transcripts, most of which have putative orthologs in B. mori (Figure 6B, in Additional file 6: Table S5). Of these, sixteen serpins have a predicted signal peptide for secretion, while the rest (Serpin2, 4, 10, 12, 16, 21) are assumed to be cytosolic (in Additional file 6: Table S5). H. armigera Serpin1-6 are highly similar to their M. sexta and B. mori orthologs, suggestive of their essential roles in regulating innate immunity . One interesting feature of M. sexta and B. mori Serpin1 and 2 is the mutual exclusive use of alternative exons to generate RSL variants on the same framework. While initial analysis suggests the same mechanism in H. armigera (data not shown), a detailed study of Serpin1 and 2 gene sequences will reveal the scope of this phenomenon. 1:1 orthologous relationship between H. armigera, B. mori, and M. sexta Serpins3, 4 and 5 allowed us to make predictions about the putative functions of these HaSerpins as summarized below [61,62]. D. melanogaster Spn27A, M. sexta serpin-3 and B. mori Serpin3 are negative regulators of the melanization cascade. In the phylogenetic tree, H. armigera Serpin3 was located in the same clade with them and also contained the conserved P2-P2′ sequence (NK-FG) (Figure 6, in Additional file 1: Figure S10). Serpin3 may negatively regulate melanization and the activation of Spätzle, while Serpin4 and 5 might inhibit hemolymph cSPs regulating immune response in general (Figure 6, in Additional file 1: Figure S10). We predict their proteolytic cleavage site based on the multiple sequence alignment of the hypervariable RSLs of the H. armigera serpins (in Additional file 1: Figure S10). Serpin10 and 14 lack a conserved carboxyl-terminus and may not fold properly with a long insert separating the partial serpin domain. Serpin1, 3 ~ 6, 8, 9, 13, 21, and 22 (with Arg or Lys located at the predicted P1 position) may inhibit trypsin-like enzymes, while serpin-2, 7, 11, 12, and 17 ~ 20 (with Phe, Tyr, or Leu located at the putative P1 site) are anticipated to inhibit chymotrypsin-like SPs. Only three serpins, Serpin12, 15, and 20, contain Val and Met at the P1 position to control elastase-like enzymes. These predicted inhibitory activities and selectivity features surely need further experimental data to validate.
The immune signaling pathways and effector gene families
The insect IMD pathway, analogous to the mammalian TNF signaling pathway, is critical for combating Gram-negative bacteria . IMD protein itself is similar to mammalian TNF receptor as both contain a homologous death domain. Activation of IMD leads to the proteolytic cleavage of a Rel family transcription activator, Relish. The released active domain of Relish is then translocated into the nucleus to exert its function . The IMD pathway components, including IMD, Dredd, FADD, TAK1, Tab2, IKK-β, IKK-γ and Relish were all identified in the H. armigera transcriptome strongly suggesting the conservation of IMD mediated immunity in cotton bollworm (in Additional file 6: Table S5). Additionally, two members of IAP-2 were also identified. In D. melanogaster, IAP-2 acts as a modulator of IMD pathway aiding the nuclear localization of Relish .Besides Toll and IMD pathways, JNK and JAK/STAT pathways are also implicated to play a role in insect immunity. Activated by TAK1 kinase, JNK pathways function through transcription factor Relish resulting in production of AMPs in both A. gambiae and D. melanogaster . One copy of JNK, Hem, Jun and Fos genes of the JNK pathway were all identified in H. armigera transcriptome (in Additional file 6: Table S5). In Drosophila, JAK/STAT pathway controls the transcription of some immune molecules including TEPs, promotes phagocytosis and participates in antiviral response . Based on homology to D. melanogaster JAK/STAT pathway members, cytokine receptor Domeless, Hopscotch, and STAT were present as single copy identified in H. armigera, indicating the conservation of this pathway (in Additional file 6: Table S5). The negative regulators of JAK/STAT pathway, SOCS and PIAS, were found in the transcriptome as well. Taken together, all these four pathways-Toll, IMD, JNK and JAK/STAT, coordinates in relaying the ‘danger’ signal ultimately stimulating the production of immune effector molecules (Figure 8).
Most AMPs are low molecular weight (<30 kDa) proteins forming amphiphilic α-helix, hairpin-like β-sheet, or structures stabilized by disulfide bonds . Induction of AMPs in fat body and hemocytes is thoroughly studied in lepidopteran insects including M. sexta and B. mori [28,31]. Analysis of the H. armigera transcriptome revealed unigenes homologous to AMPs, including 4 lysozymes, 5 cecropins, 4 moricin, 3 gloverins, 1 defensin, 1 gallerimycin, and 1 lebocin (in Additional file 6: Table S5). The number of annotated AMPs is much lower than that of B. mori  and M. sexta , probably due to the lack of genome information. While cecropins and defensins are found in various insects, moricins and gloverins are only reported to be present in Lepidoptera. Furthermore, four putative lysozyme transcripts are present in H. armigera transcriptome (in Additional file 1: Figure S11). Lysozymes are 14–16 kDa enzymes that hydrolyze peptidoglycans in bacterial cell wall . All four H. armigera lysozymes belong to the c-type lysozymes with four disulfide bridges, and lysozymes 1–3 contain the conserved active site including the Asp and Glu residues. Reactive oxygen and nitrogen species (ROS and RNS) played an important role in killing microbes and parasites . However, the overproduction of free radicals can be toxic to the host cell, therefore a tight regulation of ROS/RNS concentration is required. Superoxide dismutase (SOD), glutathione S-transferase (GST), or peroxiredoxin converts ROS/RNS to less toxic product. Here, we have annotated some of genes, including 5 GSTs, 3 SODs, and 6 peroxiredoxin transcripts (in Additional file 6: Table S5), indicating that free radicals produced by ROS and RNS are likely to be components of the defense system in H. armigera. Further experiments are needed to verify suggested role of the immune signaling pathways and effectors in H. armigera as described above.
Expression profile of immunity-related genes in the fat body after fungal and bacterial infections
Humoral defense proteins are mainly produced and secreted by fat body, an insect tissue analogous to liver and adipose tissue of mammals in terms of immunity and metabolism. For PRR genes, RNA-seq based transcriptome analysis indicated that short PGRPs (PGRP-SA1, -SA2, -SB1, -SB2, and -SD) were induced in fat body after the fungal challenge. Similarly, transcription of ScR-B1, B4, B8, TEP2, TEP3, CTL9 ~ 11, 18, 20 ~ 24 and galectins is up-regulated. This change is not limited to PRRs, some signal transducers and regulators are also induced. Cluster analysis of expression patterns revealed that cSP5, 8 ~ 10, cSPH11, 12, Serpin6, 14, 20, 21 mRNAs became much more abundant in fat body after the fungal challenge whereas higher levels of cSP2 ~ 4 and cSP6 transcripts were detected in hemocytes. The pronounced up-regulation of Toll pathway genes (e.g. Spz2, Pelle, Toll-8, Toll-9, Toll-14) occurred in fat body elicited by B. bassiana, suggesting their involvement in antifungal immunity. The gloverin1, lysozyme2, and lysozyme4 mRNA levels increased most dramatically to defend the fungal infection. Notably, most of these genes did not change sharply at transcript level in the other experiment conditions, indicating that fat body plays a major role in the antifungal response.
Bacterial challenge also resulted in the activation of a broad spectrum of immunity-related genes. These include PRRs such as CTLs (CTL2, 12, 14, 16, 17, and 26), signaling modulators (cSP3, Serpin1, 4, and 9), and effector molecules (gloverin1, gloverin2, Lebocin, and Moricin4). The mRNA abundances of FADD, Fos, JNK significantly increased in fat body elicited by E. cloacae, as evidenced by 12 genes in cluster 4, which were significantly up-regulated (fold change > 2, P < 0.001) in response to such bacterial infection. This number was apparently lower than that elicited by fungi, which is 24 (cluster 3, Figure 1C). Although more immunity-related genes were elicited by fungi (67) than bacteria (42), there is an overlap between them, including up-regulated cSP8, lysozyme3, gloverin1 and hemolin, and down-regulated βGRP3 (in Additional file 1: Figure S12).
Expression profile of immunity-related genes in hemocytes in response to fungal and bacterial challenges
Hemocyte mediated defense reactions such as phagocytosis and encapsulation also depend on the recognition of foreign pathogens to trigger downstream immune signaling pathways. The expression profiles of immunity-related genes in hemocytes elicited by the pathogens were highly different from those in fat body. To better understand the complexity of responses of hemocytes against pathogens, the expression profile of immunity-related genes elicited by fungi and bacteria was compared. Gene cohorts of cluster 1 include 21 such genes induced in B. bassiana induced hemocytes, including βGRP1, βGRP2, hemolin, cSP2, cSP6, Serpin2, Serpin5, lysozyme1, lysozyme3 and so on. Immune activation by E. cloacae is also noticeable as demonstrated by induction of PRRs (PGRP-SA2, TEP1) and immune modulators (MyD88, Toll-5, Cactus, Toll-12, Serpin3, 11, 21).
Thirty-three immunity-related genes in cluster 2 were significantly suppressed in B. bassiana elicited hemocytes. RNA-seq based transcriptome analysis indicated that long PGRPs (PGRP-LA, -LB, -LC, -LD) were strongly suppressed by B. bassiana but induced by E. cloacae. Quantitative real-time PCR also revealed that the βGRP1 mRNA level dramatically increased in hemocytes elicited by B. bassiana. In contrast, CTL mRNA levels did not change significantly in hemocytes. The transcript levels of PPOs were high in hemocytes, but significantly decreased in response to fungal challenge. Likewise, components of the IMD pathway (FADD, IAP2a, 2b, IKKβ, IKKγ, Dredd, Relish, UBc13, TAK1, IMD, Tab2), JNK and JAK-STAT pathways (Domeless, PIAS, Hopscotch, STAT, Fos, JNK, Hem, SOCS) were all dramatically suppressed in hemocytes in response to the fungal infection but moderately induced by bacterial infection, as validated by the quantitative real-time PCR assay of Domeless, JAK, and STAT transcripts. This indicated that signal transduction in hemocytes responds differently to bacterial and fungal infection. More immunity-related genes were regulated by fungi (60) than bacteria (18) in hemocytes, and most of them tended to be suppressed (in Additional file 1: Figure S12).
Although fat body and hemocytes both are important organs involved in the immune response, our transcriptome analysis revealed that more immunity-related genes are induced in fat body (clusters 3 and 4) than hemocytes (cluster 1). The expression of immune recognition and modulation molecules was mainly controlled in fat body, while genes involved in signal transduction, such as JAK-STAT and IMD, were regulated more drastically in hemocytes (Figure 9). In Toll pathway, hemocytes regulated Toll-1, Toll-5, Toll-12, Tollip, Cactus, and MyD88, but more other genes were regulatedin fat body. Additionally, quantitative real-time PCR indicated that the transcription of putative PGRP-SA1, PGRP-LB, Toll-14 and Spz2 genes were elevated in fat body upon B. bassiana infection, while the mRNA levels of defensin, moricin and gloverin1 were up-regulated in hemocytes (Figure 9), which is consistent with the RNA-seq results.
While many immunological studies are focused on B. mori and M. sexta, H. armigera, as a notorious lepidopteran pest affecting world agriculture, has not yet been systematically investigated on aspects of immunity until now. In order to develop and apply better pest control methods, we used high throughput sequencing to analyze the transcriptome of hemocytes and fat body in response to bacterial and fungal infection. The 233 immunity-related transcripts discovered are involved in recognition, signal transduction and modulation, and execution mechanisms. This number is certainly going to increase when its genomic information is available, based on the information acquired from this study, we can conclude that the repertoire of immunity-related genes is well conserved in H. armigera. Moreover, we find that these genes displayed distinct mRNA profiles in response to B. bassiana and E. cloacae infection. Their products may constitute an integrated pathway of antifungal immunity, which provide valuable clues for further functional studies in H. armigera responding to B. bassiana, and pave the way for using H. armigera as a model system to study fungal-host interaction . We have achieved another goal of using the transcriptome dataset to reveal tissue-specific gene expression patterns. The global patterns of hemocytes and fat body remarkably differ from each other after infection by the two microbes. Furthermore, it is demonstrated that not only immunity-related genes but also metabolism genes are under the regulation of pathogen infection.
Taken together, this study provided a global view of host defense gene expression profiles in response to bacterial and fungal challenges in a non-model insect. The results could benefit future in-depth study on the role of candidate genes involved in anti-fungal and -bacterial immunity, thereby helping improve the understanding of host-pathogen interactions and evolutionary history of immunogenetics from insects to mammals.
Insect rearing, bacterium and fungus culture
A colony of H. armigera was obtained from Dr. Qi-Lian Qin at Institute of Zoology, Chinese Academy of Sciences, and maintained in the laboratory. Larvae were reared on an artificial diet at 28 ± 1°C under a photoperiod of 14:10 (l: d) and 70% relative humidity. Male and female pupae were placed in separate glass dishes under the same conditions for eclosion. Adults were fed with 10% honey solution. The colony of H. armigera in the laboratory was maintained for over five generations before experiments. E. cloacae strain (strain no. 1.2022) was obtained from China General Microbiological Culture Collection Center. The bacteria at stationary phase in LB broth (A600 ≈ 2.0) were diluted four times (A600 ≈ 0.5) and used for septic injury. B. bassiana strain ARSEF 2860, was cultured on potato dextrose agar (PDA) plates at 25°C and 80% humidity .
Septic injury and RNA sample preparation
At day 2 after molting, approximately 120 fifth instar larvae from the same batch were divided into three groups. Each larva in the treated groups was inoculated in hemocoel with 2 μl of the B. bassiana suspension (1 × 107 conidia/ml), E. cloacae suspension (A600 ≈ 0.5), or sterile PBS using a microinjector. Then all the larvae were kept under the same rearing conditions. The samples were collected from E. cloacae infected larvae at 6 h and B. bassiana infected larvae at 48 h, respectively, at which time larvae started to die. Hemolymph samples were individually collected into ice-cold PBS containing 0.1% 1-phenyl-2-thiourea. Suspensions were pooled in a 1.5 ml Eppendorf tubes and centrifuged at 500 g for 5 min at 4°C to collect hemocytes. Fat body tissues were removed from larvae under a dissection microscope. Control samples from the PBS-injected larvae were prepared at 6 h and 48 h in the same way, and RNA was combined in equal amount.
Samples of hemocytes and fat body were homogenized using a motor-driven pellet pestle mixer (Kontes, USA) and lysed in 1.0 ml of TRI Reagent (Sigma, USA). Total RNA was extracted using RNeasy Mini Kit (Qiagen, USA) according to the manufacturer’s instructions. RNA concentrations were determined on a Nanodrop ND-2000 spectrophotometer (NanoDrop products, USA). RNA integrity was verified on Agilent 2100 BioAnalyzer (Agilent, USA).
Day 2 fifth instar H. armigera larvae (24 in each group) were individually challenged with B. bassiana or E. cloacae. Insects in the treatment groups and control group of buffer injection were maintained in individual containers and fed on the standard diet. Dead larvae were counted at different time intervals for calculation of mortality rate. The survival curves were compared using the Kaplan-Meier and Cox’s proportional hazards model to assess variables that affected survival . The threshold of p value was adjusted by the Bonferroni correction. GraphPad software was used in all statistical analysis. The experiment was repeated three times.
Library construction and Illumina sequencing
Ten μg of total RNA from each treatment or control group was used to enrich poly(A) mRNA using oligo(dT) magnetic beads (Invitrogen, USA). In the next step, paired-end RNA-seq libraries were prepared by following the Illumina’s library construction protocol. The challenged and control libraries were sequenced on Illumina HiSeq 2000 platform (Illumina, USA) in the Beijing Institutes of Biological Sciences (Chinese Academy of Sciences). FASTQ files of raw-reads were produced and sorted out by barcodes for further analysis.
Assembly of transcriptomes and functional classification
Prior to assembly, 2 × 100 bp paired-end raw reads from each cDNA library were processed to remove adaptors, low quality sequences (Q < 20), and reads contaminated with microbes using Trimmomatic package . The FastQC package was used to verify the quality of resulting trimmed and filtered reads. The clean reads were de novo assembled to produce contigs using Trinity, the short reads assembling program [19,20] using default parameters. One set of transcript contigs was assembled for each library. To increase the transcriptome coverage, second instar larvae were injected with the bacteria, fungi or PBS (Xiong and Zou, unpublished data). Whole body mRNA samples were prepared and sequenced in the same way. Clean reads from the nine samples were assembled into the tenth dataset, which had the most full-length CDSs and splicing variants. To reduce redundancy, CAP3  was used to link overlapping transcripts into larger contigs. The final collection of transcripts was made up from consensus cluster sequences and singletons. The final collection of reference transcripts (or unigenes) was made up of consensus cluster sequences and singletons.
Non-coding RNAs were obtained using homologous search NONCODE . For functional annotations, we searched all the assembled sequences against the NCBI non-redundant (NR) sequence database using BLASTX, with an E-value cut-off of 1e−5. An in-house Perl script was used to identify putative coding regions from the sequences, and translate them into amino acid sequences for further analysis. To determine enriched GO terms, significantly up- and down-regulated genes were subjected to gene ontology enrichment analysis through BLAST2GO against GO database [21,25]. Finally, KEGG analysis was performed to discover the enriched pathways and metabolic networks via online tool KAAS at the KEGG website . In addition, GO annotations were performed for all unigenes with BLAST2GO program according to the GO association found by a BLASTX search against the NR database. For GO enrichment analysis, the corrected p value less than 0.05 were considered as significantly enriched in DETs. For KEGG pathway enrichment analysis, pathways with Q value less than 0.05 was regarded as statistically significant.
Identification and hierarchical clustering of differentially expressed transcripts (DETs)
Transcripts abundance were quantified as follows: First, aligning the filtered raw paired-end reads from each of the six tissue-specific library to the corresponding non-redundant H. armigera unigenes using Bowtie . Next, RSEM  package was applied for calculation of normalized gene expression value FPKM (an estimate of the number of sequence reads that are originated from a given transcript, accounting for the possibility that a single read could align to multiple transcripts). Subsequently, DETs between the control and treatment libraries were calculated based on the significance level (p value < 0.001) using DEGseq package in R environment . These transcripts with expression difference more than 2 or less than 0.5 folds were considered as up- and down-regulated genes, respectively. Hierarchical clustering of DET was performed using heatmap.2 function within gplots package in R environment, where the average method was used to compute the hierarchical clustering and Pearson correlation method was adopted to calculate the distance between both rows and columns. For principal component analysis (PCA), PCs, PC variances, and PC scores were calculated using eigen function in R language .
Comparative analysis of immunity-related genes from H. armigera transcriptome
The available immunity-related gene sequences from other model insect species were used as templates to search the H. armigera unigene sequence collection. The reference insect species include D. melanogaster , A. gambiae , A. mellifera , B. mori , and T. castaneum . The potential candidates of H. armigera immunity-related genes were confirmed by searching the BLASTX algorithm against the Nr database using a cut-off E-value of 0.1. The deduced amino acid sequences were analyzed by Pfam (http://pfam.xfam.org/) and SMART (http://smart.embl-heidelberg.de/) to detect conserved domain structures required for specific functions. Signal peptide predictions were determined by SignalP4.0.
H. armigera sequences were aligned with their homologues from other representative insect species from holometabolous insect orders, including D. melanogaster, A. gambiae, A. mellifera, T. castaneum, and B. mori. The sequences were retrieved from NCBI, Immunodb, or Ensembl. Multiple sequence alignments were performed by CLUSTALX software with BLOSUM series of weight matrices . Phylogenetic trees were constructed by the neighbor-joining method with statistical analysis by the bootstrap method, using 1000 repetitions .
Quantitative real-time PCR analysis
To validate the expression profile from RNA-seq results, specific primers were designed to perform a quantitative real-time PCR analysis of selected immunity-related genes (in Additional file 9: Table S8) . The annealing temperatures of the primers were controlled at about 62°C. 2 μg of total RNA was treated with DNase I (Invitrogen, USA) to remove contaminating genomic DNA, and then used for cDNA synthesis by M-MLV reverse transcriptase (Promega, USA). Ribosomal protein S3 was used as an internal standard to normalize cDNA templates. The quantitative real-time PCR was performed on a MX3000P system (Stratagene, USA) using SYBR green PCR Master Mix (Tiangen, China), according to the manufacturer’s instructions. The thermal cycling conditions were: 94°C, 5 s; 59°C, 20 s; 72°C, 20 s. Quantitative real-time PCR data were collected and exported to EXCEL for analysis.
Sequence data submission
All the raw sequence data and the final Helicoverpa transcriptome described in the manuscript are available under the BioProject accession No. PRJNA264881. Illumina sequence reads have been deposited at NCBI SRA database under the following accession numbers (FB_Mock: Sample: SAMN03120571, Experiment: SRX735419, Reads: SRR1630741; FB_Bb: Sample: SAMN03144202, Experiment: SRX743611, Reads: SRR1630794; FB_Ec: Sample: SAMN03144203, Experiment: SRX743612, Reads: SRR1630796; HC_Mock: Sample: SAMN03144201, Experiment: SRX743613, Reads: SRR1630797; HC_Bb: Sample: SAMN03144199, Experiment: SRX743614, Reads: SRR1630798; HC_Ec: Sample: SAMN03144200, Experiment: SRX743615, Reads: SRR1630799). H. armigera unigenes are available under NCBI TSA accession number affiliated with the same BioProject.
This work was supported by Strategic Priority Research Program of the CAS (No. XDB11030600), National Basic Research Program of China (No. 2014CB138405), National Natural Science Foundation of China (No. 31472008, 31272367), Open research Fund Program of State Key Laboratory of Integrated Pest Management (Chinese IPM1304) to Z.Z., and National Institutes of Health grant GM58634 to H.J.
- Tay WT, Soria MF, Walsh T, Thomazoni D, Silvie P, Behere GT, et al. A brave new world for an old world pest: Helicoverpa armigera (Lepidoptera: Noctuidae) in Brazil. PLoS One. 2013;8(11):e80134.PubMed CentralPubMedView ArticleGoogle Scholar
- Wang CS, Feng MG. Advances in fundamental and applied studies in China of fungal biocontrol agents for use against arthropod pests. Biol Control. 2014;68:129–35.View ArticleGoogle Scholar
- Wraight SP, Ramos ME, Avery PB, Jaronski ST, Vandenberg JD. Comparative virulence of Beauveria bassiana isolates against lepidopteran pests of vegetable crops. J Invertebr Pathol. 2010;103(3):186–99.PubMedView ArticleGoogle Scholar
- Gillespie JP, Kanost MR, Trenczek T. Biological mediators of insect immunity. Annu Rev Entomol. 1997;42:611–43.PubMedView ArticleGoogle Scholar
- Lemaitre B, Hoffmann J. The host defense of Drosophila melanogaster. Annu Rev Immunol. 2007;25:697–743.PubMedView ArticleGoogle Scholar
- Ferrandon D, Imler JL, Hetru C, Hoffmann JA. The Drosophila systemic immune response: sensing and signalling during bacterial and fungal infections. Nat Rev Immunol. 2007;7(11):862–74.PubMedView ArticleGoogle Scholar
- Hetru C, Hoffmann JA. NF-kappaB in the immune response of Drosophila. Cold Spring Harb Perspect Biol. 2009;1(6):a000232.PubMed CentralPubMedView ArticleGoogle Scholar
- Kanost MR, Jiang H, Yu X-Q. Innate immune responses of a lepidopteran insect Manduca Sexta. Immunol Rev. 2004;198(1):97–105.PubMedView ArticleGoogle Scholar
- Gottar M, Gobert V, Matskevich AA, Reichhart JM, Wang C, Butt TM, et al. Dual detection of fungal infections in Drosophila via recognition of glucans and sensing of virulence factors. Cell. 2006;127(7):1425–37.PubMed CentralPubMedView ArticleGoogle Scholar
- Yassine H, Kamareddine L, Osta MA. The mosquito melanization response is implicated in defense against the entomopathogenic fungus Beauveria bassiana. PLoS Pathog. 2012;8(11):e1003029.PubMed CentralPubMedView ArticleGoogle Scholar
- Zou Z, Shin SW, Alvarez KS, Kokoza V, Raikhel AS. Distinct melanization pathways in the mosquito Aedes aegypti. Immunity. 2010;32(1):41–53.PubMedView ArticleGoogle Scholar
- Kim CH, Kim SJ, Kan H, Kwon HM, Roh KB, Jiang R, et al. A three-step proteolytic cascade mediates the activation of the peptidoglycan-induced toll pathway in an insect. J Biol Chem. 2008;283(12):7599–607.PubMedView ArticleGoogle Scholar
- Zhao Z, Wu G, Wang J, Liu C, Qiu L. Next-generation sequencing-based transcriptome analysis of Helicoverpa armigera Larvae immune-primed with Photorhabdus luminescens TT01. PLoS One. 2013;8(11):e80146.PubMed CentralPubMedView ArticleGoogle Scholar
- Wang Q, Liu Y, He HJ, Zhao XF, Wang JX. Immune responses of Helicoverpa armigera to different kinds of pathogens. BMC Immunol. 2010;11:9.PubMed CentralPubMedView ArticleGoogle Scholar
- Wang JL, Zhang Q, Tang L, Chen L, Liu XS, Wang YF. Involvement of a pattern recognition receptor C-type lectin 7 in enhancing cellular encapsulation and melanization due to its carboxyl-terminal CRD domain in the cotton bollworm, Helicoverpa armigera. Dev Comp Immunol. 2014;44(1):21–9.PubMedView ArticleGoogle Scholar
- Tian YY, Liu Y, Zhao XF, Wang JX. Characterization of a C-type lectin from the cotton bollworm, Helicoverpa armigera. Dev Comp Immunol. 2009;33(6):772–9.PubMedView ArticleGoogle Scholar
- Yang DQ, Su ZL, Qiao C, Zhang Z, Wang JL, Li F, et al. Identification and characterization of two peptidoglycan recognition proteins with zinc-dependent antibacterial activity from the cotton bollworm, Helicoverpa armigera. Dev Comp Immunol. 2013;39(4):343–51.PubMedView ArticleGoogle Scholar
- Shin SW, Bian G, Raikhel AS. A toll receptor and a cytokine, Toll5A and Spz1C, are involved in toll antifungal immune signaling in the mosquito Aedes aegypti. J Biol Chem. 2006;281(51):39388–95.PubMedView ArticleGoogle Scholar
- Haas BJ, Papanicolaou A, Yassour M, Grabherr M, Blood PD, Bowden J, et al. De novo transcript sequence reconstruction from RNA-seq using the Trinity platform for reference generation and analysis. Nat Protoc. 2013;8(8):1494–512.PubMedView ArticleGoogle Scholar
- Grabherr MG, Haas BJ, Yassour M, Levin JZ, 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.PubMed CentralPubMedView ArticleGoogle Scholar
- Aparicio G, Gotz S, Conesa A, Segrelles D, Blanquer I, Garcia JM, et al. Blast2GO goes grid: developing a grid-enabled prototype for functional genomics analysis. Studies Health Technol Informatics. 2006;120:194–204.Google Scholar
- Chen S, Yang P, Jiang F, Wei Y, Ma Z, Kang L. De novo analysis of transcriptome dynamics in the migratory locust during the development of phase traits. PLoS One. 2010;5(12):e15633.PubMed CentralPubMedView ArticleGoogle Scholar
- Wang L, Feng Z, Wang X, Wang X, Zhang X. DEGseq: an R package for identifying differentially expressed genes from RNA-seq data. Bioinformatics. 2010;26(1):136–8.PubMedView ArticleGoogle Scholar
- Ogata H, Goto S, Sato K, Fujibuchi W, Bono H, Kanehisa M. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 1999;27(1):29–34.PubMed CentralPubMedView ArticleGoogle Scholar
- Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, et al. Gene Ontology: tool for the unification of biology. Nat Genet. 2000;25(1):25–9.PubMed CentralPubMedView ArticleGoogle Scholar
- Wu J, Mao X, Cai T, Luo J, Wei L. KOBAS server: a web-based platform for automated annotation and pathway identification. Nucleic Acids Res. 2006;34(Web Server issue):W720–4.PubMed CentralPubMedView ArticleGoogle Scholar
- Zou Z, Evans JD, Lu Z, Zhao P, Williams M, Sumathipala N, et al. Comparative genomic analysis of the Tribolium immune system. Genome Biol. 2007;8(8):R177.PubMed CentralPubMedView ArticleGoogle Scholar
- Jiang HB. The biochemical basis of antimicrobial responses in Manduca sexta. Insect Sci. 2008;15(1):53–66.View ArticleGoogle Scholar
- Werner T, Liu G, Kang D, Ekengren S, Steiner H, Hultmark D. A family of peptidoglycan recognition proteins in the fruit fly Drosophila melanogaster. Proc Natl Acad Sci U S A. 2000;97(25):13772–7.PubMed CentralPubMedView ArticleGoogle Scholar
- Li L, Li YP, Song CX, Xiao M, Wang JL, Liu XS. Identification and functional characterization of a peptidoglycan recognition protein from the cotton bollworm, Helicoverpa armigera. Arch Insect Biochem Physiol. 2014;86(4):240–58.PubMedView ArticleGoogle Scholar
- Tanaka H, Ishibashi J, Fujita K, Nakajima Y, Sagisaka A, Tomimoto K, et al. A genome-wide analysis of genes and gene families involved in innate immunity of Bombyx mori. Insect Biochem Mol Biol. 2008;38(12):1087–110.PubMedView ArticleGoogle Scholar
- Chen K, Liu C, He Y, Jiang H, Lu Z. A short-type peptidoglycan recognition protein from the silkworm: expression, characterization and involvement in the prophenoloxidase activation pathway. Dev Comp Immunol. 2014;45(1):1–9.PubMedView ArticleGoogle Scholar
- Lee WJ, Lee JD, Kravchenko VV, Ulevitch RJ, Brey PT. Purification and molecular cloning of an inducible gram-negative bacteria-binding protein from the silkworm, Bombyx mori. Proc Natl Acad Sci U S A. 1996;93(15):7888–93.PubMed CentralPubMedView ArticleGoogle Scholar
- Wang Y, Sumathipala N, Rayaprolu S, Jiang H. Recognition of microbial molecular patterns and stimulation of prophenoloxidase activation by a beta-1,3-glucanase-related protein in Manduca sexta larval plasma. Insect Biochem Mol Biol. 2011;41(5):322–31.PubMed CentralPubMedView ArticleGoogle Scholar
- Fujita T, Matsushita M, Endo Y. The lectin-complement pathway–its role in innate immunity and evolution. Immunol Rev. 2004;198:185–202.PubMedView ArticleGoogle Scholar
- Vasta GR. Roles of galectins in infection. Nat Rev Microbiol. 2009;7(6):424–38.PubMed CentralPubMedView ArticleGoogle Scholar
- Chai LQ, Tian YY, Yang DT, Wang JX, Zhao XF. Molecular cloning and characterization of a C-type lectin from the cotton bollworm, Helicoverpa armigera. Dev Comp Immunol. 2008;32(1):71–83.PubMedView ArticleGoogle Scholar
- Wang JL, Liu XS, Zhang Q, Zhao HB, Wang YF. Expression profiles of six novel C-type lectins in response to bacterial and 20E injection in the cotton bollworm (Helicoverpa armigera). Dev Comp Immunol. 2012;37(2):221–32.PubMedView ArticleGoogle Scholar
- Yu XQ, Kanost MR. Immulectin-2, a lipopolysaccharide-specific lectin from an insect, Manduca sexta, is induced in response to gram-negative bacteria. J Biol Chem. 2000;275(48):37373–81.PubMedView ArticleGoogle Scholar
- Yu XQ, Kanost MR. Immulectin-2, a pattern recognition receptor that stimulates hemocyte encapsulation and melanization in the tobacco hornworm, Manduca sexta. Dev Comparative Immunol. 2004;28(9):891–900.View ArticleGoogle Scholar
- Osta MA, Christophides GK, Kafatos FC. Effects of mosquito genes on Plasmodium development. Science. 2004;303(5666):2030–2.PubMedView ArticleGoogle Scholar
- Tanji T, Ohashi-Kobayashi A, Natori S. Participation of a galactose-specific C-type lectin in Drosophila immunity. Biochem J. 2006;396:127–38.PubMed CentralPubMedView ArticleGoogle Scholar
- Blandin S, Levashina EA. Thioester-containing proteins and insect immunity. Mol Immunol. 2004;40(12):903–8.PubMedView ArticleGoogle Scholar
- Bartholomay LC, Waterhouse RM, Mayhew GF, Campbell CL, Michel K, Zou Z, et al. Pathogenomics of Culex quinquefasciatus and meta-analysis of infection responses to diverse pathogens. Science. 2010;330(6000):88–90.PubMed CentralPubMedView ArticleGoogle Scholar
- Levashina EA, Moita LF, Blandin S, Vriend G, Lagueux M, Kafatos FC. Conserved role of a complement-like protein in phagocytosis revealed by dsRNA knockout in cultured cells of the mosquito, Anopheles gambiae. Cell. 2001;104(5):709–18.PubMedView ArticleGoogle Scholar
- Strand MR. The insect cellular immune response. Insect Sci. 2008;15(1):1–14.View ArticleGoogle Scholar
- Vlisidou I, Wood W. Drosophila blood cells and their role in immune responses. FEBS J. 2015; in pressGoogle Scholar
- Li X, Ma M, Liu F, Chen Y, Lu A, Ling QZ, et al. Properties of Drosophila melanogaster prophenoloxidases expressed in Escherichia coli. Dev Comp Immunol. 2012;36(4):648–56.PubMedView ArticleGoogle Scholar
- Waterhouse RM, Kriventseva EV, Meister S, Xi Z, Alvarez KS, Bartholomay LC, et al. Evolutionary dynamics of immune-related genes and pathways in disease-vector mosquitoes. Science. 2007;316(5832):1738–43.PubMed CentralPubMedView ArticleGoogle Scholar
- Ross J, Jiang H, Kanost MR, Wang Y. Serine proteases and their homologs in the Drosophila melanogaster genome: an initial analysis of sequence conservation and phylogenetic relationships. Gene. 2003;304:117–31.PubMedView ArticleGoogle Scholar
- Ji C, Wang Y, Guo X, Hartson S, Jiang H. A pattern recognition serine proteinase triggers the prophenoloxidase activation cascade in the tobacco hornworm, Manduca sexta. J Biol Chem. 2004;279(33):34101–6.PubMedView ArticleGoogle Scholar
- Buchon N, Poidevin M, Kwon HM, Guillou A, Sottas V, Lee BL, et al. A single modular serine protease integrates signals from pattern-recognition receptors upstream of the Drosophila Toll pathway. Proc Natl Acad Sci U S A. 2009;106(30):12442–7.PubMed CentralPubMedView ArticleGoogle Scholar
- Christophides GK, Zdobnov E, Barillas-Mury C, Birney E, Blandin S, Blass C, et al. Immunity-related genes and gene families in Anopheles gambiae. Science. 2002;298(5591):159–65.PubMedView ArticleGoogle Scholar
- Zou Z, Lopez DL, Kanost MR, Evans JD, Jiang H. Comparative analysis of serine protease-related genes in the honey bee genome: possible involvement in embryonic development and innate immunity. Insect Mol Biol. 2006;15(5):603–14.PubMed CentralPubMedView ArticleGoogle Scholar
- Lee SY, Cho MY, Hyun JH, Lee KM, Homma KI, Natori S, et al. Molecular cloning of cDNA for pro-phenol-oxidase-activating factor I, a serine protease is induced by lipopolysaccharide or 1,3-beta-glucan in coleopteran insect, Holotrichia diomphalia larvae. Eur J Biochem. 1998;257(3):615–21.PubMedView ArticleGoogle Scholar
- Lee KY, Zhang R, Kim MS, Park JW, Park HY, Kawabata S, et al. A zymogen form of masquerade-like serine proteinase homologue is cleaved during pro-phenoloxidase activation by Ca2+ in coleopteran and Tenebrio molitor larvae. Eur J Biochem. 2002;269(17):4375–83.PubMedView ArticleGoogle Scholar
- Jiang H, Wang Y, Kanost MR. Pro-phenol oxidase activating proteinase from an insect, Manduca sexta: a bacteria-inducible protein similar to Drosophila easter. Proc Natl Acad Sci U S A. 1998;95(21):12220–5.PubMed CentralPubMedView ArticleGoogle Scholar
- Jiang H, Wang Y, Yu XQ, Zhu Y, Kanost M. Prophenoloxidase-activating proteinase-3 (PAP-3) from Manduca sexta hemolymph: a clip-domain serine proteinase regulated by serpin-1 J and serine proteinase homologs. Insect Biochem Mol Biol. 2003;33(10):1049–60.PubMedView ArticleGoogle Scholar
- An C, Ishibashi J, Ragan EJ, Jiang H, Kanost MR. Functions of Manduca sexta hemolymph proteinases HP6 and HP8 in two innate immune pathways. J Biol Chem. 2009;284(29):19716–26.PubMed CentralPubMedView ArticleGoogle Scholar
- Kanost MR. Serine proteinase inhibitors in arthropod immunity. Dev Comp Immunol. 1999;23(4–5):291–301.PubMedView ArticleGoogle Scholar
- Zhu Y, Wang Y, Gorman MJ, Jiang H, Kanost MR. Manduca sexta serpin-3 regulates prophenoloxidase activation in response to infection by inhibiting prophenoloxidase-activating proteinases. J Biol Chem. 2003;278(47):46556–64.PubMedView ArticleGoogle Scholar
- Zou Z, Picheng Z, Weng H, Mita K, Jiang H. A comparative analysis of serpin genes in the silkworm genome. Genomics. 2009;93(4):367–75.PubMed CentralPubMedView ArticleGoogle Scholar
- Michel T, Reichhart JM, Hoffmann JA, Royet J. Drosophila Toll is activated by Gram-positive bacteria through a circulating peptidoglycan recognition protein. Nature. 2001;414(6865):756–9.PubMedView ArticleGoogle Scholar
- Wang Y, Cheng T, Rayaprolu S, Zou Z, Xia Q, Xiang Z, et al. Proteolytic activation of pro-spatzle is required for the induced transcription of antimicrobial peptide genes in lepidopteran insects. Dev Comp Immunol. 2007;31(10):1002–12.PubMed CentralPubMedView ArticleGoogle Scholar
- Ooi JY, Yagi Y, Hu X, Ip YT. The Drosophila Toll-9 activates a constitutive antimicrobial defense. EMBO Rep. 2002;3(1):82–7.PubMed CentralPubMedView ArticleGoogle Scholar
- Kleino A, Valanne S, Ulvila J, Kallio J, Myllymaki H, Enwald H, et al. Inhibitor of apoptosis 2 and TAK1-binding protein are components of the Drosophila Imd pathway. EMBO J. 2005;24(19):3423–34.PubMed CentralPubMedView ArticleGoogle Scholar
- Osta MA, Christophides GK, Vlachou D, Kafatos FC. Innate immunity in the malaria vector Anopheles gambiae: comparative and functional genomics. J Exp Biol. 2004;207(Pt 15):2551–63.PubMedView ArticleGoogle Scholar
- Boutros M, Agaisse H, Perrimon N. Sequential activation of signaling pathways during innate immune responses in Drosophila. Dev Cell. 2002;3(5):711–22.PubMedView ArticleGoogle Scholar
- Bulet P, Hetru C, Dimarcq JL, Hoffmann D. Antimicrobial peptides in insects; structure and function. Dev Comp Immunol. 1999;23(4–5):329–44.PubMedView ArticleGoogle Scholar
- Yu KH, Kim KN, Lee JH, Lee HS, Kim SH, Cho KY, et al. Comparative study on characteristics of lysozymes from the hemolymph of three lepidopteran larvae, Galleria mellonella, Bombyx mori, Agrius convolvuli. Dev Comp Immunol. 2002;26(8):707–13.PubMedView ArticleGoogle Scholar
- Kumar S, Molina-Cruz A, Gupta L, Rodrigues J, Barillas-Mury C. A peroxidase/dual oxidase system modulates midgut epithelial immunity in Anopheles gambiae. Science. 2010;327(5973):1644–8.PubMed CentralPubMedView ArticleGoogle Scholar
- Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013;30(12):2725–9.PubMed CentralPubMedView ArticleGoogle Scholar
- Xiao G, Ying SH, Zheng P, Wang ZL, Zhang S, Xie XQ, et al. Genomic perspectives on the evolution of fungal entomopathogenicity in Beauveria bassiana. Sci Reports. 2012;2:483.Google Scholar
- Wang Y, Yang P, Cui F, Kang L. Altered immunity in crowded locust reduced fungal (Metarhizium anisopliae) pathogenesis. PLoS Pathog. 2013;9(1):e1003102.PubMed CentralPubMedView ArticleGoogle Scholar
- Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30(15):2114–20.PubMed CentralPubMedView ArticleGoogle Scholar
- Huang X, Madan A. CAP3: A DNA sequence assembly program. Genome Res. 1999;9(9):868–77.PubMed CentralPubMedView ArticleGoogle Scholar
- Bu D, Yu K, Sun S, Xie C, Skogerbo G, Miao R, et al. NONCODE v3.0: integrative annotation of long noncoding RNAs. Nucleic Acids Res. 2012;40(Database issue):D210–5.PubMed CentralPubMedView ArticleGoogle Scholar
- Langmead B, Trapnell C, Pop M, Salzberg SL. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009;10(3):R25.PubMed CentralPubMedView ArticleGoogle Scholar
- Anderson WA, Banerjee U, Drennan CL, Elgin SC, Epstein IR, Handelsman J, et al. Science education. Changing the culture of science education at research universities. Science. 2011;331(6014):152–3.PubMedView ArticleGoogle Scholar
- Evans JD, Aronstein K, Chen YP, Hetru C, Imler JL, Jiang H, et al. Immune pathways and defence mechanisms in honey bees Apis mellifera. Insect Mol Biol. 2006;15(5):645–56.PubMed CentralPubMedView ArticleGoogle Scholar
- Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG. The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res. 1997;25(24):4876–82.PubMed CentralPubMedView ArticleGoogle Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. 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.