Next Article in Journal
Ecological Features and Conservation of Urtica rupestris Guss. (Urticaceae): A Narrow Endemic Species of Sicily
Next Article in Special Issue
Preliminary Investigation of Essentially Derived Variety of Tea Tree and Development of SNP Markers
Previous Article in Journal
Kukoamine B from Lycii Radicis Cortex Protects Human Keratinocyte HaCaT Cells through Covalent Modification by Trans-2-Nonenal
Previous Article in Special Issue
Integrative Transcriptomic and Phytohormonal Analyses Provide Insights into the Cold Injury Recovery Mechanisms of Tea Leaves
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Plants
Volume 12
Issue 1
10.3390/plants12010162
plants-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Development of SNP Markers for Original Analysis and Germplasm Identification in Camellia sinensis

by
Liubin Wang
1,2,
Hanshuo Xun
1,
Shirin Aktar
1,
Rui Zhang
1,
Liyun Wu
1,
Dejiang Ni
2,*,
Kang Wei
1,* and
Liyuan Wang
1,*
1
Key Laboratory of Tea Biology and Resources Utilization, Ministry of Agriculture, Nature Center for Tea Improvement, Tea Research Institute Chinese Academy of Agricultural Sciences (TRICAAS), Hangzhou 310008, China
2
College of Horticulture and Forestry Science, Huazhong Agricultural University, Wuhan 430070, China
*
Authors to whom correspondence should be addressed.
Plants 2023, 12(1), 162; https://doi.org/10.3390/plants12010162
Submission received: 8 December 2022 / Revised: 26 December 2022 / Accepted: 27 December 2022 / Published: 29 December 2022
(This article belongs to the Special Issue Breeding and Genetic Mechanism of Tea Plants)
Browse Figures
Review Reports Versions Notes

Abstract

:
Tea plants are widely grown all over the world because they are an important economic crop. The purity and authenticity of tea varieties are frequent problems in the conservation and promotion of germplasm resources in recent years, which has brought considerable inconvenience and uncertainty to the selection of parental lines for breeding and the research and cultivation of superior varieties. However, the development of core SNP markers can quickly and accurately identify the germplasm, which plays an important role in germplasm identification and the genetic relationship analysis of tea plants. In this study, based on 179,970 SNP loci from the whole genome of the tea plant, all of 142 cultivars were clearly divided into three groups: Assam type (CSA), Chinese type (CSS), and transitional type. Most CSA cultivars are from Yunnan Province, which confirms that Yunnan Province is the primary center of CSA origin and domestication. Most CSS cultivars are distributed in east China; therefore, we deduced that east China (mainly Zhejiang and Fujian provinces) is most likely the area of origin and domestication of CSS. Moreover, 45 core markers were screened using strict criteria to 179,970 SNP loci, and we analyzed 117 well-Known tea cultivars in China with 45 core SNP markers. The results were as follows: (1) In total, 117 tea cultivars were distinguished by eight markers, which were selected to construct the DNA fingerprint, and the remaining markers were used as standby markers for germplasm identification. (2) Ten pairs of parent and offspring relationships were confirmed or identified, and among them, seven pairs were well-established pedigree relationships; the other three pairs were newly identified. In this study, the east of China (mainly Zhejiang and Fujian provinces) is most likely the area of origin and domestication of CSS. The 45 core SNP markers were developed, which provide a scientific basis at the molecular level to identify the superior tea germplasm, undertake genetic relationship analysis, and benefit subsequent breeding work.

1. Introduction

Tea is one of the most popular nonalcoholic beverages in the world. The popularity of tea is associated with its excellent taste and flavor as well as its many health benefits, such as reducing the risks of cancer, cardiovascular disease, and incidence of kidney stones; protecting teeth and bones; and its anti-obesity effect [1,2,3]. Currently, the tea plant is an important economic crop that is extensively cultivated in more than 60 countries [4]. The main tea-cultivating countries, for instance, China, India, Japan, and Kenya, have undertaken a major task involving the continuous genetic improvement of tea plants. More than 1200 tea cultivars with special traits have been developed and released for cultivation worldwide [5]. Based on this situation, a simple and accurate method is needed for the variety identification, kinship, and phylogenetic analysis of tea plants. This is very important to assure the fidelity and purity of varieties when cultivars spread and to protect the rights of breeders and consumers. For breeders, the germplasm source and pedigree information can be used to construct genetic linkage maps and to improve breeding efficiency in tea plants [6].
The tea plant (Camellia sinensis ) is diploid (2n = 30 chromosomes), with a genome size of about 3G [7,8]. It is a self-incompatible plant species, which greatly contributes to its high genetic diversity. In recent decades, hundreds of new tea cultivars have been bred. Compared to old varieties, many new cultivars show excellent agronomic traits, such as high production and improved stress resistance and quality. Hence, most new tea gardens have been recently created using improved clonal cultivars. Cultivated tea plants can mainly be divided into two varieties: Camellia sinensis var. assamica (CSA, Assam type) and Camellia sinensis var. sinensis (CSS, Chinese type). CSA is mainly cultivated in warm tropical areas and is processed into black tea or puer tea, whereas CSS is cultivated in a broader geographical location, mainly in high-latitude areas, for green tea production. Recent studies have shown that the divergence between CSS and CSA occurred about 0.38 to 1.54 million years ago [7]. However, their genetic diversity and differentiation, especially the original relationship between CSS and CSA, still remain a mystery.
Molecular marker technologies such as RAPD and SSR have been widely used for studying tea plants [9,10,11,12]. However, with advanced bioinformation technology, traditional molecular markers no longer meet the needs of researchers. Single-nucleotide polymorphism (SNP) markers as an efficient tool have become the main choice for this work. Compared with traditional RAPD and SSR molecular markers, SNP markers have the following advantages: (1) abundant polymorphism in the plant genome [13], such as in the rice genome, where the average density of SNP markers has been reported to be about one SNP per 232 bp [14]; (2) dimorphism, as SNP markers generally consist of only two alleles compared with SSR markers and distinguish polymorphisms by fragment size and can reduce the detection error rate [15,16]; and (3) high-throughput technology eliminates traditional molecular labeling and the tedious gel electrophoresis typing procedure and improves the detection rate [17,18]. In other plants, such as tomato and litchi, a small number of SNPs have been screened for genetic diversity analysis and population structure analysis [19]. Recently, high-throughput sequencing technology and SNP array have been applied to SNP genotyping to analyze the genetic relationships and population structure [16,20]. However, it is not practical to conduct high-throughput sequencing on a large number of individuals; therefore, it is particularly important to select a small number of core SNP markers for tea plant research.
In our previous study [20], we reported on the development of the first whole tea genome 200 K SNP array with 179,970 unique and informative SNPs. At the same time, the SNP information was used to analyze the population structure and phylogenetic tree of 142 tea cultivars from different provinces in China. In this study, we performed a more in-depth analysis of the population structure and phylogenetic tree of tea cultivars, and 45 core SNP loci were carefully selected from 117 leading Chinese clonal tea cultivars. The objectives of this study were (1) to reveal the original relationship between CSS and CSA cultivars and (2) to characterize a set of suitable SNP marker tools, i.e., a standard fingerprint identification map and parentage verification.

2. Results

2.1. Diversification and Geographical Distribution of Tea Cultivars

In total, 142 tea cultivars from different provinces in China were used to perform population structure and phylogenetic analyses in our previous report. The results of the two analyses were surprisingly consistent; i.e., all cultivars were clearly divided into three groups, namely CSA, CSS, and the transitional type. (The detailed information is shown in the Supplementary Materials Table S1.) In this study, we further analyzed the geographic distribution of the CSA, CSS, and transitional type varieties (Figure 1). Most CSA cultivars are found in south and west China, mainly Yunnan Province, which confirms previous reports that Yunnan Province is the primary center of CSA origin and domestication [21,22], while most CSS cultivars are distributed in east China, with a clear bias towards Zhejiang and Fujian Provinces. Our data support the concept of the polyphyletic domestication of tea plants and indicate that east China (mainly Zhejiang and Fujian Provinces) is the center of CSS domestication.

2.2. Polymorphism and Location Distribution of 4115 SNP Loci

We selected 4115 high-quality SNPs of “poly high resolution”, “no minor homozygote”, “mono high resolution”, and “hemizygous” categories from 179,970 SNP loci. In order to visually display the position of 4115 SNP loci on chromosomes, the SNPs were mapped to chromosomes of “Shuchazao” [23]. The number of markers ranged from 221 to 401 on each chromosome, and the result is shown in Figure 2a. The results of the polymorphism analysis of 4115 markers are shown in Table S2. The minor allele frequency (MAF) ranged from 0.018 to 0.50, with an average of 0.33, and more than 99.73% of SNP loci had values greater than 0.05 (Figure 2b). The polymorphism information content (PIC) was between 0.03 to 0.38, with an average of 0.33, and more than 88.58% of SNP loci showed high diversity, with PIC exceeding 0.250 (Figure 2c). These data show that the loci are evenly distributed across the tea plant genome, with very few missing loci and high diversity, providing a solid foundation for the subsequent screening of core markers.

2.3. Core SNP Loci Selection and Analyses of Polymorphism

To obtain a more convenient SNP genotyping method for tea plants in the future, we further screened 4415 SNP loci with strict screening criteria: no missing allelic data, MAF, and PIC values. In order to obtain the SNP loci by PCR amplification using the designed primers directly, 35 bp up- and downstream sequences of the SNP sites were strictly conserved within the genome of “Shuchazao”. (Primers can be designed by referring to the genome to expand the sequence in the up- and downstream regions of the markers.) Based on the criteria, all 45 SNP loci were finally selected from three SNP loci on each of the 15 chromosomes (Table 1), and the up- and downstream sequences of the 45 SNP markers are available in Table S3.
Genetic diversity analyses were performed on the 117 tea cultivars (Table S4) using 45 SNP loci (data shown in Table 2). The MAF were in the range of 0.475 to 0.500, the mean was 0.447, and the observed number of alleles for each SNP locus was 2. The effective number (Ne) of alleles ranged from 1.93 to 2, and the mean was 1.989. Shannon’s information index (I) was between 0.675 to 0.693, with an average of 0.690. The measured heterozygosity (Ho) ranged from 0.299 to 0.590, with an average of 0.464. The expected degree of heterozygosity (He) was in the range of 0.482 to 0.502, with an average of 0.499. The PIC was between 0.366 to 0.375, with an average of 0.374. These results indicated that these core markers exhibited high polymorphism.

2.4. Development and Validation Marker Tools for the Identification of Germplasm and Parentage Verification

2.4.1. Identification of Germplasm and Construction of DNA Fingerprints

The probability of identity (PI) is often used as a parameter for identifying the ability of markers to distinguish breeds [10,11]. For the screening of SNP loci, each SNP locus was located on different chromosomes; the nearest locus was ~6.3 Mb (2_SNP_15 and 3_SNP_15), and all sites were considered to be independently isolated. Based on the genotyping data of the 117 tea cultivars at the 45 SNP loci, PI values were calculated in GenAlex 6.5, and the results are shown in Table 2. For the SNP loci, the PI ranged from 0.375 to 0.384, with an average of 0.376. Moreover, PI values for increasing the number of locus combinations of the 45 loci were determined (Figure S1). The combined PI value for the eight SNP markers was 3.93 × 10−4, which indicated that 117 tea cultivars could be distinguished by using eight markers. Eight low-PI-value markers (1_SNP_2, 1_SNP_1, 3_SNP_14, 3_SNP_15, 1_SNP_13, 2_SNP_14, 1_SNP_5, and 2_SNP_4) were selected to construct the fingerprint, and all of them had unique genotypes (Figure 3). Therefore, these eight markers can be used as primary markers to identify the germplasm of tea plants. However, considering that there are at least thousands of tea plant resources [24], the remaining markers can be used as standby markers for germplasm identification. The detailed genotypes of the 117 tea cultivars based on the 45 loci are shown in Table S5, which can provide a reference for the identification and comparison of tea germplasm resources.

2.4.2. Parentage Verification

A multi-locus probability of exclusion (PE) was calculated for the 45 markers based on their SNP type as a measurement of the parentage verification performance. Three PE types were created (Figure 4): multi-locus combination types with the probability of exclusion when the genotypes of both parents are known (termed Q1) or when the genotype of only one parent is known (termed Q2) and two putative parents unknown (termed Q3). With an increase in the number of marker combinations, the PE values also increased significantly. When the number of randomly combined markers reached 40, the PE values of the three types were 0.9997, 0.9949, and 1.0, respectively. The parentage of the 117 tea cultivars was completely confirmed by the core markers. The results suggest that the markers selected in our analysis are informative and suitable for parentage verification.
The parent–child pairs among the 117 tea cultivars were calculated using Cervus 3.0 [25]. We consulted the literature to understand the existing pedigree relationships in all tea cultivars and then verified their parents and offspring using the maximum likelihood method with 95% confidence level. Methods to infer both paternal and maternal half sibship as well as full sibship in a sample of offspring and to jointly infer the parentage of the offspring when candidate father and mother samples are also available are shown in [26].
The most likely 10 pairs of parents and offspring relationships are listed in Table 3. Interestingly, we found all of the pairs of related varieties that were perfectly matched at the 45 SNP loci. Huangmeigui, Chungui, and Mingke 1 were offspring of Huangdan, and Longjing 43 and Changye Baihao were the parents of Zhongcha 108 and Foxiang 3, respectively. Our data results are consistent with the breeding records in [27]. The existing pedigree information of Mingke1 (Tieguanyin♀ × Huangdan♂) [27] and Zhongming 7 (Zhongcha 108♀ × Longjing 43♂) [28] was very clear when Tieguanyin and Zhongcha 108 were set as the mothers of Mingke 1 and Zhongming 7, respectively. A search was conducted for the candidate fathers from the 117 tea cultivars. The results showed that Huangdan was identified as the paternal parent of Mingke1. However, our data did not indicate that Longjing 43 was the paternal parent of Zhongming 7. In addition, we found that Tieguanyin was the parent of Benshan and Xingrencha, Longjing 43 was the parent of Shuchazao, and Yungui was the parent of Aifeng.

2.4.3. Analysis of the Phylogenetic Tree

To verify the ability of our core markers in the cluster analysis, a phylogenetic tree of the 117 tea cultivars was constructed based on the genetic distance. First, the groups (CSS, transitional type, and CSA) to which each cultivar belonged were clearly divided [20] and are shown together with the clustering results in Figure 5. Interestingly, individuals belonging to the CSA group were clearly differentiated from the other two groups, with only a few scattered between them. However, while some individuals were distinct between two groups, in general, CSS and the transitional type were generally not clear.

3. Discussion

3.1. New Insights into the Origin and Domestication of Chinese Teas

According to available documentation, China is the first country to cultivate and utilize tea plants [29]. Although it has been recognized that China is the origin of tea plants [30], the origin of domestication of tea plants in China has long been unilateral or controversial. A previous hypothesis stated that cultivated tea plants originated in southwest China (mainly Yunnan Province) and then spread to other parts of China and the world [12,21,31]. This hypothesis was supported by the evidence that (1) the genetic diversity of tea plants shows the highest value in original regions such as Yunnan, with a decreasing trend in regions farther away from the origin center [12]; and (2) the ancient tea plants from Yunnan Province are located at the base of the phylogenetic tree according to the resequencing data of 81 tea accessions, and strong gene flow occurs from the southwest provinces to other parts of China [23]. Meegahakumbura et al. (2016, 2018) proposed that the CSA and CSS cultivars in China were independently domesticated based on nuclear DNA microsatellite analysis [21,22]. However, the exact domestication center of CSS cultivars has not yet been confirmed. A key limitation is the lack of sufficient molecular markers for diversity evaluation [30]. The development of a 200 K SNP array largely overcomes this problem [20].
In this research, we provide rich data to support inferences about the evolutionary and domestication history of tea plants. Our data support the concept of polyphyletic domestication and indicate that southwest China (mainly Yunnan Province) is the domestication center of CSA cultivars (Figure 1). Moreover, our results also support the new idea that the formation of CSS cultivars has not likely been driven by human intervention and that east China (mainly Zhejiang and Fujian Provinces) is likely the domestication center of CSS cultivars. The genomic analyses revealed that CSA and CSS cultivars diverged around 0.38 to 1.54 million years ago, raising the suggestion that there may be a long predomestication period [7]. Furthermore, recent archaeological evidence has also shown that many well-preserved tea roots have been identified at the Tianluoshan site of Yuyao, Zhejiang, which suggests CSS tea plants were cultivated in Zhejiang about 6000 years ago [32]. Due to the environment and transportation mode at that time, it was unlikely that such artificial migration of tea plants from Yunnan to Zhejiang would be realized. Therefore, we deduce that east China (mainly Zhejiang and Fujian Provinces) is most likely the area of origin and domestication of CSS tea plants. Moreover, tea plants can clearly be divided into CSA and CSS cultivars as well as separate centers of domestication origin, which is also evidence that tea plants may have undergone adaptive evolution in contrast to other crops [33] (e.g., rice) that were subjected to strong domestication pressures.

3.2. Development of SNP Marker Tools: Identification of Germplasm, Parentage Verification, and Phylogenetic Analysis

Genetic diversity is the basis for the research of tea germplasm resources and genetic breeding. The research on the genetic diversity of tea plants has changed from morphological and biochemical markers to molecular markers [16]. Previously, the commonly used molecular markers were mainly SSR markers. However, compared with SSR markers, SNP markers have more obvious advantages as a new generation of molecular marker technology. Generally, SNPs are composed of two bases with low mutation frequency and high genetic stability; high-throughput data are more accurate, and the cost is low [18]. In other crops, such as tomato and cotton, core SNP molecular markers have been developed for genetic diversity and germplasm identification [34,35]. However, there are few reports on tea plants. Therefore, there is a need to develop high-quality core SNP markers for genetic diversity research of tea plants.
In this study, with the convenience of a practical application and a reduction in application costs, based on the 200 K gene chip technology developed, for the first time in tea plants, 45 high-quality SNP loci were screened from 177,790 high-quality SNP loci covering the whole genome of tea plants evenly with the following filter criteria: MAF, PIC, and the position on the chromosome. The polymorphism of 45 markers (A = 1.989, I = 0.690, Ho = 0.464, He = 0.499, and PIC = 0.374) is higher than the 60 SNP markers in a previous report (I = 0.512, Ho = 0.401, He = 0.341) [19] and higher than 86 SNP markers (I = 0.517, Ho = 0.370, He = 0.346) in [36]. Finally, DNA fingerprint construction, parentage verification, and the phylogenetic analysis of 117 tea varieties that are widely cultivated in China were performed. The results show that these markers are highly polymorphic and have strong fingerprinting power and kinship analysis capabilities. These results provide important reference data for genetic diversity analysis, germplasm identification, and genetic relationship analysis of tea plants

3.2.1. Identification of Germplasm

Our SNP marker tools can effectively solve the problem of identifying tea germplasm resources. With the increase in the number of tea plant varieties, it is very important to accurately identify the uniqueness of each variety and to ensure the authenticity and purity of germplasm in the promotion of planting and scientific research. PI is an important indicator to evaluate the ability of molecular markers to discriminate germplasm [19]. In this study, the PI values for discriminative combinations were 8.05 × 10−20 for forty-five SNP markers and 3.93 × 10−3 for eight (Figure S1). These results suggest that these markers can accurately distinguish the 117 accessions. Since SNPs are dimorphic, and theoretically, 256 varieties can be distinguished by eight markers (28 = 256), there is still room for improvement for germplasm identification. Moreover, taking into consideration that there are currently thousands of tea accessions worldwide [5], there is great potential for applications using these 45 SNP markers. These 45 SNP markers can be used to establish a standard fingerprint for each germplasm in the future. When an unknown germplasm is obtained, the SNP site information can be detected by biotechnology, which can be directly compared with the standard fingerprint, and then, the germplasm information can be accurately identified. In order to provide a reference for the construction of a standard fingerprint of a tea plant tree, we have published detailed information (Table S4), such as the up- and downstream DNA sequences of the 45 loci, their positions on chromosomes, and the SNP genotypes of the 117 tea varieties.

3.2.2. Parentage Verification

SNP marker tools can provide a basis for determining the genetic relationship of tea plants. Correct pedigree relationships are helpful for linkage map construction, core germplasm conservation, and breeding [10]. PE is a favorable parameter for evaluating kinship identification [37]. In this study, three types of PE (Q1, Q2, and Q3) for increasing combinations of the 45 markers had values of 100%, which showed that the parentage of the 117 tea cultivars could be accurately verified.
Seven cultivars that showed clear parentage to each other were also confirmed in this study (Table 3). For example, Zhongcha 108 was selected from the progeny of Longjing 43 that had suffered from radiation exposure [27]. Shuchazao was more likely to be an offspring of Longjing 43, and Shuchazao was bred from a local population in Anhui Province. Since Anhui and Zhejiang are adjacent, and breed exchanges are frequent, the relationship between the cultivars may be true. Our results are similar to those reported by Tan et al. (2015) [10]. Tieguanyin is an elite cultivar for making oolong tea; it has been cultivated for more than 200 years and is widely cultivated in Fujian Province [27], and it is frequently used as a crossing parent in hybridization breeding programs. Therefore, Benshan and Xingrencha were most likely the offspring of Tieguanyin through natural hybridization programs. Aifeng and Yungui were selected from the same local population at Puwen farm [27], and due to frequent natural interbreeding, they may be related.

3.2.3. Phylogenetic Analysis

A sufficient number of molecular markers is an important basis for analyzing population diversity [30]. Previously, we used high-quality molecular markers to accurately distinguish tea plants [20]. The three varieties that originated from Zhejiang were classified into groups as the transitional type (S1). We examined the existing pedigree information [27] and found that Hanlv, Zhenong 12, and Zhenong 21 were a cross between the CSS and CSA groups, which indicates our previous cluster analysis was relatively accurate. Here, the core SNP markers we screened can provide some reference for the differentiation between CSA and other groups (CSS and transitional type), but they are not ideal in the differentiation between CSS and transitional groups. Although these markers were selected using strict criteria, a comparison with previous accurate classifications was difficult to achieve. Therefore, for the phylogenetic tree analysis, more markers should be selected as much as possible. Moreover, the results showed that many varieties were selected between CSA and CSS through the artificial hybrid, which reminds that we should pay attention to the conservation and collection of tea plant resources to improve the diversity of germplasm.

4. Materials and Methods

4.1. Plant Materials

A total of 142 cultivars (lines) that originated from the main tea growing provinces (Anhui, Fujian, Guangdong, Guangxi, Guizhou, Hubei, Hunan, Jiangxi, Sichuan, Taiwan, Yunnan, and Zhejiang) in China were examined. Young shoots of those cultivars (lines) were collected from the Songyang tea germplasm resource nursery in Zhejiang and the puer tea breeding center in Yunnan, China. Detailed information of the origin and the main characteristics of each cultivar (line) are available in the Supplementary Materials Table S1.

4.2. SNP Detection and Selection from the Genome

The genome re-sequencing data of 139 tea accessions generated in previous study (Wang et al., 2020) were used for SNP discovery [38]. The methods of SNP calling and filtering were also described in previous study [38]. More than 218 million high-quality variants from the re-sequencing analysis were identified and used to design the 200 K AXIOM SNP array. The criteria used for the initial SNP quality analysis of the array included the following: one SNP marker for every 500 bp to ensure uniform distribution throughout the tea genome and the locations of SNPs on or close to genes. Finally, a set of 201,536 quality-filtered SNPs from 15 chromosomes of tea plants was selected and sent to the Affymetrix Bioinformatics Services (Santa Clara, CA, USA) for the final selection and design of the array. The quality of each SNP was further assessed using an in silico validation with proprietary software. The final 200 K AXIOM SNP array contained a total of 179,970 SNP markers. A panel of 142 diverse tea cultivars was genotyped for the initial validation of the 200 K AXIOM SNP array. The National Center for Biotechnology Information has received related information and assigned accession number GSE182082 to it [20]. Then, based on the 179,970 SNPs, information on 142 cultivars (lines) was used for the subsequent phylogenetic analysis.
In our previous report [20], we selected 4115 high-quality SNPs of “poly high resolution”, “no minor homozygote”, “mono high resolution”, and “hemizygous” categories from 179,970 SNP loci that had been used to construct the “LJ43×BHZ” linkage map. Therefore, the 4115 SNPs loci were used for the next step in developing genomic tools. The online web of bioinformatics (http://www.bioinformatics.com.cn/ accessed on 14 November 2022) was used for the visualization of SNP loci.

4.3. Core SNP Marker Selection

For convenience of scientific research and industrial applications in the future, from the 4115 SNP loci, 45 SNP markers were selected based on the following criteria: (1) SNP information was not missing in all cultivars; (2) minor allele frequency (MAF) >0.05; (3) the 35 bp up- and downstream areas of the SNP sites were completely conserved compared with the reference genome of tea plants (“Shuchazao”) [22], and the Blast tool in the Tea Plant Genomic Variations Database (TeaGVD) [39] (http://www.teaplant.top/teagvd accessed on 14 November 2022) was aligned to the reference genome “Shuchazao” [22]; (4) the polymorphism information content (PIC) was available. Finally, among the SNPs filtered based on the above criteria, three markers were selected on each chromosome based on the highest PIC index and 45 SNP markers evenly distributed on 15 chromosomes. Detailed genotype information of 117 tea cultivars on 45 SNP markers is shown in Table S2.

4.4. Data Analysis

The SNP genotype information was converted to a software readable format (such as A/A = 1/1, T/T = 2/2, and A/T = 1/2). The minor allele frequency (MAF), gene diversity, heterozygosity, PIC, and Nei’s genetic distances for markers were calculated in PowerMarker v3.25 [40,41]. The neighbor-joining method was used to construct the phylogeny tree based on genetic distance, and the phylogenetic tree was constructed using ITOL (https://itol.embl.de/ accessed on 14 November 2022). GenAlex 6.5 was used to obtain the number of effective alleles (Ne), observed heterozygosity (Ho), expected heterozygosity (He), and Shannon’s information index (I) [42].
In the same way, the probability of identity (PI) for each marker and their combinations were calculated using GenAlex 6.5 [10,36]. PI represents the average probability of two random individuals having the same genotype, and the computation formula is as follows:
PI = 2 (∑pi2) 2 − ∑pi4
where pi is the frequency of the ith allele at a locus. All sites were considered to be independently isolated. For multiple loci combinations, the PI was calculated as the product of the individual locus PI values, assuming that all loci segregated independently. In this study, we took into account that all loci should belong to independent segregations at different chromosomal locations; therefore, PIsibs was not calculated [10].
The probability of exclusion (PE) values for each marker and the 45 marker combinations were calculated using GenAlex 6.5 [26]. A multi-locus combination is the exclusion probability when the genotypes of both parents are known (Q1) or when the genotype of only one parent is known (Q2), and when there are two putative parents unknown (Q3) [43].
Cervus 3.0 [25] was used to calculate the parentage analysis with exclusion and likelihood methods. The parameters were set as follows: the parents sampled, the mistyped, and the proportions of loci typed were set as 0.5, 0.01, and 0.99, respectively; the minimum typed loci, candidate fathers, and numbers for offspring were set as 13, 117, and 1000, respectively. Based on the simulation analysis, a critical LOD value at a 95% confidence level was calculated.

5. Conclusions

In this study, our data confirmed that Yunnan Province is a primary center of CSA origin and domestication, and we deduced that east China (mainly Zhejiang and Fujian Provinces) is most likely the area of origin and domestication of CSS. Moreover, 45 core markers were screened using strict criteria to obtain information from 179,970 SNPs. We analyzed 117 well-known tea cultivars in China with 45 core SNP markers. The 45 core SNP markers can be used as tools for the identification of germplasm, parentage verification, and phylogenetic analysis of tea plants.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/plants12010162/s1, Table S1: The detailed information of main characteristics and type of population structure of the tea cultivars analyzed in this study; Table S2: the polymorphism analysis of 4115 SNP loci; Table S3: the position and the up- and downstream sequences of the 45 SNP markers; Table S4: the detail genotypes information of 45 SNP markers in 117 tea cultivars; Figure S1: Probability of identity for increasing combinations of the 45 loci.

Author Contributions

L.W. (Liubin Wang), K.W., D.N. and L.W. (Liyuan Wang) conceived and designed the experiments; H.X., S.A., R.Z. and L.W. (Liyun Wu) helped to data collection and performed the experiments; K.W. and L.W. (Liyuan Wang) contributed reagents, materials, and analysis tools; L.W. (Liubin Wang), K.W. and L.W. (Liyuan Wang) wrote the paper. All authors have read and agreed to the published version of the manuscript.

Funding

We gratefully acknowledge the Zhejiang Science and Technology Major Program on Agricultural New Variety Breeding-Tea Plant (2021C02067-7-1) and China Agriculture Research System of MOF and MARA(CARS-19).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare that they have no conflict of interest. All authors read this paper and confirm the context.

References

  1. Huang, J.; Wang, Y.; Xie, Z.; Zhou, Y.; Zhang, Y.; Wan, X. The anti-obesity effects of green tea in human intervention and basic molecular studies. Eur. J. Clin. Nutr. 2014, 68, 1075–1087. [Google Scholar] [CrossRef] [Green Version]
  2. Ng, K.W.; Cao, Z.J.; Chen, H.B.; Zhao, Z.Z.; Zhu, L.; Yi, T. Oolong tea: A critical review of processing methods, chemical composition, health effects, and risk. Crit. Rev. Food Sci. Nutr. 2018, 58, 2957–2980. [Google Scholar] [CrossRef] [PubMed]
  3. Shu, X.; Cai, H.; Xiang, Y.B.; Li, H.; Lipworth, L.; Miller, N.L.; Zheng, W.; Shu, X.O.; Hsi, R.S. Green tea intake and risk of incident kidney stones: Prospective cohort studies in middle-aged and elderly Chinese individuals. Int. J. Urol. 2019, 26, 241–246. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Drew, L. The growth of tea. Nature 2019, 566, S2–S4. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Shehasen, M.Z. Tea Plant (Camellia Sinensis) Breeding Mechanisms Role in Genetic Improvement and Production of Major Producing Countries. Int. J. Res. Stud. Sci. Eng. Technol. 2019, 6, 10–20. [Google Scholar]
  6. Tan, L.Q.; Wang, L.Y.; Xu, L.Y.; Wu, L.Y.; Peng, M.; Zhang, C.C.; Wei, K.; Bai, P.X.; Li, H.L.; Cheng, H.; et al. SSR-based genetic mapping and QTL analysis for timing of spring bud flush, young shoot color, and mature leaf size in tea plant (Camellia sinensis). Tree Genet. Genomes 2016, 12, 52. [Google Scholar] [CrossRef]
  7. Wei, C.; Yang, H.; Wang, S.; Zhao, J.; Liu, C.; Gao, L.; Xia, E.; Lu, Y.; Tai, Y.; She, G.; et al. Draft genome sequence of Camellia sinensis ar. sinensis provides insights into the evolution of the tea genome and tea quality. Proc. Natl. Acad. Sci. USA 2018, 115, E4151–E4158. [Google Scholar] [CrossRef] [Green Version]
  8. Xia, E.H.; Zhang, H.B.; Sheng, J.; Li, K.; Zhang, Q.J.; Kim, C.; Zhang, Y.; Liu, Y.; Zhu, T.; Li, W.; et al. The Tea Tree Genome Provides Insights into Tea Flavor and Independent Evolution of Caffeine Biosynthesis. Mol. Plant 2017, 10, 866–877. [Google Scholar] [CrossRef] [Green Version]
  9. Liu, S.; Liu, H.; Wu, A.; Hou, Y.; An, Y.; Wei, C. Construction of fingerprinting for tea plant (Camellia sinensis) accessions using new genomic SSR markers. Mol. Breed. 2017, 37, 93. [Google Scholar] [CrossRef]
  10. Tan, L.Q.; Peng, M.; Xu, L.Y.; Wang, L.Y.; Chen, S.X.; Zou, Y.; Qi, G.N.; Cheng, H. Fingerprinting 128 Chinese clonal tea cultivars using SSR markers provides new insights into their pedigree relationships. Tree Genet. Genomes 2015, 11, 90. [Google Scholar] [CrossRef]
  11. Waits, L.P.; Luikart, G.; Taberlet, P. Estimating the probability of identity among genotypes in natural populations: Cautions and guidelines. Mol. Ecol. 2001, 10, 249–256. [Google Scholar] [CrossRef]
  12. Yao, M.Z.; Ma, C.L.; Qiao, T.T.; Jin, J.Q.; Chen, L. Diversity distribution and population structure of tea germplasms in China revealed by EST-SSR markers. Tree Genet. Genomes 2012, 8, 205–220. [Google Scholar] [CrossRef]
  13. Rafalski, A. Applications of single nucleotide polymorphisms in crop genetics. Curr. Opin. Plant Biol. 2002, 5, 94–100. [Google Scholar] [CrossRef] [PubMed]
  14. Nasu, S.; Suzuki, J.; Ohta, R.; Hasegawa, K.; Yui, R.; Kitazawa, N.; Monna, L.; Minobe, Y. Search for and analysis of single nucleotide polymorphisms (SNPs) in rice (Oryza sativa, Oryza rufipogon) and establishment of SNP markers. DNA Res. 2002, 9, 163–171. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Korir, N.K.; Han, J.; Shangguan, L.; Wang, C.; Kayesh, E.; Zhang, Y.; Fang, J. Plant variety and cultivar identification: Advances and prospects. Crit. Rev. Biotechnol. 2013, 33, 111–125. [Google Scholar] [CrossRef]
  16. Fang, W.P.; Meinhardt, L.W.; Tan, H.W.; Zhou, L.; Mischke, S.; Zhang, D. Varietal identification of tea (Camellia sinensis) using nanofluidic array of single nucleotide polymorphism (SNP) markers. Hortic. Res. 2014, 1, 14035. [Google Scholar] [CrossRef] [Green Version]
  17. Sim, S.C.; Robbins, M.D.; Chilcott, C.; Zhu, T.; Francis, D.M. Oligonucleotide array discovery of polymorphisms in cultivated tomato (Solanum lycopersicum L.) reveals patterns of SNP variation associated with breeding. BMC Genom. 2009, 10, 466. [Google Scholar] [CrossRef] [Green Version]
  18. Huang, X.F.; Wei, Y.L.; Yuan, Y.; Yu, J.C.; Yuan, Z.Y.; Wang, W.S. Genetic diversity analysis and core collection construction of 221 cultivars (strains) of Litchi chinensis based on SNP molecular markers. J. Plant Resour. Environ. 2022, 31, 74–84. [Google Scholar]
  19. Fang, W.P.; Meinhardt, L.W.; Tan, H.; Zhou, L.; Mischke, S.; Wang, X.; Zhang, D. Identification of the varietal origin of processed loose-leaf tea based on analysis of a single leaf by SNP nanofluidic array. Crop J. 2016, 4, 304–312. [Google Scholar] [CrossRef] [Green Version]
  20. Wei, K.; Wang, X.; Hao, X.; Qian, Y.; Li, X.; Xu, L.; Ruan, L.; Wang, Y.; Zhang, Y.; Bai, P.; et al. Development of a genome-wide 200K SNP array and its application for high-density genetic mapping and origin analysis of Camellia sinensis. Plant Biotechnol. J. 2022, 20, 414–416. [Google Scholar] [CrossRef]
  21. Meegahakumbura, M.K.; Wambulwa, M.C.; Thapa, K.K.; Li, M.M.; Möller, M.; Xu, J.C.; Yang, J.B.; Liu, B.Y.; Ranjitkar, S.; Liu, J.; et al. Indications for Three Independent Domestication Events for the Tea Plant (Camellia sinensis (L.) O. Kuntze) and New Insights into the Origin of Tea Germplasm in China and India Revealed by Nuclear Microsatellites. PLoS ONE 2016, 11, E0155369. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Meegahakumbura, M.K.; Wambulwa, M.C.; Li, M.M.; Thapa, K.K.; Sun, Y.S.; Möller, M.; Xu, J.C.; Yang, J.B.; Liu, J.; Liu, B.Y.; et al. Domestication Origin and Breeding History of the Tea Plant (Camellia sinensis) in China and India Based on Nuclear Microsatellites and cpDNA Sequence Data. Front. Plant Sci. 2017, 8, 2270. [Google Scholar] [CrossRef] [PubMed]
  23. Xia, E.H.; Tong, W.; Hou, Y.; An, Y.L.; Chen, L.B.; Wu, Q.; Liu, Y.L.; Yu, J.; Li, F.D.; Li, R.P.; et al. The Reference Genome of Tea Plant and Resequencing of 81 Diverse Accessions Provide Insights into Its Genome Evolution and Adaptation. Mol. Plant 2020, 13, 1013–1026. [Google Scholar] [CrossRef] [PubMed]
  24. Yao, M.Z.; Chen, L. Tea Germplasm and Breeding in China. In Global Tea Breeding: Achievements, Challenges and Perspectives; Springer: Berlin/Heidelberg, Germany, 2012; pp. 13–68. [Google Scholar]
  25. Kalinowski, S.T.; Taper, M.L.; Marshall, T.C. Revising how the computer program CERVUS accommodates genotyping error in creases success in paternity assignment. Mol. Ecol. 2007, 16, 1099–1106. [Google Scholar] [CrossRef]
  26. Wang, J.; Santure, A.W. Parentage and Sibship Inference from Multilocus Genotype Data Under Polygamy. Genetics 2009, 181, 1579–1594. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Yang, Y.J.; Liang, Y.R. Clonal Tea Cultivars in China; Shanghai Scientific and Technical Publishers: Shanghai, China, 2014. [Google Scholar]
  28. Zhang, Y.; Wang, L.; Wei, K.; Ruan, L.; Wu, L.; He, M.; Tong, H.; Cheng, H. Differential regulatory mechanisms of secondary metabolites revealed at different leaf positions in two related tea cultivars. Sci. Hortic. 2020, 272, 109579. [Google Scholar] [CrossRef]
  29. Chen, L.; Apostolides, Z.; Chen, Z. Global Tea Breeding: Achievements, Challenges and Perspectives; Advanced Topics in Science & Technology in China: Beijing, China, 2013. [Google Scholar]
  30. Xia, E.H.; Tong, W.; Wu, Q.; Wei, S.; Zhao, J.; Zhang, Z.Z.; Wei, C.L.; Wan, X.C. Tea plant genomics: Achievements, challenges and perspectives. Hortic. Res. 2020, 7, 7. [Google Scholar] [CrossRef] [Green Version]
  31. Yang, H.; Wei, C.L.; Liu, H.W.; Wu, J.L.; Li, Z.G.; Zhang, L.; Jian, J.B.; Li, Y.Y.; Tai, Y.L.; Zhang, J.; et al. Genetic Divergence between Camellia sinensis and Its Wild Relatives Revealed via Genome-Wide SNPs from RAD Sequencing. PloS ONE 2016, 11, E0151424. [Google Scholar] [CrossRef] [Green Version]
  32. Cheng, Q.K. Understanding of 6000-year-old human planted tea tree roots discovered in Tianluoshan Site. China Tea 2016, 2, 30–33. (In Chinese) [Google Scholar]
  33. Xu, X.; Liu, X.; Ge, S.; Jensen, J.D.; Hu, F.; Li, X.; Dong, Y.; Gutenkunst, R.N.; Fang, L.; Huang, L.; et al. Resequencing 50 accessions of cultivated and wild rice yields markers for identifying agronomically important genes. Nat. Biotechnol. 2011, 30, 105–111. [Google Scholar] [CrossRef]
  34. Kuang, M.; Wei, S.j.; Wang, Y.Q.; Zhou, D.Y.; Ma, L.; Fang, D.; Yang, W.-h.; Ma, Z.Y. Development of a core set of SNP markers for the identification of upland cotton cultivars in China. J. Integr. Agric. 2016, 15, 954–962. [Google Scholar] [CrossRef]
  35. Kim, M.; Jung, J.K.; Shim, E.J.; Chung, S.M.; Park, Y.; Lee, G.P.; Sim, S.C. Genome-wide SNP discovery and core marker sets for DNA barcoding and variety identification in commercial tomato cultivars. Sci. Hortic. 2021, 276, 109734. [Google Scholar] [CrossRef]
  36. Fan, X.J.; Yu, W.T.; Cai, C.P.; Lin, Y.; Wang, Z.H.; Fang, W.P.; Zhang, J.M.; Ye, N.X. Construction of Molecular ID for Tea Cultivars by Using of Singlenucleotide Polymorphism (SNP) Markers. Sci. Agric. Sin. 2021, 54, 1751–1772. [Google Scholar]
  37. Yang, T.; Miller, M.; Forgacs, D.; Derr, J.; Stothard, P. Development of SNP-Based Genomic Tools for the Canadian Bison Industry: Parentage Verification and Subspecies Composition. Front. Genet. 2020, 11, 585999. [Google Scholar] [CrossRef] [PubMed]
  38. Wang, X.C.; Feng, H.; Chang, Y.X.; Ma, C.L.; Wang, L.Y.; Hao, X.Y.; Li, A.L.; Cheng, H.; Wang, L.; Cui, P.; et al. Population sequencing enhances understanding of tea plant evolution. Nat. Commun. 2020, 11, 4447. [Google Scholar] [CrossRef]
  39. Chen, J.D.; He, W.Z.; Chen, S.; Chen, Q.Y.; Ma, J.Q.; Jin, J.Q.; Ma, C.L.; Moon, D.G.; Ercisli, S.; Yao, M.Z.; et al. TeaGVD: A comprehensive database of genomic variations for uncovering the genetic architecture of metabolic traits in tea plants. Front. Plant Sci. 2022, 13, 1056891. [Google Scholar] [CrossRef]
  40. Liu, K.; Muse, S.V. PowerMarker: An integrated analysis environment for genetic marker analysis. Bioinformatics 2005, 21, 2128–2129. [Google Scholar] [CrossRef] [Green Version]
  41. Nei, M.; Tajima, F.; Tateno, Y. Accuracy of estimated phylogenetic trees from molecular data. II. Gene frequency data. J. Mol. Evol. 1983, 19, 153–170. [Google Scholar] [CrossRef]
  42. Peakall, R.; Smouse, P.E. GenAlEx 6.5: Genetic analysis in Excel. Population genetic software for teaching and research--an update. Bioinformatics 2012, 28, 2537–2539. [Google Scholar] [CrossRef] [Green Version]
  43. Tan, L.Q.; Liu, Q.L.; Zhou, B.; Yang, C.J.; Zou, X.; Yu, Y.Y.; Wang, Y.; Hu, J.-H.; Zou, Y.; Chen, S.-X.; et al. Paternity analysis using SSR markers reveals that the anthocyanin-rich tea cultivar ‘Ziyan’ is self-compatible. Sci. Hortic. 2019, 245, 258–262. [Google Scholar] [CrossRef]
Plants 12 00162 g001 550
Figure 1. Geographical origins of tea cultivars from different provinces in China; green, Assam type (CSA); orange, Chinese type (CSS); and yellow, transitional type. The different colored areas in the pie chart represent the proportion relative to the total number of tea cultivars. Scale bar: 200 Km.
Figure 1. Geographical origins of tea cultivars from different provinces in China; green, Assam type (CSA); orange, Chinese type (CSS); and yellow, transitional type. The different colored areas in the pie chart represent the proportion relative to the total number of tea cultivars. Scale bar: 200 Km.
Plants 12 00162 g001
Plants 12 00162 g002 550
Figure 2. Genomic location distribution and MAF and PIC values of 4115 SNP loci: (a) genomic location distribution of all SNP loci; the number represents the number of markers on the chromosome; (b,c) distribution of the respective MAF and PIC values of the SNP loci.
Figure 2. Genomic location distribution and MAF and PIC values of 4115 SNP loci: (a) genomic location distribution of all SNP loci; the number represents the number of markers on the chromosome; (b,c) distribution of the respective MAF and PIC values of the SNP loci.
Plants 12 00162 g002
Plants 12 00162 g003 550
Figure 3. SNP fingerprinting of 117 tea cultivars based on eight SNP markers.
Figure 3. SNP fingerprinting of 117 tea cultivars based on eight SNP markers.
Plants 12 00162 g003
Plants 12 00162 g004 550
Figure 4. Probability of exclusion relative to the increase in the number of combinations of the 45 SNP markers. Multi-locus combination is the exclusion probability when the genotypes of both parents are known (Q1) or when the genotype of only one parent is known (Q2) and when there are two putative parents unknown (Q3).
Figure 4. Probability of exclusion relative to the increase in the number of combinations of the 45 SNP markers. Multi-locus combination is the exclusion probability when the genotypes of both parents are known (Q1) or when the genotype of only one parent is known (Q2) and when there are two putative parents unknown (Q3).
Plants 12 00162 g004
Plants 12 00162 g005 550
Figure 5. Phylogenetic analysis of the 117 tea cultivars based on the 45 SNP markers.
Figure 5. Phylogenetic analysis of the 117 tea cultivars based on the 45 SNP markers.
Plants 12 00162 g005
Table
Table 1. The position of the 45 SNPs’ markers on the chromosome.
Table 1. The position of the 45 SNPs’ markers on the chromosome.
Name of MarkerChrPositionName of MarkerChrPosition
1_SNP_11205556811_SNP_9997146833
2_SNP_112034929322_SNP_9914194449
3_SNP_11425723563_SNP_9933065116
1_SNP_22795154021_SNP_1010120082494
2_SNP_22388222592_SNP_1010162024357
3_SNP_22686966543_SNP_101014142424
1_SNP_33215640281_SNP_1111103843763
2_SNP_331607121622_SNP_1111120074865
3_SNP_331340580463_SNP_1111115906813
1_SNP_44965652771_SNP_1212157283886
2_SNP_44148987102_SNP_121284717331
3_SNP_441708313553_SNP_1212119752353
1_SNP_5530736491_SNP_1313110373506
2_SNP_55123985692_SNP_131368674252
3_SNP_55564020233_SNP_1313103578624
1_SNP_661782684441_SNP_141433020570
2_SNP_66283142392_SNP_141452333124
3_SNP_66782653333_SNP_1414102052966
1_SNP_771750653961_SNP_151585702925
2_SNP_77289885762_SNP_1515108863635
3_SNP_77475708983_SNP_1515102549019
1_SNP_8842937619
2_SNP_88117013679
3_SNP_8843579005
Table
Table 2. Genetic diversity parameters of 45 SNP loci in 117 tea cultivars.
Table 2. Genetic diversity parameters of 45 SNP loci in 117 tea cultivars.
Name of MarkerMAFNaNeIHoPICPI
1_SNP_10.4932.0002.0000.6930.4870.3750.375
2_SNP_10.4682.0001.9910.6910.5560.3750.375
3_SNP_10.4582.0001.9910.6910.2990.3740.376
1_SNP_20.4892.0002.0000.6930.4700.3750.375
2_SNP_20.4862.0001.9930.6910.3930.3740.376
3_SNP_20.4582.0001.9790.6880.5040.3720.378
1_SNP_30.4822.0001.9980.6930.4700.3750.375
2_SNP_30.4752.0001.9760.6870.4270.3750.375
3_SNP_30.4722.0001.9950.6920.5380.3740.376
1_SNP_40.5002.0001.9990.6930.4870.3750.375
2_SNP_40.4932.0001.9990.6930.4020.3750.375
3_SNP_40.4612.0001.9540.6810.3500.3690.381
1_SNP_50.4862.0001.9990.6930.4530.3750.375
2_SNP_50.4582.0001.9320.6750.4870.3710.379
3_SNP_50.4582.0001.9950.6920.4870.3740.376
1_SNP_60.4752.0001.9960.6920.4620.3720.378
2_SNP_60.4722.0001.9980.6930.4020.3750.375
3_SNP_60.4682.0001.9540.6810.4700.3690.381
1_SNP_70.4822.0001.9850.6890.4360.3730.377
2_SNP_70.4542.0001.9930.6910.4440.3710.379
3_SNP_70.4472.0001.9540.6810.4360.3690.381
1_SNP_80.4892.0001.9930.6910.3930.3740.376
2_SNP_80.4862.0001.9930.6910.4440.3740.376
3_SNP_80.4792.0001.9930.6910.5300.3740.376
1_SNP_90.4962.0001.9990.6930.4360.3750.375
2_SNP_90.4862.0001.9990.6930.4790.3750.375
3_SNP_90.4862.0001.9950.6920.4700.3740.376
1_SNP_100.4932.0001.9990.6930.4270.3750.375
2_SNP_100.4862.0001.9980.6930.5210.3750.375
3_SNP_100.4822.0001.9960.6920.4790.3750.375
1_SNP_110.4962.0001.9950.6920.5900.3740.376
2_SNP_110.4822.0001.9850.6890.5900.3730.377
3_SNP_110.4542.0001.9950.6920.4530.3740.376
1_SNP_120.4652.0001.9990.6930.4190.3740.376
2_SNP_120.4472.0001.9980.6930.4020.3740.376
3_SNP_120.4662.0001.9930.6910.4440.3740.376
1_SNP_130.4612.0001.9720.6860.5210.3750.375
2_SNP_130.4582.0001.9720.6860.4530.3660.384
3_SNP_130.4472.0001.9630.6840.4530.3700.380
1_SNP_140.4862.0001.9980.6930.5560.3750.375
2_SNP_140.4822.0001.9990.6930.4530.3750.375
3_SNP_140.4722.0002.0000.6930.4620.3750.375
1_SNP_150.4822.0001.9950.6920.4360.3740.376
2_SNP_150.4822.0001.9980.6930.5380.3750.375
3_SNP_150.4722.0001.9990.6930.4700.3750.375
Mean0.4752.0001.9890.6900.4640.3740.376
MAF, minor allele frequency; Na, number of allele; Ne, number of effective alleles; I, Shannon’s information index; Ho, observed heterozygosity; PIC, polymorphic information content; PI, probability of identity.
Table
Table 3. Putative parentage revealed by the fingerprinting data of the 45 SNP markers.
Table 3. Putative parentage revealed by the fingerprinting data of the 45 SNP markers.
Offspring 1Parenta 1Pair Loci MismatchingPair LOD ScorebCandidate Father IDPair Loci MismatchingPair LOD Score 2Known/New 3
HuangmeiguiHuangdan09.95 *-- Known
ChunguiHuangdan08.75 *-- Known
Mingke 1Huangdan 013.5 *Tieguanyin013.53 *Known
BenshanTieguanyin014.3 *-- New
XingrenchaTieguanyin08.58 *-- New
Zhongming 7Zhongcha 10806.2 *Longjing 435−13.32Known
ShuchazaoLongjing 43011 *-- Known
Zhongcha 108Longjing 43010.4 *-- Known
Foxiang 3Changye Baihao010.2 *-- Known
Aifeng Yungui 019.4 *-- New
1 In each pair, which cultivar is the parent or offspring is inferred from known pedigrees or their breeding information, such as year of registration. 2 LOD score, the natural log of the overall likelihood ratio; a higher LOD score represents a higher chance of being the true parent according to the likelihood method. -- represents that the candidate father’s ID is unknown. * 95% level of confidence. 3 “Known” means that their relationship is consistent with breeding records, and “New” refers to newly discovered information based on parentage verification in this study.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Wang, L.; Xun, H.; Aktar, S.; Zhang, R.; Wu, L.; Ni, D.; Wei, K.; Wang, L. Development of SNP Markers for Original Analysis and Germplasm Identification in Camellia sinensis. Plants 2023, 12, 162. https://doi.org/10.3390/plants12010162

AMA Style

Wang L, Xun H, Aktar S, Zhang R, Wu L, Ni D, Wei K, Wang L. Development of SNP Markers for Original Analysis and Germplasm Identification in Camellia sinensis. Plants. 2023; 12(1):162. https://doi.org/10.3390/plants12010162

Chicago/Turabian Style

Wang, Liubin, Hanshuo Xun, Shirin Aktar, Rui Zhang, Liyun Wu, Dejiang Ni, Kang Wei, and Liyuan Wang. 2023. "Development of SNP Markers for Original Analysis and Germplasm Identification in Camellia sinensis" Plants 12, no. 1: 162. https://doi.org/10.3390/plants12010162

APA Style

Wang, L., Xun, H., Aktar, S., Zhang, R., Wu, L., Ni, D., Wei, K., & Wang, L. (2023). Development of SNP Markers for Original Analysis and Germplasm Identification in Camellia sinensis. Plants, 12(1), 162. https://doi.org/10.3390/plants12010162

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop