Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Dynamic Quantitative Trait Loci Analysis of Seed Reserve Utilization during Three Germination Stages in Rice

  • Xinxin Cheng ,

    Contributed equally to this work with: Xinxin Cheng, Jinping Cheng

    Affiliation The Laboratory of Seed Science and Technology, State Key Laboratory of Crop Genetics and Germplasm Enhancement, Nanjing Agricultural University, Nanjing, PR China

  • Jinping Cheng ,

    Contributed equally to this work with: Xinxin Cheng, Jinping Cheng

    Affiliation The Laboratory of Seed Science and Technology, State Key Laboratory of Crop Genetics and Germplasm Enhancement, Nanjing Agricultural University, Nanjing, PR China

  • Xi Huang,

    Affiliation The Laboratory of Seed Science and Technology, State Key Laboratory of Crop Genetics and Germplasm Enhancement, Nanjing Agricultural University, Nanjing, PR China

  • Yanyan Lai,

    Affiliation The Laboratory of Seed Science and Technology, State Key Laboratory of Crop Genetics and Germplasm Enhancement, Nanjing Agricultural University, Nanjing, PR China

  • Ling Wang,

    Affiliation The Laboratory of Seed Science and Technology, State Key Laboratory of Crop Genetics and Germplasm Enhancement, Nanjing Agricultural University, Nanjing, PR China

  • Wenli Du,

    Affiliation The Laboratory of Seed Science and Technology, State Key Laboratory of Crop Genetics and Germplasm Enhancement, Nanjing Agricultural University, Nanjing, PR China

  • Zhoufei Wang ,

    wzf@njau.edu.cn (ZW); hszhang@njau.edu.cn (HZ)

    Affiliation The Laboratory of Seed Science and Technology, State Key Laboratory of Crop Genetics and Germplasm Enhancement, Nanjing Agricultural University, Nanjing, PR China

  • Hongsheng Zhang

    wzf@njau.edu.cn (ZW); hszhang@njau.edu.cn (HZ)

    Affiliation The Laboratory of Seed Science and Technology, State Key Laboratory of Crop Genetics and Germplasm Enhancement, Nanjing Agricultural University, Nanjing, PR China

Abstract

In this study, one rice population of recombinant inbred lines (RILs) was used to determine the genetic characteristics of seed reserve utilization during the early (day 6), middle (day 10) and late (day 14) germination stages. The seedling dry weight (SDW) and weight of the mobilized seed reserve (WMSR) were increased, while the seed reserve utilization efficiency (SRUE) decreased, during the process of seed germination. The SDW and WMSR were affected by the seed weight, while the SRUE was not affected by the seed weight. A total of twenty unconditional and twenty-one conditional additive QTLs and eight epistatic QTLs were identified at three germination stages, and the more QTLs were expressed at the late germination stage. Among them, twelve additive and three epistatic QTLs for SDW, eight additive and three epistatic QTLs for WMSR and thirteen additive and two epistatic QTLs for SRUE were identified, respectively. The phenotypic variation explained by each additive QTL, epistatic QTL and QTL × development interaction ranged from 6.10 to 23.91%, 1.79 to 6.88% and 0.22 to 2.86%, respectively. Two major additive QTLs qWMSR7.1 and qSRUE4.3 were identified, and each QTL could explain more than 20% of the total phenotypic variance. By comparing the chromosomal positions of these additive QTLs with those previously identified, eleven QTLs might represent novel genes. The best four cross combinations of each trait for the development of RIL populations were selected. The selected RILs and the identified QTLs might be applicable to improve rice seed reserve utilization by the marker-assisted selection approach.

Introduction

Vigorous seedling growth is essential for profitable and sustainable crop production. Seedling establishment begins with the appearance of the radicle and terminates when the seedling has exhausted the seed's energy reserves and starts to carry out photosynthesis [1]. Germination and heterotrophic growth are therefore crucial steps for the establishment of crops. By definition, seed vigor is the ability of a seed lot to establish seedlings successfully under a wide range of conditions [2]. High seed vigor, i.e., the rapid, uniform and complete emergence of vigorous seedlings, leads to high competition against weeds, uniform harvest and a high grain yield potential in rice. Improving rice seed vigor is an important objective in rice breeding programs.

Seed vigor is a complex trait that is determined by the interactions between genetic and environmental factors [3]. To date, several technological solutions, such as seed selection and seed treatment, have been used to improve seed vigor [2], [4][7]. In recent years, the emphasis on improving seed vigor has shifted to the potential of genetic improvement in breeding programs. Quantitative trait loci (QTL) analysis for seed vigor had been reported in rice, tomato and lettuce [8][13]. Those QTLs can be used to improve seed vigor using a marker-assisted selection (MAS) breeding program; furthermore, it becomes possible to identify candidate genes associated with seed vigor [14].

As seeds imbibe water, the mobilization of seed reserves such as starch, proteins, and lipids from storage tissues or the endosperm begins, providing essential energy to fuel germination and seedling growth until the seedling becomes photoautotrophic [15]. Heterotrophic seedling growth (seedling dry weight; SDW) can be defined as a product of two key traits: the weight of the mobilized seed reserve (WMSR) and the conversion efficiency of the utilized seed reserve (seed reserve utilization efficiency; SRUE) into seedling tissue [16][17]. These traits of seed reserve utilization play important roles in high seed vigor to improve germination and heterotrophic growth. Several phenotypic studies of seed reserve utilization have been conducted in wheat [16], soybean [17] and chickpea [18], but no studies have been reported in rice. Furthermore, the genetic control of seed reserve utilization in rice remains uncharacterized.

In rice, QTLs for seed vigor traits such as seed germination, seedling growth, seed longevity and tolerance to stress have been identified [8], [11][13], [19][25]. In most of these reports, the phenotypic values of the final developmental stage were used for QTL analysis. These QTLs with unconditional genetic effects (defined by Sen and Churchill [26]) are effective throughout seed germination. According to the theory of developmental genetics, genes are expressed selectively at different germination stages. QTL analysis can be adapted to include the effects of developmental stages [27]. Therefore, the QTL with conditional genetic effects (defined by Sen and Churchill [26]) are the loci detected at a specific growth stage [28]. The distinct gene actions for seed vigor at different developmental stages have been ignored in previous studies. Identification of both conditional and unconditional QTL will be desirable for MAS [28]. It will be essential, therefore, to include the dynamics of gene expression when analyzing the traits of seed reserve utilization.

Recently, several reports have explored the developmental behavior of various quantitative traits in rice [29][30], maize [31], cotton [32], soybean [28] and potato [33]. In rice, the dynamic QTL analysis has been used to dissect the genetic control of several quantitative traits such as tiller number [29][30], grain-filling [34] and root morphological traits [35]. The objectives of the present study were to investigate the developmental behavior of seed reserve utilization in rice, to identify both conditional and unconditional QTL and to measure the epistasis effects on seed reserve utilization during seed germination.

Materials and Methods

Plant materials

Two rice varieties, Jiucaiqing (japonica) and IR26 (indica), and 150 recombinant inbred lines (RIL) (F2:9) derived from Jiucaiqing × IR26 were used in this study [19][20]. All materials were grown at the Jiangpu Experimental Station of Nanjing Agricultural University, and the seeds were harvested at the mature stage. After the release of seed dormancy by post-ripening after four months of storage at room temperature (22–26°C), the seeds were used for the germination test.

Evaluation of seed reserve utilization

To determine the seed water content (WC) of each RIL and their parents, three replicates of seed (∼4.0 g per replicate) were powdered and dried at 104°C for 24 h. For the evaluation of seed reserve utilization, three replicates of each RIL and their parents (thirty seeds per replicate) were used and weighed (W). According to the corresponding WC, the initial seed dry weight (W0) (ISDW; mg per seed) of each RIL and their parents were calculated as [W×(1−WC)/30].

Then, the seed germination was conducted in a growth chamber at 30±1°C in the dark for 6, 10 or 14 days according to the methods of Zhang et al. [24]. The shoots and roots of seedlings were separated from the seeds by hand after germination. The seedlings dry weight (W1) (SDW; mg per seed) and the remnant seed dry weight (W2) (RSDW; mg per seed) were obtained after oven drying at 104°C for 24 h. The weight of the mobilized seed reserve (W3) (WMSR; mg per seed) was calculated as [W0–W2]. The seed reserve utilization efficiency (SRUE; mg per mg) was calculated as [W1/W3] [16].

Data analysis

Experimental data were analyzed using the Statistical Analysis System (SAS) software, and the traits of the parents were compared using Student's t-test at the 5% and 1% levels of probability. The correlations of the traits were computed using the PROC CORR by SAS software [19][20].

QTL mapping

The linkage map has been constructed previously [19][20]. Finally, a set of 168 SSR markers covering most of the rice genome at an average interval of 15 cM was constructed. Unconditional QTL was detected using the mean phenotypic values at day 6, 10 and 14 of germination. The genetic behavior measured at time t is the result of genes expressed before time (t-1) and extra effects within the period from (t-1) to t. Conditional phenotypic values (t/t-1) at time t were given by subtracting the phenotypic means measured at time t-1 from the mean at time t [27], [36]. Conditional QTL was detected using the mean phenotypic values (t/t-1) in the intervals from day 0 to 6 (t6/t0), 6 to 10 (t10/t6) and 10 to 14 (t14/t10) of germination.

Using the phenotypic values at a single germination stage, the unconditional and conditional additive QTL was identified respectively by the inclusive composite interval mapping (ICIM) method in QTL IciMapping [37][38]. Briefly, the P values for entering variables (PIN) and removing variables (POUT) were set at 0.01 and 0.02, and the scanning step was 1 cM [19]. The LOD score of 2.5 was used as the threshold value to declare the presence of a putative additive QTL. The QTLNetwork program ver. 2.0, based on a mixed linear model [39] was used to identify the epistatic QTL for seed reserve utilization in joint analysis of the unconditional phenotypic values of three germination stages. Briefly, the values for testing window and filtration window were set at 10 cM, respectively, and the walking speed was 1 cM. The putative QTL detection was determined at 95% confidence level. The proportion of observed phenotypic variance explained by each additive and epistatic QTL and the corresponding additive effects were also estimated. The QTL nomenclature followed the method of McCouch and CGSNL [40].

Prediction for novel parental combination

The best cultivars with maximum phenotypic value and elite alleles might be used to design parental combinations for crop breeding [41]. According to the information of phenotypic values and alleles of detected loci among the RILs, the best four cross combinations of each trait were selected. Firstly, the RILs with relatively highest phenotypic values of SDW, WMSR or SRUE at three germination stages (Day 6, 10 and 14) were selected. Then, the positive alleles of the detected additive QTLs among the selected RILs were analyzed. Using Monte Carlo simulation experiments, a novel cultivar combination can be determined by maximizing the number of elite alleles [41]. Finally, the best four cross combinations of each trait were selected.

Results

Seed reserve utilization phenotypes

During the process of seed germination, the SDW and WMSR were increased among the parents and RILs, while the SRUE was decreased. There were significant differences in the SDW, WMSR and SRUE between two parents at the middle (Day 10) and late (Day 14) germination stages (Table 1). The SDW and WMSR of Jiucaiqing were significantly higher than those of IR26 at the middle and late germination stages. In contrast, IR26 had a significantly higher SRUE at the middle and late germination stages. There was a continuous frequency distribution and transgressive segregation in the SDW, WMSR and SRUE among RIL population at three germination stages.

thumbnail
Table 1. Phenotypic values of seed reserve utilization related traits among the parents and RILs population at three germination stages.

https://doi.org/10.1371/journal.pone.0080002.t001

Correlation of seed weight with seed reserve utilization related traits

There were significant positive correlations between the ISDW and SDW during three germination stages (Table 2). The correlations between the ISDW and the WMSR were significantly positive during the middle and late germination stages, but there were no significant correlations between ISDW and SRUE. These results indicated that the SDW and WMSR are affected by the seed dry weight, while not the SRUE.

thumbnail
Table 2. Correlation coefficients between seed reserve utilization related traits and seed weight at three germination stages.

https://doi.org/10.1371/journal.pone.0080002.t002

Unconditional and conditional QTL for SDW

A total of six unconditional additive QTLs were identified for the SDW at three germination stages (Table 3; Fig. 1). One QTL (qSDW1.2) and six QTLs (qSDW1.2, qSDW2, qSDW4.2, qSDW5, qSDW8.2 and qSDW9) were detected at day 10 and 14 stage, respectively. Among them, the qSDW1.2 was consistently detected in both germination stages. The phenotypic variance explained by each QTL ranged from 6.78 to 16.75%. The positive alleles of three QTLs (qSDW2, qSDW5 and qSDW9) from Jiucaiqing and three QTLs (qSDW1.2, qSDW4.2 and qSDW8.2) from IR26 contributed to the increase of SDW.

thumbnail
Figure 1. Unconditional and conditional QTL for seed reserve utilization related traits at three germination stages.

* The unconditional QTL at day 6 stage is same with the conditional QTL during day 0 to 6 stage. The circle size represents the magnitude of phenotypic variance explained by each QTL. Yellow circle indicates that the positive allele of Jiucaiqing could increase the value of traits. Green circle indicates that the positive allele of IR26 could increase the value of traits.

https://doi.org/10.1371/journal.pone.0080002.g001

thumbnail
Table 3. Unconditional and conditional QTLs for seed reserve utilization related traits at three germination stages.

https://doi.org/10.1371/journal.pone.0080002.t003

Seven conditional additive QTLs were identified for the SDW during three germination stages (Table 3; Fig. 1). One QTL (qSDW12) and six QTLs (qSDW1.1, qSDW3, qSDW4.1, qSDW8.1, qSDW8.2 and qSDW8.3) were detected during day 6 to 10 stage and day 10 to 14 stage, respectively. The phenotypic variance explained by each QTL ranged from 6.97 to 16.18%. The positive alleles of three QTLs (qSDW1.1, qSDW3 and qSDW8.1) from Jiucaiqing and four QTLs (qSDW4.1, qSDW8.2, qSDW8.3 and qSDW12) from IR26 contributed to the increase of SDW.

Unconditional and conditional QTL for WMSR

Four unconditional additive QTLs were identified for the WMSR at three germination stages (Table 3; Fig. 1). Two QTLs (qWMSR2 and qWMSR4), one QTL (qWMSR1) and two QTLs (qWMSR3 and qWMSR4) were detected at day 6, 10 and 14 stage, respectively. Among them, the qWMSR4 was consistently detected in both day 6 and 14 stages. The phenotypic variance explained by each QTL ranged from 8.40 to 14.60%. The positive alleles of one QTL (qWMSR3) from Jiucaiqing and three QTLs (qWMSR1, qWMSR2 and qWMSR4) from IR26 contributed to the increase of WMSR.

Seven conditional additive QTLs were identified for the WMSR during three germination stages (Table 3; Fig. 1). Two QTLs (qWMSR2 and qWMSR4) and five QTLs (qWMSR1, qWMSR7.1, qWMSR7.2, qWMSR10 and qWMSR12) were detected during day 0 to 6 stage and day 6 to 10 stage, respectively. The phenotypic variance explained by each QTL ranged from 8.04 to 20.51%; the major QTL qWMSR7.1 accounted for the largest amount of phenotypic variation (20.51%) during the day 6 to 10 stage. The positive alleles of two QTLs (qWMSR7.1 and qWMSR10) from Jiucaiqing and five QTLs (qWMSR1, qWMSR2, qWMSR4, qWMSR7.2 and qWMSR12) from IR26 contributed to the increase of WMSR.

Unconditional and conditional QTL for SRUE

Ten unconditional additive QTLs were identified for the SRUE at three different germination stages (Table 3; Fig. 1). Four QTLs (qSRUE2.1, qSRUE2.2, qSRUE8.2, and qSRUE12.2), four QTLs (qSRUE4.1, qSRUE4.2, qSRUE4.3 and qSRUE12.2) and three QTLs (qSRUE4.4, qSRUE4.5 and qSRUE8.1) were detected at day 6, 10 and 14 stage, respectively. Among them, the qSRUE12.2 was consistently detected in both day 6 and 10 stages. The phenotypic variance explained by each QTL ranged from 6.69 to 23.91%; the major QTL qSRUE4.3 accounted for the largest amount of phenotypic variation (23.91%) at the day 10 stage. The positive alleles of five QTLs (qSRUE2.2, qSRUE4.1, qSRUE4.2, qSRUE4.5 and qSRUE8.1) from Jiucaiqing and five QTLs (qSRUE2.1, qSRUE4.3, qSRUE4.4, qSRUE8.2 and qSRUE12.2) from IR26 contributed to the increase of SRUE.

Seven conditional additive QTLs were identified for the SRUE during three germination stages (Table 3; Fig. 1). Four QTLs (qSRUE2.1, qSRUE2.2, qSRUE8.2 and qSRUE12.2) and three QTLs (qSRUE1, qSRUE3 and qSRUE12.1) were detected during day 0 to 6 stage and day 10 to 14 stage, respectively. The phenotypic variance explained by each QTL ranged from 6.10 to 18.65%. The positive alleles of three QTLs (qSRUE2.2, qSRUE3 and qSRUE12.1) from Jiucaiqing and four QTLs (qSRUE1, qSRUE2.1, qSRUE8.2 and qSRUE12.2) from IR26 contributed to the increase of SRUE.

Epistatic QTL for seed reserve utilization

A total of eight epistatic QTLs were identified by the joint analysis of the unconditional phenotypic values at three germination stages (Day 6, 10 and 14) (Table 4; Fig. 2). Among them, three epistatic QTLs each for SDW and WMSR were identified, respectively, with an epistatic main effect. The SRUE was controlled by two epistatic QTLs: one with both an epistatic main effect and an epistasis × development interaction effect and one with only an epistatic main effect. The phenotypic variance explained by each epistatic QTL ranged from 1.79 to 6.88%, and the phenotypic variation explained by each epistatic QTL × development interaction ranged from 0.22 to 2.86%.

thumbnail
Figure 2. Location of epistatic QTLs for seed reserve utilization related traits on linkage groups in joint analysis of three germination stages.

SDW: seedling dry weight; WMSR: weight of mobilized seed reserve; SRUE: seed reserve utilization efficiency. The chromosome names are listed on the top left of the bar; the map distances (cM) are shown below the bar and the markers are shown above the bar.

https://doi.org/10.1371/journal.pone.0080002.g002

thumbnail
Table 4. The epistatic QTL for seed reserve utilization related traits in joint analysis of three germination stages.

https://doi.org/10.1371/journal.pone.0080002.t004

Prediction for novel parental combination

The best four cross combinations of each trait for the development of RIL populations were selected to improve seed reserve utilization (Table 5). Four RILs each for SDW, WMSR and SRUE were selected with a relatively high value at three germination stages respectively (Fig. 3). The selected RILs had two to eight positive alleles of the additive QTLs for each trait (Table 6). To improve the three traits (SDW, WMSR and SRUE) simultaneously, 33 elite alleles could be pyramided by the combination of RILs 1750, 1791, 1878, and 1890 (Table 6).

thumbnail
Figure 3. The phenotypes of seed reserve utilization related traits among the selected RILs and two parents at three germination stages.

https://doi.org/10.1371/journal.pone.0080002.g003

thumbnail
Table 5. Parental combinations predicted from RILs of seed reserve utilization related traits.

https://doi.org/10.1371/journal.pone.0080002.t005

thumbnail
Table 6. The related additive QTLs with positive alleles in selected RILs.

https://doi.org/10.1371/journal.pone.0080002.t006

Discussion

Seed germination is shortly followed by the mobilization of seed reserves from the seed storage organs or endosperm, providing essential energy to fuel growth until the seedling becomes photoautotrophic [15]. Because of its importance in seed germination and post-germinative growth, the mobilization of reserves such as starch, protein and lipids during germination has been extensively studied [15], [42][44]. In this study, we focused on the three traits (SDW, WMSR and SRUE) related to seed reserve utilization. The seedling establishment of the most RILs begins at day 6 of germination with the appearance of the radicle and shoot and the seed's energy reserves have mostly been exhausted at day 14 of germination. Therefore, the traits of seed reserve utilization were conducted at three germination stages (Day 6, 10 and 14). Those traits were identified in the dark to avoid the photosynthesis for energy usage in seed germination.

During seed germination, the dry weight of the growing seedling is always lower than that of the mobilized reserves due to respiration [16]. Similar results were observed in our results. Therefore, except for the WMSR, the SRUE, as the conversion efficiency of the mobilized reserve into seedling tissue, might be an important factor influencing vigorous seedling growth. In this study, we found that there were significant cultivar differences for seed reserve utilization that are in agreement with the findings of Soltani et al. [45] in wheat. The japonica variety Jiucaiqing had the higher WMSR, while the indica variety IR26 had the higher SRUE in the middle and late germination stages. Among the RIL population, the genetic difference in the SDW, WMSR and SRUE were also observed. These results indicate that genetic variation can be used in rice breeding programs for the improvement of seed reserve utilization.

In this study, we found that large seeds have the advantage of producing more vigorous seedlings, similar to previous results [16], [46][47]. The ISDW was significantly corrected with the SDW and WMSR, while no significant relationship between ISDW and SRUE. Due to the SDW affected not only the WMSR but also the SRUE in germinating seed, the SDW of small seed with higher SREU might be similar that of larger seed with low SRUE. Therefore, the lack of a significant positive relation between seed weight and seedling dry weight has also been reported in several studies [48][50].

Improving the trait of rice seed reserve utilization using conventional strategies is extremely difficult because it is a quantitative trait affected by both genetic and environmental factors. It is therefore necessary to perform QTL analysis of seed reserve utilization for marker assisted selection. Unconditional QTL mapping can explain cumulative gene actions from the beginning to the end of a specific process, and conditional mapping can identify the QTLs that act during a specific growth period but are not affected by the genes expressed in the previous stages [36]. In this study, a total of 20 unconditional QTLs were identified. Among them, 6, 6 and 11 unconditional QTLs were identified at day 6, 10 and 14 of germination stages, respectively. Similarly, a total of 21 conditional QTLs were identified. Among them, 6, 6 and 9 conditional QTLs were identified during day 0 to 6, 6 to 10 and 10 to 14 germination stages, respectively. These results showed that more QTLs were expressed during the late stage than during the early and middle stages. Combined using the conditional and unconditional QTLs, we can improve the seed reserve utilization traits not only during a specific germination period but also at the end of a specific germination stage by MAS.

We found that the additive QTLs were rarely co-localized among the three germination stages. Only the qSDW1.2 was expressed at both day 10 and 14 germination stages, the qWMSR4 was expressed at both day 6 and 14 germination stages, and the qSRUE12.2 was expressed at both day 6 and 10 germination stages. These co-located QTLs might play important roles in controlling the seed reserve utilization during different germination stages. For the three traits of seed reserve utilization, the identification of co-localized additive QTLs only occurred in the genomic region RM6950–RM5759 of chromosome 1 for qSDW1.2 and qWMSR1, the genomic region RM240–RM5404 of chromosome 2 for qWMSR2 and qSRUE2.2 and the genomic region RM519–RM7376 of chromosome 12 for qSDW12 and qWMSR12. These co-localized QTLs could be very useful in the simultaneous improvement of more than one trait, as the desirable alleles of these additive QTLs (except for qSRUE2.2) were contributed by a single parent, IR26. Additionally, the joint analysis of the multi-development phenotypic values suggested that epistatic QTLs and QTL × development interactions were important components for seed reserve utilization, even though the degree of interaction was low.

By comparing the chromosomal positions of these additive QTLs, 36 QTLs identified in previous studies were found to map near the additive QTLs identified in this study (Fig. 4). The qSDW1.2 and qWMSR1 located to a position that coincided with the region of qSH1.2 for shoot height [20], and the qSDW1.1 coincided with the region of qTDW1-2 for total seedling dry weight and qRDW1-1 for root dry weight [8]. The qSDW2 was on the similar location of qGP-2 for germination percentage [12], and the qSDW3 was similar with the location of qIR-3 for imbibition rate [12] and q7GR for germination rate [24]. The qSDW5 mapped near the region of qTDW5-2 for total dry weight and qTAA5-1 for total amylase activity [8], and the qSDW8.3 was on the similar location of qIR-8 for imbibition rate [12]. The qSDW9 identified here was near the region of qRA9-1 for root activity and qSW9-1 for seed weight [8], and the qSDW12 and qWMSR12 was near the region of qSH12.1 for shoot height [20] and qIR-12 for imbibition rate [12].

thumbnail
Figure 4. The co-located QTLs in this study with previously mapped QTLs in rice.

The mapped QTLs in this study are shown on the left (red) and previously mapped QTLs are shown on the right (green) [8], [11][13], [19][25], [51][53].

https://doi.org/10.1371/journal.pone.0080002.g004

Similarly, the qWMSR2 and qSRUE2.2 was on the similar location of qGR-2 for germination rate [11], gw2.1 for 1000-seed weight [51] and qGW2d for grain width [52]. The qWMSR3 located to a position that coincided with the region of qLTG-3 for low temperature germinability [25]. The qWMSR7.1 was on the similar location of qGR7-1 [8] and q7GR [24] for germination rate, qGI-7 for germination index [11] and qGP-7-1 for germination percentage [12]. The qWMSR10 mapped near the region of qIR-10 for imbibition rate [12].

Additionally, the qSRUE1 identified here was near the region of qTDW1-1 for total dry weight, qSDW1-1 for shoot dry weight, qSW1-1 for seed weight [8] and qSH1.3 for shoot height [20]. The qSRUE3 was near the region of qLTG-3-2 for low temperature germinability [22]. The qSRUE4.1 was on the similar location of qGP-4-1 for germination rate [12] and qLTG-4 for low temperature germinability [22], and the qSRUE4.2 was similar with the region of qSH4 for shoot height [20]. The qSRUE4.4 located on the same region of qSD4 for seed dormancy [53], qLTG-4 for low temperature germinability [25] and qGP-4 for germination percentage [11]. The qSRUE4.5 was on the similar location of qGP-4-2 for germination percentage and qIR-4 for imbibition rate [13]. The qSRUE12.2 mapped near the region of qDSW12.1 for seedling dry weight [12]. However, there were no QTLs previously reported to be close to qSDW4.1, qSDW4.2, qSDW8.1, qSDW8.2, qWMSR4, qWMSR7.2, qSRUE2.1, qSRUE4.3, qSRUE8.1, qSRUE8.2, and qSRUE12.1, which indicates that these additive QTLs might be novel (Fig. 4).

Molecular evidence indicates that genes related to reserve mobilization and endosperm weakening could possibly affect seed germination [54][55]. Recently, several major QTLs for seed germination have been identified. For example, Wang et al. [11] identified one major QTL (qGI-11) for seed germination speed. Wang et al. [12] identified one major QTL (qGR-2) for seed germination under salt stress conditions. Recently, examination of the mechanism of seed vigor in rice has made significant progress. Fujino et al. [22] mapped a major rice QTL (qLTG3-1) for low-temperature germinability on chromosome 3. Fujino et al. [23] placed qLTG3-1 in a 4.8-kb region of chromosome 3 using high-resolution mapping and revealed that the expression of qLTG3-1 was tightly associated with the vacuolation of the tissues covering the embryo. This vacuolation may cause tissue weakening, resulting in the reduction of the mechanical resistance to the growth potential of the coleoptile. Seed reserve utilization is a very complex physiological process during seed germination. In this study, two major QTLs qWMSR7.1 and qSRUE4.3 were identified, and each QTL could explain more than 20% of the total phenotypic variance. Fine mapping the two QTLs are now in progress to elucidate the molecular mechanism of seed reserve utilization using near isogenic lines (NILs). Furthermore, the linked marker with the two major QTLs might be useful in marker-assisted selection to increase the seed reserve utilization in rice.

In a hypothetical cross of two cultivars, the trait values of produced RILs can be predicted by the effects of all the detected loci. The best RIL with maximum value would represent the best cross [41]. In this study, the best four cross combinations of each trait for the development of RIL populations were selected. It was found that some RILs present repeatedly in the novel parental combinations: for example, RIL 1750 in SDW and WMSR and RIL 1890 in WMSR and SRUE, indicating high general combining ability. To improve multiple traits, all the elite alleles might be pyramided into one cultivar as far as possible [41]. In this study, a total of 33 elite alleles could be pyramided by the combination of four RILs to improve the seed reserve utilization traits. These results showed that the selected RILs as new materials will be valuable in future rice breeding programs. The identified QTLs will be applicable to improve rice seed reserve utilization by MAS.

Author Contributions

Conceived and designed the experiments: ZW HZ. Performed the experiments: XC JC. Analyzed the data: JC XH YL. Contributed reagents/materials/analysis tools: LW WD. Wrote the paper: ZW.

References

  1. 1. Ichie T, Ninomiya I, Ogino K (2001) Utilization of seed reserves during germination and early seedling growth by Dryobalanops lanceolata (Dipterocarpaceae). J Trop Ecol 17: 371–378.
  2. 2. Finch-Savage WE, Clay HA, Lynn JR, Morris K (2010) Towards a genetic understanding of seed vigour in small-seeded crops using natural variation in Brassica oleracea. Plant Sci 179: 582–589.
  3. 3. Hodgkin T, Hegarty TW (1978) Genetically determined variation in seed germination and field emergence of Brassica oleracea. Ann Appl Biol 88: 407–413.
  4. 4. Guan YJ, Hu J, Wang XJ, Shao CX (2009) Seed priming with chitosan improves maize germination and seedling growth in relation to physiological changes under low temperature stress. J Zhejiang Univ Sci B 10: 427–433.
  5. 5. Schwember AR, Bradford KJ (2010) A genetic locus and gene expression patterns associated with the priming effect on lettuce seed germination at elevated temperatures. Plant Mol Biol 73: 105–118.
  6. 6. Chen K, Arora R (2011) Dynamics of the antioxidant system during seed osmopriming, post-priming germination, and seedling establishment in Spinach (Spinacia oleracea). Plant Sci 180: 212–220.
  7. 7. Zhang M, Wang Z, Yuan L, Yin C, Cheng J, et al. (2012) Osmopriming improves tomato seed vigor under aging and salinity stress. Afr J Biotechnol 11: 6305–6311.
  8. 8. Cui K, Peng S, Xing Y, Xu C, Yu S, et al. (2002) Molecular dissection of seedling-vigor and associated physiological traits in rice. Theor Appl Genet 105: 745–753.
  9. 9. Foolad MR, Zhang LP, Subbiah P (2003) Genetics of drought tolerance during seed germination in tomato: inheritance and QTL mapping. Genome 46: 536–545.
  10. 10. Hayashi E, Aoyama N, Still DW (2008) Quantitative trait loci associated with lettuce seed germination under different temperature and light environments. Genome 51: 928–947.
  11. 11. Wang ZF, Wang JF, Bao YM, Wang FH, Zhang HS (2010) Quantitative trait loci analysis for rice seed vigor during the germination stage. J Zhejiang Univ Sci B 11: 958–964.
  12. 12. Wang Z, Wang J, Bao Y, Wu Y, Zhang H (2011) Quantitative trait loci controlling rice seed germination under salt stress. Euphytica 178: 297–307.
  13. 13. Wang Z, Wang J, Zhou R, Wang J, Zhang H (2011) Identification of quantitative trait loci for cold tolerance during the germination and seedling stages in rice (Oryza sativa L.). Euphytica 181: 405–413.
  14. 14. Alonso-Blanco C, Aarts MGM, Bentsink L, Keurentjes JJB, Reymond M, et al. (2009) What has natural variation taught us about plant development, physiology and adaptation. Plant Cell 21: 1877–1896.
  15. 15. Pritchard SL, Charlton WL, Baker A, Graham IA (2002) Germination and storage reserve mobilization are regulated independently in Arabidopsis. Plant J 31: 639–647.
  16. 16. Soltani A, Gholipoor M, Zeinali E (2006) Seed reserve utilization and seedling growth of wheat as affected by drought and salinity. Env Exp Bot 55: 195–200.
  17. 17. Mohammadi H, Soltani A, Sadeghipour HR, Zeinali E (2011) Effects of seed aging on subsequent seed reserve utilization and seedling growth in soybean. Int J Plant Prod 5: 65–70.
  18. 18. Soltani A, Galeshi S, Zeinali E, Latifi N (2002) Germination, seed reserve utilization and seedling growth of chickpea as affected by salinity and seed size. Seed Sci Technol 30: 51–60.
  19. 19. Wang Z, Chen Z, Cheng J, Lai Y, Wang J, et al. (2012) QTL analysis of Na+ and K+ concentrations in roots and shoots under different levels of NaCl stress in rice (Oryza sativa L.). PLoS ONE 7: e51202.
  20. 20. Wang Z, Cheng J, Chen Z, Huang J, Bao Y, et al. (2012) Identification of QTLs with main, epistatic and QTL × environment interaction effects for salt tolerance in rice seedlings under different salinity conditions. Theor Appl Genet 125: 807–815.
  21. 21. Miura K, Lin S, Yano M, Nagamine T (2002) Mapping quantitative trait loci controlling seed longevity in rice (Oryza sativa L.). Theor Appl Genet 104: 981–986.
  22. 22. Fujino K, Sekiguchi H, Sato T, Kiuchi H, Nonoue Y, et al. (2004) Mapping of quantitative trait loci controlling low-temperature germinability in rice (Oryza sativa L.). Theor Appl Genet 108: 794–799.
  23. 23. Fujino K, Sekiguchi H, Matsuda Y, Sugimoto K, Ono K, et al. (2008) Molecular identification of a major quantitative trait locus, qLTG3-1, controlling low-temperature germinability in rice. Proc Natl Acad Sci USA 105: 12623–12628.
  24. 24. Zhang ZH, Yu SB, Yu T, Huang Z, Zhu YG (2005) Mapping quantitative trait loci (QTLs) for seedling-vigor using recombinant inbred lines of rice (Oryza sativa L.). Field Crops Res 91: 161–170.
  25. 25. Jiang L, Liu S, Hou M, Tang J, Chen L, et al. (2006) Analysis of QTLs for seed low temperature germinability and anoxia germinability in rice (Oryza sativa L.). Field Crops Res. 98: 68–75.
  26. 26. Sen S, Churchill GA (2001) A statistical framework for quantitative trait mapping. Genetics 159: 371–387.
  27. 27. Zhu J (1995) Analysis of conditional genetic effects and variance components in developmental genetics. Genetics 141: 1633–1639.
  28. 28. Han Y, Xie D, Teng W, Zhang S, Chang W, et al. (2011) Dynamic QTL analysis of linolenic acid content in different developmental stages of soybean seed. Theor Appl Genet 122: 1481–1488.
  29. 29. Liu G, Zhu H, Liu S, Zeng R, Zhang Z, et al. (2010) Unconditional and conditional QTL mapping for the developmental behavior of tiller number in rice (Oryza sativa L.). Genetica 138: 885–893.
  30. 30. Liu G, Zhu H, Zhang G, Li L, Ye G (2012) Dynamic analysis of QTLs on tiller number in rice (Oryza sativa L.) with single segment substitution lines. Theor Appl Genet 125: 143–153.
  31. 31. Zhang Z, Liu Z, Cui Z, Hu Y, Wang B, et al. (2013) Genetic analysis of grain filling rate using conditional QTL mapping in maize. PLoS ONE 8: e56344.
  32. 32. Ye Z, Lu Z, Zhu J (2003) Genetic analysis for developmental behavior of some seed quality traits in upland cotton. Euphytica 129: 183–191.
  33. 33. Li J, Lindqvist-Kreuze H, Tian Z, Liu J, Song B, et al. (2012) Conditional QTL underlying resistance to late blight in a diploid potato population. Theor Appl Genet 124: 1339–1350.
  34. 34. Takai T, Fukuta Y, Shiraiwa T, Horie T (2005) Time-related mapping of quantitative trait loci controlling grain-filling in rice (Oryza sativa L.). J Exp Bot 56: 2107–2118.
  35. 35. Qu Y, Mu P, Zhang H, Chen CY, Gao Y, et al. (2008) Mapping QTLs of root morphological traits at different growth stages in rice. Genetica 133: 187–200.
  36. 36. Yan J, Zhu J, He C, Benmoussa M, Wu P (1998) Quantitative trait loci analysis for the developmental behavior of tiller number in rice (Oryza sativa L.). Theor Appl Genet 97: 267–274.
  37. 37. Li H, Ye G, Wang J (2007) A modified algorithm for the improvement of composite interval mapping. Genetics 175: 361–374.
  38. 38. Li H, Ribaut JM, Li Z, Wang J (2008) Inclusive composite interval mapping (ICIM) for digenic epistasis of quantitative traits in biparental populations. Theor Appl Genet 116: 243–260.
  39. 39. Yang J, Hu CC, Ye XZ, Zhu J (2005) QTLNetwork 2.0. Available: http://ibi.zju.edu.cn/software/qtlnetwork. Accessed 2013 Jun 3.
  40. 40. McCouch SR, CGSNL (Committee on Gene Symbolization, Nomenclature, Linkage, Rice Genetics Cooperative) (2008) Gene nomenclature system for rice. Rice 1: 72–84.
  41. 41. Niu Y, Xu Y, Liu XF, Yang SX, Wei SP, et al. (2013) Association mapping for seed size and shape traits in soybean cultivars. Mol Breeding 31: 785–794.
  42. 42. Penfield S, Rylott EL, Gilday AD, Graham S, Larson TR, et al. (2004) Reserve mobilization in the Arabidopsis endosperm fuels hypocotyl elongation in the dark, is independent of abscisic acid, and requires PHOSPHOENOLPYRUVATE CARBOXYKINASE1. Plant Cell 16: 2705–2718.
  43. 43. Aoki N, Scofield GN, Wang XD, Offler CE, Patrick JW, et al. (2006) Pathway of sugar transport in germinating wheat seeds. Plant Physiol 141: 1255–1263.
  44. 44. Leonova S, Grimberg A, Marttila S, Stymne S, Carlsson AS (2010) Mobilization of lipid reserves during germination of oat (Avena sativa L.), a cereal rich in endosperm oil. J Exp Bot 61: 3089–3099.
  45. 45. Soltani A, Zeinali E, Galeshi S, Latifi N (2001) Genetic variation for and interrelationships among seed vigor traits in wheat from the Caspian Sea coast of Iran. Seed Sci Technol 29: 653–662.
  46. 46. Douglas CL, Wilkins DE, Churchill DB (1994) Tillage, seed size and seed density effects on performance of soft white winter wheat. Agron J 86: 707–711.
  47. 47. Bockous WW, Shroyer JP (1996) Effect of seed size on seedling vigor and forage production of winter wheat. Can J Plant Sci 76: 101–105.
  48. 48. Carleton AE, Cooper CS (1972) Seed size effects upon seedling vigor of three forage legumes. Crop Sci 12: 183–186.
  49. 49. Johnson DR, Wax LM (1978) Relationships of soybean germination and vigor tests to field performance. Agron J 70: 273–278.
  50. 50. Shroyer JP, Cox TS (1984) Effects of cultivar, environment and their interaction on seed quality of hard red winter wheat from production fields. J Appl Seed Prod 2: 24–28.
  51. 51. Marri PR, Sarla N, Reddy LV, Siddiq EA (2005) Identification and mapping of yield and yield related QTLs from an Indian accession of Oryza rufipogon. BMC Genet 6: 33.
  52. 52. Bai X, Luo L, Yan W, Kovi MR, Zhan W, et al. (2010) Genetic dissection of rice grain shape using a recombinant inbred line population derived from two contrasting parents and fine mapping a pleiotropic quantitative trait locus qGL7. BMC Genet 11: 16.
  53. 53. Gu XY, Kianian SF, Foley ME (2004) Multiple loci and epistases control genetic variation for seed dormancy in weedy rice (Oryza sativa). Genetics 166: 1503–1516.
  54. 54. Fait A, Angelovici R, Less H, Ohad I, Urbanczyk-Wochniak E, et al. (2006) Arabidopsis seed development and germination is associated with temporally distinct metabolic switches. Plant Physiol 142: 839–854.
  55. 55. Bethke PC, Libourel IG, Aoyama N, Chung YY, Still DW, et al. (2007) The Arabidopsis aleurone layer responds to nitric oxide, gibberellin, and abscisic acid and is sufficient and necessary for seed dormancy. Plant Physiol 143: 1173–1188.