Debjani Roy Choudhury
Ramesh Kumar
Vimala Devi S2
Kuldeep Singh
Rakesh Singh- 1Division of Genomic Resources, NBPGR, New Delhi, India
- 2Division of Germplasm Conservation, NBPGR, New Delhi, India
- 3NBPGR, New Delhi, India
- 4NIPB, New Delhi, India
In India, rice (Oryza sativa L.) is cultivated under a variety of climatic conditions. Due to the fragility of the coastal ecosystem, rice farming in these areas has lagged behind. Salinity coupled with floods has added to this trend. Hence, to prevent genetic erosion, conserving and characterizing the coastal rice, is the need of the hour. This work accessed the genetic variation and population structure among 2,242 rice accessions originating from India’s east coast comprising Andhra Pradesh, Orissa, and Tamil Nadu, using 36 SNP markers, and have generated a core set (247 accessions) as well as a mini-core set (30 accessions) of rice germplasm. All the 36 SNP loci were biallelic and 72 alleles found with average two alleles per locus. The genetic relatedness of the total collection was inferred using the un-rooted neighbor-joining tree, which grouped all the genotypes (2,242) into three major clusters. Two groups were obtained with a core set and three groups obtained with a mini core set. The mean PIC value of total collection was 0.24, and those of the core collection and mini core collection were 0.27 and 0.32, respectively. The mean heterozygosity and gene diversity of the overall collection were 0.07 and 0.29, respectively, and the core set and mini core set revealed 0.12 and 0.34, 0.20 and 0.40 values, respectively, representing 99% of distinctiveness in the core and mini core sets. Population structure analysis showed maximum population at K = 4 for total collection and core collection. Accessions were distributed according to their population structure confirmed by PCoA and AMOVA analysis. The identified small and diverse core set panel will be useful in allele mining for biotic and abiotic traits and managing the genetic diversity of the coastal rice collection. Validation of the 36-plex SNP assay was done by comparing the genetic diversity parameters across two different rice core collections, i.e., east coast and northeast rice collection. The same set of SNP markers was found very effective in deciphering diversity at different genetic parameters in both the collections; hence, these marker sets can be utilized for core development and diversity analysis studies.
Introduction
Rice (Oryza sativa L.) is an important cereal crop which is a predominant food for over three billion people across the globe. Since it can adapt to an extensive spectrum of environmental conditions, it is therefore considered a varied crop species (Chang, 1976). In general, more genetically diverse crops will have an increased capacity to adapt with climatic conditions, while uniformity reduces the genetic diversity (Luan et al., 2006). There should be more significant variation in the breeding population so as to stimulate genetic gain from selection to improved yield, biotic and abiotic resistance to stress, and other traits. Genebanks with proper coordination can be explored to find genetic variations present in its accessions which can further help in finding advanced traits and selecting better parental combinations for developing lineage with maximum genetic variability (Barrett and Kidwell, 1998). Genetic diversity and population structure knowledge form the backbone in building core sets adequately representing variations found in the whole collection, and thereby making the collection small and condensed (Yan et al., 2007; Agrama et al., 2009; El Bakkali et al., 2013), thus creating standard data and calculating the potential loss of genetic diversity during conservation and management (Reif et al., 2005). Trait-specific germplasm identification from large collections is crucial to successful introduction of new diversity into the crop improvement programs (Upadhyaya et al., 2014). It is pivotal to comprehend the scope of genetic variation in crop germplasm and how to control it; results from molecular marker-based diversity studies should be used with caution, especially when it comes to germplasm conservation initiatives, the reason being adaptability, which plays a major role during the process of evolution and individuals’ survival in populations. As a result, it becomes ambiguous, determining whether plant selection is based on markers directly or on linked loci responsible for adaptive traits (Raybould et al., 1996).
Single-nucleotide polymorphisms (SNPs) are supposed to be the most bountiful variation found across the genome and therefore are excellent for high-resolution genotyping, in turn useful for analysis of genetic diversity and association mapping, linkage mapping, and marker-assisted selection (MAS) (Gonzaga et al., 2015). Simple sequence repeats (SSRs) are being replaced by SNPs at a very high degree. They have become a choice for utilization in plant breeding and genetics (McCouch et al., 2010). Using second-generation sequencing methods, scientists have been able to detect millions of SNPs in rice genome (Huang et al., 2009, 2010; Xu et al., 2011). Array-based SNP detection is currently one of the most popular high-throughput marker detection methods, allowing data to be evaluated in real time. Because of their abundance, locus specificity, low error rates, and co-dominant inheritance SNP assays are popular markers (Rafalski, 2002; Schlotterer, 2004). There are numerous multiplex and uniplex genotyping platforms are available (Syvänen, 2001; Chen and Sullivan, 2003). SNP genotyping platforms including Illumina BeadXpress (Yamamoto et al., 2010; Thomson et al., 2012), Affymetrix (Singh et al., 2015; McCouch et al., 2016), and the KASP marker system (Cheon et al., 2018; Yang et al., 2019) have been developed recently and applied to rice molecular breeding. In the year 2020, Seo et al. (2020) developed two 96-plex indica-japonica SNP genotyping sets using the Fluidigm platform; however, genotyping platforms are very expensive.
Plant Genetic Resources of native rice needs to be conserved promptly. A major conservation strategy counts on developing core collections, thereby screening gene bank accessions for important agronomic traits. To create a core collection, one needs a good construction strategy. Frankel (1984) and Brown (1989) suggested the concept of core collection; it is a small set of accessions from the entire collection with least amount of repetition and maximal genetic diversity of a species. It is usually 5% to 10% representation of the total population. Van Hintum et al. (2000), in the year 2000, outlined a comprehensive process for assembling a core collection. These can be summarized as 1) determine the overall sampling ratio; 2) partition them into groups; 3) decide on the percentage of the group that will be sampled; and 4) choose from each group entry. Several other strategies have been proposed to create a core collection such as PowerCore (Kim et al., 2007), MStrat (Gouesnard et al., 2001), stepwise clustering (Hu et al., 2000), and least distance stepwise clustering (Wang et al., 2007). These various methodologies are dependent on several parameters such as species’ genetic diversity, collection size, grouping, and data type (i.e., phenotypic or molecular data). However, if adequately characterized, a quality core collection can form a progressive collection for long-term conservation.
In the present study, we demonstrate the identification of a diverse core set panel from the east coast region of India comprising 247 core accessions and a mini core set of 30 accessions to enable elite gene mining for breeding and conservation and on a long-term perspective to form a panel for association studies. We have also done a comparative analysis of genetic diversity parameters across two different rice collections (northeast rice collection by Roy Choudhury et al., in 2014, and the coastal rice collection), thus validating the use of SNP markers for studying diversity parameters across any other rice collection.
Materials and methods
Germplasm resources
A set of 2,242 seed samples of the east coast region of India (Andhra Pradesh, Orissa, and Tamil Nadu) were procured from the Indian National Genebank, National Bureau of Plant Genetic Resources (NBPGR), New Delhi, with passport data (National ID, i.e., Indigenous Collection (IC) number) and the states to which they belong. This is given in Supplementary Table S1.
Genomic DNA isolation and molecular characterization using SNP markers
The seed was de-husked, and genomic DNA was isolated using the QIAGEN DNeasy Plant Mini Kit. A tissue lyser (Retsch, Haan, Germany) and a tissue lyser adapter set (Qiagen, Hilden, Germany) was used to grind kernels into fine powder. Working stocks with 10 ng/μl of the genomic DNA samples was prepared, and 30 µl of the diluted sample was transferred to a 96-well plate to be run on the Sequenom MassARRAY which adopts the matrix-assisted laser desorption ionization-time of flight (MALDI-TOF) mass spectrometer for most authentic detection of SNPs (www.sequenom.com). The information about the genetic location of the multiform assays designed for 36 SNPs having conserved single-copy rice genes was derived from Singh et al. (2007). Sequenom Corporation (San Diego, CA, USA) created and validated a 36-plex assay with three genes per chromosome. The assay Design 3.1 program was used to build pre-amplification primers and genotyping primers, purchased and utilized for SNP validation according to the Sequenom user manual’s methodology. There were two polymerase chain reactions (PCR). The first was a normal PCR of 45 cycles with pre-amplification primers followed by removal of unincorporated dNTPs. After adjusting genotyping primers, the second PCR was an extension PCR. The extension rate was enhanced by increasing the number of cycles to 300, thus giving the highest call rates. Call rates and extension rates were also adjusted according to genotype calling algorithms. MassARRAY Typer Analyzer 3.4 Software was used to visualize SNP calling.
Genetic diversity and phylogenetic analyses
The results of SNP data were subjected to analysis using PowerMarker (V3.25) (Liu and Muse, 2005) to calculate major allele frequency, heterozygosity, gene diversity, and PIC (polymorphic information content). The genotypes’ genetic distances (Nei et al., 1983) were computed, and a neighbor joining (NJ) tree was generated and viewed in FigTree v 1.4.3 (Rambaut, 2010). To infer historical origin, software STRUCTURE V2.3.1 was used which provides clusters of related genotypes (Pritchard et al., 2000). To infer the value of genetic cluster (K), each sample was run from K = 2 to K = 10 with the admixture model and correlated allele frequency. Each K run was replicated thrice with 100,000 burn-in period and 100,000 Monte Carlo Markov Chain replicates (Evanno et al., 2005). The analysis was carried out regardless of the accessions’ geographical origin. The dataset optimal K value was obtained using program “structure harvester” (http://taylor0.biology.ucla.edu). Additional hierarchical structure analysis was done after observing additional peaks at two different K values to unmask the groups (Ambreen et al., 2018). Accessions with membership proportions >80% were considered as pure, whereas membership proportions less than 80% were judged as admixed. The software GenAlEx V6.5 (Peakall and Smouse, 2012) was used to perform principal coordinate analysis (PCoA) and analysis of molecular variance (AMOVA) between the STRUCTURE populations.
Development of core collection
PowerCore 1.0 software (Kim et al., 2007) was used to develop the core. The data set was in simple excel format, and the first column in the data set was given the name % Accessions and the subsequent columns were named NM1, NM2, and so on. After loading the data set in the software, settings were set to heuristic search with maximum possible entries. The diversity index was also mentioned in the same window of the software. There were a number of runs carried out to extract the best possible entries in the core set. Since the number of rice accessions from each state differed ranging from 1,133 accessions from Andhra Pradesh, 378 accessions from Orissa, and 731 accessions from Tamil Nadu, therefore, to avoid individual state collections taking precedence over the core, a separate core subset was developed using the rice collections from each state. Each state accession was analyzed for maximum genetic diversity parameters forming the basis of core subsets. Finally, the core accessions from all the three states were gathered together to develop a core of coastal Indian rice of 247 accessions (126 accessions from Andhra Pradesh, 45 accessions from Orissa, and 76 accessions from Tamil Nadu) (Supplementary Table S2). A mini core comprising 30 accessions was also generated using PowerCore (Supplementary Table S3). This core collection is of core collection type I (CC-I) giving a uniform representation of the original population. Here, each entry in the core set has one or more accessions that jointly make up the whole collection (Odong et al., 2013). Shannon’s diversity index and Nei’s gene diversity index were used to evaluate the diversity captured in the core and mini core collection relative to the entire collection using PowerCore. The statistical analysis of the core set and the mini core set was done using PowerMarker (V3.25) (Liu and Muse, 2005) to calculate major allele frequency, gene diversity, heterozygosity, and PIC. This has been done to give a complete depiction of diversification in the total collection and string out the best core set which can stand out as a benchmark collection.
Kinship analysis of the core collection
Kinship analysis was used to determine shared ancestry between individuals in the core collection. Tassel v 5 (Bradbury et al., 2007) was used to create a kinship matrix. All the negative values were considered as zero, and a simple bar graph was prepared using Microsoft Excel. An interactive heat map was prepared using the online available tool https://build.ngchm.net/NGCHM-web-builder, NG-CHM BUILDER: Interactive Heat Map (Ryan et al., 2020).
Validation of 36-Plex SNP assay in east coast rice collection and northeast rice collection
Genetic validation of 36-plex SNP assay (Table1) was done by comparing the genetic parameters (e.g., PIC, gene diversity, major allele frequency, and heterozygosity) in northeast rice (Roy Choudhury et al. in 2014) and east coast rice collection. A comparative analysis of these parameters across the two populations was summarized.
TABLE 1. List of SNP primers used for genotyping of 2,242 rice accessions along with gene diversity, heterozygosity, PIC, and major allele frequency.
Results
Genetic diversity of the total rice collection
Genotyping of the 2,242 rice accessions were performed using 36 SNP markers. The markers generated 72 alleles with a mean of two alleles per locus for the entire 2,242 coastal collection (Table 1). The maximum PIC was 0.37 for markers 01-608-4_C and 03-3478-1_C, and the minimum was 0.01 for marker 04-19-4_C with an average value of 0.24. The maximum and minimum heterozygosity was 0.42 for marker 04-1801-20_C and 0.009 for marker 12-1794_C, respectively, with a mean value of 0.07. Likewise, maximum and minimum gene diversity was found to be 0.49 for marker 01-608-4_C and 0.01 for marker 04-19-4_C, respectively, with an average of 0.29. The maximum major allele frequency was 0.99 for marker 04-19-4_C, and the minimum major allele frequency was observed to be 0.53 for marker 01-608-4_C with a mean value of 0.78 (Table 1).
Phylogenetic analysis of the total rice collection
Cluster analysis of 2,242 rice accessions was performed and PowerMarker (v3.25) was used to determine the dissimilarity matrix which was used to construct an unrooted phylogenetic tree using FigTree v1.4.3 (Supplementary Figure S1). Total collection got grouped into three major groups. However, no differentiation could be made between these groups based on their geographical origin with each group displaying a heterogeneous clustering of individuals.
Genetic diversity of the rice collection of the three coastal states (Andhra Pradesh, Orissa, and Tamil Nadu)
The rice collections of coastal states belonging to the east coast of India were analyzed for the genetic diversity and population structure study. The lowest PIC recorded was 0.0 for marker 04-19-4_C from the state of Orissa, and the highest was 0.37 for marker v03-3478-1_C from the state of Andhra Pradesh. Heterozygosity was observed to be 0.0 for three markers 04-19-4_C, 11-3935_C, and 12-3200-2_C from Orissa state, and the highest value observed was 0.69 for marker 04-1801-20_C from the state of Andhra Pradesh. The lowest value of genetic diversity was found to be 0.0 for marker 04-19-4_C, and the highest was 0.49 for marker 02-267_C both from Orissa. The lowest and highest major allele frequencies recorded were 0.52 for 02-267_C and 1.0 for 04-19-4_C, respectively, from Orissa state (Table 2). For the collections from three states, phylogenetic analysis was done using the neighbor joining (NJ) method, and an unrooted NJ tree was created. Rice collection from the state of Andhra Pradesh got clustered into two major groups. Group 1 had 220 accessions, and group 2 had 913 accessions (Supplementary Figure S2). There were two major groups each in case of rice collection from Orissa (Supplementary Figure S3) and Tamil Nadu (Supplementary Figure S4). The NJ tree of rice collection of Orissa exhibited 7 accessions in group 1, and 371 accessions got clustered in group 2. The NJ tree of rice collection of Tamil Nadu exhibited 5 accessions in group 1 and 726 accessions in group 2.
TABLE 2. List of genetic diversity parameters estimated for three east coast states, core, and mini core set.
Population structure of the total rice collection
The genetic link between individual rice accessions was determined using STRUCTURE, a model-based tool. Each accessions’ membership was run from K = 2 to K = 10 (Figure 1). The ultimate number of populations was determined using Structure Harvester (http://taylor0. biology.ucla.edu). The number of populations was found to be four for the entire 2,242 east coast rice collection (Figure 2). Population 1 had 128 pure and 51 admix accessions. Population 2 showed 821 pure and 107 admixed accessions, population 3 showed 157 pure and 78 admixed accessions, and population 4 showed 700 pure and 200 admixed accessions. Most of the aromatic rice accessions got grouped in population 3. Such grouping of aromatic accessions was not seen in the NJ tree. An overview of the state-wise distribution in population shows that population 1 had 60% accessions from Orissa, while in populations 2 and 3, 72% and 93% accessions were from Andhra Pradesh, respectively, fairly exhibiting the dominance of states over population. Population 4 showed around 50% accessions from Tamil Nadu (Figure 2). Likewise, the population mean value of alpha, Fst1, Fst2, Fst3, and Fst4 generated from the model-based approach and allele-freq. divergence among populations (net nucleotide distance), computed using point estimates of population obtained using the model-based approach, is given in Supplementary Table S4 and Supplementary Table S5. Fst values showed good genetic differentiation and acceptable population structure. Venn diagrams between NJ tree and population structure was constructed to find co-linearity between them. There were 21 accessions (0.93%) overlapping between population 1 of STRUCTURE and group 2 of the NJ tree. A total of 254 accessions (11%) were found to be overlapping in population 2 of STRUCTURE and group 2 of the NJ tree. Similarly, 225 accessions (10%) were common in population 3 of STRUCTURE and group 3 of NJ tree, and 754 accessions (33%) were common in population 4 of STRUCTURE and group 3 of the NJ tree. This shows less similarity being observed between groups of the NJ tree and populations of the model-based approach (Supplementary Figures S5–S8).
FIGURE 1. Estimation of population using LnP(D) derived Δk for k from 2 to 10 of a total 2,242 east coast rice collection.
FIGURE 2. Population structure bar plot of 2,242 east coastal rice accessions.
Clustering of accessions in STRUCTURE was diverse. Apart from obtaining the highest peak at K = 4, additional smaller peaks were obtained at K = 6 and K = 8 (Figure 1), implying that there are subgroups within the four major groups. As a result, an independent STRUCTURE was run with K values from 2 to 10 for all four populations obtained above. Subclustering of population 1 identified the highest peak at K = 4, giving subpopulations named subpopulation 1a, subpopulation 1b, subpopulation 1c, and subpopulation 1d (Supplementary Figures S9, S10). Subclustering of population 2 identified the highest peak at K = 6, giving subpopulations named subpopulation 2a, subpopulation 2b, subpopulation 2c, subpopulation 2d, subpopulation 2e, and subpopulation 2f (Supplementary Figures S11, S12). Subclustering of population 3 also identified the highest peak at K = 6, giving subpopulations named subpopulation 3a, subpopulation 3b, subpopulation 3c, subpopulation 3d, subpopulation 3e, and subpopulation 3f (Supplementary Figures S13, S14). Subclustering of population 4 identified the highest peak at K = 3, giving subpopulations named subpopulation 4a, subpopulation 4b, and subpopulation 4c (Supplementary Figures S15, S16). Grouping in Structure was diverse, and no dominance of states was observed in any of the population. The allocation of accessions from different states and the expected heterozygosity of the 19 subpopulations are shown in Supplementary Table S6. In this table, there are slight variations in the values of the expected heterozygosity ranging between 0.06 and 0.40 with a mean value of 0.22, indicating good genetic diversity in the subpopulation (Luo et al., 2019).
Population structure of the rice collection of three coastal states (Andhra Pradesh, Orissa, and Tamil Nadu)
Another aspect of population structure was studied to see the clustering of accessions state-wise. Population structure analysis grouped rice collection from Andhra Pradesh into four different populations (Supplementary Figure S17) Population 1 has 150 pure and 40 admix accessions, population 2 has 260 pure and 145 admix accessions, population 3 has 150 pure and 107 admix accessions, and population 4 has 179 pure and 97 admix accessions (Supplementary Figure S18). All aromatic accessions from Andhra Pradesh got grouped in population 3. Rice collection from Orissa got grouped into three populations (Supplementary Figure S19) with population 1 having 62 pure and 24 admix accessions, population 2 having 96 pure and 6 admix accessions, and population 3 having 156 pure and 34 admix accessions (Supplementary figure S20). Rice collection from Tamil Nadu got grouped into five populations (Supplementary Figure S21), population 1 having 44 pure and 20 admix accessions, 73 pure and 41 admix accessions in population 2, 168 pure and 42 admix accessions in population 3, 95 pure and 67 admix accessions in population 4, and 104 pure and 77 admix accessions in population 5 (Supplementary Figure S22). The mean values of alpha, Fst1, Fst2, Fst3, Fst4, Fst5, and Allele-freq. divergence among populations (net nucleotide distance), computed using point estimates of population generated from the model-based approach, are shown in Supplementary Table S7 and Supplementary Table S8, respectively. The values of Fst showed standardized genetic differentiation, suggesting a good population structure. Allele frequency divergence among populations computed using the point estimates of the population also gives about the genetic variation among populations. The values are indicative of the accessions being diverged in the population structure (Luo et al., 2019).
AMOVA and PCoA of total 2,242 east coast rice collection
The distribution of genetic diversity between and within the populations obtained following STRUCTURE analysis was investigated using AMOVA for total rice accessions (2242). In the first case, four populations were assumed and AMOVA analyses revealed that 29% diversity exists among populations, 20% within individuals, and 51% among individuals of the total east coast rice collection (Supplementary Table S9), while the PCoA plot showed that out of four populations obtained population 3 was getting distinctly separated from the others (Supplementary Table S10; Supplementary Figure S23). In the second case, assuming 19 subpopulations when AMOVA was done there was 30% diversity existing among populations, 21% within individuals, and 49% among individuals of the total east coast rice collection (Supplementary Table S11), while the PCoA plot showed plots with overlapping populations (Supplementary Table S12; Supplementary Figure S24). The AMOVA and PCoA studies confirmed that the actual population number in the case of the east coast collection is only four because assumption of the population of 19 did not show any advantage.
AMOVA and PCoA of the rice collection coastal states
The AMOVA study of Andhra Pradesh revealed 23% variance within individuals and 30% and 47% variance among populations and among individuals, respectively (Supplementary Table S13; Supplementary figure S25). PCoA plot analysis revealed that population 3 is distinctly isolated from the rest of the populations, and all the aromatic samples from Andhra Pradesh got grouped in to this Supplementary figure S26. Likewise, the AMOVA study of Orissa showed 19%, 68%, and 14% among populations, among individuals and within individuals, respectively (Supplementary Table S13; Supplementary figure S27). The PCoA plot showed mixing of three populations (Supplementary figure S28). The AMOVA study of Tamil Nadu showed 24%, 53%, and 23% variance among populations, among individuals, and within individuals, respectively (Supplementary Table S13; Supplementary figure S29). The PCoA plot revealed slightly isolated population 3 (Supplementary figure S30); however, there was intermixing between populations 1, 2, 4, and 5. Principal coordinate analyses of rice collection of coastal states, with percentage of variation explained by the first three axes, are summarized in Supplementary Table S14.
Generation of the core set
Out of 2,242 rice accessions studied, a core set of 247 accessions (i.e., 126 accessions from Andhra Pradesh, 45 accessions from Orissa, and 76 accessions from Tamil Nadu) and a mini core set of 30 accessions were selected using POWERCORE (Supplementary Tables S2, S3). Thirty-six SNP markers produced nine allele types, four of which were homozygous and six were heterozygous (three transitions and two transversions). We found no C/G- or G/C-type substitutions in our research. Allele frequency was determined for all three state collections and their core sets. The study of allele frequency revealed that no alleles were lost in the resulting core set and they were 99.9% similar. The same has been plotted in line plots for the total 2,242 in the east coast collection as well as for the core sets (Figure 3). These results were also in concordance with the results obtained by PowerCore where the Shannon’s diversity index and Nei’s gene diversity showed an increasing trend starting from the entire 2,242 collection to 247 core to 30 mini core collections (Figure 4). The 247 core accessions represent 11.01% of the entire collection, and the 30 mini core set represents 12.14% of the core collection (Supplementary Figure S31). Thus, we have been able to fulfill the recommended size required for a perfect core collection, which is between 5% and 20% of the original size (Bisht et al., 1998).
FIGURE 3. Graph showing the allele frequency of total east coast collection, core collection, and mini core collection.
FIGURE 4. Graph showing average diversity indices of entire collection, core collection, and mini core collection.
Genetic diversity of the core set and the mini core set
The degree of genetic diversity of the core set was studied to find the degree of genetic diversity captured from the overall coastal collection. Major allele frequency, gene diversity, heterozygosity, and PIC were observed to be 0.75, 0.34, 0.12, and 0.27 respectively (Table 2). The values of genetic diversity of the mini core set were approximately similar to the core, as shown in Table 2. A greater value genetic diversity has been observed in case of core/mini core than that of the total east coast rice collection. Hence, definitely the east coast core collection developed has rich and diverse representatives having good diversity parameters. The NJ tree showed two distinct groups in the core set and three groups in the mini core set (Figures 5, 6).
FIGURE 5. NJ tree of the 247 east coast core collection.
FIGURE 6. NJ tree of the east coast mini core set.
Population structure of core set
The population structure grouped the core set accessions into four populations (Figures 7, 8). The STRUCTURE bar plot showed population 1 having 27 pure and 20 admixed accessions, population 2 having 55 pure and 22 admixed accessions, population 3 having 20 pure and 7 admixed accessions, and population 4 having 60 pure and 36 admixed accessions (Figure 8). The mean value of alpha (0.10) for the core collection is greater than the mean value of alpha for the total collection (0.06). An alpha value close to zero means that individuals are essentially from different populations (Evanno et al., 2005). In our case, a 0.10 value of alpha in the core collection as compared to the 0.06 value of alpha in the total collection signifies more admixed individuals in the total collection. Allele-freq. divergence among populations (net nucleotide distance), computed using point estimates of population (core collection), are given in Supplementary Table S15 and Supplementary Table S16.
FIGURE 7. Estimation of population using LnP(D) derived Δk for k from 2 to 10 of the east coast core collection.
FIGURE 8. Population structure bar plot of the east coast core collection.
AMOVA and PCoA of the core set
The AMOVA study of the core collection revealed 29% variance within individuals as well as among population and 42% variance among individuals (Supplementary Table S17; Supplementary Figure S32). The PCoA plot showed populations being scattered in different quadrants (Supplementary Figure S33; Supplementary Table S18).
Kinship analysis of the core set collection
Kinship analysis of the core set showed that more than 50% of the samples had a kinship value less than zero, and less than 10% of the accessions had kinship values between 0.5 and 0.75 (Supplementary Figure S34). The kinship index and clustered heat map showed more diversity in the core collection because clustering based on this map was more heterogeneous, which indicates that maximum unique genotypes have been selected in the core collection (Figure 9).
FIGURE 9. Structured heat map showing the relatedness of rice accessions from core collection.
Validation of SNP markers in coastal rice collection and northeast rice collection
A comparative analysis of genetic diversity parameters, i.e., PIC and gene diversity between the coastal rice and northeast rice collection (Roy Choudhury et al., 2014), was done and are shown in Figure 10. On comparing the values, in both collections it was observed that in the case of the PIC, the same markers had given the maximum and minimum values in both collections (i.e., coastal rice collection and northeast rice collection). This means that markers 01-608-4_C and 03-3478-1_C had given the maximum PIC value of 0.37 in both the collections and marker 04-19-4_C had given the minimum PIC value of 0.01 in both collections. Similarly, markers 01-608-4_C and 03-3478-1_C had displayed the highest gene diversity with a value of 0.49 across both collections and a minimum value of gene diversity 0.02 with marker 04-19-4_C across both collections. Generally, line graphs for PIC and gene diversity were overlapping for both collections except at certain points where deviations were observed. For example, marker 03-1691-1_C gave a gene diversity value of 0.32 in the coastal rice collection and 0.4 in the northeast rice collection. The same marker gave a PIC value of 0.26 in the coastal rice collection and 0.32 in the northeast rice collection. Marker 11-522-1_C gave a gene diversity value of 0.45 in the coastal collection and 0.19 across the northeast collection; also, this marker 11-522-1_C gave PIC values of 0.35 and 0.17 across the coastal rice and northeast rice collections. A subsequent analysis between major allele frequency and heterozygosity from the current study and with the northeast rice collection was also evaluated (Supplementary figure S35). The values for major allele frequency and heterozygosity were overlapping except at few points where deviation was observed. For example, marker 04-1801-20_C gave heterozygosity values of 0.42 and 0.28 in coastal collection and northeast collection, respectively depicting a small amount of deviation. Similarly, marker 11-522-1_C gave major allele frequency values of 0.64 and 0.89 in the coastal collection and northeast collection, respectively, again depicting slight deviations. The validation of the same set of SNP markers (36-plex assay) in different collections (northeast and east coast collection) shows that they are very effective in deciphering the genetic diversity parameters in both the collections.
FIGURE 10. Comparison of gene diversity and PIC of east coast rice and northeast rice collection.
Discussion
The east coast collection of rice germplasm available at the National Genebank, NBPGR, New Delhi, is a valuable collection of rice accessions for the assessment of genetic diversity and other important traits. Rice accessions from different collections have served as and continued to act as sources of genes for desired qualities, contributing to the variety developments that have been reported (Choudhury et al., 2013; Das et al., 2013). Regular floods hit coastal areas and have saline soils, and various other constraints make them a fragile ecosystem with lower productivity and slow trend in growth rate (Amanullah et al., 2007). Therefore, a coastal core collection would be an appropriate measure to conserve rice in these areas for better management studies. Despite advances in genomics, the Indian rice collection has remained uncharacterized at the molecular level in terms of genetic diversity and population structure. This has been a key stumbling block in their ability to use and develop superior cultivars. In this study, effort has been made to characterize the east coast rice collection available at the National Genebank, NBPGR, New Delhi, using 36 SNP markers to enhance genome wide studies in rice. These unlinked SNP markers, which were generated and used in diverse studies, are located on the short arm, centromeric region, and long arm of all 12 rice chromosomes. As discussed, a state-wise study of the east coast collection of rice showed interesting results. A total of 72 alleles were amplified with 2 alleles per locus. The average PIC values ranged from 0.20 for Tamil Nadu, 0.21 for Orissa, and 0.25 for Andhra Pradesh. The values observed are concurrent with those observed by Singh et al. (2013) on 375 rice varieties (0.25) and Roy Choudhury et al. (2014) on the northeast rice collection (0.23) using SNP markers. Chen et al. (2017) observed PIC values of 0.27 on Ziziphus jujuba Mill, China’s most important fruit species, and 0.29 reported by Luo et al. (2019) on Camelina sativa using SNP markers. A PIC value of 0.4 was reported by Mourad et al. (2020) while they were studying the genetic diversity, population structure, and linkage disequilibrium in the spring wheat core collection, which is higher than that reported in the present study. The PIC value is generally high when SSR markers are used as observed by Rashmi et al. (2017) in 65 rice accessions (0.38) characterized using SSR markers. Pathaichindachote et al. (2019) reported a PIC value of 0.56 in 167 Thai and exotic rice varieties using 49 SSR markers and a PIC value of 0.63 reported by Jasim Aljumaili et al. (2018) on 50 aromatic rice accessions with 32 SSR markers. A mean PIC value of 0.61 has been reported by Suvi et al. (2020) while they were accessing the genetic diversity and population structure of 54 rice accessions using 14 SSR markers. Tarang et al. (2020) reported a PIC value of 0.92 with 60 microsatellite markers of 63 rice genotypes in Central and West Asia. PIC values and expected heterozygosity (He, also called gene diversity) are both indices of genetic diversity among genotypes in breeding populations. This also reveals the evolutionary pressure on the alleles as well as the mutation rate a locus may have experienced over time (Botstein et al., 1980; Shete et al., 2000; Luo et al., 2019). In our study, the average gene diversity was observed to be 0.24, 0.26, and 0.30 for Tamil Nadu, Orissa, and Andhra Pradesh, respectively (Table 2). As a result, the overall gene diversity value was slightly higher than the PIC value, which was expected. The PIC values will always be lower than gene diversity and will become closer to gene diversity as more alleles are added and with rising evenness of allele frequencies (Shete et al., 2000). PIC values are limited to 0.5 due to the biallelic nature of the SNPs (where the two alleles have identical frequencies) (Eltaher et al., 2018) and could possibly be attributable to low mutation rates in SNPs (Coates et al., 2009; Eltaher et al., 2018; Luo et al., 2019). There are some accessions in the total list of accessions which have the same IC numbers followed by X and P; these are accessions which were collected at two different periods from the same areas, hence denoted by X or P. It has been distinctly noticed that these accessions and their original counterparts did not give the same results with SNP markers; this could be due to the high evolutionary drive during collection at different periods. (Kasso and BalaKrishnan, 2013).
The genetic distances were estimated, and the dissimilarity matrix was used to build the NJ tree. The NJ tree of the 2,242 coastal samples showed three major groups. However, nothing very captivating was observed in the clusters formed. Such widely overlapped groups in the NJ tree has also been reported by Xu et al. (2016) on indica rice.
Initially, the population structure of the overall east coast collections revealed four populations. A subsequent population structure analysis gave 19 populations altogether. A weak population structure and low relatedness as revealed in kinship analysis between the east coast rice accessions and the core accessions support the statement by Nachimuthu et al. (2015) that these are critical factors to circumventing spurious data hindering the downstream study (). In the present study, the NJ tree, Bayesian-based STRUCTURE, and AMOVA and PCoA did not show any consensus clustering that could be highlighted. Similar results were also observed by Ambreen et al. in 2018. Also, a population structure study on the rice collection of the east coast states revealed four, three, and five populations for Andhra Pradesh, Orissa, and Tamil Nadu, respectively. In case of Andhra Pradesh, population structure showed a conspicuous grouping of aromatic samples in population 3. Similar type of grouping was reported in basmati rice by Civáň et al. (2019). Admixtures were observed, which suggests that besides pure lines there are samples which are heterogeneous in nature. The mean values of alpha ranged from 0.04 for Orissa, 0.05 for Tamil Nadu, and 0.10 for Arunachal Pradesh (Supplementary Table S4). When the alpha value approaches zero, it means that the majority of individuals are from distinct populations (Li et al., 2014). The values of Fst correspond to a standardized genetic differentiation, suggesting an acceptable population structure. The STRUCTURE analysis indicated good genetic diversity among the rice accessions of the east coast collection. The presence four populations were confirmed by model-based analysis. This method has been used extensively by Edae et al. (2014) to explore association mapping.
The advanced M strategy with minimum redundancy and heuristic approach was used for east coast core collection. The minimum redundancy is required for increasing allelic richness in the core collection; hence, accessions need to be of unique allelic combinations and an unstructured population (Ambreen et al., 2018). A kinship analysis study of the east coast core collection demonstrated a low amount of genetic relatedness, meeting the key condition of an ideal core collection as well as an idealistic association panel (Kumar et al., 2020). Thus, this coastal core set qualifies all the benchmarks of a standard core set.
While validating our results with the northeast rice collection (Roy Choudhury et al., 2014), markers 01-608-4_C and 03-3478-1_C were found to give the highest PIC value of 0.37, while marker 04-19-4_C was least informative giving the lowest PIC of 0.01 in both collections. Marker 01-608-4_C is a locus of evolutionary conserved genes from the saponin family of proteins, which is involved in the sphingolipid metabolic process and active in extracellular space (Bruhn, 2005); in the present study, this marker is highly conserved yet highly polymorphic. However, marker 03-3478-1_C, which is a locus of the GRP (gibberellin-regulated protein) family and an evolutionary conserved gene involved in plants’ defense mechanism as well as in growth (Inomata, 2020), has not been observed polymorphic in both collections, which is contradictory to earlier an report by Mukesh Jain et al. (2014) in rice.
Validation of the same set of SNP markers on two collections (northeast rice and east coast rice collection) has established that the 36-plex SNP assay is sufficient and efficient for initial diversity analysis and core development. Hence, this 36-plex SNP assay can be exploited by researchers for the genetic diversity study and development of core based on their own collections, thus accelerating their breeding program.
Conclusion
This is the first study where India’s east coast rice collections were characterized using SNP markers. The genetic diversity and population structure were studied, and core and mini core collections with maximum diversity and minimum redundancy were developed. A total of 2,242 east coast rice accessions from three different states of India, i.e., Andhra Pradesh, Orissa, and Tamil Nadu, have been characterized, and a wide range of gene diversity and PIC was observed. A phylogenetic analysis of the total east coast rice collection revealed three groups, and a population structure analysis revealed four populations. The 36-SNP assay used in this study was validated by comparing the genetic diversity parameters (gene diversity, PIC, major allele frequency, and heterozygosity) across two different rice collections, i.e., east coastal rice and northeast rice collection, and it was observed the these markers were sufficient to decipher all genetic parameters very efficiently; hence, they can be effectively utilized for core development and diversity study of different rice genotypes.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material; further inquiries can be directed to the corresponding author.
Author Contributions
RS conceived and designed the experiments; DC and RK performed the experiments; RS and DC analyzed the data; VS and RS contributed reagents/materials/analysis tools; DC and RS contributed to the writing of the manuscript; and KS and NS edited the manuscript.
Funding
This work was supported by ICAR grant for the project Network Project on Functional Genomics and Genetic Modification in crops.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The reviewer BPS declared a shared affiliation with one of the authors NKS to the handling editor at the time of the review.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We are thankful to Avantika Maurya and Shantanu Das for their help in kinship analysis. We are grateful and thankful to the Director, NBPGR, New Delhi, who provided facilities for this work.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2021.726152/full#supplementary-material
Supplementary Figure S1 | Neighbor joining tree of the total 2242 east coast rice collection
Supplementary Figure S2 | Neighbor joining tree of 1133 rice collection of Andhra Pradesh
Supplementary Figure S3 | Neighbor joining tree of 378 rice collection of Orissa
Supplementary Figure S4 | Neighbor joining tree of 731 rice collection of Tamil Nadu
Supplementary Figure S5 | Venn diagram showing co-linearity between all three groups of neighbor joining tree and population 1 of population structure of total east coast rice collection
Supplementary Figure S6 | Venn diagram showing co-linearity between all three groups of neighbor joining tree and population 2 of population structure of total east coast rice collection
Supplementary Figure S7 | Venn diagram showing co-linearity between all three groups of neighbor joining tree and population 3 of population structure of total east coast rice collection
Supplementary Figure S8 | Venn diagram showing co- linearity between all three groups of neighbor joining tree and population 4 of population structure of total east coast rice collection
Supplementary Figure S9 | Hierarchical population structure analysis of subpopulation of population 1 total east coast rice collection
Supplementary Figure S10 | Bar plot of population 1 highlighting the sub populations in east coast rice collection
Supplementary Figure S11 | Hierarchical population structure analysis of subpopulation of population 2 total east coast rice collection
Supplementary Figure 12 | Bar plot of population 2 highlighting the sub populations in east coast rice collection
Supplementary Figure S13 | Hierarchical population structure analysis of subpopulation of population 3 total east coast rice collection
Supplementary Figure S14 | Bar plot of population 3 highlighting the sub populations in east coast rice collection
Supplementary Figure S15 | Hierarchical population structure analysis of subpopulation of population 4 total east coast rice collection
Supplementary Figure S16 | Bar plot of population 4 highlighting the sub populations in east coast rice collection
Supplementary Figure S17 | Estimation of population using LnP(D) derived Δk for k from 2 to10 of Andhra Pradesh rice collection
Supplementary Figure S18 | Bar plot of population structure of Andhra Pradesh rice collection
Supplementary Figure S19 | Estimation of population using LnP(D) derived Δk for k from 2 to10 of Orissa rice collection
Supplementary Figure S20 | Bar plot of population structure of Orissa rice collection
Supplementary Figure S21 | Estimation of population using LnP(D) derived Δk for k from 2 to10 of Tamil Nadu rice collection
Supplementary Figure S22 | Bar plot of population structure of Tamil Nadu rice collection
Supplementary Figure S23 | Scattered plot (PCoA) of total east coast rice accessions (2242) (among 4 populations)
Supplementary Figure S24 | Scattered plot (PCoA) of total east coast rice accessions (2242) (among 19 sub populations)
Supplementary Figure S25 | Pie chart showing percentage of molecular variance of rice collection of Andhra Pradesh
Supplementary Figure S26 | Scattered plot (PCoA) of rice collection of Andhra Pradesh
Supplementary Figure S27 | Pie chart showing percentage of molecular variance of rice collection of Orissa
Supplementary Figure S28 | Scattered plot (PCoA) of rice collection of Orissa
Supplementary figure S29 | Pie chart showing percentage of molecular variance of rice collection of Tamil Nadu
Supplementary Figure S30 | Scattered plot (PCoA) of rice collection of Tamil Nadu
Supplementary Figure S31 | Venn Diagram showing the distribution of accessions in the east coast core and the mini-core collection
Supplementary Figure S32 | Pie chart showing percentage of molecular variance of the east coast rice core set (247)
Supplementary Figure S33 | Scattered plot (PCoA) of the east coast rice core set (247)
Supplementary Figure S34 | Histogram showing the kinship status of rice accessions in the east coast core collection
Supplementary Figure S35 | Comparative analysis of heterozygosity and major allele frequency values across eat coast rice and north-east rice collection
References
Agrama, H. A., Yan, W., Lee, F., Fjellstrom, R., Chen, M.-H., Jia, M., et al. (2009). Genetic Assessment of a Mini-Core Subset Developed from the USDA Rice Genebank. Crop Sci. 49, 1336–1346. doi:10.2135/cropsci2008.06.0551
Amanullah, M. M., Natarajan, S., Vanathi, D., Ramasamy, S., and Sathyamoorthi, K. (2007). Lowland Rice in Coastal Saline Soils - A Review. Agric. Rev. 28, 235–238.
Ambreen, H., Kumar, S., Kumar, A., Agarwal, M., Jagannath, A., and Goel, S. (2018). Association Mapping for Important Agronomic Traits in Safflower (Carthamus tinctorius L.) Core Collection Using Microsatellite Markers. Front. Plant Sci. 9, 402. doi:10.3389/fpls.2018.00402
Barrett, B. A., and Kidwell, K. K. (1998). AFLP‐Based Genetic Diversity Assessment Among Wheat Cultivars from the Pacific Northwest. Crop Sci. 38, 1261–1271. doi:10.2135/cropsci1998.0011183X003800050025x
Bisht, I. S., Mahajan, R. K., Loknathan, T. R., and Agrawal, R. C. (1998). Diversity in Indian Sesame Collection and Stratification of Germplasm Accessions in Different Diversity Groups. Genet. Resour. Crop Evol. 45, 325–335. doi:10.1023/A:1008652420477
Botstein, D., White, R. L., Skolnick, M., and Davis, R. W. (1980). Construction of a Genetic Linkage Map in Man Using Restriction Fragment Length Polymorphisms. Am. J. Hum. Genet. 32, 314–331.
Bradbury, P. J., Zhang, Z., Kroon, D. E., Casstevens, T. M., Ramdoss, Y., and Buckler, E. S. (2007). TASSEL: Software for Association Mapping of Complex Traits in Diverse Samples. Bioinformatics 23, 2633–2635. doi:10.1093/bioinformatics/btm308
Brown, A. H. D. (1989). “The Case for Core Collections,” in The Use of Plant Genetic Resources. Editors A. H. D. Brown, O. H. Frankel, D. R. Marshall, and J. T. Williams (England, London: Cambridge University Press), 136–156.
Bruhn, H. (2005). A Short Guided Tour through Functional and Structural Features of Saposin-like Proteins. Biochem. J. 389, 249–257. doi:10.1042/BJ20050051
Chang, T.-T. (1976). The Origin, Evolution, Cultivation, Dissemination, and Diversification of Asian and African Rices. Euphytica 25, 425–441. doi:10.1007/bf00041576
Chen, W., Hou, L., Zhang, Z., Pang, X., and Li, Y. (2017). Genetic Diversity, Population Structure, and Linkage Disequilibrium of a Core Collection of Ziziphus Jujuba Assessed with Genome-wide SNPs Developed by Genotyping-By-Sequencing and SSR Markers. Front. Plant Sci. 8, 575–588. doi:10.3389/fpls.2017.00575
Chen, X., and Sullivan, P. F. (2003). Single Nucleotide Polymorphism Genotyping: Biochemistry, Protocol, Cost and Throughput. Pharmacogenomics J. 3, 77–96. doi:10.1038/sj.tpj.6500167
Cheon, K.-S., Baek, J., Cho, Y.-i., Jeong, Y.-M., Lee, Y.-Y., Oh, J., et al. (2018). Single Nucleotide Polymorphism (SNP) Discovery and Kompetitive Allele-specific PCR (KASP) Marker Development with Korean Japonica Rice Varieties. Plant Breed. Biotech. 6, 391–403. doi:10.9787/pbb.2018.6.4.391
Choudhury, B., Khan, M. L., and Dayanandan, S. (2013). Genetic Structure and Diversity of Indigenous rice (Oryza Sativa) Varieties in the Eastern Himalayan Region of Northeast India. SpringerPlus 2., Springer, 228–237. doi:10.1186/2193-1801-2-228
Civáň, P., Ali, S., Batista-Navarro, R., Drosou, K., Ihejieto, C., Chakraborty, D., et al. (2019). Origin of the Aromatic Group of Cultivated Rice (Oryza Sativa L.) Traced to the Indian Subcontinent. Genome Biol. Evol. 11, 832–843. doi:10.1093/gbe/evz039
Coates, B. S., Sumerford, D. V., Miller, N. J., Kim, K. S., Sappington, T. W., Siegfried, B. D., et al. (2009). Comparative Performance of Single Nucleotide Polymorphism and Microsatellite Markers for Population Genetic Analysis. J. Hered. 100, 556–564. doi:10.1093/jhered/esp028
Das, B., Sengupta, S., Parida, S. K., Roy, B., Ghosh, M., Prasad, M., et al. (2013). Genetic Diversity and Population Structure of rice Landraces from Eastern and North Eastern States of India. BMC Genet. 14, 71–85. doi:10.1186/1471-2156-14-71
Edae, E. A., Byrne, P. F., Haley, S. D., Lopes, M. S., and Reynolds, M. P. (2014). Genome-wide Association Mapping of Yield and Yield Components of spring Wheat under Contrasting Moisture Regimes. Theor. Appl. Genet. 127, 791–807. doi:10.1007/s00122-013-2257-8
El Bakkali, A., Haouane, H., Moukhli, A., Costes, E., van Damme, P., and Khadari, B. (2013). Construction of Core Collections Suitable for Association Mapping to Optimize Use of Mediterranean Olive (Olea Europaea l.) Genetic Resources. PLoS One 8, e61265. doi:10.1371/journal.pone.0061265
Eltaher, S., Sallam, A., Belamkar, V., Emara, H. A., Nower, A. A., Salem, K. F. M., et al. (2018). Genetic Diversity and Population Structure of F3:6 nebraska winter Wheat Genotypes Using Genotyping-By-Sequencing. Front. Genet. 9, 76–84. doi:10.3389/fgene.2018.00076
Evanno, G., Regnaut, S., and Goudet, J. (2005). Detecting the Number of Clusters of Individuals Using the Software STRUCTURE: a Simulation Study. Mol. Ecol. 14, 2611–2620. doi:10.1111/j.1365-294X.2005.02553.x
Frankel, O. H. (1984). “Genetic Perspective of Germplasm Conservation,” in Genetic Manipulation: Impact on Man and Society. Editors A WK. Llimensee, WL Peacock, and P Starlinger (England: Cambridge University Press), 161–170.
Gonzaga, Z. J., Aslam, K., Septiningsih, E. M., and Collard, B. C. Y. (2015). Evaluation of SSR and SNP Markers for Molecular Breeding in Rice. Plant Breed. Biotech. 3, 139–152. doi:10.9787/PBB.2015.3.2.139
Gouesnard, B., Bataillon, T. M., Decoux, G., Razole, C., Schoen, D. J., and david, J. L. (2001). MSTRAT: An Algorithm for Building Germ Plasm Core Collections by Maximizing Allelic or Phenotypic Richness. J. Hered. 92, 93–94. doi:10.1093/jhered/92.1.93
Hu, J., Zhu, J., and Xu, H. M. (2000). Methods of Constructing Core Collections by Stepwise Clustering with Three Sampling Strategies Based on the Genotypic Values of Crops. Theor. Appl. Genet. 101, 264–268. doi:10.1007/s001220051478
Huang, X., Wei, X., Sang, T., Zhao, Q., Feng, Q., Zhao, Y., et al. (2010). Genome-wide Association Studies of 14 Agronomic Traits in rice Landraces. Nat. Genet. 42, 961–967. doi:10.1038/ng.695
Huang, X., Zhao, Y., Wei, X., Li, C., Wang, A., Zhao, Q., et al. (2012). Genome-wide Association Study of Flowering Time and Grain Yield Traits in a Worldwide Collection of rice Germplasm. Nat. Genet. 44, 32–39. doi:10.1038/ng.1018
Inomata, N. (2020). Gibberellin-regulated Protein Allergy: Clinical Features and Cross-Reactivity. Allergol. Int. 69, 11–18. doi:10.1016/j.alit.2019.10.007
Jain, M., Moharana, K. C., Shankar, R., Kumari, R., and Garg, R. (2014). Genomewide Discovery of DNA Polymorphisms in rice Cultivars with Contrasting Drought and Salinity Stress Response and Their Functional Relevance. Plant Biotechnol. J. 12, 253–264. doi:10.1111/pbi.12133
Jasim Aljumaili, S., Rafii, M. Y., Latif, M. A., Sakimin, S. Z., Arolu, I. W., and Miah, G. (2018). Genetic Diversity of Aromatic Rice Germplasm Revealed by SSR Markers. Biomed. Res. Int. 2018, 1–11. doi:10.1155/2018/7658032
Kasso, M., and Balakrishnan, M. (2013). Ex Situ Conservation of Biodiversity with Particular Emphasis to Ethiopia. Hindawi Publishing Corporation, ISRN Biodiversity, 1–11. Article ID 985037, 2013.
Kim, K.-W., Chung, H.-K., Cho, G.-T., Ma, K.-H., Chandrabalan, D., Gwag, J.-G., et al. (2007). PowerCore: a Program Applying the Advanced M Strategy with a Heuristic Search for Establishing Core Sets. Bioinformatics 23, 2155–2162. doi:10.1093/bioinformatics/btm313
Kumar, A., Kumar, S., Singh, K. B. M., Prasad, M., and Thakur, J. K. (2020). Designing a Mini-Core Collection Effectively Representing 3004 Diverse rice Accessions. Plant Commun. 1, 100049. doi:10.1016/j.xplc.2020.100049
Li, F. P., Lee, Y. S., Kwon, S. W., Li, G., and Park, Y. J. (2014). Analysis of Genetic Diversity and Trait Correlations Among Korean Landrace rice (Oryza Sativa L.). Genet. Mol. Res. 13, 6316–6331. doi:10.4238/2014.april.14.12
Liu, K., and Muse, S. V. (2005). PowerMarker: an Integrated Analysis Environment for Genetic Marker Analysis. Bioinformatics 21, 2128–2129. doi:10.1093/bioinformatics/bti282
Luan, S., Chiang, T.-Y., and Gong, X. (2006). High Genetic Diversity vs. Low Genetic Differentiation in Nouelia Insignis (Asteraceae), a Narrowly Distributed and Endemic Species in China, Revealed by ISSR Fingerprinting. Ann. Bot. 98, 583–589. doi:10.1093/aob/mcl129
Luo, Z., Brock, J., Dyer, J. M., Kutchan, T., Schachtman, D., Augustin, M., et al. (2019). Genetic Diversity and Population Structure of a Camelina Sativa Spring Panel. Front. Plant Sci. 10, 184–195. doi:10.3389/fpls.2019.00184
McCouch, S. R., Wright, M. H., Tung, C.-W., Maron, L. G., McNally, K. L., Fitzgerald, M., et al. (2016). Open Access Resources for Genome-wide Association Mapping in rice. Nat. Commun. 7, 1–13. doi:10.1038/ncomms10532
McCouch, S. R., Zhao, K., Wright, M., Tung, C.-W., Ebana, K., Thomson, M., et al. (2010). Development of Genome-wide SNP Assays for rice. Breed. Sci. 60, 524–535. doi:10.1270/jsbbs.60.524
Mourad, A. M. I., Belamkar, V., and Baenziger, P. S. (2020). Molecular Genetic Analysis of spring Wheat Core Collection Using Genetic Diversity, Population Structure, and Linkage Disequilibrium. BMC Genomics 21, 434. doi:10.1186/s12864-020-06835-0
Nachimuthu, V. V., Muthurajan, R., Duraialaguraja, S., Sivakami, R., Pandian, B. A., Ponniah, G., et al. (2015). Analysis of Population Structure and Genetic Diversity in rice Germplasm Using SSR Markers: An Initiative towards Association Mapping of Agronomic Traits in Oryza Sativa. Rice 8, 30–54. doi:10.1186/s12284-015-0062-5
Nei, M., Tajima, F., and Tateno, Y. (1983). Accuracy of Estimated Phylogenetic Trees from Molecular Data. J. Mol. Evol. 19, 153–170. doi:10.1007/bf02300753
Odong, T. L., Jansen, J., van Eeuwijk, F. A., and van Hintum, T. J. L. (2013). Quality of Core Collections for Effective Utilisation of Genetic Resources Review, Discussion and Interpretation. Theor. Appl. Genet. 126, 289–305. doi:10.1007/s00122-012-1971-y
Pathaichindachote, W., Panyawut, N., Sikaewtung, K., Patarapuwadol, S., and Muangprom, A. (2019). Genetic Diversity and Allelic Frequency of Selected Thai and Exotic Rice Germplasm Using SSR Markers. Rice Sci. 26, 393–403. doi:10.1016/j.rsci.2018.11.002
Peakall, R., and Smouse, P. E. (2012). GenAlEx 6.5: Genetic Analysis in Excel. Population Genetic Software for Teaching and Research-Aan Update. Bioinformatics 28, 2537–2539. doi:10.1093/bioinformatics/bts460
Pritchard, J. K., Stephens, M., and Donnelly, P. (2000). Inference of Population Structure Using Multilocus Genotype Data. Genetics 155, 945–959. doi:10.1111/j.1471-8286.2007.01758.x10.1093/genetics/155.2.945
Rafalski, A. (2002). Applications of Single Nucleotide Polymorphisms in Crop Genetics. Curr. Opin. Plant Biol. 5, 94–100. doi:10.1016/S1369-5266(02)00240-6
Rambaut, A. (2010). FigTree v1.3.1. A Graphical Viewer of Phylogenetic Trees. Available at: http://tree.bio.ed.ac.uk/software/figtree/ (Accessed February 25, 2021).
Rashmi, D., Bisen, P., Saha, S., Loitongbam, B., Singh, S., Pallavi, P., et al. (2017). Genetic Diversity Analysis in Rice (Oryza Sativa L.) Accessions Using SSR Markers. Intern. Jour. Agricul., Environ. Biotech. 10, 457–467. doi:10.5958/2230-732X.2017.00057.2
Raybould, A. F., Mogg, R. J., and Clarke, R. T. (1996). The Genetic Structure of Beta Vulgaris Ssp. Maritima (Sea Beet) Populations: RFLPs and Isozymes Show Different Patterns of Gene Flow. Heredity 77, 245–250. doi:10.1038/hdy.1996.138
Reif, J. C., Zhang, P., Dreisigacker, S., Warburton, M. L., van Ginkel, M., Hoisington, D., et al. (2005). Wheat Genetic Diversity Trends during Domestication and Breeding. Theor. Appl. Genet. 110, 859–864. doi:10.1007/s00122-004-1881-8
Roy Choudhury, D., Singh, N., Singh, A. K., Kumar, S., Srinivasan, K., Tyagi, R. K., et al. (2014). Analysis of Genetic Diversity and Population Structure of Rice Germplasm from North-Eastern Region of India and Development of a Core Germplasm Set. PLoS One 9, e113094. doi:10.1371/journal.pone.0113094
Ryan, M. C., Stucky, M., Wakefield, C., Melott, J. M., Akbani, R., Weinstein, J. N., et al. (2020). Interactive Clustered Heat Map Builder: An Easy Web-Based Tool for Creating Sophisticated Clustered Heat Maps. F1000Res 8, 1750. doi:10.12688/f1000research.20590.2
Schlötterer, C. (2004). The Evolution of Molecular Markers - Just a Matter of Fashion. Nat. Rev. Genet. 5, 63–69. doi:10.1038/nrg1249
Seo, J., Lee, G., Jin, Z., Kim, B., ChinKoh, J. H., and Koh, H. J. (2020). Development and Application of Indica–Japonica SNP Assays Using the Fluidigm Platform for rice Genetic Analysis and Molecular Breeding. Mol. Breed. 40, 1–16. doi:10.1007/s11032-020-01123-x
Shete, S., Tiwari, H., and Elston, R. C. (2000). On Estimating the Heterozygosity and Polymorphism Information Content Value. Theor. Popul. Biol. 57, 265–271. doi:10.1006/tpbi.2000.1452
Singh, N., Choudhury, D. R., Singh, A. K., Kumar, S., Srinivasan, K., Tyagi, R. K., et al. (2013). Comparison of SSR and SNP Markers in Estimation of Genetic Diversity and Population Structure of Indian Rice Varieties. PLoS One 8, e84136. doi:10.1371/journal.pone.0084136
Singh, N., Jayaswal, P. K., Panda, K., Mandal, P., Kumar, V., Singh, B., et al. (2015). Single-copy Gene Based 50 K SNP Chip for Genetic Studies and Molecular Breeding in rice. Sci. Rep. 5, 11600. doi:10.1038/srep11600
Singh, N. K., Dalal, V., Batra, K., Singh, B. K., Chitra, G., Singh, A., et al. (2006). Single-copy Genes Define a Conserved Order between rice and Wheat for Understanding Differences Caused by Duplication, Deletion, and Transposition of Genes. Funct. Integr. Genomics. 7, 17–35. doi:10.1007/s10142-006-0033-4
Suvi, W. T., Shimelis, H., Laing, M., Mathew, I., and Shayanowako, A. I. T. (2020). Assessment of the Genetic Diversity and Population Structure of rice Genotypes Using SSR Markers. Acta Agriculturae Scand. Section B - Soil Plant Sci. 70, 76–86. doi:10.1080/09064710.2019.1670859
Syvänen, A. C. (2001). Accessing Genetics Variation: Genotyping Single Nucleotide Polymorphisms. Nat. Rev. Genet. 2, 930–942.
Tarang, A., Kordrostami, M., Shahdi Kumleh, A., Hosseini Chaleshtori, M., Forghani Saravani, A., Ghanbarzadeh, M., et al. (2020). Study of Genetic Diversity in rice (Oryza Sativa L.) Cultivars of Central and Western Asia Using Microsatellite Markers Tightly Linked to Important Quality and Yield Related Traits. Genet. Resour. Crop Evol. 67, 1537–1550. doi:10.1007/s10722-020-00927-2
Thomson, M. J., Zhao, K., WrightMcNally, M. K. L., McNally, K. L., Rey, J., Tung, C.-W., et al. (2012). High-throughput Single Nucleotide Polymorphism Genotyping for Breeding Applications in rice Using the BeadXpress Platform. Mol. Breed. 29, 875–886. doi:10.1007/s11032-011-9663-x
Upadhyaya, H. D., Dwivedi, S. L., Sharma, S., Lalitha, N., Singh, S., Varshney, R. K., et al. (2014). Enhancement of the Use and Impact of Germplasm in Crop Improvement. Plant Genet. Resour. 12, S155–S159. doi:10.1017/S1479262114000458
van Hintum, Th. J. L., Brown, A. H. D., Spillane, C., and Hodgkin, T. (2000). Core Collections of Plant Genetic Resources. IPGRI Technical Bulletin No. 3. Rome, Italy: International Plant Genetic Resources Institute.
Wang, J. C., Hu, J., Xu, H. M., and Zhang, S. (2007). A Strategy on Constructing Core Collections by Least Distance Stepwise Sampling. Theor. Appl. Genet. 115, 1–8. doi:10.1007/s00122-007-0533-1
Xu, Q., Yuan, X., Wang, S., Feng, Y., Yu, H., Wang, Y., et al. (2016). The Genetic Diversity and Structure of Indica rice in China as Detected by Single Nucleotide Polymorphism Analysis. BMC Genet. 17, 53. doi:10.1186/s12863-016-0361-x
Xu, X., Liu, X., Ge, S., Jensen, J. D., Hu, F., Li, X., et al. (2011). Resequencing 50 Accessions of Cultivated and Wild rice Yields Markers for Identifying Agronomically Important Genes. Nat. Biotechnol. 30, 105–111. doi:10.1038/nbt.2050
Yamamoto, T., Nagasaki, H., Yonemaru, J.-i., Ebana, K., Nakajima, M., Shibaya, T., et al. (2010). Fine Definition of the Pedigree Haplotypes of Closely Related rice Cultivars by Means of Genome-wide Discovery of Single-Nucleotide Polymorphisms. BMC Genomics 11, 267. doi:10.1186/1471-2164-11-267
Yan, W., Rutger, J. N., Bryant, R. J., Bockelman, H. E., Fjellstrom, R. G., Chen, M.-H., et al. (2007). Development and Evaluation of a Core Subset of the USDA rice Germplasm Collection. Crop Sci. 47, 869–876. doi:10.2135/cropsci2006.07.0444
Keywords: rice, SNP markers, genetic diversity, genotyping, SNP, coastal rice
Citation: Choudhury DR, Kumar R, S VD, Singh K, Singh NK and Singh R (2021) Identification of a Diverse Core Set Panel of Rice From the East Coast Region of India Using SNP Markers. Front. Genet. 12:726152. doi: 10.3389/fgene.2021.726152
Received: 16 June 2021; Accepted: 26 October 2021;
Published: 25 November 2021.
Edited by:
Ahmed Sallam, Assiut University, EgyptReviewed by:
Bikram Pratap Singh, Indian Council of Agricultural Research, IndiaYasser Shaaban Sayed Moursi, Fayoum University, Egypt
Shamseldeen Shehabeldin Eltaher, University of Sadat City, Egypt
Copyright © 2021 Choudhury, Kumar, S, Singh, Singh and Singh. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Rakesh Singh, rakesh.singh2@icar.gov.in