A parapatric propensity for breeding precludes the completion of speciation in common teal (Anas crecca, sensu lato)

Molecular and acoustic evidence support the species status of Anthus rubescens rubescens and Anthus [rubescens] japonicus (Passeriformes: Motacillidae)

The Buff-bellied Pipit Anthus rubescens comprises two allopatric subspecies groups: A. r. rubescens and A. r. alticola in North America and A. [r.] japonicus in north-east Asia. Despite their great morphological resemblance in breeding plumage, most individuals can be assigned to one or the other subspecies group in non-breeding plumage. Allopatric distributions, morphological differentiation and previously reported molecular divergence suggested the need for additional taxonomic study to assess the rank of these two populations. To resolve the taxonomy of the Buff-bellied Pipit species complex we analysed i) two mitochondrial DNA (mtDNA) loci and ii) nine bioacoustic parameters across 69 sound recordings (338 flight calls) recovered from public databases using principal component analysis and Euclidean distance measures. By comparing our mtDNA and call divergence measures with similar values measured between long-recognised species pairs of the genus, we show that the level of mitochondrial and acoustic divergence between the two Buff-bellied Pipit subspecies groups is typical of species-level divergence in the genus Anthus. Therefore, we recommend splitting the Buff-bellied Pipit species complex into two species: Anthus rubescens (American Pipit) and Anthus japonicus (Siberian Pipit). Our results also suggest that the Water Pipit A. spinoletta deserves taxonomic reassessment as its lineages are highly divergent in acoustics and mtDNA, while mtDNA relationships suggest paraphyly relative to the Rock Pipit A. petrosus. Our work highlights the crucial importance of integrative approaches in taxonomy and the usefulness of bioacoustics in studying cryptic diversity.

Beringia as a high-latitude engine of avian speciation

Beringia is a biogeographically dynamic region that extends from northeastern Asia into northwestern North America. This region has affected avian divergence and speciation in three important ways: (i) by serving as a route for intercontinental colonisation between Asia and the Americas; (ii) by cyclically splitting (and often reuniting) populations, subspecies, and species between these continents; and (iii) by providing isolated refugia through glacial cycles. The effects of these processes can be seen in taxonomic splits of shallow to increasing depths and in the presence of regional endemics. We review the taxa involved in the latter two processes (splitting-reuniting and isolation), with a focus on three research topics: avian diversity, time estimates of the generation of that diversity, and the regions within Beringia that might have been especially important. We find that these processes have generated substantial amounts of avian diversity, including 49 pairs of avian subspecies or species whose breeding distributions largely replace one another across the divide between the Old World and the New World in Beringia, and 103 avian species and subspecies endemic to this region. Among endemics, about one in three is recognised as a full biological species. Endemic taxa in the orders Charadriiformes (shorebirds, alcids, gulls, and terns) and Passeriformes (perching birds) are particularly well represented, although they show very different levels of diversity through evolutionary time. Endemic Beringian Charadriiformes have a 1.31:1 ratio of species to subspecies. In Passeriformes, endemic taxa have a 0.09:1 species-to-subspecies ratio, suggesting that passerine (and thus terrestrial) endemism might be more prone to long-term extinction in this region, although such 'losses' could occur through their being reconnected with wider continental populations during favourable climatic cycles (e.g. subspecies reintegration with other populations). Genetic evidence suggests that most Beringian avian taxa originated over the past 3 million years, confirming the importance of Quaternary processes. There seems to be no obvious clustering in their formation through time, although there might be temporal gaps with lower rates of diversity generation. For at least 62 species, taxonomically undifferentiated populations occupy this region, providing ample potential for future evolutionary diversification.

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

The processes leading to divergence and speciation can differ broadly among taxa with different life histories. We examine these processes in a small clade of ducks with historically uncertain relationships and species limits. The green-winged teal (Anas crecca) complex is a Holarctic species of dabbling duck currently categorized as three subspecies (Anas crecca crecca, A. c. nimia, and A. c. carolinensis) with a close relative, the yellow-billed teal (Anas flavirostris) from South America. A. c. crecca and A. c. carolinensis are seasonal migrants, while the other taxa are sedentary. We examined divergence and speciation patterns in this group, determining their phylogenetic relationships and the presence and levels of gene flow among lineages using both mitochondrial and genome-wide nuclear DNA obtained from 1,393 ultraconserved element (UCE) loci. Phylogenetic relationships using nuclear DNA among these taxa showed A. c. crecca, A. c. nimia, and A. c. carolinensis clustering together to form one polytomous clade, with A. flavirostris sister to this clade. This relationship can be summarized as (crecca, nimia, carolinensis)(flavirostris). However, whole mitogenomes revealed a different phylogeny: (crecca, nimia)(carolinensis, flavirostris). The best demographic model for key pairwise comparisons supported divergence with gene flow as the probable speciation mechanism in all three contrasts (crecca-nimia, crecca-carolinensis, and carolinensis-flavirostris). Given prior work, gene flow was expected among the Holarctic taxa, but gene flow between North American carolinensis and South American flavirostris (M ∼0.1-0.4 individuals/generation), albeit low, was not expected. Three geographically oriented modes of divergence are likely involved in the diversification of this complex: heteropatric (crecca-nimia), parapatric (crecca-carolinensis), and (mostly) allopatric (carolinensis-flavirostris). Our study shows that ultraconserved elements are a powerful tool for simultaneously studying systematics and population genomics in systems with historically uncertain relationships and species limits.

Estimating Movement Rates Between Eurasian and North American Birds That Are Vectors of Avian Influenza

Avian influenza (AI) is a zoonotic disease that will likely be involved in future pandemics. Because waterbird movements are difficult to quantify, determining the host-specific risk of Eurasian-origin AI movements into North America is challenging. We estimated relative rates of movements, based on long-term evolutionary averages of gene flow, between Eurasian and North American waterbird populations to obtain bidirectional baseline rates of the intercontinental movements of these AI hosts. We used population genomics and coalescent-based demographic models to obtain these gene-flow-based movement estimates. Inferred rates of movement between these continental populations varies greatly among species. Within dabbling ducks, gene flow, relative to effective population size, varies from ∼3 to 24 individuals/generation between Eurasian and American wigeons (Mareca penelope and Mareca americana) to ∼100-300 individuals/generation between continental populations of northern pintails (Anas acuta). These are evolutionary long-term averages and provide a solid foundation for understanding the relative risks of each of these host species in potential intercontinental AI movements. We scale these values to census size for evaluation in that context. In addition to being AI hosts, many of these bird species are also important in the subsistence diets of Alaskans, increasing the risk of direct bird-to-human exposure to Eurasian-origin AI virus. We contrast species-specific rates of intercontinental movements with the importance of each species in Alaskan diets to understand the relative risk of these taxa to humans. Assuming roughly equivalent AI infection rates among ducks, greater scaup (Aythya marila), mallard (Anas platyrhynchos), and northern pintail (Anas acuta) were the top three species presenting the highest risks for intercontinental AI movement both within the natural system and through exposure to subsistence hunters. Improved data on AI infection rates in this region could further refine these relative risk assessments. These directly comparable, species-based intercontinental movement rates and relative risk rankings should help in modeling, monitoring, and mitigating the impacts of intercontinental host and AI movements.

Avian introgression patterns are consistent with Haldane's Rule

According to Haldane’s Rule, the heterogametic sex will show the greatest fitness reduction in a hybrid cross. In birds, where sex is determined by a ZW system, female hybrids are expected to experience lower fitness compared to male hybrids. This pattern has indeed been observed in several bird groups, but it is unknown whether the generality of Haldane’s Rule also extends to the molecular level. First, given the lower fitness of female hybrids, we can expect maternally inherited loci (i.e., mitochondrial and W-linked loci) to show lower introgression rates than biparentally inherited loci (i.e., autosomal loci) in females. Second, the faster evolution of Z-linked loci compared to autosomal loci and the hemizygosity of the Z-chromosome in females might speed up the accumulation of incompatible alleles on this sex chromosome, resulting in lower introgression rates for Z-linked loci than for autosomal loci. I tested these expectations by conducting a literature review which focused on studies that directly quantified introgression rates for autosomal, sex-linked, and mitochondrial loci. Although most studies reported introgression rates in line with Haldane’s Rule, it remains important to validate these genetic patterns with estimates of hybrid fitness and supporting field observations to rule out alternative explanations. Genomic data provide exciting opportunities to obtain a more fine-grained picture of introgression rates across the genome, which can consequently be linked to ecological and behavioral observations, potentially leading to novel insights into the genetic mechanisms underpinning Haldane’s Rule.

Divergence, gene flow, and speciation in eight lineages of trans-Beringian birds

Determining how genetic diversity is structured between populations that span the divergence continuum from populations to biological species is key to understanding the generation and maintenance of biodiversity. We investigated genetic divergence and gene flow in eight lineages of birds with a trans-Beringian distribution, where Asian and North American populations have likely been split and reunited through multiple Pleistocene glacial cycles. Our study transects the speciation process, including eight pairwise comparisons in three orders (ducks, shorebirds and passerines) at population, subspecies and species levels. Using ultraconserved elements (UCEs), we found that these lineages represent conditions from slightly differentiated populations to full biological species. Although allopatric speciation is considered the predominant mode of divergence in birds, all of our best divergence models included gene flow, supporting speciation with gene flow as the predominant mode in Beringia. In our eight lineages, three were best described by a split-migration model (divergence with gene flow), three best fit a secondary contact scenario (isolation followed by gene flow), and two showed support for both models. The lineages were not evenly distributed across a divergence space defined by gene flow (M) and differentiation (F_ST ), instead forming two discontinuous groups: one with relatively shallow divergence, no fixed single nucleotide polymorphisms (SNPs), and high rates of gene flow between populations; and the second with relatively deeply divergent lineages, multiple fixed SNPs, and low gene flow. Our results highlight the important role that gene flow plays in avian divergence in Beringia.

Genetic, phenotypic and ecological differentiation suggests incipient speciation in two Charadrius plovers along the Chinese coast

Background

Speciation with gene flow is an alternative to the nascence of new taxa in strict allopatric separation. Indeed, many taxa have parapatric distributions at present. It is often unclear if these are secondary contacts, e.g. caused by past glaciation cycles or the manifestation of speciation with gene flow, which hampers our understanding of how different forces drive diversification. Here we studied genetic, phenotypic and ecological aspects of divergence in a pair of incipient shorebird species, the Kentish (Charadrius alexandrinus) and the White-faced Plovers (C. dealbatus), shorebirds with parapatric breeding ranges along the Chinese coast. We assessed divergence based on molecular markers with different modes of inheritance and quantified phenotypic and ecological divergence in aspects of morphometric, dietary and climatic niches.

Results

Our integrative analyses revealed small to moderate levels of genetic and phenotypic distinctiveness with symmetric gene flow across the contact area at the Chinese coast. The two species diverged approximately half a million years ago in dynamic isolation with secondary contact occurring due to cycling sea level changes between the Eastern and Southern China Sea in the mid-late Pleistocene. We found evidence of character displacement and ecological niche differentiation between the two species, invoking the role of selection in facilitating divergence despite gene flow.

Conclusion

These findings imply that ecology can indeed counter gene flow through divergent selection and thus contributes to incipient speciation in these plovers. Furthermore, our study highlights the importance of using integrative datasets to reveal the evolutionary history and assist the inference of mechanisms of speciation.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1449-5) contains supplementary material, which is available to authorized users.

Avian influenza

Previous introductions of highly pathogenic avian influenza virus (HPAIV) to the EU were most likely via migratory wild birds. A mathematical model has been developed which indicated that virus amplification and spread may take place when wild bird populations of sufficient size within EU become infected. Low pathogenic avian influenza virus (LPAIV) may reach similar maximum prevalence levels in wild bird populations to HPAIV but the risk of LPAIV infection of a poultry holding was estimated to be lower than that of HPAIV. Only few non-wild bird pathways were identified having a non-negligible risk of AI introduction. The transmission rate between animals within a flock is assessed to be higher for HPAIV than LPAIV. In very few cases, it could be proven that HPAI outbreaks were caused by intrinsic mutation of LPAIV to HPAIV but current knowledge does not allow a prediction as to if, and when this could occur. In gallinaceous poultry, passive surveillance through notification of suspicious clinical signs/mortality was identified as the most effective method for early detection of HPAI outbreaks. For effective surveillance in anseriform poultry, passive surveillance through notification of suspicious clinical signs/mortality needs to be accompanied by serological surveillance and/or a virological surveillance programme of birds found dead (bucket sampling). Serosurveillance is unfit for early warning of LPAI outbreaks at the individual holding level but could be effective in tracing clusters of LPAIV-infected holdings. In wild birds, passive surveillance is an appropriate method for HPAIV surveillance if the HPAIV infections are associated with mortality whereas active wild bird surveillance has a very low efficiency for detecting HPAIV. Experts estimated and emphasised the effect of implementing specific biosecurity measures on reducing the probability of AIV entering into a poultry holding. Human diligence is pivotal to select, implement and maintain specific, effective biosecurity measures.

Complex history of isolation and gene flow in hoary, Olympic, and endangered Vancouver Island marmots

Climate change resulting in a reduction of alpine habitat is believed to pose a considerable risk to alpine-dependent species, including many marmots. Hoary marmots (Marmota caligata) range throughout much of the mountainous Pacific Northwest (PNW) and Rocky Mountains while the closely related Olympic and Vancouver Island marmots (M. olympus and M. vancouverensis, respectively) are restricted to small isolated regions of the PNW. The endemic Vancouver Island marmot is currently classified as Critically Endangered and the Olympic marmot has recently experienced dramatic population declines. Previous phylogenetic studies of PNW marmot species have had limited power as they focused on resolving interspecific relationships, implicitly assumed an absence of gene flow among currently recognized species, included relatively few individuals, and relied heavily or entirely on mitochondrial DNA. We sequenced 2 mitochondrial and 4 nuclear markers from 167 hoary, 4 Vancouver Island, and 5 Olympic marmots in order to investigate phylogenetic relationships and historic gene flow among these species. We recovered 2 monophyletic (and predominantly allopatric) mitochondrial clades of hoary marmots that are not sister groups. Instead, Vancouver Island marmots formed a monophyletic mitochondrial sister clade to 1 of the hoary marmot clades. Nuclear loci did not recover the 2 mitochondrial clades of hoary marmots and suggest that Vancouver Island marmots may have experienced mitochondrial introgression from coastal mainland hoary marmots. Additionally, our nuclear results suggest possible gene flow between hoary and Olympic marmots despite different chromosomal formulas. Rather than resolving what has previously been considered a straightforward 3-taxon phylogenetic question, our findings suggest a complicated history of rapid divergence of the 3 species followed by intermittent and possibly ongoing gene flow between hoary marmots and both Olympic and Vancouver Island marmots. These results therefore have significant implications for the conservation of the latter 2 species, both of which are conservation concerns.

Understanding the mechanisms of antitropical divergence in the seabird White-faced Storm-petrel (Procellariiformes: Pelagodroma marina) using a multilocus approach

Analytical methods that apply coalescent theory to multilocus data have improved inferences of demographic parameters that are critical to understanding population divergence and speciation. In particular, at the early stages of speciation, it is important to implement models that accommodate conflicting gene trees, and benefit from the presence of shared polymorphisms. Here, we employ eleven nuclear loci and the mitochondrial control region to investigate the phylogeography and historical demography of the pelagic seabird White-faced Storm-petrel (Pelagodroma marina) by sampling subspecies across its antitropical distribution. Groups are all highly differentiated: global mitochondrial ΦST = 0.89 (P < 0.01) and global nuclear ΦST varies between 0.22 and 0.83 (all P < 0.01). The complete lineage sorting of the mitochondrial locus between hemispheres is corroborated by approximately half of the nuclear genealogies, suggesting a long-term antitropical divergence in isolation. Coalescent-based estimates of demographic parameters suggest that hemispheric divergence of P. marina occurred approximately 840 000 ya (95% HPD 582 000-1 170 000), in the absence of gene flow, and divergence within the Southern Hemisphere occurred 190 000 ya (95% HPD 96 000-600 000), both probably associated with the profound palaeo-oceanographic changes of the Pleistocene. A fledgling sampled in St Helena (tropical South Atlantic) suggests recent colonization from the Northern Hemisphere. Despite the great potential for long-distance dispersal, P. marina antitropical groups have been evolving as independent, allopatric lineages, and divergence is probably maintained by philopatry coupled with asynchronous reproductive phenology and local adaptation.

Mito-nuclear discord in six congeneric lineages of Holarctic ducks (genus Anas)

Many species have Holarctic distributions that extend across Europe, Asia and North America. Most genetics research on these species has examined only mitochondrial (mt) DNA, which has revealed wide variance in divergence between Old World (OW) and New World (NW) populations, ranging from shallow, unstructured genealogies to deeply divergent lineages. In this study, we sequenced 20 nuclear introns to test for concordant patterns of OW-NW differentiation between mtDNA and nuclear (nu) DNA for six lineages of Holarctic ducks (genus Anas). Genetic differentiation for both marker types varied widely among these lineages (idiosyncratic population histories), but mtDNA and nuDNA divergence within lineages was not significantly correlated. Moreover, compared with the association between mtDNA and nuDNA divergence observed among different species, OW-NW nuDNA differentiation was generally lower than mtDNA divergence, at least for lineages with deeply divergent mtDNA. Furthermore, coalescent estimates indicated significantly higher rates of gene flow for nuDNA than mtDNA for four of the six lineages. Thus, Holarctic ducks show prominent mito-nuclear discord between OW and NW populations, and we reject differences in sorting rates as the sole cause of the within-species discord. Male-mediated intercontinental gene flow is likely a leading contributor to this discord, although selection could also cause increased mtDNA divergence relative to weak nuDNA differentiation. The population genetics of these ducks contribute to growing evidence that mtDNA can be an unreliable indicator of stage of speciation and that more holistic approaches are needed for species delimitation.

Limited phylogeographic signal in sex-linked and autosomal loci despite geographically, ecologically, and phenotypically concordant structure of mtDNA variation in the Holarctic avian genus Eremophila

Phylogeographic studies of Holarctic birds are challenging because they involve vast geographic scale, complex glacial history, extensive phenotypic variation, and heterogeneous taxonomic treatment across countries, all of which require large sample sizes. Knowledge about the quality of phylogeographic information provided by different loci is crucial for study design. We use sequences of one mtDNA gene, one sex-linked intron, and one autosomal intron to elucidate large scale phylogeographic patterns in the Holarctic lark genus Eremophila. The mtDNA ND2 gene identified six geographically, ecologically, and phenotypically concordant clades in the Palearctic that diverged in the Early - Middle Pleistocene and suggested paraphyly of the horned lark (E. alpestris) with respect to the Temminck's lark (E. bilopha). In the Nearctic, ND2 identified five subclades which diverged in the Late Pleistocene. They overlapped geographically and were not concordant phenotypically or ecologically. Nuclear alleles provided little information on geographic structuring of genetic variation in horned larks beyond supporting the monophyly of Eremophila and paraphyly of the horned lark. Multilocus species trees based on two nuclear or all three loci provided poor support for haplogroups identified by mtDNA. The node ages calculated using mtDNA were consistent with the available paleontological data, whereas individual nuclear loci and multilocus species trees appeared to underestimate node ages. We argue that mtDNA is capable of discovering independent evolutionary units within avian taxa and can provide a reasonable phylogeographic hypothesis when geographic scale, geologic history, and phenotypic variation in the study system are too complex for proposing reasonable a priori hypotheses required for multilocus methods. Finally, we suggest splitting the currently recognized horned lark into five Palearctic and one Nearctic species.

Robust estimates of divergence times and selection with a poisson random field model: a case study of comparative phylogeographic data

Mutation frequencies can be modeled as a Poisson random field (PRF) to estimate speciation times and the degree of selection on newly arisen mutations. This approach provides a quantitative theory for comparing intraspecific polymorphism with interspecific divergence in the presence of selection and can be used to estimate population genetic parameters. Although the original PRF model has been extended to more general biological settings to make statistical inference about selection and divergence among model organisms, it has not been incorporated into phylogeographic studies that focus on estimating population genetic parameters for nonmodel organisms. Here, we modified a recently developed time-dependent PRF model to independently estimate genetic parameters from a nuclear and mitochondrial DNA data set of 22 sister pairs of birds that have diverged across a biogeographic barrier. We found that species that inhabit humid habitats had more recent divergence times and larger effective population sizes than those that inhabit drier habitats, and divergence time estimated from the PRF model were similar to estimates from a coalescent species-tree approach. Selection coefficients were higher in sister pairs that inhabited drier habitats than in those in humid habitats, but overall the mitochondrial DNA was under weak selection. Our study indicates that PRF models are useful for estimating various population genetic parameters and serve as a framework for incorporating estimates of selection into comparative phylogeographic studies.

Global lack of flyway structure in a cosmopolitan bird revealed by a genome wide survey of single nucleotide polymorphisms

Knowledge about population structure and connectivity of waterfowl species, especially mallards (Anas platyrhynchos), is a priority because of recent outbreaks of avian influenza. Ringing studies that trace large-scale movement patterns have to date been unable to detect clearly delineated mallard populations. We employed 363 single nucleotide polymorphism markers in combination with population genetics and phylogeographical approaches to conduct a population genomic test of panmixia in 801 mallards from 45 locations worldwide. Basic population genetic and phylogenetic methods suggest no or very little population structure on continental scales. Nor could individual-based structuring algorithms discern geographical structuring. Model-based coalescent analyses for testing models of population structure pointed to strong genetic connectivity among the world's mallard population. These diverse approaches all support the conclusion that there is a lack of clear population structure, suggesting that the world's mallards, perhaps with minor exceptions, form a single large, mainly interbreeding population.