Introduction
Rice (
Oryza sativa L.) is not only the most important food crop and staple food for 50 percent of the world population, but also it is a model for monocotyledonous plants. Phenotypic traits, especially quantitative ones, are among the most important and least understood traits. Genetic variability of the simple sequence repeats (SSRs) markers is valuable for understanding the phenotypic traits, including fingerprinting genotypes, analyzing genetic diversity, determining variety identity, marker-assisted breeding, phylogenetic analysis, and map-based cloning of genes (
Shen et al., 2004), as well as key component of any breeding program for broadening the gene pool of rice. Expansion of a (GA) dinucleotide at a microsatellite locus associated with domestication of rice (
Ramakrishna et al., 1998). So far, 2740 SSR markers have been genetically mapped in rice, which amounts to one SSR marker every 157 kb (
McCouch et al., 2002). Microsatellites are associated with several genes that control agronomic traits. For example, the allele of RM5 marker on chromosome 1 in
O. rufipogon was associated with an 18% increase in grain yield per plant (
Xiao et al., 1996), and RM5303 was linked with a low silicon rice (
Lsi1) gene (
Ma et al., 2006). Historically, SSRs were used as genetic markers but were functionally unimportant. However, recent research has shown that SSRs have many important functions in terms of development, gene regulation, and evolution (
Lawson and Zhang, 2006). SSRs within genes can be subjected to stronger selective pressure than other genomic regions because of their functional importance (
Li et al., 2004). Functional markers are derived from polymorphic sites with gene involved in phenotypic trait variation (
Andersen and Lübberstedt, 2003). The distribution of SSRs appears highly nonrandom and varies a great deal in different regions of the genes in the genomes (
Lawson and Zhang, 2006). To date, there exist 16 genes cloned by Chinese scientists and 33 genes fine-mapped ones (
Jiang et al., 2007).
Ghd7, encoding a CCT domain protein, isolated from an elite rice hybrid, has major effects on an array of traits in rice, including number of grains per panicle, plant height, and heading date (
Xue et al., 2008).
Yunnan is not only one of the largest centers of genetic ecological diversity in the world (
Zeng et al., 2003;
Zhang et al., 2006) and is the center of origin of Asian
O. sativa (
Zeng et al., 2007), but also the genetic differentiation center of
indica and
japonica subspecies of Asian cultivated rice (
Zeng et al., 2007). It is located in the center of the genetic diversity of Asian cultivated rice—an area that encompasses southwest China, Nepal, Bhutan, Northeast India, Myanmar, Laos, and Northern Thailand (
Chang, 1976;
Gao et al., 2002;
Zeng et al., 2003). In the present study, we examined the allele size of 20 SSR markers and the content of 8 mineral elements and 23 phenotypic traits using 629 accessions of core collection for rice landraces from Yunnan, Southwest China. Our objectives were to reveal the allele size of SSRs associated with mineral elemental contents in brown rice and phenotypic traits, especially in grain traits, aid to QTL mapping and provide useful information for study on possible SSR functions, as well as the reference evidence of markers associated for seeking the relationship between microsatellites and phenotypic traits in rice.
Materials and methods
Plant materials
The 629 accessions from core collection of Yunnan rice landraces came from 100 counties in 16 prefectures in Yunnan Province. The collection was established in our laboratory. These accessions represented over 90% phenotypic diversity of 6121 accessions in Yunnan rice landraces (
Zeng et al., 2007).
Morphological traits
Field experiment was conducted on a farm in Xinping (24ºN, 102ºE, elevation of 500 m above sea level) county of Yunnan Province. All materials were sown in midseason (April-August) of double-cropping rice zone characterized by high temperature. Five-row plots were planted with 20 plants per row with a spacing of 10 cm × 20 cm. Twenty-three phenotypic traits including six plant traits (plant height, tillers per plant, panicles per plant, days to heading, flagleaf length and width), nine grain traits (1000-grain weight, grain length and width, grain length/width ratio, grain thickness, rice length and width, rice thickness, and shattering), eight panicle traits (panicle length, 1-2 internode length, awn length, filled grains per panicle, blighted grains per panicle, total grains per panicle, seed setting rate, and grain density) were examined at appropriate stages during plant growth. Days to heading for each accession was recorded when the first panicle emerged the sheath of flagleaf, and other 22 phenotypic values of each accession were represented by the mean values from 10 individual observations according to ‘Characters to be studied and standard for rating of rice genetic resource’ (Institute of Crop Germplasm Research of Chinese Academy of Agricultural Science, 1986, unpublished literature).
Mineral element content
The experimental plots were established on red mudstone lowland rice soil. Soil analysis showed average pH 5.53, organic matter 2.96 %, available phosphorus 10.08 mg·kg-1, available nitrogen 108.47 mg·kg-1, available potassium 61.42 mg·kg-1, available calcium 1505 mg·kg-1, and available magnesium 310 mg·kg-1. Total P, N, and K content were 0.053%, 0.176%, and 1.99%, respectively. Average available micronutrient level of Cu, Zn, Fe, and Mg in the soil was 3.80 mg·kg-1, 1.45 mg·kg-1, 297 mg·kg-1, and 22.75 mg·kg-1, respectively.
The concentrations of P, K, Ca, Mg, Fe, Zn, Cu, and Mg in 629 rice seeds harvested from the experimental plots were tested with two replications in Supervision and Testing Center for Farm Products Quality, Ministry of Agriculture of the People’s Republic of China. A 0.5 g sample of each brown rice was precisely weighed and put into a beaker with 5 mL of nitric acid and 1 mL of perchloric acid, which was heated to nitrify and decompose rice sample until the solution became clear. The clear solution was continuously heated and evaporated to dryness. The residue was dissolved with 5 mL of 1 mol·L
-1 hydrochloric acid and transferred to a 50-mL graduated bottle. The sample was reduced to ashes at high temperature or decomposed with warm acid decomposition. The residue was dissolved in acid solution in order to transform the elements into inorganic ion. This solution was diluted and analyzed with an ICP-AES spectrometer (
Zeng et al., 2008)
Microsatellites
Twenty SSR markers covering all 12 chromosomes were analyzed on the basis of the published rice microsatellite framework map. The repeat motifs, primer sequences, and chromosomal position for these markers could be found in the RiceGenes database (http://www.gramene.org/microsat/RMprimers.html). Microsatellite primer pairs were obtained from Seagon Bioengineer Limited Company, China. The 30-day seedlings of each accession were used to extract the genomic DNA with a modified 1% CTAB method (
Doyle and Doyle, 1990). The PCR amplifying procedure was that of Panaud et al. (
1996) with slight modification and subsequently run on 8% denatured polyacrylamide gel at 70 W. For the same marker, all runs after the first run included not only the samples but also the checks and a standard molecular weight marker-PUC19 DNA digested by
MspI. When all the samples were run completely, checks were run with another standard molecular weight marker 10 bp DNA ladder from Invitrogen, and the molecular weight for each allele was estimated. All gels were stained with silver method (
Bassam et al., 1991). In the case of null alleles in these species, PCR amplification was repeated to exclude failed PCR reaction. More than two alleles per locus occurred for some accessions. In such cases, we amplified and run them again and selected the alleles that stably occurred in the two replications.
Statistical methods
The mean of two replications for each mineral element concentration in brown rice and each phenotypic value per accession based on 10 individual observations were calculated. Given that some accession showed its heterozygosity at a certain SSR locus, the molecular weight for that SSR marker in that accession was represented by the mean of two allele size. Correlation was determined by applying Pearson’s method; Pearson correlation coefficient was measured between the molecular weight of 20 SSR markers and eight mineral elements, as well as 23 measured phenotypic traits of 629 accessions. Statistical significance was defined as P < 0.05. Statistical analysis was performed on a personal computer using SPSS software. All analyses were performed using SPSS 10.5 software (SPSS, Inc., Chicago, USA).
Results
Number and size of alleles of microsatellites
Allele size of microsatellites and its linked traits are shown in Table 1. There was very high polymorphism distributed over 12 chromosomes. A total of 309 alleles of 20 SSR markers for 629 accessions in Yunnan rice landraces were detected, and the number of alleles per marker was RM257(40) > RM247(28) > RM253(22) > RM241(21) > RM224(18) = RM235 > RM263(17) = RM234= RM223 > RM255(15) > RM258(13) > RM81A(12) = RM225= RM18 > RM5(11) = RM249 > RM244(9) > RM221(8) > RM60(5), ranging from 5 to 40, with an average of 15.5 alleles per locus. This is distinctly higher than 5.8 alleles per locus reported so far (http://www.gramene.org/microsat/microsats.txt). The overall size of PCR products amplified using 20 SSR primer pairs ranged from 95-137 bp (RM235) to 123-210 bp (RM257) (Table 1). The molecular size difference between the smallest and the largest allele for a given SSR locus varied from 8 bp (RM60) to 88 bp (RM253).
Allele size of microsatellites and mineral element content in brown rice
Correlation analysis indicated that there was a significant correlation only in 23 (
r = 0.083
*-0.144
**) of 160 pair traits between the allele size of microsatellites and mineral content (Table 2). P content in brown rice showed a significant correlation with the allele size of RM81A, RM253, RM232, RM234, and RM244, comparatively, indicating that P content was tightly associated with panicle and grain traits. K content was significantly correlated with the allele size of RM235, RM253, RM225, RM5, RM81A, and RM244, indicating that K content was associated with yield traits (
Zeng et al., 2005). Ca content was significantly correlated with the allele size of RM247 and RM255 and thus tightly associated with plant and grains traits. Mg content was significantly correlated with the allele size of RM225, RM244, and RM18, indicating its tight association with yield traits (
Zeng et al., 2005). Fe content was significantly associated with the allele size of RM225, which was linked with nitrogen use efficiency, but Zn content had no correlation with any of the 20 SSR markers. Cu content was significantly associated with the allele size of RM81A, RM253, RM60, and RM247. Mn content was correlated with the allele size of RM60 and RM225.
Allele size of microsatellites and plant traits
Correlation analysis indicated that there was a significant correlation only in 17 (
r = 0.084
*-0.200
**) of 120 pair traits between allele size of microsatellites and plant traits (Table 3). Plant height displayed a significant correlation with the allele size of RM224, RM225, RM60, and RM263. Tillers per plant had a significant correlation with the allele size of RM257, RM18, RM5, RM224, and RM247, among which RM18 and RM224 were also linked with tillers per plant. Panicles per plant displayed significant correlations with the allele size of RM258 related with spikelet fertility (
Cui et al., 2004). Days to heading displayed a significant correlation with the allele size of RM223 (
r = 0.200
**), RM232 (0.143
**), RM224 (-0.113
**), and RM60 (-0.084
*), of which RM60 was also linked with days to heading (
Thomson et al., 2003). Flagleaf length displayed significant correlations with the allele sizes of RM257 and RM224.
Allele size of microsatellites and grain traits
Correlation analysis indicated that there was a significant correlation in 94 (
r = -0.085
*-0.424
**) of 180 pair traits between allele size of microsatellites and grain traits (Table 4). 1000-grain weight displayed a significant correlation with the allele size of six SSR markers ranging from 0.089
* (RM5) to -0.190
** (RM224). RM258 (-0.151
**) was also linked with the 1000-grain weight. There was also a significant correlation in 27 (
r = 0.089
*-0.424
**) of 40 pair traits among allele size of microsatellites and grain size (length and width), among which RM5 and RM211 were also associated with grain size and yield improvement (
Xiao et al., 1996). Grain length/width ratio was one of
indica-japonica classification traits (
Zeng et al., 2000), and it had a significant correlation with the allele size of 14 SSR markers ranging from 0.089
* (RM249) to -0.420
** (RM234). RM5 was linked with grain length/width ratio, and three SSR markers (RM18, RM247 and RM258) were associated with subspecies differentiating DNA markers (Cheng et al., 2004). There was a significant correlation in 30 (
r = -0.087
*-0.349
**) of 60 pair traits among allele size of microsatellites and rice size (length, width and thickness).
Allele size of microsatellites and panicle traits
Correlation analysis indicated that there was a significant correlation in 48 (
r = 0.082
*-0.438
**) of 160 pair traits between the allele size of microsatellites and panicle traits (Table 5). The 1-2 internode length, as one of
indica-japonica classification traits (
Zeng et al., 2000), displayed a significant correlation with the allele size of 14 SSR markers ranging from 0.094
* (RM249) to -0.438
** (RM234), and three SSR markers (RM18, RM247, and RM258) associated with subspecies differentiating DNA markers (Cheng et al., 2004). There was a significant correlation in 16 (
r = -0.082
*-0.349
**) of 80 pair traits among the allele size of microsatellites, seed setting rate, and panicle composition traits (filled grains per panicle, blighted grains per panicle, and total grains per panicle), among which four seed setting traits were significantly correlated with the allele size of RM225 which in turn was linked with nitrogen use efficiency. Seed setting rate was significantly correlated with the allele size of RM257 associated with cold tolerance.
Discussion
Size of alleles for SSR markers associated with elemental contents
The overall size of PCR products amplified using 20 markers for 629 accessions in Yunnan rice landraces ranged from 95-137 bp (RM235) to 123-210 bp (RM257) (Table 1). Mineral elements played an important role in forming genetic diversity (
Zeng et al., 2006). Twenty three pairs (
r = 0.083
*-0.200
**) of 160 pairs between the allele size and mineral element contents, such as P, Cu, and K content in brown rice displayed a significant correlation with the allele size (molecular weight) of RM81A and RM253, and meanwhile, P content was associated with allele size of other markers (RM232, RM244, and RM234), and K content associated with allele size of other markers (RM5, RM255, RM244, and RM235). Recently, there are some evidences for QTL linked to the three SSR markers (RM81A, RM253, and RM244) associated with morphological traits of plants. RM81A is linked with grain yield per plant, 1000-grain weight, and plant height (
Cui et al., 2004), RM253 is associated with percent head rice, amylase content, alkali spreading score, protein content, grain length, heading date, awn length, days to maturity and grain length/width ratio (
Xiao et al., 1996;
Linh et al., 2006), with RM244 associated with fertility restorer genes (Jing et al., 2000), RM232 associated with 1000-grain weight, panicles per plant, panicle length, grains per panicle, and grain yield per plant (
Tan et al., 2007), and RM234 associated with days to heading, days to maturity and maximum root length (
Thomson et al., 2003;
Miyata et al., 2007). There existed significant correlations (
P < 0.01) between P content in brown rice and 11 traits (protein content, rice width, panicle length, filled grains per panicle, blighted grains per panicle, seed setting rate, plant height, tillers per plant, days to heading, grain width, and thickness) and between K content in brown rice and 8 traits (amylose contents, alkali spreading score, gel consistency, protein content, plant height, tillers per plant, grain thickness, and panicle length) (
Zeng et al., 2005). The genetic diversity investigated based on eight mineral contents in brown rice from Lincang, Simao, Xishuangbanna, and Dehong Prefectures has the similar trend with the genetic diversity of rice landraces of Yunnan revealed based on morphological traits and isozyme and SSR markers (
Zeng et al., 2006). Therefore, we suggest that these eight elemental concentrations of P, K, Ca, Mg, Fe, Zn, Cu, and Mn in brown rice are associated with the allele size of rice landraces in Yunnan Province, and further understanding the relationship of mineral elements associated with origin of life is needed. As far as we know, this is the first report on genotypic variations in the eight elemental concentrations in brown rice associated with the allele size of SSR markers.
Allele size of SSR markers associated with phenotypic traits
Seventeen pairs (
r = 0.083
*-0.144
**) of 120 pairs between the allele size and plant traits, such as some markers associated with plant height (RM263, RM60, RM225 and RM224) and days to heading (RM60, RM232, RM223 and RM224), were associated with the allele size of 4 SSR markers, respectively. Tillers per plant displayed a significant correlation with the allele size (molecular weight) of 5 SSR markers (RM5, RM18, RM257, RM224 and RM247). Panicles per plant displayed a significant correlation with the allele size (molecular weight) of RM249 and RM258. In the previous study, RM224 had been found associated with plant height, tillers per plant, panicle length and root-shoot ratio (
Miyata et al., 2007), with RM60 associated with high-tillering dwarf 1, days to heading and shattering (
Septiningsih et al., 2003;
Thomson et al., 2003). Tillers per plant also displayed a significant correlation with the allele size of 5 SSR markers (RM5, RM18, RM257, RM224 and RM247) and RM249 was associated with tillers per plant, panicles per plant and grains number per plant (
Abdelkhalik et al., 2005;
Pradeep et al., 2005;
Miyata et al., 2007). Furthermore 94 pairs (
r = -0.085
*-0.424
**) of 180 pairs between allele sizes and grain traits, such as nine grain traits displayed a significant correlation with the allele size (molecular weight) of RM5 and RM224. Previous studies showed that RM5 was associated with grains per panicle, percent germination, head rice grain and grain length/width ratio, especially the yield-enhancing genes increased the yield of hybrid rice by 18% (
Xiao et al., 1996;
Thomson et al., 2003;
Liang et al., 2004). Eight grain traits displayed a significant correlation with the allele size (molecular weight) of RM257 and RM244, with RM257 associated with cold resistance and heading date (
Miyata et al., 2007). Seven grain traits were significantly correlated with the allele size (molecular weight) of RM255 and RM235, as well as RM247, with RM255 associated with panicle size, panicles per plant, panicle neck diameter, and 1000-grain weight (
Xu et al., 2004) and RM235 associated with panicles per plant (
Thomson et al., 2003). There were six grain traits significantly correlated with the allele size of four SSR markers (RM81A and RM263, RM225, and RM258). Previous studies also showed that RM263 was associated with broken rice grain, crushed rice grain, 1000-grain weight, and spotted leaf (
Septiningsih et al., 2003;
Xu et al., 2004), and RM258 was associated with the spikelet fertility, 1000-grain weight, tiller number, and heading date (
Cui et al., 2004;
Miyata et al., 2007). Five grain traits displayed a significant correlation with the allele size of RM211 and RM18. In addition, RM211 was related with yield-improving, grain width and days to heading (
Miyata et al., 2007), while RM18 was associated with frizzle panicle, grain yield, tillers per plant, grains per panicle, and heading date (
Duan et al., 2003). Grain length/width ratio and grain width were significantly correlated with the allele size of RM241 etc., which proved to be associated with panicles per plant as well as plant height (
Thomson et al., 2003;
Cui et al., 2004), with 14 SSR markers ranging from 0.089
* to -0.438
**, and three SSR markers (RM18, RM247, RM258) were associated with subspecies differentiating DNA markers (Cheng et al., 2004). Grain length/width ratio and 1-2 internode length are two of
indica-japonica classification traits (
Zeng et al., 2000). 48 pairs (
r = 0.082
*-0.438
**) of 160 between allele size and panicle traits, such as four panicle traits, displayed a significant correlation with the allele size of 5 SSR markers (RM211, RM225, RM223, RM258, and RM224), among which RM225 proved to be associated with nitrogen use efficiency and grain breadth (http://www.cropscience.org.au/icsc2004/poster/3/4/1/1296_senthilvels.htm) and RM223 associated with panicle per plant, days to heading and panicle length (
Pradeep et al., 2005;
Miyata et al., 2007). Three panicle traits displayed a significant correlation with the allele size of four SSR markers (RM81A, RM257, RM235 and RM247), among which RM247 proved to be associated with phosphorus deficiency and plant height (
Thomson et al., 2003). In general, allele size (molecular weight) of microsatellites is more associated with grain or panicle traits than plant traits or mineral element content in brown rice.
Correlation analysis indicated that there existed a significant correlation not only in 182 of 620 pairs between the allele size of 20 SSR markers and 23 phenotypic traits, not only in eight mineral elements in brown rice but also in 17 pairs between allele size of RM224 and 23 phenotypic traits and eight mineral elements in brown rice, with RM224 associated with plant height, tillers per plant, panicle length, and root–shoot ratio (
Miyata et al., 2007). Up to now, 8558 QTLs controlling various complex traits have been located on different chromosome regions in
Oryza sativa (http://www.gramene.org/qtl/index.html). However, those published have not clearly discussed the allele size (molecular weight) of microsatellites and phenotypic traits, because of the limited available information about SSR locations on chromosomes. SSRs in different positions of a gene can play an important role in determining protein function, genetic development, and regulation of gene expression (
Lawson and Zhang, 2006). These 20 SSR markers for 629 accessions in Yunnan rice landraces were linked to phenotypic traits based on previous study. SSR allele size of core collection associated with phenotypic traits in rice has not been reported so far. Therefore, the data on combined phenotypic traits and allele size (molecular weight) of microsatellites within genes may provide strong evidence on the allele size of SSRs associated with mineral elemental contents in brown rice, and the SSRs on gene regulation and phenotypic traits, especially in grain traits, so as to further seek for the relationship between the molecular weight of microsatellites and phenotypic traits in rice based on markers associated.
Higher Education Press and Springer-Verlag Berlin Heidelberg