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.