Comparison of genetic diversity estimates within and among populations of maritime pine using chloroplast simple-sequence repeat and amplified fragment length polymorphism data

We compared the genetic variation of Pinus pinaster populations using amplified fragment length polymorphism (AFLP) and chloroplast simple-sequence repeat (cpSSR) loci. Populations' levels of diversity within groups were found to be similar with AFLPs, but not with cpSSRs. The high interlocus variance associated with the AFLP loci could account for the lack of differences in the former. Although AFLPs revealed much lower genetic diversity than cpSSRs, the levels of among-population differentiation found with the two types of marker were similar, provided that loci showing fewer than four null-homozygotes, in any population, were pruned from the AFLP data. Moreover, the French and Portuguese populations were clearly differentiated from each other, with both markers. The Mantel test showed that the genetic distance matrix calculated using the AFLP data was correlated with the matrix derived from the cpSSRs. Because of the concordance found between markers we conclude that gene flow was indeed the predominant force shaping nuclear and chloroplastic genetic variation of the populations within regions, at the geographical scale studied.     

Genetic variation and island biogreography: Microsatellite and mitochondrial DNA variation in island populations of the Australian bush rat, Rattus fuscipes greyii

To understand the impact of various factors on the maintenance of genetic variation in natural populations, we need to focus on situations where at least some of these factors are removed or controlled. In this study, we used highly variable, presumably neutral, microsatellite and mtDNA markers to assess the nature of genetic variation in 14 island and two mainland populations of the Australian bush rat, where there is no migration between islands. Thus we are controlling for selection and gene flow. Both marker sets revealed low levels of diversity within the small island populations and extreme differentiation between populations. For six microsatellite loci, all of the small island populations had less genetic variation than the mainland populations; reduction in allelic diversity was more pronounced than loss of heterozygosity. Kangaroo Island, the large island population, had similar levels of diversity to the mainland populations. A 442 base pair (bp) section of the mtDNA control region was screened for variation by outgroup heteroduplex analysis/temperature gradient gel electrophoresis (OHA/TGGE). Only three of the 13 small island populations showed haplotypic diversity: Gambier (2), Waldegrave (2), and Eyere (3). The level of haplotypic diversity in the small island populations was similar to that on the mainland, most likely reflecting a recent population bottleneck on the mainland. In contrast, Kangaroo Island had 9 mtDNA haplotypes. The dominant factor influencing genetic diversity on the islands was island size. No correlation was detected between genetic diversity and the time since isolation or distance form the mainland. The combination of genetic drift within and complete isolation among the small island populations has resulted in rapid and extreme population divergence. Population pair-wise comparisons of allele frequency distributions showed significant differences for all populations for all loci (F _st = 0.11–0.84, R _st = 0.07–0.99). For the mtDNA control region, 92.6% of variation was apportioned between populations; only the Pearson islands shared a haplotype. Mantel tests of pair-wise genetic distance with pair-wise geographic distance showed no significant geographical clustering of haplotypes. However, population substructuring was detected within populations where sampling was conducted over a broader geographical range, as indicated by departures from Hardy-Weinberg equilibrium. Thus substructuring in the ancestral population cannot be ruled out. The dominant evolutionary forces on the islands, after the initial founder event, are stochastic population processes such as genetic drift and mutation. This revised version was published online in July 2006 with corrections to the Cover Date.     

Genetic conservation and management of the California endemic,Torrey pine (Pinus torreyana Parry): Implications of genetic rescue in a genetically depauperate species

Rare species present a challenge under changing environmental conditions as the genetic consequences of rarity may limit species ability to adapt to environmental change. To evaluate the evolutionary potential of a rare species, we assessed variation in traits important to plant fitness using multigenerational common garden experiments. Torrey pine, Pinus torreyana Parry, is one of the rarest pines in the world, restricted to one mainland and one island population. Morphological differentiation between island and mainland populations suggests adaptation to local environments may have contributed to trait variation. The distribution of phenotypic variances within the common garden suggests distinct population‐specific growth trajectories underlay genetic differences, with the island population exhibiting substantially reduced genetic variance for growth relative to the mainland population. Furthermore, F1 hybrids, representing a cross between mainland and island trees, exhibit increased height accumulation and fecundity relative to mainland and island parents. This may indicate genetic rescue via intraspecific hybridization could provide the necessary genetic variation to persist in environments modified as a result of climate change. Long‐term common garden experiments, such as these, provide invaluable resources to assess the distribution of genetic variance that may inform conservation strategies to preserve evolutionary potential of rare species, including genetic rescue.     

Among population differentiation at nuclear genes in native Scots pine (Pinus sylvestris L.) in Scotland

In the Scottish Highlands, Scots pine is at the north-western extreme of its wide natural distribution. Here, the remaining native populations are patchily distributed in highly variable environments, from the more continental, drier eastern Highlands to the milder, wetter Atlantic Ocean coast. As these pinewoods are the remnants of a naturally established forest, they form a valuable system for analysis of genetic and adaptive variation in heterogeneous environments. Using samples from across the Scottish population, we analysed data from nuclear and mitochondrial genes to assess patterns of within and between population genetic variation. Within population diversity levels were high, and significant genetic differentiation among pairs of Scottish populations at relatively small spatial scales was present at several nuclear loci. At these loci, no differentiation had been found among continental populations, even those separated by large geographic distances. Overall, no clear clustering of Scottish samples was found in population structure analysis suggesting that geographically distant populations with high intra-population nucleotide diversity are not strongly isolated or diverged from each other. Scottish populations lacked a mitotype that is widespread in eastern and north-eastern Europe, indicating that pines from that area may not have participated in the most recent colonisation of the British Isles.     

Genetic variation and structure of house sparrow populations: is there an island effect?

Population genetic structure and intrapopulation levels of genetic variation have important implications for population dynamics and evolutionary processes. Habitat fragmentation is one of the major threats to biodiversity. It leads to smaller population sizes and reduced gene flow between populations and will thus also affect genetic structure. We use a natural system of island and mainland populations of house sparrows along the coast of Norway to characterize the different population genetic properties of fragmented populations. We genotyped 636 individuals distributed across 14 populations at 15 microsatellite loci. The level of genetic differentiation was estimated using F‐statistics and specially designed Mantel tests were conducted to study the influence of population type (i.e. mainland or island) and geographic distance on the genetic population structure. Furthermore, the effects of population type, population size and latitude on the level of genetic variation within populations were examined. Our results suggest that genetic processes on islands and mainland differed in two important ways. First, the intrapopulation level of genetic variation tended to be lower and the occurrence of population bottlenecks more frequent on islands than the mainland. Second, although the general level of genetic differentiation was low to moderate, it was higher between island populations than between mainland populations. However, differentiation increased in mainland populations somewhat faster with geographical distance. These results suggest that population bottleneck events and genetic drift have been more important in shaping the genetic composition of island populations compared with populations on the mainland. Such knowledge is relevant for a better understanding of evolutionary processes and conservation of threatened populations.     

Mitochondrial DNA variation in the European otter (Lutra lutra) and the use of spatial autocorrelation analysis in conservation

To add genetic information to the international conservation efforts on European otters Lutra lutra, we investigated the genetic population structure in and around a known "source" population of the otter, the Oberlausitz (OL) in eastern Germany. This was complemented by a first survey of genetic variation levels in the Central European otter population. Sequence analysis of 300bp of the mitochondrial control region in 76 specimens from the eastern German study region and 53 individuals from several other European populations revealed a low level of genetic variation, with only 5 haplotypes present and nucleotide diversities within populations ranging from 0.00% to 0.17%. Apart from eastern Germany, one haplotype was by far the most abundant one, from which other, only locally occurring types, could be derived by a single point mutation. This suggests a single Pleistocene refugium from which the analyzed European regions have been reinvaded after the glaciations. Within eastern Germany, two abundant haplotypes were found. Their occurrence differed significantly among subregions of eastern Germany. The uneven distribution of a locally restricted but abundant haplotype could be explained by isolation-by-distance and might reflect emigration from the OL source population to surrounding regions. This suggests that vital local populations can indeed serve as "sources" for the invasion of surrounding areas. Given a suitable genetic marker, we suggest a spatial autocorrelation analysis to monitor the genetic effect of such an emigration from a source population.     

Large mainland populations of South Island robins retain greater genetic diversity than offshore island refuges

For conservation purposes islands are considered safe refuges for many species, particularly in regions where introduced predators form a major threat to the native fauna, but island populations are also known to possess low levels of genetic diversity. The New Zealand archipelago provides an ideal system to compare genetic diversity of large mainland populations where introduced predators are common, to that of smaller offshore islands, which serve as predator-free refuges. We assessed microsatellite variation in South Island robins (Petroica australis australis), and compared large mainland, small mainland, natural island and translocated island populations. Large mainland populations exhibited more polymorphic loci and higher number of alleles than small mainland and natural island populations. Genetic variation did not differ between natural and translocated island populations, even though one of the translocated populations was established with five individuals. Hatching failure was recorded in a subset of the populations and found to be significantly higher in translocated populations than in a large mainland population. Significant population differentiation was largely based on heterogeneity in allele frequencies (including fixation of alleles), as few unique alleles were observed. This study shows that large mainland populations retain higher levels of genetic diversity than natural and translocated island populations. It highlights the importance of protecting these mainland populations and using them as a source for new translocations. In the future, these populations may become extremely valuable for species conservation if existing island populations become adversely affected by low levels of genetic variation and do not persist.     

Recurrence of the R408W Mutation in the Phenylalanine Hydroxylase Locus in Europeans

The relative frequency of the common phenylalanine hydroxylase (PAH) mutation R408W and its associations with polymorphic RFLP, VNTR, and short-tandem-repeat (STR) sites in the PAH gene were examined in many European populations and one representative North American population of defined European descent. This mutation was found to cluster in two regions: in northwest Europe among Irish and Scottish peoples, and in eastern Europe, including the Commonwealth of Independent States. This allele was significantly less frequent in intervening populations. In eastern European populations, the R408W mutation is strongly associated with RFLP haplotype 2, the three-copy VNTR allele (VNTR 3), and the 240-bp STR allele. In northwestern European populations, it is strongly associated with RFLP haplotype 1, the VNTR allele containing eight repeats (VNTR 8), and the 244-bp STR allele. An examination of the linkage between the R408W mutation and highly polymorphic RFLP, VNTR, and STR haplotypes suggests that recurrence is the most likely mechanism to account for the two different major haplotype associations of R408W in Europe.     

Genetic variation of Pinus uncinata (Pinaceae) in the Pyrenees determined with cpSSR markers

The genetic variation within and between 13 populations (385 individuals) of Pinus uncinata was analyzed with ten chloroplast microsatellite markers. Both the infinite allele mutation and stepwise mutation model (SMM) have been applied to the analysis of the genetic structure and the geographical distribution of haplotypic variation. High level of genetic diversity and low but significant differentiation among compared population were found. Three marginal populations, Sierra de Cebollera, Margaride Mountains and Sierra de Gúdar are strongly differentiated from the rest. Mutations following SMM-like process contributed significantly to the regional differentiation. The pattern of genetic structure observed in mountain pine is common in conifers with a wide distribution range. Lack of significant genetic structuring may be a result of a recent fragmentation of a historically larger population and/or interspecific hybridization and introgression. The southernmost populations from the Sierra Cebollera and the Sierra de Gúdar are the most genetically distinct. This suggests a long period of spatial isolation and/or origin from different ancestral populations.     

Genetic diversity and population size: island populations of the common shrew, Sorex araneus

Populations of many species are currently being fragmented and reduced by human interactions. These processes will tend to reduce genetic diversity within populations and reduce individual heterozygosities because of genetic drift, inbreeding and reduced migration. Conservation biologists need to know the effect of population size on genetic diversity, as this is likely to influence a population's ability to persist. Island populations represent an ideal natural experiment with which to study this problem. In a study of common shrews (Sorex araneus) on offshore Scottish islands, 497 individuals from 13 islands of different sizes and 6 regions on the mainland were trapped and genotyped at eight microsatellite loci. Previous genetic work had revealed that most of the islands in this study were highly genetically divergent from one another and the mainland. We found that most of the islands exhibited lower genetic diversity than the mainland populations. In the island populations, mean expected heterozygosity, mean observed heterozygosity and mean allelic richness were significantly positively correlated with log island size and log population size, which were estimated using habitat population density data and application of a Geographic Information System.     

Genetic variation and biogeography of Mertensia maritima (Boraginaceae)

To elucidate the population dynamics and biogeographic connections of Mertensia maritima in the North Atlantic we studied the genetic variation within and among 16 local populations in the Oslofjord region (south-east Norway), Finnmark (northern Norway), Iceland and Svalbard, using starch gel isozyme electrophoresis. We also tested the self-compatibility experimentally. Mertensia maritima was found to be self-compatible and able to self-pollinate by autodeposition. The total level of genetic variation was extremely low. All the screened individuals were homozygous at all loci. All individuals from the populations of mainland Norway and Iceland were isozymatically identical, whereas individuals from Svalbard differed from these at three loci (AAT-1, PGM-3 and SKD-2). The genetic identity between any population at Svalbard and any population on the mainland was 0.86. These genetic patterns may have resulted from a combination of inbreeding, frequent local population extinctions, and long-distance dispersal. The isozymatic distinctness of the Svalbard populations is consistent with the arctic distribution of the morphologically based var. tenella , here raised to subspecific rank as Mertensia maritima ssp. tenella.     

Microsatellite analysis of genetic structure in the mangrove species Avicennia marina (Forsk.) Vierh. (Avicenniaceae)

The level of genetic variation throughout the entire worldwide range of the mangrove species Avicennia marina (Forsk.) Vierh. was examined using microsatellite markers. Three microsatellite loci detected high levels of allelic diversity (70 alleles in total), essential for an accurate estimation of population genetic parameters. The informativeness of the microsatellite loci tended to increase with increasing average number of repeats. The levels of heterozygosity detected for each population, over all loci, ranged from 0.0 to 0.8, with an average of 0.407, indicating that some populations had little or no genetic variation, whereas others had a large amount. Populations at the extremes of the distribution range showed reduced levels of heterozygosity, and significant levels of inbreeding. This is not unexpected as these populations may be subject to founder effects and environmental constraints. The presence of genetic structure was tested in A. marina populations using three models: (i) a single panmictic model; (ii) the discrete subpopulation model; and (iii) the isolation by distance model. The discrete subpopulations model was supported by the overall measures of population differentiation based on the infinite alleles model (F-statistics), and the stepwise mutation model (R statistics). In addition, an analysis of molecular variance (AMOVA), using both theoretical models, found that most of the variation was between populations (41-71%), and within individuals in the total population (31-49%). There was little variation among individuals within populations (0-10%). There was no significant isolation by distance. The high levels of genetic differentiation observed among populations of A. marina may be due to environmental and ecological factors, particularly past sea level and climatic changes.     

Phylogeographic structure and postglacial evolutionary history of water voles (Arvicola terrestris) in the United Kingdom

The phylogeographic pattern of mitochondrial DNA variation in water voles (Arvicola terrestris) from 57 localities across the United Kingdom and representative samples from Spain, France, Switzerland and Finland was determined from sequence variation in the central portion of the control region. Twenty-seven different haplotypes were resolved which formed two distinct phylogenetic clades. This major division separated haplotypes found in Scotland from those found in England and Wales. Nested clade analysis of haplotypes indicated that such a division was a consequence of allopatric fragmentation. The haplotypes found in Switzerland, France and Spain clustered with Scottish haplotypes, whereas the haplotype from Finland clustered with the English/Welsh haplotypes. These patterns indicate that contemporary Scottish populations are derived from an Iberian glacial refugium, whereas English and Welsh populations are derived from an eastern European refugium. As such, the postglacial recolonization of the United Kingdom must have involved two colonization events, either in different localities with no subsequent contact, or as two waves separated over time, with the second wave of colonizers displacing the first. An analysis of molecular variance (AMOVA) identified significant population genetic divergence within both the major clades, indicative of restricted gene flow and regional population isolation. The implications of both phylogeographical and population genetic structure are discussed in context with the conservation of water voles in Britain.     

Distinct Phylogeographic Structures of Wild Radish (Raphanus sativus L. var. raphanistroides Makino) in Japan

Coastal plants with simple linear distribution ranges along coastlines provide a suitable system for improving our understanding of patterns of intra-specific distributional history and genetic variation. Due to the combination of high seed longevity and high dispersibility of seeds via seawater, we hypothesized that wild radish would poorly represent phylogeographic structure at the local scale. On the other hand, we also hypothesized that wild radish populations might be geographically differentiated, as has been exhibited by their considerable phenotypic variations along the islands of Japan. We conducted nuclear DNA microsatellite loci and chloroplast DNA haplotype analyses for 486 samples and 144 samples, respectively, from 18 populations to investigate the phylogeographic structure of wild radish in Japan. Cluster analysis supported the existence of differential genetic structures between the Ryukyu Islands and mainland Japan populations. A significant strong pattern of isolation by distance and significant evidence of a recent bottleneck were detected. The chloroplast marker analysis resulted in the generation of eight haplotypes, of which two haplotypes (A and B) were broadly distributed in most wild radish populations. High levels of variation in microsatellite loci were identified, whereas cpDNA displayed low levels of genetic diversity within populations. Our results indicate that the Kuroshio Current would have contributed to the sculpting of the phylogeographic structure by shaping genetic gaps between isolated populations. In addition, the Tokara Strait would have created a geographic barrier between the Ryukyu Islands and mainland Japan. Finally, extant habitat disturbances (coastal erosion), migration patterns (linear expansion), and geographic characteristics (small islands and sea currents) have influenced the expansion and historical population dynamics of wild radish. Our study is the first to record the robust phylogeographic structure in wild radish between the Ryukyu Islands and mainland Japan, and might provide new insight into the genetic differentiation of coastal plants across islands.     

Reduced MHC class II diversity in island compared to mainland populations of the black-footed rock-wallaby (<Emphasis Type="Italic">Petrogale lateralis lateralis</Emphasis>)

Many animal populations that are endangered in mainland areas exist in stable island populations, which have the potential to act as an “ark” in case of mainland population declines. Previous studies have found neutral genetic variation in such species to be up to an order of magnitude lower in island compared to mainland populations. If low genetic variation is prevalent across fitness-related loci, this would reduce the effectiveness of island populations as a source of individuals to supplement declining mainland populations or re-establish extinct mainland populations. One such species, the black-footed rock-wallaby (Petrogale lateralis lateralis), exists within fragmented mainland populations and small island populations off Western Australia. We examined sequence variation in this species within a fitness-related locus under positive selection, the MHC class II DAB β1 locus. The mainland populations displayed greater levels of allelic diversity (4–7 alleles) than the island population, despite being small and isolated, and contained at least two DAB gene copies. The island population displayed low allelic diversity (2 alleles) and fewer alleles per individual in comparison to mainland populations, and probably possesses only one DAB gene copy. The patterns of DAB diversity suggested that the island population has a markedly lower level of genetic variation than the mainland populations, in concordance with results from microsatellites (genotyped in a previous study), but preserved unique alleles which were not found in mainland populations. Where possible, conservation actions should pool individuals from multiple populations, not only island populations, for translocation programs, and focus on preventing further declines in mainland populations.     

Population structure and colonization of Bactrocera dorsalis (Diptera: Tephritidae) in China,inferred from mtDNA COI sequences

The oriental fruit fly, Bactrocera dorsalis, is a serious pest of fruits and vegetables in South‐east Asia, and, because of quarantine restrictions, impedes international trade and economic development in the region. Revealing genetic variation in oriental fruit fly populations will provide a better understanding of the colonization process and facilitate the quarantine and management of this species. The genetic structure in 15 populations of oriental fruit fly from southern China, Laos and Myanmar in South‐east Asia was examined with a 640‐bp sequence of the mitochondrial cytochrome oxidase subunit I (COI) gene. The highest levels of genetic diversity were found in Laos and Myanmar. Low to medium levels of genetic differentiation (F_ST ≤ 0.134) were observed among populations. Pooled populations from mainland China differed from those in Laos and Myanmar (F_ST = 0.024). Genetic structure across the region did not follow the isolation‐by‐distance model. The high genetic diversity observed in Laos and Myanmar supports the South‐east Asian origin of B. dorsalis. High genetic diversity and significant differentiation between some populations within mainland China indicate B. dorsalis populations have been established in the region for an extended period of time. High levels of genetic diversity observed among the five populations from Hainan Island and similarity between the Island and Chinese mainland populations indicate that B. dorsalis was introduced to Hainan from the mainland and has been on the island for many years. High genetic diversity in the recently established population in Shanghai (Pudong) suggests multiple introductions or a larger number of founders.     

Regional heterogeneity and gene flow maintain variance in a quantitative trait within populations of lodgepole pine

Genetic variation is of fundamental importance to biological evolution, yet we still know very little about how it is maintained in nature. Because many species inhabit heterogeneous environments and have pronounced local adaptations, gene flow between differently adapted populations may be a persistent source of genetic variation within populations. If this migration-selection balance is biologically important then there should be strong correlations between genetic variance within populations and the amount of heterogeneity in the environment surrounding them. Here, we use data from a long-term study of 142 populations of lodgepole pine (Pinus contorta) to compare levels of genetic variation in growth response with measures of climatic heterogeneity in the surrounding region. We find that regional heterogeneity explains at least 20% of the variation in genetic variance, suggesting that gene flow and heterogeneous selection may play an important role in maintaining the high levels of genetic variation found within natural populations.     

Phylogeography of Scots pine in Europe and Asia based on mtDNA polymorphisms

We analyzed mitochondrial DNA polymorphisms to search for evidence of the genetic structure and patterns of admixture in 124 populations (N = 1407 trees) across the distribution of Scots pine in Europe and Asia. The markers revealed only a weak population structure in Central and Eastern Europe and suggested postglacial expansion to middle and northern latitudes from multiple sources. Major mitotype variants include the remnants of Scots pine at the north-western extreme of the distribution in the Scottish Highlands; two main variants (western and central European) that contributed to the contemporary populations in Norway and Sweden; the central-eastern European variant present in the Balkan region, Finland, and Russian Karelia; and a separate one common to most eastern European parts of Russia and western Siberia. We also observe signatures of a distinct refugium located in the northern parts of the Black Sea basin that contributed to the patterns of genetic variation observed in several populations in the Balkans, Ukraine, and western Russia. Some common haplotypes of putative ancient origin were shared among distant populations from Europe and Asia, including the most southern refugial stands that did not participate in postglacial recolonization of northern latitudes. The study indicates different genetic lineages of the species in Europe and provides a set of genetic markers for its finer-scale population history and divergence inference.     

Microsatellite Variation in Honey Bee (Apis Mellifera L.) Populations: Hierarchical Genetic Structure and Test of the Infinite Allele and Stepwise Mutation Models

Samples from nine populations belonging to three African (intermissa, scutellata and capensis) and four European (mellifera, ligustica, carnica and cecropia) Apis mellifera subspecies were scored for seven microsatellite loci. A large amount of genetic variation (between seven and 30 alleles per locus) was detected. Average heterozygosity and average number of alleles were significantly higher in African than in European subspecies, in agreement with larger effective population sizes in Africa. Microsatellite analyses confirmed that A. mellifera evolved in three distinct and deeply differentiated lineages previously detected by morphological and mitochondrial DNA studies. Dendrogram analysis of workers from a given population indicated that super-sisters cluster together when using a sufficient number of microsatellite data whereas half-sisters do not. An index of classification was derived to summarize the clustering of different taxonomic levels in large phylogenetic trees based on individual genotypes. Finally, individual population X loci data were used to test the adequacy of the two alternative mutation models, the infinite allele model (IAM) and the stepwise mutation models. The better fit overall of the IAM probably results from the majority of the microsatellites used including repeats of two or three different length motifs (compound microsatellites).     

Population genetic structure and colonisation history of the tool-using New Caledonian crow

New Caledonian crows exhibit considerable variation in tool making between populations. Here, we present the first study of the species' genetic structure over its geographical distribution. We collected feathers from crows on mainland Grande Terre, the inshore island of Toupéti, and the nearby island of Maré where it is believed birds were introduced after European colonisation. We used nine microsatellite markers to establish the genotypes of 136 crows from these islands and classical population genetic tools as well as Approximate Bayesian Computations to explore the distribution of genetic diversity. We found that New Caledonian crows most likely separate into three main distinct clusters: Grande Terre, Toupéti and Maré. Furthermore, Toupéti and Maré crows represent a subset of the genetic diversity observed on Grande Terre, confirming their mainland origin. The genetic data are compatible with a colonisation of Maré taking place after European colonisation around 1900. Importantly, we observed (1) moderate, but significant, genetic differentiation across Grande Terre, and (2) that the degree of differentiation between populations on the mainland increases with geographic distance. These data indicate that despite individual crows' potential ability to disperse over large distances, most gene flow occurs over short distances. The temporal and spatial patterns described provide a basis for further hypothesis testing and investigation of the geographical variation observed in the tool skills of these crows.