Next Article in Journal
Slow-Release Fertilizer Improves the Growth, Quality, and Nutrient Utilization of Wintering Chinese Chives (Allium tuberosum Rottler ex Spreng.)
Previous Article in Journal
Effect of Digestate on Soil Organic Carbon and Plant-Available Nutrient Content Compared to Cattle Slurry and Mineral Fertilization
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Agronomy
Volume 10
Issue 3
10.3390/agronomy10030380
agronomy-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Fine Mapping of a Grain Shape Gene from a Rice Landrace Longliheinuo-Dwarf (Oryza sativa L. ssp. japonica)

by
Fei Shang
,
Xu Chao
,
Kaiwen Meng
,
Xianghe Meng
,
Qin Li
,
Lifang Chen
and
Jianfei Wang
*
The State Key Laboratory of Crop Genetics and Germplasm Enhancement, Nanjing Agricultural University, Nanjing 210095, China
*
Author to whom correspondence should be addressed.
Agronomy 2020, 10(3), 380; https://doi.org/10.3390/agronomy10030380
Submission received: 7 February 2020 / Revised: 26 February 2020 / Accepted: 5 March 2020 / Published: 10 March 2020
Browse Figures
Review Reports Versions Notes

Abstract

:
Identification of grain shape genes can facilitate breeding of rice cultivars with optimal grain shape and appearance quality. In this study, we selected two rice germplasms, namely Longliheinuo-dwarf (LH) and N643, with different grain shape, to construct a genetic population for quantitative trait locus (QTL) analysis. A major QTL (qGS7), controlling the ratio of grain length to grain width, was mapped on the chromosome 7 in a BC1F4 line. By high-resolution linkage analysis, qGS7 was delimited to a 52.8 kb region including eight predicted genes. Through sequence alignment and real-time PCR expression analysis of these ORFs, ORF3 (LOC_Os07g42410) was selected as the candidate gene for further analysis. Single nucleotide polymorphisms (SNP) diversity analysis of ORF3 revealed that a single nucleotide deletion in the 7th exon resulted in a frameshift in parent LH and the parent in which a premature stop codon was identified. It was a rare mutation that caused grain shape difference. Real-time PCR analyses showed that the expression characteristics of ORF3 was in accordance with the development of spikelets. Of the 18 agronomic traits investigation in qGS7 near isogenic lines (NILs) showed that qGS7 not only changed grain shape but also affected plant height, panicle curvature, panicle length, the length of second leaf from the top, and chalkiness.

1. Introduction

As one of the most important cereal crops in the world, rice provides staple food for approximately half of the world population [1,2]. Grain shape is a key determinant of grain yield and quality in rice. In general, large grain size is associated with high yield but poor quality [3]. Grains with different shape are preferred by the consumers in different geographical locations worldwide [4]. Therefore, it is of great significance to identify grain shape genes and apply them to rice breeding.
To date, many genes controlling grain shape have been detected and cloned in rice [5]. Multiple signaling pathways are involved in regulating the grain shape, including the ubiquitin–proteasome pathway, protein kinase signaling pathway, phytohormones, transcriptional regulatory factors, and G-protein signaling pathways [6,7]. GS3 is known to encode a putative trans-membrane protein involved in the G protein signal pathway, which is the first major QTL regulating grain size [8,9]. GW2 and GW5 (also known as qSW5), negatively regulate rice grain width and weight. GW2 encodes a RING-type E3 ubiquitin ligase, while GW5 encodes a calmodulin-binding protein [10,11,12,13]. OsSPL13, encoded by GLW7, acts as an SBP (Squamosa promoter binding protein)-domain transcription factor and positively regulates grain length by affecting cell size [14]. GW8 encodes the transcription factor OsSPL16, which binds to the GW7 promoter and regulates grain shape by reducing GW7 expression [15,16]. GW6a, as a histone acetyltransferase, elevates its expression activity and raises grain weight and yield [17]. qGL3 encodes the protein phosphatase OsPPKL1, which participates in brassinosteroid-signaling in rice, and GL3.1 controls grain size by regulating cyclin-T1-3 [18,19,20]. GS5 and GS2/GL2 are also involved in the BR signaling pathway [21,22,23]. qTGW3 encodes the OsSK41 in the glycogen 38 synthase kinase 3/SHAGGY-like family. OsSK41 interacts with OsARF4 to negatively regulate grain size and weight [24]. TGW6 and BG1 also participate in the auxin pathway [25,26]. However, how these signaling pathways regulate the grain shape still awaits further research.
In this study, we used Longliheinuo-dwarf (LH) as the donor parent and N643 as the recurrent parent to develop different mapping populations for grain shape. A major QTL (qGS7), controlling the ratio of grain width to grain length, was mapped on the chromosome 7 in a BC1F4 line. Map-based cloning revealed that qGS7 was found to be allelic to EP2/DEP2/SRS1/cl7(t). The investigation of agronomic traits in the paired near isogenic lines (NILs) showed that qGS7 only affected the partially evaluated traits of plants and had no significant effect on the yield per plant. The results not only provide some new information for grain shape but also have an application value in appearance quality breeding in rice.

2. Materials and Methods

2.1. Plant Materials and Population Development

Longliheinuo-dwarf (LH) was obtained from Co60 irradiated Longliheinuo, and it was insensitive to gibberellin. N643 is a semi-dwarf line derived from the cross of two cytoplasmic male sterility maintenance lines, namely Longtepu and Yuefeng.
BC1F2 plants with various grain shape were selected to develop into BC1F2~4 lines. A BC1F4 line (RHL2401) was found to be separated only in grain shape and used to map the grain-shape gene. Short-grain plants of RHL2401 were back-crossed with N643 to develop a BC2F2 population for fine mapping. The heterozygotes from the BC2F2 population were selected to generate BC2F3 populations for further fine mapping (Figure 1).
In addition, heterozygous plants from the BC2 population were successively self-pollinated and kept residual heterozygous of qGS7 by marker assisted selection (MAS) to BC2F6. In the BC2F6 population, homozygotes were selected to form paired NILs for agronomic trait investigation (Figure 1).

2.2. Trait Evaluation

Rice materials were grown in the field at the Experimental Station of Nanjing Agricultural University (Jiangsu Province, China). Field management was carried out in keeping with the local standard methods [27]. Grain length and grain width were measured by electronic vernier caliper (0.01 mm) after the seeds were harvested and dried. Ten fully filled grains were randomly selected from each plant.

2.3. Molecular Marker Analysis

Simple sequence repeat (SSR) markers were designed as described on the Gramene database. Insert and deletion (InDel) markers were newly designed by comparison of DNA sequences between Nipponbare and 93-11 at the target region. The genomic DNA of each plant was extracted from fresh rice leaves according to a modified CTAB protocol. PCR was performed as follows: 95 °C for 5 min, followed by 35 cycles of 95 °C for 15 s, 55–60 °C for 20 s, 72 °C for 30 s, and a final elongation step at 72 °C for 10 min. The PCR products were analyzed on 6%–8% nondenaturing polyacrylamide gel.

2.4. Genetic Mapping

Bulked segregant analysis (BSA) was applied for the primary mapping. A total of 207 polymorphic SSR markers evenly distributed on 12 chromosomes were used to identify the genotype of individuals. The genetic map construction and QTL analysis were conducted by the composite interval mapping method using the IciMapping 3.1 software. The candidate genes were predicted using the data at http://rice.plantbiology.msu.edu/cgi-bin/gbrowse.

2.5. Sequencing Analysis and Expression Analysis

The sequence of gene were downloaded from Rice Genome Annotation Project database. Based on the sequence, primers were designed using Primer Premier 5 (Supplementary Materials Table S1). The genome sequence was amplified from the leaves using Phusion high fidelity DNA polymerase (Takara Biotech, Japan). Three independent purified PCR products were directly sequenced by Springen Biotechnology Co., Ltd. (Nanjing, China). Bio-XM software (Nanjing Agricultural University, China) was used for sequence alignment.
Total RNA was extracted from panicles at four different development stages with the RNA extraction kit (Bioteke, Wuxi, China), including pollen mother cell formation stage (young panicle 2–5 cm), pollen mother cell meiosis (young panicle 5–10 cm), pollen content filling stage (young panicle close to full length), and pollen maturation (anthers turn yellow, spikelets setting). The complementary DNA (cDNA) was synthesized with random oligo nucleotides utilizing a reverse transcription system (Vazyme Biotech, Nanjing, China). Quantitative real-time PCR (20 uL reaction volume) was conducted using 1 uL of cDNA, 0.5 uL of each primer, and 10 uL of AceQ qPCR SYBR Green Master Mix (Vazyme Biotech, Nanjing, China) in a Roche480 real-time PCR detection system. The fold change in cDNA relative to the reference gene (18s-rRNA) was determined by the comparative cycle threshold method. Relative expression levels of gene was calculated using the comparative Ct method [28].

2.6. Agronomic Traits Evaluation of Paired NILs for qGS7

Three paired qGS7 NILs were seeded in nursery trays containing nutrient matrix. After 30 days, the NILs were transplanted to experimental plots in a randomized complete block design with three replications. Each replication had forty plants in a density of 25 cm (4 lines) × 13.3 cm (10 plants). Twenty plants, excluding marginal effects, were collected to evaluated agronomy traits. The agronomy traits were surveyed included plant height, panicles per plant, stem diameter, top internode diameter, flag leaf length, flag leaf angle, length of second leaf from the top, panicle length, panicle exertion length, panicle curvature, primary branch number, secondary branch number, spikelet number per panicle, filled grain number per panicle, grain length, grain width, 1000-grain weight, and chalkiness.

3. Results

3.1. qGS7 is a Major QTL for Grain Length and Grain Width

N643 showed slender grains, and LH showed small round grains (Figure 2A). A BC1F4 line (RHL2401) was found to have separation in grain shape, and other agronomic traits were approximate to N643. The frequency distribution of grain length, grain width, and the ratio of grain length to grain width were both bimodal distributions (Figure 2B–D). The segregation of long-grain to short-grain plants in RHL2401 closely fitted the ratio of 3:1 (χ2 = 2.41, p > 0.05). The segregation ratio of narrow grain to wide grain plants also was 3:1 (χ2 = 3.57, p > 0.05). Furthermore, significant negative correlation (R = −0.776, p < 0.01) between the two traits, suggested that they were controlled by the same gene or by two closely linked genes.
The mixed DNA samples of RHL2401 plants were whole genome genotyped by 207 polymorphic markers (Supplementary Materials Figure S1). QTL mapping was performed in the heterozygous intervals with polymorphic markers. Two QTLs for grain length were identified on chromosomes 3 and 7, named qGL3 and qGL7, respectively (Table 1). The qGL3 was located between SSR markers RM7370 and RM6832, explaining 7.2% phenotypic variation with a LOD value of 8.85. This interval had a reported locus GS3, which was a major QTL for grain length and weight, and widely presented in indica rice. qGL7 was mapped between SSR markers RM6344 and RM8261, explaining 64.55% of phenotypic variation. It was a major QTL controlling grain length.
Two QTL for grain width were detected on chromosomes 5 and 7, named qGW5 and qGW7 (Table 1). The qGW5, between markers RM7444 and RM1089, explained 13.47% phenotypic variation with LOD value of 15.72, and its additive effect came from LH. In this interval existed a reported gene qSW5. qGW7 was located between markers RM6344 and RM8261, explaining 60.13% phenotypic variation with a LOD value of 35.15. Coincidentally, this interval was shared by the qGL7 and qGW7. This indicated that this fragment of LH was able to increase grain width and inhibit grain length. The corresponding N643 fragment could increase the grain length and reduce the grain width. A grain shape QTL for the ratio of grain length to width was also mapped in this interval, which explained 68.67% phenotypic variation with a LOD value of 43.94 (Table 1). It was determined that the single site, affecting grain length and grain width simultaneously, was named qGS7 (QTL for grain shape on chromosome 7).

3.2. Fine Mapping of qGS7

For fine mapping qGS7, new polymorphic InDel markers were developed (Supplementary Materials Table S1). In all, 605 BC1F2:5 plants were genotyped by those markers. Through phenotype assays of populations, qGS7 was located to an interval between markers In8 and In4 (Figure 3A). To isolate the qGS7 gene, 2089 plants from the BC2F2 population were used. The qGS7 were narrow down to 101 kb region between markers In13 and In4 (Figure 3B). High-resolution mapping was conducted using 6000 BC2F3 plants to further delimit the qGS7 locus to a 52.8 kb region that contained eight predicted genes (Figure 3C,D).

3.3. Candidate Genes Predication of qGS7

Among them, ORF2 and ORF8 encoded transposon protein. Therefore, we first cloned and sequenced the remaining six genes. ORF1 had no nucleotide variance between two parents. ORF3 had six SNPs in the coding region between two parents. SNP1 from T (LH) to A (N643) in the 7th exon caused amino acid changes from isoleucine to asparagine. SNP2 from one base deletion of LH resulted in a missense mutation, which changed the 502th and its subsequent amino acid sequence until terminated translation at position 508. SNP3 from T (LH) to G (N643) in the 7th exon caused amino acid changes from leucine to tryptophan. SNP4, SNP5, and SNP6 were synonymous mutations. ORF4 had four SNPs between two parents, but they did not cause amino acid changes. ORF5 had two SNPs between two parents. SNP1 was a synonymous mutation, and SNP2 from A (LH) to T (N643) caused premature termination of translation. ORF6 had one SNP between two parents, which caused amino acid changes from isoleucine to threonine. ORF7 had two SNPs between two parents, of which SNP1 was a synonymous mutation. SNP2 from A (LH) to G (N643) caused amino acid changes from valine to isoleucine (Table 2).
ORF1 and ORF4 did not cause amino acid difference between two parents. Therefore, the expression levels of ORF3, ORF5, ORF6, and ORF7 were examined at four different stages of panicle development in the paired near-isogenic lines (NIL-qGS7N643 and NIL-qGS7LH). ORF5 and ORF6 showed no expression in the two NILs plants. The expression level of the ORF3 and ORF7 both gradually decreased with panicle development in the NIL-qGS7N643 plants. In the NIL-qGS7LH plants, they increased at first, then decreased and were barely detectable when the panicles reached their final lengths (Figure 4). Previous studies demonstrated that ORF3 (EP2/DEP2/SRS1/cl7(t)) was the erect panicle and small round grain gene [29,30,31,32]. ORF7 encodes ribosomal protein S2. Consequently, ORF3 (LOC_Os07g42410) was selected as the candidate gene for further analysis.

3.4. SNP Diversity Analysis of qGS7

In order to analyze the SNP diversity of qGS7, nucleotide sequence analysis of LOC_Os07g42410 was further carried out in 16 rice varieties (including the parents). The sequencing fragments covered three missense mutations in the coding region (Table 3). The long and short grain were distinguished by the ratio of grain length to width. The results showed that variation in SNP1 and SNP3 were not related to grain shape, while SNP2 only occurred in the LH. Therefore, we inferred that SNP2 was the core mutation site that caused the grain shape difference in the two parents and may also be a rare mutation.

3.5. Agronomic Traits Evaluation of Paired NILs for qGS7

Of the 18 agronomic traits investigated, there were significant difference in seven agronomic traits between NIL-qGS7N643 and NIL-qGS7LH, including plant height, the length of second leaf from the top, panicle length, panicle curvature, grain length, grain width, and chalkiness. Compared to NIL-qGS7N643, NIL-qGS7LH showed that the plant height and its second leaves from top were shorter, the panicles length were shorter and erect, and the grains were shorter, wider, and had more chalkiness. These results indicated that qGS7 had pleiotropic effects. However, there was no significant difference in the yield traits consisting of panicls per plant, filled grain number per panicle, and the 1000-grain weight. These results indicated that qGS7 affected partial evaluated traits of plants and had no significant effect on yield per plant (Figure 5; Table 4).

4. Discussion

Grain size and shape are key agronomic traits and also affect grain yield and grain appearance in rice. Grain width, length, and thickness coordinately determine grain shape [33]. The physical appearance of grain dimensions (length, width, and the ratio of length to width) is essential for grain quality traits in rice. A number of QTLs are known to regulate grain shape in rice, of which several have been cloned and characterized [34]. However, our understanding remains fragmentary, and thus, further efforts to excavate elite genes/alleles are required.
In this study, we selected various grain shape plants of BC1F2 to develop a QTL mapping population. A major QTL (qGS7), controlling the ratio of grain length to grain width, was delimited to a region of 52.8 kb that contained eight predicted genes on chromosome 7 (Figure 3). Sequence alignment of predicted genes showed that the nucleotide variation of ORF3, ORF5, ORF6, and ORF7 caused the change of amino acids. The real time PCR analysis of four genes at different stages of panicle development showed that ORF5 and ORF6 had no expression in the two NILs plants. The expression level of the ORF3 was similar to ORF7 (Figure 4). Previous studies demonstrated that ORF3 (EP2/DEP2/SRS1/cl7(t)) was erect panicle and small round grain gene [29,30,31,32]. ORF7 encodes ribosomal protein S2. Therefore, we selected ORF3 as the candidate gene for further analysis.
The core mutation site between two parents in ORF3 was identified as SNP2, which a single nucleotide (T) deletion resulted in a frameshift in parent LH and premature termination of the protein. SNP diversity analysis of ORF3 revealed that SNP2 only occurred in LH, and other SNPs were not associated with grain shape (Table 3). We further compared the mutation sites of cl7 (t), EP2, DEP2, and SRS1 with qGS7. The results showed that although qGS7 and EP2/DEP2/SRS1/cl7(t) belonged to the same gene, they had different mutation site (Supplementary Materials Table S2). The grain size mutation in LH of the qGS7 candidate gene was novel compared to previous studies [29,30,31,32].
Expression levels of SRS1 were higher in the developing panicle than the fully developed one [29]. DEP2 gene expression became weaker as the panicles grew longer and was barely detectable when the panicles reached their final lengths [30]. The real time PCR analysis showed that expression level of qGS7 gradually decreased with panicle development (Figure 4). The result was consistent with previous studies.
Of the 18 agronomic traits investigated in qGS7 NILs (Figure 5; Table 4), there were significant difference of seven agronomic traits including plant height, second leaf length from the top, panicle length, panicle curvature, grain length, grain width, and chalkiness. However, there was no significant difference in the yield traits. The qGS7N643 increased grain length but decreased grain width, resulting in little effect on grain weight and enhanced effect on the ratio of grain length to width. In addition, milled rice from NIL-qGS7LH showed more chalkiness than that of NIL-qGS7N643. Although the first internode thickness from the top of NIL-qGS7LH plants was slightly thicker than NIL-qGS7 N643 in this study, they did not show significant differences. The first internode thickness from the top of EP2/DEP2 mutants is relatively thicker and contains a large number of vascular bundles, which can increase the strength of the internode [30,31,32].
To date, about 14 QTLs having large effects on grain length and width in rice have been cloned. One of them, GL7/GW7, has similar effects on grain length and width with opposite allelic directions, controlling grain shape but hardly influencing grain weight [16,35]. The other 13 genes affect grain size and weight. Four of them mainly control grain width, including GW2, GS5, qSW5/GW5, and GW8 [11,12,13,21]. Eight genes mainly control grain length, including GS2/GL2, OsLG3, qLGY3/OsLG3b, GS3, GL3.1/qGL3, GL4, TGW6, and GLW7 [34,36,37,38]. The remaining one, GW6a, has similar effects on grain length and width with the same allelic direction, and exhibits a larger impact on grain weight [17]. qGS7 and GL7/GW7, had similar effects on grain length and width, controlling grain shape without influencing grain weight. This type of QTL could be used for providing a potential gene resource for fine-tuning grain shape to modify grain appearance quality without yield penalty.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4395/10/3/380/s1, Figure S1. Graphical genotype of RHL2401 line on 12 rice chromosomes; Table S1. The primer used in this study; Table S2. Mutant site of erect panicle and small round gene.

Author Contributions

Conceptualization, X.M.; Formal analysis, Q.L. and L.C.; Investigation, X.C., K.M., and F.S.; Writing—original draft, F.S.; Funding acquisition, J.W. All authors have read and agreed to the published version of the manuscript.

Acknowledgments

This work was supported by the natural science foundation of Jiangsu Province, China (BK20151427).

Conflicts of Interest

The authors have declared that they have no conflict of interest.

References

  1. Sakamoto, T.; Matsuoka, M. Identifying and exploiting grain yield genes in rice. Curr. Opin. Plant Biol. 2008, 11, 209–214. [Google Scholar] [CrossRef] [PubMed]
  2. Wang, Z.K.; Zeng, D.D.; Qin, R.; Liu, J.L.; Shi, C.H.; Jin, X.L. A Novel and Pleiotropic Factor SLENDER GRAIN3 Is Involved in Regulating Grain Size in Rice. Rice Sci. 2018, 25, 132–141. [Google Scholar] [CrossRef]
  3. Xing, Y.Z.; Zhang, Q.F. Genetic and Molecular Bases of Rice Yield. Annu. Rev. Plant Biol. 2010, 61, 421–442. [Google Scholar] [CrossRef] [PubMed]
  4. Harberd, N.P. Shaping Taste: The Molecular Discovery of Rice Genes Improving Grain Size, Shape and Quality. J. Genet. Genom. 2015, 42, 597–599. [Google Scholar] [CrossRef]
  5. Huang, R.Y.; Jiang, L.R.; Zheng, J.S.; Wang, T.S.; Wang, H.C.; Huang, Y.M.; Hong, Z.L. Genetic bases of rice grain shape: So many genes, so little known. Trends Plant Sci. 2013, 18, 218–226. [Google Scholar] [CrossRef]
  6. Zuo, J.R.; Li, J.Y. Molecular Genetic Dissection of Quantitative Trait Loci Regulating Rice Grain Size. Annu. Rev. Genet. 2014, 48, 99–118. [Google Scholar] [CrossRef]
  7. Li, N.; Xu, R.; Li, Y. Molecular Networks of Seed Size Control in Plants. Annu. Rev. Plant Biol. 2019, 70, 435–463. [Google Scholar] [CrossRef]
  8. Fan, C.H.; Xing, Y.Z.; Mao, H.L.; Lu, T.T.; Han, B.; Xu, C.G.; Li, X.H.; Zhang, Q.F. 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. 2006, 112, 1164–1171. [Google Scholar] [CrossRef]
  9. Mao, H.L.; Sun, S.Y.; Yao, J.L.; Wang, C.R.; Yu, S.B.; Xu, C.G.; Li, X.H.; Zhang, Q.F. Linking differential domain functions of the GS3 protein to natural variation of grain size in rice. Proc. Natl. Acad. Sci. USA 2010, 107, 19579–19584. [Google Scholar] [CrossRef] [Green Version]
  10. Song, X.J.; Huang, W.; Shi, M.; Zhu, M.Z.; Lin, H.X. A QTL for rice grain width and weight encodes a previously unknown RING-type E3 ubiquitin ligase. Nat. Genet. 2007, 39, 623–630. [Google Scholar] [CrossRef]
  11. Weng, J.F.; Gu, S.H.; Wan, X.Y.; Gao, H.; Guo, T.; Su, N.; Lei, C.L.; Zhang, X.; Cheng, Z.J.; Guo, X.P.; et al. Isolation and initial characterization of GW5, a major QTL associated with rice grain width and weight. Cell Res. 2008, 18, 1199–1209. [Google Scholar] [CrossRef] [PubMed]
  12. Duan, P.G.; Xu, J.S.; Zeng, D.L.; Zhang, B.L.; Geng, M.F.; Zhang, G.Z.; Huang, K.; Huang, L.J.; Xu, R.; Ge, S.; et al. Natural Variation in the Promoter of GSE5 Contributes to Grain Size Diversity in Rice. Mol. Plant 2017, 10, 685–694. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Liu, J.F.; Chen, J.; Zheng, X.M.; Wu, F.Q.; Lin, Q.B.; Heng, Y.Q.; Tian, P.; Cheng, Z.J.; Yu, X.W.; Zhou, K.N.; et al. GW5 acts in the brassinosteroid signalling pathway to regulate grain width and weight in rice. Nat. Plants 2017, 3. [Google Scholar] [CrossRef] [PubMed]
  14. Si, L.Z.; Chen, J.Y.; Huang, X.H.; Gong, H.; Luo, J.H.; Hou, Q.Q.; Zhou, T.Y.; Lu, T.T.; Zhu, J.J.; Shangguan, Y.Y.; et al. OsSPL13 controls grain size in cultivated rice. Nat. Genet. 2016, 48, 447. [Google Scholar] [CrossRef]
  15. Wang, S.K.; Wu, K.; Yuan, Q.B.; Liu, X.Y.; Liu, Z.B.; Lin, X.Y.; Zeng, R.Z.; Zhu, H.T.; Dong, G.J.; Qian, Q.; et al. Control of grain size, shape and quality by OsSPL16 in rice. Nat. Genet. 2012, 44, 950. [Google Scholar] [CrossRef]
  16. Wang, S.K.; Li, S.; Liu, Q.; Wu, K.; Zhang, J.Q.; Wang, S.S.; Wang, Y.; Chen, X.B.; Zhang, Y.; Gao, C.X.; et al. The OsSPL16-GW7 regulatory module determines grain shape and simultaneously improves rice yield and grain quality. Nat. Genet. 2015, 47, 949. [Google Scholar] [CrossRef]
  17. Song, X.J.; Kuroha, T.; Ayano, M.; Furuta, T.; Nagai, K.; Komeda, N.; Segami, S.; Miura, K.; Ogawa, D.; Kamura, T.; et al. Rare allele of a previously unidentified histone H4 acetyltransferase enhances grain weight, yield, and plant biomass in rice. Proc. Natl. Acad. Sci. USA 2015, 112, 76–81. [Google Scholar] [CrossRef] [Green Version]
  18. Qi, P.; Lin, Y.S.; Song, X.J.; Shen, J.B.; Huang, W.; Shan, J.X.; Zhu, M.Z.; Jiang, L.W.; Gao, J.P.; Lin, H.X. The novel quantitative trait locus GL3.1 controls rice grain size and yield by regulating Cyclin-T1;3. Cell Res. 2012, 22, 1666–1680. [Google Scholar] [CrossRef] [Green Version]
  19. Zhang, X.J.; Wang, J.F.; Huang, J.; Lan, H.X.; Wang, C.L.; Yin, C.F.; Wu, Y.Y.; Tang, H.J.; Qian, Q.; Li, J.Y.; et al. Rare allele of OsPPKL1 associated with grain length causes extra-large grain and a significant yield increase in rice. Proc. Natl. Acad. Sci. USA 2012, 109, 21534–21539. [Google Scholar] [CrossRef] [Green Version]
  20. Gao, X.Y.; Zhang, J.Q.; Zhang, X.J.; Zhou, J.; Jiang, Z.S.; Huang, P.; Tang, Z.B.; Bao, Y.M.; Cheng, J.P.; Tang, H.J.; et al. Rice qGL3/OsPPKL1 Functions with the GSK3/SHAGGY-Like Kinase OsGSK3 to Modulate Brassinosteroid Signaling. Plant Cell 2019, 31, 1077–1093. [Google Scholar] [CrossRef]
  21. Li, Y.B.; Fan, C.C.; Xing, Y.Z.; Jiang, Y.H.; Luo, L.J.; Sun, L.; Shao, D.; Xu, C.J.; Li, X.H.; Xiao, J.H.; et al. Natural variation in GS5 plays an important role in regulating grain size and yield in rice. Nat. Genet. 2011, 43, 1266. [Google Scholar] [CrossRef] [PubMed]
  22. Che, R.H.; Tong, H.N.; Shi, B.H.; Liu, Y.Q.; Fang, S.R.; Liu, D.P.; Xiao, Y.H.; Hu, B.; Liu, L.C.; Wang, H.R.; et al. Control of grain size and rice yield by GL2-mediated brassinosteroid responses (vol 1, 15195, 2015). Nat. Plants 2016, 2. [Google Scholar] [CrossRef]
  23. Hu, J.; Wang, Y.X.; Fang, Y.X.; Zeng, L.J.; Xu, J.; Yu, H.P.; Shi, Z.Y.; Pan, J.J.; Zhang, D.; Kang, S.J.; et al. A Rare Allele of GS2 Enhances Grain Size and Grain Yield in Rice. Mol. Plant 2015, 8, 1455–1465. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Hu, Z.J.; Lu, S.J.; Wang, M.J.; He, H.H.; Sun, L.; Wang, H.R.; Liu, X.H.; Jiang, L.; Sun, J.L.; Xin, X.Y.; et al. A Novel QTL qTGW3 Encodes the GSK3/SHAGGY-Like Kinase OsGSK5/OsSK41 that Interacts with OsARF4 to Negatively Regulate Grain Size and Weight in Rice. Mol. Plant 2018, 11, 736–749. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Ishimaru, K.; Hirotsu, N.; Madoka, Y.; Murakami, N.; Hara, N.; Onodera, H.; Kashiwagi, T.; Ujiie, K.; Shimizu, B.; Onishi, A.; et al. Loss of function of the IAA-glucose hydrolase gene TGW6 enhances rice grain weight and increases yield. Nat. Genet. 2013, 45, 707. [Google Scholar] [CrossRef]
  26. Liu, L.C.; Tong, H.N.; Xiao, Y.H.; Che, R.H.; Xu, F.; Hu, B.; Liang, C.Z.; Chu, J.F.; Li, J.Y.; Chu, C.C. Activation of Big Grain1 significantly improves grain size by regulating auxin transport in rice. Proc. Natl. Acad. Sci. USA 2015, 112, 11102–11107. [Google Scholar] [CrossRef] [Green Version]
  27. Cheng, J.P.; Wang, L.; Du, W.L.; Lai, Y.Y.; Huang, X.; Wang, Z.F.; Zhang, H.S. Dynamic quantitative trait locus analysis of seed dormancy at three development stages in rice. Mol. Breed. 2014, 34, 501–510. [Google Scholar] [CrossRef]
  28. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2(T)(-Delta Delta C) method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef]
  29. Abe, Y.; Mieda, K.; Ando, T.; Kono, I.; Yano, M.; Kitano, H.; Iwasaki, Y. The SMALL AND ROUND SEED1 (SRS1/DEP2) gene is involved in the regulation of seed size in rice. Genes Genet. Syst. 2010, 85, 327–339. [Google Scholar] [CrossRef] [Green Version]
  30. Li, F.; Liu, W.B.; Tang, J.Y.; Chen, J.F.; Tong, H.N.; Hu, B.; Li, C.L.; Fang, J.; Chen, M.S.; Chu, C.C. Rice DENSE AND ERECT PANICLE 2 is essential for determining panicle outgrowth and elongation. Cell Res. 2010, 20, 838–849. [Google Scholar] [CrossRef]
  31. Zhu, K.M.; Tang, D.; Yan, C.J.; Chi, Z.C.; Yu, H.X.; Chen, J.M.; Liang, J.S.; Gu, M.H.; Cheng, Z.K. ERECT PANICLE2 Encodes a Novel Protein That Regulates Panicle Erectness in Indica Rice. Genetics 2010, 184, 343–354. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Ni, D.H.; Li, J.; Duan, Y.B.; Yang, Y.C.; Wei, P.C.; Xu, R.F.; Li, C.R.; Liang, D.D.; Li, H.; Song, F.S.; et al. Identification and utilization of cleistogamy gene cl7(t) in rice (Oryza sativa L.). J. Exp. Bot. 2014, 65, 2107–2117. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Huang, K.; Wang, D.K.; Duan, P.G.; Zhang, B.L.; Xu, R.; Li, N.; Li, Y.H. WIDE AND THICK GRAIN 1, which encodes an otubain-like protease with deubiquitination activity, influences grain size and shape in rice. Plant J. 2017, 91, 849–860. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Li, N.; Li, Y.H. Signaling pathways of seed size control in plants. Curr. Opin. Plant Biol. 2016, 33, 23–32. [Google Scholar] [CrossRef] [PubMed]
  35. Wang, Y.X.; Xiong, G.S.; Hu, J.; Jiang, L.; Yu, H.; Xu, J.; Fang, Y.X.; Zeng, L.J.; Xu, E.B.; Xu, J.; et al. Copy number variation at the GL7 locus contributes to grain size diversity in rice. Nat. Genet. 2015, 47, 944. [Google Scholar] [CrossRef] [PubMed]
  36. Wu, W.G.; Liu, X.Y.; Wang, M.H.; Meyer, R.S.; Luo, X.J.; Ndjiondjop, M.N.; Tan, L.B.; Zhang, J.W.; Wu, J.Z.; Cai, H.W.; et al. A single-nucleotide polymorphism causes smaller grain size and loss of seed shattering during African rice domestication. Nat. Plants 2017, 3. [Google Scholar] [CrossRef]
  37. Yu, J.P.; Xiong, H.Y.; Zhu, X.Y.; Zhang, H.L.; Li, H.H.; Miao, J.L.; Wang, W.S.; Tang, Z.S.; Zhang, Z.Y.; Yao, G.X.; et al. OsLG3 contributing to rice grain length and yield was mined by Ho-LAMap. BMC Biol. 2017, 15. [Google Scholar] [CrossRef] [Green Version]
  38. Yu, J.P.; Miao, J.L.; Zhang, Z.Y.; Xiong, H.Y.; Zhu, X.Y.; Sun, X.M.; Pan, Y.H.; Liang, Y.T.; Zhang, Q.; Rehman, R.M.A.; et al. Alternative splicing of OsLG3b controls grain length and yield in japonica rice. Plant Biotechnol. J. 2018, 16, 1667–1678. [Google Scholar] [CrossRef] [Green Version]
Agronomy 10 00380 g001 550
Figure 1. Schematic of the generation of the genetic materials used in this study.
Figure 1. Schematic of the generation of the genetic materials used in this study.
Agronomy 10 00380 g001
Agronomy 10 00380 g002 550
Figure 2. Grain shape characteristics of parents and frequency distribution of grain shape in the RHL2401 population. (A) Grain morphology of Longliheinuo-dwarf (LH), N643. Scale bar = 4 mm. Frequency distribution of grain length (B), grain width (C), and the ratio of grain length to grain width (D) in the RHL2401 population.
Figure 2. Grain shape characteristics of parents and frequency distribution of grain shape in the RHL2401 population. (A) Grain morphology of Longliheinuo-dwarf (LH), N643. Scale bar = 4 mm. Frequency distribution of grain length (B), grain width (C), and the ratio of grain length to grain width (D) in the RHL2401 population.
Agronomy 10 00380 g002
Agronomy 10 00380 g003 550
Figure 3. Map-based cloning of qGS7. (A) Physical map of the qGS7 locus using 605 BC1F2:5 plants. (B) Physical map of the qGS7 locus using 2089 BC2F2 plants. (C) High-resolution linkage analysis of the qGS7 locus using 6000 BC2F3 plants. Grey solid bars represent homozygous chromosomal segments for N643 and grille for heterozygote. (D) Predicted open reading frames in the target region of the rice genome.
Figure 3. Map-based cloning of qGS7. (A) Physical map of the qGS7 locus using 605 BC1F2:5 plants. (B) Physical map of the qGS7 locus using 2089 BC2F2 plants. (C) High-resolution linkage analysis of the qGS7 locus using 6000 BC2F3 plants. Grey solid bars represent homozygous chromosomal segments for N643 and grille for heterozygote. (D) Predicted open reading frames in the target region of the rice genome.
Agronomy 10 00380 g003
Agronomy 10 00380 g004 550
Figure 4. Real-time PCR analysis of candidate genes at four different stage with panicle development. (A) The expression level of ORF3. (B) The expression level of ORF7. Values are the means ± SD of three replicates.
Figure 4. Real-time PCR analysis of candidate genes at four different stage with panicle development. (A) The expression level of ORF3. (B) The expression level of ORF7. Values are the means ± SD of three replicates.
Agronomy 10 00380 g004
Agronomy 10 00380 g005 550
Figure 5. Agronomic traits of paired NILs for qGS7. (A) Gross morphology of qGS7N643 (left) and qGS7LH (right), bar = 25 cm. (B) Top internode diameter qGS7N643 (left) and qGS7LH (right), bar = 5 cm. (C) Flag leave length of qGS7N643 (left) and qGS7LH (right), bar = 5 cm. (D) Second leaf length from the top of qGS7N643 (up) and qGS7LH (down), bar = 5 cm. (E) Main panicle of qGS7N643 (left) and qGS7LH (right), bar = 5 cm. (F) Grain length of qGS7N643 (left) and qGS7LH (right), bar =1 cm. (G) Stem diameter of qGS7N643 (left) and qGS7LH (right), bar =1 cm. (H) Grain width of qGS7N643 (up) and qGS7LH (down), bar = 1 cm. (I) Comparison of the appearance of milled rice between the qGS7N643 (left) and qGS7LH (right), bar = 1 cm.
Figure 5. Agronomic traits of paired NILs for qGS7. (A) Gross morphology of qGS7N643 (left) and qGS7LH (right), bar = 25 cm. (B) Top internode diameter qGS7N643 (left) and qGS7LH (right), bar = 5 cm. (C) Flag leave length of qGS7N643 (left) and qGS7LH (right), bar = 5 cm. (D) Second leaf length from the top of qGS7N643 (up) and qGS7LH (down), bar = 5 cm. (E) Main panicle of qGS7N643 (left) and qGS7LH (right), bar = 5 cm. (F) Grain length of qGS7N643 (left) and qGS7LH (right), bar =1 cm. (G) Stem diameter of qGS7N643 (left) and qGS7LH (right), bar =1 cm. (H) Grain width of qGS7N643 (up) and qGS7LH (down), bar = 1 cm. (I) Comparison of the appearance of milled rice between the qGS7N643 (left) and qGS7LH (right), bar = 1 cm.
Agronomy 10 00380 g005
Table
Table 1. QTL detected in a BC1F4 line (RHL2401).
Table 1. QTL detected in a BC1F4 line (RHL2401).
TraitsChrQTLFlanking MarkersLODPVE (%)AddDom
Grain Length3qGL3RM7370-RM68328.857.200.323−0.189
7qGL7RM6344-RM826141.2664.56−0.9570.505
Grain Width5qGW5RM7444-RM108915.7213.470.1260.027
7qGW7RM6344-RM826135.1560.1360.268−0.133
Grain Length/Grain Width7qGS7RM6344-RM826143.9468.67−0.6060.268
Table
Table 2. Candidate genes for the qGS7 in the 52.8 kb region.
Table 2. Candidate genes for the qGS7 in the 52.8 kb region.
ORF Locus NameNucleotide LengthPutative FunctionCoding Sequence and Amino Acid Change
ORF1LOC_Os07g423952636DNA-directed RNA polymerase II subunit RPB9 No difference in coding region
ORF2LOC_Os07g424002229Transposon protein-
ORF3LOC_Os07g424107850Expressed protein7 SNPs, SNP1:T1145→A caused Ile382→Asn; SNP2: A base deletion at 1505 caused a frameshift and termination at 1525 position
ORF4LOC_Os07g4242056613-oxoacyl-synthase4 SNPs, no AA change
ORF5LOC_Os07g424301673Expressed protein2 SNPs, SNP1:T135→C, no AA change;
SNP2:T369→A caused a 2 AA deletion at C-terminal
ORF6LOC_Os07g424405540Glycolate oxidase1 SNP, C566→T caused Ile189→Thr
ORF7LOC_Os07g424502260Ribosomal protein S22 SNPs, SNP1:C488→T, no AA change
SNP2: A772→G caused Val258 →Ile
ORF8LOC_Os07g42460 2268Transposon protein-
Table
Table 3. SNP diversity analysis for qGS7.
Table 3. SNP diversity analysis for qGS7.
AccessionSubspGrain Length(mm)Grain Width (mm)Grain Length/Grain Width SNP1SNP2SNP3
114515051547
Longliheinuo-dwarfJaponica4.53.361.34T-G
N643Indica8.892.343.8ATT
80-4Japonica6.963.282.12TTG
Jiangdu6102Indica7.453.52.13ATT
J4195/zhen 511Japonica6.813.12.2TTG
NipponbareJaponica7.413.292.25TTG
SWWRIndica6.282.742.29ATT
Siyang 87-2566Japonica7.53.252.3TTG
Jiguang1Indica7.492.932.56ATT
93-11Indica9.842.563.84ATT
IR50Indica8.992.293.92ATT
YR303-304-1-3Indica9.832.254.37TTG
YR3030-12-6-39Indica9.642.144.5TTG
H94-1Indica9.672.074.67TTG
Table
Table 4. Agronomic traits of paired NILs for qGS7.
Table 4. Agronomic traits of paired NILs for qGS7.
TraitsLHN-643NIL-1NIL-2NIL-3
qGS7N643qGS7LHqGS7N643qGS7LHqGS7N643qGS7LH
Plant Height (cm)84.6 ± 1.5389.67 ± 3.06103.43 ± 3.9792.25 ± 3.30 ***118.95 ± 3.47105.43 ± 5.38 ***115.99 ± 5.12102.95 ± 4.94 ***
Panicles Per Plant3.00 ± 0.0013.33 ± 2.0810.67 ± 3.939.70 ± 1.779.50 ± 1.7810.50 ± 1.9610.10 ± 2.289.60 ± 1.96
Stem Diameter (mm)6.78 ± 0.175.80 ± 0.204.94 ± 0.555.36 ± 0.425.29 ± 0.795.05 ± 0.715.54 ± 0.635.53 ± 0.48
Top Internode Diameter (mm)4.56 ± 0.193.17 ± 0.173.24 ± 0.363.76 ± 0.34 *2.86 ± 0.423.19 ± 0.533.19 ± 0.493.46 ± 0.31
Flag Leaf Length (cm)22.44 ± 2.3332.00 ± 3.7728.43 ±3.4728.00 ± 9.2128.63 ± 2.0326.95 ± 1.9328.66 ± 2.2727.98 ± 3.02
Flag Leaf Angle (°)15.78 ± 0.5111.92 ± 0.1416.07 ± 4.7716.52 ± 3.4725.70 ± 3.5819.43 ± 2.8117.86 ± 2.1821.07 ± 2.57
Length of Second Leaf from the Top (cm)36.02 ± 2.7243.25 ± 3.8345.03 ± 3.1439.80 ± 1.75 **46.37 ± 2.2342.42 ± 2.37 ***46.39 ± 2.1943.48 ± 2.70 *
Panicle Length (cm)15.81 ± 0.1319.94 ± 0.9322.69 ± 2.1317.95 ± 3.14 **23.79 ± 1.2618.37 ± 0.91 **23.75 ± 1.3418.86 ± 0.77 ***
Panicle Exertion (cm)5.58 ± 0.92−3.21 ± 1.07−1.01 ± 1.70−1.56 ± 1.280.90 ± 0.652.11 ± 1.061.32 ± 0.713.10 ± 1.42
Panicle Curvature (°)4.11 ± 0.19104.36 ± 12.39120.96 ± 4.2131.98 ± 12.65 ***111.68 ± 9.8460.15 ± 9.99 ***113.83 ± 11.8961.93 ± 6.09 ***
Primary Branch Number13.11 ± 0.3813.65 ± 0.2613.93 ± 1.2514.72 ± 0.7415.18 ± 0.8714.75 ± 0.7615.34 ± 0.8815.26 ± 0.64
Secondary Branch Number29.00 ± 5.2930.07 ± 7.4231.02 ± 7.9628.79 ± 3.4534.49 ± 7.8032.47 ± 5.5835.19 ± 3.6036.87 ± 5.18
Spikelet Number Per Panicle169.33 ± 18.76171.21 ± 24.00192.21 ± 42.39185.54 ± 15.45209.90 ± 34.32191.75 ± 23.75211.27 ± 17.92211.12 ± 21.83
Filled Grain Number Per Panicle131.44 ± 12.22135.49 ± 16.59164.59 ± 36.68160.35 ± 15.54167.40 ± 29.87156.43 ± 19.56179.53 ± 18.86178.10 ± 20.71
Grain Length (mm)4.50 ± 0.078.89 ± 0.048.00 ± 0.256.80 ± 0.07 ***8.06 ± 0.137.06 ± 0.11 ***8.02 ± 0.146.97 ± 0.10 ***
Grain Width (mm)3.56 ± 0.082.34 ± 0.082.44 ± 0.052.86 ± 0.06 ***2.50 ± 0.042.87 ± 0.04 ***2.40 ± 0.032.85 ± 0.03 ***
1000-Grain Weight (g)11.77 ± 1.0519.37 ± 0.3218.98 ± 0.4718.75 ± 0.6719.81 ± 0.3019.34 ± 0.44 *18.88 ± 0.4218.97 ± 0.39
Chalkiness (%)glutinous10.4 ± 2.8811.25 ± 3.5924.28 ± 6.84 ***13.4 ± 2.4434.25 ± 6.23 ***15.65 ± 3.6532.30 ± 8.32 ***
*, ** and *** indicate significant difference at the 0.05, 0.01, and 0.001 probability level, respectively.

Share and Cite

MDPI and ACS Style

Shang, F.; Chao, X.; Meng, K.; Meng, X.; Li, Q.; Chen, L.; Wang, J. Fine Mapping of a Grain Shape Gene from a Rice Landrace Longliheinuo-Dwarf (Oryza sativa L. ssp. japonica). Agronomy 2020, 10, 380. https://doi.org/10.3390/agronomy10030380

AMA Style

Shang F, Chao X, Meng K, Meng X, Li Q, Chen L, Wang J. Fine Mapping of a Grain Shape Gene from a Rice Landrace Longliheinuo-Dwarf (Oryza sativa L. ssp. japonica). Agronomy. 2020; 10(3):380. https://doi.org/10.3390/agronomy10030380

Chicago/Turabian Style

Shang, Fei, Xu Chao, Kaiwen Meng, Xianghe Meng, Qin Li, Lifang Chen, and Jianfei Wang. 2020. "Fine Mapping of a Grain Shape Gene from a Rice Landrace Longliheinuo-Dwarf (Oryza sativa L. ssp. japonica)" Agronomy 10, no. 3: 380. https://doi.org/10.3390/agronomy10030380

APA Style

Shang, F., Chao, X., Meng, K., Meng, X., Li, Q., Chen, L., & Wang, J. (2020). Fine Mapping of a Grain Shape Gene from a Rice Landrace Longliheinuo-Dwarf (Oryza sativa L. ssp. japonica). Agronomy, 10(3), 380. https://doi.org/10.3390/agronomy10030380

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop