Reticulate sympatric speciation in Cameroonian crater lake cichlids

Museum Genomics Reveals the Hybrid Origin of an Extinct Crater Lake Endemic

Crater lake fishes are common evolutionary model systems, with recent studies suggesting a key role for gene flow in promoting rapid adaptation and speciation. However, the study of these young lakes can be complicated by human-mediated extinctions. Museum genomics approaches integrating genetic data from recently extinct species are, therefore, critical to understanding the complex evolutionary histories of these fragile systems. Here, we examine the evolutionary history of an extinct Southern Hemisphere crater lake endemic, the rainbowfish Melanotaenia eachamensis. We undertook a comprehensive sampling of extant rainbowfish populations of the Atherton Tablelands of Australia alongside historical museum material to understand the evolutionary origins of the extinct crater lake population and the dynamics of gene flow across the ecoregion. The extinct crater lake species is genetically distinct from all other nearby populations due to historic introgression between 2 proximate riverine lineages, similar to other prominent crater lake speciation systems, but this historic gene flow has not been sufficient to induce a species flock. Our results suggest that museum genomics approaches can be successfully combined with extant sampling to unravel complex speciation dynamics involving recently extinct species.

Cryptic diversity and speciation in an endemic copepod crustacean Harpacticella inopinata within Lake Baikal

Ancient lakes are hotspots of species diversity, posing challenges and opportunities for exploration of the dynamics of endemic diversification. Lake Baikal in Siberia, the oldest lake in the world, hosts a particularly rich crustacean fauna, including the largest known species flock of harpacticoid copepods with some 70 species. Here, we focused on exploring the diversity and evolution within a single nominal species, Harpacticella inopinata Sars, 1908, using molecular markers (mitochondrial COI, nuclear ITS1 and 28S rRNA) and a set of qualitative and quantitative morphological traits. Five major mitochondrial lineages were recognized, with model-corrected COI distances of 0.20-0.37. A concordant pattern was seen in the nuclear data set, and qualitative morphological traits also distinguish a part of the lineages. All this suggests the presence of several hitherto unrecognized cryptic taxa within the baikalian H. inopinata, with long independent histories. The abundances, distributions and inferred demographic histories were different among taxa. Two taxa, H. inopinata CE and H. inopinata CW, were widespread on the eastern and western coasts, respectively, and were largely allopatric. Patterns in mitochondrial variation, that is, shallow star-like haplotype networks, suggest these taxa have spread through the lake relatively recently. Three other taxa, H. inopinata RE, RW and RW2, instead were rare and had more localized distributions on either coast, but showed deeper intraspecies genealogies, suggesting older regional presence. The rare taxa were often found in sympatry with the others and occasionally introgressed by mtDNA from the common ones. The mitochondrial divergence between and within the H. inopinata lineages is still unexpectedly deep, suggesting an unusually high molecular rate. The recognition of true systematic diversity in the evaluation and management of ecosystems is important in hotspots, as it is everywhere else, while the translation of the diversity into a formal taxonomy remains a challenge.

Ecological disruptive selection acting on quantitative loci can drive sympatric speciation

The process of speciation generates biodiversity. According to the null model of speciation, barriers between populations arise in allopatry, where, prior to biology, geography imposes barriers to gene flow. On the other hand, sympatric speciation requires that the process of speciation happen in the absence of a geographical barrier, where the members of the population have no spatial, temporal barriers. Several attempts have been made to theoretically identify the conditions in which speciation can occur in sympatry. However, these efforts suffer from several limitations. We propose a model for sympatric speciation based on adaptation for resource utilization. We use a genetics-based model to investigate the relative roles of prezygotic and postzygotic barriers, from the context of ecological disruptive selection, sexual selection, and genetic architecture, in causing and maintaining sympatric speciation. Our results show that sexual selection that acts on secondary sexual traits does not play any role in the process of speciation in sympatry and that assortative mating based on an ecologically relevant trait forces the population to show an adaptive response. We also demonstrate that understanding the genetic architecture of the trait under ecological selection is very important and that it is not required for the strength of ecological disruptive selection to be very high in order for speciation to occur in sympatry. Our results provide an insight into the kind of scenarios in which sympatric speciation can be demonstrated in the lab.

Sequential hybridization may have facilitated ecological transitions in the Southwestern pinyon pine syngameon

Multispecies interbreeding networks, or syngameons, have been increasingly reported in natural systems. However, the formation, structure, and maintenance of syngameons have received little attention. Through gene flow, syngameons can increase genetic diversity, facilitate the colonization of new environments, and contribute to hybrid speciation. In this study, we evaluated the history, patterns, and consequences of hybridization in a pinyon pine syngameon using morphological and genomic data to assess genetic structure, demographic history, and geographic and climatic data to determine niche differentiation. We demonstrated that Pinus edulis, a dominant species in the Southwestern US and a barometer of climate change, is a core participant in the syngameon, involved in the formation of two drought-adapted hybrid lineages including the parapatric and taxonomically controversial fallax-type. We found that species remain morphologically and genetically distinct at range cores, maintaining species boundaries while undergoing extensive gene flow in areas of sympatry at range peripheries. Our study shows that sequential hybridization may have caused relatively rapid speciation and facilitated the colonization of different niches, resulting in the rapid formation of two new lineages. Participation in the syngameon may allow adaptive traits to be introgressed across species barriers and provide the changes needed to survive future climate scenarios.

Trophic specialization on unique resources despite limited niche divergence in a celebrated example of sympatric speciation

Trophic niche partitioning is observed in many adaptive radiations and is hypothesized to be a central process underlying species divergence. However, patterns of dietary niche partitioning are inconsistent across radiations and there are few studies of niche partitioning in putative examples of sympatric speciation. Here, we conducted the first quantitative study of dietary niche partitioning using stomach contents and stable isotope analyses in one of the most celebrated examples of sympatric speciation: the cichlid radiation from crater lake Barombi Mbo, Cameroon. We found little evidence for trophic niche partitioning among cichlids, including the nine species coexisting in the narrow littoral zone. Stable isotope analyses supported these conclusions of substantial dietary overlap. Our data, however, did reveal that five of eleven species consume rare dietary items, including freshwater sponge, terrestrial ants, and nocturnal foraging on shrimp. Stomach contents of the spongivore (Pungu maclareni) were 20% freshwater sponge, notable considering that only 0.04% of all fishes consume sponges. Overall, we conclude that cichlid species in lake Barombi Mbo overlap considerably in broad dietary niches-in part due to the large proportion of detritus in the stomach contents of all species-but there is evidence for divergence among species in their diet specializations on unique resources. We speculate that these species may utilize these additional specialized resources during periods of low resource abundance in support of Liem's paradox.

Exon-based phylogenomics and the relationships of African cichlid fishes: tackling the challenges of reconstructing phylogenies with repeated rapid radiations

African cichlids (subfamily: Pseudocrenilabrinae) are among the most diverse vertebrates, and their propensity for repeated rapid radiation has made them a celebrated model system in evolutionary research. Nonetheless, despite numerous studies, phylogenetic uncertainty persists, and riverine lineages remain comparatively underrepresented in higher-level phylogenetic studies. Heterogeneous gene histories resulting from incomplete lineage sorting (ILS) and hybridization are likely sources of uncertainty, especially during episodes of rapid speciation. We investigate relationships of Pseudocrenilabrinae and its close relatives while accounting for multiple sources of genetic discordance using species tree and hybrid network analyses with hundreds of single-copy exons. We improve sequence recovery for distant relatives, thereby extending the taxonomic reach of our probes, with a hybrid reference guided/de novo assembly approach. Our analyses provide robust hypotheses for most higher-level relationships and reveal widespread gene heterogeneity, including in riverine taxa. ILS and past hybridization are identified as sources of genetic discordance in different lineages. Sampling of various Blenniiformes (formerly Ovalentaria) adds strong phylogenomic support for convict blennies (Pholidichthyidae) as sister to Cichlidae, and points to other potentially useful protein-coding markers across the order. A reliable phylogeny with representatives from diverse environments will support ongoing taxonomic and comparative evolutionary research in the cichlid model system.

The Syngameon Enigma

Despite their evolutionary relevance, multispecies networks or syngameons are rarely reported in the literature. Discovering how syngameons form and how they are maintained can give insight into processes such as adaptive radiations, island colonizations, and the creation of new hybrid lineages. Understanding these complex hybridization networks is even more pressing with anthropogenic climate change, as syngameons may have unique synergistic properties that will allow participating species to persist. The formation of a syngameon is not insurmountable, as several ways for a syngameon to form have been proposed, depending mostly on the magnitude and frequency of gene flow events, as well as the relatedness of its participants. Episodic hybridization with small amounts of introgression may keep syngameons stable and protect their participants from any detrimental effects of gene flow. As genomic sequencing becomes cheaper and more species are included in studies, the number of known syngameons is expected to increase. Syngameons must be considered in conservation efforts as the extinction of one participating species may have detrimental effects on the survival of all other species in the network.

Intergenus F1-hybrids of African weakly electric fish (Mormyridae: Gnathonemus petersii ♂ × Campylomormyrus compressirostris ♀) are fertile

Hybridisation is an important element of adaptive radiation in fish but data are limited in weakly electric mormyrid fish in this respect. Recently, it has been shown that intragenus hybrids (Campylomormyrus) are fertile and are able to produce F2-fish. In this paper, we demonstrate that even intergenus hybrids (Gnathonemus petersii ♂ × Campylomormyrus compressirostris ♀) are fertile. Three artificial reproduction (AR) trials, with an average fertilisation rate of ca. 23%, yielded different numbers of survivals (maximally about 50%) of the F1-hybrids. The complete ontogenetic development of these hybrids is described concerning their morphology and electric organ discharge (EOD). Two EOD types emerged at the juvenile stage, which did not change up to adulthood. Type I consisted of four phases and Type II was triphasic. The minimum body length at sexual maturity was between 10 and 11 cm. Malformations, growth and mortality rates are also described.

Rapid generation of ecologically relevant behavioral novelty in experimental cichlid hybrids

The East African cichlid radiations are characterized by repeated and rapid diversification into many distinct species with different ecological specializations and by a history of hybridization events between nonsister species. Such hybridization might provide important fuel for adaptive radiation. Interspecific hybrids can have extreme trait values or novel trait combinations and such transgressive phenotypes may allow some hybrids to explore ecological niches neither of the parental species could tap into. Here, we investigate the potential of second-generation (F2) hybrids between two generalist cichlid species from Lake Malawi to exploit a resource neither parental species is specialized on: feeding by sifting sand. Some of the F2 hybrids phenotypically resembled fish of species that are specialized on sand sifting. We combined experimental behavioral and morphometric approaches to test whether the F2 hybrids are transgressive in both morphology and behavior related to sand sifting. We then performed a quantitative trait loci (QTL) analysis using RADseq markers to investigate the genetic architecture of morphological and behavioral traits. We show that transgression is present in several morphological traits, that novel trait combinations occur, and we observe transgressive trait values in sand sifting behavior in some of the F2 hybrids. Moreover, we find QTLs for morphology and for sand sifting behavior, suggesting the existence of some loci with moderate to large effects. We demonstrate that hybridization has the potential to rapidly generate novel and ecologically relevant phenotypes that may be suited to a niche neither of the parental species occupies. Interspecific hybridization may thereby contribute to the rapid generation of ecological diversity in cichlid radiations.

Intragenus F1-hybrids of African weakly electric fish (Mormyridae: Campylomormyrus tamandua ♂ × C. compressirostris ♀) are fertile

Hybridization is widespread in fish and constitutes an important mechanism in fish speciation. There is, however, little knowledge about hybridization in mormyrids. F1-interspecies hybrids between Campylomormyrus tamandua ♂ × C. compressirostris ♀ were investigated concerning: (1) fertility; (2) survival of F2-fish and (3) new gene combinations in the F2-generation concerning the structure of the electric organ and features of the electric organ discharge. These F1-hybrids achieved sexual maturity at about 12–13.5 cm total length. A breeding group comprising six males and 13 females spawned 28 times naturally proving these F1-fish to be fertile. On average 228 eggs were spawned, the average fertilization rate was 47.8%. Eggs started to hatch 70–72 h after fertilization, average hatching rate was 95.6%. Average mortality rate during embryonic development amounted to 2.3%. Average malformation rate during the free embryonic stage was 27.7%. Exogenous feeding started on day 11. In total, we raised 353 normally developed larvae all of which died consecutively, the oldest specimen reaching an age of 5 months. During survival, the activities of the larval and adult electric organs were recorded and the structure of the adult electric organ was investigated histologically.

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.

Interspecific hybridization can generate functional novelty in cichlid fish

The role of interspecific hybridization in evolution is still being debated. Interspecific hybridization has been suggested to facilitate the evolution of ecological novelty, and hence the invasion of new niches and adaptive radiation when ecological opportunity is present beyond the parental species niches. On the other hand, hybrids between two ecologically divergent species may perform less well than parental species in their respective niches because hybrids would be intermediate in performance in both niches. The evolutionary consequences of hybridization may hence be context-dependent, depending on whether ecological opportunities, beyond those of the parental species, do or do not exist. Surprisingly, these complementary predictions may never have been tested in the same experiment in animals. To do so, we investigate if hybrids between ecologically distinct cichlid species perform less well than the parental species when feeding on food either parent is adapted to, and if the same hybrids perform better than their parents when feeding on food none of the species are adapted to. We generated two first-generation hybrid crosses between species of African cichlids. In feeding efficiency experiments we measured the performance of hybrids and parental species on food types representing both parental species niches and additional 'novel' niches, not used by either of the parental species but by other species in the African cichlid radiations. We found that hybrids can have higher feeding efficiencies on the 'novel' food types but typically have lower efficiencies on parental food types when compared to parental species. This suggests that hybridization can generate functional variation that can be of ecological relevance allowing the access to resources outside of either parental species niche. Hence, we provide support for the hypothesis of ecological context-dependency of the evolutionary impact of interspecific hybridization.

Adaptive Radiation from a Chromosomal Perspective: Evidence of Chromosome Set Stability in Cichlid Fishes (Cichlidae: Teleostei) from the Barombi Mbo Lake, Cameroon

Cichlid fishes are the subject of scientific interest because of their rapid adaptive radiation, resulting in extensive ecological and taxonomic diversity. In this study, we examined 11 morphologically distinct cichlid species endemic to Barombi Mbo, the largest crater lake in western Cameroon, namely Konia eisentrauti, Konia dikume, Myaka myaka, Pungu maclareni, Sarotherodon steinbachi, Sarotherodon lohbergeri, Sarotherodon linnellii, Sarotherodon caroli, Stomatepia mariae, Stomatepia pindu, and Stomatepia mongo. These species supposedly evolved via sympatric ecological speciation from a common ancestor, which colonized the lake no earlier than one million years ago. Here we present the first comparative cytogenetic analysis of cichlid species from Barombi Mbo Lake using both conventional (Giemsa staining, C-banding, and CMA₃/DAPI staining) and molecular (fluorescence in situ hybridization with telomeric, 5S, and 28S rDNA probes) methods. We observed stability on both macro and micro-chromosomal levels. The diploid chromosome number was 2n = 44, and the karyotype was invariably composed of three pairs of meta/submetacentric and 19 pairs of subtelo/acrocentric chromosomes in all analysed species, with the same numbers of rDNA clusters and distribution of heterochromatin. The results suggest the evolutionary stability of chromosomal set; therefore, the large-scale chromosomal rearrangements seem to be unlikely associated with the sympatric speciation in Barombi Mbo.

Patterns, Mechanisms and Genetics of Speciation in Reptiles and Amphibians

In this contribution, the aspects of reptile and amphibian speciation that emerged from research performed over the past decade are reviewed. First, this study assesses how patterns and processes of speciation depend on knowing the taxonomy of the group in question, and discuss how integrative taxonomy has contributed to speciation research in these groups. This study then reviews the research on different aspects of speciation in reptiles and amphibians, including biogeography and climatic niches, ecological speciation, the relationship between speciation rates and phenotypic traits, and genetics and genomics. Further, several case studies of speciation in reptiles and amphibians that exemplify many of these themes are discussed. These include studies of integrative taxonomy and biogeography in South American lizards, ecological speciation in European salamanders, speciation and phenotypic evolution in frogs and lizards. The final case study combines genomics and biogeography in tortoises. The field of amphibian and reptile speciation research has steadily moved forward from the assessment of geographic and ecological aspects, to incorporating other dimensions of speciation, such as genetic mechanisms and evolutionary forces. A higher degree of integration among all these dimensions emerges as a goal for future research.

East African cichlid lineages (Teleostei: Cichlidae) might be older than their ancient host lakes: new divergence estimates for the east African cichlid radiation

BACKGROUND

Cichlids are a prime model system in evolutionary research and several of the most prominent examples of adaptive radiations are found in the East African Lakes Tanganyika, Malawi and Victoria, all part of the East African cichlid radiation (EAR). In the past, great effort has been invested in reconstructing the evolutionary and biogeographic history of cichlids (Teleostei: Cichlidae). In this study, we present new divergence age estimates for the major cichlid lineages with the main focus on the EAR based on a dataset encompassing representative taxa of almost all recognized cichlid tribes and ten mitochondrial protein genes. We have thoroughly re-evaluated both fossil and geological calibration points, and we included the recently described fossil †Tugenchromis pickfordi in the cichlid divergence age estimates.

RESULTS

Our results estimate the origin of the EAR to Late Eocene/Early Oligocene (28.71 Ma; 95% HPD: 24.43-33.15 Ma). More importantly divergence ages of the most recent common ancestor (MRCA) of several Tanganyika cichlid tribes were estimated to be substantially older than the oldest estimated maximum age of the Lake Tanganyika: Trematocarini (16.13 Ma, 95% HPD: 11.89-20.46 Ma), Bathybatini (20.62 Ma, 95% HPD: 16.88-25.34 Ma), Lamprologini (15.27 Ma; 95% HPD: 12.23-18.49 Ma). The divergence age of the crown haplochromine H-lineage is estimated to 22.8 Ma (95% HPD: 14.40-26.32 Ma) and of the Lake Malawi radiation to 4.07 Ma (95% HDP: 2.93-5.26 Ma). In addition, we recovered a novel lineage within the Lamprologini tribe encompassing only Lamprologus of the lower and central Congo drainage with its divergence estimated to the Late Miocene or early Pliocene. Furthermore we recovered two novel mitochondrial haplotype lineages within the Haplochromini tribe: 'Orthochromis' indermauri and 'Haplochormis' vanheusdeni.

CONCLUSIONS

Divergence time estimates of the MRCA of several Tanganyika cichlid tribes predate the age of the extant Lake Tanganyika basin, and hence are in line with the recently formulated "Melting-Pot Tanganyika" hypothesis. The radiation of the 'Lower Congo Lamprologus clade' might be linked with the Pliocene origin of the modern lower Congo rapids as has been shown for other Lower Congo cichlid assemblages. Finally, the age of origin of the Lake Malawi cichlid flock agrees well with the oldest age estimate for lacustrine conditions in Lake Malawi.

Don't throw out the sympatric speciation with the crater lake water: fine-scale investigation of introgression provides equivocal support for causal role of secondary gene flow in one of the clearest examples of sympatric speciation

Genomic data has revealed complex histories of colonization and repeated gene flow previously unrecognized in some of the most celebrated examples of sympatric speciation and radiation. However, much of the evidence for secondary gene flow into these radiations comes from summary statistics calculated from sparse genomic sampling without knowledge of which specific genomic regions introgressed. This tells us little about how gene flow potentially influenced sympatric diversification. Here, we investigated whole genomes of Barombi Mbo crater lake cichlids for fine-scale patterns of introgression with neighboring riverine cichlid populations. We found evidence of secondary gene flow into the radiation scattered across <0.24% of the genome; however, from our analyses, it is not clear if the functional diversity in these regions contributed to the ecological, sexual, and morphological diversity found in the lake. Unlike similar studies, we found no obvious candidate genes for adaptive introgression and we cannot rule out that secondary gene flow was predominantly neutral with respect to the diversification process. We also found evidence for differential assortment of ancestral polymorphisms found in riverine populations between sympatric sister species, suggesting the presence of an ancestral hybrid swarm. Although the history of gene flow and colonization is more complicated than previously assumed, the lack of compelling evidence for secondary gene flow's role in species diversification suggests that we should not yet rule out one of the most celebrated examples of sympatric speciation in nature without a more thorough investigation of the timing and functional role of each introgressed region.

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.

Ancient hybridization fuels rapid cichlid fish adaptive radiations

Understanding why some evolutionary lineages generate exceptionally high species diversity is an important goal in evolutionary biology. Haplochromine cichlid fishes of Africa's Lake Victoria region encompass >700 diverse species that all evolved in the last 150,000 years. How this 'Lake Victoria Region Superflock' could evolve on such rapid timescales is an enduring question. Here, we demonstrate that hybridization between two divergent lineages facilitated this process by providing genetic variation that subsequently became recombined and sorted into many new species. Notably, the hybridization event generated exceptional allelic variation at an opsin gene known to be involved in adaptation and speciation. More generally, differentiation between new species is accentuated around variants that were fixed differences between the parental lineages, and that now appear in many new combinations in the radiation species. We conclude that hybridization between divergent lineages, when coincident with ecological opportunity, may facilitate rapid and extensive adaptive radiation.

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.

Demographic modelling with whole-genome data reveals parallel origin of similar Pundamilia cichlid species after hybridization

Modes and mechanisms of speciation are best studied in young species pairs. In older taxa, it is increasingly difficult to distinguish what happened during speciation from what happened after speciation. Lake Victoria cichlids in the genus Pundamilia encompass a complex of young species and polymorphic populations. One Pundamilia species pair, P. pundamilia and P. nyererei, is particularly well suited to study speciation because sympatric population pairs occur with different levels of phenotypic differentiation and reproductive isolation at different rocky islands within the lake. Genetic distances between allopatric island populations of the same nominal species often exceed those between the sympatric species. It thus remained unresolved whether speciation into P. nyererei and P. pundamilia occurred once, followed by geographical range expansion and interspecific gene flow in local sympatry, or if the species pair arose repeatedly by parallel speciation. Here, we use genomic data and demographic modelling to test these alternative evolutionary scenarios. We demonstrate that gene flow plays a strong role in shaping the observed patterns of genetic similarity, including both gene flow between sympatric species and gene flow between allopatric populations, as well as recent and early gene flow. The best supported model for the origin of P. pundamilia and P. nyererei population pairs at two different islands is one where speciation happened twice, whereby the second speciation event follows shortly after introgression from an allopatric P. nyererei population that arose earlier. Our findings support the hypothesis that very similar species may arise repeatedly, potentially facilitated by introgressed genetic variation.

Cichlids of the Banc d'Arguin National Park, Mauritania: insight into the diversity of the genus Coptodon

To determine the species diversity of cichlids in the Banc d'Arguin National Park (PNBA) and their phylogenetic relationships with other species in West Africa, a morphometric and meristic and molecular phylogenetic study was conducted. Both approaches not only confirm the presence of Sarotherodon melanotheron in PNBA but also demonstrate the presence of a second species from the genus Coptodon. While morphometric characteristics match the description of the Guinean tilapia Coptodon guineensis, phylogenetic reconstructions based on three mitochondrial and one nuclear DNA fragment demonstrate that C. guineensis is paraphyletic over its range. Because different lineages of C. guineensis are allopatric, the distribution of C. guineensis should be restricted to Ghana and Côte d'Ivoire. The many other lineages of this species should be considered as C. sp. aff. guineensis.

HEXT, a software supporting tree-based screens for hybrid taxa in multilocus data sets, and an evaluation of the homoplasy excess test

The homoplasy excess test (HET) is a tree-based screen for hybrid taxa in multilocus nuclear phylogenies. Homoplasy between a hybrid taxon and the clades containing the parental taxa reduces bootstrap support in the tree. The HET is based on the expectation that excluding the hybrid taxon from the data set increases the bootstrap support for the parental clades, whereas excluding non-hybrid taxa has little effect on statistical node support. To carry out a HET, bootstrap trees are calculated with taxon-jackknife data sets, that is excluding one taxon (species, population) at a time. Excess increase in bootstrap support for certain nodes upon exclusion of a particular taxon indicates the hybrid (the excluded taxon) and its parents (the clades with increased support).We introduce a new software program, hext, which generates the taxon-jackknife data sets, runs the bootstrap tree calculations, and identifies excess bootstrap increases as outlier values in boxplot graphs. hext is written in r language and accepts binary data (0/1; e.g. AFLP) as well as co-dominant SNP and genotype data.We demonstrate the usefulness of hext in large SNP data sets containing putative hybrids and their parents. For instance, using published data of the genus Vitis (~6,000 SNP loci), hext output supports V. × champinii as a hybrid between V. rupestris and V. mustangensis.With simulated SNP and AFLP data sets, excess increases in bootstrap support were not always connected with the hybrid taxon (false positives), whereas the expected bootstrap signal failed to appear on several occasions (false negatives). Potential causes for both types of spurious results are discussed.With both empirical and simulated data sets, the taxon-jackknife output generated by hext provided additional signatures of hybrid taxa, including changes in tree topology across trees, consistent effects of exclusions of the hybrid and the parent taxa, and moderate (rather than excessive) increases in bootstrap support. hext significantly facilitates the taxon-jackknife approach to hybrid taxon detection, even though the simple test for excess bootstrap increase may not reliably identify hybrid taxa in all applications.

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.

Lake Tanganyika--a 'melting pot' of ancient and young cichlid lineages (Teleostei: Cichlidae)?

A long history of research focused on the East Africa cichlid radiations (EAR) revealed discrepancies between mtDNA and nuclear phylogenies, suggesting that interspecific hybridisation may have been significant during the radiation of these fishes. The approximately 250 cichlid species of Lake Tanganyika have their roots in a monophyletic African cichlid assemblage, but controversies remain about the precise phylogenetic origin and placement of different lineages and consequently about L. Tanganyika colonization scenarios. 3312 AFLP loci and the mitochondrial ND2 gene were genotyped for 91 species representing almost all major lacustrine and riverine haplotilapiine east African cichlid lineages with a focus on L. Tanganyika endemics. Explicitly testing for the possibility of ancient hybridisation events, a comprehensive phylogenetic network hypothesis is proposed for the origin and diversification of L. Tanganyika cichlids. Inference of discordant phylogenetic signal strongly suggests that the genomes of two endemic L. Tanganyika tribes, Eretmodini and Tropheini, are composed of an ancient mixture of riverine and lacustrine lineages. For the first time a strong monophyly signal of all non-haplochromine mouthbrooding species endemic to L. Tanganyika (“ancient mouthbrooders”) was detected. Further, in the genomes of early diverging L. Tanganyika endemics Trematocarini, Bathybatini, Hemibatini and Boulengerochromis genetic components of other lineages belonging to the East African Radiation appear to be present. In combination with recent palaeo-geological results showing that tectonic activity in the L. Tanganyika region resulted in highly dynamic and heterogeneous landscape evolution over the Neogene and Pleistocene, the novel phylogenetic data render a single lacustrine basin as the geographical cradle of the endemic L. Tanganyika cichlid lineages unlikely. Instead a scenario of a pre-rift origin of several independent L. Tanganyika precursor lineages which diversified in ancient rivers and precursor lakes and then amalgamated in the extant L. Tanganyika basin is put forward as an alternative: the 'melting pot Tanganyika' hypothesis.

Morphology, molecules, and monogenean parasites: an example of an integrative approach to cichlid biodiversity

The unparalleled biodiversity of Lake Tanganyika (Africa) has fascinated biologists for over a century; its unique cichlid communities are a preferred model for evolutionary research. Although species delineation is, in most cases, relatively straightforward, higher-order classifications were shown not to agree with monophyletic groups. Here, traditional morphological methods meet their limitations. A typical example are the tropheine cichlids currently belonging to Simochromis and Pseudosimochromis. The affiliations of these widespread and abundant cichlids are poorly understood. Molecular work suggested that genus and species boundaries should be revised. Moreover, previous morphological results indicated that intraspecific variation should be considered to delineate species in Lake Tanganyika cichlids. We review the genera Simochromis and Pseudosimochromis using an integrative approach. Besides a morphometric study and a barcoding approach, monogenean Cichlidogyrus (Platyhelminthes: Ancyrocephalidae) gill parasites, often highly species-specific, are used as complementary markers. Six new species are described. Cichlidogyrus raeymaekersi sp. nov., C. muterezii sp. nov. and C. banyankimbonai sp. nov. infect S. diagramma. Cichlidogyrus georgesmertensi sp. nov. was found on S. babaulti and S. pleurospilus, C. franswittei sp. nov. on both S. marginatus and P. curvifrons and C. frankwillemsi sp. nov. only on P. curvifrons. As relatedness between Cichlidogyrus species usually reflects relatedness between hosts, we considered Simochromis monotypic because the three Cichlidogyrus species found on S. diagramma belonged to a different morphotype than those found on the other Simochromis. The transfer of S. babaulti, S. marginatus, S. pleurospilus and S. margaretae to Pseudosimochromis was justified by the similarity of their Cichlidogyrus fauna and the intermediate morphology of S. margaretae. Finally parasite data also supported the synonymy between S. pleurospilus and S. babaulti, a species that contains a large amount of geographical morphological variation.

Hybridisation and diversification in the adaptive radiation of clownfishes

BACKGROUND

The importance of hybridisation during species diversification has long been debated among evolutionary biologists. It is increasingly recognised that hybridisation events occurred during the evolutionary history of numerous species, especially during the early stages of adaptive radiation. We study the effect of hybridisation on diversification in the clownfishes, a clade of coral reef fish that diversified through an adaptive radiation process. While two species of clownfish are likely to have been described from hybrid specimens, the occurrence and effect of hybridisation on the clade diversification is yet unknown.

RESULTS

We generate sequences of three mitochondrial genes to complete an existing dataset of nuclear sequences and document cytonuclear discordance at a node, which shows a drastic increase of diversification rate. Then, using a tree-based jack-knife method, we identify clownfish species likely stemming from hybridisation events. Finally, we use molecular cloning and identify the putative parental species of four clownfish specimens that display the morphological characteristics of hybrids.

CONCLUSIONS

Our results show that consistently with the syngameon hypothesis, hybridisation events are linked with a burst of diversification in the clownfishes. Moreover, several recently diverged clownfish lineages likely originated through hybridisation, which indicates that diversification, catalysed by hybridisation events, may still be happening.

Speciation in Freshwater Fishes

The extraordinary species richness of freshwater fishes has attracted much research on mechanisms and modes of speciation. We here review research on speciation in freshwater fishes in light of speciation theory, and place this in a context of broad-scale diversity patterns in freshwater fishes. We discuss several major repeated themes in freshwater fish speciation and the speciation mechanisms they are frequently associated with. These include transitions between marine and freshwater habitats, transitions between discrete freshwater habitats, and ecological transitions within habitats, as well as speciation without distinct niche shifts. Major research directions in the years to come include understanding the transition from extrinsic environment-dependent to intrinsic reproductive isolation and its influences on species persistence and understanding the extrinsic and intrinsic constraints to speciation and how these relate to broad-scale diversification patterns through time.

Population differentiation of zander (Sander lucioperca) across native and newly colonized ranges suggests increasing admixture in the course of an invasion

In addition to ecological factors, evolutionary processes can determine the invasion success of a species. In particular, genetic admixture has the potential to induce rapid evolutionary change, which can result from natural or human-assisted secondary contact between differentiated populations. We studied the recent range expansion of zander in Germany focusing on the interplay between invasion and genetic admixture. Historically, the rivers Elbe and Danube harboured the most north-western source populations from which a north-westward range expansion occurred. This was initiated by introducing zander outside its native range into rivers and lakes, and was fostered by migration through artificial canals and stocking from various sources. We analysed zander populations of the native and invaded ranges using nuclear and mitochondrial genetic markers. Three genetic lineages were identified, which were traced to ancestral ranges. Increased genetic diversity and admixture in the invaded region highlighted asymmetric gene flow towards this area. We suppose that the adaptive potential of the invading populations was promoted by genetic admixture, whereas competitive exclusion in the native areas provided a buffer against introgression by novel genotypes. These explanations would be in line with evidence that hybridization can drive evolutionary change under conditions when new niches can be exploited.

Multiple ITS copies reveal extensive hybridization within Rheum (Polygonaceae), a genus that has undergone rapid radiation

Background

During adaptive radiation events, characters can arise multiple times due to parallel evolution, but transfer of traits through hybridization provides an alternative explanation for the same character appearing in apparently non-sister lineages. The signature of hybridization can be detected in incongruence between phylogenies derived from different markers, or from the presence of two divergent versions of a nuclear marker such as ITS within one individual.

Methodology/Principal Findings

In this study, we cloned and sequenced ITS regions for 30 species of the genus Rheum, and compared them with a cpDNA phylogeny. Seven species contained two divergent copies of ITS that resolved in different clades from one another in each case, indicating hybridization events too recent for concerted evolution to have homogenised the ITS sequences. Hybridization was also indicated in at least two further species via incongruence in their position between ITS and cpDNA phylogenies. None of the ITS sequences present in these nine species matched those detected in any other species, which provides tentative evidence against recent introgression as an explanation. Rheum globulosum, previously indicated by cpDNA to represent an independent origin of decumbent habit, is indicated by ITS to be part of clade of decumbent species, which acquired cpDNA of another clade via hybridization. However decumbent and glasshouse morphology are confirmed to have arisen three and two times, respectively.

Conclusions

These findings suggested that hybridization among QTP species of Rheum has been extensive, and that a role of hybridization in diversification of Rheum requires investigation.

Outgroup effects on root position and tree topology in the AFLP phylogeny of a rapidly radiating lineage of cichlid fish

Phylogenetic analyses of rapid radiations are particularly challenging as short basal branches and incomplete lineage sorting complicate phylogenetic inference. Multilocus data of presence-absence polymorphisms such as obtained by AFLP genotyping overcome some of the difficulties, but also present their own intricacies. Here we analyze >1000 AFLP markers to address the evolutionary history of the Limnochromini, a cichlid fish lineage endemic to Lake Tanganyika, and to test for potential effects of outgroup composition on tree topology. The data support previous mitochondrial evidence on the tribe's taxonomy by confirming the polyphyly of the genus Limnochromis and - in contradiction to a recent taxonomic revision - nesting the genus Greenwoodochromis within the Limnochromini. Species relationships suggest that ecological segregation occurred during the rapid basal radiation of the Limnochromini. The large phylogenetic distance between candidate outgroup taxa and the Limnochromini radiation caused random outgroup effects. Bootstrap support for ingroup nodes was lower in outgroup-rooted than in midpoint-rooted trees, and root positions and ingroup tree topologies varied in response to the composition of the outgroup. These observations suggest that the predisposition for homoplastic evolution makes AFLP-based phylogenetic analyses particularly susceptible to random biases introduced by too-distant outgroup taxa.

Functional basis of ecological divergence in sympatric stickleback

Background

The evolution of ecological divergence in closely related species is a key component of adaptive radiation. However, in most examples of adaptive radiation the mechanistic basis of ecological divergence remains unclear. A classic example is seen in the young benthic and limnetic stickleback species pairs of British Columbia. In each pair the benthic species feeds on littoral macroinvertebrates whereas the limnetic feeds on pelagic zooplankton. Previous studies indicate that in both short-term feeding trials and long-term enclosure studies, benthics and limnetics exhibit enhanced performance on their own resource but fare more poorly on the other species’ resource. We examined the functional basis of ecological divergence in the stickleback species pair from Paxton Lake, BC, using biomechanical models of fish feeding applied to morphological traits. We examined the consequences of morphological differences using high speed video of feeding fish.

Results

Benthic stickleback possess morphological traits that predict high suction generation capacity, including greatly hypertrophied epaxial musculature. In contrast, limnetic stickleback possess traits thought to enhance capture of evasive planktonic prey, including greater jaw protrusion than benthics and greater displacement advantage in both the lower jaw-opening lever system and the opercular four-bar linkage. Kinematic data support the expectations from the morphological analysis that limnetic stickleback exhibit faster strikes and greater jaw protrusion than benthic fish, whereas benthics exert greater suction force on attached prey.

Conclusions

We reveal a previously unknown suite of complex morphological traits that affect rapid ecological divergence in sympatric stickleback. These results indicate that postglacial divergence in stickleback involves many functional systems and shows the value of investigating the functional consequences of phenotypic divergence in adaptive radiation.

Molecular and fossil evidence place the origin of cichlid fishes long after Gondwanan rifting

Cichlid fishes are a key model system in the study of adaptive radiation, speciation and evolutionary developmental biology. More than 1600 cichlid species inhabit freshwater and marginal marine environments across several southern landmasses. This distributional pattern, combined with parallels between cichlid phylogeny and sequences of Mesozoic continental rifting, has led to the widely accepted hypothesis that cichlids are an ancient group whose major biogeographic patterns arose from Gondwanan vicariance. Although the Early Cretaceous (ca 135 Ma) divergence of living cichlids demanded by the vicariance model now represents a key calibration for teleost molecular clocks, this putative split pre-dates the oldest cichlid fossils by nearly 90 Myr. Here, we provide independent palaeontological and relaxed-molecular-clock estimates for the time of cichlid origin that collectively reject the antiquity of the group required by the Gondwanan vicariance scenario. The distribution of cichlid fossil horizons, the age of stratigraphically consistent outgroup lineages to cichlids and relaxed-clock analysis of a DNA sequence dataset consisting of 10 nuclear genes all deliver overlapping estimates for crown cichlid origin centred on the Palaeocene (ca 65-57 Ma), substantially post-dating the tectonic fragmentation of Gondwana. Our results provide a revised macroevolutionary time scale for cichlids, imply a role for dispersal in generating the observed geographical distribution of this important model clade and add to a growing debate that questions the dominance of the vicariance paradigm of historical biogeography.

Speciation patterns and processes in the zooplankton of the ancient lakes of Sulawesi Island, Indonesia

Although studies of ancient lake fauna have provided important insights about speciation patterns and processes of organisms in heterogeneous benthic environments, evolutionary forces responsible for speciation in the relatively homogenous planktonic environment remain largely unexplored. In this study, we investigate possible mechanisms of speciation in zooplankton using the freshwater diaptomids of the ancient lakes of Sulawesi, Indonesia, as a model system. We integrate phylogenetic and population genetic analyses of mitochondrial and nuclear genes with morphological and genome size data. Overall, our results support the conclusion that colonization order and local adaptation are dominant at the large, island scale, whereas at local and intralacustrine scales, speciation processes are regulated by gene flow among genetically differentiated and locally adapted populations. In the Malili lakes, the diaptomid populations are homogenous at nuclear loci, but show two highly divergent mitochondrial clades that are geographically restricted to single lakes despite the interconnectivity of the lake systems. Our study, based on coalescent simulations and population genetic analyses, indicates that unidirectional hybridization allows gene flow across the nuclear genome, but prevents the introgression of mitochondria into downstream populations. We suggest that hybridization and introgression between young lineages is a significant evolutionary force in freshwater plankton.

Successive invasion-mediated interspecific hybridizations and population structure in the endangered cichlid Oreochromis mossambicus

Hybridization between invasive and native species accounts among the major and pernicious threats to biodiversity. The Mozambique tilapia Oreochromis mossambicus, a widely used freshwater aquaculture species, is especially imperiled by this phenomenon since it is recognized by the IUCN as an endangered taxon due to genetic admixture with O. niloticus an invasive congeneric species. The Lower Limpopo and the intermittent Changane River (Mozambique) drain large wetlands of potentially great importance for conservation of O. mossambicus, but their populations have remained unstudied until today. Therefore we aimed (1) to estimate the autochthonous diversity and population structure among genetically pure O. mossambicus populations to provide a baseline for the conservation genetics of this endangered species, (2) to quantify and describe genetic variation of the invasive populations and investigate the most likely factors influencing their spread, (3) to identify O. mossambicus populations unaffected by hybridization. Bayesian assignment tests based on 423 AFLP loci and the distribution of 36 species-specific mitochondrial haplotypes both indicate a low frequency of invasive and hybrid genotypes throughout the system, but nevertheless reveal evidence for limited expansion of two alien species (O. niloticus and O. andersonii) and their hybrids in the Lower Limpopo. O. mossambicus populations with no traces of hybridization are identified. They exhibit a significant genetic structure. This contrasts with previously published estimates and provides rather promising auspices for the conservation of O. mossambicus. Especially, parts of the Upper Changane drainage and surrounding wetlands are identified as refugial zones for O. mossambicus populations. They should therefore receive high conservation priority and could represent valuable candidates for the development of aquaculture strains based on local genetic resources.

Repeated trans-watershed hybridization among haplochromine cichlids (Cichlidae) was triggered by Neogene landscape evolution

The megadiverse haplochromine cichlid radiations of the East African lakes, famous examples of explosive speciation and adaptive radiation, are according to recent studies, introgressed by different riverine lineages. This study is based on the first comprehensive mitochondrial and nuclear DNA dataset from extensive sampling of riverine haplochromine cichlids. It includes species from the lower River Congo and Angolan (River Kwanza) drainages. Reconstruction of phylogenetic hypotheses revealed the paradox of clearly discordant phylogenetic signals. Closely related mtDNA haplotypes are distributed thousands of kilometres apart and across major African watersheds, whereas some neighbouring species carry drastically divergent mtDNA haplotypes. At shallow and deep phylogenetic layers, strong signals of hybridization are attributed to the complex Late Miocene/Early Pliocene palaeohistory of African rivers. Hybridization of multiple lineages across changing watersheds shaped each of the major haplochromine radiations in lakes Tanganyika, Victoria, Malawi and the Kalahari Palaeolakes, as well as a miniature species flock in the Congo basin (River Fwa). On the basis of our results, introgression occurred not only on a spatially restricted scale, but massively over almost the whole range of the haplochromine distribution. This provides an alternative view on the origin and exceptional high diversity of this enigmatic vertebrate group.

Weak disruptive selection and incomplete phenotypic divergence in two classic examples of sympatric speciation: cameroon crater lake cichlids

Recent documentation of a few compelling examples of sympatric speciation led to a proliferation of theoretical models. Unfortunately, plausible examples from nature have rarely been used to test model predictions, such as the initial presence of strong disruptive selection. Here I estimated the form and strength of selection in two classic examples of sympatric speciation: radiations of Cameroon cichlids restricted to Lakes Barombi Mbo and Ejagham. I measured five functional traits and relative growth rates in over 500 individuals within incipient species complexes from each lake. Disruptive selection was prevalent in both groups on single and multivariate trait axes but weak relative to stabilizing selection on other traits and most published estimates of disruptive selection. Furthermore, despite genetic structure, assortative mating, and bimodal species-diagnostic coloration, trait distributions were unimodal in both species complexes, indicating the earliest stages of speciation. Long waiting times or incomplete sympatric speciation may result when disruptive selection is initially weak. Alternatively, I present evidence of additional constraints in both species complexes, including weak linkage between coloration and morphology, reduced morphological variance aligned with nonlinear selection surfaces, and minimal ecological divergence. While other species within these radiations show complete phenotypic separation, morphological and ecological divergence in these species complexes may be slow or incomplete outside optimal parameter ranges, in contrast to rapid divergence of their sexual coloration.

Physiological adaptation along environmental gradients and replicated hybrid zone structure in swordtails (Teleostei: Xiphophorus)

Local adaptation is often invoked to explain hybrid zone structure, but empirical evidence of this is generally rare. Hybrid zones between two poeciliid fishes, Xiphophorus birchmanni and X. malinche, occur in multiple tributaries with independent replication of upstream-to-downstream gradients in morphology and allele frequencies. Ecological niche modelling revealed that temperature is a central predictive factor in the spatial distribution of pure parental species and their hybrids and explains spatial and temporal variation in the frequency of neutral genetic markers in hybrid populations. Among populations of parentals and hybrids, both thermal tolerance and heat-shock protein expression vary strongly, indicating that spatial and temporal structure is likely driven by adaptation to local thermal environments. Therefore, hybrid zone structure is strongly influenced by interspecific differences in physiological mechanisms for coping with the thermal environment.