Ancient and recent riverine gene flow contributed to the adaptive radiation of sailfin silversides in Wallace's Dreampond

While adaptive radiations significantly contribute to the worlds biodiversity, much is unknown about the genetic and ecological factors underlying these rapid successions of speciation. It has been suggested that hybridisation can facilitate the speciation process by generating genetic diversity on which diversifying selection can act. Sailfin silverside fishes (Telmatherinidae) in the Malili Lakes system in Sulawesi have diversified within the last 2 million years. We assembled and annotated a chromosome-scale reference genome of the riverine sailfin silverside Telmatherina bonti and generated whole genome sequences of all species of Telmatherina in Lake Matano, South Sulawesi, Indonesia, one of the worlds oldest and deepest lakes. We reconstructed the phylogenetic relationships within the adaptive radiation of sailfin silversides and inferred past and ongoing introgression patterns. Genome-wide tests confirmed two monophyletic clades, sharpfins and roundfins. However, within clades, we found mismatches between morphology-based taxonomic assignments and genome-wide genetic relationships. We found signs of both old and ongoing introgression between river-dwelling T. bonti and the lacustrine sharpfin group, as shown in elevated D-statistic, f4-ratio and f-branch statistic. Levels of excess allele sharing between riverine species and the three most common lacustrine species declined with increasing distance from the river-inlet, indicating ongoing introgression at the lake-river interface. This combination of past and ongoing hybridisation in a radiating species flock makes Lake Matano Telmatherina a particularly valuable system to study fundamental mechanisms driving rapid speciation under genomic exchange. The phylogenomic framework elaborated in this study provides the foundation for studies of the processes shaping this charismatic radiation.