Spatial Geographic Mosaic in an Aquatic Predator-Prey Network

The geographic mosaic theory of coevolution predicts 1) spatial variation in predatory structures as well as prey defensive traits, and 2) trait matching in some areas and trait mismatching in others mediated by gene flow. We examined gene flow and documented spatial variation in crushing resistance in the freshwater snails Mexipyrgus churinceanus, Mexithauma quadripaludium, Nymphophilus minckleyi, and its relationship to the relative frequency of the crushing morphotype in the trophically polymorphic fish Herichthys minckleyi. Crushing resistance and the frequency of the crushing morphotype did show spatial variation among 11 naturally replicated communities in the Cuatro Ciénegas valley in Mexico where these species are all endemic. The variation in crushing resistance among populations was not explained by geographic proximity or by genetic similarity in any species. We detected clear phylogeographic patterns and limited gene flow for the snails but not for the fish. Gene flow among snail populations in Cuatro Ciénegas could explain the mosaic of local divergence in shell strength and be preventing the fixation of the crushing morphotype in Herichthys minckleyi. Finally, consistent with trait matching across the mosaic, the frequency of the fish morphotype was negatively correlated with shell crushing resistance likely reflecting the relative disadvantage of the crushing morphotype in communities where the snails exhibit relatively high crushing resistance.