Evidence for a bimodal distribution of hybrid indices in a hybrid zone with high admixture

The genetic structure of a hybrid zone can provide insights into the relative roles of the various factors that maintain the zone. Here, we use a multilocus approach to characterize a hybrid zone between two subspecies of killifish (Fundulus heteroclitus, Walbaum 1792) found along the Atlantic coast of North America. We first analysed clinal variation along the Atlantic coast using a single-nucleotide polymorphism in the mitochondrial DNA (mtDNA) displacement loop (D-loop) and a panel of nine nuclear microsatellite markers. A model constraining all clines to the same width and centre was not significantly different from a model in which the clines were allowed to vary independently. Locus-by-locus analysis indicated that the majority of nuclear clines shared the same centre as the mtDNA cline, and the widths of these clines were also narrower than that predicted by a neutral model, suggesting that selection is operating to maintain the hybrid zone. However, two of the nuclear clines had widths greater than the neutral prediction and had centres that were displaced relative to the mtDNA cline centre. We also found that a marsh located near the centre of the mtDNA cline demonstrated a bimodal distribution of nuclear hybrid index values, suggesting a deficit of first-generation hybrids and back-crossed genotypes. Thus, selection against hybrid genotypes may be playing a role in maintaining this hybrid zone and the associated steep nuclear and mtDNA clines.