Molecular species delimitation of a symbiotic fig-pollinating wasp species complex reveals extreme deviation from reciprocal partner specificity
© Darwell et al.; licensee BioMed Central Ltd. 2014
Received: 22 May 2014
Accepted: 13 August 2014
Published: 18 September 2014
Symbiotic relationships have contributed to major evolutionary innovations, the maintenance of fundamental ecosystem functions, and the generation and maintenance of biodiversity. However, the exact nature of host/symbiont associations, which has important consequences for their dynamics, is often poorly known due to limited understanding of symbiont taxonomy and species diversity. Among classical symbioses, figs and their pollinating wasps constitute a highly diverse keystone resource in tropical forest and savannah environments. Historically, they were considered to exemplify extreme reciprocal partner specificity (one-to-one host-symbiont species relationships), but recent work has revealed several more complex cases. However, there is a striking lack of studies with the specific aims of assessing symbiont diversity and how this varies across the geographic range of the host.
Here, we use molecular methods to investigate cryptic diversity in the pollinating wasps of a widespread Australian fig species. Standard barcoding genes and methods were not conclusive, but incorporation of phylogenetic analyses and a recently developed nuclear barcoding gene (ITS2), gave strong support for five pollinator species. Each pollinator species was most common in a different geographic region, emphasising the importance of wide geographic sampling to uncover diversity, and the scope for divergence in coevolutionary trajectories across the host plant range. In addition, most regions had multiple coexisting pollinators, raising the question of how they coexist in apparently similar or identical resource niches.
Our study offers a striking example of extreme deviation from reciprocal partner specificity over the full geographical range of a fig-wasp system. It also suggests that superficially identical species may be able to co-exist in a mutualistic setting albeit at different frequencies in relation to their fig host¿s range. We show that comprehensive sampling and molecular taxonomic techniques may be required to uncover the true structure of cryptic biodiversity underpinning intimate ecological interactions.
KeywordsSpecies delimitation Agaonidae Ficus Pollinator Wasp
Symbiosis between disparate organisms has been responsible for key innovations in evolutionary history, such as the origin of plants . It is also pivotal in several vital ecosystem functions, such as nitrogen fixation  and pollination e.g. , and has been implicated as a major driver in the generation of biodiversity ,. The intimate nature and high partner specificity of many symbioses establishes conditions for strong interactions between the species involved, in terms of both population dynamics and coevolution ,. A basic but important question is “how many species are involved in a host/symbiont interaction?”, because species interaction strength, a major determinant of ecological and evolutionary dynamics, will differ between a system with one host and one symbiont species and one with multiple host or symbiont species . In addition, when multiple symbionts compete to utilise host resources, this may select for more selfish or virulent behaviour.
The obligate mutualism between fig trees and their pollinating wasps is a classic and much-studied example of symbiosis -. Fig trees are important components of rainforest and savannah ecosystems and provide food for many vertebrate and invertebrate animals . Figs can only be pollinated by host-specific wasp species, which, in turn, can only reproduce by laying their eggs in fig flowers, upon which their developing larvae feed. This association was long considered to be a textbook example of extreme reciprocal partner specificity, with each fig species having a unique pollinator species . However, this picture has been eroded by several reports of a single fig species hosting multiple pollinator species and a few reports of two fig species sharing pollinator species -.
It is now clear that the paradigm of one-to-one reciprocal partner specificity is no longer tenable ,, but true patterns of partner species associations remain to be revealed. Recent reviews have posited that perhaps a third or half of the >750 fig species worldwide may have multiple pollinator species ,. However, these estimates are best regarded as informed guesses for two main reasons. First, most of the evidence for multiple pollinators is the by-product of studies conducted for other reasons, such as exploring wasp reproductive behaviour , or phylogeographic histories . Thus, these lines of evidence do not come from targeted studies aiming to document the numbers of species involved in the plant/pollinator symbiosis. In particular, there is a striking lack of studies involving wide sampling of insects from across host plant geographic ranges ,.
Second, there is a ¿taxonomy gap¿ that hinders assessment of species associations -. There are only about 150 described agaonid species, but the true number seems likely to exceed 1000 . This means that investigation of the pollinators of many fig species begins with only a genus level identification and with no described wasp species already linked to the fig species in question. Moreover, even when a described pollinator is known, genetic studies often reveal further pollinator species that are either morphologically cryptic within the one named entity, morphologically distinguishable but not previously sampled, or previously sampled but unrecognised within mixed species collections, e.g. .
Molecular data have already contributed substantially to rejection of the old 1:1 paradigm of fig/pollinator specificity and should now play a key role, in tandem with morphological analysis, in establishing true patterns of pollinator diversity and variation in plant/pollinator interactions at local and regional scales . Wide geographic sampling is crucial , because: a) some species may not occur in all parts of the host range, and b) if intraspecific genetic variation is underestimated by sampling few sites, it may be harder to identify the molecular `barcoding gap¿ (or appropriate clades in phylogeny-based methods) between species ,.
Another important issue is the choice of markers for molecular taxonomy. The animal barcoding approach developed by CBOL uses a standard section of the mitochondrial COI gene ,. This has proved valuable in some studies of fig wasps, but in a recent phylogenetic study of 200 species from 19 genera , the success rate for PCR and direct sequencing was much higher with cytochrome b and these two linked mitochondrial markers tend to reveal very similar patterns . Both markers can vary substantially within fig wasp (and other) species, which may hinder correct species delimitation, in particular via over-splitting ,. Consequently, it is desirable to also use a nuclear marker and to seek congruent species delimitation between the two genomes . With nuclear markers, the challenge is to find one that both amplifies reliably across species and shows sufficient variation for discrimination between closely related species . None of the currently used nuclear markers achieve the widespread utility of favoured mtDNA markers like COI and cytb. However, an internal transcribed spacer region (ITS2) of rDNA has recently been proposed as a useful nuclear barcoding marker for animals -.
A further issue with molecular data is the choice of species delimitation method. Barcoding-type approaches use genetic distance data directly to identify a barcoding gap between the pairwise genetic distances found within and between species. These approaches often work well in practice, but threshold genetic distances vary across taxa, and well-known biological phenomena, such as introgression and selective sweeps due to Wolbachia bacteria, can confound the expected patterns. Meanwhile, phylogeny-based methods aim to identify clades that are evolutionary significant units (ESUs) and invoke the phylogenetic species concept . These have a justifiable conceptual basis, but are yet to be as widely used and are experiencing a period of relatively rapid methods development.
In this study, we investigate the diversity of pollinator wasps associated with a single widespread fig species, Ficus rubiginosa. This plant is endemic to Australia and occurs widely in diverse habitats, including eucalypt scrub and rainforest, in a roughly 2500 km coastal belt that stretches from tropical northern Queensland to temperate southern New South Wales . It belongs to the Ficus section Malvanthera, which is pollinated by wasp species in the genus Pleistodontes. Its only known pollinator species is P. imperialis, which was originally described by Saunders (1882). Following new wasp collections from several malvantheran fig species, a taxonomic revision of Pleistodontes was carried out . This led to the description of seven new Pleistodontes species overall, but no change to the conclusion that F. rubiginosa was associated with a single pollinator species (P. imperialis).
Compare the 'performance of mitochondrial (cytb and COI) and nuclear (ITS2) molecular markers, and different species delimitation methods (distance and phylogeny-based) for identifying species.
Uncover the full diversity of pollinator species associated with F. rubiginosa and the impact of sampling effort on their detection.
Explore geographic turnover (beta diversity) of pollinator species across the host plant¿s wide natural latitudinal range.
Bayesian analysis on a subset (n = 44) of individuals for the 416 bp COI region suggests six major clades (Additional file 2: Figure S2), comprising the five major cytb clades but with species 1 split into northern and southern clades. However, support values for clades 2¿5 are all 1.00 whilst the two species 1 sub-clades have support values of 0.53 and 0.80 indicating poor support for a six species hypothesis. Moreover, jMOTU analysis of COI indicates five species (Additional file 3: Figure S3). Visual investigation of pairwise distances (Additional file 4: Figure S4) also does not reveal a clear barcode gap in COI data, with intraspecific divergences of 0¿5.6% overlapping slightly with interspecific divergences (5.0-13.1%). Sub-clade congruence is absolute for individuals sequenced for both cytb and COI.
Symbioses have been implicated in the innovation of key evolutionary leaps, the development of a number of fundamental ecosystem functions, and as major drivers in the generation of biodiversity. However, the exact nature of the ecological relationships between symbiotic partners is often unknown due to the lack of targeted investigations, taxonomic impediments and insufficient sampling. Among classic examples of symbioses, the fig-wasp system provides a crucial ecological resource across tropical and savannah environments. Here we provide the first comprehensive targeted investigation to assess pollinator diversity across the geographic range of a common and widespread fig species. Our study supports the existence of at least five wasp species pollinating Ficus rubiginosa. This P. imperialis complex involves more pollinator species than any other fig species studied to date , but see , and our study adds to the growing body of evidence that the diversity of fig wasps, and insects in general, has been greatly underestimated ,,.
Our most comprehensive sampling used mitochondrial cytb and delineated up to 11 ESUs that might each constitute a species (Figure 1). This cytb diversity is far greater than previously identified and is most likely explained by a big increase in sampling effort relative to that of Haine et al. , whose analyses highlighted four largely unstructured clades. Most striking among this new diversity are the six individuals ascribed to species 5 that are distinguished from other species by extended branch lengths in the phylogeny. These most likely represent a rare fifth pollinator species not sampled previously. Moreover, we found strong sub-clade structure within some of the putative species and this generally corresponds to largely or wholly Northern versus Southern populations of the more widespread species. Such geographic isolation is suggestive of restricted gene flow and could therefore indicate distinct species, but further results argue against this.
Investigation of modelled pairwise genetic distances for cytb data suggests that species delimitation is not straightforward among these wasps due to the lack of a barcoding gap. However, jMOTU analyses of cytb data clearly supported the existence of five cryptic species and this algorithm identified a barcoding gap between 17¿30 base pairs. Additionally, GMYC analyses on cytb data give strong statistical support for the existence of five P. imperialis pollinators. This method adds valuable weight to the five species hypothesis as it statistically tests the crossover point between intra- and interspecific branch lengths derived from genetic markers.
In order to compare cytb phylogeny with the standard cox1 barcoding gene, a shortened COI fragment of only 416 bp cf 658 bp;  was employed on a subset of 44 individuals. The COI topology is very similar to that for cytb but splits species 1 into two. However, the very low support values (p = 0.53 and p = 0.80) offer no strong evidence that either of the species 1 sub-clades are more closely related to species 2 for this marker. Moreover, jMOTU analyses support five species for this marker in the same way as for cytb. The nuclear ITS2 phylogeny is similar to that for cytb in revealing the same major ESUs. However, sub-clade structure is absent for this marker and there is essentially zero variation within any of the five species, allowing straightforward delimitation of taxa. Interestingly, this means that ITS2 behaves far more like a barcode in the wider general sense of this word than mtDNA markers like COI or cytb, which typically show considerable within species variation even when there is a clear barcoding gap. This disparity in resolution between mtDNA and nDNA is most likely explained by different population genetic properties of these markers , and by multigene concerted evolution of the repeated rDNA cluster.
Our sampling reveals that each of the five species has a different geographic range (Figure 5). Whilst there is no absolute geographical replacement of species, four species either dominate or are codominant in particular regions. Species 3 is dominant in northern Queensland, while species 4 appears dominant in central and southern Queensland and in the inland region around Forty Mile Scrub. The yellow species 2 is found only in the Townsville region, where it is common and co-occurs with species 3 at similar frequencies. Finally, species 1 occupies widely disjunct regions in the far north of the host plant range and in New South Wales, where it is the only pollinator species recorded. Moreover, species 1 appears excluded from most of the intermediate regions as no individuals were recorded amongst 220 wasps from the Townsville and Forty Mile Scrub populations. Only species 5 appears to have no geographical stronghold, with only six individuals sampled overall. However, this may be an artifact of the low sample size for this rare species.
From the perspective of understanding the ecological dynamics of symbiosis this study clearly rejects the simple model of reciprocally partner-specific fig and pollinator species. Over most of the host plant range, it is associated with more than one pollinator species, but the frequency of the different species varies considerably between regions. Given the wide latitudinal range of the host plant, from the wet tropics to the temperate zone, it is possible that pollinator diversity is related to local adaptation to different climates. Whatever the causes of variation in pollinator identity and diversity between regions, the variation itself has important implications for insect/plant coevolution . First, the association appears to be a 1:1 match of a fig and pollinator species in the southern part of the range in NSW, but to typically involve two or three (but not always the same) co-occurring pollinator species at sites further North in Queensland. This could lead to different coevolutionary trajectories in different parts of the range, e.g. because multiple symbiont species may increase host/symbiont conflict as a result of antagonistic competition between the symbiont species .
Another set of questions arising from our improved understanding of pollinator species diversity and distributions surrounds the local coexistence of multiple pollinator species utilizing the same host resources. Ecological theory suggests that this is difficult to achieve without some degree of niche separation, although an intriguing possibility with fig wasps is that it might occur through density-dependent sex ratio (and therefore population growth rate) variation . Different sites have different sets of pollinators (Figure 5), creating the possibility to compare different coexistence patterns in different but overlapping sets of species.
There has been only limited investigation of the comparative ecology of multiple sympatric pollinator species associated with a given fig species ,-, and most such studies have not revealed clear differences. Two exceptions are cases where one of the pollinators is a ¿cheat¿ that does not actively pollinate its host plant ,, and cases where one species is a diurnal, and the other a nocturnal, disperser e.g. . The latter case may be relevant in the P. imperialis complex because species two is yellow rather than dark brown or black, and previous studies in Africa and Asia have shown that pale coloured wasps tend to be nocturnal dispersers, while dark ones tend to be diurnal e.g. . This is one potential axis for niche differentiation, but four of the five P. imperials complex species are brown/black wasps so are unlikely to be differentiated in this way. Further research into host use by the P. imperialis species complex could focus on genetic and morphological variation within F. rubiginosa.
Another key message from our study is the importance of a comprehensive sampling regime. Haine et al.  sampled 71 P. imperialis from across F. rubiginosa¿s range and hypothesised the existence of four cryptic species. Our sampling covered much of the same geographic region in a more thorough manner, as well as some additional locations. Despite sampling 415 wasps, we only captured six individuals of the newly identified species 5 from three coastal sites up to 1000 km apart in central and southern Queensland. It remains possible that further `P. imperialis¿ species are yet to be sampled. The sampling effort from NSW is much less comprehensive than that of Queensland, but this so far appears justified as only a single species has been found from 20 specimens from several sites.
Our study offers comprehensive evidence that mutualistic P. imperialis wasps pollinating the fig species F. rubiginosa have diversified into five distinct (including some cryptic) species across their entire host plant range along the east coast of Australia. Given the absence of diagnostic morphological and ecological differences in this sister-species complex that may offer corroborating evidence for the assessment of species status in an integrative taxonomic framework, we have shown the utility of the internal transcribed spacer region (ITS2) of rDNA as a diagnostic tool for species delimitation among problematic taxa. Furthermore, we have shown that four out of five of these identified taxa show strong patterns of regional geographic dominance that suggests an adaptive origin to localised ecological conditions and offers a likely mechanism in mediating their co-existence performing the same ecological role in the same fig species. Moreover, where we identify a pattern of geographic co-dominance there is reason to believe that variation in pigmentation may indicate niche differentiation in diurnal/nocturnal dispersal patterns. In addition to this interspecific geographic structuring, our comprehensive sampling also suggests geographic structure at the intraspecific level for three of the species. Thus, our results show that one-to-one specificity may often break down in symbioses and offer an explanation as to under what circumstances it may do so. Additionally, the increasing level of pollinator diversity at lower latitudes implies that coevolutionary trajectories between host and symbiont will vary across the range. In summary, our understanding of the structuring of biodiversity in a symbiotic context is likely to be simplistic and will require targeted studies using molecular taxonomic techniques and comprehensive sampling schemes if we are to better understand the coevolutionary dynamics underpinning intimate ecological interactions.
Most sampling was conducted in Queensland from 2007¿2009 along the eastern seaboard and immediate hinterland between Brisbane (26° 46S, 153°02E) in the south and Dimbulah (17° 01S, 145°19E) in the north. Sporadic sampling was also undertaken in New South Wales (NSW), and also from some planted trees outside their natural range in Victoria and South Australia between 2000 and 2010. Near-ripe figs were collected from trees and placed into hatching jars with meshed lids that allowed air flow, prevented overheating, and prevented wasp escape. After 48 h each fig and all its emerged wasps were placed into 70% ethanol. Alternatively, figs were placed directly into alcohol and wasps were dissected out at a later date. Since most pollinator wasps developing in a given fig are siblings, we used only one wasp per syconium for DNA extraction to maximise the independence of samples. Sampling intensity was greater in Queensland than in New South Wales, but evidence for much lower diversity in New South Wales both before  and after this study (T. Sutton, pers. comm.) justifies the lower sampling effort in NSW. Additional file 6: Table S2 gives location, markers sequenced and GenBank accession numbers for all 415 wasps used in this study.
A Chelex method was used for DNA extraction  and two mitochondrial (cytb and COI) and one nuclear (ITS2) marker were amplified. NB the COI marker is a region of cox1 shorter than the standard barcoding region, where alternate primers have been designed to avoid nuclear pseudo-genes within certain taxa: . Our widest genetic sampling of individuals focussed on cytochrome b (cytb) using mtDNA primers CB1 and CB2 . A 396 bp fragment of cytb was amplified for 415 wasps in order to assign individuals to the putative species identified by Haine et al. . Cytb was chosen as our primary marker in order to integrate our study with the earlier one by Haine et al. . We also analysed a subset of individuals using the COI region to facilitate comparison between mtDNA genes and connectivity with other datasets.
Subsets of 3¿15 individuals per cytb clade were sequenced for ITS2 (n = 54)  and a 416 bp fragment of COI (n = 44) (primers CI-J-1751 and CI-N-2191; see ) in order to clarify and confirm species delimitation. Cytb was amplified using a Techne Touchgene gradient machine with 3 min at 94°C, 30 cycles of 15 s at 95°C, 20 s at 45°C, 30 s at 72°C, and a final elongation step of 10 min at 72°C. Amplification of COI consisted of 5 min at 94°C, 30 cycles of 30 s at 94°C, 45 s at 50°C, 60 s at 72°C then 10 min at 72°C. Amplification of ITS2 used 5 min at 94°C, 35 cycles of 30 s at 94°C, 40 s at 55°C, 40 s at 72°C then 10 min at 72°C. Subsequent purification and sequencing reactions were conducted by Macrogen Inc. Purification was performed using ethanol precipitation and sequencing by BigDyeTM terminator cycling conditions and a 3730xl DNA analyser. All sample sequences have been deposited in GenBank: cytb [GenBank: KM249475-KM249848], COI [GenBank: KM249375 - KM249419], ITS2 [GenBank: KM249420 - KM249474]; see Additional file 6: Table S2.
Sequence alignment and phylogenetic analysis
Sequence quality was checked using Finch TV Version 1.4.0 and sequences were edited and aligned using BioEdit . In order to place collected samples into the clade (hereafter species) categorisations of Haine et al. , and as a template with which to align newly generated sequences, 23 P. imperialis cytb sequences with geographical location data were downloaded from GenBank/EMBI (from accession numbers AJ298439 and AY567594 - AY567638). ITS2 sequences were checked for microsatellites  and a small region of 10-16 bp `TC¿ repeats was identified. Analyses were run with and without this region, which caused no differences in phylogenetic inference, other than a 2% and 3% reduction in support values at two internal nodes. P. imperialis ITS2 sequence lengths were between 312-317 bp due to the presence of indels. Indels ranged between 0-4 bp within individual species with no indels found in species 2. Within species 1, a 2 bp indel was diagnostic of specimens from Queensland and New South Wales clades. No obvious pattern was evident among the distribution of indels for other species. Additionally, ITS2 sequence length was 474 bp in the outgroup Pleistodontes sp from F. glandifera and created indels throughout the alignment. No evidence of pseudogenes or heteroplasmy was noted for mitochondrial regions. Bayesian methods were employed to construct phylogenies using MrBayes . The model of nucleotide substitution for each gene was chosen using jModelTest2 . For mitochondrial markers the TrN + I + G and TPM1uf + G models were chosen for cytb and COI respectively whilst the TPM3uf + G model was chosen for ITS2.
Species delimitation methods
We used two main strategies, distance (barcoding-type) and phylogeny-based, to delimit species. We applied two distance methods to cytb and COI data: (i) we calculated pairwise genetic distances in PAUP  and plotted them in order to determine visually whether a barcode gap was present; (ii) we used jMOTU software  to identify the inflection point in the frequency distribution of genetic distances, which is purported to show the barcoding gap. For jMOTU analyses, only sequences without ambiguous base calls were utilized.
As a complementary phylogeny-based analysis on cytb, we identified evolutionary significant units (ESUs) by applying the statistical GMYC method of Pons et al.  to the cytb phylogeny, using the R-Splits package  in R . The GMYC method requires an ultrametric phylogenetic tree and this was created from cytb data with BEAST software . Models of evolution implemented in BEAST were chosen using Path Sampling and Stepping Stone techniques . These allow the calculation of Bayes factors, which were evaluated using the criteria of Raftery , and the HKY + G model was chosen. Two independent runs of 120 million generations were performed, under constant-clock conditions with a constant coalescent model of species evolution. Genealogies and model parameters were sampled once every 6000 iterations. Preliminary runs of the GMYC using only P. imperialis sequences yielded non-significant results due to a lack of statistical power, probably due to the limited number of species despite large numbers of individuals. Consequently, we added sequences from four other Pleistodontes species and this led to significant results.
Finally, we compared the geographic distributions of the pollinator species across eight regional populations. We used Chi-square tests in R  to test null hypotheses that each species occurs at the same frequency in each region.
Availability of supporting data
Sequence data has been deposited on Genbank whilst population frequencies are in Supporting Information.
CTD: conducted fieldwork collections, molecular genetic work, sequence alignments, phylogenetic inference, statistical analyses and drafted the manuscript. SAB: conducted fieldwork collections, molecular genetic work and sequence alignments. JMC: conducted fieldwork collections, drafted the manuscript and conceived the study. All authors read and approved the final manuscript.
We would like to thank Brent Emerson and Simon van Noort for helpful comments regarding this manuscript, Simon Segar and Derek Dunn for helping with data collection, and Jean-Yves Rasplus for taxonomic advice. This work was funded by a NERC studentship (#NE/G523912/1 - CTD). There was no involvement of the funding sources in the study design; in the collection, analysis and interpretation of data; in the writing of the report; and in the decision to submit the article for publication.
- Margulis L: Origin of Eukaryotic Cells. 1970, Yale University Press, Newhaven, ConnecticutGoogle Scholar
- Kiers ET, Rousseau RA, West SA, Denison RF: Host sanctions and the legume¿rhizobium mutualism. Nature. 2003, 425: 78-81. 10.1038/nature01931.PubMedView ArticleGoogle Scholar
- Darwin C: On the Origin of Species by Means of Natural Selection, or the Preservation of Favoured Races in the Struggle for Life. 1859, John Murray, London, UKView ArticleGoogle Scholar
- Leigh EG: The evolution of mutualism. J Evol Biol. 2010, 23 (12): 2507-2528. 10.1111/j.1420-9101.2010.02114.x.PubMedView ArticleGoogle Scholar
- Althoff DM, Segraves KA, Johnson MTJ: Testing for coevolutionary diversification: linking pattern with process. Trends Ecol Evol. 2014, 29 (2): 82-89. 10.1016/j.tree.2013.11.003.PubMedView ArticleGoogle Scholar
- Thrall PH, Hochberg ME, Burdon JJ, Bever JD: Coevolution of symbiotic mutualists and parasites in a community context. Trends Ecol Evol. 2007, 22 (3): 120-126. 10.1016/j.tree.2006.11.007.PubMedView ArticleGoogle Scholar
- Rooney N, McCann KS: Integrating food web diversity, structure and stability. Trends Ecol Evol. 2012, 27 (1): 40-46. 10.1016/j.tree.2011.09.001.PubMedView ArticleGoogle Scholar
- Herre EA, Knowlton N, Mueller UG, Rehner SA: The evolution of mutualisms: exploring the paths between conflict and cooperation. Trends Ecol Evol. 1999, 14 (2): 49-53. 10.1016/S0169-5347(98)01529-8.PubMedView ArticleGoogle Scholar
- Cook JM, Rasplus JY: Mutualists with attitude: coevolving fig wasps and figs. Trends Ecol Evol. 2003, 18 (7): 241-248. 10.1016/S0169-5347(03)00062-4.View ArticleGoogle Scholar
- Herre EA, Jander KC, Machado CA: Evolutionary ecology of figs and their associates: recent progress and outstanding puzzles. Annu Rev Ecol Evol Syst. 2008, 39: 439-458. 10.1146/annurev.ecolsys.37.091305.110232.View ArticleGoogle Scholar
- Segar ST, Cook JM: The dominant exploiters of the fig/pollinator mutualism vary across continents, but their costs fall consistently on the male reproductive function of figs. Ecol Entomol. 2012, 37: 342-349. 10.1111/j.1365-2311.2012.01370.x.View ArticleGoogle Scholar
- Kissling WD, Rahbek C, Bohning-Gaese K: Food plant diversity as broad-scale determinant of avian frugivore richness. Proc R Soc B-Biol Sci. 2007, 274 (1611): 799-808. 10.1098/rspb.2006.0311.View ArticleGoogle Scholar
- Lopez-Vaamonde C, Rasplus JY, Weiblen GD, Cook JM: Molecular phylogenies of fig wasps: partial cocladogenesis of pollinators and parasites. Mol Phylogenet Evol. 2001, 21 (1): 55-71. 10.1006/mpev.2001.0993.PubMedView ArticleGoogle Scholar
- Weiblen GD: How to be a fig wasp. Annu Rev Entomol. 2002, 47: 299-330. 10.1146/annurev.ento.47.091201.145213.PubMedView ArticleGoogle Scholar
- Molbo D, Machado CA, Sevenster JG, Keller L, Herre EA: 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 U S A. 2003, 100 (10): 5867-5872. 10.1073/pnas.0930903100.PubMedPubMed CentralView ArticleGoogle Scholar
- Molbo D, Machado CA, Herre EA, Keller L: Inbreeding and population structure in two pairs of cryptic fig wasp species. Mol Ecol. 2004, 13 (6): 1613-1623. 10.1111/j.1365-294X.2004.02158.x.PubMedView ArticleGoogle Scholar
- Haine ER, Martin J, Cook JM: Deep mtDNA divergences indicate cryptic species in a fig-pollinating wasp. BMC Evol Biol. 2006, 6: 83-10.1186/1471-2148-6-83.PubMedPubMed CentralView ArticleGoogle Scholar
- Chen Y, Compton SG, Liu M, Chen X-Y: Fig trees at the northern limit of their range: the distributions of cryptic pollinators indicate multiple glacial refugia. Mol Ecol. 2012, 21 (7): 1687-1701. 10.1111/j.1365-294X.2012.05491.x.PubMedView ArticleGoogle Scholar
- Cornille A, Underhill JG, Cruaud A, Hossaert-McKey M, Johnson SD, Tolley KA, Kjellberg F, van Noort S, Proffit M: Floral volatiles, pollinator sharing and diversification in the fig-wasp mutualism: insights from Ficus natalensis, and its two wasp pollinators (South Africa). Proc Roy Soc B: Biol Sci. 2012, 279: 1731-1739. 10.1098/rspb.2011.1972.View ArticleGoogle Scholar
- McLeish MJ, van Noort S: Codivergence and multiple host species use by fig wasp populations of the Ficus pollination mutualism.BMC Evol Biol, 12:1¿2012.,Google Scholar
- Compton SG, Grehan K, van Noort S: A fig crop pollinated by three or more species of agaonid fig wasps. African Entomol. 2009, 17 (2): 215-222. 10.4001/003.017.0212.View ArticleGoogle Scholar
- Machado CA, Robbins N, Gilbert MTP, Herre EA: Critical review of host specificity and its coevolutionary implications in the fig/fig-wasp mutualism. Proc Natl Acad Sci U S A. 2005, 102: 6558-6565. 10.1073/pnas.0501840102.PubMedPubMed CentralView ArticleGoogle Scholar
- Cook JM, Segar ST: Speciation in fig wasps. Ecol Entomol. 2010, 35: 54-66. 10.1111/j.1365-2311.2009.01148.x.View ArticleGoogle Scholar
- Kobmoo N, Hossaert-Mckey M, Rasplus JY, Kjellberg F: Ficus racemosa is pollinated by a single population of a single agaonid wasp species in continental South-East Asia. Mol Ecol. 2010, 19: 2700-2712. 10.1111/j.1365-294X.2010.04654.x.PubMedView ArticleGoogle Scholar
- Kaartinen R, Stone GN, Hearn J, Lohse K, Roslin T: Revealing secret liaisons: DNA barcoding changes our understanding of food webs. Ecol Entomol. 2010, 35 (5): 623-638. 10.1111/j.1365-2311.2010.01224.x.View ArticleGoogle Scholar
- Smith MA, Eveleigh ES, McCann KS, Merilo MT, McCarthy PC, Van Rooyen KI: Barcoding a quantified food web: crypsis, concepts, ecology and hypotheses. PLoS One. 2011, 6 (7): e14424-10.1371/journal.pone.0014424.PubMedPubMed CentralView ArticleGoogle Scholar
- Wheeler QD, Raven PH, Wilson EO: Taxonomy: impediment or expedient?. Sci (New York, NY). 2004, 303 (5656): 285-10.1126/science.303.5656.285.View ArticleGoogle Scholar
- Cook JM, West SA: Figs and fig wasps. Curr Biol. 2005, 15 (24): R978-R980. 10.1016/j.cub.2005.11.057.PubMedView ArticleGoogle Scholar
- Lopez-Vaamonde C, Dixon DJ, Cook JM, Rasplus JY: Revision of the Australian species of pleistodontes (Hymenoptera : Agaonidae) fig-pollinating wasps and their host-plant associations. Zool J Linn Soc. 2002, 136 (4): 637-683. 10.1046/j.1096-3642.2002.00040.x.View ArticleGoogle Scholar
- Vogler AP, Monaghan MT: Recent advances in DNA taxonomy. J Zool Syst Evol Res. 2007, 45 (1): 1-10. 10.1111/j.1439-0469.2006.00384.x.View ArticleGoogle Scholar
- Moritz C, Cicero C: DNA barcoding: promise and pitfalls. Plos Biol. 2004, 2 (10): 1529-1531. 10.1371/journal.pbio.0020354.View ArticleGoogle Scholar
- Will KW, Rubinoff D: Myth of the molecule: DNA barcodes for species cannot replace morphology for identification and classification. Cladistics- Int J Willi Hennig Soc. 2004, 20 (1): 47-55. 10.1111/j.1096-0031.2003.00008.x.View ArticleGoogle Scholar
- Hebert PDN, Ratnasingham S, deWaard JR: ¿Barcoding animal life: cytochrome c oxidase subunit 1 divergences among closely related species,¿. Proc R Soc London Ser B-Biological Sci. 2003, 270: S96-S99. 10.1098/rsbl.2003.0025.View ArticleGoogle Scholar
- Hebert PDN, Cywinska A, Ball SL, DeWaard JR: Biological identifications through DNA barcodes. Proc R Soc London Ser B-Biological Sci. 2003, 270 (1512): 313-321. 10.1098/rspb.2002.2218.View ArticleGoogle Scholar
- Cruaud A, Jabbour-Zahab R, Genson G, Couloux A, Peng Y-Q, Da Rong Y, Ubaidillah R, Santinelo Pereira RA, Kjellberg F, van Noort S, Kerdelhue C, Rasplus J-Y: Out of Australia and back again: the world-wide historical biogeography of non-pollinating fig wasps (Hymenoptera: Sycophaginae). J Biogeogr. 2011, 38 (2): 209-225. 10.1111/j.1365-2699.2010.02429.x.View ArticleGoogle Scholar
- Lopez-Vaamonde C, Wikstrom N, Kjer KM, Weiblen GD, Rasplus JY, Machado CA, Cook JM: Molecular dating and biogeography of fig-pollinating wasps. Mol Phylogenet Evol. 2009, 52 (3): 715-726. 10.1016/j.ympev.2009.05.028.PubMedView ArticleGoogle Scholar
- Avise J: Phylogeography: The History And Formation Of Species. 2000, Harvard University Press, Cambridge Massachusetts, London EnglandGoogle Scholar
- Wiemers M, Fiedler K: Does the DNA barcoding gap exist? - a case study in blue butterflies (Lepidoptera: Lycaenidae). Front Zool. 2007, 4: 8-10.1186/1742-9994-4-8.PubMedPubMed CentralView ArticleGoogle Scholar
- Hurst GDD, Jiggins FM: Problems with mitochondrial DNA as a marker in population, phylogeographic and phylogenetic studies: the effects of inherited symbionts. Proc R Soc B-Biol Sci. 2005, 272 (1572): 1525-1534. 10.1098/rspb.2005.3056.View ArticleGoogle Scholar
- Lohse K, Sharanowski B, Blaxter M, Nicholls JA, Stone GN: Developing EPIC markers for chalcidoid Hymenoptera from EST and genomic data. Mol Ecol Resour. 2011, 11 (3): 521-529. 10.1111/j.1755-0998.2010.02956.x.PubMedView ArticleGoogle Scholar
- Coleman AW: Is there a molecular key to the level of `biological species¿ in eukaryotes? A DNA guide. Mol Phylogenet Evol. 2009, 50 (1): 197-203. 10.1016/j.ympev.2008.10.008.PubMedView ArticleGoogle Scholar
- Li YW, Zhou X, Feng G, Hu HY, Niu LM, Hebert PDN, Huang DW: COI and ITS2 sequences delimit species, reveal cryptic taxa and host specificity of fig-associated Sycophila (Hymenoptera, Eurytomidae). Mol Ecol Resour. 2010, 10 (1): 31-40. 10.1111/j.1755-0998.2009.02671.x.PubMedView ArticleGoogle Scholar
- Muller T, Philippi N, Dandekar T, Schultz J, Wolf M: Distinguishing species. Rna-a Publ Rna Soc. 2007, 13 (9): 1469-1472. 10.1261/rna.617107.View ArticleGoogle Scholar
- Zhou MJ, Xiao JH, Bian SN, Li YW, Niu LM, Hu HY, Wu WS, Murphy RW, Huang DW: Molecular approaches identify known species, reveal cryptic species and verify host specificity of Chinese Philotrypesis (Hymenoptera: Pteromalidae). Mol Ecol Resour. 2012, 12 (4): 598-606. 10.1111/j.1755-0998.2012.03127.x.PubMedView ArticleGoogle Scholar
- Eldredge N, Cracraft J: Phylogenetic Patterns And The Evolutionary Process Method And Theory In Comparative Biology. Phylogenetic Patterns and the Evolutionary Process: Method and Theory in Comparative Biology. Viii+349p. Edited by: Eldredge N, Cracraft J. 1980, Columbia University Press, New York, N.Y., USA. Illus. MapsGoogle Scholar
- Dixon DJ, Jackes BR, Bielig LM: Figuring out the figs: the ficus obliqua-ficus rubiginosa complex (Moraceae : Urostigma sect. Malvanthera). Aust Syst Bot. 2001, 14 (1): 133-154. 10.1071/SB99029.View ArticleGoogle Scholar
- Lopez-Vaamonde C, Dixon DJ, Cook JM, Rasplus JY: Revision of the Australian species of pleistodontes (Hymenoptera : Agaonidae) fig-pollinating wasps and their host-plant associations. Zool J Linn Soc. 2002, 136 (4): 637-683. 10.1046/j.1096-3642.2002.00040.x.View ArticleGoogle Scholar
- De Queiroz K: Species concepts and species delimitation. Syst Biol. 2007, 56 (6): 879-886. 10.1080/10635150701701083.PubMedView ArticleGoogle Scholar
- Moe AM, Weiblen GD: Molecular divergence in allopatric ceratosolen (Agaonidae) pollinators of geographically widespread ficus (Moraceae) species. Ann Entomol Soc Am. 2010, 103 (6): 1025-1037. 10.1603/AN10083.View ArticleGoogle Scholar
- Bruno JF, Cardinale BJ: Cascading effects of predator richness. Front Ecol Environ. 2008, 6 (10): 539-546. 10.1890/070136.View ArticleGoogle Scholar
- Lowe Harris A, S Ashton P: Ecological Genetics. 2004, Blackwell Science limited, Carlton, Victoria, AustraliaGoogle Scholar
- Thompson JN: The Geographic Mosaic of Coevolution. 2005, The University of Chicago Press, ChicagoGoogle Scholar
- Zhang DY, Lin K, Hanski I: Coexistence of cryptic species. Ecol Lett. 2004, 7 (3): 165-169. 10.1111/j.1461-0248.2004.00569.x.View ArticleGoogle Scholar
- Weiblen GD, Bush GL: Speciation in fig pollinators and parasites. Mol Ecol. 2002, 11 (8): 1573-1578. 10.1046/j.1365-294X.2002.01529.x.PubMedView ArticleGoogle Scholar
- Compton SG, Holton KC, Rashbrook VK, Vannoort S, Vincent SL, Ware AB: Studies of Ceratosolen galili, a non-pollinating agaonid fig wasp. Biotropica. 1991, 23 (2): 188-194. 10.2307/2388305.View ArticleGoogle Scholar
- Peng Y-Q, Duan Z-B, Yang D-R, Rasplus J-Y: Co-occurrence of two Eupristina species on Ficus altissima in Xishuangbanna, SW China. Symbiosis. 2008, 45 (1¿3): 9-14.Google Scholar
- Warren M, Robertson MP, Greeff JM: A comparative approach to understanding factors limiting abundance patterns and distributions in a fig tree-fig wasp mutualism. Ecography (Cop). 2010, 33 (1): 148-158. 10.1111/j.1600-0587.2009.06041.x.View ArticleGoogle Scholar
- West SA, Cook JM, Werren JH, Godfray HCJ: Wolbachia in two insect host-parasitoid communities. Mol Ecol. 1998, 7 (11): 1457-1465. 10.1046/j.1365-294x.1998.00467.x.PubMedView ArticleGoogle Scholar
- Xiao J-H, Wang N-X, Li Y-W, Murphy RW, Wan D-G, Niu L-M, Hu H-Y, Fu Y-G, Huang D-W: Molecular approaches to identify cryptic species and polymorphic species within a complex community of fig wasps. PLoS One. 2010, 5 (11): e15067-10.1371/journal.pone.0015067.PubMedPubMed CentralView ArticleGoogle Scholar
- Jermiin LS, Crozier RH: The cytochrome-B region in the mitochondrial DNA of the ant Tetraponera rufoniger - sequence divergence in hymenoptera may be associated with nucleotide content. J Mol Evol. 1994, 38 (3): 282-294. 10.1007/BF00176090.PubMedView ArticleGoogle Scholar
- White T, Bruns T, Lee S, Taylor J: Amplification And Direct Sequencing Of Fungal Ribosomal Genes For Phylogenies. 1990, Academic, New York, N.Y.Google Scholar
- Simon C, Frati F, Beckenbach A, Crespi B, Liu H, Flook P: Evolution, weighting, and phylogenetic utility of mitochondrial gene-sequences and a compilation of conserved polymerase chain-reaction primers. Ann Entomol Soc Am. 1994, 87 (6): 651-701.View ArticleGoogle Scholar
- Hall TA: BioEdit: a user friendly biological sequence alignment editor and analysis program for Windows95/98NT. Nucleic Acids Symp. 1999, 41: 95-98.Google Scholar
- Harris DJ, Crandall KA: Intragenomic variation within ITS1 and ITS2 of freshwater crayfishes (Decapoda: Cambaridae): implications for phylogenetic and microsatellite studies. Mol Biol Evol. 2000, 17: 284-291. 10.1093/oxfordjournals.molbev.a026308.PubMedView ArticleGoogle Scholar
- Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19 (12): 1572-1574. 10.1093/bioinformatics/btg180.PubMedView ArticleGoogle Scholar
- Darriba D, Taboada GL, Doallo R, Posada D: jModelTest 2: more models, new heuristics and parallel computing. Nat Methods. 2012, 9 (8): 772-10.1038/nmeth.2109.PubMedPubMed CentralView ArticleGoogle Scholar
- Swofford DL: PAUP* Phylogenetic Analysis Using Parsimony (*and Other Methods). 2002, Sinauer Associates, Sunderland, Massachusetts, USAGoogle Scholar
- Jones M, Ghoorah A, Blaxter M: jMOTU and taxonerator: turning DNA barcode sequences into annotated operational taxonomic units. PLoS One. 2011, 6 (4): 10-10.1371/journal.pone.0019259.View ArticleGoogle Scholar
- Pons J, Barraclough TG, Gomez-Zurita J, Cardoso A, Duran DP, Hazell S, Kamoun S, Sumlin WD, Vogler AP: Sequence-based species delimitation for the DNA taxonomy of undescribed insects. Syst Biol. 2006, 55 (4): 595-609. 10.1080/10635150600852011.PubMedView ArticleGoogle Scholar
- Ezard T, Fujisawa T, Barraclough T: Species¿ Limits by Threshold Statistics. 2009Google Scholar
- R: A language and environment for statistical computing. R Found Stat Comput. 2011, 1: 409-Google Scholar
- Drummond AJ, Rambaut A: BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol Biol. 2007, 7: 214-10.1186/1471-2148-7-214.PubMedPubMed CentralView ArticleGoogle Scholar
- Baele G, Lemey P, Rambaut A, Suchard MA, Alekseyenko AV: Improving the accuracy of demographic and molecular clock model comparison while accommodating phylogenetic uncertainty.Mol Biol Evol 2012, In press.,Google Scholar
- Raftery AE: Hypothesis testing and model selection. Markov chain Monte Carlo in practice. Edited by: Gilks WR, Richardson S, Spiegelhalter DJ. 1996, Chapman & Hall, LondonGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.