Cretaceous environmental changes led to high extinction rates in a hyperdiverse beetle family
© Kergoat et al.; licensee BioMed Central Ltd. 2014
Received: 8 August 2014
Accepted: 6 October 2014
Published: 21 October 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.
KeywordsCretaceous-palaeogene mass extinction Cretaceous terrestrial revolution Dating analyses Diversification analyses Environmental changes Palaeoclimates Tenebrionidae
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
Age estimates obtained with a crown calibration scenario using either a BD or a Yule model of speciation
95.58 - 159.96
118.88 - 176.46
57.81 - 234.13
71.31 - 186.65
185.73 - 235.21
191.95 - 236.60
152.04 - 177.51
152.00 - 168.20
152.03 - 186.74
152.01 - 177.19
210.27 - 228.34
210.68 - 229.83
120.00 - 142.65
120.02 - 142.83
120.01 - 162.72
120.00 - 158.92
101.99 - 182.42
93.72 - 182.24
108.72 - 136.86
112.40 - 147.30
152.03 - 160.35
152.00 - 157.72
90.40 - 129.19
84.94 - 139.08
95.37 - 146.98
101.22 - 175.36
126.00 - 135.38
126.00 - 135.56
168.40 - 188.71
169.56 - 191.62
122.15 - 175.71
120.95 - 164.16
126.00 - 130.86
126.00 - 130.81
138.98 - 154.64
132.14 - 149.53
59.88 - 97.34
61.63 - 84.12
117.91 - 168.19
116.69 - 171.77
85.25 - 108.91
88.76 - 109.90
80.72 - 105.69
82.69 - 113.95
59.60 - 101.78
52.03 - 107.92
47.94 - 74.39
42.62 - 79.75
53.99 - 79.22
50.34 - 75.20
81.94 - 104.39
87.38 - 107.20
102.87 - 128.11
112.07 - 128.01
77.40 - 117.74
81.84 - 115.80
52.11 - 89.13
55.97 - 94.31
45.09 - 67.74
37.83 - 63.89
49.29 - 67.79
49.07 - 69.47
34.66 - 43.14
35.58 - 46.27
19.48 - 24.00
20.01 - 24.00
Impacts of cretaceous environmental changes on tenebrionid diversification dynamics
Results of the main TreePar analyses for BD and Yule crown calibration scenarios
Trees: BD crown calibration scenario
BD constant 1 shift
BD constant 2 shifts
BD constant 3 shifts
BD constant 4 shifts
shift time #1
shift time #2
shift time #3
shift time #4
Trees: Yule crown calibration scenario
BD constant 1 shift
BD constant 2 shifts
BD constant 3 shifts
BD constant 4 shifts
shift time #1
shift time #2
shift time #3
shift time #4
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 .
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.
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.
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Vakhrameev VA: Jurassic and Cretaceous floras and climates of the Earth. 1991, Cambridge University Press, CambridgeGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- Föllmi KB: Early Cretaceous life, climate and anoxia. Cretaceous Res. 2011, 35: 230-257. 10.1016/j.cretres.2011.12.005.View ArticleGoogle Scholar
- Labandeira CC: A paleobiologic perspective on plant-insect interactions. Curr Opin Plant Biol. 2013, 16: 414-421. 10.1016/j.pbi.2013.06.003.PubMedView ArticleGoogle Scholar
- Benton MJ: The origins of modern biodiversity on land. Phil Trans R Soc B. 2010, 365: 3667-3679. 10.1098/rstb.2010.0269.PubMedPubMed CentralView ArticleGoogle Scholar
- Farrell BD: “Inordinate fondness” explained: why are there so many beetles?. Science. 1998, 281: 555-559. 10.1126/science.281.5376.555.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- Grimaldi D, Engel MS: Evolution of the insects. 2005, Cambridge University Press, CambridgeGoogle Scholar
- Stadler T: Mammalian phylogeny reveals recent diversification rate shifts. Proc Natl Acad Sci U S A. 2011, 108: 6187-6192. 10.1073/pnas.1016876108.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Wang B, Zhang H, Jarzembowski EA: Early Cretaceous angiosperms and beetle evolution. Front Plant Sci. 2013, 4: 360-PubMedPubMed CentralGoogle Scholar
- Ślipiński SA, Leschen RAB, Lawrence JF: Order coleoptera linnaeus, 1758. Zootaxa. 2011, 3148: 203-208.Google Scholar
- 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.Google Scholar
- Watt JC: A revised subfamily classification of Tenebrionidae (Coleoptera). New Zealand J Zool. 1974, 1: 381-452. 10.1080/03014223.1974.9517846.View ArticleGoogle Scholar
- 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.Google Scholar
- 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.Google Scholar
- Matthews EG: Classification, phylogeny and biogeography of the genera of Adeliini (Coleoptera: Tenebrionidae). Inver Taxon. 1998, 12: 685-824. 10.1071/IT97008.View ArticleGoogle Scholar
- Matthews EG, Bouchard P: Tenebrionid Beetles of Australia: Description of Tribes, Keys to Genera, Catalogue of Species. 2008, Australian Biological Resources Study, CanberraGoogle Scholar
- McKenna DD, Farrell BD: Beetles (Coleoptera). The Timetree of Life. Edited by: Hedges SB, Kumar S. 2009, Oxford University Press, New York, 278-289.Google Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- Parrish JT: Climate of the supercontinent Pangea. J Geol. 1993, 101: 215-233. 10.1086/648217.View ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- Matthews EG: Origins of Australian arid-zone tenebrionid beetles. Invert Syst. 2000, 14: 941-951. 10.1071/IT00021.View ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.,Google Scholar
- 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).View ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- Ikeda H, Nishikawa M, Sota T: Loss of flight promotes beetle diversification. Nat Commun. 2012, 3: 648-10.1038/ncomms1659.PubMedPubMed CentralView ArticleGoogle Scholar
- Cornell HV: Is regional species diversity bounded or unbounded?. Biol Rev. 2013, 88: 140-165. 10.1111/j.1469-185X.2012.00245.x.PubMedView ArticleGoogle Scholar
- Erwin DH: Climate as a driver of evolutionary change. Curr Biol. 2009, 19: R575-R583. 10.1016/j.cub.2009.05.047.PubMedView ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- Condamine FL, Rolland J, Morlon H: Macroevolutionary perspectives to environmental change. Ecol Lett. 2013, 16 (Suppl 1): 72-85. 10.1111/ele.12062.PubMedView ArticleGoogle Scholar
- 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.,Google Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Gradstein FM, Ogg JG, Schmitz M, Ogg G: The Geologic Time Scale 2012. 2012, Elsevier, BostonGoogle Scholar
- 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.Google Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Lartillot N, Philippe H: Computing Bayes factors using thermodynamic integration. Syst Biol. 2006, 55: 195-207. 10.1080/10635150500433722.PubMedView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Hugall AF, Stuart-Fox D: Accelerated speciation in colour-polymorphic birds. Nature. 2012, 485: 631-634. 10.1038/nature11050.PubMedView ArticleGoogle Scholar
- Magallón S, Sanderson MJ: Absolute diversification rates in angiosperm clades. Evolution. 2001, 55: 1762-1780. 10.1111/j.0014-3820.2001.tb00826.x.PubMedView ArticleGoogle Scholar
- Lewis P: A likelihood approach to estimating phylogeny from discrete morphological character data. Syst Biol. 2001, 50: 913-925. 10.1080/106351501753462876.PubMedView ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- Ricklefs RE: Estimating diversification rates from phylogenetic information. Trends Ecol Evol. 2007, 22: 602-610. 10.1016/j.tree.2007.06.013.View ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Rabosky DL: Primary controls on species richness in higher taxa. Syst Biol. 2010, 59: 634-645. 10.1093/sysbio/syq060.PubMedView ArticleGoogle Scholar
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 (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.