First evidence of hybridization between golden jackal (Canis aureus) and domestic dog (Canis familiaris) as revealed by genetic markers

Simulating the efficacy of wolf-dog hybridization management with individual-based modeling

Introgressive hybridization between wolves and dogs is a conservation concern due to its potentially deleterious long-term evolutionary consequences. European legislation requires that wolf-dog hybridization be mitigated through effective management. We developed an individual-based model (IBM) to simulate the life cycle of gray wolves that incorporates aspects of wolf sociality that affect hybridization rates (e.g., the dissolution of packs after the death of one/both breeders) with the goal of informing decision-making on management of wolf-dog hybridization. We applied our model by projecting hybridization dynamics in a local wolf population under different mate choice and immigration scenarios and contrasted results of removal of admixed individuals with their sterilization and release. In several scenarios, lack of management led to complete admixture, whereas reactive management interventions effectively reduced admixture in wolf populations. Management effectiveness, however, strongly depended on mate choice and number and admixture level of individuals immigrating into the wolf population. The inclusion of anthropogenic mortality affecting parental and admixed individuals (e.g., poaching) increased the probability of pack dissolution and thus increased the probability of interbreeding with dogs or admixed individuals and boosted hybridization and introgression rates in all simulation scenarios. Recognizing the necessity of additional model refinements (appropriate parameterization, thorough sensitivity analyses, and robust model validation) to generate management recommendations applicable in real-world scenarios, we maintain confidence in our model's potential as a valuable conservation tool that can be applied to diverse situations and species facing similar threats.

Detection and complete genome characterization of a genogroup X (GX) sapovirus (family Caliciviridae) from a golden jackal (Canis aureus) in Hungary

In this study, a novel genotype of genogroup X (GX) sapovirus (family Caliciviridae) was detected in the small intestinal contents of a golden jackal (Canis aureus) in Hungary and characterised by viral metagenomics and next-generation sequencing techniques. The complete genome of the detected strain, GX/Dömsöd/DOCA-11/2020/HUN (PP105600), is 7,128 nt in length. The ORF1- and ORF2-encoded viral proteins (NSP, VP1, and VP2) have 98%, 95%, and 88% amino acid sequence identity to the corresponding proteins of genogroup GX sapoviruses from domestic pigs, but the nucleic acid sequence identity values for their genes are significantly lower (83%, 77%, and 68%). During an RT-PCR-based epidemiological investigation of additional jackal and swine samples, no other GX strains were detected, but a GXI sapovirus strain, GXI/Tótfalu/WBTF-10/2012/HUN (PP105601), was identified in a faecal sample from a wild boar (Sus scrofa). We report the detection of members of two likely underdiagnosed groups of sapoviruses (GX and GXI) in a golden jackal and, serendipitously, in a wild boar in Europe.

Protoparvovirus carnivoran 1 infection of golden jackals Canis aureus in Serbia

Parvoviruses are among the major animal pathogens that can cause considerable health disorders ranging from subclinical to lethal in domestic and wild animals. Golden jackal (Canis aureus), an expanding European species, is a reservoir of many pathogens, including vector-borne diseases and zoonoses. Given the importance of parvovirus infections in dogs and cats, this study aimed to unfold the virus prevalence and molecular characterisation in the golden jackal population in Serbia. The spleen samples from 68 hunted jackals during 2022/2023 were tested for the VP2-specific genome region of Protoparvovirus carnivoran 1 by PCR. BLAST analysis of partial VP2 sequences obtained from three animals (4.4%) revealed the highest similarity to Protoparvovirus carnivoran 1, genogroup Feline panleukopenia virus, which is the second report on FPV infection in jackals. Based on specific amino acid residues within partial VP2, the jackals' Protoparvovirus carnivoran 1 was also classified as FPV. One jackal's strain showed two synonymous mutations at positions 699 and 1167. Although species cross-transmission could not be established, jackals' health should be maintained by preventing the transmission of viruses to native species and vice versa. Although jackals are considered pests, their role as natural cleaners is of greater importance. Therefore, further monitoring of their health is needed to understand the influence of infectious diseases on population dynamics and to determine the relationship between domestic predators and jackals and the direction of cross-species transmission.

Recent selection created distinctive variability patterns on MHC class II loci in three dolphin species from the Mediterranean Sea

The major histocompatibility complex (MHC) includes highly polymorphic genes involved in antigen presentation, which is crucial for adaptive immune response. They represent fitness related genetic markers particularly informative for populations exposed to environmental challenges. Here we analyse the diversity and evolutionary traits of MHC class II DQA and DQB genes in the dolphins Stenella coeruleoalba and Grampus griseus from the Mediterranean Sea. We found substantial nucleotide and functional diversity, as well as strong evidence of balancing selection indicated by allele and supertype frequencies, Tajima's D statistics and dN/dS tests. The Risso's dolphin, considered the least abundant in the region, showed the effect of divergent allele advantage at the nucleotide and functional-peptide levels. An outstanding polymorphism was found in the striped dolphin, particularly intriguing in the DQA gene where the Ewens-Watterson test detected a selection sweep that occurred in recent history. We hypothesize that morbillivirus, which has recurrently invaded Mediterranean populations over the last decades, exerted the detected selective pressure.

White and other fur colourations and hybridization in golden jackals (Canis aureus) in the Carpathian basin

The golden jackal (Canis aureus) is a reoccurring species in the centre of the Carpathian basin, in Hungary. In total, 31 golden jackal tissue samples were collected, from 8 white-coated, 2 black-coated and one mottled animal across Hungary. Sequences and fragment length polymorphisms were studied for white colour (MC1R), and for black coat colouration (CBD103). In each white animal, the most widespread mutation causing white fur colour in dogs in homozygous form was detected. Three animals were found to carry the mutation in heterozygous form. The two black golden jackals were heterozygous for the 3 bp deletion in CBD103 that mutation for black coat colouration in dogs, and one of them also carried the mutation causing white fur. None of the white animals showed signs of hybridization, but both the black and the mottled coloured individuals were found to be hybrids based on genetic testing. Kinship was found three times, twice between white animals, and once between a white animal and an agouti animal carrying the mutation of white coat. Our results confirm the findings that golden jackal-dog hybrids may occur without human intervention, and the detected mutation causing white fur colour in golden jackals could possibly be due to an early hybridization event.

Mid-Pleistocene Transitions Forced Himalayan ibex to Evolve Independently after Split into an Allopatric Refugium

Pleistocene glaciations had profound impact on the spatial distribution and genetic makeup of species in temperate ecosystems. While the glacial period trapped several species into glacial refugia and caused abrupt decline in large populations, the interglacial period facilitated population growth and range expansion leading to allopatric speciation. Here, we analyzed 40 genomes of four species of ibex and found that Himalayan ibex in the Pamir Mountains evolved independently after splitting from its main range about 0.1 mya following the Pleistocene species pump concept. Demographic trajectories showed Himalayan ibex experienced two historic bottlenecks, one each c. 0.8-0.5 mya and c. 50-30 kya, with an intermediate large population expansion c. 0.2-0.16 mya coinciding with Mid-Pleistocene Transitions. We substantiate with multi-dimensional evidence that Himalayan ibex is an evolutionary distinct phylogenetic species of Siberian ibex which need to be prioritized as Capra himalayensis for taxonomic revision and conservation planning at a regional and global scale.

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.

Possible origins and implications of atypical morphologies and domestication-like traits in wild golden jackals (Canis aureus)

Deciphering the origins of phenotypic variations in natural animal populations is a challenging topic for evolutionary and conservation biologists. Atypical morphologies in mammals are usually attributed to interspecific hybridisation or de-novo mutations. Here we report the case of four golden jackals (Canis aureus), that were observed during a camera-trapping wildlife survey in Northern Israel, displaying anomalous morphological traits, such as white patches, an upturned tail, and long thick fur which resemble features of domesticated mammals. Another individual was culled under permit and was genetically and morphologically examined. Paternal and nuclear genetic profiles, as well as geometric morphometric data, identified this individual as a golden jackal rather than a recent dog/wolf-jackal hybrid. Its maternal haplotype suggested past introgression of African wolf (Canis lupaster) mitochondrial DNA, as previously documented in other jackals from Israel. When viewed in the context of the jackal as an overabundant species in Israel, the rural nature of the surveyed area, the abundance of anthropogenic waste, and molecular and morphological findings, the possibility of an individual presenting incipient stages of domestication should also be considered.

The computational implementation of a platform of relative identity-by-descent scores algorithm for introgressive mapping

With the development of genotyping and sequencing technology, researchers working in the area of conservation genetics are able to obtain the genotypes or even the sequences of a representative sample of individuals from the population. It is of great importance to examine the genomic variants and genes that are highly preferred or pruned during the process of adaptive introgression or long-term hybridization. To the best of our knowledge, we are the first to develop a platform with computational integration of a relative identity-by-descent (rIBD) scores algorithm for introgressive mapping. The rIBD algorithm is designed for mapping the fine-scaled genomic regions under adaptive introgression between the source breeds and the admixed breed. Our rIBD calculation platform provides compact functions including reading input information and uploading of files, rIBD calculation, and presentation of the rIBD scores. We analyzed the simulated data using the rIBD calculation platform and calculated the average IBD score of 0.061 with a standard deviation of 0.124. The rIBD scores generally follow a normal distribution, and a cut-off of 0.432 and -0.310 for both positive and negative rIBD scores is derived to enable the identification of genomic regions showing significant introgression signals from the source breed to the admixed breed. A list of genomic regions with detailed calculated rIBD scores is reported, and all the rIBD scores for each of the considered windows are presented in plots on the rIBD calculation platform. Our rIBD calculation platform provides a user-friendly tool for the calculation of fine-scaled rIBD scores for each of the genomic regions to map possible functional genomic variants due to adaptive introgression or long-term hybridization.

African or Asiatic origin? Genome analyses solved the mystery of the hybrid origin of the rescued lion cub

Species assignment of any seized material using DNA analysis has been a routine and widely accepted standard procedure in providing scientific advisory for the legal prosecution of wildlife cases. Scientific advancements and rigorous application of genetic tools have led to the development of a variety of molecular markers with their defined efficacy in wildlife forensics. However, in a few unusual cases where a hybrid needs to be identified or assignment need to be made at sub-species level, mitochondrial markers often fail or else provide biased results, which can affect the overall judgment in the court of law. Here, we report one such challenging case of lion cub rescued by the law enforcement from illegal trafficking. Phylogenetic assessment based on complete mitogenome assigned rescued lion cub with African lion (Panthera leo leo). However, the TSPY gene of the Y chromosome established that the lion cub shared its paternal lineage from Asiatic lion (Panthera leo persica). With the use of maternally and paternally inherited markers, we conclude a hybrid origin of the rescued lion cub which shared ancestry from both Asiatic as well as African lion. The present study exhibits the application of genome sequencing in thinking beyond routine identification and contributes to the operating procedures of wildlife forensics.

A reduced SNP panel to trace gene flow across southern European wolf populations and detect hybridization with other Canis taxa

Intra- and inter-specific gene flow are natural evolutionary processes. However, human-induced hybridization is a global conservation concern across taxa, and the development of discriminant genetic markers to differentiate among gene flow processes is essential. Wolves (Canis lupus) are affected by hybridization, particularly in southern Europe, where ongoing recolonization of historic ranges is augmenting gene flow among divergent populations. Our aim was to provide diagnostic canid markers focused on the long-divergent Iberian, Italian and Dinaric wolf populations, based on existing genomic resources. We used 158 canid samples to select a panel of highly informative single nucleotide polymorphisms (SNPs) to (i) distinguish wolves in the three regions from domestic dogs (C. l. familiaris) and golden jackals (C. aureus), and (ii) identify their first two hybrid generations. The resulting 192 SNPs correctly identified the five canid groups, all simulated first-generation (F1) hybrids (0.482 ≤ Q_i ≤ 0.512 between their respective parental groups) and all first backcross (BC1) individuals (0.723 ≤ Q_i ≤ 0.827 to parental groups). An assay design and test with invasive and non-invasive canid samples performed successfully for 178 SNPs. By separating natural population admixture from inter-specific hybridization, our reduced panel can help advance evolutionary research, monitoring, and timely conservation management.

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.

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.

Reliable wolf-dog hybrid detection in Europe using a reduced SNP panel developed for non-invasively collected samples

Background

Understanding the processes that lead to hybridization of wolves and dogs is of scientific and management importance, particularly over large geographical scales, as wolves can disperse great distances. However, a method to efficiently detect hybrids in routine wolf monitoring is lacking. Microsatellites offer only limited resolution due to the low number of markers showing distinctive allele frequencies between wolves and dogs. Moreover, calibration across laboratories is time-consuming and costly. In this study, we selected a panel of 96 ancestry informative markers for wolves and dogs, derived from the Illumina CanineHD Whole-Genome BeadChip (174 K). We designed very short amplicons for genotyping on a microfluidic array, thus making the method suitable also for non-invasively collected samples.

Results

Genotypes based on 93 SNPs from wolves sampled throughout Europe, purebred and non-pedigree dogs, and suspected hybrids showed that the new panel accurately identifies parental individuals, first-generation hybrids and first-generation backcrosses to wolves, while second- and third-generation backcrosses to wolves were identified as advanced hybrids in almost all cases. Our results support the hybrid identity of suspect individuals and the non-hybrid status of individuals regarded as wolves. We also show the adequacy of these markers to assess hybridization at a European-wide scale and the importance of including samples from reference populations.

Conclusions

We showed that the proposed SNP panel is an efficient tool for detecting hybrids up to the third-generation backcrosses to wolves across Europe. Notably, the proposed genotyping method is suitable for a variety of samples, including non-invasive and museum samples, making this panel useful for wolf-dog hybrid assessments and wolf monitoring at both continental and different temporal scales.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07761-5.

Positive selection on the MHC class II DLA-DQA1 gene in golden jackals (Canis aureus) from their recent expansion range in Europe and its effect on their body mass index

Background

In Europe, golden jackals (Canis aureus) have been expanding their range out of the southern and southeastern Balkans towards central Europe continually since the 1960s. Here, we investigated the level of functional diversity at the MHC class II DLA-DQA1 exon 2 in golden jackal populations from Bulgaria, Serbia, and Hungary. Specifically, we tested for positive selection on and geographic variation at that locus due to adaptation to supposedly regionally varying pathogenic landscapes. To test for potential fitness effects of different protein variants on individual body condition, we used linear modeling of individual body mass indexes (bmi) and accounted for possible age, sex, geographical, and climatic effects. The latter approach was performed, however, only on Serbian individuals with appropriate data.

Results

Only three different DLA-DQA1 alleles were detected, all coding for different amino-acid sequences. The neutrality tests revealed no significant but positive values; there was no signal of spatial structuring and no deviation from the Hardy–Weinberg equilibrium across the studied range of expansion. However, we found a signal of trans-species polymorphism and significant test results for positive selection on three codons. Our information-theory based linear modeling results indicated an effect of ambient temperature on the occurrence of individual DLA-DQA1 genotypes in individuals from across the studied expansion range, independent from geographical position. Our linear modeling results of individual bmi values indicated that yearlings homozygous for DLA-DQA1*03001 reached values typical for adults contrary to yearlings carrying other genotypes (protein combinations). This suggested better growth rates and thus a possible fitness advantage of yearlings homozygous for DLA-DQA1*03001.

Conclusions

Our results indicate a demographic (stochastic) signal of reduced DLA-DQA1 exon 2 variation, in line with the documented historical demographic bottleneck. At the same time, however, allelic variation was also affected by positive selection and adaptation to varying ambient temperature, supposedly reflecting geographic variation in the pathogenic landscape. Moreover, an allele effect on body mass index values of yearlings suggested differential fitness associated with growth rates. Overall, a combination of a stochastic effect and positive selection has shaped and is still shaping the variation at the studied MHC locus.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01856-z.

Updating the Bibliography of Interbreeding among Canis in North America

This bibliography provides a collection of references that documents the evolution of studies evidencing interbreeding among Canis species in North America. Over the past several decades, advances in biology and genomic technology greatly improved our ability to detect and characterize species interbreeding, which has significance for understanding species in a changing landscape as well as for endangered species management. This bibliography includes a discussion within each category of interbreeding, the timeline of developing evidence, and includes a review of past research conducted on experimental crosses. Research conducted in the early 20th century is rich with detailed records and photographs of hybrid offspring development and behavior. With the progression of molecular methods, studies can estimate historical demographic parameters and detect chromosomal patterns of ancestry. As these methods continue to increase in accessibility, the field will gain a deeper and richer understanding of the evolutionary history of North American Canis.

Adding Nuance to Our Understanding of Dog-Wildlife Interactions and the Need for Management

The interactions of dogs and wildlife are receiving increasing attention across the globe. Evidence suggests that dogs can negatively affect wildlife in a variety of ways, including through the risk of predation, by hybridizing with other Canis species, by acting as a reservoir or vector of pathogens, and by competing with wildlife for resources. A multitude of observations and case studies for each of these interactions has led to calls to prioritize increased management of dogs, for instance, through population reduction, vaccination, movement restrictions, and enhanced care of the dogs themselves. Here, I argue that while the risks that dogs present to wildlife are real, they vary in their importance across the globe. Furthermore, often the approaches used to address these risks are championed without a full understanding of the likelihood of success given the necessary spatial scale of management, the availability of alternative management approaches, and the role dogs play in societies. For instance, culling, vaccination, and animal husbandry approaches in reducing the impacts of dogs on wildlife sometimes fail to recognize that local human populations actively recruit replacement dogs, that vaccination often fails to reach levels necessary for herd immunity, and that enhanced dog husbandry may have indirect impacts on wildlife by requiring environmentally impactful activities. I suggest there is a need for attention not only to identify the impacts of dogs but also to determine where and when such interactions with wildlife are problematic and the likelihood of success for any proposed management approach. The impacts of dogs should be mediated in a context-specific manner that accounts for factors such as the local density of dogs, the susceptibility of local wildlife populations to the risks that dog populations may represent, and the local societal norms that underpin how dogs might be managed.

Non-invasive diagnostic PCRs for rapid detection of golden jackal, red fox, and gray wolf/domestic dog and application to validate golden jackal presence in Styria, Austria

Human-predator conflicts are frequently caused by livestock and/or game depredation. The golden jackal’s (Canis aureus) range expansion in Europe, as well as the recent re-expansion of several gray wolf (Canis lupus) populations, might increase risk of such conflicts. In Austria, golden jackal presence has been increasing since the 1990s including reports of wildlife and livestock kills, frequently occurring in the provinces Styria and Burgenland. We developed a rapid, two-step genetic screening protocol to (1) detect canid mtDNA from non-invasively collected samples like swabs from kills using diagnostic PCRs, and (2) assign this DNA to red fox (Vulpes vulpes), golden jackal, or gray wolf/dog. To monitor golden jackal presence in the region, a total of 167 signs of presence were collected over a period of 30 months throughout the Styrian province. Among these, 14 non-invasive genetic samples (13 swabs from kill sites and one scat) were screened with the developed protocol. Four of these samples revealed golden jackal mtDNA and six samples red fox mtDNA. The developed genetic screening protocol represents a quick and inexpensive method to assess canid presence, e.g., at kill sites, and therefore possesses high value for the conservation and wildlife management community.

Metabolomics shows the Australian dingo has a unique plasma profile

Dingoes occupy a wide range of the Australian mainland and play a crucial role as an apex predator with a generalist omnivorous feeding behaviour. Dingoes are ecologically, phenotypically and behaviourally distinct from modern breed dogs and have not undergone artificial selection since their arrival in Australia. In contrast, humans have selected breed dogs for novel and desirable traits. First, we examine whether the distinct evolutionary histories of dingoes and domestic dogs has lead to differences in plasma metabolomes. We study metabolite composition differences between dingoes (n = 15) and two domestic dog breeds (Basenji n = 9 and German Shepherd Dog (GSD) n = 10). Liquid chromatography mass spectrometry, type II and type III ANOVA with post-hoc tests and adjustments for multiple comparisons were used for data evaluation. After accounting for within group variation, 62 significant metabolite differences were detected between dingoes and domestic dogs, with the majority of differences in protein (n = 14) and lipid metabolites (n = 12), mostly lower in dingoes. Most differences were observed between dingoes and domestic dogs and fewest between the domestic dog breeds. Next, we collect a second set of data to investigate variation between pure dingoes (n = 10) and dingo-dog hybrids (n = 10) as hybridisation is common in regional Australia. We detected no significant metabolite differences between dingoes and dingo-dog hybrids after Bonferroni correction. However, power analysis showed that increasing the sample size to 15 could result in differences in uridine 5′-diphosphogalactose (UDPgal) levels related to galactose metabolism. We suggest this may be linked to an increase in Amylase 2B copy number in hybrids. Our study illustrates that the dingo metabolome is significantly different from domestic dog breeds and hybridisation is likely to influence carbohydrate metabolism.

High-quality carnivoran genomes from roadkill samples enable comparative species delineation in aardwolf and bat-eared fox

In a context of ongoing biodiversity erosion, obtaining genomic resources from wildlife is essential for conservation. The thousands of yearly mammalian roadkill provide a useful source material for genomic surveys. To illustrate the potential of this underexploited resource, we used roadkill samples to study the genomic diversity of the bat-eared fox (Otocyon megalotis) and the aardwolf (Proteles cristatus), both having subspecies with similar disjunct distributions in Eastern and Southern Africa. First, we obtained reference genomes with high contiguity and gene completeness by combining Nanopore long reads and Illumina short reads. Then, we showed that the two subspecies of aardwolf might warrant species status (P. cristatus and P. septentrionalis) by comparing their genome-wide genetic differentiation to pairs of well-defined species across Carnivora with a new Genetic Differentiation index (GDI) based on only a few resequenced individuals. Finally, we obtained a genome-scale Carnivora phylogeny including the new aardwolf species.

Wildlife forensics: A boon for species identification and conservation implications

Wildlife trade and fraudulence in food, artefacts and cosmetic industries had raised serious concern in protection of the wild faunal diversity. Lack of proper tools and molecular based techniques for identification of wild species are some of the major constrains faced by the judiciary and law enforcement agencies while framing charges against poachers and illicit agitator. The emergence of wildlife forensics serves as a boon in solving long pending cases of wildlife crimes. Wildlife forensics have proven to be fast, accurate and reliable criminal investigation processes with comprehensive coverage and easy accessibility. It has also helped resolving taxonomic disputes, determining spatiotemporal genetic divergence, evolutionary history, origins and even endemism. Collaboration among inter-disciplinary fields has even led to engineered signature markers and phylogenetics for several species. Development in fields of genetics, molecular and evolutionary biology and other omics techniques have further contributed in accurate identification of species. Wildlife forensics, with the support of proper international mega database units for population reference, will be fundamental in wildlife investigations through its unlimited information sharing ability. The efficient conservation of species will, however, require a collaborative approach consisting of national policy makers, local stakeholders and implementation agencies in addition to experts from the scientific communities.

The relevance of gene flow with wild relatives in understanding the domestication process

The widespread use of genomic tools has allowed for a deeper understanding of the genetics and the evolutionary dynamics of domestication. Recent studies have suggested that multiple domestications and introgression are more common than previously thought. However, the ability to correctly infer the many aspects of domestication process depends on having an adequate representation of wild relatives. Cultivated maize (Zea mays ssp. mays) is one of the most important crops in the world, with a long and a relatively well-documented history of domestication. The current consensus points towards a single domestication event from teosinte Zea mays ssp. parviglumis from the Balsas Basin in Southwestern Mexico. However, the underlying diversity of teosintes from Z. mays ssp. parviglumis and Zea mays ssp. mexicana was not taken into account in early studies. We used 32 739 single nucleotide polymorphisms (SNPs) obtained from 29 teosinte populations and 43 maize landraces to explore the relationship between wild and cultivated members of Zea. We then inferred the levels of gene flow among teosinte populations and maize, the degree of population structure of Zea mays subspecies, and the potential domestication location of maize. We confirmed a strong geographic structure within Z. mays ssp. parviglumis and documented multiple gene flow events with other members of the genus, including an event between Z. mays ssp. mexicana and maize. Our results suggest that the likely ancestor of maize may have been domesticated in Jalisco or in the southern Pacific Coast and not in the Balsas Basin as previously thought. In this context, different populations of both teosinte subspecies have contributed to modern maize's gene pool. Our results point towards a long period of domestication marked by gene flow with wild relatives, confirming domestication as long and ongoing process.

Matching STR and SNP genotyping to discriminate between wild boar, domestic pigs and their recent hybrids for forensic purposes

The genetic discrimination between phylogenetically close taxa can be challenging if their gene pools are not differentiated and there are many shared polymorphisms. The gene flow between wild boar (Sus scrofa) and domestic pig (S. s. domesticus) has never been interrupted from domestication onwards, due to non-stop natural and human-mediated crossbreeding. To date there are no individual genetic markers that are able to distinguish between the two forms, nor even to identify effectively their hybrids. We developed a combined molecular protocol based on multiplex porcine-specific STR-profiling system and new real time PCR-based assays of single polymorphisms in the NR6A1 and MC1R genes to gain high diagnostic power in the differentiation of wild boar, pig and hybrids for forensic purposes. The combined approach correctly assigned individuals to one or the other parental gene pool and identified admixed genotypes. Evidence was found for substantial reduction of false negative results by using multiple marker systems jointly, compared to their use individually. Our protocol is a powerful and cost-effective diagnostic tool that can easily be adopted by most forensic laboratories to assist authorities contrast food adulteration, assure veterinary public health and fight against wildlife crimes, like poaching and illegal detention of wild animals.

Range-wide patterns of human-mediated hybridisation in European wildcats

Hybridisation between wild taxa and their domestic congeners is a significant conservation issue. Domestic species frequently outnumber their wild relatives in population size and distribution and may therefore genetically swamp the native species. The European wildcat (Felis silvestris) has been shown to hybridise with domestic cats (Felis catus). Previously suggested spatially divergent introgression levels have not been confirmed on a European scale due to significant differences in the applied methods to assess hybridisation of the European wildcat. We analysed 926 Felis spp. samples from 13 European countries, using a set of 86 selected ancestry-informative SNPs, 14 microsatellites, and ten mitochondrial and Y-chromosome markers to study regional hybridisation and introgression patterns and population differentiation. We detected 51 hybrids (four F1 and 47 F2 or backcrosses) and 521 pure wildcats throughout Europe. The abundance of hybrids varied considerably among studied populations. All samples from Scotland were identified as F2 hybrids or backcrosses, supporting previous findings that the genetic integrity of that wildcat population has been seriously compromised. In other European populations, low to moderate levels of hybridisation were found, with the lowest levels being in Central and Southeast Europe. The occurrence of distinct maternal and paternal markers between wildcat and domestic cat suggests that there were no severe hybridisation episodes in the past. The overall low (< 1%) prevalence of F1 hybrids suggests a low risk of hybridisation for the long-term genetic integrity of the wildcat in most of Europe. However, regionally elevated introgression rates confirm that hybridisation poses a potential threat. We propose regional in-depth monitoring of hybridisation rates to identify factors driving hybridisation so as to develop effective strategies for conservation.

Fast speciations and slow genes: uncovering the root of living canids

Despite ongoing efforts relying on computationally intensive tree-building methods and large datasets, the deeper phylogenetic relationships between living canid genera remain controversial. We demonstrate that this issue arises fundamentally from the uncertainty of root placement as a consequence of the short length of the branch connecting the major canid clades, which probably resulted from a fast radiation during the early diversification of extant Canidae. Using both nuclear and mitochondrial genes, we investigate the position of the canid root and its consistency by using three rooting methods. We find that mitochondrial genomes consistently retrieve a root node separating the tribe Canini from the remaining canids, whereas nuclear data mostly recover a root that places the Urocyon foxes as the sister lineage of living canids. We demonstrate that, to resolve the canid root, the nuclear segments sequenced so far are significantly less informative than mitochondrial genomes. We also propose that short intervals between speciations obscure the place of the true root, because methods are susceptible to stochastic error in the presence of short internal branches near the root.

Dispersal history of the golden jackal (Canis aureus moreoticus Geoffroy, 1835) in Europe and possible causes of its recent population explosion

Background

Data on the historical distribution of the golden jackal in Europe and its primary habitats are scarce. There are many new data on the population explosion and the rapid spread of the in Europe. However, the main factors for this expansion, the core population and its routes of dispersal, remain controversial or insufficiently studied.

New information

This study provides a profound analysis of the history of the jackal’s (Canisaureusmoreoticus Geoffroy, 1835) occurrence in Europe, the factors limiting or those triggering its expansion on the continent. The analysis shows that the timing of the species appearance in Europe still remains unclear. Historical data show that the species is a typical inhabitant of South-Eastern Europe, with some pulsations within its core area, as well as extensions to the north and west of it in favourable periods. Nowadays, the increase of the species range in Europe is the largest documented population explosion on the continent. We argue that this expansion originates from only three core populations, the Peri-Strandja area and the Dalmatian coast in the Balkans and the east parts of Western Transcaucasia in the Caucasus. This population explosion is largely due to a unique combination of factors of an anthropogenic nature.

Interspecific Gene Flow Shaped the Evolution of the Genus Canis

The evolutionary history of the wolf-like canids of the genus Canis has been heavily debated, especially regarding the number of distinct species and their relationships at the population and species level [1-6]. We assembled a dataset of 48 resequenced genomes spanning all members of the genus Canis except the black-backed and side-striped jackals, encompassing the global diversity of seven extant canid lineages. This includes eight new genomes, including the first resequenced Ethiopian wolf (Canis simensis), one dhole (Cuon alpinus), two East African hunting dogs (Lycaon pictus), two Eurasian golden jackals (Canis aureus), and two Middle Eastern gray wolves (Canis lupus). The relationships between the Ethiopian wolf, African golden wolf, and golden jackal were resolved. We highlight the role of interspecific hybridization in the evolution of this charismatic group. Specifically, we find gene flow between the ancestors of the dhole and African hunting dog and admixture between the gray wolf, coyote (Canis latrans), golden jackal, and African golden wolf. Additionally, we report gene flow from gray and Ethiopian wolves to the African golden wolf, suggesting that the African golden wolf originated through hybridization between these species. Finally, we hypothesize that coyotes and gray wolves carry genetic material derived from a "ghost" basal canid lineage.

Human-Mediated Introgression of Haplotypes in a Modern Dairy Cattle Breed

Domestic animals can serve as model systems of adaptive introgression and their genomic signatures. In part, their usefulness as model systems is due to their well-known histories. Different breeding strategies such as introgression and artificial selection have generated numerous desirable phenotypes and superior performance in domestic animals. The modern Danish Red Dairy Cattle is studied as an example of an introgressed population. It originates from crossing the traditional Danish Red Dairy Cattle with the Holstein and Brown Swiss breeds, both known for high milk production. This crossing happened, among other things due to changes in the production system, to raise milk production and overall performance. The genomes of modern Danish Red Dairy Cattle are heavily influenced by regions introgressed from the Holstein and Brown Swiss breeds and under subsequent selection in the admixed population. The introgressed proportion of the genome was found to be highly variable across the genome. Haplotypes introgressed from Holstein and Brown Swiss contained or overlapped known genes affecting milk production, as well as protein and fat content (CD14, ZNF215, BCL2L12, and THRSP for Holstein origin and ITPR2, BCAT1, LAP3, and MED28 for Brown Swiss origin). Genomic regions with high introgression signals also contained genes and enriched QTL associated with calving traits, body confirmation, feed efficiency, carcass, and fertility traits. These introgressed signals with relative identity-by-descent scores larger than the median showing Holstein or Brown Swiss introgression are mostly significantly correlated with the corresponding test statistics from signatures of selection analyses in modern Danish Red Dairy Cattle. Meanwhile, the putative significant introgressed signals have a significant dependency with the putative significant signals from signatures of selection analyses. Artificial selection has played an important role in the genomic footprints of introgression in the genome of modern Danish Red Dairy Cattle. Our study on a modern cattle breed contributes to an understanding of genomic consequences of selective introgression by demonstrating the extent to which adaptive effects contribute to shape the specific genomic consequences of introgression.

Canis aureus (Carnivore: Canidae)

Canis aureus (Linnaeus, 1758), the golden jackal, is a medium-sized, wide spread, terrestrial carnivore. It is 1 of 7 species found in the genus Canis. It ranges from Africa to Europe, the Middle East, Central Asia, and Southeast Asia. Due to its tolerance of dry habitats and its omnivorous diet, C. aureus can live in a wide variety of habitats. It normally lives in open grassland habitat but also occurs in deserts, woodlands, mangroves, and agricultural and rural habitats in India and Bangladesh. It ranges from sea level in Eritrea to 3,500 m in the Bale Mountains of Ethiopia and 2,000 m in India. C. aureus is listed as “Least Concern” by the International Union for Conservation of Nature and Natural Resources Red List of Threatened Species version 2016.1.

Disentangling Timing of Admixture, Patterns of Introgression, and Phenotypic Indicators in a Hybridizing Wolf Population

Hybridization is a natural or anthropogenic process that can deeply affect the genetic make-up of populations, possibly decreasing individual fitness but sometimes favoring local adaptations. The population of Italian wolves (Canis lupus), after protracted demographic declines and isolation, is currently expanding in anthropic areas, with documented cases of hybridization with stray domestic dogs. However, identifying admixture patterns in deeply introgressed populations is far from trivial. In this study, we used a panel of 170,000 SNPs analyzed with multivariate, Bayesian and local ancestry reconstruction methods to identify hybrids, estimate their ancestry proportions and timing since admixture. Moreover, we carried out preliminary genotype-phenotype association analyses to identify the genetic bases of three phenotypic traits (black coat, white claws, and spur on the hind legs) putative indicators of hybridization. Results showed no sharp subdivisions between nonadmixed wolves and hybrids, indicating that recurrent hybridization and deep introgression might have started mostly at the beginning of the population reexpansion. In hybrids, we identified a number of genomic regions with excess of ancestry in one of the parental populations, and regions with excess or resistance to introgression compared with neutral expectations. The three morphological traits showed significant genotype-phenotype associations, with a single genomic region for black coats and white claws, and with multiple genomic regions for the spur. In all cases the associated haplotypes were likely derived from dogs. In conclusion, we show that the use of multiple genome-wide ancestry reconstructions allows clarifying the admixture dynamics even in highly introgressed populations, and supports their conservation management.

Eurasian golden jackal as host of canine vector-borne protists

Background

Jackals are medium-sized canids from the wolf-like clade, exhibiting a unique combination of ancestral morphotypes, broad trophic niches, and close phylogenetic relationships with the wolf and dog. Thus, they represent a potential host of several pathogens with diverse transmission routes. Recently, populations of the Eurasian golden jackal Canis aureus have expanded into the Western Palaearctic, including most of Europe. The aim of our study was to examine Eurasian golden jackals from Romania, Czech Republic and Austria for a wide spectrum of vector-borne protists and to evaluate the role of this species as a reservoir of disease for domestic dogs and/or humans.

Results

Diagnostic polymerase chain reaction (PCR) DNA amplifications revealed 70% of jackals to be positive for Hepatozoon, 12.5% positive for piroplasms, and one individual positive for Leishmania infantum. Phylogenetic analyses of partial 18S rDNA sequences invariably placed sequenced isolates of Hepatozoon into the H. canis clade. For piroplasms, both the 18S and cox1 sequences obtained confirmed the presence of Babesia canis and “Theileria annae” in 5 and 2 individuals, respectively, providing the first records of these two piroplasmids in Eurasian golden jackals. A single animal from Dolj County (Romania) was PCR-positive for L. infantum, as confirmed also by sequencing of ITS1-5.8S.

Conclusions

Apparently, expanding populations of jackals can play a significant role in spreading and maintaining new Babesia canis foci in Central Europe. The role of jackals in the epidemiology of “Theileria annae” and H. canis is probably similar to that of red foxes and should be taken into account in further research on these parasites. Also the presence of L. infantum deserves attention. Our study confirms that once established, the populations of Eurasian golden jackals constitute natural reservoirs for many canine vector-borne diseases, analogous to the role of the coyotes in North America.

Golden jackal (Canis aureus) in the Czech Republic: the first record of a live animal and its long-term persistence in the colonized habitat

A golden jackal (Canis aureus) individual was recorded ~40 km east of Prague in the Czech Republic. It is the first record of a living golden jackal in the country; up to now several individuals have been recorded but all of them were either shot dead or killed by a vehicle. The observed animal was documented by camera traps set up for research of carnivore diversity in different habitats in the study area. It was first photographed on 19 June 2015, and in total there were 57 records made by 12 traps until 24 March 2016 when the animal was still present in the area. Forty-nine of the 57 records were made in a shrubby grassland over an area of ~100 ha, 39% of sightings were during the day and 61% in the night. There were two distinct peaks in the circadian activity of the animal, from 4 to 10 a.m., and from 6 p.m. to midnight. We also review the verified records of the golden jackal in the Czech Republic, some of which were only published in local hunting magazines. However, the observation reported in this paper represents the first evidence of a long-term occurrence in Europe of the same golden jackal individual, that persisted for at least nine months and over winter, northwest of Hungarian-Austrian border where the population has been known to reproduce.

Population structure and gene flow in the global pest, Helicoverpa armigera

Helicoverpa armigera is a major agricultural pest that is distributed across Europe, Asia, Africa and Australasia. This species is hypothesized to have spread to the Americas 1.5 million years ago, founding a population that is at present, a distinct species, Helicoverpa zea. In 2013, H. armigera was confirmed to have re-entered South America via Brazil and subsequently spread. The source of the recent incursion is unknown and population structure in H. armigera is poorly resolved, but a basic understanding would highlight potential biosecurity failures and determine the recent evolutionary history of region-specific lineages. Here, we integrate several end points derived from high-throughput sequencing to assess gene flow in H. armigera and H. zea from populations across six continents. We first assemble mitochondrial genomes to demonstrate the phylogenetic relationship of H. armigera with other Heliothine species and the lack of distinction between populations. We subsequently use de novo genotyping-by-sequencing and whole-genome sequences aligned to bacterial artificial chromosomes, to assess levels of admixture. Primarily, we find that Brazilian H. armigera are derived from diverse source populations, with strong signals of gene flow from European populations, as well as prevalent signals of Asian and African ancestry. We also demonstrate a potential field-caught hybrid between H. armigera and H. zea, and are able to provide genomic support for the presence of the H. armigera conferta subspecies in Australasia. While structure among the bulk of populations remains unresolved, we present distinctions that are pertinent to future investigations as well as to the biosecurity threat posed by H. armigera.