Cretaceous environmental changes led to high extinction rates in a hyperdiverse beetle family
BMC Evolutionary Biology volume 14, Article number: 220 (2014)
As attested by the fossil record, Cretaceous environmental changes have significantly impacted the diversification dynamics of several groups of organisms. A major biome turnover that occurred during this period was the rise of angiosperms starting ca. 125 million years ago. Though there is evidence that the latter promoted the diversification of phytophagous insects, the response of other insect groups to Cretaceous environmental changes is still largely unknown. To gain novel insights on this issue, we assess the diversification dynamics of a hyperdiverse family of detritivorous beetles (Tenebrionidae) using molecular dating and diversification analyses.
Age estimates reveal an origin after the Triassic-Jurassic mass extinction (older than previously thought), followed by the diversification of major lineages during Pangaean and Gondwanan breakups. Dating analyses indicate that arid-adapted species diversified early, while most of the lineages that are adapted to more humid conditions diversified much later. Contrary to other insect groups, we found no support for a positive shift in diversification rates during the Cretaceous; instead there is evidence for an 8.5-fold increase in extinction rates that was not compensated by a joint increase in speciation rates.
We hypothesize that this pattern is better explained by the concomitant reduction of arid environments starting in the mid-Cretaceous, which likely negatively impacted the diversification of arid-adapted species that were predominant at that time.
In the fossil record of marine and terrestrial taxa, the Cretaceous Period (145–66 Million years ago (Ma)) is usually considered a time of major reorganization and modernization of ecosystems characterized by the extinction of groups that were formerly dominant and the appearance and subsequent diversification of new groups ,. The climate during this period is often considered a textbook example of a greenhouse world , because of the high average global temperatures that lasted until about 70 Ma. Still, this period had its share of both warming and cooling events, and was notably characterized by an extensive tropical ecosystem  and an important reduction of arid zones that started about 110–120 Ma ,. Of particular interest, a period called the Cretaceous terrestrial revolution (KTR; also referred as the angiosperm revolution ,), saw the explosive radiation of angiosperms which rose from 0 to 80% of floral composition between 125–90 Ma . The KTR likely provided new ecological and evolutionary opportunities for insects - and for several vertebrate clades that underwent major diversification during this period ,. Although it remains challenging to investigate the diversification dynamics of species-rich groups in deep time, ancient clades represent important avenues to assess the impacts of past environmental changes such as those we now face (e.g. see Roelants et al. ). However, the picture from the fossil record is incomplete and biased for most species-rich groups such as insects . Thankfully, recent developments of cutting-edge methods that allow estimation of variation in diversification rates among lineages , provide new tools for assessing the diversification dynamics of groups with a poor fossil record. However, very few studies have addressed the effects of Cretaceous environmental changes, in particular the KTR and the Cretaceous-Palaeogene (K-Pg) mass extinction, using these new methods (but see Lloyd et al.  and Meredith et al. ). In birds and mammals an increase of diversification rates was revealed, which can be accounted for by high speciation and low extinction rates ,,. In contrast, the diversification rates in dinosaurs remained constant through the KTR .
In insects, the explosive radiation of angiosperms is thought to have provided diversification opportunities for pollinators, leaf-mining flies, as well as butterflies and moths . Though numerous dated phylogenies are available, little information can be gleaned on the diversification dynamics of most major groups or families during the Mesozoic. The most notable exceptions are a study on Lepidoptera , which highlighted three increases in diversification rates during the Cretaceous; several studies on ants , that revealed several diversification bursts about 100 Ma; and a study on carpenter bees , which provided support for a K-Pg mass extinction event. In Coleoptera, Hunt et al.  also postulated that the diversification of Coleoptera is better explained by the persistence of old lineages (see also Wang et al. ) that diversified in multiple niches, rather than high diversification rates or the effect of the Cretaceous rise of angiosperms. That said, at the time this study was conducted, methods that explicitly incorporate birth-death models allowing diversification rates to vary were not widely available, requiring authors to rely instead on sister-clade comparisons.
To provide novel insights on insect diversification dynamics, and especially on insect groups that do not feed on live plant tissues, we investigated how Cretaceous changes impacted the diversification of a hyperdiverse beetle family, the Tenebrionidae (darkling beetles). This family currently encompasses ca. 20,000 described species worldwide, and ranks seventh in terms of species richness within the Coleoptera . Adults and larvae are mostly saprophagous, feeding on decaying vegetation, but some species are also mycetophagous and predatory larvae are known in a few tribes . Most of the extant tenebrionid species are either distributed in arid/semiarid regions (e.g., members of the large subfamily Pimeliinae) or in subtropical/tropical forested regions ,. The habitat preference is so marked in this family that shifts in their abundance and range have been used as indicators of climate change . The tenebrionid fossil record highlights an age as old as the Middle Jurassic for the family , and indicates that Mesozoic Pangaean and Gondwanan breakups may be reflected in the disjunct distribution patterns observed in several groups such as the Pimeliinae or the Adeliini ,,. However, due to its relative scarcity, the fossil record alone is not sufficient to provide an accurate temporal framework for the family.
To gain a better understanding of the age of specific divergence events, we conducted molecular dating analyses using 21 calibration points combined with conservative priors and calibration settings. We then used the resulting timeframe to conduct diversification analyses to investigate the early tempo and mode of diversification of the family in relation to Cretaceous environmental changes. Postulating that past environmental changes may have contributed to the diversification dynamics of Tenebrionidae, as more than half of extant species are adapted to arid/semi-arid conditions, we also assessed whether the rise of angiosperms also coincided with changes in the diversification dynamics of clades that encompass forest-dwelling groups.
Results and discussion
Mesozoic origin and diversification of a species-rich beetle family
Dating analyses unequivocally support an ancient origin for Tenebrionidae at the end of the Early Jurassic (ca. 180 Ma, Figure 1; Table 1), which significantly predates (about 30 Myrs) previous molecular estimates of tenebrionid ages ,. This earlier date of origin has several consequences on our understanding of tenebrionid origin and evolution against the backdrop of past environmental changes. The first lineages to diversify were a major clade of the paraphyletic subfamily Pimeliinae (ca. 9,000 species; median age of 144.29 Ma, 95% HPD: 116.69-171.77 Ma for the best-fit calibration procedure) and the Lagriinae (ca. 2,000 species; median age of 147.1 Ma, 95% HPD: 132.14-149.53 Ma). For these two groups, extant distribution patterns are consistent with potential vicariance events, as advocated for Adeliini (ca. 400 species; median age of 95.13 Ma, 95% HPD: 82.69-113.95 Ma), which occur today in Australia, Chile, New Caledonia and New Zealand ,. For Pimeliinae, our results support the hypothesis of Matthews et al. , who postulated that the highly disjunct extant distributions of several pimeliine tribes are better explained by an ancient origin and Gondwanan fragmentation. The diversification of these two subfamilies was followed by other subfamilies during the KTR (125 to 80 Ma), such as Alleculinae (median age of 126.83 Ma, 95% HPD: 120.0-130.81 Ma), Stenochiinae (median age of 100.22 Ma, 95% HPD: 88.76-109.9 Ma) and several lineages of Diaperinae and Tenebrioninae.
Contrary to the arid-dwelling Pimeliinae, the other subfamilies are currently more diverse in tropical environments (see Matthews et al.  for details). Thus if we consider that phylogenetic biome conservatism holds in many groups , this suggests that distinct tenebrionid lineages were already adapted to arid or tropical environments before the end of the KTR. Results from character optimization analyses support this hypothesis and suggest that the common ancestor of tenebrionid beetles was adapted to arid environments (Figure 2). Arid environments were extremely widespread in the Early Jurassic ,; though they progressively receded with the break-up of the Pangaea ,, they never completely disappeared, which could explain the persistence of distinct lineages of pimeliine beetles across eastern and western hemispheres . Interestingly, pimeliine beetles are currently nearly absent from Australia, where aridification began in the Miocene and deserts only appeared very recently ,. As underlined by Matthews et al. , this suggests that the group was strictly adapted to arid environments from the beginning, and unable to colonize the Australian continent before its isolation in the middle of the Cenozoic. The only exceptions are the few representatives (nine species) of two plesiotypic tribes (Cnemeplatiini and Vacronini) that were possibly present as coastal sand dune inhabitants before separation from Gondwana . Other groups such as the rain-forest specialists Adeliini  originated about 96 Ma (median age of 95.13 Ma, 95% HPD: 82.69-113.95 Ma), possibly in the warm temperate forests that were widespread in the southern hemisphere at that time . This is also the case for other tropical groups in the subfamily Tenebrioninae (tribes Heleini and Titaeini), whose most recent common ancestor diversified at the beginning of the Late Cretaceous, around 100 Ma (median age of 99.83 Ma, 95% HPD: 91.76-111.76 Ma). In some of these groups the phylogenetic biome conservatism is not irreversible, highlighted by the results of character optimizations indicating that several lineages (e.g. Blaptini in Tenebrioninae) became secondarily adapted to arid environments. These secondary shifts likely occurred multiple times during the evolutionary history of darkling beetles, as several subfamilies encompass xerophilic tribes (two tribes out of nine in Lagriinae; 11 tribes out of 29 in Tenebrioninae and five tribes out of 11 in Diaperinae).
Impacts of cretaceous environmental changes on tenebrionid diversification dynamics
The lineages-through-time plots suggested a rapid initial diversification, followed by a slowdown in the accumulation of lineages around 100–103 Ma (Figure 3A, Additional file 1: Figure S1). Similar to Meredith et al. , we used birth-death models, implemented in the TreePar approach, to analyse the time-calibrated trees for the time period spanning 50 Ma to the origin of the group (ca. 180 Ma). We found that a diversification model with varying rates through time was supported by both likelihood ratio tests (P <0.001) and corrected Akaike information criteria compared to a constant birth-death model (Table 2). The model with one shift time (two-rate diversification model) fit the trees better (Table 2). Similar trends were also obtained using randomly sampled trees from the posterior trees of the Bayesian dating analysis with the best-fit calibration procedure (Additional file 2: Tables S1-S2).
The period when diversification (=speciation - extinction) and turnovers (=extinction/speciation) changed was inferred to have been the boundary between the Early and Late Cretaceous in the middle of the KTR (shift time at ca. 103 Ma for the best-fit calibration procedure; Figure 3a, Table 1). We found that from the origin of the family to about 103 Ma, the net diversification rates of tenebrionid beetles (mean of 0.085 lineages Myr–1) were more elevated than the average net diversification rates for all coleopterans estimated in Hunt et al.  (ranging from 0.046 to 0.068 lineages Myr–1). Turnovers were not elevated before the shift time (0.305 and 0.178 for the “Yule crown” and “BD crown”, respectively). When calculating speciation and extinction using the diversification rate and turnover ratio, the speciation rate was 0.122 events/lineages Myr–1 and the extinction rate was 0.037 events/lineages Myr–1 before the shift time (Figure 3A).
After the shift time (103 Ma), turnover rates reached high levels (0.935) while net diversification rates dropped drastically (0.021 lineages/Myr). This important slowdown of net diversification rates (4-fold decrease) resulted from a marked increase of extinction rates (8.5-fold increase, from 0.037 to 0.316 events/lineages Myr–1), which was not compensated by a concurrent increase in speciation rates (2.8-fold increase, from 0.122 to 0.338 events/lineages Myr–1). Thus, we found no support for a positive upturn in diversification rates during the KTR, contrary to the patterns observed for the diversification of Lepidoptera , ants ,, mammals , and birds . This result suggests that the changes in floral communities concomitant with the radiation of Angiosperms  did not significantly foster the diversification of darkling beetles, a pattern in agreement with those observed in several other groups such as dinosaurs  and flies . It is worth highlighting that the decrease of net diversification rates we revealed cannot be attributed to the incompleteness of our sampling, because we would not have otherwise been able to capture the simultaneous increase of speciation rates at that time.
To better understand the potential contribution of major abiotic and biotic factors, we provide an overview of the possible correlates of diversification for tenebrionids on a global scale (Figure 3B). Our results allow us to rule out any role of the K-Pg event on the significant increase in extinction rates, as it postdates the shift of diversification rates by about 38 Myr. On the contrary, we postulate that the drastic reduction of arid habitats that started at the end of the Early Cretaceous specifically impacted the diversification of arid-specialist darkling beetles. In particular the subfamily Pimeliinae, currently the single most speciose tenebrionid subfamily, is very conspicuous in desert ecosystems and is found sister to all other tenebrionids, which likely indicates that pimeliines were also predominant during the long-lasting arid episode that persisted between the Late Jurassic and the Early Cretaceous. Therefore the shift toward humid conditions ca. 110–120 Ma and the contractions of desert ecosystems (Figure 3B) probably precipitated the extinction of multiple pimeliine lineages, leaving behind relict lineages with completely disjunct distributions. At least five instances of the latter are known  for tribes Caenocrypticini (Namib desert in southern Africa and coastal deserts of Peru and Chile), Cossyphodini (Africa, India and South America), Evaniosomini (Namib desert and South America), Elenophorini (western Mediterranean and southern South America) and Vacronini (North American and central Australian). It is also important to underline that the vast majority of Pimeliinae are flightless ; reduced dispersal abilities may have increased their sensitivity to the overall reduction of available habitats at a broader scale , limiting their ability to track climate change . While rapid habitat loss may have accelerated the extinction rate of this family dominated at the time by arid-adapted species, one could also imagine that their low vagility may have promoted the speciation of geographically isolated relictual populations through reduced migration between populations at a more local scale .
We may hypothesize that high environmental temperatures during the Jurassic and Cretaceous increased rates of biological processes, shortening generation time and thus possibly increasing evolutionary speed , leading to relatively high speciation rates. This role of temperature as a key driver of biodiversity is widely acknowledged for large-scale studies in both marine and terrestrial groups ,-. As tenebrionid beetles are more adapted to arid or tropical regions ,,, they may have benefited overall from the favourable warm climatic conditions that persisted during most of the KTR. After the inferred shift time, the concomitant increase in speciation rate can be attributed to the tropical-adapted lineages. These groups probably took advantage of the expansion of warm tropical and subtropical-forested biomes between the Early and the Late Cretaceous , which likely provided them with new ecological opportunities. This hypothesis is also supported by the fact that major tropical groups started their diversification shortly before (Alleculinae) or after (Lagriinae: Adeliini; Stenochiinae; Tenebrioninae: Heleini, Titaeini) the inferred time of environmental shift.
Overall the global diversification pattern that was inferred for tenebrionid beetles is quite intricate (slowdown of diversification despite an increase of speciation rates) and illustrates the importance of integrating knowledge of the biology of taxa whenever possible. Though we cannot exclude a possible role of the KTR (for the tropical-adapted species), we postulate that the reduction of arid environments that started 110–120 Ma is a key factor that explains the slowdown in diversification dynamics for the family.
We used a molecular dataset of eight gene fragments (four mitochondrial and four nuclear) encompassing 404 species (from the study of Kergoat et al. ; see also Additional file 3: Table S3). This dataset was first introduced in a study that has mostly focused on discussing tenebrionid relationships . Within the family Tenebrionidae this dataset comprises 250 species (150 of which were collected and sequenced by the authors) belonging to 98 distinct genera. The sampled subfamilies (seven out of nine) represent more than 99% of the tenebrionid species diversity. The sampled tribes (37 tribes out of 96) also encompass most (ca. 77%) of the tenebrionid generic diversity (see Additional file 4: Table S4). For calibration purposes, this dataset includes representatives of several tenebrionid genera (Bolitophagus, Gonocephalum, Isomira, Lorelus, Platydema, Pentaphyllus and Tribolium), for which several unambiguous fossil representatives are known . One hundred fifty four species (from 42 beetle families) are also used as outgroups. The rationale here is to sample representatives of deep lineages in order to maximize the number of available calibration points for molecular dating analyses.
For all dating analyses, we used a Bayesian relaxed-clock (BRC) approach as implemented in BEAST v1.7.5 . Constraints on clade ages were enforced using sixteen fossil calibrations and five geological calibrations. Fossil constraints were defined following Parham et al.  and Sauquet et al.  and used to specify minimum ages at nodes (see Additional file 5 for detailed information on each fossil). For all fossil constraints we used the youngest age of the geological stage (using the latest geological time scale from Gradstein et al. ) or the youngest age known for a specific amber deposit as minimum age for the calibration to avoid any overestimation of fossil age ,. For geological calibrations, we relied on taxa endemic from particular islands of the Canary archipelago, which are commonly used for molecular dating analyses ,. Because the geological history of these volcanic islands is well documented , it allows the specification of maximum age constraints for each island. We used conservative estimates for all islands, namely: 20 Ma for Fuerteventura; 15.5 Ma for Lanzarote; 16 Ma for Gran Canaria; 12 Ma for La Gomera; 11.6 Ma for Tenerife; 2 Ma for La Palma and 1.1 Ma for El Hierro . Given that the magmatic event represented by the Canary Islands started 24 Ma, we set a maximum age of 24.0 Myr (corresponding to the age of the two oldest islands) for the node leading to the clade encompassing the 13 Pimelia endemic to the Canary Islands. A maximum age of 16.0 Myr was set for the clade including the three subspecies of P. sparsa that are only found in Gran Canaria. A maximum age of 11.6 Myr was set for the clade encompassing P. canariensis, P. radula ascendens, P. r. radula and P. r. granulata, which are all endemic to Tenerife. A maximum age constraint of 12.0 Myr was set for the clade encompassing P. laevigata laevigata (endemic to La Palma), P. l. validipes (endemic to La Gomera) and P. l. costipennis (endemic to El Hierro). Finally a maximum age constraint of 1.1 Myr was set for the clade encompassing P. l. validipes (endemic to La Gomera) and P. l. costipennis (endemic to El Hierro). Finally, to optimize the exploration of parameter space and avoid overestimation of root age we specify a conservative maximum age of 265 million years (Myr) for the root, which significantly predates any known fossil occurrence for the clade that encompasses the sampled superfamilies .
BEAST analyses were implemented with partitioned relaxed-clock models (one uncorrelated lognormal clock per gene). In order to limit the numbers of parameter to estimate, we used a guide tree that corresponds to the best maximum likelihood tree inferred with the same dataset in the study of Kergoat et al. . Fossil constraints were either placed on crown or stem nodes. Fossil constraints put on stems are generally more conservative, especially if the sampling of the group is not enough to provide a relevant representation of the clade diversity ,. On the other hand, fossil constraints put on stems can also result in severe biases and underestimate the age of a group if the sister clade is too distantly related ,. Since our taxon sampling is far from being complete and given the difficulty to assign a calibration to a precise node, we adopt a conservative way in which all node calibrations were enforced using uniform distributions (hard-bound constraints ,,). To account for the fact that our trees describe inter-specific relationships, we also perform distinct analyses with two distinct tree speciation priors: birth-death (BD) and Yule.
For each of the four distinct calibration procedures, two distinct runs were carried out with 50 million generations and trees sampled every 5,000 generations (10,000 trees were sampled for each run). BEAST .xml files were also modified to implement the path-sampling procedure , which allows an unbiased approximation of the marginal likelihood of runs . The latter is critical to properly infer the Bayes factors that are used to sort among competing calibration procedures . We used a conservative burn-in-period of 12.5 million generations per run. Post burn-in trees from the two distinct runs (7,500 trees for each run) were further combined using the LogCombiner module of BEAST. Convergence of runs was assessed graphically under Tracer v1.5  and by examining the effective sample size (ESS) of parameters. Bayes factors  were then estimated using the log files of the four distinct calibration procedures, using scripts detailed in Baele et al. . Convergence was indicated by ESS of parameters >200 for the post burn-in trees. Dating analyses based on stem calibrations recovered several age estimates that were not consistent with the fossil record (see Additional file 6: Table S5), and were thus discarded. Out of the two crown calibrations, the “Yule crown” calibration procedure was recovered as the best-fit calibration procedure by the Bayes factor comparison (Additional file 7: Table S6). A median age of 180.05 Ma (95% HPD: 169.56-191.62; “Yule crown”) was estimated for the Tenebrionidae (178.04 Ma for the “BD crown”)(see also Additional file 8: Figure S2, Additional file 9: Figure S3).
Diversification analyses were conducted on chronograms resulting from BRC calibration procedures. Analyses were realized under the R environment software implementing the TreePar package. To visualize the tempo and mode of diversification of the group, we first reconstructed lineages-through-time plots. We then assessed whether diversification rates remained constant during the early evolutionary history of darkling beetles. Similar to the way Meredith et al.  examined the diversification rates of mammals and the impact of the K-Pg extinction event, we focused on the time period between the estimated origin of the family and 50 Ma (early Eocene). We did not examine the tempo or mode of diversification of the group after this date because our chronograms included only a small fraction of the standing diversity of crown-group tenebrionids that were likely extant in the Cenozoic. We are aware of potential biases in diversification rate analyses that use incomplete phylogenies ,. However, our aim was not to reveal a precise diversification pattern but instead to get a broad sense of the early diversification dynamics. Several previous studies on diversification rates for highly incomplete phylogenies have revealed interesting evolutionary patterns despite methodological limitations of their trait-dependent analyses (e.g., Goldberg et al.  and Hugall & Stuart-Fox ), which are even more sensitive to the missing lineages than are those used in the current study .
The TreePar package was used to assess speciation and extinction rates through time. This method relaxes the assumption of constant rates by allowing rates to change at specific points in time . Such a model allows for the detection of rapid changes in speciation and extinction rates due to environmental factors like the K-Pg mass extinction. We employed the `bd.shifts.optim’ function that allows for estimating discrete changes in speciation and extinction rates and mass extinction events in under-sampled phylogenies . Going backward in time, it estimates the maximum likelihood speciation and extinction rates together with the rate shift times t = (t1,t2,…,tn) in a phylogeny. At each time t, the rates are allowed to change and the species may undergo a mass extinction event. TreePar analyses were run with the following settings: start =50, end = crown age estimated by dating analyses, grid =1 Myr, four possible shift times were tested, and posdiv = FALSE to allow the diversification rate to be negative (i.e. allows for periods of declining diversity). As our taxon sampling did not include all lineages between the origin of the group and 50 Ma, we adjusted the sampling fraction accordingly. This value was computed using the inferred diversity at 50 Ma divided by the number of lineages we had in the phylogeny. We first estimated a whole diversification rate for Tenebrionidae using the method of moments of Magallón & Sanderson  with three relative extinction rates (ε = 0/0.5/0.9) and the current total diversity of the group (20,000 species). We then calculated the diversity at 50 Ma using the following formula:
We investigated the robustness of our results with the maximum clade credibility tree. We took into account the effect of dating uncertainties by repeating the diversification analyses on 500 and 1,000 random chronograms sampled from the Bayesian posterior distributions. The latter allowed us to estimate the mean of all free parameters and their standard errors. Given the diversification rate (=speciation – extinction) and turnover (=extinction/speciation), we inferred the speciation and extinction rate as follows:
Using the method of moments and three relative extinction rates (ε = 0/0.5/0.9), the tenebrionid diversification rates were estimated at 0.0529/0.0513/0.0434 lineages Myr–1. Given the current diversity and the diversification rates, we inferred the diversity at 50 Ma around 7556/7800/9218 species (see formula above). Knowing the number of sampled lineages, we deduced the following sampling fraction 0.0265/0.0256/0.0217 for our phylogenies at 50 Ma. Thus, we conservatively set the sampling fraction to 0.03 in the TreePar analyses.
Character trait optimizations
Character trait optimization analyses were conducted to investigate the phylogenetic biome conservatism of the sampled tenebrionid lineages. We used as a guide tree the best tree inferred under maximum likelihood from the study of Kergoat et al. . This tree was further modified by removing all outgroup taxa and retaining only one species (randomly chosen among specimens with fewest missing data) per tenebrionid genus. Habitat preferences were subsequently coded for the corresponding 98 genera (Additional file 10: Table S7), using a two-state categorization (“arid/semi-arid” and “other”). Ancestral character state estimations were further carried out under maximum likelihood using a one-parameter Markov k-state model with symmetrical rates , as implemented in Mesquite 2.75 . The support of one state over another (at a given node) was considered as significant if the difference between their log-likelihoods was greater than or equal to 2.0 .
Limitations of the taxon sampling on the analyses
Conclusions about the temporal nature of diversification are dependent upon the quality of the data . Thus there are important methodological points that may affect analyses and results presented here.
An important methodological issue might be the estimates of divergence times, which are highly suspected to introduce error into diversification analyses or distorted conclusions -. For instance, too few fossils, failure to use an appropriate prior distribution, or poor calibration choices can lead to inaccurate and illusory dating results (e.g. ). These problems may be important as all the post-tree analyses rely on the branching times of the group. Here we have assessed the divergence times of tenebrionids using a rigorous analytic pipeline with BRC analyses, all unambiguous described fossils plus external calibrations; we also repeated the analyses to check likelihood congruence and to cross validate our results. Our dating analyses provided similar divergence times for the phylogeny, and diversification rate analyses were congruent.
Incomplete taxon sampling is a potentially severe problem that is difficult to address. First, missing lineages can lead to inaccurate phylogenetic reconstructions, and anomalous branch lengths can, in turn, bias dating analyses. In our study, we tried to recover the densest taxon sampling that restricts these potential effects, and we have taken into account missing lineages in our analyses and interpretation. Second, missing taxon sampling can be a problem for selecting the best-fit model of diversification. Spurious results of diversification rates can be inferred with under-sampled phylogenies . Rabosky  explained that the main issue is not the quality of the taxon sampling, but rather whether a rate is informative if diversification is diversity-dependent. Although it is difficult to detect and incorporate diversity-dependence in the diversification analyses at our level, it could be nonetheless interesting to investigate whether subclades are closer to their diversity limits. A study on the genus Blaps (Tenebrionidae: Blaptini) found no effect of diversity-dependence on its diversification .
In addition, there are inherent limitations for reconstructing ancestral character states, and thus for inferring when a biome shift occurred. This might be less of a problem here, as the evolution of biome association is very clear and includes few shifts. As taxonomic sampling becomes more extensive in phylogenetic studies of tenebrionids, there will be opportunities for more thorough comparisons of diversification rate variation and biome evolution.
GJK designed the study, participated in field sampling, conducted analyses, and wrote the paper. PB and JLA wrote the paper. ALC participated in field sampling and performed experiments. HJ designed and participated in the sampling. RJ-Z and GG performed experiments. LS participated in field sampling and performed the identification of sampled specimens. FLC designed the study, participated in field sampling, performed experiments, conducted analyses, and wrote the paper. All authors commented on the manuscript. All authors read and approved the final manuscript.
Lloyd GT, Davis KE, Pisani D, Tarver JE, Ruta M, Sakamoto M, Hone DW, Jennings R, Benton MJ: Dinosaurs and the Cretaceous terrestrial revolution. Proc R Soc B. 2008, 275: 2483-2490. 10.1098/rspb.2008.0715.
Coiffard C, Gomez B, Daviero-Gomez V, Dilcher DL: Rise to dominance of angiosperm pioneers in European Cretaceous environments. Proc Natl Acad Sci U S A. 2012, 109: 20955-20959. 10.1073/pnas.1218633110.
Spicer RA, Chapman JL: Climate change and the evolution of high-latitude terrestrial vegetation and floras. Trends Ecol Evol. 1990, 5: 279-284. 10.1016/0169-5347(90)90081-N.
Vakhrameev VA: Jurassic and Cretaceous floras and climates of the Earth. 1991, Cambridge University Press, Cambridge
Morley RJ: Cretaceous and Tertiary Climate Change and the Past Distribution of Megathermal Rainforests. Tropical Rainforest Responses to Climatic Changes. Edited by: Bush MB, Flenley J. 2007, Praxis Publishing, Chichester, 1-31. 10.1007/978-3-540-48842-2_1.
Hallam A: Continental humid and arid zones during the Jurassic and Cretaceous. Palaeogeogr Palaeoclim Palaeoecol. 1984, 47: 195-223. 10.1016/0031-0182(84)90094-4.
Föllmi KB: Early Cretaceous life, climate and anoxia. Cretaceous Res. 2011, 35: 230-257. 10.1016/j.cretres.2011.12.005.
Labandeira CC: A paleobiologic perspective on plant-insect interactions. Curr Opin Plant Biol. 2013, 16: 414-421. 10.1016/j.pbi.2013.06.003.
Benton MJ: The origins of modern biodiversity on land. Phil Trans R Soc B. 2010, 365: 3667-3679. 10.1098/rstb.2010.0269.
Farrell BD: “Inordinate fondness” explained: why are there so many beetles?. Science. 1998, 281: 555-559. 10.1126/science.281.5376.555.
Moreau CS, Bell CD, Vila R, Archibald SB, Pierce NE: Phylogeny of the ants: diversification in the age of angiosperms. Science. 2006, 312: 101-104. 10.1126/science.1124891.
McKenna DD, Sequeira AS, Marvaldi AE, Farrell BD: Temporal lags and overlap in the diversification of weevils and flowering plants. Proc Natl Acad Sci U S A. 2009, 106: 7083-7088. 10.1073/pnas.0810618106.
Wahlberg N, Wheat CW, Peña C: Timing and patterns in the taxonomic diversification of Lepidoptera (butterflies and moths). PLoS One. 2013, 8: e80875-10.1371/journal.pone.0080875.
Meredith RW, Janečka JE, Gatesy J, Ryder OA, Fisher CA, Teeling EC, Goodbla A, Eizirik E, Simao TLL, Stadler T, Rabosky DL, Honeycutt RL, Flynn JJ, Ingram CM, Steiner C, Williams TL, Robinson TJ, Burk-Herrick A, Westerman M, Ayoub NA, Springer MS, Murphy WJ: Impacts of the Cretaceous terrestrial revolution and KPg extinction on mammal diversification. Science. 2011, 334: 521-524. 10.1126/science.1211028.
Jetz W, Thomas GH, Joy JB, Hartmann K, Mooers AO: The global diversity of birds in space and time. Nature. 2012, 491: 444-448. 10.1038/nature11631.
Roelants K, Gower DJ, Wilkinson M, Loader SP, Biju SD, Guillaume K, Moriau L, Bossuyt F: Global patterns of diversification in the history of modern amphibians. Proc Natl Acad Sci U S A. 2007, 104: 887-892. 10.1073/pnas.0608378104.
Grimaldi D, Engel MS: Evolution of the insects. 2005, Cambridge University Press, Cambridge
Stadler T: Mammalian phylogeny reveals recent diversification rate shifts. Proc Natl Acad Sci U S A. 2011, 108: 6187-6192. 10.1073/pnas.1016876108.
Pyron RA, Burbrink FT: Phylogenetic estimates of speciation and extinction rate for testing ecological and evolutionary hypotheses. Trends Ecol Evol. 2013, 28: 729-736. 10.1016/j.tree.2013.09.007.
Moreau CS, Bell CD: Testing the museum versus cradle tropical biological diversity hypothesis: phylogeny, diversification, and ancestral biogeographic range evolution of the ants. Evolution. 2013, 67: 2240-2257. 10.1111/evo.12105.
Rehan SM, Leys R, Schwarz MP: First evidence of a massive extinction event affecting bees close to the K-T boundary. PLoS One. 2013, 8: e76683-10.1371/journal.pone.0076683.
Hunt T, Bergsten J, Levkaničova Z, Papadopoulou A, St John O, Wild R, Hammond PM, Ahrens D, Balke M, Caterino MS, Gómez-Zurita J, Ribera I, Barraclough TG, Bocakova M, Bocak L, Vogler AP: A comprehensive phylogeny of beetles reveals the evolutionary origins of a superradiation. Science. 2007, 318: 1913-1916. 10.1126/science.1146954.
Wang B, Zhang H, Jarzembowski EA: Early Cretaceous angiosperms and beetle evolution. Front Plant Sci. 2013, 4: 360-
Ślipiński SA, Leschen RAB, Lawrence JF: Order coleoptera linnaeus, 1758. Zootaxa. 2011, 3148: 203-208.
Matthews EG, Lawrence JF, Bouchard P, Steiner WE, Slipiński SA: 11.14. Tenebrionidae Latreille, 1802. Handbook of Zoology. A Natural History of the Phyla of the Animal Kingdom. Volume IV - Arthropoda: Insecta. Part 38. Coleoptera, Beetles. Volume 2: Systematics (Part 2). Edited by: Leschen RAB, Beutel RG, Lawrence JF. 2010, Walter de Gruyter, Berlin, 574-659.
Watt JC: A revised subfamily classification of Tenebrionidae (Coleoptera). New Zealand J Zool. 1974, 1: 381-452. 10.1080/03014223.1974.9517846.
Geisthardt M: Tenebrionidae (Insecta, Coleoptera) as an indicator for climatic changes on the Cape Verde Islands. Spec Bull Jap Soc Coleopt. 2003, 6: 331-337.
Kirejtshuk AG, Merkl O, Kernegger F: A new species of the genus Pentaphyllus Dejean, 1821 (Coleoptera: Tenebrionidae, Diaperinae) from the Baltic amber and checklist of the fossil Tenebrionidae. Zoosyst Ross. 2008, 17: 131-137.
Matthews EG: Classification, phylogeny and biogeography of the genera of Adeliini (Coleoptera: Tenebrionidae). Inver Taxon. 1998, 12: 685-824. 10.1071/IT97008.
Matthews EG, Bouchard P: Tenebrionid Beetles of Australia: Description of Tribes, Keys to Genera, Catalogue of Species. 2008, Australian Biological Resources Study, Canberra
McKenna DD, Farrell BD: Beetles (Coleoptera). The Timetree of Life. Edited by: Hedges SB, Kumar S. 2009, Oxford University Press, New York, 278-289.
Crisp MD, Arroyo MTK, Cook LG, Gandolfo MA, Jordan GJ, McGlone MS, Weston PH, Westoby M, Wilf P, Linder HP: Phylogenetic biome conservatism on a global scale. Nature. 2009, 458: 754-756. 10.1038/nature07764.
Chandler MA, Rind D, Ruedy R: Pangaean during the Early Jurassic: GCM simulations and the sedimentary record of paleoclimate. Geol Soc Am Bull. 1992, 104: 543-559. 10.1130/0016-7606(1992)104<0543:PCDTEJ>2.3.CO;2.
Parrish JT: Climate of the supercontinent Pangea. J Geol. 1993, 101: 215-233. 10.1086/648217.
Byrne M, Yeates DK, Joseph L, Kearney M, Bowler J, Williams MAJ, Cooper S, Donnellan SC, Keogh JS, Leys R, Melville J, Murphy DJ, Porch N, Wyrwoll KH: Birth of a biome: insights into the assembly and maintenance of the Australian arid zone biota. Mol Ecol. 2008, 17: 4398-4417. 10.1111/j.1365-294X.2008.03899.x.
Byrne M, Steane DA, Joseph L, Yeates DK, Jordan GJ, Crayn D, Aplin K, Cantrill DJ, Cook LG, Crisp MD, Keogh JS, Melville J, Moritz C, Porch N, Sniderman JMK, Sunnucks P, Weston PH: Decline of a biome: evolution, contraction, fragmentation, extinction and invasion of the Australian mesic zone biota. J Biogeogr. 2011, 38: 1635-1656. 10.1111/j.1365-2699.2011.02535.x.
Matthews EG: Origins of Australian arid-zone tenebrionid beetles. Invert Syst. 2000, 14: 941-951. 10.1071/IT00021.
Zachos JC, Dickens GR, Zeebe RE: An early Cenozoic perspective on greenhouse warming and carbon-cycle dynamics. Nature. 2008, 451: 279-283. 10.1038/nature06588.
Veizer J, Godderis Y, François LM: Evidence for decoupling of atmospheric CO2 and global climate during the Phanerozoic eon. Nature. 2000, 408: 698-701. 10.1038/35047044.
Prokoph A, Shields GA, Veizer J: Compilation and time-series analysis of a marine carbonate δ18O, δ13C,87Sr/86Sr and δ34S database through Earth history.Earth-Sci Rev 2008, 87:113–133.,
Blakey RC: Gondwana Paleogeography From Assembly to Breakup – A 500 Million Year Odyssey. Resolving the Late Paleozoic Ice Age in Time and Space. Edited by: Fielding CR, Frank TD, Isbell JL. 2008, Geological Society of America Special Paper 441. Geological Society of America, Boulder, CO, 1-28. 10.1130/2008.2441(01).
Mitterboeck TF, Adamowicz SJ: Flight loss linked to faster molecular evolution in insects. Proc R Soc B. 2013, 280: 20131128-10.1098/rspb.2013.1128.
Parmesan C: Ecological and evolutionary responses to recent climate change. Annu Rev Ecol Evol Syst. 2006, 37: 637-669. 10.1146/annurev.ecolsys.37.091305.110100.
Ikeda H, Nishikawa M, Sota T: Loss of flight promotes beetle diversification. Nat Commun. 2012, 3: 648-10.1038/ncomms1659.
Cornell HV: Is regional species diversity bounded or unbounded?. Biol Rev. 2013, 88: 140-165. 10.1111/j.1469-185X.2012.00245.x.
Erwin DH: Climate as a driver of evolutionary change. Curr Biol. 2009, 19: R575-R583. 10.1016/j.cub.2009.05.047.
Mayhew PJ, Jenkins GN, Benton TG: A long-term association between global temperature and biodiversity, origination and extinction in the fossil record. Proc R Soc Lond B. 2009, 275: 47-53. 10.1098/rspb.2007.1302.
Condamine FL, Rolland J, Morlon H: Macroevolutionary perspectives to environmental change. Ecol Lett. 2013, 16 (Suppl 1): 72-85. 10.1111/ele.12062.
Kergoat GJ, Soldati L, Clamens A-L, Jourdan H, Jabbour-Zahab R, Genson G, Bouchard P, Condamine FL: Higher-level molecular phylogeny of darkling beetles (Coleoptera, Tenebrionoidea, Tenebrionidae).Syst Entomol. in press.,
Drummond AJ, Suchard MA, Xie D, Rambaut A: Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol. 2012, 29: 1969-1973. 10.1093/molbev/mss075.
Parham JF, Donoghue PC, Bell CJ, Calway TD, Head JJ, Holroyd PA, Inoue JG, Irmis RB, Joyce WG, Ksepka DT, Patané JSL, Smith ND, Tarver JE, van Tuinen M, Yang Z, Angielczyk KD, Greenwood JM, Hipsley CA, Jacobs L, Makovicky PJ, Müller J, Smith KT, Theodor JM, Warnock RCM, Benton MJ: Best practices for justifying fossil calibrations. Syst Biol. 2012, 61: 346-359. 10.1093/sysbio/syr107.
Sauquet H, Ho SYW, Gandolfo MA, Jordan GJ, Wilf P, Cantrill DJ, Bayly MJ, Bromham L, Brown GK, Carpenter RJ, Lee DM, Murphy DJ, Sniderman JMK, Udovicic F: Testing the impact of calibration on molecular divergence times using a fossil-rich group: the case of Nothofagus (Fagales). Syst Biol. 2012, 61: 289-313. 10.1093/sysbio/syr116.
Gradstein FM, Ogg JG, Schmitz M, Ogg G: The Geologic Time Scale 2012. 2012, Elsevier, Boston
Condamine FL, Soldati L, Clamens A-L, Rasplus J-Y, Kergoat GJ: Diversification patterns and processes of wingless endemic insects in the Mediterranean Basin: historical biogeography of the genus Blaps (Coleoptera: Tenebrionidae). J Biogeogr. 2013, 40: 1899-1913. 10.1111/j.1365-2699.2012.02787.x.
Kergoat GJ, Le Ru BP, Genson G, Cruaud C, Couloux A, Delobel A: Phylogenetics, species boundaries and timing of resource tracking in a highly specialized group of seed beetles (Coleoptera: Chrysomelidae: Bruchinae). Mol Phylogenet Evol. 2011, 59: 746-760. 10.1016/j.ympev.2011.03.014.
Fernández-Palacios JM, de Nascimento L, Otto R, Delgado JD, Garcia-del-Rey E, Arévalo JR, Whittaker RJ: A reconstruction of Palaeo-Macaronesia, with particular reference to the long-term biogeography of the Atlantic island laurel forest. J Biogeogr. 2011, 38: 226-246. 10.1111/j.1365-2699.2010.02427.x.
Lopez-Vaamonde C, Wikström N, Labandeira C, Godfray HCJ, Goodman SJ, Cook JM: Fossil-calibrated molecular phylogenies reveal that leaf-mining moths radiated millions of years after their host plants. J Evol Biol. 2006, 19: 1314-1326. 10.1111/j.1420-9101.2005.01070.x.
Biffin E, Hill RS, Lowe AJ: Did Kauri (Agathis: Araucariaceae) really survive the Oligocene drowning of New Zealand?. Syst Biol. 2010, 59: 594-602. 10.1093/sysbio/syq030.
Crisp MD, Cook LG: Cenozoic extinctions account for the low diversity of extant gymnosperms compared with angiosperms. New Phyt. 2011, 192: 997-1009. 10.1111/j.1469-8137.2011.03862.x.
Yang Z, Rannala B: Bayesian estimation of species divergence times under a molecular clock using multiple fossil calibrations with soft bounds. Mol Biol Evol. 2006, 23: 212-226. 10.1093/molbev/msj024.
Ho SYW, Phillips MJ: Accounting for calibration uncertainties in phylogenetic estimation of evolutionary divergence times. Syst Biol. 2009, 58: 367-380. 10.1093/sysbio/syp035.
Lartillot N, Philippe H: Computing Bayes factors using thermodynamic integration. Syst Biol. 2006, 55: 195-207. 10.1080/10635150500433722.
Baele G, Li WLS, Drummond AJ, Suchard MA, Lemey P: Accurate model selection of relaxed clocks in Bayesian phylogenetics. Mol Biol Evol. 2013, 30: 239-243. 10.1093/molbev/mss243.
Rambaut A, Drummond AJ: Tracer v1.5.., [http://beast.bio.ed.ac.uk/Tracer]
Kass RE, Raftery AE: Bayes factors. J Am Stat Assoc. 1995, 90: 773-795. 10.1080/01621459.1995.10476572.
Cusimano N, Renner SS: Slowdowns in diversification rates from real phylogenies may be not real. Syst Biol. 2010, 59: 458-464. 10.1093/sysbio/syq032.
Davies MP, Midford PE, Maddison W: Exploring power and parameter estimation of the BiSSE method for analysing species diversification. BMC Evol Biol. 2013, 13: 38-10.1186/1471-2148-13-38.
Goldberg EE, Kohn JR, Lande R, Robertson KA, Smith SA, Igic B: Species selection maintains self-incompatibility. Science. 2010, 330: 493-495. 10.1126/science.1194513.
Hugall AF, Stuart-Fox D: Accelerated speciation in colour-polymorphic birds. Nature. 2012, 485: 631-634. 10.1038/nature11050.
Magallón S, Sanderson MJ: Absolute diversification rates in angiosperm clades. Evolution. 2001, 55: 1762-1780. 10.1111/j.0014-3820.2001.tb00826.x.
Lewis P: A likelihood approach to estimating phylogeny from discrete morphological character data. Syst Biol. 2001, 50: 913-925. 10.1080/106351501753462876.
Maddison WP, Maddison DR: Mesquite: A Modular System for Evolutionary Analysis, Version 2.75. Available at . 2011., [http://mesquiteproject.org]
Schluter D, Price T, Mooers AO, Ludwig D: Likelihood of ancestor states in adaptive radiation. Evolution. 1997, 51: 1699-1711. 10.2307/2410994.
Ricklefs RE: Estimating diversification rates from phylogenetic information. Trends Ecol Evol. 2007, 22: 602-610. 10.1016/j.tree.2007.06.013.
Rabosky DL: Ecological limits and diversification rate: alternative paradigms to explain the variation in species richness among clades and regions. Ecol Lett. 2009, 12: 735-743. 10.1111/j.1461-0248.2009.01333.x.
Rabosky DL: Primary controls on species richness in higher taxa. Syst Biol. 2010, 59: 634-645. 10.1093/sysbio/syq060.
We would like to thank Alfried Vogler, Emily Crow and four anonymous reviewers for the numerous constructive comments and corrections they made on a preliminary version of the manuscript. Interesting insights on the results of the TreePar diversification analyses were kindly provided by Tanja Stadler. FLC also wishes to thank Hélène Morlon for supervision. Partial funding for this study was obtained through the program “ANR Biodiversit” of the French National Agency for Research (Project BIONEOCAL 2008–2012) and by funds from INRA for GJK. Part of the sequencing was also supported by the program “Bibliothèque du Vivant” (Project CIAM) supported by a joint CNRS, INRA and MNHN consortium. Both FLC and JLA were supported by ANR grants (Projects ECOEVOBIO-CHEX2011 and SPATEVOLEPID-10-PDOC-017-01, respectively). No conflicts of interest were discovered.
The authors declare that they have no competing interests.
Electronic supplementary material
Additional file 1: Figure S1.: Lineages-through-time plots for the four calibration scenarios used in this study. (PDF 40 KB)
Additional file 2: Table S1.: TreePar diversification analyses: means and standard errors based on the analysis of 500 posterior trees. Table S2. TreePar diversification analyses: means and standard errors based on the analysis of 1000 posterior trees. (PDF 47 KB)
Additional file 4: Table S4.: Number of genera per tribe. Sampled tribes are highlighted using bold characters. (PDF 84 KB)
Additional file 5: Additional information on fossil taxa and justifications on ages of geological formations and amber deposits.(PDF 236 KB)
Additional file 6: Table S5.: Comparison of age estimates (under either a BD or a Yule model of speciation) between crown calibration and stem calibration scenarios. For all calibration schemes the median ages and 95% HPD are reported. Bold italic fonts are used to highlight age estimates that are contradicted by known fossil evidences (for the stem calibration scenario). (PDF 70 KB)
Additional file 7: Table S6.: Bayes factor scores (2 ln BF) from Comparisons of alternative calibration procedures of BRC analyses. (PDF 146 KB)
Additional file 8: Figure S2.: Timetree corresponding to the results of the best--fit calibration procedure implemented with BEAST (`Yule crown’). Horizontal bars on nodes correspond to the 95% higher posterior probabilities (HPD) of age estimates. (PDF 277 KB)
Additional file 9: Figure S3.: Timetree corresponding to the results of the `BD crown’ calibration procedure implemented with BEAST.Horizontal bars on Nodes correspond to the 95% higher posterior probabilities (HPD) of age estimates. (PDF 240 KB)
About this article
Cite this article
Kergoat, G.J., Bouchard, P., Clamens, AL. et al. Cretaceous environmental changes led to high extinction rates in a hyperdiverse beetle family. BMC Evol Biol 14, 220 (2014). https://doi.org/10.1186/s12862-014-0220-1