Mitochondrial DNA analysis reveals Holarctic homogeneity and a distinct Mediterranean lineage in the Golden eagle (Aquila chrysaetos)

A renewed glance at the Palearctic golden eagle: Genetic variation in space and time

Anthropogenic pressures on nature have been causing population declines for centuries. Intensified persecution of apex predators, like the golden eagle, resulted in population bottlenecks during the 19th and 20th centuries. To study population genetics and demographic history of the golden eagle throughout its distribution, we collected museum samples from previously underrepresented regions, such as Russia and Central Asia. We used 12 microsatellite loci and a fragment of the mitochondrial DNA control region to re-evaluate phylogeography of Eurasian golden eagles and study the impacts of the population bottleneck. Our results revealed a north-south genetic gradient, expressed by the difference between Mediterranean and Holarctic lineages, as well as genetically distinct Northern Europe and Central Asia and Caucasus regions. Furthermore, Northern Europe exhibited the lowest, whereas Central Asia and Caucasus had the highest genetic diversity. Although golden eagles maintained relatively high genetic diversity, we detected genetic signatures of the recent bottleneck, including reduced genetic diversity and a decline in the effective female population size around the year 1975. Our study improves the knowledge of the genetic composition of Eurasian golden eagles and highlights the importance of understanding their historical population dynamics in the face of ongoing and future conservation efforts.

Lead, cadmium, and other trace elements in the liver of golden eagles and white-tailed eagles: recent data from Poland and a systematic review of previous studies

The golden eagle (Aquila chrysaetos) and the white-tailed eagle (Haliaeetus albicilla), being apex predators and facultative scavengers, can bioaccumulate different environmental contaminants, including toxic elements that may adversely affect their health. We analyzed the levels of cadmium (Cd), lead (Pb), and other metals and metalloids, including arsenic (As), barium (Ba), beryllium (Be), cobalt (Co), chromium (Cr), copper (Cu), iron (Fe), magnesium (Mg), manganese (Mn), molybdenum (Mo), selenium (Se), thorium (Th), thallium (Tl), uranium (U), vanadium (V), and zinc (Zn) in liver samples taken from three golden eagles and 36 white-tailed eagles that were found dead across Poland to verify their exposure. We also used a systematic review to summarize the available literature data on Cd, Pb, and other studied elements in the liver of both eagle species. Analyses of trace elements in the liver samples of the Polish eagles revealed interspecific differences in Cd, Cu, and Mn and differences in Co, Mn, Tl, and Zn among study regions. All elements tested except Pb were below the suggested thresholds linked with adverse health effects in birds. The hepatic Pb found in almost half of all the tested individuals suggests environmental exposure to this toxic element. One of the tested white-tailed eagles had hepatic Pb above the threshold of sublethal poisoning. Although our results seem optimistic, as previous Polish studies showed a higher prevalence of birds with hepatic Pb exceeding the toxicity threshold, they indicate that exposure to this toxic metal could still pose an additional threat to the health of Polish eagles.

Contrasting patterns of genetic diversity and lack of population structure in the lesser spotted eagle Clanga pomarina (Aves: Accipitriformes) across its breeding range

Characterizing the genetic diversity and structure of populations is essential for understanding their evolutionary history and planning species conservation. The lesser spotted eagle (Clanga pomarina) is a large migratory raptor with a relatively small breeding range concentrated in Eastern Europe. We evaluated the level of genetic diversity and population structuring by estimating the length diversity of 23 microsatellite markers in 306 individuals and sequencing 473 nucleotides from the mitochondrial pseudo-control region in 265 individuals across the distribution range. The microsatellite data suggested shallow differentiation between geographical regions and moderate genetic diversity across the range; no recent population bottlenecks were detected. Mitochondrial diversity was relatively low; however, high values were recorded at the southern edge of the distribution range. This, in combination with the star-like distribution of mitochondrial haplotypes, suggests the expansion of the European population from a single (Balkan) refugium during the late Pleistocene or early Holocene after the glacial population bottleneck. However, the Caucasian population may have survived in a separate refugium. We conclude that the lack of clear population structuring and ongoing gene flow across Europe support the treatment of the geographically restricted global population of the lesser spotted eagle as a single evolutionary and conservation unit.

The genome sequence of the European golden eagle, Aquila chrysaetos chrysaetos Linnaeus 1758

We present a genome assembly from an individual female Aquila chrysaetos chrysaetos (the European golden eagle; Chordata; Aves; Accipitridae). The genome sequence is 1.23 gigabases in span. The majority of the assembly is scaffolded into 28 chromosomal pseudomolecules, including the W and Z sex chromosomes.

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.

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.

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.

Only true pelagics mix: comparative phylogeography of deepwater bathybatine cichlids from Lake Tanganyika

In the absence of dispersal barriers, species with great dispersal ability are expected to show little, if at all, phylogeographic structure. The East African Great Lakes and their diverse fish faunas provide opportunities to test this hypothesis in pelagic fishes, which are presumed to be highly mobile and unrestricted in their movement by physical barriers. Here, we address the link between panmixis and pelagic habitat use by comparing the phylogeographic structure among four deepwater cichlid species of the tribe Bathybatini from Lake Tanganyika. We show that the mitochondrial genealogies (based on the most variable part or the control region) of the four species are very shallow (0.8-4% intraspecific divergence across entire distribution ranges) and that all species experienced recent population growth. A lack of phylogeographic structure in the two eupelagic species, Bathybates fasciatus and B. leo, was consistent with expectations and with findings in other pelagic cichlid species. Contrary to expectations, a clear phylogeographic structure was detected in the two benthopelagic species, B. graueri and Hemibates stenosoma. Differences in genetic diversity between eupelagic and benthopelagic species may be due to differences in their dispersal propensity, mediated by their respective predatory niches, rather than precipitated by external barriers to dispersal.

Nuclear and Mitochondrial DNA Analyses of Golden Eagles (Aquila chrysaetos canadensis) from Three Areas in Western North America; Initial Results and Conservation Implications

Understanding the genetics of a population is a critical component of developing conservation strategies. We used archived tissue samples from golden eagles (Aquila chrysaetos canadensis) in three geographic regions of western North America to conduct a preliminary study of the genetics of the North American subspecies, and to provide data for United States Fish and Wildlife Service (USFWS) decision-making for golden eagle management. We used a combination of mitochondrial DNA (mtDNA) D-loop sequences and 16 nuclear DNA (nDNA) microsatellite loci to investigate the extent of gene flow among our sampling areas in Idaho, California and Alaska and to determine if we could distinguish birds from the different geographic regions based on their genetic profiles. Our results indicate high genetic diversity, low genetic structure and high connectivity. Nuclear DNA Fst values between Idaho and California were low but significantly different from zero (0.026). Bayesian clustering methods indicated a single population, and we were unable to distinguish summer breeding residents from different regions. Results of the mtDNA AMOVA showed that most of the haplotype variation (97%) was within the geographic populations while 3% variation was partitioned among them. One haplotype was common to all three areas. One region-specific haplotype was detected in California and one in Idaho, but additional sampling is required to determine if these haplotypes are unique to those geographic areas or a sampling artifact. We discuss potential sources of the high gene flow for this species including natal and breeding dispersal, floaters, and changes in migratory behavior as a result of environmental factors such as climate change and habitat alteration. Our preliminary findings can help inform the USFWS in development of golden eagle management strategies and provide a basis for additional research into the complex dynamics of the North American subspecies.