Front. Ecol. Evol. Frontiers in Ecology and Evolution Front. Ecol. Evol. 2296-701X Frontiers Media S.A. 10.3389/fevo.2025.1596591 Ecology and Evolution Original Research Rapid radiations underlie most of the known diversity of life Wiens John J. 1 * Moen Daniel S. 2 1 Department of Ecology and Evolutionary Biology, University of Arizona, Tucson, AZUnited States 2 Department of Evolution, Ecology, and Organismal Biology, University of California, Riverside, Riverside, CAUnited States

Edited by: Diego Baldo, CONICET Institute of Subtropical Biology (IBS), Argentina

Reviewed by: Juan Diego Daza, Sam Houston State University, United States

Frederico Henning, Federal University of Rio de Janeiro, Brazil

*Correspondence: John J. Wiens, wiensj@arizona.edu

20 08 2025 2025 13 1596591 19 03 2025 27 06 2025 Copyright © 2025 Wiens and Moen 2025 Wiens and Moen

This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

Rapid radiations, including adaptive radiations, are of considerable interest to evolutionary biologists, in large part because they are thought to underlie much of the species diversity of life. Yet, this fundamental idea has only been tested at a limited scale, within frogs. Here, we test this idea across living organisms and within many of the largest clades (e.g. animals, plants). Specifically, we quantify how much of Earth’s species richness is contained within rapid radiations (clades with high net diversification rates). We find that among the major clades of living organisms and among land plant phyla and animal phyla, >80% of known species richness is contained within the few clades in the upper 90th percentile for diversification rates in each group. Thus, these exceptionally rapid radiations contain most of Earth’s extant species diversity. Patterns were broadly similar using smaller clades (orders, families) and in insects and vertebrates, with the majority of species generally contained within clades in the upper 75th percentile. Results were also similar using large-scale clades defined by their ages instead of taxonomic ranks. Overall, these results show for the first time that most of the known species richness of life is explained by rapid radiations. Moreover, phenotypic evidence from previous studies suggests that some of the most species-rich rapid radiations across life, animals, and plants may also qualify as adaptive radiations.

adaptive radiation biodiversity diversification macroevolution species richness section-in-acceptance Phylogenetics, Phylogenomics, and Systematics

香京julia种子在线播放

    1. <form id=HxFbUHhlv><nobr id=HxFbUHhlv></nobr></form>
      <address id=HxFbUHhlv><nobr id=HxFbUHhlv><nobr id=HxFbUHhlv></nobr></nobr></address>

      Introduction

      Adaptive radiation has become a central topic in modern evolutionary biology. Adaptive radiations are generally thought to involve rapid diversification of species (i.e. high rates of speciation minus extinction) accompanied by the evolution of ecologically relevant phenotypes (Schluter, 2000; Gillespie et al., 2020; Moen et al., 2021). One reason that adaptive radiations are of such great interest is that they are thought to underlie much of the species diversity and phenotypic diversity of life (Simpson, 1953; Schluter, 2000; Glor, 2010; Givnish, 2015). After all, if the idea of adaptive radiation applied to only a few exceptional clades with limited species diversity (e.g. African rift-lake cichlids, Caribbean anoles, Galápagos finches, sticklebacks), then it is unclear why adaptive radiation should be a broadly important topic in the field. However, with the exception of a recent study in frogs (Morinaga et al., 2023), no studies have explicitly quantified how much of Earth’s overall species richness is contained in clades resembling adaptive radiations.

      The study in frogs (Morinaga et al., 2023) developed a framework for addressing how much species diversity and phenotypic diversity is explained by clades with different dynamics of diversification and phenotypic evolution (following Moen et al., 2021). Thus, they classified clades as resembling adaptive radiations (those with rapid phenotypic change and rapid species diversification), non-adaptive radiations (those with rapid species diversification but unexceptional rates of phenotypic change), adaptive non-radiations (with rapid phenotypic change but unexceptional diversification rates), and non-adaptive, non-radiations (with relatively slow rates of phenotypic change and species diversification). They found that in frogs ~75% of both species richness and phenotypic diversity was contained in clades resembling adaptive radiations, with above-average rates of diversification and phenotypic change. Here and throughout, we refer to diversification rates as the rate of speciation minus the rate of extinction, or the rate of species accumulation over time. For brevity, we use “diversification rate” instead of “net diversification rate”.

      Here, we apply this general framework (Moen et al., 2021; Morinaga et al., 2023) to ask: how much of life’s current species diversity evolved from relatively rapid radiations? It would be very difficult to estimate comparable species-level rates of phenotypic change for hundreds of species in every major clade across life. Therefore, without phenotypic data and rates, one cannot distinguish adaptive radiations from non-adaptive radiations, adaptive non-radiations, and non-adaptive non-radiations. Nevertheless, we can quantify how much of life’s species richness is contained within relatively rapid radiations (i.e. clades with high diversification rates) as opposed to non-radiations (i.e. clades with lower diversification rates).

      Some previous analyses have found that variation in species richness of named clades of comparable rank (e.g. families, phyla) is strongly related to the diversification rates of these clades (Scholl and Wiens, 2016). However, this relationship alone says nothing about what proportion of species belong to rapidly diversifying clades as opposed to clades with more moderate diversification rates. For example, in both of the hypothetical examples in Figure 1 , diversification rates and species richness are very strongly related (Example A: r 2 = 0.98, P=0.002; Example B: r 2 = 0.93, P=0.008; details in Figure 1 , Supplementary Table S1 ). But in Example A, the most rapidly diversifying clade contains only 30% of the group’s species richness, whereas in Example B, the most rapidly diversifying clade contains 80%. Overall, the contribution of rapid radiations to life’s species diversity remains unknown, and is not addressed by the relationship between diversification rates and species richness.

      Two hypothetical examples, illustrating how the proportion of species in a group that are contained in the most rapidly diversifying clade(s) can vary. In Example A, the most rapidly diversifying clade (Clade 1) includes only 30% of the group’s overall richness. In Example B, the most rapidly diversifying clade (Clade 1) contains 80%. Our hypothesis is that, within a given taxonomic group, the majority of species will belong to the most rapidly diversifying clade (Clade 1), as in Example B This hypothesis follows from the idea that rapid radiations (such as adaptive radiations) contain most of the diversity of life. The alternative hypothesis is that species richness will be more evenly distributed among clades (Example A), such that only a minority of species are in the most rapidly diversifying clade (Clade 1). Importantly, the proportion of species in a clade that are contained within the most rapidly diversifying clade can be independent of the relationship between diversification rates and species richness among clades. Thus both of these examples show a strong relationship between species richness and diversification rates among clades (Example A: r 2 = 0.98, P=0.002; Example B: r 2 = 0.93, P=0.008), regardless of whether the most rapidly diversifying clade contains a minority (Example A) or majority (Example B) of the group’s overall species richness. For these examples, we assumed that all clades were 50 million years old and estimated diversification rates using the stem-group MS estimator with ϵ=0.5 (values in Supplementary Table S1 ).

      Example A shows five clades with species counts of 30, 25, 20, 15, and 10, respectively. Example B combines these into Clade 1 with 80 species and smaller clades with 9, 5, 4, and 2 species. Clade 1 in both examples is blue, with others in orange.

      Two factors now facilitate addressing this question. First, several recent studies have analyzed patterns of diversification rates and richness across all of life, and among major clades of land plants, animals, insects, and vertebrates (see Methods for details). Thus, many of the building blocks needed to address this question are now available. Second, recent studies have also shown that the clade-based estimator of diversification rates used in these studies (Magallón and Sanderson, 2001) yields estimates that are strongly related to true diversification rates in simulations (r 2~0.7), including simulations in which speciation and extinction rates vary within clades over time and among subclades (Meyer and Wiens, 2018; Meyer et al., 2018). Empirical analyses of randomly selected clades across life (Yu and Wiens, 2024; see also Supplementary Dataset S1, Supplementary Table S2 ) show that these clade-based estimates (Magallón and Sanderson, 2001) generally are strongly related (r 2~0.8) to estimates from a new Bayesian species-based method (ClaDS; Maliet et al., 2019; Maliet and Morlon, 2022). These two lines of evidence suggest that these clade-based estimates are accurate and a reasonable proxy for estimates from this Bayesian species-level method. These clade-based estimators can estimate diversification rates given only the age and species richness of each clade (Magallón and Sanderson, 2001). Therefore, a well-sampled, fully-resolved, time-calibrated, species-level phylogeny is not necessary to estimate diversification rates within each clade using this clade-based approach, unlike ClaDS and many other estimators based on species-level branch lengths. By using this clade-based estimator, we can conduct analyses across all of life and within many major groups that lack a comprehensive species-level phylogeny (i.e. most of them).

      In this study we quantify how much of life’s species richness is contained within relatively rapid radiations. We do this across living organisms and within some of the most species-rich groups, including land plants, animals, insects, and vertebrates. For each of these five groups, we compile data on diversification rates and species richness for named clades. We then determine what proportion of the richness of each group is contained within clades with above-average diversification rates (following Morinaga et al., 2023), and those with diversification rates in the upper 75th, 90th, and 95th percentiles. These new analyses allow us to address the long-standing question of whether most of the diversity of life is contained within rapid radiations.

      Materials and methods Data sources

      We first assembled data on species richness and diversification rates of named, ranked clades (e.g. kingdoms, phyla, classes, orders, families) from studies across all of life (Scholl and Wiens, 2016; Chen and Wiens, 2021), and among clades of land plants (Hernández-Hernández and Wiens, 2020), animals (Wiens, 2015a, Jezkova and Wiens, 2017), insects (Wiens et al., 2015), and vertebrates (Wiens, 2015b). These sources contained information on the richness, stem-group age, and estimated diversification rate for each named, ranked clade. We initially used this information directly for our calculations. We then explored the sensitivity of the results to new analyses based on alternative estimates of species richness, clade age, and the corresponding diversification rates ( Supplementary Appendix S1 ). All species numbers refer to extant, described species. The taxonomy used here (e.g. ranking of groups as phyla vs. classes) followed the taxonomy used in those papers. We performed separate analyses comparing clades of the same rank (e.g. kingdoms), across life and within selected major groups (animals, plants, insects, vertebrates). Some clades were nested inside others (families within kingdoms) but nested clades were not included in the same analysis ( Table 1 ).

      Percentage of species diversity contained within clades with high diversification rates.

      Group Rank Tree Total clades Total species Clades above average Percentage of species in high-diversification clades
      Above average 75th 90th 95th
      Life Kingdoms NA 17 2,068,138 6 99.0 91 90.5 16.7
      Life Families NA 2,545 1,409,134 941 77.4 58.5 32.4 19.2
      Plants Phyla FPU 10 306,976 5 98.0 94.8 90.6 90.6
      Plants Phyla FPC 10 306,976 5 98.0 90.8 90.6 90.6
      Plants Orders FPC 140 306,976 62 96.2 84.5 46.0 26.0
      Plants Orders FPU 140 306,976 55 92.0 69.4 49.4 26.0
      Plants Orders FPCM 140 306,976 61 96.6 84.5 59.2 39.2
      Plants Families FPC 678 306,976 262 94.0 82.6 43.7 32.8
      Plants Families FPU 678 306,976 257 91.9 74.5 54.6 30.8
      Animals Phyla Tree 1 28 1,515,954 12 98.4 97.8 86.5 1.6
      Animals Phyla Tree 2 28 1,515,954 13 99.2 97.3 86.5 1.6
      Animals Phyla Tree 3 28 1,515,954 14 99.3 97.8 89.4 1.6
      Animals Families NA 1,710 1,092,718 630 74.5 53.2 27.1 16.5
      Vertebrates Classes Tree 1 12 66,113 6 97.9 38.6 15.6 15.6
      Vertebrates Classes Tree 2 12 66,113 6 97.9 38.6 15.6 15.6
      Vertebrates Orders NA 144 58,636 63 79.2 63.9 44.6 39.8
      Vertebrates Families NA 778 51,763 299 81.1 61.4 26.4 10.4
      Insects Orders NA 31 1,063,532 15 96.2 71.0 29.9 29.6
      Insects Families NA 870 1,036,830 304 72.3 52.8 26.3 14.9

      For each group, we give the rank of the clades used, the specific tree used, the total number of clades of that rank included, the total number of species contained among those clades, the number of those clades with above-average diversification rates (based on the stem-group estimator with ϵ=0.5), and then the percentage of all species contained among clades with above-average diversification rates and among clades with diversification rates in the upper 75th, 90th, and 95th percentiles. FPC = constrained tree (Fiz-Palacios et al., 2011). FPU = unconstrained tree (Fiz-Palacios et al., 2011). FPCM = constrained tree (Fiz-Palacios et al., 2011), modified to reflect the tree of Magallón et al. (2015). NA, not applicable (only a single tree was analyzed).

      For some groups, two or more trees were used to estimate diversification rates in the original studies (see details below): different trees can lead to different clade ages and different rate estimates. We performed separate analyses on each tree, to address the sensitivity of the results to these differences in clade ages.

      Across life, we used the 17 higher-level clades (Supplementary Dataset S2) analyzed by Chen and Wiens (2021). The tree used in that study was from Parfrey et al. (2011). These clades encompass the groups usually ranked as kingdoms (bacteria, archaeans, fungi, animals, plants) and various clades interspersed among them that are usually referred to as protists and algae. These 17 clades collectively encompassed 2.068 million extant, described species (i.e. most of the 2.18 million species currently known; Bánki et al., 2024). We corrected errors in a few divergence dates reported in the supplementary materials in that study, such that the dates fully matched the methods given in that study. Note: all supplementary tables and appendices are available in the Supplementary Material , and all datasets (Supplementary Datasets S1–S27) are available on FigShare: https://figshare.com/s/09c264942aab86d1cc05.

      These 17 clades do not include the “Asgard archaea”, the recently discovered clades that are more closely related to eukaryotes than to other archaeans (Williams et al., 2013; Zaremba-Niedzwiedzka et al., 2017). We do not know of time-calibrated phylogenies that include them. However, those taxa encompass very few species (all described archaeans together include ~400 species; Bánki et al., 2024), and so their exclusion here should have minimal impact. Specifically, by excluding small clades with (presumably) low rates, it makes it less likely that most species will be in the few large clades with higher rates. The latter pattern is the pattern that we typically observed (see Results).

      We also used the 2,545 families across life (Supplementary Dataset S3) analyzed by Scholl and Wiens (2016). These clades collectively encompass only 1.409 million species, because generally only clades included in time-calibrated phylogenies could be used to estimate diversification rates. These families encompassed most species of animals and plants, and many families of archaeans, bacteria, protists, and fungi. We did not analyze classes or orders across life because their coverage is relatively incomplete in that dataset. However, we analyzed orders of land plants and insects along with classes and orders of vertebrates. Our coverage was more complete in these groups.

      For analyses across land plants, we initially used the data from Hernández-Hernández and Wiens (2020). These included all phyla (Supplementary Dataset S4), orders (Supplementary Dataset S5), and families (Supplementary Dataset S6), spanning 306,976 species. The ages were based on two time-calibrated trees, from Fiz-Palacios et al. (2011), which included all land plant families. These two trees account for the considerable uncertainty regarding the crown-group age of angiosperms, and span much of the range among most recent estimates (Sauquet et al., 2022). Specifically, in these two trees, the crown-group age of angiosperms ranges from 130–267 million years. One tree was called “Fiz-Palacios constrained” (FPC; constrained to the younger age) and the other “Fiz-Palacios unconstrained” (FPU). Because the trees had different ages for each clade, their diversification rates differed. For plant orders, we also used a third tree, which corresponded to the FPC tree but with the phylogeny in angiosperms following Magallón et al. (2015). We refer to this as the FPCM tree. We did not use the family-level version of that tree because it did not include all angiosperm families.

      For analyses across animals, we used the three phylum-level trees used by Jezkova and Wiens (2017), originally from Wiens (2015a). These trees differ in the topological and age constraints used. These trees included 28 of the ~32 known animal phyla and collectively encompassed 1,515,954 animal species (Supplementary Datasets S7–S9). The four unsampled phyla contain relatively few species (~200 in total; Wiens, 2015a).

      For animal families, we used the dataset of Scholl and Wiens (2016). This dataset included 1,710 families that encompass 1,092,718 species (Supplementary Dataset S10).

      For insect orders, we used the dataset from Wiens et al. (2015). This dataset included 31 hexapod orders, encompassing 1,063,532 species (Supplementary Dataset S11). This included most known species (1.0 million; Bánki et al., 2024). Those authors used three trees, but here we used the tree from Misof et al. (2014), which was based on the most extensive molecular dataset (in terms of markers). We refer to hexapods here as insects, although hexapods also include three relatively small non-insect clades (Collembola, Protura, and Diplura).

      For insect families, we used a portion of the family-level dataset from Scholl and Wiens (2016), which was originally from Rainford et al. (2014). We included data from 870 families (Supplementary Dataset S12), which encompassed 1,036,830 species (again, almost all known species).

      For major vertebrate clades, we initially used the dataset of Wiens (2015b), which included 12 major clades (Supplementary Dataset S13). These are largely equivalent to classes (Bánki et al., 2024), but with Reptilia separated into Lepidosauria (lizards, snakes, tuatara), Testudines (turtles), Aves (birds), and Crocodylia (crocodilians). Wiens (2015b) included two trees, which differed only in the placement of the two clades of cyclostomes (hagfishes [Myxini], lampreys [Petromyzontiformes]). These 12 clades together encompass 66,113 species. For vertebrate orders and families we used the dataset of Scholl and Wiens (2016). These taxa encompassed 58,636 species and 51,763 species, respectively.

      The data described above were used for our baseline analyses. We also performed sensitivity analyses using alternative, updated estimates of clade ages, species numbers, and diversification rates ( Supplementary Appendix S1 ). These results were generally very similar to the baseline results, strongly suggesting that the overall conclusions are robust to reasonable variation in these parameters.

      Diversification rates

      For diversification rates, we used the method-of-moments estimator (Magallón and Sanderson, 2001). We refer to this as the MS estimator hereafter. We focused on the MS estimator for stem-group ages:

      r = 1 / t *  log [ n ( 1 ϵ ) + ϵ ]

      where r is the estimated net diversification rate, t is the clade’s age, n is the clade’s species richness, and ϵ is the assumed ratio of speciation to extinction rates (see below). For a given clade (e.g. phylum), the diversification rate is calculated based on its richness and age, not based on a summary of rate estimates for included taxa within that clade. Nevertheless, previous analyses (Scholl and Wiens, 2016) show that diversification rates of higher taxa (e.g. kingdoms) are strongly related to the mean rates of the clades within them (e.g. families).

      We used the stem-group estimator because it allows clades to be included even when only one species from that clade is included in the tree, and is fully insensitive to incomplete phylogenetic sampling of species within clades (Meyer and Wiens, 2018) and largely insensitive to incomplete sampling among clades (Scholl and Wiens, 2016). By contrast, the crown-group estimator requires at least two species per clade, and is biased when taxon sampling does not encompass the crown-group age (i.e. the oldest split within the clade; Sanderson, 1996; Meyer and Wiens, 2018).

      Yu and Wiens (2024) tested if rate estimates from the MS estimator are significantly related to those from ClaDS (a Bayesian estimator based on species-level phylogenies; Maliet et al., 2019; Maliet and Morlon, 2022). They found a strong relationship between these rate estimates across and among groups, using the crown-group MS estimator (animals: r 2 = 0.73; plants: r 2 = 0.94; fungi: r 2 = 0.93; bacteria: r 2 = 0.92; archaeans: r 2 = 0.67; across life: r 2 = 0.86). We found generally similar results using the stem-group estimator, except in bacteria and archaea (animals: r 2 = 0.67; plants: r 2 = 0.89; fungi: r 2 = 0.88; bacteria: r 2 = 0.04; archaeans: r 2 = 0.37; across life: r 2 = 0.82; our Supplementary Table S2 ; Supplementary Dataset S1). Both ClaDS and the crown-group MS estimator focus on the age and relationships within the crown group, which might explain the weaker relationships between ClaDS and the stem-group estimator.

      The MS estimator includes a correction (ϵ) for the inclusion of only living clades in rate estimation (Magallón and Sanderson, 2001). This correction is the assumed ratio between extinction and speciation rates across the tree (Magallón and Sanderson, 2001). For the stem-group estimator, simulations show that the value of ϵ used does not strongly affect the accuracy of the estimated rates (despite extensive variation in speciation and extinction rates among clades; Meyer and Wiens, 2018). Similarly, empirical studies show that ϵ generally has little impact on relationships between rates and richness and other conclusions (e.g. Wiens, 2015a; Scholl and Wiens, 2016; Chen and Wiens, 2021; Morinaga et al., 2023). Therefore, for simplicity, we primarily present results using an intermediate value (ϵ=0.5). However, we also performed analyses assuming a higher value (ϵ=0.9). The results were generally similar and often identical for the question asked here.

      Some authors have criticized the MS estimator on the grounds that diversification rates estimated from fossils using this method do not necessarily predict later species richness in time series of fossils (Rabosky and Benson, 2021). Unsurprisingly, the MS estimator is unable to predict future mass extinctions and key innovations based on patterns of richness before these events occurred. However, the method is used to estimate diversification rates, not future richness. Therefore, their results do not address the usefulness of this method for estimating diversification rates. Moreover, those authors only subjected the MS estimator to this test. Thus, the MS estimator might actually be the best method according to this test, if other methods were included in the comparison. We note that the standard approach to evaluating diversification-rate estimators is instead using simulations in which the true diversification rates are known and one evaluates the correlation between the true and known rates. The MS estimator performs well by this standard criterion (i.e. with strong correlations between true and estimated rates; Meyer and Wiens, 2018; Meyer et al., 2018), and better than a widely used Bayesian species-level estimator (review in Wiens, 2024).

      Those authors (Rabosky and Benson, 2021) also implied that the MS estimator performs poorly when the true diversification rates are not constant. Yet, simulations repeatedly show that this method yields similarly strong relationships between true and estimated rates, regardless of whether rates are constant within clades, variable among clades, or variable within clades (Kozak and Wiens, 2016; Meyer and Wiens, 2018; Meyer et al., 2018). These simulations included variable rates among subclades within a clade and rates that changed over time within a clade, including linear and exponential increases and decreases in extinction rates and speciation rates. Furthermore, the MS estimator does not estimate separate speciation and extinction rates, and therefore does not suffer from the problem of the potential non-identifiability of these rates (Louca and Pennell, 2020). Again, our conclusions are largely robust to different assumptions about the ratio of extinction to speciation rates (ϵ).

      Another misconception is that strong relationships between diversification rates and richness are inevitable when using the MS estimator, because richness is used to calculate diversification rates. Simulations show that strong relationships between species richness and diversification from the MS estimator are not artifactual (Kozak and Wiens, 2016). Further, empirical analyses show that such strong relationships are not inevitable (Scholl and Wiens, 2016; Yu and Wiens, 2024). For example, relationships between richness and diversification rates are generally weak when clades are chosen randomly instead of comparing clades of the same rank (Yu and Wiens, 2024). Simulations and empirical analyses also show that faster diversification rates in younger clades can decouple richness and diversification rates (Kozak and Wiens, 2016; Scholl and Wiens, 2016). Therefore, there is no circularity associated with examining the relationships between diversification rates and richness of clades, nor with examining the percentage of species that belong to clades with high diversification rates (a percentage which can vary extensively even when there is a strong diversification-richness relationship; Figure 1 ). Finally, under standard definitions of circularity, the method determines the results (which is not true here) and the results are then used to justify the method (also not true).

      There has been debate about the extent to which diversification rates are dependent on the ages of clades (e.g. Henao-Diaz et al., 2019; Louca et al., 2022). However, the more important point is that the method used here gives strong relationships between true and estimated rates, even if there can be significant relationships between clade age and diversification rates (but not always; e.g. Yu and Wiens, 2024). Furthermore, if faster diversification rates were solely a function of clade age, then there is no reason to expect a few clades with exceptional diversification rates to contain most species richness. This latter pattern is the focus of our paper, not the relationship between diversification rates and richness. Named clades of the same rank do not have identical ages, but instead can have higher variance in clade ages than randomly selected clades, at least for the ranks examined here (families and above; Yu and Wiens, 2024).

      Statistical analyses

      For each group and each taxonomic rank (e.g. animal phyla; Table 1 ), we examined the distribution of diversification rates among clades and determined the mean and the upper 75th, 90th, and 95th percentiles. We then summed the number of species contained in the clades with above-average diversification rates, and calculated the percentage of the total number of species (all clades summed) that belonged to the above-average clades. We followed the same procedure for the upper 75th, 90th, and 95th percentiles.

      Based on past statements about the importance of adaptive radiations for species diversity (see Introduction), we predicted that the majority of species in each group would belong to clades with above-average diversification rates (i.e. rapid radiations). For example, Morinaga et al. (2023) found that ~75% of frog species belonged to clades with above-average rates.

      Note that this question is not based on a statistical comparison to a simulated null distribution. Instead, this question is about empirical patterns. Deviations in the observed richness of clades relative to simulated constant diversification rates or randomly varying diversification rates are not the question of interest here. Simulations show that variation in diversification rates can drive variation in richness among clades (Kozak and Wiens, 2016), but we do not know how much richness belongs to the most rapidly diversifying clades in empirical datasets. That is our question here.

      Similarly, there is empirical evidence that diversification rates vary across living organisms and that this variation seems to underlie variation in species richness (e.g. Scholl and Wiens, 2016). Therefore, comparison to a null model of equal rates among clades (or over time) would be pointless: we already know that rates are not constant. But the fact that rates vary among clades does not tell us how species richness will be distributed among those clades ( Figure 1 ).

      Potential methodological issues

      Here we describe several potential methodological issues, and the additional analyses performed to address them.

      Use of named clades

      The use of named clades can be controversial (Poe et al., 2021; Baker et al., 2021). For our study, it was essential to use named clades because phylogenies do not yet include all species in every group. Using named clades allowed us to assign species to clades without including every species in the tree. To address the sensitivity of the results to the particular set of clades used, we performed analyses for most groups using different taxonomic ranks (e.g. phyla vs. families). However, we do not expect the results to be identical when using higher taxonomic ranks (e.g. kingdoms, phyla) vs. lower-ranked clades (e.g. orders, families). Specifically, when using higher-ranked clades, richness may be aggregated into a smaller number of clades. When using lower-ranked clades, richness may be more evenly distributed among clades (because large clades will be divided into many subclades).

      There are also different ways to divide species among major clades. We explored the impacts of making taxonomic groups based on the ages of clades at the largest taxonomic scales (kingdoms, plant phyla, animal phyla). Thus, the clades were aggregated based on their minimum ages rather than based on taxonomy alone. We describe the details of how this was done in Supplementary Appendix S2 . Note that if we defined clades to have identical ages, then the diversification rates of clades would depend strongly on their species richness. We emphasize that these were merely alternative analyses designed to test the robustness of the results to alternative clade divisions.

      We acknowledge that using named clades of the same rank (or unranked clades of similar age) influences the results relative to using randomly selected clades. Recent analyses (Yu and Wiens, 2024) demonstrated that when using randomly selected clades, the richness of clades tends to be more strongly related to their ages than to their diversification rates. Here, we focused on comparing clades of the same rank or similar age, such that there is the potential for richness to be divided more evenly among these clades (rather than biasing our results in favor of larger disparities). For example, when clades are selected randomly, most clades are very young and have very limited richness relative to older clades (Yu and Wiens, 2024).

      Non-comparability of species

      Another potential issue is the idea that species are not comparable across life. However, across most living organisms, species are typically united by gene flow with their conspecifics (Hernández-Hernández et al., 2021). Although bacteria might appear to be an exception, recombination occurs between closely related individuals in >90% of examined bacterial species (Diop et al., 2022). Given this, bacterial species may be generally comparable to species in other groups (Bobay and Ochman, 2017). There are also asexual eukaryote species, but these solely or predominantly asexual groups appear to have limited species diversity (Chen and Wiens, 2021). Nonetheless, to partially address this issue, we performed a set of analyses using only eukaryotic species.

      We note that our study relies on previous taxonomy and species delimitation. We recognize that diverse lines of evidence can be used to recognize species, and that most presently described species have traditionally been recognized based on morphological evidence. We think that morphological data should provide a minimum estimate of species numbers, even if these morphology-based species contain many morphologically cryptic species. We addressed how large numbers of undescribed species would impact our conclusions with focused analyses (see next section). Importantly, the larger projections of undescribed species numbers explicitly incorporate morphologically cryptic species (Larsen et al., 2017).

      Described vs. undescribed species

      The most important issue for our conclusions may be that many species remain undescribed. These undescribed species could number in the millions, billions, or trillions (Mora et al., 2011; Locey and Lennon, 2016; Larsen et al., 2017; Li and Wiens, 2023; recent review in Wiens, 2023).

      Chen and Wiens (2021) estimated richness and diversification rates among 17 kingdom-level clades using projections of undescribed richness (largely from Larsen et al., 2017). We also performed analyses using these projections, including both “low” projected richness (282 million species across life) and high projected richness (2.2 billion). The lower value might still be considered unrealistically high, but these analyses were intended to address the robustness of our conclusions to species numbers that are very different from current numbers of described species. The number of bacterial species is especially uncertain (Wiens, 2023), so we included analyses both across all kingdoms and among eukaryotes only.

      Incomplete sampling of clades and species

      Our sampling of higher-level clades is nearly complete (e.g. kingdoms, plant and animal phyla), such that the sampled clades incorporate most known species within each group. However, some analyses of lower-level clades did not include every named taxon (i.e. they were not present in the same time-calibrated tree as the others) and so did not include all species (e.g. families across life [~70% of species included], animals [~67%], and vertebrates [~78%]). The impact of this sampling on our analyses is difficult to address. However, our conclusions from families are similar across groups, regardless of whether those groups have relatively complete sampling of families (e.g. plants, insects) or less complete sampling (e.g. life, animals). Therefore, it seems unlikely that this factor explains our results.

      Variation in rates within clades

      We appreciate that some readers might be skeptical of our results because diversification rates vary within clades. Our measure of diversification for each clade reflects the overall rate within that clade, which (in effect) averages over time and among subclades within each clade. Diversification rate is a continuous variable, and we expect it to vary extensively among subclades and over time. Furthermore, some clades may have an overall fast rate primarily because of an exceptional subclade with a fast rate (e.g. angiosperms within plants). But this general pattern could apply to any clade at any level. We dealt with this issue primarily by exploring patterns at different taxonomic levels (e.g. among kingdoms, phyla, orders, and families). It would be problematic to selectively include and exclude subclades based on variation in their rates. Furthermore, we note again that diversification rates of higher taxa tend to be strongly related to the mean rates of the clades within them (Scholl and Wiens, 2016).

      Results

      Among known species, we found that clades with above-average diversification rates generally encompassed >70% of the species in each major clade ( Table 1 ). Furthermore, when considering higher-level clades (kingdoms, phyla), those clades with exceptional diversification rates (upper 90th percentile, 10% of the clades) generally included >85% of the species in each group ( Figures 2 , 3 ; Table 1 ). However, the clades with the fastest diversification rates (upper 95th percentile) often contained only a minority of species. When clades were subdivided into lower ranking taxa (e.g. families), richness was more evenly distributed among clades ( Table 1 ), but the majority of species in each group generally belonged to those clades in the upper 75th percentile for diversification rates ( Figure 3 , Table 1 ). We describe these results in detail for each taxonomic group below.

      Phylogeny, divergence times, and species richness of major clades across life, animals, and plants. For each group, we show a time-calibrated phylogeny and the species richness of each higher taxon. The species richness of clades that are in the upper 90th percentile for diversification rates (rapid radiations) are indicated with blue, whereas richness for other clades are shown in orange. Many of the clades with lower diversification rates are barely visible or not visible because they have relatively low species richness (note that all rapid radiations are visible). However, not every group with relatively high richness has rapid diversification rates (e.g. fungi, chordates, mollusks, bryophytes). We did not use a log scale here, in order to better show the relative differences in richness among groups. The tree for kingdom-level clades across life is from Chen and Wiens (2021; modified from Parfrey et al., 2011), the tree for animal phyla is from Wiens (2015a; Tree 2), and the tree for land plant phyla is the constrained tree from Fiz-Palacios et al. (2011). The data on species richness, clade ages, and diversification rates used here are given in Supplementary Dataset S2 (life), Supplementary Dataset S4 (plants), and Supplementary Dataset S8 (animals, Tree 2). Images of plants and animals are from John J. Wiens.

      Phylogenetic tree diagram displaying life, animals, and plants, categorized by their evolutionary relationships and number of species. Rapid radiations are highlighted in blue, and other clades in orange. Images of a butterfly, clownfish, and plant represent different groups. Time is shown in billions of years ago on the horizontal axis.

      The percentage of species belonging to rapid radiations and to other clades, across life and across plants and animals. The top row shows that most known species richness belongs to clades that are in the upper 90th percentile for diversification rates. The bottom shows that species are more evenly distributed among families, but that the majority of species nevertheless belong to clades with relatively rapid diversification rates (upper 75th percentile). The results for plants are based on the FPU tree (unconstrained tree of Fiz-Palacios et al., 2011). Data for all trees, taxonomic ranks, and diversification cutoffs are given in Table 1 .

      Six pie charts showing the percentage of species in rapid radiations versus other clades. Top row: Kingdoms, Plant phyla, Animal phyla; 90th percentile rapid radiations in blue. Bottom row: Families, Plant families, Animal families; 75th percentile rapid radiations in blue. Orange represents other clades.

      Across 17 kingdom-level clades spanning all of life ( Figure 2 ; Table 1 ; Supplementary Dataset S2), we found that the six clades with above-average diversification rates encompassed 99% of known species, and the two clades in the upper 90th percentile (animals and land plants) included 90% of known species ( Figure 3 ). However, the most rapidly diversifying of these clades (land plants) included only 17% of known species. Among these 17 clades, land plants are much younger than animals but have higher diversification rates and less richness. When these clades were subdivided into families (Supplementary Dataset S3), the majority of species (58%) were in clades in the upper 75th percentile for diversification rates, but with only a minority of species among those clades in the upper 90th and 95th percentiles ( Table 1 ).

      Among the 10 phyla of land plants ( Figure 2 ; Table 1 ; Supplementary Dataset S4), 98% of species belonged to clades with above-average diversification rates, and 91% belonged to a single clade (angiosperms, Magnoliophyta) with the fastest diversification rate (95th percentile). Among 140 land plant orders ( Figure 3 ; Table 1 ; Supplementary Dataset S5), 92–97% of species belonged to orders with above-average diversification rates (ranges refer to different trees), whereas 69–84% belonged to those in the upper 75th percentile and 46–59% in the upper 90th percentile. Among the 678 land plant families ( Table 1 ; Supplementary Dataset S6), 92–94% of species were in families with above-average diversification rates, 74–83% were in those with high diversification rates (upper 75th percentile; Figure 2 ), and 44–55% were in those with exceptional diversification rates (upper 90th percentile).

      Among 28 phyla of animals ( Figure 2 ; Table 1 ; Supplementary Datasets S7–S9), 97–98% of species belonged to phyla with high diversification rates (75th percentile) and 86–89% were in phyla with exceptional diversification rates (90th percentile; Figure 3 ). These patterns are largely explained by the relatively rapid diversification and numerical dominance of arthropods ( Figure 2 ). Among 1,710 animal families ( Table 1 ; Supplementary Dataset S10), 75% of included species belonged to clades with above-average diversification rates and 53% belonged to those with high diversification rates (upper 75th percentile; Figure 3 ). Only a minority of species (16–27%) belonged to families with exceptional diversification rates (upper 90th–95th percentiles).

      Among 31 insect (hexapod) orders ( Table 1 ; Supplementary Dataset S11), 96% of species belonged to those with above-average diversification rates, and 71% belonged to those with high diversification rates (75th percentile; Figure 4 ). Only a minority of species belonged to orders with exceptional diversification rates (90th and 95th percentile). Among 870 insect families ( Table 1 ; Supplementary Dataset S12), 72% of species belonged to families with above-average diversification rates, and 53% belonged to families with high diversification rates (75th percentile).

      The percentage of species richness belonging to rapid radiations and other clades, among land plant orders, insect orders, and vertebrate classes. The top row shows that most known species richness belongs to clades that are in the upper 75th percentile for diversification rates for plant and insect orders, but not for vertebrate classes. The bottom row shows that most known species in these groups do belong to clades with above-average diversification rates. The results for plants are based on the FPU tree (unconstrained tree of Fiz-Palacios et al., 2011). Data for all trees, taxonomic ranks, and diversification cutoffs are given in Table 1 .

      Six pie charts compare the percentage of species in rapid radiations and other clades across plant orders, insect orders, and vertebrate classes. The upper charts show the 75th percentile with more species in clades than rapid radiations, while the lower charts depict above-average species, with most in rapid radiations. A legend indicates blue for rapid radiations and orange for other clades.

      Among the 12 major clades of vertebrates (classes or equivalent; Table 1 ; Supplementary Dataset S13), 98% of species belonged to clades with above-average diversification rates ( Figure 4 ), but only a minority of species (16–39%) belonged to clades with high and exceptional diversification rates (75th to 95th percentiles; Figure 4 ). This pattern occurred because almost half of vertebrate species belonged to the clade Actinopterygia (ray-finned fishes), which has diversification rates that are merely above-average and not high or exceptional (Supplementary Dataset S13). For 144 vertebrate orders (Supplementary Dataset S14), 79% of species were in orders with above-average rates, and 40–64% belonged to orders with high and exceptional rates (75th to 95th percentiles). Among 778 vertebrate families (Supplementary Dataset S15), 81% of species richness was in families with above-average rates, 53% was in families with high rates (75th percentile), and a minority of species belonged to families with exceptional rates (90th and 95th percentiles).

      These baseline results utilized the MS estimator with ϵ=0.5. Results were similar using ϵ=0.9 ( Supplementary Table S3 ). The most notable difference was among animal phyla: using ϵ=0.9, most species richness (83%) was in the clade with an exceptionally fast diversification rate (Arthropoda; upper 95th percentile), but using ϵ=0.5 the clade with the fastest rate (Nematoda) contained relatively little known richness (2%).

      We then revisited the higher-level analyses with new estimates of diversification based on new estimates of species richness ( Table 2 ; Supplementary Appendix S1 ; animal phyla: Supplementary Datasets S16–S18; insect orders: Supplementary Dataset S19), and new estimates of richness combined with new estimates of clade ages (vertebrate classes: Supplementary Dataset S20; land plant phyla: Supplementary Dataset S21). The results were very similar, in terms of how much richness was contained within rapidly diversifying clades. The only substantive difference was (again) that in some analyses of animal phyla with revised richness data, the clades with exceptional rates (95th percentile) included 78% of the species (for Trees 1 and 2), instead of 2%.

      Percentage of species diversity contained within clades with high diversification rates, using alternative estimates of species numbers and clade ages.

      Group Rank Tree Total clades Total species Clades above average Percentage of species in high-diversification clades
      Above average 75th 90th 95th
      Land plants Phyla new 10 299,162 6 99.4 92.3 88.4 88.4
      Animals Phyla Tree 1 28 1,523,141 12 97.8 96.1 87.8 77.8
      Animals Phyla Tree 2 28 1,523,141 13 99.1 95.2 80.6 77.8
      Animals Phyla Tree 3 28 1,523,141 14 99.2 95.9 87.8 1.2
      Vertebrates Classes new 12 71,408 6 98.0 39.7 14.9 14.9
      Insects Orders NA 31 989,369 16 95.5 71.9 34.4 34.1

      See Supplementary Appendix S1 for details about alternative estimates. For each group, we give the rank of the clades used, the specific tree used, the total number of clades of that rank included, the total number of species contained among those clades, the number of those clades with above-average diversification rates (based on the stem-group estimator with ϵ=0.5), and then the percentage of all species contained among clades with above-average diversification rates and among clades with diversification rates in the upper 75th, 90th, and 95th percentiles. NA, not applicable (only a single tree was analyzed). new=a new tree was analyzed, different from those used in the baseline analyses.

      We also examined the impact of using projected species numbers across life ( Table 3 ; Supplementary Dataset S22), instead of only described species. We found that only 22–23% of projected species belonged to clades with above-average diversification rates. This disparate pattern arose because in these projections ~77% of all species are bacteria, which have below-average diversification rates, given their ancient age relative to other clades.

      Percentage of species diversity contained within clades with high diversification rates, using projected species numbers instead of numbers of described species.

      Group Rank Total clades Total species Clades above average Percentage of species in high-diversification clades
      Above average 75th 90th 95th
      Life (low) Kingdom 17 281,967,821 6 22.7 15.4 7.4 7.2
      Eukaryotes (low) Kingdom 15 63,967,444 6 99.9 68.1 32.4 31.9
      Life (high) Kingdom 17 2,238,367,821 6 22.0 22.0 7.3 7.3
      Eukaryotes (high) Kingdom 15 492,367,444 6 99.9 99.9 33.2 33.1

      We used two projections for overall species richness, one that is relatively low and one that is relatively high. For each group, we give the rank of the clades used, the total number of clades of that rank included, the total number of species contained among those clades, the number of those clades with above-average diversification rates (based on the stem-group estimator with ϵ=0.5), and then the percentage of all species contained among clades with above-average diversification rates and among clades with diversification rates in the upper 75th, 90th, and 95th percentiles.

      Considering only eukaryotes ( Table 3 ; Supplementary Dataset S22), then among the 15 kingdom-level clades, 99.9% of projected species belonged to those with above-average diversification rates and 68.1–99.9% belonged to those with high diversification rates (75th percentile; range represents low and high richness projections). Only a minority of projected species belonged to clades with exceptional diversification rates (90th and 95th percentiles).

      Additionally we examined the impact of defining the largest-scale clades based on clade ages rather than taxonomic ranks alone ( Table 4 ). Among kingdoms, we used 14 clades that were each at least 1 billion years old (range=1.08–4.20 billion years; Supplementary Dataset S23) and then 9 clades that were at least 1.5 billion years old (range=1.55–4.20 billion years; Supplementary Dataset S24). Using this approach, almost all species richness (98–99%) was in clades with above-average or high diversification rates (75th percentile). With these 14 clades, 74% of species richness was in clades with exceptional diversification rates (90–95th percentile), but only 18% was when using the 9 clades that were >1.5 billion years old (similar to the original results). Among land plant phyla ( Table 4 ), we used 6 clades that were at least 400 million years (Myr) old (range=459–556 Myr or 431–511 Myr depending on the tree; Supplementary Dataset S25): the results were almost identical to those using 10 phyla (with 91% of richness contained among clades with exceptional diversification rates; 90th–95th percentile). Among animal phyla ( Table 4 ; Supplementary Datasets S26–S27), we used seven clades that were at least 700 Myr old (Tree 1: range=757–1,031 Myr; Tree 2: 703–876). Results were similar to those based on 28 phyla, except that 94% of richness was contained in clades with exceptional diversification rates (90th–95th percentile). In summary, using clades defined based on clade ages gave similar results to the original analyses based on named clades, but with the most rapidly diversifying clades often encompassing more richness.

      Percentage of species diversity contained within clades with high diversification rates, using age-based clade definitions.

      Group Rank Tree Total clades Total species Clades above average Percentage of species in high-diversification clades
      Above average 75th 90th 95th
      Life 1.0 Bya NA 14 2,068,138 5 99.4 99.0 73.8 73.8
      Life 1.5 Bya NA 9 2,068,138 3 99.3 98.0 17.6 17.6
      Plants 400 Mya FPU 6 306,976 3 98.0 90.9 90.9 90.9
      Plants 400 Mya FPC 6 306,976 3 98.0 90.9 90.9 90.9
      Animals 700 Mya Tree 1 7 1,515,954 3 99.4 98.5 93.5 93.5
      Animals 700 Mya Tree 2 7 1,515,954 4 99.9 98.5 93.5 93.5

      For each group, we give the age of the clades used (Bya=billion years ago; Mya=million years ago), the specific tree used, the total number of clades of that rank included, the total number of species contained among those clades, the number of those clades with above-average diversification rates (based on the stem-group estimator with ϵ=0.5), and then the percentage of all species contained among clades with above-average diversification rates and among clades with diversification rates in the upper 75th, 90th, and 95th percentiles.

      Discussion

      Rapid radiations are thought to contain most of life’s species diversity, but this idea has not been tested explicitly. Here, we find that most known species belong to a small fraction of higher-level clades with relatively rapid diversification rates. We show this pattern across the major clades of life and among the major clades of plants and animals ( Figures 2 , 3 ; Table 1 ).

      These results are not apparent from previous studies nor do they appear to be methodological artifacts. Given previous analyses showing a relationship between diversification rates and species richness across life (Scholl and Wiens, 2016), one might argue that it is inevitable that those clades with more species would have higher diversification rates. But the patterns described here depend on how species are divided among clades with different diversification rates, not the relationship between diversification rates and species richness among clades (see example in Figure 1 ). We show that at the highest taxonomic levels (kingdoms, phyla), these rapidly diversifying clades (upper 90th percentile) contain >85% of all species ( Table 1 ). This pattern is also not obvious because the clades with the very fastest diversification rates (95th percentile) need not contain the majority of species ( Table 1 ). For example, among kingdom-level clades, plants have faster diversification rates than animals, but plants contain only a small minority of known species. These patterns of disparity in richness among clades are also reduced using lower taxonomic ranks (e.g. families). When using lower taxonomic ranks, the large, rapidly diversifying clades are subdivided into smaller clades that must each contain a smaller proportion of the group’s overall species richness. Nevertheless, even when using these lower taxonomic ranks, a minority of rapidly diversifying clades (upper 75th percentile) typically contains the majority of a group’s species richness ( Table 1 ; Figure 2 ). These overall patterns were generally robust to alternative estimates of clade ages, richness, and diversification rates ( Tables 2 , 3 ; Supplementary Table S3 ), and the results for families were similar regardless of whether sampling of families in the group is largely complete or spanned ~67–78% of the group’s species ( Table 1 ). We obtained similar results using ages to define clades, rather than using named ranked taxa ( Table 4 ). Thus, our results are not simply contingent on how humans divide the Tree of Life into named clades (see also Table 1 , showing results from alternative divisions of species into clades). Furthermore, the estimator of diversification rates used here gives estimates that are strongly related to the true rates (in simulations) and to a Bayesian species-based estimator (ClaDS) in empirical analyses. Thus, there is no basis for claiming that this method is problematic. Our goal was not to test the overall relationship between diversification rates and richness, but even if it were, such relationships are neither inevitable nor artifactual (Kozak and Wiens, 2016) nor are they circular (see Methods). Indeed, other empirical studies have suggested that diversification rates often do not underlie richness patterns (McPeek and Brown, 2007; Rabosky et al., 2012; Hedges et al., 2015; Yu and Wiens, 2024), further showing that our results are not obvious. On the other hand, our goal was to empirically test a long-standing idea in the macroevolutionary literature, and so our results (which support that idea) are not shocking from that perspective.

      The most important sensitivity of the results related to the inclusion of projected, undescribed bacterial species ( Table 3 ). At present, there are only ~10,000 described species of bacteria (Bánki et al., 2024). Among recent estimates of actual bacterial richness, even lower estimates are in the millions (Louca et al., 2019) and upper estimates are in the billions or trillions (Locey and Lennon, 2016). Because bacteria are very old relative to other major clades, they would still have relatively low overall diversification rates even if they contained billions or trillions of species (see calculations in Scholl and Wiens, 2016). Thus, if actual bacterial richness really is much higher than described richness for other groups, then a clade with low diversification rates would contain the majority of species across life (unlike our other results). Therefore, we caution that our results apply primarily to known species diversity. Nevertheless, if we restrict ourselves to eukaryotes, then the majority of species belong to clades in the upper 75th percentile for diversification rates, even if there are tens of millions of undescribed eukaryotic species. We also acknowledge that these projections of undescribed richness could be incorrect, which is why we treat our analyses of known diversity as the baseline.

      Our focus here has been on how much species richness is contained within rapid radiations, but how many of these rapid radiations might be adaptive radiations? There is considerable debate about the definition of adaptive radiation (e.g. Glor, 2010; Givnish, 2015; Hernández-Hernández, 2019; Gillespie et al., 2020; Moen et al., 2021). Some standard definitions require relatively rapid speciation within a clade, along with evidence for correlations between phenotypic traits and environments among species and the functional utility of the traits (Schluter, 2000; Gillespie et al., 2020). Many groups analyzed here may meet this broad criterion. Previous analyses have identified specific phenotypic traits and/or key innovations that seem to drive rapid diversification in these groups. Across life, a key trait that seems to explain rapid diversification rates is multicellularity (Chen and Wiens, 2021), which may have paved the way for much of the phenotypic variability (e.g. different cells and tissue types) in plants, animals, and fungi (Carroll, 2001). Within land plants, insect pollination seems to be a key trait that helps explain the rapid diversification of angiosperms relative to other land plant clades (Hernández-Hernández and Wiens, 2020), as suggested by many previous authors (e.g. Raven, 1977; Stebbins, 1981; Niklas, 2016). Insect pollination may drive considerable phenotypic diversification in flowers (Grant, 1949; Harder and Barrett, 2006) and much plant speciation (Grant, 1949; Sargent, 2004; Kay et al., 2006; van der Niet et al., 2014). In animals and vertebrates, terrestriality seems to be an important driver of diversification among major clades (Wiens, 2015a; 2015b; Jezkova and Wiens, 2017), and the invasion of land may represent a key innovation that allowed access to many open ecological niches (May, 1994; Benton, 2001; Vermeij and Grosberg, 2010). Analyses of diversification rates among insect orders suggest that herbivory may represent a key innovation in some insect clades (Wiens et al., 2015), a longstanding idea in evolutionary ecology (Ehrlich and Raven, 1964; Mitter et al., 1988; Farrell, 1998; Mayhew, 2007; Futuyma and Agrawal, 2009). Host shifts may be a particularly important driver of herbivory-related diversification in insects (e.g. Hardy and Otto, 2014; Forbes et al., 2017). In summary, previous studies identified key traits associated with many of these rapid radiations, and so these particular large clades might represent adaptive radiations (relative to comparable clades in these same groups that did not rapidly diversify). Conversely, the absence of these key traits may help why the majority of clades did not rapidly diversify in each of these groups.

      On the other hand, these rapidly diversifying clades may not fit every definition of adaptive radiation. Some definitions require exceptional rates of diversification and phenotypic evolution relative to other clades (e.g. Poe et al., 2018; Moen et al., 2021), but this may be too restrictive. Thus, Morinaga et al. (2023) calculated the relative richness of clades that were “adaptive-radiation-like”, having only above-average rates of diversification and phenotypic evolution. In our study, most species belonged to clades with above-average diversification rates, but not necessarily those in the upper 95th percentile ( Table 1 ). Furthermore, to better test for adaptive radiation, future analyses might also incorporate rates of phenotypic evolution. However, it may be challenging to estimate phenotypic rates from species-level phylogenies from many species across life.

      Our results cover all of life and several major groups (plants, animals, insects, vertebrates), but the pattern documented here may also occur within smaller groups. For example, many vertebrate groups are numerically dominated by a single relatively recent clade, such as teleosts within actinopterygian fishes, frogs within amphibians, placentals within mammals, neognaths within birds, and colubroids within snakes (trees from Irisarri et al., 2017; richness data from Bánki et al., 2024). Given that these clades are relatively recent and have high richness, they should have relatively high diversification rates. Similarly, within angiosperms, most species belong to the clade uniting the “core angiosperms” and not the species-poor clades at the base of the angiosperm tree (Magallón et al., 2015; Zuntini et al., 2024).

      Similar patterns may also occur in groups that we did not examine in detail here. For example, among the 154,534 currently described species of extant fungi (Bánki et al., 2024), most species (98%) belong to the phyla Ascomycota (98,334 species) and Basidiomycota (52,979 species), and not the other 9 fungal phyla. Phylogenies of fungi show these two phyla as sister taxa that are relatively young among higher-level fungal clades (James et al., 2006; Li et al., 2021). Given their relatively young age and high species richness, these two phyla would almost certainly represent rapid radiations relative to other fungal phyla (higher rates were shown in Scholl and Wiens (2016) but relatively few fungal phyla were included). We acknowledge that we have not confirmed the patterns that we found here within every group across life and at every phylogenetic scale, but we have shown these patterns at the largest phylogenetic scales and within some of the largest groups (e.g. animals, plants). This was our initial question. If this pattern does not hold within most genera or families (for example), it would not overturn our conclusions, which are specifically focused at larger scales.

      We emphasize that our study is focused on patterns of extant species richness. We recognize that the relative richness of many clades has changed dramatically over geological time. For example, angiosperms began as a single species, and initially diversified when many other plant clades were already extant and relatively diverse (Knoll, 1986; Niklas, 2016). Moreover, many groups with low diversity today had higher diversity in the past (e.g. coelacanths; Torino et al., 2021), and this might also help explain some of the disparities in richness observed among clades today (especially why some ancient groups are species poor). These observations do not make the patterns described here incorrect. Within a group, the fact that some clades have declined does not invalidate other clades as rapid radiations: diversification includes extinction, and “rapid” is relative and pertains to the timescale of any given comparison. Instead, we caution that these patterns may apply primarily to present day richness patterns. Although it would be interesting to study these patterns in the fossil record for some exemplar groups, it would be difficult to conduct similar analyses across life based on paleontological data.

      Many authors have suggested that adaptive radiations are characterized by declining diversification rates over time. We find that much of life (at multiple scales) is numerically dominated by relatively young clades with rapid diversification rates (e.g. animals among kingdoms, arthropods and angiosperms among animal and plant phyla). We favor the idea that adaptive radiations have relatively rapid diversification rates relative to other clades, and not simply declining rates over time within clades. The pattern of declining rates could have many explanations that are unrelated to adaptive radiation (Moen and Morlon, 2014).

      Conclusions

      In summary, we show here that most of life’s known, extant species richness belongs to relatively rapid radiations. Furthermore, many of these rapid radiations seem to be related to specific phenotypic traits, and thus may qualify as adaptive radiations (under some definitions). Overall, these results support the emphasis in evolutionary biology on finding the drivers of adaptive radiations and other rapid radiations. These results also suggest that much of life’s species diversity might be the products of a nested series of rapid radiations, with each radiation potentially triggered by a different trait (e.g. first multicellularity, then terrestriality, and then herbivory [within insects]). We speculate that such a pattern might help explain the paradox of how rapid radiation is sustained over long timescales (Martin and Richards, 2019). An exciting area for future research will be to see how general these nested radiations are across scales, and to understand how and why certain traits seem to trigger them.

      Data availability statement

      The original contributions presented in the study are publicly available. This data can be found here: FigShare, https://doi.org/10.6084/m9.figshare.25954669.v1.

      Author contributions

      JW: Investigation, Formal Analysis, Writing – review & editing, Writing – original draft, Data curation, Methodology, Visualization, Conceptualization, Supervision. DM: Writing – review & editing, Conceptualization.

      Funding

      The author(s) declare that no financial support was received for the research and/or publication of this article.

      Conflict of interest

      The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

      Generative AI statement

      The author(s) declare that no Generative AI was used in the creation of this manuscript.

      Publisher’s note

      All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

      Supplementary material

      The Supplementary Material for this article can be found online at: /articles/10.3389/fevo.2025.1596591/full#supplementary-material

      References Baker J. Meade A. Pagel M. Venditti C. (2021). Nothing wrong with the analysis of clades in comparative evolutionary studies: a reply to Poe et al. Syst. Biol. 70, 197201. doi: 10.1093/sysbio/syaa067, PMID: 32845334 Bánki O. Roskov Y. Doring M. Ower G. Hernández Robles D. R. Plata Corredor C. A. . (2024). Catalogue of life checklist (Version 2021-08-25) (Catalogue of Life). Available online at: http://www.catalogueoflife.org. Benton M. J. (2001). Biodiversity on land and in the sea. Geol. J. 36, 211230. doi: 10.1002/gj.877 Bobay L.-M. Ochman H. (2017). Biological species are universal across life’s domains. Genome Biol. Evol. 9, 491501. doi: 10.1093/gbe/evx026, PMID: 28186559 Carroll S. B. (2001). Chance and necessity: the evolution of morphological complexity and diversity. Nature 409, 11021109. doi: 10.1038/35059227, PMID: 11234024 Chen L. Wiens J. J. (2021). Multicellularity and sex helped shape the Tree of Life. Proc. R. Soc. Lond. B 288, 20211265. doi: 10.1098/rspb.2021.1265, PMID: 34315265 Diop A. Torrance E. L. Stott C. M. Bobay L.-M. (2022). Gene flow and introgression are pervasive forces shaping the evolution of bacterial species. Genome Biol. 23, 239. doi: 10.1186/s13059-022-02809-5, PMID: 36357919 Ehrlich P. R. Raven P. H. (1964). Butterflies and plants: a study in coevolution. Evolution 18, 586608. doi: 10.2307/2406212 Farrell B. D. (1998). ‘Inordinate fondness’ explained: why are there so many beetles? Science 281, 555559. doi: 10.1126/science.281.5376.555, PMID: 9677197 Fiz-Palacios O. Schneider H. Heinrichs J. Savolainen V. (2011). Diversification of land plants: insights from a family-level phylogenetic analysis. BMC Evol. Biol. 11, 341. doi: 10.1186/1471-2148-11-341, PMID: 22103931 Forbes A. A. Devine S. N. Hippee A. C. Tvedte E. S. Ward A. K. Widmayer H. A. . (2017). Revisiting the particular role of host shifts in initiating insect speciation. Evolution 71, 11261137. doi: 10.1111/evo.13164, PMID: 28052326 Futuyma D. J. Agrawal A. A. (2009). Macroevolution and the biological diversity of plants and herbivores. Proc. Natl. Acad. Sci. U.S.A. 106, 1805418061. doi: 10.1073/pnas.0904106106, PMID: 19815508 Gillespie R. G. Bennett G. M. De Meester L. Feder J. L. Fleischer R. C. Harmon L. J. . (2020). Comparing adaptive radiations across space, time, and taxa. J. Hered. 111, 120. doi: 10.1093/jhered/esz064, PMID: 31958131 Givnish T. J. (2015). Adaptive radiation versus ‘radiation’ and ‘explosive diversification’: why conceptual distinctions are fundamental to understanding evolution. New Phytol. 207, 297303. doi: 10.1111/nph.13482, PMID: 26032979 Glor R. E. (2010). Phylogenetic insights on adaptive radiation. Annu. Rev. Ecol. Evol. Syst. 41, 251270. doi: 10.1146/annurev.ecolsys.39.110707.173447 Grant V. (1949). Pollination systems as isolating mechanisms in angiosperms. Evolution 3, 8297. doi: 10.1111/j.1558-5646.1949.tb00007.x, PMID: 18115119 Harder L. D. Barrett S. C. H. (2006). Ecology and evolution of flowers (Oxford: Oxford University Press). Hardy N. Otto S. P. (2014). Specialization and generalization in the diversification of phytophagous insects: tests of the musical chairs and oscillation hypotheses. Proc. R. Soc. Lond. 281, 20132960. doi: 10.1098/rspb.2013.2960, PMID: 25274368 Hedges S. B. Marin J. Suleski M. Paymer M. Kumar S. (2015). Tree of Life reveals clock-like speciation and diversification. Mol. Biol. Evol. 32, 835845. doi: 10.1093/molbev/msv037, PMID: 25739733 Henao-Diaz L. F. Harmon L. J. Sugawara M. T. C. Miller E. T. Pennell M. W. (2019). Macroevolutionary diversification rates show time dependency. Proc. Natl. Acad. Sci. U.S.A. 116, 74037408. doi: 10.1073/pnas.1818058116, PMID: 30910958 Hernández-Hernández T. (2019). Evolutionary rates and adaptive radiations. Biol. Philos. 34, 41. doi: 10.1007/s10539-019-9694-y Hernández-Hernández T. Miller E. C. Román-Palacios C. Wiens J. J. (2021). Speciation across the Tree of Life. Biol. Rev. 96, 12051242. doi: 10.1111/brv.12698, PMID: 33768723 Hernández-Hernández T. Wiens J. J. (2020). Why are there so many flowering plants? A multi-scale analysis of plant diversification. Am. Nat. 195, 948963. doi: 10.1086/708273, PMID: 32469653 Irisarri I. Baurain D. Brinkmann H. Delsuc F. Sire J. Y. Kupfer A. . (2017). Phylotranscriptomic consolidation of the jawed vertebrate timetree. Nat. Ecol. Evol. 1, 13701378. doi: 10.1038/s41559-017-0240-5, PMID: 28890940 James T. Y. Kauff F. Schoch C. Matheny P. B. Hofstetter V. Cox C. . (2006). Reconstructing the early evolution of the fungi using a six gene phylogeny. Nature 443, 818822. doi: 10.1038/nature05110, PMID: 17051209 Jezkova T. Wiens J. J. (2017). What explains patterns of diversification and richness among animal phyla? Am. Nat. 189, 201212. doi: 10.1086/690194, PMID: 28221832 Kay K. M. Voelckel C. Yang J. Y. Hufford K. M. Kaska D. D. Hodges S. A. (2006). “Floral characters and species diversification,” in Ecology and evolution of flowers. Eds. Harder L. Barret S. C. H. (Oxford University Press, New York), 311325. Knoll A. H. (1986). “Patterns of change in plant communities through geological time,” in Community ecology. Eds. Diamond J. Case T. J. (Harper and Row, New York), 126141. Kozak K. H. Wiens J. J. (2016). Testing the relationships between diversification, species richness, and trait evolution. Syst. Biol. 65, 975988. doi: 10.1093/sysbio/syw029, PMID: 27048703 Larsen B. B. Miller E. C. Rhodes M. K. Wiens J. J. (2017). Inordinate fondness multiplied and redistributed: the number of species on Earth and the new Pie of Life. Q. Rev. Biol. 92, 229265. doi: 10.1086/693564 Li Y. Steenwyk J. L. Chang Y. Wang Y. James T. Y. Stajich J. E. . (2021). A genome-scale phylogeny of the kingdom Fungi. Curr. Biol. 31, 16531665. doi: 10.1016/j.cub.2021.01.074, PMID: 33607033 Li X. Wiens J. J. (2023). Estimating global biodiversity: the role of cryptic insect species. Syst. Biol. 72, 391403. doi: 10.1093/sysbio/syac069, PMID: 36301156 Locey K. J. Lennon J. T. (2016). Scaling laws predict global microbial diversity. Proc. Natl. Acad. Sci. U.S.A. 113, 59705975. doi: 10.1073/pnas, PMID: 35121666 Louca S. Henao-Diaz L. F. Pennell M. (2022). The scaling of diversification rates with age is likely explained by sampling bias. Evolution 76, 16251637. doi: 10.1111/evo.14515, PMID: 35567800 Louca S. Mazel F. Doebeli M. Parfrey L. W. (2019). A census-based estimate of Earth’s bacterial and archaeal diversity. PLoS Biol. 17, e3000106. doi: 10.1371/journal.pbio.3000106, PMID: 30716065 Louca S. Pennell M. W. (2020). Extant timetrees are consistent with a myriad of diversification histories. Nature 580, 502505. doi: 10.1038/s41586-020-2176-1, PMID: 32322065 Magallón S. Gomez-Acevedo S. Sanchez-Reyes L. L. Hernández-Hernández T. (2015). A metacalibrated time-tree documents the early rise of flowering plant phylogenetic diversity. New Phytol. 207, 437453. doi: 10.1111/nph.13264, PMID: 25615647 Magallón S. Sanderson M. J. (2001). Absolute diversification rates in angiosperm clades. Evolution 55, 17621780. doi: 10.1111/j.0014-3820.2001.tb00826.x, PMID: 11681732 Maliet O. Hartig F. Morlon H. (2019). A model with many small shifts for estimating species-specific diversification rates. Nat. Ecol. Evol. 3, 10861092. doi: 10.1038/s41559-019-0908-0, PMID: 31160736 Maliet O. Morlon H. (2022). Fast and accurate estimation of species-specific diversification rates using data augmentation. Syst. Biol. 71, 353366. doi: 10.1093/sysbio/syab055, PMID: 34228799 Martin C. H. Richards E. J. (2019). The paradox behind the pattern of adaptive radiation: how can the speciation process sustain itself through an early burst? Ann. Rev. Ecol. Evol. Syst. 50, 569594. doi: 10.1146/annurev-ecolsys-110617-062443, PMID: 36237480 May R. M. (1994). Biological diversity: differences between land and sea. Phil. Trans. R. Soc Lond. B 343, 105111. doi: 10.1098/rstb.1994.0014 Mayhew P. J. (2007). Why are there so many insect species? Perspectives from fossils and phylogenies. Biol. Rev. 82, 425454. doi: 10.1111/j.1469-185X.2007.00018.x, PMID: 17624962 McPeek M. A. Brown J. M. (2007). Clade age and not diversification rate explains species richness among animal taxa. Am. Nat. 169, E97106. doi: 10.1086/512135, PMID: 17427118 Meyer A. L. S. Román-Palacios C. Wiens J. J. (2018). BAMM gives misleading rate estimates in simulated and empirical datasets. Evolution 72, 22572266. doi: 10.1111/evo.13574, PMID: 30101971 Meyer A. L. S. Wiens J. J. (2018). Estimating diversification rates for higher taxa: BAMM can give problematic estimates of rates and rate shifts. Evolution 72, 3953. doi: 10.1111/evo.13378, PMID: 29055133 Misof B. S. Liu S. Meusemann K. Peters R. S. Donath A. Mayer C. . (2014). Phylogenomics resolves the timing and pattern of insect evolution. Science 346, 763767. doi: 10.1126/science.1257570, PMID: 25378627 Mitter C. Farrell B. Wiegmann B. (1988). The phylogenetic study of adaptive zones: has phytophagy promoted insect diversification? Am. Nat. 132, 107128. doi: 10.1086/284840 Moen D. S. Morlon H. (2014). Why does diversification slow down? Trends Ecol. Evol. 29, 190197. doi: 10.1016/j.tree.2014.01.010, PMID: 24612774 Moen D. S. Ravelojaona R. N. Hutter C. R. Wiens J. J. (2021). Testing for adaptive radiation: a new approach applied to Madagascar frogs. Evolution 75, 30083025. doi: 10.1111/evo.14328, PMID: 34396527 Mora C. Tittensor D. P. Adl S. Simpson A. G. B. Worm B. (2011). How many species are there on Earth and in the Ocean? PLoS Biol. 9, e1001127. doi: 10.1371/journal.pbio.1001127, PMID: 21886479 Morinaga G. Wiens J. J. Moen D. S. (2023). The radiation continuum and the evolution of frog diversity. Nat. Commun. 14, 7100. doi: 10.1038/s41467-023-42745-x, PMID: 37925440 Niklas K. J. (2016). Plant evolution: an introduction to the history of life (Chicago: University of Chicago Press). Parfrey L. W. Lahr D. J. G. Knoll A. H. Katz L. A. (2011). Estimating the timing of early eukaryotic diversification with multigene molecular clocks. Proc. Natl. Acad. Sci. U.S.A. 108, 1362413629. doi: 10.1073/pnas.1110633108, PMID: 21810989 Poe S. Anderson C. Barnett J. (2021). On the selection and use of clades in comparative evolutionary analyses. Syst. Biol. 70, 190196. doi: 10.1093/sysbio/syaa022, PMID: 32196114 Poe S. Nieto-Montes de Oca A. Torres-Carvajal O. de Queiroz K. Velasco J. A. Truett B. . (2018). Comparative evolution of an archetypal adaptive radiation: innovation and opportunity in Anolis lizards. Am. Nat. 191, E185E194. doi: 10.1086/697223, PMID: 29750558 Rabosky D. L. Benson R. B. J. (2021). Ecological and biogeographic drivers of biodiversity cannot be resolved using clade age-richness data. Nat. Commun. 12, 2945. doi: 10.1038/s41467-021-23307-5, PMID: 34011982 Rabosky D. L. Slater G. J. Alfaro M. E. (2012). Clade age and species richness are decoupled across the eukaryotic tree of life. PLoS Biol. 10, e1001381. doi: 10.1371/journal.pbio.1001381, PMID: 22969411 Rainford J. L. Hofreiter M. Nicholson D. B. Mayhew P. J. (2014). Phylogenetic distribution of extant richness suggests metamorphosis is a key innovation driving diversification in insects. PloS One 9, e109085. doi: 10.1371/journal.pone.0109085, PMID: 25275450 Raven P. H. (1977). A suggestion concerning the Cretaceous rise to dominance of the angiosperms. Evolution 31, 451452. doi: 10.1111/j.1558-5646.1977.tb01029.x, PMID: 28563229 Sanderson M. J. (1996). How many taxa must be sampled to identify the root node of a large clade? Syst. Biol. 45, 168173. doi: 10.1093/sysbio/45.2.168 Sargent R. D. (2004). Floral symmetry affects speciation rates in angiosperms. Proc. R. Soc Lond. B. 271, 603608. doi: 10.1098/rspb.2003.2644, PMID: 15156918 Sauquet H. Ramírez-Barahona S. Magallón S. (2022). What is the age of flowering plants? J. Exp. Bot. 73, 38403853. doi: 10.1093/jxb/erac130, PMID: 35438718 Schluter D. (2000). The ecology of adaptive radiation (Oxford, UK: Oxford University Press). Scholl J. P. Wiens J. J. (2016). Diversification rates and species richness across the Tree of Life. Proc. R. Soc Lond. B 283, 20161335. doi: 10.1098/rspb.2016.1334, PMID: 27605507 Simpson G. G. (1953). The major features of evolution (New York, NY: Columbia Univ. Press). Stebbins G. L. (1981). Why are there so many species of flowering plants? BioScience 31, 573577. doi: 10.2307/1308218 Torino P. Soto M. Perea D. (2021). A comprehensive phylogenetic analysis of coelacanth fishes (Sarcopterygii, Actinistia) with comments on the composition of the Mawsoniidae and Latimeriidae: evaluating old and new methodological challenges and constraints. Hist. Biol. 33, 34233443. doi: 10.1080/08912963.2020.1867982 Van der Niet T. Peakall R. Johnson S. D. (2014). Pollinator driven ecological speciation in plants: new evidence and future perspectives. Ann. Bot. 113, 199212. doi: 10.1093/aob/mct290, PMID: 24418954 Vermeij G. J. Grosberg R. K. (2010). The great divergence: when did diversity on land exceed that in the sea? Int. Comp. Biol. 50, 675682. doi: 10.1093/icb/icq078, PMID: 21558232 Wiens J. J. (2015a). Faster diversification on land than sea helps explain global biodiversity patterns among habitats and animal phyla. Ecol. Lett. 18, 12341241. doi: 10.1111/ele.12503, PMID: 26346782 Wiens J. J. (2015b). Explaining large-scale patterns of vertebrate diversity. Biol. Lett. 11, 20150506. doi: 10.1098/rsbl.2015.0506, PMID: 26202428 Wiens J. J. (2023). The number of species on Earth: progress and problems. PLoS Biol. 21, e300238. doi: 10.1371/journal.pbio.3002388, PMID: 37983223 Wiens J. J. (2024). Speciation across life and the origins of biodiversity patterns. Evol. J. Linn. Soc 3, kzae025. doi: 10.1093/evolinnean/kzae025 Wiens J. J. Lapoint R. T. Whiteman N. K. (2015). Herbivory increases diversification across insect clades. Nat. Commun. 6, 8370. doi: 10.1038/ncomms9370, PMID: 26399434 Williams T. A. Foster P. G. Cox C. J. Embley T. M. (2013). An archaeal origin of eukaryotes supports only two primary domains of life. Nature 504, 231236. doi: 10.1038/nature12779, PMID: 24336283 Yu D. Wiens J. J. (2024). The causes of species richness patterns among clades. Proc. R. Soc Lond. B 291, 20232436. doi: 10.1098/rspb.2023.2436, PMID: 38262607 Zaremba-Niedzwiedzka K. Caceres E. F. Saw J. H. Backstrom D. Juzokaite L. Vancaester E. . (2017). Asgard archaea illuminate the origin of eukaryotic cellular complexity. Nature 541, 353358. doi: 10.1038/nature21031, PMID: 28077874 Zuntini A. R. Carruthers T. Maurin O. Bailey P. C. Leempoel K. Brewer G. G. . (2024). Phylogenomics and the rise of the angiosperms. Nature 629, 843850. doi: 10.1038/s41586-024-07324-0, PMID: 38658746
      ‘Oh, my dear Thomas, you haven’t heard the terrible news then?’ she said. ‘I thought you would be sure to have seen it placarded somewhere. Alice went straight to her room, and I haven’t seen her since, though I repeatedly knocked at the door, which she has locked on the inside, and I’m sure it’s most unnatural of her not to let her own mother comfort her. It all happened in a moment: I have always said those great motor-cars shouldn’t be allowed to career about the streets, especially when they are all paved with cobbles as they are at Easton Haven, which are{331} so slippery when it’s wet. He slipped, and it went over him in a moment.’ My thanks were few and awkward, for there still hung to the missive a basting thread, and it was as warm as a nestling bird. I bent low--everybody was emotional in those days--kissed the fragrant thing, thrust it into my bosom, and blushed worse than Camille. "What, the Corner House victim? Is that really a fact?" "My dear child, I don't look upon it in that light at all. The child gave our picturesque friend a certain distinction--'My husband is dead, and this is my only child,' and all that sort of thing. It pays in society." leave them on the steps of a foundling asylum in order to insure [See larger version] Interoffice guff says you're planning definite moves on your own, J. O., and against some opposition. Is the Colonel so poor or so grasping—or what? Albert could not speak, for he felt as if his brains and teeth were rattling about inside his head. The rest of[Pg 188] the family hunched together by the door, the boys gaping idiotically, the girls in tears. "Now you're married." The host was called in, and unlocked a drawer in which they were deposited. The galleyman, with visible reluctance, arrayed himself in the garments, and he was observed to shudder more than once during the investiture of the dead man's apparel. HoME香京julia种子在线播放 ENTER NUMBET 0016jieohm.com.cn
      hunlizhe.com.cn
      www.hnzz666.com.cn
      www.jyyczz.org.cn
      jynknp.com.cn
      wbit.org.cn
      www.qy0.net.cn
      qxkpoo.com.cn
      www.stackcube.com.cn
      wypiano.com.cn
      处女被大鸡巴操 强奸乱伦小说图片 俄罗斯美女爱爱图 调教强奸学生 亚洲女的穴 夜来香图片大全 美女性强奸电影 手机版色中阁 男性人体艺术素描图 16p成人 欧美性爱360 电影区 亚洲电影 欧美电影 经典三级 偷拍自拍 动漫电影 乱伦电影 变态另类 全部电 类似狠狠鲁的网站 黑吊操白逼图片 韩国黄片种子下载 操逼逼逼逼逼 人妻 小说 p 偷拍10幼女自慰 极品淫水很多 黄色做i爱 日本女人人体电影快播看 大福国小 我爱肏屄美女 mmcrwcom 欧美多人性交图片 肥臀乱伦老头舔阴帝 d09a4343000019c5 西欧人体艺术b xxoo激情短片 未成年人的 插泰国人夭图片 第770弾み1 24p 日本美女性 交动态 eee色播 yantasythunder 操无毛少女屄 亚洲图片你懂的女人 鸡巴插姨娘 特级黄 色大片播 左耳影音先锋 冢本友希全集 日本人体艺术绿色 我爱被舔逼 内射 幼 美阴图 喷水妹子高潮迭起 和后妈 操逼 美女吞鸡巴 鸭个自慰 中国女裸名单 操逼肥臀出水换妻 色站裸体义术 中国行上的漏毛美女叫什么 亚洲妹性交图 欧美美女人裸体人艺照 成人色妹妹直播 WWW_JXCT_COM r日本女人性淫乱 大胆人艺体艺图片 女同接吻av 碰碰哥免费自拍打炮 艳舞写真duppid1 88电影街拍视频 日本自拍做爱qvod 实拍美女性爱组图 少女高清av 浙江真实乱伦迅雷 台湾luanlunxiaoshuo 洛克王国宠物排行榜 皇瑟电影yy频道大全 红孩儿连连看 阴毛摄影 大胆美女写真人体艺术摄影 和风骚三个媳妇在家做爱 性爱办公室高清 18p2p木耳 大波撸影音 大鸡巴插嫩穴小说 一剧不超两个黑人 阿姨诱惑我快播 幼香阁千叶县小学生 少女妇女被狗强奸 曰人体妹妹 十二岁性感幼女 超级乱伦qvod 97爱蜜桃ccc336 日本淫妇阴液 av海量资源999 凤凰影视成仁 辰溪四中艳照门照片 先锋模特裸体展示影片 成人片免费看 自拍百度云 肥白老妇女 女爱人体图片 妈妈一女穴 星野美夏 日本少女dachidu 妹子私处人体图片 yinmindahuitang 舔无毛逼影片快播 田莹疑的裸体照片 三级电影影音先锋02222 妻子被外国老头操 观月雏乃泥鳅 韩国成人偷拍自拍图片 强奸5一9岁幼女小说 汤姆影院av图片 妹妹人艺体图 美女大驱 和女友做爱图片自拍p 绫川まどか在线先锋 那么嫩的逼很少见了 小女孩做爱 处女好逼连连看图图 性感美女在家做爱 近距离抽插骚逼逼 黑屌肏金毛屄 日韩av美少女 看喝尿尿小姐日逼色色色网图片 欧美肛交新视频 美女吃逼逼 av30线上免费 伊人在线三级经典 新视觉影院t6090影院 最新淫色电影网址 天龙影院远古手机版 搞老太影院 插进美女的大屁股里 私人影院加盟费用 www258dd 求一部电影里面有一个二猛哥 深肛交 日本萌妹子人体艺术写真图片 插入屄眼 美女的木奶 中文字幕黄色网址影视先锋 九号女神裸 和骚人妻偷情 和潘晓婷做爱 国模大尺度蜜桃 欧美大逼50p 西西人体成人 李宗瑞继母做爱原图物处理 nianhuawang 男鸡巴的视屏 � 97免费色伦电影 好色网成人 大姨子先锋 淫荡巨乳美女教师妈妈 性nuexiaoshuo WWW36YYYCOM 长春继续给力进屋就操小女儿套干破内射对白淫荡 农夫激情社区 日韩无码bt 欧美美女手掰嫩穴图片 日本援交偷拍自拍 入侵者日本在线播放 亚洲白虎偷拍自拍 常州高见泽日屄 寂寞少妇自卫视频 人体露逼图片 多毛外国老太 变态乱轮手机在线 淫荡妈妈和儿子操逼 伦理片大奶少女 看片神器最新登入地址sqvheqi345com账号群 麻美学姐无头 圣诞老人射小妞和强奸小妞动话片 亚洲AV女老师 先锋影音欧美成人资源 33344iucoom zV天堂电影网 宾馆美女打炮视频 色五月丁香五月magnet 嫂子淫乱小说 张歆艺的老公 吃奶男人视频在线播放 欧美色图男女乱伦 avtt2014ccvom 性插色欲香影院 青青草撸死你青青草 99热久久第一时间 激情套图卡通动漫 幼女裸聊做爱口交 日本女人被强奸乱伦 草榴社区快播 2kkk正在播放兽骑 啊不要人家小穴都湿了 www猎奇影视 A片www245vvcomwwwchnrwhmhzcn 搜索宜春院av wwwsee78co 逼奶鸡巴插 好吊日AV在线视频19gancom 熟女伦乱图片小说 日本免费av无码片在线开苞 鲁大妈撸到爆 裸聊官网 德国熟女xxx 新不夜城论坛首页手机 女虐男网址 男女做爱视频华为网盘 激情午夜天亚洲色图 内裤哥mangent 吉沢明歩制服丝袜WWWHHH710COM 屌逼在线试看 人体艺体阿娇艳照 推荐一个可以免费看片的网站如果被QQ拦截请复制链接在其它浏览器打开xxxyyy5comintr2a2cb551573a2b2e 欧美360精品粉红鲍鱼 教师调教第一页 聚美屋精品图 中韩淫乱群交 俄罗斯撸撸片 把鸡巴插进小姨子的阴道 干干AV成人网 aolasoohpnbcn www84ytom 高清大量潮喷www27dyycom 宝贝开心成人 freefronvideos人母 嫩穴成人网gggg29com 逼着舅妈给我口交肛交彩漫画 欧美色色aV88wwwgangguanscom 老太太操逼自拍视频 777亚洲手机在线播放 有没有夫妻3p小说 色列漫画淫女 午间色站导航 欧美成人处女色大图 童颜巨乳亚洲综合 桃色性欲草 色眯眯射逼 无码中文字幕塞外青楼这是一个 狂日美女老师人妻 爱碰网官网 亚洲图片雅蠛蝶 快播35怎么搜片 2000XXXX电影 新谷露性家庭影院 深深候dvd播放 幼齿用英语怎么说 不雅伦理无需播放器 国外淫荡图片 国外网站幼幼嫩网址 成年人就去色色视频快播 我鲁日日鲁老老老我爱 caoshaonvbi 人体艺术avav 性感性色导航 韩国黄色哥来嫖网站 成人网站美逼 淫荡熟妇自拍 欧美色惰图片 北京空姐透明照 狼堡免费av视频 www776eom 亚洲无码av欧美天堂网男人天堂 欧美激情爆操 a片kk266co 色尼姑成人极速在线视频 国语家庭系列 蒋雯雯 越南伦理 色CC伦理影院手机版 99jbbcom 大鸡巴舅妈 国产偷拍自拍淫荡对话视频 少妇春梦射精 开心激动网 自拍偷牌成人 色桃隐 撸狗网性交视频 淫荡的三位老师 伦理电影wwwqiuxia6commqiuxia6com 怡春院分站 丝袜超短裙露脸迅雷下载 色制服电影院 97超碰好吊色男人 yy6080理论在线宅男日韩福利大全 大嫂丝袜 500人群交手机在线 5sav 偷拍熟女吧 口述我和妹妹的欲望 50p电脑版 wwwavtttcon 3p3com 伦理无码片在线看 欧美成人电影图片岛国性爱伦理电影 先锋影音AV成人欧美 我爱好色 淫电影网 WWW19MMCOM 玛丽罗斯3d同人动画h在线看 动漫女孩裸体 超级丝袜美腿乱伦 1919gogo欣赏 大色逼淫色 www就是撸 激情文学网好骚 A级黄片免费 xedd5com 国内的b是黑的 快播美国成年人片黄 av高跟丝袜视频 上原保奈美巨乳女教师在线观看 校园春色都市激情fefegancom 偷窥自拍XXOO 搜索看马操美女 人本女优视频 日日吧淫淫 人妻巨乳影院 美国女子性爱学校 大肥屁股重口味 啪啪啪啊啊啊不要 操碰 japanfreevideoshome国产 亚州淫荡老熟女人体 伦奸毛片免费在线看 天天影视se 樱桃做爱视频 亚卅av在线视频 x奸小说下载 亚洲色图图片在线 217av天堂网 东方在线撸撸-百度 幼幼丝袜集 灰姑娘的姐姐 青青草在线视频观看对华 86papa路con 亚洲1AV 综合图片2区亚洲 美国美女大逼电影 010插插av成人网站 www色comwww821kxwcom 播乐子成人网免费视频在线观看 大炮撸在线影院 ,www4KkKcom 野花鲁最近30部 wwwCC213wapwww2233ww2download 三客优最新地址 母亲让儿子爽的无码视频 全国黄色片子 欧美色图美国十次 超碰在线直播 性感妖娆操 亚洲肉感熟女色图 a片A毛片管看视频 8vaa褋芯屑 333kk 川岛和津实视频 在线母子乱伦对白 妹妹肥逼五月 亚洲美女自拍 老婆在我面前小说 韩国空姐堪比情趣内衣 干小姐综合 淫妻色五月 添骚穴 WM62COM 23456影视播放器 成人午夜剧场 尼姑福利网 AV区亚洲AV欧美AV512qucomwwwc5508com 经典欧美骚妇 震动棒露出 日韩丝袜美臀巨乳在线 av无限吧看 就去干少妇 色艺无间正面是哪集 校园春色我和老师做爱 漫画夜色 天海丽白色吊带 黄色淫荡性虐小说 午夜高清播放器 文20岁女性荫道口图片 热国产热无码热有码 2015小明发布看看算你色 百度云播影视 美女肏屄屄乱轮小说 家族舔阴AV影片 邪恶在线av有码 父女之交 关于处女破处的三级片 极品护士91在线 欧美虐待女人视频的网站 享受老太太的丝袜 aaazhibuo 8dfvodcom成人 真实自拍足交 群交男女猛插逼 妓女爱爱动态 lin35com是什么网站 abp159 亚洲色图偷拍自拍乱伦熟女抠逼自慰 朝国三级篇 淫三国幻想 免费的av小电影网站 日本阿v视频免费按摩师 av750c0m 黄色片操一下 巨乳少女车震在线观看 操逼 免费 囗述情感一乱伦岳母和女婿 WWW_FAMITSU_COM 偷拍中国少妇在公车被操视频 花也真衣论理电影 大鸡鸡插p洞 新片欧美十八岁美少 进击的巨人神thunderftp 西方美女15p 深圳哪里易找到老女人玩视频 在线成人有声小说 365rrr 女尿图片 我和淫荡的小姨做爱 � 做爱技术体照 淫妇性爱 大学生私拍b 第四射狠狠射小说 色中色成人av社区 和小姨子乱伦肛交 wwwppp62com 俄罗斯巨乳人体艺术 骚逼阿娇 汤芳人体图片大胆 大胆人体艺术bb私处 性感大胸骚货 哪个网站幼女的片多 日本美女本子把 色 五月天 婷婷 快播 美女 美穴艺术 色百合电影导航 大鸡巴用力 孙悟空操美少女战士 狠狠撸美女手掰穴图片 古代女子与兽类交 沙耶香套图 激情成人网区 暴风影音av播放 动漫女孩怎么插第3个 mmmpp44 黑木麻衣无码ed2k 淫荡学姐少妇 乱伦操少女屄 高中性爱故事 骚妹妹爱爱图网 韩国模特剪长发 大鸡巴把我逼日了 中国张柏芝做爱片中国张柏芝做爱片中国张柏芝做爱片中国张柏芝做爱片中国张柏芝做爱片 大胆女人下体艺术图片 789sss 影音先锋在线国内情侣野外性事自拍普通话对白 群撸图库 闪现君打阿乐 ady 小说 插入表妹嫩穴小说 推荐成人资源 网络播放器 成人台 149大胆人体艺术 大屌图片 骚美女成人av 春暖花开春色性吧 女亭婷五月 我上了同桌的姐姐 恋夜秀场主播自慰视频 yzppp 屄茎 操屄女图 美女鲍鱼大特写 淫乱的日本人妻山口玲子 偷拍射精图 性感美女人体艺木图片 种马小说完本 免费电影院 骑士福利导航导航网站 骚老婆足交 国产性爱一级电影 欧美免费成人花花性都 欧美大肥妞性爱视频 家庭乱伦网站快播 偷拍自拍国产毛片 金发美女也用大吊来开包 缔D杏那 yentiyishu人体艺术ytys WWWUUKKMCOM 女人露奶 � 苍井空露逼 老荡妇高跟丝袜足交 偷偷和女友的朋友做爱迅雷 做爱七十二尺 朱丹人体合成 麻腾由纪妃 帅哥撸播种子图 鸡巴插逼动态图片 羙国十次啦中文 WWW137AVCOM 神斗片欧美版华语 有气质女人人休艺术 由美老师放屁电影 欧美女人肉肏图片 白虎种子快播 国产自拍90后女孩 美女在床上疯狂嫩b 饭岛爱最后之作 幼幼强奸摸奶 色97成人动漫 两性性爱打鸡巴插逼 新视觉影院4080青苹果影院 嗯好爽插死我了 阴口艺术照 李宗瑞电影qvod38 爆操舅母 亚洲色图七七影院 被大鸡巴操菊花 怡红院肿么了 成人极品影院删除 欧美性爱大图色图强奸乱 欧美女子与狗随便性交 苍井空的bt种子无码 熟女乱伦长篇小说 大色虫 兽交幼女影音先锋播放 44aad be0ca93900121f9b 先锋天耗ばさ无码 欧毛毛女三级黄色片图 干女人黑木耳照 日本美女少妇嫩逼人体艺术 sesechangchang 色屄屄网 久久撸app下载 色图色噜 美女鸡巴大奶 好吊日在线视频在线观看 透明丝袜脚偷拍自拍 中山怡红院菜单 wcwwwcom下载 骑嫂子 亚洲大色妣 成人故事365ahnet 丝袜家庭教mp4 幼交肛交 妹妹撸撸大妈 日本毛爽 caoprom超碰在email 关于中国古代偷窥的黄片 第一会所老熟女下载 wwwhuangsecome 狼人干综合新地址HD播放 变态儿子强奸乱伦图 强奸电影名字 2wwwer37com 日本毛片基地一亚洲AVmzddcxcn 暗黑圣经仙桃影院 37tpcocn 持月真由xfplay 好吊日在线视频三级网 我爱背入李丽珍 电影师傅床戏在线观看 96插妹妹sexsex88com 豪放家庭在线播放 桃花宝典极夜著豆瓜网 安卓系统播放神器 美美网丝袜诱惑 人人干全免费视频xulawyercn av无插件一本道 全国色五月 操逼电影小说网 good在线wwwyuyuelvcom www18avmmd 撸波波影视无插件 伊人幼女成人电影 会看射的图片 小明插看看 全裸美女扒开粉嫩b 国人自拍性交网站 萝莉白丝足交本子 七草ちとせ巨乳视频 摇摇晃晃的成人电影 兰桂坊成社人区小说www68kqcom 舔阴论坛 久撸客一撸客色国内外成人激情在线 明星门 欧美大胆嫩肉穴爽大片 www牛逼插 性吧星云 少妇性奴的屁眼 人体艺术大胆mscbaidu1imgcn 最新久久色色成人版 l女同在线 小泽玛利亚高潮图片搜索 女性裸b图 肛交bt种子 最热门有声小说 人间添春色 春色猜谜字 樱井莉亚钢管舞视频 小泽玛利亚直美6p 能用的h网 还能看的h网 bl动漫h网 开心五月激 东京热401 男色女色第四色酒色网 怎么下载黄色小说 黄色小说小栽 和谐图城 乐乐影院 色哥导航 特色导航 依依社区 爱窝窝在线 色狼谷成人 91porn 包要你射电影 色色3A丝袜 丝袜妹妹淫网 爱色导航(荐) 好男人激情影院 坏哥哥 第七色 色久久 人格分裂 急先锋 撸撸射中文网 第一会所综合社区 91影院老师机 东方成人激情 怼莪影院吹潮 老鸭窝伊人无码不卡无码一本道 av女柳晶电影 91天生爱风流作品 深爱激情小说私房婷婷网 擼奶av 567pao 里番3d一家人野外 上原在线电影 水岛津实透明丝袜 1314酒色 网旧网俺也去 0855影院 在线无码私人影院 搜索 国产自拍 神马dy888午夜伦理达达兔 农民工黄晓婷 日韩裸体黑丝御姐 屈臣氏的燕窝面膜怎么样つぼみ晶エリーの早漏チ○ポ强化合宿 老熟女人性视频 影音先锋 三上悠亚ol 妹妹影院福利片 hhhhhhhhsxo 午夜天堂热的国产 强奸剧场 全裸香蕉视频无码 亚欧伦理视频 秋霞为什么给封了 日本在线视频空天使 日韩成人aⅴ在线 日本日屌日屄导航视频 在线福利视频 日本推油无码av magnet 在线免费视频 樱井梨吮东 日本一本道在线无码DVD 日本性感诱惑美女做爱阴道流水视频 日本一级av 汤姆avtom在线视频 台湾佬中文娱乐线20 阿v播播下载 橙色影院 奴隶少女护士cg视频 汤姆在线影院无码 偷拍宾馆 业面紧急生级访问 色和尚有线 厕所偷拍一族 av女l 公交色狼优酷视频 裸体视频AV 人与兽肉肉网 董美香ol 花井美纱链接 magnet 西瓜影音 亚洲 自拍 日韩女优欧美激情偷拍自拍 亚洲成年人免费视频 荷兰免费成人电影 深喉呕吐XXⅩX 操石榴在线视频 天天色成人免费视频 314hu四虎 涩久免费视频在线观看 成人电影迅雷下载 能看见整个奶子的香蕉影院 水菜丽百度影音 gwaz079百度云 噜死你们资源站 主播走光视频合集迅雷下载 thumbzilla jappen 精品Av 古川伊织star598在线 假面女皇vip在线视频播放 国产自拍迷情校园 啪啪啪公寓漫画 日本阿AV 黄色手机电影 欧美在线Av影院 华裔电击女神91在线 亚洲欧美专区 1日本1000部免费视频 开放90后 波多野结衣 东方 影院av 页面升级紧急访问每天正常更新 4438Xchengeren 老炮色 a k福利电影 色欲影视色天天视频 高老庄aV 259LUXU-683 magnet 手机在线电影 国产区 欧美激情人人操网 国产 偷拍 直播 日韩 国内外激情在线视频网给 站长统计一本道人妻 光棍影院被封 紫竹铃取汁 ftp 狂插空姐嫩 xfplay 丈夫面前 穿靴子伪街 XXOO视频在线免费 大香蕉道久在线播放 电棒漏电嗨过头 充气娃能看下毛和洞吗 夫妻牲交 福利云点墦 yukun瑟妃 疯狂交换女友 国产自拍26页 腐女资源 百度云 日本DVD高清无码视频 偷拍,自拍AV伦理电影 A片小视频福利站。 大奶肥婆自拍偷拍图片 交配伊甸园 超碰在线视频自拍偷拍国产 小热巴91大神 rctd 045 类似于A片 超美大奶大学生美女直播被男友操 男友问 你的衣服怎么脱掉的 亚洲女与黑人群交视频一 在线黄涩 木内美保步兵番号 鸡巴插入欧美美女的b舒服 激情在线国产自拍日韩欧美 国语福利小视频在线观看 作爱小视颍 潮喷合集丝袜无码mp4 做爱的无码高清视频 牛牛精品 伊aⅤ在线观看 savk12 哥哥搞在线播放 在线电一本道影 一级谍片 250pp亚洲情艺中心,88 欧美一本道九色在线一 wwwseavbacom色av吧 cos美女在线 欧美17,18ⅹⅹⅹ视频 自拍嫩逼 小电影在线观看网站 筱田优 贼 水电工 5358x视频 日本69式视频有码 b雪福利导航 韩国女主播19tvclub在线 操逼清晰视频 丝袜美女国产视频网址导航 水菜丽颜射房间 台湾妹中文娱乐网 风吟岛视频 口交 伦理 日本熟妇色五十路免费视频 A级片互舔 川村真矢Av在线观看 亚洲日韩av 色和尚国产自拍 sea8 mp4 aV天堂2018手机在线 免费版国产偷拍a在线播放 狠狠 婷婷 丁香 小视频福利在线观看平台 思妍白衣小仙女被邻居强上 萝莉自拍有水 4484新视觉 永久发布页 977成人影视在线观看 小清新影院在线观 小鸟酱后丝后入百度云 旋风魅影四级 香蕉影院小黄片免费看 性爱直播磁力链接 小骚逼第一色影院 性交流的视频 小雪小视频bd 小视频TV禁看视频 迷奸AV在线看 nba直播 任你在干线 汤姆影院在线视频国产 624u在线播放 成人 一级a做爰片就在线看狐狸视频 小香蕉AV视频 www182、com 腿模简小育 学生做爱视频 秘密搜查官 快播 成人福利网午夜 一级黄色夫妻录像片 直接看的gav久久播放器 国产自拍400首页 sm老爹影院 谁知道隔壁老王网址在线 综合网 123西瓜影音 米奇丁香 人人澡人人漠大学生 色久悠 夜色视频你今天寂寞了吗? 菲菲影视城美国 被抄的影院 变态另类 欧美 成人 国产偷拍自拍在线小说 不用下载安装就能看的吃男人鸡巴视频 插屄视频 大贯杏里播放 wwwhhh50 233若菜奈央 伦理片天海翼秘密搜查官 大香蕉在线万色屋视频 那种漫画小说你懂的 祥仔电影合集一区 那里可以看澳门皇冠酒店a片 色自啪 亚洲aV电影天堂 谷露影院ar toupaizaixian sexbj。com 毕业生 zaixian mianfei 朝桐光视频 成人短视频在线直接观看 陈美霖 沈阳音乐学院 导航女 www26yjjcom 1大尺度视频 开平虐女视频 菅野雪松协和影视在线视频 华人play在线视频bbb 鸡吧操屄视频 多啪啪免费视频 悠草影院 金兰策划网 (969) 橘佑金短视频 国内一极刺激自拍片 日本制服番号大全magnet 成人动漫母系 电脑怎么清理内存 黄色福利1000 dy88午夜 偷拍中学生洗澡磁力链接 花椒相机福利美女视频 站长推荐磁力下载 mp4 三洞轮流插视频 玉兔miki热舞视频 夜生活小视频 爆乳人妖小视频 国内网红主播自拍福利迅雷下载 不用app的裸裸体美女操逼视频 变态SM影片在线观看 草溜影院元气吧 - 百度 - 百度 波推全套视频 国产双飞集合ftp 日本在线AV网 笔国毛片 神马影院女主播是我的邻居 影音资源 激情乱伦电影 799pao 亚洲第一色第一影院 av视频大香蕉 老梁故事汇希斯莱杰 水中人体磁力链接 下载 大香蕉黄片免费看 济南谭崔 避开屏蔽的岛a片 草破福利 要看大鸡巴操小骚逼的人的视频 黑丝少妇影音先锋 欧美巨乳熟女磁力链接 美国黄网站色大全 伦蕉在线久播 极品女厕沟 激情五月bd韩国电影 混血美女自摸和男友激情啪啪自拍诱人呻吟福利视频 人人摸人人妻做人人看 44kknn 娸娸原网 伊人欧美 恋夜影院视频列表安卓青青 57k影院 如果电话亭 avi 插爆骚女精品自拍 青青草在线免费视频1769TV 令人惹火的邻家美眉 影音先锋 真人妹子被捅动态图 男人女人做完爱视频15 表姐合租两人共处一室晚上她竟爬上了我的床 性爱教学视频 北条麻妃bd在线播放版 国产老师和师生 magnet wwwcctv1024 女神自慰 ftp 女同性恋做激情视频 欧美大胆露阴视频 欧美无码影视 好女色在线观看 后入肥臀18p 百度影视屏福利 厕所超碰视频 强奸mp magnet 欧美妹aⅴ免费线上看 2016年妞干网视频 5手机在线福利 超在线最视频 800av:cOm magnet 欧美性爱免播放器在线播放 91大款肥汤的性感美乳90后邻家美眉趴着窗台后入啪啪 秋霞日本毛片网站 cheng ren 在线视频 上原亚衣肛门无码解禁影音先锋 美脚家庭教师在线播放 尤酷伦理片 熟女性生活视频在线观看 欧美av在线播放喷潮 194avav 凤凰AV成人 - 百度 kbb9999 AV片AV在线AV无码 爱爱视频高清免费观看 黄色男女操b视频 观看 18AV清纯视频在线播放平台 成人性爱视频久久操 女性真人生殖系统双性人视频 下身插入b射精视频 明星潜规测视频 mp4 免賛a片直播绪 国内 自己 偷拍 在线 国内真实偷拍 手机在线 国产主播户外勾在线 三桥杏奈高清无码迅雷下载 2五福电影院凸凹频频 男主拿鱼打女主,高宝宝 色哥午夜影院 川村まや痴汉 草溜影院费全过程免费 淫小弟影院在线视频 laohantuiche 啪啪啪喷潮XXOO视频 青娱乐成人国产 蓝沢润 一本道 亚洲青涩中文欧美 神马影院线理论 米娅卡莉法的av 在线福利65535 欧美粉色在线 欧美性受群交视频1在线播放 极品喷奶熟妇在线播放 变态另类无码福利影院92 天津小姐被偷拍 磁力下载 台湾三级电髟全部 丝袜美腿偷拍自拍 偷拍女生性行为图 妻子的乱伦 白虎少妇 肏婶骚屄 外国大妈会阴照片 美少女操屄图片 妹妹自慰11p 操老熟女的b 361美女人体 360电影院樱桃 爱色妹妹亚洲色图 性交卖淫姿势高清图片一级 欧美一黑对二白 大色网无毛一线天 射小妹网站 寂寞穴 西西人体模特苍井空 操的大白逼吧 骚穴让我操 拉好友干女朋友3p