Abstract Anisakis pegreffii and A. simplex (s.s.) are the two zoonotic anisakids infecting cetaceans as well as pelagic/demersal fish and squids. In European waters, A. pegreffii prevails in the Mediterranean Sea, while A. simplex (s.s.) in the NE Atlantic Ocean. Abiotic conditions likely play a significant role in shaping their geographical distribution. The Iberian Atlantic and Alboran Sea waters are sympatric areas of the two species. A total of 429 adults and L3 stage from both sympatric and allopatric areas were studied by a wide nuclear genotyping approach (including newly and previously found diagnostic single nucleotide polymorphisms (SNPs) at nuclear DNA (nDNA) and microsatellite DNA loci) and sequenced at mitochondrial DNA (mtDNA) cox 2. Admixture between the two species was detected in the sympatric areas studied by STRUCTURE Bayesian analysis; NEWHYBRIDS revealed different categories of hybridization between the two species, representing approximately 5%. A tendency for F1 female hybrids to interbreed with the parental species at the geographical distribution limits of both species was observed. This finding suggests that hybridization occurs when the two parental species significantly differ in abundance. Mitochondrial introgression of A. simplex (s.s.) in A. pegreffii from Mediterranean waters was also detected, likely as a result of past and/or paleo-introgression events. The high level of genetic differentiation between the two species and their backcrosses indicates that, despite current hybridization, reproductive isolation which maintains evolutionary boundaries between the two species, exists. Possible causes of hybridization phenomena are attempted, as well as their evolutionary and ecological implications, also considering a sea warming scenario in European waters.