Close banner

2022-05-14 20:38:19 By : Ms. Rosa Song

Thank you for visiting You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

Communications Biology volume  4, Article number: 1172 (2021 ) Cite this article

Diverse lines of geological and geochemical evidence indicate that the Eocene-Oligocene transition (EOT) marked the onset of a global cooling phase, rapid growth of the Antarctic ice sheet, and a worldwide drop in sea level. Paleontologists have established that shifts in mammalian community structure in Europe and Asia were broadly coincident with these events, but the potential impact of early Oligocene climate change on the mammalian communities of Afro-Arabia has long been unclear. Here we employ dated phylogenies of multiple endemic Afro-Arabian mammal clades (anomaluroid and hystricognath rodents, anthropoid and strepsirrhine primates, and carnivorous hyaenodonts) to investigate lineage diversification and loss since the early Eocene. These analyses provide evidence for widespread mammalian extinction in the early Oligocene of Afro-Arabia, with almost two-thirds of peak late Eocene diversity lost in these clades by ~30 Ma. Using homology-free dental topographic metrics, we further demonstrate that the loss of Afro-Arabian rodent and primate lineages was associated with a major reduction in molar occlusal topographic disparity, suggesting a correlated loss of dietary diversity. These results raise new questions about the relative importance of global versus local influences in shaping the evolutionary trajectories of Afro-Arabia’s endemic mammals during the Oligocene.

Deep sea cores consistently preserve evidence for a strong oxygen isotope shift (Oi-1, ~33.5 Ma (million years ago)) just after the Eocene-Oligocene boundary (EOB, 33.9 Ma), the timing of which closely correlates with the expansion of Antarctic glaciation and drops in sea level, calcite compensation depth, and sea surface temperatures1,2,3,4. These abrupt environmental changes in the marine realm roughly coincide with marked shifts in mammalian community structure on some landmasses, most notably Europe5,6 and Asia7,8,9. Overall, however, terrestrial records from the early Oligocene provide geographically mixed evidence for environmental change relative to late Eocene conditions10,11. The nature of the Eocene-Oligocene transition (EOT) on the Afro-Arabian landmass has long been enigmatic, but recent paleontological sampling has improved our understanding of this phase12,13,14,15,16,17; this enriched context has inspired proposals of increased geographic provincialism of mammalian communities in the early Oligocene of Africa13, and led to the identification of multiple strepsirrhine primate extinctions in northern Africa near the EOB18.

The inconsistent signal for early Oligocene cooling across the planet suggests that some areas might have experienced gradual environmental change, if any10,11. Detection of a long-term biotic response to such gradual environmental change can be complicated by sampling gaps in the fossil record, but Bayesian tip-dating (BTD) phylogenetic analyses19 that incorporate rates of morphological evolution to estimate divergence times among living and extinct taxa20,21,22,23,24,25,26,27 can partially overcome this problem by allowing the identification of “ghost” lineages that cross temporal and geographic sampling gaps. Here we employ these methods to study lineage diversity in Afro-Arabian primates and rodents, the most common small mammals in Paleogene fossil deposits of that landmass.

Two clades of rodents dispersed into Afro-Arabia during the Eocene—the now largely arboreal but ancestrally terrestrial Anomaluroidea (the scaly-tails or anomalures), which probably arrived in the early Eocene27,28, and terrestrial Hystricognathi, which likely dispersed to Afro-Arabia from Asia in the middle Eocene25,29. Arboreal primates mirror this pattern of staggered arrival times, with strepsirrhines likely dispersing from Europe or Asia to Afro-Arabia in the early Eocene, and anthropoids appearing in the middle Eocene following a dispersal from Asia22,26. These four clades include a great diversity of living and extinct species that have lived in a wide range of arboreal and terrestrial habitats that could have been impacted by environmental change during the early Oligocene.

We constructed a near-comprehensive composite time-scaled phylogeny that includes living and extinct members of all four clades, with divergences among extinct taxa based largely on BTD analyses (see “Methods”). The resulting composite tree includes 317 tips, 62% of which (198) are extinct (Fig. 1a and Supplementary Fig. 1). A lineage-through-time (LTT) plot derived from the composite tree (Fig. 1b) reveals three distinct diversity peaks (latest Eocene (35 Ma), earliest Miocene (20–20.5 Ma), and Recent), and two major diversity troughs (mid-Oligocene (28–29.5 Ma) and mid-Miocene (13 Ma)).

a Composite time-scaled tree of Afro-Arabian primates and rodents used for calculation of lineage diversity shifts from the Eocene to Recent. See Supplementary Fig. 1 for full tree. b Black line is the lineage through time (LTT) plot based on the composite tree in (a); light blue histogram shows number of tips in the composite tree in 1 Ma long bins; top, high-resolution astronomically dated benthic carbon (δ13C‰) and oxygen (δ18O‰) isotope curves from the early Eocene to Recent, from Westerhold et al.4. δ18O‰ curve provides a proxy for mean temperature difference from the present day. c LTT plot for the composite tree compared to 10,000 random trees with the same tip ages as the composite tree, but random topologies and branch lengths. d LTT plot for Afro-Arabian Hyaenodonta, a clade of carnivorous placentals that were likely the primary mammalian predators of the primates and rodents in the composite tree. Dan. = Danian, Se. = Selandian, Th. = Thanetian, Bart. = Bartonian, Rup. = Rupelian, Chat. = Chattian, Aq. or Aqui. = Aquitanian, Burd. = Burdigalian, L. or Lan. = Langhian, S. or Ser. = Serravallian, Tort. = Tortonian, M or Mes. = Messian, UD = undivided.

The mid-Oligocene and mid-Miocene diversity troughs bottom out near sampling gaps in the Afro-Arabian terrestrial vertebrate fossil record (Fig. 1b), leaving open the possibility that the LTT curves are artificially low during these time periods simply due to a lack of sampling. This explanation can be rejected because the LTT curve shows high lineage diversification during all other periods for which there is little or no paleontological sampling (Fig. 1b). Furthermore, the major losses of lineage diversity that precede the mid-Oligocene and mid-Miocene diversity troughs occur during time periods that are particularly well-sampled14,16,30.

We tested whether the magnitude of extinction required by the composite tree could be observed in a randomized sample by generating 10,000 trees that had the same tip ages as those in the composite tree, but random topologies and branch lengths (Fig. 1c). On average, the random trees require a ~17% loss of lineages from peak Eocene diversity to the deepest point in the early Oligocene trough, whereas the composite tree requires a loss of ~61% of lineages across the same interval. The magnitude of the lineage loss observed in the composite tree therefore cannot be explained as an artifact of sampling alone; rather, it is highly dependent on the non-random phylogenetic structure of the composite tree, which reflects a strong lack of phylogenetic continuity (i.e., few ghost lineages) linking the similarly diverse late Eocene and early Miocene primate and rodent communities.

The lineage diversification pattern of Afro-Arabian primates and rodents closely matches that of Afro-Arabian members of Hyaenodonta, a clade of carnivorous mammals that is known to have gone extinct in the late Miocene23 (Fig. 1d). Prior to 15 Ma, Afro-Arabian hyaenodonts show the same LTT peaks as are seen in the composite tree of primates and rodents (latest Eocene and early Miocene), with a distinct mid-Oligocene diversity trough (Fig. 1d). Because Hyaenodonta goes extinct in the late Miocene, the veracity of the middle Miocene diversity drop that is also seen in primates and rodents is not open to debate. When hyaenodont lineage numbers are added to those in the composite tree of primates and rodents, the magnitude of total lineage loss from the late Eocene into the mid-Oligocene rises to ~63% (Supplementary Fig. 2). This value closely approximates the estimated worldwide loss of ~67% of all species during the transition from the late Eocene to the early Oligocene31.

The vast majority of the extinct species included in the composite tree of primates and rodents are only known from dental remains, and often only isolated teeth. Fortunately tooth shape reflects diet, and homology-free dental topographic metrics (DTMs) have proven useful for discriminating among primate and rodent species with different dietary preferences32,33,34. More generally, primates and rodents must be able to process food to maintain basic biological functions, so molar dental topography has a strong and unambiguous link to survival. We take disparity in DTM values through time to broadly reflect temporal shifts in dietary breadth. We calculated three DTMs (Dirichlet normal energy (ariaDNE)35, orientation patch count rotated (OPCR)34, and relief index (RFI)36) on lower second molars (M2) across a large subset (89%) of the extinct species older than 15 Ma for which M2 is known. Lineage counts between 55 and 15 Ma (taken at 0.5 Ma intervals) in the reduced sample on which we calculated DTMs closely track lineage counts in the composite tree (r2 = 0.97). We calculated ancestral states throughout the tree for the first two principal components (PCs) derived from PC analysis of all three DTMs, and calculated disparity within each major clade at 0.5 Ma intervals (via interpolation) using four different disparity metrics (square root of 2D hull area (SR2DHA), minimum spanning tree length (MST), sum of ranges (SoR), and sum of variances (SoV)). As our estimates of disparity are based on a principal components analysis (PCA) of DTM values, they are only capturing “functional richness” of M2 topography within each clade (sensu Villeger et al.)37.

Among the disparity metrics used here, SR2DHA is the least impacted by outliers, and allows us to sum values for each of the four clades considered independently (Fig. 2). A comparison of lineage number with sumSR2DHA and sumMST between 55 and 15 Ma shows very similar temporal trends (slope 0.98–1.02; r2 = 0.86–0.88). Extinctions in the early Oligocene resulted in a substantial loss of sumSR2DHA and sumMST, with values at 30.5 Ma being about one half of peak late Eocene values. Over the subsequent 10 Myr, sumSR2DHA and sumMST gradually increase, reaching an early Miocene peak at ~95% of late Eocene values.

The blue line is the LTT plot for the primate and rodent species for which M2 dental topographic metrics were calculated; the orange line shows temporal changes in the sum of the square roots of 2D hull areas for all four clades, with 2DHA calculated for each clade (Anomaluroidea, Anthropoidea, Hystricognathi, and Strepsirrhini) independently (sumSR2DHA); the green line shows temporal changes in the length of the minimum spanning tree (MST), calculated for each clade (Anomaluroidea, Anthropoidea, Hystricognathi, and Strepsirrhini) independently (sumMST). sumSR2DHA and sumMST show a tight correlation with lineage number, as shown by the r2 values in the insets (n = 80). Bart. = Bartonian, Aqui. = Aquitanian, Burd. = Burdigalian, Lan. = Langhian.

Rodents show a drop in dental topographic disparity (DTD) at the EOB regardless of which disparity metric is used (Supplementary Fig. 3), and regardless of whether Anomaluroidea and Hystricognathi are analyzed together (Supplementary Fig. 3), or independently (Supplementary Figs. 4 and 5). The drop in DTD is most precipitous among hystricognaths (Supplementary Fig. 4), whereas anomaluroids show a very gradual decline in DTD from the EOB into the early Miocene (Supplementary Fig. 3). Strepsirrhines show a small DTD dip from the latest Eocene to the EOB, but the overall trend through the Oligocene is for DTD to increase (Supplementary Figs. 6 and 7); however, this pattern is driven largely by the early Miocene stem daubentoniid Propotto, whose relatively flat molars are radically different from those of other Afro-Arabian strepsirrhines22. When Propotto is excluded, strepsirrhine DTD drops sharply just before the EOB, and then drops again at 32 Ma, with a slow subsequent increase that peaks in the early Miocene, but only to a maximum of half the DTD observed in the late Eocene (Supplementary Fig. 7). Anthropoidea is the only clade that shows an increase in DTD across the EOB, followed by large drops at 32 and 30 Ma, when DTD returns to its lowest level since the middle Eocene (i.e., shortly after Anthropoidea appeared in Afro-Arabia) (Supplementary Fig. 8). DTD then increases steadily through the late Oligocene into the early Miocene, but never reaches the peak seen at 32 Ma, even though lineage numbers are higher in the early Miocene than at any point in the Paleogene.

When all of the ancestral states for primate and rodent DTMs are pooled, long-term shifts toward relatively high OPCR values and relatively low RFI and ariaDNE values are evident in the mid-Oligocene, following periods of increased variability for both metrics in the late Eocene (Fig. 3). RFI values remain relatively low through the early Miocene, while ariaDNE values show increased variability from the late Oligocene into the early Miocene.

Violin plots of all ancestral reconstructions for ariaDNE, OPCR, and RFI, collated at 0.5 Ma intervals.

The centroid of the hystricognath dental topographic morphospace shifted farther and faster than that of any other clade considered here (movement of 0.28/Ma for hystricognaths versus 0.10–0.13/Ma for the other clades), and moved rapidly toward the ancestral anomaluroid morphospace during the middle and late Eocene (Fig. 4a). This shift suggests that early Afro-Arabian hystricognaths were evolving dental adaptations similar to those of early terrestrial stem anomaluroids such as zegdoumyids and nementchamyids, and so might have been exploiting similar dietary resources in terrestrial settings (Fig. 4b). The anomaluroid centroid shows a strong shift away from the group’s ancestral morphospace around 47–48 Ma, prior to the estimated ~44 Ma arrival time of Hystricognathi. The timing of this shift roughly coincides with the origin of crown Anomaluroidea27, members of which are either fully arboreal (Anomaluridae)38 or at least partially arboreal (Zenkerellidae)39, and so might reflect exploitation of novel dietary resources within an arboreal milieu. The evolution of arboreality presumably allowed crown anomaluroids to avoid direct competition with terrestrial hystricognaths for the remainder of the Cenozoic, but the same cannot be said of terrestrial or semi-terrestrial stem and basal crown anomaluroids, which were extirpated during the late Eocene interval when hystricognaths were undergoing a major adaptive radiation in Afro-Arabia.

a Plot of Euclidean movement of the centroid of each clade’s dental topographic morphospace from their inferred arrival time in Afro-Arabia (first divergence) to 15 Ma, and overall path length, along the first two principal components derived from PCA of the three dental topographic metrics (ariaDNE, OPCR, and RFI). b Comparison of Euclidean movement of morphospace centroids within each order (Primates and Rodentia), showing relative shifts through time; red and orange circles indicate the position of the Eocene-Oligocene and Oligocene-Miocene boundaries, respectively.

Prior to the EOB, Afro-Arabian anthropoids show little movement of their dental topographic (DT) morphospace (total Eocene Euclidean path length of 0.59, Fig. 4a, b), whereas that of strepsirrhines shifts rather dramatically over the same interval (Eocene path length of 1.97) (Fig. 4a, b). This pattern suggests the possibility of competitive displacement of middle and late Eocene strepsirrhines by anthropoids within their shared arboreal habitats. However, in the early Oligocene, anthropoid feeding ecology undergoes a major change, as reflected by a strong shift in the position of that clade’s DT morphospace (i.e., away from that of strepsirrhines) between 32 and 30 Ma, likely due in part to the extinction of smaller-bodied anthropoid lineages (Fig. 4a; see also Kirk and Simons40). With an apparent release of ecological competition between anthropoids and strepsirrhines during the Oligocene, the latter then undergo a shift back toward the ancestral anthropoid DT morphospace from which they were displaced in the Eocene (Fig. 4a).

Our data strongly suggest that Afro-Arabian primate and rodent clades lost most of their peak Eocene lineage diversity and DTD over the course of the first four million years of the Oligocene. All of the clades considered here (including Hyaenodonta) show a drop in lineage diversity within 500 Kyr of the onset of the early Oligocene global cooling event. There are no reliable high-resolution temperature proxies from continental Afro-Arabia across the EOT, so the possible influence of global cooling on this pattern of lineage loss is currently an open question; sea surface temperatures off the coast of Nigeria showed a ~2 °C drop across the EOT41, but estimates of mean annual temperatures based on poorly constrained equatorial terrestrial pollen records show no change during the same interval11. Afro-Arabia was largely isolated from northern continents at this time, so these extinctions cannot be attributed to competition with immigrant terrestrial taxa12,42. The same cannot be said of the middle Miocene diversity trough, by which time numerous Eurasian rodent clades had appeared in Africa43.

Anthropoidea, Hystricognathi, and Hyaenodonta all show a recovery after the EOB, followed by a second major loss of lineage diversity that bottoms out around 29.5–30 Ma. This time period is roughly coincident with land mammal age boundaries in North America (Whitneyan-Arikareean boundary, ~30 Ma)44, Europe (Suevian-Arvernian boundary, ~29.2 Ma)44, and South America (Tinguirirican-Deseadan boundary, ~30 Ma)45. All of these boundaries predate the Mid-Oligocene Glacial Interval (28–26.3), but are close to a major spike in atmospheric CO2 that peaks at ~29.5 Ma, followed by a steep drop through the rest of the Oligocene46.

The final extirpation of numerous Afro-Arabian mammal lineages around 30 Ma also roughly coincides with the onset of a series of major geological events that permanently and profoundly transformed the landscape of eastern Africa and the southwestern part of the Arabian peninsula—the uplift of the Ethiopian plateau47, volcanic super eruptions in Ethiopia at 31.1 and 30.8 Ma48, a massive outpouring of flood basalts of the Afro-Arabian Large Igneous Province (AALIP)48 (Fig. 2), and the opening of the Red Sea and Gulf of Aden49. Basalts of the AALIP are estimated to be ~900,000 km3 in volume47 and up to 3 km thick in some areas48. Unlike the effects of early Oligocene global cooling, which might have been buffered at lower latitudes, it is indisputable that these volcanic events must have had a devastating impact on the plant and animal communities that inhabited the AALIP, while also having a major influence on large-scale patterns of atmospheric moisture transport50 and water drainage51 in eastern Africa. It has recently been estimated that the release of sulfur and halogen molecules from the AALIP would have been comparable in volume to that which has been estimated for the Deccan and Siberian Traps48, both of which are temporally associated with mass extinction events (end-Cretaceous and end-Permian, respectively)52. Primates and rodents of early Oligocene age have been found around the periphery of the AALIP, indicating that this region was likely a major center of diversification for these clades. The only productive African fossil sites that fall within the mid-Oligocene diversity trough (the Chilga region of northwest Ethiopia53, which is geographically situated within the AALIP) have yielded numerous paenungulate species, but no rodents or primates, though the absence of small mammals at these localities might be due to taphonomic bias53. Nevertheless, the small number of late Oligocene sites that have been found in Afro-Arabia have thus far shown low primate and rodent diversity54,55, consistent with a scenario of early Oligocene extinction and late Oligocene reassembly of primate and rodent communities.

All of the foregoing suggests that explanations for the Oligocene extinctions that our analyses have revealed are likely complex and multifactorial, as there were global influences (such as cooling across the EOT, and changes in atmospheric CO2 levels) and local influences (such as AALIP volcanism) that could have caused major environmental changes in Afro-Arabia between ~34 and ~30 Ma. Any such scenarios must also take into account the varied evolutionary responses of other endemic clades during this turbulent interval—for instance, proboscideans diversified morphologically while embrithopods persisted into the late Oligocene relatively unchanged, at least dentally56; some hyracoid lineages that were present in the early Oligocene extended into the Miocene, but several went extinct in the early Oligocene57; and crown macroscelideans started to diversify in the Oligocene, but probably far to the southwest of the AALIP58. Crown bathyergoid rodents are estimated to have appeared in the mid-Oligocene in Afro-Arabia59, indicating the existence of selection pressures favoring fossoriality early in the Oligocene. It is tempting to suggest that hyaenodont extinctions in the early Oligocene might have followed from the loss of lineages that produced their preferred prey species; hyaenodont lineage diversity more closely tracks that of the relatively common anthropoids and hystricognaths (from 44 to 15 Ma; r2 = 0.81) than that of the relatively rare strepsirrhines and anomaluroids (from 54 to 15 Ma; r2 = 0.32). The lineage diversity and DTD patterns revealed here can be tested via similar analyses of these, and other, endemic African clades, but ultimately the more urgent need is improved paleontological sampling of the Afro-Arabian Oligocene. The environmental context for these pulses of extinction and diversification will only be clarified through a dedicated effort to extract paleotemperature estimates, stable isotopes, and paleobotanical data from diverse settings across Afro-Arabia that preserve terrestrial sediments of Oligocene age.

Our application of these phylogenetic and dental topographic methods to study long-term trends in lineage and dietary diversity is novel, and as DTMs are homology-free, this approach can be applied to any mammalian clade whose members retain functional molar teeth. Similar analyses of other mammalian clades evolving on other landmasses will therefore produce data that are directly comparable to those we present here, and we anticipate that such studies could provide significant new insights into the terrestrial mammalian response to global climate change through the Cenozoic.

The composite near-comprehensive tree of Afro-Arabian primates, anomaluroid rodents, and hystricognathous rodents combines information from multiple BTD analyses of character matrices published by Gunnell et al.22, Rasmussen et al.60, Dembo et al.61, DeSilva et al.62, Sallam and Seiffert25, and Marivaux et al.16. The tree is supplemented by molecular phylogenies of extant taxa59,63, and grafting of several extinct taxa that were not included in any of the aforementioned phylogenetic character matrices. All of the new BTD analyses employed the fossilized birth–death (=FBD) prior, flat beta priors for the extinctionpr and fossilizationpr parameters, exp(10) for the speciationpr parameter, and used the R code provided by Gunnell et al.22 to calculate an appropriate clockratepr prior (which uses as input the allcompat tree from a nonclock Bayesian analysis and the midpoint of each species’ age prior). Convergence of runs was assessed by tracking minimum effective sample sizes (=minESS) of all parameters, and average standard deviations of split frequencies (=ASDSF). Additional details of all phylogenetic analyses, and justifications for placements of grafted species, are available in the Supplementary Methods. The BTD analysis of hyaenodonts is based on a modified matrix that augments that of Borths and Stevens64 and includes a new hyaenodont genus and species from the terminal Eocene Fayum Locality 41. Phylogenetic data matrices, settings, and output are available in the Dryad repository associated with this study (

The composite time-scaled phylogeny includes 317 tips (198 fossil, 119 extant), all of which descend from unquestionably Afro-Arabian clades. The root age of this tree (representing the divergence between Primates and Rodentia) is ~76.16 Ma, based on the average of four different estimates for the origin of Euarchontoglires recovered by Springer et al.65. To quantify the number of lineages (i.e., branches) present in the tree at any given time, we used the function ltt.plot.coords() from the R package ape66 for continuous-age lineage through time (LTT) summarization.

With the set of 317 taxa, the function rtree() from the R package ape was used to create 10,000 random phylogenies. Branch lengths were stripped from these trees leaving only topological relationships. Each topology was then converted to a set of hard constraints using the function createMrBayesConstraints() from the R package paleotree67. Next, base R functions were used to create 10,000 Nexus character-taxon matrices, each containing the 317 taxa and 1,000 standard characters where all matrix cells were randomly coded as 0 or 2. After doing so, each set of topological constraints was paired with a unique character-taxon matrix. We used scripting with MrBayes v3.2.768 to conduct Bayesian time-scaled phylogenetic analyses. Three constraint types were applied to each of these calls: (1) a fixed topology (see above), (2) fixed ages for taxa in the matrix (extracted from the starting composite tree), and (3) a fixed root age. MCMC settings for all analyses specified 1 run, 8 chains, and sampling for 1 M generations in increments of 1000. Fixed root ages were paired with approximately flat clock rate priors [norm(1,100)]. All characters were treated as ordered for branch length estimation from the random character matrices (i.e., binary standard characters treated as unordered will less effectively inform branch lengths in MrBayes analyses). Proposals that include Tau moves were disabled because these are incompatible with the application of hard topological constraints for the full tree. Posterior distributions were summarized using the option for 50% majority rule plus all compatible groups. The resultant set contains 10,000 trees where topology, branch lengths, and node ages are random—but tip times preserve the ages of their respective taxa. Lastly, LTT summarization was applied to each tree in the randomized set (see Fig. 1c). Both the composite tree and randomized tree-set are included in the Dryad repository associated with this study:

A total of 329 second lower molar specimens of 134 species were scanned using a Nikon XTH 225 ST scanner (110–121 kVp, 116–120 mA, 0.013–0.034 voxel size), a Scanco vivaCT40 scanner (75kVp, 177 mA, 0.0125–0.035 voxel size), or a GE Phoenix Nanotom M scanner (90 kVp, 180 mA, 0.012 voxel size). Original fossil specimens or epoxy casts of fossil specimens were scanned. Surface files were created and oriented into occlusal view in Avizo version 8.0. Digital surface files were cropped at the lowest point of occlusal basin for calculation of ariaDNE and OPCR. All digital surfaces were also cropped along the base of the crown for the calculation of RFI. Minor deformities such as fine cracks in the enamel or small bubbles introduced during the molding process were manually removed in Geomagic. Surface files that had been cropped at the lowest point of the occlusal basin were simplified to 8000 triangles and smoothed for 30 iterations in Avizo (used for the calculation of ariaDNE and OPCR). Surfaces that had been cropped along the base of the crown were simplified to 10,000 triangles and smoothed for 30 iterations in Avizo (used for the calculation of RFI).

The variables ariaDNE (ε = 0.1), RFI, and OPCR were used to quantify dental shape. Dental topography was chosen to quantify dental shape as these GIS-based methods summarize the three-dimensional shape of a tooth in various functionally distinct metrics: complexity34, relief36, and bending energy69.

Complexity, or OPCR34,70 captures the number of “breakage sites”, or “patches” on a crown surface and is calculated as the OPCR. It is calculated by summing the number of “patches” that are formed by continuous points on the digital surface that face the same direction and have the same orientation34,70. The counting of patches is repeated eight times, every time rotating the crown surface slightly (5.625°). This “rotated” patch count decreases the error that is introduced by manually orientating the tooth into occlusal view. The OPCR value is the average of the eight Orientation Patch Counts. The more “patches”, the more “tools” the tooth has for breaking down food and thus the higher the complexity. OPCR was calculated for surfaces using MorphoTester71. Relief (or RFI) is measured as the ratio between 3D crown surface area and the 2D crown surface area36. RFI was calculated in Excel using the surface area and outline area output from MorphoTester following the same equation used by Boyer (2008):

Bending energy of a surface was calculated as the DNE of a surface, or ariaDNE35,69 and, simply put, captures the “curvedness” of the surface. ariaDNE was calculated using a the Matlab script provided by Shan et al.35 in Matlab version R2018a. The ε-value was set to 0.1, as we intended to capture large-scale features of molar shape but remained within the recommended ε-value range by Shan et al.35.

The topographic variables capture different aspects of the crown shape. Simply put, complexity reflects number of features on a surface, or “tools” (such as cusps, Evans et al., 2007)34, whereas curvature and relief reflect the shape of features (such as a sharp cusp compared to a blunt cusp, Bunn et al., 2011)69. However, in practice, all dental topographic variables are affected by both the number and shape of features of teeth, but they are affected to different degrees (see Winchester (2016)71 for a further discussion of how these different topographic variables correlate with each other or are affected by simulated differences in 3D shapes).

Dental topographic values were averaged per species. On average, each species was represented by 2.5 specimens. Specimens with extreme wear were excluded following criteria outlined by Allen et al.72. A list of all specimens used in the DTM analyses is provided in the Dryad repository associated with this study:

A PCA was applied to the species averages of the complete sample using the three variables ln(OPCR), RFI, and ariaDNE. The natural log of the OPCR values was taken as count-data, and cannot be treated as continuous data73. The PCA was performed in R using the prcomp() function from the stats package and setting “scale=TRUE”, which standardizes all variables prior to the analysis. PC1 captured 46.9% of the variation, PC2 38.7%, and PC3 14.4%.

Ancestral state reconstruction (ASR) analyses were run independently for the three DTMs (ariaDNE, OPCR, and RFI), and jointly for the first two PCs from the PCA of these metrics. Again, together PC1 and PC2 explain ~85.5% of the observed variation. Before ASR, the data for each DTM were z-transformed (i.e., centered and scaled with means and standard deviations). In total, the trait data include 134 euarchontan taxa—47 anthropoid primates (43 fossil, 4 extant), 24 strepsirrhine primates (19 fossil, 5 extant), 47 hystricognath rodents (43 fossil, 4 extant), and 16 anomaluroid rodents (12 fossil, 4 extant). A time-scaled tree that corresponds to these taxa was created by pruning tips from the larger composite phylogeny leaving only terminals with data affirmative DTMs.

We tested the application of multiple Brownian motion (BM) rate matrices to explore if DTM evolution among the four focus clades (anthropoids, strepsirrhines, hystricognaths, and anomaluroids) could be better modeled in a multivariate context. To do so, we defined each group as its own “regime” and then created all possible combinations thereof. These combinations include: one version with a single BM rate (=regime), seven versions with two rates, six versions with three rates, and one version with four rates. For the two-PC trait set, we also created three variants for each version—a variant where the two trait rate matrices are estimated independently (i.e., unique variances and shared covariance; unconstrained), a variant where the rate (=sigma) for each PC is estimated independently but covariance is disallowed (i.e., a “diagonal” constraint), and a variant where the two PCs share a single rate (i.e., an “equal” constraint, without covariance). In total, this strategy produced 15 combinatorial versions for each DTM data set and 45 combinatorial versions/variants for the two-PC data set.

The function mvBM() (multivariate Brownian motion) from the R package mvMORPH74 was used to obtain AICc decision-criterion scores for all combinatorial versions (and their variants). The AICc calculation includes penalty terms for the invocation of additional model parameters which allows for optimal model choice, without overfitting. For each of the four trait-sets (ariaDNE, OPCR, RFI, and two-PC), model tests on the version-suites were repeated ten times and the models producing the lowest AICc scores were selected as preferred.

For the tests using two-PC trait data, all ten replicates identified the best-fit model as including two BM rate matrices (1 = anomaluroids, 2 = anthropoids + strepsirrhines + hystricognaths) treated as unconstrained (i.e., unique rates for each PC with shared covariance). For ariaDNE, the ten replicates produced mixed results, with a near tie between two models—one that includes four BM rates (i.e., one rate per clade) and one that includes three BM rates (1 = hystricognaths, 2 = anomaluroids, 3 = anthropoids + strepsirrhines). The ten OPCR model tests also produced mixed results with optimal AICc scores corresponding to a three-rate model (1 = anthropoids, 2 = hystricognaths, 3 = strepsirrhines + anomaluroids) and a two-rate model (1 = anthropoids + hystricognaths, 2 = strepsirrhines + anomaluroids). For RFI, the ten tests found two best-fit models which incorporate four rates (i.e., one rate per clade) and three rates (1 = hystricognaths, 2 = anomaluroids, 3 = anthropoids + strepsirrhines). Using the phylogenetic tree, the four trait data sets, and these model test results, we used the function estim() from the R package mvMORPH to estimate continuous-trait values for all internal tree nodes. For each of the four trait data sets, we repeated these ASR estimations ten times (i.e., once per model test result). Final nodal estimates were calculated as the arithmetic means from the ten rounds of ASR.

Within a phylogenetic tree, each edge (i.e., branch) is defined by (and connected by) a unique pair of nodes, the parent and child. Tree tips are considered terminal tree nodes for which there are no descending children. For a time-scaled phylogenetic tree, each parent–child pair has an estimated start and end age. If the same tree has been used in an ASR analysis, each parent–child pair will also be associated with start and end values for the reconstructed trait. Thus, a simple linear rate of change can be calculated for each tree edge and the trait’s value can be interpolated at all positions along any edge in the tree.

Further, at any given age contained in a time-scaled tree, all edges (i.e., lineages) that cross that time-slice can be counted. This is discrete-age lineage through time (LTT) summarization. Similarly, at any given age we can collect a distribution of interpolated trait values from the set of edges that cross the time-slice. If this strategy is repeated at small age increments, the distributions of trait values from multiple time-slices can characterize (nearly continuously) the evolution of variation in the trait through time.

We produced novel R code to implement time-slice trait value interpolation. Input for the code includes a time-scaled phylogenetic tree (ours is discussed above), observed trait values for tree tips (i.e., input for the ASR analyses), and estimated trait values for internal tree nodes (i.e., output from the ASR analyses). The R code also includes an input argument that specifies age increment sampling. We selected increments at 0.5 Myr which samples our tree at 153 time-slices ranging from present (0 Ma) to the late Cretaceous (76 Ma). The output generated by our code is a simple table where column names are time-slice ages and row names are tree edge numbers (following the standard edge numbering scheme used for class “phylo” objects in R). Table cells contain the interpolated values of the trait for tree edges at their corresponding time-slice. Where edges do not cross a time-slice, table cells are populated with not applicable (NA). We applied this code to all three DTMs and both PCs.

The phylogeny used for ASR analyses contains 12 tree edges that cannot (firmly) be associated with residency on the Afro-Arabian landmass. These edges are the stem branches leading to the clades (1) primates, (2) anthropoids, (3) strepsirrhines, (4) rodents, (5) hystricognaths, and (6) anomaluroids. Further, the branches (7,8) leading to the terminals for Lemur catta and Nycticebus coucang (extant taxa from Madagascar and Indonesia, respectively) almost certainly do not represent occupancy on the African continent. We also identified the stem branches (9,10,11,12) associated with the clades (a) Azibiidae + Djebelemuridae, (b) Adapidae, and (c) Nosmips + crown Strepsirrhini as possibly non-African. With the intent to limit our characterizations to lineages that can be confidently referred to Afro-Arabia, we elected to exclude these tree edges from the results—prior to summarization—by removing their corresponding rows from the trait value interpolation tables. We assigned each of the remaining tree edges to one of the focus clades creating four primary edge subsets (anthropoids, strepsirrhines, hystricognaths, and anomaluroids), and then regrouped two of these subsets into a “rodents” set, two into a “primates” set, and all four back into an “all clades” set. One additional exploratory subset removed the Propotto edge from the strepsirrhine group.

We used the “all clades” edge set to summarize ariaDNE, OPCR, and RFI for time-slice samples. These distributions were characterized using kernel density estimation (KDE) and plotted as violins (Fig. 3). Summarization of PC1 + PC2 was applied independently to the four focus subsets (i.e., clades), to the strepsirrhine (less Propotto) subset, to the recombined “rodents” set, and to the recombined “primates” set (Supplementary Figs. 3–8). Each of these summarizations begins with the oldest data affirmative time-slice and ends at present day (0 Ma). We calculated four measures of variation for each time-slice: SoV, SoR, SR2DHA, and MST. SoV is simply the variance of PC1 added to the variance of PC2. Similarly, SoR is the range of PC1 (difference between max and min) added to the range of PC2. SR2DHA treats PC1 + PC2 coordinates in each time-slice distribution as an XY point cloud and a hull is estimated around the points. The area of the resulting polygon quantifies the distribution’s occupancy of morphospace. To find SR2DHA, we used the function concaveman() from the R package concaveman75 with the concavity argument set to 10,000—a convex setting for our data which ranges approximately from −3 to 3 on both axes. Then, area was calculated using the function Polygon from the R package sp76, followed by square root transformation. MST point associations were found using the function mst() from the R package ape. The lengths of each of these linkages were calculated as standard Euclidean distance and all lengths in the network were summed to yield the MST metric.

Additionally, we computed per time-slice point cloud centroids (where X and Y are the arithmetic means of PC1 and PC2, respectively) and calculated the Euclidean distance between point cloud centroids of adjacent time-slices (Fig. 4a, b). Finally, we counted the number of lineages present in each time-slice sample. As a total measure of PC1 + PC2 variation, we summed the independent SR2DHAs of the four focus clades at each time-slice (Fig. 2).

Original data, model tests, ASR results, interpolation tables, and summarizations are available in the Dryad repository associated with this study ( We also provide all novel R scripts (and functions) as well as a step-by-step R workflow that reproduces this portion of the research.

The principal components analysis was performed in R using the prcomp() function from the stats package and setting “scale=TRUE”, which standardizes all variables prior to the analysis. Correlation analyses were performed in R. All input files for phylogenetic analyses, ancestral state reconstructions, time-slice trait value interpolation, and time-slice summarization and calculation of disparity through time are available on the Dryad Digital Repository ( and include all of the settings required to reproduce our results.

Further information on research design is available in the Nature Research Reporting Summary linked to this article.

Input data files, settings, code and results from phylogenetic, ASR, and disparity analyses are available on the Dryad Digital Repository ( Digital surface models for all of the figured fossil specimens are available on MorphoSource (, except those housed at the National Museums of Kenya.

All novel R code (R version 4.0.2) used in this study is available on the Dryad Digital Data Repository (

Liu, Z. et al. Global cooling during the Eocene-Oligocene climate transition. Science 323, 1187–1190 (2009).

CAS  PubMed  Article  Google Scholar 

Coxall, H. K., Wilson, P. A., Pälike, H., Lear, C. H. & Backman, J. Rapid stepwise onset of Antarctic glaciation and deeper calcite compensation in the Pacific Ocean. Nature 433, 53–57 (2005).

CAS  PubMed  Article  Google Scholar 

Zachos, J., Pagani, M., Sloan, L., Thomas, E. & Billups, K. Trends, rhythms, and aberrations in global climate 65 Ma to present. Science 292, 686–693 (2001).

CAS  PubMed  Article  PubMed Central  Google Scholar 

Westerhold, T. et al. An astronomically dated record of Earth’s climate and its predictability over the last 66 million years. Science 369, 1383–1387 (2020).

CAS  PubMed  Article  Google Scholar 

Stehlin, H. G. Remarques sur les faunules de mammifère des couches éocènes et Oligocènes du Bassin de Paris. Bull. Soc. Géol. Fr. 19, 488–520 (1909).

Escarguel, G. & Legendre, S. New methods for analysing deep-time meta-community dynamics and their application to the Paleogene mammals from the Quercy and Limagne area (Massif Central, France). Strata 13, 245–273 (2006).

Meng, J. & McKenna, M. C. Faunal turnovers of Palaeogene mammals from the Mongolian Plateau. Nature 394, 364–367 (1998).

Sun, J. et al. Synchronous turnover of flora, fauna, and climate at the Eocene–Oligocene boundary in Asia. Sci. Rep. 4, 7463 (2014).

CAS  PubMed  PubMed Central  Article  Google Scholar 

Ni, X., Li, Q., Li, L. & Beard, K. C. Oligocene primates from China reveal divergence between African and Asian primate evolution. Science 352, 673–677 (2016).

CAS  PubMed  Article  Google Scholar 

Sheldon, N. D. et al. Coupling of marine and continental oxygen isotope records during the Eocene-Oligocene transition. Geol. Soc. Am. Bull. 128, 502–510 (2016).

Pound, M. J. & Salzmann, U. Heterogeneity in global vegetation and terrestrial climate change during the late Eocene to early Oligocene transition. Sci. Rep. 7, 43386 (2017).

CAS  PubMed  PubMed Central  Article  Google Scholar 

Seiffert, E. R. Early primate evolution in Afro-Arabia. Evol. Anthropol. 21, 239–253 (2012).

Beard, K. C., Coster, P. M. C., Salem, M. J., Chaimanee, Y. & Jaeger, J. J. Biogeographic provincialism shown by Afro-Arabian mammals during the middle Cenozoic: climate change, Red Sea rifting and global eustasy. In Human Interaction with the Environment in the Red Sea (eds. Agius, D. A., Khalil, E., Scerri, E. M. L. & Williams, A.) 48–68 (Brill, 2017).

Coster, P. M. C. et al. A new early Oligocene mammal fauna from the Sirt Basin, central Libya: biostratigraphic and paleobiogeographic implications. J. Afr. Earth Sci. 104, 43–55 (2015).

Jaeger, J.-J. et al. New rodent assemblages from the Eocene Dur at-Talah escarpment (Sahara of central Libya): systematic, biochronological, and palaeobiogeographical implications. Zool. J. Linn. Soc. 160, 195–213 (2010).

Marivaux, L., Adnet, S., Benammi, M. & Tabuce, R. Anomaluroid rodents from the earliest Oligocene of Dakhla, Morocco, reveal the long-lived and morphologically conservative pattern of the Anomaluridae and Nonanomaluridae during the Tertiary in Africa. J. Syst. Palaeontol. 15, 539–569 (2017).

Marivaux, L., Adnet, S., Benammi, M., Tabuce, R. & Yans, J. Earliest Oligocene hystricognathous rodents from the Atlantic margin of northwestern Saharan Africa (Dakhla, Morocco): systematic, paleobiogeographical, and paleoenvironmental implications. J. Vertebr. Paleontol. 37, e1357567 (2017).

Seiffert, E. R. Evolution and extinction of Afro-Arabian primates near the Eocene-Oligocene boundary. Folia Primatol. 78, 314–327 (2007).

Ronquist, F. et al. A total-evidence approach to dating with fossils, applied to the early radiation of the Hymenoptera. Syst. Biol. 61, 973–999 (2012).

PubMed  PubMed Central  Article  Google Scholar 

Beck, R. M. D. & Lee, M. S. Y. Ancient dates or accelerated rates? Morphological clocks and the antiquity of placental mammals. Proc. R. Soc. B-Biol. Sci. 281, 10 (2014).

Ni, X., Flynn, J. J., Wyss, A. R. & Zhang, C. Cranial endocast of a stem platyrrhine primate and ancestral brain conditions in anthropoids. Sci. Adv. 5, eaav7913 (2019).

PubMed  PubMed Central  Article  Google Scholar 

Gunnell, G. F. et al. Fossil lemurs from Egypt and Kenya suggest an African origin for Madagascar’s aye-aye. Nat. Commun. 9, 3193 (2018).

PubMed  PubMed Central  Article  CAS  Google Scholar 

Borths, M. R., Holroyd, P. A. & Seiffert, E. R. Hyainailourine and teratodontine cranial material from the late Eocene of Egypt and the application of parsimony and Bayesian methods to the phylogeny and biogeography of Hyaenodonta (Placentalia, Mammalia). PeerJ 4, e2639 (2016).

PubMed  PubMed Central  Article  Google Scholar 

Sallam, H. & Seiffert, E. New phiomorph rodents from the latest Eocene of Egypt, and the impact of Bayesian “clock”-based phylogenetic methods on estimates of basal hystricognath relationships and biochronology. PeerJ 4, e1717 (2016).

Sallam, H. M. & Seiffert, E. R. Revision of Oligocene ‘Paraphiomys’ and an origin for crown Thryonomyoidea (Rodentia: Hystricognathi: Phiomorpha) near the Oligocene-Miocene boundary in Africa. Zool. J. Linn. Soc. 190, 352–371 (2020).

Seiffert, E. R. et al. A parapithecid stem anthropoid of African origin in the Paleogene of South America. Science 368, 194–197 (2020).

CAS  PubMed  Article  Google Scholar 

Heritage, S. et al. Ancient phylogenetic divergence of the enigmatic African rodent Zenkerella and the origin of anomalurid gliding. PeerJ 4, e2320 (2016).

PubMed  PubMed Central  Article  Google Scholar 

Fabre, P. H. et al. Flightless scaly-tailed squirrels never learned how to fly: a reappraisal of Anomaluridae phylogeny. Zoologica Scr. 47, 404–417 (2018).

Marivaux, L. & Boivin, M. Emergence of hystricognathous rodents: Palaeogene fossil record, phylogeny, dental evolution and historical biogeography. Zool. J. Linn. Soc. 187, 929–964 (2019).

Simons, E. Egyptian Oligocene primates: a review. Yrbk. Phys. Anthropol. 38, 199–238 (1995).

Hoyal Cuthill, J. F., Guttenberg, N. & Budd, G. E. Impacts of speciation and extinction measured by an evolutionary decay clock. Nature 588, 636–641 (2020).

Seiffert, E. R. et al. A fossil primate of uncertain affinities from the earliest late Eocene of Egypt. Proc. Natl. Acad. Sci. USA 107, 9712–9717 (2010).

CAS  PubMed  PubMed Central  Article  Google Scholar 

Winchester, J. M. et al. Dental topography of platyrrhines and prosimians: convergence and contrasts. Am. J. Phys. Anthropol. 153, 29–44 (2014).

Evans, A. R., Wilson, G. P., Fortelius, M. & Jernvall, J. High-level similarity of dentitions in carnivorans and rodents. Nature 445, 78–81 (2007).

CAS  PubMed  Article  Google Scholar 

Shan, S., Kovalsky, S. Z., Winchester, J. M., Boyer, D. M. & Daubechies, I. ariaDNE: a robustly implemented algorithm for Dirichlet energy of the normal. Methods Ecol. Evolution 10, 541–552 (2019).

Boyer, D. M. Relief index of second mandibular molars is a correlate of diet among prosimian primates and other euarchontan mammals. J. Hum. Evol. 55, 1118–1137 (2008).

Villéger, S., Mason, N. W. H. & Mouillot, D. New multidimensional functional diversity indices for a multifaceted framework in functional ecology. Ecology 89, 2290–2301 (2008).

Kingdon, J. Rodents, hares and rabbits. In Mammals of Africa, Vol. III (ed. Happold, D. C. D.) 602–617 (Bloomsbury Publishing, 2013).

Dinets, V. First observations on the behavior of the flightless anomalure (Zenkerella insignis). Zoology 123, 121–123 (2017).

Kirk, E. C. & Simons, E. L. Diets of fossil primates from the Fayum Depression of Egypt: a quantitative analysis of molar shearing. J. Hum. Evolution 40, 203–229 (2001).

Cramwinckel, M. J. et al. Synchronous tropical and polar temperature evolution in the Eocene. Nature 559, 382–386 (2018).

CAS  PubMed  Article  Google Scholar 

Seiffert, E. R. Chronology of Paleogene mammal localities. In Cenozoic Mammals of Africa (eds. Werdelin, L. & Sanders, W. J.) 19–26 (University of California Press, 2010).

Winkler, A. J. Neogene paleobiogeography and East African paleoenvironments: contributions from the Tugen Hills rodents and lagomorphs. J. Hum. Evol. 42, 237–256 (2002).

Gradstein, F., Ogg, J. G., Schmitz, M. & Ogg, G. The Geological Time Scale 2012 (Elsevier, 2012).

Dunn, R. E. et al. A new chronology for middle Eocene–early Miocene South American Land Mammal Ages. GSA Bull. 125, 539–555 (2013).

Zhang, Y. G., Pagani, M., Liu, Z., Bohaty, S. M. & DeConto, R. A 40-million-year history of atmospheric CO2. Philos. Trans. A Math. Phys. Eng. Sci. A 371, 20130096 (2013).

Sembroni, A., Faccenna, C., Becker, T. W., Molin, P. & Abebe, B. Longterm, deep-mantle support of the Ethiopia-Yemen Plateau. Tectonics 35, 469–488 (2016).

Prave, A. R. et al. Geology and geochronology of the Tana Basin, Ethiopia: LIP volcanism, super eruptions and Eocene–Oligocene environmental change. Earth Planet. Sci. Lett. 443, 1–8 (2016).

Wolfenden, E., Ebinger, C. J., Yirgu, G., Renne, P. & Kelley, S. P. Evolution of the southern Red Sea rift: Birth of a magmatic margin. Geol. Soc. Am. Bull. 117, 846–864 (2005).

Sepulchre, P. et al. Tectonic uplift and eastern Africa aridification. Science 313, 1419–1423 (2006).

CAS  PubMed  Article  Google Scholar 

Pik, R., Marty, B., Carignan, J. & Lavé, J. Stability of the Upper Nile drainage network (Ethiopia) deduced from (U-Th)/He thermochronometry: implications for uplift and erosion of the Afar plume dome. Earth Planet. Sci. Lett. 215, 73–88 (2003).

Ernst, R. E. & Youbi, N. How large igneous provinces affect global climate, sometimes cause mass extinctions, and represent natural markers in the geological record. Palaeogeogr. Palaeoclimatol. Palaeoecol. 478, 30–52 (2017).

Kappelman, J. et al. Oligocene mammals from Ethiopia and faunal exchange between Afro-Arabia and Eurasia. Nature 426, 549–552 (2003).

CAS  PubMed  Article  Google Scholar 

Stevens, N. et al. Palaeontological evidence for an Oligocene divergence between Old World monkeys and apes. Nature 497, 611–614 (2013).

CAS  PubMed  Article  Google Scholar 

Rasmussen, D. T. & Gutierrez, M. A mammalian fauna from the late Oligocene of northwestern Kenya. Palaeontogr. Abt. A 288, 1–52 (2009).

Sanders, W. J., Rasmussen, D. T. & Kappelman, J. New large-bodied mammals from the late Oligocene site of Chilga. Ethiop. Acta Palaeont. Polonica 49, 365–392 (2004).

Rasmussen, D. T. & Gutiérrez, M. Hyracoidea. In Cenozoic Mammals of Africa (eds Sanders, W. J. & Werdelin, L.) 123–145 (University of California Press, 2010).

Heritage, S., Rayaleh, H., Awaleh, D. G. & Rathbun, G. B. New records of a lost species and a geographic range expansion for sengis in the Horn of Africa. PeerJ 8, e9652 (2020).

PubMed  PubMed Central  Article  CAS  Google Scholar 

Patterson, B. & Upham, N. A newly recognized family from the Horn of Africa, the Heterocephalidae (Rodentia: Ctenohystrica). Zool. J. Linn. Soc. 172, 942–963 (2014).

Rasmussen, D. et al. Primitive Old World monkey from the earliest Miocene of Kenya and the evolution of cercopithecoid bilophodonty. Proc. Natl. Acad. Sci. USA 116, 6051–6056 (2019).

CAS  PubMed  PubMed Central  Article  Google Scholar 

Dembo, M. et al. The evolutionary relationships and age of Homo naledi: an assessment using dated Bayesian phylogenetic methods. J. Hum. Evolution 97, 17–26 (2016).

DeSilva, J., McNutt, E., Benoit, J. & Zipfel, B. One small step: a review of Plio-Pleistocene hominin foot evolution. Am. J. Phys. Anthropol. 168, 63–140 (2019).

Springer, M. et al. Macroevolutionary dynamics and historical biogeography of primate diversification inferred from a species supermatrix. PLoS ONE 7, e49521 (2012).

Borths, M. & Stevens, N. Simbakubwa kutokaafrika, gen. et sp. nov. (Hyainailourinae, Hyaenodonta, ‘Creodonta,’ Mammalia), a gigantic carnivore from the earliest Miocene of Kenya. J. Vertebr. Paleontol. 39, e1570222 (2019).

Springer, M. et al. Waking the undead: implications of a soft explosive model for the timing of placental mammal diversification. Mol. Phylogenet. Evol. 106, 86–102 (2017).

Paradis, E., Claude, J. & Strimmer, K. APE: analyses of phylogenetics and evolution in R language. Bioinformatics 20, 289–290 (2004).

CAS  PubMed  Article  Google Scholar 

Bapst, D. paleotree: an R package for paleontological and phylogenetic analyses of evolution. Methods Ecol. Evol. 3, 803–807 (2012).

Ronquist, F. et al. MrBayes 3.2: efficient Bayesian phylogenetic inference and model choice across a large model space. Syst. Biol. 61, 539–542 (2012).

PubMed  PubMed Central  Article  Google Scholar 

Bunn, J. et al. Comparing Dirichlet normal surface energy of tooth crowns, a new technique of molar shape quantification for dietary inference, with previous methods in isolation and in combination. Am. J. Phys. Anthropol. 145, 247–261 (2011).

Evans, A. & Jernvall, J. Patterns and constraints in carnivoran and rodent dental complexity and relative tooth size. J. Vertebr. Paleontol. 29, 92A–92A (2009).

Winchester, J. MorphoTester: an open source application for morphological topographic analysis. PLoS ONE 11, e0147649 (2016).

Allen, K., Cooke, S., Gonzales, L. & Kay, R. Dietary inference from upper and lower molar morphology in platyrrhine primates. PLoS ONE 10, e0118732 (2015).

Sokal, R. R. & Rohlf, F. J. Biometry: the Principles and Practice of Statistics in Biological Research 4th edn (WH Freeman, 2012).

Clavel, J., Escarguel, G. & Merceron, G. mvMORPH: an R package for fitting multivariate evolutionary models to morphometric data. Methods Ecol. Evol. 6, 1311–1319 (2015).

Gombin, J., Vaidyanathan, R., Agafonkin, V. Package ‘concaveman’: A Very Fast 2D Concave Hull Algorithm, v. 1.1.0 (2020).

Pebesma, E. et al. Package ‘sp’: Classes and Methods for Spatial Data. Version 1.4-5 (2021).

de Vries, D. et al. Widespread loss of mammalian lineage and dietary diversity in the early Oligocene of Afro-Arabia (Dryad,) 2021.

D.D.’s work has been supported by funds from The Leakey Foundation, the U.S. National Science Foundation (BSC-1824745), and the Natural Environment Research Council (NERC NE/T000341/1). M.R.B. was supported by NSF grant DBI-1612062. Field work in the Fayum Depression, Egypt, and digital curation of Fayum fossils has been supported by the U.S. National Science Foundation (BCS-0416164, BCS-0819186, and BCS-1231288) as well as Gordon and Ann Getty and The Leakey Foundation. Micro-CT scanning was partially supported by NSF grants DBI-1458192 and DBI-2023087, and IMLS grant MA-245704-OMS-20. K.C. Beard, E. Gheerbrant, F.K. Manthi, L. Marivaux, and M. Pickford kindly provided access to original fossils or casts. We thank Jeroen Smaers for his comments on an earlier version of this work, and for help with R code development for D.D.’s dissertation analyses. We also thank L. Werdelin, M. Orliac, and two anonymous reviewers for their helpful comments. This is Duke Lemur Center publication #1489.

Ecosystems and Environment Research Centre, School of Science, Engineering and Environment, University of Salford, Manchester, UK

Interdepartmental Doctoral Program in Anthropological Sciences, Stony Brook University, Stony Brook, NY, 11794, USA

Dorien de Vries & Steven Heritage

Duke Lemur Center Museum of Natural History, Durham, NC, 27705, USA

Steven Heritage, Matthew R. Borths, Hesham M. Sallam & Erik R. Seiffert

Mansoura University Vertebrate Paleontology, Department of Geology, Faculty of Science, Mansoura, Egypt

Institute of Global Health and Human Ecology (I-GHHE), School of Sciences and Engineering, American University in Cairo, New Cairo, Egypt

Department of Integrative Anatomical Sciences, Keck School of Medicine of USC, University of Southern California, Los Angeles, CA, USA

Department of Mammalogy, Natural History Museum of Los Angeles County, Los Angeles, CA, 90007, USA

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

D.D. collected, processed, and analyzed all dental topographic data; M.R.B., H.M.S., and E.R.S. created phylogenetic character matrices; S.H. and E.R.S. ran phylogenetic analyses; E.R.S. created the composite tree; S.H. created novel R code for, and ran, randomization, ASR, and disparity analyses; E.R.S. created figures with assistance from S.H.; E.R.S. wrote the first draft of the main text, and all authors contributed edits; all authors contributed text to the “Methods” section.

Correspondence to Erik R. Seiffert.

The authors declare no competing interests.

Peer review information Communications Biology thanks Lars Werdelin, Fabien Condamine and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Primary handling editor: Brooke LaFlamme.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit

de Vries, D., Heritage, S., Borths, M.R. et al. Widespread loss of mammalian lineage and dietary diversity in the early Oligocene of Afro-Arabia. Commun Biol 4, 1172 (2021).


Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.

Communications Biology (Commun Biol) ISSN 2399-3642 (online)

Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.