Skip to main content

Molecular phylogeny of bark and ambrosia beetles reveals multiple origins of fungus farming during periods of global warming



Fungus farming is an unusual life style in insects that has evolved many times in the wood boring weevils named ‘ambrosia beetles’. Multiple occurrences of this behaviour allow for a detailed comparison of the different origins of fungus farming through time, its directionality, and possible ancestral states. We tested these hypotheses with a phylogeny representing the largest data set to date, nearly 4 kb of nucleotides from COI, EF-1α, CAD, ArgK, 28S, and 200 scolytine taxa.


Phylogenetic analyses using Bayesian or parsimony approaches placed the root of Scolytinae close to the tribe Scolytini and Microborus, but otherwise indicated low resolution at older nodes. More recent clades were well resolved, including ten origins of fungus farming. There were no subsequent reversals to bark or phloem feeding in the fungus farming clades. The oldest origin of fungus farming was estimated near 50 Ma, long after the origin of Scolytinae (100-120 Ma). Younger origins included the species rich Xyleborini, dated to 21 Ma. Sister group comparisons and test of independence between traits indicated that neither gregarious larval feeding nor regular inbreeding by sibling mating was strongly correlated with the origin of fungus farming.


Origins of fungus farming corresponded mainly with two periods of global warming in the Cenozoic era, which were characterised by broadly distributed tropical forests. Hence, it seems likely that warm climates and expanding tropical angiosperm forests played critical roles in the successful radiation of diverse fungus farming groups. However, further investigation will likely reveal additional biological factors that promote fungus farming.


Bark and ambrosia beetles in the weevil subfamily Scolytinae are remarkably diverse in their ecological adaptations to a variety of habitats [1]. With more than 6,000 species currently recognised, they comprise about 10 percent of the total weevil diversity, and thus constitute one of the greatest species radiations on earth [2]. Scolytine beetles are generally enormously abundant and dominate forest insect communities associated with recently dead wood and other lignified plant material [35]. Their ubiquitous presence in nearly all woody material indicates that these beetles are the most important organisms in the early stages of forest decomposition. In these habitats they excavate diverse tunnel systems reflecting different reproductive modes and variable diets in different plant tissues. A large scale phylogeny of Scolytinae will therefore illuminate many aspects of key evolutionary traits in this group of beetles.

One of the most successful ecological adaptations in bark and ambrosia beetles is their symbiotic relationship with microbes. This relationship provides nutritional enrichment of dead plant material in exchange for transmission between plant resources. These microbes include bacteria, yeast and mycelial fungi. While the exact role of the first two is not known [6, 7], the multicellular fungi contribute an important but variable component of a wood boring beetle’s diet [8, 9]. In about 2,000 species of Scolytinae, Ophiostomales and Microascales fungi are the sole source of food for both larvae and adults and are actively cultivated by the beetle. This extraordinary agricultural system shares many similarities with fungus farming in ants and termites [1012]. These are the ambrosia beetles and they include common pest species such as the striped ambrosia beetle Trypodendron lineatum and the red bay ambrosia beetle Xyleborus glabratus. In addition to scolytine ambrosia beetles, fungus farming also occurs in another weevil subfamily, the Platypodinae [13, 14]. Fungi are transmitted as spores in specialised cuticular pockets in ambrosia beetles, and are inoculated in the wood during the excavation of a new tunnel system in their host plant. They grow fine mycelia in the wood during egg laying [8]. After hatching, the larvae graze on a dense carpet of conidia which covers the tunnel walls. Ambrosia beetles are both functionally (food) and physiologically (hormones) dependent on these fungi. It has been documented that moulting and metamorphosis does not occur in the absence of fungal steroids [15]. It is therefore expected that a reversal to a wood based diet is not very likely once fungus farming has evolved.

The remaining species of Scolytinae mainly feed on bark and phloem, occasionally on seeds, and they take most of their nutrients from dead plant tissue excavated during tunnel construction [1]. These beetles are also associated with Ophiostomales fungi, which may provide food enrichment. Although it is not unusual that bark beetles have mycangia to facilitate transportation of fungal spores [8], they cannot make fungal gardens and have only a facultative association with fungi; hence, they can complete their life cycle without the addition of fungus to their diets. The fitness of certain bark beetle species (e.g. Dendroctonus) is nevertheless increased by addition of such fungi in their phloem diet [16, 17]. This is a selective advantage for effective transmission. One can therefore readily imagine that a transition from a primarily bark- and phloem-based diet to a nutritional dependence on fungi is not a particularly difficult evolutionary event. Fungus farming has evolved multiple times in Scolytinae because of this nutritional advantage [13, 18].

Multiple origins of fungus farming enable a meaningful comparison of the underlying ecological circumstances that may have spurred the transition from phloem to fungal feeding. Fungus farming beetles are found in wet tropical forests (with a few boreal exceptions), which is perhaps indicative of the ideal climate for ambrosia beetle species radiations. However, the ecological or climatic conditions under which fungus farming may have evolved and the timing of these feeding transitions have not been explored in detail.

Ambrosia beetles are also characterised by a generally gregarious feeding behaviour in the larvae, a behaviour that seems pre-adaptive for close inbreeding by sibling mating [19, 20]. Regular inbreeding is the rule in several unrelated scolytine lineages and is characterised by female-biased offspring sex ratios. Fungus farming is potentially advantageous to developing sibling mating by congregating brood members in areas of fungal growth, which facilitates efficient mating between siblings. The direction of evolution from gregarious inbreeding to fungus farming or vice versa, or the assessment of whether or not these are completely unlinked phenomena, has not yet been examined in a phylogenetic context but see [1, 8, 19, 21].

A detailed phylogenetic analysis of Scolytinae will enable us to test the sequence of the origin of fungus farming, gregarious feeding modes, and sibling mating. Hence, we have reconstructed the most comprehensively sampled phylogeny of Scolytinae to date, based on five genes for 200 Scolytinae species. This represents 123 of 250 genera in 24 of the currently 27 recognized tribes, including all known fungus farming lineages. Dates of evolutionary origin were estimated by calibrating molecular divergence rates using a significantly updated fossil record that set the minimum age of both Scolytinae and Curculioninae to more than 100 Ma [2224].


All phylogenetic analyses revealed limited phylogenetic structure at deeper nodes, with most of the strongly supported nodes occurring at tribal and generic levels. Differences between the various Bayesian and parsimony analyses were not strongly supported and we have therefore used the Bayesian analysis based on seven partitions to illustrate the major topological findings in this report (Figure 1). Incongruent signal in the five independent markers was revealed by an overall negative support from EF-1α or CAD in the partitioned Bremer analyses (Additional file 1: Figure S1). However, re-analyses with exclusion of one or more partitions did not improve overall node support. Analyses of amino acid translated data from the four protein encoding genes combined resulted in minimal resolution, mainly with the ingroup separate from the outgroups, and with Scolytini subtending the most basal node in Scolytinae.

Figure 1
figure 1

Molecular phylogeny of Scolytinae. Topology based on Bayesian analysis of 3,694 nucleotides from five unlinked gene fragments, using seven unlinked data partitions (mtDNA positions 1, 2, and 3, nucDNA protein encoding gene positions 1, 2, and 3, 28S rRNA). Posterior probability marked on nodes by ** = 100, * = > 95. Obligate fungus farming is marked in blue, with blue hatch mark signifying ambrosia symbiosis in a single species. Red dots indicate the inferred origin of regular sibling mating.

Phylogenetic patterns

Eleven of the 24 included tribes (Additional file 2: Table S1) were monophyletic in all or some analyses, including Bothrosternini, Corthylini, Crypturgini, Phloeotribini, Phrixosomatini, Scolytoplatypodini, Scolytini, Xyleborini and Xyloterini (Hyorrhynchini and Cactopinini were represented by single species). Micracidini included Cactopinini and these taxa were monophyletic in some analyses (Figure 1). Dryocoetini was paraphyletic with respect to Xyleborini and Ipini was paraphyletic with respect to Premnobiini; these two clades formed well supported sister groups. Scolytoplatypodini was nested within parts of Hexacolini (genera Scolytodes and Gymnochilus). Several well-sampled tribes such as Hypoborini and Polygraphini were monophyletic with the exclusion of one genus each. The majority of conifer associated Hylurgini (previously Tomicini) formed two main clades – one southern hemisphere clade consisting of Araucaria associated species, and one northern clade associated with boreal Pinaceae. Although only weakly supported, Scolytini was monophyletic and subtended the Hexacolini genus Microborus and all other Scolytinae species.

Timing and directionality of ambrosia feeding

Ten independent origins of fungus cultivation could be traced on the various topologies under accelerated parsimony optimisation (Figure 1), with eleven origins under delayed transformation. Two of these origins were the recent origins of single species in the genera Hypothenemus and Scolytodes. Hyorrhynchini was represented only by a single species (Sueus niisimai). The remaining seven groups of fungus farming beetles were monophyletic in all analyses, with the exception of Camptocerus that was sometimes paraphyletic with respect to Cnemonyx.

Minimum (crown) and maximum (stem) age estimates for ambrosia beetle clades are listed in Table 1. For the seven clades where crown age could be assessed, dates ranged from 8-48 Ma in analysis A (Figure 2; Additional file 3: Figure S2), with Corthylina as the likely earliest clade of fungus farming species. Xyloterini, Scolytoplatypodini and Camptocerus followed with ages of 40, 35 and 33 Ma. The origins of the Bothrosternus + Eupagiocerus and Xyleborini clades were younger, estimated to 20 and 21 Ma. The two different analyses provided very similar time estimates, with analysis B providing time estimates roughly 1-3 myr older (Table 1).

Table 1 Estimates of clade ages
Figure 2
figure 2

Dates of origins of key evolutionary traits in Scolytinae. Age of gregarious larval feeding (G), regular inbreeding by sibling mating (thin line), and ambrosia fungus feeding (grey box) in Scolytinae, based on crown ages as estimated in Beast (Analysis A). Stippled line and light grey indicate uncertainties associated with the lack of observations (Bothrosternini) or lack of phylogenetic sampling (Hyorrhynchini). Only those larvae that move freely and feed gregariously are considered truly ‘gregarious.’ Above, the Zachos curve showing variation in global temperature, and the posterior 95% distribution of crown age for Corthylina, Xyloterini, Scolytoplatypus, Camptocerus, Xyleborini, and the Bothrosternus-Eupagiocerus clade.

Each fungus farming lineage was sister to a bark feeding beetle lineage, with the exception of Bothrosternus-Eupagiocerus where the sister group Cnesinus consisted of pith feeders, and, in part, Xyleborini where the sister group Coccotrypes contained some seed feeding species (Table 2).

Table 2 Sister group contrasts of the fungus farming clade and its inferred closest relative

The relative timing of gregarious feeding and inbreeding

The relative timing of ambrosia fungus cultivation, regular inbreeding and gregarious feeding modes are depicted in Figure 2. The oldest origins of regular inbreeding were estimated to be 50.1 Ma (analysis A) for Hypothenemus + Ptilopodius, and 29.3 Ma for the inbreeding Dryocoetini (Ozopemon, Coccotrypes, Dryocoetiops) plus Xyleborini. In clades where both inbreeding and ambrosia fungus farming occur (n = 6), inbreedingevolved first twice, fungus farming three times, and simultaneously once. Gregarious larval feeding occured in six clades and always occurred before or simultaneously with fungus farming (n = 4), and inbreeding (n = 6). Independence between fungus farming and gregarious feeding mode, and between fungus farming and regular inbreeding by sibling mating, was not rejected, while independence between gregarious feeding and sibling mating was rejected (Table 3). Transition rates indicated that fungus farming is more likely to evolve in groups with regular inbreeding or gregarious feeding mode, with the highest transition rates to fungus farming found in lineages with regular inbreeding.

Table 3 Correlation of fungus farming, gregariousness and regular inbreeding


Fungus farming has evolved at least ten times in Scolytinae, in contrast to the single origin of fungus farming in attine ants and macrotermitine termites [27, 28]. Although limited resolution in tree topology was found in all types of analyses, the fungus farming taxa and sister lineages were well resolved. We found that all origins of fungus farming in Scolytinae were derived and a reversal to a non-fungal diet could not be traced on any of the tree topologies examined. It is noteworthy that although these findings are concordant with previous studies [18], our data are more complete in terms of taxon sampling, inclusive of a higher number of fungus farming taxa per clade. The only possible reversal indicated by some of the Bayesian topologies related to Camptocerus with respect to Cnemonyx. However, a more complete taxon sampling for these two genera showed that Camptocerus is indeed monophyletic [29]; Smith and Cognato, unpublished molecular data]. Our data thus corroborate the hypothesis that fungus farming is indeed a non-reversible evolutionary transition.

The many origins of fungus farming did not correlate strongly with some of the biological factors that benefit from a symbiotic relationship between fungi and beetles. Although we did observe a trend in fungus farming evolving more often in lineages with close inbreeding, the reverse transition rate (from outbreeding to inbreeding in fungus farming lineages) were negligibly low and the association between specific reproductive modes and fungus farming was not significantly correlated (see Table 3). This is perhaps most clearly illustrated by the complete lack of regular inbreeding in four ambrosia beetle lineages (Scolytodes unipunctatus, the genera Camptocerus and Scolytoplatypus, and the entire subtribe Corthylina). Furthermore is regular inbreeding the norm in true bark beetles such as in some Dendroctonus and Araptus [see e.g. [19]. Fungus farming is therefore a trait that at least sometimes evolves relatively independent of reproductive biology. Repeated origins of fungus farming must therefore be explained by additional ecological factors such as the frequent facultative association between bark beetles and fungi that grow in the phloem and bark of the host trees. Based on this perspective, it is notable that termites have only evolved fungus farming on one occasion even though fungus is an important food component for many other termite groups [27].

Compared with the timing of the origin of Scolytinae, more than 100 Ma, the development of obligate symbiotic fungus farming occurred relatively late. In all groups where an estimate of crown age was reliable, they revealed origins younger than 50 Ma, with 95% confidence interval ± 12 myr (Table 1). Xyloterini and Scolytoplatypodini had stem ages older than the crown age for Corthylina, which could potentially indicate a slightly older origin of fungus farming in these groups. However, it is equally likely that close relatives of these taxa were not included in our study which would overestimate the age of these fungus farmers. Regardless of these uncertainties, the Ophiostomales fungi have certainly existed much longer than the ambrosia beetles as shown by the multiple independent origins of the symbiotic fungi [18, 30], and thus have likely been nutritionally advantageous to the early lineages of bark beetles that preceded the first ambrosia gardeners. In light of the ubiquitous presence of ambrosia beetles in pantropical forests, and the likely early availability of ambrosia fungi, one may wonder why such a successful adaptation should have taken so long to evolve. There are two particularly relevant factors that may not have been optimal at the earliest stage of bark beetle evolution – tropical forest diversity and climate.

About 98 percent of the known ambrosia beetle fauna is tropical or subtropical [25], which emphasizes that fungal symbiosis is largely dependent on moist conditions in warm climates [8, 31]. Thus, the timing of modern moist tropical forests expansion may be relevant to the origin of fungus farming beetles. Elements of angiosperm-dominated tropical forests developed during the mid-Cretaceous, but did not radiate extensively until the Palaeocene or early Eocene era [3235]. This time period experienced a thermal maximum (PETM) of some 5-8 degrees warmer climate from 58 to 45 Ma [36]. Several groups of animals and plants showed increased diversification associated with the increasing angiosperm dominance [37], in particular during or just after PETM [33, 3841] when tropical elements dominated floras and faunas from the equator to mid-latitudes e.g. [35, 42, 43]. Corthylina, Xyloterini, Scolytoplatypus and Camptocerus originated during or immediately after PETM and had likely taken advantage of the large tropical angiosperm forests emerging during this time period.

The only group of fungus cultivating insects that may have occurred in the Cretaceous period is a related group of weevils in the subfamily Platypodinae. Recent studies are inconsistent about the phylogenetic position of these beetles, but they are definitely part of the advanced weevil radiation [13, 44, 45]. Although the timing of this group seems problematic as a consequence of a generally higher substitution rate at independent genetic loci see [13], a late Cretaceous origin at 100-80 Ma seems realistic based on molecular data [13, 45] and a fossil from Burmese amber (Grimaldi, pers. comm). Climate during this time period is less well understood, but was probably quite warm, dominated by Magnoliales and the early expanding Malphigiales [35, 46]. However, the greatest part of the platypodine radiation took place much later, with more than 90 percent of the diversity originating in the Eocene and later time periods [13].

It is interesting that fungus farming in ants and termites have similarly late origins as in most Scolytinae beetles. Attine ants first originated around 50 Ma [12], similar to Corthylina beetles. However, the major radiation of these ants occurred later, around 20 Ma, which corresponds to our estimates for Bothrosternus-Eupagiocerus and the great Xyleborini radiation. During this intermediate ‘Antarctic thawing’ period [36], which lasted some 10 million years, tropical climates again dominated near mid-latitudes [35]. This is also the time period when the fungus gardening termites (Macrotermitinae) diversified [47], after their origin in tropical rainforests of Africa [48].

The late origin of the greatest ambrosia beetle radiation in Xyleborini is well supported by our data. Stem age was only 23 Ma (Additional file 3: Figure S2) for a clade that is closely related to bark beetles in the genera Coccotrypes, Dryocoetiops and Ozopemon[4952]. The species diversity in Xyleborini is therefore unparalleled by any other ambrosia beetle lineage or other scolytine lineage. A recent origin of Xyleborini fits well with their absence from Dominican amber, a fossil source otherwise rich on older ambrosia beetle groups such as Corthylina and Platypodinae [53]. The great diversity of Xyleborini stands in contrast to a relatively modest diversity in the other clades of scolytine ambrosia beetles, particularly so in perspective of time. The reason for their great diversity is unclear. There are at least nine other scolytine clades of ambrosia beetles and only three of these are marginally more diverse than their sister group (see Table 2). Xyleborini are also characterized by regular inbreeding by sibling mating which is generally a great success factor in scolytine evolution, including Hypothenemus and related genera in Cryphalini. However, among the seven origins of regular inbreeding, only Xyleborini and inbreeding Cryphalini are more diverse than their sister group, so evidently there is no direct connection between inbreeding and diversification. There is therefore nothing overtly unique with this group of beetles compared with other ambrosia beetle lineages. To conclude, Xyleborini is most likely diverse because of chance effects, evolving at the right time (global warming in the ‘Antarctic thawing’ period) in the most productive parts of the globe, in the tropical regions.


This study has provided the most accurate time estimate to date for the origin of fungus farming in bark and ambrosia beetles. Many origins occurred at different time periods, which resulted in a highly variable diversity for each lineage per time unit. However, they all shared a relatively late origin, which most likely occurred in periods of global warming and expanding wet tropical forests. Our taxonomic coverage was generally broad, which resulted in precise estimates for most origins of fungus farming. However, these estimates were less certain in Scolytoplatypodini, Xyloterini and Hyorrhynchini in which sister group relations were more distant. Additional research must aim at closing these taxonomic gaps.

The present study has furthermore demonstrated the challenge connected with resolving highly diverse taxon groups. Almost 4000 nucleotides were included from five independent molecular markers, but only the more derived clades were resolved, including all origins of fungus farming. Scolytinae is a very interesting group for many other kinds of evolutionary studies and detailed hypotheses may be tested if also deeper phylogenetic nodes obtain better resolution. As it currently stands, more protein encoding genes must be optimised to obtain such goals – a project of high priority in our laboratories.


Taxon sampling and data acquisition

Individuals used for DNA extraction and PCR are listed in Additional file 2: Table S1, inclusive of country of origin and genbank accession numbers. Three closely related outgroups were included. Platypodinae and Coptonotus were excluded because of their uncertain phylogenetic placement and the fact that Platypodinae has generally very long branch lengths which may influence the rooting at the base of Scolytinae [13].

Ambrosia beetles in the subfamily Scolytinae are defined as beetles actively cultivating fungi, and which differ from other bark and seed beetles by their exclusively fungal diet [8]. Inbreeding is defined as regular mating between siblings, with strongly female biased offspring sex ratio [19]. Gregarious feeding is defined by clustered siblings feeding and congregated development (without independent larval tunnels and pupal cradles), often forming variously rounded or cave-like structures in the wood. The character states for each included taxon are listed in Additional file 4: Table S2.

DNA extractions, PCR reactions, purifications and DNA sequencing reactions followed the protocols given in Jordal et al. [13].

Alignment and phylogenetic analyses

All protein encoded gene sequences were aligned by eye, with introns identified by insertions demarcated by the initiating GT and terminating AG motif. These were excised before further analyses. Alignments of rDNA sequences from the D2-D3 domains of the large ribosomal subunit 28S were initially aligned by the MUSCLE software [54] using default parameters, then re-aligned after pruning long expansion segments in Phloeoborus sp, Diamerus curvifer, Dolurgocleptes punctifer and Dactylipalpus grouvellei. The resulting Muscle alignment was adjusted slightly according to a secondary structure model for Scolytinae [55].

Phylogenetic trees were reconstructed in a Bayesian framework using Mr Bayes 3.1.2 [56] and by parsimony searches in PAUP* [57]. For the Bayesian analyses we selected the best model for each partition using Mr Model test [58]. We partitioned the nucleotide data by (a) genomes (2: mtDNA vs. nDNA) (b) genes (5: COI, EF-1α, CAD, ArgK, 28S), (c) positions per genome (7: 1st, 2nd and 3rd positions in mtDNA vs. 1st, 2nd and 3rd positions in the combined protein coding DNA data, and 28S) and (d) gene specific positions (13: each codon position for each protein coding gene, and 28S; see Additional file 5: Table S3). Most partitions had a GTR + I + Γ model selected by AIC, with only COI 3rd positions optimizing a GTR + Γ model (Additional file 5: Table S3). Amino acid translated data from the four protein encoding genes were subject to estimation of gene specific models (mixed evolutionary models). Fifty million generations were run on a Titan cluster of 8 CPUs hosted at, with sampling every 1000 generation. The level of convergence from two parallel runs was inspected in live views of likelihoods for the two runs.

Parsimony analyses consisted of 2,000 heuristic searches with 30 random additions and TBR swapping for each search. Node support was estimated by 200 bootstrap replicates of 20 random addition replicates each. Gene specific contribution to each node was measured by partitioned Bremer support [59]. Ecological character transformation was traced in Mesquite [60] using the parsimony criterion, with test of correlation between fungus farming, gregariousness and inbreeding using Pagel’s [26] modified test of independence between traits. Statistical distributions were based on 20 maximum likelihood searches of 100 data simulations.

Dating of nodes

We estimated divergence times in the software BEAST [61], with input files generated in the BEAST module BEAUti (Additional file 6: Table S4, Additional file 7: Table S5). Data were divided into five partitions consisting of mtDNA1 + 2, mtDNA3, nuc1 + 2, nuc3, 28S. A Yule speciation process birth rate was implemented with a uniform distribution between 1 and 1,000. The Beast tree was calibrated by dating nodes with a relatively precise fossil date, hence using a normally distributed age for these calibration points (see below). Each analysis ran for 20 million generations, with a total of 4,000 trees sampled, deleting the first 2,000 trees as burn-in (Additional file 6: Table S4, Additional file 7: Table S5). The analyses were replicated once with the settings changed according to suggestions from BEAST: scale factor = 0.822 (default 0.75); window-size = 2.0 (default 1.0).

The oldest known Scolytinae fossil Cylindrobrotus pectinatus is possibly from early Aptian Lebanese amber [23]. This species cannot be assigned to any scolytine genus or tribe and likely is indicative of the most ancient form of Scolytinae, predating Scolytinae as currently defined [1]. Burmese amber includes the slightly younger (100 Ma) scolytine fossil in the extant genus Microborus[22] and thus indicates a more exact minimum age for Scolytinae. The fossil fauna of the advanced weevils is relatively young in this respect. The oldest fossil currently known is Ararioerhinus (Anthonomini) from the mid-Cretaceous Santana formation about 112-116 Ma [24]. Previous phylogenetic analyses of weevils has demonstrated an extremely narrow time window for the radiation of the advanced weevils, with Brachycerinae and Entiminae (here represented by Polydrusus) marginally older than the more advanced weevils such as Molytinae, Baridinae and Scolytinae [13, 45, 62]. Thus it is highly unlikely that Scolytinae and other advanced weevils are much older than the 120 myr indicated by the lower Aptian scolytine fossil.

Because there is a conflict between the age of the oldest scolytine and other advanced weevil fossils, two different analyses were made. In the first analysis (A) we used two fossil calibrations, including Scolytinae at 100 Ma (st. dev. = 1) and the advanced weevils at 112 Ma, excluding the Entiminae taxon Polydrusus. In the second analysis (B) we allowed for an older age of Scolytinae corresponding to the lower Aptian age of Cylindrobrotus at 120 Ma, the age of 100 Ma for the oldest node that includes Microborus (excluding the tribe Scolytini according to the amino acid coded data and supported by a previous study [13]), and we included the oldest known dryocoetine fossil that is likely ancestral to Dryocoetes (at the node subtending D. alni). The latter fossil was taken from Dominican amber with an approximate age of 30 Ma as one of very few reliable fossils that fits with our taxon sampling (see also Additional file 6: Table S4, Additional file 7: Table S5).


  1. Jordal BH: Scolytinae. Handbook of Zoology, Band IV Arthropoda: Insecta Part 38: Coleoptera, Beetles, V. Edited by: Beutel R, Leschen R. 2012, Berlin: deGruyter Press, in press

    Google Scholar 

  2. Oberprieler RG, Marvaldi AE, Anderson RS: Weevils, weevils, weevils everywhere. Zootaxa. 2007, 1668: 491-520.

    Google Scholar 

  3. Browne FG: Some aspects of host selection among ambrosia beetles in the humid tropics of south-east Asia. Malayan Forestry. 1958, 21: 164-182.

    Google Scholar 

  4. Browne FG: The biology of Malayan Scolytidae and Platypodidae. Malayan Forest Records. 1961, 22: 1-255.

    Google Scholar 

  5. Schedl KE: Breeding habits of arboricole insects in Central Africa. Xth International Congress of Entomology, Montreal, Proceedings. 1956, 183-197.

    Google Scholar 

  6. Scott JJ, Oh D-C, Yuceer MC, Klepzig KD, Clardy J, Currie CR: Bacterial protection of beetle-fungus mutualism. Science. 2008, 322 (5898): 63-10.1126/science.1160423.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  7. Hulcr J, Adams A, Raffa K, Hofstetter R, Klepzig K, Currie C: Presence and diversity of streptomyces in dendroctonus and sympatric bark beetle galleries across north america. Microbial Ecology. 2011, 61 (4): 759-768. 10.1007/s00248-010-9797-0.

    Article  PubMed  Google Scholar 

  8. Beaver RA: Insect - fungus relationships in the bark and ambrosia beetles. Insect - fungus interactions. Edited by: Wilding N, Collins NM, Hammond PM, Webber JF. 1989, London: Academic, 121-143.

    Chapter  Google Scholar 

  9. Hulcr J, Dunn RR: The sudden emergence of pathogenicity in insect–fungus symbioses threatens naive forest ecosystems. Proceedings of the Royal Society B: Biological Sciences. 2011, 278: 2866-2873. 10.1098/rspb.2011.1130.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Aanen DK, Eggleton P, Rouland-Lefevre C, Guldberg-Frà , slev T, Rosendahl S, Boomsma JJ: The evolution of fungus-growing termites and their mutualistic fungal symbionts. Proceedings of the National Academy of Sciences. 2002, 99 (23): 14887-14892. 10.1073/pnas.222313099.

    Article  CAS  Google Scholar 

  11. Mueller UG, Rehner SA, Schultz TR: The evolution of agriculture in ants. Science. 1998, 281: 2034-2038.

    Article  PubMed  CAS  Google Scholar 

  12. Schultz TR, Brady SG: Major evolutionary transitions in ant agriculture. Proceedings of the National Academy of Sciences. 2008, 105 (14): 5435-5440. 10.1073/pnas.0711024105.

    Article  CAS  Google Scholar 

  13. Jordal BH, Sequeira AS, Cognato AI: The age and phylogeny of wood boring weevils and the origin of subsociality. Molecular Phylogenetics and Evolution. 2011, 59: 708-724. 10.1016/j.ympev.2011.03.016.

    Article  PubMed  Google Scholar 

  14. Jordal BH: Platypodinae. Handbook of Zoology, Band IV Arthropoda: Insecta Part 38: Coleoptera, Beetles, V. Edited by: Beutel R, Leschen R. 2012, Berlin: deGruyter Press, in press

    Google Scholar 

  15. Kok LT, Norris DM, Chu HM: Sterol metabolism as a basis for a mutualistic symbiosis. Nature. 1970, 225: 661-662.

    Article  PubMed  CAS  Google Scholar 

  16. Bentz BJ, Six DL: Ergosterol content of fungi associated with Dendroctonus ponderosae and Dendroctonus rufipennis (Coleoptera: Curculionidae, Scolytinae). Annals of the Entomological Society of America. 2006, 99 (2): 189-194. 10.1603/0013-8746(2006)099[0189:ECOFAW]2.0.CO;2.

    Article  CAS  Google Scholar 

  17. Six DL, Paine TD: Effects of mycangial fungi and host tree species on progeny survival and emergence of Dendroctonus ponderosae (Coleoptera: Scolytidae). Annals of the Entomological Society of America. 1998, 27 (6): 1393-1401.

    Google Scholar 

  18. Farrell BD, Sequeira A, O’Meara B, Normark BB, Chung J, Jordal B: The evolution of agriculture in beetles (Curculionidae: Scolytinae and Platypodinae). Evolution. 2001, 55 (10): 2011-2027.

    Article  PubMed  CAS  Google Scholar 

  19. Kirkendall LR: Ecology and evolution of biased sex ratios in bark and ambrosia beetles. Evolution and diversity of sex ratio: insects and mites. Edited by: Wrensch DL, Ebbert MA. 1993, New York: Chapman and Hall, 235-345.

    Chapter  Google Scholar 

  20. Kirkendall LR, Kent DS, Raffa KA: Interactions between males, females and offspring in bark and ambrosia beetles: the significance of living in tunnels for the evolution of social behaviour. Social behaviour in insects and arachnids. Edited by: Crespi BJ, Choe JC. 1997, Cambridge: Cambridge University Press, 181-215.

    Google Scholar 

  21. Kirkendall LR: The evolution of mating systems in bark and ambrosia beetles (Coleoptera: Scolytidae and Platypotidae). Zoological Journal of the Linnean Society. 1983, 77: 293-352. 10.1111/j.1096-3642.1983.tb00858.x.

    Article  Google Scholar 

  22. Cognato AI, Grimaldi D: 100 million years of morphological conservation in bark beetles (Coleoptera: Curculionidae: Scolytinae). Systematic Entomology. 2009, 34: 93-100. 10.1111/j.1365-3113.2008.00441.x.

    Article  Google Scholar 

  23. Kirejtshuk AG, Azar D, Beaver RA, Mandelshtam MY, Nel A: The most ancient bark beetle known: a new tribe, genus and species from Lebanese amber (Coleoptera, Curculionidae, Scolytinae). Systematic Entomology. 2009, 34: 101-112. 10.1111/j.1365-3113.2008.00442.x.

    Article  Google Scholar 

  24. Santos MFDEA, Mermudes JRM, Fonseca VMMD: A specimen of Curculioninae (Curculionidae, Coleoptera) from the Lower Cretaceous, Araripe Basin, north-eastern Brazil. Palaeontology. 2011, 54 (4): 807-814. 10.1111/j.1475-4983.2011.01057.x.

    Article  Google Scholar 

  25. Wood SL, Bright D: A catalog of Scolytidae and Platypodidae (Coleoptera). Part 2: taxonomic index. Great Basin Naturalist Memoirs. 1992, 13: 1-1553.

    Google Scholar 

  26. Pagel M: Detecting correlated evolution on phylogenies: a general method for the comparative analysis of discrete characters. Proceedings of the Royal Society of London series B. 1994, 255: 37-45. 10.1098/rspb.1994.0006.

    Article  Google Scholar 

  27. Mueller UG, Gerardo N: Fungus farming insects: multiple origins and diverse evolutionary histories. Proceedings of the National Acedemy of Science USA. 2002, 99 (24): 15247-15249. 10.1073/pnas.242594799.

    Article  CAS  Google Scholar 

  28. Mueller UG, Gerardo NM, Aanen DK, Six DL, Schultz TR: The evolution of agriculture in insects. Annual Review of Ecology, Evolution, and Systematics. 2005, 36 (1): 563-595. 10.1146/annurev.ecolsys.36.102003.152626.

    Article  Google Scholar 

  29. Smith SM, Cognato AI: A taxonomic revision of Camptocerus Dejean (Coleoptera: Curculionidae: Scolytinae). Insecta Mundi. 2010, 148: 1-88.

    Google Scholar 

  30. Alamouti SM, Tsui CKM, Breuil C: Multigene phylogeny of filamentous ambrosia fungi associated with ambrosia and bark beetles. Mycological Research. 2009, 113: 822-835. 10.1016/j.mycres.2009.03.003.

    Article  CAS  Google Scholar 

  31. Beaver RA: Host specificity of temperate and tropical animals. Nature. 1979, 281: 139-141. 10.1038/281139a0.

    Article  Google Scholar 

  32. Couvreur T, Forest F, Baker W: Origin and global diversification patterns of tropical rain forests: inferences from a complete genus-level phylogeny of palms. BMC Biology. 2011, 9 (1): 1-12. 10.1186/1741-7007-9-1.

    Article  Google Scholar 

  33. Merckx V, Chatrou L, Lemaire B, Sainge M, Huysmans S, Smets E: Diversification of myco-heterotrophic angiosperms: evidence from Burmanniaceae. BMC Evolutionary Biology. 2008, 8 (1): 1-16. 10.1186/1471-2148-8-1.

    Article  Google Scholar 

  34. Graham A: The age and diversification of terrestrial new world ecosystems through Cretaceous and Cenozoic time. American Journal of Botany. 2011, 98 (3): 336-351. 10.3732/ajb.1000353.

    Article  PubMed  Google Scholar 

  35. Morley RJ: Origin and Evolution of Tropical Rain Forests. 2000, New York: Wiley

    Google Scholar 

  36. Zachos J, Pagani M, Sloan L, Thomas E, Billups K: Trends, rhythms, and aberrations in global climate 65 Ma to present. Science. 2001, 292 (5517): 686-693. 10.1126/science.1059412.

    Article  PubMed  CAS  Google Scholar 

  37. Moreau CS, Bell CD, Vila R, Archibald SB, Pierce NE: Phylogeny of the Ants: Diversification in the Age of Angiosperms. Science. 2006, 312 (5770): 101-104. 10.1126/science.1124891.

    Article  PubMed  CAS  Google Scholar 

  38. Vieites DR, Min M-S, Wake DB: Rapid diversification and dispersal during periods of global warming by plethodontid salamanders. Proceedings of the National Academy of Sciences. 2007, 104 (50): 19903-19907. 10.1073/pnas.0705056104.

    Article  CAS  Google Scholar 

  39. Bininda-Emonds ORP, Cardillo M, Jones KE, MacPhee RDE, Beck RMD, Grenyer R, Price SA, Vos RA, Gittleman JL, Purvis A: The delayed rise of present-day mammals. Nature. 2007, 446 (7135): 507-512. 10.1038/nature05634.

    Article  PubMed  CAS  Google Scholar 

  40. Ericson PGP, Anderson CL, Britton T, Elzanowski A, Johansson US, Källersjö M, Ohlson JI, Parsons TJ, Zuccon D, Mayr G: Diversification of Neoaves: integration of molecular sequence data and fossils. Biology Letters. 2006, 2 (4): 543-547. 10.1098/rsbl.2006.0523.

    Article  PubMed  PubMed Central  Google Scholar 

  41. Ward PS, Brady SG, Fisher BL, Schultz TR: Phylogeny and biogeography of dolichoderine ants: effects of data partitioning and relict taxa on historical inference. Systematic Biology. 2010, 59 (3): 342-362. 10.1093/sysbio/syq012.

    Article  PubMed  CAS  Google Scholar 

  42. Archibald SB, Bossert WH, Greenwood DR, Farrell BD: Seasonality, the latitudinal gradient of diversity, and eocene insects. Paleobiology. 2010, 36 (3): 374-398. 10.1666/09021.1.

    Article  Google Scholar 

  43. Harrington GJ, Eberle J, Le-Page BA, Dawson M, Hutchison JH: Arctic plant diversity in the early eocene greenhouse. Proceedings of the Royal Society B: Biological Sciences. 2011, 279: 1515-1521.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Kuschel G, Leschen RAB, Zimmerman EC: Platypodidae under scrutiny. Invertebrate Taxonomy. 2000, 14 (6): 771-805. 10.1071/IT00024.

    Article  Google Scholar 

  45. McKenna DD, Sequeira AS, Marvaldi AE, Farrell BD: Temporal lags and overlap in the diversification of weevils and flowering plants. Proceedings of the National Academy of Science USA. 2009, 106 (17): 7083-7088. 10.1073/pnas.0810618106.

    Article  CAS  Google Scholar 

  46. Davis Charles C, Webb Campbell O, Wurdack Kenneth J, Jaramillo Carlos A, Donoghue Michael J: Explosive radiation of malpighiales supports a mid‐cretaceous origin of modern tropical rain forests. The American Naturalist. 2005, 65 (3): E36-E65.

    Article  Google Scholar 

  47. Nobre T, Koné NA, Konaté S, Linsenmair KE, Aanen DK: Dating the fungus-growing termites’ mutualism shows a mixture between ancient codiversification and recent symbiont dispersal across divergent hosts. Molecular Ecology. 2011, 20 (12): 2619-2627. 10.1111/j.1365-294X.2011.05090.x.

    Article  PubMed  CAS  Google Scholar 

  48. Aanen DK, Eggleton P: Fungus-growing termites originated in african rain forest. Current biology. 2005, 15 (9): 851-855. 10.1016/j.cub.2005.03.043.

    Article  PubMed  CAS  Google Scholar 

  49. Jordal BH: Elongation Factor 1α resolves the monophyly of the haplodiploid ambrosia beetles Xyleborini (Coleoptera: Curculionidae). Insect Molecular Biology. 2002, 11: 453-465. 10.1046/j.1365-2583.2002.00354.x.

    Article  PubMed  CAS  Google Scholar 

  50. Jordal BH, Beaver RA, Normark BB, Farrell BD: Extraordinary sex ratios and the evolution of male neoteny in sib-mating Ozopemon beetles. Biological Journal of the Linnean Society. 2002, 75: 353-360. 10.1111/j.1095-8312.2002.tb02076.x.

    Article  Google Scholar 

  51. Jordal BH, Normark BB, Farrell BD: Evolutionary radiation of an inbreeding haplodiploid beetle lineage (Curculionidae, Scolytinae). Biological Journal of the Linnean Society. 2000, 71 (3): 483-499. 10.1111/j.1095-8312.2000.tb01270.x.

    Article  Google Scholar 

  52. Normark BB, Jordal BH, Farrell BD: Origin of a haplodiploid beetle lineage. Proceedings of the Royal Society of London serie B. 1999, 266: 2253-2259. 10.1098/rspb.1999.0916.

    Article  Google Scholar 

  53. Bright DE, Poinar GO: Scolytidae and Platypodidae (Coleoptera) from Dominican Republic amber. Annals of the Entomological Society of America. 1994, 87 (2): 170-194.

    Article  Google Scholar 

  54. Edgar RC: MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Research. 2004, 32 (5): 1792-1797. 10.1093/nar/gkh340.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  55. Jordal BH, Gillespie JJ, Cognato AI: Secondary structure alignment and direct optimization of 28S rDNA sequences provide limited phylogenetic resolution in bark and ambrosia beetles (Curculionidae: Scolytinae). Zoologica Scripta. 2008, 37: 1-14.

    Google Scholar 

  56. Huelsenbeck JP, Ronquist F: MRBAYES: Bayesian inference of phylogeny. Biometrics. 2001, 17 (8): 754-755.

    CAS  Google Scholar 

  57. Swofford D: PAUP*. Phylogenetic Analysis Using Parsimony (*and other methods) version 4. 2002, Massachusetts: Sinauer Associates, Sunderland

    Google Scholar 

  58. Nylander JAA: Mr Model test. Version 2. Computer program distributed by the author. 2004, Uppsala University: Evolutionary Biology Centre, 2

    Google Scholar 

  59. Baker RH, DeSalle R: Multiple sources of character information and the phylogeny of Hawaiian drosophilids. Systematic Biology. 1997, 46: 645-673.

    Google Scholar 

  60. Maddison WP, Maddison DR: Mesquite: a modular system for evolutionary analysis. Version 2.74. Computer program distributed by the author, 2011,

  61. Drummond AJ, Rambaut A: BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evolutionary Biology. 2007, 7: 214-214. 10.1186/1471-2148-7-214.

    Article  PubMed  PubMed Central  Google Scholar 

  62. Marvaldi AE, Sequeira AS, O’Brien CW, Farrell BD: Molecular and morphological phylogenetics of weevils (Coleoptera, Curculionoidea): do niche shifts accompany diversification?. Systematic Biology. 2002, 51 (5): 761-785. 10.1080/10635150290102465.

    Article  PubMed  Google Scholar 

Download references


We thank Debra and Sarah Smith for reviewing a previous version of this manuscript. Funded in part by the National Science Foundation (DEB-0328920) and the National Geographic Society (7922-05) to AIC, and by the Norwegian Research Council no. 214232/F20 to BHJ.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Bjarte H Jordal.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

BHJ designed the study, collected and assembled beetle samples, carried out the molecular phylogenetic analyses and drafted the manuscript. AIC participated in the overall design of the study, collected some of the samples, and contributed to drafting the manuscript. Both authors have read and approved the final manuscript.

Electronic supplementary material


Additional file 1: Figure S1. Partitioned Bremer support values for each of the five gene fragments included. Their relative support was estimated on two tree topologies: the 84-taxa parsimony topology and the Bayesian 7-partitions topology pruned to 84 taxa. (PDF 271 KB)


Additional file 2: Table S1. Taxon sampling and their respective accession numbers in GenBank (missing data denoted by ’-’). (DOCX 67 KB)


Additional file 3: Figure S2. Beast tree topology showing each origin of fungus farming in thick blue branches, and each origin of regular inbreeding marked by a black dot. Below is the Zachos curve of temperature variation during the Cenozoic. (PDF 357 KB)


Additional file 4: Table S2. Character states for the taxa sampled, with the absence (0) or presence (1) of fungus farming, gregarious larval feeding and regular inbreeding. (DOCX 41 KB)

Additional file 5: Table S3. Properties of each partition, estimated from 84 taxa with no missing data. (DOCX 18 KB)

Additional file 6: Table S4. Xml file generated by Beauti and analyzed in Beast (analysis A). (TXT 798 KB)

Additional file 7: Table S5. Xml file generated by Beauti and analyzed in Beast (analysis B). (TXT 799 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and Permissions

About this article

Cite this article

Jordal, B.H., Cognato, A.I. Molecular phylogeny of bark and ambrosia beetles reveals multiple origins of fungus farming during periods of global warming. BMC Evol Biol 12, 133 (2012).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Agriculture
  • Symbiosis
  • Molecular phylogeny
  • Bark and ambrosia beetles
  • Evolutionary origins