Next Article in Journal
Histological Observation of Primary and Secondary Aerenchyma Formation in Adventitious Roots of Syzygium kunstleri (King) Bahadur and R.C.Gaur Grown in Hypoxic Medium
Previous Article in Journal
Impact Assessment of Skidding Extraction: Effects on Physical and Chemical Properties of Forest Soils and on Maple Seedling Growing along the Skid Trail
Previous Article in Special Issue
Infection Levels of the Microsporidium Larssoniella duplicati in Populations of the Invasive Bark Beetle Ips duplicatus: From Native to New Outbreak Areas
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Shallow Genetic Structure among the European Populations of the Six-Toothed Bark Beetle Ips sexdentatus (Coleoptera, Curculionidae, Scolytinae)

1
Forest Research Institute, Hellenic Agricultural Organization Demeter, Vassilika, 57006 Thessaloniki, Greece
2
Institute of Silviculture and Forest Protection, University of Sopron, H-9400 Sopron, Hungary
3
Department of Ecology, University of Alicante, 3690 San Vicente del Raspeig, Spain
4
Sanidad Agricola Econex S.L., 30149 Siscar–Santomera, Murcia, Spain
5
Croatian Forest Research Institute, Cvjetno naselje 41, 10450 Jastrebarsko, Croatia
6
Department of Agronomy, Food, Natural Resources, Animals and Environment (DAFNAE), University of Padova, 35122 Padova, Italy
7
Department of Forest and Soil Sciences, Boku, University of Natural Resources and Life Sciences, 1190 Vienna, Austria
*
Author to whom correspondence should be addressed.
Forests 2019, 10(2), 136; https://doi.org/10.3390/f10020136
Submission received: 3 January 2019 / Revised: 5 February 2019 / Accepted: 6 February 2019 / Published: 7 February 2019

Abstract

:
The six-toothed bark beetle, Ips sexdentatus, is one of the most abundant scolytid species of the central and southern European countries. It mostly feeds on Pinus sp., whereas during population outbreaks it can also attack Picea sp. In spite of its broad distribution, its phylogeography has never been studied before. To do that, we employed an mtDNA marker on 489 individuals that covered most of its native range in Europe. Geographic distribution of the 86 haplotypes showed that at least three glacial refugia have played a significant role in shaping the currently observed pattern of genetic divergence in Europe, without excluding the contribution of minor refugial areas that acted in a similar manner. The revealed shallow structure can be considered an artifact of factors that reduced intraspecific diversity, at the same time favoring gene flow. As such, biological traits of the species itself (flying ability and host preference) and even human-mediated transport of wood seem to be the most prevailing and probable reasons that gave rise to the observed pattern.

1. Introduction

Comparative phylogeographic studies allow the indication of the common refugial areas that species used during ice ages and the migration routes they followed in the postglacial recolonization [1,2]. For Europe, glacial refugia were located in the southeastern parts of Europe where climatic conditions were milder but recently cryptic glacial refugia within the northern parts of Europe were also revealed [3,4,5]. Consequently, the contemporary patterns of diversity observed for many organisms emerged through the amalgamation of lineages that diverged both in major and cryptic refugia [6,7,8,9]. Forest insects are of particular interest as phylogeographic analyses can be used to understand how climatic changes affected their genetic structure and migration behavior in the past. After doing that, we can then assess their behavior within the context of climate change and even under the influence of human-mediated migration, as both these factors shape species distribution today [10].
The Curculionidae family includes the subfamily Scolytinae, including some of the most notorious forest pests affecting the environment in many ways [11,12,13]. Among those, Dendroctonus ponderosae H. for North America and Ips typographus L. for Europe are responsible for enormous economic loss and tremendous ecological damage [14,15,16]. Numerous studies have attempted to decipher the biology, behavior and distribution of these species, even employing phylogeographic approaches [17,18]. For example, the studies on I. typographus [17,19] facilitated the accurate delimitation of the postglacial colonization routes and possible glacial refugial areas in Europe, knowledge that was later further enhanced through investigations on other Scolytinae species [20,21,22,23]. However, a species that is relatively less investigated despite its broad distribution is the six-toothed bark beetle, Ips sexdentatus B. This bark beetle feeds predominantly on pine trees (Pinus spp.) but during outbreaks may attack even spruce (Picea spp.) [24,25], which can be found in most of the Palearctic regions [26], although in Europe is spread mostly in central and southern countries [27]. Even though I. sexdentatus is generally not considered an aggressive species, this scolytid may cause strong damage in weakened trees particularly after fire or in association with other more aggressive species such as Tomicus spp. or other Ips spp. [28,29,30]; at high population levels it can even attack and kill healthy trees [31,32].
As previous studies on European bark beetles have shown, glacial refugia have left a detectable imprint on their genetic diversity pattern, either in a subtler [19,23] or more emphatic manner [33]. The rate of influence is largely determined by the biological traits in concert with the duration of evolutionary history of the species themselves [19]. Here we present the phylogeography of I. sexdentatus of European populations based on an mtDNA marker, with the aim to assess the influence of glacial refugia in the observed pattern of intraspecific diversity within its native range.

2. Materials and Methods

2.1. Sampling—DNA Analysis

From 2011 to 2015, Ips sexdentatus adults were collected by pheromone traps and shipped to the Forest Research Institute (Hellenic Agricultural Organization Demeter, Thessaloniki, Greece). In total 50 localities from 17 countries were sampled, covering almost the whole European distribution of the species (Supplementary Table S1). Specimens were sent in absolute ethanol and stored at −20 °C in the laboratory. DNA of 489 individuals was extracted using the PureLink® Genomic DNA kit (Invitrogen, Carlsbad, CA, USA) following the manufacturer’s instructions, with the only modification that the initial tissue grinding was performed using stainless steel grinding balls. Polymerase Chain Reaction (PCR) was run in 25 μL volumes, with primers 2183 and 3041 targeting an 850 bp-fragment of mtDNA’s cytochrome oxidase 1 gene [34]. Each reaction contained 8 μL of the extracted DNA, 10.6 μL of double distilled water, 5 μL of Red Bioline buffer (provided with the Taq), 0.5 μL of each primer and 0.4 μL of MyTaq (Red Bioline), adding up to a total of 25 μL. PCR procedure was as follows: 3 min at 94 °C, followed by 40 cycles of 30 s at 94 °C (denaturation), 30 s at 45 °C (annealing) and 1.5 min at 72 °C (extension). The final extension period was carried out at 72 °C over the course of 7 min. Purification of PCR products was performed using the PureLink® PCR Purification Kit (Invitrogen, Carlsbad, CA, USA) following the manufacturer’s protocol. Sequencing took place at CEMIA SA (Larissa, Greece) using an ABI 3730XL.

2.2. Data Analysis

Chromas Lite® was used to visualize the sequences which then were aligned and checked using Clustal X (version 2) [35]. Haplotypes that were represented by a single individual were verified by additional sequencing of an independent amplicon, in order to exclude cases of base misincorporation as a result of a PCR error [36]. After this double check, haplotypes were described and deposited in National Center for Biotechnology Information (NCBI) GenBank.
Patterns of molecular diversity were assessed by estimating haplotype (Hd), nucleotide diversity (p) [37] and the average number of nucleotide differences (k) [38] for every population using MEGA v.6 [39]. The ratio of haplotypes/individuals was then compared with previous mtDNA studies on bark beetles [40]. Bayesian phylogeny was constructed with MrBayes version 3.1.1 [41] using the nucleotide substitution model that best fits the data under the hierarchical likelihood–ratio test (hLRT) and Akaike criteria determined by MrModeltest version 2.1 [42]. The most appropriate model was found to be Hasegawa-Kishino-Yano with invariable sites and gamma distributed heterogeneity (HKY + I + G) [43], subsequently integrated in the Bayesian analysis. The initial run lasted 2 × 106 generations, with a sampling frequency of 100 generations in order to determine the number of trees to be discarded in the final run. Calculation of standard deviation showed that after 5 × 105 generations a plateau has been reached, thus the first 500 trees were discarded.
Population dynamics were assessed using Tajima’s D [44], Fu and Li’s F [45] statistics and a mismatch distribution analysis which were calculated by MEGA 6 [39], whereas the statistical parsimony network was estimated using Automated Nested Clade Analysis (ANECA) [46] that relies on TCS [47] and GeoDis [48] to infer the demographic processes that produced the observed pattern of divergence.
Clustering of populations was initially assessed with Spatial Analysis of Molecular Variance (SAMOVA 1.0) [49], to unveil the geographic clustering pattern and the minimum number of groups that produces FCT (among-groups variation) values higher than FST (among-populations/within groups variation), both with and without taking into account the coordinates of populations. This estimation would allow identification of the point at which variation is mostly shaped by geographic groups rather than single populations. As the impact of isolation by distance [50] is very often heavily imprinted on intraspecific diversity [51], the relationship between genetic and geographic distances was assessed by the Mantel test [52] as this is implemented in Alleles in Space (AIS) [53] with 1000 permutations to assess the statistical significance of the coefficient. AIS was also used to visualize the spatial pattern of genetic diversity by running a Landscape Shape Interpolation (LSI) analysis (distance weighting parameter a = 1/50 × 50 grid). The 3-d surface plot produced exhibited the geographic coordinates (X- and Y-axis) and the genetic distances (Z-axis), with peaks representing areas of high genetic distance among individuals (restricted gene flow) and troughs point out areas of low genetic distance (unimpeded gene flow). Visualization of species diversity distribution was also performed using the Genetic Landscape Toolbox [54] under ArcMap 9.3 (ESRI, Redlands, CA, USA) according to author’s instructions. In both these approaches to visualize the distribution of genetic divergence, populations with few individuals (less than 3) were excluded.

3. Results

Analysis of the 674-long region of the mitochondrial Cytochrome Oxidase One (COI) gene of 489 I. sexdentatus individuals yielded 86 haplotypes. The mean number of substitutions was 1.9 nucleotides (=0.28%) whereas the maximum was 15 nucleotides (=2.2%). Among those 86 haplotypes, seven had more than ten individuals while 48 haplotypes contained a single individual. No amino acid change could be found translating the DNA of the haplotypes (Supplementary Table S1). Previous genetic studies on bark beetles [40] showed that between the number of individuals analyzed and the respective haplotypes retrieved there is a positive correlation that is statistically significant (Pearson’s correlation coefficient r = 0.6181, p < 0.05). Integrating the number of haplotypes of the current investigation produced a similar coefficient and statistical significance (r = 0.6276, p < 0.05) (Supplementary Figure S1).
A Bayesian tree was calculated (Figure 1) and only three clusters supported by high (>85%) posterior probabilities indicated by colors could be identified. These clusters were mapped onto a European map (Figure 2). The green clade with six haplotypes only included individuals from Italy and the blue clade with two haplotypes that included only individuals from Greece. The red clade encompassed samples from several central–north eastern countries (Germany, Austria, Slovenia, Hungary and Ukraine) (Figure 2). All other haplotypes were marked with pale red and thus represented the biggest parts of the pie charts. HT25 was the widest distributed haplotype, occurring in the majority of populations. HT20 was found in populations occurring in the southeastern/central parts of Europe, with its mutationally related haplotypes (HT17, HT18, HT19, HT21, HT51, HT52 and HT53) being geographically limited in southeast Europe. Something similar was observed with HT02, which was also found in southeastern/central populations, whereas the closest related haplotypes are found in the Balkans (Supplementary Figure S2).
In terms of population dynamics, both neutrality indices were negative (D = −2.30, p < 0.01/F = −5.09, p < 0.02), showing an excess of rare alleles that can be observed in a population expanding after a recent bottleneck event. This was also verified by the unimodal distribution shown by the mismatched distribution (Supplementary Figure S2). Automated Nested Clade Analysis (ANECA) did not reach a conclusion for the total cladogram, likely a consequence of the shallow structure among populations; for the higher level clusters that corresponded to the identified clades (2-3, 3-1, 4-1), the conclusion was restricted gene flow with isolation-by-distance, suggesting that certain areas became isolated, favoring the genetic divergence of individuals therein (Supplementary Figure S3).
For both options with and without coordinates, the number of clusters that produced an FCT value (=0.39473) higher than FST (=0.38497) is achieved only when they were separated into many groups (N = 15); however, 2/3 of these groups contained a single population, mainly located in southeastern Europe (data not shown). A similar image was concluded by Alleles in Space that allows the spatial visualization of the areas that exhibit high (peaks) and low (troughs) diversity. Southern populations exhibited overall higher diversity than central and northern ones, with numerous peaks located in the southern longitudes (Supplementary Figure S4). ArcMap 9.1 visualized these data where red color (points of high haplotype diversity) was mostly found in Italy, the northern parts of the Iberian Peninsula and the Balkans, indicating the effect of multiple areas that presumably might have acted as sources of high diversity (Supplementary Figure S5).

4. Discussion

The first genetic analysis of European I. sexdentatus populations identified 86 haplotypes analyzing 489 individuals from 17 countries, that covered most of the natural range of this species in Europe. This haplotype diversity falls within the range of previous phylogenetic investigations on bark beetle species worldwide [40] (Supplementary Figure S1). Based on the statistically supported clades observed in the topology of the phylogenetic tree, it can be deduced that at least three refugial areas have shaped the genetic structure of I. sexdentatus: one located on the Italian Peninsula, one located in the southern part of the Balkans and one likely in the Dinaric region. The distinct clustering and separation of haplotypes originating from these regions suggest that a fraction of diversity was maintained there at least during the latest glacial periods, a pattern similar to the ones retrieved for other conifer-feeding bark beetles in Europe like Tomicus piniperda L. [22,23], Ips typographus [8,55] and Pityogenes chalcographus L. [8,33,56]. However, the occurrence of additional haplotype clusters with common geographic origins (clusters 2-2 and 2-3; Supplementary Figure S3) argues for a possible impact of additional minor refugia beside the above mentioned major refugial areas on the genetic divergence among the European I. sexdentatus populations. This finding is congruent with recent investigations [6] that support the concept of “refugia within refugia” [57]. Southeastern parts of Europe in particular seem to have been such a region, with the diverse landscape favoring the emergence of intraspecific divergence [58], something that was visualized by AIS and ArcMap (Supplementary Figures S4 and S5, respectively). Both approaches demonstrated not only that southern populations outperformed central and northern ones in terms of diversity (Supplementary Figure S4), but also that multiple areas acted as diversity hot spots for I. sexdentatus populations (Supplementary Figure S5).
However, this pattern seems to be blurred by factors that acted upon the initial differentiation. Such a shallow genetic structure may arise due to recent lineage divergence [59] or the life history traits of the species itself [60]. The latter particularly seems to have played an important role for I. sexdentatus; as it has been recently observed, the flight ability of beetles largely determines their diversification rate, with flightless species being in favor of allopatric speciation not only by retaining higher genetic differentiation among distant populations but also demonstrating a higher number of genetically distinct lineages than flight-capable species [61]. With the flight capability of I. sexdentatus ranging between 5 km (98% of the samples) and 45 km (10%) [62], gene flow among populations cannot be overlooked. It thus seems possible that gene flow effectively reduced most of the divergence that could derive from different refugial areas resulting in this shallow pattern of intraspecific differentiation [2,63]. Moreover, in line with the hypothesis expressed previously [64], host specialization of the bark beetle species strongly determines the degree of intraspecific genetic diversity. Specialist species being more prone to higher divergence than generalist species as host selection cannot exert any kind of selection that could facilitate the emergence of differentiation for the latter. As a generalist, I. sexdentatus can attack spruce (Picea abies (L.) Karst) beside its main host (Pinus sp.), particularly during mass outbreaks, and thus the observed low intraspecific divergence could be attributed to this trait. To them, another factor should also be added: The unimpeded trade of wood and wood packaging material among the European countries. Recent studies [65] have shown that I. sexdentatus is the second most common scolytid species in international ports and in the associated wood waste landfills. As expected, human-mediated movement of individuals over long distances smoothened the genetic landscape, conforming populations to panmixia as has been already demonstrated [66,67].

5. Conclusions

This study verified the impact of at least three glacial refugia on the European populations of I. sexdentatus, with indications of additional minor refugial areas that contributed to the currently observed pattern of diversity. Nevertheless, the shallow structure evinced the effect of biotic (biological traits of the species) and abiotic (human-mediated transport) factors that smoothened divergence. Taking this into account, future studies should aim at a refugia-oriented sampling coupled with a comprehensive set of genome-wide molecular markers using ddRADSeq [33], something that could considerably increase the resolution capacity of the approach, allowing the conclusion even of subtle structuring, as is the case for I. sexdentatus.

Supplementary Materials

The following are available online at https://www.mdpi.com/1999-4907/10/2/136/s1, Supplementary Figure S1. Correlation of haplotypes retrieved against the number of individuals analyzed in bark beetle species phylogenetic studies. Supplementary Figure S2. Mismatch distribution diagram inferred from the mitochondrial DNA sequences of I. sexdentatus. Supplementary Figure S3. Cladogram with respective conclusions inferred with ANECA. Colored clusters represent statistically supported (>85% posterior probability) clades of the phylogenetic tree (colors identical with Figure 1). Supplementary Figure S4. Landscape of genetic diversity created with Alleles in Space; peaks represent areas of high genetic distance among individuals and troughs point out areas of low genetic distance. Supplementary Figure S5. Distribution of haplotype diversity (based on Tamura-Nei distance) over the populations sampled. Red color indicates points of high and blue color shows areas of low diversity. Supplementary Table S1. Individuals analyzed per population, and their haplotype assignment.

Author Contributions

D.N.A., M.F., C.S. and F.L. conceived the idea; D.N.A. conducted the experiments; R.W., M.P., D.G., F.L. and M.F. organized the thorough sampling; D.N.A., F.L. and C.S. analyzed the data. All authors contributed to the writing.

Funding

Austrian Science Fund FWF (project number P26749-B25).

Acknowledgments

We thank the Austrian Science Fund FWF for financial support of this study (project number P26749-B25).

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Taberlet, P.; Fumagalli, L.; Wust-Saucy, A.-G.; Cosson, J.-F. Comparative phylogeography and postglacial colonization routes in Europe. Mol. Ecol. 1998, 7, 453–464. [Google Scholar] [CrossRef] [PubMed]
  2. Hewitt, G.M. The genetic legacy of the Quaternary ice ages. Nature 2000, 405, 907–913. [Google Scholar] [CrossRef]
  3. Stewart, J.R.; Lister, A.M.; Barnes, I.; Dalén, L. Refugia revisited: Individualistic responses of species in space and time. Proc. R. Soc. Lond. B Biol. Sci. 2010, 277, 661–671. [Google Scholar] [CrossRef] [PubMed]
  4. Tzedakis, P.C.; Emerson, B.C.; Hewitt, G.M. Cryptic or mystic? Glacial tree refugia in northern Europe. Trends Ecol. Evol. 2013, 28, 696–704. [Google Scholar] [CrossRef] [PubMed]
  5. Kajtoch, Ł.; Cieślak, E.; Varga, Z.; Paul, W.; Mazur, M.A.; Sramkó, G.; Kubisz, D. Phylogeographic patterns of steppe species in Eastern Central Europe: A review and the implications for conservation. Biodivers. Conserv. 2016, 25, 2309–2339. [Google Scholar] [CrossRef]
  6. Médail, F.; Diadema, K. Glacial refugia influence plant diversity patterns in the Mediterranean Basin. J. Biogeogr. 2009, 36, 1333–1345. [Google Scholar] [CrossRef]
  7. Theissinger, K.; Balint, M.; Feldheim, K.A.; Haase, P.; Johannesen, J.; Laube, I.; Pauls, S.U. Glacial survival and post-glacial recolonization of an arctic-alpine freshwater insect (Arcynopteryx dichroa, Plecoptera, Perlodidae) in Europe. J. Biogeogr. 2012, 40, 236–248. [Google Scholar] [CrossRef]
  8. Bertheau, C.; Bankhead-Dronnet, S.; Martin, C.; Lieutier, F.; Roux-Morabito, G. Lack of genetic differentiation after host range extension argues for the generalist nature of Pityogenes chalcographus (Curculionidae: Scolytinae). Ann. For. Sci. 2012, 69, 313–323. [Google Scholar] [CrossRef]
  9. Daneck, H.; Fer, T.; Marhold, K. Glacial survival in northern refugia? Phylogeography of the temperate shrub Rosa pendulina L. (Rosaceae): AFLP vs. chloroplast DNA variation. Biol. J. Linn. Soc. 2016, 119, 704–718. [Google Scholar] [CrossRef]
  10. Rousselet, J.; Zhao, R.; Argal, D.; Simonato, M.; Battisti, A.; Roques, A.; Kerdelhue, C. The role of topography in structuring the demographic history of the pine processionary moth, Thaumetopoea pityocampa (Lepidoptera, Notodontidae). J. Biogeogr. 2010, 37, 1478–1490. [Google Scholar] [CrossRef]
  11. Brockerhoff, E.G.; Bain, J.; Kimberley, M.; Knižek, M. Interception frequency of exotic ark and ambrosia beetles (Coleoptera: Scolytinae) and relationship with establishment in New Zealand and worldwide. Can. J. For. Res. 2006, 36, 289–298. [Google Scholar] [CrossRef]
  12. Fettig, C.J.; Klepzig, K.D.; Billings, R.F.; Munson, A.S.; Nebeker, T.E.; Negrón, J.F.; Nowak, J.T. The effectiveness of vegetation management practices for prevention and control of bark beetle infestations in coniferous forests of the western and southern United States. For. Ecol. Manag. 2007, 238, 24–53. [Google Scholar] [CrossRef]
  13. Abdulah, H.; Darvishzadeh, R.; Skidmore, A.K.; Groen, T.A.; Heurich, M. European spruce bark beetle (Ips typographus L.) green attack affects foliar reflectance and biochemical properties. Int. J. Appl. Earth Obs. Geoinf. 2018, 64, 199–209. [Google Scholar] [CrossRef]
  14. Kurz, W.A.; Dymond, C.C.; Stinson, G. Mountain pine beetle and forest carbon: Feedback to climate change. Nature 2008, 454, 987–990. [Google Scholar] [CrossRef] [PubMed]
  15. Waring, K.M.; Reboletti, D.M.; Mork, L.A.; Huang, C.-H.; Hofstetter, R.W.; Garcia, A.M.; Fulé, P.Z.; Davis, T.S. Modeling the impacts of two bark beetle species under a warming climate in the southwestern USA: Ecological and economic consequences. Environ. Manag. 2009, 44, 824–835. [Google Scholar] [CrossRef] [PubMed]
  16. Vega, F.E.; Hofstetter, R.W. Bark Beetles: Biology and Ecology of Native and Invasive Species, 1st ed.; Elsevier Academic Press: Oxford, UK, 2015. [Google Scholar]
  17. Sallé, A.; Arthofer, W.; Lieutier, F.; Stauffer, C.; Kerdelhué, C. Phylogeography of a host-specific insect: Genetic structure of Ips typographus in Europe does not reflect past fragmentation of its host. Biol. J. Linn. Soc. 2007, 90, 239–246. [Google Scholar] [CrossRef]
  18. Cullingham, C.I.; Roe, A.D.; Sperling, F.A.H.; Coltman, D.W. Phylogeographic insights into an irruptive pest outbreak. Ecol. Evol. 2012, 2, 908–919. [Google Scholar] [CrossRef] [PubMed]
  19. Bertheau, C.; Schuler, H.; Arthofer, W.; Avtzis, D.N.; Mayer, F.; Krumböck, S.; Moodley, Y.; Stauffer, C. Divergent evolutionary histories of two sympatric spruce bark beetle species. Mol. Ecol. 2013, 22, 3318–3332. [Google Scholar] [CrossRef] [PubMed]
  20. Kerdelhué, C.; Roux-Morabito, G.; Forichon, J.; Chambon, J.-M.; Robert, A.; Lieutier, F. Population genetic structure of Tomicus piniperda L. (Curculionidae: Scolytinae) on different pines species and validation of T. destruens. Mol. Ecol. 2002, 11, 483–494. [Google Scholar]
  21. Gallego, D.; Cánovas, F.; Esteve, M.A.; Galián, J. Descriptive biogeography of Tomicus (Coleoptera: Scolytidae) species in Spain. J. Biogeogr. 2004, 31, 2011–2014. [Google Scholar] [CrossRef]
  22. Faccoli, M.; Piscedda, A.; Salvato, P.; Simonato, M.; Masutti, L.; Battisti, A. Genetic structure and phylogeography of pine shoot beetle populations (Tomicus destruens and T. piniperda, Coleoptera, Scolytidae) in Italy. Ann. For. Sci. 2005, 62, 361–368. [Google Scholar] [CrossRef]
  23. Horn, A.; Stauffer, C.; Lieutier, F.; Kerdelhué, C. Complex postglacial history of the temperate bark beetle Tomicus piniperda L. (Coleoptera, Scolytinae). Heredity 2009, 103, 238–247. [Google Scholar] [CrossRef] [PubMed]
  24. Rener, I.; Maja, J. The six-toothed bark beetle Ips sexdentatus (Boerner, 1767) (Coleoptera: Scolytidae) in the fire site Kojnik. J. For. Sci. 2001, 2, 154–155. [Google Scholar]
  25. Ozcan, G.E.; Eroglu, M.; Akinci, H.A. Use of pheromone-baited traps for monitoring Ips sexdentatus (Boerner) (Coleoptera: Curculionidae) in oriental spruce stands. Afr. J. Biotechnol. 2011, 10, 16351–16360. [Google Scholar]
  26. Pfeffer, A. Zentral- und westpaläarktische Borken- und Kernkäfer (Coleoptera: Scolytidae, Platypodidae); Pro Entomologica, Naturhistorisches Museum: Basel, Switzerland, 1995. [Google Scholar]
  27. Alonso-Zarazaga, M.; Barrios, H.; Borovec, R.; Bouchard, P.; Caldara, R.; Colonnelli, E.; Gültekin, L.; Hlaváč, P.; Korotyaev, B.; Lyal, C.; et al. Cooperative Catalogue of Palaearctic Coleoptera Curculionoidea; Sociedad Entomológica Aragonesa (S.E.A.): Zaragoza, Spain, 2017. [Google Scholar]
  28. Paiva, M.R.; Pessoa, M.F.; Vité, J.P. Reduction in the pheromone attractant response of Orthotomicus erosus (Woll.) and Ips sexdentatus Boern. (Col., Scolytidae). J. Appl. Entom. 1988, 106, 198–200. [Google Scholar] [CrossRef]
  29. Lombardero, M.J.; Ayres, M.P. Factors influencing bark beetle outbreaks after forest fires on the Iberian Peninsula. Environ. Entomol. 2011, 40, 1007–1018. [Google Scholar] [CrossRef] [PubMed]
  30. Faccoli, M. European Bark and Ambrosia Beetles: Types, Characteristics and Identification of Mating Systems; World Biodiversity Association Handbooks 5: Verona, Italy, 2015. [Google Scholar]
  31. Rossi, J.P.; Samalens, J.-C.; Guyon, D.; van Halder, I.; Jactel, H.; Menassieu, P.; Piou, D. Multiscale spatial variation of the bark beetle Ips sexdentatus damage in a pine plantation forest (Landers de Gascogne, southwestern France). For. Ecol. Manag. 2009, 257, 1551–1557. [Google Scholar] [CrossRef]
  32. Pineau, X.; Bourguignon, M.; Jactel, H.; Lieutier, F.; Sallé, A. Pyrrhic victory for bark beetles: Successful standing tree colonization triggers strong intraspecific competition for offspring of Ips sexdentatus. For. Ecol. Manag. 2017, 399, 188–196. [Google Scholar] [CrossRef]
  33. Schebeck, M.; Dowle, E.J.; Schuler, H.; Avtzis, D.N.; Bertheau, C.; Feder, J.L.; Ragland, G.J.; Stauffer, C. Pleistocene climate cycling and host plant association shaped the demographic history of the bark beetle Pityogenes chalcographus. Sci. Rep. 2018, 8, 14207. [Google Scholar] [CrossRef] [PubMed]
  34. Simon, C.; Frati, F.; Beckenbach, A.; Crespi, B.; Liu, H.; Flook, P. Evolution, weighting, and phylogenetic utility of mitochondrial gene sequences and a compilation of conserved polymerase chain reaction primers. Ann. Entomol. Soc. Am. 1994, 87, 651–701. [Google Scholar] [CrossRef]
  35. Thompson, J.D.; Gibson, T.J.; Plewniak, F.; Jeanmougin, F.; Higgind, D.G. The CLUSTAL_X windows interface: Flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res. 1997, 25, 4876–4882. [Google Scholar] [CrossRef] [PubMed]
  36. Kobayashi, N.; Tamura, K.; Aotsuka, T. PCR error and molecular population genetics. Biochem. Gen. 1999, 37, 317–321. [Google Scholar] [CrossRef]
  37. Nei, M. Molecular Evolutionary Genetics; Columbia University Press: New York, NY, USA, 1987. [Google Scholar]
  38. Tajima, F. Evolutionary relationship of DNA sequences in finite populations. Genetics 1983, 105, 437–460. [Google Scholar] [PubMed]
  39. Tamura, K.; Stecher, G.; Peterson, D.; Filipski, A.; Kumar, S. MEGA6: Molecular evolutionary genetics analysis version 6.0. Mol. Biol. Evol. 2013, 30, 2725–2729. [Google Scholar] [CrossRef] [PubMed]
  40. Avtzis, D.N.; Bertheau, C.; Stauffer, C. What is next in bark beetle phylogeography? Insects 2012, 3, 453–472. [Google Scholar] [CrossRef] [PubMed]
  41. Ronquist, F.; Huelsenbeck, J. MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics 2003, 19, 1572–1574. [Google Scholar] [CrossRef] [PubMed]
  42. Nylander, J.; Ronquist, F.; Huelsenbeck, J.; Nieves-Aldrey, J. Bayesian phylogenetic analysis of combined data. Syst. Biol. 2004, 53, 47–67. [Google Scholar] [CrossRef]
  43. Hasegawa, M.; Kishino, K.; Yano, T. Dating the human-ape splitting by a molecular clock of mitochondrial DNA. J. Mol. Evol. 1985, 22, 160–174. [Google Scholar] [CrossRef]
  44. Tajima, F. Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics 1989, 123, 585–595. [Google Scholar]
  45. Fu, Y.-X.; Li, W.H. Statistical test of neutrality of mutations. Genetics 1993, 133, 693–709. [Google Scholar]
  46. Panchal, M. The automation of Nested Clade Phylogeographic Analysis. Bioinformatics 2007, 23, 509–510. [Google Scholar] [CrossRef] [PubMed]
  47. Clement, M.; Posada, D.; Crandall, K.A. TCS: A computer program to estimate gene genealogies. Mol. Ecol. 2000, 9, 1657–1659. [Google Scholar] [CrossRef] [PubMed]
  48. Posada, D.; Crandall, K.A.; Templeton, A.R. GeoDis: A program for the cladistic nested analysis of the geographical distribution of genetic haplotypes. Mol. Ecol. 2000, 9, 487–488. [Google Scholar] [CrossRef]
  49. Dupanloup, I.; Schneider, S.; Excoffier, L. A simulated annealing approach to define the genetic structure of populations. Mol. Ecol. 2002, 11, 2571–2581. [Google Scholar] [CrossRef] [PubMed]
  50. Wright, S. Isolation by distance. Genetics 1943, 28, 114–138. [Google Scholar] [PubMed]
  51. Safner, T.; Miller, M.P.; McRae, B.H.; Fortin, M.J.; Manel, S. Comparison of Bayesian clustering and edge detection methods for inferring boundaries in landscape genetics. Int. J. Mol. Sci. 2011, 12, 865–889. [Google Scholar] [CrossRef] [PubMed]
  52. Mantel, N. The detection of disease clustering and a generalized regression approach. Cancer Res. 1967, 27, 209–220. [Google Scholar] [PubMed]
  53. Miller, M.P. Alleles in space: Computer software for the joint analysis of interindividual spatial and genetic information. J. Hered. 2005, 96, 722–724. [Google Scholar] [CrossRef]
  54. Vandergast, A.G.; Perry, W.M.; Lugo, R.V.; Hathaway, S.A. Genetic landscapes GIS Toolbox: Tools to map patterns of genetic divergence and diversity. Mol. Ecol. Resour. 2011, 11, 158–161. [Google Scholar] [CrossRef]
  55. Stauffer, C.; Lakatos, F.; Hewitt, G.M. Phylogeography and postglacial colonization routes of Ips typographus (Coleoptera, Scolytidae). Mol. Ecol. 1999, 8, 763–773. [Google Scholar] [CrossRef]
  56. Avtzis, D.N.; Arthofer, W.; Stauffer, C. Sympatric occurrence of diverged mtDNA lineages of Pityogenes chalcographus (Coleoptera, Scolytinae) in Europe. Biol. J. Linn. Soc. 2008, 94, 331–340. [Google Scholar] [CrossRef]
  57. Gómez, A.; Lunt, D.H. Refugia within Refugia: Patterns of Phylogeographic Concordance in the Iberian Peninsula. In Phylogeography of Southern European Refugia; Weiss, S., Ferrand, N., Eds.; Springer: Dordrecht, The Netherland, 2007; pp. 155–188. [Google Scholar]
  58. Poulakakis, N.; Kapli, P.; Lymberakis, P.; Trichas, A.; Vardinoyiannis, K.; Sfenthourakis, S.; Mylonas, M. A review of phylogeographic analyses of animal taxa from the Aegean and surrounding regions. J. Zool. Syst. Evol. Res. 2014, 53, 18–32. [Google Scholar] [CrossRef]
  59. Fontanella, F.M. Mitochondrial DNA panmixia in the little brown seastar Echinaster spinulosus suggests a recent population expansion. DNA Seq. 2016, 27, 1–6. [Google Scholar] [CrossRef] [PubMed]
  60. Neethling, M.; Matthee, C.A.; Bowie, R.C.K.; von der Heyden, S. Evidence of panmixia despite barriers to gene flow in the southern African endemic, Caffrogobius caffer (Teleostei: Gobiidae). BMC Evol. Biol. 2008, 8, 325. [Google Scholar] [CrossRef] [PubMed]
  61. Ikeda, H.; Nishikawa, M.; Sota, T. Loss of flight promotes beetle diversification. Nat. Commun. 2012, 3, 648. [Google Scholar] [CrossRef] [PubMed]
  62. Jactel, H. Individual variability of the flight potential of Ips sexdentatus Boern. (Coleoptera: Scolytidae) in relation to day of emergence, sex, size and lipid content. Can. Entomol. 1993, 125, 919–930. [Google Scholar] [CrossRef]
  63. Schmitt, T. Molecular biogeography of Europe: Pleistocene cycles and postglacial trends. Front. Zool. 2007, 4, 11. [Google Scholar] [CrossRef]
  64. Kelley, S.T.; Farrell, B.D.; Mitton, J.B. Effects of specialization on genetic differentiation in sister species of bark beetles. Heredity 2000, 84, 218–227. [Google Scholar] [CrossRef]
  65. Rassati, D.; Faccoli, M.; Marini, L.; Haack, R.A.; Battisti, A.; Toffolo, E.P. Exploring the role of wood waste landfills in early detection of non-native wood-boring beetles. J. Pest Sci. 2015, 88, 563–572. [Google Scholar] [CrossRef]
  66. Robinet, C.; Imbert, C.-E.; Rousselet, J.; Sauvard, D.; Garcia, J.; Goussard, F.; Roques, A. Human-mediated long-distance jumps of the pine processionary moth in Europe. Biol. Invasions 2012, 14, 1557–1569. [Google Scholar] [CrossRef]
  67. Lacković, N.; Pernek, M.; Bertheau, C.; Franjević, D.; Stauffer, C.; Avtzis, D.N. Limited genetic structure of gypsy moth populations reflecting a recent history in Europe. Insects 2018, 9, 143. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Phylogenetic tree of the 86 haplotypes created with MrBayes version 3.1.1 and inferred using the nucleotide substitution model (HKY + I + G) suggested by MrModeltest version 2.1.
Figure 1. Phylogenetic tree of the 86 haplotypes created with MrBayes version 3.1.1 and inferred using the nucleotide substitution model (HKY + I + G) suggested by MrModeltest version 2.1.
Forests 10 00136 g001
Figure 2. Distribution of mtDNA clades. Colors are similar to Figure 1. Pie charts are proportional to the number of individuals analyzed.
Figure 2. Distribution of mtDNA clades. Colors are similar to Figure 1. Pie charts are proportional to the number of individuals analyzed.
Forests 10 00136 g002

Share and Cite

MDPI and ACS Style

Avtzis, D.N.; Lakatos, F.; Gallego, D.; Pernek, M.; Faccoli, M.; Wegensteiner, R.; Stauffer, C. Shallow Genetic Structure among the European Populations of the Six-Toothed Bark Beetle Ips sexdentatus (Coleoptera, Curculionidae, Scolytinae). Forests 2019, 10, 136. https://doi.org/10.3390/f10020136

AMA Style

Avtzis DN, Lakatos F, Gallego D, Pernek M, Faccoli M, Wegensteiner R, Stauffer C. Shallow Genetic Structure among the European Populations of the Six-Toothed Bark Beetle Ips sexdentatus (Coleoptera, Curculionidae, Scolytinae). Forests. 2019; 10(2):136. https://doi.org/10.3390/f10020136

Chicago/Turabian Style

Avtzis, Dimitrios N., Ferenc Lakatos, Diego Gallego, Milan Pernek, Massimo Faccoli, Rudolf Wegensteiner, and Christian Stauffer. 2019. "Shallow Genetic Structure among the European Populations of the Six-Toothed Bark Beetle Ips sexdentatus (Coleoptera, Curculionidae, Scolytinae)" Forests 10, no. 2: 136. https://doi.org/10.3390/f10020136

APA Style

Avtzis, D. N., Lakatos, F., Gallego, D., Pernek, M., Faccoli, M., Wegensteiner, R., & Stauffer, C. (2019). Shallow Genetic Structure among the European Populations of the Six-Toothed Bark Beetle Ips sexdentatus (Coleoptera, Curculionidae, Scolytinae). Forests, 10(2), 136. https://doi.org/10.3390/f10020136

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop