Maybe So, Maybe Not: Canis lepophagus at Hagerman Fossil Beds National Monument, Idaho, USA

A canid dentary is described from the Pliocene Glenns Ferry Formation at Hagerman Fossil Beds National Monument, south-central Idaho, USA. The specimen possesses traits in alliance with and measurements falling within or exceeding those of Canis lepophagus. The dentary, along with a tarsal IV (cuboid) and an exploded canine come from the base of the fossiliferous Sahara complex within the monument. Improved geochronologic control provided by new tephrochronologic mapping by the U.S. Geological Survey-National Park Service Hagerman Paleontology, Environments, and Tephrochronology Project supports an interpolated age of approximately 3.9 Ma, placing it in the early Blancan North American Land Mammal Age. It is conservatively referred to herein as Canis aff. C. lepophagus with the caveat that it is an early and robust example of that species. A smaller canid, initially assigned to Canis lepophagus and then to Canis ferox, is also known from Hagerman. Most specimens of Canis ferox, including the holotype, were recently reassigned to Eucyon ferox, but specimens from the Hagerman and Rexroad faunas were left as Canis sp. and possibly attributed to C. lepophagus. We agree that these smaller canids belong in Canis and not Eucyon but reject placing them within C. lepophagus; we refer to them here as Hagerman-Rexroad Canis. This study confirms the presence of two approximately coyote-sized canids at Hagerman and adds to the growing list of carnivorans now known from these fossil beds.


Canidae
Canidae include the Hesperocyoninae, Borophaginae, and Caninae. Hesperocyonines were small and omnivorous and persisted from the late Eocene through early Miocene in North America (Wang 1994). The small, omnivorous borophagines of the Oligocene diversified in the Miocene and Pliocene to include the massive-jawed genus Borophagus (Wang et al. 1999. These later "bone crushing dogs" may have had an ecological role similar to that of Old World hyaenids (Biknevicius and Ruff 1992). Borophagines spread from North America into Central America and went extinct by the early Pleistocene (Wang et al. 1999. The Caninae appeared in the early Oligocene of North America (Tedford et al. 2009) and remained small, fox-like generalists until the late Miocene, when they quickly diversified and spread into Eurasia (Montoya et al. 2009;Sotnikova and Rook 2010) and Africa (Morales et al. 2005;de Bonis et al. 2007;Geraads 2011), reaching South America by the late Pliocene (Prevosti et al. 2009). Only the Caninae (true canids) persist today.

Pliocene Caninae of North America
The basal canid genus Eucyon (Tedford and Qiu 1996;Tedford et al. 2009) dates to the Miocene late Clarendonian and Miocene-Pliocene Hemphillian Land Mammal Ages of North America (Bell et al. 2004) and late Miocene localities of Eurasia and Africa (Rook 1992(Rook , 2009Morales et al. 2005;Montoya et al. 2009;Sotnikova and Rook 2010;Werdelin et al. 2015). Eucyon has been considered a paraphyletic taxon of questionable validity, with some specimens subsequently referred to the Vulpini or Amphicyonidae (Werdelin et al. 2015). Eucyon davisi is considered basal to tribe Canini, while other members of this genus appear to align well with the wolf-like subtribe Canina (Zrzavý et al. 2018). Eucyon is generally thought to have given rise to Canis (Wang and Tedford 2007).
Eucyon and early Canis share many synapomorphies (Tedford et al. 2009) and several early species of Canis have been moved to Eucyon (Rook 1992;Spassov and Rook 2006). This includes the recent reassignment by Lucenti and Rook (2020) of Canis ferox to Eucyon ferox. This late Hemphillian to late early Blancan canid was first described by Miller and Carranza-Castañeda (1998) as intermediate in morphology to the smaller, more fox-like, and gracile Eucyon species and the typically larger, more robust Blancan Canis lepophagus. Lucenti and Rook (2020) left specimens identified as Canis ferox from the Hagerman Local Fauna of Idaho (Bjork 1970) and Rexroad Local Fauna of Kansas (Hibbard 1941) as Canis sp. and suggested its placement in Canis lepophagus or Canis sp. nov. due to a lack of key diagnostic features attributable to Eucyon.
The remaining definitive early member of Canis is the coyote-sized Canis lepophagus. Its mix of basal and derived traits places it intermediately between Eucyon and larger, later members of the genus Canis (Tedford et al. 2009). It is restricted to the (primarily late) Blancan with one early and one late early Blancan locality in Nevada and Nebraska, respectively (Tedford et al. 2009). Previous studies placed C. lepophagus as the direct or intermediary ancestor of the coyote, Canis latrans (Johnston 1938;Giles 1960;Bjork 1970;Kurtén 1974;Nowak 1979;Kurtén and Anderson 1980;Anderson 1984;Bever 2005) or wolf, Canis lupus (Kurtén 1974;Nowak 1979Nowak , 2002Nowak , 2003Kurtén and Anderson 1980). However, those extant canids share several traits absent in C. lepophagus (Tedford et al. 2009), and the Canis lineage leading to wolves and coyotes is now considered of Eurasian origin (Perri et al. 2021). Tedford et al. (2009) placed later Blancan C. lepophagus as ancestral to the Canis thooides-Canis feneus-Canis cedazoensis jackal-like lineage. Zrzavý et al. (2018) placed C. lepophagus in a basal and likely paraphyletic group of stem Canina that included Canis ferox and Canis thooides leading to Canis arnensis, but with unresolved polytomies, leaving its monospecificity in question. Systematics of this and other early canines remain largely unresolved.

Hagerman Fossil Beds National Monument (HAFO)
The Pliocene fluvial-lacustrine Glenns Ferry Formation is a stratigraphic unit that extends across southwestern Idaho. The most fossiliferous area of this formation comprises the nearly 4,300 acres of Hagerman Fossil Beds National Monument (HAFO), National Park Service, near the town of Hagerman, Idaho (Fig. 1). The Glenns Ferry dates from at least 4.2 to 3.07 Ma at Hagerman (Neville et al. 1979;Hart and Brueseke 1999;Hart unpublished data). Pleistocene sediments of Tuana Gravel (~ 1.9 Ma) and Yahoo Clay (poorly constrained, Malde 1982) occur as isolated outcrops in some areas of the monument, but these units have produced little to no fossil remains, respectively, and do not occur within the study areas discussed herein.
The Hagerman strata preserve a diverse Blancan fauna of nearly 200 species including mammals, birds, reptiles, amphibians, fishes, and invertebrates (McDonald 1993). The Carnivora are well represented with six families and, with this new canid, 19 species (Table 1). Notable carnivoran taxa include the enigmatic meline, Ferinestrix vorax (Bjork 1970), a Blancan-aged remnant of the Hemphillian ursid, Agriotherium (Samuels et al. 2009), and the oldest New World river otter, Lontra weiri (Prassack 2016). Bjork (1970) described two canids from the Hagerman fossil beds: Canis lepophagus, a canine (see also Nowak 1979) and Borophagus hilli, a borophagine. Tedford et al. (2009) reassigned most early Blancan specimens of C. lepophagus, including those from Hagerman and Rexroad (Hibbard 1941) to Canis ferox. As noted, C. ferox is now invalid with the holotype and most other specimens moved to Eucyon by Lucenti and Rook (2020). We concur with Lucenti and Rook (2020) that the smaller canids of both Hagerman and Rexroad do not belong in Eucyon and are more closely aligned with the genus Canis, but not with Canis 1 3 lepophagus. We herein refer to the Hagerman and Rexroad specimens as Hagerman-Rexroad Canis to differentiate those specimens from other canids discussed here.

Site Description
The primary specimen described herein is a robust canid dentary (HAFO 21175) collected within the Sahara complex in 2009 by a former HAFO paleontologist and a group of students and park staff. The Sahara complex is an east-facing bluff (Fig. 2) with a series of steep (50-75°) exposures primarily along its eastern and southern flanks. The complex is named for a large sandy exposure colloquially referred to by park staff as the "Sahara blowout" that extends across a large swath of the major eroded eastern knob of the bluff.
Field notes from the day HAFO 21175 was collected do not indicate the depositional context or if the specimen was found in-situ, but coordinates collected with the specimen place it within an erosive sandy zone near the base of the complex (X in Figs. 2 and 3). This is along the western edge of HAFO locality 042, a fossiliferous eight-meter exposed section of cross-bedded sand, silt, and clay. There, in-situ remains of Mammut americanum occur within a fining upward sequence of brownish gray coarse to medium lithic sand overlain by a thin lens of silty sand that contain a variety of microfauna (e.g., rodents, fish, frogs, and gastropods). Unfortunately, coordinates for HAFO 21175 are within GPS error of a wash that runs down the slope from just below the Sahara blowout to the collection site (Fig. 3). This initially prevented our confident placement of the specimen in this lower deposit.
Immediately above HAFO locality 042 are alternating angular blocky layers of beige and black clay and silty clay. Some fossils, primarily turtle shell fragments, have eroded out of these clay units, but a thick layer of overburden largely limits erosion. Continuing up the wall of the complex, an overlying covered interval of approximately 10 m is topped by an extensive bed of paper shales. These shales preserve abundant reed-like residues and carbon films and contain Fig. 1 The location of Hagerman Fossil Beds National Monument in south-central Idaho, USA. Inset shows the outline of the monument boundary along the western bank of the Snake River, west of the town of Hagerman, with a star denoting the general area of the main fossil localities discussed in this paper Table 1 List of Carnivora known from the Hagerman Fossil Beds (Bjork 1970) a Identified as Borophagus sp. by Bjork (1970) and assigned to B. hilli by McDonald (2020) b Identified as Canis lepophagus by Bjork (1970), Canis ferox by Tedford et al. (2009), Canis sp. by Lucenti and Rook (2020) and referred to here as Hagerman-Rexroad Canis c Specimens described in this paper d Werdelin (1985)  several volcanic tephra layers, including a basaltic tephra informally named bed G (Malde and Powers 1962) that is dated at 3.84 ± 0.05 Ma ( 40 Ar/ 39 Ar, recalculated) (Hart and Brueseke 1999). Overlying the paper shales, the Sahara blowout is comprised of thick beds of poorly consolidated sand interspersed with thin layers of clay. A wide range of fossil taxa have come from the eastern flank of this upper sandy unit while comparatively few fossils erode out of its steeper southern face that lies directly above HAFO locality 042. Another basaltic tephra, bed J, overlies these sands and directly underlies the 3.44 ± 0.04 Ma ( 40 Ar/ 39 Ar, recalculated) Deer Gulch lava flow (Hart and Brueseke 1999) in the northern end of the monument. If HAFO 21175 originated from these upper deposits, rather than the stratigraphic horizon at which it was collected, then its age range would be 3.84 to 3.44 Ma instead of > 3.84 Ma (Fig. 2).
The Mammut-bearing sand-silt-clay outcrop of HAFO locality 042 extends 40 m eastward along the base of the next major knob of the Sahara complex to HAFO locality 535 (Y, Fig. 2) where additional fragments of mastodon and other faunal elements have eroded out. In 2019, dry screening of sediments at HAFO locality 535 produced a canid tarsal IV (cuboid) (HAFO 26531) and an exploded canine (HAFO 26532, not described). Similar elevation, lithology, and additional surface occurrences of M. americanum support HAFO localities 042 and 535 as a synchronous depositional event. Canid material is relatively rare at Hagerman. The presence of additional canid material, from a deposit synchronous with that from where the coordinates placed HAFO 21175, strongly implies that the material belongs to a single canid. The distribution across 40 m can be explained as landscape scatter by various taphonomic agents prior to burial such as a scavenger removing the hind limb and transporting it away from the carcass (Hill 1979;Behrensmeyer 1983).

Site Geochronology
A monument-wide tephrostratigraphic study by the U.S. Geological Survey (USGS)-National Park Service (NPS) Fig. 2 The Sahara complex. Six tephra layers are directly mapped with a seventh projected from an adjacent bluff. Locations of HAFO 21175 (left, HAFO locality 042) and HAFO 26531 (right, HAFO locality 535) are marked by an X and Y, respectively. The canid icon shown in the right stratigraphic column represents the stratigraphic position of both fossil localities. The age shown for bed G is recal-culated from Hart and Brueseke (1999) using a Fish Canyon sanidine fluence monitor age of 28.201 Ma. The age shown for bed J is based on the stratigraphic position of bed J immediately below the 3.44 ± 0.04 Ma Deer Gulch lava flow in the northern portion of the monument. The age of the Deer Gulch lava flow is similarly recalculated from Hart and Brueseke (1999) Hagerman Paleontology, Environments, and Tephrochronology (PET) Project is expanding on previous work by early researchers at Hagerman (Malde and Powers 1962;Malde 1972;Neville et al. 1979;Hart and Brueseke 1999) and has mapped and characterized several volcanic tephra layers in and adjacent to the Sahara complex (Fig. 2). The location where these canid specimens were collected is approximately 25 m stratigraphically below the 3.84 ± 0.05 Ma bed G. We have further mapped and geochemically correlated bed F, another basaltic tephra, to approximately five meters stratigraphically below HAFO localities 042 and 535. The age of bed F is estimated at ca. 3.9 Ma based on previous paleomagnetic studies (Neville et al. 1979, updated to GPTS calibration of Ogg 2012. Thus, an early Blancan age of ca. 3.9 Ma is interpolated HAFO localities 042 and 535 and, therefore, HAFO 21175, HAFO 26531, and HAFO 26532. At another monument locality, an isolated lower molar (HAFO 23808) occurs higher up in the sequence, approximately 17 m above bed G at that location. An age of ca. 3.6 Ma is extrapolated for this stratigraphic interval. These data support an estimated age range of approximately 3.9 to 3.6 Ma for this larger canid at Hagerman.

Excavation Techniques
The slopes above and adjacent to HAFO locality 042 were scraped clean of overburden where possible and matrix sediments were screened with a hand-held 3.2 mm mesh sieve to determine whether those strata contained any fossil material that could have eroded out and been transported downslope to the elevation of HAFO locality 042. Sediments were also screened from several small excavation pits dug in the blocky clay and silty clay units immediately above the sandsilt-clay backwall of HAFO locality 042, along the southern face of the upper Sahara blowout, and at HAFO localities 042 and 535. The covered intervals between the lower exposures and the paper shales ( Fig. 2) were not screened as the heavy vegetation likely prevents any fossils potentially present in those layers from eroding out to the surface.

Fig. 3
The base of the Sahara complex on the day the canid dentary (HAFO 21175) was collected. Note the base of the wash to the left, where the only set of coordinates (X) used for this and other specimens collected that day was taken Table 2 Measurements of fossil canid dentaries See text for institutional abbreviations. Measurement abbreviations as follows: DDm1, dentary depth taken from the base of the m1 paraconid; DPW, dentary width taken perpendicular to the horizontal ramus along the molar arcade; DWm1, dentary width taken parallel to horizontal ramus at the point of the m1 protoconid; Sp3, shape of p3 (width/length); Sp4, shape of p4 (width/ length); Lm1, maximum anteroposterior length of m1; Wm1, maximum width of m1; Hm1, heigh of m1 taken at from the base of the paraconid; Sm1, m1 shape (Wm1/Lm1); Hm1:DD1, height of the m1 in relation to the dentary depth taken at the base of the m1 paraconid; Lm2, maximum anteroposterior length of m2; Wm2, maximum width of m2; RBL, relative blade length measured as the trigonid length in relation to the total length of the m1); and RGA , relative grinding area measured as the square root of the summed area for the talonid and m2 divided by the trigonid length; Sm2, m2 shape (Wm2/Lm2). Mean and standard deviation provided for taxonomic groups with an NISP > 2

Tephrochronology
Tephra deposits throughout the monument were mapped, described, and sampled. Samples were then laboratory processed, petrographically characterized, and analyzed for major and minor elemental concentrations via electron microprobe (EMP) and then geochemically correlated utilizing standard tephrochronology methods (Sarna-Wojcicki 1976, 2000Westgate and Gorton 1981;Sarna-Wojcicki and Davis 1991;Lowe 2011). An initial tephrochronological framework was developed and progressively refined using supplementary field mapping data collected over several years. Volcanic glass compositions acquired by EMP analysis provided geochronologic constraints by geochemically correlating our new samples to previously identified samples in the USGS Tephrochronology Project's reference databases and to published analytical results of tephra samples collected and analyzed by Hart and Brueseke (1999) and Malde and Powers (1962). All tephra names are informal and were established by Malde and Powers (1962) and subsequent authors, except for bed F2 and the basilisk ash layer, both of which were identified for the first time by the USGS-NPS PET Project. These tephra layers are also informally named in keeping with the historic nomenclature put forth by Malde and Powers (1962).
Dates of tephra layers were recalculated from the 40 Ar/ 39 Ar data first reported by Hart and Brueseke (1999) relative to the astronomically tuned Fish Canyon sanidine (FCs) age of 28.201 Ma (Kuiper et al. 2008 using decay constants of Min et al. 2000) and are quoted with a 2-sigma error, although the original dates were reported with a 1-sigma error (written communication, WK Hart). Bed G was directly dated via 40 Ar/ 39 Ar methods (3.84 ± 0.05 Ma recalculated, Hart and Brueseke 1999). The estimated age of bed J is based on its occurrence immediately below the 3.44 ± 0.04 Ma Deer Gulch lava flow in the northern portion of the monument. Therefore, the inferred age of bed J is ≥ 3.44 ± 0.04 Ma (Hart and Brueseke 1999). Paleomagnetism data support an age of ca. 3.9 Ma for bed F based on its position elsewhere in the monument where it lies approximately 120 m stratigraphically above sediments that preserve the Cochiti Reversed-Polarity Subchron of the Gilbert  (Neville et al. 1979;Hart and Brueseke 1999).

Dental Nomenclature and Measurements
See Online Resource 1 for the list of comparative materials. Dental nomenclature follows Tedford et al. (2009). Measurements (Table 2) follow Van Valkenburgh and Wayne (1994) and Tedford et al. (2009). Measurements were acquired using an electronic digital caliper (accuracy to 0.01 mm) with serial data output capacity and rounded to the nearest 0.01 mm. Each measurement was repeated three times with the mean used for each specimen.
Abbreviations are as follows: DDm1, dentary depth taken from the base of the m1 paraconid; DPW, dentary width taken perpendicular to the horizontal ramus along the molar arcade; DWm1, dentary width taken parallel to the horizontal ramus at the point of the m1 protoconid; Sp3, shape of p3 (width/length); Sp4, shape of p4 (width/length); Lm1, maximum anteroposterior length of m1; Wm1, maximum width of m1; Hm1, height of m1 taken from the base of the paraconid; Sm1, m1 shape (Wm1/Lm1); Hm1:DD1, height of the m1 in relation to the dentary depth taken at the base medial, c. distal, and d. anterior view. Abbreviations as follows: cal, calcaneal facet; cun, cuneiform facet; nav, navicular facet; ten, tendineal sulcus and groove for tendon of fibularis; tub, plantar tuberance; IV, facet for digit IV; V, facet for digit V. Scale bar equals 10 mm of the m1 paraconid; Lm2, maximum anteroposterior length of m2; Wm2, maximum width of m2; RBL, relative blade length measured as the trigonid length in relation to the total length of the m1; RGA , relative grinding area measured as the square root of the summed area for the talonid and m2 divided by the trigonid length; and Sm2, m2 shape (Wm2/ Lm2).
We used a variance-covariance matrix that disregards groups in Principal Component Analysis (PCA) to maximize variance across groups. We did not standardize data because it was commensurable. Substitution of mean (mean imputation) was used to address missing data when measurements for a given variable were available for at least 50% of a given group. Estimating missing data by substitution of mean can lead to a decrease in variance and correlativity but does not affect the overall species grouping patterns in small sample sets (Strauss et al. 2003). A lack of suitable measurements across groups led to the removal of RGA, Sm2, SP3, and SP4 from further analysis. Preliminary test runs found the loadings for those measurements had little effect on the PCA. Four PCA analyses were conducted: all taxonomic groups, larger canids, smaller canids, and for the m1 to allow for inclusion of HAFO 23808. Principle component eigenvalues greater than 1.0 and character loading values of 0.5 or greater are considered. We did not perform a MANOVA due to the small sample sizes. Site locality Hagerman Fossil Beds National Monument, Twin Falls County, Idaho. Fossil location coordinates are protected by federal law (Paleontological Resource Protection Act, 16 U.S.C. 470aaa-aaa-11). Exact location information will be shared with qualified researchers on request through the NPS Research Permit and Reporting System.
Age Geological age range of specimens approximately 3.9-3.6 Ma, Early Blancan North American Land Mammal Age.
Diagnosis HAFO 21175 (Fig. 4a-c) is a left dentary broken anterior to the p3 with p3-m3 intact and the base of the masseteric fossa preserved. Wear is mild on the p3-4 and m1 trigonid such that no features are lost and the maximum heights of cuspulids and conids are measurable. Wear is moderate on the m1 talonid and m2, and moderately heavy on the m3 leaving the margins of some features difficult to distinguish. This coyote-sized canid is confidently assigned to Canis and further referred to Canis aff. C. lepophagus based on the following morphological traits assigned to C. lepophagus by Tedford et al. (2009): deep masseteric fossa flanked by a strong anterior margin and thick masseteric line; relatively deep horizontal ramus; robust, posteriorly expanded premolars with large, bulbous principle accessory cuspulids and a small secondary accessory cuspulid; p4 crown superior to m1 paraconid; robust m1; entoconulid connecting the metaconid and entoconid; hypoconid and entoconid worn but connected by a weak transverse cristid with the hypoconid broader and the entoconid slightly displaced posteriorly; small hypoconulid presented as a second basin in the talonid; notable anterior expansion of the m2; retention of two trigonid cusps (metaconid and protoconid); m2 expanded anteriorly and narrows posteriorly; m3 ovoid.

HAFO 21175
The mandibular corpus is deep and wide through to the masseteric fossa and exceeds measurements, particularly below and posterior to the carnassial (DDm1, DPW, and DWm1, Table 2), of the largest specimens of Canis lepophagus. The masseteric fossa is shallow and flanked by robust anterior and ventral margins. Tedford et al. (2009) note two mental foramina below the p1-p2 diastema in C. lepophagus; the ramus of HAFO 21175 is broken anterior to the p3 but it lacks the foramen ventral to the p3 that is seen in Eucyon. The premolars are uncrowded with slight diastemata following the p2, p3, and p4, as are observed in most specimens of C. lepophagus. Eucyon davisi, E. ferox, and Hagerman/Rexroad Canis exhibit notable diastemata posterior to the p1-p3, but not between the p4-m1. Tooth crowding can vary intraspecifically (Ameen et al. 2017) and may not be a reliable diagnostic trait.
HAFO 21175 has bulbous and posteriorly expanded premolars (Sp3 and Sp4, Table 2) as in Canis arnensis (Lucenti and Rook 2016) that differ from the elongate and bilaterally tapered premolars of E. davisi, E. ferox, Hagerman-Rexroad Canis, and C. lepophagus. The only other early Blancan specimen of C. lepophagus (F:AM 49295) lacks mandibular material for comparison, but upper premolars of F:AM 49295 are elongate as seen in other specimens of C. lepophagus ( fig. 41 in Tedford et al. 2009). In Canis thooides, the premolars are also robust, but less elongate and taper posteriorly, while in Canis edwardii the p4 is reduced relative to the m1 paraconid.
The principal cusp of the p4 extends just above the m1 paraconid and is labially expanded but posteriorly constrained. The p4 primary accessory cuspulid is worn but lobed and more offset labially compared to the accessory cuspulid of the p3. In Eucyon, including E. ferox, a large secondary accessory cuspulid is distinct from the cingulid (Lucenti and Rook 2020). An ancestral trait, the secondary accessory cuspulid is typically lost in later Canis, usually with the primary accessory cuspulid projecting posteriorly to meet an enlarged cingulid ridge, but it is retained in some later taxa, including Canis (Xenocyon) lycaonoides and some specimens of Canis lupus (Tedford et al. 2009) and Canis anthus (USNM 476034). We found the secondary accessory cuspulid in E. davisi (F:AM 49294 and F:AM 63009B) and E. ferox (F:AM 63035, 63060, and 63061) confirms as observed by Lucenti and Rook (2020). In Hagerman-Rexroad Canis and C. edwardii, a greatly reduced and trenchant secondary accessory cuspulid emerges from the base separate from the posterior accessory cuspulid. It is typically reduced in early C. lepophagus (see also Tedford et al. 2009), and in a few specimens (CWT 2617 and LACM 1246) the primary accessory cuspulid extends posteriorly to the cingulid such that the secondary accessory cuspulid is lost. In HAFO 21175, C. thooides (F:AM 63092), and Canis aureus thai (USNM 236632) the secondary accessory cuspulid is reduced but retains prominence and emerges off the medial edge of the primary accessory cusp. Canis arnensis also exhibits two distinct accessory cuspulids (Lucenti and Rook 2016).
The m1 is short in length and height relative to dentary depth compared to most specimens of C. lepophagus but with greater breadth expansion across the protoconid creating a broad carnassial (Wm1, Hm1:DD1, and Sm1, Table 2). Much of this length is attributed to the trigonid, with a moderately trenchant talonid basin. In C. lepophagus the m1 paraconid's anterior face is typically nearly vertical, but it can also present itself as slightly oblique as seen in WTUC 560 (fig. 42 e in Tedford et al. 2009), LACM 1343, USNM 26116, and in HAFO 21175. The metaconid and entoconid are connected by an entoconulid. Severe wear along the labial side of the talonid prevents confirmation of similar cristids between the hypoconid and protoconid. The hypoconid and entoconid are worn, but the wear facet for the hypoconid is twice that of the entoconid with the latter displaced more posterior. Small cristids occur along the cingulid, just posterior to the hypoconid and entoconid, but the cingulid itself is poorly developed. These cristids produce a posteriorly opened hypoconulid shelf that is mesially constrained and separated from the talonid basin by a weak transverse crest. This crest is used by Tedford et al. (2009) to distinguish C. lepophagus from the more primitive Eucyon (and Hagerman-Rexroad Canis) where the talonid basin lies open and continuous with the hypoconulid shelf (see also Lucenti and Rook 2020). In Hagerman-Rexroad Canis the hypoconid and entoconid merge, separating the hypoconulid shelf from the talonid basin, as seen in later canids. It is deeper and more triangular but also open due to a poorly defined cingulid ridge, compared to HAFO 21175.
The m2 is sub-ovoid in occlusal view with expansion across the protoconid and metaconid, labial contraction just posterior to the protoconid, and expansion across the hypoconid and entoconid reduced only relative to the trigonid. The m2 exhibits heavy, uneven wear with the protoconid worn to the incipient anterolabial cingulid. The metaconid is also worn and only slightly offset posterior to where the protoconid is presumed to have been pre-wear.
The m3 is ovoid with slight posterior constriction and is heavily worn anteriorly with an unworn hypoconulid shelf.  Table 3 for summary and loadings  Table 4 for summary and loadings

HAFO 23808
HAFO 23808 (Fig. 4d-f) is an unworn, isolated right m1. The specimen aligns with C. lepophagus and not Eucyon or Hagerman-Rexroad Canis based on the following: low series of cristids connects the metaconid and entoconid and the protoconid and hypoconid; strong transverse crest connects the hypoconid and entoconid, closing the talonid basin posterior-labially; entoconid lies posterior to the hypoconid; hypoconulid shelf with a small cusp creates a second distinct basin to the talonid; anterior face of paraconid semi-oblique.
Unlike HAFO 21175, the transverse crest is a distinct series of cristids with a blade-like edge to both the anterior and posterior labial edge of the protoconid. The molar length is shorter than seen in C. lepophagus, but other molar tooth measurements (Wm1, Hm1, and Sm1, Table 2) align it with C. lepophagus. It differs from HAFO 21175 in being longer, taller, and narrower.

HAFO 26531
HAFO 26531 (Fig. 5) is a left tarsal IV (cuboid) assigned to Canis based on comparison to modern Canis latrans. It falls within the size range of C. latrans but is more gracile. The navicular facet is poorly preserved, but the cuneiform facet extends further antero-laterally than in C. latrans (Fig. 5a). On the plantar face, a shortened neck and greatly reduced plantar tuberosity produce an almost shear face in lateral view (Fig. 5b). This differs from the long neck, bulbous tuberosity, and gradual slope in C. latrans. In dorsal view, the anterior face is of approximate width, but the tendineal sulcus is reduced (Fig. 5c). The anterior facet for digit V is slightly damaged but appears planar while the anterior facet for digit IV is concave and narrows ventrally, as in C. latrans, but with a stronger medial curve (Fig. 5d). Anteriorly, the calcaneal facet is narrower and more rectangular than squared (Fig. 5e). There was no tarsal IV specimen of C. lepophagus for comparison, but its Canis-like attributes and physical location within nearby geologically contemporaneous sediments, relative to HAFO 21175, leaves us confident that it belongs to the same individual as HAFO 21175.

Taxonomic Assignment
We assign HAFO 21175 to Canis aff. C. lepophagus but not without acknowledging some notable differences. Features in the lower dentition best align it with C. lepophagus, but the depth and thickness of the mandibular corpus, particularly below and posterior to the carnassial, exceed that of this mostly late Blancan species. The only other earlier Blancan specimen of C. lepophagus (F:AM 49295) lacks a preserved dentary for direct comparison, but cranial measurements place that specimen within the range for later C. lepophagus (Tedford et al. 2009). Principal Component Analysis supports the distinctive morphospace of HAFO 21175 compared to other canid groups (Figs. 6,7,and 9), while HAFO 23808 falls within C. lepophagus (Fig. 9, Tables 2 and 6). Measurements for SP3 and SP4 were not used in the PCA, but their shape further distinguishes HAFO 21175 from C. lepophagus.
Mandibular thickening below the m1 and m2 in HAFO 21175 points to greater force and either increased durophagy and/or reliance on larger prey (Van Valkenburgh 1991;Wroe et al. 2005). Wider and rounder premolars also support greater durophagy (Van Valkenburgh 1991;Van Valkenburgh et al. 2004) and further differentiates it from C. lepophagus. Heavy wear along the grinding surface (RGA) may lend additional support for hard part consumption and/or be an indicator of age. The enlarged premolars and reduced grinding area, coupled with a more robust mandibular corpus, could suggest a robust male (Van Valkenburgh and Sacco 2002;Tedford et al. 2009;O'Keefe et al. 2013;Brannick et al. 2015), but HAFO 21175 would still be an anomalously robust jawed, early specimen of C. lepophagus. It may instead represent an early branching event, potentially leaving C. lepophagus polyphyletic. Given the limited material and high variability seen in this taxon, the conservative approach is to refer it to C. aff. C. lepophagus. HAFO 23808 better aligns with C. lepophagus, but as an isolated m1 is also left as C. aff. C. lepophagus.

Discussion
The study of extant canid relationships, particularly among the wolf-and jackal-like canids, incorporates morphological, behavioral, and genetic data, yet our understanding of those canids remains in flux (Nowak and Federof 2002;Schwartz and Vucetich 2009;Rueness et al. 2011;Gaubert et al. 2012;Viranta et al. 2017;Waples et al. 2018;Joshi et al. 2020;Stoyanov 2020). Recently diverged sister taxa can overlap substantially morphologically, and species introgression is common today (Fredrickson and Hedrick 2006;Galov et al. 2015;Gopalakrishnan et al. 2018) and likely in the past. Allometric and phenotypic plasticity in extant canids (Machado and Teta 2020), particularly of facial traits, can leave larger or younger conspecifics as outliers and potentially lead to misidentification as distinct species (Wayne 1986;Prevosti et al. 2013;Gopalakrishnan et al. 2018;Machado et al. 2018;Chemisquy et al. 2019;Machado and Teta 2020 Sansalone et al. 2015;Flower 2016). These factors have led to difficulties in distinguishing extant canid species (Prevosti et al. 2013;von Holdt et al. 2016;Chemisquy et al. 2019). Not surprisingly, similar issues arise in the study of fossil canid accumulations, but with even fewer types of data to aid us in reconstructing canid relationships. In the late Pleistocene, larger accumulations of canids have revealed clinal shifts in canid body size and morphology (Leonard et al. 2007;Slater et al. 2009;Meachen and Samuels 2012;O'Keefe et al. 2013;Sansalone et al. 2015;Flower 2016) and changes to in-situ populations over time due to changes in climate (Brannick et al. 2015) and competition (Meachen et al. 2014), but even then, phylogenetic relationships remain in flux (Perri et al. 2021). Pre-Pleistocene North American canid samples, such as those from Hagerman, are currently of insufficient size for population variance studies.
The limited material and often-ambiguous morphological traits of Canis leaves the taxonomic status of Hagerman's canids as tentative, with the larger assigned to the group of canids currently referred to Canis lepophagus and a more gracile species known from both Hagerman and Rexroad that is best left, for now, as Canis sp. Whether Canis lepophagus gave rise to the jackal-like lineage of C. thooides and, possibly, to C. arnensis is also unresolved but possible given some overall similarities in dental traits.
If Canis lepophagus remain monospecific, then this species was widespread across much of what is now the southwestern, central, and northwestern United States, and long-lived with a known range of 4.5 Ma (Panaca Formation, Nevada, Lindsay et al. 2002) to 2.5 Ma (Grandview Fauna, Idaho, Repenning et al. 1995). This same land area today supports several contemporaneous species (and subspecies) of Canis including taxa that are only genetically identifiable as taxonomically distinct (Wilson et al. 2000;Schwartz and Vucetich 2009;Heppenheimer et al. 2018;Machado and Teta 2020).

Conclusion
Canid systematics is in a state of flux as new discoveries and a broader range of applied techniques continue to change our understanding of extant and extinct canid relationships. For extant canids, the application of genetic, isotopic, and ecological data has deemed some morphologically distinct populations conspecific while other morphologically indistinct populations were found to differ significantly in diet, behavior, and/or genetics, thus leading to taxonomic reassignment. For extinct canids, especially in the North American Pliocene, a mix of synapomorphies and apparent intraspecific variance coupled with small, incomplete sample sizes have led to multiple morphological-based taxonomic reassignments of specimens at both the species and genus level.
The canids at Hagerman Fossil Beds National Monument highlight the ambiguities of canid morphological traits in assigning taxonomic affinity. Hagerman's smaller canid was published as Canis lepophagus, Canis ferox, and most recently relegated to Canis sp., with C. ferox no longer considered a valid taxon. Newly described Hagerman specimens of a larger coyote-sized canid morphologically align best with and are assigned to Canis aff C. lepophagus, but with an unusually robust dentary that adds to the variation already evident in that species. Canis lepophagus spans two million years and much of the continental United States and may represent multiple speciation events or phenotypic variance related to geographic and ecological differences that additional fossil material or the application of new techniques may help to shed light on.