Female mate choice has often been proposed to play an important role in cases of rapid speciation, in particular in the explosively evolved haplochromine cichlid species flocks of the Great Lakes of East Africa. Little, if anything, is known in cichlid radiations about the heritability of female mating preferences. Entirely sympatric distribution, large ecological overlap and conspicuous differences in male nuptial coloration, and female preferences for these, make the sister species Pundamilia pundamilia and P. nyererei from Lake Victoria an ideally suited species pair to test assumptions on the genetics of mating preferences made in models of sympatric speciation. Female mate choice is necessary and sufficient to maintain reproductive isolation between these species, and it is perhaps not unlikely therefore, that female mate choice has been important during speciation. A prerequisite for this, which had remained untested in African cichlid fish, is that variation in female mating preferences is heritable. We investigated mating preferences of females of these sister species and their hybrids to test this assumption of most sympatric speciation models, and to further test the assumption of some models of sympatric speciation by sexual selection that female preference is a single-gene trait. We find that the differences in female mating preferences between the sister species are heritable, possibly with quite high heritabilities, and that few but probably more than one genetic loci contribute to this behavioural speciation trait with no apparent dominance. We discuss these results in the light of speciation models and the debate about the explosive radiation of cichlid fishes in Lake Victoria.

No abstract

Major environmental events that fragment populations among multiple island habitats have potential to drive large-scale episodes of speciation and adaptive radiation. A recent palaeolimnological study of sediment cores indicated that Lake Malawi underwent major climate-driven desiccation events 75,000-135,000 years ago that lowered the water level to at least 580 m below the present state and severely reduced surface area. After this period, lake levels rose and stabilized, creating multiple discontinuous littoral rocky habitats. Here, we present evidence supporting the hypothesis that establishment and expansion of isolated philopatric rock cichlid populations occurred after this rise and stabilization of lake level. We studied the Pseudotropheus (Maylandia) species complex, a group with both allopatric and sympatric populations that differ in male nuptial colour traits and tend to mate assortatively. Using coalescent analyses based on mitochondrial DNA, we found evidence that populations throughout the lake started to expand and accumulate genetic diversity after the lake level rise. Moreover, most haplotypes were geographically restricted, and the greatest genetic similarities were typically among sympatric or neighbouring populations. This is indicative of limited dispersal and establishment of assortative mating among populations following the lake level rise. Together, this evidence is compatible with a single large-scale environmental event being central to evolution of spatial patterns of genetic and species diversity in P. (Maylandia) and perhaps other Lake Malawi rock cichlids. Equivalent climate-driven pulses of habitat formation and fragmentation may similarly have contributed to observed rapid and punctuated cladogenesis in other adaptive radiations.

Genetic linkage maps are essential for comparative genomics, high quality genome sequence assembly and fine scale quantitative trait locus (QTL) mapping. In the present study we identified and genotyped markers via restriction-site associated DNA (RAD) sequencing and constructed a genetic linkage map based on 1,597 SNP markers of an interspecific F2 cross of two closely related Lake Victoria cichlids (Pundamilia pundamilia and P. sp. ‘red head’). The SNP markers were distributed on 22 linkage groups and the total map size was 1,594 cM with an average marker distance of 1.01 cM. This high-resolution genetic linkage map was used to anchor the scaffolds of the Pundamilia genome and estimate recombination rates along the genome. Via QTL mapping we identified a major QTL for sex in a ∼1.9 Mb region on Pun-LG10, which is homologous to Oreochromis niloticus LG 23 (Ore-LG23) and includes a well-known vertebrate sex-determination gene (amh).