A 3D geometric morphometric analysis of the bovid distal humerus, with special reference to Rusingoryx atopocranion (Pleistocene, Eastern Africa).

The family Bovidae [Mammalia: Artiodactyla] is speciose and has extant representatives on every continent, forming key components of mammal communities. For these reasons, bovids are ideal candidates for studies of ecomorphology. In particular, the morphology of the bovid humerus has been identified as highly related to functional variables such as body mass and habitat. This study investigates the functional morphology of the bovid distal humerus in isolation due to its increased likelihood of preservation in the fossil record, and the resulting opportunity for a better understanding of the ecomorphology of extinct bovids. A landmark scheme of 30 landmarks was used to capture the 3D distal humerus morphology in 111 extant bovid specimens. We find that the distal humerus has identifiable morphologies associated with body mass, habitat preference and tribe affiliation and that some characteristics are shared between high body mass bovids and those living on hard, flat terrain which is likely due to the high stress on the bone in both cases. We directly apply our findings regarding extant bovids to the extinct alcelaphine bovid, Rusingoryx atopocranion from the mid to late Pleistocene (>33-45 ka) Lake Victoria region of Kenya. This species is known for some peculiar morphologies including a domed cranium with hollow nasal crests, and having small hooves for a bovid of its size. Another interesting aspect of Rusingoryx's skeletal morphology which has not been addressed is an unusual protrusion on the lateral epicondyle of the distal humerus. Despite considerable individual variation in the Rusingoryx specimens, we find evidence to support its historical assignment to the tribe Alcelaphini, and that it likely preferred open grassland habitats, which is consistent with independent reconstructions of the palaeoenvironment. We also provide the most accurate body mass estimate for Rusingoryx to date, based on distal humerus centroid size. Overall, we are able to conclude that the distal humerus in extant bovids is highly informative regarding body mass, habitat preference and tribe, and that this can be applied directly to a fossil taxon with promising results.

Behavioral and phylogenetic correlates of limb length proportions in extant apes and monkeys: Implications for interpreting hominin fossils.

The body proportions of extant animals help inform inferences about the behaviors of their extinct relatives, but relationships between body proportions, behavior, and phylogeny in extant primates remain unclear. Advances in behavioral data, molecular phylogenies, and multivariate analytical tools make it an opportune time to perform comprehensive comparative analyses of primate traditional limb length proportions (e.g., intermembral, humerofemoral, brachial, and crural indices), body size-adjusted long bone proportions, and principal components. In this study we used a mix of newly-collected and published data to investigate whether and how the limb length proportions of a diverse sample of primates, including monkeys, apes, and modern humans, are influenced by behavior and phylogeny. We reconfirm that the intermembral index, followed by the first principal component of traditional limb length proportions, is the single most effective variable distinguishing hominoids and other anthropoids. Combined limb length proportions and positional behaviors are strongly correlated in extant anthropoid groups, but phylogeny is a better predictor of limb length proportion variation than of behavior. We confirm convergences between members of the Atelidae and extant apes (especially Pan), members of the Hylobatidae and Pongo, and a potential divergence of Presbytis limb proportions from some other cercopithecoids, which correlate with adaptations for forelimb-dominated behaviors in some colobines. Collectively, these results substantiate hypotheses indicating that extinct hominins and other hominoid taxa can be distinguished by analyzing combinations of their limb length proportions at different taxonomic levels. From these results, we hypothesize that fossil skeletons characterized by notably disparate limb length proportions are unlikely to have exhibited similar behavioral patterns.

Relative abundance of grazing and browsing herbivores is not a direct reflection of vegetation structure: Implications for hominin paleoenvironmental reconstruction.

The diet of fossil herbivores inferred from enamel stable carbon isotopes is often used to make paleoenvironmental reconstructions. While many studies have focused on using environmental indicator taxa to make paleoenvironmental reconstructions, community-based approaches are considered to provide a more complete picture of paleolandscapes. These studies assume that the diet and relative abundance of herbivores are related to the areal extent of different vegetation types on the landscape. Here, we quantitatively test this assumption in 16 modern ecosystems in eastern and southern Africa with a wide range of woody vegetation cover. We conducted a landscape-level spatial analysis of vegetation patterns using a published land cover data set and computed landscape metrics. We compiled data on relative abundance and diet of herbivores inferred from carbon isotope studies for all large herbivores in these ecosystems. We found that despite differences in the total areal extent of different vegetation types, numerous sizable patches of each vegetation type are available in most ecosystems. However, despite variation across the ecosystems examined, grazers are typically the most abundant herbivores even in sites that have a higher proportion of forest and shrub cover. This indicates that the diet and relative abundance of herbivores is not a simple reflection of the total areal extent of vegetation types available on the landscape. The higher proportion of grazers observed in these ecosystems is a result of multiple factors including habitat heterogeneity, differences in biomass turnover rate between grasses and woody vegetation, resource partitioning, and the advantages of group living in open environments. Comparison of diet and relative abundance of herbivores in modern ecosystems to fossil herbivore assemblages shows that very different vegetation regimes can support similar herbivore assemblages. This study has significant implications for paleolandscape reconstructions and cautions against a simplistic wooded vs. grassland paleoenvironmental interpretations based on fossil herbivore assemblages.

No sustained increase in zooarchaeological evidence for carnivory after the appearance of Homo erectus.

The appearance of Homo erectus shortly after 2.0 Ma is widely considered a turning point in human dietary evolution, with increased consumption of animal tissues driving the evolution of larger brain and body size and a reorganization of the gut. An increase in the size and number of zooarchaeological assemblages after the appearance of H. erectus is often offered as a central piece of archaeological evidence for increased carnivory in this species, but this characterization has yet to be subject to detailed scrutiny. Any widespread dietary shift leading to the acquisition of key traits in H. erectus should be persistent in the zooarchaeological record through time and can only be convincingly demonstrated by a broad-scale analysis that transcends individual sites or localities. Here, we present a quantitative synthesis of the zooarchaeological record of eastern Africa from 2.6 to 1.2 Ma. We show that several proxies for the prevalence of hominin carnivory are all strongly related to how well the fossil record has been sampled, which constrains the zooarchaeological visibility of hominin carnivory. When correcting for sampling effort, there is no sustained increase in the amount of evidence for hominin carnivory between 2.6 and 1.2 Ma. Our observations undercut evolutionary narratives linking anatomical and behavioral traits to increased meat consumption in H. erectus, suggesting that other factors are likely responsible for the appearance of its human-like traits.

Intrataxonomic trends in herbivore enamel δ13C are decoupled from ecosystem woody cover.

Analysis of enamel stable carbon isotopes (δ13Cenamel) of fossil herbivores is an important tool for making inferences about Plio-Pleistocene vegetation structure in Africa and the environmental context of hominin evolution. Many palaeoecological studies implicitly or explicitly assume that individual variation in C3-C4 plant consumption among fossil herbivores directly reflects the abundance of C3 (trees, shrubs) or C4 (low-altitude tropical grasses) vegetation. However, a strong link between δ13Cenamel of herbivores and ecosystem vegetation structure has not been rigorously established. Here we combine δ13Cenamel data from a large dataset (n = 1,643) with multidecadal Landsat estimates of C3 woody cover across 30 African ecosystems to show that there is little relationship between intrataxonomic variation in δ13Cenamel and vegetation structure. This is especially true when removing forested ecosystems (>80% woody cover)-which numerous lines of evidence suggest are rare in the Plio-Pleistocene fossil record of eastern Africa-from our analyses. Our findings stand in contrast with the common assumption that variation in herbivore δ13Cenamel values reflects changes in the relative abundance of C3-C4 vegetation. We conclude that analyses using herbivore δ13Cenamel data to shed light on the environmental context of hominin evolution should look to explicitly community-level approaches for making vegetation inferences.

The environments of Australopithecus anamensis at Allia Bay, Kenya: A multiproxy analysis of early Pliocene Bovidae.

Australopithecus anamensis, among the earliest fully bipedal hominin species, lived in eastern Africa around 4 Ma. Much of what is currently known about the paleoecology of A. anamensis comes from the type locality, Kanapoi, Kenya. Here, we extend knowledge of the range of environments occupied by A. anamensis by presenting the first multiproxy paleoecological analysis focusing on Bovidae excavated from another important locality where A. anamensis was recovered, locality 261-1 (ca. 3.97 Ma) at Allia Bay, East Turkana, Kenya. Paleoenvironments are reconstructed using astragalar ecomorphology, mesowear, hypsodonty index, and oxygen and carbon isotopes from dental enamel. We compare our results to those obtained from Kanapoi. Our results show that the bovid community composition is similar between the two fossil assemblages. Allia Bay and Kanapoi bovid astragalar ecomorphology spans the spectrum of modern morphologies indicative of grassland, woodland, and even forest-adapted forms. Dietary reconstructions based on stable isotopes, mesowear, and hypsodonty reveal that these bovids' diet encompassed the full C3 to C4 dietary spectrum and overlap in the two data sets. Our results allow us to confidently extend our reconstructions of the paleoenvironments of A. anamensis at Kanapoi to Allia Bay, where this pivotal hominin species is associated with heterogeneous settings including habitats with varying degrees of tree cover, including grasslands, bushlands, and woodlands.

Investigating Biotic Interactions in Deep Time.

Recent renewed interest in using fossil data to understand how biotic interactions have shaped the evolution of life is challenging the widely held assumption that long-term climate changes are the primary drivers of biodiversity change. New approaches go beyond traditional richness and co-occurrence studies to explicitly model biotic interactions using data on fossil and modern biodiversity. Important developments in three primary areas of research include analysis of (i) macroevolutionary rates, (ii) the impacts of and recovery from extinction events, and (iii) how humans (Homo sapiens) affected interactions among non-human species. We present multiple lines of evidence for an important and measurable role of biotic interactions in shaping the evolution of communities and lineages on long timescales.

Mammal functional diversity and habitat heterogeneity: Implications for hominin habitat reconstruction.

Hominin habitats are frequently described as 'mosaic' based on interpretations of fossil assemblages comprising taxa with divergent functional adaptations (e.g., both grazers and browsers). This interpretation rests on an assumption that mammal functional diversity is positively associated with habitat heterogeneity. We test this assumption using modern mammal data for 141 sites in Africa. Species average body mass and locomotor and dietary information was compiled for all species >500 g. The functional diversity of each species assemblage was measured using five metrics: locomotor richness, trophic richness, functional richness, functional divergence, and functional evenness. We used a high-resolution woody cover estimate for sub-Saharan Africa to compute the coefficient of variation of percentage of woody cover for each site. We used a published land cover classification to compute the number of habitat patches and the number of distinct habitat types at each site. Multiple regressions were conducted at 9 different spatial resolutions (pixel size ranging from 100 m to 30 km) to explore the relationship between functional diversity and habitat heterogeneity metrics. The overall number of species found at a site is strongly positively associated with functional richness, locomotor richness, and trophic richness at all spatial resolutions. The number of habitat types at a site and the density of habitat patches show a modest positive relationship with most functional diversity metrics at most spatial resolutions. The coefficient of variation in woody cover is a very poor predictor of functional diversity. The locomotor and trophic richness of mammal communities are positively associated with habitat heterogeneity and may be useful for reconstructing aspects of heterogeneity at hominin sites. However, the overall number of species is an important confounding variable, which in turn depends on the sample size of the overall fossil sample. Researchers should carefully consider impacts of sample size on faunal reconstructions of habitat heterogeneity.

Reorganization of surviving mammal communities after the end-Pleistocene megafaunal extinction.

Large mammals are at high risk of extinction globally. To understand the consequences of their demise for community assembly, we tracked community structure through the end-Pleistocene megafaunal extinction in North America. We decomposed the effects of biotic and abiotic factors by analyzing co-occurrence within the mutual ranges of species pairs. Although shifting climate drove an increase in niche overlap, co-occurrence decreased, signaling shifts in biotic interactions. Furthermore, the effect of abiotic factors on co-occurrence remained constant over time while the effect of biotic factors decreased. Biotic factors apparently played a key role in continental-scale community assembly before the extinctions. Specifically, large mammals likely promoted co-occurrence in the Pleistocene, and their loss contributed to the modern assembly pattern in which co-occurrence frequently falls below random expectations.

Comparative isotopic evidence from East Turkana supports a dietary shift within the genus Homo.

It has been suggested that a shift in diet is one of the key adaptations that distinguishes the genus Homo from earlier hominins, but recent stable isotopic analyses of fossils attributed to Homo in the Turkana Basin show an increase in the consumption of C4 resources circa 1.65 million years ago, significantly after the earliest evidence for Homo in the eastern African fossil record. These data are consistent with ingesting more C4 plants, more animal tissues of C4 herbivores, or both, but it is also possible that this change reflects factors unrelated to changes in the palaeobiology of the genus Homo. Here we use new and published carbon and oxygen isotopic data (n = 999) taken from large-bodied fossil mammals, and pedogenic carbonates in fossil soils, from East Turkana in northern Kenya to investigate the context of this change in the isotope signal within Homo. By targeting taxa and temporal intervals unrepresented or undersampled in previous analyses, we were able to conduct the first comprehensive analysis of the ecological context of hominin diet at East Turkana during a period crucial for detecting any dietary and related behavioural differences between early Homo (H. habilis and/or H. rudolfensis) and Homo erectus. Our analyses suggest that the genus Homo underwent a dietary shift (as indicated by δ13Cena and δ18Oena values) that is (1) unrelated to changes in the East Turkana vegetation community and (2) unlike patterns found in other East Turkana large mammals, including Paranthropus and Theropithecus. These data suggest that within the Turkana Basin a dietary shift occurred well after we see the first evidence of early Homo in the region.

Phylogenetic signal in tooth wear dietary niche proxies.

In the absence of independent observational data, ecologists and paleoecologists use proxies for the Eltonian niches of species (i.e., the resource or dietary axes of the niche). Some dietary proxies exploit the fact that mammalian teeth experience wear during mastication, due to both tooth-on-tooth and food-on-tooth interactions. The distribution and types of wear detectible at micro- and macroscales are highly correlated with the resource preferences of individuals and, in turn, species. Because methods that quantify the distribution of tooth wear (i.e., analytical tooth wear methods) do so by direct observation of facets and marks on the teeth of individual animals, dietary inferences derived from them are thought to be independent of the clade to which individuals belong. However, an assumption of clade or phylogenetic independence when making species-level dietary inferences may be misleading if phylogenetic niche conservatism is widespread among mammals. Herein, we test for phylogenetic signal in data from numerous analytical tooth wear studies, incorporating macrowear (i.e., mesowear) and microwear (i.e., low-magnification microwear and dental microwear texture analysis). Using two measures of phylogenetic signal, heritability (H2) and Pagel's λ, we find that analytical tooth wear data are not independent of phylogeny and failing to account for such nonindependence leads to overestimation of discriminability among species with different dietary preferences. We suggest that morphological traits inherited from ancestral clades (e.g., tooth shape) influence the ways in which the teeth wear during mastication and constrain the foods individuals of a species can effectively exploit. We do not suggest that tooth wear is simply phylogeny in disguise; the tooth wear of individuals and species likely varies within some range that is set by morphological constraints. We therefore recommend the use of phylogenetic comparative methods in studies of mammalian tooth wear, whenever possible.

Digital data collection in paleoanthropology.

Understanding patterns of human evolution across space and time requires synthesizing data collected by independent research teams, and this effort is part of a larger trend to develop cyber infrastructure and e-science initiatives. At present, paleoanthropology cannot easily answer basic questions about the total number of fossils and artifacts that have been discovered, or exactly how those items were collected. In this paper, we examine the methodological challenges to data integration, with the hope that mitigating the technical obstacles will further promote data sharing. At a minimum, data integration efforts must document what data exist and how the data were collected (discovery), after which we can begin standardizing data collection practices with the aim of achieving combined analyses (synthesis). This paper outlines a digital data collection system for paleoanthropology. We review the relevant data management principles for a general audience and supplement this with technical details drawn from over 15 years of paleontological and archeological field experience in Africa and Europe. The system outlined here emphasizes free open-source software (FOSS) solutions that work on multiple computer platforms; it builds on recent advances in open-source geospatial software and mobile computing.

Taphonomy of fossils from the hominin-bearing deposits at Dikika, Ethiopia.

Two fossil specimens from the DIK-55 locality in the Hadar Formation at Dikika, Ethiopia, are contemporaneous with the earliest documented stone tools, and they collectively bear twelve marks interpreted to be characteristic of stone tool butchery damage. An alternative interpretation of the marks has been that they were caused by trampling animals and do not provide evidence of stone tool use or large ungulate exploitation by Australopithecus-grade hominins. Thus, resolving which agents created marks on fossils in deposits from Dikika is an essential step in understanding the ecological and taphonomic contexts of the hominin-bearing deposits in this region and establishing their relevance for investigations of the earliest stone tool use. This paper presents results of microscopic scrutiny of all non-hominin fossils collected from the Hadar Formation at Dikika, including additional fossils from DIK-55, and describes in detail seven assemblages from sieved surface sediment samples. The study is the first taphonomic description of Pliocene fossil assemblages from open-air deposits in Africa that were collected without using only methods that emphasize the selective retention of taxonomically-informative specimens. The sieved assemblages show distinctive differences in faunal representation and taphonomic modifications that suggest they sample a range of depositional environments in the Pliocene Hadar Lake Basin, and have implications for how landscape-based taphonomy can be used to infer past microhabitats. The surface modification data show that no marks on any other fossils resemble in size or shape those on the two specimens from DIK-55 that were interpreted to bear stone tool inflicted damage. A large sample of marks from the sieved collections has characteristics that match modern trampling damage, but these marks are significantly smaller than those on the DIK-55 specimens and have different suites of characteristics. Most are not visible without magnification. The data show that the DIK-55 marks are outliers amongst bone surface damage in the Dikika area, and that trampling is not the most parsimonious interpretation of their origin.

Paleoenvironments of the Shungura Formation (Plio-Pleistocene: Ethiopia) based on ecomorphology of the bovid astragalus.

The Shungura Formation in southwestern Ethiopia preserves a long and relatively continuous record of eastern African mammalian and hominin evolution. This study reconstructs habitat preferences of the Shungura Formation bovids from ca. 3.4-1.9 Ma (million years ago), based on the ecomorphology of fossil bovid astragali. Habitat predictions are made using a Discriminant Function Analysis informed by functional analysis of the astragalus that controls for body size and phylogenetic signal. The high abundance of astragali in the Shungura record allows for habitat reconstructions on sub-unit timescales, rather than aggregating samples at the level of geological member. During much of the time period examined, the lower Omo Valley was dominated by relatively mesic habitats in close proximity to the ancestral Omo River. Previous research has suggested that environmental change at ∼ 2.85 Ma caused significant habitat change at Shungura, but the astragalar data do not support major habitat change at this time. However, there are significant fluctuations in inferred habitats, with more open habitats indicated in sub-units C-08 and C-09 (∼ 2.56 Ma), and in unit F-01 (just after 2.36 Ma). The cause of sub-unit fluctuations in habitat signal is difficult to determine, but local factors, including spatial habitat heterogeneity around the river channel and in its floodplain, may play a contributing role. In the light of faunal evidence from elsewhere in the greater Omo-Turkana Basin, hominins in the area likely had access to a variety of heterogeneous habitats.

Functional morphology of the bovid astragalus in relation to habitat: controlling phylogenetic signal in ecomorphology.

Bovid astragali are one of the most commonly preserved bones in the fossil record. Accordingly, astragali are an important target for studies seeking to predict the habitat preferences of fossil bovids based on bony anatomy. However, previous work has not tested functional hypotheses linking astragalar morphology with habitat while controlling for body size and phylogenetic signal. This article presents a functional framework relating the morphology of the bovid astragalus to habitat-specific locomotor ecology and tests four hypotheses emanating from this framework. Highly cursorial bovids living in structurally open habitats are hypothesized to differ from their less cursorial closed-habitat dwelling relatives in having (1) relatively short astragali to maintain rotational speed throughout the camming motion of the rotating astragalus, (2) a greater range of angular excursion at the hock, (3) relatively larger joint surface areas, and (4) a more pronounced "spline-and-groove" morphology promoting lateral joint stability. A diverse sample of 181 astragali from 50 extant species was scanned using a Next Engine laser scanner. Species were assigned to one of four habitat categories based on the published ecological literature. A series of 11 linear measurements and three joint surface areas were measured on each astragalus. A geometric mean body size proxy was used to size-correct the measurement data. Phylogenetic generalized least squares (PGLS) was used to test for differences between habitat categories while controlling for body size differences and phylogenetic signal. Statistically significant PGLS results support Hypotheses 1 and 2 (which are not mutually exclusive) as well as Hypothesis 3. No support was found for Hypothesis 4. These findings confirm that the morphology of the bovid astragalus is related to habitat-specific locomotor ecology, and that this relationship is statistically significant after controlling for body size and phylogeny. Thus, this study validates the use of this bone as an ecomorphological indicator.

Ecomorphology and phylogenetic risk: Implications for habitat reconstruction using fossil bovids.

Reconstructions of paleohabitats are necessary aids in understanding hominin evolution. The morphology of species from relevant sites, understood in terms of functional relationships to habitat (termed ecomorphology), offers a direct link to habitat. Bovids are a speciose radiation that includes many habitat specialists and are abundant in the fossil record. Thus, bovids are extremely common in ecomorphological analyses. However, bovid phylogeny and habitat preference are related, which raises the possibility that analyses linking habitat with morphology are not 'taxon free' but 'taxon-dependent.' Here we analyze eight relative dimensions and one shape index of the metatarsal for a sample of 72 bovid species and one antilocaprid. The selected variables have been previously shown to have strong associations with habitat and to have functional explanations for these associations. Phylogenetic generalized least squares analyses of these variables, including habitat and size, resulted in estimates for the parameter lambda (used to model phylogenetic signal) varying from zero to one. Thus, while phylogeny, morphology, and habitat all march together among the bovids, the odds that phylogeny confounds ecomorphological analyses may vary depending on particular morphological characteristics. While large values of lambda do not necessarily indicate that habitat differences are unimportant drivers of morphology, we consider the low value of lambda for relative metatarsal width suggestive that conclusions about habitat built on observations of this particular morphology carry with them less 'phylogenetic risk.' We suggest that the way forward for ecomorphology is grounded in functionally relevant observations and careful consideration of phylogeny designed to bracket probable habitat preferences appropriately. Separate consideration of different morphological variables may help to determine the level of 'phylogenetic risk' attached to conclusions linking habitat and morphology.

Phylogenetic comparative methods complement discriminant function analysis in ecomorphology.

In ecomorphology, Discriminant Function Analysis (DFA) has been used as evidence for the presence of functional links between morphometric variables and ecological categories. Here we conduct simulations of characters containing phylogenetic signal to explore the performance of DFA under a variety of conditions. Characters were simulated using a phylogeny of extant antelope species from known habitats. Characters were modeled with no biomechanical relationship to the habitat category; the only sources of variation were body mass, phylogenetic signal, or random "noise." DFA on the discriminability of habitat categories was performed using subsets of the simulated characters, and Phylogenetic Generalized Least Squares (PGLS) was performed for each character. Analyses were repeated with randomized habitat assignments. When simulated characters lacked phylogenetic signal and/or habitat assignments were random, <5.6% of DFAs and <8.26% of PGLS analyses were significant. When characters contained phylogenetic signal and actual habitats were used, 33.27 to 45.07% of DFAs and <13.09% of PGLS analyses were significant. False Discovery Rate (FDR) corrections for multiple PGLS analyses reduced the rate of significance to <4.64%. In all cases using actual habitats and characters with phylogenetic signal, correct classification rates of DFAs exceeded random chance. In simulations involving phylogenetic signal in both predictor variables and predicted categories, PGLS with FDR was rarely significant, while DFA often was. In short, DFA offered no indication that differences between categories might be explained by phylogenetic signal, while PGLS did. As such, PGLS provides a valuable tool for testing the functional hypotheses at the heart of ecomorphology.