Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Stridulations Reveal Cryptic Speciation in Neotropical Sympatric Ants

Abstract

The taxonomic challenge posed by cryptic species underlines the importance of using multiple criteria in species delimitation. In the current paper we tested the use of acoustic analysis as a tool to assess the real diversity in a cryptic species complex of Neotropical ants. In order to understand the potential of acoustics and to improve consistency in the conclusions by comparing different approaches, phylogenetic relationships of all the morphs considered were assessed by the analysis of a fragment of the mitochondrial DNA cytochrome b. We observed that each of the cryptic morph studied presents a morphologically distinct stridulatory organ and that all sympatric morphs produce distinctive stridulations. This is the first evidence of such a degree of specialization in the acoustic organ and signals in ants, which suggests that stridulations may be among the cues used by these ants during inter-specific interactions. Mitochondrial DNA variation corroborated the acoustic differences observed, confirming acoustics as a helpful tool to determine cryptic species in this group of ants, and possibly in stridulating ants in general. Congruent morphological, acoustic and genetic results constitute sufficient evidence to propose each morph studied here as a valid new species, suggesting that P. apicalis is a complex of at least 6 to 9 species, even if they present different levels of divergence. Finally, our results highlight that ant stridulations may be much more informative than hitherto thought, as much for ant communication as for integrative taxonomists.

Introduction

The Tropics are home to nearly two-thirds of the World's known biodiversity, but also to a large amount of species that have remained unnoticed [1], partly due to the occurrence of cryptic species, i.e. of two or more distinct species that are erroneously classified (and hidden) under a single species name, due to their very similar morphology [2]. Such species cause great problems for taxonomists as they cannot be readily or reliably distinguished only on a morphological basis, and the taxonomic challenge they pose underlines the importance of using multiple criteria in species delimitation. Indeed, concordant changes in several characteristics of an organism, and corroboration from independent data constitute better evidence for separating species [2], [3]. As accurate species identification is crucial both to research in all areas of biology and to biodiversity conservation [2], [4], there is an urgent need to coalesce effective tools which allow the clarification of these taxonomic issues in order to estimate the real biodiversity in cryptic species groups.

Acoustic signals may differ between cryptic species and were often the first clue of a hidden diversity in many insect groups, like hemipterans [5], lacewings [6], orthopterans [7][9] and flies [10], [11]. By delimitating species phenotypic variability, acoustic analyses have revealed many unsuspected species and solved several confusing taxonomic problems. Thus, acoustic descriptions should be systematically included in species diagnosis [12]. The use of DNA sequences (“DNA barcoding”) in extensive phylogenetic studies has also revealed an effectively high level of hidden biodiversity [13][16], suggesting that molecular data should be incorporated by taxonomists as a matter of routine [2].

In ants, stridulatory sound production is known since the late nineteenth century [17][21]. This faculty seems to have evolved several times independently and can be found in a great number of species from five subfamilies (Myrmicinae, Pseudomyrmecinae, Ponerinae, Ectatomminae, Nothomyrmecinae) [22], [23]. Stridulations are produced during dorso-ventral movements of the gaster by the rubbing of the distal border of the 3rd segment of the abdomen that acts like a scraper on the stridulatory file, made up of perfectly parallel and rectilinear tegument ridges on the mid-dorsal edge of the fourth abdominal segment [24], [25]. They consist in chirps made by train of pulses, in which each pulse corresponds to the rubbing movement of the scraper on one ridge on the stridulatory file [26]. These signals are essentially transmitted by the substrate [27], [28], but evidence for a perception of air-transmitted sounds, at least over distances of a few centimetres, is also available [29]. These signals are generally barely audible without amplification [22], [29] but in some Ponerinae species, like Pachycondyla apicalis, the intensity of the airborne sound emitted by a single ant at a distance of 1 cm can reach more than 93 dB [30]. The frequency of the signal can vary from a few hertz, like in Solenopsis richteri [29] and Myrmica spp [31], [32] to up to 84 kHz in P. apicalis [30]. Thus far, stridulations have been shown to be produced in several behavioural contexts depending on the species examined, like food recruitment [27], [33][34], trophallaxis [35], nest emigration [36], intra- and inter-specific conflicts [37], [38], and mating [39], [40]. Some ants can also respond to stridulations produced by their mutualists [41][44] or even their parasites [31][32].

However, even if stridulations are common events in ant societies, it is still probably the least understood mode of communication and detailed studies on the acoustic characteristics of these signals are scarce [45]. Furthermore, up to now, most work on ant stridulations refers only to the subfamily Myrmicinae, and almost nothing is known about primitive ants like the Ponerinae. In this ant subfamily, the genus Pachycondyla is one of the most ancient still living [26], [46]. It presents over 60 described species just for the Neotropics [47], and after some assessment studies [48][50], it appears to be a good bioindicator of the myrmecological diversity and quality of Neotropical ecosystems. Recent studies [51][53] demonstrated that many Neotropical Pachycondyla actually consist of cryptic species complexes, indicating that the diversity in this genus can be really underestimated.

The Pachycondyla apicalis species complex is a good example of this problem and its taxonomy remains unsettled. Ants in this group are large, conspicuous insects found in Neotropical forests from southern Mexico to Paraguay, overall presenting similar ecological and biological features [54], [55]. Furthermore, ants in this complex also share the same general morphology and coloration (the “morph apicalis”) which is probably under stabilizing selection as it is also copied by other arthropods like spiders [56][58]. These sympatric species could thus form a massive mimicry ring which would represent an advantage against predators [15], [57], [59]. Wild [54] revised the taxonomy of this group using several morphological and biometrical criteria and recognised three broadly sympatric species: P. apicalis [60], P. obscuricornis [61] and P. verenae [62] instead of only two as thought before [63]. Although he observed a considerable morphological variation for P. verenae and P. apicalis across their distribution ranges he did not consider further division of the complex. Later, in a morphological study combining cytogenetical and ecological data when available, Delabie et al. [59] demonstrated that, given the stability of differences through multiple cases of sympatry, the variability pointed by Wild [54] could in fact refer to a species mosaic rather than a geographic cline. They recognized seven distinct taxa within this group: four for P. apicalis, two for P. verenae and one for P. obscuricornis, and despite supporting a really intricate diversity inside the P. apicalis species complex, the authors were not able to conclude about the validity and the taxonomic status (e.g. species, subspecies, ecotypes) of the different morphs within each species without further investigations.

Here for the first time, we test acoustics as a tool to assess the real diversity in a cryptic species complex of Neotropical ants. We measured the overall structure of the stridulatory file and performed a detailed acoustic analysis of the stridulations produced in five morphs of the P. apicalis species complex. To improve the consistency in conclusions by comparing different approaches, we also analysed a fragment of the mitochondrial DNA cytochrome b for all morphs considered.

Methods

Ants

Colonies of the P. apicalis species complex were collected at Petit Saut, French Guiana (n = 10) and Los Tuxlas, Mexico (n = 1). Colonies were reared in the laboratory in artificial plastered nests. The nests were maintained at 25±1°C, with approximately 65±10% relative humidity, and a 12L∶12D photoperiod. All colonies were provided with an identical diet (honey/apple mixture and crickets) twice a week. Ant collection, husbandry and experimental procedures used in this study fulfilled all the legal requirements concerning insect experimentation of France.

Ants were classified into morphs within the currently named species according to Delabie et al. [59] classification. Three different morphs could be identified, one for P. verenae (PVE) and two for P. apicalis (PAP), respectively: PVE Morph 1 (5 colonies), PAP Morph 3 (1 colony) and PAP Morph 4 (2 colonies) (see supplementary Figures S1, S2 and S3). Among our collected colonies, we did not find the P. verenae Morph 2 or P. apicalis Morphs 1 or 2. Moreover, three colonies of P. apicalis from Petit Saut did not fit any described morph in Delabie et al. [59]. We treated them here as 2 new morphs: PAP Morph 5 (2 colonies) and PAP Morph 6 (1 colony), as they also presented subtle distinctive morphological traits. PAP Morph 5 is moderate sized, with an emarginated petiole and a very finely striated cuticle and hairy head and body (see supplementary Figure S4). PAP Morph 6 presents a rounder petiole and is clearly bigger than all other P. apicalis examined (see supplementary Figure S5). PVE Morph 1 and PAP Morphs 4, 5, and 6 were found in sympatry in the same site at Petit Saut, French Guyana. PAP Morph 3 is allopatric to all other morphs and occurs only in Mexico [59]. Vouchers of each morph were deposited in the CPDC collection of the Laboratório de Mirmecologia, Cocoa Research Center at Itabuna (Bahia, Brazil).

Morphometric Study of the Stridulatory Organ

A total of 41 workers from 8 colonies (Table 1) were dissected and the segments of the gaster containing the stridulatory file were cleaned in an ultrasonic-wave bath, air-dried and then placed on aluminium stubs. The samples were coated with a mixture of 80% gold/20% palladium and examined with a Leica Stereoscan 440 scanning electron microscope (SEM). Measurements of the stridulatory organ were obtained from the digitalized SEM images. Eight variables were measured from each stridulatory file: Length, maximal width, 1st, 2nd and 3rd quartile widths, number of ridges and inter-ridge distance in the medial and distal portions of the stridulatory file (Figure 1). For these two last variables, 5 measures were taken from each worker and the mean value was computed. Also, as an estimate of the ants' size, the thorax length of 42 other workers was measured with a Zeiss Stereo Microscope at a magnification of 10x.

thumbnail
Figure 1. Scanning electron micrographs of the stridulatory file of a Pachycondyla apicalis species complex worker.

(A) General view of the stridulatory file, (B) detail of the ridges showing the interface between the fine structure of the file and the surrounding cuticle and (C) detail in the medial portion of the stridulatory file, showing the inter-ridge distances.

https://doi.org/10.1371/journal.pone.0015363.g001

thumbnail
Table 1. Morphometric characteristics (mean±SE) of the stridulatory file in five morphs of the Pachycondyla apicalis species complex.

https://doi.org/10.1371/journal.pone.0015363.t001

Stridulation recording and analysis

Ants from the P. apicalis species complex produce stridulations which result in both airborne sound and substrate vibrations. Airborne sound presents audible and ultrasound components, and in the following we refer to the audible component only. A total of 40 workers from 5 colonies (Table 2) were recorded. All recordings were carried out in a low-noise room where the ambient temperature was kept at 25±1°C and the relative humidity at 65±10%. The recording setup consisted of an omnidirectional Sennheiser K6 microphone (frequency response: 30 to 20000 Hz±1 dB) connected to a Marantz PMD 671 digital recorder with sampling frequency at 48 kHz. We did not consider frequencies superior to 20 KHz, due to technical limitations of the microphone. Ants were held with forceps 1 cm from the microphone during recording. The following temporal parameters were analysed using the software Avisoft-SASLab Pro, version 4.40 [64]: the chirp duration, the inter-chirp interval, and for each chirp we measured the number of pulses, the pulse repetition rate, the mean inter-pulse interval as well as the inter-pulse interval in the 1st, 2nd and 3rd thirds of the chirp. The frequency parameters considered for each chirp were: the dominant frequency, the frequencies at 25, 50 and 75% of the signal energy, and the percentage of energy below 14 kHz. The maximal and minimal intra-pulse frequencies were also calculated for each individual, by the zero-crossing method [65]. For each chirp analysed, we calculated the maximal and minimal intra-pulse frequencies for 10 pulses and the mean value was computed (Figure 2). A series of ten chirps was analyzed for each ant, and the mean value was computed.

thumbnail
Figure 2. Stridulatory signal of a Pachycondyla apicalis species complex worker.

(A) Oscillogram of a series of chirps, showing the chirp duration (CD), the inter-chirp interval (ICI), the number of pulses (NP) in a chirp and the inter-pulse interval (IPI). (B) Spectrogram of a chirp, showing the dominant frequency (DM) and the quartiles of frequencies at 25, 50 and 75% of the signal energy (Q25, Q50 and Q75 respectively).

https://doi.org/10.1371/journal.pone.0015363.g002

thumbnail
Table 2. Acoustic characteristics (mean±SE) of the stridulatory signals in five morphs of the Pachycondyla apicalis species complex.

https://doi.org/10.1371/journal.pone.0015363.t002

DNA extraction, amplification and sequencing

The DNA of 22 workers from 11 colonies was extracted from ethanol-preserved tissues (head and thorax) using a DNeasy Blood & Tissue kit (Qiagen Inc., Valencia, CA) following the manufacturer's protocol. Mitochondrial DNA variation was assayed by the amplification of a portion of the mtDNA cytochrome b (cyt b, ∼700 bp) using primers CB1 (5′-TAT-GTA-CTA-CCA-TGA-GGA-CAA-ATA-TC-3′) and tRS (5′-TAT-TTC-TTT-ATT-ATG-TTT-TCA-AAA-C-3′) from Simon et al. [66]. Each PCR was carried out in a 50-µL volume according to a standard protocol using a T1 thermal cycler (Biometra). The thermal cycle profile was as follows: 2 min at 94°C; 35 cycles at [30 s at 94°C/60 s at 50°C/60 s at 72°C]; 5 min at 72°C. Amplified products were sequenced using the same primers as used for the amplification by Genoscreen (Lille, France) using an ABI 3730XL automatic sequencer (Applied Biosystems).

Sequences analyses were edited and aligned using the default settings of Clustal X [67] and checked by eye. To generate phylogenetic trees, we used pairwise distances (Neighbor Joining algorithm, NJ). As Pachycondyla is a paraphyletic genus (C. Schmidt, pers. com.), we used two other divergent Pachycondyla species as outgroups: P. villosa and P. goeldi, respectively. Average intermorph genetic divergence was calculated using the Kimura 2-parameter model [68] using the MEGA4 program [69]. Other models (Jukes–Cantor or Tamura-Nei distance) when applied, resulted in similar results. NJ tree was constructed with the MEGA4 program. The robustness of the tree was tested with 1000 bootstrap replications. The equality of evolutionary rate between all sequences was tested using Tajima' relative rate test [70] in MEGA4.

Statistical analysis

Workers used for the morphometric measurements came from two different colonies for PVE Morph 1 and PAP Morphs 4 and 5. As these colonies did not differ for the structure of the stridulatory organ (MANOVA “colonies nested within morphs”, Wilk's λ = 0.249, F24,47 = 1.205, p = 0,2863), their data were pooled for further comparisons between morphs. Discriminant function analyses (DFA) were performed to identify potential differences of the stridulatory file morphometry and the acoustic characteristics of stridulations between morphs and to calculate the success rate of individual discrimination (correct classification rate) using Mahalanobis distances. As the temporal parameter of stridulations “inter-chirp interval” showed a very high coefficient of variation within morphs (ranging from 13.31 to 45.50%), this parameter was not included in the DFA. Additionally, we compared each morphometric and acoustic parameter using one-way ANOVA followed by Unequal N HSD (Honest Significance Difference) post hoc tests [71] to understand how these parameters varied between morphs. Paired Student's t-tests were used to establish the differences between the medial and distal inter-ridge distances within groups. When necessary, Box-Cox transformation was applied to achieve normality on some parameters [72]. All results are stated as mean±SE. The significance level was taken at α≤0.05 to assess differences. All analyses were conducted using Statistica v8.0 [73].

Results

Morphometry of the stridulatory organ

The discriminant function analysis of all morphometric parameters considered for the stridulatory file clearly separates all studied morphs (Wilk's λ = 0.005, F32,108 = 10.712, p<0.001; Figure 3, Table 1). The file (see supplementary Figure S6) can vary in length (Figures 4 and 5.a), maximal and quartile widths (Figure 4), medial and distal inter-ridge distances (Figure 6) and number of ridges between morphs (Figure 5.b). These differences are combined in a distinctive way within each morph, which allows a 97.5% correct classification rate of all individuals based on the overall pattern of the stridulatory file.

thumbnail
Figure 3. Discriminant function analysis of the stridulatory file morphometry of the P. apicalis species complex.

PVE: P. verenae, PAP: P. apicalis. Ellipses are 95% confidence intervals around centroids.

https://doi.org/10.1371/journal.pone.0015363.g003

thumbnail
Figure 4. Pattern of the stridulatory file of five morphs from the P. apicalis species complex.

First (Q1), second (Q2) and third (Q3) quartile widths of the file of ridges in relation with the length (L). ANOVAQ1, F4,36 = 18.74, ANOVAQ2, F4,36 = 20.25, ANOVAQ3, F4,36 = 25.88 and ANOVALength, F4,36 = 9.82. All p<0.001. PVE: P. verenae, PAP: P. apicalis.

https://doi.org/10.1371/journal.pone.0015363.g004

thumbnail
Figure 5. Ant size and stridulatory file features.

Relationship between the stridulatory file length and (A) the ant size (represented by thorax length) and (B) the number of rigdes of five morphs from the P. apicalis species complex. (ANOVAThorax, F4,37 = 57.82, p<0.001, ANOVAFile, F4,36 = 9.83, p<0.001 and ANOVARidges, F4,36 = 9.66, p<0.001). When different, capital letters indicate significant differences for parameters on the X axes, and small letters indicate significant differences b for parameters on the Y axes tests, p<0.05. PVE: P. verenae, PAP: P. apicalis.

https://doi.org/10.1371/journal.pone.0015363.g005

thumbnail
Figure 6. Inter-ridge distances of the stridulatory file in the P. apicalis species complex.

Comparison of the inter-ridge distances in the medial and distal portions of the stridulatory file within morphs (Paired Student's t-tests within morphs, tPVE Morph 1 = 4.28, tPAP Morph 3 = −3.11, tPAP Morph 4 = 3.86, tPAP Morph 5 = −2.93, tPAP Morph 6 = 7.18; * p<0.05 and ** p<0.01) and between morphs (ANOVAMedial, F4,36 = 15.17, p<0.001, ANOVADistal, F4,36 = 21.61, p<0.001. When different, small letters indicate significant differences of inter-ridge distances on the medial portion of the stridulatory file and capital letters indicate significant differences of inter-ridge distances on the distal portion of the stridulatory file (Unequal N HSD tests, p<0.05)). PVE: P. verenae, PAP: P. apicalis.

https://doi.org/10.1371/journal.pone.0015363.g006

Differences observed in the pattern of the stridulatory file are not due only to allometric differences between individuals. Morphs with different body size present the same length for the stridulatory file (Figure 5.a). Moreover, stridulatory files of the same length can comprise different number of ridges (Figure 5.b), which directly affect the inter-ridges distance. This latter parameter presents complex variation patterns along the stridulatory file depending on the morph. It can vary between the medial (ANOVA, F4,36 = 15.17, p<0.001) and distal (ANOVA, F4,36 = 21.61, p<0.001) regions of the file and the pattern is not the same for all morphs: PVE Morph 1 and PAP Morphs 4 and 6 present ridges more spaced in the middle than in the distal region of the file, while PAP Morph 3 and 5 present the opposite pattern (Figure 6).

Acoustic analyses of stridulations

The DFA of all temporal and frequency parameters (except the inter-chirp interval) reveals that each sympatric morph from French Guiana produce distinctive sounds (Wilk's λ: 0.01987, F56,87: 2.724, p<0,0001; Figure 7, see supplementary Audio S1, Audio S2, Audio S3, Audio S4 and Audio S5). Here, the differences in the acoustic characteristics of stridulations allowed a correct classification rate of 95% of all individuals. PVE Morph 1 produces shorter chirps than PAP Morphs 4 and 6 (ANOVAChirps, F4,35 = 5.24, p<0.01, Unequal N HSD, p<0.01 and p<0.05 respectively, Table 2), and with fewer pulses than PAP Morph 6 (ANOVANPulses, F4,35 = 4.18, p<0.01, Unequal N HSD, p<0.01, Table 2), notwithstanding the similarity in the number of ridges in the stridulatory file (Figure 5.b). This result is due to the smaller inter-ridge distances presented by PVE Morph 1 (Figure 6), but also to the smaller number of ridges that are rubbed during stridulatory movements (Table 2). In spite of presenting a high inter-individual variation, the inter-chirp interval in PAP Morph 6 is bigger than all the other PAP morphs (ANOVAICI, F4,35 = 4.84, p<0.01, Unequal N HSD, p<0.05 for PAP Morph 3 and p<0.01 for PAP Morphs 4 and 5, Table 2), which means this morph produces fewer stridulations in a given time. Nevertheless, this morph seems to rub the stridulatory file faster during stridulation production, as suggested by the tendencies for a higher pulse repetition rate and a lower mean inter-pulse interval together with larger inter-ridge distances (Figure 6, Tables 1 and 2).

thumbnail
Figure 7. Discriminant function analysis of the stridulations of the P. apicalis species complex.

All temporal and frequency parameters measured (except the inter-chirp interval) are considered in the model. PVE: P. verenae, PAP: P. apicalis. Ellipses are 95% confidence intervals around centroids.

https://doi.org/10.1371/journal.pone.0015363.g007

Concerning frequency parameters, PAP Morph 4 stridulations have higher dominant frequency than PVE Morph 1 and PAP Morph 5 (ANOVADF, F4,35 = 4.21, p<0.01, Unequal N HSD, p<0.05 for both morphs, Table 2), higher frequencies at 25% and 50% of the signal energy than PVE Morph 1 (ANOVAF25%, F4,35 = 2.82, p<0.05, Unequal N HSD, p<0.05 and ANOVAF50%, F4,35 = 4.05, p<0.01, Unequal N HSD, p<0.05, Table 2), and higher intra-pulse minimal frequency than all other morphs (ANOVAMFP, F4,35 = 8.17, p<0.01, Unequal N HSD, p<0.01 for all morphs). The intra-pulse maximal frequency is also significantly lower in PVE Morph 1 compared to PAP Morphs 4 and 6 (ANOVALFP, F4,35 = 5.62, p<0.01, Unequal N HSD, p<0.01 for PAP Morph 4 and p<0.05 for PAP Morph 6). Last, the stridulations of the allopatric PAP Morph 3 from Mexico did not significantly differ only from those of PAP Morph 5 from French Guiana (Figure 7) even if these morphs present distinctive stridulatory organs (Figures 3, 4, 5, 6).

Nucleotide composition and sequence variation

A total of 690 base pairs were analysed for cyt b. Accession numbers range from HM770106 to HM770124 in Genbank. No pseudogenes, no insertions, deletions nor any rearrangements were detected. As in other hymenopteran mitochondrial genome [74], there is an A-T bias in the base composition of cyt b: we obtained on average 33.9% of A, 44.1% of T, 14.7% of C, and 7.3% of G. Over the 690 base pairs analysed, the number of variable characters was 259 among which 156 were found to be parsimony-informative. The transition/transversion rate ratios are k1 = 1.207 (purines) and k2 = 2.974 (pyrimidines), with an overall transition/transversion bias of 0.584. NJ tree of haplotypes obtained for the different morphs is presented in Figure 8, rooted with P. villosa and P. goeldi as outgroups. Tajima's relative tests are significant (X2 = 4.25, p = 0.038 to X2 = 8, p = 0.005, df = 1) when morphs of the clade constituted by P. apicalis are compared to P. verenae or PAP Morph 5, rejecting as a consequence the molecular clock hypothesis for this set of sequences.

thumbnail
Figure 8. Neighbour-joining phylogenetic tree of mtDNA haplotypes of the Pachycondyla apicalis species complex.

The tree is rooted using P. villosa and P. goeldi sequences as outgroups 1 and 2, respectively. Numbers indicate bootsrap values. PVE: P. verenae, PAP: P. apicalis.

https://doi.org/10.1371/journal.pone.0015363.g008

We could observe 3 main groups, in which PAP Morph 5 clearly differs from all other morphs (about 15%), and presents a basal position to the P. apicalis and P. verenae groups (Figure 8).

Within P. verenae Morph 1, we obtain two clades but the level of divergence is low (mean distance of 1.5±0.3%). Within the group P. apicalis, the genetic distance between Morphs 3, 4 and 6 is considerable. When compared to the other morphs of P. apicalis, PAP Morph 3 varies in 10% of the base pairs considered, which nearly corresponds of the level of divergence between P. verenae and the clade P. apicalis. Between the sympatric P. apicalis Morphs 4 and Morph 6 the genetic variation is about 7% (Figure 8).

Discussion

Our results demonstrate that each studied cryptic morph of the P. apicalis complex presents a morphologically distinct stridulatory organ, and that all sympatric morphs make distinctive sounds. The differentiation observed for the stridulatory organs is not only due to allometric differences between individuals but also to intrinsic morph features. The distinct acoustic signals produced, in their turn, are the result of this morphological specificity together with ant behaviour for stridulation production. Indeed, inter-specific competitive interactions in sympatry may have led to divergent selection acting in contrasting directions between morphs. In contrast, the similarity observed between the acoustic signals produced by PAP Morphs 3 and 5 in spite of the distinctiveness of their stridulatory files can be due to the fact that the allopatric PAP Morph 3 is not subjected to the same selective pressure as the sympatric morphs.

To our knowledge, we provide the first record of such a degree of acoustic specialization in closely related ants, both at the level of the production organ and of the produced signal. For a long time, stridulations have been considered as a mere generalist alarm signal that has no selective advantage to be species-specific and that could rarely be found specific even at the genus level [22]. In the few groups of sympatric and/or related species studied to date [three sympatric species of Pogonomyrmex ssp [39], four sympatric Messor species [38], [75], four Neotropical Ectatomma [30] and five Myrmica spp. [31], [32] no prominent inter-specific differences could be observed in the structures of the stridulatory file and in the characteristics of the signal. The only specificity already demonstrated in ant stridulations was at the intra-specific level, between different castes: for major and minor workers in Atta cephalotes [76] and for gynes and males or workers in four sympatric Messor species [38]. Recently, Barbero et al [31], [32] demonstrated that acoustic signals carry information about the caste and the status of a colony member in Myrmica schencki and trigger distinct behavioural responses by workers as a function of the identity of the emitter.

In ants, the most important channel for communication involves chemical and, to a smaller degree, tactile cues [26], [77]. However, the clear differentiation and specificity in the stridulatory file and signals observed in this study in a group of species with considerable morphological stasis suggests that acoustic communication may have a more important role than generally thought during interspecific relationships in these ants. We can also expect that, in this group of cryptic ants, acoustic signals might modulate chemical and tactile cues in different ways [31], and that a synergy between distinct-source information might improve communication in different behavioural contexts.

The variations observed here for some acoustic parameters would have likely been decreased and the differences between morphs would have likely been built up by studying individuals coming from different colonies for each morph. Nevertheless, mitochondrial DNA variation estimated for cytochrome b well supports the acoustic differences observed for all five morphs and confirms the taxonomic potential of acoustics for this group of cryptic ants and possibly for other stridulating ant species. We thus found that each studied morph in the P. apicalis species presents a real genetic identity. The most surprising result is the basal position of PAP Morph 5 and its high level of genetic divergence compared to other sympatric P. apicalis. The genetic isolation of PAP Morph 5 is mirrored by some morphological characters that could have been neglected, leaving this morph hidden until now. Moreover, despite the non-divergence of PAP Morph 5 stridulations from those of the allopatric PAP Morph 3, the corroboration of the distinct morphometry of the stridulatory file, genetics and morphology are enough evidence to separate these two taxa. Similarly, the genetic distances observed between the other two morphs of P. apicalis, even if lower than between PAP Morphs 3 and 5, are all well above the usual inter-specific values found for some cryptic insects groups, i.e. the 3% sequence divergence threshold typically used in the barcoding studies [78].

The congruence of the genetic data with acoustic and morphological information leads us to propose each morph studied here as a valid new species. In fact, the subtle morphological variations observed for these ants confirm not to be random, and what Wild [54] previously thought to be only intra-specific morphological variation, is verified by acoustic analysis and genetics to be distinct inter-specific traits, as supposed by Delabie et al. [59]. Given that acoustics matched the mitochondrial DNA divergence in the P. apicalis species complex and showed a high potential for species diagnosis in this group of cryptic species, further studies should consider including this tool to investigate cryptic diversity in stridulating ants. Indeed, with all the specific acoustic and morphological traits evidenced by the present study, the cryptic species in the P. apicalis complex can be considered from now on as pseudo-cryptic species, i.e. cryptic species for which after detailed comparisons of morphological and non-morphological features key characters can be established for their identification [3], [79]. In addition, further work should try to compare the type specimens of the described species, seeking to verify which morphs are the true P. apicalis [60] and P. verenae [62], and even P. flavicornis [60], a related P. apicalis species described by Latreille at the same time as P. apicalis, but synonymised later by Brown [63].

Speciation is a very complex process which is affected by many different factors (genetic, ecological, developmental, environmental, etc.) interacting in nonlinear ways [80]. The integrative approach undertaken in this paper was thus essential for the recognition of the real biodiversity in our sample of morphs of the P. apicalis species complex. Yet, the biodiversity inside this complex could even be more important than shown here. In this study, we had access to only three of seven morphs described in Delabie et al. [59] (P. verenae Morph 1 and P. apicalis Morphs 3 and 4) and to two new sympatric morphs of P. apicalis that we identified for French Guiana (P. apicalis Morphs 5 and 6). When taking into consideration that within the morphs non studied here there are still the ‘rare’ P. obscuricornis, a second morph of P. verenae and two additional P. apicalis morphs [59], the diversity inside the P. apicalis complex is expected to reach nine cryptic species, and a complete survey over the whole distribution range of the complex might uncover many more species. Additional research on the biogeography, ecology and behaviour of these ants could also reveal species idiosyncrasies, and help us to better understand the speciation process within this Neotropical species complex. For example, one might uncover the selective forces that have driven this high diversification in the P. apicalis clade and not in the P. verenae, or elucidate if this group diverged in sympatry or if their actual distribution is derived. Finally, as evidenced in skipper butterflies [15], fig-pollinating wasps [81], cerambycid beetles [14], pseudoscorpions [13] and ants [51][53], our results add support to the hypothesis of a higher incidence of cryptic species in the tropics. They thus highlight the importance of large-scale studies and the necessity of testing new complementary conclusive tools to correctly quantify Neotropical biological diversity. Such research endeavours are certainly overwhelming, but they are essential for the understanding of the world's true diversity of life and the first step to assure its conservation.

Supporting Information

Figure S1.

Worker specimen of Pachycondyla verenae (PVE) Morph 1 from Petit Saut, French Guiana. Lateral view (A), Full-face view (B) and petiole, oblique lateral view (C).

https://doi.org/10.1371/journal.pone.0015363.s001

(TIF)

Figure S2.

Worker specimen of Pachycondyla apicalis (PAP) Morph 3 from Los Tuxlas, Mexico. Lateral view (A), Full-face view (B) and petiole, oblique lateral view (C).

https://doi.org/10.1371/journal.pone.0015363.s002

(TIF)

Figure S3.

Worker specimen of Pachycondyla apicalis (PAP) Morph 4 from Petit Saut, French Guiana. Lateral view (A), Full-face view (B) and petiole, oblique lateral view (C).

https://doi.org/10.1371/journal.pone.0015363.s003

(TIF)

Figure S4.

Worker specimen of Pachycondyla apicalis (PAP) Morph 5 from Petit Saut, French Guiana. Lateral view (A), Full-face view (B) and petiole, oblique lateral view (C).

https://doi.org/10.1371/journal.pone.0015363.s004

(TIF)

Figure S5.

Worker specimen of Pachycondyla apicalis (PAP) Morph 6 from Petit Saut, French Guiana. Lateral view (A), Full-face view (B) and petiole, oblique lateral view (C).

https://doi.org/10.1371/journal.pone.0015363.s005

(TIF)

Figure S6.

Scanning electron micrographs of the stridulatory file of five morphs from the Pachycondyla apicalis species complex. PVE Morph 1 (A), PAP Morph 3 (B), PAP Morph 4 (C), PAP Morph 5 (D) and PAP Morph 6 (E). PVE Morph 1 and PAP Morphs 4, 5, and 6 from Petit Saut, French Guiana and PAP Morph 3 from Los Tuxlas, Mexico. PVE: P. verenae, PAP: P. apicalis.

https://doi.org/10.1371/journal.pone.0015363.s006

(TIF)

Audio S1.

Series of chirps from a Pachycondyla verenae (PVE) Morph 1 worker.

https://doi.org/10.1371/journal.pone.0015363.s007

(WAV)

Audio S2.

Series of chirps from a Pachycondyla apicalis (PAP) Morph 3 worker.

https://doi.org/10.1371/journal.pone.0015363.s008

(WAV)

Audio S3.

Series of chirps from a Pachycondyla apicalis (PAP) Morph 4 worker.

https://doi.org/10.1371/journal.pone.0015363.s009

(WAV)

Audio S4.

Series of chirps from a Pachycondyla apicalis (PAP) Morph 5 worker.

https://doi.org/10.1371/journal.pone.0015363.s010

(WAV)

Audio S5.

Series of chirps from a Pachycondyla apicalis (PAP) Morph 6 worker.

https://doi.org/10.1371/journal.pone.0015363.s011

(WAV)

Acknowledgments

We gratefully acknowledge Marie-Claire Malherbe, Paul Devienne and Thierry Aubin for their great help with the SEM images, ant collection and acoustic analyses respectively. Thanks are also due to Nicolas Châline for his critical comments on early versions of the manuscript and to Philippe Cerdan (HYDRECO, French Guiana) for logistical support during field work. We are very grateful to the anonymous referee for revision comments and helpful suggestions.

Author Contributions

Conceived and designed the experiments: RSF CP DF FR. Performed the experiments: RSF CP FR. Analyzed the data: RSF CP JHCD DF FR. Wrote the paper: RSF CP JHCD DF FR.

References

  1. 1. Willig MR, Kaufman DM, Stevens RD (2003) Latitudinal gradients of biodiversity: Pattern, process, scale, and synthesis. Annu Rev Ecol Evol Syst 34: 273–309.
  2. 2. Bickford D, Lohman DJ, Sodhi NS, Ng PKL, Meier R, et al. (2007) Cryptic species as a window on diversity and conservation. Trends Ecol Evol 22: 148–155.
  3. 3. Schlick-Steiner BC, Steiner FM, Seifert B, Stauffer C, Christian E, et al. (2010) Integrative Taxonomy: A Multisource Approach to Exploring Biodiversity. Annu Rev Entomol 55: 421–438.
  4. 4. Balakrishnan R (2005) Species Concepts, Species Boundaries and Species Identification: A View from the Tropics. Syst Biol 54: 689–693.
  5. 5. Claridge MF, Dawah HA, Wilson MR (1997) Species in insect herbivores and parasitoids - sibling species, host races and biotypes. In: Claridge MF, Dawah HA, Wilson MR, editors. Species: The units of biodiversity. Chapman & Hall. pp. 247–272.
  6. 6. Henry CS (1994) Singing and cryptic speciation in insects. Trends Ecol Evol 9: 388–392.
  7. 7. Broza M, Blondheim Y, Nevo E (1998) New species of mole crickets of the Gryllotalpa gryllotalpa group (Orthoptera: Gryllotalpidae) from Israel, based on morphology, song recordings, chromosomes and cuticular hydrocarbons, with comments on the distribution of the group in Europe and the Mediterranean region. Syst Entomol 23: 125–135.
  8. 8. Cade WH, Otte D (2000) Gryllus texensis n. sp.: A widely studied field cricket (Orthoptera; Gryllidae) from the southern United States. Trans Am Entomol Soc 126: 117–123.
  9. 9. Walker TJ, Forrest TG, Spooner JD (2003) The rotundifolia complex of the genus Amblycorypha (Orthoptera: Tettigoniidae): Songs reveal new species. Ann Entomol Soc Am 96: 433–447.
  10. 10. Ritchie MG, Gleason JM (1995) Rapid evolution of courtship song pattern in Drosophila willistoni sibling species. J Evol Biol 8: 463–479.
  11. 11. Noor MAF, Aquadro CF (1998) Courtship songs of Drosophila pseudoobscura and D. persimilis: analysis of variation. Anim Behav 56: 115–125.
  12. 12. Sueur J (2006) Insect species and their songs. In: Drosopoulos S, Claridge M, editors. Insect Sounds and Communication: Physiology, Behaviour, Ecology, and Evolution. Boca Raton, FL: Taylor & Francis. pp. 207–217.
  13. 13. Wilcox TP, Hugg L, Zeh JA, Zeh DW (1997) Mitochondrial DNA sequencing reveals extreme genetic differentiation in a cryptic species complex of neotropical pseudoscorpions. Mol Phylogenet Evol 7: 208–216.
  14. 14. Berkov A (2002) The impact of redefined species limits in Palame (Coleoptera: Cerambycidae: Lamiinae: Acanthocinini) on assessments of host, seasonal, and stratum specificity. Biol J Linn Soc 76: 195–209.
  15. 15. Hebert PDN, Penton EH, Burns JM, Janzen DH, Hallwachs W (2004) Ten species in one: DNA barcoding reveals cryptic species in the neotropical skipper butterfly Astraptes fulgerator. Proc Natl Acad Sci U S A 101: 14812–14817.
  16. 16. Smith MA, Rodriguez JJ, Whitfield JB, Deans AR, Janzen DH, et al. (2008) Extreme diversity of tropical parasitoid wasps exposed by iterative integration of natural history, DNA barcoding, morphology, and collections. Proc Natl Acad Sci USA 105: 12359–12364.
  17. 17. Landois H (1874) Thierstimmen. Freiburg i. Br. Die Sprache der Ameisen 129–134.
  18. 18. Lubbock J (1877) On some points in the anatomy of ants. Mon Microsc J. 18. : 121–142, pls. 189-192.
  19. 19. Janet C (1893) Note sur la production des sons chez les fourmis et sur les organes qui les produisent. Ann Soc Entomol France 62: 159–167.
  20. 20. Janet C (1894) Étude sur les fourmis. 6e note. Sur l'appareil de stridulation de Myrmica rubra L. Ann Soc Entomol France 63: 109–117.
  21. 21. Sharp D (1893) On stridulation in ants. Trans Entomol Soc London 1893: 199–213.
  22. 22. Markl H (1973) The evolution of stridulatory communication in ants. In:Proceedings IUSSI VIIth International Congress, London, 10-15 September, 1973. Southampton: University of Southampton. pp. 258–265.
  23. 23. Taylor RW (1978) Nothomyrmecia macrops: a living-fossil ant rediscovered. Science 201: 979–985.
  24. 24. Spangler HG (1967) Ant stridulations and their synchronization with abdominal movement. Science 155: 1687–1689.
  25. 25. Kermarrec A, Mauléon H, Antun AA (1976) La stridulation de Acromyrmex octospinosus Reich. (Formicidae, Attini): biométrie de l'appareil stridulateur et analyse du signal produit. Insect Soc 23: 29–47.
  26. 26. Hölldobler B, Wilson EO (1990) The ants. Cambridge, Mass: Harvard University Press. 732 p.
  27. 27. Roces F, Tautz J, Hölldobler B (1993) Stridulation in leaf-cutting ants: short-range recruitment through plant borne vibrations. Naturwissenschaften 80: 521–524.
  28. 28. Roces F, Tautz J (2001) “Ants are deaf” J Acoust Soc Amer. 109. : 3080–3082.
  29. 29. Hickling R, Brown RL (2000) Analysis of acoustic communication by ants. J Acoust Soc Amer 108: 1920–1929.
  30. 30. Pavan G, Priano M, De Carli P, Fanfani A, Giovannotti M (1997) Stridulatory organ and ultrasonic emission in certain species of ponerine ants (genus: Ectatomma and Pachycondyla, Hymenoptera, Formicidae). Bioacoustics 8: 209–221.
  31. 31. Barbero F, Bonelli S, Thomas JA, Balletto E, Schonrogge K (2009) Acoustical mimicry in a predatory social parasite of ants. J Exp Biol 212: 4084–4090.
  32. 32. Barbero F, Thomas JA, Bonelli S, Balletto E, Schönrogge K (2009) Queen ants make distinctive sounds that are mimicked by a butterfly social parasite. Science 323: 782–785.
  33. 33. Hölldobler B, Stanton RC, Markl H (1978) Recruitment and food-retrieving behavior in Novomessor (Formicidae, Hymenoptera). I. Chemical signals. Behav Ecol Sociobiol 4: 163–181.
  34. 34. Baroni Urbani C, Buser MW, Schilliger E (1988) Substrate vibration during recruitment in ant social organization. Insect Soc 35: 241–250.
  35. 35. Stuart RJ, Bell PD (1980) Stridulation by workers of the ant, Leptothorax muscorum (Nylander) (Hymenoptera, Formicidae). Psyche 87: 199–210.
  36. 36. Maschwitz U, Schönegge P (1983) Forage communication, nest moving recruitment, and prey specialization in the oriental ponerine Leptogenys chinensis. Oecologia 57: 175–182.
  37. 37. Markl H (1965) Stridulation in leaf-cutting ants. Science 149: 1392–1393.
  38. 38. Grasso DA, Priano M, Pavan G, Mori A, Le Moli F (2000) Stridulation in four species of Messor ants (Hymenoptera, Formicidae). Ital J Zool 67: 281–283.
  39. 39. Markl H, Hölldobler B, Hölldobler T (1977) Mating behavior and sound production in harvester ants (Pogonomyrmex, Formicidae). Insect Soc 24: 191–212.
  40. 40. Mercier JL, Lenoir JC, Eberhardt A, Frohschammer S, Williams C, et al. (2007) Hammering, mauling, and kissing: stereotyped courtship behavior in Cardiocondyla ants. Insect Soc 54: 403–411.
  41. 41. DeVries PJ (1990) Enhancement of symbioses between butterfly caterpillars and ants by vibrational communication. Science 248: 1104–1106.
  42. 42. DeVries PJ (1991) Call production by myrmecophilous riodinid and lycaenid butterfly caterpillars (Lepidoptera): morphological, acoustical, functional, and evolutionary patterns. Am Mus Novit 3025: 1–23.
  43. 43. Travassos MA, Pierce NE (2000) Acoustics, context and function of vibrational signalling in a lycaenid butterfly-ant mutualism. Anim Behav 60: 13–26.
  44. 44. Morales MA, Barone JL, Henry CS (2008) Acoustic alarm signalling facilitates predator protection of treehoppers by mutualist ant bodyguards. Proc R Soc London B Biol Sci 275: 1935–1941.
  45. 45. Ferreira RS, Fresneau D (2009) Stridulation: Le chant méconnu des fourmis. In: Megherbi H, Fourment MC, editors. Langages Cahiers de l'Infantile. Paris: L'Harmattan. pp. 71–91.
  46. 46. Bolton B (1995) A new general catalogue of the ants of the world. Cambridge, Mass: Harvard University Press. 504 p.
  47. 47. Fernández F, Sendoya S (2004) List of Neotropical ants (Hymenoptera: Formicidae). Biota Colomb 5: 3–93.
  48. 48. Majer JD, Delabie JHC (1993) An evaluation of Brazilian cocoa farm ants as potential biological control agents. J Plant Prot Tropics 10: 43–49.
  49. 49. Majer JD, Delabie JHC (1999) Impact of tree isolation on arboreal and ground ant communities in cleared pasture in the Atlantic rain forest region of Bahia, Brazil. Insect Soc 46: 281–290.
  50. 50. Delabie JHC, Fisher BL, Majer J, Wright IW (2000) Sampling effort and choice of methods. In: Agosti D, Majer DJ, Alonso L, Schultz T, editors. Measuring and monitoring biodiversity: standard methods for ground-dwelling ants. Washington D.C.: Smithsonian Institution Press. pp. 145–154.
  51. 51. Mariano CSF, Pompolo SG, Delabie JHC (2000) Citogenética das espécies gêmeas e simpátricas Pachycondyla villosa e Pachycondyla sp ‘inversa’ (Ponerinae). Naturalia 24: 215–217.
  52. 52. Lucas C, Fresneau D, Kolmer K, Heinze J, Delabie JHC, et al. (2002) A multidisciplinary approach to discriminating different taxa in the species complex Pachycondyla villosa (Formicidae). Biol J Linn Soc 75: 249–259.
  53. 53. Mariano CSF, Pompolo SG, Borges DS, Delabie JHC (2006) Are the Neotropical ants Pachycondyla crenata (Roger, 1861) and Pachycondyla mesonotalis (Santschi, 1923) (Hymenoptera: Formicidae: Ponerinae) good species? A cytogenetic approach. Myrmecol Nachr 8: 277–280.
  54. 54. Wild AL (2005) Taxonomic revision of the Pachycondyla apicalis species complex (Hymenoptera: Formicidae). Zootaxa 834: 1–25.
  55. 55. Fresneau D (1994) Biologie et comportement social d'une fourmi ponerine neotropicale (Pachycondyla apicalis): Ph.D. Thesis. Paris XIII: Université Villetaneuse. 331 p.
  56. 56. Reiskind J (1977) Ant-mimicry in Panamanian clubionid and salticid spiders (Araneae: Clubioinidae, Salticidae). Biotropica 9: 1–8.
  57. 57. McIver JD, Stonedahl G (1993) Myrmecomorphy: morphological and behavioral mimicry of ants. Annu Rev Entomol 38: 351–379.
  58. 58. Cushing PE (1997) Myrmecomorphy and myrmecophily in spiders: A review. Fla Entomol 80: 165–193.
  59. 59. Delabie JHC, Mariano CSF, Mendes LF, Pompolo SG, Fresneau D (2008) Problemas apontados por estudos morfológicos, ecológicos e citogenéticos no gênero Pachycondyla na região Neotropical: o caso do complexo apicalis. In: Vilela EF, Santos IA, Schoereder JH, Neto JL, Serrão JE, Campos LAO, editors. Insetos Sociais: Da Biologia à Aplicação. Viçosa: Editora UVF. pp. 197–222.
  60. 60. Latreille PA (1802) Histoire naturelle des fourmis, et recueil de mémoires et d'observations sur les abeilles, les araignées, les faucheurs, et autres insectes. Paris: De l'impr. de Crapelet (chez T. Barrois). 445 p.
  61. 61. Emery C (1890) Voyage de M. E. Simon au Venezuela (Décembre 1887 - Avril 1888). (6): Formicides. Ann Soc Entomol France. pp. 55–76.
  62. 62. Forel A (1922) Glanures myrmécologiques en 1922. Rev Suisse Zool 30: 87–102.
  63. 63. Brown WL (1957) Ants from Laguna Ocotal (Hymenoptera: Formicidae). Bull Mus Comp Zool 116: 228–237.
  64. 64. Specht R (2008) Avisoft-SASLab version 4.40. Berlin: Avisoft Bioacoustics. Available from http://www.avisoft.com.
  65. 65. Mbu Nyamsi RG, Aubin T, Bremond JC (1994) On the extraction of some time dependent parameters of an acoustic signal by means of the analytic signal concept. Its application to animal sound study. Bioacoustics 5: 187–203.
  66. 66. Simon C, Frati F, Beckenbach A, Crespi B, Liu H, et al. (1994) Evolution, weighting, and phylogenetic utility of mitochondrial gene-sequences and a compilation of conserved polymerase chain-reaction primers. Ann Entomol Soc Am 87: 651–701.
  67. 67. Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG (1997) The Clustal X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res 25: 4876–4882.
  68. 68. Kimura M (1980) A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. J Mol Evol 16: 111–120.
  69. 69. Tamura K, Dudley J, Nei M, Kumar S (2007) MEGA4: Molecular evolutionary genetics analysis (MEGA) software version 4.0. Mol Biol Evol 24: 1596–1599.
  70. 70. Tajima F (1993) Simple methods for testing molecular clock hypothesis. Genetics 135: 599–607.
  71. 71. Spjotvol E, Stoline MR (1973) Extension of t-method of multiple comparison to include cases with unequal sample sizes. J Am Stat Assoc 68: 975–978.
  72. 72. Box GEP, Cox DR (1964) An analysis of transformations. J R Stat Soc Ser B-Stat Methodol 26: 211–252.
  73. 73. StatSoft Bremond (2007) STATISTICA (logiciel d'analyses de données), version 8.0. www.statsoft.fr.
  74. 74. Crozier RH, Crozier YC (1993) The mitochondrial genome of the honeybee Apis mellifera - complete sequence and genome organization. Genetics 133: 97–117.
  75. 75. Grasso DA, Mori A, Le Moli F, Giovannotti M, Fanfani A (1998) The stridulatory organ of four Messor ant species (Hymenoptera, Formicidae). Ital J Zool 65: 167–174.
  76. 76. Markl H (1968) Die Verständigung durch Stridulationssignale bei Blattschneiderameisen, II: Erzeugung und Eigenschaften der Signale. Z vgl Physiol 60: 103–150.
  77. 77. Hölldobler B, Wilson EO (2009) The Superorganism: the beauty, elegance and strangeness of insect societies. New York: W.W. Norton and Company. 576 p.
  78. 78. Song H, Buhay JE, Whiting MF, Crandall KA (2008) Many species in one: DNA barcoding overestimates the number of species when nuclear mitochondrial pseudogenes are coamplified. Proc Natl Acad Sci USA 105: 13486–13491.
  79. 79. Saez AG, Lozano E (2005) Body doubles. Nature 433: 111–111.
  80. 80. Gavrilets S (2003) Perspective: Models of speciation: What have we learned in 40 years? Evolution 57: 2197–2215.
  81. 81. Molbo D, Machado CA, Sevenster JG, Keller L, Herre EA (2003) Cryptic species of fig-pollinating wasps: Implications for the evolution of the fig-wasp mutualism, sex allocation, and precision of adaptation. Proc Natl Acad Sci USA 100: 5867–5872.