Hui Yu1

Enwei Tian1

Linna Zheng1

Xiaoxia Deng1

Yufen Cheng1

Lianfu Chen1

Wei Wu1

Wattana Tanming2

Dayong Zhang3

Stephen G Compton4

Finn Kjellberg5

Hui Yu
email: yuhui@scib.ac.cn

Multiple parapatric pollinators have radiated across a continental fig tree displaying clinal genetic variation

Keywords: Co-evolution, dispersal, gene flow, mutualism, plant growth form, speciation

The ways that plant-feeding insects have diversified are central to our understanding of terrestrial ecosystems. Obligate nursery pollination mutualisms provide highly relevant model systems of how plants and their insect associates have diversified and the over 800 species of fig trees (Ficus) allow comparative studies. Fig trees can have one or more pollinating fig wasp species (Agaonidae) that breed within their figs, but factors influencing their number remain to be established. In some widely distributed fig trees, the plants form populations isolated by large swathes of sea, and the different populations are pollinated by different wasp species.

Other Ficus species with continuous distributions may present genetic signatures of isolation by distance, suggesting more limited pollinator dispersal, which may also facilitate pollinator speciation. We tested the hypothesis that Ficus hirta, a species for which preliminary data showed genetic isolation by distance, would support numerous pollinator species across its range. Our results show that across its range F. hirta displays clinal genetic variation and is pollinated by nine parapatric species of Valisia. This is the highest number of pollinators reported to date for any Ficus species, and it is the first demonstration of the occurrence of parapatric pollinator species on a fig host displaying continuous genetic structure. Future comparative studies across Ficus species should be able to establish the plant traits that have driven the evolution of pollinator dispersal behaviour, pollinator speciation and host plant spatial genetic structure.

INTRODUCTION

Insects feeding on plants represent one of the most common forms of trophic interaction seen in terrestrial ecosystems, and the relationship between plants and insects has generated the most species -rich animal taxa ever seen on the planet (Price, 1980[START_REF] Price | Species interactions and the evolution of biodiversity[END_REF]. These insects are typically associated with one or a few related host plants, so different plants host different insects [START_REF] Strong | Insects on plants: Community patterns and mechanisms[END_REF]. Processes such as local adaptation, the geography of speciation and the consequences of secondary sympatry are key issues in understanding the dynamics of community diversification [START_REF] Mittelbach | Ecological and evolutionary perspectives on community assembly[END_REF]. However, almost no studies have explored how the insect community on a plant changes across its geographic range [START_REF] Lawton | Patterns of diversity for the insect herbivores on bracken[END_REF][START_REF] Leather | Insect species richness of the British Rosaceae: The importance of host range, plant architecture, age of establishment, taxonomic isolation and species -area relationships[END_REF][START_REF] Lewinsohn | Four ways towards tropical herbivore megadiversity[END_REF]. A study across Papuan rainforests found little variation in herbivorous insect communities over distances of up to 1,000 km [START_REF] Craft | Population genetics of ecological communities with DNA barcodes: An example from New Guinea Lepidoptera[END_REF][START_REF] Novotny | Low beta diversity of herbivorous insects in tropical forests[END_REF], while a study of nonpollinating fig wasps on Ficus rubiginosa Desf. ex Ventenat in Australia detected, within a set of seven widespread morphospecies, four pairs of parapatric cryptic species, along a 2,000km transect [START_REF] Darwell | Cryptic diversity in a fig wasp community-Morphologically differentiated species are sympatric but cryptic species exhibit competitive exclusion[END_REF].

In nursery pollination mutualisms, pollinating insects breed in floral structures of the plants they pollinate [START_REF] Dufaÿ | Conflicts between plants and pollinators that reproduce within inflorescences: Evolutionary variations on a theme[END_REF]. Insects involved in nursery pollination mutualisms constitute particularly favourable biological models for investigating the geographic variation of diversity on a host plant as the insect's prevalence on a host is high and their presence is easy to detect. Furthermore, some of these mutualisms are species rich allowing comparative studies. Nursery pollination systems, such as those involving Yucca and Yucca moths [START_REF] Pellmyr | Yuccas, yucca moths, and coevolution: A review[END_REF], Glochidion and Epicephala moths [START_REF] Kawakita | Evolution of obligate pollination mutualism in the tribe Phyllantheae (Phyllanthaceae)[END_REF] and fig trees (Ficus) and fig wasps (Agaonidae) [START_REF] Cook | Mutualists with attitude: Coevolving fig wasps and figs[END_REF], were initially envisioned as systems in which a single species of pollinating insect breeds in and pollinates a single host plant species. This simple pattern is now largely rejected because of accumulating examples where one insect species pollinates several host species [START_REF] Cornille | Floral volatiles, pollinator sharing and diversification in the fig -wasp mutualism: insights from Ficus natalensis, and its two wasp pollinators (South Africa)[END_REF][START_REF] Hembry | Does biological intimacy shape ecological network structure -A test using a brood pollination mutualism on continental and oceanic islands[END_REF][START_REF] Pellmyr | Systematic revision of the yucca moths in the Tegeticula yuccasella complex (Lepidoptera: Prodoxidae) north of Mexico[END_REF][START_REF] Wachi | Genome -wide sequence data suggest the possibility of pollinator sharing by host shift in dioecious figs (Moraceae, Ficus)[END_REF][START_REF] Wang | Pollinator sharing and gene flow among closely related sympatric dioecious fig taxa[END_REF], or cases where several insect species pollinate a single host [START_REF] Darwell | Molecular species delimitation of a symbiotic fig -pollinating wasp species complex reveals extreme deviation from reciprocal partner specificity[END_REF][START_REF] Li | Four new species of Epicephala Meyrick, 1880 (Lepidoptera, Gracillariidae) associated with two species of Glochidion (Phyllanthaceae) from Hainan Island in China[END_REF][START_REF] Yang | The incidence and pattern of copollinator diversification in dioecious and monoecious figs[END_REF]. A general model of plant-insect diversification in these nursery pollination mutualisms must therefore take into account the factors influencing their relative rates of speciation and explain why some of the plants have a single pollinator when others have two or more.

Fig trees (Ficus) represent the most species-rich group of plants offering nursery pollination rewards and the most diversified lineage within the Moraceae (Brunn -Lund, Verstraete, [START_REF] Bruun -Lund | Rush hour at the Museum -Diversification patterns provide new clues for the success of figs (Ficus L., Moraceae)[END_REF]. only transfer pollen from those plants in which their larvae can develop. This is a more direct link between suitability for offspring development and subsequent pollen transfer than in any other brood pollination mutualism and may favour co-adaptation [START_REF] Anstett | Figs and fig pollinators: Evolutionary conflicts in a coevolved mutualism[END_REF]. However, a more rapid rate of speciation among fig wasps, compared with their hosts, is to be expected given their much shorter generation times [START_REF] Petit | Some evolutionary consequences of being a tree[END_REF][START_REF] Thomas | A generation time effect on the rate of molecular evolution in invertebrates[END_REF]. From this perspective, the long -held assumption that each fig tree was pollinated by a single species of pollinator was problematic.

More extensive sampling and molecular techniques that have facilitated separation of morphologically close ("cryptic") species have revealed an increasing number of Ficus species with several pollinators [START_REF] Darwell | Molecular species delimitation of a symbiotic fig -pollinating wasp species complex reveals extreme deviation from reciprocal partner specificity[END_REF]. However, we still know little about how genetic diversity is structured within Ficus species and how this is related to the distributions and genetic diversity among the fig wasp pollinators that it supports [START_REF] Bain | Geographic structuring into vicariant species -pairs in a wide -ranging, high -dispersal plant -insect mutualism: The case of Ficus racemosa and its pollinating wasps[END_REF][START_REF] Rodriguez | Diversification and spatial structuring in the mutualism between Ficus septica and its pollinating wasps in insular South East Asia[END_REF][START_REF] Wachi | Genome -wide sequence data suggest the possibility of pollinator sharing by host shift in dioecious figs (Moraceae, Ficus)[END_REF]. Where several pollinator species share a single Ficus host, they are generally sister species [START_REF] Yang | The incidence and pattern of copollinator diversification in dioecious and monoecious figs[END_REF]. The recorded exceptions, where nonsister pollinator species occupy the same host species, have almost all been within two Ficus clades that have diversified recently, and this appears to have favoured host shifts [START_REF] Cruaud | An extreme case of plant-insect codiversification: Figs and figpollinating wasps[END_REF][START_REF] Jousselin | One fig to bind them all: Host conservatism in a fig wasp community unravelled by cospeciation analyses among pollinating and nonpollinating fig wasps[END_REF][START_REF] Machado | Critical review of host specificity and its coevolutionary implications in the fig fig -wasp mutualism[END_REF][START_REF] Yang | The incidence and pattern of copollinator diversification in dioecious and monoecious figs[END_REF]. The main modes of fig wasp speciation therefore appear to vary across lineages.

Widely distributed plants grow in areas with a wide range of environmental conditions, which will vary in suitability for their associated insects, including their pollinators. Reflecting this, plant -eating insects usually only occupy a subset of the ranges of their hosts [START_REF] Strong | Insects on plants: Community patterns and mechanisms[END_REF], but fig trees require an effective and specific pollinator to be present wherever they grow. Widely distributed Ficus species would therefore be predicted to support additional species of pollinator compared to species with narrow, more climatically homogeneous, ranges. Only a few widely distributed Ficus species have been sampled at multiple sites across their geographic ranges. Ficus racemosa L. is pollinated by a single fig wasp across India and by another in China and SE Asia, with additional sister species of pollinators in the island of Borneo and Australia [START_REF] Bain | Geographic structuring into vicariant species -pairs in a wide -ranging, high -dispersal plant -insect mutualism: The case of Ficus racemosa and its pollinating wasps[END_REF][START_REF] Kobmoo | Ficus racemosa is pollinated by a single population of a single agaonid wasp species in continental South -East Asia[END_REF]. More localized differentiation is evident among island populations of Ficus septica Burm. f., in the Philippines [START_REF] Conchou | 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[END_REF][START_REF] Lin | The lack of pollinator specificity in a dioecious fig tree: Sympatric figpollinating wasps of Ficus septica in southern Taiwan[END_REF][START_REF] Rodriguez | Diversification and spatial structuring in the mutualism between Ficus septica and its pollinating wasps in insular South East Asia[END_REF], and along elevation gradients [START_REF] Segar | Speciation in a keystone plant genus is driven by elevation: A case study in New Guinean Ficus[END_REF]SoutoVilarós et al., 2018), where several different pollinators are present.

The absence of divergence among the pollinators of F. racemosa across large swathes of its continental range is likely to reflect their ability to disperse over large distances. Spatial genetic structuring across the range of a plant species provides an indication of the geographic extent of their gene flow, and in the case of Ficus, it has revealed a widely varying extent of dispersal among their pollinators. Paternity analyses have shown that some fig wasps transport pollen between large fig trees growing over 100 km apart [START_REF] Ahmed | Windborne insects mediate directional pollen transfer between desert fig trees 160 kilometers apart[END_REF], whereas population structuring shows that other, smaller, species have much more localized gene flow [START_REF] Chen | Genetic diversity and differentiation of the extremely dwarf Ficus tikoua in Southwestern China[END_REF][START_REF] Liu | Movements of genes between populations: Are pollinators more effective at transferring their own or plant genetic markers[END_REF]. These differences in dispersal behaviour among the pollinators of different trees appear highly likely to influence differentiation and speciation of both their own populations as well as those of the host plants they pollinate.

Previous data have shown that Ficus hirta Vahl presents a pattern of spatial genetic structure suggesting genetic isolation by distance across continental South -East Asia [START_REF] Yu | Nuclear and chloroplast DNA phylogeography of Ficus hirta: Obligate pollination mutualism and constraints on range expansion in response to climate change[END_REF]. Here, based on extensive sampling, we describe the population genetic structure of F.

MATERIALS AND METHODS

Study species

Like all Ficus species, F. hirta (family Moraceae, subgenus Ficus, section Eriosycea, [START_REF] Berg | Flora Malesiana precursor for the treatment of Moraceae 3: Ficus subgenus Ficus[END_REF] In contrast, the figs on female trees develop seeds and no pollinator offspring. Ficus hirta is a species of secondary, disturbed habitats and typically grows as a shrub or treelet. It can reach a maximum height of about 5 m, but most mature individuals are much shorter [START_REF] Berg | Moraceae (Ficus)[END_REF]. Its figs are small and produced in the leaf axils. Female figs ripen to a bright red and are mainly eaten by birds [START_REF] Corlett | Figs (Ficus, Moraceae) in Urban Hong Kong, South China[END_REF]. As many as 50 figs can be present on a tree, but their development is asynchronous, which reduces the peak numbers of figs available for pollination or dispersal at any one time and can also allow pollinator cycling between figs on the same male tree [START_REF] Jia | jMOTU and taxonerator: Turning DNA barcode sequences into annotated operational taxonomic units[END_REF].

The distribution of F. hirta extends northwards from the island of Java (Indonesia) in the south to China in the north and westwards into northeast India, Nepal and Sikkim [START_REF] Berg | Moraceae (Ficus)[END_REF]. This extensive geographic range (covering over 30° of latitude and 30° of longitude) encompasses a range of tropical and subtropical biomes. [START_REF] Berg | Precursory taxonomic studies on Ficus (Moraceae) for the Flora of Thailand[END_REF] restricted F. hirta to what was previously known as F. h. subsp. hirta. In the following, we will use the name F. hirta for F. hirta sensu [START_REF] Berg | Precursory taxonomic studies on Ficus (Moraceae) for the Flora of Thailand[END_REF]. Following that definition, F. hirta has a disjunct distribution. South of Thailand, it is absent from Peninsular Malaysia and it is only present in South Sumatra and North Java [START_REF] Berg | Moraceae (Ficus)[END_REF]. Population genetic studies, mainly from China, detected some spatial genetic structuring across F. hirta populations and concluded that low nuclear differentiation, combined with high interpopulation differentiation and geographic structuring of chloroplast variation, indicated that gene flow via seed dispersal was more limited than via dispersal of pollen [START_REF] Yu | Nuclear and chloroplast DNA phylogeography of Ficus hirta: Obligate pollination mutualism and constraints on range expansion in response to climate change[END_REF]. Genetic differentiation in nuclear genes between populations of F. hirta on the Chinese mainland and Hainan Island has also been detected, and corresponding differences were also present between the populations of the tree's pollinators [START_REF] Tian | Lack of genetic isolation by distance, similar genetic structuring but different demographic histories in a fig -pollinating wasp mutualism[END_REF]. Two morphologically distinguishable pollinator fig wasps have been described from the figs of F. hirta. Valisia javana javana Mayr was reared from F. hirta figs collected in Java, and a second subspecies, Valisia j. hilli Wiebes, was described from figs collected in Hong Kong [START_REF] Wiebes | Agaonidae (Hymenoptera Chalcidoidea) and Ficus (Moraceae): Fig wasps and their figs, xi (Blastophaga) s.l[END_REF].

Sampling

Between 2006 and 2014, we sampled the leaves of F. hirta (31 locations) and its pollinating fig wasps (32 locations) across mainland South -East Asia and south to the island of Java (Supporting Information Table S1). The samples were up to 4,100 km apart (north-south) and reached the northern, eastern and southern limits of the range of F. hirta.

Within locations, F. hirta individuals were typically sampled at intervals of 3-5 m, with no individuals located more than 200 m apart. Their leaves were collected and dried in plastic bags containing silica gel.

Pollinating wasps 2.3.1 Pollinating wasp DNA extraction, amplification and sequencing

The mitochondrial genetic marker mtCOI was sequenced from an average of 9.9 individuals per location (range 1-28, total 330), and microsatellite loci were genotyped from an average of 18.9 other individuals from the same locations (range 2-32, total 568). To complement the data, the nuclear ITS2 nuclear gene was amplified for 201 individuals.

Genomic DNA was extracted from the whole body of each fig wasp using the EasyPure Genomic DNA Extraction Kit (TransGen, Beijing, China). A 681 bp fragment of the mtCOI gene was then sequenced following the protocol used in previous studies [START_REF] Tian | Lack of genetic isolation by distance, similar genetic structuring but different demographic histories in a fig -pollinating wasp mutualism[END_REF]. A 689 bp fragment of the ITS2 gene was amplified in 201 individuals using the universal primer pair (ITSR -5 --CGCCTGATCTGAGGTCGTGA -3 --ITSF -5 -ATTCCCGCACCACGCCTGGCTGA -3 -; Lopez -Vaamonde, Rasplus, Weiblen, & [START_REF] Lopez -Vaamonde | Molecular phylogenies of fig wasps: Partial cocladogenesis of pollinators and parasites[END_REF] and the same PCR amplification reaction volume as for the COI gene. The reaction was optimized and programmed on a MJ Thermal Cycler (PTC 200) as one cycle of denaturation at 94°C for 5 min, 35 cycles of 30 -s denaturation at 94°C, 30 s at a 55°C annealing temperature and 30 -s extension at 72°C, followed by 8 -min extension at 72°C. All amplified PCR products were purified using QIAquick spin columns (Qiagen) and were sequenced in an ABI 3730xl capillary sequencer using bigdye terminator V 3.1 chemistry (Applied Biosystems).

Previously published protocols were also used to genotype individuals at nine unlinked microsatellite loci (1-78, 1-141, A34, A80, A99, B30, C25, F17 and H33) that had been previously developed for Valisia j. hilli [START_REF] Tian | Development of microsatellite loci for Blastophaga javana (Agaonidae), the pollinating wasp of Ficus hirta (Moraceae)[END_REF]. The alleles used in the present study were the same as those used previously to analyse southeast China and Hainan populations [START_REF] Tian | Lack of genetic isolation by distance, similar genetic structuring but different demographic histories in a fig -pollinating wasp mutualism[END_REF].

Pollinating fig wasp sequence analysis

We did not detect any indications of pseudogenes, such as multiple peaks in chromatograms, stop codons or frameshift mutations [START_REF] Song | Many species in one: DNA barcoding overestimates the number of species when nuclear mitochondrial pseudogenes are coamplified[END_REF]. Sequences were aligned using MUSCLE [START_REF] Edgar | MUSCLE: Multiple sequence alignment with high accuracy and high throughput[END_REF] implemented in mega 6.0 [START_REF] Tamura | MEGA6: Molecular evolutionary genetics analysis version 6.0[END_REF] with manual corrections.

Previous work has shown that the fig wasp genus Valisia, including samples collected from F. hirta, is monophyletic [START_REF] Cruaud | Laying the foundations for a new classification of Agaonidae (Hymenoptera: Chalcidoidea), a multilocus phylogenetic approach[END_REF]. We downloaded from GenBank the COI sequences of the eight Valisia species available including one sequence (FJ619191) of V. javana. Published phylogenies have shown that the genera Ceratosolen and Kradibia constitute an outgroup relative to all other pollinating fig wasps [START_REF] Cruaud | Laying the foundations for a new classification of Agaonidae (Hymenoptera: Chalcidoidea), a multilocus phylogenetic approach[END_REF][START_REF] Cruaud | An extreme case of plant-insect codiversification: Figs and figpollinating wasps[END_REF]. We included representatives of these two genera (two species of Ceratosolen and 11 species of Kradibia) in the phylogenetic analysis.

Dated phylogenetic trees that included the downloaded sequences and all our sequenced haplotypes were estimated using Bayesian methods. The best -fit model, GTR+I+G, was selected by hierarchical likelihood ratio tests in the program modeltest 3.7 [START_REF] Posada | MODELTEST: Testing the model of DNA substitution[END_REF], where GTR+I+G was favoured. We ran beast 2.3.1 [START_REF] Bouckaert | BEAST 2: A software platform for Bayesian evolutionary analysis[END_REF] to explore the best combinations of substitution, clock and population models, with 10 million generations. Twice the difference in Ln harmonic mean of the likelihood of each model combination (LnBF; [START_REF] Nicholls | Concordant phylogeography and cryptic speciation in two Western Palaearctic oak gall parasitoid species complexes[END_REF] was calculated with tracer 1.6 to assess the preferred models, based on the LnBF table [START_REF] Kass | Bayes factors[END_REF]. Using the GTR+I+G model, we ran beast under a strict clock, an uncorrelated exponential relaxed clock or lognormal relaxed clock with each population model (constant size, exponential growth, yule process, birth-death process). Bayes factors indicated that the combination of an uncorrelated exponential clock and the exponential population growth model was the best (with LnBFs from 0.47 to 198.19 and most >15).

Two independent runs of 30 million iterations were performed, with genealogies and model parameters sampled every 1,000 iterations. The chain convergence was checked based on ESSs (effective sample sizes) viewed in tracer 1.6. All ESSs for each parameter had to be larger than 200. The phylogenetic tree was summarized by treeannotator 1.8.1 and then viewed by figtree 1.3.1 [START_REF] Rambaut | FigTree v1.3.1[END_REF]. The dates of the most recent common ancestor were scaled by a mutation rate of 1 because no suitable fossil record was available to calibrate node ages. We used the published 34.8 (46.5-24.2) Ma age for the crown group of the genus Valisia that had been estimated previously for a very similar set of Valisia species that used numerous genetic markers and a comprehensive set of species within Agaonidae [START_REF] Cruaud | An extreme case of plant-insect codiversification: Figs and figpollinating wasps[END_REF]. We used the phylogenetic tree to detect candidate species characterized by low within -clade genetic distances and much larger between -clade genetic distances.

We calculated Kimura 2 -parameter (K2P) distances within and between clades for COI haplotypes using mega 6.06. Tajima's D, Fu and Li's D, and Fu and Li's F (Fu & Li, 1993) were used to detect signatures of population expansion or selection using dnasp 5.0 [START_REF] Librado | DnaSP v5: A software for comprehensive analysis of DNA polymorphism data[END_REF]. We then explored the relationships of COI haplotypes within each clade using phylogenetic networks built separately for each clade with tcs 1.21 [START_REF] Clement | TCS: A computer program to estimate gene genealogies[END_REF] using the 95% statistical parsimony criterion as a connection limit, with loops in the network resolved following Crandall and Templeton's methods (1993). Variation in COI sequences among populations within clades and within populations was partitioned using hierarchical analysis of molecular variance (AMOVA) implemented in arlequin 3.0 [START_REF] Excoffier | Arlequin ver. 3.0: An integrated software package for population genetics data analysis[END_REF]. Significance tests were based on 10,000 permutations.

We used jMOTU as a complementary method to detect molecular operational taxonomic units (MOTUs) using a range of threshold differences. All the COI sequences used in the phylogenetic analysis were analysed using jMOTU at cut -offs from 0 to 21 bases (Jones, Ghoorah, & Blaxter, 2011). The aggregation parameter was 97%, and the input minimum sequence percentage was set to 95%.

As ITS2 evolves more slowly than COI, we applied a simplified analysis of the data. We calculated K2P distances within and between clades as for COI haplotypes using mega 6.06, and we used a maximum -likelihood tree to reconstruct the phylogenetic relationships based on all ITS haplotypes. The ML tree was reconstructed using mega 6.06 [START_REF] Tamura | MEGA6: Molecular evolutionary genetics analysis version 6.0[END_REF], and node supports were assessed based on 2,000 bootstrap replicates.

| -Pollinating fig wasp microsatellite data analysis

Classical indices of genetic diversity were estimated using genalex 6.1 [START_REF] Peakall | GENALEX 6: Genetic analysis in Excel. Population genetic software for teaching and research[END_REF]. In order to represent the global data, we performed a factorial correspondence analysis following [START_REF] Nenadi | Correspondence analysis in R, with two -and threedimensional graphics: The ca package[END_REF] as implemented in Genetix (Belkhir, Borsa, Chikhi, Raufaste, & Bonhomme, 1996[START_REF] Edgar | MUSCLE: Multiple sequence alignment with high accuracy and high throughput[END_REF]. Because multiple locus population genetics data should not be analysed using phylogenetic trees, we constructed an unweighted neighbour -joining tree of multilocus microsatellite genotypes using darwin v6 [START_REF] Perrier | DARwin software[END_REF], thus grouping genotypes according to the numbers of shared alleles, without taking into account gene evolution. Bootstrap number was set at 2,000. Finally, we used Bayesian clustering to assign multilocus microsatellite genotypes to clusters using structure 2.2 [START_REF] Pritchard | Inference of population structure using multi -locus genotype data[END_REF]. The admixture ancestry and correlated allele frequencies model were used with five independent runs each of 500,000 MCMC iterations and 500,000 burn -in steps. We ran structure with varying K values (the number of clusters) from 2 to 30 (the maximum number of populations available for microsatellite analysis) to explore how different values of K resulted in the assignment of species and populations to different clusters. The resulting assignment patterns were explored to detect co -occurrence of different clades within a sampling location. Results obtained from the factorial correspondence analysis, neighbour -joining tree and Bayesian clustering were compared.

The presence of a relationship between (log) geographic distance and genetic differentiation FST/(1 -FST) values (isolation by distance) was evaluated with genalex 6.1 within the two clades detected in the above analysis and for which data from more than three locations were available, using a Mantel test with 10,000 permutations. When structuring was detected within clades, the distribution of pairwise FST values was explored to detect potential structuring into subclades.

Ficus hirta 2.4.1 Ficus hirta DNA extraction and amplification

We used previously published methodological details for cpDNA and nuclear microsatellite sequencing and genotyping [START_REF] Yu | Nuclear and chloroplast DNA phylogeography of Ficus hirta: Obligate pollination mutualism and constraints on range expansion in response to climate change[END_REF], and re -used all the previously published cytoplasmic data. The set of analysed microsatellite loci was reduced to avoid some amplification/reading problems, leaving seven microsatellite loci available: FS4 -11, Frub38, Frub398 and Frub436 [START_REF] Yu | Nuclear and chloroplast DNA phylogeography of Ficus hirta: Obligate pollination mutualism and constraints on range expansion in response to climate change[END_REF], and FH3, FH10 and FH47 [START_REF] Zheng | Development and characterization of microsatellite loci for Ficus hirta (Moraceae)[END_REF]. Using these, we analysed plant microsatellite data from two locations in southeast China and two locations in Hainan, and cytoplasmic data from 14 locations in southeast China and Hainan. This was in addition to the microsatellite data and cytoplasmic data from 17 locations elsewhere within the range of the plant, so we had a total of 21 locations with microsatellite data and 31 locations with cytoplasmic data.

Ficus hirta cytoplasmic DNA analysis

The sequences of the two cpDNA regions were concatenated (1,367 bp) and then aligned using MUSCLE in mega 6.06 with manual corrections. A matrix of combined sequences for trnL -trnF and trnS -trnG was constructed, and haplotypes were distinguished on the basis of nucleotide and insertion/deletion differences. The same procedures as for wasp COI data were used to determine polymorphism indices, detect signatures of population expansion and establish haplotype networks.

We used AMOVA to test the significance of cpDNA differentiation among populations (999 permutations) and tested for genetic isolation by distance by conducting a Mantel test of the correlation between FST/(1 -FST) and log -transformed geographic distance for all population pairs (999 permutations).

Ficus hirta microsatellite analyses

Parameters of genetic diversity were calculated as for pollinator microsatellite data. As with the pollinating wasp microsatellite data, we performed a factorial correspondence analysis, we build a neighbour -joining tree, and we used Bayesian clustering to assign multilocus microsatellite genotypes to clusters. The number of clusters was set to vary from 1 to 21.

The relationship between FST and geographic distance was plotted to visualize genetic differentiation. The significance of the correlation between FST/(1 -FST) and the log geographic distance was evaluated as with the insects.

RESULTS

Pollinating fig wasps 3.1.1 Pollinating fig wasp gene sequences

After alignment and exclusion of sites with gaps, 502 bp long COI sequences were obtained for phylogenetic analysis. A total of 162 haplotypes were obtained and 267 COI polymorphic sites were identified from the 330 sequenced fig wasps (Supporting Information Table S2). Haplotype sequences have been deposited in GenBank under Accession nos. KR873011-KR873047 and MF472722-MF472846. A 689 bp fragment of the ITS2 gene was amplified in 201 individuals. A total of 18 haplotypes were obtained. The sequences have been deposited in GenBank under Accession nos. MF467418-MF467426 and MF467428-MF467436.

The Bayesian COI tree separated the F. hirta pollinating fig wasps into nine clades (Figure 1). Gene sequence differences were weak within clades (Kimura -2 -parameter = 0.001-0.014) and high between clades (0.064-0.272) (Figure 1, Supporting Information Figure S1, Supporting Information Table S3). The between -clade distances are similar to those between currently recognized Agaonidae species [START_REF] Chen | Fig trees at the northern limit of their range: The distributions of cryptic pollinators indicate multiple glacial refugia[END_REF][START_REF] Yang | The incidence and pattern of copollinator diversification in dioecious and monoecious figs[END_REF].

The ITS2 maximum -likelihood tree separated the F. hirta pollinating wasps into seven clades (Supporting Information Figure S2, Supporting Information Table S4). ITS2 data were obtained for several individuals from each of 30 of the 31 locations for which COI data were obtained (Supporting Information Table S1). The ITS2 sequences grouped the genotypes according to sampling location (Supporting Information Table S5) in the same way as the S2), except that the ITS2 sequences provided less resolution and pooled the locations that gave COI clade 2 and COI clade 9 and pooled the locations that gave COI clade 6 and COI clade 7. Hence, cytoplasmic and nuclear sequence data provided consistent location -level assignment to clades, although the ITS2 sequences provided less resolution.

The program jMOTU groups sequences into clades (molecular taxonomic units) so that the minimum difference, expressed in numbers of bases, between two different MOTUs is larger than the chosen cut -off value. With increasing cut -off values, the number of inferred MOTUs drops rapidly and then stabilizes when withinspecies variation is included within a single MOTU, before dropping again when closely related species are pooled. The number of clades inferred from the COI sequences dropped to 33 MOTUs for a cut -off value of 9 and remained constant up to a cut -off value of 11 (Supporting Information Figure S3). For these cut -off values, all the clades recognized in the COI tree were separated into single MOTUs except for clade 7, which was separated into three MOTUs and clade 6 that was also separated into 3 MOTUs. Starting at the cut -off value of 17 bases (which is the 2.5% cut -off threshold value proposed as a general rule -of -thumb for species discrimination with COI; Jones et al., 2011) and up to a cut -off value of 21 bases, jMOTU gave 29 MOTUs. Each of the nine clades recognized from the phylogenetic tree (Figure 1) was recognized as a single MOTU, and the two sequences of fig wasps (AY842415 and HM802690) collected from Ficus langkokensis grouped into a single MOTU. Hence, the jMOTU analysis supports the proposition that each clade recognized in the COI tree is a separate operational taxonomic unit, that is a species.

Overall, the different analyses of sequence data support the conclusion that clades 1-9 can be considered as distinct species, and this conclusion is applied from here on (as sp1-sp9).

The pollinators collected in Singapore (sp5) formed part of a distinct clade that had the largest gene sequence differences relative to the other species pollinating F. hirta (for cytoplasmic COI, Figure 1, Supporting Information Table S3; for nuclear ITS2, Supporting Information Figure S2, Supporting Information Table S4). This clade included a COI GenBank sequence (AY642456) of a fig wasp collected from F. androchaete Corner, a Ficus species endemic to Borneo, where F. hirta is absent [START_REF] Berg | Moraceae (Ficus)[END_REF]. The other major clade comprised sp1-4 and sp6-9 reared from F. hirta, plus the GenBank sequence for Valisia esquiroliana Chen & Chou. The dated tree obtained with beast (Supporting Information Figure S1) suggests that separation of sp5 from the remaining eight species associated with F. hirta began about 16 (11-22) Ma and the two tightly knit subgroups of species (sp2 + sp9) and (sp4 + sp6 + sp7) diverged from each other about 8 Ma (Supporting Information Figure S1). These two species groups were also recovered in the ITS2 phylogeny (Supporting Information Figure S2). The estimated dates of divergence among the taxa within these subgroups that we are recognizing as distinct species varied between 4.6 (sp2 + sp9 and sp4 + sp6 + sp7) and 6.0 (sp1 + sp3) Ma.

A striking feature was the limited overlap in the distribution of species. We only observed examples of a pair of species occurring together in two locations, location CS (with sp2 and a low frequency of sp3) in north Thailand and Wu in northeast Thailand (with sp4 together with small numbers of sp7; Figure 2a). In addition, GenBank Accession no. FJ619191 belonged to sp3 and was recorded from XI in SW China, an area where we only collected sp2. Where species were collected from several sites, their distributions were generally geographically coherent, but sp7 has a disjoint distribution. It is present in east Thailand and in Java, but absent from peninsular Malaysia where its host plant is also absent and is apparently replaced in peninsular Thailand by sp8 (Figure 2a). AMOVA revealed considerable haplotype variation within species, as well as within and among locations (Supporting Information Table S6). The exception was sp5, where haplotype diversity was strikingly low, with only two almost identical haplotypes represented in the 17 individuals (Supporting Information Figure S4). In the two more extensively sampled species, we observed negative deviations from neutral assumptions in both sp1 (Tajima's D value = -1.51, p = 0.055; Fu and Li's D = -2.75, p < 0.05; and Fu and Li's F = -2.69, p < 0.05) and sp2 (Tajima's D = -1.88, p < 0.05; Fu and Li's D = -3.49, p < 0.01; and Fu and Li's F = -3.40, p < 0.01). Deviations were not significant in the other species. The COI haplotype networks (Supporting Information Figure S4) revealed a subdivision in sp1 into a southeast China group of haplotypes, a Hainan Island group of haplotypes and a Vietnamese group of haplotypes, with some haplotype exchange between southeast China and Hainan Island.

Pollinating fig wasps microsatellite data

Diversity indices for microsatellite markers are given in Supporting Information Table S7. Pairwise genetic differentiation between location values (FST) is given in Supporting Information Table S8.

The first eight components of the discriminant analysis separated the nine species (Supporting Information Figure S5). The analysis also shows differentiation within sp1 according to geographic origin (China, Hainan and North Vietnam), as suggested by the COI analysis (Supporting Information Figure S4). While ITS2 sequences did not allow separation of sp2 from sp9 and sp6 from sp7, the discriminant analysis of microsatellite data separates them without ambiguity on component 2. Hence, the discriminant analysis confirms the presence of nine nuclear genetic clusters corresponding to the species detected using the cytoplasmic COI sequences.

In the neighbour -joining tree (Figure 3a), the different species are again separated, except for sp5, which seems to be close to sp1 from Vietnam, a result that is not compatible with the nuclear and cytoplasmic sequence data. The subdivision of sp1 into three geographic entities is supported. The close proximity between sp4 and sp6 is also supported, as is the slight separation within sp7 depending on geographic origin.

The bar plots of assignments of individuals to clusters using structure with variable numbers of clusters are presented in Supporting Information Figure S6a. Irrespective of the different values of K, assignments to clusters are generally strong, except for the south China locations for sp1 for which mixed assignments suggest the presence of some variation among locations. From K = 6 to 22, the populations of sp1 from Hainan island and Vietnam are each assigned to their own cluster. For K = 11, all the species are separated into different clusters except for sp4 and sp6 (Figure 4a). For K = 20 and above, sp4 and sp6 are separated into distinct clusters. For K = 22 and above, the individuals of sp7 collected in Thailand were perfectly separated from those collected in Java. Consistently, the separation of sp3 to sp9 is somewhat unstable across the different values of K. As with the COI data, and despite different individuals being used for COI and microsatellite analyses, the microsatellite data suggested assignment of one individual to sp3 in location CS and another individual to sp7 in location Wu.

Genetic isolation by distance could only be explored for sp1 and sp2 as the other species were only sampled in 1-3 locations. There was significant isolation by distance, as shown by the relationship of FST/(1 -FST) with distance (Mantel test) in sp1 (R2 = 0.13, p = 0.008), but not in sp2 (R2 = 0.023, p = 0.459). However, the signature of isolation by distance in sp1 was generated by the differentiation of sp1 into the three genetically differentiated geographic entities already detected in the COI, in the discriminant and in the structure analyses (located in southeast China, Hainan and Vietnam, Figure 5), despite obtaining a single ITS2 haplotype for sp1 (with 33, 11 and 9 individuals genotyped, respectively). There was no genetic isolation by distance within any of these three geographic entities.

Ficus hirta 3.2.1 Ficus hirta cpDNA analysis

We detected 24 haplotypes (H1-H24). Sequences of the trnL -trnF and trnS -trnG regions have been deposited in GenBank under Accession nos. GQ452019-GQ452032 and MF467405-MF467416. We did not detect any signature of selection or of population expansion on cpDNA sequences. The cpDNA network and the distribution of these haplotypes among populations are presented in Figure 2b and Supporting Information Table S9. The haplotype network is centred on haplotype H8, and no haplotype differed from H8 by more than six mutations (Figure 2b). Haplotype H8 was observed throughout the range of F. hirta. Haplotype diversity was high in the north -central part of our sampling zone, where sampling density was highest. AMOVA confirmed significant differentiation among populations (FST = 0.799; p < 0.001). The colours correspond to different branches of the network, and their geographic distribution suggests spatial structuring. The regression of population pairwise FST on the natural logarithm of geographic distance was also significant (R2 = 0.0472, p = 0.003, Mantel test), confirming that there was spatial genetic structure among populations of the plant.

Ficus hirta microsatellite analysis

Genetic diversity parameters are given in Supporting Information Table S10. Pairwise genetic differentiation between location values (FST) is given in Supporting Information Table S11.

The first four components of the discriminant analysis organized the data along a north-south axis, except for the genotypes from Singapore, which were placed close to the northern genotypes (Supporting Information Figure S7). The discriminant analysis suggests much more continuous genetic variation in F. hirta than in its pollinators (compare Supporting Information Figures S5 andS7).

As with the discriminant analysis, the neighbour -joining tree separates the data along a north-south axis (Figure 3b). The Singapore samples and the Trang samples each form highly homogeneous clusters. As in the discriminants analysis, the Trang samples are placed with southern locations while the Singapore samples are placed with northern locations.

The bar plots of assignments of individuals to clusters of varying total number using structure are presented in Supporting Information Figure S6b. For K = 2, genotypes from all northern locations plus Singapore were mainly assigned to cluster one and genotypes from southern locations were mainly assigned to cluster two (Supporting Information Figure S6b). The geographic distribution of assignments (Figure 2a) showed: (a) progressive geographic genetic variation of F. hirta along a north-south axis; (b) some geographic east-west geographic genetic structure; and (c) a placement of the Singapore samples that does not correspond to its geographic location. For K = 4, genotypes from the northern locations had mixed assignments to two clusters, supporting an east-west spatial genetic structure, genotypes from Singapore still clustered with genotypes from the north, and genotypes from Trang formed a distinct group (Figure 4b).

We also explored spatial genetic structure by plotting genetic differentiation between locations against geographic distance (Figure 6). Comparisons between locations showed a general pattern of genetic isolation by distance (R2 = 0.204, permutation test, p = <10 -3). Singapore was an exception, as it was most similar to the geographically distant northern locations.

Taken together, the factorial analysis, the structure assignments to clusters and the neighbour -joining tree all suggest that the genetic diversity of F. hirta is organized according to a pattern of genetic

DISCUSSION

Species status of the pollinators

Although different pollinating wasp individuals were genotyped using COI, ITS2 and microsatellites, the parapatric distribution of the wasp clades allows direct comparison of results obtained with the different markers. Cytoplasmic and nuclear markers subdivide the pollinating wasps into the same nine clades. The COI divergence between clades (ranging from 5.4% to 28%, with most values above 10%) is larger than the divergence previously reported between sister species of fig pollinating wasps that share the same host (2.4%-7.4%; [START_REF] Yang | The incidence and pattern of copollinator diversification in dioecious and monoecious figs[END_REF]. Furthermore, the nine clades can be separated morphologically [START_REF] Wiebes | Agaonidae (Hymenoptera Chalcidoidea) and Ficus (Moraceae): Fig wasps and their figs, xi (Blastophaga) s.l[END_REF]J.Y. Rasplus, com. pers.) and we therefore conclude that F. hirta is pollinated by nine species of fig wasps in our study zone, which covered most but not all of the range of the plant.

Sp5 is not closely related to the other pollinators of F. hirta. It was recorded from Singapore, where our data confirm an earlier proposition that F. hirta is not native [START_REF] Berg | Moraceae (Ficus)[END_REF]. We suggest that sp5 is the regular pollinator of another Ficus species that has colonized F. hirta in Singapore. [START_REF] Cook | Speciation in fig wasps[END_REF][START_REF] Wang | Pollinator sharing and gene flow among closely related sympatric dioecious fig taxa[END_REF]. Sequences from GenBank suggest that Valisia esquiroliana Chen & Chou is nested within the group of Valisia species associated with F. hirta. This fig wasp was described from Ficus triloba Buch. -Ham. ex Voigt, a species closely related to F. hirta [START_REF] Berg | Precursory taxonomic studies on Ficus (Moraceae) for the Flora of Thailand[END_REF][START_REF] Berg | Moraceae (Ficus)[END_REF]. Genetic analyses incorporating these two Ficus species will be required to clarify their evolutionary relationships.

Genetic differentiation was detected within the distribution of pollinator sp1. Its populations are divided into three genetic groups located in continental China, Hainan Island and Vietnam. Such differentiation could be a first step towards speciation. Experimental data are needed to determine whether these genetic groups are interfertile. Reproductive isolation could arise rapidly in fig pollinating wasps as they display a high prevalence of Wolbachia that could cause cytoplasmic incompatibilities [START_REF] Haine | Convergent incidences of Wolbachia infection in fig wasp communities from two continents[END_REF] and they display systematic assortative mating due to mating in their natal figs before dispersal [START_REF] Anstett | Figs and fig pollinators: Evolutionary conflicts in a coevolved mutualism[END_REF].

The biogeography of Ficus and their pollinators

Ficus hirta displays a signature of genetic isolation by distance with both nuclear microsatellites and the cytoplasmic haplotypes, with no marked subdivision into distinct gene pools. It has been proposed that speciation in nursery pollinators and their host plants may be decoupled [START_REF] Hembry | Diversification and coevolution in brood pollination mutualisms: Windows into the role of biotic interactions in generating biological diversity[END_REF] and our results provide an example of this phenomenon.

The disjunct distribution of sp7, with a 2,000 -km gap and with three species of pollinating wasps (sp5, sp8 and sp9) located in between, suggests that pollinators can successfully colonize distant host plants and establish populations there. A similar example is known from Australia, where one pollinator of F. rubiginosa is present in the cool south and at a higher altitude in the warmer north, resulting in a 1,500 -km distribution gap [START_REF] Sutton | One step ahead: A parasitoid disperses farther and forms a wider geographic population than its fig wasp host[END_REF]. This suggests that pollinating fig wasp distributions are plastic and do not necessarily reflect a classical biogeographic history. Indeed, the biogeographic history of F. racemosa populations differs from that of its pollinator species, despite their identical separation into large gene pools [START_REF] Bain | Geographic structuring into vicariant species -pairs in a wide -ranging, high -dispersal plant -insect mutualism: The case of Ficus racemosa and its pollinating wasps[END_REF]. Similarly, in F. septica, an ecologically divergent pollinating fig wasp has colonized a 2,000km range previously only occupied by three parapatric pollinating wasp species [START_REF] Rodriguez | Diversification and spatial structuring in the mutualism between Ficus septica and its pollinating wasps in insular South East Asia[END_REF]. Thus, the distributions of fig pollinating wasp species seem to be dynamic and reveal cases of successful establishment of long -distance migrants. The longrange dispersal abilities of at least some fig wasps should facilitate this [START_REF] Ahmed | Windborne insects mediate directional pollen transfer between desert fig trees 160 kilometers apart[END_REF].

The mismatch distributions of COI sequences in the two bestsampled species (sp1 and sp2) provide signatures of either selective sweeps or population expansion. This is similar to what has been found in the pollinators of F. pumila [START_REF] Chen | Fig trees at the northern limit of their range: The distributions of cryptic pollinators indicate multiple glacial refugia[END_REF] and in a pollinator of F. septica [START_REF] Lin | Fig trees at the northern limit of their range: The distributions of cryptic pollinators indicate multiple glacial refugia[END_REF]. It has been proposed that these signatures are due to recurrent regional extinctions of pollinating wasps during extreme climatic events, but not of their host plants, followed by population expansions [START_REF] Tian | Lack of genetic isolation by distance, similar genetic structuring but different demographic histories in a fig -pollinating wasp mutualism[END_REF]. Contemporary examples of this phenomenon have been documented in Florida (hurricane), southern France (frost) and the north of Borneo Island (El Nino -associated draught; [START_REF] Tian | Lack of genetic isolation by distance, similar genetic structuring but different demographic histories in a fig -pollinating wasp mutualism[END_REF]. Climatic episodes such as these make fig wasp species ranges intrinsically dynamic.

The establishment of new disjunct populations by long -distance migrants depends on a combination of suitable environmental and biological conditions in their new locations. Founder populations will also usually be small, which can lead to Allee effects. These result from (a) genetic inbreeding and loss of heterozygosity; (b) demographic stochasticity (including sex -ratio fluctuations); and (c) a reduction in cooperative interactions when individuals are at low density [START_REF] Courchamp | Inverse density dependence and the Allee effect[END_REF] [START_REF] Molbo | Inbreeding and population structure in two pairs of cryptic fig wasp species[END_REF], and consequently, we expect founder populations of fig wasps to be little affected by inbreeding. Furthermore, although offspring sex ratios are heavily female -biased, foundress females typically begin by laying male eggs in a clutch, so female offspring rarely have difficulty finding mates [START_REF] Raja | The mechanism of sex ratio adjustment in a pollinating fig wasp[END_REF]. Finally, foundress females compete within figs for oviposition sites, so no cooperation between them is expected [START_REF] Herre | Coevolution and reproductive characteristics in 12 species of New World figs and their pollinator wasps[END_REF]. Consequently, pollinating fig wasps are probably not sensitive to Allee effects and they should be much less dependent on progressive geographic range expansion than many other organisms. Within this context, the parapatric distributions of many pollinating fig wasps suggest an important role for interspecific competition in shaping their distributions, as has been proposed for nonpollinating fig wasps [START_REF] Darwell | Cryptic diversity in a fig wasp community-Morphologically differentiated species are sympatric but cryptic species exhibit competitive exclusion[END_REF].

Diversification in plant-insect mutualisms

The total of nine pollinator species found on F. hirta represents the highest ever -reported number of pollinator species for a Ficus species. Almost as high pollinator diversity has been documented for some other Ficus species (five for F. rubiginosa, [START_REF] Darwell | Molecular species delimitation of a symbiotic fig -pollinating wasp species complex reveals extreme deviation from reciprocal partner specificity[END_REF]; five for F. racemosa, [START_REF] Bain | Geographic structuring into vicariant species -pairs in a wide -ranging, high -dispersal plant -insect mutualism: The case of Ficus racemosa and its pollinating wasps[END_REF]; four for F. septica, [START_REF] Rodriguez | Diversification and spatial structuring in the mutualism between Ficus septica and its pollinating wasps in insular South East Asia[END_REF]. This is strikingly different from the comprehensively known Yucca moth-Yucca association, where from a total of 24 pollinating moth species [START_REF] Pellmyr | Phylogeny of the pollinating yucca moths, with revision of Mexican species (Tegeticula and Parategeticula; Lepidoptera, Prodoxidae)[END_REF] there are only five instances of two copollinators sharing on a host [START_REF] Smith | Host specificity and reproductive success of yucca moths (Tegeticula spp. Lepidoptera: Prodoxidae) mirror patterns of gene flow between host plant varieties of the Joshua tree (Yucca brevifolia: Agavaceae)[END_REF]. Available data on the less comprehensively investigated Epicephala-Phyllantheae association suggest an intermediate situation with up to four pollinator species associated with one host [START_REF] Li | Four new species of Epicephala Meyrick, 1880 (Lepidoptera, Gracillariidae) associated with two species of Glochidion (Phyllanthaceae) from Hainan Island in China[END_REF]. The comparatively low pollinator diversity observed in Yucca moths could be related to several factors. First, Yuccas have much more restricted distribution ranges than Ficus species, with only the distribution ranges of the two most northern Yucca species reaching 2,000 km [START_REF] Althoff | Geographic isolation trumps coevolution as a driver of yucca and yucca moth diversification[END_REF]. Second, Yucca moths can exhibit prolonged diapause, a factor that will slow down the dynamics of population divergence [START_REF] Powel | Interrelationships of yuccas and yucca moths[END_REF]. Third, assortative mating is not systematic, unlike in fig wasps [START_REF] Powel | Interrelationships of yuccas and yucca moths[END_REF].

In addition to their contrasting rates of differentiation, the population genetics of F. hirta and its pollinators are strikingly different. Whereas F. hirta exhibits clinal genetic variation with a signature of genetic isolation by distance, its pollinator species sp1 and sp2 lack any such signature. The same pattern has been observed for F. pumila L. and one of its pollinating wasps in southeast China [START_REF] Liu | Movements of genes between populations: Are pollinators more effective at transferring their own or plant genetic markers[END_REF]. Hence, although pollinating fig wasps diversify faster than their host figs, they also display signatures of stronger gene flow. This apparent contradiction suggests that factors which facilitate speciation, such as a short generation time, infection by Wolbachia and assortative mating, can compensate for extensive gene flow and be decisive factors facilitating pollinating fig wasp speciation. Some nonpollinating fig wasps (NPFW, belonging to several families of Chalcidoidea) may disperse as widely as pollinator fig wasps [START_REF] Kjellberg | Tracking the elusive history of diversification in plantherbivorous insect-parasitoid food webs: Insights from figs and fig wasps[END_REF][START_REF] Sutton | One step ahead: A parasitoid disperses farther and forms a wider geographic population than its fig wasp host[END_REF]. NPFW have a largely similar biology to the pollinators, but exhibit varying intensities of local mate competition due to varying patterns of oviposition and mating sites. Obligatory assortative mating among NPFW ranges from minimal (among species that mate outside natal figs), through intermediate (NPFW with offspring dispersed across numerous figs that mate in the cavity of their natal figs) to being highly similar to that of pollinating wasps (with offspring aggregated in a single or small number of figs and mating inside these natal figs; [START_REF] Cook | Mutualists with attitude: Coevolving fig wasps and figs[END_REF]. If assortative mating indeed facilitates speciation, then we would expect a correlation between NPFW reproductive behaviour and speciation rates in different NPFW clades.

The mismatch between genetic structuring of host plant and fig wasps shows that, as in Yucca moths [START_REF] Godsoe | Absence of population -level phenotype matching in an obligate pollination mutualism[END_REF] and Epicephala moths [START_REF] Hembry | Does biological intimacy shape ecological network structure -A test using a brood pollination mutualism on continental and oceanic islands[END_REF], fine population -level phenotypic matching between mutualists is not required to allow population persistence. Nevertheless, on Hainan Island, the populations of F. hirta and its pollinator present the same signature of marked genetic differentiation from continental populations. Hence, the mismatch in population genetic structure is not consistent over the whole range of the association. Such differences in genetic costructuring are also encountered among associations in other species -specific plant-insect interactions. For instance, in a specialized association between two ants and an ant plant, with populations arranged as stepping stones, the ants and the plant presented very similar genetic signatures of population geographic expansion [START_REF] Léotard | Range expansion drives dispersal evolution in an Equatorial three -species symbiosis[END_REF]. In a similar association showing a history of restriction to refugia followed by expansion leading to more continuous populations, the plants and insects had the same pattern of genetic geographic structuring, but different histories [START_REF] Blatrix | Comparative analysis of spatial genetic structure in an ant-plant symbiosis reveals a tension zone and highlights speciation processes in tropical Africa[END_REF]. The latter study also showed that the trees had evolved stronger reproductive isolation than their ant symbionts. Hence, available case studies suggest that patterns of diversification in plant-insect mutualistic interactions are not uniform among representatives of the same type of interaction. As more case studies become available, exploring how variation is explained by the biology of individual species is the next challenge.

Relationships between tree characteristics and diversity of pollinators

Pollinator diversity has now been explored throughout large parts of the distribution ranges of four widely distributed Asian and Australasian Ficus species (F. hirta, F. septica, F. racemosa and F. rubiginosa; [START_REF] Bain | Geographic structuring into vicariant species -pairs in a wide -ranging, high -dispersal plant -insect mutualism: The case of Ficus racemosa and its pollinating wasps[END_REF][START_REF] Darwell | Molecular species delimitation of a symbiotic fig -pollinating wasp species complex reveals extreme deviation from reciprocal partner specificity[END_REF][START_REF] Rodriguez | Diversification and spatial structuring in the mutualism between Ficus septica and its pollinating wasps in insular South East Asia[END_REF]. Their pollinator fig wasps each form old monophyletic species complexes (dating back >10 Ma) associated with a single host [START_REF] Bain | Geographic structuring into vicariant species -pairs in a wide -ranging, high -dispersal plant -insect mutualism: The case of Ficus racemosa and its pollinating wasps[END_REF][START_REF] Cruaud | An extreme case of plant-insect codiversification: Figs and figpollinating wasps[END_REF], but the species groups display varying geographic scales of pollinator species diversification. The host figs also display varying scales of genetic differentiation. Further studies are required to investigate whether plant traits such as crop size, crop synchrony and growth form may combine to select for different pollinator dispersal behaviour and may result in pollinator and fig genetic diversity being expressed at different geographic scales.

  Fig trees are pollinated by female fig wasps that enter Ficus inflorescences to lay their eggs. Because fig wasps carry pollen from their natal plants, they

  hirta and its pollinating fig wasps across most of their continental range. We addressed the following questions: (a) To what extent is the population structure of the plant co -incident with that of its pollinator(s)-have they diversified at the same spatial scales -(b) Even in the absence of major geographic barriers, is this widespread host tree pollinated by a diverse assemblage of fig wasps within its range -And (c) if so, what are the ecological consequences-do individual plants have the chance to receive pollination services from more than one fig wasp species -This is the first study comparing spatial genetic structuring in a widespread fig species and its pollinating wasps that encompasses most of their ranges. It is also the first broad -scale study of joint genetic structuring between pollinator fig wasps and a fig tree species displaying genetic isolation by distance. We discuss the results in terms of the factors that may have resulted in the patterns of diversification we observed and what they suggest about modes of speciation.

  has unique protogynous inflorescences called figs. These are hollow and lined internally with numerous tiny male and female florets. Pollination can only be achieved when pollen is carried by female fig wasps from their natal fig into receptive figs through a narrow aperture called the ostiole. Like about half of all Ficus species, F. hirta is functionally dioecious, with figs that differ in floral structure on "male" and "female" trees. Figs of the former support the development of the fig wasp offspring that when mature can transport pollen to receptive figs.

  At each location, 10-30 figs containing mature fig wasp offspring were also removed from the plants and placed individually in fine -mesh bags, where the fig wasps were allowed to emerge. The pollinators were then separated from other species of fig wasps and preserved in 95% ethanol that was stored at -20°C until DNA extraction. A single female fig wasp per fig was used for genetic analyses.

FIGURE 1 -

 1 FIGURE 1 -COI Bayesian phylogenetic tree of the Valisia fig wasp pollinators associated with Ficus hirta, including all GenBank sequences of pollinating sp5 wasps reared from Ficus subg. Ficus sect. Eriosycea Miq. and some sequences of Ceratosolen and Kradibia used as outgroups. Posterior probabilities of the nodes are indicated as percentages

FIGURE 3 -

 3 FIGURE 3 -Neighbour -joining microsatellite genotype networks for (a) the pollinating wasps and (b) Ficus of F. hirta. Note the disjoint distribution of sp7, the allopatric distribution of the closely related (according to COI data) sp2 and sp9 and the parapatric distribution of the closely related (according to COI data) sp4, sp6 and sp7. For cpDNA, colours are used to show the localized distribution of haplotype H1 and of the different branches of the haplotype network hirta.

FIGURE 4 -

 4 FIGURE 4 -Bar plots of membership probabilities of individuals to the different clusters (as vertical bars) from the structure analysis. (a) Genetic structure of the pollinating wasps obtained at K = 11; (b) genetic structure of the plant host, Ficus hirta, obtained at K = 4. For the pollinator, sp1 genotypes from continental China present mixed assignments to two clusters suggesting genetic isolation by distance, while genotypes from Hainan and genotypes from Vietnam are each assigned to a particular cluster. For F. hirta, there is no clear global pattern of

FIGURE 5 -

 5 FIGURE 5 -Spatial genetic structuring of the pollinator fig wasp sp1 based on nuclear microsatellite data. Genetic differentiation between locations is plotted according to geographic distance. The colours indicate the geographic origin of the locations being compared. Dark blue: comparison between two locations from southeast China; red: comparison between a location in southeast China and one in Hainan Island; green: comparison between a location from southeast China and one from Vietnam; violet: comparison between a location from Hainan Island and one from Vietnam; light blue, comparison between two locations from Vietnam; orange: comparison between the two Hainan Island locations. The species is structured into three genetic groups: southeast China, Hainan and Vietnam, the same structure as revealed by the cytoplasmic data (Supporting Information Figure S2)

FIGURE 6 -

 6 FIGURE 6 -Pairwise genetic differentiation between locations in Ficus hirta, as a function of distance, based on microsatellite data. Red points: comparisons involving the Singapore location.

  Fig trees planted outside their native range are known to similarly be capable of supporting successful development by fig wasps that usually breed in the figs of another Ficus species

  . As a fig is often colonized by a single fig wasp and given that their offspring mate in their natal fig before dispersal, inbreeding among fig wasps is routine

ACKNOWLEDGEMENTS

We thank Prof. Xia Nianhe, Peng Yanqiong, Vu Quang Nam (Vietnam National Foundation for Science and Technology Development under grant numbers 106.11 -2012.82), Tien Tran, Nong Van Duy, Zhang Xianchun (collecting permission number in Indonesia: No. 1799IPH.1.02/KS.01/XI/2014), Harry Wiriadinata, Arief Hidayat, Wei Ran, Chen Wenli, Mashhor Mansor, Nora Sikin, Serena Lee (provided Singapore samples), Khin Me Aung, Chuck Cannon, Thitima Tharawoot and Liang Dan for sample collection and help in the field. We thank Mei Qiming, Liu Qing, Kong Hanghui and Tu Tieyao for their help with cluster computers. We also thank Liu Yunxiao for her help with figures. We thank Simon Segar, Tatiana Giraud and the anonymous reviewers for helpful comments on the manuscript. This work was supported by the National Natural Science Foundation of China (31670395; 31630008; 31370409), the Guangzhou Science and Technology Plan Project (201707010351), the External Cooperation Program of Bureau of International Co -operation Chinese Academy of Sciences (No. GJHZ201321) and the CAS President's International Fellowship Initiative (2017VBA0039).

Funding information

National Natural Science Foundation of China, Grant/Award Number: 31670395, 31630008 and 31370409; Guangzhou Science and Technology Plan Project, Grant/Award Number: 201707010351; External Cooperation Program of Bureau of International Co-operation Chinese Academy of Sciences, Grant/Award Number: GJHZ201321; CAS President's International Fellowship Initiative, Grant/Award Number: 2018VBA0040; Vietnam National Foundation for Science and Technology Development, Grant/Award Number: 106.11 -2012.and cowrote the manuscript.

CONFLICT OF INTEREST

The authors declared that they have no competing interests.