Genomic islands of speciation harbor genes underlying coloration differences in a pair of Neotropical seedeaters

Incomplete speciation can be leveraged to associate phenotypes with genotypes, thus providing insights into the traits relevant to the reproductive isolation of diverging taxa. We investigate the genetic underpinnings of the phenotypic differences between Sporophila plumbea and S. beltoni. S. beltoni has only recently been described based, most notably, on differences in bill coloration (yellow vs. black in S. plumbea). Both species are indistinguishable through mtDNA or reduced-representation genomic data, and even whole-genome sequencing revealed low genetic differentiation. Demographic reconstructions attribute this genetic homogeneity to gene flow, despite divergence in the order of millions of generations. We found a narrow hybrid zone in southern Brazil where genetically, yet not phenotypically, admixed individuals appear to be prevalent. Despite the overall low genetic differentiation, we identified three narrow peaks along the genome with highly differentiated SNPs. These regions harbor six genes, one of which is involved in pigmentation (EDN3) and is a candidate for controlling bill color. Within the outlier peaks we found signatures of resistance to gene flow, as expected for islands of speciation. Our study shows how genes related to coloration traits are likely involved in generating prezygotic isolation and establishing species boundaries early in speciation.