versión On-line ISSN 0718-1620
Cienc. Inv. Agr. vol.38 no.3 Santiago dic. 2011
Cien. Inv. Agr. 38(3):441-452. 2011
PLANT BREEDING, GENETICS AND GENETIC RESOURCES
Heterogeneous genetic structure in a natural population of Raulí(Nothofagus nervosa)
Heterogeneidad de la estructura genética, en una población natural de Raulí (Nothofagus nervosa)
Basilio Carrasco1, Lafayette Eaton2, Luis Letelier3, Carole Díaz1, and Rolando García-Gonzáles3
1Facultad de Agronomía e Ingeniería Forestal. Pontificia Universidad Católica de Chile. Vicuña Mackenna 4860, Macul. Casilla 306, Santiago, Chile.
2Facultad de Ciencias, Departamento de Ecología, Universidad de Chile. Las Palmeras 3425, Ñuñoa Santiago, Chile.
3Laboratorio de Biotecnología de los Recursos Naturales, Facultad de Ciencias Agrarias y Forestales. Universidad Católica del Maule, Avenida San Miguel 3650, Talca, Chile.
Heterozygote deficiencies in natural populations of outbreeding tree species are common and thought to be due mainly to biparental inbreeding. Inbreeding is believed to be caused by family structure within populations, a product of limited seed dispersal and probably limited pollen dispersal. Although both theory and simulation studies predict that structure should be apparent where trees are isolated by distance, most studies of structure in natural populations have detected only a weak spatial genetic structuring. In this contribution, we compare the use of spatial autocorrelation methodology and F statistics with the concept of relatedness to examine the spatial genetic structure in the natural population of a native southern beech and to explore the discrepancy between theory and observations. Autocorrelation detected structure in only a few of the nine enzyme loci tested in an estimated patch size of approximately 10 m. By successively eliminating the largest distances in the Gabriel map, the population was separated into groups or patches of neighbors, which were then tested for relatedness. Three groups of relatives were found interspersed with seven groups of unrelated individuals. The F statistics for these groups also showed weak genetic structure. We suggest that heterogeneity of family structure within natural populations may be one reason why more spatial genetic structure has not been detected.
Key words: Genetic structure, spatial structure, Raulí, spatial autocorrelation.
Es común encontrar deficiencias de heterocigotos en poblaciones naturales de especies forestales alógamas, lo cual se explicaría por un aumento del nivel de consanguinidad. Se ha postulado que la consanguinidad es causada por la presencia de estructuras de individuos emparentados dentro de las poblaciones, como producto de la limitada dispersión de polen y semillas. Si bien, antecedentes teóricos y estudios de simulaciones predicen que la estructura debería ser aparente cuando existe aislamiento por distancia, muchos estudios de estructura en poblaciones naturales, han detectado sólo un débil estructuración espacial. En este artículo se compara el uso de la metodología de autocorrelación espacial, estadísticos F y estimadores de parentesco. El objetivo es examinar la estructura genética espacial de una especial nativa de Nothofagus y explorar la discrepancia entre teoría y observaciones experimentales. El análisis de autocorrelación detectó estructura sólo para algunas de las nueve enzimas analizadas. Se estimaron tamaños de vecindades de alrededor de 10 m. Al eliminar sucesivamente las distancias más grandes, en el mapa de Gabriel, las poblaciones se separaron en vecindades, las que fueron analizadas para su nivel de parentesco. Se detectaron tres grupos de individuos emparentados, mezclados con siete grupos de individuos relativamente no emparentados. El estadístico F para los grupos identificados también mostró una débil estructura genética. Se sugiere que la heterogeneidad de la estructura familiar dentro de las poblaciones naturales puede ser una de las razones que explica la escasa estructura genética espacial observada en Nothofagus nervosa.
Palabras clave: Autocorrelación espacial, estructura especial, estructura genética, raulí.
Most temperate and tropical trees have reproductive systems that favor outcrossing (Seltman et al, 2009; Muona, 1990; Bawa et al, 1985). Surprisingly, many tree species with these characteristics show a deficiency of heterozygotes for allozyme loci in their natural populations (Brown, 1979). Hypotheses to explain this phenomenon, which make up half of Brown's (1979) "heterozygosity paradox," postulate a genetic structure within populations caused by a limited dispersion of pollen and/or seeds.
The limited movement of pollen restricts neighborhood size, while limited seed dispersal produces a family or kinship structure (Wright, 1943; Malecót, 1948; Levin and Kerster, 1974). Especially if both processes operate, the continuous model of isolation by distance (Wright, 1943; Malecót, 1948) predicts that neighboring individuals should tend to be relatives and should be more similar genetically than a random pair from the population.
It is generally accepted that most plant species do not form panmictic populations (Levin and Kerster, 1974; Brown, 1979) even when their distribution is continuous. There has been considerable interest in methods to document and measure genetic structure within populations, concentrated in spatial autocorrelation analysis (Sokal and Oden, 1991; Hyewood, 1991; Epperson, 1993; Epperson, 2000; Epperson, 2003) but also including F statistics (Slatkin and Arter, 1991; Merzeau et al., 1994; Hardy and Veckemans, 2002, Suarez et al., 2008) and multiple regressions (Furnier et al., 1987). Although there has been some controversy about the efficiency of spatial autocorrelation for detecting structure (Slatkin and Arter, 1991; Sokal and Oden, 1991; Bruno et al., 2008), it has been indicated that both join-counts and Moran's I should accurately detect small-scale structure under an isolation-by-distance model (Sokal and Jacquez, 1978; Epperson and Li, 1996, 1997; Bruno et al., 2008).
Quite a few northern-hemisphere species of the Fagaceae have been examined for spatial structure. In most cases, although the studied populations present a deficiency of heterozygotes, only a weak spatial genetic structure has been detected. A number of explanations have been proposed for the lack of structure, including the effects of neutral and selected loci (Bacilieri et al., 1994; Gapare and Aitken, 2005), "low heterozygosity alleles" (Leonardi and Menozzi, 1996; Belletti et al., 2005), null alleles and allele frequencies (Bacilieri et al, 1994; Aldrich et al., 2005), a small number of generations, stochastic effects (Bacilieri et al., 1994; Lopez-Aljorna et al., 2007), local differences in flowering times (Gregorius et al., 1986; Sampson et al, 1990; Hendry and Day, 2005; Sampson and Byrne, 2008) and seed dispersal (Berg and Hamrick, 1994; Kuss et al., 2008; Sethi et al., 2009).
Raulí (Nothofagus nervosa) show a considerable heterozygote deficit within its natural populations, more so than the majority of the species of Fagaceae (Carrasco, 1998; Marchelli and Gallo, 2000). Thus, we should expect spatial genetic structure to be particularly evident within the populations of this southern beech. In this contribution, we use spatial autocorrelation, F statistics and relat-edness coefficients to examine the nature of the spatial genetic structure of a small population of N. nervosa with an inbreeding coefficient of almost 0.2. We will argue that the heterogeneity of family structure within populations may be an important reason why more spatial genetic structure is not detected in natural populations of plant species.
Materials and methods
Nothofagus nervosa grows in Chile on the lower slopes of the Andes from 36°S to 40°S and in the Coast Range from 38°S to 41°S at altitudes from 100 m to 1200 m (Donoso, 1995). We mapped the locations of all 110 adult individuals (dbh > 25 cm, height > 5 cm) of a population of N. nervosa located 4 km west of Recinto, in the Bio-Bio Region of Chile (71° 41'W, 36° 49'S, 750 m altitude). In contrast to its Northern Hemisphere counterpart Fagus sylvatica L., Raulí is never found in pure stands; this population was mixed with Roble [N. obliqua (Mirb.) Oerst.], Geviuna avellana Molina and Sophora microphylla Meyen (Pollman, 2003, 2005). The stand is second growth; we estimated that it was cut 30-40 years ago. The majority of the individuals originated from stump sprouts, but perhaps one third were less than 40 years of age.
We collected a terminal lateral branch including at least 6 leaves from each individual. The bases of the branches were wrapped in wet paper towels and kept on ice in a cooler until returned to the laboratory, after which the branches were kept at 4°C with the bases in water until processed.
We ground 3-5 of the newest leaves that had been stripped of their main veins, along with scrapings of the cortex, in approximately 2 ml of an extraction buffer, pH 7.5, containing Trizma base 6.5 g, citric acid 1.5 g, cisteine 1.0 g, ascorbic acid 1.0 g, polyethylene glycol 5.0 g, 2-mercaptoethanol 2 drops, soluble Polyvinylpyrrolidone (PVP)13 g in 1 L of water. The resultant slurry was centrifuged at 5000 rpm for 4 min, and then the supernatant was poured into Eppendorf tubes in duplicate and stored at -80°C until used for electrophoresis.
The electrophoretic and staining procedures followed the methods of Conkle et al. (1982). We used their System D ("morpholine citrate") for malic dehydrogenase (MDH, EC 188.8.131.52), menadione reductase (MNR, EC 184.108.40.206) and shikimate dehydrogenase (SKDH, EC.1.LL25); their System A (Lithium Borate) for alanine amino peptidase (AAP, EC 220.127.116.11), alanine amino transferase (AAT, EC 18.104.22.168), acid phosphatase (ACP, EC 22.214.171.124) and phosphoglucoisomerase (PGI, EC 126.96.36.199); and histidine, pH 8.0, for leucine amino peptidase (LAP, EC 188.8.131.52), peroxidase (PER, EC 184.108.40.206) and fluorescent esterase (Fl-EST, EC 220.127.116.11). These chemical products were purchased from SIGMA.
The zymograms were given a Mendelian interpretation, and the loci were called "variable isozymes locf" instead of "putative loci." Adequate results were not obtained for 6 of the trees or three of the staining systems; these were eliminated from consideration, and the analyses were based upon 104 individuals and seven systems that provided nine variable loci.
Autocorrelation was performed using the methods and formulas of Sokal and Oden (1978). Nominal data were tested using join-counts based upon Gabriel maps. Because each locus proved to have one allele with a frequency higher than 0.6, we compared homozygotes for the common allele to other genotypes. We also tested the total number of unlike joins. We quantified the number of the most common alleles as 1, 0.5 or 0 depending upon whether 2, 1 or 0 copies of this allele were present in an individual respectively (Dewey and Heywood, 1988). For construction of the correlograms, the distances between all pairs of individuals were calculated, and these were divided into 20 size classes. Moran's I was calculated for each class.
Moran's I is calculated using standardized variables, Zi = Xi -.X. If p is the frequency of the most common allele (A1) at a locus, then under our quantification scheme, it is also the mean of X, and the values of Z1 will be 1 - p for A1A1,½ - p for A1A1, and -p for A1A1, where j represents any of the other alleles. In the present case, where P>0.6 for each locus, the homozygotes for the most common allele will have positive values, while the other genotypes will have negative values. We summed the values of Zi over the loci for each individual, producing an index that will have positive values for the genotypes that are mostly homozygous for the most common alleles and negative values for the less common genotypes. We then calculated Moran's I for all the pairs of individuals with these values.
We calculated F statistics according to the methods of Weir and Cockerham (1984) and Hamilton's kinship coefficients according to equation 6 from the work of Queller and Goodnight (1989). The variances of the estimators were estimated by jack-knifing over the loci and over the groups. We prefer to call the kinship coefficients "Hamilton's relatedness," to avoid possible confusion with Malecot's coefficient of kinship and Wright's coefficient of relationship. The critical values for relatedness were estimated by randomly assigning the individuals to groups, and calculating r values for 10,000 replicates. Genetic diversity and spatial structure were determined using GeneAlex 6.0 (Peakall and Smouse, 2006) and GS+TM5 (Gamma Design Software, Plainwell, MI, USA) (Chen et al., 2008) and SPAGeDi 1.3 (Hardy and Vekemans, 2002).
The nine resolved loci all presented three or four alleles (Table 1). In each case, one allele had a frequency of at least 0.6. The values of F were positive and significant (P ≤ 0.05) for all the loci except PER4, for which we found no homozygotes for the uncommon alleles. The average of F over the loci was 0.197. Fast-running allele 1 of PER2 is found in other populations of Raulí (Table 1).
Nominal data. The physical distribution of the 104 analyzed individuals is shown in Figure 1. The numbers assigned to the individuals in the field run roughly from the bottom to the top of Figure 1. Because each locus had one high-frequency allele, the homozygotes for the most common allele were designated as A and compared with the other genotypes (B) using join-counts of the Gabriel map. A few individuals could not be scored for all loci, so we constructed the corresponding Gabriel map for each locus according to the individuals scored.
|Figure 1. Representation of the locations of 104 individuals of Nothofagus nervosa in a population located 4 km west of Recinto, Bío-Bío Region, Chile. The scales are in meters (m). Dotted and solid circles are groups of neighboring trees determined by successively eliminating the largest distances in the Gabriel map. The shaded areas indicate trees with significant relatedness (Hamilton's r) values.|
Five of the nine analyzed loci showed an excess of like join-counts, but only PGI had a standard normal greater than 2.0 (Table 2). The total number of unlike joins was fewer than expected for five loci but was significant only for PGI (Table 2). Thus, join-counts did not provide clear evidence for the spatial genetic structure, contrary to our original expectations.
Interval data. The value of Moran's I, calculated using pairs of individuals connected in the Gabriel maps, was positive for the most common allele in five of the nine loci, but only two were significant, PGI and ACP (data not shown). The correlograms using all the pairwise distances divided into 20 groups showed significant (positive) values of I for the smallest distance class (Figure 2) only for MDH, PER2 and PGI, and for the second smallest class for PER3 and SKDH. This indicates that there is some tendency for neighboring trees to be more similar, but this shows only a weak genetic structure.
|Figure 2. Correlograms for two enzymes and a summary of individuals in a population of Nothofagus nervosa located 4 km west of Recinto, Bío-Bío Region, Chile. The physical distances between all pairs of individuals were divided into 20 sizes classes, and Moran's I was calculated for each class. PGI, phosphoglucoisomerase; MDH, malic dehydrogenase.|
Moran's I for the Z sums was significant and positive for the first two distance classes (Figure 2). This overall measure confirms that there is some genetic structure in the population and indicates a tendency for the more common genotypes to be neighbors, and the same is true for the less common genotypes.
The range of values in the smallest distance class was 0.8-5.1 m and for the second class approximately 5-12 m (the ranges varied somewhat for the different loci because of the different numbers of individuals). Using Sokal's criterion of the distance at which a correlogram first intercepts the abscissa (Figure 2), we should conclude that for N. nervosa, more similar individuals tend to be found within approximately 10 m of one another.
Relatedness of groups within the population
To use F statistics or relatedness, one must first define the groups to be compared. Previous studies have used agglomerative methods, defining quadrants that are then combined. A divisive method to delimit the subgroups within the population was used by successively eliminating the largest distances in the Gabriel map of the 104 individuals. This procedure first separated the last individual sampled (at the top of Figure 1), which was not further considered. Next, the largest distances were eliminated, producing 2, 3 and up to 10 groups. The r values measure the relatedness of the individuals in a group relative to the whole population. Note that the r value of a group does not change when other groups are subdivided. We also calculated F statistics for 2, 4, 6 and 10 groups.
The first division separates a group of 24 more-isolated individuals, including the two circles at the upper left of Figure 1, which are relatively unrelated compared with the larger group of 79 (Figure 3). Further division of the smaller group also revealed unrelated individuals. Subdivision of the larger group (Figure 3) first produced two groups whose relatedness values, although positive, were not significant. Further subdivisions revealed three physical groups of significantly related individuals interspersed with groups of unrelated individuals (Figures 1 and 3).
|Figure 3. Schematic representation of neighboring groups of individuals by successive elimination of the largest distances in the Gabriel map in a Nothofagus nervosa population. The numbers inside the circles indicate the number of individuals in a group. The numbers to the right of the circles are the values of Hamilton's (relatedness) for each group, relative to the other groups on the same horizontal line, for 2, 4, 6, and 9 groups. The asterisks indicate relatedness values significant at the 0.05 level using a randomization test.|
To determine the significance of the observed relatedness values, we divided the individuals randomly into groups of the observed sizes and calculated r for each group. Table 3 presents the results of 10,000 replicates of this procedure, for six groups and 10 groups. Note that standard deviations are appropriated because the observed values are not means. Because there was an allele at high frequency for each locus, a number of genotypes presented few differences, thus critical values are large for the small groups. The 5% level for a group of 6 (0.342, Table 3) is greater than the expected value for a group of half sibs; for this reason, we did not continue to produce larger numbers of smaller groups.
The main increase in FST was produced between the values for two groups and for four groups (Table 4). Specifically, the group of 36 trees (the dotted circle at the bottom of Figure 1) has allele frequencies somewhat different than the rest of the population. The Fst values were significant for four or more groups (Table 4). We used jack-knife values over the loci because jack-knifing over the groups produced a larger variance. The FIS values for the groups slowly decreased as the number of groups increased; however, the Fis for the 10 groups is only 9% less than for the entire population (Table 4). The F statistics would allow us to conclude that there is some genetic structure in the population but only discovers the group of 36 as different. Note that this group includes two groups of related individuals and two groups of unrelated individuals (Figures 1 and 3).
What is known of the reproductive biology of Nothofagus nervosa led us to expect that within the population, genetic structure should be detected. Although Raulí is monoecious, anemochorous and anemophilous, most seeds fall near the mother tree, and pollen production is rather limited (Riveros et al., 1995). While Raulí has not been tested for self-compatibility, other species of the genus are self-incompatible or almost completely so (Riveros et al., 1995), as is the case for Fagus sylvatica (Merzeau et al, 1994). The average F for 20 populations of Raulí was 0.186 (Carrasco, 1998; Carrasco and Eaton, 2002), whereas the F for this population was 0.197, which represents a considerable level of inbreeding and suggests that there might be some selfing.
Under these conditions, it is rather surprising that spatial autocorrelation methods did not detect more structure. Autocorrelation statistics revealed a weak spatial genetic structure in the population, which was better detected by Moran's I than by join-counts. However, only three loci had significant I values for the two smallest distance classes, and two were significant using the Gabriel map.
Some simulation studies of isolation by distance have "sampled" every individual in the virtual populations and concur that both join-counts (Epperson and Li, 1997; Geng et al, 2009) and Moran's I for the smallest distance classes (Epperson and Li, 1996; Chung et al, 2000) should accurately detect spatial genetic structure as long as the neighborhood size is not too large. In marked contrast, a number of field studies using species of Fagaceae have found only a weak structure using spatial autocorrelation statistics (Sork et al, 1993; Houston and Houston. 1994; Merzeau et al., 1994; Leonardi and Manozzi. 1996; Tero et al, 2005). It is worthwhile to consider the possible causes for this discrepancy.
One possibility is that the assumptions of the model of isolation by distance are not met in the natural populations, in particular where pollen and/ or seed dispersion is not very limited (Tero et al., 2005; Epperson, 2005). There is some evidence that this may be the case in Silene tatarica, where seed dispersal seems to be the limiting factor (Tero et al, 2005), and in Quercus for both seeds (Houston and Houston 1994) and pollen (Dow and Ashley, 1996). For C. tepejilote, nonrandom genetic distribution among the nearest neighbors was detected even from small spatial values, and this distribution seems to be consistent with the neighborhood size (of approximately 300 individuals). For this species, seed dispersal, mortality among life cycle stages, overlapping generations and contrasting traits of mating and reproduction could have influenced the standing spatial genetic structure within populations (Luns et al., 2007).
The model of isolation by distance also requires that enough time has elapsed so that a "quasi-equilibrium" state has been reached (Epperson, 1993; Epperson, 2005). In simulation studies, approximately 50 generations after a random start are needed to reach this stage, at which point, estimations such as Moran's I have stabilized (Epperson and Li, 1997; Epperson, 2005; Chen et al, 2008). Although this assumption may be reasonable for many species, it may not be valid for long-lived temperate forest trees whose ranges have changed considerably due to the last glacial cycle (Yacine and Lumaret, 1989; Leonardi and Menozzi, 1995; Villagrán, 1991; Aldrich et al, 2005). The current distribution range for N. nervosa on the Chilean side is considered a plant refugia in the last glacial period (Villagran, 1991; Carrasco and Eaton, 2002; Pastorino and Gallo, 2002; Premoli et al, 2002). Long-living tree species in the southern temperate forests in Chile show a genetic fragmentation that could be a consequence of their geographic fragmentation during that ice period (Torres-Diaz et al, 2007) and the re-colonization process after the glacial periods (Marchelli and Gallo, 2004). For N. nervosa, it was suggested that genetic diversity was different between the eastern populations located in Argentina and the western populations located on the Chilean side (Marchelli and Gallo, 2004), and these differences could be related to the existence of glacial refugia. However, it could also reflect low levels of gene flow given the eastward unidirectional winds that limit pollen dispersal or by the low seed dispersal alongside the eastern populations.
The simulations of isolation by distance have all used a "saturated" environment, a 100 x 100 matrix always filled with "individuals." This simulates species that form monospecific stands, such as Fagus sylvatica, but may not be appropriate to judge species such as Nothofagus nervosa that occur in multi-specific assemblages and whose density is variable within a stand.
Finally, both the model of isolation by distance and the simulations assume that dispersion of gametes and offspring is homogeneous over the whole population, which eventually produces a given level of structure. This kind of spatial genetic structure, which consists of patches of homozygotes with relatively few heterozy-gotes in between, is well measured by spatial autocorrelation methods. These "patches" or neighborhoods were smaller than the 25 and 28 m estimated for Quercus laevis Walter (Berg and Hamrick, 1994) and Glycine soja (Zhao et al., 2009) and the 40 m estimated for Dipteryx alata (Suarez et al., 2008). However, they were similar in size to the 10 m estimated for Quercus petraea (Matt.) Liebl. (Bacilieri et al., 1994) and mangrove (Geng et al., 2009) and larger than the 5 m estimated for Fumana thymifolia (Jump et al., 2008). However, a heterogeneous genetic structure in a population would not give such clear results because neighbors would be relatives only in parts of the population.
These findings show that the studied population of N. nervosa is a mosaic composed of groups of related individuals mixed with groups of unrelated individuals. A heterogeneous structure is not expected under the model of isolation by distance and must be accounted by other models. Seed dispersion and/ or seedling establishment must have been irregular in this population, very limited at some times and much less limited in others because it is hard to imagine a model of selection at the microhabitat level that would produce adjacent groups of related and unrelated individuals. As in many temperate trees, N. nervosa tends to reproduce in clearings or gaps (Veblin and Donoso, 1987; Carrasco and Eatton 2002; Pollman, 2003, 2005), which together with the limited seed dispersion may produce "family patches" (Houston and Houston, 1994). Seed dispersion may not be exponential, thus colonization may occur injumps rather than an advancing wave (Nichols and Hewitt, 1994). Marchelli and Gallo (2004) found that the genetic diversity of N. nervosa did not relate to the geographic distance. This supports the idea that the heterogeneous composition of individuals should be related to seed dispersion and/or recruitment inside the stands.
We suggest that an equilibrium theory of population structure could not be an appropriate model for N. nervosa nor for many species whose populations do not persist for many generations or have not occupied their present locations for a long time. We propose that many species may have heterogeneous genetic structure within their populations, and this may explain why spatial autocorrelation methods have not detected more structure. More simulations of non-equilibrium populations are needed to predict the conditions under which structure will be produced, what forms it will take and how long it will last.
Mitton (1992) predicted that the coefficient of relationship would "allow direct examination of the degree relatedness within stands of forest trees." Relatedness makes use of data for all the variable loci, which are summarized in a single estimator r. Because kinship structuring affects all the neutral loci equally and independently (Malecót, 1948; Lewontin and Krakauer, 1973), relatedness should be an adequate estimator to detect structure within populations. However, in more or less continuous populations, we need criteria to determine the groups of individuals that will be compared. Successively removing the largest distances of the Gabriel map produced adequate results in the present case, but other exploratory methods could detect related groups. In those cases in which enough family groups may be determined within a population, it should be possible to test whether a locus is subject to selection because its pattern of variation should differ only from subject to kinship structuring.
We thank Dionisia Sepúlveda for her electrophoresis skill and Oscar Chandía for able assistance in the field and laboratory. This work was financed by European Community Project N° CI-CT 930043 and Fondo de Investigación del Bosque Nativo de CONAF, Project 025/2010.
Aldrich, P.R., J.C Glaubitz, G.R Parker, O.E Rhodes, and C.H Michler. 2005. Genetic structure inside a declining red oak community in old-growth forest. Journal of Heredity 96: 627-634. [ Links ]
Bacilieri, R., T Labbe, and A. Kremer. 1994. Intra-specific genetic structure in a mixed population of Quercus petraea (Matt.) Leibl. and Q. robur L. Journal of Heredity 73: 130-141. [ Links ]
Bawa, K.S., D.R. Perry, and J.H. Beach. 1985. Reproductive biology of tropical lowland rain forest trees. I. Sexual systems and incompatibility mechanisms. Am Journal of Botany 72: 331-345. [ Links ]
Belletti, P., S. Leonardi, I. Monteleone, and R. Piovani. 2005. Allozyme variation in different species of deciduous oaks from northwestern Italy. Silvae Genetica 54: 9-16. [ Links ]
Berg, E., and J. Hamrick. 1994. Spatial and genetic structure of two sandhills oaks: Quercus laevis and Quercus margaretta (Fagaceae). Am Journal of Botany 81: 1-7. [ Links ]
Bruno, C., R. Macchiavelli, and M. Balzarini. 2008. Non-parametric smoothing of multivariate genetic distances in the analysis of spatial population structure at fine scale. Theoretical and Applied Genetics 117: 435-447. [ Links ]
Carrasco, B. 1998. Patrones de variabilidad alo-enzimática en raulí (Nothofagus alpinanervosa (Poepp. et Endl.) Oerst.): Efecto de su sistema reproductivo, distribución geográfica y última glaciación. MSc. Thesis, University of Chile. [ Links ]
Carrasco, B., and L. Eaton. 2002. Natural history and genetic structure of raulí (Nothofagus nervosa (Phil.) Dim. et Mil). Forest Genetics 9:275-284 [ Links ]
Chen, X.W., B.L. Li, and X.S. Zhang. 2008. Using spatial analysis to monitor tree diversity at a large scale: a case study in Northeast China Transect. Journal of Plant Ecology-UK 1: 137-141. [ Links ]
Chung, M.G., M.Y. Chung, G.S. Oh, and B. K. Epperson. 2000. Spatial genetic structure in a Neolitsea sericea population (Lauraceae). Journal of Heredity 85: 490-497. [ Links ]
Conkle, M.T., P. Hodgkiss, L. Nunnaly, and S. Hunter. 1982. Starch gel electrophoresis of conifer seeds: A laboratory manual. USDA, U.S. Forest Service. [ Links ]
Dewey, S.E., and J. S. Heywood. 1988. Spatial genetic structure in a population of Psycotria nervosa. I. Distribution of genotypes. Evolution 42: 834-838. [ Links ]
Donoso, C. 1995. Bosques templados de Chile y Argentina: Variación, Estructura y Dinámica. 3rd Edition. Editorial Universitaria, Santiago, Chile. [ Links ]
Dow, B.D., and M.V. Ashley. 1996. Microsatellite analysis of seed dispersal and parentage of saplings in bur oak, Quercus macrocarpa. Molecular Ecology 5: 615-627. [ Links ]
Epperson, B.K. 1993. Recent advances in correlation studies of spatial patterns of genetic variation. Evolutionary Biology 27: 95-155. [ Links ]
Epperson, B.K. 2000. Spatial and spacetime correlations in ecological models. Ecological Modelling 132: 63-76. [ Links ]
Epperson, B.K. 2003. Covariances among join-count spatial autocorrelation measures. Theoretical Population Biology 64: 81-87. [ Links ]
Epperson, B.K. 2005. Estimating dispersal from short distance spatial autocorrelation. Journal of Heredity 95: 7-15. [ Links ]
Epperson, B.K., and T.Q. Li. 1996. Measurement of genetic structure within populations using Moran's spatial autocorrelation statistics. P.N.A.S. (USA) 93: 10528-10532. [ Links ]
Epperson, B.K., and T.Q. Li. 1997. Gene dispersal and spatial genetic structure. Evolution 51: 672-681. [ Links ]
Furnier, G.R., P. Knowles, M.A. Clyde, and B.P. Dancik. 1987. Effects of avian seed dispersal on the genetic structure of whitebark pine populations. Evolution 41: 607-612. [ Links ]
Gapare, W. J., and S.N. Aitken. 2005. Strong spatial genetic structure in peripheral but not core populations of Sitka spruce <Picea sitchensis (Bong.) Carr.>. Molecular Ecology 14: 2659-2667. [ Links ]
Geng, Y.P., S.Q. Tang, T.R. Tashi, Z.P. Song, G.R. Zhang, L.Y. Zeng, J.Y. Zhao, L. Wang, J. Shi, J.K. Chen, and Y. Zhong. 2009. Fine- and landscape-scale spatial genetic structure of cushion rockjasmine, Androsace tapete (Primulaceae), across southern Qinghai-Tibetan Plateau. Genetica 135: 419-427. [ Links ]
Gregorius, H., J. Krauhausen, and G. Müller-Starck. 1986. Spatial and temporal genetic differentiation among the seed in a stand of Fagus sylvatica L. Journal of Heredity 57: 255-262. [ Links ]
Hardy, O.J., and Vekemans, X. 2002. SPAGEDi: a versatile computer program to analyze spatial genetic structure at the individual or population levels. Molecular Ecology Notes 2: 618-620. [ Links ]
Hendry, A.P., and T. Day. 2005. Population structure attributable to reproductive time: isolation by time and adaptation by time. Molecular Ecology 14: 901-916. [ Links ]
Heywood, J.S. 1991. Spatial analysis of genetic variation in plant populations. Annual Review Ecology Evolution, and Systematics 22: 335-355. [ Links ]
Houston, D.B., and D. R. Houston. 1994. Variation in American beech (Fagus grandifolia Ehrh.): Isozyme analysis of genetic structure in selected stands. Silvae Genetica 43: 277-284. [ Links ]
Jump, A.S., L. Rico, F. Lloret, and J. Peñuelas. 2008. Microspatial population genetic structure of the Mediterranean shrub Fumana thymifolia. Plant Biology 11: 152-160. [ Links ]
Kuss, P., A.R. Pluess, H.H. Aegisdottir, and J. Stocklin. 2008. Spatial isolation and genetic differentiation in naturally fragmented plant populations of the Swiss Alpes. Journal of Plant Ecology-UK 1: 149-159. [ Links ]
Leonardi, S., and P. Menozzi. 1995. Genetic variability of Fagus sylvatica L. in Italy: The role of postglacial recolonization. Journal of Heredity 75: 35-44. [ Links ]
Leonardi, S., and P. Menozzi. 1996. Spatial structure of genetic variability in natural stands of Fagus sylvatica L (beech) in Italy. Journal of Heredity 77: 359-368. [ Links ]
Levin, D., and H. Kerster. 1974. Gene flow in seed plants. Evolutionary Biology 17: 139-220. [ Links ]
Lewontin, R.C., and J. Krakauer. 1973. Distribution of gene frequency as a test of the theory of the selective neutrality of polymorphisms. Genetics 74: 175-195. [ Links ]
Lopez-Aljorna, A., M.A. Bueno, I. Aguinagalde, and J.P. Martin. 2007. Fingerprinting and genetic variability in cork oak (Quercus suber L.) elite trees using ISSR and SSR markers. Annals of Forest Science 64: 773-779. [ Links ]
Luns, R., B.K. Epperson, and K. Oyama. 2007. Spatial genetic structure of two sympatric neotropical palms with contrasting life histories. Journal of Heredity 95: 298-305 [ Links ]
Malecót, G. 1948. Les Mathématiques de l'Hérédité. Masson, Paris. [ Links ]
Marchelli, P., and L.A. Gallo. 2000. Genetic analysis of isozyme variants in open pollinated families of southern beech Nothofagus nervosa (PHIL.) DIM. et MIL. Silvae Genetica 49: 90-98. [ Links ]
Marchelli, P., and L.A. Gallo. 2004. The combined role of glaciation and hybridization in shaping the distribution of genetic variation in a Patagonian southern beech. Journal of Biogeography 31: 451-460. [ Links ]
Merzeau, D., B. Comps, B. Thiébaut, J. Cuguen, and J. Letouzey. 1994. Genetic structure of natural stands of Fagus sylvatica L. (beech). Journal of Heredity 72: 269-277. [ Links ]
Mitton, J.B. 1992. The dynamic mating system of conifers. New Forests 6: 197-216. [ Links ]
Muona, O. 1990. Population genetics in forest tree improvement. In A.H.D. Brown, M.T. Clegg, A.L. Kahler, and B.S. Weir (eds.). Sinauer, Sunderland, Massachusetts. Plant Population Genetics, Breeding and Genetic Resources. [ Links ]
Nichols, R.A., and G.M. Hewitt 1994. The genectic consequences of long-distance dispersal during colonization. Journal of Heredity 72:312-317. [ Links ]
Pastorino, M.J., and L.A. Gallo. 2002. Quaternary evolutionary history of Austrocedrus chilensis, a cypress native to the Andean-Patagonian forest. Journal of Biogeography 29: 1167-1178. [ Links ]
Peakall, R., and P.E. Smouse. 2006. GENALEX 6: genetic analysis in Excel. Population genetic software for teaching and research. Molecular Ecology Notes 6: 288-295. [ Links ]
Pollmann, W. 2003. Stand structure and dendroecology of an old-growth Nothofagus forest in Conguillio National Park, south Chile. Forest Ecology and Management 176: 87-103. [ Links ]
Pollmann, W. 2005. A long-term record of Nothofagus dominance in the southern Andes, Chile. Austral Ecology 30: 91-102. [ Links ]
Premoli, A.C., C.P. Souto, A.E. Rovere, T.R. Allnut, and A.C. Newton. 2002. Patterns of isozyme variation as indicators of biogeographic history in Pilgerodendron uviferum (D. Don) Florin. Diversity and Distributions 8: 57-66. [ Links ]
Queller, D.C., and K.F. Goodknight. 1989. Estimation of relatedness using genetic markers. Evolution 43: 258-275. [ Links ]
Riveros, M., M. Paredes, M. Rosas, E. Cárdenas, J. Armesto, M. Arroyo, and B. Palma. 1995. Reproductive biology in species of the genus Nothofagus. Environmental and Experimental Botany 35: 519-524. [ Links ]
Sampson, J., S. Hopper, and S. James. 1990. Temporal variation in allele frequencies in the pollen pool of Eucalyptus rhodantha. Journal of Heredity 65: 189-199. [ Links ]
Sampson, J.F., and M. Byrne. 2008. Outcrossing between an agroforestry plantation and remnant native populations of Eucalyptus loxophleba. Molecular Ecology 17: 2769-2781. [ Links ]
Seltmann, P., A. Cocucci, D. Renison, A. Cierjacks, and I. Hensen. 2009. Mating system, outcrossing distance effects and pollen availability in the wind-pollinated treeline species Polylepis australis BITT (Rosaceae). Basic and Applied Ecology 10: 52-60. [ Links ]
Sethi, P., and H.F. Howe. 2009. Recruitment of Hornbill-Dispersed Trees in Hunted and Logged Forests of the Indian Eastern Himalaya. Conservation Biology 23: 710-718. [ Links ]
Slatkin, M., and H.E. Arter. 1991. Spatial autocorrelation methods in population genetics. The American Naturalist 138: 499-517. [ Links ]
Sokal, R.R., and G.M. Jacquez. 1978. Testing inferences about microevolutionary processes by means of spatial autocorrelation analysis. Evolution 45: 152-168. [ Links ]
Sokal, R.R., and M.L. Oden. 1991. Spatial autocorrelation analysis as an inferential tool in population genetics. The American Naturalist 138: 518-521. [ Links ]
Sokal,R.R., and N. L. Oden. 1978. Spatial autocorrelation in biology. 1. Methodology. Biol. J. Linn. Soc. 10 199-228. [ Links ]
Sork, V., S. Huang, and E. Wiener. 1993. Macrogeographic and fine scale genetic structure in North American oak species, Quercus rubra L. Annals of Forest Science 50: 261-270. [ Links ]
Suarez, A., G. Castillo, and M.I. Chacon. 2008. Genetic diversity and spatial genetic structure within a population of an aromatic shrub, Lippia origanoides (Verbenaceae), in the Chicamocha Canyon, northeastern Colombia. Genetics Research 90: 455-465. [ Links ]
Tero, N., J. Aspi, P. Siikamaki, and A. Jakalaniemi. 2005. Local genetic population structure in an endangered plant species, Silene tatarica (Caryophyllaceae). Journal of Heredity 94: 478-487. [ Links ]
Torres-Diaz, C., E. Ruiz, F. González, G. Fuentes, and L.A. Cavieres. 2007. Genetic diversity in Nothofagus alessandrii (Fagaceae), an endangered endemic tree species of the coastal maulino forest of central Chile. Annals of Botany 100: 75-82. [ Links ]
Veblin, T.T., and C. Donoso. 1987. Alteración y dinámica regenerativa de las especies chilenas de Nothofagus de la Región de los Lagos. Bosque 8: 133-142. [ Links ]
Villagrán, C. 1991. Historia de los bosques templados del sur de Chile durante el tardiglacial y postglacial. Revista Chilena de Historia Natural 64: 447-460. [ Links ]
Weir, B., and C. Cockerham. 1984. Estimating F-statistics for the analysis of population structure. Evolution 38: 1358-1370. [ Links ]
Wright, S. 1943. Isolation by distance. Genetics 28: 114-138. [ Links ]
Yacine, A., and R. Lumaret. 1989. Genetic diversity in holm-oak (Quercus ilex L.): Insight from several enzyme markers. Silvae Genetica 38: 140-148. [ Links ]
Zhao, R., H. Xia, and B.R. Lu. 2009. Fine-scale genetic structure enhances biparental inbreeding by promoting mating events between more related individuals in wild soybean (Glycine soja; Fabaceae) populations. American Journal of Botany 96: 1138-1147. [ Links ]
Received June 29, 2010. Accepted January 15, 2011.
Corresponding author: email@example.com