Global Population Genetic Structure of Caenorhabditis remanei Reveals Incipient Speciation

Mating system transitions dramatically alter the evolutionary trajectories of genomes that can be revealed by contrasts of species with disparate modes of reproduction. For such transitions in Caenorhabditis nematodes, some major causes of genome variation in selfing species have been discerned. And yet, we have only limited understanding of species-wide population genetic processes for their outcrossing relatives, which represent the reproductive state of the progenitors of selfing species. Multilocus-multipopulation sequence polymorphism data provide a powerful means to uncover the historical demography and evolutionary processes that shape genomes. Here we survey nucleotide polymorphism across the X chromosome for three populations of the outcrossing nematode Caenorhabditis remanei and demonstrate its divergence from a fourth population describing a closely related new species from China, C. sp. 23. We find high genetic variation globally and within each local population sample. Despite geographic barriers and moderate genetic differentiation between Europe and North America, considerable gene flow connects C. remanei populations. We discovered C. sp. 23 while investigating C. remanei, observing strong genetic differentiation characteristic of reproductive isolation that was confirmed by substantial F(2) hybrid breakdown in interspecific crosses. That C. sp. 23 represents a distinct biological species provides a cautionary example of how standard practice can fail for mating tests of species identity in this group. This species pair permits full application of divergence population genetic methods to obligately outcrossing species of Caenorhabditis and also presents a new focus for interrogation of the genetics and evolution of speciation with the Caenorhabditis model system.     

Inbreeding and outbreeding depression in Stylidium hispidum: implications for mixing seed sources for ecological restoration

The benefits of composite rather than local seed provenances for ecological restoration have recently been argued, largely on the basis of maximizing evolutionary potential. However, these arguments have downplayed the potentially negative consequences of outbreeding depression once mixed provenances interbreed. In this study, we compared intraspecific F1 hybrid performance and molecular marker differentiation among four populations of Stylidium hispidum, a species endemic to Southwestern Australia. Multivariate ordination of 134 AFLP markers analyzed genetic structure and detected two clusters of paired sites that diverged significantly for marker variation along a latitudinal boundary. To test for outbreeding depression and to determine the consequences of molecular population divergence for hybrid fitness, we conducted controlled pollinations and studied germination and survival for three cross categories (within‐population crosses, short‐ and long‐distance F1 hybrids) for paired sites distributed within and between the two genetically differentiated regions. We found evidence of outbreeding depression in long‐distance hybrids (111–124 km), and inbreeding depression among progeny of within‐population crosses, relative to short‐distance (3–10 km) hybrids, suggesting an intermediate optimal outcrossing distance in this species. These results are discussed in light of the evolutionary consequences of mixing seed sources for biodiversity restoration.     

Population genetic divergence corresponds with species-level biodiversity patterns in the large genus Begonia

Begonia is one of the largest angiosperm genera, containing over 1500 species. Some aspects of the distribution of biodiversity in the genus, such as the geographical restrictions of monophyletic groups, the rarity and morphological variability of widespread species, and a preponderance of narrow endemics, suggest that restricted gene flow may have been a factor in the formation of so many species. In order to investigate whether this inference based on large-scale patterns is supported by data at the population level, we examined the distribution of genetic variation within Begonia sutherlandii in the indigenous forests of Kwazulu-Natal, South Africa, using microsatellite markers. Despite the species being predominantly outbreeding, we found high and significant levels of population structure (standardized      = 0.896). Even within individual populations, there was evidence for clear differentiation of subpopulations. There is thus congruence in evolutionary patterns ranging from interspecific phylogeny, the distribution of individual species, to the levels of population differentiation. Despite this species-rich genus showing a pan-tropical distribution, these combined observations suggest that differentiation occurs over very local scales. Although strongly selected allelic variants can maintain species cohesion with only low levels of gene flow, we hypothesize that in Begonia , gene flow levels are often so low, that divergence in allopatry is likely to be a frequent occurrence, and the lack of widespread species may in part be attributable to a lack of a mechanism for holding them together.     

Molecular evolution and quantitative variation for chemosensory behaviour in the nematode genus Caenorhabditis

Caenorhabditis elegans is a model organism in biology, yet despite the tremendous information generated from genetic, genomic and functional analyses, C. elegans has rarely been used to address questions in ecological genetics. Here, we analyse genetic variation for chemosensory behaviour, an ecologically important trait that is also genetically well characterized, at both the phenotypic and molecular levels within three species of the genus Caenorhabditis. We show that the G-protein ODR-3 plays an important role in chemosensory avoidance behaviour and identify orthologues of odr-3 in C. briggsae and C. remanei. Both quantitative genetic analysis of chemosensory behaviour and molecular population genetic analysis of odr-3 show that there is little genetic variation among a worldwide collection of isolates of the primarily selfing C. elegans, whereas there is substantially more variation within a single population of the outcrossing C. remanei. Although there are a large number of substitutions at silent sites within odr-3 among the three species, molecular evolution at the protein level is extremely conserved, suggesting that odr-3 plays an important role in cell signalling during chemosensation and/or neuronal cilia development in C. remanei and in C. briggsae as it does in C. elegans. Our results suggest that C. remanei may be a more suitable subject for ecological and evolutionary genetic studies than C. elegans.     

Specialist versus generalist life histories and nucleotide diversity in Caenorhabditis nematodes

Species with broad ecological amplitudes with respect to a key focal resource, niche generalists, should maintain larger and more connected populations than niche specialists, leading to the prediction that nucleotide diversity will be lower and more subdivided in specialists relative to their generalist relatives. This logic describes the specialist-generalist variation hypothesis (SGVH). Some outbreeding species of Caenorhabditis nematodes use a variety of invertebrate dispersal vectors and have high molecular diversity. By contrast, Caenorhabditis japonica lives in a strict association and synchronized life cycle with its dispersal host, the shield bug Parastrachia japonensis, itself a diet specialist. Here, we characterize sequence variation for 20 nuclear loci to investigate how C. japonica''s life history shapes nucleotide diversity. We find that C. japonica has more than threefold lower polymorphism than other outbreeding Caenorhabditis species, but that local populations are not genetically disconnected. Coupled with its restricted range, we propose that its specialist host association contributes to a smaller effective population size and lower genetic variation than host generalist Caenorhabditis species with outbreeding reproductive modes. A literature survey of diverse organisms provides broader support for the SGVH. These findings encourage further testing of ecological and evolutionary hypotheses with comparative population genetics in Caenorhabditis and other taxa.     

Microevolutionary studies in nematodes: a beginning

Comparisons between related species often allow the detailed genetic analysis of evolutionary processes. Here we advocate the use of the nematode Caenorhabditis elegans (and several other rhabditid species) as model systems for microevolutionary studies. Compared to Drosophila species, which have been a mainstay of such studies, C. elegans has a self-fertilizing mode of reproduction, a shorter life cycle and a convenient cell-level analysis of phenotypic variation. Data concerning its population genetics and ecology are still scarce, however. We review molecular, behavioral and developmental intraspecific polymorphisms for populations of C. elegans, Oscheius sp. 1 and Pristionchus pacificus. Focusing on vulval development, which has been well characterized in several species, we discuss relationships between patterns of variations: (1) for a given genotype (developmental variants), (2) after mutagenesis (mutability), (3) in different populations of the same species (polymorphisms) and (4) between closely related species. These studies have revealed that evolutionary variations between sister species affect those characters that show phenotypic developmental variants, that are mutable and that are polymorphic within species.     

Population genetics of Caenorhabditis elegans: the paradox of low polymorphism in a widespread species

Caenorhabditis elegans has become one of the most widely used model research organisms, yet we have little information on evolutionary processes and recent evolutionary history of this widespread species. We examined patterns of variation at 20 microsatellite loci in a sample of 23 natural isolates of C. elegans from various parts of the world. One-half of the loci were monomorphic among all strains, and overall genetic variation at microsatellite loci was low, relative to most other species. Some population structure was detected, but there was no association between the genetic and geographic distances among different natural isolates. Thus, despite the nearly worldwide occurrence of C. elegans, little evidence was found for local adaptation in strains derived from different parts of the world. The low levels of genetic variation within and among populations suggest that recent colonization and population expansion might have occurred. However, the patterns of variation are not consistent with population expansion. A possible explanation for the observed patterns is the action of background selection to reduce polymorphism, coupled with ongoing gene flow among populations worldwide.     

Reintroduction of the endangered and endemic plant species Cochlearia bavarica—Implications from conservation genetics

Population reintroduction is a common practice in conservation, but often fails, also due to the effects of inbreeding or outbreeding depression. Cochlearia bavarica is a strongly endangered plant species endemic to Bavaria in Germany, constantly declining since the late 1980s. Therefore, population reintroduction is intended. In this study, we analyzed genetic diversity within and genetic differentiation between all 32 remnant populations of the species in Swabia and Upper Bavaria using amplified fragment length polymorphisms. Our aim was to increase reintroduction success by providing data to avoid negative effects of inbreeding and outbreeding and to preserve the natural genetic pattern of the species. Genetic diversity within populations was low but similar to other rare and endemic species and varied strongly between populations but did not depend on population size. Our analysis revealed a strong geographic pattern of genetic variation. Genetic differentiation was strongest between Swabia and Upper Bavaria and at the population level, whereas differentiation between subpopulations was comparatively low. Isolation by distance and genetic differentiation was stronger among populations from Upper Bavaria than from Swabia. From the results of our study, we derived recommendations for a successful reintroduction of the species. We suggest using rather genetically variable than large populations as reintroduction sources. Moreover, the exchange of plant material between Swabia and Upper Bavaria should be completely avoided. Within these regions, plant material from genetically similar populations should preferably be used for reintroduction, whereas the exchange among subpopulations seems to be possible without a negative impact on genetic variation due to natural gene flow.     

Heterosis and outbreeding depression in crosses between natural populations of Arabidopsis thaliana

Understanding the causes and architecture of genetic differentiation between natural populations is of central importance in evolutionary biology. Crosses between natural populations can result in heterosis if recessive or nearly recessive deleterious mutations have become fixed within populations because of genetic drift. Divergence between populations can also result in outbreeding depression because of genetic incompatibilities. The net fitness consequences of between-population crosses will be a balance between heterosis and outbreeding depression. We estimated the magnitude of heterosis and outbreeding depression in the highly selfing model plant Arabidopsis thaliana, by crossing replicate line pairs from two sets of natural populations (C↔R, B↔S) separated by similar geographic distances (Italy↔Sweden). We examined the contribution of different modes of gene action to overall differences in estimates of lifetime fitness and fitness components using joint scaling tests with parental, reciprocal F₁ and F₂, and backcross lines. One of these population pairs (C↔R) was previously demonstrated to be locally adapted, but locally maladaptive quantitative trait loci were also found, suggesting a role for genetic drift in shaping adaptive variation. We found markedly different genetic architectures for fitness and fitness components in the two sets of populations. In one (C↔R), there were consistently positive effects of dominance, indicating the masking of recessive or nearly recessive deleterious mutations that had become fixed by genetic drift. The other set (B↔S) exhibited outbreeding depression because of negative dominance effects. Additional studies are needed to explore the molecular genetic basis of heterosis and outbreeding depression, and how their magnitudes vary across environments.     

Exploring the impact of neutral evolution on intrapopulation genetic differentiation in functional traits in a long-lived plant

Most plant species, particularly long-lived plants, harbor a large amount of genetic variation within populations. A central issue in evolutionary ecology is to explore levels of genetic variation and understand the mechanisms that influence them. In this study, our goals were to examine the impact of neutral evolutionary processes on the genetic variance and functional diversity within three populations of a long-lived plant (Quercus suber L.). For this purpose, we genotyped the progeny of 45 open-pollinated mother trees from three populations originating from Spain, Portugal, and Morocco using six microsatellite markers. Seedlings were planted in a common garden trial and were phenotypically characterized by seven leaf functional traits. Molecular analyses revealed weak genetic differences between Iberian and Moroccan populations. Nevertheless, high genetic differentiation was observed among maternal families within populations. Differentiation between particular maternal families from the same population reached values of 29.2 %, which far exceeds the values reported between the most genetically distant populations for this species (11.7 %). Maternal families differed also in phenology, leaf size, and shape traits. In the Moroccan population, there were correlations among matrices of distances for molecular markers, leaf shape traits (e.g., leaf circularity index), and phenology, indicating that maternal families with contrasting phenologies were genetically and functionally distinct. This, together with the moderate heritability for phenology in Moroccan population, suggests that besides selective forces, neutral evolutionary processes have promoted intrapopulation genetic divergence and contribute to maintain high levels of genetic variation within this population. Overall, our results reinforce the importance of intrapopulation studies in long-lived plants under an evolutionary context.     

THE QUANTITATIVE AND MOLECULAR GENETIC ARCHITECTURE OF A SUBDIVIDED SPECIES

In an effort to elucidate the evolutionary mechanisms that determine the genetic architecture of a species, we have analyzed 17 populations of the microcrustacean Daphnia pulex for levels of genetic variation at the level of life-history characters and molecular markers in the nuclear and mitochondrial genomes. This species is highly subdivided, with approximately 30% of the variation for nuclear molecular markers and 50% of the variation for mitochondrial markers being distributed among populations. The average level of genetic subdivision for quantitative traits is essentially the same as that for nuclear markers, which superficially suggests that the life-history characters are diverging at the neutral rate. However, the existence of a strong correlation between the levels of population subdivision and broadsense heritabilities of individual traits argues against this interpretation, suggesting instead that the among-population divergence of some quantitative traits (most notably body size) is being driven by local adaptation to different environments. The fact that the mean phenotypes of the individual populations are also strongly correlated with local levels of homozygosity indicates that variation in local inbreeding plays a role in population differentiation. Rather than being a passive consequence of local founder effects, levels of homozygosity may be selected for directly for their effects on the phenotype (adaptive inbreeding depression). There is no relationship between the levels of variation within populations for molecular markers and quantitative characters, and this is explained by the fact that the average standing genetic variation for life-history characters in this species is equivalent to only 33 generations of variation generated by mutation.     

Oscheius tipulae, a widespread hermaphroditic soil nematode, displays a higher genetic diversity and geographical structure than Caenorhabditis elegans

The nematode Oscheius tipulae belongs to the same family (Rhabditidae) as the model species Caenorhabditis elegans . Both species reproduce through self-fertilizing hermaphrodites and facultative males. Recent studies have shown that the self-fertile C. elegans and C. briggsae displayed a 20-fold lower genetic diversity than the male–female species C. remanei . Several explanations have been put forward to account for this difference, including their mode of reproduction and dynamic population structure. Here, we present the results of extensive worldwide sampling of O. tipulae , which we previously used as a laboratory organism for developmental genetics. We found that O. tipulae is much more widespread and common in soil throughout the world than Caenorhabditis species. We analysed 63 O. tipulae isolates from several continents using amplified fragment length polymorphism (AFLP). We found that O. tipulae harbours a 5-fold higher genetic diversity than C. elegans and C. briggsae . As in C. elegans , a high proportion of this diversity was found locally. Yet, we detected significant geographical differentiation, both at the worldwide scale with a latitudinal structure and between three localities in France. In summary, O. tipulae exhibited significantly higher levels of genetic diversity and large-scale geographical structure than C. elegans , despite their shared mode of reproduction. This species difference in genetic diversity may be explained by a number of other differences, such as population size, distribution, migration and dynamics. Due to its widespread occurrence and relatively high genetic diversity, O. tipulae may be a promising study species for evolutionary studies.     

Spatial and ecological population genetic structures within two island‐endemic Aeonium species of different niche width

The Crassulacean genus Aeonium is a well‐known example for plant species radiation on oceanic archipelagos. However, while allopatric speciation among islands is documented for this genus, the role of intra‐island speciation due to population divergence by topographical isolation or ecological heterogeneity has not yet been addressed. The aim of this study was to investigate intraspecific genetic structures and to identify spatial and ecological drivers of genetic population differentiation on the island scale. We analyzed inter simple sequence repeat variation within two island‐endemic Aeonium species of La Palma: one widespread generalist that covers a large variety of different habitat types (Ae. davidbramwellii) and one narrow ecological specialist (Ae. nobile), in order to assess evolutionary potentials on this island. Gene pool differentiation and genetic diversity patterns were associated with major landscape structures in both species, with phylogeographic implications. However, overall levels of genetic differentiation were low. For the generalist species, outlier loci detection and loci–environment correlation approaches indicated moderate signatures of divergent selection pressures linked to temperature and precipitation variables, while the specialist species missed such patterns. Our data point to incipient differentiation among populations, emphasizing that ecological heterogeneity and topographical structuring within the small scales of an island can foster evolutionary processes. Very likely, such processes have contributed to the radiation of Aeonium on the Canary Islands. There is also support for different evolutionary mechanisms between generalist and specialist species.     

Caenorhabditis evolution: if they all look alike, you aren't looking hard enough

Caenorhabditis elegans is widely known as a model organism for cell, molecular, developmental and neural biology, but it is also being used for evolutionary studies. A recent meeting of researchers in Portugal covered topics as diverse as phylogenetics, genetic mapping of quantitative and qualitative intraspecific variation, evolutionary developmental biology and population genetics. Here, we summarize the main findings of the meeting, which marks the formal birth of a research community dedicated to Caenorhabditis species evolution.     

OUTBREEDING DEPRESSION WITH LOW GENETIC VARIATION IN SELFING CAENORHABDITIS NEMATODES

Theory and empirical study produce clear links between mating system evolution and inbreeding depression. The connections between mating systems and outbreeding depression, whereby fitness is reduced in crosses of less related individuals, however, are less well defined. Here we investigate inbreeding and outbreeding depression in self‐fertile androdioecious nematodes, focusing on Caenorhabditis sp. 11. We quantify nucleotide polymorphism for nine nuclear loci for strains throughout its tropical range, and find some evidence of genetic differentiation despite the lowest sequence diversity observed in this genus. Controlled crosses between strains from geographically separated regions show strong outbreeding depression, with reproductive output of F1s reduced by 36% on average. Outbreeding depression is therefore common in self‐fertilizing Caenorhabditis species, each of which evolved androdioecious selfing hermaphroditism independently, but appears strongest in C. sp. 11. Moreover, the poor mating efficiency of androdioecious males extends to C. sp. 11. We propose that self‐fertilization is a key driver of outbreeding depression, but that it need not evolve as a direct result of local adaptation per se. Our verbal model of this process highlights the need for formal theory, and C. sp. 11 provides a convenient system for testing the genetic mechanisms that cause outbreeding depression, negative epistasis, and incipient speciation.     

Polymorphism and evolution of vulval precursor cell lineages within two nematode genera, Caenorhabditis and Oscheius

BACKGROUND: The cell lineage of nematodes is mostly invariant for a given species, but varies between species. One can thus wonder how a cell lineage varies during evolution. We have started a microevolutionary approach within two genera by observing lineage variations of vulval precursor cells in different natural nematode populations of the same and closely related species. RESULTS: In Caenorhabditis elegans, the P3.p cell lineage is variable within a genetically homogeneous population and polymorphic between wild strains. Irrespective of its division pattern, P3.p is competent to form vulval tissue in different C. elegans strains, whereas it is not competent in C. briggsae. In Oscheius sp. 1, P4.p and P8.p lineages are strongly polymorphic. Within each genus, these intraspecies polymorphisms in cell lineages are amplified between closely related species. In Oscheius sp. 1, the large polymorphisms in P4.p and P8.p lineages allowed us to undertake a genetic analysis of the variation between two pairs of strains. Multiple loci are involved in cell lineage differences, and variation at one locus appears to have a relatively strong effect. In addition to these large lineage variations in cells that do not normally contribute to the vulva, we find minor variations (errors) in vulval lineages, which represent the precision level of the vulval-patterning process and point to a selection pressure for maintenance of a large vulval equivalence group. CONCLUSIONS: Polymorphisms in vulval cell lineage are found within a given nematode species, and could be instrumental in explaining evolutionary variations between closely related species.     

Herbivore host-associated genetic differentiation depends on the scale of plant genetic variation examined

Herbivore adaptation to plant genetic variation can lead to reproductive isolation and the formation of host-associated lineages (host-associated differentiation, or HAD). Plant genetic variation exists along a scale, ranging from variation among individual plant genotypes to variation among plant species. Along this scale, herbivores may adapt and diverge at any level, yet few studies have examined whether herbivore differentiation exhibits scaling with respect to host variation (e.g., from genotypes to species). Determining at which level(s) herbivore differentiation occurs can provide insight into the importance of plant genetic variation on herbivore evolution. Previous studies have found strong genetic differentiation in the eriophyid mite, Aceria parapopuli, between hybrid Populus hosts and parental Populus species, but minimal neutral-locus differentiation among individual trees of the same species. We tested whether genetic differentiation in A. parapopuli scales with genetic variation in its Populus hosts. Using mite ITS1 sequence data collected among host species and among host populations, two key patterns emerged. (1) We found strong differentiation of A. parapopuli among Populus species, supporting the hypothesis that plant species differences drive reproductive isolation and HAD. (2) We did not find evidence of host-driven genetic differentiation in mites at the level of plant populations, suggesting that this level of plant variation is insufficiently strong to drive differentiation at a neutral locus. In combination with previous studies, we found that HAD occurs at the higher levels of plant genetic variation, but not at lower levels, and conclude that HAD depends on the scale of plant genetic variation examined.     

Extensive clonality and strong differentiation in the insular pacific tree Santalum insulare: implications for its conservation

BACKGROUND AND AIMS: The impact of evolutionary forces on insular systems is particularly exacerbated by the remoteness of islands, strong founder effects, small population size and the influence of biotic and abiotic factors. Patterns of molecular diversity were analysed in an island system with Santalum insulare, a sandalwood species endemic to eastern Polynesia. The aims were to evaluate clonality and to study the genetic diversity and structure of this species, in order to understand the evolutionary process and to define a conservation strategy. METHODS: Eight nuclear microsatellites were used to investigate clonality, genetic variation and structure of the French Polynesian sandalwood populations found on ten islands distributed over three archipelagos. KEY RESULTS: It was found that 58 % of the 384 trees analysed were clones. The real size of the populations is thus dramatically reduced, with sometimes only one genet producing ramets by root suckering. The diversity parameters were low for islands (n(A) = 1.5-5.0; H(E) = 0.28-0.49). No departure from Hardy-Weinberg proportion was observed except within Tahiti island, where a significant excess of homozygotes was noted in the highland population. Genetic structure was characterized by high levels of differentiation between archipelagos (27 % of the total variation) and islands (F(ST) = 0.50). The neighbour-joining tree did not discriminate the three archipelagos but separated the Society archipelago from the other two. CONCLUSIONS: This study shows that clonality is a frequent phenomenon in S. insulare. The genetic diversity within populations is lower than the values assessed in species distributed on the mainland, as a consequence of insularity. But this can also be explained by the overexploitation of sandalwood. The differentiation between archipelagos and islands within archipelagos is very high because of the limited gene flow due to oceanic barriers. Delineation of evolutionary significant units and principles for population management are proposed based on this molecular analysis.     

Comparative Study of Seed Germination and Genetic Variation of Rare Endemic Rhinanthus osiliensis and Related Widespread Congener R. rumelicus (Orobanchaceae)

Genetic diversity is affected by breeding systems, life history traits and the distribution of species. Generally, inbred species with restricted geographic distribution have lower genetic variation compared to common outbreeding plant species. Rhinanthus osiliensis is a narrow endemic growing in calcareous spring fens on the island Saaremaa, Estonia. Presumably the closest congener of R. osiliensis is R. rumelicus, which is widely distributed in Eastern and Central Europe and the Balkan Peninsula. Isozymes were used to describe the mating system, levels of genetic variation and differentiation between the populations of R. osiliensis and R. rumelicus. Genetic diversity was found to be higher in widespread R. rumelicus compared to the endemic R. osiliensis. A significantly higher inbreeding coefficient with a low germination rate indicated substantial autogamy in R. osiliensis, whereas R. rumelicus displayed a mixed mating mode with a moderate inbreeding coefficient and a high germination rate. The low genetic differentiation between R. osiliensis and R. rumelicus confirmed their close affinity and the short evolutionary period of both taxa after the postglacial colonization. Considering the restricted geographic range and inbreeding in populations of R. osiliensis, this species may face a loss of genetic variation in the future.     

Partitioning of genetic variation in the weedy clonal herb Aegopodium podagraria (Apiaceae) in Sweden

Although sexual reproduction appears to be rare in many clonal plant species, the majority of clonal species show similar levels of genetic variation to non-clonal plant species. Aegopodium podagraria is a noxious garden weed and has become a successful invader in many natural habitats. Earlier studies have demonstrated population differentiation in life-history traits between different A. podagraria populations in Sweden. In this study, we used three methods to assess genetic variation in this species. Using analyses of molecular markers (chloroplast DNA PCR-RFLP and allozyme electrophoresis) we did not detect any variation between different A. podagraria clones. However, a multivariate analysis of leaflet shape in five populations from central Sweden revealed considerable variation both within and between populations. The variation found in leaflet shape is suggested to be genetically based because the leaves were collected from plants grown from seed under similar conditions in a common garden environment. These relatively high levels of variation within and between populations indicate that there is likely to be repeated seedling recruitment within established populations. ¹