PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of plosonePLoS OneView this ArticleSubmit to PLoSGet E-mail AlertsContact UsPublic Library of Science (PLoS)
 
PLoS One. 2010; 5(12): e14405.
Published online 2010 December 21. doi:  10.1371/journal.pone.0014405
PMCID: PMC3014355

Cuticular Compounds Bring New Insight in the Post-Glacial Recolonization of a Pyrenean Area: Deutonura deficiens Deharveng, 1979 Complex, a Case Study

Simon Joly, Editor

Abstract

Background

In most Arthropod groups, the study of systematics and evolution rely mostly on neutral characters, in this context cuticular compounds, as non-neutral characters, represent an underexplored but potentially informative type of characters at the infraspecific level as they have been routinely proven to be involved in sexual attraction.

Methods and Findings

The collembolan species complex Deutonura deficiens was chosen as a model in order to test the utility of these characters for delineating four infraspecific entities of this group. Specimens were collected for three subspecies (D. d. deficiens, D. d. meridionalis, D. d. sylvatica) and two morphotypes (D. d. sylvatica morphoype A and B) of the complex; an additional species D. monticola was added. Cuticular compounds were extracted and separated by gas chromatography for each individual. Our results demonstrate that cuticular compounds succeeded in separating the different elements of this complex. Those data allowed also the reconstruction of the phylogenetic relationships among them.

Conclusions

The discriminating power of cuticular compounds is directly related to their involvement in sexual attraction and mate recognition. These findings allowed a discussion on the potential involvement of intrinsic and paleoclimatic factors in the origin and the diversification of this complex in the Pyrenean zone. This character type brings the first advance from pattern to process concerning the origin of this species complex.

Introduction

Species and subspecies are taxonomic levels generally considered as objective descriptors of biological reality, a mere fact at the origin of the numerous different studies proposed to date. Whereas long diverged lineages are generally easy to separate with various sets of characters, recently diverged lineages are usually much more difficult to address. In those cases, classical morphological approaches sometimes fail to provide consistent discrimination or hypotheses about evolutionary history and relationships, and benefit from the inclusion of other markers with more appropriate and objectively measurable rates of divergence [1]. Besides the most commonly used neutral genetic markers such as ISSR [2], microsatellites [3] nuclear and mitochondrial sequences [4.5], allozymes [6] and karyotypes [7], [8], alternative sets of characters, such as cuticular compounds, have been proposed and tested with promising results [9], [10]. These compounds are lipids present on the most external cuticular layer of all terrestrial arthropods: the epicuticular layer. This layer is composed of various categories of lipids. In most cases, the majority of these compounds are hydrocarbons but other types of compounds can be present such as fatty acids, ketones, alcohols and esters [11], [12]. As the benefit of using the maximum number of compounds for population level studies has already been emphasized in the literature [10], in the present study we consider the total composition of the cuticular profile regardless of families to which the chemical compounds belong. In insects, the different types of cuticular compounds are functionally involved, particularly in sexual communication [13], [14]. In Collembola, intraspecific interactions such as reproductive attraction and recognition signalling [15], [16], aggregation [17], [18], [19], [20], [21] and alarm [22], [23], [24] rely heavily on olfactory interactions through chemical compounds.

As they also play a major role in sexual recognition [25], [26] and thus in premating isolation mechanisms [27], cuticular compounds, as a type of character, deserve a special attention in the study of groups of subspecific entities. This is especially relevant for Collembola since cuticular compounds were proven to play a major role in sexual recognition [16] suggesting that they could be potential key drivers of reproductive isolation, making them good candidates for infraspecific taxon discrimination in this group.

A previous population study of the wingless grasshopper Chorthippus parallelus from the Pyrenees revealed significant differences between subspecies cuticular profiles and suggested their major implication in specific sexual recognition [28]. A quantitative sexual dimorphism in cuticular blends was found in this grasshopper species [10], [29] and subsequently those compounds were proven to be involved in reproductive isolation mechanisms [30]. A similar sexual dimorphism was also demonstrated in Collembola [31].

In the present study, we tested the discriminating power of such characters at the infraspecific level for Collembola using the Deutonura deficiens complex as a model. This thoroughly studied group displays a very complex pattern of closely related subspecies and morphotypes living in parapatry within a relatively limited geographic range. Narrow overlapping areas and hybrid zones were documented in the Pyrenees region [32], [33]. This complex includes three subspecies: D. d. deficiens, D. d. meridionalis and D. d. sylvatica. Each one splits into several morphotypes that are geographically well delimited. These infraspecific taxonomic categories defined by Deharveng [32] are based on their distribution and morphological similarity within populations (cephalic tuberculisation among subspecies and abdominal setae and body pigmentation for morphotypes within subspecies).

The extensive knowledge of the repartition areas of these different infraspecific entities allowed an overall assessment of the role that ecological factors played in producing the current biogeographical pattern and the position of the contact zones [32]. None of the ecological factors examined (vegetation, altitude, precipitation, temperature or hydrographic network) were found to be significantly related to the species boundaries within the mosaic pattern [33]. Concerning the origin of the pattern, these results suggest the involvement of inherent characteristics of the organisms and/or an historical inheritance.

In this context we address the utility of cuticular compounds as a set of characters for (1) separating infraspecific lineages in Collembola and (2) establishing phylogenetic relationships between them. We also discuss the origin of the parapatric pattern in the context of European glacial history.

Materials and Methods

Samples

Specimens were collected in the field, either by direct capture on substrate or by extraction from litter and decayed wood with Berlese funnels. They were freeze killed and stored at -28°C. In order to avoid the potential bias of variation of the cuticular profile with age [34], only adult individuals were used. In this study, 21 individuals from the D. d. deficiens subspecies, 13 from D. d. meridionalis and 21 from D. d. sylvatica (with 12 individuals from morphotype C and 9 from morphotype A) were analysed (Table 1). In order to infer phylogenetic relationships among these infraspecific entities, we used an external species, Deutonura monticola (17 individuals), to root the evolutionary hypothesis. Identifications to subspecies and morphotypes were carried out through morphological examination of specimens after microscopic preparation of all individuals.

Table 1
Size and locations of sampled populations.

Gas Chromatography

Cuticular lipids were extracted by soaking each individual for 5 minutes in 20 µl of hexane. Separation of compounds in the extracts was performed on a Hewlett Packard Chromatograph model 5890 equipped with a Flame Ionization Detector (320°C) and an automatic sampler (HP 6890). For each individual, 5 µl of extract was injected in an apolar capillary column (HT-5) of 25m long. Helium was used as the carrier gas (pressure 16 psi - injector temperature 280°C). The temperature program started from 150°C and reached 320°C with a rate of 5°C/min. Then a 20 minutes 320°C isotherm was used to clean up the column between two samples. Each compound was characterized with its Kováts index (KI) [35] calculated using the coelution of samples with a standard solution of linear alkanes (C16 to C40) (Table S1). This homology hypothesis relying on KI brings limitations because compounds were not structurally identified. Some compounds even if structurally different can have the same KI. In some cases, this could lead to an underestimation of the variability of the profiles.

Statistical Analysis

The data were recovered with the Millennium 2.15 software (Waters). For each individual, the relative surfaces of each peak in the chromatogram matched the concentration of the corresponding cuticular compound in the blend. Those relative surfaces were calculated and the data were exported to Systat 10.0 for statistical treatment. A stepwise forward discriminant analysis was performed to explore infraspecific variation, and the “F-to-remove” value was retained for each discriminant peak. This value is a good indicator of peak importance for distinction between groups. But a discriminant analysis alone can eliminate highly correlated values and thus give a partial view of the data set separating different classes. For this reason we also performed Pearson Correlation Matrix with a Bonferroni correction (P<0.005) to recover the correlated peaks (Table S1 - Supplementary material). A jackknife resampling method (default for discriminant analysis in the SYSTAT 10.0 software) was used to evaluate the discriminatory power of the cuticular compounds.

Phylogenetic Reconstruction

Distance relations were established between the different infraspecific entities through a Neighbor Joining calculation with MEGA 4 [36]. The Mahalanobis distances matrix obtained from the previous discriminant analysis (i.e. distances between centroïds of the different infra-specific entities) was used for this purpose (Table 2). The external species D. monticola was used to root the phylogenetic hypothesis.

Table 2
Mahalanobis distances matrix obtained via the discriminant analysis.

Results

A total of 128 distinct compounds were found for D. deficiens subspecies, morphotypes, and the external species Deutonura monticola. After the statistical analysis, we obtained a significant discrimination (P<0.00001, F = 27.6126) between the four infra-specific taxa and D. monticola (Fig. 1). The Jacknife classification value was 100% for all taxa but D. d. sylvatica morphotype C with one individual misplaced (Table 3). Thirty-five cuticular compounds were designated as discriminant. Eighty-two other compounds were recovered using the correlation matrix. Eventually, 117 products allowed separation of the groups (Table S1 - Supplementary material). All variations were quantitative. Most compounds are implicated in the discrimination, indicating profile proximity between the different morphotypes and subspecies with numerous quantitative variations at a fine scale.

Figure 1
Graphic representation of the discriminant analysis following the two first axis accounting for 92.6% of the variability.
Table 3
Jackniffe assignation values of the discriminant analysis for each infraspecific entity.

Phylogenetic relationships were established and rooted with Deutonura monticola (Fig. 2). The subspecies D. d. sylvatica appeared as the first derived form of the lineage with its two morphotypes aggregated in basal position. Then the two other subspecies, D. d. deficiens and D. d. meridionalis, branched together suggesting their later individualisation.

Figure 2
Phylogenetic relationships between subspecies and morphotypes of Deutonura deficiens build from the Mahalanobis distances rooted with D. monticola.

Discussion

During the last twenty years, cuticular compounds have been investigated in chemotaxonomy for species specificity [37], [38], [39], [40], [41]. But the systematic use of those characters for taxonomy was challenged by some findings in termites: if cuticular compounds were proven useful in this group for species delimitation [40], [42], a strong variation of the profiles was found among several colonies of the same species [43], [44], [45]. But this is the peculiar case of social insect in which conspecific colonies are competing with each other and thus have a strong pressure for chemical recognition cues among colonies [13], [46], [47] even if a ‘dear enemy’ effect was demonstrated in termites [48]. It was suggested that a part of the cuticular compounds can be incorporated from some of the substances found in the alimentation as a possible explanation for the high variability of profiles among conspecific colonies of termites [48], such a phenomenon was already proven to occur in ants [49].

By contrast there is not such a competitive pressure among individuals in solitary animals like collembolans and thus a selection pressure toward a strong individualization of their cuticular profiles cannot be expected. Furthermore the repartition areas of the infraspecific entities of the species complex studied in this paper are very restricted and so the variations of feeding possibilities are expected to be quite narrow on such a range where the types of habitats are similar. The low variability of the collembolans cuticular profile within and among different populations of the species was documented in this respect [50] which allows to push the test of the taxonomic resolution further here.

In this study, cuticular compounds successfully discriminated all of the infraspecific taxa assigned by morphology within the Deutonura deficiens complex. The congruence between the two types of characters supports the biological reality of those entities. In addition, they brought further understanding of the phylogenetic relationships among these taxa. Those data confirmed that in Collembola, this character type is phylogenetically informative at low taxonomical levels [50].

The phylogenetic relationships established from cuticular profiles support two hypotheses: (a) an early individualization of D. d. sylvatica morphotypes A and C, which is consistent with the current sympatry of some of their populations exhibiting a high density of both parental phenotypes and an absence of hybrids; such a pattern and their basal position in the chemical phylogeny support the hypothesis of an achieved speciation phase with pre-reproductive isolation completed as hypothesized by Deharveng [32], (b) a more recent appearance of D. d. deficiens and D. d. meridionalis; the strict parapatry and the occurrence of morphologically intermediate specimens (assumed to be hybrids) among mixed populations in narrow contact zones between those two subspecies provides further support for their recent individualization.

As these chemicals play a crucial role in sexual attraction and mate recognition in Collembola [16], the discrepancies found among the cuticular profiles of all the infraspecific entities of the D. deficiens complex strongly suggest a status of an ongoing process of reproductive isolation involving those compounds. Clear divergence of cuticular profiles between sexes was demonstrated for other Collembola species [31]. Such sexual dimorphism exists also in D. deficiens, and all but one of the cuticular compounds involved in the sexual variation were implicated in the subspecies discrimination (Table S2 - Supplementary material). This indicates that sexual selection is mostly responsible for the variation of these traits, implying a major role of the cuticular compounds in sexual and specific recognition, which fits with the last findings in this field [16]. Then the possibility of a reinforcement phenomenon weighting on these characters as described by Higgie and Blows [51] could be considered.

The elements brought by the analysis of cuticular profiles suggest that several events of allopatry could have led to different levels of reproductive isolation responsible for the appearance of the distinct intraspecific entities and the origin of the current parapatric pattern of the D. deficiens complex. Pleistocene glacial events could be at the origin of the occurrence of allopatric phases for several populations of this species. A comparable mosaic distribution of the subspecies of another Collembola species of the same subfamily (Monobella grassei) was uncovered in the same region [52]. One of the hypothetical causes examined for such a pattern was the major historical event of the Pleistocene glaciations. These key paleoclimatic episodes are mainly accountable for the current repartition of animal and plant species populations: a correlation was established between variation in vegetal cover and climatic osrelationships were established cillations in the last 135 000 years [53], and Hewitt [54] proposed three European models for post-glacial recolonization in the Pleistocene based on plant and animal cases. The knowledge of the geographic distribution of the taxa studied here allows a comparison with Hewitt's recolonization models [54] in order to discuss the origin of the actual pattern (described in [33]) in a historical perspective. One of the postglacial recolonization models described by Hewitt [54] matches the existing distribution of the D. deficiens subspecies and the location of the contact zones (the recent discovery of a D. d. sylvatica population in the Italian Peninsula strengthens this support). It implies a faster advance of the populations from eastern refuges (Balkan, southern Carpathians and Caucasus) by contrast with the populations from the Iberic and Italian Peninsula refuges. As a result, the northern zones of these southern refuges were colonized by eastern populations blocking the advance of the southern ones and leading to the occurrence of contact zones. The current repartition and the contact zones position observed in the D. deficiens complex is congruent with the scenario of post-glacial recolonization exhibited by several taxa [54]. This brings support to the involvement of the last Pleistocene glacial events in the establishment of this distribution pattern.

Though we acknowledge some limitations in the understanding of some aspects of the complex wider geographic repartition and history: the possible occurrence of populations belonging to D. deficiens complex in eastern refuges, or the role played by Atlantic and Mediterranean microhabitats near actual contact zones which could have interfered with this model at fine scales as described in the ‘refugia within refugia’ concept developed for the Balkan Peninsula [55].

The successful use of a type of character directly involved in the mate choice on genetically distinct but closely related entities in a parapatric context, underlies their predominance in the origin and the maintenance of this parapatric pattern. The complexity of the whole pattern calls for a comparison of other markes types and a more exhaustive sampling through existing and potential repartition areas of the complex. The elucidation of the origin and the maintenance of such a pattern is crucial to understand what has driven the repartition and diversity of understudied soil microarthropods in this critical European zone.

Supporting Information

Table S1

Discriminant analysis results for the 4 infraspecific entities and the external species D. monticola.

(0.24 MB DOC)

Table S2

Discriminant analysis results between sexes in D. d. meridionalis (P<0.00001, F = 193,3723; Jackknife 100%).

(0.03 MB DOC)

Acknowledgments

The authors want to thank here André Fabres for his advice on statistical analyses, Sarah Adamowicz and Rodolphe Rougerie for their useful comments on the manuscript.

Footnotes

Competing Interests: The authors have declared that no competing interests exist.

Funding: This research was funded by the Laboratoire d'Ecologie Terrestre (Université Paul Sabatier - CNRS). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

References

1. Audzijonyte A, Ovcarenko I, Bastrop R, Vainola R. Two cryptic species of the Hediste diversicolor group (Polychaeta, Nereididae) in the Baltic Sea, with mitochondrial signatures of different population histories. Mar Biol. 2008;155:599–612.
2. Gonzalez J, Wink M, Garcia-Del-Rey E, Castro GD. Evidence from DNA nucleotide sequences and ISSR profiles indicates paraphyly in subspecies of the Southern Grey Shrike (Lanius meridionalis). J Ornithol. 2008;149:495–506.
3. Gava A, de Freitas TRO. Microsatellite analysis of a hybrid zone between chromosomally divergent populations of Ctenomys minutus from southern brazil (Rodentia: Ctenomyidae). J Mammal. 2004;85:1201–1206.
4. Hwang DS, Ki JS, Jeong DH, Kim BH, Lee BK, et al. A comprehensive analysis of three Asiatic black bear mitochondrial genomes (subspecies ussuricus, formosanus and mupinensis), with emphasis on the complete mtDNA sequence of Ursus thibetanus ussuricus (Ursidae). Mitochondr. DNA. 2008;19:418–429. [PubMed]
5. Godinho R, Crespo EG, Ferrand N. The limits of mtDNA phylogeography: complex patterns of population history in a highly structured Iberian lizard are only revealed by the use of nuclear markers. Mol Ecol. 2008;17:4670–4683. [PubMed]
6. Driscoll DA, Roberts JD. A hybrid zone defined by allozymes and ventral colour in Geocrinia rosea (Anura : Myobatrachidae). Aust J Zool. 2007;55:371–376.
7. de Freitas TRO. Cytogenetics status of four Ctenomys species in the south of Brazil. Genetica. 2006;126:227–235. [PubMed]
8. Volobouev V, Auffray JC, Debat V, Denys C, Gautun JC, et al. Species delimitation in the Acomys cahirinus dimidiatus complex (Rodentia, Muridae) inferred from chromosomal and morphological analyses. Biol J Linnean Soc. 2007;91:203–214.
9. Tregenza T, Pritchard VL, Butlin RK. Patterns of trait divergence between populations of the meadow grasshopper, Chorthippus parallelus. Evolution. 2000c;54:574–585. [PubMed]
10. Buckley SH, Tregenza T, Butlin RZ. Transitions in cuticular composition across a hybrid zone: historical accident or environmental adaptation? Biol J Linnean Soc. 2003;78:193–201.
11. Lockey KH. Lipids of the insect cuticle: origin, composition and function. Comp Biochem Physiol. 1988;89B:595–645.
12. Nelson DR, Charlet LD. Cuticular hydrocarbons of the sunflower beetle, Zygogramma exclamationis. Comp Biochem Physiol B. 2003;135:273–284. [PubMed]
13. Singer TL. Role of hydrocarbons in the recognition systems of insects. American Zoologist. 1998;38:394–405.
14. Nojima S, Shimomura K, Honda H, Yamamoto I, Ohsawa K. Contact sex pheromone components of the cowpea weevil, Callosobruchus maculatus. J Chem Ecol. 2007;33:923–933. [PubMed]
15. Waldorf ES. Sex pheromone in the Springtail, Sinella curviseta. Environ Entomol. 1974;3:916–918.
16. Porco D, Deharveng L, Skarżyńskic D. Sex pheromone in Tetrodontophora bielanensis (Waga, 1842) (Collembola: Onychiuridae): Indirect reproduction mediated by cuticular compounds. Pedobiologia. 2009;53:59–63.
17. Mertens J, Bourgoignie R. Aggregation Pheromone in Hypogastrura viatica (Collembola). Behav Ecol Sociobiol. 1977;2:41–48.
18. Joosse ENG, Koelman TACM. Evidence for the presence of aggregation pheromones in Onychiurus armatus Collembola a pest insect in sugar beet. Entomol Exp Appl. 1979;26:197–201.
19. Verhoef HA, Nagelkerke CJ, Joose ENG. Aggregation pheromones in Collembola. J Insect Physiol. 1977;23:1009–1013.
20. Manica A, McMeechan FK, Foster WA. An aggregation pheromone in the intertidal collembolan Anurida maritima. Entomol Exp Appl. 2001;99:393–395.
21. Benoit JB, Elnitsky MA, Schulte GG, Lee RE, Denlinger DL. Antarctic Collembolans use chemical signals to promote aggregation and egg laying. J Insect Behav. 2009;22:121–133.
22. Purrington FF, Kendall PA, Bater JE, Stinner BR. Alarm pheromone in a gregarious Poduromorph collembolan (Collembola, Hypogastruridae). Gt Lakes Entomol. 1991;24:75–78.
23. Messer C, Dettner K, Schulz S, Francke W. Phenolic compounds in Neanura muscorum (Collembola, Neanuridae) and the role of 1,3-dimethoxybenzene as alarming substance. Pedobiologia. 1999;43:174–182.
24. Porco D, Deharveng L. 1,3-Dimethoxybenzene, a chemotaxonomic marker for the Neanurinae subfamily (Collembola). Biochem Syst Ecol. 2007;35:160–161.
25. Sugeno W, Hori M, Matsuda K. Identification of the contact sex pheromone of Gastrophysa atrocyanea (Coleoptera: Chrysomelidae). Appl Entomol Zool. 2006;41:269–276.
26. Peterson MA, Dobler S, Larson EL, Juárez D, Schlarbaum T, et al. Profiles of cuticular hydrocarbons mediate male mate choice and sexual isolation between hybridising Chrysochus (Coleoptera: Chrysomelidae). Chemoecology. 2007;17:87–96.
27. Smadja C, Butlin RK. On the scent of speciation: the chemosensory system and its role in premating isolation. Heredity. 2009;102:77–97. [PubMed]
28. Neems RM, Butlin RK. Divergence in cuticular hydrocarbons between parapatric subspecies of the meadow grasshopper, Chorthippus parallelus (Orthoptera, Acrididae). Biol J Linnean Soc. 1995;54:139–149.
29. Tregenza T, Buckley SH, Pritchard VL, Butlin RK. Inter and intrapopulation effects of sex and age on epicuticular composition of meadow grasshopper, Chorthippus parallelus. J Chem Ecol. 2000a;26:257–278.
30. Tregenza T, Pritchard VL, Butlin RK. The origins of the premating reproductive isolation: testing hypotheses in the grasshopper Chorthippus parallelus. Evolution. 2000b;54:1–13. [PubMed]
31. Porco D, Deharveng L, Gers C. Sexual discrimination with cuticular lipids in Schoettella ununguiculata (Tullberg, 1869) (Collembola: Hypogastruridae). Pedobiologia. 2004;48:581–583.
32. Deharveng L. Polymorphisme et polytypisme morphologiques chez quelques Neanurinae européens. Rev Ecol Biol Sol. 1984;21:533–562.
33. Rougerie R, Deharveng L, Porco D, Bedos A. The distribution of Deutonura deficiens morphotypes (Collembola: Neanurinae): geographical barriers, ecological gradients and historical factors. Pedobiologia. 2004;50:83–93.
34. Mpuru S, Blomquist G, Schal C, Roux M, Kuenzli M, et al. Effect of age and sex on the production of internal and external hydrocarbons and pheromones in the housefly, Musca domestica. Insect Biochem Mol Biol. 2001;31:139–155. [PubMed]
35. Kováts E. Gas chromatographic characterization of organic substances in the retention index system. Adv Chromatogr. 1965;1:229–247.
36. Tamura K, Dudley J, Nei M, Kumar S. MEGA4: Molecular evolutionary genetics analysis (MEGA) software version 4.0. Mol Biol Evol. 2007;24:1596–1599. [PubMed]
37. Carlson DA, Reinert JF, Bernier UR, Sutton BD, Seawright JA. Analysis of the cuticular hydrocarbons among species of the Anopheles quadrimaculatus complex (Diptera : Culicidae). J Am Mosq Control Assoc. 1997;13(Suppl):103–111. [PubMed]
38. Page M, Nelson LJ, Blomquist GJ, Seybold SJ. Cuticular hydrocarbons as chemotaxonomic characters of pine engraver beetles (Ips ssp.) in the grandicollis subgeneric group. J Chem Ecol. 1997;23:1053–1099.
39. Chapman RF, Espelie KE, Sword GA. Use of cuticular lipids in grasshopper taxonomy: a study of variation in Shcistocerca shoshone (Thomas). Biochem Sys Ecol. 1995;23:383–398.
40. Haverty MI, Woodrow RJ, Nelson LJ, Grace JK. Identification of termite species by the hydrocarbons in their feces. J Chem Ecol. 2005;31:2119–2151. [PubMed]
41. Baracchi D, Dapporto L, Teseo S, Hashim R, Turillazzi S. Medium molecular weight polar substances of the cuticle as tools in the study of the taxonomy, systematics and chemical ecology of tropical hover wasps (Hymenoptera: Stenogastrinae). J Zool Sys Evol Res. 2010;48:109–114.
42. Austin JW, Bagneres AG, Szalanski AL, Scheffrahn RH, Heintschel BP, et al. Reticulitermes malletei (Isoptera: Rhinotermitidae): a valid nearctic subterranean termite from eastern North America. Zootaxa. 2007;1554:1–26.
43. Bagine RKN, Brandl R, Kaib M. Species delimitation in Macrotermes (Isoptera, Macrotermitidae) - evidence from epicuticular hydrocarbons, morphology, and ecology. Ann Entomol Soc Am. 1994;87:498–506.
44. Haverty MI, Nelson LJ, Forschler BT. New cuticular hydrocarbon phenotypes of Reticulitermes (Isoptera : Rhinotermitidae) from the United States. Sociobiology. 1999;34:1–21.
45. Marten A, Kaib M, Brandl R. Cuticular hydrocarbon phenotypes do not indicate cryptic species in fungus-growing termites (Isoptera: Macrotermitinae). J Chem Ecol. 2009;35:572–579. [PubMed]
46. Clement JL, Bagneres AG. Nestmate recognition in termites. Pheromone communication in social insects: ants, wasps, bees and termites. In: Vander Meer RK, Breed MD, Espelie KE, Winston ML, editors. Westview Press; 1998. pp. 126–155.
47. Lahav S, Soroker V, Hefetz A, Vander Meer RK. Direct behavioral evidence for hydrocarbons as ant recognition discriminators. Naturwissenschaften. 1999;86:246–249.
48. Kaib M, Franke S, Francke W, Brandl R. Cuticular hydrocarbons in a termite: phenotypes and a neighbour-stranger effect. Physiol Entomol. 2002;27:189–198.
49. Buczkowski G, Silverman J. Geographical variation in Argentine ant aggression behaviour mediated by environmentally derived nestmate recognition cues. Anim Behav. 2006;71:327–335.
50. Porco D, Deharveng L. Phylogeny of Collembola based on cuticular compounds: inherent usefulness and limitation of a character type. Naturwissenschaften. 2009;96:943–954. [PubMed]
51. Higgie M, Blows MW. The evolution of reproductive character displacement conflicts with how sexual selection operates within a species. Evolution. 2008;62:1192–1203. [PubMed]
52. Deharveng L, Bedos A, Gisclard C. Environmental factors, microgeographic patterns of endemism and hybrid zones in Monobella grassei (Insecta : Collembola : Neanuridae). Biol J Linnean Soc. 1998;64:527–554.
53. Guiot J. Back to the last interglacial. Nature. 1997;388:25–27.
54. Hewitt GM. Post-glacial re-colonization of European biota. Biol J Linnean Soc. 1999;68:87–112.
55. Previsic A, Walton C, Kucinic M, Mitrikeski PT, Kerovec M. Pleistocene divergence of Dinaric Drusus endemics (Trichoptera, Limnephilidae) in multiple microrefugia within the Balkan Peninsula. Mol Ecol. 2009;18:634–647. [PubMed]

Articles from PLoS ONE are provided here courtesy of Public Library of Science