Dataset documenting prevalence and counts of pine processionary moth tents on local host trees in 3 regions of France with different climatic environments

The pine processionary moth Thaumetopoea pityocampa is a defoliating lepidopter that develops during winter. The larvae are gregarious and bear urticating setae that are harmful to humans and vertebrates. They shelter in conspicuous silk tents that are easy to detect. We here present a dataset comprising tree characterization and tent counts from 3 agglomerations in France located in regions with different climatic environments. The studied trees belong to various conifer species that are potential hosts for the caterpillars. In each site, we defined clusters as one target tree and its 10-62 nearest neighbors, and surveyed each tree within the clusters by informing: tree species, coordinates, size, number of tents. We characterized a total of 3690 trees, including 2009 trees in Orléans (grouped in 68 clusters), 359 trees in La Baule (18 clusters) and 1322 trees in Montpellier (52 clusters). We provide the raw data characterizing each individual tree, graphs showing the prevalence and mean number of tents for the tree species included in the survey, and maps allowing to locate each tree. This dataset brings information about host preference of the pine processionary moth and will be useful as a baseline to study spatio-temporal variability of host-insect relationships. It can also be informative for decision-makers and managers of urban greenings to avoid trees that are likely to be heavily infested for plantation in proximity to vulnerable people.

Intraspecific niche models for the invasive ambrosia beetle Xylosandrus crassiusculus suggest contrasted responses to climate change

Xylosandrus crassiusculus is an invasive ambrosia beetle comprising two differentiated genetic lineages, named cluster 1 and cluster 2. These lineages invaded different parts of the world at different periods of time. We tested whether they exhibited different climatic niches using Schoener's D and Hellinger's I indices and modeled their current potential geographical ranges using the Maxent algorithm. The resulting models were projected according to future and recent past climate datasets for Europe and the Mediterranean region. The future projections were performed for the periods 2041-2070 and 2071-2100 using 3 SSPs and 5 GCMs. The genetic lineages exhibited different climate niches. Parts of Europe, the Americas, Sub-Saharan Africa, Asia, and Oceania were evaluated as suitable for cluster 1. Parts of Europe, South America, Central and South Africa, Asia, and Oceania were considered as suitable for cluster 2. Models projection under future climate scenarios indicated a decrease in climate suitability in Southern Europe and an increase in North Eastern Europe in 2071-2100. Most of Southern and Western Europe was evaluated as already suitable for both clusters in the early twentieth century. Our results show that large climatically suitable regions still remain uncolonized and that climate change will affect the geographical distribution of climatically suitable areas. Climate conditions in Europe were favorable in the twentieth century, suggesting that the recent colonization of Europe is rather due to an increase in propagule pressure via international trade than to recent environmental changes.

Sequence, assembly and count datasets of viruses associated to the pine processionary moth Thaumetopoea pityocampa (Denis & Schiffermüller) (Lepidoptera, Notodontidae) identified from transcriptomic high-throughput sequencing

The pine processionary moth Thaumetopoea pityocampa is a Lepidopteran pest species occurring in the Western Mediterranean. It causes heavy pine defoliations and it is a public and animal health concern because of its urticating caterpillars. Very little is known about the viruses associated to this species, as only two viruses were described so far. We here present a dataset corresponding to 34 viral transcripts, among which 27 could be confidently assigned to 9 RNA and DNA viral families (Iflaviridae, Reoviridae, Partitiviridae, Permutotetraviridae, Flaviviridae, Rhabdoviridae, Parvoviridae, Baculoviridae and PolyDNAviridae). These transcripts were identified from an original transcriptome assembled for the insect host, using both blast search and phylogenetic approaches. The data were acquired from 2 populations in Portugal and 2 populations in Italy. The transcripts were de novo assembled and used to identify viral sequences by homology searches. We also provide information about the populations and life stages in which each virus was identified. The data produced will allow to enrich the virus taxonomy in Lepidopteran hosts, and to develop PCR-based diagnostic tools to screen colonies across the range and determine the distribution and prevalence of the identified viral species.

The worldwide invasion history of a pest ambrosia beetle inferred using population genomics

Xylosandrus crassiusculus, a fungus-farming wood borer native to Southeastern Asia, is the most rapidly spreading invasive ambrosia species worldwide. Previous studies focusing on its genetic structure suggested the existence of cryptic genetic variation in this species. Yet, these studies used different genetic markers, focused on different geographical areas and did not include Europe. Our first goal was to determine the worldwide genetic structure of this species based on both mitochondrial and genomic markers. Our second goal was to study X. crassiusculus' invasion history on a global level and identify the origins of the invasion in Europe. We used a COI and RAD sequencing design to characterize 188 and 206 specimens worldwide, building the most comprehensive genetic data set for any ambrosia beetle to date. The results were largely consistent between markers. Two differentiated genetic clusters were invasive, albeit in different regions of the world. The markers were inconsistent only for a few specimens found exclusively in Japan. Mainland USA could have acted as a source for further expansion to Canada and Argentina through stepping stone expansion and bridgehead events. We showed that Europe was only colonized by Cluster II through a complex invasion history including several arrivals from multiple origins in the native area, and possibly including bridgehead from the United States. Our results also suggested that Spain was colonized directly from Italy through intracontinental dispersion. It is unclear whether the mutually exclusive allopatric distribution of the two clusters is due to neutral effects or due to different ecological requirements.

Complete mitogenome data from a summer population specimen of the urticating pine defoliator Thaumetopoea pityocampa (Lepidoptera, Notodontidae, Thaumetopoeinae, Thaumetopoea)

We present a de novo mitogenome assembly obtained from specimens sampled in the so-called summer population (SP) of Thaumetopoea pityocampa (Denis and Schiffermüller, 1775) in Portugal. Contrary to the typical larval development occurring in winter in this species, the larvae of this unique population develop during summer. The sequencing data used were obtained from genomic libraries originally generated to assemble the nuclear genome of T. pityocampa[1]. We also provide a complete annotation and a phylogenetic representation which positions the Portuguese summer population of T. pityocampa and an Italian typical individual of the same species among the Notodontidae family and more distant Noctuoidea species. This data represents a valuable new resource for an expanding and urticating insect pest.

Role of Genes in Regulating Host Plants Expansion in Tephritid Fruit Flies (Diptera) and Potential for RNAi-Based Control

Host plant expansion is an important survival strategy for tephritids as they expand their range. Successful host expansion requires tephritids to adapt to the chemical and nonchemical properties of a novel host fruit, such as fruit color, phenology, and phytochemicals. These plant properties trigger a series of processes in tephritids, with each process having its own genetic basis, which means that various genes are involved in regulating host plant expansion by tephritids. This review summarizes current knowledge on the categories and roles of genes involved in host plant expansion in several important tephritid species, including genes related to chemoreception (olfactory and gustation), vision, digestion, detoxification, development, ribosomal and energy metabolism. Chemoreception- and detoxification- and digestion-related genes are stimulated by volatile chemicals and secondary chemicals of different hosts, respectively, which are involved in the regulation of nervous signal transduction that triggers behavioral, physical, and chemical responses to the novel host fruit. Vision-, nerve-, and development-related genes and metabolism-associated genes are activated in response to nonchemical stimuli from different hosts, such as color and phenology, to regulate a comprehensive adaptation of the extending host for tephritids. The chemical and nonchemical signals of hosts activate ribosomal and energy-related genes that result in the basic regulation of many processes of host expansion, including detoxification and development. These genes do not regulate novel host use individually, but multiple genes regulate multilevel adaptation to novel host fruits via multiple mechanisms. These genes may also be potential target genes for RNAi-based control of tephritid pests.

Monitoring of adult emergence in the pine processionary moth between 1970 and 1984 in Mont Ventoux, France

Background

The current climate change has marked impacts on the phenology of species, i.e. the timing of the various stages of their life cycle. Yet, to fully understand how phenological patterns can be modified according to changes in temperature regimes, it is of prime importance to rely on high quality historical data. Here, we propose a very valuable dataset including individual monitoring from pupation to adult emergence of 46 479 individuals of pine processionary moth (Thaumetopoea pityocampa) surveyed between 1970 and 1984 in southern France along an altitudinal gradient. As optional prolonged diapause occurs in this species, i.e. some individuals experience one or more years of diapause before emerging, the caterpillars sampled in any given year were monitored during up to 5 years. The goal was to give precise information about phenology in this species to further analyse its temporal patterns of variation.

New information

This dataset is unique by its richness and the type of data it contains. Phenology in the pine processionary moth is often monitored by the use of pheromone traps in the field, which does not provide all the necessary information, because it is then not possible to trace back the exact origin of the moth trapped, nor to characterise other steps of the life cycle. Moreover, as it corresponds to historical data dating back to the 70s and the 80s, the dataset provides a historical baseline of trends in the pre-warming period.

From refugia to contact: Pine processionary moth hybrid zone in a complex biogeographic setting

Contact zones occur at the crossroad between specific dispersal routes and are facilitated by biogeographic discontinuities. Here, we focused on two Lepidoptera sister species that come in contact near the Turkish Straits System (TSS). We aimed to infer their phylogeographic histories in the Eastern Mediterranean and finely analyze their co-occurrence and hybridization patterns in this biogeographic context. We used molecular mitochondrial and nuclear markers to study 224 individuals from 42 localities. We used discordances between markers and complementary assignment methods to identify and map hybrids and parental individuals. We confirmed the parapatric distribution of Thaumetopoea pityocampa (Lepidoptera: Notodontidae) in the west and Thaumetopoea wilkinsoni in the east and identified a narrow contact zone. We identified several glacial refugia of T. wilkinsoni in southern Turkey with a strong east-west differentiation in this species. Unexpectedly, T. pityocampa crossed the TSS and occur in northern Aegean Turkey and some eastern Greek islands. We found robust evidence of introgression between the two species in a restricted zone in northwestern Turkey, but we did not identify any F1 individuals. The identified hybrid zone was mostly bimodal. The distributions and genetic patterns of the studied species were strongly influenced both by the Quaternary climatic oscillations and the complex geological history of the Aegean region. T. pityocampa and T. wilkinsoni survived the last glacial maximum in disjoint refugia and met in western Turkey at the edge of the recolonization routes. Expanding population of T. wilkinsoni constrained T. pityocampa to the western Turkish shore. Additionally, we found evidence of recurrent introgression by T. wilkinsoni males in several T. pityocampa populations. Our results suggest that some prezygotic isolation mechanisms, such as differences in timing of the adult emergences, might be a driver of the isolation between the sister species.

Shifted phenology in the pine processionary moth affects the outcome of tree-insect interaction

In the Mediterranean and temperate regions, an increase in the frequency and intensity of drought events has been recorded, probably due to climate change. In consequence, trees will more frequently experience hydric stress, a condition that can be expected to affect insect-tree interactions, while adaptation mechanisms may be further in course. The effect of tree water stress on the performance of two allochronic populations of Thaumetopoea pityocampa was here studied. Namely, we compared a unique population of this insect, in which the larvae develop in the summer (SP), with the typical population having winter larval development (WP), to test the adaptation hypothesis to host plant status. Larvae of each population were fed on needles of young potted Pinus pinaster plants under two water supply regimes: (i) well-watered (control) and (ii) subjected to 3 months of drought stress. Compared to control, stressed plants had higher amounts of soluble sugars, phenols, and higher C/N ratio, whereas water content and chlorophylls concentrations were lower. In general, T. pityocampa larvae had lower performances on water-stressed plants, as shown by lower survival rates, lower needle consumption, and longer development times. Yet, the detrimental effects of tree stress were only significant for the WP larvae, while SP larvae were able to overcome such conditions. Results demonstrate that tree water stress can negatively affect T. pityocampa populations. Furthermore, the evidence is also provided that responses to the physiological condition of the host trees may occur at the population level, as a result of adaptation mechanisms driven by climate change.

Draft genome and reference transcriptomic resources for the urticating pine defoliator Thaumetopoea pityocampa (Lepidoptera: Notodontidae)

The pine processionary moth Thaumetopoea pityocampa (Lepidoptera: Notodontidae) is the main pine defoliator in the Mediterranean region. Its urticating larvae cause severe human and animal health concerns in the invaded areas. This species shows a high phenotypic variability for various traits, such as phenology, fecundity and tolerance to extreme temperatures. This study presents the construction and analysis of extensive genomic and transcriptomic resources, which are an obligate prerequisite to understand their underlying genetic architecture. Using a well-studied population from Portugal with peculiar phenological characteristics, the karyotype was first determined and a first draft genome of 537 Mb total length was assembled into 68,292 scaffolds (N50 = 164 kb). From this genome assembly, 29,415 coding genes were predicted. To circumvent some limitations for fine-scale physical mapping of genomic regions of interest, a 3X coverage BAC library was also developed. In particular, 11 BACs from this library were individually sequenced to assess the assembly quality. Additionally, de novo transcriptomic resources were generated from various developmental stages sequenced with HiSeq and MiSeq Illumina technologies. The reads were de novo assembled into 62,376 and 63,175 transcripts, respectively. Then, a robust subset of the genome-predicted coding genes, the de novo transcriptome assemblies and previously published 454/Sanger data were clustered to obtain a high-quality and comprehensive reference transcriptome consisting of 29,701 bona fide unigenes. These sequences covered 99% of the cegma and 88% of the busco highly conserved eukaryotic genes and 84% of the busco arthropod gene set. Moreover, 90% of these transcripts could be localized on the draft genome. The described information is available via a genome annotation portal (http://bipaa.genouest.org/sp/thaumetopoea_pityocampa/).

De novo genome and transcriptome resources of the Adzuki bean borer Ostrinia scapulalis (Lepidoptera: Crambidae)

We present a draft genome assembly with a de novo prediction and automated functional annotation of coding genes, and a reference transcriptome of the Adzuki bean borer, Ostrinia scapulalis, based on RNA sequencing of various tissues and developmental stages. The genome assembly spans 419 Mb, has a GC content of 37.4% and includes 26,120 predicted coding genes. The reference transcriptome holds 33,080 unigenes and contains a high proportion of a set of genes conserved in eukaryotes and arthropods, used as quality assessment of the reconstructed transcripts. The new genomic and transcriptomic data presented here significantly enrich the public sequence databases for the Crambidae and Lepidoptera, and represent useful resources for future researches related to the evolution and the adaptation of phytophagous moths. The genome and transcriptome assemblies have been deposited and made accessible via a NCBI BioProject (id PRJNA390510) and the LepidoDB database (http://bipaa.genouest.org/sp/ostrinia_scapulalis/).

Proteome Analysis of Urticating Setae From Thaumetopoea pityocampa (Lepidoptera: Notodontidae)

Thaumetopoea pityocampa (Denis & Schiffermüller) (Lepidoptera: Notodontidae) is harmful to conifer trees because of defoliation and to public health because of the release of urticating setae from the caterpillars. Contact with setae by humans and domestic animals induces dermatitis, usually localized to the exposed areas. Recent studies demonstrated the presence of a complex urticating mechanism where proteins present in the setae may play a role as activators of immune responses. Yet, limited information is available at present about the proteins occurring in the setae of T. pityocampa. Using a refined method for protein extraction from the setae, and a combination of liquid chromatography tandem-mass spectrometry (LC-MS/MS), de novo assembly of transcriptomic data, and sequence similarity searches, an extensive data set of 353 proteins was obtained. These were further categorized by molecular function, biological process, and cellular location. All the 353 proteins identified were found to match through BLAST search with at least one Lepidoptera sequence available in databases. We found the previously known allergens Tha p 1 and Tha p 2 described from T. pityocampa, as well as enzymes involved in chitin biosynthesis, one of the principal components of the setae, and serine proteases that were responsible for inflammatory and allergic reactions in other urticating Lepidoptera. This new proteomic database may allow for a better understanding of the complexity of allergenic reactions due to T. pityocampa and to other Lepidoptera sharing similar defense systems.

Effect of heat waves on embryo mortality in the pine processionary moth

Extreme climate events such as heat waves are predicted to become more frequent with climate change, representing a challenge for many organisms. The pine processionary moth Thaumetopoea pityocampa is a Mediterranean pine defoliator, which typically lays eggs during the summer. We evaluated the effects of heat waves on egg mortality of three populations with different phenologies: a Portuguese population with a classical life cycle (eggs laid in summer), an allochronic Portuguese population reproducing in spring, and a Tunisian population from the extreme southern limit of T. pityocampa distribution range, in which eggs are laid in fall. We tested the influence of three consecutive hot days on egg survival and development time, using either constant (CT) or daily cycling temperatures (DT) with equivalent mean temperatures. Maximum temperatures (T max) used in the experiment ranged from 36 to 48°C for DT and from 30 to 42°C for CT. Heat waves had a severe negative effect on egg survival when T max reached 42°C for all populations. No embryo survived above this threshold. At high mean temperatures (40°C), significant differences were observed between populations and between DT and CT regimes. Heat waves further increased embryo development time. The knowledge we gained about the upper lethal temperature to embryos of this species will permit better prediction of the potential expansion of this insect under different climate warming scenarios.

Experimental evidence for heritable reproductive time in 2 allochronic populations of pine processionary moth

Phenology allows organisms to overcome seasonally variable conditions through life-cycle adjustment. Changes in phenology can drastically modify the evolutionary trajectory of a population, while a shift in the reproductive time may cause allochronic differentiation. The hypothesis of heritable reproductive time was experimentally tested, by studying a unique population of the pine processionary moth Thaumetopoea pityocampa (Den. & Schiff.) which has a shifted phenology, and however co-occurs with the typical population following the classical life cycle. When populations of both types were reared under controlled conditions, the reproductive time was maintained asynchronous, as observed in the field. The shifted population was manipulated in the laboratory to reproduce later than usual, yet the offspring emerged in the next year at the expected dates thus "coming back" to the usual cycle. Hybrids from crosses performed between the 2 populations showed an intermediate phenology. From the emergence times of parents and offspring, a high heritability of the reproductive time (h = 0.76) was observed. The offspring obtained from each type of cross was genetically characterized using microsatellite markers. Bayesian clustering analysis confirmed that hybrids can be successfully identified and separated from the parental genetic classes by genotyping. Findings support the hypothesis that, for this particular population, incipient allochronic speciation is due to a heritable shift in the reproductive time that further causes assortative mating and might eventually cause ecological adaptation/maladaptation in response to environmental changes.

Genetic differentiation of the pine processionary moth at the southern edge of its range: contrasting patterns between mitochondrial and nuclear markers

The pine processionary moth (Thaumetopoea pityocampa) is an important pest of coniferous forests at the southern edge of its range in Maghreb. Based on mitochondrial markers, a strong genetic differentiation was previously found in this species between western (pityocampa clade) and eastern Maghreb populations (ENA clade), with the contact zone between the clades located in Algeria. We focused on the moth range in Algeria, using both mitochondrial (a 648 bp fragment of the tRNA‐cox2) and nuclear (11 microsatellite loci) markers. A further analysis using a shorter mtDNA fragment and the same microsatellite loci was carried out on a transect in the contact zone between the mitochondrial clades. Mitochondrial diversity showed a strong geographical structure and a well‐defined contact zone between the two clades. In particular, in the pityocampa clade, two inner subclades were found whereas ENA did not show any further structure. Microsatellite analysis outlined a different pattern of differentiation, with two main groups not overlapping with the mitochondrial clades. The inconsistency between mitochondrial and nuclear markers is probably explained by sex‐biased dispersal and recent afforestation efforts that have bridged isolated populations.

Genetic structure and colonization history of the fruit fly Bactrocera tau (Diptera: Tephritidae) in China and Southeast Asia

Bactrocera tau (Walker), a major invasive pest worldwide, was first described in Fujian (China) in 1849 and has dispersed to tropical and subtropical Asia and the South Pacific region. Few data are available on its colonization history and expansion processes. This pilot study attempted to reconstruct the colonization history and pathways of this pest in China and neighboring Southeast Asian countries based on mitochondrial DNA. Results of the study showed six genetic groups corresponding to geographical characteristics, although the pattern was relatively weak. Homogeneous genetic patterns were observed within southern and central China, and northern Vietnam. Continuous colonization from the coast of southern China to inland regions of China and northern Vietnam was suggested. Strong genetic structure was observed in western China, Thailand, and Laos. The isolation of four of the six groups was most probably attributable to major topographical barriers of western China. Yunnan acted as a contact zone of B. tau in China and neighboring Southeast Asia. The absence of isolation by distance and the overall low phylogeographic structure of B. tau suggested that long distance dispersal events and human activities could play a major role in the colonization and expansion patterns of B. tau. By analyzing the genetic diversity, gene flow, haplotype phylogeny, and demographic history of 23 fly populations, we hypothesized that B. tau could have been introduced long ago in southern China, from which it further expanded or that southern China could correspond to the native range of this species.

Contrasted invasion processes imprint the genetic structure of an invasive scale insect across southern Europe

Deciphering the colonization processes by which introduced pests invade new areas is essential to limit the risk of further expansion and/or multiple introductions. We here studied the invasion history of the maritime pine bast scale Matsucoccus feytaudi. This host-specific insect does not cause any damage in its native area, but it devastated maritime pine forests of South-Eastern France where it was detected in the 1960s, and since then reached Italy and Corsica. We used population genetic approaches to infer the populations' recent evolutionary history from microsatellite markers and Approximate Bayesian Computation. Consistent with previous mitochondrial data, we showed that the native range is geographically strongly structured, which is probably due to the patchy distribution of the obligate host and the limited dispersal capacity of the scale. Our results show that the invasion history can be described in three successive steps involving different colonization and dispersal processes. During the mid-XXth century, massive introductions occurred from the Landes planted forest to South-Eastern France, probably due to transportation of infested wood material after World War II. Stepping-stone expansion, consistent with natural dispersal, then allowed M. feytaudi to reach the maritime pine forests of Liguria and Tuscany in Italy. The island of Corsica was accidentally colonized in the 1990s, and the most plausible scenario involves the introduction of a limited number of migrants from the forests of South-Eastern France and Liguria, which is consistent with an aerial dispersal due to the dominant winds that blow in spring in this region.

Comparative analysis of two phenologically divergent populations of the pine processionary moth (Thaumetopoea pityocampa) by de novo transcriptome sequencing

The pine processionary moth Thaumetopoea pityocampa is a Mediterranean lepidopteran defoliator that experiences a rapid range expansion towards higher latitudes and altitudes due to the current climate warming. Its phenology - the time of sexual reproduction - is certainly a key trait for the local adaptation of the processionary moth to climatic conditions. Moreover, an exceptional case of allochronic differentiation was discovered ca. 15 years ago in this species. A population with a shifted phenology (the summer population, SP) co-exists near Leiria, Portugal, with a population following the classical cycle (the winter population, WP). The existence of this population is an outstanding opportunity to decipher the genetic bases of phenology. No genomic resources were so far available for T. pityocampa. We developed a high-throughput sequencing approach to build a first reference transcriptome, and to proceed with comparative analyses of the sympatric SP and WP. We pooled RNA extracted from whole individuals of various developmental stages, and performed a transcriptome characterisation for both populations combining Roche 454-FLX and traditional Sanger data. The obtained sequences were clustered into ca. 12,000 transcripts corresponding to 9265 unigenes. The mean transcript coverage was 21.9 reads per bp. Almost 70% of the de novo assembled transcripts displayed significant similarity to previously published proteins and around 50% of the transcripts contained a full-length coding region. Comparative analyses of the population transcriptomes allowed to investigate genes specifically expressed in one of the studied populations only, and to identify the most divergent homologous SP/WP transcripts. The most divergent pairs of transcripts did not correspond to obvious phenology-related candidate genes, and 43% could not be functionally annotated. This study provides the first comprehensive genome-wide resource for the target species T. pityocampa. Many of the assembled genes are orthologs of published Lepidoptera genes, which allows carrying out gene-specific re-sequencing. Data mining has allowed the identification of SNP loci that will be useful for population genomic approaches and genome-wide scans of population differentiation to identify signatures of selection.

Phenotypic divergence in reproductive traits of a moth population experiencing a phenological shift

Allochrony that is reproductive isolation by time may further lead to divergence of reproductive adaptive traits in response to different environmental pressures over time. A unique "summer" population of the pine processionary moth Thaumetopoea pityocampa, reproductively isolated from the typical winter populations by allochronic differentiation, is here analyzed. This allochronically shifted population reproduces in the spring and develops in the summer, whereas "winter" populations reproduce in the late summer and have winter larval development. Both summer and winter populations coexist in the same pine stands, yet they face different climatic pressures as their active stages are present in different seasons. The occurrence of significant differences between the reproductive traits of the summer population and the typical winter populations (either sympatric or allopatric) is thus hypothesized. Female fecundity, egg size, egg covering, and egg parasitism were analyzed showing that the egg load was lower and that egg size was higher in the summer population than in all the studied winter populations. The scales that cover the egg batches of T. pityocampa differed significantly between populations in shape and color, resulting in a looser and darker covering in the summer population. The single specialist egg parasitoid species of this moth was almost missing in the summer population, and the overall parasitism rates were lower than in the winter population. Results suggest the occurrence of phenotypic differentiation between the summer population and the typical T. pityocampa winter populations for the life-history traits studied. This work provides an insight into how ecological divergence may follow the process of allochronic reproductive isolation.

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.

Host and phenology shifts in the evolution of the social moth genus Thaumetopoea

The genus Thaumetopoea contains the processionary moths, a group of lepidopteran associated with forest trees, well known for the social behaviour of the larvae and for carrying urticating setae. The taxonomy of the genus is partly unresolved and a phylogenetic approach is lacking. The goal of this work is to produce a phylogeny for Thaumetopoea and to identify the main traits driving the evolution of this group. Eighteen mitochondrial and three nuclear genes were fully/partly sequenced. Markers were aligned and analysed singularly or in various combinations. Phylogenetic analyses were performed according to maximum likelihood and Bayesian inference methods. Trees obtained from largest data sets provided identical topologies that received strong statistical support. Three main clades were identified within Thaumetopoea and were further supported by several signatures located in the mitochondrial tRNAs and intergenic spacers. The reference topology was used to investigate the evolution of life history traits related to biogeography, host plant, ecology, and morphology. A multigenic approach allowed to produce a robust phylogenetic analysis of the genus Thaumetopoea, with the identification of three major clades linked to different ecological and life history traits. The first clade is associated with Angiosperm host plants and has a fast spring development of larvae on young foliage. The other clades have originated by one event of host plant shift to Gymnosperm Pinaceae, which implied a longer larval developmental time due to the lower nutritional quality of leaves. These clades showed different adaptations to such a constraint, the first with a switch of larval feeding to cold season (winter pine processionary moths), and the second with a retraction to high altitude and latitude and a development cycle extended over two years (summer pine processionary moths). Recent global warming is affecting all species and seems able to further shape the evolution of the group.

Host and phenology shifts in the evolution of the social moth genus Thaumetopoea

The genus Thaumetopoea contains the processionary moths, a group of lepidopteran associated with forest trees, well known for the social behaviour of the larvae and for carrying urticating setae. The taxonomy of the genus is partly unresolved and a phylogenetic approach is lacking. The goal of this work is to produce a phylogeny for Thaumetopoea and to identify the main traits driving the evolution of this group. Eighteen mitochondrial and three nuclear genes were fully/partly sequenced. Markers were aligned and analysed singularly or in various combinations. Phylogenetic analyses were performed according to maximum likelihood and Bayesian inference methods. Trees obtained from largest data sets provided identical topologies that received strong statistical support. Three main clades were identified within Thaumetopoea and were further supported by several signatures located in the mitochondrial tRNAs and intergenic spacers. The reference topology was used to investigate the evolution of life history traits related to biogeography, host plant, ecology, and morphology. A multigenic approach allowed to produce a robust phylogenetic analysis of the genus Thaumetopoea, with the identification of three major clades linked to different ecological and life history traits. The first clade is associated with Angiosperm host plants and has a fast spring development of larvae on young foliage. The other clades have originated by one event of host plant shift to Gymnosperm Pinaceae, which implied a longer larval developmental time due to the lower nutritional quality of leaves. These clades showed different adaptations to such a constraint, the first with a switch of larval feeding to cold season (winter pine processionary moths), and the second with a retraction to high altitude and latitude and a development cycle extended over two years (summer pine processionary moths). Recent global warming is affecting all species and seems able to further shape the evolution of the group.

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.].

Microsatellite and mitochondrial data provide evidence for a single major introduction for the Neartic leafhopper Scaphoideus titanus in Europe

Scaphoideus titanus, a leafhopper native to North America and invasive in Europe, is the vector of the Flavescence dorée phytoplasma, the causal agent of the most important form of grapevine yellows in European vineyards. We studied 10 polymorphic microsatellite loci and a 623 bp fragment of the mitochondrial cytochrome oxidase II gene in native S. titanus from north-eastern America and introduced European populations, to elucidate the colonization scenario. Consistent with their recent history, invasive European populations were less genetically diverse than American populations for both types of markers, suggesting a recent bottleneck. Significant isolation by distance was detected between American populations but not between European populations. None of the European mitochondrial haplotypes was found in the American vineyards, from which they are assumed to have originated. The precise source of the invasive S. titanus populations therefore remains unclear. Nevertheless, the high heterozygosity of North-East American populations (which contained 92% of the observed alleles) suggests that this region is part of the native range of S. titanus. Clustering population genetics analyses with microsatellite and mitochondrial data suggested that European populations originated from a single introduction event. Most of the introduced populations clustered with populations from Long Island, the Atlantic Coast winegrowing region in which Vitis aestivalis occurs.

Genetic structure and inferences on potential source areas for Bactrocera dorsalis (Hendel) based on mitochondrial and microsatellite markers

Bactrocera dorsalis (Diptera: Tephritidae) is mainly distributed in tropical and subtropical Asia and in the Pacific region. Despite its economic importance, very few studies have addressed the question of the wide genetic structure and potential source area of this species. This pilot study attempts to infer the native region of this pest and its colonization pathways in Asia. Combining mitochondrial and microsatellite markers, we evaluated the level of genetic diversity, genetic structure, and the gene flow among fly populations collected across Southeast Asia and China. A complex and significant genetic structure corresponding to the geographic pattern was found with both types of molecular markers. However, the genetic structure found was rather weak in both cases, and no pattern of isolation by distance was identified. Multiple long-distance dispersal events and miscellaneous host selection by this species may explain the results. These complex patterns may have been influenced by human-mediated transportation of the pest from one area to another and the complex topography of the study region. For both mitochondrial and microsatellite data, no signs of bottleneck or founder events could be identified. Nonetheless, maximal genetic diversity was observed in Myanmar, Vietnam and Guangdong (China) and asymmetric migration patterns were found. These results provide indirect evidence that the tropical regions of Southeast Asia and southern coast of China may be considered as the native range of the species and the population expansion is northward. Yunnan (China) is a contact zone that has been colonized from different sources. Regions along the southern coast of Vietnam and China probably served to colonize mainly the southern region of China. Southern coastal regions of China may also have colonized central parts of China and of central Yunnan.

Permanent genetic resources added to Molecular Ecology Resources Database 1 August 2011-30 September 2011

This article documents the addition of 299 microsatellite marker loci and nine pairs of single-nucleotide polymorphism (SNP) EPIC primers to the Molecular Ecology Resources (MER) Database. Loci were developed for the following species: Alosa pseudoharengus, Alosa aestivalis, Aphis spiraecola, Argopecten purpuratus, Coreoleuciscus splendidus, Garra gotyla, Hippodamia convergens, Linnaea borealis, Menippe mercenaria, Menippe adina, Parus major, Pinus densiflora, Portunus trituberculatus, Procontarinia mangiferae, Rhynchophorus ferrugineus, Schizothorax richardsonii, Scophthalmus rhombus, Tetraponera aethiops, Thaumetopoea pityocampa, Tuta absoluta and Ugni molinae. These loci were cross-tested on the following species: Barilius bendelisis, Chiromantes haematocheir, Eriocheir sinensis, Eucalyptus camaldulensis, Eucalyptus cladocalix, Eucalyptus globulus, Garra litaninsis vishwanath, Garra para lissorhynchus, Guindilla trinervis, Hemigrapsus sanguineus, Luma chequen. Guayaba, Myrceugenia colchagüensis, Myrceugenia correifolia, Myrceugenia exsucca, Parasesarma plicatum, Parus major, Portunus pelagicus, Psidium guayaba, Schizothorax richardsonii, Scophthalmus maximus, Tetraponera latifrons, Thaumetopoea bonjeani, Thaumetopoea ispartensis, Thaumetopoea libanotica, Thaumetopoea pinivora, Thaumetopoea pityocampa ena clade, Thaumetopoea solitaria, Thaumetopoea wilkinsoni and Tor putitora. This article also documents the addition of nine EPIC primer pairs for Euphaea decorata, Euphaea formosa, Euphaea ornata and Euphaea yayeyamana.

Phylogeny and evolution of life-history strategies in the Sycophaginae non-pollinating fig wasps (Hymenoptera, Chalcidoidea)

BACKGROUND

Non-pollinating Sycophaginae (Hymenoptera, Chalcidoidea) form small communities within Urostigma and Sycomorus fig trees. The species show differences in galling habits and exhibit apterous, winged or dimorphic males. The large gall inducers oviposit early in syconium development and lay few eggs; the small gall inducers lay more eggs soon after pollination; the ostiolar gall-inducers enter the syconium to oviposit and the cleptoparasites oviposit in galls induced by other fig wasps. The systematics of the group remains unclear and only one phylogeny based on limited sampling has been published to date. Here we present an expanded phylogeny for sycophagine fig wasps including about 1.5 times the number of described species. We sequenced mitochondrial and nuclear markers (4.2 kb) on 73 species and 145 individuals and conducted maximum likelihood and Bayesian phylogenetic analyses. We then used this phylogeny to reconstruct the evolution of Sycophaginae life-history strategies and test if the presence of winged males and small brood size may be correlated.

RESULTS

The resulting trees are well resolved and strongly supported. With the exception of Apocrytophagus, which is paraphyletic with respect to Sycophaga, all genera are monophyletic. The Sycophaginae are divided into three clades: (i) Eukoebelea; (ii) Pseudidarnes, Anidarnes and Conidarnes and (iii) Apocryptophagus, Sycophaga and Idarnes. The ancestral states for galling habits and male morphology remain ambiguous and our reconstructions show that the two traits are evolutionary labile.

CONCLUSIONS

The three main clades could be considered as tribes and we list some morphological characters that define them. The same biologies re-evolved several times independently, which make Sycophaginae an interesting model to test predictions on what factors will canalize the evolution of a particular biology. The ostiolar gall-inducers are the only monophyletic group. In 15 Myr, they evolved several morphological adaptations to enter the syconia that make them strongly divergent from their sister taxa. Sycophaginae appears to be another example where sexual selection on male mating opportunities favored winged males in species with small broods and wingless males in species with large broods. However, some species are exceptional in that they lay few eggs but exhibit apterous males, which we hypothesize could be due to other selective pressures selecting against the re-appearance of winged morphs.

Quaternary history and contemporary patterns in a currently expanding species

Background

Quaternary climatic oscillations had dramatic effects on species evolution. In northern latitudes, populations had to survive the coldest periods in refugial areas and recurrently colonized northern regions during interglacials. Such a history usually results in a loss of genetic diversity. Populations that did not experience glaciations, in contrast, probably maintained most of their ancestral genetic diversity. These characteristics dramatically affected the present-day distribution of genetic diversity and may influence the ability of species to cope with the current global changes. We conducted a range-wide study of mitochondrial genetic diversity in the pine processionary moth (Thaumetopoea pityocampa/T. wilkinsoni complex, Notodontidae), a forest pest occurring around the Mediterranean Basin and in southern Europe. This species is responding to the current climate change by rapid natural range expansion and can also be accidentally transported by humans. Our aim was to assess if Quaternary climatic oscillations had a different effect across the species' range and to determine if genetic footprints of contemporary processes can be identified in areas of recent introduction.

Results

We identified three main clades that were spatially structured. In most of Europe, the genetic diversity pattern was typical for species that experienced marked glaciation cycles. Except in refugia, European populations were characterized by the occurrence of one main haplotype and by a strong reduction in genetic diversity, which is expected in regions that were rapidly re-colonized when climatic conditions improved. In contrast, all other sub-clades around the Mediterranean Basin occurred in limited parts of the range and were strongly structured in space, as is expected in regions in which the impact of glaciations was limited. In such places, genetic diversity was retained in most populations, and almost all haplotypes were endemic. This pattern was extreme on remote Mediterranean islands (Crete, Cyprus, Corsica) where highly differentiated, endemic haplotypes were found. Recent introductions were typified by the existence of closely-related haplotypes in geographically distant populations, which is difficult to detect in most of Europe because of a lack of overall genetic structure.

Conclusion

In regions that were not prone to marked glaciations, recent moth introductions/expansions could be detected due to the existence of a strong spatial genetic structure. In contrast, in regions that experienced the most intense Quaternary climatic oscillations, the natural populations are not genetically structured, and contemporary patterns of population expansion remain undetected.

Complex postglacial history of the temperate bark beetle Tomicus piniperda L. (Coleoptera, Scolytinae)

Tomicus piniperda is an economically important pine bark beetle infesting European Pinus spp. stands. We sequenced and analyzed 797 bp of the mitochondrial genome from individuals obtained from 34 populations sampled throughout the European range. We obtained 36 haplotypes, from which a haplotype network was constructed. In the Iberian Peninsula, high-genetic variability was detected with numerous endemic haplotypes. In contrast, the other European populations were less diverse with a single haplotype predominating from the Pyrenees to Scandinavia. Nevertheless, even within Europe, a few populations showed significant amounts of diversity. Four groups were obtained by Spatial Analysis of Molecular Variance, illustrating the regional characteristics of the species. T. piniperda had multiple fragmented refugia in the Iberian Peninsula. These currently isolated populations only partly contributed to postglacial re-colonizations of Northern Europe during interglacials. Nevertheless, few long-range migration events up to Northern Europe were detected, mostly originating from the Pyrenees. In the rest of Europe, the phylogeographical patterns were unclear, because of repeated cycles of contraction and expansion. The genetic analysis showed one glacial refugium in North-Central Europe, whereas other refugia most likely occurred in the Southern Alps, Apennine and the Balkans. The phylogeographical pattern depicted here reflects partly the postglacial history of the beetles' main host tree P. sylvestris.

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.

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.

Comparative population genetic study of two oligophagous insects associated with the same hosts

A parallel study of the genetic structure of two oligophagous species associated with the same hosts was conducted to determine the main factors shaping the distribution of genetic diversity. The bark beetle Tomicus piniperda and the pine processionary moth (PPM) Thaumetopoea pityocampa are both associated with the genus Pinus and belong to different guilds (xylophagous vs defoliating species). The PPM is an ectophagous species that feeds on the needles of living trees, whereas T. piniperda is endophagous and bores galleries in the inner bark of weakened trees. Both species were sampled in the main regions of France, and their genetic structure was assessed after genotyping with five microsatellite markers. Populations of the PPM were significantly structured. A pattern of isolation by distance was found when distances were calculated as bypassing the Massif Central, whereas no such pattern could be found with raw geographic distances. On the contrary, most populations of T. piniperda were not differentiated. No effect of host species could be detected in either of the two species. We conclude that the two taxa have contrasting effective dispersal rates per generation, and we hypothesize that this reflects the different selection pressures acting on individual fitness via different strategies of host use.

Genetic study of the forest pest Tomicus piniperda (Col., Scolytinae) in Yunnan province (China) compared to Europe: new insights for the systematics and evolution of the genus Tomicus

The pine shoot beetle Tomicus piniperda is present throughout Eurasia. In Europe, it is considered as a secondary pest that rarely causes tree mortality, while heavy damage is observed in Yunnan Province (China) where it exhibits a novel aggregative behaviour during shoot attack. To understand why the ecological characteristics of the European and Chinese populations differ so strongly, we conducted an analysis of population genetic structure on 12 populations in Yunnan and one in JiLin using mitochondrial (COI-COII) and nuclear (ITS2 and 28S rDNA) DNA sequences, and compared the results to those obtained in France. We showed that the Yunnan populations differed markedly from French and JiLin populations. For all three markers, the genetic distances measured between the Tomicus from Yunnan and those from France were similar to distances previously observed between species. Similar distances were found between Yunnan and JiLin populations. Conversely, the distances between French and JiLin individuals were substantially lower, falling in the intraspecific range. We concluded that the individuals sampled in Yunnan belong to a new, undescribed species (Tomicus sp. nov.). We also showed that some individuals belong to the species T. brevipilosus that had never been recorded from this region before. Evolution of the genus Tomicus is discussed in the light of these new results.

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.

Species and recombination effects on DNA variability in the tomato genus

Population genetics theory predicts that strong selection for rare, beneficial mutations or against frequent, deleterious mutations reduces polymorphism at linked neutral (or weakly selected) sites. The reduction of genetic variation is expected to be more severe when recombination rates are lower. In outbreeding species, low recombination rates are usually confined to certain chromosomal regions, such as centromeres and telomeres. In contrast, in predominantly selfing species, the rarity of double heterozygotes leads to a reduced effective recombination rate in the whole genome. We investigated the effects of restricted recombination on DNA polymorphism in these two cases, analyzing five Lycopersicon species with contrasting mating systems: L. chilense, L. hirsutum, L. peruvianum, L. chmielewskii, and L. pimpinellifolium, of which only the first three species have self-incompatibility alleles. In each species, we determined DNA sequence variation of five single-copy genes located in chromosomal regions with either high or low recombination rate. We found that the mating system has a highly significant effect on the level of polymorphism, whereas recombination has only a weak influence. The effect of recombination on levels of polymorphism in Lycopersicon is much weaker than in other well-studied species, including Drosophila. To explain these observations, we discuss a number of hypotheses, invoking selection, recombination, and demographic factors associated with the mating system. We also provide evidence that L. peruvianum, showing a level of polymorphism (almost 3%) that is comparable to the level of divergence in the whole genus, is the ancestral species from which the other species of the genus Lycopersicon have originated relatively recently.

Molecular phylogeny of fig wasps Agaonidae are not monophyletic

According to the present classification, the family Agaonidae contains all fig pollinators as well as five subfamilies of non-pollinating fig wasps. The molecular phylogeny of the family was reconstructed using partial sequences of the 28S rRNA (D1 and D2 domains). Our results show that the family Agaonidae is not monophyletic. As a consequence, we restrict the famiy to the pollinator clade, and assign the non-pollinating subfamilies to various chalcid families. Sycoecinae, Otitesellinae and Sycoryctinae are included in Pteromalidae, whereas Sycophaginae and Epichrysomallinae are left unclassified and will require more in-depth morphological studies. Moreover, we proved that the fig pollination syndrome evolved only once, early in group history. The resource due to the fig-pollinator mutualism has secondarily been colonized independently by different Chalcid lineages.