<figure><img src="image_1.png" /></figure>

HAL
open science

# Diversification and spatial structuring in the mutualism between Ficus septica and its pollinating wasps in insular South East Asia

Lillian Jennifer Rodriguez, Anthony Bain, Lien-Siang Chou, Lucie Conchou, Astrid Cruaud, Regielene Gonzales, Martine Hossaert-Mckey, Jean-Yves Rasplus, Hsy-Yu Tzeng, Finn Kjellberg

<figure><img src="image_2.png" /></figure>

## To cite this version:

Lillian Jennifer Rodriguez, Anthony Bain, Lien-Siang Chou, Lucie Conchou, Astrid Cruaud, et al. "Diversification and spatial structuring in the mutualism between Ficus septica and its pollinating wasps in insular South East Asia." BMC Evolutionary Biology, 2017, 17, pp.207. 10.1186/s12862-017-1034-8. doi:10.1186/s12862-017-1034-8

HAL Id: hal-01589491

https://hal.sorbonne-universite.fr/hal-01589491

Submitted on 18 Sep 2017

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_3.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Diversification and spatial structuring in the mutualism between *Ficus septica* and its pollinating wasps in insular South East Asia

<figure><img src="image_4.png" /></figure>

CrossMark

Lillian Jennifer Rodriguez<fnref n="1" />,<fnref n="2" />,³*, Anthony Bain⁴, Lien-Siang Chou⁴, Lucie Conchou<fnref n="2" />,⁵, Astrid Cruaud³, Regielene Gonzales<fnref n="1" />, Martine Hossaert-McKey<fnref n="2" />, Jean-Yves Rasplus³, Hsy-Yu Tzeng⁶ and Finn Kjellberg<fnref n="2" />

**Abstract**

**Background:** Interspecific interactions have long been assumed to play an important role in diversification. Mutualistic interactions, such as nursery pollination mutualisms, have been proposed as good candidates for diversification through co-speciation because of their intricate nature. However, little is known about how speciation and diversification proceeds in emblematic nursery pollination systems such as figs and fig wasps. Here, we analyse diversification in connection with spatial structuring in the obligate mutualistic association between *Ficus septica* and its pollinating wasps throughout the Philippines and Taiwan.

**Results:** *Ceratosolen* wasps pollinating *F. septica* are structured into a set of three vicariant black coloured species, and a fourth yellow coloured species whose distribution overlaps with those of the black species. However, two black pollinator species were found to co-occur on Lanyu island. Microsatellite data on *F. septica* indicates the presence of three gene pools that broadly mirrors the distribution of the three black clades. Moreover, receptive fig odours, the specific message used by pollinating wasps to locate their host tree, varied among locations.

**Conclusions:** *F. septica* and its black pollinator clades exhibited similar geographic structuring. This could be due originally to geographic barriers leading to isolation, local adaptation, and finally co-structuring. Nevertheless, the co-occurrence of two black pollinator species on Lanyu island suggests that the parapatric distribution of the black clades is now maintained by the inability of migrating individuals of black pollinators to establish populations outside their range. On the other hand, the distribution of the yellow clade strongly suggests an initial case of character displacement followed by subsequent range extension: in our study system, phenotypic or microevolutionary plasticity has allowed the yellow clade to colonise hosts presenting distinct odours. Hence, while variation in receptive fig odours allows specificity in the interaction, this variation does not necessarily lead to coevolutionary plant-insect diversification. Globally, our results evidence evolutionary plasticity in the fig-fig wasp mutualism. This is the first documentation of the presence of two distinct processes in pollinating fig wasp diversification on a host species: the formation of vicariant species and the co-occurrence of other species over large parts of their ranges probably made possible by character displacement.

**Keywords:** Biogeography, *Ceratosolen*, Mutualism, Philippines, Speciation

<footnote>* Correspondence: ljvrodriguez@gmail.com</footnote>

<footnote n="1">¹Institute of Biology, University of the Philippines, Diliman, Quezon City, Philippines</footnote>

<footnote n="2">²CEFE UMR 5175, CNRS—Université de Montpellier—Université Paul-Valéry Montpellier—EPHE, Montpellier, France</footnote>

<footnote>Full list of author information is available at the end of the article</footnote>
and pollen but no seeds, while female trees bear figs that produce seeds but neither wasp offspring nor pollen [18]. Two polymorphic species (formally considered as subspecies) pollinate *F. septica*, namely *Ceratosolen bisulcatus bisulcatus* (Mayr) and *C. b. jucundus* Grandi [21]. Recent molecular analysis has demonstrated that these species are species complexes of closely related species [17], two of which were shown to co-occur in South Taiwan and in Luzon, Philippines, even at the fig-level [19, 22].

## General overview of sampling strategy

To get a broad image of the biogeographic distribution of *F. septica* pollinating wasp clades, we collected, sorted to morpho-species and sequenced pollinating wasps from all over the Philippines, from Taiwan and from Okinawa, and we examined published molecular data from the same region (Table 1). Genetic structuring within *F. septica* was investigated on a more limited scale, with four sampling points in the Philippines (one location in Mindanao, one in Negros, two in Luzon), one in Lanyu Island (60 km east of South Taiwan), and two in Taiwan (Table 1). The aim was to 1) confirm plant genetic homogeneity within islands and to detect potential genetic differentiation among islands and 2) to compare spatial genetic structuring within *F. septica* with the distribution of the different pollinator clades. Composition of the odours released by receptive figs (i.e. figs ready to be visited by pollinators), was investigated in four locations, one in Mindanao, one in Negros, one in Luzon and one in Taiwan, so as to represent the zones investigated for plant genetic structuring, and for which we also had information on pollinating wasps (Table 1). The aim was to detect whether differentiation in receptive fig

odours, a highly selected trait of the plants, correlated or not with genetic differentiation in neutral markers of the plants. A complementary aim was to detect whether a given wasp clade could be associated with plant populations exhibiting divergent receptive fig odours, i.e. to detect whether geographic variation in receptive fig odours constrained the distribution of the different wasp clades.

## Structuring into clades, distribution and phylogeny of *Ficus septica* pollinators

We sampled pollinator wasps from a large series of sites (Table 1, Fig. 1) from March 2011 to May 2014. Nearly ripe male figs were bagged and wasps were collected from the bags after their emergence from the figs. The samples were stored in vials with 70% ethanol. For each location, wasps were sorted to morpho-species and two to five wasp individuals of each morpho-species were sequenced for cytochrome *c* oxidase subunit I (*COI*), cytochrome *b* (*cyt b*) and elongation factor 1α (*EF1α*) (following [23, 24]). Extractions and amplifications were carried out at CBGP-INRA, Montferrier-sur-Lez, France.

## *Ficus septica* microsatellite analysis

To investigate genetic differentiation of *F. septica* among sites, we sampled leaves from seven study sites between April 2010 and September 2011. One leaf was collected from 17 to 37 trees per site in the Philippines and Taiwan (Table 1). Two pairs of collection sites were located on the same islands: on Luzon, they were about 100 km apart and on Taiwan, they were about 350 km apart. The leaves were oven-dried for two to three days and preserved dry. DNA was extracted from the leaf samples, amplified and

<table><caption>**Table 1** Collection sites for pollinator, leaf and fig odour samples</caption><thead><tr><th>Country</th><th>Locality</th><th><i>F. septica</i> pollinators</th><th><i>F. septica</i> leaves</th><th><i>F. septica</i> receptive fig odour</th></tr></thead><tbody><tr><td>Japan</td><td>Okinawa</td><td>[55]</td><td></td><td></td></tr><tr><td>Japan</td><td>Okinawa</td><td>This study</td><td></td><td></td></tr><tr><td>Taiwan</td><td>North</td><td>This study, [17, 23]</td><td>This study</td><td>This study</td></tr><tr><td>Taiwan</td><td>Throughout the country</td><td>[56]</td><td></td><td></td></tr><tr><td>Taiwan</td><td>South</td><td>This study, [17, 22]</td><td>This study</td><td></td></tr><tr><td>Taiwan</td><td>Lanyu Island</td><td>[17, 22, 56]</td><td>This study</td><td></td></tr><tr><td>Philippines</td><td>Northern Luzon</td><td>This study</td><td></td><td></td></tr><tr><td>Philippines</td><td>Central Luzon</td><td>This study</td><td>This study</td><td>This study</td></tr><tr><td>Philippines</td><td>Southern Luzon</td><td>This study, [57]</td><td>This study</td><td></td></tr><tr><td>Philippines</td><td>Panay Island</td><td>This study</td><td></td><td></td></tr><tr><td>Philippines</td><td>Negros Island</td><td>This study</td><td>This study</td><td>This study</td></tr><tr><td>Philippines</td><td>Palawan Island</td><td>This study</td><td></td><td></td></tr><tr><td>Philippines</td><td>Camiguin Island</td><td>This study</td><td></td><td></td></tr><tr><td>Philippines</td><td>Mindanao Island</td><td>This study</td><td>This study</td><td>This study</td></tr></tbody></table>

Table 1 Legend: Sites are arranged in a north-to-south orientation. Also indicated are the references of the GenBank sequences (Additional file 8) that were blasted on our sequences for clade assignation. Pollinator genes analysed - [17, 50–52]: COI; [48]: COI, 28S; [49]: 18S, 28S, COI, cyt b, Wg; this study: COI, cyt b, EF1α. The COI sequences allow the assignment of every sequenced individual to clade and subclade recognised in this study
<figure><img src="image_5.png" /><figcaption>**Fig. 1** Geographic distribution of the *Ficus septica* gene pools and the different clades (and subclades) of the pollinators of *F. septica*</figcaption></figure>

genotyped (following [25]). Fourteen nuclear microsatellite markers specifically designed for the genus *Ficus* were used [25]. They were selected because of their high amplification rates and polymorphism (Additional file 1). Genotyping was carried out using an ABI 3730 Sequencer (Applied Biosystems).

**Ficus septica** receptive fig odour collection and analysis

To collect odours from receptive figs, 5–20 *F. septica* trees were sampled per study site from March to September 2011. Scent extraction was carried out from 10.00 to 14.00 when scent emission is maximum [19]. As much as possible, male and female trees were equally represented in
the sampling. Moreover, in the Philippines, fig scent collection was performed twice per site, at different seasons. Odours were collected using the headspace technique. Thirty to seventy receptive figs from each of the selected trees were placed inside a polyethylene terephthalate bag. The bag was closed for 30 min to allow the scent to accumulate. After the accumulation period, air was pumped out of the bag through two tubes, each with a chromatoprobe filter, for 5 min at a rate of 158 mL min⁻¹. The volatile compounds were adsorbed into these filters, which were previously injected with known amounts of standard compounds (nonane and dodecane) to allow quantification of the volatiles. Empty bags were used as controls. The filters were placed in vials and stored at -18 °C until further analysis.

To identify and quantify the compounds present in the scent samples, filters were injected in a gas chromatograph (CP-3800, Varian Inc., Palo Alto, CA) coupled with a mass spectrometer (Saturn 2000, Varian; for identification of compounds) and equipped with a flame ion detector (for quantification of compounds). Samples were injected using a 1079 programmed temperature injector with a chromatoprobe thermodesorption kit (Varian). Temperature program and other conditions followed [19]. Component identification was based on the gas chromatograph (GC) retention indices of the compounds compared to those reported in literature [26] and through computer matching of mass spectra (MS) with NIST'98 MS library. These analyses were performed at the *Plate-Forme d'Analyses Chimique en Ecologie* of the LabEx CeMEB.

## Data analyses

For the phylogenetic analysis of pollinator wasp sequences, gene sequences were cleaned and edited in Geneious v5.5.9 and checked for presence of pseudogenes in Mega v6.06. We ran maximum likelihood (ML) analyses on the mitochondrial DNA dataset (COI - 1527 bp + cyt b - 734 bp) and on the nuclear DNA dataset (EF1α - 517 bp) using a GTR + Γ4 model with 1000 bootstrap replicates. We had 54 COI, 65 cyt b and 58 EF1α in-group sequences. We also used sequences of pollinator wasps of *F. botryocarpa*, *F. malayana* and *F. uncinata* as outgroups. In total, we had 13 COI, 10 cyt b and 8 EF1α outgroup sequences.

To quantify genetic variation between the *F. septica* populations, heterozygosity, mean number of alleles, and FST values were computed using *GENALEX* v6.5 [27] and *SPAGEDI* v1.4 [28]. Exact tests of inter-population genetic differentiation were carried out using Fisher's method in *GENEPOP* v4.2 [29]. Finally, the microsatellite data was subjected to Bayesian clustering analysis using *STRUCTURE* v2.3.4 [30] in order to determine the different genetic groups. The parameter set was a burn-in period of 100,000 steps and a run of 1,000,000 steps with no admixture, repeated five times. The number of individual clusters was determined using the Evanno method implemented in *STRUCTURE* HARVESTER [31].

We used the program *STRUCTURE* PLOT to visualize the Bayesian clusters [32].

Significant differences in total quantity of volatiles were tested using One-way ANOVA. For the subsequent analyses, relative abundance data was used. Fig scent composition was analyzed by constructing several models of MANOVA using geography, sex and season (or combinations of these) as explanatory variables. Significance values of each response variable were evaluated using the false discovery rate [33]. Spread of the variation in the data was visualized in a principal components analysis (PCA) performed using the package *FACTOMINER* v1.31.5. A dissimilarity matrix using the Bray-Curtis index was also generated. Inter-population comparison of the blends was done by subjecting the Bray-Curtis indices to permutational analysis of variance (permanova) tests (permutations = 100,000) in the package *VEGAN* v2.3-3. All analyses were performed in R v3.2.1.

## Results

### *Ficus septica* microsatellite analysis

Genetic parameters averaged over 14 *Ficus septica* microsatellite loci used in this study are shown in Additional files 1 and 2. Expected heterozygosity ($H_e$) and number of alleles per locus (Na) exhibited minimal variation in the seven sites (0.313–0.405) while private alleles (PA) varied more between islands. Exact tests of genotypic differentiation showed that all the population pairs, except the two from Luzon Island, were differentiated.

$F_{ST}$ values showed that the least differentiated populations were those from Luzon and those from Taiwan and Lanyu (Additional file 3). An intermediate homogeneity between the Luzon populations and the Negros population was evidenced by intermediate $F_{ST}$ values demonstrating some genetic differentiation. The Bayesian clustering gave an optimal number of clusters, K, equal to 3 (Fig. 2). From South to North, the clusters gathered 1) the lone Mindanao population, 2) the Negros and Luzon populations, 3) the Lanyu Island population and the two Taiwan populations (Fig. 2).

### *Ficus septica* receptive odours

A total of 39 compounds were found in the scent of *F. septica*, with a little less than half classified as sesquiterpenes (Additional file 4). However, the two dominant compounds, (E)-β-ocimene and (Z)-3-hexen-1-ol, were not sesquiterpenes but an acyclic monoterpene and a simple aliphatic, respectively. The most abundant fig odour component was α-pinene in Mindanao samples, (Z)-3-hexen-1-ol in Negros samples, (E)-β-ocimene in Central Luzon samples, and linalool in North Taiwan samples. Total amount of volatiles emitted did not differ significantly between locations, seasons and sexes (One-way ANOVA, $F_{10,45} = 1.853$, *P* = 0.078). Of the 39 compounds, 18 had a mean amount of at least 1% of the blend (Additional file 4). To allow
<figure><img src="image_6.png" /><figcaption>**Fig. 2** Plant genetics. Bar plots of the Bayesian clustering runs at the best K (3). Each column displays the probability of each individual to belong to each of the clusters. The three groups corresponded to i) Taiwan including Lanyu Island, ii) Luzon plus Negros Island, iii) Mindanao Island</figcaption></figure>

significance testing, only the abundances of these 18 compounds were used for subsequent analyses.

Using relative abundance values of the 18 compounds, we found that geography/location, season and their interaction factor were significant in explaining the variation in odour blend of receptive figs from Philippines (permanova, geography: *F* = 6.93, *P* < 0.0001; season: *F* = 2.55, *P* = 0.016; geography*season: *F* = 1.896, *P* = 0.026). Sex did not have a significant effect on scent variation (permanova, sex: *F* = 0.49, *P* = 0.868). When the Philippine dataset was partitioned into the three different sites, season had a significant effect in the scent variation in both Negros and Central Luzon, but not in Mindanao (permanova: Negros – *F* = 2.60, *P* = 0.012; Central Luzon – *F* = 2.89, *P* = 0.021; Mindanao – *F* = 1.20, *P* = 0.30; Additional file 5). In an analysis of all samples, including Taiwan, geography/location still showed a significant effect and sex still presented no effect on the blend (MANOVA, geography: approx. *F* = 5.45, *P* < 0.0001; sex: approx. *F* = 0.44, *P* = 0.967). Through pairwise comparisons of samples from different locations, it was found that the odour blend of the four *F. septica* populations significantly differed from each other (Additional file 6). This result was supported by the findings that the main compounds in the blends of the four sites differed (Additional file 7). PCA also reflected this differentiation (Fig. 3). The first two PCs explain about 24% of the variation in the data. The first axis separated the Mindanao samples from the other samples. The second axis separated the remaining three groups. In this representation, the barycenters of fig odours from Central Luzon and Negros were located next to each other.

In summary, data on receptive fig odour show that Mindanao, Negros, Central Luzon, and North Taiwan each have their own distinct odour profiles. These also show that

geography/location is the main factor affecting fig scent variation, with the sex of the tree contributing very little, if at all, to the variation seen. Lastly, in the Philippines, season also accounted to some extent for scent variation in Negros and Central Luzon but not in Mindanao.

## Structuring into clades and phylogeny of *Ficus septica* pollinators

Number of individuals sequenced, location of collection and GenBank accession numbers are given in Additional file 8. The phylogenetic trees derived from nuDNA (Fig. 4A) and mtDNA (Fig. 4B) show the presence of four major clades associated with *F. septica*. In our samples, three major pollinator clades (1–3) corresponded to black coloured morphs and the last major clade (clade 4) was composed of more or less yellow-coloured morphs, either dark backed-yellow sided (subclades 4A and 4B) or totally yellow (subclades 4C and 4D). Clades 1 to 3 followed a south to north distribution: clade 1 was present only in Mindanao, clade 2 ranged from Negros Island to Luzon (plus, according to COI sequences, one specimen collected on Lanyu Island and analyzed in [22]), and clade 3 occurred from Lanyu Island to Okinawa Island (Fig. 1). On the other hand, clade 4 showed less geographic structuring and was present throughout the Philippines up to South Taiwan (Fig. 1). Subclade 4A occurred in Mindanao and neighboring Camiguin Island, subclade 4B stretched from Palawan to the southern tip of Luzon. Subclade 4C ranged from Central to North Luzon and subclade 4D was present in Lanyu Island and South Taiwan. Clade 1 was found in sympatry with subclade 4A. Clade 2 coexisted with subclade 4C (Luzon) and with subclade 4D and clade 3 (Lanyu Island). Lastly, subclade 4B was not found in sympatry with any other clade. All previously published sequences from our sampling region fitted
<figure><img src="image_7.png" /><figcaption>**Fig. 3** Plant odours. Principal Components Analysis of relative abundance of each compound in receptive fig odour showing groupings according to location. (North Taiwan – Yellow, Central Luzon – Red, Negros Island – Blue, Mindanao Island – Green). Odour profiles vary significantly between each of the four study sites (MANOVA, geography: approx. F = 5.45, P < 0.0001)</figcaption></figure>

into the picture presented here. Samples from Indonesia and New Guinea [17] belonged to different clades from those recovered in this study.

## Discussion

### Genetic and odour structuring in *Ficus septica*

*F. septica* was spatially structured into three separate gene pools (from South to North): 1) the Mindanao population, 2) the Negros and Luzon populations, 3) the Lanyu Island population and the two Taiwan populations (Fig. 2). However, this genetic structuring was not reflected by the receptive fig odours because Mindanao, Negros, Central Luzon, and North Taiwan each had their own distinct odour profiles (Fig. 3). Hence receptive fig odours seem to evolve differently from neutral markers, suggesting selection. Beyond the case of *F. septica*, geographic variation in receptive fig odours has also been demonstrated in two fig species, *F. hispida* and *F. racemosa*. At least in *F. racemosa*, receptive fig odours varied between plants belonging to different gene pools and pollinated by different wasp species, while receptive fig odours did not differ between plants growing 800 km apart but belonging to the same gene pool and pollinated by the same wasp species [14, 34]. This similarity

a

<figure><img src="image_8.png" /><figcaption>**Fig. 4** Insect phylogeny. Maximum likelihood (ML) trees from (a) nuclear DNA (EF1α gene sequences) and (b) mitochondrial DNA (COI and cyt b gene sequences) showing the four pollinator clades (clades 1 to 4) locally-associated with *Ficus septica*. Clade 4 in (b) is further subdivided into subclades 4A to 4D</figcaption></figure>

b

<figure><img src="image_9.png" /></figure>
in fig odours suggests stabilising selection at work. Further, in the exceptional situation where two *Ficus* species freely share pollinators, their receptive fig odours are not distinguishable, suggesting selection by the wasps for evolutionary convergence of receptive fig odours [35, 36]. While we do not provide a direct demonstration that the difference in receptive fig odour reported here is meaningful for pollinating wasp attraction, more rapid evolution than neutral markers among islands is strikingly different from the stability and convergence of receptive fig odours documented in other *Ficus* species. Therefore we hypothesise that the differentiation of receptive fig odours in our study system is meaningful in terms of pollinator attraction.

The importance of receptive fig odours in pollinating wasp attraction has been largely demonstrated [37–39]. In contrast, variation within fig species [34] and among closely related species has been much less investigated [35, 36]. We document here seasonal variation in receptive fig odours in Luzon and Negros, but not in Mindanao, a difference that may be explained by limited seasonality in Mindanao comparatively to the two other islands [40]. This fits the very limited available data showing seasonal variation in receptive fig odours, measured within or between sexes depending on the species’ flowering phenology [20, 41]. Nevertheless, in all documented cases, differences between seasons remain limited, more limited than differences among species (except for *Ficus* species sharing pollinators [35, 36]) or than differences among islands for *F. septica*.

## Spatial structuring in *Ficus septica* and its pollinator clades

In *F. septica*, we show that gene pools do not neatly costructure with receptive fig odours, as one gene pool (Luzon-Negros) corresponded to two distinct scent profiles (Luzon and Negros). This difference provides us with a new perspective to analyse pollinator structuring. Will pollinator species/populations reflect the geographic structuring of *F. septica* gene pools? Or will pollinators co-structure with receptive odour profiles, which are the main wasp attracting agents?

The pollinators of *Ficus septica* formed a monophyletic entity, structured into four major genetic clades that can be recognized morphologically, and four subclades within clade 4 that could not be distinguished morphologically, to the exception of colouration differences. According to available data on genetic distances within and among fig pollinating wasp species, the four major clades should be considered as distinct species while the subclades should be conservatively interpreted as within species variation [42]. The type of *Ceratosolen jucundus* was sampled on Mount Makiling [43], where we collected subclade 4C. Within this framework, *Ceratosolen bisulcatus jucundus* corresponds to clade 4 and is upgraded to species rank *Ceratosolen jucundus* (stat. nov.), while clades 1 to 3 correspond to three

undescribed species that all appear different from *C. bisulcatus* (stat. nov.).

The three black coloured clades were allopatric to parapatric, presenting a geographic South-North zonation (Fig. 1). The black clades presented a zonation similar to that of the *F. septica* gene pools (Fig. 2). While sampling was not strictly comparable preventing formal statistical comparison, the similarity in spatial structure of clade distributions was striking. We hypothesize that this co-structuring is due to the presence of geographic barriers between the populations leading to isolation and subsequent diversification. *F. septica* and the black pollinator clades might have colonized the islands along the same South-to-North direction. Upon establishment in the different islands, populations must have been geographically isolated, subsequently facilitating local adaptation and diversification, such that the co-structuring between the fig and the wasps is evident now. Co-structuring, in this case, is thus driven more by geographic (and accompanying climatic) differences rather than by strict plant-insect reciprocal adaptation. Indeed, according to the Köppen-Geiger climate classification, the climates of Luzon and Negros are similar while the climates of Mindanao and Taiwan are more different [40]. This climatic differentiation is in broad concordance with the distributions of pollinator clades 1 to 3 and *F. septica* gene pools.

How do these patterns relate to receptive odour structuring? If an olfaction barrier maintains the specificity of the fig-wasp interactions, then we would expect each pollinator clade to recognize a single receptive odour profile. Nevertheless, geographic variation in receptive fig odours in *F. septica* has not precluded the presence of the following three situations: (1) coexistence of clades 2 and 3 on Lanyu Island, (2) presence of clade 2 in Luzon and Negros despite differing fig scent, and if our scenario is correct, (3) geographic expansion of clade 4 throughout the Philippines and up to Taiwan, bridging the different receptive fig odours produced by *F. septica*. These suggest either that wasps can rapidly evolve detection of new odours or that the wasps can still recognize receptive figs because the odour differences are limited. Hence, while receptive fig odour is a cue wasps use to locate their host figs [39], its differentiation among fig populations or among related *Ficus* species may allow phenotypic or micro-evolutionary plasticity in wasp odour recognition rather than being a constraint that would lead to strict insect-plant coevolutionary diversification. We may conjecture, for instance, that response to receptive fig odours, and maybe even perception of odours in wasps of clade 2 may vary between wasps originating from Negros, from Luzon and from Lanyu Island. Future experimental studies are needed to establish possible differential attraction to the various receptive odour compounds.

We therefore demonstrate here that structuring in *F. septica* and their pollinator clades is driven mainly by
geographic (and climatic) factors rather than a species-specific olfaction barrier.

## Pollinator clade distribution and diversification

It should be noted that clades 2 and 3 are both present in Lanyu Island, which could be considered as their contact zone. Although these two clades are distributed over large stretches of sea presenting only few small islands, they still retain within-clade genetic homogeneity. This suggests sufficient effective numbers of migrants and hence sufficient dispersal. How then is the co-structuring maintained? What prevents expansion of clade 2 into South Taiwan (a distance of about 65 km from Lanyu Island) when it is able to hop along a series of islands up to 100 km apart from Luzon to Lanyu? Our working hypothesis is that migrant pollinators are unable to establish viable populations outside their realized range. This might be because local pollinators are more adapted to local conditions than the migrants, leading to competitive exclusion of migrating species. This pattern was also exhibited by the non-pollinating fig wasps associated with *Ficus rubiginosa* in Australia, where parapatric species pairs had smaller ranges than expected probably due to competition [44]. Therefore, the parapatric distribution of the black clades might be maintained by the inability of migrating species to establish populations outside their range.

Within this context, the distribution and incipient diversification of clade 4 is of particular interest. Clade 4 is present throughout the Philippines and extends 60 km northwards into South Taiwan, bridging almost the whole set of clients. However, the distribution of subclades is structured geographically with a South-to-North zonation. Also, subclades 4C and 4D are characterized by completely yellow colour in all the specimens we collected and their distributions totally overlap with the distributions of black coloured clades. Yellow colour is exceptional in diurnal fig pollinator species [19]. Individuals of subclade 4C have been demonstrated to exhibit very short survival times outside figs comparatively to co-occurring individuals of the black coloured clade 2 [19]. The authors interpreted the yellow colour of the clade as part of an adaptive syndrome based on short survival but high competitiveness, a “live fast, die young” strategy. Hence, clade 4 seems to have evolved a modified life strategy, which has enabled colonization of areas where *F. septica* is also pollinated by a black coloured clade and the subsequent stable coexistence of the two wasp species. This suggests recent range expansion into regions inhabited by clades 1 to 3 and ongoing differentiation into distinct gene pools. Such a scenario would involve initial ecological character displacement in a region of co-occurrence followed by subsequent colonization of the regions inhabited by the other black clades. In this sense, the initial clade 4 colonization of habitats already occupied by a black clade

would qualify as ecological character displacement sensu stricto, while subsequent range extension would be made possible by differences evolved in allopatry, in a process of species sorting [5].

Interestingly, our results show the widespread co-occurrence of yellow *C. jucundus* wasps with the black pollinator clades. This has been previously reported, but on a much smaller scale [19, 22, 45]. Co-existence of reproductively isolated yellow and black pollinator clades were also documented in *F. rubiginosa* collected from Townsville, Australia [42, 46]. Aside from these studies and a few more [47], co-occurrence of sympatric, actively pollinating species over large areas has been rarely documented and their population genetics barely characterized. Our study fills this gap and we have proposed mechanisms on how this co-existence was achieved.

## Biogeography and adaptation in relation to other fig-wasp systems

How do our results compare with what is known about biogeography and adaptation in the fig-pollinating wasp system? The situation in *F. racemosa* resembles the situation of *F. septica* and its black pollinator clades, but on a broader scale. Indeed, the pollinators are structured into a set of wide-ranging vicariant species and the plant is structured along the same geographic lines into different gene pools, the limit of which seem to correspond to geographic barriers limiting gene flow [14]. Similarly, results on *F. hirta* and its pollinator in South East China show overall strong gene flow over continental expanses in both plant and insect, with genetically divergent populations for both plant and insect on Hainan island [48]. The authors suggested that the geographic barrier in itself was not sufficient to explain the data, and that isolation by adaptation, facilitated by a decrease in dispersal across a geographic barrier, had to be invoked [48]. In the data presented here, the presence of clade 2 in Lanyu Island supports this vision, as it shows that fig pollinating wasp species can achieve continuity across almost 400 km of sea, hopping from island to island with some over 100 km apart: the 30 km separating Hainan island from the continent can hardly be considered by themselves a major obstacle for fig wasp dispersal. Moreover, data on *F. pumila*, a creeping fig, and its pollinators in Eastern China showed some genetic isolation by distance in the plant [49], and separation of the pollinators into two wide ranging vicariant species. South of Shanghai, in the Zhoushan archipelago, within a context of some geographic obstacles to dispersal, a localized third pollinator species is present in combination with one of the wide-ranging species [50]. The two co-occurring species differ in precocity of emergence from figs after over-wintering, suggesting differences in temperature thresholds. This ecological differentiation may facilitate the coexistence of competing species [50], while the parapatric distribution of the two
wide-ranging species is probably best explained by competitive exclusion in a context of strong gene flow limiting local adaptation in widespread species [51]. Finally, *F. rubiginosa* provides the most complex pattern of distribution of diverse pollinator species documented to date, a feature that may result from the range of the host tree encompassing a combination of two climatic gradients, a progressive North-South gradient and a steep East-West one [42, 46, 52]. Hence, the biogeographic patterns observed in fig-pollinating wasp associations tend to support a global vision in which the pollinators of wide ranging fig species, but not necessarily the fig species, are often structured into broad, geographically vicariant, genetically homogeneous entities. New lineages of pollinating wasp species may colonize such systems following a local event of ecological character displacement. This can be followed by subsequent colonization of parts of the range of the fig occupied by other pollinators through a process of species sorting resulting in the regional coexistence of two pollinating wasps on the same host. We were able to demonstrate these general emergent patterns for pollinating fig wasp species diversification in a relatively simple fig-wasp system. This leads us to suggest that within *Ficus* species, in an insular setting, fig-fig wasp diversification and co-structuring may be frequent. These patterns do not need to result from coevolutionary diversification [1]. Indeed, in our fig-fig wasp system, geographic barriers seem to have been the main driving factor leading to plant-insect co-structuring, a feature that does not necessarily entail parallel cladogenesis [14]. Interestingly, the co-occurrence of clades 2 and 3 in Lanyu Island suggests that barriers that might arise from coevolutionary diversification can be overcome, limiting its role in the process of joint plant-insect diversification.

The Philippines constitute a model island archipelago to investigate evolutionary processes of diversification [53]. Much focus has been on limited dispersal species, for which routes of colonization can be reconstructed using phylogenies. The results presented here on *Ficus septica* allow a different perspective. What biogeographic structure do we get in strong dispersal model species and what does it tell us about evolutionary processes at work in lower dispersal species? Our study may shed light on results such as those on shrews [54] as it supports a role for geographic barriers, but probably also for ecological character displacement followed by species sorting in diversification processes. Our results also emphasize the importance of climatic variation among islands as a force structuring diversity and diversification.

## Conclusions

We have shown that *F. septica* and its black pollinator clades exhibited similar geographic structuring. We hypothesize that this could be due originally to strong geographic (and climatic) barriers leading to isolation, local

adaptation, and finally plant-insect co-structuring. We also show the co-existence of pollinating species in an area, which probably had been brought about by initial character displacement and subsequent range expansion. This is the first documentation of the presence of two distinct processes in pollinating fig wasp diversification on a host species: the formation of vicariant species and the co-occurrence of other species over large parts of their ranges probably made possible by character displacement. Lastly, we provide evidence that receptive fig odours do not constrain the specificity of the fig-pollinator wasp interactions.

## Additional files

**Additional file 1:** Characteristics of 14 microsatellite loci used in the study across all *Ficus septica* individuals. $H_e$ - observed heterozygosity, $H_e$ - expected heterozygosity, Na - number of alleles per locus, NA - Not Applicable. (DOCX 61 kb)

**Additional file 2:** Genetic parameters from 14 microsatellite loci of *Ficus septica*. n - number of leaves sampled, $H_e$ - expected heterozygosity, Na - mean number of alleles per locus, PA - mean number of private alleles. (DOCX 51 kb)

**Additional file 3:** Genetic differentiation in *Ficus septica* populations as shown from $F_{ST}$ values computed from 14 microsatellite loci. All $F_{ST}$ values (besides between the two populations from Luzon) are significant at *P* < 0.05 using the individual permutation tests. TW - Taiwan; PH - Philippines. (DOCX 62 kb)

**Additional file 4:** Volatile organic compounds emitted by *Ficus septica* from the Philippines and North Taiwan. Compounds with more than 1% mean relative abundance are shown in boldface. (DOCX 79 kb)

**Additional file 5:** Plot of the mean relative abundance of the 18 most abundant compounds found in *Ficus septica* odours between two seasons (dry and rainy) in two Philippine sites: Central Luzon and Negros Island. (PDF 9 kb)

**Additional file 6:** Permutational analysis of variance (permanova) comparisons of the fig odour profiles (relative abundance) of *Ficus septica* populations from North Taiwan and three Philippine sites (Central Luzon, Negros Island, Mindanao Island) (*-significant difference at *P* < 0.05). Significant *P*-values for all comparisons show that each of the 4 sites has its own distinct odour profile. (DOCX 56 kb)

**Additional file 7:** Plot of the mean relative abundance of the 18 most abundant compounds found in *Ficus septica* odours from North Taiwan and three sites in the Philippines (Central Luzon, Negros Island, Mindanao Island): the most abundant compound is different in each of the sites. (PDF 9 kb)

**Additional file 8:** Sampling table for *Ceratosolen bisulcatus* wasps indicating collection locality and GenBank accession number for each sequence used for the phylogenetic analysis. (DOCX 108 kb)

## Abbreviations

ANOVA: analysis of variance; COI: cytochrome c oxidase subunit I; cyt b: cytochrome b; EF1a: elongation factor 1α; GC: gas chromatograph; GTR: general time reversible; $H_e$: expected heterozygosity; MANOVA: multivariate analysis of variance; ML: maximum likelihood; MS: mass spectra; mtDNA: mitochondrial DNA; Na: number of alleles per locus; nuDNA: nuclear DNA; PA: private alleles; PC: principal component; PCA: principal components analysis; Permanova: permutational analysis of variance

## Acknowledgements

We are grateful to Jean-Marie Bessiere, Marjorie Garcia, Bruno Buatois, Florian Massip, Rian Bernal, Nonoy Dulman, Luciano Palmieri, Ian Fontanilla, Dino Ramos, Yoyi Rodriguez and Reynand Canoy for helping us in the identification of compounds, laboratory work and field collection.

## Funding

This work was supported by the University of the Philippines Office of the Vice-Chancellor for Research and Development, Commission on Higher
Education (Philippines), ANR-09-BLAN-0392-CSD 7, NSC 99-2923-B-002-001-MY3, Bibliothèque du Vivant funded by the CNRS, Muséum National d'Histoire Naturelle and the Institut National de la Recherche Agronomique.

**Availability of data and materials**

Microsatellite genotyping data, odour profiles and gene alignments are available from the corresponding author on request. Wasp DNA sequences are deposited in GenBank (accession numbers are listed in the Additional file 8).

**Authors' contributions**

MHM, LSC, JYR and FK conceived the study. LJR, AB, LC, AC, RG, HYT, JYR and FK performed the experiments and data collection. LJR, AB, LC, AC, JYR and FK analysed and interpreted the data. LJR and FK wrote the manuscript. All authors commented on the manuscript and gave final approval for publication.

**Ethics approval and consent to participate**

We confirm that the leaf and insect collection presented here were conducted in accordance with the Wildlife Resources Conservation and Protection Act (Republic Act 9147) set forth by the Philippine government, the Department of Environment and Natural Resources and the University of the Philippines.

**Consent for publication**

Not applicable

**Competing interests**

The authors declare that they have no competing interests.

**Publisher's Note**

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Author details**

¹Institute of Biology, University of the Philippines, Diliman, Quezon City, Philippines. ²CEFE UMR 5175, CNRS—Université de Montpellier—Université Paul-Valéry Montpellier—EPHE, Montpellier, France. ³INRA, UMR 1062 CBGP, Montferrier-sur-Lez, France. ⁴Institute of Ecology and Evolutionary Biology, National Taiwan University, Taipei, Taiwan. ⁵Institut d'Ecologie et des Sciences de l'Environnement de Paris—ECOSENS, INRA-UPMC, Versailles, France. ⁶Department of Forestry and Natural Resources, National Chung-Hsing University, Taichung, Taiwan.

Received: 17 March 2017 Accepted: 4 August 2017

Published online: 29 August 2017

**References**

1. Hembry DH, Yoder JB, Goodman KR. Coevolution and the diversification of life. Am Nat. 2014;184:425–38.

2. Stuart YE, Campbell TS, Hohenlohe PA, Reynolds RG, Revell LJ, Losos JB. Rapid evolution of a native species following invasion by a congener. Science. 2014;346:463–6.

3. Grant PR, Grant BR. How and why species multiply: the radiation of Darwin's finches. New Jersey: Princeton University Press; 2008.

4. Schluter D, McPhail JD. Ecological character displacement and speciation in sticklebacks. Am Nat. 1992;140:85–108.

5. Stuart YE, Losos JB. Ecological character displacement: glass half full or half empty? Trends Ecol Evol. 2013;28:402–8.

6. Thompson JN. The coevolutionary process. Chicago: The University of Chicago Press; 1994.

7. Ramirez SR, Eltz T, Fujiwara MK, Gerlach G, Goldman-Huertas B, Tsutsui ND, et al. Asynchronous diversification in a specialized plant-pollinator mutualism. Science. 2011;333:1742–6.

8. Yokoyama J. Cospeciation of figs and fig-wasps: a case study of endemic species pairs in the Ogasawara Islands. Popul Ecol. 2003;45:249–56.

9. Kawakita A. Evolution of obligate pollination mutualism in the tribe Phyllanthetheae (Phyllanthaceae). Plant Spec Biol. 2010;25:3–19.

10. Cruaud A, Ránsted N, Chantarasuwan B, Chou LS, Clement WL, Couloux A, et al. An extreme case of plant-insect codiversification: figs and fig-pollinating wasps. Syst Biol. 2012;61:1029–47.

11. Cook JM, Segar ST. Speciation in fig wasps. Ecol Entomol. 2010;35:54–66.

12. Kiester AR, Lande R, Schemske DW. Models of coevolution and speciation in plants and their pollinators. Am Nat. 1984;124:220–43.

13. Yoder JB, Smith CI, Rowley DJ, Flatz R, Godsoe W, Drummond C, et al. Effects of gene flow on phenotype matching between two varieties of Joshua tree (Yucca brevifolia; Agavaceae) and their pollinators. J Evolution Biol. 2013;26:1220–33.

14. Bain A, Borges RM, Chevallier MH, Vignes H, Kobmoo N, Peng YQ, et al. Geographic structuring into vicariant species-pairs in a wide-ranging, high-dispersal plant-insect mutualism: the case of *Ficus racemosa* and its pollinating wasps. Evol Ecol. 2016;30:663–84. doi:10.1007/s10682-016-9836-5.

15. Hembry DH, Okamoto T, Gillespie RG. Repeated colonization of remote islands by specialized mutualists. Biol Lett. 2011;8:258–61.

16. Cruaud A, Rasplus JY, Jabbour-Zahab R, Genson G, Ungricht S. Testing the emergence of New Caledonia: fig wasp mutualism as a case study and a review of evidence. PLoS One. 2012;7:e30941.

17. Moe AM, Weiblen GD. Molecular divergence in allopatric *Ceratosolen* (Agaonidae) pollinators of geographically widespread *Ficus* (Moraceae) species. Ann Entomol Soc Am. 2010;103:1025–37.

18. Berg CC, Corner EJH. Moraceae (*Ficus*). In: HPNooteboom, editor. Flora Malesiana. Netherlands: Nationaal Herbarium Nederland, Universiteit Leiden branch. 2005. p. 1–730.

19. Conchou L, Cabioch L, Rodriguez LJV, Kjellberg F. Daily rhythm of mutualistic pollinator activity and scent emission in *Ficus septica*: ecological differentiation between co-occurring pollinators and potential consequences for chemical communication and facilitation of host speciation. PLoS One. 2014;9:e103581.

20. Hossaert-McKey M, Proffit M, Soler C, Chen C, Bessière JM, Schatz B, et al. How to be a dioecious fig: chemical mimicry between sexes matters only when both sexes flower synchronously. Sci Rep. 2016;6:21236.

21. Wiebes JT. Taxonomy and host plant preference of indo-Australian fig wasps of the genus *Ceratosolen* (Agaonidae). Tijdschr Entomol. 1963;106:1–112.

22. Lin RC, Yeung CKL, Fong JJ, Tzeng HY, Li SH. The lack of pollinator specificity in a dioecious fig tree: sympatric fig-pollinating wasps of *Ficus septica* in southern Taiwan. Biotropica. 2011;43:200–7.

23. Cruaud A, Jabbour-Zahab R, Genson G, Cruaud C, Couloux A, Kjellberg F, et al. Laying the foundations for a new classification of Agaonidae (hymenoptera: Chalcidoidea), a multilocus phylogenetic approach. Cladistics. 2010;26:359–87.

24. Cruaud A, Jabbour-Zahab R, Genson G, Couloux A, Peng YQ, Yang DR, et al. Out of Australia and back again: the world-wide historical biogeography of non-pollinating fig wasps (hymenoptera: Sycophaginae). J Biogeogr. 2011;38:209–25.

25. Garcia M, Bain A, Tzeng HY, Peng YQ, Chou LS, Kjellberg F. Portable microsatellite primers for *Ficus* (Moraceae). Am J Bot. 2012;99:e187–92.

26. Adams RP. Identification of essential oil components by gas chromatography/mass spectroscopy. 4th ed. USA: Allured Publishing Corporation; 2007.

27. Peakall R, Smouse PE. GenAlEx 6.5: genetic analysis in excel. Population genetic software for teaching and research—an update. Bioinformatics. 2012;28:2537–9.

28. Hardy OJ, Vekemans X. SPAGEDI: a versatile computer program to analyse spatial genetic structure at the individual or population levels. Mol Ecol Notes. 2002;2:618–20.

29. Rousset F. GENEPOP'007: a complete re-implementation of the GENEPOP software for windows and Linux. Mol Ecol Resour. 2008;8:103–6.

30. Pritchard JK, Stephens M, Donnelly P. Inference of population structure using multilocus genotype data. Genetics. 2000;155:945–59.

31. Earl DA. vonHoldt BM. STRUCTURE HARVESTER: a website and program for visualizing STRUCTURE output and implementing the Evanno method. Conserv Genet Res. 2011;4:359–61.

32. Ramasamy RK, Ramasamy S, Bindroo BB, Naik VG. STRUCTURE PLOT: a program for drawing elegant STRUCTURE bar plots in user friendly interface. SpringerPlus. 2014;3:431.

33. Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J Roy Stat Soc B Met. 1995;57:289–300.

34. Soler C, Hossaert-McKey M, Buatois B, Bessière JM, Schatz B, Proffit M. Geographic variation of floral scent in a highly specialized pollination mutualism. Phytochemistry. 2011;72:74–81.

35. Cornille A, Underhill JG, Cruaud A, Hossaert-McKey M, Johnson SD, Tolley KA, et al. Floral volatiles, pollinator sharing and diversification in the figwasp mutualism: insights from *Ficus natalensis*, and its two wasp pollinators (South Africa). Proc R Soc B. 2012;279:1731–9.

36. Wang G, Cannon CH, Chen J. Pollinator sharing and gene flow among closely related sympatric dioecious fig taxa. Proc R Soc B. 2016;283:20152963.
37. Hossaert-McKey M, Gibernau M, Frey JE. Chemosensory attraction of fig wasps to substances produced by receptive figs. Entomol Exp Appl. 1994;70:185-91.

38. Grison-Pigé L, Bessière JM, Hossaert-McKey M. Specific attraction of fig-pollinating wasps: role of volatile compounds released by tropical figs. J Chem Ecol. 2002;28:267-79.

39. Chen C, Song Q, Proffit M, Bessière JM, Li Z, Hossaert-McKey M. Private channel: a single unusual compound assures specific pollinator attraction in *Ficus semicordata*. *Funct Ecol.* 2009;23:941-50.

40. Peel MC, Finlayson BL, McMahon TA. Updated world map of the Köppen-Geiger climate classification. Hydrol Earth Syst Sci. 2007;11:1633-44.

41. Soler CC, Proffit M, Bessiere JM, Hossaert-McKey M, Schatz B. Evidence for intersexual chemical mimicry in a dioecious plant. Ecol Lett. 2012;15:978-85.

42. Darwell CT. Al-Beid S, cook JM. Molecular species delimitation of a symbiotic fig-pollinating wasp species complex reveals extreme deviation from reciprocal partner specificity. BMC Evol Biol. 2014;14:189.

43. Grandi G. Hyménoptères sycophiles récoltés aux îles Philippines par C.F. Baker, I. Agaonini. Philipp J Sci. 1927;33:309-29.

44. Darwell CT, Cook JM. Cryptic diversity in a fig wasp community – morphologically differentiated species are sympatric but cryptic species are allopatric. Mol Ecol. 2017;26:937-50.

45. Bain A, Tzeng HY, Wu WJ, Chou LS. *Ficus* (Moraceae) and fig wasps (*Hymenoptera: Chalcidoidea*) in Taiwan. Bot Stud. 2015;56:11.

46. Sutton TL, DeGabriel JL, Riegler M, Cook JM. Local coexistence and genetic isolation of three pollinator species on the same fig tree species. Heredity. 2017; doi:10.1038/hdy.2016.125.

47. Kerdelhué C, Hochberg ME, Rasplus JY. Active pollination of *Ficus sur* by two sympatric fig wasp species in West Africa. Biotropica. 1997;29:69-75.

48. Tian E, Nason JD, Machado C, Zheng L, Yu H, Kjellberg F. Lack of genetic isolation by distance, similar genetic structuring but different demographic histories in a fig-pollinating wasp mutualism. Mol Ecol. 2015;24:5976-91.

49. Liu M, Zhang JI, Chen Y, Compton SG, Chen XY. Contrasting genetic responses to population fragmentation in a coevolving fig and fig wasp across a mainland-island archipelago. Mol Ecol. 2013;22:4384-96.

50. Liu M, Rui Zhao R, Chen Y, Zhang J, Compton SG, Chen XY. Competitive exclusion among fig wasps achieved via entrainment of host plant flowering phenology. PLoS One. 2014;9:e97783. doi:10.1371/journal.pone.0097783.

51. Kjellberg F, Proffit M. Tracking the elusive history of diversification in plant-herbivorous insect-parasitoid food webs: insights from figs and fig-wasps. Mol Ecol. 2016;25:843-5.

52. Sutton TL, Riegler M, Cook JM. One step ahead: a parasitoid disperses farther and forms a wider geographic population than its fig wasp host. Mol Ecol. 2016;25:882-94.

53. Brown RM, Siler CD, Oliveros CH, Esselstyn JA, Diesmos AC, Hosner PA, et al. Evolutionary processes of diversification in a model island archipelago. Annu Rev Ecol Evol Syst. 2013;44:411-35.

54. Esselstyn JA, Maher SP, Brown RM. Species interactions during diversification and community assembly in an island radiation of shrews. PLoS One. 2011; 6:e21885.

55. Azuma H, Harrison RD, Nakamura K, Su ZH. Molecular phylogenies of figs and fig-pollinating wasps in the Ryukyu and Bonin (Ogasawara) islands. Japan Genes Genet Syst. 2010;85:177-92.

56. Lin RC, Yeng CKL, Li SH. Drastic post-LGM expansion and lack of historical genetic structure of a subtropical fig-pollinating wasp (*Ceratosolen* sp. 1) of *Ficus septica* in Taiwan. Mol Ecol. 2008;17:5008-22.

57. Machado CA, Jousselin E, Kjellberg F, Compton SG, Herre EA. Phylogenetic relationships, historical biogeography and character evolution of fig-pollinating wasps. Proc R Soc Lond B. 2001;268:685-94.

Submit your next manuscript to BioMed Central
and we will help you at every step:

* We accept pre-submission inquiries

* Our selector tool helps you to find the most relevant journal

* We provide round the clock customer support

* Convenient online submission

* Thorough peer review

* Inclusion in PubMed and all major indexing services

* Maximum visibility for your research

Submit your manuscript at
www.biomedcentral.com/submit

BioMed Central