The genetic dissection of naturally occurring phenotypes sheds light on many fundamental and longstanding questions in speciation and adaptation and is a central research topic in evolutionary biology. Until recently, forward-genetic approaches were virtually impossible to apply to nonmodel organisms, but the development of next-generation sequencing techniques eases this difficulty. Here, we use the ddRAD-seq method to map a colour trait with a known adaptive function in cichlid fishes, well-known textbook examples for rapid rates of speciation and astonishing phenotypic diversification. A suite of phenotypic key innovations is related to speciation and adaptation in cichlids, among which body coloration features prominently. The focal trait of this study, horizontal stripes, evolved in parallel in several cichlid radiations and is associated with piscivorous foraging behaviour. We conducted interspecific crosses between Haplochromis sauvagei and H. nyererei and constructed a linkage map with 867 SNP markers distributed on 22 linkage groups and total size of 1130.63 cM. Lateral stripes are inherited as a Mendelian trait and map to a single genomic interval that harbours a paralog of a gene with known function in stripe patterning. Dorsolateral and mid-lateral stripes were always coinherited and are thus under the same genetic control. Additionally, we directly quantify the genotyping error rates in RAD markers and offer guidelines for identifying and dealing with errors. Uncritical marker selection was found to severely impact linkage map construction. Fortunately, by applying appropriate quality control steps, a genotyping accuracy of >99.9% can be reached, thus allowing for efficient linkage mapping of evolutionarily relevant traits.

Individuals of the scale-eating cichlid fish, Perissodus microlepis, from Lake Tanganyika tend to have remarkably asymmetric heads that are either left-bending or right-bending. The 'left' morph opens its mouth markedly towards the left and preferentially feeds on the scales from the right-hand side of its victim fish, and the 'right' morph bites scales from the victims' left-hand side. This striking dimorphism made these fish a textbook example of their astonishing degree of ecological specialization and as one of the few known incidences of negative frequency-dependent selection acting on an asymmetric morphological trait, where left and right forms are equally frequent within a species. We investigated the degree and the shape of the frequency distribution of head asymmetry in P. microlepis to test whether the variation conforms to a discrete dimorphism, as generally assumed. In both adult and juvenile fish, mouth asymmetry appeared to be continuously and unimodally distributed with no clear evidence for a discrete dimorphism. Mixture analyses did not reveal evidence of a discrete or even strong dimorphism. These results raise doubts about previous claims, as reported in textbooks, that head variation in P. microlepis represents a discrete dimorphism of left-and right-bending forms. Based on extensive field sampling that excluded ambiguous (i.e. symmetric or weakly asymmetric) individual adults, we found that left and right morphs occur in equal abundance in five populations. Moreover, mate pairing for 51 wild-caught pairs was random with regard to head laterality, calling into question reports that this laterality is maintained through disassortative mating.

Scale-eating cichlid fish, Perissodus microlepis, from Lake Tanganyika display handed (lateralized) foraging behavior, where an asymmetric ‘left’ mouth morph preferentially feeds on the scales of the right side of its victim fish and a ‘right’ morph bites the scales of the left side. This species has therefore become a textbook example of the astonishing degree of ecological specialization and negative frequency-dependent selection. We investigated the strength of handedness of foraging behavior as well as its interaction with morphological mouth laterality in P. microlepis. In wild-caught adult fish we found that mouth laterality is, as expected, a strong predictor of their preferred attack orientation. Also laboratory-reared juvenile fish exhibited a strong laterality in behavioral preference to feed on scales, even at an early age, although the initial level of mouth asymmetry appeared to be small. This suggests that pronounced mouth asymmetry is not a prerequisite for handed foraging behavior in juvenile scale-eating cichlid fish and might suggest that behavioral preference to attack a particular side of the prey plays a role in facilitating morphological asymmetry of this species.

Population genetic analyses were conducted to investigate whether random mating occurs between left and right-mouth morphs of the dimorphic scale-eating cichlid fish Perissodus microlepis from two geographical sites in southern Lake Tanganyika. The mitochondrial and nuclear DNA markers (13 microsatellite loci) revealed no genetic differentiation between left and right morphs (i.e. widespread interbreeding). The observed lack of genetic divergence between the different morphs allowed for the exclusion of the possibility of assortative mating between same morph types. The microsatellite data showed no significant departures of heterozygosity from Hardy-Weinberg equilibrium, suggesting purely random mating between the morphs. Overall, this study indicated no genetic evidence for either assortative or disassortative mating, but it did provide support for the random mating hypothesis. Highly significant, albeit weak, spatial population structure was also found when samples of different morphs were pooled according to geographical sites. An additional analysis of two microsatellite loci that were recently suggested to be putatively linked to the genetic locus that determines the laterality of these mouth morphs did not show any such association.

We investigated the effect of development mode on the spatial and temporal population genetic structure of four littorinid gastropod species. Snails were collected from the same three sites on the west coast of Vancouver Island, Canada in 1997 and again in 2007. DNA sequences were obtained for one mitochondrial gene, cytochrome b (Cyt b), and for up to two nuclear genes, heat shock cognate 70 (HSC70) and aminopeptidase N intron (APN54). We found that the mean level of genetic diversity and long-term effective population sizes (N(e)) were significantly greater for two species, Littorina scutulata and L. plena, that had a planktotrophic larval stage than for two species, Littorina sitkana and L. subrotundata, that laid benthic egg masses which hatched directly into crawl-away juveniles. Predictably, two poorly dispersing species, L. sitkana and L. subrotundata, showed significant spatial genetic structure at an 11- to 65-km geographical scale that was not observed in the two planktotrophic species. Conversely, the two planktotrophic species had more temporal genetic structure over a 10-year interval than did the two direct-developing species and showed highly significant temporal structure for spatially pooled samples. The greater temporal genetic variation of the two planktotrophic species may have been caused by their high fecundity, high larval dispersal, and low but spatially correlated early survivorship. The sweepstakes-like reproductive success of the planktotrophic species could allow a few related females to populate hundreds of kilometres of coastline and may explain their substantially larger temporal genetic variance but lower spatial genetic variance relative to the direct-developing species.

The introduction of invasive Nile tilapia (Oreochromis niloticus), and the rapacious predator Nile perch (Lates niloticus), into Lake Victoria resulted in a decline in population sizes, genetic diversity and even extirpation of native species which were previously the mainstay of local fisheries. However, remnant populations of native fish species, including tilapia, still persist in satellite lakes around Lake Victoria where they may coexist with O. niloticus. In this study we assessed population genetic structure, diversity, and integrity of the native critically endangered Singidia tilapia (O. esculentus) in its refugial populations in the Yala swamp, Kenya, and contrasted this diversity with populations of the invasive tilapia O. niloticus in satellite lakes (Kanyaboli, Namboyo and Sare) and Lake Victoria. Based on mtDNA control region sequences and eight nuclear microsatellite loci, we did not detect any mtDNA introgression between the native and the invasive species in Lakes Kanyaboli and Namboyo, but did find low levels of nuclear admixture, primarily from O. niloticus to O. esculentus. Some genetic signal of O. esculentus in O. niloticus was found in Lake Sare, where O. esculentus is not found, suggesting it has recently been extirpated by the O. niloticus invasion. In both species, populations in the satellite lakes are significantly genetically isolated from each other, with private mtDNA haplotypes and microsatellite alleles. For O. niloticus, genetic diversity in satellite lakes was similar to that found in Lake Victoria. Our data imply a low frequency of immigration exchange between the two populations of O. esculentus and we suggest that the populations of this endangered species and important fisheries resource should be conserved separately in Lakes Kanyaboli and Namboyo and with high priority.

The scale‐eating cichlid fish, Perissodus microlepis, from Lake Tanganyika are a well‐known example of an asymmetry dimorphism because the mouth/head is either left‐bending or right‐bending. However, how strongly its pronounced morphological laterality is affected by genetic and environmental factors remains unclear. Using quantitative assessments of mouth asymmetry, we investigated its origin by estimating narrow‐sense heritability (h 2) using midparent–offspring regression. The heritability estimates [field estimate: h 2 = 0.22 ± 0.06, P = 0.013; laboratory estimate: h 2 = 0.18 ± 0.05, P = 0.004] suggest that although variation in laterality has some additive genetic component, it is strongly environmentally influenced. Family‐level association analyses of a putative microsatellite marker that was claimed to be linked to gene(s) for laterality revealed no association of this locus with laterality. Moreover, the observed phenotype frequencies in offspring from parents of different phenotype combinations were not consistent with a previously suggested single‐locus two‐allele model, but they neither were able to reject with confidence a random asymmetry model. These results reconcile the disputed mechanisms for this textbook case of mouth asymmetry where both genetic and environmental factors contribute to this remarkable case of morphological asymmetry.

The present population structure of a species reflects the influence of population history as well as contemporary processes. To examine the relative importance of these factors in shaping the current population structure of Littorina keenae, we sequenced 762 base pairs of the mitochondrial ND6 and cytochrome b genes in 584 snails from 13 sites along the northeastern Pacific coast. Haplotype network analysis revealed a 'star-like' genealogy indicative of a recent population expansion. Nested clade and mismatch analyses also supported the hypothesis of sudden population expansion following a population bottleneck during the Last Glacial Maximum. Analysis of molecular variance and pairwise Phi(ST) showed no significant spatial population differentiation from Mexico to Oregon - not even across the recognized biogeographic boundary at Point Conception. This is probably due to high contemporary gene flow during the free-swimming larval stage of this snail. Surprisingly, we found a highly significant temporal population differentiation between a San Pedro sample from 1996 and one from 2005, which gave an estimate of effective population size (N(e)) of only 135. Nearly statistically significant changes in the frequency of a particular haplotype in three other populations over 2-3 years further support Hedgecock's 'sweepstakes' hypothesis. When by chance nearly all of the progeny from an aggregation of highly fecund sisters that possess a rare haplotype successfully recruit to become the next generation, the rare haplotype can become temporarily common across the entire species' range. This modification of the sweepstakes hypothesis can explain why the temporal variation that we observed was much greater than the spatial variation.