Construction of SNP fingerprint and population genetic analysis of honeysuckle germplasm resources in China

Introduction The flower buds of Lonicera japonica Thunb. are widely used in Chinese medicine for their anti-inflammatory properties, and they have played an important role in the fight against SARS COVID-19 and other major epidemics. However, due to the lack of scientific and accurate variety identification methods and national unified standards, scattered and non-standardized management in flower bud production has led to mixed varieties that have caused significant difficulties in the cataloging and preservation of germplasm resources and the identification, promotion, and application of new L. japonica varieties. Methods In this study, we evaluated the population structure, genetic relationships, and genetic fingerprints of 39 germplasm resources of Lonicera in China using simplified genome sequencing technology. Results A total of 13,143,268 single nucleotide polymorphisms (SNPs) were identified. Thirty-nine samples of Lonicera were divided into four subgroups, and the population structure and genetic relationships among existing Lonicera germplasm resources were determined using principal component analysis, population structure analysis, and phylogenetic tree analysis. Through several stringent selection criteria, 15 additional streamlined, high-quality DNA fingerprints were filtered out of the validated 50 SNP loci and verified as being able to effectively identify the 39 Lonicera varieties. Discussion To our knowledge, this is the first comprehensive study measuring the diversity and population structure of a large collection of Lonicera varieties in China. These results have greatly broadened our understanding of the diversity, phylogeny, and population structure of Lonicera. The results may enhance the future analysis of genetic diversity, species identification, property rights disputes, and molecular breeding by providing a scientific basis and reference data for these efforts.


Introduction
Lonicera japonica Thunb. a perennial semi-evergreen twining shrub of the Caprifoliaceae family, is native to China and is primarily distributed in temperate regions of the northern hemisphere, with smaller numbers in Japan and Korea. The honeysuckle flower is a dry flower bud or open flower of L. japonica and is rich in organic acids, flavonoids , volatile oils, and other components. It can reduce fever and remove toxins without causing gastrointestinal issues (Li, 2020). Especially in China, honeysuckle flower has played an important role in the treatment of SARS COVID-19 and other major epidemics (Ren et al., 2022).
Before 2014, the National Pharmacopoeia Commission of China regarded Flos Lonicerae as a classification under L. japonica . The 2015 edition of the Chinese Pharmacopoeia detailed the differences between L. japonica and Flos Lonicerae regarding medicinal history, sources, characteristics, chemical components, and other aspects. The medicinal source of honeysuckle flower is L. japonica, and the medicinal sources of Flos Lonicerae include Lonicera macranthoides Hand.-Mazz., Lonicera hypoglauca Miq., Lonicera confusa DC., and Lonicera fulvotomentosa Hsu et S.C.Cheng. (Zhu et al., 2021). However, it is worth noting that compared with the three Flos Lonicerae, L. hypoglauca, L. confusa, and L. fulvotomentosa, L. macranthoides currently has more mature cultivation varieties with the largest number of plantings and the widest planting area, which is more recognized and accepted by the public . In the research on Flos Lonicerae, people usually choose L. macranthoides as the research object rather than L. hypoglauca, L. confuse, and L. fulvotomentosa (Yao et al., 2022).
With increasing market recognition of L. japonica and L. macranthoides germplasm and the continuous optimization and improvement of breeding technology, the number of Lonicera varieties is increasing . A total of 39 varieties of L. japonica and L. macranthoides were authorized or registered in China in April 2022 (Lan, 2017). Most of these varieties are propagated with cuttings and have strong adaptability with no strict soil or climate needs. Sandy loam with a thick soil layer is the best substrate, as it provides an extremely cold-resistant medium . However, due to the lack of standardized management regulations and decentralized management, the Lonicera germplasm resources are poorly organized, and there are unclear genetic relationships among varieties. These problems not only lead to intellectual property disputes between varieties but also make it very difficult to catalog and preserve the germplasm resources of Lonicera and cultivate new varieties for a wide range of applications. Therefore, accurate and efficient variety identification technology is needed for breeding new high-quality Lonicera varieties, and this makes variety identification increasingly important in the breeding industry. At present, the identification of Lonicera germplasm is limited to morphological and chemical fingerprint analyses. For example, methods s uch as high-performance liquid chromatography (HPLC) (Liu et al., 2019) and Fourier transform infrared spectroscopy (FTIR) (Nikzad-Langerodi et al., 2017) have been used to analyze Lonicera index components in most studies. DNA molecular marker technology is widely used in the identification of plant varieties due to its advantages of being free from environmental and species restrictions, its simple operation, and its ability to identify uniform and abundant loci as well as strong polymorphism. In plant biology, simple sequence repeat (SSR) and single nucleotide polymorphism (SNP) markers have been used to construct fingerprint databases, identify genetic relationships, construct genetic maps, and isolate and clone diverse genes of Nicotiana tabacum L. (Wang et al., 2021a), Brassica oleracea var. italica (Shen et al., 2021), Lagenaria siceraria (Mol.) Standl. (Wang et al., 2021b), Triticum aestivum L. (Li et al., 2018), Camellia sinensis L. (Karunarathna et al., 2021), Beta vulgaris L. (Schneider et al., 2007), and Ipomoea batatas (L.) Lam. , and other species. Shao et al. (2020) performed genetic detection on 113 olive germplasm resources using eight pairs of SSR fluorescence markers. Zhang et al. (2021) conducted molecular identification and DNA barcode construction of Dracaena germplasm resources from Liliaceae. (Wang et al., 2021a, b) constructed core SNPs to identify calabash germplasm resources using SNP markers, thereby providing data support for the molecular marker-assisted breeding of calabash varieties. SSR and SNP molecular markers have the advantages of rich marker polymorphism, good experimental repeatability, easy standardization of data, clear distribution of marker sites, and mature technology; they are the only recommended markers used to construct a DNA fingerprint database in the Biological and Molecular Marker Technology (BMT) Molecular Test Guidelines of the International Union for the Protection of New Plant Varieties (UPOV) and the General Principles of DNA Fingerprinting Methods for Plant Variety Identification in China (NY/T 2594-2016) (Hayward et al., 2015;Li et al., 2020).
Compared with SSR markers, SNP markers are abundant in number and display greater polymorphism; they are easily detected and convenient for statistical analysis, and they can be identified via high-throughput automated detection (Pei et al., 2015). With the continuous improvement of next-generation sequencing (NGS) technology, the development and detection of SNP sites have become simple and efficient. However, in the field of Lonicera variety identification, only Peng et al. (2010) designed a pair of allele-specific primers for PCR analysis of nrDNA ITS sequences and successfully distinguished the authenticity of five L. japonica species. Xu et al. (2015) used SSR technology to construct SSR fingerprints in Excel format for six different types of L. japonica, providing the earliest data support for variety identification. In addition, there have been no studies on the identification of Lonicera varieties based on DNA molecular markers (Jiang et al., 2013). The SNP molecular marker technology used in this study refers to the DNA sequence polymorphism caused by the variation of a single base in the genomic DNA sequence (Gazendam et al., 2022). In recent years, an increasing number of studies have used SNP molecular marker technology to identify varieties. DNA fingerprinting technology, with the advantages of being fast and accurate, is a powerful tool for identifying varieties and strains and has been widely used in the diversity and purity identification of many crops (Freixas-Coutin et al., 2019). Therefore, SNP molecular markers are regarded as the most important and promising DNA molecular markers in plant biological research (Azizi et al., 2021).
In this study, 39 Lonicera germplasm resources were used for Illumina NovaSeq sequencing. The results were compared with the reference genome to screen and mine SNP core markers and interpret their genetic relationships, genetic diversity, and population structure. In addition, a DNA fingerprint of Lonicera was constructed to efficiently and cheaply distinguish different species and varieties of Lonicera. The results provide a scientific basis and a data reference for genetic diversity analysis, variety identification, and molecular breeding of Lonicera (Supplementary Figure 1).

Plant materials
Considering that L. macranthoides in Flos Lonicerae has many cultivated varieties, a wide planting area, and more stable and representative genetic characters. In contrast, L. hypoglauca, L. confuse, and L. fulvotomentosa currently have no cultivated varieties, and there are few wild varieties with low representativeness and other factors. In this experiment, L. macranthoides was selected as the representative of Flos Lonicerae. The experiment was conducted with 39 wild and cultivated varieties of L. japonica and L. macranthoides existing in China. There were 35 varieties of L. japonica and four varieties of L. macranthoides ( Figure 1). The specific variety information is shown in Table 1. For each germplasm accession, seven vigorous plants were randomly selected at the seedling stage, and young leaf samples totaling 5 g were collected on tinfoil, frozen immediately in liquid nitrogen, and transported to the laboratory, where they were stored at −80°C until they were used for DNA extraction.

DNA extraction and library construction
After grinding the leaf samples, genomic DNA was extracted using a Plant Genomic DNA Kit (TIANGEN, China), and the quality and concentration of DNA were measured using a NanoDrop2000 UV spectrophotometer (Thermo Fisher, Waltham, MA, USA). Then, a paired-end library with a length range of 300-500 bp was constructed via double-digested (ddRAD) library construction of qualified sample DNA. First, 500 ng of genomic DNA was reacted with 0.6U EcoRI (NEB), T4 DNA ligase (NEB), ATP (NEB), and EcoRI connectors (including index sequences of differentiated samples) at 37°C for 3 h and annealed at 65°C for 1 h. The restriction enzyme NlaIII (NEB) and the NlaIII connector were then added and allowed to react for 3 h at 37°C . After the reaction, the endonuclease was inactivated at 65°C for 30 min in a polymerase chain reaction (PCR) amplifier. Digested products of 400-600 bp were recovered via agarose gel electrophoresis and quantified using Qubit 3.0 (Life Technology). After mixing 39 samples in equal quantities, an Illumina TruSeq kit was used to construct a DNA library of the mixed products (Hanania et al., 2004).

Simplified genome sequencing and reference genome alignment
An Illumina NovaSeq 6000 PE150 was used for sample sequencing after library construction, and data quality control was performed on the original sequenced reads (Dellaporta et al., 2016). Fastp software (version: 0.20.0) was used to remove reads with an unknown base number N <5, reads with a length of bases <50%, quality value <5, connector sequences, and other low-quality sequences to obtain clean data. The detailed parameters were set to -q5 -n5. Burrows-Wheeler Aligner (BWA, 0.7.17-R1188) was then used to compare the sequenced reads with the reference Simplified genome sequencing of 39 wild and cultivated varieties of L. japonica and L. macranthoides in China. Groups G1, G3, and G4 are L. japonica and Group G2 is L. macranthoides. Blue represents cultivar, orange represents wild, and brown represents crossbreed.  (Byers et al., 2012). The reference genome used was the GWHAAZE000,000,000 genome Fasta (L. Lonicera) (d ownlo ad web site: h t t p s : / / ng d c . c nc b . a c . c n/ s e a r c h / ? dbId=gwh&q=SAMC097356), and the parameters were set as -M -R. The insert size and coverage depth of each sample were counted, and variation was detected by comparing the positions of clean reads on the reference genome (Chao et al., 2009). Then, the same files generated by the comparison were converted to bam format using Samtools software (version: 1.9). Finally, Picard MarkDuplicates (version 2.21.2) was used to detect duplicate tags, and high-quality reads were retained for subsequent analysis.

Fingerprint construction
DNA fingerprint construction was performed based on the high-quality SNPs obtained (Xu et al., 2017). As few markers as possible were used to identify as many varieties as possible to achieve the purposes of simplicity, efficiency, and economy (Wang et al., 2009). Fifty core markers with a high detection rate and significant amounts of polymorphism that could distinguish all varieties were screened out based on the size of the polymorphic information content (PIC) value and distribution frequency (Pavana et al., 2012), and then a DNA fingerprint was constructed. Furthermore, the 15 leanest SNP combinations were further screened out (Jones and Mackay, 2015) to identify Lonicera varieties with lower costs and faster speeds.

Verification of SNP locus authenticity
Primers for the 15 SNP loci were designed using Primer Premier 5 software (Table 2), and the DNA of the validation test sample was extracted to conduct PCR. The parameters were set as follows: length, 18-30 bp; Tm value, 55-65 degrees; and GC content, 40%-70%. The primers were synthesized by Sangon Bioengineering Co., Ltd. (Shanghai, China). Four samples were randomly selected from 39 known samples, and three unknown samples were added (Table 1). A total of seven verification DNA samples were selected as templates for PCR amplification. The total volume of the PCR mixture was 25 ml, containing 1 ml DNA template (100 ng/ ml), 1 ml each of Primer F and Primer R (10 pmol/L), 0.2 ml Taq Plus DNA polymerase (5 U/ml), 2.5 ml 10× PCR buffer with Mg 2+ , and 1 ml dNTP (10 mM). The PCR reaction program involved predenaturation at 95°C for 5 min, denaturation at 94°C for 30 s, annealing for 30 s to 63°C, and extension at 72°C for a total of 30 s; this was repeated for 30 cycles, after which there was a final extension at 72°C for 10 min before storage at 4°C. After PCR amplification, 5 ml of PCR products were taken and subjected to 1% agarose gel electrophoresis at 150 V and 100 mA. After 10-20 min of observation, the target PCR band was cut from the gel, recovered, and then sequenced using 3730XL. The sequencing data were searched in the results group, and sequence analysis software SnapGene was used for analysis. Finally, SeqMan software was used to observe the peak map to test the reproducibility, genetic stability, and specificity of 15 SNPs.

Results
Simplified genome sequencing and reference genome alignment The sequencing of 39 Lonicera samples on the Illumina NovaSeq 6000 PE150 yielded a total of 84.88 Gb of clean data, and Q30 reached more than 93.25% per sample (Table 3). The obtained clean reads were mapped to the reference genome, and the mapping efficiency reached 99.38% (Table 4).  Read Number, the total number of paired-end reads in the clean data; Base Number, clean data; Q30 (%), the percentage of bases whose clean data quality value is greater than or equal to 30; Q20 (%), the percentage of bases whose clean data quality value is greater than or equal to 20; Average Q, average quality value. Total Reads refers to the total number of reads for sequencing, and a pair of paired-end reads refers to two counts of reads; Genome Size is the genome size; Cover Size is the size of covering genome; Cover Bases is the base number of the sequencing coverage; Mapped Reads is the number of reads in the map.

Selection and identification of high-quality SNPs
After sequencing, GATK software detected many SNP variants in 39 Lonicera samples; sequencing generated a total of 13,143,268 SNPs (Figure 2A). A total of 3,374,929 filtered SNPs were obtained, distributed on nine chromosomes ( Figure 2B). Finally, the SNP distribution map on each chromosome was drawn according to the number and density of SNPs ( Figure 2C). There were about 400,000 SNPs on chromosome GWHAAZE000000000 1 and 600,000 SNPs on chromosome GWHAAZE000000000 2. For chromosomes GWHAAZE000000000 3 to 9, the numbers of SNPs on these six chromosomes were similar, and there were about 300,000 SNPs on each chromosome.

Genetic relationships and population structure analysis
The PCA of high-quality SNPs screened from 39 Lonicera samples was conducted using GCTA software. L. macranthoides samples from Hunan were clearly clustered on one side ( Figure 3A and Table 1), while L. japonica samples from the same species were clustered on the other side ( Figure 3B and Table 1). For the two clusters classified by species, detailed clustering classification was further conducted according to a variety of characteristics of different Lonicera samples. For example, honeysuckle varieties with red flower buds were clustered together. Crossbred samples were completely clustered together. Admixture software was used to analyze the population structure of the 39 Lonicera samples. The cross-validation error rate was the lowest when K = 4, a result that was consistent with the variety characteristics and source classification of the 39 Lonicera samples ( Figure 3C and Table 1). The 39 Lonicera germplasm resources were divided into four subgroups (Figures 3D, E).
Finally, FastTree software was used to construct an evolutionary tree, and the 39 Lonicera samples were clustered into four groups (Figure 4), consistent with the results of the PCA and the population structure analysis. In addition, through further analysis, we found that the four clusters also had high similarity in phenotype, species, and origin. Each group was supported by a high bootstrap value.

Construction of a DNA fingerprint
Based on the results of simplified sequencing, SNP loci with a PIC value greater than 0.30 that were uniformly distributed on nine chromosomes were screened as the core loci for the construction of  Phylogenetic tree of 39 Lonicera varieties. G1 is a L. japonica variety with red and light red flower buds; outgroup G2 is L. macranthoides; G3 is mainly wild L. japonica from Shandong and Henan; and G4 is a mature and highly recognized L. japonica in the current market. The text annotation in the figure is: species name (species type, origin). Li et al. 10.3389/fpls.2023.1080691 Frontiers in Plant Science frontiersin.org different loci between each pair. To quickly and cheaply distinguish Lonicera varieties, in this study we screened out the combination with the minimum number of SNPs that could distinguish 39 Lonicera varieties from the obtained 50 core loci; we obtained the most simplified SNP combination that contained 15 SNPs distributed on eight chromosomes ( Figure 5B). When using the combination of 15 SNPs to distinguish samples, there was at least one different SNP between each of the two samples that could effectively distinguish the 39 Lonicera varieties. Some Lonicera samples could only be distinguished based on one SNP. For example, "Longhua" (SPL022), "Yincuilei" (SPL019), and "Jincuilei" (SPL020) could only be distinguished by Marker 7; "Yujin 2" (SPL010) and "Yate 1" (SPL004) could only be distinguished by Marker 6. "Yujin 5" (SPL014) and "Yujin 5 1-2" (SPL028) could only be distinguished by Marker 11. These results also illustrate that it is difficult to develop fingerprints against the background of complex and chaotic genetic relationships in L. japonica germplasm.

Verification of SNP loci and genetic stability of honeysuckle
We randomly selected four samples, "Yujin 1," "Yujin 2," "Mixian Damaohua," and "Telei 1," from a set of 39 samples to verify the reproducibility of SNPs. Fifteen SNPs had the same site specificity as the previous test results (Table 5), with a clear and clutter-free locus peak map, 100% reproducibility, and ideal results. Next, we tested the genetic stability of 2-, 5-, 7-, 8-, 9-, and 10-year-old plants of "Yujin No. 1." The results showed that the genetic stability of honeysuckle of the same variety in different years was as high as 97.33% (Table 6). Finally, three new varieties, "Yateliben," "Bainong 2," and "Weizi," were added to verify the specificity of 15 SNP loci (Table 7). The results showed that 15 SNP loci could effectively distinguish the genotypes of the newly added varieties (Supplementary Figure 2).

SNP-based genetic relationships among the Lonicera
As the original and major producer of Lonicera, China has three main production areas: Fengqiu in Henan Province, Julu in Hebei Province, and Pingyi in Shandong Province (Song, 2020). The planting area of L. japonica in 2022 reached 1,066.7 km 2 (Lv, 2020;Ma et al., 2022), with an annual output value of nearly 10 billion CNY (Liu, 2021;Liu, 2022). At present, many Lonicera varieties have been widely planted, but there has been a lack of standardized and decentralized management. Therefore, it was necessary to clarify the genetic relationship and population structure of existing Lonicera varieties.
The 2015 Edition of the Chinese Pharmacopoeia details the differences between L. japonica and Flos Lonicerae in medicinal history, sources, characteristics, chemical components, and other aspects, showing that Flos Lonicerae is classified as L. macranthoides, L. hypoglauca, L. confusa, and L. fulvotomentosa. L. macranthoides is more recognized because of its easy cultivation, large number of plants, and wide planting area. Lu et al. (2017) established the fingerprints of five batches of L. japonica and three batches of L. macranthoides using proton nuclear magnetic resonance (1H-NMR) that provided a basis for the quality control analysis of L. japonica and L. macranthoides. Ren et al. (2022) used microscopy and thin-layer chromatography (TLC) to evaluate the

Marker6 T/T T/T T/T T/T T/T T/T T/T T/T T/T T/T
Marker is the name of the SNP site; REF is the reference genome genotype; ALT is the verification result. quality of L. japonica varieties in Henan, Hebei, Shandong, and other areas and concluded that the index components of Shandong L. japonica were the highest. This work provided a scientific basis for the selection of raw materials for honeysuckle decoctions. Deng et al. (2021) conducted a comparative study on tree-shaped and wild-type L. japonica using HPLC fingerprints and found that wild L. japonica contained more chemical components, while treeshaped L. japonica was more consistent. Li et al. (2020) conducted a comparative analysis of the characteristics and pharmacological effects of L. japonica and L. macranthoides, providing data support for their application in traditional Chinese medicine. The related studies focused on Lonicera in the early stages; most of the studies on the germplasm focused on morphological analysis, index component differences, and quality evaluation among varieties, while there was almost no research on the genetic relationships and population structure of Lonicera species and varieties. In this study, 39 existing Lonicera germplasm resources in China were collected and divided into four subgroups via simplified sequencing, PCA, population structure analysis, and evolutionary tree construction. Subgroup I contains the red variety L. acuminata of L. japonica. Subgroup II contained L. macranthoides. Subgroup III comprises the main, wild, and hybrid varieties of L. japonica from Shandong and Henan. Subgroup IV contained cultivars of L. japonica that comprised varieties with high yields, strong resistance to pests and diseases, and good characteristics that have been recognized on the current market, for example, SPL009, SPL002, SPL003, and other varieties. Based on these results, it should be clear that the lineages of wild Lonicera species were extremely complex, while the lineages of systematically bred L. japonica were relatively simple. For example, "Mixian xianhua" (SPL006), "Mixian wild" (SPL007), and "Mixian Damaohua" (SPL039) are wild varieties of L. japonica with extremely complex genetic relationships. In the population structure analysis, they all had ancestral lines from the four subgroups. However, the genetic relationship between "Yujin 1" (SPL009) and "Yate liangzhong" (SPL035) was simple and derived from one pedigree. After further analysis, this study also clarified the controversial genetic relationship of "Longyao" (SPL023), a variety produced in Hunan, as belonging to L. macranthoides. However, our results confirmed that "Longyao" (SPL023) was a cross-species hybrid of a female parent L. japonica and a male parent L. macranthoides, which explained why it simultaneously had the pedigree of two subgroups in the PCA and population structure analysis. The genetic relationship of "Yujin 6" (SPL027) was further clarified, and the results showed that it was crossbred from the female parent "Telei 1" (SPL011) and the male parent "Yujin 1" (SPL009). This finding well explained its position in Group IV of the evolutionary tree. In addition, four other hybrids were used in this study. Population structure analysis showed that all hybrids have genetic genes from their parents, which was largely consistent with expectations. At present, this research group is conducting a more in-depth analysis of hybrid characteristics.    , the present study focused on aspects including analyzing the population structure of L. japonica, excavating the differences among varieties of L. japonica, and using SNP to distinguish varieties of L. japonica. This work can not only explain the population relationships of existing Lonicera germplasm resources but can also further clarify the genetic relationships between Lonicera varieties. These results provide a favorable reference and a basis for the management of Lonicera germplasm resources and have important significance for the breeding of new Lonicera varieties in the future.
In this study, SNP molecular marker technology and DNA fingerprinting technology were used; 3,374,929 SNPs were obtained from 39 samples of Lonicera; 50 core SNPs with a strong ability to identify Lonicera germplasm were screened; and DNA fingerprints were constructed. A fingerprint map constructed with 50 core SNP markers effectively identified and distinguished differences among different species and varieties of Lonicera. In addition, this study further filtered out 15 SNPs as the most compact site combination. The DNA fingerprints constructed based on the 15 SNP loci were used to detect different Lonicera samples, and ≥1 different locus was found in each Lonicera variety that could be used to distinguish different species and varieties rapidly, accurately, and at a low cost.
Considering the extremely complex genetic background of Lonicera germplasm, this study conducted further research on varieties that had several different loci ≤1 in the fingerprints of 15 SNP loci. After further analysis, it was found that these varieties not only had high repeatability in DNA fingerprints but also had high similarity in origin, species, and characteristics. For example, "Mihua 1" (SPL017) and "Fengjin 1" (SPL012) were introduced and domesticated from "Mixian Damaohua" (SPL008). The varieties "Longhua" (SPL022), "Jincuilei" (SPL020), and "Yincuilei" (SPL019) originated from Hunan Province and were L. macranthoides "Yujin 2" (SPL010) and "Yate 1" (SPL004) had purplish red flower buds and were red varieties of L. japonica. "Yujin 5" (SPL014) and "Yujin 5 1-2" (SPL028) were both hybrids, and their parents were the same varieties (Table 1). Therefore, it was concluded that when the difference in the DNA fingerprint between two varieties is ≤1 locus, it means that the two varieties have certain similarity in species, origin, parents, flower color, etc., which is very likely to be two identical or similar varieties.

Reliability and stability of Lonicera SNP fingerprints
To verify the reliability of the DNA fingerprints comprising 15 SNP loci, four samples were randomly selected from the original 39 samples, and 15 SNP loci were analyzed via PCR and Sanger sequencing. The 15 SNP verification results were completely consistent with the previous test results. Subsequently, genetic analysis was conducted on 10 honeysuckle samples of the same variety from different years. The results showed that the genetic stability of these honeysuckle samples was extremely strong, reaching 97.33%, a level that was consistent with the characteristics of honeysuckle as a vegetatively propagated plant. However, it is worth noting that the genotypes of the two-year-old samples "J01," "YJ02," and "YJ03" and the five-year-old sample "YJ07" changed from A/A to A/G in Marker 10, a result that may have been related to the second allele of the SNP. However, as a vegetatively propagated plant, the genetic traits of honeysuckle are very stable, and the probability of genetic mutation between the same varieties is very low; thus, this phenomenon was most likely caused by the duality of SNP sites. Finally, three new varieties verified the specificity of SNPs, further demonstrating the reliability of the SNPs and DNA fingerprints developed in this study.
In general, the verification results of the 15 SNP loci and the genetic stability of honeysuckle were in line with expectations; the analysis not only confirmed that the 15 SNP loci constituted the authenticity and reliability of DNA fingerprints but also verified the high genetic stability of honeysuckle, providing a data reference and theoretical support for the identification of honeysuckle germplasm and the SNP fingerprint identification method.

Data availability statement
The original contributions presented in the study are publicly available. This data can be found here: NCBI, PRJNA895242.

Author contributions
JL conceived of the project. JL and XC designed the experiments. XC performed most of the experiments. QH assisted in preparing seedlings for experiments. PL and XZ helped overcome technical difficulties in SNP site validation. CC, XZ, FL, and YW assisted in analyzing the data. JL, XC, and CC interpreted the data generated in the figures and wrote the manuscript. All authors contributed to the article and approved the submitted version.

Funding
This work was supported by the Science and Technology Major Project of Xinxiang City (ZD2020002),the 2022 Provincial Financial Forest Germplasm Resources Construction Project Honeysuckle Germplasm Resources Construction, the Key project at central government level: The ability establishment of sustainable use for valuable Chinese medicine resources (2060302), and the Postgraduate Education Reform and Quality Improvement Project of Henan Province (YJS2023JD33).