Journal of Heredity Advance Access originally published online on February 24, 2005
Journal of Heredity 2005 96(4):381-387; doi:10.1093/jhered/esi040
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Genetic Differentiation Among Populations of the Salt Marsh Beetle Pogonus littoralis (Coleoptera: Carabidae): A Comparison Between Atlantic and Mediterranean Populations
From the Department of Entomology, Royal Belgian Institute of Natural Sciences, Vautierstreet 29, B-1000 Brussels, Belgium (Dhuyvetter, Gaublomme, Verdyck, and Desender); and University of Antwerp, Groenenborgerlaan 171, 2020 Antwerp, Belgium (Verdyck)
Address correspondence to Hilde Dhuyvetter at the address above, or e-mail: hildedhuyvetter{at}yahoo.com.
| Abstract |
|---|
|
|
|---|
Genetic differentiation was studied among populations of the salt marsh beetle Pogonus littoralis (Coleoptera: Carabidae), comparing Atlantic and Mediterranean populations over a large part of its range. The genetic structure of this highly mobile beetle was investigated by studying allozyme polymorphism at nine enzyme loci in 13 populations. Mediterranean P. littoralis were highly significantly differentiated from Atlantic populations. Moreover, more isolated Atlantic populations showed increased differentiation and decreased genetic diversity compared to less fragmented Mediterranean populations.
The genetic structure of a species is influenced by its past history and by current gene flow. In Europe, the present-day patterns of phenotypic and genetic variation of many organisms have their origin in Pleistocene events (Hewitt 1999; Taberlet et al. 1998). After the retreat of the ice sheets since the last glacial maximum (about 18,000 years ago), northern Europe was subsequently reinvaded from several refuges in southern Europe. Under this scenario, populations in southern Europe should be genetically more diverse than populations that recently recolonized northern Europe, thus producing south-north clines of genetic variability (Eber and Brandl 1997). In theory, latitudinal differences in variability can be expected, as the colonizing populations are more prone to founder events and serial bottlenecking stemming from the greater demographic instability of these marginal populations. Moreover, not only past history, but also the current balance between gene flow, genetic drift, and selection shapes the observed genetic patterns of species. In small populations at the edge of a species range, genetic drift may lead to reduced within-population genetic diversity and increased between-population differentiation, while populations admixture can lead to increased genetic diversity (Walter and Epperson 2001).
Another important factor in recent history is anthropogenic habitat fragmentation. Many areas of formerly continuous natural habitat have been subdivided into relatively small habitat islands surrounded by human-altered environments. Fragmentation generally results in reduced genetic diversity (Avise and Hamrick 1996; Frankham et al. 2002), but may also increase genetic differentiation between populations as a result of reduced gene flow (Slatkin 1994).
An excellent example of fragmented habitats are salt marshes. European salt marshes are relatively recent habitats, in many cases well documented historically. Present coastlines have only been in their current position for at most the last 6000 years (Adam 1990; Houthuys et al. 1993). Since their origin, the area of salt marshes in Europe has been considerably reduced due to geomorphological evolution and, in recent centuries, human activity (Schubel and Hirschberg 1978). Estimates all over western Europe show a recent and dramatic decrease in the surface area and habitat quality of salt marshes (Dijkema et al. 1984).
To investigate if an observed genetic differentiation results from fragmentation of the habitat or merely from a distance effect, the genetic population structure should be compared between continuous and more fragmented areas (e.g., Gibbs 1998; Knutsen et al. 2000; Peacock and Smith 1997; Van Dongen et al. 1998). These comparisons allow assessment of the impact of spatial habitat structure on population dynamics within species with specific life-history characteristics. However, these studies typically focus on species thought to have relatively low levels of vagility, possibly because fragmentation is more likely to disrupt gene flow in those species. Williams et al. (2003) were the first to identify genetic patterns consistent with recent habitat fragmentation in a wide-ranging, high gene flow butterfly on a large geographical scale.
Ground beetles (Coleoptera, Carabidae) appear to be ideal model organisms for such studies. The carabid beetle of this study, Pogonus littoralis, is a highly specialized halobiontic species, limited in its occurrence to salt marsh microhabitats, where it can reach relatively high densities (Desender 1989). This beetle has to be highly mobile because it is regularly forced to move between temporarily dry salt marsh ponds or creeks during its life cycle.
In this study we analyze the genetic diversity and differentiation of P. littoralis from the north to the south of its geographical distribution area using allozyme electrophoresis. The three main goals were (1) to investigate the distribution of genetic variation both within (genetic diversity) and among (genetic differentiation) Atlantic and Mediterranean populations; (2) to compare genetic diversity and differentiation among populations in continuous and more fragmented areas on a comparable geographical scale; and (3) to investigate genetic differentiation and variability in populations from two regions.
| Materials and Methods |
|---|
|
|
|---|
Studied Species and Sampling
Pogonus littoralis is a halobiont and highly mobile ground beetle, with a specialized microhabitat preference for unvegetated, temporarily dry salt marsh ponds or creeks, where it lives between cracks in humid sea clay. The geographical area of the species extends from the Mediterranean to western European Atlantic coast salt marshes, north to the British Isles and The Netherlands (Turin 2000). This constantly macropterous species shows maximally developed wings and functional flight musculature and is active during the day (Desender 1989).
Since 1993 we have collected P. littoralis beetles from five Atlantic (region with discontinuous or more fragmented distribution) and seven Mediterranean salt marshes (region with more continuous distribution) varying in size and isolation in France. We also sampled one population in Thessaloniki (Greece). All study sites are shown in Figure 1. In some areas, different populations could be studied. Beetles were collected by standardized hand-catches (unit of effort), mostly by using an aspirator, or by flotation whenever brackish water was available near a sampling site. Beetles were transported alive, and identified and counted with a binocular dissecting microscope. They were then frozen in liquid nitrogen or in an ultrafreezer until subjected to electrophoresis.
|
|
Allozyme Electrophoresis
Individuals were prepared for allozyme electrophoresis by homogenizing part of the body (head and thorax) in 40 µl of distilled water on ice. Overall, more than 600 beetles, at least 40 per population (Table 1) were processed for this study. Continuous cellulose acetate electrophoresis was carried out using standard methods (Hebert and Beaton 1989). Two buffer systems were used: Tris-glycine 10% (pH 8.5; Hebert and Beaton 1989) and Tris-maleate (pH 7.8; Richardson et al. 1986). Genotypes were scored for nine putative loci coding for seven enzymes: aldehyde oxidase (AO; EC 1.2.3.1 [EC] ), peptidase-A (dipeptide substrate: leucyl glycine, PEP-A1 and PEP-A2; EC 3.4..), peptidase-Z (tripeptide substrate: leucyl glycine glycine, Pep-Z; EC 3.4..), and phosphoglucomutase (PGM; EC 2.7.5.1 [EC] ) on a Tris-glycine buffer, and glucose-6-phosphate isomerase (GPI; EC 5.3.1.9 [EC] ), malate dehydrogenase (MDH; EC 1.1.1.37 [EC] ), and isocitrate dehydrogenase 1 and 2 (IDH1, IDH2; EC 1.1.1.42 [EC] ) on a Tris-maleate buffer.
|
Data Analysis
Tests for Hardy-Weinberg equilibrium and linkage disequilibrium were performed with GENEPOP (Raymond and Rousset 1995). Genetic diversity within populations was estimated using POPGENE (version 1.31; Yeh et al. 1997) in three ways: (1) the mean number of alleles over all loci, (2) the percentage of polymorphic loci, and (3) gene diversity: Nei's (1978) unbiased expected heterozygosity (He), averaged across loci. The mean number of alleles was standardized to samples of 40 individuals using the rarefaction technique (Hurlbert 1971). Differences between means were tested using STATISTICA (StatSoft Inc., Tulsa, OK). Genetic structure and differentiation were analyzed by evaluating among-population and among-region differentiation. Total genetic variance was therefore partitioned among regions, among populations within regions, and within populations by carrying out a hierarchical analysis of molecular variance (AMOVA) using ARLEQUIN (version 2.000; Schneider et al. 2000). GENETIX (Belkhir et al. 19961998) yielded FST estimates (Weir and Cockerham 1984), tested for their significant departure from the null hypothesis (no differentiation) by means of a permutation procedure. Wright's FST was also calculated (Wright 1969). Genetic distances between populations were visualized in an unweighted pair group method with arithmetic mean (UPGMA) dendrogram based on Nei's (1978) unbiased genetic distance.
| Results |
|---|
|
|
|---|
Genetic Variability
A total of 33 alleles were scored at nine loci across 13 populations. The number of alleles detected at each locus ranged from two (IDH1 and PEP-Z) to six (AO). Allele frequency tables for all loci are available upon request from the authors. No significant deviations from Hardy-Weinberg equilibrium were observed for any of the populations or loci. No significant linkage disequilibrium was found after applying the sequential Bonferroni correction. Therefore the studied loci can be interpreted as independent markers.
Within-Population and Within-Region Genetic Diversity
Genetic diversity varied a lot within regions, but there was a clear tendency for higher mean values in the Mediterranean area (Table 1): a Mann-Whitney U-test showed significantly higher values of mean expected heterozygosity (U = 0.0, P = .004) as well as mean alleles per locus (U = 0.0, P = .004) and mean percentage of polymorphic loci (U = 0.0, P = .004) for Mediterranean populations than for Atlantic ones. This higher variability is only to some degree caused by private alleles (three unique Atlantic alleles as compared to six for the Mediterranean; all in very low frequencies), but results especially from several nonunique alleles occurring at more equal frequencies and thus increasing heterozygosity in Mediterranean populations. Inspection of the allele frequencies shows that a higher diversity in the Mediterranean samples is especially visible at the PEP-A1 locus (cf. Figure 1). The population from Greece has a much lower overall genetic diversity compared to the other Mediterranean samples from France.
Among-Population and Among-Region Genetic Differentiation
The overall FST value among all populations was 0.165 (P < .0001; Weir and Cockerham 1984). The Atlantic region FST was 0.023 (P < .0030) and the Mediterranean FST was 0.013 (P < .0001; sample of Greece not included). We also computed Wright's FST (Wright 1969) and obtained an overall FST of 0.161. Genetic differentiation estimates clearly increase at a larger geographic scale.
Based on Nei's (1978) genetic distances, we constructed a UPGMA dendrogram (Figure 2). All samples from Mediterranean France cluster together, as well as those from Atlantic France, whereas the sample from Thessaloniki forms a third group.
The major part of genetic variation was found within populations (75.36%), with 23.55% among regions and only 1.10% among populations within regions (AMOVA; Table 2). Exact tests showed a strong genetic differentiation among the three regions (all three comparisons P < .001). To test whether the Mediterranean France and Atlantic France groups really represent two differentiated regions, we excluded the Greek population from the AMOVA. Variance among regions dropped to 18.74%, but was still highly significant (Table 2). When the two regions (Mediterranean France and Atlantic France) were analyzed separately, most of the variation was again detected within populations, with up to 98.86% of the variance found within populations in Mediterranean France, but variance among populations was still highly significant (Table 2).
|
| Discussion |
|---|
|
|
|---|
Allozyme Variability in Perspective
Significant genetic structuring (allozymes) has been reported for many insects, including beetles (Assmann and Weber 1997; De Jong et al. 2001; Desender et al. 1998; Desender and Serrano 1999; Desender and Verdyck 2001; Dhuyvetter et al. 2004; Epps et al. 2000; Hsiao 1989; Knoll and Rowell-Rahier 1998; Knoll et al. 1996; Knutzen et al. 2000). Estimates of overall genetic differentiation in the present study (FST = 0.161) are higher than the mean values (Wright's FST = 0.103) obtained for 30 other beetle species (Hsiao 1989), which are known to be among the highest recorded for insects (Ward et al. 1992). Next, FST values of P. littoralis are compared to those obtained in other salt marsh beetles. We only included comparable studies using both Atlantic and Mediterranean population samples to calculate an overall FST. Overall FST in the wing polymorphic Pogonus chalceus (see Desender and Serrano 1999; Dhuyvetter et al. 2004) appears to be larger, as expected, than in the related highly mobile P. littoralis. In P. chalceus, nearly 25% of the total observed genetic variation is due to differentiation among populations (Wright's FST; Desender and Serrano 1999). Theories that relate variation in FST to variation in rates of gene flow indeed predict that species with a high dispersal power should show less population structuring (Waples 1998; Ward et al. 1992). In Bembidion normannum, another highly mobile salt marsh beetle, the overall FST is much lower than in P. littoralis (Wright's FST = 0.059; Desender and Verdyck 2001). The higher FST in P. littoralis is possibly due to the specific biology and microhabitat preference of this species (see further).
Population Differentiation
Mediterranean P. littoralis are highly significantly differentiated from Atlantic populations (Table 2 and Figure 2). Atlantic beetles of two other salt marsh species, P. chalceus and B. normannum, were also genetically distinct from Mediterranean populations (Desender and Serrano 1999; Desender and Verdyck 2001; Dhuyvetter et al. 2004). The highest gene diversity estimates of P. littoralis were found in the Mediterranean region and the lowest values were found in the Atlantic region (Table 1). This was also the case for B. normannum (Desender and Verdyck 2001), and is consistent with a much higher incidence of P. littoralis in Mediterranean salt marshes. The high genetic diversity observed within Mediterranean populations gives a strong indication that the evolutionary origin of P. littoralis lies in the Mediterranean area and/or that this area has served as glacial refugium, which was also suggested for B. normannum (Desender and Verdyck 2001).
We found that P. littoralis populations were genetically differentiated in both areas. However, the level of genetic differentiation among populations, as estimated by FST, was considerably greater in the Atlantic area, where the species shows a much more discontinuous distribution than in the Mediterranean area, where P. littoralis has a more continuous distribution. A similar result was obtained in another study comparing a fragmented area with a more continuous one in the highly mobile butterfly Speyeria idalia (Williams et al. 2003).
In the Mediterranean populations, we also observed a high level of genetic diversity within populations. This suggests that all populations are continuously exposed to gene flow. In view of the microhabitat preference of P. littoralis, populations in the Mediterranean area are probably to be viewed as metapopulations made up of ephemeral local populations. The resulting regular emigration (or local extinction) versus colonization dynamics is expected to enhance gene flow and consequently to continuously break down genetic differentiation, at least at a smaller geographic scale. In this area, the lowest within-population diversity was observed in the sample from Thessaloniki (Greece). This could be the result of genetic drift (random fixation of alleles) in less dense populations at the eastern edge of the species distribution. The population from Thessaloniki showed a three to eight times lower relative density compared to other Mediterranean samples. In metapopulation systems, edge populations have fewer neighbors than core populations. Fewer neighbors mean fewer immigrants, reduced chances of demographic rescue, and less gene flow (Segelbacher and Storch 2002). To understand this situation better in P. littoralis, we need to study some populations between France and Greece; for example, from Italy.
For the Atlantic populations, we suggest that P. littoralis has lost genetic variation due to founder events, population bottlenecks, or both. This interpretation is supported by the finding that populations sampled in southern Europe are more variable than northern European populations, a finding consistent with the general tendency for populations in areas glaciated during the Pleistocene to exhibit reduced genetic variability (e.g., Merilä et al. 1996). Five of the nine studied enzymes are fixed or almost fixed in Atlantic populations compared to Mediterranean ones. Moreover, the observed elevated genetic differentiation is probably a direct consequence of genetic drift in conjunction with reduced or unbalanced dispersal following habitat fragmentation. Latitudinal differences in variability could be expected as the colonizing populations are prone to founder events and serial bottlenecking stemming from the greater demographic instability of these marginal populations. This is suggested by the allele frequencies of the PEP-A1 locus, where allele two is fixed in the northernmost population from the Authie estuary (Figure 1). We intend to study these patterns and suggested processes further in the near future, if possible by means of additional and more powerful genetic markers, such as recently developed microsatellite markers (Dhuyvetter and Desender 2003) or DNA sequencing.
Although the observed elevated genetic differentiation in the Atlantic region might in itself have limited consequences, the genome-wide implications are reduced population sizes and loss of possible advantageous alleles, as well as fixation of disadvantageous ones. This could therefore ultimately increase the chances of extinction, even for a species with high dispersal power. Distribution data for P. littoralis in several western European countries indicate a significant regression of this salt marsh beetle during the past century (e.g., Desender et al. 1995; Hyman and Parsons 1992), as would be expected from the genetic data. It is, however, not yet clear whether habitat loss and fragmentation of Atlantic salt marshes or an overall decrease in their habitat quality is most responsible for this regression. In this context, we therefore plan to study some additional populations with more powerful markers.
| Acknowledgments |
|---|
This work was supported by the Entomology Department of the Royal Belgian Institute of Natural Sciences and the Federale Diensten voor Wetenschappelijke, Technische en Culturele Aangelegenheden (DWTC; project MO 36/006). It is carried out within the framework of the Flemish Research Network (FWO.010.97N, "Ecological genetics: patterns and processes of genetic variation in natural populations"; and FWO.017.02N, "Ecological genetics: a new approach"). A. Drumont assisted with part of the laboratory work. V. Versteirt and J.-Y. Rasplus helped in field sampling. T. Backeljau provided many helpful comments on the manuscript.
| Footnotes |
|---|
Corresponding Editor: Ron C. Woodruff
Received April 22, 2004
Accepted November 6, 2004
| References |
|---|
|
|
|---|
-
Adam P, 1990. Saltmarsh ecology Cambridge: Cambridge University Press.
Assmann T and Weber F, 1997. On the allozyme differentiation of Carabus punctatoauratus Germar (Coleoptera, Carabidae). J Zool Syst Evol Res 35:3343.
Avise JC and Hamrick JL, 1996. Conservation genetics New York: Chapman & Hall.
Belkhir K, Borsa P, Goudet J, Chikhi L, and Bonhomme F, 19961998. GENETIX, logiciel sous WindowsTM pour la génétique des populations Laboratoire Génome et Populations, CNRS UPR 9060. Montpellier, France: Université de Montpellier II.
De Jong PW, De Vos H, and Nielsen K, 2001. Demic structure and its relation with the distribution of an adaptive trait in Danish flea beetles. Mol Ecol 10:13121332.
Desender K, 1989. Dispersievermogen en ecologie van loopkevers (Coleoptera, Carabidae) in België: een evolutionaire benadering Studiedocumenten van het K.B.I.N.
Desender K, Backeljau T, Delahaye K, and De Meester L, 1998. Age and size of European saltmarshes and the population genetic consequences for ground beetles. Oecologia 114:503513.[CrossRef]
Desender K, Maes D, Maelfait J-P, and Van Kerckvoorde M, 1995. Een gedocumenteerde Rode lijst van de zandloopkevers en loopkevers van Vlaanderen Mededelingen van het Instituut voor Natuurbehoud.
Desender K and Serrano J, 1999. A genetic comparison of Atlantic and Mediterranean populations of a saltmarsh beetle. Belg J Zool 129:8394.
Desender K and Verdyck P, 2001. Geographic scaling and genetic differentiation in two highly mobile European saltmarsh beetles. Belg J Zool 131:3142.
Dhuyvetter H and Desender K, 2003. Isolation and characterisation of microsatellite loci in the saltmarsh beetle Pogonus chalceus (Coleoptera, Carabidae). Mol Ecol Notes 3:460461.[CrossRef]
Dhuyvetter H, Gaublomme E, and Desender K, 2004. Genetic differentiation and local adaptation in the saltmarsh beetle Pogonus chalceus: a comparison between allozyme and microsatellite loci. Mol Ecol 13:10651074.[Medline]
Dijkema KS, Beeftink WG, Doody JP, Gehu JM, Heydemann B, and Rivas Martinez S, 1984. La végétation halophile en Europe (prés salés) Strasbourg: Council of Europe.
Eber S and Brandl R, 1997. Genetic differentiation of the tephritid fly Urophora cardui in Europe as evidence for its biogeographical history. Mol Ecol 6:651660.[CrossRef]
Epps TM, Britten HB, and Rust RW, 2000. Allozyme variability, isolation, and dispersal of Eusattus muricatus (Coleoptera: Tenebrionidae) within silver state dune complex great basin, western North America. West North Am Nat 60:281290.
Frankham R, Ballou JD, and Briscoe DA, 2002. Introduction to conservation genetics Cambridge: Cambridge University Press.
Gibbs JP, 1998. Genetic structure of redback salamander Plethodon cinereus populations in continuous and fragmented forests. Biol Conserv 86:7781.[CrossRef]
Hebert PDN and Beaton MJ, 1989. Methodologies for allozyme analysis using cellulose acetate electrophoresis Beaumont, TX: Helena Laboratories.
Hewitt GM, 1999. Post-glacial recolonization of European biota. Biol J Linn Soc 68:87112.[CrossRef][ISI]
Houthuys R, De Moor G, and Sommé J, 1993. The shaping of the French-Belgian North Sea coast throughout recent geology and history In: Coastlines of the southern North Sea (Hillen R and Verhagen HJ, eds). New York: American Society of Civil Engineers.
Hsiao TH, 1989. Estimation of genetic variability amongst Coleoptera In: Electrophoretic studies on agricultural pests (Loxdale HD and den Hollander J, eds). Oxford: Clarendon; 143180.
Hurlbert SH, 1971. The nonconcept of species diversity: a critique and alternative parameters. Ecology 52:577586.[CrossRef][ISI]
Hyman PS and Parsons MS, 1992. Review of the scare and threatened Coleoptera of Great Britain Peterborough: Joint Nature Conservation Committee.
Knoll S and Rowell-Rahier M, 1998. Distribution of genetic variance and isolation by distance in two leaf beetle species: Oreina cacaliae and Oreina speciosissima. Heredity 81:412421.
Knoll S, Rowell-Rahier M, Mardulyn P, and Pasteels JM, 1996. Spatial genetic structure of leaf beetle species with special emphasis on alpine populations In: Chrysomelidae biology, vol 1: The classification, phylogeny and genetics (Jolivet PHA and Cox ML, eds). Amsterdam: SPB Academic Publishing; 379388.
Knutsen H, Rukke BA, Jorde PE, and Ims RA, 2000. Genetic differentiation among populations of the beetle Bolitophagus reticulatus (Coleoptera: Tenebrionidae) in a fragmented and a continuous landscape. Heredity 84:667676.
Mantel N, 1967. The detection of disease clustering and a generalized regression approach. Cancer Res 27:209220.[ISI][Medline]
Merilä J, Björklund M, and Baker AJ, 1996. Genetic population structure and gradual northward decline of genetic variability in the greenfinch (Carduelis chloris). Evolution 50:25482557.[CrossRef]
Nei M, 1978. Estimation of average heterozygosity and genetic distance from a small number of individuals. Genetics 89:583590.
Peacock MM and Smith AT, 1997. The effect of habitat fragmentation on dispersal patterns, mating behaviour, and genetic variation in a pika (Ochotona princeps) metapopulation. Oecologia 112:524533.[CrossRef]
Raymond M and Rousset F, 1995. GENEPOP version 1.2: a population genetics software for exact tests and ecumenicism. J Hered 86:248249.
Richardson BJ, Baverstock PR, and Adams M, 1986. Allozyme electrophoresis. A handbook for animal systematics and population studies San Diego: Academic Press.
Schneider S, Roessli D, and Excoffier L, 2000. ARLEQUIN version 2000: a software for population genetics data analysis Geneva, Switzerland: University of Geneva.
Schubel JR and Hirschberg DJ, 1978. Estuarine graveyards, climatic change and the importance of the estuarine environment In: Estuarine interactions (Wiley ML, ed). New York: Academic Press; 285303.
Segelbacher G and Storch I, 2002. Capercaillie in the Alps: genetic evidence of metapopulation structure and population decline. Mol Ecol 11:16691677.[Medline]
Slatkin M, 1994. Gene flow and population structure In: Ecological genetics (Real LA, ed). Princeton, NJ: Princeton University Press; 285303.
Taberlet P, Fumagalli L, Wust-Saucy A-G, and Cosson J-F, 1998. Comparative phylogeography and postglacial colonization routes in Europe. Mol Ecol 7:453464.[CrossRef][Medline]
Turin H, 2000. De Nederlandse loopkevers, verspreiding en oecologie (Coleoptera: Carabidae) Nederlandse Fauna3. Nationaal Natuurhistorisch Museum Naturalis, KNNV Uitgevery & EIS Nederland, Leiden.
Van Dongen S, Backeljau T, Matthysen E, and Dhondt AA, 1998. Genetic population structure of the winter moth (Operophtera brumata L.) (Lepidoptera, Geometridae) in a fragmented landscape. Heredity 80:92100.[CrossRef]
Walter R and Epperson BK, 2001. Geographic pattern of genetic variation in Pinus resinosa: area of greatest diversity is not the origin of postglacial populations. Mol Ecol 10:103111.[Medline]
Ward RD, Skibinski DOF, and Woodwark M, 1992. Protein heterozygosity, protein structure, and taxonomic differentiation In: Evolutionary biology (Hecht MK, Wallace B, and Macintyre RJ, eds). New York: Plenum; 73159.
Weir BS and Cockerham CC, 1984. Estimating F-statistics for the analysis of population structure. Evolution 43:13491368.
Williams BL, Brawn JD, and Paige KN, 2003. Landscape scale genetic effects of habitat fragmentation on a high gene flow species: Speyeria idalia (Nymphalidae). Mol Ecol 12:1120.[CrossRef][Medline]
Wright S, 1969. Evolution and the genetics of populations, vol 2: The theory of gene frequencies Chicago: University of Chicago Press.
Yeh F, Yang RC, and Boyle T, 1997. POPGENE. The user-friendly shareware for population genetic analysis Calgary, Alberta, Canada: University of Alberta, Center for International Forestry Research.
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||

