Asymmetric reproductive isolation between terminal forms of the salamander ring species Ensatina eschscholtzii revealed by fine-scale genetic analysis of a hybrid zone

<p>Abstract</p> <p>Background</p> <p>Ring species, exemplified by salamanders of the <it>Ensatina eschscholtzii </it>complex, represent a special window into the speciation process because they allow the history of species formation to be traced back in time through the geographically differentiated forms connecting the two terminal forms of the ring. Of particular interest is the nature and extent of reproductive isolation between the geographically terminal forms, in this case <it>E. e. eschscholtzii </it>and <it>E. e. klauberi</it>. Previous studies have documented infrequent hybridization at the end of the ring. Here, we report the first fine-scale genetic analysis of a hybrid zone between the terminal forms in southern California using individual-based Bayesian analyses of multilocus genetic data to estimate levels and direction of hybridization and maximum-likelihood analysis of linkage disequilibrium and cline shape to make inferences about migration and selection in the hybrid zone.</p> <p>Results</p> <p>The center of the hybrid zone has a high proportion of hybrids, about half of which were classified as F1s. Clines are narrow with respect to dispersal, and there are significant deviations from Hardy-Weinberg equilibrium as well as nonrandom associations (linkage disequilibria) between alleles characteristic of each parental type. There is cytonuclear discordance, both in terms of introgression and the geographic position of mitochondrial versus nuclear clines. Genetic disequilibrium is concentrated on the <it>eschscholtzii </it>side of the zone. Nearly all hybrids possess <it>klauberi </it>mtDNA, indicating that most hybrids are formed from female <it>klauberi </it>mating with male <it>eschscholtzii </it>or male hybrids (but not vice versa).</p> <p>Conclusions</p> <p>Our results are consistent with a tension zone trapped at an ecotone, with gene combinations characteristic of <it>klauberi </it>showing up on the <it>eschscholtzii </it>side of the zone due to asymmetric hybridization. We suggest that the observed asymmetry is best explained by increased discriminatory power of <it>eschscholtzii </it>females, or asymmetric postzygotic isolation. The relatively high frequency of hybrids, particularly F1s, contrasts with other contacts between the terminal forms, and with other contacts between other divergent <it>Ensatina </it>lineages, highlighting the diverse outcomes of secondary contact within a single species complex.</p