Genetic structure of natural Eugenia dysenterica DC (Myrtaceae) populations in northeastern Goias, Brazil, accessed by morphological traits
http://www.100md.com
《遗传学和分子生物学》
IEmbrapa Trigo, Goiania, Go, Brazil
IIUniversidade Federal de Goias, Escola de Agronomia e Engenharia de Alimentos, Goiania, Go, Brazil
ABSTRACT
The Brazilian Cerrado (Savannah) is a biome with great biodiversity. Many of the species found in this biome have unlimited economic potential and are used by the local population for various purposes. With the opening of the agricultural frontiers in this region, much of the native flora has been replaced by extensive areas of monocultures. Thirteen populations of the Eugenia dysenterica DC species were sampled in the northeast of the state of Goias, Brazil, to analyze their genetic structure. Morphological traits were measured for all the progeny obtained and some of them were analyzed by using RAPD markers. The natural populations of the species showed high levels of variability for the morphological traits. The phenotypic, genotypic and genetic variabilities as structured in these populations were highly significant among populations, although most of the variability was detected within populations. The phenotypic and genetic variabilities were structured spatially as shown by the positive and highly significant correlation coefficients between the phenotypic and genetic distance matrixes and the geographic distance matrix. The ideal conservation strategy for the species in the region requires sampling of a high number of individuals of each population and also a significant number of populations to ensure suitable effective population sizes.
Key words: Cerrado, cagaiteira, population genetics, genetic structure and conservation.
Introduction
The cagaiteira (Eugenia dysenterica DC) is a species that occurs naturally in stricto sensu Cerrado areas. This species shows flowering strategy of the big-bang type, in which bloom synchronization occurs in a very short time period. Bees are the main pollinators including species of the Bombus genus type, among others. The cagaiteira as well as other species of the Myrtaceae family are dispersed by zoocory. However, the Cerrado has recently undergone a process of human disturbance and extensive areas of monoculture have replaced much of the native flora. Many of the natural species in this biome have economic potential and are used by the local population for various purposes. The domestication and incorporation of these species in the regional production systems and the design of conservation strategies depend fundamentally on knowledge of the genetic structure of their natural populations.
Different classes of markers have been used in the genetic analysis of natural populations. Among these classes the morphological traits have been used extensively in genetic analysis because they are easy to handle. Many studies have been carried out in species from the Cerrado using this type of trait (Silva et al., 2001). More recently, the application of technologies based on DNA polymorphism has provided a wealth of information on the variability of wild and cultivated populations. The molecular markers used mostly for determining the genetic variability at the DNA level for natural population analysis are RAPD markers and in some specific cases microsatellites and RFLP. These techniques have been applied to wild species to answer specific questions, mainly regarding the crossing system (Hayward and Hamilton, 1997). Some studies based on microsatellites have been employed in the study of natural populations (Collevatti et al., 2001; Zucchi et al., 2003) but most of the studies are performed with RAPD markers in plants (Chalmers et al., 1992; Wolf and Rijn, 1993, Bucci and Menozzi, 1995; Yeh et al., 1995; Maguire and Sedgley 1997; Gillies et al., 1997, Palacios and Gonzalez-Candelas, 1997; Gallois et al., 1998, Cattan-Toupance et al., 1998, Skabo et al., 1998), and in other organisms (Peever and Milgroom, 1994, Lou et al., 1998).
Eugenia dysenterica DC populations from the southeast region of Goias were studied using morphological traits (Silva et al., 2001), isoenzymatic (Telles et al., 2001 a and b) and RAPD and SSR markers (Zucchi et al., 2003). The authors detected stochastic processes of genetic divergence among populations, such as the stepping stone model. In addition to the aspects related to the biology of the species, an important factor in explaining the spatial pattern obtained was the advanced fragmentation process mainly caused by the expansion of the agricultural frontier in the Cerrado region. This process reduces the size of the cagaiteira populations and of their pollinators and dispersers, and is most certainly a contributing factor in the reduction of the genetic flow among the populations. As the fragmentation process in these areas is relatively recent, the divergence among these populations may become even more intense in the future.
The objective of the present study was to characterize the genetic structure of natural Eugenia dysenterica DC populations in the northeastern region of Goias state, Brazil, using morphological traits and RAPD markers to direct future conservation, domestication and breeding programs.
Materials and Methods
The material for the present study was obtained by field collection in the northeastern region of Goias, including 13 areas of stricto sensu Cerrado. Twelve mother plants were sampled in each area by collecting a sufficient number of fruits for characterization and to obtain progenies. In Table 1, we show the location of the counties where the samples were obtained with their respective geographic coordinates.
Morphological trait assessment
The fruits collected from each plant were characterized by measuring the fruit weight (FW), seed weight (SW) and number of seeds (NS). Later the seeds extracted from the fruits were sown in a greenhouse with 80% shading to study the variables related to the initial development of the progenies. A completely randomized block design was adopted with three replications and each plot consisted of five plants. Data were collected on percentage emergence (EP), number of days to emergence (NDE), plant height (PH90, PH178) and plant diameter (PD90, PD178) at 90 and 178 days after sowing, respectively.
RAPD analysis
Three field populations collected (1, 2 and 4) were used for the genetic characterization based on RAPD markers. The protocol for DNA extraction was adapted from Ferreira and Grattapaglia (1996). Two preliminary experiments were carried out for the qualitative and quantitative selection of the primers used. Eighty decamer primers from Operon technologies were analyzed and among them six primers of the OPA series, which showed a more polymorphic pattern, were selected for this analysis (Table 2). Amplification reactions were processed with the selected primers from the DNA samples. To compose each amplification reaction, a total volume of 25 μL of the mixture, containing distilled and deionized H2O (11.5 μL), 10 mM Tris-HCl, pH 8.3 (2.5 μL), 1.25 mM of each dNTP (2.5 μL), 0.25 μM of a decamer primer (Operon Technologies, Alameda, CA, USA), 10 ng/μL (3.0 μL), 5 mM MgCl2 (1.5 μL), 1 ud Taq DNA polymerase (1.0 μL) and 10 ng/μL of genomic DNA (3.0 μL) was used. The enzymatic amplification was performed in a thermocycler programmed initially for 5 min at 94 °C followed by 48 cycles of: 30 s at 92 °C for DNA denaturation, 1.5 min at 37 °C for annealing of the primer to the DNA template and 6 min at 72 °C for primer extension. The fragments generated by amplification were separated by electrophoresis in horizontal 1.4% agarose gel at a concentration of 1% in TBE buffer (90 mM Tris-borate and 2 mM EDTA). After separation the results were visualized by treatment with 10 μL ethidium bromide (10 mg/mL) and exposure to ultraviolet light. After photographing, the gel was interpreted for presence (1) or absence (0) of each band, for all the primers used. For the statistical analysis only the strong and reproductive bands were selected.
Analysis of the data obtained from morphological traits
Data collected from the fruits and young plants were utilized in structuring of the analysis of variance and the partitioning of variance in components. A hierarchical model which considered the effects of population, progeny within population and fruit or seedling within progeny was used. The data also allowed the construction of a phenotypic (fruit data) and genotypic (seedling data) distance matrix using the Mahalanobis generalized distance (Mahalanobis, 1936) obtained by a multivariate analysis for quantitative variables implemented by the GENES program (Cruz, 1997).
Analysis of the data obtained from RAPD markers
The coefficient of similarity by simple matching was calculated from the binary data among all the individuals to construct a dissimilarity matrix (genetic distance). The decomposition of the total variation in its components among and within populations was performed by the analysis of molecular variance (Excoffier et al., 1992) from the genetic distance matrix, using the Arlequin software (Schneider et al., 2000). The significance of the fixation indexes was tested using 1,000 permutations.
Analysis of the spatial pattern of divergence among populations
The spatial variation patterns were analyzed by the Pearson coefficient of correlation (r) between the phenotypic or genotypic average distance matrixes among populations, and the geographic distance matrix among them, using the NTSYS 1.7 software. The pattern of this variation was studied by the construction of dendrograms from the phenotypic and genotypic distance matrixes, using the UPGMA clustering method. The significances of the dendrograms were tested by cophenetic correlation with 5,000 permutations.
Results and Discussion
The results of the analysis of variance and their components for fruit and progeny data are presented in Tables 3 and 4. Fruit weight contributed most to phenotypic divergence among populations. On average, the fruit weight value was 15.34 g, indicating a relatively small fruit size and the seeds were also small with a mean of 2.19 g. The mean number of seeds per fruit was 1.46. The mean NDE values (38.12 days) and EP (81.71%) indicate a high percentage of emergence in a relatively short time period compared with other fruit tree species native to the Cerrado.
The phenotypic and genotypic structuring of the populations showed that most of the variability is found within populations for most of the variables. The occurrence of significant genetic variability among progeny within populations is expected and corresponds to 1/4 of the additive variance of the population, with the assumption of random crossing and many individuals per progeny. Deviations from random crossing due to a restricted number of individuals lead generally to an increase among progeny variation. The variation among populations occurs when there is some restriction in gene flow and also arises through the effect of genetic drift and specific adaptations to particular environments via selection. In the present case, the proportion of the genetic variance among populations, compared to the total variance was very high when compared with what is normally detected with selectively neutral markers in cross-pollinating species. This fact probably indicates an effect of selection in the differentiation among populations. However, the influence of maternal effects, which include non-genetic effects via seed vigor, cannot be disregarded.
In Table 5, we show the results of the structuring of the genetic variability with RAPD markers. The analysis of molecular variance indicated a total variance of 0.43456. This value is similar to that obtained by Gillies et al. (1997) who analyzed Cedrela odorata populations in Costa Rica, with the same type of marker. High levels of total variance are also reported for temperate climate species (Yeh et al., 1995; Liao and Hsiao, 1998). Although most of the genetic variance is found within populations there is a highly significant quantity among populations (about 8% of total variability), indicating gene flow restriction among populations. Liao and Hsiao (1998) used RAPD markers to study the Acorus gramineus Soland (Aracae) species that occurs in Southeast Asia and detected levels of differentiation (16.88%) among populations. According to these authors, the presence of specific bands at each sampling site may be the result of the separate evolution of the populations by restriction to the gene flow. Lower levels of differentiation among populations were obtained by Yeh et al. (1995) in natural populations of the Populus tremuloides Michx species. The differentiation among populations, although small (2.6%), was highly significant. This level of differentiation among populations was also observed in other studies with the same species using isoenzymes markers. These results suggest that neutral forces are a better explanation of the structure observed for selection. When the differentiation values among populations for the morphological traits are compared, either at phenotypic or genotypic level, the lower value observed for RAPD markers represents the role of genetic drift in the differentiation of the population, because these markers are selectively neutral they are not specific targets of the natural selection process. The differentiation values observed with progeny data quantify the role of drift and the selection that may explain the greater proportions of variability observed for most of the phenotypic variables. The analysis of data obtained from fruit variables showed that besides these two sources (drift and selection) the environment exercises a strong influence.
Zucchi et al. (2003) analyzed Eugenia dysenterica DC populations from southeastern Goias with the same type of marker and obtained a value of fST = 0.2703 and a genetic structure that concentrated 27.03% of the genetic variability among populations. The lower values detected in the present study might have occurred because of the greater proximity among the populations studied and also because of the smaller degree of fragmentation existing in the northeastern region that so far has not been exploited for agricultural purposes. In this context, these populations may have been in a much slower population differentiation process than those found in the southeastern region of the state. Analyzing the same populations with microsatellite markers (SSR), Zucchi et al. (2003) ascertained that the species came close to alogamy and the allelic frequencies of the populations were detected in the proportions of Hardy-Weinberg equilibrium for the generation sampled, in contrast with the results obtained with isoenzymes by Telles et al. (2001 a and b) suggesting panmixy deviations. However, according to Zucchi et al. (2003) there is congruence among the population parameters obtained with different types of markers (isoenzymes, RAPD and SSR markers) despite some discrepancies. The results obtained by these studies in the southeastern region of the state show that in spite of the species having a tendency to alogamy, the gene flow between populations is small and considered restrictive, probably because of the human disturbance of the Cerrado. The restricted flow allowed the drift to increase the differentiation among populations.
The comparison of population parameters obtained with different markers should be interpreted with caution. One restriction that should be considered in this type of comparison is the sampling of the genome itself. RAPD and SSR markers affect coding as well as non-coding regions of the genome and the isoenzymes are restricted to the coding region. Due to their dominant nature, RAPD markers can lead to an underestimation of the recessive allele frequency in a population causing a bias in the estimates of the genetic parameters.
The dendrogram patterns (Figure 1) indicate an almost continuous variation of the phenotypic and genotypic distances since no characteristic clusters were formed. This performance, reflected in the dendrograms, is expected under the null hypothesis. Those populations that are closer geographically are more similar phenotypically and genotypically and there is a clinal variation pattern. This pattern, which reflects the evolutionary model of isolation by distance, was observed by Telles et al. (2001 a and b) and can also be visualized in Figure 2. In Figure 2, we show that there was a clinal variation pattern because small geographic distances result in small phenotypic and genotypic distances. However, for greater geographic distances there are both small and large phenotypic and genotypic variations (envelope effect), indicating that the differentiation process becomes random. Therefore there should be a point (distance in km) at which the phenotypic and genotypic variation becomes independent from the geographic distances (indicating isolation by distance). In the present case, no evidence was observed of stabilization in the variability among populations, indicating that greater phenotypic and genotypic variations could be obtained if collections were made in wider regions. Telles et al. (2001a) confirmed this clinal variation pattern by analysis of spatial autocorrelation in populations of the species in the southeastern part of the state.
The genetic variability pattern showed a space structuring by the positive and highly significant coefficient of correlation between the phenotypic or genotypic variation matrixes with the geographic distance matrix (Figure 3). This fact reinforces the supposition that the genetic variability in these populations is structured in geographic space, confirming that this variability is best explained by a neutral evolution model (Telles et al., 2001b).
To establish in situ conservation or sampling areas for ex situ conservation the concepts of representativeness (Vencovsky and Crossa, 1999) and the minimum viable population (Menges, 2000a and b) should be considered. These approaches are important for the sustainability of the populations being preserved. The levels of variability among populations observed in this study indicate the need for a large number of areas for in situ conservation and an adequate number of individuals from each population to ensure appropriate effective sizes.
Acknowledgements
The authors wish to thank the Graduate Program of the Agronomy School of the Federal University of Goias, FINEP and CAPES for financial support.
References
Bucci G and Menozzi P (1995) Genetic variation of RAPD markers in a Picea abies Karst. Population. Heredity 75:188-197.
Cattan-Toupance I, Michalakis Y and Neema C (1998) Genetic structure of wild bean population in their South-Andean centre of origin. Theor Appl Genet 96:844-851.
Chalmers KJ, Waugh R, Sprent JL, Simons AJ and Powell W (1992) Detection of genetic variation between and within populations of Gliricidia sepium and G. maculata using RAPD markers. Heredity 69:465-472.
Collevatti RG, Gratapaglia D and Hay JD (2001) Population genetic structure of the endangered tropical tree species Caryocar brasiliense, based on variability at microsatellite loci. Molecular Ecology 10:349-356.
Cruz CD (1997) Programa GENES. Aplicativo Computacional em Genetica e Estatistica. Editora UFV, Viosa, 590 pp.
Excoffier L, Smouses PE and Quattro JM (1992) Analysis of molecular variance inferred from metric distances among DNA Haplotypes: Application to human mitocondrial DNA restriction data. Genetics 131:479-491.
Ferreira ME and Grattapaglia D (1996) Introduo ao uso de marcadores moleculares em analise genetica. 2nd edition. Embrapa, Brasilia, 220 pp.
Gallois A, Audran JC and Burrus M (1998) Assessment of genetic relationship and population discrimination among fagus sylvatica L. by RAPD. Theor Appl Genet 97:211-219.
Gillies ACM, Cornelius JP, Newton AC, Navarro C and Hernandez M (1997) Genetic variation in Costa Rican populations of the tropical timber species Cedrela odorata L. assessed using RAPDs. Mol Ecol 6:1133-1145.
Hayward MD and Hamilton NRS (1997) Genetic diversity, population structure and conservation. In: Callow JA, Ford-Loyd BV and Newbury HJ (eds) Biotechnology and Plant Genetic Resources - Conservation and Use. Biotechnology in agriculture series n. 19, CAB International, Wallingford, Oxon, UK, pp 49-76.
Lou KF, Weiss MJ, Bruckner PL, Morril WL, Tabert LE and Martin JM (1998) RAPD variation within and among geographic populations of wheat stem sawfly (Cephus cinctus Norton). The American Genetic Association 89:329-335.
Liao LC and Hsiao JY (1998) Relationship between population genetic structure and riparian habitat as revealed by RAPD analysis of the rheophyte Acorus gramineus Soland. (Araceae) in Taiwan. Mol Ecol 7:1275-1281.
Maguire TL and Sedgley M (1997) Genetic diversity in Banksia and Dryandra (Protaceae) with emphasis on Banksia cuneata, a rare and endangered species. Heredity 79:394-401.
Mahalanobis PC (1936) On the generalized distance in statistics. Proc Natl Inst Sci 2:49-55.
Menges ES (2000a) Population viability analyses in plants. Challenges and opportunities. Tree 15:51-56.
Menges ES (2000b) Applications of population viability analyses in plant. Ecological Bulletins 48:73-84.
Palacios C and Gonzalez-Candelas F (1997) Analysis of population genetic structure and variability using RAPD markers in the endemic and endangered Limonium dufourii (Plumbaginaceae). Mol Ecol 6:1107-1121.
Peever TL and Milgroom MG (1994) Genetic structure of Pyrenophora teres population determined with random amplified polymorphic DNA markers. Can J Bot 72:915-923.
Scheneider S, Roessli D and Excoffier L (2000) Alerquin ver. 2000. A Software for Population Data Analysis. Genetic and Biometry Laboratory, University of Geneva, Switzerland.
Silva RSM, Chaves LJ and Naves RV (2001) Caracterizao de frutos e arvores de cagaita (Eugenia dysenterica DC) no sudeste de Goias, Brasil. Revista Brasileira de Fruticultura 23:320-334.
Skabo S, Vaillancourt RE and Potts, BM (1998) Fine-scale genetic structure of Eucalyptus globulus Spp. Globulus forest revealed by RAPDs. Aust J Bot 46:583-594.
Telles MPC, Diniz-Filho JAF, Coelho ASG and Chaves LJ (2001a) Autocorrelao espacial das freqvências alelicas em subpopulaes de cagaiteira (Eugenia dysenterica DC, Myrtaceae) no sudeste de Goias. Revista Brasileira de Botanica 24:145-154.
Telles MPC, Silva RSM, Chaves LJ, Coelho ASG and Diniz-Filho JAF (2001b) Divergência entre subpopulaes de cagaiteira (Eugenia dysenterica) em resposta a padres edaficos e distribuio espacial. Pesquisa Agropecuaria Brasileira 36:1387-1394.
Vencovsky R and Crossa J (1999) Measures of representativeness. In: Borem A, Giudice MP and Sakiyama NS (eds) Plant Breeding in the Turn of the Millennium. Editora da UFV, Viosa, pp 335-354.
Zucchi MI, Brondani RPV, Pinheiro JB, Chaves LJ, Coelho ASG and Vencovsky R (2003) Genetic structure and gene flow in Eugenia dysenterica DC in the Brazilian Cerrado utilizing SSR makers. Genetics and Molecular Biology 26:449-457.
Wolf FC and Rijn JPV (1993) Rapid detection of genetic variability in Chrysanthemum (Dendrantema grandiflora Tzvelev) using random primers. Heredity 71:335-341.
Yeh FC, Chong KX and Yang RC (1995) RAPD variation within and among natural populations of Trembling Aspen (Populus tremuloides Michx) from Alberta. Journal of Heredity 86:454-460.(Maria da Gloria Trindade;)
IIUniversidade Federal de Goias, Escola de Agronomia e Engenharia de Alimentos, Goiania, Go, Brazil
ABSTRACT
The Brazilian Cerrado (Savannah) is a biome with great biodiversity. Many of the species found in this biome have unlimited economic potential and are used by the local population for various purposes. With the opening of the agricultural frontiers in this region, much of the native flora has been replaced by extensive areas of monocultures. Thirteen populations of the Eugenia dysenterica DC species were sampled in the northeast of the state of Goias, Brazil, to analyze their genetic structure. Morphological traits were measured for all the progeny obtained and some of them were analyzed by using RAPD markers. The natural populations of the species showed high levels of variability for the morphological traits. The phenotypic, genotypic and genetic variabilities as structured in these populations were highly significant among populations, although most of the variability was detected within populations. The phenotypic and genetic variabilities were structured spatially as shown by the positive and highly significant correlation coefficients between the phenotypic and genetic distance matrixes and the geographic distance matrix. The ideal conservation strategy for the species in the region requires sampling of a high number of individuals of each population and also a significant number of populations to ensure suitable effective population sizes.
Key words: Cerrado, cagaiteira, population genetics, genetic structure and conservation.
Introduction
The cagaiteira (Eugenia dysenterica DC) is a species that occurs naturally in stricto sensu Cerrado areas. This species shows flowering strategy of the big-bang type, in which bloom synchronization occurs in a very short time period. Bees are the main pollinators including species of the Bombus genus type, among others. The cagaiteira as well as other species of the Myrtaceae family are dispersed by zoocory. However, the Cerrado has recently undergone a process of human disturbance and extensive areas of monoculture have replaced much of the native flora. Many of the natural species in this biome have economic potential and are used by the local population for various purposes. The domestication and incorporation of these species in the regional production systems and the design of conservation strategies depend fundamentally on knowledge of the genetic structure of their natural populations.
Different classes of markers have been used in the genetic analysis of natural populations. Among these classes the morphological traits have been used extensively in genetic analysis because they are easy to handle. Many studies have been carried out in species from the Cerrado using this type of trait (Silva et al., 2001). More recently, the application of technologies based on DNA polymorphism has provided a wealth of information on the variability of wild and cultivated populations. The molecular markers used mostly for determining the genetic variability at the DNA level for natural population analysis are RAPD markers and in some specific cases microsatellites and RFLP. These techniques have been applied to wild species to answer specific questions, mainly regarding the crossing system (Hayward and Hamilton, 1997). Some studies based on microsatellites have been employed in the study of natural populations (Collevatti et al., 2001; Zucchi et al., 2003) but most of the studies are performed with RAPD markers in plants (Chalmers et al., 1992; Wolf and Rijn, 1993, Bucci and Menozzi, 1995; Yeh et al., 1995; Maguire and Sedgley 1997; Gillies et al., 1997, Palacios and Gonzalez-Candelas, 1997; Gallois et al., 1998, Cattan-Toupance et al., 1998, Skabo et al., 1998), and in other organisms (Peever and Milgroom, 1994, Lou et al., 1998).
Eugenia dysenterica DC populations from the southeast region of Goias were studied using morphological traits (Silva et al., 2001), isoenzymatic (Telles et al., 2001 a and b) and RAPD and SSR markers (Zucchi et al., 2003). The authors detected stochastic processes of genetic divergence among populations, such as the stepping stone model. In addition to the aspects related to the biology of the species, an important factor in explaining the spatial pattern obtained was the advanced fragmentation process mainly caused by the expansion of the agricultural frontier in the Cerrado region. This process reduces the size of the cagaiteira populations and of their pollinators and dispersers, and is most certainly a contributing factor in the reduction of the genetic flow among the populations. As the fragmentation process in these areas is relatively recent, the divergence among these populations may become even more intense in the future.
The objective of the present study was to characterize the genetic structure of natural Eugenia dysenterica DC populations in the northeastern region of Goias state, Brazil, using morphological traits and RAPD markers to direct future conservation, domestication and breeding programs.
Materials and Methods
The material for the present study was obtained by field collection in the northeastern region of Goias, including 13 areas of stricto sensu Cerrado. Twelve mother plants were sampled in each area by collecting a sufficient number of fruits for characterization and to obtain progenies. In Table 1, we show the location of the counties where the samples were obtained with their respective geographic coordinates.
Morphological trait assessment
The fruits collected from each plant were characterized by measuring the fruit weight (FW), seed weight (SW) and number of seeds (NS). Later the seeds extracted from the fruits were sown in a greenhouse with 80% shading to study the variables related to the initial development of the progenies. A completely randomized block design was adopted with three replications and each plot consisted of five plants. Data were collected on percentage emergence (EP), number of days to emergence (NDE), plant height (PH90, PH178) and plant diameter (PD90, PD178) at 90 and 178 days after sowing, respectively.
RAPD analysis
Three field populations collected (1, 2 and 4) were used for the genetic characterization based on RAPD markers. The protocol for DNA extraction was adapted from Ferreira and Grattapaglia (1996). Two preliminary experiments were carried out for the qualitative and quantitative selection of the primers used. Eighty decamer primers from Operon technologies were analyzed and among them six primers of the OPA series, which showed a more polymorphic pattern, were selected for this analysis (Table 2). Amplification reactions were processed with the selected primers from the DNA samples. To compose each amplification reaction, a total volume of 25 μL of the mixture, containing distilled and deionized H2O (11.5 μL), 10 mM Tris-HCl, pH 8.3 (2.5 μL), 1.25 mM of each dNTP (2.5 μL), 0.25 μM of a decamer primer (Operon Technologies, Alameda, CA, USA), 10 ng/μL (3.0 μL), 5 mM MgCl2 (1.5 μL), 1 ud Taq DNA polymerase (1.0 μL) and 10 ng/μL of genomic DNA (3.0 μL) was used. The enzymatic amplification was performed in a thermocycler programmed initially for 5 min at 94 °C followed by 48 cycles of: 30 s at 92 °C for DNA denaturation, 1.5 min at 37 °C for annealing of the primer to the DNA template and 6 min at 72 °C for primer extension. The fragments generated by amplification were separated by electrophoresis in horizontal 1.4% agarose gel at a concentration of 1% in TBE buffer (90 mM Tris-borate and 2 mM EDTA). After separation the results were visualized by treatment with 10 μL ethidium bromide (10 mg/mL) and exposure to ultraviolet light. After photographing, the gel was interpreted for presence (1) or absence (0) of each band, for all the primers used. For the statistical analysis only the strong and reproductive bands were selected.
Analysis of the data obtained from morphological traits
Data collected from the fruits and young plants were utilized in structuring of the analysis of variance and the partitioning of variance in components. A hierarchical model which considered the effects of population, progeny within population and fruit or seedling within progeny was used. The data also allowed the construction of a phenotypic (fruit data) and genotypic (seedling data) distance matrix using the Mahalanobis generalized distance (Mahalanobis, 1936) obtained by a multivariate analysis for quantitative variables implemented by the GENES program (Cruz, 1997).
Analysis of the data obtained from RAPD markers
The coefficient of similarity by simple matching was calculated from the binary data among all the individuals to construct a dissimilarity matrix (genetic distance). The decomposition of the total variation in its components among and within populations was performed by the analysis of molecular variance (Excoffier et al., 1992) from the genetic distance matrix, using the Arlequin software (Schneider et al., 2000). The significance of the fixation indexes was tested using 1,000 permutations.
Analysis of the spatial pattern of divergence among populations
The spatial variation patterns were analyzed by the Pearson coefficient of correlation (r) between the phenotypic or genotypic average distance matrixes among populations, and the geographic distance matrix among them, using the NTSYS 1.7 software. The pattern of this variation was studied by the construction of dendrograms from the phenotypic and genotypic distance matrixes, using the UPGMA clustering method. The significances of the dendrograms were tested by cophenetic correlation with 5,000 permutations.
Results and Discussion
The results of the analysis of variance and their components for fruit and progeny data are presented in Tables 3 and 4. Fruit weight contributed most to phenotypic divergence among populations. On average, the fruit weight value was 15.34 g, indicating a relatively small fruit size and the seeds were also small with a mean of 2.19 g. The mean number of seeds per fruit was 1.46. The mean NDE values (38.12 days) and EP (81.71%) indicate a high percentage of emergence in a relatively short time period compared with other fruit tree species native to the Cerrado.
The phenotypic and genotypic structuring of the populations showed that most of the variability is found within populations for most of the variables. The occurrence of significant genetic variability among progeny within populations is expected and corresponds to 1/4 of the additive variance of the population, with the assumption of random crossing and many individuals per progeny. Deviations from random crossing due to a restricted number of individuals lead generally to an increase among progeny variation. The variation among populations occurs when there is some restriction in gene flow and also arises through the effect of genetic drift and specific adaptations to particular environments via selection. In the present case, the proportion of the genetic variance among populations, compared to the total variance was very high when compared with what is normally detected with selectively neutral markers in cross-pollinating species. This fact probably indicates an effect of selection in the differentiation among populations. However, the influence of maternal effects, which include non-genetic effects via seed vigor, cannot be disregarded.
In Table 5, we show the results of the structuring of the genetic variability with RAPD markers. The analysis of molecular variance indicated a total variance of 0.43456. This value is similar to that obtained by Gillies et al. (1997) who analyzed Cedrela odorata populations in Costa Rica, with the same type of marker. High levels of total variance are also reported for temperate climate species (Yeh et al., 1995; Liao and Hsiao, 1998). Although most of the genetic variance is found within populations there is a highly significant quantity among populations (about 8% of total variability), indicating gene flow restriction among populations. Liao and Hsiao (1998) used RAPD markers to study the Acorus gramineus Soland (Aracae) species that occurs in Southeast Asia and detected levels of differentiation (16.88%) among populations. According to these authors, the presence of specific bands at each sampling site may be the result of the separate evolution of the populations by restriction to the gene flow. Lower levels of differentiation among populations were obtained by Yeh et al. (1995) in natural populations of the Populus tremuloides Michx species. The differentiation among populations, although small (2.6%), was highly significant. This level of differentiation among populations was also observed in other studies with the same species using isoenzymes markers. These results suggest that neutral forces are a better explanation of the structure observed for selection. When the differentiation values among populations for the morphological traits are compared, either at phenotypic or genotypic level, the lower value observed for RAPD markers represents the role of genetic drift in the differentiation of the population, because these markers are selectively neutral they are not specific targets of the natural selection process. The differentiation values observed with progeny data quantify the role of drift and the selection that may explain the greater proportions of variability observed for most of the phenotypic variables. The analysis of data obtained from fruit variables showed that besides these two sources (drift and selection) the environment exercises a strong influence.
Zucchi et al. (2003) analyzed Eugenia dysenterica DC populations from southeastern Goias with the same type of marker and obtained a value of fST = 0.2703 and a genetic structure that concentrated 27.03% of the genetic variability among populations. The lower values detected in the present study might have occurred because of the greater proximity among the populations studied and also because of the smaller degree of fragmentation existing in the northeastern region that so far has not been exploited for agricultural purposes. In this context, these populations may have been in a much slower population differentiation process than those found in the southeastern region of the state. Analyzing the same populations with microsatellite markers (SSR), Zucchi et al. (2003) ascertained that the species came close to alogamy and the allelic frequencies of the populations were detected in the proportions of Hardy-Weinberg equilibrium for the generation sampled, in contrast with the results obtained with isoenzymes by Telles et al. (2001 a and b) suggesting panmixy deviations. However, according to Zucchi et al. (2003) there is congruence among the population parameters obtained with different types of markers (isoenzymes, RAPD and SSR markers) despite some discrepancies. The results obtained by these studies in the southeastern region of the state show that in spite of the species having a tendency to alogamy, the gene flow between populations is small and considered restrictive, probably because of the human disturbance of the Cerrado. The restricted flow allowed the drift to increase the differentiation among populations.
The comparison of population parameters obtained with different markers should be interpreted with caution. One restriction that should be considered in this type of comparison is the sampling of the genome itself. RAPD and SSR markers affect coding as well as non-coding regions of the genome and the isoenzymes are restricted to the coding region. Due to their dominant nature, RAPD markers can lead to an underestimation of the recessive allele frequency in a population causing a bias in the estimates of the genetic parameters.
The dendrogram patterns (Figure 1) indicate an almost continuous variation of the phenotypic and genotypic distances since no characteristic clusters were formed. This performance, reflected in the dendrograms, is expected under the null hypothesis. Those populations that are closer geographically are more similar phenotypically and genotypically and there is a clinal variation pattern. This pattern, which reflects the evolutionary model of isolation by distance, was observed by Telles et al. (2001 a and b) and can also be visualized in Figure 2. In Figure 2, we show that there was a clinal variation pattern because small geographic distances result in small phenotypic and genotypic distances. However, for greater geographic distances there are both small and large phenotypic and genotypic variations (envelope effect), indicating that the differentiation process becomes random. Therefore there should be a point (distance in km) at which the phenotypic and genotypic variation becomes independent from the geographic distances (indicating isolation by distance). In the present case, no evidence was observed of stabilization in the variability among populations, indicating that greater phenotypic and genotypic variations could be obtained if collections were made in wider regions. Telles et al. (2001a) confirmed this clinal variation pattern by analysis of spatial autocorrelation in populations of the species in the southeastern part of the state.
The genetic variability pattern showed a space structuring by the positive and highly significant coefficient of correlation between the phenotypic or genotypic variation matrixes with the geographic distance matrix (Figure 3). This fact reinforces the supposition that the genetic variability in these populations is structured in geographic space, confirming that this variability is best explained by a neutral evolution model (Telles et al., 2001b).
To establish in situ conservation or sampling areas for ex situ conservation the concepts of representativeness (Vencovsky and Crossa, 1999) and the minimum viable population (Menges, 2000a and b) should be considered. These approaches are important for the sustainability of the populations being preserved. The levels of variability among populations observed in this study indicate the need for a large number of areas for in situ conservation and an adequate number of individuals from each population to ensure appropriate effective sizes.
Acknowledgements
The authors wish to thank the Graduate Program of the Agronomy School of the Federal University of Goias, FINEP and CAPES for financial support.
References
Bucci G and Menozzi P (1995) Genetic variation of RAPD markers in a Picea abies Karst. Population. Heredity 75:188-197.
Cattan-Toupance I, Michalakis Y and Neema C (1998) Genetic structure of wild bean population in their South-Andean centre of origin. Theor Appl Genet 96:844-851.
Chalmers KJ, Waugh R, Sprent JL, Simons AJ and Powell W (1992) Detection of genetic variation between and within populations of Gliricidia sepium and G. maculata using RAPD markers. Heredity 69:465-472.
Collevatti RG, Gratapaglia D and Hay JD (2001) Population genetic structure of the endangered tropical tree species Caryocar brasiliense, based on variability at microsatellite loci. Molecular Ecology 10:349-356.
Cruz CD (1997) Programa GENES. Aplicativo Computacional em Genetica e Estatistica. Editora UFV, Viosa, 590 pp.
Excoffier L, Smouses PE and Quattro JM (1992) Analysis of molecular variance inferred from metric distances among DNA Haplotypes: Application to human mitocondrial DNA restriction data. Genetics 131:479-491.
Ferreira ME and Grattapaglia D (1996) Introduo ao uso de marcadores moleculares em analise genetica. 2nd edition. Embrapa, Brasilia, 220 pp.
Gallois A, Audran JC and Burrus M (1998) Assessment of genetic relationship and population discrimination among fagus sylvatica L. by RAPD. Theor Appl Genet 97:211-219.
Gillies ACM, Cornelius JP, Newton AC, Navarro C and Hernandez M (1997) Genetic variation in Costa Rican populations of the tropical timber species Cedrela odorata L. assessed using RAPDs. Mol Ecol 6:1133-1145.
Hayward MD and Hamilton NRS (1997) Genetic diversity, population structure and conservation. In: Callow JA, Ford-Loyd BV and Newbury HJ (eds) Biotechnology and Plant Genetic Resources - Conservation and Use. Biotechnology in agriculture series n. 19, CAB International, Wallingford, Oxon, UK, pp 49-76.
Lou KF, Weiss MJ, Bruckner PL, Morril WL, Tabert LE and Martin JM (1998) RAPD variation within and among geographic populations of wheat stem sawfly (Cephus cinctus Norton). The American Genetic Association 89:329-335.
Liao LC and Hsiao JY (1998) Relationship between population genetic structure and riparian habitat as revealed by RAPD analysis of the rheophyte Acorus gramineus Soland. (Araceae) in Taiwan. Mol Ecol 7:1275-1281.
Maguire TL and Sedgley M (1997) Genetic diversity in Banksia and Dryandra (Protaceae) with emphasis on Banksia cuneata, a rare and endangered species. Heredity 79:394-401.
Mahalanobis PC (1936) On the generalized distance in statistics. Proc Natl Inst Sci 2:49-55.
Menges ES (2000a) Population viability analyses in plants. Challenges and opportunities. Tree 15:51-56.
Menges ES (2000b) Applications of population viability analyses in plant. Ecological Bulletins 48:73-84.
Palacios C and Gonzalez-Candelas F (1997) Analysis of population genetic structure and variability using RAPD markers in the endemic and endangered Limonium dufourii (Plumbaginaceae). Mol Ecol 6:1107-1121.
Peever TL and Milgroom MG (1994) Genetic structure of Pyrenophora teres population determined with random amplified polymorphic DNA markers. Can J Bot 72:915-923.
Scheneider S, Roessli D and Excoffier L (2000) Alerquin ver. 2000. A Software for Population Data Analysis. Genetic and Biometry Laboratory, University of Geneva, Switzerland.
Silva RSM, Chaves LJ and Naves RV (2001) Caracterizao de frutos e arvores de cagaita (Eugenia dysenterica DC) no sudeste de Goias, Brasil. Revista Brasileira de Fruticultura 23:320-334.
Skabo S, Vaillancourt RE and Potts, BM (1998) Fine-scale genetic structure of Eucalyptus globulus Spp. Globulus forest revealed by RAPDs. Aust J Bot 46:583-594.
Telles MPC, Diniz-Filho JAF, Coelho ASG and Chaves LJ (2001a) Autocorrelao espacial das freqvências alelicas em subpopulaes de cagaiteira (Eugenia dysenterica DC, Myrtaceae) no sudeste de Goias. Revista Brasileira de Botanica 24:145-154.
Telles MPC, Silva RSM, Chaves LJ, Coelho ASG and Diniz-Filho JAF (2001b) Divergência entre subpopulaes de cagaiteira (Eugenia dysenterica) em resposta a padres edaficos e distribuio espacial. Pesquisa Agropecuaria Brasileira 36:1387-1394.
Vencovsky R and Crossa J (1999) Measures of representativeness. In: Borem A, Giudice MP and Sakiyama NS (eds) Plant Breeding in the Turn of the Millennium. Editora da UFV, Viosa, pp 335-354.
Zucchi MI, Brondani RPV, Pinheiro JB, Chaves LJ, Coelho ASG and Vencovsky R (2003) Genetic structure and gene flow in Eugenia dysenterica DC in the Brazilian Cerrado utilizing SSR makers. Genetics and Molecular Biology 26:449-457.
Wolf FC and Rijn JPV (1993) Rapid detection of genetic variability in Chrysanthemum (Dendrantema grandiflora Tzvelev) using random primers. Heredity 71:335-341.
Yeh FC, Chong KX and Yang RC (1995) RAPD variation within and among natural populations of Trembling Aspen (Populus tremuloides Michx) from Alberta. Journal of Heredity 86:454-460.(Maria da Gloria Trindade;)