A continuous fish fossil record reveals key insights into adaptive radiation.

Adaptive radiations have been instrumental in generating a considerable amount of life's diversity. Ecological opportunity is thought to be a prerequisite for adaptive radiation1, but little is known about the relative importance of species' ecological versatility versus effects of arrival order in determining which lineage radiates2. Palaeontological records that could help answer this are scarce. In Lake Victoria, a large adaptive radiation of cichlid fishes evolved in an exceptionally short and recent time interval3. We present a rich continuous fossil record extracted from a series of long sediment cores along an onshore-offshore gradient. We reconstruct the temporal sequence of events in the assembly of the fish community from thousands of tooth fossils. We reveal arrival order, relative abundance and habitat occupation of all major fish lineages in the system. We show that all major taxa arrived simultaneously as soon as the modern lake began to form. There is no evidence of the radiating haplochromine cichlid lineage arriving before others, nor of their numerical dominance upon colonization; therefore, there is no support for ecological priority effects. However, although many taxa colonized the lake early and several became abundant, only cichlids persisted in the new deep and open-water habitats once these emerged. Because these habitat gradients are also known to have played a major role in speciation, our findings are consistent with the hypothesis that ecological versatility was key to adaptive radiation, not priority by arrival order nor initial numerical dominance.

The ecological and genomic basis of explosive adaptive radiation.

Speciation rates vary considerably among lineages, and our understanding of what drives the rapid succession of speciation events within young adaptive radiations remains incomplete1-11. The cichlid fish family provides a notable example of such variation, with many slowly speciating lineages as well as several exceptionally large and rapid radiations12. Here, by reconstructing a large phylogeny of all currently described cichlid species, we show that explosive speciation is solely concentrated in species flocks of several large young lakes. Increases in the speciation rate are associated with the absence of top predators; however, this does not sufficiently explain explosive speciation. Across lake radiations, we observe a positive relationship between the speciation rate and enrichment of large insertion or deletion polymorphisms. Assembly of 100 cichlid genomes within the most rapidly speciating cichlid radiation, which is found in Lake Victoria, reveals exceptional 'genomic potential'-hundreds of ancient haplotypes bear insertion or deletion polymorphisms, many of which are associated with specific ecologies and shared with ecologically similar species from other older radiations elsewhere in Africa. Network analysis reveals fundamentally non-treelike evolution through recombining old haplotypes, and the origins of ecological guilds are concentrated early in the radiation. Our results suggest that the combination of ecological opportunity, sexual selection and exceptional genomic potential is the key to understanding explosive adaptive radiation.

Testing alternative hypotheses for the decline of cichlid fish in Lake Victoria using fish tooth time series from sediment cores.

Lake Victoria is well known for its high diversity of endemic fish species and provides livelihoods for millions of people. The lake garnered widespread attention during the twentieth century as major environmental and ecological changes modified the fish community with the extinction of approximately 40% of endemic cichlid species by the 1980s. Suggested causal factors include anthropogenic eutrophication, fishing, and introduced non-native species but their relative importance remains unresolved, partly because monitoring data started in the 1970s when changes were already underway. Here, for the first time, we reconstruct two time series, covering the last approximately 200 years, of fish assemblage using fish teeth preserved in lake sediments. Two sediment cores from the Mwanza Gulf of Lake Victoria, were subsampled continuously at an intra-decadal resolution, and teeth were identified to major taxa: Cyprinoidea, Haplochromini, Mochokidae and Oreochromini. None of the fossils could be confidently assigned to non-native Nile perch. Our data show significant decreases in haplochromine and oreochromine cichlid fish abundances that began long before the arrival of Nile perch. Cyprinoids, on the other hand, have generally been increasing. Our study is the first to reconstruct a time series of any fish assemblage in Lake Victoria extending deeper back in time than the past 50 years, helping shed light on the processes underlying Lake Victoria's biodiversity loss.

Ancient DNA is preserved in fish fossils from tropical lake sediments.

Tropical freshwater lakes are well known for their high biodiversity, and particularly the East African Great Lakes are renowned for their adaptive radiation of cichlid fishes. While comparative phylogenetic analyses of extant species flocks have revealed patterns and processes of their diversification, little is known about evolutionary trajectories within lineages, the impacts of environmental drivers, or the scope and nature of now-extinct diversity. Time-structured palaeodata from geologically young fossil records, such as fossil counts and particularly ancient DNA (aDNA) data, would help fill this large knowledge gap. High ambient temperatures can be detrimental to the preservation of DNA, but refined methodology now allows data generation even from very poorly preserved samples. Here, we show for the first time that fish fossils from tropical lake sediments yield endogenous aDNA. Despite generally low endogenous content and high sample dropout, the application of high-throughput sequencing and, in some cases, sequence capture allowed taxonomic assignment and phylogenetic placement of 17% of analysed fish fossils to family or tribe level, including remains which are up to 2700 years old or weigh less than 1 mg. The relationship between aDNA degradation and the thermal age of samples is similar to that described for terrestrial samples from cold environments when adjusted for elevated temperature. Success rates and aDNA preservation differed between the investigated lakes Chala, Kivu and Victoria, possibly caused by differences in bottom water oxygenation. Our study demonstrates that the sediment records of tropical lakes can preserve genetic information on rapidly diversifying fish taxa over time scales of millennia.

The Genetic Population Structure of Lake Tanganyika's Lates Species Flock, an Endemic Radiation of Pelagic Top Predators.

Understanding genetic connectivity plays a crucial role in species conservation decisions, and genetic connectivity is an important component of modern fisheries management. In this study, we investigated the population genetics of four endemic Lates species of Lake Tanganyika (Lates stappersii, L. microlepis, L. mariae, and L. angustifrons) using reduced-representation genomic sequencing methods. We find the four species to be strongly differentiated from one another (mean interspecific FST = 0.665), with no evidence for contemporary admixture. We also find evidence for strong genetic structure within L. mariae, with the majority of individuals from the most southern sampling site forming a genetic group that is distinct from the individuals at other sampling sites. We find evidence for much weaker structure within the other three species (L. stappersii, L. microlepis, and L. angustifrons). Our ability to detect this weak structure despite small and unbalanced sample sizes and imprecise geographic sampling locations suggests the possibility for further structure undetected in our study. We call for further research into the origins of the genetic differentiation in these four species-particularly that of L. mariae-which may be important for conservation and management of this culturally and economically important clade of fishes.

Structural genomic variation leads to genetic differentiation in Lake Tanganyika's sardines.

Identifying patterns in genetic structure and the genetic basis of ecological adaptation is a core goal of evolutionary biology and can inform the management and conservation of species that are vulnerable to population declines exacerbated by climate change. We used reduced-representation genomic sequencing methods to gain a better understanding of genetic structure among and within populations of Lake Tanganyika's two sardine species, Limnothrissa miodon and Stolothrissa tanganicae. Samples of these ecologically and economically important species were collected across the length of Lake Tanganyika, as well as from nearby Lake Kivu, where L. miodon was introduced in 1959. Our results reveal differentiation within both S. tanganicae and L. miodon that is not explained by geography. Instead, this genetic differentiation is due to the presence of large sex-specific regions in the genomes of both species, but involving different polymorphic sites in each species. Our results therefore indicate rapidly evolving XY sex determination in the two species. Additionally, we found evidence of a large chromosomal rearrangement in L. miodon, creating two homokaryotypes and one heterokaryotype. We found all karyotypes throughout Lake Tanganyika, but the frequencies vary along a north-south gradient and differ substantially in the introduced Lake Kivu population. We do not find evidence for significant isolation by distance, even over the hundreds of kilometres covered by our sampling, but we do find shallow population structure.

Movement of transposable elements contributes to cichlid diversity.

African cichlid fishes are a prime model for studying speciation mechanisms. Despite the development of extensive genomic resources, it has been difficult to determine which sources of genetic variation are responsible for cichlid phenotypic variation. One of their most variable phenotypes is visual sensitivity, with some of the largest spectral shifts among vertebrates. These shifts arise primarily from differential expression of seven cone opsin genes. By mapping expression quantitative trait loci (eQTL) in intergeneric crosses of Lake Malawi cichlids, we previously identified four causative genetic variants that correspond to indels in the promoters of either key transcription factors or an opsin gene. In this comprehensive study, we show that these indels are the result of the movement of transposable elements (TEs) that correlate with opsin expression variation across the Malawi flock. In tracking the evolutionary history of these particular indels, we found they are endemic to Lake Malawi, suggesting that these TEs are recently active and are segregating within the Malawi cichlid lineage. However, an independent indel has arisen at a similar genomic location in one locus outside of the Malawi flock. The convergence in TE movement suggests these loci are primed for TE insertion and subsequent deletions. Increased TE mobility may be associated with interspecific hybridization, which disrupts mechanisms of TE suppression. This might provide a link between cichlid hybridization and accelerated regulatory variation. Overall, our study suggests that TEs may be an important driver of key regulatory changes, facilitating rapid phenotypic change and possibly speciation in African cichlids.

The genomic substrate for adaptive radiation in African cichlid fish.

Cichlid fishes are famous for large, diverse and replicated adaptive radiations in the Great Lakes of East Africa. To understand the molecular mechanisms underlying cichlid phenotypic diversity, we sequenced the genomes and transcriptomes of five lineages of African cichlids: the Nile tilapia (Oreochromis niloticus), an ancestral lineage with low diversity; and four members of the East African lineage: Neolamprologus brichardi/pulcher (older radiation, Lake Tanganyika), Metriaclima zebra (recent radiation, Lake Malawi), Pundamilia nyererei (very recent radiation, Lake Victoria), and Astatotilapia burtoni (riverine species around Lake Tanganyika). We found an excess of gene duplications in the East African lineage compared to tilapia and other teleosts, an abundance of non-coding element divergence, accelerated coding sequence evolution, expression divergence associated with transposable element insertions, and regulation by novel microRNAs. In addition, we analysed sequence data from sixty individuals representing six closely related species from Lake Victoria, and show genome-wide diversifying selection on coding and regulatory variants, some of which were recruited from ancient polymorphisms. We conclude that a number of molecular mechanisms shaped East African cichlid genomes, and that amassing of standing variation during periods of relaxed purifying selection may have been important in facilitating subsequent evolutionary diversification.

Genomics of Rapid Incipient Speciation in Sympatric Threespine Stickleback.

Ecological speciation is the process by which reproductively isolated populations emerge as a consequence of divergent natural or ecologically-mediated sexual selection. Most genomic studies of ecological speciation have investigated allopatric populations, making it difficult to infer reproductive isolation. The few studies on sympatric ecotypes have focused on advanced stages of the speciation process after thousands of generations of divergence. As a consequence, we still do not know what genomic signatures of the early onset of ecological speciation look like. Here, we examined genomic differentiation among migratory lake and resident stream ecotypes of threespine stickleback reproducing in sympatry in one stream, and in parapatry in another stream. Importantly, these ecotypes started diverging less than 150 years ago. We obtained 34,756 SNPs with restriction-site associated DNA sequencing and identified genomic islands of differentiation using a Hidden Markov Model approach. Consistent with incipient ecological speciation, we found significant genomic differentiation between ecotypes both in sympatry and parapatry. Of 19 islands of differentiation resisting gene flow in sympatry, all were also differentiated in parapatry and were thus likely driven by divergent selection among habitats. These islands clustered in quantitative trait loci controlling divergent traits among the ecotypes, many of them concentrated in one region with low to intermediate recombination. Our findings suggest that adaptive genomic differentiation at many genetic loci can arise and persist in sympatry at the very early stage of ecotype divergence, and that the genomic architecture of adaptation may facilitate this.

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

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

Ecological disturbance reduces genomic diversity across an Alpine whitefish adaptive radiation.

Genomic diversity is associated with the adaptive potential of a population and thereby impacts the extinction risk of a species during environmental change. However, empirical data on genomic diversity of populations before environmental perturbations are rare and hence our understanding of the impact of perturbation on diversity is often limited. We here assess genomic diversity utilising whole-genome resequencing data from all four species of the Lake Constance Alpine whitefish radiation. Our data covers a period of strong but transient anthropogenic environmental change and permits us to track changes in genomic diversity in all species over time. Genomic diversity became strongly reduced during the period of anthropogenic disturbance and has not recovered yet. The decrease in genomic diversity varies between 18% and 30%, depending on the species. Interspecific allele frequency differences of SNPs located in potentially ecologically relevant genes were homogenized over time. This suggests that in addition to the reduction of genome-wide genetic variation, the differentiation that evolved in the process of adaptation to alternative ecologies between species might have been lost during the ecological disturbance. The erosion of substantial amounts of genomic variation within just a few generations in combination with the loss of potentially adaptive genomic differentiation, both of which had evolved over thousands of years, demonstrates the sensitivity of biodiversity in evolutionary young adaptive radiations towards environmental disturbance. Natural history collections, such as the one used for this study, are instrumental in the assessment of genomic consequences of anthropogenic environmental change. Historical samples enable us to document biodiversity loss against the shifting baseline syndrome and advance our understanding of the need for efficient biodiversity conservation on a global scale.

Anthropogenic Eutrophication Drives Major Food Web Changes in Mwanza Gulf, Lake Victoria.

Discerning ecosystem change and food web dynamics underlying anthropogenic eutrophication and the introduction of non-native species is necessary for ensuring the long-term sustainability of fisheries and lake biodiversity. Previous studies of eutrophication in Lake Victoria, eastern Africa, have focused on the loss of endemic fish biodiversity over the past several decades, but changes in the plankton communities over this same time remain unclear. To fill this gap, we examined sediment cores from a eutrophic embayment, Mwanza Gulf, to determine the timing and magnitude of changes in the phytoplankton and zooplankton assemblages over the past century. Biogeochemical proxies indicate nutrient enrichment began around ~ 1920 CE and led to rapid increases in primary production, and our analysis of photosynthetic pigments revealed three zones: pre-eutrophication (prior to 1920 CE), onset of eutrophication with increases in all pigments (1920-1990 CE), and sustained eutrophication with cyanobacterial dominance (1990 CE-present). Cladoceran remains indicate an abrupt decline in biomass in ~ 1960 CE, in response to the cumulative effects of eutrophication and lake-level rise, preceding the collapse of haplochromine cichlids in the 1980s. Alona and Chydorus, typically benthic littoral taxa, have remained at relatively low abundances since the 1960s, whereas the abundance of Bosmina, typically a planktonic taxon, increased in the 1990s concurrently with the biomass recovery of haplochromine cichlid fishes. Overall, our results demonstrate substantial changes over the past century in the biomass structure and taxonomic composition of Mwanza Gulf phytoplankton and zooplankton communities, providing a historical food web perspective that can help understand the recent changes and inform future resource management decisions in the Lake Victoria ecosystem. Supplementary Information: The online version contains supplementary material available at 10.1007/s10021-024-00908-x.

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.

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.

Strong species structure but weak geographical structure in demersal Lake Victoria cichlids.

Studying phenotypic and genetic differentiation between very young species can be very informative with regard to learning about processes of speciation. Identifying and characterizing genetic species structure and distinguishing it from spatial genetic structure within a species is a prerequisite for this and is often not given sufficient attention. Young radiations of cichlid fish are classical speciation study systems. However, it is only during the past decade that population genomics based on next-generation sequencing has begun to provide the power to resolve species and distinguish speciation from spatial population structure for the youngest of these radiations. The Lake Victoria haplochromine cichlids constitute the youngest large cichlid fish radiation, probably <20,000 years old. Earlier work showed that communities of rocky reef cichlids are composed of many reciprocally monophyletic species despite their very recent origins. Here, we build on this work by studying assemblages of offshore demersal cichlids, adding analyses of within-species spatial structure to the sympatric species structure. We sampled seven multispecies communities along a 6-km-long transect from one side of the Mwanza Gulf to the other side. We investigated whether phenotypically diagnosed putative species are reciprocally monophyletic and whether such monophyly is stable across species geographic ranges. We show that all species are genetically strongly differentiated in sympatry, that they are reciprocally monophyletic, and that monophyly is stable across distribution ranges. We found significant differentiation between geographically distinct populations in two species, but no or weak isolation by distance. We further found subtle but significant morphological differences between all species and a linear relationship between genomic and morphological distance which suggests that differences in morphology begin to accumulate after speciation has already affected genome-wide restrictions of gene flow.

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.

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.

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.

Population genomic tests of models of adaptive radiation in Lake Victoria region cichlid fish.

Adaptive radiation is usually thought to be associated with speciation, but the evolution of intraspecific polymorphisms without speciation is also possible. The radiation of cichlid fish in Lake Victoria (LV) is perhaps the most impressive example of a recent rapid adaptive radiation, with 600+ very young species. Key questions about its origin remain poorly characterized, such as the importance of speciation versus polymorphism, whether species persist on evolutionary time scales, and if speciation happens more commonly in small isolated or in large connected populations. We used 320 individuals from 105 putative species from Lakes Victoria, Edward, Kivu, Albert, Nabugabo and Saka, in a radiation-wide amplified fragment length polymorphism (AFLP) genome scan to address some of these questions. We demonstrate pervasive signatures of speciation supporting the classical model of adaptive radiation associated with speciation. A positive relationship between the age of lakes and the average genomic differentiation of their species, and a significant fraction of molecular variance explained by above-species level taxonomy suggest the persistence of species on evolutionary time scales, with radiation through sequential speciation rather than a single starburst. Finally the large gene diversity retained from colonization to individual species in every radiation suggests large effective population sizes and makes speciation in small geographical isolates unlikely.