Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Genetic evidence that Chinese chestnut cultivars in Japan are derived from two divergent genetic structures that originated in China

  • Sogo Nishio ,

    Roles Data curation, Formal analysis, Funding acquisition, Investigation, Resources, Writing – original draft

    nishios@affrc.go.jp

    Affiliation Institute of Fruit Tree and Tea Science, NARO, Tsukuba, Ibaraki, Japan

  • Shuan Ruan,

    Roles Resources

    Affiliation Liaoning Institute of Economic Forestry, Dalian, Liaoning, China

  • Yutaka Sawamura,

    Roles Investigation

    Affiliation Institute of Fruit Tree and Tea Science, NARO, Tsukuba, Ibaraki, Japan

  • Shingo Terakami,

    Roles Data curation, Formal analysis

    Affiliation Institute of Fruit Tree and Tea Science, NARO, Tsukuba, Ibaraki, Japan

  • Norio Takada,

    Roles Resources

    Affiliation Institute of Fruit Tree and Tea Science, NARO, Tsukuba, Ibaraki, Japan

  • Yukie Takeuchi,

    Roles Resources

    Affiliation Institute of Fruit Tree and Tea Science, NARO, Tsukuba, Ibaraki, Japan

  • Toshihiro Saito,

    Roles Supervision, Validation, Writing – review & editing

    Affiliation Institute of Fruit Tree and Tea Science, NARO, Tsukuba, Ibaraki, Japan

  • Eiich Inoue

    Roles Conceptualization, Data curation, Resources, Supervision, Validation, Writing – review & editing

    Affiliation College of Agriculture, Ibaraki University, Ami-machi, Ibaraki, Japan

Abstract

The Chinese chestnut (Castanea mollissima Bl.) was introduced into Japan about 100 years ago. Since then, a number of Chinese chestnut cultivars and Japanese–Chinese hybrid cultivars have been selected by farmers and plant breeders, but little information has been available about their origins and genetic relationships. A classification based on simple sequence repeat markers was conducted using 230 cultivars including Japanese chestnut (Castanea crenata Sieb. et Zucc.) cultivars originated in Japan, Japanese–Chinese hybrid cultivars, and Chinese chestnut cultivars originated in both Japan and China. First, a search for synonyms (cultivars with identical genotypes) revealed 23 synonym groups among the Chinese chestnut cultivars, and all but one cultivar from each synonym group was omitted from further analyses. Second, genetic structure analysis showed a clear division between Japanese and Chinese chestnut, and most of the Japanese and Chinese cultivars had a simple genetic structure corresponding to the expected species. On the other hand, most Japanese–Chinese hybrid cultivars had admixed genetic structure. Through a combination of parentage and chloroplast haplotype analyses, 16 of the 18 hybrid cultivars in this study were inferred to have parent–offspring relationships with other cultivars originated in Japan. Finally, Bayesian clustering and chloroplast haplotype analysis showed that the 116 Chinese chestnut cultivars could be divided into two groups: one originated in the Hebei region of China and the other originated in the Jiangsu and Anhui regions of China. The Chinese chestnut cultivars selected in Japan showed various patterns of genetic structure including Hebei origin, Jiangsu or Anhui origin, and admixed. The chestnut cultivar genetic classifications obtained in this study will be useful for both Japanese and Chinese chestnut breeding programs.

Introduction

Japanese chestnut (Castanea crenata Sieb. et Zucc.) and Chinese chestnut (C. mollissima Bl.) are naturally distributed throughout Japan and China, respectively. Both species are reported to have been prehistorically domesticated [13] and are still economically important for production of edible nuts. These two species have large differences in both genetic and morphological properties. Japanese chestnut has large fruit size, adherent pellicle, and comparatively good yield, whereas Chinese chestnut has small nut size, easy-to-peel pellicle, and low yield [4]. In addition, the yellow brown shoot and pubescence at the nut tip are characters specific to Chinese chestnut and can be used to distinguish between the two species.

Interspecific hybridization has been used in chestnut breeding programs all over the world [4,5]. The breeding objectives of such programs have varied depending on the species. In the United States, the chestnut blight resistance gene from Chinese chestnut has been introduced into American chestnut (C. dentata [Marsh.] Borkh.) [6,7]. To improve European chestnut (C. sativa Mill.), quantitative trait loci for resistance to Phytophthora cinnamomic and chestnut gall wasp were introduced from Japanese chestnut [8,9]. In Japan, Chinese chestnut has been applied in Japanese chestnut breeding programs to improve nut quality and pellicle peelability [10,11]. Thus, it is important to clarify the genetic relationships among chestnut cultivars of different species.

According to Tanaka [12] and Isaki [13], Chinese chestnut was introduced into Japan about 100 years ago and attempts were made to cultivate it in Japan. However, this species was not suited to the Japanese climate [12] and had high susceptibility to the chestnut gall wasp [13]. Some local farmers in Japan planted nuts of Chinese chestnut and selected for adaptation to the local climate and conditions, leading to the development of Chinese chestnut cultivars such as ‘Houji’, ‘Miyagawa’, ‘Hinoharu’, and ‘Aioi’ in the early 20th century [13] Subsequently, the Japanese–Chinese hybrid cultivar ‘Riheiguri’ was selected by a local farmer and became one of the leading cultivars in Japan, accounting for 7% of total chestnut cultivation area in Japan in 2016. Even though this cultivar does not give high yields, its nut quality is highly valued by consumers and producers.

A number of Chinese chestnut and Japanese–Chinese hybrid cultivars have been collected and preserved at the NARO (National Agriculture and Food Organization) Genebank (www.gene.affrc.go.jp) and can be readily used for research and breeding purposes. However, information on the origin and genetic structure of these cultivars has been unclear or limited, since most of them were developed by local farmers and local agricultural experiment stations more than 50 years ago. Also, it is quite difficult to distinguish between pure Chinese chestnut and Japanese–Chinese hybrids by nut appearance. Although Japanese and Chinese chestnut readily produce interspecific hybrids, the extent of admixing from other species has not been clarified by molecular marker analyses.

The Chinese chestnut has the widest distribution among the chestnut species [14]. The northern range reaches 41°N, following the ranges of the Yanshan Mountains, and the southern range extends to 18°N on the Wuzhi Mountain of Hainan Island. Due to its large distribution, the Chinese chestnut is presumed to have higher genetic diversity than other chestnut species [15]. Several reports have suggested that central China, possibly the Shengnongia region near the Chang Jiang River, is the center of genetic diversity of Chinese chestnut and is one of the candidates for a refugium of this species [1517]. On the other hand, Zhang and Liu [18] suggested that the southwest area of China is the center of diversity. So far, no obvious genetic structure corresponding to geographical location has been found by clustering analyses, suggesting that human-mediated transportation might have affected the wild chestnut population structure [17]. While a large number of chloroplast haplotypes were identified within wild chestnut populations, only two haplotypes were identified in a collection of cultivars [19]. Ovesná et al. [20] also suggested that the genetic variability of Chinese chestnut cultivars was less than that of wild populations. Thus, several studies have examined the genetic diversity of wild chestnut and the genetic relationship between wild populations and cultivar collections. However, classifications based on clustering analyses using a large number of chestnut cultivars collected from diverse regions have not previously been conducted. In addition, the genetic relationship between Chinese chestnut cultivars within and outside of China has been unclear.

The classification and clustering of chestnut cultivars in Europe and Japan have been assessed using simple sequence repeats (SSRs) [2125]. Because SSRs are highly reliable markers, they have been used to conduct Bayesian structure analyses and to identify synonyms (differently named cultivars with identical genotypes), homonyms (different genotypes with the same cultivar name), and parent–offspring relationships. The Bayesian clustering analyses usually correspond to region-based classification or prevalent nut use. On the other hand, humans have carried scions between different locations, resulting in some cultivars showing unexpected cultivar origin and population structure [25]. In fact, many synonym groups contain local cultivars from different regions [21,24] and parentage analyses have clarified that traditional cultivars have contributed to the appearance of many local cultivars [23,25].

The main objective of the present study was to clarify the origin and genetic characteristics of Chinese chestnut cultivars in Japan. We used a three-step strategy for cultivar classification. The first step was to identify synonym groups and eliminate duplicate genotypes prior to further analyses. The second was to clarify the genetic relationship between Japanese–Chinese hybrids and Chinese chestnut cultivars in Japan by using highly reliable materials as references, i.e., Japanese and Chinese chestnut cultivars originated in Japan and China, respectively. The third was to clarify the genetic relationship between Chinese chestnut cultivars selected in Japan and those originated in China. Identification of synonyms, parentage, and genetic relationships among cultivars would be useful for chestnut breeding programs and further genetic classification studies.

Materials and methods

Plant materials and DNA extraction

The 12 groups (230 cultivars) used in this study are shown in Tables 1 and S1. The Japanese chestnut, Japanese–Chinese hybrid, and Chinese chestnut cultivars from Japan are preserved at the NARO Genebank (www.gene.affrc.go.jp). These include Japanese chestnut local cultivars that originated in the Kanto region of Japan (designated J_KA), the Tanba region of Japan (J_TA), other regions of Japan (J_OJ), and Korea (KOR); Japanese–Chinese chestnut hybrids developed in Japan (HYB); and Chinese chestnut cultivars selected in Japan (C_SJ) and introduced from other countries (C_IO). The Chinese chestnut cultivars originated in China were provided by the Liaoning Province Economic Forest Research Institute, Shandong Institute of Pomology, and Hebei Agriculture and Forestry Academy of Sciences Changli Guoshu Institute. Chinese chestnut cultivar groups in China were defined by region of origin, i.e., Hebei (C_HE), Shandong (C_SH), Anhui (C_AN), Jiangsu (C_JI), and other regions of China (C_OR). To avoid using duplicate genotypes in the analysis, only one cultivar was used from each synonym group identified by Nishio et al. [24]. Genomic DNA was extracted from young leaves with a DNeasy Plant Mini Kit (Qiagen, Germany) according to the manufacturer’s instructions.

thumbnail
Table 1. Names, accession numbers, and genotype information for the 230 cultivars used in this study.

https://doi.org/10.1371/journal.pone.0235354.t001

Genetic markers

The 230 cultivars were genotyped for 31 nuclear SSRs [2628] (S2 Table) and 5 chloroplast SSRs (cpSSRs) (Cmcs1–3, Cmcs5, and Cmcs7) [29]. PCR amplification was performed in a 10-μL solution containing 5 μL of 2× Green GoTaq G2 Hot Start Master Mix (0.4 mM each dNTP, Taq DNA polymerase, and 4 mM MgCl2, pH 8.5; Promega, Madison, WI, USA), 20 pmol of each forward primer labeled with a fluorescent dye (5-FAM or 5-HEX) and unlabeled reverse primer, and 2.5 ng of genomic DNA. Amplification was performed in 35 cycles of 94°C for 1 min, 55°C for 1 min, and 72°C for 2 min. PCR products were separated and detected with a 3130xl Genetic Analyzer (Life Technologies, Carlsbad, CA, USA). The size of each amplified band was determined by comparison with a set of internal-standard DNA fragments (400HD ROX, Life Technologies) in GeneMapper software v. 5.0 (Life Technologies).

Data analysis

Prior to analysis of population structure and parent–offspring relationships, synonym groups were identified using the 31 nuclear SSR markers to analyze the 230 cultivars. After eliminating duplicate genotypes, a set of 200 unique cultivars (see Results) was used for further analyses. The probability of identity (PI) for each locus and for the whole SSR set (Cumulative PI) was calculated using the software Gimlet v1.3.3 [30] to check the power of discrimination. Chloroplast haplotypes were determined by using the 5 cpSSRs (S1 Table). Cultivars that had the identical combination of alleles for all 5 cpSSRs were considered to carry the same haplotype. The positions of the 5 cpSSRs within the complete Castanea mollissima chloroplast genome are shown in S1 Table.

Bayesian statistical inference on the population structure was performed by using Structure 2.3.4 software [31] with the independent model for allele frequency, without any prior information about the origin of each cultivar. First, the 200 cultivars representing unique genotypes were used to clarify the genetic structure of Japanese and Chinese chestnut and their hybrids. Next, only the pure Chinese chestnut cultivars were used in a second structure analysis to examine genetic structures within that species. The analysis was run 10 times for each value of K (number of inferred ancestral populations) from 2 to 8 for 1,000,000 iterations after a burn-in period of 100,000 iterations. Evanno et al.’s [32] criterion of ΔK was used to estimate the appropriate K value. Principal coordinate analysis (PCoA) was performed in GenAlEx 6.5 from the pairwise genetic distances obtained with the covariance-standardized method [33].

For the Japanese–Chinese hybrid cultivar group, putative parent–offspring relationships were calculated with the parent calculation program MARCO [34], which identifies possible parents from among the genotypes in a set of cultivars. Genotypes were considered to have a parent–offspring relationship if they shared at least one allele per SSR locus, with the exception that a discrepancy at a single SSR locus was accepted to allow for possible genotyping errors, presence of null alleles, or mutation, as previously proposed [25,3538]. To determine whether a parent was the seed parent or pollen parent, chloroplast haplotype data were used.

For the Chinese chestnut groups that had more than 8 cultivars each (C_SJ, C_HE, C_SH, C_AN, C_JI, and C_OR), the observed heterozygosity (HO), expected heterozygosity (HE), and inbreeding coefficient (F) were calculated using GenAlEx v. 6.5 software [33], and allelic richness (AR, n = 9) was calculated using the R package Hierfstat [39].

Evaluation of phenotypic traits

Nut harvest date and nut weight were recorded in 2001–2003 for cultivars preserved at the NARO Genebank, which include Japanese chestnut cultivars, Japanese–Chinese hybrid cultivars, and Chinese chestnut cultivars preserved in Japan (S3 Table). To evaluate nut harvest date and nut weight, each bur was harvested when it had changed from green to brown and had begun to split open or had dropped. The nuts were removed from the burs and the number of nuts harvested on a given day was recorded for each tree. Burs and nuts were harvested every three or four days from late August to October. The harvest date for each nut was expressed as the number of days after July 31 (i.e., August 1 = day 1), and the average value of nut harvest date for each tree was used as its score for this trait. Nut weight (g) per nut was measured on a digital scale on each harvest date. The average nut weight was calculated as the total nut weight divided by the total number of intact nuts.

For Chinese chestnut cultivars originated in China, nut harvest date and nut weight could not be obtained because only DNA samples were available.

Results

Identification of synonym groups

We could differentiate all but 53 of the 230 chestnut cultivars with the 31 nuclear SSR markers. These 53 cultivars were divided into 23 synonym groups, each consisting of 2 to 5 cultivars having the same genotypes at all SSR loci (Table 2). Consequently, we identified 200 unique genotypes from the 230 cultivars. The probability of identity (PI) for each locus ranged from 0.016 for PRG79 to 0.269 for PRD83 (mean = 0.099), whereas the total PI was 6.56× 10−35. Synonym groups were found only within Chinese chestnut cultivars originated in China: none of the synonym groups contained cultivars originated in Japan. Some of the cultivars showing identical genotypes originated in different regions of China; however, most of them originated in regions adjacent to one another. Examples of synonym groups include ‘Yanshanduanzhi’ from Hebei and ‘Laiyangduanzhi’, ‘Laizhouduanzhi’, and ‘Shuheduanzhi’ from Shandong (Syn-10); ‘Yanshanzaofeng’ from Hebei and ‘Xinzhuang 2’ from Beijing (Syn-6); ‘Mifengqiu’ from Anhui and ‘Paoche 2’ from Jiangsu (Syn-2); and ‘Hefeichushuhong’ from Jiangsu and ‘Shuangjili’ from Anhui (Syn-1).

thumbnail
Table 2. Chinese chestnut cultivar groups with identical SSR genotypes for all 31 SSR markers.

https://doi.org/10.1371/journal.pone.0235354.t002

Chloroplast haplotype frequency

In total, 6 chloroplast haplotypes were identified among the 200 unique cultivars by using 5 cpSSRs (Table 3). The three Japanese chestnut cultivar groups (J_KA, J_TA, and J_OJ) carried only HAP1, whereas the Chinese chestnut cultivar groups mainly carried HAP3 and HAP4. HAP2 was only found in the Korean cultivar ‘Pochun B-1’ (ID#66). HAP5 and HAP6 were only found in Chinese chestnut cultivars ‘Kousei 2’ (ID#103) and ‘Duanzhatou’ (ID#175), respectively. Among the Chinese chestnut cultivars, HAP3 was the only haplotype found in C_JI and C_AN (East central regions), whereas HAP4 dominated in C_HE (Northern region). The cultivars originated in Shandong (C_SH), which is located between Hebei and Jiangsu, carried both HAP3 and HAP4. Chinese chestnut cultivars from other regions of China (C_OR) and those from Japan (C_SJ) and other world areas (C_IO) also carried both HAP3 and HAP4.

thumbnail
Table 3. Haplotype frequencies in the set of 200 cultivars.

https://doi.org/10.1371/journal.pone.0235354.t003

Genetic relationship between Japanese and Chinese chestnut cultivars and their hybrids

To clarify the genetic relationships between Japanese and Chinese chestnut cultivars and their hybrids, we performed Bayesian clustering analyses (Fig 1). The values of Δ(K) were much higher at K = 2 than at K = 3 to K = 8, so we constructed bar plot diagrams at K = 2. The “red” and “light green” shading indicate the Japanese and Chinese chestnut clusters, respectively. Most of the Japanese and Chinese chestnut cultivars had membership in a single cluster corresponding to their respective species. On the other hand, most Japanese–Chinese hybrid cultivars (HYB) had admixed structure, i.e., had approximately equal membership in both the “red” and “light green” clusters. Among the HYB cultivars, only ‘Ikaba’ (ID#67) had membership predominantly in the “light green” (Chinese) cluster. In both the Japanese and Chinese chestnut cultivar groups, a few cultivars had admixed structures (ratio of red:light green = 0.78:0.22 for ‘Obiwase’, 0.24:0.76 for ‘Hamjung 3’, and 0.16:0.84; for ‘Miyagawa 18’). Some Chinese chestnut cultivars had very low membership in the “red” (Japanese) cluster, but the amount was too small to declare that they were hybrids (for example, red:light green = 0.04:0.96 for ‘Tokuganri A’ and 0.03:0.97 for ‘2399’).

thumbnail
Fig 1. Detailed bar plot diagram for K = 2 in the independent model using 66 Japanese chestnut cultivars, 18 Japanese–Chinese hybrids, and 116 Chinese chestnut cultivars.

The first number under each bar represents the individual accession ID number (1–200); the second number (in parentheses) represents the group number (1–12). ID numbers and groups are defined in Table 1.

https://doi.org/10.1371/journal.pone.0235354.g001

Putative parentage of Japanese–Chinese hybrid cultivars

Because some Japanese–Chinese hybrid cultivars had records indicating that they had been selected from crosses between Japanese and Chinese chestnut cultivars in Japan, we used MARCO software to perform parent–offspring relationship analysis using cultivars originated in Japan (Table 4). Out of the 18 Japanese–Chinese cultivars in this study, we could infer both parents for 6 cultivars (‘Ikaba’, ‘Hayashi 1’, ‘Hayashi 3’, ‘Hayashi amaguri’, ‘Hyogo 493’, and ‘Yamewase’; ID#67–71, 84) and one parent for 10 cultivars (‘Kurakata amaguri’, ‘Nishiharima’, ‘Omatsuguri’, ‘Senshu amaguri’, ‘Shimaki 1’, ‘Shimaki 2’, ‘Shimaki 3’, ‘Shimaki 4’, ‘Shimaki 5’, and ‘Shimaki 6’; ID#72–74, 76–82). Out of the six cultivars for which we were able to infer both parents, five cultivars were F1 hybrids between Japanese and Chinese chestnut and one cultivar, ‘Ikaba’ (ID#67), was presumed to be an offspring between Japanese–Chinese hybrid cultivar ‘Riheiguri’ (ID#75) and Chinese chestnut cultivar ‘Gifu 1’ (ID#87).

thumbnail
Table 4. Putative parent–offspring relationships for Japanese–Chinese hybrid cultivars.

https://doi.org/10.1371/journal.pone.0235354.t004

Phenotypic trait evaluation of Japanese chestnut cultivars

Nut harvest date ranged from August 19 [‘Hassaku’ (ID#4)] to October 15 [‘Shimokatsugi’ (ID#32)] (average, September 18) for Japanese chestnut cultivars and from August 31 [‘Hayashi 3’ (ID#69)] to October 1 [‘Shimaki 1’ (ID#77)] (average, September 17) for Japanese–Chinese chestnut cultivars (S3 Table). The average nut harvest date of the Chinese chestnut cultivars preserved in Japan (September 29) was later than those of the Japanese and Japanese–Chinese chestnut hybrid cultivars and ranged from September 10 [‘Kahoku 10 (ID#98) and’‘Miyagawa 84’ (ID#106)] to October 17 [‘Konan 22’ (ID#100) and ‘Kousei 2’ (ID#103)]).

The average nut weight was large (24.7 g) in Japanese chestnut cultivars, intermediate (20.0 g) in Japanese–Chinese chestnut hybrids, and small (10.9 g) in Chinese chestnut cultivars preserved in Japan. The range was 11.9–40.9 g for Japanese chestnut cultivars, 11.9–29.5 g for Japanese–Chinese chestnut hybrids, and 5.2–17.6 g for Chinese chestnut cultivars in Japan. The average nut weight of Chinese chestnut cultivars that carried HAP3 (12.6 g) was larger than that of those carrying HAP4 (7.5g).

Genetic relationships among Chinese chestnut cultivars in Japan and China

To determine the genetic relationships among Chinese chestnut cultivars from different areas, Bayesian clustering analysis was performed using only the 116 Chinese chestnut cultivars (Fig 2). The values of Δ(K) were highest at K = 2: values of Δ(K) at K = 3 to K = 8 were less than one-thousandth that at K = 2. The two Bayesian clusters strongly corresponded to the results of chloroplast haplotype analysis (Fig 3). Membership in the “red” cluster was dominant in cultivars from northern China (C_HE), all but one of which were HAP4, whereas membership in the “light green” cluster was dominant in cultivars from east central China (C_JI and C_AN), all of which were HAP3. The C_SJ, C_SH, and C_OR cultivars had membership in both the “red” and “light green” clusters. Although most of the C_SH cultivars were admixed, many of the C_SJ cultivars had membership in a single cluster (ID#92–98, 105–109, 111). In C_SJ, cultivars that carried HAP3 (S1 Table) belonged to the “light green” or admixed clusters (Fig 2; ID#85–95, 97, 99–102, 108–110), while most of the cultivars that carried HAP4 belonged to the “red” cluster (ID#98, 105–107, 111). To validate the results of the Bayesian clustering analysis, PCoA was conducted after classifying each cultivar as predominantly part of the “red” cluster (ratio of “red” > 0.8), the “light green” cluster (ratio of “light green” > 0.8), or admixed (0.8 > ratio of “red” > 0.2). The results of PCoA were similar to those of the Structure analysis, and the first informative PCo component corresponded to the separation between the “red” and “green” clusters (S1 Fig).

thumbnail
Fig 2. Detailed bar plot diagram for K = 2 in the independent model using 117 Chinese chestnut cultivars.

The first number under each bar represents the individual accession ID number (1–200); the second number (in parentheses) represents the group number (6–12). ID numbers and groups are defined in Table 1. Red and green are used to indicate each of the K = 2 populations in this analysis and do not correspond to the same colors in Fig 1.

https://doi.org/10.1371/journal.pone.0235354.g002

thumbnail
Fig 3. Geographic locations of Chinese chestnut cultivar groups and genetic structures.

Donut and pie charts indicate composition of chloroplast haplotypes and clusters identified in STRUCTURE, respectively. The green, orange, and blue in the donut charts indicate HAP3, HAP4, and HAP5, respectively. The colors for the pie charts are based on the results shown in Fig 2.

https://doi.org/10.1371/journal.pone.0235354.g003

Genetic parameters were calculated to clarify the genetic diversity of Chinese chestnut cultivar groups that had more than 8 cultivars (Table 5). HO was the lowest for C_HE (0.457), whereas HO for C_AN (0.642) was slightly higher than for the other groups. HE and AR for C_HE were both lower than for the other groups. On the other hand, AR in C_SJ, C_SH, and C_OR, all of which had both HAP3 and HAP4 cultivars (Table 3), was higher than in the other groups. The inbreeding coefficient (F) was highest for C_SJ (0.085), whereas those for C_JI and C_AN were negative (−0.118 and −0.138, respectively).

thumbnail
Table 5. Genetic characteristics of Chinese chestnut cultivar groups analyzed using 31 SSRs.

https://doi.org/10.1371/journal.pone.0235354.t005

Discussion

Our set of 31 nuclear SSR loci proved to have an high discriminative power (total probability of identity: 6.56× 10−35) for the 200 unique cultivars. This value of the total probability of identity is quite low compared to those in other studies related to identification of synonyms (3.73× 10−12–2.99× 10−8) [4042]. It is highly unlikely to detect false synonyms with the 31 nuclear SSR markers. We identified 23 synonym groups among the Chinese chestnut cultivars originated in China. A previous study using Japanese and Chinese chestnut cultivars in Japan had already identified two synonym groups for Chinese chestnut cultivars selected in Japan and one for Japanese–Chinese hybrid cultivars [24]. Synonyms have been commonly identified in chestnut cultivar collections [2123]. One reason that chestnut cultivars may have many synonyms is that nut appearance is quite similar among cultivars. For most major fruit crops, fruit color would be a good characteristic to distinguish cultivars, but color differences are not helpful for distinguishing chestnut cultivars in most cases. Interestingly, synonyms were found only among cultivars originated in China, not between cultivars from Japan and China. Most Chinese chestnut genotypes introduced into Japan were assumed to be from seeds from China or Korea [13]. Because seeds can be produced via outcrossing, their progeny would be genetically diverse. Clonal propagation would have been difficult compared with seed propagation because of graft incompatibility between Japanese and Chinese chestnut. In addition, some of the Chinese chestnut cultivars might not have been suited to Japanese climates, reducing the chance of introducing the Chinese chestnut cultivars by clonal propagation. The cultivars in C_SJ would have been selected either directly from seed introduced from China and Korea or from successive generations of hybridization among the introduced genotypes.

Genetic relationships among Chinese chestnut, Japanese chestnut, and Japanese–Chinese hybrid cultivars were clarified by Bayesian clustering analysis. Most of the Chinese and Japanese chestnut cultivars had simple genetic structure (Fig 1), indicating that those cultivars were pure Chinese chestnut or Japanese chestnut, respectively. The cultivar collections from Japan and China have been present in those nations for a long time, limiting the chances for interspecific hybridization. According to Isaki [13], introduction of nuts of Chinese chestnut into Japan began in the 20th century. On the other hand, Korean native chestnut had been considered as an intermediate between Japanese and Chinese chestnut on the basis of its morphological characteristics [43]. Consistent with that previous report, our study showed that the Korean cultivar ‘Hamjung 3’ (ID#62) was an admixture between Chinese and Japanese chestnut (Fig 1).

In addition to ‘Hamjung 3’, definitive introgressions of Chinese chestnut into Japanese chestnut and vice versa were identified in Japanese chestnut ‘Obiwase’ (ID#47) and Chinese chestnut ‘Miyagawa 18’ (ID#105). In a previous study, however, ‘Obiwase’ was presumed to have a genetic structure derived from a wild chestnut population distributed on Kyusyu Island; its structure was different from that of other native cultivars but it was still considered to be a pure Japanese chestnut cultivar [25]. Since this study did not include wild populations, this cultivar may have been miscategorized. Likewise, ‘Miyagawa 18’ carried both Chinese and Japanese chestnut genetic structure. The percentage of Japanese chestnut genetic structure was about 16%; thus, this cultivar would be a first backcross (BC1) or a second backcross (BC2), not an F1 hybrid. For Japanese–Chinese chestnut hybrid cultivars, both parentage and chloroplast haplotype analyses were conducted to determine the putative seed and parent cultivars. Out of 18 Japanese–Chinese cultivars, we were able to presume both parents for six cultivars and one parent for 10 cultivars (Table 4). Because Chinese chestnut cultivars were relatively uncommon in Japan and the hybridizations had been done relatively recently (within 100 years), many parent–offspring relationships were identified. The results of the analyses were compared with those of Isaki [13], who reported the parentages of some of the same cultivars. ‘Shimaki 1’, ‘Shimaki 2’, ‘Shimaki 3’, ‘Shimaki 4’, and ‘Shimaki 5’ (ID#77–81) were reported by Isaki [13] to be selected from seedlings derived from ‘Ganne’ and Chinese chestnut accessions. Here, the parentage of ‘Shimaki 1’, ‘Shimaki 2’, and ‘Shimaki 5’ was reconfirmed. Also, the parentage of ‘Ikaba’, listed in a plant variety protection database in Japan (http://www.hinshu2.maff.go.jp/en/en_top.html), matched our results. On the other hand, ‘Hayashi 1’, ‘Hayashi 3’, and ‘Hayashi amaguri’ were reported by Isaki [13] to be offspring of ‘Kasashi 1’, an F1 hybrid between ‘Kasaharawase’ (ID#43) and a Chinese chestnut accession; however, in the present study they were presumed to be offspring of ‘Kasaharawase’ or ‘Kanotsume’ (ID#24).

The Japanese–Chinese chestnut hybrid cultivars showed nut size intermediate (20.0 g) between those of Japanese (24.7 g) and Chinese chestnut cultivars (10.9 g). The average nut harvest date of the Japanese–Chinese chestnut hybrid cultivars was the similar to that of Japanese chestnut and earlier than that of Chinese chestnut (S3 Table). Thus, for Chinese chestnut breeding, Japanese chestnut cultivars have the potential to increase nut weight and shorten the time to harvest. On the other hand, introducing Chinese chestnut cultivars into Japanese chestnut breeding programs would be not an effective way to shorten the time to harvest or increase nut weight. Today, some Japanese–Chinese cultivars are highly valued by Japanese farmers and consumers because of the high nut quality and moderate nut pellicle peelability. QTL analyses of interspecific backcross populations are necessary to identify interesting species-specific genes and accelerate breeding programs.

The origin of hybrid cultivar ‘Riheiguri’ (ID#75), one of the major cultivars in Japan, is an interesting example of interspecies hybridization giving rise to a desirable new cultivar. ‘Riheiguri’ was assumed to be an F1 hybrid between Japanese and Chinese chestnut since it had been developed in a Tsuchida orchard that contained both species [13]. Although we were unable to identify its parents by parentage analysis, Bayesian structure analysis revealed it had approximately equal membership in both the “red” and “light green” clusters, supporting the assumption that it was an F1 hybrid. In addition, ‘Riheiguri’ had HAP3, which originated in Jiangsu and Anhui. This cultivar might have been selected from seeds of other cultivars originated in these regions. It has relatively large nut size like Japanese chestnut (23.9 g; S3 Table), mealy texture like Chinese chestnut [13], and moderate pellicle peelability [44], which was likely inherited from Chinese chestnut. Because of its good nut quality, ‘Riheiguri’ and its relatives have been used in breeding programs, resulting in the release of new hybrid cultivars such as ‘Shuhou’ and ‘Mikuri’. On the other hand, trials to release Japanese chestnut cultivars with the easy-peeling pellicle trait from Chinese chestnut have not yet succeeded [4].

We identified four chloroplast haplotypes (HAP3–HAP6) among the Chinese chestnut cultivars. HAP4 was mainly found in cultivars from Hebei, Shandong, and Japan, whereas HAP3 dominated in most groups except for C_HE. HAP5 and HAP6 were each identified in only one cultivar. At least 38 chloroplast haplotypes were identified from wild populations by Chen and Huang [19] and Liu et al. [17] using cpSSRs. However, only two haplotypes [19] and four haplotypes (this study) were identified from Chinese chestnut cultivar collections, suggesting that cultivars have limited genetic diversity compared to wild populations. Because the numbers and types of markers were quite similar in those studies and ours, it is reasonable to compare the results. On the other hand, it is possible that more chloroplast haplotypes would be found if we detected a larger number of polymorphisms by sequencing the whole chloroplast genome of several cultivars. The finding that cultivars showed less genetic diversity than wild populations corresponds to the suggestions of Mattioni et al. [45] and Ovesná et al. [20], i.e., that because traits and genes useful for chestnut cultivation were artificially selected, the domestication process would typically reduce genetic diversity.

Both chloroplast haplotype and Bayesian clustering analyses showed that the Chinese chestnut cultivars used in the present study could be divided into two groups (Fig 3): one originated in the Hebei region and the other originated in Jiangsu and Anhui. Most of the cultivars from Shandong had admixed genetic structure (Fig 2), whereas cultivars selected in other regions of China had various patterns of genetic structure. The cultivar groups that had admixed structure showed higher HE and AR, which is not unexpected because hybridization between cultivars from different clusters would increase genetic diversity. Since we had no information about cultivars selected in other regions of China, it was quite difficult for us to clarify the breeding history of these cultivars. Some cultivars from Shandong might have been selected from crosses between cultivars derived from Hebei and from Jiangsu or Anhui. Alternatively, cultivars might have been selected from wild chestnuts with an admixed genetic structure growing in Shandong. Although HE and AR were low in C_HE, C_JI, and C_AN, the values of F were positive in C_HE but negative in C_JI and C_AN, indicating that artificial selection pressure was higher in the Hebei region. According to Kikuchi [46], nuts of cultivars from Hebei are small and sweet, whereas cultivars from central China have comparatively large nut size and low flavor, and cultivars from Shandong have intermediate characteristics. Consistent with this previous report, the average nut weight of Chinese chestnut cultivars that carried HAP3, the only haplotype found in the cultivars from Jiangsu and Anhui, was larger than those carried HAP4, the most common haplotype found in the cultivars from Hebei (Tables 3 and S3).Thus, breeders could use cultivars showing differences in genetic structure according to the objectives of their breeding programs.

The Chinese chestnut cultivars selected in Japan showed various patterns of genetic structure. More than half of those cultivars showed a “light green”-dominated structure presumed to have originated in Jiangsu or Anhui, while some cultivars showed admixed structure or “red”-dominated structure presumed to have originated in Hebei. The cultivars selected by Houji in Kochi prefecture (‘Houji 354’, ‘Houji 445’, ‘Houji 350’, ‘Houji 446’, and ‘Houji 480’; ID#91–95) had high membership in the “light green” cluster, and the cultivars selected by Miyagawa in Yamanashi prefecture (‘Miyagawa 100’, ‘Miyagawa 18’, ‘Miyagawa 84’, and ‘Miyagawa 85’; ID#104–107) had high membership in the “red” cluster. According to PCoA analysis, the cultivars selected by Houji were closest to those originated in Jiangsu or Anhui, while the cultivars selected by Miyagawa were closest to those originated in Hebei (S1 Fig), suggesting that Houji cultivars were selected from seeds originated in Jiangsu or Anhui and that Miyagawa cultivars were selected from seeds originated in Hebei. Geographical data support this hypothesis because Hebei is located in the northern part of China, and Hakushu in Yamanashi, where Miyagawa cultivars were selected, is located in a cold, high-altitude part of Japan. According to one account, Miyagawa cultivars were introduced from Songchong in North Korea, which is relatively close to Hebei [13]. The nut sizes of Miyagawa cultivars (which had genetic structure originated in Hebei) were 5.2–9.2 g, whereas those of cultivars that had genetic structure originated in Jiangsu or Anhui averaged around 12 g (S3 Table), corresponding to previous reports that cultivars in northern China have smaller nut size than those in central China [13]. In this study, we analyzed only Chinese chestnuts from Japan and from the seaside of China, not from western or southern China. Fortunately, chloroplast haplotype analysis clarified that the Chinese chestnut materials from Japan and China were similar (Fig 3), suggesting that the Chinese chestnut cultivars collected in this study were sufficient basic materials to classify Chinese chestnut cultivars in Japan. However, we are sure that there are many other genetic structures in Chinese chestnuts in other regions. Further analyses, including cultivars and wild populations from all over China, would clarify the detailed genetic relationships and domestication process of Chinese chestnut.

Conclusions

We used SSRs to genotype chestnut cultivars preserved in both Japan and China and to determine the genetic structure of Chinese and Japanese chestnut cultivars. The synonym groups and putative parentages of some Chinese chestnut cultivars were identified here for the first time. Most of the Chinese and Japanese chestnut cultivars had a simple genetic structure corresponding to their respective species, whereas Japanese–Chinese hybrid cultivars had admixed structures. The Chinese chestnut cultivars could be divided into two groups: one that originated in Hebei and one that originated in Jiangsu and Anhui. The Chinese chestnut cultivars selected in Japan also carried the genetic structures originated in these two divergent regions, suggesting that their ancestral genotypes originated in those two groups. The information obtained in this study will be useful for population genetic studies for those species and for Japanese and Chinese chestnut breeding programs.

Supporting information

S1 Table. Names, accession numbers, and genotype information for the 230 cultivars used in this study.

https://doi.org/10.1371/journal.pone.0235354.s001

(PDF)

S2 Table. List of the 31 SSR nuclear markers used in the present study.

The linkage group (LG) and position (pos) of each marker are based on an integrated map of a Kunimi × 709–34 population (Kx709CP; Nishio et al., 2018).

https://doi.org/10.1371/journal.pone.0235354.s002

(XLSX)

S1 Fig. Principal coordinate analysis (PCoA) plot generated from genetic distance calculations among the 116 Chinese chestnut cultivars in GenAlEx software.

Cultivars were classified as predominantly part of the “red” cluster (ratio of “red” > 0.8), the “light green” cluster (ratio of “light green” > 0.8), or admixed (0.8 > ratio of “red” > 0.2; shown here in black).

https://doi.org/10.1371/journal.pone.0235354.s004

(TIF)

Acknowledgments

The authors thank Dr. Hiroyuki Iketani for his valuable suggestions and discussions.

References

  1. 1. Kitagawa J, Yasuda Y. The influence of climatic change on chestnut and horse chestnut preservation around Jomon sites in Northeastern Japan with special reference to the Sannai-Maruyama and Kamegaoka sites. Quaternary International. 2004;123:89–103.
  2. 2. Liu L, Zhou J. Some considerations on chestnut development in the 21st century in China. Acta Hortic. 1999;494:85–8.
  3. 3. Minaki M. Enlargement of the fruit size of Castanea crenata Sieb. et Zucc. since the Jomon Period (since ca. 10000 yBP) in Japan. Jpn J Histor Bot. 1994;2:3–10.
  4. 4. Pereira-Lorenzo S, Ballester A, Corredoira E, Vieitez AM, Agnanostakis S, Costa R, et al. Chestnut. In: Badenes ML, Byrne DH, editors. Fruit Breeding. Springer; 2012. pp. 729–69.
  5. 5. Pereira-Lorenzo S, Costa R, Anagnostakis S, Serdar U, Yamamoto T, Saito , et al. Interspecific hybridization of chestnut. In: Mason AS, editor. Polyploidy and hybridization for crop improvement. Boca Raton; 2016. pp 377–407.
  6. 6. Hebard FV. The backcross breeding program of the American chestnut foundation. J Amer Chestnut Found. 2006;19:55–77.
  7. 7. Wheeler N, Sederoff R. Role of genomics in the potential restoration of the American chestnut. Tree Genet Genomes 2009;5(1):181–7.
  8. 8. Santos C, Nelson CD, Zhebentyayeva T, Machado H, Gomes-Laranjo J, Costa RL. First interspecific genetic linkage map for Castanea sativa x Castanea crenata revealed QTLs for resistance to Phytophthora cinnamomi. PLoS One. 2017;12(9):e0184381. pmid:28880954
  9. 9. Torello-Marinoni D, Nishio S, Portis E, Valentini N, Sartor C, Dini F, et al., editors. Development of a genetic linkage map for molecular breeding of chestnut. Acta Hortic. 2018;1220:23–8.
  10. 10. Sawamura Y. Chestnut. In: The Japanese Society for Horticultural Science, editors. Horticulture in Japan 2003. Nakanishi Printing; 2006. pp. 82–5.
  11. 11. Pereira-Lorenzo S, Costa R, Anagnostakis S, Serdar U, Yamamoto T, Saito , et al. Interspecific hybridization of chestnut. In: Mason AS, editor. Polyploidy and hybridization for crop improvement. Boca Raton; 2016. pp 377–407.
  12. 12. Tanaka Y. Cultivation of chestnut. Tokyo: Meibundou; 1933.
  13. 13. Isaki M. Theory and actual condition of chestnut cultivation. Tokyo: Hakuyusya; 1978.
  14. 14. Fei S, Liang L, Paillet FL, Steiner KC, Fang J, Shen Z, et al. Modelling chestnut biogeography for American chestnut restoration. Divers Distrib. 2012;18(8):754–68.
  15. 15. Huang H, Dane F, Norton J. Allozyme diversity in Chinese, Seguin and American chestnut (Castanea spp.). Theor Appl Genet. 1994;88(8):981–5. pmid:24186251
  16. 16. Lang P, Huang H. Genetic variation and population structure of three endemic Castanea species in China. Acta Hortic. 1999;494:269–76.
  17. 17. Liu W, Kang M, Tian H, Huang HW. A range wide geographic pattern of genetic diversity and population structure of Castanea mollissima populations inferred from nuclear and chloroplast microsatellites. Tree Genet Genomes. 2013;9(4):975–87.
  18. 18. Zhang H, Liu L. The genetic diversity of Castanea mollissima and the effect of artificial selection. Acta Bot Yunnan. 1998;20(1):81–8.
  19. 19. Cheng L, Huang W, editors. Diversity of chloroplast DNA SSRs in wild and cultivated C. mollissima: evidence of cultivated C. mollissima origin. Acta Hortic. 2010;866:243–9.
  20. 20. Ovesná J, Kučera L, Jiang L, Vagnerová D. Characterisation of Chinese elite cultivars and genetic resources of chestnut by AFLP. Biol Plant. 2005;49(1):125–7.
  21. 21. Gobbin D, Hohl L, Conza L, Jermini M, Gessler C, Conedera M. Microsatellite-based characterization of the Castanea sativa cultivar heritage of southern Switzerland. Genome. 2007;50(12):1089–103. pmid:18059537
  22. 22. Marinoni DT, Akkak A, Beltramo C, Guaraldo P, Boccacci P, Bounous G, et al. Genetic and morphological characterization of chestnut (Castanea sativa Mill.) germplasm in Piedmont (north-western Italy). Tree Genet Genomes. 2013;9(4):1017–30.
  23. 23. Pereira-Lorenzo S, Costa R, Ramos-Cabrer AM, Ciordia-Ara M, Ribeiro CAM, Borges O, et al. Chestnut cultivar diversification process in the Iberian Peninsula, Canary Islands, and Azores. Genome. 2011;54(4):301–15. pmid:21491973
  24. 24. Nishio S, Yamamoto T, Terakami S, Sawamura Y, Takada N, Saito T. Genetic diversity of Japanese chestnut cultivars assessed by SSR markers. Breeding Sci. 2011;61(2):109–20.
  25. 25. Nishio S, Iketani H, Fujii H, Yamamoto T, Terakami S, Takada N, et al. Use of population structure and parentage analyses to elucidate the spread of native cultivars of Japanese chestnut. Tree Genet Genomes. 2014;10(5):1171–80.
  26. 26. Buck EJ, Hadonou M, James CJ, Blakesley D, Russell K. Isolation and characterization of polymorphic microsatellites in European chestnut (Castanea sativa Mill.). Mol Ecol Notes. 2003;3(2):239–41.
  27. 27. Kubisiak TL, Nelson CD, Staton ME, Zhebentyayeva T, Smith C, Olukolu BA, et al. A transcriptome-based genetic map of Chinese chestnut (Castanea mollissima) and identification of regions of segmental homology with peach (Prunus persica). Tree Genet Genomes. 2013;9(2):557–71.
  28. 28. Nishio S, Yamamoto T, Terakami S, Sawamura Y, Takada N, Nishitani C, et al. Novel genomic and EST-derived SSR markers in Japanese chestnuts. Scientia Horticulturae. 2011;130(4):838–46.
  29. 29. Sebastiani F, Carnevale S, Vendramin G. A new set of mono‐and dinucleotide chloroplast microsatellites in Fagaceae. Mol Ecol Notes. 2004;4(2):259–61.
  30. 30. Valière N. GIMLET: a computer program for analysing genetic individual identification data. Mol Ecol Notes. 2002;2(3):377–9.
  31. 31. Pritchard JK, Stephens M, Donnelly P. Inference of population structure using multilocus genotype data. Genetics. 2000;155(2):945–59. pmid:10835412
  32. 32. Evanno G, Regnaut S, Goudet J. Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol Ecol. 2005;14(8):2611–20. pmid:15969739
  33. 33. Peakall R, Smouse PE. GENALEX 6: genetic analysis in Excel. Population genetic software for teaching and research. Mol Ecol Notes. 2006;6(1):288–95.
  34. 34. Fujii H, Yamashita H, Hosaka F, Terakami S, Yamamoto T. Development of a software to presume the parent-child relationship using the result of DNA marker typing. Hort Res (Japan). 2010;9 (SUPPL. 1):34.
  35. 35. Di Vecchi-Staraz M, Bandinelli R, Boselli M, This P, Boursiquot JM, Laucou V, et al. Genetic structuring and parentage analysis for evolutionary studies in grapevine: Kin group and origin of the cultivar Sangiovese revealed. J Am Soc Hortic Sci. 2007;132(4):514–24.
  36. 36. Boursiquot JM, Lacombe T, Laucou V, Julliard S, Perrin FX, Lanier N, et al. Parentage of Merlot and related winegrape cultivars of southwestern France: discovery of the missing link. Australian Journal of Grape and Wine Research. 2009;15(2):144–55. PubMed PMID: ISI:000266170700004.
  37. 37. Cipriani G, Marrazzo MT, Peterlunger E. Molecular characterization of the autochthonous grape cultivars of the region Friuli Venezia Giulia—North-Eastern Italy. Vitis. 2010;49(1):29–38.
  38. 38. Lacombe T, Boursiquot JM, Laucou V, Di Vecchi-Staraz M, Peros JP, This P. Large-scale parentage analysis in an extended set of grapevine cultivars (Vitis vinifera L.). Theor Appl Genet. 2013;126(2):401–14. pmid:23015217
  39. 39. Goudet J. Hierfstat, a package for R to compute and test hierarchical F‐statistics. Mol Ecol Notes. 2005;5(1):184–6.
  40. 40. Mellano MG, Beccaro GL, Donno D, Marinoni DT, Boccacci P, Canterino S, et al. Castanea spp. biodiversity conservation: collection and characterization of the genetic diversity of an endangered species. Genet Resour Crop Evol. 2012;59(8):1727–41.
  41. 41. Marinoni DT, Akkak A, Beltramo C, Guaraldo P, Boccacci P, Bounous G, et al. Genetic and morphological characterization of chestnut (Castanea sativa Mill.) germplasm in Piedmont (north-western Italy). Tree Genet Genomes 2013;9(4):1017–30.
  42. 42. Bacchetta L, Rovira M, Tronci C, Aramini M, Drogoudi P, Silva A, et al. A multidisciplinary approach to enhance the conservation and use of hazelnut Corylus avellana L. genetic resources. Genet Resour Crop Evol. 2015;62(5):649–63.
  43. 43. Kim MJ, Lee U, Kim SC, Hwang MS, Lee MH, editors. Comparison of nut characteristics between Korean native chestnut accessions and prevailing cultivars cultivated in Korea. Acta Hortic. 2005;693:299–304.
  44. 44. Sato A, Tanaka K, Takada N, Sawamura Y, Hirabayashi T. Comparison of phenolic content of easily removed pellicle of Japanese chestnut 'Porotan' with other Japanese and Chinese chestnut cultivars. J Jpn Soc Hortic Sci. 2010;79(3):258–62.
  45. 45. Mattioni C, Cherubini M, Micheli E, Villani F, Bucci G. Role of domestication in shaping Castanea sativa genetic variation in Europe. Tree Genet Genomes. 2008;4(3):563–74.
  46. 46. Kikuchi A. Horticulture of fruit trees, vol 1. Tokyo: Yokendo; 1948.