
==== Front
Sci Rep
Sci Rep
Scientific Reports
2045-2322
Nature Publishing Group UK London

39294146
65729
10.1038/s41598-024-65729-3
Article
Phenotypic diversity in early Australian dingoes revealed by traditional and 3D geometric morphometric analysis
Koungoulos Loukas G. loukas.koungoulos@anu.edu.au

127
Hulme-Beaman Ardern 34
Fillios Melanie 5
Willandra Lakes Region World Heritage Aboriginal Advisory GroupKelly Mick 6
Slade Joan 6
Ellis Coral 6
Reyland Maureen 6
Winch Patricia 6
Kelly Daniel 6
Clark Warren 6
Mitchell Ernest 6
Johnston Ivan 6

1 https://ror.org/0384j8v12 grid.1013.3 0000 0004 1936 834X Department of Archaeology, School of Humanities, The University of Sydney, Sydney, Australia
2 https://ror.org/019wvm592 grid.1001.0 0000 0001 2180 7477 Archaeology and Natural History, College of Asia and the Pacific, School of Culture, History and Language, The Australian National University, Canberra, Australia
3 https://ror.org/04xs57h96 grid.10025.36 0000 0004 1936 8470 Department of Veterinary Anatomy, Physiology and Pathology, Institute of Infection, Veterinary and Ecological Sciences, University of Liverpool, Liverpool, UK
4 https://ror.org/04xs57h96 grid.10025.36 0000 0004 1936 8470 Department of Archaeology, Classics and Egyptology, School of Histories, Languages and Cultures, University of Liverpool, Liverpool, UK
5 https://ror.org/04r659a56 grid.1020.3 0000 0004 1936 7371 Department of Archaeology, School of Humanities, Arts and Social Sciences, The University of New England, Armidale, Australia
6 Department of Planning and Environment, Buronga, NSW Australia
7 https://ror.org/02zv4ka60 grid.438303.f 0000 0004 0470 8815 Australian Museum Research Institute, Australian Museum, Sydney, Australia
18 9 2024
18 9 2024
2024
14 2122822 12 2023
24 6 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by-nc-nd/4.0/ Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.
The dingo is a wild dog endemic to Australia with enigmatic origins. Dingoes are one of two remaining unadmixed populations of an early East Asian dog lineage, the other being wild dogs from the New Guinea highlands, but morphological connections between these canid groups have long proved elusive. Here, we investigate this issue through a morphometric study of ancient dingo remains found at Lake Mungo and Lake Milkengay, in western New South Wales. Direct accelerated mass spectrometry (AMS) radiocarbon dates from an ancient Lake Mungo dingo demonstrate that dingoes with a considerably smaller build than the predominant modern morphotype were present in semi-arid southeastern Australia c.3000–3300 calBP. 3D geometric morphometric analysis of a near-complete Mungo cranium finds closest links to East Asian and New Guinean dogs, providing the first morphological evidence of links between early dingoes and their northern relatives. This ancient type is no longer extant within the range of modern dingo variability, but populations from nearby southeastern Australia show a closer resemblance than those to the north and west. Our results reaffirm prior characterisations of regional variability in dingo phenotype as not exclusively derived from recent domestic dog hybridisation but as having an earlier precedent, and suggest further that the dingo’s phenotype has changed over time.

Subject terms

Archaeology
Palaeontology
Animal physiology
http://dx.doi.org/10.13039/501100001774 University of Sydney Australian Government Research Training Program Koungoulos Loukas G. issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

The dingo is a canid endemic to Australia, with enigmatic origins and a heavily-debated specific taxonomic identity. Dingoes are closer genetic relatives of East Asian dogs than wolves or other wild canids, but simultaneously exhibit a number of behavioural and phenotypic traits comparable to wild canids that are absent or otherwise unusual in domestic dogs1,2. Dingoes are “dogs” which maintain stable wild populations entirely independent of reliance on human-derived resources and/or recruitment from domestic populations3, a uniquely naturalised status shared only by closely-related dogs found in the highlands of New Guinea4,5. As such, dingoes are considered by some to be amongst the earliest-established naturalised populations of any formerly domestic mammal, and certainly the oldest known example for Canis familiaris6.

Molecular dating estimates of the divergence of dingoes from New Guinean and Asian dogs suggests they arrived in Australia prior to 5000BP7, but the earliest directly-dated skeletal remains are less than 3500 years old, following closely the earliest remains of dogs in Indonesia and Near Oceania8. Nevertheless, genetic assessments of both ancient and modern canids widely agree that dingoes represent a surviving branch of the earliest southern foray of C. familiaris into the region of Australasia9. As such, knowledge of the dingo’s origins and relationship to non-Australian canids is of key importance to understanding the early history of dogs in the Asia–Pacific region, and by extension, the movement and relations of its’ resident human cultures10.

Before the era of wide-scale genetic and genomic applications to understanding the phylogenetic relations of dogs and their relatives, morphometric approaches struggled to find close matches to the dingo amongst the native dogs of the nearby Asia–Pacific region. Gollan11, in an examination of both modern and ancient samples from a wide area spanning from South Asia to Polynesia, suggested that dingoes were most closely and probably directly derived from Indian pariah dogs, rejecting any connection with Southeast Asian dogs. Later, Corbett12 suggested that Thai pariah dogs formed an intermediary stage between the early dingo and modern domestic dog. More recently, Gonzalez13, argued that dingoes derived from a particular type of pariah dog still extant in small numbers in India, which subsequently drifted towards an ancestral form (Canis lupus pallipes) during isolation and natural selection in Australia. These morphometric studies all rejected any morphological relation of dingo to the wild dogs found in the New Guinea highlands, on the basis of major differences in size and overall form between the two.

Conversely, genetic and genomic research unanimously rejects all of the above suppositions, clearly indicating that the dingo’s closest living relatives are New Guinean canids: the Singing Dogs (NGSD) and recently-described Highland Wild Dog (HWD) of Irian Jaya14–20. More broadly, the dingo, NGSD and HWD (grouped henceforth as NGC = New Guinea Canid) belong to a lineage of East Asian dogs reaching back at least as far as the early Neolithic of eastern China20, but are unique in having avoided admixture from other dog lineages until the introduction of domestic dogs by European colonists21. Members of this lineage previously existing in Mainland and Island Asia appear to have been replaced or assimilated by more recent East Asian and then West Eurasian lineages16, which constitute the present-day dog populations of these regions.

Recent phylogenies emphasise a split between dingoes found in northwestern and southeastern regions of Australia, which may be linked to separate introduction events15,16,20. Interestingly, they also place NGSDs within the dingo clade altogether, as a sister branch to southeastern dingo. Preliminary application of 3D geometric morphometric approaches to dingo cranial morphology have supported this phylogenetic division, finding that the cranial form of southeastern dingoes also bears a closer resemblance to NGCs than their northwestern conspecifics22. However, as some of the similarities between dingoes and NGCs are also shared with domestic dogs, the strength of these conclusions is somewhat limited by the possibility of recent domestic dog hybridisation resulting in equifinality of cranial forms between southeastern dingoes and NGCs.

Accordingly, there is currently a major disjunction between the impressions of dingo ancestry and population history given by genetic/genomic and morphometric approaches respectively, suggesting that the early history of canids in the southern Asia–Pacific region involved shifts in phenotype both in Australia and the broader Asia–Pacific. This issue could be better informed by study of ancient dingoes, but this has thus far been very limited as available specimens are very rare, generally highly incomplete and/or fragmentary, and often represent subadults. The only comparisons of reasonably intact crania in particular were performed several decades ago by Gollan11, using traditional caliper-based measurements, and mainly concern individuals from recent pre-contact times (< 1000BP).

Here, we address this gap by presenting results of dating and morphological analysis of the remains of several ancient dingoes from lunette contexts in southwestern New South Wales (Table 1) at Lake Mungo and Lake Milkengay. These individuals were discovered and excavated from 1979 to 1981 as their positions within older (Pleistocene) strata were revealed by ongoing erosion and deflation of the dried lakes’ lunettes23. Some of these—MD1, MD2 and MLK1—were briefly studied by Gollan11, but the others were discovered after his work. One such individual (MD4) comprises a near-complete skeleton with an intact cranium (Fig. 1), which offered the first opportunity for 3D geometric morphometric comparison of an ancient dingo with present-day dingoes and their international relatives. Table 1 Ancient dingoes from western NSW sampled in this study.

Specimen	Location	Sedimentary context	Description	
MD1	Lake Mungo	Eroding, brown sandy fill in wombat burrow fill intrusive to Zanci unit	Incomplete skeleton including fragments of cranium, mandible and postcranial elements	
MD2	Lake Mungo	Eroding surface of Zanci unit	Fragments of mandible and postcranial elements	
MD3	Lake Mungo	Eroding surface of Zanci unit	Mandible	
MD4	Lake Mungo	Exposed collapsed burrow in Arumpo unit, underlaid by brown sandy sediment	Mostly complete but fragmentary skeleton with adhering concreted sediment	
MD4b	Lake Mungo	Exposed wombat burrow in Arumpo unit	Proximal femur, found nearby to MD4 but from a separate individual	
MD5	Lake Mungo	Adjacent to wombat in eroding channel fill in Zanci unit	Calcaneus	
MD6	Lake Mungo	Exposed wombat burrow in Zanci 100 m east of Dog 3 findspot	Complete humerus and vertebrae	
MLK1	Lake Milkengay	Surface of deflating lunette	Highly fragmented cranium and anterior portion of postcranial skeleton	

Figure 1 Cranium of MD4. (a) left lateral view; (b) right lateral view; (c) dorsal view; (d) ventral view.

Results

AMS dating and dietary stable isotopes

AMS dates of 3035 ± 18BP (3329–3071 calBP) and 2878 ± 18BP (3060–2866 calBP) were obtained from MD1 and MLK1 respectively (Table 2). The calibrated age ranges for MD1 and MLK1, at approximately 3000BP, are both considerably lower than their sedimentary contexts might have initially suggested. It is likely that a similar situation applies to the undated MD4. MD1 is near-contemporaneous with other early, directly-dated specimens from Australia, also in semi-arid environments but over 1500 km away: Madura Cave (3069 ± 27BP8; and Koonalda Cave (3031 ± 34BP)24. Like the Mungo and Milkengay specimens, the Madura and Koonalda fossil specimens appear to have accumulated in the caves naturally rather than having been domestically associated with humans. Together, these dates suggest that wild dingoes were numerous enough by 3000BP to be frequently integrated into the archaeological record. Table 2 AMS dates and dietary isotopes from ancient Lake Mungo and Lake Milkengay dingoes.

Sample ID	Specimen	Date (BP)	δ13C	δ15N	Calibrated age (BP)	
WK-50336	MD1	3035 ± 18 BP	− 16.30	11.23	3330–3070 (95.4%)	
WK-50341	MLK1	2878 ± 18 BP	–	–	3040–2860 (95.4%)	

Dietary δ13C and δ15N stable isotopes were also obtained for MD1. In absence of directly comparable data from this area of Australia, these were compared with 99 published isotopes from dingoes and ancient dogs from elsewhere (Fig. 2). Of these, MD1 is very similar to the younger Madura Cave dingo8, but the older Madura individual is differentiated by very high δ15N comparable to Cook Islands and Pacific Northwest dogs heavily reliant on marine resources25–27. There is no indication of similarity to Neolithic Chinese dogs or the Matja Kuru 2 (Timor) village dog that were provisioned on agricultural produce13,28–30. When compared to recent dingoes from the tropics31, MD1’s δ13C score appears consistent with this sample, but it scores amongst the highest δ15N. This is consistent with most of Morrant et al.’s sample31 reflecting dingoes foraging in mixed/forest environs (C3 plants common) with a smaller proportion from open areas, as opposed to the mixed grassland/heath/Mallee scrub at Lake Mungo. The latter biome is rich in C4 plants, which results in comparably enriched δ15N for local foragers32.Figure 2 Carbon and nitrogen dietary stable isotopes for MD1 compared to modern and ancient dingoes and non-Australian archaeological dogs in the Pacific region.

Previous research has demonstrated that Late Holocene human remains from riverine sites to the southwest of Lake Mungo exhibit higher mean δ15N scores further inland (13.4 ± 1.2) than closer to the coast (10.1 ± 1.1). This is because the former region afforded greater access to terrestrial fauna that subsist on semi-arid C4-rich grasses32. Both ancient riverine populations however exhibit near identical mean δ13C (− 20 ± 0.8 and − 20.1 ± 1.2 respectively), reflecting C3 derived from additional heavy aquatic fish and shellfish consumption32. These observations contextualise MD1’s dietary isotopic results and indicate they are likely to reflect the normal foods of a wild dingo in the semi-arid zone—large marsupials, reptiles and birds33.

3D geometric morphometrics

The relationship between the cranial morphology of MD4, modern dingoes and other Asia–Pacific dogs is visualised in a neighbour-joining tree denoting Procrustes distances (Fig. 3a). Surprisingly, MD4 falls within a cluster of East Asian dogs; the nearest branches observed are two Japanese dogs and the next two are Chinese dogs, with the portion of New Guinean dogs representing Southern Highlands NGSD forming a sister branch to this group. Timorese village dogs cluster with a number of unverified dogs from the Eastern Highlands of New Guinea11, and, interestingly, with a number of dingoes from the southeast coast of NSW. A broad split in the tree between most northwestern and most southeastern dingoes is apparent, with the latter more closely associated with the branches containing New Guinean and Asian dogs, and MD4. Accordingly, modern dingoes near in cranial form to MD4 are exclusively from southeastern Australia, but MD4 is closer to several East Asian and New Guinean dogs than to any sampled dingo.Figure 3 3D geometric morphometric distance-based comparisons of MD4, dingoes and Asia–Pacific dogs. (a) Neighbour-joining tree; (b) group-mean neighbour-joining tree; (c) population mean cranial morphs.

The population-level relationships between dingoes, New Guinean and East Asian dogs, and MD4 are further elucidated in a grouped-mean neighbour-joining tree (Fig. 3b). Modern dingoes cluster and separate in a northwest-southeast fashion, with southeastern dingoes intermediate to New Guinean and Asian dogs. MD4’s branch is closest to East Asian and New Guinean dogs, but maintains a substantial distance from either. This is because these two populations are groupings of convenience, being internally heterogeneous with higher in-group variability (dispersion) than any sampled dingo population (with southeastern dingoes intermediate to these dogs and northwestern dingoes; Table S2).

The traits evident in MD4 responsible for these associations are, compared to modern dingoes: a smaller size; overall shorter and broader shape; a dorsally-rotated (airorhynchous) palate with slightly curved toothrow; a prominent and anterior-rotated frontal-orbital (facial) area; and an elevated cranial height. These traits are evident in the population-mean cranial morphs for dingo and non-dingo groups (Fig. 3c). However, MD4 has an exceptionally well-developed and prominent sagittal crest and large auditory bullae, which are generally rare in East Asian and New Guinean dogs, but strongly characteristic of dingo11. Hence, MD4 displays mixed traits of dingoes and their northern relatives.

The East Asian population includes breeds and landraces from several different regions, whilst variability within New Guinean dogs has been previously noted22 and may reflect different introductions of dog to New Guinea. MD4 is close in conformation to only some of the members of these groups, as illustrated by the lowest Procrustes distances to this individual (Table 3). Of the ten nearest crania, four are Japanese dogs, and five are New Guinean dogs, from the original Southern Highlands breeding group4 and a 19th century Central Province specimen originally classified as “Canis papuensis”34. MD4 is conversely far-removed from the gracile Timorese village dogs and similarly gracile wild dogs from the Mt Wilhelm-Eastern Highlands collection areas in PNG11. Table 3 Ten lowest Procrustes distances to MD4 from sampled canids.

Rank	Specimen ID	Population	Geographic origin	Distance	
1	PRICT-226	Asian dog	Kagoshima, Japan	0.0619	
2	M7989	New Guinea canid	Captive (Taronga Zoo)	0.0621	
3	PRICT-178	Asian dog	Japan	0.0627	
4	PRICT-179	Asian dog	Japan	0.0653	
5	A3652	New Guinea canid	Central Province, PNG	0.0658	
6	PRICT-222	Asian dog	Aichi, Japan	0.0671	
7	M9135	New Guinea canid	Captive (Taronga Zoo)	0.0675	
8	M9185	New Guinea canid	Captive (Taronga Zoo)	0.0677	
9	UNSW-572	Alpine Dingo	Gippsland, Australia	0.0689	
10	BBM-X 157421	New Guinea canid	Captive (Honolulu Zoo)	0.0691	

When the major axes of variation are examined using PCA MD4 also falls outside of the range of modern dingoes from either northwestern or southeastern Australia22, and within the range of Asian and New Guinean dogs (Fig. 4a). MD4’s centroid size tends towards the lowest end of modern dingoes, and is greater than all New Guinean canids (Fig. 4b), but very similar to some Asian dogs (in particular to a dog from Guangzhou, southern China). Linear regression of size on shape determined that the morphological trajectories described by PC1 (20.84%) and PC2 (8.74%) have an allometric component, with size contributing 10.90% of overall form, and 42.38% of PC1 in particular compared to just 3.72% of PC2 (all results significant, p < 0.0001). A secondary (Fig. 4c) PCA using the residuals of this regression to correct for influence of size reveals the association of MD4 with NGCs and Asian dogs, whilst a tertiary PCA (Fig. 4d) with the allometric components of form removed reveals that MD4’s association with Asian dogs specifically is stronger when only shape is considered. These results indicate some of the similarity between MD4 and New Guinean dogs is due to convergence in their overall sizes, which occur within a range outside that observed within most modern dingoes. These observations of similarity between the cranial form of MD4 and Asian dogs were tested through predictive classification using a linear discriminant analysis optimised for unequal sample sizes (sensu lato Evin et al.35; code available Hulme-Beaman36). This returned a predominantly Asian Dog identification followed by NGC, and lastly with a few consecutive combinations of principal components occasionally returning southeastern Australia dingo (see SI for full breakdown of all returned classifications from all consecutive combinations of principal components, and distributions of classification probability based on all resampling iterations). The outstanding position of MD4 relative to modern dingoes cannot be attributed to demographic factors. The specimen is categorically adult as indicated by closed and fused sutures throughout the posterior cranium and its fully erupted (and worn) adult teeth. The very well-developed sagittal crest in particular indicates that it is not a young adult. Although there is some sexual dimorphism in the cranial form of dingoes, this is very minor compared to variation on the basis of geography22. Our modern sample includes numerous individuals of both sex from several geographic populations. There are no signs of pathology or developmental issues in MD4 which might have affected its cranial form. Similarly, there are no indications that the specimen has been deformed by taphonomic processes; although the zygomatic arch was broken and has been digitally reattached, the areas of the cranium which are important in distinguishing MD4 from modern dingoes are located elsewhere in the intact regions of the cranium.Figure 4 Comparisons of MD4 3D cranial morphology with dingoes and Asia–Pacific dogs. (a) PCA of form including size; (b) violin plot of centroid size; (c) PCA of size-adjusted form; (d) PCA of shape, with allometric component removed.

Size of body and carnassial teeth

The estimated shoulder height of ancient Lake Mungo dingoes ranged between approximately 40–47.5 cm (Fig. 5a). MD4 and MD1 are shorter than any of the 117 modern dingoes from around Australia for which height was reconstructed using the same method (46.37–61.50 cm, x̄ = 54.22 ± 29.87 cm), whilst the tallest individual, MD6, is shorter than 97% of these. MD4 and MD6 are taller than all 12 New Guinean canids for which limb bones were available (33.24–41.76 cm, x¯  = 37.5 ± 29.87 cm), though MD1 is very close in height to the largest three: two NGSD from the Southern Highlands and one wild dog of uncertain affiliations from the Central Province mountains4. According to prior research, the intermediate space between dingoes and New Guinean dogs is occupied by various Southeast Asian dogs11–13. Particularly interesting is the closeness of MD1 to the Matja Kuru 2 (Timor) village dog dated to 2967 ± 58BP, at 39.3 cm tall13; and that of MD4 to the Hoekgrot (Java) pariah-type dog from a deposit aged 2655 ± 60 to 3265 ± 55 BP, which was 45 cm tall37.Figure 5 Size measurements from ancient dingoes. (a) Reconstructed statures of ancient Mungo canids compared to overall mean for sampled modern dingo; (b) size of maxillary carnassial; (c) size of mandibular carnassial. Convex hulls have been drawn around limits of modern NW and SE dingo ranges.

Estimated body masses of MD4 and MLK1 were 14.53 kg and 14.27 kg respectively, falling within the range and close to the mean estimates for 95 sampled dingo crania (11.39–22.57 kg, x¯ = 14.25 ± 1.69), also above almost all of the 20 New Guinean dogs sampled (2.49–14.28 kg, x¯ = 8.58 ± 2.90 kg) but comparable to many Southeast Asian dogs11,13,38. Considering MD4’s mass together with its substantially shorter height suggests a stockier build compared to the usual lithe, lean frame of most modern dingoes. As MD4’s palate is unusually broad for its length, a proportion which strongly deviates from modern dingo, there is a possibility it has resulted in an inflated mass estimate.

Ancient Mungo and Milkengay dingoes exhibit maxillary (Fig. 5b) and mandibular (Fig. 5c) carnassial lengths near the minimum of the range of modern dingoes and within the range of many Asian and most New Guinean dogs, but generally maintain a breadth consistent with modern dingo. The exception is MD1 which possesses an exceptionally slim upper carnassial tooth, with equivalents in the study sample only amongst small New Guinean wild dogs of uncertain classification22, and Timorese village dogs. The differences between the dated specimens MD1 and MLK1 suggest substantial variability in the body size of dingoes from the same area within a period of 100–200 years.

Discussion

The dating and morphometric characterisation of ancient dingoes from western NSW informs several key aspects of the dingo’s early history long obscured by prior lack of materials suitable for study. Firstly, our results demonstrate that early dingoes were often smaller than modern dingoes, some intermediate in terms of both body size and dentition to these and New Guinean and Asian dogs. Reduced size for early dingoes was at one time suggested by the smaller teeth of archaeological specimens, but later retracted based on a revised assessment of individual age39. The validity of these prior suspicions are demonstrated through adult specimens of similar or greater antiquity in our study.

3D geometric morphometric characterisation of the MD4 Lake Mungo dingo cranium demonstrates the first clear links between ancient dingoes and dogs to the north of Australia; in particular, with East Asian and New Guinean dogs for which morphological connections were not previously evident. Prior efforts to elucidate the relationships of dingoes to non-Australian dogs by means of morphology have struggled largely due to an assumption that the ancestral dingoes must necessarily closely resemble their present-day descendants. Our results demonstrate that this is not the case, though there is no exact match for MD4 amongst modern-day East Asian or New Guinean dogs. MD4’s cranial morphology might reflect a form common to the ancestors of both dingoes and NGSDs, and by extension, that of an early lineage of East Asian dog which once existed across a large portion of the southern Asia–Pacific region and Australia20.

The divergence of modern dingoes from those described from Lake Mungo and Lake Milkengay suggests that dingoes have increased in size in Australia over the last 3000 years, a change accompanied by a drift in cranial form. An alternative possibility is that the ancient dingoes studied here represent an otherwise extinct lineage which did not successfully adapt, but was absorbed or replaced by an expansion of today’s dominant lineages. Modern dingoes sampled from the Mallee region, which is geographically and ecologically close to the Willandra and Anabranch lake systems, indeed appear as morphological intermediates to the ancient dingoes of this area and modern dingoes from the north and west, in terms of both cranial form and carnassial tooth size.

Finally, our results reveal the deep roots of regional variation in modern dingo morphology along a broadly northwestern-southeastern cline, in accordance with differential maternal and paternal haplotype distributions between these areas15,16. Features such as reduced height, mass and carnassial tooth size observed in modern southeastern dingoes have been explained as resulting from of recent admixture from European breed dogs40,41(Daniels and Corbett 2003; Radford et al. 2012), but have a precedent in dingoes which lived ~ 3000 years before this took place. Other “dog-like” aspects of cranial form in dingoes from southeastern Australia22 also have a precedent in the apparently non-Australian conformation of MD4. Though this of course does not entirely preclude a role of recent European dog admixture, some phenotypic differences observed in contemporary southeastern dingo may be of ancient origin and therefore unrelated.

Ongoing geometric morphometric investigations will further clarify the nature and direction of morphological drift between ancient and modern dingoes, how it relates to the wild and domestic relatives of the dingo, and how it may vary in different biogeographic areas of Australia. Of particular interest to continuing research is the degree to which modern dingoes resemble Asian wolves42, and which specific populations of ancient Asian dogs (a diverse group, like their present-day counterparts) bear the closest resemblance to ancient dingoes. Preliminary expeditionary studies of wild dogs in the Irian Jaya portion of the New Guinea highlands5, and the size of the lone lower carnassial tooth from this area sampled in this study suggests the local presence of an as-yet unstudied and more robust morphotype, which may prove to bear a closer resemblance to ancient dingoes such as MD4 than wild dogs from eastern and southeastern New Guinea.

Materials and methods

Discovery and excavation of ancient canid remains

Canid remains from Lake Mungo were originally identified from 1979 to 1985 in exposed positions eroding from the surface of the lake’s southern lunette, and were retrieved from these locations with minimal excavation. All were found eroding from the Zanci unit (24.5–16 ka), apart from MD4 which is from the older Arumpo unit (30–24.5 ka)43. Wombats became locally extinct before the twentieth century (potentially before European settlement began in the 1830s) but their burrows are common throughout the Mungo lunette. Many feature the bones of other terrestrial species which repurposed the burrows, including the dingoes reported here. The collapse of these burrows preserved their bones within the surrounding sediments, giving the impression of genuine in-situ provenance in Pleistocene-age strata23. At least one of the Mungo dingo skeletons (MD1) exhibits carnivore damage, suggesting it was already deceased for some time when the burrows collapsed.

The Milkengay dingo (MLK1) was discovered during archaeological survey in an eroded portion of the Lake Milkengay lunette containing lithic artefacts and shells, in an area with hearths dated to 3390 ± 500BP and 4590 ± 210BP. As such, it was considered one of, if not the oldest dingo remains in Australia by Gollan11. However, the actual association of the dingo with these archaeological dates was dubious considering that the ages of the dated charcoal samples were stratigraphically inverted, and that the dingo was laterally displaced from them by ~ 10 m11. Due to its fragmented and incomplete nature, study of MLK1 was limited to dental measurements and reconstruction of body mass.

Despite their provenance from lunette areas containing Aboriginal artefacts, more broadly identified as archaeological sites23, there is no indication of direct human association for the ancient dingo remains studied here11. The Mungo dingoes in particular appear to have been exhibiting denning behaviour at the time of their deaths within wombat burrows. The common frequency of skeletal remains of other species within collapsed/filled burrows or hollows in Willandra lunette archaeological sites, particularly bettongs and wombats, further suggests the deposition of these individuals was independent of human activity. However, the possibility that these individuals were tamed in young age before returning to the wild44, or lived in some commensal association with nearby people45, cannot be ruled out.

Community engagement and research design

Archaeological dingoes studied here come from the traditional lands of the Barkandji (Paakantyi), Mutthi Mutthi, and Ngiyampaa peoples. Representatives of these First Nations communities comprising the Willandra Lakes Region World Heritage Aboriginal Advisory Group were consulted during research design, data collection, analysis, results and interpretation, through in-person and digital meetings and progress updates from 2018 to 2021. With community feedback and comments considered at these stages, a final research proposal and submission manuscript was presented to and approved by the WLRWH AAG in 2023.

Dating and isotopes

Elements from each skeleton were subjected to FTIR (Fourier Transform Infrared Spectroscopy) analysis prior to sampling in order to identify those which most likely contain well-preserved collagen and thus yield radiocarbon dates. Based on the results of this analysis, seven samples were taken: from the astragalus of MD1; the mandible of MD3; the mandible and tibia of MD4; the femur of MD4b (a different individual from the same site); the calcaneus of MD5; and the canine tooth of MLK1. Approximately 1–2 g of bone from each specimen was collected for the sample using a Dremel drill.

AMS dating and dietary isotopic work of MD1 and MLK1 was performed at the University of Waikato Radiocarbon Dating Laboratory, Auckland, New Zealand (see Figures S5 and S6. for raw AMS plots). Despite promising FTIR results and good visual condition of sampled elements, the collagen preservation in all Mungo and Milkengay specimens was determined during AMS processing as very poor to non-existent. The only successful samples were from MD1 and MLK1; those submitted from MD3 and MD5 failed, as did two separate samples from the mandible and tibia of MD4. Calibrated ranges for successful samples were produced using the updated SHCal20 calibration curve for the southern hemisphere46.

Morphological analyses

MD4’s cranium was 3D scanned and landmarked at 45 points along the midline and left side of the cranium, following the methodology of Koungoulos22, which uses a Creaform 700 Handyscan to digitise crania and Stratovan Checkpoint47 to apply landmarks (Fig. S2). As the left zygomatic arch of MD4 was broken prior to excavation, this anatomical feature was digitally reassembled in Blender before landmarking following the protocols outlined by Lautenschlager48. MD4 was compared to a dataset of 475 dingo and dog crania scanned following Koungoulos22, excepting a portion of the East Asian sample which were obtained as open-source CT scans via KUPRI (Table S1). Modern dingo crania were assigned to 14 subpopulations based on the biogeographic area in which they were originally collected (Fig. S7). The dog populations were divided into dogs from New Guinea (including Irian Jaya and New Britain), and dogs from Asia.

Specimens were aligned using a Generalised Procrustes Analysis (GPA using Morpho49, and then between known modern group differences (i.e. removing the archaeological specimen) were tested using a permuted ANOVA using geomorph50. Dispersions (i.e. group morphological diversity) was tested using a permuted ANOVA on distances from the group centroids using geomorph50. The archaeological specimen was separated from the modern specimens; the raw data of the known specimens were resampled to equal sample size, and for each resampling iteration a new GPA was carried out and followed by a principal components analysis (PCA) and then a linear discriminant analysis on a reduced number of PCs. The archaeological specimen was projected into this GPA space based around equal sample size by aligning the archaeological specimen to the new GPA mean shape using Morpho functions49and then the specimen was projected using the predict function into PCA space, then LDA space for each iteration. Classification results vary depending on the number of PCs used, therefore the process was repeated for all consecutive combinations of PCs to the smallest known group sample size minus 1 (i.e. PCs 1–15, by testing each combination of PCs 1–2, then 1–3, 1–4 and so on). The posterior probability for the archaeological specimen was recorded in every iteration and PC combination and plotted as probability distributions for each PC combination. See SI R markdown for full work through of the data with all corresponding tests, and plots generated.

All 3D geometric morphometric analyses were performed in Rstudio (v.2022.07.2+576) using the packages ape51, geomorph50, GeoOrigins52, MASS53, Morpho49, and OptLDA36. A full record of all packages and versions used is located in the Supplementary Information. Unrooted neighbour-joining trees (BIONJ) based both on individual scores and population means were constructed to more clearly illustrate MD4’s relations in terms of cranial morphology54,55. Measurements of centroid size, and PCAs based on raw data and residuals from a regression of shape on size were performed to further explore these relationships and determine the role of size in characterisations of MD4’s nearest morphological relations.

Secondary metric analyses involved comparisons of three elements of body size for MD4 and other archaeological dingoes. Height at the shoulder was reconstructed from length of limb bones using Harcourt’s method56 for MD1, MD4 and MD6. The body mass of MD4 and MLK1 was estimated using measurements of the palate length and breadth, taken digitally from the 3D scan models, following the Lucas et al.57 dingo-specific method. These were compared with estimates for modern dingoes and New Guinean dogs obtained in the same manner from the crania sampled in this research. Measurements of carnassial teeth length and breadth were taken digitally from the 3D scans.

Supplementary Information

Supplementary Information.

Supplementary Information.

Supplementary Information

The online version contains supplementary material available at 10.1038/s41598-024-65729-3.

Acknowledgements

Acknowledgements are owed to the original excavators and collectors of the ancient dingo remains that form the focus of this paper: Peter Clark, Peter Brown, Klim Gollan, Jeanette Hope and Mike McIntyre, and to Allison Dejanovic of the Australian Museum who facilitated extended study access to the specimens. We thank the Barkandji (Paakantji), Mutthi Mutthi and Ngiyampaa communities of western NSW for their endorsement of this research, and are grateful to Dan Rosendahl, Tressna Martin, Aaron Fogel and Pat Faulkner for their invaluable assistance in facilitating community engagements. For facilitating data collection from the comparative specimens of modern dingoes and dogs, we thank Sandy Ingleby and Harry Parnaby (Australian Museum), Karen Roberts (Museum Victoria), Denise Donlon (Shellshear Museum), Leo Joseph and Chris Wilson (CSIRO National Wildlife Collection), Sofia Samper-Carro (Australian National University), Heather Janetzki (Queensland Museum), Gavin Dally (Museum and Art Gallery Northern Territory), Jude Philp (Macleay Museum), David Stemmer (South Australian Museum), Mike Letnic (University of New South Wales), Kenny Travouillon (Western Australian Museum) and Molly Hagemann (Bernice Bishop Museum). We also acknowledge the Kyoto University Primatology Research Institute’s Digital Morphology Museum (KUPRI DMM) for providing 3D-scan data, and the private collectors and institutions (KUPRI) who provided the original specimens; we are additionally grateful for Naoko Egi’s assistance in determining their affiliations. We extend thanks to Julie Sebanc-Butler for performing the preliminary FTIR analysis of the Mungo and Milkengay skeletal remains; to Fiona Petchey and Lorena Becerra-Valdivia for assistance with dating and isotopic analysis; Nicola Stern, Nathan Jankowski and Jim Bowler for expert insights into the sedimentology and stratigraphy of Lake Mungo’s southern lunette; and to Harvey Johnston and Peter Brown for first-hand information regarding the discovery and provenance of many of the Mungo dingoes. We are also grateful to Damian Morrant for access and permission to use his isotopic data from modern dingoes. This research was funded by an Australian Government Research Training Program scholarship, and the Carlyle-Greenwell Research Scholarship from the University of Sydney. Assistance with publication expenses was also provided by the Ben Sandford Cullen Scholarship from the University of Sydney and the Deputy Vice Chancellor, Research at the University of New England.

Author contributions

L.G.K. conceived the study under guidelines of and consultation with the WLRWHAAG. L.G.K. and M.F. designed and conducted sampling procedures for AMS and isotopic dating. L.G.K. collected morphometric data used in the study. L.G.K. and A.H-B. conducted the geometric morphometric statistical analyses and prepared corresponding figures together. L.G.K. wrote main text and supplementary information file. All authors contributed to text of Discussion, reviewed and approved the final manuscript.

Data availability

Data is provided within the manuscript or supplementary information files.

Competing interests

The authors declare no competing interests.

Publisher's note

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

A list of authors and their affiliations appears at the end of the paper.
==== Refs
References

1. Jackson SM The Dogma of Dingoes—Taxonomic status of the dingo: A reply to Smith et al. Zootaxa 2019 4564 1 198 212 10.11646/zootaxa.4564.1.7
Jackson, S. M. et al. The Dogma of Dingoes—Taxonomic status of the dingo: A reply to Smith et al.. Zootaxa 4564(1), 198–212. 10.11646/zootaxa.4564.1.7 (2019).
2. Smith B Taxonomic status of the Australian dingo: The case for Canis dingo Meyer, 1793 Zootaxa 2019 4564 1 173 197 10.11646/zootaxa.4564.1.6
Smith, B. et al. Taxonomic status of the Australian dingo: The case for Canis dingo Meyer, 1793. Zootaxa 4564(1), 173–197. 10.11646/zootaxa.4564.1.6 (2019).
3. Boitani L Ciucci P Comparative social ecology of feral dogs and wolves Ethol. Ecol. Evol. 1995 7 1 49 72 10.1080/08927014.1995.9522969
Boitani, L. & Ciucci, P. Comparative social ecology of feral dogs and wolves. Ethol. Ecol. Evol. 7(1), 49–72. 10.1080/08927014.1995.9522969 (1995).
4. Koler-Matznick J Brisbin IL Feinstein M Bulmer S An updated description of the New Guinea singing dog (Canis hallstromi, Troughton 1957) J. Zool. 2003 261 2 109 118 10.1017/S0952836903004060
Koler-Matznick, J., Brisbin, I. L., Feinstein, M. & Bulmer, S. An updated description of the New Guinea singing dog (Canis hallstromi, Troughton 1957). J. Zool. 261(2), 109–118. 10.1017/S0952836903004060 (2003).
5. McIntyre JK Wolf LL Sacks BN Koibur J Brisbin IL Jr A population of free-living highland wild dogs in Indonesian Papua Austral. Mammal. 2019 42 2 160 166 10.1071/AM18039
McIntyre, J. K., Wolf, L. L., Sacks, B. N., Koibur, J. & Brisbin, I. L. Jr. A population of free-living highland wild dogs in Indonesian Papua. Austral. Mammal. 42(2), 160–166. 10.1071/AM18039 (2019).
6. Gering E Incorvaia D Henriksen R Conner J Getty T Wright D Getting back to nature: Feralization in animals and plants Trends Ecol. Evol. 2019 34 12 1 15 10.1016/j.tree.2019.07.018
Gering, E. et al. Getting back to nature: Feralization in animals and plants. Trends Ecol. Evol. 34(12), 1–15. 10.1016/j.tree.2019.07.018 (2019).
7. Koungoulos L Fillios M Hunting dogs down under? On the Aboriginal use of tame dingoes in dietary game acquisition and its relevance to Australian prehistory J. Anthropol. Archaeol. 2020 58 101146 10.1016/j.jaa.2020.101146
Koungoulos, L. & Fillios, M. Hunting dogs down under? On the Aboriginal use of tame dingoes in dietary game acquisition and its relevance to Australian prehistory. J. Anthropol. Archaeol. 58, 101146. 10.1016/j.jaa.2020.101146 (2020).
8. Balme J O’Connor S Fallon S New dates on dingo bones from Madura Cave provide oldest firm evidence for arrival of the species in Australia Sci. Rep. 2018 8 9933 10.1038/s41598-018-28324-x 30026564
Balme, J., O’Connor, S. & Fallon, S. New dates on dingo bones from Madura Cave provide oldest firm evidence for arrival of the species in Australia. Sci. Rep. 8, 9933. 10.1038/s41598-018-28324-x (2018).30026564
9. Greig K Complex history of dog (Canis familiaris) origins and translocations in the Pacific revealed by ancient mitogenomes Sci. Rep. 2018 8 9130 10.1038/s41598-018-27363-8 29904060
Greig, K. et al. Complex history of dog (Canis familiaris) origins and translocations in the Pacific revealed by ancient mitogenomes. Sci. Rep. 8, 9130. 10.1038/s41598-018-27363-8 (2018).29904060
10. Fillios MA Taçon PSC Who let the dogs in? A review of the recent genetic evidence for the introduction of the dingo to Australia and implications for the movement of people J. Archaeol. Sci. Rep. 2016 7 782 792 10.1016/j.jasrep.2016.03.001
Fillios, M. A. & Taçon, P. S. C. Who let the dogs in? A review of the recent genetic evidence for the introduction of the dingo to Australia and implications for the movement of people. J. Archaeol. Sci. Rep. 7, 782–792. 10.1016/j.jasrep.2016.03.001 (2016).
11. Gollan, K. Prehistoric Dingo (Doctor of Philosophy Thesis, Australian National University, 1982)10.25911/5d6cfce36fe39.
12. Corbett LK Morphological comparisons of Australian and Thai dingoes: A reappraisal of dingo status, distribution and ancestry Proc. Ecol. Soc. Austral. 1985 13 277 291
Corbett, L. K. Morphological comparisons of Australian and Thai dingoes: A reappraisal of dingo status, distribution and ancestry. Proc. Ecol. Soc. Austral. 13, 277–291 (1985).
13. Gonzalez A Clark G O’Connor S Matisoo-Smith L A 3000 year old dog burial in Timor-Leste Austral. Archaeol. 2013 76 1 13 20 10.1080/03122417.2013.11681961
Gonzalez, A., Clark, G., O’Connor, S. & Matisoo-Smith, L. A 3000 year old dog burial in Timor-Leste. Austral. Archaeol. 76(1), 13–20. 10.1080/03122417.2013.11681961 (2013).
14. Ardalan A Narrow genetic basis for the Australian dingo confirmed through analysis of paternal ancestry Genetica 2012 140 1–3 65 73 10.1007/s10709-012-9658-5 22618967
Ardalan, A. et al. Narrow genetic basis for the Australian dingo confirmed through analysis of paternal ancestry. Genetica 140(1–3), 65–73. 10.1007/s10709-012-9658-5 (2012).22618967
15. Cairns JM Brown SK Sacks BN Ballard JWO Conservation implications for dingoes from the maternal and paternal genome: Multiple populations, dog introgression, and demography Ecol. Evol. 2017 7 22 9787 9807 10.1002/ece3.3487 29188009
Cairns, J. M., Brown, S. K., Sacks, B. N. & Ballard, J. W. O. Conservation implications for dingoes from the maternal and paternal genome: Multiple populations, dog introgression, and demography. Ecol. Evol. 7(22), 9787–9807. 10.1002/ece3.3487 (2017).29188009
16. Cairns KM Shannon LM Koler-Matznick J Ballard JWO Boyko AR Elucidating biogeographical patterns in Australian native canids using genome wide SNPs PLoS ONE 2018 13 e0198754 10.1371/journal.pone.0198754 29889854
Cairns, K. M., Shannon, L. M., Koler-Matznick, J., Ballard, J. W. O. & Boyko, A. R. Elucidating biogeographical patterns in Australian native canids using genome wide SNPs. PLoS ONE 13, e0198754. 10.1371/journal.pone.0198754 (2018).29889854
17. Oskarsson MCR Mitochondrial DNA data indicate an introduction through Mainland Southeast Asia for Australian dingoes and Polynesian domestic dogs Proc. R. Soc. B Biol. Sci. 2011 279 1730 967 974 10.1098/rspb.2011.1395
Oskarsson, M. C. R. et al. Mitochondrial DNA data indicate an introduction through Mainland Southeast Asia for Australian dingoes and Polynesian domestic dogs. Proc. R. Soc. B Biol. Sci. 279(1730), 967–974. 10.1098/rspb.2011.1395 (2011).
18. Savolainen P Leitner T Wilton AN Matisoo-Smith E Lundeberg J A detailed picture of the origin of the Australian dingo, obtained from the study of mitochondrial DNA Proc. Natl. Acad. Sci. USA 2004 101 33 12387 12390 10.1073/pnas.0401814101 15299143
Savolainen, P., Leitner, T., Wilton, A. N., Matisoo-Smith, E. & Lundeberg, J. A detailed picture of the origin of the Australian dingo, obtained from the study of mitochondrial DNA. Proc. Natl. Acad. Sci. USA 101(33), 12387–12390. 10.1073/pnas.0401814101 (2004).15299143
19. Surbakti S New Guinea highland wild dogs are the original New Guinea singing dogs Proc. Natl. Acad. Sci. 2020 10.1073/pnas.2007242117 32868416
Surbakti, S. et al. New Guinea highland wild dogs are the original New Guinea singing dogs. Proc. Natl. Acad. Sci.10.1073/pnas.2007242117 (2020).32868416
20. Zhang M Ancient DNA evidence from China reveals the expansion of Pacific dogs Mol. Biol. Evol. 2020 37 5 1462 1469 10.1093/molbev/msz311 31913480
Zhang, M. et al. Ancient DNA evidence from China reveals the expansion of Pacific dogs. Mol. Biol. Evol. 37(5), 1462–1469. 10.1093/molbev/msz311 (2020).31913480
21. Bergström A Origins and genetic legacy of prehistoric dogs Science 2020 370 6516 557 564 10.1126/science.aba9572 33122379
Bergström, A. et al. Origins and genetic legacy of prehistoric dogs. Science 370(6516), 557–564. 10.1126/science.aba9572 (2020).33122379
22. Koungoulos L Old dogs, new tricks: 3D geometric analysis of cranial morphology supports ancient population substructure in the Australian dingo Zoomorphology 2020 139 2 263 275 10.1007/s00435-019-00475-z
Koungoulos, L. Old dogs, new tricks: 3D geometric analysis of cranial morphology supports ancient population substructure in the Australian dingo. Zoomorphology 139(2), 263–275. 10.1007/s00435-019-00475-z (2020).
23. Clark, P.M. Willandra Lakes World Heritage Area Archaeological Resource Study (NSW Department of Environment and Planning, Buronga, 1987).
24. Walshe K Koonalda Cave, Nullarbor Plain, South Australia—issues in optical and radiometric dating of deep karst caves Geochronometria 2017 44 1 366 373 10.1515/geochr-2015-0081
Walshe, K. Koonalda Cave, Nullarbor Plain, South Australia—issues in optical and radiometric dating of deep karst caves. Geochronometria 44(1), 366–373. 10.1515/geochr-2015-0081 (2017).
25. Ames KM Stable isotope and ancient DNA analysis of dog remains from Cathlapotle (45CL1), a contact-era site on the Lower Columbia River J. Archaeol. Sci. 2015 57 268 282 10.1016/j.jas.2015.02.038
Ames, K. M. et al. Stable isotope and ancient DNA analysis of dog remains from Cathlapotle (45CL1), a contact-era site on the Lower Columbia River. J. Archaeol. Sci. 57, 268–282. 10.1016/j.jas.2015.02.038 (2015).
26. Craig, J.A. Stable isotope analysis of prehistoric human and commensal diet on Aitutaki, southern Cook Islands (Doctor of Philosophy Thesis, University of Auckland, 2009).
27. Hillis D McKechnie I Guiry E St. Claire DE Darimont CT Ancient dog diets on the Pacific Northwest Coast: Zooarchaeological and stable isotope modelling evidence from Tseshaht territory and beyond Sci. Rep. 2020 10 15630 10.1038/s41598-020-71574-x 33004834
Hillis, D., McKechnie, I., Guiry, E., St. Claire, D. E. & Darimont, C. T. Ancient dog diets on the Pacific Northwest Coast: Zooarchaeological and stable isotope modelling evidence from Tseshaht territory and beyond. Sci. Rep. 10, 15630. 10.1038/s41598-020-71574-x (2020).33004834
28. Barton L Agricultural origins and the isotopic identity of domestication in northern China Proc. Natl. Acad. Sci. 2009 106 14 5523 5528 10.1073/pnas.080996010 19307567
Barton, L. et al. Agricultural origins and the isotopic identity of domestication in northern China. Proc. Natl. Acad. Sci. 106(14), 5523–5528. 10.1073/pnas.080996010 (2009).19307567
29. Chen XL Raising practices of neolithic livestock evidenced by stable isotope analysis in the Wei River Valley, North China Int. J. Osteoarchaeol. 2016 26 1 42 52 10.1002/oa.2393
Chen, X. L. et al. Raising practices of neolithic livestock evidenced by stable isotope analysis in the Wei River Valley, North China. Int. J. Osteoarchaeol. 26(1), 42–52. 10.1002/oa.2393 (2016).
30. Pechenkina EA Ambrose SH Xiaolin M Benfer RA Reconstructing northern Chinese Neolithic subsistence practices by isotopic analysis J. Archaeol. Sci. 2005 32 8 1176 1189 10.1016/j.jas.2005.02.015
Pechenkina, E. A., Ambrose, S. H., Xiaolin, M. & Benfer, R. A. Reconstructing northern Chinese Neolithic subsistence practices by isotopic analysis. J. Archaeol. Sci. 32(8), 1176–1189. 10.1016/j.jas.2005.02.015 (2005).
31. Morrant DS Wurster CM Johnson CN Butler JRA Congdon BC Prey use by dingoes in a contested landscape: Ecosystem service provider or biodiversity threat? Ecol. Evol. 2017 7 21 8927 8935 10.1002/ece3.3345 29152188
Morrant, D. S., Wurster, C. M., Johnson, C. N., Butler, J. R. A. & Congdon, B. C. Prey use by dingoes in a contested landscape: Ecosystem service provider or biodiversity threat?. Ecol. Evol. 7(21), 8927–8935. 10.1002/ece3.3345 (2017).29152188
32. Pate FD Owen TD Stable carbon and nitrogen isotopes as indicators of sedentism and territoriality in late Holocene South Australia Archaeol. Oceania 2014 49 2 86 94 10.1002/arco.5019
Pate, F. D. & Owen, T. D. Stable carbon and nitrogen isotopes as indicators of sedentism and territoriality in late Holocene South Australia. Archaeol. Oceania 49(2), 86–94. 10.1002/arco.5019 (2014).
33. Doherty TS Continental patterns in the diet of a top predator: Australia's dingo Mamm. Rev. 2018 49 1 31 44 10.1111/mam.12139
Doherty, T. S. et al. Continental patterns in the diet of a top predator: Australia’s dingo. Mamm. Rev. 49(1), 31–44. 10.1111/mam.12139 (2018).
34. Parnaby HE Ingleby S Divljan A Type specimens of non-fossil mammals in the Australian Museum, Sydney Rec. Austral. Museum 2017 69 5 277 420 10.3853/j.2201-4349.69.2017.1653
Parnaby, H. E., Ingleby, S. & Divljan, A. Type specimens of non-fossil mammals in the Australian Museum, Sydney. Rec. Austral. Museum 69(5), 277–420. 10.3853/j.2201-4349.69.2017.1653 (2017).
35. Evin A Unravelling the complexity of domestication: A case study using morphometrics and ancient DNA analyses of archaeological pigs from Romania Philos. Trans. R. Soc. B Biol. Sci. 2015 370 166 20130616 10.1098/rstb.2013.0616
Evin, A. et al. Unravelling the complexity of domestication: A case study using morphometrics and ancient DNA analyses of archaeological pigs from Romania. Philos. Trans. R. Soc. B Biol. Sci. 370(166), 20130616. 10.1098/rstb.2013.0616 (2015).
36. Hulme-Beaman, A. OPTLDA: LDA optimisation with resampling to equal sample size. https://github.com/ArdernHB/OptLDA (2022).
37. Storm P Resten van een hond uit de vindplaats Hoekgrot (Java) Cranium 2001 18 1 31 44
Storm, P. Resten van een hond uit de vindplaats Hoekgrot (Java). Cranium 18(1), 31–44 (2001).
38. Corbett LK The Dingo in Australia and Asia 1995 Cornell University Press
Corbett, L. K. The Dingo in Australia and Asia (Cornell University Press, 1995).
39. Macintosh NWG Fox MW 'The origin of the dingo: An enigma The Wild Canids: Their Systematics, Behavioural Ecology and Evolution 1975 Van Nostrand Reinhold Company 87 106
Macintosh, N. W. G. ’The origin of the dingo: An enigma. In The Wild Canids: Their Systematics, Behavioural Ecology and Evolution (ed. Fox, M. W.) 87–106 (Van Nostrand Reinhold Company, 1975).
40. Daniels MJ Corbett LK Redefining introgressed protected mammals: When is a wildcat a wild cat and a dingo a wild dog? Wild. Res. 2003 30 3 213 218 10.1071/WR02045
Daniels, M. J. & Corbett, L. K. Redefining introgressed protected mammals: When is a wildcat a wild cat and a dingo a wild dog?. Wild. Res. 30(3), 213–218. 10.1071/WR02045 (2003).
41. Radford CG Letnic M Fillios M Crowther MS An assessment of the taxonomic status of wild canids in south-eastern New South Wales: Phenotypic variation in dingoes Austral. J. Zool. 2012 60 2 73 80 10.1071/ZO12006
Radford, C. G., Letnic, M., Fillios, M. & Crowther, M. S. An assessment of the taxonomic status of wild canids in south-eastern New South Wales: Phenotypic variation in dingoes. Austral. J. Zool. 60(2), 73–80. 10.1071/ZO12006 (2012).
42. Gonzalez, A. The Pariah Case: Some Comments on the Origin and Evolution of Primitive Dogs and on the Taxonomy of Related Species (Doctor of Philosophy Thesis, The Australian National University, 2012)10.25911/5d5e766e646dc.
43. Jankowski NR Stern N Lachlan TJ Jacobs Z A high-resolution late Quaternary depositional history and chronology for the southern portion of the Lake Mungo lunette, semi-arid Australia Quat. Sci. Rev. 2020 233 106224 10.1016/j.quascirev.2020.106224
Jankowski, N. R., Stern, N., Lachlan, T. J. & Jacobs, Z. A high-resolution late Quaternary depositional history and chronology for the southern portion of the Lake Mungo lunette, semi-arid Australia. Quat. Sci. Rev. 233, 106224. 10.1016/j.quascirev.2020.106224 (2020).
44. Smith BP Litchfield CA A review of the relationship between Indigenous Australians, dingoes (Canis dingo) and domestic dogs (Canis familiaris) Anthrozoös 2013 22 2 111 128 10.2752/175303709X434149
Smith, B. P. & Litchfield, C. A. A review of the relationship between Indigenous Australians, dingoes (Canis dingo) and domestic dogs (Canis familiaris). Anthrozoös 22(2), 111–128. 10.2752/175303709X434149 (2013).
45. Brumm A Dingoes and domestication Archaeol. Oceania 2021 56 1 17 31 10.1002/arco.5226
Brumm, A. Dingoes and domestication. Archaeol. Oceania 56(1), 17–31. 10.1002/arco.5226 (2021).
46. Hogg AG SHCal20 Southern Hemisphere Calibration, 0–55,000 Years cal BP Radiocarbon 2020 62 4 759 778 10.1017/RDC.2020.59
Hogg, A. G. et al. SHCal20 Southern Hemisphere Calibration, 0–55,000 Years cal BP. Radiocarbon 62(4), 759–778. 10.1017/RDC.2020.59 (2020).
47. Stratovan Corporation. Stratovan Checkpoint. 2018.08.07 edn.(2018).
48. Lautenschlager S Reconstructing the past: Methods and techniques for the digital restoration of fossils R. Soc. Open Sci. 2016 3 10 160342 10.1098/rsos.160342 27853548
Lautenschlager, S. Reconstructing the past: Methods and techniques for the digital restoration of fossils. R. Soc. Open Sci. 3(10), 160342. 10.1098/rsos.160342 (2016).27853548
49. Schlager S Zheng G Li S Székely G Chapter 9 – Morpho and Rvcg – shape analysis in R: R-packages for geometric morphometrics, shape analysis and surface manipulations Statistical Shape and Deformation Analysis: Methods, Implementation and Applications 2017 Academic Press 217 256
Schlager, S. Chapter 9 – Morpho and Rvcg – shape analysis in R: R-packages for geometric morphometrics, shape analysis and surface manipulations. In Statistical Shape and Deformation Analysis: Methods, Implementation and Applications (eds Zheng, G. et al.) 217–256 (Academic Press, 2017). 10.1016/B978-0-12-810493-4.00011-0.
50. Adams DC Otárola-Castillo geomorph: An r package for the collection and analysis of geometric morphometric shape data Methods Ecol. Evol. 2013 4 4 393 399 10.1111/2041-210X.12035
Adams, D. C. & Otárola-Castillo,. geomorph: An r package for the collection and analysis of geometric morphometric shape data. Methods Ecol. Evol. 4(4), 393–399. 10.1111/2041-210X.12035 (2013).
51. Paradis E Schliep K ape 5.0: An environment for modern phylogenetics and evolutionary analyses in R Bioinformatics 2019 35 3 536 528 10.1093/bioinformatics/bty633
Paradis, E. & Schliep, K. ape 5.0: An environment for modern phylogenetics and evolutionary analyses in R. Bioinformatics 35(3), 536–528. 10.1093/bioinformatics/bty633 (2019).
52. Hulme-Beaman A geoorigins: A new method and r package for trait mapping and geographic provenancing of specimens without categorical constraints Methods Ecol. Evol. 2020 11 10 1247 1257 10.1111/2041-210X.13444
Hulme-Beaman, A. et al. geoorigins: A new method and r package for trait mapping and geographic provenancing of specimens without categorical constraints. Methods Ecol. Evol. 11(10), 1247–1257. 10.1111/2041-210X.13444 (2020).
53. Venables WN Ripley BD Modern Applied Statistics with S 2002 4 Springer
Venables, W. N. & Ripley, B. D. Modern Applied Statistics with S 4th edn. (Springer, 2002).
54. Gascuel O BIONJ: An improved version of the NJ algorithm based on a simple model of sequence data Mol. Biol. Evol. 1997 14 7 685 695 10.1093/oxfordjournals.molbev.a025808 9254330
Gascuel, O. BIONJ: An improved version of the NJ algorithm based on a simple model of sequence data. Mol. Biol. Evol. 14(7), 685–695. 10.1093/oxfordjournals.molbev.a025808 (1997).9254330
55. Saitou N Nei M The neighbor-joining method: A new method for reconstructing phylogenetic trees Mol. Biol. Evol. 1987 4 4 406 425 10.1093/oxfordjournals.molbev.a040454 3447015
Saitou, N. & Nei, M. The neighbor-joining method: A new method for reconstructing phylogenetic trees. Mol. Biol. Evol. 4(4), 406–425. 10.1093/oxfordjournals.molbev.a040454 (1987).3447015
56. Harcourt RA The dog in prehistoric and early historic Britain J. Archaeol. Sci. 1974 1 2 151 175 10.1016/0305-4403(74)90040-5
Harcourt, R. A. The dog in prehistoric and early historic Britain. J. Archaeol. Sci. 1(2), 151–175. 10.1016/0305-4403(74)90040-5 (1974).
57. Lucas T Smith BP Norris RM Henneberg M Reconstructing body mass of the Australian dingo (Canis dingo) from two simple measurements of the hard palate J. Archaeol. Sci. Rep. 2019 23 534 539 10.1016/j.jasrep.2018.11.018
Lucas, T., Smith, B. P., Norris, R. M. & Henneberg, M. Reconstructing body mass of the Australian dingo (Canis dingo) from two simple measurements of the hard palate. J. Archaeol. Sci. Rep. 23, 534–539. 10.1016/j.jasrep.2018.11.018 (2019).
