Rice is considered to be one of the most important cereal crops, as it is the staple food for almost half of the world’s population. In recent years, there has been a growing interest from many rice-producing countries on improving the cooking and eating quality, two of the most important components of grain quality of rice (Juliano 1985; Unnevehr et al., 1992; Ge et al., 2005).
In the best-quality grain, the kernel elongates but its width changes little during cooking (Khush et al., 1979, which has been associated with the Basmati type (Sood et al., 1979; 1983; Ahn et al., 1993), an important commercial commodity mainly found in India, Pakistan, and Thailand (Khush et al., 1979).
For cooked rice grain quality, greater CRE was associated with high-quality rice varieties (Tang et al., 1987). Some previous studies had performed genetic analyses andgene mapping of these cooked rice properties. By studying a collection of upland rice varieties, Hussain et al (1987) found that CRE had rather high broad-sense heritability, suggesting it is more difficult to improve CRE genetically than for either milked rice length or cooked rice length. Bao et al(2001) confirmed that CRE phenotypes are dependent on a small number of major genes.Ahn et al mapped a major CRE QTL on chromosome 8 in Oryza sativa‘Basmati 370’, and consequently Li et al (2004) detected a second CRE QTL on chromosome 3. Chromosome 2, 5, 6, and 7 were found to exert further effects on CRE (Shen et al., 2005), and a CRE QTL closely associated with Wx gene on chromosome 6 had also been described by researchers (He et al., 2003; Zhang et al., 2004; Ge et al., 2005; Tian et al., 2005). In addition, Amarawathi et al (2008) used a RIL population to detectone QTL for CRE on chromosome 11. With these populations, such as F2, F2/F3, BC1, DH, and RIL, it was difficult to differentiate QTLs,especially ones that have small and/or interacting effects from background noise (Eshed and Zamir, 1995).
Eshed and Zamir (1995) proposed to apply introgression line (IL)populations so that a high-resolution mapping of QTLs can be achieved. Several permanent mapping populations in rice, such as CSSLs and backcross inbred lines (BILs), had been developed (Aida et al., 1997; Kubo et al., 1999; Xiao et al., 2005). For CSSLs, each chromosome segment substitution line (CSSL) carried a single or fewer donor segments in the background of a recurrent genotype. Interactions between donor alleles were limited to those genes on homozygous substituted tracts, thus reducing the effects of interferences from genetic background (Howell et al., 1996; Yano, 2001). Several CSSLs in rice had been developed and a lot of QTLs for traits with biological and economic interest had been detected through these CSSLs by various researchers (Kubo et al., 2002; Ebitani et al., 2005; Mei et al., 2006; Takai et al., 2007; Ye et al., 2010; Li et al., 2011; Ujiie et al., 2012; Marzougui et al., 2012). However, few studies have performed genetic analyses for CRE using CSSLs. Therefore, the objective of this study is to find specific QTLs for CRE using CSSLs, These QTLs should be helpful in marker-assisted selection (MAS) of high eating quality rice varieties and map-based cloning of desirable QTLs.
1 Results
1.1 Characteristics of the CSSLs
The substituted chromosome segments in the CSSLs covered most of the 12 chromosomes. Each line contained a substituted segment of a particular chromosomal region and additional small segments in non-target regions. If we assume that recombination occurred midway between two adjacent markers, the estimated length of the substituted segments in 103 CSSLs ranged from a minimum 1.4 cM to a maximum 79.2 cM, and had an average of 25.2 cM. The total length of the substituted segments in the CSSLs population was 2590.6 cM, which was 1.7 times that of rice genome (Zhu et al., 2009).
1.2 Traits of MRL, CRL, and CRE in parents
As illustrated in Table 1, significant phenotypic differences were detected between the parents for MRL, CRL, and CRE (t-test, P<0.01), and 93-11 showed higher level in three traits compared with Nipponbare, which indicated the approach used to evaluate the three traits in this study was valid.
Table 1 Differences of the parents among MRL, CRL, and CRE
|
1.3 Analysis of Quantitative trait loci
In total, three QTLs for MRL, designated as qMRL-3-1, qMRL-3-2, qMRL-8, were detected based on the difference oft-testbetween the mean of each CSSL and 93-11 and the additive effect ranged from -0.41 to -0.14. The QTL qMRL-3-1, with the largest effect, was located in the interval between RM6832 and RM3513 on chromosome 3. For all the three QTLs, the alleles from Nipponbare had negative effects (Table 2; Figure 1).
Table 2 QTLs and their additive effect on the substituted segments in CSSLs for MRL, CRL, and CRE |
Figure 1 Genetic map and chromosomal location of QTLs for MRL, CRL and CRE in rice. |
Four QTLs for CRL, qCRL-3-2, qCRL-6, qCRL-8, and qCRL-9, were detected based on the difference oft-testbetween the mean of each CSSL and 93-11, ranging from -0.41 to -0.16 for the additive effect. The QTL qCRL-6, with the largest effect, was located in the interval between RM508 and RM510 on chromosome 6. Similarly, the alleles from Nipponbare had negative effects for all four QTLs (Table 2; Figure 1).
Five QTLs for CRE, designated as qCRE-4, qCRE-6, qCRE-9, qCRE-10,and qCRE-11, were detected based on the difference oft-testbetween the mean of each CSSL and 93-11, and the rangeof the additive effect was from -5.80 to -2.82. The QTL qCRE-6with the largest effect was in the same region with Wx locus and located in the interval between RM508 and RM510 on chromosome 6. In the same way, for all five QTLs, the alleles from Nipponbare also had negative effects (Table 2; Figure 1).
1.4 Substitutional mapping QTLs
To map these QTLs, the substitution mapping was performed according to the method described by Wissuwa et al (2002). Using 5 CSSLs, qMRL-3-1 was mapped in the marker interval RM5551-RM6931, which spanned 8.5 cM in genetic distance on chromosome 3, while qMRL-3-2 and qCRL-3-2 were simultaneously located in the marker interval RM6832-RM3513, with the genetic distance of 5.3 cM (Figure 2A). Using CSSL43 and CSSL41, qCRE-4 was mapped in the marker interval RM5709-RM280 on chromosome 4 at the genetic distance of 9.5 cM ( Figure 2B). Using CSSL53 and CSSL52, qCRL-6 and qCRE-6 were both located in the marker interval RM508-RM510 on chromosome 6, and the genetic distance was 4.6 cM (Figure 2C). Using CSSL72, CSSL73, and CSSL74, qMRL-8 and qCRL-8 were located in the same marker interval RM4085-RM3395 on chromosome 8, which spanned 8.6 cM in genetic distance (Figure 2D). Using CSSL81 and CSSL83, qCRL-9 and qCRE-9 were simultaneously mapped in the marker interval RM566-RM242 on rice chromosome 9 at the genetic distance of 10.7 cM (Figure 2E). Using 5 CSSLs, qCRE-10 and qCRE-11 were respectively located in the marker interval RM8201-RM1108 on chromosome 10 and RM5349- RM2136 on chromosome 11, with the genetic distance of 14.3 cM and 27.9 cM (Figure 2F; Figure 2G).
Figure 2 Substitutional mapping QTLs for MRL, CRL and CRE in rice
|
2 Discussion
2.1 QTLs detected through CSSLs
In this study, we used CSSLs and the substitution mapping method to analyze additive effects of QTLs on MRL, CRL, and CRE in rice. Since each CSSL used in this study contains only one substituted segment from a donor in 93-11 genetic background, all the genetic variation between each CSSL and 93-11 was contributed by the substituted segment, thus providing more reliable QTL detection and estimation of QTL effects by minimizing the background genetic effects and. In addition, the smaller substituted chromosome segments in the CSSLs will facilitate the detection of QTLs in rice. Further experiments need to be conducted to detect QTLs elsewhere in the genome and/or determine whether some segments contained more than one QTL, as the substituted chromosome segments used in this study did not include the whole genome.
2.2 Three QTLs controlling CRE might be novel
QTLs affecting CRE were detected on all chromosomes except chromosome 9 in previous studies. We did not detect a major QTL on chromosome 8 associated with cooked kernel elongation of rice in this study as previously reported by Ahn et al., which was possibly due to the presence of a genomic segment introgression from Basmati 370 (a famous variety with more than 100% elongation in length and with little increase in width after cooking) in the germplasm used in their study. Zhang et al (2004) used a DH population to identify five QTLs for CRE on chromosome 1, 3, 5, 6, and 10, respectively. One QTL for CRE was mapped to the same region on chromosome 2 using different populations (Ge et al., 2005; Liu et al., 2008), while another CRE QTL was located in the marker interval RM301-RM29 on chromosome 2 using a DH population (Tian et al., 2005). Liu et al (2008) detected three CRE QTLs on chromosome 4, 5, and 12, respectively, and the qCRE-4 on chromosome 4 near qER-4 was detected in this study. Li et al (2004) mapped a CRE QTL on chromosome 3, with the favorable allele obtained from the African rice O. glaberrima. Using a RIL population, Wang et al (2007) identified four CRE QTLs on chromosome 3, 6, 7, and 8, respectively. Interestingly, a CRE QTL on chromosome 6 closely linked to the Wx gene had been simultaneously described by previous researches (He et al., 2003;Zhang et al., 2004; Ge et al., 2005; Tian et al., 2005). In this paper, the qCRE-6 with the largest effect was located in the near region with the Waxy gene on chromosome 6 (Figure 3), and the Waxy gene region played an important role in CRE.Recently, Amarawathi et al (2008) detected one QTL for CRE on chromosome 11 using a RIL population. Two QTLs on chromosomes 4 and 6 detected here were likely in common with QTLs detected in previous studies, but other QTLs (qCRE-9, qCRE-10, and qCRE-11) detected here had not been reported (Zhang et al., 2004; Tian et al., 2005; Liu et al., 2008; Amarawathi et al., 2008).
Figure 3 Comparison between the results in this study and in previous researches |
2.3 Pleiotropic QTLs in same location
Some QTLs affecting MRL, CRL, and CRE were mapped to the same genome regions (Figure 1). For example, the QTL clusters corresponding to the Wx locus simultaneously controlled CRE and CRL, and the locus controlling MRL and CRL simultaneously were both detected on chromosomes 3 and 8. Co-localization of these QTLs, affected by pleiotropic or close linkage may explain the genetic basis for correlations among various quality traits. Therefore, further genetic studies, such as fine mapping of these QTLs, which are required to determine the actual number of distinct QTLs in these regions and to elucidate the genetic mechanism of the co-localization are currently underway in secondary F2 populations derived from crosses between target CSSLs and 93-11.
2.4 Genetic Improvement of cooked rice elongation
Our results have clearly demonstrated that the region of waxy gene plays a major role in determining the CRE. The presence of closely linked markers on the side of the Wx locus will greatly facilitate the precise replacement of the alleles of the poor-quality parent using MAS. Wx alleles in donor parent a specialty Japonica with low MRL, low CRL, and low CRE, are used to replace alleles from the recipient parent 93-11, a specialty Indica rice cultivar in China with high MRL, high CRL, and high CRE (Table 1). Cooked recombinant rice stayed tender and slender, suggesting93-11 could be a desirable parent for improving the quality of rice by introducing its Wx alleles into the breeding lines.
3 Materials and methods
3.1 Developing the CSSLs
The strategy for the development of the CSSLs is shown in Figure 4. Rice cultivar 93-11 was crossed with Nipponbare, and a resultant F1 plant was backcrossed to 93-11 to produce BC1F1, which were again backcrossed with 93-11 to produce BC2F1. Subsequently, 135 BC3F1 lines were obtained from 1500 BC4F1 individuals.In addition, 230 SSR markers selected from dense rice microsatellite maps (McCouch et al., 2002) were used to detect the polymorphism between 93-11 and Nipponbare. A whole-genome survey of the 1500 BC4F1 using the polymorphic markers identified 139 plants with most of their genomic regions homozygous of 93-11, and these plants were self-crossed to produce BC4F2 segregation populations. Finally, a total of 103 CSSLs were selected in 6192 BC4F2 plants with MAS (Zhu et al.,2008).
Figure 4 Breeding strategy for developing CSSLs carrying Nipponbare segments in the genetic background of 93-11
|
3.2 Plant materials used in this research
The population used in this study consisted of 103 CSSLs, each line and both parent plants were field-grown in Jiangsu Academy of Agricultural Sciences (Nanjing, China) in 2008. Field management essentially followed normal agricultural practice. Before conducting quality analysis, all seeds samples were stored at room temperature for a period of at least 3 months after being harvested.
3.3 Assay of cooked rice elongation
The analysis of CRE was performed according to the method described by Wang et al(2007) with slight modifications. Ten intact milled grains were selected randomly from each CSSL and measured for MRL. The grains were soaked in distilled water for 30 min, which were then placed between two pieces of wet ï¬lter paper in a petridish filled with an appropriate amount of water. The dish was then placed in a covered container and the grains were cooked by steaming over boiling water for 10 min and simmering for 10 min (with power off). The cooked rice grain was transferred onto a piece of dry ï¬lter paper at the bottome of a fresh petridish, which was placed in a desiccator with a constant temperature (19℃). Then the CRL was measured. For repeated sample of eachCSSL, the two traits above were measured in duplicate, with values averaged over replications. TheCREwas calculated as follows: CRE = (CRL –MRL)/ MRL · 100%.
3.4 Determining the length of substituted segments in CSSLs
Rice microsatellite map was used to locate and determine the substituted segments in CSSLs (McCouch et al, 2002), the length of which was then estimated based on graphical genotypes (Young and Tanksley, 1989; Xi et al, 2006). 100% donor genotype was defined as a chromosome segment flanked by two markers of donor genotype (DD), whereas 0% donor genotype segment was flanked by 2 recipient genotype (RR) markers, and 50% donor genotype was flanked by one recipient marker (DR). Thus, the length of the substituted segments was estimated as the sum of the length of DD plus 1 half-length DR on each side.
The substituted segment was determined based on its location on rice microsatellite map. The length of substituted segments in CSSLs was estimated based on graphical genotypes. A chromosome segment flanked by two markers of donor genotype (DD) was considered as 100% donor genotype, a chromosome segment flanked by two markers of recipient genotype (RR) was considered as 0% donor genotype. A chromosome segment flanked by one marker of donor type and one marker of recipient type (DR) was considered as 50% donor type. In other words, the estimated length of a substituted segment was a sum of the length of DD and two of 1/2 length of DR both sides.
3.5 Substitutional mapping QTLs
QTLs were detected based on the difference oft-test between the mean of each CSSL and 93-11. A probability level of 0.01 was used as the threshold for the detection of a putative QTL. QTLs were named according to McCouch et al. (1997).We detected the exact QTL position by using the substitution mapping approach which was performed according to the method described by Wissuwa et al (2002). Moreover, the additive effect was estimated in accordance with the method described by Eshed et al (1995): The additive effect=phenotypic value of CSSLï¼phenotypic value of 93-11; The additive effect percentage=(the additive effect/phenotypic value of 93-11) ×100%.
Acknowledgement
This work was supported by the Key Support Program of Jiangsu Science and Technology (BE2008354) and Jiangsu Agriculture Science and Technology Innovation Fund (CX(12)1003).
Ahn S.N., Bollich C.N., McClung A.M., and Tanksley S.D., 1993, RFLP analysis of genomic regions associated with cooked-kernel elongation in rice, Theoretical and Applied Genetics, 87: 27-32
Aida Y., Tsunematsu H., Doi K., and Yosimura A., 1997, Development of a series of introgression lines of Japonica in the background of Indica rice, Rice Genetics Newsletters,14: 41-43
Amarawathi Y., Singh R., Singh A.K., Singh V.P., Mohapatra T., Sharma T.R., and Singh N.K., 2008, Mapping of quantitative trait loci for basmati quality traits in rice (Oryza sativa L.), Molecular Breeding, 21: 49-65
Bao J.S., Xie J.K., and Xia Y.W., 2001, Genetic analysis of cooked rice elongation in Indica rice (Oryza sative L). Acta Agronomica Sinica, 27: 489-492 (in chinese)
Ebitani T., Takeuchi Y., Nonoue Y., Yamamoto T., Takeuchi K., and Yano M., 2005, Construction and evaluation of chromosome segment substitution lines carrying overlapping chromosome segments of Indica rice cultivar ‘Kasalath’ in a genetic background of Japonica elite cultivar ‘Koshihikari’, Breeding Science, 55: 65-73
Eshed Y., and Zamir D., 1995, An introgression line population of Lycopersicon pennellii in the cultivated tomato enables the identification and fine mapping of yield-associated QTL, Genetics, 141: 1147-1162
Ge X.J., Xing Y.Z., Xu C.G., and He Y. Q., 2005, QTL analysis of cooked rice grain elongation, volume expansion, and water absorption using a recombinant inbred population, Plant Breeding, 124: 121-126
He Y.Q., Xing Y.Z., Ge X.J., Li X.H., and Xu C.G., 2003, Gene mapping for elongation index related traits on cooked rice grain quality, Molecular Plant Breeding, 1: 613-619 (in chinese)
Howell P.M., Lydiate D.J., Marshall D.F., 199, Towards developing intervarietal substitution lines in Brassica napus using marker-assisted selection, Genome, 39: 348-358
Hussain A.A., Maurya D.M., and Vaish C.P., 1987, Study on quality status of indigenous upland rice (Oryza. sativa L), Indian Journal of Genetics and Plant Breeding, 47: 145-152
Juliano B.O., 1985, Rice hull and rice straw, Rice Chemistry and Technology, 2d ed. American Association of Cereal Chemists. St. Paul, Minnesota, USA
Khush G.S., Paule C.M., and DelaCruz N.M., 1979, Rice grain quality evaluation and improvement at IRRI. In: Chemical Aspects of Rice Grain Quality, International Rice Research Institute, Los Banos, Philippines
Kubo T., Nakamura K., and Yoshimura A., 1999, Development of a series of indica chromosome segment substitution lines in japonica background of rice, Rice Genetics Newsletters, 16: 104-105
Kubo T., Aida Y., Nakamura K., Tsunematsu H., Doi K., and Yoshimura A., 2002, Reciprocal chromosome segment substitution series derived
from Japonica and Indica cross of rice (Oryza sativa L.), Breeding Science, 52: 319-325
Li J., Xiao J., Grandillo S., Jiang L., Wan Y., Deng Q., Yuan L., and McCouch S.R., 2004, QTL detection for rice grain quality traits using an interspecific backcross population derived from cultivated Asian (Oryza sative L.) and African (Oryza glaberrima S.) rice, Genome, 47: 697-704
Li M.,Sun P.L., Zhou H.J., Chen S., Yu S.B., 2011, Identification of quantitative trait loci associated with germination using chromosome segment substitution lines of rice (Oryza sativa L.), Theoretical and Applied Genetics, 123: 411-420
Liu L.L., Yan X.Y., Jiang L., Zhang W.W., Wang M.Q., Zhou S.R., Shen Y., Shen Y.Y., Liu S.J., Chen L.M., Wang J.K., and Wan J.M., 2008, Identification of stably expressed quantitative trait loci for cooked rice elongation in non-Basmati varieties, Genome, 51: 104-112
Marzougui S.,Sugimoto K.,Yamanouchi U.,Shimono M.,Hoshino T., Hori K., Kobayashi M.,Ishiyama K., and Yano M., 2012, Mapping and characterization of seed dormancy QTLsusing chromosome segment substitution lines in rice,Theoretical and Applied Genetics,124: 893-902
McCouch S.R., Cho Y.G., Yano M., Paul E., Blinstrub M., Morishima H., and Kinoshita T., 1997, Report on QTL nomenclature, Rice Genetics Newsletters, 14: 11-13
McCouch S.R., Teytelman L., Xu Y.B., Lobos K.B., Clare K., Walton M., Fu B.Y., Maghirang R., Li Z.K., Xing Y.Z., Zhang Q.F., Kono I., Yano M., Fjellstrom R., DeClerrck G., Schneider D., Cartinhour S., Ware D., and Stein L., 2002, Development and mapping of 2240 new SSR markers for rice (Oryza sativa L.), DNA Research, 9: 199-207
Mei H.W., Xu J.L., Li Z.K., Yu X.Q., Guo L.B., Wang Y.P., Ying C.S., and Luo L.J., 2006, QTLs influencing panicle size detected in two reciprocal introgressive line (IL) populations in rice (Oryza sativa L.), Theoretical and Applied Genetics, 112: 648-656
Shen S.Q., Zhuang J.Y., Wang S.Z., Su Q.Y., Bao J.S., and Xia Y.W., 2005, Analysis on the QTLs with main, epistasis and genotype-environmental interaction effects for cooked rice elongation, Chinese Journal of Rice Science, 19: 319-322 (in chinese)
Sood B.C., Siddiq E.A., and Zaman F.U., 1979, The mechanism of kernel elongation in rice, Indian Journal of Genetics and Plant Breeding, 39: 457-460
Sood B.C., Siddiq E.A., and Zaman F.U., 1983, Genetic analysis of kernel elongation in rice, Indian Journal of Genetics and Plant Breeding, 43: 40-43
Takai T., Nonoue Y., Yamamoto S., Yamanouchi U., Matsubara K., Liang Z.W., Lin H.X., Ono N., Uga Y., and Yano M., 2007, Development of chromosome segment substitution lines derived from backcross between Indica donor cultivar ‘Nona bokra’ and Japonica recipient cultivar ‘Koshihikari’, Breeding Science, 57: 257-261
Tang S., 1987, A study on the cooking and eating quality of Chinese hybrid rice (Oryza sativa L.), Scientia Agriculture Sinica, 20: 17-22
Tian R., Jiang G.H., Shen L.H., Wang L.Q., and He Y.Q., 2005, Mapping quantitative trait loci underlying the cooking and eating quality of rice using a DH population, Molecular Breeding, 15: 117-124
Ujiie K., Kashiwagi T., and Ishimaru K., 2012, Identification and functional analysis of alleles for productivity in two sets of chromosome segment substitution lines of rice, Euphytica, 187: 325-337
Unnevehr L.J., Duff B., and Juliano B.O., 1992, Consumer demand for rice grain quality. International Rice Research Institute, Manila, The Philippines and International Development Research Centre, Ottawa
Wang L.Q., Liu W.J., Xu Y., He Y.Q., Luo L.J., Xing Y.Z., XuC.G., and Zhang Q.F., 2007, Genetic basis of 17 traits and viscosity parameters characterizing the eating and cooking quality of rice grain, Theoretical and Applied Genetics, 115: 463-476
Wissuwa M., Wegner J., Ae N., and Yano M., 2002, Substitution mapping of Pupl: a major QTL increasing phosphorus uptake of rice from a phosphorus-deficient soil, Theoretical and Applied Genetics,105: 890-897
Xiao Y.H., Liu L.L., Jiang L., Lu C.G., Yu C.Y., Zhang W.W., Zuo J.S., Zhai H.Q., and Wan J.M., 2005, Development of the chromosome segment substitution lines (CSSLs) derived from a hybrid rice cross, Pai64S/93-11 with super high yield potential, Rice Genetics Newsletters, 22: 17-20
Xi Z.Y., He F.H., Zeng R.Z., Zhang Z.M., Ding X.H., Li W.T., and Zhang G.Q., 2006. Development of a wide population of chromosome single-segment substitution lines in the genetic background of an elite cultivar of rice (Oryza sativa L.), Genome, 49: 476-484
Yano M., 2001,Naturally occurring allelic variations as a new resource for functional genomics in rice, Rice genetics IV (Proceeding of the fourth international rice genetics symposium 2000), PP.227-238
Ye G.Y.,Liang S.S., and Wan J.M., 2010, QTL mapping of protein content in rice using single chromosome segment substitution lines,Theoretical and Applied Genetics, 121: 741-750
Young N.D., and Tanksley S.D., 1989, Restriction fragment length polymorphism maps and the concept of graphical genotypes. Theoretical and Applied Genetics, 77: 95-101
Zhang G.H., Zeng D.L., Guo L.B., Qian Q., Zhang G.P., Teng S., and Zhu L.H., 2004, Genetic dissection of cooked rice elongation in rice (Oryza sativaL.), Hereditas, 26: 887-892
Zhu W.Y., Lin J., Yang D.W., Zhao L., Zhang Y.D., Zhu Z., Chen T., and Wang C.L., 2009, Development of chromosome segment substitution lines derived from backcross between two sequenced rice cultivars, Indica recipient 93-11 and Japonica donor Nipponbare. Plant Molecular Biology Reporter, 27: 126-131