Kyu-Won Kim
Bhagwat Nawade
Jungrye Nam1
Sang-Ho Chu
Jungmin Ha
Yong-Jin Park- 1Center for Crop Breeding on Omics and Artificial Intelligence, Kongju National University, Yesan, South Korea
- 2Department of Plant Science, Gangneung-Wonju National University, Gangneung, South Korea
- 3Department of Plant Resources, College of Industrial Sciences, Kongju National University, Yesan, South Korea
Rice is a globally cultivated crop and is primarily a staple food source for more than half of the world’s population. Various single-nucleotide polymorphism (SNP) arrays have been developed and utilized as standard genotyping methods for rice breeding research. Considering the importance of SNP arrays with more inclusive genetic information for GWAS and genomic selection, we integrated SNPs from eight different data resources: resequencing data from the Korean World Rice Collection (KRICE) of 475 accessions, 3,000 rice genome project (3 K-RGP) data, 700 K high-density rice array, Affymetrix 44 K SNP array, QTARO, Reactome, and plastid and GMO information. The collected SNPs were filtered and selected based on the breeder’s interest, covering all key traits or research areas to develop an integrated array system representing inclusive genomic polymorphisms. A total of 581,006 high-quality SNPs were synthesized with an average distance of 200 bp between adjacent SNPs, generating a 580 K Axiom Rice Genotyping Chip (580 K _ KNU chip). Further validation of this array on 4,720 genotypes revealed robust and highly efficient genotyping. This has also been demonstrated in genome-wide association studies (GWAS) and genomic selection (GS) of three traits: clum length, heading date, and panicle length. Several SNPs significantly associated with cut-off, −log10 p-value >7.0, were detected in GWAS, and the GS predictabilities for the three traits were more than 0.5, in both rrBLUP and convolutional neural network (CNN) models. The Axiom 580 K Genotyping array will provide a cost-effective genotyping platform and accelerate rice GWAS and GS studies.
1 Introduction
Rice (Oryza sativa) is a staple food source for more than half of the global population (Wing et al., 2018). Rice production is expected to increase by 50–70% by 2050, with improved quality, reliability, and sustainability of global food demand (Zhao et al., 2011; Seck et al., 2012). However, sustainable production with fewer resources will require the efficient utilization of high-throughput and intensive systems in increasingly variable environments (Rasheed et al., 2017; Yu et al., 2022).
With advances in high-throughput sequencing technologies, -omics-based studies on rice have progressed considerably, enabling the efficient identification of a large number of single nucleotide polymorphisms (SNPs) (Wing et al., 2018; Nguyen et al., 2019). In addition to being highly prevalent, biallelic, codominant, and stable SNPs also play a significant role in phenotypic variation. SNPs are the most effective and highly informative genetic markers used to unravel functional variants underlying traits for crop improvement (Yu et al., 2022). Next-generation sequencing technologies enable accurate detection of SNPs from various genomic backgrounds. With the availability of several million SNPs, the challenge is efficient and economical genotyping of these SNPs (Rasheed et al., 2017).
High-throughput SNP genotyping is an attractive genotyping tool for identifying sequence polymorphisms (Rasheed et al., 2017; Guo et al., 2021). It is typically accomplished using SNP arrays or ‘chips’ or genotyping-by-sequencing (GBS). SNP arrays and GBS are cost-effective for genotyping thousands to millions of SNPs, whereas PCR-based genotyping requires hundreds to a few thousand SNPs, and is laborious, time-consuming, and suitable for small-scale genotyping. Although GBS has low ‘set-up’ and per-sample costs and performs SNP discovery and genotyping simultaneously, its experimental operation and data analysis are beyond the reach of average breeders (Verlouw et al., 2021). In contrast, high-throughput genotyping arrays can be used repeatedly to genotype different populations in a short period of time with straightforward data analysis (Rasheed et al., 2017).
Several genotyping platforms, including Illumina BeadXpress (Chen et al., 2011; Thomson et al., 2012), Fluidigm platform (Seo et al., 2020), Illumina Infinium (Yu et al., 2014; Thomson et al., 2017; Morales et al., 2020), and Affymetrix (Zhao et al., 2011; Singh et al., 2015; McCouch et al., 2016) have been developed and utilized in rice molecular breeding. RiceSNP50 was designed based on over 10M SNP loci from the resequencing data of 801 rice varieties (Chen et al., 2014). OsSNPnks include 50 K high-quality non-redundant SNPs (Singh et al., 2015). Because they mainly consist of SNPs within single-copy genes, SNP information has been widely applied in evolutionary and domestication-related studies of the Oryza genus. McCouch et al. constructed a high-density rice array consisting of 700 K SNPs surpassing the largest publicly available genotyping platform for any crop species (McCouch et al., 2016). Cornell_ 6 K _Array_Infinium_Rice (C6AIR) was designed and developed to be polymorphic within and between target germplasm groups and to map populations of interest (Thomson et al., 2017). Therefore, C6AIR provides a highly informative dataset to Cornell University and the IRRI, indicating the importance of data resources in designing SNP arrays. C6AIR was updated to C7AIR, covering polymorphisms between and within O. sativa, O. glaberrima, O. rufipogon, and O. nivara (Morales et al., 2020). Seo et al. developed two 96-plex indica-japonica SNP genotyping assays for particular target populations containing functional SNPs associated with agronomic traits for efficient genotyping (Seo et al., 2020). High-throughput genotyping platforms play critical roles in genetic diversity, gene mapping, germplasm resource analysis, genome-wide association study (GWAS), evolution analysis, and genomic selection (Xu et al., 2021).
However, most of these arrays included whole-genome random SNPs but were not inclusive of SNPs related to key traits or research interests. Taking advantage of the accrued rice genomic sequence data, we collected SNPs from eight different highly informative datasets and selected high-throughput SNPs across the breeder’s research interests (Figure 1) to develop a large-scale genotyping array on the Affymetrix platform. Further validation of this array using a large set of accessions for GWAS analysis and genomic selection for different traits has demonstrated its usefulness in the global rice community.
Figure 1 Schematic of 580K _ KNU chip design. SNPs were collected from eight different datasets, including 1) Resequencing data of 475 KRICE, 2) 3K-RGP from IRRI, 3) 700K rice array from Cornell University, 4) Affymetrix 44K Rice Chip, 5) QTARO database, 6) Plant Reactome Gramene database, 7) Chloroplast and mitochondria genomic sequence and 8) GMO information.
2 Materials and methods
2.1 Sequence resources
We utilized different databases for designing the Axiom rice genotyping chip:1) Resequencing data of Korean World Rice Collection (KRICE) of 475 accessions composed of 417 cultivated and 58 wild accessions (Phitaktansakul et al., 2021), 2) Rice genome project data of 3,000 accessions (3 K-RGP) from the International Rice Research Institute (IRRI) (https://snp-seek.irri.org), 3) High-Density Rice Array (HDRA, 700K) array reported from Cornell university (McCouch et al., 2016), 4) Affymetrix 44 K Rice Chip (Affy44K) (Zhao et al., 2011), 5) QTARO database (http://qtaro.abr.affrc.go.jp/) (Yonemaru et al., 2010), 6) Plant Reactome Gramene Pathways database (https://plantreactome.gramene.org/) 7) Chloroplast (Tong et al., 2016; Cheng et al., 2019) and mitochondria (Tong et al., 2017) genomic sequences of japonica and rufipogon and 8) GMO, we used transgenic plants and genomes of various microorganisms as references for GMO marker design. The known GMO events consisted of host and insert regions; we utilized NCBI Primer-BLAST tool to get the full-length products and modified them into Axiom probes for 155 events. Besides, GMO markers were also developed from ds, tDNA, and tos17 as these insertion elements have been employed for generating large-scale mutant pools in different crops. We performed BLAST queries against rice reference genome sequences and unmapped sequences were selected for GMO marker design. Further, selected SNPs were aligned with different reference genomes including, japonica (ftp://ftp.ensemblgenomes.org/pub/plants/release-36/fasta/oryza_sativa/dna) (Kawahara et al., 2013), indica (http://rice.hzau.edu.cn/rice/download_ext/MH63RS2.LNNK00000000.fasta.gz) (Zhang et al., 2016), and O. rufipogon (https://www.ncbi.nlm.nih.gov/nuccore/NC_013816.1?report=fasta) (Fujii et al., 2010).
2.2 SNP filtering and integration
A VCF file was created from 3,475 rice accessions (3 K-RGP and KRICE), and SNP/indel sites with MAF<0.05, and missing rate> 0.1 were removed. Obtained SNPs were further enriched with Affymetrix (Affy44K), High-Density Rice Array (HDRA, 700K), and SNPs of selected genes from QTARO and Reactome databases, which resulted in a total of 7,682,442 markers. These markers were categorized into different classes viz., japonica, indica, and rufipogon specific, based on their sources (Table S1). The selected SNPs (called tag-SNPs) and the corresponding flanking sequences were submitted to Affymetrix (Axiom® BioFx Services) service for initial probe screening. The priority was given as 0, 1, and 2, with 2 being the highest priority and 0 being the lowest. We assigned priority 2 to agronomically important genes, that have known pathways, and exist exclusively in either japonica or indica. After removing the tag-SNPs with a design score (pconvert)<0.6, a total of 3,204,347 SNPs met the Affymetrix probe designing criteria (Table S1).
2.3 Selection of SNPs for array development
Among technical suitable variants, we allocated SNPs based on genetic diversity and breeder’s interest, which were divided into five divisions (Figure 2A). In the genetic diversity division, we selected 40 K high-quality SNPs from KRICE (Minimum allele frequency, MAF< 0.01, non-missing) (Phitaktansakul et al., 2021), 220,135 SNPs from the Affymetrix and Cornell chip, and SNPs specific to the indica/japonica group from the 3 K-RGP and KRICE data. We also selected chloroplasts (Tong et al., 2016; Cheng et al., 2019) and mitochondrial genomes (Tong et al., 2017) and filtered high-quality SNPs/indels (Figure 2A).
Figure 2 Development of integrated array system. (A) SNPs collected from genomic databases were prioritized based on breeders’ interest and key research divisions. (B) Summary of distribution of SNPs selected for Axiom 580K Array.
Furthermore, in the key research division, SNPs detected from 11 different genome-wide association studies (GWAS) for candidate genes were filtered with MAF< 0.01 (Figure 2). The evaluation of these studies resulted in 400 K polymorphisms, including 300 K SNPs and 100 K indels, which are associated with different traits.
In the breeder division, the key traits and genes of rice breeder interest were chosen, including starch synthesis, blast resistance, and yield-related genes (Figure 2A). A total of 1800 genes/QTL regions from the QTARO database were used for the SNP collection. In addition, for metabolite division, a set of genes was selected from the plant reactome (https://plantreactome.gramene.org) to retrieve SNPs involved in different pathways or processes in rice (Figure 2A, Table S2).
To prioritize probe sets for polymorphisms and select the final set of SNPs for array design, we used following criteria: (1) 71mer sequence (35 bp flanking sequence of target SNP) was to score the marker; (2) A marker was “not recommended” having one or more polymorphisms within 24 bases; (3) If a marker has the same recommendation for each strand, we recommended tilling the one with the highest pconvert value; (4) A marker/strand was recommended if: pconvert > 0.6, there are no wobbles, and poly count = 0; (5) A marker/strand was not_recommended if: duplicate count > 0, or poly count > 0, or pconvert< 0.4, or wobble distance< 21, or wobble count >= 3; (4) A marker was considered not_possible on a given strand if we cannot build a probe to interrogate the SNP in that direction.
A list of 620,852 candidate SNP loci was sent to Affymetrix Bioinformatics Services (Santa Clara, CA, USA) for array design. The quality of each SNP was assessed again and designated as ‘recommended’, ‘neutral’, ‘not recommended’, and ‘not_possible’ using in silico validation with proprietary software. We retained one SNP marker every 200 bp to ensure a uniform distribution and high density of SNPs throughout the rice genome. The final 580 K Axiom Rice Genotyping Chip (580 K _ KNU chip) contained a total of 581,006 SNP markers (Figure 2B).
2.4 Plant material and phenotyping
A total of 4,720 genotypes were genotyped for initial validation of the 580 K _KNU chip. This 4,720 genotype set was composed of different genotypes, including RILs (1,821), backcross inbred lines (209), backcross lines (BC1F3, 96), F1-F7 (209), breed (1,123), landraces (252), weedy (488), Rufipogon (96), wild type (245), and mutant lines (181) (Table S3).
The traits clum or culm length (cm), heading date, and panicle length (cm) were evaluated for GWAS, and amylose content (%), panicle length (cm), number of grains/panicle, heading date, number of panicles/hill, and 100-grain weight were used in genomic selection analysis. Phenotyping was conducted at the Kongju National University, Yesan, South Korea, during the dry seasons of 2017 and 2018. Five plants were randomly selected from the middle row of the plot, and each parameter was recorded using the Standard Evaluation System (SES) (IRRI, 2002). Days to heading were defined as the time when half of the plants in each accession showed panicles. Amylose content was determined using the iodine colorimetric method at 620 nm absorbance on a UV-1800 spectrophotometer (Shimadzu Co., Kyoto, Japan). Frequency distributions of phenotypic data were tested for normality using the Shapiro–Wilk function in R environment (Royston, 1995).
2.5 Genotyping
Genomic DNA was extracted from young green leaf tissue using a Qiagen plant DNeasy kit (Qiagen, Germantown, MD, USA) and quantified using a NanoDrop spectrophotometer (Thermo Scientific, USA). The DNA quality was checked using a 1% agarose gel. Genomic DNA (200 ng) from all lines was hybridized into arrays using the Affymetrix GeneTitan system, according to the manufacturer’s instructions. SNP genotyping, quality control (QC), and SNP filtering were performed according to the Axiom Genotyping Solution Data Analysis User Guide (http://www.affymetrix.com/). Briefly, genotype calling and QC metrics were performed using Affymetrix Genotyping Console™ (GTC) v.4.2. Samples with a development quality check (DQC) value<0.83 and call rate<0.97 were excluded from further analysis. GTC results were post-processed using the SNPolisher R package (v.3.0). The Ps_Metrics function was used to generate SNP metrics, and the Ps_Classification function with the default setting classified SNPs into six categories: PolyHighResolution (SNPs had good cluster resolution and at least two examples of the minor allele), MonoHighResolution (SNPs had good SNP clustering but less than two samples had the minor allele), Off-Target Variant, CallRateBelow-Threshold (SNPs had call rates CR below the threshold, but the other properties were above the threshold), NoMinorHom (SNPs had good cluster resolution but no samples had the minor allele), and Other (more than one cluster property was below the threshold) (Gao et al., 2014b). Furthermore, SNP QC metrics, including call rate (CR, ≥97%), Fisher’s linear discriminant (FLD, ≥3.6), heterozygous strength offset (HetSO, ≥_0.1), and homozygote ratio offset (HomRO, ≥0.3) values, were applied to assess SNPs. The remaining SNPs retained for further analysis were annotated using an in-house Python script. SNP distribution and count across the rice genome were analyzed using a 100 kb sliding window with an R package.
2.6 Genome-wide association studies
High-quality SNPs obtained from the Affymetrix chip were used in a genome-wide analysis of associations. A GWAS was performed on three phenotypic traits: clum length, heading date, and panicle length. Association analyses were conducted using the Genome Association and Prediction Integrated Tool (GAPIT) (Lipka et al., 2012) and TASSEL 5.0 (Trait Analysis by Association, Evolution, and Linkage) (Bradbury et al., 2007). GWAS analysis was performed using the mixed linear model (MLM) of GAPIT (Yu et al., 2006) to predict the association between each SNP and the phenotypic data. The kinship (K) matrix represents the variance-covariance matrix between individuals. The R package qqman (https://cran.r-project.org/web/packages/qqman/index.html) was used to draw Manhattan plots. A P value of 3.16 × 10−7 was used to consider marker-trait association (MTA) as significant.
2.7 GS Models
2.7.1 Penalized Regression Model
The ridge regression best linear unbiased prediction (rr-BLUP) model was implemented using the R package rrBLUP (Endelman, 2011). The model is described as follows.
where y is an N × 1 vector of adjusted means for all genotypes, μ is the overall mean, Z is an N × M matrix of markers, u is a vector of marker effects as u ~ N(0, Iσ2u), and e is the residual error with e ~ N(0, Iσ2e).
GS was performed with fourfold cross-validation by including 80% of the samples in the training population and predicting the genomic estimated breeding values (GEBVs) of the remaining 20% of the samples. For the accuracy assessment, two 50 replication sets were performed, with each replicate consisting of five iterations.
2.7.2 Convolutional Neural Networks
Convolutional neural network (CNN) is a deep learning model that accommodates inputs distributed along with space patterns (Pérez-Enciso and Zingaretti, 2019). In a CNN, the input data first passes through a convolutional layer, followed by a pooling layer, dropout layer, fully connected dense layer, batch normalization layer, and finally to the output layer containing one node with the predicted trait value. During each convolutional layer, the CNN applies kernels and filters, and performs the convolution operation with a predefined width and strides, providing the same weights for all SNP marker windows. The filter moves for the same window size across the input SNP markers, and the CNN obtains a locally weighted sum (Sandhu et al., 2021). The earlyStopping function in Keras (https://keras.io/callbacks/#earlystopping) was applied to control model overfitting (Zingaretti et al., 2020). A pooling layer is added after each convolutional layer for dimensionality reduction, and the filters are invariant to small changes in the input. Finally, pooling results in a smoothed representation and merging of the kernel output of the previous convolutional layer by taking the minimum, mean, and maximum (Bellot et al., 2018).
3 Result
3.1 Alignment, SNP selection, and Axiom Array design
A detailed description of the detection, filtering, and final selection of SNPs included in the array is provided in the Methods section and Figure 1. SNPs collected from eight different datasets were selected based on the breeder’s interests and key research divisions (Figure 2).
Alignment of whole-genome resequencing of 1) 475 Korean World Rice Collection (KRICE) accessions and 2) 3,000 rice genome project (3 K-RGP) data accessions against both indica and japonica rice reference genomes were performed to identify sequence variations (SNPs and indels) (Kawahara et al., 2013; Zhang et al., 2016). The alignment resulted in the identification of over 3.1 million SNPs in KRICE, of which 40 K high-quality SNPs (MAF< 0.01 and zero missing rates), all exonic SNPs/indels and sub-species specific SNPs were selected. Furthermore, SNPs from the high-density rice array assay from 3) Cornell University, 4) Affymetrix 44 K Rice Chip, 5) QTARO database, and 6) Plant Reactome were mapped and aligned to select potential SNPs. SNPs from (7) chloroplasts (Tong et al., 2016; Cheng et al., 2019) and mitochondrial genomes (Tong et al., 2017) of japonica and rufipogon rice were filtered and 3,449 and 3,329 SNPs/indels were selected, respectively. In GMO markers, 20,895 markers were derived from binary vectors, and 2,089, 13,697, and 746 markers derived from ds, tDNA, and tos17, respectively. As described in the ‘Methods’ section, after applying different criteria to identified SNPs data, a set of 620,852 high-quality SNPs were selected. Finally, 581,006 SNPs were tiled on the 580 K _ KNU chip SNP array that includes 3,329 from rufipogon, 34,276 from indica, and 542,432 from japonica (3,449 plastid and 538,983 nuclear). Among japonica specific 538,983 SNPs/indels (500,725 SNPs and 38,258 indels) 126,076 SNPs were for agronomic traits, 162,315 for intergenic, and 94,393 for genic SNPs (Figure 2B). Among the 34,276 indica-specific SNPs/indels, 22,820 were SNPs, and the remaining 11,456 were indels. In the case of 8) GMO, 969 regions were selected from transgenic genes and vectors, and a total of 40,715 probes were selected as candidate SNPs (Figure 2B).
Regarding the distribution of 500,725 japonica-specific SNPs in different parts of the genes, 82,666 (16.51%) SNPs were present in exons, 78,963 (15.77%) in introns, 27,139 (5.42%) in the UTR, and 26,3539 (52.63%) in the intergenic region (Figure 3A). Most indels were detected within an intron of japonica (Figure 3B). In the case of 22,820 indica-specific SNPs, no intergenic SNP were detected, whereas only two indica-specific indels were observed in the intergenic region (Figure 3C, D). Of the indica-specific SNPs, 53.6% (12,234), 27.7% (6,317) and 06.0% (1,367) were distributed within introns, exons, and UTR, respectively (Figure 3C). SNPs were located along each of the 12 rice chromosomes, with an average density of 154 SNPs/100 K and a median density of 130 SNPs/100 K (Figure 4). The average gap between two adjacent SNPs was 200 bp, and gaps between more than 90% were less than 2 kb (Table S4).
Figure 3 Genomic position of selected SNPs and indels. (A) Distribution of japonica-specific SNPs in different genomic regions. (B) Distribution of japonica-specific indels in different genomic regions. (C) Distribution of genomic regions indica-specific SNPs in different genomic regions (D) Distribution of indica-specific indels in different genomic regions.
Figure 4 Distribution of the converted SNPs on the array in 100 Kb windows along the rice chromosomes. (A) Density of japonica-specific SNPs across chromosomes. (B) Density of japonica-specific indels across chromosomes (C) Density of indica-specific SNPs across chromosomes (D) Density of indica-specific indels across chromosomes.
3.2 Genotyping performance of 580K _ KNU chip
The performance of the 580 K _ KNU Axiom Array was evaluated by genotyping eight different sets of genotypes on an integrated Affymetrix GeneTitan® platform. All samples passed the quality assessment with a high DQC value (>0.89) and call rate (>95%), and duplicate samples showed 99% SNP reproducibility. Thus, the genotyping results of 4,720 genotypes validated the chip performance with both a high sample success rate and genotyping call. The SNP genotyping results from the 4,720 genotypes were classified into six categories based on the Affymetrix quality control metrics (Table 1; Figure 5). Based on the filtering parameters, Fisher’s linear discriminant (FLD), HetSO, HomRO, and CR ≥95% as the filtering options, approximately 79.2% (n = 18,087) of the total indica-specific array SNPs were converted. Of the 38,258 designed japonica-specific indels, 23,759 passed the bead representation and decoding quality metrics. In japonica-specific SNPs, 53.13% SNPs were categorized in the ‘PolyHighResolution’ category, whereas 15.94% SNPs were found in ‘OTV’ and 6.08% in the ‘MonoHighResolution’ category (Figure 5A). A total of 7,735 japonica-specific indels (20.2%) were classified ‘PolyHighResolution’ and 23.0% as ‘OTV’ (Figure 5B). While in the case of indica-specific SNPs, 24.8% of SNPs were found in the ‘PolyHighResolution’ category, but the highest SNPs (42.6%) were grouped fall in the ‘other’ category (Figure 5C). A similar trend was found for indica-specific indels with the highest 50.8% in ‘other’ category, followed by 22.6% of SNPs in the ‘PolyHighResolution’ (Figure 5D).
Table 1 Classification of SNPs in the 580K SNP Chip after genotyping on 4720 accessions.
Figure 5 Summary of SNP genotyping data in 4720 accessions using the 580K _ KNU chip. (A) Classification of japonica-specific SNPs (B) Classification of japonica-specific indels (C) Classification of indica-specific SNPs (D) Classification of indica-specific indels. NoMinorHom: No minor homozygote, MonoHighResolution: Monomorphic high-resolution, PolyHighResolution: Polymorphic high-resolution, CallRateBelowThreshold: Call rate below threshold, OTV: Off-target variant.
3.3 Genome-wide association studies and genomic selection
To evaluate the practicality of this array in GWAS, we used a set of 1,288 lines, representing a subset of the 4,720 collection that were phenotyped for clum length, heading date, and panicle length. The frequency distribution curves revealed continuous variation and displayed a normal distribution for all the traits (Figure 6). A total of 60 SNPs were detected for clum length with a threshold above 7 (−log10 p-value) having the highest −log10 p-value of 22 for a marker on chromosome 7 (chr07_ 29494004) (Figure 7A; Table S5). In contrast, 277 markers on chromosomes 6 and 7 were significantly associated with heading date with a 7< −log10 p-value (Figure 7B; Table S6). We obtained a total of 84 significant SNPs exceeding the threshold –log10 p-values of 7 for panicle length (Figure 7C; Table S7). From the GWAS results, we observed several significantly associated SNPs with three traits located on chromosomes 6, 7, and 9 based on the genome-wide significance cut-off, −log10 p-values 7 (Figure 7). For panicle length, three SNPs from Os09g0456100 (OsLP1; LONG PANICLE 1) showed significant association with −log10 p-value of 7<. In addition, significant SNPs mapped on the 2.4 Mbp region of chromosome 9 were belongs to 15 genes that include DENSE AND ERECT PANICLE1 (DEP1; Os09g0441900), lecithine cholesterol acyltransferase (Os09g0444200), Serine carboxypeptidase 42 (Os09g0462875), gibberellin receptor (Os09g0455900), alpha-amylase isozyme 3A precursor (Os09g0457400), etc were significantly associated with panicle length (Table 2). Similarly, annotation of significant SNP regions with -log10 p-values > 7 predicted a total of 42 genes on chromosomes 6 and 7 for heading date and 12 genes for clum length (Tables S8).
Figure 6 Frequency distribution for all the traits used for the genome-wide association study. (A) Clum length, (B) Heading date, and (C) Panicle length. The histograms with purple colors are normal as expected while other colors denote difference from normal distribution. The observed and expected normal distributions fitting for data were represented with dashed and solid red lines, respectively.
Figure 7 Manhattan plots of genome-wide association studies using the 580K _ KNU chip (A) Manhattan plot for Clum length, (B) Manhattan plot for Heading date, and (C) Manhattan plot for Panicle length.
Table 2 List of genes identified from significantly associated SNPs for phenotypic traits.
To demonstrate the performance of the 580 K _KNU chip rice array in quantitative phenotype prediction, we utilized the rrBLUP and CNN statistical models. The predictability of genomic selection was evaluated using four-fold cross-validation, where the sample was randomly partitioned into four parts to estimate the parameters. Finally, all parts were predicted once and used ten times to estimate the parameters. A total of 80, 7 genotypes, including having bred and weedy rice and corresponding 121,208 markers from chip data, were used to evaluate the accuracy of prediction for amylose content (Table S9). Both models showed the highest predictability for amylose content, and values of predictability for all traits with rrBLUP were higher than those of CNN (Figure 8). The predictabilities of the three traits were more than 0.5 in both models. The CNN model was unable to predict the grain number per panicle; however, rrBLUP had a lower value (Figure 8).
Figure 8 Genomic selection predictive ability (r) for different traits by using rrBLUP (A), and CNN (B) model.
4 Discussion
Rice was the first crop species to be fully sequenced (International Rice Genome Sequencing Project, 2005). Thousands of varieties have been re-sequenced (The 3,000 rice genomes project, 2014; Duitama et al., 2015), and de novo assemblies have been performed for several subspecies (Yu et al., 2002; Kawahara et al., 2013; Schatz et al., 2014). Based on accrued genomic information, QTL studies and GWAS have been good strategies for understanding the genetic basis underlying complex traits in rice (Yano et al., 2019; Wang et al., 2020; Yuan et al., 2020), and major associations have been identified as important traits for the last few decades. In current rice breeding programs, a major limitation in the genetic dissection of agronomically important traits is the tight linkage between undesirable loci and preferable loci, for example, two loci within a linkage disequilibrium (LD) block (Xiao et al., 2021). To break tightly linked loci, breeding experts either have to enlarge the population size or advance generations until the block is dissected and the marker interval is minimized. High-throughput genotyping arrays can genotype hundreds of thousands of markers over a large number of samples in a short timeframe. Owing to this technology, the breeding system enables the handling of the population required to dissect LD blocks in a relatively short time. Through association studies and linkage mapping, the rate of development of trait-linked DNA markers can be accelerated, and the breeding cycle can be dramatically reduced, even for tightly linked traits within LD blocks. In the current rice breeding program, where major associations have been detected for important traits, SNP array chips must be designed based on breeders’ interests to address challenging problems in rice breeding, along with efficient and fast genotyping.
In this study, the 580 K _KNU chip array was developed based on inclusive genomic polymorphisms targeting breeders’ interests, covering all key traits and research areas. To construct an integrated array system, in addition to the SNPs from 475 KRICE, SNPs from major rice chips that had been published previously were added after strict filtration (Figure 1). To make the 580 K _KNU chip more informative for the current rice breeding program, we selected candidate SNPs associated with important traits to reflect breeders’ interests (Figure 2). Hence, along with a dense marker interval of 200 bp, the 580 K _KNU chip consisting of 581,006 high-quality SNPs representing the genomic polymorphisms in rice is an excellent platform for genetic dissection of agronomic traits and is highly informative, especially for current research topics in rice breeding programs. The development and cultivation of GM crops have expanded over time, and microarray technology is a flexible method of detecting GMO varieties. In this study, we incorporated GMO-specific loci covering 40,715 markers that will be helpful in rice breeding for preliminary screening of GMOs contaminations. Previously such microarray-based GMO detection has been applied at a small scale in detection of GM soya, rice and maize lines (Kim et al., 2010; Turkec et al., 2016; Chen et al., 2021; Kutateladze et al., 2021).
Using the 580 K _KNU chip, we evaluated the associations between the three phenotypes and the 580 K SNP genotypes in different lines. For club length and heading date, the most significant signals were located in the candidate gene Os07g0695100 (Pseudo-Response Regulator37; OsPRR37), which was identified as being responsible for the early heading7-2 (EH7-2)/heading date2 (Hd2) QTL (Koo et al., 2013; Yan et al., 2013; Gao et al., 2014a). In the case of panicle length, the most significant signal (chr09_16415391 with −log10 p-value 21) was observed in the OsDEP1 gene (Os09g0441900) (Table 2), which controls the erect panicle (EP) architecture, which is a typical characteristic of super rice utilized in rice breeding for nearly a century owing to its high yield, lodging tolerance with strong stems, reasonable population structure, and high nitrogen use efficiency (Huang et al., 2009; Zhou et al., 2009; Xu et al., 2016; Sibo et al., 2021). The same marker chr09_16415391 (Os09g0441900) was found to be the second most significant signal (−log10 p-value 22) for clum length (Table 2). Os09g0456100 (OsLP1; LONG PANICLE 1) that encodes Remorin_C-containing proteins showed significant association with panicle length. Previously, two SNPs from the third and fifth exons of LP1 were reported to reduce panicle length (Liu et al., 2016). The clusters of significant SNPs associated with heading date were detected around Os07g0695100 covering an extensive region of 0.15 Mb, followed by a region (~1-0.89 Mb) on chromosome 6 containing Os06g0285900 (embryogenesis transmembrane protein), Os06g0285100 (OsFLA23; fasciclin domain-containing protein) Os06g0286310 (Oxidoreductase-like protein), Os06g0286400 (nitrate induced protein), Os06g0289200 (UDP-glucuronosyl/UDP-glucosyltransferase family protein), Os06g0286351 (Armadillo-type fold domain containing protein), Os06g0285400 (Serine/threonine-specific kinase-like protein), Os06g0296700 (C-type lectin conserved site domain-containing protein) which have not previously been shown to be associated with any trait in rice (Table 2). Similarly, for heading date, 34 and 13 QTLs were detected near Hd1 and Hd17, RFT1, and Hd3a genes, respectively, on chromosome 6, while 10 QTLs were detected near the OsPRR37 gene on chromosome 7 (Ebana et al., 2011; Hori et al., 2015). Therefore, identified candidate loci could be targeted for fine mapping in order to determine the exact genes/alleles underlying these GWAS signals pertaining to agronomic traits.
GBLUP is the most robust method and the most commonly used tool in rice because it provides high predictability (Xu et al., 2014; Spindel et al., 2015; Xu et al., 2021). In GS selection using the 580K_KNU chip, rrBLUP performed better than CNN, but both models showed > 0.5 prediction values (Figure 8). GS predictability is influenced by various factors, including heritability, relatedness between populations, sample size, marker density, genetic architecture, statistical model, and factors. GS accuracy in rice breeding populations has been performed for various quantitative traits varied by trait, population, and model and moderate to high predictive ability has been reported (Xu et al., 2021). Onogi et al., 2015 recorded the predictive ability of heading date (0.8), clum length (0.75), panicle length (0.6), panicle number (0.4), and grain length (0.4) in a population of 110 Asian rice cultivars using GBLUP (Onogi et al., 2015). The genotyping of 413 rice inbred lines with a 44 K chip showed the predictive abilities for florets per panicle (~0.6), flowering time (~0.63), plant height (~0.7), and protein content ~0.44 (Isidro et al., 2015), while in a panel of 363 elite breeding lines predictive abilities of 0.31, 0.34, and 0.63 for grain yield, plant height and flowering time, respectively were reported (Spindel et al., 2015). The predictive abilities for 1000 grain weight were 0.82–0.83 in 210 recombinant inbred lines and 278 hybrids (Xu et al., 2014), and 0.54 in 1495 hybrids derived from incomplete NC II design (Cui et al., 2020). The GS predictive values reported in this study were in accordance with previous reports on rice (Xu et al., 2014; Grenier et al., 2015; Spindel et al., 2015; Xu et al., 2018; Cui et al., 2020). Compared with other arrays, where several major genes have not been fully integrated, our 580 K _KNU chip system might be better suited for GS. Similarly, the gene–coding sequence–haplotype (gcHap)-based GS showed higher predictive ability than SNP-based GS because the gcHap dataset represents the diversity of 45 963 rice genes in 3010 rice accessions (Zhang et al., 2021). GS-specific SNP arrays could further improve rice breeding accuracy, intensity, and efficiency as well as reduce cost and time.
The current manuscript mainly focused on the development of an inclusive SNP array system that will help rice breeders for multiple breeding programs. Also, we have conducted validation analyses for the GWAS and genomic selection. In collaboration with DNA Link, Inc. (Korea), we provide microarray analysis services using our 580 K _KNU chip. Scientists and breeders in Korea are using it for a variety of applications, and we expect to see the results of studies on the use of our SNP array in the near future. As an accurate high-density genotyping tool, the 580 K _KNU chip is an excellent platform for GWAS, QTL mapping, evolutionary studies, genetic diversity, and genomic selection, especially genetic dissection of important traits of breeders’ interest that have not been fully identified. Hence, it will play a pivotal role in rice breeding applications.
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/s.
Author contributions
K-WK, BN and JH wrote the manuscript. K-WK, JN, BN and S-HC collected and analyzed the data. JH and Y-JP supervised and revised the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This study was supported by a National Research Foundation of Korea (NRF) grant funded by the Korean government (MSIT) (No. 2022R1A4A1030348). This work was supported by the Korea Institute of Planning and Evaluation for Technology in Food, Agriculture and Forestry (IPET) through the Digital Breeding Transformation Technology Development Program funded by the Ministry of Agriculture, Food and Rural Affairs (MAFRA) (322060031HD020). This work was carried out with the support of “Cooperative Research Program for Agriculture Science and Technology Development (Project No. PJ015935), Rural Development Administration, Republic of Korea. This work was supported by a research grant from the Kongju National University in 2022.
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.
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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2022.1036177/full#supplementary-material
References
Bellot, P., de los Campos, G., Pérez-Enciso, M. (2018). Can deep learning improve genomic prediction of complex human traits? Genetics 210, 809–819. doi: 10.1534/genetics.118.301298
Bradbury, P. J., Zhang, Z., Kroon, D. E., Casstevens, T. M., Ramdoss, Y., Buckler, E. S. (2007). TASSEL: software for association mapping of complex traits in diverse samples. Bioinformatics 23, 2633–2635. doi: 10.1093/bioinformatics/btm308
Cheng, L., Nam, J., Chu, S.-H., Rungnapa, P., Min, M., Cao, Y., et al. (2019). Signatures of differential selection in chloroplast genome between japonica and indica. Rice 12, 65. doi: 10.1186/s12284-019-0322-x
Chen, H., He, H., Zou, Y., Chen, W., Yu, R., Liu, X., et al. (2011). Development and application of a set of breeder-friendly SNP markers for genetic analyses and molecular breeding of rice (Oryza sativa l.). Theor. Appl. Genet. 123, 869. doi: 10.1007/s00122-011-1633-5
Chen, H., Xie, W., He, H., Yu, H., Chen, W., Li, J., et al. (2014). A high-density SNP genotyping array for rice biology and molecular breeding. Mol. Plant 7, 541–553. doi: 10.1093/mp/sst135
Chen, L., Zhou, J., Li, T., Fang, Z., Li, L., Huang, G., et al. (2021). GmoDetector: An accurate and efficient GMO identification approach and its applications. Food Res. Int. 149, 110662. doi: 10.1016/j.foodres.2021.110662
Cui, Y., Li, R., Li, G., Zhang, F., Zhu, T., Zhang, Q., et al. (2020). Hybrid breeding of rice via genomic selection. Plant Biotechnol. J. 18, 57–67. doi: 10.1111/pbi.13170
Duitama, J., Silva, A., Sanabria, Y., Cruz, D. F., Quintero, C., Ballen, C., et al. (2015). Whole genome sequencing of elite rice cultivars as a comprehensive information resource for marker assisted selection. PloS One 10, e0124617. doi: 10.1371/journal.pone.0124617
Ebana, K., Shibaya, T., Wu, J., Matsubara, K., Kanamori, H., Yamane, H., et al. (2011). Uncovering of major genetic factors generating naturally occurring variation in heading date among Asian rice cultivars. Theor. Appl. Genet. 122, 1199–1210. doi: 10.1007/s00122-010-1524-1
Endelman, J. B. (2011). Ridge regression and other kernels for genomic selection with r package rrBLUP. Plant Genome 4:250–255. doi: 10.3835/plantgenome2011.08.0024
Fujii, S., Kazama, T., Yamada, M., Toriyama, K. (2010). Discovery of global genomic re-organization based on comparison of two newly sequenced rice mitochondrial genomes with cytoplasmic male sterility-related genes. BMC Genomics 11, 209. doi: 10.1186/1471-2164-11-209
Gao, H., Jin, M., Zheng, X.-M., Chen, J., Yuan, D., Xin, Y., et al. (2014a). Days to heading 7, a major quantitative locus determining photoperiod sensitivity and regional adaptation in rice. Proc. Natl. Acad. Sci. 111, 16337–16342. doi: 10.1073/pnas.1418204111
Gao, H., Pirani, A., Webster, T., Shen, M. M. (2014b)Systems and methods for SNP characterization and identifying off target variants. In: US 20140274749 (United States). Available at: https://patents.google.com/patent/US11302417B2/en (Accessed August 12, 2022).
Grenier, C., Cao, T.-V., Ospina, Y., Quintero, C., Châtel, M. H., Tohme, J., et al. (2015). Accuracy of genomic selection in a rice synthetic population developed for recurrent selection breeding. PloS One 10, e0136594. doi: 10.1371/journal.pone.0136594
Guo, Z., Yang, Q., Huang, F., Zheng, H., Sang, Z., Xu, Y., et al. (2021). Development of high-resolution multiple-SNP arrays for genetic analyses and molecular breeding through genotyping by target sequencing and liquid chip. Plant Commun. 2, 100230. doi: 10.1016/j.xplc.2021.100230
Hori, K., Nonoue, Y., Ono, N., Shibaya, T., Ebana, K., Matsubara, K., et al. (2015). Genetic architecture of variation in heading date among Asian rice accessions. BMC Plant Biol. 15, 115. doi: 10.1186/s12870-015-0501-x
Huang, X., Qian, Q., Liu, Z., Sun, H., He, S., Luo, D., et al. (2009). Natural variation at the DEP1 locus enhances grain yield in rice. Nat. Genet. 41, 494–497. doi: 10.1038/ng.352
International Rice Genome Sequencing Project (2005). The map-based sequence of the rice genome. Nature 436, 793–800. doi: 10.1038/nature03895
IRRI (2022). Standard Evaluation System for Rice (SES), 4th Ed. (Los Banos, Philippines:International Rice Research Institute (IRRI))
Isidro, J., Jannink, J.-L., Akdemir, D., Poland, J., Heslot, N., Sorrells, M. E. (2015). Training set optimization under population structure in genomic selection. Theor. Appl. Genet. 128, 145–158. doi: 10.1007/s00122-014-2418-4
Kawahara, Y., de la Bastide, M., Hamilton, J. P., Kanamori, H., McCombie, W. R., Ouyang, S., et al. (2013). Improvement of the oryza sativa nipponbare reference genome using next generation sequence and optical map data. Rice 6, 4. doi: 10.1186/1939-8433-6-4
Kim, J.-H., Kim, S.-Y., Lee, H., Kim, Y.-R., Kim, H.-Y. (2010). An event-specific DNA microarray to identify genetically modified organisms in processed foods. J. Agric. Food Chem. 58, 6018–6026. doi: 10.1021/jf100351x
Koo, B.-H., Yoo, S.-C., Park, J.-W., Kwon, C.-T., Lee, B.-D., An, G., et al. (2013). Natural variation in OsPRR37 regulates heading date and contributes to rice cultivation at a wide range of latitudes. Mol. Plant 6, 1877–1888. doi: 10.1093/mp/sst088
Kutateladze, T., Bitskinashvili, K., Sapojnikova, N., Kartvelishvili, T., Asatiani, N., Vishnepolsky, B., et al. (2021). Development of multiplex PCR coupled DNA chip technology for assessment of endogenous and exogenous allergens in GM soybean. Biosensors 11, 481. doi: 10.3390/bios11120481
Lipka, A. E., Tian, F., Wang, Q., Peiffer, J., Li, M., Bradbury, P. J., et al. (2012). GAPIT: genome association and prediction integrated tool. Bioinformatics 28, 2397–2399. doi: 10.1093/bioinformatics/bts444
Liu, E., Liu, Y., Wu, G., Zeng, S., Tran Thi, T. G., Liang, L., et al. (2016). Identification of a candidate gene for panicle length in rice (Oryza sativa l.) Via association and linkage analysis. Front. Plant Sci. 7. doi: 10.3389/fpls.2016.00596
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, 10532. doi: 10.1038/ncomms10532
Morales, K. Y., Singh, N., Perez, F. A., Ignacio, J. C., Thapa, R., Arbelaez, J. D., et al. (2020). An improved 7K SNP array, the C7AIR, provides a wealth of validated SNP markers for rice breeding and genetics studies. PloS One 15, e0232479. doi: 10.1371/journal.pone.0232479
Nguyen, K. L., Grondin, A., Courtois, B., Gantet, P. (2019). Next-generation sequencing accelerates crop gene discovery. Trends Plant Sci. 24, 263–274. doi: 10.1016/j.tplants.2018.11.008
Onogi, A., Ideta, O., Inoshita, Y., Ebana, K., Yoshioka, T., Yamasaki, M., et al. (2015). Exploring the areas of applicability of whole-genome prediction methods for Asian rice (Oryza sativa l.). Theor. Appl. Genet. 128, 41–53. doi: 10.1007/s00122-014-2411-y
Pérez-Enciso, M., Zingaretti, L. M. (2019). A guide on deep learning for complex trait genomic prediction. Genes 10, 553. doi: 10.3390/genes10070553
Phitaktansakul, R., Kim, K.-W., Myo Aung, K., Zin Maung, T., Min, M.-H., Somsri, A., et al. (2021). Multi-omics analysis reveals the genetic basis of rice fragrance mediated by betaine aldehyde dehydrogenase 2. J. Advanced Res. doi: 10.1016/j.jare.2021.12.004
Rasheed, A., Hao, Y., Xia, X., Khan, A., Xu, Y., Varshney, R. K., et al. (2017). Crop breeding chips and genotyping platforms: Progress, challenges, and perspectives. Mol. Plant 10, 1047–1064. doi: 10.1016/j.molp.2017.06.008
Royston, P. (1995). Remark AS R94: A remark on algorithm AS 181: The W-test for normality. J. R. Stat. Society. Ser. C (Applied Statistics) 44, 547–551. doi: 10.2307/2986146
Sandhu, K. S., Lozada, D. N., Zhang, Z., Pumphrey, M. O., Carter, A. H. (2021). Deep learning for predicting complex traits in spring wheat breeding program. Front. Plant Sci. 11. doi: 10.3389/fpls.2020.613325
Schatz, M. C., Maron, L. G., Stein, J. C., Wences, A. H., Gurtowski, J., Biggers, E., et al. (2014). Whole genome de novo assemblies of three divergent strains of rice, oryza sativa, document novel gene space of aus and indica. Genome Biol. 15, 506. doi: 10.1186/s13059-014-0506-z
Seck, P. A., Diagne, A., Mohanty, S., Wopereis, M.C.S (2012).Crops that feed the world 7: Rice. C. S. Food Sec. 4, 7–24. doi: 10.1007/s12571-012-0168-1
Seo, J., Lee, G., Jin, Z., Kim, B., Chin, J. H., 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, 39. doi: 10.1007/s11032-020-01123-x
Sibo, C., Liang, T., Jian, S., Quan, X., Zhengjin, X., Wenfu, C. (2021). Contribution and prospect of erect panicle type to japonica super rice. Rice Sci. 28, 431–441. doi: 10.1016/j.rsci.2021.07.004
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
Spindel, J., Begum, H., Akdemir, D., Virk, P., Collard, B., Redoña, E., et al. (2015). Genomic selection and association mapping in rice (Oryza sativa): Effect of trait genetic architecture, training population composition, marker number and statistical model on accuracy of rice genomic selection in elite, tropical rice breeding lines. PloS Genet. 11, e1004982. doi: 10.1371/journal.pgen.1004982
The 3,000 rice genomes project (2014). The 3,000 rice genomes project. GigaScience 3, 7. doi: 10.1186/2047-217X-3-7
Thomson, M. J., Singh, N., Dwiyanti, M. S., Wang, D. R., Wright, M. H., Perez, F. A., et al. (2017). Large-Scale deployment of a rice 6 K SNP array for genetics and breeding applications. Rice 10, 40. doi: 10.1186/s12284-017-0181-2
Thomson, M. J., Zhao, K., Wright, M., 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
Tong, W., He, Q., Park, Y.-J. (2017). Genetic variation architecture of mitochondrial genome reveals the differentiation in Korean landrace and weedy rice. Sci. Rep. 7, 43327. doi: 10.1038/srep43327
Tong, W., Kim, T.-S., Park, Y.-J. (2016). Rice chloroplast genome variation architecture and phylogenetic dissection in diverse oryza species assessed by whole-genome resequencing. Rice 9, 57. doi: 10.1186/s12284-016-0129-y
Turkec, A., Lucas, S. J., Karacanli, B., Baykut, A., Yuksel, H. (2016). Assessment of a direct hybridization microarray strategy for comprehensive monitoring of genetically modified organisms (GMOs). Food Chem. 194, 399–409. doi: 10.1016/j.foodchem.2015.08.030
Verlouw, J. A. M., Clemens, E., de Vries, J. H., Zolk, O., Verkerk, A. J. M. H., am Zehnhoff-Dinnesen, A., et al. (2021). A comparison of genotyping arrays. Eur. J. Hum. Genet. 29, 1611–1624. doi: 10.1038/s41431-021-00917-7
Wang, Q., Tang, J., Han, B., Huang, X. (2020). Advances in genome-wide association studies of complex traits in rice. Theor. Appl. Genet. 133, 1415–1425. doi: 10.1007/s00122-019-03473-3
Wing, R. A., Purugganan, M. D., Zhang, Q. (2018). The rice genome revolution: from an ancient grain to green super rice. Nat. Rev. Genet. 19, 505–517. doi: 10.1038/s41576-018-0024-z
Xiao, N., Pan, C., Li, Y., Wu, Y., Cai, Y., Lu, Y., et al. (2021). Genomic insight into balancing high yield, good quality, and blast resistance of japonica rice. Genome Biol. 22, 283. doi: 10.1186/s13059-021-02488-8
Xu, Y., Ma, K., Zhao, Y., Wang, X., Zhou, K., Yu, G., et al. (2021). Genomic selection: A breakthrough technology in rice breeding. Crop J. 9, 669–677. doi: 10.1016/j.cj.2021.03.008
Xu, Y., Wang, X., Ding, X., Zheng, X., Yang, Z., Xu, C., et al. (2018). Genomic selection of agronomic traits in hybrid rice using an NCII population. Rice 11, 32. doi: 10.1186/s12284-018-0223-4
Xu, H., Zhao, M., Zhang, Q., Xu, Z., Xu, Q. (2016). The DENSE AND ERECT PANICLE 1 (DEP1) gene offering the potential in the breeding of high-yielding rice. Breed Sci. 66, 659–667. doi: 10.1270/jsbbs.16120
Xu, S., Zhu, D., Zhang, Q. (2014). Predicting hybrid performance in rice using genomic best linear unbiased prediction. Proc. Natl. Acad. Sci. 111, 12456–12461. doi: 10.1073/pnas.1413750111
Yan, W., Liu, H., Zhou, X., Li, Q., Zhang, J., Lu, L., et al. (2013). Natural variation in Ghd7.1 plays an important role in grain yield and adaptation in rice. Cell Res. 23, 969–971. doi: 10.1038/cr.2013.43
Yano, K., Morinaka, Y., Wang, F., Huang, P., Takehara, S., Hirai, T., et al. (2019). GWAS with principal component analysis identifies a gene comprehensively controlling rice architecture. Proc. Natl. Acad. Sci. 116, 21262–21267. doi: 10.1073/pnas.1904964116
Yonemaru, J.-i., Yamamoto, T., Fukuoka, S., Uga, Y., Hori, K., Yano, M. (2010). Q-TARO: QTL annotation rice online database. Rice 3, 194–203. doi: 10.1007/s12284-010-9041-z
Yu, S., Ali, J., Zhou, S., Ren, G., Xie, H., Xu, J., et al. (2022). From green super rice to green agriculture: Reaping the promise of functional genomics research. Mol. Plant 15, 9–26. doi: 10.1016/j.molp.2021.12.001
Yuan, J., Wang, X., Zhao, Y., Khan, N. U., Zhao, Z., Zhang, Y., et al. (2020). Genetic basis and identification of candidate genes for salt tolerance in rice by GWAS. Sci. Rep. 10, 1–9. doi: 10.1038/s41598-020-66604-7
Yu, J., Hu, S., Wang, J., Wong, G. K.-S., Li, S., Liu, B., et al. (2002). A draft sequence of the rice genome (Oryza sativa l. ssp. indica). Science 296, 79–92. doi: 10.1126/science.1068037
Yu, J., Pressoir, G., Briggs, W. H., Vroh Bi, I., Yamasaki, M., Doebley, J. F., et al. (2006). A unified mixed-model method for association mapping that accounts for multiple levels of relatedness. Nat. Genet. 38, 203–208. doi: 10.1038/ng1702
Yu, H., Xie, W., Li, J., Zhou, F., Zhang, Q. (2014). A whole-genome SNP array (RICE6K) for genomic breeding in rice. Plant Biotechnol. J. 12, 28–37. doi: 10.1111/pbi.12113
Zhang, J., Chen, L.-L., Sun, S., Kudrna, D., Copetti, D., Li, W., et al. (2016). Building two indica rice reference genomes with PacBio long-read and illumina paired-end sequencing data. Sci. Data 3, 160076. doi: 10.1038/sdata.2016.76
Zhang, F., Wang, C., Li, M., Cui, Y., Shi, Y., Wu, Z., et al. (2021). The landscape of gene–CDS–haplotype diversity in rice: Properties, population organization, footprints of domestication and breeding, and implications for genetic improvement. Mol. Plant 14, 787–804. doi: 10.1016/j.molp.2021.02.003
Zhao, K., Tung, C.-W., Eizenga, G. C., Wright, M. H., Ali, M. L., Price, A. H., et al. (2011). Genome-wide association mapping reveals a rich genetic architecture of complex traits in oryza sativa. Nat. Commun. 2, 467. doi: 10.1038/ncomms1467
Zhou, Y., Zhu, J., Li, Z., Yi, C., Liu, J., Zhang, H., et al. (2009). Deletion in a quantitative trait gene qPE9-1 associated with panicle erectness improves plant architecture during rice domestication. Genetics 183, 315–324. doi: 10.1534/genetics.109.102681
Keywords: SNP array, genotyping, genome selection, genome assisted breeding, GWAS
Citation: Kim K-W, Nawade B, Nam J, Chu S-H, Ha J and Park Y-J (2022) Development of an inclusive 580K SNP array and its application for genomic selection and genome-wide association studies in rice. Front. Plant Sci. 13:1036177. doi: 10.3389/fpls.2022.1036177
Received: 04 September 2022; Accepted: 30 September 2022;
Published: 24 October 2022.
Edited by:
Sung Don Lim, Sangji University, South KoreaReviewed by:
Bo-Keun Ha, Chonnam National University, South KoreaSun-Goo Hwang, Sangji University, South Korea
Copyright © 2022 Kim, Nawade, Nam, Chu, Ha and Park. 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: Yong-Jin Park, eWpwYXJrQGtvbmdqdS5hYy5rcg==; Jungmin Ha, ai5oYUBnd251LmFjLmty
†These authors have contributed equally to this work