Inference of Raja miraletus population structure using mitochondrial and nuclear DNA: comparing the resolution power of molecular markers in exploring species and population boundaries

This study focused on the role of oceanographic discontinuities and the presence of transitional areas in shaping the population structure and the phylogeography of the Raja miraletus species complex, coupled with the test of the effective occurrence of past speciation events. The comparisons between the Atlantic African and the North-Eastern Atlantic-Mediterranean geographic populations were unravelled using both Cytochrome Oxidase I and eight microsatellite loci. This approach guaranteed a robust dataset for the identification of a speciation event between the Atlantic African clade, corresponding to the ex Raja ocellifera nominal species, and the NE Atlantic-Mediterranean R. miraletus clade. As a matter of fact, the origin of the Atlantic Africa and the NE Atlantic-Mediterranean deep split dated about 11.74MYA and was likely due to the synergic influence currents and two upwelling areas crossing the Western African Waters. Within the Mediterranean Sea, particular attention was also paid to the transitional area represented by Adventura and Maltese Bank, that might have contributed in sustaining the connectivity of the Western and the Eastern Mediterranean geographical populations. Furthermore, the geology of the easternmost part of Sicily and the geo-morphological depression of the Calabrian Arc could have driven the differentiation of the Eastern Mediterranean Sea. Although bathymetric and oceanographic discontinuity could represent barriers to dispersal and migration between Eastern and Western Mediterranean samples, a clear and complete genetic separation among them was not detected. Results produced by this work identified a speciation event defining Raja ocellifera and R. miraletus as two different species, and describing the R. miraletus species complex as the most ancient cryptic speciation event in the family Rajidae, representing another example of how strictly connected the environment, the behavioural habits and the evolutionary and ecologic drivers are.

DNA barcoding taxonomy and diversity of the Mediterranean sharks and rays

Elasmobranchs are an important by-catch of commercial fisheries targeting bony fishes. Fisheries targeting sharks are rare, but usually almost all specimen bycatched are marketed. They risk extinction if current fishing pressure continues (Ferretti et al., 2008). Accurate species identification is critical for the design of sustainable fisheries and appropriate management plans, especially since not all species are equally sensitive to fishing pressure (Walker & Hislop 1998). The identification of species constitutes the first basic step for biodiversity monitoring and conservation (Dayrat B et al., 2005). More recently, mtDNA sequencing has also been used for species identification and its use has become widespread under the DNA Barcode initiative (e.g. Hebert et al. 2003a, 2003b; Ward et al. 2005, 2008a; Moura et al 2008; Steinke et al. 2009). The aims of this work were: 1) identify sharks and skates species using DNA barcode; 2) compare species of different provenance; 3) use DNA barcode for misidentified species. Using DNA barcode 15 species of sharks (Alopias vulpinus, Centrophorus granulosus, Cetorhinus maximus, Dalatias licha, Etmopterus spinax, Galeorhinus galeus, Galeus melastomus, Heptranchias perlo, Hexanchus griseus, Mustelus mustelus, Mustelus punctulatus, Oxynotus centrina, Scyliorhinus canicula Squalus acanthias, Squalus blainville), 1 species of chimaera (Chimaera monstrosa) and 21 species of rays/skayes (Dasyatis centroura, Dasyatis pastinaca, Dasyatis sp., Dipturus nidarosiensis, Dipturus oxyrinchus, Leucoraja circularis, Leucoraja melitensis, Myliobatis aquila, Pteromylaeus bovinus, Pteroplatytrygon violacea, Raja asterias, Raja brachyura, Raja clavata, Raja miraletus, Raja montagui, Raja radula, Raja polystigma, Raja undulata, Rostroraja alba, Torpedo marmorata, Torpedo nobiliana, Torpedo torpedo) was identified.

Evolution and historical biogeography of the Eastern Atlantic skates

In my thesis, I tested the hypothesis that the diversification of the Eastern Atlantic skate faunas arose through vicariance rather than dispersal, using combined approach of molecular phylogeny reconstruction and zoogeography (namely historical biogeography). This analyses have been carried out independently on four Rajidae genera belonging to two different tribes: Rajini (Raja and Dipturus) and Amblyrajini (Rajella and Leucoraja). These taxa were selected because they displayed high species diversity and richness of endemic species in the Eastern Atlantic and Mediterranean. The verification of this hypothesis was carried out by reconstructing the best phylogenetic relationships among four genera and 26 species (including several endemism) based on mtDNA and nuDNA gene variation and several statistical approaches. Divergence times of taxa have been estimated based on molecular clock and fossil calibration to explain evolutionary patterns in the context of geological framework. Main issues are (i) the evidence that Eastern Atlantic skate evolution and displacement of species diversity occurred from pulsed geographical speciation (i.e. repeated series of parallel and independent speciation events) started in the Late Eocene-Early Miocene and they have occurred prevalently during Miocene; (ii) such relatively ancient origin of diversification has been allowed the sympatric displacement and evolution of several congeneric taxa likely because they have accumulated huge differences in the genomic and physiological/behavioural phenotypic traits; (iii) recently diverged sister species and taxa showed allopatric or parapatric evolution by the presence of oceanographic or hydrogeographical barriers which likely prevent large mixing between parapatric sister species.

Molecular phylogeny and historical zoogeography of four genera of Eastern Atlantic skates (Rajiformes, Rajidae)

The aim of this study was to reconstruct a solid phylogeny of four genera of the Rajidae family (Chondrichthyans: Batoidea) using a concatenated alignment of mtDNA genes. Then use the resultant tree to estimate divergence time between taxa based on molecular clock and fossil calibration and conduct biogeographic analysis. The intent was to prove that the actual distribution of species of Eastern Atlantic and Mediterranean skates is due to a series of vicariant events. The species considered belongs to two different tribe: Rajini (Raja and Dipturus) and Amblyrajini (Leucoraja and Rajella). The choice of this genera is due to their high presence in the area of interest and to the richness of endemic species. The results show that despite the ancient origin of Rajidae (97 MYA), the Eastern Atlantic and Mediterranean faunas originated more recently, during Middle Miocene-Late Pliocene, after the closure of connection between these areas and the Indo-Pacific ocean (15 MYA). The endemic species of the Mediterranean (Raja asterias, R. radula, R. polystigma and Leucoraja melitensis) originated after the Messinian salinity crisis (7-5 MYA), when the recolonization of the basin occurred, and are still maintained in allopatric distribution by the presence of biogeographic barriers. Moreover from 4 to 2.6 MYA we can observe the formation of sister species for Raja, Leucoraja and Rajella, one of which has a Northern distribution, and the other has a Southern distribution (R. clavata vs R. straeleni, L. wallacei vs L. naevus, R. fyllae vs R. caudaspinosa and R. kukujevi vs R. leopardus + R. barnardi). The Quaternary and present oceanographic discontinuities that occur along the western African continental shelf (e.g., Cape Blanc and the Angola–Benguela Front) might contribute to the maintenance of low or null levels of gene flow between these closely related siblings species. Also sympatric speciation must be invoked to explain the evolution of skates, for example for the division between R. leopardus and R. barnardi. The speciation processes followed a south-to-north pathways for Dipturus and a north-to-south pathways for Raja, Leucoraja and Rajella underling that the evolution of the genera occurred independently. In the end, it is conceivable that the evolutionary pathways of the tribes followed the costal line during the gondwana fragmentation. The results demonstrate that the evolution of this family is characterized by a series of parallel and independent speciation events, strictly correlated to the tectonic movement of continental masses and paleogeographic and paleoclimatic events and so can be explained by a panbiogeographical (vicariance) model.

Was there a second adaptive radiation of giant tortoises in the Indian Ocean? Using mitochondrial DNA to investigate speciation and biogeography of Aldabrachelys (Reptilia, Testudinidae)

A radiation of five species of giant tortoises (Cylindraspis ) existed in the southwest Indian Ocean, on the Mascarene islands, and another (of Aldabrachelys ) has been postulated on small islands north of Madagascar, from where at least eight nominal species have been named and up to five have been recently recognized. Of 37 specimens of Madagascan and small-island Aldabrachelys investigated by us, 23 yielded significant portions of a 428-base-pair (bp) fragment of mitochondrial (cytochrome b and tRNA-Glu), including type material of seven nominal species (A. arnoldi, A. dussumieri, A. hololissa, A. daudinii, A. sumierei, A. ponderosa and A. gouffei ). These and nearly all the remaining specimens, including 15 additional captive individuals sequenced previously, show little variation. Thirty-three exhibit no differences and the remainder diverge by only 1-4 bp (0.23-0.93%). This contrasts with more widely accepted tortoise species which show much greater inter- and intraspecific differences. The non-Madagascan material examined may therefore only represent a single species and all specimens may come from Aldabra where the common haplotype is known to occur. The present study provides no evidence against the Madagascan origin for Aldabra tortoises suggested by a previous molecular phylogenetic analysis, the direction of marine currents and phylogeography of other reptiles in the area. Ancient mitochondrial DNA from the extinct subfossil A. grandidieri of Madagascar differs at 25 sites (5.8%) from all other Aldabrachelys samples examined here.

Speciation and phylogeography in Caridina indistincta, a complex of freshwater shrimps from Australian heathland streams

The Caridina indistincta complex is a group of closely related atyid shrimps that inhabit coastal freshwater streams throughout north-eastern Australia. Using mitochondrial DNA sequence data (cytochrome oxidase 1, CO1), we (1) inferred the timing of speciation in the C. indistincta group and (2) examined the intraspecific phylogeographic patterns within the group. Assuming a shrimp-specific rate of CO1 evolution, the level of sequence divergence among species suggests that speciation took place during the Miocene epoch. Within one widespread mainland species, phylogeographic patterns suggest strong geographic 'regionalisation' of mtDNA lineages that are most likely of Pleistocene origin. By contrast, another species comprises two highly divergent mtDNA lineages that occur in sympatry. We suggest that although Pleistocene sea-level regressions appear important in generating population-level phylogeographic patterns, these events were largely unimportant in the formation of species in this group.

Inferring the degree of incipient speciation in secondary contact zones of closely related lineages of Palearctic green toads (Bufo viridis subgroup).

Reproductive isolation between lineages is expected to accumulate with divergence time, but the time taken to speciate may strongly vary between different groups of organisms. In anuran amphibians, laboratory crosses can still produce viable hybrid offspring >20 My after separation, but the speed of speciation in closely related anuran lineages under natural conditions is poorly studied. Palearctic green toads (Bufo viridis subgroup) offer an excellent system to address this question, comprising several lineages that arose at different times and form secondary contact zones. Using mitochondrial and nuclear markers, we previously demonstrated that in Sicily, B. siculus and B. balearicus developed advanced reproductive isolation after Plio-Pleistocene divergence (2.6 My, 3.3-1.9), with limited historic mtDNA introgression, scarce nuclear admixture, but low, if any, current gene flow. Here, we study genetic interactions between younger lineages of early Pleistocene divergence (1.9 My, 2.5-1.3) in northeastern Italy (B. balearicus, B. viridis). We find significantly more, asymmetric nuclear and wider, differential mtDNA introgression. The population structure seems to be molded by geographic distance and barriers (rivers), much more than by intrinsic genomic incompatibilities. These differences of hybridization between zones may be partly explained by differences in the duration of previous isolation. Scattered research on other anurans suggests that wide hybrid zones with strong introgression may develop when secondary contacts occur <2 My after divergence, whereas narrower zones with restricted gene flow form when divergence exceeds 3 My. Our study strengthens support for this rule of thumb by comparing lineages with different divergence times within the same radiation.

Rapid radiation in spiny lobsters (Palinurus spp) as revealed by classic and ABC methods using mtDNA and microsatellite data

Background: Molecular tools may help to uncover closely related and still diverging species from a wide variety of taxa and provide insight into the mechanisms, pace and geography of marine speciation. There is a certain controversy on the phylogeography and speciation modes of species-groups with an Eastern Atlantic-Western Indian Ocean distribution, with previous studies suggesting that older events (Miocene) and/or more recent (Pleistocene) oceanographic processes could have influenced the phylogeny of marine taxa. The spiny lobster genus Palinurus allows for testing among speciation hypotheses, since it has a particular distribution with two groups of three species each in the Northeastern Atlantic (P. elephas, P. mauritanicus and P. charlestoni) and Southeastern Atlantic and Southwestern Indian Oceans (P. gilchristi, P. delagoae and P. barbarae). In the present study, we obtain a more complete understanding of the phylogenetic relationships among these species through a combined dataset with both nuclear and mitochondrial markers, by testing alternative hypotheses on both the mutation rate and tree topology under the recently developed approximate Bayesian computation (ABC) methods. Results Our analyses support a North-to-South speciation pattern in Palinurus with all the South-African species forming a monophyletic clade nested within the Northern Hemisphere species. Coalescent-based ABC methods allowed us to reject the previously proposed hypothesis of a Middle Miocene speciation event related with the closure of the Tethyan Seaway. Instead, divergence times obtained for Palinurus species using the combined mtDNA-microsatellite dataset and standard mutation rates for mtDNA agree with known glaciation-related processes occurring during the last 2 my. Conclusion The Palinurus speciation pattern is a typical example of a series of rapid speciation events occurring within a group, with very short branches separating different species. Our results support the hypothesis that recent climate change-related oceanographic processes have influenced the phylogeny of marine taxa, with most Palinurus species originating during the last two million years. The present study highlights the value of new coalescent-based statistical methods such as ABC for testing different speciation hypotheses using molecular data.

Sex-linked genetic markers in the shrews of the Sorex araneus group : from free gene flow to speciation

SUMMARY: The shrews of the Sorex araneus group are morphologically .very similar, but have undergone a spectacular chromosomal evolution. Altogether, the shrews of this group present a complete array of every possible level of chromosomal and genetic differentiation. In South-Western Europe, four species are recognised: S. antiriorii, S. araneus, S. coronatus and S. granarius, which differ essentially by the amount and the composition of Robertsonian metacentric chromosomés. Additionally, several chromosome races of S. araneus are also present in the same region (i.e. Bretolet, Carlit, Cordon, Jura and Vaud). The objective of this thesis was to examine the genetic relationships between populations, races and /or species of the Sorex araneus group with a special emphasis onsex-specific markers (mtDNA and Y chromosome). We first investigate the evolutionary history of the shrews of the Sorex araneus group distributed in the South-Western Europe. The results of. these analyses confirmed the difficulty to draw a single dichotomic tree within this group. Incongruent mtDNA and Y chromosome phylogenies suggest further that genetic and chromosomal evolution are in this group partially independent processes and that the evolutionary history of the south-western European populations of the S. araneus group can only be understood if we consider secondary contacts between taxa, after their divergence (with genetic exchanges by means of hybridization and / or introgression). Using one male-inherited, one female inherited and eight biparentally inherited markers, we investigate the population genetic structure of the Valais shrew (Sorex antinorii). Overall there results suggest that two already well-differentiated genetic lineages colonized the Swiss Alps after the last glacial period and came into contact in the Rhône Valley. After the Valais shrew (Sorex antinorii) reached the Swiss Alps, it came into contact with the common shrew (Sorex araneus). When two species come into contact and hybridize, endogenous counter-selection of hybrids is usually first expressed as a reduced fertility or viability in hybrids of the heterogametic sex, a mechanism know as Haldane's rule (Haldane 1922). We first evaluated the extent of introgression for Y chromosome, mtDNA and autosomal markers in a hybrid zone between S. antinoriii and S. araneus. The overall level of genetic and karyotypic differentiation between the two species must be strong .enough to allow the detection asymmetric introgression. Secondly, we compared the levels of gene flow between chromosome common to both species and chromosome differently rearranged in each of them. We detected a significantly stronger genetic structure in rearranged chromosomes. Over a 10-year period, we even observed a decrease of genetic structure for common chromosomes. These results strongly support the role of chromosomal rearrangements in the reproductive barrier between S. araneus and S. anfinorii. Overall, this thesis underlines the need to use different inherited (paternally, maternally and / or biparentally) and chromosomally located (on common vs. on rearranged chromosomes) markers to obtain more accurate pictures of genetic relationships between populations or species. RÉSUMÉ: Les musaraignes du groupe Sorex araneus sont morphologiquement très proches, mais ont connu une spectaculaire évolution chromosomique. Prises dans leur ensemble, les musaraignes de ce groupe présentent tous les nivaux possibles de différenciation génétique et chromosomique. Dans le sud-ouest de l'Europe, quatre espèces appartenant à ce groupe sont présentes : S. antinorii, S. araneus, S. coronatus et S. granarius. Celles-ci diffèrent essentiellement par leur caryotype dont la variabilité est principalement due à des fusions Robertsoniennes. De plus, plusieurs races chromosomiques appartenant à S. araneus sont aussi présentes dans la même région (i.e. les races Bretolet, Carlit, Cordon, Jura et Vaud). L'objectif de cette thèse était d'examiner les relations génétiques entre populations, races et/ou espèces du groupe S. araneus, en utilisant particulièrement des marqueurs liés aux sexes (ADN mitochondrial et Chromosome Y). Nous avons dans un premier temps retracé l'histoire évolutive des musaraignes de ce groupe dans le sud-ouest de l'Europe. Les résultats dé ces analyses confirment qu'il est difficile de tracer un simple arbre dichotomique au sein de ce groupe. Les arbres phylogénétiques obtenus sur l'ADN mitochondrial et le chromosome Y sont incongruents et suggèrent de plus que l'évolution génétique et chromosomique sont des processus indépendants. L'histoire évolutive -des populations de ce groupe ne peut. être comprise qu'en considérant des contacts secondaires entre taxa postérieure à leur divergence et induisant des échanges génétiques par hybridation et/ou introgression. Par la suite, nous avons examiné la structure génétique des populations de la musaraigne du Valais, S. antinorii, en utilisant un marqueur transmis par les mâles, un marqueur transmis par les femelles et huit marqueurs transmis par les 2 sexes. Nos résultats suggèrent que deux lignées génétiquement bien différenciées aient colonisé les Alpes Suisses, après les dernières glaciations et entrent en contact dans là Vallée du Rhône. Après avoir franchi les Alpes Suisses, la musaraigne du Valais est entrée en contact avec là musaraigne commune (S. araneus). Lorsque deux espèces entrent en contact et s'hybrident, la sélection contre les hybrides implique habituellement une baisse de fertilité ou de viabilité des hybrides du sexe hétérogamétique (i.e. les mâles XY chez les mammifères). Ce mécanisme est connu sous le nom de règle de Haldane (Haldane 1922) et implique une plus forte structuration génétique de marqueurs males - spécifiques que des marqueurs femelles spécifiques. Nous avons donc évalué le degré d'introgression des marqueurs situés sur le chromosome Y, sur l'ADN mitochondrial et sur des autosomes dans une zone hybride entre S. araneus et S. antinorii. Le niveau de différenciation chromosomique et génétique entre les 2 espèces doit être suffisamment fort pour ne pas permettre la détection d'une introgression asymétrique entre les sexes. Dans un second temps, nous avons comparé les niveaux de flux de gênes mesurés à l'échelle du chromosome, pour des chromosomes communs aux deux espèces et pour des chromosomes différemment arrangées dans chacune des deux espèces. Nous avons détecté une structure génétique significativement plus forte sur les chromosomes réarrangés et comme la zone hybride a été étudiée à dix années d'intervalle, nous observons même une diminution de la structure génétique pour les chromosomes communs au cours du temps.. Ces résultats soutiennent fortement l'hypothèse d'un rôle des réarrangements chromosomiques dans l'établissement d'une barrière reproductive entre S. araneus et S. antinorii. Ainsi cette thèse souligne l'utilité d'utiliser des marqueurs génétiques avec différents modes de transmission. (par les mâles, par les femelles et/ou par les 2 sexes) ou localisés au niveau du chromosome (chromosomes communs vs chromosomes réarrangés) afin d'obtenir une image plus juste ou du moins plus complète des relations génétiques entre populations ou espèces.

Rapid radiation in spiny lobsters (Palinurus spp) as revealed by classic and ABC methods using mtDNA and microsatellite data

Background: Molecular tools may help to uncover closely related and still diverging species from a wide variety of taxa and provide insight into the mechanisms, pace and geography of marine speciation. There is a certain controversy on the phylogeography and speciation modes of species-groups with an Eastern Atlantic-Western Indian Ocean distribution, with previous studies suggesting that older events (Miocene) and/or more recent (Pleistocene) oceanographic processes could have influenced the phylogeny of marine taxa. The spiny lobster genus Palinurus allows for testing among speciation hypotheses, since it has a particular distribution with two groups of three species each in the Northeastern Atlantic (P. elephas, P. mauritanicus and P. charlestoni) and Southeastern Atlantic and Southwestern Indian Oceans (P. gilchristi, P. delagoae and P. barbarae). In the present study, we obtain a more complete understanding of the phylogenetic relationships among these species through a combined dataset with both nuclear and mitochondrial markers, by testing alternative hypotheses on both the mutation rate and tree topology under the recently developed approximate Bayesian computation (ABC) methods. Results: Our analyses support a North-to-South speciation pattern in Palinurus with all the South-African species forming a monophyletic clade nested within the Northern Hemisphere species. Coalescent-based ABC methods allowed us to reject the previously proposed hypothesis of a Middle Miocene speciation event related with the closure of the Tethyan Seaway. Instead, divergence times obtained for Palinurus species using the combined mtDNA-microsatellite dataset and standard mutation rates for mtDNA agree with known glaciation-related processes occurring during the last 2 my. Conclusion: The Palinurus speciation pattern is a typical example of a series of rapid speciation events occurring within a group, with very short branches separating different species. Our results support the hypothesis that recent climate change-related oceanographic processes have influenced the phylogeny of marine taxa, with most Palinurus species originating during the last two million years. The present study highlights the value of new coalescent-based statistical methods such as ABC for testing different speciation hypotheses using molecular data.

Comparisons of host mitochondrial, nuclear and endosymbiont bacterial genes reveal cryptic fig wasp species and the effects of Wolbachia on host mtDNA evolution and diversity

Background Figs and fig-pollinating wasp species usually display a highly specific one-to-one association. However, more and more studies have revealed that the "one-to-one" rule has been broken. Co-pollinators have been reported, but we do not yet know how they evolve. They may evolve from insect speciation induced or facilitated by Wolbachia which can manipulate host reproduction and induce reproductive isolation. In addition, Wolbachia can affect host mitochondrial DNA evolution, because of the linkage between Wolbachia and associated mitochondrial haplotypes, and thus confound host phylogeny based on mtDNA. Previous research has shown that fig wasps have the highest incidence of Wolbachia infection in all insect taxa, and Wolbachia may have great influence on fig wasp biology. Therefore, we look forward to understanding the influence of Wolbachia on mitochondrial DNA evolution and speciation in fig wasps. Results We surveyed 76 pollinator wasp specimens from nine Ficus microcarpa trees each growing at a different location in Hainan and Fujian Provinces, China. We found that all wasps were morphologically identified as Eupristina verticillata, but diverged into three clades with 4.22-5.28% mtDNA divergence and 2.29-20.72% nuclear gene divergence. We also found very strong concordance between E. verticillata clades and Wolbachia infection status, and the predicted effects of Wolbachia on both mtDNA diversity and evolution by decreasing mitochondrial haplotypes. Conclusions Our study reveals that the pollinating wasp E. verticillata on F. microcarpa has diverged into three cryptic species, and Wolbachia may have a role in this divergence. The results also indicate that Wolbachia strains infecting E. verticillata have likely resulted in selective sweeps on host mitochondrial DNA.

Molecular characterization of fire ants, Solenopsis spp., from Brazil based on analysis of mtDNA gene cytochrome oxidase I

Species from the Solenopsis saevissima (Smith) (Hymenoptera: Formicidae) species group are native to South America and have a cosmopolitan distribution because they have been accidentally introduced in many countries around the world. In Brazil, they have a wide distribution, including urban areas. The present study was conducted to investigate the characterization of Solenopsis genus populations associated with urban/human interference sites in Brazil by analyzing the mitochondrial gene cytochrome oxidase I and estimating the degree of relatedness of these populations to make inferences about their phylogeny and also observe the patterns of mitochondrial haplotype (mitotype) distribution across their range. The results revealed complete geographical coherence and polyphyly for the Solenopsis invicta Buren and Solenopsis saevissima species groups, which confirms the diversity of the genera. It also suggests the possibility that reproductively-isolated populations occur, resulting in the evolutionary process of speciation. No predominant haplotype was found in the populations analyzed, but some were more prevalent.

Structure and dynamics of hybrid zones at different stages of speciation in the common vole (Microtus arvalis)

The genetic structure and dynamics of hybrid zones provide crucial information for understanding the processes and mechanisms of evolutionary divergence and speciation. In general, higher levels of evolutionary divergence between taxa are more likely to be associated with reproductive isolation and may result in suppressed or strongly restricted hybridization. In this study, we examined two secondary contact zones between three deep evolutionary lineages in the common vole (Microtus arvalis). Differences in divergence times between the lineages can shed light on different stages of reproductive isolation and thus provide information on the ongoing speciation process in M. arvalis. We examined more than 800 individuals for mitochondrial (mtDNA), Y-chromosome and autosomal markers and used assignment and cline analysis methods to characterize the extent and direction of gene flow in the contact zones. Introgression of both autosomal and mtDNA markers in a relatively broad area of admixture indicates selectively neutral hybridization between the least-divergent lineages (Central and Eastern) without evidence for partial reproductive isolation. In contrast, a very narrow area of hybridization, shifts in marker clines and the quasi-absence of Y-chromosome introgression support a moving hybrid zone and unidirectional selection against male hybrids between the lineages with older divergence (Central and Western). Data from a replicate transect further support non-neutral processes in this hybrid zone and also suggest a role for landscape history in the movement and shaping of geneflow profiles.