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.

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.

Climate change impact on the potential geographical distribution of two invading Xylosandrus ambrosia beetles.

Xylosandrus compactus and X. crassiusculus are two polyphagous ambrosia beetles originating from Asia and invasive in circumtropical regions worldwide. Both species were recently reported in Italy and further invaded several other European countries in the following years. We used the MaxEnt algorithm to estimate the suitable areas worldwide for both species under the current climate. We also made future projections for years 2050 and 2070 using 11 different General Circulation Models, for 4 Representative Concentration Pathways (2.6, 4.5, 6.0 and 8.5). Our analyses showed that X. compactus has not been reported in all potentially suitable areas yet. Its current distribution in Europe is localised, whereas our results predicted that most of the periphery of the Mediterranean Sea and most of the Atlantic coast of France could be suitable. Outside Europe, our results also predicted Central America, all islands in Southeast Asia and some Oceanian coasts as suitable. Even though our results when modelling its potential distribution under future climates were more variable, the models predicted an increase in suitability poleward and more uncertainty in the circumtropical regions. For X. crassiusculus, the same method only yielded poor results, and the models thus could not be used for predictions. We discuss here these results and propose advice about risk prevention and invasion management of both 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.

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/).

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/).

Deciphering the demographic history of allochronic differentiation in the pine processionary moth Thaumetopoea pityocampa.

Understanding the processes of adaptive divergence, which may ultimately lead to speciation, is a major question in evolutionary biology. Allochronic differentiation refers to a particular situation where gene flow is primarily impeded by temporal isolation between early and late reproducers. This process has been suggested to occur in a large array of organisms, even though it is still overlooked in the literature. We here focused on a well-documented case of incipient allochronic speciation in the winter pine processionary moth Thaumetopoea pityocampa. This species typically reproduces in summer and larval development occurs throughout autumn and winter. A unique, phenologically shifted population (SP) was discovered in 1997 in Portugal. It was proved to be strongly differentiated from the sympatric "winter population" (WP), but its evolutionary history could only now be explored. We took advantage of the recent assembly of a draft genome and of the development of pan-genomic RAD-seq markers to decipher the demographic history of the differentiating populations and develop genome scans of adaptive differentiation. We showed that the SP diverged relatively recently, that is, few hundred years ago, and went through two successive bottlenecks followed by population size expansions, while the sympatric WP is currently experiencing a population decline. We identified outlier SNPs that were mapped onto the genome, but none were associated with the phenological shift or with subsequent adaptations. The strong genetic drift that occurred along the SP lineage certainly challenged our capacity to reveal functionally important loci.

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.

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.

Inferences on population history of a seed chalcid wasp: invasion success despite a severe founder effect from an unexpected source population.

Most invasive species established in Europe originate from either Asia or North America, but little is currently known about the potential of the Anatolian Peninsula (Asia Minor) and/or the Near East to constitute invasion sources. Mediterranean forests are generally fragile ecosystems that can be threatened by invasive organisms coming from different regions of the Mediterranean Basin, but for which historical data are difficult to gather and the phylogeographic patterns are still poorly understood for most terrestrial organisms. In this study, we characterized the genetic structure of Megastigmus schimitscheki, an invasive seed-feeding insect species originating from the Near East, and elucidated its invasion route in South-eastern France in the mid 1990s. To disentangle the evolutionary history of this introduction, we gathered samples from the main native regions (Taurus Mountains in Turkey, Lebanon and Cyprus) and from the invaded region that we genotyped using five microsatellite markers and for which we sequenced the mitochondrial Cytochrome Oxidase I gene. We applied a set of population genetic statistics and methods, including approximate Bayesian computation. We proposed a detailed phylogeographic pattern for the Near East populations, and we unambiguously showed that the French invasive populations originated from Cyprus, although the available historical data strongly suggested that Turkey could be the most plausible source area. Interestingly, we could show that the introduced populations were founded from an extremely restricted number of individuals that realized a host switch from Cedrus brevifolia to C. atlantica. Evolutionary hypotheses are discussed to account for this unlikely scenario.

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.

Temperature niche shift observed in a Lepidoptera population under allochronic divergence.

A process of adaptive divergence for tolerance to high temperatures was identified using a rare model system, consisting of two sympatric populations of a Lepidoptera (Thaumetopoea pityocampa) with different life cycle timings, a 'mutant' population with summer larval development, Leiria SP, and the founder natural population, having winter larval development, Leiria WP. A third, allopatric population (Bordeaux WP) was also studied. First and second instar larvae were experimentally exposed to daily-cycles of heat treatment reaching maximum values of 36, 38, 40 and 42 °C; control groups placed at 25 °C. A lethal temperature effect was only significant at 42 °C, for Leiria SP, whereas all temperatures tested had a significant negative effect upon Leiria WP, thus indicating an upper threshold of survival c.a. 6 °C above that of the WP. Cox regression model, for pooled heat treatments, predicted mortality hazard to increase for Leiria WP (+108%) and Bordeaux WP (+78%) in contrast to Leiria SP; to increase by 24% for each additional °C; and to decrease by 53% from first to second instar larvae. High variability among individuals was observed, a population characteristic that may favour selection and consequent adaptation. Present findings provide an example of ecological differentiation, following a process of allochronic divergence. Results further contribute to a better understanding of the implications of climate change for ecological genetics.

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.

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.

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.

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.

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.