QTL Identification for Cooking and Eating Quality in indica Rice Using Multi-Parent Advanced Generation Intercross (MAGIC) Population

Association mapping using a multi-parent advanced generation intercross (MAGIC) population provides a promising tool in genetic dissection of rice cooking and eating quality (CEQ). In this study, QTLs were identified for ten physicochemical properties related to CEQ using 508 F₆ MAGIC lines. The whole population and eight founder lines were genotyped with 6K Illumina Infinium HD Assay. All traits had high heritability estimates and showed a large genetic variation in the MAGIC population. Highly significant phenotypic correlations were present between traits. AC was significantly positively correlated with PKT, TV, FV, SBV, PKT, and RT but significantly negatively correlated with GC and BDV. Seventeen QTLs were identified for all traits. GBSSI locus was hosted or closely to nine QTLs, qAC6, qGC6.1, qPKT6.1, qPKV6, qBDV6.1, qTV6.1, qFV6, qSBV6, and qRT6, suggesting that GBSSI impacts the overall CEQ. Another locus closed to SSIIa, located at 6.99 Mb, affects five traits, GC, PKT, BDV, SBV, and PT. The identified QTLs revealed small to modest effects where the highest percentage of phenotypic variance explained was 17.18%. These QTLs are directly relevant and useful in breeding for CEQ in indica rice. These results also confirmed that QTL mapping via association mapping using a MAGIC population is a powerful method in genetic analysis of complex traits.

Introduction

Grain quality is one of the most important quantitative traits in rice. It is a multi-faceted trait involving physical and biochemical aspects relating to milling, appearance, cooking and eating quality (CEQ), and nutrition (He et al., 1999; Fitzgerald et al., 2009; Bao, 2014). Among these, CEQ is the most important aspect especially in consumers’ perspective serving as the basis of choice for quality rice. CEQ mainly affects the easiness of cooking, firmness, and cohesiveness of cooked rice and CEQ is directly related to four physicochemical properties, the amylose content (AC), gel consistency (GC), gelatinization temperature (GT), and pasting viscosity (PV) (Dela Cruz and Khush, 2000; Fitzgerald et al., 2009; Bao, 2014). Hence, most of the genetics studies use these traits as rice quality parameters.

Amylose content is the most important chemical property affecting CEQ, which is considered as a good indicator of the appearance and texture of cooked rice (Lanceras et al., 2000). GC is a measure of firmness and stickiness of cooked rice upon cooling and therefore, a measure of strength of the gel. It is used as one of the key parameters to determine the textural property of cooked rice (Cuevas and Fitzgerald, 2012). It is also used to differentiate rice varieties having similar AC class, since GC is highly associated with AC (i.e., rice variety with low AC tend to have high GC or soft gel upon cooling) (Dela Cruz and Khush, 2000; Cuevas and Fitzgerald, 2012). Thus, GC is negatively correlated with AC. Pasting viscosity is widely used to assess starch properties in rice (Bao and Xia, 1999; Yan et al., 2005; Allahgholipour et al., 2006). It is measured using rapid visco analyzer (RVA), a rotational viscometer, which mimics the cooking process (Bao, 2014). To some extent, PV is interchangeably used with gelatinization but the two parameters are different. Gelatinization specifically refers to molecular disintegration of starch molecules, whereas PV refers to starch viscosity development. This suggests that PV occurs after gelatinization.

Several studies have reported major genes/QTLs for CEQ in rice using DH population (He et al., 1999; Bao et al., 2000), RILs (Bao et al., 2003), and CSSLs (Zhang et al., 2013) via bi-parental linkage mapping. Among these genes identified, it is well known that GBSSI gene is the major gene responsible for CEQ (He et al., 1999; Tran et al., 2011). Most of these QTLs detected had very high phenotypic variance explained, however, these QTLs were often not transferable to other genetic backgrounds since the estimated effects are limited to the population studied (Holland, 2007). Hence, development of rice varieties with superior CEQ is confined to very few genetic backgrounds.

Genetic analysis via association mapping (AM) offers a great opportunity to identify significant marker-trait associations (MTAs) (Huang et al., 2012). AM is a practical application of linkage disequilibrium (LD) mapping since its principle lies in LD. Hence, LD greatly influence the resolution of mapping QTL associated with the trait under study. The use of natural accessions for AM results in high mapping resolution due to its great LD decay over time brought by years of numerous recombination. However, there is an extensive population structure and cryptic kinship in natural accessions which increases the chance of discovering false positive MTAs (Zhu et al., 2008; Huang et al., 2015). Multi-parent advanced generation intercross (MAGIC) population, an extension of advanced intercrossed lines (Darvasi and Soller, 1995), offers an alternative type of mapping population that overcomes the limitations of family-based mapping using populations derived from biparental cross and AM using natural accessions. It involves the use of multiple genetically diverse parents to come up with a population of high allelic and phenotypic diversity in combination with high levels of recombination events brought about by several cycles of inter-mating (Mackay and Powell, 2007; Cavanagh et al., 2008; Huang et al., 2015). Increased recombination events reduce the extent of LD and thus increase genetic mapping resolution (Huang et al., 2015). Furthermore, since the founder lines contributing to the allelic and phenotypic diversity are well known, the population structure and kinship are highly controlled. This suggests that MAGIC population offers a great potential for mapping complex quantitative traits by preventing spurious MTAs. The effective use of MAGIC population in QTL identification via AM has been conducted in several cereal crops such as rice (Bandillo et al., 2013; Meng et al., 2016, 2017), wheat (Mackay et al., 2014), and barley (Sannemann et al., 2015). Thus, the use of MAGIC population for genetic studies in crops has risen in the past years.

In this study, we aimed to identify genes/QTLs underlying CEQ in rice using MAGIC population using 6K Illumina Infinium HD Assay. Association analysis was conducted by employing mixed linear model (MLM) to identify the loci associated with CEQ. The results will be valuable in breeding indica rice to improved CEQ and for further use of MAGIC population to illuminate genetic basis of other grain quality traits.

Materials and Methods

Plant Material

The MAGIC population was developed at IRRI using eight diverse founder lines as described in detail by Meng et al. (2016, 2017). The eight founder lines including SAGC-08 (A), HHZ5-SAL9-Y3-Y1 (B), BP1976B-2-3-7-TB-1-1 (C), PR 33282-B-8-1-1-1-1-1 (D), FFZ1 (E), CT 16658-5-2-2SR-2-3-6MP (F), IR 68 (G), and IR 02A127 (H) were selected based on genetic relationship and relevance to rice breeding in terms of yield, grain quality, and biotic or abiotic stress tolerance (Supplementary Table 1). A total of 532 F₆ lines randomly chosen from the population were used in this study.

Phenotyping

The field trials were conducted at the headquarters of International Rice Research Institute (IRRI), Los Baños, Laguna, Philippines in 2014 wet season. Trials design and field management were described by Meng et al. (2016). One hundred and six lines (One-fifth of 532 F₆ lines) and the eight parental lines were replicated three times. Seeds were sown in the seedling nursery on 26 May, 2014. Twenty-day old seedlings were transplanted into field of six-by-six plots with a total of 36 plants in each plot and with a spacing distance of 20 cm × 20 cm. An augmented randomized complete block design was used to layout the trials. A total of 16 plants from middle rows were harvested. Finally, 508 lines yielded enough paddies were selected for further grain quality analysis. Freshly harvested paddy was dried to moisture content of 12–14% and equilibrated in paper bags at room temperature for 2 months. Physicochemical traits of parents and MAGIC population were evaluated at Grain Quality and Nutrition Services Laboratory, IRRI. AC was measured by standard iodine colorimetry method described in ISO 6647-2-2011. GC was measured based on the method described by Cagampang et al. (1973). Pasting viscosity (PV) parameters, the peak viscosity (PKV), trough viscosity (TV), breakdown viscosity (BDV), final viscosity (FV), setback viscosity (SBV), peak time (PKT), pasting temperature (PT), and retrogradation (RT) were measured by RVA to evaluate rheological properties of starch structure (RVA 4500, Perten, Sweden) based on the manufacturer’s instructions. AC, GC, and PV parameters were measured in three replicates.

DNA Extraction and SNP Genotyping

High-throughput DNA extraction technique was carried out at the IRRI Genotyping Service Laboratory using sbeadex^® mini plant kit. High quality DNA samples were quantified with Quant-iT^TM PicoGreen^® dsDNA Assay Kit. SNP Genotyping was performed using 6K Illumina Infinium HD Assay (Illumina Inc., San Diego, CA, United States) following manufacturer’s protocol. Automated calling of SNP array data was done using Illumina GenomeStudio software. A total of 4695 SNP markers yielded good quality genotypic calls. A stringent filtering strategy was further conducted to choose high quality SNPs for association and linkage mapping. All heterozygous markers were set to missing. Markers with minor allele frequency (<0.03) and with missing values (≥0.1) were removed. Highly correlated markers (r² > 0.95) were also excluded from the SNP data set using custom R scripts.

Statistical Analysis

Adjusted trait values for each RIL were obtained using the PBTools software developed by IRRI¹. The broad sense heritability ($H_{\text{B}}^2$) of each trait was calculated using the following equation, $H_{\text{B}}^2$ = σ_G/(σ_G + σ_E/1.2), where σ_G,σ_E and 1.2 represent estimate genotypic variance, experimental error variances, and the number of replication, respectively. Pairwise correlation coefficients for all the traits were estimated using the STAR software (see footnote 1).

Analysis of population structure, LD, and MTA mapping were carried out as described by Meng et al. (2016, 2017). MLM with population structure and familial relatedness between genotypes was conducted using TASSEL V5.0 (Bradbury et al., 2007). Significant MTAs were identified based on their significant association threshold (positive false discovery rate, q^FDR< 0.05) (Benjamini and Hochberg, 1995). The q-values were calculated using the QVALUE in R (Storey and Tibshirani, 2003). Peaks exhibiting the significance threshold level within a physical distance of 1.25 Mb were considered as a single QTL (Meng et al., 2016, 2017). R² was represented explained phenotypic variation using TASSEL V5.0. The plots were visualized using qqman package in R.

Results

Phenotypic Variation

PR 33282-B-8-1-1-1-1-1 had the lowest AC (0.8%) among the founder lines, therefore considered as glutinous parent (Table 1). The rest of the founder lines were non-waxy. SAGC-08, FFZ1, and IR 02A127 were under low AC class (10–20%) while CT 16658-5-2-2SR-2-3-6MP, BP1976B-2-3-7-TB-1-1, HHZ 5-SAL 9-Y 3-Y 1, and IR 68 were under intermediate AC class (20–25%). There were 232, 235, and 12 MAGIC lines with low, intermediate, and high AC classes, respectively (Supplementary Table 2). In addition, 29 lines were glutinous. The average AC of MAGIC population was 17.9%.

TABLE 1

TABLE 1. Characteristics and broad sense heritability for cooking and eating quality traits of parents and MAGIC population.

The GC of PR 33282-B-8-1-1-1-1-1 was the highest (Table 1) among founder lines suggesting its soft textural property. Similarly, BP1976B-2-3-7-TB-1-1 had the soft textural property with a mean value of (61.5 mm). HHZ 5-SAL 9-Y 3-Y 1, SAGC-08, FFZ1, CT 16658-5-2-2SR-2-3-6MP, and IR 02A127 were under medium GC class, therefore had flaky textural property. IR 68 had a hard GC class suggesting very flaky textural property. A total of 190, 145, and 173 MAGIC lines with high, intermediate, and low GC, respectively (Supplementary Table 2). The average GC of the whole population was 50.2 mm.

The RVA was used to mimic the cooking process in rice. In the process, viscosity increases reaching its maximum in the first phase called the PKV [6, 29]. This phase is characterized by swelling and paste formation as the rice starch is being heated up to 65°C with water. PKV of PR 33282-B-8-1-1-1-1-1 was the lowest among founder lines (Table 1), since this line was glutinous with glossy sticky firm and soft textural property. The other founder lines had almost the same PKV values and were much higher than PR 33282-B-8-1-1-1-1-1. This suggests rice with soft textural property will tend to swell and burst faster than those with hard and very flaky textural property. Average PKV of the MAGIC population was 2910. PKT is the time required to reach PKV hence low PKV is indicative of low PKT. The founder line PR 33282-B-8-1-1-1-1-1, with the lowest PKV, had also the lowest PKT values (Table 1). The PKT values of MAGIC population ranged from 3.7 to 7 min with a mean of 6.1 min. Similar with PKT, TV, FV, SBV, PT, and RT of PR 33282-B-8-1-1-1-1-1 are the lowest among founder lines (Table 1). On the contrary, except for BDV, the remains of other seven RVA properties of IR 68 are the highest among eight parents (Table 1).

Broad sense heritability was estimated for all the traits (Table 1). Results showed that all the physicochemical traits had high heritability estimates ranging from 0.81 to 0.95. Notably, AC had the highest heritability and therefore the trait could be easily passed onto the progeny. The high heritability estimates for the other traits were not surprising since AC mainly influences the overall CEQ traits.

Trait Correlations

Based on Pearson’s correlation analysis between physicochemical parameters, almost all the traits were significantly correlated (Table 2). AC was significantly positively correlated with PKT, TV, FV, SBV, PKT, and RT but significantly negative correlated with GC and BDV, which indicate that disruption of starch molecules become harder as AC increase. These findings supported the fact that indeed AC is the most important factor affecting cooking and eating properties. In addition, the correlation coefficient between BDV and other traits analyzed was strongly significant as well. BDV had positive correlation with GC, PKV, and PT, while it was negatively correlated with other traits. It is also noteworthy that PT only shows significant correlation with three traits (AC, PKV, and BDV) with relatively low coefficient.

TABLE 2

TABLE 2. Pairwise correlation values among ten physicochemical properties of MAGIC population.

QTLs Identified by Association Analysis

A total of 62 significant MTAs were identified for the ten physicochemical traits related to CEQ. The Manhattan plots for all the traits were presented in Figure 1. The details of MTAs were given in Supplementary Table 3. A total of 17 QTLs were identified by delineating significant MTAs within a physical distance of 2.5 Mb into a single QTL (Table 3).

FIGURE 1

FIGURE 1. Association analysis for (A) AC, (B) GC, (C) PKT, (D) PKV, (E) TV, (F) BDV, (G) FV, (H) SBV, (I) PT, and (J) RT. X-axis and Y-axis represent chromosome and negative log (q-value), respectively. Horizontal solid line indicates the significant threshold P < 0.05.

TABLE 3

TABLE 3. Putative QTLs for cooking and eating quality.

A total of seven significant MTAs were identified for AC, all of which were located on chromosome 6 (Figure 1A and Supplementary Table 3). These MTAs were located on the interval 1.38–2.13 Mb and with peak SNP at an A/G polymorphism. The allele G was associated with decreased AC, accounting 13.49% of total phenotypic variance (Table 3). The SNP marker, fd8, an A/C polymorphism was just located on GBSSI and explained 6.24% of the total phenotypic variance.

Three QTLs qGC5, qGC6.1, and qGC6.2 were identified for GC from a total of nine significant MTAs (Figure 1B and Supplementary Table 3). These were located on chromosomes 5 and 6, and explained 3.14, 5.96, and 5.24% of the phenotypic variation, respectively (Table 3). The peaked SNPs for qGC5 and qGC6.1 were both A/G polymorphisms located at 8.38 and 2.13 Mb on chromosome 5 and 6, respectively, with the allele G increasing GC. The allele C in C/T polymorphism of qGC6.2 increases GC. An A/C polymorphism was on GBSSI gene explained 3.63% of the total phenotypic variation. This result was consistent with previous study that GBSSI was a major gene affecting GC.

Two QTLs were identified for PKT from a total of three significant MTAs (Table 3, Figure 1C and Supplementary Table 3). The qPKT6.1 and qPKT6.2 were located on chromosome 6 within 1.38 to 1.50 Mb and at 6.99 Mb, explained 4.03 and 4.95% of the phenotypic variations, respectively. The qPKT6.1 peaked at A/G polymorphism whereas qPKT6.2 peaked at C/T polymorphism (Table 3).

One QTL (qPKV6) was identified for PKV (Figure 1D and Supplementary Table 3). This QTL was located on chromosome 6 at 1.77 Mb and peaked at A/C polymorphism (Table 3). The allele C was associated with decreased PKV, explained about 8.5% of phenotypic variation. Moreover, the peaked SNP corresponded with GBSSI gene.

A total of six significant MTAs were identified for TV on chromosome 6, which were delineated into one QTL, qTV6 (Figure 1E and Supplementary Table 3). The qTV6 was located within 1.38 to 2.87 Mb and peaked at A/G polymorphism with allele G being associated with decreased TV, explained 7.69% of the phenotypic variation (Table 3).

Two QTLs, qBDV6.1 and qBDV6.2, both located on chromosome 6 were identified for BDV from a total of eight significant MTAs (Figure 1F and Supplementary Table 3). The qBDV6.1 located within 1.38 to 2.15 Mb peaked at A/G polymorphism (Table 3), which co-localized with the GBSSI gene. The qBDV6.2 peaked at C/T polymorphisms located at 6.99 Mb. The alleles G on qBDV6.1 and allele C on qBDV6.2 were associated with decreased BDV, accounting for 6.01 and 4.78% of phenotypic variation, respectively.

Six significant MTAs were identified for FV (Figure 1G and Supplementary Table 3). These MTAs were delineated into a single QTL, the qFV6 (Table 3). The QTL was located in the interval 1.38 to 2.15 Mb on chromosome 6. Moreover, this peaked at A/G polymorphism and explained 10.34% of the total phenotypic variance. The allele G in this SNP was associated with decreased FV.

A total of two QTLs qSBV6.1 and qSBV6.2 were identified for SBV from eight significant MTAs (Figure 1H and Supplementary Table 3). The qSBV6.1 located within 1.38 to 2.13 Mb on chromosome 6 peaked at A/G polymorphism. The qSBV6.2 on the same chromosome peaked at C/T polymorphism and was located on the same chromosome at 6.99 Mb (Table 3). The alleles G and C of the two peaked SNPs were associated with decreased SBV.

Eight significant MTAs were identified for PT and delineated into three QTILs, the qPT6.1, qPT6.2, and qPT6.3 (Figure 1I and Supplementary Table 3). The qPT6.1 located within 5.71 to 6.40 Mb peaked at the A/G polymorphism and with allele A associated with increased PT, explained 4.33% of the total phenotypic variance. The qPT6.2 located within 6.99 to 7.43 Mb peaked at the C/T polymorphism and with allele C associated with decreased PT, explained 17.18% of the total phenotypic variance. The qPT6.3 peaked at A/G polymorphism, and explained 3.14% of the phenotypic variance (Table 3).

Six significant MTAs were identified for RT and delineated into a single QTL qRT6 (Figure 1J and Supplementary Table 3). The QTL was located on chromosome 6 within 1.38 to 2.13 Mb, accounting for 6.85% phenotypic variation. The peaked SNP, A/G polymorphism, in the region was located at 1.38 Mb and with allele G associated with decreased RT.

Discussion

Pleiotropic Effects and Significant Correlations Among Traits

Pleiotropy is common phenomena in plant genetics. In this study, one QTL locus was hosted or closely to nine QTLs, qAC6, qGC6.1, qPKT6.1, qPKV6, qBDV6.1, qTV6.1, qFV6, qSBV6, and qRT6, suggesting that this locus impact mostly characteristics of CEQ in rice (Figure 1 and Table 3). Specifically, the SNP marker fd8, an A/C polymorphism, located at a physical distance of 1.77 Mb corresponded to GBSSI gene. Another QTL locus, closed to SSIIa, anchored at 6.99 Mb affecting five traits, GC, PKT, BDV, SBV, and PT. These results were coincident with other studies (Lanceras et al., 2000; Wan et al., 2004; Ebadi et al., 2013; Bao, 2014), indicating pleiotropic effects are very common traits related to rice quality. Pleiotropic effects are important in crop improvement especially when QTL/gene affects multiple unrelated traits (i.e., desirable and undesirable). In addition, AC was negatively correlated with GC and BDV (Table 2). Positive correlation of AC with all the RVA parameters except BDV indicates that high AC increase disruption of starch molecules resulting to harder cooked rice gel. Negative correlation of BDV with all RVA parameters especially with SBV indicates that starch remains undisrupted when BDV is low. Above results suggest glossy sticky and soft rice texture tend to undergo fast breakdown of starch molecules (Yan et al., 2005; Wang et al., 2007). These findings supported the fact that indeed AC is the most important factor affecting CEQ (Dela Cruz and Khush, 2000; Fitzgerald et al., 2009; Ni et al., 2011). In the present study, the QTLs exhibiting pleiotropy could be directly utilized for improving CEQ traits in indica rice, since the traits it affects were highly correlated.

QTLs With Small Effects

In the present study, the best MTAs revealed only modest R² value suggesting low variance predicted by each SNP marker. The qPT6.2 revealed only 17.18% of the phenotypic variance whereas the same QTL explained 88% variation in QTL mapping using RILs derived from indica-japonica cross (Wang et al., 2007). Moreover, candidate gene association analysis showed higher phenotypic variance explained for AC (Kharabian-Masouleh et al., 2012; Yang et al., 2014; Zhao et al., 2015). Many association analyses in plants reported low to modest R² value. Bandillo et al. (2013) and Qiu et al. (2015) reported R² values ranging from 11 to 20% and 6.6 to 22.2% for the several quantitative traits in rice, respectively. These corresponded with our results where R² values ranged from 2.3 to 17.2% (Table 3). Possible reasons for these includes insufficient marker coverage used in the study where the causal variant is not in perfect LD with the genotype thereby reducing the association power and variation explained by each SNP marker. This problem could be circumvented with the use of populations showing greater levels of LD (Sorkheh et al., 2008) and more marker coverage. Covering the genome with high marker density is now possible with increasing availability and great reduction in the cost of high throughput SNP genotyping platforms such as genotype-by-sequencing (GBS) (Elshire et al., 2011), Rice 44K (Zhao et al., 2011), and Rice 700K platforms (McCouch et al., 2016). The population size is another consideration to increase the variation explained by the SNP marker. In our study, the population size (508 lines) used is enough to capture high mapping power and resolution since it was suggested that at least 50 individuals per founder lines should be the size of MAGIC population (Cavanagh et al., 2008). However, it was not able to significantly increase the variation explained by best MTAs due to low marker coverage used in the study.

Potential of MAGIC Population for Gene Discovery and Breeding

The use of MAGIC population for gene discovery and breeding are becoming more attractive for geneticists and breeders. Genetic diversity is paramount to crop improvement. Notably, MAGIC population covers high level of recombination which allows for wide allelic and phenotypic diversity and increased precision of QTL mapping (Cavanagh et al., 2008). More recombination events allow for the development of alien diversity resource for use in conventional breeding and conduit for genetic analysis studies. Moreover high throughput SNP genotyping platforms are becoming cheap and statistical approaches for such populations are now available. The MAGIC population used in this study showed wider phenotypic variation in almost all traits than the eight founder lines (Table 1), which suggested the formation of transgressive segregants on both directions. This offers breeders a wide variability to select superior lines to further develop varieties with improved CEQ. The MAGIC lines were presented in a matrix based on AC and GC as these two parameters were the most important factors determining CEQ (Supplementary Table 2). It is interesting that the population can be classified into ten subsets. Compared with eight parents, only four groups, the low AC-soft GC, low AC-hard GC, high AC-medium GC, and high AC-hard GC, were created during three times of intermating and several times of selfing. It is indicating that it is possible to select specific recombinant of AC and GC for different CEQ using MAGIC population. Currently, the MAGIC population we used in this study was part of population improvement program via recurrent selection in IRRI to develop elite varieties of good performance. Interestingly, MAGIC lines were also included in multi-environment yield trials in South Asia, South-east Asia, and Africa. A total of 21 high yielding lines with satisfactory agronomic performance were selected by IRRI breeders for further testing together with 700 lines developed from other IRRI breeding programs. Continuous evaluation of MAGIC lines for biotic and abiotic stresses would allow the development of outstanding lines with high yield and satisfactory agronomic performance under stress conditions. At the present situation, some of the MAGIC lines used in this study were part of salinity tolerance and nutrient deficiency experiments in IRRI.

Conclusion

Overall our results provided details on the chances of QTL detection via association mapping in indica rice using MAGIC population. A total of 17 QTLs were detected for all the physicochemical traits using Infinium 6K genotyping platform, of which 9 and 5 QTLs corresponded to GBSSI and SSIIa genes, respectively. This supports previous studies that GBSSI and SSIIa genes mainly influence the overall CEQ. The QTLs identified revealed small to modest phenotypic variance ranging from 2.89 to 17.18%. These QTLs and the markers highly associated with the traits under study may be useful for breeding indica rice to improve CEQ. Moreover, these findings provided a direct support on the use of MAGIC population as valuable resource for detecting QTL via association mapping and for practical breeding.

Author Contributions

XZ and GY designed the experiment. KP and XZ performed all the phenotypic evaluation, data analysis, and drafted the manuscript. All authors revised the paper and approved the final version to be published.

Funding

This research was financially supported by Key Research and Development Plan of Zhejiang Province (2018C02055), Discipline Construction Plan of ZAAS to XZ’s team and the CAAS Innovative Team Award to GY’s Team.

Conflict of Interest Statement

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.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2018.00868/full#supplementary-material

Footnotes

1. ^ bbi.irri.org

References

Allahgholipour, M., Ali, A. J., Alinia, F., Nagamine, T., and Kojima, Y. (2006). Relationship between rice grain amylose and pasting properties for breeding better quality rice varieties. Plant Breed. 125, 357–362. doi: 10.1111/j.1439-0523.2006.01252.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Bandillo, N., Raghavan, C., Muyco, P. A., Sevilla, M. A. L., Lobina, I. T., Dilla-Ermita, C. J., et al. (2013). Multi-parent advanced generation inter-cross (MAGIC) populations in rice: progress and potential for genetics research and breeding. Rice 6:E11. doi: 10.1186/1939-8433-6-11

PubMed Abstract | CrossRef Full Text | Google Scholar

Bao, J. S. (2014). “Genes and QTLs for rice grain quality improvement,” in Rice - Germplasm Genetics and Improvement, ed. W. Yan (Rijeka: InTech), 239–278.

Google Scholar

Bao, J. S., Corke, H., He, P., and Zhu, L. H. (2003). Analysis of quantitative trait loci for starch properties of rice based on an RIL population. Act. Bot. Sin. 45, 986–994.

Google Scholar

Bao, J. S., and Xia, Y. W. (1999). Genetic control of paste viscosity characteristics in indica rice (Oryza sativa L.). Theor. Appl. Genet. 98, 1120–1124. doi: 10.1007/s001220051175

CrossRef Full Text | Google Scholar

Bao, J. S., Zheng, X. W., Xia, Y. W., He, P., Shu, Q. Y., Lu, Y., et al. (2000). QTL mapping for the paste viscosity in rice (Oryza sativa L.). Theor. Appl. Genet. 100, 280–284. doi: 10.1007/s001220050037

CrossRef Full Text

Benjamini, Y., and Hochberg, Y. (1995). Controlling the false discovery rate: a practical and powerful approach to multiple testing. J. R. Stat. Soc. 57, 289–300.

Google Scholar

Bradbury, P. J., Zhang, Z., Kroon, D. E., Casstevens, T. M., Ramdoss, Y., and Buckler, E. S. (2007). TASSEL: software for association mapping of complex traits in diverse samples. Bioinformatics 23, 2633–2635. doi: 10.1093/bioinformatics/btm308

PubMed Abstract | CrossRef Full Text | Google Scholar

Cagampang, B. G., Perez, C. M., and Juliano, B. O. (1973). A gel consistency test for eating quality of rice. J. Sci. Food. Agric. 24, 1589–1594. doi: 10.1002/jsfa.2740241214

PubMed Abstract | CrossRef Full Text | Google Scholar

Cavanagh, C., Morell, M., Mackay, I., and Powell, M. (2008). From mutations to MAGIC: resources for gene discovery validation and delivery in crop plants. Curr. Opin. Biol. 11, 215–221. doi: 10.1016/j.pbi.2008.01.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Cuevas, R. P., and Fitzgerald, M. A. (2012). “Genetic diversity of rice grain quality,” in Genetic Diversity in Plants, ed. M. Caliskan (Rijeka: InTech), 286–310.

Google Scholar

Darvasi, A., and Soller, M. (1995). Advanced intercross lines an experimental population for fine genetic mapping. Genetics 141, 1199–1207.

Google Scholar

Dela Cruz, N., and Khush, G. S. (2000). “Rice Grain Quality Evaluation Procedures,” in Aromatic Rices, eds R. K. Singh, U. S. Singh, and G. S. Khush (New Delhi: Mohan Primlani for Oxford & IBH Publishing Co), 16–28.

Google Scholar

Ebadi, A., Farshadfar, E., and Rabiei, B. (2013). Mapping QTLs controlling cooking and eating quality indicators of Iranian rice using RILs across three years. Aust. J. Crop Sci. 7, 1494–1502.

Google Scholar

Elshire, R. J., Glaubitz, J. C., Sun, Q., Poland, J. A., Kawamoto, K., Buckler, E. S., et al. (2011). A robust simple genotyping-by-sequencing (GBS) approach for high diversity species. PLoS One 6:e19379. doi: 10.1371/journal.pone.0019379

PubMed Abstract | CrossRef Full Text | Google Scholar

Fitzgerald, M. A., McCouch, S. R., and Hall, R. D. (2009). Not just a grain of rice: the quest of quality. Trends Plant Sci. 14, 133–139. doi: 10.1016/j.tplants.2008.12.004

PubMed Abstract | CrossRef Full Text | Google Scholar

He, P., Li, S. G., Qian, Q., Ma, Y. Q., Li, J. Z., Wang, W. M., et al. (1999). Genetic analysis of rice grain quality. Theor. Appl. Genet. 98, 502–508. doi: 10.1007/s001220051098

CrossRef Full Text | Google Scholar

Holland, J. (2007). Genetic architecture of complex traits in plants. Curr. Opin. Plant Biol. 10, 156–161. doi: 10.1016/j.pbi.2007.01.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, B. E., Verbyla, K. L., Verbyla, A. P., Raghavan, C., Singh, V. K., Gaur, P., et al. (2015). MAGIC populations in crops: current status and future prospects. Theor. Appl. Genet. 128, 999–1017. doi: 10.1007/s00122-015-2506-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, X., Zhao, Y., Wei, X., Li, C., Wang, A., Zhao, Q., et al. (2012). Genome-wide association study of flowering time and grain yield traits in a worldwide collection of rice germplasm. Nat. Genet. 44, 32–39. doi: 10.1038/ng.1018

PubMed Abstract | CrossRef Full Text | Google Scholar

Kharabian-Masouleh, A., Waters, D. L. E., Reinke, R. F., Ward, R., and Henry, R. J. (2012). SNP in starch biosynthesis genes associated with nutritional and functional properties of rice. Sci. Rep. 2:E557. doi: 10.1038/srep00557

PubMed Abstract | CrossRef Full Text | Google Scholar

Lanceras, J. C., Huang, Z. L., Naivikul, O., Vanavichit, A., Ruanjaichon, V., and Tragoonrung, S. (2000). Mapping of genes for cooking and eating quality in Thai jasmine rice (KDML105). DNA Res. 7, 93–101. doi: 10.1093/dnares/7.2.93

PubMed Abstract | CrossRef Full Text | Google Scholar

Mackay, I., Bansept-Basler, J. P., Barber, T., Bentley, A. R., Cockram, J., Gosman, N., et al. (2014). An eight-parent multiparent advanced generation inter-cross population for winter-sown wheat: creation properties and validation. G3 4, 1603–1610. doi: 10.1534/g3.114.012963

PubMed Abstract | CrossRef Full Text | Google Scholar

Mackay, I., and Powell, W. (2007). Methods for linkage disequilibrium mapping in crops. Trends Plant Sci. 12, 1360–1385. doi: 10.1016/j.tplants.2006.12.001

PubMed Abstract | CrossRef Full Text | Google Scholar

McCouch, S., Wright, M., Tung, C.-W., Maron, L., McNally, K., and Fitzgerald, M. (2016). Open access resources for genome wide association mapping in rice. Nature Comm. 7:10532. doi: 10.1038/ncomms10532

PubMed Abstract | CrossRef Full Text | Google Scholar

Meng, L., Wang, B., Zhao, X., Ponce, K., Qian, Q., and Ye, G. (2017). Association mapping of ferrous, zinc, and aluminum tolerance at the seedling stage in indica rice using MAGIC populations. Front. Plant Sci. 8:1822. doi: 10.3389/fpls.2017.01822

PubMed Abstract | CrossRef Full Text | Google Scholar

Meng, L., Zhao, X., Ponce, K., Ye, G., and Leung, H. (2016). QTL mapping for agronomic traits using multi-parent advanced generation inter-cross (MAGIC) populations derived from diverse elite indica rice lines. Field Crops Res. 189, 19–42. doi: 10.1016/j.fcr.2016.02.004

CrossRef Full Text | Google Scholar

Ni, D., Zhang, S., Cheng, S., Xu, Y., Li, L., Li, H., et al. (2011). Improving cooking and eating quality of Xieyou57 an elite indica hybrid rice by marker-assisted selection of the Wx locus. Euphytica 175, 355–362. doi: 10.1007/s10681-011-0377-2

CrossRef Full Text | Google Scholar

Qiu, X., Pang, Y., Yuan, Z., Xing, D., Xu, J., Dingkuhn, M., et al. (2015). Genome-wide association study of grain appearance and milling quality in a worldwide collection of indica rice germplasm. PLoS One 10:E0145577. doi: 10.1371/journal.pone.0145577

PubMed Abstract | CrossRef Full Text | Google Scholar

Sannemann, W., Huang, B. E., Mathew, B., and Leon, J. (2015). Multi-parent advanced generation inter-cross in barley: high-resolution quantitative trait locus mapping for flowering time as a proof of concept. Mol. Breed. 35, 1–16. doi: 10.1007/s11032-015-0284-7

CrossRef Full Text | Google Scholar

Sorkheh, K., Malysheva-Otto, L. V., Wirthensohn, M. G., Tarkesh-Esfahani, S., and Martínez-Gómez, P. (2008). Linkage disequilibrium, genetic association mapping and gene localization in crop plants. Genet. Mol. Biol. 31, 805–814. doi: 10.1590/S1415-47572008005000005

CrossRef Full Text | Google Scholar

Storey, J. D., and Tibshirani, R. (2003). Statistical significance for genomewide studies. Proc. Natl. Acad. Sci. U.S.A. 100, 9440–9445. doi: 10.1073/pnas.1530509100

PubMed Abstract | CrossRef Full Text | Google Scholar

Tran, N. A., Daygon, V., Resurreccion, A. P., Cuevas, R. P., Corpuz, H. M., and Fitzgerald, M. A. (2011). A single nucleotide polymorphism in the Waxy gene explains a significant component of gel consistency. Theor. Appl. Genet. 123, 519–525. doi: 10.1007/s00122-011-1604-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Wan, X. Y., Wan, J. M., Su, C. C., Wang, C. M., Shen, W. B., Shen, J. M., et al. (2004). QTL detection for eating quality of cooked rice in a population of chromosome segment substitution lines. Theor. Appl. Genet. 110, 71–79. doi: 10.1007/s00122-004-1744-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, L. Q., Liu, W. J., Xu, Y., He, Y. Q., Luo, L. J., Xing, Y. Z., et al. (2007). Genetic basis of 17 traits and viscosity parameters characterizing the eating and cooking quality of rice grain. Theor. Appl. Genet. 115, 463–476. doi: 10.1007/s00122-007-0580-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Yan, C., Li, X., Zhang, R., Sui, J., Liang, G., Shen, X., et al. (2005). Performance and inheritance of rice starch RVA profile characteristics. Rice Sci. 2005, 39–47.

Google Scholar

Yang, F., Chen, Y., Tong, C., Huang, Y., Xu, F., Li, K., et al. (2014). Association mapping of starch physicochemical properties with starch synthesis-related gene markers in nonwaxy rice (Oryza sativa L.). Mol. Breed. 34, 1747–1763. doi: 10.1007/s11032-014-0135-y

CrossRef Full Text | Google Scholar

Zhang, C., Hu, B., Zhu, K., Zhang, H., Leng, Y., Tang, S., et al. (2013). QTL mapping for rice RVA properties using high-throughput re-sequenced chromosome segment substitution lines. Rice Sci. 20, 407–414. doi: 10.1016/S1672-6308(13)60131-6

CrossRef Full Text | Google Scholar

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:E467. doi: 10.1038/ncomms1467

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, X., Zhou, L., Ponce, K., and Ye, G. (2015). The usefulness of known genes/QTLs for grain quality in an indica population of diverse breeding lines tested using association analysis. Rice 8:E29. doi: 10.1186/s12284-015-0064-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, C., Gore, M., Buckler, E. S., and Yu, J. (2008). Status and prospects of association mapping in plants. Plant Genome. 1, 5–20. doi: 10.3835/plantgenome2008.02.0089

CrossRef Full Text | Google Scholar