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.     

Microsatellite variation and population structure of the moor frog (Rana arvalis) in Scandinavia

Several recent studies have found amphibian populations to be genetically highly structured over rather short geographical distances, and that the rate of genetically effective dispersal may differ between the sexes. However, apart from the common frog ( Rana temporaria ) little is known about the genetic structuring and sex-biased dispersal in northern European amphibians. We investigated the patterns of genetic diversity and differentiation within and among Scandinavian populations of the moor frog ( Rana arvalis ) using microsatellite markers. The genetic diversity within local R. arvalis populations was not a simple linear negative function of latitude but a convex one: genetic diversity peaked in mid-latitude populations, and declined thereafter dramatically towards the north. The average degree of genetic differentiation among populations ( F _ST = 0.14) was lower than that observed for the common frog ( F _ST = 0.21), though the pattern of isolation by distance was similar for both species. Contrary to common frogs, no evidence for female-biased dispersal was found. The results reinforce the view that amphibian populations are—in general—highly structured over relatively small geographical distances, even in comparatively recently colonized areas.     

Microsatellite analysis of genetic population structure in an endangered salmonid: the coastal cutthroat trout (Oncorhynchus clarki clarki)

The genetic population structure of coastal cutthroat trout ( Oncorhynchus clarki clarki ) in Washington state was investigated by analysis of variation in allele frequencies at six highly polymorphic microsatellite loci for 13 anadromous populations, along with one outgroup population from the Yellowstone subspecies ( O. clarki bouvieri) (mean heterozygosity = 67%; average number of alleles per locus = 24). Tests for genetic differentiation revealed highly significant differences in genotypic frequencies for pairwise comparisons between all populations within geographical regions and overall population subdivision was substantial ( F _ST = 0.121, R _ST = 0.093), with 44.6% and 55.4% of the among-population diversity being attributable to differences between streams ( F _SR = 0.054) and between regions ( F _RT = 0.067), respectively. Analysis of genetic distances and geographical distances did not support a simple model of isolation by distance for these populations. With one exception, neighbour-joining dendrograms from the Cavalli-Sforza and Edwards' chord distances and maximum likelihood algorithms clustered populations by physiogeographic region, although overall bootstrap support was relatively low (53%). Our results suggest that coastal cutthroat trout populations are ultimately structured genetically at the level of individual streams. It appears that the dynamic balance between gene flow and genetic drift in the subspecies favours a high degree of genetic differentiation and population subdivision with the simultaneous maintenance of high heterozygosity levels within local populations. Results are discussed in terms of coastal cutthroat trout ecology along with implications for the designation of evolutionarily significant units pursuant to the US Endangered Species Act of 1973 and analogous conservation units.     

Extremely low levels of genetic diversity in the terrestrial orchid Epipactis thunbergii (Orchidaceae) in South Korea: implications for conservation

Levels of allozyme variation, population genetic structure, and fine-scale genetic structure (FSGS) of the rare, both sexually and clonally reproducing terrestrial orchid Epipactis thunbergii were examined for eight ( N  = 734) populations in a 20 × 20-km area in South Korea. Twenty-three putative allozyme loci resolved from 15 enzyme systems were used. Extremely low levels of allozyme variation were found within populations: the mean frequency of polymorphic loci was 3.8% [isocitrate dehydrogenase ( Idh-2 ) with two alleles was polymorphic across populations], the mean number of alleles per locus was 1.04, and the mean expected heterozygosity was 0.013. The overall fixation index was not significantly different from zero ( F _IS = 0.069), although the species is self-compatible. However, a significantly high degree of population differentiation was found between populations at Idh-2 ( F _ST = 0.388) in the studied area. Furthermore, spatial autocorrelation analyses revealed a significant FSGS (up to 3 m) within populations. These observations suggest that the main explanatory factors for the extremely low levels of genetic diversity and the shaping of the population genetic structure of E. thunbergii are genetic drift as a result of a small effective population size, a restricted gene flow, and the isolation of populations. Considering the current genetic structure of E. thunbergii , three guidelines are suggested for the development of conservation strategies for the species in South Korea: (1) protection of habitats of standing populations; (2) prohibition by law of any collection of E. thunbergii ; and (3) protection of nearby pollinator populations, given the fact that fruit set in natural habitats is very low.  © 2007 The Linnean Society of London, Botanical Journal of the Linnean Society , 2007, 155 , 161–169.     

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.     

Population genetic patterns suggest a behavioural change in wild common frogs (Rana temporaria) following disease outbreaks (Ranavirus)

We use 14 microsatellite loci to investigate the impact of a viral disease ( Ranavirus ) on the population genetic structure of wild common frogs ( Rana temporaria ). Populations with a history of Ranavirus mortalities (and 83% declines in the number of frogs) were compared with populations with no history of infection. Infected ponds showed significantly elevated F _IS (homozygote excess), significantly reduced relatedness, and no detectable effect on allelic richness. We hypothesize that the elevated F _IS and reduced relatedness are consequences of assortative mating, and that allelic richness is maintained by immigration from nearby populations. Simulations indicate that the elevated F _IS cannot be explained by population size reductions, but can indeed be explained by assortative mating (even if a mate choice locus is unlinked to the genetic markers). While the majority of studies consider demographic outcomes following disease outbreaks, our results indicate that emerging infectious diseases could also result in behavioural changes.     

Microsatellite analysis reveals high genetic diversity but low genetic structure in Ethiopian indigenous cattle populations

Ethiopian cattle are under threat from uncontrolled mating practices and are at high risk of becoming genetically homogeneous. Therefore, to evaluate genetic diversity, population structure and degree of admixture, 30 microsatellite markers were genotyped using 351 DNA samples from 10 Ethiopian cattle populations and the Holstein breed. The mean number of alleles per cattle population ranged from 6.93 ± 2.12 in Sheko to 7.50 ± 2.35 in Adwa. The mean observed and expected heterozygosities were 0.674 ± 0.015 and 0.726 ± 0.019 respectively. Ethiopian cattle populations have maintained a high level of within-population genetic differentiation (98.7%), the remainder being accounted for by differentiation among populations (1.3%). A highly significant deficiency in heterozygotes was detected within populations ( F _IS = 0.071; P  <   0.001) and total inbreeding ( F _IT = 0.083; P  <   0.001). The study populations were highly admixed but distinct from pure Bos taurus and Bos indicus breeds. The various levels of admixture and high genetic diversity make Ethiopian cattle populations suitable for future genetic improvement and utilization under a wide range of agro-ecologies in Ethiopia.     

Genetic structure within and between island populations of the flightless cormorant (Phalacrocorax harrisi)

We assessed colony- and island-level genetic differentiation for the flightless cormorant ( Phalacrocorax harrisi ), an endangered Galápagos endemic that has one of the most limited geographical distributions of any seabird, consisting of only two adjacent islands. We screened 223 individuals from both islands and nine colonies at five microsatellite loci, recovering 23 alleles. We found highly significant genetic differentiation throughout the flightless cormorant's range on Fernandina and Isabela Islands (global F _ST = 0.097; P  < 0.0003) both between islands (supported by Bayesian analyses, F _ST and R _ST values) and within islands (supported only by F _ST and R _ST values). An overall pattern of isolation-by-distance was evident throughout the sampled range ( r =  0.4169, one-sided P  ≤ 0.02) and partial Mantel tests of this relationship confirmed that ocean is a dispersal barrier ( r =  0.500, one-sided P  ≤ 0.003), especially across the 5-km gap between the two islands. The degree of detected genetic differentiation among colonies is surprising, given the flightless cormorant's limited range, and suggests a role for low vagility, behavioural philopatry, or both to limit dispersal where physical barriers are absent. We argue that this population should be managed as at least two genetic populations to better preserve the species-level genetic diversity, but, for demographic reasons, advocate the continued conservation of all breeding colonies.     

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.     

Understanding population structure and historical demography in a conservation context: population genetics of an endangered fern

The construction of the world's largest hydroelectric scheme across the Yangtze River, the Three Gorges Dams (TGD), in the centre of a southern-central Chinese biodiversity hot spot, the Three Gorges Reservoir Area (TGRA), has attracted international concern and conservation action. To examine whether landscape changes to date have impacted regional flora, and to establish long-term monitoring baselines, we assessed the distribution and dynamics of an endangered and TGRA endemic fern, Adiantum reniforme var. sinense . For eight nuclear microsatellites, high levels of genetic diversity ( H _E = 0.653–0.781) and slightly elevated inbreeding ( F _IS = 0.077–0.197) were found across 13 surveyed populations. The population history of this fern is characterized by a balance of gene flow and genetic drift, where historical dispersal, inferred from coalescent ( F =  0.129) and genetic differentiation ( F _ST = 0.094 and R _ST = 0.180) approaches, is moderate, reflecting an isolation by distance relationship. Importantly, most populations exhibited mutation-drift disequilibrium, suggesting a recent population decline, which is congruent with the known demographic history of the species following dam-related activities. Based on these results, populations of A. reniforme var. sinense are expected to lose genetic diversity and increase genetic structure as dam-related activities decrease size and increase genetic isolation of remnants.     

Genetic structuring and transfer of marine dinoflagellate Cochlodinium polykrikoides in Japanese and Korean coastal waters revealed by microsatellites

To determine the process of population expansion and ascertain the origin of the Sea of Japan population, in a noxious red tide forming dinoflagellate Cochlodinium polykrikoides , 13 samples, isolated from 11 different localities in Japanese and Korean coasts, were analysed using 10 polymorphic microsatellites. Analyses by nonmetric multidimensional scaling plots of pairwise F _ST, global amova , and genetic admixture analysis identified three clusters — the Sea of Japan populations, Yatsushiro Sea (Kumamoto Pref.) populations, and other populations — indicating genetic structuring of the 13 samples into three distinct populations. In the proportion of shared alleles by pairwise individuals ( P _SAxy) analyses between the Sea of Japan and the other samples, P _SAxy was extremely low compared with that among the Sea of Japan or among other samples, indicating that a large genetic barrier has occurred between the populations. No significant relationship of isolation-by-distance patterns and almost no genetic distance were detected between pairwise samples of the Sea of Japan, although there is a maximal distance of > 600 km between samples. In addition, P _SAxy data among the samples were extremely high compared with those among other samples, clearly showing that a large-scale transfer from west to east has occurred via the Tushima Warm Current. In the P _SAxy data of the Seto Inland Sea and Pacific samples, individuals showing relatively high P _SAxy were concentrated in the three areas of Nagasaki, Harima, and Mie, suggesting that frequent transfer may have occurred by human-assisted dispersal, although Nagasaki and Mie are separated by a distance of approximately 700 km.     

Population genetic structure and effective population size of ayu (Plecoglossus altivelis), an amphidromous fish

The temporal and spatial population genetic structure of ayu Plecoglossus altivelis (Salmoniformes: Plecoglossidae), an amphidromous fish, was examined using analysis of variation at six microsatellite DNA loci. Intracohort genetic diversities, as measured by the number of alleles and heterozygosity, were similar among six cohorts (2001–2006) within a population (Nezugaseki River), with the mean number of alleles per cohort ranging from 11·0 to 12·5 and the expected heterozygosity ranging from 0·74 to 0·77. Intrapopulational genetic diversities were also similar across the three studied populations along the 50 km coast, with the mean number of alleles and the expected heterozygosity ranging from 11·33 to 11·67 and from 0·75 to 0·76, respectively. The authors observed only one significant difference in pair-wise population differentiation ( F _ST-value) between the cohorts within a population and among three populations. Estimates of the effective population size ( N _e) based on maximum-likelihood method yielded small values (ranging from 94·8 to 135·5), whereas census population size ranged from c. 4800 to 24 000. As a result, the ratio of annual effective population sizes to census population size ( N _e/ N ) ranged from 0·004 to 0·023. These estimates of N _e/ N agree more closely with estimates for marine fishes than that of the larger estimates for freshwater fishes. The present study suggests that ayu which is highly fecund and shows low survival during the early life stages is also characterized by having low value of N _e/ N , similar to marine species with a pelagic life cycle.     

Genetic variability in the European bison (Bison bonasus) population from Białowieża forest over 50 years

The aim of this study was to assess potential post-bottleneck temporal genetic differentiation following the reintroduction of the species into the Białowieża Forest. Variability of 12 polymorphic microsatellite markers was analysed for 178 individuals born between 1955 and 2005, divided by birth year into five temporal groups. Low overall allelic richness (AR) per locus (AR = 2.0) and a low overall expected heterozygosity (H_E = 0.28) were observed. The overall F _IS was not significantly different from zero. The mean F _IS values were, however, significantly different from zero for individuals born between 1955 and 1965 ( F _IS = 0.19). A Bayesian computation was used to estimate effective population size (N_e) for each temporal group. We observed relatively small values of N_e ranging from 7.9 to 28.4. The low N_e values confirm that, despite a rapid growth of the bison population following the founder event, N_e increased only slowly.  © 2009 The Linnean Society of London, Biological Journal of the Linnean Society , 2009, 97 , 801–809.     

Genetic differentiation amongst fragrant orchids (Gymnadenia conopsea s.l.) in the British Isles

Genetic variation was examined in five microsatellite loci to seek evidence of genetic differentiation and restricted gene flow that would support the taxonomic division of Gymnadenia into three species ( G. borealis , G. conopsea , and G. densiflora ). A total of 107 alleles was detected in 17 populations from England, Scotland, and Ireland. The mean expected heterozygosities ranged from 0.48 to 0.81. The differentiation in allele frequencies amongst populations that had been assigned to each taxon on the basis of morphology was sufficiently large to support the taxa as distinct species. Phylogenetic trees based on microsatellite allele frequencies, as well as assignment tests, supported the existence of three distinct groups with at least partial restriction of gene flow between them. There was substantial homozygote excess, leading to high F _IS estimates, for most loci in most populations. This is unlikely to have been a result of widespread null alleles, and more probably reflects a high level of inbreeding in G. conopsea . This inference requires further investigation. The implications of the results of this and other taxonomic studies for the conservation of Gymnadenia in Britain are discussed.  © 2007 The Linnean Society of London, Botanical Journal of the Linnean Society , 2007, 155 , 349–360.     

Population structure of Moroccan water frogs: genetic cohesion despite a fragmented distribution

This study analyses the allozymic variation of 20 presumptive loci in eight populations of Rana saharica from Morocco. Populations were collected from the very different climatic zones of this country: the Rif area, the Atlas mountains and the desert. Moroccan water-frog populations are genetically well differentiated from the geographically closed Algerian populations. Thus, to check if such a differentiation process is taking place within Moroccan water frogs, we attempted to analyse the genetic structure and patterns of gene flow of Moroccan populations, by means of estimates of Hardy-Weinberg equilibrium, F-statistics and indirect measures of gene flow. F_st(0.250) and F_is(0.254) values were similar, which means that both intra and interpopulation differentiation contribute equally to the amount of genetic divergence revealed. F_is values indicated some degree of structure within ponds, which is possibly related to the homing behaviour of some amphibians. On the other hand, F_st and genetic distances between populations were not very high. Despite the low levels of gene flow estimated, together with the homing behaviour revealed and the spatially discontinuous distribution, it was found that genetic differentiation among populations was not as high as expected. The likelihood of genetic homogeneity being the consequence of continuous population extinction and recolonization events is discussed.     

Fine‐scale genetic structure overrides macro‐scale structure in a marine snail: nonrandom recruitment,demographic events or selection?

The planktotrophic littorinid species Littoraria flava occupies a continuous habitat on rocky shores close to brackish and freshwater sources. Previous studies of this species have shown a moderate genetic structure over a broad geographical scale, with high deviations from Hardy–Weinberg expectations in many allozymic loci. Local-scale subdivision in marine species with a long dispersal phase is unexpected, but occasionally found. Using a horizontal transect at three locations, we examined whether microscale and short-term subdivision also occurred in L. flava populations and, if so, whether this could explain the Hardy–Weinberg deviations. Littoraria flava showed even more structuring on a microgeographical scale (4–300 m) than on a large-scale (> 200 km). The Ewens–Watterson neutrality test showed that 18% of the tests deviated significantly from the neutrality model. A homogeneity test for each locus across samples within transects showed homogeneous and high F _IS values in many loci. These results and the apparent genetic patchiness within transects suggest that asynchronous spawning associated with recurrent colonizations in L. flava can explain the local differentiation without a recognizable pattern. In addition, there could be a balance between these factors and diversifying selection acting on different loci at different times and localities. © 2007 The Linnean Society of London, Biological Journal of the Linnean Society , 2007, 91 , 23–36.     

Contrasting levels of genetic diversity between the common, self-compatible Liparis kumokiri and rare, self-incompatible Liparis makinoana (Orchidaceae) in South Korea

Levels of allozyme variation and intrapopulation spatial genetic structure of the two terrestrial clonal orchids Liparis kumokiri , a self-compatible relatively common species, and L. makinoana , a self-incompatible rare species, were examined for 17 ( N  = 1875) and four ( N  = 425) populations, respectively, in South Korea. Populations of L. makinoana harboured high levels of genetic variation ( H _e = 0.319) across 15 loci. In contrast, L. kumokiri exhibited a complete lack of allozyme variation ( H _e = 0.000). Considering the lack of genetic variability, it is suggested that current populations of L. kumokiri in South Korea originated from a genetically depauperate ancestral population. For L. makinoana , a significant deficit of heterozygosity (mean F _IS = 0.198) was found in population samples excluding clonal ramets, suggesting that pollen dispersal is localized, generating biparental inbreeding. The significant fine-scale genetic structuring (≤ 2 m) found in a previous study, in addition to the moderate levels of population differentiation ( F _ST = 0.107) and the significant relationship between genetic and geographical distances ( r  = 0.680) found here, suggests a leptokurtic distribution of seed dispersal for L. makinoana . Although populations of L. makinoana harbour high levels of genetic variation, they are affected by a recent genetic bottleneck. This information suggests that genetic drift and limited gene flow could be the main evolutionary forces for speciation of a species-rich genus such as Liparis .  © 2007 The Linnean Society of London, Botanical Journal of the Linnean Society , 2007, 153 , 41–48.     

Population differentiation in an annual legume: genetic architecture

Abstract. The presence or absence of epistasis, or gene interaction, is explicitly assumed in many evolutionary models. Although many empirical studies have documented a role of epistasis in population divergence under laboratory conditions, there have been very few attempts at quantifying epistasis in the native environment where natural selection is expected to act. In addition, we have little understanding of the frequency with which epistasis contributes to the evolution of natural populations. In this study we used a quantitative genetic design to quantify the contribution of epistasis to population divergence for fitness components of a native annual legume, Chamaecrista fasciculata . The design incorporated the contrast of performance of F₂ and F₃ segregating progeny of 18 interpopulation crosses with the F₁ and their parents. Crosses were conducted between populations from 100 m to 2000 km apart. All generations were grown for two seasons in the natural environment of one of the parents. The F₁ often outperformed the parents. This F₁ heterosis reveals population structure and suggests that drift is a major contributor to population differentiation. The F₂ generation demonstrated that combining genes from different populations can sometimes have unexpected positive effects. However, the F₃ performance indicated that combining genes from different populations decreased vigor beyond that due to the expected loss of heterozygosity. Combined with previous data, our results suggest that both selection and drift contribute to population differentiation that is based on epistatic genetic divergence. Because only the F₃ consistently expressed hybrid breakdown, we conclude that the epistasis documented in our study reflects interactions among linked loci.     

High Genetic Diversity and Differentiation in Relict Lowland Populations of Gentianella austriaca (A. and J. Kern.) Holub (Gentianaceae)

Abstract: The biennial Gentianella austriaca (A. & J. Kern.) Holub, representing a nutrient-poor grassland taxon of low competition power, is becoming rare in the lowlands of eastern Austria due to changes in land use. To estimate effects of isolation and decreasing population sizes, as well as evolutionary relationships, we investigated variation in isozymes and morphological characters within and between seven populations from the mountains, foothills, and lowlands. Additionally, data on reproduction, habitat, germination and population sizes were collected to examine possible causes of variation and differentiation. We found highest genetic diversity (v_a, v_go) in the lowland and foothill populations, and highest genetic differentiation (D_ja, D_jgo) (i.e., lowest genetic identity: Nei's I) in the lowland populations. The low diversity of the mountain populations might indicate that they are derived from lowland populations. Surprisingly, highest genetic diversity was found in the smallest population. This indicates that in small remnant populations of taxa with a mixed mating system, genetic diversity might be maintained even after many generations after reducing population size dramatically. We found some relationship between genetic diversity and high fitness (germination success) and (inversely) with seed size. Plant size and reproductive success are negatively correlated with altitude, whereas flower size and seed size seem to be subject to other forces of selection. Combining all morphometric, reproductive and genetic traits, the lowland populations are most strongly differentiated and therefore of highest conservation priority.     

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.