Ecosystem size matters: The dimensionality of intralacustrine diversification in Icelandic stickleback is predicted by lake size

Cases of evolutionary diversification can be characterized along a continuum from weak to strong genetic and phenotypic differentiation. Several factors may facilitate or constrain the differentiation process. Comparative analyses of replicates of the same taxon at different stages of differentiation can be useful to identify these factors. We estimated the number of distinct phenotypic groups in threespine stickleback populations from nine lakes in Iceland and in one marine population. Using the inferred number of phenotypic groups in each lake, genetic divergence from the marine population, and physical lake and landscape variables, we tested if ecosystem size, approximated by lake size and depth, or isolation from the ancestral marine gene pool predict the occurrence and the extent of phenotypic and genetic diversification within lakes. We find intralacustrine phenotypic diversification to be the rule rather than the exception, occurring in all but the youngest lake population and being manifest in ecologically important phenotypic traits. Neutral genetic data further indicates non-random mating in four out of nine studied lakes, and restricted gene flow between sympatric phenotypic groups in two. Although neither the phenotypic variation nor the number of intralacustrine phenotypic groups were associated with any of our environmental variables, the number of phenotypic traits that were differentiated was significantly positively related to lake size, and evidence for restricted gene flow between sympatric phenotypic groups was only found in the largest lakes where trait specific phenotypic differentiation was highest.