Patterns of genetic, phenotypic, and acoustic variation across a chiffchaff (Phylloscopus collybita abietinus/tristis) hybrid zone

Trait Variation and Spatiotemporal Dynamics across Avian Secondary Contact Zones

A secondary contact zone (SCZ) is an area where incipient species or divergent populations may meet, mate, and hybridize. Due to the diverse patterns of interspecific hybridization, SCZs function as field labs for illuminating the on-going evolutionary processes of speciation and the establishment of reproductive isolation. Interspecific hybridization is widely present in avian populations, making them an ideal system for SCZ studies. This review exhaustively summarizes the variations in unique traits within avian SCZs (vocalization, plumage, beak, and migratory traits) and the various movement patterns of SCZs observed in previous publications. It also highlights several potential future research directions in the genomic era, such as the relationship between phenotypic and genomic differentiation in SCZs, the genomic basis of trait differentiation, SCZs shared by multiple species, and accurate predictive models for forecasting future movements under climate change and human disturbances. This review aims to provide a more comprehensive understanding of speciation processes and offers a theoretical foundation for species conservation.

Multi-character approach reveals a new mangrove population of the Yellow Warbler complex, Setophaga petechia, on Cozumel Island, Mexico

The Setophaga petechia complex includes 43 subspecies distributed within the new world, of which some are migratory and others are resident, with only two resident subspecies in the Mexican Caribbean: Setophaga petechia bryanti a mangrove subspecies belonging to the erithachorides group resident on the mainland of the Yucatan Peninsula and Setophaga petechia rufivertex endemic to Cozumel Island and belonging to the petechia group. Recently, a new population of individuals presenting intermediate phenotypic traits and living in mangrove ecosystems was discovered and reported for Cozumel Island. In this study, we used a multi-character approach including genetic (five ISSR genetic markers), morphometric (eight traits), phenotypic (four characteristics of males), and acoustic dataset (11 parameters) to understand the process of differentiation and the status of these new island individuals in relation to the two well-established subspecies using a total of 60 individuals (20 for each group). Through multivariate analyses based on different dataset used in our study, we show how the new population is related to the endemic island subspecies, S. p. rufivertex and to the mainland subspecies, S. p. bryanti while demonstrating finite differences. We conclude that the new population of S. petechia on Cozumel Island is a well-established population with high level of differentiation.

Bioacoustic differentiation of calls in the chiffchaff complex

The chiffchaff complex is a group of common forest bird species, notorious for the number of cryptic taxa recently discovered, being a great example of speciation in action. Vocalizations have been crucial to unveil its hidden diversity. In this study we quantitatively analyze the acoustic characteristics of their calls with permutational analysis of variance, canonical variate analysis and a self-organizing map, to determine their variability and differences. We related these differences with the geographical and genetic distances between taxonomic groups, by means of Pearson correlations. We used recordings from Xeno-canto, an open database of bird vocalizations. Inter-taxa distances based on call traits were broadly consistent with geographic distances but not correlated with genetic distances. The Iberian Chiffchaff (Phylloscopus ibericus), presumably the most ancient lineage, was the most central in the variation space, while the Siberian Chiffchaff (P. collybita tristis) was the most peripheric and also very uniform, in contrast with the Canarian Chiffchaff (P. canariensis) highly variable, as expected by the "character release hypothesis" on islands. Calls proved to be an excellent tool, especially amenable for non-biased mathematical analyses which, combined with the wide availability of records in Xeno-canto, greatly facilitates the widespread use of this methodology in a wide range of species and geographical areas.

Introgression at the emerging secondary contact zone of magpie Pica pica subspecies (Aves: Corvidae): integrating data on nuclear and mitochondrial markers, vocalizations, and field observations

Zones of secondary contact provide a good opportunity to investigate the origin and dynamics of reproductive isolation between related populations. We analyzed genetic and phenotypic patterns and gene flow between two subspecies of the Eurasian magpie Pica pica s.l. which recently came into contact after presumably long periods of isolation. We describe the distribution of subspecies in a young contact zone at Argun’ river basin in southern Siberia where populations occur in parapatry and an older hybrid population in eastern Mongolia. Based on genome-wide SNP data, we analyzed patterns and strength of gene flow between the subspecies. Our results indicate occasional hybridization with backcrossing and asymmetric introgression along a wide range in Transbaikalia and locally in eastern Mongolia. Males of P. p. jankowskii apparently exhibit higher dispersal ability towards the west compared to P. p. leucoptera (towards the east). The former occasionally migrates to eastern Mongolia and Transbaikalia where introgression of nuclear, but not mitochondrial DNA was evident. Bioacoustic investigations showed differences between the subspecies in speed and structure of vocalization. We discovered intermediate calls of hybrid magpies and bilingual birds alternating calls that are typical for the two taxa. Furthermore, we found dramatically decreased reproductive success in hybridogeneous populations. By complementing our results with established phylogeographic patterns of P. pica s.l. based on a mitochondrial marker sequence, and considering indications of sterility of hybrids in the contact zone, we propose to elevate the two corresponding subspecies to species level: P. pica for the western form and P. serica for the eastern form.

Asymmetric allelic introgression across a hybrid zone of the coal tit (Periparus ater) in the central Himalayas

In the Himalayas, a number of secondary contact zones have been described for vicariant vertebrate taxa. However, analyses of genetic divergence and admixture are missing for most of these examples. In this study, we provide a population genetic analysis for the coal tit (Periparus ater) hybrid zone in Nepal. Intermediate phenotypes between the distinctive western "spot-winged tit" (P. a. melanolophus) and Eastern Himalayan coal tits (P. a. aemodius) occur across a narrow range of <100 km in western Nepal. As a peculiarity, another distinctive cinnamon-bellied form is known from a single population so far. Genetic admixture of western and eastern mitochondrial lineages was restricted to the narrow zone of phenotypically intermediate populations. The cline width was estimated 46 km only with a center close to the population of the cinnamon-bellied phenotype. In contrast, allelic introgression of microsatellite loci was asymmetrical from eastern P. a. aemodius into far western populations of phenotypic P. a. melanolophus but not vice versa. Accordingly, the microsatellite cline was about 3.7 times wider than the mitochondrial one.

Landscape resistance constrains hybridization across contact zones in a reproductively and morphologically polymorphic salamander

Explicitly accounting for phenotypic differentiation together with environmental heterogeneity is crucial to understand the evolutionary dynamics in hybrid zones. Species showing intra-specific variation in phenotypic traits that meet across environmentally heterogeneous regions constitute excellent natural settings to study the role of phenotypic differentiation and environmental factors in shaping the spatial extent and patterns of admixture in hybrid zones. We studied three environmentally distinct contact zones where morphologically and reproductively divergent subspecies of Salamandra salamandra co-occur: the pueriparous S. s. bernardezi that is mostly parapatric to its three larviparous subspecies neighbours. We used a landscape genetics framework to: (i) characterise the spatial location and extent of each contact zone; (ii) assess patterns of introgression and hybridization between subspecies pairs; and (iii) examine the role of environmental heterogeneity in the evolutionary dynamics of hybrid zones. We found high levels of introgression between parity modes, and between distinct phenotypes, thus demonstrating the evolution to pueriparity alone or morphological differentiation do not lead to reproductive isolation between these highly divergent S. salamandra morphotypes. However, we detected substantial variation in patterns of hybridization across contact zones, being lower in the contact zone located on a topographically complex area. We highlight the importance of accounting for spatial environmental heterogeneity when studying evolutionary dynamics of hybrid zones.

Geographic patterns of mtDNA and Z-linked sequence variation in the Common Chiffchaff and the 'chiffchaff complex'

The Common Chiffchaff Phylloscopus collybita is an abundant, polytypic Palearctic bird. Validity of some of its subspecies is controversial and birds from some parts of the species range remain unclassified taxonomically. The relationships among populations from different geographic areas have not been sufficiently explored with molecular data. In this study we analyzed the relationships among the four species in the ‘chiffchaff complex’ (Common Chiffchaff, Iberian Chiffchaff P. ibericus, Canary Islands Chiffchaff P. canariensis and Mountain Chiffchaff P. sindianus), and the patterns of intraspecific geographic variation in the mtDNA ND2 gene and intron 9 of the Z-linked aconitase gene (ACO1I9) across the Common Chiffchaff range, including a recently discovered population breeding on Mt. Hermon (Anti-Lebanon mountains). Our data supported the monophyly of the chiffchaff complex and its current systematics at the species level. Within the Common Chiffchaff, the Siberian race P. c. tristis was the most differentiated subspecies and may represent a separate or incipient species. Other Common Chiffchaff subspecies also were differentiated in their mtDNA, however, lineages of neighboring subspecies formed wide zones of introgression. The Mt. Hermon population was of mixed genetic origin but contained some birds with novel unique lineage that could not be assigned to known subspecies. All Common Chiffchaff lineages diverged at the end of the Ionian stage of Pleistocene. Lineage sorting of ACO1I9 alleles was not as complete as that of mtDNA. Chiffchaff species were mostly distinct at ACO1I9, except the Common and Canary Islands Chiffchaffs that shared multiple alleles. An AMOVA identified geographic structure in Common Chiffchaff ACO1I9 variation that was broadly consistent with that of mtDNA ND2 gene. The genetic and other data suggest the chiffchaff complex to be a group of evolutionarily young taxa that represent a paradigm of ‘species evolution in action’ from intergrading subspecies through to apparently complete biological speciation.

Heterogeneous Patterns of Genetic Diversity and Differentiation in European and Siberian Chiffchaff (Phylloscopus collybita abietinus/P. tristis)

Identification of candidate genes for trait variation in diverging lineages and characterization of mechanistic underpinnings of genome differentiation are key steps toward understanding the processes underlying the formation of new species. Hybrid zones provide a valuable resource for such investigations, since they allow us to study how genomes evolve as species exchange genetic material and to associate particular genetic regions with phenotypic traits of interest. Here, we use whole-genome resequencing of both allopatric and hybridizing populations of the European (Phylloscopus collybita abietinus) and the Siberian chiffchaff (P. tristis)—two recently diverged species which differ in morphology, plumage, song, habitat, and migration—to quantify the regional variation in genome-wide genetic diversity and differentiation, and to identify candidate regions for trait variation. We find that the levels of diversity, differentiation, and divergence are highly heterogeneous, with significantly reduced global differentiation, and more pronounced differentiation peaks in sympatry than in allopatry. This pattern is consistent with regional differences in effective population size and recurrent background selection or selective sweeps reducing the genetic diversity in specific regions prior to lineage divergence, but the data also suggest that postdivergence selection has resulted in increased differentiation and fixed differences in specific regions. We find that hybridization and backcrossing is common in sympatry, and that phenotype is a poor predictor of the genomic composition of sympatric birds. The combination of a differentiation scan approach with identification of fixed differences pinpoint a handful of candidate regions that might be important for trait variation between the two species.

Patterns of genetic, phenotypic, and acoustic variation across a chiffchaff (Phylloscopus collybita abietinus/tristis) hybrid zone

Characterizing patterns of evolution of genetic and phenotypic divergence between incipient species is essential to understand how evolution of reproductive isolation proceeds. Hybrid zones are excellent for studying such processes, as they provide opportunities to assess trait variation in individuals with mixed genetic background and to quantify gene flow across different genomic regions. Here, we combine plumage, song, mtDNA and whole‐genome sequence data and analyze variation across a sympatric zone between the European and the Siberian chiffchaff (Phylloscopus collybita abietinus/tristis) to study how gene exchange between the lineages affects trait variation. Our results show that chiffchaff within the sympatric region show more extensive trait variation than allopatric birds, with a large proportion of individuals exhibiting intermediate phenotypic characters. The genomic differentiation between the subspecies is lower in sympatry than in allopatry and sympatric birds have a mix of genetic ancestry indicating extensive ongoing and past gene flow. Patterns of phenotypic and genetic variation also vary between regions within the hybrid zone, potentially reflecting differences in population densities, age of secondary contact, or differences in mate recognition or mate preference. The genomic data support the presence of two distinct genetic clades corresponding to allopatric abietinus and tristis and that genetic admixture is the force underlying trait variation in the sympatric region—the previously described subspecies (“fulvescens”) from the region is therefore not likely a distinct taxon. In addition, we conclude that subspecies identification based on appearance is uncertain as an individual with an apparently distinct phenotype can have a considerable proportion of the genome composed of mixed alleles, or even a major part of the genome introgressed from the other subspecies. Our results provide insights into the dynamics of admixture across subspecies boundaries and have implications for understanding speciation processes and for the identification of specific chiffchaff individuals based on phenotypic characters.