Being cosmopolitan: evolutionary history and phylogeography of a specialized raptor, the Osprey Pandion haliaetus

The evolutionary history of "suboptimal" migration routes

Migratoriness in birds is evolutionary labile, with many examples of increasing or decreasing migration distances on the timescale of modern ornithology. In contrast, shifts of migration to more nearby wintering grounds seem to be a slow process. We examine the history of how Palearctic migratory landbirds have expanded their wintering ranges to include both tropical Africa and Asia, a process that has involved major shifts in migratory routes. We found that species with shorter migration distances and with resident populations in the Palearctic more often winter in both Africa and Asia. Our results suggest that changes in wintering grounds are not by long-distance migrant populations per se, but through historic intermediate populations that were less migratory from which long-distance migration evolved secondarily. The failure of long-distance migrants to shift migration direction to more nearby winter quarters indicates that major modifications to the migratory program may be difficult to evolve.

DNA metabarcoding suggests adaptive seasonal variation of individual trophic traits in a critically endangered fish

Dietary studies are critical for understanding foraging strategies and have important applications in conservation and habitat management. We applied a robust metabarcoding protocol to characterize the diet of the critically endangered freshwater fish Zingel asper (the Rhone streber). We conducted modelling and simulation analyses to identify and characterize some of the drivers of individual trophic trait variation in this species. We found that population density and ontogeny had minor effects on the trophic niche of Z. asper. Instead, our results suggest that the majority of trophic niche variation was driven by seasonal variation in ecological opportunity. The total trophic niche width of Z. asper seasonally expanded to include a broader range of prey. Furthermore, null model simulations revealed that the increase of between-individual variation in autumn indicates that Z. asper become more opportunistic relative to summer and spring, rather than being associated with a seasonal specialization of individuals. Overall, our results suggest an adaptive variation of individual trophic traits in Z. asper: the species mainly consumes a few ephemeropteran taxa (Baetis fuscatus and Ecdyonurus) but seems to be capable of adapting its foraging strategy to maintain its body condition. This study illustrates how metabarcoding data obtained from faeces can be validated and combined with individual-based modelling and simulation approaches to explore inter- and intrapopulational individual trophic traits variation and to test hypotheses in the conventional analytic framework of trophic ecology.

Common patterns in the molecular phylogeography of western palearctic birds: a comprehensive review

A plethora of studies have offered crucial insights in the phylogeographic status of Western Palearctic bird species. However, an overview integrating all this information and analyzing the combined results is still missing. In this study, we compiled all published peer-reviewed and grey literature available on the phylogeography of Western Palearctic bird species. Our literature review indicates a total number of 198 studies, with the overwhelming majority published as journal articles (n = 186). In total, these literature items offer information on 145 bird species. 85 of these species are characterized by low genetic differentiation, 46 species indicate genetic variation but no geographic structuring i.e. panmixia, while 14 species show geographically distinct lineages and haplotypes. Majority of bird species inhabiting the Western Palearctic display genetic admixture. The glaciation cycles in the past few million years were pivotal factors in shaping this situation: during warm periods many species expanded their distribution range to the north over wide areas of Eurasia; whereas, during ice ages most areas were no longer suitable and species retreated to refugia, where lineages mixed.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s10336-021-01893-x.

Thermoregulation and heat exchange in ospreys (Pandion haliaetus)

The osprey (Pandion haliaetus) is a cosmopolitan and long-distant migrant, found at all thermal extremes ranging from polar to tropical climates. Since ospreys may have an unusually flexible thermal physiology due to their migration over, and use of, a wide range of habitats, they represent an interesting study system to explore thermoregulatory adaptations in a raptor. In this study, we investigated the efficiency of heat exchange between body and environment in ospreys using micro-computed tomography (μ-CT), infrared thermography and behavioral observations. μ-CT revealed that the osprey bill has its largest potential for heat exchange at the proximal bill region, where arteries are situated most closely under the surface. However, thermal images of 10 juvenile ospreys showed that the bill contributes to only 0.3% of the bird's total heat exchange. The long legs and protruding claws played a more prominent role as heat dissipation areas with a contribution of 6% and 7%, respectively. Operative thresholds, i.e. the ambient temperature below which heat is lost, were high (>38.5 °C) in these body parts. However, we found no indication of active regulation of heat exchange. Instead we observed multiple behavioral adaptations starting at relatively low ambient temperatures. At 26.3 °C ospreys had a 50% probability of showing panting behavior and above 27.9 °C they additionally spread their wings to enable heat dissipation from the less insulated ventral side. The thermal images revealed that at an ambient temperature of 32.1 °C ospreys had a 50% probability of developing a ≥2 °C and up to 7.5 °C colder stripe on the head, which was likely caused by cutaneous evaporation. Our observations suggest that ospreys more strongly rely on behavioral mechanisms than on active thermal windows to cope with heat stress. This study not only improves our understanding of the role of different body parts in ospreys' total heat exchange with the environment but further provides an insight about additional adaptations of this raptor to cope with heat stress.

Population history of the golden eagle inferred from whole-genome sequencing of three of its subspecies

The application of evolutionary genetic research to investigate the potential for endangered species to adapt to changing environments is important for conservation biology. Effective population size (Ne) is informative for understanding adaptive potential as it refers to the genetic variation in breeding individuals who have contributed to contemporary and historic population diversity. We reconstruct fluctuations in Ne in three golden eagle subspecies (Japanese, Scottish, North American) using the pairwise sequential Markovian coalescent (PSMC) model based on whole-genome sequence data. Our results indicate the timing of subspeciation events and suggest significant ongoing demographic reductions since the start of the Last Glacial Period. Importantly, we find evidence for gene flow from continental populations into the ancestral Japanese population resulting in a short, sharp recovery in genetic diversity. Timing agrees with the palaeogeographic estimates of land bridge connections between the Japanese archipelago and Asian continent and matches a similar Ne spike in the Scottish population, but not in the North American population. Given contemporary declines in isolated Japanese and UK island populations, our study highlights a concerning loss of local genetic diversity, but also indicates the likely response of populations to genetic reinforcement from neighbouring subspecies, increasing management options and encouraging a range-wide species conservation approach.

Genetic Structure and Diversity of Two Populations of the Eastern Buzzard (Buteo japonicus japonicus and B. j. toyoshimai) in Japan

The nominotypical subspecies of the Eastern buzzard (Buteo japonicus japonicus; BJJ) is a common raptor inhabiting East Asia and Japan. Another subspecies, B. j. toyoshimai (BJT), inhabits only the Bonin Islands of the Ogasawara Islands, where there are only an estimated 85 breeding pairs. Because of this low population size, this subspecies is classified as endangered (class IB) in Japan. The aims of the present study were to examine genetic differences between BJJ and BJT, determine the genetic structure of the Eastern Buzzard, and assess genetic diversity within each subspecies. We sequenced 1526 bp within the control region of the mtDNA of 10 BJJ individuals during the breeding season in four sites; similarly, we sequenced 23 BJJ individuals during winter in three sites. We detected 24 haplotypes among the 33 individuals. In a similar analysis performed with 12 BJT individuals, three haplotypes were detected. The phylogenetic analysis showed that BJJ and BJT have diverged into distinct clades, supporting the genetic differentiation between the subspecies. Network and mismatch distribution analyses indicated that BJJ may have experienced population expansion. In addition, comparisons with other raptors revealed a high degree of genetic diversity in the BJJ population. In contrast, the genetic diversity of the BJT population is lower than that in other raptors. Our results indicated that it is necessary to protect BJT to prevent the reduction in its genetic diversity.

New insights into population structure of the European golden eagle (Aquila chrysaetos) revealed by microsatellite analysis

Connectivity between golden eagle (Aquila chrysaetos) populations is poorly understood. Field studies exploring natal dispersal suggest that this raptor is a philopatric species, but with the ability to roam far. However, little is known about the population structure of the species in Europe. Our study is based on 14 microsatellite loci and is complemented by new and previously published mitochondrial control region DNA data. The present dataset includes 121 eagles from Scotland, Norway, Finland, Estonia, the Mediterranean and Alpine regions. Our sampling focused on the Alpine and Mediterranean populations because both mitochondrial DNA (mtDNA) lineages found in golden eagles, the Holarctic and the Mediterranean, are known to co-occur there. Cluster analyses of nuclear DNA support a shallow split into northern and southern populations in Europe, similar to the distribution of the two mtDNA lineages, with the Holarctic lineage occurring in the north and the Mediterranean lineage predominating in the south. Additionally, Scotland shows significant differentiation and low relative migration levels that indicate isolation from the mainland populations. Alpine and Mediterranean golden eagles do not show nuclear structure corresponding to divergent mtDNA lineages. This indicates that the presence of northern Holarctic mitochondrial haplotypes in the Alps and the Mediterranean is attributable to past admixture rather than recent long-distance dispersal.

Migrating ospreys use thermal uplift over the open sea

Most large raptors on migration avoid crossing the sea because of the lack of atmospheric convection over temperate seas. The osprey Pandion haliaetus is an exception among raptors, since it can fly over several hundred kilometres of open water. We equipped five juvenile ospreys with GPS-Accelerometer-Magnetometer loggers. All birds were able to find and use thermal uplift while crossing the Mediterranean Sea, on average 7.5 times per 100 km, and could reach altitudes of 900 m above the sea surface. Their climb rate was 1.6 times slower than over land, and birds kept flapping most of the time while circling in the thermals, indicating that convections cells were weaker than over land. The frequency of thermal soaring was correlated with the difference between the sea surface and air temperature, indicating that atmospheric convection occurred when surface waters were warmer than the overlaying air. These observations help explain the transoceanic cosmopolitan distribution of osprey, and question the widely held assumption that water bodies represent strict barriers for large raptors.

One-locus-several-primers: A strategy to improve the taxonomic and haplotypic coverage in diet metabarcoding studies

In diet metabarcoding analyses, insufficient taxonomic coverage of PCR primer sets generates false negatives that may dramatically distort biodiversity estimates. In this paper, we investigated the taxonomic coverage and complementarity of three cytochrome c oxidase subunit I gene (COI) primer sets based on in silico analyses and we conducted an in vivo evaluation using fecal and spider web samples from different invertivores, environments, and geographic locations. Our results underline the lack of predictability of both the coverage and complementarity of individual primer sets: (a) sharp discrepancies exist observed between in silico and in vivo analyses (to the detriment of in silico analyses); (b) both coverage and complementarity depend greatly on the predator and on the taxonomic level at which preys are considered; (c) primer sets' complementarity is the greatest at fine taxonomic levels (molecular operational taxonomic units [MOTUs] and variants). We then formalized the "one-locus-several-primer-sets" (OLSP) strategy, that is, the use of several primer sets that target the same locus (here the first part of the COI gene) and the same group of taxa (here invertebrates). The proximal aim of the OLSP strategy is to minimize false negatives by increasing total coverage through multiple primer sets. We illustrate that the OLSP strategy is especially relevant from this perspective since distinct variants within the same MOTUs were not equally detected across all primer sets. Furthermore, the OLSP strategy produces largely overlapping and comparable sequences, which cannot be achieved when targeting different loci. This facilitates the use of haplotypic diversity information contained within metabarcoding datasets, for example, for phylogeography and finer analyses of prey-predator interactions.

Unravelling population processes over the Late Pleistocene driving contemporary genetic divergence in Palearctic buzzards

Population range expansions and contractions as a response to climate and habitat change throughout the Quaternary are known to have contributed to complex phylogenetic and population genetic events. Speciation patterns and processes in Palearctic buzzards (genus Buteo) are a long-standing example of morphological and genetic data incongruence, attributed to panmixia, habitat range shifts, contact zones, and climate change. Here we assess the systematics, phylogeography and population genetic structure of three nominal species of Palearctic buzzards, Buteo buteo (including B. b. vulpinus), B. rufinus (including B. r. cirtensis) and B. hemilasius. Phylogenetic analyses inferred from mitochondrial data recover B. hemilasius as sister to the sister clades B. r. rufinus and B. buteo complex (B. b. buteo, B. b. vulpinus, but also including B. r. cirtensis). In contrast, we find an unresolved genetic delimitation inferred from four nuclear loci, suggesting an ancestral genetic pool for all species. Time-trees suggest population contractions and expansions throughout the Pleistocene, which likely reflect habitat change and contrasting ecological niche requirements between species. Microsatellite-based extended Bayesian skyline plots reveal relatively constant population sizes for B. hemilasius, B. r. rufinus, and B. b. vulpinus, in contrast to a dramatic population expansion in B. r. cirtensis within the last 3 kya. Overall, our study illustrates how complex population processes over the Late Pleistocene have shaped the patterns of genetic divergence in Palearctic buzzards, due to the joint effects of shared ancestral polymorphisms, population expansions and contractions, with hybridization at contact zones leading to admixture and introgression.

Analysis of the Genetic Diversity and Structure of the Grey-Faced Buzzard (Butastur indicus) in Japan, Based on mtDNA

The grey-faced buzzard (Butastur indicus) is a raptor that inhabits East Asia, including Japan. Because the number of individuals has decreased by 75% over the last 40 years, this species is classified as vulnerable (VU) in Japan. In the present study, wesought to reveal the genetic structure of the Japanese grey-faced buzzard population at several breeding sites, and to assess the levels of genetic diversity within the Japanese population. We sequenced 555 bp of the mitochondrial DNA of 96 individuals sampled during the breeding season at 18 sites, and 11 individuals sampled during the winter season at one site. In total, 21 variable sites were found in the control region, and we detected 26 haplotypes among the 107 individuals. Fukuoka represented the core breeding area for grey-faced buzzards, as half of all haplotypes were detected there. Four unique haplotypes were detected in the overwintering area. The results of the network and mismatch distribution analyses indicated that the grey-faced buzzard has not experienced a genetic bottleneck in the past, but did experience recent population expansion. In addition, comparisons with other raptors revealed rich genetic diversity in the grey-faced buzzard population. Our results indicate that conservation of both breeding and wintering areas is important for the protection of the grey-faced buzzard.

Comprehensive molecular phylogeny of barn owls and relatives (Family: Tytonidae), and their six major Pleistocene radiations

The owl family Tytonidae comprises two genera: Phodilus, limited to the forests of central Africa and South-East Asia, and the ubiquitous Tyto. The genus Tyto is majorly represented by the cosmopolitan Common Barn Owl group, with more than 30 subspecies worldwide. Discrete differences in body size and plumage colouration have led to the classification of this family into many species and subspecies, but the taxonomic status and phylogenetic relationships between taxa remain unclear, and in some groups controversial. Although several previous studies attempted to resolve this problem, they have been limited in their taxonomic and geographical coverage, or have relied on restricted molecular evidence and low sample sizes. Based on the most comprehensive sampling to date (16 out of 17 Tyto species, and one out of three Phodilus species), a multi-locus approach using seven mitochondrial and two nuclear markers, and taking advantage of field data and museum collections available worldwide, our main questions in this study were: (1) what are the phylogenetic relationships and classification status of the whole family; (2) when and where did the most important speciation events occur? We confirm that the Common Barn Owl, Tyto alba is divided into three main evolutionary units: the American Barn Owl, T. furcata; the Western Barn Owl, T. alba; and the Eastern Barn Owl, T. javanica, and suggest a Late Miocene (ca. 6 mya) Australasian and African origin of the group. Our results are supported by fossil age information, given that the most recent common ancestor between the Tytonidae genera Phodilus and Tyto was probably from the Oligocene (ca. 28 mya) of Australasia. We finally reveal six major Pleistocene radiations of Tyto, all resulting in wide-range distributions.

A population genomics approach shows widespread geographical distribution of cryptic genomic forms of the symbiotic fungus Rhizophagus irregularis

Arbuscular mycorrhizal fungi (AMF; phylum Gomeromycota) associate with plants forming one of the most successful microbe-plant associations. The fungi promote plant diversity and have a potentially important role in global agriculture. Plant growth depends on both inter- and intra-specific variation in AMF. It was recently reported that an unusually large number of AMF taxa have an intercontinental distribution, suggesting long-distance gene flow for many AMF species, facilitated by either long-distance natural dispersal mechanisms or human-assisted dispersal. However, the intercontinental distribution of AMF species has been questioned because the use of very low-resolution markers may be unsuitable to detect genetic differences among geographically separated AMF, as seen with some other fungi. This has been untestable because of the lack of population genomic data, with high resolution, for any AMF taxa. Here we use phylogenetics and population genomics to test for intra-specific variation in Rhizophagus irregularis, an AMF species for which genome sequence information already exists. We used ddRAD sequencing to obtain thousands of markers distributed across the genomes of 81 R. irregularis isolates and related species. Based on 6 888 variable positions, we observed significant genetic divergence into four main genetic groups within R. irregularis, highlighting that previous studies have not captured underlying genetic variation. Despite considerable genetic divergence, surprisingly, the variation could not be explained by geographical origin, thus also supporting the hypothesis for at least one AMF species of widely dispersed AMF genotypes at an intercontinental scale. Such information is crucial for understanding AMF ecology, and how these fungi can be used in an environmentally safe way in distant locations.

A from-benchtop-to-desktop workflow for validating HTS data and for taxonomic identification in diet metabarcoding studies

The main objective of this work was to develop and validate a robust and reliable "from-benchtop-to-desktop" metabarcoding workflow to investigate the diet of invertebrate-eaters. We applied our workflow to faecal DNA samples of an invertebrate-eating fish species. A fragment of the cytochrome c oxidase I (COI) gene was amplified by combining two minibarcoding primer sets to maximize the taxonomic coverage. Amplicons were sequenced by an Illumina MiSeq platform. We developed a filtering approach based on a series of nonarbitrary thresholds established from control samples and from molecular replicates to address the elimination of cross-contamination, PCR/sequencing errors and mistagging artefacts. This resulted in a conservative and informative metabarcoding data set. We developed a taxonomic assignment procedure that combines different approaches and that allowed the identification of ~75% of invertebrate COI variants to the species level. Moreover, based on the diversity of the variants, we introduced a semiquantitative statistic in our diet study, the minimum number of individuals, which is based on the number of distinct variants in each sample. The metabarcoding approach described in this article may guide future diet studies that aim to produce robust data sets associated with a fine and accurate identification of prey items.