Journal of Heredity Advance Access originally published online on June 26, 2008
Journal of Heredity 2008 99(5):476-482; doi:10.1093/jhered/esn055
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Genetic Structure and Mating System of Euterpe edulis Mart. Populations: A Comparative Analysis Using Microsatellite and Allozyme Markers
Monsanto do Brasil, Uberlandia Research Station, Rodovia Uberlandia/Araxa, BR 452, Km 149, 38405-232, Uberlandia, MG—Brazil, (Conte); Universidade Federal de Santa Catarina, Departamento de Fitotecnia, Rodovia Admar Gonzaga, 1346, CP 476, 88040-900, Florianópolis, SC, Brazil (Sedrez dos Reis); Universidade do Estado de Santa Catarina, Departamento de Engenharia Florestal, Av. Luiz de Camões, 2090, 88520-000, Lages, SC, Brazil (Mantovani); and Universidade de São Paulo, Escola Superior de Agricultura "Luiz de Queiroz", Departamento de Genética, Av. Padua Dias, 11, CP 83, 13400-970, Piracicaba, SP, Brazil (Vencovsky)
Address correspondence to Rudimar Conte at the address above, or e-mail: rudimar.conte{at}monsanto.com.
A comparative study between microsatellite and allozyme markers was conducted on the genetic structure and mating system in natural populations of Euterpe edulis Mart. Three cohorts, including seedlings, saplings, and adults, were examined in 4 populations using 10 allozyme loci and 10 microsatellite loci. As expected, microsatellite markers had a much higher degree of polymorphism than allozymes, but estimates of multilocus outcrossing rate (
= 1.00), as well as estimates of genetic structure (FIS, GST), were similar for the 2 sets of markers. Estimates of RST, for microsatellites, were higher than those of GST, but results of both statistics revealed a close agreement for the genetic structure of the species. This study provides support for the important conclusion that allozymes are still useful and reliable markers to estimate population genetic parameters. Effects of sample size on estimates from hypervariable loci are also discussed in this paper.
In the last 4 decades, the main genetic markers used in population studies in most species have been allozymes (Hedrick 1999). Although the number of molecular markers (DNA markers) has increased significantly in recent years, allozymes still continue to be widely used in population genetic studies (Sebbenn et al. 2000; Mariot et al. 2002; Lumaret et al. 2004; Cuevas et al. 2006; Haag et al. 2006; Mantovani et al. 2006). An important aspect of these markers is related to their neutrality. Although most studies using allozyme markers assume selective neutrality, some studies have pointed out the existence of adaptive value associated with protein polymorphism (Eanes 1999; Conte et al. 2003; Borrell et al. 2004; Araki et al. 2005; Volis et al. 2005). This fact has generated some controversy about the usefulness of these markers in population genetics because natural selection may influence the maintenance and population organization of genetic polymorphism. However, there is evidence that only a few enzymes are actually under selection and that the rest must be neutral or nearly neutral (Eanes 1999; Thelen and Allendorf 2001; Volis et al. 2005).
The recent development of techniques for the examination of hypervariable regions of the genome has significantly raised the number of applications of genetic markers (Haig 1998; Zane et al. 2002). Currently, microsatellites or simple sequence repeats (SSRs) are a popular tool in population and conservation genetics (Collevatti et al. 2001; Gao et al. 2002; Gaiotto et al. 2003; Dutech et al. 2004). Because most microsatellites are located in noncoding regions and no direct evidence exists for a general role of SSRs in eukaryotic genomes, they are more likely to be neutral than functional markers such as allozymes (Queller et al. 1993; Oliveira et al. 2006). However, evidence for nonneutrality at particular microsatellite loci is accumulating (Larsson et al. 2007). The use of microsatellites has made a great contribution toward the improvement of genetic parameter estimates, especially in populations possessing insufficient diversity to be assessed with other techniques (Oliveira et al. 2006).
The most important difference between the 2 categories of markers is their level of polymorphism. Microsatellites are much more polymorphic than allozymes in natural populations, with more alleles, and heterozygosities are often greater than 0.5 (Estoup et al. 1998). Mutation rates at microsatellite loci have been estimated to be in the order of 5 x 10–5 to 5 x 10–3 (Amos et al. 1996), that is 2–4 orders of magnitude higher than at allozyme loci (Estoup et al. 2002). Comparative analyses of microsatellite and allozyme markers in natural populations have revealed some divergences between the 2 markers, and some advantages have been pointed out for microsatellites (Estoup et al. 1998; Streiff et al. 1998; Gao et al. 2002). In particular, the higher level of variation at microsatellite loci suggests that they should be more powerful than allozymes to detect changes in effective population size, genetic structure, and outcrossing rates. However, according to Hedrick (1999), data from highly variable loci must be carefully evaluated because the information given may be quite different than that obtained from less variable markers.
Heart of palm (Euterpe edulis Mart.; Arecaceae) is a native shade-tolerant palm of the Brazilian Atlantic forest which is economically valuable because it is the main species used for palm heart extraction. Genetic studies using microsatellite and allozyme markers in natural populations have demonstrated high values of outcrossing rates and gene flow for the species together with high levels of genetic variability within populations and low interpopulation genetic divergence (Reis et al. 2000c; Conte et al. 2003; Gaiotto et al. 2003). In natural populations, the species presents a high density and a J-shaped age structure (Reis et al. 1996, 2000a). It is a monoecious species with strongly protandrous flowers pollinated by insects (Mantovani and Morellato 2000). Currently, natural populations are intensely fragmented and reduced in its area of occurrence (Reis et al. 2000b). Although a sustainable management system has been proposed for the species (Reis et al. 2000a), the intensive exploitation occurring since the 1960s has resulted in the elimination of several populations and serious alterations to remaining populations (Galetti and Fernandez 1998).
In this study, we combined analyses of microsatellite and allozyme loci in order to compare the genetic structure and mating system of E. edulis Mart. populations. Our working hypothesis was based on the assumption that, as the SSR loci are highly multiallelic, they improve resolution over allozymes for the estimation of population genetic parameters.
| Material and Methods |
|---|
|
|
|---|
Populations and Sampling
Four natural populations of E. edulis were surveyed in the districts of São Pedro de Alcântara and Ibirama, Santa Catarina, Brazil (Figure 1). The vegetation of the study sites is defined as "Evergreen Atlantic Tropical Forest." In spite of the current fragmentation and degradation of the populations of E. edulis, many remnant fragments exist that are well conserved, as is the case of the populations here investigated. Also, the fragments still retain connectivity due to their mutual proximity, such that disturbed populations will benefit from the influence of the undisturbed populations in the surrounding area.
|
To study variability and genetic structure, leaf material from 3 age cohorts was sampled randomly in the 4 populations as follows: 1) Seedlings—insertion height of the youngest leaf less than 10 cm, corresponding to the individuals of the last reproductive period; 2) Saplings—insertion height from 30 cm to 1 m, absence of exposed stem and the presence of 4–5 mature leaves; and 3) Adults—plants usually higher than 10 m of height, presence of open inflorescence and/or abscission scars on the stem from previous events of reproduction. A total of 50 individuals were sampled per cohort in each population. For the mating system study, open-pollinated seeds were collected from 15 randomly selected adult trees in São Pedro (Population 4). Seeds were germinated and grown in a greenhouse under controlled temperature (20–35 °C) and natural photoperiod. Progeny arrays of 13 individuals per open-pollinated family were used to determine outcrossing rates.
Genetic Analysis
In total, 10 enzymatic systems were used:
-esterase (
-EST, EC 3.1.1.1
[EC]
), shikimic dehydrogenase (SKDH, EC 1.1.1.25
[EC]
), peroxidase (PRX, EC 1.11.1.7
[EC]
), 6-phosphogluconate dehydrogenase (PGDH, EC 1.1.1.44
[EC]
), nicotinamide adenine dinucleotide dehydrogenase (NADHDH, EC 1.6.99.3
[EC]
), malate dehydrogenase (MDH, EC 1.1.1.37
[EC]
), phosphoglucoisomerase (PGI, EC 5.3.1.9
[EC]
), isocitrate dehydrogenase (IDH, EC 1.1.1.42
[EC]
), phosphoglucomutase (PGM, EC 5.4.2.2
[EC]
), and glucose-6-phosphate dehydrogenase (G6PDH, EC 1.1.1.49
[EC]
). Among the 10 systems employed, 16 zones with enzymatic activity were observed, revealing 10 presumed polymorphic loci. Each zone that showed a behavior apparently independent of the others and with Mendelian segregation was considered an allozyme locus, taking into account the possible quaternary structure of each enzyme in question (Kephart 1990). Only the 10 polymorphic loci were used in the genetic analyses.
A group of 10 SSR loci, previously developed and optimized for E. edulis (Gaiotto et al. 2001), were used to genotype all sampled individuals. The total reaction volume used in polymerase chain reaction was 13 µl, containing 7.5 ng of genomic DNA; 50 mM KCl, 20 mM Tris–HCl pH 8.4; 10% dimethyl sulfoxide; 1.5 mM MgCl2; 250 µM of each deoxynucleoside triphosphates; 0.3 µM of each primer (forward and reverse); 1.0 unit of Taq DNA polymerase (Ludwig Biotecnologia Porto Alegre, Brazil); and ultrapure water to complete the final volume. Amplifications were performed using a PTC-100 thermal controller (MJ Research, Waltham, MA) with the following conditions: 94 °C for 5 min, 30 cycles of 94 °C for 1 min, the primer-specific annealing Ta temperature (Gaiotto et al. 2001) for 1 min, 72 °C for 1 min, ending with 72 °C for 7 min. The amplified fragments were separated in 4% polyacrylamide gel in a run with 1x Tris/Borate/EDTA at 60 W for 1 h, using 10-bp ladder size standard, and stained with silver nitrate (Creste et al. 2001).
Data Analysis
Molecular genetic variability was described in terms of average number of alleles per locus (A), observed heterozygosity (Ho), Nei's (1978) gene diversity (He), and fixation index (f) using the GDA program (Lewis and Zaykin 2001). The fixation index (f or FIS) was estimated for each population by analysis of variance (ANOVA), and the significance was tested by bootstrapping over loci using 10 000 replicates in order to estimate interlocus variation. Genotypic frequencies of each locus were subjected to a goodness-of-fit test (Fisher's exact test) to expected Hardy–Weinberg proportions, as defined by Weir (1996), using the GDA program (Lewis and Zaykin 2001). Such tests were performed by the conventional Monte Carlo method using 10 batches with 1000 permutations per batch.
Molecular genetic structure among populations was investigated by the analysis of gene diversity in subdivided populations (Nei 1973). Parameters HS, HT, DST, and GST were estimated using the FSTAT program (Goudet 2002), and the null hypothesis was tested by bootstrapping over loci using 10 000 replicates. As most mutations within microsatellites involve the addition or subtraction of a small number of repeat units, the mutation process is not in line with expectations under an infinite allele model with low mutation rates. Therefore, an analogue of the FST statistics, namely the RST parameter (Slatkin 1995), developed for microsatellite data, was also used to quantify genetic structure. ANOVA of allele size followed Goodman (1997) using the FSTAT program (Goudet 2002).
The mating system of E. edulis was analyzed under the mixed mating model described by Ritland and Jain (1981) using the MLTR program (Ritland 1997). This procedure estimates the multilocus outcrossing rate (tm), the average single-locus outcrossing rate (ts), the outcrossing rate between related individuals [(tm – ts)
0], the average single locus inbreeding coefficient of maternal parents (f), the correlation of outcrossing rate within progeny arrays (rt), the correlation of outcrossed paternity within progeny arrays (or the probability that randomly chosen progeny pairs from the same array are full sibs; rp), and the pollen (p) and ovule (o) gene population frequencies. Standard errors (SEs) were estimated by bootstrapping over progenies using 10 000 replicates.
| Results |
|---|
|
|
|---|
Number of Alleles at Allozyme and Microsatellite Loci
From the allozyme analysis, a total of 41 alleles at 16 loci were identified in the 3 cohorts of the 4 populations. From these loci, 6 were monomorphic (
-Est-1,
-Est-2, Pgdh-1, Nadhdh-1, Mdh-2, and Idh-1) and 10 were polymorphic (Skdh-1, Prx-2, Prx-3, Prx-4, Prx-5, Pgdh-2, Mdh-1, Pgi-2, Pgm-1, and G6pdh-1). The locus Prx-5 showed a very low level of polymorphism, with the most frequent allele close to fixation in all cohorts studied. However, some loci showed high levels of polymorphism (Prx-3, Prx-4, Pgm-1, with 4 alleles each, and Pgdh-2, with 5 alleles). All 10 microsatellite loci displayed high levels of polymorphism and a total of 161 alleles were described. The most variable loci (EE48 and EE47) displayed 24 and 26 alleles, respectively, whereas the least variable (EE5) displayed 13 alleles.
Genetic Variability
The estimates of genetic variation within populations for allozymes and microsatellite markers are presented in Table 1. An average of 52 individuals was analyzed per cohort in each population. As expected, the values of A, Ho, and He were much higher for microsatellites than allozymes. For allozymes, the following diversity estimates were obtained among the cohorts:
: 3.05–3.15;
e: 0.416–0.431;
o: 0.378–0.403. For microsatellites, the estimates were as follows:
: 14.12–14.72;
e: 0.781–0.785;
o: 0.678–0.709.
|
Although considerable variation was observed among populations, both markers displayed similar trend within populations (Table 1). In general, both markers displayed lower
Departures from Hardy–Weinberg equilibrium were observed in a greater extent for microsatellites than allozymes, using Fisher's exact test. Whereas allozymes presented an average of 70% of the loci at equilibrium for most populations, the average for microsatellites was 50%; some populations, however, presented values which were even lower for seedlings and saplings (37.5%). On the other hand, as Fishers exact test was performed for several loci, the number of loci where the null hypothesis was rejected decreased when the Bonferroni correction was applied (Weir 1996).
Genetic Structure
The analysis of gene diversity in subdivided populations (Nei 1973) showed a high congruence between the 2 markers, both of which revealed low molecular genetic divergence among populations. On average, the estimated values of GST for allozyme loci were low (0.023, 0.017, and 0.011 for seedlings, saplings, and adults, respectively), but significantly different from 0. Similarly, for microsatellites estimates of GST were also low (0.024, 0.021, and 0.028, respectively) but significantly different than 0. In the pairwise comparisons, the estimated values of GST varied from 0.013 to 0.038 for microsatellites and were significantly different than 0 in all comparisons performed with this marker. For allozymes, more variation of GST was observed among the different comparisons (0.008–0.041), and the majority were not significantly different than 0, especially for populations of the same region.
The estimates of RST over all populations ranged from 0.033 to 0.048 among the 3 cohorts. In the pairwise comparisons, these estimates ranged from 0.004 to 0.091 and mostly were significantly different from 0. Despite the fact that the estimates of RST were higher than those of GST, results of both statistics revealed a close agreement, indicating that the divergence was mostly among regions.
Mating system
SEs of pollen and ovule gene frequencies were low for both markers (data not shown), indicating that the sample size used was adequate for the estimates of mating system parameters. Most differences in allele frequencies between pollen and ovules were not significant on the basis of SEs, indicating homogeneity of the pollen pool by each maternal tree. A great congruence between the 2 markers was observed in the estimates of mating system parameters (Table 2). The multilocus outcrossing rate (
) was equal to 0.996 (0.006) for microsatellites and 1.000 (0.000) for allozymes, revealing preferentially allogamic behavior of E. edulis. The single locus outcrossing rate (
) was high (0.987 for allozymes and 0.929 for microsatellites), but significantly different from the multilocus estimates, especially for microsatellites, suggesting that crosses among related individuals had occurred in the population. The fixation index for maternal genotypes (
) was null in both markers, indicating the absence of inbreeding. The probability that a randomly chosen pair of progenies from the same array were full sibs
was equal to 0.109 (0.000) for allozymes and 0.103 (0.012) for microsatellites, suggesting that the effective number of pollen parents (Nep = 1/rp) is somewhat less than 10 individuals. The normalized variation of outcrossing rate among progeny arrays
was low and nonsignificant for both markers, which is in agreement with the floral biology of this strongly protandrous species.
|
| Discussion |
|---|
|
|
|---|
The 10 allozyme and 10 SSR loci used for genetic analysis in this study detected high levels of molecular genetic variation, confirming the results obtained in previous studies with E. edulis populations with allozymes (Reis et al. 2000c; Conte et al. 2003) and with microsatellites (Gaiotto et al. 2003). As expected, microsatellite markers showed a higher degree of polymorphism than allozymes for all cohorts examined. Differences in within-population genetic diversity (He) were detected for microsatellites, as observed in some studies (Estoup et al. 1998; Streiff et al. 1998), but not detected in others (e.g., Freville et al. 2001).
Differences in estimates of He among cohorts were of the same magnitude for the 2 sets of markers. Therefore, both markers expressed similar inference on the genetic diversity of the species.
Our study also showed departures from Hardy–Weinberg equilibrium for both markers, mainly for microsatellite loci. However, observed departures from Hardy–Weinberg at some loci may derive from statistical sampling (Weir 1996; Collevatti et al. 2001) beacause the microsatellite loci used in this work are highly polymorphic and the number of individuals genotyped per cohort per population (
50) is relatively limited for the sampling of all possible genotypes at a locus. This aspect is particularly important in population genetic studies once the use of hypervariable loci might give raise to departures in both genotype and allele frequencies, which may substantially affect the accuracy of the estimates of genetic parameters. In spite of the well-known low contribution of rare alleles to estimated f values, the sample size used in our study might have contributed to the significant inbreeding values among seedlings and samplings detected with microsatellites. Besides, we cannot discard the possible occurrence of null alleles, which might overestimate the homozygosity values, affecting the inbreeding estimates, as discussed by Ally and Ritland (2007).
As discussed by Sun et al. (2001), GST is a reliable statistic for the comparison between allozymes and microsatellite once it is based on the relative distribution of variability within and among populations. In our study, on average, the estimated values of GST for allozymes and microsatellites between all analyzed cohorts were very similar. This result is congruent to those reported in other studies (Estoup et al. 1998; Sun et al. 2001). According to Estoup et al. (1998), the similarity on FST estimates suggests that selection is not operating strongly on the genetic differentiation at both loci.
The results from both markers revealed that most of the molecular genetic variability is distributed within populations, confirming the results obtained in other studies with this species (Reis et al. 2000c; Conte et al. 2003; Gaiotto et al. 2003). The small amount of population genetic differentiation (
ST,
ST) is also consistent with the high rate of outcrossing of this species and the localized distribution of the study populations. Moreover, our results are in agreement with other studies on tropical species, where species with large geographic ranges, outcrossing breeding systems, and high rates of gene flow, have more genetic variability within populations and consequently low divergence among populations (Hamrick and Godt 1990).
Hedrick (1999) mentions that because microsatellite loci often have very high heterozygosity within populations, the magnitude of differentiation may be quite small. However, similar magnitude of genetic divergence was also observed using allozyme markers, which demonstrates the consistency of the genetic information obtained with these 2 markers. Similarly, an examination of the estimates of molecular genetic variability, such as the number of alleles per locus and the gene diversity, revealed that all 4 populations exhibited similar levels of variability for the 3 categories studied.
The multilocus outcrossing rate (tm) in the present study (
= 1.000 for allozymes and
= 0.996 for microsatellites) confirmed the results obtained in previous studies of this species ranging from 0.940 to 1.040 for allozymes (Reis et al. 2000c) and from 0.910 to 0.980 for microsatellites (Gaiotto et al. 2003). These results are all consistent with the conclusion that E. edulis is an obligately outcrossing species. The estimates of tm are also congruent with the floral biology and pollination ecology of the species (Mantovani and Morellato 2000).
Despite the high rate of outcrossing found in this species, nonrandom matings have been observed by the differences between the multilocus (
) and single locus (
) estimates for both markers, indicating that crosses among related individuals had occurred in the population. The occurrence of nonrandom matings was also observed in the estimates of correlation of outcrossed paternity within progeny arrays (
), revealing a substantial number of full sibs originating from a limited pool of pollen, a possible result of limited pollen dispersal and/or nonsynchronous flowering of all individuals. The larger SEs observed for microsatellites is the result of less accurate estimates typical for highly multiallelic loci.
In summary, SSR loci showed a much higher degree of absolute polymorphism than allozymes, but relative measures of the mating system and genetic structure were quite similar for the 2 sets of markers. Both markers displayed few pronounced genetic structure of E. edulis, with the exception of a small increase in the inbreeding levels among seedlings and saplings of the populations 1 and 2 (private areas).
In spite of limited number of estimates provided in this study, the present results give support for the important prediction that allozymes are still useful and reliable markers to estimate population genetic parameters. Regardless of the limitation of sample size for the sampling of all possible genotypes, microsatellites also confirmed the high content of genetic information for studies of population genetic parameters.
| Funding |
|---|
|
|
|---|
Coordenação de Aperfeiçoamento de Pessoal de Ensino Superior and Conselho Nacional de Desenvolvimento Científico e Tecnológico.
| Acknowledgments |
|---|
We thank the researchers of the Núcleo de Pesquisas em Florestas Tropicais/Universidade Federal de Santa Catarina for their assistance in collecting and organizing the samples and also the Laboratório de Fisiologia do Desenvolvimento e Genética Vegetal/Universidade Federal de Santa Catarina for the infrastructure provided in support of this work.
| Footnotes |
|---|
Corresponding Editor: Jim Hamrick
Received August 7, 2007
Accepted May 22, 2008
| References |
|---|
|
|
|---|
-
Ally D, Ritland K. A case study: looking at the effects of fragmentation on genetic structure indifferent life history stages of Old-Growth Mountain Hemlock (Tsuga mertensiana). J Hered (2007) 98(1):73–78.
Amos W, Sawcer SJ, Feakes RW, Rubinzstein DC. Microsatellites show directional bias and heterozygote instability. Nat Genet (1996) 13:390–391.[CrossRef][Web of Science][Medline]
Araki H, Yoshizumi S, Inomata N, Yamazaki T. Genetic coadaptation of the Amylase gene system in Drosophila melanogaster: evidence for the selective advantage of the lowest AMY activity and of its epistatic genetic background. J Hered (2005) 96:388–395.
Borrell YJ, Pineda H, McCarthy I, Vázquez E, Sánchez JA, Lizana GB. Correlation between fitness and heterozygosity at allozyme and microsatellite loci in the Atlantic salmon, Salmo salar L. Heredity (2004) 92:585–593.[CrossRef][Web of Science][Medline]
Collevatti RG, Grattapaglia D, Hay JD. Population genetic structure of the endangered tropical tree species Caryocar brasiliense, based on variability at microsatellite loci. Mol Ecol (2001) 10:349–356.[CrossRef][Medline]
Conte R, Nodari RO, Vencovsky R, Reis MS. Genetic diversity and recruitment of the tropical palm, Euterpe edulis Mart. in a natural population from the Brazilian Atlantic Forest. Heredity (2003) 91:401–406.[CrossRef][Web of Science][Medline]
Creste C, Tulmann-Neto A, Figueira A. Detection of single sequence repeat polymorphisms in denaturing polyacrylamide sequencing gels by silver staining. Plant Mol Biol Rep (2001) 19:299–306.[CrossRef][Web of Science]
Cuevas E, Arias DM, Domínguez CA, Castillo RA, Molina-Freaner F. The genetic structure of the gynodioecious Kallstroemia grandiflora (Zygophyllaceae): the role of male sterility and colonization history. Heredity (2006) 97:269–274.[CrossRef][Web of Science][Medline]
Dutech C, Joly HI, Jarne P. Gene flow, historical population dynamics and genetic diversity within French Guianan populations of a rainforest tree species, Vouacapoua americana. Heredity (2004) 92:69–77.[CrossRef][Web of Science][Medline]
Eanes WF. Analysis of selection on enzyme polymorphisms. Annu Rev Ecol Syst (1999) 30:301–326.[CrossRef][Web of Science]
Estoup A, Jarne P, Cornuet JM. Homoplasy and mutation model at microsatellite loci and their consequences for population genetics analysis. Mol Ecol (2002) 11:1591–1604.[CrossRef][Medline]
Estoup A, Rousset F, Michalakis Y, Cornuet JM, Adriamanga M, Guyomard R. Comparative analysis of microsatellite and allozyme markers: a case study investigating microgeographic differentiation in brown trout (Salmo trutta). Mol Ecol (1998) 7:339–353.[CrossRef][Medline]
Freville H, Justy F, Olivieri I. Comparative allozyme and microsatellite population structure in a narrow endemic plant species, Centaurea corymbosa Pourret (Asteraceae). Mol Ecol (2001) 10:879–889.[CrossRef][Medline]
Gaiotto FA, Brondani RPV, Grattapaglia D. Microsatellite markers for heart of palm—Euterpe edulis and E. oleracea Mart. (Arecaceae). Mol Ecol (2001) 1:86–88.
Gaiotto FA, Grattapaglia D, Vencovsky R. Genetic structure, mating system, and long-distance gene flow in heart of palm (Euterpe edulis Mart). J Hered (2003) 94:399–406.
Galetti M, Fernandez JC. Palm heart harvesting in the Brazilian Atlantic forest: changes in industry structure and the illegal trade. J Appl Ecol (1998) 35:294–301.[CrossRef]
Gao LZ, Schaal BA, Zhang CH, Jia JZ, Dong YS. Assessment of population genetic structure in common wild rice Oryza rufipogon Griff. using microssatellite and allozyme markers. Theor Appl Genet (2002) 106:173–180.[Web of Science][Medline]
Goodman SJ. RST Calc. A collection of computer program for calculating estimates of genetic differentiation from microsatellite and determining their significance. Mol Ecol (1997) 6:881–885.[CrossRef]
Goudet J. FSTAT. A computer program to estimate and test gene diversities and fixation indices (Version 2.9.3.2). (2002) Available from: http://www2.unil.ch/popgen/softwares/fstat.htm.
Haag CR, Riek M, Hottinger JW, Pajunen VI, Ebert D. Founder events as determinants of within-island and among-island genetic structure of Daphnia metapopulations. Heredity (2006) 96:150–158.[CrossRef][Web of Science][Medline]
Haig SM. Molecular contributions to conservation. Ecology (1998) 79:413–425.[CrossRef][Web of Science]
Hamrick JL, Godt MJW. Allozyme diversity in plant species. In: Plant population genetics, breeding and genetic resources—Brown AHD, Clegg MT, Kahler AL, Weir BS, eds. (1990) Sunderland (MA): Sinauer Associates. 43–63.
Hedrick PW. Highly variable loci and their interpretation in evolution and conservation. Evolution (1999) 53:313–318.[CrossRef][Web of Science]
Kephart SR. Starch gel electrophoresis of plant isozymes: a comparative analysis of techniques. Am J Bot (1990) 77:693–712.[CrossRef][Web of Science]
Larsson LC, Laikre L, Palm S, André C, Carvalho GR, Ryman N. Concordance of allozyme and microsatellite differentiation in a marine fish, but evidence of selection at a microsatellite locus. Mol Ecol (2007) 16:1135–1147.[CrossRef][Medline]
Lewis P, Zaykin D. GDA. A computer program for the analysis of allelic data (Version 1.0d16c). (2001) Available from: [cited July 15, 2007] http://hydrodictyon.eeb.uconn.edu/people/plewis/software.php.
Lumaret R, Ouazzani N, Michaud H, Vivier G, Deguilloux M-F, Giusto FD. Allozyme variation of oleaster populations (wild olive tree) (Olea europaea L.) in the Mediterranean Basin. Heredity (2004) 92:343–351.[CrossRef][Web of Science][Medline]
Mantovani A, Morellato LPC. Fenologia da floração, frutificação, mudança foliar e aspectos da biologia floral. Sellowia (2000) 49–52:23–38.
Mantovani A, Morellato LPC, Reis MS. Internal genetic structure and outcrossing rate in natural population of Araucaria angustifolia (Bert.) O. Kuntze. J Hered (2006) 97:466–472.
Mariot A, Di Stasi LC, Reis MS. Genetic diversity in natural populations of Piper cernuum. J Hered (2002) 93:365–369.
Nei M. Analysis of gene diversity in subdivided populations. Proc Natl Acad Sci USA (1973) 70:3321–3323.
Nei M. Estimation of average heterozygosity and genetic distance from a small number of individuals. Genetics (1978) 89:583–590.
Oliveira EJ, Pádua JG, Zucchi MI, Vencovsky R, Vieira MLC. Origin, evolution and genome distribution of microsatellites. Genet Mol Biol (2006) 29:294–307.
Queller DC, Strassmann JE, Hughes CR. Microsatellites and kinship Tree. Trends in Ecol and Evol. (1993) 8:285–288.[CrossRef]
Reis MS, Fantini AC, Nodari RO, Reis A, Guerra MP, Mantovani A. Management and conservation of natural populations in Atlantic Rain forest: the case study of palm heart (Euterpe edulis Martius). Biotropica (2000a) 32:894–902.[Web of Science]
Reis MS, Guerra MP, Nodari RO, Reis A, Ribeiro RJ. Distribuição geográfica e situação atual das populações na área de ocorrência de Euterpe edulis Martius. Sellowia (2000b) 49–52:324–335.
Reis MS, Kageyama PY, Guimarães E, Nodari RO, Fantini AC, Mantovani A, Vencovsky R. Variação genética em populações naturais de Euterpe edulis Martius na Floresta Ombrófila Densa. Sellowia (2000c) 49–52:131–149.
Reis A, Kageyama PY, Reis MS, Fantini AC. Demografia de Euterpe edulis Martius (Arecaceae) em uma Floresta Ombrófila Densa Montana, em Blumenau (SC). Sellowia (1996) 45–48:13–45.
Ritland K. MLTR multilocus mating system program (version 1.1). (1997) Available from: http//genetics.forestry.ubc.ca/ritland/programs.html.
Ritland K, Jain SK. A model for the estimation of outcrossing rate and gene frequencies using independent loci. Heredity (1981) 47:35–52.[CrossRef][Web of Science]
Sebbenn AM, Seoane CES, Kageyama PY, Vencovsky R. Effects of the management on the genetic structure of caixeta (Tabebuia cassinoides) populations. Sci Forest (2000) 58:127–143.
Slatkin M. A measure of population subdivision based on microsatellite allele frequencies. Genetics (1995) 139:457–462.[Web of Science][Medline]
Streiff R, Labbé T, Bacilieri R, Steinkellner H, Gloessl J, Kremer A. Within-population genetic structure in Quercus robur L. and Quercus petraea (Matt.) Liebl. assessed with isozymes and microsatellites. Mol Ecol (1998) 7:317–328.[CrossRef]
Sun GL, Díaz O, Salomon B, von Botmer R. Genetic diversity and structure in a natural Elymus caninus population from Denmark based on microsatellite and isozyme analyses. Plant Syst Evol (2001) 227:235–244.[CrossRef]
Thelen GC, Allendorf FW. Heterozygosity-fitness correlations in rainbow trout: effects of allozyme loci or associative overdominance? Evolution (2001) 55:1180–1187.[CrossRef][Web of Science][Medline]
Volis S, Yakubov B, Shulgina I, Ward D, Mendlinger S. Distinguishing adaptive from nonadaptive genetic differentiation: comparison of QST and FST at two spatial scales. Heredity (2005) 95:466–475.[CrossRef][Web of Science][Medline]
Weir BS. Genetic data analysis II: methods for discrete population genetic data (1996) Sunderland (MA): Sinauer Associates.
Zane L, Bargelloni L, Patarnello T. Strategies for microsatellite isolation: a review. Mol Ecol (2002) 11:1–16.[CrossRef][Medline]
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
