Swainson's Warbler (Limnothlypis swainsonii) is a species of conservation concern because of its small wintering range in the Caribbean Basin, relatively low population densities, and habitat fragmentation in its core breeding range in the southeastern United States. We investigated microsatellite DNA variation among 11 breeding populations from eastern Texas to Virginia and two populations from wintering areas in Jamaica and Mexico. Analyses of six polymorphic loci indicated a moderate level of gene flow among breeding populations, relatively small effective population sizes (<200 individuals in each sampled population), and subtle population variation. We detected no evidence of population bottlenecks in breeding or wintering populations. Bayesian assignment tests suggested that substantial mixing of breeding populations may occur in wintering areas. Genetic differences between the Mexican and Jamaican populations indicate they may be drawn from different subsets of breeding populations. Patterns of genetic variation among breeding and wintering populations suggest a network of local and regional conservation programs may be necessary to maintain genetic diversity in Swainson's Warbler.
Genetic differentiation in migratory species is often associated with migratory divides where populations that winter in different regions meet during the breeding season in parapatric contact zones (Salomonsen 1955). Studies of Sylvia warblers (Sylviidae) in Europe have demonstrated that migratory behaviors can have a strong genetic basis and can evolve rapidly (Berthold and Querner 1981; Helbig 1991, 1996; Berthold 2003; Bearhop et al. 2005). Experimental crosses between populations that winter in different regions produced offspring that exhibited intermediate degrees of migratory orientation and restlessness. Helbig (1991) suggested that hybrids between populations of Eurasian Blackcap (Sylvia atricapilla) from opposite sides of the migratory divide would be selected against, because they would probably attempt to migrate across the Mediterranean and end in the uninhabitable expanses of the Sahara Desert. This example suggests a plausible mechanism for maintenance of genetic structure in breeding populations of migratory species that winter in different areas (allohiemy).
Swainson's Warbler (Limnothlypis swainsonii) is a sparsely-distributed wood warbler (Parulidae) that breeds in the unglaciated southeastern United States and winters in the Caribbean Basin (Meanley 1971; Brown and Dickson 1994; Graves 2001, 2002). Geographic variation in plumage color is subtle and the species is currently regarded as monotypic (Brown and Dickson 1994). Although locally common, Swainson's Warbler has been ranked as one of the most vulnerable breeding songbirds in the southeastern United States because of habitat fragmentation, relatively low population densities, and a small disjunct wintering range (Morse 1989, Terborgh 1989, Hunter et al. 1993, Rappole 1995). Most contemporary breeding populations occur in second-growth forest on alluvial soils in the Mississippi Valley and on the coastal plain from eastern Texas to southeastern Virginia (Graves 1998, 2001, 2002). The nonbreeding distribution is poorly known (Brown and Dickson 1994, Graves 1996), but the primary wintering areas appear to be in Cuba (Kirkconnell et al. 1996), Jamaica (Graves 1996), and in humid forests of mainland Middle America from Veracruz to Tabasco in Mexico (Winker et al. 1992, Winker et al. 1999b) and in Belize (K. Winker, unpubl. data).
Winker et al. (2000) compared allozyme variation at 26 loci among five breeding populations of Swainson's Warbler. Allele frequencies at five loci indicated modest population structure, chiefly between samples in the Mississippi Valley from Arkansas and populations on the coastal plain from Louisiana to Virginia; patterns of genetic variation were inconsistent with an isolation-by-distance model. It was hypothesized the observed population structure was due to genetic drift due to the absence of barriers to gene flow in the contemporary landscape. Alternatively, genetic structuring could be influenced by the warbler's divided wintering range (Greater Antilles vs. mainland Middle America). If migratory behavior is under genetic control in this nocturnally migrating species, interbreeding between populations that winter on the mainland and those that winter in the Greater Antilles might produce offspring that end their migrations over open water in the Caribbean or in the Gulf of Mexico.
We investigated allelic variation at six polymorphic microsatellite loci in 11 breeding populations of Swainson's Warbler sampled from eastern Texas to Virginia and two populations from wintering areas in Jamaica and Mexico. We had three principal objectives: (1) characterize microsatellite diversity within breeding and wintering populations, (2) examine the evidence for mixing of breeding populations on wintering areas, and (3) contrast microsatellite and allozyme variation in the subset of breeding populations studied by Winker et al. (2000).
METHODS
Data Collection
Swainson's Warblers are territorial in breeding (Meanley 1971) and wintering areas (Graves 1996); they are also monogamous, although polygyny may occur (Graves 1992), and pairing occurs in breeding areas. We obtained tissue samples from 205 territorial males from breeding locations between 27 April and 31 May (1986–1996) in the southeastern United States (Fig. 1) under state and national permits: Sulphur River, Texas (TX1: n = 17); Sam Houston National Forest, Texas (TX2: n = 20); Atchafalaya River, Louisiana (LA1: n = 22); Tensas River, Louisiana (LA2: n = 10); Homochitto River, Mississippi (MS1: n = 20); Little Sunflower River, Mississippi (MS2: n = 10); White and Mississippi rivers, Arkansas (AR: n = 20); Apalachicola River, Florida (FL: n = 24); Ocmulgee River, Georgia (GA: n = 20); Cooper and Santee rivers, South Carolina (SC: n = 21); and the Great Dismal Swamp, Virginia, and Chowan River, North Carolina (VA: n = 21). These samples nearly span the latitudinal and longitudinal distribution of the core breeding range (Graves 2002). Individuals from five population samples (LA1, AR, FL, SC, and VA) were previously assayed in the allozyme study of Winker et al. (2000). Genetic samples from the wintering range were obtained in Veracruz (n = 25) and Tabasco (n = 1) in southern Mexico (MEX) and in the Blue Mountains of Jamaica (JAM: n = 19; Fig. 1). Voucher specimens (except for Jamaican samples, for which only blood was obtained) are housed in the National Museum of Natural History (Smithsonian Institution), Bell Museum of Natural History (University of Minnesota), and the Coleccion Nacional de Aves (Universidad Nacional Autonoma de Mexico).
Genomic DNA was extracted from body tissues or from small pieces of museum skins using a diatomaceous earth/guanidine thiocyanate extraction protocol (Carter and Milton 1993) and diluted to 20 ng/μL for amplification of microsatellite loci developed for this species by Winker et al. (1999a). Extractions and amplifications were conducted in a PCR-free laboratory. All amplification batches included negative controls. Individuals were randomized across extractions, amplifications, and gel runs. Fragments were generated using dye-labeled primers and visualized on an ABI 373A automated sequencer. Fragment size was measured using an internal size standard (350 TAMRA) and GeneScan software (both from Applied Biosystems Inc., Foster City, CA, USA). All electropherograms were examined manually, allele scoring was done without reference to source population, and results were not sorted to population until completion of the study.
Analyses
Basic statistics of allelic frequencies, genetic distance measures between populations, allelic diversity, and expected and observed heterozygosities were calculated using the computer programs BIOSYS-1 and GDA version 1.0 (Swofford and Selander 1981, Lewis and Zaykin 1999). Tests for Hardy-Weinberg (H-W) equilibrium, allelic heterogeneity, levels of population structure, and differences between population pairs were performed using GENEPOP version 3.1d (Raymond and Rousset 1995). Levels of population structure are given using the θ of Weir and Cockerham (1984), which is an unbiased estimator of traditional FST. We used CERVUS version 1.0 (Marshall et al. 1998) to infer frequencies of null alleles. We used analysis of covariance (ANCOVA) to evaluate population variability in allelic diversity. The relationship between genetic and geographic distances among populations was assessed with Mantel and permutation tests using NTSYS-pc version 1.50 (Rohlf 1988). We examined geographic subsets of the population matrix with Spearman rank correlation coefficients (rs). Pairwise comparisons of population parameters are not independent, because each population is subjected to multiple contrasts. Therefore, we emphasize the relative strength of the correlative relationship (rs) between genetic and geographic distances for these populations rather than P values.
We used STRUCTURE, version 2 (Pritchard et al. 2000, Falush et al. 2003) to examine how well the predefined populations corresponded to genetic groups (K). Individual genotypes in this Bayesian clustering approach are assigned to clusters with Hardy-Weinberg equilibrium and linkage equilibrium achieved within each cluster. A Markov Chain Monte Carlo approach was used to identify the number of clusters (K) that are most likely given the observed genotypes. STRUCTURE was used twice for each defined K (1–11) after a burn-in of 105 iterations, followed by an additional 106 iterations on the full breeding area data set (11 populations). No prior information (e.g., on the population of origin of each individual) was used. Admixture and correlations models were used in which individuals can have mixed ancestry, and where the allele frequencies of closely related populations may be correlated. We further refined our analysis of genetic groups (K) using the approach of Evanno et al. (2005), which examines the second order rate of change of the log probability of data with respect to the number of clusters. These analyses were based on 20 independent STRUCTURE runs for each K under the same conditions of a 105 burn-in, plus another 106 iterations. A bootstrapped (100 replicate), neighbor-joining tree was developed using SEQBOOT, GENDIST, NEIGHBOR, and CONSENSE in the software package PHYLIP (Felsenstein 1993).
We tested for population bottlenecks using the computer program M (Garza and Williamson 2001). This test calculates the ratio (M) of total alleles to the range of allele sizes within populations. Reduced populations are expected to have a smaller M ratio than populations in mutation-drift equilibrium (Garza and Williamson 2001). Simulations estimating the probability of the observed M ratio used Garza and Williamson's (2001) suggested values for the proportion of one-step mutations (90%) and the average size of non-one-step mutations (Δg = 3.5). We used 10,000 replicates for simulations and set θ (4Neμ) to 1, 10, and 25, values corresponding to equilibrium effective population sizes (Ne) before a bottleneck of 500, 5,000, and 12,500, respectively, at a mutation rate (μ) of 5 × 10−4 (Goldstein and Schlotterer 1999). Population sizes of Swainson's Warbler are unknown, but values of θ well above 1 seem reasonable given our field experience.
Estimates of 4Neμ were obtained using MIGRATE (Beerli and Felsenstein 2001), where Ne is effective population size and μ is the microsatellite mutation rate. Our estimates of long-term effective population sizes (Ne) were made using a mutation rate (μ) of 5 × 10−4 (Goldstein and Schlötterer 1999). Estimates of levels of gene flow among breeding populations were made using the rare alleles method (Slatkin 1985, Barton and Slatkin 1986, Slatkin and Barton 1989) as implemented in GENEPOP version 3.1d (Raymond and Rousset 1995) and assignment tests; the latter (Cornuet et al. 1999) enable a more direct estimate of gene flow by calculating the likelihood that an individual genotype originated from the population where it was obtained. These tests remove an individual from the sample and compare it with the remainder (Rannala and Mountain 1997, Cornuet et al. 1999) to identify the most likely breeding population of origin for wintering individuals. We used an α of 0.01, a stringency level shown to have high levels of accuracy in assigning populations of origin for dispersing individuals (Berry et al. 2004). We used the Dunn-Sidák correction (Sokal and Rohlf 1995) where appropriate to adjust probabilities for simultaneous statistical tests.
RESULTS
Variability of Microsatellite DNA
The six microsatellite loci exhibited 5–16 alleles per locus (Table 1). Sixty-seven alleles were detected across all loci, ranging from a minimum of 32 detected in the MS2 population (n = 10 individuals) to a maximum of 45 in the LA1 population (n = 22 individuals). Ten alleles were unique to a single population (FL = 3; JAM = 2; LA2 = 1; MEX = 2; TX2 = 2). The number of alleles detected per population was correlated (r2 = 0.54; P < 0.005) with sample size. Populations from localities that drain into the Atlantic Ocean (VA, SC, GA) had significantly lower allelic diversity (ANCOVA, F1,8 = 9.57; P = 0.015) than populations from drainages that empty into the Gulf of Mexico after factoring out the effects of sample size (FL, MS1, MS2, LA1, LA2, AR, TX1, TX2; Table 2). An overall heterozygote deficiency (P = 0.011) was caused by loci Lswμ5B and Lswμ19 (Table 3). Heterozygote deficiency can be caused by a number of factors, including inbreeding, the Wahlund effect, and null alleles. Inferred frequencies for the latter possibility suggested that null alleles were not responsible for the H-W disequilibrium (Table 3).
Population Structure
Breeding populations had significant but weak structure (θ = 0.0068, P < 0.00005; Table 4) where θ is the unbiased estimator for FST; this was driven by heterogeneous distributions of alleles in loci Lswμ5B (P < 0.00005) and Lswμ18 (P = 0.039; Table 4). The Bayesian assessment of population clusters reflected this weak structure with the highest probabilities of K occurring for 1–7 populations (Fig. 2A). Further testing to identify the true number of genetic clusters (K), performed separately for both breeding and wintering populations, suggested that 2–4 populations are involved (2–4 among breeding populations and 2–3 in the wintering samples), although the peak height of delta K in both sets of analyses suggested a lack of strong signal in the data set (Fig. 2B, C). A bootstrapped distance tree (10,000 replicates) provided little support for any breeding population associations; just two populations formed one supported clade (FL and MS1 supported by 60% of bootstrap replicates; not shown).
Mantel and permutation tests (NTSYS-pc) of the pooled breeding population data set indicated that genetic distance was at best only weakly associated with geographic distance (Z = 0.27, P = 0.11; 10,000 random permutations). However, when population contrasts were partitioned by geographic region, a clearer view of population differentiation emerged (Figs. 1, 3). Pairwise comparisons between Gulf populations (n = 28 contrasts) indicated a slightly negative relationship between genetic and geographic distance (rs = −0.39). The relationship between genetic and geographic distance was substantially weaker (rs = −0.23) when pairwise comparison was limited to the seven populations from the Mississippi and Sabine river drainages (n = 21 contrasts) of the western Gulf. In contrast, pairwise contrasts (n = 24) between populations from Atlantic drainages with those from Gulf drainages exhibited a markedly positive correlation (rs = 0.49), suggesting a weak longitudinal trend toward genetic isolation-by-distance. The most geographically isolated population in our survey (VA) also exhibited the greatest average pairwise genetic distance from other populations. The most centrally located population (FL) also showed the lowest average genetic distance between populations.
Wintering populations exhibited significant but weak structure (θ = 0.0111, P = 0.044; Table 4). This was driven largely by locus Lswμ9 (P = 0.004; other loci P > 0.098; Table 4). The loci responsible for breeding and wintering population structure were different (Table 4).
Allozyme and Microsatellite Comparisons
We compared genetic distance matrices derived from allozymes and the microsatellite loci examined in this study for five breeding populations (LA1, AR, FL, SC, VA). The null hypothesis that the two matrices were not associated could not be rejected with Mantel and permutations tests (Z = 0.18, P = 0.30; 10,000 random permutations) indicating the two nuclear genetic marker systems examined for this subset of breeding populations exhibited discordant patterns of differentiation.
Gene Flow
The rare alleles method, an indirect technique to estimate the number of migrants per generation, inferred Nm = 9.16 among breeding populations after correction for sample size (where mean sample size was 18.6 and mean frequency of private alleles P[1] = 0.026). Assignment tests, which are more reflective of contemporary gene flow, indicated that three to seven individuals (15– 70%) from each breeding population had genotypes that were significantly unlikely to have originated there (averaging 4.4 individuals or 25% per population; Table 5). These values ranged from two to five individuals per population (averaging 3.2 individuals or 18% per population; Table 5) when adjusted for type-I error (Dunn-Sidák test). These numbers were similar in the two wintering populations (averaging 4.5 individuals per population, or 3.5 with the Dunn-Sidák correction; Table 5). The genotypes of three individuals from Florida and single individuals from Mississippi (MS1), South Carolina, and Texas (TX2) could not be assigned to any breeding population (all with P < 0.00005). Similarly, two individuals from Jamaica and one from Mexico could not be assigned to any breeding or wintering population. Both the rare alleles method and assignment tests suggested a moderate amount of gene flow among breeding populations (Nm = 3.2–9.2).
Long-term Population Size
We obtained estimates of 4Neμ, where Ne is effective population size and μ is the microsatellite mutation rate. These analyses suggested uniformly small long-term effective population sizes (Ne) ranging from a low of 44 (LA2) to a high of 174 (LA1; Table 6). There was no evidence for population genetic bottlenecks in breeding or wintering populations at any of the modeled values of θ (1, 10, and 25; P > 0.1), despite these rather low estimates of long-term effective population size.
Evidence of Winter Mixing
We compared the allelic frequencies of wintering populations from Mexico and Jamaica with those of each of the 11 breeding populations to investigate the genetic relationships between breeding and wintering populations. No genetic differences were found between any pairwise combination of wintering and breeding population when α was adjusted for multiple tests (Table 7). Evidence of geographic segregation of breeding populations in wintering areas remained weak under a less conservative approach (no adjustment of α). Allelic frequencies of the Mexican wintering population were significantly different from those of breeding populations from Arkansas (AR) in the Mississippi Valley and Virginia (VA) and South Carolina (SC) on the Atlantic coastal plain (Table 7). Similarly, allelic frequencies of the Jamaican wintering population were significantly different from those of breeding populations from Louisiana (LA2), Arkansas (AR), Texas 1 (TX1), and Virginia (VA; all P < 0.02; Table 7).
We assigned wintering individuals to their most likely breeding population by choosing the highest probability P-value from the array of breeding populations. All but two individuals from Mexico and two from Jamaica could be assigned to one of the 11 breeding populations (P < 0.01). Of the 24 assignable individuals from Mexico, only Arkansas (AR) and Georgia (GA) were eliminated as likely breeding populations; and among the 16 assignable individuals from Jamaica, only Louisiana (LA2), Mississippi (MS2), and Virginia (VA) were unlikely breeding sources. This suggests that individuals from a substantial number of geographically distinct breeding populations mix in wintering areas in Mexico and Jamaica.
DISCUSSION
The breeding distribution of Swainson's Warbler has decreased significantly during the past century because of deforestation of bottomlands, flood mitigation projects, and forestry practices that minimize early successional habitats (Twedt and Loesch 1999; Graves 2001, 2002). Several populations near the northern periphery of the warbler's breeding range have disappeared in recent decades and many of the remaining populations are isolated. We detected no evidence of genetic bottlenecks in breeding or wintering populations despite relatively small effective population sizes and fragmentation of the warbler's geographic range. However, the lower allelic diversity exhibited among Atlantic-drainage populations and the significant correlation between genetic and geographic distances observed in pairwise comparisons of Gulf-drainage populations with those sampled in Atlantic drainages suggests that subtle population structure exists among breeding populations. Clustering tests further suggested that breeding populations sampled in this study are composed of two to four genetic populations. Whether differentiation is caused by genetic drift, stochastic sampling effects, selection associated with the warbler's disjunct wintering range, or other factors is unknown.
The moderate levels of gene flow among breeding populations revealed by microsatellite data (Nm = 3.2–9.2) were consistent with the results of an earlier survey of allozymes (Winker et al. 2000) in a subset of five populations (Nm = 1.5–11.7). However, spatial patterns of genetic variation revealed by microsatellites and allozymes were discordant, a not uncommon event between molecular marker systems (Allendorf and Seeb 2000). Significant population structure was caused by two microsatellite loci in the present study and three of 16 loci examined in our earlier allozyme study.
Mixing of breeding populations in wintering areas is believed to be a common phenomenon because breeding ranges of most Nearctic-Neotropic migratory songbirds are considerably larger than their wintering ranges (Terborgh 1989). The signal of geographic structure in our microsatellite data is weak, but Bayesian assignment tests imply that at least six breeding populations (TX1, TX2, LA1, MS1, FL, and SC) could have contributed wintering individuals to both the Mexican and Jamaican wintering populations. However, genetic differences between the Mexican and Jamaican populations suggest they are comprised of different subsets of breeding populations.
CONSERVATION IMPLICATIONS
Should breeding populations of Swainson's Warbler be managed as a single conservation unit? Patterns of microsatellite variation among breeding and wintering populations suggest that multi-regional and international efforts will likely be required to maintain the current level of genetic diversity in the species. Implementation of a viable management plan for the species throughout its global range, however, is contingent upon the development of additional genetic markers that can more fully resolve the genetic structure of breeding populations and their distributions during the nonbreeding season. Future field efforts should focus first on obtaining genetic samples from breeding populations in the Appalachian and Ozark mountains and from wintering populations in Cuba and Belize.
Acknowledgments
The Legacy Resource Management Program of the U.S. Department of Defense (DAMD17-93-J-3073, DACA87-94-H-0012) and the Research Opportunities Fund of the Smithsonian Institution supported this study. Numerous agencies kindly granted permits for the research: Jamaican Conservation and Development Trust; National Resources Conservation Authority (Kingston); USDA Forest Service; U.S. Fish and Wildlife Service; U.S. Department of Defense; Arkansas Game and Fish Commission; Florida Fish and Wildlife Conservation Commission; Georgia Department of Natural Resources; Louisiana Department of Wildlife and Fisheries; Mississippi Department of Wildlife, Fisheries, and Parks; North Carolina Wildlife Resources Commission; South Carolina Department of Natural Resources; Texas Parks and Wildlife Department; and Virginia Department of Game and Inland Fisheries. R. T. Brumfield and T. C. Glenn helped with marker development. R. M. Zink permitted us to sample study skins in the Bell Museum of Natural History. Patricia Escalante P provided additional samples from Mexico. B. K. Schmidt, C. M. Milensky, J. P. Angle, and J. P. Dean assisted with specimen preparation. Francois Rousset provided helpful correspondence on making multiple comparisons. M. J. Braun, R. T. Brumfield, P. O. Lewis, and C. L. Pruett provided helpful comments on an earlier draft.
LITERATURE CITED
TABLE 1. Allele frequencies for six microsatellite loci among 11 breeding and two wintering populations of Swainson's Warbler (n = 250).
TABLE 1. Continued.
TABLE 2. Sample sizes (n), average allelic diver sity (A), and expected and observed heterozygosities (He and Ho) for six microsatellite loci in breeding and wintering populations of Swainson's Warbler
TABLE 3. Exact tests for Hardy-Weinberg equi librium (P-values; GENEPOP, Raymond and Rousset 1995) and inferred frequencies of null alleles (CER VUS, Marshall et al. 1998)
TABLE 4. Allelic heterogeneity and levels of pop ulation structure (θ and associated P-values; Weir and Cockerham 1984) among 11 breeding and between two wintering populations of Swainson's Warbler (GENEPOP, Raymond and Rousset 1995)
TABLE 5. Assignment test results comparing ge notypes of individuals with genotypic characteristics of the breeding populations from which they were sampled (GeneClass, Cornuet et al. 1999)
TABLE 6. Estimates of 4Neμ and the 95% confidence interval (following Beerli and Felsenstein 2001) where Ne is effective population size and μ is the mutation rate of microsatellite loci. Subsequent estimates of long-term effective population sizes (Ne est.) and the corresponding 95% confidence intervals are based on a mutation rate of 5 × 10−4 (Goldstein and Schlötterer 1999) and are well below contemporary census sizes of these populations
TABLE 7. Pairwise comparisons showing genetic similarities and differences between two wintering and 11 breeding populations of Swainson's Warbler (θ and associated P-values; GENEPOP)