1. Introduction
As one of the main food crops, rice (Oryza sativa L.) feeds billions of people worldwide [1]. With increasing living standards, the demand of consumers has shifted the goal of rice breeding toward the improvement of rice grain quality [2]. Rice quality is a complex trait that includes milling quality, appearance quality, nutrition quality, eating and cooking quality (ECQ), and hygiene quality [3]. Among these, appearance quality (which is mainly determined by grain shape and chalkiness [4]) and ECQ (which is mainly determined by the amylose content) are the primary factors affecting the market value of rice [3].
Grain shape, which affects rice yield and quality simultaneously, is determined according to three traits: grain length, grain width, and grain thickness. Although a great number of major genes regulating grain shape have been cloned in recent years, there are still many minor genes to be discovered [5].
Chalkiness can be divided into white-belly, white-core, and white-back based on its distribution on the rice grain [6]. Chalkiness is an undesirable trait for consumers that negative affects rice grain appearance quality, milling quality, and ECQ [7]. The formation of chalkiness is regulated by both genetic factors and environmental factors [8]. Over past few decades, many chalkiness QTLs (quantitative trait loci) have been identified [9], but only a major chalkiness gene Chalk5 and a major white-core gene WCR1 have been cloned [10,11]. Chalk5 encodes a vacuolar H+-translocating pyrophosphatase (V-ppase) which disrupts the pH homeostasis of the developing endomembrane transport system, thereby affecting the biogenesis of protein bodies and forming air gaps between storage component in the endosperm, ultimately leading to chalkiness [10]. WCR1 encodes an F-box protein. OsDOF17 binds to the WCR1 promoter to increase its transcription level to promote the transcription of MT2b and inhibit 26S proteasome-mediated MT2b degradation, promoting the elimination of excess ROS, maintaining redox homeostasis, and delaying programmed cell death (PCD) in the endosperm, ultimately leading to a decrease in the white-core rate [11]. Thus, discovering more chalkiness genes is conducive to deepening our understanding of rice chalkiness.
Eating and cooking quality is reflected in multiple indices, including amylose content (AC), gel consistency (GC), gelatinization temperature (GT), and aroma and taste value [3]. AC is mainly regulated by a gene specifically expressed by the endosperm, Wx/GBSSI, which fluctuates between 0% to 30% according to the different allele of Wx [12]. However, there are also many minor genes that affect AC, such as Osbzip58, OsBP5, and Du1 [13,14,15].
In this study, two back-cross populations derived from two rice cultivars (Cypress and KY131) were developed and used for QTL mapping of the rice grain shape, chalkiness, and AC. As a result, 18 grain shape QTLs, 18 chalkiness QTLs, and 1 AC QTL were detected. The genetic effects of three grain shape QTLs—qGL7.2, qLWR7.2, and qGW3—and three chalkiness QTLs—qBR2, qWCR3, and qWCR11—were then confirmed in a random population. Using the progeny test method, qWCR3 was fine-mapped to a 100 kb region of chromosome 3.
2. Results
2.1. Phenotype Variation in BC3F2 and BC3F2:3 Populations
The grain length (GL), grain width (GW), length-to-width ratio (LWR), and amylose content (AC) had consistently higher values in Cypress than in KY131, while an opposite trend was observed for white-core rate (WCR), white-belly rate (WBR), white-back rate (BR), and chalkiness rate (CR) in both the BC3F2 (Hainan 2017) and BC3F2:3 populations (Wuhan 2018; Table 1). The GL, GW, and LWR in BC3F2 and BC3F2:3 populations showed a normal distribution (Figure 1A–F), as well as WBR, WCR, and AC in the BC3F2:3 population (Figure 1H,I,K); however, BR and CR showed skewed distributions in the BC3F2:3 population (Figure 1G,J).
Correlation analysis was conducted for all traits in the BC3F2:3 population. Among grain-shape-related traits, an extremely significant positive correlation was detected between GL and LWR, while an extremely significant negative correlation was found between GW and LWR. Among the chalkiness-related traits, BR was positively correlated to CR with high significance, and there were also significant positive correlations among WCR, WBR, and CR (Table 2). Moreover, GL was significantly positively correlated to WBR, while GW was negatively correlated to WBR and WCR with high significance.
2.2. Construction of Genetic Linkage Map
After screening nearly 1000 markers, 151 SSR (simple sequence repeat) and InDel (insertion and deletion) polymorphic makers between KY131 and Cypress were selected for genetic linkage map construction (Table S1). As a result, the genetic map covered a 948.6 cM genome region with an average interval size of 6.3 cM (Figure 2).
2.3. 37 QTLs Were Detected in BC3F2 and BC3F2:3 Populations
Using the genetic linkage analysis method, a total of 37 QTLs were detected on 11 chromosomes (excluding chromosome 10) in the BC3F2 and BC3F2:3 populations (Table 3).
A total of 18 grain-shape-related QTLs were found, of which there were seven for GL (qGL3, qGL4, qGL6, qGL7.1, qGL7.2, qGL7.3, and qGL8); the phenotypic variance explained by each QTL for GL ranged from 4.55% to 20.60%, and qGL7.1 was detected repeatedly both years. There were six QTLs for GW (qGW1, qGW3, qGW5.1, qGW5.2, qGW7.1, and qGW7.2), with qGW3 and qGW7.2 located in the same intervals as qGL3 and qGL7.2, respectively. We also detected five QTLs for LWR (qLWR3, qLWR5, qLWR7.1, qLWR7.2, and qLWR9), of which qLWR5 and qLWR7.1 were detected repeatedly both years, and qLWR7.1 was located in the same interval as qGW7.1.
A total of 18 QTLs for chalkiness were found in BC3F2:3 population, including four QTLs for BR (qBR1, qBR2, qBR9, and qBR11), four QTLs for WBR (qWBR1.1, qWBR1.2, qWBR5, and qWBR6), five QTLs for WCR (qWCR3, qWCR5.1, qWCR5.2, qWCR11, and qWCR12), and five QTLs for CR (qCR1.1, qCR1.2, qCR2, qCR9, and qCR11). Of these, qBR2 and qBR11 were located in common intervals with qCR2 and qWCR11, respectively.
For AC, only one QTL (qAC4) was detected in the BC3F2:3 population, explaining 4.01% of the phenotype variation.
2.4. Genetic Effect Validation of Seven QTLs
Based on the results of QTL mapping, a number of QTLs with large effects were selected for genetic effect validation. The seeds of BC3F2:3 plants with heterozygous genotype in each QTL interval were selected to construct near-isogenic lines for genetic effect validation. As shown in Table 4, three grain-shape-related QTLs (qGW3, qGL7.2, and qLWR7.2) and four chalkiness-related QTLs (qWBR1.2, qBR2, qWCR3, and qWCR11) were chosen for genetic effect validation. In particular, qWBR1.2, qGL7.2, and qGW3 were validated in one environment, while qLWR7.2, qBR2, qWCR3, and qWCR11 were validated in two environments.
2.5. qWCR3 Was Fine-Mapped to a 100 kb Region
We selected a stable QTL, qWCR3, for further research. For the fine mapping of qWCR3, a BC3F4 population consisting of 200 individuals was developed, and eight recombinants between D3-6 and C3.26.6 were identified. All recombinants were used for progeny testing of qWCR3, and 10 InDel makers were developed to genotype these recombinants according to the sequence variation of the two parents (Figure 3A). The progeny test involves comparing the phenotypes of two homozygous genotype plants in the inbred line of each recombinant and, if there is a significant difference, the candidate gene is considered to be located in the heterozygous fragment of the recombinant; otherwise, it is considered to be in the homozygous fragment. Progeny test results are shown in Figure 3B, from which it can be seen that significant WCR differences were observed between homozygous plants of lines 5–1 and 5–2, while no significant WCR differences were detected in other recombinant lines, suggesting that the candidate gene should be located in the 100 kb region between makers 25.49 and 25.59, and that qWCR3 co-segregated with maker 25.52.
2.6. LOC_Os03g45210 Could Be the Candidate Genes of qWCR3
The target 100 kb region of qWCR3 contains 16 annotated genes (Table 5), according to the annotation information of the japonica variety Nipponbare from the RGAP website (Rice Genome Annotation Project;
3. Discussion
Rice grain quality reflects the demands of consumers, sellers, and producers, consists of complex traits, and is one of the most important goals of rice breeding. Mapping-based cloning can help us identify QTLs (quantitative trait loci) that regulate rice grain quality traits, thus providing a theoretical basis for rice breeding. In this study, we identified 37 QTLs for grain shape, chalkiness, and amylose content in back-cross populations between Cypress and KY131.
3.1. Cloned Genes in the QTL Mapping Intervals
Among these QTLs, qGL3 and qGW3 were located in the same interval, which contains two cloned genes: PGL1 and BG1 [16,17]. PGL1 encodes a bHLH protein which positively regulates grain length through regulation of cell length in the lemma [16]. BG1 encodes a membrane localization protein, which plays a role in regulating the transport and distribution of auxin, thus positively regulating grain size [17]. The interval of qGL6 contains a cloned gene, BU1, which regulates BR signaling pathways through OsBRI1 and RGA1 and controls the degree of curvature of rice leaf nodes [18]. The interval of qGL8 contains a cloned gene, GAD1, which encodes a small secretory signal peptide leading to shorter grain length [19]. The interval of qGW5.1 contains a cloned gene, GS5, which encodes a serine carboxypeptidase and positively regulates grain width [20]. qGW5.2 and qLWR5 are located in the same interval which contains two genes, GSK2 and OSCYP51G3, that can regulate grain shape [21,22]. The intervals of qGW5.1, qGW5.2, and qLWR5 share partial overlaps, with the overlapping area containing a grain width-controlling gene, GW5. A deletion mutation in GW5 led to down-regulation of its expression level and to an increase in grain width. In addition, the GW5 protein can inhibit the activity of GSK2, thereby regulating the positive expression of genes in the downstream BR signaling pathway and positively regulating grain width [23]. The qLWR3 interval contains a cloned gene, GS3. Through competitively binding with the rice G protein β subunit RGB1, GS3 inhibits the signal pathway of the G protein, thus regulating grain length [24,25]. The overlapping regions of qGL7.2, qGW7.2, qLWR7.1, and qLWR7.2 contain five cloned genes, GLW7, GL7, OsBZR1, BG2, and SRS1/DEP2, which are related to grain shape [26,27,28,29,30]. These QTLs containing cloned genes may be used in rice breeding for the improvement of rice quality.
3.2. Four Chalkiness QTLs Were Newly Found
The formation of rice chalkiness is regulated by complex genetic networks and environments. Over the past few decades, researchers have detected a large number of quantitative trait loci controlling rice chalkiness [9]; however, few QTLs/genes have been fine-mapped or cloned. In this study, we detected 18 QTLs for chalkiness and verified four QTLs (qWBR1.2, qBR2, qWCR3, and qWCR11) repeatedly. The Cypress allele at qWCR3 effectively increased the white-core rate (Table 3 and Table 4), while the Cypress alleles at qWBR1.2, qBR2, and qWCR11 were found to effectively decrease the white-belly rate, white-back rate, and white-core rate, respectively. These results could serve as a basis for further fine mapping of the four QTLs and could be used in breeding programs to reduce chalkiness.
3.3. LOC_Os03g45210 Could Be a New Gene for Rice Chalkiness
qWCR3 was fine-mapped to a 100 kb region, where only LOC_Os03g45210 has an 11 bp InDel variation on the second exon between the two parents, leading to frameshift mutation and premature protein termination. LOC_Os03g45210 encodes a protein belonging to the plant cysteine oxidase (PCOs) family, which act as O2-sensing enzymes in plants. The PCOs family in Arabidopsis controls the stability of ERF-VIIs through the N/Arg degradation pathway, first through oxidization of the N-terminal cysteine of ERF-VIIs proteins, then by carrying out Argyl-tRNA-protein transferase 1 (ATE1)-mediated arginine installation to control the stability of ERF-VIIs, thus regulating oxygen homeostasis [31]. When plants are submerged, the oxygen concentration level decreases and the activity level of PCOs decreases accordingly, leading to enhanced stability of ERF VIIs and accumulation of ethylene, allowing plants to adapt to the transcription of genes related to hypoxic conditions [32]. At present, the cloned white-core rate gene—WCR1—is known to be related to oxygen homeostasis in rice endosperm cells [11]. Taken together, the candidate gene of qWCR3 may function as an O2 sensor in rice, controlling the formation of chalkiness through regulation of oxygen homeostasis in the endosperm.
4. Materials and Methods
4.1. Plant Materials and Field Experiment
A BC3F1 population consisting of 313 lines was derived from a cross between a temperate japonica cultivar KY131 (the recurrent parent) and a tropical japonica cultivar Cypress (the donor parent). The BC3F2 population consisted of 286 lines, with each line planted in one row including 12 plants. The BC3F2 population and their parents were planted in the experimental field of Huazhong Agricultural University in Lingshui city, Hainan province, in 2017. When the seeds of the BC3F2 plant matured, we mix-harvested the seeds of the BC3F2 lines to produce the BC3F2:3 lines, ensuring that the seeds of each of the BC3F2 lines corresponded strictly to a line of the BC3F2:3 population, i.e., each BC3F2:3 line originated from a mixed seed of a specific BC3F2 line. All the BC3F2:3 lines and their parents were planted in the experimental field of Huazhong Agricultural University in Wuhan city, Hubei province, in 2018. The phenotypes of rice quality-related traits in the BC3F2 and BC3F2:3 populations were collected for QTL mapping analysis.
In 2019, 2020, 2021, and 2022, the BC3F2 and BC3F2:3 lines with double heterozygous genotypes in the target QTL interval were planted for evaluation of genetic effects in the experimental field of Huazhong Agricultural University in Wuhan and Hainan. For the fine mapping of qWCR3, eight recombinants between D3-6 and C3.26.6 were identified from a BC3F4 population derived from a BC3F3 plant with a double heterozygous qWCR3 segment, where all recombinants were planted in 8 rows with 12 plants per row.
The 30-day-old progeny seedlings of these populations were transplanted into eight rows with 12 seedlings each, having 16.5 cm spacing in single-row plots in the field, while the rows were 26.4 cm apart. Field management followed local practices. Ten plants from the middle of each row were harvested individually for the measurement of traits.
4.2. Phenotyping and Statistical Analysis
Before measurement, harvested grains were air-dried at room temperature for at least 3 months to maintain a consistent water balance. The grain length (GL), grain width (GW), and grain length-to-width ratio (LWR) were measured after threshing using a BenQ scanner and analyzed with SmartGrain software [33]. Grain chalkiness, including white-belly (WBR), white-core (WCR), and white-back (BR) phenotypes, was assessed visually. The percentage of chalky grains in the total number of dehulled grains was used for measurement of the grain chalkiness rate. The chalkiness rate (CR) was calculated as the sum of the three types of chalkiness. Flour ground from milled grain was used to measure the amylose content using the iodine staining method [34].
The grain shape, chalkiness, and amylose content of two different homozygous genotypes in the same QTL were compared through a t-test. If the p-value was less than 0.05, the candidate gene was considered to be located in the heterozygous fragment, otherwise it was considered to be located in the homozygous fragment.
4.3. Genotyping and QTL Analysis
Leaves from the parent varieties KY131 and Cypress were used to extract genome DNA using the CTAB method [35], and the population’s leaves were used to extract genome DNA using the TPS method. The parent genomes were sequenced and assembled according to the rice reference genome (Rice Genome Annotation Project;
4.4. Progeny Testing Analysis
Progeny testing was conducted in the progeny segregation population of each recombinant. In the progeny line of each recombinant, a student’s t test was performed to compare the grain size differences between plants with a homozygous KY131 allele and a homozygous Cypress allele. If the p-value was less than 0.05, the candidate gene was considered to be located in the heterozygous fragment, otherwise it was considered to be located in the homozygous fragment.
4.5. Genetic Map Construction and QTL Analysis
Mapmaker/Exp3.0 was used to construct the genetic linkage map [37], while the composite interval mapping method included in the WinQTLCart software version 2.5 was used to perform QTL scanning on population phenotype data and genetic maps, as well as to screen QTLs with a logarithm of odds (LOD) threshold of 2.0 [38].
5. Conclusions
Using back-cross populations of Cypress and KY131, a genetic linkage map with 151 polymorphic markers was constructed and 37 QTLs for grain shape, chalkiness, and amylose content were detected through QTL mapping analysis. Among them, three QTLs for grain shape (qGW3, qGL7.2, and qLWR7.1) and four QTLs for chalkiness (qWBR1.2, qBR2, qWCR3, and qWCR11) presented repeatable and detectable genetic effects in different generation populations. Through progeny testing, we fine-mapped qWCR3 to a 100 kb region that contains 16 annotated genes, of which only LOC_Os03g45210 is specifically expressed in seed and endosperm, with an 11 bp InDel variation on the second exon between two parents leading to a frameshift mutation and premature protein termination. As a result, LOC_Os03g45210 could be the candidate gene for qWCR3. The results of this study lay the foundation for map-based cloning of these QTLs.
Conceptualization, Y.H.; methodology, J.L., G.L. and H.Z.; investigation, J.L., H.Z., Y.W., E.L., G.G., Q.Z. and H.S.; resources, G.J. and Y.H.; data curation, J.L. and H.Z.; writing, J.L.; supervision, Y.H.; funding acquisition, Y.H. All authors have read and agreed to the published version of the manuscript.
Not applicable.
Not applicable.
The datasets generated during the current study are available from the corresponding author upon reasonable request.
The authors extend their appreciation for the support received from the National Key Laboratory of Crop Genetic Improvement, Huazhong Agricultural University.
The authors declare that there are no conflicts of interest.
Footnotes
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Enlarge this image.Figure 1. Distribution of phenotypic values among the BC3F2 and BC3F2:3 populations derived from KY131/Cypress. Grain shape-related traits in 2017 (A–C), grain shape-related traits in 2018 (D–F), chalkiness-related traits (G–J), and amylose content (K) in 2018. GL, grain length; GW, grain width; LWR, length-to-width ratio; WBaR, white-back rate; WBeR, white-belly rate; WCR, white-core rate; CR, chalkiness rate; AC, amylose content.
Enlarge this image.Figure 2. Genetic linkage map based on 151 SSR markers. RM primers were derived from public databases, while the other InDel primers were designed based on parental sequence variations.
Enlarge this image.Figure 3. Fine mapping of qWCR3: (A) markers used in the fine mapping of qWCR3. (B) Progeny test of the eight recombinants narrowed qWCR3 to a 100 kb region. Black, white, and grey blocks represent the genotype of homozygous CYPRESS, homozygous KY131, and heterozygote, respectively.
Phenotypes of grain traits in the BC3F2 and BC3F2:3 populations.
| Trait | Parents | 2017HN BC3F2 | 2018WH BC3F2:3 | |||||||
|---|---|---|---|---|---|---|---|---|---|---|
| KY131 | Cypress | Mean | SD | MIN | MAX | Mean | SD | MIN | MAX | |
| GL (mm) | 6.81 | 8.67 | 6.99 | 0.29 | 5.37 | 7.87 | 6.94 | 0.01 | 6.14 | 7.67 |
| GW (mm) | 2.88 | 2.52 | 3.50 | 0.13 | 3.05 | 3.84 | 3.00 | 0.01 | 2.47 | 3.38 |
| LWR | 2.39 | 3.46 | 2.01 | 0.10 | 1.79 | 2.46 | 2.33 | 0.01 | 2.11 | 3.08 |
| BR (%) | 22.06 | 0.99 | 97.48 | 3.32 | 79.21 | 100 | ||||
| WCR (%) | 3.82 | 3.16 | 41.38 | 15.08 | 7.66 | 80.72 | ||||
| WBR (%) | 2.96 | 1.56 | 24.91 | 12.27 | 3.72 | 74.89 | ||||
| CR (%) | 24.42 | 4.98 | 98.41 | 2.24 | 82.67 | 100 | ||||
| AC (%) | 16.67 | 25.06 | 15.81 | 4.49 | 5.63 | 24.51 | ||||
Note: GL, grain length; GW, grain width; LWR, length-to-width ratio; BR, white-back rate; WBR, white-belly rate; WCR, white-core rate; CR, chalkiness rate; AC, amylose content. BC3F2 was planted in Hainan (HN) in 2017. BC3F2:3 was planted in Wuhan (WH) in 2018.
Correlation analysis of traits in the BC3F2:3 (KY131/Cypress) population.
| GL (mm) | GW (mm) | LWR | BR (%) | WCR (%) | WBR (%) | CR (%) | |
|---|---|---|---|---|---|---|---|
| GW (mm) | 0.20 ** | ||||||
| LWR | 0.50 ** | −0.74 ** | |||||
| BR (%) | 0.06 | 0.01 | 0.04 | ||||
| WCR (%) | −0.04 | 0.28 ** | −0.28 ** | 0.14 * | |||
| WBR (%) | −0.12 * | 0.24 ** | −0.29 ** | −0.01 | 0.55 ** | ||
| CR (%) | 0.03 | 0.08 | −0.04 | 0.93 ** | 0.22 ** | 0.13 * | |
| AC (%) | −0.04 | −0.01 | −0.02 | 0.03 | 0.24 ** | 0.01 | 0.04 |
Note: GL, grain length; GW, grain width; LWR, length-to-width ratio; BR, white-back rate; WBR, white-belly rate; WCR, white-core rate; CR, chalkiness rate; AC, amylose content. ** Correlation is significant at p < 0.01; * Correlation is significant at p < 0.05.
QTL mapping of KY131/Cypress populations.
| Traits | QTL | Chr | Interval | BC3F2 | BC3F2:3 | ||||
|---|---|---|---|---|---|---|---|---|---|
| LOD | ADD | V (%) | LOD | ADD | V (%) | ||||
| GL | qGL3 | 3 | T3-1-RM545 | 11.46 | 1.291 | 20.60 | |||
| qGL4 | 4 | Z4-20.6-Z4-26.5 | 3.06 | −0.188 | 4.85 | ||||
| qGL6 | 6 | Z6-6.0-TT6-1 | 4.81 | −0.140 | 5.88 | ||||
| qGL7.1 | 7 | RM481-Z7-5.4 | 3.69 | 0.125 | 4.55 | ||||
| qGL7.2 | 7 | T7-3-RM478 | 6.72 | −0.342 | 9.12 | 8.74 | −0.316 | 12.82 | |
| qGL7.3 | 7 | RM478-Z7-28.3 | 4.03 | −0.252 | 7.07 | ||||
| qGL8 | 8 | RM515-Z8-25.3 | 4.44 | −0.178 | 6.61 | ||||
| GW | qGW1 | 1 | RM486-A1.40.39 | 2.96 | 0.097 | 4.58 | |||
| qGW3 | 3 | T3-1-RM545 | 7.87 | 0.608 | 17.52 | ||||
| qGW5.1 | 5 | RM405-TT5-1 | 3.49 | 0.156 | 10.03 | ||||
| qGW5.2 | 5 | Z5-5.1-Z5-7.2 | 6.35 | 0.139 | 13.03 | ||||
| qGW7.1 | 7 | RM346-RM478 | 5.83 | 0.196 | 12.50 | ||||
| qGW7.2 | 7 | T7-3-RM478 | 5.36 | 0.140 | 8.94 | ||||
| LWR | qLWR3 | 3 | Z3-14.8-Z3-17.6 | 3.83 | −0.087 | 4.10 | |||
| qLWR5 | 5 | Z5-5.1-Z5-7.2 | 6.20 | −0.095 | 8.82 | 10.56 | −0.179 | 21.13 | |
| qLWR7.1 | 7 | RM346-RM478 | 25.80 | −0.237 | 33.26 | 17.73 | −0.246 | 26.86 | |
| qLWR7.2 | 7 | T7-3-RM346 | 10.19 | −0.203 | 13.20 | ||||
| qLWR9 | 9 | 9.8.0-Z9-10.4 | 3.72 | −0.275 | 14.78 | ||||
| BR | qBR1 | 1 | A1.40.39-RM104 | 2.69 | 6.398 | 4.25 | |||
| qBR2 | 2 | T2-1-T2-2 | 4.37 | 2.790 | 8.94 | ||||
| qBR9 | 9 | 9.2.0-9.8.0 | 3.03 | 2.680 | 3.79 | ||||
| qBR11 | 11 | C11-3.5-11-6.6 | 3.77 | −2.135 | 5.05 | ||||
| WBR | qWBR1.1 | 1 | C1-23.3-RM5 | 3.24 | 6.771 | 4.32 | |||
| qWBR1.2 | 1 | D1-8-B1.40.61 | 4.52 | −9.274 | 5.46 | ||||
| qWBR5 | 5 | Z5-5.1-Z5-7.2 | 4.54 | 8.925 | 6.73 | ||||
| qWBR6 | 6 | RM586-Z6-6.0 | 9.00 | −19.770 | 18.17 | ||||
| WCR | qWCR3 | 3 | D3-6-C3.26.6 | 2.54 | −10.607 | 4.88 | |||
| qWCR5.1 | 5 | Z5-5.1-TT5-1 | 2.98 | 10.873 | 5.75 | ||||
| qWCR5.2 | 5 | TT5-1-TT5-2 | 2.81 | 7.307 | 3.99 | ||||
| qWCR11 | 11 | C11-3.5-11.6.6 | 2.98 | −1.310 | 4.47 | ||||
| qWCR12 | 12 | Z12-19.9-Z12-25.2 | 4.38 | −10.997 | 7.02 | ||||
| CR | qCR1.1 | 1 | Z1-18.0-C1-23.3 | 3.25 | 3.393 | 26.30 | |||
| qCR1.2 | 1 | C1-23.3-RM5 | 4.04 | 1.448 | 5.59 | ||||
| qCR2 | 2 | T2-1-T2-2 | 3.43 | 1.675 | 7.22 | ||||
| qCR9 | 9 | H9.4.7-9.8.0 | 3.48 | 1.880 | 4.57 | ||||
| qCR11 | 11 | 11-6.6-Z11-16.9 | 2.88 | −1.026 | 3.76 | ||||
| AC | qAC4 | 4 | Z4-6.0-Z4-14.3 | 2.52 | 3.687 | 4.01 | |||
Note: QTL, quantitative trait locus; GL, grain length; GW, grain width; LWR, length-to-width ratio; BR, white-back rate; WBR, white-belly rate; WCR, white-core rate; CR, chalkiness rate; AC, amylose content; Chr, chromosome; LOD, logarithm of odds; ADD, additive effect, where positive values indicate that alleles from Cypress increase the trait scores; V, variance explained by the QTL.
Validation of QTLs in BC3F3 and BC3F4 populations.
| Year | QTL | Cypress Genotype | KY131 Genotype | p Value | ADD |
|---|---|---|---|---|---|
| 2019WH | qLWR7.2 | 2.89 ± 0.17 | 2.58 ± 0.20 | 0.0016 | 0.15 |
| qBR2 (%) | 40.05 ± 20.93 | 61.97 ± 18.97 | 0.033 | −10.96 | |
| qWBR1.2 (%) | 25.65 ± 11.52 | 8.66 ± 7.50 | 0.0019 | 8.50 | |
| 2020WH | qGL7.2 (mm) | 8.42 ± 0.39 | 7.70 ± 0.24 | 0.0059 | 0.36 |
| qGW3 (mm) | 3.15 ± 0.15 | 3.31 ± 0.08 | 0.0069 | −0.08 | |
| qLWR7.2 | 2.80 ± 0.07 | 2.50 ± 0.06 | 5.58 × 10−9 | 0.15 | |
| 2021WH | qBR2 (%) | 29.06 ± 6.19 | 41.81 ± 9.02 | 0.0026 | −6.37 |
| 2021HN | qWCR3(%) | 14.47 ± 5.90 | 2.46 ± 2.99 | 0.001 | 6.01 |
| qWCR11(%) | 3.91 ± 1.26 | 14.86 ± 3.86 | 6.14 × 10−7 | −5.47 | |
| 2022WH | qWCR3(%) | 8.71 ± 4.29 | 1.16 ± 0.75 | 1.09 × 10−7 | 3.78 |
| qWCR11(%) | 6.08 ± 4.24 | 12.02 ± 6.46 | 0.0070 | −2.97 |
Note: QTL, quantitative trait locus. WH, Wuhan. HN, Hainan. GL, grain length. GW, grain width. LWR, length-to-width ratio. BR, white-back rate. WBR, white-belly rate. WCR, white-core rate. ADD, additive effect; positive values indicate that alleles from Cypress increase the trait scores.
The ORFs in the fine-mapping region of qWCR3.
| Gene | Gene Product Name |
|---|---|
| LOC_Os03g45150 | LTP family protein precursor |
| LOC_Os03g45160 | Hypothetical protein |
| LOC_Os03g45170 | Amino acid permease |
| LOC_Os03g45180 | Expressed protein |
| LOC_Os03g45194 | Oxidoreductase |
| LOC_Os03g45210 | Plant Cysteine oxidase-3 |
| LOC_Os03g45220 | Expressed protein |
| LOC_Os03g45230 | Expressed protein |
| LOC_Os03g45250 | Plant Cysteine oxidase-2 |
| LOC_Os03g45260 | Vesicle transport v-SNARE protein |
| LOC_Os03g45270 | CS domain containing protein |
| LOC_Os03g45280 | Dehydrin |
| LOC_Os03g45290 | Ankyrin repeat domain-containing protein 50 |
| LOC_Os03g45300 | Transposon protein |
| LOC_Os03g45310 | Hypothetical protein |
| LOC_Os03g45320 | Dehydrogenase |
Supplementary Materials
The supporting information can be downloaded at:
References
1. Wing, R.; Purugganan, M.; Zhang, Q. The rice genome revolution: From an ancient grain to Green Super Rice. Nat. Rev. Genet.; 2018; 19, pp. 505-517. [DOI: https://dx.doi.org/10.1038/s41576-018-0024-z] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/29872215]
2. Tian, Z.; Qian, Q.; Liu, Q.; Yan, M.; Liu, X.; Yan, C.; Liu, G.; Gao, Z.; Tang, S.; Zeng, D. et al. Allelic diversities in rice starch biosynthesis lead to a diverse array of rice eating and cooking qualities. Proc. Natl. Acad. Sci. USA; 2009; 106, pp. 21760-21765. [DOI: https://dx.doi.org/10.1073/pnas.0912396106] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/20018713]
3. Zhou, H.; Xia, D.; He, Y. Rice grain quality-traditional traits for high quality rice and health-plus substances. Mol. Breed.; 2020; 40, 1. [DOI: https://dx.doi.org/10.1007/s11032-019-1080-6]
4. Fitzgerald, M.; McCouch, S.; Hall, R. Not just a grain of rice: The quest for quality. Trends Plant Sci.; 2009; 14, pp. 133-139. [DOI: https://dx.doi.org/10.1016/j.tplants.2008.12.004] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/19230745]
5. Zhao, D.; Zhang, C.; Li, Q.; Liu, Q. Genetic control of grain appearance quality in rice. Biotechnol. Adv.; 2022; 60, 108014. [DOI: https://dx.doi.org/10.1016/j.biotechadv.2022.108014] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/35777622]
6. Tan, Y.; Xing, Y.; Li, J.; Yu, S.; Xu, C.; Zhang, Q. Genetic bases of appearance quality of rice grains in Shanyou 63, an elite rice hybrid. Theor. Appl. Genet.; 2000; 101, pp. 823-829. [DOI: https://dx.doi.org/10.1007/s001220051549]
7. Chun, A.; Song, J.; Kim, K.; Lee, H. Quality of head and chalky rice and deterioration of eating quality by chalky rice. J. Crop Sci. Biotechnol.; 2009; 12, pp. 239-244. [DOI: https://dx.doi.org/10.1007/s12892-009-0142-4]
8. Tang, S.; Zhang, H.; Liu, W.; Dou, Z.; Zhou, Q.; Chen, W.; Wang, S.; Ding, Y. Nitrogen fertilizer at heading stage effectively compensates for the deterioration of rice quality by affecting the starch-related properties under elevated temperatures. Food Chem.; 2019; 277, pp. 455-462. [DOI: https://dx.doi.org/10.1016/j.foodchem.2018.10.137] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30502170]
9. Wu, B.; Xia, D.; Zhou, H.; Cheng, S.; Wang, Y.; Li, M.; Gao, G.; Zhang, Q.; Li, X.; He, Y. Fine mapping of qWCR7, a grain chalkiness QTL in rice. Mol. Breed.; 2021; 41, 68. [DOI: https://dx.doi.org/10.1007/s11032-021-01260-x] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/37309362]
10. Li, Y.; Fan, C.; Xing, Y.; Yun, P.; Luo, L.; Yan, B.; Peng, B.; Xie, W.; Wang, G.; Li, X. et al. Chalk5 encodes a vacuolar H(+)-translocating pyrophosphatase influencing grain chalkiness in rice. Nat. Genet.; 2014; 46, pp. 398-404. [DOI: https://dx.doi.org/10.1038/ng.2923]
11. Wu, B.; Yun, P.; Zhou, H.; Xia, D.; Gu, Y.; Li, P.; Yao, J.; Zhou, Z.; Chen, J.; Liu, R. et al. Natural variation in WHITE-CORE RATE 1 regulates redox homeostasis in rice endosperm to affect grain quality. Plant Cell; 2022; 34, pp. 1912-1932. [DOI: https://dx.doi.org/10.1093/plcell/koac057] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/35171272]
12. Mikami, I.; Uwatoko, N.; Ikeda, Y.; Yamaguchi, J.; Hirano, H.; Suzuki, Y.; Sano, Y. Allelic diversification at the wx locus in landraces of Asian rice. Theor. Appl. Genet.; 2008; 116, pp. 979-989. [DOI: https://dx.doi.org/10.1007/s00122-008-0729-z] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/18305920]
13. Bello, B.; Hou, Y.; Zhao, J.; Jiao, G.; Wu, Y.; Li, Z.; Wang, Y.; Tong, X.; Wang, W.; Yuan, W. et al. NF-YB1-YC12-bHLH144 complex directly activates Wx to regulate grain quality in rice (Oryza sativa L.). Plant Biotechnol. J.; 2019; 17, pp. 1222-1235. [DOI: https://dx.doi.org/10.1111/pbi.13048] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30552799]
14. Zhu, Y.; Cai, X.; Wang, Z.; Hong, M. An interaction between a MYC protein and an EREBP protein is involved in transcriptional regulation of the rice Wx gene. J. Biol. Chem.; 2003; 278, pp. 47803-47811. [DOI: https://dx.doi.org/10.1074/jbc.M302806200] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/12947109]
15. Zeng, D.; Yan, M.; Wang, Y.; Liu, X.; Qian, Q.; Li, J. Du1, encoding a novel Prp1 protein, regulates starch biosynthesis through affecting the splicing of Wxb pre-mRNAs in rice (Oryza sativa L.). Plant Mol. Biol.; 2007; 65, pp. 501-509. [DOI: https://dx.doi.org/10.1007/s11103-007-9186-3]
16. Heang, D.; Sassa, H. Antagonistic actions of HLH/bHLH proteins are involved in grain length and weight in rice. PLoS ONE; 2012; 7, e31325. [DOI: https://dx.doi.org/10.1371/journal.pone.0031325] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/22363621]
17. Liu, L.; Tong, H.; Xiao, Y.; Che, R.; Xu, F.; Hu, B.; Liang, C.; Chu, J.; Li, J.; Chu, C. Activation of Big Grain1 significantly improves grain size by regulating auxin transport in rice. Proc. Natl. Acad. Sci. USA; 2015; 112, pp. 11102-11107. [DOI: https://dx.doi.org/10.1073/pnas.1512748112] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/26283354]
18. Tanaka, A.; Nakagawa, H.; Tomita, C.; Shimatani, Z.; Ohtake, M.; Nomura, T.; Jiang, C.; Dubouzet, J.; Kikuchi, S.; Sekimoto, H. et al. BRASSINOSTEROID UPREGULATED1, encoding a helix-loop-helix protein, is a novel gene involved in brassinosteroid signaling and controls bending of the lamina joint in rice. Plant Physiol.; 2009; 151, pp. 669-680. [DOI: https://dx.doi.org/10.1104/pp.109.140806] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/19648232]
19. Jin, J.; Hua, L.; Zhu, Z.; Tan, L.; Zhao, X.; Zhang, W.; Liu, F.; Fu, Y.; Cai, H.; Sun, X. et al. GAD1 Encodes a Secreted Peptide That Regulates Grain Number, Grain Length, and Awn Development in Rice Domestication. Plant Cell; 2016; 28, pp. 2453-2463. [DOI: https://dx.doi.org/10.1105/tpc.16.00379] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/27634315]
20. Li, Y.; Fan, C.; Xing, Y.; Jiang, Y.; Luo, L.; Sun, L.; Shao, D.; Xu, C.; Li, X.; Xiao, J. et al. Natural variation in GS5 plays an important role in regulating grain size and yield in rice. Nat. Genet.; 2011; 43, pp. 1266-1269. [DOI: https://dx.doi.org/10.1038/ng.977]
21. Tong, H.; Liu, L.; Jin, Y.; Du, L.; Yin, Y.; Qian, Q.; Zhu, L.; Chu, C. DWARF AND LOW-TILLERING acts as a direct downstream target of a GSK3/SHAGGY-like kinase to mediate brassinosteroid responses in rice. Plant Cell; 2012; 24, pp. 2562-2577. [DOI: https://dx.doi.org/10.1105/tpc.112.097394] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/22685166]
22. Xia, K.; Ou, X.; Tang, H.; Wang, R.; Wu, P.; Jia, Y.; Wei, X.; Xu, X.; Kang, S.; Kim, S. et al. Rice microRNA osa-miR1848 targets the obtusifoliol 14α-demethylase gene OsCYP51G3 and mediates the biosynthesis of phytosterols and brassinosteroids during development and in response to stress. New Phytol.; 2015; 208, pp. 790-802. [DOI: https://dx.doi.org/10.1111/nph.13513] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/26083975]
23. Liu, J.; Chen, J.; Zheng, X.; Wu, F.; Lin, Q.; Heng, Y.; Tian, P.; Cheng, Z.; Yu, X.; Zhou, K. et al. GW5 acts in the brassinosteroid signalling pathway to regulate grain width and weight in rice. Nat. Plants; 2017; 3, 17043. [DOI: https://dx.doi.org/10.1038/nplants.2017.43] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/28394310]
24. Mao, H.; Sun, S.; Yao, J.; Wang, C.; Yu, S.; Xu, C.; Li, X.; Zhang, Q. Linking differential domain functions of the GS3 protein to natural variation of grain size in rice. Proc. Natl. Acad. Sci. USA; 2010; 107, pp. 19579-19584. [DOI: https://dx.doi.org/10.1073/pnas.1014419107] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/20974950]
25. Sun, S.; Wang, L.; Mao, H.; Shao, L.; Li, X.; Xiao, J.; Ouyang, Y.; Zhang, Q. A G-protein pathway determines grain size in rice. Nat. Commun.; 2018; 9, 851. [DOI: https://dx.doi.org/10.1038/s41467-018-03141-y] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/29487318]
26. Wang, Y.; Xiong, G.; Hu, J.; Jiang, L.; Yu, H.; Xu, J.; Fang, Y.; Zeng, L.; Xu, E.; Xu, J. et al. Copy number variation at the GL7 locus contributes to grain size diversity in rice. Nat. Genet.; 2015; 47, pp. 944-948. [DOI: https://dx.doi.org/10.1038/ng.3346] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/26147619]
27. Bai, M.; Zhang, L.; Gampala, S.; Zhu, S.; Song, W.; Chong, K.; Wang, Z. Functions of OsBZR1 and 14-3-3 proteins in brassinosteroid signaling in rice. Proc. Natl. Acad. Sci. USA; 2007; 104, pp. 13839-13844. [DOI: https://dx.doi.org/10.1073/pnas.0706386104] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/17699623]
28. Si, L.; Chen, J.; Huang, X.; Gong, H.; Luo, J.; Hou, Q.; Zhou, T.; Lu, T.; Zhu, J.; Shangguan, Y. et al. OsSPL13 controls grain size in cultivated rice. Nat. Genet.; 2016; 48, pp. 447-456. [DOI: https://dx.doi.org/10.1038/ng.3518] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/26950093]
29. Xu, F.; Fang, J.; Ou, S.; Gao, S.; Zhang, F.; Du, L.; Xiao, Y.; Wang, H.; Sun, X.; Chu, J. et al. Variations in CYP78A13 coding region influence grain size and yield in rice. Plant Cell Environ.; 2015; 38, pp. 800-811. [DOI: https://dx.doi.org/10.1111/pce.12452] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/25255828]
30. 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, pp. 327-339. [DOI: https://dx.doi.org/10.1266/ggs.85.327] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/21317545]
31. Taylor-Kearney, L.; Flashman, E. Targeting plant cysteine oxidase activity for improved submergence tolerance. Plant J.; 2022; 109, pp. 779-788. [DOI: https://dx.doi.org/10.1111/tpj.15605]
32. Licausi, F.; Kosmacz, M.; Weits, D.; Giuntoli, B.; Giorgi, F.; Voesenek, L.; Perata, P.; van Dongen, J. Oxygen sensing in plants is mediated by an N-end rule pathway for protein destabilization. Nature; 2011; 479, pp. 419-422. [DOI: https://dx.doi.org/10.1038/nature10536] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/22020282]
33. Tanabata, T.; Shibaya, T.; Hori, K.; Ebana, K.; Yano, M. SmartGrain: High-throughput phenotyping software for measuring seed shape through image analysis. Plant Physiol.; 2012; 160, pp. 1871-1880. [DOI: https://dx.doi.org/10.1104/pp.112.205120] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/23054566]
34. Bao, J.; Corke, H.; Sun, M. Nucleotide diversity in starch synthase IIa and validation of single nucleotide polymorphisms in relation to starch gelatinization temperature and other physicochemical properties in rice (Oryza sativa L.). Theor. Appl. Genet.; 2006; 113, pp. 1171-1183. [DOI: https://dx.doi.org/10.1007/s00122-006-0355-6] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/16850313]
35. Miura, K.; Agetsuma, M.; Kitano, H.; Yoshimura, A.; Matsuoka, M.; Jacobsen, S.E.; Ashikari, M. A metastable DWARF1 epigenetic mutant affecting plant stature in rice. Proc. Natl. Acad. Sci. USA; 2009; 106, pp. 11218-11223. [DOI: https://dx.doi.org/10.1073/pnas.0901942106] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/19541604]
36. Li, H.; Durbin, R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics; 2009; 25, pp. 1754-1760. [DOI: https://dx.doi.org/10.1093/bioinformatics/btp324] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/19451168]
37. Lander, E.; Green, P.; Abrahamson, J.; Barlow, A.; Daly, M.; Lincoln, S.; Newberg, L. MAPMAKER: An interactive computer package for constructing primary genetic linkage maps of experimental and natural populations. Genomics; 1987; 1, pp. 174-181. [DOI: https://dx.doi.org/10.1016/0888-7543(87)90010-3] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/3692487]
38. Wang, S.; Basten, C.; Zeng, Z. Windows QTL Cartographer 2.5. Department of Statistics, North Carolina State. 2006; Available online: http://statgen.ncsu.edu/qtlcart/WQTLCart.htm (accessed on 8 April 2024).
© 2024 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Details
; Gao, Guanjun 2 ; Zhang, Qinglu 2 ; Lou, Guangming 2 ; Jiang, Gonghao 3 ; He, Yuqing 2
1 National Key Laboratory of Crop Genetic Improvement and National Center of Plant Gene Research, Hubei Hongshan Laboratory, Huazhong Agricultural University, Wuhan 430070, China;
2 National Key Laboratory of Crop Genetic Improvement and National Center of Plant Gene Research, Hubei Hongshan Laboratory, Huazhong Agricultural University, Wuhan 430070, China;
3 College of Life Science, Heilongjiang University, Harbin 150080, China