The Use of Genetics for the Management of a Recovering Population: Temporal Assessment of Migratory Peregrine Falcons in North America

PLOS ONE, Nov 2010

Background Our ability to monitor populations or species that were once threatened or endangered and in the process of recovery is enhanced by using genetic methods to assess overall population stability and size over time. This can be accomplished most directly by obtaining genetic measures from temporally-spaced samples that reflect the overall stability of the population as given by changes in genetic diversity levels (allelic richness and heterozygosity), degree of population differentiation (FST and DEST), and effective population size (Ne). The primary goal of any recovery effort is to produce a long-term self-sustaining population, and these genetic measures provide a metric by which we can gauge our progress and help make important management decisions. Methodology/Principal Findings The peregrine falcon in North America (Falco peregrinus tundrius and anatum) was delisted in 1994 and 1999, respectively, and its abundance will be monitored by the species Recovery Team every three years until 2015. Although the United States Fish and Wildlife Service makes a distinction between tundrius and anatum subspecies, our genetic results based on eleven microsatellite loci suggest limited differentiation that can be attributed to an isolation by distance relationship and warrant no delineation of these two subspecies in its northern latitudinal distribution from Alaska through Canada into Greenland. Using temporal samples collected at Padre Island, Texas during migration (seven temporal time periods between 1985–2007), no significant differences in genetic diversity or significant population differentiation in allele frequencies between time periods were observed and were indistinguishable from those obtained from tundrius/anatum breeding locations throughout their northern distribution. Estimates of harmonic mean Ne were variable and imprecise, but always greater than 500 when employing multiple temporal genetic methods. Conclusions/Significance These results, including those from simulations to assess the power of each method to estimate Ne, suggest a stable or growing population, which is consistent with ongoing field-based monitoring surveys. Therefore, historic and continuing efforts to prevent the extinction of the peregrine falcon in North America appear successful with no indication of recent decline, at least from the northern latitude range-wide perspective. The results also further highlight the importance of archiving samples and their use for continual assessment of population recovery and long-term viability.

A PDF file should load here. If you do not see its contents the file may be temporarily unavailable at the journal website or you do not have a PDF plug-in installed and enabled in your browser.

Alternatively, you can download the file locally and open with any standalone PDF reader:

The Use of Genetics for the Management of a Recovering Population: Temporal Assessment of Migratory Peregrine Falcons in North America

et al. (2010) The Use of Genetics for the Management of a Recovering Population: Temporal Assessment of Migratory Peregrine Falcons in North America. PLoS ONE 5(11): e14042. doi:10.1371/journal.pone.0014042 The Use of Genetics for the Management of a Recovering Population: Temporal Assessment of Migratory Peregrine Falcons in North America Jeff A. Johnson 0 Sandra L. Talbot 0 George K. Sage 0 Kurt K. Burnham 0 Joseph W. Brown 0 Tom L. Maechtle 0 William S. Seegar 0 Michael A. Yates 0 Bud Anderson 0 David P. Mindell 0 Robert C. Fleischer, Smithsonian Institution National Zoological Park, United States of America 0 1 Department of Biological Sciences, Institute of Applied Sciences, University of North Texas, Denton, Texas, United States of America, 2 Alaska Science Center, U.S. Geological Survey, Anchorage, Alaska, United States of America, 3 High Arctic Institute, Orion, Illinois, United States of America, 4 Department of Biological Sciences, University of Idaho, Moscow, Idaho, United States of America , 5 Sheridan , Wyoming, United States of America, 6 Edgewood Research Development and Engineering Center, Department of Army, Aberdeen Proving Ground, Maryland, United States of America, 7 Raptor Research Center, Boise State University , Minden , Nevada, United States of America, 8 Falcon Research Group, Bow, Washington, United States of America, 9 California Academy of Sciences , San Francisco, California , United States of America Background: Our ability to monitor populations or species that were once threatened or endangered and in the process of recovery is enhanced by using genetic methods to assess overall population stability and size over time. This can be accomplished most directly by obtaining genetic measures from temporally-spaced samples that reflect the overall stability of the population as given by changes in genetic diversity levels (allelic richness and heterozygosity), degree of population differentiation (FST and DEST), and effective population size (Ne). The primary goal of any recovery effort is to produce a longterm self-sustaining population, and these genetic measures provide a metric by which we can gauge our progress and help make important management decisions. Methodology/Principal Findings: The peregrine falcon in North America (Falco peregrinus tundrius and anatum) was delisted in 1994 and 1999, respectively, and its abundance will be monitored by the species Recovery Team every three years until 2015. Although the United States Fish and Wildlife Service makes a distinction between tundrius and anatum subspecies, our genetic results based on eleven microsatellite loci suggest limited differentiation that can be attributed to an isolation by distance relationship and warrant no delineation of these two subspecies in its northern latitudinal distribution from Alaska through Canada into Greenland. Using temporal samples collected at Padre Island, Texas during migration (seven temporal time periods between 1985-2007), no significant differences in genetic diversity or significant population differentiation in allele frequencies between time periods were observed and were indistinguishable from those obtained from tundrius/anatum breeding locations throughout their northern distribution. Estimates of harmonic mean Ne were variable and imprecise, but always greater than 500 when employing multiple temporal genetic methods. Conclusions/Significance: These results, including those from simulations to assess the power of each method to estimate Ne, suggest a stable or growing population, which is consistent with ongoing field-based monitoring surveys. Therefore, historic and continuing efforts to prevent the extinction of the peregrine falcon in North America appear successful with no indication of recent decline, at least from the northern latitude range-wide perspective. The results also further highlight the importance of archiving samples and their use for continual assessment of population recovery and long-term viability. - Funding: Funding for field work and sample collection was generously provided by the U.S. Armys Edgewood Research Development and Engineering Center, The Wolf Creek Charitable Foundation, The Grassland Charitable Foundation, Ruth and Brian Mutch, the North American Falconers Association, The Offield Family Foundation, and The Peregrine Fund. Funding for the genetic portion of this study was provided by The Peregrine Fund and the U. S. Geological Survey Alaska Science Center. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. Competing Interests: The authors have declared that no competing interests exist. In cases where populations or species have a recent history of decline followed by increase, the use of genetic data can be a powerful tool for monitoring progress in conservation efforts [14]. For example, estimates of the genetic diversity (allelic and heterozygosity), effective population size (Ne), gene flow or dispersal, and population admixture can provide information useful for making future management decisions to prevent further population decline and extinction. For species previously considered threatened or endangered, section 4(g)(1) of the Endangered Species Act requires the U.S. Fish and Wildlife Service in cooperation with the States to monitor a species for a minimum of five years after being removed from the List of Endangered and Threatened Wildlife and Plants to ensure they maintain nonthreatened status. The incorporation of genetic monitoring into such programs can provide information on the progress made in creating and maintaining a self-sustaining population, regardless if genetic measures were addressed in the original rulemaking or listing process [e.g., 57]. This approach can be particularly important with populations or species that have wide geographic distributions over challenging terrain (e.g., mountainous) where accurate demographic measures from the field are costly or difficult to obtain. The peregrine falcon (Falco peregrinus) provides an example species recovery plan that could benefit from ongoing genetic monitoring. Globally, the peregrine falcon consists of nineteen subspecies and is found on every continent with the exception of Antarctica [8,9]. In North America, three subspecies are currently recognized [10]. F. p. pealei is a year-round resident of the Pacific Northwest from the coasts of northern Washington and British Columbia extending to the Aleutian Islands in Alaska. F. p. tundrius breeds throughout the Arctic tundra of Alaska, Canada and western Greenland, and F. p. anatum breeds south of the tundra to northern Mexico, except in coastal areas in the Pacific Northwest. Both F. p. tundrius and anatum are migratory with tundrius wintering as far south as central Argentina and Chile in South America [11,12; see also 13]. Historic estimates of peregrine falcon abundance range from 400 to 500 pairs in Greenland [14], 1,000 to 3,500 [15] to 7,548 pairs [16] in the Arctic, and 7,000 to 10,000 total pairs in North America [17]. Following the late 1940s, many peregrine falcon populations suffered a steady decline due primarily to exposure from organochlorines, including DDT (1,1,1-trichloro-2,2-bis[pchlorophenyl]-ethane) and its principle metabolite DDE (1,1dichloro-2, 2-bis[p-chlorophenyl]-ethylene), which caused direct mortality or adversely affected their reproduction and egg production [1821]. By 1964, peregrine falcons nesting east of the Rocky Mountains south of the boreal forests in Canada (F. p. anatum) were essentially extirpated [22,23]. To the west in the Rocky Mountains, the number of peregrine falcons (F. p. anatum) was also significantly reduced with only 15 (29%) of 51 known historic nest sites occupied in 1964 [24], and by 1979 after a much more extensive survey of historic nest site locations, only 12 (7.5%) out of 160 were occupied covering an area of four million km2 [25]. Declines in the Arctic (F. p. tundrius) were less severe; however, it is estimated that their abundance was reduced approximately 50 to 60% by 1975 [20,26]. For the Peales falcon (F. p. pealei), abundance remained relatively stable during this time period [26], presumably due to their specialized diet feeding predominately on sea birds (e.g., alcids; [27,28] as opposed to terrestrial avian prey that were more likely to be exposed to DDT. In the past two decades, peregrine falcons in the U.S. and Canada, including Europe [8,29], have made a remarkable recovery due to the ban of DDT in 1969 and 1972 in Canada and the U.S., respectively, the inclusion of F. p. anatum and F. p. tundrius on the precursor of the federal Endangered Species List in 1970 (35 FR 16047), and extensive propagation and release efforts made by many conservation groups [30]. As a result, the species was delisted in 1994 (F. p. tundrius; 59 FR 50796), and 1999 (F. p. anatum; 64 FR 46541-46558) in the U. S., and reassigned Special Concern status in April 2007 in Canada. Current estimate of total breeding population size for both F. p. tundrius and F. p. anatum is between 4,300 and 10,400 [31]; considering immature and floater (non-breeders) individuals, the total population size could be between 40,000 to 50,000 [13]. These estimates are based on both direct and indirect counts, including projections made from available potential nesting sites. Because our current estimates of population size are imprecise, significant uncertainty exists for making decisions for management purposes. The primary aim of this study was to determine the overall stability and effective size of high latitude peregrine falcon populations in North American and Greenland, F. p. anatum and tundrius. This was accomplished using temporally segregated samples from both migratory and breeding peregrine falcons sampled during fall and spring migration through Padre Island, Texas over a 21 year period (,7 generations) and from southwestern Greenland over a 1214 year period (,4 generations), respectively. We assume that the sampled migratory individuals from Padre Island used in this study possess high latitude breeding distributions throughout Alaska, Canada and Greenland as supported by band recovery [31], satellite telemetry studies [12,32,33] and genetics (this study; see below), and therefore represent the northern peregrine falcon breeding population in North America. Our temporal sampling allowed the investigation of overall population stability by assessing allele frequency change over time and the estimation of Ne for both the migrant population and a focal breeding population in Greenland. These results are useful for conservation monitoring purposes and for making decisions that may influence future peregrine falcon population viability in North America. Sampling and DNA extraction All samples were obtained from birds caught and bled under government permits, and all birds were released after processing. A total of 292 peregrine falcons were sampled for genetic analyses during migration through Padre Island, Texas. Blood samples were taken during both autumn and spring migration for each of the following temporal subsets: 1985/86, 1988/89, and 2006/07, with additional samples from spring 2001 migration period (see Table 1 for sample sizes). All samples were kept frozen, and DNA extractions were performed using methods described elsewhere [34]. An additional 349 samples were obtained from the contemporary northern breeding distribution of peregrine falcons throughout Alaska, Canada and western Greenland, of which 168 samples were used in a previous study [35]. The samples collected from breeding territories were included in this study to verify the degree of population subdivision throughout their northern breeding distribution from which the Padre Island migrants likely originated. For analysis purposes and adequate sample size considerations, individuals were grouped into geographic sampling regions for each of the subspecies (see Fig. 1; Table 1). Two additional subspecies, F. p. cassini (n = 25) from South America and F. p. macropus (n = 15) from Australia, were included in the analyses for comparative purposes. Genotyping Eleven microsatellite loci originally developed for the peregrine falcon (Fp5, Fp13, Fp31, Fp46-1, Fp54, Fp79-4, Fp82-2, Fp86-2, Fp89, Fp92-1, Fp107; [36]) were used for the microsatellite analyses. All microsatellite loci were dinucleotide repeats, and protocols used for PCR amplification have been described elsewhere [34,35,37]. Genotypic data generated in different laboratories using different procedures were calibrated using a subset of samples (n$4) for each of the eleven microsatellite loci. No ambiguities were observed across all loci after calibration. Statistical analyses Genetic diversity. linkage equilibrium Microsatellite genotypes were tested for and departure from Hardy-Weinberg Greenland 2001-043 F. p. tundrius and anatum migrants (Padre Island, TX) equilibrium within each population at each locus using the computer program GDA [38]. Sequential Bonferroni corrections were applied to correct for multiple simultaneous comparisons [39]. Mean number of alleles per locus (A) and observed (Ho) and expected (He) heterozygosity values were calculated using GDA. Measures of allelic richness (AR) were calculated using the program FSTAT version [40]. AR estimates control for uneven sample sizes among populations [41]. Differences in microsatellite genetic diversity estimates between sample locations and time periods were tested for significance using a Wilcoxon signed-rank test. Measures of FIS and its significance for each sampled population was calculated using Fishers exact test within Genepop v. 4.0.10 ([42,43]; after adjusting the p-value to account for multiple simultaneous comparisons [39]. Population subdivision. The degree of population subdivision between sample locations and temporal sampling periods was investigated using the Bayesian method of Pritchard et al. [44] and Falush et al. [45], implemented in the program STRUCTURE version 2.1. The number of genetically distinct clusters (K), or populations, was identified based on allele frequencies across loci while minimizing linkage and violations to Hardy-Weinberg equilibrium. The most likely value of K is determined by comparing the likelihood of the data for different values of K. To determine the number of clusters, we also calculated the rate of change in the log probability of the data between successive K values (DK) plotted against K following Evanno et al. [46]. Analyses were performed with all samples and subspecies, including additional analysis without F. p. pealei, F. p. cassini and F. p. macropus samples to determine the influence the latter three subspecies have on the overall approximation of K. Calculations were conducted with a burn-in period of 105 iterations followed by an additional 106 iterations. Each simulation from K = 1 to 8 was performed four times using an ancestry model allowing admixture where individual a was inferred from the data for each cluster (alpha .1 means that most individuals are admixed; [45]), and a model of correlated allele frequencies that did not include prior information on sampling origin. Final results from STRUCTURE were visualized using the program DISTRUCT [47]. The degree of population subdivision was also explored as implemented in the software TESS [48]. This latter approach determines the number of groups similar to STRUCTURE, but differs by taking into account the spatial organization of individuals and incorporates a regularization procedure that helps facilitate the choice of K [4850]. The method implemented in TESS is also less influenced by Isolation by Distance compared to methods such as STRUCTURE when identifying the number of distinct clusters when clinal variation exists [50,51]. The MCMC algorithm was run under Figure 1. Peregrine falcon population sampling locations in North America. Samples sizes for each area are given in Table 1. doi:10.1371/journal.pone.0014042.g001 admixture model with interaction parameter Y = 0.7, with 10,000 burn-in and 50,000 sweeps. Twenty independent iterations were run for K = 27 and after identifying the value of K that produced the highest likelihood, this was run 100 times and the 20 highest likelihood runs for Kmax were averaged using CLUMPP version 1.1.2 [52] applying the Full Search algorithm and the G pairwise matrix similarity statistics. Estimates of genetic differentiation based on pairwise FST and Dest values were also obtained to further investigate the overall stability between sampling locations and temporal sampling periods. FST values were calculated following Weir & Cockerham [53] as implemented in ARLEQUIN version 3.11 [54], and Dest values [55] were calculated using SPADE [56] and bootstrap proportions for estimates of 95% confidence intervals (CI) were based on 1,000 permutations. The program Isolation By Distance Web Service (IBDWS), version 3.16 [57] was used to perform a Mantel test with 10,000 randomizations to examine the correlation between matrices of genetic distance (pairwise FST and Dest) and geographic distance of breeding sample locations of anatum and tundrius. A second set of analyses that included pealei samples were also performed. Geographic distance between sample locations was measured as the euclidean distance (km) using the ruler implemented in Google Earth version Effective population size. Genetic estimates of effective population size (Ne) were calculated using different methodological approaches to assess the robustness of our results. Two general approaches to estimating Ne were explored: 1) analyzing single time period population samples, and 2) analyzing multiple temporal samples from the same population. Although the temporal approach has been shown to outperform single-sample methods when assuming a closed population [1,58,59], recent analytical developments have improved both the precision and accuracy of Ne estimates from individual population samples [e.g., 60,61]. We used a method originally based on linkage disequilibrium (LD; [62]) to estimate contemporary Ne from individual sampling periods. The LD method included a bias correction [60] as implemented in the program LDNE [63], which has been shown to improve performance even with non-ideal populations (e.g. skewed sex ratios or non-random variance in reproductive success; [60]). Estimates of Ne were obtained for each of the spring migratory sampling periods collected from Padre Island, TX. We do not use this method to estimate Ne for any of the breeding populations because we cannot assume a closed population [64,65]. A jackknife method was used to obtain 95% confidence intervals (CI) on loci, and estimates were calculated assuming random mating and excluded all alleles #0.01 [63]. To estimate Ne based on multiple sampling periods for the migrant population from Padre Island, we employed three methods that are based on the premise that temporal variance in neutral genetic allele frequencies is inversely proportional to Ne due to the effects of genetic drift in the absence of migration and mutation [66,67]. The first method is based on the standardized variance of change in allele frequencies (Fk) between at least two sampling periods (equation 11, [68]; see also [69]). Because bias can exist with this method when estimates are based on small sample sizes and skewed allele frequencies [70], we used the weighing scheme of Jorde & Ryman [71] to provide an estimate of Ne with our dataset. Using sampling plan I [68], estimates of Ne and 95% confidence intervals were obtained using the program TEMPOFS [71]. Sampling plan I (i.e. nondestructive sampling) requires an estimate of population census size (N) to calculate Ne. Because we do not have a precise estimate of N for the peregrine falcon population, we calculated Ne using a range of values from 1,000 to 100,000 individuals to determine if uncertainty in N influences our estimate of Ne. Two additional estimates of Ne were obtained for the Padre Island migrant population using a coalescent-based method as employed in the program TM3 [72] and a pseudo-likelihood method implemented in the program MLNE 2.3 [65,73]. Both methods were used to calculate Ne of the migrant population while assuming an Ne-MAX of 10,000, no immigration, and a generation time of three years. To estimate Ne based on multiple temporal sampling periods for the Greenland population, we used MLNE 2.3 while assuming 1) a closed population (NeCLOSED), and 2) accounting for immigration from a potential source population (NeOPEN). Estimates of Ne from temporal data when mistakenly assuming a closed population are likely to be incorrect because, in addition to the effects of genetic drift, immigration will influence allele frequencies of the population to an extent that is related to the amount of differentiation between populations [59,64,65,74]. Therefore, when immigration is present, methods that account for this effect should be employed to generate accurate estimates of Ne. Using MLNE, two temporally spaced datasets from Greenland (1990 and 2001-04; four generations) were used to estimate both NeCLOSED and NeOPEN to assess the potential influence of immigration on our estimate of Ne and for comparative purposes with our estimates from the migrant Padre Island population Ne. Similar to the migrant dataset, 10,000 was used as our NeMAX. We used the pooled allele frequencies from contemporary tundrius and anatum breeding locations as the potential source population for immigrants into Greenland for estimating NeOPEN. Additional estimates were also calculated using the spring migrants from Padre Island, TX as the source population to evaluate the choice of source population on NeOPEN. Simulations. To assess the utility and precision of methods used to estimate Ne, we used simulated data representing multiple populations of specified size. This was done primarily to determine our ability to estimate Ne in populations of large size where drift is not likely to play a strong role influencing allele frequency change over a short time periods, e.g., seven generations (as with this study). Using the empirical data from spring 1986 as our initial sampling period (T0), we used the program BottleSim [75] to simulate seven generations (T7) at population sizes Ne = 50, 100, 200, 300, 500, 1000, 2000, and 5000. For each simulation based on 1000 iterations, we used the settings for maximum generation overlap (100%), random mating, three years for age at first breeding, 12-year longevity, and equal sex ratios [13]. Estimates of Ne were calculated similar to the empirical data, and their deviations from the specified Ne were then determined and directly compared to results obtained using the migrant Padre, TX temporal dataset. We also assessed levels of differentiation among ten populations from each of the simulated datasets of known size using similar samples sizes (n = 46). This was done to investigate the development of genetic differentiation relative to Ne after seven generations had passed similar to our empirical dataset. Genetic diversity measures Eight of the eleven microsatellite loci were polymorphic in all peregrine falcon sampling locations. Locus Fp5 was monomorphic in F. p. tundrius Northwest and Ungava Bay, F. p. anatum Alberta and Northeast, F. p. cassini, and Padre Island spring 2001 migrant sampling locations. Loci Fp54 and Fp92-1were both monomorphic for F. p. macropus. After adjusting for multiple comparisons, significant departures from Hardy-Weinberg equilibrium in the form of heterozygote deficiencies were observed in one locus (Fp92-1) among four sampling locations (Northwest and 1990 Greenland, F. p. tundrius; Northwest, F. p. anatum; Padre Island fall migrants 1985). Similarly, significant FIS values were observed with five sampled locations (Table 1; heterozygote deficit), and three (Greenland 1990 & 2001-04, F. p. tundrius; Northwest, F. p. anatum) of the five remained significant after excluding locus Fp921 from the analysis. No pairwise comparisons testing for linkage disequilibrium were significant after correcting for multiple comparisons. The majority of microsatellite genetic diversity estimates do not differ significantly (Wilcoxon signed-rank test, p.0.05) between geographic sampling locations in North America or between temporal sampling periods from Padre Island, Texas or Greenland (Table 1). Allelic richness (AR) varied from 3.760.5 (6 s.e.) alleles per locus in F. p. pealei to 4.660.7 alleles per locus in F. p. tundrius from western Greenland. The few cases for F. p. pealei possessed significantly lower AR compared to the tundrius populations in Nunavut (Z = 22.667; p = 0.008) and Greenland (Z = 22.667; p = 0.008) and the anatum population in Ontario (Z = 22.134; p = 0.033). Expected heterozygosity (He) ranged from 0.48560.092 in Ungava Bay F. p. tundrius from northeastern Canada to 0.55760.083 in F. p. anatum from eastern Canada. A significant difference in He was observed between F. p. tundrius populations in Ungava Bay and Nunavut (Z = 22.223; p = 0.026). Genetic diversity estimates for F. p. cassini and macropus from Argentina and Australia, respectively, were significantly lower in all comparisons of AR when compared to all sampling locations for the North American subspecies anatum, tundrius, and pealei (p,0.02); Table 1). F. p. macropus He was significantly lower than He estimates from all anatum, tundrius, and pealei populations (p,0.04), with the exception of Ungava Bay F. p. tundrius (Z = 21.778; p = 0.075). F. p. cassini He was not significantly different (p.0.06) from any of the He estimates from anatum, tundrius, and pealei populations. Population structure The posterior probability values for each value of K with STRUCTURE, while using the complete dataset that included F. p. cassini and macropus, plateau at K = 3 to K = 5 with Ln P(D) values from multiple runs at K = 3 (SD = 2.2) and K = 4 (SD = 3.6) being more consistent across runs compared to K = 5 (SD = 18.2; Fig. 2). When we used DK to infer the number of clusters, K = 3 was clearly inferred for the complete dataset using all sampled subspecies. Results from STRUCTURE when using only data from North American subspecies (F. p. pealei, tundrius, anatum) indicated the highest posterior probability values for K = 2 (Ln P(D) = 215843.5), while K = 1 and K = 3 had lower posterior probability values (215886.7 and 215937.6, respectively). We are unable to evaluate between K = 1 and K = 2 using the DK method (see Evanno et al. 2005). For K = 2 in this second analysis, F. p. tundrius and anatum sample locations possessed a relatively high proportion of membership to the same inferred cluster (0.870 to 0.959 and 0.759 to 0.922, respectively) and the majority of the pealei samples were assigned to the second cluster at a lower proportion (proportion of membership = 0.661). Samples collected from Padre Island across all years clustered with high support (.0.924; data not shown) with samples collected from breeding grounds identified as F. p. tundrius/anatum. The results from TESS corroborated those found with STRUCTURE. The minimum DIC value (30669) was achieved with Figure 2. Results from STRUCTURE analysis for all sampled peregrine falcon populations. (A) Assignment of individuals to K = 3 to 5 inferred clusters based on 11 microsatellite loci. Colors indicate different inferred clusters and their magnitude represents the posterior probability that the individual belongs to a particular cluster. (B) Estimated log probability values [Ln (P(D)] for each run for K = 1 to 8. The box highlighting K = 3 to 5 indicates the lowest values of K with the highest likelihood values. These results remained similar after excluding F. p. cassini and macropus from the analysis (data not shown). doi:10.1371/journal.pone.0014042.g002 K = 4 (K = 2, DIC value = 31736; K = 3, DIC value = 31163) and plateaus with higher values of K. The average value for the 20% best runs when K = 4 showed that the tundrius and anatum individuals clustered in a single group, which also included all of the migrant samples collected from Padre Island, whereas the remaining three subspecies, pealei, macropus and cassini, were in separate clusters each with high support (data not shown). Estimates of FST and Dest largely agreed with the results from STRUCTURE and TESS, indicating strong genetic differentiation between subspecies with the exception of those comparisons between F. p. tundrius and anatum which showed much lower levels of differentiation (Table S1). Although 16 of the 20 pairwise FST comparisons were significant (p,0.001) after correction for multiple comparisons between F. p. tundrius and anatum sampled breeding territories, the values were low (FST = 0.006 to 0.050) and similar in magnitude to those obtained from within subspecies comparisons (FST = 0.007 to 0.026 and 0.016 to 0.033 for tundrius and anatum, respectively) of which some were also significant (p,0.001; see Table S1). Estimates of Dest between F. p. tundrius and anatum were also low (Dest = 0.000 to 0.032), with only one comparison being significantly different from zero (Ontario & Greenland_1990). Dest values between F. p. tundrius and anatum were similar to those obtained from pairwise comparisons between sample locations within F. p. anatum (Dest = 0.003 to 0.029), while pairwise comparisons between tundrius subspecies locations were consistently low (Dest = 0.000 to 0.012) and not significantly different from zero (Table S1). After excluding samples from Ungava Bay due to low population sample size (n = 15), significant isolation by distance was observed among breeding sample locations of F. p. tundrius and anatum using both pairwise FST (r = 0.663; Mantel test P = 0.003) and Dest (r = 0.649; Mantel test P = 0.001) measures. Similarly, significant isolation by distance was observed among F. p. tundrius and anatum sample locations when southwest Ontario samples were excluded from the analyses (FST, r = 0.668, Mantel test P = 0.012; DEST, r = 0.579, Mantel test P = 0.035). However, isolation by distance was not supported when we included F. p. pealei genetic distance measures with tundrius and anatum breeding locations (FST, r = 0.092, Mantel test P = 0.324; Dest, r = 0.113, Mantel test P = 0.288). When comparing migrant peregrines from Padre Island with samples collected on breeding territories, far fewer pairwise FST and Dest comparisons were significant with those made with F. p. tundrius (23% and 0% out of 35 comparisons, respectively) than with anatum samples (71% and 14% out of 28 comparisons, respectively; Table S1). No pairwise FST or Dest comparisons between migrant Padre Island temporal samples were significant across sampling periods, indicating stable allele frequency distributions over a 22-year period. The simulated datasets of known size at Ne of #300 following seven generations were all significantly different from each other based on allele frequency distributions after sequential Bonferroni correction; whereas at Ne of 500, eight of forty-five comparisons were significant and at Ne of $1000, none of the pairwise comparisons were significant (data not shown). Effective population size Point estimates of Ne for the migrant peregrine falcon population varied depending on the choice of method, but in all cases, the values were high and ranged from 509 to .10,000 breeding individuals (Table 2). Reported 95% confidence levels around each point estimates were wide, with all of cases, regardless of method, extending to infinity, or at least the maximum allowable value (.10,000) used in each of the analyses (i.e., Ne-MAX). The choice of the census population size (N, 1,000 to 100,000; see Methods) used with the method implemented in the program TempoFs did not substantially influence our estimate of Ne. For example, Ne was 450 (117-infinity) and 516 (121-infinity) using an N of 1,000 or 100,000, respectively. Estimates of Ne from the method LDNe using spring migratory single time periods ranged from 187.8 in the 2001 dataset and .10,000 in the remaining three periods (1986, 1989, and 2007) with 95% CIs ranging between 59.5 to infinity. MLNE analyses on the Greenland population while allowing for immigration (see Methods) produced an NeOPEN estimate of 122.7 (95% CI 55.3-590.7) with a joint migration estimate of 0.103 (95% CI 0.024-0.234). The choice of source population did not substantially affect our estimates of NeOPEN. When we defined the source population for potential immigrants as Padre Spring 2007, our estimates of NeOPEN was 102.5 (95% CI 44.7-470.7) with a joint migration estimate of 0.221 (95% CI 0.056-0.734). Estimate of NeCLOSED for the Greenland peregrine falcon population was 158.9 (95% CI 75.7-804.5). Results from our simulated datasets of known size ranging from 50 to 5,000 breeding individuals further supported our empirical data suggesting that the migrant population is of large size. As population size increased across simulations, the accuracy and precision of each method for calculating Ne decreased (Fig. 3; Table 3). Above Ne of 500, for example, point estimates from all methods differed from the actual simulated size by more than 200 individuals with wide 95% confidence intervals (Fig. 3). Interestingly, point count estimates of Ne at levels below 500 were more often overestimated, while $500 tended to be underestimated. The one obvious exception was with LDNe where estimates of Ne$1,000 were overestimated (Fig. 3; Table 3). These results suggest that our estimate of the migrant peregrine falcon population Ne is at least 500 and possibly .1000 breeding individuals. Method used to estimate Ne MLNE 95% confidence intervals are provided in parentheses below each point estimate. Values indicated as infinity represent .10,000 breeding individuals (Ne-MAX of 10,000). 1Estimates of Ne from LDNe are based on a single time period (i.e., the 7th generation at the particular population size). doi:10.1371/journal.pone.0014042.t002 Figure 3. Estimates of Ne from simulated populations of known size. Three different temporal methods (TempoFs, TM3, MLNE) and a fourth method (LDNe) based on a single time period were used to estimate Ne. Similar to our empirical data from Padre Island, TX, temporal estimates were based on seven generations (T0 - T7), while the single time period estimate is from the simulated 7th generation (see methods). Bars represent deviated estimates of Ne from simulated population size and the dotted vertical lines reflect 95% confidence intervals. = values beyond the range of y-axis (see also Table 3 for point estimates used in this figure). doi:10.1371/journal.pone.0014042.g003 Table 3. Estimates of Ne from simulated populations of known size. Method used to estimate Ne 95% confidence intervals are provided in parentheses below each point estimate. Values indicated as infinity represent .10,000 breeding individuals (Ne-MAX of 10,000). 1Estimates of Ne from LDNe are based on a single time period (i.e., the 7th generation at the particular population size). doi:10.1371/journal.pone.0014042.t003 Population genetic data are a valuable tool for monitoring populations and species [e.g., 24,7], particularly those in the process of recovery. The peregrine falcon in North America is one such species that has required extensive monitoring to assess its progress toward achieving a sustainable population as defined by a federally mandated monitoring plan [76] subsequent to the species delisting from endangered status in 1999 (64 FR 4654146558). Here, we have utilized genetic data to assess the stability of the migrant peregrine falcon population in North America. Our results based on multiple methodological approaches indicate that this species is stable from a population genetics perspective as documented for breeding territories sampled throughout its entire northern latitudinal distribution including its migratory population at Padre Island, Texas. Assessment of population structure F. p. tundrius and anatum To support our use of peregrine falcon samples collected during migration at Padre Island to assess the species northern high latitude population genetic stability, it was important to document levels of genetic differentiation among all presumed northern breeding areas that contribute to the migrant population. In a previous study investigating the genetic structure of peregrine falcon sampling locations across Canada, Brown et al. [35] documented no differentiation based on FST and STRUCTURE analyses between F. p. tundrius and anatum using samples collected prior to their population decline in the 1950s. These results suggested that a continuous phenotypic cline existed based on the subtle morphological and behavioral characters used originally to describe the separation of the two subspecies (e.g., [77]). In contrast, however, Brown et al. [35] did document low but significant differentiation between the two subspecies when using contemporary samples. This was particularly the case with sample locations from southeastern Ontario where stronger allele frequency fluctuations and allelic introgression from captive released birds were more likely to have occurred as compared to their sampling locations further north and to the west (see also [78,79]; Fig. 2; Table S1). Nearly 7,000 peregrine falcons were released in the United States and Canada between the years 1974 and 1999 [13]. Many of these birds originated from breeding stock that included nonnative peregrine subspecies such as peregrinus and brookei from Europe, cassini from South America, and macropus from Australia, but to a lesser extent than individuals with either anatum, tundrius or pealei pure or mixed-ancestry [7981]. Because peregrine falcons were extirpated throughout southeastern Canada and the United States east of the Rocky Mountains, the peregrine falcons that now reside in these areas are largely the result of the release programs [e.g., 82]. Therefore, it is not surprising that contemporary low levels of genetic differentiation were observed between peregrine falcon populations sampled in geographic areas that have been reestablished with captive-bred individuals, while no genetic differentiation was observed prior to the decline. Other release and supplementation programs have also documented similar effects [e.g., 8386]. For example, Jacobsen et al. [29] reported significant microsatellite allele frequency changes before and after the reintroduction project of peregrine falcons (F. p. peregrinus) in southern Scandinavia [87]. Using F. p. tundrius and anatum samples from Brown et al. ([35]; n = 140) with additional sampling (n = 176) of both subspecies from previously unsampled geographic areas in Alaska (tundrius and anatum) and Greenland (tundrius), our results indicate little if no genetic population differentiation (STRUCTURE, TESS, FST and Dest) among sampling locations of these two subspecies. The significant values that were observed were those comparisons that included F. p. anatum samples from eastern Canada or geographically distant locations, including three other subspecies (Fig. 2; Table S1). A significant Mantel correlation between genetic (FST) and geographic distance matrices suggests an isolation by distance model of population differentiation (or regional equilibrium; [88]) throughout the high-latitude breeding distribution of tundrius and anatum. These results were also supported by a new measure (i.e., Dest, [55]) that is useful for investigating population differentiation between geographic locations when using highly polymorphic markers such as microsatellite DNA [see also 8991]. In this case, only a few pairwise tundrius/anatum or anatum/anatum comparisons were significant, all of which were those comparisons made with the Ontario sampled location (Table S1), and similar to FST, an isolation by distance model of population structure could not be rejected based on Dest. This is consistent with there being considerable gene flow among the populations considered as tundrius or anatum in Alaska, Canada and Greenland, and suggests that enough time may have passed for localized gene flow and genetic drift to stabilize and produce a pattern characteristic of migration-drift equilibrium throughout their high-latitude distribution (e.g., [88]). In contrast, all pairwise comparisons (FST and Dest) that included other subspecies pealei, cassini and macropus were significant. Although the degree of differentiation between pealei and tundrius/anatum was not as strong with STRUCTURE and TESS compared to FST and Dest, the lower proportion of membership for a few individuals from the two groups was likely due to pealei also being used in the captive breeding program for which some post release individuals subsequently possessed mixed ancestry [7981; see Figure 2]. Similarly, a few of the individuals collected as pealei may have been misidentified in the field where the subspecies distribution overlaps with anatum. Although, F. p. pealei are nonmigratory and they occupy coastal territories in the Pacific Northwest, extensive plumage variation in immature individuals exists, particularly in the southern portion of their range [13] where our pealei samples originated. No significant isolation by distance was observed after including the F. p. pealei sample location with F. p. tundrius/anatum breeding territories further suggesting no contemporary gene flow between pealei and tundrius/ anatum populations. Therefore, the STRUCTURE results were likely influenced by introgression of alleles from individuals used in the release program [e.g., 29] and/or the lack of power to detect population structure in cases with relatively low levels of differentiation (e.g., FST,0.03; [92,93]). With the exception of those comparisons made between high latitude populations of F. p. tundrius (see below), we remain cautious and resist equating reported population differentiation estimates with actual relative degree of differentiation between populations [e.g., 55,9496]. Peregrine falcons in North America, particularly those on the East Coast and Midwest in both Canada and the U.S., have experienced dramatic changes in population abundance over the past half-century [13,17,30]. These lower latitude populations are unlikely to have reached the equilibrium conditions necessary for a direct interpretation of measures of differentiation. Additional work is necessary to explore how these measures are influenced by data from populations with differing demographic histories, particularly those that experienced recent decline followed by a rapid recovery. Our results from STRUCTURE, TESS, FST and Dest do indicate, however, no genetic differentiation between F. p. tundrius and high latitude anatum sample locations in Alaska and northwest Canada. The high latitude sampling locations used in this study are all geographic regions that have not experienced any direct supplementation from translocated or captive bred individuals [97]. Similarly, the contemporary sampling locations from Alaska, northern Canada and Greenland do not differ in allele frequencies or levels of genetic diversity when compared to pre-decline historic samples from Canada ([35]; data not shown) suggesting genetic stability over time at least in their northern distribution. This conclusion is further supported by a significant isolation by distance relationship among breeding sample locations of these two subspecies, even after excluding samples from southern Ontario where introgression from captive-released birds is likely to have occurred (see [35]). A few studies have documented reduced levels of genetic variation in declining populations of birds of prey (e.g., [98100]). However, populations that have experienced recent declines yet have recovered in abundance can retain levels of genetic diversity when certain demographic conditions exist. For example, Hailer et al. [101] documented high levels of genetic diversity in recovering white-tailed eagle (Haliaeetus albicilla) populations suggesting that their long life span (,17 years) has helped buffer against the effects of genetic drift and loss of genetic variability [see also 102104]. In comparison, maximum longevity of banded wild peregrine falcons ranges from 16 to 20 years [13], suggesting that this factor may have helped populations maintain high levels of genetic variability in their northern distribution in North America, despite significant declines in abundance. After adjusting significance levels to account for multiple simultaneous comparisons, five sample locations possessed significant FIS values, or heterozygote deficiencies (Table 1). Possible explanations for these results include null alleles, sampled multiple populations (Wahlund effect), or nonrandom mating within sampled locations. Two of the populations were no longer significant after excluding locus Fp92-1 from the analysis, suggesting possible null alleles with this locus. For the three remaining sampled locations, despite observing little population genetic differentiation (FST, DEST, STRUCTURE, TESS) between breeding areas across their North American high latitude distribution, the more plausible explanation for significant FIS values is due to sampling multiple subpopulations within the three areas (1990 & 2001-04 Greenland F. p. tundrius and Northwest F. p. anatum). For example, the samples for Greenland 2001-04 were obtained from both Thule (n = 15) and Kangerlussuaq (n = 27), two geographic areas separated by .1,000 km, and the samples from Northwest F. p. anatum were obtained from a large geographic area throughout Alaska and northwest Canada [see Figure 1]. Although fewer samples were obtained from Thule, significant FIS values were observed, but not with Kangerlussuaq after separating the two datasets (data not shown). More work is required to investigate whether fine scale geographic structure may exist in areas on the periphery of the species distribution and whether recent expansion or growth may influence these results (e.g., [105]) because peregrine falcons have recently expanded northward into areas such as Thule [33]. We currently do not possess an adequate number of samples collected in the United States to determine if peregrine falcons (i.e., anatum) show a similar lack of population breeding structure south of the Canadian and U.S. border (e.g., continental U.S.) and whether contemporary gene flow exists throughout their latitudinal distribution. We feel that this deficiency in sampling, however, does not negate the utility of the current analysis because we were primarily interested in determining overall genetic stability and effective population size of migrant passage peregrine falcon population. It has been shown that peregrine falcons that breed south of the U.S./Canada border possess reduced migratory behavior than those further north [11,13,32] and, therefore, less likely sampled in Padre Island, TX. Although, anecdotal evidence suggests that individuals with mixed tundrius ancestry that were released in eastern U.S. may migrate further south than predecline individuals (i.e., anatum) from the same area [81], other work in the Midwestern U.S. has documented an increasing number of urban-nesting peregrine falcons overwintering consistently in or near their breeding territories [106]. The most recent USFWS Monitoring Results for F. p. anatum in 2003 [107] reported that while 92% of recorded nest sites throughout five of the six defined regions for monitoring purposes in continental U.S. were located on natural substrates (e.g., cliffs), 68% in the Midwestern/Northeast region were on human-built structures in urban settings such as tall buildings and bridges [see also 79,106,108,109]. Whether introgression of non-native genes into the breeding population in the U.S. has had any negative effects on fitness (e.g., outbreeding depression; [110]) or changes in population dispersal patterns remains to be shown ([79]; see also [111]). Additional work is required specifically to address the genetic stability of populations in the contiguous U.S., whereas our study is primarily focused on high latitude populations, which are more likely the source of migrant birds passing through Padre Island. Migrant population genetic stability Peregrine falcons sampled during migration at Padre Island, Texas clustered with high support with individuals sampled throughout their northern breeding distribution in Alaska, Canada and Greenland (F. p. tundrius/anatum; Fig. 2). In contrast to peregrine falcons sampled in southern Canada (anatum), the Padre Island samples had consistently higher posterior probability assignment values similar to peregrine falcons identified as tundrius/anatum and possessed little if no signal indicating admixture from the other sampled subspecies (Fig. 2). These results, along with results from FST and DEST pairwise comparison (Table S1), therefore suggest that peregrine falcons passing through Padre Island were likely individuals with breeding territories located further north than southeastern Canada. Previous studies investigating migratory patterns of peregrine falcons in North America using banding records or satellite telemetry have also documented peregrines passing through Padre Island that originate or finalize their migration in northern high latitude areas rather than further south in southern Canada and continental U.S. ([11,32,33]; see also Fig. 2 in [31]). Fuller et al. [12], for example, identified a wide distribution of breeding territories across northern latitudes for passage peregrines (n = 54), including those surveyed migrating through Padre Island. Across migratory seasons (fall and spring) and seven time periods sampled over a twenty-two year period, no significant changes in levels of microsatellite diversity were observed at Padre Island, and diversity levels were similar to those obtained from high latitude breeding peregrine falcon sample locations (Table 1). These results suggest that the population is large enough in size to offset the negative effects of drift, with adequate levels of gene flow between areas (see also [112115]). Further, no significant levels of differentiation were observed between each of the sampled time periods (Fig. 2; Table S1). The latter result is important because monitoring changes in population differentiation (e.g., FST or DEST) is often a more sensitive indicator of population decline than is the loss of allelic diversity [116,117]. Similarly, no significant change in diversity levels (Table 1) or population differentiation (Table S1) was observed between the two temporal sample periods (1990 and 2001-04) in western Greenland. Multiple studies have documented significant allele frequency change associated with increased population differentiation in small or declining populations [74,118121]. With our simulated datasets of known size, for example, the development of significant population differentiation (FST) was observed in all pairwise comparisons at Ne of #300 in as little as seven generations (the time period between our samples collected from Padre). However, at Ne of 500 only eight of the forty-five comparisons were significant (18%) and none of the comparisons at Ne of $1000 were significant (0%). These results, along with those from the USFWS nationwide monitoring efforts [107], suggest that the higher latitude migratory and breeding peregrine falcon population is stable with no indication of decline. In fact, monitoring efforts in the field, including multiple long-term migration watchsites (e.g., Cape May Bird Observatory, Hawk Mountain Sanctuary, Hawk Ridge Bird Observatory) suggest that this species continues to increase in abundance [13,31,122124]. Similarly, levels of organochlorine pesticides continue to decline in peregrine falcon migrants returning from Central and South America. Henny et al. [125] reported a 9697% decline in blood DDE concentrations in female peregrine falcons sampled between 1994 (n = 45) and 2004 (n = 27) at Padre Island. Out of the 27 adult peregrine falcons sampled in 2004, DDE concentrations were below detectible levels (,0.02 mg/g) in 20 birds (77%), while in contrast only two of the 156 adult samples (1%) between 1978 and 1994 were below the detectible limit [125]. These are definitely reassuring signs that the peregrine falcon population is moving toward full recovery. Our inability to obtain a precise estimate of Ne for the migratory population of peregrine falcons also suggests a large population [see also 126,127]. Estimates of Ne ranged from 509 to infinity (.10,000; Ne-MAX), with extremely wide 95% confidence intervals (Table 2). The power to estimate Ne using genetic data is dependent on multiple factors. When populations are of small size (,500 breeding individuals), a variety of methods, some of which were employed in this study, do exceptionally well in inferring how strong genetic drift was or how large the Ne of the population must be to cause the observed change in allele frequencies over time when assuming no mutation, selection and migration during the sampled time period [58,64,127]. When populations are of large size, however, allele frequencies are less likely to change due to drift and our ability to estimate Ne becomes much more difficult. As was observed with our estimates of population differentiation measures (FST and DEST; Table S1) between temporal Padre Island sampling periods, no significant changes in allele frequencies were identified over the 22-year period (,7 generations), further supporting that the breeding population of high latitude peregrine falcons is of large size. This conclusion is also supported by our estimates of Ne from the local Greenland peregrine falcon population (Ne = ,120) while accounting for immigration, suggesting that a much larger migratory population must exist when extrapolated to the remainder of its breeding distribution across North America. There are ways to improve the precision of our estimate of Ne. Either we increase the sample size of each time period and/or increase the number of loci characterizing allele frequency change over time. In a recent study investigating Ne of an Australian tiger prawn (Penaeus esculentus) population, Ovenden et al. [128] determined using simulation that they would require 2,000 samples taken one generation apart to reliably estimate Ne of about 8,000 breeding individuals with eight microsatellite loci using similar temporal methods as this study. By increasing the time between sampling periods to four generations, their simulations suggested that the same sample size would produce accurate estimates for Ne of 10,000; however, by decreasing the sample size to 1,000 individuals, their ability to obtain finite estimates dropped from 100% to 65% [128]. In this study, our samples sizes were 46 individuals for the majority of temporal periods. We could possibly double the sample size to approximately 100 individuals per time period, but sampling beyond that number is unrealistic given the difficulty in trapping migratory falcons. Palstra & Ruzzante [64] recommend that the S/Ne ratio (where S is the sample size) be approximately 0.10 for adequate sampling for the temporal approach for estimating Ne with genetic data. With a sample size of 46 individuals, we should be able to provide reliable estimates of Ne,460. Therefore, the wide 95% confidence intervals, or low precision, we obtained from all of the temporal methods for estimating Ne in this study suggests that the actual Ne is of larger size. Alternatively, increasing the number of loci would also increase the power to estimate Ne from large populations using genetic data; however, the extent of its improvement will depend on the variability of the loci with increased polymorphism required [64]. In a recent review, Leberg [1] commented that increasing the number of loci sampled will increase the precision of estimates to a greater extent than increasing the number of individuals sampled. (but see [127]). Therefore, because of difficulties obtaining large sample sizes, it may prove worthwhile to explore additional markers, such as SNPs [e.g., 129131] and additional microsatellite loci [e.g., 132,133] to obtain a more precise estimate of Ne that can be used to make future management decisions for this species. Conservation implications The above results can inform future management decisions impacting the full recovery of peregrine falcons breeding in North America. Recently, the U.S. Fish and Wildlife Service allowed the take in 2009 of up to 36 first-year (FY) autumn migrant passage peregrine falcons east of 100uW longitude for use in falconry [31,134]. The plan made a distinction between peregrines with natal sites south of 54uN latitude (F. p. anatum) and those further north, which includes both F. p. tundrius and northern F. p. anatum subspecies. Peregrine falcons (F. p. anatum) south of 54uN latitude possess reduced migratory behavior compared to those further north, producing a leap-frog breeding/migratory peregrine distribution [11,13,31,32]. This distinction between northern migratory birds and those further south was important because the two geographic groups differ in their current census estimates. Approximately 2,700 to 8,000 pairs, which include both tundrius and anatum subspecies, are estimated for the northern population ($54uN latitude) based on non-genetic methods. The southern populations (,54 uN latitude) east of 100uW longitude (F. p. anatum) is estimated at ,450 pairs, while west of 100uW longitude (F. p. anatum and F. p. pealei) consist of ,1,400 to 1,800 pairs [31]. Therefore, the USFWS specifically targets the take of migratory individuals that breed $54uN latitude. Based on our genetic results, including census estimates from the field, and because the proposal specified autumn first-year individuals, the removal of 36 FY autumn migrants from the population is unlikely to adversely affect their continued recovery [see also 31]. First-year survivorship in the wild is estimated at 40 50%, while breeding adult survival for migrants is likely between 8085% [13]. If the peregrine falcon population was of small size (Ne,500), the methods that we employed should have had sufficient power to provide precise estimates of Ne; yet, our inability to obtain such an estimate suggests this population is of larger size. The estimate of Ne as measured in this study (with the exception of Ne from LDNe; see [1,73]) reflects the harmonic mean effective size of the migratory high latitude peregrine falcon population over the sampled seven-generation time period (e.g., 19892007). It roughly approximates to the number of breeding individuals that produce offspring that live to reproductive age during the sampled time period, which is typically a value much smaller than the actual census size (N) of the population (Ne/N,0.11; see [135]). Genetic measures of Ne incorporate all demographic effects in their estimate, such as fluctuating population size, unequal sex ratio, and variance in reproductive success [1,136138], which all decrease Ne relative to N. Therefore, even if we assume that Ne of 500 is correct, the actual census size is likely an order of magnitude larger; however, more work is required to determine the actual ratio between the two measures specific to peregrine falcon populations. The number of non-breeders (floaters) can outnumber the actual breeders in some areas by severalfold [13,139], and the numbers of migrants recorded at specific monitoring sites in the U.S. are high. For example, between 1999 and 2004, the mean annual fall migration peregrine falcon count at Curry Hammock Florida State Park in the middle Keys alone was 1,908 individuals (1,4322,858, minmax; [122]). Taking into account the relatively high first year mortality for migrant peregrine falcons ($5060%, [13]), the Ne/N ratio is further reduced due to a large number of individuals not producing offspring relative to fall census estimates. Although we were unable to provide a precise point estimate of Ne in this study, we can conclude that the Ne for the migrant peregrine falcon population is unlikely to be smaller than 500. This agrees with field data suggesting a much larger breeding population size [13,31]. Reducing organochlorine pesticides and other contaminants (e.g., mercury; [140,141]) in their environment is of greater importance for securing the long-term viability of peregrine falcon populations, and recent results from Henny et al. [125] suggest that these conditions are improving. Consequences of the illegal take of peregrine falcons in their wintering distribution in South America [142] also deserves more attention, and certainly, local monitoring of specific areas within the species breeding distribution should continue. With the exception of our results from Greenland, the analyses using migrant samples from Padre Island provide a coarse description of population genetic stability over time for the high-latitude breeding population of this species, while more local demographic perturbations associated with specific breeding locations (see [143]) would not necessarily be reflected in our results using migrant samples alone. Supporting Information Table S1 Pairwise estimates of FST (below diagonal) and Dest (above diagonal) based on 11 microsatellite loci between regional peregrine falcon sample locations. Found at: doi:10.1371/journal.pone.0014042.s001 (0.09 MB PDF) We thank the managers and staff of Laguna Atascosa National Wildlife Refuge, managers and staff of USGS Bird Banding Laboratory, Dr. William G. Mattox (Greenland Peregrine Falcon Survey), Dr. Mark R. Fuller, and all the field researchers who captured peregrines and collected, processed and archived samples. Marit Nesje and Jan Lifjeld graciously provided a portion of the data and samples from Australia, and samples from Alaska and Canada were collected by numerous field biologists and provided by Jon Longmire. We thank The Danish Polar Center and Greenland Home Rule Government for providing permits and allowing us to work in Greenland. Mention of trade names or commercial products does not constitute endorsement or recommendation for use. Conceived and designed the experiments: JAJ. Performed the experiments: JAJ ST GKS JWB. Analyzed the data: JAJ. Contributed reagents/ materials/analysis tools: JAJ ST GKS KKB JWB TLM WSS MAY BA DPM. Wrote the paper: JAJ ST KKB JWB DPM. 1. Leberg P ( 2005 ) Genetic approaches for estimating the effective size of populations . Journal of Wildlife Management 69 : 1385 - 1399 . 2. Schwartz MK , Luikart G , Waples RS ( 2007 ) Genetic monitoring as a promising tool for conservation and management . Trends in Ecology and Evolution 22 : 25 - 33 . 3. Laikre L , Nilsson T , Primmer CR , Ryman N , Allendorf FW ( 2009 ) Importance of genetics in the interpretation of favourable conservation status . Conservation Biology 23 : 1378 - 1381 . 4. De Barba M , Waits LP , Garton EO , Genovesi P , Randi E , et al. ( 2010 ) The power of genetic monitoring for studying demography, ecology and genetics of a reintroduced brown bear population . Molecular Ecology 19 : 3938 - 3951 . 5. Hall D ( 2005 ) Memo: Policy on Genetics in Endangered Species Activity , January 27 , 2005 . US Fish & Wildlife Service. Available: fws/05_24_5_hallmemo.pdf via the Internet. Accessed 30 Apr 2010 . 6. Morgenweck R ( 2005 ) Letter from Ralph Morgenweck to Dale Hall , March 11 , 2005 . Available: pdf via the Internet. Accessed 30 Apr 2010 . 7. Laikre L , Allendorf FW , Aroner LC , Baker CS , Gregovich DP , et al. ( 2010 ) Neglect of genetic diversity in implemetation of the Convention on Biological Diversity . Conservation Biology 24 : 86 - 88 . 8. Ratcliffe D ( 1993 ) The Peregrine Falcon , 2nd edn. London: T & AD Poyser. 9. Ferguson-Lees J , Christie DA ( 2001 ) Raptors of the World . London: Christopher Helm. 10. White CM , Boyce DA , Jr. ( 1988 ) An overview of peregrine falcon subspecies . In: Cade TJ, Enderson JH , Thelander CG , White CM, eds. Peregrine Falcon Populations: Their Management and Recovery . Boise, Idaho: The Peregrine Fund . pp 789 - 810 . 11. Schmutz JK , Fyfe RW , Banasch U , Armbruster H ( 1991 ) Routes and timing of migration of falcons banded in Canada . Wilson Bulletin 103 : 44 - 58 . 12. Fuller MR , Seegar WS , Schueck LS ( 1998 ) Routes and travel rates of migrating Peregrine Falcons Falco peregrinus and Swainson's Hawks Buteo swainsoni in the Western Hemisphere . Journal of Avian Biology 29 : 433 - 440 . 13. White CM , Clum NJ , Cade TJ , Hunt WG ( 2002 ) Peregrine Falcon (Falco peregrinus) . In: Poole A, Gill F, eds. The Birds of North America, No. 660. Philadelphia: The Birds of North America, Inc. 14. Falk K , Mller S ( 1988 ) Status of the peregrine falcon in South Greenland: population density and reproduction . In: Cade TJ, Enderson JH , Thelander CG , White CM, eds. Peregrine Falcon Populations: Their Management and Recovery . Boise, Idaho: The Peregrine Fund . pp 37 - 43 . 15. Cade TJ ( 1982 ) The Falcons of the World . Ithaca, NY: Cornell University Press. 16. Fyfe RW ( 1969 ) The Peregrine Falcon in northern Canada . In: Hickey JJ, ed. Peregrine Falcon Populations: Their Biology and Decline . Madison, WI: University of Wisconsin Press. pp 101 - 114 . 17. Kiff LF ( 1988 ) Commentary - Changes in the status of the Peregrine in North America: an overview . In: Cade TJ, Enderson JH , Thelander CG , White CM, eds. Peregrine Falcon Populations: Their Management and Recovery . Boise, Idaho: The Peregrine Fund . pp 123 - 139 . 18. Jefferies DJ , Prestt I ( 1966 ) Post-mortems of peregrines and lanners with particular reference to organochlorine residues . British Birds 59 : 49 - 64 . 19. Ratcliffe D ( 1967 ) Decreases in eggshell weight in certain birds of prey . Nature 215 : 208 - 210 . 20. Cade TJ , Lincer JL , White CM , Roseneau DG , Swartz LG ( 1971 ) DDE residues and eggshell changes in Alaska falcons and hawkbills. Science 172 : 955 - 957 . 21. Peakall DB , Reynolds LM , French MC ( 1976 ) DDE in eggs of the peregrine falcon . Bird Study 23 : 183 - 186 . 22. Berger DD , Sidebar CR , Jr., Gamble KE ( 1969 ) The status of breeding Peregrines in the eastern United States . In: Hickey JJ, ed. Peregrine Falcon Populations: Their Biology and Decline . Madison, WI: University of Wisconsin Press. pp 165 - 173 . 23. Berry RB ( 2003 ) History and Extinction of the Appalachian Peregrine . In: Cade TJ, Burnham W, eds. Return of the Peregrine . Boise, ID: The Peregrine Fund . pp 35 - 55 . 24. Enderson JH ( 1969 ) Population trends among Peregrine Falcons in the Rocky Mountain region . In: Hickey JJ, ed. Peregrine Falcon Populations: Their Biology and Decline . Madison, WI: University of Wisconsin Press. pp 73 - 79 . 25. Burnham W ( 2003 ) Peregrine Falcon Restoration in the Rocky Mountains/ Northwest . In: Cade TJ, Burnham W, eds. Return of the Peregrine . Boise, ID: The Peregrine Fund . pp 127 - 153 . 26. Fyfe RW , Temple SA , Cade TJ ( 1976 ) The 1975 North American peregrine falcon survey . Canadian Field Naturalist 90 : 228 - 273 . 27. White CM , Emison WB , Williamson FL ( 1973 ) DDE in a resident Aleutian Island peregrine population . The Condor 75 : 306 - 311 . 28. Nelson RW ( 1990 ) Status of the Peregrine Falcon, Falco peregrinus pealei , on Langara Island, Queen Charlotte Islands, British Columbia , 1968 - 1989 . The Canadian Field-Naturalist 104 : 193 - 199 . 29. Jacobsen F , Nesje M , Bachmann L , Lifjeld JT ( 2008 ) Significant genetic admixture after reintroduction of peregrine falcon (Falco peregrinus) in southern Scandinavia . Conservation Genetics 9 : 581 - 591 . 30. Cade TJ , Burnham W ( 2003 ) Return of the Peregrine. A North American Saga of Tenacity and Teamwork . Boise, Idaho: The Peregrine Fund. 31. United States Fish and Wildlife Service ( 2008 ) Final environmental assessment and management plan: take of migrant peregrine falcons from the wild for use in falconry, and reallocation of nestling/fledgling take . Arlington, Virginia: U.S. Fish and Wildlife Service, Division of Migratory Bird Management. 32. McGrady MJ , Maechtle TL , Vargas JJ , Seegar WS , Porras Pen a MC ( 2002 ) Migration and ranging of peregrine falcons wintering on the Gulf of Mexico coast , Tamaulipas, Mexico. The Condor 104 : 39 - 48 . 33. Burnham KK , Burnham WA , Newton I , Johnson JA , Gosler AG ( 2010 ) Recent population expansion of Peregrine Falcons in northwest Greenland as a result of climate change . Meddelelser om Grnland, BioScience . In press. 34. Johnson JA , Burnham KK , Burnham WA , Mindell DP ( 2007 ) Genetic structure among continental and island populations of gyrfalcons . Molecular Ecology 16 : 3145 - 3160 . 35. Brown JW , Van Coeverden de Groot PJ , Birt TP , Seutin G , Boag PT , et al. ( 2007 ) Appraisal of the consequences of the DDT-induced bottleneck on the level and geographic distribution of neutral genetic variation in Canadian peregrine falcons, Falco peregrinus . Molecular Ecology 16 : 327 - 343 . 36. Nesje M , Red KH , Lifjeld JT , Lindberg P , Steen OF ( 2000 ) Genetic relationships in the peregrine falcon (Falco peregrinus) analysed by microsatellite DNA markers . Molecular Ecology 9 : 53 - 60 . 37. Talbot SL , Palmer AG , Sage GK , Sonthagen SA , Swem T , et al. ( 2010 ) Lack of genetic polymorphism among peregrine falcons of Fiji. Avian Biology . In review. 38. Lewis PO , Zaykin D ( 2001 ) Genetic Data Analysis: Computer Program for the Analysis of Allelic Data , Version 1.1. Storrs , Connecticut: University of Connecticut. 39. Rice WR ( 1989 ) Analyzing tables of statistical tests . Evolution 43 : 223 - 225 . 40. Goudet J ( 1995 ) FSTAT (version 1.2): a computer program to calculate Fstatistics . Journal of Heredity 86 : 485 - 486 . 41. Leberg P ( 2002 ) Estimating allelic richness: effects of sample size and bottlenecks . Molecular Ecology 11 : 2445 - 2449 . 42. Raymond M , Rousset F ( 1995 ) GENEPOP (version 1.2): population genetics software for exact tests and ecumenicism . Journal of Heredity 86 : 248 - 249 . 43. Rousset F ( 2008 ) Genepop'007: a complete reimplementation of the Genepop software for Windows and Linux . Molecular Ecology Resources 8 : 103 - 106 . 44. Pritchard J , Stephens M , Donnelly P ( 2000 ) Inference of population structure using multilocus genotype data . Genetics 155 : 945 - 959 . 45. Falush D , Stephens M , Pritchard J ( 2003 ) Inference of population structure using multilocus genotype data: linked loci and correlated allele frequencies . Genetics 164 : 1567 - 1587 . 46. Evanno G , Regnaut S , Goudet J ( 2005 ) Detecting the number of clusters of individuals using the software Structure: a simulation study . Molecular Ecology 14 : 2611 - 2620 . 47. Rosenberg N ( 2004 ) DISTRUCT: a program for the graphical display of population structure . Molecular Ecology Notes 4 : 137 - 138 . 48. Chen C , Durand E , Forbes F , Francois O ( 2007 ) Bayesian clustering algorithms ascertaining special population structure: a new computer program and a comparison study . Molecular Ecology Notes 7 : 747 - 756 . 49. Francois O , Ancelet S , Guillot G ( 2006 ) Bayesian Clustering Using Hidden Markov Random Fields in Spatial Population Genetics . Genetics 174 : 805 - 816 . 50. Durand E , Jay F , Gaggiotti OE , Francois O ( 2009 ) Spatial inference of admixture proportions and secondary contact zones . Molecular Biology and Evolution 26 : 1963 - 1973 . 51. Francois O , Blum MGB , Jakobsson M , Rosenberg NA ( 2008 ) Demographic history of European populations of Arabidopsis thaliana . PLoS Genetics 4 ( 5 ): e1000075 . 52. Jakobsson M , Rosenberg NA ( 2007 ) CLUMPP: a cluster matching and permutation program for dealing with label switching and multimodality in analysis of population structure . Bioinformatics 23 : 1801 - 1806 . 53. Weir BS , Cockerham CC ( 1984 ) Estimating F-statistics for the analysis of population structure . Evolution 38 : 1358 - 1370 . 54. Excoffier L , Laval G , Schneider S ( 2005 ) ARLEQUIN ver . 3 .0: An integrated software package for population genetics data analysis . Evolutionary Bioinformatics Online 1 : 47 - 50 . 55. Jost L ( 2008 ) GST and its relatives do not measure differentiation . Molecular Ecology 17 : 4015 - 4026 . 56. Chao A , Shen T-J SPADE . Available: softwareCE.html via the Internet. Accessed 25 Dec 2009 . 57. Jensen JL , Bohonak AJ , Kelly ST Isolation by distance, web service . BMC Genetics 6 : 13 . [v.3.16]. 58. Wang J ( 2005 ) Estimation of effective population sizes from data on genetic markers . Philosophical Transactions of the Royal Society Series B 360 : 1395 - 1409 . 59. Fraser DJ , Hansen MM , stergaard S , Tessier N , Legault M , et al. ( 2007 ) Comparative estimation of effective population size and temporal gene flow in two contrasting population systems . Molecular Ecology 16 : 3866 - 3889 . 60. Waples RS ( 2006 ) A bias correction for estimates of effective population size based on linkage disequilibrium at unlinked gene loci . Conservation Genetics 7 : 167 - 184 . 61. Tallmon DA , Koyuk A , Luikart G , Beaumont MA ( 2008 ) ONeSAMP: a program to estimate effective population size using approximate Bayesian computation . Molecular Ecology Resources 8 : 299 - 301 . 62. Hill WG ( 1981 ) Estimation of effective population size from data on linkage disequilibrium . Genetical Research 38 : 209 - 216 . 63. Waples RS , Do C ( 2008 ) LDNE: a program for estimating effective population size from data on linkage disequilibrium . Molecular Ecology Resources 8 : 753 - 756 . 64. Palstra FP , Ruzzante DE ( 2008 ) Genetic estimates of contemporary effective population size: what can they tell us about the importance of genetic stochasticity for wild population persistence? Molecular Ecology 17 : 3428 - 3447 . 65. Wang J , Whitlock MC ( 2003 ) Estimating effective population size and migration rates from genetic samples over space and time . Genetics 163 : 429 - 446 . 66. Krimbas CB , Tsakas S ( 1971 ) The genetics of Dacus oleae V. Changes of esterase polymorphism in a natural population following insecticide control: selection or drift? Evolution 25 : 454 - 460 . 67. Pollak E ( 1983 ) A new method for estimating the effective population size from allele frequency changes . Genetics 104 : 531 - 548 . 68. Waples RS ( 1989 ) A generalized approach for estimating effective population size from temporal changes in allele frequency . Genetics 121 : 379 - 391 . 69. Waples RS ( 1990 ) Conservation genetics of Pacific salmon: III, estimating effective population size . Journal of Heredity 81 : 277 - 289 . 70. Turner TF , Salter LA , Gold JR ( 2001 ) Temporal-method estimates of Ne from highly polymorphic loci . Conservation Genetics 2 : 297 - 308 . 71. Jorde PE , Ryman N ( 2007 ) Unbiased estimator for genetic drift and effective population size . Genetics 177 : 927 - 935 . 72. Berthier P , Beaumont MA , Cornuet J-M , Luikart G ( 2002 ) Likelihood estimation of the effective population size using temporal changes in allele frequencies: a genealogical approach . Genetics 160 : 741 - 751 . 73. Wang J ( 2001 ) A pseudo-likelihood method for estimating effective population size from temporally placed samples . Genetical Research 78 : 243 - 257 . 74. Johnson JA , Bellinger MA , Toepfer JE , Dunn PO ( 2004 ) Temporal changes in allele frequencies and low effective population size in greater prairie-chickens . Molecular Ecology 13 : 2617 - 2630 . 75. Kuo CH , Janzen FJ ( 2003 ) BottleSim: a bottleneck simulation program for long-lived species with overlapping generations . Molecular Ecology Notes 3 : 669 - 673 . 76. United States Fish and Wildlife Service ( 2003 ) Monitoring plan for the American Peregrine Falcon, a species recovered under the Endangered Species Act . PortlandOregon: U.S. Fish and Wildlife Service, Divisions of Endangered of Endangered Species and Migratory Birds and State Programs, Pacific Region. 77. White CM ( 1968 ) Diagnosis and relationships of the North American tundrainhabiting peregrine falcons . Auk 85 : 179 - 191 . 78. Nesje M , Red KH , Bell DA , Lindberg P , Lifjeld JT ( 2000 ) Microsatellite analysis of population structure and genetic variability in peregrine falcons (Falco peregrinus) . Animal Conservation 3 : 267 - 275 . 79. Tordoff HB , Redig PT ( 2001 ) Role of genetic background in success of reintroduced Peregrine Falcons . Conservation Biology 15 : 528 - 532 . 80. Temple SA ( 1988 ) Future goals and needs for the management and conservation of the Peregrine Falcon . In: Cade TJ, Enderson JH , Thelander CG , White CM, eds. Peregrine Falcon Populations: Their Management and Recovery . BoiseIdaho : The Peregrine Fund . pp 843 - 848 . 81. Cade TJ , Tordoff HB , Barclay JH ( 2000 ) Re-introduction of peregrines in the eastern United States: an evaluation . Re-Introduction News 19 : 19 - 21 . 82. Tordoff HB , Redig PT ( 1997 ) Midwest peregrine falcon demography , 1982 - 1995 . Journal of Raptor Research 31 : 339 - 346 . 83. DeYoung RW , Demarias S , Honeycutt RL , Rooney AP , Gonzales RA , et al. ( 2003 ) Genetic consequences of white-tailed deer (Odocoileus virginianus) restoration in Mississippi . Molecular Ecology 12 : 3237 - 3252 . 84. Vernesi C , Crestanello B , Pecchioli E , Tartari D , Caramelli D , et al. ( 2003 ) The genetic impact of demographic decline and reintroduction in the wild boar (Sus scrofa): a microsatellite analysis . Molecular Ecology 12 : 585 - 595 . 85. Eldridge WH , Naish KA ( 2007 ) Long term effects of translocation and release numbers on fine scale population structure among coho salmon (Onchorhynchus kisutch) . Molecular Ecology 16 : 2407 - 2421 . 86. Fitzpatrick BM , Johnson JR , Kump DK , Shaffer HB , Smith JJ , et al. ( 2009 ) Rapid fixation of non-native alleles revealed by genome-wide SNP analysis of hybrid tiger salamanders . BMC Evolutionary Biology 9 : 176 . 87. Lindberg P ( 1988 ) Reintroducing the peregrine falcon in Sweden . In: Cade TJ, Enderson JH , Thelander CG , White CM, eds. Peregrine Falcon Populations: Their Management and Recovery . Boise, Idaho: The Peregrine Fund . pp 619 - 628 . 88. Hutchison DW , Templeton AR ( 1999 ) Correlation of pairwise genetic and geographic distance measures: inferring the relative influences of gene flow and drift on the distribution of genetic variability . Evolution 53 : 1898 - 1914 . 89. Jost L ( 2009 ) D vs . GST: Response to Heller and Siegismund (2009) and Ryman and Leimar ( 2009 ). Molecular Ecology 18 : 2088 - 2091 . 90. Heller R , Siegismund H ( 2009 ) Relationship between three measures of genetic differentiation GST,DEST, and GST': how wrong have we been? Molecular Ecology 18 : 2080 - 2083 . 91. Gerlach G , Jueterbock A , Kraemer P , Deppermann J , Harmand P ( 2010 ) Calculations of population differentiation based on GST and D: forget GST but not all of statistics! Molecular Ecology 19 : 3845 - 3852 . 92. Latch EK , Dharmarajan G , Glaubitz JC , Rhodes OE , Jr. ( 2006 ) Rela- tive performance of Bayesian clustering software for inferring population substructure and individual assignment at low levels of population differentiation . Conservation Genetics 7 : 295 - 302 . 93. Waples RS , Gaggiotti O ( 2006 ) What is a population? An empirical evaluation of some genetic methods for identifying the number of gene pools and their degree of connectivity . Molecular Ecology 15 : 1419 - 1439 . 94. Balloux F , Bru nner H , Lugon-Moulin N , Hausser J , Goudet J ( 2000 ) Microsatellites can be misleading: an empirical and simulation study . Evolution 54 : 1414 - 1422 . 95. Balloux F , Lugon-Moulin N ( 2002 ) The estimation of population differentiation with microsatellite markers . Molecular Ecology 11 : 155 - 165 . 96. Ryman N , Leimar O ( 2009 ) GST is still a useful measure of genetic differentiation - a comment on Jost's D. Molecular Ecology 18 : 2084 - 2087 . 97. Ambrose RE , Ritchie RJ ( 2003 ) Peregrine falcons , naturalists, scientists, and biologists in Alaska, 1900 to 2000 . In: Cade TJ, Burnham W, eds. Return of the Peregrine . Boise, ID: The Peregrine Fund . pp 189 - 197 . 98. Groombridge JJ , Jones CG , Bruford MW , Nichols RA ( 2000 ) ' Ghost' alleles of the Mauritius kestrel . Nature 403 : 616 . 99. Ewing SR , Nager RG , Nicoll MA , Aumjaud A , Jones CG , et al. ( 2008 ) Inbreeding and loss of genetic variation in a reintroduced population of Mauritius kestrel . Conservation Biology 22 : 395 - 404 . 100. Gautschi B , Muller JP , Schmid B , Shykoff JA ( 2003 ) Effective number of breeders and maintenance of genetic diversity in the captive bearded vulture population . Heredity 91 : 9 - 16 . 101. Hailer F , Helander B , Folkestad AO , Ganusevich SA , Garstad S , et al. ( 2006 ) Bottlenecked but long-lived: high genetic diversity retained in white-tailed eagles upon recovery from population decline . Biology Letters 2 : 316 - 319 . 102. Halbert ND , Derr JN ( 2008 ) Patterns of genetic variation in US federal bison herds . Molecular Ecology 17 : 4963 - 4977 . 103. Johnson JA , Gilbert M , Virani MZ , Asim M , Mindell DP ( 2008 ) Temporal genetic analysis of the critically endangered oriental white-backed vulture in Pakistan . Biological Conservation 141 : 2403 - 2409 . 104. Lerner HRL , Johnson JA , Lindsay AR , Kiff LF , Mindell DP ( 2009 ) It's not too late for the harpy eagle (Harpia harpyja): high levels of genetic diversity and differentiation can fuel conservation programs . PLoS ONE 4 : e7336 . (doi:10.1371/journal.pone. 0007336 ). 105. Eckert CG , Samis KE , Lougheed SC ( 2008 ) Genetic variation across species' geographical ranges: the central-marginal hypothesis and beyond . Molecular Ecology 17 : 1170 - 1188 . 106. Septon G ( 2000 ) Overwintering by urban-nesting peregrine falcons Falco peregrinus in Midwestern North America . In: Chancellor RD, Meyburg B-U, eds. Raptors at Risk. Surrey, B.C.: Hancock House Publications , Ltd. pp 455 - 461 . 107. Green MG , Swem T , Morin M , Mesta R , Klee M , et al. ( 2006 ) Monitoring results for breeding American Peregrine Falcon (Falco peregrinus anatum) , 2003 . Washington DC: U.S. Department of Interior, Fish and Wildlife Service, Biological Technical Publication FWS/BTP-R1005-2006. 108. Cade TJ , Martell M , Redig P , Septon GA , Tordoff HB ( 1996 ) Peregrine falcons in urban North America . In: Bird DM, Varland DE , Negro JJ, eds. Raptors in Human Landscapes . London: Academic Press. pp 3 - 14 . 109. Wakamiya SM , Roy CL ( 2009 ) Use of monitoring data and population viability analysis to inform reintroduction decisions: peregrine falcons in the Midwestern United States . Biological Conservation 142 : 1767 - 1776 . 110. Edmands S ( 2007 ) Between a rock and a hard place: evaluating the relative risks of inbreeding and outbreeding for conservation and management . Molecular Ecology 16 : 463 - 475 . 111. Temple SA , Cade TJ ( 1988 ) Genetic issues associated with recovery efforts for three endangered raptors . In: Garcelon DK, Roemer GW, eds. Proceedings of the International Symposium on Raptor Reintroduction , 1985 . Arcata, CA: Institute of Wildlife Studies . pp 17 - 29 . 112. Hansson B , Bensch S , Hasselquist D , Lillandt BG , Wennerberg L , et al. ( 2000 ) Increase of genetic variation over time in a recently founded population of great reed warblers (Acrocephalus arundinaceus) revealed by microsatellites and DNA fingerprinting . Molecular Ecology 9 : 1529 - 1538 . 113. Keller LF , Jeffery KJ , Arcese P , Beaumont MA , Hochachka WM , et al. ( 2001 ) Immigration and the ephemerality of a natural population bottleneck: evidence from molecular markers . Proceedings of the Royal Society of London Series BBiological Sciences 268 : 1387 - 1394 . 114. Vila C , Sundqvist AK , Flagstad O , Seddon J , Bjornerfeldt S , et al. ( 2003 ) Rescue of a severely bottlenecked wolf (Canis lupus) population by a single immigrant . Proceedings of the Royal Society of London Series B-Biological Sciences 270 : 91 - 97 . 115. Ortego J , Aparicio JM , Calabuig G , Cordero PJ ( 2007 ) Increase of heterozygosity in a growing population of lesser kestrels . Biology Letters 3 : 585 - 588 . 116. Luikart G , Cornuet J-M , Allendorf FW ( 1999 ) Temporal changes in allele frequencies provide estimates of population bottleneck size . Conservation Biology 13 : 523 - 530 . 117. Spencer CC , Neigel JE , Leberg PL ( 2000 ) Experimental evaluation of the usefulness of microsatellite DNA for detecting demographic bottlenecks . Molecular Ecology 9 : 1517 - 1528 . 118. Pertoldi C , Hansen MM , Loeschcke V , Madsen AB , Jacobsen L , Baagoe H ( 2001 ) Genetic consequences of population decline in the European otter (Lutra lutra): an assessment of microsatellite DNA variation in Danish otters from 1883 to 1993 . Proceedings of the Royal Society Series B 268 : 1775 - 1781 . 119. Martnez-Cruz B , Godoy JA , Negro JJ ( 2007 ) Population fragmentation leads to spatial and temporal genetic structure in the endangered Spanish imperial eagle . Molecular Ecology 16 : 477 - 486 . 120. Segelbacher G , Manel S , Tomiuk J ( 2008 ) Temporal and spatial analyses disclose consequences of habitat fragmentation on the genetic diversity in capercaillie (Tetrao urogallus) . Molecular Ecology 17 : 2356 - 2367 . 121. Demandt MH ( 2010 ) Temporal changes in genetic diversity of isolated populations of perch and roach . Conservation Genetics 11 : 249 - 255 . 122. Lott CA ( 2006 ) A new raptor migration monitoring site in the Florida Keys: counts from 1999-2004 . Journal of Raptor Research 40 : 200 - 209 . 123. Farmer CJ , Goodrich LJ , Ruelas Inzunza E , Smith JP ( 2008 ) Conservation status of North America's birds of prey . In: Bildstein KL, Smith JP , Ruelas Inzunza E , Veit RR, eds. State of North America's Birds of Prey. Series in Ornithology 3. Cambridge, MA: Nuttall Ornithological Club: Washington DC: American Ornithologists Union . pp 303 - 415 . 124. Smith JP , Farmer CJ , Hoffman SW , Lott CA , Goodrich LJ , et al. ( 2008 ) Trends in autumn counts of migratory raptors around the Gulf of Mexico , 1995 - 2005 . In: Bildstein KL, Smith JP , Ruelas Inzunza E , Veit RR, eds. State of North America's Birds of Prey. Series in Ornithology 3. CambridgeMA: Nuttall Ornithological Club: Washington DC: American Ornithologists Union . pp 253 - 277 . 125. Henny CJ , Yates MA , Seegar WS ( 2009 ) Dramatic declines of DDE and other organochlorines in spring migrant peregrine falcons from Padre Island , Texas, 1978 - 2004 . Journal of Raptor Research 43 : 37 - 42 . 126. Poulsen NA , Nielsen EE , Schierup MH , Loeschcke V , Grnkjaer P ( 2006 ) Long-term stability and effective population size in North Sea and Baltic Sea cod (Gadus morhua) . Molecular Ecology 15 : 321 - 331 . 127. Antao T , Perez-Figueroa A , Luikart G ( 2010 ) Evolutionary Applications (in press) Early detection of population declines: high power of genetic monitoring using effective population size estimators . DOI: 10.1111/j .1752- 4571 . 2010 .00150.x. 128. Ovenden JR , Peel D , Street R , Courtney AJ , Hoyle SD , et al. ( 2007 ) The genetic effective and adult census size of an Australian population of tiger prawns (Penaeus esculentus) . Molecular Ecology 16 : 127 - 138 . 129. Morin PA , Luikart G , Wayne RK , SNP workshop group (2004) SNPs in ecology, evolution and conservation . Trends in Ecology and Evolution 19 : 208 - 216 . 130. Morin PA , Martien KK , Taylor BL ( 2009 ) Assessing statistical power of SNPs for population structure and conservation studies . Molecular Ecology Resources 9 : 66 - 73 . 131. Backstrom N , Fagerberg S , Ellegren H ( 2008 ) Genomics of natural bird populations: a gene-based set of reference markers evenly spread across the avian genome . Molecular Ecology 17 : 964 - 980 . 132. Caballero IC , Ashley MV , Bates JM ( 2006 ) Return of the peregrine: what microsatellites can tell us about population genetics dynamics in reintroduced peregrine falcons . Forth North American Ornithological Conference, Veracruz, Mexico. 133. Dawnay N , Ogden R , Wetton JH , Thorpe RS , McEwing R ( 2009 ) Genetic data from 28 STR loci for forensic individual identification and parentage analyses in 6 bird of prey species . Forensic Science International: Genetics 3 : e63 - e69 . 134. United States Fish and Wildlife Service ( 2008 ) Take of migrant peregrine falcons in the United States for use in falconry . Federal Register 73 ( 236 ): 74508 - 9 . 135. Frankham R ( 1995 ) Effective population-size adult-population size ratios in wildlife - a review . Genetics Research 66 : 95 - 107 . 136. Caballero A ( 1994 ) Developments in the prediction of effective population size . Heredity 73 : 657 - 679 . 137. Wang J , Caballero A ( 1999 ) Developments in predicting the effective size of subdivided populations . Heredity 82 : 212 - 226 . 138. Turner TF , Osborne MJ , Moyer GR , Benavides MA , Alo D ( 2006 ) Life history and environmental variation interact to determine effective population to census size ratio . Proceedings of the Royal Society Series B 273 : 3065 - 3073 . 139. Newton I ( 1988 ) Commentary - Population regulation in Peregrines: an overview . In: Cade TJ, Enderson JH , Thelander CG , White CM, eds. Peregrine Falcon Populations: Their Management and Recovery . BoiseIdaho : The Peregrine Fund . pp 761 - 770 . 140. Mora M , Skiles R , McKinney B , Paredes M , Buckler D , et al. ( 2002 ) Environmental contaminants in prey and tissues of the peregrine falcon in the Big Bend Region , Texas, USA. Environmental Pollution 116 : 169 - 176 . 141. Mora M , Skiles RS , Paredes M ( 2007 ) Further assessment of environmental contaminants in avian prey of the peregrine falcon in Big Bend National Park , Texas. The Southwestern Naturalist 52 : 54 - 59 . 142. Maseveu OB ( 2009 ) Is illegal take affecting the peregrine falcon population in central Peru: why and how ? In: Rich TD, Arizmendi C , Demarest DW , Thompson C, eds. Tundra to Tropics : Connecting Birds, Habitats and People. Proceedings of the Fourth International Partners in Flight Conference , 13 - 16 February 2008 . McAllen, TX: Partners in Flight . pp 171 - 178 . 143. Soule ME , Estes JA , Miller B , Honnold DL ( 2005 ) Strongly interacting species: conservation policy , management, and ethics. BioScience 55 : 168 - 176 .

This is a preview of a remote PDF:

Jeff A. Johnson, Sandra L. Talbot, George K. Sage, Kurt K. Burnham, Joseph W. Brown, Tom L. Maechtle, William S. Seegar, Michael A. Yates, Bud Anderson, David P. Mindell. The Use of Genetics for the Management of a Recovering Population: Temporal Assessment of Migratory Peregrine Falcons in North America, PLOS ONE, 2010, DOI: 10.1371/journal.pone.0014042