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.

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.

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.

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.

Adaptive radiations are an important source of biodiversity and are often characterized by many speciation events in very short succession. It has been proposed that the high speciation rates in these radiations may be fuelled by novel genetic combinations produced in episodes of hybridization among the young species. The role of such hybridization events in the evolutionary history of a group can be investigated by comparing the genealogical relationships inferred from different subsets of loci, but such studies have thus far often been hampered by shallow genetic divergences, especially in young adaptive radiations, and the lack of genome-scale molecular data. Here, we use a genome-wide sampling of SNPs identified within restriction site-associated DNA (RAD) tags to investigate the genomic consistency of patterns of shared ancestry and adaptive divergence among five sympatric cichlid species of two genera, Pundamilia and Mbipia, which form part of the massive adaptive radiation of cichlids in the East African Lake Victoria. Species pairs differ along several axes: male nuptial colouration, feeding ecology, depth distribution, as well as the morphological traits that distinguish the two genera and more subtle morphological differences. Using outlier scan approaches, we identify signals of divergent selection between all species pairs with a number of loci showing parallel patterns in replicated contrasts either between genera or between male colour types. We then create SNP subsets that we expect to be characterized to different extents by selection history and neutral processes and describe phylogenetic and population genetic patterns across these subsets. These analyses reveal very different evolutionary histories for different regions of the genome. To explain these results, we propose at least two intergeneric hybridization events (between Mbipia spp. and Pundamilia spp.) in the evolutionary history of these five species that would have lead to the evolution of novel trait combinations and new species.

Speciation rates vary tremendously among lineages, and our understanding of whatfuels the rapid succession of speciation events within young adaptive radiations remains incomplete [1][2][3][4][5][6][7][8][9][10][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 radiations 12 . By reconstructing a large phylogeny of all described cichlid species, we show that explosive speciation is solely concentrated in several large young lake species flocks. Speciation rate increases are associated with absence of top predators but this is not a sufficient explanation for explosive speciation. Across lake radiations we observe a positive relationship between speciation rate and enrichment with large indel polymorphisms. Assembly of one hundred cichlid genomes within the most rapidly speciating cichlid radiation, found in Lake Victoria, reveals exceptional 'genomic potential' -hundreds of ancient haplotypes bearing indel polymorphisms, many associated with specific ecologies and shared with ecologically similar species from other older radiations elsewhere in Africa. Network analysis reveals fundamentally nontreelike evolution through recombining old haplotypes, with origins of ecological guilds 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.

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.

Man-made habitat fragmentation is a major concern in river ecology and is expected to have particularly detrimental effects on aquatic species with limited dispersal abilities, like the bullhead (Cottus gobio). We used ten microsatellite markers to investigate small-scale patterns of gene flow, current dispersal and neutral genetic diversity in a morphologically diverse river where fragmented and unfragmented sections could be compared. We found high genetic differentiation between sampling sites with a maximum F ST of 0.32 between sites separated by only 35 km. A significant increase of genetic differentiation with geographical distance was observed in the continuous river section as well as in the full dataset which included headwater populations isolated by anthropogenic barriers. Several lines of evidence are consistent with the hypothesis that such barriers completely block upstream movement while downstream dispersal may be little affected. In the unfragmented habitat, dispersal rates were also higher in the direction of water flow than against it. The resulting asymmetry in gene flow likely contributes to the decrease of genetic variation observed from the lower reaches towards the headwaters, which is particularly pronounced in physically isolated populations. Our findings suggest that headwater populations, due to their isolation and low genetic variation, may be particularly vulnerable to extinction.