Genetic analysis of a contact zone between two lineages of the ocellated lizard (Lacerta lepida Daudin 1802) in south‐eastern Iberia reveal a steep and narrow hybrid zone

Measuring the diffusion of genes between diverging taxa through zones of secondary contact is an essential step to understand the extent and nature of the reproductive isolation that has been achieved. Previous studies have shown that the ocellated lizard (Lacerta lepida Daudin, 1802) has endured repeated range fragmentation associated with the climatic oscillations of the Plio‐Pleistocene that promoted diversification of many different evolutionary units within the species. However, the oldest divergence within the group is estimated to have occurred much earlier, during the Miocene, around 9 Ma and corresponds to the split between the subspecies Lacerta lepida nevadensis Buchholz (1963) and Lacerta lepida lepida Daudin (1802). Although these two evolutionary units have documented genetic and morphological differentiation, most probably accumulated during periods of allopatry, little is known about patterns of gene flow between them. In this study, we performed a population genetic analysis of a putative area of secondary contact between these two taxa, using mtDNA and microsatellite data. We assessed levels of gene flow across the contact zone to clarify to what extent gene flow may be occurring. Hybridization between the subspecies was observed by the presence of genetically introgressed individuals. However, the overall coincidence of mitochondrial and multilocus nuclear clines and generally steep clines support the idea that this contact zone is acting as a barrier to gene flow. Taken together, these results suggest that L. l. lepida and L. l. nevadensis are in independent evolutionary trajectories and should be considered as two different species.