Skip to main content

Detection of Novel QTLs Regulating Grain Size in Extra-Large Grain Rice (Oryza sativa L.) Lines

Abstract

Background

Grain size is an important trait that affects rice yield. Although many genes that contribute to grain size have been cloned from mutants or by quantitative trait locus (QTL) analysis based on bi-parental mapping, the molecular mechanisms underlying grain-size determination remain poorly understood. In this study, we identified the lines with the largest grain size and detected novel QTLs affecting the grain size.

Results

We screened the National Institute for Agrobiological Sciences Genebank database and identified two rice lines, BG23 with the widest grain and LG10 with the longest grain. Using these two lines, we performed QTL analysis for grain size. Eight QTLs were detected during the QTL analyses using F2 populations derived from crosses between the large-grain lines BG23 or LG10 and the middle-size grain cultivars Nipponbare and Kasalath. Both BG23 and LG10 possessed large-grain alleles of four major QTLs: GW2, GS3, qSW5/GW5, and GW8. Other three minor QTLs were derived from BG23. However, these QTLs did not explain the differences in grain size between these two lines. Additionally, four QTLs for grain length or width were detected in an F2 population derived from a cross between BG23 and LG10; this population lacked the strong effects of the four major QTLs shared by both parent plants. Of these newly detected QTLs, the effects of two QTLs, GL3b and GL6, were confirmed by progeny testing. Comparison of the length of inner epidermal cells in plants homozygous for BG23 and LG10 alleles indicated that GL3b and GL6 genes regulate cell elongation and cell division, respectively.

Conclusions

In this study, we detected 12 loci including 14 QTLs regulating grain size from two lines with largest grains available in Japanese stock. Of these loci, we confirmed the effect of two gene loci and mapped their candidate region. Identification of novel genes regulating grain size will contribute to our understanding of the molecular mechanisms controlling grain size.

Background

Grain size is an important trait in breeding, serving as a factor determining rice yield. Generally, genes regulating grain size are classified into two groups, those controlling cell elongation and those controlling cell division. A short-grain phenotype due to shortened cell length in the loss-of-function mutants of brassinosteroid (BR)-related gene D61 implies that these genes control cell elongation (Segami et al. 2012). Furthermore, SRS3 and SRS5 are also involved in cell elongation. SRS3 and SRS5 encode a kinesin-13 protein, and an α-tubulin, respectively (Kitagawa et al. 2010, Segami et al. 2012). These srs mutants commonly show short-grain phenotypes due to shortened cell lengths, similar to the BR-related mutants. However, the Srs5 mutant is seemingly not involved in BR signal transduction. D1 and TUD1 are involved in cell division (Izawa et al. 2010; Hu et al. 2013). The mutants d1 and tud1 commonly show a short-grain and dwarf phenotype owing to reduced cell numbers. D1 and TUD1 encode the α-subunit of a heterotrimeric G-protein and U-box ubiquitin ligase, respectively, and d1 is epistatic to the tud1 mutation (Ashikari et al. 1999; Fujisawa et al. 1999; Hu et al. 2013). Of these mutants, we succeeded in classifying srs3, Srs5, d1, and d61 into cell elongation- or division-type mutants by comparing the cell length and cell number of the inner-epidermal cells of lemma (Izawa et al. 2010; Kitagawa et al. 2010; Segami et al. 2012). Quantitative trait locus (QTL) analysis also revealed other genes involved in grain size regulation. GS3 controls grain length and weight (Fan et al. 2006; Mao et al. 2010; Takano-Kai et al. 2009); GW2 and qSW5/GW5 control grain width and weight (Shomura et al. 2008; Song et al. 2007; Weng et al. 2008); GS5 controls grain width, filling, and weight (Li et al. 2011); GW8 controls grain width (Wang et al. 2012); GL3.1 (also named OsPPKL1) controls grain length (Qi et al. 2012; Zhang et al. 2012); TGW6 and GW6a control grain weight (Ishimaru et al. 2013; Song et al. 2015); and GL7/GW7 controls grain length and width (Wang et al. 2015a, b). Interestingly, among these genes cloned by QTL analysis, only GL7/GW7 was involved in cell elongation, whereas the others were involved in cell division. It is noteworthy that these genes encode a variety of different proteins involved in various signal transduction pathways. GS3 encodes the γ-subunit 3 of a heterotrimeric G-protein (Li et al. 2012); GW2 encodes a RING-type domain with E3 ubiquitin ligase activity (Song et al. 2007); qSW5/GW5 encodes a novel nucleoprotein that interacts with polyubiquitin in a yeast two-hybrid experiment (Weng et al. 2008); GS5 encodes a putative serine carboxypeptidase and acts as a positive regulatory factor in the cell cycle (Li et al. 2011); GW8 encodes the transcription factor Squamosa promoter-binding protein-like 16 containing the miR156-targeted site (Wang et al. 2012); GL3.1 encodes a Ser/Thr phosphatase of the protein phosphatase of the kelch-like family (Qi et al. 2012; Zhang et al. 2012); TGW6 encodes a protein with indole-3-acetic acid-glucose hydrolase activity (Ishimaru et al. 2013); GW6a encodes a functional GCN5-related N-acetyl-transferase-like protein that harbors intrinsic histone acetyl-transferase activity (Song et al. 2015); and GL7/GW7 encodes a homolog of the Arabidopsis thaliana TONNEAU1 recruiting motif (TRM) protein (Wang et al. 2015a, b). Although many genes were cloned and identified, the molecular mechanisms underlying grain size determination are not well understood, because most of these genes are not involved in the same signal transduction pathways. To clarify these mechanisms, it is necessary to identify and classify additional genes that contribute to grain size, whether through cell elongation or cell division, by unified methodology. Therefore, the aim of this study was to identify, from the National Institute for Agrobiological Sciences (NIAS) Genebank database (http://www.gene.affrc.go.jp), the rice lines with the largest grain, and using these lines, detect novel QTLs and elucidate how these QTLs control grain size.

Results

BG23 and LG10 Display the Extra-Large Grain Phenotype

The parental lines used in QTL analysis for grain size were selected by screening for the largest-grain rice lines in the public database NIAS Genebank. The grain width and length of screened lines ranged from 1.7 to 4.6 mm and 3.4 to 14.3 mm, respectively (Additional file 1: Figure S1A and B). Of these screened lines in NIAS Genebank, BG23 (grain width 4.6 mm) and LG10 (grain length 14.3 mm), which produced seeds with the largest width and length, respectively, were selected as parental lines for QTL analysis. Because the grain length of Nipponbare and Kasalath was 6.97 ± 0.28 mm and 7.99 ± 0.33 mm, respectively, and near to the median value of the NIAS Genebank database, we used these lines as middle-size lines for grain length and width comparison (Fig. 1, Table 1). The grain length of BG23 and LG10 was 11.24 ± 0.42 mm and 14.25 ± 0.46 mm, respectively (Fig. 1, Table 1). The grain width of BG23 and LG10 was 4.43 ± 0.22 mm and 4.06 ± 0.36 mm, respectively, whereas that of Nipponbare and Kasalath was 3.27 ± 0.27 mm and 2.55 ± 0.21 mm, respectively (Fig. 1, Table 1). Because the grains of BG23 and LG10 were significantly larger compared to those of the middle-size grain cultivars, these lines were expected to possess different composition of genes regulating grain size. Clarifying the compositions and effects of each gene will provide information about epistasis and effect of interactions of these genes. Therefore, this analysis was expected to provide novel and useful genes with the potential to increase grain yield by controlling grain size.

Fig. 1
figure 1

Grain size and morphology of middle-size grain cultivars Nipponbare and Kasalath and extra-large grain lines BG23 and LG10. Bar = 5 mm

Table 1 Comparisons of grain length and width in rice lines Nipponbare, Kasalath, BG23, and LG10

BG23 and LG10 Shared Four Major QTLs for Grain Size

To identify the genes regulating grain length and width in BG23 and LG10, four combinations of F2 populations were derived from the following crosses: Kasalath × BG23, Kasalath × LG10, Nipponbare × BG23, and Nipponbare × LG10. A plot of grain length and width of the F2 individuals in each population showed a continuous distribution (Fig. 2a–h), indicating that the differences in grain traits between the extra-large grain lines and middle-size grain lines were regulated by QTLs. QTL analysis for grain length and width detected two loci on Chr2 and Chr3 that were shared by all four F2 populations (Fig. 2i–l, Table 2, Additional file 2: Figure S2A–H) and one locus on Chr5 shared by two F2 populations derived from Kasalath crosses (Fig. 2i and j, Table 2). Unshared QTLs were detected on Chr1, Chr4, Chr10, and Chr8 (Fig. 2i–k, Table 2). In total, we detected eight QTLs including seven loci in this analysis (Fig. 2i–k, Table 2).

Fig. 2
figure 2

Frequency distribution of grain length and width and quantitative trait loci (QTLs) detected in F2 populations. Frequency distribution of grain length is shown in a, c, e, and g. Frequency distribution of grain width is shown in b, d, f, and h. F2 populations derived from Kasalath × BG23 (a and b), Kasalath × LG10 (c and d), Nipponbare × BG23 (e and f), and Nipponbare × LG10 (g and h) crosses. Vertical white and black arrowheads indicate the means of the middle-size grain cultivars and extra-large grain lines, respectively. Genetic linkage maps of F2 populations derived from Kasalath × BG23 (i), Kasalath × LG10 (j), Nipponbare × BG23 (k), and Nipponbare × LG10 (l). Boxes indicate positions of regions harboring QTLs with the logarithm of the odds > 3.0 for grain length (black) and width (white). Horizontal black arrowheads indicate the position of previously cloned QTLs for grain size. Double-headed arrows indicate the position of previously detected QTLs for grain size. Bars indicate marker positions

Table 2 Quantitative trait loci for grain length and width detected in F2 populations derived from the crosses Kasalath × BG23, Kasalath × LG10, Nipponbare × BG23, and Nipponbare × LG10

As several QTLs were detected around previously cloned gene loci, we performed marker analysis for the genes GW2, GS3, qSW5/GW5, and GW8. Both BG23 and LG10 possessed all of the alleles responsible for large grain size in all four genes (Additional file 3: Figure S3A–H). The QTLs for grain width and length, detected on Chr2 and Chr3, respectively, in all tested populations correspond to GW2 and GS3, whereas the QTLs for grain width detected on Chr5 and Chr8 correspond to qSW5/GW5 and GW8.

The QTLs for grain length on Chr4 and Chr10 detected from the Kasalath × LG10 cross and Chr1 detected from the Nipponbare × BG23 cross include qGL1, qGL4b, and qGL10; loci have been previously reported (Kato et al. 2011; Wang et al. 2012; Yan et al. 2014) (Fig. 2j and k). In total, we detected seven loci including four cloned and three non-cloned QTLs for grain size. However, as GW8 was detected only in the F2 population derived from the Kasalath × BG23 cross, these minor QTLs appeared to be masked by the effects of the major QTLs. Here, we examined whether the QTL analysis for grain length using F2 population derived from BG23 and LG10 cross is possible to detect the same loci on Chr1, Chr4, and Chr10 detected in other combinations of F2 populations and whether these loci can fully explain the differences between BG23 and LG10. Hence, the QTLs for grain length and width were investigated using the F2 population derived from the BG23 × LG10 cross, thus avoiding the strong effects of the four QTLs shared by the parents.

Novel QTLs Detected from the BG23 × LG10 Cross

Screening for simple sequence repeat (SSR) markers for QTL analysis that could indicate the differences in genotypes between BG23 and LG10 was unsuccessful, because these two parental lines exhibited a highly similar banding pattern on the gel. To construct the linkage map, we performed re-sequencing of the parents and designed BeadArray markers using detected single nucleotide polymorphisms (SNPs).

The histograms of grain length and width in F2 individuals showed a continuous distribution, indicating that the differences between BG23 and LG10 with respect to these parameters were likely regulated by multiple QTLs (Fig. 3a and b). QTL analysis of the F2 population from the BG23 × LG10 cross revealed four QTLs for grain length on Chr1, Chr3, the short arm of Chr6 (Chr6S), and the long arm of Chr6 (Chr6L) and two QTLs for grain width on Chr1 and Chr6S (Fig. 3c, Table 3, and Additional file 3: Figure S3I and J). Additionally, the QTL for grain length on Chr6L included GW6a and GW6b (Song et al. 2015) (Fig. 3c). Overall, we detected six QTLs including five loci (Fig. 3c, Table 3). The QTL for grain width on Chr1 and that for grain length on Chr3 and Chr6S represent novel loci regulating grain size. A two-dimensional scan revealed no significant epistatic interactions among these QTLs in this population (Additional file 4: Figure S4).

Fig. 3
figure 3

Frequency distributions of grain length and width and quantitative trait loci (QTLs) detected in F2 populations derived from the BG23 × LG10 cross. Frequency distributions of grain length (a) and width (b). Vertical white and black arrowheads indicate the means for BG23 and LG10, respectively. (c) Boxes indicate the positions of regions harboring QTLs with the logarithm of the odds > 3.0 for grain length (black) and width (white). Horizontal black arrowheads indicate the positions of previously cloned QTLs for grain size. Double-headed arrows indicate the positions of previously detected QTLs for grain size. Bars indicate marker position

Table 3 Quantitative trait loci for grain length and width detected in F2 population derived from the BG23 × LG10 cross

Phenotypic Evaluation of Novel QTLs

To confirm the effect of detected QTLs on grain length, we screened F3 lines including one segregating locus and three other fixed loci. Of the four loci, we screened two segregation lines for QTLs detected on Chr3 and Chr6L (Additional file 5: Figure S5) and from these segregation lines selected plants homozygous for BG23 and LG10 alleles on Chr3 (Chr3-BG23 and Chr3-LG10) and Chr6L (Chr6L-BG23 and Chr6L-LG10). The grain length of Chr3-LG10 and Chr6L-LG10 were significantly longer compared to those of Chr3-BG23 and Chr6L-BG23, respectively (Fig. 4a, c). These results clearly indicate that these two loci include genes responsible for grain length. LG10 alleles on Chr3 and Chr6L have a positive effect on grain length. We named these two genes GL3b and GL6. To determine whether these grain-size regulating genes are involved in cell division or cell elongation, we compared the length of the inner epidermal cell in lemma by using scanning electron microscopy (SEM) (Fig. 4b). The length of the inner epidermal cells was significantly longer in Chr3-LG10 than in Chr3-BG23, whereas no significant difference was observed between Chr6L-BG23 and Chr6L-LG10 (Fig. 4d). In contrast, the estimated number of cells in Chr6L-LG10 was significantly higher than that in Chr6L-BG23 (Fig. 4e). These results indicate that GL3b and GL6 affect cell elongation and cell division, respectively.

Fig. 4
figure 4

Comparison of grain length and inner-epidermal cell length in segregation lines. a Grain morphology of Chr3-BG23, Chr3-LG10, Chr6-BG23, and Chr6-LG10. Bar = 5 mm. b Inner-epidermal cells in the lemma of Chr3-BG23, Chr3-LG10, Chr6-BG23, and Chr6-LG10. Bars = 100 μm. c Grain length of Chr3-BG23, Chr3-LG10, Chr6-BG23, and Chr6-LG10. d Cell length of Chr3-BG23, Chr3-LG10, Chr6-BG23, and Chr6-LG10. e Estimated cell number of Chr3-BG23, Chr3-LG10, Chr6-BG23, and Chr6-LG10. Numbers in c–d indicate the averages. Student’s t-test was used to generate the P values. * indicates statistical significance from the neighboring value (P < 0.05). n.s. indicates no significance (P > 0.05)

Genetic Mapping of GL3b and GL6

To narrow down the candidate region of GL3b and GL6, we performed genetic mapping using recombinant lines of each locus. Recombinant screening of these two loci using F3 population of segregating lines resulted in four and three recombinants in GL3b and GL6 regions, respectively (Fig. 5a and b). Because the F3 plants possess heterozygous chromosomes, we selected F4 plants homozygous for BG23 or LG10 and recombinant chromosome of each GL3b and GL6 region to compare the phenotype of grain length. Using the data on grain length in F5 seeds obtained from selected F4 plants from each line, we performed mapping analysis for GL3b and GL6 gene loci.

Fig. 5
figure 5

Genetic mapping of GL3b and GL6. a Recombinant plants possessing breakpoint at the Chr3 locus were screened from the F3 population of segregation lines by using PCR markers RM7000 and BGLG_3-37.2_Hinf. We found four F3 recombinant plants. Because the F3 plants were heterozygous, we selected F4 plants homozygous for BG23 or LG10 chromosome and recombinant chromosome to compare the phenotype for grain length. At least three homozygous plants were selected for each genotype. The grain length of these homozygous F4 plants was measured by using 50 F5 seeds. Two types of genotype selected among homozygous plants from each line are indicated to the left of the F2 line number. The grain length of each line is indicated to the right of the F2 line number. The results of the t-test for grain length between the plants with or without recombination are indicated to the right of the grain length of F5 grains. GL3b gene locus was mapped into a 1.2-Mb region between BGLG_3-36.1_Hinf and BGLG_3-37.2_Hinf. b Recombinant plants possessing breakpoint at the Chr6L locus were screened from the F3 population of segregation lines by using PCR markers BGLG_6-26.0_Hinf and BGLG_6-28.12_Hae. We found three F3 recombinant plants. The grain length of homozygous F4 plants was measured by using 50 F5 seeds. GL6 gene locus was mapped into a 1.3-Mb region between BGLG_6-26.79_Hinf and BGLG_6-28.12_Hae. White and black boxes indicate BG23 and LG10 chromosomes, respectively. Student’s t-test was used to generate the P values. * indicates statistical significance from the neighboring value (P < 0.05). n.s. indicates no significance (P > 0.05). The numbers under the marker name indicate marker’s position (Mb)

At the GL3b locus, we screened recombinant plants with a recombination breakpoint between the PCR markers RM7000 and BGLG_3-37.2_Hinf and obtained four F3 recombinant plants (Fig. 5a). These plants were grown to obtain F4 seeds, which were then sawn and the resulting plants were screened homozygous for BG23 or LG10 chromosome as well as for recombinant chromosomes and compared the grain length of the resulting F5 seeds. The recombinant plant selected from F4 line no. 1 with a recombination event between 35.20 Mb and 35.35 Mb showed a significantly longer grain compared to that of sibling plants homozygous for BG23, indicating that the causal gene of GL3b is located on the right from 35.20 Mb (Fig. 5a). Similarly, the F4 line no. 7 also indicates that GL3b is located to the right from 35.20 Mb because grain length of the recombinant plant is the same as that in plants homozygous for BG23. Similarly, the right border of GL3b was indicated at 36.37 Mb according to F4 lines 16 and 18 (Fig. 5a). Taken together, we concluded that GL3b locus is located in a region 1.2 Mb long, between 35.20 Mb and 36.37 Mb on Chromosome 3.

At the GL6 locus, we screened three recombinant plants possessing recombination breakpoint between the PCR markers BGLG_6-26.79_Hinf and BGLG_6-28.12_Hae (Fig. 5b). Because the recombinant plant selected from the F4 line no. 1 with a recombination event between 26.79 Mb and 27.09 Mb shows same overlap in grain length with sibling plants homozygous for BG23, we conclude that the causal gene GL6 locates to the right from 26.79 Mb (Fig. 5b). Similarly, the right border of GL6 was indicated as 28.12 Mb by F4 lines 36 and 20 (Fig. 5b). Taken together, we concluded that GL6 locus locates in a 1.3-Mb long region, between 26.79 Mb and 28.12 Mb on Chromosome 6.

Discussion

In this study, we detected eight QTLs on seven loci from four populations derived from crosses between large-grain lines and middle-size grain lines. Furthermore, we newly detected six QTLs on five loci from the population derived from the cross between large grain lines. In total, we detected 12 loci including four cloned QTL genes. Both BG23 and LG10 contained four major grain-size regulating genes: GW2, GS3, qSW5/GW5, and GW8 (Fig. 2i-l). The results presented herein indicate that the additive effect of GW2 locus for grain width, and of GS3 locus for grain length, detected from four combinations of F2 population were distributed from 0.2 to 0.3 and from 0.44 to 0.7, respectively (Fig. 2i-l, Table 2). These additive effects are consistent with previously reported effects of GW2 and GS3 (Fan et al. 2006; Song et al. 2007). The dominant GS3 locus detected in our four populations was distributed from −0.42 to −0.59. This indicates that plants heterozygous for GS3 locus show equal phenotype as plants homozygous for functional GS3 gene. Thus, the GS3 locus detected in our study behaves in a recessive manner, which is consistent with previous report (Fan et al. 2006). The locus qSW5/GW5 was detected in the F2 population derived from Kasalath × BG23 and Kasalath × LG10 crosses (Fig. 2i and j, Table 2). Originally, qSW5 was detected from descendants derived from the Kasalath × Nipponbare cross (Shomura et al. 2008), with the Nipponbare allele of qSW5 being recessive. Our results corroborate this feature of qSW5. GW8 was detected only in the F2 population derived from the cross between BG23 × Kasalath, whereas LG10 also contained large grain allele (Fig. 2i, Additional file 3: Figure S3D). The GW8 locus detected in the BG23 × Kasalath cross showed additive effect of 0.15 and dominance effect of −0.16 (Table 2). This result indicates that the large grain allele of GW8 locus behaves in a recessive manner, although previously GW8 was reported as semi-dominant gene (Wang et al. 2012). Because the effects of GW8 locus for grain size regulation are comparatively weaker compared to the effects of the other three genes and not detected in the BG23× Nipponbare cross, the effects of GW8 seemed to be masked or distorted by other three major QTLs. Therefore, the effects of minor QTLs may not be detectable. Although it is difficult to compare the effects of QTLs detected from different populations because the genetic background affects the effect of each QTL, we assumed that there should be some alternative QTLs explain the differences between BG23 and LG10. Thus, we tried to detect QTLs for grain size between BG23 and LG10 while avoiding the four major QTLs. As we expected, we could detect additional four QTLs in the F2 population derived from the cross between BG23 and LG10. Of these four QTLs, two loci located on Chr3 and Chr6L and named GL3b and GL6 were confirmed to be involved in cell elongation and cell division, respectively (Fig. 4a-e). Observation of the inner epidermal cells in lemma revealed that GL3b promotes grain length by promoting cell elongation (Fig. 4b-e). Using the same method, the genes responsible for cell elongation had been identified from a BR-related mutant (d61) and microtubule-related mutants (srs3 and Srs5) (Kitagawa et al. 2010; Segami et al. 2012). TRM-containing protein (GL7/GW7) cloned in the QTL analysis also regulates cell elongation in rice (Wang et al. 2015b). Because these studies cloned only those genes involved in cell elongation, genetic diversity of the grain size in rice was considered to have developed mainly owing to the differences in the genes responsible for cell division. Therefore, clarifying the molecular functions of GL3b via cloning and genetic analysis of these genes will be important in elucidating the complex genetic mechanism underlying grain size regulation. GL6 is involved in cell division (Fig. 4b-e). Since it was detected in the F2 population derived from the crossing of BG23 and LG10, both of which share the long-grain alleles of GS3 and GW8 (Fig. 3c, Additional file 3: Figure S3), GL6 promotes grain length in an additive manner with GS3 and GW8. This indicates that GL6 belongs to a genetically separated signal transduction pathway from those of GS3 and GW8.

In contrast, the grain of BG23 is wider than the grain of LG10 (Fig. 1, Table 1). Nevertheless, LG10 provided positive additive effects for all detected QTLs for grain width in the F2 population derived from the crossing of BG23 and LG10. This is explained by the grain length of LG10, which is too long in the longitudinal direction to fill grain hull. Therefore, starch input is not sufficient to fill the volume of the grain hull. This starch shortage leads to unfilled and thinner grains than predicted by the genetic potential of LG10.

Conclusion

In this study, we detected 12 loci that are involved in the regulation of grain size, using rice lines with largest grain size in Japanese stock (Figs. 2i-l, 3c). Among these loci, we identified two loci GL3b and GL6 that promote grain length through increased cell length and cell division, respectively (Fig. 4b-e). Although we confirmed the effects of the two loci detected in this study, the effects of six loci remain unknown (Fig. 2i-l, and Fig. 3c). The results indicate that the molecular mechanisms underlying grain size determination are complex. To clarify the mechanisms and genetic architecture, we need to isolate the genes and test the genetic interactions of each gene in future studies.

Methods

Plant Materials

Two rice (Oryza sativa L.) cultivars, Nipponbare (japonica) and Kasalath (indica), and two extra-large grain rice lines, BG23 (japonica) and LG10 (japonica), were used in this study. BG23 (accession number: 8966) and LG10 (accession number: 84410) were obtained from NIAS Genebank after screening for lines with the widest and the longest grain, omitting the data with input errors (containing wrong trait annotations, low reliability of the reads, etc.).

For the QTL analysis, five F2 populations were produced from the following parental crosses: Kasalath × BG23, Kasalath × LG10, Nipponbare × BG23, Nipponbare × LG10, and BG23 × LG10. Phenotypic evaluations for QTL analysis using these F2 populations were performed in the 2012 growing season. Phenotypic evaluations for F2, F3, and F4 populations derived from BG23 × LG10 cross were performed in the 2013, 2014, and 2015 growing seasons, respectively. Rice seeds were sown at the beginning of April, and at the beginning of May the seedlings were transplanted into a paddy field at the Research Center for Bioresources Development, Fukui Prefectural University, Fukui, Japan, and grown under natural conditions.

Phenotypic Evaluation

Harvested rice seeds were air-dried and stored in a seed dryer for at least 2 weeks. Fully filled grains were then evaluated for grain length and width. Fifty grains were randomly selected from each plant and scanned with a scanner (GT8300, EPSON, Nagano, Japan) to obtain images that were then used to measure each size using Smart Grain software (Tanabata et al. 2012).

The inner-epidermal cells of lemmas in mature grain hull were observed by SEM (TM3030, Hitachi, Tokyo, Japan). We measured cell length of the central cell layer of the hull (n > 20).

QTL Analysis

Genomic DNA was extracted from each F2 progeny. The number of individuals and markers used from each F2 population for each genotype were as follows: 94 individuals and 92 markers in Kasalath × BG23 F2 population, 94 individuals and 90 markers in Kasalath × LG10 F2 population, 94 individuals and 66 markers in Nipponbare × BG23 F2 population, 94 individuals and 105 markers in Nipponbare × LG10 F2 population, and 318 individuals and 106 markers in BG23 × LG10 F2 population (Additional file 6: Table S1 and Additional file 7: Table S2). PCR-based markers, including SSR markers (McCouch et al. 2002; Ware et al. 2002; http://www.gramene.org/) and cleaved amplified polymorphic sequence (CAPS) markers (Konieczny and Ausubel 1993), were used to construct linkage maps. An Illumina GoldenGate Assay with VeraCode (San Diego, CA, USA) was used according to the manufacturer’s instructions.

QTL analysis was performed using the R/qtl software (http://www.rqtl.org; Broman et al. 2003). Putative QTLs were detected by using the interval mapping function. Two QTLs of the same chromosome were confirmed by composite interval mapping function. The QTLs were assumed present when the logarithm of the odds score was 3.0 or higher. The additive and dominant effects and the phenotypic variance explained by each QTL at the maximum LOD score were estimated by using the sim.geno, makeqtl, fitqtl, and effectplot functions in the R/qtl software (Broman et al. 2003).

Marker Analysis

GW2-specific CAPS markers were amplified using the primers GW2-Hinf1-F and GW2-Hinf1-R. GS3 marker analysis was carried out according to the protocol described by Fan et al. (2009) using primers SF28-U and SF28-L. qSW5/GW5 marker analysis was conducted according to the method presented in Weng et al. (2008) using primers Indel2-F and Indel2-R, and GW8 marker analysis was carried out using indel markers to detect a 10-bp deletion in the promoter region (Additional file 8: Table S3).

Resequencing Analysis of the Extra-Large Grain Lines

Construction of the shotgun library and sequencing were conducted with a Genome Analyzer-IIx (GA-IIx; Illumina, USA). Rice genomic DNA was fragmented by Covaris (Woburn, MA, USA). The size-selected DNA (average length 300 bp) was purified with an AmpureXP (Beckman Coulter, Inc., Fullerton, CA, USA) and gel extraction. The standard short-read library was built using a TruSeq DNA Sample Prep Kit v2-SetA (Illumina, USA) and TruSeq SBS v5 (Illumina, USA) according to manufacturer’s instructions for sequencing runs at 2 × 100 bp total. After sequencing, GA-IIx real time analysis 1.13.48.0 and CASAVA 1.8.2 (Illumina, USA) were utilized for base calling.

The rice genome IRGSP sequence (IRGSP Build5 Pseudomolecules of the Rice Genome, http://rgp.dna.affrc.go.jp/E/IRGSP/Build5/build5.html) was used for reference mapping. The pair-ended reads were used for reference mapping with BWA ver.0.6.1 (Li and Durbin 2010) and SAMtools ver. 0.1.18 (Li et al. 2009), converting the mapped data to a SAM file and then generating a pileup of the sequence data. SNP, deletions, and insertions were detected using CLC Genomics Workbench v6.5.1 (CLCbio, Qiagen, Aarhus, Denmark). After selecting the nucleotide variations with frequency equal to 100.00, we identified 252,343 SNPs and INDELs between BG23 and Nipponbare and 326,351 SNPs and INDELs between LG10 and Nipponbare. SNPs included in the center of single-copied 60-bp sequence were selected to design a marker by running a BLAST search of the Rice Annotation Project Database (http://rapdb.dna.affrc.go.jp/).

Abbreviations

BS:

Brassinosteroid

CAPS:

Cleaved amplified polymorphic sequence

Chr6L:

Long arm of Chr6

Chr6S:

Short arm of Chr6

QTL:

Quantitative trait locus

SEM:

Scanning electron microscopy

SNP:

Single nucleotide polymorphism

SSR:

Simple sequence repeat

TRM:

TONNEAU1 recruiting motif

References

  • Ashikari M, Wu J, Yano M, Sasaki T, Yoshimura A (1999) Rice gibberellin-insensitive dwarf mutant gene Dwarf 1 encodes the alpha-subunit of GTP-binding protein. Proc Natl Acad Sci U S A 96:10284–10289

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  • Broman KW, Wu H, Sen Åš, Churchill GA (2003) R/qtl: QTL mapping in experimental crosses. Bioinformatics 19:889–890

    Article  CAS  PubMed  Google Scholar 

  • Fan C, Xing Y, Mao H, Lu T, Han B, Xu C, Li X, Zhang Q (2006) GS3, a major QTL for grain length and weight and minor QTL for grain width and thickness in rice, encodes a putative transmembrane protein. Theor Appl Genet 112:1164–1171

    Article  CAS  PubMed  Google Scholar 

  • Fan C, Yu S, Wang C, Xing Y (2009) A causal C-A mutation in the second exon of GS3 highly associated with rice grain length and validated as a functional marker. Theor Appl Genet 118:465–472

    Article  CAS  PubMed  Google Scholar 

  • Fujisawa Y, Kato T, Ohki S, Ishikawa A, Kitano H, Sasaki T, Asahi T, Iwasaki Y (1999) Suppression of the heterotrimeric G protein causes abnormal morphology, including dwarfism, in rice. Proc Natl Acad Sci U S A 96:7575–7580

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  • Hu X, Qian Q, Xu T, Zhang Y, Dong G, Gao T, Xie Q, Xue Y (2013) The U-Box E3 Ubiquitin Ligase TUD1 Functions with a Heterotrimeric G α Subunit to Regulate Brassinosteroid-Mediated Growth in Rice. PLoS Genet 9:e1003391

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  • Ishimaru K, Hirotsu N, Madoka Y, Murakami N, Hara N, Onodera H, Kashiwagi T, Ujiie K, Shimizu B, Onishi A et al (2013) Loss of function of the IAA-glucose hydrolase gene TGW6 enhances rice grain weight and increases yield. Nat Genet 45:707–711

    Article  CAS  PubMed  Google Scholar 

  • Izawa Y, Takayanagi Y, Inaba N, Abe Y, Minami M, Fujisawa Y (2010) Function and Expression Pattern of the α Subunit of the Heterotrimeric G Protein in Rice. Plant Cell Physiol 51:271–281

    Article  CAS  PubMed  Google Scholar 

  • Kato T, Segami S, Toriyama M, Kono I, Ando T, Yano M, Kitano H, Miura K, Iwasaki Y (2011) Detection of QTLs for grain length from large grain rice (Oryza sativa L.). Breed Sci 61:269–274

    Article  CAS  Google Scholar 

  • Kitagawa K, Kurinami S, Oki K, Abe Y, Ando T, Kono I, Yano M, Kitano H, Iwasaki Y (2010) A novel kinesin 13 protein regulating rice seed length. Plant Cell Physiol 51:1315–1329

    Article  CAS  PubMed  Google Scholar 

  • Konieczny A, Ausubel F (1993) A procedure for mapping Arabidopsis mutations using co-dominant ecotype-specific PCR-based markers. Plant J 4:403–410

    Article  CAS  PubMed  Google Scholar 

  • Li H, Durbin R (2010) Fast and accurate long-read alignment with Burrows–Wheeler transform. Bioinformatics 26:589–595

    Article  PubMed  PubMed Central  Google Scholar 

  • Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R (2009) The Sequence Alignment/Map format and SAMtools. Bioinformatics 25:2078–2079

    Article  PubMed  PubMed Central  Google Scholar 

  • Li Y, Fan C, Xing Y, Jiang Y, Luo L, Sun L, Shao D, Xu C, Li X, Xiao J et al (2011) Natural variation in GS5 plays an important role in regulating grain size and yield in rice. Nat Genet 43:1266–1269

    Article  CAS  PubMed  Google Scholar 

  • Li S, Liu W, Zhang X, Liu Y, Li N, Li Y (2012) Roles of the Arabidopsis G protein γ subunit AGG3 and its rice homologs GS3 and DEP1 in seed and organ size control. Plant Signal Behav 7:1357–1359

    Article  PubMed  PubMed Central  Google Scholar 

  • Mao H, Sun S, Yao J, Wang C, Yu S, Xu C, Li X, Zhang Q (2010) Linking differential domain functions of the GS3 protein to natural variation of grain size in rice. Proc Natl Acad Sci U S A 107:19579–19584

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  • McCouch SR, Teytelman L, Xu Y, Lobos KB, Clare K, Walton M, Fu B, Maghirang R, Li Z, Xing Y et al (2002) Development and mapping of 2240 new SSR markers for rice (Oryza sativa L.). DNA Res 9:199–207

    Article  CAS  PubMed  Google Scholar 

  • Qi P, Lin Y, Song X, Shen J, Huang W, Shan J, Zhu M, Jiang L, Gao J, Lin H (2012) The novel quantitative trait locus GL3.1 controls rice grain size and yield by regulating Cyclin-T1;3. Cell Res 22:1666–1680

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  • Segami S, Kono I, Ando T, Yano M, Kitano H, Miura K, Iwasaki Y (2012) Small and round seed 5 gene encodes alpha-tubulin regulating seed cell elongation in rice. Rice 5:4

    Article  PubMed  PubMed Central  Google Scholar 

  • Shomura A, Izawa T, Ebana K, Ebitani T, Kanegae H, Konishi S, Yano M (2008) Deletion in a gene associated with grain size increased yields during rice domestication. Nat Genet 40:1023–1028

    Article  CAS  PubMed  Google Scholar 

  • Song XJ, Huang W, Shi M, Zhu M, Lin H (2007) A QTL for rice grain width and weight encodes a previously unknown RING-type E3 ubiquitin ligase. Nat Genet 39:623–630

    Article  CAS  PubMed  Google Scholar 

  • Song XJ, Kuroha T, Ayano M, Furuta T, Nagai K, Komeda N, Segami S, Miura K, Ogawa D, Kamura T et al (2015) Rare allele of a previously unidentified histone H4 acetyltransferase enhances grain weight, yield, and plant biomass in rice. Proc Natl Acad Sci U S A 112:76–81

    Article  CAS  PubMed  Google Scholar 

  • Takano-Kai N, Jiang H, Kubo T, Sweeney M, Matsumoto T, Kanamori H, Padhukasahasram B, Bustamante C, Yoshimura A, Doi K et al (2009) Evolutionary history of GS3, a gene conferring grain length in rice. Genet 182:1323–1334

    Article  CAS  Google Scholar 

  • Tanabata T, Shibaya T, Hori K, Ebana K, Yano M (2012) SmartGrain: High-throughput phenotyping software for measuring seed shape through image analysis. Plant Physiol 160:1871–1880

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  • Wang S, Wu K, Yuan Q, Liu X, Liu Z, Lin X, Zeng R, Zhu H, Dong G, Qian Q et al (2012) Control of grain size, shape and quality by OsSPL16 in rice. Nat Genet 44:950–954

    Article  CAS  PubMed  Google Scholar 

  • Wang S, Li S, Liu Q, Wu K, Zhang J, Wang S, Wang Y, Chen X, Zhang Y, Gao C, Wang F, Huang H, Fu X (2015a) The OsSPL16-GW7 regulatory module determines grain shape and simultaneously improves rice yield and grain quality. Nat Genet 47:949–954

    Article  CAS  PubMed  Google Scholar 

  • Wang Y, Xiong G, Hu J, Jiang L, Yu H, Xu J, Fang Y, Zeng L, Xu E, Xu J, Ye W, Meng X, Liu R, Chen H, Wang Y, Zhu X, Li J, Qian Q (2015b) Copy number variation at the GL7 locus contributes to grain size diversity in rice. Nat Genet 47:944–948

    Article  CAS  PubMed  Google Scholar 

  • Ware DH, Jaiswal P, Ni J, Yap IV, Pan X, Clark KY, Teytelman L, Schmidt SC, Zhao W, Chang K et al (2002) Gramene, a tool for grass genomics. Plant Physiol 130:1606–1613

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  • Weng J, Gu S, Wan X, Gao H, Guo T, Su N, Lei C, Zhang X, Cheng Z, Guo X et al (2008) Isolation and initial characterization of GW5, a major QTL associated with rice grain width and weight. Cell Res 18:1199–1209

    Article  CAS  PubMed  Google Scholar 

  • Yan B, Liu R, Li Y, Wang Y, Gao G, Zhang Q, Liu X, Jiang G, He Y (2014) QTL analysis on rice grain appearance quality as exemplifying the typical events of transgenic or backcrossing breeding. Breed Sci 64:231–239

    Article  PubMed  PubMed Central  Google Scholar 

  • Zhang X, Wang J, Huang J, Lan H, Wang C, Yin C, Wu Y, Tang H, Qian Q, Li J et al (2012) Rare allele of OsPPKL1 associated with grain length causes extra-large grain and a significant yield increase in rice. Proc Natl Acad Sci U S A 109:21534–21539

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgments

We thank Drs. K. Takeda, I. Takamure, and T. Takita for informative suggestions about large grain rice lines.

This study was supported by grants from the Ministry of Agriculture, Forestry, and Fisheries of Japan (Genomics-based Technology for Agricultural Improvement, RBS-2001, IPG-0002) and JSPS KAKENHI Grant Number 25•10703.

Authors’ contribution

SS carried out the QTL analysis and drafted the manuscript. TY and KO conducted the QTL analysis. TN and MA performed the GoldenGate Assay. HK, HS, SM, and YK performed re-sequencing. HK developed all of the F2 populations. YI and KM conceived and participated in the design and coordination of the study and helped to draft the manuscript. All authors read and approved the final version of the manuscript.

Competing interests

The authors declare that they have no competing interests.

Author information

Authors and Affiliations

Authors

Corresponding authors

Correspondence to Yukimoto Iwasaki or Kotaro Miura.

Additional files

Additional file 1: Figure S1.

Frequency distribution of grain length and width in rice public dataset at NIAS Genebank. Frequency distribution of grain length (A) and grain width (B). Arrowheads indicate the trait data of BG23 and LG10 with measured place and year. (JPEG 612 kb)

Additional file 2: Figure S2.

Logarithm of odds (LOD) curves of five quantitative trait locus (QTL) analyses. LOD curves of QTL analyses for grain length are shown in A, C, E, G, and I. LOD curves of QTL analyses for grain length are shown in B, D, F, H, and J. F2 populations derived from Kasalath × BG23 (A and B), Kasalath × LG10 (C and D), Nipponbare × BG23 (E and F), Nipponbare × LG10 (G and H), and BG23 × LG10 (I and J) crosses. Dashed lines indicate LOD = 3.0 as threshold value. (JPG 541 kb)

Additional file 3: Figure S3.

Marker analysis of GW2, GS3, qSW5/GW5, and GW8 alleles in Nipponbare, Kasalath, BG23, and LG10 lines. The results of marker analysis of GW2 (A), GS3 (B), qSW5/GW5 (C), and GW8 (D) alleles. Arrows indicate the larger grain (L) and smaller grain (S) alleles. (JPG 1337 kb)

Additional file 4: Figure S4.

Heat map for a two-dimensional genome scan with two-quantitative trait locus (QTL) models in F2 population derived from the BG23 × LG10 cross. The heat map of the maximum logarithm of odds (LOD) score. Upper left triangle: interaction (Full–Add) model. Lower right triangle: full model. Color-coded scales indicate the values on the left for the interaction model (LOD threshold = 5.35) and on the right for the full model (LOD threshold = 8.13). (JPG 879 kb)

Additional file 5: Figure S5.

Selection of segregation lines to evaluate the effect of quantitative trait loci (QTLs) detected on Chr3 and Chr6L. F3 lines with one segregating locus for Chr3 or Chr6L and three other fixed loci were selected. From these two lines, F4 plants homozygous for BG23 or LG10 were selected to evaluate the effect of QTLs on Chr3 and Chr6L. White and black boxes indicate BG23 and LG10 chromosomes, respectively. (JPG 287 kb)

Additional file 6: Table S1.

Primer sequences used for the quantitative trait locus (QTL) analysis of four populations. (XLSX 51 kb)

Additional file 7: Table S2.

Primer sequences employed in the quantitative trait locus (QTL) analysis and genetic mapping using the population from the BG23 × LG10 cross. (XLSX 43 kb)

Additional file 8: Table S3.

Primer sequences used in marker analysis. (XLSX 29 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Segami, S., Yamamoto, T., Oki, K. et al. Detection of Novel QTLs Regulating Grain Size in Extra-Large Grain Rice (Oryza sativa L.) Lines. Rice 9, 34 (2016). https://doi.org/10.1186/s12284-016-0109-2

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12284-016-0109-2

Keywords