For well over 50 years, the telome theory of Walter Zimmermann has been extremely influential in interpreting the evolutionary history of land plant architecture. Using the “telome/mesome” distinction, and the concept of universal “elementary processes” underlying the change in form in all plants, the theory was an ambitious synthesis based on the proposition that evolutionary change might be understood by a simple set of developmental or evolutionary rules. However, a major problem resides in deciding exactly how assertions of change are to span both developmental and evolutionary domains simultaneously, and, we argue, the theory critically fails testability as a scientific theory. Thus, despite continued popularity for the descriptive terms derived from the theory in evolutionary studies of early land plants, time has come to replace it with a more explicit, testable approach. Presented here is an attempt to clarify perhaps the most important issue raised by the telome theory—whether simple changes in basic developmental processes suffice to describe much of early land plant evolution. Considering the morphology of Silurian–Devonian fossil members, it is hypothesized that early land plants possessed a common set of developmental processes centered on primary growth of shoot apical meristems. Among these were (1) the capacity to monitor and act upon internal physiological status here modeled as “apex strength,” (2) a mechanism for allocation of apex strength in a context-dependent way at each point of branching, (3) a rule for context-dependent apex angle for branches, (4) a largely invariant phyllotaxis unrelated to physiological status, and (5) a simple switch for terminating primary growth, based in part on genetics. Implemented as a set of developmental rules within a simple L-system model, these aspects of primary development in plants determine a sizable range of resultant morphologies, some of which are highly reminiscent of the early fossils. Thus, some support is found, perhaps, for Zimmermann's intuition. However, traditional concepts of growth patterns in plants, including the contrast between epidogenesis and apoxogenesis, require updating. In our reformulation, developmental processes, stated as rules of developmental dynamics, together constitute what we term the plant's developmental state. Using a hypothetico-deductive format, one may hypothesize intrinsic (or genetic) developmental processes that play out as realized developmental activity in specific spatial/temporal contexts, as modified by multiple context factors. The resultant plant morphology is highly dependent on multiple and simultaneous pathway ontogenetic trajectories. Within a likely set of developmental rules reasonably inferred from plant development, some of Zimmermann's elementary processes are perhaps recognizable whereas others are not. Progressively “overtopped” morphologies are easily produced by modifying intrinsic branch allocation. However, even so, the other developmental rules have a profound effect on final architectures. Planate architectures and circination vernation, often treated as special cases by plant morphologists, are perhaps better understood in terms of recurrent or iterative developmental relationships. Much analytic work remains before a completely specified system of rules will emerge. A well-articulated relationship between ontogeny and phylogeny remains fundamentally important in assessing evolutionary change. Fossil and living plants make it abundantly clear that current evolutionary concepts involving modification of a single ontogenetic trajectory from ancestor to descendant need to be greatly expanded into consideration of the entire logical geometry of causation in development. A mechanism for testing is also required that need not wait for complete elucidation at the molecular level. The relative simplicity of plant development, combined with an outstanding fossil record of early members, offers unique opportunities along these lines.

A number of natural events can cause ozone depletion, including asteroid and comet impacts, large-scale volcanism involving the stratospheric injection of chlorine, and close cosmic events such as supernovae. These events have previously been postulated to have been sole or contributory causes of mass extinctions. Following such events, UV-B radiation would have been elevated at the surface of the earth. The possibilities for detecting elevated UV-B as a kill mechanism in the fossil record are discussed. In the case of impact events and large-scale volcanism, the taxa affected by increases in UV-B radiation are likely to be similar to those affected by cooling and by the initial drop in irradiance caused by stratospheric dust injection. Thus UV-B may synergistically exacerbate the effects of these other environmental changes and contribute to stress in the biosphere, although UV-B alone is unlikely to cause a mass extinction. By the same token, however, this similarity in affected taxa is likely to make delineating the involvement of UV-B radiation in the fossil record more difficult. Cosmic events such as supernovae may produce smaller extinction events, but ones that are “cleaner” UV catastrophes without the involvement of other environmental changes.

Spinal ossified tendons are a defining character for Ornithischia, one of the two major clades of dinosaurs. The function of these bony rods has remained a mystery since their first detailed description in 1886. Qualitative approaches to understand ossified tendon function have resulted in different ecological and behavioral interpretations for ornithopod dinosaurs. To evaluate ossified tendon function, this study constructed finite element models of the vertebral column for two ornithopod taxa: Tenontosaurus, which shows the plesiomorphic condition of longitudinally arrayed tendons along the spinous processes, and Brachylophosaurus, which exhibits a lattice of tendons along the spinous processes. Both models predict that ossified tendons stiffened the vertebral column, especially the tail, but the derived lattice of ossified tendons in iguanodontoidean dinosaurs, like Brachylophosaurus, increased spinal stiffness more than the plesiomorphic condition. Caudofemoral muscles that retracted the hindlimb during locomotion attached the femur to the tail in ornithopods. Increased tail stiffness caused by intratendinous ossification may have influenced locomotion by rigidly anchoring M. caudofemoralis longus to the tail, thereby allowing a more forceful retraction of the hindlimb by reducing ventral flexion of the tail during muscle contraction. Ossified tendons may also have been important for storing elastic energy throughout the gait cycle. Moreover, the lattice of ossified tendons stiffened the trunk and tail nearly equally in Brachylophosaurus, indicating the evolution of a postural function by passively supporting the epaxial musculature in maintaining a horizontal vertebral column.

Decay experiments are becoming a more widespread tool in evaluating the fidelity of the fossil record. Character interpretations of fossil specimens stand to benefit from an understanding of how decay can result in changes in morphology and, potentially, total character loss. We performed a decay experiment for the Class Enteropneusta to test the validity of anatomical interpretations of the Burgess Shale enteropneust Spartobranchus tenuis and to determine how the preservation of morphological features compares with the sequence of character decay in extant analogues. We used three species of enteropneust (Saccoglossus pusillus, Harrimania planktophilus, and Balanoglossus occidentalis) representing the two major families of Enteropneusta. Comparisons between decay sequences suggest that morphological characters decay in a consistent and predictable manner within Enteropneusta, and do not support the hypothesis of stemward slippage. The gill bars and nuchal skeleton were the most decay resistant, whereas the gill pores and pre-oral ciliary organ were unequivocally the most decay prone. Decay patterns support the identification of the nuchal skeleton, gill bars, esophageal organ, trunk, and proboscis in Spartobranchus tenuis and corroborate a harrimaniid affinity. Bias due to the taphonomic loss of taxonomically informative characters is unlikely. The morphologically simple harrimaniid body plan can be seen, therefore, to be plesiomorphic within the enteropneusts. Discrepancies between the sequence of decay in a laboratory setting and fossil preservation also exist. These discrepancies are highlighted not to discredit the use of modern decay studies but rather to underline their non-actualistic nature. Paleoenvironmental variables besides decay, such as the timeframe between death and early diagenesis as well as postmortem transport, are discussed relative to decay data. These experiments reinforce the strength of a comprehensive understanding of decay sequences as a benchmark against which to describe fossil taxa and understand the conditions leading to fossilization.

Truth, goes an old proverb, is the daughter of time. Fifty years ago, G. G. Simpson (1944) brought paleontology into the Neodarwinian fold, arguing that evolutionary tempo can be documented in the geological record and used to inform debate about evolutionary mode. Today, increasingly sophisticated paleontological investigations of rate—be it diversification, extinction, migration, morphological change, or divergence in macromolecular sequence—require calibration of the geological time scale with a precision far greater than Simpson could have anticipated. Expanding research on the relationships between environmental history and evolution also requires unprecedented resolution in correlation and geochronometry.

Several groups of siliceous microorganisms possess a skeletal latticework of interconnected rods. Skeletal configurations of one of these groups, the silicoflagellates, can be produced by a simple mathematical model that minimizes the apical surface area for a given basal area and internal volume. A similar model that minimizes the total length of the skeletal elements, and thus the silica utilization and skeletal weight, produces configurations that are generally less common in silicoflagellates. The diversity of silicoflagellate skeletal morphologies suggests that both the minimization of apical surface area and the conservation of skeletal material may be important factors in skeletal design. The two most important morphologies found in modern oceans, the four-sided Dictyocha and the six-sided Distephanus, can co-occur in an environment where both factors have some relative importance. However, these models do not explain the range of silicoflagellate skeletal morphology found in nature.

The hypothesis is proposed that mean level of heterozygosity is functionally related to rate of speciation in evolutionary phylads. Under this hypothesis, phylads which speciate more rapidly do so because of increased level of within-species genetic variability which is then available to conversion to species differences under appropriate ecological or environmental conditions. An important corollary is that rate of speciation could be limited in phylads with low genetic variability, irrespective of environmental considerations.

This hypothesis has been tested with respect to electrophoretically detectable variation in products of structural genes in two families of North American fishes characterized by grossly different rates of speciation. Totals of 69 species of the highly speciose Cyprinidae, and 19 species of the relatively depauperate Centrarchidae, were assayed for mean level of heterozygosity at 11–24 genetic loci. Since Cyprinidae and Centrarchidae exhibit on the average nearly identical levels of genic variation (Ĥ = 0.052 ± 0.004, and Ĥ = 0.049 ± 0.009, respectively), the hypothesis that level of heterozygosity affects rate of speciation in these fishes is not supported.

Nonetheless, the amount of genic variability in both Cyprinidae and Centrarchidae is large, comparable to mean levels in previously studied vertebrates. The great wealth of genome variability, reflected in the electrophoretic variation present in virtually all outcrossing organisms, apparently can accommodate considerable flexibility in rate and pattern of evolutionary response to the various environmental regimes challenging organisms.

It has been suggested that a sudden event at the end of the Cretaceous period caused a major extinction that was felt disproportionately by creatures in the water column. It also has been argued that benthic deposit feeders, being relatively independent of abundance of organic particles in the water column, should have survived the crisis more readily than suspension feeders, which depended more upon feeding upon phytoplankton. I argue that the hypothesis of relative immunity of deposit feeders is insufficient, because deposit feeders by and large depend upon a supply of organic matter from the water column and would have succumbed to food shortage nearly as rapidly as suspension feeders, possibly within a maximum of three to six months. This near simultaneity of extinction would have been especially true of continental shelf environments. Even in some parts of the deep sea, it is likely that a dependence upon the water column above might have caused deep-sea deposit feeders to succumb rapidly. Therefore, deposit feeders would not necessarily have outsurvived suspension feeders during a crisis of depleted water-column phytoplankton, increased shading by inert particles, or poisoning of the water column and killing of phytoplankton. The relatively lower rate of extinction of nuculoid bivalves may relate instead to their presence in deeper-water refuge habitats, their apparent relative ability to diversify in higher latitudes, or their resistance to factors other than food shortage.

Body size is one of the most studied phenotypic attributes because it is biologically important and easily measured. Despite a long history of study, however, the pattern of body-size change in diverse higher taxa over the Phanerozoic remains largely unknown because few relevant data sets span more than a single geological period or provide comprehensive, global coverage. In this study, we measured representative specimens of 3414 brachiopod genera illustrated in the Treatise on Invertebrate Paleontology. We applied these size data to stage-resolved stratigraphic ranges from the Treatise and the Paleobiology Database to develop a Phanerozoic record of trends in brachiopod size. Using a model comparison approach, we find that temporal variation in brachiopod size exhibits two distinct modes—a Paleozoic mode of size increase and a post-Paleozoic mode indistinguishable from a random walk. This transition reflects a change in the identities of the most diverse brachiopod orders rather than a shift in mode within any given order. Paleozoic size increase reflects a small, persistent bias toward the origination of new genera larger than those surviving from the previous stage and is identifiable as a statistically supported trend in three orders representing both Class Strophomenata (Order Productida) and Class Rhynchonellata (orders Atrypida and Spiriferida). Extinction exhibits no consistent bias with respect to size. The shift in evolutionary mode across the end-Permian mass extinction adds to long-standing evidence from studies of diversity and abundance that this biotic catastrophe suddenly and permanently altered the evolutionary history of what was, until that time, the most diverse animal phylum on Earth.

Ecologists and paleontologists alike are increasingly using the fossil record as a spatial data set, in particular to study the dynamics and distribution of geographic range sizes among fossil taxa. However, no attempts have been made to establish how accurately range sizes and range-size dynamics can be preserved. Two fundamental questions are: Can common paleo range-size reconstruction methods accurately reproduce known species’ ranges from locality (i.e., point) data? And, are some reconstruction methods more reliable than others? Here, we develop a methodological framework for testing the accuracy of commonly used paleo range-size reconstruction methods (maximum latitudinal range, maximum great-circle distance, convex hull, and alpha convex hull) in different extinction-related biogeographic scenarios. We use the current distribution of surface water bodies as a proxy for “preservable area,” in which to test the performance of the four methods. We find that maximum great-circle distance and convex-hull methods most reliably capture changes in range size at low numbers of fossil sites, whereas convex hull performs best at predicting the distribution of “victims” and “survivors” in hypothetical extinction scenarios. Our results suggest that macroevolutionary and macroecological patterns in the relatively recent past can be studied reliably using only a few fossil occurrence sites. The accuracy of range-size reconstruction undoubtedly changes through time with the distribution and area of fossiliferous sediments; however, our approach provides the opportunity to systematically calibrate the quality of the spatial fossil record in specific environments and time intervals, and to delineate the conditions under which paleobiologists can reconstruct paleobiogeographical, macroecological, and macroevolutionary patterns over critical intervals in Earth history.

Parasitic trematode worms leave characteristic pits in their bivalve mollusk hosts and represent an ideal system for analyzing parasite-host interactions through space and time with statistically meaningful sample sizes. Previous work in Late Pleistocene–Holocene sequences from the Po plain revealed significant long-term fluctuations in trematode prevalence values: higher prevalence in retrogradational environments (TST) and negligible prevalence in progradational environments (HST). Here we expand upon this work by investigating traces of parasitism, kleptoparasitism, and predation on mollusk death assemblages from two domains along the northern Adriatic coastline. The domain north of the Po delta (TST-like) and the southern domain (including the Po delta; HST-like) comprise environments comparable to those recovered in late Holocene (<6 Kyr) subsurface progradational deposits. We collected 17,299 specimens representing 111 species from 11 locations on the northern Adriatic coast of Italy. Our results reveal high predation pressure, a high diversity of host taxa, and widespread presence of trematode infestation in starved, oligotrophic, environmentally more stable (i.e., TST-like) settings north of the Po delta. Immediately south of the Po delta, in settings with strong and variable sedimentary input, almost no infestation is recorded. The reappearance of infestation is evident in the southern portion of the study area (i.e., Cattolica-Montemarciano), relatively far from the highly stressed environments south of the Po River. There is no significant difference in trematode prevalence values between fossil and modern samples. The distribution of spionid traces (an indicator of stressed environments) was nearly the opposite of that displayed by trematodes. Drilling frequency is highest in TST-like environments and is not correlated with diversity indices. These results suggest that temporal trends of trematode prevalence (and possibly also other biotic interactions) in sedimentary successions are controlled by environmental changes driven by glacio-eustatic dynamics, and reaffirm the importance of interpreting temporal trends in the context of spatial variation.

The Tulip Beds locality on Mount Stephen (Yoho National Park, British Columbia) yields one of the most abundant and diverse (~10,000 specimens in 110 taxa) Burgess Shale fossil assemblages in the Canadian Rockies. Detailed semi quantitative and quantitative analyses of this assemblage suggest strong similarities with the Walcott Quarry on Fossil Ridge. Both assemblages are dominated by epibenthic, sessile, and suspension feeding taxa, mostly represented by arthropods and sponges and have comparable diversity patterns, despite sharing only about half the genera. However, the Tulip Beds has a higher relative abundance of suspension feeders and taxa of unknown affinity compared to the Walcott Quarry. These biotic variations are probably largely attributable to ecological and evolutionary differences between the two temporally distinct communities that adapted to similar, but not identical, environmental settings. For instance, the Tulip Beds is farther away from the Cathedral Escarpment than the Walcott Quarry. The Tulip Beds and Walcott Quarry assemblages are more similar to each other than either one is to the assemblages of the Chengjiang biota, although the relative diversity of major taxonomic groups and ecological patterns are similar in all assemblages. The conserved diversity patterns and ecological structures among sites suggest that the ecological composition of Cambrian Burgess Shale-type communities was relatively stable across wide geographic and temporal scales.

We describe three giant palaeobatrachid fossil tadpoles of the genus Palaeobatrachus (Nieuwkoop-Faber [NF] stages 60–64) from the Miocene of Randecker Maar, Germany. The largest was 150 mm at the beginning of metamorphosis (stage 60), whereas the smallest was 100 mm and approaching the end of metamorphosis (stage 64). In contrast, normal palaeobatrachid tadpoles and their pipid relatives, both extinct and extant, rarely exceed 60 mm in length. We review here both ecological and pathological conditions that are conducive to the development of gigantism in tadpoles. Tadpoles that lack a thyroid gland become exceptionally large and arrest development at early hindlimb stages (NF stages 53–56). However, the advanced metamorphic stages of the giant Palaeobatrachus tadpoles indicate that they were able to metamorphose, and thus were not athyroid. Environmental factors—pond size and permanence, predators, duration of the growing season—may all contribute to tadpole gigantism in certain extant anuran species. We identify suites of ecological features that distinguish extant anurans with large tadpoles from high-latitude and high-altitude permanent lakes in temperate regions (e.g., certain Rana and Telmatobius) from tropical species, such as Pseudis paradoxa, whose tadpoles normally achieve large size in temporary seasonal ponds. The paleoecology of Randecker Maar suggests that Palaeobatrachus tadpoles lived in a permanent semitropical lake, but one with few predators.

Stephen Jay Gould made impressive contributions to macroevolutionary theory; one of the topics in this area that particularly interested him was how to define and recognize species selection. Here we explore how and why Gould's ideas on concepts related to species selection evolved over 30 years, from the punctuated equilibria paper of 1972 to his “Structure of Evolutionary Theory” magnum opus published in 2002. Throughout his career his ideas on species selection shifted between three phases. Initially, Gould favored a definition of species selection that was more descriptive. Later, he came to distinguish between species sorting, which he called species selection in the broad sense, and true species selection, which is tied to the concept of species-level aptations. Finally, he came to view species selection in a broader, more inclusive way, effectively merging the two earlier viewpoints. His ideas on species selection changed over the years because he was trying to square his views on complex concepts like adaptation, natural selection, emergence, and the independence of macroevolutionary theory. Gould's thoughts on species selection not only help to define the history of debate on the concept but also help set a course for the future.

Large-scale trends in planktonic foraminiferal diversity have so far been based on utilization of synoptic biostratigraphic range charts. Although this approach ensures the taxonomic consistency and quality of the data being used, it takes no formal account of any sampling biases that might exist in the fossil record. We demonstrate that the occurrence data of planktonic foraminifera, as recorded in the primary literature, are strongly biased by sampling. We do this by demonstrating that raw diversity curves derived from the land-based and deep-sea records are strikingly different, but that they each correlate with the intensity of sampling in their respective environments, and thus are ultimately controlled by the structure of the geological record in each setting. Because sampling of the Mesozoic record is best in our land record whereas sampling of the Cenozoic is best in our deep-sea record, we combine the two to generate the best-supported estimates of species and genus diversity over time from these data. We correct for sampling bias using shareholder quorum subsampling and a modeling approach. The data are then transformed to generate a range-through plot of species richness that is compared with two earlier estimates of the diversity history where comparable species-in-bin data can be recovered. No robust statistical correlation is found among the three estimates. Although differences in amplitude are to be expected, differences in the actual shape of the curve are surprising. We conclude that these differences stem from the nature of the data themselves, namely the taxonomic scheme adopted and the taxonomic coverage used.

The fossil record has been used to show that in some geologic intervals certain traits of taxa may increase their survivability, and therefore that the risk of extinction is not randomly distributed among taxa. It has also been suggested that traits that buffer against extinction in background times do not confer the same resistance during mass extinction events. An open question is whether at any time in geologic history extinction probabilities were randomly distributed among taxa. Here we use a method for detecting random extinction to demonstrate that during both background and mass extinction times, extinction of marine invertebrate genera has been nonrandom with respect to species richness categories of genera. A possible cause for this nonrandom extinction is selective clustering of extinctions in genera consisting of species which possess extinction-biasing traits. Other potential causes considered here include geographic selectivity, increased extinction susceptibility for species in species-rich genera, or biases related to taxonomic practice and/or sampling heterogeneity. An important theoretical result is that extinction selectivity at the species level cannot be smoothly extrapolated upward to genera; the appearance of random genus extinction with respect to species richness of genera results when extinction has been highly selective at the species level.

Phyletic speciation arising from a shift between adaptive peaks in the selective surface is modeled. At a locus initially monomorphic, a new allele is introduced that triggers phenotypic evolution by making a new adaptive peak available. A phenotypic model is used to describe the evolution of a quantitative character that is important for fitness but not directly affected by the new allele. Recursion equations are derived for both the allele frequency and the mean of the quantitative character. It is shown that the dynamics of the quantitative trait fit a “punctuated equilibrium” pattern of evolution. Phyletic speciation is seen to occur over a time span so short as to make observation of the transition in the fossil record unlikely.

Attempts to use Fourier methods to quantify ammonoid suture shape have failed to yield robust, repeatable results because sutures are complex curves that violate the assumptions of Fourier mathematics. In particular, sampled sutures are artificially truncated such that the first and last sampled points are not equivalent, the folds are non-stationary, and a position along a horizontal reference axis may map to multiple amplitudes along the suture path. Here I introduce an alternative Fourier method—the windowed short-time Fourier transform (STFT)—that accommodates these characteristics of complex curves. For each suture, digitized landmarks were parameterized using a tangent angle function and then smoothed by convolving with an apodization function. Piece-wise Fourier transforms were then calculated and averaged, resulting in a robust, unique quantification of line morphology. STFT coefficients and estimated power spectra describing the relative weights of harmonics were generated for 576 Paleozoic-basal Triassic ammonoid genera, representative of the range of suture morphotypes. While insensitive to major episodes of taxonomic turnover (Frasnian/Famennian, end-Devonian, and Permian/Triassic extinctions), the summed power data support the previously observed trend toward increasing suture complexity through time. Moreover, partitioning the summed power statistic into harmonic ranges allows novel insight into Paleozoic suture evolution. In particular, the data show significant shifts in the dominant morphotypes during periods of rebound and radiation and suggest that basal Triassic ammonoids possessed unique suture morphotypes when compared with those of the Paleozoic.

Macroevolutionary and macroecological studies must account for biases in the fossil record, especially when questions concern the relative abundance and diversity of taxa that differ in preservation and sampling potential. Using Cenozoic marine mollusks from a temperate setting (New Zealand), we find that much of the long-term temporal variation in gastropod versus bivalve occurrences is correlated with the stage-level sampling probabilities of aragonitic versus calcitic taxa. Average sampling probabilities are higher for calcitic species, but this contrast is time-varying in a predictable way, being concentrated in stages with widespread carbonate deposition.

To understand these results fully, we link them with analyses at the level of individual point occurrences. Doing so reveals that aragonite bias is effectively absent in terrigenous clastic sediments. In limestones, by contrast, calcitic species have at least twice the odds of sampling as aragonitic species. This result is most pronounced during times of widespread carbonate deposition, where the difference in the per-collection odds of sampling species is a factor of eight. During carbonate-rich intervals, calcitic taxa also have higher odds of sampling in clastics. At first glance this result may suggest simple preservational bias against aragonite. However, comparing relative odds of aragonitic versus calcitic sampling with absolute sampling rates shows that the positive calcite bias during carbonate-rich times reflects higher than average occurrence rates for calcitic taxa (rather than lower rates for aragonitic taxa) and that the negative aragonite bias in limestones reflects lower than average occurrence rates for aragonitic taxa (rather than higher rates for calcitic taxa).

Our results therefore indicate a time-varying interplay of two main factors: (1) taphonomic loss of aragonitic species in carbonate sediments, with no substantial bias in terrigenous clastics; and (2) an ecological preference of calcitic taxa for environments characteristic of periods with pervasive carbonate deposition, irrespective of lithology per se.

Organismal metabolic rates reflect the interaction of environmental and physiological factors. Thus, calcifying organisms that record growth history can provide insight into both the ancient environments in which they lived and their own physiology and life history. However, interpreting them requires understanding which environmental factors have the greatest influence on growth rate and the extent to which evolutionary history constrains growth rates across lineages. We integrated satellite measurements of sea-surface temperature and chlorophyll-a concentration with a database of growth coefficients, body sizes, and life spans for 692 populations of living marine bivalves in 195 species, set within the context of a new maximum-likelihood phylogeny of bivalves. We find that environmental predictors overall explain only a small proportion of variation in growth coefficient across all species; temperature is a better predictor of growth coefficient than food supply, and growth coefficient is somewhat more variable at higher summer temperatures. Growth coefficients exhibit moderate phylogenetic signal, and taxonomic membership is a stronger predictor of growth coefficient than any environmental predictor, but phylogenetic inertia cannot fully explain the disjunction between our findings and the extensive body of work demonstrating strong environmental control on growth rates within taxa. Accounting for evolutionary history is critical when considering shells as historical archives. The weak relationship between variation in food supply and variation in growth coefficient in our data set is inconsistent with the hypothesis that the increase in mean body size through the Phanerozoic was driven by increasing productivity enabling faster growth rates.