Genetic differentiation in an endangered and strongly philopatric, migrant shorebird

Background Populations living in fragmented habitats may suffer from loss of genetic variation and reduced between-patch dispersal, which are processes that can result in genetic differentiation. This occurs frequently in species with reduced mobility, whereas genetic differentiation is less common among mobile species such as migratory birds. The high dispersal capacity in the latter species usually allows for gene flow even in fragmented landscapes. However, strongly philopatric behaviour can reinforce relative isolation and the degree of genetic differentiation. The Southern Dunlin (Calidris alpina schinzii) is a philopatric, long-distance migratory shorebird and shows reduced dispersal between isolated breeding patches. The endangered population of the Southern Dunlin breeding at the Baltic Sea has suffered from habitat deterioration and fragmentation of coastal meadows. We sampled DNA across the entire population and used 12 polymorphic microsatellite loci to examine whether the environmental changes have resulted in genetic structuring and loss of variation. Results We found a pattern of isolation-by-distance across the whole Baltic population and genetic differentiation between local populations, even within the southern Baltic. Observed heterozygosity was lower than expected throughout the range and internal relatedness values were positive indicating inbreeding. Conclusions Our results provide long-term, empirical evidence for the theoretically expected links between habitat fragmentation, population subdivision, and gene flow. They also demonstrate a rare case of genetic differentiation between populations of a long-distance migratory species. The Baltic Southern Dunlin differs from many related shorebird species that show near panmixia, reflecting its philopatric life history and the reduced connectivity of its breeding patches. The results have important implications as they suggest that reduced connectivity of breeding habitats can threaten even long-distance migrants if they show strong philopatry during breeding. The Baltic Southern Dunlin warrants urgent conservation efforts that increase functional connectivity and gene flow between breeding areas. Supplementary Information The online version contains supplementary material available at 10.1186/s12862-021-01855-0.

depression, and in the long-term genetic stochasticity due to drift, poses a substantial threat in isolated and rapidly declining populations (e.g. [1]), as shown by many case studies (e.g. [6][7][8]). For such populations, gene flow and connectivity to other populations are crucial for maintaining genetic variation [9,10].
The effects of a small population size may be reinforced by life-history characteristics such as poor dispersal ability or philopatry. Strong natal and adult philopatry suggests that there must be benefits from returning to the natal or previous breeding site [11]. However, philopatry can be a detrimental strategy at the population level, when fragmentation-driven reduction in connectivity has reduced gene flow [12]. Strong philopatry is expected to lead to genetic structuring, isolation-by-distance (IBD) and increased inbreeding due to decreased opportunities for individuals to mate outside of kin [2]. Multiple examples of these processes exist in species with reduced mobility, but they are much less common in birds, especially in long-distance migratory species, because of their higher dispersal rates [13,14].
The Southern Dunlin (Calidris alpina schinzii) is a small, migratory shorebird species breeding exclusively in short-vegetated and wet grasslands [15,16]. The Baltic population, which breeds mainly on coastal grasslands, has suffered a dramatic decline (about 60% during the last 20 years) and is one of the most endangered shorebird populations in Europe [17,18] despite the species-level assessment of "Least Concern" [18]. In the beginning of the twentieth century, the Southern Dunlin was common and widespread in most parts of the Baltic [19]. By the start of the twenty-first century, the population size of the Baltic Southern Dunlin had declined to 1110-1360 breeding pairs [17]. Since then, several local populations have gone extinct, and the number of pairs is currently closer to 500 pairs [18]. While the populations have suffered from high nest predation pressure across the range [16,20,21] and possibly decreased adult survival [22], the initial reasons for the decline were likely large-scale agricultural changes and eutrophication which led to overgrowth, habitat loss, and fragmentation of previously connected meadow systems [18]. It is important to note that the population declines have continued despite breeding habitats being available in many areas around the Baltic (see [7]).
Habitat fragmentation can be expected to lead to genetic effects in the Southern Dunlin because both adults [23][24][25] and juveniles [26] are highly site-faithful to their breeding and natal sites. Yet, an analysis using data collected mostly before the drastic declines occurring after the turn of the twenty-first century, detected no genetic structuring or signs of genetic impoverishment in Baltic Southern Dunlins [27]. However, because reduced structural connectivity of patches has been shown to lead to decreased between-patch movements in this species [26], the current situation differs from the historical situation when breeding sites of the Southern Dunlin around the Baltic Sea were better connected and environmental predictability, favoring philopatry, was probably higher. Therefore, habitat fragmentation has since likely reduced the movements of individuals between populations. The resulting reduction in gene flow has already led to inbreeding, substantially increasing the likelihood of extinction of some populations [7]. Given these findings and the incessant decline of the Baltic population, we expected an overall reduction in genetic variation as well as increased genetic differentiation -even in this longdistance migratory species.
Building on extensive sampling from the entire Baltic Southern Dunlin population and genetic analyses based on polymorphic microsatellite markers, we (A) examine whether the populations at the Baltic are genetically differentiated, and (B) estimate levels of genetic variation of the Baltic populations in order to assess if the observed decline, increased isolation, and lack of connectivity of breeding sites have had an effect since previous studies. Importantly, the genetic connectedness of the Bothnian Bay population at the northernmost location of the Baltic Sea in Finland [28] to the schinzii populations in the southern Baltic and the alpina population in Lapland has never been studied. Therefore, we further (C) examine whether there is indication of gene flow between the nominate subspecies C. a. alpina and Baltic Southern Dunlin populations with special emphasis on the population at Bothnian Bay. Finally, we (D) discuss the evolutionary implications of our study and suggest relevant conservation measures.

Genetic diversity
Genotyping error rate was low, with a mean of 1.6% when genotyped twice. When the ambiguous samples were again genotyped twice, the error rate dropped to 0.1% (only one locus in four individuals remained unclearthis data were excluded from the analyses). In two populations, the presence of null alleles was suggested for loci CAS23 and Cme1 (see Additional file 1: Table S1 for information on the loci used), and these loci were also suspected to show stuttering in other populations. For CAS23, however, this was likely caused by the marker being sex-linked. Since the error rate was low and the suggested genotyping errors were not consistent among populations, all loci were used in the analyses. Linkage disequilibrium was present only in a few locus pairs, seemingly randomly in different populations, suggesting no strong linkage between the loci.
Allelic richness (A), corrected for the difference in sample size, was similar in all schinzii populations (varying from 3.61 in Pori to 3.76 in Jurmo; Table 1). The highest A (3.93) was found in the alpina subspecies. Observed heterozygosity was lower than expected (i.e. F IS values were positive) in every population (Table 1). F IS was significant in Pori and Jurmo, and especially in Estonia. The lowest F IS was in Bothnian Bay (0.001). The mean internal relatedness (IR) values per population were highest in Pori (0.231), Denmark (0.133) and Estonia (0.093), and lowest in Bothnian Bay (0.020; Table 1).

Isolation-by-distance
We found that the kinship coefficient slowly decreased with increasing distance between individuals, indicating an isolation-by-distance pattern ( Fig. 1) with regression coefficient per one km being -0.00002 (SE = 0.000004; p = 0.021, number of individuals = 374). The intra-group (IG) class and the first distance class (mean distance 405 km) showed significantly positive kinship coefficients (p = 0.000 and p = 0.021, respectively), whereas the third and fourth distance classes (mean distances 716 and 1129 km, respectively) were negative, with the last class significantly so (p = 0.028).

Genetic structure
The clustering analysis with program Structure resulted in the highest ΔK value at K = 4 (Additional file 2: Figure S1 and Table S2). Birds from the alpina population formed their own cluster ( (Table 2), the overall among-population fixation index being highly significant (F ST = 0.011, p < 0.001). The Bothnian Bay population differed from all Swedish and the Danish population, the Danish population from western and southern Sweden, and alpina from every other population except Pori, Jurmo and Estonia, which did not differ significantly from any population. The largest pairwise differences arose in pairwise D est comparisons: they were up to an order of magnitude higher than, but still concordant with the F ST values (R 2 = 0.747, p < 0.001), varying from − 0.007 to 0.334 (Table 2).

Discussion
We show that the Southern Dunlin populations breeding on the Baltic coastal meadows are genetically differentiated, providing a rare example of genetic structuring on a regional scale of a long-distance migratory species. This finding is in strong contrast to observations from many other northern shorebird species, which show limited genetic differentiation and high levels of gene flow on a geographic scale similar to our study [30] and even across much larger spatial scales (e.g. [31][32][33][34][35]). Genetic differentiation has mainly been found in island populations of otherwise panmictic species (e.g. [36], but see also [37]). Indeed, the Southern Dunlin is a habitat specialist, and the fragmented habitat patches resemble islands amidst a matrix of unsuitable breeding habitat. Accordingly, genetic differentiation is partly linked to the distance between breeding sites, as indicated by the significant IBD pattern. While these results are in line with the continued fragmentation of its breeding range, the philopatric life history of the species may be the process that prevents movement of this otherwise mobile species and thus further reduces functional connectivity of the populations, eventually affecting allele frequencies and resulting into genetic differentiation.
The IBD pattern showed a decrease of genetic similarity with increasing geographic distance across the Baltic region. In general, individuals breeding up to 400 km away have significantly positive kinship coefficients (i.e. they are genetically more similar to each other than if chosen by random), which is consistent with known longdistance movements between Pori and Bothnian Bay [28]. Genetic similarity becomes negative (i.e. individuals are genetically less similar than if chosen by random) when the distance extends over 600 km. Accordingly, the highest F ST values were observed between Bothnian Bay and the southwestern populations (Denmark and the Swedish populations) that are furthest apart. The Estonian population, which is located somewhat in the middle between the Bothnian Bay and the southwestern populations, showed mixed ancestry in the Structure analysis, mostly with the Bothnian Bay and western Swedish populations. This indicates that the Estonian population receives immigrants from other parts of the range. DAPC suggested the best number of clusters would be the same four as suggested by the Structure analysis. However, the clustering was not as clearly in accordance with the geographic origins of the samples as in Structure. Based on the results of these two analyses, it seems that there is a gradual change in allele frequencies along the Baltic Sea from the southwestern areas through Estonia towards the Bothnian Bay, resembling the cline reported by Marthinsen et al. [38] but on a much smaller spatial scale. Importantly, we also found differentiation (both with fixation indices and clustering analyses) within a smaller scale of 200 to 500 km between the populations in the southern Baltic. This result fits earlier reports that show movements to be rare in these local populations [7,24,26]. Despite the relatively short distances to other populations in the southwestern Baltic, western Sweden does not receive immigrants [7]. Therefore, even though movement generally decreases with increasing distance between breeding sites [26], unknown environmental factors operating together with distance seem to create ecological barriers and differentiation also on smaller scales. It is, for example, possible that the breeding sites in western Sweden are located too far west from the migration route that follows the coast of the Baltic Sea [39].
We identified four genetic clusters within our data: three of them were formed by the Baltic Southern Dunlin individuals and the fourth by C. a. alpina individuals (with the exception of Danish Southern Dunlins that clustered with alpina at small K-values; Fig. 2). This contradicts earlier results by Marthinsen et al. [38], who could not distinguish between Dunlin subspecies with clustering analysis using seven microsatellites. Given our results, it seems unlikely that C. a. alpina mix with the Southern Dunlin populations as a result of migratory short-stopping, as indicated in other subspecies [40].
We found lower levels of heterozygosity than expected (i.e. positive F IS values) and positive estimates of internal relatedness throughout the range, which is consistent with inbreeding. Indeed, previous work has documented severe inbreeding in the very small population in western Sweden [7] and a significantly positive inbreeding coefficient in the island of Öland, Eastern Sweden [27]. In the present study, the F IS values of these populations were positive though not statistically significant (Table 1), which may be a sampling effect. On the other hand, the highly significant F IS estimated for Estonia, one of the largest local populations, may in fact be explained by within-population genetic structure (Wahlund effect), which is further supported by the mixed ancestry in Estonia and weak differentiation from the other populations. The Bothnian Bay population showed the lowest F IS and IR values and highest observed heterozygosity. These attributes are consistent with annually occurring immigration to the Bothnian Bay [28] and relatively high movement rates between breeding patches [26].

Conclusions
Our results highlight that philopatry and reduced structural connectivity can result in population differentiation and IBD at a small spatial scale even in a mobile species. However, the pronounced site fidelity in the declining and fragmented Baltic Southern Dunlin population results in inbreeding ( [7], this study) and may lead into an evolutionary trap [7,41]. Indeed, the ongoing population decline suggests that the local populations across the Baltic are becoming smaller and more vulnerable to stochastic processes and eventually extinction. In addition to maintaining existing breeding sites, the Baltic Southern Dunlin therefore needs urgent conservation efforts that increase structural connectivity among present and potential breeding sites. For isolated and inbred populations such as the one in western Sweden, translocation of individuals might be the only way to ensure the exchange of genetic material. This could be achieved by captive breeding and introduction, perhaps from the genetically and geographically closest populations in eastern and southern Sweden.

Study species
The Southern Dunlin is one of several subspecies of Dunlin that breed across arctic and subarctic tundra, alpine wetlands, and wet grasslands in the temperate zone [43]. In addition to the Baltic region, Southern Dunlins breed in Iceland, British Isles, Faroe Islands and southeastern Greenland with an estimated 970 000-990 000 individuals [18,44]. Their main autumn migration route follows the Atlantic coast of Europe and continues south to the main wintering areas in northwest and northern Africa, where the wintering grounds are shared with the northern subspecies alpina and arctica [39,43,45]. Based on previous genetic analyses, the Southern Dunlin is part of a larger phylogenetic group including its geographically closest subspecies alpina and arctica. Even though these subspecies are genetically very similar, some genetic differences have been found. These differences mainly reflect a gradual change in allele frequencies throughout their breeding ranges [38,46,47].
The largest numbers of breeding Southern Dunlins in the Baltic area are found in Denmark and Estonia, where the populations are estimated to 170 and 180-230 pairs, respectively [18,42,48]. Sweden holds about 75 pairs [49] divided in three populations in western, eastern and southern Sweden, 250-300 km apart. Finland holds around 40-50 pairs (own observations). Most of them breed at Bothnian Bay; two other populations exist in southern Finland in Pori and Jurmo, over 400 km away (Fig. 4). The populations in Poland and Lithuania have disappeared and there is only one breeding site left in Germany [18].

Sampling
Blood samples of 344 adult Southern Dunlins were collected from eight populations around the Baltic Sea during the breeding seasons in 1997-2016 (Fig. 4, Additional file 3: Table S3). The length of the sampling period varied among the populations; ranging from only 1 year in some to 17 years in western Sweden (Additional file 3: Table S3). However, even the longest sampling periods represent only a few generations given that the Dunlin has an average longevity of 7 years [18].
The samples were obtained by puncturing the brachial vein. In addition, four eggs with embryos from four different deserted nests were obtained from the Danish population. The only sample from Gotland was combined with the eastern Swedish population due to their proximity (< 150 km).
Samples (feathers) of the northern subspecies alpina were obtained from 26 breeding adults from Finnish Lapland, collected during 2010-2017. The total number of sampled Dunlins was thus 374.

DNA extraction and microsatellite amplification
DNA was extracted from blood and tissue samples either with the standard phenol-chloroform method [50], or using the UltraClean ® BloodSpin ™ Kit or UltraClean ® Tissue & Cells DNA Isolation Kit (MoBio Laboratories), and from feathers using the method described in Rönkä et al. [34]. Individuals were genotyped for 12 microsatellite loci, which were amplified in 10 µl volumes containing 20-100 ng of template DNA, 0.1 µM of each primer, 0.8-1 mM MgCl 2 , 0.2 mM of dNTPs, 1 µl of 10 × PCR-Buffer and 0.l U of DNA-polymerase (Biotools). The amplification profile was 94 °C for 1 min followed by 35 cycles of 94 °C for 30 s, 52-58 °C for 45 s (see Additional file 1: Table S1), 72 °C for 45 s and 72 °C for 10 min for final extension. The PCR reactions were run with ABI 3730 sequencer using GS500-Liz size standard (Applied Biosystems) and the loci were scored with GeneMapper v. 4.0. (Applied Biosystems), except for the Swedish samples, which were scored with CEQ TM 8000 Genetic Analysis System (Beckman Coulter) using the Fragment Analysis Module v. 8.0.52. Due to possible differences between the allele sizes defined by the two sequencers, samples were calibrated by genotyping five Swedish individuals with both sequencers. Genotyping error rate was calculated by amplifying most individuals twice. If differences were found between the two runs, samples were genotyped twice more.

Basic statistics
The microsatellite data were checked for potential genotyping errors (stutter bands, null-alleles and large allele dropouts) with MicroChecker v. 2.2.3. [51] for all populations except Denmark, Pori and Jurmo, which contained too few samples. GenePop v. 4.2. [52,53] was used to test linkage disequilibrium (LD) and deviation from the Hardy-Weinberg (HW) equilibrium for each locus and population, as well as calculating population-wise F IS (i.e. deviation from random mating). FSTAT v. 2.9.3.2. [54] was used to calculate allelic richness (individuals with any missing data excluded and using rarefraction to account for different sample sizes) and an Excel macro [55] to calculate internal relatedness of individuals. Observed and expected heterozygosity was calculated with Arlequin v. 3.5.1.2. [56] for each population. One locus (CAS23) was found to be sex-linked as females had always only one allele, whereas males were often heterozygous. This locus was excluded from the calculations of heterozygosity, F IS and allelic richness. For the Structure run (see below), all females and the individuals whose sex remained undetermined were coded as both alleles missing regarding this locus.

Isolation-by-distance
Isolation-by-distance within the Southern Dunlin was tested with program SpaGeDi v. 1.4. [57]. The program calculates the genetic relatedness of all possible pairs of individuals, compares the relatedness to the corresponding geographic distances and tests whether there is a correlation. The number and intervals of the distance classes were decided after fine-tuning them so that they met the 'rule of thumb' of SpaGeDi (# pairs > 500, % partic > 50% and CV partic < 1 per distance class; see SpaGeDi manual for further information). This resulted in three distance classes: birds breeding up to 650 km, 651-850 km and 851-1400 km from each other, respectively. In addition, an intra-group (i.e. population-wise) Fig. 4 Locations of the sampled populations (in blue) of the Baltic Southern Dunlin (Calidris alpina schinzii) and Dunlin (C. a. alpina) with sample sizes (N) and the non-sampled (in red) breeding distribution of the Baltic Dunlin [18]. The breeding distribution was drawn according to [18] and Rönkä et al. (unpublished data). See text and [18] for details on the current distribution