Elucidating biogeographical patterns in Australian native canids using genome wide SNPs

A global assessment of large terrestrial carnivore kill rates

Through killing and instilling fear in their prey, large terrestrial carnivores shape the structure and function of ecosystems globally. Most large carnivore species have experienced severe range and population declines due to human activities, and many are now threatened with extinction. Consequently, the impacts of these predators on food webs have been diminished or lost completely from many ecosystems. Kill rates provide a fundamental metric for understanding large carnivore ecology and assessing and comparing predation within and across ecological communities. Our systematic review of large terrestrial mammalian carnivore kill rates reveals significant positive geographic (North America, Europe, and Africa) and taxonomic (grey wolf Canis lupus, puma Puma concolor, lion Panthera leo, and Eurasian lynx Lynx lynx) bias, with most studies apparently motivated by human-carnivore conflict over access to ungulate prey and wildlife management objectives. Our current understanding of the behaviour and functional roles of many large carnivore species and populations thus remains limited. By synthesising and comparing kill rates, we show that solitary carnivores (e.g. brown bears Ursus arctos and most felids) exhibit higher per capita kill rates than social carnivores. However, ungulate predation by bears is typically limited to predation of neonates during a short period. Lower per capita kill rates by social carnivores suggests group living significantly reduces energetic demands, or, alternatively, that group-living carnivores defend and consume a greater proportion of large prey carcasses, or may acquire more food through other means (e.g. scavenging, kleptoparasitism) than solitary hunters. Kill and consumption rates for Canidae - measured as kilograms of prey per kilogram of carnivore per day - are positively correlated with body mass, consistent with increasing energy costs associated with a cursorial hunting strategy. By contrast, ambush predators such as felids show an opposite trend, and thus the potential energetic advantage of an ambush hunting strategy for carnivores as body mass increases. Additionally, ungulate kill rates remain relatively constant across solitary felid body sizes, indicative of energetic constraints and optimal foraging. Kill rate estimates also reveal potential insights into trophic structuring within carnivore guilds, with subordinate carnivores often killing more than their larger counterparts, which may be indicative of having to cope with food losses to scavengers and dominant competitors. Subordinate carnivores may thus serve an important role in provisioning food to other trophic levels within their respective ecosystems. Importantly, kill rates also clarify misconceptions around the predatory behaviour of carnivores (e.g. spotted hyaenas Crocuta crocuta and wolverines Gulo gulo are often considered scavengers rather than the capable hunters that they are) and thus the potential impacts of various carnivore species on their ecological communities. Despite the importance of kill rates in understanding predator-prey interactions, their utility is not widely recognised, and insufficient research limits our ability to fully appreciate and predict the consequences of modified predation regimes, justify current management actions affecting carnivores, or inform effective conservation measures. Together with other important research on predator-prey interactions, robust kill rate studies that address the research deficiencies we highlight will provide a deeper understanding of the foraging behaviours and potential ecosystem impacts of many of the world's carnivores, thus aiding effective conservation and management actions.

Phenotypic diversity in early Australian dingoes revealed by traditional and 3D geometric morphometric analysis

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.

Ancient genomes reveal over two thousand years of dingo population structure

Dingoes are culturally and ecologically important free-living canids whose ancestors arrived in Australia over 3,000 B.P., likely transported by seafaring people. However, the early history of dingoes in Australia-including the number of founding populations and their routes of introduction-remains uncertain. This uncertainty arises partly from the complex and poorly understood relationship between modern dingoes and New Guinea singing dogs, and suspicions that post-Colonial hybridization has introduced recent domestic dog ancestry into the genomes of many wild dingo populations. In this study, we analyzed genome-wide data from nine ancient dingo specimens ranging in age from 400 to 2,746 y old, predating the introduction of domestic dogs to Australia by European colonists. We uncovered evidence that the continent-wide population structure observed in modern dingo populations had already emerged several thousand years ago. We also detected excess allele sharing between New Guinea singing dogs and ancient dingoes from coastal New South Wales (NSW) compared to ancient dingoes from southern Australia, irrespective of any post-Colonial hybrid ancestry in the genomes of modern individuals. Our results are consistent with several demographic scenarios, including a scenario where the ancestry of dingoes from the east coast of Australia results from at least two waves of migration from source populations with varying affinities to New Guinea singing dogs. We also contribute to the growing body of evidence that modern dingoes derive little genomic ancestry from post-Colonial hybridization with other domestic dog lineages, instead descending primarily from ancient canids introduced to Sahul thousands of years ago.

Genomic Consequences of Isolation and Inbreeding in an Island Dingo Population

Dingoes come from an ancient canid lineage that originated in East Asia around 8,000 to 11,000 years BP. As Australia's largest terrestrial predator, dingoes play an important ecological role. A small, protected population exists on a world heritage listed offshore island, K'gari (formerly Fraser Island). Concern regarding the persistence of dingoes on K'gari has risen due to their low genetic diversity and elevated inbreeding levels. However, whole-genome sequence data is lacking from this population. Here, we include five new whole-genome sequences of K'gari dingoes. We analyze a total of 18 whole-genome sequences of dingoes sampled from mainland Australia and K'gari to assess the genomic consequences of their demographic histories. Long (>1 Mb) runs of homozygosity (ROHs)-indicators of inbreeding-are elevated in all sampled dingoes. However, K'gari dingoes showed significantly higher levels of very long ROH (>5 Mb), providing genomic evidence for small population size, isolation, inbreeding, and a strong founder effect. Our results suggest that, despite current levels of inbreeding, the K'gari population is purging strongly deleterious mutations, which, in the absence of further reductions in population size, may facilitate the persistence of small populations despite low genetic diversity and isolation. However, there may be little to no purging of mildly deleterious alleles, which may have important long-term consequences, and should be considered by conservation and management programs.

Genomic consequences of isolation and inbreeding in an island dingo population

Dingoes come from an ancient canid lineage that originated in East Asia around 8000-11,000 years BP. As Australia's largest terrestrial predator, dingoes play an important ecological role. A small, protected population exists on a world heritage listed offshore island, K'gari (formerly Fraser Island). Concern regarding the persistence of dingoes on K'gari has risen due to their low genetic diversity and elevated inbreeding levels. However, whole-genome sequencing data is lacking from this population. Here, we include five new whole-genome sequences of K'gari dingoes. We analyze a total of 18 whole genome sequences of dingoes sampled from mainland Australia and K'gari to assess the genomic consequences of their demographic histories. Long (>1 Mb) runs of homozygosity (ROH) - indicators of inbreeding - are elevated in all sampled dingoes. However, K'gari dingoes showed significantly higher levels of very long ROH (>5 Mb), providing genomic evidence for small population size, isolation, inbreeding, and a strong founder effect. Our results suggest that, despite current levels of inbreeding, the K'gari population is purging strongly deleterious mutations, which, in the absence of further reductions in population size, may facilitate the persistence of small populations despite low genetic diversity and isolation. However, there may be little to no purging of mildly deleterious alleles, which may have important long-term consequences, and should be considered by conservation and management programs.

Dingoes, companions in life and death: The significance of archaeological canid burial practices in Australia

The dingo, also known as the Australian native dog, was introduced in the late Holocene. Dingoes were primarily wild animals but a number resided in Aboriginal people's camps. Traditionally, these individuals were taken from wild litters before weaning and raised by Aboriginal people. It is generally believed that these dingoes were not directly provided for, and upon sexual maturity, returned to reproduce in the wild. However, some died while in the company of people and, were buried in occupation sites. This Australian practice parallels the burial of domestic dogs in many regions of the Asia-Pacific and beyond but has attracted very little research. We explore the historical and archaeological evidence for dingo burial, examining its different forms, chronological and geographic distribution, and cultural significance. Dingoes were usually buried in the same manner as Aboriginal community members and often in areas used for human burial, sometimes alongside people. This practice probably occurred from the time of their introduction until soon after European colonisation. We present a case study of dingo burials from Curracurrang Rockshelter (NSW) which provides insights into the lives of ancient tame dingoes, and suggests that domestication and genetic continuity between successive camp-dwelling generations may have occurred prior to European contact.

Genome-wide variant analyses reveal new patterns of admixture and population structure in Australian dingoes

Admixture between species is a cause for concern in wildlife management. Canids are particularly vulnerable to interspecific hybridisation, and genetic admixture has shaped their evolutionary history. Microsatellite DNA testing, relying on a small number of genetic markers and geographically restricted reference populations, has identified extensive domestic dog admixture in Australian dingoes and driven conservation management policy. But there exists a concern that geographic variation in dingo genotypes could confound ancestry analyses that use a small number of genetic markers. Here, we apply genome-wide single-nucleotide polymorphism (SNP) genotyping to a set of 402 wild and captive dingoes collected from across Australia and then carry out comparisons to domestic dogs. We then perform ancestry modelling and biogeographic analyses to characterise population structure in dingoes and investigate the extent of admixture between dingoes and dogs in different regions of the continent. We show that there are at least five distinct dingo populations across Australia. We observed limited evidence of dog admixture in wild dingoes. Our work challenges previous reports regarding the occurrence and extent of dog admixture in dingoes, as our ancestry analyses show that previous assessments severely overestimate the degree of domestic dog admixture in dingo populations, particularly in south-eastern Australia. These findings strongly support the use of genome-wide SNP genotyping as a refined method for wildlife managers and policymakers to assess and inform dingo management policy and legislation moving forwards.

Small mammals reduce activity during high moon illumination under risk of predation by introduced predators

Predation influences prey survival and drives evolution of anti-predator behaviour. Anti-predator strategies by prey are stimulated by direct encounters with predators, but also by exposure to indicators of risk such as moonlight illumination and vegetation cover. Many prey species will suffer increased risk on moonlit nights, but risk may be reduced by the presence of dense vegetation. Determining the role of vegetation in reducing perceived risk is important, especially given predictions of increased global wildfire, which consumes vegetation and increases predation. We used remote cameras in southeastern Australia to compare support for the predation risk and habitat-mediated predation risk hypotheses. We examined the influence of moonlight and understorey cover on seven 20-2500 g mammalian prey species and two introduced predators, red foxes and feral cats. Activity of all prey species reduced by 40-70% with increasing moonlight, while one species (bush rat) reduced activity in response to increasing moonlight more sharply in low compared to high understorey cover. Neither predator responded to moonlight. Our findings supported the predation risk hypothesis and provided limited support for the habitat-mediated predation risk hypothesis. For prey, perceived costs of increased predation risk on moonlit nights outweighed any benefits of a brighter foraging environment.

The Australasian dingo archetype: de novo chromosome-length genome assembly, DNA methylome, and cranial morphology

BACKGROUND

One difficulty in testing the hypothesis that the Australasian dingo is a functional intermediate between wild wolves and domesticated breed dogs is that there is no reference specimen. Here we link a high-quality de novo long-read chromosomal assembly with epigenetic footprints and morphology to describe the Alpine dingo female named Cooinda. It was critical to establish an Alpine dingo reference because this ecotype occurs throughout coastal eastern Australia where the first drawings and descriptions were completed.

FINDINGS

We generated a high-quality chromosome-level reference genome assembly (Canfam_ADS) using a combination of Pacific Bioscience, Oxford Nanopore, 10X Genomics, Bionano, and Hi-C technologies. Compared to the previously published Desert dingo assembly, there are large structural rearrangements on chromosomes 11, 16, 25, and 26. Phylogenetic analyses of chromosomal data from Cooinda the Alpine dingo and 9 previously published de novo canine assemblies show dingoes are monophyletic and basal to domestic dogs. Network analyses show that the mitochondrial DNA genome clusters within the southeastern lineage, as expected for an Alpine dingo. Comparison of regulatory regions identified 2 differentially methylated regions within glucagon receptor GCGR and histone deacetylase HDAC4 genes that are unmethylated in the Alpine dingo genome but hypermethylated in the Desert dingo. Morphologic data, comprising geometric morphometric assessment of cranial morphology, place dingo Cooinda within population-level variation for Alpine dingoes. Magnetic resonance imaging of brain tissue shows she had a larger cranial capacity than a similar-sized domestic dog.

CONCLUSIONS

These combined data support the hypothesis that the dingo Cooinda fits the spectrum of genetic and morphologic characteristics typical of the Alpine ecotype. We propose that she be considered the archetype specimen for future research investigating the evolutionary history, morphology, physiology, and ecology of dingoes. The female has been taxidermically prepared and is now at the Australian Museum, Sydney.

The Australasian dingo archetype: De novo chromosome-length genome assembly, DNA methylome, and cranial morphology

Background

One difficulty in testing the hypothesis that the Australasian dingo is a functional intermediate between wild wolves and domesticated breed dogs is that there is no reference specimen. Here we link a high-quality de novo long read chromosomal assembly with epigenetic footprints and morphology to describe the Alpine dingo female named Cooinda. It was critical to establish an Alpine dingo reference because this ecotype occurs throughout coastal eastern Australia where the first drawings and descriptions were completed.

Findings

We generated a high-quality chromosome-level reference genome assembly (Canfam_ADS) using a combination of Pacific Bioscience, Oxford Nanopore, 10X Genomics, Bionano, and Hi-C technologies. Compared to the previously published Desert dingo assembly, there are large structural rearrangements on Chromosomes 11, 16, 25 and 26. Phylogenetic analyses of chromosomal data from Cooinda the Alpine dingo and nine previously published de novo canine assemblies show dingoes are monophyletic and basal to domestic dogs. Network analyses show that the mtDNA genome clusters within the southeastern lineage, as expected for an Alpine dingo. Comparison of regulatory regions identified two differentially methylated regions within glucagon receptor GCGR and histone deacetylase HDAC4 genes that are unmethylated in the Alpine dingo genome but hypermethylated in the Desert dingo. Morphological data, comprising geometric morphometric assessment of cranial morphology place dingo Cooinda within population-level variation for Alpine dingoes. Magnetic resonance imaging of brain tissue show she had a larger cranial capacity than a similar-sized domestic dog.

Conclusions

These combined data support the hypothesis that the dingo Cooinda fits the spectrum of genetic and morphological characteristics typical of the Alpine ecotype. We propose that she be considered the archetype specimen for future research investigating the evolutionary history, morphology, physiology, and ecology of dingoes. The female has been taxidermically prepared and is now at the Australian Museum, Sydney.

Genetic Structure of the Ca Rater Mallorquí Dog Breed Inferred by Microsatellite Markers

Ca Rater Mallorquí is a dog breed from the Island of Mallorca (Spain) traditionally used as a hunting and ratting dog to prevent disease spread and economic losses related to rodent activities on farms. However, the census data shows a population decline that should be addressed by implementing a conservation program. The first step to implementing a conservation plan is knowing the genetic situation of the Ca Rater Mallorquí population. Therefore, we aimed to genetically characterise the breed in our study. We analysed 33 microsatellites recommended by the International Society of Animal Genetics (ISAG) in 77 samples. Data were obtained from 13 samples of Balearic, Spanish, and international dog breeds to study the genetic diversity among breeds. The population did not significantly deviate from the Hardy-Weinberg equilibrium with heterozygosity (Ho) of 0.655 and expected heterozygosity (He) of 0.685. The Wright's fixation indices, the Factorial Correspondence Analysis (FCA), a dendrogram representing Reynolds genetic distance between populations, and the pairwise F_ST values establish the Ca Rater Mallorquí as an independent breed distinct from the Balearic, Spanish, and international breeds.

The myth of wild dogs in Australia: are there any out there?

Hybridisation between wild and domestic canids is a global conservation and management issue. In Australia, dingoes are a distinct lineage of wild-living canid with a controversial domestication status. They are mainland Australia’s apex terrestrial predator. There is ongoing concern that the identity of dingoes has been threatened from breeding with domestic dogs, and that feral dogs have established populations in rural Australia. We collate the results of microsatellite DNA testing from 5039 wild canids to explore patterns of domestic dog ancestry in dingoes and observations of feral domestic dogs across the continent. Only 31 feral dogs were detected, challenging the perception that feral dogs are widespread in Australia. First generation dingo × dog hybrids were similarly rare, with only 27 individuals identified. Spatial patterns of genetic ancestry across Australia identified that dingo populations in northern, western and central Australia were largely free from domestic dog introgression. Our findings challenge the perception that dingoes are virtually extinct in the wild and that feral dogs are common. A shift in terminology from wild dog to dingo would better reflect the identity of these wild canids and allow more nuanced debate about the balance between conservation and management of dingoes in Australia.

Genome-Wide SNPs Detect Hybridisation of Marsupial Gliders (Petaurus breviceps breviceps × Petaurus norfolcensis) in the Wild

Petaurus breviceps and Petaurus norfolcensis have produced hybrids in captivity, however there are no reported cases of Petaurus hybridisation in the wild. This study uses morphological data, mitochondrial DNA, and nuclear genome-wide SNP markers to confirm P. breviceps breviceps × P. norfolcensis hybridisation within their natural range on the central coast of New South Wales, Australia. Morphological data identified a potential hybrid that was confirmed with next-generation sequencing technology and 10,111 genome-wide SNPs. Both STRUCTURE and NewHybrid analyses identified the hybrid as a P. norfolcensis backcross, which suggests an initial F1 hybrid was fertile. The mitochondrial DNA matched that of a P. b. breviceps, indicating that a P. b. breviceps female initially mated with a P. norfolcensis male to produce a fertile female offspring. Our study is an important example of how genome-wide SNPs can be used to identify hybrids where the distribution of congeners overlaps. Hybridisation between congeners is likely to become more frequent as climate changes and habitats fragment, resulting in increased interactions and competition for resources and mates.

Conservation concerns associated with low genetic diversity for K'gari-Fraser Island dingoes

The dingo population on world heritage-listed K'gari-Fraser Island (K'gari) is amongst the most well-known in Australia. However, an absence of population genetic data limits capacity for informed conservation management. We used 9 microsatellite loci to compare the levels of genetic diversity and genetic structure of 175 K'gari dingo tissue samples with 264 samples from adjacent mainland regions. Our results demonstrated that the K'gari population has significantly lower genetic diversity than mainland dingoes (A_R, H_E, P_AR; p < 0.05) with a fourfold reduction in effective population size (N_e = 25.7 vs 103.8). There is also strong evidence of genetic differentiation between the island and mainland populations. These results are in accordance with genetic theory for small, isolated, island populations, and most likely the result of low initial diversity and founder effects such as bottlenecks leading to decreased diversity and drift. As the first study to incorporate a large sample set of K'gari dingoes, this provides invaluable baseline data for future research, which should incorporate genetic and demographic monitoring to ensure long-term persistence. Given that human-associated activities will continue to result in dingo mortality, it is critical that genetic factors are considered in conservation management decisions to avoid deleterious consequences for this iconic dingo population.

New Guinea highland wild dogs are the original New Guinea singing dogs

New Guinea singing dogs (NGSD) are identifiable by their namesake vocalizations, which are unlike any other canid population. Their novel behaviors and potential singular origin during dog domestication make them an attractive, but elusive, subject for evolutionary and conservation study. Although once plentiful on the island of New Guinea (NG), they were presumed to currently exist only in captivity. This conclusion was based on the lack of sightings in the lowlands of the island and the concurrent expansion of European- and Asian-derived dogs. We have analyzed the first nuclear genomes from a canid population discovered during a recent expedition to the highlands of NG. The extreme altitude (>4,000 m) of the highland wild dogs' (HWD) observed range and confirmed vocalizations indicate their potential to be a wild NGSD population. Comparison of single-nucleotide polymorphism genotypes shows strong similarity between HWD and the homogeneous captive NGSD, with the HWD showing significantly higher genetic diversity. Admixture analyses and estimation of shared haplotypes with phylogenetically diverse populations also indicates the HWD is a novel population within the distinct evolutionary lineage of Oceanic canids. Taken together, these data indicate the HWD possesses a distinct potential to aid in the conservation of NGSD both in the wild and under human care.

Genomic regions under selection in the feralization of the dingoes

Dingoes are wild canids living in Australia, originating from domestic dogs. They have lived isolated from both the wild and the domestic ancestor, making them a unique model for studying feralization. Here, we sequence the genomes of 10 dingoes and 2 New Guinea Singing Dogs. Phylogenetic and demographic analyses show that dingoes originate from dogs in southern East Asia, which migrated via Island Southeast Asia to reach Australia around 8300 years ago, and subsequently diverged into a genetically distinct population. Selection analysis identifies 50 positively selected genes enriched in digestion and metabolism, indicating a diet change during feralization of dingoes. Thirteen of these genes have shifted allele frequencies compared to dogs but not compared to wolves. Functional assays show that an A-to-G mutation in ARHGEF7 decreases the endogenous expression, suggesting behavioral adaptations related to the transitions in environment. Our results indicate that the feralization of the dingo induced positive selection on genomic regions correlated to neurodevelopment, metabolism and reproduction, in adaptation to a wild environment.

The Australian dingo: untamed or feral?

BACKGROUND

The Australian dingo continues to cause debate amongst Aboriginal people, pastoralists, scientists and the government in Australia. A lingering controversy is whether the dingo has been tamed and has now reverted to its ancestral wild state or whether its ancestors were domesticated and it now resides on the continent as a feral dog. The goal of this article is to place the discussion onto a theoretical framework, highlight what is currently known about dingo origins and taxonomy and then make a series of experimentally testable organismal, cellular and biochemical predictions that we propose can focus future research.

DISCUSSION

We consider a canid that has been unconsciously selected as a tamed animal and the endpoint of methodical or what we now call artificial selection as a domesticated animal. We consider wild animals that were formerly tamed as untamed and those wild animals that were formerly domesticated as feralized. Untamed canids are predicted to be marked by a signature of unconscious selection whereas feral animals are hypothesized to be marked by signatures of both unconscious and artificial selection. First, we review the movement of dingo ancestors into Australia. We then discuss how differences between taming and domestication may influence the organismal traits of skull morphometrics, brain and size, seasonal breeding, and sociability. Finally, we consider cellular and molecular level traits including hypotheses concerning the phylogenetic position of dingoes, metabolic genes that appear to be under positive selection and the potential for micronutrient compensation by the gut microbiome.

CONCLUSIONS

Western Australian Government policy is currently being revised to allow the widespread killing of the Australian dingo. These policies are based on an incomplete understanding of the evolutionary history of the canid and assume the dingo is feralized. However, accumulated evidence does not definitively show that the dingo was ever domesticated and additional focused research is required. We suggest that incorporating ancient DNA data into the debate concerning dingo origins will be pivotal to understanding the evolutionary history of the canid. Further, we advocate that future morphological, behavioural and genetic studies should focus on including genetically pure Alpine and Desert dingoes and not dingo-dog hybrids. Finally, we propose that future studies critically examine genes under selection in the dingo and employ the genome from a wild canid for comparison.