Gene flow among wild and cultivated common walnut (Juglans regia) trees in the Qinling Mountains revealed by microsatellite markers

respectively. Pollination of mother trees was not by nearest neighbors, and a paternity of 60.7% of offspring evaluated could not be assigned to a local, sampled tree. Pollen flow from cultivated trees to wild trees was infrequent (5.4%), and selfing rates ranged from zero to 25.0%. Male parents were located from 0 to 1005 m from their female partners, with an average pollination distance of 285.1 m. These results are discussed in light of the cultivated species’ diversity, outlining the frequent spontaneous genetic contributions from the wild to the cultivated compartment. In addition, the pollen flow parameters provide useful information about the dynamics of pollen movement within J. regia populations.


Introduction
Mating patterns, including pollen dispersal, gene flow, and the movement of genes among lineages, play an important role in the evolution of organisms by shuffling the genetic diversity and structure within species (Jones et al. 2008). Genetic invasion between wild species and cultivars by pollen dispersal, seed dispersal, and hybridization has played an important role in plant evolutionary history, especially the domestication of trees as crops (Miller and Gross 2011;Delplancke et al. 2013). Historically, wild species have been used in breeding programs as a source of morphological and genetic variation for plant improvement. In places where crops are grown near their wild ancestors, the reverse movement of genes, from cultivated to wild plants, has important evolutionary consequences for local conspecifics and their relatives (McGranahan et al. 1988; Barbour et al. 2002;Delplancke et al. 2012;Cornille et al. 2013). Because cultivated genotypes can be highly competitive, they Abstract Patterns of gene flow and gene introgression can be used to assess the risk of genetic pollution of wild forest trees from widespread cultivated trees. A comprehensive understanding of the genetic relationships and levels of gene flow among wild and cultivated common walnut (Juglans regia) has become an urgent issue. Using twelve microsatellite markers, we investigated the genetic diversity and gene flow between cultivated and wild trees of J. regia in the Qinling Mountains, China. A high level of genetic variation was detected in both cultivated and wild trees. The mean number of alleles per locus was 17.5. Observed heterozygosity (H O ) and expected heterozygosity (H E ) were 0.777 and 0.800, can contribute to the local exclusion of wild conspecifics (Trucco et al. 2009). In other cases, cultivated plant species have been selected for non-competitiveness or are less robust because of genetic bottlenecks caused by long-term cultivation (Miller and Gross 2011;Delplancke et al. 2012), but these phenotypes and domestication bottlenecks are often much less pronounced in perennial tree crops.
Hybridization between cultivated and wild trees is facilitated by the frequent absence of an intrinsic reproductive barrier (Gepts and Papa 2003;Cornille et al. 2013). Cultivated to wild gene flow can reduce biodiversity, cause a decline in reproductive fitness, and lead to loss of genetic distinctiveness of native populations (Wolf et al. 2001;Sánchez-de León and Johnson-Maynard 2009). Knowledge of pollen dispersal patterns within forest landscapes is important for managing the integrity of native gene pools by limiting the potential for genetic pollution from non-native plantings.
Walnuts (most Juglans species) belong to a widespread genus of wind-pollinated heterodichogamous species (Bai et al. 2006(Bai et al. , 2007Kimura et al. 2012;Zhao et al. 2018). The common walnut (also called the English or Persian walnut, J. regia L.) is a monoecious wind-pollinated deciduous tree native to Eurasia but now grown in many countries with a temperate climate (Manning 1978;Gunn et al. 2010;Aradhya et al. 2017). It has been cultivated for its timber and nuts for at least 6800 years (Beer et al. 2008;Feng et al. 2018). In many parts of Asia it is difficult to know if autochthonous trees are wild or feral (Pollegioni et al. 2015;Aradhya et al. 2017;Bernard et al. 2018) as there are no obvious phenotypic indicators of domestication. Nuts from apparently wild trees are harvested by villagers who transport them to markets where they may be sold alongside locally grown improved varieties. Because common walnut is long-lived and anemophilous, gene flow between wild trees and modern cultivated genotypes or their seedlings is likely, even in areas that may be centers of diversity.
We investigated gene flow between (apparently) wild and cultivated common walnut trees in the Qinling Mountains of China, a region believed to contain an ancient population of J. regia (Cheng et al. 2013;Guo et al. 2013). Our goals were to: (1) Characterize the genetic diversity of common walnuts (J. regia) in the Qinling Mountains; (2) Investigate the relative genetic contribution of cultivated J. regia to the seeds of wild trees in the Qinling Mountains; and, (3) Determine the pollen dispersal distance at research locations.

Sampling and plant materials
Common walnut is valued for both nuts and timber. To characterize the genetic diversity of cultivated common walnuts, we sampled seeds from trees planted by farmers in 16 locations around the Qinling Mountains, Shaanxi province, China (Table 1, Fig. 1). For clarity, these cultivated trees are described as "cultivars" even though they were not grafted but instead chance seedlings or seedlings of locally desirable genotypes. At sites with cultivated trees, a sample of genotypes were collected that were traditionally cultivated by local farmers who had no access to modern, commercial cultivars. To investigate the relative genetic contribution of cultivated J. regia to the seeds of wild trees, we sampled 275 seeds from 19 mother trees at a single location (Table 1; Fig. 1). Both the seeds and the mother trees were genotyped at 12 loci and a parentage analysis used to determine the percent of offspring derived from self-pollinations, pollinations from nearby mother trees, pollinations from trees outside the study area that were not genotyped, and pollinations from the closest cultivated population (HH and II; Fig. 1). To determine the pollen dispersal distance (Table 1, Fig. 1), the sample locations were mapped using GIS software ArcGIS (version 10.0; ESRI 2010).

SSR genotyping
DNA was extracted from leaf and seed samples using the methods of Doyle and Doyle (1987) and Zhao and Woeste (2011). Twelve microsatellite loci (Woeste et al. 2002;Victory et al. 2006;Ahmed et al. 2012;Pollegioni et al. 2015) ( Table 2) were used for genotyping. Each forward primer was marked with fluorescent dye (FAM, TAMRM, HEX and ROX). The PCR products were sent to Shagon Biotech Life Science Products and Services for analysis. The genetic analysis process was similar to the methods in Feng et al. 2018.

Microsatellite data analysis
Genemapper software (Applied Biosystem, Foster City, CA, USA) was used to score all genotype data. After removing 18 samples of seeds of cultivated individuals (seeds) and 36 samples of seeds of wild individuals (seeds) that had genotypes that did not match their presumed maternal parent based on CERVUS v3.0 (Kalinowski et al. 2007), a total of 317 individuals were retained for analysis. The exclusion probabilities for one parent and two parents was 80%. We tested for evidence of null alleles using MICRO-CHECKER (Van Oosterhout et al. 2004).
Deviations per locus and per mother tree from Hardy-Weinberg equilibrium (HWE) were tested using GENEPOP (Raymond and Rousset 1995) and Arlequin 3.5 (Excoffier and Lischer 2010). CERVUS v3.0 software (Kalinowski et al. 2007) and POPGEN version 1.32 (Yeh and Boyle 1996) were used to calculate the genetic diversity parameters, including percent polymorphic loci (95% criterion), the total number of alleles (Na), effective number of alleles (Ne), allelic richness (Ar), gene diversity (H S ), and overall gene diversity (H T ). The observed (H O ) and expected (H E ) heterozygosities were calculated based on the allele frequencies using POPGENE (Yeh and Boyle 1996).
Genetic differentiation between groups (F ST ) and linkage disequilibrium (LD) between loci were tested using the program GENEPOP v4.0 (Raymond and Rousset 1995). The significance of difference between genetic groups, measured as F ST , was determined by permutation tests (10,000) using Arlequin v3.5 (Excoffier and Lischer 2010). R ST was calculated using SPAGEDI (Hardy and Vekemans 2002). We performed a two-level (among trees and with trees) analysis of molecular variance (AMOVA) based on 239 seed from 19 wild mother trees using ARLEQUIN v3.5 (Excoffier and Lischer 2010).
The overall genetic variation among cultivated and wild trees was represented using principal coordinate analysis (PCoA) using GenALEx v6.5 (Peakall and Smouse 2012). To evaluate the admixture and genetic populations within wild trees only, a genetic structure analysis of 239 wild seeds was performed using a model-based Bayesian clustering of individuals, as implemented in STRU CTU RE v2.3.4 (Pritchard et al. 2000). We also add the cultivated seeds to the STRU CTU RE analysis to check the genetic mixture. This software uses a Markov Chain Monte Carlo (MCMC) framework in which the algorithm explores a parameter space considering individual admixture proportion, locus-specific ancestries, population allele frequencies, and the expected admixture of the data set, assuming a user-defined K number of groups. STRU CTU RE was run using 100,000 burn-in MCMC iterations, with run lengths of 1,000,000 iterations, and ten replicates per run for K = 1-15 clusters, with admixture (Pritchard et al. 2000). STRU CTU RE HARVESTER was used to calculate the optimal value of K (Earl and Vonholdt. 2012) using the deltaK criterion (Evanno et al. 2005), and the inferred clusters were drawn as colored box-plots using program DISTRUCT (Rosenberg 2004).

Maternity and paternity analysis
A total of 275 seeds from 19 adult wild walnut trees were analysed for maternity using CERVUS v3.0 (Kalinowski et al. 2007), a software based on maximum-likelihood methods. The simulation parameters were: 10,000 cycles, 19 candidate mother trees, 19/275 as the proportion of candidate mothers sampled, 0.975 as the proportion of loci typed, 0.01 as the proportion of mistyped, 0.01 error rate in likelihood calculations, 0.80 for the relaxed confidence level and 0.95 for the strict confidence level. After verification of maternal parentage for each offspring, the number of candidate male parents specified in the simulation was the approximate size of the  Table 1). The number of the different scales on the map. QH = Qinghai Province, GS = Gansu Province, SX = Shaanxi Province, SXI = Shanxi Province, HN = Henan Province, SC = Sichuan Province, HB = Hubei Province, CS = Chongqing Direct-controlled Municipality male gene pool for the cultivated and wild trees in this study, which included both sampled and unsampled male parents. The total of number of identified male parents was 371 and the proportion of male parents sampled were adjusted to 0.80, 0.05 as the proportion of mistyped, 0.05 error rate in likelihood calculations, 0.65 for the relaxed confidence level and 0.95 for the strict confidence level.

Mating pattern and pollen dispersal
The relatedness among all 19 adult wild trees was calculated using software RELATEDNESS v5.0.5 (Goodnight and Queller 1999). Observed distances of pollen migration were used to construct a dispersal curve to describe pollen movement within the population. Measurements consisted of the distance between the mother and the putative father when the most likely father was assigned. In addition, the distances between each mother tree and all other candidate fathers were also measured to determine whether the observed inter-mate distance was influenced by the spatial arrangement of the adult trees. Gene flow (Nm) was estimated using Arlequin v3.11 (Excoffier and Lischer 2010) and MIGRATE v3.6 (Beerli 2006).

Genetic variation and structure of wild walnut seeds and cultivars as revealed by microsatellite data
Estimates of genetic diversity varied among microsatellite loci and among the wild offspring cohorts and adult (wild) trees [N = 18 after removal of wild parent Tree-N ( Fig. 1) due to poor DNA quality]. Among wild samples (N = 239), alleles per locus (Na) ranged from 9 to 34, with a total of 192 alleles across the 12 loci. The mean number of alleles per locus was 16.0. The mean polymorphic information content (PIC) was 0.774 (Table 2). Allelic richness (Ar) estimated among all trees ranged from 4.00 to 9.06 with an average of 6.85. Observed heterozygosity (H O ) among loci was 0.803, with a range of 0.548-0.927. Expected heterozygosity (H E ) among loci was 0.796, with a range of 0.563-0.957 (Table 2). Among samples from cultivated trees (78 seeds from 16 locations), the range of alleles per locus (Na) was from 7 to 29, with a total of 187 alleles across the 12 loci. The mean number of alleles per locus was 15.58. The mean polymorphic information content (PIC) was 0.779 (Table 3).

Genetic structure
Bayesian clustering of 239 progenies from 18 wild mother trees identified three populations as most likely (Fig. 4). These populations are distinguishable despite relatively high gene flow among the wild mother trees (Nm = 1.534). Visual examination of the spatial arrangement of the populations within the site from which the wild seeds were collected   (Fig. 1). The wild populations were genetically distinct from the seeds from cultivated trees (Fig. 3), which all belonged to a common population.

Maternity and paternity analysis
Paternity analysis of 371 seeds revealed that the male parent of 94 (25.3%) samples could be identified based on the genotypes of the seeds and their presumed female parent. Of these, 36 samples (15.1%) were seeds from cultivated trees, 61 (25.5%) were from wild adult mother trees with 80% confidence (Table 4). Among wild trees, mean distance between successful males and their female mates was 285.1 m, ranging from 61.2 to 1005.5 m, not including self-pollinations. Most successful pollinizers (75%) were from 200 to 600 m from their mate, and only three trees were more than 1250 m from their mate (Fig. 5). The rate of gene flow (0.425) between wild and cultivated J. regia was estimated using F ST, which was 0.088, ranged from 0.046 (GG-LL) and from 0.147 (FF-NN), which was lower than the average F ST among 18 progenies of wild trees (F ST = 0.094), ranged from 0.053 (J-K) and from 0.137 (C-I). The genetic differentiation (F ST ) was 0.042 between wild and cultivated groups; AMOVA revealed that most of the genetic variation (96%) was within wild or cultivated trees (Table 5). Paternity analysis showed that none of the wild mother trees was father of a seed from a cultivated tree, although coalescent-based estimations of gene flow (Nm = 1.48 among wild mother trees; Nm = 2.58 among cultivated trees, Nm = 1.76 between wild and cultivated trees) indicated that cultivated J. regia was expected to contributed about 13 migrants (5.4%) to the seeds of the wild progenies (Fig. 6). showing the locations of eighteen wild mother trees. Small circles and pie graphs indicate the profiles of genotyped progeny arrays from each mother tree for seed lot collection from 2011. Large letters around the pie indicate the wild mother tree's name. Dark gray, white, light gray, and black wedges represent selfing, male parent was a cultivated plant, two parents within study site, and unassigned seeds, respectively (details see Table 5)

Mating system and pollen dispersal
Of 239 progeny samples, 36 (15.1%) were determined to have arisen from self-pollination. Among wild mother trees, the selfing rate varied between zero to 16.7% (Table 6; Fig. 6). The percentage of seeds with two wild parents within the study site varied among mother trees from zero to 38.5% (Table 6). In addition, 145 progenies (60.7%) were not assigned pollen parents (paternity analysis did not match any sampled wild or cultivated genotype). These were most likely the result of immigrant pollen from outside the study area or were sired by un-sampled trees within the study site. Contamination rates for individual wild mother trees ranged from 50.0% to 91.7% (Table 6).

Performance of microsatellite markers in genetic variation and genetic structure of walnut
Microsatellite markers provided robust and highly informative signals based on high levels of polymorphism (Woeste et al. 2002;Dangl et al. 2005;Victory et al. 2006;Pollegioni et al. 2015). In general, it is expected that the genetic diversity of wild trees would be higher than cultivated trees. Our investigation revealed genetic variation in progenies of wild and cultivated J. regia in the Qinling Mountains were similar, although there were some important differences, and wild and cultivated trees belonged to distinct genetic groups (Fig. 3). Cultivated and wild trees were similar in terms of number of alleles (Na), expected heterozygosity (H E ), and polymorphic information content (PIC) (Tables 2 and 3). In contrast, cultivated trees and wild trees differed in terms of effective number of alleles (Ne) (1.93 versus 3.41) and observed heterozygosity (H O ) (0.672 versus 0.803) (Tables 2  and 3). These results indicate that the cultivated trees and wild trees differed considerably in two important measures of diversity, possibly because of human selection, although it is likely that most cultivated trees in the Qinling Mountains are only one or two generations removed from a wild population.

Contributions of wild walnuts to cultivated trees
These results are in accordance with studies of other plants that revealed substantial genetic exchange between wild and cultivated varieties of the same species. For example, patterns with a similar trend were observed for maize (Matsuoka et al. 2002), apples (Coart et al. 2006;Cornille et al. 2012), cultivated olive trees (Breton et al. 2008), and grape cultivars (Myles et al. 2011). Wild walnut trees have long been used as a source of genetic novelty in breeding programs (Pollegioni et al. 2015). On the other hand, these results raised a question of whether wild-to-cultivated exchanges remained exclusively spontaneous or whether traditional practices could have facilitated the fixation of introgressed wild genes into cultivated walnut trees. For Juglans, several authors reported the direct use of wild species as rootstocks (McGranahan and Leslie 2009). Additionally, germplasm resources are rich and the varieties numerous in J. regia which may be caused by a long cultivation history, and many are biologically and economically important because of their high-quality timer and nutritious nuts. Traditional practices in walnuts breeding remain largely unknown s (Rogers 2004;Beer et al. 2008).
Nuclear microsatellites (SSRs) showed that genetic exchanges between cultivated J. regia and wild J. regia were directional (from cultivar to wild), and outlined substantial cultivated-to-wild gene flow (Nm = 0.4252), confirming that J. regia genes could spontaneously be introgressed into their wild resources. And this genetic transfer can have significant evolutionary consequences, especially if the inserted transgene is adaptive under natural conditions (Felber et al. 2007;Sun et al. 2019). For example, traits such as enhanced fertility (Zhao and Woeste 2011), late flowering (Cornille et al. 2012), andantifungal (e.g., S. clavigignenti-juglandacearum, Ross-Davis et al. 2008) could favor the emergence of highly competitive phenotypes if transferred into wild species. Our results also highlighted the importance when inferring the origins of walnut domestication using molecular genetic data, because gene flow from wild walnuts can have a bias on distance-based ancestry inferences (Cornille et al. 2012;De Andrés et al. 2012). Finally, this study revealed that cultivated-to-wild gene flow occurred commonly among the walnuts located in the Qinling Mountains. There are no barriers preventing gene flow between wild and cultivated walnut trees at the current time in the study site. Further, these results suggest that transgenes could potentially introgress from cultivated into wild Persian walnut populations, demonstrating the  Table 6 The wild progeny samples mating system, out crossing rates, and pollen dispersal genotyped using 12 SSR loci Name need for detailed characterization of cultivated to wild gene flow within the Juglans species complex (Yuan et al. 2018;Dang et al. 2019).

Mating system, out crossing rates, and pollen dispersal
The mating system of walnuts is characterized by a protogynous-protandrous dimorphism which is a form of heterodichogamy (Bai et al. 2007;Kimura et al. 2012) intended to avoid self-pollination. A genetic model of this system predicts a stable polymorphism with both mating types in equal proportion in a randomly mating population (Härdling and Bergsten 2006). The proportion of progeny derived from self-pollinations in our study was 15.1%, less than three times as outside trees. This selfing rate was considerably higher than rates reported in other Juglans species (Bai et al. 2007) and in other wind-pollinated species (Kimura et al. 2012). The rate of self-pollination can vary over years and among populations and may be strongly affected by local factors (Kery et al. 2003). Bai et al. (2007) showed that self-mating was rare in a high-density, morph-ratiounbiased population of the heterodichogamous Manchurian walnut, Juglans mandshurica Maxim. Kimura et al. (2012) reported that the selfing rate was 13.8% of Japanese walnut (J. ailantifolia Carr.) in a 2.5 ha study site. The population was wild, low density, and morph-ratio-biased. In our study, selfing rate varied between zero to 25.0% with a mean of 15.1% (Fig. 6); the highest rates of selfing were observed in tree-M and tree-H; Fig. 6). Moreover, the un-assigned wild trees impact the selfing rate caused by the missing data in this study. The simplest explanation for variance in selfing rates among trees is difference in the availability of mating partners in the neighbourhood, but factors related to phenology likely play a role, especially within-plant bloom overlap. Therefore, a relatively lengthy intra-morph bloom overlap may have affected the extent of assortative mating in this study. The inbreeding coefficient for the seeds of the wild parents was 0.077. Averaged across all seeds, the mean distance between parents (d) was 285.1 m for wild-mother trees. Immigration rates were approximately 60.9%. These findings indicate that pollination by distant males may be relatively common in this population (Fig. 5). Kimura et al. (2012) showed that distance to pollinizers in a wild population of Juglans ailantifolia was 392 m and 233 m for protogynous and protandrous mother trees, respectively. Bai et al. (2007) showed that pollen dispersal in J. mandshurica was predominantly short-distance, and local density of reproductive trees affected the patterns of pollen dispersal in forest trees. Mating patterns were thus also affected by both density and tree size (Jolivet et al. 2012).
Our study revealed that in this population of J. regia in the Qinling Mountains, the longest pollination distance was 1.0 km. Large-scale studies are needed to disentangle relative influences of these factors on the mating system and pollination success (Jolivet et al. 2012;Davies et al. 2015).

Conclusions
Pollination of mother trees was outcrossed, often by neighbors, and pollinizers outside the immediate sampled area accounted for about 60.7% of the sampled seeds. In a few cases, pollination of wild trees was by nearby cultivated trees (5.4%). The selfing rate among the wild trees varied from zero to 25.0%. Pollinizers ranged from zero (self-pollinations) to 1005 m from their respective female partners, with an average distance of 285.1 m. This description of pollen flow provides useful information about the dynamics of pollen movement within wild J. regia populations.
Acknowledgements The authors wish to thank Hailong Xia, Qi Xu, Jia Yang, and Tao Zhou for sample collection. Mention of a trademark, proprietary product, or vendor does not constitute a guarantee or warranty of the product by the U.S. Department of Agriculture and does not imply its approval to the exclusion of other products or vendors that also may be suitable.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.