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.

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.

Conservation genomics and systematics of a near-extinct island radiation.

Conservation benefits from incorporating genomics to explore the impacts of population declines, inbreeding, loss of genetic variation and hybridization. Here we use the near-extinct Mariana Islands reedwarbler radiation to showcase how ancient DNA approaches can allow insights into the population dynamics of extinct species and threatened populations for which historical museum specimens or material with low DNA yield (e.g., scats, feathers) are the only sources for DNA. Despite their having paraphyletic mitochondrial DNA (mtDNA), nuclear single nucleotide polymorphisms (SNPs) support the distinctiveness of critically endangered Acrocephalus hiwae and the other three species in the radiation that went extinct between the 1960s and 1990s. Two extinct species, A. yamashinae and A. luscinius, were deeply divergent from each other and from a third less differentiated lineage containing A. hiwae and extinct A. nijoi. Both mtDNA and SNPs suggest that the two isolated populations of A. hiwae from Saipan and Alamagan Islands are sufficiently distinct to warrant subspecies recognition and separate conservation management. We detected no significant differences in genetic diversity or inbreeding between Saipan and Alamagan, nor strong signatures of geographical structuring within either island. However, the implications of possible signatures of inbreeding in both Saipan and Alamagan, and long-term population declines in A. hiwae that pre-date modern anthropogenic threats require further study with denser population sampling. Our study highlights the value that conservation genomics studies of island radiations have as windows onto the possible future for the world's biota as climate change and habitat destruction increasingly fragment their ranges and contribute to rapid declines in population abundances.

Cryptic Speciation of the Oriental Greenfinch Chloris sinica on Oceanic Islands.

The Oriental greenfinch, Chloris sinica, is a small seed-eating finch that breeds in the eastern Palearctic region, an area that spans from Russia in the east to China, Korea, and Japan in the south and southwest. Several subspecies have been described based on subtle morphological characteristics, although the taxonomy varies among different authors. Although many ecological studies have been performed, there has been no phylogenetic study that encompasses the species' entire geographical range. We used four regions of mitochondrial DNA to analyze the intraspecies genetic phylogeny and diversity of the Oriental greenfinch. In addition, we performed morphometric analyses using museum specimens. Genetic analysis identified two clades that diverged approximately 1.06 million years ago. These were a population from the Ogasawara Islands, Japan (subspecies kittlitzi, Clade B), and the other populations (Clade A, which could not be subdivided according to geographic context). Morphometric analyses showed that the population on the Kuril Islands (subspecies kawarahiba) had the longest mean wing length, whereas C. s. kittlitzi had the shortest wings. Chloris s. kittlitzi also had the longest mean bill length, probably because it has adapted to feeding on the Ogasawara Islands. Based on molecular phylogeny and morphology analyses, we recommend that C. s. kittlitzi should be treated as a completely distinct species, called the Ogasawara greenfinch, Chloris kittlitzi. It is critically endangered and needs to be specially protected.

Explosive radiation and spatial expansion across the cold environments of the Old World in an avian family.

Our objective was to elucidate the biogeography and speciation patterns in an entire avian family, which shows a complex pattern of overlapping and nonoverlapping geographical distributions, and much variation in plumage, but less in size and structure. We estimated the phylogeny and divergence times for all of the world's species of Prunella based on multiple genetic loci, and analyzed morphometric divergence and biogeographical history. The common ancestor of Prunella was present in the Sino-Himalayan Mountains or these mountains and Central Asia-Mongolia more than 9 million years ago (mya), but a burst of speciations took place during the mid-Pliocene to early Pleistocene. The relationships among the six primary lineages resulting from that differentiation are unresolved, probably because of the rapid radiation. A general increase in sympatry with increasing time since divergence is evident. With one exception, species in clades younger than c. 3.7 my are allopatric. Species that are widely sympatric, including the most recently diverged (2.4 mya) sympatric sisters, are generally more divergent in size/structure than allo-/parapatric close relatives. The distributional pattern and inferred ages suggest divergence in allopatry and substantial waiting time until secondary contact, likely due to competitive exclusion. All sympatrically breeding species are ecologically segregated, as suggested by differences in size/structure and habitat. Colonizations of new areas were facilitated during glacial periods, followed by fragmentation during interglacials-contrary to the usual view that glacial periods resulted mainly in fragmentations.

DNA barcoding reveals 24 distinct lineages as cryptic bird species candidates in and around the Japanese Archipelago.

DNA barcoding using a partial region (648 bp) of the cytochrome c oxidase I (COI) gene is a powerful tool for species identification and has revealed many cryptic species in various animal taxa. In birds, cryptic species are likely to occur in insular regions like the Japanese Archipelago due to the prevention of gene flow by sea barriers. Using COI sequences of 234 of the 251 Japanese-breeding bird species, we established a DNA barcoding library for species identification and estimated the number of cryptic species candidates. A total of 226 species (96.6%) had unique COI sequences with large genetic divergence among the closest species based on neighbour-joining clusters, genetic distance criterion and diagnostic substitutions. Eleven cryptic species candidates were detected, with distinct intraspecific deep genetic divergences, nine lineages of which were geographically separated by islands and straits within the Japanese Archipelago. To identify Japan-specific cryptic species from trans-Paleartic birds, we investigated the genetic structure of 142 shared species over an extended region covering Japan and Eurasia; 19 of these species formed two or more clades with high bootstrap values. Excluding six duplicated species from the total of 11 species within the Japanese Archipelago and 19 trans-Paleartic species, we identified 24 species that were cryptic species candidates within and surrounding the Japanese Archipelago. Repeated sea level changes during the glacial and interglacial periods may be responsible for the deep genetic divergences of Japanese birds in this insular region, which has led to inconsistencies in traditional taxonomies based on morphology.

Distribution of two distinctive mitochondrial DNA lineages of the Japanese robin Luscinia akahige across its breeding range around the Japanese islands.

The genetic structure of the Japanese robin Luscinia akahige across its breeding range was investigated based on mitochondrial cytochrome b and a control region sequence to reveal the phylogeographic history of this species, which is endemic to the islands around Japan. We obtained the 2,230-bp combined sequences of 77 robins from eight locations and identified 33 haplotypes. Of these, two phylogenetic clades, Ak and Tn, were evident in both maximum-likelihood and Bayesian analyses. Clade Ak was found throughout the range, whereas clade Tn was endemic to the Izu Islands, resulting in the mixture of both clades at this location. These two mtDNA lineages were inferred to have originated in different refugia during the mid-Pleistocene glacial period, probably one in southwestern Japan and the other in the Izu Islands. The population in the Izu Islands has been described as an independent subspecies, tanensis, and is known to be distinguished phenotypically from the nominotypical akahige. Therefore, the distribution pattern of the two mtDNA lineages may be better explained by past introgression of mitochondrial DNA rather than by assuming continuous hybridization with immigrants from the mainland. Future nuclear DNA analysis is needed to evaluate the effect of this introgression event on the genetic bases of the current Izu Islands population. Within clade Ak, the haplotype distribution pattern was inferred to be primarily a result of contiguous range expansion, reflecting the history of both range and population expansion after the last glacial period, with some effects of allopatric fragmentation in its southwestern range.

Old divergences in a boreal bird supports long-term survival through the Ice Ages.

BACKGROUND: Unlike northern Europe and most of northern North America, the Eastern Palearctic and the northwesternmost tip of North America are believed to have been almost unglaciated during the Quarternary glacial periods. This could have facilitated long-term survival of many organisms in that area. To evaluate this, we studied the phylogeography in east Asia and Alaska of a boreal migratory passerine bird, the Arctic Warbler Phylloscopus borealis, and compared our results with published data on especially North American species. RESULTS: In a sample of 113 individuals from 18 populations we identified 42 haplotypes of the mitochondrial cytochrome b gene, which separated into three clades: A--Alaska and mainland Eurasia (except Kamchatka); B--Kamchatka, Sakhalin and Hokkaido; and C--Honshu, Shikoku and Kyushu (i.e. Japan except Hokkaido). The oldest split among these clades, between A/B and C, is estimated to have taken place sometime between the mid Pliocene and early Pleistocene, and the second divergence, between clades A and B, in the early to mid Pleistocene. Within all of the three main clades, there are signs of population expansion. CONCLUSIONS: The Arctic Warbler separated into three main clades in close succession around the Pliocene/Pleistocene border, with the two northern clades diverging last. All three clades probably experienced population bottlenecks during the Pleistocene as a result of range shifts and contractions, but nevertheless survived and maintained their integrities. Several other clades of Northeastern Palearctic birds are noted to have diversified during the Pliocene. In contrast, avian species or phylogroups presently occupying formerly glaciated North American ground are generally younger. The differences between these regions could be due to slower speciation rates in the Eastern Palearctic due to less fragmentation of forest habitats during glacial periods, or to longer survival of Eastern Palearctic clades as a result of less severe conditions in that region compared to northern North America. Several other Palearctic organisms show concordant biogeographical patterns to that of the Arctic Warbler, indicating common causes of their diversifications.