Abstract
Walnuts have substantial economic value and are of significant interest being a wild-cultivated species. The study has re-sequenced the entire genome of the wild walnut, aligning it with the walnut reference genome, to identify 2,021,717 single nucleotide polymorphisms (SNPs). These were used to examine the genetics of 130 wild walnut samples collected from three countries. Utilizing structural and principal component analysis, the walnut samples from Central Asia were classified into four populations: Ili ah in Xinjiang (I), Dushanbe region in Tajikistan (II), Sary-Chelek, Arslanbob in Kara-Alma regions of Kyrgyzstan (III), and Kok-Tundy region of Kyrgyzstan (IV). The 4 groups showed large differences in nucleotide diversity, population differentiation, and linkage disequilibrium decay, as well as gene flow among them. The present geographic distribution of these populations does not align with the genetic distribution pattern as the populations of Central Asian wild walnuts have experienced similar population dynamics in the past, i.e., the highest effective population size at ca. 6 Ma, two sharp population declines at 6 and 0.2 Ma, and convergence at ca. 0.2 Ma. The genetic distribution patterns are better explained by human activity, notably through archaeological findings of walnut use and the influence of the Silk Road, rather than by current geographic distributions.
1 Introduction
The walnut (Juglans regia L.) belongs to genus Juglans (Juglandaceae) and it is also known as the Persian walnut, and it is a wind-pollinated, monoecious, long-lived, perennial tree that has been domesticated and cultivated worldwide (, , ). The walnut tree is highly valued for its nutritional and medicinal values found in the nuts, as well as for the quality of its wood (; ). The natural habitats of walnut wood are confined to the mountainous regions of the Asian continent, stretching from Xinjiang Province in western China to the Caucasus (, ; ). The relationship between the ancient centers of diversity for walnuts and their domestication and subsequent spread by humans remains a subject of considerable debate (; ; ). Recent studies of plant domestication history have mostly utilized agronomy, genetics, and archaeobotany (; ). Quantitative trait locus localization, genome-wide association studies, and genome-wide resequencing studies have identified genes associated with the initial domestication and further diversification of crops (; ). Additionally, the integration of phylogeny, archaeology, and genomics provides a broader insight into the domestication process for various plant species (; ; ). The combination of whole genome resequencing data with fossils and archaeological evidence offers a comprehensive method to review the domestication patterns in cultivated plants (; ), as the case here studied.
Juglans regia is thought to be a Neoproterozoic relict of Tertiary forests in Eurasia, where fossil pollen evidence suggests that walnuts have been present since the Pleistocene () and that his record persists from the mid-Paleoproterozoic to the Neoproterozoic (). However, the evolutionary history of this species during the Late Tertiary climatic deterioration and Quaternary glaciation is replete with mass extinctions, range reductions, fragmentation, and bottlenecks (; ). In addition, the isolation of walnut groups from each other was accelerated by the gradual desertification of Central Asia during the Holocene (). Recent studies have proposed that J. regia originated as an ancient hybrid between American and Asian walnut lineages in the late Pliocene (). Although, man-made walnut is native to the mountains of southeastern Europe and western Central Asia (), recent research suggests that Xinjiang Province in China could also be the origin of J. regia (). The relationship between the ancient centers of diversity for walnuts and their domestication and subsequent spread by humans remains a subject of considerable debate (). Human activities have had a profound impact on the distribution and range expansion of walnuts throughout history (). There is a substantial historical evidence of humans’ close association with walnuts. Dried walnut seeds have been found in the Near East (6230–5790 yr. BP) (), Central Asia (5149 yr. BP) (), and along the Yellow River Basin in northeastern China (7300 yr. BP) (). The consumption of walnuts, valued for their edibility and religious significance, dates back to Persia (7,000 BCE) (), and was facilitated by trade networks such as the Silk Road (). Despite the long history of use, the domestication of walnuts is relatively recent, with the improvement of asexual varieties of Persian walnuts since the middle of the 20th century (; ). Plant cultivation and domestication is a spatially and temporarily dynamic multi-stage process, and taking populations of wild plants to cultivated forms requires human intervention to ensure survival (; ; ). Generally, some reproductive methods or fruit characteristics of perennial plants are changed in this process, but this is not evident in walnuts, due to the fact that the reproductive practices of walnuts under artificial conditions do not differ from those of wild walnuts (; ; ; ; ).
The advancement of molecular markers, including Randomly Amplified Polymorphic DNA (RAPD), simple sequence repeats (SSR), amplified fragment length polymorphism (AFLP), expressed sequence marker SSR (EST-SSR), and genomic SSR markers (gSSR), has significantly enhanced the study of plant genetic diversity and structure (; ; ; ; ). The study of population genetic structure is essential to analyze the adaptive evolution and genetic relationships of walnut populations across various regions (). Single nucleotide polymorphisms (SNPs), defined as any single base substitution/indel in the genome of an individual (), represent next-generation sequencing technologies. SNPs have gained popularity due to the reduced cost of next-generation sequencing, their abundance in the genome, and their ability to express a rich amount of genetic information (; ). The walnut genome sequenced, SNP markers can be utilized for rapid and highly automated genotyping, making them ideal for genetic diversity studies and genomic selection analyses.
Central Asia, a region encompassing both the distribution area of wild walnuts and a historically crucial segment of the ancient Silk Road, is of great significance for the study of wild walnut genetic diversity which is important for the conservation and utilization of the walnut species. In this study, we sequenced 130 wild walnuts samples from Central Asia, including Xinjiang China, the Wild Fruit Forest Reserve of Kyrgyzstan, and the Dushanbe region of Tajikistan, using whole genome sequencing technology and walnut reference genomes to elucidate the genetic relationships and genetic diversity in these regions. Further, our aim was to elucidate the genetic relationships and genetic diversity within these regions and to describe the spatial distribution of walnut genetics in Central Asia, considering both geo-historical variability and the impact of human activities throughout history.
2 Materials and methods
2.1 Plant materials and DNA extraction
Samples were collected from large natural walnut forests in the mountains or along forest roads. A minimum of 20 meters was maintained between sampled trees, and a total of 130 fresh young leaf samples from healthy J. regia (Table 1) were collected. Of the 130 samples, 27 came from Yili, Xinjiang, China (Group Y). 83 samples were gathered from four key wild fruit forest reserves in Kyrgyzstan: 21 samples were collected from the first area, Sary-Cheek (Group J); 21 samples were collected from the second area, Arslanbob (Group K); 26 samples were collected from the Kok-Tundy Area C (Group L), and the last 15 samples came from the Kara-Alma area (Group M). Additionally, 20 samples were obtained from Dushanbe, Uzbekistan (Group D). The fresh leaves were dried in silica gel and stored for future analysis. The geographical distribution of all sampling points is plotted by ArcGIS v10.2 (ESRI Inc., Redlands, California, USA) (Figure 1).
Table 1
| Sampling Site | Country | ID | N | Latitude | Longitude | Elevation (m) |
|---|---|---|---|---|---|---|
| Yili, Xinjiang | China | Y | 27 | 43.217701 | 82.152833 | 1533 |
| Dushanbe | Tajikistan | D | 20 | 38.837793 | 68.959988 | 1752 |
| Sary-Chelek | Kyrgyzstan | J | 21 | 41.85481 | 71.964071 | 1778 |
| Arslanbob | Kyrgyzstan | K | 21 | 41.339877 | 72.889697 | 1652 |
| Kok-Tundy | Kyrgyzstan | L | 26 | 40.853421 | 73.66128 | 1319 |
| Kara-Alma | Kyrgyzstan | M | 15 | 41.173196 | 73.311407 | 1264 |
The 130 walnut samples from three regions of Central Asia.
Figure 1
DNA extraction was performed using the improved CTAB method (). Since high-quality gDNA is required in the RAD genotyping technique, we evaluated the concentration and purity of the extracted DNA. These assessments were determined by the ratios of absorbance at 260/280 nm and at 260/230 nm, using a NanoDrop ND-2000 spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA). This crucial procedure ensured high-quality samples and comparable DNA concentrations for subsequent analysis.
2.2 Sequencing and sequencing quality control
For library preparation, DNA was isolated from fresh leaves with 2% cetyltrimethylammonium bromide (CTAB) according to a previously published protocol. A TruSeq library was prepared using the KAPA Hyper Prep Kit (Illumina® platforms, KAPA BIOSYSTEMS, Boston, MA, USA. Cat No. KK8504) following the manufacturers’ manual. All the libraries were sequenced by the HiSeq X® platform (Illumina®, San Diego, CA, USA) according to the manufacturer’s instructions. The genome samples were cleaved using the restriction endonuclease EcoRI, and then each one was physically fragmented to construct DNA libraries, and we selected a library of 200–400 bp fragments for insertion. The raw number is obtained by reading 150 bp from the opposite end. The HiSeq X-Ten platform (Illumina, Inc., San Diego, CA, United States) was used for adaptor ligation and DNA cluster preparation.
The original sequencing data contained joint information, low-quality bases, and undetected bases, which would greatly interfere with subsequent information analysis. To obtain clean reads and prevent further interference, the following filters were used: (1) remove reads containing splice sequences. After the first round of quality control using FastQC (version 0.11.5), DNA sequencing reads were trimmed and filtered using Trimmomatic (version 0.36.5). (2) Remove pairs of paired reads whose N content in a single-end sequencing read was more than 10% of the length (3) remove pair of paired reads when a single-end sequencing read contained more than 50% low-quality (<= 5) bases. Perform the above operations based on FastQC to ensure the quality of the trimmed data.
2.3 SNP calling and filtering
We used Burrows-Wheeler alignment (BWA) V0.7.17 software () (parameters: mem -t 4 -k 32 -M) to map clean reads to the walnut reference genome (Jregia.genome_v1.0) (). Use the SAMtools48 (version: 0.1.18) software to convert mapping results to BAM format and filter unmapped and non-unique readings. SAMtools and Picard are used for sorting and removing duplicates, respectively. We removed duplicates of the comparison results with the SAMTOOL () software (Parameters: rmdup), and duplicate reads were filtered using the Picard package (picard.sourceforge.net, version: 1.87).To ensure high quality data, we utilized SAMTOOL () to perform strict filtering of SNP calls based on the following criteria: (1) a minimum SNP support number (coverage depth) of more than 3, (2) a minimum allele frequency of more than 0.05, (3) a deletion rate of less than 10%. Variant identification, programed to detect SNPs, was conducted with the Unified Genotyper application of the Genome Analysis Toolkit (). GATK was used to construct the general variant call file and merge the gVCF files to form a single variant call file. In order to obtain high-quality SNPs, we first filtered the merged VCF data using the GATK hard filter. Using the ANNOVAR package for genome-based annotation, SNPs were categorized into exonic regions (overlapping coding exons), splice sites (within 2 bp of splice junctions), 5′UTRs and 3′UTRs, intronic regions (overlapping introns), upstream and downstream regions (upstream or downstream of the transcriptional start site by 1 kb), and intergenic regions. SNPs in coding exons were further categorized into synonymous SNPs (which do not cause amino acid changes) and non-synonymous SNPs (which cause amino acid changes; mutations that cause stop-gain and stop-loss are also included in this group), and indels in exonic regions were classified according to the presence or absence of frameshift (3 bp insertion or deletion) mutations.
2.4 Structure analysis and phylogenetic tree linkage disequilibrium analysis
To eliminate the potential effect of physical linkage between variants, we selected high-quality SNPs with no more than 20% missing data, and based on this, we refined the loci to ensure that no two loci were located within the same 2000 bp region. To investigate genetic structure among the walnut populations, we conducted a principal component analysis (PCA) using PLINK v1.09 () The principal components were also plotted against each other using R3.4 to visualize patterns of genetic variation, and the first three eigenvectors were plotted in 3D space. Additionally, we utilized a hierarchical Bayesian modeling using the STRUCTURE software; this involved a burn-in of 10,000 iterations and running 100,000 iterations with 10 replicates. Subsequently, we clustered the individuals according to K = 1–8 using STRUCTURE v. 2.3.4 (); this process utilized an admixture model with associated allele frequencies. The optimal value of K was determined using STRUCTURE HARVESTER v.0.6.94 () according to the Delta K method of Evanno (), Ln (D|K), and the final posterior probability of K (). IQTREE (version 1.6.9) () was used to perform phylogenetic analyses to construct maximum likelihood (ML)-based phylogenetic trees based on the GTR + F + R5 model, and 1000 fast bootstrap replicates were performed to determine branching confidence values. After evaluating 130 DNA models, ModelFinder was used to estimate the best-fit model, and finally an ML phylogenetic tree was constructed based on intergenic region SNPs. We used PHYLIP software v 3.696 to produce a cladogram based on the genetic distance matrix using the p-distance formula (). The phylogenetic tree mapping was done by Figtree using the Neighborhood Joining (NJ) method.
2.5 LD analysis and genetic diversity analysis
We calculated the pair-wise LD as r2 for each population and within chromosomes using Plink v1.09 following the formula proposed by Hill and Robertson (). The LD attenuation calculation was based on the r2 between two SNPs and the distance between two SNPs, and the LD decay was similarly determined based on the r2 between two SNPs and the distance between the two SNPs. The r2 statistics was obtained using the PLINK software (version 1.07) (). For each SNP pair, bins of 100 kb were created based on pairwise physical distance. The nucleotide diversity (θπ) of the six groups and genetic distance (FST) were calculated in stepping windows 40 kb in size by VCFtools v0.1.13 (). The horizontal axis represents the position of the genome, and the vertical axis represents the θπ value. The population θπ value reflects the base diversity of the population genome.
2.6 Estimated effective population size
To reconstruct the historical demography of the Asian walnuts, we used PSMC (). We filtered sequencing data by the following conditions: 1) a mean genome coverage of 18, 2) a per-site filter of 10 reads, and 3) no more than 25% of missing data were followed (). The parameters in PSMC were set to a mass adjustment of 50, a minimum localization mass adjustment of 20, a minimum depth adjustment of one-third of the average depth genome coverage, and a maximum depth adjustment of 2 times. In line with existing literature, the time from birth to offspring was set at 20 years (). The annual mutation rate (u) was estimated to be 5.6×10–10 nucleotide substitutions per year using Ka. This converts to 1.12×10–8 per generation (20 years).
3 Results
3.1 Sequencing and SNP variation
After filtering the sequencing data, we obtained high quality clean data, including sequencing data yield, sequencing error rate, Q20 content, Q30 content, GC content and other data. A total of 130 samples were sequenced in this batch, yielding a high sequence quality and volume of data (611.42 Gb). The lowest raw data output of a single sample was 3.46 G, and the highest raw data output was 6.38 G, while the average Effective Rate(%) was 99.83. The mean value of the proportion of bases with mass values >= 20 (Q20) was 97.18%, and the average of the proportion of bases with mass values >= 30 (Q30) was 92.38%. The average GC Content(%) was 37.19. The lowest Effective Rate of a single sample was 97.80, and the highest was 99.95. The lowest Q20 (%) for individual samples was 95.07, and the highest was 97.64. The lowest Q30 (%) for individual samples was 88.34, highest was 93.38. The lowest GC Content (%) for individual samples was 36.39, and the highest was 38.98. GC Distribution is normal (Supplementary Table 1). The average depth for the 130 samples was 7.99, with an average Coverage of 1X of 94.53% and an average Coverage of 4X of 78.81%. Over 90% of sequenced reads from the walnuts aligned with the reference genome (Supplementary Table 2). The sample comparison rate reflects the similarity between the sample sequencing data and the reference genome, and the coverage depth and coverage can directly reflect the homogeneity of the sequencing data and the homology with the reference sequence. The results show that their similarity to the reference genome meets the requirements for resequencing analysis.
Our SNP detection analysis revealed a total of 2021717 SNP variant loci. Of these, 73090 (3.62%) SNPs in Upstream, 84110 SNPs in Downstream (4.16%), and 1492191 SNPs in Intergenic (73.81%). Additionally, 5304 (0.26%) SNPs were found in 1 Kb region upstream of the gene and 1 Kb downstream of the other gene (Upstream/Downstream). There were 88,425 (4.37%) variants located in the exonic region, of which 1,226 (0.06%) were in Stop gain, 198 (0.01%) in Stop loss, 40,447 (2%) in Synonymous, and 46,554 (2.3%) in Non-synonymous. There were 278,186 SNPs located in Intronic, 367 (0.02%) in Splicing, and 149,191 (73.81%) in Intergenic (Figure 2A; Supplementary Table 3). The distribution of SNP loci on the chromosome is shown in Figure 3A, where the denser distribution of loci is displayed with darker shades, and less dense regions are shown with lighter shades (Figure 2B). ANNOVAR functionally annotates genetic variants detected in multiple genomes, enabling subsequent correlation analyses based on the availability of data on the chromosome where the variant is located, the start site, the stop site, the reference nucleotide, and the variant nucleotide.
Figure 2
Figure 3
3.2 Phylogenetic relationships and population structure of walnut
The STRUCTURE analysis showed that at K = 2 the whole panel could be divided into two groups: the Ili walnut ditch area formed one, and the remaining areas clustered together as the other group (groups D, K, J, M, L). When K = 3, the Dushanbe walnut group emerged as one group (group D), while the walnuts of the Kyrgyzstan area were still the other (groups L, M, J, K, groups). When K=4, the wild walnuts of the Kyrgyzstan region (groups L, M, J, K) were separated from the wild walnuts of the Kok-Tundy region. When K=5, the distinction between the three remaining Kirghizstan regions is not obvious. (Figure 3A).
The PCA corroborated our findings regarding the population structure of wild walnuts, and based on the results of STRUCTURE, it showed a clear structure separation of the panel. We utilized the PCA to analyze the relationship between groups in different places as a complementary method to visualize their clusters. This revealed four well-separated clusters among 130 samples of wild walnuts. The largest cluster contains 51 J. regia accessions, and the second cluster mainly includes samples of Ili wild walnuts. The third contains samples of wild walnuts from the Kok-Tundy region of Kyrgyzstan, and the fourth major contains samples of wild walnuts from the Dushanbe region. The findings of both STRUCTURE and PCA revealed that the wild walnut population in Xinjiang was strikingly distinct and unique from other populations. Furthermore, that data indicated obvious gene flow between the walnuts in the three regions of M, J, and K (Figure 3C.).
In our study, we assessed the genetic relationships and calculated pairwise genetic distances among over two million Single Nucleotide Polymorphisms (SNPs) employing neighbor-joining tree analysis. Utilizing the genetic information from all 130 samples, we constructed a phylogenetic tree that was categorized into two principal groups. The first group encompassed wild walnuts from the Ili and Dushanbe regions, while the second group was comprised of wild walnuts originating from four distinct wild fruit forest reserves in Kyrgyzstan (Figure 3B.). The phylogenetic trees we developed provide a vivid illustration of the genetic affinities among these six geographically diverse regions. A striking observation was that the wild walnuts from Xinjiang demonstrated a closer genetic affinity to those from the Dushanbe region, despite their significant geographic separation. This result is congruent with the classification outcomes derived from both STRUCTURE and PCA analyses, suggesting a robust genetic relationship that transcends geographical barriers.
3.3 Genetic diversity and linkage disequilibrium in the genomes of walnuts
The observed diversity patterns revealed genetic disparities between walnuts from different regions. Nucleotide diversity (θπ) of walnut reveals genetic diversity in different regions. Xinjiang had the lowest nucleotide diversity (θπ = 0.0065) and Sary-Cheek had the highest nucleotide diversity, (θπ = 0.001076). Further, the results showed that the lowest genetic diversity was found in the wild walnut population from Xinjiang, and the highest genetic diversity was found in the Sary-Cheek population (Figure 4A.). The fixation index (FST) is an effective measurement of differences across populations by reflecting the allelic heterozygosity level within a population. The FST index of Y & K is the highest at 0.569014 and also the FST index of Xinjiang, China (Group Y) and other groups show high values all above 0.5. However, the FST index for K & J is the lowest at 0.10313, and all the FST indexes between the Kok-Tundy Area C (Group L), Arslanbob (Group K), Sary-Cheek (Group J), Kara-Alma area (Group M) groups show lower values. Genome level comparisons show that the wild walnut population from Xinjiang has the highest genetic differentiation as compared to the other regions while the lowest genetic differentiation was between the Arslanbob (Group K) and Sary-Cheek (Group J) groups (Figure 4A.).
Figure 4
In plant-forward genetics studies, determining the linkage disequilibrium (LD) pattern is essential. Factors such as the population size of the species, fertility (selfing or hybridization), selection pressure, and the rate of recombination all affect the LD level. We estimated the LD decay for each population as a function of physical distance. Then, we sorted SNP pairs in 100 kb-bins based on the distance between pairs and estimated the average r 2 values for each bin. The walnut tree populations from different localities exhibited differing linkage disequilibrium distances (LD) and linkage disequilibrium attenuation (Figure 4B.). Wild walnuts from the Dushanbe region reached 0.2 r2 at about 100 kb, wild walnuts from the Sary-Cheek (Group J), Arslanbob (Group K), and Kara-Alma area (Group M) regions reached 0.2 r2 at 300–400 kb, and wild walnuts from the Kara-Alma area (Group M) region reached 0.2 r2 at 500 kb. Wild walnuts from the Ili region were more heterogeneous, with their r2 values failing to reach half of the value up to 500 kb. Moreover, the r2 values were higher, showing a stronger linkage imbalance. Notably, the wild walnut population from the Dushanbe region of Tajikistan showed a much shorter average distance for LD to decay to 50% of the maximum compared to endemic walnuts. In contrast, the slowest LD decay was in the wild walnut population from Ili, China.
3.4 Demographic history
By employing the pairwise sequential Markovian coalescence, PSMC, we inferred the demographic histories of four representative individuals from each of the six groups in our sampling sites. Twenty-four individuals from these six sampling sites yielded similar inferred population histories, although the number of effective groups prior to 0.2 Ma was not the same. Our analysis revealed that the effective population size (Ne) was higher around 6 Ma and underwent two large declines at a later time. The first decline occurred at about 2 Ma, followed by a period of stability, and then a more dramatic decline at 0.2 Ma. Interestingly, the Persian walnut Ne curves converged at 0.2Ma; this convergence is likely reflecting the time when all sources of J. regia last shared a common ancestor (Figure 5).
Figure 5
4 Discussion
The publication of a high-quality walnut genome advances research on walnut and other related valuable species (). This includes structural analysis, genetic map construction, and molecular marker-assisted breeding (; ; ). However, until an updated reference genome is published, there is limited information available about the genome (; ). Sequencing the walnut genome has enabled us to identify SNPs efficiently. This study identified 2021717 SNPs, and these variants provide useful genomic resources for future genetic differentiation studies. Genetic diversity is crucial for the long-term survival of species. The higher the genetic diversity (), the greater the adaptation of the species to the environment, and conversely, the greater the risk of extinction as adaptation declines ().This study found relatively low genetic diversity in walnuts from the Ili and Dushanbe regions and generally high genetic diversity in walnuts from the four Kyrgyz regions of ary-Cheek, Arslanbob, Kok-Tundy, and Kara-Alma. This difference in genetic diversity levels suggests distinct origins of walnuts from the Ili or Dushanbe regions and that of the four regions of Kirghizstan. This hypothesis is consistent with the findings of other recent studies (; ).
Genetic diversity is influenced by a multitude of factors, such as reproductive mechanisms, genetic mutations, natural selection, and the isolation of populations, which can be exacerbated by environmental shifts and anthropogenic interference (; ; ). Wind pollination, coupled with self-incompatibility, can inadvertently enhance the likelihood of inbreeding, potentially diminishing the heterozygosity of alleles within a population (). This, in turn, may result in a lack of sufficient heterozygosity to maintain high genetic diversity (). Our sampling map indicates that the four regions of Kyrgyzstan were sampled in close proximity to each other. This reduced geographic isolation could be a contributing factor to their higher genetic diversity. The results also show large differences in prior differentiation among populations. The highest index of genetic differentiation was found between the wild walnut populations in Ili, China while lower indices of genetic differentiation were observed among the four regions of ary-Cheek, Arslanbob, Kok-Tundy, and Kara-Alma. Genetic differentiation between species is influenced by a variety of factors, such as gene flow, reproductive practices, and different natural environments. Fragmentation and isolation can lead to restricted gene flow and increased inbreeding in many tree species ().
As a relict plant of the Tertiary period, extensive fossil research suggests that the ancestor of the walnut was widely distributed across Asia and Europe from the Miocene era. This evidence extends from the Late Tertiary to the Quaternary period, with distribution sites extending from Europe to South Asia in the east (; ; ; ; ; ; ; ; ). Based on other studies, it is presumed the Persian walnut has survived in LGMs in western Europe, Xinjiang Province, and northeastern to southwestern China over a broad geographic area from 30°N to 45°N (; ; ). During the Pleistocene period, population differentiation may have been influenced by repeated extinctions and colonization during the Quaternary Climate Oscillation (). Thus, the Late Miocene and Pleistocene climatic changes were distributed in broadly similar but fragmented environments across Asia. This fragmented distribution has likely contributed to the walnut’s regionally adaptability, with its effective population size continuously limited by climatic fluctuations and its own long generation time (; ; ; ). This result is also reflected in the effective population size dynamics in the PSMC plot. The Himalayan movement, marked by the collision and convergence of the Indian plate with the Eurasian continent, significantly influenced the formation of the Tien Shan mountain basin structure (; ; ). The formation of the Tibetan Plateau and the Tien Shan’s geographic structure both affect the local climate pattern. As a result, the geographical segregation of the walnut populations along the Tien Shan has become increasingly pronounced (; ). Based on the results of our Structure, PCA, and phylogenetic tree, these 130 wild walnut materials from Central Asia were classified into four major taxa. The wild walnuts from the Ili region were the first to be separated in Structure, whereas the four regions of Kyrgyzstan, Sary-Cheek, Arslanbob, Kok-Tundy, Kara-Alma were not as distinctly separated. The Structure results also indicated that walnuts from four regions, Sary-Cheek, Arslanbob, Kok-Tundy, and Kara-Alma, showed heterozygous genetic information, suggesting gradual genetic infiltration between walnuts from the four regions.
PCA and Phylogenetic Tree analysis indicate that the walnut populations from the Ili region and Dushanbe region are more closely related while the genetic composition of walnut in the Kirghizia regions are more complex. Thus, the genetic spatial structure can’t be explained by using purely geographic space. Given the extensive history of human utilization of walnuts, and the clear historical record of human dispersal of walnuts (; ), it is essential that our structural analysis of the genetic structure of walnuts incorporates the influence of anthropogenic factors. Walnut shells were found in a grave site in the South Caucasus (4,230–3,790 B.C., Late Neolithic) () and in the vicinity of wheat and barley grains in the Kashmir region (2,700–2,000 B.C.) (). Further, nutshells were found in the remains of the 3200-year-old Barikot in the Swat Valley of Pakistan deal (). Research suggests that walnuts were first domesticated in the Iranian-Anatolian region and then spread westward and eastward by humans (; , ; ).
Archaeological evidence from the Xinjiang region suggests the utilization of walnuts in Xinjiang can be dated back to the Early Iron Age (seventh to third centuries B.C.) A small number of peach kernels, apricot kernels, and other relics were found in the disturbed soil layer of M1 of the Alagou cemetery in Urumqi (510 B.C.-390.B.C.) (). Emperor Wu of the Han Dynasty came to Central Asia to fight the Xiongnu by connecting the ancient routes in 138 B.C. and 119 B.C. This act opened the first step of the northern route of the Silk Road (). The rise of the Silk Road (; ) facilitated the widespread trade in walnuts across Asia and Europe (; ). A number of sites containing walnuts have been discovered in the silk road following the emergence of the Han Dynasty and in the Tarim Basin side of the North-South-East edge (; ; ; ; ; ; ; ). This evidence indicates that walnuts were widely planted in the Xinjiang region during the Han Dynasty.
According to the above information, the use of walnuts in the Xinjiang region was limited before the Han Dynasty as few walnut artifacts were unearthed from that time. However, a large number of walnut artifacts were unearthed since the Han Dynasty, and walnuts should have been domesticated by human beings before the Silk Road (). Thus, with the emergence of the Silk Road, walnuts were widely cultivated and utilized in Xinjiang. Although walnuts are crops, they do not meet the criteria of being domesticated as they can grow normally even when left out of human influence. Over the course of thousands of years of history, humans selectively bred walnuts by choosing suitable seedlings from different geographical locations (; ; ; ; ; ). Therefore, we must consider human impact on walnut genetic space as well as the influence of pre-existing geographic differentiation.
Reviewing paleoecological data in southern Europe, Syria, Kyrgyzstan, China, and the Himalayas, suggests that Juglans may have survived the last ice age (; ; ; ; ; ; ). According to results from the PSMC, during the Quaternary glacial period (1–2 Ma), the effective population sizes of Juglans declined due to glacial-induced conditions and have maintained a decline since then. Notably, they converged around 0.2 Ma, potentially sharing a common ancestor. However, the formation of genetic spatial patterns observed in this study may be more influenced by human activities than natural factors. Results from Structure, PCA, and phylogenetic tree analyses for walnuts from four regions of Kirghizia show that although the four sites were collected in close proximity to each other, walnuts from Kok-Tundy could be isolated independently. In contrast, those from the other three sites could not. Similar findings have been reported in studies using SSR ().
The wild walnut populations in Kirghizstan existed for at most 2,000 years (). Walnuts from different regions met in Kirghizstan, and the Silk Road further influenced the genetic space of wild walnuts there. Anthropogenic impacts, likely stemming from the extensive cultivation of common walnuts for nut production and for reforestation and afforestation initiatives, have been documented since the 5th century AD (). Although the results of the developmental trees indicate that walnuts from the Ili and Dushanbe regions are closely related, there is no clear evidence that the walnuts from Ili and Dushanbe have existed here since the Tertiary period. This uncertainty is partly attributable to the limited availability of fossil walnut records in Tajikistan (; ).
One thing is certain: despite human activities, the walnut species is still evolving (). Our results suggest that geographic isolation plays a significant role in the differences between walnut populations. In particular, the wild walnut population in the Ili region has low genetic diversity and faces under high selection pressure, compared to high divergence results in other regions. In contrast, the wild walnut population in the Dushanbe region, which is closely related to Ili, was not subjected to strong selection pressure. Given the characteristics of walnut reproduction, including wind-borne and seed propagation, and their long-lived generations (), this selection pressure may have existed for an extended period. Therefore, we need to consider a deeper perspective to conserve the existing wild walnut forests. We could use existing natural conditions to select specific adaptive traits, but the strong anthropogenic influence on walnuts in Central Asia necessitates a different approach. It may be more informative to compare these wild walnuts with modern cultivars to select traits that have been domesticated by humans (GWAS analysis) to elucidate the reasons for the current distribution pattern of wild walnuts in Central Asia.
5 Conclusions
This study analyzed the genetic structure of wild walnuts in Central Asia by utilizing SNP markers. Our results revealed high levels of genetic diversity in Central Asia. Specifically, we found that wild walnuts from Ili and Dushanbe are more closely related to each other and likely share a similar origin. These populations differ in origin from walnuts from the four regions of Kyrgyzstan, yet there is evidence of gene flow with each other. Our findings also support the previous view that wild walnut was geographically isolated by Quaternary autochthonous and climatic changes. Furthermore, the widespread use of the wild walnut in Central Asia by humans, particularly after the emergence of the Silk Road, led to greater anthropogenic disturbances and contributes to the present genetic pattern.
Statements
Data availability statement
The data presented in the study are deposited in the National Center for Biotechnology Information (NCBI) repository, accession number PRJNA1072832 and PRJNA1107910.
Author contributions
XL: Data curation, Methodology, Resources, Visualization, Writing – original draft, Writing – review & editing. XW: Formal Analysis, Investigation, Writing – original draft. DZ: Project administration, Supervision, Writing – original draft. JH: Software, Writing – review & editing. WS: Writing – original draft, Writing – review & editing, Conceptualization, Data curation, Funding acquisition, Methodology, Resources, Supervision. JW: Writing – review & editing, Conceptualization, Formal Analysis, Funding acquisition, Project administration, Software.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This research was funded by the Strategic Biological Resources Capacity Building Project, CAS (grant number KFJ-BRP-017-72); Western Young Scholar, CAS (grant number 2021-XBQNXZ-010); Chinese Academy of Sciences (CAS) Talent Program (grant number 2021-XBQNXZ-013) and Xinjiang Forestry and Grassland Bureau (grant number DZXJZB202203S).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2024.1394409/full#supplementary-material
Supplementary Table 1130 samples sequenced data output.
Supplementary Table 2130 Sequencing Sample Comparison and Coverage Statistics.
Supplementary Table 3Statistics of SNP detection results.
References
1
AbboS.vanOssR. P.GopherA.SarangaY.OfnerR.PelegZ. (2014). Plant domestication versus crop evolution: a conceptual framework for cereals and grain legumes. Trends Plant Sci.19, 351–360. doi: 10.1016/j.tplants.2013.12.002
2
AkçaY.SenS.M. (1994). “Studies on selection of walnut (Juglans regia L.) in Gürün,” in Progress in temperate fruit breeding (Springer, Berlin/Heidelberg, Germany), 179–181.
3
AradhyaM.VelascoD.IbrahimovZ.ToktoralievB.MaghradzeD.MusayevM.et al. (2017). Genetic and ecological insights into glacial refugia of walnut (Juglans regia L.). PloS One12, 1–28. doi: 10.1371/journal.pone.0185974
4
AradhyaM.WoesteK.VelascoD. (2009). “Genetic diversity, structure and differentiation in cultivated walnut (Juglans regia L.),” in VI International Walnut Symposium, Vol. 861, 127–132.
5
Archaeology (1995). New harvest of cultural relics and archaeology in Xinjiang 1979–1989 (China, Xinjiang, Urumqi: Xinjiang People’s Publishing House).
6
Archaeology and Museum (1997). New Harvest of Cultural Relics and Archaeology in Xinjiang Continued 1990–1996 (China, Xinjiang,Urumqi: Xinjiang Fine Arts Photography Press).
7
AubaileF. (2012). Pathways of diffusion of some plants and animals between Asia and the Mediterranean region Vol. 1) (Paris, France: Revue d’ethnoécologie), Vol. 1. doi: 10.4000/ethnoecologie
8
BaiW. N.YanP. C.ZhangB. W.WoesteK. E.LinK.ZhangD. Y. (2018). Demographically idiosyncratic responses to climate change and rapid Pleistocene diversification of the walnut genus Juglans (Juglandaceae) revealed by whole-genome sequences. New Phytol.217, 1726–1736. doi: 10.1111/nph.14917
9
BayazitS.KazanK.GülbittiS.CevikV.AyanoğluH.ErgülA. (2007). AFLP analysis of genetic diversity in low chill requiring walnut (Juglans regia L.) genotypes from Hatay, Turkey. Sci. Hortic.111, 394–398. doi: 10.1016/j.scienta.2006.11.006
10
BeerR.KaiserF.SchmidtK.ArnmannB.CarraroG.GrisaE.et al. (2008). Vegetation history of the walnut forests in Kyrgyzstan (Central Asia): natural or anthropogenic origin? Quaternary Sci. Rev.27, 621–632. doi: 10.1016/j.quascirev.2007.11.012
11
BernardA.LheureuxF.DirlewangerE. (2018). Walnut: past and future of genetic improvement. Tree Genet. Genom.14, 1. doi: 10.1007/s11295-017-1214-0
12
BerryE. (1912). Notes on the geological history of the walnuts and hickories. Plant World15, 225–240.
13
BottemaS. (1980). On the history of the walnut (Juglans regia L.) in Southeastern Europe. Acta Bot. Neerl.29, 343–349. doi: 10.1111/plb.1980.29.issue-5-6
14
CarrionJ. S.Sanchez-GomezP. (1992). Palynological data in support of the survival of walnut (Juglans regia L.) in the western mediterranean area during last glacial times. J. Biogeogr.19, 623–630. doi: 10.2307/2845705
15
ChenF.ChenJ.WangZ.ZhangJ.LiX.LinM.et al. (2019). Genomics: cracking the mysteries of walnuts. J. Genet.98, 33. doi: 10.1007/s12041-019-1084-3
16
ChengW. (1994). Atlas of Chinese Agricultural Archaeology (China, Jiangxi,Nanchang: Jiangxi Science and Technology Press).
17
ChomickiG.SchaeferH.RennerS. S. (2020). Origin and domestication of Cucurbitaceae crops: insights from phylogenies, genomics and archaeology. New Phytol.226, 1240–1255. doi: 10.1111/nph.16015
18
ChristianD. (2000). Silk roads or steppe roads? silk roads World history J. World history1, 1–26. doi: 10.1484/srs-eb
19
CostaJ.CarrapatosoI.OliveiraM. B. P. P.MafraI. (2014). Walnut allergens: molecular characterization, detection and clinical relevance. Clin. Exp. Allergy44, 319–341. doi: 10.1111/cea.12267
20
De CandolleA. (1855). Géographie botanique raisonnée ou exposition des faits principaux et des lois concernant la distribution géographique des plantes de l’époque actuelle. (French, Paris: V. Masson), Vol. 2, 30–40.
21
DePristoM. A.BanksE.PoplinR.GarimellaK. V.MaguireJ. R.HartlC.et al. (2011). A framework for variation discovery and genotyping using next-generation DNA sequencing data. Nat. Genet.43, 491–498. doi: 10.1038/ng.806
22
DingY.CaoY.ZhangW.ChenJ.LiuJ.LiP.et al. (2022). Population-genomic analyses reveal bottlenecks and asymmetric introgression from Persian into iron walnut during domestication. Genome Biol.23, 1. doi: 10.1186/s13059-022-02720-z
23
DodonovA. E.BaiguzinaL. L. (1995). Loess stratigraphy of Central Asia: Palaeoclimatic and palaeoenvironmental aspects. Quaternary Sci. Rev.14, 707–720. doi: 10.1016/0277-3791(95)00054-2
24
DodonovA. E.SadchikovaT. A.SedovS. N.SimakovaA. N.ZhouL. P. (2006). Multidisciplinary approach for paleoenvironmental reconstruction in loess-paleosol studies of the Darai Kalon section, Southern Tajikistan. Quaternary Int.152, 48–58. doi: 10.1016/j.quaint.2005.12.001
25
DongQ.CaoX.YangG.YuH.NicholasK. K.WangC.et al. (2010). Discovery and characterization of SNPs in Vitis vinifera and genetic assessment of some grapevine cultivars. Sci. Hortic.125, 233–238. doi: 10.1016/j.scienta.2010.03.023
26
EarlD. A.VonHoldtB. M. (2012). STRUCTURE HARVESTER: a website and program for visualizing STRUCTURE output and implementing the Evanno method. Conserv. Genet. Resour.4, 359–361. doi: 10.1007/s12686-011-9548-7
27
ErdeiB.MagyariE. (2011). Late Miocene plant remains from Bükkábrány, Hungary. Studia Botanica Hungarica42, 135–151.
28
EvannoG.RegnautS.GoudetJ. (2005). Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol. Ecol.14, 2611–2620. doi: 10.1111/j.1365-294X.2005.02553.x
29
FengX.ZhouH.ZulfiqarS.LuoX.HuY.FengL.et al. (2018). The phytogeographic history of common walnut in China. Front. Plant Sci.9, 1–19. doi: 10.3389/fpls.2018.01399
30
Filipova MarinovaM. V.KvavadzeE. V.ConnorS. E.SjogrenP. (2010). Estimating absolute pollen productivity for some European Tertiary-relict taxa. Vegetation History Archaeobotany19, 351–364. doi: 10.1007/s00334-010-0257-z
31
FjellstromR. G.ParfittD. E. (1994). Walnut (Juglans spp) genetic diversity determined by restriction-fragment-length-polymorphisms. Genome37, 690–700. doi: 10.1139/g94-097
32
FrancisG.SpogliF.FunahashiN.MukaiG.MohnotA.SekiguchiR. (2006). Along the silk road, stanford program on international and cross-cultural education (SPICE) and the silk road project. (New York, USA: Oxford University Press).
33
FrenchD. (1998). Pre-and early-Roman roads of Asia Minor: The Persian royal road. Iran36, 15–43. doi: 10.2307/4299973
34
FullerD. Q. (2007). Contrasting patterns in crop domestication and domestication rates: Recent archaeobotanical insights from the old world. Ann. Bot.100, 903–924. doi: 10.1093/aob/mcm048
35
FullerD. Q.AllabyR. G.StevensC. (2010). Domestication as innovation: the entanglement of techniques, technology and chance in the domestication of cereal crops. World Archaeol.42, 13–28. doi: 10.1080/00438240903429680
36
FullerD. Q.DenhamT.Arroyo-KalinM.LucasL.StevensC. J.QinL.et al. (2014). Convergent evolution and parallelism in plant domestication revealed by an expanding archaeological record. Proc. Natl. Acad. Sci. United States America111, 6147–6152. doi: 10.1073/pnas.1308937110
37
GaisbergerH.LegayS.AndreC.LooJ.AzimovR.AalievS.et al. (2020). Diversity under threat: connecting genetic diversity and threat mapping to set conservation priorities for juglans regia L. Populations in central Asia. Front. Ecol. Evol.8, 1–18. doi: 10.3389/fevo.2020.00171
38
GautB. S.SeymourD. K.LiuQ.ZhouY. (2018). Demography and its effects on genomic variation in crop domestication. Nat. Plants4, 512–520. doi: 10.1038/s41477-018-0210-1
39
GunnB. F.AradhyaM.SalickJ. M.MillerA. J.YongpingY.LinL.et al. (2010). Genetic variation in walnuts (Juglans regia and J. sigillata; Juglandaceae): species distinctions, human impacts, and the conservation of agrobiodiversity in Yunnan, China. Am. J. Bot.97, 660–671. doi: 10.3732/ajb.0900114
40
GuochaoZ. U. O.ZuohengZ.ZhiliangW.MinL.LongshengW. (2008). Tectonic division, stratigraphical stystem and the evolution of western Tianshan mountains, Xinjiang. Geological Rev.54, 748–767. doi: 10.16509/j.georeview.2008.06.007
41
GuptaA. K. (2004). Origin of agriculture and domestication of plants and animals linked to early Holocene climate amelioration. Curr. Sci.87, 54–59.
42
HarrisD. R. (2001). Domestication of plants in the old world: The origin and spread of cultivated plants in west Asia, Europe and the Nile valley. Agr. Hist. Rev.49, 226–227.
43
HemeryG. E.SavillP. S.ThakurA. (2005). Height growth and flushing in common walnut (Juglans regia L.): 5-year results from provenance trials in Great Britain. Forestry78, 121–133. doi: 10.1093/forestry/cpi012
44
HillW.RobertsonA. (1968). Linkage disequilibrium in finite populations. Theor. Appl. Genet.38, 226–231. doi: 10.1007/BF01245622
45
IvanovD.UtescherT.AshrafA. R.MosbruggerV.BozukovV.DjorgovaN.et al. (2012). Late miocene palaeoclimate and ecosystem dynamics in Southwestern Bulgaria - A study based on pollen data from the gotse-delchev basin. Turkish J. Earth Sci.21, 187–211. doi: 10.3906/yer-1004-45
46
JanickJ. (2002). “History of Asian horticultural technology,” in 26th International Horticultural Congress, Toronto, Canada. 19–32.
47
JanickJ. (2003). History of Asian horticultural technology. Acta Hortic.620, 19–32. doi: 10.17660/ActaHortic.2003.620.1
48
JinQ.GaoS.MoR.ShengF.ZhangQ.WuC.et al. (2022). A preliminary study for identifying genes associated with pellicle development in Xinjiang walnut (Juglans regia L.). Horticulturae8, 1–15. doi: 10.3390/horticulturae8090784
49
JinyiL. I.KezhuoW.YapingL. I.GuihuaS. U. N.ChunhuaC. H. U.LiqunL. I.et al. (2006). Geomorphological features, crustal composition and geological evolution of the Tianshan Mountains. Geological Bull. China25, 895–909.
50
KhanazarovA. A.ChernovaG. M.VenglovskiyB. I.AshimovK. S.KenjabaevS. K. (2003). Man-made green monuments of Central Asia: some examples of Uzbekistan. Schweizerische Z. fur Forstwesen154, 226–228. doi: 10.3188/szf.2003.0226
51
Kovar EderJ.KvacekZ.MartinettoE.RoironP. (2006). Late Miocene to Early Pliocene vegetation of southern Europe (7–4Ma) as reflected in the megafossil plant record. Palaeogeogr. Palaeoclimatol. Palaeoecol.238, 321–339.
52
KvacekZ.HablyL.SzakmányG. (1994). Additions to the pliocene flora of gérce (western Hungary). Földtani Közlöny124, 69–87.
53
Lam-TungN.SchmidtH. A.von HaeselerA.Bui QuangM. (2015). IQ-TREE: A fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol. Biol. Evol.32, 268–274. doi: 10.1093/molbev/msu300
54
LiH. (2011). A statistical framework for SNP calling, mutation discovery, association mapping and population genetical parameter estimation from sequencing data. Bioinformatics27, 2987–2993. doi: 10.1093/bioinformatics/btr509
55
LiH.DurbinR. (2009). Fast and accurate short read alignment with Burrows–Wheeler transform. Bioinformatics25, 1754–1760. doi: 10.1093/bioinformatics/btp324
56
LiH.DurbinR. (2011). Inference of human population history from individual whole-genome sequences. Nature475, 493–496. doi: 10.1038/nature10231
57
LiH.HandsakerB.WysokerA.FennellT.RuanJ.HomerN.et al. (2009). The sequence alignment/map format and SAMtools. Bioinformatics25, 2078–2079. doi: 10.1093/bioinformatics/btp352
58
LiY.LiX.CaoH.LiC.JiangH.LiC. (2013). Grain remains from archaeological sites and development of oasis agriculture in Turpan, Xinjiang. Sci. Bull.58, 40–45. doi: 10.1360/tb-2013-suppl017
59
LiangM.DongQ.ZhangX.LiuY.LiH.GuoS.et al. (2023). Metabolomics and transcriptomics analyses reveals the molecular regulatory mechanisms of walnut (Juglans regia L.) embryos in response to shade treatment. Int. J. Mol. Sci.24, 1–26. doi: 10.3390/ijms241310871
60
LiuX. (2005). A study of the crustal tectonic evolution in basin-mountain areas of northwest China in the light of Xinjiang geotransect. Acta Geoscientia Sin.26, 105–112.
61
LuY.ZhangC.LiX.LiangY.WangY.LiW. (2020). Development of EST-SSR markers and their application in the analysis of the genetic diversity of Sophora japonica Linn. Trees-Structure Funct.34, 1147–1156. doi: 10.1007/s00468-020-01985-w
62
LuoX.ZhouH.CaoD.YanF.ChenP.WangJ.et al. (2022). Domestication and selection footprints in Persian walnuts (Juglans regia). PloS Genet.18, 1–20. doi: 10.1371/journal.pgen.1010513
63
ManiB. (2008). Kashmir Neolithic and early Harappan: a linkage. Pragdhara18, 229–247.
64
MannersV.KumariaS.TandonP. (2013). SPAR methods revealed high genetic diversity within populations and high gene flow of Vanda coerulea Griff ex Lindl (Blue Vanda), an endangered orchid species. Gene519, 91–97. doi: 10.1016/j.gene.2013.01.037
65
ManningW. E. (1978). The classification within the Juglandaceae. Ann. Missouri botanical garden65, 1058–1087. doi: 10.2307/2398782
66
McGranahanG.LeslieC. (1991). Genetic Resources of Temperate Fruit and Nut Crops (Wageningen, The Netherlands: International Society of Horticultural Sciences).
67
MercuriA. M.MarignaniM.SadoriL. (2013). 2013 Palynology: the bridge between palaeoecology and ecology for the understanding of human-induced global changes in the mediterranean area. Annali di Botanica3, 107–113. doi: 10.4462/annbotrm-10332
68
MeyerR. S.PuruggananM. D. (2013). Evolution of crop species: genetics of domestication and diversification. Nat. Rev. Genet.14, 840–852. doi: 10.1038/nrg3605
69
MillerA. J.GrossB. L. (2011a). From forest to field: perennial fruit crop domestication. Am. J. Bot.98, 1389–1414. doi: 10.3732/ajb.1000522
70
MillerA. J.GrossB. L. (2011b). vFrom forest to field: Perennial fruit Crop Domestication. Am. J. Bot.98, 1389–1414. doi: 10.3732/ajb.1000522
71
MishraR. K. (2020). The ‘Silk road’: historical perspectives and modern constructions. Indian Historical Rev.47, 21–39. doi: 10.1177/0376983620922431
72
MolnarT.ZaurovD.CapikJ.EisenmanS.FordT.NikolyiL.et al. (2011). Persian walnuts (Juglans regia L.) in Central Asia. Annu. Rep. North Nut Grow Assoc.101, 56–69.
73
Nadachowska-BrzyskaK.BurriR.SmedsL.EllegrenH. (2016). PSMC analysis of effective population sizes in molecular ecology and its application to black-and-white Ficedula flycatchers. Mol. Ecol.25, 1058–1072. doi: 10.1111/mec.13540
74
NiceseF. P.HormazaJ. I.McGranahanG. H. (1998). Molecular characterization and genetic relatedness among walnut (Juglans regia L.) genotypes based on RAPD markers. Euphytica101, 199–206. doi: 10.1023/A:1018390120142
75
PecoraroC.BabbucciM.VillamorA.FranchR.PapettiC.LeroyB.et al. (2016). Methodological assessment of 2b-RAD genotyping technique for population structure inferences in yellowfin tuna (Thunnus albacares). Mar. Genomics25, 43–48. doi: 10.1016/j.margen.2015.12.002
76
PerrinoE. V.WagensommerR. P. (2022). Crop wild relatives (CWRs) threatened and endemic to Italy: urgent actions for protection and use. Biology-Basel11 (2). doi: 10.3390/biology11020193
77
PokhariaA. K.ManiB. R.SpateM.BettsA.SrivastavaA. (2018). Early neolithic agriculture, (2700–2000 B.C.) and Kushan period developments (AD 100–300): macrobotanical evidence from Kanispur in Kashmir, India. Vegetation History Archaeobotany27, 477–491. doi: 10.1007/s00334-017-0645-8
78
PollegioniP.WoesteK.ChiocchiniF.Del LungoS.CiolfiM.OlimpieriI.et al. (2017). Rethinking the history of common walnut (Juglans regia L.) in Europe: Its origins and human interactions. PloS One12, 1–24. doi: 10.1371/journal.pone.0172541
79
PollegioniP.WoesteK. E.ChiocchiniF.Del LungoS.OlimpieriI.TortolanoV.et al. (2015). Ancient humans influenced the current spatial genetic structure of common walnut populations in Asia. PloS One10, 1–16. doi: 10.1371/journal.pone.0135980
80
PollegioniP.WoesteK. E.ChiocchiniF.OlimpieriI.TortolanoV.ClarkJ.et al. (2014). Landscape genetics of Persian walnut (Juglans regia L.) across its Asian range. Tree Genet. Genom.10, 1027–1043. doi: 10.1007/s11295-014-0740-2
81
PopovM. (1929). Wild growing fruit trees and shrubs of Asia Minor. Bull. Appl. Bot. Pl Breed22, 241–483.
82
PotterD.GaoF. Y.AielloG.LeslieC.McGranahanG. (2002). Intersimple sequence repeat markers for fingerprinting and determining genetic relationships of walnut (Juglans regia) cultivars. J. Am. Soc Hortic. Sci.127, 75–81. doi: 10.21273/JASHS.127.1.75
83
PrimmerC.BorgeT.LindellJ.SætreG. P. (2002). Single-nucleotide polymorphism characterization in species with limited available sequence information: high nucleotide diversity revealed in the avian genome. Mol. Ecol.11, 603–612. doi: 10.1046/j.0962-1083.2001.01452.x
84
PritchardJ. K.StephensM.DonnellyP. (2000). Inference of population structure using multilocus genotype data. Genetics155, 945–959. doi: 10.1093/genetics/155.2.945
85
PurcellS.NealeB.Todd-BrownK.ThomasL.FerreiraM. A. R.BenderD.et al. (2007). PLINK: A tool set for whole-genome association and population-based linkage analyses. Am. J. Hum. Genet.81, 559–575. doi: 10.1086/519795
86
RenG.ZhangX.LiY.RidoutK.Serrano-SerranoM. L.YangY.et al. (2021). Large-scale whole-genome resequencing unravels the domestication history of Cannabis sativa. Sci. Adv.7, 9. doi: 10.1126/sciadv.abg2286
87
Renault MiskovskyJ.Bui Thi MaiM.GirardM. (1984). A propos de l’indigenat ou de l’introduction de Juglans et Platanus dans l’ouest de l’Europe au Quaternaire. Rev. paleobiologievol. spécial, 155–178.
88
RetiefJ. D. (2000). Phylogenetic analysis using PHYLIP. Methods Mol. Biol. (Clifton N.J.)132, 243–258. doi: 10.1385/1-59259-192-2:243
89
RoironP. (1991). The upper Miocene macroflora from diatomites of Murat (Cantal, France) paleoclimatic implications. Paleontogr Abt B223, 169–203.
90
Rong-TingX. (1989). “Discussion on the origin of walnut in China,” in I International Symposium on Walnut Production, Vol. 284, 353–362.
91
RoorW.KonradH.MamadjanovD.GeburekT. (2017). Population Differentiation in Common Walnut (Juglans regia L.) across Major Parts of Its Native Range-Insights from Molecular and Morphometric Data. J. Hered.108, 391–404. doi: 10.1093/jhered/esw122
92
ShatilovaI.RukhadzeL.KokolashviliI. (2014). The history of genus Juglans L. on the territory of Georgia. Bull. Georg. Natl. Acad. Sci.8, 1–7.
93
SpenglerR. N.IIITangL.NayakA.BoivinN.OlivieriL. M. (2021). The southern Central Asian mountains as an ancient agricultural mixing zone: new archaeobotanical data from Barikot in the Swat valley of Pakistan. Vegetation History Archaeobotany30, 463–476. doi: 10.1007/s00334-020-00798-8
94
TorokeldievN.ZieheM.GailingO.FinkeldeyR. (2019). Genetic diversity and structure of natural Juglans regia L. populations in the southern Kyrgyz Republic revealed by nuclear SSR and EST-SSR markers. Tree Genet. Genom.15, 1. doi: 10.1007/s11295-018-1311-8
95
UNESCO (2024) Silk Road Project. Available online at: https://zh.unesco.org/silkroad/guanyusichouzhilu.
96
VahdatiK. (2013). “Traditions and folks for walnut growing around the silk road,” in 1st International Symposium on Fruit Culture and its Traditional Knowledge along Silk Road Countries, Yerevan, Armenia, 19–24. doi: 10.17660/ActaHortic.2014.1032.1
97
WangB. (1983). Overview of agricultural archaeology in Xinjiang. Agric. Archaeol.01), 102–117 + 102.
98
WangL.HeF.HuangY.HeJ.YangS.ZengJ.et al. (2018). Genome of wild mandarin and domestication history of Mandarin. Mol. Plant11, 1024–1037. doi: 10.1016/j.molp.2018.06.001
99
WangL.JiangL.ZhangG.JiangH. e. (2020). Study on the remains of horticultural plants in tuyog grottoes of Turpan. Agric. Archaeol.06), 15–26.
100
WangZ.ZhangH.TongB.HanB.LiuD.ZhangP.et al. (2023). EST-SSR marker-based investigation on genetic diversity and genetic structure of Juglans mandshurica Maxim. in Shandong Province of China. Genet. Resour. Crop Evol.70, 981–991. doi: 10.1007/s10722-022-01482-8
101
WenY.HanW.WuS. (2010). Plant genetic diversity and its influencing factors. J. Cent. South Univ. For. Technol.30, 80–87. doi: 10.14067/j.cnki.1673-923x.2010.12.002
102
WilkinsonK. N.GasparianB.PinhasiR.AvetisyanP.HovsepyanR.ZardaryanD.et al. (2012). Areni-1 Cave, Armenia: A Chalcolithic-Early Bronze Age settlement and ritual site in the southern Caucasus. J. Field Archaeol.37, 20–33. doi: 10.1179/0093469011Z.0000000002
103
WuY.TuohutiW.LI (2002). Excavation of the Yingpan cemetery in Yuli County, Xinjiang. Archaeology06), 58–74 + 102–104 + 102.
104
XiaohongT.YongW.LifuA. (2014). Urumqi yu’ergou site and alagou (Alwighul) cemetery. Archaeology04, 19–35 + 12.
105
XuX.KleidonA.MillerL.WangS.WangL.DongG. (2010). Late Quaternary glaciation in the Tianshan and implications for palaeoclimatic change: a review. Boreas39, 215–232. doi: 10.1111/j.1502-3885.2009.00118.x
106
YangY.ForsytheE. S.DingY.ZhangD.BaiW. (2023). Genomic analysis of plastid-nuclear interactions and differential evolution rates in coevolved genes across juglandaceae species. Genome Biol. Evol.15, 8. doi: 10.1093/gbe/evad145
107
YaoY.-F.BruchA. A.MosbruggerV.LiC.-S. (2011). Quantitative reconstruction of Miocene climate patterns and evolution in Southern China based on plant fossils. Palaeogeogr. Palaeoclimatol. Palaeoecol.304, 291–307. doi: 10.1016/j.palaeo.2010.04.012
108
YuanX.SunY.BaiX.DangM.FengX.ZulfiqarS.et al. (2018). Population structure, genetic diversity, and gene introgression of two closely related walnuts (Juglans regia and J. sigillata) in Southwestern China revealed by EST-SSR markers. Forests9, 10.
109
ZhangY. (1983). An introduction to ancient crops unearthed in Xinjiang. Agric. Archaeol.01), 122–126 + 118–121.
110
ZhangD. F.FengquanL.JianminB. (2000). Eco-environmental effects of the Qinghai-Tibet Plateau uplift during the Quaternary in China. Environ. Geol.39, 1352–1358. doi: 10.1007/s002540000174
111
ZhangZ.SunJ. (2011). Palynological evidence for Neogene environmental change in the foreland basin of the southern Tianshan range, northwestern China. Global Planet. Change75, 56–66. doi: 10.1016/j.gloplacha.2010.10.006
112
ZhangB.-W.XuL.-L.LiN.YanP.-C.JiangX.-H.WoesteK. E.et al. (2019). Phylogenomics reveals an ancient hybrid origin of the persian walnut. Mol. Biol. Evol.36, 2451–2461. doi: 10.1093/molbev/msz112
113
ZhangJ.ZhangW.JiF.QiuJ.SongX.BuD.et al. (2020). A high-quality walnut genome assembly reveals extensive gene expression divergences after whole-genome duplication. Plant Biotechnol. J.9), 1848–1850. doi: 10.1111/pbi.13350
114
ZhaoP.WoesteK. E. (2011). DNA markers identify hybrids between butternut (Juglans cinerea L.) and Japanese walnut (Juglans ailantifolia Carr.). Tree Genet. Genom.7, 511–533. doi: 10.1007/s11295-010-0352-4
115
ZhouH.YanF.HaoF.YeH.YueM.WoesteK.et al. (2023). Pan-genome and transcriptome analyses provide insights into genomic variation and differential gene expression profiles related to disease resistance and fatty acid biosynthesis in eastern black walnut (Juglans nigra). Horticulture Res.10, 1. doi: 10.1093/hr/uhad015
Summary
Keywords
Juglans regia, SNP, population genomics, gene flow, genetic differentiation
Citation
Li X, Wang X, Zhang D, Huang J, Shi W and Wang J (2024) Historical spread routes of wild walnuts in Central Asia shaped by man-made and nature. Front. Plant Sci. 15:1394409. doi: 10.3389/fpls.2024.1394409
Received
01 March 2024
Accepted
17 May 2024
Published
06 June 2024
Volume
15 - 2024
Edited by
Robert Philipp Wagensommer, Free University of Bozen-Bolzano, Italy
Reviewed by
Elisa Vendramin, Council for Agricultural Research and Economics, Italy
Enrico Vito Perrino, International Centre for Advanced Mediterranean Agronomic Studies, Italy
Updates
Copyright
© 2024 Li, Wang, Zhang, Huang, Shi and Wang.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Wei Shi, shiwei@ms.xjb.ac.cn; Jiancheng Wang, www-1256@ms.xjb.ac.cn
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.