- Research article
- Open Access
A comprehensive transcriptome and immune-gene repertoire of the lepidopteran model host Galleria mellonella
© Vogel et al; licensee BioMed Central Ltd. 2011
Received: 25 February 2011
Accepted: 11 June 2011
Published: 11 June 2011
The larvae of the greater wax moth Galleria mellonella are increasingly used (i) as mini-hosts to study pathogenesis and virulence factors of prominent bacterial and fungal human pathogens, (ii) as a whole-animal high throughput infection system for testing pathogen mutant libraries, and (iii) as a reliable host model to evaluate the efficacy of antibiotics against human pathogens. In order to compensate for the lack of genomic information in Galleria, we subjected the transcriptome of different developmental stages and immune-challenged larvae to next generation sequencing.
We performed a Galleria transcriptome characterization on the Roche 454-FLX platform combined with traditional Sanger sequencing to obtain a comprehensive transcriptome. To maximize sequence diversity, we pooled RNA extracted from different developmental stages, larval tissues including hemocytes, and from immune-challenged larvae and normalized the cDNA pool. We generated a total of 789,105 pyrosequencing and 12,032 high-quality Sanger EST sequences which clustered into 18,690 contigs with an average length of 1,132 bases. Approximately 40% of the ESTs were significantly similar (E ≤ e-03) to proteins of other insects, of which 45% have a reported function. We identified a large number of genes encoding proteins with established functions in immunity related sensing of microbial signatures and signaling, as well as effector molecules such as antimicrobial peptides and inhibitors of microbial proteinases. In addition, we found genes known as mediators of melanization or contributing to stress responses. Using the transcriptomic data, we identified hemolymph peptides and proteins induced upon immune challenge by 2D-gelelectrophoresis combined with mass spectrometric analysis.
Here, we have developed extensive transcriptomic resources for Galleria. The data obtained is rich in gene transcripts related to immunity, expanding remarkably our knowledge about immune and stress-inducible genes in Galleria and providing the complete sequences of genes whose primary structure have only partially been characterized using proteomic methods. The generated data provide for the first time access to the genetic architecture of immunity in this model host, allowing us to elucidate the molecular mechanisms underlying pathogen and parasite response and detailed analyses of both its immune responses against human pathogens, and its coevolution with entomopathogens.
The introduction of novel high through-put sequencing technologies provides insight into the genetic architecture of an increasing number of non-model organisms including insects. Next-generation (NextGen) pyrosequencing has become an important tool in transcriptomic studies and allows targeted identification of genes which are (differentially) expressed in distinct tissues or cells, during development, or upon activation of immune responses. This technology has been used, for example, to characterize both the midgut-specific and the immunity-related transcriptome of Manduca sexta, which has emerged as a model in lepidopteran biochemistry and physiology [1, 2]. In this study, we subjected the immunity-related transcriptome of the greater wax moth Galleria mellonella to a combination of Sanger and NextGen sequence analysis. Our study was motivated by two reasons. Firstly, Galleria is suited to identify ancient features of innate immunity in lepidopterans because it belongs to the family Pyralidae which has been placed in a basal phylogenetic position within the Lepidoptera. Secondly, Galleria represents a powerful, reliable and proven model system for innate immunity studies. It is currently used as a host system to reconstruct rapid reciprocal adaptations during host-parasite coevolution  and as a an alternative model host for testing human pathogens, which is ethically better acceptable than mammalian hosts such as mice, rats and rabbits [4, 5]. Galleria caterpillars prosper world-wide in use as alternative mini-hosts because they combine advantages shared with other invertebrate host models with benefits that are unique to this lepidopteran. The advantages of the nematode Caenorhabditis elegans and the fruit fly Drosophila melanogaster are complete, well-annotated genomes and that microarrays, RNA interference libraries and mutant strains are available which allow analysis of host-pathogen interactions at the molecular level . However, the larger size of Galleria caterpillars enables precise injection of antibiotics or a number of pathogens, easy manipulation and collection of tissue and hemolymph samples to study pathophysiology with, for example, proteomic approaches.
Further advantages of Galleria are (i) the low overall costs of breeding large numbers, providing an inexpensive whole-animal high throughput infection assay system , (ii) their worldwide commercial availability, e.g., they are sold as bait for fishermen or as food for pets (reptiles), (iii) the positive correlation between the pathogenicity of bacteria and fungi when evaluated in Galleria and mice , (iiii) and that this heterologous insect host can be adapted in the laboratory to human physiological temperature (37°C). This is essential in order to mimic the physiological conditions in mammals because human pathogens are adapted to the physiological temperature of their host which is often required for the synthesis and the release of their pathogenic or virulence factors [4, 5]. These advantages have convinced an increasing number of researchers to favor Galleria as a mini-host model for prominent pathogenic bacteria and fungi that are responsible for severe human diseases such as Bacillus cereus , Enterococcus faecalis , Listeria monocytogenes , Pseudomonas aeruginosa , Staphylococcus aureus , Candida albicans  and Cryptococcus neoformans . In addition, a number of antimicrobial peptides and inhibitors of microbial virulence factors have been discovered during the past decade in Galleria whose therapeutic potential in medicine and plant protection is presently being explored [16, 17].
The major disadvantage of Galleria as a heterologous host system is that neither genome nor transcriptome sequence data are available and, therefore, important information about the immunity and stress related genes and their expression are lacking. Consequently, this study was designed to fill this gap and to provide a data set which enables more detailed studies, for example microarray or proteomic analysis, in the future. In order to induce expression of immunity-related genes in this lepidopteran species we injected a bacterial lipopolysaccharide (LPS) preparation into last instar larvae which has been proven as a potent elicitor of immune responses in Galleria  and other insect species [19, 20]. Normalized larval dscDNA was sequenced using Roche 454 FLX and Sanger (directional long reads) methods. The combining of both technologies provided deep sequencing coverage of the expressed genes relevant to this research project.
Because of the large hemolymph sample volumes that can be obtained from Galleria caterpillars, their host response to pathogens can easily be studied at the peptide and protein level . To test the correlation between transcriptomic and proteomic data we collected hemolymph samples from untreated and LPS-injected larvae. In order to identify peptides and proteins that are secreted within the hemolymph upon activation of innate immune responses, we used 2D-gelelectrophoresis combined with mass spectrometric analysis of spots that appear or are enhanced upon injection of LPS. Complementary proteomic analysis of hemolymph samples confirmed induced expression and release into the hemolymph of proteins known to mediate recognition of microbes, immunity-related signaling or killing of microbes.
Results and Discussion
Transcriptome assembly and functional analyses using Gene Ontologies
Summary statistics for Galleria mellonella expressed sequence tag (EST) analysis
454 FLX technology
Sanger + 454 FLX
Total number of reads
Average length of read (bases)
Total number of reads in the assembly
Total number of contigs
Total number of singlets
Number of unique sequences
Average/Largest contig size
The patterns of GO category associations nonetheless differed between these two insect species in a few categories, with relatively high abundance of Multicellular organismal development, Anatomical structure morphogenesis and Cellular developmental process in Bombyx and Biosynthetic process and Macromolecule metabolic process and response to stress being more abundant in Galleria. In addition to this, several categories were only present in Galleria (Oxidation reduction and Cell death). Differences in GO category associations between Bombyx and Galleria might be attributed to the fact that ESTs of the latter originate predominantly from larvae.
Recognition of pathogen or damage associated molecular pattern genes
Gram negative bacteria binding proteins (GNBPs) and β-1,3-glucan recognition proteins (βGRPs) have been extensively studied as pattern recognition proteins in Lepidoptera [26–28]. Most of these proteins are produced in the fat body and secreted into the caterpillar's hemolymph. Some are constitutively present whereas others are induced upon microbial infection. We have identified five different ßGRPs in the Galleria EST data collection, including one most similar to the midgut-specific short ßGRP with glucanase activity as previously described . To further examine the relationships among βGRP proteins across insects and the ßGRPs identified in Galleria, a total of 45 sequences from 24 species, including many proteins that had previously been found in insect hemolymph, were collected and used to construct a Bayesian phylogeny (Additional File 3). The phylogenetic analysis revealed that these sequences clustered in two distinct clades. One of these clades is clearly separated from the other clades by a high posterior probability and contains the Helicoverpa armigera Glucanase-1 protein (described in Pauchet et al. ), and sequences from cDNA libraries made from midgut tissue of different Lepidoptera species, including one Galleria sequence. This phylogeny suggests that Galleria does have all of the ßGRPs found in more derived Lepidopteran species, including the gene coding for a protein with glucanase activity. This supports the idea of an ancient ßGRP duplication event in Lepidoptera, leading to paralogues that have different functions.
Immunity related signaling
In insects, cell signaling against fungal and bacterial pathogens occurs through the Toll, Imd, and Jak-STAT pathways . These pathways are quite similar to the vertebrate (e.g. TNF) signaling pathways, and induce the expression of antimicrobial peptides and other molecules through interaction with NFkB factors. The major signaling pathways Toll and Imd are represented by central receptors such as toll, toll-like, 18 wheeler and related LRR repeat-containing G-protein coupled receptors. We have identified at least three different toll or toll-like receptor transcripts in the Galleria dataset. The exact number of different toll receptors is not easy to evaluate, as some of the transcripts are incomplete and the predicted amino acid sequences do not always overlap. In addition to toll, we identified two different 18 wheeler partial transcripts with homology to Spodoptera frugiperda 18 wheeler (Genbank entry ADV41489: Li, S: A 18 wheeler toll receptor gene from S. frugiperda cell is in response to LPS and Saccharomyces cerevisiae stimulation). However, a critical evaluation of the role of 18 wheeler in Drosophila has put its postulated function as a pattern recognition receptor for Gram negative bacteria into question . Interestingly, we identified transcripts encoding for the transcription factors NFkB and relish which function as obligate dimmers. Relish regulates downstream of the IMD pathway expression of antimicrobial peptides in Drosophila . While the signaling pathways that stimulate immune gene expression have been well characterized by genetic analysis in Drosophila, they are far from being well understood in most other insect species. However, several proteins involved in these pathways have recently been characterized in Lepidoptera. One such pathway involves proteolytic activation of a cytokine called Spaetzle, which functions in dorsal-ventral patterning during early embryonic development and in the antimicrobial immune response in larvae and adults. Most interestingly, it could be shown that injection of Spaetzle into M. sexta larvae stimulated expression of several immune-related peptides and proteins, including cecropin, attacin, moricin and lysozyme . We have identified a Spaetzle homolog in Galleria. The Gme-Spaetzle cDNA encodes a polypeptide with 29%, 42% and 44% identity to N. vitripennis, B. mori and M. sexta, respectively (Additional file 4A-alignment of Galleria, Bombyx, Manduca Spaetzle).
In addition to major immune signaling proteins, we identified a calreticulin sequence in the immune-induced Galleria transcriptome data. Calreticulin is involved in signal transduction events associated with innate immunity, cell adhesion, angiogenesis and apoptosis in mammals. The level of calreticulin on the surface of human dentritic cells and polymorphonuclear phagocytes correlates with their phagocytotic ability . Induction of calreticulin upon LPS challenge has recently been determined in other invertebrates such as the planarian Schmidtea mediterranea which is suggestive for its evolutionarily conserved roles in innate immunity . LPS-challenge also induced expression of tetraspanins whose role in modulating immune signal complexes in vertebrates is well established . Its induced expression upon LPS-injection has also been documented in ancient insects such as the firebrat . Similarly, an ankyrin repeat domain containing protein was found both in this basal insect and in Galleria.
Antimicrobial peptides and proteins
Cecropins represent another group of linear and amphipathic peptides with a-helical structure. The first member of this peptide family exhibiting antibacterial and antifungal activity was discovered in and isolated from the hemolymph of the silk moth Hyalophora cecropia and has therefore been named cecropin . The cecropin-like peptide from Galleria is synthesized as a propeptide, with a putative 22-residue signal peptide, a 4-residue propeptide and a 39-residue mature peptide with a mass of 4.3 kDa. Like cecropins from other insects it exhibits potent activity against both Gram-positive and Gram-negative bacteria . We have identified four different cecropins in the Galleria transcriptome dataset, including a more diverged D-type cecropin. This surprisingly large number of different cecropins (Additional file 4B) covers a larger fraction of the amino acid diversity encountered when comparing cecropins from across the Lepidoptera.
We determined both cysteine-rich peptides reported from Galleria which exclusively inhibit growth of filamentous fungi, the defensin-like antifungal peptides galiomicin  and gallerimycin . At least the latter contributes to innate immune responses mediating resistance of G. mellonella larvae against normally lethal infection by the human pathogenic yeast C. albicans . Transgenic expression of gallerimycin has been shown to confer resistance to fungal diseases on crops . A homologue of spodoptericin, the third defensin-like peptide discovered in Lepidoptera , is also present in our Galleria transcriptome.
In a previous study, we used the suppression subtractive hybridization method to screen for genes that are induced in Galleria upon challenge with LPS . This approach resulted in the discovery of novel peptides and protein families which were also found in this extended transcriptomic study. For example, we discovered a cobatoxin-like molecule and a protein which was named Gall-6-tox due to its six conserved tandem repeats of cysteine-stabilized alpha beta motifs (CS-αβ), the structural scaffold characteristic of invertebrate defensins and scorpion toxins. Homologues of Gal-6-tox differing in the number of tandem repeats of the CS-αβ motif were later found in other lepidopterans such as Bombyx mori and Spodoptera exigua. It turned out that they belong to a novel family of atypical defensin-derived immune-related proteins, which is specific to Lepidoptera and which is now called X-tox . Moreover, our study confirmed the induced expression of tenascin-like proteins in Galleria upon LPS-challenge , which represent immune effector molecules known from vertebrates. However, using RACE-PCR we obtained the full-length cDNA which is considerably shorter than vertebrate tenascins and lacks characteristic tenascin domains such as fibronectin type-3-like repeats. These findings make the relation of the identified sequences to tenascins unlikely.
In addition, we identified a full-length cDNA sequence which is identical to the deduced protein sequence of a Galleria proline-rich peptide  and almost identical to two protein fragments identified in a previous study analyzing hemolymph peptide fragments in Galleria . Finally, our transcriptomic analysis confirmed the presence of genes encoding cobatoxin-like peptides [18, 45] (Additional file 4C).
Inhibitors of microbial proteases
Several induced transcripts encode for transferrin which represents a multifunctional and evolutionarily conserved player in innate immunity. Its role in binding and removing available free iron ions, thus creating unfavorable environments for bacteria has first been reported in vertebrates . A recent study using B. mori confirmed both the induced expression of transferrin upon LPS-challenge and its contribution to antibacterial iron-withholding strategy in Lepidoptera B. mori .
Stress response genes
In line with our previous studies in which we used LPS-challenge to screen for inducible immunity-related genes in insects and other invertebrates [18–20, 29, 61] we determined induced expression of genes involved in detoxification and stress adaptation such as apolipoprotein D, cytochrome P450s, gluthathione S-transferases, and a number of heat shock proteins which further supports our hypothesis that interdependencies between immune and stress responses are evolutionarily conserved in insects [18–20, 29, 61]. Glutathione S-transferases (GSTs) are a large and diverse family of detoxification enzymes found in most organisms. GSTs help to protect cells from oxidative stress, but they also play a central role in the detoxification of both endogenous and xenobiotic compounds (e.g. plant secondary metabolites or insecticides) and are involved in intracellular transport and biosynthesis of hormones. Eukaryotes contain multiple GSTs belonging to different GST classes and with differing enzyme activities to accommodate the wide range of functions of this enzyme family. The insect GST supergene family encodes a group of proteins that have been assigned to at least six classes: Delta, Epsilon, Omega, Sigma, Theta and Zeta [62, 63]. The Delta and Epsilon classes, both specific to insects, are the largest classes and are often involved in xenobiotic metabolism whereas the Omega, Sigma, Theta and Zeta classes have a much wider taxonomic distribution and likely play essential housekeeping roles [62, 63]. Herbivorous insects have to cope with toxic plant metabolites taken up with their diet and GSTs can play an important role in their detoxification [64–66].
We identified a total of 19 different GSTs in Galleria larval ESTs out of which 2 were microsomal GSTs. Five out of the six classes identified in other insect species are represented and most of the Galleria GSTs belong to the insect-specific Delta and Epsilon classes with 4 and 6 members each, respectively. However, in contrast to a comparable larval EST dataset of the generalist plant herbivore lepidopteran H. armigera  both the total number of GSTs identified and the strong overrepresentation of the insect-specific GSTs is much lower in Galleria (Additional file 5). The insect-specific Delta and Epsilon GST classes are often involved in detoxification of xenobiotics and the limited number of GSTs from those classes may point at the unique ecological niche and highly specialized diet of Galleria which is devoid of any (potentially toxic) plant secondary metabolites.
We have generated a comprehensive larval transcriptome map of the phylogenetically ancient lepidopteran Galleria mellonella. This data set complements and massively expands the known spectrum of immunity and stress related genes of this model host which have been found in previous studies using peptidomic  or SSH-based transcriptomic approaches . Besides genes encoding proteins that mediate recognition of microbial signatures such as GNBPs, βGRPs, PGRPs and Toll or immunity-related signalling, we determined a broad spectrum of defence related effectors such as antimicrobial peptides and proteins among which moricins and gloverins are restricted to Lepidoptera. In line with other studies, the spectrum of genes which is up-regulated in response to injected LPS includes proteins involved in detoxification (apolipoprotein D, cytochrome P450s, gluthathione S-transferases) and stress response (e.g. heat shock proteins). The secretion of induced immunity-and stress-related peptides and proteins into the hemolymph has been confirmed by comparative proteomic analysis of hemolymph samples from untreated and immunized larvae. Importantly, the spectrum of immunity-related genes identified in this study shares high similarity with that found in another lepidopteran species, the tobacco hornworm M. sexta, whose killed bacteria-induced transcriptome has previously been analyzed by pyrosequencing . Furthermore, except for attacins, we identified in Galleria members of all families of antimicrobial peptides which are predicted from the complete genome sequence of B. mori, . Taken together we postulate that all effector molecule families contributing to lepidopteran innate immunity are present in the phylogenetically basal family Pyralidae to which Galleria belongs. The entity of generated data provide a valuable platform for more detailed analyses of immune responses in Galleria and, therefore, improve the suitability of this lepidopteran both as a model host for human pathogens and for studies addressing coevolution with entomopathogens.
Galleria mellonella individuals used here were obtained from the laboratory culture which has been used in our previous studies. Galleria caterpillars were reared on an artificial diet (22% maize meal, 22% wheat germ, 11% dry yeast, 17.5% bee wax, 11% honey, and 11% glycerin) at 31°C in darkness. Last-instar larvae, each weighing between 250 and 350 mg, were used for immunization using 10 mg/ml LPS dissolved in water (Sigma, Taufkirchen, Germany). Ten microliters of sample volume per caterpillar was injected dorsolaterally into the hemocoel using 1-ml disposable syringes and 0.4-by 20-mm needles mounted on a microapplicator. Larvae were homogenized at 8 h postinjection for RNA isolation or bled at 24 h postinjection to obtain hemolymph samples.
RNA extraction, cDNA normalization and Next Generation Sequencing
Total RNA was extracted from different life stages, from hemocytes, and from immune-challenged larvae (injections) using TRIZOL and mRNA was subsequently isolated from total RNA using the MN-NucleoTrap mRNA kit according to the manufacturers' instructions (Macherey & Nagel). cDNAs were generated from 1 μg of poly(A)+ mRNA using the SMART PCR cDNA synthesis kit (BD Clontech) following the manufacturer's protocol. Reverse transcription was performed with the SMART KIT reverse transcriptase (Takara) for 60 min at 42°C. In order to prevent over-representation of the most common transcripts, the resulting single-stranded cDNAs were normalized using the Kamchatka crab duplex-specific nuclease method (Trimmer cDNA normalization kit, Evrogen) . Subsequently, SMART kit components and Triple-Taq enzyme with proof-reading activity were used to generate full-length enriched double-stranded long cDNAs. Each step of the normalization procedure was carefully monitored to avoid the generation of artefacts and overcycling. The optimal condition for ds-cDNA synthesis was empirically determined by subjecting the cDNA to a range of thermocycling numbers and their products checked by electrophoresis. The optimal cycle number was defined as the maximum number of PCR cycles without any signs of overcycling. The resulting normalized cDNA library was used for 454 pyrosequencing  using the Roche 454 FLX machine and Sanger sequencing using an ABI 3730 × l capillary sequencer. The 454 sequence reads were assembled using the newbler assembler with standard settings and using the CLC Genomics Workbench as an alternative assembly method. Before assembly, obtained reads were preprocessed by masking PolyA tails and removing SMART adapters using custom written Perl scripts. We compared the resulting contigs to the refseq protein database containing all information on coding sequences so far obtained (March 2010). Furthermore, we set up species specific databases from Drosophila, Bombyx, and human in order to find species specific similarities.
Sanger Sequencing and Generation of EST Databases
A fraction of the dscDNAs was cloned in the pGEM-T-easy vector. Ligations were transformed into E. coli ELECTROMAX DH5α-E electro-competent cells (Invitrogen). Plasmid minipreparation from bacterial colonies grown in 96 deep-well plates was performed using the 96 well robot plasmid isolation kit (NextTec) on a Tecan Evo Freedom 150 robotic platform (Tecan). Sequencing of both the 5' and 3' termini of cDNA library clones was carried out on an ABI 3730 xl automatic DNA sequencer (PE Applied Biosystems). Vector clipping, quality trimming and sequence assembly using stringent conditions (e.g. high quality sequence trimming parameters, 95% sequence identity cutoff, 25 bp overlap) was done with the Lasergene software package (DNAStar Inc.).
Blast homology searches and sequence annotation
We set up individual searchable databases for the complete sequence dataset and used this to identify the genes we describe in more detail in the text. Blast searches were conducted on a local server using the National Center for Biotechnology Information (NCBI) blastall program. Homology searches (BLASTx and BLASTn) of unique sequences and functional annotation by gene ontology terms (GO; http://www.geneontology.org), InterPro terms (InterProScan, EBI), enzyme classification codes (EC), and metabolic pathways (KEGG, Kyoto Encyclopedia of Genes and Genomes) were determined using the BLAST2GO software suite v2.3.1 http://www.blast2go.de. Homology searches were performed remotely on the NCBI server through QBLAST, and followed a sequential strategy. First, Sequences were searched against the NCBI non-redundant (nr) protein database using an E-value cut-off of 10-3, with predicted polypeptides of a minimum length of 15 amino acids. Second, sequences retrieving no BLASTx hit were searched again by BLASTn, against an NCBI nr nucleotide database using an E-value cut-off of 10-10. The GO data presented represent the level 3 analysis, illustrating general functional categories. Enzyme classification codes, and KEGG metabolic pathway annotations, were generated from the direct mapping of GO terms to their enzyme code equivalents. Finally, InterPro searches were performed remotely from BLAST2GO via the InterProEBI web server. In order to obtain a rough transcriptome coverage estimate for the Galleria larval cDNA library, we went through a series of search steps in order to i) obtain all hits against the conserved KEGG pathway database, and ii) estimate genome coverage by identifying the complete ribosomal protein dataset as compared to the full B. mori set. Based on these findings we estimate the theoretical transcriptome coverage to be close to 90% (e.g. 77/79 B. mori ribosomal proteins were found). Nucleotide sequences were analyzed in more detail using the commercial Lasergene Software package and the freeware BioEdit program. Genes were aligned by their amino acid sequences using the ClustalW function  or the MAFFT program. If necessary, alignments were then corrected by eye and reverted back to the nucleotide sequences for the phylogenetic analyses and in order to remove redundant contigs.
We have deposited the EST (Sanger) and short read (454 Roche) data with the following accession numbers: ERP000555 (SRA) and JG394435-JG406465 (dbEST). Phylogenetic data was deposited at TreeBASE with submission ID 11389. All of the predicted protein sequences used for alignments and phylogenies can be found in additional file 6. Note that the names of the validated proteins are made from the letters Gme followed by the number of the contig from the assembly. An assembly of the Galleria data with contig consensus sequences, Blast2GO hits against nr database, hit accessions, and annotations including InterPro scans can be found in Additional file 7.
The phylogenetic reconstruction implemented for the analysis of several proteins was performed using two different methods. For the Neighbour-Joining (NJ) method we implemented the TREECON program. Amino acid sequences were aligned by MAFFT http://mafft.cbrc.jp/alignment/server/index.html and each visually inspected for regions of high quality alignment. The NJ consensus tree was generated with TREECON. Distance calculations were performed after Tajima & Nei and bootstrap analysis, running 1000 bootstrap samples. Conserved residues in the alignments were highlighted with BOXSHADE 3.21 http://www.ch.embnet.org/software/BOX_form.html. In addition to the Neighbour-Joining method, for some gene trees the phylogenetic reconstruction was done by Bayesian inference using Mr. Bayes 3.1. The prior was set for the amino acid models to mix, thereby allowing model jumping between fixed-rate amino acid models. Markov Chain Monte Carlo runs were carried out for 10,000,000 generations after which log likelihood values showed that equilibrium had been reached after the first 5000 generations in all cases, and those data were discarded from each run and considered as 'burnin'. Two runs were conducted for the dataset showing agreement in topology and likelihood scores. The Neighbour-joining and the Bayesian tree topologies including their general subfamily relationships and node supports were in agreement. The gene trees were visualized and optimized with the MEGA4 software package .
Two-dimensional gel electrophoresis of hemolymph proteins
Proteomic analysis of immune hemolymph has been performed as described previously . In brief, hemolymph samples from 10 larvae 24 h post immune challenge and from 10 untreated larvae used as controls were collected directly into 1.5 ml pre-cooled plastic tubes containing traces of phenylthiourea to prevent melanisation reactions. Hemocytes were removed by brief centrifugation step and cell-free hemolymph was precipitated by the addition of 3 volumes of 100% acetone and 0.4 volumes of 100% trichloroacetic acid and incubation at 20°C for 1 h. After centrifugation at 20,000 × g for 10 min, the pellet was washed three times with 100% acetone and resolved under agitation in 8 M urea at 22°C for 16 h. Protein concentrations were determined using a Micro BC assay kit (Uptima, Montlucon, France). Two-dimensional gel electrophoresis was done with the Ettan IPGphor II system and the Ettan DALTsix electrophoresis unit (Amersham Biosciences, Uppsala, Sweden) according to the instructions of the manufacturer. Briefly, 1 mg of protein was mixed with immobilized pH gradient (IPG) buffer (pH 3 to 11 nonlinear gradient [NL]) and applied on an IPG strip (24 cm; pH 3 to 11 NL). Isoelectric focusing was performed at 20°C and 75 μA per IPG strip as follows: swelling for 24 h and isoelectric focusing for 1 h at 500 V, 8-h gradient to reach 1,000 V, 3-h gradient to reach 8,000 V, and isoelectric focusing for 4 h at 8,000 V. Prior to Tris-Tricine-sodium dodecyl sulfate (SDS)-polyacrylamide gel electrophoresis (5) with 26-by 20-cm 15% gels, the strips were equilibrated with 6 M urea, 30% glycerin, 2% SDS, and 50 mM Tris-HCl at pH 8 for 30 min. After electrophoresis at 20°C, the gels were stained using colloidal Coomassie brilliant blue (Carl Roth). For image analysis, the gels were scanned using an Umax PowerLook II scanner and analyzed with Delta2D software (Decodon, Greifswald, Germany). Spot identity has been determined by comparing spots with our recent study  combined with additional peptide mass fingerprinting analyses. In brief, spots excised from the gel were carbamidomethylated and in-gel digested using mass spectrometry grade trypsin (Promega) in 0.025 M NH4HCO3. The mass spectra of the resulting tryptic peptides were recorded using an Ultraflex TOF/TOF mass spectrometer (Bruker Daltonik, Bremen, Germany) operating under FlexControl 2.4 (Bruker) in the positive-ion reflectron mode, with dihydroxy benzoic acid as the matrix. Peptide mass profiles were analyzed with local Mascot http://www.matrixscience.com, using deduced protein sequences from our present Galleria transcriptome analysis database.
We acknowledge the help of M. Linder (University of Giessen) with Maldi-tof MS analysis, A. Berisha (University of Giessen) with local Mascot analysis and M. Fischer (University of Giessen) for technical assistance. The authors thank the German Research Foundation for a Heisenberg Fellowship to BA (AL902/4-1) and a grant within the DFG priority program 1399 (Host-parasite Coevolution-rapid reciprocal adaptation and its genetic base) to AV (VI 219/3-1). This work was supported by the Max Planck Society.
- Zou Z, Najar F, Wang Y, Roe B, Jiang H: Pyrosequence analysis of expressed sequence tags for Manduca sexta hemolymph proteins involved in immune responses. Insect Biochem Mol Biol. 2008, 38: 677-82. 10.1016/j.ibmb.2008.03.009.View ArticlePubMedPubMed CentralGoogle Scholar
- Pauchet Y, Wilkinson P, Vogel H, Nelson DR, Reynolds SE, Heckel DG, ffrench-Constant RH: Pyrosequencing the Manduca sexta larval midgut transcriptome: messages for digestion, detoxification and defence. Insect Mol Biol. 2010, 19: 61-75. 10.1111/j.1365-2583.2009.00936.x.View ArticlePubMedGoogle Scholar
- Vilcinskas A: Coevolution between pathogen-derived proteinases and proteinase inhibitors of host insects. Virulence. 2010, 1: 206-14. 10.4161/viru.1.3.12072.View ArticlePubMedGoogle Scholar
- Scully L, Bidochka M: Developing insects as models for current and emerging human pathogens. FEMS Microbiol Lett. 2006, 263: 1-9. 10.1111/j.1574-6968.2006.00388.x.View ArticlePubMedGoogle Scholar
- Fuchs B, Mylonakis E: Using non-mammalian hosts to study fungal virulence and host defense. Curr Opin Microbiol. 2006, 9: 346-51. 10.1016/j.mib.2006.06.004.View ArticlePubMedGoogle Scholar
- Kavanagh K, Reeves EP: Exploiting the potential of insects for the in vivo pathogenicity testing of microbial pathogens. FEMS Microbiol. 2004, 28: 101-12. 10.1016/j.femsre.2003.09.002.View ArticleGoogle Scholar
- Mylonakis E, Casadevall A, Ausubel FM: Exploiting amoeboid and non-vertebrate animal model systems to study virulence of human pathogenic fungi. Plos Path. 2007, 3: 3101-View ArticleGoogle Scholar
- Jander G, Rahme L, Ausubel FM, Drenkard E: Positive correlation between virulence of Pseudomonas aeroginosa mutants in mice and insects. J Bacteriol. 2000, 182: 3843-45. 10.1128/JB.182.13.3843-3845.2000.View ArticlePubMedPubMed CentralGoogle Scholar
- Fedhila S, Daou N, Lereclus D, Nielsen-LeRoux C: Identification of Bacillus cereus internalin and other candidate virulence genes specifically induced during oral infection in insects. Mol Microbiol. 2006, 62: 339-55. 10.1111/j.1365-2958.2006.05362.x.View ArticlePubMedGoogle Scholar
- Park S, Kim KM, Lee JH, Seo SJ, Lee IH: Extracellular gelatinase of Enterococcus faecalis destroys a defense system in insect hemolymph and human serum. Infect Immun. 2007, 75: 1861-9. 10.1128/IAI.01473-06.View ArticlePubMedPubMed CentralGoogle Scholar
- Mukherjee K, Altincicek B, Hain T, Domann E, Vilcinskas A, Chakraborty T: Galleria mellonella as model system to study Listeria pathogenesis. Appl Environ Microbiol. 2010, 76: 310-7. 10.1128/AEM.01301-09.View ArticlePubMedGoogle Scholar
- Miyata S, Casey M, Frank D, Ausubel F, Drenkard E: Use of the Galleria mellonella caterpillar as a model host to study the role of type III secretion system in Pseudomonas aeroginosa pathogenesis. Infect Immun. 2003, 71: 2404-13. 10.1128/IAI.71.5.2404-2413.2003.View ArticlePubMedPubMed CentralGoogle Scholar
- Garcia-Lara J, Needham A, Foster S: Invertebrates as animal models for Staphylococcus aureus pathogenesis: A window into host-pathogen interaction. FEMS Immunol Med Microbiol. 2005, 43: 311-23. 10.1016/j.femsim.2004.11.003.View ArticlePubMedGoogle Scholar
- Bergin D, Murphy L, Keenan J, Clynes M, Kavanagh K: Pre-exposure to yeast protects larvae of Galleria mellonella from a subsequent lethal infection by Candida albicans and is mediated by the increased expression of antimicrobial peptides. Microbes Infect. 2006, 8: 2105-12. 10.1016/j.micinf.2006.03.005.View ArticlePubMedGoogle Scholar
- Mylonakis E, Moreno R, El Khoury J, Idnurm A, Heitman J, Calderwood SB, Ausubel FM, Diener A: Galleria mellonella as a model system to study Cryptococcus neoformans pathogenesis. Infect Immun. 2005, 73: 3842-50. 10.1128/IAI.73.7.3842-3850.2005.View ArticlePubMedPubMed CentralGoogle Scholar
- Clermont A, Wedde M, Seitz V, Podsiadlowski L, Hummel M, Vilcinskas A: Cloning and expression of an inhibitor against microbial metalloproteinases from insects (IMPI) contributing to innate immunity. Biochem J. 2004, 382: 315-22. 10.1042/BJ20031923.View ArticlePubMedPubMed CentralGoogle Scholar
- Vilcinskas A: Anti-Infective therapeutics from the lepidopteran model host Galleria mellonella. Curr Pharm Des. 2011, 17 (13): 1240-5.View ArticlePubMedGoogle Scholar
- Seitz V, Clermont A, Wedde M, Hummel M, Vilcinskas A, Schlatterer K, Podsiadlowski L: Identification of immunorelevant genes from greater wax moth (Galleria mellonella) by a subtractive hybridization approach. Dev Comp Immunol. 2003, 27: 207-215. 10.1016/S0145-305X(02)00097-6.View ArticlePubMedGoogle Scholar
- Altincicek B, Vilcinskas A: Identification of immune-related genes from an apterygote insect, the firebrat Thermobia domestica. Insect Biochem Mol Biol. 2007, 37: 726-31. 10.1016/j.ibmb.2007.03.012.View ArticlePubMedGoogle Scholar
- Altincicek B, Vilcinskas A: Analysis of the immune-inducible transcriptome from microbial stress resistant, rat-tailed maggots of the drone fly Eristalis tenax. BMC Genomics. 2007, 8: 326-10.1186/1471-2164-8-326.View ArticlePubMedPubMed CentralGoogle Scholar
- Altincicek B, Lindner M, Linder D, Preissner K, Vilcinskas A: Microbial metalloproteinases mediate sensing of invading pathogens and activate innate immune responses in the lepidopteran model host Galleria mellonella. Infect Immun. 2007, 75: 175-83. 10.1128/IAI.01385-06.View ArticlePubMedGoogle Scholar
- Götz S, García-Gómez JM, Terol J, Williams TD, Nueda MJ, Robles M, Talón M, Dopazo J, Conesa A: High-throughput functional annotation and data mining with the Blast2GO suite. Nucleic Acids Res. 2008, 36: 3420-35. 10.1093/nar/gkn176.View ArticlePubMedPubMed CentralGoogle Scholar
- Beldade P, McMillan WO, Papanicolaou A: Butterfly genomics eclosing. Heredity. 2008, 100: 150-7. 10.1038/sj.hdy.6800934.View ArticlePubMedGoogle Scholar
- Whitten MM, Tew IF, Lee BL, Ratcliffe NA: A novel role for an insect apolipoprotein (apolipophorin III) in beta-1,3-glucan pattern recognition and cellular encapsulation reactions. J Immunol. 2004, 172 (4): 2177-85.View ArticlePubMedGoogle Scholar
- Altincicek B, Stötzel S, Wygrecka M, Preissner KT, Vilcinskas A: Host-derived extracellular nucleic acids enhance innate immune responses, induce coagulation, and prolong survival upon infection in insects. J Immunol. 2008, 181 (4): 2705-12.View ArticlePubMedGoogle Scholar
- Pauchet Y, Muck A, Svatos A, Heckel DG: Chromatographic and electrophoretic resolution of proteins and protein complexes from the larval midgut microvilli of Manduca sexta. Insect Biochem Mol Biol. 2009, 39: 467-74. 10.1016/j.ibmb.2009.05.001.View ArticlePubMedGoogle Scholar
- Pauchet Y, Freitak D, Heidel-Fischer HM, Heckel DG, Vogel H: Immunity or digestion: glucanase activity in a glucan-binding protein family from Lepidoptera. J Biol Chem. 2009, 284 (4): 2214-24.View ArticlePubMedGoogle 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 USA. 1996, 93: 7888-93. 10.1073/pnas.93.15.7888.View ArticlePubMedPubMed CentralGoogle Scholar
- Lemaitre B, Hoffmann JA: The host defense of Drosophila melanogaster. Annu Rev Immunol. 2007, 25: 697-743. 10.1146/annurev.immunol.25.022106.141615.View ArticlePubMedGoogle Scholar
- Ligoxygakis P, Bulet P, Reichhart JM: Critical evaluation of the role of the Toll-like receptor 18-Wheeler in the host defense of Drosophila. EMBO Rep. 2002, 3: 666-73. 10.1093/embo-reports/kvf130.View ArticlePubMedPubMed CentralGoogle Scholar
- Tanji T, Yun EY, Ip YT: Heterodimers of NF-kappaB transcription factors DIF and Relish regulate antimicrobial peptide genes in Drosophila. Proc Natl Acad Sci USA. 2010, 107: 14715-20. 10.1073/pnas.1009473107.View ArticlePubMedPubMed CentralGoogle Scholar
- An C, Jiang H, Kanost MR: Proteolytic activation and function of the cytokine Spätzle in the innate immune response of a lepidopteran insect, Manduca sexta. FEBS J. 2010, 277: 148-62. 10.1111/j.1742-4658.2009.07465.x.View ArticlePubMedGoogle Scholar
- Johnson S, Michalak M, Opas M, Eggleton P: The ins and outs of calreticulin: from the ER lumen to the extracellular space. Trends Cell Biol. 2001, 11: 122-9. 10.1016/S0962-8924(01)01926-2.View ArticlePubMedGoogle Scholar
- Altincicek B, Vilcinskas A: Comparative analysis of septic injury-inducible genes in phylogenetically distant model organisms of regeneration and stem cell research, the planarian Schmidtea mediterranea and the cnidarian Hydra vulgaris. Front Zool. 2008, 5: 6-10.1186/1742-9994-5-6.View ArticlePubMedPubMed CentralGoogle Scholar
- Levy S, Shoham T: The tetraspanin web modulates immune-signalling complexes. Nat Rev Immunol. 2005, 5: 136-48. 10.1038/nri1548.View ArticlePubMedGoogle Scholar
- Brown SE, Howard A, Kasprzak AB, Gordon KH, East PD: The discovery and analysis of a divergend family of novel antifungal moricin-like peptides in the wax moth Galleria mellonella. Insect Biochem Mol Biol. 2008, 38: 201-12. 10.1016/j.ibmb.2007.10.009.View ArticlePubMedGoogle Scholar
- Hara S, Yamakawa M: Moricin, a novel type of antibacterial peptide isolated from the silkworm, Bombyx mori. J Biol Chem. 1995, 270: 29923-7. 10.1074/jbc.270.50.29923.View ArticlePubMedGoogle Scholar
- Axén A, Carlsson A, Engström A, Bennich H: Gloverin, an antibacterial protein from the immune hemolymph of Hyalophora pupae. Eur J Biochem. 1997, 247: 614-9. 10.1111/j.1432-1033.1997.00614.x.View ArticlePubMedGoogle Scholar
- Cheng T, Zhao P, Liu C, Xu P, Gao Z, Xia Q, Xiang Z: Structures, regulatory regions, and inductive expression patterns of antimicrobial peptide genes in the silkworm Bombyx mori. Genomics. 2006, 87: 356-65. 10.1016/j.ygeno.2005.11.018.View ArticlePubMedGoogle Scholar
- Steiner H, Hultmark D, Engström A, Bennich H, Boman HG: Sequence and specificity of two antibacterial proteins involved in insect immunity. Nature. 1981, 292: 246-8. 10.1038/292246a0.View ArticlePubMedGoogle Scholar
- Kim CH, Lee JH, Kim I, Seo SJ, Son SM, Lee KY, Lee IH: Purification and cDNA cloning of a cecropin-like peptide from the great wax moth, Galleria mellonella. Mol Cells. 2004, 17: 262-6.PubMedGoogle Scholar
- Lee YS, Yun EK, Jang WS, Kim I, Lee JH, Park SY, Ryu KS, Seo SJ, Kim CH, Lee IH: Purification, cDNA cloning and expression of an insect defensin from the great wax moth, Galleria mellonella. Insect Mol Biol. 2004, 13: 65-72. 10.1111/j.1365-2583.2004.00462.x.View ArticlePubMedGoogle Scholar
- Langen G, Imani J, Altincicek B, Kieseritzky G, Kogel KH, Vilcinskas A: Transgenic expression of gallerimycin, a novel antifungal insect defensin from the greater wax moth Galleria mellonella, confers resistance to pathogenic fungi in tobacco. Biol Chem. 2006, 387: 549-57. 10.1515/BC.2006.071.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Girard PA, Boublik Y, Wheat CW, Volkoff AN, Cousserans F, Brehélin M, Escoubas JM: X-tox: an atypical defensin derived family of immune-related proteins specific to Lepidoptera. Dev Comp Immunol. 2008, 32: 575-84. 10.1016/j.dci.2007.09.004.View ArticlePubMedGoogle Scholar
- Mohrig W, Messner B: Immunreaktionen bei Insekten. I. Lysozyme als grundlegender antibakterieller Faktor im humoralen Abwehrgeschehen. Biol Zentralbl. 1968, 87: 439-47.Google Scholar
- Jolles J, Schoentgen F, Croizier G, et al: Insect lysozymes from three species of Lepidoptera: Their structural relatedness to the C (chicken) type lysozyme. J Mol Evol. 1979, 14: 267-71. 10.1007/BF01732494.View ArticlePubMedGoogle Scholar
- Yu KH, Kim KN, Lee JH, Lee HS, Kim SH, Cho KY, Nam MH, Lee IH: Comparative study and characteristics of lysozymes from the hemolymph of three lepidopteran larvae, Galleria mellonella, Bombyx mori, Agrius convolvuli. Dev Comp Immunol. 2002, 26: 707-13. 10.1016/S0145-305X(02)00027-7.View ArticlePubMedGoogle Scholar
- Vilcinskas A, Matha V: Effect of the entomopathogenic fungus Beauveria bassiana on humoral immune response of Galleria mellonella larvae (Lepidoptera: Pyralidae). Eur J Entomol. 1997, 94: 461-72.Google Scholar
- Vilcinskas A, Götz P: Parasitic fungi and their interactions with the insect immune system. Adv Parasitol. 1999, 43: 267-313.View ArticleGoogle Scholar
- Nilsen IW, Overbø K, Sandsdalen E, Sandaker E, Sletten K, Myrnes B: Protein purification and gene isolation of chlamysin, a cold-active lysozyme-like enzyme with antibacterial activity. FEBS Lett. 1999, 464: 153-8. 10.1016/S0014-5793(99)01693-2.View ArticlePubMedGoogle Scholar
- Brown SE, Howard A, Kasprzak AB, Gordon KH, East PD: A peptidomic study reveals the impressive antimicrobial peptide arsenal of the wax moth Galleria mellonella. Insect Biochem Mol Biol. 2009, 39: 792-800. 10.1016/j.ibmb.2009.09.004.View ArticlePubMedGoogle Scholar
- Cytrynska M, Zdybicka-Barabas A, Suder P, Jakubowicz T: Purification and characterization of eight peptides from Galleria mellonella immune hemolymph. Peptides. 2007, 28: 533-546. 10.1016/j.peptides.2006.11.010.View ArticlePubMedGoogle Scholar
- Fröbius AC, Kanost MR, Götz P, Vilcinskas A: Isolation and characterization of novel inducible serine protease inhibitors from larval hemolymph of the greater wax moth Galleria mellonella. Eur J Biochem. 2000, 267: 2046-53. 10.1046/j.1432-1327.2000.01207.x.View ArticlePubMedGoogle Scholar
- Vilcinskas A, Wedde M: Insect inhibitors of metalloproteinases. IUBMB Life. 54: 339-43.Google Scholar
- Adekoya O, Sylte I: The thermolysin family (M4) of enzymes: therapeutic and biotechnological potential. Chem Biol Drug Des. 2009, 73: 7-16. 10.1111/j.1747-0285.2008.00757.x.View ArticlePubMedGoogle Scholar
- Wedde M, Weise C, Kopacek P, Franke P, Vilcinskas A: Purification and characterization of an inducible metalloprotease inhibitor from the hemolymph of greater wax moth larvae, Galleria mellonella. Eur J Biochem. 1998, 255: 534-43.View ArticleGoogle Scholar
- Wedde M, Weise C, Nuck C, Altincicek B, Vilcinskas A: The insect metalloproteinase inhibitor gene of the lepidopteran Galleria mellonella encodes two distinct inhibitors. Biol Chem. 2007, 388: 119-27. 10.1515/BC.2007.013.View ArticlePubMedGoogle Scholar
- Kumagai T, Awai M, Okada S: Mobilization of iron and iron-related proteins in rat spleen after intravenous injection of lipopolysaccharides (LPS). Pathol Res Pract. 1992, 188: 931-41.View ArticlePubMedGoogle Scholar
- Yun EY, Lee JK, Kwon OY, Hwang JS, Kim I, Kang SW, Lee WJ, Ding JL, You KH, Goo TW: Bombyx mori transferrin: genomic structure, expression and antimicrobial activity of recombinant protein. Dev Comp Immunol. 2009, 33: 1064-9. 10.1016/j.dci.2009.05.008.View ArticlePubMedGoogle Scholar
- Altincicek B, Vilcinskas A: Analysis of the immune-related transcriptome of a lophotrochozoan model, the marine annelid Platynereis dumerilii. Front Zool. 2007, 4: 18-10.1186/1742-9994-4-18.View ArticlePubMedPubMed CentralGoogle Scholar
- Ranson H, Claudianos C, Ortelli F, Abgall C, Hemingway J, Sharakhova MV, Unger MF, Collins FH, Feyereisen R: Evolution of supergene families associated with insecticide resistance. Science. 2002, 298 (5591): 179-81. 10.1126/science.1076781.View ArticlePubMedGoogle Scholar
- Ding Y, Ortelli F, Rossiter LC, Hemingway J, Ranson H: The Anopheles gambiae glutathione transferase supergene family: annotation, phylogeny and expression profiles. BMC Genomics. 2003, 4: 35-10.1186/1471-2164-4-35.View ArticlePubMedPubMed CentralGoogle Scholar
- Caccuri AM, Antonini G, Board PG, Flanagan J, Parker MW, Paolesse R, Turella P, Chelvanayagam G, Ricci G: Identification, characterization, and crystal structure of the omega class glutathione transferases. J Biol Chem. 2000, 275 (32): 24798-806. 10.1074/jbc.M001706200.View ArticleGoogle Scholar
- Yu SJ: Insect glutathione S-transferases. Zoological Studies. 1996, 35 (1): 9-19.Google Scholar
- Francis F, Vanhaelen N, Haubruge E: Glutathione S-transferases in the adaptation to plant secondary metabolites in the Myzus persicae aphid. Arch Insect Biochem Physiol. 2005, 58: 166-174. 10.1002/arch.20049.View ArticlePubMedGoogle Scholar
- Dabrowska P, Freitak D, Vogel H, Heckel DG, Boland W: The phytohormone precursor OPDA is isomerized in the insect gut by a single, specific glutathione transferase. Proc Natl Acad Sci USA. 2009, 106 (38): 16304-9. 10.1073/pnas.0906942106.View ArticlePubMedPubMed CentralGoogle Scholar
- Brown SE, Howard A, Kasprzak AB, Gordon KH, East PD: A peptidomic study reveals the impressive antimicrobial peptide arsenal of the wax moth Galleria mellonella. Insect Biochem Mol Biol. 2009, 39: 792-800. 10.1016/j.ibmb.2009.09.004.View ArticlePubMedGoogle Scholar
- Zhulidov PA, Bogdanova EA, Shcheglov AS, Vagner LL, Khaspekov GL, Kozhemyako VB, Matz MV, Meleshkevitch E, Moroz LL, Lukyanov SA, Shagin DA: Simple cDNA normalization using kamchatka crab duplex-specific nuclease. Nucleic Acids Res. 2004, 32: e37-10.1093/nar/gnh031.View ArticlePubMedPubMed CentralGoogle Scholar
- Margulies M, Egholm M, Altman WE, Attiya S, Bader JS, Bemben LA, Berka J, Braverman MS, Chen YJ, Chen Z, Dewell SB, Du L, Fierro JM, Gomes XV, Godwin BC, He W, Helgesen S, Ho CH, Irzyk GP, Jando SC, Alenquer ML, Jarvie TP, Jirage KB, Kim JB, Knight JR, Lanza JR, Leamon JH, Lefkowitz SM, Lei M, Li J, et al: Genome sequencing in microfabricated high-density picolitre reactors. Nature. 2005, 437: 376-80.PubMedPubMed CentralGoogle Scholar
- Conesa A, Götz S: Blast2GO: A Comprehensive Suite for Functional Analysis in plant genomics. Int J Plant Genom. 2008, 619832-Google Scholar
- Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19: 1572-4. 10.1093/bioinformatics/btg180.View ArticlePubMedGoogle Scholar
- Tamura K, Dudley J, Nei M, Kumar S: MEGA4: molecular evolutionary genetics analysis (MEGA) software version 4.0. Mol Biol Evol. 2007, 24: 1596-9. 10.1093/molbev/msm092.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.