BACKGROUND
natural hybridization and subsequent introgression provide opportunities for alleles to travel between species within a genus  <cit> .
the presence of new genetic variants causes an increase of genetic diversity and may also bring novel adaptations in the recipient species  <cit> .
most of the studies on plant natural hybridization focused on one pair of species and only in rare cases the occurrence of this phenomenon was investigated in a three or more species context.
the hybridization patterns can be more accurately estimated when another related species occurring in the proximity of a two-species hybrid zone are considered  <cit> .
the availability of highly polymorphic genetic markers  and powerful statistical procedures  has increased our ability to identify first and later-generation hybrids in a wide range of plant species.
among the long-lived species of the temperate zone, poplars  and oaks  <cit>  have received much of the attention.
genetic exchange within the genus quercus was suspected for a long time  <cit>  and confirmed by numerous recent studies which used both chloroplast and nuclear genetic markers  <cit> .
hybridization appears to be influenced by the spatial structure of the population  <cit>  and relative frequency of the parental species  <cit> .
however, a recent study that used a set of microsatellite loci to study pairs of populations distributed across the natural range suggested the existence of reproductive isolation rather than high rates of gene flow between q. robur and q. petraea  <cit> .
paternity analysis studies have shown that wind-pollinated tree species such as oaks can spread their pollen over long distances  <cit> .
despite the high levels of distant gene flow reported, matings occur commonly among neighboring trees determining significant differences in the genetic composition of the effective pollen pools received by each maternal tree in a single year  <cit>  but also by the same maternal trees among reproductive years  <cit> .
paternity analyses based on microsatellite markers were used for estimating hybridization rates at a fine scale in mixed oak stands  <cit> .
pollination experiments conducted over several years revealed that hybridization between the most common oak species in europe, q. robur and q. petraea  <cit> , and between a first-generation hybrid and these two parental species  <cit>  is possible.
pollination of q. robur seed trees by q. petraea pollen is usually more successful than the reciprocal cross  <cit> .
here, we investigate the patterns of actual hybridization in a mixed natural forest comprising all white oak species which occur in the region.
the specific objectives of this study were to:  find any evidence of interspecific matings  determine the directionality of gene flow between species  quantify the genetic differences between the effective pollen pools and to correlate these differences with the physical distances between the maternal parents  compare the proportion of hybrids present in acorn collections with the proportion of hybrids among the adult trees of the forest.
methods
species sampled and study area
oaks are monoecious species, i.e.
they have male and female flowers on the same tree.
their annual regularity of flower production contrasts with their irregular fructification, which may range from abundant crops  in some years to poor or no crops in others  <cit> .
the distance between abundant crops may reach 6- <dig> or even more years for some oak species  <cit>  which makes acorn collection very difficult.
the four closely related species included in this study, q. robur, q. petraea, q. pubescens and q. frainetto, respectively, are deciduous trees and belong to section quercus s.s. .
the study plot is a part of the bejan oak reserve  which is situated in the carpathian basin, west-central romania.
more details about the study area are given elsewhere  <cit> .
adult sampling
a total of  <dig> adult trees were morphologically and genetically analyzed in a previous study  <cit> .
the trees are about 120- <dig> yr old and originated from natural regeneration.
most of the individual trees are located within the core plot, an area where the sampling of all white oak trees was exhaustive.
as the sample size for q. frainetto within the core plot was very small  the sampling for this species was extended to the nearby area .
based on pubescence and leaf descriptors the vast majority of individual trees  were classified as belonging to a morphological species  <cit> .
furthermore, based on a set of  <dig> microsatellite and  <dig> isozyme markers and using a bayesian admixture analysis implemented in the program structure version  <dig>   <cit>  without any a priori information on morphology we assigned adult individuals to pure species and hybrids  <cit> .
there was a very good correspondence of each morphological species  with one of the inferred genetic clusters .
in this study we classified each tree with a high probability to belong to one of the four genetic clusters  as pure species and q <  <dig>  as hybrid .
however, as in the recent study of lepais and his colleagues  <cit>  individuals with q <  <dig>  for one genetic cluster , and q <  <dig>  for the remaining three clusters  were considered as pure species.
offspring sampling
we genotyped  <dig> acorns collected in autumn  <dig> from eight seed parents: four q. petraea, two q. frainetto, one q. robur, and one q. pubescens, respectively.
according to the genetic assignment of the adult tree generation,  <dig> out of  <dig> seed trees were classified as pure species  <cit> .
only one seed parent  was morphologically classified as q. frainetto but genetically was considered as hybrid  between 'frainetto' and 'pubescens' clusters.
the sampling design was affected by the nearly complete lack of fructification for two species, q. robur and q. pubescens, in  <dig> at our study site.
all maternal trees are located in the central part of the study plot .
the acorns were collected on the ground near the trunk of each seed parent after shooting the crown with a gun.
dna isolation and microsatellite typing
details of the method used for dna extraction, pcr amplification and microsatellite genotyping of the adult trees are given elsewhere  <cit> .
using the same procedure the genotypes of the offspring were scored at six microsatellite loci.
loci ssrqpzag1/ <dig>  ssrqpzag <dig>  ssrqpzag <dig> and ssrqpzag <dig> were originally developed for q. petraea  <cit> , ssrqrzag <dig> was developed for q. robur  <cit>  and msq <dig> for q. macrocarpa  <cit> .
statistical data analysis
paternity analysis  was performed using the likelihood method implemented in the program famoz  <cit> .
a simulation assuming both parents inside and outside the study plot was done using  <dig>   <dig> offspring to estimate the threshold for paternity .
since mistyping is very likely to occur when scoring microsatellite alleles, the error rate was set to  <dig>  - both in the simulation and in the assignment of the most likely father.
the genotype with the highest log  of the odds ratio  score above the threshold value was considered as the most likely pollen parent.
the lod score represents the likelihood of a particular genotype being the father compared to all other genotypes.
the male gametic contribution  was inferred by subtracting the female contribution from the offspring genotype.
in ambiguous cases, in which the seed parent and its offspring were heterozygotes and shared the same pair of alleles at one locus, the pollen profile was constructed using the maximum likelihood estimator of pollen parent contribution .
we assumed that an unknown pollen parent from outside the study plot  was successful, if none of the investigated trees could produce a pollen with the inferred multilocus paternal haplotype, for example due to the occurrence of an allele in the progeny which was not observed for any adult tree.
successful pollination by external pollen cannot be ruled out even if one or more potential pollen parents within the study plot were identified and the most likely pollen parent was selected based on its lod score.
thus, the reported proportions of external pollen are minimum estimates.
however, strong deviations from these estimates to the 'true' rate of successful external pollen are unlikely due to high exclusion probabilities.
to quantify the genetic differentiation among the effective pollen pools of the seed trees, an analysis of molecular variance  of male haplotypes was performed using the software package genalex version  <dig>   <cit> .
the significance of the Φpt, an analogous of fst, was tested by  <dig> random permutations.
this analysis was conducted for the total effective pollen pools and for the separate effective pollen pools originating from inside and outside of the study plot.
by using the same software we then tested the isolation by distance hypothesis.
the association between the two matrices, pairwise Φpt between the effective pollen pools of single tree progenies and pairwise physical distances between the maternal parents, was analyzed using the mantel test statistics  <cit> .
statistical significance was tested by means of  <dig> permutations.
RESULTS
paternity analysis
the six microsatellite markers were highly polymorphic, resulting in high exclusion probabilities for paternity analysis .
only a few alleles that were not present in the adult trees were found among the offspring .
strong evidence for a null  allele was found at locus ssrqpzag1/ <dig> in q. petraea.
mismatches between one seed parent  and its offspring were observed at this locus whereas the same offspring shared at least one allele with the seed parent at five remaining loci.
moreover, the seed parent as well as all the mismatching offspring appeared to be homozygotes, supporting the hypothesis that both seed parent and mismatching offspring are heterozygous for a null allele.
in accordance with the presence of null alleles, the greatest deficit of heterozygotes  amongst the six microsatellite loci in the four oak species was observed at the locus ssrqpzag1/ <dig> in q. petraea adult trees  <cit> .
number of alleles in the adult tree generation , number of new alleles detected among the offspring , size range of alleles and exclusion probability for six microsatellite loci.
the pollen donors that were assigned to the offspring are shown exemplarily in figure 2a-d for four maternal trees, one of each oak species.
among  <dig> offspring having their pollen donor within the study plot only three appeared to be self-pollinated.
averaged across all offspring and seed parents, selfing rate was very low  and due to selfing in one single q. frainetto seed parent .
mean pollen dispersal distances are given for each seed parent in table  <dig>  the observed values differ considerably among seed parents ranging from  <dig>  m to  <dig>  m. spatial distribution of each oak species within the plot, e.g.
the presence or absence of conspecific trees in the neighborhood of each seed parent seems to affect markedly the dispersal distances .
the gene flow from outside the study plot was estimated at 48% when averaged across all seed parents.
however, this rate varied markedly among the seed parents although they are located relatively close to each other .
for instance, relatively low rates of apparent gene flow from outside were observed for the two q. frainetto  and one q. pubescens  maternal parents compared to the q. robur maternal parent .
the proportion of gene flow originating from outside substantially varied across q. petraea seed parents .
species names were abbreviated as follows: rob - q. robur; pet - q. petraea; pub - q. pubescens; fra - q. frainetto.
hybridization between oak species
the majority of acorns gathered from the seven pure species seed parents were fertilized by conspecific pollen .
eleven offspring  appear to be first-generation hybrids since their parents belong to different pure species .
seven were f <dig> hybrids between q. petraea  and q. pubescens , three between q. robur and q. petraea, and one between q. frainetto and q. pubescens, respectively.
even though offspring from four pure q. petraea mother trees were investigated, all first-generation hybrids between q. petraea and q. pubescens were detected among the acorns collected from one q. petraea tree, pet- <dig> .
the presence of putative hybrids  among effective pollen donors was interpreted as introgressive hybridization in  <dig> out of  <dig> mating events .
putative hybrids between q. petraea and q. pubescens  and one hybrid between q. pubescens and q. frainetto  fertilized q. petraea mother trees.
one hybrid between q. petraea and q. pubescens  fertilized four ovules of the q. pubescens mother tree.
moreover, one hybrid between q. robur and q. pubescens fertilized seven ovules of the same q. pubescens tree.
three different hybrids between q. frainetto and q. pubescens  were observed among the pollen donors of the q. frainetto seed tree, fra- <dig>  none of the q. frainetto trees situated outside the core plot  was found among the effective pollen donors of the q. frainetto seed tree .
one putative hybrid between q. petraea and q. robur  fertilized one ovule of the q. robur mother tree.
among  <dig> offspring of the single mother tree, fra- <dig>  which was considered as hybrid between q. frainetto and q. pubescens  in spite its q. frainetto physical appearance, we detected as pollen donors both pure parental species  and putative hybrids  between this pair of species.
the total pollen flow between pure species as well as between pure species and their hybrids, detected in our sample, is illustrated in figure  <dig>
differentiation among the effective pollen pools
the genetic differentiation  among the eight effective pollen pools averaged over all loci was significant .
pairwise Φpt values between the effective pollen pools of eight single-tree progenies varied from  <dig>  to  <dig>  .
no significant differences were detected among the effective pollen pools of three q. petraea seed parents.
only the effective pollen pool of the highly hybridizing q. petraea tree, pet- <dig>  differed significantly from the effective pollen pools of other q. petraea and q. pubescens seed trees.
interestingly, there was no significant difference between the effective pollen pools received by the q. robur seed parent, rob- <dig>  and one q. petraea tree, pet- <dig>  there was almost no genetic difference between the effective pollen pools of the q. frainetto seed parent and the putative hybrid between q. frainetto and q. pubescens  .
species names were abbreviated as follows: rob - q. robur; pet - q. petraea; pub - q. pubescens; fra - q. frainetto.
ns p >  <dig> ; * p <  <dig> ; ** p <  <dig> ; *** p <  <dig> .
p values were estimated after  <dig> permutations.
when considering the bulked data over eight seed parents the genetic differences between the effective pollen pools from inside and outside the study plot were significant across all loci .
however, at two loci, ssrqpzag1/ <dig> and ssrqpzag <dig>  the differences were not significant .
our data provided no support for a strong relationship between the pairwise Φpt values and spatial distances .
even with regard to the limited sample size this result suggests that spatial distance between the seed parents had no significant impact on the differentiation of the effective pollen pools received by each seed parent.
discussion
the present study provides new insights into pollen mediated-gene flow and hybridization within the european white oak species complex.
mismatches between putative parents and offspring due to the presence of a suspected null allele were reported in other studies  and seem to be relatively common when employing microsatellite markers in parentage analyses .
selfing was low, as expected in monoecious wind-pollinated tree species of the temperate zone such as oaks .
the selfed offspring originate from one q. frainetto seed parent, the species with the lowest density in the area  <cit> .
according to streiff and her colleagues  <cit>  the effective pollen dispersal in oaks involves two components: local dispersal and long-distance transport.
in our study, the portion of local pollen dispersal was quite variable among the seed parents.
for instance, the local component of dispersal seems to be prevalent in q. frainetto and q. pubescens seed parents, the species with the lowest relative abundance within the whole forest reserve  <cit> .
in contrast, the amount of pollen originating from outside the study plot was much higher for the three q. petraea seed parents  which is consistent with other results obtained via paternity analyses in different oak species  <cit> .
according to the forest records, q. petraea is the predominant species in the reserve and in the nearby forest which may explain the larger amount of pollen coming from outside the study plot for this species.
by using paternity analyses we studied pollination patterns in a very complex situation involving all species which occur not only in the mixed stand but also in the nearby area on tens of kilometers.
the importance of considering all species when studying hybridization was very well highlighted in france where hybrids were detected in considerable proportions in a mixed stand in the absence of parental species  <cit> .
our findings suggest the existence of reproductive barriers between oak species since the majority of the effective pollen originates from the same species.
however, the hybridization rate was estimated at  <dig> % in the seed generation.
we found strong evidence for the production of first generation hybrids  because both parents were phenotypically and genetically pure species.
we observed f <dig> hybrids between three out of twelve possible species combinations, including reciprocal crossings, which suggest that hybridization at this site is not restricted to one species pair .
this result is strongly influenced by our unbalanced sample of mother trees of different species.
previously, in a two-oak-species stand streiff and her colleagues  <cit>  reported hybridization in  <dig> out of  <dig> offspring .
however, this rate of interspecific gene flow was estimated considering only the physical appearance of oak trees.
as in our study the hybridization rate varied markedly among maternal parents.
those oaks showing high levels of hybridization were mostly surrounded by individuals belonging to other species.
in such contact zones between species, hybridization may occur due to physical proximity  <cit> .
hybridization was not restricted to the formation of first-generation hybrids since  <dig> % of the pollen donors were putative introgressed forms between the oak species.
nearly 3/ <dig> of the total number of hybrids were later-generation hybrids.
the predominance of backcrossed individuals in the detriment of first generation hybrids was reported for the same species complex in western europe  <cit> .
by comparing the overall proportion of hybrids in the adult generation  with the proportion of hybrids in the offspring  we found nearly 80% more hybrids in the seed generation.
the decrease in the number of hybrids from offspring to adults as well as the maintenance of species identity may be the result of the selection against hybrids  <cit>  which acts in the ecological context of this forest.
environmental selection against hybrids was invoked to explain the overestimation of hybridization rates when analyzing seeds  <cit> .
another explanation for this situation would be the stand history.
the uneven-aged structure of the adult population might indicate that some of the trees were established when less than four species were found at this site.
the higher proportion of hybrid acorns compared to hybrids adult trees may also be caused by past introgression events which have reduced the strength of the reproductive barriers over time  <cit> .
our results indicate that hybrids participate at reproduction both as pollen donor and receptor .
similar situations were reported in a mixed stand with q. petraea, q. pubescens and intermediate individuals between both species in italy  <cit> .
moreover, hybrids between one pair of species prefer to mate with one of their parental species.
only in two instances  we found hybrids which mate with other species than the parental ones.
these preferential matings indicate weaker reproductive barriers between the parental species and their hybrids as compared to other related species and hybrid combinations.
the asymmetric gene flow from q. petraea to q. robur revealed by paternity analysis in our sample is consistent with the outcome of the bayesian analysis in the adult generation  <cit> .
a possible explanation is the contrasting relative abundance of the two species at bejan: q. petraea is very common in the study area while q. robur is frequent only in a small marginal part of the reserve.
higher success rates of the combination q. robur  × q. petraea  were also observed in artificial crosses  <cit> .
our finding supports the model of asymmetrical hybridization between q. petraea and q. robur  <cit> .
moreover, the gene exchange between q. pubescens and q. petraea was mainly from q. pubescens towards q. petraea  in our sample.
the opposite direction was observed between the two species in a hybrid zone in italy  <cit> .
the ecological conditions and forest composition which differ between the two geographical areas may have influenced the direction of introgression  <cit> .
these discrepancies between hybridization patterns described throughout oaks' distribution range stress the considerable role played by ecological factors and the necessity for comparative studies.
the level of hybridization may be higher for one individual or another because of an overlap in flowering time with most of the individuals of other species.
almost no differences in flowering time between oak species have been reported in different mixed forests  <cit> .
for instance, large differences in flowering time between the single q. robur seed parent and the majority of q. robur adult trees and an overlap of the flowering time with most of the q. petraea individuals can be advocated as one of the causes to explain the high propensity for hybridization of this seed tree.
differences in flowering time may lead to a preponderance of pollen coming from q. petraea relative to the pollen coming from its conspecific individuals.
phenological observations at this site are needed to validate this hypothesis.
furthermore, the present results are based on only one flowering season and the patterns of pollen movement may certainly vary from year to year  <cit> .
hence, the physical barrier represented by q. cerris trees might have been also hindered the pollen coming from the eastern part of the study plot, where q. robur is predominant .
the low amount of genetic differentiation among the effective pollen pools for seed parents which belong to the same species may be due to the limited number of sampled maternal parents.
genetic heterogeneity of effective pollen pools was reported in other studies based on larger sample sizes  <cit> .
genetic heterogeneity of effective pollen pools may be caused by uneven reproductive contribution of local pollen donors  <cit> .
differences in flowering phenology have also been mentioned as a possible cause  <cit> .
however, the differentiation of effective pollen pools for seed parents which belong to different species was highly significant, except for the effective pollen pool received by the single q. robur seed parent which was very similar to the effective pollen pool of q. petraea trees.
the low amount of genetic differentiation between the effective pollen pool received by the q. robur seed parent and the three other q. petraea neighboring trees is consistent with the result of the paternity analysis that indicated a high rate of hybridization with q. petraea trees.
in contrast to the study of streiff and her colleagues  <cit>  a significant differentiation between the effective pollen pools from outside and inside the study plot was detected.
this genetic heterogeneity may be caused by the small number of pollen donors identified within the study plot that were not representative for the genetic composition of the whole oak forest.
our data did not follow an isolation by distance pattern like in other tree species  <cit> .
the high number of oak species which are present among the sampled maternal trees correlated with their preference for conspecific pollen rather than heterospecific pollen may explain this pattern.
CONCLUSIONS
by using paternity analysis and a set of six highly polymorphic microsatellite loci we have traced patterns of contemporary gene flow within and between oak species in a natural reserve.
our results confirm that introgressive hybridization between different pairs of species is quite common and play a significant evolutionary role within the white oak species complex.
selection against hybrids over years might contribute to the maintenance of species identity despite considerable amounts of interspecific gene flow.
however, the differences between the proportion of hybrid acorns and the proportion of hybrid adults may also be explained by the stand history or past introgression.
further phenological investigations and a larger sample size are required for a better understanding of the hybridization events taking place in this four-oak-species stand.
authors' contributions
alc conceived the idea, collected the seed samples, did the laboratory work, performed the paternity analysis, and wrote the majority of the text.
og helped to analyze the data and to draft the manuscript.
rf participated in the design and coordination of the study, data interpretation and made substantial contributions to writing this paper.

 natural hybridization and subsequent introgression provide opportunities for alleles to travel between species within a genus. the presence of new genetic variants causes an increase of genetic diversity and may also bring novel adaptations in the recipient species. most of the studies on plant natural hybridization focused on one pair of species and only in rare cases the occurrence of this phenomenon was investigated in a three or more species context. the hybridization patterns can be more accurately estimated when another related species occurring in the proximity of a two-species hybrid zone are considered. the availability of highly polymorphic genetic markers and powerful statistical procedures has increased our ability to identify first and later-generation hybrids in a wide range of plant species. among the long-lived species of the temperate zone, poplars and oaks have received much of the attention. genetic exchange within the genus quercus was suspected for a long time and confirmed by numerous recent studies which used both chloroplast and nuclear genetic markers. hybridization appears to be influenced by the spatial structure of the population and relative frequency of the parental species. however, a recent study that used a set of microsatellite loci to study pairs of populations distributed across the natural range suggested the existence of reproductive isolation rather than high rates of gene flow between q. robur and q. petraea. paternity analysis studies have shown that wind-pollinated tree species such as oaks can spread their pollen over long distances. despite the high levels of distant gene flow reported, matings occur commonly among neighboring trees determining significant differences in the genetic composition of the effective pollen pools received by each maternal tree in a single year but also by the same maternal trees among reproductive years. paternity analyses based on microsatellite markers were used for estimating hybridization rates at a fine scale in mixed oak stands. pollination experiments conducted over several years revealed that hybridization between the most common oak species in europe, q. robur and q. petraea, and between a first-generation hybrid and these two parental species is possible. pollination of q. robur seed trees by q. petraea pollen is usually more successful than the reciprocal cross. here, we investigate the patterns of actual hybridization in a mixed natural forest comprising all white oak species which occur in the region. the specific objectives of this study were to: find any evidence of interspecific matings determine the directionality of gene flow between species quantify the genetic differences between the effective pollen pools and to correlate these differences with the physical distances between the maternal parents compare the proportion of hybrids present in acorn collections with the proportion of hybrids among the adult trees of the forest. methods species sampled and study area oaks are monoecious species, i.e. they have male and female flowers on the same tree. their annual regularity of flower production contrasts with their irregular fructification, which may range from abundant crops in some years to poor or no crops in others. the distance between abundant crops may reach 6- or even more years for some oak species which makes acorn collection very difficult. the four closely related species included in this study, q. robur, q. petraea, q. pubescens and q. frainetto, respectively, are deciduous trees and belong to section quercus s.s.. the study plot is a part of the bejan oak reserve which is situated in the carpathian basin, west-central romania. more details about the study area are given elsewhere. adult sampling a total of adult trees were morphologically and genetically analyzed in a previous study. the trees are about 120- yr old and originated from natural regeneration. most of the individual trees are located within the core plot, an area where the sampling of all white oak trees was exhaustive. as the sample size for q. frainetto within the core plot was very small the sampling for this species was extended to the nearby area. based on pubescence and leaf descriptors the vast majority of individual trees were classified as belonging to a morphological species. furthermore, based on a set of microsatellite and isozyme markers and using a bayesian admixture analysis implemented in the program structure version without any a priori information on morphology we assigned adult individuals to pure species and hybrids. there was a very good correspondence of each morphological species with one of the inferred genetic clusters. in this study we classified each tree with a high probability to belong to one of the four genetic clusters as pure species and q < as hybrid. however, as in the recent study of lepais and his colleagues individuals with q < for one genetic cluster, and q < for the remaining three clusters were considered as pure species. offspring sampling we genotyped acorns collected in autumn from eight seed parents: four q. petraea, two q. frainetto, one q. robur, and one q. pubescens, respectively. according to the genetic assignment of the adult tree generation, out of seed trees were classified as pure species. only one seed parent was morphologically classified as q. frainetto but genetically was considered as hybrid between'frainetto' and'pubescens' clusters. the sampling design was affected by the nearly complete lack of fructification for two species, q. robur and q. pubescens, in at our study site. all maternal trees are located in the central part of the study plot. the acorns were collected on the ground near the trunk of each seed parent after shooting the crown with a gun. dna isolation and microsatellite typing details of the method used for dna extraction, pcr amplification and microsatellite genotyping of the adult trees are given elsewhere. using the same procedure the genotypes of the offspring were scored at six microsatellite loci. loci ssrqpzag1/ ssrqpzag ssrqpzag and ssrqpzag were originally developed for q. petraea, ssrqrzag was developed for q. robur and msq for q. macrocarpa. statistical data analysis paternity analysis was performed using the likelihood method implemented in the program famoz. a simulation assuming both parents inside and outside the study plot was done using offspring to estimate the threshold for paternity. since mistyping is very likely to occur when scoring microsatellite alleles, the error rate was set to - both in the simulation and in the assignment of the most likely father. the genotype with the highest log of the odds ratio score above the threshold value was considered as the most likely pollen parent. the lod score represents the likelihood of a particular genotype being the father compared to all other genotypes. the male gametic contribution was inferred by subtracting the female contribution from the offspring genotype. in ambiguous cases, in which the seed parent and its offspring were heterozygotes and shared the same pair of alleles at one locus, the pollen profile was constructed using the maximum likelihood estimator of pollen parent contribution. we assumed that an unknown pollen parent from outside the study plot was successful, if none of the investigated trees could produce a pollen with the inferred multilocus paternal haplotype, for example due to the occurrence of an allele in the progeny which was not observed for any adult tree. successful pollination by external pollen cannot be ruled out even if one or more potential pollen parents within the study plot were identified and the most likely pollen parent was selected based on its lod score. thus, the reported proportions of external pollen are minimum estimates. however, strong deviations from these estimates to the'true' rate of successful external pollen are unlikely due to high exclusion probabilities. to quantify the genetic differentiation among the effective pollen pools of the seed trees, an analysis of molecular variance of male haplotypes was performed using the software package genalex version. the significance of the Φpt, an analogous of fst, was tested by random permutations. this analysis was conducted for the total effective pollen pools and for the separate effective pollen pools originating from inside and outside of the study plot. by using the same software we then tested the isolation by distance hypothesis. the association between the two matrices, pairwise Φpt between the effective pollen pools of single tree progenies and pairwise physical distances between the maternal parents, was analyzed using the mantel test statistics. statistical significance was tested by means of permutations.  paternity analysis the six microsatellite markers were highly polymorphic, resulting in high exclusion probabilities for paternity analysis. only a few alleles that were not present in the adult trees were found among the offspring. strong evidence for a null allele was found at locus ssrqpzag1/ in q. petraea. mismatches between one seed parent and its offspring were observed at this locus whereas the same offspring shared at least one allele with the seed parent at five remaining loci. moreover, the seed parent as well as all the mismatching offspring appeared to be homozygotes, supporting the hypothesis that both seed parent and mismatching offspring are heterozygous for a null allele. in accordance with the presence of null alleles, the greatest deficit of heterozygotes amongst the six microsatellite loci in the four oak species was observed at the locus ssrqpzag1/ in q. petraea adult trees. number of alleles in the adult tree generation, number of new alleles detected among the offspring, size range of alleles and exclusion probability for six microsatellite loci. the pollen donors that were assigned to the offspring are shown exemplarily in figure 2a-d for four maternal trees, one of each oak species. among offspring having their pollen donor within the study plot only three appeared to be self-pollinated. averaged across all offspring and seed parents, selfing rate was very low and due to selfing in one single q. frainetto seed parent. mean pollen dispersal distances are given for each seed parent in table the observed values differ considerably among seed parents ranging from m to m. spatial distribution of each oak species within the plot, e.g. the presence or absence of conspecific trees in the neighborhood of each seed parent seems to affect markedly the dispersal distances. the gene flow from outside the study plot was estimated at 48% when averaged across all seed parents. however, this rate varied markedly among the seed parents although they are located relatively close to each other. for instance, relatively low rates of apparent gene flow from outside were observed for the two q. frainetto and one q. pubescens maternal parents compared to the q. robur maternal parent. the proportion of gene flow originating from outside substantially varied across q. petraea seed parents. species names were abbreviated as follows: rob - q. robur; pet - q. petraea; pub - q. pubescens; fra - q. frainetto. hybridization between oak species the majority of acorns gathered from the seven pure species seed parents were fertilized by conspecific pollen. eleven offspring appear to be first-generation hybrids since their parents belong to different pure species. seven were f hybrids between q. petraea and q. pubescens, three between q. robur and q. petraea, and one between q. frainetto and q. pubescens, respectively. even though offspring from four pure q. petraea mother trees were investigated, all first-generation hybrids between q. petraea and q. pubescens were detected among the acorns collected from one q. petraea tree, pet-. the presence of putative hybrids among effective pollen donors was interpreted as introgressive hybridization in out of mating events. putative hybrids between q. petraea and q. pubescens and one hybrid between q. pubescens and q. frainetto fertilized q. petraea mother trees. one hybrid between q. petraea and q. pubescens fertilized four ovules of the q. pubescens mother tree. moreover, one hybrid between q. robur and q. pubescens fertilized seven ovules of the same q. pubescens tree. three different hybrids between q. frainetto and q. pubescens were observed among the pollen donors of the q. frainetto seed tree, fra- none of the q. frainetto trees situated outside the core plot was found among the effective pollen donors of the q. frainetto seed tree. one putative hybrid between q. petraea and q. robur fertilized one ovule of the q. robur mother tree. among offspring of the single mother tree, fra- which was considered as hybrid between q. frainetto and q. pubescens in spite its q. frainetto physical appearance, we detected as pollen donors both pure parental species and putative hybrids between this pair of species. the total pollen flow between pure species as well as between pure species and their hybrids, detected in our sample, is illustrated in figure differentiation among the effective pollen pools the genetic differentiation among the eight effective pollen pools averaged over all loci was significant. pairwise Φpt values between the effective pollen pools of eight single-tree progenies varied from to. no significant differences were detected among the effective pollen pools of three q. petraea seed parents. only the effective pollen pool of the highly hybridizing q. petraea tree, pet- differed significantly from the effective pollen pools of other q. petraea and q. pubescens seed trees. interestingly, there was no significant difference between the effective pollen pools received by the q. robur seed parent, rob- and one q. petraea tree, pet- there was almost no genetic difference between the effective pollen pools of the q. frainetto seed parent and the putative hybrid between q. frainetto and q. pubescens. species names were abbreviated as follows: rob - q. robur; pet - q. petraea; pub - q. pubescens; fra - q. frainetto. ns p >; * p <; ** p <; *** p <. p values were estimated after permutations. when considering the bulked data over eight seed parents the genetic differences between the effective pollen pools from inside and outside the study plot were significant across all loci. however, at two loci, ssrqpzag1/ and ssrqpzag the differences were not significant. our data provided no support for a strong relationship between the pairwise Φpt values and spatial distances. even with regard to the limited sample size this result suggests that spatial distance between the seed parents had no significant impact on the differentiation of the effective pollen pools received by each seed parent. discussion the present study provides new insights into pollen mediated-gene flow and hybridization within the european white oak species complex. mismatches between putative parents and offspring due to the presence of a suspected null allele were reported in other studies and seem to be relatively common when employing microsatellite markers in parentage analyses. selfing was low, as expected in monoecious wind-pollinated tree species of the temperate zone such as oaks. the selfed offspring originate from one q. frainetto seed parent, the species with the lowest density in the area. according to streiff and her colleagues the effective pollen dispersal in oaks involves two components: local dispersal and long-distance transport. in our study, the portion of local pollen dispersal was quite variable among the seed parents. for instance, the local component of dispersal seems to be prevalent in q. frainetto and q. pubescens seed parents, the species with the lowest relative abundance within the whole forest reserve. in contrast, the amount of pollen originating from outside the study plot was much higher for the three q. petraea seed parents which is consistent with other results obtained via paternity analyses in different oak species. according to the forest records, q. petraea is the predominant species in the reserve and in the nearby forest which may explain the larger amount of pollen coming from outside the study plot for this species. by using paternity analyses we studied pollination patterns in a very complex situation involving all species which occur not only in the mixed stand but also in the nearby area on tens of kilometers. the importance of considering all species when studying hybridization was very well highlighted in france where hybrids were detected in considerable proportions in a mixed stand in the absence of parental species. our findings suggest the existence of reproductive barriers between oak species since the majority of the effective pollen originates from the same species. however, the hybridization rate was estimated at % in the seed generation. we found strong evidence for the production of first generation hybrids because both parents were phenotypically and genetically pure species. we observed f hybrids between three out of twelve possible species combinations, including reciprocal crossings, which suggest that hybridization at this site is not restricted to one species pair. this result is strongly influenced by our unbalanced sample of mother trees of different species. previously, in a two-oak-species stand streiff and her colleagues reported hybridization in out of offspring. however, this rate of interspecific gene flow was estimated considering only the physical appearance of oak trees. as in our study the hybridization rate varied markedly among maternal parents. those oaks showing high levels of hybridization were mostly surrounded by individuals belonging to other species. in such contact zones between species, hybridization may occur due to physical proximity. hybridization was not restricted to the formation of first-generation hybrids since % of the pollen donors were putative introgressed forms between the oak species. nearly 3/ of the total number of hybrids were later-generation hybrids. the predominance of backcrossed individuals in the detriment of first generation hybrids was reported for the same species complex in western europe. by comparing the overall proportion of hybrids in the adult generation with the proportion of hybrids in the offspring we found nearly 80% more hybrids in the seed generation. the decrease in the number of hybrids from offspring to adults as well as the maintenance of species identity may be the result of the selection against hybrids which acts in the ecological context of this forest. environmental selection against hybrids was invoked to explain the overestimation of hybridization rates when analyzing seeds. another explanation for this situation would be the stand history. the uneven-aged structure of the adult population might indicate that some of the trees were established when less than four species were found at this site. the higher proportion of hybrid acorns compared to hybrids adult trees may also be caused by past introgression events which have reduced the strength of the reproductive barriers over time. our results indicate that hybrids participate at reproduction both as pollen donor and receptor. similar situations were reported in a mixed stand with q. petraea, q. pubescens and intermediate individuals between both species in italy. moreover, hybrids between one pair of species prefer to mate with one of their parental species. only in two instances we found hybrids which mate with other species than the parental ones. these preferential matings indicate weaker reproductive barriers between the parental species and their hybrids as compared to other related species and hybrid combinations. the asymmetric gene flow from q. petraea to q. robur revealed by paternity analysis in our sample is consistent with the outcome of the bayesian analysis in the adult generation. a possible explanation is the contrasting relative abundance of the two species at bejan: q. petraea is very common in the study area while q. robur is frequent only in a small marginal part of the reserve. higher success rates of the combination q. robur × q. petraea were also observed in artificial crosses. our finding supports the model of asymmetrical hybridization between q. petraea and q. robur. moreover, the gene exchange between q. pubescens and q. petraea was mainly from q. pubescens towards q. petraea in our sample. the opposite direction was observed between the two species in a hybrid zone in italy. the ecological conditions and forest composition which differ between the two geographical areas may have influenced the direction of introgression. these discrepancies between hybridization patterns described throughout oaks' distribution range stress the considerable role played by ecological factors and the necessity for comparative studies. the level of hybridization may be higher for one individual or another because of an overlap in flowering time with most of the individuals of other species. almost no differences in flowering time between oak species have been reported in different mixed forests. for instance, large differences in flowering time between the single q. robur seed parent and the majority of q. robur adult trees and an overlap of the flowering time with most of the q. petraea individuals can be advocated as one of the causes to explain the high propensity for hybridization of this seed tree. differences in flowering time may lead to a preponderance of pollen coming from q. petraea relative to the pollen coming from its conspecific individuals. phenological observations at this site are needed to validate this hypothesis. furthermore, the present results are based on only one flowering season and the patterns of pollen movement may certainly vary from year to year. hence, the physical barrier represented by q. cerris trees might have been also hindered the pollen coming from the eastern part of the study plot, where q. robur is predominant. the low amount of genetic differentiation among the effective pollen pools for seed parents which belong to the same species may be due to the limited number of sampled maternal parents. genetic heterogeneity of effective pollen pools was reported in other studies based on larger sample sizes. genetic heterogeneity of effective pollen pools may be caused by uneven reproductive contribution of local pollen donors. differences in flowering phenology have also been mentioned as a possible cause. however, the differentiation of effective pollen pools for seed parents which belong to different species was highly significant, except for the effective pollen pool received by the single q. robur seed parent which was very similar to the effective pollen pool of q. petraea trees. the low amount of genetic differentiation between the effective pollen pool received by the q. robur seed parent and the three other q. petraea neighboring trees is consistent with the result of the paternity analysis that indicated a high rate of hybridization with q. petraea trees. in contrast to the study of streiff and her colleagues a significant differentiation between the effective pollen pools from outside and inside the study plot was detected. this genetic heterogeneity may be caused by the small number of pollen donors identified within the study plot that were not representative for the genetic composition of the whole oak forest. our data did not follow an isolation by distance pattern like in other tree species. the high number of oak species which are present among the sampled maternal trees correlated with their preference for conspecific pollen rather than heterospecific pollen may explain this pattern.  by using paternity analysis and a set of six highly polymorphic microsatellite loci we have traced patterns of contemporary gene flow within and between oak species in a natural reserve. our results confirm that introgressive hybridization between different pairs of species is quite common and play a significant evolutionary role within the white oak species complex. selection against hybrids over years might contribute to the maintenance of species identity despite considerable amounts of interspecific gene flow. however, the differences between the proportion of hybrid acorns and the proportion of hybrid adults may also be explained by the stand history or past introgression. further phenological investigations and a larger sample size are required for a better understanding of the hybridization events taking place in this four-oak-species stand. authors' contributions alc conceived the idea, collected the seed samples, did the laboratory work, performed the paternity analysis, and wrote the majority of the text. og helped to analyze the data and to draft the manuscript. rf participated in the design and coordination of the study, data interpretation and made substantial contributions to writing this paper. 
the minimum amount of effective pollen originating from outside the plot varied markedly among the seed parents, ranging from to. only two out of eight seed parents, belonging to q. petraea and q. robur, showed a high propensity to hybridize with q. pubescens and q. petraea, respectively. however, the isolation was not complete since we found strong evidence for both first-generation and later-generation hybrids. the majority of the successful matings occurred between conspecific individuals indicating the existence of reproductive barriers between oak species. first-generation hybrids were observed only between three out of six possible species combinations. to investigate the contemporary pattern of hybridization we sampled and genotyped offspring from a natural mixed forest comprising four species of the european white oak complex: quercus robur, q. petraea, q. pubescens, and q. frainetto. the observation of first and later-generation hybrids in higher frequency in acorns than in adults might be explained by selection against hybrid genotypes, the history of this uneven-aged forest or past introgression between species. significant structure of the effective pollen pools was detected in our sample. however, no support was found for the isolation by distance hypothesis. the proportion of hybrids was much higher in the seed generation when compared to the adult tree generation. a total of offspring were assigned unambiguously to one of the pollen donors within the study plot. hybrids between one pair of species preferred to mate with one of their parental species.
 few studies address the issue of hybridization in a more than two-species context. the species-rich quercus complex is one of the systems which can offer such an opportunity. to investigate the contemporary pattern of hybridization we sampled and genotyped offspring from a natural mixed forest comprising four species of the european white oak complex: quercus robur, q. petraea, q. pubescens, and q. frainetto.  a total of offspring were assigned unambiguously to one of the pollen donors within the study plot. the minimum amount of effective pollen originating from outside the plot varied markedly among the seed parents, ranging from to. the majority of the successful matings occurred between conspecific individuals indicating the existence of reproductive barriers between oak species. however, the isolation was not complete since we found strong evidence for both first-generation and later-generation hybrids. only two out of eight seed parents, belonging to q. petraea and q. robur, showed a high propensity to hybridize with q. pubescens and q. petraea, respectively. significant structure of the effective pollen pools was detected in our sample. however, no support was found for the isolation by distance hypothesis. the proportion of hybrids was much higher in the seed generation when compared to the adult tree generation.  first-generation hybrids were observed only between three out of six possible species combinations. hybrids between one pair of species preferred to mate with one of their parental species. the observation of first and later-generation hybrids in higher frequency in acorns than in adults might be explained by selection against hybrid genotypes, the history of this uneven-aged forest or past introgression between species. 