#PAGE_PARAMS# #ADS_HEAD_SCRIPTS# #MICRODATA#

Korean soybean core collection: Genotypic and phenotypic diversity population structure and genome-wide association study


Authors: Namhee Jeong aff001;  Ki-Seung Kim aff002;  Seongmun Jeong aff003;  Jae-Yoon Kim aff003;  Soo-Kwon Park aff001;  Ju Seok Lee aff005;  Soon-Chun Jeong aff005;  Sung-Taeg Kang aff006;  Bo-Keun Ha aff007;  Dool-Yi Kim aff001;  Namshin Kim aff003;  Jung-Kyung Moon aff008;  Man Soo Choi aff001
Authors place of work: National Institute of Crop Science, Rural Development Administration, Wanju-gun, Jeollabuk-do, Republic of Korea aff001;  FarmHannong, Ltd., Daejeon, Republic of Korea aff002;  Genome Editing Research Center, Korea Research Institute of Bioscience and Biotechnology, Daejeon, Republic of Korea aff003;  Department of Bioinformatics, KRIBB School of Bioscience, Korea University of Science and Technology, Daejeon, Republic of Korea aff004;  Bio-Evaluation Center, Korea Research Institute of Bioscience and Biotechnology, Cheongju, Chungcheongbuk-do, Republic of Korea aff005;  Department of Crop Science & Biotechnology, Dankook University, Cheonan, Chungcheongnam-do, Republic of Korea aff006;  Division of Plant Biotechnology, College of Agriculture and Life Sciences, Chonnam National University, Gwangju, Republic of Korea aff007;  National Institute of Agricultural Sciences, Rural Development Administration, Jeonju, Jeollabuk-do, Republic of Korea aff008
Published in the journal: PLoS ONE 14(10)
Category: Research Article
doi: https://doi.org/10.1371/journal.pone.0224074

Summary

A core collection is a subset that represents genetic diversity of the total collection. Soybean (Glycine max (L.) Merr.) is one of major food and feed crops. It is the world’s most cultivated annual herbaceous legume. Constructing a core collection for soybean could play a pivotal role in conserving and utilizing its genetic variability for research and breeding programs. To construct and evaluate a Korean soybean core collection, genotypic and phenotypic data as well as population structure, were analyzed. The Korean soybean core collection consisted of 430 accessions selected from 2,872 collections based on Affymetrix Axiom® 180k SoyaSNP array data. The core collection represented 99% of genotypic diversity of the total collection. Analysis of population structure clustered the core collection into five subpopulations. Accessions from South Korea and North Korea were distributed across five subpopulations. Analysis of molecular variance indicated that only 2.01% of genetic variation could be explained by geographic origins while 16.18% of genetic variation was accounted for by subpopulations. Genome-wide association study (GWAS) for days to flowering, flower color, pubescent color, and growth habit confirmed that the core collection had the same genetic diversity for tested traits as the total collection. The Korean soybean core collection was constructed based on genotypic information of the 180k SNP data. Size and phenotypic diversity of the core collection accounted for approximately 14.9% and 18.1% of the total collection, respectively. GWAS of core and total collections successfully confirmed loci associated with tested traits. Consequently, the present study showed that the Korean soybean core collection could provide fundamental and practical material and information for both soybean genetic research and breeding programs.

Keywords:

Phenotypes – Population genetics – Genome-wide association studies – Crops – Soybean – Plant growth and development – Genetic polymorphism – Crop genetics

Introduction

Genetic diversity is fundamental in crop breeding and research programs. Although genetic diversity could be a barrier to progress, it is essential for crop improvement. Monoculture of a few improved cultivars has decreased the genetic diversity across crop species [1]. One of important ways to achieve successful crop improvement is through continuous supply of genetic diversity, including new or improved variability in target traits [2]. Therefore, managing and utilizing large and diverse germplasm collection are important and valuable challenges for successful biological research and crop improvement programs.

Across the world, germplasm conservation centers have been established to preserve genetic diversity of target crop species. Although large numbers of collections are desirable to preserve genetic variability, their usefulness and accessibility might be inversely related to their size [3]. An increase in the number of accessions without detailed information could decrease the efficiency and value of these collections [4, 5]. To overcome the size issue, core collections have been constructed for important crop species. A core collection was first defined as “a limited set of accessions” representing genetic diversity with minimum repetitiveness in the total collection [1, 6, 7]. Presently, core collections for approximately 30 species have been constructed, including rice (Oryza sativa L.) [8], wheat (Triticum aestivum L.) [9], cotton (Gossypiumbar badense L.) [10], peanut (Arachis hypogea L.) [11], pepper (Capsicum spp.) [12], jujube (Ziziphus jujube Mill.) [13], melon (Cucumis melo L.) [14], medicago (Medicago spp.) [15], and soybean (Glycine max L. Merr.) [16].

Soybean is one of the most important crop species as both food and feed sources in the world. Soybean genetic resources are largely preserved in China, the USA, Japan, and Korea [6, 17, 18, 19, 20]. In Korea, approximately 20,000 accessions of soybean are maintained in the National Agrobiodiversity Center of Rural Development Administration (http://genebank.rda.go.kr/). China has the largest soybean collection. There are approximately 26,000 accessions in the Institute of Crop Germplasm Resources of the Chinese Academy of Agricultural Science [19, 20]. The United States Department of Agriculture (USDA) has nearly 17,000 soybean accessions [6]. In Japan, around 11,000 soybean accessions are maintained at the National Institute of Agrobiological Sciences (NIAS) Genebank [17].

Many studies have evaluated various strategies to construct core collections of many crop species [6, 17, 18, 21, 22]. These strategies include random, proportion, constant, logarithmic, and genetic diversity-based methods. Among them, genetic diversity-based method has been considered the simplest and the most efficient [21]. Through genetic diversity-based strategy using 48 SNPs, Lee et al. [12] have elected 240 accessions from 3,821 capsicum collections. Xu et al. [13] have selected a core collection consisting of 150 accessions from 947 jujube collections using 24 SSR markers. However, considering genome sizes of crop species, the above studies used relatively few genetic markers to select accessions for core collections.

Although a core collection of soybean was first constructed in 1987 [23], no follow-up study using this collection has been reported. A Chinese soybean core collection has been established from 23,857 collections based on agronomic traits and SSR markers [16]. Cho et al. [24] have reported a core collection of Korean landraces selected from around 7,000 collections. Oliveira et al. [6] have developed a soybean core collection of the USDA Soybean Germplasm Collection by using several sampling strategies. A Japanese mini-core collection consisting of 96 accessions has been selected by using 191 SNP markers from 1,603 collections [17]. Priolli et al. [25] have developed a core collection with 31 accessions from 435 Brazilian soybean cultivars using 27 SSR markers. Kuroda et al. [22] have developed two core collections from 1,359 wild and cultivated soybean collections using information of 20 SSR markers. Likewise, genetic diversity analyses have been performed to construct core collections in soybean by using phenotypic information or a small number of SSR markers. In recent years, genetic markers have been increasingly developed and applied to many crop species, including barley [26], wheat [27], common bean [28], maize [29], sunflower [30], apple [31], tomato [32], and soybean [33].

The objectives of this study were to: 1) construct a Korean soybean core collection from 2,872 soybean collections, 2) evaluate genetic and phenotypic diversity of core and total collections, and 3) conduct genome-wide association studies (GWAS) for several important phenotypic traits in both core and total collections and compare GWAS results.

Material and methods

Plant materials

A total of 2,872 soybean collections (total collection) maintained in the National Agrobiodiversity Center in the Rural Development Administration (Jeonju, S. Korea) were used for this study. These soybean collections originated from China, Japan, the USA, N. Korea, S. Korea, and other countries (Canada, France, Sweden, and unknown). Among these collections, 2,872 cultivated soybeans originated from China, Japan, the USA, N. Korea, S. Korea, and others (76, 41, 79, 108, 2,556, and 12 accessions, respectively).

DNA extraction and SNP genotype analysis

Young trifoliate leaves from a single plant of each collection were collected for DNA extraction. Genomic DNA was isolated using the CTAB (cetyltrimethylammonium bromide) method [34] with a minor modification. Axiom® 180k (180,961) SoyaSNP array was used for SNP genotyping [35]. Genomic DNAs from the total collection were hybridized to arrays using Affymetrix GeneTitan system based on the manufacturer’s instructions. SNP genotyping was conducted based on Axiom® Genotyping Solution Data Analysis User Guide [35].

Genetic diversity analysis of the total collection

Phylogenetic analysis of the total collection was performed using genotype data from the Axiom® 180k SoyaSNP array. Phylogenetic tree and principal component analysis (PCA) were conducted by using APE R package [36]. A distance matrix was calculated from pairwise DNA sequences of samples. Phylogenetic tree was constructed based on the distance matrix using a neighbor-joining algorithm. The distance matrix was also used for PCA. To estimate population structure of a core collection, fastSTRUCTURE [37] was used. This program has been widely used to calculate posterior inference based on the Bayesian framework [38, 39]. SNP genotype data set from 2,872 soybean accessions were converted to variant call format (VCF) format and fastSTRUCTURE was executed for various numbers of populations (K = 1,…,10). Numbers of subpopulations were defined using the marginal likelihood function. Arlequin v 3.5 [40] program was used to compare molecular diversity between and within subpopulations and geographic regional groups using analysis of molecular variance (AMOVA) [41, 42].

Construction of a core collection

To select accessions for a core collection, genotype data of the total collection were analyzed using GenoCore [43], a fast and consistent method for large datasets. GenoCore offers two options: coverage and delta. The coverage option provides a percentage value representing how many of the selected subsets of samples reflect the diversity of the total collection while the delta value represents increasing ratio in the coverage. In the present study, 99% of the coverage and 0.001% of the delta value were applied.

Evaluation of phenotypic traits

Five qualitative traits and two quantitative traits were evaluated for both core and total collections in field tests. The total collection was tested in an experimental field at the National Institute of Crop Science (Suwon, S. Korea, 37°15’29.1” N, 126°58’34.3” E) in 2013 and the core collection was grown in an experimental field at the National Institute of Crop Science (Wanju, S. Korea, 35°50’35.2” N, 127°02’45.9” E) in 2016. The experimental design was a randomized block design (RBD) with two replicates. The field plot was 3 m long with 0.75 m of row spacing and plant spacing within rows of 0.15 m. Two seeds were sown per spot (or hole).

Five qualitative traits, including hypocotyl color, leaflet shape, flower color, growth habit, and pubescent color, were investigated at V1, R1, R2, R3, and R7 stages, respectively. Flower color was recorded as purple or white. Hypocotyl color was purple or green. Pubescent color was tan or gray. Growth habit was recorded as either determinate or indeterminate type. Leaflet shape was classified as ovate or narrow. In the case of quantitative traits, days to flowering (R1) were recorded as numbers of days from planting to beginning to bloom. Each accession was individually bulk-harvested after its full maturity stage (R8). A sample of 100 cleaned seeds from each collection was randomly selected and weighed.

Genome-wide association study

To evaluate whether the core collection could sufficiently represent genetic diversity of the total collection, four phenotypic traits (days to flowering, flower color, pubescent color, and growth habit) with known controlling genes were tested through GWAS. After filtration of SNPs with < 1% minor allele frequency (MAF) and > 5% heterozygosity, 131,620 SNPs were selected and used for GWAS using a compressed mixed linear model [44]. All analyses were conducted using GAPIT [45] and GenABEL package [46] in R Project (http://www.r-project.org/). GWAS was conducted for both core and total collections to compare their genetic diversities for tested phenotypic traits.

Genes controlling flower color (W1), pubescent color (T), and growth habit (TFL1) of soybean have been mapped to chromosome 13, 6, and 19, respectively [47, 48, 49]. Quantitative trait loci (QTL) associated with days to flowering, named E1, E2, and E3, have been genetically mapped to chromosome 6, 10, and 19 in soybean, respectively [50, 51, 52, 53]. Above-listed traits were analyzed for both core and total collections.

Results

Soybean germplasm genotyping

To identify genetic diversity in the total collection, 170,223 high-quality SNPs from the Axiom® 180k SoyaSNP array were used. In the total collection, average sample call rate, homozygosity rate, and heterozygosity rate were 99.41%, 98.82%, and 0.59%, respectively. In the core collection, average sample call rate, homozygosity rate, and heterozygosity rate were 98.83%, 97.21%, and 1.62%, respectively (Table 1). Soybean accessions with call rates less than 97% were excluded from future analysis (Affymetrix, Analysis Guide Axiom® Genotyping solution Data Analysis Guide).

Tab. 1. Genotypic diversity of 2,872 accessions in the total collection and 430 accessions in the core collection by Axiom® 180K SoyaSNP array.
Genotypic diversity of 2,872 accessions in the total collection and 430 accessions in the core collection by Axiom<sup>®</sup> 180K SoyaSNP array.

Construction of the core collection

The 170,223 high-quality SNPs were also used to construct the core collection from the total collection (S1 Table). Based on genotyping data, accessions with more than 99.9% similarity or with SNPs ≤ 1% of MAF were excluded from further analysis.

We first selected 407 accessions reflecting 99% of genetic diversity of total accessions. Second, as additional sources of genetic variations, 24 cultivated soybeans used for whole-genome re-sequencing were included with the first selection. This addition produced a core collection consisting of 431 soybean accessions. However, we found that one collection (Danbaek) was included in both the first selection and the re-sequencing accessions. Finally, a Korean soybean core collection consisting of 430 accessions was constructed (Table 2 and S2 Table).

Tab. 2. Construction of the core collection by using Axiom® 180K SoyaSNP array and GenoCore.
Construction of the core collection by using Axiom<sup>®</sup> 180K SoyaSNP array and GenoCore.

PCA was conducted to confirm whether the core collection represented the genetic diversity of the total soybean collection. Results showed that the distribution of collections by genetic diversity showed similar trends for both core and total accessions (Fig 1). As shown in the scree plot, genetic variations of 72 accessions reflected more than 95% of genetic variations in the total collection (Fig 2). In addition, 407 accessions reflected more than 99% of genetic variations of the total collection (S2 Table). The number of accessions in the core collection accounted for about 15% of the total collection. Accessions in the core collection originated from South Korea (319, 74%), North Korea (41, 9%), China (42, 10%), Japan (13, 3%), the USA (12, 3%), and other countries (3, 1%) (Table 2).

Fig. 1.
Principal component analysis of the core collection (red triangle) and the total collection (black circle).
Principal component analysis model statistics for the core collection.
Fig. 2. Principal component analysis model statistics for the core collection.
The core collection was made of 407 accessions, reflecting more than 99% of the genetic diversity of the total collection. Specifically, 72 accessions reflected more than 95% of the genetic variations of the total collection (black arrow).

Population structure and phylogenetic analysis

SNP information was used to evaluate the population structure of the core collection. Marginal likelihood was found to be the highest when K was 5 (Fig 3A). The number of accessions belonging to each subpopulation was as follows: subpopulation I, 84; subpopulation II, 89; subpopulation III, 111; subpopulation IV, 93; and subpopulation V, 53. AMOVA indicated that only 2.01% of the molecular variance was explained by the origin while 16.18% was contributed by the subpopulation. The variation among individuals within populations was the highest (81.82%, Table 3). Accessions from South Korea and North Korea were distributed across five subpopulations while 30 of 42 collections from China belonged to subpopulation V. Most of the improved lines belonged to subpopulations I and II while most of landraces belonged to subpopulations III and IV.

Population structure of the core collection.
Fig. 3. Population structure of the core collection.
(A) The number of populations (K = 5) was selected using marginal likelihood provided by fastSTRUCTURE. (B) Five subpopulation clusters inferred by fastSTRUCTURE are indicated by different colors.
Tab. 3. Analysis of molecular variance of geographic groups and five subpopulations in the Korean core collection.
Analysis of molecular variance of geographic groups and five subpopulations in the Korean core collection.

The population structure of the core collection was doubled-checked using both phylogenetic tree and PCA. In the phylogenetic tree, most accessions with more than 50% admixture fractions for a genome were distributed over the five subpopulations while admixed accessions were placed between subpopulation clades (Fig 4). Scree plot showed that five subpopulations were optimal (Fig 5). PCA plots indicated that accessions from the South and North Korea were evenly distributed and 28 out of 42 Chinese accessions (67%) were grouped together (Fig 6). As a result, although most of the resources in the Korean soybean core collection were Korean soybean resources and evenly distributed, the variation of Korean resources did not cover variation in Chinese resources.

Phylogenetic tree of the core collection.
Fig. 4. Phylogenetic tree of the core collection.
The core collection was divided into five subpopulation clades (I-V).
Scree plot of principal component analysis in the core collection between variance and number of principal components.
Fig. 5. Scree plot of principal component analysis in the core collection between variance and number of principal components.
Principal component analysis of the Korean core collection.
Fig. 6. Principal component analysis of the Korean core collection.
CA: Canada; CN: China; JP: Japan; NK: North Korea; SD: Sweden; SK: South Korea; USA: United States of America.

Phenotypic variations in the core collection

Important phenotypic traits of soybean in both core and total collections were evaluated. In the core collection, days to flowering ranged from 27 to 68 days (Fig 7B). Two hundred ninety-seven accessions had purple color hypocotyls and flowers while 133 had green hypocotyls and white flowers. For pubescent color, 210 accessions were tan and 220 were gray. Regarding growth habit, 305 accessions presented a determinate growth type while 121 accessions had an indeterminate type. Growth habits of four accessions were not noted. Most accessions (412) had ovate leaflet shapes while only 18 had narrow leaflet shapes (Fig 7A and S3 Table). The 100-seed weight of these accessions ranged from 7.3 to 59.4 g (Fig 7C).

Morphological trait measurements for the core collection and the total collection.
Fig. 7. Morphological trait measurements for the core collection and the total collection.
(A) Flower color, hypocotyl color, pubescent color, growth habit, and leaflet shape in the core collection and total collection. (B) Days to flowering in the core collection and total collection. (C) The 100-seed weight in the core collection and the total collection.

In the present study, we investigated and compared five phenotypes in core and total collections. The proportion of each trait (flower color, hypocotyl color, pubescent color, growth habit, and leaflet shape) in the core collection to that in the total collection was approximately 18.1%. Specifically, proportions of these qualitative traits in the core collection to that in the total collection were around 15.8% in South Korea, 38% in North Korea, 51.4% in China, 27.7% in Japan, 16.6% in the USA, and 20.2% in other countries (Fig 7A and S3 Table). These results indicated that phenotypic proportions of the tested phenotypic traits in the core collection were similar to those in the total collection, although the core collection was selected based on SNP genotypes.

In the case of quantitative traits, the proportion of days to flowering in the core collection to that in the total collection was around 26.1%. The proportion of days to flowering in the core collection to that in the total collection was around 12.1% in S. Korea, 48.9% in N. Korea, 68.6% in China, 47.8% in Japan, 19.7% in the USA, and 39.3% in other countries (Fig 7B and S4 Table). The proportion of the 100-seed weight in the core collection to that in the total collection was about 25.2%, 23.7% in S. Korea, 33.3% in N. Korea, 77.9% in China, 37.5% in Japan, 24.1% in the USA, and 100% in other countries (Fig 7C and S5 Table). The 100-seed weight in the total collection ranged from 6.2 g to 59.4g. Similarly, the 100-seed weight in the core collection ranged from 7.3 g to 59.4g. Days to flowering ranged from 21 to 78 days in the total collection and from 27 to 68 days in the core collection.

We also compared the distribution of accessions between total and core collections by their origins (S1, S2, and S3 Figs). The distribution trend of the accessions by origin in the core collection was similar to that in the total collection. In particular, the percentage of the core collection extracted from the total collection of Chinese resources was the highest (qualitative traits, 51.4%; days to flowering trait, 68.6%; 100-seed weight trait, 77.9% for Chinese resources extracted from the total collection) (S1C, S2C, and S3C Figs). These results might indicate that the genetic diversity of Chinese collections was higher than that of collections from other countries.

Genome-wide association study

Results of GWAS for the four important phenotypic traits in core and total collections were compared. For flower color (W1), pubescent color (T), and growth habit (TFL1), GWAS confirmed that both collections had significant SNPs at gene regions of the same positions. QTLs associated with days to flowering genes E1, E2, and E3 were positioned on the same chromosomes in both core and total collections (Fig 8).

GWAS for days to flowering in the core collection and the total collection.
Fig. 8. GWAS for days to flowering in the core collection and the total collection.
(A) Manhattan plot of the core collection. (B) Manhattan plot of the total collection. (C) Quantile-Quantile plot for days to flowering in the core collection. (D) Quantile-Quantile plot for days to flowering in the total collection.

Discussion

In the present study, a Korean soybean core collection was constructed based on SNP information from the Axiom® 180k SoyaSNP array. Genotypic, phenotypic diversity, and population structure of the population were then analyzed. As mentioned previously, earlier studies have tested various selection methods to construct core collections in diverse crop species. To date, most of core collections for crop species have been developed based on geographic origin, morphological, and phenotypic traits. Ultimate goals of a core collection are not only to have the maximum possible genetic diversity of the total collection, but also to have minimum repetition. However, it has been reported that phenotypic data cannot perfectly reflect the genetic diversity of the total collection since most of quantitative agronomic traits are easily affected by environmental factors [54]. In contrast, molecular markers can directly reflect the genetic diversity in a total collection at DNA sequence level [55].

Although several studies have previously established core collections based on genotypic data, relatively few molecular markers have been used to determine genetic diversity in total collections and select core collections. Kuroda et al. [22] have used 20 SSR markers to develop two core collections of wild and cultivated soybean. To develop a core set from 435 Brazilian soybean cultivars, Priolli et al. [25] have used 130 alleles from 27 SSR loci. Xu et al. [13] have used 24 SSR markers to construct the Chinese jujube (Ziziphus jujube Mill) core collection. For Capsicum germplasm, Lee et al. [12] have constructed a core collection ‘CC240’ consisting of 240 accessions by using 48 SNP markers and 32 phenotypic/morphological traits. The number of SNP markers (170,223) used to construct the core collection in the present study was significantly greater than that used in previously reported studies. It sufficiently covered all genomic regions of soybean. In addition to high-density molecular markers, the use of precisely tested phenotypic and agronomic traits can improve the quality of selection strategies used to construct core collections in crop species. After the development of a core collection, accurate evaluation of important and target phenotypes of the core collection in multiple environments will accelerate the course of gene discovery and marker development.

Compared to STRUCTURE analysis, the phylogenetic tree did not present clear separation according to geographical origin. In particular, most of the Chinese collections belonged to subpopulation V (S2 Table). However, collections from other countries were evenly distributed across subpopulations in both the phylogenetic tree and the population structure analysis. AMOVA indicated that genetic variability of populations within groups was higher than that between groups. Although few accessions from other countries were included in the core collection and the number of accessions from Korea was larger than that from other countries, the genetic diversity of Korean accessions was significantly higher than that of accessions from other countries. In other words, there is genetic distance between resources from Korean and other countries. If a new core collection was constructed by integrating the core collection developed from China, Japan, USA, and Korea, genetic distances between the collections from each origin would be clearly elucidated.

The ultimate aim of constructing the Korean soybean core collection in this study was to improve the utilization of useful breeding materials for genome prediction or selection in soybean breeding programs as well as to explore candidate genes related to particular traits through GWAS. In the present study, GWAS was performed to confirm whether the core collection represented genetic diversity of the total collection. Monogenic traits including flower color (W1), pubescent color (T), and growth habit (TFL1) were identically detected in both collections. However, the GWAS for polygenic traits presented different results for the two collections. In the case of E1 (major gene for days to flowering), GWAS results were identical for the two collections. The p-value of the SNP associated with E1 was highly significant in the two collections: p = 1.83E-08 in the core collection and p = 1.21E-11 in the total collection. In the case of minor genes (E2 and E3), GWAS results for E3 were similar for the two collections. The p-value of the SNP representing E3 was also highly significant in the two populations: p = 3.67E-07 in the core collection and p = 5.40E-10 in the total collection. However, the p-value of the SNP representing E2 was found to be highly significant (p = 3.81E-12) in the total collection, but less significant (p = 0.0021) in the core collection (S6 Table). These results indicated that the size of the populations might be an important factor in GWAS for minor QTLs and that the Korean core collection was suitable for GWAS.

Approximately 10% of the total collection has been recommended for the size of a core collection [1]. Most of the core collections surveyed by Spillane et al. (unpublished) were also 5–20% of the size of the total collection [56]. A small-sized population might not identify QTLs with minor effects. However, the size of the core collections could be determined by genetic diversity of the total population. If the genetic diversity of the total population is low, the size of the core collection does not need to be large. In contrast, the size of the extracted core collection should be large if the genetic diversity of the total population is high. Therefore, we strongly recommend extracting core collections from large accessions with extensive genetic diversity. After developing core collections for crop species, utilization of the collections will become more and more important.

In conclusion, in the work, we constructed a Korean soybean core collection representing genetic and phenotypic diversity of the total collection based on the 180K SNP genotypes. Several genetic studies and breeding programs using core collection are currently ongoing to identify tolerances to biotic and/or abiotic stresses, find candidate genes for important agricultural traits, develop molecular markers for these traits, and develop marker sets for genome selection. The Korean soybean core collection can provide useful genetic material for soybean diversity and genetic studies as well as effective information for soybean breeding programs. Detailed information on soybean accessions included in this Korean core collection is available from the corresponding authors. It can also be obtained from the National Agrobiodiversity Center, Rural Development Administration through its website at http://genebank.rda.go.kr.

Supporting information

S1 Fig [tif]
Comparison of the five morphological traits in the core collection and the total collection by origin.

S2 Fig [tif]
Comparison of days to flowering between the core collection and the total collection by origin.

S3 Fig [tif]
Comparison of the 100-seed weight between the core collection and the total collection by origin.

S1 Table [xlsx]
List of the total collection.

S2 Table [xlsx]
List of the core collection.

S3 Table [xlsx]
Morphological traits investigated in the core collection and total collection.

S4 Table [xlsx]
Investigation of days to flowering in the core collection and the total collection.

S5 Table [xlsx]
Investigation of 100-seed weight in the core collection and the total collection.

S6 Table [xlsx]
Information on SNPs linked to days to flowering (, , and ) from GWAS.


Zdroje

1. Brown AHD. Core collections: a practical approach to genetic resources management. Genome. 1989; 31: 818–24.

2. Upadhyaya HD. Establishing core collections for enhanced use of germplasm in crop improvement. Ekin J Crop Breed and Gen. 2015; 1–1: 1–12.

3. Bhandari HR, Bhanu AN, Srivastava K, Singh MN, Shreya, Hemantaranjan A. Assessment of genetic diversity in crop plants–An overview. Advances in Plants & Agriculture Research. 2017; 7(3): 00255.

4. Guo Y, Li Y, Hong H, Qiu LJ. Establishment of the integrated applied core collection and its comparison with mini core collection in soybean (Glycine Max). The Crop Journal. 2014. pp. 38–45.

5. Marshall DR. Limitations to the use of germplasm collections. In: Brown AHD, Frankel OH, Marshall ER, Williams JT, editors. The use of plant genetic resources. New York: Cambridge University Press; 1989. pp. 105–120.

6. Oliveira MF, Nelson RL, Geraldi IO, Cruz CD, Toledo JFF. Establishing a soybean germplasm core collection. Field Crops Research. 2010; 119: 227–289.

7. Odong TL, Jansen J, Eeuwijk FA, Hintum TJL. Quality of core collection for effective utilization of genetic resources review, discussion and interpretation. Theor Appl Genet. 2013; 126: 289–305. doi: 10.1007/s00122-012-1971-y 22983567

8. Li ZC, Zhang HI, Cao YS, Qiu ZE, Wei XH, Tang SX, et al. Studies on the sampling strategy for primary core collection of Chinese ingenious rice. Acta Agron Sin. 2003; 29: 20–24.

9. Dong YC, Cao YS, Zhang SC, Wang LF, You GX, Pang BS, et al. Establishment of candidate core collections in Chinese common wheat germplasm. J Plant Genet Resour. 2003; 4: 1–8.

10. Xu H, Mei Y, Hu J, Zhu J, Gong P. Sampling a core collection of island cotton (Gossypiumbar badense L.) based on the genotypic values of fiver traits. Resour Crop Evol. 2006; 53: 515–521.

11. Holbook CC, Anderson WF. Evaluation of a core collection to identify resistance to late leafspot in peanut. Crop Sci. 1995; 35: 1700–1702.

12. Lee HY, Ro NY, Jeong HJ, Kwon JK, Jo J, Ha Y, et al. Genetic diversity and population structure analysis to construct a core collection from a large Capsicum germplasm. BMC Genetics. 2016; 17: 142. doi: 10.1186/s12863-016-0452-8 27842492

13. Xu C, Gao J, Du Z, Li D, Wang Z, Li Y, et al. Identifying the genetic diversity, genetic structure and a core collection of Ziziphus jujube Mill. Var. jujube accessions using microsatellite markers. Nature Scientific Reports. 2016; 6: 31503.

14. Hu J, Wang P, Su Y, Wang R, Li Q, Sun K. Microsatellite diversity, population structure, and core collection formation in Melon germplasm. Plant Mol Biol Rep. 2015; 33: 439–447.

15. Diwan N, McIntosh MS, Bauchan GR. Methods of developing a core collection of annual Medicago species. Theor Appl Genet. 1995; 90: 755–761. doi: 10.1007/BF00222008 24172915

16. Wang LX, Guan Y, Guan RX, Li YH, Ma YS, Dong ZM, et al. Establishment of Chinese soybean (Glycine max) core collections with agronomic traits and SSR markers. Euphytica. 2006; 151: 215–223.

17. Kaga A, Shimizu T, Watanabe S, Tsubokura Y, Katoyose Y, Harada K, et al. Evaluation of soybean germplasm conserved in NIAS genebank and development of mini core collections. Breeding Science. 2012; 61: 566–592. doi: 10.1270/jsbbs.61.566 23136496

18. Qiu LJ, Xing LL, Guo Y, Wang J, Jackson SA, Change RZ. A platform for soybean molecular breeding: the utilization of core collections for food security. Plant mol. Biol. 2013; 83: 41–50. doi: 10.1007/s11103-013-0076-6 23708950

19. Chang RZ, Qiu J, Sun J, Chen Y, Li X, Xu Z. Collection and conservation of soybean germplasm in China. In: Proc. World Soybean Research Conference VI, Chicago, IL, 4–7, August 1999. National Soybean Research Lab, Urbana, 1999. pp. 172–176.

20. Carter TE, Nelson RL, Sneller CH, Cui Z. Genetic diversity in soybean. In: Boerma HR, Specht JE. (Eds), Soybeans: Improvement, Production and Uses, vol. 16, 3rd ed. American Society of Agronomy, Madison, 2004. pp. 303–416.

21. Zhao L, Dong Y, Liu B, Hao S, Wang K, Li X. Establishment of a core collection for the Chinese annual wild soybean (Glycine soja). Chinese Science Bulletin. 2005; 50: 989–996.

22. Kuroda Y, Tomooka N, Kaga A, Wanigadeva SMSW, Vaughan DA. Genetic diversity of wild soybean (Glycine soja Sieb. Et Zucc.) and Japanese cultivated soybeans [G. max (L.) Merr.] based on microsatellite (SSR) analysis and the selection of a core collection. Cenet Resour Crop Evol. 2009; 56: 1045–1055.

23. Brown AHD, Grace JP, Speer SS. Designation of a core collection of perennial Glycine. Soybean Genetics Newsletter. 1987; 14: 59–70.

24. Cho GT, Yoon MS, Lee J, Baek HJ, Kang JH, Kim TS, et al. Development of a core set of Korean soybean landraces [Glycine max (L.) Merr.]. J Crop Sci Biotech. 2008; 11: 157–162.

25. Priolli RHG, Wysmierski PT, Cunha CP, Pinheiro JB, Vello NA. Genetic structure and a selected core set of Brazilian soybean cultivars. Genetics and Molecular Biology. 2013; 36: 382–390. doi: 10.1590/S1415-47572013005000034 24130446

26. Close TJ, Bhat PR, Lonardi S, Wu Y, Rostoks N, Ramsay L. Development and implementation of high-throughput SNP genotyping in barley. BMC Genomics. 2009; 10: 582. doi: 10.1186/1471-2164-10-582 19961604

27. Akhunov E, Nicolet C, Dvorak J. Single nucleotide polymorphism genotyping in polyploidy wheat with the Illumina GoldenGate assay. Theor Appl Genet. 2009; 119: 507–517. doi: 10.1007/s00122-009-1059-5 19449174

28. Hyten DL, Song QJ, Fickus EW, Quigley CV, Lim JS, et al. High throughput SNP discovery and assay development in common bean. BMC Genomics. 2010; 11: 475. doi: 10.1186/1471-2164-11-475 20712881

29. Ganal MW, Durstewitz G, Polley A, Be´rard A, Buckler ES, et al. A large maize (Zea mays L.) SNP genotyping array: Development and germplasm genotyping, and genetic mapping to compare with the B73 reference genome. PLoS ONE. 2011; 6: e28334. doi: 10.1371/journal.pone.0028334 22174790

30. Bachlava E, Taylor CA, Tang S, Bowers JE, Mandel JR, et al. SNP discovery and development of a high-density genotyping array for sunflower. PLoS ONE. 2012; 7: e29814. doi: 10.1371/journal.pone.0029814 22238659

31. Chagne´ D, Crowhurst RN, Troggio M, Davey MW, Gilmore B, et al. Genome-wide SNP detection, validation, and development of an 8K SNP array for apple. PLoS ONE. 2012; 7: e31745. doi: 10.1371/journal.pone.0031745 22363718

32. Sim SC, Durstewitz G, Plieske J, Wieseke R, Ganal MW, et al. Development of a large SNP genotyping array and generation of high-density genetic maps in tomato. PLoS ONE. 2012; 7: e40563. doi: 10.1371/journal.pone.0040563 22802968

33. Song Q, Hyten DL, Jia G, Quigley CV, Fickus EW, Nelson RL, et al. Development and evaluation of SoySNP50K, a high density genotyping array for soybean. PLoS ONE. 2013; 8: e54985. doi: 10.1371/journal.pone.0054985 23372807

34. Doyle JJ, Doyle JL. Isolation of plant DNA from fresh tissue. Focus. 1990; 12: 13–15.

35. Lee YG, Jeong N, Kim JH, Lee K, Kim KH, Pirani A, et al. Development, validation and genetic analysis of a large soybean SNP genotyping array. Plant J. 2015; 81: 625–636. doi: 10.1111/tpj.12755 25641104

36. Paradis E, Claude J, Strimmer K. APE: analyses of phylogenetics and evolution in R language. Bioinformatics. 2004; 20: 289–290. doi: 10.1093/bioinformatics/btg412 14734327

37. Raj A, Stephens M, Pritchard JK. fastSTRUCTURE: Variational inference of population structure in large SNP data sets. Genetics. 2014; 197: 573–589. doi: 10.1534/genetics.114.164350 24700103

38. Logsdon BA, Hoffman GE, Mezey JG. A variational bayes algorithm for fast and accurate multiple locus genome-wide association analysis. BMC Bioinformatics. 2010; 11: 58. doi: 10.1186/1471-2105-11-58 20105321

39. Carbonetto P, Stephens M. Scalable variational inference for Bayesian variable selection in regression, and its accuracy in genetic association studies. Bayesian Anal. 2012; 7: 73–108.

40. Excoffier L, Lischer HEL. Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under linux and windows. Molecular Ecology Resources. 2010; 10: 564–567. doi: 10.1111/j.1755-0998.2010.02847.x 21565059

41. Michalakis Y, Excoffier L. A generic estimation of population subdivision using distances between alleles with special reference for microsatellite loci. Genetics. 1996; 142:1061–1064. 8849912

42. Excoffier L, Smouse PE, Quattro JM. Analysis of molecular variance inferred from metric distances among DNA haplotypes: Application to human mitochondrial DNA restriction data. Genetics. 1992; 131: 479–491. 1644282

43. Jeong S, Kim JY, Jeong SC, Kang ST, Moon JK, Kim N. GenoCore: A simple and fast algorithm for core subset selection from large genotype datasets. PLoS ONE. 2017; 12: e0181420. doi: 10.1371/journal.pone.0181420 28727806

44. Zhang Z, Ersoz E, Lai CQ, Todhunter RJ, Tiwari HK, Gore MA, et al. Mixed linear model approach adapted for genome-wide association studies. Nat Genet. 2010; 42: 355–360. doi: 10.1038/ng.546 20208535

45. Lipka AE, Tian F, Wang Q, Peiffer J, Li M, Bradbury PJ, et al. GAPIT: genome association and prediction integrated tool. Bioinformatics. 2012; 28: 2397–2399. doi: 10.1093/bioinformatics/bts444 22796960

46. Aulchenko YS, Ripke S, Isaacs A, Duijn CM. GenABEL: an R library for genome-wide association analysis. Bioinformatics. 2007; 23: 1294–1296. doi: 10.1093/bioinformatics/btm108 17384015

47. Takahashi R, Benitez ER, Funatsuki H, Ohnishi S. Soybean maturity and pubescence color genes improve chilling tolerance. Crop Sci. 2005; 45: 1387–1393.

48. Tian Z, Wang X, Lee R, Li Y, Specht JE, Nelson RL, et al. Artificial selection for determinate growth habit in soybean. Proc Natl Acad Sci. 2010; 107: 8563–8568. doi: 10.1073/pnas.1000088107 20421496

49. Yang K, Jeong N, Moon JK, Lee YH, Lee SH, Kim HM, et al. Genetic analysis of genes controlling natural variation of seed coat and flower colors in soybean. Journal of Heredity. 2010; 101: 757–768. doi: 10.1093/jhered/esq078 20584753

50. Liu B, Kanazawa A, Matsumura H, Takahashi R, Harada K, Abe J. Genetic redundancy in soybean photoresponses associated with duplication of the phytochrome A gene. Genetics. 2008; 180: 995–1007. doi: 10.1534/genetics.108.092742 18780733

51. Watanabe S, Hideshima R, Xia Z, Tsubokura Y, Sato S, Nakamoto Y, et al. Map-based cloning of the gene associated with the soybean maturity locus E3. Genetics. 2009; 182: 1251–1262. doi: 10.1534/genetics.108.098772 19474204

52. Watanabe S, Xia Z, Hideshima R, Tsubokura Y, Sato S, Yamanaka N, et al. A map-based cloning strategy employing a residual heterozygous line reveals that the GIGANTEA gene is involved in soybean maturity and flowering. Genetics. 2011; 188:395–407. doi: 10.1534/genetics.110.125062 21406680

53. Xu M, Yamagishi N, Zhao C, Takeshima R, Kasai M, Watanabe S, et al. The soybean-specific maturity gene E1 family of floral repressors controls night-break responses through down-regulation of FLOWERING LOCUS T orthologs. Plant Physiology. 2015; 168:1735–1746. doi: 10.1104/pp.15.00763 26134161

54. Hu J, Zhu J, Xu HM. Methods of constructing core collection by stepwise clustering with three sampling strategies based on the genotypic values of crops. Theor Appl Genet. 2000; 101: 264–268.

55. Liu XB, Li J, Yang ZL. Genetic diversity and structure of core collection of winter mushroom (Flammulina velutipes) developed by genomic SSR markers. Hereditas. 2018; 155: 3. doi: 10.1186/s41065-017-0038-0 28690478

56. van Hintum ThJL, Brown AHD, Spillane C, Hodgkin T. Core collections of plant genetic resources. IPGRI Technical Bulletin No. 3. International Plant Genetic Resources Institute, Rome, Italy. 2000.


Článek vyšel v časopise

PLOS One


2019 Číslo 10
Nejčtenější tento týden
Nejčtenější v tomto čísle
Kurzy Podcasty Doporučená témata Časopisy
Přihlášení
Zapomenuté heslo

Zadejte e-mailovou adresu, se kterou jste vytvářel(a) účet, budou Vám na ni zaslány informace k nastavení nového hesla.

Přihlášení

Nemáte účet?  Registrujte se

#ADS_BOTTOM_SCRIPTS#