Variation in morphological traits in a recent hybrid zone between closely related Quercus liaotungensis and Q. mongolica (Fagaceae)

Aims Hybridization usually leads to gene introgression between related species in hybrid zones, associated with complex patterns of morphological variation. Nevertheless, previous studies have tended to ignore the effects of geographic variation in hybridization rates on species taxonomy. This study aims to investigate the variation of morphological traits between two sympatric and taxonomically confused oak species, Quercus liaotungensis and Q. mongolica, and reveal the effects of hybridization rates on morphological traits and the taxonomic boundary.Methods We used seven microsatellite loci to evaluate species status and measured 15 morphological traits in 26 trees in the recent hybrid zone between Q. liaotungensis and Q. mongolica, and we characterized the differences between the two oak species and their hybrids for the investigated traits.Important findings Molecular analyses indicated that 74% of 78 sampled maternal trees were hybrids between Q. liaotungensis and Q. mongolica although the observed morphological variation suggested that they had remained distinct species. Across all of the differentiated leaf and reproductive traits, the hybrids expressed patterns similar to Q. liaotungensis, which may suggest dominant expression of parental characters. These results are consistent with our expectation that hybrids will be difficult to distinguish from parental species in a recent hybrid zone.     

Natural hybridization and hybrid zones between Quercus crassifolia and Quercus crassipes (Fagaceae) in Mexico: morphological and molecular evidence

Hybrid zones provide interesting systems to study genetic and ecological interaction between different species. The correct identification of hybrids is necessary to understand the evolutionary process involved in hybridization. An oak species complex occurring in Mexico formed by two parental species, Quercus crassifolia H. & B. and Q. crassipes H. & B., and their putative hybrid species, Q. dysophylla, was analyzed with molecular markers (random amplified polymorphic DNA [RAPDs]) and morphological tools in seven hybrid zones (10 trees per taxa in each hybrid zone) and two pure sites for each parental species (20 trees per site). We tested whether geographic proximity of hybrid plants to the allopatric site of a parental species increases its morphological and genetic similarity with its parent. Seventeen morphological traits were measured in 8700 leaves from 290 trees. Total DNA of 250 individuals was analyzed with six diagnostic RAPD primers. Quercus crassifolia differed significantly from Q. crassipes in all the examined characters. Molecular markers and morphological characters were highly coincident and support the hypothesis of hybridization in this complex, although both species remain distinct in mixed stands. Clusters and a hybrid index (for molecular and morphological data) showed that individuals from the same parental species were more similar among themselves than to putative hybrids, indicating occasional hybridization with segregation in hybrid types or backcrossing to parents. Evidence does not indicate a unidirectional pattern of gene flow.     

Evidence for shared ancestral polymorphism rather than recurrent gene flow at microsatellite loci differentiating two hybridizing oaks (Quercus spp.)

Quercus petraea and Quercus robur are two closely related oak species, considered to hybridize. Genetic markers, however, indicate that despite sharing most alleles, the two species remain separate genetic units. Analysis of 20 microsatellite loci in multiple populations from both species suggested a genome-wide differentiation. Thus, the allele sharing between both species could be explained either by low rates of gene flow or shared ancestral variation. We performed further analyses of population differentiation in a biogeographical setting and an admixture analysis in mixed oak stands to distinguish between both hypotheses. Based on our results we propose that the low genetic differentiation among these species results from shared ancestry rather than high rates of gene flow.     

Contrasting ectomycorrhizal fungal communities on the roots of co-occurring oaks (Quercus spp.) in a California woodland

Plant host species is considered an important factor influencing ectomycorrhizal (EM) communities. To gain insights into the role of host species in structuring EM communities, EM communities on sympatric oak (Quercus) species were compared in the Sierra Nevada foothills of California. Using molecular methods (polymerase chain reaction, cloning, restriction fragment length polymorphism and DNA sequencing), EM fungi on roots of deciduous Quercus douglasii and evergreen Quercus wislizeni trees were identified from 64 soil cores. The total EM species richness was 140, of which 40 taxa were detected on both oak hosts. Greater diversity and frequency of EM fungi with epigeous fruiting habit were found on Q. wislizeni, while taxa in the Ascomycota were more frequent and diverse on Q. douglasii. Using ordination, it was determined that both soil extractable phosphorus and oak host species explained a significant proportion of the variation in EM species distribution. These results indicate that plant host species can be an important factor influencing EM fungal community composition, even within congeneric trees.     

Molecular characterization and genetic structure of Quercus acutissima germplasm in China using microsatellites

Quercus acutissima is native to eastern Asia. It has a wide distribution in China and China is an important component in understanding the ecology and genetic structure of this species. Q. acutissima attained high economic value for hardwood product and can be managed as an energy tree species. To investigate the genetic variation of Q. acutissima provenances, 12 microsatellite primer pairs were used to analyze 672 trees sampled from 28 provenances of Q. acutissima in China. All of the tested microsatellite loci proved to be effective for the studied Q. acutissima provenances. The results revealed that allele numbers varied from 5 to 13 per locus, with an average of 8 alleles per locus. The mean observed heterozygosity and expected heterozygosity were 0.4927 and 0.7023, respectively. The relatedness of the provenances was studied using the arithmetic mean algorithm based on Nei’s genetic distance and principal coordinates analysis. Interestingly, both approaches revealed two main groups: one consisted of the eastern Chinese provenances, and the other comprised of the western Chinese provenances. An analysis of molecular variance indicated that most genetic variation was contained within populations (84 %). The two microsatellite markers developed in this study may be employed for genetic characterization of other oak species. Considering the management or breeding programs of Q. acutissima provenances in China, we should treat each main group as a single gene resource.     

Demography and recruitment limitations of three oak species in California

We review published studies on the demography and recruitment of California oak trees and focus on the widespread dominant species of the foothill woodlands, Quercus douglasii, Q. lobata, and Q. agrifolia, to ascertain the nature and strength of evidence for a decline in populations of these species. The vast majority of studies have been of short duration (less than three years), focused on the acorn and seedling life stages, and conducted at few locations within each species geographic range. We summarize the extensive body of research that has been conducted on the biological and physical factors that limit natural seedling recruitment of oaks. The oak "regeneration problem" has largely been inferred from current stand structure rather than by demographic analyses, which in part reflects the short-term nature of most oak research. When viewed over longer periods of time usingfield surveys or historical photos, the evidence for a regeneration problem in foothill oaks is mixed. Q. douglasii shows very limited seedling or sapling recruitment at present, but longer term studies do not suggest a decline in tree density, presumably because rare recruitment is sufficient to offset low rates of mortality of overstory individuals. Q. agrifolia appears to be stable or increasing in some areas, but decreasing in areas recently impacted by the disease Phytophthora ramorum. Evidence from the few available studies is more consistent in suggesting long-term declines in foothill populations of Q. lobata. Long-term monitoring, age structure analysis, and population models are needed to resolve the current uncertainty over the sustainability of oak woodlands in California.     

Genetic variation and differentiation in population of Japanese emperor oak (Quercus dentata Thunb.) and Mongolian oak (quercus mongolica fisch. ex ledeb.) in the south of the Russian far east

Allozyme variation of Japanese emperor oak Quercus dentata Thunb. and Mongolian oak Quercus mongolica Fisch. ex Lebed.) was examined in 11 populations of these species from southern Primorye. Analysis of 18 loci showed that in these populations, 66.7% of genes of Mongolian oak are polymorphic, the number of alleles per locus being 2.28. In Japanese emperor oak, these parameters were respectively 59.8 and 2.36%. The observed heterozygosity in both species did not differ from the expected values, constituting 0.198 in Japanese emperor oak and 0.161 in Maongolian oak. The Nei's distance between Quercus mongolica and Q. dentate was 0.053. The Japanese emperor oak and Mongolian oak populations formed separate clusters in an UPGMA dendrogram. Taking into account morphological and ecological isolation of the species, we suggest that the differences in gene frequencies may reflect selection acting in favor of adaptation to different growth conditions and counteracting between-species hybridization.     

The clonal structure of Quercus geminata revealed by conserved microsatellite loci

The scrub oak communities of the southeastern USA may have existed at their present locations for thousands of years. These oaks form suckers, and excavations of root systems suggest that clones may occupy very large areas. Resolution of the clonal nature of scrub oaks is important both to manage the tracts of this ecosystem that remain, and in conducting long-term ecological studies, where the study area must substantially exceed the area occupied by any single clone. Microsatellites were used to determine the genetic diversity of a dominant oak species within a 2-ha long-term experimental site on Merritt Island at the Kennedy Space Center. This area contains a long-term study of the effects of elevated CO2 on the ecosystem. Conservation of seven microsatellite loci, previously identified in the sessile oak, Quercus petraea, was tested in two Florida scrub oak species, Q. geminata and Q. myrtifolia. Sequence analysis revealed that all seven microsatellite loci were conserved in Q. geminata and five loci were conserved in Q. myrtifolia. Six microsatellite loci were polymorphic in Q. geminata and these were subsequently used to investigate the clonal structure of the Q. geminata population. Twenty-one unique combinations of microsatellites, or haplotypes, occurred only once, whereas the remaining 26 individuals belonged to a total of seven different haplotypes. Trees with identical haplotypes were in close proximity, supporting the interpretation that they were clones. The results showed that there is significant genetic diversity within the 2-ha experimental site. Microsatellites provided a powerful and noninvasive tool for distinguishing individual genotypes and determining an adequate area for long-term ecosystem studies.     

Genetic admixing of two evergreen oaks,Quercus acuta and Q. sessilifolia (subgenus Cyclobalanopsis), is the result of interspecific introgressive hybridization

In forests worldwide, Quercus is a major genus; however, the boundaries between the constituent species are relatively weak, and hybridization is reported frequently. In this study, we examined Quercus acuta and Quercus sessilifolia (subgenus Cyclobalanopsis), which have a putative hybrid—Q. x takaoyamensis. We investigated leaf morphological traits and microsatellites of Q. acuta and Q. sessilifolia in the area where the two species are both found. Although the leaf traits overlapped, the two species could be distinguished morphologically as demonstrated by principal component analysis based on a range of these traits. They were also genetically differentiated, with F _ST?=?0.104. However, they shared most of the alleles at all eight loci examined, and considerable genetic admixing was detected. Admixture analysis demonstrated that Q. acuta and Q. sessilifolia, respectively, contained 11 and 24 % of individuals with a probability of less than 0.9 of being correctly assigned to their species. Model-based testing showed that this admixing was created by not only shared ancestral polymorphism but also by hybridization. Effective population size and migration rate were estimated using the coalescent approach. We estimated 8.843 and 71.98 effective numbers of migrants per generation to Q. acuta and Q. sessilifolia, respectively. Theoretically, one to ten migrants per generation are required to prevent complete genetic differentiation. Based on the results of this study, it appears that genetic admixing, with sharing of most alleles, is probably common in the two species and is maintained by interspecific introgressive hybridization.     

Natural hybridisation between Quercus petraea (Matt.) Liebl. and Quercus pubescens Willd. within an Italian stand as revealed by microsatellite fingerprinting

Interspecific gene flow is frequently reported in the genus Quercus . However, interfertile oak species often seem to remain distinct, even within areas of sympatry. This study employed molecular markers to verify, at a fine scale, the presence of interspecific gene flow in a natural population of Quercus petraea and Quercus pubescens . Within a delimited area of 6 ha, all adult trees belonging to the studied oak complex and seeds from a subsample of such trees were collected and analysed using molecular microsatellite markers. A low interspecific genetic differentiation and a high level of interspecific genetic admixture suggested past hybridisation. Paternity inference of seeds allowed the estimation of pollination frequencies from the three groups of pollen donors ( Q. petraea, Q. pubescens , intermediate). We also assayed pollen viability and germinability of each species group. We observed natural hybridisation between Q. petraea and Q. pubescens, with a predominant component in the direction Q. petraea → Q. pubescens : Q. pubescens displayed a higher level of heterospecific pollination by Q. petraea (25.8%) and intermediate morphotypes (14.7%), compared to Q. petraea acting as pollen receptor (with less than 5% heterospecific pollinations). Intermediate 'mother trees' were pollinated in similar proportions by Q. petraea (23.1%), Q. pubescens (37.8%) and intermediate morphotypes (39.1%). The asymmetrical introgression observed for the studied generation may be caused, among other factors, by the relative abundance of trees from each species group in the studied area.     

Mating patterns of a subdivided population of the andean oak (Quercus humboldtii Bonpl., Fagaceae)

Mating patterns play a critical role in the maintenance of genetic variation. We analyzed the mating system in a recently fragmented population of the Andean oak (Quercus humboldtii) using four microsatellite loci. Five fragments in northeastern Colombia, South America, were selected consisting of 30.4 trees on average. We sampled about 30 seeds from three target trees in each fragment and genotyped them with four microsatellite loci with a total of 40 alleles across loci. Progenies were analyzed under the mixed mating system model (MLTR program) and the TwoGener pollen pool structure analyses. The number of unshared pollen donors per family (Nu) was estimated using gametotypic counts with the program HAPLOTYPES. Low selfing (3%) is occurring at the population and fragment level. Biparental inbreeding is present (4.9%), but reduced, in the largest fragment. The average pollen neighborhood size (Nep = 5.1 to 6.1) appears comparable to other oak species in sparse landscapes. In contrast, Nu consists of 2 to 4 main donors, although up to 12 are possible, with the mode positively correlated with fragment size. The Andean oak appears to be a resilient species capable of tolerating population subdivision, provided landscapes include large fragments.     

Effect of a red oak species gradient on genetic structure and diversity of Quercus castanea (Fagaceae) in Mexico

Incipient reproductive barriers are a common characteristic of oak species. Disruption of these barriers promotes changes in diversity and genetic structure of the species involved. Quercus castanea is a red oak with wide geographic distribution in Mexico, which presents atypically high morphological variability when it occurs in sympatry with other red oak species, suggesting that hybridization may explain the observed variation. We tested if the genetic structure and diversity levels of Q. castanea are related to the number of red oak species growing in sympatry. In total, 14 microsatellite (SSRs) primers (six nSSRs and eight cpSSRs) were used in 120 Q. castanea individuals (20/site) belonging to six populations, where the number of red oak species associated varied from zero to five. Results showed a positive and significant relationship between the genetic diversity of Q. castanea and the number of red oak species growing in sympatry, regardless of the marker type or the parameter of genetic diversity analyzed. Also, we found a higher genetic differentiation of Q. castanea populations using cpSSRs in comparison with nSSRs. Our results suggest that temperate forests with high red oaks species richness co-dominated by Q. castanea promote the increase in this species genetic diversity. From a conservation perspective, high genetic diversity levels of foundation species such as Q. castanea may have positive cascade effects extending to other species in the community.     

Hybridization and divergence in multi‐species oak (Quercus) communities

Oaks (Quercus: Fagaceae) commonly interbreed yet retain their morphological, genetic and ecological distinctiveness. Post‐zygotic isolation mechanisms, such as ecologically dependent selection on adaptive loci, may therefore limit introgression. To test this hypothesis, we quantified hybridization and genetic divergence across the contact zone of four red oaks (Quercus section Lobatae) in the Great Lakes region of North America using a suite of 259 amplified fragment length polymorphisms and 27 genic and genomic microsatellite markers. First, we identified hybrids using genetic structure analysis and confirmed the reliability of our assignments via simulations. Then, we identified candidate loci for species maintenance with three complementary tests for selection and obtained partial gene sequences linked to an outlier locus and three other loci. We detected evidence of recent hybridization among all species and considerable gene flow between Q. ellipsoidalis and Q. velutina. Overall, c. 20% of Q. velutina had recent ancestry from Q. ellipsoidalis, whereas nearly 30% of Q. ellipsoidalis had a Q. velutina ancestor. Most loci were negligibly to weakly differentiated among species, but two gene‐linked microsatellites deviated significantly from neutral expectations in multiple, complementary outlier tests. Both outlier loci were located in the same 15‐cM bin on an existing Q. robur linkage map, a region under divergent selection in other oak species. Adaptive loci in this highly differentiated genomic region may contribute to ecological divergence among species and limit introgression.     

Coincidence of small-scale spatial discontinuities in leaf morphology and nuclear microsatellite variation of Quercus petraea and Q. robur in a mixed forest

BACKGROUND AND AIMS: The taxon complex comprising Quercus petraea and Q. robur shows distinct morphologies and ecological preferences, but mostly low differentiation in various types of molecular markers at a broad spatial range. Local, spatially explicit analyses may reveal patterns induced by microevolutionary processes operating mainly over short distances. However, no attempts have been made to date to explore the potential of spatial analyses combining morphological and genetic data of these oaks. METHODS: A mixed oak stand was studied to elucidate the small-scale population genetic structure. All adult individuals were classified and putative hybrids were identified using multivariate discrimination analysis of leaf morphological characters. Likewise, all trees were genotyped with five nuclear microsatellites, and a Bayesian assignment method was applied based on maximum likelihood of multilocus genotypes for taxon and putative hybrid classification. KEY RESULTS: Multivariate analyses of leaf morphological data recognized two groups with few individuals as putative hybrids. These groups were significantly differentiated at the five microsatellites, and genetic taxon assignment coincided well with morphological classification. Furthermore, most putative hybrids were assigned to the taxon found in their spatial neighbourhood. When grouping trees into clusters according to their spatial positions, these clusters were clearly dominated by one taxon. Discontinuities in morphological and genetic distance matrices among these clusters showed high congruence. CONCLUSIONS: The spatial-genetic analyses and the available literature led to the assumption that reproductive barriers, assortative mating, limited seed dispersal and microsite-induced selection in favour of the locally adapted taxon at the juvenile stage may reinforce taxon-specific spatial aggregation that fosters species separation. Thus, the results tend to support the hypothesis that Q. petraea and Q. robur are distinct taxa which share a recent common ancestry. Occasional hybrids are rarely found in adults owing to selection during establishment of juveniles.     

Gene flow and fine-scale genetic structure in a wind-pollinated tree species, Quercus lobata (Fagaceaee)

California Valley oak (Quercus lobata), one of the state's most distinctive oak species, has experienced serious demographic attrition since the 19th century, due to human activities. Recent estimates of pollen dispersal suggest a small reproductive neighborhood. Whether small neighborhood size is a recent phenomenon, a consequence of reduced gene flow caused by demographic changes, or whether it has been historically restricted, remains unclear. To examine this question, we have characterized the spatial genetic structure of N = 191 Q. lobata individuals, spread over an area of 230 ha, using eight microsatellite loci. The observed autocorrelogram suggests an historical standard deviation of gene flow distance of about 350 m per generation, higher than contemporary pollen dispersal estimates. To determine whether our estimates were affected by strong prevailing winds from the west-northwest, we developed and utilized a novel anisotropic autocorrelation analysis. We detected no more than a hint of anisotropy, and we concluded that adult spatial structure is indicative of strong historical signature of "isolation by distance." This historical estimate provides a useful reference value against which to gauge the future gene flow consequences of ongoing anthropogenic disturbance.     

Genetic variation and differentiation in populations of Japanese emperor oak <Emphasis Type="Italic">Quercus dentata</Emphasis> Thunb. and Mongolian oak <Emphasis Type="Italic">Quercus mongolica</Emphasis> Fisch. ex Ledeb.) in the south of the Russian Far East

Allozyme variation of Japanese emperor oak Quercus dentata Thunb. and Mongolian oak Quercus mongolica Fisch. ex Ledeb.) was examined in 11 populations of these species from southern Primorye. Analysis of 18 loci showed that in these populations, 66.7% of genes of Mongolian oak are polymorphic, the number of alleles per locus being 2.28. In Japanese emperor oak, these parameters were respectively 59.8 and 2.36%. The observed heterozygosity in both species did not differ from the expected values, constituting 0.198 in Japanese emperor oak and 0.161 in Mongolian oak. The Nei’s distance between Quercus mongolica and Q. dentata was 0.053. The main contribution to the species differentiation is made by loci Lap-1, Idh-1, Pgm-1, and Ndh-1, although none of them is diagnostic. The Japanese emperor oak and Mongolian oak populations formed separate clusters in an UPGMA dendrogram. Taking into account morphological and ecological isolation of the species, we suggest that the differences in gene frequencies may reflect selection acting in favor of adaptation to different growth conditions and counteracting between-species hybridization.     

Molecular, morphological, and ecological niche differentiation of sympatric sister oak species, Quercus virginiana and Q. geminata (Fagaceae)

The genus Quercus (the oaks) is notorious for interspecific hybrization, generating questions about the mechanisms that permit coexistence of closely related species. Two sister oak species, Quercus virginiana and Q. geminata, occur in sympatry in Florida and throughout the southeastern United States. In 11 sites from northern and southeastern regions of Florida, we used a leaf-based morphological index to identify individuals to species. Eleven nuclear microsatellite markers significantly differentiated between the species with a high correspondence between molecular and morphological typing of specimens. Nevertheless, Bayesian clustering analysis indicates interspecific gene flow, and six of 109 individuals had mixed ancestry. The identity of several individuals also was mismatched using molecular markers and morphological characters. In a common environment, the two species performed differently in terms of photosynthetic performance and growth, corresponding to their divergent ecological niches with respect to soil moisture and other edaphic properties. Our data support earlier hypotheses that divergence in flowering time causes assortative mating, allowing these ecologically distinct sister species to occur in sympatry. Limited gene flow that permits ecological differentiation helps to explain the overdispersion of oak species in local communities.     

Morphological and RAPD analysis of hybridization between Quercus affinis and Q. laurina (fagaceae), two Mexican red oaks

Quercus affinis and Q. laurina are two closely related Mexican red oaks with partially overlapping distributions. Within the area of overlap, there are localities where morphological intergradation occurs. A previous hypothesis explained this pattern as a result of secondary contact between the two species, followed by hybridization and introgression. This possibility was analyzed here by examining foliar and genetic variation in 16 localities situated along a macrogeographic gradient, which included morphologically representative populations of both species and populations from within the area of overlap. Maximum-likelihood hybrid index scores calculated from nine semi-diagnostic RAPD markers indicated a shift in the genetic composition of populations from one species to the other along the macrogeographic gradient, with genetically intermediate populations situated in the area of overlap. Foliar variation followed a partially congruent pattern, but Q. laurina-like morphology predominated in some of the genetically intermediate populations. There were several instances of correlated frequency changeovers of single RAPD markers and morphological characters along the macrogeographic gradient and a few cases of markedly parallel patterns between markers. The results were interpreted as consistent with a hypothesis of secondary contact between the two oak species that has resulted in some differential introgression among markers.     

The Description of Clones of Quercus robur L. and Q. petraea (Matt.) Liebl. with Microsatellites and AFLP in an Ancient Woodland

Abstract: In densely populated areas autochthonous Quercus robur L. (pedunculate oak) and Q. petraea (Matt.) Liebl. (sessile oak) (Fagaceae) populations have been maintained as ancient devastated woodlands. The continuous cutting, grazing and resprouting of such woodlands has enabled the development of clonal structures. For conservation purposes, an analysis of the actual number, size and spatial distribution of clones is necessary, especially when there is an interest in genetic variation of the population. This study describes for the first time - based on microsatellite and AFLP™ analysis - clones in an autochthonous mixed Q. robur and Q. petraea population that has been coppiced and grazed for several centuries. Based on six microsatellite loci and 69 polymorphic AFLP markers, only 14 unique genotypes were detected in a plot that consisted of 80 trees. Clones were observed for both Q. robur and Q. petraea. The largest clone diameters were observed for Q. robur, with distances up to 5.8 m. The observed clone sizes may indicate the age of the trees.     

栓皮栎天然群体SSR遗传多样性研究

利用微卫星(SSR)标记对我国4个省内的5个栓皮栎（Quercus variabilis Bl.）天然群体的遗传多样性进行了研究。16对SSR标记揭示了栓皮栎丰富的遗传多样性：等位基因数（A）平均8.4375个,有效等位基因数(Ne)平均为5.9512个,平均期望杂合度(He)0.8059,Nei多样性指数(h)为0.8041。栓皮栎自然分布区中心地带的群体具有较高的遗传多样性,而人为对森林的破坏将降低林木群体的遗传多样性。栓皮栎群体的变异主要来源于群体内,群体间分化较小,遗传分化系数仅为0.0455。此外,栓皮栎群体间的遗传距离与地理距离之间存在显著的正相关。这些遗传信息为栓皮栎遗传多样性的保护和利用提供了一定依据。Abstract: Genetic diversity of five Quercus variabilis natural populations in four provinces of China was studied with microsatellite (SSR) markers. A relatively high level of genetic diversity was detected in Q. variabilis species with 16 polymorphic microsatellite loci. Average number of alleles (A) and effective number of alleles (Ne) were 8.4375 and 5.9512 respectively. The mean expected heterozygosity (He) was 0.8059 and Nei diversity index (h) was 0.8041. Higher diversity was found with the populations from the central range of the species in contrast to those from peripheral areas and human activities might decrease the genetic diversity of populations. The majority of genetic variation occurred within populations, which could be concluded from the low coefficient of genetic differentiation (Fst=0.0455). In addition, significant correlation was found between geographical distance and genetic distance. All these results present a basis to the conservation and utilization of genetic diversity of Quercus variabilis.