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

BACKGROUND: Resolving the phylogeny of rapidly radiating lineages presents a challenge when building the Tree of Life. An Old World avian family Prunellidae (Accentors) comprises twelve species that rapidly diversified at the Pliocene-Pleistocene boundary. RESULTS: Here we investigate the phylogenetic relationships of all species of Prunellidae using a chromosome-level de novo assembly of Prunella strophiata and 36 high-coverage resequenced genomes. We use homologous alignments of thousands of exonic and intronic loci to build the coalescent and concatenated phylogenies and recover four different species trees. Topology tests show a large degree of gene tree-species tree discordance but only 40-54% of intronic gene trees and 36-75% of exonic genic trees can be explained by incomplete lineage sorting and gene tree estimation errors. Estimated branch lengths for three successive internal branches in the inferred species trees suggest the existence of an empirical anomaly zone. The most common topology recovered for species in this anomaly zone was not similar to any coalescent or concatenated inference phylogenies, suggesting presence of anomalous gene trees. However, this interpretation is complicated by the presence of gene flow because extensive introgression was detected among these species. When exploring tree topology distributions, introgression, and regional variation in recombination rate, we find that many autosomal regions contain signatures of introgression and thus may mislead phylogenetic inference. Conversely, the phylogenetic signal is concentrated to regions with low-recombination rate, such as the Z chromosome, which are also more resistant to interspecific introgression. CONCLUSIONS: Collectively, our results suggest that phylogenomic inference should consider the underlying genomic architecture to maximize the consistency of phylogenomic signal.

Limited Song Mixing Without Genomic Gene Flow in a Contact Zone Between Two Songbird Species.

Song is considered to play an important role in the maintenance of prezygotic reproductive isolation between closely related songbird species. Therefore, song mixing in a contact zone between closely related species is often considered as evidence of hybridization. The Sichuan Leaf Warbler Phylloscopus forresti and the Gansu Leaf Warbler Phylloscopus kansuensis, which diverged 2 million years ago, have formed a contact zone in the south of the Gansu Province of China, where mixed songs have been observed. In this study, we investigated the potential causes and consequences of song mixing by integrating bioacoustic, morphological, mitochondrial, and genomic data with field ecological observations. We found that the two species display no apparent morphological differences, whereas their songs differ dramatically. We demonstrated that ∼11% of the males in the contact zone sang mixed songs. Two males singing mixed song were genotyped, and both were found to be P. kansuensis. Despite the presence of mixed singers, population genomic analyses detected no signs of recent gene flow between the two species, although two possible cases of mitochondrial introgression were identified. We conclude that the rather limited song mixing does not lead to, or result from, hybridization, and hence does not result in the breakdown of reproductive barriers between these cryptic species.

When colors mislead: Genomics and bioacoustics prompt re-classification of Asian flycatcher radiation (Aves: Niltavinae).

Traditional classification of many animals, including birds, has been highly dependent on external morphological characters like plumage coloration. However, both bioacoustics and genetic or genomic data have revolutionized our understanding of the relationships of certain lineages and led to sweeping taxonomic re-organizations. In this study, we present a case of erroneous delimitation of genus boundaries in the species-rich flycatcher subfamily Niltavinae. Genera within this subfamily have historically been delineated based on blue versus brown male body plumage until recent studies based on a few mitochondrial and nuclear loci unearthed several cases of generic misclassification. Here we use extensive bioacoustic data from 43 species and genomic data from 28 species for a fundamental reclassification of species in the Niltavinae. Our study reveals that song is an important trait to classify these birds even at the genus level, whereas plumage traits exhibit ample convergence and have led to numerous historic misattributions. Our taxonomic re-organization leads to new biogeographic limits of major genera, such that the genus Cyornis now only extends as far east as the islands of Sulawesi, Sula, and Banggai, whereas Eumyias is redefined to extend far beyond Wallace's Line to the islands of Seram and Timor. Our conclusions advise against an over-reliance on morphological traits and underscore the importance of integrative datasets.

Introgression underlies phylogenetic uncertainty but not parallel plumage evolution in a recent songbird radiation.

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 five species in the songbird genus Motacilla: one species, Motacilla alba, shows wide intra-specific plumage variation, while the four others form two pairs of very similar-looking species (M. aguimp + M. samveasnae and M. grandis + M. maderaspatensis, respectively). However, the two 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 investigation of genome-wide patterns of gene tree heterogeneity to help understanding the mechanisms underlying phenotypic evolution.

Phylogenetic conflict between species tree and maternally inherited gene trees in a clade of Emberiza buntings (Aves: Emberizidae).

Different genomic regions may reflect conflicting phylogenetic topologies primarily due to incomplete lineage sorting and/or gene flow. Genomic data are necessary to reconstruct the true species tree and explore potential causes of phylogenetic conflict. Here, we investigate the phylogenetic relationships of four Emberiza species (Aves: Emberizidae) and discuss the potential causes of the observed mitochondrial non-monophyly of Emberiza godlewskii (Godlewski's bunting) using phylogenomic analyses based on whole genome resequencing data from 41 birds. Analyses based on both the whole mitochondrial genome and ~39 kilobases from the non-recombining W chromosome reveal sister relationships between each the northern and southern populations of E. godlewskii with E. cioides and E. cia, respectively. In contrast, the monophyly of E. godlewskii is reflected by the phylogenetic signal of autosomal and Z chromosomal sequence data as well as demographic inference analyses, which - in combination - support the following tree topology: (((E. godlewskii, E. cia), E. cioides), E. jankowskii). Using D-statistics, we detected multiple gene flow events among different lineages, indicating pervasive introgressive hybridization within this clade. Introgression from an unsampled lineage that is sister to E. cioides or introgression from an unsampled mitochondrial + W chromosomal lineage of E. cioides into northern E. godlewskii may explain the phylogenetic conflict between the species tree estimated from genome-wide data versus mtDNA/W tree topologies. These results underscore the importance of using genomic data for phylogenetic reconstruction and species delimitation.

Genomic and phenotypic changes associated with alterations of migratory behaviour in a songbird.

The seasonal migration of birds is a fascinating natural wonder. Avian migratory behaviour changes are common and are probably a polygenic process as avian migration is governed by multiple correlated components with a variable genetic basis. However, the genetic and phenotypic changes involving migration changes are poorly studied. Using one annotated near-chromosomal level de novo genome assembly, 50 resequenced genomes, hundreds of morphometric data and species distribution information, we investigated population structure and genomic and phenotypic differences associated with differences in migratory behaviour in a songbird species, Yellow-throated Bunting Emberiza elegans (Aves: Emberizidae). Population genomic analyses reveal extensive gene flow between the southern resident and the northern migratory populations of this species. The hand-wing index is significantly lower in the resident populations than in the migratory populations, indicating reduced flight efficiency of the resident populations. Here, we discuss the possibility that nonmigratory populations may have originated from migratory populations though migration loss. We further infer that the alterations of genes related to energy metabolism, nervous system and circadian rhythm may have played major roles in regulating migration change. Our study sheds light on phenotypic and polygenic changes involving migration change.

A near-complete and time-calibrated phylogeny of the Old World flycatchers, robins and chats (Aves, Muscicapidae).

The Old World flycatchers, robins and chats (Aves, Muscicapidae) are a diverse songbird family with over three hundred species. Despite continuous efforts over the past two decades, there is still no comprehensive and well-resolved species-level phylogeny for Muscicapidae. Here we present a supermatrix phylogeny that includes all 50 currently recognized genera and ca. 92% of all the species, built using data from up to 15 mitochondrial and 13 nuclear loci. In addition to assembling nucleotide sequences available in public databases, we also extracted sequences from the genome assemblies and raw sequencing reads from GenBank and included a few unpublished sequences. Our analyses resolved the phylogenetic position for several previously unsampled taxa, for example, the Grand Comoro Flycatcher Humblotia flavirostris, the Collared Palm Thrush Cichladusa arquata, and the Taiwan Whistling-Thrush Myophonus insularis, etc. We also provide taxonomic recommendations for genera that exhibit paraphyly or polyphyly. Our results suggest that Muscicapidae diverged from Turdidae (thrushes and allies) in the early Miocene, and the most recent common ancestors for the four subfamilies (Muscicapinae, Niltavinae, Cossyphinae and Saxicolinae) all arose around the middle Miocene.

Cryptic species in a colorful genus: Integrative taxonomy of the bush robins (Aves, Muscicapidae, Tarsiger) suggests two overlooked species.

Several cryptic avian species have been validated by recent integrative taxonomic efforts in the Sino-Himalayan mountains, indicating that avian diversity in this global biodiversity hotspot may be underestimated. In the present study, we investigated species limits in the genus Tarsiger, the bush robins, a group of montane forest specialists with high species richness in the Sino-Himalayan region. Based on comprehensive sampling of all 11 subspecies of the six currently recognized species, we applied an integrative taxonomic approach by combining multilocus, acoustic, plumage and morphometric analyses. Our results reveal that the isolated north-central Chinese populations of Tarsiger cyanurus, described as the subspecies albocoeruleus but usually considered invalid, is distinctive in genetics and vocalisation, but only marginally differentiated in morphology. We also found the Taiwan endemic T. indicus formosanus to be distinctive in genetics, song and morphology from T. i. indicus and T. i. yunnanensis of the Sino-Himalayan mountains. Moreover, Bayesian species delimitation using BPP suggests that both albocoeruleus and formosanus merit full species status. We propose their treatment as 'Qilian Bluetail' T. albocoeruleus and 'Taiwan Bush Robin' T. formosanus, respectively.

Most Genomic Loci Misrepresent the Phylogeny of an Avian Radiation Because of Ancient Gene Flow.

Phylogenetic trees based on genome-wide sequence data may not always represent the true evolutionary history for a variety of reasons. One process that can lead to incorrect reconstruction of species phylogenies is gene flow, especially if interspecific gene flow has affected large parts of the genome. We investigated phylogenetic relationships within a clade comprising eight species of passerine birds (Phylloscopidae, Phylloscopus, leaf warblers) using one de novo genome assembly and 78 resequenced genomes. On the basis of hypothesis-exclusion trials based on D-statistics, phylogenetic network analysis, and demographic inference analysis, we identified ancient gene flow affecting large parts of the genome between one species and the ancestral lineage of a sister species pair. This ancient gene flow consistently caused erroneous reconstruction of the phylogeny when using large amounts of genome-wide sequence data. In contrast, the true relationships were captured when smaller parts of the genome were analyzed, showing that the "winner-takes-all democratic majority tree" is not necessarily the true species tree. Under this condition, smaller amounts of data may sometimes avoid the effects of gene flow due to stochastic sampling, as hidden reticulation histories are more likely to emerge from the use of larger data sets, especially whole-genome data sets. In addition, we also found that genomic regions affected by ancient gene flow generally exhibited higher genomic differentiation but a lower recombination rate and nucleotide diversity. Our study highlights the importance of considering reticulation in phylogenetic reconstructions in the genomic era.[Bifurcation; introgression; recombination; reticulation; Phylloscopus.].

Earth history and the passerine superradiation.

Avian diversification has been influenced by global climate change, plate tectonic movements, and mass extinction events. However, the impact of these factors on the diversification of the hyperdiverse perching birds (passerines) is unclear because family level relationships are unresolved and the timing of splitting events among lineages is uncertain. We analyzed DNA data from 4,060 nuclear loci and 137 passerine families using concatenation and coalescent approaches to infer a comprehensive phylogenetic hypothesis that clarifies relationships among all passerine families. Then, we calibrated this phylogeny using 13 fossils to examine the effects of different events in Earth history on the timing and rate of passerine diversification. Our analyses reconcile passerine diversification with the fossil and geological records; suggest that passerines originated on the Australian landmass ∼47 Ma; and show that subsequent dispersal and diversification of passerines was affected by a number of climatological and geological events, such as Oligocene glaciation and inundation of the New Zealand landmass. Although passerine diversification rates fluctuated throughout the Cenozoic, we find no link between the rate of passerine diversification and Cenozoic global temperature, and our analyses show that the increases in passerine diversification rate we observe are disconnected from the colonization of new continents. Taken together, these results suggest more complex mechanisms than temperature change or ecological opportunity have controlled macroscale patterns of passerine speciation.

Multiple species delimitation approaches applied to the avian lark genus Alaudala.

Species delimitation has advanced from a purely phenotypic exercise to a branch of science that integrates multiple sources of data to identify independently evolving lineages that can be treated as species. We here test species limits in the avian Lesser Short-toed Lark Alaudala rufesens-Sand Lark A. raytal complex, which has an intricate taxonomic history, ranging from a single to three recognised species, with different inclusiveness in different treatments. Our integrative taxonomic approach is based on a combination of DNA sequences, plumage, biometrics, songs, song-flights, geographical distributions, habitat, and bioclimatic data, and using various methods including a species delimitation program (STACEY) based on the multispecies coalescent model. We propose that four species should be recognised: Lesser Short-toed Lark A. rufescens (sensu stricto), Heine's Short-toed Lark A. heinei, Asian Short-toed Lark A. cheleensis and Sand Lark A. raytal. There is also some evidence suggesting lineage separation within A. cheleensis and A. raytal, but additional data are required to evaluate this. The species delimitation based on STACEY agrees well with the non-genetic data. Although computer-based species delimitation programs can be useful in identifying independently evolving lineages, we stress that whenever possible, species hypotheses proposed by these programs should be tested by independent, non-genetic data. Our results highlight the difficulty and subjectivity of delimiting lineages and species, especially at early stages in the speciation process.

"Ghost Introgression" As a Cause of Deep Mitochondrial Divergence in a Bird Species Complex.

In the absence of nuclear-genomic differentiation between two populations, deep mitochondrial divergence (DMD) is a form of mito-nuclear discordance. Such instances of DMD are rare and might variably be explained by unusual cases of female-linked selection, by male-biased dispersal, by "speciation reversal" or by mitochondrial capture through genetic introgression. Here, we analyze DMD in an Asian Phylloscopus leaf warbler (Aves: Phylloscopidae) complex. Bioacoustic, morphological, and genomic data demonstrate close similarity between the taxa affinis and occisinensis, even though DMD previously led to their classification as two distinct species. Using population genomic and comparative genomic methods on 45 whole genomes, including historical reconstructions of effective population size, genomic peaks of differentiation and genomic linkage, we infer that the form affinis is likely the product of a westward expansion in which it replaced a now-extinct congener that was the donor of its mtDNA and small portions of its nuclear genome. This study provides strong evidence of "ghost introgression" as the cause of DMD, and we suggest that "ghost introgression" may be a widely overlooked phenomenon in nature.

A comprehensive phylogeny and taxonomic evaluation of the waxbills (Aves: Estrildidae).

We present a revised taxonomy of Estrildidae based on the first time-calibrated phylogeny of the family Estrildidae estimated from a data set including the majority of the species, and all genera except the monospecific Paludipasser, using two mitochondrial and five nuclear markers. We find that most differences in current taxonomy reflect alternative opinions among authors regarding inclusiveness of genera, which are usually not in conflict with the phylogeny. The most notable exception is the current circumscriptions of the genera Neochmia, Nesocharis and Taeniopygia, which are incompatible with the phylogeny. Estrildidae is subdivided into six well supported subclades, which we propose be recognized as the subfamilies Amandavinae, Erythrurinae, Estrildinae, Lagonostictinae, Lonchurinae and Poephilinae.

Niche filling slows the diversification of Himalayan songbirds.

Speciation generally involves a three-step process--range expansion, range fragmentation and the development of reproductive isolation between spatially separated populations. Speciation relies on cycling through these three steps and each may limit the rate at which new species form. We estimate phylogenetic relationships among all Himalayan songbirds to ask whether the development of reproductive isolation and ecological competition, both factors that limit range expansions, set an ultimate limit on speciation. Based on a phylogeny for all 358 species distributed along the eastern elevational gradient, here we show that body size and shape differences evolved early in the radiation, with the elevational band occupied by a species evolving later. These results are consistent with competition for niche space limiting species accumulation. Even the elevation dimension seems to be approaching ecological saturation, because the closest relatives both inside the assemblage and elsewhere in the Himalayas are on average separated by more than five million years, which is longer than it generally takes for reproductive isolation to be completed; also, elevational distributions are well explained by resource availability, notably the abundance of arthropods, and not by differences in diversification rates in different elevational zones. Our results imply that speciation rate is ultimately set by niche filling (that is, ecological competition for resources), rather than by the rate of acquisition of reproductive isolation.

From the Himalayas to a continental Island: Integrative species delimitation in the Brownish-flanked Bush Warbler Horornis fortipes complex.

As species serve as basic units of study in many fields of biology, assessments of species limits are fundamental for such studies. Here, we used a multilocus dataset and different coalescent-based methods to analyze species delimitation and phylogenetic relationships in the Brownish-flanked Bush Warbler Horornis fortipes complex, which is widespread in the Sino-Himalayan region. We also examined the vocal and morphometric divergence within this complex. Our genetic results suggested that Horornis fortipes is composed of at least three independently evolving lineages, which diverged 1.1-1.8 million years ago. However, these lineages have hardly diverged in song or morphometrics and only very slightly in plumage. Our result indicate that there are three incipient species in Horonis fortipes complex diverged in central Himalayas and Hengduan Mountains, but not between the continent and Taiwan island.

Near-complete phylogeny and taxonomic revision of the world's babblers (Aves: Passeriformes).

The babblers are a diverse group of passerine birds comprising 452 species. The group was long regarded as a "scrap basket" in taxonomic classification schemes. Although several studies have assessed the phylogenetic relationships for subsets of babblers during the past two decades, a comprehensive phylogeny of this group has been lacking. In this study, we used five mitochondrial and seven nuclear loci to generate a dated phylogeny for babblers. This phylogeny includes 402 species (ca. 89% of the overall clade) from 75 genera (97%) and all five currently recognized families, providing a robust basis for taxonomic revision. Our phylogeny supports seven major clades and reveals several non-monophyletic genera. Divergence time estimates indicate that the seven major clades diverged around the same time (18-20 million years ago, Ma) in the early Miocene. We use the phylogeny in a consistent way to propose a new taxonomy, with seven families and 64 genera of babblers, and a new linear sequence of names.

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

Generally, genotypes and phenotypes are expected to be spatially congruent; however, in widespread species complexes with few barriers to dispersal, multiple contact zones, and limited reproductive isolation, discordance between phenotypes and phylogeographic groups is more probable. Wagtails (Motacilla) are a genus of birds with striking plumage pattern variation across the Old World. Up to 13 subspecies are recognized within a single species, yet previous studies using mitochondrial DNA have supported polyphyletic phylogeographic groups that are inconsistent with subspecies plumage characteristics. In this study, we investigate the link between phenotypes and genotype by taking a phylogenetic approach. We use genome-wide SNPs, nuclear introns, and mitochondrial DNA to estimate population structure, isolation by distance, and species relationships. Together, our genetic sampling includes complete species-level sampling and comprehensive coverage of the three most phenotypically diverse Palearctic species. Our study provides strong evidence for species-level patterns of differentiation, however population-level differentiation is less pronounced. SNPs provide a robust estimate of species-level relationships, which are mostly corroborated by a combined analysis of mtDNA and nuclear introns (the first time-calibrated species tree for the genus). However, the mtDNA tree is strongly incongruent and is considered to misrepresent the species phylogeny. The extant wagtail lineages originated during the Pliocene and the Eurasian lineage underwent rapid diversification during the Pleistocene. Three of four widespread Eurasian species exhibit an east-west divide that contradicts both subspecies taxonomy and phenotypic variation. Indeed, SNPs fail to distinguish between phenotypically distinct subspecies within the M. alba and M. flava complexes, and instead support geographical regions, each of which is home to two or more different looking subspecies. This is a major step towards our understanding of wagtail phylogeny compared to previous analyses of fewer species and considerably less sequence data.

Comprehensive molecular phylogeny of the grassbirds and allies (Locustellidae) reveals extensive non-monophyly of traditional genera, and a proposal for a new classification.

The widespread Old World avian family Locustellidae ('grassbirds and allies') comprises 62 extant species in 11 genera. In the present study, we used one mitochondrial and, for most species, four nuclear loci to infer the phylogeny of this family. We analysed 59 species, including the five previously unsampled genera plus two genera that had not before been analysed in a densely sampled dataset. This study revealed extensive disagreement with current taxonomy; the genera Bradypterus, Locustella, Megalurus, Megalurulus and Schoenicola were all found to be non-monophyletic. Non-monophyly was particularly pronounced for Megalurus, which was widely scattered across the tree. Three of the five monotypic genera (Amphilais, Buettikoferella and Malia) were nested within other genera; one monotypic genus (Chaetornis) formed a clade with one of the two species of Schoenicola; whereas the position of the fifth monotypic genus (Elaphrornis) was unresolved. Robsonius was confirmed as sister to the other genera. We propose a phylogenetically informed revision of genus-level taxonomy, including one new generic name. Finally, we highlight several non-monophyletic species complexes and deep intra-species divergences that point to conflict in taxonomy and suggest an underestimation of current species diversity in this group.

Complete species-level phylogeny of the leaf warbler (Aves: Phylloscopidae) radiation.

The leaf warbler radiation (Aves: Phylloscopidae) has undergone a c. 50% increase in the number of recognised species over the last three decades, mainly as a result of analyses of vocalisations and DNA. Using a multilocus dataset for all of the species in this family, and multispecies coalescent-based as well as concatenation methods, we provide the first complete species-level phylogeny for this important group, as well as an estimate of the timing of diversification. The most recent common ancestor for the family was dated at 11.7 million years ago (mya) (95% highest posterior density 9.8-13.7 mya), and divergence times between sister species ranged from 0.5 mya (0.3-0.8 mya) to 6.1 mya (4.8-7.5 mya). Based on our results, we support synonymising Seicercus with Phylloscopus, which results in a monogeneric Phylloscopidae. We discuss the pros and cons of this treatment,and we argue againstproliferation of taxonomic names,and conclude that a large monogeneric Phylloscopidae leads to the fewest taxonomic changes compared to traditional classifications. We briefly discuss morphological evolution in the light of the phylogeny. The time calibrated phylogeny is a major improvement compared to previous studies based on a smaller number of species and loci and can provide a basis for future studies of other aspects of phylloscopid evolution.

Genomic differentiation and patterns of gene flow between two long-tailed tit species (Aegithalos).

Patterns of heterogeneous genomic differentiation have been well documented between closely related species, with some highly differentiated genomic regions ("genomic differentiation islands") spread throughout the genome. Differential levels of gene flow are proposed to account for this pattern, as genomic differentiation islands are suggested to be resistant to gene flow. Recent studies have also suggested that genomic differentiation islands could be explained by linked selection acting on genomic regions with low recombination rates. Here, we investigate genomic differentiation and gene-flow patterns for autosomes using RAD-seq data between two closely related species of long-tailed tits (Aegithalos bonvaloti and A. fuliginosus) in both allopatric and contact zone populations. The results confirm recent or ongoing gene flow between these two species. However, there is little evidence that the genomic regions that were found to be highly differentiated between the contact zone populations are resistant to gene flow, suggesting that differential levels of gene flow is not the cause of the heterogeneous genomic differentiation. Linked selection may be the cause of genomic differentiation islands between the allopatric populations with no or very limited gene flow, but this could not account for the heterogeneous genomic differentiation between the contact zone populations, which show evidence of recent or ongoing gene flow.