Molecular phylogenetics at the population/species interface in cave spiders of the southern Appalachians (Araneae:Nesticidae:Nesticus)

This paper focuses on the relationship between population genetic structure and speciation mechanisms in a monophyletic species group of Appalachian cave spiders (Nesticus). Using mtDNA sequence data gathered from 256 individuals, I analyzed patterns of genetic variation within and between populations for three pairs of closely related sister species. Each sister-pair comparison involves taxa with differing distributional and ecological attributes; if these ecological attributes are reflected in basic demographic differences, then speciation might proceed differently across these sister taxa comparisons. Both frequency-based and gene tree analyses reveal that the genetic structure of the Nesticus species studied is characterized by similar and essentially complete population subdivision, regardless of differences in general ecology. These findings contrast with results of prior genetic studies of cave-dwelling arthropods that have typically revealed variation in population structure corresponding to differences in general ecology. Species fragmentation through both extrinsic and intrinsic evolutionary forces has resulted in discrete, perhaps independent, populations within morphologically defined species. Large sequence divergence values observed between populations suggest that this independence may extend well into the past. These patterns of mtDNA genealogical structure and divergence imply that species as morphological lineages are currently more inclusive than basal evolutionary or phylogenetic units, a suggestion that has important implications for the study of speciation mechanisms.      

Do assemblages of Coregonus (Teleostei: Salmoniformes) in the Central Alpine region of Europe represent species flocks?

To examine models of evolution for Coregonus from the Central Alpine region of Europe, 20 populations from nine lakes were assessed for variation at six microsatellite DNA loci. Patterns of variation were tested against three evolutionary models: phenotypic plasticity, multiple invasions of lakes by divergent forms, and within-lake radiation of species flocks. All sympatric and all but one allopatric pairs of populations were significantly divergent in allele frequencies. Pairwise F -statistics indicated reduced gene flow among phenotypically divergent sympatric populations. These results reject the hypothesis that within-lake morphological and ecological diversity reflects phenotypic plasticity within a single gene pool. Genetic similarity was higher among forms within lakes than between populations of the same form in different lakes. Among-lake divergence was primarily a product of allele size differences. Mantel tests contrasting patterns of genetic divergence against patterns predicted from the multiple invasions and species flocks models indicated that the latter is the best explanation of the observed genetic variation. Thus, reproductively isolated species diverged within lakes, with similar patterns repeatedly emerging among lakes. While this study argues for a particular mode of evolution in Central Alpine Coregonus , the taxonomy of these forms remains unresolved.     

CLADES: A classification‐based machine learning method for species delimitation from population genetic data

Species are considered to be the basic unit of ecological and evolutionary studies. As multilocus genomic data are increasingly available, there have been considerable interests in the use of DNA sequence data to delimit species. In this study, we show that machine learning can be used for species delimitation. Our method treats the species delimitation problem as a classification problem for identifying the category of a new observation on the basis of training data. Extensive simulation is first conducted over a broad range of evolutionary parameters for training purposes. Each pair of known populations is combined to form training samples with a label of “same species” or “different species”. We use support vector machine (SVM) to train a classifier using a set of summary statistics computed from training samples as features. The trained classifier can classify a test sample to two outcomes: “same species” or “different species”. Given multilocus genomic data of multiple related organisms or populations, our method (called CLADES) performs species delimitation by first classifying pairs of populations. CLADES then delimits species by maximizing the likelihood of species assignment for multiple populations. CLADES is evaluated through extensive simulation and also tested on real genetic data. We show that CLADES is both accurate and efficient for species delimitation when compared with existing methods. CLADES can be useful especially when existing methods have difficulty in delimitation, for example with short species divergence time and gene flow.     

Speciation and dispersal along continental coastlines and island arcs in the Indo-West Pacific turbinid gastropod genus Lunella

Species trees were produced for the Indo-West Pacific (IWP) gastropod genus Lunella using MrBayes, BEAST, and *BEAST with sequence data from four genes. Three fossil records were used to calibrate a molecular clock. Eight cryptic species were recognized using statistical methods for species delimitation in combination with morphological differences. However, our results suggest caution in interpreting ESUs defined solely by the general mixed Yule Coalescent model in genera like Lunella, with lower dispersal abilities. Four almost entirely allopatric species groups were recovered that differ in ecology and distribution. Three groups occur predominantly along continental coastlines and one occurs on island arrays. Sympatric species occur only in the torquata and coronata groups along coastlines, whereas species in the cinerea group, distributed in two-dimensional island arrays, occur in complete allopatry. Dispersal along island arcs has been important in the maintenance of species distributions and gene flow among populations in the cinerea group. The emergence of new islands and their eventual subsidence over geological time has had important consequences for the isolation of populations and the eventual rise of new species in Lunella.     

The subspecies concept in butterflies: has its application in taxonomy and conservation biology outlived its usefulness?

Subspecies lie at the interface between systematics and population genetics, and represent a unit of biological organization in zoology that is widely used in the disciplines of taxonomy and conservation biology. In this review, we explore the utility of subspecies in relation to their application in systematics and biodiversity conservation, and briefly summarize species concepts and criteria for their diagnosis, particularly from an invertebrate perspective. The subspecies concept was originally conceived as a formal means of documenting geographical variation within species based on morphological characters; however, the utility of subspecies is hampered by inconsistencies by which they are defined conceptually, a lack of objective criteria or properties that serve to delimit their boundaries, and their frequent failure to reflect distinct evolutionary units according to population genetic structure. Moreover, the concept has been applied to populations largely comprising different components of genetic diversity reflecting contrasting evolutionary processes. We recommend that, under the general lineage (unified) species concept, the definition of subspecies be restricted to extant animal groups that comprise evolving populations representing partially isolated lineages of a species that are allopatric, phenotypically distinct, and have at least one fixed diagnosable character state, and that these character differences are (or are assumed to be) correlated with evolutionary independence according to population genetic structure. Phenotypic character types include colour pattern, morphology, and behaviour or ecology. Under these criteria, allopatric subspecies are a type of evolutionarily significant unit within species in that they show both neutral divergence through the effects of genetic drift and adaptive divergence under natural selection, and provide an historical context for identifying biodiversity units for conservation. Conservation of the adaptedness and adaptability of gene pools, however, may require additional approaches. Recent studies of Australian butterflies exemplify these points. © 2012 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, ?? , ??–??.     

Phylogenetic relationships among holarctic populations of Chironomus entis and Chironomus plumosus in view of possible horizontal transfer of mitochondrial genes

In eight Holarctic populations of two typical chironomid sibling species of the plumosus group, Chrionomus entis and Chironomus plumosus, nucleotides sequences of mitochondrial (cytb) and nuclear (gb2b) gene regions were examined. The phylogenetic trees reflecting the evolutionary histories of the nuclear and mitochondrial markers exhibited significant differences. On the tree based on the nuclear gene sequences the populations clustered according to their species affiliation, whereas on the tree based on the mitochondrial gene sequences the populations were grouped according to their geographic position. This discrepancy is probably explained by mitochondrial gene flow between sympatric species with incomplete reproductive isolation (sibling species). Based on our results together with the earlier data on nuclear and mitochondrial gene sequences of some other species from the phylogenetic group plumosus, a scheme of phylogenetic relationships within this group is proposed. This scheme is in many ways different from the traditional view on the evolutionary relationships among species of the plumosus group.     

Phylogeography and Population Genetic Structure of the Ornate Dragon Lizard, Ctenophorus ornatus

Species inhabiting ancient, geologically stable landscapes that have been impacted by agriculture and urbanisation are expected to have complex patterns of genetic subdivision due to the influence of both historical and contemporary gene flow. Here, we investigate genetic differences among populations of the granite outcrop-dwelling lizard Ctenophorus ornatus, a phenotypically variable species with a wide geographical distribution across the south-west of Western Australia. Phylogenetic analysis of mitochondrial DNA sequence data revealed two distinct evolutionary lineages that have been isolated for more than four million years within the C. ornatus complex. This evolutionary split is associated with a change in dorsal colouration of the lizards from deep brown or black to reddish-pink. In addition, analysis of microsatellite data revealed high levels of genetic structuring within each lineage, as well as strong isolation by distance at multiple spatial scales. Among the 50 outcrop populations’ analysed, non-hierarchical Bayesian clustering analysis revealed the presence of 23 distinct genetic groups, with outcrop populations less than 4 km apart usually forming a single genetic group. When a hierarchical analysis was carried out, almost every outcrop was assigned to a different genetic group. Our results show there are multiple levels of genetic structuring in C. ornatus, reflecting the influence of both historical and contemporary evolutionary processes. They also highlight the need to recognise the presence of two evolutionarily distinct lineages when making conservation management decisions on this species.     

Does genetic variation and gene flow vary with rarity in obligate seeding Persoonia species (Proteaceae)?

Theory predicts that genetic variation is a determinant of persistence, and that the abundance and distribution of variation is strongly dependent on genetic drift and gene flow. Small, isolated populations are expected to be less diverse and more differentiated than large, inter-connected populations. Thus rare species may be more at risk of extinction. We used 389 putative AFLP loci to compare genetic variation and structuring in two pairs of closely-related common (large populations geographically widespread) and rare (small populations spatially restricted) Persoonia species. We genotyped 15–22 adult plants, from four populations, covering the geographic range of each species. Although genetic diversity was low for all four species (for long-lived outcrossing perennials), we found significantly more diversity within populations of the rare species than within those of the common species. AMOVA revealed significant levels of structure both among species (21%) and populations (15%). The proportion of inter-population variation within species did not vary consistently with rarity (Pair 1 rare 21.1% versus common 16.5%; Pair 2 rare 15.8% versus common 20.6%). However populations of the rare species were more differentiated than common species with similar geographic separation, suggesting greater gene flow between populations of the common species. Therefore, even relatively small genetically isolated populations of rare Persoonia species were more diverse than large populations of common Persoonia species. We hypothesise that common Persoonia species have undergone a rapid range expansion from a narrow gene pool, while genetic diversity is maintained in the soil seed-bank of rare remnants.     

Butterfly populations in two forest fragments at the Kenya coast

Species richness, diversity and composition of butterflies in two Kenya coastal forest remnants, Muhaka and Mrima hill, were investigated. Sixty‐three species were recorded from each forest remnant from a total of 1329 individuals. Species accumulation curves for both forests did not reach an asymptote. High species similarity was recorded between the forest interior and the surrounding matrix, primarily due to invasion of the forest interior clearings by the savanna species. Despite their small sizes, these forest remnants were found to maintain viable populations of true forest butterflies. However, the number of species was less than half that recorded from the larger forest reserve of Arabuko‐Sokoke, located in the same geographical area. Records from Muhaka forest show species unique to it, not found in the larger forest reserves, underscoring the importance of small remnants in the preservation of forest biodiversity. The high species similarity between the forest remnants implied that if habitat corridors were created, gene flow between these remnants and other larger forest reserves would be possible. This would reduce the isolation of true forest butterfly populations within the remnants and potential local extinction.     

Evolutionary history of an adaptive radiation in species pairs of threespine sticklebacks (Gasterosteus): insights from mitochondrial DNA

Species pairs of threespine stickleback, Gasterosteus aculeatus, co-exist in several lakes in the Strait of Georgia, southwestern British Columbia. One species, ‘benthics’ is robust-bodied and is morphologically and behaviourally specialized for benthivory. The other species, ‘limnetics’ is specialized for planktivory in open-water habitats of the lakes. We examined mitochondrial DNA restriction site variation in benthic and limnetic sticklebacks as well as in solitary freshwater, anadromous (sea-run), and marine populations to test: (i) if benthic and limnetic pairs have evolved only once or multiple times (parallel evolution) and (ii) if the species have evolved sympatrically, or allopatrically from ‘double invasions’ of lakes by ancestral anadromous/marine sticklebacks. Stickleback mtDNA comprised a single clade with a low (mean = 0.40%) degree of sequence divergence among the 77 haplotypes resolved. Most nucleotide diversity (97%) was found within (rather than among) populations of anadromous/marine sticklebacks whereas most diversity (77%) was found among populations in freshwater sticklebacks. Significant differences in haplotype frequencies were found between benthics and limnetics in three of the four species pair lakes examined, but in all cases the pairs within lakes were characterized by unique assemblages of closely related haplotypes. Hierarchical clustering of divergence estimates suggested that comparable species from different lakes have originated independently in all lakes because in no case did comparable species from different lakes cluster together. Divergent species within lakes tended to be more closely related to one another than to species in other lakes and there were two cases were benthics and limnetics within a particular lake were monophyletic. In two of the four two-species lakes, limnetics were less divergent from putative ancestral anadromous/marine stickleback as predicted by the double invasion hypothesis, but in the two other lakes benthics were less divergent. Our data argue strongly that the species pairs have evolved independently in each lake were they now co-exist. Further, in two lakes our data are consistent with the species having evolved by sympatric divergence, but allopatric divergence followed by introgression of mtDNA that has obscured ancestral relationships cannot be discounted completely. Finally, despite remaining uncertainty about the geography of speciation, the species appear to have evolved in the face of gene flow arguing that natural selection acting on trophic ecology has been a major component of ecological speciation in sticklebacks.     

A MULTILOCUS TEST OF SIMULTANEOUS DIVERGENCE ACROSS THE ISTHMUS OF PANAMA USING SNAPPING SHRIMP IN THE GENUS ALPHEUS

The completion of the Panamanian Isthmus is one of the greatest natural experiments in evolution, sending multiple species pairs from a broad range of taxonomic groups on independent evolutionary trajectories. The resulting transisthmian sister species have been used as model systems for examining consequences that accompany cessation of gene flow in formerly panmictic populations. However, variance in pairwise genetic distances of these "geminates" often exceeds expectations, seemingly conflicting with the assumption that separation of populations was contemporaneous with the final closure of the Isthmus. Multilocus datasets and coalescent-based analytical methods can be used to estimate divergence times while accounting for variance in gene divergence that predates isolation, thus removing the need to invoke unequal divergence times. Here we present results from Bayesian analyses of sequence data from seven nuclear and one mitochondrial marker in eight transisthmian species pairs in the snapping shrimp genus Alpheus . Divergence times in two species pairs were shown to occur much earlier than the Isthmus final closure, but much of the variance in pairwise genetic distances from cytochrome oxidase I (COI) was explained when ancestral polymorphisms were accounted for. Results illustrate how coalescent approaches may be more appropriate for dating recent divergences than for estimating ancient speciation events.     

Directional genetic differentiation and relative migration

Understanding the population structure and patterns of gene flow within species is of fundamental importance to the study of evolution. In the fields of population and evolutionary genetics, measures of genetic differentiation are commonly used to gather this information. One potential caveat is that these measures assume gene flow to be symmetric. However, asymmetric gene flow is common in nature, especially in systems driven by physical processes such as wind or water currents. As information about levels of asymmetric gene flow among populations is essential for the correct interpretation of the distribution of contemporary genetic diversity within species, this should not be overlooked. To obtain information on asymmetric migration patterns from genetic data, complex models based on maximum‐likelihood or Bayesian approaches generally need to be employed, often at great computational cost. Here, a new simpler and more efficient approach for understanding gene flow patterns is presented. This approach allows the estimation of directional components of genetic divergence between pairs of populations at low computational effort, using any of the classical or modern measures of genetic differentiation. These directional measures of genetic differentiation can further be used to calculate directional relative migration and to detect asymmetries in gene flow patterns. This can be done in a user‐friendly web application called divMigrate‐online introduced in this study. Using simulated data sets with known gene flow regimes, we demonstrate that the method is capable of resolving complex migration patterns under a range of study designs.     

Habitat type as a determinant of species range sizes: the example of lotic–lentic differences in aquatic Goleoptera

Species differ in the size of their geographical ranges, but it is unclear how this is affected by the intrinsic properties of various habitat types. Using data on range sizes for 490 species of aquatic Coleoptera from the Iberian Peninsula we show that running-water (lotic) species have much smaller distributional ranges than those occurring in standing water (lentic). This robust association of habitat type and range size has independently arisen in at least four monophyletic coleopteran lineages, in Hydradephaga, Hydrophiloidea, Hydraenidae and Byrrhoidea, and several more times within these main groups. We propose that this pattern is due to different evolutionary dynamics of both habitat types: stagnant water bodies are more likely to completely disappear, requiring frequent migration of resident populations. Rivers and streams, on the contrary, have more temporal and spatial continuity, and therefore permit the long-term persistence of local populations. In less permanent habitats species will require a greater geographical mobility, which indirectly results in a larger size range. The less dispersive populations of running water should also have reduced gene flow, increasing the probability of allopatric speciation, and thus reducing the average range of more widespread ancestral species. These differences in population parameters, and the frequency of transitions between the two habitat types, may have strong macroevolutionary consequences, in particular regarding speciation rates and possible morphological specializations.     

How tree species fill geographic and ecological space in eastern North America

Background and Aims Ecologists broadly accept that the number of species present within a region balances regional processes of immigration and speciation against competitive and other interactions between populations that limit distribution and constrain diversity. Although ecological theory has, for a long time, addressed the premise that ecological space can be filled to ‘capacity’ with species, only with the availability of time-calibrated phylogenies has it been possible to test the hypothesis that diversification slows as the number of species in a region increases. Focusing on the deciduous trees of eastern North America, this study tested predictions from competition theory concerning the distribution and abundance of species.Methods Local assemblages of trees tabulated in a previous study published in 1950 were analysed. Assemblages were ordinated with respect to species composition by non-metric multidimensional scaling (NMS). Distributions of trees were analysed by taxonomically nested analysis of variance, discriminant analysis based on NMS scores, and canonical correlation analysis of NMS scores and Bioclim climate variables.Key Results Most of the variance in species abundance and distribution was concentrated among closely related (i.e. congeneric) species, indicating evolutionary lability. Species distribution and abundance were unrelated to the number of close relatives, suggesting that competitive effects are diffuse. Distances between pairs of congeneric species in NMS space did not differ significantly from distances between more distantly related species, in contrast to the predictions of both competitive habitat partitioning and ecological sorting of species.Conclusions Eastern deciduous forests of North America do not appear to be saturated with species. The distributions and abundances of individual species provide little evidence of being shaped by competition from related (i.e. ecologically similar) species and, by inference, that diversification is constrained by interspecific competition.     

Formation of a fluvial non-parasitic population of Lethenteron camtschaticum as the first step in petromyzontid speciation

To elucidate the petromyzontid speciation process, the genetic independence of the fluvial non-parasitic populations within the anadromous parasitic Lethenteron camtschaticum was estimated by using polymorphic microsatellite loci. Abundant gene flow was revealed in multitemporal scales between potentially sympatric populations, suggesting ongoing gene flow resulting from imperfect size-assortative mating between them and plastic determination of life histories. On the contrary, landlocked fluvial non-parasitic populations in the upper region of dams were genetically divergent from anadromous parasitic populations. The temporal heterogeneity of gene flow, i.e. contemporary little gene flow but significant gene flow over the long-term between the landlocked fluvial non-parasitic and anadromous parasitic populations was elucidated. In addition, the divergence time of isolation of the landlocked populations from the ancestral anadromous parasitic population was estimated to have occurred 17.9-428.2 years ago, which includes the construction times of an initial dam c. 90 years ago. These instances indicate that the landlocked populations should have very recently been established, and subsequent accumulation of divergence and development of reproductive isolation are predicted. The present landlocked fluvial non-parasitic populations should be analogous to the founder populations in terms of petromyzontid speciation. The data also strongly support the hypothesis of multitemporal and multispatial speciation in the petromyzontid stem-satellite species complex.     

Phylogenetic Relationships among Holarctic Populations of Chironomus entis and Chironomus plumosus in View of Possible Horisontal Transfer of Mitochondrial Genes

In eight Holarctic populations of two typical chironomid sibling species of the plumosus group, Chironomus entisandChironomus plumosus, nucleotide sequences of mitochondrial (cytb) and nuclear (gb2b) gene regions were examined. The phylogenetic trees reflecting the evolutionary histories of the nuclear and mitochondrial markers exhibited significant differences. On the tree based on the nuclear gene sequences the populations clustered according to their species affiliation, whereas on the tree based on the mitochondrial gene sequences the populations were grouped according to their geographic position. This discrepancy is probably explained by mitochondrial gene flow between sympatric species with incomplete reproductive isolation (sibling species). Based on our results together with the earlier data on nuclear and mitochondrial gene sequences of some other species from the phylogenetic group plumosus, a scheme of phylogenetic relationships within this group is proposed. This scheme is in many ways different from the traditional view on the evolutionary relationships among species of the plumosus group.     

Identifying the rooted species tree from the distribution of unrooted gene trees under the coalescent

Gene trees are evolutionary trees representing the ancestry of genes sampled from multiple populations. Species trees represent populations of individuals—each with many genes—splitting into new populations or species. The coalescent process, which models ancestry of gene copies within populations, is often used to model the probability distribution of gene trees given a fixed species tree. This multispecies coalescent model provides a framework for phylogeneticists to infer species trees from gene trees using maximum likelihood or Bayesian approaches. Because the coalescent models a branching process over time, all trees are typically assumed to be rooted in this setting. Often, however, gene trees inferred by traditional phylogenetic methods are unrooted. We investigate probabilities of unrooted gene trees under the multispecies coalescent model. We show that when there are four species with one gene sampled per species, the distribution of unrooted gene tree topologies identifies the unrooted species tree topology and some, but not all, information in the species tree edges (branch lengths). The location of the root on the species tree is not identifiable in this situation. However, for 5 or more species with one gene sampled per species, we show that the distribution of unrooted gene tree topologies identifies the rooted species tree topology and all its internal branch lengths. The length of any pendant branch leading to a leaf of the species tree is also identifiable for any species from which more than one gene is sampled.     

GENE GENEALOGIES WITHIN THE ORGANISMAL PEDIGREES OF RANDOM-MATING POPULATIONS

Using computer simulations, we generated and analyzed genetic distances among selectively neutral haplotypes transmitted through gene genealogies with random-mating organismal pedigrees. Constraints and possible biases on haplotype distances due to correlated ancestry were evaluated by comparing observed distributions of distances to those predicted from an inbreeding theory that assumes independence among haplotype pairs. Results suggest that: 1) mean time to common ancestry of neutral haplotypes can be a reasonably good predictor of evolutionary effective population size; 2) the nonindependence of haplotype paths of descent within a given gene genealogy typically produces significant departures from the theoretical probability distributions of haplotype distances; 3) frequency distributions of distances between haplotypes drawn from “replicate” organismal pedigrees or from multiple unlinked loci within an organismal pedigree exhibit very close agreement with the theory for independent haplotypes. These results are relevant to interpretations of current molecular data on genetic distances among nonrecombining haplotypes at either nuclear or cytoplasmic loci.     

Why the Phylogenetic Species Concept?—Elementary

Although species play a number of unique and necessary roles in biology, none are more important than as the elements of phylogeny, nomenclature, and biodiversity study. Species are not divisible into any smaller units among which shared derived characters can be recognized with fidelity. Biodiversity inventory, assessment, and conservation are dependent upon a uniformly applicable species concept. Species are the fundamental units in formal Linnaean classification and zoological nomenclature. The Biological Species Concept, long given nominal support by most zoologists, forced an essentialy taxonomic problem (what are species?) into a population genetics framework (why are there species?). Early efforts at a phylogenetic species concept focused on correcting problems in the Biological Species Concept associated with ancestral populations, then applying phylogenetic logic to species themselves. Subsequently, Eldredge and Cracraft, and Nelson and Platnick, each proposed essentially identical and truly phylogenetic species concepts that permitted the rigorous recognition of species prior to and for the purposes of phylogenetic analysis, yet maintained the integrity of the Phylogenetic Species Concept outside of cladistic analysis. Such phylogenetic elements have many benefits, including giving to biology a unit species concept applicable across all kinds of living things including sexual and asexual forms. This is possible because the Phylogenetic Species Concept is based on patterns of character distributions and is therefore consistent with the full range of possible evolutionary processes that contribute to species formation, including both biotic and abiotic (even random) factors.     

Patterns of genetic diversity and biogeographical history of the tropical wetland tree, Pterocarpus officinalis (Jacq.), in the Caribbean basin

Studies examining intraspecific variation in plant species with widespread distributions and disjunct populations have mainly concentrated on temperate species. Here, we determined the genetic structure of a broadly distributed wetland tropical tree, Pterocarpus officinalis (Jacq.), from eight Neotropical populations using amplified length fragment polymorphisms (AFLP). AFLPs proved highly variable with almost half (48%) of the genetic variation at these loci occurring among individuals within populations. Nonetheless, there was a strong geographical pattern in the distribution of AFLP variation within P. officinalis. Caribbean and continental populations fell into two well-defined genetic clusters supported by the presence of a number of unique AFLP bands. Within these two regions, there were also strong genetic differences among populations, caused mainly by frequency differences in AFLP bands, making it difficult to determine the evolutionary relationships among populations. In addition, our analysis of P. officinalis revealed striking differences in the levels of AFLP variation among the eight populations sampled. In general, Caribbean populations had lower genetic diversity than continental populations. Moreover, there was a clear loss in AFLP diversity with distance from the continent among Caribbean populations. The overall genetic pattern within P. officinalis suggests that past colonization history, coupled with genetic drift within local populations, rather than contemporary gene flow are the major forces shaping variation within this species.