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.

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.

Recent Northward Expansion of a Passerine Bird Species, Brownish-Flanked Bush Warbler (Horornis fortipes).

Northward expansions of bird distributions have been commonly observed in the Northern Hemisphere, likely as a result of climate change. The causes and ecological impacts of such range shifts have received extensive attention, but studies on the process of range shifts are still relatively scarce. The Brownish-flanked Bush Warbler (Horornis fortipes) has expanded northward from 35° N to 40° N during the past decade. In this study, we collated 77 records of the species beyond its traditional distribution during the past ten years from citizen science data. Most of the new records were from northeast of its traditional distribution, including the North China Plain, Taihang Mountains, and Taishan Mountain, and a few records from the northern margin of the Qinling Mountains and Qinghai-Tibet Plateau. We concluded that the Brownish-flanked Bush Warbler has bred in this new area in at least six sites. The newly established populations are assumed to belong to the subspecies H. f. davidianus, which can be divided into eastern and western dialect groups based on differences in songs. Song recordings from 10 males from Beijing and its adjacent areas were collected. Bayesian analysis based on the acoustic traits indicated that these males were most likely from the western dialect area, with a posterior probability of 99.975%. Combining topographical data with the habitat preference of the species, we inferred that these individuals spread northeastward from the Qinling Mountains to Taihang Mountains, and further along the Yanshan Mountains. This study is a case study of the distribution expansion of a bird species, which reflects the dynamics of a species in the early stage of its northward expansion.

The formation of the Indo-Pacific montane avifauna.

The processes generating the earth's montane biodiversity remain a matter of debate. Two contrasting hypotheses have been advanced to explain how montane populations form: via direct colonization from other mountains, or, alternatively, via upslope range shifts from adjacent lowland areas. We seek to reconcile these apparently conflicting hypotheses by asking whether a species' ancestral geographic origin determines its mode of mountain colonization. Island-dwelling passerine birds at the faunal crossroads between Eurasia and Australo-Papua provide an ideal study system. We recover the phylogenetic relationships of the region's montane species and reconstruct their ancestral geographic ranges, elevational ranges, and migratory behavior. We also perform genomic population studies of three super-dispersive montane species/clades with broad island distributions. Eurasian-origin species populated archipelagos via direct colonization between mountains. This mode of colonization appears related to ancestral adaptations to cold and seasonal climates, specifically short-distance migration. Australo-Papuan-origin mountain populations, by contrast, evolved from lowland ancestors, and highland distribution mostly precludes their further colonization of island mountains. Our study explains much of the distributional variation within a complex biological system, and provides a synthesis of two seemingly discordant hypotheses for montane community formation.

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.

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.

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.

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.

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.

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.].

Asymmetric introgression reveals the genetic architecture of a plumage trait.

Genome-wide variation in introgression rates across hybrid zones offers a powerful opportunity for studying population differentiation. One poorly understood pattern of introgression is the geographic displacement of a trait implicated in lineage divergence from genome-wide population boundaries. While difficult to interpret, this pattern can facilitate the dissection of trait genetic architecture because traits become uncoupled from their ancestral genomic background. We studied an example of trait displacement generated by the introgression of head plumage coloration from personata to alba subspecies of the white wagtail. A previous study of their hybrid zone in Siberia revealed that the geographic transition in this sexual signal that mediates assortative mating was offset from other traits and genetic markers. Here we show that head plumage is associated with two small genetic regions. Despite having a simple genetic architecture, head plumage inheritance is consistent with partial dominance and epistasis, which could contribute to its asymmetric introgression.

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.

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.

"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.

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.

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.

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.

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.