GENE FLOW AND MITOCHONDRIAL DNA VARIATION IN THE KILLIFISH, FUNDULUS HETEROCLITUS

Two subspecies of killifish, Fundulus heteroclitus, inhabit the Atlantic coast of the United States. The contact zones between them are typified by morphological, physiological, nuclear gene, and mitochondrial DNA clines. Considerable debate exists about the importance of restricted gene flow and natural selection in maintaining these clines. To evaluate the relative importance of these two evolutionary forces we employed analysis of mtDNA as an independent measure of gene flow. Solutions of equilibrium diffusion equations suggested that killifish dispersal must be less than 2 km to maintain previously observed allozyme frequency clines in the absence of strong selection. To determine whether populations separated by distances greater than 2 km accumulate significant genetic differences we examined a total of N = 480 individuals from five killifish populations spanning a continuous shoreline distance of 8.4 km. Distribution of the 25 detected haplotypes indicated that most of the variation was within rather than between sampling locations. No evidence of geographic structure was detected, nor were there any significant genetic differences between killifish populations. With these data, gene flow was evaluated by estimating effective migration rate (N_e m) between the populations. N_e m was estimated from G_st (N_e m = 49.4), from F_st (N_e m = 24.1), by the private alleles approach (N_e m = 18.5) and via four phylogenetic analyses (N_e m ranged from 11.4 to 16.9). Regardless of the analysis, N_e m was greater than one; the threshold level at which gene flow is strong enough to prevent differentiation due to genetic drift. This suggests that while the characteristics that distinguish Fundulus subspecies may have arisen in isolation, the current clinal distributions exhibited along the Atlantic coast are most likely maintained by selection.