Antimicrobial peptide-like genes in Nasonia vitripennis: a genomic perspective
BMC Genomics volume 11, Article number: 187 (2010)
Antimicrobial peptides (AMPs) are an essential component of innate immunity which can rapidly respond to diverse microbial pathogens. Insects, as a rich source of AMPs, attract great attention of scientists in both understanding of the basic biology of the immune system and searching molecular templates for anti-infective drug design. Despite a large number of AMPs have been identified from different insect species, little information in terms of these peptides is available from parasitic insects.
By using integrated computational approaches to systemically mining the Hymenopteran parasitic wasp Nasonia vitripennis genome, we establish the first AMP repertoire whose members exhibit extensive sequence and structural diversity and can be distinguished into multiple molecular types, including insect and fungal defensin-like peptides (DLPs) with the cysteine-stabilized α-helical and β-sheet (CSαβ) fold; Pro- or Gly-rich abaecins and hymenoptaecins; horseshoe crab tachystatin-type AMPs with the inhibitor cystine knot (ICK) fold; and a linear α-helical peptide. Inducible expression pattern of seven N. vitripennis AMP genes were verified, and two representative peptides were synthesized and functionally identified to be antibacterial. In comparison with Apis mellifera (Hymenoptera) and several non-Hymenopteran model insects, N. vitripennis has evolved a complex antimicrobial immune system with more genes and larger protein precursors. Three classical strategies that are likely responsible for the complexity increase have been recognized: 1) Gene duplication; 2) Exon duplication; and 3) Exon-shuffling.
The present study established the N. vitripennis peptidome associated with antimicrobial immunity by using a combined computational and experimental strategy. As the first AMP repertoire of a parasitic wasp, our results offer a basic platform for further studying the immunological and evolutionary significances of these newly discovered AMP-like genes in this class of insects.
AMPs constitute essential components of innate immunity to rapidly respond to diverse microbial pathogens [1, 2]. As key effectors, AMPs directly kill invaders by acting either as pore-formers or metabolic inhibitors . Most AMPs are cationic polypeptides, usually smaller than 100 amino acids in length, with enormous sequence diversity. Based on structural characteristics, these molecules are roughly divided into three different groups : 1) Linear peptides free of cysteines, often forming an α-helical conformation with amphiphilic surface, such as insect cecropins and mellitins, amphibian magainins and arachnidian meucins [5, 6]; 2) Cysteine-rich peptides. The most representative of this group are insect-derived defensins, which belong to the CSαβ polypeptide superfamily and are conserved across the kingdoms of eukaryotes [7, 8], even in bacteria [9, 10]; 3) Peptides with unusual bias in certain amino acids, such as glycine-rich tenecin3 , proline-rich apidaecin , histidine-rich histatin  and tryptophan-rich indolicidin . AMPs in this group are unstructured and often enter into bacterial cells to inhibit different metabolic targets [3, 15].
Insects, as a rich source of AMPs, attract great attention of scientists in both understanding the basic biology of immune system and searching molecular templates for anti-infective drug design. So far, about 200 such peptides have been identified from insects. In Drosophila melanogaster, there are 20 AMPs characterized . In the past ten years, a number of genome sequencing projects have promoted application of computational approaches for the discovery of components involved in innate immunity of a given species. Some examples include D. melanogaster, Anopheles gambiae, A. mellifera, Tribolium castaneum and Bombyx mori[17–21]. This will undoubtedly extend our knowledge in basic biology of insect innate immunity system. Furthermore, these data also provide new clues for elucidation of immunological adaptation of insects to environmental changes. For example, analysis of the complete immune system of A. gambiae results in the discovery of a marked deficit of orthologues due to gene loss and excessive expansion of some specific genes, which is significantly different from that of D. melanogaster. This possibly reflects differential selective pressures from different pathogens between these two species. These results also facilitate further in-depth analysis of the mosquito immune system's impact on the malaria parasite [17, 22].
Parasitoids (Hymenoptera, Insecta) are a group of parasitic insects. Some of them attach significant vectors of human disease, such as house flies, roaches and ticks, and some of them are extremely important regulators of agricultural pests . To study molecular mechanism of the innate immunity of parasitic wasps, the basic biology of their AMPs is needed to be established. As the first parasitic hymenopteran insect to have its genome sequenced , N. vitripennis provides a new resource for the identification of AMP genes in this class of insects in a genomic scale [25, 26]. Here, we report systemic characterization of the AMP repertoire from the N. vitripennis genome, which provide a perspective for the understanding of a possible relationship between immunity and parasitism. In comparison with A. mellifera and several non-Hymenopteran insects, N. vitripennis has developed a more complex antimicrobial immune system through genetic duplication and exon-shuffling.
Computational identification of AMPs in the N. Vitripennis genome
To identify N. vitripennis AMP repertoire, we employed three complementary approaches to database search of the N. vitripennis genome (see additional file 1). Firstly, we chose A. mellifera AMPs (e.g. mellitin, apidaecin, apisimin, abaecin, hymenoptaecin, defensins) as queries to recognize their orthologues in N. vitripennis by BLASTP and TBLASTN. Secondly, we performed pattern search by the Scanprosite program based on the cysteine arrangement pattern of CSαβ-type defensins . This strategy allowed identifying a large family of DLPs with a typical CXXXC and CXC motif and an amino-terminal signal peptide, which can be divided into three distinct subfamilies, including the known navidefensins  and two new subfamilies named nasonins and navitricins, respectively. Thirdly, we screened putative secreted peptides of < 150 amino acids from the N. vitripennis proteome by recognizing an amino-terminal signal sequence, from which AMP-like peptides rich in specific amino acids were identified. All AMPs found by the above approaches were again used as query to carry out BLASTP and TBLASTN iteratively until no new hits appeared. As a result, we identified a total of 44 AMPs in the genome of N. vitripennis, as shown in additional file 2, in which only nahymenoptaecin-1 and five defensins are recently reported [25, 26].
All these newly discovered peptides display remarkable AMP characteristics, as reflected by their small size and net positive charges at pH 7.0. Most peptides described here are smaller than 150 amino acids in length, except for nahymenoptaecin-1, nahymenoptaecin-2, nasonin-2 and nasonin-6, which all have undergone internal duplication (see additional file 2). N. vitripennis AMPs exhibit extensive sequence and structural diversity and can be distinguished into multiple molecular types including insect and fungal DLPs with the CSαβ fold; Pro- or Gly-rich abaecins and hymenoptaecins; horseshoe crab tachystatin-type AMPs with the ICK fold; and a linear α-helical peptide.
To confirm the reliability of our computational prediction, we undertook the N. vitripennis EST database search and found more than 60% peptides have corresponding transcripts, suggesting the genes encoding these peptides are expressed at the transcriptional level. For remaining 40% peptides whose cDNA sequences were not found in the EST database, one possibility is that their expression depends upon suitable microbial challenges  or is associated with different developmental stages of N. vitripennis, as observed in the D. melanogaster antifungal drosomycin , although pseudogene or wrong prediction resulted from computational methods cannot be completely excluded.
Classification of N. Vitripennis AMPs
Classical insect-type defensins (CITDs)
Defensins with the CSαβ structure are crucial effectors of innate immunity  which have been found in many insect species, the majority of them coming from different orders of the subclass Neoptera (Diptera, Coleoptera, Lepidopera, Hemiptera and Hymenoptera). These AMPs are primarily active against Gram-positive bacteria likely by forming voltage-dependent channels in bacterial membrane leading to a loss of cytoplasmic potassium [29, 30]. Their protective roles have been well documented by in vivo targeted disruption of defensin gene which resulted in the death of A. gambiae after Gram-positive bacterial infection . Insect defensins consist of 30-50 amino acids linked by three or four disulfide bridges. They share a common CSαβ structural motif with some functionally related scorpion neurotoxins targeting various ion channels, protease inhibitors, and even a plant sweet taste peptide [8, 32].
In A. mellifera, there are only two defensins named defensin1 and defensin2 . Amdefensin2 is a CITD with similar sequence, gene organization and molecular size to those from other insects. However, amdefensin1, also known as royalisin, has a lineage-specific carboxyl-terminal extension. Defensins with clear orthologous relationship with these two defensins have been identified in the N. vitripennis genome based on sequence similarity and identical precursor organization (signal peptide, propeptide and mature peptide) by our group . These two defensins have undergone gene duplication to form a multigene family containing five members (navidefensin1-1, 1-2, 2-1, 2-2 and 2-3). A detailed analysis of their exon-intron organization allowed identifying N. vitripennis-specific intron loss and gain (Figure 1). For example, navidefensin1-1 has a phase 2 intron at the residue lysine following the CXC motif, which is conserved with amdefensin1, while navidefensin1-2 lacks this intron. 3'-intron loss or gain has been frequently observed in a large number of genes by a possible reverse transcription-mediated mechanism . Whether dynamic evolution of introns is associated with the carboxyl-terminal extension in defensin1 remains an open question .
Three-dimensional (3D) models of navidefensin1-1 and navidefensin2-1 were constructed by comparative modeling with sapecin's structure (PDB: 1L4V)  as a template, which confirmed their typical CSαβ folding, as identified by an α-helix and a two-stranded β-sheet. Residues spanning S15-S23 in navidefensin1-1 and S18-A26 in navidefensin2-1 form an α-helix with separate hydrophobic (A16, A17 and L22 in navidefensin1-1 or A19, A21, V22, L25, and A26 in navidefensin2-1) and hydrophilic (S15, N20, and S23 in navidefensin1-1 or S18 and R23 in navidefensin2-1) surfaces. Two β-stands (residues G28 to C31 and C36 to R39 in navidefensin1-1 or G31 to C34 and C39 to R42 in navidefensin2-1) constitute an antiparallel sheet with a γ-core region which has been considered as a functional determinant responsible for antimicrobial activity of some defensins . Because of no suitable templates to build a model of the extended carboxyl-terminus of navidefensin1-1, alternatively an ab initio prediction method was employed. In both methods, a conserved CSαβ core was successfully predicted (Verify3D values > 0.2). For navidefensin1-1, the ab initio method predicted an α-helix in its carboxyl-terminus (residues K44 to R50) (Figure 2B). The existence of an extra glycine in the carboxyl-terminus of navidefensin1-1 suggests that Phe51 is likely amidated to provide a force for the helical stability, as predicted in bee defensin1 [33, 38].
In addition to the defensins mentioned above, we identified a new DLP subfamily without a propeptide, which contains fourteen members (named nasonin-1 to -14) (Figure 1). Sequence analysis revealed that nasonins are more closely related to two known non-classical defensins than CITDs in precursor organization, n-loop size and sequence similarity, one being antifungal termicin isolated from termites [39, 40]; another being defensin1 from the Mexican scorpion Centruroides limpidus limpidus. Intriguingly, such similarity was also observed between nasonins and a scorpion K+ channel toxin cobatoxin [42, 43] (Figure 1). These observations suggest that nasonins could have diverse functional features. For nasonin-2 and nasonin-6, their defensin unit has undergone 4 to 5 internal repeats, as described in fungal DLPs . What's the potential function of such repeats and how these repeats occur will be further discussed below.
We also predicted the structure of nasonin-1 by comparative modeling using scorpion neurotoxin cobatoxin (PDB: 1PJV)  as a template. Although having a shorter n-loop than CITDs, nasonins can also adopt a typical CSαβ architecture (Figure 2B). The reliability of this model was evaluated by Verify3D with a score of 0.248. Overall, nasonin-1 more resembles structurally termicin (PDB: 1MM0)  and cobatoxin (Figure 2B), but they differ in molecular surface charge distribution (data not shown), which could explain functional diversity of these peptides, as observed between human beta-defensin2 and the snake toxin crotamine .
Navidefensin3 is a unique non-classical defensin because it has identical precursor organization and gene structure to navidefensin1 and 2 but similar n-loop size and amino acid sequence to nasonins. In fact, several non-classical insect defensins with a short n-loop were also found in Lepidoptera insects (e.g. gallerimycin, spodoptericin and the B. mori defensins) [46–48]. In the tree presented in Figure 2A, navidefensin3 is clustered together with nasonins rather than navidefensins, suggesting that this DLP is an evolutionary link between these two multigene subfamilies of defensins. Its lack in A. mellifera is consistent with the absence of the nasonin subfamily in this species.
We also found two unique DLPs (named navitricin-1 and -2) in the N. vitripennis genome, which possess 8 cysteines with an identical alignment pattern to the family III members of fungal DLPs . Analysis of genomic sequences identified a phase 1 intron in the mature peptide-coding region of navitricins (Figure 3A).
As mentioned previously, navitricins have two additional cysteines which could form a fourth disulfide bridge. To confirm this, comparative modeling was again used to build the structure of navitricin-1 by choosing different templates, e.g. CITDs (sapecin and phormicin), the scorpion toxin BmTX1, and the plant sweet taste protein brazzein. In all these cases no acceptable models were generated and we thus tried to use the I-TASSER algorithm for structure simulation which produced a suitable model, as evaluated by Verify 3D (Figure 3B). Overall, this molecule is similar to the defensins described above including three identical disulfide bridges. Importantly, this model correctly predicts the fourth disulfide bridge linking the n- and c-loops together which makes the whole structure more compact. In fact, computational prediction was also applied to reveal an additional disulfide bridge in several insulin-like proteins .
(4) Biological activity
Previous studies have highlighted a crucial role of α-helical and γ-core regions of insect defensins in inhibiting the growth of Gram-positive bacteria, whereas their n-loops could mediate antibacterial specificity [37, 50, 51]. Interestingly, all nasonins identified here have a shorter n-loop relative to CITDs. To confirm whether these novel DLPs have different antibacterial spectrum from CITDs, we chose nasonin-1 as a representative. Reduced nasonin-1 was chemically synthesized and its in vitro refolded was performed by using reduced and oxidized glutathione (GSH/GSSG) under an alkaline condition (pH = 8.0). Refolded peptide was further purified by reverse-phase high performance liquid chromatography (RP-HPLC) (Figure 4A). In comparison with the reduced form, refolded nasonin-1 has a slightly earlier retention time (18 min vs. 19 min). The molecular weight (MW) of refolded nasonin-1 detected by matrix-assisted laser desorption/ionization time of flight mass spectrum (MALDI-TOF MS) is 3553.48 Da, about 7 Da smaller than the reduced form (3560.8 Da), in favor of the formation of three disulfide bridges (Figure 4B).
Subsequently, we evaluated antibacterial and antifungal activities of nasonin-1 using classical inhibition zone assays  and found that nasonin-1 only displayed activity against two Gram-negative bacteria at micromolar concentrations with a lethal concentrations (CL) of 1.52 μM for Stenotrophomonus sp. LZ-1 and of 15.4 μM for Escherichia coli ATCC 25922 (Figure 4C). No effect was observed on Gram-positive bacteria and fungi. This is the first example describing a DLP mainly active on Gram-negative bacteria. Recombinant navidefensin2-2 was recently reported to be active on Gram-positive bacteria . Whether such remarkable difference in target selectivity is associated with the change of n-loop length in these two Nasonia defensins needs further investigation.
AMPs rich in specific amino acids
AMPs rich in specific amino acids (e.g. glycine, proline, histidine, tryptophan or arginine) represent an additional class of AMPs acting on diverse microorganisms (e.g. bacteria, fungi, and virus) [53, 54]. The majority of these molecules are unstructured and they often inhibit microbial growth by entering into cells and interacting with proteins involved in key metabolic processes . Eighteen such AMPs have been identified here, including two known AMP families (abaecin and hymenoptaecin) [55, 56].
Abaecin is an inducible Pro-rich antibacterial peptide, firstly identified from the hemolymph of A. mellifera after bacterial challenge, and known to be regulated by the immune deficiency (IMD) pathway [55, 57]. They can act on both Gram-positive and Gram-negative bacteria with a moderate inhibitory concentration when compared to apidaecin which was also from A. mellifera. Subsequently, several abaecin-like genes were identified from the bee A. cerana and Bombus ignites, and the wasp Pteromalus puparum[58, 59]. Searching for the N. vitripennis genome identified three genes encoding protein precursors with clear sequence similarity to abaecins in their amino-terminal part of mature peptides. We named these abaecin-like genes nabaecin-1 to -3. Analysis of these precursor sequences allowed recognizing a dibasic cleavage site (RR) for putative post-translational processing by which an amino-terminal part (called N-terminal abaecin unit, NtAU) corresponding to the bee abaecin will be released (Figure 5). The NtAU of nabaecin-3 shares the highest sequence identity (50%) with A. mellifera abaecin. To understand how nabaecins extended their carboxyl-termini, we further searched the N. vitripennis database and found that a peptide (named navitripenicin) possesses highly identical amino acid sequence to the carboxyl-termini of nabaecins (herein named C-terminal navitripenicin unit, CtNU), but the region corresponding to abaecin is replaced by an acidic propeptide with a dibasic cleavage site. High content of glycines together with positively charged characteristics suggest that CtNU and navitripenicin belong to putative AMPs. A conserved phase 0 intron located at a nearly identical position among all these peptides (Figure 5) supports their evolutionary relatedness.
Hymenoptaecin is an antibacterial peptide rich in glycine, which was identified only in Hymenopteran insects [56, 59, 60]. The mature hymenoptaecin is released by cleaving off a signal peptide and an acidic propeptide after bacterial infection. Expression of the hymenoptaecin gene was also regulated by the IMD pathway . It inhibits the viability of Gram-negative and Gram-positive bacteria, including several human pathogens . Similarity search of the N. vitripennis genome identified two hymenoptaecin-like peptides (named nahymenoptaecin-1 and nahymenoptaecin-2) with a unique amino-terminal AMP-like region corresponding to the propeptide of bee hymenoptaecins in position (Figure 6A). Functional characterization and evolutionary significance of nahymenoptaecin-1 has been described recently .
Nahymenoptaecin-2 presents a more interesting structural feature in the mature peptide region, whose amino- and carboxyl-termini are respectively a Pro-rich peptide (named pronavicin) and a hymenoptaecin-like peptide with four internal repeats. Overall, these repeats can be well aligned with hymenoptaecin except a deletion of 31 amino acids in each repeat. Similarly, a conserved phase 0 intron is also present in the amino-termini of bee and parasitoid hymenoptaecins, suggesting that these peptides diverged from a common ancestor after speciation.
Although hymenoptaecins are believed to be only present in Hymenoptera, our study reveals that similar peptides were also evolved by Dipteran Drosophila (Figure 6A). These fly-derived peptides display a chimeric characteristic because they have an acidic propeptides, as bee hymenoptaecins, but their mature peptides more resemble the repeat of nahymenoptaecin-2 in size. Moreover, Drosophila hymenoptaecins have lost their introns, consistent with the lineage-specific intron loss in these species . All Drosophila hymenoptaecins clustering together in the tree (Figure 6B) indicates their monophyletic origin after separation of Diptera from Hymenoptera.
As a unique region of nahymenoptaecin-2, the amino-terminal pronavicin shares detectable sequence identity (25-41.2%) to other known Pro-rich AMPs from insects such as drosocin , apidaecin  and formaecin  (Figure 6C). Thus we synthesized this peptide to verify its antibacterial activity. Results showed that pronavicin was able to inhibit the growth of all Gram-positive and Gram-negative bacteria used here with differential potency. It displayed the strongest inhibitory effect on the Gram-positive bacterium Bacillus megaterium (C L = 3.11 μM) (Table 1). At high micromolar concentrations this peptide did not affect fungi and yeasts.
ICK fold is an evolutionarily conserved structural motif shared by a large group of polypeptides with diverse sequences and bioactivities. For example, venomous animals (e.g. spider and scorpion) usually employed this fold to develop their neurotoxins affecting channels' functions [64, 65] while horseshoe crabs, insects and some plants use ICK peptides as a part of their antibacterial defense system [66–69]. Peptides belonging to this fold type exhibit antibacterial (plant ICK peptides), antifungal (tachystatin and Alo-3), and chitin-binding (tachystatin-A) activities. Three peptides (named naickin-1, -2 and -3) with the ICK structure were characterized for the first time from the N. vitripennis genome that all share sequence similarity to an ICK peptide from A. mellifera (named amickin-1). All four peptides have a typical ICK cysteine arrangement pattern with a unique extended carboxyl-terminus rich in proline and positively charged residues. A phase 0 intron preceding the second cysteine is conserved between naickins and amickin-1 (Figure 7A), suggesting their homologous relationship. A typical cystine knot stabilized by three disulfide bridges can be observed in all these peptides, as revealed by comparative modeling using tachystatin-B1  as a template. Figure 7B presents the model structure of naickin-1 (residues T2 to T40) which is composed of an antiparallel β-sheet well fitted with tachystatin-B1.
Linear α-helical peptides
Linear peptides with an α-helical conformation are a large class of AMPs involved in immune response of arthropods, amphibians and vertebrates. Since cecropin, the first insect-derived linear AMP, was isolated from Cecropia moth , the number of such peptides dramatically increases in recent decades [72, 73]. Several linear α-helical AMPs have also been characterized in social wasps [74–76]. From the N. vitripennis genome, we identified one AMP-like peptide (called nahelixin) belonging to this class (Figure 7C). The precursor of nahelixin is 118 residues in length with an amino-terminal signal peptide and a putative carboxyl-terminal propeptide. The predicted mature peptide of 22 amino acids shares about 40% sequence identity with four known antibacterial peptides from frogs . Secondary structure prediction indicated it could adopt an α-helical structure (data not shown), which was further confirmed by ab initio structural prediction. An amphiphilic surface, with hydrophilic and hydrophobic residues separately arranged at two directions of the helical axis, is present in the model, as seen at the helical wheel projection (Figure 7C). Such structural feature represents a prerequisite for antibacterial activity of linear AMPs .
Inducible expression profiles of N. Vitripennis AMP-like genes
The up-regulation of AMP expression after infection is a common defense strategy by hosts to rapidly destroy invaders . To evaluate whether N. vitripennis AMPs have such expression feature, we chose twelve genes that belong to 4 different categories (i.e. nabaecins, nahymenoptaecins, navidefensins, and nasonins) as representatives. By using RT-PCR and DNA sequencing, we found that six genes are constitutively expressed because their cDNAs can be amplified from both infected and non-infected N. vitripennis adults. These genes are respectively nasonin-1, nasonin-3, nasonin-4, nabaecin-2, together with navidefensin1-1 and 1-2, two recently cloned genes from the non-infected N. vitripennis adults . Subsequently, we studied the expression profiles of all twelve genes after bacterial challenge by semi-quantitative RT-PCR. Results showed that only two genes (nasonin-3 and nasonin-4) whose expression was not induced by the bacteria. Four genes (navidefensin1-1, 1-2, nasonin-1 and nabaecin-2) weakly expressed in the non-infected wasp were significantly up-regulated by bacterial challenge. Transcripts of nabaecin-1, nahymenoptaecin-1 and nahymenoptaecin-2 were only detected in the infected wasp, indicating that their expression is initiated by infection (Figure 8). Differential expression patterns among different AMPs may represent a combinational strategy to form a defense network against diverse microbial pathogens. Although three navidefensin2 genes were not detected to be expressed in the non-infected and infected adults, their defense role in different developmental stages of Nasonia is possible .
AMPs in the parasitic wasp
With more and more genome sequences released, identifying new bioactive peptides/proteins by computational genomics approaches is becoming quite valuable. Some examples include hormones from nematode and mosquito, odorant binding-like proteins from honey bee, opossum immune genome, and reptile venom genes [78–82]. However, most strategies used are primarily based on sequence similarity that is useful in finding orthologues of known genes. An additional computational search strategy depends on the presence of conserved structural motifs within a peptide superfamily. Such examples include the identification of fungal DLPs and mammalian β-defensins [32, 83]. Here, we identified the N. vitripennis AMPs in a genomic scale by using an integrated strategy which combines similarity search, pattern recognition and AMP characteristics (see additional file 1).
Our previous and current studies indicate that at least 44 AMPs are present in N. vitripennis. In comparison with A. mellifera AMPs [18, 25, 26], we found an obvious complexity increase of immune peptides in the parasitic wasp, both in number and protein precursor size (Figure 9). For example, N. vitripennis has three abaecins, six defensins and fourteen nasonins whereas A. mellifera has only 1 abaecin, 2 defensins and lacks nasonins. Besides changes in number, N. vitripennis also undergoes wide sequence extension. For instance, members in the nabaecin family from N. vitripennis have an extended carboxyl-terminus of about 65 residues relative to A. mellifera while nahymenoptaecins underwent duplication in their amino- or carboxyl-terminus (Figures 5 and 6, also see ) which largely enlarges their molecular size. In addition, the same type AMPs (e.g. defensins) also considerably change in several aspects: 1) The length variation in the n-loop (i.e. nasonins); 2) The development of a wrapper disulfide bridge  (i.e. navitricins); 3) Carboxyl-terminal extension (i.e. navidefensin1s) (Figure 10); and 4) For nasonin-2 and -6, the nasonin-like domain (NLD) repeats 4-5 times and thus represents another type of complexity increase (Figure 11).
In comparison with several non-Hymenopteran insects whose genome sequences have been released (e.g. D. melanogaster and A. gambiae from Diptera, B. mori from Lepidoptera, and T. castaneum from Coleoptera) [17–21], N. vitripennis still have evolved more AMP genes both in kinds and numbers (Figure 9A). Based on the phylogenetic relationship of holometabolous insects , we elucidate several evolutionary events (e.g. gene expansion, terminal extension and tandem repeat) which might have occurred along the parasitoid branch (Figure 9B) to result in the lineage-specific complexity increase in N. vitripennis AMPs for synergetic defense in parasitic condition .
Evolution of the N. Vitripennis AMPs
Three different strategies have been recognized to shape the complex antimicrobial immune system in N. vitripennis:
Adjacent chromosome location together with gene structure conservation strongly supports tandem gene duplication occurring in at least three AMPs in the N. vitripennis genome, including navidefensins, nabaecins and nasonins (Figure 10A). An emerging trend in the study of AMPs is that their evolution is often driven by Darwin positive selection following gene duplication . However, we only detected weak positive selection signal in the amino-terminal abaecin unit (Table 2). No such signal was found in the gene families of defensins and nasonins, which is likely due to the sparse samplings of species used in this analysis because when more defensin sequences are included, accelerated amino acid substitutions were detected in the Nasonia lineage . An additional change after gene duplication is structural diversification of N. vitripennis AMPs (Figure 10B). These structural variations may result in evolutionary novelty by changing structural basis related to function. For example, n-loop has been proven to be functionally important in mosquito defensin  and its change in size will probably influence the activity and antibacterial spectra of these defensins, as observed in navidefensin2-2  and nasonin-1. There is also increasing awareness that alteration of the highly exposed wrapper disulfide bridge linkage pattern can lead to functional switch of scorpion sodium channel toxins via adjusting the conformation of key residues associated with toxin function . In this regard, navitricins with such a disulfide bridge are examples for further research. Finally, a typical amphipathic surface is present in the carboxyl-terminal α-helical region of the navidefensin1 subfamily members, suggesting a new antimicrobial unit has evolved on the classical defensin scaffold.
In N. vitripennis, we identified two unusual defensin-like genes (i.e. nasonin-2 and nasonin-6) that encode proteins carrying 4-5 repeats of NLD without protease processing signals (Figure 11). Although we do not know how these genes originated, the gene structure of nasonin-6 gives a glimpse into its history. As shown in Figure 11, each NLD is encoded by a single exon, which is evidence for exon duplication . From an evolutionary perspective, exon duplication mediated internal repeats could be viewed in the same light as gene duplication, which is usually considered to be the principal path for complexity increase of a protein. For example, the leucin-rich repeats (LRRs) in ribonuclease inhibitor (RI) are displayed in 1-30 tandem copies within a protein and such high LRR unit number could be critical for the evolution of both broad specificity and high affinity by RI . The presence of multiple defensin repeats in nasonin-2 and -6 could have some structural significance if we consider the functional importance of forming oligomerization for membrane permeabilization of β-defensins . Alternatively, new functional features were also evolved in some repeats, such as protease inhibitory activity  which may be useful in preventing from degradation of the whole molecule. The last domain of nasonin-2 has high sequence identity (86%) and a conserved phase 1 intron after signal peptide to nasonin-1, suggesting it originated from one ancestral copy of nasonin-1 by exon duplication and subsequent intron loss (Figure 1).
Exon-shuffling represents a major mechanism by which a limited set of modules can be shuffled by intronic recombination. Such mechanism has significantly contributed to the emergence of novel complex multidomain proteins, novel functions and increased organismic complexity. Gene structure conservation in multidomain immune molecules suggests that exon shuffling may be involved in the evolution of the immune system . However, in most case tracing such event is rather difficult due to accelerated sequence changes. Nabaecins, chimeric genes of abaecin and navitripenicin as mentioned previously (Figure 5), provides some clues to see a history of immune-related genes, in which exon-shuffling is likely involved. At the sequence level, the NtAU of nabaecin-3 is highly similar to A. mellifera abaecin (41.2% (protein) or 45.1% (DNA)), while its CtNU shares 41.4% (protein) and 56.2% (DNA) identity to navitripenicin. And the sequence identity extends to their 3' UTRs (Figure 12A). At the genomic level, all these three genes retain a conserved phase 0 intron at identical positions (Figure 5), indicating that they are evolutionarily related. In all these genes, CtNU and navitripenicin are encoded by a single exon. Therefore, we assume that after the wasp diverged from the bee lineage, an exon-shuffling, mediated by the conserved phase 0 intron, occurred between the ancestral ortholog of the A. mellifera abaecin and a duplicated navitripenicin gene, which resulted in the emergence of the first nabaecin gene. In subsequent evolution, a multigene family of three members formed by gene duplication (Figure 12B).
By using a combined computational and experimental strategy, we established for the first time the N. vitripennis peptidome associated with innate immunity. Three basic evolutionary scenarios are recognized in the generation of a complex antimicrobial system in N. vitripennis. Our work presented here will offer a basic platform for further studying the immunological and evolutionary significances of these newly discovered AMPs in parasitic insects. Whether these AMP-like genes with a complicate structure or much more copies are involved in parasitism remains an open question which constitutes our next research direction. In particular, our efficient approaches for finding new AMPs in a genomic scale can easily be applied to other model organisms including humans, which will accelerate our understandings on AMP-mediated innate immunity at a defense network level. Further functional characterization of these parasitic wasp-derived AMPs will also help explore new-type of drug leads for anti-infective therapy.
Strategies for gene discovery used here are provided in the supplementary information (see additional file 1). BLASTP and TBLASTN programs were used to characterize orthologues of known AMPs from A. mellifera against the database of N. vitripennis. When searching for hymenoptaecin in flies, BLAST was carried out in Flybase . The program WinGene  was employed to predict and translate a complete open reading frame from a selected nucleotide sequence. The N. vitripennis protein sequences downloaded from GenBank  and the FTP site  were applied to perform ScanProsite on the ExPaSy sever .
All sequences we identified as AMP-like peptides were submitted to SignalP 3.0 server  for signal peptide prediction. Peptide characteristics identified here were further analyzed on the ExPaSy server  for estimation of molecular weights and calculation of net charges at Protein Calculator . Secondary structure prediction was done on the ExPaSy server . To predict amphiphilic structure of peptides, helical wheels models were calculated on the server .
Sequence alignment and structural analysis
Amino acid sequences were aligned using Clustal X  with fine adjustment by hands and phylogenetic trees were constructed by MEGA 4.0 . Structures were modeled from ab initio by I-TASSER on-line [105–107] or on the Structural Bioinformatics server by comparative modeling method . Structure evaluation was carried out by Verify 3D [109, 110]. The MultiProt server was used to do structural superimposition and calculation of root mean square deviation (RMSD) . All model structures predicted here have been despoited in Protein Model DataBase  under id number of PM0075686-PM0075709.
Maximum likelihood analysis
Codon-substitution models were selected to estimate the nonsynonymous-to-synonymous rate ratio (ω = dN/dS) using the CODEML program of the PAML software package . Four models recommended by Yang make two likelihood ratio tests (LRTs) by M1a/M2a and M7/M8: M1a (nearly neutral model) constraints a proportion p0 of conserved sites with 0 < ω < 1, while a proportion p1 = 1-p0 of neutral sites with ω1 = 1; M2a (positive selection model) adds an extra class of sites with the proportion p2 = 1-p0 -p1 and with ω estimated from the data. M7 (β distribution model) does not allow for positively selected sites and M8(β and ω model) adds an extra class of sites to M7, allowing for ω > 1, which means the presence of positively selected sites. In addition, Model 0 (M0) assuming one ω for all sites that does not allow the existence of positive selection was chosen for negative control [114, 115]. Upon detection of the positively selected signals, the calculation of posterior probabilities was completed using the Bayes Empirical Bayes (BEB) method [115, 116].
Chemical synthesis and oxidative refolding
Pronavicin and reduced nasonin-1 were chemically synthesized by Xi'an Huachen Bio-Technology Co., Ltd. (Xi'an, China). The cyclization reaction to form disulfide bridges in the reduced nasonin-1 molecule was carried out in 0.1 M Tris-HCl (pH = 8.0) with 2 mM GSH and 0.2 mM GSSG. Refolded peptide was purified to be homogeneity by C18-RP-HPLC. Its MW was determined by MALDI-TOF MS on a Kratos PC Axima CFR plus (Shimazu Co. LTD, Kyoto).
Inhibition zone assays were carried out to evaluate antimicrobial activities of peptides [28, 52, 117]. Briefly, 50 μl bacteria or fungal spores with OD600 = 0.5 were inoculated to pre-heated 6 ml Luria-Bertaini's medium (1% bactotryptone, 0.5% bactoyeast extract and 0.5% NaCl) for bacteria or MEA medium (1% Malt Extract, 0.1% Peptone, 2% Glucose) for fungi containing 0.8% agar. The mixture was then spread on a 9-cm Petri dish, giving a depth of 1 mm. After settling, 2-mm wells were punched in the plate and then 2 μl of peptide under different concentrations was added to each well. Diameters of the inhibition zones were measured 16 hours after incubation at 30°C. C L was calculated by the Hultmark's method . The microorganisms used in this assay were listed in additional file 3.
Establishment of the infection model
N. vitripennis was maintained as a laboratory culture on pupae of the housefly Musca domestica. They were reared under a 14:10 light-dark cycle at 25°C in glass containers fed on 20% (v/v) honey solution. A standard method was applied to establish an infection model  of N. vitripennis. Briefly, a bacterial mixture of E. coli ATCC 25922 and M. luteus which had been separately incubated overnight was injected into abdomen of N. vitripennis adults with a micro-injector. About 50 wasps were challenged in total. Non-injected wasps were used as a control.
Preparation of total RNA
For isolating N. vitripennis total RNA, 50 wasps challenged for 16 hours by bacteria were grounded into fine powder in liquid nitrogen. The Trizol reagent (SBS Genetech, Beijing) was used to prepare total RNA according to the supplier's instructions. This method was also applied to prepare non-challenged total RNA .
Reverse transcription of total RNAs prepared from non-infected or infected N. vitripennis was performed using RT PreMix kit (TransGen, Beijing) and a universal oligo (dT) - containing adaptor primer dT3AP. The first-strand cDNA amplification in both non-infected and infected wasps was performed by using an AMP-like gene-specific primer and the universal primer 3AP according to the method described . PCR product was cloned into pGEM-T (TIANGEN, Beijing), and transformed into E. coli DH5α competent cells (TransGen, Beijing). Their sequences were confirmed by DNA sequencing. Primers used in this study are all listed in additional file 4.
For semi-quantitative RT-PCR, an AMP-like gene were amplified from non-infected or infected cDNAs by using a gene-specific primer and 3AP and PCR products obtained from different cycles (i.e. 25, 30, 35, 40, 45) were taken for comparison of the amounts of these products by electrophoresis on 1.5% agarose gels. N. vitripennis ribosome protein RP49 was chosen as an internal control which was amplified by using primers NvRp49-F2/3AP and the same cDNA templates.
Ezekowitz RA, Hoffmann JA: Innate immunity. 2003, Totowa, New Jersey: Human press
Hale JD, Hancock RE: Alternative mechanisms of action of cationic antimicrobial peptides on bacteria. Expert Rev Anti Infect Ther. 2007, 5: 951-959. 10.1586/14787220.127.116.111.
Brogden KA: Antimicrobial peptides: pore formers or metabolic inhibitors in bacteria?. Nat Rev Microbiol. 2005, 3: 238-250. 10.1038/nrmicro1098.
Jenssen H, Hamill P, Hancock RE: Peptide antimicrobial agents. Clin Microbiol Rev. 2006, 19: 491-511. 10.1128/CMR.00056-05.
Brown KL, Hancock RE: Cationic host defense (antimicrobial) peptides. Curr Opin Immunol. 2006, 18: 24-30. 10.1016/j.coi.2005.11.004.
Gao B, Sherman P, Luo L, Bowie J, Zhu S: Structural and functional characterization of two genetically related meucin peptides highlights evolutionary divergence and convergence in antimicrobial peptides. FASEB J. 2009, 23: 1230-1245. 10.1096/fj.08-122317.
Bulet P, Stocklin R: Insect antimicrobial peptides: structures, properties and gene regulation. Protein Pept Lett. 2005, 12: 3-11. 10.2174/0929866053406011.
Zhu S, Gao B, Tytgat J: Phylogenetic distribution, functional epitopes and evolution of the CSalphabeta superfamily. Cell Mol Life Sci. 2005, 62: 2257-2269. 10.1007/s00018-005-5200-6.
Gao B, Rodriguez M, Lanz-Mendoza H, Zhu S: AdDLP, a bacterial defensin-like peptide exhibits anti-Plasmodium activity. Biochem Biophys Res Commun. 2009, 387: 393-398. 10.1016/j.bbrc.2009.07.043.
Zhu S: Evidence for myxobacterial origin of eukaryotic defensins. Immunogenetics. 2007, 59: 949-954. 10.1007/s00251-007-0259-x.
Kim DH, Lee DG, Kim KL, Lee Y: Internalization of tenecin 3 by a fungal cellular process is essential for its fungicidal effect on Candida albicans. Eur J Biochem. 2001, 268: 4449-4458. 10.1046/j.1432-1327.2001.02364.x.
Casteels P, Ampe C, Jacobs F, Vaeck M, Tempst P: Apidaecins, antibacterial peptides from honeybees. EMBO J. 1989, 8: 2387-2391.
Oppenheim FG, Xu T, McMillian FM, Levitz SM, Diamond RD, Offner GD, Troxler RF: Histatins, a novel family of histidine-rich proteins in human parotid secretion. Isolation, characterization, primary structure, and fungistatic effects on Candida albicans. J Biol Chem. 1988, 263: 7472-7477.
Selsted ME, Novotny MJ, Morris WL, Tang YQ, Smith W, Cullor JS: Indolicidin, a novel bactericidal tridecapeptide amide from neutrophils. J Biol Chem. 1992, 267: 4292-4295.
Otvos L: The short proline-rich antibacterial peptide family. Cell Mol Life Sci. 2002, 59: 1138-1150. 10.1007/s00018-002-8493-8.
Lemaitre B, Hoffmann JA: The host defense of Drosophila melanogaster. Annu Rev Immunol. 2007, 25: 697-743. 10.1146/annurev.immunol.25.022106.141615.
Christophides GK, Zdobno EV, Barillas-Mury C, Birney E, Blandin S, Blass C, Brey PT, Collins FH, Danielli A, Dimopoulos G, Hetru C, Hoa NT, Hoffmann JA, Kanzok SM, Letunic I, Levashina EA, Loukeris TG, Lycett G, Meister S, Michel K, Moita LF, Müller HM, Osta MA, Paskewitz SM, Reichhart JM, Rzhetsky A, Troxler L, Vernick KD, Vlachou D, Volz J, von Mering C, Xu J, Zheng L, Bork P, Kafatos FC: Immunity-related genes and gene families in Anopheles gambiae. Science. 2002, 298: 159-165. 10.1126/science.1077136.
Evans JD, Aronstein K, Chen YP, Hetru C, Imler JL, Jiang H, Kanost M, Thompson GJ, Zou Z, Hultmark D: Immune pathways and defence mechanisms in honey bees Apis mellifera. Insect Mol Biol. 2006, 15: 645-656. 10.1111/j.1365-2583.2006.00682.x.
Sackton TB, Lazzaro BP, Schlenke TA, Evans JD, Hultmark D, Clark AG: Dynamic evolution of the innate immune system in Drosophila. Nat Genet. 2007, 39: 1461-1468. 10.1038/ng.2007.60.
Tanaka H, Ishibashi J, Fujita K, Nakajima Y, Sagisaka A, Tomimoto K, Suzuki N, Yoshiyama M, Kaneko Y, Iwasaki T, Sunaqawa T, Asaoka A, Mita K, Yamakawa M: A genome-wide analysis of genes and gene families involved in innate immunity of Bombyx mori. Insect Biochem Mol Biol. 2008, 38: 1087-1110. 10.1016/j.ibmb.2008.09.001.
Zou Z, Evans JD, Lu Z, Zhao P, Williams M, Sumathipala N, Hetru C, Hultmark D, Jiang H: Comparative genomic analysis of the Tribolium immune system. Genome Biol. 2007, 8: R177-10.1186/gb-2007-8-8-r177.
Waterhouse RM, Kriventseva EV, Meister S, Xi Z, Alvarez KS, Bartholomay LC, Barillas-Mury C, Bian G, Blandin S, Christensen BM, Dong Y, Jiang H, Kanost MR, Koutsos AC, Levashina EA, Li J, Ligoxygakis P, Maccallum RM, Mayhew GF, Mendes A, Michel K, Osta MA, Paskewitz S, Shin SW, Vlachou D, Wang L, Wei W, Zheng L, Zou Z, Severson DW, Raikhel AS, Kafatos FC, Dimopoulos G, Zdobnov EM, Christophides GK: Evolutionary dynamics of immune-related genes and pathways in disease-vector mosquitoes. Science. 2007, 316: 1738-1743. 10.1126/science.1139862.
Gadau J, Niehuis O, Peire A, Werren JH, Baudry E, Beukeboom L: The jewel wasp - Nasonia. Genome mapping and genomics in animals. Edited by: Hunter W, Kole C. 2008, Springer Berlin Heidelberg, 1: 27-41. full_text.
Werren JH, Richards S, Desjardins CA, Niehuis O, Gadau J, Colbourne JK, Beukeboom LW, Desplan C, Elsik CG, Grimmelikhuijzen CJP, Kitts P, Lynch J, Murphy T, Oliveira D, Smith CD, Zande van de L, Worley KC, Zdobnov EM: Functional and evolutionary insights from the genomes of three parasitoid Nasonia species. Science. 2010, 327: 343-348. 10.1126/science.1178028.
Gao B, Zhu S: Identification and characterization of the parasitic wasp Nasonia defensins: positive selection targeting the functional region?. Dev Comp Immunol. 2010,
Gao B, Zhu S: Characterization of a unique hymenoptaecin-like antimicrobial peptide in the parasitoid wasp Nasonia vitripennis. Process Biochem. 2009, 45: 139-146. 10.1016/j.procbio.2009.08.017.
Ferrandon D, Imler JL, Hetru C, Hoffmann JA: The Drosophila systemic immune response: sensing and signaling during bacterial and fungal infections. Nat Rev Immunol. 2008, 7: 862-874. 10.1038/nri2194.
Tian C, Gao B, Rodriguez Mdel C, Lanz-Mendoza H, Ma B, Zhu S: Gene expression, antiparasitic activity, and functional evolution of the drosomycin family. Mol Immunol. 2008, 45: 3909-3916. 10.1016/j.molimm.2008.06.025.
Dimarcq JL, Bulet P, Hetru C, Hoffmann J: Cysteine-rich antimicrobial peptides in invertebrates. Biopolymers. 1998, 47: 465-477. 10.1002/(SICI)1097-0282(1998)47:6<465::AID-BIP5>3.0.CO;2-#.
Takeuchi K, Takahashi H, Sugai M, Iwai H, Kohno T, Sekimizu K, Natori S, Shimada I: Channel-forming membrane permeabilization by an antibacterial protein, sapecin: determination of membrane-buried and oligomerization surfaces by NMR. J Biol Chem. 2004, 279: 4981-4987. 10.1074/jbc.M307815200.
Blandin S, Moita LF, Kocher T, Wilm M, Kafatos FC, Levashina EA: Reverse genetics in the mosquito Anopheles gambiae: targeted disruption of the Defensin gene. EMBO Rep. 2002, 3: 852-856. 10.1093/embo-reports/kvf180.
Zhu S: Discovery of six families of fungal defensin-like peptides provides insights into origin and evolution of the CSalphabeta defensins. Mol Immunol. 2008, 45: 828-838. 10.1016/j.molimm.2007.06.354.
Klaudiny J, Albert S, Bachanova K, Kopernicky J, Simuth J: Two structurally different defensin genes, one of them encoding a novel defensin isoform, are expressed in honeybee Apis mellifera. Insect Biochem Mol Biol. 2005, 35: 11-22. 10.1016/j.ibmb.2004.09.007.
Sverdlov AV, Babenko VN, Rogozin IB, Koonin EV: Preferential loss and gain of introns in 30 portions of genes suggests a reverse-transcription mechanism of intron insertion. Gene. 2004, 338: 85-91. 10.1016/j.gene.2004.05.027.
Zhu S: Did cathelicidins, a family of multifunctional host-defense peptides, arise from a cysteine protease inhibitor?. Trends Microbiol. 2008, 16: 353-360. 10.1016/j.tim.2008.05.007.
Hanzawa H, Shimada I, Kuzuhara T, Komano H, Kohda D, Inagaki F, Natori S, Arata Y: 1H nuclear magnetic resonance study of the solution conformation of an antibacterial protein, sapecin. FEBS Lett. 1990, 269: 413-420. 10.1016/0014-5793(90)81206-4.
Yount NY, Yeaman MR: Multidimensional signatures in antimicrobial peptides. Proc Natl Acad Sci. 2004, 101: 7363-7368. 10.1073/pnas.0401567101.
Casteels-Josson , Zhang KW, Capaci T, Casteels P, Tempst P: Acute transcriptional response of the honeybee peptide-antibiotics gene repertoire and required post-translational conversion of the precursor structures. J Biol Chem. 1994, 269: 28569-28575.
Bulmer MS, Crozier RH: Duplication and diversifying selection among termite antifungal peptides. Mol Biol Evol. 2004, 21: 2256-2264. 10.1093/molbev/msh236.
Lamberty M, Zachary D, Lanot R, Bordereau C, Robert A, Hoffmann JA, Bulet P: Insect immunity. Constitutive expression of a cysteine-rich antifungal and a linear antibacterial peptide in a termite insect. J Biol Chem. 2001, 276: 4085-4092. 10.1074/jbc.M002998200.
Rodriguez de la Vega RC, Garcia BI, D'Ambrosio C, Diego-Garcia E, Scaloni A, Possani LD: Antimicrobial peptide induction in the haemolymph of the Mexican scorpion Centruroides limpidus limpidus in response to septic injury. Cell Mol Life Sci. 2004, 61: 1507-1519. 10.1007/s00018-004-4096-x.
Jouirou B, Mosbah A, Visan V, Grissmer S, M'Barek S, Fajloun Z, Van Rietschoten J, Devaux C, Rochat H, Lippens G, El Ayeb M, De Waard M, Mabrouk K, Sabatier JM: Cobatoxin 1 from Centruroides noxius scorpion venom: chemical synthesis, three-dimensional structure in solution, pharmacology and docking on K+ channels. Biochem J. 2004, 377: 37-49. 10.1042/BJ20030977.
Selisko B, Garcia C, Becerril B, Gomez-Lagunas F, Garay C, Possani LD: Cobatoxins 1 and 2 from Centruroides noxius Hoffmann constitute a subfamily of potassium-channel-blocking scorpion toxins. Eur J Biochem. 1998, 254: 468-479. 10.1046/j.1432-1327.1998.2540468.x.
Da Silva P, Jouvensal L, Lamberty M, Bulet P, Caille A, Vovelle F: Solution structure of termicin, an antimicrobial peptide from the termite Pseudacanthotermes spiniger. Protein Sci. 2003, 12: 438-446. 10.1110/ps.0228303.
Yount NY, Kupferwasser D, Spisni A, Dutz SM, Ramjan ZH, Sharma S, Waring AJ, Yeaman MR: Selective reciprocity in antimicrobial activity versus cytotoxicity of hBD-2 and crotamine. Proc Natl Acad Sci. 2009, 106: 14972-14977. 10.1073/pnas.0904465106.
Kaneko Y, Tanaka H, Ishibashi J, Iwasaki T, Yamakawa M: Gene expression of a novel defensin antimicrobial peptide in the silkworm, Bombyx mori. Biosci Biotechnol Biochem. 2008, 72: 2353-2361. 10.1271/bbb.80263.
Volkoff AN, Rocher J, d'Alencon E, Bouton M, Landais I, Quesada-Moraga E, Vey A, Fournier P, Mita K, Devauchelle G: Characterization and transcriptional profiles of three Spodoptera frugiperda genes encoding cysteine-rich peptides. A new class of defensin-like genes from lepidopteran insects?. Gene. 2003, 319: 43-53. 10.1016/S0378-1119(03)00789-3.
Wen H, Lan X, Cheng T, He N, Shiomi K, Kajiura Z, Zhou Z, Xia Q, Xiang Z, Nakagaki M: Sequence structure and expression pattern of a novel anionic defensin-like gene from silkworm (Bombyx mori). Mol Biol Rep. 2009, 36: 711-716. 10.1007/s11033-008-9233-4.
Duret L, Guex N, Peitsch MC, Bairoch A: New insulin-like proteins with atypical disulfide bond pattern characterized in Caenorhabditis elegans by comparative sequence analysis and homology modeling. Genome Res. 1998, 8: 348-353.
Landon C, Barbault F, Leqrain M, Guenneuques M, Vovelle F: Rational design of peptides active against the gram positive bacteria Staphylococcus aureus. Proteins. 2008, 72: 229-239. 10.1002/prot.21912.
Yamada K, Natori S: Characterization of the antimicrobial peptide derived from sapecin B, an antibacterial protein of Sarcophaga peregrina (flesh fly). Biochem J. 1994, 298: 623-628.
Hultmark D: Quantification of antimicrobial activity, using the inhibition-zone assay. Techniques in Insect Immunology FL. Edited by: Dumphy AG, Marmaras VJ. 1998, Lavoisier Press, 103-107.
Iijima R, Kurata S, Natori S: Purification, characterization, and cDNA cloning of an antifungal protein from the hemolymph of Sarcophaga peregrina (flesh fly) larvae. J Biol Chem. 1993, 268: 12055-12061.
Robinson WE, McDougall B, Tran D, Selsted ME: Anti-HIV-1 activity of indolicidin, an antimicrobial peptide from neutrophils. J Leukoc Biol. 1998, 63: 94-100.
Casteels P, Ampe C, Riviere L, Van Damme J, Elicone C, Fleming M, Jacobs F, Tempst P: Isolation and characterization of abaecin, a major antibacterial response peptide in the honeybee (Apis mellifera). Eur J Biochem. 1990, 187: 381-386. 10.1111/j.1432-1033.1990.tb15315.x.
Casteels P, Ampe C, Jacobs F, Tempst P: Functional and chemical characterization of Hymenoptaecin, an antibacterial polypeptide that is infection-inducible in the honeybee (Apis mellifera). J Biol Chem. 1993, 268: 7044-7054.
Schluns H, Crozier RH: Relish regulates expression of antimicrobial peptide genes in the honeybee, Apis mellifera, shown by RNA interference. Insect Mol Biol. 2007, 16: 753-759.
Cheng X, Liu G, Ye G, Wang H, Shen X, Wu K, Xie J, Altosaar I: Screening and cloning of antimicrobial DNA sequences using a vital staining method. Gene. 2009, 430: 132-139. 10.1016/j.gene.2008.10.017.
Xu P, Shi M, Chen XX: Antimicrobial peptide evolution in the Asiatic honey bee Apis cerana. PLoS ONE. 2009, 4: e4239-10.1371/journal.pone.0004239.
Choi YS, Choo YM, Lee KS, Yoon HJ, Kim I, Je YH, Sohn HD, Jin BR: Cloning and expression profiling of four antibacterial peptide genes from the bumblebee Bombus ignitus. Comp Biochem Physiol B Biochem Mol Biol. 2008, 150: 141-146. 10.1016/j.cbpb.2008.02.007.
Coulombe-Huntington J, Majewski J: Intron loss and gain in Drosophila. Mol Biol Evol. 2007, 24: 2842-2850. 10.1093/molbev/msm235.
Bulet P, Dimarcq JL, Hetru C, Lagueux M, Charlet M, Hegy G, Van Dorsselaer A, Hoffmann JA: A novel inducible antibacterial peptide of Drosophila carries an O-glycosylated substitution. J Biol Chem. 1993, 268: 14893-14897.
Mackintosh JA, Veal DA, Beattie AJ, Gooley AA: Isolation from an ant Myrmecia gulosa of two inducible O-glycosylated proline-rich antibacterial peptides. J Biol Chem. 1998, 273: 6139-6143. 10.1074/jbc.273.11.6139.
Zhu S, Darbon H, Dyason K, Verdonck F, Tytgat J: Evolutionary origin of inhibitor cystine knot peptides. FASEB J. 2003, 17: 1765-1767. 10.1096/fj.02-0515com.
Billen B, Bosmans F, Tytgat J: Animal peptides targeting voltage-activated sodium channels. Curr Pharm Des. 2008, 14: 2492-2502. 10.2174/138161208785777423.
Barbault F, Landon C, Guenneugues M, Meyer JP, Schott V, Dimarcq JL, Vovelle F: Solution structure of Alo-3: a new knottin-type antifungal peptide from the insect Acrocinus longimanus. Biochemistry. 2003, 42: 14434-14442. 10.1021/bi035400o.
Fujitani N, Kawabata S, Osaki T, Kumiki Y, Demura M, Nitta K, Kawano K: Structure of the antimicrobial peptide tachystatin A. J Biol Chem. 2002, 277: 23654-23657. 10.1074/jbc.M111120200.
Liu Y, Luo J, Ren F, Peng C, Wu G, Zhao J: Purification, characterization, and molecular cloning of gene of a seed-specific antimicrobial protein from pokeweed. Plant Physiol. 2000, 122: 1015-1024. 10.1104/pp.122.4.1015.
Osaki T, Omotezako M, Nagayama R, Hirata M, Iwanaga S, Kasahara J, Hattori J, Ito I, Sugiyama H, Kawabata S: Horseshoe crab hemocyte-derived antimicrobial polypeptides, tachystatins, with sequence similarity to spider neurotoxins. J Biol Chem. 1999, 274: 26172-26178. 10.1074/jbc.274.37.26172.
Fujitani N, Kouno T, Nakahara T, Takaya K, Osaki T, Kawabata S, Mizuguchi M, Aizawa T, Demura M, Nishimura S, Kawano K: The solution structure of horseshoe crab antimicrobial peptide tachystatin B with an inhibitory cystine-knot motif. J Pept Sci. 2007, 13: 269-279. 10.1002/psc.846.
Hultmark D, Steiner H, Rasmuson T, Boman HG: Insect immunity. Purification and properties of three inducible bactericidal proteins from hemolymph of immunized pupae of Hyalophora cecropia. Eur J Biochem. 1980, 106: 7-16.
Zasloff M: Magainins, a class of antimicrobial peptides from Xenopus skin: isolation, characterization of two active forms, and partial cDNA sequence of a precursor. Proc Natl Acad Sci. 1987, 84: 5449-5453. 10.1073/pnas.84.15.5449.
Li J, Xu X, Xu C, Zhou W, Zhang K, Yu H, Zhang Y, Zheng Y, Rees HH, Lai R, Yang D, Wu J: Anti-infection peptidomics of amphibian skin. Mol Cell Proteomics. 2007, 6: 882-894. 10.1074/mcp.M600327-MCP200.
Ifrah D, Doisy X, Ryge TS, Hansen PR: Structure-activity relationship study of anoplin. J Pept Sci. 2005, 11: 113-121. 10.1002/psc.598.
Konno K, Rangel M, Oliveira JS, Dos Santos Cabrera MP, Fontana R, Hirata IY, Hide I, Nakata Y, Mori K, Kawano M, Fuchino H, Sekita S, Neto JR: Decoralin, a novel linear cationic alpha-helical peptide from the venom of the solitary eumenine wasp Oreumenes decoratus. Peptides. 2007, 28: 2320-2327. 10.1016/j.peptides.2007.09.017.
Xu X, Li J, Lu Q, Yang H, Zhang Y, Lai R: Two families of antimicrobial peptides from wasp (Vespa magnifica) venom. Toxicon. 2006, 47: 249-253. 10.1016/j.toxicon.2005.10.015.
Olson L, Soto AM, Knoop FC, Conlon JM: Pseudin-2: an antimicrobial peptide with low hemolytic activity from the skin of the paradoxical frog. Biochem Biophys Res Commun. 2001, 288: 1001-1005. 10.1006/bbrc.2001.5884.
Belov K, Sanderson CE, Deakin JE, Wong ES, Assange D, McColl KA, Gout A, de Bono B, Barrow AD, Speed TP, Trowsdale J, Papenfuss AT: Characterization of the opossum immune genome provides insights into the evolution of the mammalian immune system. Genome Res. 2007, 17: 982-991. 10.1101/gr.6121807.
Forêt S, Maleszka R: Function and evolution of a gene family encoding odorant binding-like proteins in a social insect, the honey bee (Apis mellifera). Genome Res. 2006, 16: 1404-1413. 10.1101/gr.5075706.
Nathoo AN, Moeller RA, Westlund BA, Hart AC: Identification of neuropeptide-like protein gene families in Caenorhabditis elegans. Proc Natl Acad Sci. 2001, 98: 14000-14005. 10.1073/pnas.241231298.
Riehle MA, Garczynski SF, Crim JW, Hill CA, Brown MR: Neuropeptides and peptide hormones in Anopheles gambiae. Science. 2002, 298: 172-175. 10.1126/science.1076827.
Whittington CM, Papenfuss AT, Bansal P, Torres AM, Wong ES, Deakin JE, Graves T, Alsop A, Schatzkamer K, Kremitzki C, Ponting CP, Temple-Smith P, Warren WC, Kuchel PW, Belov K: Defensins and the convergent evolution of platypus and reptile venom genes. Genome Res. 2008, 18: 986-994. 10.1101/gr.7149808.
Schutte BC, Mitros JP, Bartlett JA, Walters JD, Jia HP, Welsh MJ, Casavant TL, McCray PB: Discovery of five conserved beta -defensin gene clusters using a computational search strategy. Proc Natl Acad Sci. 2002, 99: 2129-2133. 10.1073/pnas.042692699.
Zhu S, Gao B: Molecular characterization of a new scorpion venom lipolysis activating peptide: Evidence for disulfide bridge-mediated functional switch of peptides. FEBS Lett. 2006, 580: 6825-6836. 10.1016/j.febslet.2006.11.040.
Savard J, Tautz D, Richard S, Weinstock GM, Gibbs RA, Werren JH, Tettelin H, Lercher MJ: Phylogenomic analysis reveals bees and wasps (Hymenoptera) at the base of the radiation of Holometabolous insects. Genome Res. 2006, 16: 1334-1338. 10.1101/gr.5204306.
Khush RS, Lemaitre B: Genes that fight infection. What the Drosophila genome says about animal immunity. Trends Genet. 2000, 16: 442-449. 10.1016/S0168-9525(00)02095-3.
Tennessen JA: Molecular evolution of animal antimicrobial peptides: widespread moderate positive selection. J Evol Biol. 2005, 18: 1387-1394. 10.1111/j.1420-9101.2005.00925.x.
Street TO, Rose GD, Barrick D: The role of introns in repeat protein gene formation. J Mol Biol. 2006, 360: 258-266. 10.1016/j.jmb.2006.05.024.
Haigis MC, Haag ES, Raines RT: Evolution of ribonuclease inhibitor by exon duplication. Mol Biol Evol. 2002, 19: 959-963.
Hoover DM, Rajashankar KR, Blumenthal R, Puri A, Oppenheim JJ, Chertov O, Lubkowski J: The structure of human beta-defensin-2 shows evidence of higher order oligomerization. J Biol Chem. 2000, 275: 32911-32918. 10.1074/jbc.M006098200.
Carvalho AO, Gomes VM: Plant defensins-Prospects for the biological functions and biotechnological properties. Peptides. 2009, 30: 1007-1020. 10.1016/j.peptides.2009.01.018.
Betts MJ, Guiqó R, Aqarwal P, Russell RB: Exon structure conservation despite low sequence similarity: a relic of dramatic events in evolution?. EMBO J. 2001, 20: 5354-5360. 10.1093/emboj/20.19.5354.
BLAST program. [http://blast.ncbi.nlm.nih.gov/]
Hennig L: WinGene/WinPep: user-friendly software for the analysis of amino acid sequences. Biotechniques. 1999, 26: 1170-1172.
GenBank Nasonia vitripennis Genome FTP site. [ftp://ftp.ncbi.nih.gov/genomes/Nasonia_vitripennis/]
SignalP 3.0 Server. [http://www.cbs.dtu.dk/services/SignalP]
Protein Calculator. [http://www.scripps.edu/~cdputnam/protcalc.html]
Analysis of individual sequence. [http://www.bioinfo.org.cn/lectures/index-7.html#ex2]
Clustal X. [ftp://ftp-igbmc.u-strasbg.fr/pub/]
MEGA 4.0 Server. [http://www.megasoftware.net]
Wu S, Skolnick J, Zhang Y: Ab initio modeling of small proteins by iterative TASSER simulations. BMC Biol. 2007, 5: 17-10.1186/1741-7007-5-17.
Zhang Y: Template-based modeling and free modeling by I-TASSER in CASP7. Proteins. 2007, 69 (Suppl 8): 108-117. 10.1002/prot.21702.
Zhang Y: I-TASSER server for protein 3D structure prediction. BMC Bioinformatics. 2008, 9: 40-10.1186/1471-2105-9-40.
Douguet D, Labesse G: Easier threading through web-based comparisons and cross-validations. Bioinformatics. 2001, 17: 752-753. 10.1093/bioinformatics/17.8.752.
Bowie JU, Lüthy R, Eisenberg D: A method to identify protein sequences that fold into a known three-dimensional structure. Science. 1991, 253: 164-170. 10.1126/science.1853201.
Lüthy R, Bowie JU, Eisenberg D: Assessment of protein models with three-dimensional profiles. Nature. 1992, 356: 83-85. 10.1038/356083a0.
Shatsky M, Nussinov R, Wolfson HJ: A method for simultaneous alignment of multiple protein structures. Proteins. 2004, 56: 143-156. 10.1002/prot.10628.
Protein Model database. [http://mi.caspur.it/PMDB/]
Yang Z, Nielsen R, Goldman N, Pedersen AM: Codon-substitution models for heterogeneous selection pressure at amino acid sites. Genetics. 2000, 155: 431-449.
Yang Z, Wong WS, Nielsen R: Bayes empirical Bayes inference of amino acid sites under positive selection. Mol Biol Evol. 2005, 22: 1107-1118. 10.1093/molbev/msi097.
Zhang J, Nielsen R, Yang Z: Evaluation of an improved branch-site likelihood method for detecting positive selection at the molecular level. Mol Biol Evol. 2005, 22: 2472-2479. 10.1093/molbev/msi237.
Gao B, Zhu S: Differential potency of drosomycin to Neurospora crassa and its mutant: implications for evolutionary relationship between defensins from insects and plants. Insect Mol Biol. 2008, 17: 405-411. 10.1111/j.1365-2583.2008.00810.x.
Tzou P, Meister M, Lemaitre B: Methods for studying infection and immunity in Drosophila. Methods in microbiology. Edited by: Sansonetti PJ, Zychlinsky A. 2002, Academic Press, 507-530. full_text.
Gao B, Tian C, Zhu S: Inducible antibacterial response of scorpion venom gland. Peptides. 2007, 28: 2299-2305. 10.1016/j.peptides.2007.10.004.
We thank Human Genome Sequencing Center, Baylor College of Medicine http://www.hgsc.bcm.tmc.edu/ for providing N. vitripennis sequence data for public use and Prof. Chuanling Qiao (Institute of Zoology, Chinese Academy of Sciences, China), Prof. Fengyan Bai (Institute of Microbiology, Chinese Academy of Sciences, China) and Dr. Jianguo Zhou (Oak Ridge National Laboratory, USA) for providing microbial strains used in this study. We also thank Prof. Fengqin He (Institute of Zoology, Chinese Academy of Sciences, China) for providing M. domestica. This work was supported by grants from the National Natural Science Foundation of China (30730015, 30621003 and 90608009) and the 973 Program from the Ministry of Science and Technology of China (2010CB945304).
The authors declare that they have no competing interests.
SZ designed the experiments. CT, SZ and GY wrote the paper. CT conducted database search. CT, BG and QF performed the experiments. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Strategy of database searches of putative Nasonia vitripennis antimicrobial peptides. A diagram for the computational identification of putative AMPs. (DOC 64 KB)
Additional file 2: Summary of N. vitripennis antimicrobial peptides. A table for the summary of N. vitripennis antimicrobial peptides. (DOC 126 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Tian, C., Gao, B., Fang, Q. et al. Antimicrobial peptide-like genes in Nasonia vitripennis: a genomic perspective. BMC Genomics 11, 187 (2010). https://doi.org/10.1186/1471-2164-11-187