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

Ancient standing genetic variation facilitated the adaptive radiation of Lake Victoria cichlids

Cichlid fishes are textbook examples of explosive speciation and adaptive radiation, providing a great opportunity to understand how the genomic substrate yields extraordinary species diversity. Recently, we performed comparative genomic analyses of three Lake Victoria cichlids to reveal the genomic substrates underlying their rapid speciation and adaptation. We found that long divergent haplotypes derived from large-scale standing genetic variation, which originated before the adaptive radiation of Lake Victoria cichlids, may have contributed to their rapid diversification. In addition, the present study on genomic data from other East African cichlids suggested the reuse of alleles that may have originated in the ancestral lineages of Lake Tanganyika cichlids during cichlid evolution. Therefore, our results highlight that the primary factor that could drive repeated adaptive radiation across East African cichlids was allelic reuse from standing genetic variation to adapt to their own specific environment. In this report, we summarize the main results and discuss the evolutionary mechanisms of cichlids, based on our latest findings.

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.

Genomic fingerprints of palaeogeographic history: The tempo and mode of rift tectonics across tropical Africa has shaped the diversification of the killifish genus Nothobranchius (Teleostei: Cyprinodontiformes)

This paper reports a phylogeny of the African killifishes (Genus Nothobranchius, Order Cyprinodontiformes) informed by five genetic markers (three nuclear, two mitochondrial) of 80 taxa (seven undescribed and 73 of the 92 recognized species). These short-lived annual fishes occupy seasonally wet habitats in central and eastern Africa, and their distribution coincides largely with the East African Rift System (EARS). The fossil dates of sister clades used to constrain a chronometric tree of all sampled Nothobranchius recovered the origin of the genus at ~13.27 Mya. It was followed by the radiations of six principal clades through the Neogene. An ancestral area estimation tested competing biogeographical hypotheses to constrain the ancestral origin of the genus to the Nilo-Sudan Ecoregion, which seeded a mid-Miocene dispersal event into the Coastal ecoregion, followed closely (~10 Mya) by dispersals southward across the Mozambique coastal plain into the Limpopo Ecoregion. Extending westwards across the Tanzanian plateau, a pulse of radiations through the Pliocene were associated with dispersals and fragmentation of wetlands across the Kalahari and Uganda Ecoregions. We interpret this congruence of drainage rearrangements with dispersals and cladogenic events of Nothobranchius to reflect congruent responses to recurrent uplift and rifting. The coevolution of these freshwater fishes and wetlands is attributed to ultimate control by tectonics, as the EARS extended southwards during the Neogene. Geobiological consilience of the combined evidence supports a tectonic hypothesis for the evolution of Nothobranchius.

New haplochromine cichlid from the upper Miocene (9-10 MYA) of Central Kenya

BACKGROUND

The diversification process known as the Lake Tanganyika Radiation has given rise to the most speciose clade of African cichlids. Almost all cichlid species found in the lakes Tanganyika, Malawi and Victoria, comprising a total of 12-16 tribes, belong to this clade. Strikingly, all the species in the latter two lakes are members of the tribe Haplochromini, whose origin remains unclear. The 'out of Tanganyika' hypothesis argues that the Haplochromini emerged simultaneously with other cichlid tribes and lineages in Lake Tanganyika, presumably about 5-6 million years ago (MYA), and that their presence in the lakes Malawi and Victoria and elsewhere in Africa today is due to later migrations. In contrast, the 'melting pot Tanganyika hypothesis' postulates that Haplochromini emerged in Africa prior to the formation of Lake Tanganyika, and that their divergence could have begun about 17 MYA. Haplochromine fossils could potentially resolve this debate, but such fossils are extremely rare.

RESULTS

Here we present a new fossil haplochromine from the upper Miocene site Waril (9-10 million years) in Central Kenya. Comparative morphology, supported by Micro-CT imaging, reveals that it bears a unique combination of characters relating to dentition, cranial bones, caudal skeleton and meristic traits. Its most prominent feature is the presence of exclusively unicuspid teeth, with canines in the outer tooth row. †Warilochromis unicuspidatus gen. et sp. nov. shares this combination of characters solely with members of the Haplochromini and its lacrimal morphology indicates a possible relation to the riverine genus Pseudocrenilabrus. Due to its fang-like dentition and non-fusiform body, †W. unicuspidatus gen. et sp. nov. might have employed either a sit-and-pursue or sit-and-wait hunting strategy, which has not been reported for any other fossil haplochromine cichlid.

CONCLUSIONS

The age of the fossil (9-10 MYA) is incompatible with the 'out of Tanganyika' hypothesis, which postulates that the divergence of the Haplochromini began only 5-6 MYA. The presence of this fossil in an upper Miocene palaeolake in the Central Kenya Rift, as well as its predatory lifestyle, indicate that Haplochromini were already an important component of freshwater drainages in East Africa at that time.

The coincidence of ecological opportunity with hybridization explains rapid adaptive radiation in Lake Mweru cichlid fishes

The process of adaptive radiation was classically hypothesized to require isolation of a lineage from its source (no gene flow) and from related species (no competition). Alternatively, hybridization between species may generate genetic variation that facilitates adaptive radiation. Here we study haplochromine cichlid assemblages in two African Great Lakes to test these hypotheses. Greater biotic isolation (fewer lineages) predicts fewer constraints by competition and hence more ecological opportunity in Lake Bangweulu, whereas opportunity for hybridization predicts increased genetic potential in Lake Mweru. In Lake Bangweulu, we find no evidence for hybridization but also no adaptive radiation. We show that the Bangweulu lineages also colonized Lake Mweru, where they hybridized with Congolese lineages and then underwent multiple adaptive radiations that are strikingly complementary in ecology and morphology. Our data suggest that the presence of several related lineages does not necessarily prevent adaptive radiation, although it constrains the trajectories of morphological diversification. It might instead facilitate adaptive radiation when hybridization generates genetic variation, without which radiation may start much later, progress more slowly or never occur.

New fossil cichlid from the middle Miocene of East Africa revealed as oldest known member of the Oreochromini

A new genus and species of fossil cichlid fishes of middle Miocene age (12.5 Ma) is described from the Ngorora fish Lagerstätte (Tugen Hills, Kenya) in the East African Rift Valley. Parsimony analysis of morphological characters using published phylogenetic frameworks for extant cichlids combined with the application of a comprehensive best-fit approach based on morphology was employed to place the new fossil taxon in the phylogenetic context of the African cichlids. The data reveal that the fossil specimens can be assigned to the tribe Oreochromini within the haplotilapiines. †Oreochromimos kabchorensis gen. et sp. nov. shows a mosaic set of characters bearing many similarities to the almost pan-African Oreochromis and the East African lake-endemic Alcolapia. As the striking diversity of present-day African cichlids, with 1100 recognised species, has remained largely invisible in the fossil record, the material described here adds significantly to our knowledge of the Miocene diversity of the group. It effectively doubles the age of a fossil calibration point, which has hitherto been used to calibrate divergence times of the East African cichlids in molecular phylogenetic investigations. Furthermore, the comparative dataset derived from extant cichlids presented here will greatly facilitate the classification of fossil cichlids in future studies.

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.

Whole-genome sequences of Malawi cichlids reveal multiple radiations interconnected by gene flow

The hundreds of cichlid fish species in Lake Malawi constitute the most extensive recent vertebrate adaptive radiation. Here we characterize its genomic diversity by sequencing 134 individuals covering 73 species across all major lineages. The average sequence divergence between species pairs is only 0.1-0.25%. These divergence values overlap diversity within species, with 82% of heterozygosity shared between species. Phylogenetic analyses suggest that diversification initially proceeded by serial branching from a generalist Astatotilapia-like ancestor. However, no single species tree adequately represents all species relationships, with evidence for substantial gene flow at multiple times. Common signatures of selection on visual and oxygen transport genes shared by distantly related deep-water species point to both adaptive introgression and independent selection. These findings enhance our understanding of genomic processes underlying rapid species diversification, and provide a platform for future genetic analysis of the Malawi radiation.

Phylogenomics uncovers early hybridization and adaptive loci shaping the radiation of Lake Tanganyika cichlid fishes

Lake Tanganyika is the oldest and phenotypically most diverse of the three East African cichlid fish adaptive radiations. It is also the cradle for the younger parallel haplochromine cichlid radiations in Lakes Malawi and Victoria. Despite its evolutionary significance, the relationships among the main Lake Tanganyika lineages remained unresolved, as did the general timescale of cichlid evolution. Here, we disentangle the deep phylogenetic structure of the Lake Tanganyika radiation using anchored phylogenomics and uncover hybridization at its base, as well as early in the haplochromine radiation. This suggests that hybridization might have facilitated these speciation bursts. Time-calibrated trees support that the radiation of Tanganyika cichlids coincided with lake formation and that Gondwanan vicariance concurred with the earliest splits in the cichlid family tree. Genes linked to key innovations show signals of introgression or positive selection following colonization of lake habitats and species' dietary adaptations are revealed as major drivers of colour vision evolution. These findings shed light onto the processes shaping the evolution of adaptive radiations.

Disentangling Incomplete Lineage Sorting and Introgression to Refine Species-Tree Estimates for Lake Tanganyika Cichlid Fishes

Adaptive radiation is thought to be responsible for the evolution of a great portion of the past and present diversity of life. Instances of adaptive radiation, characterized by the rapid emergence of an array of species as a consequence to their adaptation to distinct ecological niches, are important study systems in evolutionary biology. However, because of the rapid lineage formation in these groups, and occasional gene flow between the participating species, it is often difficult to reconstruct the phylogenetic history of species that underwent an adaptive radiation. In this study, we present a novel approach for species-tree estimation in rapidly diversifying lineages, where introgression is known to occur, and apply it to a multimarker data set containing up to 16 specimens per species for a set of 45 species of East African cichlid fishes (522 individuals in total), with a main focus on the cichlid species flock of Lake Tanganyika. We first identified, using age distributions of most recent common ancestors in individual gene trees, those lineages in our data set that show strong signatures of past introgression. This led us to formulate three hypotheses of introgression between different lineages of Tanganyika cichlids: the ancestor of Boulengerochromini (or of Boulengerochromini and Bathybatini) received genomic material from the derived H-lineage; the common ancestor of Cyprichromini and Perissodini experienced, in turn, introgression from Boulengerochromini and/or Bathybatini; and the Lake Tanganyika Haplochromini and closely related riverine lineages received genetic material from Cyphotilapiini. We then applied the multispecies coalescent model to estimate the species tree of Lake Tanganyika cichlids, but excluded the lineages involved in these introgression events, as the multispecies coalescent model does not incorporate introgression. This resulted in a robust species tree, in which the Lamprologini were placed as sister lineage to the H-lineage (including the Eretmodini), and we identify a series of rapid splitting events at the base of the H-lineage. Divergence ages estimated with the multispecies coalescent model were substantially younger than age estimates based on concatenation, and agree with the geological history of the Great Lakes of East Africa. Finally, we formally tested the three hypotheses of introgression using a likelihood framework, and find strong support for introgression between some of the cichlid tribes of Lake Tanganyika.

Relic populations of Fukomys mole-rats in Tanzania: description of two new species F. livingstoni sp. nov. and F. hanangensis sp. nov

Previous studies of African mole-rats of the genera Heliophobius and Fukomys (Bathyergidae) in the regions of East and south central Africa have revealed a diversity of species and vicariant populations, with patterns of distribution having been influenced by the geological process of rifting and changing patterns of drainage of major river systems. This has resulted in most of the extant members of the genus Fukomys being distributed west of the main Rift Valley. However, a small number of isolated populations are known to occur east of the African Rift Valley in Tanzania, where Heliophobius is the most common bathyergid rodent. We conducted morphological, craniometric and phylogenetic analysis of mitochondrial cytochrome b (cyt b) sequences of two allopatric populations of Tanzanian mole-rats (genus Fukomys) at Ujiji and around Mount Hanang, in comparison with both geographically adjacent and more distant populations of Fukomys. Our results reveal two distinct evolutionary lineages, forming clades that constitute previously unnamed species. Here, we formally describe and designate these new species F. livingstoni and F. hanangensis respectively. Molecular clock-based estimates of divergence times, together with maximum likelihood inference of biogeographic range evolution, offers strong support for the hypothesis that vicariance in the Western Rift Valley and the drainage patterns of major river systems has subdivided populations of mole-rats. More recent climatic changes and tectonic activity in the “Mbeya triple junction” and Rungwe volcanic province between Lakes Rukwa and Nyasa have played a role in further isolation of these extra-limital populations of Fukomys in Tanzania.

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.

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: A Host of Opportunities for Evolutionary Parasitology

Thanks to high species diversity and a broad range of speciation mechanisms, cichlid fishes represent a textbook model in evolutionary biology. They are also of substantial economic value. Despite this importance, cichlid parasites remain understudied, although some are more diverse than their hosts. They may offer important insights into cichlid evolution and the evolution of host-parasite interactions. We review five major lines of research conducted on cichlid parasites so far: the study of parasite diversity and speciation; the role of parasites in cichlid diversification; the evolutionary ecology of host specificity; historical biogeography; and biological invasions. We call for more research in these areas and suggest approaches to valorise the potential that cichlid parasites hold for the study of evolutionary parasitology.

Deep-water parasite diversity in Lake Tanganyika: description of two new monogenean species from benthopelagic cichlid fishes

Background

Lake Tanganyika is the world’s second deepest lake. Its diverse cichlid assemblage offers a unique opportunity for studying a deep-water host-parasite model in freshwater. Low host specificity and a broad host range including representatives of the Bathybatini tribe in the only monogenean parasite described from this habitat, Cichlidogyrus casuarinus Pariselle, Muterezi Bukinga & Vanhove, 2015 suggest a link between lower specificity and lower host density. Conversely, high host specificity and species richness are reported for monogeneans of the lake’s littoral cichlids. We further investigated whether the deep-water environment in Lake Tanganyika is really monogenean species-depauperate by investigating the monogenean fauna of Trematocara unimaculatum (a representative of the tribe Trematocarini, the sister lineage of the Bathybatini) and Benthochromis horii, a member of the tribe Benthochromini, found in the same deep-water habitat as the already known hosts of C. casuarinus.

Methods

Sclerotised structures of the collected monogenean individuals were characterised morphologically using light microscopy and morphometrics.

Results

Both examined cichlid species are infected by a single monogenean species each, which are new to science. They are described as Cichlidogyrus brunnensis n. sp., infecting T. unimaculatum, and Cichlidogyrus attenboroughi n. sp., parasitising on B. horii. Diagnostic characteristics include the distal bifurcation of the accessory piece in C. brunnensis n. sp. and the combination of long auricles and no heel in C. attenboroughi n. sp. In addition C. brunnensis n. sp. does not resemble C. casuarinus, the only species of Cichlidogyrus thus far reported from the Bathybatini. Also Cichlidogyrus attenboroughi n. sp. does not resemble any of the monogenean species documented from the pelagic zone of the lake and is among the few described species of Cichlidogyrus without heel.

Conclusions

As two new and non-resembling Cichlidogyrus species are described from T. unimaculatum and B. horii, colonisation of the deep-water habitat by more than one morphotype of Cichlidogyrus is evident. Based on morphological comparisons with previously described monogenean species, parasite transfers with the littoral zone are possible. Therefore, parasites of pelagic cichlids in the lake do not seem to only mirror host phylogeny and the evolutionary history of this host-parasite system merits further attention.

A robust phylogeny among major lineages of the East African cichlids

The huge monophyletic group of the East African cichlid radiations (EAR) consists of thousands of species belonging to 12-14 tribes; the number of tribes differs among studies. Many studies have inferred phylogenies of EAR tribes using various genetic markers. However, these phylogenies partly contradict one another and can have weak statistic support. In this study, we conducted maximum-likelihood (ML) phylogenetic analyses using restriction site-associated DNA (RAD) sequences and propose a new robust phylogenetic hypothesis among Lake Tanganyika cichlid fishes, which cover most EAR tribes. Data matrices can vary in size and contents depending on the strategies used to process RAD sequences. Therefore, we prepared 23 data matrices with various processing strategies. The ML phylogenies inferred from 15 large matrices (2.0×10(6) to 1.1×10(7) base pairs) resolved every tribe as a monophyletic group with 100% bootstrap support and shared the same topology regarding relationships among the tribes. Most nodes among the tribes were supported by 100% bootstrap values, and the bootstrap support for the other node varied among the 15 ML trees from 70% to 100%. These robust ML trees differ partly in topology from those in earlier studies, and these phylogenetic relationships have important implications for the tribal classification of EAR.

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.

Molecular phylogenetics reveals convergent evolution in lower Congo River spiny eels

Background

The lower Congo River (LCR) is a region of exceptional species diversity and endemism in the Congo basin, including numerous species of spiny eels (genus Mastacembelus). Four of these exhibit distinctive phenotypes characterized by greatly reduced optic globes deeply embedded into the head (cryptophthalmia) and reduced (or absent) melanin pigmentation, among other characteristics. A strikingly similar cryptophthalmic phenotype is also found in members of a number of unrelated fish families, strongly suggesting the possibility of convergent evolution. However, little is known about the evolutionary processes that shaped diversification in LCR Mastacembelus, their biogeographic origins, or when colonization of the LCR occurred.

Methods

We sequenced mitochondrial and nuclear genes from Mastacembelus species collected in the lower Congo River, and compared them with other African species and Asian representatives as outgroups. We analyzed the sequence data using Maximum Likelihood and Bayesian phylogenetic inference.

Results

Bayesian and Maximum Likelihood phylogenetic analyses, and Bayesian coalescent methods for species tree reconstruction, reveal that endemic LCR spiny eels derive from two independent origins, clearly demonstrating convergent evolution of the cryptophthalmic phenotype. Mastacembelus crassus, M. aviceps, and M. simbi form a clade, allied to species found in southern, eastern and central Africa. Unexpectedly, M. brichardi and brachyrhinus fall within a clade otherwise endemic to Lake Tanganikya (LT) ca. 1500 km east of the LCR. Divergence dating suggests the ages of these two clades of LCR endemics differ markedly. The age of the crassus group is estimated at ~4 Myr while colonization of the LCR by the brichardi-brachyrhinus progenitor was considerably more recent, dated at ~0.5 Myr.

Conclusions

The phylogenetic framework of spiny eels presented here, the first to include LCR species, demonstrates that cryptophthalmia and associated traits evolved at least twice in Mastacembelus: once in M. brichardi and at least once in the M. crassus clade. Timing of diversification is broadly consistent with the onset of modern high-energy flow conditions in the LCR and with previous studies of endemic cichlids. The close genetic relationship between M. brichardi and M. brachyrhinus is particularly notable given the extreme difference in phenotype between these species, and additional work is needed to better understand the evolutionary history of diversification in this clade. The findings presented here demonstrate strong, multi-trait convergence in LCR spiny eels, suggesting that extreme selective pressures have shaped numerous phenotypic attributes of the endemic species of this region.

Electronic supplementary material

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

Molecular phylogenetics reveals convergent evolution in lower Congo River spiny eels

BACKGROUND

The lower Congo River (LCR) is a region of exceptional species diversity and endemism in the Congo basin, including numerous species of spiny eels (genus Mastacembelus). Four of these exhibit distinctive phenotypes characterized by greatly reduced optic globes deeply embedded into the head (cryptophthalmia) and reduced (or absent) melanin pigmentation, among other characteristics. A strikingly similar cryptophthalmic phenotype is also found in members of a number of unrelated fish families, strongly suggesting the possibility of convergent evolution. However, little is known about the evolutionary processes that shaped diversification in LCR Mastacembelus, their biogeographic origins, or when colonization of the LCR occurred.

METHODS

We sequenced mitochondrial and nuclear genes from Mastacembelus species collected in the lower Congo River, and compared them with other African species and Asian representatives as outgroups. We analyzed the sequence data using Maximum Likelihood and Bayesian phylogenetic inference.

RESULTS

Bayesian and Maximum Likelihood phylogenetic analyses, and Bayesian coalescent methods for species tree reconstruction, reveal that endemic LCR spiny eels derive from two independent origins, clearly demonstrating convergent evolution of the cryptophthalmic phenotype. Mastacembelus crassus, M. aviceps, and M. simbi form a clade, allied to species found in southern, eastern and central Africa. Unexpectedly, M. brichardi and brachyrhinus fall within a clade otherwise endemic to Lake Tanganikya (LT) ca. 1500 km east of the LCR. Divergence dating suggests the ages of these two clades of LCR endemics differ markedly. The age of the crassus group is estimated at ~4 Myr while colonization of the LCR by the brichardi-brachyrhinus progenitor was considerably more recent, dated at ~0.5 Myr.

CONCLUSIONS

The phylogenetic framework of spiny eels presented here, the first to include LCR species, demonstrates that cryptophthalmia and associated traits evolved at least twice in Mastacembelus: once in M. brichardi and at least once in the M. crassus clade. Timing of diversification is broadly consistent with the onset of modern high-energy flow conditions in the LCR and with previous studies of endemic cichlids. The close genetic relationship between M. brichardi and M. brachyrhinus is particularly notable given the extreme difference in phenotype between these species, and additional work is needed to better understand the evolutionary history of diversification in this clade. The findings presented here demonstrate strong, multi-trait convergence in LCR spiny eels, suggesting that extreme selective pressures have shaped numerous phenotypic attributes of the endemic species of this region.