POPULATION GENOMIC TESTS OF MODELS OF ADAPTIVE RADIATION IN LAKE VICTORIA REGION CICHLID FISH

Severe Bottleneck Impacted the Genomic Structure of Egg-Eating Cichlids in Lake Victoria

Within 15,000 years, the explosive adaptive radiation of haplochromine cichlids in Lake Victoria, East Africa, generated 500 endemic species. In the 1980s, the upsurge of Nile perch, a carnivorous fish artificially introduced to the lake, drove the extinction of more than 200 endemic cichlids. The Nile perch predation particularly harmed piscivorous cichlids, including paedophages, cichlids eat eggs and fries, which is an example of the unique trophic adaptation seen in African cichlids. Here, aiming to investigate past demographic events possibly triggered by the invasion of Nile perch and the subsequent impacts on the genetic structure of cichlids, we conducted large-scale comparative genomics. We discovered evidence of recent bottleneck events in 4 species, including 2 paedophages, which began during the 1970s to 1980s, and population size rebounded during the 1990s to 2000s. The timing of the bottleneck corresponded to the historical records of endemic haplochromines" disappearance and later resurgence, which is likely associated with the introduction of Nile perch by commercial demand to Lake Victoria in the 1950s. Interestingly, among the 4 species that likely experienced bottleneck, Haplochromis sp. "matumbi hunter," a paedophagous cichlid, showed the most severe bottleneck signatures. The components of shared ancestry inferred by ADMIXTURE suggested a high genetic differentiation between matumbi hunter and other species. In contrast, our phylogenetic analyses highly supported the monophyly of the 5 paedophages, consistent with the results of previous studies. We conclude that high genetic differentiation of matumbi hunter occurred due to the loss of shared genetic components among haplochromines in Lake Victoria caused by the recent severe bottleneck.

Rapid clade divergence and phyletic gradualism in an interacting particle model of sympatric speciation

The coexistence of cladogenesis, i.e., the branching of lineages along an evolutionary tree as observed in the fossil record, and anagenesis, which is the progressive evolution within populations, lacks a clear explanation. In this study, we examine a simple model that simulates the evolutionary changes occurring within populations inhabiting the same environment in sympatry, and driven by ecological competition. Our model characterizes populations through a set of evolving morphological traits represented by mathematical points within a two-dimensional morphospace. Such points may reproduce or die due to overcrowding, implying competition in morphospace as suggested by the ecological phenomenon of character displacement. By focusing on the morphospace rather than physical space, the model effectively captures the simultaneous evolution of coexisting populations. Central to the model is the delicate balance between the range of competition and the range of reproduction within the morphospace. Interesting patterns emerge when the ratio between the competition to reproducetion ranges, referred to as CR ratio, changes from values slightly smaller to significantly larger than unity. When competition acts over short distances relative to the reproduction range (low CR), the phylogenetic tree takes on a nearly uniform appearance, gradually transforming into a more bush-like structure for slightly higher CR values. With further increases in CR, evolutionary lineages become more discernible, and the morphogenetic pattern shifts from a bush-like shape to a more tree-like arrangement and few branches for very large CRs. At specific time sections, the synthetic phylogenetic tree appears as an assembly of clusters of individuals within the morphospace. These clusters, interpretable as simulated models of species, exhibit distinct separation within the morphospace and are subject to dynamic inter-cluster repulsion. Notably, clusters tend to be resistant to change. They maintain relatively constant abundances while gradually shifting their positions within the morphospace-a phase that aligns with the concept of phyletic gradualism. However, this predictable pattern is occasionally upset by the abrupt divisions into multiple groups, interpreted as cladogenesis events. The intricacies of the splitting process are explored, revealing that in scenarios with large CR values, the splitting can emerge much more rapidly than phyletic changes. This accelerated process of splitting is initiated by one or few individuals at the fringes of a cluster, where competition is minimal. The newly generated cluster then undergoes deformation, swiftly followed by divergence and splitting (seen as branching in the synthetic phylogenetic tree), as if an inherent "repulsion" triggered the division between species. The simple rules implied in the interacting-particle model may provide insight into the coexistence of gradualism and cladogenesis along lineages, illustrating the capacity for rapid shifts during cladogenesis and the more gradual process of anagenesis.

Cycles of fusion and fission enabled rapid parallel adaptive radiations in African cichlids

Although some lineages of animals and plants have made impressive adaptive radiations when provided with ecological opportunity, the propensities to radiate vary profoundly among lineages for unknown reasons. In Africa's Lake Victoria region, one cichlid lineage radiated in every lake, with the largest radiation taking place in a lake less than 16,000 years old. We show that all of its ecological guilds evolved in situ. Cycles of lineage fusion through admixture and lineage fission through speciation characterize the history of the radiation. It was jump-started when several swamp-dwelling refugial populations, each of which were of older hybrid descent, met in the newly forming lake, where they fused into a single population, resuspending old admixture variation. Each population contributed a different set of ancient alleles from which a new adaptive radiation assembled in record time, involving additional fusion-fission cycles. We argue that repeated fusion-fission cycles in the history of a lineage make adaptive radiation fast and predictable.

Mitogenomic Characterization of Cameroonian Endemic Coptodon camerunensis (Cichliformes: Cichlidae) and Matrilineal Phylogeny of Old-World Cichlids

The mitogenomic evolution of old-world cichlids is still largely incomplete in Western Africa. In this present study, the complete mitogenome of the Cameroon endemic cichlid, Coptodon camerunensis, was determined by next-generation sequencing. The mitogenome was 16,557 bp long and encoded with 37 genes (13 protein-coding genes, two ribosomal RNA genes, 22 transfer RNA genes, and a control region). The C. camerunensis mitogenome is AT-biased (52.63%), as exhibited in its congener, Coptodon zillii (52.76% and 53.04%). The majority of PCGs start with an ATG initiation codon, except COI, which starts with a GTG codon and five PCGs and ends with the TAA termination codon and except seven PCGs with an incomplete termination codon. In C. camerunensis mitogenome, most tRNAs showed classical cloverleaf secondary structures, except tRNA-serine with a lack of DHU stem. Comparative analyses of the conserved blocks of two Coptodonini species control regions revealed that the CSB-II block was longer than other blocks and contained highly variable sites. Using 13 concatenated PCGs, the mitogenome-based Bayesian phylogeny easily distinguished all the examined old-world cichlids. Except for Oreochromini and Coptodinini tribe members, the majority of the taxa exhibited monophyletic clustering within their respective lineages. C. camerunensis clustered closely with Heterotilapia buttikoferi (tribe Heterotilapiini) and had paraphyletic clustering with its congener, C. zillii. The Oreochromini species also displayed paraphyletic grouping, and the genus Oreochromis showed a close relationship with Coptodinini and Heterotilapiini species. In addition, illustrating the known distribution patterns of old-world cichlids, the present study is congruent with the previous hypothesis and proclaims that prehistoric geological evolution plays a key role in the hydroclimate of the African continent during Mesozoic, which simultaneously disperses and/or colonizes cichlids in different ichthyological provinces and Rift Lake systems in Africa. The present study suggests that further mitogenomes of cichlid species are required, especially from western Africa, to understand their unique evolution and adaptation.

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

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

Identification of a novel sex determining chromosome in cichlid fishes that acts as XY or ZW in different lineages

Sex determination systems are highly conserved among most vertebrates with genetic sex determination, but can be variable and evolve rapidly in some. Here, we study sex determination in a clade with exceptionally high sex chromosome turnover rates. We identify the sex determining chromosomes in three interspecific crosses of haplochromine cichlid fishes from Lakes Victoria and Malawi. We find evidence for different sex determiners in each cross. In the Malawi cross and one Victoria cross the same chromosome is sex-linked but while females are the heterogametic sex in the Malawi species, males are the heterogametic sex in the Victoria species. This chromosome has not previously been reported to be sex determining in cichlids, increasing the number of different chromosomes shown to be sex determining in cichlids to 12. All Lake Victoria species of our crosses are less than 15,000 years divergent, and we identified different sex determiners among them. Our study provides further evidence for the diversity and evolutionary flexibility of sex determination in cichlids, factors which might contribute to their rapid adaptive radiations.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s10750-021-04560-7.

The propagation of admixture-derived adaptive radiation potential

Adaptive radiations (ARs) frequently show remarkable repeatability where single lineages undergo multiple independent episodes of AR in distant places and long-separate time points. Genetic variation generated through hybridization between distantly related lineages can promote AR. This mechanism, however, requires rare coincidence in space and time between a hybridization event and opening of ecological opportunity, because hybridization generates large genetic variation only locally and it will persist only for a short period. Hence, hybridization seems unlikely to explain recurrent AR in the same lineage. Contrary to these expectations, our evolutionary computer simulations demonstrate that admixture variation can geographically spread and persist for long periods if the hybrid population becomes separated into isolated sub-lineages. Subsequent secondary hybridization of some of these can reestablish genetic polymorphisms from the ancestral hybridization in places far from the birthplace of the hybrid clade and long after the ancestral hybridization event. Consequently, simulations revealed conditions where exceptional genetic variation, once generated through a rare hybridization event, can facilitate multiple ARs exploiting ecological opportunities available at distant points in time and space.

Similar ecology, different morphology: Three new species of oral-mollusc shellers from Lake Edward

Lake Edward, East Africa, harbours a largely understudied assemblage of haplochromine cichlids that displays a range of adaptions to various specialised trophic niches. In this system, we discovered specimens of Haplochromis with morphologies similar to those of oral-mollusc shellers from Lake Victoria. These morphologies are characterised by short oral jaws with stout teeth that are used either to crush molluscs or to grab the soft bodies of snails and wrench them out of their shells. A morphometric study on 47 specimens from Lake Edward revealed the presence of three new species with an oral-shelling morphology: Haplochromis concilians sp. nov., H. erutus sp. nov. and H. planus sp. nov. All three species are formally described. Stomach-content observations confirmed an opportunistic oral-shelling ecology for H. concilians sp. nov. and H. erutus sp. nov. Within H. planus sp. nov., only large specimens displayed a specialised oral-shelling morphology, but their stomachs were nearly empty, while small specimens consumed mainly Ostracoda and Hydrachnidia. Remarkably, the three species differed considerably in morphology from each other, but they each resembled oral-sheller species from Lake Victoria.

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.

Divergent parasite infections in sympatric cichlid species in Lake Victoria

Parasitism has been proposed as a factor in host speciation, as an agent affecting coexistence of host species in species-rich communities and as a driver of post-speciation diversification. Young adaptive radiations of closely related host species of varying ecological and genomic differentiation provide interesting opportunities to explore interactions between patterns of parasitism, divergence and coexistence of sympatric host species. Here, we explored patterns in ectoparasitism in a community of 16 fully sympatric cichlid species at Makobe Island in Lake Victoria, a model system of vertebrate adaptive radiation. We asked whether host niche, host abundance or host genetic differentiation explains variation in infection patterns. We found significant differences in infections, the magnitude of which was weakly correlated with the extent of genomic divergence between the host species, but more strongly with the main ecological gradient, water depth. These effects were most evident with infections of Cichlidogyrus monogeneans, whereas the only host species with a strictly crevice-dwelling niche, Pundamilia pundamilia, deviated from the general negative relationship between depth and parasitism. In accordance with the Janzen-Connell hypothesis, we also found that host abundance tended to be positively associated with infections in some parasite taxa. Data on the Pundamilia sister species pairs from three other islands with variable degrees of habitat (crevice) specialization suggested that the lower parasite abundance of P. pundamilia at Makobe could result from both habitat specialization and the evolution of specific resistance. Our results support influences of host genetic differentiation and host ecology in determining infections in this diverse community of sympatric cichlid species.

Pattern of genetic differentiation of an incipient speciation process: The case of the high Andean killifish Orestias

During the Pleistocene and Holocene, the southwest Andean Altiplano (17°-22°S) was affected by repeated fluctuations in water levels, high volcanic activity and major tectonic movements. In the early Holocene the humid Tauca phase shifted to the arid conditions that have lasted until the present, producing endorheic rivers, lakes, lagoons and wetlands. The endemic fish Orestias (Cyprinodontidae) represents a good model to observe the genetic differentiation that characterizes an incipient speciation process in allopatry since the morphospecies described inhabit a restricted geographic area, with present habitat fragmentation. The genetic diversity and population structure of four endemic morphospecies of Orestias (Cyprinodontidae) found in the Lauca National Park (LNP) analyzed with mitochondrial markers (Control Region) and eight microsatellites, revealed the existence of genetic groups that matches the fragmentation of these systems. High values of genetic and phylogeographic differentiation indices were observed between Chungará Lake and Piacota lagoon. The group composed of the Lauca River, Copapujo and Chuviri wetlands sampling sites showed a clear signal of expansion, with a star-like haplotype network. Levels of genetic differentiation were lower than in Chungará and Piacota, suggesting that these localities would have differentiated after the bottlenecks linked to the collapse of Parinacota volcano. The Parinacota sample showed a population signal that differed from the other localities revealing greater genetic diversity and a disperse network, presenting haplotypes shared with other LNP localities. A mixing pattern of the different genetic groups was evident using the microsatellite markers. The chronology of the vicariance events in LNP may indicate that the partition process of the Orestias populations was gradual. Considering this, and in view of the genetic results, we may conclude that the morphospecies from LNP are populations in ongoing differentiation process.

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.

Evaluating genomic divergence and parallelism in replicate ecomorphs from young and old cichlid adaptive radiations

Comparative genomic studies of closely related species typically focus on single species pairs at one given stage of divergence. That makes it difficult to infer the continuum of evolutionary process during speciation and beyond. Here, we use whole-genome resequencing to examine genomic patterns of divergence in three sympatric cichlid species pairs with very similar functional and ecological differentiation, but different ages. We find a strong signature of increasing genomic divergence with time in both the mitochondrial genome and the nuclear genome. In contrast to many other systems, we find that in these cichlids, regions of elevated relative differentiation also exhibit increased absolute differentiation. We detect a signature of convergent evolution in a comparison of outlier regions across all three species pair comparisons, but the extent of it is modest, and regions that are strongly divergent in any one pair tend to be only slightly elevated in the other pairs, consistent with a repeatable but polygenic basis of traits that characterize the ecomorphs. Our results suggest that strong functional phenotypic differentiation, as seen in all three species pairs, is generally associated with a clear signature of genomic divergence, even in the youngest species pair.

Process and pattern in cichlid radiations - inferences for understanding unusually high rates of evolutionary diversification

The cichlid fish radiations in the African Great Lakes differ from all other known cases of rapid speciation in vertebrates by their spectacular trophic diversity and richness of sympatric species, comparable to the most rapid angiosperm radiations. I review factors that may have facilitated these radiations and compare these with insights from recent work on plant radiations. Work to date suggests that it was a coincidence of ecological opportunity, intrinsic ecological versatility and genomic flexibility, rapidly evolving behavioral mate choice and large amounts of standing genetic variation that permitted these spectacular fish radiations. I propose that spatially orthogonal gradients in the fit of phenotypes to the environment facilitate speciation because they allow colonization of alternative fitness peaks during clinal speciation despite local disruptive selection. Such gradients are manifold in lakes because of the interaction of water depth as an omnipresent third spatial dimension with other fitness-relevant variables. I introduce a conceptual model of adaptive radiation that integrates these elements and discuss its applicability to, and predictions for, plant radiations.

Turnover and accumulation of genetic diversity across large time-scale cycles of isolation and connection of populations

Major climatic and geological events but also population history (secondary contacts) have generated cycles of population isolation and connection of long and short periods. Recent empirical and theoretical studies suggest that fast evolutionary processes might be triggered by such events, as commonly illustrated in ecology by the adaptive radiation of cichlid fishes (isolation and reconnection of lakes and watersheds) and in epidemiology by the fast adaptation of the influenza virus (isolation and reconnection in hosts). We test whether cyclic population isolation and connection provide the raw material (standing genetic variation) for species evolution and diversification. Our analytical results demonstrate that population isolation and connection can provide, to populations, a high excess of genetic diversity compared with what is expected at equilibrium. This excess is either cyclic (high allele turnover) or cumulates with time depending on the duration of the isolation and the connection periods and the mutation rate. We show that diversification rates of animal clades are associated with specific periods of climatic cycles in the Quaternary. We finally discuss the importance of our results for macroevolutionary patterns and for the inference of population history from genomic data.

Limited genetic divergence among Australian alpine Poa tussock grasses coupled with regional structuring points to ongoing gene flow and taxonomic challenges

BACKGROUND AND AIMS

While molecular approaches can often accurately reconstruct species relationships, taxa that are incompletely differentiated pose a challenge even with extensive data. Such taxa are functionally differentiated, but may be genetically differentiated only at small and/or patchy regions of the genome. This issue is considered here in Poa tussock grass species that dominate grassland and herbfields in the Australian alpine zone.

METHODS

Previously reported tetraploidy was confirmed in all species by sequencing seven nuclear regions and five microsatellite markers. A Bayesian approach was used to co-estimate nuclear and chloroplast gene trees with an overall dated species tree. The resulting species tree was used to examine species structure and recent hybridization, and intertaxon fertility was tested by experimental crosses.

KEY RESULTS

Species tree estimation revealed Poa gunnii, a Tasmanian endemic species, as sister to the rest of the Australian alpine Poa. The taxa have radiated in the last 0·5-1·2 million years and the non-gunnii taxa are not supported as genetically distinct. Recent hybridization following past species divergence was also not supported. Ongoing gene flow is suggested, with some broad-scale geographic structure within the group.

CONCLUSIONS

The Australian alpine Poa species are not genetically distinct despite being distinguishable phenotypically, suggesting recent adaptive divergence with ongoing intertaxon gene flow. This highlights challenges in using conventional molecular taxonomy to infer species relationships in recent, rapid radiations.

Genetic structure of pelagic and littoral cichlid fishes from Lake Victoria

The approximately 700 species of cichlids found in Lake Victoria in East Africa are thought to have evolved over a short period of time, and they represent one of the largest known examples of adaptive radiation. To understand the processes that are driving this spectacular radiation, we must determine the present genetic structure of these species and elucidate how this structure relates to the ecological conditions that caused their adaptation. We analyzed the genetic structure of two pelagic and seven littoral species sampled from the southeast area of Lake Victoria using sequences from the mtDNA control region and 12 microsatellite loci as markers. Using a Bayesian model-based clustering method to analyze the microsatellite data, we separated these nine species into four groups: one group composed of pelagic species and another three groups composed mainly of rocky-shore species. Furthermore, we found significant levels of genetic variation between species within each group at both marker loci using analysis of molecular variance (AMOVA), although the nine species often shared mtDNA haplotypes. We also found significant levels of genetic variation between populations within species. These results suggest that initial groupings, some of which appear to have been related to habitat differences, as well as divergence between species within groups took place among the cichlid species of Lake Victoria.

How many marker loci are necessary? Analysis of dominant marker data sets using two popular population genetic algorithms

The number of marker loci required to answer a given research question satisfactorily is especially important for dominant markers since they have a lower information content than co-dominant marker systems. In this study, we used simulated dominant marker data sets to determine the number of dominant marker loci needed to obtain satisfactory results from two popular population genetic analyses: STRUCTURE and AMOVA (analysis of molecular variance). Factors such as migration, level of population differentiation, and unequal sampling were varied in the data sets to mirror a range of realistic research scenarios. AMOVA performed well under all scenarios with a modest quantity of markers while STRUCTURE required a greater number, especially when populations were closely related. The popular ΔK method of determining the number of genetically distinct groups worked well when sampling was balanced, but underestimated the true number of groups with unbalanced sampling. These results provide a window through which to interpret previous work with dominant markers and we provide a protocol for determining the number of markers needed for future dominant marker studies.

Peak and persistent excess of genetic diversity following an abrupt migration increase

Genetic diversity is essential for population survival and adaptation to changing environments. Demographic processes (e.g., bottleneck and expansion) and spatial structure (e.g., migration, number, and size of populations) are known to shape the patterns of the genetic diversity of populations. However, the impact of temporal changes in migration on genetic diversity has seldom been considered, although such events might be the norm. Indeed, during the millions of years of a species' lifetime, repeated isolation and reconnection of populations occur. Geological and climatic events alternately isolate and reconnect habitats. We analytically document the dynamics of genetic diversity after an abrupt change in migration given the mutation rate and the number and sizes of the populations. We demonstrate that during transient dynamics, genetic diversity can reach unexpectedly high values that can be maintained over thousands of generations. We discuss the consequences of such processes for the evolution of species based on standing genetic variation and how they can affect the reconstruction of a population's demographic and evolutionary history from genetic data. Our results also provide guidelines for the use of genetic data for the conservation of natural populations.

Genome-wide RAD sequence data provide unprecedented resolution of species boundaries and relationships in the Lake Victoria cichlid adaptive radiation

Although population genomic studies using next generation sequencing (NGS) data are becoming increasingly common, studies focusing on phylogenetic inference using these data are in their infancy. Here, we use NGS data generated from reduced representation genomic libraries of restriction-site-associated DNA (RAD) markers to infer phylogenetic relationships among 16 species of cichlid fishes from a single rocky island community within Lake Victoria's cichlid adaptive radiation. Previous attempts at sequence-based phylogenetic analyses in Victoria cichlids have shown extensive sharing of genetic variation among species and no resolution of species or higher-level relationships. These patterns have generally been attributed to the very recent origin (<15,000 years) of the radiation, and ongoing hybridization between species. We show that as we increase the amount of sequence data used in phylogenetic analyses, we produce phylogenetic trees with unprecedented resolution for this group. In trees derived from our largest data supermatrices (3 to >5.8 million base pairs in width), species are reciprocally monophyletic with high bootstrap support, and the majority of internal branches on the tree have high support. Given the difficulty of the phylogenetic problem that the Lake Victoria cichlid adaptive radiation represents, these results are striking. The strict interpretation of the topologies we present here warrants caution because many questions remain about phylogenetic inference with very large genomic data set and because we can with the current analysis not distinguish between effects of shared ancestry and post-speciation gene flow. However, these results provide the first conclusive evidence for the monophyly of species in the Lake Victoria cichlid radiation and demonstrate the power that NGS data sets hold to resolve even the most difficult of phylogenetic challenges.

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.

The impact of the geologic history and paleoclimate on the diversification of East african cichlids

The cichlid fishes of the East African Great Lakes are the largest extant vertebrate radiation identified to date. These lakes and their surrounding waters support over 2,000 species of cichlid fish, many of which are descended from a single common ancestor within the past 10 Ma. The extraordinary East African cichlid diversity is intricately linked to the highly variable geologic and paleoclimatic history of this region. Greater than 10 Ma, the western arm of the East African rift system began to separate, thereby creating a series of rift basins that would come to contain several water bodies, including the extremely deep Lakes Tanganyika and Malawi. Uplifting associated with this rifting backponded many rivers and created the extremely large, but shallow Lake Victoria. Since their creation, the size, shape, and existence of these lakes have changed dramatically which has, in turn, significantly influenced the evolutionary history of the lakes' cichlids. This paper reviews the geologic history and paleoclimate of the East African Great Lakes and the impact of these forces on the region's endemic cichlid flocks.