AN ABSTRACT OF THE THESIS OF
Leah R. Gorman for the degree of Master of Science in Wildlife Science presented on
June 6 2000. Title: Population Differentiation among Snowy Plovers (Charadrius
in North America.
Abstract approved: Susan M. Haig
Two North American subspecies of Snowy Plovers (Charadrius alexandrinus) have been described: the Cuban
Snowy Plover (C. a. tenuirostris) and the Western Snowy Plover (C. a. nivosus). Coastal populations of the Western Snowy Plover
are listed as Threatened under the U.S. Endangered Species Act, while populations ofthe Cuban Snowy Plover in the southeastern
U.S. and Puerto Rico show evidence of decline and are being considered for listing. To clarify the relationships among populations,
we examined variation in mitochondrial (mtDNA) sequences and at inter-simple sequence repeat (ISSR) loci among 8 populations
distributed across North America and the Caribbean. MtDNA d-loop region sequences (322 bp,n = 126) revealed 8 haplotypes,
with one haplotype unique to Puerto Rico, and indicated that Puerto Rico is significantly differentiated from all mainland
populations. Conversely, neighbor-joining analysis of 16 ISSR loci suggests that Puerto Rico is nested within a cluster ofpopulations
from eastern continental North America. Evidence for structure among continental populations was weak among mtDNA hapiotypes
(cD 5 t = 0.025, p = 0.178), but analysis ofISSR markers supported subdivision into groups east and west of the Rocky Mountains
('DST = 0.445, p <0.001). A significant relationship between genetic distance (pairwise FST) and geographic distance was
observed (Mantel r = 0.85, p c 0.01). Overall,
two management units in continental North America (divided by the Rocky Mountains) and a separate unit for the Puerto Rican
population may be warranted. We recommend that both genetic and demographic considerations be weighed in making policy decisions
regarding the status and listing of Snowy Plover population segments.
Population Differentiation among Snowy Plovers
in North America
Leah R. Gorman
A Thesis Submitted
Oregon State University
In Partial Fulfillment of the
requirements for the degree of Master of Science
Presented June 6, 2000-Commencement June 2001
Master of Science thesis of Leah R. Gorman presented on June 6,2000.
M. Haig, representing Wildlife Science
Chair of Department of Fisheries and Wildlife
Dean of the Graduate School
I understand that my thesis will become part of the permanent collection of Oregon State University libraries.
My signature below authorizes release of my thesis to any reader upon request.
Leah R. Gorman,
I would like to express my gratitude to my major professor
Sue Haig for her constant support in the design, implementation, analysis, and writing of this thesis. I am also grateful
for the valuable input provided by my committee members Joe Beatty, Aaron Liston, and Dan Roby. Thanks to C. Hallett, P. Sanzenbacher,
and 0. Taft for assistance in the field. Additional samples were provided by A. Archuleta, R. Estelle, A. Powell, M. Stern,
R. Swift, J. Whittier, and the Louisiana Museum o fNatural Sciences. The field component of the project benefitted from the
advice of R. Boyd, J. Collazo, J. Gore, P. Paton, J. Plissner, J. Rupert, F. Vilella, and C. Zonick. Logistical support was
provided by the USGS Forest and Rangeland Ecosystem Science Center, Cabo Rojo National Wildlife Refuge, the Florida Game and
Fresh Water Fish Commission, the Florida Division of Recreation and Parks, Gulf Island National Seashore, Eglin Air Force
Base, Texas Parks and Wildlife Division, Lower Rio Grande Valley National Wildlife Refuge, Quivira National Wildlife Refuge,
Inland Sea Shorebird Reserve, Bureau of Land Management Salt Lake District Office, Utah Division of Wildlife Resources, and
Summer Lake Wildlife Area. T. Mullins and S. Wagner offered guidance in the laboratory. J. Beatty, D. Kesler, A. Liston, D.
Roby, and S. Wagner provided helpful comments on the manuscript. The work was funded by grants from the USGS Forest and Rangeland
Ecosystem Science Center, USGS Biological Resources Division Species-at-Risk program, and the U.S. Fish and Wildlife Service.
Table of Contents
and Methods.................................... 4
Tissue Sampling.......................................... 4
DNA Amplification and Sequencing
Amplification, Screening, and Scoring of ISSR Markers.8
Data Analysis ...........................................10
MtDNA Sequence Data .....................................14
ISSR Markers ............................................22
Genetic Variation within Populations ....................30
Genetic Structure and Phylogeographic Patterns.......... 31
Relationships among Subspecies
Unit Designation ...........................36
List of Figures
||Breeding range of the Snowy Plover in North America, Caribbean, and Bahamas (modified
from Page eta!. 1995a) and sampling site locations
||Minimum spanning network reflecting the evolutionary relationships of 8 maternal
haplotypes of Snowy Plover
||Distribution of Snowy Plover haplotypes. Pie charts represent haplotype frequencies
||Neighbor-joining tree of North American and Caribbean Snowy Plover populations
based on Manhattan distances of 16 variable ISSR loci
||Non-metric multidimensional scaling based on Jaccard's distances for Snowy Plovers
(n=104) using 16 variable ISSR loci
||Scatterplot of FST estimates based on 16 variable ISSR loci versus geographic distance for all pairwise
comparisons of populations
List of Tables
||Variable nucleotide sites from the aligned 322 bp of the mitochondrial control
region for North American, Puerto Rican, and Peruvian Snowy Plover populations
|Haplotype (/i) and nucleotide (it) diversity for North American and Puerto Rican
Snowy Plovers based on mtDNA control region sequence
||Analysis of molecular variance for Snowy Plover populations based on mtDNA sequence
variation and 16 variable ISSR loci
||Probability of non-differentiation between populations based on mtDNA sequence
variation (below diagonal) and 16 variable ISSR loci (above diagonal)
||Estimates of ISSR dominant marker allele frequencies (÷)using Lynch and Milligan's
(1994) Taylor Expansion for 8 Snowy Plover breeding areas across North America and the Caribbean
||Genetic variability within populations of Snowy Plovers averaged over 16 ISSR loci
Population Differentiation among Snowy Plovers
alexandrinus) in North America
Snowy Plovers, or Kentish Plovers (Charadrius alexandrinus), are small sequentially polygamous shorebirds
with a global distribution that have suffered population declines due to human disturbance, habitat loss, and habitat degradation
(Page ci at 1995a). At least five subspecies of Kentish Plovers have been described based on morphological variation, although
the taxonomy is controversial (Hayman ci at 1986). In North America, clarity regarding subspecies definition and the relationship
among regional populations is important in determining where conservation efforts should be focused and which populations
are listed for protection under the U.S. Endangered Species Act.
There are two subspecies of Snowy
Plovers recognized in North America. The Western Snowy Plover (C. a. nivosus) is thought to comprise resident populations
on the Pacific Coast of the United States, Texas, and Mexico, as well as migratory populations in the Great Basin and Great
Plains (American Ornithologists Union 1957; Figure 1). The Cuban Snowy Plover (C. a. tenuirostris) is generally thought to
encompass resident populations on the United States Gulf Coast east of Louisiana, Greater Antilles, southern Bahama Islands,
Yucatan peninsula of Mexico, and islands off Venezuela. Often, C. a. tenuirostris is not mentioned in avian taxonomy discussions,
and all North American Snowy Plovers are considered nivosus (e.g., Hayman ci cii. 1986). No genetic or morphological examination
has been conducted to resolve these issues, aside from describing more or less pale plumage (Conover 1945, AOU 1957, Binford
1989, Page et al. 1995a). A third subspecies in the Western Hemisphere, C. a. occidentauis, is found in coastal Peru and Chile.
This subspecies is distinguished from the North American subspecies by a wider eye stripe and less sexual dimorphism in breeding
adults. Together, the three form a group known as the 'Snowy Plover', which has been distinguished from other races of Kentish
Plover by the all white lores of their breeding plumage (Hayman et at 1986).
The Western Snowy
Plover is one of a few species listed as a distinct population segment under the U.S. Endangered Species Act (Federal Register
1996). Coastal populations of this subspecies are federally listed as Threatened, while others are not (Federal Register 1992).
Many of the same pressures that affect populations of the Western Snowy Plover have impacted populations in the southeastern
United States (Chase & Gore 1989, Gore 1996, Sprandel et al. 1997). As a result, the Cuban Snowy Plover has been listed
as threatened by the state ofFlorida and endangered by the state of Alabama. Moreover, in Puerto Rico, the breeding population
only consists of approximately 20 pairs (Lee 1989). Should the U.S. Fish and Wildlife Service proceed with a federal listing
of these populations, it will be important to understand whether the Cuban Snowy Plover represents a distinct genetic unit
within the species, and should be managed separately from the Western Snowy Plover. Additionally, it will be important to
determine the distinctiveness of individual populations of the Snowy Plover in the United States, as threats to the viability
of populations vary geographically.
In this study, we describe the population structure and phylogeography
of Snowy Plovers in North America by examining variation in the mitochondrial DNA (mtDNA) control region and at inter-simple
sequence repeat (ISSR) loci. We analyze the relationship between the geographic distribution of genetic variation and gene
flow, genetic drift, and historic processes. We also discuss the implications for defining conservation units and conserving
genetic diversity within and among populations.
Materials and Methods
Following protocols outlined by the American Ornithologists
Union (Oring et at 1988), blood samples were collected from breeding populations in Puerto Rico, Florida, Kansas, Texas, Utah,
eastern Oregon, and California (Figure 1). Once Snowy Plover nests were located, standard shorebird nest traps were used to
catch adult birds as they returned to their nests (see Hill & Talent 1990). In order to minimize nest desertion, attempts
to trap birds were not made unless at least two eggs had been laid. Birds were marked with a U.S. Fish and Wildlife Service
aluminum band to insure that individuals were not recaptured. Blood samples were taken from the brachial vein using a 27 gauge
needle, collected in a 0.25 ml heparinizedcapillary tube, and transferred to a 1 niL cryogenic vial containing a buffer solution
(100 mM Tris HCI, pH 8.0; 100mM EDTA, pH 8.0; 10 mM NaCl; 0.5% SDS). Samples were stored at room temperature until they could
be put into permanent storage at -80°C.
In Oregon and Oklahoma, heart, pectoralis, and embryonic
tissue samples were collected from carcasses of individuals that died of natural causes. In addition, membranes were collected
from eggs in Colorado. Finally, tissue samples from six individuals of C. a. occidentalis from Peru were provided by the Louisiana
Museum of Natural Sciences (LMUMZB24O19, B-24029, B-24030, B-24047, B-24107, B-24l57). The dates of collection for the samples
from Peru are unknown.
Figure 1. Breeding range of the Snowy Plover in North America, Caribbean,
and Bahamas (modified from Page et al. 1995a) and sampling site locations. Dashed line indicates currently accepted subspecies
boundary between the Western Snowy Plover (C. a. tenuirostris) and Cuban Snowy Plover (C. a. nivosus). Sites sampled for C.
a. tenuirostris include: PUERTO Rico (PR): a. Cabo Rojo (n = 8); FLOPJDA (FL): b. Cayo Costa State Park (n = 2), c. St. George
Island State Park (n = 4), d. Philips Inlet (n = 1), e. Eglin Air Force Base (n = 7), f. Gulf Island National Seashore (n
= 5). Sites sampled for C. a. nivosus include: TEX.ks (TX): g. Lower Rio Grande Valley National Wildlife Refuge (n = 7); EASTERN
GREAT PLAINS (GP): h. Great Salt Plains National Wildlife Refuge, Oklahoma (n = 7); i. Quivira National Wildlife Refuge, Kansas
(n = 20); COLORADO (CO): j. Arkansas River (n = 7); UTAH (UT): k. Inland Sea Shorebird Reserve (n=12), I. Saltwell Flats (n
= 3); EASTERN OREGON (EO): m. Summer Lake (n = 6), n. Lake Abert (n = 2); OREGoN COAST (DC): o. Baker Beach, Lane County (n
= 1), p. Coos County (n = 11); CMiFoI~NL&(CA) q. San Diego County (n = 20).
High molecular weight DNA was obtained for analysis by a standard phenol/chloroform extraction from blood
and tissue samples. Briefly, 10ml ofblood or a 1 mm section of tissue was digested in 400 ml of extraction buffer A (50mM
Tris, pH 8.0; 10mM EDTA, 200mM NaCI, 2% SDS) and 30, ml Proteinase K (20mg/mi). Samples were vortexed and incubated overnight
(-18 h) at 50°C. If blood clots or tissue did not fully disperse, a second aliquot of Proteinase K was added, and samples
were incubated for an additional 2 h. Samples were extracted with an equal volume of phenol (saturated with 10mM Tris, pH
8.0) followed by an equal volume of chloroforin/isoamyl alcohol (25:1). DNA was precipitated by adding a 1/10 volume of 3M
sodium acetate, two volumes of cold 95% ethanol, and pelleted at approximately 15,000 xG for20 minutes. Pellets were rinsed
once with 70% ethanol, dried under vacuum, and resuspended in 30j4 sterile water. DNA concentration of samples was quantified
with a DyNA Quant 200 fluorometer (Hoefer).
DNA Amp1~ficationand Sequencing
The polymerase chain reaction (PCR) was used to amplify a 322 bp fragment of the niitochondrial control
region, or d-loop, using the primer pair T5437L and TS778H (Wenink ci al. 1994). Each reaction was carried out using 100 ng
of genomic DNA in a 50 ml volume. Fragments were amplified using a PTC 100 thermal cycler (MJ Research). An initial denaturation
of 10 miii at 93°C was followed by 44 cycles of 1 mm denaturation at 92°C,1 mm annealing at 45°C, and 2 mm extension at 72°C,
with a single final extension at 72°Cfor 10 mm. The reactions were purified electrophoretically; bands of approximately 350
bp were extracted from a 1% agarose gel using an ultra-free-mc 0.45 filter (Millipore) and purified over a micron-50 filter
(Millipore). Sequence data were generated with an automated sequencer (Applied Biosystems model 373A) located in the Central
Services Laboratory at Oregon State University. The sequencer uses a modified Sanger method (Sanger 1977) with fluorescently
labeled dideoxy terminators. The region was bidirectionally sequenced to insure sequences of high quality.
PCR Amplification, Screening, and Scoring of ISSR Markers
The ISSR method utilizes a single primer anchored to a simple sequence repeat (SSR), or microsatellite,
to amplify the region between adjacent repeats (Gupta ci aL 1994). Amplifications were carried out in lox buffer (20mM Tris-C
2 H 3 0 2 , pH 9,0; 10 mMNH 4 SO 4 sulfate; 75mM KC 2 H 3 0 2 0.05% Tween® 20); 1.5 mMMgCl 2 5 mg/mL acetamide; 100MM dNTPs;
0.15 MM primer; SOng of genomie DNA; fluorescent dUTPs (0.2 mlvi [Ri 10]dUTP or 0.2mM [R6G]dUTP or 0.8 mM dUTP [TAMIRA]);
and 1.9 units of Taq polymerase (Promega) per 25 1 iL reaction. A PTC- 100 thermal cycler (MJ Research) was programmed for
34 cycles ofthe following parameters: denaturation for 30 s at 93°C,annealling for 30s at 45°C,and elongation for 2 mm at
72°C. A single 10 mm elongation at 72°Cfollowed the cycles, and reactions were held at 4°Cuntil removed from the thermal cycler.
Two methods of electrophoresis were utilized to visualize bands of different sizes. First, we loaded 15 giL of each reaction
into a 2% agarose gel stained with EtEr, and electrophoresed for 4 h in TBE (90 mM Tris base; 90 mlvi Boric Acid; 2 mM EDTA,
pH 8.0). A 1 Kb ladder (GibcoBRL) served as a molecular size standard. Amplification products between 1000 and 1400 bp were
scored as present (1) or absent (0) from the agarose gel. Further, PCR products were electrophoresed on an ABJ Prism 377 automated
sequencer (Applied Biosystems) at the Central Services Laboratory at Oregon State University. Briefly, amplification products
were loaded in a 5% LongRanger (6 M Urea) denaturing gel, and run for 5 h with Filter Wheel A, and Genescan 2500 [Rox] as
an internal lane standard. Gels were analyzed utilizing GeneScan Analysis 2.1 software (Applied Biosystems). Fragments between
400 and 1000 bp were scored as present (1) or absent (0) using Genotyper software (Applied Biosystems). Ninety primers from
the University of British Columbia Biotechnology Laboratory (IJBCBL) primer set #9 were initially screened for PCR amplification
on agarose gels. Of these, 44 with distinct, well separated bands were screened for variability using
4 individuals from across the North American range (Cabo Rojo, Puerto Rico; St. George Island, Florida; Lower Rio Grande Valley
NWR, Texas; San Diego County, California). We chose 16 bright, reproducible bands from eight primers for scoring in the final
analysis. The primers utilized in the final analysis were UBC826: (AC)8 T; UBCS46: (CA)8 RT; UBCS4S:(CA)3 RG; UBC849: (GT)8
YA; IJBC85O: (GT)8 YC; UBC857 (AC)8 YG; UBC859: (TG)3 RC; and UBC886: VDV(CT)7 . DNA was amplified from single extractions
to minimize variation in DNA quality and concentration. Negative controls were run to prevent scoring artifact bands. Low
quality DNA samples were eliminated from the analyses; thus, the set of samples analyzed does not completely overlap the mtDNA
data set, and the Colorado population is not represented in ISSR analyses.
Mitochondrial DNA Sequences Snowy Plover sequences were compared to control region sequences from Ruddy
Turnstone (Arenaria interpres) obtained from Genbank (Accession number L2OI 36) and Piping Plover (Charadrius melodus; S.
Haig & T. Mullins, unpublished data). Alignments were constructed by hand using the Genetic Data Environment (ODE) software
package (Smith et al. 1992). Estimates of nucleotide diversity (it) and haplotype diversity (h) were derived using the equations
of Nei (1987). Genetic distance between haplotypes was estimated using the Kimura 2-parameter method (Kimura 1980). We estimated
population subdivision via analysis of molecular variance using haplotype frequencies (AM0VA; Excoffier ci al. 1992). In this
analysis, FST, the correlation of random genotypes within a population relative to the species, is analogous to Wright's FST
statistic (Wright 1951). To test significance of the F statistics, values were compared to a null distribution obtained by
generating 1000 random permutations of the haplotypes among populations. To further test the null hypothesis that the distribution
of haplotypes among populations was random, an exact test of population differentiation was conducted (Raymond & Rousset
1995). All of the analyses listed above were carried out in the program ARLEQUIN 1.1 (Schneider et al. 1997).
ISSR markers have a dominant or codominant mode of inheritance (e.g., Gupta ci al. 1994, Wolfe et
al. 1998), and thus were scored as phenotypes of "present" or "absent". Absence of a band is assumed to indicate the loss
of a priming site through divergence, deletion, or chromosomal rearrangement (Wolfe & Liston 1998). Scored loci were assumed
to be non-allelic and in Hardy-Weinberg equilibrium. To correct for bias associated with dominance, band frequencies, percent
polymorphic loci (P, 95% criteria), and expected heterozygosity (He, Nei 1978) were calculated using Lynch and Milligan's
(1994) Taylor expansion in TFPGA (Miller 1998b).
In addition, the effective number of alleles per
locus (As) was calculated using the program POPGENE (Yeh ci al. 1997). Phylogenetic relationships among populations were examined
using the neighbor-joining method (Saitou & Nei 1987). Distance trees were constructed by calculating a matrix of Manhattan
distances (Prevosti distance in Wright 1978) of ISSR phenotypes in the program RAPDDIST (Black 1998b). One hundred bootstrap
matrices were calculated in RAPDDIST, and used to construct a consensus bootstrap tree with NEIGHBOR and CONSENSE programs
in PHYLIP v 3.5C (Felsenstein 1993). We used non-metric multidimensional scaling, a multivariate approach based on ranked
differences between populations, to qualitatively assess differences between populations (Kruskal 1964, Mather 1976). Distance
between populations was estimated using Jaccard's coefficient (Jaccard 1908 in Sneath & Sokal 1973) because of its utility
for two-state data (±). Final stress (opposite of goodness of fit) was plotted against number of dimensions to determine the
minimum number of dimensions that accurately assessed minimum stress; three dimensions were chosen for the final analysis.
A plot of number of iterations versus stress determined 30 iterations were sufficient to obtain a stable final solution, and
100 iterations were run in the final analysis. Monte Carlo simulations were conducted to estimate the probability that the
final stress could have been determined by chance. All calculations were carried out with PC-ORD software (McCune & Mefford
1999). Genetic differentiation among population groups suggested by geographic barriers were tested using analysis of molecular
variance (AMOVA; Excoffier ci al. 1992). First, a matrix of squared Euclidean distances between populations was prepared with
AMOVA-Prep software (Miller 1998a), where E = n [1 - (nxy/n)], and nxy is the number of bands shared by two individuals (Excoffier
ci al. 1992). We then calculated within and among group variance components and FST, in WINAMOVA (Exeoffier 1993). In addition, the level of heterogeneity within these groups was compared to that expected by chance using multi-response
permutation procedure (MRPP), carried out on the Euclidean distance matrix in PC-ORD (McCune & Mefford 1999). MRPP is
a non-parametric approach, which does not require multi variate normality and homogeneity of variances (Mielke 1984). Pairwise
comparisons of genetic differentiation among populations were made with an exact test of population differentiation (Raymond
& Rousset 1995), using 1000
permutations in TFPGA (Miller 1998b). In order to assess the relative
historical influence of gene flow and genetic drift on population structure, we compared the relationship between genetic
distance and geographic distance to that expected under a stepping-stone model of population structure (Kimura 1953, Malécot
1955, Kimura & Weiss 1964). Pairwise FST values were calculated as an approximation of genetic distance (Reynolds ct al.
1983), using a second order Taylor expansion (Lynch & Milligan 1994) in RAPDDIST (Black 1998). Following Hutchinson and
Templeton (1999), we constructed scatterplots of pairwise F 5 T values versus geographic distance to see if genetic distance
and scatter increased with geographic distance. Residuals from a linear regression of pairwise F 51 values versus geographic
distance were plotted against geographic distance. A Mantel test, which allows for lack of independence among points, was
used to examine correlation between the geographic distance matrix and matrices of pairwise FST values and residual values
(Hutchinson & Templeton 1999). Significance of the resulting correlation coefficients (r-values) was assessed with a Monte
Carlo simulation (1000 iterations) using PC-ORD software (McCune & Mefford 1999). Lynch and Milligan's (1994) F~1, statistic,
and the associated estimate of gene flow (Nm; Wright 1951), were calculated in the program RAPDFST (Black 1998a).
MIDNA Sequence Data
In the 322 bp of sequence obtained from 126 individuals, we observed 6 variable nucleotide sites comprising
8 haplotypes (Table 1). Of these 8 substitutions, 7 were transitions and 1 was a transversion. No haplotypes were separated
from other haplotypes by more than a single base pair change, as indicated in a hand-generated parsimony netivork (Figure
2). Haplotypes varied greatly in their degree of geographic dispersion (Figure 3). The most common haplotypes, C and D, were
quite widespread and were identified in all populations in the continental United States. These haplotypes accounted for 83.3%
of individuals. Haplotype C was also present in two individuals from Peru. In contrast, haplotype A was unique to Puerto Rico
and was the only haplotype represented in that population. Haplotypes B and G were present at low frequencies in the Texas
and the Great Plains populations. Haplotype H was also present in one individual from Kansas and one from Colorado. Haplotype
F was unique to California, but appeared only in a single individual. The final haplotype was identified in three individuals
from Peru and one from Eastern Oregon. Overall, nucleotide diversity for mtDNA sequences was low within all populations, averaging
0.003, and ranging from 0 to 0.003 (Table 2). Haplotype diversity, on the other hand, ranged from 0 in Puerto Rico and Colorado
populations to 0.7 14 in Texas, with an overall value of 0.617 ±0.031 (SE) for continental North America and Puerto Rico.
Haplotype diversity was lower in Florida than in western populations, which all reflected haplotype diversity values between
0.485 and 0.679 (Table 2).
Analysis of molecular variance (AMOVA) revealed that the majority of variance in haplotypic
diversity was attributable to the within population component, regardless of which populations were considered (Table 3).
When only continental populations were considered, analysis of molecular variance provided no evidence of structure (FST =
0.025, p = 0.178). However, when all North American populations were considered (including Puerto Rico), significant differentiation
between populations was observed, indicating that population structure was due to differentiation between Puerto Rico and
mainland populations (Table 3). Hierarchical analysis revealed that birds with different migratory patterns (east and west
of the Rocky Mountains) were not structured with respect to mtDNA haplotypes. The exact test of population differentiation
revealed low probability of non-differentiation between the Puerto Rican population and all continental populations (Table
4). Among eastern populations, Texas appears to be significantly differentiated from Florida. Texas is also significantly
differentiated from California and coastal Oregon. There was no support for differentiation of Florida from any population
west of the Great Plains, norevidence that the Great Plains population was differentiated from any other continental population
at mtDNA haplotypes. Furthermore, the exact test provided no support for population differentiation in the western United
States. Results indicate that the probability of non-differentiation between all pairs of western populations (coastal Oregon,
Eastern Oregon, California, and Utah) is greater than 0.45, and the probability of non-differentiation between Oregon and
California populations is 1.0.
Screening for genetic variability resulted in 16 scorable JSSR loci (Table 5), ranging in size from 450
to 1400 bp (Table 5). Across all populations, the effective number of alleles was 1.45 ± 0.30, the percentage of polymorphic
loci was 87.5, and level of expected heterozygosity was 0.28 (Table 6). Puerto Rico had the lowest genetic diversity of all
populations, while western populations had the highest diversity. Neighbor-joining analysis revealed a cluster of populations
east of the Rocky Mountains (Figure 4). In contrast to mitochondrial analysis, Puerto Rico was not isolated from continental
populations. In fact, Puerto Rico fell within a cluster of eastern continental populations, and grouped most closely with
Florida, although bootstrap support for the Florida - Puerto Rico relationship was weak. Populations from Texas and the Great
Plains clustered with Florida and Puerto Rico, which suggests that the nivosus subspecies may be aparaphyletic group. Non-metric
multidimensional scaling of 16 loci further supported differences between populations east and west of the Rocky Mountains.
Two dimensions (axes) represented most of the variance among individuals, while the third accounted for the remainder (axis
1: R 2 = 0.15, axis 2: R 2 = 0.26, axis 3: R 2 = 0.08, cumulative R 2 = 0.48). Additional dimensions did not significantly
improve the model. Monte Carlo simulations indicated that the probability of obtaining a lower final stress by chance was
very low (final stress = 15.95, p <0.01). A plot of the final solution reveals clustering of populations east and west
of the Rocky Mountains, with a good deal of overlap among individuals within these two groups (Figure 5). Of the populations
east of the Rocky Mountains, the Great Plains population displayed the greatest overlap with western populations.
Analysis of molecular variance (AMOVA) of ISSR markers provided
further support for differences between eastern and western groups. When nested analyses were conducted, a greater proportion
of the variation was attributable to groups (east and west of the Rocky Mountains) than among populations within groups, or
within groups (Table 3). The eastern group appeared to be more highly structured than the western group, with a greater proportion
of the variance attributable to the among population component. Likewise, multi-response permutation procedures resulted in
significant within group agreement among all populations (R = 0.33, p <0.05) and for the eastern and western groups (R
= 0.22; p <0.05). Differentiation among all pairs of populations was compared with an exact test of population differentiation
(Table 4). As with mtDNA results, Puerto Rico had a high probability of differentiation from all continental populations.
Results provide evidence that Florida is differentiated from the Great Plains; however, Texas was not significantly differentiated
from either Florida or the Great Plains. All populations in the eastern group (Puerto Rico, Florida, Texas, Great Plains)
appear to be significantly differentiated from the western group (California, coastal Oregon, and eastern Oregon, and Utah).
Members of the western group were not significantly differentiated from one another.
however, had a much lower probability of differentiation from other western populations. The correlation between pairwise
F 5 T values and geographic distance was examined to illustrate the relative influence of genetic drift and gene flow on population
structure. The null hypothesis of equilibrium can be rejected if the relationship between geographic distance and genetic
distance is not significant, or if a scatterplot of geographic distance does not reveal a positive relationship between genetic
and geographic distance and a pattern of increasing variance (Hutchinson & Templeton 1999). Comparison of pairwise F 5
T values and geographic distance reveals a positive linear relationship, starting near the origin, as predicted under an equilibrium
model (Figure 6). A Mantel test revealed a highly significant relationship between pairwise F 5 T and geographic distance
(r = 0.85; p <0.01). However, we did not observe a significant positive correlation between the residuals and geographic
distance (r = 0.00; p = 0.47), as expected under equilibrium conditions. As equilibrium conditions were not met, it was inappropriate
to estimate gene flow across all populations. Assuming equilibrium conditions within groups, the eastern group (F 5~ = 0.30
± 0.24; N,,, = 0.6) appears to be more highly structured and has lower levels of gene flow than the western group (F 5 T =
0.16 ± 0.13; N~= 1.3). However, if Puerto Rico is excluded from the eastern group, gene flow between populations is higher
east of the Rocky Mountains (FST = 0.06 ± 0.07; Nm = 3.9). Gene flow between the eastern and western groups appears to be
limited (F 5 T = 0.25 ± 0.22; N,,, = 0.7).
Genetic Variation within Populations
Few studies have estimated within population genetic variation in shorebirds (Baker & Strauch 1988,
Haig & Oring 1988a, Baker 1992, Haig et aT 1997). Allozyme analyses revealed low variability within populations of Piping
Plovers (Haig & Oring 1988a). Detailed phylogeographic analyses have been completed on only one shorebird species, the
Dunlin (Calidris alpina; Wenink et al. 1993, 1994, 1996; Wenink and Baker 1996). Interestingly, estimates of within population
haplotypic diversity in Dunlin for Region II of the mtDNA control region were similar to those of Snowy Plover populations
(Dunlin: h = 0.00 - 0.72). When data from cytochrome b and control Region I were considered, however, higher levels of within
population diversity were observed, ranging from 0.63 to 1.00 (Wenink 1994). This illustrates the need to exercise caution
when interpreting estimates of genetic diversity based on a single region, as different domains of the control region appear
to have different mutation rates (Baker & Marshall 1997). Overall, our estimates of haplotypic and nucleotide diversity
are low, but within the range observed among other avian species (Gill et al. 1993, Seutin a al. 1994, Gill & Slikas 1992).
Populations of Snowy Plover varied in degree of genetic diversity, with the Puerto Rico samples exhibiting the lowest level
at both mtDNA and ISSR markers. Given the small population size of the resident Puerto Rican population (n = 40; Lee 1989),
a local population bottleneck may have led to a reduction in genetic variability. Low variability may also be attributable
to founder effect, given the insularity of the Puerto Rican population. Comparison of the Puerto Rican population with other
populations in the Caribbean will be necessary in order to distinguish between these hypotheses. Within the continental United
States, patterns of genetic diversity at ISSR loci were not concordant with measures o fhaplotypic diversity. The Texas population
had the highest mtDNA haplotypic diversity, while western populations had the highest diversity at
Genetic Structure and Phylogeographic Patterns
MtDNA structure among North American and Caribbean Snowy Plovers is due primarily to differentiation of
Puerto Rican Snowy Plovers from continental populations. While ISSR markers reveal significant differentiation between Puerto
Rico and mainland populations, Puerto Rico appears to be paraphyletic with populations in the Eastern United States. Thus,
lack of overlap in haplotype composition between Puerto Rico and continental populations is consistent with a recent population
bottleneck, versus a prolonged period of evolutionary separation from mainland populations.
analyses support further subdivision of continental populations of Snowy Plovers in North America into two groups: east and
west of the Rocky Mountains. While we observed no statistical differences in mtDNA haplotype composition between groups, three
haplotypes common to populations in the central United States were not observed west of the Rocky Mountains. Groups east and
west of the Rocky Mountains are concordant with known migratory patterns: eastern birds have been observed to winter along
the coast of the Gulf of Mexico, while western populations winter along the Pacific coast and the Gulf of California (Page
a aT 1995b; W. Conway, pers. comm.).
Within the eastern United States, pairwise comparisons of
populations suggest some structure to populations. However, the results are not consistent between mtDNA and ISSR markers,
except that Texas and Great Plains populations do not appear to be differentiated. Qualitative examination of mtDNA haplotype
distribution reveals two haplotypes identified exclusively among Texas, Great Plains, and Colorado populations. Presence of
these shared haplotypes might reflect either a pattern of gene flow between resident Texas and migratory Great Plains birds
or a common population history. For example, the Great Plains breeding population could represent a historic breeding range
expansion of the Texas population. Theory on the evolution of migratory behavior in birds predicts that partial migration
is a stage in the evolution of fully migratory populations (e.g., Cox 1985).
Lack of differentiation
between coastal Oregon, California, and Eastern Oregon is not unexpected because birds have been observed dispersing between
breeding sites in the region (Stenzel a al. 1994). Within a single breeding season, female birds from the central California
coast were observed as far south as San Diego County (615 km) and inland at Summer and Abert Lakes, Oregon, 660 km north in
the Great Basin. Between seasons, movements as far north as Washington (1140 km) have been documented (Stenzel a al. 1994).
The low magnitude of structure observed among western populations might also reflect high rates of natal dispersal. Rates
of natal philopatry have not been examined for North American Snowy Plovers, but anecdotal evidence suggests that dispersal
of hundreds of kilometers from the natal site is not uncommon (Page et at. 1995 a).
At ISSR loci,
the Great Salt Lake population has a much lower probability of non-differentiation from birds at other Great Basin sites or
Pacific coastal sites, although the populations are not differentiated at mtDNA alleles. One might predict differentiation
between these populations because movements between the Great Salt Lake and other Great Basin sites have not been documented
despite banding efforts, and high breeding site fidelity at Great Salt Lake sites has been reported (Paton & Edwards 1996).
In addition, birds from the Great Salt Lake apparently do not mix as extensively on wintering grounds with coastal California
birds as birds from the Western Great Basin (Page a al. 1995b). While mixing of populations on the wintering grounds is not
as reliable an indicator of gene flow as natal philopatry and breeding site fidelity (Oring and Lank 1984), in this case,
the breeding range overlaps the wintering range. Because pair bond formation occurs on the breeding grounds and may occur
up to two months before egg laying (Warmer et al. 1986, Chase & Gore 1989), wintering together might increase the likelihood
of pair bond formation and decrease the likelihood of a bird to migrate.
Overall, mtDNA markers
revealed less structure within continental North America than ISSR makers. Lack of concordance between mONA and ISSR markers
might reflect inadequate time since separation of populations for mONA differentiation to become apparent (Neigel & Avise
1986). The relatively low variation observed consistently across all populations supports this hypothesis. An alternative
hypothesis might be that differences between mtDNA and nuclear genes reflect differential dispersal patterns between males
and females. For organisms with sex-biased dispersal patterns, maternally inherited mtDNA may not accurately represent population
structure for biparentally transmitted genes or the species as a whole (Avise a' al. 1992, Melnick & Hoelzer 1992, Bowen
a at 1992). Shorebirds, on average, have higher rates of female than male dispersal, although this varies depending on the
mating system (Oring et at 1983, Oring & Lank 1984, Gratto et al. 1985, Haig & Oring 1988b, Colwell & Oring 1988,
Reed & Oring 1993). The Snowy Plover mating system is sequentially polygamous and females have lower breeding site fidelity
than males (Page et al 1983, Warriner et al 1986, Stenzel et al 1994, Paton 1994).
female dispersal is greater than male dispersal, we would expect to see weak population structure in both mtDNA haplotypes
and autosomal genes, as female mediated gene flow homogenizes both. Therefore, differences between male and female dispersal
alone do not adequately explain the observed pattern in Snowy Plovers. Because the effective population size for mtDNA is
generally smaller than for nuclear genes, mtDNA markers are expected to sort more rapidly, and often reveal structure not
observed at nuclear genes (Neigel & Avise 1986). However, the effective population size for mitochondrial genes may not
always be smaller than autosomal genes if the assumption of random mating is not met (Chesser & Baker 1996, Piertney et
at!. 2000). In Snowy Plovers, deviations from the assumption of random mating, including sex-biased dispersal rates and a
polygamous mating system, may explain greater structure at ISSR loci.
Given the putative barriers
to gene flow and complex dispersal patterns of Snowy Plovers, it is not surprising that the null hypothesis of equilibrium
was rejected at ISSR loci. Snowy Plovers, as partial migrants, may not fit the assumptions of the steppingstone model that
dispersal is limited to the closest populations. Although we found a highly significant relationship between geographic distance
and genetic distance, the scatter of points did not reflect a pattern of increasing variance expected under equilibrium (Hutchinson
& Templeton 1999). From the observed pattern, it is difficult to distinguish between a hypothesis of isolation by distance
and a more complex structure where individual regions are influenced differentially by gene flow and genetic drift. MtDNA
evidence for the isolation of the Puerto Rican population and the likelihood of the Rocky Mountains as a barrier suggest that
the latter might be the case.
Lack of equilibrium between gene flow and genetic drift indicates
that it would not be appropriate to estimate gene flow across all populations (Hutchinson & Templeton 1999). If we assume
equilibrium within regional groups, comparison of relative estimates reveals greater gene flow among western populations than
among eastern populations only when Puerto Rico is included in the eastern group. Migration of a few individuals per generation
has been shown to be adequate to prevent populations from becoming fixed for different alleles (Wright 1951). Our estimates
suggest that gene flow is higher than this threshold for continental populations of Snowy Plovers.
Our results call for a reexamination of subspecies boundaries for Snowy Plovers. At mtDNA alleles, the Florida population appears to be much more closely related to other continental populations
of C. a. nivosus than to the Puerto Rican population of C. a. tenuirostris. Conover (1945) similarly found no morphological
differences between birds from Florida and California, but lacked specimens to draw conclusions about the identity of Caribbean
populations. At ISSR loci, the Puerto Rican population is nested within a group of populations from the eastern United States.
However, this group also does not correspond to previously defined subspecies, as Great Plains and Texas birds appear to be
more closely related to eastern than western populations.
This study provides evidence of overlap
in haplotype composition of the South American subspecies (C. a. occidentalis) and the Western Snowy Plover (C. a. nivosus).
Haplotype C, one of the more common North American haplotypes, was also observed in two of five individuals from Peru. In
addition, a single copy of haplotype E was detected in eastern Oregon, which was otherwise found only in the Peruvian population.
The observed pattern might reflect incomplete lineage sorting in populations that diverged relatively recently on an evolutionary
time scale. Additional sampling of C. a. occidentalis would be required to detect statistical differences in allele frequencies
between the two subspecies. Overall, an examination of the global population structure of Kentish Plovers would be prudent
before drawing conclusions about subspecies boundaries.
Conservation Unit Designation
Preservation of biological diversity below the species level has become increasingly recognized as a priority
in endangered species conservation. One approach to this problem has been the identification of intraspecific conservation
units, which rely to a varying degree on genetic criteria (e.g., evolutionary significant unit, management unit, distinct
population segment; Ryder 1986; Waples 1991, 1995; Moritz 1994, Moritz et at 1995). Based on differences at ISSR loci, populations
of Snowy Plovers east and west of the Rocky Mountains could be considered separate management units according to Moritz's
(1994) criteria, requiring differences in allele frequencies at nuclear loci. Because the populations do not exhibit reciprocal
monophyly at mtDNA alleles, these groups could not be considered separate evolutionary significant units (ESU; Moritz 1994;
Moritz et. al. 1995). However, both ISSR and mtNDA data suggest a separate management unit (MU) for the Puerto Rican population
may also be warranted. While the Puerto Rican population appears to be closely related to other eastern populations, Puerto
Rico has significant differences in mtDNA and nuclear gene frequencies from all mainland populations. Additional sampling
of other Caribbean populations will clarify whether a separate or overlapping management unit should be considered for the
eastern United States and Caribbean.
Within management units, distinguishing whether populations
are separate genetic and demographic units is important because threats to the viability of these breeding populations differ. Coastal and inland populations of Snowy Plovers in the western United States
are currently being managed separately; coastal populations are protected as a Distinct Population Segment under the U.S.
Endangered Species Act, while inland populations are not listed. Our study provides no evidence of genetic differentiation
between coastal and inland populations. Likewise,
decisions must be made about the status and listing of the resident Florida population (n = 400 individuals; Chase & Gore
1989), the resident coastal Texas population (n = 500 individuals; Page et al. 1995a), and the migratory Great Plains population
(n = 1900-2150 individuals; Page et al. 1995a).
Considered together, Texas and Great Plains populations
contain unique mtDNA lineages within North American Snowy Plovers which should be recognized to conserve genetic diversity.
Nonetheless, evidence for genetic differentiation among populations in the eastern continental United States is inadequate
to support designation of additional management units, according to Moritz's criteria.
policy decisions about whether to manage populations separately, genetic data need to be evaluated in light of demographic
considerations. Evidence of differentiation at mtDNA alleles provides a strong indication of demographic autonomy (Avise 1995).
However, in some cases levels of gene flow may be high enough to prevent genetic differentiation when populations are demographically
independent (Avise 1995, Taylor & Dizon 1996). In the context of conservation, a population experiencing declines may
have enough dispersal to homogenize populations genetically, but inadequate recruitment to prevent extirpation (Taylor &
Dizon 1996, 1999). If the policy goal is to halt a decline over an extensive part of a species range, a conservative approach
may be to manage populations in smaller units based on demographic considerations.
to the U.S. Fish and Wildlife Service policy, the Eastern United States, the Western United States, and the Puerto Rican populations
could all be considered Distinct Population Segments due to the genetic discontinuity between populations, corresponding to
the three management units described above (Federal Register 1996). However, additional groups may also be considered Distinct
Population Segments, as the policy allows for other criteria to be used to judge distinctiveness of populations (Federal Register
1996). We recommend that both genetic and demographic considerations
be weighed in making policy decisions regarding the status and listing of Snowy Plover population segments.
American Ornithologists Union (1957) Check-list of North American Birds. Lord Baltimore Press, Baltimore.
Avise JC (1995) Mitochondrial DNA polymorphism and a connection between genetics and demography of relevance to
conservation. Conservation Biology, 9, 686-690.
Avise JC, Nelson WS, Ankney CD (1992) Matriarchal
population genetic structure in an avian species with female natal philopatry. Evolution, 46, 1084-1096.
AJ (1992) Molecular genetics of Calidris with special reference to Knots. Wader Study Group Bulletin, 64 (suppl.), 29-35.
Baker AJ, Marshall BID (1997) Mitochondrial control region sequences as tools for understanding evolution. In:
Avian Molecular Evolution and Systematics (ed. Mindell D), pp. 5 1-82. Academic Press, San Diego, California.
Baker AJ, Strauch JG (1988) Genetic variation and differentiation in shorebirds. Acta Xix Congressus Internationalis
Ornithologic. Pp. 1639-1645.
Binford LC (1989) A distributional survey of the birds of the Mexican
state of Oaxaca.
Ornithological Monographs, 43.
Black WC (I 998a) RAPDFST
- FORTRAN programs for analysis of genetic relationships among individuals using RAPD-PCR markers. Colorado State University,
Fort Collins, Colorado (ftp:!! larnar.colostate.edu).
Black WC (1998b) RAPDDIST 1.0. Department
of Microbiology, Colorado State University, Fort Collins, Colorado. (ftp:f/lamar.colostate.edu!pub/wcb4!).
Bowen BW, Richardson II, Meylan AB, Margaritoulis D, Hopkins Murphy R, Avise JC (1992) Global population structure
and natural history of the green turtle (Chelonia mydas) in terms of matriarchal phylogeny. Evolution, 46, 865-881.
Castro G, Myers JP (1988). Snowy Plover (Charadrius alexandrinus) records from Panama. American Birds, 42, 374.
Chase C, Gore J (1989) Snowy PloverBreeding Distribution. Final Performance Report. Nongame Wildlife Division,
Florida Game and Freshwater Fish Commsssion.
Chesser RK, Baker Ri (1996) Effective sizes and dynamics
of uniparentally and diparentally inherited genes. Genetics, 144, 1225-1235.
Colwell MA, OringLW
(1988) Return rates of prairie shorebirds: sex and species differences. Wader Study Group Bulletin, 55, 21-24.
Conover B (1945) Notes on some North American shorebirds. Condor, 57, 211-214.
(1985) The evolution of avian migration systems between temperate and tropical regions of the New World. American Naturalist,
Excoffier L (1993) W]INAMOVA. Genetics and Biometry Laboratory, University of Geneva,
Carouge, Switzerland (http:l/anthropologie.unige.chfftpfcomplwin/amova).
Excoffier L, Smouse P,
Quattro J (1992) Analysis of molecular variance inferred from metric distances among DNA haplotypes: application to human
mitochondrial DNA restriction data. Genetics, 131, 479-491.
Federal Register (1992) Endangered
and threatened wildlife and plants: determination of
threatened status for the Pacific coast population
of the Western Snowy Plover.
Federal Register 58, 12864-12874.
Register (1996) Policy regarding the recognition of distinct vertebrate population segments under the Endangered Species Act.
Federal Register61, 4722.
Felsenstein J (1993) PHYLIP (Phylogeny Inference Package) Version 3.5c.
Department of Genetics, University of Washington, Seattle, Washington.
Gi1IFB, Slikas B (1992) Patterns of mitochondrial DNA divergence in North American crested titmice. Condor 94,
Gill FB, Mostrom AM, Mack AL(1993) Speciation of North American chickadees. I. Patterns
of mtDNA genetic divergence. Evolution, 47, 195-212.
Gore J (1996) Cuban Snowy Plover. In: Rare
and Endangered Biota of Florida. V Birds (eds.Rodgers H, Kale W, Smith KT), pp. 6 1-72. University Press of Florida, Gainesville,
Gratto CL, Morrison RIG, Cooke F (1985) Philopatry, site tenacity, and mate fidelity in
the Semipalmated Sandpiper. Auk, 102, 16-24.
Gupta M, Chyi Y-S, Romero-Severson I, Owen IL (1994)
Amplification of DNA markers from evolutionary diverse genomes using single primers of simple sequence repeats. Theoretical
and Applied Genetics, 89, 998-1006.
Haig SM, Oring LW (1988a) Genetic differentiation of Piping
Plovers across North America. Auk, 105, 260-267.
Haig SM, Oring LW (1988b) Distribution and dispersal
in the Piping Plover. Auk, 105, 630-638.
Haig SM, Gratto-Trevor CL, Mullins TD, CoIwell, MA (1997)
Population identification of western hemisphere shorebirds throughout the annual cycle. Molecular Ecology, 6, 413-427.
Hayman P, Marchant J, Prater T (1986) Shorebirds: an Identification Guide. Houghton Mifflin, Boston.
Hill LA, Talent LO (1990) Effects of capture, handling, banding, and radio-marking on breeding Least Terns and
Snowy Plovers. Journal ofField Ornithology, 61, 310-319.
Hutchinson DW, Templeton AR (1999) Correlation
of pairwise genetic and geographic distance measures: inferring the relative influences of gene flow and drift on distribution
of genetic variability. Evolution, 53, 1898-1914.
Jaccard P (1908) Nouvelles recherches sur la
distribution florale. Bulletin Societe Vaud. Sciences Naturates, 44, 223-270.
Kimura M (1953) Stepping
stone model of population. Animal Report of the National Institute of Genetics, 3, 62-63.
M (1980) A simple method for estimating evolutionary rate of base substitution through comparative studies of nucleotide sequences.
Journal of Molecular Evolution, 16, 111-120.
Kimura M, Weiss GH (1964) The stepping stone model
of population structure and the decrease of genetic correlation with distance. Genetics, 49, 561-576.
JB (1964) Non-metric multidimensional scaling: a numerical method. Psychoenetrika, 29, 1-27.
G (1989) Breeding ecology and habitat use patterns of Snowy and Wilson's Plovers at the Cabo Rojo Salt Flats, Puerto Rico.
Master's Thesis, Clemson University, Clemson, South Carolina.
Lynch M, Milligan BG (1994) Analysis
of population genetic structure with RAPD markers. Molecular Ecology, 3, 91-99.
Malécot 0 (1955)
The decrease of relationship with distance. Cold Spring Harbor Symposium of Quantitative Biology, 20, 52-53.
Mather PM (1976) computational Methods of Multivariate Analysis in Physical Geography. J. Wiley and Sons, London.
McCune B, Mefford MI (1999) PC-ORD-Multivariate analysis of ecological data. Version 4,28 beta. MjM Software, Oleneden
Melnick DJ, HoelzerGA (1992) Differences in male and female dispersal lead to contrasting
distributions of nuclear and mitochondrial variation. International Journal of Primatology, 13, 379-393.
PW Jr (1984) Meteorological applications of permutation techniques based on distance functions. In: Handbook of Statistics,
Vol.4 (eds., Krishnaiah PR and Sen PK), pp. 813-830. Elsevier Science Publishers.
Miller MP (1998a)
AMOVA-PREP. Department of Biological Sciences, Northern Arizona University, Flagstaff, Arizona. (http://herb.bio.nau.edul-miller/amovaprp.htm).
Miller MP (1998b) Tools for population genetic analysis (TFPGA). Department of Biological Sciences, Northern Arizona
University, Flagstaff, Arizona. (http:tlherb.bio.nau.edul-millerI).
Moritz C (1994) Defining "Evolutionary
Significant Units" for conservation. Trends in Ecology and Evolution, 9, 373-375.
Moritz C, Lavery
5, Slade R (1995) Using allele frequency and phylogeny to define units for conservation and management. American Fisheries
Society Symposium, 17, 249-262.
Nei M (1978) Estimation of average heterozygosities and genetic
distance from a small number of individuals. Genetics, 41, 583-590.
Nei M (1987) MolecularEvolutionary
Genetics. Columbia University Press, New York.
Neigel SE, Avise JC (1986) Phylogenetic relationships
ofmitochondrial DNA under various demographic models of speciation. In: Evolutionary Processes and Theory (eds., Nevo E and
Karlin 5), pp. 5 15-534. Academic Press, New York.
Oring LW, Lank DB (1984). Breeding area fidelity,
natal philopatry, and the social systems of sandpipers. In: Shorebirds: Breeding Behavior and Populations (eds. Burger J,
Olla BL), pp. 125-147. Plenum, New York.
Oring LW, Lank DB, Maxson SJ (1983) Population studies
of the polyandrous Spotted Sandpiper. Auk, 100, 272-285.
Oring LW, Able KP, Anderson DW, Baptista
LF, Barlow IC, Gaunt AS, Gill AS, Gill FB, Wingfield JC, (1988) Guidelines for the use of wild birds in research. Au/c, 105
Page OW, Winkler DW, Swarth CW (1983) Spacing out at Mono Lake: breeding
success, nest density, and predation in the Snowy Plover. Auk, 100, 13-24.
Page OW, WarrinerJS,
WarrinerlC, Paton PWC (1995a) The Snowy Plover (Charadrius alexandrinus). In: The Birds of North America (eds. Poole A, Gill
F), no. 154. The Academy of Natural Sciences, Philadelphia, and the American Ornithologists Union, Washington D.C.
Page GW, Stem MA, Paton PW (1995b) Differences in wintering areas of Snowy Plovers from inland breeding sites in
western North America. Condor, 97, 258-262.
Paton PW (1994) Survival estimates for Snowy Plovers
breeding at Great Salt Lake, Utah. Condor, 96, 1106-1109.
Paton PWC, Edwards TC (1996) Factors
affecting interannual movements of Snowy Plovers. Au/c, 113, 534-543.
Piertney SB, MacCoil ADC,
Bacon PJ, Racey PA, Lambin X, Dallas JF (2000) Matrilineal genetic structure and female-mediated gene flow in Red Grouse (Lagopus
scoticus): and analysis using mitochondrial DNA. Evolution, 54, 279-289.
Raymond M, Rousset F (1995)
An exact test for population differentiation. Evolution, 49, 1280-1283.
Reed SM, Oring LW (1993)
Philopatry, site fidelity, dispersal, and survival of Spotted Sandpipers. Auk, 110, 541-551.
J, Weir ES, Cockerham CC (1983) Estimation for the coancestry coefficient: basis for short-term genetic distance. Genetics,
Ryder OA (1986) Species conservation and systematics: the dilemma of subspecis. Trends
in Ecology and Evolution, 1, 9-10.
Sanger F, Nicklen 5, Coulson AR (1977) DNA sequencing with chain
terminating inhibitors. Proceedings of the National Academy of Sciences USA, 74, 5463-5467.
N, Nei M (1987) The neighbor-joining method: a new method for reconstructing phylogenetic trees. MolecularBiology and Evolution,
Schneider 5, Kueffer J-M, Roessli D, Excoffier L (1997) Arlequin V. 1.1. A software
for population genetic analysis. (http://anthropologie.unige.chlarlequinl).
Seutin 0, Klein NK,
Ricklefs RE, Bermingham E (1994) Historical biogeography of the bananaquit (Coerebaflaveola) in the Caribbean region: a mtDNA
assessment. Evolution, 48, 1041-1061.
Smith SW, Wang C, Gillevet PM, Gilbert W (1992) Genetic Data
Environment and Harvard Genome Database. Genome Mapping and Sequencing. Cold Spring Harbor Laboratory. (http://fastlink.nih.gov/gde_sw.html).
Sneath PH, Sokal RR (1973) Numerical Taxonomy: The Principles and Practices of Numerical Taxonomy. W.H. Freeman
and Co., San Francisco.
Sprandel GL, Gore JA, Cobb DT (1997) Winter Shorebird Survey. Final Performance
Report. Florida Game and Freshwater Fish Commission, Talahassee, Florida.
Stenzel L, WarrinerJC,
WarrinerlS, Wilson KS, Bidstrup FC, Page OW (1994) Long-distance breeding dispersal of Snowy Plovers in western North America.
Journal of Animal Ecology, 63, 887-902.
Taylor BL, Dizon AE (1996) The need to estimate power to
link genetics and demography for conservation. Conservation Biology, 10, 66 1-664.
Taylor BL, Dizon
AE (1999) First policy then science: why a managment unit based solely on genetic criteria cannot work. Molecular Ecology,
Waples RS (1991) Pacific Salmon, Onchorynchus spp., and the definition of"Species"
under the Endangered Species Act. Marine Fisheries Review, 53, 11-22. Waples RS (1995) Evolutionary significant units and
the conservation of biological diversity under the U.S. Endangered Species Act. American Fisheries Symposium, 17, 8-27.
Warriner IS, Warriner JC, Page OW, Stenzel LW (1986) Mating system and reproductive success of a small population
of polygamous Snowy Plovers. Wilson Bulletin, 98, 15- 37.
Wenink PW, Baker AJ (1996) Mitochondrial
DNA lineages in composite flocks of migratory and wintering Dunlins (Calidris alpina) Auk, 113, 744-756.
PW, Baker AJ, Tilanus MGI (1993) Hypervariable control region sequences reveal global population structuring in a long-distance
migrant shorebird, the Dunlin (Calidris alpina). Proceedings of the National Academy of Science, USA, 90, 94-98.
Wenink PW, Baker AJ, Tilanus MGI (1994) Mitochondrial control-region sequences in two shorebird species, the Turnstone
and the Dunlin, and their utility in population genetic studies. Molecular Biology and Evolution, 11,22-31.
Wenink PW, Baker AJ, Rosne H-U, Tilanus MGJ (1996) Global mitochondrial DNA phylogeny of Holarctic breeding Dunlins
(Calidris alpina). Evolution, 50, 3 18-330.
Whitlock MC, McCauley DE (1999) Indirect measures of
gene flow and migration: FST=1/(4Nm+l). Heredity, 82, 117-125.
Wolfe AD, Liston A (1998) Contributions
of PCR-based methods to plant systematics and evolutionary biology. In: Molecular Systematics of Plants II: DNA Sequencing
(eds Soltis DE, Soltis PS, Doyle JJ), pp.43-So. Chapman Hall, New York.
Wolfe AD, Xiang Q-Y, Kephart
SR (1998) Assessing hybridization in natural populations of Penstemon (Scrophulariaceae) using hyervariable intersimple sequence
repeat bands. Molecular Ecology, 7, 1107-1125.
Wright S (1951) The genetical structure of populations.
Annals of Eugenetics, 15,323- 354.
Wright S (1978) Evolution and the Genetics of Populations. VoL
IL Variability Within and Among Natural Populations. University of Chicago Press, Chicago.
FC, Yang R-C, Boyle T, Ye Z-H, Mao JX (1997) POPGENE, the user friendly shareware for population genetic analysis. Molecular
Biology and Biotechnology Centre, University of Alberta, Canada. (http://www.ualberta.cal-fyeh).