Differentiation between natural and cultivated populations of Medicago sativa (Leguminosae) from Spain: analysis with random amplified polymorphic DNA (RAPD) markers and comparison to allozymes

The conservation of a crop's wild relatives as genetic resources requires an understanding of the way genetic diversity is maintained in their populations, notably the effect of crop-to-wild gene flow. In this study, the amount of differentiation between natural and cultivated populations of Medicago sativa was analysed using random amplified polymorphic DNA (RAPD) markers and an extension of the AMOVA procedure adapted to autotetraploid organisms. Simulations of structured populations were performed to test whether AMOVA provides estimates of population structure in autotetraploids that can be directly compared to those obtained for allozyme data. Simulations showed that straight phi-statistics allow a good estimation of population differentiation when unbiased allelic frequencies are used to correct the conditional expectations of squared genetic distances. But such unbiased estimates can not be practically guaranteed, and population structure is notably overestimated when some populations are fixed for the presence of amplified fragments. However, removing fixed loci from the data set improves the statistical power of the test for population structure. The genetic variation of 15 natural and six cultivated populations of M. sativa was analysed at 25 RAPD loci and compared to estimates computed with allozymes on the same set of populations. Although RAPD markers revealed less within-population genetic diversity than allozymes, the quantitative and qualitative patterns of population structure were in full agreement with allozymes. This confirmed the conclusions drawn from the allozymic survey: crop-to-wild gene flow occurred in many locations, but some other mechanisms opposed cultivated traits to be maintained into natural populations.     

The effects of sample size on population genetic diversity estimates in song sparrows Melospiza melodia

To empirically determine the effects of sample size on commonly used measures of average genetic diversity, we genotyped 200 song sparrows Melospiza melodia from two populations, one genetically depauperate (n=100) and the other genetically diverse (n=100), using eight microsatellite loci. These genotypes were used to randomly create 10,000 datasets of differing sizes (5 to 50) for each population to determine what the effects of sample size might be on several estimates of genetic diversity (number of alleles per locus, average observed heterozygosity, and unbiased average expected heterozygosity) in natural populations of conservation concern. We found that at small sample sizes of 5 to 10 individuals, estimates of unbiased heterozygosity outperformed those based on observed heterozygosity or allelic diversity for both low- and high-diversity populations. We also found that when comparing across populations in which different numbers of individuals were sampled, rarefaction provided a useful way to compare estimates of allelic diversity. We recommend that standard errors should be reported for all diversity estimators, especially when sample sizes are small. We also recommend that at least 20 to 30 individuals be sampled in microsatellite studies that assess genetic diversity when working in a population that has an unknown level of diversity. However, research on critically endangered populations (where large sample sizes are impossible or extremely difficult to obtain) should include measures of genetic diversity even if sample sizes are less than ideal. These estimates can be useful in assessing the genetic diversity of the population.     

Assessing population genetic structure and variability with RAPD data: Application to Vaccinium macrocarpon (American Cranberry)

A method for estimating and comparing population genetic variation using random amplified polymorphic DNA (RAPD) profiling is presented. An analysis of molecular variance (AMOVA) is extended to accomodate phenotypic molecular data in diploid populations in Hardy-Weinberg equilibrium or with an assumed degree of selfing. We present a two step strategy: 1) Estimate RAPD site frequencies without preliminary assumptions on the unknown population structure, then perform significance testing for population substructuring. 2) If population structure is evident from the first step, use this data to calculate better estimates for RAPD site frequencies and sub-population variance components. A nonparametric test for the homogeneity of molecular variance (HOMOVA) is also presented. This test was designed to statistically test for differences in intrapopulational molecular variances (heteroscedasticity among populations). These theoretical developments are applied to a RAPD data set in Vaccinium macrocarpon (American cranberry) using small sample sizes, where a gradient of molecular diversity is found between central and marginal populations. The AMOVA and HOMOVA methods provide flexible population analysis tools when using data from RAPD or other DNA methods that provide many polymorphic markers with or without direct allelic data.     

Evaluation of factors affecting individual assignment precision using microsatellite data from horse breeds and simulated breed crosses

Assignment tests have been utilized to investigate population classification, measure genetic diversity and to solve forensic questions. Using microsatellite data from 26 loci genotyped in eight horse breeds we examined how population differentiation, number of scored loci, number of scored animals per breed and loci variability affected individual assignment precision applying log likelihood methods. We found that both genetic differentiation and number of scored loci were highly important for recognizing the breed of origin. When comparing two and two breeds, a proportion of 95% of the most differentiated breeds (0.200 < or = FST < or = 0.259) could be identified scoring only three loci, while the corresponding number was six for the least differentiated breeds (0.080 < or = FST < or = 0.139). An identical proportion of simulated breed crosses, differentiated from their parental breeds by FST estimates in the range 0.050-0.069, was identified when scoring 12 loci. This level of source identification was not obtained for the less differentiated breed crosses. The current data further suggested that population sample size and locus variability were not critical for the assignment precision as long as moderately large sample sizes (> or = 20 animals per population) and fairly variable loci were used.     

内蒙古典型草原羊草种群遗传分化的RAPD分析

运用 RAPD技术对内蒙古典型草原不同生境 8个羊草种群进行分析。采用 2 4个随机引物 (10 nt)在 8个种群中共检测到2 2 4个扩增片断 ,其中多态性片断 173个 ,总的多态位点百分率达 77.2 % ,特异性片断 2 2个 ,占 9.82 % ,平均每个引物扩增的DNA带数为 9.3 3条。利用 Nei指数和 Shannon指数估算了 8个种群的遗传多样性 ,并计算种群相似系数和遗传距离 ,运用UPGMA法进行聚类分析。结果表明 :羊草大部分的遗传变异存在于种群内 ,只有少部分的遗传变异存在于种群间 ,Nei指数和Shannon指数计算结果分别为 85.4%和 66.8% ;羊草不同种群的遗传多样性存在差异 ;8个羊草种群平均遗传距离为 0 .2 3 16,变异范围为 0 .1587～ 0 .2 70 0 ,说明 8个羊草种群间的遗传变异不大 ,即 :在较小地理范围内羊草的遗传分化程度较小 ;8个种群可聚为 3个类群 ,聚类结果显示生境相似的种群能够聚在一起 ,而地理距离最近的种群不一定归为一类 ,说明小范围内羊草种群间的遗传分化与地理距离不存在相关性 ,而与其生境间的相似度相关。影响遗传相似性的不是单一因子而是各种因子的综合作用 ,较小地理范围内羊草种群间的遗传分化主要是由环境的异质性所引起的     

History vs. current demography: explaining the genetic population structure of the common frog (Rana temporaria)

The amount of genetic variability at neutral marker loci is expected to decrease, and the degree of genetic differentiation among populations to increase, as a negative function of effective population size. We assessed the patterns of genetic variability and differentiation at seven microsatellite loci in the common frog (Rana temporaria) in a hierarchical sampling scheme involving three regions (208-885 km apart), three subregions within regions and nine populations (5-20 km apart) within subregions, and related the variability and differentiation estimates to variation in local population size estimates. Genetic variability within local populations decreased significantly with increasing latitude, as well as with decreasing population size and regional site occupancy (proportion of censured localities occupied). The positive relationship between population size and genetic variability estimates was evident also when the effect of latitude (cf. colonization history) was accounted for. Significant genetic differentiation was found at all hierarchical levels, and the degree of population differentiation tended to increase with increasing latitude. Isolation by distance was evident especially at the regional sampling level, and its strength increased significantly towards the north in concordance with decreasing census and marker-based neighbourhood size estimates. These results are in line with the conjecture that the influence of current demographic factors can override the influence of historical factors on species population genetic structure. Further, the observed reductions in genetic variability and increased degree of population differentiation towards the north are in line with theoretical and empirical treatments suggesting that effective population sizes decline towards the periphery of a species' range.     

Genetic variation at allozyme and RAPD loci in sessile oak Quercus petraea (Matt.) Liebl.: the role of history and geography

The nuclear genetic variation within and among 21 populations of sessile oak was estimated at 31 RAPD loci in conjunction with previous estimates of variation at eight allozyme loci. The aim of the study was to assess the relative role of isolation-by-distance and postglacial history on patterns of nuclear variation. Because of its small effective population size and maternal transmission, the chloroplast genome is a good marker of population history. Both kinds of nuclear variation (RAPD and allozyme) were therefore compared, first, to the geographical distances among populations and, secondly, to chloroplast DNA restriction polymorphism in the same populations. Multiple Mantel tests were used for this purpose. Although RAPDs revealed less genetic diversity than allozymes, levels of genetic differentiation ( G _ST) were identical. The standard genetic distance calculated at all RAPD loci was correlated with geographical distances but not with the genetic distance calculated from chloroplast DNA data. Conversely, allozyme variation was correlated with chloroplast DNA variation, but not with geography. Possibly, divergent selection at two allozyme loci during the glacial period could explain this pattern. Because of its greater number of loci assayed, RAPDs probably provided a less biased picture of the relative role of geography and history.     

Population genetics of the yellow fever mosquito in Trinidad: comparisons of amplified fragment length polymorphism (AFLP) and restriction fragment length polymorphism (RFLP) markers

Recent development of DNA markers provides powerful tools for population genetic analyses. Amplified fragment length polymorphism (AFLP) markers result from a polymerase chain reaction (PCR)-based DNA fingerprinting technique that can detect multiple restriction fragments in a single polyacrylamide gel, and thus are potentially useful for population genetic studies. Because AFLP markers have to be analysed as dominant loci in order to estimate population genetic diversity and genetic structure parameters, one must assume that dominant (amplified) alleles are identical in state, recessive (unamplified) alleles are identical in state, AFLP fragments segregate according to Mendelian expectations and that the genotypes of an AFLP locus are in Hardy-Weinberg equilibrium (HWE). The HWE assumption is untestable for natural populations using dominant markers. Restriction fragment length polymorphism (RFLP) markers segregate as codominant alleles, and can therefore be used to test the HWE assumption that is critical for analysing AFLP data. This study examined whether the dominant AFLP markers could provide accurate estimates of genetic variability for the Aedes aegypti mosquito populations of Trinidad, West Indies, by comparing genetic structure parameters using AFLP and RFLP markers. For AFLP markers, we tested a total of five primer combinations and scored 137 putative loci. For RFLP, we examined a total of eight mapped markers that provide a broad coverage of mosquito genome. The estimated average heterozygosity with AFLP markers was similar among the populations (0.39), and the observed average heterozygosity with RFLP markers varied from 0.44 to 0.58. The average FST (standardized among-population genetic variance) estimates were 0.033 for AFLP and 0.063 for RFLP markers. The genotypes at several RFLP loci were not in HWE, suggesting that the assumption critical for analysing AFLP data was invalid for some loci of the mosquito populations in Trinidad. Therefore, the results suggest that, compared with dominant molecular markers, codominant DNA markers provide better estimates of population genetic variability, and offer more statistical power for detecting population genetic structure.     

Using molecular markers with high mutation rates to obtain estimates of relative population size and to distinguish the effects of gene flow and mutation: a demonstration using data from endemic Mauritian skinks

We propose a method of analysing genetic data to obtain separate estimates of the size (N(p)) and migration rate (m(p)) for the sampled populations, without precise prior knowledge of mutation rates at each locus ( micro(L)). The effects of migration and mutation can be distinguished because high migration has the effect of reducing genetic differentiation across all loci, whereas a high mutation rate will only affect the locus in question. The method also takes account of any differences between the spectra of immigrant alleles and of new mutant alleles. If the genetic data come from a range of population sizes, and the loci have a range of mutation rates, it is possible to estimate the relative sizes of the different N(p) values, and likewise the m(p) and the micro(L). Microsatellite loci may also be particularly appropriate because loci with a high mutation rate can reach mutation-drift-migration equilibrium more quickly, and because the spectra of mutants arriving in a population can be particularly distinct from the immigrants. We demonstrate this principle using a microsatellite data set from Mauritian skinks. The method identifies low gene flow between a putative new species and populations of its sister species, whereas the differentiation of two other populations is attributed to small population size. These distinct interpretations were not readily apparent from conventional measures of genetic differentiation and gene diversity. When the method is evaluated using simulated data sets, it correctly distinguishes low gene flow from small population size. Loci that are not at mutation-migration-drift equilibrium can distort the parameter estimates slightly. We discuss strategies for detecting and overcoming this effect.     

Detection and analysis of genetic variation in Hordeum spontaneum populations from Israel using RAPD markers

Randomly amplified polymorphic DNA (RAPD) markers were used to analyse genetic diversity within and between Hordeum spontaneum populations sampled from Israel. Nei's index of genetic differentiation was used to partition diversity into within and between population components. Fifty-seven per cent of the variation detected was partitioned within 10 H. spontaneum populations. Using principal component and multiple regression analysis, part of the variation detected between populations was seen to be associated with certain ecogeographical factors. Fifty-eight per cent of the distribution of the phenotypic frequencies of three RAPD phenotypes detected using a single primer in 20 H. spontaneum populations could be accounted for by four ecogeographical variables, suggesting adaptive variation at certain RAPD loci.     

羊草种群遗传分化的RAPD分析Ⅱ.RAPD数据的统计分析

对松嫩草原上分布的灰绿型和黄绿型羊草9个种群进行了15个引物的RAPD分析,统计结果表明,两类种群的扩增片段数和多态位点比率明显不同,黄绿型种群低于灰绿型,其值分别＜90与＞100,＜50％与＞70％,比较了7种不同统计方法据RAPD表型或基因型频率估算的种群遗传多样性,几种统计结果都揭示,黄绿型种群低于灰绿型种群,用F1s值矫正种群对Hardy-Weinberg平衡的偏离后,估算等位基因频率,通过Shannon指数和Nei指数估计羊草种群间分化分别为37．6％和35．7％,高于等位酶的分析,讨论比较了等位酶和RAPD分析结果的异同。     

Sampling Issues Affecting Accuracy of Likelihood-based Classification Using Genetical Data

We demonstrate the effectiveness of a genetic algorithm for discovering multi-locus combinations that provide accurate individual assignment decisions and estimates of mixture composition based on likelihood classification. Using simulated data representing different levels of inter-population differentiation (F_st～ 0.01 and 0.10), genetic diversities (four or eight alleles per locus), and population sizes (20, 40, 100 individuals in baseline populations), we show that subsets of loci can be identified that provide comparable levels of accuracy in classification decisions relative to entire multi-locus data sets, where 5, 10, or 20 loci were considered. Microsatellite data sets from hatchery strains of lake trout, Salvelinus namaycush, representing a comparable range of inter-population levels of differentiation in allele frequencies confirmed simulation results. For both simulated and empirical data sets, assignment accuracy was achieved using fewer loci (e.g., three or four loci out of eight for empirical lake trout studies). Simulation results were used to investigate properties of the ‘leave-one-out’ (L1O) method for estimating assignment error rates. Accuracy of population assignments based on L1O methods should be viewed with caution under certain conditions, particularly when baseline population sample sizes are low (<50).     

Power for detecting genetic divergence: differences between statistical methods and marker loci

Information on statistical power is critical when planning investigations and evaluating empirical data, but actual power estimates are rarely presented in population genetic studies. We used computer simulations to assess and evaluate power when testing for genetic differentiation at multiple loci through combining test statistics or P values obtained by four different statistical approaches, viz. Pearson's chi-square, the log-likelihood ratio G-test, Fisher's exact test, and an F(ST)-based permutation test. Factors considered in the comparisons include the number of samples, their size, and the number and type of genetic marker loci. It is shown that power for detecting divergence may be substantial for frequently used sample sizes and sets of markers, also at quite low levels of differentiation. The choice of statistical method may be critical, though. For multi-allelic loci such as microsatellites, combining exact P values using Fisher's method is robust and generally provides a high resolving power. In contrast, for few-allele loci (e.g. allozymes and single nucleotide polymorphisms) and when making pairwise sample comparisons, this approach may yield a remarkably low power. In such situations chi-square typically represents a better alternative. The G-test without Williams's correction frequently tends to provide an unduly high proportion of false significances, and results from this test should be interpreted with great care. Our results are not confined to population genetic analyses but applicable to contingency testing in general.     

Genetic variation of wild and hatchery populations of the catla Indian major carp (Catla catla Hamilton 1822: Cypriniformes, Cyprinidae) revealed by RAPD markers

Genetic variation is a key component for improving a stock through selective breeding programs. Randomly amplified polymorphic DNA (RAPD) markers were used to assess genetic variation in three wild population of the catla carp (Catla catla Hamilton 1822) in the Halda, Jamuna and Padma rivers and one hatchery population in Bangladesh. Five decamer random primers were used to amplify RAPD markers from 30 fish from each population. Thirty of the 55 scorable bands were polymorphic, indicating some degree of genetic variation in all the populations. The proportion of polymorphic loci and gene diversity values reflected a relatively higher level of genetic variation in the Halda population. Sixteen of the 30 polymorphic loci showed a significant (p < 0.05, p < 0.01, p < 0.001) departure from homogeneity and the F(ST) values in the different populations indicated some degree of genetic differentiation in the population pairs. Estimated genetic distances between populations were directly correlated with geographical distances. The unweighted pair group method with averages (UPGMA) dendrogram showed two clusters, the Halda population forming one cluster and the other populations the second cluster. Genetic variation of C. catla is a useful trait for developing a good management strategy for maintaining genetic quality of the species.     

大别山山核桃天然群体遗传结构的初步分析

利用RAPD分子标记技术检测了大别山山核桃3个天然居群的遗传多样性和遗传结构。20条10 bp随机引物共检测到238个扩增位点,其中多态性位点162个,多态位点百分率(PPB)为68.1%。居群水平Shannon’s多态性信息指数(I)介于0.2651~0.2801之间;居群水平Ne i’s基因多样性指数(H)介于0.1789~0.1890之间。遗传变异计算显示大别山山核桃居群间基因分化系数(Gst)为0.4063,分子方差分析(AMOVA)表明居群间基因分化水平为0.4177,居群间基因流(Nm)为0.7306,说明大别山山核桃大部分变异存在于居群内,居群间基因交流相对较少。这一结果符合大别山山核桃风媒、异交的繁育系统特点,但其居群间基因分化程度明显高于异交植物的平均水平(Gst=0.1930)。地理隔离、居群内近交及居群间基因流受阻可能是形成目前大别山山核桃天然群体遗传结构的主要因素。     

Evolutionary implications of allozyme and RAPD variation in diploid populations of dioecious buffalograss Buchloë dactyloides

Buffalograss, Buchloë dactyloides, is widely distributed throughout the Great Plains of North America, where it is an important species for rangeland forage and soil conservation. The species consists of two widespread polyploid races, with narrowly endemic diploid populations known from two regions: central Mexico and Gulf Coast Texas. We describe and compare the patterns of allozyme and RAPD variation in the two diploid races, using a set of 48 individuals from Texas and Mexico (four population samples of 12 individuals each). Twelve of 22 allozyme loci were polymorphic, exhibiting 35 alleles, while seven 10-mer RAPD primers revealed 98 polymorphic bands. Strong regional differences were detected in the extent of allozyme polymorphism: Mexican populations exhibited more internal gene diversity (H_e= 0.20, 0.19) than did the Texan populations (H_e= 0.08, 0.06), although the number of RAPD bands in Texas (n= 62) was only marginally smaller than in Mexico (n= 68). F-statistics for the allozyme data, averaged over loci, revealed strong regional differentiation (mean F_RT=+ 0.30), as well as some differentiation among populations within regions (mean F_PR=+ 0.09). In order to describe and compare the partitioning of genetic variation for multiple allozyme and RAPD loci, we performed an Analysis of Molecular Variance (AMOVA). AMOVA for both allozyme and RAPD data revealed similar qualitative patterns: large regional differences and smaller (but significant) population differences within regions. RAPDs revealed greater variation among regions (58.4% of total variance) than allozymes (45.2%), but less variation among individuals within populations (31.9% for RAPDs vs. 45.2% for allozymes); the proportion of genetic variance among populations within regions was similar (9.7% for RAPDs vs. 9.6% for allozymes). Despite this large-scale concordance of allozyme and RAPD variation patterns, multiple correlation Mantel techniques revealed that the correlations were low on an individual by individual basis. Our findings of strong regional differences among the diploid races will facilitate further study of polyploid evolution in buffalograss.     

Estimating genetic drift and effective population size from temporal shifts in dominant gene marker frequencies

Measurement of temporal change in allele frequencies represents an indirect method for estimating the genetically effective size of populations. When allele frequencies are estimated for gene markers that display dominant gene expression, such as, e.g. random amplified polymorphic DNA (RAPD) and amplified fragment length polymorphism (AFLP) markers, the estimates can be seriously biased. We quantify bias for previous allele frequency estimators and present a new expression that is generally less biased and provides a more precise assessment of temporal allele frequency change. We further develop an estimator for effective population size that is appropriate when dealing with dominant gene markers. Comparison with estimates based on codominantly expressed genes, such as allozymes or microsatellites, indicates that about twice as many loci or sampled individuals are required when using dominant markers to achieve the same precision.     

Neighbour joining trees, dominant markers and population genetic structure

Population genetic theory for 'traditional' codominant loci showing low levels of allelic diversity (eg allozymes) has been well characterised and evaluated. In contrast, appropriate methods for the analysis of data from more recently developed marker systems are still being refined. For multilocus dominant markers such as amplified fragment length polymorphisms (AFLPs) and randomly amplified polymorphic DNA (RAPDs), the methods of data analysis can be split into two main categories. In population-based approaches, population allele frequencies are compared to obtain some measure of the partitioning of genetic diversity into within- and between-population components. In contrast, individual-based approaches use individual multilocus genotypes as the unit of analysis. Inferences on population processes such as gene flow are based on inter-relationships among individual samples as visualised on phenetic diagrams such as neighbour joining trees. Using a simulation approach coupled with neighbour joining analyses, we show that while the underlying population genetic structure is an important determinant of tree shape in the analysis of dominant data, the number of loci examined also affects the topology. At low levels of population differentiation (eg FST=0.07), mutually exclusive clustering of individuals into their respective populations can occur when sufficiently large numbers of loci are scored (eg 250 loci, typical of many AFLP studies). In contrast, unresolved star-shaped topologies can be recovered at higher levels of population differentiation (FST= >0.15) when lower numbers of loci are employed (eg 50 loci, typical of many RAPD studies). Thus, the relationship between tree topology and the extent of genetic structuring of populations is contingent upon the number of dominant loci scored. The consequences of these findings for the biological interpretation of individual-based analysis of dominant data sets are discussed.     

Assessment of the population structure of five Finnish dog breeds with microsatellites

Genetic variabilities within and between Finnish populations of Golden Retrievers, German Shepherds, Wirehaired Dachshunds, Pembroke Welsh Corgis and Bedlington Terriers were quantified with microsatellite allele numbers, observed heterozygosities, expected heterozygosities, FIS estimates, FST estimates and DS distances. In a sample of 50 individuals from each, breed and ten polymorphic loci, the highest genetic diversity was exhibited in the Wirehaired Dachshunds and lowest in the Bedlington Terriers. Although statistically significant deviations from the Hardy-Weinberg (H-W) equilibrium were observed, they occurred at an unexpectedly low frequency. Most strikingly, the extremely small Bedlington Terrier population displayed genotypes in H-W proportions in all investigated loci. The H-W deviations always occurred with positive FIS estimates, which, on average, were not larger than values reported for free-living canids. Genetic differentiation between the breeds was very large. As a comparison, present estimates were, on average, over two times higher than previously observed between breeds of sheep, and over two times higher than the highest estimates reported between human populations. Moreover, the highest DS distances were only slightly lower than the lowest values inferred between humans and chimpanzees. Severe bottlenecks in the recent past of the examined breeds were not statistically supported. The presented data imply genetic isolation and intense artificial selection in the history of these breeds of dogs.     

DNA EXTRACTION AND RAPD MARKERS TO ASSESS THE GENETIC SIMILARITY AMONG GELIDIUM SESQUIPEDALE (RHODOPHYTA) POPULATIONS1

A simple method developed for genomic DNA isolation from fungus was tested on the red alga, Gelidium sesquipedale (Clem.) Born et Thur., which is commercially exploited for its high sulfated polysaccharide (agar) content. This method is faster, cheaper, and less toxic than conventional phenol/chloroform methods. Random amplified polymorphic DNA (RAPD) amplifications were performed successfully without the necessity of purifying the DNA. RAPD markers were used to investigate the genetic similarity among three natural populations of G. sesquipedale from southern Portugal. Bulked-genomic DNA samples of 15 different individuals were made in each population. These can be conceived of as a sample of the population DNA. Of the 62 primers screened, 41 produced bands and 22 revealed polymorphisms. Genetic similarities among populations were high. Populations that are further away from each other have the lowest similarity coefficients, whereas the intermediate Ingrina population, located on the south coast, showed higher genetic similarity with the Odeceixe population located on the southwest coast, than with the São Rafael southern population. This suggests a higher genetic flow between Odeceixe and Ingrina or the result may be a founder effect in the sense that the species has propagated from the east coast to the south coast of Portugal. We conclude that the use of this isolation method with RAPD analysis is appropriate to characterize the genetic variability of this commercial species along its geographical distribution. Large sample sizes can be screened at a relatively low cost. Finding genetic markers for commercial populations of G. sesquipedale may be of industrial interest.