Skip to main content

QTL mapping for starch paste viscosity of rice (Oryza sativa L.) using chromosome segment substitution lines derived from two sequenced cultivars with the same Wx allele



The eating and cooking qualities (ECQs) of rice (Oryza sativa L.) are key characteristics affecting variety adoption and market value. Starch viscosity profiles tested by a rapid visco analyzer (RVA) offer a direct measure of ECQs and represent the changes in viscosity associated with starch gelatinization. RVA profiles of rice are controlled by a complex genetic system and are also affected by the environment. Although Waxy (Wx) is the major gene controlling amylose content (AC) and ECQs, there are still other unknown genetic factors that affect ECQs.


Quantitative trait loci (QTLs) for starch paste viscosity in rice were analyzed using chromosome segment substitution lines (CSSLs) developed from the two cultivars 9311 and Nipponbare, which have same Wx-b allele. Thus, the effect of the major locus Wx was eliminated and the other locus associated with the RVA profile could be identified. QTLs for seven parameters of the starch RVA profile were tested over four years in Nanjing, China. A total of 310 QTLs were identified (from 1 to 55 QTLs per trait) and 136 QTLs were identified in more than one year. Among them, 6 QTLs were stalely detected in four years and 26 QTLs were detected in at least three years including 13 pleiotropic loci, controlling 2 to 6 RVA properties simultaneously. These stable QTL hotspots were co-located with several known starch synthesis-related genes (SSRGs). Sequence alignments showed that nucleotide and amino acid sequences of most SSRGs were different between the two parents. Finally, we detected stable QTLs associated with multiple starch viscosity traits near Wx itself, supporting the notion that additional QTLs near Wx control multiple characteristic values of starch viscosity.


By eliminating the contribution from the major locus Wx, multiple QTLs associated with the RVA profile of rice were identified, several of which were stably detected over four years. The complexity of the genetic basis of rice starch viscosity traits might be due to their pleiotropic effects and the multiple QTL hot spots. Minor QTLs controlling starch viscosity traits were identified by using the chromosome segment substitution strategy. Allele polymorphism might be the reason that QTLs controlling RVA profile characteristics were detected in some known SSRG regions.


Rice (Oryza sativa L.) is one of the most important cereal crops worldwide, with about half of the world’s population consuming it as a main source of dietary calories [1]. Grain quality is a key factor affecting rice production and market value. Current breeding programs aim to improve rice quality particularly with respect to eating and cooking qualities (ECQs).

As the main chemical component of grains, starch and its fine structure determine rice ECQs [2]. The starch viscosity profile is typically tested with a rapid visco analyzer (RVA) and offers a direct measure of ECQs by characterizing the pasting behavior and measuring the changes in viscosity associated with starch gelatinization during cooking. In recent years, starch RVA profiles have become an increasingly popular measure of rice quality, as they are easy to perform and only small samples are required [3]. Indeed, starch RVA profiles have been used as a novel tool to estimate ECQs, and some models were developed successfully in breeding programs [2, 4, 5].

Starch RVA profiles of rice grains are controlled by a complex genetic system that involves multiple quantitative trait loci (QTLs), several with major effects on the trait, and many starch synthesis-related genes (SSRGs) [6]. Linkage mapping using different populations identified a major QTL in the genomic region that contains the gene Waxy (Wx) on chromosome 6. Wx encodes the granule-bound starch synthase (GBSS) that is mainly responsible for longer amylose chains [7,8,9]. Many additional loci associated with starch RVA profiles were also shown to affect the ECQs of rice [10]. The alkali degeneration gene locus (ALK) on chromosome 6 encoding starch synthase II a (SSII-3), is the major factor responsible for gelatinization temperature (GT) and the distribution of amylopectin chain length [11,12,13]. Other SSRGs, such as genes encoding debranching enzyme (DBE), isoamylases (ISA), starch branching enzyme (SBE), soluble starch synthase (SSS), and pullulanase (PUL), play minor effects on starch RVA profiles [14]. These SSRGs have widely been proposed to control amylose content (AC), GT, or certain other rice quality traits [15,16,17].

The heredity of starch RVA profile characters is complex and could easily be effected by environmental factors, especially the temperature during the grain filling stage [18, 19]. This brings more difficulty in mapping QTLs of consistent RVA traits. Over 200 starch RVA traits-related QTLs have been identified in various mapping populations [3, 7, 18,19,20,21,22,23,24,25]. A few QTLs have been cloned or finely mapped (, Certain varieties carrying the same Wx allele exhibit different ECQs, indicating that other minor QTLs affect ECQs. Identifying these minor QTLs will be critical to precisely modulate rice grain quality in the future. However, because the effects of minor QTLs might be covered by major QTLs such as Wx, fine mapping and the cloning of the underlying locus can be challenging. Populations derived from two parents with similar AC or with the same Wx allele have been developed to eliminate the effects of Wx and detect minor QTLs for starch RVA and ECQs of rice [23, 24, 26, 27].

The genetic basis of complex traits in rice and the contribution of single-nucleotide polymorphisms have been explored by genome-wide association studies, allowing the generation of testable hypotheses relating allele variation and starch RVA traits [28]. Allelic diversity of SSRGs can explain the genetic basis for the observed phenotypic diversity in starch physicochemical properties across germplasm [14, 17]. However, the distribution of functional alleles is strongly correlated with population structure, which can lead to false results if the population is small [29]. As one of the secondly populations, chromosome segment substitution lines (CSSLs) are very useful for precisely mapping QTLs and dissecting the genetic basis of complex traits [30]. A set of CSSLs had been developed which derived from two sequenced rice cultivars, the indica variety 9311 (recipient) and the japonica variety Nipponbare (donor) [31,32,33]. Since both cultivars carry the same Wx-b allele [22], this CSSLs could be used as an excellent mapping population to detect new QTLs that might be covered by the major Wx locus.

In this study, we aimed to detect stable QTLs in the absence of variation at the Wx locus for parameters of starch RVA profiles using CSSLs. Our results establish the foundation for fine mapping and subsequent cloning of genes responsible for these QTLs, which will help to improve ECQs in rice by molecular breeding.


Performance of starch RVA profiles in parents and CSSLs in different years

The two parents 9311 and Nipponbare showed similar apparent amylose content (AAC), 15.5 % ± 0.2 and 14.6 % ± 0.2 %, respectively, in Nanjing from 2017 to 2018. The AAC of CSSLs ranged from 13.1 to 17.2 %, with an average AAC of 14.92 % in 2017. The starch RVA profiles of 9311 and Nipponbare differed significantly over four years of the study. Most RVA parameters varied greatly over four years, with the exception of peak time (PeT) and pasting temperature (PaT), which were less affected by environment (Table 1). For indica rice 9311, almost all parameters were much smaller than those in Nipponbare, with the exception of setback viscosity (SB) and PeT.

Table 1 Phenotypic variations of starch RNA characteristics among 9311, Nipponbare, and their CSSL populations

The starch RVA profiles showed a continuous distribution in the CSSL population over the different years, with the exception of PaT, which exhibited a double-peak distribution (Fig. 1). The maximum value of peak viscosity (PV) was greater than the other parameters as consequence of greater kurtosis. The starch RVA profiles of CSSLs showed a partial separation for some parameters, such as PV, final viscosity (FV), and PaT, in 2016 and 2017. Among the CSSLs, the phenotypic values of PV, FV, and PaT exhibited a greater range of variation than other traits over four years, whereas breakdown viscosity (BD), SB, and PeT varied little. The mean values for the various RVA parameters of the CSSLs were near the mid-parent value, but some values were not (Table 1). In addition, the phenotypic values for all starch RVA profiles showed bidirectional ultra-parental genetic types in the CSSL population, consistent with polygenic control of those quantitative traits.

Fig. 1
figure 1

Distribution of starch rapid visco analyzer profile characteristics in the Nipponbare/9311 chromosome segment substitution line population in 2014. Note: The mean values for 9311 and Nipponbare are indicated by arrows labeled 9 and N, respectively

Correlation analysis using average values for each year as a variable revealed that PV was positively correlated with trough viscosity (TV), BD, and FV, but negatively correlated with SB and PaT. TV was positively correlated with FV and PeT. FV was positively correlated with PeT and PaT but negatively correlated with SB. BD was negatively correlated with all other RVA parameters. SB was positively correlated with PeT and PaT. These correlations were highly significant (P < 0.01). Over four years, the relationship between SB and TV or FV changed (Supplementary Tables 1, 2, 3). SB was positively correlated with FV except in 2013. TV was positively correlated with SB in 2016 and 2017.

Microclimate analysis

Every year, 9311 and Nipponbare bloomed on July 20th and August 25th, respectively. The heading time of CSSLs ranged from August 15th to the 25th. Although the maximum temperature of July and August was more than 40 °C in 2017, CSSLs and parents showed normal seed sets across the four years. The mean daily maximum temperature during the recorded periods ranged from 28.7 to 34.3 °C over the four years of this study. The daily average temperatures during observational periods in 2013, 2016, and 2017 were about 28.0 °C, higher than the average temperature of 25.1 °C measured in 2014 (Fig. 2). The average daily maximum temperatures showed the same tendency as the daily average temperatures over all four years. The highest recorded temperature of the four years was in 2017. The relative humidity (RH) varied extensively over four years independently of temperature. Average daily RH during the observational periods was higher in 2014 and 2017 at about 91 % lower than those in 2013 and 2016 at 80.5 and 82.5 %.

Fig. 2
figure 2

Daily maximum temperature recorded during the observational periods over the four years of this study

QTL analysis

To elucidate the effects of environmental factors on starch viscosity among CSSLs, we mapped QTLs that influence starch RVA profile characteristics each year. In total, 310 QTLs were detected for seven RVA characteristics. Of them, 44 QTLs were related to PV, 42 to TV, and 43 to FV; 43, 49, 34, and 55 QTLs were identified for BD, SB, PeT, and PaT, respectively.

22 loci for PV, 15 for TV, 25 for FV, and 11 for BD were stably detected across multiple years. 28, 8, and 27 QTLs were mapped in different years for SB, PeT and PaT, respectively. Totally, 81 QTLs were mapped to the same genomic region across two years, and 50 QTLs were identified across three years. Five stable QTLs controlling TV, FV and PaT were identified across all four years (Table 2). The remaining 174 QTLs appeared in only one year. These results suggested that, although most traits showed large variation across the years, several QTLs related to starch RVA profile characters were robust enough to be consistently identified each year. For example, we repeatedly detected a QTL around 16.9–19 Mb of chromosome 2 that affects both FV and TV each year. Similarly, qPaT7 and qTV7 located near marker RM432 on chromosome 7 and qFV9 located near RM219 on chromosome 9 were also mapped every year. A QTL near marker RM3827 on chromosome 6 for SB and a QTL near marker RM1812 on chromosome 11 for BD have not been previously described to be associated with these traits at the genomic locations (Table 3). Six QTLs were not affected by environmental factors in our research, and they are likely to substantially contribute to the RVA characteristics of rice starch.

Table 2 Quantitative trait locus (QTL) analysis for properties of the starch rapid visco analyzer profile
Table 3 Stable rapid visco analyzer (RVA) quantitative trait loci (QTLs) detected more than three times

Several of the QTLs detected here exhibited pleiotropic effects, as they influenced multiple starch RVA profile characteristics. In total, 26 genomic intervals corresponding to 55 QTLs associated with starch RVA profiles were identified in at least three years. Of those, 13 chromosomal intervals showed pleiotropy by affecting more than one starch RVA profile characteristic (Table 3). In particular, the region near marker RM219 on chromosome 9 contained one QTL associated with all starch RVA profile characteristics with the exception of TV. Moreover, the interval between markers RM469 and RM587 on chromosome 6 was associated with five parameters, while the interval near marker RM3795 on chromosome 2, markers RM289 and RM178 on chromosome 5, and marker 10-1.63 on chromosome 10 were each associated with four RVA parameters.

Several reports have described QTLs affecting rice starch RVA profiles [18, 21,22,23, 34]. The 26 intervals detected more than three times were compared with previously mapped QTL. With the exception of two intervals, all other intervals co-located with known QTLs associated with starch RVA profiles from different mapping populations and different environments (Table 3). Specifically, a number of certain hot spots were identified, such as the region between markers RM469 and RM587 on chromosome 6 and the region around marker RM1375 on chromosome 10. In the genomic interval from markers RM469 to RM587 on chromosome 6, we mapped QTLs for SB, PV, BD, PeT, and PaT, as the previous studies did for the traits PeT, TV, FV, BD, and SB [1, 3, 7, 21, 24]. qBD10, located near RM1375 (15.9–18 Mb) on chromosome 10, overlapped with the mapping interval described for certain starch RVA profile QTLs related to FV, SB, PeT, TV, PaT, and PV [3, 24, 34].

We noticed that several of our stable QTL hotspots co-located with some known SSRGs, such as Wx, SSIII-1, SSIV-2, and SBE3 (Table 4). Besides the QTL in the interval 0.56–2.86 Mb on chromosome 6, which co-located with Wx, the region spanning 3.24–5.38 Mb and 30.45–32.65 Mb on chromosome 4, 24.59–26.37 Mb on chromosome 5, and 16.9–20.95 Mb on chromosome 2 contained PUL, SSIII-1, SSIV-2, and SBE3, respectively. The loci ADPlar (Adenosine diphosphate glucose pyrophosphorylase large subunit 1), ADPsma (Adenosine diphosphate glucose pyrophosphorylase small subunit), SSII-1, and ISA were near our starch RVA QTLs. In the case of ADPlar, ADPiso (Adenosine diphosphate glucose pyrophosphorylase large subunit 2), GBSSII, SSII-2, SSIII-2, SSIV-1, and PUL, the nearby or overlapping QTLs were identified in only one or two years, suggesting that the underlying QTLs are easily influenced by the environment. No QTLs associated with starch RVA traits was detected in regions that overlap with the genomic location of SSI, SSII-3, SBE1, or SBE4. The region between markers RM6748 and RM5473 on chromosome 4 contained a QTL for SB and overlapped with the genomic location of SSIII-2, which is related to SB. SSIII-2 is a critical gene controlling Pat and PeT [6]. The effect of this region on Pat and PeT was not tested in our research.

Table 4 Co-location of known SSRGs with the starch rapid visco analyzer (RVA) quantitative trait loci (QTLs) detected in this study

The QTL hotspots near SSII-1 and SSII-3 were previously reported to have major effects on PT and minor effects on gel consistency (GC), AC, PV, CPV, BD, and SB [2]. Here, we found that the QTL near SSII-1 is associated with BD. A QTL hotspot near SBE3 on chromosome 2 was also reported, in which four major QTLs associated with HPV, consistency viscosity, viscosity at 95 °C, and BD over two years [27]. We also detected QTL hotspots that co-located with SBE3 and are associated with FV, TV, SB, and PaT.

Sequence alignment of SSRGs between parents 9311 and Nipponbare

The stable QTL hotspots were co-located with several known SSRGs. To explore the potential sequence variation between the two parental cultivars at these SSRGs, we aligned the nucleotide and predicted protein sequences from 9311 and Nipponbare, which revealed limited allelic diversity between the two parents, with some exceptions (Table 4). For example, the parental alleles for SSI and SSII-3 were identical in sequence in 9311 and Nipponbare, possibly explaining why we failed to detect a QTL overlapping with either of these intervals in this study. Likewise, we detected no QTL near SBE1 and SBE4, whose sequence identity between the two parents was 99.71 and 99.87 %, respectively.

The two parents also shared the same allele at ADPlar and SSIV-2, whose genomic coordinates did not overlap with any of our QTL intervals. However, our QTLs related to RVA profile parameters overlapped with or in close proximity to 12 SSRGs. Although the gene sequences of Wx are identical between 9311 and Nipponbare, we identified QTLs associated with multiple RVA traits in the Wx region. In addition, the nucleotide sequences and encoded protein sequences for ADPsma, SSII-1, SSII-2, SSIII-1, SSIII-2, SSIV-1, SBE3, ISA, and PUL were different between the two parents, suggesting that sequence polymorphisms at these loci contribute to the observed QTLs controlling RVA profile characteristics.


Recent studies have confirmed that changes in the environment largely affect rice starch PV, TV, and FV, while AC and other starch RVA parameters are mainly influenced by genotypes [1, 11]. In this study, starch RVA profiles of both parents and CSSLs varied greatly across four years, especially in 2016 and 2017. One possible reason is the high temperature of the grain filling stage during the summer of these years [38,39,40]. High temperature during flowering will make the spikelet sterile and reduce the seed set. Air temperature during grain filling significantly affects grain composition, as well as starch structure and properties [41, 42]. High temperatures result in lower AC, poor ECQs, and higher pasting properties by decreasing the activity of starch synthesis enzymes. Under such conditions, rice starch contained fewer long-chain amylose and short-chain amylopectin, but more intermediate- and long- chain amylopectin than under normal conditions, as well as higher crystallinity and gelatinization properties [42].

We observed that PV is positively correlated with TV, BD, and FV but negatively correlated with SB and PaT. Similar correlations were previously reported [1, 3, 34]. We also established that PV, FV, TV, BD, PaT, and PeT show consistent cross-correlations regardless of their great variations across the years, indicating that these traits are mainly affected by genotypic variation. There were main effect loci that controlled these starch RVA profile characteristics. The interaction between genotype and environment exhibited strongly for SB, as the correlation coefficients between SB and other parameters (except PV) differed for each year.

The confounding effects of the environment and the major effect gene Wx make it difficult to identify more minor QTLs for starch RVA profiles. Only a few genes and QTLs related to starch RVA profiles have been cloned, such as qAC2, qGC6, ALK, Chalk5, and Dul (Dull) [6]. In this study, we specifically selected two cultivars with the same Wx allele, Wx-b, to generate CSSLs as the mapping population [22]. SSII-3, the major gene controlling GT, showed no differences in its genomic sequence and only a 2-amino acid difference in its protein sequence between 9311 and Nipponbare. We therefore hypothesize that 9311 and Nipponbare carry the same allele at SSII-3. New QTLs that exert minor effects on starch RVA profiles and GT could be identified without the influence of major Wx and SSII-3. Zhang developed a CSSL population consisting of 38 lines with the same parents selected here. Following the same mapping strategy, 10 stable QTLs for RVA properties and two minor stable QTLs for GT were identified over two years and in two environments [22, 23].

The QTLs for the RVA properties they mapped also displayed pleiotropy but didn’t map to the same or similar chromosomal regions with our QTLs, which might not be surprising in considering that the two sets of CSSL lines harbored different chromosome segments from Nipponbare. Liu et al. (2011) also performed an extensive QTL mapping analysis for 16 rice quality traits across eight environments using a set of CSSLs developed from two rice varieties with similar AC, and detected 10 stable RVA profiles cross four environments [19]. However, the number of molecular markers was limited, raising the possibility of missing small introgression segments and thus lowering the accuracy of QTL detection. With more lines and smaller substitution segments, our CSSLs are adequate to identify minor QTLs for starch RVA profiles and GT without the effect of SSII-3 and Wx.

In our research, although 310 QTLs were detected for seven paste viscosity properties for rice starch RVA profile traits, only 136 QTLs were mapped repeatedly. These results indicated that the RVA profile characteristics are largely affected by the environment, consistent with previous reports [18, 19]. The novel loci we detected for all starch RVA parameters, and in particular the six QTLs we identified every year, are important and require further validation. In addition, we propose that new and stable QTLs, such as qSB6 and qBD11, which had not been previously reported, may be useful for research on marker-assisted selection of ECQs.

Based on our mapping, the clustering of multiple QTLs controlling RVA traits indicate that pleiotropic effects and QTL hotspots are key factors affecting starch RVA traits in rice. Similar observations of pleiotropy were reported previously (, [3, 24]). QTLs with high correlations are often grouped in the same or adjacent marker intervals on a chromosome [43]. Starch RVA traits showed significant correlations that confirmed the linkage or pleiotropy of the corresponding loci.

Overlap between the intervals of our mapped QTLs and known SSRGs was common in this research (Table 4). Most known SSRGs showed different alleles between the parents, 9311, and Nipponbare. This result suggested that alleles of SSRGs contributed to our QTLs controlling RVA profile characteristics under the same major gene Wx, which should be taken into consideration in rice quality breeding. The effect of SSRG allele combinations on starch quality and ECQs should be investigated more in the future. We failed to identify QTLs associated with RVA profiles near some SSRGs, such as SSII-3 and SSI, likely due to the lack of allele polymorphism between the two parents. The 5ʹ untranslated region (UTR), genomic and protein sequences of Wx are identical between 9311 and Nipponbare, although we noted several RVA traits associated with the Wx genomic region. Xu also reported several SNPs close to Wx that were significantly associated with RVA parameters in subpopulations with the same Wx allele [15]. These loci will be high-priority candidates for future characterization.

The highly consistent genotypes between SSRGs and QTLs reported here confirmed the accuracy of our mapping. The stable or newly developed QTLs located in the region that no known SSRGs in are meaningful to identify new genes controlling starch RVA profiles in the future, which will facilitate further research into the genetic mechanism regulating RVA profiles and ECQ of rice.


We mapped QTLs associated with starch viscosity profile, one of the most important factors contributing to ECQs in rice. The effect of major locus Wx could be eliminated by generating CSSLs developed from 9311 and Nipponbare with the same Wx-b allele, which allowed the identification of new QTLs associated with RVA profile traits. We analyzed the genetic basis of variation in RVA profile over four years, and identified 136 repeated QTLs. Among them, 6 stable QTLs were detected every year, which are therefore likely to be very important for the RVA characteristics. In addition, we uncovered 13 intervals detected for 3 of the four years that showed pleiotropy with respect to controlling two to six starch RVA profile properties simultaneously. Finally, we highlighted four intervals, such as the interval between markers RM469 and RM587 on chromosome 6, that are associated with more than four RVA parameters for a given year. Pleiotropic effects and QTL hotspots appear to be key factors affecting starch RVA traits in rice. There was high consistency between allelic diversity at known SSRGs and the QTLs reported here. Future research will explore these stable QTLs and hotspots in more detail.


Plant materials and field planting

An advanced backcross population was developed by our lab using the indica variety ‘9311’ (recipient) and the japonica variety ‘Nipponbare’ (donor). 9311 was obtained from its breeder, Yangzhou Institute of Agricultural Sciences, Jiangsu Academy of Agricultural Sciences. Nipponbare was obtained from Jiangsu Provincial Platform for Conservation and Utilization of Agricultural Germplasm.

The backcrossed population consisted of 119 BC4F2 lines. Backcrossing and simple sequence repeat marker selection were performed as described in detail by Zhu and Zhao [20, 31]. Each introgression line in the population contained one to seven segments originating from Nipponbare. The 119 lines contained 318 substituted segments with an average of two to seven segments per line and covered 84.0 % of the whole rice genome [31].

The CSSLs and two parents were planted in fields at the Jiangsu Academy of Agricultural Sciences (32°02’N, 118°52’E; elev. 10 m) in Nanjing in 2013, 2014, 2016, and 2017. Each year, all seeds were planted within two blocks on May 15th and transplanted on June 20th. Each line was planted in three rows with a row-to-row distance of 30 cm and plant-to-plant distance of 13.3 cm. The seeds of five plants were sampled from each CSSL line and dried naturally. Milled rice was grinded into powder and then passed through a 100-mesh sieve. After drying at 4 °C in an oven, each powder sample was balanced for 2 days at room temperature and kept at 4 °C for 3 months. Then, paste viscosity was measured as described below. For each line, five samples were used to determine the starch RVA profile.


The climate parameters at the field site (air temperature, RH) were measured using a data logger (Thermo Recorder TR-72U, T & D Corp, Japan). The sensor was placed at a height of 170 cm; air temperature and RH were collected every 10 min. The climate parameters were collected from the beginning of flowering to full maturity each year.

Starch RVA profile

Starch paste viscosity was measured with a Rapid Visco Analyser (Tecmaster, Perten, Sweden) according to the American Association of Cereal Chemists Standard Method (AACC 61 − 02) with TCW software 3 (Thermal Cycle for Windows) [44]. 3 g rice flour (accounting for 12 % moisture basis) with 25 mL distilled water were used. The heat 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 over 3.75 min; (3) the temperature was held at 95 °C over 2.5 min; (4) the temperature was ramped down linearly to 50 °C for 3.75 min; (5) the temperature was held at 50 °C for 1.4 min. The RVA paddle speed was set to 960 rpm for the first 10 s of the test, after which the speed was 160 rpm.

Starch paste viscosity characteristics are characterized by five parameters : PV, TV, FV, PeT, and PaT. BD and SB were calculated as: BD = PV – TV, and SB = FV – PV [36]. Correlations between the RVA parameters for each year were analyzed by IBM SPSS Statistics v22.

The AAC was determined using the iodine staining method described in the European Standard EN ISO 6647-2-2015. The absorbance of the solution was measured at 620 nm against the blank solution using a continuous flow analyzer (Seal Analytical AA3, GER). The AAC was calculated using a standard curve made from four rice samples with known AAC.

QTL mapping

Genotype data for 250 polymorphic loci, including 211 simple sequence repeat and 39 sequence tag site markers, were used for QTL detection. Molecular linkages were established using composite interval mapping with version 3.3 of QTL IciMapping software [45].

QTL were detected according to the method described by Eshed and Zamir [46]. The significance of each QTL was determined by comparing the mean RVA profile values of a CSSL line with the recipient parent 9311 using analysis of variance and Dunnett’s test. A QTL was considered as present when a CSSL line exhibited a significant difference compared to 9311 with corresponding probability value P < 0.05. If more than three CSSLs showed differences, then the QTL was estimated as being located within the chromosomal region shared by those CSSLs [22]. QTL nomenclature followed as that of McCouch et al. [47].

Sequence alignment of SSRGs between parents 9311 and Nipponbare

The genomic and predicted protein sequences of known SSRGs that co-locate with starch RVA QTLs were downloaded from online resources, as 9311 and Nipponbare genomes are sequenced. Nipponbare sequences were downloaded from Gramene ( The genomic sequences for 9311 were obtained from the Rice Genome Project (Beijing Genomics Institute) [32]. The predicted protein sequences were obtained from Gramene or NCBI ( Sequence alignment was performed with DNAMAN (version 6).

Availability of data and materials

All materials analyzed during the current study are available from the Jiangsu Provincial Platform for Conservation and Utilization of Agricultural Germplasm ( The stock numbers for 9311 and Nipponbare are M3A00407620 and SD_NJAU_B149, respectively. The results of correlation analysis among RVA parameters in different years are in the Additional file (1) The stock numbers, raw phenotype data and genotype data for all CSSLs are in Additional file (2) The details about the sequence alignments are in the additional file 3.



Eating and cooking qualities


Rapid Visco Analyzer


Amylose content


Quantitative trait loci


Chromosome segment substitution lines


Starch synthesis related gene


Granule-bound starch synthase


Alkali degeneration


Starch synthase


Gelatinization temperature


Debranching enzyme




Starch branching enzyme


Soluble starch synthase




Apparent amylose content


Pasting temperature


Peak time


Setback viscosity


Peak viscosity


Final viscosity


Breakdown viscosity


Trough viscosity


Relative humidity


Adenosine diphosphate glucose pyrophosphorylase large subunit 1


Adenosine diphosphate glucose pyrophosphorylase small subunit


Adenosine diphosphate glucose pyrophosphorylase large subunit 2


Gel consistency




Untranslated region


  1. Tong C, Chen Y, Tang F, Xu F, Hung Y, Cheng H, et al. Genetic diversity of amylose content and RVA pasting parameters in 20 rice accessions grown in Hainan, China. Food Chem. 2014;161:239–45.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  2. Cozzolino D. The use of the rapid visco analyser (RVA) in breeding and selection of cereals. J Cereal Sci. 2016;70:282–90.

    CAS  Article  Google Scholar 

  3. Yang Y, Ni D, Song F, Li L, Lu X, Li Z, et al. Identification of QTLs for rice starch RVA profile properties in different ecological sites. Acta Agronomica Sinica. 2012;38(2):264–74 (in Chinese with English abstract).

    CAS  Article  Google Scholar 

  4. Vk BR, Ns C. Deploying viscosity and starch polymer properties to predict cooking and eating quality models: a novel breeding tool to predict texture. Carbohyd Polym. 2021;260:117766.

    Article  CAS  Google Scholar 

  5. Mao T, Zhang Z, Ni S, Zhao Y, Li X, Zhang L, et al. Assisted selection of eating quality progeny of indica (O. sativa L. ssp. indica) and japonica (O. sativa L. ssp. japonica) hybrids using rice starch properties. Genet Resour Crop Evol. 2021;68:411–20.

    CAS  Article  Google Scholar 

  6. Yan C, Tian Z, Fang Y, Yang Y, Li J, Zeng S, et al. Genetic analysis of starch paste viscosity parameters in glutinous rice Oryza sativa L. Theor Appl Genet. 2011;122:63–76.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  7. Bao J, Zheng X, Xia Y, He P, Shu Q, Lu X, et al. QTL mapping for the paste viscosity characteristics in rice (Oryza sativa L.). Theor Appl Genet. 2000;100:280–4.

    CAS  Article  Google Scholar 

  8. Lapitan V, Redona E, Abe T, Brar D. Mapping of quantitative trait loci using a double–haploid population from the cross of indica and japonica cultivars of rice. Crop Sci. 2009;49:1620–8.

    CAS  Article  Google Scholar 

  9. Zhang C, Chen S, Ren X, Lu Y, Liu D, Cai X, et al. Molecular structure and physicochemical properties of starches from rice with different amylose contents resulting from modification of OsGBSI activity. J Agr Food Chem. 2017;65(10):2222–32.

    CAS  Article  Google Scholar 

  10. Su Y, Rao Y, Hu S, Yang Y, Gao Z, Zhang G, et al. Map-based cloning proves qGC-6, a major QTL for gel consistency of japonica/indica cross, responds by Waxy in rice (Oryza sativa L.). Theor Appl Genet. 2011;123:859–67.

    PubMed  Article  PubMed Central  Google Scholar 

  11. Bao J, Sun M, Zhu L. Corke H. Analysis of quantitative trait locus for some starch properties in rice (Oryza sativa L.), thermal properties, gel texture, swelling volume. J Cereal Sci. 2004;39:379–85.

    CAS  Article  Google Scholar 

  12. Fan C, Yu X, Xing Y, Xu C, Luo L, Zhang Q. The main effects, epistatic effects and environmental interactions of QTLs on the cooking and eating quality of rice in a doubled–haploid line population. Theor Appl Genet. 2005;110:1445–52.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  13. Umemoto T, Yano M, Satoh H, Shomura A, Nakamura Y. Mapping of a gene responsible for the difference in amylopectin structure between japonica–type and indica–type rice varieties. Theor Appl Genet. 2002;104:1–8.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  14. Tian Z, Quian Q, Liu Q, Yan M, Liu X, Yan C, et al. Allelic diversities in rice starch biosynthesis lead to a diverse array of rice eating and cooking qualities. Proc Natl Acad Sci USA. 2009;106:21760–5.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  15. Xu F, Bao J, He Q, Park Y. Genome–wide association study of eating and cooking qualities in different subpopulations of rice (Oryza sativa L.). BMC Genom. 2016;17:663.

    Article  Google Scholar 

  16. Xiang X, Kang C, Xu S, Yang B. Combined effects of Wx and SSIIa haplotypes on rice starch physicochemical properties. J Sci Food Agri. 2017;97:1229–34.

    CAS  Article  Google Scholar 

  17. Yang F, Chen Y, Tong C, Huang Y, Xu F, Li K, ,et al. Association mapping of starch physicochemical properties with starch synthesis–related gene markers in nonwaxy rice (Oryza sativa L.). Mol Breeding. 2014;34:1747–63.

    CAS  Article  Google Scholar 

  18. Zhang Q, Zhan Y, Zhu Z, Zhao L, Zhao Q, Xu L, et al. Inheritance analysis and QTL mapping of rice starch viscosity (Rapid Visco Analyzer profile) characteristics. Rice Sci. 2008;15:186–94.

    CAS  Article  Google Scholar 

  19. Liu X, Wan X, Ma X, Wan J. Dissecting the genetic basis for the effect of rice chalkiness, amylose content, protein content, and rapid viscosity analyzer profile characteristics on the eating quality of cooked rice using the chromosome segment substitution line population across eight environments. Genome. 2010;54:64–80.

    Article  CAS  Google Scholar 

  20. Zhu W, Lin J, Yang D, Zhao L, Zhang Y, Zhu Z, et al. Development of chromosome segment substitution lines derived from backcross between two sequenced rice cultivars, indica recipient 93 – 11 and japonica donor Nipponbare. Plant Mol Biol Rep. 2009;27:126–31.

    CAS  Article  Google Scholar 

  21. Zhang Y, Jiang L, Liu X, Liu S, Chen L, Zhai H, et al. Analysis of QTLs for starch RVA profile properties in the superior rice cultivar Koshihikari. China J Rice Sci. 2010;24(2):137–44. (in Chinese with English abstract).

    Google Scholar 

  22. Zhang C, Bing H, Zhu K, Zhang H, Leng Y, Tang S, et al. QTL mapping for rice RVA properties using high-throughput resequenced chromosome segment substitution lines. Rice Sci. 2013;20(6):407–14.

    Article  Google Scholar 

  23. Liu X, Zhu K, Zhang C, Hong R, Sun P, Tang S, et al. Mapping of minor QTLs for rice gelatinization temperature using chromosome segment substitution lines from indica 9311 in the japonica background. Acta Agronomica Sinica. 2014;40(10):1740–7. (in Chinese with English abstract).

    CAS  Article  Google Scholar 

  24. Yao X, Wang J, Liu J, Zhang J, Ren C, Ma D, et al. Mapping quantitative trait loci associated with starch paste viscosity in rice (Oryza sativa L.) under different environmental conditions. Plant Breed. 2017;136:591–602.

    CAS  Article  Google Scholar 

  25. Shar T, Sheng Z, Ali U, Fiaz S, Wei X, Xie L, et al. Mapping quantitative trait loci associated with starch paste viscosity attributes by using double haploid populations of rice (Oryza sativa L.). J Integr Agric. 2020;19(7):1691–703.

    CAS  Article  Google Scholar 

  26. Bao J, Wu Y, Hu B, Wu P, Cui H, Shu Q. QTL for rice grain quality based on a DH population derived from parents with similar apparent amylose content. Euphytica. 2002;128:317–24.

    CAS  Article  Google Scholar 

  27. Yan B, Tondi Y, Chen J, Wang Y, Gao G, Zhang Q, et al. 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 Breed. 2014;34:2151–63.

    CAS  Article  Google Scholar 

  28. Lu F, Park Y. Sequence variations in OsAGPase significantly associated with amylose content and viscosity properties in rice (Oryza sativa L.). Genet Res. 2012;94:179–89.

    CAS  Article  Google Scholar 

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

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  30. Xu J, Zhao Q, Du P, Xu C, Wang B, Feng Q, et al. Developing high throughput genotyped chromosome segment substitution lines based on population whole–genome re–sequencing in rice (Oryza sativa L.). BMC Genom. 2010;11:656.

    CAS  Article  Google Scholar 

  31. Zhao C, Zhou L, Yu X, Zhao Q, Chen T, Yao S, et al. Identification of tiller angle quantitative trait loci based on chromosome segment substituted lines and high-density physical map in rice. Bull Botany. 2012;47(6):1–8. (in Chinese with English abstract).

    Google Scholar 

  32. Yu J, Hu SN, Wang J, Wong G, Li S, Liu B, et al. A draft sequence of the rice genome (Oryza sativa L. ssp. indica). Science. 2002;296:79–92.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  33. Goff SA, Ricke D, Lan TH, Presting G, Wang R, Dunn M, et al. A draft sequence of the rice genome (Oryza sativa L. ssp. japonica). Science. 2002;296:92–100.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  34. Zhang J, Zheng L, Cai Y, You X, Kong F, Wang G, et al. Correlation analysis and QTL mapping for starch RVA profile properties and amylose and protein contents in rice. China J Rice Sci. 2017;31(1):31–9. (in Chinese with English abstract).

    Google Scholar 

  35. Shao X, Chen Z, Liu Z, Ding C, Tang S, Li G, et al. QTL mapping for rice RVA properties using chromosome segment substitution lines derived from a cross between Koshihikari and 9311. Mol Plant Breed. 2015;2:261–8. (in Chinese with English abstract).

    Google Scholar 

  36. Chen L, Yan X, Wang L, Gao W, Yang J, Chen S, et al. Pasting properties, grain-filling characteristics and allelic variation linked to the grain quality in diverse rice. Euphytica. 2017;213:29–44.

    Article  CAS  Google Scholar 

  37. Hsu Y, Tseng M, Wu Y, Lin M, Wei F, Hwu K, et al. Genetic factors responsible for eating and cooking qualities of rice grains in a recombinant inbred population of an inter-subspecific cross. Mol Breed. 2014;34:655–73.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  38. Mei M, Jiang Y, Wan Z, Liu L, Ye D, Wang Y, et al. National Climate Centre. Climatic characteristics and major meteorological events over China in 2016. Meteorological Monthly. 2017;43(4):468–76. (in Chinese with English abstract).

    Google Scholar 

  39. Feng A, Zeng H, Yin Y, Song Y, Liu Y, Wang Y, et al. Climatic characteristics and major meteorological events over China in 2017. Meteorological Monthly. 2018;44(4):548–55. (in Chinese with English abstract).

    Google Scholar 

  40. Jiang H, Dian W, Wu P. Effect of high temperature on fine structure of amylopectin in rice endosperm by reducing the activity of the starch branching enzyme. Phytochemistry. 2003;63:53–9.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  41. Lin CJ, Li CY, Lin SK, Yang FH, Huang JJ, Liu YH, Lur HS. Influence of high temperature during grain filling on he accumulation of storage proteins and grain quality in rice (Oryza sativa L.). J Agric Food Chem. 2010;58:10545–52.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  42. Zhang C, Zhou L, Zhu Z, Lu H, Zhou X, Qian Y, et al. Characterization of grain quality and starch fine structure of two japonica rice (Oryza Sativa) cultivars with good sensory properties at different temperatures during the filling stage. J Agric Food Chem. 2016;64:4048–57.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  43. Hittalmani S, Huang N, Courtois B, Venuprasad R, Shashidhar H, Zhuang J, et al. Identification of QTL for growth and grain yield–related traits in rice across nine locations of Asia. Theor Appl Genet. 2003;107:679–90.

    PubMed  Article  PubMed Central  Google Scholar 

  44. Bao J, Xia Y. Genetic control of the paste viscosity characteristics in indica rice (Oryza sativa L.). Theor Appl Genet. 1999;98:1120–4.

    Article  Google Scholar 

  45. Wang J, Li H, Zhang L, Meng L. 2013. Users’ Manual of QTL IciMapping. The Quantitative genetics group, Institute of crop science, Chinese Academy of Agricultural Sciences (CAAS), Beijing 100081, China, and Genetic Resources Program, International Maize and Wheat Improvement Center (CIMMYT), Mexico.

  46. Eshed Y, Zamir D. An introgression line population of Lycopersican pennelli in the cultivated tomato enables the identification and fine mapping of yield associated QTL. Genetics. 1995;141(3):1147–62.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  47. McCouch S, Cho Y, Yano M, Paul E, Blinstrub M, Morishima H, et al. Report on QTL nomenclature. Rice Genetics Newslett. 1997;14:11–3.

    Google Scholar 

Download references


We thank Dr. Wunyin Zhu and Miss Jin Lin for their works during the development of CSSLs.


This work was supported by grants from the Natural Science Foundation of Jiangsu Province of China (BK20180302), the Agriculture Research System of China (CARS-01-62), Jiangsu Agriculture Science and Technology Innovation Fund (CX(18)1001), and the Provincial Key Laboratory of Agrobiology, Jiangsu Academy of Agricultural Sciences, China. The funding bodies played no role in the design of the study, collection, analysis, and interpretation of data and in writing the manuscript.

Author information

Authors and Affiliations



LZ, CZ and CW designed the experiments. LZ and SY prepared powers of the grain and measured the RVA profile. LZ, TC, QZ, ZZ, YZ developed the CSSLs populations and CZ identified the CSSLs by molecular markers. TC and ZZ managed the planting of materials in field. QZ performed analyses. LZ wrote the manuscript which was modified by YZ. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Yadong Zhang or Cailin Wang.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Table S1.

Correlation analysis of seven RVA parameters in 2013 and 2014. Table S2. Correlation analysis of seven RVA parameters in 2016 and 2017. Table S3. Correlation analysis of starch RVA parameters.

Additional file 2: Table S4.

Genotype and phynotype of CSSLs.

Additional file 3.

Sequence Alignments of SSRGs between 9311 and Nipponbare.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Zhao, L., Zhao, C., Zhou, L. et al. QTL mapping for starch paste viscosity of rice (Oryza sativa L.) using chromosome segment substitution lines derived from two sequenced cultivars with the same Wx allele. BMC Genomics 22, 596 (2021).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Rice (Oryza sativa L.)
  • Starch paste viscosity
  • Quantitative trait locus
  • Chromosome segment substitution lines