Early evolution of beetles regulated by the end-Permian deforestation

  1. Xianye Zhao
  2. Yilun Yu
  3. Matthew E Clapham
  4. Evgeny Yan
  5. Jun Chen
  6. Edmund A Jarzembowski
  7. Xiangdong Zhao
  8. Bo Wang  Is a corresponding author
  1. State Key Laboratory of Palaeobiology and Stratigraphy, Nanjing Institute of Geology and Palaeontology and Center for Excellence in Life and Paleoenvironment, Chinese Academy of Sciences, China
  2. University of Chinese Academy of Sciences, China
  3. Institute of Vertebrate Paleontology and Paleoanthropology, Chinese Academy of Sciences, China
  4. Department of Earth and Planetary Sciences, University of California, Santa Cruz, United States
  5. Palaeontological Institute, Russian Academy of Sciences, Russian Federation
  6. Institute of Geology and Paleontology, Linyi University, China
  7. Department of Earth Sciences, Natural History Museum, United Kingdom

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.sa0

Introduction

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.

Examples of Permian beetles.

(A and B) Tshekardocoleidae, Moravocoleus permianus Kukalová, 1969, photograph and reconstruction. (C and D) Permocupedinae, Permocupes sojanensis Ponomarenko, 1969, photograph and reconstruction. (E) Tshekardocoleidae, Sylvacoleus richteri Ponomarenko, 1963, elytra photograph. (F) Taldycupedinae, Taldycupes reticulatus Ponomarenko, 1969, elytra photograph. Scale bars represent 1 mm.

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.

Figure 2 with 1 supplement see all
Diversity of Coleoptera from the Early Permian to Middle Triassic.

Natural taxa and mixed taxa (natural taxa and formal taxa) are counted at family, genus, and species levels separately. (A) Family-level diversity of natural taxa. (B) Family-level diversity of mixed taxa. (C) Genus-level diversity of natural taxa. (D) Genus-level diversity of mixed taxa. (E) Species-level diversity of natural taxa. (F) Species-level diversity of mixed taxa. Abbreviations: P1, Early Permian; P2, Middle Permian; P3, Late Permian; T1, Early Triassic; T2, Middle Triassic.

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).

Figure 3 with 1 supplement see all
Ecological shifts of Coleoptera from the Early Permian to Middle Triassic.

(A) Simplified phylogeny of Coleoptera from the Early Permian to Middle Triassic. Thick lines indicate the known extent of the fossil record. The branches representing stem groups are shown in red. The ‘dead clade walking’ pattern is symbolized by the dashed line. For details of the phylogenetic analysis, see Figure 3—figure supplement 1. (B) Genus percentage of xylophagous groups from the Early Permian to Middle Triassic. Yellow graded band represents the ‘coal gap’.

Figure 3—source data 1

Character state matrix for the phylogenetic analysis.

https://cdn.elifesciences.org/articles/72692/elife-72692-fig3-data1-v2.xlsx

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 18).

Figure 4 with 8 supplements see all
Morphospace comparisons of Coleoptera from the Early Permian to Middle Triassic.

(A) Morphospace three-dimensional (3D) plot ordinated by principal coordinates analysis (PcoA), maximum observable rescale distance (MORD) matrices, based on species-level dataset. (B and C) Disparity comparisons ordinated by PcoA, MORD matrices, based on species-level dataset, proxy by pov and sov. Abbreviations: pov, product of variance; sov, sum of variance; P1, Early Permian; P2, Middle Permian; P3, Late Permian; T1, Early Triassic; T2, Middle Triassic.

Figure 4—source data 1

Fossil character matrix for the morphospace analysis.

https://cdn.elifesciences.org/articles/72692/elife-72692-fig4-data1-v2.xlsx
Figure 4—source data 2

Result of permutation test for morphological disparity.

https://cdn.elifesciences.org/articles/72692/elife-72692-fig4-data2-v2.xlsx

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 protocol

We 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 protocol

In 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 protocol

We 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 47).

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 47) 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 15). 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 15). 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.

The following data sets were generated
    1. Zhao X
    2. Yu Y
    3. Clapham M
    4. Yan E
    5. Chen J
    6. Jarzembowski E
    7. Zhao X
    8. Wang B
    (2021) Dryad Digital Repository
    Dataset of Early evolution of beetles regulated by the end-Permian deforestation.
    https://doi.org/10.5061/dryad.7m0cfxpvd

References

  1. Book
    1. Grimaldi D
    2. Engel MS
    (2005)
    Evolution of the Insects
    Cambridge University Press.
    1. Ponomarenko AG
    (1969)
    Historical development of archostematan beetles
    Trudy Paleontologicheskogo Instituta AN SSSR 125:1–240.
    1. Ponomarenko AG
    (2003)
    Ecological evolution of beetles (Insecta: Coleoptera)
    Acta Zoologica Cracoviensia 46:319–328.
    1. Ponomarenko AG
    (2004)
    Beetles (Insecta, Coleoptera) of the late Permian and early Triassic
    Paleontological Journal 38:S185–S196.
    1. Ulyshen MD
    (2016) Wood decomposition as influenced by invertebrates
    Biological Reviews of the Cambridge Philosophical Society 91:70–85.
    https://doi.org/10.1111/brv.12158

Article and author information

Author details

  1. Xianye Zhao

    1. State Key Laboratory of Palaeobiology and Stratigraphy, Nanjing Institute of Geology and Palaeontology and Center for Excellence in Life and Paleoenvironment, Chinese Academy of Sciences, Nanjing, China
    2. University of Chinese Academy of Sciences, Beijing, China
    Contribution
    Data curation, Formal analysis, Investigation, Methodology, Software, Validation, Visualization, Writing - original draft, Writing - review and editing
    Competing interests
    No competing interests declared
  2. Yilun Yu

    1. University of Chinese Academy of Sciences, Beijing, China
    2. Institute of Vertebrate Paleontology and Paleoanthropology, Chinese Academy of Sciences, Beijing, China
    Contribution
    Formal analysis, Investigation, Methodology, Software, Validation, Visualization, Writing - original draft, Writing - review and editing
    Competing interests
    No competing interests declared
  3. Matthew E Clapham

    Department of Earth and Planetary Sciences, University of California, Santa Cruz, Santa Cruz, United States
    Contribution
    Investigation, Validation, Writing - review and editing
    Competing interests
    No competing interests declared
  4. Evgeny Yan

    Palaeontological Institute, Russian Academy of Sciences, Moscow, Russian Federation
    Contribution
    Investigation, Visualization, Writing - review and editing
    Competing interests
    No competing interests declared
  5. Jun Chen

    1. State Key Laboratory of Palaeobiology and Stratigraphy, Nanjing Institute of Geology and Palaeontology and Center for Excellence in Life and Paleoenvironment, Chinese Academy of Sciences, Nanjing, China
    2. Institute of Geology and Paleontology, Linyi University, Linyi, China
    Contribution
    Formal analysis, Investigation, Methodology, Writing - review and editing
    Competing interests
    No competing interests declared
  6. Edmund A Jarzembowski

    1. State Key Laboratory of Palaeobiology and Stratigraphy, Nanjing Institute of Geology and Palaeontology and Center for Excellence in Life and Paleoenvironment, Chinese Academy of Sciences, Nanjing, China
    2. Department of Earth Sciences, Natural History Museum, London, United Kingdom
    Contribution
    Investigation, Writing - review and editing
    Competing interests
    No competing interests declared
  7. Xiangdong Zhao

    1. State Key Laboratory of Palaeobiology and Stratigraphy, Nanjing Institute of Geology and Palaeontology and Center for Excellence in Life and Paleoenvironment, Chinese Academy of Sciences, Nanjing, China
    2. University of Chinese Academy of Sciences, Beijing, China
    Contribution
    Investigation, Writing - review and editing
    Competing interests
    No competing interests declared
  8. Bo Wang

    State Key Laboratory of Palaeobiology and Stratigraphy, Nanjing Institute of Geology and Palaeontology and Center for Excellence in Life and Paleoenvironment, Chinese Academy of Sciences, Nanjing, China
    Contribution
    Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Validation, Visualization, Writing - original draft, Writing - review and editing
    For correspondence
    bowang@nigpas.ac.cn
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-8001-9937

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).

Version history

  1. Received: August 1, 2021
  2. Preprint posted: October 13, 2021 (view preprint)
  3. Accepted: November 3, 2021
  4. Accepted Manuscript published: November 8, 2021 (version 1)
  5. Version of Record published: November 11, 2021 (version 2)

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

  • 1,797
    views
  • 286
    downloads
  • 18
    citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

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)

  1. Xianye Zhao
  2. Yilun Yu
  3. Matthew E Clapham
  4. Evgeny Yan
  5. Jun Chen
  6. Edmund A Jarzembowski
  7. Xiangdong Zhao
  8. Bo Wang
(2021)
Early evolution of beetles regulated by the end-Permian deforestation
eLife 10:e72692.
https://doi.org/10.7554/eLife.72692

Share this article

https://doi.org/10.7554/eLife.72692

Further reading

    1. Evolutionary Biology
    2. Immunology and Inflammation
    Mark S Lee, Peter J Tuohy ... Michael S Kuhns
    Research Advance

    CD4+ T cell activation is driven by five-module receptor complexes. The T cell receptor (TCR) is the receptor module that binds composite surfaces of peptide antigens embedded within MHCII molecules (pMHCII). It associates with three signaling modules (CD3γε, CD3δε, and CD3ζζ) to form TCR-CD3 complexes. CD4 is the coreceptor module. It reciprocally associates with TCR-CD3-pMHCII assemblies on the outside of a CD4+ T cells and with the Src kinase, LCK, on the inside. Previously, we reported that the CD4 transmembrane GGXXG and cytoplasmic juxtamembrane (C/F)CV+C motifs found in eutherian (placental mammal) CD4 have constituent residues that evolved under purifying selection (Lee et al., 2022). Expressing mutants of these motifs together in T cell hybridomas increased CD4-LCK association but reduced CD3ζ, ZAP70, and PLCγ1 phosphorylation levels, as well as IL-2 production, in response to agonist pMHCII. Because these mutants preferentially localized CD4-LCK pairs to non-raft membrane fractions, one explanation for our results was that they impaired proximal signaling by sequestering LCK away from TCR-CD3. An alternative hypothesis is that the mutations directly impacted signaling because the motifs normally play an LCK-independent role in signaling. The goal of this study was to discriminate between these possibilities. Using T cell hybridomas, our results indicate that: intracellular CD4-LCK interactions are not necessary for pMHCII-specific signal initiation; the GGXXG and (C/F)CV+C motifs are key determinants of CD4-mediated pMHCII-specific signal amplification; the GGXXG and (C/F)CV+C motifs exert their functions independently of direct CD4-LCK association. These data provide a mechanistic explanation for why residues within these motifs are under purifying selection in jawed vertebrates. The results are also important to consider for biomimetic engineering of synthetic receptors.

    1. Evolutionary Biology
    Robert Horvath, Nikolaos Minadakis ... Anne C Roulin
    Research Article

    Understanding how plants adapt to changing environments and the potential contribution of transposable elements (TEs) to this process is a key question in evolutionary genomics. While TEs have recently been put forward as active players in the context of adaptation, few studies have thoroughly investigated their precise role in plant evolution. Here, we used the wild Mediterranean grass Brachypodium distachyon as a model species to identify and quantify the forces acting on TEs during the adaptation of this species to various conditions, across its entire geographic range. Using sequencing data from more than 320 natural B. distachyon accessions and a suite of population genomics approaches, we reveal that putatively adaptive TE polymorphisms are rare in wild B. distachyon populations. After accounting for changes in past TE activity, we show that only a small proportion of TE polymorphisms evolved neutrally (<10%), while the vast majority of them are under moderate purifying selection regardless of their distance to genes. TE polymorphisms should not be ignored when conducting evolutionary studies, as they can be linked to adaptation. However, our study clearly shows that while they have a large potential to cause phenotypic variation in B. distachyon, they are not favored during evolution and adaptation over other types of mutations (such as point mutations) in this species.