To understand mechanisms structuring diversity in young adaptive radiations, quantitative and unbiased information about genetic and phenotypic diversity is much needed. Here, we present the first in-depth investigation of whitefish diversity in a Swiss lake, with continuous spawning habitat sampling in both time and space. Our results show a clear cline like pattern in genetics and morphology of populations sampled along an ecological depth gradient in Lake Neuchâtel. Divergent natural selection appears to be involved in shaping this cline given that trait specific P(ST)-values are significantly higher than F(ST)-values when comparing populations caught at different depths. These differences also tend to increase with increasing differences in depth, indicating adaptive divergence along a depth gradient, which persists despite considerable gene flow between adjacent demes. It however remains unclear, whether the observed pattern is a result of currently stable selection-gene flow balance, incipient speciation, or reverse speciation due to anthropogenic habitat alteration causing two formerly divergent species to collapse into a single gene pool.