A Phylogenomic Assessment of Processes Underpinning Convergent Evolution in Open-Habitat Chats.

Insights into the processes underpinning convergent evolution advance our understanding of the contributions of ancestral, introgressed, and novel genetic variation to phenotypic evolution. Phylogenomic analyses characterizing genome-wide gene tree heterogeneity can provide first clues about the extent of ILS and of introgression and thereby into the potential of these processes or (in their absence) the need to invoke novel mutations to underpin convergent evolution. Here, we were interested in understanding the processes involved in convergent evolution in open-habitat chats (wheatears of the genus Oenanthe and their relatives). To this end, based on whole-genome resequencing data from 50 taxa of 44 species, we established the species tree, characterized gene tree heterogeneity, and investigated the footprints of ILS and introgression within the latter. The species tree corroborates the pattern of abundant convergent evolution, especially in wheatears. The high levels of gene tree heterogeneity in wheatears are explained by ILS alone only for 30% of internal branches. For multiple branches with high gene tree heterogeneity, D-statistics and phylogenetic networks identified footprints of introgression. Finally, long branches without extensive ILS between clades sporting similar phenotypes provide suggestive evidence for the role of novel mutations in the evolution of these phenotypes. Together, our results suggest that convergent evolution in open-habitat chats involved diverse processes and highlight that phenotypic diversification is often complex and best depicted as a network of interacting lineages.

Genome-wide evidence supports mitochondrial relationships and pervasive parallel phenotypic evolution in open-habitat chats.

In wheatears and related species ('open-habitat chats'), molecular phylogenetics has led to a comprehensively revised understanding of species relationships and species diversity. Phylogenetic analyses have suggested that, in many cases, phenotypic similarities do not reflect species' relationships, revealing traditionally defined genera as non-monophyletic. This led to the suggestion of pervasive parallel evolution of open-habitat chats' plumage coloration and ecological phenotypes. However, to date, the molecular evidence for the phylogenetic relationships among open-habitat chats is mainly limited to mitochondrial DNA. Here, we assessed whether the mitochondrial relationships are supported by genome-wide data. To this end, we reconstructed the species tree among 14 open-habitat chat taxa using multi-species coalescent analyses based on ~1'300 SNPs. Our results confirm previous ones based chiefly on mitochondrial DNA; notably the paraphyly of the Oenanthe lugens complex and the clustering of individual species formerly placed in the genera Cercomela and Myrmecocichla within Oenanthe. Since several variable morphological and ecological characteristics occur in multiple places across the open-habitat chat phylogeny, our study consolidates the evidence for pervasive parallel evolution in the plumage coloration and ecology of open-habitat chats.

Parallel plumage colour evolution and introgressive hybridization in wheatears.

Genetic and phenotypic mosaics, in which various phenotypes and different genomic regions show discordant patterns of species or population divergence, offer unique opportunities to study the role of ancestral and introgressed genetic variation in phenotypic evolution. Here, we investigated the evolution of discordant phenotypic and genetic divergence in a monophyletic clade of four songbird taxa-pied wheatear (O. pleschanka), Cyprus wheatear (Oenanthe cypriaca), and western and eastern subspecies of black-eared wheatear (O. h. hispanica and O. h. melanoleuca). Phenotypically, black back and neck sides distinguish pied and Cyprus wheatears from the white-backed/necked black-eared wheatears. Meanwhile, mitochondrial variation only distinguishes western black-eared wheatear. In the absence of nuclear genetic data, and given frequent hybridization among eastern black-eared and pied wheatear, it remains unclear whether introgression is responsible for discordance between mitochondrial divergence patterns and phenotypic similarities, or whether plumage coloration evolved in parallel. Multispecies coalescent analyses of about 20,000 SNPs obtained from RAD data mapped to a draft genome assembly resolve the species tree, provide evidence for the parallel evolution of colour phenotypes and establish western and eastern black-eared wheatears as independent taxa that should be recognized as full species. The presence of the entire admixture spectrum in the Iranian hybrid zone and the detection of footprints of introgression from pied into eastern black-eared wheatear beyond the hybrid zone despite strong geographic structure of ancestry proportions furthermore suggest a potential role for introgression in parallel plumage colour evolution. Our results support the importance of standing heterospecific and/or ancestral variation in phenotypic evolution.

An annotated chromosome-scale reference genome for Eastern black-eared wheatear (Oenanthe melanoleuca).

Pervasive convergent evolution and in part high incidences of hybridization distinguish wheatears (songbirds of the genus Oenanthe) as a versatile system to address questions at the forefront of research on the molecular bases of phenotypic and species diversification. To prepare the genomic resources for this venture, we here generated and annotated a chromosome-scale assembly of the Eastern black-eared wheatear (Oenanthe melanoleuca). This species is part of the Oenanthe hispanica complex that is characterized by convergent evolution of plumage coloration and high rates of hybridization. The long-read-based male nuclear genome assembly comprises 1.04 Gb in 32 autosomes, the Z chromosome, and the mitogenome. The assembly is highly contiguous (contig N50, 12.6 Mb; scaffold N50, 70 Mb), with 96% of the genome assembled at the chromosome level and 95.5% benchmarking universal single-copy orthologs (BUSCO) completeness. The nuclear genome was annotated with 18,143 protein-coding genes and 31,333 mRNAs (annotation BUSCO completeness, 98.0%), and about 10% of the genome consists of repetitive DNA. The annotated chromosome-scale reference genome of Eastern black-eared wheatear provides a crucial resource for research into the genomics of adaptation and speciation in an intriguing group of passerines.