Globally distributed Xyleborus species reveal recurrent intercontinental dispersal in a landscape of ancient worldwide distributions
© Gohli et al. 2016
Received: 11 December 2015
Accepted: 4 February 2016
Published: 15 February 2016
Invasive species can have devastating effects on native ecosystems and therefore impose a significant threat to human welfare. The introduction rate of invasive species has accelerated dramatically in recent times due to human activity (anthropogenic effects), with a steadily growing pool of widespread tramp species. We present an in-depth analysis of four pantropical species of Xyleborus ambrosia beetles (Xyleborus volvulus, Xyleborus perforans, Xyleborus ferrugineus, and Xyleborus affinis) with similar ecology (fungus cultivation in dead wood), reproductive biology (permanent inbreeding) and genetic system (haplodiploidy). The unique combination of reproductive traits and broad host plant usage pre-adapts these beetles for colonizing of new areas.
We found that all four species were broadly distributed long before human-assisted dispersal became common, and that the impact of anthropogenic effects varied among the species. For X. volvulus, X. perforans, and X. affinis there was evidence of ancient establishment in numerous regions, but also of abundant recent introductions into previously colonized areas. For X. ferrugineus, we found clear biogeographical structuring of old clades, but little evidence for recent successful introductions.
Our results indicate that current human-aided transoceanic dispersal has strongly affected the genetic makeup of three of the species in this study. However, current biogeographical patterns of all four species are equally, if not more strongly, influenced by ancient establishment on different continents.
Invasive species constitute a major threat to native ecosystems worldwide [1, 2]. Many different types of organisms can be successful colonizers, and particularly those easily transported by humans over long distances . Anthropogenic dispersal is therefore regarded as the most significant factor in explaining threats to ecosystems [2, 4] and agricultural crops . However, long-range dispersal is not always human-aided and the relative importance of natural vs human-aided dispersal is often not obvious.
Even though many species are introduced to new continents or oceans every year, only a few become permanently established . Serendipity plays a role in propagule establishment, but flexible habitat preferences and high propagule pressure are certainly important [7, 8]. Certain organisms have biological traits that likely facilitate establishment after dispersal [pre-adapted colonizers; 6, 9, 10]. Among invasive insects, mites, annelids, and plants, species that reproduce by sibling mating (or selfing) or by clonal reproduction through parthenogenesis are overrepresented [11–14]. For instance, the bark and ambrosia beetle fauna (Scolytinae) on remote oceanic islands is generally represented by two or three times more sib-mating species compared to the mainland , and nearly three-quarters of the established alien Scolytinae in North America are brother-sister maters .
With both inbreeding and parthenogenesis, a single mated female is sufficient for a species to become established in a new area [10, 16, 17]. Reproductive assurance is advantageous to insects such as bark beetles that depend on locating scattered, ephemeral resources, and enables long-distance colonization and establishment. For normally outbreeding species, a large number of introduced individuals  or repeated introductions  can ameliorate the effects of inbreeding depression. Permanent inbreeders , on the other hand, are expected to have purged strongly deleterious alleles early in a lineage's history, and both parthenogens and inbreeders have been under strong selection for combinations of genes which work well together. Regularly inbreeding species are thus expected to be largely immune to the effects of inbreeding that result from bottleneck events such as colonisations , and are therefore relatively unaffected by many ecological and genetic mechanisms normally causing Allee effects .
In both the North American and European wood boring fauna, exotic species are disproportionately ambrosia beetles that cultivate and feed upon mutualistic fungi in tunnels they excavate in dead trees [15, 23]. The fungal symbionts are host generalists, which enable the beetles to breed in many host plants families [24, 25], pre-adapting them for successful colonization of distant regions [4, 10–13]. Further, in contrast to the native faunas, exotic ambrosia beetles are almost all inbreeders [26, 27].
The ambrosia beetle genus Xyleborus contains many abundant and widespread invasive species [28–32]. All inbreed by sibling mating and all are haplodiploid (males are haploid and produced by unfertilized eggs). By mating with siblings,having the potential for re-mating with clonally produced sons in the absence of mates, and the ability to grow fungi for food in almost any kind of wood, Xyleborus and similar ambrosia beetles are exceptionally efficient at colonizing and establishing in new areas.
Human activity results in a constant stream of introductions of bark and ambrosia beetles to new areas [26, 33–36]. For wood boring beetles such as these, the most common mechanism for intercontinental dispersal is transport in timber and wooden packing materials [26, 34]. They survive such transport particularly well, being ensconced in material which provides food, a buffered microenvironment, and protection from most natural enemies. Hence, a wide variety of exotic bark and ambrosia beetle species are regularly trapped near harbours around the world [26, 27, 34, 37–39]. While anthropogenic effects strongly influence the recent spread of these beetles, many species of Xyleborus are also known from the earliest examinations of tropical scolytine fauna on different continents [40–42]. Their presence over several centuries may well be due to anthropogenic factors, but could also indicate prehistorical dispersal. It is therefore uncertain if human transport alone can explain the pantropical distribution of these beetles.
In order to test the relative importance of ancient vs. modern dispersal for current geographical distributions, we reconstructed the biogeographical history of four Xyleborus species. The species were chosen because they are largely pantropical in distribution and are among the most numerous ambrosia beetles wherever they are found [30, 43–45]. Due to their nearly global distribution, it has not been possible to determine their geographical origin. We envision three alternative scenarios for the biogeographical history of these species: (i) If Xyleborus beetles became widely established before the dawn of human influence on species distributions, we expect evidence for ancient establishment in multiple areas. Ancient distributions would be reflected in multiple divergent clades that are largely restricted to specific regions. High haplotypic diversity within several geographical regions would furthermore imply stable populations since prehistorical time. (ii) Contrarily, if a current pantropical distribution is due largely to multiple human-aided dispersal events, we expect the source area to be represented by a single genetically diverse clade. Very recent colonisations of new regions would then be evident from young derived clades nested within the aforementioned clade. (iii) A third possible result is a combination of the features from (i) to (ii), which would result from ancient dispersal to multiple regions followed by recent human-aided introductions to the same and other areas.
Females of the four focal species were sampled from four continents to 20 countries (Additional file 1: Table S1). Beetles were collected by hand, or in ethanol-baited flight intercept traps, or light traps, and preserved in 95 % ethanol. All specimens were identified based on morphology (BHJ and LRK) before extraction. We extracted DNA from 53 Xyleborus ferrugineus, 62 X. affinis, 25 X. volvulus, and 30 X. perforans individuals. We included sequences from Coccotrypes cyperi, C. advena, Ozopemon brownei, and Xylosandrus morigerus as outgroups , and we included data from nine related Xyleborus species to test for monophyly in our focal species (Genbank acc. no and metadata in Additional file 1: Table S3).
We sequenced one mitochondrial gene, cytochrome oxidase I (COI), and one nuclear gene, Elongation Factor 1-alpha (EF1α). The rapidly evolving mitochondrial gene is suited for examining population structure, while the nuclear gene may be more informative when attempting to resolve older relationships among species. A subset of individuals representing all identified COI haplo-groups was selected for EF1α sequencing (Additional file 1: Table S1). Before combining the two loci in an analytical framework, we assessed their phylogenetic concordance with a Congruence Among Distance Matrices (CADM) test . Using the function tanglegram from the R package dendextend, we visualized separate phylogenetic reconstructions for each locus—and their concordance—as a tanglegram.
We extracted genomic DNA using E.Z.N.A® Tissue DNA Kit (Omega Bio-Tek, Atlanta, GA, USA) following the manufacturer’s protocol. 610 base pairs (bp) of COI gene and 1043 bp of the EF1α nuclear gene were amplified and sequenced using primers from Normark et al. [48; Additional file 1: Table S2]. PCR reactions were performed in 25 μL volume with 10 μM of each forward and reverse primer, 0.25 mM dNTPS, 0.625 U HotStar Taq DNA Polymerase, 17.4 μL ddH2,0, 2.5 μL of 10x PCR buffer, 25 mM MgCl2 and 1.0 μL of DNA template using the following PCR cycle: 95°, 15 min; 35 cycles [94°, 30s; 48°, 45 s; 72°, 60s]; 72°, 7 min. PCR reactions for EF1α were performed in 26 μL volume containing similar ingredients and concentrations as for CO1 gene amplifications, but with 2.0 μL of DNA. A touchdown PCR profile was employed for this nuclear gene: 95°, 15 min; 13 cycles TD[94°, 30s; 58–44°, 45 s; 72°, 30s]; 26 cycles [94°, 30s; 44°, 45 s; 72°, 30s]; 72°, 6 min. We purified amplified DNA using ExoSAP (Exonuclease I—Shrimp Alkaline Phosphatase) and sequenced in both directions using standard protocol for ABI BigDye® Terminator v3.1 Cycle sequencing kit (Applied Biosystems).
Sequence contigs were assembled and edited using SEQMAN II, a contig assembly module from DNASTAR Lasergene6®, and aligned using ClustalW  in BioEdit 7.0.1 . The alignments were further refined using amino acid translation. We found no indications of nuclear insertions (NUMTs) of the COI gene in our sequenced data (all single-peak chromatograms, no stop codons or indels). COI (acc.no: KP941137-KP941327) and EF1α (acc.no: KP941328-KP941418) sequences are deposited in Genbank.
The phylogenetic and biogeographical reconstruction included both genes and all sequence data (194 individuals: the four focal species, outgroup, and ingroup), and seven geographical regions (Afrotropical mainland, Afrotropical islands (Indian Ocean), Neotropical mainland, Neotropical islands (Pacific Ocean), Indo-Malaysia, Australasia, Australasian small islands (Pacific Ocean). Haplotype networks were constructed independently for both COI and EF1α for this dataset using haploNet from the pegas R package . 2) Ancestral reconstruction of geographical distribution was done by defining geographical location as discrete traits in BEAST, using a symmetric trait substitution model. Social network was inferred with BSSVS  in BEAST. The biogeographical reconstruction was performed using a conditional reference prior . Since we were interested in identifying supported diffusions (dispersal/migration routes) for each species, we performed analyses on each individual species set using the same approach as outlined above. The results from these latter analyses were further analysed in SPREAD  (a tool developed for epidemiology studies, which can also be used to visualize biogeographical histories in general, e.g. [55, 56]). A Bayes factor cutoff of 10 was implemented so that only highly supported diffusions were returned by the software.
The best-fit evolutionary models of sequence evolution were determined with jModeltest 2.1.3 [57, 58] (Additional file 1: Table S4 and Additional file 1: Table S5), and phylogenetic and biogeographical reconstructions were performed with BEAST v1.8.2 . The Yule speciation process tree prior was used in the reconstruction which included all species and the constant size coalescence tree prior was used for reconstructions of individual species datasets. The larger dataset was run for 70 M and the smaller datasets for 20 M generations. Mean values and effective sample size (ESS values) for all parameters were obtained using Tracer  with 10 % burn-in (Additional file 1: Table S4 and Additional file 1: Table S5). Maximum clade credibility trees were obtained using TreeAnnotator from the BEAST package.
In addition to phylogenetic and biogeographic analyses, we performed a molecular dating analysis using COI and EF1α data. For time calibration, we referred to an unpublished phylogenetic analysis of Scolytinae (317 spp.) which include four fossil calibration points. The sister group of Xyleborus in this large analysis was used as the outgroup in the calibration analysis (Additional file 1: Table S3). In the original higher level analysis, Xyleborus was represented by X. alluaudi and X. affinis. We thus included X. alluadi in order to evaluate the age estimates for the MRCA for these two species in both analyses (as a qualitative test). The stem age of Xyleborus was set to 17.3 Ma based on our global Scolytinae dating, with a normal distribution prior (stdev=0.173 myr). Since our system is relatively young and contain similar species, we opted for a strict clock rate prior with rate estimation for both markers and a Yule speciation prior in BEAST.
Tests of molecular variance (AMOVA) were performed in Arlequin v3.5 . We performed two AMOVA tests based on different sets of geographical regions. In the first test, we clustered populations into the same seven regions as used in the large phylogenetic reconstruction (Afrotropical mainland, Afrotropical islands (Indian Ocean), Neotropical mainland, Neotropical islands (Pacific Ocean), Indo-Malaysia, Australasia, Australasian small islands (Pacific Ocean)); in the second test, we merged the oceanic regions with their continental counterparts giving a total of four regions. COI and EF1α were analysed separately. We also estimated genetic diversity in terms of segregating sites, nucleotide and haplotype diversity (S, π and H) along with Tajima's D, which tests for neutrality and recent population expansion or contraction, using DNAsp . The association between genetic and geographical distance was evaluated with mantel tests performed in R using mantel.rtest . To evaluate the distribution of genetic distances among haplotypes within species, mismatch distribution plots for COI were created in R  using the function MMD from the pegas package .
Each species contained distinct intraspecific clades separated by a large number of substitutions, especially in the COI data (Fig. 1). Considerable geographical mixing was observed within subclades in three of the four species (X. perforans, X affinis and X. volvulus) except for one X. affinis clade that was purely Afrotropical. Xyleborus ferrugineus differed from the other species, containing numerous old clades that were generally restricted to distinct geographical regions. The differences noted between X. ferrugineus and the other three species were also apparent in the COI mismatch distribution plots (Additional file 1: Figure S3) showing large genetic distances between haplotypes in X. ferrugineus, while shorter genetic distances predominate in the mismatch distributions for the other three species.
The biogeographical reconstruction revealed that X. ferrugineus and X. volvulus most likely originated in the Neotropics (pie charts on phylogeny in Fig. 1). Neotropical and Afrotropical origin were both highly likely for X. affinis. For X. perforans, an Australasian origin had the highest probability, but Neotropical and Indo-Malaysian origins also had substantial probability scores.
COI haplotype networks were largely congruent with the phylogenies (Fig. 1). We found large genetic distances among clusters in all four species, indicating ancient splits between intra-specific lineages. Haplotype clusters were not geographically homogeneous for X. affinis, X. volvulus and X. perforans; clusters contained individuals predominantly from one geographical region together with one or a few individuals from a separate region. A few single haplotypes were shared by individuals from different geographical regions. Sequences of X. ferrugineus showed stronger geographical structuring, with geographically homogeneous haplotype clusters. We observed only one example of a very recent dispersal event (between the Neotropics and Afrotropics) in this species. There was also a tendency towards higher haplotypic diversity within geographical regions in X. ferrugineus, in particular the Afrotropical and Neotropical regions. The EF1α haplotype networks revealed little additional information to those of COI, but rather supported the general patterns seen in Fig. 1 (Additional file 1: Figure S4). As expected, there was significantly less fine scale resolution in the EF1α haplotype network compared to the COI network, which is due to the lower level of genetic variation in EF1α (Table 2).
Tests of molecular variance (AMOVA) for the four species and both genetic markers. The seven regions are Afrotropics, Afrotropical islands (Indian Ocean), Neotropics, Neotropical islands (Pacific Ocean), Indo-Malaysia, Australasia, Australasia small islands (Pacific Ocean). In the four regions the oceanic regions were merged with their continental counterparts
Source of variation
Sum of squares
Sum of squares
ФCT : 0.59
ФCT : 0.34
Among populations within regions
ФSC : 0.28
ФSC : 0.55
ФST : 0.71
ФST : 0.71
ФCT : 0.59
ФCT : 0.18
Among populations within regions
ФSC : 0.34
ФSC : 0.66
ФST : 0.73
ФST : 0.71
ФCT : 0.45
ФCT : 0.10
Among populations within regions
ФSC : 0.17
ФSC : 0.48
ФST : 0.55
ФST : 0.53
Among populations within regions
ФSC : 0.58
ФSC : 0.45
ФST : 0.44
ФST : 0.44
ФCT : 0.09
Among populations within regions
ФSC : 0.49
ФSC : 0.54
ФST : 0.54
ФST : 0.53
Among populations within regions
ФSC : 0.50
ФSC : 0.13
ФST : 0.01
ФST : 0.02
ФCT : 0.21
ФCT : 0.01
Among populations within regions
ФSC : 0.24
ФSC : 0.36
ФST : 0.40
ФST : 0.36
ФCT : 0.24
Among populations within regions
ФSC : 0.79
ФST : 0.84
Three measures of genetic diversity and Tajima’s D for the four focal species
Number of segregating sites (S)
Total number of mutations
Nucleotide diversity (π)
Standard deviation of π
Total number of haplotypes
Standard deviation of Hd
For all four Xyleborus species, we found conclusive evidence for colonization of separate continents millions of years ago (see Fig. 2). With the exception of X. ferrugineus, these species also show evidence of accelerated rates of colonization in recent times (Figs. 1 and 2), which likely stems from increased introductions through anthropogenic effects. However, the geographically broad distribution of these species, which was discovered during the earliest insect inventories in tropical Africa, South America and SE-Asia [40–42], is not primarily a result of introductions just prior to these first inventories. The results from our dating analysis and biogeographically structured COI haplotype clusters clearly indicate ancient dispersal and demonstrates that anthropogenic effects are not a necessity for the global spread of these or similar species.
The emerging picture of a clade with many widespread species dispersing since the time of their origin is a relatively novel finding and distinguishes this group from many other invasive organisms such as certain termites  and ants [66, 67] and several other bark and ambrosia beetles. The coffee berry borer (Hypothenemus hampei), for instance, has considerable genetic variation only in its native African region, with recent dispersal across continents following the cultivation of coffee . Termites are particularly comparable in terms of ecology and includes the dry-wood inhabiting termites; despite dozens of widespread tramp species, none of these are ancient globetrotters, which is surprising given a flexible reproductive system and capacity for rafting in draft wood . In this regard, termites show many similarities with a group of ambrosia beetles in the weevil subfamily Platypodinae, where a handful of recent world expansions is the exception among otherwise highly endemic beetles . In this perspective it is truly remarkable that several species of Xyleborus have such a genuinely wide distribution pattern. However, they do possess a truly remarkable mode of reproduction that possibly holds the key for understanding their wide distribution. They are permanent inbreeders, hence inseminated before leaving their nest, and can even mate with sons produced from unfertilized haploid eggs . Although there are many other species of Xyleborini that have very narrow distributions, species from this tribe are, together with another group of permanent inbreeders in the genus Hypothenemus, clearly overrepresented among the most widespread forest insects, indicative of their colonisation ability [70, 71].
Inferring the geographical origin of widespread Xyleborus is not straightforward given their complex patterns of repeated dispersal across and between continents. A close evolutionary relationship between a large clade of endemic Neotropical genera (Sampsonius, Dryocoetoides, Coptoborus and Theoborus) and the clade of widespread Xyleborus species may indicate that the Neotropical region is important in their evolutionary history . However, some of the closest relatives of our focal species are endemic to the Afrotropical region (e.g. X. alluaudi and related species) , which complicates the inference of biogeographical history from such higher level taxonomy, and highlights that dispersal, not vicariance, is likely the primary driver of biogeographical structure in this group. Our biogeographical analyses support the Neotropical association in three of the focal species (ambiguous in X. perforans). It should be pointed out, however, that the reconstructed region for X. affinis seems to be strongly affected by the underlying node. Since the topology at this node has poor support (posterior probability=0.87), it may be unwise to emphasize this state when determining the ancestral region. This, together with the high haplotypic diversity in the Afrotropics, suggests an Afrotropical origin isfor X. affinis.
Given that ancient long-range dispersal has occurred repeatedly in this clade of Xyleborus, it may not be unrealistic to expect cryptic speciation in the absence of sexual selection [70, 73]. Several previous synonyms of X. ferrugineus has recently been resurrected to full species status [71, 74], which makes us alert to this possibility. We note that the New Caledonia population, which likely constitutes a very old dispersal event from the Neotropics, was separate from the other populations in both data sets (COI and EF1α; Additional file 1: Figure S2), but without any morphological differentiation. Xyleborus bispinatus and X. impressus, which were previously considered as X. ferrugineus, can on the other hand be readily distinguished genetically and morphologically and hence they are more clearly separated than the New Caledonia population. It is furthermore possible that broader geographical sampling of X. ferrugineus also will reduce some of the genetic gaps observed between populations of this species. Questions have been raised about the validity of species status of the near identical X. volvulus and X. perforans [45, 75]. Our data nevertheless solved these issues confidently and we conclude that X. volvulus and X. perforans represent two individual species—both with regard to phylogenetic monophyly, and to their separate biogeographical origins. However, they are both globally widespread, and not restricted to the Neotropics and the Old World, respectively, as previously assumed .
We have tested the hypothesis that the global distributions in many Xyleborus species were fostered by human transportation. We found ample evidence for recent colonisations of new regions in three of the four examined Xyleborus species, which is in line with this hypothesis. However, the genetic data also strongly suggest that all four species were present across most of their current distributional range several million years ago, which was clearly before humans crossing the oceans. The combination of a pre-adapted biology for propagule establishment far from the source population, and the frequent transportation of wood material by humans in modern times, has radically increased migration rates in these beetles. It is therefore not surprising that many species of Xyleborus are among the most frequently collected species in any tropical or subtropical region.
Availability of supporting data
Supporting figures and tables are found in the additional file. All analysed DNA sequences are available in Genbank. COI acc.no: KP941137-KP941327), EF1α acc.no: KP941328-KP941418, accession numbers for ingroup and outgroup sequences are listed in (Additional file 1: Table S3).
We thank Jiri Hulcr, Carlos Ruiz and Anthony Cognato for providing samples and the Research Council of Norway for funding (grant no.: 170565/V40).
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Walker B, Steffen W. An overview of the implications of global change for natural and managed terrestrial ecosystems. Conserv Ecol. 1997;1(2):2.Google Scholar
- Gurevitch J, Padilla DK. Are invasive species a major cause of extinctions? TREE. 2004;19(9):470–4.PubMedGoogle Scholar
- McNeely JA. Human dimensions of invasive alien species. Washington: Island Press; 2005.Google Scholar
- Cronk QCB, Fuller JL. Plant invaders: the threat to natural ecosystems. London: Chapman & Hall; 1995.Google Scholar
- Pimentel D, Zuniga R, Morrison D. Update on the environmental and economic costs associated with alien-invasive species in the united states. Ecol Econ. 2005;52(3):273–88.View ArticleGoogle Scholar
- Williamson M, Fitter A. The varying success of invaders. Ecology. 1996;77(6):1661–6.View ArticleGoogle Scholar
- Kolar CS, Lodge DM. Progress in invasion biology: predicting invaders. TREE. 2001;16(4):199–204.PubMedGoogle Scholar
- Simberloff D. The role of propagule pressure in biological invasions. Annu Rev Ecol Evol Syst. 2009;40:81–102.View ArticleGoogle Scholar
- Simberloff D. Introduced insects: a biogeographic and systematic perspective. In: Mooney HA, Drake JA, editors. Ecology of ecological invasions of north america and hawaii. New York: Springer; 1986. p. 3–26.View ArticleGoogle Scholar
- Baker HG. Self-compatibility and establishment after ‘long-distance’ dispersal. Evolution. 1955;9:347–9.View ArticleGoogle Scholar
- Jordal BH, Beaver RA, Kirkendall LR. Breaking taboos in the tropics: inbreeding promotes colonization by wood-boring beetles. Global Ecol Biogeogr. 2001;10(4):345–58.View ArticleGoogle Scholar
- Allendorf FW, Lundquist LL. Introduction: population biology, evolution, and control of invasive species. Conserv Biol. 2003;17(1):24–30.View ArticleGoogle Scholar
- Price SC, Jain SK. Are inbreeders better colonizers? Oecologia. 1981;49(2):283–6.View ArticleGoogle Scholar
- Terhivuo J, Saura A. Dispersal and clonal diversity of north-European parthenogenetic earthworms. Biol Invasions. 2006;8(6):1205–18.View ArticleGoogle Scholar
- Haack RA, Rabaglia RJ. Exotic bark and ambrosia beetles in the USA: potential and current invaders. In: Peña JE, editor. Potential invasive pests of agricultural crops. Wallingford: CAB International; 2013. p. 48–74.View ArticleGoogle Scholar
- Andersen HF, Jordal BH, Kambestad M, Kirkendall LR. Improbable but true: the invasive inbreeding ambrosia beetle xylosandrus morigerus has generalist genotypes. Ecol Evol. 2012;2(1):247–57.PubMed CentralView ArticlePubMedGoogle Scholar
- Darwin C. The effects of cross and self fertilisation in the vegetable kingdom. London: John Murray; 1876.View ArticleGoogle Scholar
- Blackburn TM, Lockwood JL, Cassey P. The influence of numbers on invasion success. Mol Ecol. 2015;24(9):1942–53.View ArticlePubMedGoogle Scholar
- Bock DG, Caseys C, Cousens RD, Hahn MA, Heredia SM, Hübner S, Turner KG, Whitney KD, Rieseberg LH. What we still don't know about invasion genetics. Mol Ecol. 2015;24(9):2277–97.Google Scholar
- Avilés L, Purcell J. The evolution of inbred social systems in spiders and other organisms: from short-term gains to long-term evolutionary dead ends. Adv Study Behav. 2012;44:99–133.View ArticleGoogle Scholar
- Crnokrak P, Barrett SCH. Perspective: purging the genetic load: a review of the experimental evidence. Evolution. 2002;56(12):2347–58.View ArticlePubMedGoogle Scholar
- Courchamp F, Clutton-Brock T, Grenfell B. Inverse density dependence and the allee effect. TREE. 1999;14:405–10.PubMedGoogle Scholar
- Marini L, Haack RA, Rabaglia RJ, Toffolo EP, Battisti A, Faccoli M. Exploring associations between international trade and environmental factors with establishment patterns of exotic scolytinae. Biol Invasions. 2011;13(10):2275–88.View ArticleGoogle Scholar
- Beaver RA. Insect—fungus relationships in the bark and ambrosia beetles. In: Wilding N, Collins NM, Hammond PM, Webber JF, editors. Insect—fungus interactions. London: Academic; 1989. p. 121–43.Google Scholar
- Hulcr J, Mogia M, Isua B, Novotny V. Host specificity of ambrosia and bark beetles (Col., curculionidae: scolytinae and platypodinae) in a New guinea rainforest. Ecol Entomol. 2007;32(6):762–72.View ArticleGoogle Scholar
- Haack RA. Exotic bark—and wood-boring coleoptera in the United States: recent establishments and interceptions. Can J For Res-Rev Can Rech For. 2006;36(2):269–88.View ArticleGoogle Scholar
- Sauvard D, Branco M, Lakatos F, Faccoli M, Kirkendall L. Weevils and bark beetles (coleoptera, curculionoidea). chapter 8.2. BioRisk. 2010;4:219–66.View ArticleGoogle Scholar
- Beeson CFC. The ecology and control of the forest insects of India and the neighbouring countries. 2nd ed. Government of India: New Delhi; 1961.Google Scholar
- Browne FG. Some aspects of host selection among ambrosia beetles in the humid tropics of south-east Asia. Malayan Forestry. 1958;21:164–82.Google Scholar
- Browne FG. The biology of malayan scolytidae and platypodidae. Malayan Forest Records. 1961;22:1–255.Google Scholar
- Kalshoven LGE. Studies on the biology of indonesian scolytoidea. 4. Data on the habits of scolytidae. Tijdschr Entomol. 1958;101:157–80.Google Scholar
- Schedl KE. Breeding habits of arboricole insects in central Africa. In: Xth international congress of entomology: august 17–25 1956. Montreal. 1956;1:183–97.Google Scholar
- Bain J. Overseas wood-and bark-boring insects intercepted at New Zealand ports. Forest Research Institute, New Zealand Forest Service. 1977;63:28.Google Scholar
- Haack RA. Intercepted scolytidae (coleoptera) at U.S. Ports of entry: 1985–2000. Integr Pest Manage Rev. 2001;6:253–82.View ArticleGoogle Scholar
- Brockerhoff EG, Bain J, Kimberly M, Knizek M. Interception frequency of exotic bark and ambrosia beetles (coleoptera: scolytinae) and relationship with establishment in New Zealand and worldwide. Can J For Res-Rev Can Rech For. 2006;36:289–98.View ArticleGoogle Scholar
- Brockerhoff EG, Jones DC, Kimberley MO, Suckling DM, Donaldson T. Nationwide survey for invasive wood-boring and bark beetles (coleoptera) using traps baited with pheromones and kairomones. For Ecol Manage. 2006;228(1–3):234–40.View ArticleGoogle Scholar
- Kirkendall L, Faccoli M. Bark beetles and pinhole borers (curculionidae, scolytinae, platypodinae) alien to Europe. ZooKeys. 2010;56:227–51.View ArticleGoogle Scholar
- Haack RA, Petrice TR. Bark—and wood-borer colonization of logs and lumber after heat treatment to ISPM 15 specifications: the role of residual bark. J Econ Entomol. 2009;102(3):1075–84.View ArticlePubMedGoogle Scholar
- Brockerhoff EG, Kimberley M, Liebhold AM, Haack RA, Cavey JF. Predicting how altering propagule pressure changes establishment rates of biological invaders across species pools. Ecology. 2014;95(3):594–601.View ArticlePubMedGoogle Scholar
- Schedl KE. New bark beetles and ambrosia beetles. Stylops. 1935;4(12):273.Google Scholar
- Schreiner J. Neue tomicinen von der goldkuste Afrikas. Deutsche Entomologische Zeitschrift. 1882;26:246–8.Google Scholar
- Blandford WFH. Family scolytidae. Biologia Centrali-Americana, Coleoptera. 1897;4(6):145–84.Google Scholar
- Wood SL. Bark and ambrosia beetles of south America (coleoptera, scolytidae). Provo, Utah: Brigham Young University; 2007.Google Scholar
- Wood SL. The bark and ambrosia beetles of north and central America (coleoptera: scolytidae), a taxonomic monograph. Gt Basin Nat Mem. 1982;6:1–1359.Google Scholar
- Hulcr J, Cognato AI. Xyleborini of New Guinea, a taxonomic monograph (coleoptera: curculionidae: scolytinae). Lanham, Maryland: Thomas Say Publications in Entomology, Entomological Society of America; 2013.Google Scholar
- Jordal BH. Elongation factor 1α resolves the monophyly of the haplodiploid ambrosia beetles xyleborini (coleoptera: curculionidae). Insect Mol Biol. 2002;11:453–65.View ArticlePubMedGoogle Scholar
- Campbell V, Legendre P, Lapointe F-J. The performance of the congruence among distance matrices (CADM) test in phylogenetic analysis. BMC Evol Biol. 2011;11(1):64.PubMed CentralView ArticlePubMedGoogle Scholar
- Normark BB, Jordal BH, Farrell BD. Origin of a haplodiploid beetle lineage. Proc R Soc B. 1999;266:2253–9.PubMed CentralView ArticleGoogle Scholar
- Thompson JD, Higgins DG, Gibson TJ. Clustal-W—improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific Gap penalties and weight matrix choice. Nucleic Acids Res. 1994;22(22):4673–80.PubMed CentralView ArticlePubMedGoogle Scholar
- Hall TA. BioEdit: a user-friendly biological sequence alignment editor and analysis program for windows 95/98/NT. Nucleic Acids Symp Ser. 1999;41:95–8.Google Scholar
- Paradis E. Pegas: an R package for population genetics with an integrated-modular approach. Bioinformatics. 2010;26(3):419–20.View ArticlePubMedGoogle Scholar
- Lemey P, Rambaut A, Drummond AJ, Suchard MA. Bayesian phylogeography finds its roots. PLoS Comput Biol. 2009;5(9), e1000520.PubMed CentralView ArticlePubMedGoogle Scholar
- Ferreira MAR, Suchard MA. Bayesian analysis of elapsed times in continuous-time markov chains. Canadian Journal of Statistics-Revue Canadienne De Statistique. 2008;36(3):355–68.View ArticleGoogle Scholar
- Bielejec F, Rambaut A, Suchard MA, Lemey P. SPREAD: spatial phylogenetic reconstruction of evolutionary dynamics. Bioinformatics. 2011;27(20):2910–2.PubMed CentralView ArticlePubMedGoogle Scholar
- Wang Q, Liu J, Allen GA, Ma Y, Yue W, Marr KL, Abbott RJ. Arctic plant origins and early formation of circumarctic distributions: a case study of the mountain sorrel. Oxyria digyna New Phytol. 2016;209(1):343–53.Google Scholar
- Guarnizo CE, Werneck FP, Giugliano LG, Santos MG, Fenker J, Sousa L, D’Angiolella AB, dos Santos AR, Strüssmann C, Rodrigues MT et al. Cryptic lineages and diversification of an endemic anole lizard (squamata, dactyloidae) of the cerrado hotspot. Mol Phylogenet Evol. 2016;94(Part A):279–89.Google Scholar
- Darriba D, Taboada GL, Doallo R, Posada D. JModelTest 2: more models, new heuristics and parallel computing. Nat Methods. 2012;9(8):772–2.Google Scholar
- Guindon S, Gascuel O. A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003;52(5):696–704.View ArticlePubMedGoogle Scholar
- Drummond AJ, Suchard MA, Xie D, Rambaut A. Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol. 2012;29(8):1969–73.PubMed CentralView ArticlePubMedGoogle Scholar
- Rambaut A, Suchard MA, Xie D, Drummond AJ. Tracer v1.6, Available from http://beast.bio.ed.ac.uk/Tracer. 2014.
- Excoffier L, Laval G, Schneider S. Arlequin (version 3.0): an integrated software package for population genetics data analysis. Evol Bioinforma. 2005;1:47–50.Google Scholar
- Librado P, Rozas J. DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics. 2009;25(11):1451–2.View ArticlePubMedGoogle Scholar
- Dray S, Dufour A-B. The ade4 package: implementing the duality diagram for ecologists. J Stat Softw. 2007;22(4):1–20.View ArticleGoogle Scholar
- R Development Core Team. R: a language and environment for statistical computing. In. Vienna, Austria, http://www.R-project.org; 2014.
- Evans TA, Forschler BT, Grace JK. Biology of invasive termites: a worldwide review. Annu Rev Entomol. 2013;58(1):455–74.View ArticlePubMedGoogle Scholar
- Suarez AV, Holway DA, Case TJ. Patterns of spread in biological invasions dominated by long-distance jump dispersal: insights from argentine ants. Proc Natl Acad Sci. 2001;98(3):1095–100.PubMed CentralView ArticlePubMedGoogle Scholar
- Foucaud J, Orivel J, Loiseau A, Delabie JH, Jourdan H, Konghouleux D, Vonshak M, Tindo M, Mercier JL, Fresneau D. Worldwide invasion by the little fire ant: routes of introduction and eco‐evolutionary pathways. Evol Appl. 2010;3(4):363–74.Google Scholar
- Gauthier N. Multiple cryptic genetic units in hypothenemus hampei (coleoptera: scolytinae): evidence from microsatellite and mitochondrial DNA sequence data. Biol J Linn Soc. 2010;101(1):113–29.View ArticleGoogle Scholar
- Jordal BH. Molecular phylogeny and biogeography of the weevil subfamily platypodinae reveals evolutionarily conserved range patterns. Mol Phylogenet Evol. 2015;92:294–307.View ArticlePubMedGoogle Scholar
- Kirkendall LR, Biedermann PHW, Jordal BH. Evolution and diversity of bark and ambrosia beetles. In: Vega FE, Hofstetter RW, editors. Bark beetles: biology and ecology of native and invasive species. Waltham, Massachusetts: Academic; 2014.Google Scholar
- Kirkendall LR, Jordal BH. The bark and ambrosia beetles (curculionidae, scolytinae) of cocos island, costa rica and the role of mating systems in island zoogeography. Biol J Linn Soc. 2006;89(4):729–43.View ArticleGoogle Scholar
- Cognato AI, Hulcr J, Dole SA, Jordal BH. Phylogeny of haplo-diploid, fungus-growing ambrosia beetles (coleoptera: curculionidae: scolytinae: xyleborini) inferred from molecular and morphological data. Zool Scr. 2011;40(2):174–86.Google Scholar
- Kirkendall LR. Ecology and evolution of biased sex ratios in bark and ambrosia beetles (scolytidae). In: Wrensch DL, Ebbert MA, editors. Evolution and diversity of sex ratio: Insects and mites. New York: Chapman and Hall; 1993. p. 235–345.View ArticleGoogle Scholar
- Rabaglia RJ. The validity of xyleborus impressus eichhoff (coleoptera : curculionidae : scolytinae) as distinct from xyleborus ferrugineus (fabricius). Coleopt Bull. 2005;59(2):261–6.View ArticleGoogle Scholar
- Wood SL, Bright D. A catalog of scolytidae and platypodidae (coleoptera). part 2: taxonomic index. Gt Basin Nat Mem. 1992;13:1–1553.Google Scholar