Inferences of selection and migration in the Danish house mouse hybrid zone

We analysed the patterns of allele frequency change for ten diagnostic autosomal allozyme loci in the hybrid zone between the house mouse subspecies Mus musculus domesticus and M. m. musculus in central Jutland. After determining the general orientation of the clines of allele frequencies, we analysed the cline shapes along the direction of maximum gradient. Eight of the ten clines are best described by steep central steps with coincident positions and an average width of 8.9 km (support limits 7.6–12.4) flanked by tails of introgression, indicating the existence of a barrier to gene flow and only weak selection on the loci studied. We derived estimates of migration from linkage disequilibrium in the centre of the zone, and by applying isolation by distance methods to microsatellite data from some of these populations. These give concordant estimates of σ =  0.5–0.8 km generation     . The barrier to gene flow is of the order of 20 km (support limits 14–28), and could be explained by selection of a few per cent at 43–120 underdominant loci that reduces the mean fitness in the central populations to 0.45. Some of the clines appear symmetrical, whereas others are strongly asymmetrical, and two loci appear to have escaped the central barrier to gene flow, reflecting the differential action of selection on different parts of the genome. Asymmetry is always in the direction of more introgression into musculus , indicating either a general progression of domesticus into the musculus territory, possibly mediated by differential behaviour, or past movement of the hybrid zone in the opposite direction, impeded by potential geographical barriers to migration in domesticus territory.  © 2005 The Linnean Society of London, Biological Journal of the Linnean Society , 2005, 84 , 593–616.     

Inference of selection and stochastic effects in the house mouse hybrid zone

We explored the transition of 13 X-linked markers across two separate portions of the house mouse hybrid zone, asking whether such a comparison can distinguish the effects of selection from random factors. A heuristic search in the likelihood landscape revealed more complex likelihood profiles for data sampled in two-dimensional (2D) space relative to data sampled along a linear transect. Randomized resampling of localities analyzed for individual loci showed that deletion of sites away from the zone center can decrease cline width estimates whereas deletion of sites close to the center can significantly increase the width estimates. Deleting localities for all loci resulted in wider clines if the number of samples from the center was limited. The results suggest that, given the great variation in width estimates resulting from inclusion/exclusion of sampling sites, the geographic sampling design is important in hybrid zone studies and that our inferences should take into account measures of uncertainty such as support intervals. The comparison of the two transects indicates cline widths are narrower for loci in the central part of the X chromosome, suggesting selection is stronger in this region and genetic incompatibilities may have at least partly common architecture in the house mouse hybrid zone.     

Testing for selection on the androgen-binding protein in the Danish mouse hybrid zone

The three peripheral subspecies of the house mouse Mus musculus have fixed specific variants of the androgen-binding protein (ABP) that have been proposed to be part of a recognition mechanism that could participate in sexual isolation between the subspecies. We tested for selection on Abpa by characterizing the pattern of Abpa introgression across a transect of the hybrid zone between M. m. musculus and M. m. domesticus in Jutland. On the musculus side, the cline for Abpa resembled that of a nearly neutral allozyme more than that of strongly selected X and Y chromosome markers. However, the high central step which displaces the tail of introgression of Abpa to higher frequencies was best accounted for by linkage to a locus under strong selection against hybrids. Still, we cannot exclude that this pattern results from weak selection against Abpa in the tail of introgression, which would be compatible with an assortative choice mechanism. On the domesticus side there was little introgression close to the hybrid zone, presumably due to a geographical barrier to migration. However, substantial frequencies of musculus alleles occurred further away, suggesting mixed colonization patterns as well as flow across the hybrid zone.  © 2005 The Linnean Society of London, Biological Journal of the Linnean Society , 2005, 84 , 447–459.     

Counterselection on sex chromosomes in the Mus musculus European hybrid zone

The extent to which alleles can disperse across a hybrid zone depends on the selection they are subjected to in the hybrid genetic background or, for those that are selectively neutral, on their ability to escape from the unfavourable environment by recombination. Three markers spanning a 45 cM segment in the center of the X chromosome were used to investigate the degree to which selection against X chromosome linked genes helps to maintain the barrier to gene flow in the hybrid zone between Mus musculus domesticus and M. m. musculus in Denmark. The introgression of all the sex chromosome specific markers was more limited than that of the autosomal enzymes (Idh1, Amy, Gpd1, Pgm1, Es1, Es2, Mpi, Np1, Es10, Sod1) and the mitochondrial DNA. The cline for DXPas2, which is in the center of the X chromosome, is extremely steep and shows that certain genes located in this region are strongly selected against in the hybrid background. The clines of the other two X-linked markers, Hprt and DXPas1, and of the Y chromosome are not as abrupt and all three have similar asymmetric introgression patterns. Although the musculus variants appear to behave in much the same way as those of the autosomal genes, the domesticus variants do not introgress. The results show that X-linked and to a lesser extent Y-linked genes are more strongly selected against in the hybrid genome than the mitochondrial genome or the different autosomal loci. This suggests that co-adapted gene systems involving the sex chromosomes may play an important role in the hybrid breakdown between the two subspecies.     

Haldane's rule in an avian system: using cline theory and divergence population genetics to test for differential introgression of mitochondrial, autosomal, and sex-linked loci across the Passerina bunting hybrid zone

Using cline fitting and divergence population genetics, we tested a prediction of Haldane's rule: autosomal alleles should introgress more than z-linked alleles or mitochondrial haplotypes across the Passerina amoena/Passerina cyanea (Aves: Cardinalidae) hybrid zone. We screened 222 individuals collected along a transect in the Great Plains of North America that spans the contact zone for mitochondrial (two genes), autosomal (four loci) and z-linked (two loci) markers. Maximum-likelihood cline widths estimated from the mitochondrial (223 km) and z-linked (309 km) datasets were significantly narrower on average than the autosomal cline widths (466 km). We also found that mean coalescent-based estimates of introgression were larger for the autosomal loci (0.63 genes/generation, scaled to the mutation rate mu) than for both the mitochondrial (0.27) and z-linked loci (0.59). These patterns are consistent with Haldane's rule, but the among-locus variation also suggests many independently segregating loci are required to investigate introgression patterns across the genome. These results provide the first comprehensive comparison of mitochondrial, sex-linked, and autosomal loci across an avian hybrid zone and add to the body of evidence suggesting that sex chromosomes play an important role in the formation and maintenance of reproductive isolation between closely related species.     

Chromosomal introgression in house mice from the hybrid zone between M. m. domesticus and M. m. musculus in Denmark

The recent discovery of Robertsonian (Rb) translocations in Danish mice from the hybrid zone between Mus musculus musculus and M. m. domesticus stimulated the chromosomal analysis of populations along a north-south transect through this zone. G-Banding identified the Rb fusions as Rb(3.8), Rb(2.5) and Rb(6.9). The cytogenetic results show that there is a gradual decrease in the number of fusions as one proceeds north, the translocations abruptly ending in populations from the centre of the hybrid zone determined by seven diagnostic allozymic markers. These results indicate that Rb fusions are present only in domesticus or predominantly domesticus-genotype mice and that they do not introgress into M. m. musculus . To test if genie incompatibilities between the musculus genetic background and Rb fusions were involved in the systematic elimination of the latter, predominantly musculus mice from the hybrid zone were crossed with Rb domesticus mice carrying Rb(3.8). The karyotypic analysis of the progeny showed no distortion of the transmission ratio of this fusion.
The chromosomal and allozymic analysis of these mice further indicates that (i) recombination is not suppressed between metacentrics and their acrocentric homologues and (ii) specific domesticus chromosomal segments are tolerated in the musculus genomes whereas the Rb centromeres are not.     

Transgressive segregation in a behavioural trait? Explorative strategies in two house mouse subspecies and their hybrids

Hybrid zones between genetically diverged populations are widespread among animals and plants. Their dynamics usually depend on selection against admixture and dispersal of parental forms in the zone. Although indirect estimates of selection have been the target of many studies, dispersal has been neglected. In this study we carried out open field experiments to test whether males of two house mouse subspecies, Mus musculus musculus and Mus musculus domesticus, differ in their propensity to disperse and in their character of exploration. We tested wild‐caught males and males of two wild‐derived inbred strains. In addition, we examined reciprocal F₁ crosses to test the prediction that these hybrids display intermediate behaviours. We revealed that M. m. musculus males were less hesitant to enter the experimental arena than were M. m. domesticus males, but once inside the arena their movements were more timid. F₁ males differed from both parental strains, with longer latencies to enter the arena, but explored the arena in a similar fashion as the M. m. domesticus males, thus displaying transgressive behavioural phenotypes. These results contribute to our knowledge of behavioural divergence between the mouse subspecies, and add a new facet to the study of speciation. © 2012 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, ●●, ●●–●●.     

Speciation in Passerina buntings: introgression patterns of sex-linked loci identify a candidate gene region for reproductive isolation

Sex-chromosomes are thought to play an important role in speciation, but few studies of non-model organisms have investigated the relative influence of multiple sex-linked markers on reproductive isolation. We collected 222 individuals along a geographical transect spanning the hybrid zone between Passerina amoena and P. cyanea (Aves: Cardinalidae). Using maximum-likelihood cline fitting methods, we estimated locus-specific introgression rates for 10 z-linked markers. Although the cline width estimates ranged from 2.8 to 584 km, eight of 10 loci had cline widths between 224 and 271 km. We also used coalescent-based estimates of locus-specific divergence times between P. amoena and P. cyanea to test a recently proposed hypothesis of an inverse relationship between divergence time and cline width but did not find a significant association. The narrow width (2.8 km) of the cline estimated from the VLDLR9 locus indicates strong selection retarding introgression of alleles at this locus across the hybrid zone. Interestingly, a mutation in the very low density lipoprotein receptor ( VLDLR ) gene, in which VLDLR9 is an intron, is known to reduce the egg-laying ability of some chickens, suggesting a possible link between this gene region and reproductive isolation between P. amoena and P. cyanea . These results underscore the importance of sampling multiple loci to investigate introgression patterns across a chromosome or genome and support previous findings of the importance of sex-linked genes in speciation.     

ABRUPT CLINE FOR SEX CHROMOSOMES IN A HYBRID ZONE BETWEEN TWO SPECIES OF MICE

We compared the patterns of movement of sex chromosomal and autosomal loci along a 160 km transect across a zone of hybridization between M. domesticus and M. musculus in southern Germany and western Austria using seven genetic markers. These included one Y-specific DNA sequence (YB10), two X-specific loci (DXWas68 and DXWas31), and four autosomal isozyme loci (Es-10, Es-1, Mpi-1, and Np-1). Random effects logistic regression analysis enabled us to examine the relationship between M. domesticus allele frequency and geographic distance from the western edge of the hybrid zone and allowed statistical evaluation of differences in cline midpoint and width among loci. More limited movement was observed for all three sex chromosomal markers across the zone compared with three of the four autosomal markers. If differential movement reflects fitness differences of specific alleles (or alleles at closely linked loci) on a hybrid background, then alleles that move to a limited extent across a hybrid zone may contribute to hybrid breakdown between two species. The limited flow of both X- and Y-specific alleles suggest that sex chromosomes have played an important role in Mus speciation.     

THE PONTIA DAPLIDICE-ED USA HYBRID ZONE IN NORTHWESTERN ITALY

The pierid butterflies Pontia daplidice and P. edusa, parapatrically distributed in southern Europe, have very similar morphologies and life histories, but show fixed differences at four allozyme markers. We sampled these allozymes in a 28-population transect north of Genoa in Italy, through the hybrid zone where these taxa meet. We used the numerical techniques developed for hybrid zone analysis to study the patterns of genetic differentiation and their underlying evolutionary causes. The hybrid zone is characterized by a very short and steep central region, flanked by broad tails of introgression extended up to 100 km in either direction. From mean two-locus disequilibium of D = 0.148 (maximum-likelihood two-unit support limits 0.139-0.153), and after accounting for minor differences in the center locations of the single-locus clines, which act to bias the dispersal estimate, we estimated a dispersal rate of σ = 4.4 (3.7-5.5) km/gen^1/2. The effective selection needed to maintain the steep central portion is strong, 0.47 < s? < 0.64, when combined over potential intrinsic (genetic background) and extrinsic (ecological) sources of selection. The clines in allozyme loci showed variation that was significantly different between the most divergent shapes, and the differences are attributable to different degrees of introgression on the edusa side of the zone. The average selection acting on individual allozyme loci was high at s_???e  1.5%, but because of the narrowness of the central region of the cline, we suspect that this estimate is somewhat biased by selection on loci closely linked to the allozyme markers. A common question for taxa that show fixed allozyme differences in parapatry is whether or not they are genetically isolated. A fairly general measure of genetic isolation across hybrid zones is the time, T, that it takes a neutral allele to cross the hybrid zone and recombine into the opposite genetic background, given by T = (β/σ)², where β is the barrier strength of the hybrid zone. Genetic isolation in the Pontia zone is weak, with T  25 generations for most allozyme markers. By this measure, populations of daplidice and edusa on opposite sides of the hybrid zone share more identical-by-descent alleles than do populations of phenotypically pure daplidice in, say, France and Morocco. Accordingly, we think it best for systematists to consider edusa as a well-marked subspecies of P. daplidice.     

Inferences from a rapidly moving hybrid zone

Anartia fatima and Anartia amathea (Lepidoptera: Nymphalidae) are sister taxa whose ranges abut in a narrow hybrid zone in eastern Panama. At the center of the zone, hybrids are abundant, although deviations from Hardy-Weinberg and linkage disequilibria are strong, due in part to assortative mating. We measured differences across the zone in four wing color-pattern characters, three allozyme loci, and mitochondrial haplotype. Wing pattern, allozyme, and mitochondrial clines were coincident (i.e., had the same positions) and concordant (i.e.. all markers had similar cline shapes, about 28 km wide). Repeated samples demonstrated that the hybrid zone has been moving eastwards at an average rate of 2.5 km/year over the past 20 years, accompanied by an equivalent movement of the mtDNA cline. No introgression of mtDNA haplotypes were found in the 'wake" of the moving cline, as might be expected for a neutral marker. The concordance of morphological and mtDNA clines between 1994 and 2000, in spite of hybrid zone movement, suggests strong epistasis between the mitochondrial genome and nuclear loci. Cline movement is achieved mainly by pure fatima immigrating into amathea populations; hybrids had little effect, and were presumably outcompeted by fitter pure fatima genotypes. This movement can be explained if random dispersal of 7-19 km x gen(1/2) is coupled with a competitive advantage to A. fatima genomes of 2-5%. Hybrid zone motion is equivalent to Phase III of Wright's shifting balance. Hybrid zone movement has rarely been considered likely in the past, but our results show that it may be more important in biogeography and evolution than generally realized.     

Differential patterns of introgression across the X chromosome in a hybrid zone between two species of house mice

A complete understanding of the speciation process requires the identification of genomic regions and genes that confer reproductive barriers between species. Empirical and theoretical research has revealed two important patterns in the evolution of reproductive isolation in animals: isolation typically arises as a result of disrupted epistatic interactions between multiple loci and these disruptions map disproportionately to the X chromosome. These patterns suggest that a targeted examination of natural gene flow between closely related species at X-linked markers with known positions would provide insight into the genetic basis of speciation. We take advantage of the existence of genomic data and a well-documented European zone of hybridization between two species of house mice, Mus domesticus and M. musculus, to conduct such a survey. We evaluate patterns of introgression across the hybrid zone for 13 diagnostic X-linked loci with known chromosomal positions using a maximum likelihood model. Interlocus comparisons clearly identify one locus with reduced introgression across the center of the hybrid zone, pinpointing a candidate region for reproductive isolation. Results also reveal one locus with high frequencies of M. domesticus alleles in populations on the M. musculus side of the zone, suggesting the possibility that positive selection may act to drive the spread of alleles from one species on to the genomic background of the other species. Finally, cline width and cline center are strongly positively correlated across the X chromosome, indicating that gene flow of the X chromosome may be asymmetrical. This study highlights the utility of natural populations of hybrids for mapping speciation genes and suggests that the middle of the X chromosome may be important for reproductive isolation between species of house mice.     

Little and large: body size and genetic clines in a New Zealand gecko (Woodworthia maculata) along a coastal transect

Clinal variation can result from primary differentiation or secondary contact and determining which of these two processes is responsible for the existence of a cline is not a trivial problem. Samples from a coastal transect of New Zealand geckos (Woodworthia maculata) identified for the first time a body size cline 7-10 km wide. The larger geckos are almost twice the mass of the small adult geckos. Clines in allele and haplotype frequency were found at two of the four genetic loci examined. Estimated width of the morphological cline was concordant with neither the narrower mtDNA cline (3-7 m) nor the wider nuclear cline (RAG-2; 34-42 km), and cline centers were not coincident. Although the body size cline is narrow compared to the entire range of the species, it is 2-3 orders of magnitude greater than estimates of dispersal distance per generation for these geckos. No evidence of assortative mating, nor of hybrid disadvantage was identified, thus there is little evidence to infer that endogenous selection is maintaining a hybrid zone. We cannot distinguish secondary contact from primary origin of this body size cline but conclude that secondary contact is likely due to the occurrence of mtDNA haplotypes from three distinct clades within the coastal transect and the presence of two frequency clines within this region.     

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.     

Asymmetric size and shape variation in the Central European transect across the house mouse hybrid zone

We studied asymmetric variation of the mandible in the Central European portion of the hybrid zone between two house mouse subspecies, Mus musculus musculus and Mus musculus domesticus. Within introgression classes, defined by the share of diagnostic allozymes, we quantified the directional and fluctuating component of asymmetric variation, as well as skewness and kurtosis of individual asymmetry distributions. Furthermore, in the same manner we re‐analysed asymmetric variation of the ventral side of the skull. According to the quadratic polynomial model, the mandible shape‐fluctuating asymmetry, but not size‐fluctuating asymmetry, was significantly decreased in the centre of the hybrid zone (with a minimum predicted for a hybrid index of 0.41). On the contrary, the skull shape‐fluctuating asymmetry non‐monotonically increased towards the musculus side of the hybrid zone (with a peak predicted for a hybrid index of 0.86). Thus, the impact of hybridization on fluctuating asymmetry is trait‐specific in this portion of the house mouse hybrid zone. The only general feature of asymmetric variation we observed was the shift towards the platykurtosis of asymmetry distributions in the centre of the hybrid zone. Taken together, we suggest genetic variability for right–left asymmetries to be generally increased, but the developmental instability of mandible shape to be decreased, by hybridization. We hypothesize the decrease of developmental instability to be caused by overdominant effects on developmental dynamics rather than by increased heterozygosity. © 2010 The Linnean Society of London, Biological Journal of the Linnean Society, 2010, 101 , 13–27.     

An analysis of signatures of selective sweeps in natural populations of the house mouse

Population and locus-specific reduction of variability of polymorphic loci could be an indication of positive selection at a linked site (selective sweep) and therefore point toward genes that have been involved in recent adaptations. Analysis of microsatellite variability offers a way to identify such regions and to ask whether they occur more often than expected by chance. We studied four populations of the house mouse (Mus musculus) to assess the frequency of such signatures of selective sweeps under natural conditions. Three samples represent the subspecies Mus m. dometicus [corrected] and came from Germany, France, and Cameroon. One sample came from Kazakhstan and constitutes a population of the subspecies Mus m. [corrected] musculus. Mitochondrial D-loop sequences from all animals confirm their respective assignments. Approximately 200 microsatellite loci were typed for up to 60 unrelated individuals from each population and evaluated for signs of selective sweeps on the basis of Schl?tterer's ln RV and ln RH statistics. Our data suggest that there are slightly more signs of selective sweeps than would have been expected by chance alone in each of the populations and also highlights some of the statistical challenges faced in genome scans for detecting selection. Single-nucleotide polymorphism typing of one sweep signature in the M. m. domesticus populations around the beta-defensin 6 locus confirms a lowered nucleotide diversity in this region and limits the potential sweep region to about 20 kb. However, no amino acid exchange has occurred in the coding region when compared to M. m. musculus. If this sweep signature is due to a recent adaptation, it is expected that a regulatory change would have caused it. Our data provide a framework for conducting a systematic whole genome scan for signatures of selective sweeps in the mouse genome.     

There is no heterotic effect upon developmental stability in the ventral side of the skull within the house mouse hybrid zone

Fluctuating asymmetry (FA) of a body is expected to be related to an organism's developmental instability. We studied patterns of FA in the ventral side of the skull along a transect across the central-European portion of the hybrid zone between two house mouse subspecies, Mus musculus musculus and Mus musculus domesticus. We found that FA was not significantly different between introgression classes relative to differences between localities within those classes and that the within-class differences were much larger than differences between individual-signed asymmetries within localities. However, if year was added to the same analysis as another factor, FA was not significantly different among localities within the same introgression class. When individual asymmetries were plotted against individual hybrid indices, hybrids appeared more asymmetric than individuals from outside of the zone. Thus contrary to previous studies, we did not find lower FA indicating heterotic effect in hybrids for the traits studied. It is suggested that the impact of hybridization on FA in the ventral side of the mouse skull is negligible or overwhelmed by other factors.     

Estimates of Selection and Gene Flow from Measures of Cline Width and Linkage Disequilibrium in Heliconius Hybrid Zones

Hybrid zones can yield estimates of natural selection and gene flow. The width of a cline in gene frequency is approximately proportional to gene flow (sigma) divided by the square root of per-locus selection (square root of s). Gene flow also causes gametic correlations (linkage disequilibria) between genes that differ across hybrid zones. Correlations are stronger when the hybrid zone is narrow, and rise to a maximum roughly equal to s. Thus cline width and gametic correlations combine to give estimates of gene flow and selection. These indirect measures of sigma and s are especially useful because they can be made from collections, and require no field experiments. The method was applied to hybrid zones between color pattern races in a pair of Peruvian Heliconius butterfly species. The species are Müllerian mimics of one another, and both show the same changes in warning color pattern across their respective hybrid zones. The expectations of cline width and gametic correlation were generated using simulations of clines stabilized by strong frequency-dependent selection. In the hybrid zone in Heliconius erato, clines at three major color pattern loci were between 8.5 and 10.2 km wide, and the pairwise gametic correlations peaked at R approximately 0.35. These measures suggest that s approximately 0.23 per locus, and that sigma approximately 2.6 km. In erato, the shapes of the clines agreed with that expected on the basis of dominance. Heliconius melpomene has a nearly coincident hybrid zone. In this species, cline widths at four major color pattern loci varied between 11.7 and 13.4 km. Pairwise gametic correlations peaked near R approximately 1.00 for tightly linked genes, and at R approximately 0.40 for unlinked genes, giving s approximately 0.25 per locus and sigma approximately 3.7 km. In melpomene, cline shapes did not perfectly fit theoretical shapes based on dominance; this deviation might be explained by long-distance migration and/or strong epistasis. Compared with erato, sample sizes in melpomene are lower and the genetics of its color patterns are less well understood. In spite of these problems, selection and gene flow are clearly of the same order of magnitude in the two species. The relatively high per locus selection coefficients agree with "major gene" theories for the evolution of Müllerian mimicry, but the genetic architecture of the color patterns does not. These results show that the genetics and evolution of mimicry are still only sketchily understood.     

Assessing multilocus introgression patterns: a case study on the mouse X chromosome in central Europe

Multilocus hybrid zone (HZ) studies predate genomics by decades. The power of early methods is becoming apparent and now large datasets are commonplace. Relating introgression along a chromosome to evolutionary process is challenging: although reduced introgression regions may indicate speciation genes, this pattern may be obscured by asymmetric introgression of linked invasive genes. Further, HZ movement may form salients and leave islands in its wake. Barton's concordance was proposed 24 years ago for assessing introgression where geographic patterns are complex. The geographic axis of introgression is replaced with the hybrid index. We compare this, a recently proposed genomic clines approach, and two-dimensional (2D) geographic analyses, for 24 X chromosome loci of 2873 mice from the central-European house mouse HZ. In 2D, 14 loci show linear contact, seven precisely matching previous studies. Four show introgression islands to the east of the zone, suggesting past westward zone movement, two show westward salients. Barton's concordance both recovers and refines this information. A region of reduced introgression on the central X is supported, despite X centromere-proximal male-biased westward introgression matching a westward 2D geographic salient. Genomic clines results are consistent regarding introgression asymmetries, but otherwise more difficult to interpret. Evidence for genetic conflict is discussed.     

Estimation of the dispersal of a major pest of maize by cline analysis of a temporary contact zone between two invasive outbreaks

Dispersal is a key factor in invasion and in the persistence and evolution of species. Despite the importance of estimates of dispersal distance, dispersal measurement remains a real methodological challenge. In this study, we characterized dispersal by exploiting a specific case of biological invasion, in which multiple introductions in disconnected areas lead to secondary contact between two differentiated expanding outbreaks. By applying cline theory to this ecological setting, we estimated σ, the standard deviation of the parent–offspring distance distribution, of the western corn rootworm, Diabrotica virgifera virgifera, one of the most destructive pests of maize. This species is currently invading Europe, and the two largest invasive outbreaks, in northern Italy and Central Europe, have recently formed a secondary contact zone in northern Italy. We identified vanishing clines at 12 microsatellite loci throughout the contact zone. By analysing both the rate of change of cline slope and the spatial variation of linkage disequilibrium at these markers, we obtained two σ estimates of about 20 km/generation^1/2. Simulations indicated that these estimates were robust to changes in dispersal kernels and differences in population density between the two outbreaks, despite a systematic weak bias. These estimates are consistent with the results of direct methods for measuring dispersal applied to the same species. We conclude that secondary contact resulting from multiple introductions is very useful for the inference of dispersal parameters and should be more widely used in other species.