Population fragmentation leads to spatial and temporal genetic structure in the endangered Spanish imperial eagle

The fragmentation of a population may have important consequences for population genetic diversity and structure due to the effects of genetic drift and reduced gene flow. We studied the genetic consequences of the fragmentation of the Spanish imperial eagle (Aquila adalberti) population into small patches through a temporal analysis. Thirty‐four museum individuals representing the population predating the fragmentation were analysed for a 345‐bp segment of the mitochondrial control region and a set of 10 nuclear microsatellite loci. Data from a previous study on the current population (N = 79) were re‐analysed for this subset of 10 microsatellite markers and results compared to those obtained from the historical sample. Three shared mitochondrial haplotypes were found in both populations, although fluctuations in haplotype frequencies and the occurrence of a fourth haplotype in the historical population resulted in lower current levels of haplotype and nucleotide diversity. However, microsatellite markers revealed undiminished levels of nuclear diversity. No evidence for genetic structure was observed for the historical Spanish imperial eagle population, suggesting that the current pattern of structure is the direct consequence of population fragmentation. Temporal fluctuations in mitochondrial and microsatellite allelic frequencies were found between the historical and the current population as well as for each pairwise comparison between historical and current Centro and historical and current Parque Nacional de Doñana nuclei. Our results indicate an ancestral panmictic situation for the species that management policies should aim to restore. A historical analysis like the one taken here provides the baseline upon which the relative role of recent drift in shaping current genetic patterns in endangered species can be evaluated and this knowledge is used to guide conservation actions.     

Eighteen polymorphic microsatellite markers for the highly endangered Spanish imperial eagle (Aquila adalberti) and related species

Here we describe the development of 18 polymorphic microsatellite markers for the endangered Spanish imperial eagle (Aquila adalberti). Microsatellites were tested in five other raptor species. These markers were revealed as good molecular tools for genetic population studies, individual identification and parentage assessment in Spanish imperial eagle and closely related species.     

Tetranucleotide microsatellites for aquila and haliaeetus eagles

A unique community of four syntopic eagle species exists in north‐central Kazakhstan. Questions about behaviour and genetics in these four species would benefit from the development of microsatellite markers. We isolated eight polymorphic microsatellite repeats (AAAG)_n from the eastern imperial eagle (Aquila heliaca) genome using a hybridization enrichment technique. These loci revealed moderate diversity in a local population of eastern imperial eagles (observed heterozygosity 0.26–0.78), and were also polymorphic in steppe eagles (A. nipalensis) and white‐tailed sea eagles (Haliaeetus albicilla). These primers may be polymorphic in other species of Aquila and Haliaeetus eagles.     

Behavioral evidence of hunting and foraging techniques by a top predator suggests the importance of scavenging for preadults

Scavenging may be a regular feeding behavior for some facultative raptor species occupying low quality habitats and/or with little experience in hunting techniques. However, its importance has been largely underestimated due to methodological limitations in identifying the real proportion in the diet. Here, through direct observations, we assessed the hunting and foraging success of the threatened Spanish imperial eagle Aquila adalberti determining the influence of age, sex, breeding status, habitat quality, prey type, and landscape characteristics. From 465 observations, Spanish imperial eagles used hunting in flight (42%), scavenging (30%), hunting from a perch (16%) and kleptoparasitism (12%). Our model suggests that Prey size and Prey type best explain hunting success, followed by Landscape and Sex. Our findings suggest that Spanish imperial eagles increase hunting success with age, with scavenging and kleptoparasitism regularly used as juveniles. The absence of relationships with any of the variables considered suggests that kleptoparasitism is an opportunistic behavior used sporadically. Scavenging is also independent of habitat quality and landscape characteristics. Accordingly, low prey density is not a driver of carrion use for preadult individuals, suggesting that a lack of hunting ability obliges this age‐class to use this alternative feeding technique regularly. As a result, the threatened Spanish imperial eagle population is also prone to mortality related to the illegal use of poison baits and, potentially, veterinary drugs (i.e., diclofenac).     

Genetics at the verge of extinction: insights from the Iberian lynx

Population viability might become compromised by the loss of genetic diversity and the accumulation of inbreeding resulting from population decline and fragmentation. The Iberian lynx (Lynx pardinus) provides a paradigmatic example of a species at the verge of extinction, and because of the well‐documented and different demographic histories of the two remaining populations (Doñana and Andújar), it provides the opportunity to evaluate the performance of analytical methods commonly applied to recently declined populations. We used mitochondrial sequences and 36 microsatellite markers to evaluate the current genetic status of the species and to assess the genetic signatures of its past history. Mitochondrial diversity was extremely low with only two haplotypes, alternatively fixed in each population. Both remnant populations have low levels of genetic diversity at microsatellite markers, particularly the population from Doñana, and genetic differentiation between the two populations is high. Bayesian coalescent‐based methods suggest an earlier decline starting hundreds of years ago, while heterozygosity excess and M‐ratio tests did not provide conclusive and consistent evidence for recent bottlenecks. Also, a model of gene flow received overwhelming support over a model of pure drift. Results that are in conflict with the known recent demography of the species call for caution in the use of these methods, especially when no information on previous demographic history is available. Overall, our results suggest that current genetic patterns in the Iberian lynx are mainly the result of its recent decline and fragmentation and alerts on possible genetic risks for its persistence. Conservation strategies should explicitly consider this threat and incorporate an integrated genetic management of wild, captive and re‐introduced populations, including genetic restoration through translocations.     

Primary moult of continental Black-tailed Godwits Limosa limosa limosa in the Doñana wetlands,Spain

Capsule: Most Continental Black-tailed Godwits Limosa limosa limosa using the Doñana wetlands during post-breeding migration appear to begin moult before they arrive and suspend moult before they migrate onwards to West Africa.

Aims: We aim to describe the primary moult strategies and patterns in the Continental Black-tailed Godwits using the Doñana wetlands, a major passage and wintering area for waterbirds in southern Spain.

Methods: Individual godwits were captured, marked and their primary moult was scored in Doñana during the non-breeding season (June–March) in 2011 and 2012. Data from resightings of colour-marked godwits and birds equipped with satellite transmitters were used to estimate stopover duration during post-breeding migration (June–September) to determine if godwits move to West Africa before completing their primary moult.

Results: Average primary moult duration was estimated to be 84 days?±?9 se, during 29 June–21 September and did not differ between sexes. Only 2% of individuals were observed with suspended moult. We estimated stopover duration in Doñana to be 13 days?±?2 se before migrating to West Africa.

Conclusions: Most godwits stage for about two weeks in the Doñana wetlands during southward migration, moult their primaries and appear to suspend moult before crossing the Sahara. Others may complete their primary moult in Doñana, or elsewhere in Europe and overwinter in Doñana where increasing numbers of godwits have been detected in recent years. A few individuals may finish the moult in Doñana and migrate to West Africa late in the post-breeding season.     

NADH diaphorase polymorphism in European fallow deer

NADH diaphorase polymorphism is reported for European fallow deer (D. d. dama) from Coto de Doñana National Park (Spain). Allele frequencies ofp=0.682 (anodic variant) andq=0.318 (cathodic variant) in 110 specimens suggest population genetic usefulness of this biallelic system in an otherwise fairly monomorphic cervid species.     

Solving man-induced large-scale conservation problems: the Spanish imperial eagle and power lines

Background

Man-induced mortality of birds caused by electrocution with poorly-designed pylons and power lines has been reported to be an important mortality factor that could become a major cause of population decline of one of the world rarest raptors, the Spanish imperial eagle (Aquila adalberti). Consequently it has resulted in an increasing awareness of this problem amongst land managers and the public at large, as well as increased research into the distribution of electrocution events and likely mitigation measures.

Methodology/Principal Findings

We provide information of how mitigation measures implemented on a regional level under the conservation program of the Spanish imperial eagle have resulted in a positive shift of demographic trends in Spain. A 35 years temporal data set (1974–2009) on mortality of Spanish imperial eagle was recorded, including population censuses, and data on electrocution and non-electrocution of birds. Additional information was obtained from 32 radio-tracked young eagles and specific field surveys. Data were divided into two periods, before and after the approval of a regional regulation of power line design in 1990 which established mandatory rules aimed at minimizing or eliminating the negative impacts of power lines facilities on avian populations. Our results show how population size and the average annual percentage of population change have increased between the two periods, whereas the number of electrocuted birds has been reduced in spite of the continuous growing of the wiring network.

Conclusions

Our results demonstrate that solving bird electrocution is an affordable problem if political interest is shown and financial investment is made. The combination of an adequate spatial planning with a sustainable development of human infrastructures will contribute positively to the conservation of the Spanish imperial eagle and may underpin population growth and range expansion, with positive side effects on other endangered species.     

Population genetic inferences using immune gene SNPs mirror patterns inferred by microsatellites

Single nucleotide polymorphisms (SNPs) are replacing microsatellites for population genetic analyses, but it is not apparent how many SNPs are needed or how well SNPs correlate with microsatellites. We used data from the gopher tortoise, Gopherus polyphemus—a species with small populations, to compare SNPs and microsatellites to estimate population genetic parameters. Specifically, we compared one SNP data set (16 tortoises from four populations sequenced at 17 901 SNPs) to two microsatellite data sets, a full data set of 101 tortoises and a partial data set of 16 tortoises previously genotyped at 10 microsatellites. For the full microsatellite data set, observed heterozygosity, expected heterozygosity and F_ST were correlated between SNPs and microsatellites; however, allelic richness was not. The same was true for the partial microsatellite data set, except that allelic richness, but not observed heterozygosity, was correlated. The number of clusters estimated by structure differed for each data set (SNPs = 2; partial microsatellite = 3; full microsatellite = 4). Principle component analyses (PCA) showed four clusters for all data sets. More than 800 SNPs were needed to correlate with allelic richness, observed heterozygosity and expected heterozygosity, but only 100 were needed for F_ST. The number of SNPs typically obtained from next‐generation sequencing (NGS) far exceeds the number needed to correlate with microsatellite parameter estimates. Our study illustrates that diversity, F_ST and PCA results from microsatellites can mirror those obtained with SNPs. These results may be generally applicable to small populations, a defining feature of endangered and threatened species, because theory predicts that genetic drift will tend to outweigh selection in small populations.     

Rabbit (Oryctolagus cuniculus) abundance and protected areas in central-southern Spain: why they do not match?

European rabbits (Oryctolagus cuniculus), a keystone species in the Iberian Mediterranean ecosystem, are the staple prey of the Iberian lynx (Lynx pardinus) and the Spanish imperial eagle (Aquila adalberti). These predators require medium to high rabbit densities and a low degree of human disturbance. We compared rabbit abundances in areas of central-southern Spain under three levels of protection and management: protected areas, intensively managed (nonprotected) hunting estates, and other nonprotected areas. We used pellet abundance indices to estimate rabbit density in 118 surveys conducted during the summers of 2002 and 2003. We observed greater rabbit abundance in intensively managed hunting estates compared to protected areas and other nonprotected areas, perhaps because policy makers did not consider rabbit numbers when selecting priority areas. Alternatively, differences in game management practices (e.g., predator control or habitat management) may explain the higher rabbit densities observed in managed hunting estates. Our results suggest that the best feeding conditions for the Iberian lynx and the Spanish imperial eagle occur in intensively managed hunting areas, where such predators are frequently persecuted. The conservation of these endangered predators may require efforts to increase rabbit densities in protected areas.     

Modelling habitat use and distribution of golden eagles <Emphasis Type="Italic">Aquila chrysaetos</Emphasis> in a low-density area of the Iberian Peninsula

We analyse the current situation of the Golden eagle (Aquila chrysaetos) in the region of Galicia in NW Spain. At present, the entire Galician population (five pairs) is located within an area of about 2000 km² in the province of Ourense. To identify high-priority areas for golden eagle conservation, we derived predictive models of habitat suitability using logistic regression and a Geographic Information System (GIS). Specifically, to model the distribution of the breeding population we considered topographic features, land use and degree of humanization, using a 10 × 10 km grid. Presence/absence of golden eagle nests was used as the dependent variable; analyses were performed both considering current nesting areas and considering old nesting areas (1960s and 70s). At the spatial scale considered, the best predictors of habitat suitability for breeding were topographical variables indicative of rugged relief. For current nesting areas the most parsimonious model included maximum altitude. We consider that the predictive models obtained may be of use for the monitoring and conservation management of the golden eagle population in this region. Conservation problems associated with habitat constraints such as food supply, availability of nesting sites, changes in land use and human disturbance are discussed.     

Novel tetranucleotide microsatellite DNA markers for the wood frog,Rana sylvatica

Fifteen tetranucleotide microsatellite loci were identified and characterized for wood frogs (Rana sylvatica) collected from three vernal pools in the southeastern US. These markers revealed a high degree of genetic diversity (nine to 34 alleles per locus), heterozygosity (30.6–92.3%) and allelic heterogeneity (69% of comparisons were statistically significant). Considerable differentiation among populations was observed as genetic distances (chord) ranged between 0.40 and 0.55 and all F_ST values (0.02–0.05) were statistically significant. Genotypic assignment tests correctly classified 103 of 113 individuals to their respective collection. These markers should prove useful for investigating fine‐scale population structure and metapopulation dynamics.     

Evaluating mortality rates and causalities in a critically endangered felid across its whole distribution range

The conservation of endangered species requires accurate data, and knowledge of cause-specific mortality rates is one of the most important issues. In recent years, conservation programs for the critically endangered Iberian lynx Lynx pardinus have been developed on the basis of mortality data derived 30 years ago from the small Doñana population. Thus, there is an urgent need for an update of mortality rates and causes in both populations (Sierra Morena and Doñana). Here we use radio-tracking information from the whole range of the Iberian lynx to quantify mortality rates and identify their causes. Between 2006 and 2011, we radio-tagged 78 Iberian lynxes from its two remaining populations (39 from Sierra Morena and 39 from Doñana). Mortality events were evaluated to identify causes, and cause-specific annual mortality rates (AMR) were obtained using the nonparametric cumulative incidence function estimator. Overall, AMR was estimated at 0.16?±?0.05 (0.19?±?0.09 in Sierra Morena and 0.12?±?0.07 in Doñana). Disease was the main cause of mortality both for the whole population and the Doñana population. Poaching was the main cause of mortality in Sierra Morena. Our results suggest that the best strategy for conserving this species is to focus action on decreasing the fatal effect of disease and poaching. Given the possible existence of an underlying inbreeding-mediated immunosuppression, genetic management aimed at increasing the genetic diversity of this population is also recommended.     

GENETIC STRUCTURE AND MALE-MEDIATED GENE FLOW IN THE GHOST BAT (MACRODERMA GIGAS)

The Australian ghost bat is a large, opportunistic carnivorous species that has undergone a marked range contraction toward more mesic, tropical sites over the past century. Comparison of mitochondrial DNA (mtDNA) control region sequences and six nuclear microsatellite loci in 217 ghost bats from nine populations across subtropical and tropical Australia revealed strong population subdivision (mtDNA φ_ST = 0.80; microsatellites UR_ST = 0.337). Low-latitude (tropical) populations had higher heterozygosity and less marked phylogeographic structure and lower subdivision among sites within regions (within Northern Territory [NT] and within North Queensland [NQ]) than did populations at higher latitudes (subtropical sites; central Queensland [CQ]), although sampling of geographically proximal breeding sites is unavoidably restricted for the latter. Gene flow among populations within each of the northern regions appears to be male biased in that the difference in population subdivision for mtDNA and microsatellites (NT φ_ST = 0.39, UR_ST = 0.02; NQ φ_ST = 0.60, UR_ST = ?0.03) is greater than expected from differences in the effective population size of haploid versus diploid loci. The high level of population subdivision across the range of the ghost bat contrasts with evidence for high gene flow in other chiropteran species and may be due to narrow physiological tolerances and consequent limited availability of roosts for ghost bats, particularly across the subtropical and relatively arid regions. This observation is consistent with the hypothesis that the contraction of the species' range is associated with late Holocene climate change. The extreme isolation among higher-latitude populations may predispose them to additional local extinctions if the processes responsible for the range contraction continue to operate.     

Isolation and characterization of novel tetranucleotide microsatellite DNA markers for the spotted salamander,Ambystoma maculatum

Fifteen tetranucleotide microsatellite loci were identified and characterized for spotted salamanders (Ambystoma maculatum) collected from three vernal pools in the south‐eastern USA. These markers revealed a high degree of genetic diversity (7–32 alleles per locus), heterozygosity (31.6–86.3%) and allelic heterogeneity (91% of comparisons were statistically significant). Considerable differentiation among populations was observed as genetic distances (chord) ranged between 0.50 and 0.65, and all F_ST values (0.08–0.14) were statistically significant. Moreover, genotypic assignment tests correctly classified all individuals to their respective collection. These markers should prove useful for investigating fine‐scale population structure and mating system.     

Disentangling the impact of demographic factors on population differentiation of an endangered freshwater crayfish (Austropotamobius pallipes) using population density and microsatellite data

1. Habitat fragmentation of stream ecosystems often results in decreased connectivity between populations and lower population sizes. Hence, understanding how habitat fragmentation affects genetic erosion is important for the preservation of freshwater biodiversity, in particular, as small populations suffer from loss of genetic diversity through genetic drift and loss of fitness because of inbreeding, increasing the risk of extinction. 2. Here, we assess the impact of demographic factors on population differentiation in the endangered freshwater crayfish Austropotamobius pallipes by analysing population genetic structure, estimating effective population sizes and comparing levels of polymorphism at five microsatellite loci with density estimates of 10 populations within a small French catchment that has become progressively confined to headwaters over the last six decades. 3. Levels of expected heterozygosity and allelic richness per population were relatively low (0.214–0.396 and 1.6–2.6, respectively). We found strong genetic differentiation between these geographically close populations (F_ST = 0.283), with weak statistical evidence for a pattern of isolation by distance. Estimates of effective population size were low (<150) in most populations, but potentially reached several thousands in three populations. 4. Population density and allelic richness were strongly positively correlated. A robust relationship between population density and heterozygosity values was also noted, but only after discarding two populations for which significant genetic signatures of a recent bottleneck were found; these two populations displayed high expected heterozygosity compared with a very low density. Populations with the highest densities of individuals had the highest effective population size estimates and vice versa. 5. Our results clearly show the importance of demographic factors and genetic drift on A. pallipes populations. Furthermore, analysis of genetic and population density data is a pragmatic and efficient approach to corroborate inferences from genetic data and can be particularly useful in the identification of populations experiencing a bottleneck and therefore in conservation genetics studies aiming at identifying priority populations for conservation.     

Ongoing anthropogenic eutrophication of the catchment area threatens the Doñana World Heritage Site (South-west Spain)

Eutrophication is a major cause of wetland degradation worldwide. In recent decades, reductions in nutrient inputs have led to improvements in water quality in many rivers and lakes in central and northern Europe, but long-term trends are less clear in southern Europe. We conducted the first comprehensive study of water quality in Doñana (SW Spain), one of the most important wetland complexes in Europe and the Mediterranean region. The core area of Doñana is a large shallow, seasonal marsh (UNESCO World Heritage Site—WHS) that floods during rainy, cool winter months, then dries out during the summer. The marsh is fed by three main streams whose catchments are within a Biosphere Reserve but are impacted by greenhouses (for cultivating fruit, irrigated with groundwater), poorly treated urban wastewaters and tourism. From 2013 to 2016, we monitored nutrient (Total P, Total N, PO₄^?3, NH₄⁺, NO₃^?, NO₂^?) and phytoplankton chlorophyll-a (chla) concentrations in surface waters of the Doñana marsh and the three main streams. We quantified changes in greenhouse cover since 1995 using satellite images. Nutrient concentrations in streams were consistently higher than in the marsh, particularly in the Partido and Rocina streams that regularly reached concentrations equivalent to a “bad physico-chemical status” under the EU Water Framework Directive (WFD), and whose catchments suffered a fivefold expansion of greenhouses from 1995 to 2016. The Partido was the most polluted stream, and the most affected by wastewater effluents, and had particularly high concentrations of NH₄⁺ (Geometric Mean?=?0.3 mg L^?1) and NO₂^? (GM?=?0.52 mg L^?1) across seasons. Patterns in chla concentrations were less consistent, but streams (GM?=?6.78 µg L^?1) generally had higher concentrations than the marsh (GM?=?4.27 µg L^?1). Nutrient concentrations in spot samples within the marsh largely depended on a combination of evaporation (as revealed by higher stable isotope δ²H values in the water column) and spatial processes (concentrations increase close to stream entry points, where conductivity is lower). Anthropogenic nutrient pollution in entry streams is a serious problem in Doñana, with extensive stretches too toxic for fish. Reinforcement of policies aimed at reducing nutrient inputs to Doñana are urgently required to meet the biodiversity conservation and environmental objectives for the WHS and WFD, respectively. Paradoxically, the marsh is currently relied upon to purify the water entering from streams.

     

A non-invasive genetic evaluation of population size,natal philopatry,and roosting behavior of non-breeding eastern imperial eagles (<Emphasis Type="Italic">Aquila heliaca</Emphasis>) in central Asia

Roughly one-third of all eagle species are considered to be threatened or endangered, but the ecology of most eagles remains poorly understood. While the pronounced territoriality of breeding adults facilitates behavioral studies, the demography of non-breeding individuals (pre-adults and non-territorial floaters) is almost completely unknown. Traditionally, limited data on pre-adult and floater movement come from wing-tagging and/or telemetry studies. As an alternative to these methods we used genetic analyses of non-invasively collected feathers to investigate the population biology of non-breeding eastern imperial eagles (Aquila heliaca) in Kazakhstan. Microsatellite profiles of shed feathers indicate that eastern imperial eagles roost communally with other raptors. Furthermore, roosts are large and dynamic: 287 non-breeding eastern imperial eagles were detected in our sample, and a mark-recapture analysis estimated the total number to be 308 ± 8. The natal origins of these individuals were investigated by comparing their microsatellite profiles to those available for >90% of the eastern imperial eagle chicks hatched at the study site over the six previous breeding seasons. Only 4% of the individuals genetically matched a chick, suggesting that the reserve may serve as a critical refugium for pre-adults and itinerant floaters. Feathers have long been recognized as a suitable source of DNA, but few studies have used wide-scale, non-invasive collections of feathers (>1,000 samples) to investigate fundamental aspects of avian biology. Our research demonstrates that non-invasive genetic analyses of feathers can be used to evaluate population size, natal philopatry, and local movements of birds that are difficult to study using traditional means.     

Population Genetic Structure of Pacific White Shrimp (<Emphasis Type="Italic">Litopenaeus vannamei</Emphasis>) from Mexico to Panama: Microsatellite DNA Variation

Genetic variation and population structure of wild white shrimp (Litopenaeus vannamei) from 4 geographic locations from Mexico to Panama were investigated using 5 microsatellite DNA loci. The genetic diversity between populations was indicated by the mean number of alleles per locus and mean observed heterozygosity, which ranged from 7.4 to 8.6 and from 0.241 to 0.388, respectively. Significant departures from Hardy-Weinberg equilibrium were found at most locations at each locus, with the exception Guatemala at Pvan0013, and were caused by high heterozygote deficiencies. Genetic differences between localities were detected by pairwise comparison based on allelic and genotypic frequencies, with the exception of locus Pvan1003. Significant pairwise F _ST values between locations and total F _ST showed that the white shrimp population is structured into subpopulations. However, population differentiation does not follow an isolation-by-distance model. Knowledge of the genetic diversity and structure of L.vannamei populations will be of interest for aquaculture and fisheries management to utilize and preserve aquatic biodiversity.     

Genetic Structure of Wild Rice Oryza Glumaepatula Populations in Three Brazilian Biomes Using Microsatellite Markers

The existence of Oryza glumaepatula is threatened by devastation and, thus, the implementation of conservation strategies is extremely relevant. This study aimed to characterize the genetic variability and estimate population parameters of 30 O. glumaepatula populations from three Brazilian biomes using 10 microsatellite markers. The levels of allelic variability for the SSR loci presented a mean of 10.3 alleles per locus and a value of 0.10 for the average allelic frequency value. The expected total heterozygosity (H_e) ranged from 0.63 to 0.86. For the 30 populations tested, the mean observed (H_o) and expected heterozygosities (H_e) were 0.03 and 0.11within population, respectively, indicating an excess of homozygotes resulting from the preferentially self-pollinating reproduction habit. The estimated fixation index ( _IS ) was 0.79 that differed significantly from zero, indicating high inbreeding within each O. glumaepatula population. The total inbreeding of the species (_IT ) was 0.98 and the genetic diversity indexes among populations, _ST and _ST, were 0.85 and 0.90, respectively, indicating high genetic variability among them. Thus, especially for populations located in regions threatened with devastation, it is urgent that in situ preservation conditions should be created or that collections be made for ex situ preservation to prevent loss of the species genetic variability.