Unveiling the origin of Quercus serrata subsp. mongolicoides found in Honshu,Japan, by using genetic and morphological analyses

Quercus mongolica is a tree found in temperate deciduous forests in east Asia. In Japan, Q. mongolica var. crispula is commonly found; moreover, an oak whose morphology is similar to that of Q. mongolica var. mongolica of the Asian continent has been found in certain areas of Honshu and Hokkaido. Recently, the oak found in Honshu was described as Q. serrata subsp. mongolicoides (QSM). However, genetic comparison between this oak and Q. mongolica var. mongolica of the Asian continent has not been performed; the origin of QSM is thus unknown. Therefore, we aimed to determine the origin of QSM by conducting nuclear microsatellite (nSSR), chloroplast DNA (cpDNA) and leaf morphology analyses for the three taxa, as well as other congeners. The cpDNA variation overlapped among the three taxa, suggesting low discrimination ability for these taxa. Although morphological congruency was found between QSM and Q. mongolica var. mongolica, results of nSSR analyses showed that QSM contained a genetic admixture of Q. mongolica var. mongolica of the Asian continent and Q. mongolica var. crispula of Japan, bolstering an admixture hypothesis. The nSSR and cpDNA analyses suggested that Q. mongolica var. crispula can be the progenitor of Q. mongolica var. mongolica and harbors the ancestral cpDNA haplotypes. Therefore, we concluded that QSM might have been created by an admixture that likely occurred within Japan between Q. mongolica var. crispula and putative relict Q. mongolica var. mongolica, which might have diverged in or around Japan from Q. mongolica var. crispula during the late Pleistocene.     

Geographic patterns of genetic variation in nuclear and chloroplast genomes of two related oaks (<Emphasis Type="Italic">Quercus aliena</Emphasis> and <Emphasis Type="Italic">Q. serrata</Emphasis>) in Japan: implications for seed and seedling transfer

In this study, we assessed geographic patterns of genetic variations in nuclear and chloroplast genomes of two related native oaks in Japan, Quercus aliena and Q. serrata, in order to facilitate development of genetic guidelines for transfer of planting stocks for each species. A total of 12 populations of Q. aliena and 44 populations of Q. serrata were analyzed in this study. Genotyping of nuclear microsatellites in Q. aliena was done with only nine populations (n = 212) due to limited numbers of individuals in two populations, while all 12 populations (n = 89) were used in sequencing chloroplast DNA (cpDNA). In Q. serrata, 43 populations (n = 1032) were genotyped by nuclear microsatellite markers, while cpDNA of 44 populations (n = 350) was sequenced. As anticipated, geographic patterns detected in the variations of Q. aliena’s nuclear genome and its chloroplast haplotype distribution clearly distinguished northern and southern groups of populations. However, those of Q. serrata were inconsistent. The geographic distribution of its chloroplast haplotypes tends to show the predicted differentiation between northern and southern lineages, but geographic signals in the genetic structure of its nuclear microsatellites are weak. Therefore, treating northern and southern regions of Japan as genetically distinct transferrable zones for planting stocks is highly warranted for Q. aliena. For Q. serrata, the strong NE-SW geographic structure of cpDNA should be considered.     

Geographic distribution and origin of the chloroplast T/C‐type in Quercus mongolica var. crispula in northeastern Japan

In order to investigate the historical expansion of common oaks (Quercus sect. Prinus) in the northeastern part of Japan, the relation between the chloroplast haplotypes (I and II) in Quercus mongolica var. crispula and the chloroplast types (T‐ and C‐type) were examined. Complete linkage between haplotype II and chloroplast C‐type was found. The chloroplasts examined in the oak species collected from Sakhalin and Primorski Krai, Russia, and Harbin, China were all T‐type. This suggests that the T to C mutation had occurred in haplotype II in Japan. Neither of haplotype I nor haplotype II was found outside Japan, suggesting both occurred in Japan independently from an ancestral haplotype VI. Haplotype I, which has been known only in Mount Hayachine within Honshu, is distributed southward to the Kanto district, where refugia might have occurred during glaciations.     

Isolation and characterization of microsatellite loci from Quercus mongolica var. crispula

Microsatellites were isolated and characterized for the Japanese oak species, Quercus mongolica var. crispula, distributed in temperate deciduous forests of Japan. Eleven of the 48 primer pairs designed successfully amplified unambiguous and polymorphic single loci among 67 Q. mongolica var. crispula individuals within a plot in southwestern Japan. The observed and expected heterozygosities of the 11 microsatellite markers ranged from 0.522 to 0.896 and from 0.536 to 0.882, respectively. These polymorphic microsatellite markers are useful for estimating pollen‐mediated gene flow in Q. mongolica var. crispula.     

The genetic structure of <Emphasis Type="Italic">Quercus crispula</Emphasis> in northeastern Japan as revealed by nuclear simple sequence repeat loci

Previous studies have reached different discussions about the genetic variation and genetic structure of Quercus crispula populations in northeastern Japan. This is a common oak species in Eastern Asia. Some studies have suggested that the populations in northeastern Japan were derived from those remaining in the southwest after the last glacial maximum (LGM), whilst other studies have found evidence that populations persisted in northeastern Japan during the LGM. Using seven highly polymorphic nuclear simple sequence repeat loci, we investigated the genetic structure of 16 Q. crispula populations along a latitudinal gradient in northeastern Japan (northern Honshu and Hokkaido), spanning about half of the species’ biogeographic range in the country. Although the level of population differentiation was low (F _ST = 0.021; G^￠_\textST G^{\prime}_{\text{ST}}  = 0.090), two geographically differentiated clusters were detected by STRUCTURE analysis. The first cluster included most of the populations in Hokkaido, and may indicate continued survival throughout past glacial periods. We found a significant decrease in allelic richness with latitude, so the second cluster may represent an expansion of the lineage from Honshu during the post-glacial period. These results should enhance our understanding of historical north–south migrations of this species in northeastern Japan.     

A set of primers for the amplification of chloroplast microsatellites in Quercus

The increase in demand for the certification of oak seed lots, as well as control of the geographical origin of oak wood, has led us to develop powerful genetic markers permitting us to discriminate among provenance regions. With the aim of detecting new chloroplast variants, we have identified 17 potential cpSSRs motifs from available oak sequences and tested their variability among French oak populations. Six loci were polymorphic at the intraspecific level in Quercus petraea and Q. robur. Moreover, conservation of the primer pairs was checked on a set of 21 forest tree species and they were all shown to work well on several Quercus species, and even within Fagacaea.     

Population differentiation of sessile oak at the altitudinal front of migration in the French Pyrenees

To assess the effects of altitude on the level and structure of genetic diversity, a genetic survey was conducted in 12 populations of sessile oak (Quercus petraea) located between 130 and 1660 m in two parallel valleys on the northern side of the Pyrenees Mountains. Genetic diversity was monitored at 16 nuclear microsatellite loci and 5 chloroplast DNA (cpDNA) markers. The cpDNA survey suggested that extant populations in both valleys shared the same source populations from the plain. There was no visible trend of nuclear genetic diversity along altitude, even if indirect estimates of effective population sizes revealed a consistent reduction at higher altitudes. Population differentiation, although low, was mostly present among populations of the same valleys and reached similar levels than differentiation across the range of distribution of sessile oak. Contribution to the overall differentiation in the valleys was mostly due to the genetic divergence of the highest populations and the altitudinal variation of allelic frequencies at a few loci. Bayesian inference of migration between groups of populations showed that gene flow is preferentially unidirectional from lower altitudes in one valley to other groups of populations. Finally, we found evidence of clonal reproduction in high altitude populations. The introgression of Quercus robur and Quercus pubescens was also more frequent at the altitudinal margin suggesting that this mechanism may have contributed to the present migration and adaptation of Q. petraea and may also facilitate its future upslope shift in the context of climate change.     

Influence of late Quaternary climate change on present patterns of genetic variation in valley oak,Quercus lobata Née

Phylogeography and ecological niche models (ENMs) suggest that late Quaternary glacial cycles have played a prominent role in shaping present population genetic structure and diversity, but have not applied quantitative methods to dissect the relative contribution of past and present climate vs. other forces. We integrate multilocus phylogeography, climate‐based ENMs and multivariate statistical approaches to infer the effects of late Quaternary climate change on contemporary genetic variation of valley oak (Quercus lobata Née). ENMs indicated that valley oak maintained a stable distribution with local migration from the last interglacial period (~120 ka) to the Last Glacial Maximum (~21 ka, LGM) to the present compared with large‐scale range shifts for an eastern North American white oak (Quercus alba L.). Coast Range and Sierra Nevada foothill populations diverged in the late Pleistocene before the LGM [104 ka (28–1622)] and have occupied somewhat distinct climate niches, according to ENMs and coalescent analyses of divergence time. In accordance with neutral expectations for stable populations, nuclear microsatellite diversity positively correlated with niche stability from the LGM to present. Most strikingly, nuclear and chloroplast microsatellite variation significantly correlated with LGM climate, even after controlling for associations with geographic location and present climate using partial redundancy analyses. Variance partitioning showed that LGM climate uniquely explains a similar proportion of genetic variance as present climate (16% vs. 11–18%), and together, past and present climate explains more than geography (19%). Climate can influence local expansion–contraction dynamics, flowering phenology and thus gene flow, and/or impose selective pressures. These results highlight the lingering effect of past climate on genetic variation in species with stable distributions.     

Differential effects of host plant hybridization on herbivore community structure and grazing pressure on forest canopies

Interspecific hybridization in plants is known to have ecological effects on associated organisms. We examined the differences in insect herbivore community structure and grazing pressure on tree canopy leaves among natural hybrids and their parental oak species. We measured leaf traits, herbivore community structure, and grazing pressure on leaves of two oak species, Quercus crispula and Q. dentata, and their hybrids. The concentration of nitrogen in canopy leaves was greater in hybrids and in Q. dentata than in Q. crispula. The concentration of total phenolics was lower in hybrids than in Q. crispula. The concentration of condensed tannin was greater in hybrids than in Q. crispula. Relative herbivore abundance and species richness were greater on oak hybrids than on either parental species; herbivore species diversity and composition on hybrids were close to those on Q. crispula. Herbivore grazing pressure was lower on hybrids and Q. dentata than on Q. crispula. There was a negative correlation between herbivore grazing pressure and leaf nitrogen, suggesting that interspecific variation among oak taxa in herbivore pressure may be explained by leaf nitrogen; variation in herbivore community structure among oak taxa is likely to be controlled by polygenic leaf traits. Differing responses of (1) herbivore community structure and (2) herbivore grazing pressure to host plant hybridization may play important roles in regulating herbivore biodiversity in cool‐temperate forest canopies.     

Molecular phylogeny reveals the presence of cryptic speciation within Erysiphe japonica (≡ Typhulochaeta japonica), a powdery mildew on Quercus spp.

Molecular phylogenetic analyses based on 28S, ITS and IGS rDNA sequences indicate that Erysiphe japonica (≡ Typhulochaeta japonica) consists of two different genetic groups, one group on Quercus aliena, Q. robur and Q. serrata, and another group on Q. crispula var. crispula and Q. crispula var. horikawae. As morphological difference between the two groups are not yet marked distinctly, we suppose that the process of speciation has not yet been finished and propose a new variety, E. japonica var. crispulae, for the latter group based on the difference in host range and the distinct genetic segregation. Quercus robur (pedunculate oak) is a new host of E. japonica.     

High genetic diversity and stable Pleistocene distributional ranges in the widespread Mexican red oak Quercus castanea Née (1801) (Fagaceae)

The Mexican highlands are areas of high biological complexity where taxa of Nearctic and Neotropical origin and different population histories are found. To gain a more detailed view of the evolution of the biota in these regions, it is necessary to evaluate the effects of historical tectonic and climate events on species. Here, we analyzed the phylogeographic structure, historical demographic processes, and the contemporary period, Last Glacial Maximum (LGM) and Last Interglacial (LIG) ecological niche models of Quercus castanea, to infer the historical population dynamics of this oak distributed in the Mexican highlands. A total of 36 populations of Q. castanea were genotyped with seven chloroplast microsatellite loci in four recognized biogeographic provinces of Mexico: the Sierra Madre Occidental (western mountain range), the Central Plateau, the Trans‐Mexican Volcanic Belt (TMVB, mountain range crossing central Mexico from west to east) and the Sierra Madre del Sur (SMS, southern mountain range). We obtained standard statistics of genetic diversity and structure and tested for signals of historical demographic expansions. A total of 90 haplotypes were identified, and 29 of these haplotypes were restricted to single populations. The within‐population genetic diversity was high (mean h_S = 0.72), and among‐population genetic differentiation showed a strong phylogeographic structure (N_ST = 0.630 > G_ST = 0.266; p < .001). Signals of demographic expansion were identified in the TMVB and the SMS. The ecological niche models suggested a considerable percentage of stable distribution area for the species during the LGM and connectivity between the TMVB and the SMS. High genetic diversity, strong phylogeographic structure, and ecological niche models suggest in situ permanence of Q. castanea populations with large effective population sizes. The complex geological and climatic histories of the TMVB help to explain the origin and maintenance of a large proportion of the genetic diversity in this oak species.     

The effects of herbivory and soil fertility on the growth patterns of <Emphasis Type="Italic">Quercus serrata</Emphasis> and <Emphasis Type="Italic">Q. crispula</Emphasis> saplings at the shoot and individual levels

We investigated the combined effects of herbivore damage and soil fertility on shoot growth patterns of Quercus serrata and Q. crispula saplings at both the shoot and individual levels. Saplings were grown in herbivore-damaged or undamaged areas of the greenhouse with the two fertilization treatment levels, low or high. We measured the leaf area loss, number of flushes, length of extension units (EUs; the first vs the higher), number of leaves on each individual, and number of EUs. At the shoot level, the leaf area loss at high soil fertility was significantly greater than that at low soil fertility among the highest EUs of Q. serrata, while this difference was not significant in Q. crispula, suggesting that effect of soil fertility on leaf area loss is species-specific. Furthermore, herbivore damage was associated with a significant increase in the number of EUs and a reduction in the length of the higher EUs under both soil fertility treatments, although saplings had a tendency to produce significantly more flushes and longer individual EUs under the high soil fertility. At the individual level, herbivore-damaged saplings exhibited a significant increase in leaf numbers; however, the total length of the EUs in Q. serrata or Q. crispula was not significantly affected by herbivore damage, regardless of soil fertility. These results suggest that Q. serrata and Q. crispula saplings produce shorter EUs in response to herbivore damage in order to reduce the cost of mechanical support and spread the risk for any subsequent herbivore damage.     

Phylogeographical variation of chloroplast DNA in holm oak (Quercus ilex L.)

Variation in the lengths of restriction fragments (RFLPs) of the whole chloroplast DNA molecule was studied in 174 populations of Quercus ilex L. sampled over the entire distribution of this evergreen and mainly Mediterranean oak species. By using five endonucleases, 323 distinct fragments were obtained. From the 29 and 17 cpDNA changes identified as site and length mutations, respectively, 25 distinct chlorotypes were distinguished, mapped and treated cladistically with a parsimony analysis, using as an outgroup Q. alnifolia Poech, a closely related evergreen oak species endemic to Cyprus where Q. ilex does not grow. The predominant role of Q. ilex as maternal parent in hybridization with other species was reflected by the occurrence of a single very specific lineage of related chlorotypes, the most ancestral and recent ones being located in the southeastern and in the northwestern parts of the species’ geographical distribution, respectively. The lineage was constituted of two clusters of chlorotypes observed in the ‘ilex’ morphotyped populations of the Balkan and Italian Peninsulas (including the contiguous French Riviera), respectively. A third cluster was divided into two subclusters identified in the ‘rotundifolia’ morphotyped populations of North Africa, and of Iberia and the adjacent French regions, respectively. Postglacial colonization probably started from three distinct southerly refugia located in each of the three European peninsulas, and a contact area between the Italian and the Iberian migration routes was identified in the Rhône valley (France). Chlorotypes identical or related to those of the Iberian cluster were identified in the populations from Catalonia and the French Languedoc region, which showed intermediate morphotypes, and in the French Atlantic populations which possessed the ‘ilex’ morphotype, suggesting the occurrence of adaptive morphological changes in the northern part of the species’ distribution.     

Phylogeographical structure in Zelkova serrata in Japan and phylogeny in the genus Zelkova using the polymorphisms of chloroplast DNA

The genetic differentiation inherent in Zelkova serrata in Japan and the southern portion of the Korean Peninsula was examined by comparing a chloroplast DNA (cpDNA) sequence over a 16?k baselength in 40 individual samples collected from an area covering the natural distribution range of Z. serrata in Japan and the southern portion of the Korean Peninsula. We detected over 50 single-nucleotide polymorphisms in the protein-coding and intergenic regions, and over 30 insertions/deletions in the intergenic region. From the polymorphisms detected in the cpDNA, 14 haplotypes were identified. These 14 haplotypes had cluster-like structures and genetic differentiation between the clusters was large. Closely related haplotypes existed in adjacent regions. One haplotype existed in both Japan and the Korean Peninsula. By comparison with other Zelkova species, Z. serrata is apparently distinct from European and East Asian Zelkova species and Z. serrata is closest to the Ulmus species in the genus Zelkova. The effects of the analyzed length of the cpDNA sequence on the detection of polymorphisms were analyzed by re-sampling simulation.     

Fungal Community Analyses of Endophytic Fungi from Two Oak Species,Quercus mongolica and Quercus serrata,in Korea

Fungal endophytes have been recorded in various plant species with a richness of diversity, and their presence plays an essential role in host plant protection against biotic and abiotic stresses. This study applied the Illumina MiSeq sequencing platform based on the amplification of fungal ribosomal ITS2 region to analyze fungal endophytic communities of two oak species (Quercus mongolica and Q. serrata) with different oak wilt disease susceptibilities in Korea. The results showed a total of 230,768 sequencing reads were obtained and clustered at a 97% similarity threshold into 709 operational taxonomic units (OTUs). The OTUs of Q. serrata were higher than that of Q. mongolica with the number of 617 OTUs and 512 OTUs, respectively. Shannon index also showed that Q. serrata had a significantly higher level of fungal diversity than Q. mongolica. Total of OTUs were assigned into 5 fungal phyla, 17 classes, 60 orders, 133 families, 195 genera, and 280 species. Ascomycota was the dominant phylum with 75.11% relative abundance, followed by Basidiomycota with 5.28%. Leptosillia, Aureobasidium and Acanthostigma were the most abundant genera detected in Q. serrata with the average relative abundance of 2.85, 2.76, and 2.19%, respectively. On the other hand, Peltaster, Cladosporium and Monochaetia were the most common genera detected in Q. mongolica with the average relative abundance of 4.83, 3.03, and 2.87%, respectively. Our results indicated that fungal endophytic communities were significantly different between two oak species and these differences could influence responses of host trees to oak wilt disease caused by Raffaelea quercus-mongolicae.     

Identification of glacial refugia in south-eastern North America by phylogeographical analyses of a forest understorey plant, Trillium cuneatum

Aim We examine several hypotheses emerging from biogeographical and fossil records regarding glacial refugia of a southern thermophilic plant species. Specifically, we investigated the glacial history and post‐glacial colonization of a forest understorey species, Trillium cuneatum. We focused on the following questions: (1) Did T. cuneatum survive the Last Glacial Maximum (LGM) in multiple refugia, and (if so) where were they located, and is the modern genetic structure congruent with the fossil record‐based reconstruction of refugia for mesic deciduous forests? (2) What are the post‐glacial colonization patterns in the present geographical range? Location South‐eastern North America. Methods We sampled 45 populations of T. cuneatum throughout its current range. We conducted phylogeographical analyses based on maternally inherited chloroplast DNA (cpDNA haplotypes) and used TCS software to reconstruct intraspecific phylogeny. Results We detected six cpDNA haplotypes, geographically highly structured into non‐overlapping areas. With one exception, none of the populations had mixed haplotype composition. TCS analysis resulted in two intraspecific cpDNA lineages, with one clade subdivided further by shallower diversification. Main conclusions Our investigation revealed that T. cuneatum survived the LGM in multiple refugia, belonging to two (western, eastern) genealogical lineages geographically structured across south‐eastern North America. The western clade is confined to the south‐western corner of T. cuneatum’s modern range along the Lower Mississippi Valley, where fossil records document a major refugium of mesic deciduous forest. For the eastern clade, modern patterns of cpDNA haplotype distribution suggest cryptic vicariance, in the form of forest contractions and subsequent expansions associated with Pleistocene glacial cycles, rather than simple southern survival and subsequent northward colonization. The north–south partitioning of cpDNA haplotypes was unexpected, suggesting that populations of this rather southern thermophilic species may have survived in more northern locations than initially expected based on LGM climate reconstruction, and that the Appalachian Mountains functioned as a barrier to the dispersal of propagules originating in more southern refugia. Furthermore, our results reveal south‐west to north‐east directionality in historical migration through the Valley and Ridge region of north‐west Georgia.     

Leaf morphological analyses in four European oak species (Quercus) and their hybrids: A comparison of traditional and geometric morphometric methods

Abstract

In this study, leaf morphology was assessed in a mixed oak stand (western France) using two geometric morphometric (landmark and outline) datasets and one dataset of 19 leaf measures. Adult oaks (817 oaks), comprising four white oak species (Quercus petraea, Q. robur, Q. pubescens and Q. pyrenaica), were sampled for DNA extraction and genetic analysis (nuclear microsatellites). Leaf morphology was assessed on 336 oaks, comprising pure species and hybrids as determined by genetic assignment. This comparative study of oak leaf morphology, based on the use of two free size geometric morphometric methods and a set of leaf measurements, combined with the genetic assignment of individuals to pure species or hybrids, provided information about the differences among species and the intermediate leaf morphology of their hybrids.     

Early seedling growth of pine (Pinus densiflora) and oaks (Quercus serrata,Q. mongolica,Q. variabilis) in response to light intensity and soil moisture

In climatic chambers seed germination and seedling growth of Pinus densiflora Sieb. et Zucc., Quercus serrata Thunb., Quercus mongolica Fisch. ex Turcz. and Quercus variablilis Bl. were investigated as functions of light intensity and soil moisture. In Korea these tree species occur widely and form mixed forests with different distributions. Species clearly differed in the pattern of germination and early seedling growth between light and soil treatments. The germination of pine did not differ between the experimental treatments until the breaking of the primary buds. After that, light intensity was the deciding factor for further development. In the most moist treatment (approx. field capacity) growth of the pine seedlings was strongly inhibited. For the three oak species, differences between experimental treatments first occurred after complete formation of primary leaves. Seed development strongly correlated with the weight of the acorn. The large seeded Q. variabilis (acorns with mean weight of 4.7 g) developed faster and reached larger dimensions towards the end of the experiment than Q. mongolica (2.8 g per acorn) and Q. serrata (0.9 g per acorn). Regarding height and biomass growth, the oak species showed a higher shade tolerance than pine. The proleptic shoot growth was clearly influenced by the light intensity. Root formation was favoured by a high exposure to light. In case of the oak species reduction of soil moisture increased the length of primary roots and the number of secondary roots.     

Genetic variation and differentiation within a natural community of five oak species (Quercus spp.)

Chloroplast DNA and two categories of nuclear markers - isozymes and microsatellites - were used to examine a very rich natural community of oaks (Quercus spp.) situated in west-central Romania. The community consists of five oak species: Q. robur, Q. petraea, Q. pubescens, and Q. frainetto - that are closely related -, and Q. cerris. A total of five chloroplast haplotypes was identified. Q. cerris was fixed for a single haplotype. The other four species shared the two most common haplotypes. One haplotype was confined to Q. robur and a very rare one was restricted to Q. petraea. Both types of nuclear markers revealed a larger genetic variation for Q. pubescens and Q. petraea than for Q. frainetto and Q. robur, although the differences between species are in most cases not significant. At the nuclear level, Q. cerris could be clearly separated from the other four oak species confirming the taxonomic classification. Regardless of the estimate used, the levels of polymorphism revealed by microsatellites were much higher than those based on isozymes. For the four closely related species the overall genetic differentiation was significant at both categories of nuclear markers. Several loci, such as Acp-C for isozymes, and ssrQpZAG36 and ssrQrZAG96 for microsatellites were very useful to discriminate among species. However, the level of differentiation varied markedly between pairs of species. The genetic affinities among the species may reflect different phylogenetic distances and/or different rates of recurrent gene flow at this site.     

Assessment of within-population genetic structure in <Emphasis Type="Italic">Quercus crispula</Emphasis> and <Emphasis Type="Italic">Q. dentata</Emphasis> by amplified fragment length polymorphism analysis

The spatial genetic structure was assessed by amplified fragment length polymorphism (AFLP) technique for Quercus crispula Blume (37 individuals) and Q.dentata Thunberg (54 individuals) growing along a 550-m transect in Ishikari, Hokkaido, northern Japan. The simple Mantel test revealed significant negative correlation between genetic similarity and geographic distance in Q.dentata and negative but insignificant correlation in Q.crispula. Spatial autocorrelation analysis revealed that Mantels r generally decreased from positive to negative values with the increase of spatial distance in both oak species with significant deviation from zero for the 50–100-m (positive) and 250–300-m classes (negative) in Q.dentata. Thus, significant spatial genetic structure was revealed at least in Q.dentata