Abstract:Instances of parallel phenotypic evolution offer great opportunities to understand the evolutionary processes underlying phenotypic changes. However, confirming parallel phenotypic evolution and studying its causes requires a robust phylogenetic framework. One such example is the "black-and-white wagtails," a group of 5 species in the songbird genus Motacilla: 1 species, Motacilla alba, shows wide intra-specific plumage variation, while the 4r others form 2 pairs of very similar-looking species (M. aguimp + M. samveasnae and M. grandis + M. maderaspatensis, respectively). However, the 2 species in each of these pairs were not recovered as sisters in previous phylogenetic inferences. Their relationships varied depending on the markers used, suggesting that gene tree heterogeneity might have hampered accurate phylogenetic inference. Here, we use whole genome resequencing data to explore the phylogenetic relationships within this group, with a special emphasis on characterizing the extent of gene tree heterogeneity and its underlying causes. We first used multispecies coalescent methods to generate a "complete evidence" phylogenetic hypothesis based on genome-wide variants, while accounting for incomplete lineage sorting (ILS) and introgression. We then investigated the variation in phylogenetic signal across the genome to quantify the extent of discordance across genomic regions and test its underlying causes. We found that wagtail genomes are mosaics of regions supporting variable genealogies, because of ILS and inter-specific introgression. The most common topology across the genome, supporting M. alba and M. aguimp as sister species, appears to be influenced by ancient introgression. Additionally, we inferred another ancient introgression event, between M. alba and M. grandis. By combining results from multiple analyses, we propose a phylogenetic network for the black-and-white wagtails that confirms that similar phenotypes evolved in non-sister lineages, supporting parallel plumage evolution. Furthermore, the inferred reticulations do not connect species with similar plumage coloration, suggesting that introgression does not underlie parallel plumage evolution in this group. Our results demonstrate the importance of investing genome-wide patterns of gene tree heterogeneity to help understand the mechanisms underlying phenotypic evolution. [Gene tree heterogeneity; incomplete lineage sorting; introgression; parallel evolution; phylogenomics; plumage evolution; wagtails.].

Introgression Underlies Phylogenetic Uncertainty But Not Parallel Plumage Evolution in a Recent Songbird Radiation

Discordance between genomic divergence and phenotypic variation in a rapidly evolving avian genus (Motacilla)

Asymmetric introgression reveals the genetic architecture of a plumage trait

A late burst of colour evolution in a radiation of songbirds (Passeriformes: Parulidae) suggests secondary contact drives signal divergence

Genetic basis and evolution of structural color polymorphism in an Australian songbird

Selection on a small genomic region underpins differentiation in multiple color traits between two warbler species

Introgression across narrow contact zones shapes the genomic landscape of phylogenetic variation in an African bird clade

Mutations in different pigmentation genes are associated with parallel melanism in island flycatchers

Mosaic genome evolution in a recent and rapid avian radiation

Demographic histories and genome-wide patterns of divergence in incipient species of shorebirds

Introgression of regulatory alleles and a missense coding mutation drive plumage pattern diversity in the rock pigeon

Complexity of Avian Evolution Revealed by Family-Level Genomes

Comparative Genomics and Transcriptomics of Chrysolophus Provide Insights into the Evolution of Complex Plumage Coloration

Incomplete Lineage Sorting and Introgression in the Diversification of Chinese Spot-Billed Ducks and Mallards

Genome‐wide SNPs Confirm Plumage Polymorphism and Hybridisation Within a Cyornis Flycatcher Species Complex

Gene flow and an anomaly zone complicate phylogenomic inference in a rapidly radiated avian family (Prunellidae)

'Ghost introgression' as a cause of deep mitochondrial divergence in a bird species complex.

Whole-genome analyses resolve early branches in the tree of life of modern birds

Divergent Selection in Low Recombination Regions Shapes the Genomic Islands in Two Incipient Shorebird Species.

Population genomics indicate three different modes of divergence and speciation with gene flow in the green-winged teal duck complex

Comparative analyses identify genomic features potentially involved in the evolution of birds-of-paradise