Early evolution of beetles regulated by the end-Permian deforestation
Abstract
The end-Permian mass extinction (EPME) led to a severe terrestrial ecosystem collapse. However, the ecological response of insects—the most diverse group of organisms on Earth—to the EPME remains poorly understood. Here, we analyse beetle evolutionary history based on taxonomic diversity, morphological disparity, phylogeny, and ecological shifts from the Early Permian to Middle Triassic, using a comprehensive new dataset. Permian beetles were dominated by xylophagous stem groups with high diversity and disparity, which probably played an underappreciated role in the Permian carbon cycle. Our suite of analyses shows that Permian xylophagous beetles suffered a severe extinction during the EPME largely due to the collapse of forest ecosystems, resulting in an Early Triassic gap of xylophagous beetles. New xylophagous beetles appeared widely in the early Middle Triassic, which is consistent with the restoration of forest ecosystems. Our results highlight the ecological significance of insects in deep-time terrestrial ecosystems.
Editor's evaluation
The study proposes a new evolutionary-ecological scenario for Late Paleozoic and early Mesozoic beetles, supported by the summary of all available knowledge about early beetle fossils, including analyses of their taxon and morphological diversity and phylogenetic relationships. The effects of xylophagous beetles during the Paleozoic may have played a fundamental role in global biochemical cycles. The results advance our understanding of the evolutionary success of beetles and the many ways in which large environmental changes may affect biodiversity in general.
https://doi.org/10.7554/eLife.72692.sa0Introduction
The end-Permian mass extinction (EPME; approximately 252 million years ago) was the most severe extinction event in the Phanerozoic (Benton and Newell, 2014). The EPME was primarily caused by the eruption of the Siberian flood basalts (Burgess and Bowring, 2015; Fielding et al., 2019), which generated excessive emissions of thermogenic methane, CO2, and SO2 that cascaded rapid global warming (Wu et al., 2021; Black et al., 2018), oceanic acidification and anoxia/euxinia (Schobben et al., 2020), aridification and other shifts in the hydrological cycle (Sun et al., 2012), acid rain (Black et al., 2014), wildfires (Shen et al., 2011), and ozone destruction (Benca et al., 2018). The response of terrestrial ecosystems to the EPME is quite heterogeneous, probably due to biotic, topographic, and latitudinal differences (Fielding et al., 2019; Zhao et al., 2020; Dal Corso et al., 2020). Moreover, how terrestrial ecosystems were affected during the EPME is still highly controversial (Benton and Newell, 2014; Gastaldo, 2019). Terrestrial tetrapods and plants are considered to have been severely affected by the EPME mostly based on diversity and taxonomic composition (Benton and Newell, 2014; Viglietti et al., 2021); however, such mass extinction was questioned by a more comprehensive dataset of plant macro- and microfossils (Gastaldo, 2019; Nowak et al., 2019). Similarly, Permian insects are thought to have suffered a significant extinction (Labandeira and Sepkoski, 1993; Béthoux et al., 2005; Labandeira, 2005; Condamine et al., 2020; Condamine et al., 2016), but this was not supported by other molecular phylogenetic and fossil record analyses (Ponomarenko, 2016; Dmitriev et al., 2018; Montagna et al., 2019; Schachat et al., 2019). In addition, the ecological response of insects to the EPME remains poorly understood (Benton and Newell, 2014; Schachat et al., 2021).
Beetles (Coleoptera) are the most speciose group of extant insects (Stork, 2018), with a stratigraphic range dating back to at least the lowest Permian (Ponomarenko, 2016; Kirejtshuk et al., 2013). They have a rich fossil record since the Permian and display a wide array of lifestyles (Figure 1; Ponomarenko, 1969; Ponomarenko, 2003). Their fossil record thus offers a unique and complementary perspective for studying the ecological response of insects to the EPME. The evolutionary history of Coleoptera has been widely investigated through molecular phylogenetic analyses (Condamine et al., 2016; McKenna et al., 2019; Zhang et al., 2018), morphological phylogenetic analyses (Beutel et al., 2008; Beutel et al., 2019), and fossil record analyses (Ponomarenko, 2003; Ponomarenko, 2016; Smith and Marcot, 2015). Although a long-term Palaeozoic-Mesozoic turnover of beetle assemblages is supported by almost all analyses, the detailed ecological response to the EPME and its explanatory mechanisms remain unclear. Most of the Permian and Triassic beetles belong to stem groups (extinct suborders or families; Figure 1), and thus they show character combinations and evolutionary history that cannot be inferred or predicted from phylogenetic analysis of modern beetles. In particular, two problems were ignored by previous analyses. First, phylogenetic relationships of some key fossils remain poorly resolved, particularly in their evolutionary relationships to modern taxa. Second, there are two complementary taxonomic systems for Permian and Triassic beetles: one is artificial formal taxa (based on isolated elytra that cannot be definitely classified into any natural group), and the other is the natural taxonomy (commonly based on more complete fossils including bodies and elytra). The formal taxa, like trace fossil taxa, lack comprehensive phylogenetic data (Ponomarenko, 2004), and thus they cannot be used unreservedly for biodiversity and phylogenetic analyses, but can be helpful in the morphospace analysis. These issues cloud the temporal resolution of coleopteran biodiversity in deep time and complicate the evolutionary trajectory of beetles but can be overcome through a combination of multiple analytical methods. Therefore, taxonomic diversity, morphological disparity, and ecological shifts are best evaluated jointly to better understand how the EPME has shaped the evolutionary history of beetles.
Here, we compile an updated database of beetles from the Early Permian to Middle Triassic based on the taxonomic revision of fossils (including formal taxa). We analyse the evolution of taxonomic diversity, morphological disparity, and palaeoecological shifts of beetles from the Early Permian to Middle Triassic through phylogenetic and palaeoecological reconstructions and morphospace analyses of fossil material. Our results suggest that xylophagous (feeding on or in wood) beetles probably played a key and underappreciated role in the Permian carbon cycle and that the EPME had a profound ecological influence on beetle evolution. These results provide new insights into the ecological role of insects in deep-time terrestrial ecosystems and the ecological response of insects to deforestation and global warming.
Results
Taxonomic diversity
We compiled an updated database of beetles (21 families, 125 genera, and 299 species) from the Early Permian to Middle Triassic based on the taxonomic revision of natural and formal taxa (Figure 2—source data 1). Our database contains 18 families, 109 genera, and 220 species of natural taxa. There is a steady increase of families from the Early Permian to Middle Triassic, which is consistent with the result of Smith and Marcot, 2015, whose analyses were only conducted at the family level. The diversity of natural taxa displays almost the same trajectory at both species and genus levels (Figure 2C and E). The diversity is roughly stable in the Early Permian (Cisuralian), mainly represented by Tshekardocoleidae (Figure 1), increases rapidly in the Middle Permian (Guadalupian) and Late Permian (Lopingian), with the rise of the major clades Permocupedidae (Permocupedinae and Taldycupedinae) and Rhombocoleidae. Subsequently, it plunges in the Early Triassic and recovered gradually from the Anisian (early Middle Triassic). In the Ladinian (late Middle Triassic), the diversity clearly exceeds that of the Late Permian (Figure 2C and E). From the Middle Triassic, the Permian coleopteran assemblage characterized by Tshekardocoleidae, Permocupedidae, and Rhombocoleidae is completely replaced by a Triassic assemblage dominated by Cupedidae, Phoroschizidae and Triaplidae.
Our database also contains 3 families, 17 genera, and 79 species of formal taxa. A considerable proportion of Permian beetles belong to such taxa (Permosynidae, Schizocoleidae, and Asiocoleidae). These formal taxa mostly belong to stem groups, but some should probably be attributed to the two extant suborders Adephaga and Polyphaga. Both species and genus-level diversities of formal taxa show a gradual increase from the Middle to Late Permian, but decrease distinctly from the Triassic (Figure 2—figure supplement 1). The mixed taxa diversity (combining natural and formal taxa) displays the same trajectory to that of natural taxa at both species and genus levels (Figure 2D and F).
Phylogeny
We carried out a phylogenetic analysis based on 93 adult and larval characters across 15 natural taxa representing all natural families of Coleoptera from the Early Permian to Middle Triassic (Figure 3—source data 1). Our parsimony analysis result is consistent with a previous analysis (Beutel et al., 2008), and confirms that Tshekardocoleidae, Permocupedidae (Permocupedinae and Taldycupedinae), and Rhombocolediae are the stem group of Coleoptera (Figure 3A, Figure 3—figure supplement 1).
Morphological disparity
We chose beetle elytra—hardened forewings primarily serving as protective covers for the hindwings and body underneath—to perform the morphological disparity analysis for three reasons: (1) elytra are the most commonly preserved fossils of Palaeozoic and Mesozoic beetles, and they are easily accessible in the literature and in online databases; (2) Permian and Triassic elytra display complex morphological structure (Ponomarenko, 1969); (3) elytra morphology has long been studied in relation to taxonomic diversity of living and extinct beetles (Ponomarenko, 2004; Tong et al., 2021).
We assembled two discrete character matrices (at species and genus levels) based on 35 characters of 197 genera and 346 species (including undetermined species and unnamed specimens) for morphological disparity analyses (Figure 4—source data 1). The taxa were ordinated into a multivariate morphospace using both principal coordinates analysis (PcoA) and non-metric multidimensional scaling (NMDS) with two distance metrics, including the generalized Euclidean distance (GED) and maximum observable rescale distance (MORD). We chose both sum of variance (sov) and product of variance (pov) as the proxy for morphological disparity due to their robustness in sample size (Simões et al., 2020). The use of discrete characters produces results that have non-metric properties, but this approach can be used to elucidate broad patterns of similarities and clustering within multidimensional space (Lloyd, 2016; Deline et al., 2018).
The patterns of morphospace occupation of beetles in different time-bins are shown in three-dimensional plots delimited by combinations of the first three axes of the PcoA and NMDS results based on MORD metrics (Figure 4A, Figure 4—figure supplement 1). The morphological disparity results of two ordination methods within the MORD and GED matrix shows the same trajectory at both genus and species levels. The evolutionary pattern of morphological disparity is robust in different disparity metrics. Disparity is low in the Early Permian, with a significant increase during the Middle Permian. It is roughly stable in the Middle and Late Permian, subsequently showing a distinct plunge in the Early Triassic, slightly recovering in the Middle Triassic but is still significantly lower than in the Middle and Late Permian (Figure 4, Figure 4—figure supplements 1–8).
Discussion
Our results demonstrate that beetles display a steady accumulation of taxonomic diversity throughout the Permian (Figure 2). The earliest definite beetles are Tshekardocoleidae (including the genus Coleopsis) from the Early Permian of Germany, Czech Republic, USA, and Russia (Figure 1), although the origin of Coleoptera is dated to the Carboniferous by molecular phylogenetic analysis (McKenna et al., 2019; Zhang et al., 2018; Toussaint et al., 2016). The coleopteran diversity radiation in the Middle Permian is consistent with an expansion of morphological disparity, corresponding to the appearance of multiveined, smooth, and striate elytra as well as some other patterns (Figure 4, Figure 4—figure supplement 1). Taxonomic diversity and morphological disparity were decoupled during the Late Permian when taxonomic diversity increased but morphological disparity was almost stable (Figures 2 and 4). The abrupt Middle Permian increase of coleopteran morphological disparity conforms to the early burst model of clade disparity commonly arising early in radiations (Simões et al., 2020; Hughes et al., 2013).
The Permian coleopteran assemblage was dominated by stem groups including Tshekardocoleidae, Permocupedidae (Permocupedinae and Taldycupedinae), and Rhombocolediae in terms of richness and abundance. These ancient beetles were most likely xylophagous because they display a prognathous head, a characteristic elytral pattern with window punctures, a cuticular surface with tubercles (or scales), and a plesiomorphic pattern of ventral sclerites, very similar to the extant wood associated archostematans (Figure 1; Kirejtshuk et al., 2013; Ponomarenko, 1969). Moreover, Permian trace fossils showing wood boring provide convincing evidence for the xylophagous habit of these ancient beetles (Naugolnykh and Ponomarenko, 2010; Feng et al., 2019). Aquatic or semi-aquatic beetles including Phoroschizidae and Ademosynidae, belonging to the suborder Archostemata, first appeared in the Middle Permian and diversified in the Late Permian (Ponomarenko, 2003). The three other suborders of Coleoptera, comprising Polyphaga, Adephaga, and Myxophaga, most likely evolved by the Late Permian, but definite fossils are rare at this time.
Permian beetles probably played an important ecological role in forest ecosystems because most Permian beetles were most likely xylophagous insects that consumed living and dead woody stems (Figure 3). Some Permian xylophagous beetles fed on living wood tissues (Feng et al., 2017; Feng et al., 2019), which likely reduced tree productivity and could have caused extensive tree mortality. Insect-mediated tree mortality is known to result in large transfers of carbon from biomass to dead organic matter (Seidl et al., 2018; Fei et al., 2019). The other Permian xylophagous beetles were likely saproxylic (feeding on dead wood) (Ponomarenko, 2003), and they could also impact terrestrial carbon dynamics by accelerating wood decomposition (Ulyshen, 2018). Saproxylic animals first appeared in the Devonian and are mainly represented by small invertebrates such as oribatid mites, until the Permian (Labandeira et al., 1997; Labandeira, 2007). Whereas grazing by micro- and meso-invertebrates (nematodes, collembolans, enchytraeids and oribatid mites) did not significantly affect wood decomposition, consumption by macro-invertebrates (dominated by saproxylic beetles and termites in modern ecosystems) significantly sped up wood decomposition (Tapanila and Roberts, 2012). In addition to those that directly facilitated decomposition by consuming wood, Permian saproxylic beetles are likely to have had a variety of indirect effects on decomposition, including creating tunnels that facilitate the movement of fungi into wood (Naugolnykh and Ponomarenko, 2010; Feng et al., 2017), and vectoring fungi and other decay organisms on or within their bodies, like their extant counterparts (Ulyshen, 2016). In conclusion, Permian beetles that feed on living and dead wood probably could impact terrestrial carbon dynamics by reducing forests’ carbon sequestration capacity, and by converting live materials to dead organic matter and subsequent decomposition.
The oxygen concentration of the atmosphere began to rise in the early Palaeozoic, probably with a peak in the Carboniferous and large decline from the beginning of the Permian (Dahl et al., 2010; Berner, 2009; Krause et al., 2018). The reason for this plunge was attributed to a tectonic- or climate-driven reduction in the extent of coal swamps (Berner and Canfield, 1989) or to the evolution of lignin-consuming fungi (Floudas et al., 2012). However, global recoverable coal is only equivalent to a few percent of the oxygen budget in the atmosphere, and thus cannot account for the large drop of atmospheric oxygen (Nelsen et al., 2016). Furthermore, lignin-consuming fungi may have been present before the Carboniferous (Nelsen et al., 2016). Recently, a new geochemical model proposed that the development of Permian terrestrial herbivores may have limited transport and long-term burial of terrestrial organic compounds in marine sediments, resulting in less organic carbon burial and attendant declines in atmospheric oxygen (Laakso et al., 2020). Herbivorous insects and amniotes are thought to be the major herbivorous animals during the Permian and thus are considered to be the most important drivers of the Permian change in biogeochemical cycles of carbon (Laakso et al., 2020). However, Permian herbivorous amniotes mainly fed on leaves, stems, roots, and rhizomes (Sues and Reisz, 1998; Pearson et al., 2013) and could normally digest cellulose by fermentation but could not consume lignin, as in extant herbivorous vertebrates (Pearson et al., 2013). The majority of terrestrial plant biomass is stored in forest woody tissue consisting of decay-resistant lignin (Hibbett et al., 2016; Bar-On et al., 2018). In modern forests, the total carbon stock in woody tissue (including living and dead wood) is approximately 340 Pg carbon, much more than 72 Pg carbon in roots (below ground), 43 Pg carbon in foliage, and 43 Pg carbon in litter (Reich et al., 2014; Pan et al., 2011). In extant forest ecosystems, insects may account for 29% of the total carbon flux from dead wood and thus they have a functional importance in the decomposition of dead wood and the carbon cycle (Seibold et al., 2021). During the Permian, beetles were probably the dominant consumers of woody tissue, while a few other insect groups may have sometimes fed on dead wood (such as stem dictyopterans and protelytropterans) (Grimaldi and Engel, 2005). Permian beetles had probably evolved close interactions with various microorganisms, especially lignin-consuming fungi (Nelsen et al., 2016), which also accelerated the decomposition of dead wood. The Early Permian onset of the decrease in oxygen concentrations is consistent with the origin and radiation of the xylophagous beetles in the fossil record. Therefore, we propose that Permian xylophagous beetles could have been responsible for at least part of the change in Permian biogeochemical cycles in Laakso’s model (Laakso et al., 2020).
As the most taxonomically and functionally diverse group of living organisms on Earth (Stork, 2018), extant insects have significant effects on terrestrial carbon and nutrient cycling by modulating the quality and quantity of resources that enter the detrital food web (Belovsky and Slade, 2000; Kurz et al., 2008; Yang and Gratton, 2014; Seibold et al., 2021). However, the effects of insects on terrestrial ecosystems in deep time have been viewed as unimportant or overlooked (Doughty, 2017). Permian beetles were among the principal degraders of wood and played a fundamental role in deep-time carbon and nutrient cycling and niche creation. Insects may have been one of the major regulating factors of forest ecosystems at least from the Permian.
Our results show that both the taxonomic diversity and morphological disparity dropped dramatically during the Early Triassic (Figures 2 and 4). Combined with the phylogenetic results (Figure 3), our suite of analyses yields a clear ecological signal from beetles across the Permian/Triassic boundary: all xylophagous stem-group beetles become extinct near the Permian-Triassic boundary or abruptly decreased in the Early Triassic (a pattern called ‘dead clade walking’; Barnes et al., 2021), while aquatic phoroschizid and ademosynid lineages crossed the Permian/Triassic boundary and diversified in the Middle Triassic. Coleoptera recovered in taxonomic diversity during the Middle Triassic by the rise of new predatory and herbivorous groups, synchronized with the recovery of terrestrial ecosystems (Zhao et al., 2020). However, the morphological disparity is significantly lower than that of the Middle and Late Permian due to the lack of stem-group beetles that possess complex elytra structures (Figure 4, Figure 4—figure supplement 1). Polyphagan groups increased in taxonomic diversity during the Middle Triassic, which is a transitory epoch from a Palaeozoic stem-group beetle assemblage to a Mesozoic polyphagan-dominated assemblage.
Xylophagous groups are absent or rare in Early Triassic coleopteran assemblages, becoming widespread again from the Middle Triassic, mainly represented by more derived archostematans (such as Cupedidae) and polyphagans (Ponomarenko, 2003). This gap in xylophagous beetles coincided chronologically with the gap in coal deposition (‘coal gap’), a time during which peat-forming forests were rare or absent (Figure 3B), extending across at least the entire Early Triassic (Benton and Newell, 2014; Retallack et al., 1996; Nowak et al., 2020; Zhao et al., 2020). During the latest Permian and earliest Triassic, gymnosperm-dominated forests abruptly collapsed (Vajda et al., 2020) and were replaced by other biomes (such as isoetalean-dominated herbaceous heathlands; Feng et al., 2020) in most areas due to extreme conditions including aridity (Sun et al., 2012), wildfires (Shen et al., 2011), and ozone destruction (Benca et al., 2018). In some regions, the plant extinction was less severe, or the recovery was rapid (Hochuli et al., 2010), or there may have been multiple crises during the Early Triassic (Schneebeli-Hermann et al., 2017), but even short-term ecosystem disruption could have led to extinctions among xylophagous beetles. Previous studies have not provided a clear picture of insect evolution in response to possible environmental stresses, nor of their response to the EPME (Benton and Newell, 2014). Our results show for the first time that the demise of most forests (deforestation event; Vajda et al., 2020) most likely resulted in the extinction of most Palaeozoic xylophagous beetles, analogous to the extinction of tree-dwelling birds and mammals resulting from end-Cretaceous deforestation (Field et al., 2018; Hughes et al., 2021).
Our results reveal an Early Triassic gap in xylophagous beetles, suggesting that early archaic beetles experienced the severe ecological consequences of end-Permian deforestation. Extant insects are suffering from dramatic declines in abundance and diversity largely due to the anthropogenic deforestation and global warming (van Klink et al., 2020; Wagner et al., 2021). However, xylophagous insects have been largely neglected in studies of the current extinction crisis (van Klink et al., 2020). In particular, the diversity and abundance of xylophagous beetles are extremely sensitive to climate change and can also entail forest collapse and carbon cycle disturbance (Kurz et al., 2008; Fei et al., 2019; Šamonil et al., 2020). Our findings may help to better understand future changes in insect diversity and abundance and its consequences faced with global environmental change.
Materials and methods
Diversity analysis
Request a detailed protocolWe compiled an updated database of all coleopteran species from the Early Permian to Middle Triassic from published literature in the early 1800s through to early 2020. In addition, we incorporated data from other open access database projects, including the Fossil Insect Database (EDNA) and Paleobiology Database (PBDB). We re-examined all published occurrences and taxonomy of Coleoptera from the Early Permian to Middle Triassic (Figure 2—source data 1). We standardized and corrected for nomenclatural consistency of all taxa using a classification of extinct beetle taxa above the genus rank (Bouchard et al., 2011). The data were filtered and cleaned by removing or reassigning illegitimate, questionable, and synonymous taxa and converting local to global chronostratigraphic units (Supplementary file 1).
We allocated fossil species into 12 stage-level time-bins covering the Early Permian-Middle Triassic interval (from the Asselian to Ladinian, 298–237 Ma). Considering the short duration of the Induan and Olenekian stages, we combined both stages into one time-bin. The formal taxa were erected based only on isolated elytra that cannot be classified definitely into any natural group. Thus, we separately counted the diversity of natural, formal, and mixed groups (Figure 2, Figure 2—figure supplement 1). The species Coleopsis archaica was attributed to Tschekardocoleidae by Kirejtshuk et al., 2013, but was later elevated to a new family Coleopseidae by Kirejtshuk, 2020. We followed the former opinion because it is premature to erect a family without a detailed cladistic analysis. We determined the stratigraphical ranges of families, genera, and species as the maximum and minimum ages in stage-level time-bins. All diversity was calculated using the range-through method.
Phylogenetic analysis
Request a detailed protocolIn light of the new taxa and characters available for further testing the phylogenetic status of ancient stem-group beetles, we reconstructed the phylogenetic relationships among the stem groups by incorporating the presently described new taxa and revised characters coding into the previous dataset (Beutel et al., 2008). The morphological characters used for phylogenetic analysis comprise 93 adult and larval characters (Figure 3—source data 1). Unknown characters were coded as ‘?’. The taxon sampling contains two megalopterans as outgroups (Sialis and Chauliodes) and 13 coleopteran ingroup taxa (five extant and eight extinct) representing all four coleopteran extant suborders and their stem groups (Supplementary file 2). Compared to previous character matrices, we added the subfamily Taldycupedinae and three new characters (Supplementary file 3). The matrix was analysed in TNT version 1.1, through parsimony analysis and using traditional search (Goloboff et al., 2008). All characters were equally weighed and unordered (1000 replicates and 1000 trees saved per replication). Bootstrap values, consistency index, and retention index were provided (Figure 3—figure supplement 1).
Morphospace analysis
Request a detailed protocolWe performed morphospace analyses with our newly assembled discrete character matrices (Figure 4—source data 1; Supplementary file 4). The analyses were performed using the free software R. 4.0.4. Both the MORD matrix and GED matrix were calculated based on two discrete character matrices (Lloyd, 2016; Wills, 1998). Recent research has revealed that the GED matrix creates a systematic bias in cases with a high percentage of missing data, and the MORD matrix can provide greater fidelity under these circumstances (Lloyd, 2016; Lehmann et al., 2019). We then ordinated all taxa into a multivariate morphospace with both PcoA and NMDS.
For PcoA, we used the function ‘ordinate cladistic matrix’ in package ‘Claddis’ with a cailliez method to correct the negative eigenvalues (Lloyd, 2016). Two disparity matrices were used to evaluate the volume of the morphospace, including the sum and product of the variances (Wills et al., 2016). The product metrics was normalized by taking the nth root (n equals the number of axes used for calculating disparity metrics). We used the scores on all axes that together comprise 90% of total variance to calculate those disparity metrics. We chose a permutation test (two-tailed) to test the null hypothesis of no difference between insect disparity of different time-bins. Each test run used 5000 replications. The test statistic was obtained by using the disparity metric of an older time-bin to minus that of a younger time-bin. If the proportion in the null distribution greater than the observed value of the test statistic is smaller than 0.025, the insect disparity of an older time-bin was considered significantly larger than that of a younger time-bin, and if the proportion was greater than 0.975, the insect disparity of a younger time-bin was considered significantly larger than that of an older time-bin (Figure 4—source data 2).
We also performed permutation tests with sample size corrected. For two design groups with different sample sizes, we first performed subsampling of the group with more samples to obtain equal sample sizes. Based on the newly obtained two groups with equal sample sizes, we calculated the observed value of the test statistic. Then we randomly permutated those species into different groups once and calculated a test statistic. We repeated this procedure (subsampling and permutation) 10,000 times and obtained a null distribution plus a group of observed values. We then calculated a set of proportions greater than the observed values in the null distribution. By analogy, if the median of the proportions is equal to or smaller than 0.025, the insect disparity of an older time-bin is significantly larger than that of a younger time-bin. If the median proportion is considered significantly greater than 0.975, the insect disparity of a younger time-bin is larger than that of an older time-bin (Figure 4—figure supplements 4–7).
For NMDS, we used the function ‘metaMDS’ in package ‘vegan’ with the number of dimension settings to 3 (Dixon, 2003). Both non-metric fit and linear fit were very high (larger than 0.90; Figure 4—figure supplements 4–7) and the stresses were smaller than 0.2, which implies that the ordinations are relatively good. Then we repeated all the previous analyses with this NMDS morphospace and acquired a similar result. The two disparity metrics of each time-bin were calculated based on two different distance matrices and two different ordination methods (Figure 4, Figure 4—figure supplements 1–5). The distribution was simulated under 500 bootstraps. Thirty-one undetermined specimens from the Grès à Voltzia Formation (Lower/Middle Triassic boundary, France) were included in our database and their age was attributed to the early Middle Triassic in our analysis (Figure 4, Figure 4—figure supplements 1–5). Considering that the age of these specimens is controversial, we repeated all our analyses assuming that the age of these specimens is Early Triassic; the result is consistent with the previous one (Figure 4—figure supplements 6 and 7).
Data availability
All source data are available at https://doi.org/10.5061/dryad.7m0cfxpvd. In addition, the source data files (Supplementary Data 1-4) have been provided for figures 2-4 and appendix figures 1-10.
-
Dryad Digital RepositoryDataset of Early evolution of beetles regulated by the end-Permian deforestation.https://doi.org/10.5061/dryad.7m0cfxpvd
References
-
Impacts of global warming on Permo-Triassic terrestrial ecosystemsGondwana Research 25:1308–1337.https://doi.org/10.1016/j.gr.2012.12.010
-
A new model for atmospheric oxygen over Phanerozoic timeAmerican Journal of Science 289:333–361.https://doi.org/10.2475/ajs.289.4.333
-
Phanerozoic atmospheric oxygen: new results using the GEOCARBSULF modelAmerican Journal of Science 309:603–606.https://doi.org/10.2475/07.2009.03
-
The Triassic radiation of the entomofaunaComptes Rendus Palevol 4:609–621.https://doi.org/10.1016/j.crpv.2005.06.005
-
Sulfur isotopic evidence for sources of volatiles in Siberian Traps magmasEarth and Planetary Science Letters 394:58–69.https://doi.org/10.1016/j.epsl.2014.02.057
-
VEGAN, a package of R functions for community ecologyJournal of Vegetation Science 14:927–930.https://doi.org/10.1111/j.1654-1103.2003.tb02228.x
-
Insect diversity from the Carboniferous to RecentPaleontological Journal 52:610–619.https://doi.org/10.1134/S0031030118060047
-
Herbivores increase the global availability of nutrients over millions of yearsNature Ecology & Evolution 1:1820–1827.https://doi.org/10.1038/s41559-017-0341-1
-
Climate, decay, and the death of the coal forestsCurrent Biology 26:R563–R567.https://doi.org/10.1016/j.cub.2016.01.014
-
Rapid demise and recovery of plant ecosystems across the end-Permian extinction eventGlobal and Planetary Change 74:144–155.https://doi.org/10.1016/j.gloplacha.2010.10.004
-
Evolution of the elytral venation and structural adaptations in the oldest Palaeozoic beetles (Insecta: Coleoptera: Tshekardocoleidae)Journal of Systematic Palaeontology 12:575–600.https://doi.org/10.1080/14772019.2013.821530
-
Stepwise oxygenation of the Paleozoic atmosphereNature Communications 9:4081.https://doi.org/10.1038/s41467-018-06383-y
-
The fossil record of insect extinction: new approaches and future directionsAmerican Entomologist 51:14–29.https://doi.org/10.1093/ae/51.1.14
-
Estimating morphological diversity and tempo with discrete character-taxon matrices: implementation, challenges, progress, and future directionsBiological Journal of the Linnean Society 118:131–151.https://doi.org/10.1111/bij.12746
-
No mass extinction for land plants at the Permian-Triassic transitionNature Communications 10:384.https://doi.org/10.1038/s41467-018-07945-w
-
Palaeophytogeographical Patterns Across the Permian–Triassic BoundaryFrontiers in Earth Science 8:613350.https://doi.org/10.3389/feart.2020.613350
-
Reconstructing the diversity of early terrestrial herbivorous tetrapodsPalaeogeography, Palaeoclimatology, Palaeoecology 372:42–49.https://doi.org/10.1016/j.palaeo.2012.11.008
-
Historical development of archostematan beetlesTrudy Paleontologicheskogo Instituta AN SSSR 125:1–240.
-
Ecological evolution of beetles (Insecta: Coleoptera)Acta Zoologica Cracoviensia 46:319–328.
-
Beetles (Insecta, Coleoptera) of the late Permian and early TriassicPaleontological Journal 38:S185–S196.
-
Insects during the time around the Permian—Triassic crisisPaleontological Journal 50:174–186.https://doi.org/10.1134/S0031030116020052
-
Global coal gap between Permian–Triassic extinction and Middle Triassic recovery of peat-forming plantsGeological Society of America Bulletin 108:195–207.https://doi.org/10.1130/0016-7606(1996)108<0195:GCGBPT>2.3.CO;2
-
Indirect biogeomorphic and soil evolutionary effects of spruce bark beetleGlobal and Planetary Change 195:103317.https://doi.org/10.1016/j.gloplacha.2020.103317
-
A Cretaceous peak in family-level insect diversity estimated with mark-recapture methodologyProceedings. Biological Sciences 286:2054.https://doi.org/10.1098/rspb.2019.2054
-
Are insects heading toward their first mass extinction? Distinguishing turnover from crises in their fossil recordAnnals of the Entomological Society of America 114:99–118.https://doi.org/10.1093/aesa/saaa042
-
Palynofloral associations before and after the Permian–Triassic mass extinction, Kap Stosch, East GreenlandGlobal and Planetary Change 155:178–195.https://doi.org/10.1016/j.gloplacha.2017.06.009
-
Invasive alien pests threaten the carbon stored in Europe’s forestsNature Communications 9:1626.https://doi.org/10.1038/s41467-018-04096-w
-
Calibrating the end-Permian mass extinctionScience 334:1367–1372.https://doi.org/10.1126/science.1213454
-
The fossil record and macroevolutionary history of the beetlesProceedings. Biological Sciences 282:20150060.https://doi.org/10.1098/rspb.2015.0060
-
How Many Species of Insects and Other Terrestrial Arthropods Are There on Earth?Annual Review of Entomology 63:31–45.https://doi.org/10.1146/annurev-ento-020117-043348
-
Origins and early evolution of herbivory in tetrapodsTrends in Ecology & Evolution 13:141–145.https://doi.org/10.1016/S0169-5347(97)01257-3
-
Wood decomposition as influenced by invertebratesBiological Reviews of the Cambridge Philosophical Society 91:70–85.https://doi.org/10.1111/brv.12158
-
End-Permian (252 Mya) deforestation, wildfires and flooding—an ancient biotic crisis with lessons for the presentEarth and Planetary Science Letters 529:115875.https://doi.org/10.1016/j.epsl.2019.115875
-
Crustacean disparity through the Phanerozoic: comparing morphological and stratigraphic dataBiological Journal of the Linnean Society 65:455–500.https://doi.org/10.1111/j.1095-8312.1998.tb01149.x
-
Insects as drivers of ecosystem processesCurrent Opinion in Insect Science 2:26–32.https://doi.org/10.1016/j.cois.2014.06.004
Article and author information
Author details
Funding
Chinese Academy of Sciences (XDA19050101 XDB26000000)
- Bo Wang
National Natural Science Foundation of China (42125201 41688103)
- Bo Wang
Natural Scientific Foundation of Shandong Province (ZR2020YQ27)
- Jun Chen
Russian Science Foundation (21-14-00284)
- Evgeny Yan
Chinese Academy of Sciences (2020VCA0020)
- Edmund A Jarzembowski
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We are grateful to George Perry, Martin Fikáček, and two anonymous reviewers for invaluable comments that improved this manuscript. We thank AG Ponomarenko, RG Beutel, AP Rasnitsyn, AG Kirejtshuk, J Xue, H Xu, B Huang, and H Zeng for helpful discussion, and D Yang for reconstructions. BW thanks members of the palaeoentomological laboratory of the Palaeontological Institute (Russian Academy of Sciences) for their help during his visit to Moscow (2010).
Copyright
© 2021, Zhao et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 2,012
- views
-
- 321
- downloads
-
- 22
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Evolutionary Biology
Eyespot patterns have evolved in many prey species. These patterns were traditionally explained by the eye mimicry hypothesis, which proposes that eyespots resembling vertebrate eyes function as predator avoidance. However, it is possible that eyespots do not mimic eyes: according to the conspicuousness hypothesis, eyespots are just one form of vivid signals where only conspicuousness matters. They might work simply through neophobia or unfamiliarity, without necessarily implying aposematism or the unprofitability to potential predators. To test these hypotheses and explore factors influencing predators’ responses, we conducted a meta-analysis with 33 empirical papers that focused on bird responses to both real lepidopterans and artificial targets with conspicuous patterns (i.e. eyespots and non-eyespots). Supporting the latter hypothesis, the results showed no clear difference in predator avoidance efficacy between eyespots and non-eyespots. When comparing geometric pattern characteristics, bigger pattern sizes and smaller numbers of patterns were more effective in preventing avian predation. This finding indicates that single concentric patterns have stronger deterring effects than paired ones. Taken together, our study supports the conspicuousness hypothesis more than the eye mimicry hypothesis. Due to the number and species coverage of published studies so far, the generalisability of our conclusion may be limited. The findings highlight that pattern conspicuousness is key to eliciting avian avoidance responses, shedding a different light on this classic example of signal evolution.
-
- Evolutionary Biology
Studying the fecal microbiota of wild baboons helps provide new insight into the factors that influence biological aging.