Testing the effect of the rock record on diversity: a multidisciplinary approach to elucidating the generic richness of sauropodomorph dinosaurs through time

The accurate reconstruction of palaeobiodiversity patterns is central to a detailed understanding of the macroevolutionary history of a group of organisms. However, there is increasing evidence that diversity patterns observed directly from the fossil record are strongly influenced by fluctuations in the quality of our sampling of the rock record; thus, any patterns we see may reflect sampling biases, rather than genuine biological signals. Previous dinosaur diversity studies have suggested that fluctuations in sauropodomorph palaeobiodiversity reflect genuine biological signals, in comparison to theropods and ornithischians whose diversity seems to be largely controlled by the rock record. Most previous diversity analyses that have attempted to take into account the effects of sampling biases have used only a single method or proxy: here we use a number of techniques in order to elucidate diversity. A global database of all known sauropodomorph body fossil occurrences (2024) was constructed. A taxic diversity curve for all valid sauropodomorph genera was extracted from this database and compared statistically with several sampling proxies (rock outcrop area and dinosaur‐bearing formations and collections), each of which captures a different aspect of fossil record sampling. Phylogenetic diversity estimates, residuals and sample‐based rarefaction (including the first attempt to capture ‘cryptic’ diversity in dinosaurs) were implemented to investigate further the effects of sampling. After ‘removal’ of biases, sauropodomorph diversity appears to be genuinely high in the Norian, Pliensbachian–Toarcian, Bathonian–Callovian and Kimmeridgian–Tithonian (with a small peak in the Aptian), whereas low diversity levels are recorded for the Oxfordian and Berriasian–Barremian, with the Jurassic/Cretaceous boundary seemingly representing a real diversity trough. Observed diversity in the remaining Triassic–Jurassic stages appears to be largely driven by sampling effort. Late Cretaceous diversity is difficult to elucidate and it is possible that this interval remains relatively under‐sampled. Despite its distortion by sampling biases, much of sauropodomorph palaeobiodiversity can be interpreted as a reflection of genuine biological signals, and fluctuations in sea level may account for some of these diversity patterns.     

Phanerozoic marine diversity: rock record modelling provides an independent test of large-scale trends

Sampling bias created by a heterogeneous rock record can seriously distort estimates of marine diversity and makes a direct reading of the fossil record unreliable. Here we compare two independent estimates of Phanerozoic marine diversity that explicitly take account of variation in sampling—a subsampling approach that standardizes for differences in fossil collection intensity, and a rock area modelling approach that takes account of differences in rock availability. Using the fossil records of North America and Western Europe, we demonstrate that a modelling approach applied to the combined data produces results that are significantly correlated with those derived from subsampling. This concordance between independent approaches argues strongly for the reality of the large-scale trends in diversity we identify from both approaches.     

Mesozoic marine tetrapod diversity: mass extinctions and temporal heterogeneity in geological megabiases affecting vertebrates

The fossil record is our only direct means for evaluating shifts in biodiversity through Earth''s history. However, analyses of fossil marine invertebrates have demonstrated that geological megabiases profoundly influence fossil preservation and discovery, obscuring true diversity signals. Comparable studies of vertebrate palaeodiversity patterns remain in their infancy. A new species-level dataset of Mesozoic marine tetrapod occurrences was compared with a proxy for temporal variation in the volume and facies diversity of fossiliferous rock (number of marine fossiliferous formations: FMF). A strong correlation between taxic diversity and FMF is present during the Cretaceous. Weak or no correlation of Jurassic data suggests a qualitatively different sampling regime resulting from five apparent peaks in Triassic–Jurassic diversity. These correspond to a small number of European formations that have been the subject of intensive collecting, and represent ‘Lagerstätten effects’. Consideration of sampling biases allows re-evaluation of proposed mass extinction events. Marine tetrapod diversity declined during the Carnian or Norian. However, the proposed end-Triassic extinction event cannot be recognized with confidence. Some evidence supports an extinction event near the Jurassic/Cretaceous boundary, but the proposed end-Cenomanian extinction is probably an artefact of poor sampling. Marine tetrapod diversity underwent a long-term decline prior to the Cretaceous–Palaeogene extinction.     

Preservational bias controls the fossil record of pterosaurs

Pterosaurs, a Mesozoic group of flying archosaurs, have become a focal point for debates pertaining to the impact of sampling biases on our reading of the fossil record, as well as the utility of sampling proxies in palaeo‐diversity reconstructions. The completeness of the pterosaur fossil specimens themselves potentially provides additional information that is not captured in existing sampling proxies, and might shed new light on the group's evolutionary history. Here we assess the quality of the pterosaur fossil record via a character completeness metric based on the number of phylogenetic characters that can be scored for all known skeletons of 172 valid species, with averaged completeness values calculated for each geological stage. The fossil record of pterosaurs is observed to be strongly influenced by the occurrence and distribution of Lagerstätten. Peaks in completeness correlate with Lagerstätten deposits, and a recovered correlation between completeness and observed diversity is rendered non‐significant when Lagerstätten species are excluded. Intervals previously regarded as potential extinction events are shown to lack Lagerstätten and exhibit low completeness values: as such, the apparent low diversity in these intervals might be at least partly the result of poor fossil record quality. A positive correlation between temporal patterns in completeness of Cretaceous pterosaurs and birds further demonstrates the prominent role that Lagerstätten deposits have on the preservation of smaller bodied organisms, contrasting with a lack of correlation with the completeness of large‐bodied sauropodomorphs. However, we unexpectedly find a strong correlation between sauropodomorph and pterosaur completeness within the Triassic–Jurassic, but not the Cretaceous, potentially relating to a shared shift in environmental preference and thus preservation style through time. This study highlights the importance of understanding the relationship between various taphonomic controls when correcting for sampling bias, and provides additional evidence for the prominent role of sampling on observed patterns in pterosaur macroevolution.     

Severe extinction and rapid recovery of mammals across the Cretaceous–Palaeogene boundary,and the effects of rarity on patterns of extinction and recovery

The end‐Cretaceous mass extinction ranks among the most severe extinctions of all time; however, patterns of extinction and recovery remain incompletely understood. In particular, it is unclear how severe the extinction was, how rapid the recovery was and how sampling biases might affect our understanding of these processes. To better understand terrestrial extinction and recovery and how sampling influences these patterns, we collected data on the occurrence and abundance of fossil mammals to examine mammalian diversity across the K‐Pg boundary in North America. Our data show that the extinction was more severe and the recovery more rapid than previously thought. Extinction rates are markedly higher than previously estimated: of 59 species, four survived (93% species extinction, 86% of genera). Survival is correlated with geographic range size and abundance, with widespread, common species tending to survive. This creates a sampling artefact in which rare species are both more vulnerable to extinction and less likely to be recovered, such that the fossil record is inherently biased towards the survivors. The recovery was remarkably rapid. Within 300 000 years, local diversity recovered and regional diversity rose to twice Cretaceous levels, driven by increased endemicity; morphological disparity increased above levels observed in the Cretaceous. The speed of the recovery tends to be obscured by sampling effects; faunas show increased endemicity, such that a rapid, regional increase in diversity and disparity is not seen in geographically restricted studies. Sampling biases that operate against rare taxa appear to obscure the severity of extinction and the pace of recovery across the K‐Pg boundary, and similar biases may operate during other extinction events.     

A long-term association between global temperature and biodiversity, origination and extinction in the fossil record

The past relationship between global temperature and levels of biological diversity is of increasing concern due to anthropogenic climate warming. However, no consistent link between these variables has yet been demonstrated. We analysed the fossil record for the last 520 Myr against estimates of low latitude sea surface temperature for the same period. We found that global biodiversity (the richness of families and genera) is related to temperature and has been relatively low during warm 'greenhouse' phases, while during the same phases extinction and origination rates of taxonomic lineages have been relatively high. These findings are consistent for terrestrial and marine environments and are robust to a number of alternative assumptions and potential biases. Our results provide the first clear evidence that global climate may explain substantial variation in the fossil record in a simple and consistent manner. Our findings may have implications for extinction and biodiversity change under future climate warming.     

Spatiotemporal sampling patterns in the 230 million year fossil record of terrestrial crocodylomorphs and their impact on diversity

The 24 extant crocodylian species are the remnants of a once much more diverse and widespread clade. Crocodylomorpha has an approximately 230 million year evolutionary history, punctuated by a series of radiations and extinctions. However, the group's fossil record is biased. Previous studies have reconstructed temporal patterns in subsampled crocodylomorph palaeobiodiversity, but have not explicitly examined variation in spatial sampling, nor the quality of this record. We compiled a dataset of all taxonomically diagnosable non‐marine crocodylomorph species (393). Based on the number of phylogenetic characters that can be scored for all published fossils of each species, we calculated a completeness value for each taxon. Mean average species completeness (56%) is largely consistent within subgroups and for different body size classes, suggesting no significant biases across the crocodylomorph tree. In general, average completeness values are highest in the Mesozoic, with an overall trend of decreasing completeness through time. Many extant taxa are identified in the fossil record from very incomplete remains, but this might be because their provenance closely matches the species’ present‐day distribution, rather than through autapomorphies. Our understanding of nearly all crocodylomorph macroevolutionary ‘events’ is essentially driven by regional patterns, with no global sampling signal. Palaeotropical sampling is especially poor for most of the group's history. Spatiotemporal sampling bias impedes our understanding of several Mesozoic radiations, whereas molecular divergence times for Crocodylia are generally in close agreement with the fossil record. However, the latter might merely be fortuitous, i.e. divergences happened to occur during our ephemeral spatiotemporal sampling windows.     

Skeletal completeness of the non‐avian theropod dinosaur fossil record

Non‐avian theropods were a highly successful clade of bipedal, predominantly carnivorous, dinosaurs. Their diversity and macroevolutionary patterns have been the subject of many studies. Changes in fossil specimen completeness through time and space can bias our understanding of macroevolution. Here, we quantify the completeness of 455 non‐avian theropod species using the skeletal completeness metric (SCM), which calculates the proportion of a complete skeleton preserved for a specimen. Temporal patterns of theropod skeletal completeness show peaks in the Carnian, Oxfordian–Kimmeridgian and Barremian–Aptian, and lows in the Berriasian and Hauterivian. Lagerstätten primarily drive the peaks in completeness and observed taxonomic diversity in the Oxfordian–Kimmeridgian and the Barremian–Aptian. Theropods have a significantly lower distribution of completeness scores than contemporary sauropodomorph dinosaurs but change in completeness through time for the two groups shows a significant correlation when conservation Lagerstätten are excluded, possibly indicating that both records are primarily driven by geology and sampling availability. Our results reveal relatively weak temporal sampling biases acting on the theropod record but relatively strong spatial and environmental biases. Asia has a significantly more complete record than any other continent, the mid northern latitudes have the highest abundance of finds, and most complete theropod skeletons come from lacustrine and aeolian environments. We suggest that these patterns result from historical research focus, modern climate dynamics, and depositional transportation energy plus association with conservation Lagerstätten, respectively. Furthermore, we find possible ecological biases acting on different theropod subgroups, but body size does not influence theropod completeness on a global scale.     

Dinosaur extinction: closing the '3 m gap'

Modern debate regarding the extinction of non-avian dinosaurs was ignited by the publication of the Cretaceous-Tertiary (K-T) asteroid impact theory and has seen 30 years of dispute over the position of the stratigraphically youngest in situ dinosaur. A zone devoid of dinosaur fossils reported from the last 3 m of the Upper Cretaceous, coined the '3 m gap', has helped drive controversy. Here, we report the discovery of the stratigraphically youngest in situ dinosaur specimen: a ceratopsian brow horn found in a poorly rooted, silty, mudstone floodplain deposit located no more than 13 cm below the palynologically defined boundary. The K-T boundary is identified using three criteria: (i) decrease in Cretaceous palynomorphs without subsequent recovery, (ii) the existence of a 'fern spike', and (iii) correlation to a nearby stratigraphic section where primary extraterrestrial impact markers are present (e.g. iridium anomaly, spherules, shocked quartz). The in situ specimen demonstrates that a gap devoid of non-avian dinosaur fossils does not exist and is inconsistent with the hypothesis that non-avian dinosaurs were extinct prior to the K-T boundary impact event.     

Ten more years of discovery: revisiting the quality of the sauropodomorph dinosaur fossil record

Spatiotemporal changes in fossil specimen completeness can bias our understanding of a group's evolutionary history. The quality of the sauropodomorph fossil record was assessed a decade ago, but the number of valid species has since increased by 60%, and 17% of the taxa from that study have since undergone taxonomic revision. Here, we assess how 10 years of additional research has changed our outlook on the group's fossil record. We quantified the completeness of all 307 sauropodomorph species currently considered valid using the skeletal completeness metric, which calculates the proportion of a complete skeleton preserved for each taxon. Taxonomic and stratigraphic age revisions, rather than new species, are the drivers of the most significant differences between the current results and those of the previous assessment. No statistical differences appeared when we use our new dataset to generate temporal completeness curves based solely on taxa known in 2009 or 1999. We now observe a severe drop in mean completeness values across the Jurassic–Cretaceous boundary that never recovers to pre-Cretaceous levels. Explaining this pattern is difficult, as we find no convincing evidence that it is related to environmental preferences or body size changes. Instead, it might result from: (1) reduction of terrestrial fossil preservation space due to sea level rise; (2) ecological specificities and relatively high diagnosability of Cretaceous species; and/or (3) increased sampling of newly explored sites with many previously unknown taxa. Revisiting patterns in this manner allows us to test the longevity of conclusions made in previous quantitative studies.     

A refined modelling approach to assess the influence of sampling on palaeobiodiversity curves: new support for declining Cretaceous dinosaur richness

Modelling has been underdeveloped with respect to constructing palaeobiodiversity curves, but it offers an additional tool for removing sampling from their estimation. Here, an alternative to subsampling approaches, which often require large sample sizes, is explored by the extension and refinement of a pre-existing modelling technique that uses a geological proxy for sampling. Application of the model to the three main clades of dinosaurs suggests that much of their diversity fluctuations cannot be explained by sampling alone. Furthermore, there is new support for a long-term decline in their diversity leading up to the Cretaceous–Paleogene (K–Pg) extinction event. At present, use of this method with data that includes either Lagerstätten or ‘Pull of the Recent’ biases is inappropriate, although partial solutions are offered.     

The fossil record of ichthyosaurs,completeness metrics and sampling biases

Ichthyosaurs were highly successful marine reptiles with an abundant and well‐studied fossil record. However, their occurrences through geological time and space are sporadic, and it is important to understand whether times of apparent species richness and rarity are real or the result of sampling bias. Here, we explore the skeletal completeness of 351 dated and identified ichthyosaur specimens, belonging to all 102 species, the first time that such a study has been carried out on vertebrates from the marine realm. No correlations were found between time series of different skeletal metrics and ichthyosaur diversity. There is a significant geographical variation in completeness, with the well‐studied northern hemisphere producing fossils of much higher quality than the southern hemisphere. Medium‐sized ichthyosaurs are significantly more complete than small or large taxa: the incompleteness of small specimens was expected, but it was a surprise that larger specimens were also relatively incomplete. Completeness varies greatly between facies, with fine‐grained, siliciclastic sediments preserving the most complete specimens. These findings may explain why the ichthyosaur diversity record is low at times, corresponding to facies of poor preservation potential, such as in the Early Cretaceous. Unexpectedly, we find a strong negative correlation between skeletal completeness and sea level, meaning the most complete specimens occurred at times of global low sea level, and vice versa. Completeness metrics, however, do not replicate the sampling signal and have limited use as a global‐scale sampling proxy.     

Disparity vs. diversity in Stegosauria (Dinosauria,Ornithischia): cranial and post-cranial sub-dataset provide different signals

The Stegosauria represents an iconic group of ornithischian dinosaurs, with a fossil record spanning the Middle Jurassic to the Late Cretaceous. In this contribution I present the first detailed analysis of the relationship between disparity and diversity through the evolutionary history of the group. The analysis has been performed on a recently published cladistic dataset, allowing the separate study of the signals deriving from discrete characters and from continuous morphometric characters. Whereas the disparity as sum of variance is decoupled with respect to diversity, the sum of ranges provides a signal fairly consistent with the trend in the number of taxa. Both sub-data sets show that evolution of stegosaurs can be considered essentially as symmetrical, i.e. the maximum exploration of the possible morphospace takes place about half way through the history of the group, with subsequent significant decline until extinction in the Upper Cretaceous. An interesting result is a decoupling of the tempo and mode of evolution of the cranium and postcranium in stegosaurs. Specifically, the evolutionary radiation with maximum saturation of morphospace is anticipated in the cranial skeleton, with maximum peak in the Oxfordian; in contrast, the postcranium explores the largest number of morphotypes subsequently during the Kimmeridgian.     

Taxonomic identification bias does not drive patterns of abundance and diversity in theropod dinosaurs

The ability of palaeontologists to correctly diagnose and classify new fossil species from incomplete morphological data is fundamental to our understanding of evolution. Different parts of the vertebrate skeleton have different likelihoods of fossil preservation and varying amounts of taxonomic information, which could bias our interpretations of fossil material. Substantial previous research has focused on the diversity and macroevolution of non-avian theropod dinosaurs. Theropods provide a rich dataset for analysis of the interactions between taxonomic diagnosability and fossil preservation. We use specimen data and formal taxonomic diagnoses to create a new metric, the Likelihood of Diagnosis, which quantifies the diagnostic likelihood of fossil species in relation to bone preservation potential. We use this to assess whether a taxonomic identification bias impacts the non-avian theropod fossil record. We find that the patterns of differential species abundance and clade diversity are not a consequence of their relative diagnosability. Although there are other factors that bias the theropod fossil record that are not investigated here, our results suggest that patterns of relative abundance and diversity for theropods might be more representative of Mesozoic ecology than often considered.     

Fossil biogeography: a new model to infer dispersal,extinction and sampling from palaeontological data

Methods in historical biogeography have revolutionized our ability to infer the evolution of ancestral geographical ranges from phylogenies of extant taxa, the rates of dispersals, and biotic connectivity among areas. However, extant taxa are likely to provide limited and potentially biased information about past biogeographic processes, due to extinction, asymmetrical dispersals and variable connectivity among areas. Fossil data hold considerable information about past distribution of lineages, but suffer from largely incomplete sampling. Here we present a new dispersal–extinction–sampling (DES) model, which estimates biogeographic parameters using fossil occurrences instead of phylogenetic trees. The model estimates dispersal and extinction rates while explicitly accounting for the incompleteness of the fossil record. Rates can vary between areas and through time, thus providing the opportunity to assess complex scenarios of biogeographic evolution. We implement the DES model in a Bayesian framework and demonstrate through simulations that it can accurately infer all the relevant parameters. We demonstrate the use of our model by analysing the Cenozoic fossil record of land plants and inferring dispersal and extinction rates across Eurasia and North America. Our results show that biogeographic range evolution is not a time-homogeneous process, as assumed in most phylogenetic analyses, but varies through time and between areas. In our empirical assessment, this is shown by the striking predominance of plant dispersals from Eurasia into North America during the Eocene climatic cooling, followed by a shift in the opposite direction, and finally, a balance in biotic interchange since the middle Miocene. We conclude by discussing the potential of fossil-based analyses to test biogeographic hypotheses and improve phylogenetic methods in historical biogeography.     

The fossil record and macroevolutionary history of the beetles

Coleoptera (beetles) is the most species-rich metazoan order, with approximately 380 000 species. To understand how they came to be such a diverse group, we compile a database of global fossil beetle occurrences to study their macroevolutionary history. Our database includes 5553 beetle occurrences from 221 fossil localities. Amber and lacustrine deposits preserve most of the beetle diversity and abundance. All four extant suborders are found in the fossil record, with 69% of all beetle families and 63% of extant beetle families preserved. Considerable focus has been placed on beetle diversification overall, however, for much of their evolutionary history it is the clade Polyphaga that is most responsible for their taxonomic richness. Polyphaga had an increase in diversification rate in the Early Cretaceous, but instead of being due to the radiation of the angiosperms, this was probably due to the first occurrences of beetle-bearing amber deposits in the record. Perhaps, most significant is that polyphagan beetles had a family-level extinction rate of zero for most of their evolutionary history, including across the Cretaceous–Palaeogene boundary. Therefore, focusing on the factors that have inhibited beetle extinction, as opposed to solely studying mechanisms that may promote speciation, should be examined as important determinants of their great diversity today.     

How diversification rates and diversity limits combine to create large-scale species-area relationships

Species-area relationships (SARs) have mostly been treated from an ecological perspective, focusing on immigration, local extinction and resource-based limits to species coexistence. However, a full understanding across large regions is impossible without also considering speciation and global extinction. Rates of both speciation and extinction are known to be strongly affected by area and thus should contribute to spatial patterns of diversity. Here, we explore how variation in diversification rates and ecologically mediated diversity limits among regions of different sizes can result in the formation of SARs. We explain how this area-related variation in diversification can be caused by either the direct effects of area or the effects of factors that are highly correlated with area, such as habitat diversity and population size. We also review environmental, clade-specific and historical factors that affect diversification and diversity limits but are not highly correlated with region area, and thus are likely to cause scatter in observed SARs. We present new analyses using data on the distributions, ages and traits of mammalian species to illustrate these mechanisms; in doing so we provide an integrated perspective on the evolutionary processes shaping SARs.     

Fossil evidence for key innovations in the evolution of insect diversity

Explaining the taxonomic richness of the insects, comprising over half of all described species, is a major challenge in evolutionary biology. Previously, several evolutionary novelties (key innovations) have been posited to contribute to that richness, including the insect bauplan, wings, wing folding and complete metamorphosis, but evidence over their relative importance and modes of action is sparse and equivocal. Here, a new dataset on the first and last occurrences of fossil hexapod (insects and close relatives) families is used to show that basal families of winged insects (Palaeoptera, e.g. dragonflies) show higher origination and extinction rates in the fossil record than basal wingless groups (Apterygota, e.g. silverfish). Origination and extinction rates were maintained at levels similar to Palaeoptera in the more derived Polyneoptera (e.g. cockroaches) and Paraneoptera (e.g. true bugs), but extinction rates subsequently reduced in the very rich group of insects with complete metamorphosis (Holometabola, e.g. beetles). Holometabola show evidence of a recent slow-down in their high net diversification rate, whereas other winged taxa continue to diversify at constant but low rates. These data suggest that wings and complete metamorphosis have had the most effect on family-level insect macroevolution, and point to specific mechanisms by which they have influenced insect diversity through time.     

The fossil record of extant elasmobranchs

Sharks and their relatives (Elasmobranchii) are highly threatened with extinction due to various anthropogenic pressures. The abundant fossil record of fossil taxa has allowed the tracing of the evolutionary history of modern elasmobranchs to at least 250 MYA; nonetheless, exactly how far back the fossil record of living taxa goes has never been collectively surveyed. In this study, the authors assess the representation and extent of the fossil record of elasmobranchs currently living in our oceans by collecting their oldest records and quantifying first appearance dates at different taxonomic levels (i.e., orders, families, genera and species), ecological traits (e.g., body size, habitat and feeding mechanism) and extinction risks (i.e., threatened, not threatened and data deficient). The results of this study confirm the robust representation of higher taxonomic ranks, with all orders, most of the families and over half of the extant genera having a fossil record. Further, they reveal that 10% of the current global species diversity is represented in the geological past. Sharks are better represented and extend deeper in time than rays and skates. While the fossil record of extant genera (e.g., the six gill sharks, Hexanchus) goes as far back as c. 190 MYA, the fossil record of extant species (e.g., the sand shark, Carcharias taurus Rafinesque 1810) extends c. 66 MYA. Although no significant differences were found in the extent of the fossil record between ecological traits, it was found that the currently threatened species have a significantly older fossil record than the not threatened species. This study demonstrate that the fossil record of extant elasmobranchs extends deep into the geologic time, especially in the case of threatened sharks. As such, the elasmobranch geological history has great potential to advance the understanding of how species currently facing extinction have responded to different stressors in the past, thereby providing a deep-time perspective to conservation.     

Developmental bias,macroevolution, and the fossil record

A fuller understanding of the role of developmental bias in shaping large‐scale evolutionary patterns requires integrating bias (the probability distribution of variation accessible to an ancestral phenotype) with clade dynamics (the differential survival and production of species and evolutionary lineages). This synthesis could proceed as a two‐way exchange between the developmental data available to neontologists and the strictly phenotypic but richly historical and dynamic data available to paleontologists. Analyses starting in extant populations could aim to predict macroevolution in the fossil record from observed developmental bias, while analyses starting in the fossil record, particularly the record of extant species and lineages, could aim to predict developmental bias from macroevolutionary patterns, including the broad range of extinct phenotypes. Analyses in multivariate morphospaces are especially effective when coupled with phylogeny, theoretical and developmental models, and diversity–disparity plots. This research program will also require assessing the “heritability” of an ancestral bias across phylogeny, and the tendency for bias change in strength and orientation over evolutionary time. Such analyses will help find a set of general rules for the macroevolutionary effects of developmental bias, including its impact on and interactions with the other intrinsic and extrinsic factors governing the movement, expansion, and contraction of clades in morphospace.