The genomics of speciation: investigating the molecular correlates of X chromosome introgression across the hybrid zone between Mus domesticus and Mus musculus

Understanding the genetic details of reproductive isolation is a key goal in the study of speciation. Hybrid zones, geographical regions where two species meet and exchange genes, can provide insight into the genetic basis of reproductive isolation. This is especially true in species with mapped molecular markers because patterns of gene flow can be compared among different genomic regions. Even greater insight can be obtained in species with complete genome sequences because gene identity, gene number and other features of interest can be assessed for genomic regions with different patterns of introgression. Here, we review recent studies on the well-characterized hybrid zone between Mus domesticus and M. musculus , including a detailed survey of patterns of introgression for 13 markers on the X chromosome. We then compare levels of introgression for these 13 regions to a number of genomic attributes inferred from the complete sequence of the X chromosome, with two purposes. First, we identify candidate genes for reproductive isolation by finding genes that map to an X-linked region of reduced introgression and that are only expressed in the male germ line or that show high rates of protein evolution in comparison with rat. Second, we ask whether patterns of gene flow are correlated with recombination rate, gene density, base composition, CpG island density, mutation rate and the rate of protein evolution, as might be expected if many genes contribute to reproductive isolation. We identify seven candidate genes for reproductive isolation between M. domesticus and M. musculus , and our analyses reveal no general correlations between levels of introgression and other measured sequence characteristics. These results underline the utility of the house mouse as a model system for the study of speciation.  © 2005 The Linnean Society of London, Biological Journal of the Linnean Society , 2005, 84 , 523–534.     

Genetic analysis of the hybridization zone between two subspecies Mus musculus domesticus and Mus musculus musculus in Bulgaria

The hybrid zone between the two subspecies of mice Mus musculus domesticus and Mus musculus musculus, which has been studied extensively in Denmark, crosses Europe to the Black Sea through the Alps and the Balkans. Two hundred and seventy-nine animals were captured in 22 localities along a transect across the Balkans. The animals were characterized for seven diagnostic nuclear loci by protein electrophoresis and by restriction pattern analysis of their mitochondrial DNA. The nuclear data show a sharp transition between the two subspecies, most of the variations in allele frequencies (from 0.9 to 0.1) occurring within a 36-km section of the transect. The introgression varies from one locus to the other and is more pronounced, in terms of distance, in M. m. musculus territory. Mitochondrial DNA introgression is important but occurs in one direction only, i.e. from M. m. musculus to M. m. domesticus, while a cytoplasmic transfer from M. m. domesticus to M. m. musculus has been reported. A previous study showed that no Y chromosome introgression occurs. The different behaviour of these three types of markers could be due to the interaction between selection against hybrid genomes and meiotic recombination. Objectively, it would appear that the genes that can introgress are neutral or nearly so and have been separated from deleterious genes they were linked to by recombination. This could explain the differential introgression between autosomal loci. The mitochondrial and Y chromosomes undergo no or very little recombination and each is transmitted as a whole. Their degree of introgression is thus indicative of the intensity of selection resulting from the amount of functional differentiation between the two taxa, which seems to be strong for the Y chromosome and weak for mitochondrial DNA. We propose that the asymmetry of nuclear introgression is due to different population structures. As M. m. musculus is relatively less structured, the rapid spreading of introgressed genes would be favoured. Such a scheme, however, can hardly account for the unidirectionality of the mitochondrial flow, which could be due to sex-dependent behaviour.     

Reproductive isolation between chromosomal races of the house mouse Mus musculus domesticus in a parapatric contact area revealed by an analysis of multiple unlinked loci

The house mouse, Mus musculus domesticus, exhibits a high level of chromosomal polymorphism because of the occurrence and fast fixation of Robertsonian fusions between telocentric chromosomes. For this reason, it has been considered a classical speciation model to analyse the role of the chromosomal changes in reproductive isolation. In this study, we analysed a parapatric contact area between two metacentric races in central Italy, the Cittaducale race (CD: 2n = 22) and the Ancarano race (ACR: 2n = 24), to estimate gene flow at the boundary. Hybrids between these two races show high levels of structural heterozygosity and are expected to be highly infertile. A sample of 88 mice from 14 sites was used. The mice were genotyped by means of eight microsatellite loci mapped in four different autosomal arms. The results show clear genetic differentiation between the CD and ACR races, as revealed by differences in allele frequencies, factorial correspondence analysis and indexes of genetic population (e.g. F(ST) and R(ST)) along the contact zone. The genetic differentiation between the races was further highlighted by assignation and clustering analyses, in which all the individuals were correctly assigned by their genotypes to the source chromosomal race. This result is particularly interesting in view of the absence of any geographical or ecological barrier in the parapatric contact zone, which occurs within a village. In these conditions, the observed genetic separation suggests an absence of gene flow between the races. The CD-ACR contact area is a rare example of a final stage of speciation between chromosomal races of rodents because of their chromosomal incompatibility.     

Hybrid zone formation and contrasting outcomes of secondary contact over transects in common toads

Much progress in speciation research stems from documenting patterns of morphological and genetic variation in hybrid zones. Contrasting patterns of marker introgression in different sections of the contact can provide valuable insights on the relative importance of various evolutionary mechanisms maintaining species differences in the face of hybridization and gene flow and on hybrid zone temporal and spatial dynamics. We studied species interactions in the common toads Bufo bufo and B. spinosus in France and northwestern Italy using morphological and molecular data from the mitochondrial and nuclear genomes in an extensive survey, including two independent transects west and east of the Alps. At both, we found sharp, coincident and concordant nuclear genetic transitions. However, morphological clines were wider or absent and mtDNA introgression was asymmetric. We discuss alternative, nonexclusive hypotheses about evolutionary processes generating these patterns, including drift, selection, long‐distance dispersal and spatial shifts in hybrid zone location and structure. The distribution of intraspecific mtDNA lineages supports a scenario in which B. bufo held a local refugium during the last glacial maximum. Present‐day genetic profiles are best explained by an advance of B. spinosus from a nearby Iberian refugium, largely superseding the local B. bufo population, followed by an advance of B. bufo from the Balkans, with prongs north and south of the Alps, driving B. spinosus southwards. A pendulum moving hybrid zone, first northwards and then southwards, explains the wide areas of introgression at either side of the current position of the contact zones.     

Genetic exchange across a hybrid zone within the Iberian endemic golden-striped salamander, Chioglossa lusitanica

The study of hybrid zones resulting from Pleistocene vicariance is central in examining the potential of genetically diverged evolutionary units either to introgress and merge or to proceed with further isolation. The hybrid zone between two mitochondrial lineages of Chioglossa lusitanica is located near the Mondego River in Central Portugal. We used mitochondrial and nuclear diagnostic markers to conduct a formal statistical analysis of the Chioglossa hybrid zone in the context of tension zone theory. Key results are: (i) cline centres are not coincident for all markers, with average widths of ca. 2-15 km; (ii) heterozygote deficit was not observed across loci near the transect centre; (iii) associations of parental allele combinations ('linkage disequilibrium'R) were not detected either across loci or across the transect. These observations suggest that the Chioglossa hybrid zone is not a tension zone with strong selection against hybrids but instead one shaped mostly by neutral mixing. The patterns uncovered suggest a complex history of populations over a small scale that may be common in southern Pleistocene refugia.     

Mitochondrial DNA variation in house mice, Mus domesticus (Rutty)

This article summarizes knowledge of house mouse diversity based on restriction analysis of mtDNA from 202 individuals representing 83 localities in Europe, North Africa, the Middle East, and the Americas. It begins by describing the variation among 34 newly collected European mice and relating the 15 new types of mtDNA found in them to the 37 types known from previous work. None of the new types represent deep new branches in the tree. Moreover, the order of branching in the tree reinforces the view that mice north of the Alps have southern origins. Two possible time scales are under consideration for the process of colonization of north-western Europe. According to one hypothesis, which assumes that the mtDNA clock in mice ticks at the standard mammalian rate, commensal mice colonized Europe north of the Alps roughly 30 000–70 000 years ago, perhaps in concert with the spread of anatomically modern Homo sapiens . Another hypothesis, which requires accelerated mtDNA evolution in mice, would have the colonization take place 5000–6000 years ago, or even more recently, following the spread of agriculture. The study also shows that mtDNAs closely resembling those M. domesticus mtDNAs that introgressed into the Danish and Swedish populations of M. musculus are widespread in M. domesticus populations of the region near the hybrid zone between these two species in Germany.     

Non-metric morphological divergence in the western house mouse, Mus musculus domesticus, from the Barcelona chromosomal hybrid zone

The effect of hybridization on morphological variation was investigated in 120 western house mice, Mus musculus domesticus , from the hybrid zone between the Barcelona and standard chromosomal races. The incidence of 37 non-metric cranial traits was calculated for standard mice (2 n  = 40) and Barcelona-standard hybrids (2 n  = 27–39). Subsequent analyses were conducted on several karyological subgroups, established by grouping the animals according to either their diploid number or their degree of chromosomal heterozygosity. Results revealed no significant difference by sex, asymmetry, or geographical distance. Significant phenetic divergences were found between the karyotypes studied in relation to several variants. Differences were especially substantial between the standard race and hybrid mice, even with respect to those hybrids with karyotypes close to that of the standard race. Within the hybrids, the maximum divergence corresponded to the 28-chromosome homozygotes, chromosomally close to the Barcelona race, and to the heterozygotes with more than two fusions. Since differences in non-metric trait frequencies are generally considered a measure of genetic divergence, the results suggest the occurrence of a barrier to gene flow in the Barcelona hybrid zone. The decrease of genetic exchange between the chromosomally differentiated mice might be due to reduced fertility in hybrids, associated with chromosomal heterozygosity.  © 2003 The Linnean Society of London, Biological Journal of the Linnean Society , 2003, 80 , 313–322.     

Steep clines within a highly permeable genome across a hybrid zone between two subspecies of the European rabbit

Maintenance of genetic distinction in the face of gene flow is an important aspect of the speciation process. Here, we provide a detailed spatial and genetic characterization of a hybrid zone between two subspecies of the European rabbit. We examined patterns of allele frequency change for 22 markers located on the autosomes, X‐chromosome, Y‐chromosome and mtDNA in 1078 individuals sampled across the hybrid zone. While some loci revealed extremely wide clines (w ≥ 300 km) relative to an estimated dispersal of 1.95–4.22 km/generation, others showed abrupt transitions (w ≈ 10 km), indicating localized genomic regions of strong selection against introgression. The subset of loci showing steep clines had largely coincident centers and stepped changes in allele frequency that did not co‐localize with any physical barrier or ecotone, suggesting that the rabbit hybrid zone is a tension zone. The steepest clines were for X‐ and Y‐chromosome markers. Our results are consistent with previous inference based on DNA sequence variation of individuals sampled in allopatry in suggesting that a large proportion of each genome has escaped the overall barrier to gene flow in the middle of the hybrid zone. These results imply an old history of hybridization and high effective gene flow and anticipate that isolation factors should often localize to small genomic regions.     

Geographical pattern of genetic variation in the Robertsonian system of Mus musculus domesticus in central Italy

The Robertsonian (Rb) system of Mus musculus domesticus in central Italy consists of three Rb races (CD, 2 n  = 22; CB, 2 n  = 22; ACR, 2 n  = 24), and additional Rb populations with new metacentrics have recently been reported. The aim of the present paper is to provide insight into the process of house mouse raciation in this area. Here we present new data concerning the geographical distribution of mtDNA haplotypes (control region sequence) in the three Rb races and in the surrounding standard populations. The presence of exclusive clades in all the three Rb races supports the idea of a past fragmentation event and a subsequent period of isolation. However, comparison of mtDNA polymorphism between Rb and standard races, and the shape of the mismatch distribution within the races, does not reveal any clear sign of a strong reduction in population size, at this geographical scale, related to race formation.  © 2005 The Linnean Society of London, Biological Journal of the Linnean Society , 2005, 84 , 395–405.     

Discordant patterns of introgression across a narrow hybrid zone between two cryptic lineages of an Iberian endemic newt

The study of natural hybrid zones can illuminate aspects of lineage divergence and speciation in morphologically cryptic taxa. We studied a hybrid zone between two highly divergent but morphologically similar lineages (south‐western and south‐eastern) of the Iberian endemic Bosca's newt (Lissotriton boscai) in SW Iberia with a multilocus dataset (microsatellites, nuclear and mitochondrial genes). STRUCTURE and NEWHYBRIDS analyses retrieved few admixed individuals, which classified as backcrosses involving parental individuals of the south‐western lineage. Our results show asymmetric introgression of mtDNA beyond the contact from this lineage into the south‐eastern lineage. Analysis of nongeographic introgression patterns revealed asymmetries in the direction of introgression, but except for mtDNA, we did not find evidence for nonconcordant introgression patterns across nuclear loci. Analysis of a 150‐km transect across the hybrid zone showed broadly coincident cline widths (ca. 3.2–27.9 km), and concordant cline centres across all markers, except for mtDNA that is displaced ca. 60 km northward. Results from ecological niche modelling show that the hybrid zone is in a climatically homogenous area with suitable habitat for the species, suggesting that contact between the two lineages is unlikely to occur further south as their distributions are currently separated by an extensive area of unfavourable habitat. Taken together, our findings suggest the genetic structure of this hybrid zone results from the interplay of historical (biogeographic) and population‐level processes. The narrowness and coincidence of genetic clines can be explained by weak selection against hybrids and reflect a degree of reproductive isolation that is consistent with cryptic speciation.     

The role of salivary androgen-binding protein in reproductive isolation between two subspecies of house mouse: Mus musculus musculus and Mus musculus domesticus

Previous behavioural studies using inbred lines have suggested that the gene ( Abpa ) for the alpha subunit of salivary androgen-binding protein (ABP) plays a role in prezygotic isolation between house mouse Mus musculus subspecies. We tested this hypothesis in animals from wild allopatric (121 individuals from four samples) and parapatric (320 animals from 15 samples) populations sampled on the Czech–Bavarian transect across the hybrid zone between M. m. domesticus and M. m. musculus . The study did not reveal a consistent statistically significant trend of homosubspecific preferences in individual allopatric and parapatric populations. Nonetheless, the whole pattern of preference was skewed toward homosubspecific preference mostly on the M. m. musculus side of the hybrid zone. The pattern of homosubspecific preferences was stronger for the time spent sniffing than it was for the first choice of the signal (the ratio of homosubspecific vs. heterosubspecific preferences for both sexes was 6 : 2 in allopatric and 21 : 9 in parapatric populations, while the same rates were 4 : 4 and 16 : 14 for the first choice). To the extent that Y-maze tests reflect preference under wild conditions, we suggest that this slight preference may not in itself be sufficient to impede gene flow between the two subspecies and thus act as a reproductive barrier. ABP most probably participates in a complex system of subspecies-specific recognition in the hybrid zone, but the picture is far too complex at this time to allow a conclusive evaluation of the importance of this role.  © 2005 The Linnean Society of London, Biological Journal of the Linnean Society , 2005, 84 , 349–361.     

Genetic analysis of a chromosomal hybrid zone in the Australian morabine grasshoppers (Vandiemenella, viatica species group)

Whether chromosomal rearrangements promote speciation by providing barriers to gene exchange between populations is one of the long-standing debates in evolutionary biology. This question can be addressed by studying patterns of gene flow and selection in hybrid zones between chromosomally diverse taxa. Here we present results of the first study of the genetic structure of a hybrid zone between chromosomal races of morabine grasshoppers Vandiemenella viatica , P24(XY) and viatica 17, on Kangaroo Island, Australia. Chromosomal and 11 nuclear markers revealed a narrow hybrid zone with strong linkage disequilibrium and heterozygote deficits, most likely maintained by a balance between dispersal and selection. Widths and positions of clines for these markers are concordant and coincident, suggesting that selection is unlikely to be concentrated on a few chromosomes. In contrast, a mitochondrial marker showed a significantly wider cline with centre offset toward the P24(XY) side. We argue that the discordance between the mitochondrial and nuclear/chromosomal clines and overall asymmetry of the clines suggest a secondary origin of the contact zone and potential movement of the zone after contact. Genome-wide scans using many genetic markers and chromosomal mapping of these markers are needed to investigate whether chromosomal differences directly reduce gene flow after secondary contact.     

Strong selection on male plumage in a hybrid zone between a hybrid bird species and one of its parents

Homoploid hybrid speciation (HHS) requires reproductive barriers between hybrid and parent species, despite incomplete reproductive isolation (RI) between the parents. Novel secondary sexual trait values in hybrids may cause prezygotic isolation from both parents, whereas signals inherited by the hybrid from one parent species may cause prezygotic isolation with the other. Here we investigate whether differences in male plumage function as a premating barrier between the hybrid Italian sparrow and one of its parent species, the house sparrow, in a narrow Alpine hybrid zone. Italian sparrow male plumage is a composite mosaic of the parental traits, with its head plumage most similar to its other parent, the Spanish sparrow. We use geographical cline analysis to examine selection on three plumage traits, 75 nuclear single nucleotide polymorphisms (SNPs) and hybrid indices based on these SNPs. Several SNPs showed evidence of restricted introgression in the Alps, supporting earlier findings. Crown colour exhibited the narrowest plumage cline, representing a 37% (range 4–65%) drop in fitness. The cline was too narrow to be due to neutral introgression. Only crown colour was significantly bimodal in the hybrid zone. Bimodality may be due to RI or a major QTL, although fitness estimates suggest that selection contributes to the pattern. We discuss the implications with respect to HHS and the species status of the Italian sparrow.     

Evolution of the chromosomal races of Mus musculus domesticus in the Rhaetian Alps: the roles of whole-arm reciprocal translocation and zonal raciation

The evolution of ten Robertsonian (Rb) races of the house mouse (Mus musculus domesticus) in the Rhaetian Alps of northern Italy and southern Switzerland is reconsidered. The mechanisms of centric fusion, zonal raciation and, for the first time, whole-arm reciprocal translocation (WART), are used in this non-mathematical approach to produce a phylogenetic tree (using chromosome fusions as characters) with the smallest number of steps. The shortest tree that we found (16 steps) is at least two to nine mutations shorter than previously published models. Three other trees (17 or 18 steps) are also considered, since they are geographically more sensible. In general, these four scenarios correspond more closely to the present distributions of the ten Rb races than previous trees. Our results suggest that zonal raciation and WARTs play an important role in the evolution of Rb races of the house mouse.     

Evolutionary modularity of the mouse mandible: dissecting the effect of chromosomal reorganizations and isolation by distance in a Robertsonian system of Mus musculus domesticus

The mouse mandible consists of several morphogenetic units that are usually grouped into two main modules: the alveolar region and the ascending ramus. The genetic/ontogenetic modularity of the two regions implies that they might evolve independently to some extent. In particular, evolutionary modularity in quantitative traits could arise during chromosomal speciation due to lower gene flow in rearranged chromosomes. With the aim of uncovering the autonomous evolution of the mandible modules, the form variation of each of them was assessed in the house mouse Robertsonian system from Barcelona, in which chromosomal variation and geographical distance may act as isolation factors. The association between these factors and morphological changes was analysed to determine their contribution to the differentiation of each module. Although size changes in the two modules were highly correlated, shape changes were not, and their association with karyotype differences, but not geographical distance, was dependent on the module. The results support the existence of two evolutionary modules and highlight the importance of size in morphological integration of the mandible. They also suggest that geographical distance and chromosomal reorganizations reduce gene flow between karyotypically divergent populations, but although geographical distance represents a global barrier to gene flow, the isolation produced by a set of chromosomal reorganizations only affects particular modules, probably depending on the number and location of loci with effects on a particular morphological region.     

Mitochondrial function and bioenergetic trade‐offs during lactation in the house mouse (Mus musculus)

Energy allocation theory predicts that a lactating female should alter the energetic demands of its organ systems in a manner that maximizes nutrient allocation to reproduction while reducing nutrient use for tasks that are not vital to immediate survival. We posit that organ‐specific plasticity in the function of mitochondria plays a key role in mediating these energetic trade‐offs. The goal of this project was to evaluate mitochondrial changes that occur in response to lactation in two of the most energetically demanding organs in the body of a rodent, the liver and skeletal muscle. This work was conducted in wild‐derived house mice (Mus musculus) kept in seminatural enclosures that allow the mice to maintain a natural social structure and move within a home range size typical of wild mice. Tissues were collected from females at peak lactation and from age‐matched nonreproductive females. Mitochondrial respiration, oxidative damage, antioxidant, PGC‐1α, and uncoupling protein levels were compared between lactating and nonreproductive females. Our findings suggest that both liver and skeletal muscle downregulate specific antioxidant proteins during lactation. The liver, but not skeletal muscle, of lactating females displayed higher oxidative damage than nonreproductive females. The liver mass of lactating females increased, but the liver displayed no change in mitochondrial respiratory control ratio. Skeletal muscle mass and mitochondrial respiratory control ratio were not different between groups. However, the respiratory function of skeletal muscle did vary among lactating females as a function of stage of concurrent pregnancy, litter size, and mass of the mammary glands. The observed changes are predicted to increase the efficiency of skeletal muscle mitochondria, reducing the substrate demands of skeletal muscle during lactation. Differences between our results and prior studies highlight the role that an animals’ social and physical environment could play in how it adapts to the energetic demands of reproduction.     

Studies on the Robertsonian chromosomal variation of Mus musculus domesticus (Rodentia, Muridae) in Greece

The karyotype of the house mouse, Mus musculus domesticus , was examined in 282 specimens from 44 localities, in an effort to gain better understanding of the Robertsonian (Rb) variation known to exist in Greece. We consider that an Rb system exists in Peloponnisos, southern Greece, distributed in an area that is substantially larger than previously known. It consists of at least three Rb races with 2 n  = 30, 2 n  = 24 and 2 n  = 28, respectively, the last being reported for the first time in this paper and carrying Rb(3.6), Rb(8.12), Rb(10.14), Rb(13.15), Rb(9.16) and Rb(11.17) in a homozygous state. Additional instances of variation in this Rb system include individuals with 2 n  = 31 and 32 of variable Rb constitution and hybrids between the Rb races with 2 n  = 30 and 2 n  = 24. In southern Peloponnisos, Rb(10.14) was found in either a homozygous or a heterozygous state (2 n  = 38 or 39). The relationships among the Rb populations of Peloponnisos are discussed and hypotheses for their evolution are proposed. Rb variation was also recorded in two new locations of eastern Sterea Ellas (2 n  = 28 and 29) and one in Ipiros, north-west Greece (2 n  = 38). These findings corroborate the existence of two separate Rb systems in those two areas. Moreover, among a number of islands surveyed, Rb variation was only found in Kythira island, with Rb(10.14) in a heterozygous state (2 n  = 39). Finally, the typical all-acrocentric karyotype (2 n  = 40) was found in 51 of the animals studied from 13 localities. © 2005 The Linnean Society of London, Biological Journal of the Linnean Society , 2005, 84 , 503–513.     

Assessment of gene flow across a hybrid zone in red-tailed chipmunks (Tamias ruficaudus)

The role of hybridization in animal speciation is controversial and recent research has challenged the long-standing criterion of complete reproductive isolation to attain species status. The speciation-with-gene-flow model posits that the genome is semi-permeable and hybridization may be a phase in the process of divergence. Here, we apply these concepts to a previously identified zone of mtDNA introgression between the two strongly morphologically differentiated subspecies of red-tailed chipmunk ( Tamias ruficaudus ) in the US Inland Northwest. Using multilocus genotype data from the southern, older contact zone, we demonstrate that neutral gene flow is unusually low between the subspecies across the Lochsa River. This is geographically congruent with the discontinuity in bacular morphology, indicating that the cline of mitochondrial DNA (mtDNA) haplotypes is displaced. Furthermore, we elucidate the evolutionary forces responsible by testing hypotheses of lineage sorting and hybridization. We determined that introgressive hybridization is the cause of mtDNA/morphology incongruence because there are non-zero levels of migration and gene flow. Although our estimate of the age of the hybrid zone has wide credibility intervals, the hybridization events occurred in the Late Pleistocene and the divergence occurred in the Middle Pleistocene. Finally, we assessed substructure within and adjacent to the hybrid zone and found that the hybrid zone constitutes a set of populations that are genetically differentiated from parental sets of populations; therefore, hybridization in this system is not likely an evolutionary sink, but has generated novel combinations of genotypes.     

Mitochondrial ghost lineages blur phylogeography and taxonomy of Natrix helvetica and N. natrix in Italy and Corsica

Grass snakes are widely distributed across the Western Palearctic. Recent phylogeographic studies provided evidence that three distinct parapatric species exist. Two of these occur in Italy, Natrix helvetica and N. natrix, and a contact zone between both taxa has been suggested for north-eastern Italy. Moreover, previous investigations revealed for the Italian Peninsula a complex phylogeographic structure. Using mtDNA sequences and microsatellite loci, we examined the situation for mainland Italy, Sicily, Sardinia, and Corsica. Our study confirmed the occurrence of N. natrix in north-eastern Italy. Cline analyses revealed limited gene flow between N. helvetica and N. natrix across a narrow hybrid zone. Within N. helvetica, conflicting patterns of mitochondrial and nuclear genomic differentiation were revealed. Three nuclear genomic clusters were found; one of them corresponded to no fewer than five distinct and, in part, deeply divergent and ancient mitochondrial lineages from mainland Italy and Sicily. This cluster was paraphyletic with respect to the two remaining mitochondrial lineages, each of which matched with another nuclear genomic cluster (one from Corsica plus Sardinia and another one from western Europe north of the Alps). This unexpected pattern most likely results from mainly male-mediated gene flow and female philopatry combined with population-density-dependent processes such as ‘high-density blocking’. With respect to taxonomy, we propose to synonymize N. h. lanzai Kramer, 1970 with N. h. sicula (Cuvier, 1829), acknowledging their lacking nuclear genomic differentiation. The studied hybrid zone of N. h. helvetica and N. h. sicula in Italy is wide, with a smooth cline for nuclear markers, supporting their subspecies status. We found no evidence for the distinctiveness of the two subspecies from Corsica (N. h. corsa) and Sardinia (N. h. cetti), suggesting their synonymy, but refrain from taxonomic conclusions because of small sample sizes and the endangered status of the Sardinian taxon.