European wildcat populations are subdivided into five main biogeographic groups: consequences of Pleistocene climate changes or recent anthropogenic fragmentation?

Extant populations of the European wildcat are fragmented across the continent, the likely consequence of recent extirpations due to habitat loss and over‐hunting. However, their underlying phylogeographic history has never been reconstructed. For testing the hypothesis that the European wildcat survived the Ice Age fragmented in Mediterranean refuges, we assayed the genetic variation at 31 microsatellites in 668 presumptive European wildcats sampled in 15 European countries. Moreover, to evaluate the extent of subspecies/population divergence and identify eventual wild × domestic cat hybrids, we genotyped 26 African wildcats from Sardinia and North Africa and 294 random‐bred domestic cats. Results of multivariate analyses and Bayesian clustering confirmed that the European wild and the domestic cats (plus the African wildcats) belong to two well‐differentiated clusters (average Ф_ST = 0.159, R_st  = 0.392, P > 0.001; Analysis of molecular variance [AMOVA]). We identified from c. 5% to 10% cryptic hybrids in southern and central European populations. In contrast, wild‐living cats in Hungary and Scotland showed deep signatures of genetic admixture and introgression with domestic cats. The European wildcats are subdivided into five main genetic clusters (average Ф_ST = 0.103, R_st  = 0.143, P > 0.001; AMOVA) corresponding to five biogeographic groups, respectively, distributed in the Iberian Peninsula, central Europe, central Germany, Italian Peninsula and the island of Sicily, and in north‐eastern Italy and northern Balkan regions (Dinaric Alps). Approximate Bayesian Computation simulations supported late Pleistocene–early Holocene population splittings (from c. 60 k to 10 k years ago), contemporary to the last Ice Age climatic changes. These results provide evidences for wildcat Mediterranean refuges in southwestern Europe, but the evolution history of eastern wildcat populations remains to be clarified. Historical genetic subdivisions suggest conservation strategies aimed at enhancing gene flow through the restoration of ecological corridors within each biogeographic units. Concomitantly, the risk of hybridization with free‐ranging domestic cats along corridor edges should be carefully monitored.     

How reliable are morphological and anatomical characters to distinguish European wildcats,domestic cats and their hybrids in France?

Phenotypic variation in hybridizing species or subspecies is a prerequisite for allowing conservation ecologists and wildlife managers to identify parental populations and their hybrids in the field. We assessed the reliability of a set of eight morphological (body size and pelage characters) and four anatomical criteria (skull and intestine morphometric measurements) to distinguish between 302 French specimens classified as wildcat, domestic cat or hybrid on the basis of a Bayesian analysis (STRUCTURE) of their multilocus microsatellite genotypes. This aim was achieved by performing a set of multivariate analyses on morphological, anatomical and genetic data sets (Hill and Smith's analysis, co‐inertia analysis and discriminant analysis of principal components). Wildcats and domestic cats were very satisfactorily distinguished, even when using simple non‐invasive morphological criteria easily usable in the field like the morphology of the tail, dorsal line or flank stripes. Using anatomical instead of morphological characters slightly increased the discriminating power. Many more difficulties arose when we tried to distinguish hybrid specimens from both wildcat and domestic ones. Anatomical characters performed better than morphological ones in recognizing hybrids, but the assignment success rate remained very low, about 31.6% and 1.5%, respectively. Overall, the most discriminating characters were two continuous, derived anatomical characters: the cranial index followed by the intestinal index. Classification of specimens in three classes based on their microsatellite genotypes appeared to be inadequate for identifying hybrid specimens, as hybrid specimens seemed to be distributed along an anatomical continuum. With this observation in mind, we assessed the linear relationships between a proxy of the individual level of hybridization (q_ik) and the cranial and intestinal indices, respectively. Both relationships were highly significant. The greatest correlation was found with the cranial index (R² = 60.4%). Altogether, our results suggest that future work should be geared towards enhancing the measure of hybridization using more discriminating molecular markers and improving morphometric skull measurements through the use of modern geometric morphometric methods, using landmarks rather than skull dimension.