Dire wolves were the last of an ancient New World canid lineage.

Dire wolves are considered to be one of the most common and widespread large carnivores in Pleistocene America1, yet relatively little is known about their evolution or extinction. Here, to reconstruct the evolutionary history of dire wolves, we sequenced five genomes from sub-fossil remains dating from 13,000 to more than 50,000 years ago. Our results indicate that although they were similar morphologically to the extant grey wolf, dire wolves were a highly divergent lineage that split from living canids around 5.7 million years ago. In contrast to numerous examples of hybridization across Canidae2,3, there is no evidence for gene flow between dire wolves and either North American grey wolves or coyotes. This suggests that dire wolves evolved in isolation from the Pleistocene ancestors of these species. Our results also support an early New World origin of dire wolves, while the ancestors of grey wolves, coyotes and dholes evolved in Eurasia and colonized North America only relatively recently.

On Information Rank Deficiency in Phenotypic Covariance Matrices.

This article investigates a form of rank deficiency in phenotypic covariance matrices derived from geometric morphometric data, and its impact on measures of phenotypic integration. We first define a type of rank deficiency based on information theory then demonstrate that this deficiency impairs the performance of phenotypic integration metrics in a model system. Lastly, we propose methods to treat for this information rank deficiency. Our first goal is to establish how the rank of a typical geometric morphometric covariance matrix relates to the information entropy of its eigenvalue spectrum. This requires clear definitions of matrix rank, of which we define three: the full matrix rank (equal to the number of input variables), the mathematical rank (the number of nonzero eigenvalues), and the information rank or "effective rank" (equal to the number of nonredundant eigenvalues). We demonstrate that effective rank deficiency arises from a combination of methodological factors-Generalized Procrustes analysis, use of the correlation matrix, and insufficient sample size-as well as phenotypic covariance. Secondly, we use dire wolf jaws to document how differences in effective rank deficiency bias two metrics used to measure phenotypic integration. The eigenvalue variance characterizes the integration change incorrectly, and the standardized generalized variance lacks the sensitivity needed to detect subtle changes in integration. Both metrics are impacted by the inclusion of many small, but nonzero, eigenvalues arising from a lack of information in the covariance matrix, a problem that usually becomes more pronounced as the number of landmarks increases. We propose a new metric for phenotypic integration that combines the standardized generalized variance with information entropy. This metric is equivalent to the standardized generalized variance but calculated only from those eigenvalues that carry nonredundant information. It is the standardized generalized variance scaled to the effective rank of the eigenvalue spectrum. We demonstrate that this metric successfully detects the shift of integration in our dire wolf sample. Our third goal is to generalize the new metric to compare data sets with different sample sizes and numbers of variables. We develop a standardization for matrix information based on data permutation then demonstrate that Smilodon jaws are more integrated than dire wolf jaws. Finally, we describe how our information entropy-based measure allows phenotypic integration to be compared in dense semilandmark data sets without bias, allowing characterization of the information content of any given shape, a quantity we term "latent dispersion". [Canis dirus; Dire wolf; effective dispersion; effective rank; geometric morphometrics; information entropy; latent dispersion; modularity and integration; phenotypic integration; relative dispersion.].

Lions and brown bears colonized North America in multiple synchronous waves of dispersal across the Bering Land Bridge.

The Bering Land Bridge connecting North America and Eurasia was periodically exposed and inundated by oscillating sea levels during the Pleistocene glacial cycles. This land connection allowed the intermittent dispersal of animals, including humans, between Western Beringia (far northeast Asia) and Eastern Beringia (northwest North America), changing the faunal community composition of both continents. The Pleistocene glacial cycles also had profound impacts on temperature, precipitation and vegetation, impacting faunal community structure and demography. While these palaeoenvironmental impacts have been studied in many large herbivores from Beringia (e.g., bison, mammoths, horses), the Pleistocene population dynamics of the diverse guild of carnivorans present in the region are less well understood, due to their lower abundances. In this study, we analyse mitochondrial genome data from ancient brown bears (Ursus arctos; n = 103) and lions (Panthera spp.; n = 39), two megafaunal carnivorans that dispersed into North America during the Pleistocene. Our results reveal striking synchronicity in the population dynamics of Beringian lions and brown bears, with multiple waves of dispersal across the Bering Land Bridge coinciding with glacial periods of low sea levels, as well as synchronous local extinctions in Eastern Beringia during Marine Isotope Stage 3. The evolutionary histories of these two taxa underline the crucial biogeographical role of the Bering Land Bridge in the distribution, turnover and maintenance of megafaunal populations in North America.

Muscular anatomy of the forelimb of tiger (Panthera tigris).

Dissection reports of large cats (family Felidae) have been published since the late 19th century. These reports generally describe the findings in words, show drawings of the dissection, and usually include some masses of muscles, but often neglect to provide muscle maps showing the precise location of bony origins and insertions. Although these early reports can be highly useful, the absence of visual depictions of muscle attachment sites makes it difficult to compare muscle origins and insertions in living taxa and especially to reconstruct muscle attachments in fossil taxa. Recently, more muscle maps have been published in the primary literature, but those for large cats are still limited. Here, we describe the muscular anatomy of the forelimb of the tiger (Panthera tigris), and compare muscle origins, insertions, and relative muscle masses to other felids to identify differences that may reflect functional adaptations. Our results reiterate the conservative nature of felid anatomy across body sizes and behavioral categories. We find that pantherines have relatively smaller shoulder muscle masses, and relatively larger muscles of the caudal brachium, pronators, and supinators than felines. The muscular anatomy of the tiger shows several modifications that may reflect an adaptation to terrestrial locomotion and a preference for large prey. These include in general a relatively large m. supraspinatus (shoulder flexion), an expanded origin for m. triceps brachii caput longum, and relatively large m. triceps brachii caput laterale (elbow extension), as well as relatively large mm. brachioradialis, abductor digiti I longus, and abductor digiti V. Muscle groups that are well developed in scansorial taxa are not well developed in the tiger, including muscles of the cranial compartment of the brachium and antebrachium, and m. anconeus. Overall, the musculature of the tiger strongly resembles that of the lion (Panthera leo), another large-bodied terrestrial large-prey specialist.

Pleistocene origins, western ghost lineages, and the emerging phylogeographic history of the red wolf and coyote.

The red wolf (Canis rufus) of the eastern US was driven to near-extinction by colonial-era persecution and habitat conversion, which facilitated coyote (C. latrans) range expansion and widespread hybridization with red wolves. The observation of some grey wolf (C. lupus) ancestry within red wolves sparked controversy over whether it was historically a subspecies of grey wolf with its predominant "coyote-like" ancestry obtained from post-colonial coyote hybridization (2-species hypothesis) versus a distinct species closely related to the coyote that hybridized with grey wolf (3-species hypothesis). We analysed mitogenomes sourced from before the 20th century bottleneck and coyote invasion, along with hundreds of modern amplicons, which led us to reject the 2-species model and to investigate a broader phylogeographic 3-species model suggested by the fossil record. Our findings broadly support this model, in which red wolves ranged the width of the American continent prior to arrival of the grey wolf to the mid-continent 60-80 ka; red wolves subsequently disappeared from the mid-continent, relegated to California and the eastern forests, which ushered in emergence of the coyote in their place (50-30 ka); by the early Holocene (12-10 ka), coyotes had expanded into California, where they admixed with and phenotypically replaced western red wolves in a process analogous to the 20th century coyote invasion of the eastern forests. Findings indicate that the red wolf pre-dated not only European colonization but human, and possibly coyote, presence in North America. These findings highlight the urgency of expanding conservation efforts for the red wolf.

Regional variation in the cheetah (Acinonyx jubatus) revisited: Morphology of wild and captive populations.

The cheetah (Acinonyx jubatus) is listed as a vulnerable species by the International union for the conservation of nature (IUCN), including two critically endangered subspecies, the Saharan cheetah, and the Iranian cheetah, so it is imperative that we understand variation in cheetah morphology to make good decisions regarding the conservation of this species. Here, we aim to determine whether northeastern African cheetahs have smaller body sizes than southern African cheetahs. This study also adds to our knowledge of cheetah morphology from two cheetah populations that do not yet have comprehensive published data: Kenya, and northeastern Africa, including captive individuals. We calculated means and standard deviations on cranial and body measurements of live or in few cases, freshly dead, cheetahs from the aforementioned populations, plus previously published data on Namibian and Botswanan cheetahs and compared them to one another using multivariate analysis of variance. Results show that northeastern African cheetahs have smaller body sizes than southern and eastern African populations. We also found that captive cheetahs retain the morphological characteristics of their ancestral population- captive cheetahs from southern Africa have similar body sizes to wild southern African cheetahs and larger body sizes than captives from northeastern Africa. Other analyses regarding cheetah growth agree with previous studies on Namibian and Botswanan cheetah populations rates. As such, this study can serve as a baseline for the care of captive cheetah populations to maintain healthy weights and body proportions.

Postcranial diversity and recent ecomorphic impoverishment of North American gray wolves.

Recent advances in genomics and palaeontology have begun to unravel the complex evolutionary history of the gray wolf, Canis lupus Still, much of their phenotypic variation across time and space remains to be documented. We examined the limb morphology of the fossil and modern North American gray wolves from the late Quaternary (

The fossil record of phenotypic integration and modularity: A deep-time perspective on developmental and evolutionary dynamics.

Variation is the raw material for natural selection, but the factors shaping variation are still poorly understood. Genetic and developmental interactions can direct variation, but there has been little synthesis of these effects with the extrinsic factors that can shape biodiversity over large scales. The study of phenotypic integration and modularity has the capacity to unify these aspects of evolutionary study by estimating genetic and developmental interactions through the quantitative analysis of morphology, allowing for combined assessment of intrinsic and extrinsic effects. Data from the fossil record in particular are central to our understanding of phenotypic integration and modularity because they provide the only information on deep-time developmental and evolutionary dynamics, including trends in trait relationships and their role in shaping organismal diversity. Here, we demonstrate the important perspective on phenotypic integration provided by the fossil record with a study of Smilodon fatalis (saber-toothed cats) and Canis dirus (dire wolves). We quantified temporal trends in size, variance, phenotypic integration, and direct developmental integration (fluctuating asymmetry) through 27,000 y of Late Pleistocene climate change. Both S. fatalis and C. dirus showed a gradual decrease in magnitude of phenotypic integration and an increase in variance and the correlation between fluctuating asymmetry and overall integration through time, suggesting that developmental integration mediated morphological response to environmental change in the later populations of these species. These results are consistent with experimental studies and represent, to our knowledge, the first deep-time validation of the importance of developmental integration in stabilizing morphological evolution through periods of environmental change.

Evolution in coyotes (Canis latrans) in response to the megafaunal extinctions.

Living coyotes modify their behavior in the presence of larger carnivores, such as wolves. However, little is known about the effects of competitor presence or absence on morphological change in coyotes or wolves over long periods of time. We examined the evolution of coyotes and wolves through time from the late Pleistocene, during which many large carnivorous species coexisted as predators and competitors, to the Recent; this allowed us to investigate evolutionary changes in these species in response to climate change and megafaunal extinctions at the end of the Pleistocene. We measured postcranial skeletal morphologies of wolves (Canis lupus) and coyotes (C. latrans) from Pleistocene-aged tar deposits, as well as early, mid, and recent Holocene populations of both. We found few morphological differences between Pleistocene and Holocene wolf populations. Conversely, we found many differences in coyotes: Pleistocene coyotes were larger and more robust than Holocene populations. However, within 1,000 y of the megafaunal extinctions, coyotes are morphologically indistinguishable from modern populations. We cannot attribute these differences directly to climate change because modern coyotes do not follow Bergmann's rule, which states body size increases with decreasing temperature. Instead, we suggest that Pleistocene coyotes may have been larger and more robust in response to larger competitors and a larger-bodied prey base. Although we cannot separate competition from predator-prey interactions, this study indicates that the effects of biotic interactions can be detected in the fossil record.

A novel multiscale assessment of community assembly across time, space, and functional niche.

A basic ecological tenet is that organisms in a community occupy different niches and have different traits, but how consistently competition, selection, and phylogenetic effects structure communities remains uncertain. Are all communities created equal? We examine how mammalian carnivoran communities are assembled with regard to mass, diet, and locomotion. Here, we use a multivariate nearest-neighbor framework to examine multiple North American localities spanning 3 million years to determine whether community assembly is consistent through time and four modern localities around the world to assess the effects of habitat. Additionally, we examined how trait patterns differ among families and how family-level evolutionary effects affect them. We found some broadly consistent patterns, although differences are more pronounced than similarities. Diet is more affected by evolutionary constraints than by time or place. Locomotion is most affected by habitat, and the ability to partition niches is related to habitat heterogeneity. Mass is influenced by family, but also by habitat and the mass-selective extinction events at the end-Pleistocene. These findings indicate that assembly patterns are not largely determined by within-community interactions but instead show that each community is a product of its independent variables.

Pre-Younger Dryas megafaunal extirpation at Rancho La Brea linked to fire-driven state shift.

The cause, or causes, of the Pleistocene megafaunal extinctions have been difficult to establish, in part because poor spatiotemporal resolution in the fossil record hinders alignment of species disappearances with archeological and environmental data. We obtained 172 new radiocarbon dates on megafauna from Rancho La Brea in California spanning 15.6 to 10.0 thousand calendar years before present (ka). Seven species of extinct megafauna disappeared by 12.9 ka, before the onset of the Younger Dryas. Comparison with high-resolution regional datasets revealed that these disappearances coincided with an ecological state shift that followed aridification and vegetation changes during the Bølling-Allerød (14.69 to 12.89 ka). Time-series modeling implicates large-scale fires as the primary cause of the extirpations, and the catalyst of this state shift may have been mounting human impacts in a drying, warming, and increasingly fire-prone ecosystem.

A mummified Pleistocene gray wolf pup.

In July 2016, a mummified carcass of an ancient wolf (Canis lupus) pup (specimen YG 648.1) was discovered in thawing permafrost in the Klondike goldfields, near Dawson City, Yukon, Canada (Figure 1A). The wolf pup mummy was recovered along a small tributary of Last Chance Creek during hydraulic thawing that exposed the permafrost sediment in which it was preserved. This mummified wolf pup is important to the local Tr'ondëk Hwëch'in people, who named it Zhùr, meaning 'wolf' in the Hän language of their community. Here, we report detailed morphometric, isotopic, and genetic analyses of Zhùr that reveal details of her appearance, evolutionary relationships to other wolves and short life-history and ecology. Zhùr is the most complete wolf mummy known. She lived approximately 57,000 years ago and died in her den during a collapse of the sediments. During her short life, she ate aquatic resources, and is related to ancient Beringian and Russian gray wolves and her clade is basal to all living gray wolves. VIDEO ABSTRACT.

Locomotor correlates of the scapholunar of living and extinct carnivorans.

The relationship of carpal morphology to ecology and habitat is under studied in carnivorans and more generally in mammals. Here, we use 3D-scanning techniques to assess the usefulness of a carpal bone, the scapholunar, in carnivorans to reflect ecology and habitat, and to reconstruct the ecology of five extinct carnivorans from two fossil sites: Rancho La Brea and Natural Trap Cave. We 3D-scanned scapholunars and measured articular surface areas and angles between articular facets using GeoMagic and Rhino 3D-software. We analyzed the difference in these metrics using multivariate analysis of variance and discriminant function analysis. Results show that the scapholunar reflects ecological signal, with clear groupings of cursorial carnivorans and grappling/climbing carnivorans; however, phylogenetic signal was also present in the results with hyaenids, canids, and large felids in distinct morphospaces. Extinct species Miracinonyx trumani (American cheetah) and Smilodon fatalis (sabertooth cat) showed surprising results with M. trumani grouping with pantherines instead of Acinonyx or Puma, suggesting it runs but still retains the ability to grapple prey. S. fatalis groups with pantherines, but also shows some unique adaptations, suggesting it had a different range of wrist motion than living cats. Overall, the scapholunar is a good indicator of ecology and functional morphology and can be another tool to use in modern and fossil carnivorans to reconstruct extinct ecologies and locomotor behaviors.

Correction: Low-cost cross-taxon enrichment of mitochondrial DNA using in-house synthesised RNA probes.

[This corrects the article DOI: 10.1371/journal.pone.0209499.].

Low-cost cross-taxon enrichment of mitochondrial DNA using in-house synthesised RNA probes.

Hybridization capture with in-solution oligonucleotide probes has quickly become the preferred method for enriching specific DNA loci from degraded or ancient samples prior to high-throughput sequencing (HTS). Several companies synthesize sets of probes for in-solution hybridization capture, but these commercial reagents are usually expensive. Methods for economical in-house probe synthesis have been described, but they do not directly address one of the major advantages of commercially synthesised probes: that probe sequences matching many species can be synthesised in parallel and pooled. The ability to make "phylogenetically diverse" probes increases the cost-effectiveness of commercial probe sets, as they can be used across multiple projects (or for projects involving multiple species). However, it is labour-intensive to replicate this with in-house methods, as template molecules must first be generated for each species of interest. While it has been observed that probes can be used to enrich for phylogenetically distant targets, the ability of this effect to compensate for the lack of phylogenetically diverse probes in in-house synthesised probe sets has not been tested. In this study, we present a refined protocol for in-house RNA probe synthesis and evaluated the ability of probes generated using this method from a single species to successfully enrich for the target locus in phylogenetically distant species. We demonstrated that probes synthesized using long-range PCR products from a placental mammal mitochondrion (Bison spp.) could be used to enrich for mitochondrial DNA in birds and marsupials (but not plants). Importantly, our results were obtained for approximately a third of the cost of similar commercially available reagents.

Ancient DNA: Saber-Toothed Cats Are the Same Beasts After All.

Ancient DNA from the saber-toothed cat Homotherium reveals that the late Pleistocene species from Europe and North America were the same. Homotherium turns out to be only distantly related to the well-known saber-toothed Smilodon.

Extinct Beringian wolf morphotype found in the continental U.S. has implications for wolf migration and evolution.

Pleistocene diversity was much higher than today, for example there were three distinct wolf morphotypes (dire, gray, Beringian) in North America versus one today (gray). Previous fossil evidence suggested that these three groups overlapped ecologically, but split the landscape geographically. The Natural Trap Cave (NTC) fossil site in Wyoming, USA is an ideally placed late Pleistocene site to study the geographical movement of species from northern to middle North America before, during, and after the last glacial maximum. Until now, it has been unclear what type of wolf was present at NTC. We analyzed morphometrics of three wolf groups (dire, extant North American gray, Alaskan Beringian) to determine which wolves were present at NTC and what this indicates about wolf diversity and migration in Pleistocene North America. Results show NTC wolves group with Alaskan Beringian wolves. This provides the first morphological evidence for Beringian wolves in mid-continental North America. Their location at NTC and their radiocarbon ages suggest that they followed a temporary channel through the glaciers. Results suggest high levels of competition and diversity in Pleistocene North American wolves. The presence of mid-continental Beringian morphotypes adds important data for untangling the history of immigration and evolution of Canis in North America.

Measures of Relative Dentary Strength in Rancho La Brea Smilodon fatalis over Time.

The late Pleistocene megafaunal extinction of approximately 12,000 years ago, included the demise of Smilodon fatalis, a hypercarnivore from the Rancho La Brea deposits, which has been studied across time by looking at different deposits or pits to determine morphological size and shape changes and trends during this time. To better understand functional aspects of these changes, this study focused on a measure of jaw strength over time, which can give an indication of morphological changes within the jaw that cannot be seen using surface morphometrics. By radiographing dentaries, cortical bone can be seen, which provides an estimate of resistance to bending forces while biting, and can be measured and used as an indicator of jaw strength. Measurements were taken at repeatable locations on the dentary of the depth of the cortical bone, and of a standardized measure of cortical bone, which allows for the comparison between different individuals. Specimens included those of five different pits ranging from about 37 Kybp to 13 Kybp (just before the extinction of S. fatalis). No significant difference was found in the depth of jaws at any of the measurement points from any of the pits. However, significant differences were found in both the actual thickness of cortical bone, and the standardized thickness of cortical bone at the lower P4 between pit 13 (which had the lowest amount of bone) and pit 61/67 (which had the highest). These conclusions support other studies that have shown that individuals in pit 13 were under physiological and perhaps dietary stress, which may be reflected in the deposition of cortical bone, while the opposite trend is seen in the individuals in pit 61/67. Our results further support findings suggesting Smilodon did not appear to be morphologically most vulnerable right before its extinction.

Ecological changes in Coyotes (Canis latrans) in response to the ice age megafaunal extinctions.

Coyotes (Canis latrans) are an important species in human-inhabited areas. They control pests and are the apex predators in many ecosystems. Because of their importance it is imperative to understand how environmental change will affect this species. The end of the Pleistocene Ice Age brought with it many ecological changes for coyotes and here we statistically determine the changes that occurred in coyotes, when these changes occurred, and what the ecological consequences were of these changes. We examined the mandibles of three coyote populations: Pleistocene Rancho La Brean (13-29 Ka), earliest Holocene Rancho La Brean (8-10 Ka), and Recent from North America, using 2D geometric morphometrics to determine the morphological differences among them. Our results show that these three populations were morphologically distinct. The Pleistocene coyotes had an overall robust mandible with an increased shearing arcade and a decreased grinding arcade, adapted for carnivory and killing larger prey; whereas the modern populations show a gracile morphology with a tendency toward omnivory or grinding. The earliest Holocene populations are intermediate in morphology and smallest in size. These findings indicate that a niche shift occurred in coyotes at the Pleistocene/Holocene boundary - from a hunter of large prey to a small prey/more omnivorous animal. Species interactions between Canis were the most likely cause of this transition. This study shows that the Pleistocene extinction event affected species that did not go extinct as well as those that did.