Genetic characterization of hybridization and introgression between anadromous rainbow trout (Oncorhynchus mykiss irideus) and coastal cutthroat trout (O. clarki clarki)

Interspecific hybridization represents a dynamic evolutionary phenomenon and major conservation problem in salmonid fishes. In this study we used amplified fragment length polymorphisms (AFLP) and mitochondrial DNA (mtDNA) markers to describe the extent and characterize the pattern of hybridization and introgression between coastal rainbow trout (Oncorhynchus mykiss irideus) and coastal cutthroat trout (O. clarki clarki). Hybrid individuals were initially identified using principle coordinate analysis of 133 polymorphic AFLP markers. Subsequent analysis using 23 diagnostic AFLP markers revealed the presence of F1, rainbow trout backcross, cutthroat trout backcross and later-generation hybrids. mtDNA analysis demonstrated equal numbers of F1 hybrids with rainbow and cutthroat trout mtDNA indicating reciprocal mating of the parental types. In contrast, rainbow and cutthroat trout backcross hybrids always exhibited the mtDNA from the recurrent parent, indicating a male hybrid mating with a pure female. This study illustrates the usefulness of the AFLP technique for generating large numbers of species diagnostic markers. The pattern of hybridization raises many questions concerning the existence and action of reproductive isolating mechanisms between these two species. Our findings are consistent with the hypothesis that introgression between anadromous populations of coastal rainbow and coastal cutthroat trout is limited by an environment-dependent reduction in hybrid fitness.     

Population subdivision in westslope cutthroat trout (Oncorhynchus clarki lewisi) at the northern periphery of its range: evolutionary inferences and conservation implications

Westslope cutthroat trout (Oncorhynchus clarki lewisi, Salmonidae) are native to the upper Columbia, Missouri, and South Saskatchewan river drainages of western North America and are at the northern periphery of their range in southeastern British Columbia, Canada. We examined geographical variation in allele frequencies at eight microsatellite loci in 36 samples of westslope cutthroat trout from British Columbia to assess levels of population subdivision and to test the hypothesis that different habitat types (principally mainstem vs. above migration barrier habitats) would influence levels of genetic diversity, genetic divergence among populations, and attainment of equilibrium between gene flow and genetic drift. Across all samples, the mean number of alleles per locus was 3.9 and mean expected heterozygosity was 0.56. Population subdivision was extensive with an overall Fst (theta) of 0.32. Populations sampled above migration barriers had significantly fewer alleles, lower expected heterozygosity, but greater average pairwise Fst than populations sampled from mainstem localities. We found evidence for isolation-by-distance from a significant correlation between genetic distance and geographical distance (r = 0.31), but the pattern was much stronger (r = 0.51) when above barrier populations and a population that may have been involved in headwater exchanges were removed. By contrast, isolation-by-distance was not observed when only above barrier populations were tested among themselves. Our data support the maintenance of separate demographic management strategies for westslope cutthroat trout inhabiting different river systems and illustrate how differing habitat structure (e.g. presence of migration barriers) may influence patterns of biodiversity and gene flow-drift equilibrium.     

Ten species specific microsatellite loci for Lahontan cutthroat trout,Oncorhynchus clarki henshawi

Ten polymorphic microsatellite loci (containing tri and tetra‐nucleotide repeats) were developed for the Lahontan cutthroat trout (Oncorhynchus clarki henshawi), a subspecies of cutthroat trout listed as threatened under the United States Endangered Species Act. Polymorphism was assessed for 445 individuals from 12 populations representing eight watersheds spread throughout the three Distinct Population Segments defined for this subspecies. All loci were polymorphic (X? = 19, range 7–26 alleles). All loci were in Hardy–Weinberg equilibrium (HWE) except for one locus (OCH 9) in a single population (P < 0.00014 after Bonferroni correction for multiple tests).     

Spatial partitioning and asymmetric hybridization among sympatric coastal steelhead trout (Oncorhynchus mykiss irideus), coastal cutthroat trout (O. clarki clarki) and interspecific hybrids

Hybridization between sympatric species provides unique opportunities to examine the contrast between mechanisms that promote hybridization and maintain species integrity. We surveyed hybridization between sympatric coastal steelhead (Oncorhynchus mykiss irideus) and coastal cutthroat trout (O. clarki clarki) from two streams in Washington State, Olsen Creek (256 individuals sampled) and Jansen Creek (431 individuals sampled), over a 3-year period. We applied 11 O. mykiss-specific nuclear markers, 11 O. c. clarki-specific nuclear markers and a mitochondrial DNA marker to assess spatial partitioning among species and hybrids and determine the directionality of hybridization. F1 and post-F1 hybrids, respectively, composed an average of 1.2% and 33.6% of the population sampled in Jansen Creek, and 5.9% and 30.4% of the population sampled in Olsen Creek. A modest level of habitat partitioning among species and hybrids was detected. Mitochondrial DNA analysis indicated that all F1 hybrids (15 from Olsen Creek and five from Jansen Creek) arose from matings between steelhead females and cutthroat males implicating a sneak spawning behaviour by cutthroat males. First-generation cutthroat backcrosses contained O. c. clarki mtDNA more often than expected suggesting natural selection against F1 hybrids. More hybrids were backcrossed toward cutthroat than steelhead and our results indicate recurrent hybridization within these creeks. Age analysis demonstrated that hybrids were between 1 and 4 years old. These results suggest that within sympatric salmonid hybrid zones, exogenous processes (environmentally dependent factors) help to maintain the distinction between parental types through reduced fitness of hybrids within parental environments while divergent natural selection promotes parental types through distinct adaptive advantages of parental phenotypes.     

A molecular analysis of hybridization between native westslope cutthroat trout and introduced rainbow trout in southeastern British Columbia, Canada

Restriction site variation in the Ikaros gene intron was used to assess the incidence of westslope cutthroat trout ( Oncorhynchus clarki lewisi ), rainbow trout ( O. mykiss ) and interspecific hybrids at 11 localities among eight streams tributary to the upper Kootenay River system in south-eastern British Columbia, Canada. Out of 356 fish assayed by this technique, hybrids ( n =16) were found at seven of the 11 sites across five different streams. Rainbow trout ( n =6) were found at two of the 11 sites. Analysis of hybrids with a second genetic marker (heat shock 71 intron) indicated that most represented either backcrosses to both westslope cutthroat and rainbow trout, or post F₁ hybrids. Mitochondrial DNA analysis indicated that hybrid matings occur between male rainbow trout and female westslope cutthroat trout and vice versa. Comparison of present hybridization in five tributaries relative to an allozyme-based analysis in the mid-1980s, that documented hybrids in only a single tributary of seven that were common to the two studies, suggests that hybridization and introgression has increased in upper Kootenay River tributaries. The present analysis is a conservative estimate of genetic interaction between the species because introgression was not tested in the majority of samples. Identification of genetically pure westslope cutthroat trout populations, and why they might be resistant to introgression from rainbow trout, are crucial conservation priorities for this unique subspecies of cutthroat trout.     

Subtle population structuring within a highly vagile marine invertebrate, the veined squid Loligo forbesi, demonstrated with microsatellite DNA markers

Microsatellite DNA markers were applied for the first time in a population genetic study of a cephalopod and compared with previous estimates of genetic differentiation obtained using allozyme and mitochondrial DNA (mtDNA) markers. Levels of genetic variation detected with microsatellites were much higher than found with previous markers (mean number of alleles per locus=10.6, mean expected heterozygosity ( H _E)=0.79; allozyme H _E=0.08; mtDNA restriction fragment length polymorphism (RFLP) H _E=0.16). In agreement with previous studies, microsatellites demonstrated genetic uniformity across the population occupying the European shelf seas of the North East Atlantic, and extreme genetic differentiation of the Azores population ( R _ST/ F _ST=0.252/0.245; allozyme F _ST=0.536; mtDNA F _ST=0.789). In contrast to other markers, microsatellites detected more subtle, and significant, levels of differentiation between the populations of the North East Atlantic offshore banks (Rockall and Faroes) and the shelf population ( R _ST=0.048 and 0.057). Breakdown of extensive gene flow among these populations is indicated, with hydrographic (water depth) and hydrodynamic (isolating current regimes) factors suggested as possible barriers to migration. The demonstration of genetic subdivision in an abundant, highly mobile marine invertebrate has implications for the interpretation of dispersal and population dynamics, and consequent management, of such a commercially exploited species. Relative levels of differentiation indicated by the three different marker systems, and the use of measures of differentiation (assuming different mutation models), are discussed.     

Molecular and quantitative trait variation within and among populations of the intertidal copepod Tigriopus californicus

Abstract While molecular and quantitative trait variation may be theoretically correlated, empirical studies using both approaches frequently reveal discordant patterns, and these discrepancies can contribute to our understanding of evolutionary processes. Here, we assessed genetic variation in six populations of the copepod Tigriopus californicus. Molecular variation was estimated using five polymorphic microsatellite loci, and quantitative variation was measured using 22-life history and morphometric characters. Within populations, no correlation was found between the levels of molecular variation (heterozygosity) and quantitative variation (heritability). Between populations, quantitative subdivision ( Q _ST) was correlated with molecular subdivision when measured as F _ST but not when measured as R _ST. Unlike most taxa studied to date, the overall level of molecular subdivision exceeded the level of quantitative subdivision ( F _ST= 0.80, R_ST = 0.89, Q _ST = 0.30). Factors that could contribute to this pattern include stabilizing or fluctuating selection on quantitative traits or accelerated rates of molecular evolution.     

Erosion of interspecific reproductive barriers resulting from hatchery supplementation of rainbow trout sympatric with cutthroat trout

The frequency of hybridization between cutthroat (Onchorhynchus clarki clarki) and rainbow (O. mykiss irideus) trout from coastal habitats in British Columbia, Canada, was examined in seven populations where the two species are sympatric with no history of rainbow trout stocking and compared with areas where native rainbow trout populations have been supplemented with hatchery fish (three populations). Four nuclear markers were used to identify each species and interspecific hybrids and one mitochondrial marker showed the direction of gene exchange between species. The frequency of hybrids was significantly higher (Fisher exact test, P < 0.001) in river systems where hatchery rainbow trout have been introduced (50.6% hybrids) than in populations where the two species naturally co-occur without supplementation (9.9% hybrids).     

Genetic Integrity and Microgeographic Population Structure of Westslope Cutthroat Trout, Oncorhynchus clarki lewisi, in the Pend Oreille Basin in Washington

Five microsatellite DNA loci (Ots-101 ^*,Ots-107 ^*,Oki-10 ^*, Ogo-3 ^*, and FGT-3 ^*) were screened to evaluate the genetic characteristics and population structure for cutthroat trout from eight tributaries of the Pend Oreille River in northeastern Washington and to compare these collections with two hatchery stocks of westslope cutthroat trout, Oncorhynchus clarki lewisi, Yellowstone cutthroat trout, Oncorhynchus clarki bouvieri and a hatchery rainbow trout, Oncorhynchus mykiss, strain that have been stocked in northeastern Washington. Relatively high levels of variation (numbers of alleles and heterozygosity) were observed in all collections and allele frequencies were quite variable among collections. Evidence of limited introgression by rainbow and/or Yellowstone cutthroat was found at several locations. Both F_ST values and tests of genetic differentiation indicated the existence of numerous, reproductively isolated populations. The population in Slate Creek was very similar to the Kings Lake Hatchery strain, and we conclude that this similarity is the result of historical introductions of this hatchery strain into what was presumably a stream without a native cutthroat population. In one stream, differences in introgression and allele frequencies were found above and below a barrier falls. Because of the substantial level of population differentiation observed among the various collections, we recommend that management and conservation actions be focused at the level of individual streams in order to maintain the productivity and genetic character of the existing populations of cutthroat trout.     

A suite of twelve single nucleotide polymorphism markers for detecting introgression between cutthroat and rainbow trout

A suite of 12 subspecies and species-specific single nucleotide polymorphism (species-specific SNP) markers was developed to distinguish rainbow trout (RT) Oncorhynchus mykiss from the four major subspecies of cutthroat trout: westslope cutthroat trout (WCT) Oncorhynchus clarki lewisi, Yellowstone cutthroat trout (YCT) Oncorhynchus clarki bouvieri, coastal cutthroat trout (CCT) Oncorhynchus clarki clarki, Lahontan cutthroat trout (LCT) Oncorhynchus clarki henshawi, and their hybrids. Several of the markers were linked to help strengthen hybrid determinations, and sex-specific species-specific SNP assays were also developed.     

Allozyme Diversity in Populations of Cymbidium goeringii (Orchidaceae)

Abstract: Using 14 allozyme loci, we investigated levels of genetic diversity within populations, and degree of genetic divergence among 24 populations of Cymbidium goeringii (Orchidaceae) in Korea and Japan. Cymbidium goeringii maintains high levels of genetic diversity both at population (mean expected heterozygosity, H _e = 0.238) and species levels (0.260). Means of H _e found in 24 populations were not significantly different from each other. About 90 % of the total variation in the species is common to all populations (mean G _ST = 0.108). No unique allele was found in any population. The indirect estimate of gene flow based on the mean G _ST was high ( Nm = 2.06). Nei's genetic identities for pairs of populations had high values (mean = 0.974 [SD = 0.013]). The Mantel-Z test showed a significant correlation between genetic distance and geographic distance. However, the mean G _ST value between 17 populations in Korea and seven Japanese populations was relatively low (0.029), even though the land connection between the southern Korean peninsula and southern Japanese archipelagos has not existed since the middle Pleistocene. Large numbers of small seeds of C. goeringii might travel long distances by wind from populations to populations both in Korea and Japan, increasing genetic diversity within populations and maintaining low genetic differentiation among populations.     

Development and evaluation of 200 novel SNP assays for population genetic studies of westslope cutthroat trout and genetic identification of related taxa

DNA sequence data were collected and screened for single nucleotide polymorphisms (SNPs) in westslope cutthroat trout (Oncorhynchus clarki lewisi) and also for substitutions that could be used to genetically discriminate rainbow trout (O. mykiss) and cutthroat trout, as well as several cutthroat trout subspecies. In total, 260 expressed sequence tag‐derived loci were sequenced and allelic discrimination genotyping assays developed from 217 of the variable sites. Another 50 putative SNPs in westslope cutthroat trout were identified by restriction‐site‐associated DNA sequencing, and seven of these were developed into assays. Twelve O. mykiss SNP assays that were variable within westslope cutthroat trout and 12 previously published SNP assays were also included in downstream testing. A total of 241 assays were tested on six westslope cutthroat trout populations (N = 32 per population), as well as collections of four other cutthroat trout subspecies and a population of rainbow trout. All assays were evaluated for reliability and deviation from Hardy–Weinberg and linkage equilibria. Poorly performing and duplicate assays were removed from the data set, and the remaining 200 assays were used in tests of population differentiation. The remaining markers easily distinguished the various subspecies tested, as evidenced by mean G_ST of 0.74. A smaller subset of the markers (N = 86; average G_ST = 0.40) was useful for distinguishing the six populations of westslope cutthroat trout. This study increases by an order of magnitude the number of genetic markers available for the study of westslope cutthroat trout and closely related taxa and includes many markers in genes (developed from ESTs).     

Clinal genetic variation and isolation by distance in the European eel Anguilla anguilla (L.)

The genetic variability and structure of the European eel ( Anguilla anguilla L.) in populations throughout Europe was reassessed using 15 allozymic loci, seven of which were polymorphic. Seven sites were sampled on a latitudinal gradient across the natural continental range, extending from southern France to southern Norway. Heterozygosity ( H _e = 0.05) and level of polymorphism (P = 0.43) were comparable to other marine fish. Populations were poorly differentiated ( G _ST = 0.014, F _ST = 0.002), which is not surprising considering the high dispersal capability of the European eel. However, a significant geographical cline was detected at two alleles ( IDH-1 * 100 and GPI-1 * 110 ), and genetic distances ( D _CE) were concordant with geographical coastal distances. Mantel tests, pairwise F _ST's and multidimensional scaling analyses identify three distinct groups: Northern Europe, Western Europe and the Mediterranean Sea. We propose that the clinal genetic structure in the European eel may be due to (1) isolation by distance (as recently detected with microsatellites), (2) temporal reproductive separation, (3) post-larval selective forces, (4) contact between formerly separated groups or (5) some combination thereof.  © 2002 The Linnean Society of London, Biological Journal of the Linnean Society , 2002, 77 , 509–521     

Widespread gene flow and high genetic variability in populations of water voles Arvicola terrestris in patchy habitats

Theory predicts that the impact of gene flow on the genetic structure of populations in patchy habitats depends on its scale and the demographic attributes of demes (e.g. local colony sizes and timing of reproduction), but empirical evidence is scarce. We inferred the impact of gene flow on genetic structure among populations of water voles Arvicola terrestris that differed in average colony sizes, population turnover and degree of patchiness. Colonies typically consisted of few reproducing adults and several juveniles. Twelve polymorphic microsatellite DNA loci were examined. Levels of individual genetic variability in all areas were high ( H _O= 0.69–0.78). Assignments of juveniles to parents revealed frequent dispersal over long distances. The populations showed negative F _IS values among juveniles, F _IS values around zero among adults, high F _ST values among colonies for juveniles, and moderate, often insignificant, F _ST values for parents. We inferred that excess heterozygosity within colonies reflected the few individuals dispersing from a large area to form discrete breeding colonies. Thus pre-breeding dispersal followed by rapid reproduction results in a seasonal increase in differentiation due to local family groups. Genetic variation was as high in low-density populations in patchy habitats as in populations in continuous habitats used for comparison. In contrast to most theoretical predictions, we found that populations living in patchy habitats can maintain high levels of genetic variability when only a few adults contribute to breeding in each colony, when the variance of reproductive success among colonies is likely to be low, and when dispersal between colonies exceeds nearest-neighbour distances.     

Genetic variation and bill size dimorphism in a passerine bird, the reed bunting Emberiza schoeniclus

In passerine birds morphological differentiation in bill size within species is not commonly observed. Bill size is usually associated with a trophic niche, and strong differences in it may reflect the process of genetic differentiation and, possibly, speciation. We used both mitochondrial DNA (mtDNA) and nuclear microsatellites to study genetic variation between two subspecies of reed bunting, Emberiza schoeniclus schoeniclus and E.s. intermedia , along their distributional boundary in western Europe. These two subspecies are characterized by a high dimorphism in bill size and, although breeding populations of the two subspecies are found very close to each other in northern Italy, apparently no interbreeding occurs. The observed morphological pattern between the two subspecies may be maintained by geographically varying selective forces or, alternatively, may be the result of a long geographical separation followed by a secondary contact. MtDNA sequences of cytochrome b and ND5 (515 bp) showed little variation and did not discriminate between the two subspecies, indicating a divergence time of less than 500 000 years. The analysis of four microsatellite loci suggested a clear, although weak, degree of genetic differentiation in the large- and small-billed populations, as indicated by F _ST and R _ST values and genetic distances. The correlation between bill size and genetic distance between populations remained significant after accounting for the geographical distances between sampling localities. Altogether, these results indicate a very recent genetic differentiation between the two bill morphs and suggest that a strong selection for large bills in the southern part of the breeding range is probably involved in maintaining the geographical differentiation of this species.     

Characterization of 13 microsatellites for Lahontan cutthroat trout (Oncorhynchus clarki henshawi) and cross-amplification in six other salmonids

Thirteen newly developed tri- and tetranucleotide repeat microsatellite markers were developed for Lahontan cutthroat trout (Oncorhynchus clarki henshawi), a threatened subspecies endemic to the Lahontan hydrographic basin in the western USA. These loci are highly polymorphic with five to 30 alleles per locus and observed heterozygosities ranging from 0.4 to 0.7. Cross-species amplification of these markers was most successful in the closely related rainbow trout, Oncorhynchus mykiss, with only three loci amplifying in brown trout, Salmo trutta. Nonoverlapping allelic distributions for many of these loci among the six salmonid species screened suggest these markers may be useful for hybrid determination.     

Population divergence in the amphicarpic species Amphicarpaea edgeworthii Benth. (Fabaceae): microsatellite markers and leaf morphology

Comparative analyses of the genetic differentiation in microsatellite markers ( F _ST) and leaf morphology characters ( Q _ST) of Amphicarpaea edgeworthii Benth. were conducted to gain insight into the roles of random processes and natural selection in the population divergence. Simple sequence repeat analyses on 498 individuals of 19 natural populations demonstrate that a significant genetic differentiation occurs among populations (mean F _ST = 0.578), and A. edgeworthii is a highly self-fertilized species (mean selfing rate s  = 0.989). The distribution pattern of genetic diversity in this species shows that central populations possess high genetic diversity (e.g. population WL with H _E = 0.673 and population JG with H _E = 0.663), whereas peripheral ones have a low H _E as in population JD (0.011). The morphological divergence of leaf shape was estimated by the elliptical Fourier analysis on the data from 11 natural and four common garden populations. Leaf morphology analyses indicate the morphological divergence does not show strong correlation with the genetic differentiation ( R  = 0.260, P  = 0.069). By comparing the 95% confidence interval of Q _ST with that of F _ST, Q _ST values for five out of 12 quantitative traits are significantly higher than the average F _ST value over eight microsatellite loci. The comparison of F _ST and Q _ST suggests that two kinds of traits can be driven by different evolutionary forces, and the population divergence in leaf morphology is shaped by local selections.  © 2009 The Linnean Society of London, Biological Journal of the Linnean Society , 2009, 96 , 505–516.     

Population genetic structure of the butterfly Melitaea didyma (Nymphalidae) along a northern distribution range border

The population genetic structure of the butterfly Melitaea didyma was studied along the northern distribution range border in Central Germany by means of allozyme electrophoresis. Individuals were sampled from a total of 21 habitat patches from four regions, and two provinces. Sampling was designed to estimate local vs. regional differentiation. High levels of variability were found, H _e= 0.14–0.21. The mean expected sample heterozygosity from one region, Mosel, was significantly lower than from the Hammelburg region, H _e= 0.17 and 0.19, respectively. Two hierarchical levels of genetic differentiation were found. Within regions individuals sampled from different patches behaved as belonging to one population with high levels of gene flow (Hammelburg F _ST= 0.015, Mosel F _ST= 0.044), though local isolation barriers did create a substructuring of these populations. The inbreeding coefficients, F _IS, were constant over all sample levels, suggesting a similar distribution of habitat patches within regions. Between regions gene flow was limited. An isolation by distance analysis indicated that the hierarchical structure, at the provincial level, may be breaking down due to isolation of regional populations. A more general observation was that the sampling design may greatly have influenced the estimation of genetic differentiation. Depending on which samples were included, overall F _ST estimates ranged from 0.059–0.090.     

Comparison of quantitative and molecular genetic variation of native vs. invasive populations of purple loosestrife (Lythrum salicaria L., Lythraceae)

Study of adaptive evolutionary changes in populations of invasive species can be advanced through the joint application of quantitative and population genetic methods. Using purple loosestrife as a model system, we investigated the relative roles of natural selection, genetic drift and gene flow in the invasive process by contrasting phenotypical and neutral genetic differentiation among native European and invasive North American populations ( Q _ST −  F _ST analysis). Our results indicate that invasive and native populations harbour comparable levels of amplified fragment length polymorphism variation, a pattern consistent with multiple independent introductions from a diverse European gene pool. However, it was observed that the genetic variation reduced during subsequent invasion, perhaps by founder effects and genetic drift. Comparison of genetically based quantitative trait differentiation ( Q _ST) with its expectation under neutrality ( F _ST) revealed no evidence of disruptive selection ( Q _ST >  F _ST) or stabilizing selection ( Q _ST <  F _ST). One exception was found for only one trait (the number of stems) showing significant sign of stabilizing selection across all populations. This suggests that there are difficulties in distinguishing the effects of nonadaptive population processes and natural selection. Multiple introductions of purple loosestrife may have created a genetic mixture from diverse source populations and increased population genetic diversity, but its link to the adaptive differentiation of invasive North American populations needs further research.     

Genetic differentiation within and between isolated Algerian subpopulations of Barbary macaques (Macaca sylvanus): evidence from microsatellites

This study of wild-living Algerian Barbary macaques ( Macaca sylvanus ) was designed to examine genetic variability in subpopulations isolated in residual forest patches, in an attempt to obtain data on the effects of habitat fragmentation. The wild population of this species (estimated at a maximum of 15 000) is vulnerable and this study therefore has direct relevance for conservation measures. Data from five microsatellite loci were analysed for 159 individuals from nine different groups living in four isolates in Algeria. Genetic polymorphism was found to be relatively high (4–12 alleles per locus) compared with other genetic markers used in previous studies of this species; mean expected heterozygosity was 65%. The four isolates are all genetically distinct ( F _ST = 0; P < 0.001). Indeed, the results suggest that dispersal is limited even between some social groups within a single isolate. Genetic distances based on models not assuming stepwise mutation ( F _ST and chord distance) gave very similar results and are highly correlated with geographical distances within one isolate but not between isolates. This may indicate that isolation by distance exerts a significant influence within an isolate but that genetic drift prevails between the four isolates. After allowing for variation in sample size, we found no evidence of reduced allelic diversity within small isolates that may have been separated for 250 years or more. The surviving population of Algerian Barbary macaques taken as a whole still shows marked variability in microsatellite alleles, but maintenance of genetic variability over the long term will surely require effective protection of all isolates.