Carole Kerdelhué

An Extreme Case of Plant-Insect Codiversification: Figs and Fig-Pollinating Wasps

It is thought that speciation in phytophagous insects is often due to colonization of novel host plants, because radiations of plant and insect lineages are typically asynchronous. Recent phylogenetic comparisons have supported this model of diversification for both insect herbivores and specialized pollinators. An exceptional case where contemporaneous plant-insect diversification might be expected is the obligate mutualism between fig trees (Ficus species, Moraceae) and their pollinating wasps (Agaonidae, Hymenoptera). The ubiquity and ecological significance of this mutualism in tropical and subtropical ecosystems has long intrigued biologists, but the systematic challenge posed by >750 interacting species pairs has hindered progress toward understanding its evolutionary history. In particular, taxon sampling and analytical tools have been insufficient for large-scale cophylogenetic analyses. Here, we sampled nearly 200 interacting pairs of fig and wasp species from across the globe. Two supermatrices were assembled: on an average, wasps had sequences from 77% of 6 genes (5.6 kb), figs had sequences from 60% of 5 genes (5.5 kb), and overall 850 new DNA sequences were generated for this study. We also developed a new analytical tool, Jane 2, for event-based phylogenetic reconciliation analysis of very large data sets. Separate Bayesian phylogenetic analyses for figs and fig wasps under relaxed molecular clock assumptions indicate Cretaceous diversification of crown groups and contemporaneous divergence for nearly half of all fig and pollinator lineages. Event-based cophylogenetic analyses further support the codiversification hypothesis. Biogeographic analyses indicate that the present-day distribution of fig and pollinator lineages is consistent with a Eurasian origin and subsequent dispersal, rather than with Gondwanan vicariance. Overall, our findings indicate that the fig-pollinator mutualism represents an extreme case among plant-insect interactions of coordinated dispersal and long-term codiversification. [Biogeography; coevolution; cospeciation; host switching; long-branch attraction; phylogeny.].

The effect of RAD allele dropout on the estimation of genetic variation within and between populations

Inexpensive short‐read sequencing technologies applied to reduced representation genomes is revolutionizing genetic research, especially population genetics analysis, by allowing the genotyping of massive numbers of single‐nucleotide polymorphisms (SNP) for large numbers of individuals and populations. Restriction site–associated DNA (RAD) sequencing is a recent technique based on the characterization of genomic regions flanking restriction sites. One of its potential drawbacks is the presence of polymorphism within the restriction site, which makes it impossible to observe the associated SNP allele (i.e. allele dropout, ADO). To investigate the effect of ADO on genetic variation estimated from RAD markers, we first mathematically derived measures of the effect of ADO on allele frequencies as a function of different parameters within a single population. We then used RAD data sets simulated using a coalescence model to investigate the magnitude of biases induced by ADO on the estimation of expected heterozygosity and FST under a simple demographic model of divergence between two populations. We found that ADO tends to overestimate genetic variation both within and between populations. Assuming a mutation rate per nucleotide between 10−9 and 10−8, this bias remained low for most studied combinations of divergence time and effective population size, except for large effective population sizes. Averaging FST values over multiple SNPs, for example, by sliding window analysis, did not correct ADO biases. We briefly discuss possible solutions to filter the most problematic cases of ADO using read coverage to detect markers with a large excess of null alleles.

Non-Pollinating Afrotropical Fig Wasps Affect the Fig-Pollinator Mutualism in Ficus within the Subgenus Sycomorus

Ficus are species-specifically pollinated by chalcidoid wasps (Agaonidae), and the fig tree-fig pollinator mutualism has long been studied. A diversity of other chalcids, both gall-makers or parasitoids, co-occur in each monoecious Ficus species, and use the fig resource by developing in ovaries that they transform into galls. The oviposition sequence and impact on the mutualism of these non-pollinating fig wasps were studied on two Ficus species of the subgenus Sycomorus in the Ivory Coast. Field observations, fig measurements and counts of wasps and seeds from mature figs were conducted. Four groups of fig wasp species oviposit in the ovaries at different stages of syconial development, and were demonstrated to lay eggs in the internal ovary layers; the most external flowers seem to be protected against oviposition and thus develop into seeds. Non-pollinating wasp species, by parasitizing pollinator larvae or by competing for oviposition sites, have a significant negative impact on the pollinator population and dynamics, and thereby on Ficus male function (i.e. pollen dispersal). In contrast, but for one species, they do not affect the seed production, that is, the tree female function.

Estimation of population allele frequencies from next-generation sequencing data: pool-versus individual-based genotyping

Molecular markers produced by next‐generation sequencing (NGS) technologies are revolutionizing genetic research. However, the costs of analysing large numbers of individual genomes remain prohibitive for most population genetics studies. Here, we present results based on mathematical derivations showing that, under many realistic experimental designs, NGS of DNA pools from diploid individuals allows to estimate the allele frequencies at single nucleotide polymorphisms (SNPs) with at least the same accuracy as individual‐based analyses, for considerably lower library construction and sequencing efforts. These findings remain true when taking into account the possibility of substantially unequal contributions of each individual to the final pool of sequence reads. We propose the intuitive notion of effective pool size to account for unequal pooling and derive a Bayesian hierarchical model to estimate this parameter directly from the data. We provide a user‐friendly application assessing the accuracy of allele frequency estimation from both pool‐ and individual‐based NGS population data under various sampling, sequencing depth and experimental error designs. We illustrate our findings with theoretical examples and real data sets corresponding to SNP loci obtained using restriction site–associated DNA (RAD) sequencing in pool‐ and individual‐based experiments carried out on the same population of the pine processionary moth (Thaumetopoea pityocampa). NGS of DNA pools might not be optimal for all types of studies but provides a cost‐effective approach for estimating allele frequencies for very large numbers of SNPs. It thus allows comparison of genome‐wide patterns of genetic variation for large numbers of individuals in multiple populations.

Phylogeographic structure and past history of the circum-Mediterranean species Tomicus destruens Woll. (Coleoptera: Scolytinae).

Phylogeographic studies are often focused on temperate European species with relict footholds in the Mediterranean region. Past climatic oscillations usually induced range contractions and expansions from refugial areas located in southern Europe, and spatial distribution of genetic diversity show that northward expansions were usually pioneer‐like. Actually, few studies have focused on circum‐Mediterranean species, which probably were not influenced in the same way by climatic oscillations. We present the phylogeography of the bark beetle Tomicus destruens, which is restricted to the whole Mediterranean basin and the Atlantic coasts of North Africa and Portugal. We systematically sequenced 617 bp of the mitochondrial genes COI and COII for 42 populations (N = 219). Analysis revealed 53 haplotypes geographically structured in two clades, namely eastern and western clades, that diverged during the Pleistocene. A contact zone was identified along the Adriatic coast of Italy. Interestingly, we found contrasting levels of genetic structure within each clade. The eastern group was characterized by a significant phylogeographic pattern and low levels of gene flow, whereas the western group barely showed a spatial structure in haplotype distribution. Moreover, the main pine hosts were different between groups, with the Aleppo‐brutia complex in the east and the maritime pine in the west. Potential roles of host species, climatic parameters and geographical barriers are discussed and the phylogeographic patterns are compared to classical models of postglacial recolonization in Europe.

Population genetic structure of Tomicus piniperda L. (Curculionidae: Scolytinae) on different pine species and validation of T. destruens (Woll.)

Genetic diversity and population structure of Tomicus piniperda was assessed using mitochondrial sequences on 16 populations sampled on 6 pine species in France. Amplifications of Internal transcribed space 1 (ITS1) were also performed. Our goals were to determine the taxonomic status of the Mediterranean ecotype T. piniperda destruens, and to test for host plant or geographical isolation effect on population genetic structure. We showed that T. piniperda clusters in two mtDNA haplotypic groups. Clade A corresponds to insects sampled in continental France on Pinus sylvestris, P. pinaster and P. uncinata, whereas clade B gathers the individuals sampled in Corsica on P. pinaster and P. radiata and in continental France on P. pinea and P. halepensis. Insects belonging to clade A and clade B also consistently differ in the length of ITS1. Individuals belonging to both clades were found once in sympatry on P. pinaster. Genetic distances between clades are similar to those measured between distinct species of Tomicus. We concluded that clade B actually corresponds to the destruens ecotype and forms a good species, T. destruens. Analyses of molecular variance (amova) were conducted separately on T. destruens and T. piniperda to test for an effect of either geographical isolation or host species. Interestingly, the effect of host plant was significant for T. piniperda only, while the effect of geographical isolation was not. Pine species therefore seems to act as a significant barrier to gene flow, even if host race formation is not observed. These results still need to be confirmed by nuclear markers.

Genetic isolation through time: allochronic differentiation of a phenologically atypical population of the pine processionary moth.

Allochronic speciation refers to a mode of sympatric speciation in which the differentiation of populations is primarily due to a phenological shift without habitat or host change. However, it has been so far rarely documented. The present paper reports on a plausible case of allochronic differentiation between sympatric populations of the pine processionary moth (PPM), Thaumetopoea pityocampa. The PPM is a Mediterranean insect with winter larval development. A phenologically atypical population with early adult activity and summer larval development was detected 10 years ago in Portugal. Mitochondrial and nuclear sequences strongly suggest that the ‘summer’ individuals are closely related to the sympatric winter population, while microsatellite data show a reduction in allelic richness, a distortion of allelic frequencies and significant genetic differentiation. Moreover, monitoring of adult flights suggests that reproductive activity does not overlap between the summer and winter populations. We postulate that the summer population appeared after a sudden phenological shift of some individuals of the sympatric winter population, leading to a founder effect and complete reproductive isolation. Given that the individuals showing this new phenology are subject to different selection pressures, the observed allochronic differentiation may rapidly lead to deeper divergence.

Population Genetics of the Oriental Fruit Fly, Bactrocera dorsalis (Diptera: Tephritidae), in Yunnan (China) Based on Mitochondrial DNA Sequences

Abstract The oriental fruit fly, Bactrocera dorsalis Hendel, is one of the most destructive pest insects of tropical and subtropical fruits and vegetables. It is thought to be an introduced species in Yunnan Province, China, where it causes severe damage. Depending on the latitude, the fly occurs year-round or only during the warm season. To assess the genetic diversity of the fly and to understand the relative isolation of its populations in this mountainous region, we conducted an analysis of population genetic structure using mitochondrial cytochrome oxidase (COI) gene sequences. Twenty-eight haplotypes were detected among 37 individuals with up to 13 mutations between haplotypes. Within-population diversity was high, and genetic distances between haplotypes reached 2.2%. The haplotype network showed that many haplotypes were missing in the sampled populations. Intraspecific variability in Bactrocera dorsalis was thus high in Yunnan. The data suggested either a longer residence of the fly in Yunnan than recognized previously or a recurrent colonization process from different origins. One population, namely Ruili, was significantly isolated from the others, probably because of geographic barriers to gene flow. This population seemed to be in a contact zone with flies originating from surrounding regions. In contrast, some populations separated by >300 km were not significantly structured. We suggest that the insects engage in long range dispersal, most probably taking advantage of prevailing air currents. The data also suggested that the region of Kunming, where the fly only occurs seasonally, is recolonized each year by migrating flies from several southern regions.

Phylogeography of the pine processionary moth Thaumetopoea wilkinsoni in the Near East

Phylogeographic structure of the eastern pine processionary moth Thaumetopoea wilkinsoni was explored in this study by means of nested clade phylogeographic analyses of COI and COII sequences of mitochondrial DNA and Bayesian estimates of divergence times. Intraspecific relationships were inferred and hypotheses tested to understand historical spread patterns and spatial distribution of genetic variation. Analyses revealed that all T. wilkinsoni sequences were structured in three clades, which were associated with two major biogeographic events, the colonization of the island of Cyprus and the separation of southwestern and southeastern Anatolia during the Pleistocene. Genetic variation in populations of T. wilkinsoni was also investigated using amplified fragment length polymorphisms and four microsatellite loci. Contrasting nuclear with mitochondrial data revealed recurrent gene flow between Cyprus and the mainland, related to the long‐distance male dispersal. In addition, a reduction in genetic variability was observed at both mitochondrial and nuclear markers at the expanding boundary of the range, consistent with a recent origin of these populations, founded by few individuals expanding from nearby localities. In contrast, several populations fixed for one single mitochondrial haplotype showed no reduction in nuclear variability, a pattern that can be explained by recurrent male gene flow or selective sweeps at the mitochondrial level. The use of both mitochondrial and nuclear markers was essential in understanding the spread patterns and the population genetic structure of T. wilkinsoni, and is recommended to study colonizing species characterized by sex‐biased dispersal.

Incipient allochronic speciation in the pine processionary moth (Thaumetopoea pityocampa, Lepidoptera, Notodontidae).

A plausible case of allochronic differentiation, where barrier to gene flow is primarily attributable to a phenological shift, was recently discovered in Portugal for the pine processionary moth Thaumetopoea pityocampa. Previous results suggested that the observed ‘summer population’ (SP) originated from the sympatric winter population (WP). Our objectives were to finely analyse these patterns and test their stability in time, through field monitoring and genetic analyses of larvae and adults across different years. Reproductive activity never overlapped between SP and WP. Microsatellites showed a clear differentiation of the SP, consistent with a strong reduction in gene flow owing to the phenological shift. Assignment tests suggested that some individuals shift from the SP to the WP phenology, causing some hybridization. We discuss these patterns and their maintenance over time. This could be a first stage of allochronic speciation, and SP should be considered as a distinct phenological race.