Complex patterns of genetic and phenotypic divergence in populations of the Lake Malawi cichlid Maylandia zebra

Ecological speciation, which relies heavily on selection driving the emergence of new species, has become the primary paradigm through which rapid species radiations are understood. In this way, selection, particularly ecological selection, is assumed to be the driver of most species radiations. However, in many radiations, such as the radiation of Lake Malawi’s cichlids, the assumption of selection as the driver of speciation has rarely been explicitly tested, and drift, often, has completely been ignored as potential factor. In order to understand the forces driving the divergence of Lake Malawi cichlids at the microevolutionary level, we studied the genetic and phenotypic divergence of ten allopatric populations of Maylandia zebra. We estimated effective population sizes as proxy for drift. Further, we compared neutral genetic differentiation to divergence in three phenotypic traits: body size, body shape, and melanophore counts. We found small, yet significant, population differentiation in all the studied traits across most populations. Population sizes were small rendering the potential for drift to be high. However, phenotypic differentiation exceeded neutral expectations for all traits suggesting divergent local selection. Our data suggest that natural, and potentially also sexual, selection may be the dominant force driving population differentiation in Lake Malawi’s rock-dwelling cichlids, despite the potential for drift in small populations.