Population Genetics and Gene Flow in Cyphotilapia frontosa and Cyphotilapia gibberosa Along the East Coast of Lake Tanganyika

The radiation of cichlid species in the East African Great Lakes is remarkable and rapid. The population genetics of two deep-water Cyphotilapia species along the east coast of Lake Tanganyika from Burundi to southern Tanzania was determined using ddRAD-seq. A combination of ADMIXTURE, PCA, genome polarization, and 2D site frequency spectrum analyses confirmed the presence of two species, C. frontosa in the north and C. gibberosa in the south, as documented in other studies. We also found evidence of a potential hybrid zone connecting the two species at a sharp genetic cline centered in the middle of the lake and apparent introgression in both directions, but predominantly from ‘gibberosa’ into ‘frontosa’. The highest proportion of introgressed ‘gibberosa’ ancestry was present in the southernmost populations of C. frontosa collected near Karilani Island and Cape Kabogo. At the intra-specific level, there was support for between 1 and 3 populations of C. frontosa, whereas the results indicated only a single homogeneous population of C. gibberosa. The presence of different morphs in the lake despite the low levels of heterozygosity suggests that a small number of loci may be involved in the morphological variation and/or that there is a more complex interplay between genetics and the environment in different locations.