Rapid sympatric ecological differentiation of crater lake cichlid fishes within historic times

Environmental correlates of adaptive diversification in postglacial freshwater fishes

Determining how environmental conditions contribute to divergence among populations and drive speciation is fundamental to resolving mechanisms and understanding outcomes in evolutionary biology. Postglacial freshwater fish species in the Northern Hemisphere are ideal biological systems to explore the effects of environment on diversification in morphology, ecology, and genetics (ecomorph divergences) within lakes. To date, various environmental factors have been implicated in the presence of multiple ecomorphs within particular lakes or regions. However, concerted evidence for generalizable patterns in environmental variables associated with speciation across geographical regions and across species and genera has been lacking. Here, we aimed to identify key biotic and abiotic factors associated with ecological divergence of postglacial freshwater fish species into multiple sympatric ecomorphs, focusing on species in the well-studied, widespread, and co-distributed genera Gasterosteus, Salvelinus, and Coregonus (stickleback, charr, and whitefish, respectively). We found that the presence of multiple sympatric ecomorphs tended to be associated with increasing lake surface area, maximum depth, and nutrient availability. In addition, predation, competition, and prey availability were suggested to play a role in divergence into multiple ecomorphs, but the effects of biotic factors require further study. Although we identified several environmental factors correlated with the presence of multiple ecomorphs, there were substantial data gaps across species and regions. An improved understanding of these systems may provide insight into both generalizable environmental factors involved in speciation in other systems, and potential ecological and evolutionary responses of species complexes when these variables are altered by environmental change.

Phenotypic rate and state are decoupled in response to river-to-lake transitions in cichlid fishes

Geographic access to isolated ecosystems is an important catalyst of adaptive radiation. Cichlid fishes repeatedly colonized rift, crater, and volcanic lakes from surrounding rivers. We test the "lake effect" on the phenotypic rate and state across 253 cichlid species. The rate of evolution was consistently higher (~10-fold) in lakes, and consistent across different dimensions of the phenotype. Rate shifts tended to occur coincident with or immediately following river-to-lake transitions, generally resulting in 2- to 5-fold faster rates than in the founding riverine lineage. By contrast, river- and lake-dwelling cichlids exhibit considerable overlap in phenotypes, generally with less disparity in lakes, but often different evolutionary optima. Taken together, these results suggest that lake radiations rapidly expand into niches largely already represented by ancestral riverine lineages, albeit in different frequencies. Lakes may provide ecological opportunity via ecological release (e.g., from predators/competitors) but need not be coupled with access to novel ecological niches.

Early stages of sympatric homoploid hybrid speciation in crater lake cichlid fishes

Homoploid hybrid speciation (i.e., hybrid speciation without a change in ploidy) has traditionally been considered to be rare in animals. Only few accepted empirical examples of homoploid hybrid speciation in nature exist, and in only one previous case (insects) was it convincingly shown that this process occurred in complete sympatry. Here, we report an instance of sympatric homoploid hybrid speciation in Midas cichlid fishes in Crater Lake Xiloá, Nicaragua. The hybrid lineage, albeit at an early stage of speciation, has genomically and phenotypically diverged from both of its two parental species. Together with a distinct stable isotope signature this suggests that this hybrid lineages occupies a different trophic niche compared to the other sympatric Midas cichlid species in Crater Lake Xiloá.

Here, the authors present a rare example of sympatric homoploid hybrid speciation, without ploidy changes, in the Midas cichlid fishes from Nicaragua. Midas cichlid hybrids occupy a different ecological niche, likely facilitated by body shape adaptations.

Divergent and non-parallel evolution of MHC IIB in the Neotropical Midas cichlid species complex

Background

Ecological diversification is the result of divergent natural selection by contrasting habitat characteristics that favours the evolution of distinct phenotypes. This process can happen in sympatry and in allopatry. Habitat-specific parasite communities have the potential to drive diversification among host populations by imposing selective pressures on their host's immune system. In particular, the hyperdiverse genes of the major histocompatibility complex (MHC) are implicated in parasite-mediated host divergence. Here, we studied the extent of divergence at MHC, and discuss how it may have contributed to the Nicaraguan Midas cichlid species complex diversification, one of the most convincing examples of rapid sympatric parallel speciation.

Results

We genotyped the MHC IIB for individuals from six sympatric Midas cichlid assemblages, each containing species that have adapted to exploit similar habitats. We recovered large allelic and functional diversity within the species complex. While most alleles were rare, functional groups of alleles (supertypes) were common, suggesting that they are key to survival and that they were maintained during colonization and subsequent radiations. We identified lake-specific and habitat-specific signatures for both allelic and functional diversity, but no clear pattern of parallel divergence among ecomorphologically similar phenotypes.

Conclusions

Colonization and demographic effects of the fish could have contributed to MHC evolution in the Midas cichlid in conjunction with habitat-specific selective pressures, such as parasites associated to alternative preys or environmental features. Additional ecological data will help evaluating the role of host–parasite interactions in the Midas cichlid radiations and aid in elucidating the potential role of non-parallel features differentiating crater lake species assemblages.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-01997-9.

The macroparasite fauna of cichlid fish from Nicaraguan lakes, a model system for understanding host-parasite diversification and speciation

The Nicaraguan lakes represent an ideal continent-island-like setting to study the colonization patterns of both fish and their parasites. The dominant fish fauna are cichlids, particularly the Midas cichlid species complex Amphilophus spp., a well-studied model for recent sympatric speciation. Here, we characterized the Midas cichlid macroparasite diversity in Nicaraguan lakes. We evaluated patterns of parasite diversity across host populations. Morphological and molecular analyses were conducted, revealing a macroparasite fauna composed by 37 taxa, including platyhelminths, nematodes, copepods, branchiurans, hirudineans and oribatids. Three invasive species are reported for the first time. The Midas cichlid was infected by 22 parasite taxa, 18 shared with other cichlids. Eight taxa conformed the core parasite fauna of the Midas cichlid. The large lakes had higher parasite diversity than the smaller and isolated crater lakes, although parasite infracommunity diversity was lower. Environmental factors along with the differential distribution of intermediate hosts, the potential resistance gained by their hosts after colonization of new lakes, competitive exclusion among parasites, or the introduction of exotic fish, may determine the observed pattern of parasite heterogeneous distribution. Our study provides a ground to explore the evolutionary history of both, hosts and parasites within the context of speciation and diversification processes.

Rapid Parallel Morphological and Mechanical Diversification of South American Pike Cichlids (Crenicichla)

Explosive bouts of diversification are one of the most conspicuous features of the tree of life. When such bursts are repeated in similar environments it suggests some degree of predictability in the evolutionary process. We assess parallel adaptive radiation of South American pike cichlids (Crenicichla) using phylogenomics and phylogenetic comparative methods. We find that species flocks in the Uruguay and Iguazú River basins rapidly diversified into the same set of ecomorphs that reflect feeding ecology. Both adaptive radiations involve expansion of functional morphology, resulting in unique jaw phenotypes. Yet, form and function were decoupled such that most ecomorphs share similar mechanical properties of the jaws (i.e., jaw motion during a feeding strike). Prey mobility explained six to nine-fold differences in the rate of morphological evolution, but had no effect on the rate of mechanical evolution. We find no evidence of gene flow between species flocks or with surrounding coastal lineages that may explain their rapid diversification. When compared to cichlids of the East African Great Lakes and other prominent adaptive radiations, pike cichlids share many themes, including rapid expansion of phenotypic diversity, specialization along the benthic-to-pelagic habitat and soft-to-hard prey axes, and the evolution of conspicuous functional innovations. Yet, decoupled evolution of form and function and the absence of hybridization as a catalyzing force are departures from patterns observed in other adaptive radiations. Many-to-one mapping of morphology to mechanical properties is a mechanism by which pike cichlids exhibit a diversity of feeding ecologies while avoiding exacerbating underlying mechanical trade-offs.

Phylogeography of high Andean killifishes Orestias (Teleostei: Cyprinodontidae) in Caquena and Lauca sub-basins of the Altiplano (Chile): mitochondrial and nuclear analysis of an endangered fish

From the early Miocene, the uplift of the Andes Mountains, intense volcanic activity and the occurrence of successive periods of dryness and humidity would have differentially influenced the modification of Altiplano watersheds, and consequently the evolutionary history of the taxa that live there. We analyzed Orestias populations from the Caquena and Lauca Altiplanic sub-basins of northern Chile to determine their genetic differentiation and relationship to their geographical distribution using mitochondrial (D-loop) and nuclear (microsatellite) molecular markers and to reconstruct its biogeographic history on these sub-basins. The results allowed reconstructing and reevaluating the evolutionary history of the genus in the area; genic diversity and differentiation together with different founding genetic groups suggest that Orestias have been spread homogeneously in the study area and would have experienced local disturbances that promoted isolation and diversification in restricted zones of their distribution.

Pervasive admixture and the spread of a large-lipped form in a cichlid fish radiation

Adaptive radiations have proven important for understanding the mechanisms and processes underlying biological diversity. The convergence of form and function, as well as admixture and adaptive introgression, are common in adaptive radiations. However, distinguishing between these two scenarios remains a challenge for evolutionary research. The Midas cichlid species complex (Amphilophus spp.) is a prime example of adaptive radiation, with phenotypic diversification occurring at various stages of genetic differentiation. One species, A. labiatus, has large fleshy lips, is associated with rocky lake substrates, and occurs patchily within Lakes Nicaragua and Managua. By contrast, the similar, but thin-lipped, congener, A. citrinellus, is more common and widespread. We investigated the evolutionary history of the large-lipped form, specifically regarding whether the trait has evolved independently in both lakes from ancestral thin-lipped populations, or via dispersal and/or admixture events. We collected samples from distinct locations in both lakes, and assessed differences in morphology and ecology. Using RAD-seq, we genotyped thousands of SNPs to measure population structure and divergence, demographic history, and admixture. We found significant between-species differences in ecology and morphology, local intraspecific differences in body shape and trophic traits, but only limited intraspecific variation in lip shape. Despite clear ecological differences, our genomic approach uncovered pervasive admixture between the species and low genomic differentiation, with species within lakes being genetically more similar than species between lakes. Taken together, our results suggest a single origin of large-lips, followed by pervasive admixture and adaptive introgression, with morphology being driven by local ecological opportunities, despite ongoing gene-flow.

The comparative genomic landscape of adaptive radiation in crater lake cichlid fishes

Factors ranging from ecological opportunity to genome composition might explain why only some lineages form adaptive radiations. While being rare, particular systems can provide natural experiments within an identical ecological setting where species numbers and phenotypic divergence in two closely related lineages are notably different. We investigated one such natural experiment using two de novo assembled and 40 resequenced genomes and asked why two closely related Neotropical cichlid fish lineages, the Amphilophus citrinellus species complex (Midas cichlids; radiating) and Archocentrus centrarchus (Flyer cichlid; nonradiating), have resulted in such disparate evolutionary outcomes. Although both lineages inhabit many of the same Nicaraguan lakes, whole-genome inferred demography suggests that priority effects are not likely to be the cause of the dissimilarities. Also, genome-wide levels of selection, transposable element dynamics, gene family expansion, major chromosomal rearrangements and the number of genes under positive selection were not markedly different between the two lineages. To more finely investigate particular subsets of the genome that have undergone adaptive divergence in Midas cichlids, we also examined if there was evidence for 'molecular pre-adaptation' in regions identified by QTL mapping of repeatedly diverging adaptive traits. Although most of our analyses failed to pinpoint substantial genomic differences, we did identify functional categories containing many genes under positive selection that provide candidates for future studies on the propensity of Midas cichlids to radiate. Our results point to a disproportionate role of local, rather than genome-wide factors underlying the propensity for these cichlid fishes to adaptively radiate.

Intraspecific variation and structuring of phenotype in a lake-dwelling species are driven by lake size and elevation

The fragmented, heterogeneous and relatively depauperate ecosystems of recently glaciated lakes present contrasting ecological opportunities for resident fish. Across a species, local adaptation may induce diverse and distinct phenotypic responses to various selection pressures. We tested for intraspecific phenotypic structuring by population in a common native lake-dwelling fish species across a medium-scale geographic region with considerable variation in lake types. We investigated potential lake-characteristic drivers of trophic morphology. Using geometric morphometric techniques, we quantified the head shapes of 759 adult brown trout (Salmo trutta L.) from 28 lakes and reservoirs across Scotland. Multivariate statistical analyses showed that almost all populations differed from one another. Trout from larger and deeper lakes had deeper, but shorter heads, and smaller eyes. Higher elevation lakes were associated with fish with shorter heads and jaws. Within-population shape variation also differed by population, and was positively correlated with lake surface area and depth. Trout within reservoirs differed subtly from those in natural lakes, having larger eyes, shorter jaws and greater variability. This study documents an extraordinary morphological variation between and within populations of brown trout, and demonstrates the role of the extrinsic environment in driving phenotypic structuring over a medium-scale and varied geographic area.

Ecological speciation in European whitefish is driven by a large-gaped predator

Lake‐dwelling fish that form species pairs/flocks characterized by body size divergence are important model systems for speciation research. Although several sources of divergent selection have been identified in these systems, their importance for driving the speciation process remains elusive. A major problem is that in retrospect, we cannot distinguish selection pressures that initiated divergence from those acting later in the process. To address this issue, we studied the initial stages of speciation in European whitefish (Coregonus lavaretus) using data from 358 populations of varying age (26–10,000 years). We find that whitefish speciation is driven by a large‐growing predator, the northern pike (Esox lucius). Pike initiates divergence by causing a largely plastic differentiation into benthic giants and pelagic dwarfs: ecotypes that will subsequently develop partial reproductive isolation and heritable differences in gill raker number. Using an eco‐evolutionary model, we demonstrate how pike's habitat specificity and large gape size are critical for imposing a between‐habitat trade‐off, causing prey to mature in a safer place or at a safer size. Thereby, we propose a novel mechanism for how predators may cause dwarf/giant speciation in lake‐dwelling fish species.

Molecular clocks, biogeography and species diversity in Herichthys with evaluation of the role of Punta del Morro as a vicariant brake along the Mexican Transition Zone in the context of local and global time frame of cichlid diversification

Using molecular dated phylogenies and biogeographic reconstructions, the species diversity, biogeography and time frame of evolution of the genus Herichthys were evaluated. In particular, we test the role of Punta del Morro (PdM) as a vicariant brake along the Mexican Transition Zone in the context of local and global time frame of cichlid diversification using several sets of calibrations. Species diversity in Herichthys is complex and the here employed dating methods suggest young age and rapid divergence for many species while species delimitation methods did not resolve these young species including both sympatric species pairs. Based on our molecular clock dating analyses, Herichthys has colonized its present distribution area significantly prior to the suggested vicariance by PdM (10-17.1 Ma vs. 5 to 7.5 Ma). The PdM constraint is in conflict with all other paleogeographic and fossil constraints including novel ones introduced in this study that are, however, congruent among each other. Our study demonstrates that any cichlid datings significantly older or younger than the bounds presented by our analyses and discussion have to be taken as highly questionable from the point of view of Middle American paleogeography and cichlid biogeography unless we allow the option that cichlid biogeography is completely independent from ecological and geological constraints.

Parallelism in eco-morphology and gene expression despite variable evolutionary and genomic backgrounds in a Holarctic fish

Understanding the extent to which ecological divergence is repeatable is essential for predicting responses of biodiversity to environmental change. Here we test the predictability of evolution, from genotype to phenotype, by studying parallel evolution in a salmonid fish, Arctic charr (Salvelinus alpinus), across eleven replicate sympatric ecotype pairs (benthivorous-planktivorous and planktivorous-piscivorous) and two evolutionary lineages. We found considerable variability in eco-morphological divergence, with several traits related to foraging (eye diameter, pectoral fin length) being highly parallel even across lineages. This suggests repeated and predictable adaptation to environment. Consistent with ancestral genetic variation, hundreds of loci were associated with ecotype divergence within lineages of which eight were shared across lineages. This shared genetic variation was maintained despite variation in evolutionary histories, ranging from postglacial divergence in sympatry (ca. 10-15kya) to pre-glacial divergence (ca. 20-40kya) with postglacial secondary contact. Transcriptome-wide gene expression (44,102 genes) was highly parallel across replicates, involved biological processes characteristic of ecotype morphology and physiology, and revealed parallelism at the level of regulatory networks. This expression divergence was not only plastic but in part genetically controlled by parallel cis-eQTL. Lastly, we found that the magnitude of phenotypic divergence was largely correlated with the genetic differentiation and gene expression divergence. In contrast, the direction of phenotypic change was mostly determined by the interplay of adaptive genetic variation, gene expression, and ecosystem size. Ecosystem size further explained variation in putatively adaptive, ecotype-associated genomic patterns within and across lineages, highlighting the role of environmental variation and stochasticity in parallel evolution. Together, our findings demonstrate the parallel evolution of eco-morphology and gene expression within and across evolutionary lineages, which is controlled by the interplay of environmental stochasticity and evolutionary contingencies, largely overcoming variable evolutionary histories and genomic backgrounds.

Implementing Large Genomic Single Nucleotide Polymorphism Data Sets in Phylogenetic Network Reconstructions: A Case Study of Particularly Rapid Radiations of Cichlid Fish

The Midas cichlids of the Amphilophus citrinellus spp. species complex from Nicaragua (13 species) are an extraordinary example of adaptive and rapid radiation ($<$24,000 years old). These cichlids are a very challenging group to infer its evolutionary history in phylogenetic analyses, due to the apparent prevalence of incomplete lineage sorting (ILS), as well as past and current gene flow. Assuming solely a vertical transfer of genetic material from an ancestral lineage to new lineages is not appropriate in many cases of genes transferred horizontally in nature. Recently developed methods to infer phylogenetic networks under such circumstances might be able to circumvent these problems. These models accommodate not just ILS, but also gene flow, under the multispecies network coalescent (MSNC) model, processes that are at work in young, hybridizing, and/or rapidly diversifying lineages. There are currently only a few programs available that implement MSNC for estimating phylogenetic networks. Here, we present a novel way to incorporate single nucleotide polymorphism (SNP) data into the currently available PhyloNetworks program. Based on simulations, we demonstrate that SNPs can provide enough power to recover the true phylogenetic network. We also show that it can accurately infer the true network more often than other similar SNP-based programs (PhyloNet and HyDe). Moreover, our approach results in a faster algorithm compared to the original pipeline in PhyloNetworks, without losing power. We also applied our new approach to infer the phylogenetic network of Midas cichlid radiation. We implemented the most comprehensive genomic data set to date (RADseq data set of 679 individuals and $>$37K SNPs from 19 ingroup lineages) and present estimated phylogenetic networks for this extremely young and fast-evolving radiation of cichlid fish. We demonstrate that the MSNC is more appropriate than the multispecies coalescent alone for the analysis of this rapid radiation. [Genomics; multispecies network coalescent; phylogenetic networks; phylogenomics; RADseq; SNPs.]

MicroRNA Gene Regulation in Extremely Young and Parallel Adaptive Radiations of Crater Lake Cichlid Fish

Cichlid fishes provide textbook examples of explosive phenotypic diversification and sympatric speciation, thereby making them ideal systems for studying the molecular mechanisms underlying rapid lineage divergence. Despite the fact that gene regulation provides a critical link between diversification in gene function and speciation, many genomic regulatory mechanisms such as microRNAs (miRNAs) have received little attention in these rapidly diversifying groups. Therefore, we investigated the posttranscriptional regulatory role of miRNAs in the repeated sympatric divergence of Midas cichlids (Amphilophus spp.) from Nicaraguan crater lakes. Using miRNA and mRNA sequencing of embryos from five Midas species, we first identified miRNA binding sites in mRNAs and highlighted the presences of a surprising number of novel miRNAs in these adaptively radiating species. Then, through analyses of expression levels, we identified putative miRNA/gene target pairs with negatively correlated expression level that were consistent with the role of miRNA in downregulating mRNA. Furthermore, we determined that several miRNA/gene pairs show convergent expression patterns associated with the repeated benthic/limnetic sympatric species divergence implicating these miRNAs as potential molecular mechanisms underlying replicated sympatric divergence. Finally, as these candidate miRNA/gene pairs may play a central role in phenotypic diversification in these cichlids, we characterized the expression domains of selected miRNAs and their target genes via in situ hybridization, providing further evidence that miRNA regulation likely plays a role in the Midas cichlid adaptive radiation. These results provide support for the hypothesis that extremely quickly evolving miRNA regulation can contribute to rapid evolutionary divergence even in the presence of gene flow.

Searching for Sympatric Speciation in the Genomic Era

Sympatric speciation illustrates how natural and sexual selection may create new species in isolation without geographic barriers. However, recent genomic reanalyses of classic examples of sympatric speciation reveal complex histories of secondary gene flow from outgroups into the radiation. In contrast, the rich theoretical literature on this process distinguishes among a diverse range of models based on simple genetic histories and different types of reproductive isolating barriers. Thus, there is a need to revisit how to connect theoretical models of sympatric speciation and their predictions to empirical case studies in the face of widespread gene flow. Here, theoretical differences among different types of sympatric speciation and speciation-with-gene-flow models are reviewed and summarized, and genomic analyses are proposed for distinguishing which models apply to case studies based on the timing and function of adaptive introgression. Investigating whether secondary gene flow contributed to reproductive isolation is necessary to test whether predictions of theory are ultimately borne out in nature.

The onset of ecological diversification 50 years after colonization of a crater lake by haplochromine cichlid fishes

Adaptive radiation research typically relies on the study of evolution in retrospective, leaving the predictive value of the concept hard to evaluate. Several radiations, including the cichlid fishes in the East African Great Lakes, have been studied extensively, yet no study has investigated the onset of the intraspecific processes of niche expansion and differentiation shortly after colonization of an adaptive zone by cichlids. Haplochromine cichlids of one of the two lineages that seeded the Lake Victoria radiation recently arrived in Lake Chala, a lake perfectly suited for within-lake cichlid speciation. Here, we infer the colonization and demographic history, quantify phenotypic, ecological and genomic diversity and diversification, and investigate the selection regime to ask if the population shows signs of diversification resembling the onset of adaptive radiation. We find that since their arrival in the lake, haplochromines have colonized a wide range of depth habitats associated with ecological and morphological expansion and the beginning of phenotypic differentiation and potentially nascent speciation, consistent with the very early onset of an adaptive radiation process. Moreover, we demonstrate evidence of rugged phenotypic fitness surfaces, indicating that current ecological selection may contribute to the phenotypic diversification.

Convergent phenotypic evolution of the visual system via different molecular routes: How Neotropical cichlid fishes adapt to novel light environments

How predictable is evolution? This remains a fundamental but contested issue in evolutionary biology. When independent lineages colonize the same environment, we are presented with a natural experiment that allows us to ask if genetic and ecological differences promote species-specific evolutionary outcomes or whether species phenotypically evolve in a convergent manner in response to shared selection pressures. If so, are the molecular mechanisms underlying phenotypic convergence the same? In Nicaragua, seven species of cichlid fishes concurrently colonized two novel photic environments. Hence, their visual system represents a compelling model to address these questions, particularly since the adaptive value of phenotypic changes is well-understood. By analyzing retinal transcriptomes, we found that differential expression of genes responsible for color vision (cone opsins and cyp27c1) produced rapid and mostly convergent changes of predicted visual sensitivities. Notably, these changes occurred in the same direction in all species although there were differences in underlying gene expression patterns illustrating nonconvergence at the molecular level. Adaptive phenotypes evolved deterministically, even when species differ substantially in ecology and genetic variation. This provides strong evidence that phenotypic evolution of the visual system occurred in response to similar selective forces of the photic environment.

Lessons from a natural experiment: Allopatric morphological divergence and sympatric diversification in the Midas cichlid species complex are largely influenced by ecology in a deterministic way

Explaining why some lineages diversify while others do not and how are key objectives in evolutionary biology. Young radiations of closely related species derived from the same source population provide an excellent opportunity to disentangle the relative contributions of possible drivers of diversification. In these settings, lineage-specific effects are shared and can be ruled out. Moreover, the relevant demographic and ecological parameters can be estimated accurately. Midas cichlid fish in Nicaragua have repeatedly colonized several crater lakes, diverged from the same source populations, and, interestingly, diversified in some of them but not others. Here, using the most comprehensive molecular and geometric morphometric data set on Midas cichlids to date (∼20,000 SNPs, 12 landmarks, ∼700 individuals), we aim to understand why and how crater lake populations diverge and why some of them are more prone to diversify in sympatry than others. Taking ancestor-descendant relationships into account, we find that Midas cichlids diverged in parallel from their source population mostly-but not exclusively-by evolving more slender body shapes in all six investigated crater lakes. Admixture among crater lakes has possibly facilitated this process in one case, but overall, admixture and secondary waves of colonization cannot predict morphological divergence and intralacustrine diversification. Instead, morphological divergence is larger the more dissimilar a crater lake is compared to the source lake and happens rapidly after colonization followed by a slow-down with time. Our data also provide some evidence that founder effects may positively contribute to divergence. The depth of a crater lake is positively associated with variation in body shapes (and number of species), presumably by providing more ecological opportunities. In conclusion, we find that parallel morphological divergence in allopatry and the propensity for diversification in sympatry across the entire Midas cichlid fish radiation is partly predictable and mostly driven by ecology.

Gene(s) and individual feeding behavior: Exploring eco-evolutionary dynamics underlying left-right asymmetry in the scale-eating cichlid fish Perissodus microlepis

The scale‐eating cichlid fish Perissodus microlepis is a textbook example of bilateral asymmetry due to its left or right‐bending heads and of negative frequency‐dependent selection, which is proposed to maintain this stable polymorphism. The mechanisms that underlie this asymmetry remain elusive. Several studies had initially postulated a simple genetic basis for this trait, but this explanation has been questioned, particularly by reports observing a unimodal distribution of mouth shapes. We hypothesize that this unimodal distribution might be due to a combination of genetic and phenotypically plastic components. Here, we expanded on previous work by investigating a formerly identified candidate SNP associated to mouth laterality, documenting inter‐individual variation in feeding preference using stable isotope analyses, and testing their association with mouth asymmetry. Our results suggest that this polymorphism is influenced by both a polygenic basis and inter‐individual non‐genetic variation, possibly due to feeding experience, individual specialization, and intraspecific competition. We introduce a hypothesis potentially explaining the simultaneous maintenance of left, right, asymmetric and symmetric mouth phenotypes due to the interaction between diverse eco‐evolutionary dynamics including niche construction and balancing selection. Future studies will have to further tease apart the relative contribution of genetic and environmental factors and their interactions in an integrated fashion.

Sympatric Speciation in the Genomic Era

Sympatric speciation has been of key interest to biologists investigating how natural and sexual selection drive speciation without the confounding variable of geographic isolation. The advent of the genomic era has provided a more nuanced and quantitative understanding of the different and often complex modes of speciation by which sympatric sister taxa arose, and a reassessment of some of the most compelling empirical case studies of sympatric speciation. However, I argue that genomic studies based on contemporary populations may never be able to provide unequivocal evidence of true primary sympatric speciation, and there is a need to incorporate palaeogenomic studies into this field. This inability to robustly distinguish cases of primary and secondary 'divergence with gene flow' may be inconsequential, as both are useful for understanding the role of large effect barrier loci in the progression from localised genic isolation to genome-wide reproductive isolation. I argue that they can be of equivalent interest due to shared underlying mechanisms driving divergence and potentially leaving similar patterns of coalescence.

Accelerated Evolution and Functional Divergence of the Dim Light Visual Pigment Accompanies Cichlid Colonization of Central America

Cichlids encompass one of the most diverse groups of fishes in South and Central America, and show extensive variation in life history, morphology, and colouration. While studies of visual system evolution in cichlids have focussed largely on the African rift lake species flocks, Neotropical cichlids offer a unique opportunity to investigate visual system evolution at broader temporal and geographic scales. South American cichlid colonization of Central America has likely promoted accelerated rates of morphological evolution in Central American lineages as they encountered reduced competition, renewed ecological opportunity, and novel aquatic habitats. To investigate whether such transitions have influenced molecular evolution of vision in Central American cichlids, we sequenced the dim-light rhodopsin gene in 101 Neotropical cichlid species, spanning the diversity of the clade. We find strong evidence for increased rates of evolution in Central American cichlid rhodopsin relative to South American lineages, and identify several sites under positive selection in rhodopsin that likely contribute to adaptation to different photic environments. We expressed a Neotropical cichlid rhodopsin protein invitro for the first time, and found that while its spectral tuning properties were characteristic of typical vertebrate rhodopsin pigments, the rate of decay of its active signalling form was much slower, consistent with dim light adaptation in other vertebrate rhodopsins. Using site-directed mutagenesis combined with spectroscopic assays, we found that a key amino acid substitution present in some Central American cichlids accelerates the rate of decay of active rhodopsin, which may mediate adaptation to clear water habitats.

Phylogenomics of pike cichlids (Cichlidae: Crenicichla): the rapid ecological speciation of an incipient species flock

The rapid rise of phenotypic and ecological diversity in independent lake-dwelling groups of cichlids is emblematic of the East African Great Lakes. In this study, we show that similar ecologically based diversification has occurred in pike cichlids (Crenicichla) throughout the Uruguay River drainage of South America. We collected genomic data from nearly 500 ultraconserved element (UCEs) loci and >260 000 base pairs across 33 species, to obtain a phylogenetic hypothesis for the major species groups and to evaluate the relationships and genetic structure among five closely related, endemic, co-occurring species (the Uruguay River species flock; URSF). Additionally, we evaluated ecological divergence of the URSF based on body and lower pharyngeal jaw (LPJ) shape and gut contents. Across the genus, we recovered novel relationships among the species groups. We found strong support for the monophyly of the URSF; however, relationships among these species remain problematic, likely because of the rapid and recent evolution of this clade. Clustered co-ancestry analysis recovered most species as well delimited genetic groups. The URSF species exhibit species-specific body and LPJ shapes associated with specialized trophic roles. Collectively, our results suggest that the URSF consists of incipient species that arose via ecological speciation associated with the exploration of novel trophic roles.

Ecosystem size predicts eco-morphological variability in a postglacial diversification

Identifying the processes by which new phenotypes and species emerge has been a long‐standing effort in evolutionary biology. Young adaptive radiations provide a model to study patterns of morphological and ecological diversification in environmental context. Here, we use the recent radiation (ca. 12k years old) of the freshwater fish Arctic charr (Salvelinus alpinus) to identify abiotic and biotic environmental factors associated with adaptive morphological variation. Arctic charr are exceptionally diverse, and in postglacial lakes there is strong evidence of repeated parallel evolution of similar morphologies associated with foraging. We measured head depth (a trait reflecting general eco‐morphology and foraging ecology) of 1,091 individuals across 30 lake populations to test whether fish morphological variation was associated with lake bathymetry and/or ecological parameters. Across populations, we found a significant relationship between the variation in head depth of the charr and abiotic environmental characteristics: positively with ecosystem size (i.e., lake volume, surface area, depth) and negatively with the amount of littoral zone. In addition, extremely robust‐headed phenotypes tended to be associated with larger and deeper lakes. We identified no influence of co‐existing biotic community on Arctic charr trophic morphology. This study evidences the role of the extrinsic environment as a facilitator of rapid eco‐morphological diversification.

Incipient speciation driven by hypertrophied lips in Midas cichlid fishes?

Sympatric speciation has been debated in evolutionary biology for decades. Although it has gained in acceptance recently, still only a handful of empirical examples are seen as valid (e.g. crater lake cichlids). In this study, we disentangle the role of hypertrophied lips in the repeated adaptive radiations of Nicaraguan crater lake cichlid fish. We assessed the role of disruptive selection and assortative mating during the early stages of divergence and found a functional trade-off in feeding behaviour between thick- and thin-lipped ecotypes, suggesting that this trait is a target of disruptive selection. Thick-lipped fish perform better on nonevasive prey at the cost of a poorer performance on evasive prey. Using enclosures in the wild, we found that thick-lipped fish perform significantly better in rocky than in sandy habitats. We found almost no mixed pairs during two breeding seasons and hence significant assortative mating. Genetic differentiation between ecotypes seems to be related to the time since colonization, being subtle in L. Masaya (1600 generations ago) and absent in the younger L. Apoyeque (<600 generations ago). Genome-wide differentiation between ecotypes was higher in the old source lakes than in the young crater lakes. Our results suggest that hypertrophied lips might be promoting incipient sympatric speciation through divergent selection (ecological divergence in feeding performance) and nonrandom mating (assortative mating) in the young Nicaraguan crater lakes. Nonetheless, further manipulative experiments are needed in order to confirm the role of hypertrophied lips as the main cue for assortative mating.

Molecular characterization of MHC class IIB genes of sympatric Neotropical cichlids

BACKGROUND

The Major Histocompatibility Complex (MHC) is a key component of the adaptive immune system of all vertebrates and consists of the most polymorphic genes known to date. Due to this complexity, however, MHC remains to be characterized in many species including any Neotropical cichlid fish. Neotropical crater lake cichlids are ideal models to study evolutionary processes as they display one of the most convincing examples of sympatric and repeated parallel radiation events within and among isolated crater lakes.

RESULTS

Here, we characterized the genes of MHC class IIB chain of the Midas cichlid species complex (Amphilophus cf. citrinellus) including fish from five lakes in Nicaragua. We designed 19 new specific primers anchored in a stepwise fashion in order to detect all alleles present. We obtained 866 genomic DNA (gDNA) sequences from thirteen individuals and 756 additional sequences from complementary DNA (cDNA) of seven of those individuals. We identified 69 distinct alleles with up to 25 alleles per individual. We also found considerable intron length variation and mismatches of alleles detected in cDNA and gDNA suggesting that some loci have undergone pseudogenization. Lastly, we created a model of protein structure homology for each allele and identified their key structural components.

CONCLUSIONS

Overall, the Midas cichlid has one of the most diverse repertoires of MHC class IIB genes known, which could serve as a powerful tool to elucidate the process of divergent radiations, colonization and speciation in sympatry.

Niche divergence facilitated by fine-scale ecological partitioning in a recent cichlid fish adaptive radiation

Ecomorphological differentiation is a key feature of adaptive radiations, with a general trend for specialization and niche expansion following divergence. Ecological opportunity afforded by invasion of a new habitat is thought to act as an ecological release, facilitating divergence, and speciation. Here, we investigate trophic adaptive morphology and ecology of an endemic clade of oreochromine cichlid fishes (Alcolapia) that radiated along a herbivorous trophic axis following colonization of an isolated lacustrine environment, and demonstrate phenotype‐environment correlation. Ecological and morphological divergence of the Alcolapia species flock are examined in a phylogenomic context, to infer ecological niche occupation within the radiation. Species divergence is observed in both ecology and morphology, supporting the importance of ecological speciation within the radiation. Comparison with an outgroup taxon reveals large‐scale ecomorphological divergence but shallow genomic differentiation within the Alcolapia adaptive radiation. Ancestral morphological reconstruction suggests lake colonization by a generalist oreochromine phenotype that diverged in Lake Natron to varied herbivorous morphologies akin to specialist herbivores in Lakes Tanganyika and Malawi.

Continental cichlid radiations: functional diversity reveals the role of changing ecological opportunity in the Neotropics

Adaptive radiations have been hypothesized to contribute broadly to the diversity of organisms. Models of adaptive radiation predict that ecological opportunity and ecological release, the availability of empty ecological niches and the response by adapting lineages to occupy them, respectively, drive patterns of phenotypic and lineage diversification. Adaptive radiations driven by 'ecological opportunity' are well established in island systems; it is less clear if ecological opportunity influences continent-wide diversification. We use Neotropical cichlid fishes to test if variation in rates of functional evolution is consistent with changing ecological opportunity. Across a functional morphological axis associated with ram-suction feeding traits, evolutionary rates declined through time as lineages diversified in South America. Evolutionary rates of ram-suction functional morphology also appear to have accelerated as cichlids colonized Central America and encountered renewed opportunity. Our results suggest that ecological opportunity may play an important role in shaping patterns of morphological diversity of even broadly distributed lineages like Neotropical cichlids.

Correlated evolution of personality, morphology and performance

Evolutionary change in one trait can elicit evolutionary changes in other traits due to genetic correlations. This constrains the independent evolution of traits and can lead to unpredicted ecological and evolutionary outcomes. Animals might frequently exhibit genetic associations among behavioural and morphological-physiological traits, because the physiological mechanisms behind animal personality can have broad multitrait effects and because many selective agents influence the evolution of multiple types of traits. However, we currently know little about genetic correlations between animal personalities and nonbehavioural traits. We tested for associations between personality, morphology and locomotor performance by comparing zebrafish (Danio rerio) collected from the wild and then selectively bred for either a proactive or reactive stress coping style ('bold' or 'shy' phenotypes). Based on adaptive hypotheses of correlational selection in the wild, we predicted that artificial selection for boldness would produce correlated evolutionary responses of larger caudal regions and higher fast-start escape performance (and the opposite for shyness). After four to seven generations, morphology and locomotor performance differed between personality lines: bold zebrafish exhibited a larger caudal region and higher fast-start performance than fish in the shy line, matching predictions. Individual-level phenotypic correlations suggested that pleiotropy or physical gene linkage likely explained the correlated response of locomotor performance, while the correlated response of body shape may have reflected linkage disequilibrium, which is breaking down each generation in the laboratory. Our results indicate that evolution of personality can result in concomitant changes in morphology and whole-organism performance, and vice versa.

Incipient sympatric speciation in Midas cichlid fish from the youngest and one of the smallest crater lakes in Nicaragua due to differential use of the benthic and limnetic habitats?

Understanding how speciation can occur without geographic isolation remains a central objective in evolutionary biology. Generally, some form of disruptive selection and assortative mating are necessary for sympatric speciation to occur. Disruptive selection can arise from intraspecific competition for resources. If this competition leads to the differential use of habitats and variation in relevant traits is genetically determined, then assortative mating can be an automatic consequence (i.e., habitat isolation). In this study, we caught Midas cichlid fish from the limnetic (middle of the lake) and benthic (shore) habitats of Crater Lake Asososca Managua to test whether some of the necessary conditions for sympatric speciation due to intraspecific competition and habitat isolation are given. Lake As. Managua is very small (<900 m in diameter), extremely young (maximally 1245 years of age), and completely isolated. It is inhabited by, probably, only a single endemic species of Midas cichlids, Amphilophus tolteca. We found that fish from the limnetic habitat were more elongated than fish collected from the benthic habitat, as would be predicted from ecomorphological considerations. Stable isotope analyses confirmed that the former also exhibit a more limnetic lifestyle than the latter. Furthermore, split‐brood design experiments in the laboratory suggest that phenotypic plasticity is unlikely to explain much of the observed differences in body elongation that we observed in the field. Yet, neutral markers (microsatellites) did not reveal any genetic clustering in the population. Interestingly, demographic inferences based on RAD‐seq data suggest that the apparent lack of genetic differentiation at neutral markers could simply be due to a lack of time, as intraspecific competition may only have begun a few hundred generations ago.

Multispecies Outcomes of Sympatric Speciation after Admixture with the Source Population in Two Radiations of Nicaraguan Crater Lake Cichlids

The formation of species in the absence of geographic barriers (i.e. sympatric speciation) remains one of the most controversial topics in evolutionary biology. While theoretical models have shown that this most extreme case of primary divergence-with-gene-flow is possible, only a handful of accepted empirical examples exist. And even for the most convincing examples uncertainties remain; complex histories of isolation and secondary contact can make species falsely appear to have originated by sympatric speciation. This alternative scenario is notoriously difficult to rule out. Midas cichlids inhabiting small and remote crater lakes in Nicaragua are traditionally considered to be one of the best examples of sympatric speciation and lend themselves to test the different evolutionary scenarios that could lead to apparent sympatric speciation since the system is relatively small and the source populations known. Here we reconstruct the evolutionary history of two small-scale radiations of Midas cichlids inhabiting crater lakes Apoyo and Xiloá through a comprehensive genomic data set. We find no signs of differential admixture of any of the sympatric species in the respective radiations. Together with coalescent simulations of different demographic models our results support a scenario of speciation that was initiated in sympatry and does not result from secondary contact of already partly diverged populations. Furthermore, several species seem to have diverged simultaneously, making Midas cichlids an empirical example of multispecies outcomes of sympatric speciation. Importantly, however, the demographic models strongly support an admixture event from the source population into both crater lakes shortly before the onset of the radiations within the lakes. This opens the possibility that the formation of reproductive barriers involved in sympatric speciation was facilitated by genetic variants that evolved in a period of isolation between the initial founding population and the secondary migrants that came from the same source population. Thus, the exact mechanisms by which these species arose might be different from what had been thought before.

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.

Fine-scale behavioural differences distinguish resource use by ecomorphs in a closed ecosystem

Temporal differences in habitat use and foraging specialisms between ecomorphs represent aspects of behavioural phenotype that are poorly understood with regard to the origin and maintenance of ecological diversity. We tested the role of behaviour in resource use divergence of two Arctic charr (Salvelinus alpinus) phenotypes, a slim, putatively pelagic-dwelling morph and a robust, putatively littoral-dwelling generalist morph, over an annual cycle, using biotelemetry and stable isotopes. Pelagic morph charr exhibited significantly greater δC¹³ depletion, concordant with increased zooplanktivory, than for the Littoral morph. Although three-dimensional space-use of the morphs strongly overlapped, on average, the Littoral morph used that habitat 19.3% more than the Pelagic morph. Pelagic morph fish were significantly more active, further from the lake bed and at greater depth than Littoral fish (annual means respectively, Pelagic, 0.069BLs⁻¹, 8.21 m and 14.11 m; Littoral, 0.047BLs⁻¹, 5.87 m and 10.47 m). Patterns of habitat use differed between ecomorphs at key times, such as during autumn and at ice break, likely related to spawning and resumption of intensive foraging respectively. Extensive space-use overlap, but fine-scale differences in habitat use between charr ecomorphs, suggests the importance of competition for generating and maintaining polymorphism, and its potential for promoting reproductive isolation and evolution in sympatry.

Tunicates: exploring the sea shores and roaming the open ocean. A tribute to Thomas Huxley

This review is a tribute to the remarkable contributions of Thomas Huxley to the biology of tunicates, the likely sister group of vertebrates. In 1851, the great biologist and philosopher published two landmark papers on pelagic tunicates in the Philosophical Transactions of the Royal Society. They were dedicated to the description of the adult anatomy and life cycle of thaliaceans and appendicularians, the pelagic relatives of ascidians. In the first part of this review, we discuss the novel anatomical observations and evolutionary hypotheses made by Huxley, which would have a lasting influence on tunicate biology. We also briefly comment on the more philosophical reflections of Huxley on individuality. In the second part, we stress the originality and relevance of past and future studies of tunicates in the resolution of major biological issues. In particular, we focus on the complex relationship between genotype and phenotype and the phenomenon of developmental system drift. We propose that more than 150 years after Huxley's papers, tunicate embryos are still worth studying in their own right, independently of their evolutionary proximity to vertebrates, as they provide original and crucial insights into the process of animal evolution. Tunicates are still at the forefront of biological research.

Descriptions of five new species of the Neotropical cichlid genus Gymnogeophagus Miranda Ribeiro, 1918 (Teleostei: Cichliformes) from the rio Uruguay drainage

ABSTRACT Gymnogeophagus is a Neotropical cichlid genus distributed in the río Paraguay, rio Paraná and rio Uruguay drainages and also in the coastal rivers of Uruguay and southern Brazil. Its monophyly is supported by two derived features: the absence of supraneurals and the presence of a forward spine in the first dorsal-fin pterygiophore. Herein, five new species of Gymnogeophagus are described from middle to upper tributaries of the rio Uruguay drainage and from the rio Negro. All these new species belong to a clade, which includes G. gymnogenys , easily recognized by sharing two synapomorphies: the absence of an oblique bar between the dorsal border of the eye and the nape, and the possession of a black bar originating in the dorsal contour near the dorsal-fin origin and directed downward and backward on the dorsum. A key to all species of Gymnogeophagus is provided.

Ecological diversification associated with the pharyngeal jaw diversity of Neotropical cichlid fishes

Innovations can facilitate bursts of diversification by increasing access to novel resources and the attainment of novel functional designs. Pharyngognathy, exhibited by highly diverse groups such as wrasses and cichlid fishes, is hypothesized to increase foraging capacity and efficiency. Here, I test the hypothesis that pharyngeal jaw shape and tooth morphology are adaptive in an ecologically diverse radiation of Neotropical cichlid fishes that spans North, Central and South America. I partitioned species into generalized trophic guilds using published stomach content analyses and quantified shape variation of the lower pharyngeal jaw (LPJ) using geometric morphometrics. Additionally, I tested for convergence in LPJ shape and trophic guild by mapping the phylogeny onto the principal components and testing for shifts towards similar evolutionary regimes. Major LPJ shape variation included the length and orientation (i.e. narrow or wide) of the lateral processes and length of the medial process, which varied based on the proportion of fishes and plants consumed. Pharyngeal tooth number, diversity and the frequency of tooth types were not evenly distributed among trophic guilds. There were seven distinct evolutionary regimes that converged upon four optima. Pharyngeal jaw diversification is associated with the exploitation of novel resources among Neotropical cichlids such that pharyngeal specialization has increased access to otherwise poorly accessible resources, such as resources that are difficult to crush (e.g. hard-shelled organisms) and assimilate (e.g. algae).

Sympatric ecological divergence associated with a color polymorphism

Background

Color polymorphisms are a conspicuous feature of many species and a way to address broad ecological and evolutionary questions. Three potential major evolutionary fates of color polymorphisms are conceivable over time: maintenance, loss, or speciation. However, the understanding of color polymorphisms and their evolutionary implications is frequently impaired by sex-linkage of coloration, unknown inheritance patterns, difficulties in phenotypic characterization, and a lack of evolutionary replicates. Hence, the role of color polymorphisms in promoting ecological and evolutionary diversification remains poorly understood. In this context, we assessed the ecological and evolutionary consequences of a color polymorphic study system that is not hampered by these restrictions: the repeated adaptive radiations of the gold/dark Midas cichlid fishes (the Amphilophus citrinellus species complex) from the great lakes and crater lakes of Nicaragua, Central America.

Results

We conducted multi-trait morphological and ecological analyses from ten populations of this young adaptive radiation (<6,000 years old), which revealed sympatric ecological differentiation associated with the conspicuous binary (gold/dark) color polymorphism. Varying degrees of intraspecific ecological divergence were observed across the ten color morph pairs, but most pairs exhibited a consistently parallel ecological and evolutionary trajectory across populations. Specifically, gold Midas cichlids are frequently deeper-bodied, have more robust pharyngeal jaws, and feed at a lower trophic level compared to conspecific, sympatric dark individuals. A common garden experiment suggests there is a genetic correlation of color and eco-morphological traits.

Conclusions

We demonstrate unprecedented ecological and evolutionary consequences of color polymorphism in this adaptive radiation. Across the species complex, sympatric conspecific individuals differed in eco-morphology depending on color morph (gold/dark) and the axis of differentiation tended to be consistent across replicates. The consistent divergence across wild populations and the common garden experiment suggests that color is genetically correlated to ecology. Because Midas cichlids are known to mate color assortatively, the putative genetic correlation of this color polymorphism with an eco-morphological divergence suggests an innate potential to promote ecological and evolutionary divergence across this species complex. However, there are to date no examples of speciation based on color in this radiation, suggesting long-term maintenance of this color polymorphism.

Electronic supplementary material

The online version of this article (doi:10.1186/s12915-015-0192-7) contains supplementary material, which is available to authorized users.

High levels of interspecific gene flow in an endemic cichlid fish adaptive radiation from an extreme lake environment

Studying recent adaptive radiations in isolated insular systems avoids complicating causal events and thus may offer clearer insight into mechanisms generating biological diversity. Here, we investigate evolutionary relationships and genomic differentiation within the recent radiation of Alcolapia cichlid fish that exhibit extensive phenotypic diversification, and which are confined to the extreme soda lakes Magadi and Natron in East Africa. We generated an extensive RAD data set of 96 individuals from multiple sampling sites and found evidence for genetic admixture between species within Lake Natron, with the highest levels of admixture between sympatric populations of the most recently diverged species. Despite considerable environmental separation, populations within Lake Natron do not exhibit isolation by distance, indicating panmixia within the lake, although individuals within lineages clustered by population in phylogenomic analysis. Our results indicate exceptionally low genetic differentiation across the radiation despite considerable phenotypic trophic variation, supporting previous findings from smaller data sets; however, with the increased power of densely sampled SNPs, we identify genomic peaks of differentiation (FST outliers) between Alcolapia species. While evidence of ongoing gene flow and interspecies hybridization in certain populations suggests that Alcolapia species are incompletely reproductively isolated, the identification of outlier SNPs under diversifying selection indicates the radiation is undergoing adaptive divergence.

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.