Understanding the origin, maintenance and significance of phenotypic variation is one of the central issues in evolutionary biology. An ongoing discussion focuses on the relative roles of isolation and selection as being at the heart of genetically based spatial variation. We address this issue in a representative of a taxon group in which isolation is unlikely: a marine broadcast spawning invertebrate. During the free-swimming larval phase, dispersal is potentially very large. For such taxa, small-scale population genetic structuring in neutral molecular markers tends to be limited, conform expectations. Small-scale differentiation of selective traits is expected to be hindered by the putatively high gene flow. We determined the geographical distribution of molecular markers and of variation in a shell shape measure, globosity, for the bivalve Macoma balthica (L.) in the western Dutch Wadden Sea and adjacent North Sea in three subsequent years, and found that shells of this clam are more globose in the Wadden Sea. By rearing clams in a common garden in the laboratory starting from the gamete phase, we show that the ecotypes are genetically different; heritability is estimated at 23%. The proportion of total genetic variation that is between sites is much larger for the morphological additive genetic variation (Q(ST) = 0.416) than for allozyme (F-ST = 0.000-0.022) and mitochondrial DNA cytochrome-c-oxidase-1 sequence variation (Phi(ST) = 0.017). Divergent selection must be involved and intraspecific spatial genetic differentiation in marine broadcast spawners is apparently not constrained by low levels of isolation.