Complex histories of repeated gene flow in Cameroon crater lake cichlids cast doubt on one of the clearest examples of sympatric speciation

One of the most celebrated examples of sympatric speciation in nature are monophyletic radiations of cichlid fishes endemic to Cameroon crater lakes. However, phylogenetic inference of monophyly may not detect complex colonization histories involving some allopatric isolation, such as double invasions obscured by genome‐wide gene flow. Population genomic approaches are better suited to test hypotheses of sympatric speciation in these cases. Here, we use comprehensive sampling from all four sympatric crater lake cichlid radiations in Cameroon and outgroups across Africa combined with next‐generation sequencing to genotype tens of thousands of SNPs. We find considerable evidence of gene flow between all four radiations and neighboring riverine populations after initial colonization. In a few cases, some sympatric species are more closely related to outgroups than others, consistent with secondary gene flow facilitating their speciation. Our results do not rule out sympatric speciation in Cameroon cichlids, but rather reveal a complex history of speciation with gene flow, including allopatric and sympatric phases, resulting in both reproductively isolated species and incipient species complexes. The best remaining non‐cichlid examples of sympatric speciation all involve assortative mating within microhabitats. We speculate that this feature may be necessary to complete the process of sympatric speciation in nature.     

Genomic architecture of ecologically divergent body shape in a pair of sympatric crater lake cichlid fishes

Determining the genetic bases of adaptations and their roles in speciation is a prominent issue in evolutionary biology. Cichlid fish species flocks are a prime example of recent rapid radiations, often associated with adaptive phenotypic divergence from a common ancestor within a short period of time. In several radiations of freshwater fishes, divergence in ecomorphological traits — including body shape, colour, lips and jaws — is thought to underlie their ecological differentiation, specialization and, ultimately, speciation. The Midas cichlid species complex (Amphilophus spp.) of Nicaragua provides one of the few known examples of sympatric speciation where species have rapidly evolved different but parallel morphologies in young crater lakes. This study identified significant QTL for body shape using SNPs generated via ddRAD sequencing and geometric morphometric analyses of a cross between two ecologically and morphologically divergent, sympatric cichlid species endemic to crater Lake Apoyo: an elongated limnetic species (Amphilophus zaliosus) and a high‐bodied benthic species (Amphilophus astorquii). A total of 453 genome‐wide informative SNPs were identified in 240 F₂ hybrids. These markers were used to construct a genetic map in which 25 linkage groups were resolved. Seventy‐two segregating SNPs were linked to 11 QTL. By annotating the two most highly supported QTL‐linked genomic regions, genes that might contribute to divergence in body shape along the benthic–limnetic axis in Midas cichlid sympatric adaptive radiations were identified. These results suggest that few genomic regions of large effect contribute to early stage divergence in Midas cichlids.     

A phylogenetic analysis and taxonomic revision of ricefishes,Oryzias and relatives (Beloniformes,Adrianichthyidae)

Ricefishes, known best by the model organism, the medaka, Oryzias latipes Temminck & Schlegel, 1846, comprise the family Adrianichthyidae, which ranges broadly throughout fresh and brackish waters of Central, South and Southeast Asia and the Indo‐Malay‐Philippines Archipelago as far east as Timor. Twenty‐eight Recent species are recognized here in two monophyletic genera, Adrianichthys and Oryzias. Xenopoecilus and Horaichthys are placed in synonymy of Oryzias for the first time. Adrianichthys comprises four species from Lake Poso, Sulawesi, Indonesia. Oryzias comprises 24 species that live throughout the range of the family. A fossil genus and species, ?Lithopoecilus brouweri from the Miocene of central Sulawesi, is included tentatively in the Adrianichthyidae. Evidence for the sister group relationship of adrianichthyids and exocoetoids is reviewed briefly and that relationship corroborated. Monophyly of adrianichthyids is likewise strongly supported here. Species groups within Oryzias are diagnosed as monophyletic largely based on osteology, colour pattern and meristic variation. They correspond only in part to species groups previously recognized based on chromosome constitution. Miniature species do not comprise a monophyletic group; disjunct absolute size in close relatives has evolved repeatedly. Oryzias latipes is a member of a species complex that includes O. luzonensis, O. curvinotus and the miniatures O. sinensis and O. mekongensis. A new species, Oryzias bonneorum sp. nov. , is described from Lake Lindu, Sulawesi, Indonesia. Lectotypes are designated for Haplochilus celebensis Weber, 1894 and Haplochilus timorensis Weber & de Beaufort, 1922. No claim to original US Government works. Journal compilation © 2008 The Linnean Society of London, Zoological Journal of the Linnean Society, 2008, 154 , 494–610.     

The shape of things to come in the study of the origin of species?

Perhaps Darwin would agree that speciation is no longer the mystery of mysteries that it used to be. It is now generally accepted that evolution by natural selection can contribute to ecological adaptation, resulting in the evolution of reproductive barriers and, hence, to the evolution of new species (Schluter & Conte 2009 ; Meyer 2011 ; Nosil 2012 ). From genes that encode silencing proteins that cause infertility in hybrid mice (Mihola et al. 2009 ), to segregation distorters linked to speciation in fruit flies (Phadnis & Orr 2009 ), or pollinator‐mediated selection on flower colour alleles driving reinforcement in Texan wildflowers (Hopkins & Rausher 2012 ), characterization of the genes that drive speciation is providing clues to the origin of species (Nosil & Schluter 2011 ). It is becoming apparent that, while recent work continues to overturn historical ideas about sympatric speciation (e.g. Barluenga et al. 2006 ), ecological circumstances strongly influence patterns of genomic divergence, and ultimately the establishment of reproductive isolation when gene flow is present (Elmer & Meyer 2011 ). Less clear, however, are the genetic mechanisms that cause speciation, particularly when ongoing gene flow is occurring. Now, in this issue, Franchini et al. ( 2014 ) employ a classic genetic mapping approach augmented with new genomic tools to elucidate the genomic architecture of ecologically divergent body shapes in a pair of sympatric crater lake cichlid fishes. From over 450 segregating SNPs in an F2 cross, 72 SNPs were linked to 11 QTL associated with external morphology measured by means of traditional and geometric morphometrics. Annotation of two highly supported QTL further pointed to genes that might contribute to ecological divergence in body shape in Midas cichlids, overall supporting the hypothesis that genomic regions of large phenotypic effect may be contributing to early‐stage divergence in Midas cichlids.     

Are sympatrically speciating Midas cichlid fish special? Patterns of morphological and genetic variation in the closely related species Archocentrus centrarchus

Established empirical cases of sympatric speciation are scarce, although there is an increasing consensus that sympatric speciation might be more common than previously thought. Midas cichlid fish are one of the few substantiated cases of sympatric speciation, and they formed repeated radiations in crater lakes. In contrast, in the same environment, such radiation patterns have not been observed in other species of cichlids and other families of fish. We analyze morphological and genetic variation in a cichlid species (Archocentrus centrarchus) that co‐inhabits several crater lakes with the Midas species complex. In particular, we analyze variation in body and pharyngeal jaw shape (two ecologically important traits in sympatrically divergent Midas cichlids) and relate that to genetic variation in mitochondrial control region and microsatellites. Using these four datasets, we analyze variation between and within two Nicaraguan lakes: a crater lake where multiple Midas cichlids have been described and a lake where the source population lives. We do not observe any within‐lake clustering consistent across morphological traits and genetic markers, suggesting the absence of sympatric divergence in A. centrarchus. Genetic differentiation between lakes was low and morphological divergence absent. Such morphological similarity between lakes is found not only in average morphology, but also when analyzing covariation between traits and degree of morphospace occupation. A combined analysis of the mitochondrial control region in A. centrarchus and Midas cichlids suggests that a difference between lineages in the timing of crater lake colonization cannot be invoked as an explanation for the difference in their levels of diversification. In light of our results, A. centrarchus represents the ideal candidate to study the genomic differences between these two lineages that might explain why some lineages are more likely to speciate and diverge in sympatry than others.     

Genomic signatures of divergent selection and speciation patterns in a ‘natural experiment’, the young parallel radiations of Nicaraguan crater lake cichlid fishes

Divergent selection is the main driving force in sympatric ecological speciation and may also play a strong role in divergence between allopatric populations. Characterizing the genome‐wide impact of divergent selection often constitutes a first step in unravelling the genetic bases underlying adaptation and ecological speciation. The Midas cichlid fish (Amphilophus citrinellus) species complex in Nicaragua is a powerful system for studying evolutionary processes. Independent colonizations of isolated young crater lakes by Midas cichlid populations from the older and great lakes of Nicaragua resulted in the repeated evolution of adaptive radiations by intralacustrine sympatric speciation. In this study we performed genome scans on two repeated radiations of crater lake species and their great lake source populations (1030 polymorphic AFLPs, n ～ 30 individuals per species). We detected regions under divergent selection (0.3% in the crater lake Xiloá flock and 1.7% in the older crater lake Apoyo radiation) that might be responsible for the sympatric diversifications. We find no evidence that the same genomic regions have been involved in the repeated evolution of parallel adaptations across crater lake flocks. However, there is some genetic parallelism apparent (seven out of 51 crater lake to great lake outlier loci are shared; 13.7%) that is associated with the allopatric divergence of both crater lake flocks. Interestingly, our results suggest that the number of outlier loci involved in sympatric and allopatric divergence increases over time. A phylogeny based on the AFLP data clearly supports the monophyly of both crater lake species flocks and indicates a parallel branching order with a primary split along the limnetic‐benthic axis in both radiations.     

Population genomics of local adaptation versus speciation in coral reef fishes (Hypoplectrus spp,Serranidae)

Are the population genomic patterns underlying local adaptation and the early stages of speciation similar? Addressing this question requires a system in which (i) local adaptation and the early stages of speciation can be clearly identified and distinguished, (ii) the amount of genetic divergence driven by the two processes is similar, and (iii) comparisons can be repeated both taxonomically (for local adaptation) and geographically (for speciation). Here, we report just such a situation in the hamlets (Hypoplectrus spp), brightly colored reef fishes from the wider Caribbean. Close to 100,000 SNPs genotyped in 126 individuals from three sympatric species sampled in three repeated populations provide genome‐wide levels of divergence that are comparable among allopatric populations (F_st estimate = 0.0042) and sympatric species (F_st estimate = 0.0038). Population genetic, clustering, and phylogenetic analyses reveal very similar patterns for local adaptation and speciation, with a large fraction of the genome undifferentiated (F_st estimate ≈ 0), a very small proportion of F_st outlier loci (0.05–0.07%), and remarkably few repeated outliers (1–3). Nevertheless, different loci appear to be involved in the two processes in Hypoplectrus, with only 7% of the most differentiated SNPs and outliers shared between populations and species comparisons. In particular, a tropomyosin (Tpm4) and a previously identified hox (HoxCa) locus emerge as candidate loci (repeated outliers) for local adaptation and speciation, respectively. We conclude that marine populations may be locally adapted notwithstanding shallow levels of genetic divergence, and that from a population genomic perspective, this process does not appear to differ fundamentally from the early stages of speciation.     

Deciphering the demographic history of allochronic differentiation in the pine processionary moth Thaumetopoea pityocampa

Understanding the processes of adaptive divergence, which may ultimately lead to speciation, is a major question in evolutionary biology. Allochronic differentiation refers to a particular situation where gene flow is primarily impeded by temporal isolation between early and late reproducers. This process has been suggested to occur in a large array of organisms, even though it is still overlooked in the literature. We here focused on a well‐documented case of incipient allochronic speciation in the winter pine processionary moth Thaumetopoea pityocampa. This species typically reproduces in summer and larval development occurs throughout autumn and winter. A unique, phenologically shifted population (SP) was discovered in 1997 in Portugal. It was proved to be strongly differentiated from the sympatric “winter population” (WP), but its evolutionary history could only now be explored. We took advantage of the recent assembly of a draft genome and of the development of pan‐genomic RAD‐seq markers to decipher the demographic history of the differentiating populations and develop genome scans of adaptive differentiation. We showed that the SP diverged relatively recently, that is, few hundred years ago, and went through two successive bottlenecks followed by population size expansions, while the sympatric WP is currently experiencing a population decline. We identified outlier SNPs that were mapped onto the genome, but none were associated with the phenological shift or with subsequent adaptations. The strong genetic drift that occurred along the SP lineage certainly challenged our capacity to reveal functionally important loci.     

Visual adaptation could aid sympatric speciation in a deep crater lake

Allopatric speciation was originally suggested to be the primary mechanism of animal speciation (Mayr, 1942; Figure 1). During allopatric speciation, populations diverge when gene flow is reduced across significant biogeographic barriers. Sympatric speciation, where species diverge while inhabiting the same location, was thought to be essentially impossible. However, the advent of theoretical models followed by new experimental evidence made sympatric speciation more plausible (Via, 2001). The cichlid fishes of Barombi Mbo, a small crater lake in western Cameroon, became one of the most widely accepted examples of sympatric speciation (Schliewen, Tautz, & Paabo, 1994). Although the phylogenetic history of this clade is not quite as simple as originally thought, it remains one of the best examples of sympatric speciation (Richards, Poelstra, & Martin, 2018). However, little is known about the molecular mechanisms contributing to the splitting of these species in situ. In a From the Cover article in this issue of Molecular Ecology, Musilova et al. (2019) focus on the diversity of visual systems among these fishes. They identify genetic changes associated with several aspects of visual adaptation that may have contributed to the ecological specialization and sympatric speciation of cichlids in this lake.     

A test for within‐lake niche differentiation in the nine‐spined sticklebacks (Pungitius pungitius)

Specialization for the use of different resources can lead to ecological speciation. Accordingly, there are numerous examples of ecologically specialized pairs of fish “species” in postglacial lakes. Using a polymorphic panel of single nucleotide variants, we tested for genetic footprints of within‐lake population stratification in nine‐spined sticklebacks (Pungitius pungitius) collected from three habitats (viz. littoral, benthic, and pelagic) within a northern Swedish lake. Analyses of admixture, population structure, and relatedness all supported the conclusion that the fish from this lake form a single interbreeding unit.     

Crater lake cichlids individually specialize along the benthic–limnetic axis

A common pattern of adaptive diversification in freshwater fishes is the repeated evolution of elongated open water (limnetic) species and high‐bodied shore (benthic) species from generalist ancestors. Studies on phenotype‐diet correlations have suggested that population‐wide individual specialization occurs at an early evolutionary and ecological stage of divergence and niche partitioning. This variable restricted niche use across individuals can provide the raw material for earliest stages of sympatric divergence. We investigated variation in morphology and diet as well as their correlations along the benthic‐limnetic axis in an extremely young Midas cichlid species, Amphilophus tolteca, endemic to the Nicaraguan crater lake Asososca Managua. We found that A. tolteca varied continuously in ecologically relevant traits such as body shape and lower pharyngeal jaw morphology. The correlation of these phenotypes with niche suggested that individuals are specialized along the benthic‐limnetic axis. No genetic differentiation within the crater lake was detected based on genotypes from 13 microsatellite loci. Overall, we found that individual specialization in this young crater lake species encompasses the limnetic‐ as well as the benthic macro‐habitat. Yet there is no evidence for any diversification within the species, making this a candidate system for studying what might be the early stages preceding sympatric divergence.     

The youngest split in sympatric schizothoracine fish (Cyprinidae) is shaped by ecological adaptations in a Tibetan Plateau glacier lake

Although new empirical evidence shows that sympatric speciation has occurred in some species, there are few indisputable model organisms for this process of speciation. The two subspecies ( Gymnocypris eckloni eckloni and G . e . scoliostomus ) of the schizothoracine Gymnocypris fish species complex from a small glacier lake in the Tibetan Plateau, Lake Sunmcuo, fit several of the key characteristics of the sympatric speciation model. We used combined mitochondrial control region sequences and the cytochrome b gene (1894 bp) to address the phylogenetics and population genetics of 232 specimens of G . e . eckloni and G . e . scoliostomus , as well as all of its closely related sister species. We found that: (i) a total of four old lineages were uncovered in the widespread G . e . eckloni , of which only one was shown to be shared with all G . e . scoliostomus individuals and (ii) the new subspecies ( G . e . scoliostomus ) evolved in Lake Sunmcuo from the ancestral G . e . eckloni population within approximately 0.057 Ma. These two taxa of the species complex are morphologically distinct, and reproductive isolation is further suggested. Ecological disruptive selection based on morphological traits (e.g. mouth cleft characters) and food utilization may be a mechanism of incipient speciation of two sympatric populations within Lake Sunmcuo. This study provides the first genetic evidence for sympatric speciation in the schizothoracine fish.     

Sympatric speciation in a genus of marine reef fishes

Sympatric speciation has been contentious since its inception, yet is increasingly recognized as important based on accumulating theoretical and empirical support. Here, we present a compelling case of sympatric speciation in a taxon of marine reef fishes using a comparative and mechanistic approach. Hexagrammos otakii and H. agrammus occur in sympatry throughout their ranges. Molecular sequence data from six loci, with complete sampling of the genus, support monophyly of these sister species. Although hybridization occurrs frequently with an allopatric congener in an area of slight distributional overlap, we found no F₁ hybrids between the focal sympatric taxa throughout their coextensive ranges. We present genetic evidence for complete reproductive isolation based on SNP analysis of 382 individuals indicating fixed polymorphisms, with no shared haplotypes or genotypes, between sympatric species. To address questions of speciation, we take a mechanistic approach and directly compare aspects of reproductive isolation between allopatric and sympatric taxa both in nature and in the laboratory. We conclude that the buildup of reproductive isolation is strikingly different in sympatric vs. allopatric taxa, consistent with theoretical predictions. Lab reared hybrids from allopatric species crosses exhibit severe fitness effects in the F₁ or backcross generation. No intrinsic fitness effects are observed in F₁ hybrids from sympatric species pairs, however these treatments exhibited reduced fertilization success and complete pre‐mating isolation is implied in nature because F₁ hybrid adults do not occur. Our study addresses limitations of previous studies and supports new criteria for inferring sympatric speciation.     

Spatial and host related genomic variation in partially sympatric cactophagous moth species

Surveys of patterns of genetic variation in natural sympatric and allopatric populations of recently diverged species are necessary to understand the processes driving intra- and interspecific diversification. The South American moths Cactoblastis cactorum, Cactoblastis doddi and Cactoblastis bucyrus are specialized in the use of cacti as host plants. These species have partially different geographic ranges and differ in patterns of host plant use. However, there are areas that overlap, particularly, in northwestern Argentina, where they are sympatric. Using a combination of genome-wide SNPs and mitochondrial data we assessed intra and interspecific genetic variation and investigated the relative roles of geography and host plants on genetic divergence. We also searched for genetic footprints of hybridization between species. We identified three well delimited species and detected signs of hybridization in the area of sympatry. Our results supported a hypothetical scenario of allopatric speciation in the generalist C. cactorum and genetic interchange during secondary geographic contact with the pair of specialists C. bucyrus and C. doddi that probably speciated sympatrically. In both cases, adaptation to new host plants probably played an important role in speciation. The results also suggested the interplay of geography and host plant use as drivers of divergence and limiting gene flow at intra and interspecific levels.     

THE GENETIC ARCHITECTURE OF REPRODUCTIVE ISOLATION DURING SPECIATION‐WITH‐GENE‐FLOW IN LAKE WHITEFISH SPECIES PAIRS ASSESSED BY RAD SEQUENCING

During speciation‐with‐gene‐flow, effective migration varies across the genome as a function of several factors, including proximity of selected loci, recombination rate, strength of selection, and number of selected loci. Genome scans may provide better empirical understanding of the genome‐wide patterns of genetic differentiation, especially if the variance due to the previously mentioned factors is partitioned. In North American lake whitefish (Coregonus clupeaformis), glacial lineages that diverged in allopatry about 60,000 years ago and came into contact 12,000 years ago have independently evolved in several lakes into two sympatric species pairs (a normal benthic and a dwarf limnetic). Variable degrees of reproductive isolation between species pairs across lakes offer a continuum of genetic and phenotypic divergence associated with adaptation to distinct ecological niches. To disentangle the complex array of genetically based barriers that locally reduce the effective migration rate between whitefish species pairs, we compared genome‐wide patterns of divergence across five lakes distributed along this divergence continuum. Using restriction site associated DNA (RAD) sequencing, we combined genetic mapping and population genetics approaches to identify genomic regions resistant to introgression and derive empirical measures of the barrier strength as a function of recombination distance. We found that the size of the genomic islands of differentiation was influenced by the joint effects of linkage disequilibrium maintained by selection on many loci, the strength of ecological niche divergence, as well as demographic characteristics unique to each lake. Partial parallelism in divergent genomic regions likely reflected the combined effects of polygenic adaptation from standing variation and independent changes in the genetic architecture of postzygotic isolation. This study illustrates how integrating genetic mapping and population genomics of multiple sympatric species pairs provide a window on the speciation‐with‐gene‐flow mechanism.     

ADAPTIVE RADIATION AND ECOLOGICAL DIVERSIFICATION OF SULAWESI'S ANCIENT LAKE SHRIMPS

Ancient lakes are natural laboratories for the study of adaptive radiation. Recently, two lake systems on the Indonesian island Sulawesi have emerged as promising new model systems. A species flock of atyid freshwater shrimps in the Malili lake system comprises 15 colorful endemic taxa. Mitochondrial DNA data suggest two independent colonizations by riverine ancestors. Only one colonization event led to subsequent radiation into 14 species, while the second clade comprises just one species. The vast majority of species (n= 12) are habitat specialists, which are confined to the larger Malili clade and include all taxa with species‐specific color patterns and a restricted distribution within the five connected Malili lakes. Morphological, genetic, and ecological data are consistent with the existence of an adaptive radiation in the Malili lakes, involving the habitat‐specific diversification of trophic morphology. In addition to testing criteria for the recognition of an adaptive radiation, an ancestral state reconstruction reveals an equal probability for either a riverine generalist or a lacustrine specialist as ancestor of the large Malili clade, which is interpreted as indicative of an early stage of habitat specialization within this radiation. Finally, our results suggest that species diversification may have been primarily driven by ecological specialization and allopatric speciation.     

Pattern,process and geographic modes of speciation

The tradition of classifying cases of speciation into discrete geographic categories (allopatric, parapatric and sympatric) fuelled decades of fruitful research and debate. Not surprisingly, as the science has become more sophisticated, this simplistic taxonomy has become increasingly obsolete. Geographic patterns are now reasonably well understood. Sister species are rarely sympatric, implying that sympatric speciation, it its most general sense, is rare. However, sympatric speciation, even in its most restricted population genetic sense, is possible. Several case studies have demonstrated that divergence has occurred in nature without geographic barriers to gene flow. Obviously, different sets of criteria for sympatric speciation will lead to different numbers of qualifying cases. But changing the rules of nomenclature to make ‘sympatric speciation’ more or less common does not constitute scientific progress. Advances in the study of speciation have come from studies of the processes that constrain or promote divergence, and how they are affected by geography.     

Introgression and selection shaped the evolutionary history of sympatric sister‐species of coral reef fishes (genus: Haemulon)

Closely related marine species with large overlapping ranges provide opportunities to study mechanisms of speciation, particularly when there is evidence of gene flow between such lineages. Here, we focus on a case of hybridization between the sympatric sister‐species Haemulon maculicauda and H. flaviguttatum, using Sanger sequencing of mitochondrial and nuclear loci, as well as 2422 single nucleotide polymorphisms (SNPs) obtained via restriction site‐associated DNA sequencing (RADSeq). Mitochondrial markers revealed a shared haplotype for COI and low divergence for CytB and CR between the sister‐species. On the other hand, complete lineage sorting was observed at the nuclear loci and most of the SNPs. Under neutral expectations, the smaller effective population size of mtDNA should lead to fixation of mutations faster than nDNA. Thus, these results suggest that hybridization in the recent past (0.174–0.263 Ma) led to introgression of the mtDNA, with little effect on the nuclear genome. Analyses of the SNP data revealed 28 loci potentially under divergent selection between the two species. The combination of mtDNA introgression and limited nuclear DNA introgression provides a mechanism for the evolution of independent lineages despite recurrent hybridization events. This study adds to the growing body of research that exemplifies how genetic divergence can be maintained in the presence of gene flow between closely related species.     

Adaptive sympatric speciation of polychromatic "roundfin" sailfin silverside fish in Lake Matano (Sulawesi)

The significance of sympatric speciation is one of the most controversial topics in evolutionary biology. Theory suggests that different factors can lead to speciation in full geographical contact, including selection and nonrandom mating. Strict criteria have been established for assessing sympatric speciation, which have been met in only a very few cases. Here, we investigate differentiation among sympatric morphospecies and color morphs of "roundfin" sailfin silversides (Telmatherinidae), small freshwater fish endemic to ancient Lake Matano in Central Sulawesi (Indonesia). Morphospecies are distinct according to body shape (geometric morphometrics), population structure (population-level amplified fragment length polymorphism [AFLP] markers), ecology, and mating behavior (habitat transects, stomach contents). Explorative genome scans based on AFLPs indicate that divergent selection affects only 1.3-4.2% of the analyzed loci, suggesting an early stage of speciation. Transect data demonstrate strong assortative mating and adaptive niche differentiation. However, we find no restrictions in gene flow among the conspicuous male color morphs. In summary, our data are consistent with a sympatric mode of divergence among three morphospecies under conditions effectively ruling out allopatric scenarios. Substantial, but incomplete, reproductive isolation suggests an early stage of speciation, most likely due to ecological selection pressure.