Skip to main content

Genome-wide association study of eating and cooking qualities in different subpopulations of rice (Oryza sativa L.)



Starch and protein are two major components of polished rice, and the amylose and protein contents affect eating and cooking qualities (ECQs). In the present study, genome-wide association study with high-quality re-sequencing data was performed for 10 ECQs in a panel of 227 non-glutinous rice accessions and four derived panels.


Population structure accounted for high phenotypic variation in three routine panels and had minor effects on subspecies-based panels. Using the mixed linear model method based on the P + K model, we detected 29, 24, 16, 17, and 29 loci that were significant for ECQ parameters in each of the five panels. Some of these loci were close to starch synthesis-related genes. Two quantitative trait loci (QTLs) (chr.9: 15417525 ~ 15474876; 17538294 ~ 18443016) for several starch paste viscosity properties detected in four panels were close to the isoamylase 3 gene, one QTL (chr.1: 30627943 ~ 31668474) for consistency detected in three panels was close to the starch synthase IV-1 gene. The QTL (chr.7: 1118122 ~ 1967247) for breakdown (BD), detected in the whole panel and japonica panel, and one QTL (chr.7: 25312126 ~ 26540950) for BD and setback (SB), detected in the whole panel and indica panel, may be specific gene alleles in japonica or indica panels. One previously detected QTL (chr.11: 22240707 ~ 22563596) for protein content and one new QTL (chr.5: 7756614 ~ 8042699) for many ECQ traits detected in more than two panels, may represent valuable targets for future cloning of the underlying genes.


This study detected minor-effect QTLs affecting ECQs, and may increase our understanding of the genetic differences regulating the formation of ECQ between indica and japonica varieties.


Improving grain quality and grain yield are among the most important objectives in rice breeding programs. Apparent amylose content (AAC) is known to affect the eating and cooking quality (ECQ) of rice [1]. However, the rapid viscosity analyser (RVA) profile can distinguish eating quality rice varieties with similar AACs by evaluating the breakdown (BD), setback (SB), and consistency (CS) viscosities [2]. In general, rice varieties with good eating quality have a BD of more than 100 rapid visco units (RVU) and an SB of less than 25 RVU [3]. Therefore, starch paste viscosities play an essential role in estimating the eating, cooking, and processing qualities of rice, and have received more and more attention in rice breeding programs [4, 5].

The Wx gene encodes granule bound starch synthase I (GBSSI), which is responsible for amylose synthesis. Starch paste viscosity properties are predominantly controlled by the Wx gene, as well as minor-effect genes with various additive effects [68]. The starch paste viscosity of glutinous rice differs from that of non-glutinous rice due to its lack of amylose [1]. Since GBSSI is inactive in glutinous varieties, the starch paste properties are thought to be controlled by other genes. The pullulanase (PUL) gene plays an important role in the control of peak viscosity (PV), hot paste viscosity (HPV), cold paste viscosity (CPV), BD, peak time (PeT), and pasting temperature (PTemp) in glutinous rice [4]. Soluble starch synthase IIa (SSIIa) is responsible for HPV, CPV, BD, SB, PeT, and PTemp, as well as the thermal properties of glutinous rice [9]. Polymorphisms in both SBE1 and SBE3 loci account for 70 % of the observed variations in HPV and CPV, and account for 40 % in both PV and CS in waxy rice [10]. Many populations derived from bi-parental crosses with similar AACs were developed to identify QTLs for starch paste viscosities and eliminate the Wx effects in non-glutinous accessions [7, 11]. A QTL cluster close to SSIIa that had major effects on the alkali spreading value (ASV), pasting time, and minor effects on gel consistency, AAC, PV, CPV, BD, and SB was identified [7]. The branching enzyme IIb (BEIIb or SBE3) gene cluster is responsible for PV, HPV, CPV, BD, CS, and gel consistency [7] . Another method to avoid the effects of the Wx gene on starch paste viscosity was performed to detect more QTLs within sub-populations of the same Wx background [8, 12].

Association mapping is a powerful method to connect structural genomics and phenomics in plants, provided that information on population structure and linkage disequilibrium (LD) are available [13]. Association mapping, especially candidate-gene association mapping, is commonly used to detect QTLs for starch qualities [4, 9, 14, 15]. Tian et al. [14] found that the starch synthesis-related genes (SSRGs) cooperated with each other to form a regulatory network that controls the ECQs and defines correlations among AAC, ASV, and gel consistency. Yan et al. [4] genotyped waxy rice with 17 SSRGs using 43 gene-tagged molecular markers, and found that 10 SSRGs are involved in controlling the RVA profile parameters. However, candidate-gene association mapping will miss other unknown loci, and genome-wide association studies (GWAS) are comprehensive approaches to systemically search the genome for causal genetic variations [16]. The distribution of functional alleles is strongly correlated with population structure since LD can be caused by an admixture of subpopulations, which leads to false-positive results, particularly in populations with a small number of samples [17]. An ideal sample with a subtle population structure and familial relatedness is the best panel to control for false-negatives [17]. Moreover, subdividing the diversity panel to analyze subpopulations independently using the mixed model provides a reasonable solution to lower false-positive or -negative rates [18]. For example, Huang et al. [19] chose only the indica landrace population to conduct GWAS for 14 agronomic traits. Huang et al. [20] collected phenotypic data of flowering time and grain-related traits for a GWAS on the indica or japonica subpopulations and the whole population, respectively.

Two subspecies of rice (Oryza sativa L.), indica and japonica, differ in several morphological and physiological characteristics including grain shape, shattering, apiculus hair length, leaf color, seed dormancy, cold tolerance, phenol, and sensitivity to potassium chlorate [2123]. Moreover, japonica and indica rice varieties vary in starch physicochemical properties, which strongly influences rice cooking and processing methods for food and industrial applications [24]. The Wx a allele is predominant in non-waxy indica cultivars with high AACs (22–29 %), whereas the Wx b allele, with low AACs (12–19 %), is common to the non-waxy japonica variety [2527]. The splice donor mutation is prevalent in temperate japonica rice but rare or absent in indica and tropical japonica rice [28]. The indica type SSIIa can convert S-type amylopectin to L-type amylopectin by elongating short chains of DP (degree of polymerization) ≤10 to form chains of DP 13–25 [24]. The japonica S-type amylopectin has enriched short branch chains of DP ≤10 and depleted intermediate chains of 12 ≤DP ≤24 compared with the indica L-type amylopectin, whereas the proportion of cluster interconnecting long chains of DP ≥25 are comparable in the two varieties [29]. Bao et al. [30] classified 56 glutinous rice varieties into two groups according to the gelatinization temperature (GT), corresponding to high GT and low GT, and found that all 15 high GT accessions belonged to indica rice. Despite the two genes (Wx and SSIIa) being well-studied, the exact genetic effects of other SSRGs in shaping grain quality and starch paste viscosities in different rice subpopulations remain unclear [4].

In the present study, we divided the 227 non-glutinous rice accessions into two panels and two ecotype based panels, and explored the genetics of 10 ECQ parameters in different panels and ecotypes based on genome-wide association mapping with high-quality re-sequencing data. These results will increase our understanding of the mechanism of ECQ properties in different ecotypes and how minor-effect QTLs function on ECQs.


Phenotypic variation and contribution of population structure to phenotypic variation

Population structure was evaluated using the ADMIXTURE program, which estimates individual ancestry and admixture proportions assuming that K populations exist based on a maximum-likelihood method [31]. Because rice (Oryza sativa L.) has two subspecies, it was not surprising that indica and japonica subpopulations could be identified. However, the indica subpopulation could be further divided into indica, aus, as well as some admixtures (admix), while the japonica subpopulation could be divided into temperate japonica, tropical japonica, and some admixtures (Additional file 1: Figure S1). Principal component analysis (PCA) plots of the first two components (Fig. 1) clearly differentiated rice accessions into different subpopulations, which corresponded with the population structure results (Fig. 1; Additional file 1: Figure S1). All 10 phenotypic variations in the admixture group showed moderate levels. The aus group had the highest AAC, PeT, SB, and CS, but the lowest PV, HPV, CPV, BD, and PTemp. The indica group had a wide phenotypic variation range with regards to AAC, PV, HPV, CPV, BD, SB, and CS compared to the other groups, and its AAC, PV, CPV, SB, and CS mid-values were higher than those of the japonica group. Temperate japonica had the lowest AAC, SB, CS, and PeT, and the highest PTemp. The phenotypic variation of tropical japonica was between that of indica and temperate japonica. However, there was no significant difference in protein content (PRO) among groups (Fig. 2).

Fig. 1

Principal component analysis (first two components) in 5 panels

Fig. 2

Phenotypic variation among the different ecotype groups

With respect to the whole panel, a wide range in all traits was observed for AAC (10.2 ~ 30.7 %), PRO (1.6 ~ 14.1 %), PV (32.7 ~ 347.8RVU), HPV (24.2 ~ 245.3 RVU), CPV (59.6 ~ 371.5 RVU), BD (8.4 ~ 179.5 RVU), SB (−117.1 ~ 97.2 RVU), CS (28.6 ~ 170.6 RVU), PeT (3.7 ~ 6.9 min), and PTemp (68.1 ~ 81.5 °C) (Table 1). When comparing the indica panel with the japonica panel, the AAC, CPV, SB, CS, and PTemp of the indica panel were significantly higher than those of the japonica panel. The PRO, HPV, and PeT of the indica panel were significantly lower than those of the japonica panel, whereas PV and BD were similar between the two panels (Table 1). Furthermore, the AAC, PRO, SB, CS, and PTemp of Panel 1 were significantly higher than those of Panel 2. Panel 1 had lower PV, HPV, BD, and PeT than Panel 2 (Additional file 1: Table S1).

Table 1 Phenotypic variation in different panels and ecotype (indica and japonica panels) contribution to phenotypic variation

Given that population structure is a main factor affecting GWAS, we evaluated phenotypic variation among different subspecies groups (indica/japonica) in the whole panel (Fig. 2), as well as the population structure contribution for phenotypic variation in each panel (Table 1 and Additional file 1: Table S1).

In the whole panel, population structure explained 42.9 (AAC), 10.0 (PV), 12.2 (HPV), 8.4 (CPV), 7.8 (SB), 36.5 (CS), 27.7 (PeT), and 22.4 % (PTemp) of phenotypic variation (P <0.01) (Table 1). However, population structure affected only two traits in the indica panel (PV and HPV) and two traits in the japonica panel (AAC and CS) at the significance level of P <0.01, and affected 6 traits (AAC, PV, HPV, CPV, PeT, and PTemp) in indica panel and 3 traits (AAC, CS, and PTemp) in japonica panel at significance level of P <0.05 (Table 1). In Panel 1, except for BD, population structure accounted for large phenotypic variation explanation (PVE), ranging from 13.8 % (PTemp) to 57.4 % (AAC) at a significance level of P <0.001. In Panel 2, population structure explained 37.5 (PV), 46.0 (BD), 25.7 (SB), 22.6 (CS), and 25.3 % (PeT) phenotypic variation (P <0.01) (Additional file 1: Table S1).

GWAS of 10 ECQ parameters

There were 29, 24, 16, 17, and 29 loci identified in the whole panel, Panel 1, Panel 2, indica panel, and japonica panel, respectively (Tables 2, 3, and 4, Fig. 3).

Table 2 Loci detected for eating and cooking qualities in the whole panel
Table 3 Loci detected for eating and cooking qualities in Panel 1 and Panel 2
Table 4 Loci detected for eating and cooking qualities in indica and japonica panels
Fig. 3

Distribution of the significant loci for ECQs in five panels (The locus in the red circle indicates the Wx gene for AAC, PV, HPV, CPV, BD, SB, and CS in whole panel and Panel 1, for AAC and BD in panel 2; for AAC, PV, HPV, CPV, BD, and SB in indica panel, for AAC, PV, BD, and CS in japonica panel. Loci with white box, unmarked, underlines, black box, and brackets were detected in the whole panel, Panel 1, Panel 2, indica and japonica panels, respectively.)

Whole panel


Only one locus (chr.6: 1765761) for AAC was detected on chromosome 6. Two loci (chr.5: 13312843; chr.11: 22562237) for PRO were detected on chromosomes 5 and 11, accounting for 15.3 and 12.4 % PVE, respectively (Table 2).

PV, HPV, CPV, BD, SB, and CS

Two, four, three, three, four, and three loci were detected for PV, HPV, CPV, BD, SB, and CS, respectively. The major locus (chr.6: 1765761) was detected for PV, BD, SB, and CS, and the nearby locus (chr.6: 1750498) was the most significant site that linked with HPV and CPV. One locus (chr.5: 7756614) for PV was close to the locus (chr.5: 7944999) for HPV. One locus (chr.1: 31640680) for CPV was close to the locus (chr.1: 30634461) for CS. One locus (chr.9: 17538294) for HPV was close to one locus (chr.9: 17539188) for CPV. One locus (chr.7: 1967427) was detected for both BD and SB. One (chr.7: 26540950) for BD and one (chr.7: 25326697) for SB were close to each other (Table 2).

PeT and PTemp

A total of five loci located on chromosomes 1, 2, 5, 6, and 8, were detected for PeT. Two loci (chr.6: 7002095; chr.12: 17937512) for PTemp on chromosomes 6 and 12, accounted for 6.5 and 4.1 % of PVE, respectively (Table 2).

Panel 1 and Panel 2


One and two loci for AAC were detected in Panel 1 and Panel 2, respectively. The locus (chr.6: 1765761) and locus (chr.6: 1585864) were major QTLs for AAC in Panel1 and Panel 2, respectively. Another locus (chr.9:19310549) for AAC was detected only in Panel 2. Two loci for PRO, located on chromosomes 5 and 11, were identified in Panel 1. No locus was detected for PRO in Panel 2 (Table 3).

PV, HPV, CPV, BD, SB, and CS

In Panel 1, GWAS detected 2, 3, 4, 1, 1, and 4 loci for PV, HPV, CPV, BD, SB, and CS, respectively. Similar to the whole panel GWAS results, the locus (chr.6: 1765761) had major effects on BD, and SB, and CS, the locus (chr.6: 1750498) was detected for HPV and CPV, whereas the locus (chr.6: 1666386) was the most significant site for PV. The locus (chr.2: 29640874) for both PV and HPV accounted for 11.8 and 12.1 % PVE, respectively. One locus (chr.9: 15417525) for CPV and one (chr.9: 15474876) for CS was close to each other (Table 3).

In Panel 2, we detected 3, 4, 2, 2, and 1 locus for PV, HPV, CPV, BD, and SB, respectively. No locus was detected for CS. One locus (chr.5: 7944999) was significant for both PV and HPV. One locus (chr.1: 3541593) was detected for both HPV and CPV, while the locus (chr.7: 16503569) for PV was close to one locus (chr.7: 16506178) for both HPV and CPV on chromosome 7 (Table 3).

PeT and PTemp

In Panel 1, we detected four loci for PeT on chromosomes 2, 3, 6, and 10, while two loci (chr.6: 6757255 and chr.6: 21712093) for PTemp were detected on chromosome 6, accounting for 21.7 and 15.5 % PVE, respectively (Table 3).

In Panel 2, only one locus (chr.2: 24259382) for PeT and one locus (chr.8: 9459681) for PTemp were detected on chromosomes 2 and 8, respectively (Table 3).

indica and japonica panels


One locus (chr.6: 1529682) and one (chr.6: 1765761) for AAC were detected in indica and japonica panels, respectively. Three and six loci for PRO were detected in indica and japonica panels, respectively. One locus (chr.11: 22240707) for PRO in indica panel was close to the locus (chr.11: 22563596) detected in japonica panel (Table 4).

PV, HPV, CPV, BD, SB, and CS

In indica panel, one locus for PV, three for HPV, one for CPV, three for BD, and three for SB were detected and no locus was detected for CS. The locus (chr.6: 1750498) was a major QTL for both HPV and CPV, whereas loci for PV (chr.6: 1665483), BD (chr.6: 1721792), and SB (chr.6: 1617418) were close to each other in the indica panel. One locus (chr.9: 18729438) for HPV was close to the locus (chr.9: 18443016) for CPV in indica panel. One locus (chr.7: 26540950) for BD was close to the locus (chr.7: 25312126) for SB. The QTL (chr.11: 3735389) for BD was close to one QTL (chr.11: 3781133) for SB in the indica panel (Table 4).

Three loci for PV, four for HPV, one for CPV, two for BD, three for SB, and two for CS were identified in the japonica panel. The QTL (chr.6: 1765761) for PV and SB was close to the locus (chr.6: 1759451) for CS. One locus (chr.9: 18431810) for HPV was close to the locus (chr.9: 18443016) for CPV, one QTL (chr.1: 31774927) for SB was close to one QTL (chr.1: 31668474) for CS in japonica panel (Table 4).

PeT and PTemp

One locus (chr.1: 23829740) for PeT and one locus (chr.6: 6722646) for PTemp were detected in the indica panel, accounting for 32.3 and 29.5 % PVE, respectively.

Three loci for PeT and four for PTemp were detected in the japonica panel. The locus (chr.5: 7996572) for PeT were close to the locus detected in the whole panel. The locus (chr.6: 6747298) for PTemp was close to the loci in the whole panel, Panel 1, and the indica panel (Table 4).


Effect of population structure on GWAS

Given that linkage disequilibrium can be caused by the admixture of subpopulations and may result in false-positives if not correctly controlled in the statistical analysis, estimation of population structure must be included in the association analysis [17]. In this study, to determine whether different networks of alleles were associated with ECQ trait variation in the different panels, we performed GWAS on each of the four panels independently and in the panel as a whole. The present study allowed us to compare GWAS results derived from different panels.

The effect of population structure on phenotypic variation was estimated based on analysis of variance (ANOVA) (Table 1 and Additional file 1: Table S1). For the whole panel, eight of 10 traits were strongly affected by population structure. For Panel 1 and 2, nine and five traits were strongly affected by population structure (P <0.01), respectively, indicating that Panel 2 (breeding lines) experienced artificial selection, which was less structured than Panel 1. However, only two traits were significantly affected by population structure in the indica and japonica panels at significance level of P <0.01, suggesting that association mapping in the subpopulations, e.g. indica or japonica subpopulations, with simple population structures can lower false-positives. indica and japonica varietal groups should be properly treated for association analysis, which may explain why transgressive offspring occur in divergent subpopulation crosses [18] .

QTL comparison among different panels

For convenience, we regarded those SNPs that were detected in different panels with the physical distance less than 2.0 Mb as the same QTL, and these SNPs were regarded as alleles.

Similar to previous studies, the Wx gene is responsible for AAC in these five panels. The single nucleotide polymorphism (SNP) (chr.6: 1765761) for AAC is the Wx a/Wx b allele (a SNP at the splice site of intron 1), which was present in the whole panel, Panel 1, and japonica panel. The Wx gene was also the major QTL for AAC in the other two panels, but the Wx a/Wx b allele was not the most significant association site, suggested that different traits maybe preferred linked with different SNPs. Besides the Wx gene, one minor effect locus (chr.9: 19310549) for AAC was detected in Panel 2 (consisted of breeding lines); this locus was close to the QTL flanked by BPH and RM160 [32] (Fig. 3).

One QTL (chr.5: 12559930 ~ 13312843) for PRO detected in the whole panel and Panel 1 was close to the protein content QTL qPRO5, which was detected by Xu et al. [8] in the M201-Wx subpopulation. The other QTL (chr.11: 22562237 ~ 22563596) for PRO detected in the whole, Panel 2, indica, and japonica panels was close to pro11, which is flanked by RM209 and RM229 [33].

For starch paste viscosity parameters, a specific SNP (chr.6: 1765761), recognized as the Wx a/Wx b allele, was responsible for 4 (PV, BD, SB, and CS), 3 (BD, SB, and CS), and 2 (PV and SB) in the whole panel, Panel 1, and japonica panel, respectively. There were also some SNPs close to Wx gene that significantly associated with RVA parameters. However, only one Wx gene linked SNP was detected for BD in Panel 2, that may due to the reason that most Korean breeding lines were low-AAC varieties with less diversity. The GWAS results in our study may increase our understanding of the genetic control of starch properties in South Korea japonica rice varieties with low AACs and the same Wx allele. Besides the Wx gene, 12 QTLs were identified in two panels, 3 QTLs were detected in 3 panels, and 2 QTLs were detected in 4 panels (Fig. 3). The QTL (chr.5: 7756614 ~ 8042699) was detected for PV in the whole panel and Panel 2, for HPV in the whole panel, Panel 2 and japonica panel, and for PeT in the whole panel and japonica panel. This locus has not been detected in previous studies, and is a potentially novel locus specifically present in japonica rice varieties. The QTL (chr.9: 17538294 ~ 18443016) for HPV in the whole panel, indica, and japonica panels, and for CPV in the whole panel and japonica panel, was close to qPV9 for PV, HPV, and CPV [8], which is close to the isoamylase 3 gene (Os09g0469400, chr.9: 17851431 ~ 17868668). The QTL (chr.7: 1118122 ~ 1967427) for BD in the whole panel and japonica panel and SB in the whole panel, was close to qSB7 [7]. The QTL (chr.7: 2532697 ~ 26540950) for BD and SB in the whole panel and indica panel was close to the QTL (linked with RM6420) for PV, HPV, CPV, and SB [34], and close to the QTL for SB and CS [8]. The QTL (chr.4: 24056499 ~ 24293919) for SB, detected in the whole panel and Panel 2, was close to the QTL for FV (final viscosity, CPV), PKT, SBV, and TV (HPV), and is linked with RM3785 [35]. The QTL (chr.1: 30627943 ~ 31774927) for CS in the whole panel, Panel 1, and japonica panel, and three loci for SB (japonica panel), PeT (whole panel), and CPV (whole panel) was 0.21 Mb from the SSIV-1 gene (Os01g0720600, chr.1:30032428 ~ 30041425), indicating that SSIV-1 may influence starch qualities through affecting CPV, SB, CS, and PeT. The QTL (chr.9: 15417525 ~ 15474876) for CS in the whole panel and Panel 1, and for CPV in Panel 1, was close to qCS9 described by Yan et al. [7]. The QTL (chr.3: 33737410 ~ 33893179) for PeT was detected in Panel 1 and the japonica panel, and was close to Btemp (peak temperature) flanked by RM203 and RM422 [15]. The QTL (chr.6: 6722646 ~ 7002095) for PTemp, detected in the whole panel, Panel 1, indica, and japonica panels, was close to the ALK gene (Os06g0229800, chr.6: 6748398 ~ 6753302) [36], which encodes starch synthase IIa (SSIIa).

QTLs detected in specific panels

Panel 1

The locus (chr.1: 27287504) for HPV, detected only in the Panel 1, was about 1.9 Mb away from the ADP-glucose pyrophosphorylase large subunit 2 (AGPL2) gene (Os01g0633100, chr.1:25353805 ~ 25361883). The locus (chr.2: 29640874) for PV and CPV was about 1.6 Mb from the starch synthase II-2 (SSII-2 or SSIIb) gene (Os02g0744700, chr.2:31233292 ~ 31238210). The locus (chr.11: 24255295) for CS was close to qSB11b (Yan et al. 2014). The locus (chr.4: 27772775) for CPV was close to the QTL for PV and FV, which is linked with RM255 [35].

Panel 2

The QTL (chr.7: 16503569 ~ 16506178) for PV, HPV, and CPV was close to qHPV7 and qCPV7, which are linked with RM418 [37]. The locus (chr.1:3541593) for HPV and CPV, and locus (chr.1: 5102566) for PV were close to the QTL for fa (DP <12) and HPV, which are linked with RM7278 [35]. The locus (chr.2: 8596549) for BD was close to qBD2 [7]. The locus (chr.2: 24259382) for PeT was close to the qBTemp2 for peak temperature [7].

indica panel

The locus (chr.3: 26688241) for HPV was close to the qHPV3.2 [8].

japonica panel

The locus (chr.4: 6536977) for PRO was close to the qPC-4 [38]. The locus (chr.8: 2885700) for PV was close to the qPV8 [39] and qPKV8 [7]. The QTL (chr.2: 4054094) for PeT was close to the qPeT2 [34].


Overall, we investigated the genetics of 10 rice ECQ parameters by GWAS mapping in the whole panel and four derived panels. GWAS within subspecies panels (indica/japonica) cannot only lower false-positives influenced by population structure, but can also increase our understanding of ECQ differences between indica and japonica subspecies. Comparison of GWAS results indicated that some loci were common in different panels, while some were in specific panels. Loci detected in more than one panel were potentially reliable QTLs for ECQs, and these QTLs should be properly used in molecular breeding to improve rice ECQs, particularly in the Korean japonica rice with similar AACs and the same Wx allele.



A total of 227 non-glutinous rice accessions from 295 rice accessions (Additional file 2: Table S2) and their whole genome re-sequencing data [40] were used in this study. Here, we defined the 227 accessions as a whole panel for GWAS, among which the core collection of 110 diverse rice accessions was defined as Panel 1, which comprises 57 Korean accessions and 53 rice accessions collected form 27 countries worldwide. The remaining 117 accessions were considered as Panel 2, which consists of breeding lines, one USDA rice accession, three Japan breeding lines, and 113 Korean breeding lines released between 1967 and 2011. Furthermore, all 227 accessions were divided into the indica panel (59) and japonica panel (168), since indica/japonica species contribute the most to the population structure. The whole panel was planted in 2014 during the growing season from early May to September at Kongju National University farm, Yeasan Gun, South Korea. A randomized block design with two replications for each line was used in field experiments.

After being air-dried and stored at room temperature for 3 months, all rice samples were de-hulled to brown rice, then milled into white rice, afterwards ground and passed through a 100-mesh sieve to obtain milled-grain flour. All samples were stored in cooling rooms at 4 °C until use.

Phenotypic data evaluation

AAC and protein content

Determination of AAC was performed using the iodine staining method [41]. The absorbance of the solution was measured at 620 nm against the blank solution with a spectrophotometer. AAC was calculated using a standard curve made from rice samples with known AACs. PRO was measured using the Kjeldahl method with a Kjeltec-Foss 70 2400 Auto analyzer using 5.95 as the nitrogen-to-protein conversion factor.

RVA profile parameters

The RVA profile of rice flour was determined using a Rapid Visco Analyser (RVA-3, Newport Scientific, Warrie wood, Australia). According to the method of AACC61-02, 3.0 g of rice flour (12 % m. b.) was placed in an aluminum canister, after which 25.0 g of distilled water was added. The heating profile was set as follows: (1) the temperature was held at 50 °C for 1 min; (2) the temperature was linearly ramped up to 95 °C until 4.8 min; (3) the temperature was held at 95 °C until 7.3 min; (4) the temperature was linearly ramped down to 50 °C at 11.1 min; and (5) held at 50 °C until 12.5 min.

Four primary parameters were obtained from the pasting curve: PV, HPV, CPV, and PeT. Three secondary parameters were calculated from the primary parameters: BD (PV minus HPV), SB (CPV minus PV), and CS (CPV minus HPV). The viscosity parameters were measured in RVU. PTemp was determined according to the method proposed by Bao [42].

Statistical analyses of phenotypic data

Analysis of variance (ANOVA) was performed on the different panels to determine the phenotypic variation and population structure influence on phenotypic variation using the general linear model procedure (PROC GLM) by the SAS program (version 9.3, SAS Institute Inc., Cary, NC). In the model, population structure was set as the fixed independent variable, and 10 phenotypic data points were set as the dependent variables.

Population structure evaluation and principal component analysis

The genome re-sequencing data of the 227 accessions were collected from the 295 rice whole genome re-sequencing data, which were previously obtained using HiSeq 2000 and HiSeq 2500 [40, 43]. After removing SNPs with a minimum frequency <0.05 by filter alignment using TASSEL 5.0, there were 2,071,925, 1,229,263, 1,563,570, 1,178,226, and 574,169 high-quality SNPs remaining for Panel 1, Panel 2, the whole panel, the indica panel, and the japonica panel, respectively. Population structure was evaluated by the ADMIXTURE program, which estimates individual ancestry and admixture proportions assuming that K populations exist based on a maximum-likelihood method [31]. We ran the program with the script:

admixture --cv genotype.bed k

k ranged from 1 to 5. Principal component analysis (PCA) of each panel was calculated by TASSEL 5.0 with default parameters by using high-quality SNPs.

Genome-wide association mapping of five panels

To control for false-positives, the P (PCA) + K model, based on a mixed linear model (MLM), was used for genome-wide association mapping for the five panels by GAPIT software [44]. Here the principal components (PCs) were used to control for population structure in GWAS. We used the default parameters for running the GAPIT software: GAPIT < − myGAPIT(Y = myY, G = myG, = 3, Model.selection = TRUE). myY is the phenotype matrix and myG is genotype matrix. SNP sites with the lowest P value in the peak region (P <10−4) were considered significant SNPs (or loci) for phenotypic variation. Otherwise, sites with regions having only one SNP significant were considered false-positives. Loci with a physical distance less than 2.0 Mb were considered the same locus.


AAC, apparent amylose content; BD, breakdown; CPV, cold paste viscosity; CS, consistency; ECQ, eating and cooking quality; GWAS, genome-wide association study; HPV, hot paste viscosity; PeT, peak time; PRO, protein content; PTemp, pasting temperature; PV, peak viscosity; PVE, phenotypic variation explanation; QTL, quantitative trait locus; RVA, rapid viscosity analyser; RVU, rapid visco units; SB, setback; SNP, single nucleotide polymorphism; SSRGs, starch synthesis related genes.


  1. 1.

    Bao JS. Toward understanding the genetic and molecular bases of the eating and cooking qualities of rice. Cereal Food World. 2012;57(4):148–56.

    CAS  Article  Google Scholar 

  2. 2.

    Bao JS. Progress in studies on inheritance and improvement of rice starch quality. Mol Plant Breed. 2007;5(6):1–20.

    CAS  Google Scholar 

  3. 3.

    Shu Q, Wu D, Xia Y, Gao M. Relationship between RVA profile character and eating quality in Oryza sativa L. Scientia Agricultura Sinica. 1998;31(3):25–9.

    CAS  Google Scholar 

  4. 4.

    Yan CJ, Tian ZX, Fang YW, Yang YC, Li JA, Zeng SY, Gu SL, Xu CW, Tang SZ, Gu MH. Genetic analysis of starch paste viscosity parameters in glutinous rice (Oryza sativa L.). Theor Appl Genet. 2011;122(1):63–76.

    CAS  Article  PubMed  Google Scholar 

  5. 5.

    He Y, Han Y, Jiang L, Xu C, Lu J, Xu M. Functional analysis of starch-synthesis genes in determining rice eating and cooking qualities. Mol Breeding. 2006;18(4):277–90.

    CAS  Article  Google Scholar 

  6. 6.

    Bao JS, Zheng XW, Xia YW, He P, Shu QY, Lu X, Chen Y, Zhu LH. QTL mapping for the paste viscosity characteristics in rice (Oryza sativa L.). Theor Appl Genet. 2000;100(2):280–4.

    CAS  Article  Google Scholar 

  7. 7.

    Yan B, Tondi Yacouba N, Chen J, Wang Y, Gao G, Zhang Q, Liu X, He Y. Analysis of minor quantitative trait loci for eating and cooking quality traits in rice using a recombinant inbred line population derived from two indica cultivars with similar amylose content. Mol Breeding. 2014;34(4):2151–63.

    CAS  Article  Google Scholar 

  8. 8.

    Xu F, Sun C, Huang Y, Chen Y, Tong C, Bao JS. QTL mapping for rice grain quality: a strategy to detect more QTLs within sub-populations. Mol Breeding. 2015;35(4):1–11.

    Article  Google Scholar 

  9. 9.

    Xu F, Zhang G, Tong C, Sun X, Corke H, Sun M, Bao JS. Association mapping of starch physicochemical properties with starch biosynthesizing genes in waxy rice (Oryza sativa L.). J Agr Food Chem. 2013;61(42):10110–7.

    CAS  Article  Google Scholar 

  10. 10.

    Han YP, Xu ML, Liu XY, Yan CJ, Korban SS, Chen XL, Gu MH. Genes coding for starch branching enzymes are major contributors to starch viscosity characteristics in waxy rice (Oryza sativa L.). Plant Sci. 2004;166(2):357–64.

    CAS  Article  Google Scholar 

  11. 11.

    Bao JS, Wu YR, Hu B, Wu P, Cui HR, Shu QY. QTL for rice grain quality based on a DH population derived from parents with similar apparent amylose content. Euphytica. 2002;128(3):317–24.

    CAS  Article  Google Scholar 

  12. 12.

    Yang F, Chen Y, Tong C, Huang Y, Xu F, Li K, Corke H, Sun M, Bao JS. Association mapping of starch physicochemical properties with starch synthesis-related gene markers in nonwaxy rice (Oryza sativa L.). Mol Breeding. 2014;34(4):1747–63.

    CAS  Article  Google Scholar 

  13. 13.

    Jin L, Lu Y, Xiao P, Sun M, Corke H, Bao JS. Genetic diversity and population structure of a diverse set of rice germplasm for association mapping. Theor Appl Genet. 2010;121(3):475–87.

    Article  PubMed  Google Scholar 

  14. 14.

    Tian Z, Qian Q, Liu Q, Yan M, Liu X, Yan C, Liu G, Gao Z, Tang S, Zeng D, et al. Allelic diversities in rice starch biosynthesis lead to a diverse array of rice eating and cooking qualities. Proc Natl Acad Sci U S A. 2009;106(51):21760–5.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  15. 15.

    Wang L, Liu W, Xu Y, He Y, Luo L, Xing Y, Xu C, Zhang Q. Genetic basis of 17 traits and viscosity parameters characterizing the eating and cooking quality of rice grain. Theor Appl Genet. 2007;115(4):463–76.

    CAS  Article  PubMed  Google Scholar 

  16. 16.

    Zhu C, Gore M, Buckler ES, Yu J. Status and prospects of association mapping in plants. Plant Genome. 2008;1(1):5–20.

    CAS  Article  Google Scholar 

  17. 17.

    Yu J, Buckler ES. Genetic association mapping and genome organization of maize. Curr Opin Chem Biol. 2006;17(2):155–60.

    CAS  Google Scholar 

  18. 18.

    Zhao K, Tung C-W, Eizenga GC, Wright MH, Ali ML, Price AH, Norton GJ, Islam MR, Reynolds A, Mezey J, et al. Genome-wide association mapping reveals a rich genetic architecture of complex traits in Oryza sativa. Nat Commun. 2011;2:467.

    Article  PubMed  PubMed Central  Google Scholar 

  19. 19.

    Huang X, Wei X, Sang T, Zhao Q, Feng Q, Zhao Y, Li C, Zhu C, Lu T, Zhang Z, et al. Genome-wide association studies of 14 agronomic traits in rice landraces. Nat Genet. 2010;42(11):961–U976.

    CAS  Article  PubMed  Google Scholar 

  20. 20.

    Huang X, Zhao Y, Wei X, Li C, Wang A, Zhao Q, Li W, Guo Y, Deng L, Zhu C, et al. Genome-wide association study of flowering time and grain yield traits in a worldwide collection of rice germplasm. Nat Genet. 2012;44(1):32–9.

    Article  Google Scholar 

  21. 21.

    Lin Z, Griffith ME, Li X, Zhu Z, Tan L, Fu Y, Zhang W, Wang X, Xie D, Sun C. Origin of seed shattering in rice (Oryza sativa L.). Planta. 2007;226(1):11–20.

    CAS  Article  PubMed  Google Scholar 

  22. 22.

    Kovach MJ, Sweeney MT, McCouch SR. New insights into the history of rice domestication. Trends Genet. 2007;23(11):578–87.

    CAS  Article  PubMed  Google Scholar 

  23. 23.

    Oka HI, Morishima H. Phylogenetic differentiation of cultivated rice.23. potentiality of wild progenitors to evolve the indica and japonica types of rice cultivars. Euphytica. 1982;31(1):41–50.

    Article  Google Scholar 

  24. 24.

    Nakamura Y, Francisco PB, Hosaka Y, Sato A, Sawada T, Kubo A, Fujita N. Essential amino acids of starch synthase IIa differentiate amylopectin structure and starch quality between japonica and indica rice varieties. Plant Mol Biol. 2005;58(2):213–27.

    CAS  Article  PubMed  Google Scholar 

  25. 25.

    Yamanaka S, Nakamura I, Watanabe KN, Sato YI. Identification of SNPs in the waxy gene among glutinous rice cultivars and their evolutionary significance during the domestication process of rice. Theor Appl Genet. 2004;108(7):1200–4.

    CAS  Article  PubMed  Google Scholar 

  26. 26.

    Hirano HY, Eiguchi M, Sano Y. A single base change altered the regulation of the Waxy gene at the posttranscriptional level during the domestication of rice. Genome Biol Evol. 1998;15(8):978–87.

    CAS  Article  Google Scholar 

  27. 27.

    Wang ZY, Zheng FQ, Shen GZ, Gao JP, Snustad DP, Li MG, Zhang JL, Hong MM. The amylose content in rice endosperm is related to the posttranscriptional regulation of the Waxy gene. Plant J. 1995;7(4):613–22.

    CAS  Article  PubMed  Google Scholar 

  28. 28.

    Olsen KM, Caicedo AL, Polato N, McClung A, McCouch S, Purugganan MD. Selection under domestication: Evidence for a sweep in the rice Waxy genomic region. Genetics. 2006;173(2):975–83.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  29. 29.

    Nakamura Y, Sakurai A, Inaba Y, Kimura K, Iwasawa N, Nagamine T. The fine structure of amylopectin in endosperm from Asian cultivated rice can be largely classified into two classes. Starch-Starke. 2002;54(3–4):117–31.

    CAS  Article  Google Scholar 

  30. 30.

    Bao JS, Corke H, Sun M. Genetic diversity in the physicochemical properties of waxy rice (Oryza sativa L) starch. J Sci Food Agr. 2004;84(11):1299–306.

    CAS  Article  Google Scholar 

  31. 31.

    Alexander DH, Novembre J, Lange K. Fast model-based estimation of ancestry in unrelated individuals. Genome Res. 2009;19(9):1655–64.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  32. 32.

    Li X, Ma W, Wang Z, Guo Y, Liu F, Wu H. Study on integration of QTLs related to amylose content in rice. Journal of Northeast Agricultural University. 2014;45(3):8–14.

    Google Scholar 

  33. 33.

    Aluko G, Martinez C, Tohme J, Castano C, Bergman C, Oard J. QTL mapping of grain quality traits from the interspecific cross Oryza sativa × O. glaberrima. Theor Appl Genet. 2004;109(3):630–9.

    CAS  Article  PubMed  Google Scholar 

  34. 34.

    Hsu YC, Tseng MC, Wu YP, Lin MY, Wei FJ, Hwu KK, Hsing YI, Lin YR. Genetic factors responsible for eating and cooking qualities of rice grains in a recombinant inbred population of an inter-subspecific cross. Mol Breeding. 2014;34(2):655–73.

    CAS  Article  Google Scholar 

  35. 35.

    Zhao WG, Chung JW, Kwon SW, Lee JH, Ma KH, Park YJ. Association analysis of physicochemical traits on eating quality in rice (Oryza sativa L.). Euphytica. 2013;191(1):9–21.

    Article  Google Scholar 

  36. 36.

    Gao Z, Zeng D, Cui X, Zhou Y, Yan M, Huang D, Li J, Qian Q. Map-based cloning of the ALK gene, which controls the gelatinization temperature of rice. Sci China Ser C. 2003;46(6):661–8.

    CAS  Article  Google Scholar 

  37. 37.

    Cho YC, Suh JP, Yoon MR, Baek MK, Won YJ, Lee JH, Park HS, Baek SH, Lee JH. QTL mapping for paste viscosity characteristics related to eating quality and QTL-NIL development in japonica rice (Oryza sativa L.). Plant Breed Biot. 2013;1(4):333–46.

    Article  Google Scholar 

  38. 38.

    Zheng L, Zhang W, Chen X, Ma J, Chen W, Zhao Z, Zhai H, Wan J. Dynamic QTL analysis of rice protein content and protein index using recombinant inbred lines. J Plant Biol. 2011;54(5):321–8.

    CAS  Article  Google Scholar 

  39. 39.

    Kobayashi A, Tomita K. QTL detection for stickiness of cooked rice using recombinant inbred lines derived from crosses between japonica rice cultivars. Breeding Sci. 2008;58(4):419–26.

    CAS  Article  Google Scholar 

  40. 40.

    He Q, Yu J, Kim TS, Cho YH, Lee YS, Park YJ. Resequencing reveals different domestication rate for BADH1 and BADH2 in rice (Oryza sativa). PLoS One. 2015;10(8):e0134801.

    Article  PubMed  PubMed Central  Google Scholar 

  41. 41.

    Bao JS, Shen S, Sun M, Corke H. Analysis of genotypic diversity in the starch physicochemical properties of nonwaxy rice: apparent amylose content, pasting viscosity and gel texture. Starch‐Stärke. 2006;58(6):259–67.

    CAS  Article  Google Scholar 

  42. 42.

    Bao JS. Accurate measurement of pasting temperature by the rapid visco-analyser: a case study using rice flour. Rice Sci. 2008;15(1):69–72.

    Article  Google Scholar 

  43. 43.

    Kim TS, He Q, Kim KW, Yoon MY, Ra WH, Li FP, Tong W, Yu J, Win HO. Choi. et al. Genome-wide resequencing of KRICE_CORE reveals their potential for future breeding, as well as functional and evolutionary studies in the post-genomic era. BMC Genomics. 2016;17(1):1.

    CAS  Google Scholar 

  44. 44.

    Lipka AE, Tian F, Wang Q, Peiffer J, Li M, Bradbury PJ, Gore MA, Buckler ES, Zhang Z. GAPIT: genome association and prediction integrated tool. Bioinformatics. 2012;28(18):2397–9.

    CAS  Article  Google Scholar 

Download references


We thank Dr. Kyu-Won Kim for his technical supports on data analysis.


This work was carried out with the support of “Cooperative Research Program for Agriculture Science & Technology Development (Project No. PJ01116101)” Rural Development Administration, Republic of Korea. This research was supported by Bio-industry Technology Development Program, Ministry of Agriculture, Food and Rural Affairs (115078–2).

Availability of data and material

The datasets generated during and/or analysed during the current study are available from the corresponding author on reasonable requests. Additional file 2: Table S2 was added to show the source of the rice accessions used in this study, and the phylogenetic data was uploaded to the Tree base (

Authors’ contributions

JB and YP conceived and designed the study; FX measured phenotype data; QH performed data analysis; FX, JB and YP drafted the manuscript. All authors critically reviewed and approved the final manuscript.

Competing interests

The authors declare that they have no competing interests.

Consent for publication

Not applicable

Ethics approval and consent to participate

Not applicable

Author information



Corresponding author

Correspondence to Yong-Jin Park.

Additional files

Additional file 1:

Table S1. Phenotypic variation in different panels and contribution of population structure to phenotypic variation. Figure S1 Population structure of the whole panel. (DOCX 6297 kb)

Additional file 2:

Table S2. The sources of the rice varieties used in this study. (DOC 32 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Xu, F., Bao, J., He, Q. et al. Genome-wide association study of eating and cooking qualities in different subpopulations of rice (Oryza sativa L.). BMC Genomics 17, 663 (2016).

Download citation


  • Rice
  • Eating and cooking quality
  • Subpopulations
  • GWAS