Genetic diversity, population structure, and DNA fingerprinting of Ailanthus altissima var. erythrocarpa based on EST-SSR markers

Ailanthus altissima var. erythrocarpa is an A. altissima variety with high economic, ecological and ornamental value, but there have been no reports on the development of SSR primers for it. According to the SSR primer information provided by the transcriptome of A. altissima var. erythrocarpa, 120 individuals with different redness levels were used to screen polymorphic primers. Transcriptomic analysis revealed 10,681 SSR loci, of which mononucleotide repeats were dominant (58.3%), followed by dinucleotide and trinucleotide repeats (16.6%, 15.1%) and pentanucleotide repeats (0.2%). Among 140 pairs of randomly selected primers, nineteen pairs of core primers with high polymorphism were obtained. The average number of alleles (Na), average number of effective alleles (Ne), average Shannon’s diversity index (I), average observed heterozygosity (Ho), average expected heterozygosity (He), fixation index (F) and polymorphic information content (PIC) were 11.623, 4.098, 1.626, 0.516, 0.696, 0.232 and 0.671, respectively. Nineteen EST-SSR markers were used to study the genetic diversity and population structure of A. altissima var. erythrocarpa. The phylogenetic tree, PCoA, and structure analysis all divided the tested resources into two categories, clearly showing the genetic variation between individuals. The population showed high genetic diversity, mainly derived from intraspecific variation. Among nineteen pairs of primers, 4 pairs (p33, p15, p46, p92) could effectively distinguish and be used for fingerprinting of the tested materials. This study is of great significance for genetic diversity analysis and molecular-assisted breeding of A. altissima var. erythrocarpa.

Ailanthus altissima (Mill) Swingle, belonging to Simaroubaceae and Ailanthus, is widely distributed on all continents except Antarctica and is a pioneer tree used in greening and barren mountain afforestation.It has high economic, medicinal, ecological and ornamental value [1][2][3][4][5][6][7] .Wild related species often become important resources for innovative plant breeding due to their higher resistance and richer genetic characteristics than those of this species.Ailanthus altissima var.erythrocarpa is a variety of A. altissima (Mill) Swingle; it is usually recognized as a single population or multiple very small populations in the northern region of China. A. altissima var.erythrocarpa has significantly better physical and chemical properties 8 and medicinal and ornamental value than A. altissima (Mill) Swingle due to its higher contents of flavonoids and other substances.The fruiting period of A. altissima var.erythrocarpa is from June to September.Because the fruit is red and persists for months, it has high ornamental value and a high profile 8 .There are significant differences in the appearance time, duration and chroma of fruit color (data not yet published).Genetic composition and the environment can lead to stability or variability in plant morphology 9 , which necessitates scientific identification and evaluation of germplasm of A. altissima var.erythrocarpa.
Expressed Sequence Tag-SSRs (EST-SSRs) are derived from mRNA transcription sequences and are applicable to genetic linkage map construction, genetic diversity analysis, population genetic structure research, and studies of species evolution and variation in plants with or without reference genomes [10][11][12] .Compared with g-SSRs,

Number and distribution of EST-SSR loci in A. altissima var. erythrocarpa
By using the Illumina HiSeq™ high-throughput sequencing platform, transcriptome sequencing of A. altissima var.erythrocarpa was performed, and an SSR search of the spliced and assembled unigenes revealed 9828 accurate SSRs and 853 compound SSRs, totaling 10,681 SSR loci.The EST-SSR types of A. altissima var.erythrocarpa were abundant, and their distribution showed (Fig. 1a) trinucleotide repeats, dinucleotide repeats and mononucleotide repeats in the largest numbers, among which 6228 mononucleotides accounted for the highest proportion (58.31%); 1777 dinucleotides accounted for 16.64%; 1621 trinucleotides accounted for 15.18%; and tetra-hexanucleotide repeats were rare, accounting for 0.99%, 0.24% and 0.66%, respectively.Other complex nucleotide repeats accounted for 7.99% of the total SSRs, indicating that mono-trinucleotides were the main repeat motifs of EST-SSRs in A. altissima var.erythrocarpa.
According to EST-SSR repeat motif type analysis, there were 166 exact SSR motifs, and the numbers of mononucleotide to hexanucleotide repeats were 2, 6, 30, 39, 23 and 66, respectively (Table 1).Mononucleotides were mainly the A/T repeat type (6214), accounting for 58.18% of the total number of mononucleotides.Dinucleotides were mainly the AT/AT (33.54%) and AG/CT (32.58%) repeat motifs, accounting for 5.58% and 5.42% of the total SSRs, respectively.The specific percentages of the most common trinucleotides, AGA/TCT, AAG/CTT and CTT/AAG, were 1.40%, 1.35% and 1.33%, respectively.The AAAT/ATTT and AAAG/CTTT tetranucleotides accounted for 0.18% and 0.08%, respectively.The penta-to hexanucleotide repeat motifs appeared in many forms, but their frequency was relatively low.Among all repeat types (Fig. 1b), the number of repeats of SSR loci motifs was 5 or more, and the frequency of nucleotide repeats decreased gradually with an increase in the number of repeats.The types with more than 10 repetitions accounted for 63.16% of the total, the largest percentage.This was followed by those with 5 and 6, accounting for 9.90% and 9.14% of the total, respectively.The lowest percentage was observed for those with 9 repetitions, accounting for only 2.32%.The number of mononucleotide repeats was 10 or more, that of dinucleotides was 6 or more, and the highest repetition rate was 6 times.The trinucleotide

Primer polymorphism analysis
Based on capillary electrophoresis, correlation analysis of 7 natural populations of A. altissima var.erythrocarpa was carried out with microsatellite sites showing polymorphism (Table 3 and Fig. 2).There were 221 alleles identified among the nineteen EST-SSR markers.All 120 individuals could be uniquely genotyped using these 221 alleles, demonstrating the high discrimination capacity of the nineteen EST-SSR markers.The markers showed considerable variation, with the number of alleles (Na) ranging from 7 to 20 (mean 11.632), number of effective alleles (Ne) from 1.550 to 11.942 (mean 4.098), observed heterozygosity (Ho) from 0.191 to 0.931 (mean 0.516), expected heterozygosity (He) from 0.355 to 0.916 (mean 0.696), unbiased expected heterozygosity (uHe) from 0.356 to 0.920 (mean 0.699), Shannon's information index (I) from 0.797 to 2.675 (mean 1.626), and polymorphism information content (PIC) from 0.340 to 0.912 (mean 0.671).The marker with the most alleles was p33 (20), and the markers with the fewest alleles were p30, p40 and p95 (7).The PIC values of the nineteen EST-SSR markers ranged from 0.340 to 0.912, with an average of 0.671.All markers exhibited high polymorphism, with the exception of p30 and p49 (PIC > 0.5).
As shown in Table 4, the average inbreeding coefficients Fis and Fit were − 0.018 and 0.199, respectively, indicating that the number of homozygous individuals was greater than that of heterozygous individuals in the seven natural populations of A. altissima var.erythrocarpa.Among the seven natural populations, the Fst differentiation coefficient of the nineteen EST-SSRs showed a variation range of 0.084-0.456,with an average value of 0.237, which is between 0.15 and 0.25.The variation range of gene flow at polymorphic sites was 0.298-2.737,with an average value of 1.056.

Development of DNA fingerprints for A. altissima var. erythrocarpa
According to the genotype statistics of nineteen pairs of SSR primers, the primers with higher PIC values and more genotypes were selected to distinguish 120 A. altissima var.erythrocarpa accessions.According to the PIC values of the SSR primers, the test materials were added successively from high to low and identified gradually until they were completely separated.Table 5 shows that SSR primers p33 and p15 could distinguish 92.5% of the materials (111).SSR primers p33, p15 and p46 could distinguish 99.17% of the materials (119).Four pairs of SSR primers, p33, p15, p46 and p92, could be used to distinguish the 120 test materials.To facilitate promotion and application, the fingerprint information and names of 120 individuals of A. altissima var.erythrocarpa were imported into QR code software to generate fingerprint QR codes, wherein the primers p33, p15, p46 and p92 were represented by A, B, C and D, respectively (Supplementary Fig. 1).Please see the appendix for all fingerprints.In this study, nineteen pairs of core primers were screened to completely separate 120 test samples of A. altissima var.erythrocarpa, and the obtained DNA fingerprints were used to evaluate and identify the genotypes of this variety at the molecular level.www.nature.com/scientificreports/

Genetic diversity analysis of A. altissima var. erythrocarpa
There were significant differences in genetic diversity among A. altissima var.erythrocarpa from 7 different populations (Table 6 and Fig. 3).Compared with those of A. altissima var.erythrocarpa in the other six populations, the values of Ne, I, He, uHe and F in the HB population were the highest, and the values of Ho were the lowest.There were differences in indicators of genetic diversity between populations: HB (Ne = 4.439) > NX (Ne = 3.330) > SHX (Ne = 3.224) > HN (Ne = 3.202) > BJ (Ne = 2.829) > TJ and SHD (Ne = 1.421).In contrast, the level of genetic diversity in SHD and TJ was low, which may be due to the small number of individuals.The fixation index was less than 0 in populations TJ, BJ and SHD, indicating that the number of heterozygotes was greater than that of homozygotes, while the number of homozygotes was greater than that of heterozygotes in the other four populations.The average paired F ST coefficient among the 7 populations of A. altissima var.erythrocarpa was 0.163.The differentiation between TJ and SHD was the highest (0.600), which may be related to too few samples.The differentiation between BJ and SHD was the second highest (0.342), and the differentiation coefficient between NX and SHX (0.016) was the lowest, but their Nm was the highest (15.401), indicating that the genetic differentiation between the NX and SHX populations was small and that gene exchange is frequent (Fig. 4 and Supplementary Table 2).The range of genetic distances between the seven natural populations of A. altissima var.erythrocarpa was 0.042-1.386,and the variation range of genetic similarity was 0.250-0.938(Table 7).Among the populations, NX and SHX had the highest genetic similarity, and the genetic distance was relatively short, indicating that they were more genetically similar than other pairs.TJ and SHD had the highest genetic similarity and the greatest genetic distance, making them the most distant relatives.
According to the genetic similarity coefficients among the test materials of A. altissima var.erythrocarpa, using MEGA7 software, the UPGMA algorithm was used to constructed a dendrogram for 120 A. altissima var.erythrocarpa materials (Fig. 5a), which could be divided into two groups (A and B) through cluster analysis.Group A was mainly composed of 24 materials from the Handan and Pingshan areas of Hebei Province and Shandong Province, accounting for 20% of the total test materials.There were 96 materials from the 7 natural populations of A. altissima var.erythrocarpa in Group B, accounting for 80% of the total materials tested.In Group B, the test materials from Henan Province were scattered in each subgroup, while those in Hebei Province were more clustered, which may be related to the genetic background of the samples and the environment at the sample collection sites.Using GenAlEx 6.502 software to carry out principal coordinate analysis (Fig. 5b), the variation contribution rates of the first and second principal components were 23.22% and 5.4%, respectively, and the classification results were consistent with those of the cluster analysis.This further proves that the primers  4d) the value of LnP(D) continued to increase, there was no inflection point, and the Delta K values could be used to determine the optimal K value (Supplementary Table 3, Fig. 5c,d).STRU CTU RE analysis showed that the seven natural populations of A. altissima var.erythrocarpa could be genetically divided into two different subpopulations using the ΔK method (Fig. 5c), and the clustering pattern was the same as that in the dendrogram.
When K was equal to 3-5, some individuals in subgroup 1 were divided into several subgroups, indicating that these individuals were hybrids and that the other individuals had a pure pedigree (Fig. 5e).Individuals HX01, WX04, and ZHL04 in subgroup 1 consistently originated from two subgroups.The individuals in subgroup 2 originated from a single ancestral group, except that P15 and P16 were divided among three subgroups (K = 3-5).

Discussion
In recent years, SSR molecular markers have been widely used in studies of species genetic diversity due to their advantages, such as high species specificity and cross-related species transfer, and primer development is their primary use.To date, only JosphatK.Saina et al. 18,33     were obtained based on transcriptome sequencing results of tissue-cultured seedlings from young leaves of A. altissima var.erythrocarpa, which was far lower than that reported by JosphatK.Saina et al. 18 (33,084 potential SSR loci), and the EST-SSR distribution characteristics were similar.The SSR markers developed in this study had a high effective amplification rate, which is helpful for analyzing the genetic background and performing auxiliary breeding of A. altissima var.erythrocarpa.Most plants show mononucleotides, dinucleotides and trinucleotides as dominant repeat motifs, such as Stephanandra incisa 23 , Bougainvillea cultivars 21 , Cocos nucifera L. 34 , A. altissima (Mill.)Swingle 18 , Pinus bungeana 11 , Ulmus pumila L. 35 , and A. altissima var.erythrocarpa.The results of the transcriptome analysis showed the same pattern, with these three types accounting for 90.12% of the total SSRs.AG/CT is the most frequent dinucleotide repeat motif in most plants, and the homopurine-homopyrimidine extension is often found in the 5′ untranslated region, which plays an important role in gene expression and nucleic acid metabolism regulation in plants.AAG/CTT is the main nucleotide repeat type in dicotyledons 18,23,24,33,[36][37][38][39][40][41] .In this study, the dinucleotide AT/AT repeat type had the highest frequency, which was consistent with the results for Pinus bungeana 11,42 , Bougainvillea cultivars 21 .The majority of the trinucleotide repeats were AGA/TCT, followed by AAG/CTT.Among the 140 pairs of primers screened and synthesized, 106 yielded amplification bands, and the effective amplification rate reached 75.71%, which was similar to that in A. altissima (Mill.)Swingle (80%) 18 and Cinnamomum chago (70.59%) 43 .The polymorphism rate of the primers was 13.57%, which was significantly lower than that reported by Saina et al. 18 (65%), which may be caused by the use of more randomly selected primers, small genetic differences between samples and different environmental selection pressure.Therefore, EST-SSRs are still the most economical and effective tool for developing a large number of primers at the same time using transcriptome data.The chloroplast genome of A. altissima is a circular molecule with a size of 160,815 base pairs and a tetrad structure.The chloroplast genome size of A. altissima was similar to that of Arabidopsis thaliana (L.) Heynh.(120-217 kb) and Populus trichocarpa (155-159 kb), so it was relatively easy to develop cpSSR markers for A. altissima.Compared with cpSSRs, EST-SSRs exist in the gene regions of the nuclear genome and www.nature.com/scientificreports/are highly polymorphic 18,44 .In this study, EST-SSRs were combined with molecular markers of A. altissima var.erythrocarpa to study the point mutations in A. altissima var.erythrocarpa, providing important resources for subsequent germplasm resource protection and evolutionary research in populations.Phenotypic traits combined with molecular markers could be used to analyze the correlations of traits, providing guidance for the utilization of germplasm resources and breeding of new varieties.PIC can measure the polymorphism of primers: PIC > 0.5 indicates highly polymorphic, 0.25 < PIC < 0.5 indicates moderately polymorphic, and PIC < 0.25 indicates weakly polymorphic 45 .In this study, 19 pairs of clearly polymorphic primers with the expected size were screened, and the average PIC value was 0.671, lower than the average PIC value of A. altissima (0.819) 18 .PIC > 0.5 indicated that the identified SSR loci had high genetic resolution, making them conducive to population genetic analysis and revealing that the studied A. altissima var.erythrocarpa accessions had high genetic diversity.This may be related to the complex living environment of wild A. altissima var.erythrocarpa.The materials of wild A. altissima var.erythrocarpa used in this experiment are widely distributed in Ningxia, Shanxi, Shandong, Henan, Tianjin, Beijing, Hebei and other places.The sampling points were mostly on hillsides, on roadsides and in other similar places.The diversity of the geographical environment and the characteristics of climate variability jointly determine the diversity of the living environment of A. altissima var.erythrocarpa.Therefore, this study further reveals the genetic diversity and structure of this species.
The average expected heterozygosity is less affected by sampling than other metrics.The higher the value is, the lower the genetic consistency and the higher the level of genetic diversity.The average expected heterozygosity of A. altissima var.erythrocarpa from 7 natural populations was 0.509, similar to the results of Saina et al. 18 , which may be related to the wide distribution and complex living environment of wild A. altissima var.erythrocarpa.The mean He was slightly lower than that reported in previous studies 18,28,30 , which may be related to the fact that A. altissima var.erythrocarpa is a variety of A. altissima and originated in the same place.In the study of Saina et al. 18 , the mean Ne of A. altissima (Mill.)Swingle is 3.757.In this study, the mean Ne of A. altissima var.erythrocarpa is 2.838.This may be due to the fact that the population size ( 14) and sample size (165) of A. altissima (Mill.)Swingle are larger than that of this study (7 and 120), and the distribution range (99-120° N and 29-37° E) is wider than that of A. altissima var.erythrocarpa (105-117° N and 33-41° E).Among the seven natural populations, HB (Ne = 4.439, I = 1.625) had the highest genetic diversity, TJ (Ne = 1.421,I = 0.292) and SHD (Ne = 1.421,I = 0.292) had the lowest genetic diversity.This may be because each group has different site conditions (city streets, mountains, rocks), and its growing environment is different; Samples from city streets, there is human interference.These all influence the measurement of genetic diversity.The F ST comparison between nineteen pairs of EST-SSRs and 7 natural populations (Fst mean = 0.163, Nm mean = 3.783) showed a high degree of genetic differentiation (Table 4 and Fig. 4) and frequent gene exchange, indicating that there was strong gene flow (Nm > 1) among the seven natural populations, which inhibited genetic differentiation caused by genetic drift, with low differentiation level and high intra-population differentiation level.This is supported by PCoA (Fig. 5b) and dendrogram based on the UPGMA algorithm (Fig. 5a).
According to genetic distance and genetic similarity analysis, NX and SHX had the closest genetic relationship.To effectively analyze A. altissima var.erythrocarpa germplasm resources, information about natural habitat, genetic evolution, genetic diversity, and the degree of species mixing is required 46 .For principal coordinate analysis, structure analysis, and cluster analysis based on genetic distance as determined by UPGMA, the seven natural populations of red-fruited Ailanthus could be divided into two subgroups; one was subgroup A dominated by materials from Handan and Pingshan County, Hebei Province, and the other was subgroup B (Fig. 5a,b,e) dominated by materials other than those from SHD, which may be because the material for transcriptome sequencing was taken from Handan, and subgroup A was highly represented among the samples.Structurebased clustering is a reliable tool for estimating the degree of species mixing, and the kinship of a small number of samples in subpopulation A and subpopulation B was complex (Fig. 5e), indicating that during the evolution of the species, frequent gene exchange and recombination between individuals in the region led to a richer genetic structure.Some germplasm resources are introduced to different places via wind and animal and human activities.NX, SHX, HN, TJ, and BJ were clustered in subgroup A, which may not be geographically isolated, and subgroup A and subgroup B included samples of the HB natural population.When K = 3-5, individuals HX01, WX04, and ZHL04 in subgroup 1 have complex lineages and are always related to both subpopulations.Individuals in subgroup 2, except for P15 and P16, which can be divided into three subgroups, all originated from a single ancestral group, these results can explain that there is no geographical separation between the 7 natural groups, and the kinship and complexity of each individual can be clearly known.The resources of the same area were almost always clustered together (Fig. 5a,b,e).This finding is similar to the classification results of Astragali Radix 47 , indicating high genetic similarity and frequent gene exchange.
DNA fingerprinting technology analyzes the genetic material itself, avoiding the influences of environmental factors, sample morphology, material source and other factors, and identifies the germplasm genotypes as indicators at the molecular level.It has been widely used in grain, oil, vegetable and fruit crops, flowers, shrubs, arbors and other plants 16,20,21,48 .DNA fingerprinting has high simplicity and stability, variability, and multiple loci and can distinguish members of the same family, such as different species, varieties, and forms and even the same strain and subtle variation between individuals.It can be used in germplasm genetic relationship inference, genetic breeding, population genetic structure determination, ecological and evolutionary studies, classification and other valuable genetic marker methods [48][49][50] .However, there have been no reports on the fingerprint information of A. altissima var.erythrocarpa.In this study, four of nineteen EST-SSR core primers (p33, p15, p46 and p92) were screened to distinguish 120 tested A. altissima var.erythrocarpa materials and construct a unique molecular ID card.The establishment of fingerprinting for A. altissima var.erythrocarpa is of great significance for its subsequent breeding, as the resource can be used to identify varieties and excellent genes as well as their variation and transmission in offspring, and lays a foundation for the establishment and improvement www.nature.com/scientificreports/ of fingerprinting databases of A. altissima var.erythrocarpa variety resources, genetic resource management and intellectual property protection.However, in this study, there were few test materials, the sample collection area was only in northern China, and the fingerprint background of the "core germplasms" of A. altissima var.erythrocarpa lacked in-depth research.As a result, further verification of the reliability of the results is required.

Conclusion
In this study, nineteen pairs of EST-SSRs revealed high genetic diversity and significant genetic structure in A. altissima var.erythrocarpa.Using four pairs of SSR core primers, a total of 120 wild individuals from 7 regions of A. altissima var.erythrocarpa were 100% distinguished and assigned unique identification IDs.The experimental resources of A. altissima var.erythrocarpa have rich genetic variations, mainly derived from intraspecific variations, and do not exhibit geographical variation.The research results presented in this article provide a public resource of SSR primers in A. altissima var.erythrocarpa, in turn providing a theoretical basis for effective management, protection, and utilization of its germplasms.The excellent varieties studied according to the SSR primers of A. altissima var.erythrocarpa can be introduced by biotechnological means and can also be preserved ex situ or in situ, which would be conducive to subsequent research.The findings also provide genetic resources for an in-depth understanding of population genetics and evolutionary history, joint genomics research, the selection of dominant parents for the cultivation of new varieties, and intellectual property protection in A. altissima var.erythrocarpa.

Plant materials
A total of 120 wild individuals of Ailanthus altissima var.erythrocarpa selected based on four indicators were collected from 7 provinces in China (Ningxia, Shanxi, Henan, Tianjin, Beijing, Hebei, and Shandong).The geographical locations of the sample spanned 106.00-117.72 east longitude, 33.76-41.12north latitude, and 34.7-1858.2m altitude, and the local soil was mostly red soil.When sampling, the latitude, longitude and altitude of the sample were recorded with a GPS device.The HN, TJ, BJ and SHD sampling sites were distributed discontinuously, with 6-10 plants per hectare.The number of individuals was small., the distribution was scattered, and the distance between the trees collected was not less than 50 m to prevent half-sib relationship between the samples.(Supplementary Table 1, Figs. 6 and 7).

DNA extraction
Fresh and tender Ailanthus altissima var.erythrocarpa leaves were selected and brought back into the laboratory in an ice box.Genomic DNA of Ailanthus altissima var.erythrocarpa was extracted by the improved CTAB method and an Omega kit.DNA concentration and quality were determined by a microspectrophotometer (model: Nano300).After quality determination, the DNA was diluted to 30 ng/μl in each tube.A total of 5 tubes (150 μl/tube) were used as working liquid.The original liquid and working liquid were stored in a freezer at − 20 °C for later use.

Development of SSR primers
Due to the lack of a reference genome sequence, leaves of tissue-cultured plantlets were sent to Wuhan Matwell Biotechnology Co., Ltd., for transcriptome sequencing and primer design.SSR primer design data were provided by Professor Dukejiu of Hebei Agricultural University.A total of 140 pairs of primers were randomly selected Redness Type Fruiting period Measurement

SSR primer screening and PCR amplification
First, the mixtures of 8 samples were used as templates to screen out primers with clear and bright bands that matched the size of the product.Second, primers with clear bands and good polymorphism were selected for amplification in 6 samples.Finally, the core primers were screened from 120 samples, which were verified with the primers.

Data analysis
The fluorescence capillary electrophoresis detection results were viewed using GeneMarker V2.4.0 software and collated in Excel.GenAlEx 6.502 software was used to calculate the average number of alleles (Na), average effective number of alleles (Ne), Shannon's information index (I), average observed heterozygosity (Ho), average expected heterozygosity (He), and fixation index (F).Hardy-Weinberg equilibrium tests were performed 18,51 .The polymorphism information content (PIC) was calculated using Cervus 307 software 18,51   www.nature.com/scientificreports/tree was drawn using MEGA7 52 .Mixed analysis based on Bayesian models was used for population structure analysis in Structure v2.3.4 53 .Then, the number of Markov chain Monte Carlo (MCMC) iterations after burn-in was set to 100,000, with a run length of 100,000 and a genetic uniform cluster number (K value) between 1 and 10.Each analysis consisted of 10 repeated runs.The optimal K value was determined using the best-K method in StructureHarvester 54 .Then used CLUMPP Windows.1.1.2bsoftware 55 to found the real K value, replaced the K value in STRU CTU RE, and generated a genetic structure map.

Statement
Permission to collect plant samples was obtained for this study.Voucher specimens were identified by Alfred Rehder and deposited in the Chinese Field Herbarium (CFH).This study complied with local and national guidelines. Vol

Figure 4 .
Figure 4. Fst and Nm values detected for population pairs of A. altissima var.erythrocarpa.Fst: above diagonal; Nm: below diagonal.

Figure 5 .
Figure 5. Structure analysis of 120 A. altissima var.erythrocarpa individuals.(a) Dendrogram based on the genetic distance of individuals in different populations (UPGMA); (b) principal coordinate analysis plot; (c) delta K distribution based on the rate of change in L(K) between consecutive K values.(d) Estimated average likelihood L(K) distribution(mean ± SD) from 2 to 10 possible clusters (K).(e) histogram of the structure assignment when K = 2-5.

Figure 6 .
Figure 6.Screening index system for 120 A. altissima var.erythrocarpa individuals.Redness: degree of redness of the samara; type: full red and gradually red; fruiting period: early red (June-July) and late red (July August); measurement: large fruit and small fruit.

Table 1 .
Types of SSR motifs in A. altissima var.erythrocarpa.

Table 3 .
Genetic diversity parameters for A. altissima var.erythrocarpa individuals at the nineteen microsatellite markers.

Table 4 .
F-statistics of different microsatellite markers in A. altissima var.erythrocarpa.

Table 5 .
Differentiation of 120 A. altissima var.erythrocarpa accessions by adding SSR marker primer combinations one by one.

Table 6 .
Genetic diversity parameters in seven natural populations of A. altissima var.erythrocarpa.
Used Populations v1.2.32 to calculate genetic distances, based on genetic similarity coefficients, a phylogenetic