Email updates

Keep up to date with the latest news and content from BMC Evolutionary Biology and BioMed Central.

Open Access Research article

A proteomic approach for studying insect phylogeny: CAPA peptides of ancient insect taxa (Dictyoptera, Blattoptera) as a test case

Steffen Roth13, Bastian Fromm1, Gerd Gäde2 and Reinhard Predel1*

Author Affiliations

1 Institute of Zoology, University of Jena, Erbertstrasse 1, D-07743 Jena, Germany

2 Zoology Department, University of Cape Town, Rondebosch 7701, South Africa

3 Institute of Biology, University of Bergen, Bergen N-5020, Norway

For all author emails, please log on.

BMC Evolutionary Biology 2009, 9:50  doi:10.1186/1471-2148-9-50


The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2148/9/50


Received:6 October 2008
Accepted:3 March 2009
Published:3 March 2009

© 2009 Roth et al; licensee 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.

Abstract

Background

Neuropeptide ligands have to fit exactly into their respective receptors and thus the evolution of the coding regions of their genes is constrained and may be strongly conserved. As such, they may be suitable for the reconstruction of phylogenetic relationships within higher taxa. CAPA peptides of major lineages of cockroaches (Blaberidae, Blattellidae, Blattidae, Polyphagidae, Cryptocercidae) and of the termite Mastotermes darwiniensis were chosen to test the above hypothesis. The phylogenetic relationships within various groups of the taxon Dictyoptera (praying mantids, termites and cockroaches) are still highly disputed.

Results

Tandem mass spectrometry of neuropeptides from perisympathetic organs was used to obtain sequence data of CAPA peptides from single specimens; the data were analysed by Maximum Parsimony and Bayesian Interference. The resulting cladograms, taking 61 species into account, show a topology which is in general agreement with recent molecular and morphological phylogenetic analyses, including the recent phylogenetic arrangement placing termites within the cockroaches. When sequence data sets from other neuropeptides, viz. adipokinetic hormones and sulfakinins, were included, the general topology of the cladogram did not change but bootstrap values increased considerably.

Conclusion

This study represents the first comprehensive survey of neuropeptides of insects for solely phylogenetic purposes and concludes that sequences of short neuropeptides are suitable to complement molecular biological and morphological data for the reconstruction of phylogenetic relationships.

Background

Peptides are short proteins, whose power to resolve phylogenetic questions have already been recognized (e.g. [1,2]). Peptide mass fingerprints support species recognition in many cases, particularly in organisms that exhibit few morphological differences such as microorganisms [3]. A specific group of peptides are the neuropeptides, structurally diverse messenger molecules, which influence a wide-range of physiological processes [4]. Due to their role as ligands, which have to fit into the respective receptors, neuropeptides are under considerable evolutionary constraint. Consequently, the regions of neuropeptide genes encoding for mature peptides may be highly conserved and suitable for the reconstruction of deep level phylogenetic relationships within higher taxa. However, only few attempts have been made to use these substances for phylogenetic purposes. Gäde [5] first introduced this approach for neuropeptides belonging to adipokinetic/hypertrehalosaemic hormones. The few sequence variations of these hormones within insects, however, do not contain sufficient information for a detailed analysis of phylogenetic relationships, although grouping of certain taxa is possible [6,7]. Other peptide families with multiple forms such as allatostatins [8] have both conserved and fast-evolving peptide sequences and are certainly more significant in this context but less extensively studied. The conserved sequences may be suitable for the reconstruction of phylogenetic relationships within higher taxa and the fast-evolving sequences may be more suitable for the reconstruction of tip-level phylogenetic relationships within closely related taxa.

Conducting a phylogenetic analysis of the genes encoding neuropeptides is not an easy task. In most cases, only small portions of these genes have been highly conserved, specifically the regions encoding for mature peptides, which interact with their receptors. Thus, primers successfully used for the identification of neuropeptide genes in a certain insect species may fail to recognize the orthologous gene in a related species (Derst, Roth, Predel; unpublished). Recent developments in mass spectrometric techniques [9], however, have paved the way for a rapid identification of mature neuropeptides from single insect specimens [10-13], thereby circumventing the genomic approach.

In the present study, tandem mass spectrometry was used for the first time to perform an extensive phylogenetic study on neuropeptides of insects, focusing on CAPA peptides of Dictyoptera. CAPA peptides were first identified from the American cockroach, Periplaneta americana [14-16]. CAPA-genes are known from a number of holometabolous insects (e.g. Drosophila melanogaster: [17], Anopheles gambiae: [18], Apis mellifera: [19], Tribolium castaneum: [20]). These genes encode for up to four peptides, which belong to CAPA-periviscerokinins (PVKs) and CAPA-pyrokinins (PKs). Both groups of CAPA peptides bind to different receptor types [21,22]. Besides their expression in a few interneurons, CAPA peptides are always part of the neuroendocrine system of the abdominal ventral nerve cord and are likely released into the haemolymph via abdominal perisympathetic organs (PSOs). Direct mass spectrometric screening of these organs (see [10,23]) allowed the unambiguous identification of the CAPA peptides from single specimens and cleared the way for a large-scale screening of these neuropeptides in the taxon Dictyoptera.

The taxon Dictyoptera includes praying mantids (Mantodea), termites (Isoptera), and cockroaches (Blattoptera) (e.g. [24]), and members are among the oldest pterygote insects known. Both morphological and molecular data support a monophyly of Mantodea and Isoptera (see [25]). The relationships of Mantodea, Isoptera, and Blattoptera, the monophyly of Blattoptera and the relationships among several cockroach lineages are, however, a topic of conflicting conclusions (e.g. [26-38]). In particular, the position of the genus Cryptocercus within the Blattoptera and its relationship with Isoptera has been the focus of numerous phylogenetic studies. Grandcolas (analysis of morpho-anatomical data: [39,40]) and Gäde et al. (analysis of adipokinetic hormones: [6]) placed these wood-feeding cockroaches in the Polyphagidae. Molecular data, however, suggest a sister-relationship between termites and Cryptocercus [41-44], a historical position [45] that is supported by Deitz et al. [46], and Klass & Meier [47] based on morpho-anatomical data. Inward et al. [44] presented convincing data to suggest that Isoptera nest within Blattoptera. The monophyly of several cockroach taxa and subgroups of these taxa is, however, doubtful. In a recent analysis of five gene loci, Inward et al. [44] found no support for the monophyly of the Blattellidae and subordinated taxa within the Blaberidae. In some of these taxa, further data acquisition of conventional molecular and morphological characters and more species may provide sufficient information to resolve more precisely the phylogenetic relationship of certain taxa within the Blattoptera (see [44,47]). In cases where these attempts result in conflicting hypotheses about the placement and monophyly of different taxa, additional characters (e.g. sequences of neuropeptides) may be required to test the robustness of the different analyses.

To test the phylogenetic information of neuropeptides in general, we used a stepwise approach by analysing the topology and stability of the phylogenetic trees, starting with the CAPA peptide data set followed by repeated analyses with additional neuropeptide sequences, namely adipokinetic hormone (AKH-1) and sulfakinins (SKs).

The cladograms obtained from these peptide sequences confirmed that certain neuropeptide sequences of insects are able to complement molecular, biological and morphological data for the reconstruction of phylogenetic relationships.

Results

Data acquisition and alignment

Direct mass spectrometric analysis of abdominal PSO preparations of single specimens (examples given in Figures 1, 2) revealed complete sequences of CAPA peptides from 61 cockroach/termite species. The species list covers major taxa of cockroaches (Blattidae, Polyphagidae, Cryptocercidae, Blaberidae, Blattellidae) and the termite Mastotermes darwiniensis. From most species, three CAPA-periviscerokinins (PVKs), and a single CAPA-pyrokinin (PK) were sequenced. Cryptocercus and the blattellid cockroaches Symploce pallens and Loboptera decipiens express only two different PVKs. A fourth PVK (designated PVK-4) was found in the Madagascan Blaberidae and the Table Mountain cockroach Aptera fusca (for sequences see [48]. These PVK-4 peptides, whose sequences suggest an internal gene duplication of PVK-1 (Elliptorhina, Gromphadorhina, Princisia) or PVK-2 (Aptera), did not influence the topology of phylogenetic trees and were not included in the final alignments. The average size of the PVKs was 11 amino acids (aa) and that of the PK was 17 aa. Sequences of the CAPA peptides were combined for each species and aligned (Table 1). With the inclusion of gaps and sequences of the outgroup species (Locusta migratoria and Drosophila melanogaster), the alignment resulted in 58 characters. Thirteen characters were constant, 12 variable characters were parsimony-uninformative, and 33 variable characters were parsimony-informative. The sequence of PVK-2 was found to be highly conserved and did not contain phylogenetically informative substitutions.

thumbnailFigure 1. Comparison of MALDI-TOF mass spectra (mass range 1000–2000 Da) of single abdominal PSO preparations of Blaberus craniifer and Bantua robusta (= peptide hormone fingerprint). Only few abundant substances are detectable. Underlying sequences were used for phylogenetic analyses.

thumbnailFigure 2. CID spectrum (ESI-QTOF MS) of Blattella germanica PVK-1 at [M+2H]2+ 536.9 ([M+H]+:1072.6). The y- and b-type fragment ions are labelled. Fragments were analyzed manually and the resulting sequence is given in the inset.

Table 1. Sequences of CAPA peptides aligned with Clustal X

Sequence variation of CAPA peptides within and among populations

We did not observe a single sequence variation of CAPA peptides from males, females, and larvae within any of the cockroach populations investigated. The PSOs of the American cockroach, P. americana, served as control in most mass spectrometric analyses (n > 400), and there was a lack of variability of neuropeptides at the individual level. We compared the CAPA peptides for a number of species (Diploptera punctata, Loboptera decipiens, Blaberus craniifer) that had been raised in a culture for multiple generations with specimens collected in the field. In addition, three South African populations of Bantua robusta that were collected in the rainforest (Tsitsikamma), fynbos (Cape Town), and Karroo vegetation (Kamieskroon) were investigated but no sequence variations were found (data not shown).

Analysis of phylogenetic relationships by means of CAPA sequences

Due to the high level of conservation in the sequences of PVK-2 as well as in the C-termini of the other CAPA peptides, only 33 amino acid positions contained phylogenetically informative characters. It was intriguing to see that the Maximum Parsimony (MP) analysis (Figure 3) obtained from these data was generally in agreement with recent molecular [44] and morphological [47] analyses, although the bootstrap values were relatively low. Significant support (bootstrapping, posterior probabilities of Bayesian analysis) was found for the monophyly of Blaberoidea (Blattellidae + Blaberidae) and Blattidae. The cladograms also support sister-group relationships between Blaberoidea and Blattoidea, Blattellidae and Blaberidae, and Blattidae and Polyphagidae + Cryptocercidae + Mastotermes. Within the latter clade, the three polyphagid species (Polyphaga aegyptiaca, Ergaula capucina, Therea petiveriana) appear as a monophyletic group separated from an unsolved sister-group containing Cryptocercus kyebangensis and Mastotermes darwiniensis. A Bayesian consensus tree (see additional file 1): Phylogenetic relationships based on neuropeptide sequences represented by a Bayesian majority rules consensus tree) yielded almost identical topologies with those that were obtained from Maximum Parsimony.

Additional file 1. Bayesian analysis of CAPA peptides sequences. Phylogenetic relationships of cockroaches based on CAPA peptides sequences represented by a Bayesian majority rules consensus tree. Numbers on the nodes indicate posterior probability values (≥ 0.49) (proportion of the 18205 sampled trees that contain the node).

Format: PDF Size: 17KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

thumbnailFigure 3. Phylogenetic relationships of cockroaches based on CAPA peptide sequences represented by a maximum parsimony (MP) 50% majority rules consensus tree. Numbers on the branches indicate bootstrap values (≥ 50) for MP. Italic numbers on the nodes indicate posterior probability values (≥ 0.5) (proportion of the 18205 sampled trees that contain the node). Tree length = 142, Consistency index (CI) = 0.768, Homoplasy Index (HI) = 0.232, Retention index (RI) = 0.907, Rescaled consistency index (RC) = 0.696.

Although the relationships within the Blaberidae (members of 9 of 11 subfamilies were included in this study) were poorly resolved, the different clades comprised, with few exceptions, only members of specific subfamilies. This was found for Blaberinae (Blaberus, Eublaberus, Lucihormetica, Archimandrita and Blaptica), Oxyhaloinae (Madagascan genera Princisia, Elliptorhina, Gromphadorhina as well as Rhyparobia), Panchlorinae (Panchlora species), and Perisphaeriinae (Southern African Cyrtotria, Perisphaeria, Bantua, Hostilia, and Pilema). In contrast, monophyly was not supported for some genera which are currently grouped in the Perisphaeriinae (see [49]). The Australian genus Laxta and African genus Derocalymma, both containing extremely flattened cockroaches which are adapted for living under bark, did not show close relationships with each other or with the remaining Perisphaeriinae. Instead, Derocalymma was found in a clade also containing Gyna and the Table Mountain cockroach, Aptera fusca. Blepharodera discoidalis, which was removed from the Perisphaeriinae by Grandcolas [49], contained CAPA peptides typical of Perisphaeriinae.

To test if the topology of the phylogenetic trees remains stable, the phylogenetic analysis was repeated with additional neuropeptide sequences, namely adipokinetic hormone (AKH-1) and sulfakinins (SKs) (see [50]). These peptides are stored in the corpora cardiaca, and mass fingerprints from these organs were sufficient for the correct assignment of the group-specific sequences in all cases. The resulting cladograms confirmed the topology of the former analysis, and increased the bootstrap values (Figure 4 and see additional file 2): Phylogenetic relationships based on peptide sequences represented by a Bayesian majority rules consensus tree).

Additional file 2. Bayesian analysis of CAPA peptides, AKH-1 and sulfakinin sequences. Phylogenetic relationships of cockroaches based on CAPA peptides, AKH-1, and sulfakinin sequences represented by a Bayesian majority rules consensus tree. Numbers on the nodes indicate posterior probability values (≥ 0.49) (proportion of the 20206 sampled trees that contain the node).

Format: PDF Size: 16KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

thumbnailFigure 4. Phylogenetic relationships of cockroaches based on CAPA peptides, AKH-1 and sulfakinin sequences represented by a maximum parsimony (MP) 50% majority rules consensus tree. Numbers on the branches indicate bootstrap values (≥ 50) for MP. Italic numbers on the nodes indicate posterior probability values (≥ 0.5) (proportion of the 20206 sampled trees that contain the node). Tree length = 181, Consistency index (CI) = 0.796, Homoplasy Index (HI) = 0.204, Retention index (RI) = 0.917, Rescaled consistency index (RC) = 0.729.

Discussion

The current investigation represents the first comprehensive survey of neuropeptides of insects for entirely phylogenetic purposes. Although the introduction of novel characters is consistently requested to corroborate existing hypotheses on phylogenetic relationships in insects (see [51]), such new character sets and methods have to compete with well-established methods. In order for our methodological approach to be acceptable by systematists using established methods, we developed techniques that allowed us to sample sufficient taxa and perform the analysis quickly.

In recent years, MALDI-TOF mass spectrometric analysis has been routinely used for studying the peptidome of the neuroendocrine system of insects [11-13,51-54]. The power of modern mass spectrometry means that only a few specimens of insects as small as the red flour beetle Tribolium castaneum are necessary to confirm the expression of more than 60 neuropeptides when genome information is available [20]. In the present study, however, genome information was not available, and the homologous peptides of the different species had to be de-novo sequenced. This approach posed a bioanalytical challenge and required a decision about the neuropeptide species to be included before extensive taxon sampling. The decision to select CAPA peptides first (see [48]) was made because these peptides fulfil certain criteria for a successful reconstruction of phylogenetic relationships. First, these peptides occur at high concentrations in neurohaemal tissues (abdominal PSOs), which are fairly easy to dissect, do not contain other neuropeptides at high concentrations and, thus, allow sequence elucidation from PSOs from a single specimen. Moreover, the detection of specific neuropeptide gene products, such as CAPA peptides, from defined neurohaemal organs usually excludes the alignment of peptides with sequence similarities that result from convergent evolution (homoplasy). Second, multiple members of related peptides encoded by single genes exist in insects. If the number of these often closely related peptide paralogues differs between related species, alignments may become difficult. Hence, it is more convenient to use a peptide family that contains the same number of peptide forms in the taxa of interest. In such a case, the storage organ as well as the conserved sequences of the peptide hormones can be used to assign the homologous peptides. Several peptide families were initially included in preliminary experiments; the CAPA peptides met the aforementioned criteria best and were thus used for this phylogenetic study. Since the sequence information from these peptides spans a length of 50 amino acids only, the resulting phylogenetic tree shows low posterior probabilities and low bootstrap levels.

In a subsequent and very rapid experimental approach, we used mass fingerprint data to include further neuropeptide sequences from relatively conserved peptides (AKH-1 and sulfakinins) in the phylogenetic analyses. The resulting topology of the cladograms did not change but the bootstrap values increased considerably. Since the additional neuropeptides did not differ very much between closely related taxa or did not differ at all, bootstrap levels of higher taxa were higher than those within lower taxa. This supports the hypothesis that, as a result of the decelerated co-evolution of neuropeptides and their receptors, neuropeptide sequences may be particularly suitable for the reconstruction of phylogenetic relationships within higher taxa.

The cladograms in Figures 3 and 4 show a topology that is in general agreement with recent molecular [44] and morphological phylogenetic analyses [47], including the recent phylogenetic arrangement placing termites within the cockroaches. Questions arising from the current data are: how can we solve existing polytomies, how can we enhance bootstrap supports for existing clades, and how can we possibly extend the analysis to higher or lower taxa? Sampling more taxa and only analysing CAPA peptides, AKHs and sulfakinins is unlikely to provide sufficient data to solve all of these questions. A combination of well chosen taxa sampling (including the outgroup taxa) and other neuropeptides will be needed to solve the relationship among the major lineages of Dictyoptera.

At a lower taxonomic level, however, a higher number of analyzed species in well-defined groups (e.g. Perisphaeriinae) may provide sufficient information to re-assess the generic composition of that group. Our data regarding the Perisphaeriinae differ, in part, from the suggestions made by Grandcolas [49], who analyzed head morphology and genitalia. The data do not support the removal of Blepharodera from this subfamily (see also [55]), and do not verify a close relatedness of Derocalymma and Laxta with the other genera of Perisphaeriinae. Indeed, we found six genera of Perisphaeriinae with completely identical neuropeptide sequences (Perisphaeria, Blepharodera, Pilema, Hostilia, Bantua, Cyrtotria) and these are exactly the genera which were placed in a single tribe (Perisphaeriini) by Roth [56].

We did not test how the choice of outgroup and ingroup taxa affects tree topology but further taxon sampling seems to be essential in termites and blattellid cockroaches. For the latter taxon, we have already obtained partial sequences from further species (unpublished data), which support the para- or polyphyletic origin of this group. In most cockroach groups (e.g. Blaberidae), however, even a more representative and comprehensive incorporation of further taxa is unlikely to provide novel insights into phylogenetic relationships. In these cases, further peptide families have to be included for phylogenetic analyses. In the present initial attempt, seven homologous neuropeptides of 61 species of Blattoptera were tested.

From a single cockroach, P. americana, roughly 80 neuropeptides have been elucidated by biochemical methods in recent years. Today, most of these peptides can be identified by mass spectrometric techniques as described in this manuscript, which makes these peptides generally suitable for phylogenetic studies. Fast evolving neuropeptides such as FMRFamides [57] can provide phylogenetic information at the generic level (see Figure 5) [58] but are not suitable for studying the deep level relationships of higher taxa within an insect order because the homology of such peptide copies among far related taxa can be difficult to assess. Other peptide families with multiple members, such as tachykinin-related peptides, pyrokinins, and allatostatins [59] are likely to be most suitable for the incorporation in phylogenetic analyses. These peptide families are represented by more than 30 paralogues in P. americana. Previous experiments have already shown that members of the Blattidae, which cannot be further separated from each other by the analysis of CAPA peptides, AKHs, and sulfakinins, are clearly distinguishable if species-specific pyrokinin sequences are identified [23]. These findings confirm that even short neuropeptide sequences of insects are suitable to complement molecular biological and morphological data for the reconstruction of phylogenetic relationships.

thumbnailFigure 5. MALDI-TOF mass spectra (neuropeptide mass fingerprints) from single thoracic PSO preparations of three Blaberus/Eublaberus species, representing FMRFamide related peptides which accumulate in the neurohaemal organs of insects (see Predel et al. 2004). The selected species were not distinguishable by screening the CAPA peptides from abdominal PSOs. All ion signals different from those of B. craniifer are marked. Such fingerprint data exist from all neurohaemal organs of all cockroach species investigated.

Conclusion

The phylogenetic relationships within the major lineages of cockroaches (Blaberidae, Blattellidae, Blattidae, Polyphagidae, Cryptocercidae) and their relationship to termites (Isoptera) were reconstructed by using the first comprehensive survey of neuropeptides of insects for solely phylogenetic purposes. The cladograms resulting from the analysis of peptide sequences of 61 Blattoptera species show a topology which is in general agreement with recent molecular and morphological phylogenetic analyses and also confirm the grouping of Isoptera within Blattoptera. Regarding other hypotheses about cockroach phylogeny, our data support the monophyly of Blaberoidea (Blattellidae + Blaberidae) and Blattidae. The cladograms also support sister-group relationships between Blaberoidea and a monophylum of the remaining cockroaches (including Isoptera), paraphyletic Blattellidae and Blaberidae, and Blattidae and Polyphagidae + Cryptocercidae + Mastotermes. This study verified that sequences of several neuropeptide families can complement molecular biological and morphological data for the reconstruction of phylogenetic relationships.

Methods

(a) Insects

In total, 61 species of Dictyoptera, representing the five cockroach taxa Polyphagidae, Cryptocercidae, Blattidae, Blattellidae, Blaberidae, and the termite Mastotermes darwiniensis were analyzed. Locusta migratoria (Orthoptera) and Drosophila melanogaster (Diptera) were used as outgroup species; the CAPA peptides of these species were identified by Predel & Gäde [60], Clynen et al. [61] and Kean et al. [17]. The names and places of collection (or sources of cockroach/termite cultures) of all species examined in this study, as well as the SWISSPROT accession numbers for peptide sequences are given in additional file 3. For most of the species, a mass fingerprint which represented about 40 peptide hormones was obtained from the major hormone release sites (corpora cardiaca, thoracic and abdominal perisympathetic organs). The respective fingerprints are typical of very closely related species (see Figure 5) and may be species-specific (see [62,63]). Remains of the insects as well as the fingerprint data can be obtained from the corresponding author.

Additional file 3. Additional information about studied species and accession numbers of peptides. Information about the species used in this study, including collecting sites/source, accession numbers of CAPA-peptides, AKH-1, sulfakinin-1 to UniProt. The sequence of sulfakinin-2 (Uni-Prot P67802) was identical in all species, except Loboptera decipiens, Symploce pallens and Blattella germanica (sequences not elucidated). Bold accession numbers correspond to sequences identified in this study. For Drosophila melanogaster peptides, the Gene-bank accession numbers are given.

Format: DOC Size: 112KB Download file

This file can be viewed with: Microsoft Word ViewerOpen Data

(b) Mass spectrometry

The dissection of the neurohaemal organs (abdominal perisympathetic organs, corpora cardiaca) as well as the sample preparation for MALDI-TOF MS (matrix-assisted laser desorption ionization time-of-flight mass spectrometry) and ESI-QTOF MS (electrospray ionization time-of-flight mass spectrometry) were performed as previously described [11,57]. MALDI-TOF MS: Mass spectra were obtained using an ABI 4700 proteomics analyzer (Applied Biosystems, Framingham, MA). To determine the sequences of the peptides, tandem MS experiments with a CID (collision induced dissociation) acceleration of 1 kV were performed. An unambiguous assignment of internal Leu/Ile was achieved by means of CID under high gas pressure that revealed unique patterns for the side chains of Leu and Ile (see [64]). Samples with CAPA peptides that contained Lys/Gln ambiguities were analysed again after dissolving the respective abdominal PSO preparations in acetic anhydride (2:1 methanol/acetic anhydride) which results in rapid acetylation of the ε-amino group of Lys. ESI-QTOF MS: In a few cases, data obtained from MALDI-TOF MS did not contain sufficient information to reveal the complete sequences of CAPA peptides. To fill the respective sequence gaps, nanoelectrospray mass spectra were acquired in the positive-ion mode using the API Qstar Pulsar (Applied Biosystems, Applera Deutschland GmbH, Darmstadt, Germany) fitted with a Protana (Odense, Denmark) nanoelectrospray source. Samples were purified using a homemade spin column and analyzed as described in Predel et al. [57].

(c) Sequence alignments and phylogenetic analysis

Homologous peptides were aligned using the Clustal × program package separately (parameter setting: gap penalty = 1; Protein Weight Matrix = BLOSUM), in contrast to aligning the whole data set simultaneously. There was no variability in the alignment results. Assignment of homologous gene products was facilitated due to their storage in specific neurohaemal organs and very similar C-terminal sequences. Phylogenetic analyses of peptides were performed under maximum-parsimony (MP) and Bayesian inference (BI) using PAUP4.0b10 [65] and MrBayes 3.1.2 [66], respectively. In the MP analysis, the heuristic search option with the tree-bisection-reconnection (TBR) branch swapping and 100 stepwise random additions of taxa was used. Gaps corresponding to missing data of few peptides were treated as missing characters, all other gaps as 21st amino acid. Levels of branch support were assessed using bootstrap resampling [67] with 1000 replicates to evaluate the reliability of the inferred topology. In the MP analysis, we tested the different data sets, i.e. CAPA peptides, adipokinetic hormone and sulfakinins, both separately and simultaneously following the total evidence approach. Because the topology of trees was similar (results not shown), we only present the results for our main data set (CAPA peptides) and overall data set. We tested the consistency by calculating the consistency index (CI), retention index (RI), and homoplasy index (HI) (see Figure 3 and 4).

For BI, we analysed the CAPA peptides and complete data set separately by using the fixed rates model test as default in MrBayes. Model free analysis of the peptide data set, however, did not change the topology of the trees (results not shown). A Markov Chain Monte Carlo (MCMC) sampling was run for 1 × 106 generations and trees were saved every 100 generations (with the first 1000 trees being discarded as "burn-in"). Gaps and missing characters were treated as missing data. Posterior probabilities with values greater than 49% are presented.

Authors' contributions

The strategy of the paper was mainly developed and coordinated by RP and to some degree by SR. RP and SR have written the manuscript. SR, BF and RP carried out insect dissection, sample preparation, and mass spectrometry; RP was responsible for species identification. BF and SR generated the phylogenetic analysis. Parts of the present study are incorporated within BF's diploma thesis. GG participated in the design and coordination of the study and helped to draft and improve the manuscript. All authors read and approved the final manuscript.

Acknowledgements

We acknowledge the assistance of Stefan Richter (Rostock, Germany), Olaf Bininda-Emonds (Jena, Germany), Dorit Liebers (Stralsund, Germany), Annett Kocum (Hiddensee, Germany), and Christian Derst (Berlin, Germany) during the phylogenetic analyses. We thank Christian Fischer (Göttingen, Germany) for helpful comments on an earlier version of the manuscript, and Susanne Neupert (Jena, Germany) as well as William K. Russell (College Station, Texas) for assistance in mass spectrometric analyses. We also thank Horst Bohn (Munich, Germany) for identification of some species of Blattellidae and Blattidae. The authors are grateful to S. Kambhampati (Manhattan, Kansas, USA), Yung Chul Park (Seoul, South Korea), Rüdiger Plarre (Berlin, Germany), Andreas Brune (Marburg, Germany), and Roland Dusi (frunol delicia, Delitzsch, Germany) for supporting the study with living or frozen specimens (see Appendix 1). We are grateful to Cathy Jenks (Bergen, Norway) for correcting the English language. The work was supported by grants from the Deutsche Forschungsgemeinschaft (Predel 595/6-1,2,4), the NRF (Pretoria, RSA; gun number 2053806 and FA2007021300002), and the Rotary Club Rondebosch (RSA).

References

  1. Rodríguez de la Vega RC, Possani LD: Overview of scorpion toxins specific for Na+ channels and related peptides: biodiversity, structure-function relationships and evolution.

    Toxicon 2005, 46:831-44. PubMed Abstract | Publisher Full Text OpenURL

  2. Olivera BM: Conus peptides: biodiversity-based discovery and exogenomics.

    J Biol Chem 2006, 281:31173-7. PubMed Abstract | Publisher Full Text OpenURL

  3. Claydon MA, Davey SN, Edwards-Jones V, Gordon DB: The rapid identification of intact microorganisms using mass spectrometry.

    Nat Biotechnol 1996, 14:1584-6. PubMed Abstract | Publisher Full Text OpenURL

  4. Strand FL: Neuropeptides: Regulators of Physiological Processes. Cambridge MA: The MIT Press; 1999. OpenURL

  5. Gäde G: The hypertrehalosaemic peptides of cockroaches: a phylogenetic study.

    Gen Comp Endocrinol 1989, 75:287-300. PubMed Abstract | Publisher Full Text OpenURL

  6. Gäde G, Grandcolas P, Kellner R: Structural data on hypertrehalosaemic neuropeptides from Cryptocercus punctulatus and Therea petiveriana: how do they fit into the phylogeny of cockroaches?

    Proc R Soc Lond B 1997, 264:763-768. Publisher Full Text OpenURL

  7. Gäde G, Marco HG, Desutter-Grandcolas L: A phylogenetic analysis of the adipokinetic neuropeptides of Ensifera.

    Physiol Entomology 2003, 28:283-289. Publisher Full Text OpenURL

  8. Bellés X, Graham LA, Bendena WG, Ding Q, Edwards JP, Weaver RJ, Tobe SS: The molecular evolution of the allatostatin precursor in cockroaches.

    Peptides 1999, 20:11-22. PubMed Abstract | Publisher Full Text OpenURL

  9. Aebersold R, Mann M: Mass spectrometry-based proteomics.

    Nature 2003, 422:198-207. PubMed Abstract | Publisher Full Text OpenURL

  10. Predel R, Kellner R, Baggerman G, Steinmetzer T, Schoofs L: Identification of novel periviscerokinins from single neurohaemal release sites in insects – MS/MS fragmentation complemented by Edman degradation.

    Eur J Biochem 2000, 267:3869-3873. PubMed Abstract | Publisher Full Text OpenURL

  11. Predel R: Peptidergic neurohemal system of an insect: Mass spectrometric morphology.

    J Comp Neurol 2001, 436:363-375. PubMed Abstract | Publisher Full Text OpenURL

  12. Clynen E, Baggerman G, Veelaert D, Cerstiaens A, Horst D, Harthoorn L, Derua R, Waelkens E, De Loof A, Schoofs L: Peptidomics of the pars intercerebralis-corpus cardiacum complex of the migratory locust, Locusta migratoria.

    Eur J Biochem 2001, 268:1929-1939. PubMed Abstract | Publisher Full Text OpenURL

  13. Neupert S, Johard HAD, Nassel DR, Predel R: Single-cell peptidomics of Drosophila melanogaster neurons identified by Gal4-driven fluorescence.

    Anal Chem 2007, 79:3690-3694. PubMed Abstract | Publisher Full Text OpenURL

  14. Predel R, Linde D, Rapus J, Vettermann S, Penzlin H: Periviscerokinin: a novel myotropic neuropeptide from the perisympathetic organs of the American cockroach.

    Peptides 1995, 16:61-66. PubMed Abstract | Publisher Full Text OpenURL

  15. Predel R, Rapus J, Holman GM, Nachman RJ, Wang Y, Penzlin H: Isolation of periviscerokinin-2 from the abdominal perisympathetic organs of the American cockroach.

    Peptides 1998, 19:801-809. PubMed Abstract | Publisher Full Text OpenURL

  16. Predel R, Kellner R, Nachman RJ, Rapus J, Gäde G: Differential distribution of pyrokinin-isoforms in cerebral and abdominal neurohemal organs of the American cockroach.

    Insect Biochem Molec Biol 1999, 29:139-144. Publisher Full Text OpenURL

  17. Kean L, Cazenave W, Costes L, Broderick KE, Graham S, Pollock VP, Davies SA, Veenstra JA, Dow JAT: Two nitridergic peptides are encoded by the gene capability in Drosophila melanogaster.

    Am J Physiol Regul Integr Comp Physiol 2002, 282:R1297-1307. PubMed Abstract | Publisher Full Text OpenURL

  18. Riehle MA, Garczynski SF, Crim JW, Hill CA, Brown MR: Neuropeptides and peptide hormones in Anopheles gambiae.

    Science 2002, 298:172-175. PubMed Abstract | Publisher Full Text OpenURL

  19. Hummon AB, Richmond TA, Verleyen P, Baggerman G, Huybrechts J, Ewing MA, Vierstraete E, Rodriguez-Zas SL, Schoofs L, Robinson GE, Sweedler JV: From the genome to the proteome: uncovering peptides in the Apis brain.

    Science 2006, 314:578-579. PubMed Abstract | Publisher Full Text OpenURL

  20. Li B, Predel R, Neupert S, Hauser F, Tanaka Y, Cazzamali G, Williamson M, Verleyen P, Schoofs L, Schachtner J, Grimmelikhuijzen CJP, Park Y: Genomics, Transcriptomics, and Peptidomics of Neuropeptides and Protein Hormones in the Red Flour Beetle Tribolium castaneum.

    Genome Res 2008, 18:113-122. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  21. Iversen A, Cazzamali G, Williamson M, Hauser F, Grimmelikhuijzen CJP: Molecular cloning and functional expression of a Drosophila receptor for the neuropeptides CAPA-1 and -2.

    Biochem Biophys Res Commun 2002, 299:628-633. PubMed Abstract | Publisher Full Text OpenURL

  22. Cazzamali G, Torp M, Hauser F, Williamson M, Grimmelikhuijzen C: The Drosophila gene CG9918 codes for a pyrokinin-1 receptor.

    Biochem Biophys Res Commun 2005, 335:14-19. PubMed Abstract | Publisher Full Text OpenURL

  23. Predel R, Gäde G: Peptidomics of neurohemal organs from species of the cockroach family Blattidae: how do neuropeptides of closely related species differ?

    Peptides 2005, 26:3-9. PubMed Abstract | Publisher Full Text OpenURL

  24. Kristensen NP: The phylogeny of hexapod 'orders'. A critical review of recent accounts.

    Z Zool Syst Evol 1975, 13:1-44. OpenURL

  25. Grimaldi D, Engel M: Evolution of the Insects. Cambridge University Press; 2005. OpenURL

  26. Thorne BL, Carpenter JM: Phylogeny of the Dictyoptera.

    Syst Entomol 1992, 17:253-268. Publisher Full Text OpenURL

  27. Kristensen NP: Forty years' insect phylogenetic systematics.

    Zool Beitr 1995, 36:83-124. OpenURL

  28. Bandi C, Sironi M, Damiani G, Magrassi L, Nalepa CA, Laudani U, Sacchi L: The establishment of intracellular symbiosis in an ancestor of cockroaches and termites.

    Proc R Soc Lond B 1995, 259:293-299. Publisher Full Text OpenURL

  29. Kambhampati S: A phylogeny of cockroaches and related insects based on DNA sequence of mitochondrial ribosomal RNA genes.

    Proc Natl Acad Sci USA 1995, 92:2017-2020. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  30. Kambhampati S: Phylogenetic Relationship among. Cockroach Families Inferred from Mitochondrial 12S rRNA Gene Sequence.

    Syst Entomol 1996, 21:89-98. OpenURL

  31. Maekawa K, Lo N, Rose HA, Matsumoto T: The evolution of soil-burrowing cockroaches (Blattaria: Blaberidae) from wood-burrowing ancestors following an invasion of the latter from Asia into Australia.

    Proc R Soc Lond B 2003, 270:1301-1307. Publisher Full Text OpenURL

  32. Grandcolas P: The phylogeny of cockroach families: A cladistic appraisal of morpho-anatomical data.

    Can J Zool 1996, 74:508-527. Publisher Full Text OpenURL

  33. Grandcolas P, D'Haese C: The phylogeny of cockroach families: Is the current molecular hypothesis robust?

    Cladistics 2001, 17:48-55. Publisher Full Text OpenURL

  34. Kjer KM: Aligned 18S and insect phylogeny.

    Syst Biol 2004, 53:506-514. PubMed Abstract | Publisher Full Text OpenURL

  35. Klass KD: The external male genitalia and the phylogeny of Blattaria and Mantodea.

    Bonner Zool Monogr 1997, 42:1-341. OpenURL

  36. Klass KD: The Proventriculus of the Dicondylia, with Comments on Evolution and Phylogeny in Dictyoptera and Odonata (Insecta).

    Zool Anz 1998, 237:15-42. OpenURL

  37. Nalepa CA, Lenz M: The ootheca of Mastotermes darwiniensis Froggatt (Isoptera: Mastotermitidae): homology with cockroach oothecae.

    Proc R Soc Lond B 2000, 267:1809-1813. Publisher Full Text OpenURL

  38. Maekawa K, Matsumoto T: Molecular phylogeny of cockroaches (Blattaria) based on mitochondrial COII gene sequences.

    Syst Entomol 2000, 25:511-519. Publisher Full Text OpenURL

  39. Grandcolas P: Phylogenetic systematics of the subfamily polyphaginae, with the assignment of Cryptocercus Scudder, 1862 to this taxon (Blattaria, Blaberoidea, Polyphagidae).

    Syst Entomol 1994, 15:145-158. OpenURL

  40. Grandcolas P: Reconstructing the Past of Cryptocercus (Blattaria: Polyphagidae): Phylogenetic Histories and Stories.

    Ann Entomol Soc Amer 1999, 92:303-307. OpenURL

  41. Lo N, Tokuda G, Watanabe H, Rose H, Slaytor M, Maekawa K, Bandi C, Noda H: Evidence from multiple gene sequences indicates that termites evolved from wood-feeding cockroaches.

    Current Biology 2000, 10:801-804. PubMed Abstract | Publisher Full Text OpenURL

  42. Lo N, Bandi C, Watanabe H, Nalepa C, Beninati T: Evidence for cocladogenesis between diverse dictyopteran lineages and their intracellular endosymbionts.

    Mol Bio Evol 2003, 20:907-913. Publisher Full Text OpenURL

  43. Maekawa K, Kitade O, Matsumoto T: Molecular phylogeny of orthopteroid insects based on the mitochondrial cytochrome oxidase II gene.

    Zool Sci 1999, 16:175-184. Publisher Full Text OpenURL

  44. Inward D, Beccaloni G, Eggleton P: Death of an order: a comprehensive molecular phylogenetic study confirms that termites are eusocial cockroaches.

    Biol Lett 2007, 3:331-335. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  45. Cleveland LR, Hall SR, Sanders EP, Collier J: The wood-feeding roach Cryptocercus, its protozoa, and the symbiosis between protozoa and roach.

    Mem Am Acad Arts Sci 1934, 17:185-342. OpenURL

  46. Deitz LL, Nalepa CD, Klass KD: Phylogeny of the Dictyoptera Re-examined (Insecta).

    Entomol Abh 2003, 61:69-91. OpenURL

  47. Klass KD, Meier R: A phylogenetic analysis of Dictyoptera (Insecta) based on morphological characters.

    Entomol Abh 2006, 63:3-50. OpenURL

  48. Predel R, Wegener C: Biology of the CAPA peptides in insects.

    Cell Mol Life Sci 2006, 63:2477-2490. PubMed Abstract | Publisher Full Text OpenURL

  49. Grandcolas P: The monophyly of the subfamily Perisphaeriinae (Dictyoptera: Blattaria: Blaberidae).

    Syst Entomol 1997, 22:123-130. Publisher Full Text OpenURL

  50. Gäde G, Hoffmann KH, Spring JH: Hormonal Regulation in Insects: Facts, Gaps, and Future Directions.

    Physiol Rev 1997, 77:963-1032. PubMed Abstract | Publisher Full Text OpenURL

  51. Grandcolas P, Deleporte P, Desutter-Grandcolas L, Daugeron C: Phylogenetics and Ecology: as many characters as possible should be included in the cladistic analysis.

    Cladistics 2001, 17:104-110. OpenURL

  52. Fricker LD, Lim JY, Pan H, Che FY: Peptidomics: identification and quantification of endogenous peptides in neuroendocrine tissues.

    Mass Spectrom Rev 2006, 25:327-344. PubMed Abstract | Publisher Full Text OpenURL

  53. Predel R, Eckert M, Pollak E, Molnar L, Scheibner O, Neupert S: Peptidomics of identified neurons demonstrates a highly differentiated expression pattern of FXPRLamides in the neuroendocrine system of an insect.

    J Comp Neurol 2007, 500:498-512. PubMed Abstract | Publisher Full Text OpenURL

  54. Audsley N, Weaver RJ: A comparison of the neuropeptides from the retrocerebral complex of adult male and female Manduca sexta using MALDI-TOF mass spectrometry.

    Regul Pept 2003, 116:127-137. PubMed Abstract | Publisher Full Text OpenURL

  55. Predel R, Roth S: Neuropeptide Evolution and the Analysis of Phylogenetic Relationships in Blattaria (Hexapoda).

    Arthropod Syst Phylogeny 2007, 65:3-6. OpenURL

  56. Roth LM: The male genitalia of Blattaria. XI. Perisphaeriinae.

    Psyche 1973, 80:305-348. Publisher Full Text OpenURL

  57. Predel R, Neupert S, Wicher D, Gundel M, Roth S, Derst C: Unique accumulation of neuropeptides in an insect: FMRFamide-related peptides in the cockroach, Periplaneta americana.

    Eur J Neurosc 2004, 20:1499-1513. Publisher Full Text OpenURL

  58. Rahman MM, Fromm B, Neupert S, Kreusch S, Predel R: Extended FMRFamides in dipteran insects: conservative expression in the neuroendocrine system is accompanied by rapid sequence evolution.

    Gen Comp Endocrinol 2008. OpenURL

  59. Kastin AJ: Handbook of Biologically Active Peptides. Amsterdam: Academic Press/Elsevier Science & Technology; 2006. OpenURL

  60. Predel R, Gäde G: Identification of the abundant neuropeptide from abdominal perisympathetic organs of locusts.

    Peptides 2002, 23:621-627. PubMed Abstract | Publisher Full Text OpenURL

  61. Clynen E, Huybrechts J, De Loof A, Schoofs L: Mass spectrometric analysis of the perisympathetic organs in locusts: identification of novel periviscerokinins.

    Biochem Biophys Res Commun 2003, 300:422-428. PubMed Abstract | Publisher Full Text OpenURL

  62. Predel R, Roth S, Neupert S, Picker M: New insect order Mantophasmatodea: species differentiation by mass fingerprints of peptide hormones?

    J Zool Syst Evol Res 2005, 43:149-156. Publisher Full Text OpenURL

  63. Roth S, Köhler G, Reinhardt K, Predel R: A discrete neuropeptide difference between two hybridizing grasshopper subspecies.

    Biol J Linn Soc 2007, 91:541-548. Publisher Full Text OpenURL

  64. Nachman RJ, Russell WK, Coast GM, Russell DH, Predel R: Mass spectrometric assignment of Leu/Ile in neuropeptides from single neurohemal organ preparations of insects.

    Peptides 2005, 26:2151-2156. PubMed Abstract | Publisher Full Text OpenURL

  65. Swofford DL: PAUP*: Phylogentic Analysis Using Parsimony (* and other methods). Vers. 4. Sunderland, MA: Sinauer Associates Press; 2002. OpenURL

  66. Ronquist F, Huelsenbeck JP: MrBAYES 3: Bayesian phylogenetic inference under mixed models.

    Bioinformatics 2003, 19:1572-1574. PubMed Abstract | Publisher Full Text OpenURL

  67. Felsenstein J: Confidence Limits on Phylogenies: An Approach Using the Bootstrap.

    Evolution 1985, 39:783-791. Publisher Full Text OpenURL