Bayesian Total-Evidence Dating Revisits Sloth Phylogeny and Biogeography: A Cautionary Tale on Morphological Clock Analyses

Combining morphological and molecular characters through Bayesian total-evidence dating allows inferring the phylogenetic and timescale framework of both extant and fossil taxa, while accounting for the stochasticity and incompleteness of the fossil record. Such an integrative approach is particularly needed when dealing with clades such as sloths (Mammalia: Folivora), for which developmental and biomechanical studies have shown high levels of morphological convergence whereas molecular data can only account for a limited percentage of their total species richness. Here, we propose an alternative hypothesis of sloth evolution that emphasizes the pervasiveness of morphological convergence and the importance of considering the fossil record and an adequate taxon sampling in both phylogenetic and biogeographic inferences. Regardless of different clock models and morphological datasets, the extant sloth Bradypus is consistently recovered as a megatherioid, and Choloepus as a mylodontoid, in agreement with molecular-only analyses. The recently extinct Caribbean sloths (Megalocnoidea) are found to be a monophyletic sister-clade of Megatherioidea, in contrast to previous phylogenetic hypotheses. Our results contradict previous morphological analyses and further support the polyphyly of "Megalonychidae," whose members were found in five different clades. Regardless of taxon sampling and clock models, the Caribbean colonization of sloths is compatible with the exhumation of islands along Aves Ridge and its geological time frame. Overall, our total-evidence analysis illustrates the difficulty of positioning highly incomplete fossils, although a robust phylogenetic framework was recovered by an a posteriori removal of taxa with high percentages of missing characters. Elimination of these taxa improved topological resolution by reducing polytomies and increasing node support. However, it introduced a systematic and geographic bias because most of these incomplete specimens are from northern South America. This is evident in biogeographic reconstructions, which suggest Patagonia as the area of origin of many clades when taxa are underrepresented, but Amazonia and/or Central and Southern Andes when all taxa are included. More generally, our analyses demonstrate the instability of topology and divergence time estimates when using different morphological datasets and clock models and thus caution against making macroevolutionary inferences when node support is weak or when uncertainties in the fossil record are not considered.

Bayesian tip-dated phylogeny and biogeography of Cissampelideae (Menispermaceae): Mitigating the effects of homoplastic morphological characters

The integration of morphological and molecular data is essential to understand the affinities of fossil taxa and spatio-temporal evolutionary processes of organisms. However, homoplastic morphological characters can mislead the placement of fossil taxa and impact downstream analyses. Here, we provide an example of how to mitigate effectively the effect of morphological homoplasy on the placement of fossil taxa and biogeographic inferences of Cissampelideae. We assembled three data types, morphological data only, morphological data with a molecular scaffold and combined morphological and molecular data. By removing high-level homoplastic morphological data or reweighting the morphological characters, we conducted 15 parsimony, 12 undated Bayesian and four dated Bayesian analyses. Our results show that the 14 selected Cissampelideae fossil taxa are placed poorly when based only on morphological data, but the addition of molecular scaffold and combination of morphological and molecular data greatly improve the resolution of fossil nodes. We raise the monotypic Stephania subg. Botryodiscia to generic status and discover that three fossils previously assigned to Stephania should be members of Diploclisia. The Bayesian tip-dated tree recovered by removing homoplastic morphological characters with a Rescaled Consistency Index <0.25 has the highest stratigraphic fit and consequently generates more reasonable biogeographic reconstruction for Cissampelideae. Cissampelideae began to diversify in Asia in the latest Cretaceous and subsequently dispersed to South America around the Cretaceous-Palaeogene boundary. Two dispersal events from Asia to Africa occurred in the Early Eocene and the Late Eocene-Late Oligocene, respectively. These findings provide guidelines and practical methods for mitigating the effects of homoplastic morphological characters on fossil placements and Bayesian tip-dating, as well as insights into the past tropical floristic exchanges among different continents.

A macroevolutionary common-garden experiment reveals differentially evolvable bone organization levels in slow arboreal mammals

Eco-morphological convergence, i.e., similar phenotypes evolved in ecologically convergent taxa, naturally reproduces a common-garden experiment since it allows researchers to keep ecological factors constant, studying intrinsic evolutionary drivers. The latter may result in differential evolvability that, among individual anatomical parts, causes mosaic evolution. Reconstructing the evolutionary morphology of the humerus and femur of slow arboreal mammals, we addressed mosaicism at different bone anatomical spatial scales. We compared convergence strength, using it as indicator of evolvability, between bone external shape and inner structure, with the former expected to be less evolvable and less involved in convergent evolution, due to anatomical constraints. We identify several convergent inner structural traits, while external shape only loosely follows this trend, and we find confirmation for our assumption in measures of convergence magnitude. We suggest that future macroevolutionary reconstructions based on bone morphology should include structural traits to better detect ecological effects on vertebrate diversification.

A timescale for placental mammal diversification based on Bayesian modeling of the fossil record

The timing of the placental mammal radiation has been the focus of debate over the efficacy of competing methods for establishing evolutionary timescales. Molecular clock analyses estimate that placental mammals originated before the Cretaceous-Paleogene (K-Pg) mass extinction, anywhere from the Late Cretaceous to the Jurassic. However, the absence of definitive fossils of placentals before the K-Pg boundary is compatible with a post-Cretaceous origin. Nevertheless, lineage divergence must occur before it can be manifest phenotypically in descendent lineages. This, combined with the non-uniformity of the rock and fossil records, requires the fossil record to be interpreted rather than read literally. To achieve this, we introduce an extended Bayesian Brownian bridge model that estimates the age of origination and, where applicable, extinction through a probabilistic interpretation of the fossil record. The model estimates the origination of placentals in the Late Cretaceous, with ordinal crown groups originating at or after the K-Pg boundary. The results reduce the plausible interval for placental mammal origination to the younger range of molecular clock estimates. Our findings support both the Long Fuse and Soft Explosive models of placental mammal diversification, indicating that the placentals originated shortly prior to the K-Pg mass extinction. The origination of many modern mammal lineages overlapped with and followed the K-Pg mass extinction.

Cranial osteology of a juvenile specimen of Acratocnus ye (Mammalia, Xenarthra, Folivora) and its ontogenetic and phylogenetic implications

The present study comprises a description of the skull and jaw anatomy of a juvenile specimen of the Antillean sloth Acratocnus ye, from the Holocene of Haiti. Detailed descriptions and illustrations are provided of the skull bones and their sutural connections, which normally fuse in adults. Descriptions are also provided for the mandible and ear ossicles, as well as endocranial surfaces and sinuses exposed by breaks. The anatomy of our juvenile A. ye is compared to that of adult A. ye to assess ontogenetic changes in the skull. Several of these ontogenetic features are significant new observations that impact the relationships within Xenarthra as a whole, or between Xenarthrans and other placental mammals, most notably, the presence of a separate mesethmoid element, the presence of alveoli for a lower deciduous canine and anterior incisor, and the presence of separate rostral and caudal entotympanic elements. A full list of such changes are provided. In addition, the specimen provides information on phylogenetically relevant characters, including features unique to the genus Acratocnus, and features of the clade Choloepodini, including Acratocnus, the smaller extinct Antillean sloth Neocnus, and the extant two-toed sloth Choloepus. Contrary to previous studies, Acratocnus shares as many features with Choloepus as it does with its fellow Antillean form Neocnus in the present study, which is consistent with current morphology-based phylogenetic hypotheses regarding the relationships within Choloepodini. The current study highlights the need for further anatomical and phylogenetic investigations of Antillean sloths (Megalocnidae/Megalonychidae), and juvenile sloths in general.

Taxonomic revision of maned sloths, subgenus Bradypus (Scaeopus), Pilosa, Bradypodidae, with revalidation of Bradypus crinitus Gray, 1850

We present a taxonomic revision of maned sloths, subgenus Bradypus (Scaeopus), a taxon endemic to the Brazilian Atlantic Forest and currently composed of a single species, the vulnerable Bradypus torquatus. Our review is based on coalescent species delimitation analyses using mitochondrial and nuclear DNA, morphological analyses, and field observations. Our integrative approach demonstrates that two species of maned sloth can be recognized: the northern maned sloth (Bradypus torquatusIlliger, 1811) occurring in the Brazilian states of Bahia and Sergipe, and the southern maned sloth (Bradypus crinitusGray, 1850), occurring in Rio de Janeiro and Espirito Santo states. The two species diverged in the Early Pliocene and are allopatrically distributed. We discuss the biogeographic pattern of the two maned sloth species, comparing it with other Atlantic Forest mammals. We also suggest that the conservation status of both maned sloths needs to be reassessed after this taxonomic rearrangement.

Revisiting proboscidean phylogeny and evolution through total evidence and palaeogenetic analyses including Notiomastodon ancient DNA

The extinct Gomphotheriidae is the only proboscidean family that colonized South America. The phylogenetic position of the endemic taxa has been through several revisions using morphological comparisons. Morphological studies are enhanced by paleogenetic analyses, a powerful tool to resolve phylogenetic relationships; however, ancient DNA (aDNA) preservation decreases in warmer regions. Despite the poor preservation conditions for aDNA in humid, sub-tropical climates, we recovered ∼3,000 bp of mtDNA of Notiomastodon platensis from the Arroyo del Vizcaíno site, Uruguay. Our calibrated phylogeny places Notiomastodon as a sister taxon to Elephantidae, with a divergence time of ∼13.5 Ma. Additionally, a total evidence analysis combining morphological and paleogenetic data shows that the three most diverse clades within Proboscidea diverged during the early Miocene, coinciding with the formation of a land passage between Africa and Eurasia. Our results are a further step toward aDNA analyses on Pleistocene samples from subtropical regions and provide a framework for proboscidean evolution.

Integrative Approach Uncovers New Patterns of Ecomorphological Convergence in Slow Arboreal Xenarthrans

Identifying ecomorphological convergence examples is a central focus in evolutionary biology. In xenarthrans, slow arboreality independently arose at least three times, in the two genera of ‘tree sloths’, Bradypus and Choloepus, and the silky anteater, Cyclopes. This specialized locomotor ecology is expectedly reflected by distinctive morpho-functional convergences. Cyclopes, although sharing several ecological features with ‘tree sloths’, do not fully mirror the latter in their outstandingly similar suspensory slow arboreal locomotion. We hypothesized that the morphology of Cyclopes is closer to ‘tree sloths’ than to anteaters, but yet distinct, entailing that slow arboreal xenarthrans evolved through ‘incomplete’ convergence. In a multivariate trait space, slow arboreal xenarthrans are hence expected to depart from their sister taxa evolving toward the same area, but not showing extensive phenotypical overlap, due to the distinct position of Cyclopes. Conversely, a pattern of ‘complete’ convergence (i.e., widely overlapping morphologies) is hypothesized for ‘tree sloths’. Through phylogenetic comparative methods, we quantified humeral and femoral convergence in slow arboreal xenarthrans, including a sample of extant and extinct non-slow arboreal xenarthrans. Through 3D geometric morphometrics, cross-sectional properties (CSP) and trabecular architecture, we integratively quantified external shape, diaphyseal anatomy and internal epiphyseal structure. Several traits converged in slow arboreal xenarthrans, especially those pertaining to CSP. Phylomorphospaces and quantitative convergence analyses substantiated the expected patterns of ‘incomplete’ and ‘complete’ convergence for slow arboreal xenarthrans and ‘tree sloths’, respectively. This work, highlighting previously unidentified convergence patterns, emphasizes the value of an integrative multi-pronged quantitative approach to cope with complex mechanisms underlying ecomorphological convergence.

Ossification pattern of the unusual pisiform in two-toed (Choloepus) and three-toed sloths (Bradypus)

Two-toed (Choloepus sp.) and three-toed (Bradypus sp.) sloths possess short, rounded pisiforms that are rare among mammals and differ from other members of Xenarthra like the giant anteater (Myrmecophaga tridactyla) which retain elongated, rod-like pisiforms in common with most mammals. Using photographs, radiographs, and μCT, we assessed ossification patterns in the pisiform and the paralogous tarsal, the calcaneus, for two-toed sloths, three-toed sloths, and giant anteaters to determine the process by which pisiform reduction occurs in sloths and compare it to other previously studied examples of pisiform reduction in humans and orangutans. Both extant sloth genera achieve pisiform reduction through the loss of a secondary ossification center and the likely disruption of the associated growth plate based on an unusually porous subchondral surface. This represents a third unique mechanism of pisiform reduction among mammals, along with primary ossification center loss in humans and retention of two ossification centers with likely reduced growth periods in orangutans. Given the remarkable similarities between two-toed and three-toed sloth pisiform ossification patterns and the presence of pisiform reduction in fossil sloths, extant sloth pisiform morphology does not appear to represent a recent convergent adaptation to suspensory locomotion, but instead is likely to be an ancestral trait of Folivora that emerged early in the radiation of extant and fossil sloths.

Dasycerine rove beetles: Cretaceous diversification, phylogeny and historical biogeography (Coleoptera: Staphylinidae: Dasycerinae)

Within the hyperdiverse beetle family Staphylinidae, Dasycerinae is one of the smallest and most cryptic subfamilies, comprising a sole extant genus characterized by a latridiid beetle-like body form. Little has been known about their early diversification, character evolution, phylogeny and historical biogeography because of limited fossil material and lack of a phylogeny integrating extant and extinct representatives. Here we report an unexpectedly diverse dasycerine fauna from the mid-Cretaceous of northern Myanmar, including a new genus and four new species. To reconstruct the early evolutionary history of Dasycerinae, we present a phylogenetic framework of the subfamily based on a dataset integrating all extant and extinct taxa using parsimony, maximum-likelihood and Bayesian methods. †Cedasyrus gen. n., characterized by distinct sexual dimorphism in antennal and elytral lengths, is recovered as the basal-most lineage, sister to the remaining two extinct genera and all living Dasycerus species. †Vetudasycerus is recovered as sister to †Protodasycerus + Dasycerus. Among all extinct taxa, †Protodasycerus bears distinctly longer elytra, and appears to represent a transitional form from †Vetudasycerus to Dasycerus. Phylogenetic inferences and ancestral distribution reconstruction support an "Out-of-Orient" model for Dasycerinae. Either the Bering- or North Atlantic Land Bridge may have served as a passageway for dasycerine dispersal between Eurasian and North American continents. An elevation-reconstruction analysis indicated that the ancestor of the extant Dasycerus probably lived at a high altitude and stayed at this elevation through the end of the Miocene. We propose that the extinction of dasycerine ancestors living on the Tethyan islands at low altitude was likely caused by sea-level rise and climatic warming during the Late Cretaceous. The high-altitude areas might have played the role of refugia that harboured subalpine derivatives which eventually gave rise to the extant Dasycerus.

Bayesian Morphological Clock versus Parsimony: An insight into the relationships and dispersal events of postvacuum Cricetidae (Rodentia, Mammalia)

Establishing an evolutionary timeline is fundamental for tackling a great variety of topics in evolutionary biology, including the reconstruction of patterns of historical biogeography, coevolution and diversification. However, the tree of life is pruned by extinction and molecular data cannot be gathered for extinct lineages. Until recently methodological challenges have prevented the application of tip-dating Bayesian approaches in morphology-based fossil-only datasets. Herein, we present a morphological dataset for a group of cricetid rodents to which we apply an array of methods fairly new in palaeontology that can be used by palaeontologists for the analysis of entirely extinct clades. We compare the tree topologies obtained by traditional parsimony, time-calibrated and non-calibrated Bayesian inference phylogenetic approaches and calculate stratigraphic congruence indices for each. Bayesian tip-dated clock methods outperform parsimony in the case of our dataset, which includes highly homoplastic morphological characters. Regardless, all three topologies support the monophyly of Megacricetodontinae, Democricetodontinae and Cricetodontinae. Dispersal and speciation events inferred through Bayesian Binary Markov chain Monte Carlo and biodiversity analyses provide evidence for a correlation between biogeographic events, climatic changes and diversification in cricetids.

Phylogenomics and Historical Biogeography of Seahorses, Dragonets, Goatfishes, and Allies (Teleostei: Syngnatharia): Assessing Factors Driving Uncertainty in Biogeographic Inferences

The charismatic trumpetfishes, goatfishes, dragonets, flying gurnards, seahorses, and pipefishes encompass a recently defined yet extraordinarily diverse clade of percomorph fishes-the series Syngnatharia. This group is widely distributed in tropical and warm-temperate regions, with a great proportion of its extant diversity occurring in the Indo-Pacific. Because most syngnatharians feature long-range dispersal capabilities, tracing their biogeographic origins is challenging. Here, we applied an integrative phylogenomic approach to elucidate the evolutionary biogeography of syngnatharians. We built upon a recently published phylogenomic study that examined ultraconserved elements by adding 62 species (total 169 species) and one family (Draconettidae), to cover ca. 25% of the species diversity and all 10 families in the group. We inferred a set of time-calibrated trees and conducted ancestral range estimations. We also examined the sensitivity of these analyses to phylogenetic uncertainty (estimated from multiple genomic subsets), area delimitation, and biogeographic models that include or exclude the jump-dispersal parameter (j). Of the three factors examined, we found that the j parameter has the strongest effect in ancestral range estimates, followed by number of areas defined, and tree topology and divergence times. After accounting for these uncertainties, our results reveal that syngnatharians originated in the ancient Tethys Sea ca. 87 Ma (84-94 Ma; Late Cretaceous) and subsequently occupied the Indo-Pacific. Throughout syngnatharian history, multiple independent lineages colonized the eastern Pacific (6-8 times) and the Atlantic (6-14 times) from their center of origin, with most events taking place following an east-to-west route prior to the closure of the Tethys Seaway ca. 12-18 Ma. Ultimately, our study highlights the importance of accounting for different factors generating uncertainty in macroevolutionary and biogeographic inferences.

Evolution of bone cortical compactness in slow arboreal mammals

Convergent evolution is a major topic in evolutionary biology. Low bone cortical compactness (CC, a measure of porosity of cortical bone) in the extant genera of "tree sloths," has been linked to their convergent slow arboreal ecology. This proposed relationship of low CC with a slow arboreal lifestyle suggests potential convergent evolution of this trait in other slow arboreal mammals. Femoral and humeral CC were analyzed in "tree sloths," lorisids, koala, and extinct palaeopropithecids and Megaladapis, in comparison to closely related but ecologically distinct taxa, in a phylogenetic framework. Low CC in "tree sloths" is unparalleled by any analyzed clade and the high CC in extinct sloths suggests the recent convergence of low CC in "tree sloths." A tendency for low CC was found in Palaeopropithecus and Megaladapis. However, lorisids and the koala yielded unexpected CC patterns, preventing the recognition of a straightforward convergence of low CC in slow arboreal mammals. This study uncovers a complex relationship between CC and convergent evolution of slow arboreality, highlighting the multifactorial specificity of bone microstructure.

Taxon incompleteness and discrete time bins affect character change rates in simulated data

Estimating how fast or slow morphology evolves through time (phenotypic change rate, PR) has become common in macroevolutionary studies and has been important for clarifying key evolutionary events. However, the inclusion of incompletely scored taxa (e.g. fossils) and variable lengths of discrete arbitrary time bins could affect PR estimates and potentially mask real PR patterns. Here, the impact of taxon incompleteness (unscored data) on PR estimates is assessed in simulated data. Three different time bin series were likewise evaluated: bins evenly spanning the tree length (i), a shorter middle bin and longer first and third bins (ii), and a longer middle bin and shorter first and third bins (iii). The results indicate that PR values decrease as taxon incompleteness increases. Statistically significant PR values, and the dispersion among PR values, depended on the time bins. These outcomes imply that taxon incompleteness can undermine our capacity to infer morphology evolutionary dynamics and that these estimates are also influenced by our choice of discrete time bins. More importantly, the present results stress the need for a better approach to deal with taxon incompleteness and arbitrary discrete time bins.

Gregariousness in the giant sloth Lestodon (Xenarthra): multi-proxy approach of a bonebed from the Last Maximum Glacial of Argentine Pampas

Megamammals constituted an important component in the Pleistocene faunal communities of South America. Paleobiological and paleoecological studies involving different megamammal taxa have increased significantly in the last years, but there are still several poorly-known issues of its life history. In this work, we analyze an assemblage composed of 13 individuals of different ontogenetic stages, and possibly different sex, belonging to the giant ground sloth Lestodon armatus (Xenarthra, Folivora), recovered from Playa del Barco site (Pampean Region, Argentina). A dating of 19,849 years Cal BP allows assigning this assemblage to a period of the MIS (Marine Isotope Stage) 2 related to the end of the Last Glacial Maximum. Based on multiple lines of research (e.g. taphonomy, paleopathology, osteohistology, isotopy), we interpret the origin of the assemblage and diverse paleobiological and paleoecological aspects (e.g. social behavior, ontogenetic changes, sexual dimorphism, diseases, resource and habitat use, trophic relationships) of L. armatus. Evidence suggests that the assemblage was formed by a local single event of catastrophic mortality, which affected different members of a social group. This record represents the first accurate evidence of gregariousness for this ground sloth, providing new data on a poorly-known behavior among extinct Folivora.

A novel method for quantifying femoral neck anteversion: A case study in extinct and extant sloths

Extinct sloths represent a wide range of morphological, locomotor, and body size variation. Researchers have examined femoral neck angle in two dimensions to hypothesize locomotor behaviors in this group; however, this measure does not account for femoral neck anteversion. Here, we present a new method for quantifying femoral neck anteversion angle, in addition to femoral neck angle, to capture the 3D position of the femoral head/neck. Femora of extant (n = 21; Bradypus and Choloepus) and extinct (n = 49; Acratocnus, Megalocnus, Neocnus, and Parocnus) sloths were surface scanned and their surface models used to calculate three angles of femoral neck anteversion and femoral neck angle. Femoral neck anteversion was calculated as the angle between the femoral neck axis and the geometric axis of the femoral condyles (GA), the 35% cross section axis, and a trochanter axis. Femoral neck angle was calculated as the angle between the femoral neck and shaft axes. Genera were compared using ANOVAs with post hoc multiple comparisons for each angle. Femoral neck angle and femoral neck anteversion relative to the cross section were also analyzed. Significant differences among genera exist for all angles, (p < .001) but not all angles separate all genera. Femoral neck and anteversion angles typically yield different results, demonstrating the utility of analyzing both angles. The GA and cross section angles are highly correlated in sloths, with the exception of comparisons among Megalocnus, Parocnus, and Neocnus, suggesting morphological variation in the distal femur. While this method was applied to sloths, it has broad applicability to mammalian groups.

Palaeoproteomics resolves sloth relationships

The living tree sloths Choloepus and Bradypus are the only remaining members of Folivora, a major xenarthran radiation that occupied a wide range of habitats in many parts of the western hemisphere during the Cenozoic, including both continents and the West Indies. Ancient DNA evidence has played only a minor role in folivoran systematics, as most sloths lived in places not conducive to genomic preservation. Here we utilize collagen sequence information, both separately and in combination with published mitochondrial DNA evidence, to assess the relationships of tree sloths and their extinct relatives. Results from phylogenetic analysis of these datasets differ substantially from morphology-based concepts: Choloepus groups with Mylodontidae, not Megalonychidae; Bradypus and Megalonyx pair together as megatherioids, while monophyletic Antillean sloths may be sister to all other folivorans. Divergence estimates are consistent with fossil evidence for mid-Cenozoic presence of sloths in the West Indies and an early Miocene radiation in South America.

Ancient Mitogenomes Reveal the Evolutionary History and Biogeography of Sloths

Living sloths represent two distinct lineages of small-sized mammals that independently evolved arboreality from terrestrial ancestors. The six extant species are the survivors of an evolutionary radiation marked by the extinction of large terrestrial forms at the end of the Quaternary. Until now, sloth evolutionary history has mainly been reconstructed from phylogenetic analyses of morphological characters. Here, we used ancient DNA methods to successfully sequence 10 extinct sloth mitogenomes encompassing all major lineages. This includes the iconic continental ground sloths Megatherium, Megalonyx, Mylodon, and Nothrotheriops and the smaller endemic Caribbean sloths Parocnus and Acratocnus. Phylogenetic analyses identify eight distinct lineages grouped in three well-supported clades, whose interrelationships are markedly incongruent with the currently accepted morphological topology. We show that recently extinct Caribbean sloths have a single origin but comprise two highly divergent lineages that are not directly related to living two-fingered sloths, which instead group with Mylodon. Moreover, living three-fingered sloths do not represent the sister group to all other sloths but are nested within a clade of extinct ground sloths including Megatherium, Megalonyx, and Nothrotheriops. Molecular dating also reveals that the eight newly recognized sloth families all originated between 36 and 28 million years ago (mya). The early divergence of recently extinct Caribbean sloths around 35 mya is consistent with the debated GAARlandia hypothesis postulating the existence at that time of a biogeographic connection between northern South America and the Greater Antilles. This new molecular phylogeny has major implications for reinterpreting sloth morphological evolution, biogeography, and diversification history.