Qing Xu1Jianhua Jiang1Chunyu Jing1,2Changmin Hu1,2Mengyuan Zhang1,2Xinru Li1,2Jiaming Shen1,2Mei Hai1,2Ying Zhang1Dezheng Wang1
Xiaojing Dang1*- 1Institute of Rice Research, Anhui Academy of Agricultural Sciences, Hefei, China
- 2College of Agronomy, Anhui Agricultural University, Hefei, China
Grain chalkiness directly affects the commercial value of rice. Genes related to chalkiness reported thus far have been discovered in mutants, but it has not been identified whether these genes can be used to improve rice quality by breeding. Therefore, discovering more quantitative trait loci (QTLs) or genes related to chalkiness in the rice germplasm is necessary. This study entails a genome-wide association study on the degree of endosperm chalkiness (DEC) and percentage of grains with chalkiness (PGWC) by combining 1.2 million single-nucleotide polymorphisms (SNPs) with the phenotypic data of 173 rice accessions. Thirteen QTLs for DEC and nine for PGWC were identified, of which four were detected simultaneously for both DEC and PGWC; further, qDEC11/qPGWC11 was identified as the major QTL. By combining linkage disequilibrium analysis and SNP information, LOC_Os11g10170 was identified as the candidate gene for DEC. There were significant differences among the haplotypes of LOC_Os11g10170, and the Hap 1 of LOC_Os11g10170 was observed to reduce the DEC by 6.19%. The qRT-PCR results showed that the gene expression levels in accessions with high DEC values were significantly higher than those in accessions with low DEC values during days 21–42 after flowering, with a maximum at 28 days. These results provide molecular markers and germplasm resources for genetic improvement of the chalkiness-related traits in rice.
Introduction
Chalkiness of rice refers to a white opaque part formed by loose deposition of starch and protein grains in the endosperm during the grain filling stage; it is negatively correlated with the appearance quality of rice and holds a certain relevance to the milling, eating, and cooking quality of rice (Chen et al., 2018). Chalkiness has always played an important role in the development of the rice industry. The degree of endosperm chalkiness (DEC) and percentage of grains with chalkiness (PGWC) are measures that directly affect the commercial value of rice; therefore, these are not only important indexes for evaluating rice chalkiness but also the main indexes for measuring the rice quality evaluation grade standard. High-chalkiness rice has obviously inferior mechanical strength than normal rice and easily breaks during processing, seriously affecting the quality, increasing the production cost, and reducing the market competitiveness of rice (Yoshioka et al., 2007). The taste of high-chalkiness rice after cooking can be significantly different from that of normal rice, in addition to its decreased viscosity and elasticity. Therefore, reducing rice chalkiness to improve rice quality is one of the important problems when breeding rice for high yield and good quality.
Chalkiness is a quantitative trait that is controlled by multiple genes with an obvious additive effect, and its inheritance is affected by the cytoplasmic genotype (Zheng et al., 2012). It also has maternal and endosperm effects (Shi et al., 2002). To the best of our knowledge, at least 46 quantitative trait loci (QTLs) controlling the DEC have been reported in literature (Tan et al., 2000; Zeng et al., 2002; Li et al., 2004; Wan et al., 2005; Liu et al., 2007; Chen et al., 2022), which were distributed on all 12 chromosomes of rice (Supplementary Table S1). Additionally, at least 84 QTLs controlling the PGWC have been reported (Li et al., 2003; Wan et al., 2005; Liu et al., 2007; Zhou et al., 2009; Wang et al., 2011; Zheng et al., 2012; Peng et al., 2014; Mirsa et al., 2021) (Supplementary Table S1), which were also distributed on all 12 chromosomes; of these 84 QTLs, 16 were found to be located on chromosome 6. Among these QTLs, only the Chalk5 QTL has been finely mapped further; its function has also been confirmed to encode a vacuolar H+-translocating pyrophosphatase that influences the grain chalkiness of rice (Li et al., 2014).
In recent years, several related studies on the flo-series chalkiness mutants have revealed numerous possible chalkiness production pathways and mechanisms. These genes are divided into three main categories. The first group is related to energy production and distribution of mutant genes; the Flo4 (Kang et al., 2005), Ogr1 (Kim et al., 2009), LOC_Os03g19890 (Jung et al., 2015), Flo10 (Wu et al., 2019), Flo12 (Zhong et al., 2019), Flo13 (Hu et al., 2018), Flo14 (Xue et al., 2019), Flo16 (Teng et al., 2019), Flo18 (Yu et al., 2021), and Flo19 (Lei et al., 2022) genes can directly or indirectly affect normal energy production in the grains. The Gif1 (Wang et al., 2008) and Nf-YB1 (Bai et al., 2016; Xu et al., 2016) genes can affect the source–sink relationship of nutrient transport directly or indirectly. The second category entails mutant genes related to starch synthesis; the Sbe (Yano et al., 1985), Wx (Mikami et al., 1999), SbeIIb (Tanaka et al., 2004), Flo5 (Ryoo et al., 2007), Pho1 (Satoh et al., 2008), Rsr1 (Fu and Xue, 2010), Agpl2 (Zhang et al., 2012), bZIP58 (Wang et al., 2013), Flo6 (Peng et al., 2014), Flo8 (Long et al., 2017), Flo20 (Yan et al., 2022), Flo22 (Yang et al., 2023), and Flo23 (Chen et al., 2023) genes can affect the normal starch synthesis in grains directly or indirectly. The third category includes mutant genes related to starch structure formation; the Rab5a (Wang et al., 2010), Pdil1-1 (Han et al., 2012), Gpa3 (Ren et al., 2014), Vps9a (Wen et al., 2015), and Bip (Yang et al., 2022) genes have been shown to induce a chalky phenotype by affecting the proteosome and vesicle systems, resulting in abnormal filling and accumulation of intracellular starch. There are a few other mutant genes related to chalkiness, such as Flo2 (Hikaru and Takeshi, 1981; She et al., 2010), Flo7 (Zhang et al., 2016), Flo11 (Zhu et al., 2018; Tabassum et al., 2020), and Flo15 (You et al., 2019), whose formation mechanisms have not been revealed fully. The genes related to chalkiness reported so far have been discovered from mutants, but it is not clear whether these genes can be used to improve rice quality by breeding. Therefore, discovering more QTLs or genes related to chalkiness in the rice germplasm is expected to provide not only the molecular basis but also the material basis for improving rice quality.
In this study, the phenotype values of the DEC and PGWC traits were investigated in 173 rice materials across two environments (E1 and E2). By combining the resequence data, a genome-wide association study (GWAS) was performed to identify the QTLs associated with DEC/PGWC and to predict the candidate genes. These results are expected to provide new gene resources and insights into the molecular mechanisms of DEC and PGWC as well as the means to genetically improve rice quality.
Materials and methods
Plant materials and field planting
A total of 173 rice accessions reported previously by Hu et al. (2022) were selected and subjected to phenotype measurements. These accession names, origins, subpopulations, and sequence read archive (SRA) accession numbers from the NCBI database are listed in Supplementary Table S2. These experimental varieties were planted at the Experimental Station of Anhui Academy of Agricultural Sciences (Hefei, Anhui Province, China) in 2021 and 2022. The field management followed routine operations. Each variety was planted twice repeatedly following a completely random design; there were 36 plants for each variety at each repetition.
Trait investigation
The grains were harvested 45 days after heading and air-dried under natural conditions for 3 months before carrying out the related experiments. Each experiment was performed twice, and the mean value of each variety from the two replicates was used to calculate its DEC and PGWC. For each repetition, about 50 whole milled rice samples were selected from each variety and placed on the glass plate of a MRS-9600 TFU2L Microtek scanner (Shanghai Zhongjing Technology Co., Ltd., Shanghai, China). These samples were then scanned using the Microtek ScanWizard EZ software. Using a rice appearance quality detection software, the proportion of the opaque parts was calculated to obtain the DEC and PGWC through the following formulas:
Genotypic data obtained
The sequence data of the 173 accessions were downloaded from the NCBI SRA with the accession number PRJNA554986 (Dang et al., 2020). The software Bowtie 2 (Li and Durbin, 2010) was then used to align all the paired-end sequences, for which the sequence of Os-Nipponbare-IRGSP 1.0 was used as the reference sequence. More than 95% of the reads were mapped to the reference genome with a mapping score exceeding 60. The HaplotypeCaller from GATK 3.8-0 was used to retrieve the single-nucleotide polymorphisms (SNPs). The software Beagle v4.1 (Browning and Browning, 2009) was used to complete the missing genotype data, and a total of 1,224,254 SNPs with minor allele frequencies of more than 5% and missing rates below 20% were selected to perform the GWAS.
Genome-wide association study
The compressed mixed linear model (MLM) was employed in the GWAS using the software GAPIT v.2.12 (Lipka et al., 2012). The “qqman” package of R software was used to draw the Manhattan plot (Turner, 2014). The significance threshold was set at 1.0 × 10-5, and the significant associated SNPs were determined on the basis of the correction method proposed by Benjamini and Hochberg (1995). The linkage disequilibrium (LD) analysis was carried out with the software Haploview 4.2, and the LD blocks surrounding the significant SNPs were detected (Barrett et al., 2005). The LD heatmap was then constructed using the “LDheatmap” program of R package (Shin et al., 2006). According to the method reported by Huang et al. (2021), there were more than three significant SNPs that exceeded the p-value threshold with clear peak-like signals within a 200-kb region of the leading SNP (having the smallest p-value), and this region was considered one QTL.
Candidate gene analysis
Using the Rice Genome Annotation Project MSU7 database (http://rice.plantbiology.Msu.edu), the candidate genes within the 200-kb genomic region were predicted. The causal genes were then determined based on gene annotation and expression analysis.
Quantitative real-time polymerase chain reaction (qRT-PCR) analysis of the candidate genes
After heading, the flowered spikelets of each day were marked. The spikelets marked for 7, 14, 21, 28, 35, and 42 days after flowering were sampled from three samples each having the highest and lowest DEC values. The TianGen Pure Plant Plus Kit (TianGen Biotech Co., Ltd., Beijing, China) was used to extract the total RNA. According to the instructions, the HisScript II Reverse Transcriptase system (Vazyme Biotech Co., Ltd., Nanjing, China) was used to synthesize the first-strand cDNA; the UBQ rRNA gene was used as the internal control. The qRT-PCR was then performed using SYBR Green (Vazyme Biotech Co., Ltd., Nanjing, China) based on the 96-well thermocycler (Roche Applied Science LightCycler 480, https://lifescience.roche.com/). The PCR procedures were performed in accordance with the description reported by Hu et al. (2022). For each sample, three biological repetitions were carried out, and the primers for the qRT-PCR are listed in Supplementary Table S3. The transcript levels of the gene expressions were calculated by 2−ΔΔCt, where ΔCt = Cttarget gene - CtUBQ rRNA (Livak and Schmittgen, 2001).
Haplotype analysis
The haplotypes of the candidate genes were determined in accordance with the RiceVarMap (http://ricevarmap.ncpgr.cn/) and China Rice Data Center (https://www.ricedata.cn/) databases. Each haplotype was confirmed by at least 20 accessions.
Results
Phenotypic analysis of chalkiness-related traits
For the DEC trait, the mean values were 12.33% ± 26.42% in E1 and 21.21% ± 26.85% in E2 (Figure 1A), with the mean coefficient of variation (CV) being 170.41% across the two environments (Figure 1A). For the PGWC trait, the maximum value was 100% while the minimum value was 0 between the two environments (Figure 1A). Figure 1B shows the distributions of the DEC and PGWC among the 173 accessions. Figure 1C displays the milled rice morphologies of the materials with different DEC values, where the accession Nuohangu shows the maximum DEC (100%) and accession Sihao4040 shows the minimum DEC (0.15%). Figure 1D shows the milled rice morphologies of the materials with different PGWC values, where the accession Jianongnuo2hao shows the maximum PGWC (100%) and accession Huajing6hao shows the minimum PGWC (0.59%). There were abundant variations in the chalkiness-related traits of the 173 rice accessions, which provided the material basis for mining the elite allele variations.
Figure 1. Phenotypic value descriptions of the degree of endosperm chalkiness (DEC) and percentage of grains with chalkiness (PGWC) traits among 173 rice accessions from two environments. (A) Basic statistics of the DEC and PGWC traits in the two environments. (B) Phenotypic value distributions of the DEC and PGWC traits in the two environments. The box edges indicate the upper and lower quantiles, and the lines in the middle of the boxes indicate the median values. The vertical lines indicate the data from the lowest quantile to the top quantile. (C) DEC and (D) PGWC performances in cultivars. Scale bar, 1 cm.
Identification of QTLs for chalkiness-related traits by GWAS
The GWAS was conducted with the 1,224,254 SNPs to investigate possible natural variations in the chalkiness-related traits, and Manhattan plots were constructed to show the significant SNP loci. For the DEC trait, 13 QTLs were identified, which were located on chromosomes 1, 3, 4, 6, 10, 11, and 12 (Figures 2A, B); among these, four QTLs were detected on chromosome 4. The QTL qDEC11 has the largest phenotypic variation explained (PVE) (24.89%). For the PGWC trait, a total of 9 QTLs were detected, which were located on chromosomes 3, 4, 6, 10, 11, and 12, with the PVE ranging from 11.89% to 18.18% (Figures 2C, D). The QTL qPGWC11 has the largest PVE (18.18%).
Figure 2. Quantitative trait loci (QTLs) identified in rice for the DEC and PGWC by the genome-wide association study (GWAS). (A) Manhattan plots for the DEC of the entire population of 173 rice accessions. (B) Information on the identified QTLs for DEC. (C) Manhattan plots for the PGWC of the entire population of 173 rice accessions. (D) Information on the identified QTLs for PGWC. Negative-log10-transformed p-values are plotted on the vertical axes, and the points above the red dashed lines show the significant single-nucleotide polymorphisms (SNPs) in the QTL regions. The QTLs identified are shown by black arrows. The QTLs shown in red indicate those that are overlapped for the DEC and PGWC traits.
By comparing the QTLs, five of them were identified to be common between the DEC and PGWC; these were qDEC4.1-qPGWC4.1, qDEC4.4-qPGWC4.4, qDEC6.2-qPGWC6, qDEC10-qPGWC10, and qDEC11-qPGWC11. Among these, both qDEC11 and qPGWC11 had the largest PVEs and largest numbers of significant SNP loci. Therefore, we consider qDEC11 to be the major QTL for DEC and further analyze it.
Identification of candidate genes for the DEC trait
To identify the candidate genes for the DEC trait, GWAS and LD analysis were conducted on chromosome 11 (Figure 3A); there were 23 candidate genes in the chromosome region of 5,388,827–5,588,827 containing qDEC11 (Figure 3B). Based on LD analysis, the LD block region was identified as 5,409,941–5,516,836, which contained 11 candidate genes (Supplementary Table S4). From the SNP information, five of the 11 genes contained non-synonymous SNPs, among which two genes located within the locus LOC_Os11g10170 were found by the GWAS to be significantly associated with DEC (Supplementary Tables S4, S5). The full length of the LOC_Os11g10170 locus is 1,716 bp, including four exons and three introns (Figure 3C). The LOC_Os11g10170 gene encodes a 387-amino acid protein. The SNPs occurred in the upstream, while the exons and introns were in the downstream sequence of the gene, resulting in the identification of three haplotypes (Figure 3C). The haplotype Hap 1 was associated with a smaller DEC, while the haplotypes Hap 2 and Hap 3 were associated with larger DEC values (Figure 3D). The SNP site at 5,515,531 changed from base T to base C at nt 649 in the cDNA sequence, resulting in a change from tryptophan (Trp) to arginine (Arg) at the amino acid 217. The SNP site at 5,516,142 changed from base G to base A at nt 1006 in the cDNA sequence, resulting in a change from aspartic acid (Asp) to asparagine (Asn) at the amino acid 336. The average DEC value of 106 accessions carrying Hap 1 was 6.19% ± 5.01%, while the corresponding average values for carrying Hap 2 and Hap 3 were 40.36% ± 37.82% and 42.14% ± 39.78%, respectively. There was no significant difference between Hap 2 and Hap 3 at p < 0.05. The average DEC values of Hap 2 and Hap 3 were significantly different from that of Hap 1 at p < 0.01 (Figure 3D). The average PGWC value of 101 accessions carrying Hap 1 was 21.82% ± 4.25%, while the corresponding average values for carrying Hap 2 and Hap 3 were 36.21% ± 12.58% and 36.60% ± 15.67%, respectively. There was no significant difference between Hap 2 and Hap 3 at p < 0.05. The average PGWC values of Hap 2 and Hap 3 were significantly different from that of Hap 1 at p < 0.01 (Figure 3E).
Figure 3. Identification of the candidate genes for the QTL qDEC11 in rice for DEC as determined from the GWAS and gene expression analysis data. (A) Identification of the linkage disequilibrium (LD) block of qDEC11. The Y axis shows the negative-log10-transformed p-values, and the points above the red line are the significant SNPs in the QTL region. The pairwise LD was calculated with the r2 value (determination coefficient between the SNP states). (B) Identification of the candidate genes in the QTL region of qDEC11. (C) Haplotypes of the LOC_Os11g10170 gene associated with DEC in rice. (D) Box plots of the DEC in accessions containing different haplotypes. (E) Box plots of the PGWC in accessions containing different haplotypes. (F) Expression analysis of the candidate gene LOC_Os11g10170 in the six materials with different AC values at different periods after flowering. The relative expressions were normalized to that of the rice UBQ gene. The error bars indicate the standard deviation, and the asterisks indicate the significant differences based on Student’s t-test (**p < 0.01).
The qRT-PCR results showed that there were no significant expression differences between the three accessions each with high and low DEC values from days 7 to 14 after flowering (Figure 3F); however, significant differences existed between the three accessions each with high and low DEC values from days 21 to 42 after flowering (Figure 3F). The expressions of LOC_Os11g10170 in the accessions with high DEC values showed gradual increases from days 21 to 28, reaching a maximum at 28 days. After 28 days, the expressions of LOC_Os11g10170 started declining gradually (Figure 3F). These results indicate that decreasing the expression of LOC_Os11g10170 could reduce the DEC.
Given that gene LOC_Os11g10140 contains the highest p-value SNP and LOC_Os11g10090 contains the premature stop codon SNP, their haplotype analyses were carried out. For gene LOC_Os11g10140, the transcript length was 1,017 bp, including four coding exons and three introns (Supplementary Figure S1A) as well as three haplotypes (Supplementary Figure S1A). The average DEC value of 83 accessions carrying Hap 1 was 14.70% ± 5.67%, while the corresponding average values for carrying Hap 2 and Hap 3 were 15.69% and 16.25%, respectively. There were no significant differences among the three haplotypes (Supplementary Figure S1B). The average PGWC value of 70 accessions carrying Hap 1 was 28.52% ± 9.83%, while the corresponding average values for carrying Hap 2 and Hap 3 were 26.59% and 28.90%, respectively. There were no significant differences among the three haplotypes (Supplementary Figure S1C).
For the gene LOC_Os11g10090, the transcript length was 2,046 bp, including four coding exons and six introns (Supplementary Figure S2A) as well as three identified haplotypes (Supplementary Figure S2A). The average DEC value of 76 accessions carrying Hap 1 was 14.97% ± 9.74%, while the corresponding average values for carrying Hap 2 and Hap 3 were 13.22% and 14.64%, respectively. There were no significant differences among the three haplotypes (Supplementary Figure S2B). The average PGWC value of 73 accessions carrying Hap 1 was 28.89% ± 7.48%, while the corresponding average values for carrying Hap 2 and Hap 3 were 26.44% and 28.43%, respectively. There were no significant differences among the three haplotypes (Supplementary Figure S2C).
Discussion
Rice quality is a complex trait, and endosperm opacity is one of the main attributes that determines the appearance quality of rice grains. Previous studies have shown that chalkiness is the result of interactions between the genetic background and environmental factors (Liao and Zhang, 2015; Zhou et al., 2016). The main environmental factors include the ecological environment, cultivation practices, and climate conditions (Xie et al., 2017). To reduce the influence of the environmental factors on phenotype, the chalkiness-related traits of 173 accessions were investigated in two environments, which showed great genetic variations. The CV ranged from 96.51% to 214.18%. These results provide a material basis for breeding to improve the chalkiness-related traits.
In this study, a total of 22 QTLs were detected, of which 13 controlled the DEC and nine controlled the PGWC traits. Among these QTLs, four common QTLs were identified for both DEC and PGWC, which were located on chromosomes 4, 6, 10, and 11 (Figure 2). By combining the information reported at the Gramene (http://www.gramene.org/markers/), BLAST (http://blast.ncbi.nlm.nih.gov/Blast.cgi), and China Rice Data Center database (http://www.ricedata.cn/gene/list/1499.htm) websites, the QTLs identified in this study were compared with those reported previously for controlling the DEC/PGWC. We compared the results identified in this study with those reported previously by Misra et al. (2019, 2021) and Sachdeva et al. (2024) and found that none of the identified loci were common, which could be related to the large differences in the geographical origins and characteristics of the varieties used. The varieties used by Misra et al. (2019, 2021) and Sachdeva et al. (2024) included indica, aus, tropical japonica from the IRRI South Asia Regional Center, Varanasi, Uttar Pradesh, India, whereas the materials used in this study included indica mainly from South China, East China, Central China, Vietnam, Philippines, and Indonesia as well as temperate japonica rice from East China, North China, and Northeast China. The results obtained in this study were also compared with those reported previously by Xin et al. (2022), and none of the loci identified were common; this may be attributed to the fact that the genetic backgrounds of the materials used were very different. The position ranges of six of the QTLs overlapped with the flanking regions of 12 of the QTLs reported previously, and the remaining 16 QTLs were newly identified in this study (Supplementary Table S6). Among the 16 QTLs newly identified herein, the QTL qDEC11 that is colocated with qPGWC11 had the largest PVE (24.89%) and highest number of significant SNP loci (225).
One GWAS signal that is significantly associated with DEC was detected to nearly single-gene resolution. In the LD block region 5,409,941–5,516,836 of chromosome 11, LOC_Os11g10170 is a newly identified candidate gene (Figure 3). The full length of LOC_Os11g10170 is 1,716 bp, including four exons and three introns; this gene encodes a 387 amino acid protein. The base T-to-C and G-to-A non-synonymous mutations in the cDNA sequence of LOC_Os11g10170 result in the high DEC phenotype (Figure 3). The qRT-PCR results showed that significant expression differences existed between the three accessions each with the highest and lowest DEC values from days 21 to 42 after flowering, reaching a maximum at 28 days (Figure 3F). However, we also suspected that the differences in these gene expression levels could be due to changes in the amino acids, leading to functional differences. Tryptophan produces indole-3-pyruvate in the presence of tryptophan transaminase, and indole-3-pyruvate produces auxin in the presence of YUCs. The change from tryptophan to arginine caused by the change in the base prevents the latter reaction, which affects the difference in gene expression level. The results of haplotype and qRT-PCR analyses indicate that LOC_Os11g10170 is the candidate gene for the DEC.
Auxin is an important plant hormone that plays a vital role in plant growth and development, while also participating in the formation and development of plant organs (Cheng et al., 2006, 2007). Indole-3-acetic acid (IAA) is the major auxin in plants, and the indole-3-pyruvate acid (IPA) pathway is the major tryptophan-dependent pathway for auxin biosynthesis in plants. For this pathway, the conversion of IPA to IAA by flavin monooxygenase encoded by the Yucca gene is the rate-limiting step (Hofmann et al., 2011; Stepanova et al., 2011). According to the website of the Rice Genome Annotation Project, the product encoded by LOC_Os11g10170 is flavin monooxygenase. Mashiguchi et al. (2011) demonstrated that Yucca-encoded flavin monooxygenase can catalyze the direct conversion of IPA to IAA based on gene function analysis, in vitro experiments, and isotope tracing observations. In the Arabidopsis thaliana genome, there are 11 Yucca genes that are mainly expressed in the apical meristems, young primordia, vascular tissues, and seeds (Cheng et al., 2006). In maize, the mutation of the Yucca1 gene was shown to result in sparse inflorescence and fewer branches (Gallavotti et al., 2008). In the rice genome, there are 14 Yucca homologous genes (Stepanova et al., 2011). OsYUCCA1 is located on chromosome 1 and affects auxin biosynthesis (Yamamoto et al., 2007). The OsFMO(t) gene located on chromosome 3 regulates IAA biosynthesis possibly locally and plays an important role in the formation of local IAA concentrations that are critical for regulating normal growth and development in rice (Yi et al., 2013). The OsFMO1 gene located on chromosome 9 can affect auxin synthesis, signal transduction, polar transport, and growth and development in plants (Liu et al., 2019). Grain filling is an important physiological process in rice growth; therefore, we hypothesize that LOC_Os11g10170 encoding flavin monooxygenase may inhibit or accelerate the transport of photosynthetic products through the phloem to the endosperm cells in the form of sucrose during grain filling to affect starch synthesis, resulting in different degrees of chalkiness among the accessions. This provides a possible scientific basis for further research on the LOC_Os11g10170 gene with regard to the developmental growth and quality of rice.
Data availability statement
The data presented in the study are deposited in the Sequence Read Archive (SRA), NCBI accession number PRJNA554986. The names of the repositories and accession numbers can be found in the article/Supplementary Material.
Author contributions
QX: formal analysis, investigation, software, and writing–original draft. JJ: data curation, formal analysis, and writing–original draft. CJ: investigation and writing–original draft. CH: investigation and writing–original draft. MZ: investigation and writing–original draft. XL: investigation and writing–original draft. JS: investigation and writing–original draft. MH: investigation and writing–original draft. YZ: investigation and writing–original draft. DW: investigation, resources, and writing–original draft. XD: funding acquisition, project administration, software, and writing–review and editing.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This research was supported by the Open Research Fund Program of Anhui Province Key Laboratory of Rice Germplasm Innovation and Molecular Improvement (grant nos. SDKF-2022-01 and SDKF-2023-01), Talent Program of Anhui Academy of Agricultural Sciences (grant nos. QNYC-201903 and QNYC-202103), and Natural Science Foundation of Anhui Province (grant no. 2108085MC100).
Acknowledgments
The authors would like to thank the Platform of National Crop Germplasm Resources of China for providing the Oryza sativa germplasms.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations or those of the publisher, editors, and reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2024.1423648/full#supplementary-material
References
Bai, A. N., Lu, X. D., Li, D. Q., Liu, J. X., and Liu, C. M. (2016). NF-YB1-regulated expression of sucrose transporters in aleurone facilitates sugar loading to rice endosperm. Cell Res. 26, 384–388. doi:10.1038/cr.2015.116
Barrett, J. C., Fry, B., Maller, J., and Daly, M. J. (2005). Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics 21, 263–265. doi:10.1093/bioinformatics/bth457
Benjamini, Y., and Hochberg, Y. (1995). Controlling the false discovery rate: a practical and powerful approach to multiple testing. J. R. Stat. Soc. 57, 289–300. doi:10.1111/j.2517-6161.1995.tb02031.x
Browning, B. L., and Browning, S. R. (2009). A unified approach to genotype imputation and haplotype-phase inference for large data sets of trios and unrelated individuals. Am. J. Hum. Genet. 84, 210–223. doi:10.1016/j.ajhg.2009.01.005
Chen, F., Wang, H. F., Gao, J., Zhou, X. B., Xin, C. Y., Zhang, S. Y., et al. (2018). Research progress on rice grain chalkiness. Barl. Cer. Sci. 35, 7–10. doi:10.14069/j.cnki.32-1769/s.2018.02.002
Chen, L., Ma, J., Liu, W., Sun, J. C., and Han, L. Z. (2022). QTL mapping for qCG5 of chalky trait in rice. J. Plant Genet. Res. 23, 1458–1464. doi:10.13430/j.cnki.jpgr.20220330003
Chen, X. L., Ji, Y., Zhao, W. Y., Niu, H. Y., Yang, X., Jiang, X. K., et al. (2023). Fructose-6-phosphate-2-kinase/fructose-2,6-bisphosphatase regulates energy metabolism and synthesis of storage products in developing rice endosperm. Plant Sci. 326, 111503. doi:10.1016/j.plantsci.2022.111503
Cheng, Y. F., Dai, X. H., and Zhao, Y. D. (2006). Auxin biosynthesis by the YUCCA flavin monooxygenases controls the formation of floral organs and vascular tissues in Arabidopsis. Gene Dev. 20, 1790–1799. doi:10.1101/gad.1415106
Cheng, Y. F., Dai, X. H., and Zhao, Y. D. (2007). Auxin synthesized by the YUCCA flavin monooxygenases is essential for embryogenesis and leaf formation in Arabidopsis. Plant Cell 19, 2430–2439. doi:10.1105/tpc.107.053009
Dang, X. J., Yang, Y., Zhang, Y. Q., Chen, X. G., Fan, Z. L., Liu, Q. M., et al. (2020). OsSYL2AA, an allele identified by gene-based association, increases style length in rice (Oryza sativa L). Plant J. 104, 1491–1503. doi:10.1111/tpj.15013
Fu, F. F., and Xue, H. W. (2010). Coexpression analysis identifies rice starch regulator1, a rice AP2/EREBP family transcription factor, as a novel rice starch biosynthesis regulator. Plant Physiol. 154, 927–938. doi:10.1104/pp.110.159517
Gallavotti, A., Barazesh, S., Malcomber, S., Hall, D., Jackson, D., Schmidt, R. J., et al. (2008). Sparse inflorescence1 encodes a monocot-specific YUCCA-like gene required for vegetative and reproductive development in maize. Pro. Natl. Acad. Sci. 105, 15196–15201. doi:10.1073/pnas.0805596105
Han, X. H., Wang, Y. H., Liu, X., Jiang, L., Ren, Y. L., Liu, F., et al. (2012). The failure to express a protein disulphide isomerase-like protein results in a floury endosperm and an endoplasmic reticulum stress response in rice. J. Exper. Bot. 63, 121–130. doi:10.1093/jxb/err262
Hikaru, S., and Takeshi, O. (1981). New endosperm mutations induced by chemical mutagens in rice, Oryza sativa L. Jap. J. Breed. 31, 316–326. doi:10.1270/jsbbs1951.31.316
Hofmann, N. R. (2011). YUC and TAA1/TAR proteins function in the same pathway for auxin biosynthesis. Plant Cell 23, 3869. doi:10.1105/tpc.111.231112
Hu, C. M., Jiang, J. H., Li, Y. L., Song, S. J., Zou, Y., Jing, C. Y., et al. (2022). QTL mapping and identification of candidate genes using a genome-wide association study for heat tolerance at anthesis in rice (Oryza sativa L.). Front. Genet. 13, 983525. doi:10.3389/fgene.2022.983525
Hu, T. T., Tian, Y. L., Zhu, J. P., Wang, Y. L., Jing, R. N., Lei, J., et al. (2018). OsNDUFA9 encoding a mitochondrial complex I subunit is essential for embryo development and starch synthesis in rice. Plant Cell Rep. 37, 1667–1679. doi:10.1007/s00299-018-2338-x
Huang, Z. B., Ying, J. F., Peng, L. L., Sun, S., Huang, C. W., Li, C., et al. (2021). A genome-wide association study reveals that the cytochrome b5 involved in seed reserve mobilization during seed germination in rice. Theor. Appl. Genet. 134, 4067–4076. doi:10.1007/s00122-021-03948-2
Jung, K. H., Kim, S. R., Giong, H. K., Nguyen, M. X., Koh, H. J., and An, G. (2015). Genome-wide identification and functional analysis of genes expressed ubiquitously in rice. Mol. Plant 8, 276–289. doi:10.1016/j.molp.2014.10.013
Kang, H. G., Park, A., Matsuoka, M., and An, G. (2005). White-core endosperm floury endosperm-4 in rice is generated by knockout mutations in the C-type pyruvate orthophosphate dikinase gene (OsPPDKB). Plant J. 42, 901–911. doi:10.1111/j.1365-313x.2005.02423.x
Kim, S. R., Yang, J. I., Moon, S., Ryu, C. H., An, K., Kim, K. M., et al. (2009). Rice OGR1 encodes a pentatricopeptide repeat-DYW protein and is essential for RNA editing in mitochondria. Plant J 59, 738–749. doi:10.1111/j.1365-313X.2009.03909.x
Lei, J., Teng, X., Wang, Y. F., Jiang, X. K., Zhao, H. H., Zheng, X. M., et al. (2022). Plastidic pyruvate dehydrogenase complex E1 component subunit Alpha1 is involved in galactolipid biosynthesis required for amyloplast development in rice. Plant Biotechnol. J. 20, 437–453. doi:10.1111/pbi.13727
Li, H., and Durbin, R. (2010). Fast and accurate long-read alignment with Burrows-Wheeler transform. Bioinformatics 26, 589–595. doi:10.1093/bioinformatics/btp698
Li, J. M., Xiao, J. H., Grandillo, S., Jiang, L. Y., Wan, Y. Z., Deng, Q. Y., et al. (2004). QTL detection for rice grain quality traits using an interspecific backcross population derived from cultivated Asian (O. sativa L.) and African (O. glaberrima S.) rice. Genome 47, 697–704. doi:10.1139/G04-029
Li, Y. B., Fan, C. C., Xing, Y. Z., Yun, P., Luo, L. J., Yan, B., et al. (2014). Chalk5 encodes a vacuolar H+- translocating pyrophosphatase influencing grain chalkiness in rice. Nat. Genet. 46, 398–404. doi:10.1038/ng.2923
Li, Z. F., Wan, J. M., Xia, J. F., and Zhai, H. Q. (2003). Mapping quantitative trait loci underlying appearance quality of rice grains (Oryza sativa L.). Acta Genet. Sin. 30, 251–259. doi:10.3969/j.issn.0379-4172(2003)03-0251-09
Liao, B., and Zhang, G. L. (2015). Research progress on rice chalkiness. Crop Res. 1, 77–83. doi:10.3969/j.issn.1001-5280.2015.01.20
Lipka, A. E., Tian, F., Wang, Q. S., Peiffer, J., Li, M., Bradbury, P. J., et al. (2012). GAPIT: genome association and prediction integrated tool. Bioinformatics 28, 2397–2399. doi:10.1093/bioinformatics/bts444
Liu, J. F., Kui, L. M., Zhu, Z. F., Tan, L. B., Wang, G. J., Li, Q. W., et al. (2007). Identification of QTLs association with processing quality and appearance quality of common wild rice. J. Agicult. Biotechnol. 15, 90–96.
Liu, S. Q., Wang, K. B., Fu, Y. X., Zhu, Q., Xiong, H. B., Nadir, S., et al. (2019). The expression characteristics and function of the new gene OsFMO1 in rice (Oryza sativa L.) flavine monooxygenase. Mol. Plant Breed. 17, 1039–1046. doi:10.13271/j.mpb.017.001039
Livak, K. J., and Schmittgen, T. D. (2001). Analysis of relative gene expression data using real-time quantitative PCR and the 2-ΔΔCt method. Methods 25, 402–408. doi:10.1006/meth.2001.1262
Long, W. H., Dong, B. N., Wang, Y. H., Pan, P. Y., Wang, Y. L., Liu, L. L., et al. (2017). Floury endosperm8, encoding the UDP-glucose pyrophosphorylase 1, affects the synthesis and structure of starch in rice endosperm. J. Plant Biol. 60, 513–522. doi:10.1007/s12374-017-0066-3
Mashiguchi, K., Tanaka, K., Sakai, T., Suqawara, S., Kawaide, H., Natsume, M., et al. (2011). The main auxin biosynthesis pathway in Arabidopsis. Proc. Natl. Acad. Sci. 108, 18512–18517. doi:10.1073/pnas.1108434108
Mikami, I., Aikawa, M., Hirano, H. Y., and Sano, Y. (1999). Altered tissue specific expression at the Wx gene of the opaque mutants in rice. Euphytica 105, 91–97. doi:10.1023/A:1003457209225
Misra, G., Anacleto, R., Badoni, S., Butardo, V. M., Molina, L., Graner, A., et al. (2019). Dissecting the genome-wide genetic variants of milling and appearance quality traits in rice. J. Exp. Bot. 19, 5115–5130. doi:10.1093/jxb/erz256
Misra, G., Badoni, S., Parween, S., Singh, R. K., Leung, H., Ladejobi, O., et al. (2021). Genome-wide association coupled gene to gene interaction studies unveil novel epistatic targets among major effect loci impacting rice grain chalkiness. Plant Biotechnol. J. 19, 910–925. doi:10.1111/pbi.13516
Peng, C., Wang, Y. H., Liu, F., Ren, Y. L., Zhou, K. N., Lv, J., et al. (2014). Floury endosperm6 encodes a CBM48 domain-containing protein involved in compound granule formation and starch synthesis in rice endosperm. Plant J. 77, 917–930. doi:10.1111/tpj.12444
Ren, Y. L., Wang, Y. H., Liu, F., Zhou, K. N., Wan, J. M., et al. (2014). Glutelin precursor accumulation3 encodes a regulator of post-golgi vesicular traffic essential for vacuolar protein sorting in rice endosperm. Plant Cell 26, 410–425. doi:10.1105/tpc.113.121376
Ryoo, N., Yu, C., Park, C. S., Baik, M. Y., Park, I. M., Cho, M. H., et al. (2007). Knockout of a starch synthase gene OsSSIIIa/Flo5 causes white-core floury endosperm in rice (Oryza sativa L.). Plant Cell Rep. 26, 1083–1095. doi:10.1007/s00299-007-0309-8
Sachdeva, S., Singh, R., Maurya, A., Singh, V. K., Singh, U. M., Kumar, A., et al. (2024). Multi-model genome-wide association studies for appearance quality in rice. Front. Plant Sci. 14, 1304388. doi:10.3389/fpls.2023.1304388
Satoh, H., Shibahara, K., Tokunaga, T., Nishi, A., Tasaki, M., Hwang, S. K., et al. (2008). Mutation of the plastidial alpha-glucan phosphory lase gene in rice affects the synthesis and structure of starch in the endosperm. Plant Cell 20, 1833–1849. doi:10.1105/tpc.107.054007
She, K. C., Kusano, H., Koizumi, K., Yamakawa, H., Hakata, M., Imamura, T., et al. (2010). A novel factor floury endosperm 2 is involved in regulation of rice grain size and starch quality. Plant Cell 22, 3280–3294. doi:10.1105/tpc.109.070821
Shi, C. H., Wu, J. G., Lou, X. B., Zhu, J., and Wu, P. (2002). Genetic analysis of transparency and chalkiness area at different filling stage of rice (Oryza sativa L.). Field Crop Res. 76, 1–9. doi:10.1016/S0378-4290(02)00011-4
Shin, J. H., Blay, S., McNeney, B., and Graham, J. (2006). LDheatmap: an R function for graphical display of pairwise linkage disequilibria between single nucleotide polymorphisms. J. Stat. Softw. 16, 1–10. doi:10.18637/jss.v016.c03
Stepanova, A. N., Yun, J., Robles, L. M., Novak, O., He, W., Guo, H., et al. (2011). The Arabidopsis YUCCA1 flavin monooxygenase functions in the indole-3-pyruvic acid branch of auxin biosynthesis. Plant Cell 23, 3961–3973. doi:10.1105/tpc.111.088047
Tabassum, R., Dosaka, T., Ichida, H., Morita, R., Ding, Y. F., Abe, T., et al. (2020). FLOURY ENDOSPERM11-2 encodes plastid HSP70-2 involved with the temperature-dependent chalkiness of rice (Oryza sativa L.) grains. Plant J. 103, 604–616. doi:10.1111/tpj.14752
Tan, Y. F., Xing, Y. Z., Li, J. X., Yu, S. B., Xu, C. G., and Zhang, Q. F. (2000). Genetic bases of appearance quality of rice grains in Shanyou 63, an elite rice hybrid. Theor. Appl. Genet. 101, 823–829. doi:10.1007/s001220051549
Tanaka, N., Fujita, N., Nishi, A., Satoh, H., Hosaka, Y., Ugaki, M., et al. (2004). The structure of starch can be manipulated by changing the expression levels of starch branching enzyme IIb in rice endosperm. Plant Biotechnol. J. 2, 507–516. doi:10.1111/j.1467-7652.2004.00097.x
Teng, X., Zhong, M. S., Zhu, X. P., Wang, C. M., Ren, Y. L., Wang, Y. L., et al. (2019). Floury endosperm16 encoding a NAD-dependent cytosolic malate dehydrogenase plays an important role in starch synthesis and seed development in rice. Plant Biotechnol. J. 17, 1914–1927. doi:10.1111/pbi.13108
Turner, S. D. (2014). qqman: an R package for visualizing GWAS results using QQ and Manhattan plots. BioRxiv. doi:10.1101/005165
Wan, X. Y., Wan, J. M., Weng, J. F., Jiang, L., Bi, J. C., Wang, C. M., et al. (2005). Stability of QTLs for rice grain dimension and endosperm chalkiness characteristics across eight environments. Theor. Appl. Genet. 110, 1334–1346. doi:10.1007/s00122-005-1976-x
Wang, E. T., Wang, J. J., Zhu, X. D., Hao, W., Wang, L. Y., Li, Q., et al. (2008). Control of rice grain filling and yield by a gene with a potential signature of domestication. Nat. Genet. 40, 1370–1374. doi:10.1038/ng.220
Wang, J. C., Xu, H., Zhu, Y., Liu, Q. Q., and Cai, X. L. (2013). OsbZIP58, a basic leucine zipper transcription factor, regulates starch biosynthesis in rice endosperm. J. Exper. Bot. 64, 3453–3466. doi:10.1093/jxb/ert187
Wang, Q. Y., Jiang, L. R., Zheng, J. S., Huang, R. Y., Wang, H. C., and Huang, Y. M. (2011). Major QTLs mapping and analysis for rice grain chalkiness and grain shape traits. Mol. Plant Breed. 5, 547–553. doi:10.3969/mpb.009.000547
Wang, Y. H., Ren, Y. L., Liu, X., Jiang, L., Chen, L. M., Han, X. H., et al. (2010). OsRab5a regulates endomembrane organization and storage protein trafficking in rice endosperm cells. Plant J. 64, 812–824. doi:10.1111/j.1365-313X.2010.04370.x
Wen, L. Y., Fukuda, M., Sunada, M., Ishino, S., Ishino, Y., Okita, T. W., et al. (2015). Guanine nucleotide exchange factor 2 for Rab5 proteins coordinated with GLUP6/GEF regulates the intracellular transport of the proglutelin from the Golgi apparatus to the protein storage vacuole in rice endosperm. J. Exper. Bot. 20, 6137–6147. doi:10.1093/jxb/erv325
Wu, M. M., Ren, Y. L., Cai, M. H., Wang, Y. L., Zhu, S. S., Zhu, J. P., et al. (2019). Rice floury endosperm 10 encodes a pentatricopeptide repeat protein that is essential for the trans-splicing of mitochondrial nad1 intron 1 and endosperm development. New Phytol. 223, 736–750. doi:10.1111/nph.15814
Xie, C. L., Tang, J. P., Yao, Y., and Lu, P. L. (2017). Advances in research on the effects of cultural practices on rice quality. China Rice 23, 13–18. doi:10.3969/j.issn.1006-8082.2017.06.003
Xin, X., Ye, J. H., Yang, Y. Y., Li, R. S., Li, Z., Wang, S., et al. (2022). Genetic diversity analysis and GWAS reveal the adaptive loci of milling and appearance quality of japonica rice (Oryza sativa L.) in Northeast China. J. Integrat. Agr. 21, 1539–1550. doi:10.1016/S2095-3119(21)63701-2
Xu, J. J., Zhang, X. F., and Xue, H. W. (2016). Rice aleurone layer specific OsNF-YB1 regulates grain filling and endosperm development by interacting with an ERF transcription factor. J. Exp. Bot. 67, 6399–6411. doi:10.1093/jxb/erw409
Xue, M. Y., Liu, L. L., Yu, Y. F., Zhu, J. P., Gao, H., Wang, Y. H., et al. (2019). Lose-of-function of a rice nucleolus-localized pentatricopeptide repeat protein is responsible for the floury endosperm14 mutant phenotypes. Rice 12, 100. doi:10.1186/s12284-019-0359-x
Yamamoto, Y., Kamiya, N., Morinaka, Y., Matsuoka, M., and Sazuka, T. (2007). Auxin biosynthesis by the YUCCA genes in rice. Plant Physiol. 143, 1362–1371. doi:10.1104/PP.106.091561
Yan, M. Y., Pan, T., Zhu, Y., Jiang, X. K., Yu, M. Z., Wang, R. Q., et al. (2022). FLOURY ENDOSPERM20 encoding SHMT4 is required for rice endosperm development. Plant Biotechnol. J. 20, 1438–1440. doi:10.1111/pbi.13858
Yang, H., Wang, Y. L., Tian, Y. H., Teng, X., Lv, Z. H., Lei, J., et al. (2023). Rice FLOURY ENDOSPERM22, encoding a pentatricopeptide repeat protein, is involved in both mitochondrial RNA splicing and editing and is crucial for endosperm development. J. Integr. Plant Biol. 65, 755–771. doi:10.1111/jipb.13402
Yang, W. P., Xu, P. K., Zhang, J. C., Zhang, S., Li, Z. W., Yang, K., et al. (2022). OsbZIP60-mediated unfolded protein response regulates grain chalkiness in rice. J. Genet. Genom. 49, 414–426. doi:10.1016/j.jgg.2022.02.002
Yano, M., Okuno, K., Kawak, J., Satoh, H., and Omura, T. (1985). High amylose mutants of rice, Oryza sativa L. Theor. Appl. Genet. 69, 253–257. doi:10.1007/BF00662436
Yi, J., Liu, L., Cao, Y., Li, J., and Mei, M. (2013). Cloning characterization and expression of OsFMO(t) in rice encoding a flavin monooxygenase. J. Genet. 92, 471–480. doi:10.1007/s12041-013-0297-0
Yoshioka, Y., Iwata, H., Tabata, M., Ninomiya, S., and Ohsawa, R. (2007). Chalkiness in rice: potential for evaluation with image analysis. Crop Sci. 47, 2113–2120. doi:10.2135/cropsci2006.10.0631sc
You, X. M., Zhang, W. W., Hu, J. L., Jing, R. N., Cai, Y., Feng, Z. M., et al. (2019). Floury endosperm 15 encodes a glyoxalase I involved in compound granule formation and starch synthesis in rice endosperm. Plant Cell Rep. 38, 345–359. doi:10.1007/s00299-019-02370-9
Yu, M. Z., Wu, M. M., Ren, Y. L., Wang, Y. H., Li, J. F., Lei, C. L., et al. (2021). Rice floury endosperm 18 encodes a pentatricopeptide repeat protein required for 5’ processing of mitochondrial nad5 messenger RNA and endosperm development. J. Integr. Plant Biol. 63, 834–847. doi:10.1111/jipb.13049
Zeng, D. L., Qian, Q., Ruan, L. Q., Teng, S., Kunihiro, Y., Fujimoto, H., et al. (2002). QTL analysis of chalkiness size in three dimensions. Chin. J. Rice Sci. 16, 11–14. doi:10.16819/j.1001-7216.2002.01.003
Zhang, D. P., Wu, J. G., Zhang, Y. J., and Shi, C. H. (2012). Phenotypic and candidate gene analysis of a new floury endosperm mutant (osagpl2-3) in rice. Plant Mol. Biol. Rep. 30, 1303–1312. doi:10.1007/s11105-012-0435-5
Zhang, L., Ren, Y. L., Lu, B. Y., Yang, C. Y., Feng, Z. M., Liu, Z., et al. (2016). Floury endosperm7 encodes a regulator of starch synthesis and amyloplast development essential for peripheral endosperm development in rice. J. Exper. Bot. 67, 633–647. doi:10.1093/jxb/erv469
Zheng, L. N., Zhang, W. W., Liu, S. J., Chen, L. M., Liu, X., Chen, X. G., et al. (2012). Genetic relationship between grain chalkiness, protein content, and paste viscosity properties in a backcross inbred population of rice. J. Cereal Sci. 56, 153–160. doi:10.1016/j.jcs.2012.05.003
Zhong, M. S., Liu, X., Liu, F., Ren, Y. L., Wang, Y. L., Zhu, J. P., et al. (2019). Floury endosperm12 encoding alanine aminotransferase 1 regulates carbon and nitrogen metabolism in rice. J. Plant Biol. 62, 61–73. doi:10.1007/s12374-018-0288-z
Zhou, L. H., Li, X. F., Wang, B. L., and Wang, Y. H. (2016). Characteristics of microclimate in paddy field and its effects on yield and quality of rice. Jiangsu Agric. Sci. 44, 81–85. doi:10.15889/j.issn.1002-1302.2016.05.022
Zhou, L. J., Liu, X., Jiang, L., Zheng, L. N., Chen, L. M., Liu, S. J., et al. (2009). Analysis of QTL and GE effects on PGWC in rice (Oryza sativa L.) using CSSL and BIL populations. Sci. Agric. Sin. 42, 1129–1135. doi:10.3864/j.issn.0578-1752.2009.04.001
Keywords: candidate gene, chalkiness, genome-wide association study, molecular marker, quantitative trait loci, single-nucleotide polymorphism
Citation: Xu Q, Jiang J, Jing C, Hu C, Zhang M, Li X, Shen J, Hai M, Zhang Y, Wang D and Dang X (2024) Genome-wide association mapping of quantitative trait loci for chalkiness-related traits in rice (Oryza sativa L.). Front. Genet. 15:1423648. doi: 10.3389/fgene.2024.1423648
Received: 26 April 2024; Accepted: 07 June 2024;
Published: 10 July 2024.
Edited by:
Joong Hyoun Chin, Sejong University, Republic of KoreaReviewed by:
Tian Li, Chinese Academy of Agricultural Sciences, ChinaSung-Ryul Kim, International Rice Research Institute (IRRI), Philippines
Copyright © 2024 Xu, Jiang, Jing, Hu, Zhang, Li, Shen, Hai, Zhang, Wang and Dang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Xiaojing Dang, dangxj@aaas.org.cn