Pollinator and host sharing lead to hybridization and introgression in Panamanian free-standing figs, but not in their pollinator wasps.

Obligate pollination mutualisms, in which plant and pollinator lineages depend on each other for reproduction, often exhibit high levels of species specificity. However, cases in which two or more pollinator species share a single host species (host sharing), or two or more host species share a single pollinator species (pollinator sharing), are known to occur in current ecological time. Further, evidence for host switching in evolutionary time is increasingly being recognized in these systems. The degree to which departures from strict specificity differentially affect the potential for hybridization and introgression in the associated host or pollinator is unclear. We addressed this question using genome-wide sequence data from five sympatric Panamanian free-standing fig species (Ficus subgenus Pharmacosycea, section Pharmacosycea) and their six associated fig-pollinator wasp species (Tetrapus). Two of the five fig species, F. glabrata and F. maxima, were found to regularly share pollinators. In these species, ongoing hybridization was demonstrated by the detection of several first-generation (F1) hybrid individuals, and historical introgression was indicated by phylogenetic network analysis. By contrast, although two of the pollinator species regularly share hosts, all six species were genetically distinct and deeply divergent, with no evidence for either hybridization or introgression. This pattern is consistent with results from other obligate pollination mutualisms, suggesting that, in contrast to their host plants, pollinators appear to be reproductively isolated, even when different species of pollinators mate in shared hosts.

Extraordinarily precise nematode sex ratios: adaptive responses to vanishingly rare mating opportunities.

Sex ratio theory predicts both mean sex ratio and variance under a range of population structures. Here, we compare two genera of phoretic nematodes (Parasitodiplogaster and Ficophagus spp.) associated with 12 fig pollinating wasp species in Panama. The host wasps exhibit classic local mate competition: only inseminated females disperse from natal figs, and their offspring form mating pools that consist of scores of the adult offspring contributed by one or a few foundress mothers. By contrast, in both nematode genera, only sexually undifferentiated juveniles disperse and their mating pools routinely consist of 10 or fewer adults. Across all mating pool sizes, the sex ratios observed in both nematode genera are consistently female-biased (approx. 0.34 males), but markedly less female-biased than is often observed in the host wasps (approx. 0.10 males). In further contrast with their hosts, variances in nematode sex ratios are also consistently precise (significantly less than binomial). The constraints associated with predictably small mating pools within highly subdivided populations appear to select for precise sex ratios that contribute both to the reproductive success of individual nematodes, and to the evolutionary persistence of nematode species. We suggest that some form of environmental sex determination underlies these precise sex ratios.

Genome-wide sequence data show no evidence of hybridization and introgression among pollinator wasps associated with a community of Panamanian strangler figs.

The specificity of pollinator host choice influences opportunities for reproductive isolation in their host plants. Similarly, host plants can influence opportunities for reproductive isolation in their pollinators. For example, in the fig and fig wasp mutualism, offspring of fig pollinator wasps mate inside the inflorescence that the mothers pollinate. Although often host specific, multiple fig pollinator species are sometimes associated with the same fig species, potentially enabling hybridization between wasp species. Here, we study the 19 pollinator species (Pegoscapus spp.) associated with an entire community of 16 Panamanian strangler fig species (Ficus subgenus Urostigma, section Americanae) to determine whether the previously documented history of pollinator host switching and current host sharing predicts genetic admixture among the pollinator species, as has been observed in their host figs. Specifically, we use genome-wide ultraconserved element (UCE) loci to estimate phylogenetic relationships and test for hybridization and introgression among the pollinator species. In all cases, we recover well-delimited pollinator species that contain high interspecific divergence. Even among pairs of pollinator species that currently reproduce within syconia of shared host fig species, we found no evidence of hybridization or introgression. This is in contrast to their host figs, where hybridization and introgression have been detected within this community, and more generally, within figs worldwide. Consistent with general patterns recovered among other obligate pollination mutualisms (e.g. yucca moths and yuccas), our results suggest that while hybridization and introgression are processes operating within the host plants, these processes are relatively unimportant within their associated insect pollinators.

Differential effects of nematode infection on pollinating and non-pollinating fig wasps: Can shared antagonism provide net benefits to a mutualism?

Species pairs that form mutualistic associations are also components of broader organismal community networks. These interaction networks have shaped the evolution of individual mutualisms through interspecific interactions ranging from secondarily mutualistic to intensely antagonistic. Our understanding of this complex context remains limited because characterizing the impacts of species interacting with focal mutualists is often difficult. How is the fitness of mutualists impacted by the co-occurring interactive network of community associates? We investigated this context using a model interaction network comprised of a fig and fig wasp mutualist, eight non-pollinating fig wasp (NPFW) antagonists/commensals and a nematode previously believed to be associated only with the pollinator wasp mutualist. Through repeated sampling and field observations, we characterized the ecological roles of these mutualist-associated organisms to identify key antagonists. We then investigated how potential nematode infection of NPFWs could impact wasp survival across key life stages and, in turn, inferred how this influences the fitness of the fig-pollinator mutualists. Unexpectedly, we found all Ficus petiolaris-associated NPFWs to be the targets for nematode infection, with infection levels sometimes exceeding that of pollinators. Experimental data collected for the most abundant NPFW species suggest that nematode infection significantly reduces their longevity. Further, comparisons of nematode loads for emerging and successfully arriving NPFWs suggest that infection severely limits their dispersal ability. Through these observations, we conclude that this infection could impact NPFWs more severely than either mutualistic partner, suggesting a novel role of density-dependent facultative mutualism between figs, pollinator wasps and the nematode. This antagonist-mediated suppression of other network antagonists may present an ecologically common mechanism through which antagonists can present net benefits for mutualists' fitness.

Inferring processes of coevolutionary diversification in a community of Panamanian strangler figs and associated pollinating wasps.

The fig and pollinator wasp obligate mutualism is diverse (∼750 described species), ecologically important, and ancient (∼80 Ma). Once thought to be an example of strict one-to-one cospeciation, current thinking suggests genera of pollinator wasps codiversify with corresponding sections of figs, but the degree to which cospeciation or other processes contribute to the association at finer scales is unclear. Here, we use genome-wide sequence data from a community of Panamanian strangler figs and associated wasp pollinators to estimate the relative contributions of four evolutionary processes generating cophylogenetic patterns in this mutualism: cospeciation, host switching, pollinator speciation, and pollinator extinction. Using a model-based approach adapted from the study of gene family evolution, our results demonstrate the importance of host switching of pollinator wasps at this fine phylogenetic and regional scale. Although we estimate a modest amount of cospeciation, simulations reveal the number of putative cospeciation events to be consistent with what would be expected by chance. Additionally, model selection tests identify host switching as a critical parameter for explaining cophylogenetic patterns in this system. Our study demonstrates a promising approach through which the history of evolutionary association between interacting lineages can be rigorously modeled and tested in a probabilistic phylogenetic framework.

A phylogenetic comparative method for evaluating trait coevolution across two phylogenies for sets of interacting species.

Evaluating trait correlations across species within a lineage via phylogenetic regression is fundamental to comparative evolutionary biology, but when traits of interest are derived from two sets of lineages that coevolve with one another, methods for evaluating such patterns in a dual-phylogenetic context remain underdeveloped. Here, we extend multivariate permutation-based phylogenetic regression to evaluate trait correlations in two sets of interacting species while accounting for their respective phylogenies. This extension is appropriate for both univariate and multivariate response data, and may use one or more independent variables, including environmental covariates. Imperfect correspondence between species in the interacting lineages can also be accommodated, such as when species in one lineage associate with multiple species in the other, or when there are unmatched taxa in one or both lineages. For both univariate and multivariate data, the method displays appropriate type I error, and statistical power increases with the strength of the trait covariation and the number of species in the phylogeny. These properties are retained even when there is not a 1:1 correspondence between lineages. Finally, we demonstrate the approach by evaluating the evolutionary correlation between traits in fig species and traits in their agaonid wasp pollinators. R computer code is provided.

Lack of genetic isolation by distance, similar genetic structuring but different demographic histories in a fig-pollinating wasp mutualism.

Historical abiotic factors such as climatic oscillations and extreme climatic events as well as biotic factors have shaped the structuring of species' genetic diversity. In obligate species-specific mutualisms, the biogeographic histories of the interacting species are tightly linked. This could be particularly true for nuclear genes in the Ficus-pollinating wasp mutualistic association as the insects disperse pollen from their natal tree. In this study, we compare spatial genetic structure of plant and pollinator for the Ficus hirta-Valisia javana association throughout southeast China including Hainan Island, for both nuclear and cytoplasmic markers. We show that dispersal of the insect leads to plant and insect presenting similar signatures of lack of genetic isolation by distance for nuclear genes on the continent over a distance of 1000 km. But we also show that the demographic histories of plant and insect are strikingly different. This is in agreement with extreme climatic events leading to transient regional extinctions of the insects, associated with local survival of the plants. We also observe evidence of genetic differentiation for both wasps and fig-tree between the continent and Hainan Island, although the Qiongzhou Strait is only on average 30 km wide, suggesting that geographic isolation by itself has not been sufficient to generate this differentiation. Hence, our results suggest that in highly dispersive mutualistic systems, isolation-by-dispersal limitation across a geographic barrier could be supplemented by isolation by adaptation, and maybe by coevolution, allowing further genetic divergence. In such systems, species may frequently be composed of a single population.

Transcriptome-facilitated development of SNPs for the Sonoran Desert rock fig, Ficus petiolaris (Moraceae).

PREMISE OF THE STUDY: Single-nucleotide polymorphism (SNP) primers were developed for a native North American desert fig, Ficus petiolaris (Moraceae), to provide markers for population genetic studies designed to quantify patterns of gene flow across a complex landscape. METHODS AND RESULTS: Transcriptome sequencing and bioinformatic protocols were implemented to discover SNPs in single-copy protein-coding genes. Multiplexes of 30 nuclear and 24 organellar (chloroplast and mitochondrial) SNPs were selected for primer development and genotyping on the Sequenom MASSArray System. Of these 54 loci, 49 reliably amplified across a panel of 96 F. petiolaris individuals. CONCLUSIONS: This study has provided SNP primers that can be applied in future studies investigating population genetics of F. petiolaris and its coevolution with associated pollinating and nonpollinating fig wasps.

Development and characterization of microsatellite loci for Ficus hirta (Moraceae).

PREMISE OF THE STUDY: Microsatellite primers were developed to investigate population genetic structure in Ficus hirta (Moraceae). METHODS AND RESULTS: Sixteen microsatellite primers were developed and optimized for F. hirta using Illumina paired-end sequencing of pre-receptive and receptive developmental-phase female flowers. Out of 16 primers, nine were found to be polymorphic in four populations of F. hirta. Alleles per locus ranged from two to 15 across the 94 F. hirta individuals, while within-population observed and expected heterozygosity per locus ranged from 0.000 to 0.955 and from 0.000 to 0.882, respectively. In addition, the 16 primers were tested in 29 additional Ficus species, with all found to amplify in at least 11 of these species and with most amplifying in a majority of the species. CONCLUSIONS: This set of microsatellite primers is the first specifically developed for F. hirta and will facilitate studies of genetic diversity within and genetic differentiation among populations of Ficus species.

Trade-offs and coexistence in fluctuating environments: evidence for a key dispersal-fecundity trade-off in five nonpollinating fig wasps.

The ecological principle of competitive exclusion states that species competing for identical resources cannot coexist, but this principle is paradoxical because ecologically similar competitors are regularly observed. Coexistence is possible under some conditions if a fluctuating environment changes the competitive dominance of species. This change in competitive dominance implies the existence of trade-offs underlying species' competitive abilities in different environments. Theory shows that fluctuating distance between resource patches can facilitate coexistence in ephemeral patch competitors, given a functional trade-off between species dispersal ability and fecundity. We find evidence supporting this trade-off in a guild of five ecologically similar nonpollinating fig wasps and subsequently predict local among-patch species densities. We also introduce a novel colonization index to estimate relative dispersal ability among ephemeral patch competitors. We suggest that a dispersal ability-fecundity trade-off and spatiotemporally fluctuating resource availability commonly co-occur to drive population dynamics and facilitate coexistence in ephemeral patch communities.

Trade-offs and coexistence: a lottery model applied to fig wasp communities.

Ecological communities in which organisms complete their life cycles on discrete ephemeral patches are common and often support an unusually large number of species. Explaining this diversity is challenging for communities of ecologically similar species undergoing preemptive competition, where classic coexistence mechanisms may not readily apply. We use nonpollinating fig wasps as a model community characterized by high diversity and preemptive competition to show how subadditive population growth and a trade-off between competitor fecundity and dispersal ability can lead to coexistence. Because nonpollinator species are often closely related, have similar life histories, and compete for the same discrete resources, understanding their coexistence is challenging given competitive exclusion is expected. Empirical observations suggest that nonpollinating fig wasp species may face a trade-off between egg loads and dispersal abilities. We model a lottery in which a species' competitive ability is determined by a trade-off between fecundity and dispersal ability. Variation in interpatch distance between figs generates temporal variability in the relative benefit of fecundity versus dispersal. We show that the temporal storage effect leads to coexistence for a range of biologically realistic parameter values. We further use individual-based modeling to show that when species' traits evolve, coexistence is less likely but trait divergence can result. We discuss the implications of this coexistence mechanism for ephemeral patch systems wherein competition is strongly preemptive.

Ecological coassociations influence species' responses to past climatic change: an example from a Sonoran Desert bark beetle.

Ecologically interacting species may have phylogeographical histories that are shaped both by features of their abiotic landscape and by biotic constraints imposed by their coassociation. The Baja California peninsula provides an excellent opportunity to examine the influence of abiotic vs. biotic factors on patterns of diversity in plant-insect species.This is because past climatic and geological changes impacted the genetic structure of plants quite differently to that of codistributed free-living animals (e.g. herpetofauna and small mammals). Thus, 'plant-like' patterns should be discernible in host-specific insect herbivores. Here, we investigate the population history of a monophagous bark beetle, Araptus attenuatus, and consider drivers of phylogeographical patterns in the light of previous work on its host plant, Euphorbia lomelii. Using a combination of phylogenetic, coalescent-simulation-based and exploratory analyses of mitochondrial DNA sequences and nuclear genotypic data, we found that the evolutionary history of A. attenuatus exhibits similarities to its host plant that are attributable to both biotic and abiotic processes. Southward range expansion and recent colonization of continental Sonora from the Baja peninsula appear to be unique to this taxon pair and probably reflect influences of the host plant. On the other hand, abiotic factors with landscape level influences on a diverse suite of codistributed arid-adapted taxa, such as Plio- and Pleistocene-aged marine incursions in the region, also left genetic signatures in beetle and host plant populations. Superimposed on these similarities, bark beetle-specific patterns and processes were also evident: our data revealed two secondarily sympatric,reproductively isolated genetic lineages, as well as a previously unrecognized mid peninsular warm desert refuge. Taken together, this work illustrates that the evolutionary history of species-specific insect herbivores may represent a mosaic of influences,including­but not limited to­those imposed by the host plant.

Nuclear and chloroplast DNA phylogeography of Ficus hirta: obligate pollination mutualism and constraints on range expansion in response to climate change.

This study uses a phylogeographic approach to investigate how interspecific interactions in an obligate pollination mutualism enhance or constrain dispersal and the range distributions of species through time. Fifteen populations of Ficus hirta, a bird-dispersed fig pollinated by a species-specific fig wasp, were sampled from Thailand to the northern limits of the tropical forest in China. These populations were assayed for six nuclear microsatellite loci and two intergenic chloroplast DNA sequences. Analyses of range expansion and genetic clustering indicated a relatively slow rate of range expansion from two or more southern glacial refugia. Low nuclear differentiation, combined with high interpopulation differentiation, and phylogeographic structuring of chloroplast variation indicated that seed dispersal has had a greater constraint than obligate interactions with fig wasps on the rate of post-glacial range expansion. This study is the first to investigate the phylogeographic history of a widely distributed southeast Asian tropical plant whose distribution extends to the northern limits of tropical forest habitat in China. It is also the first study of Ficus utilizing molecular data to evaluate whether species-specific pollination is a limitation or an aid to range expansion in response to climate change.

Flowering asynchrony and mating system effects on reproductive assurance and mutualism persistence in fragmented fig-fig wasp populations.

PREMISE OF THE STUDY: Plants and animals may experience reproductive Allee effects in fragmented populations, and obligate pollination mutualisms may be especially sensitive to extinction risk via this density-dependent process. In this study we examine how a shift from within-crown reproductive synchrony to asynchrony influences reproductive assurance through contributions to selfing and outcrossing in small, spatially isolated populations of Ficus. METHODS: The research focuses on the monoecious fig F. petiolaris and consists of phenological analyses and genetic assessments of selfing and outcrossing for populations located in Baja California's Sonoran Desert. KEY RESULTS: Phenological censuses of eight populations revealed within-crown asynchrony in 44% of reproducing trees, with 16% having sufficient overlap of male and female flowering phases to permit selfing via the cycling of pollinating fig wasps within natal trees. In mating system analyses of two of these populations, however, multilocus outcrossing rates (t(m)) were indistinguishable from 1. This result, combined with low levels of inbreeding, indicates selfing to be absent or at best a minor contributor to reproductive assurance. CONCLUSIONS: The results indicate that the fitness benefits of within-crown asynchrony lie not with selfing, as commonly asserted, but with increased opportunities for outcross pollen transmission and receipt, changing our understanding of the mechanisms by which reproduction is facilitated and extinction risk minimized in naturally fragmented Ficus populations. Given the role of fig fruit as a keystone food resource in many tropical environments, trait variation leading to reproductive assurance in figs, such as within-crown asynchrony, has broader ecosystem-level implications.

Microsatellite loci for Gossypium davidsonii (Malvaceae) and other D-genome, Sonoran Desert endemic cotton species.

PREMISE OF THE STUDY: Microsatellite primers previously developed for domesticated cotton (Gossypium hirsutum; tetraploid) were screened for their utility in investigating genetic structure and gene flow within G. davidsonii and five other wild, Mexican, D-genome cotton species (all diploid). METHODS AND RESULTS: We screened 50 microsatellite primer pairs from the Cotton Marker Database, identifying 10 loci as polymorphic within G. davidsonii. In genotyping approximately 200 individuals from four populations, we found that the number of alleles per locus ranged from seven to 17, and mean observed and expected heterozygosity ranged from 0.145 to 0.492 and from 0.436 to 0.734, respectively. We genotyped six to 20 individuals from each of the remaining species, finding these 10 loci to cross-amplify in all cases and in most cases to be polymorphic. CONCLUSIONS: These markers may be useful for further investigation of population genetics of G. davidsonii and other wild D-genome cotton species.

Significant demographic and fine-scale genetic structure in expanding and senescing populations of the terrestrial orchid Cymbidium goeringii (Orchidaceae).

PREMISE OF THE STUDY: Fine-scale genetic structure (FSGS) in plants is influenced by variation in spatial and temporal demographic processes. To determine how demographic structure and FSGS change with stages of population succession, we studied replicate expanding and senescing populations of the Asian terrestrial orchid Cymbidium goeringii. METHODS: We used spatial autocorrelation methods (O-ring and kinship statistics) to quantify spatial demographic structure and FSGS in two expanding and two senescing populations, also measuring genetic diversity and inbreeding in each. KEY RESULTS: All populations exhibited significant aggregation of individuals and FSGS at short spatial scales. In expanding populations, this finding was associated with high recruitment rates, suggesting restricted seed dispersal. In senescing populations, recruitment was minimal, suggesting alternative mechanisms of aggregation, perhaps including spatial associations with mycorrhizal fungi. All populations had significant evidence of genetic bottlenecks, and inbreeding levels were consistently high. CONCLUSIONS: Our results indicate that different successional stages can generate similar patterns of spatial demographic and genetic structure, but as a consequence of different processes. These results contrast with the only other study of senescence effects on population genetic structure in an herbaceous perennial, which found little to no FSGS in senescing populations. With the exception of populations subject to mass collection by orchid sellers, significant FSGS is characteristic of the 16 terrestrial orchid species examined to date. From a conservation perspective, this result suggests that inference of orchid population history will benefit from analyses of both FSGS and demographic structure in combination with other ecological field data.

Development of microsatellite loci for Blastophaga javana (Agaonidae), the pollinating wasp of Ficus hirta (Moraceae).

PREMISE OF THE STUDY: To investigate genetic structure of the pollinator and, indirectly, gene flow in its plant host, Blastophaga javana microsatellite primers were developed. METHODS AND RESULTS: Ten polymorphic microsatellite loci were developed using the Fast Isolation by AFLP of Sequences Containing (FIASCO) repeats protocol. Numbers of alleles per locus ranged from 2 to 19, with observed and expected heterozygosities ranging from 0.000 to 0.800 and from 0.000 to 0.925, respectively. CONCLUSIONS: These markers may be useful for further investigation of population genetics of Blastophaga javana and other congeneric species.

Slatkin's Paradox: when direct observation and realized gene flow disagree. A case study in Ficus.

In flying insects, there is frequently a lack of congruence between empirical estimates of local demographic parameters and the prediction that differentiation between populations should decrease with increasing dispersal, a puzzling phenomenon known as Slatkin's Paradox. Here, we generalize Slatkin's Paradox to other taxa, drawing from available information on dispersal to predict the relative importance of pollen vs. seed migration in structuring broad-scale patterns of genetic variation in Ficus hirta, a dioecious fig whose pollen is dispersed by minute, species-specific fig wasps and whose seeds are disseminated by a variety of vagile vertebrates (especially bats and birds). Local-scale observational and genetic studies of dioecious understory figs suggest comparable rates of pollen and seed migration. In contrast, we found unusually low nuclear differentiation (F(ST) =0.037, R(ST) =0.074) and high chloroplast differentiation (G(ST) =0.729, N(ST) =0.798) among populations separated by up to 2850km, leading us to reject the hypothesis of equal pollen and seed migration rates and to obtain an equilibrium estimate of the ratio of pollen to seed migration of r=16.2-36.3. We reconcile this example of Slatkin's Paradox with previously published data for dioecious figs and relate it to the picture of exceptionally long-distance wasp-mediated pollen dispersal that is emerging for large monoecious fig trees. More generally, we argue that Slatkin's Paradox is a general phenomenon and suggest it may be common in plants and animals.

Landscape modelling of gene flow: improved power using conditional genetic distance derived from the topology of population networks.

Landscape genetics is a burgeoning field of interest that focuses on how site-specific factors influence the distribution of genetic variation and the genetic connectivity of individuals and populations. In this manuscript, we focus on two methodological extensions for landscape genetic analyses: the use of conditional genetic distance (cGD) derived from population networks and the utility of extracting potentially confounding effects caused by correlations between phylogeographic history and contemporary ecological factors. Individual-based simulations show that when describing the spatial distribution of genetic variation, cGD consistently outperforms the traditional genetic distance measure of linearized F(ST) under both 1- and 2-dimensional stepping stone models and Cavalli-Sforza and Edward's chord distance D(c) in 1-dimensional landscapes. To show how to identify and extract the effects of phylogeographic history prior to embarking on landscape genetic analyses, we use nuclear genotypic data from the Sonoran desert succulent Euphorbia lomelii (Euphrobiaceae), for which a detailed phylogeographic history has previously been determined. For E. lomelii, removing the effect of phylogeographic history significantly influences our ability to infer both the identity and the relative importance of spatial and bio-climatic variables in subsequent landscape genetic analyses. We close by discussing the utility of cGD in landscape genetic analyses.

Comparison of quantitative and molecular genetic variation of native vs. invasive populations of purple loosestrife (Lythrum salicaria L., Lythraceae).

Study of adaptive evolutionary changes in populations of invasive species can be advanced through the joint application of quantitative and population genetic methods. Using purple loosestrife as a model system, we investigated the relative roles of natural selection, genetic drift and gene flow in the invasive process by contrasting phenotypical and neutral genetic differentiation among native European and invasive North American populations (Q(ST) - F(ST) analysis). Our results indicate that invasive and native populations harbour comparable levels of amplified fragment length polymorphism variation, a pattern consistent with multiple independent introductions from a diverse European gene pool. However, it was observed that the genetic variation reduced during subsequent invasion, perhaps by founder effects and genetic drift. Comparison of genetically based quantitative trait differentiation (Q(ST)) with its expectation under neutrality (F(ST)) revealed no evidence of disruptive selection (Q(ST) > F(ST)) or stabilizing selection (Q(ST) < F(ST)). One exception was found for only one trait (the number of stems) showing significant sign of stabilizing selection across all populations. This suggests that there are difficulties in distinguishing the effects of nonadaptive population processes and natural selection. Multiple introductions of purple loosestrife may have created a genetic mixture from diverse source populations and increased population genetic diversity, but its link to the adaptive differentiation of invasive North American populations needs further research.