Cryptic diversity in a fig wasp community-morphologically differentiated species are sympatric but cryptic species are parapatric

Cryptic diversity begets challenges and opportunities in biodiversity research

How many species of life are there on Earth? This is a question that we want to know but cannot yet answer. Some scholars speculate that the number of species may reach 2.2 billion when considering cryptic diversity and that each morphology-based insect species may contain an average of 3.1 cryptic species. With nearly two million described species, such high estimates of cryptic diversity would suggest that cryptic species are widespread. The development of molecular species delimitation has led to the discovery of a large number of cryptic species, and cryptic biodiversity has gradually entered our field of vision and attracted more attention. This paper introduces the concept of cryptic species, how they evolve, and methods by which they may be discovered and confirmed, and provides theoretical and methodological guidance for the study of hidden species. A workflow of how to confirm cryptic species is provided. In addition, the importance and reliability of multi-evidence-based integrated taxonomy are reaffirmed as a way to better standardize decision-making processes. Special focus on cryptic diversity and increased funding for taxonomy is needed to ensure that cryptic species in hyperdiverse groups are discoverable and described. An increased focus on cryptic species in the future will naturally arise as more difficult groups are studied, and thereby, we may finally better understand the rules governing the evolution and maintenance of cryptic biodiversity.

Extreme specificity in obligate mutualism-A role for competition?

Obligate mutualisms, reciprocally obligate beneficial interactions, are some of the most important mutualisms on the planet, providing the basis for the evolution of the eukaryotic cell, the formation and persistence of terrestrial ecosystems and the establishment and expansion of coral reefs. In addition, these mutualisms can also lead to the diversification of interacting partner species. Accompanying this diversification is a general pattern of a high degree of specificity among interacting partner species. A survey of obligate mutualisms demonstrates that greater than half of these systems have only one or two mutualist species on each side of the interaction. This is in stark contrast to facultative mutualisms that can have dozens of interacting mutualist species. We posit that the high degree of specificity in obligate mutualisms is driven by competition within obligate mutualist guilds that limits species richness. Competition may be particularly potent in these mutualisms because mutualistic partners are totally dependent on each other's fitness gains, which may fuel interspecific competition. Theory and the limited number of empirical studies testing for the role of competition in determining specificity suggest that competition may be an important force that fuels the high degree of specificity. Further empirical research is needed to dissect the relative roles of trait complementarity, mutualism regulation, and competition among mutualist guild members in determining mutualism specificity at local scales.

Ecological Speciation without Morphological Differentiation? A New Cryptic Species of Diodontus Curtis (Hymenoptera, Pemphredonidae) from the Centre of Europe

Upon exploring the mitotype diversity of the aphid-hunting wasp, Diodontus tristis, we revealed specimens with highly divergent mitotypes from two localities in Lithuania and nesting in clayey substrate, while the specimens with typical mitotypes were found nesting in sandy sites. The comparison of inter- and intra-specific distances and application of delimitation algorithms supported the species status of the clay-nesting populations. Using a set of DNA markers that included complete or partial sequences of six mitochondrial genes, three markers of ribosomal operon, two homeobox genes, and four other nuclear genes, we clarified the phylogenetic relationships of the new cryptic species. The endosymbiotic bacteria infestation was checked, considering the option that the divergent populations may represent clades isolated by Wolbachia infection; however, it did not demonstrate any specificity. We found only subtle morphological differences in the new clay-nesting species, D. argillicola sp. nov.; the discriminant analysis of morphometric measurements did not reliably segregate it as well. Thus, we provide the molecular characters of the cryptic species, which allow confident identification, its phylogenetic position within the genus, and an updated identification key for the D. tristis species group.

New species based on the biological species concept within the complex of Lariophagus distinguendus (Hymenoptera, Chalcidoidea, Pteromalidae), a parasitoid of household pests

The pteromalid parasitoid Lariophagus distinguendus (Foerster) belongs to the Hymenoptera, a megadiverse insect order with high cryptic diversity. It attacks stored product pest beetles in human storage facilities. Recently, it has been shown to consist of two separate species. To further study its cryptic diversity, strains were collected to compare their relatedness using barcoding and nuclear genes. Nuclear genes identified two clusters which agree with the known two species, whereas the barcode fragment determined an additional third Clade. Total reproductive isolation (RI) according to the biological species concept (BSC) was investigated in crossing experiments within and between clusters using representative strains. Sexual isolation exists between all studied pairs, increasing from slight to strong with genetic distance. Postzygotic barriers mostly affected hybrid males, pointing to Haldane's rule. Hybrid females were only affected by unidirectional Spiroplasma-induced cytoplasmic incompatibility and behavioural sterility, each in one specific strain combination. RI was virtually absent between strains separated by up to 2.8% COI difference, but strong or complete in three pairs from one Clade each, separated by at least 7.2%. Apparently, each of these clusters represents one separate species according to the BSC, highlighting cryptic diversity in direct vicinity to humans. In addition, these results challenge the recent 'turbo-taxonomy' practice of using 2% COI differences to delimitate species, especially within parasitic Hymenoptera. The gradual increase in number and strength of reproductive barriers between strains with increasing genetic distance also sheds light on the emergence of barriers during the speciation process in L. distinguendus.

Intraspecies Genomic Divergence of a Fig Wasp Species Is Due to Geographical Barrier and Adaptation

Understanding how intraspecies divergence results in speciation has great importance for our knowledge of evolutionary biology. Here we applied population genomics approaches to a fig wasp species (Valisia javana complex sp 1) to reveal its intraspecies differentiation and the underlying evolutionary dynamics. With re-sequencing data, we prove the Hainan Island population (DA) of sp1 genetically differ from the continental ones, then reveal the differed divergence pattern. DA has reduced SNP diversity but a higher proportion of population-specific structural variations (SVs), implying a restricted gene exchange. Based on SNPs, 32 differentiated islands containing 204 genes were detected, along with 1,532 population-specific SVs of DA overlapping 4,141 genes. The gene ontology (GO) enrichment analysis performed on differentiated islands linked to three significant GO terms on a basic metabolism process, with most of the genes failing to enrich. In contrast, population-specific SVs contributed more to the adaptation than the SNPs by linking to 59 terms that are crucial for wasp speciation, such as host reorganization and development regulation. In addition, the generalized dissimilarity modeling confirms the importance of environment difference on the genetic divergence within sp1. Hence, we assume the genetic divergence between DA and the continent due to not only the strait as a geographic barrier, but also adaptation. We reconstruct the demographic history within sp1. DA shares a similar population history with the nearby continental population, suggesting an incomplete divergence. Summarily, our results reveal how geographic barriers and adaptation both influence the genetic divergence at population-level, thereby increasing our knowledge on the potential speciation of non-model organisms.

Source-sink dynamics assists the maintenance of a pollinating wasp

Dispersal that unites spatially subdivided populations into a metapopulation with source-sink dynamics is crucial for species persistence in fragmented landscapes. Understanding such dynamics for pollinators is particularly urgent owing to the ongoing global pollination crisis. Here, we investigated the population structure and source-sink dynamics of a pollinating wasp (Wiebesia sp. 3) of Ficus pumila in the Zhoushan Archipelago of China. We found significant asymmetry in the pairwise migrant numbers for 22 of 28 cases on the historical timescale, but only two on the contemporary timescale. Despite a small population size, the sole island not colonized by a superior competitor wasp (Wiebesia sp. 1) consistently behaved as a net exporter of migrants, supplying large sinks. Comparable levels of genetic diversity, with few private alleles and low genetic differentiation (total F_st : 0.03; pairwise F_st : 0.0005-0.0791), were revealed among all the islands. There was a significant isolation-by-distance pattern caused mainly by migration between the competition-free island and other islands, otherwise the pattern was negligible. The clustering analysis failed to detect multiple gene pools for the whole region. Thus, the sinks were most probably organized into a patchy population. Moreover, the estimates of effective population sizes were comparable between the two timescales. Thus the source-sink dynamics embedded within a well-connected population network may allow Wiebesia sp. 3 to persist at a competitive disadvantage. This study provides evidence that metapopulations in the real world may be complicated and changeable over time, highlighting the necessity to study such metapopulations in detail.

Low host specificity and broad geographical ranges in a community of parasitic non-pollinating fig wasps (Sycoryctinae; Chalcidoidea)

Plants, phytophagous insects and their parasitoids form the most diverse assemblages of macroscopic organisms on earth. Enclosed assemblages in particular represent a tractable system for studying community assembly and diversification. Communities associated with widespread plant species are especially suitable as they facilitate a comparative approach. Pantropical fig-wasp communities represent a remarkably well-replicated system, ideal for studying these historical processes. We expect high dispersal ability in non-pollinating fig wasps to result in lower geographical turnover in comparison to pollinating fig wasps. The ability of non-pollinating wasps to utilise a number of hosts (low host specificity) is a key determinant of overall geographical range, with intraspecific competition becoming a constraining factor should diet breadth overlap among species. Finally, we expect conserved community structure throughout the host range. We aim to test these expectations, derived from population genetic and community studies, using the multi-trophic insect community associated with Ficus hirta throughout its 3,500 km range across continental and insular Asia. We collect molecular evidence from one coding mitochondrial gene, one non-coding nuclear gene and multiple microsatellites across 25 geographical sites. Using these data, we establish species boundaries, determine levels of host specificity among non-pollinating fig wasps and quantify geographical variation in community composition. We find low host specificity in two genera of non-pollinating fig wasps. Functional community structure is largely conserved across the range of the host fig, despite limited correspondence between the ranges of non-pollinator and pollinator species. While nine pollinators are associated with Ficus hirta, the two non-pollinator tribes developing in its figs each contained only four species. Contrary to predictions, we find stronger isolation by distance in non-pollinators than pollinators. Long-lived non-pollinators may disperse more gradually and be less reliant on infrequent long-distance dispersal by wind currents. Segregation among non-pollinating species across their range is suggestive of competitive exclusion and we propose that this may be a result of increased levels of local adaptation and moderate, but regular, rates of dispersal. Our findings provide one more example of lack of strict codiversification in the geographical diversification of plant-associated insect communities.

A Comparative Test for Divergent Adaptation: Inferring Speciation Drivers from Functional Trait Divergence

AbstractEcological differentiation between lineages is widely considered to be an important driver of speciation, but support for this hypothesis is mainly derived from the detailed study of a select set of model species pairs. Mounting evidence from nonmodel taxa, meanwhile, suggests that speciation often occurs with minimal differentiation in ecology or ecomorphology, calling into question the true contribution of divergent adaptation to species richness in nature. To better understand divergent ecological adaptation and its role in speciation generally, researchers require a comparative approach that can distinguish its signature from alternative processes, such as drift and parallel selection, in data sets containing many species pairs. Here we introduce new statistical models of divergent adaptation in the continuous traits of paired lineages. In these models, ecomorphological characters diverge as two lineages adapt toward alternative phenotypic optima following their departure from a common ancestor. The absolute distance between optima measures the extent of divergent selection and provides a basis for interpretation. We encode the models in the new R package diverge and extend them to allow the distance between optima to vary across continuous and categorical variables. We test model performance using simulation and demonstrate model application using published data sets of trait divergence in birds and mammals. Our framework provides the first explicit test for signatures of divergent selection in trait divergence data sets, and it will enable empiricists from a wide range of fields to better understand the dynamics of divergent adaptation and its prevalence in nature beyond just our best-studied model systems.

A morphometric assessment of species boundaries in a widespread anole lizard (Squamata: Dactyloidae)

Cryptic species – genetically distinct species that are morphologically difficult to distinguish – present challenges to systematists. Operationally, cryptic species are very difficult to identify and sole use of genetic data or morphological data can fail to recognize evolutionarily isolated lineages. We use morphometric data to test species boundaries hypothesized with genetic data in the North Caribbean bark anole (Anolis distichus), a suspected species complex. We use univariate and multivariate analyses to test if candidate species based on genetic data can be accurately diagnosed. We also test alternative species delimitation scenarios with a model fitting approach that evaluates normal mixture models capable of identifying morphological clusters. Our analyses reject the hypothesis that the candidate species are diagnosable. Neither uni- nor multivariate morphometric data distinguish candidate species. The best-supported model included two morphological clusters; however, these clusters were uneven and did not align with a plausible species divergence scenario. After removing two related traits driving this result, only one cluster was supported. Despite substantial differentiation revealed by genetic data, we recover no new evidence to delimit species and refrain from taxonomic revision. This study highlights the importance of considering other types of data along with molecular data when delimiting species.

Sky islands as foci for divergence of fig trees and their pollinators in southwest China

The dynamics of populations and their divergence over time have shaped current levels of biodiversity and in the case of the "sky islands" of mountainous southwest (SW) China have resulted in an area of exceptional botanical diversity. Ficus tikoua is a prostrate fig tree subendemic to the area that displays unique intraspecific diversity, producing figs typical of different pollination modes in different parts of its range. By combining climate models, genetic variation in populations of the tree's obligate fig wasp pollinators and distributions of the different plant phenotypes, we examined how this unusual situation may have developed. We identified three genetically distinct groups of a single Ceratosolen pollinator species that have largely parapatric distributions. The complex topography of the region contributed to genetic divergence among the pollinators by facilitating geographical isolation and providing refugia. Migration along elevations in response to climate oscillations further enhanced genetic differentiation of the three pollinator groups. Their distributions loosely correspond to the distributions of the functionally significant morphological differences in the male figs of their host plants, but postglacial expansion of one group has not been matched by spread of its associated plant phenotype, possibly due to a major river barrier. The results highlight how interplay between the complex topography of the "sky island" complex and climate change has shaped intraspecies differentiation and relationships between the plant and its pollinator. Similar processes may explain the exceptional botanical diversity of SW China.

Patterns of morphological diversification in giant Berberomeloe blister beetles (Coleoptera: Meloidae) reveal an unexpected taxonomic diversity concordant with mtDNA phylogenetic structure

Delimiting species boundaries is a complex challenge usually hindered by overlooked morphological diversification or misinterpretation of geographically structured phenotypic variability. Independent molecular data are extremely useful to characterize and understand such morphological diversity. Morphological and molecular variability of the non-phoretic and apterous, widely distributed, giant blister beetles of the genus Berberomeloe, were investigated within and between lineages across most of the distributional range of the genus. We used two mtDNA gene fragments to characterize genetic variability and to produce a time-calibrated phylogeny of the genus. Our results reveal several mitochondrial lineages, allopatrically, parapatrically and sympatrically distributed. Most clades are not distinguishable between each other based on morphometrics. However, no morphometric overlap is observed between two closely related clades, one of them occurring in sympatry with a distantly congeneric species (B. insignis), suggesting that sympatry could trigger morphological diversification. Although most species share a morphometric space, they can be morphologically identified by a combination of easily observed characteristic qualitative features. Based on the concordance between mtDNA clades and morphological units, we describe six new species of Berberomeloe (B. castuo sp. nov., B. comunero sp. nov., B. indalo sp. nov, B. yebli sp. nov., B. payoyo sp. nov. and B. tenebrosus sp. nov.), revalidate two taxa (B. maculifrons comb. nov. and B. laevigatus comb. nov.) and redefine B. majalis.

Cuticular hydrocarbons as potential mediators of cryptic species divergence in a mutualistic ant association

Upon advances in sequencing techniques, more and more morphologically identical organisms are identified as cryptic species. Often, mutualistic interactions are proposed as drivers of diversification. Species of the neotropical parabiotic ant association between Crematogaster levior and Camponotus femoratus are known for highly diverse cuticular hydrocarbon (CHC) profiles, which in insects serve as desiccation barrier but also as communication cues. In the present study, we investigated the association of the ants' CHC profiles with genotypes and morphological traits, and discovered cryptic species pairs in both genera. To assess putative niche differentiation between the cryptic species, we conducted an environmental association study that included various climate variables, canopy cover, and mutualistic plant species. Although mostly sympatric, the two Camponotus species seem to prefer different climate niches. However in the two Crematogaster species, we could not detect any differences in niche preference. The strong differentiation in the CHC profiles may thus suggest a possible role during speciation itself either by inducing assortative mating or by reinforcing sexual selection after the speciation event. We did not detect any further niche differences in the environmental parameters tested. Thus, it remains open how the cryptic species avoid competitive exclusion, with scope for further investigations.

Loss of top-down biotic interactions changes the relative benefits for obligate mutualists

The collapse of mutualisms owing to anthropogenic changes is contributing to losses of biodiversity. Top predators can regulate biotic interactions between species at lower trophic levels and may contribute to the stability of such mutualisms, but they are particularly likely to be lost after disturbance of communities. We focused on the mutualism between the fig tree Ficus microcarpa and its host-specific pollinator fig wasp and compared the benefits accrued by the mutualists in natural and translocated areas of distribution. Parasitoids of the pollinator were rare or absent outside the natural range of the mutualists, where the relative benefits the mutualists gained from their interaction were changed significantly away from the plant's natural range owing to reduced seed production rather than increased numbers of pollinator offspring. Furthermore, in the absence of the negative effects of its parasitoids, we detected an oviposition range expansion by the pollinator, with the use of a wider range of ovules that could otherwise have generated seeds. Loss of top-down control has therefore resulted in a change in the balance of reciprocal benefits that underpins this obligate mutualism, emphasizing the value of maintaining food web complexity in the Anthropocene.

Two's company, three's a crowd: co-occurring pollinators and parasite species in Breynia oblongifolia (Phyllanthaceae)

BACKGROUND

Obligate pollination mutualisms (OPMs) are specialized interactions in which female pollinators transport pollen between the male and female flowers of a single plant species and then lay eggs into those same flowers. The pollinator offspring hatch and feed upon some or all of the developing ovules pollinated by their mothers. Strong trait matching between plants and their pollinators in OPMs is expected to result in reciprocal partner specificity i.e., a single pollinator species using a single plant species and vice versa, and strict co-speciation. These issues have been studied extensively in figs and fig wasps, but little in the more recently discovered co-diversification of Epicephala moths and their Phyllanthaceae hosts. OPMs involving Epicephala moths are believed occur in approximately 500 species of Phyllanthaceae, making it the second largest OPM group after the Ficus radiation (> 750 species). In this study, we used a mixture of DNA barcoding, genital morphology and behavioral observations to determine the number of Epicephala moth species inhabiting the fruits of Breynia oblongifolia, their geographic distribution, pollinating behavior and phylogenetic relationships.

RESULTS

We found that B. oblongifolia hosts two species of pollinator that co-occurred at all study sites, violating the assumption of reciprocal specificity. Male and female genital morphologies both differed considerably between the two moth species. In particular, females differed in the shape of their ovipositors, eggs and oviposition sites. Phylogenetic analyses indicated that the two Epicephala spp. on B. oblongifolia likely co-exist due to a host switch. In addition, we discovered that Breynia fruits are also often inhabited by a third moth, an undescribed species of Herpystis, which is a non-pollinating seed parasite.

CONCLUSIONS

Our study reveals new complexity in interactions between Phyllantheae and Epicephala pollinators and highlights that host switching, co-speciation and non-pollinating seed parasites can shape species interactions in OPMs. Our finding that co-occurring Epicephala species have contrasting oviposition modes parallels other studies and suggests that such traits are important in Epicephala species coexistence.

Genomic data reject the hypothesis of sympatric ecological speciation in a clade of Desmognathus salamanders

Closely related taxa with dissimilar morphologies are often considered to have diverged via natural selection favoring different phenotypes. However, some studies have found these scenarios to be paired with limited or no genetic differentiation. Desmognathus quadramaculatus and D. marmoratus are sympatric salamander species thought to represent a case of ecological speciation based on distinct morphologies, but the results of previous studies have not resolved corresponding patterns of lineage divergence. Here, we use genome-wide data to test this hypothesis of ecological speciation. Population structure analyses partitioned individuals geographically, but not morphologically, into two adjacent regions of western North Carolina: Pisgah and Nantahala. Phylogenetic analyses confirmed the nominal species are nonmonophyletic and resolved deep divergence between the two geographic clusters. Model-testing overwhelmingly supported the hypothesis that lineage divergence followed geography. Finally, ecological niche modeling showed that Pisgah and Nantahala individuals occupy different climatic niches, and geographic boundaries for the two lineages correspond to differences in precipitation regimes across southern Appalachia. Overall, we reject the previous hypothesis of ecological speciation based on microhabitat partitioning. Instead, our results suggest that there are two cryptic lineages, each containing the same pair of morphotypes.

Plant geographic phenotypic variation drives diversification in its associated community of a phytophagous insect and its parasitoids

BACKGROUND

While the communities constituted by phytophageous insects and their parasites may represent half of all terrestrial animal species, understanding their diversification remains a major challenge. A neglected idea is that geographic phenotypic variation in a host plant may lead to heterogeneous evolutionary responses of the different members of the associated communities. This could result in diversification on a host plant by ecological speciation in some species, leading to geographic variation in community composition. In this study we investigated geographic variation of inflorescence receptacle size in a plant, Ficus hirta, and how the hymenopteran community feeding in the inflorescences has responded. Our predictions were: 1) Inflorescence size variation affects wasp species differently depending on how they access oviposition sites. 2) In some affected lineages of wasps, we may observe vicariant, parapatric species adapted to different inflorescence sizes.

RESULTS

We show that fig (the enclosed inflorescence of Ficus) wall thickness varies geographically. The fig-entering pollinating wasp was not affected, while the parasites ovipositing through the fig wall were. Two parapatric species of Philotrypesis, exhibiting strikingly different ovipositor lengths, were recorded. One species of Sycoscapter was also present, and it was restricted, like the shorter-ovipositor Philotrypesis, to the geographic zone where fig walls were thinner.

CONCLUSIONS

Previous work on fig wasps suggested that parapatric geographic ranges among congenerics were due to adaptation to variation in abiotic factors, complemented by interspecific competition. Our results show that parapatric ranges may also result from adaptation to variation in biotic factors. Within an insect community, differences among species in their response to geographic phenotypic variation of their host plant may result in geographically heterogeneous community structure. Such heterogeneity leads to heterogeneous interaction networks among sites. Our results support the hypothesis that plant geographic phenotypic variation can be a driver of diversification in associated insect communities, and can complement other diversification processes.

Conserved community structure and simultaneous divergence events in the fig wasps associated with Ficus benjamina in Australia and China

Background

Localised patterns of species diversity can be influenced by many factors, including regional species pools, biogeographic features and interspecific interactions. Despite recognition of these issues, we still know surprisingly little about how invertebrate biodiversity is structured across geographic scales. In particular, there have been few studies of how insect communities vary geographically while using the same plant host. We compared the composition (species, genera) and functional structure (guilds) of the chalcid wasp communities associated with the widespread fig tree, Ficus benjamina, towards the northern (Hainan province, China) and southern (Queensland, Australia) edges of its natural range. Sequence data were generated for nuclear and mtDNA markers and used to delimit species, and Bayesian divergence analyses were used to test patterns of community cohesion through evolutionary time.

Results

Both communities host at least 14 fig wasp species, but no species are shared across continents. Community composition is similar at the genus level, with six genera shared although some differ in species diversity between China and Australia; a further three genera occur in only China or Australia. Community functional structure remains very similar in terms of numbers of species in each ecological guild despite community composition differing a little (genera) or a lot (species), depending on taxonomic level. Bayesian clustering analyses favour a single community divergence event across continents over multiple events for different ecological guilds. Molecular dating estimates of lineage splits between nearest inter-continental species pairs are broadly consistent with a scenario of synchronous community divergence from a shared “ancestral community”.

Conclusions

Fig wasp community structure and genus-level composition are largely conserved in a wide geographic comparison between China and Australia. Moreover, dating analyses suggest that the functional community structure has remained stable for long periods during historic range expansions. This suggests that ecological interactions between species may play a persistent role in shaping these communities, in contrast to findings in some comparable temperate systems.

Electronic supplementary material

The online version of this article (10.1186/s12898-018-0167-y) contains supplementary material, which is available to authorized users.

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

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.