Genomic analyses of gray fox lineages suggest ancient divergence and secondary contact in the southern Great Plains

Abstract The gray fox (Urocyon cinereoargenteus) lineage diverged from all other extant canids at their most basal node and is restricted to the Americas. Previous mitochondrial analysis from coastal populations identified deeply divergent (up to 1 Mya) eastern and western lineages that predate most intraspecific splits in carnivores. We conducted genotyping by sequencing and mitochondrial analysis on gray foxes sampled across North America to determine geographic concordance between nuclear and mitochondrial contact zones and divergence times. We also estimated the admixture within the contact zone between eastern and western gray foxes based on nuclear DNA. Both datasets confirmed that eastern and western lineages met in the southern Great Plains (i.e. Texas and Oklahoma), where they maintained high differentiation. Admixture was generally low, with the majority of admixed individuals carrying