Deletion of the sugar importer gene OsSWEET1b accelerates sugar starvation-promoted leaf senescence in rice.

Leaf senescence is a combined response of plant cells stimulated by internal and external signals. Sugars acting as signaling molecules or energy metabolites can influence the progression of leaf senescence. Both sugar starvation and accumulation can promote leaf senescence with diverse mechanisms that are reported in different species. Sugars Will Eventually be Exported Transporters (SWEETs) are proposed to play essential roles in sugar transport, but whether they have roles in senescence and the corresponding mechanism are unclear. Here, we functionally characterized a sugar transporter, OsSWEET1b, which transports sugar and promotes senescence in rice (Oryza sativa L.). OsSWEET1b could import glucose and galactose when heterologously expressed in Xenopus oocytes and translocate glucose and galactose from the extracellular apoplast into the intracellular cytosol in rice. Loss of function of OsSWEET1b decreased glucose and galactose accumulation in leaves. ossweet1b mutants showed accelerated leaf senescence under natural and dark-induced conditions. Exogenous application of glucose and galactose complemented the defect of OsSWEET1b deletion-promoted senescence. Moreover, the senescence-activated transcription factor OsWRKY53, acting as a transcriptional repressor, genetically functions upstream of OsSWEET1b to suppress its expression. OsWRKY53-overexpressing plants had attenuated sugar accumulation, exhibiting a similar phenotype as the ossweet1b mutants. Our findings demonstrate that OsWRKY53 downregulates OsSWEET1b to impair its influx transport activity, leading to compromised sugar accumulation in the cytosol of rice leaves where sugar starvation promotes leaf senescence.

Fujian cytoplasmic male sterility and the fertility restorer gene OsRf19 provide a promising breeding system for hybrid rice.

Cytoplasmic male sterility (CMS) determined by mitochondrial genes and restorer of fertility (Rf) controlled by nuclear-encoded genes provide the breeding systems of many hybrid crops for the utilization of heterosis. Although several CMS/Rf systems have been widely exploited in rice, hybrid breeding using these systems has encountered difficulties due to either fertility instability or complications of two-locus inheritance or both. In this work, we characterized a type of CMS, Fujian Abortive cytoplasmic male sterility (CMS-FA), with stable sporophytic male sterility and a nuclear restorer gene that completely restores hybrid fertility. CMS is caused by the chimeric open reading frame FA182 that specifically occurs in the mitochondrial genome of CMS-FA rice. The restorer gene OsRf19 encodes a pentatricopeptide repeat (PPR) protein targeted to mitochondria, where it mediates the cleavage of FA182 transcripts, thus restoring male fertility. Comparative sequence analysis revealed that OsRf19 originated through a recent duplication in wild rice relatives, sharing a common ancestor with OsRf1a/OsRf5, a fertility restorer gene for Boro II and Hong-Lian CMS. We developed six restorer lines by introgressing OsRf19 into parental lines of elite CMS-WA hybrids; hybrids produced from these lines showed equivalent or better agronomic performance relative to their counterparts based on the CMS-WA system. These results demonstrate that CMS-FA/OsRf19 provides a highly promising system for future hybrid rice breeding.

Natural variation of WBR7 confers rice high yield and quality by modulating sucrose supply in sink organs.

Grain chalkiness is an undesirable trait that negatively regulates grain yield and quality in rice. However, the regulatory mechanism underlying chalkiness is complex and remains unclear. We identified a positive regulator of white-belly rate (WBR). The WBR7 gene encodes sucrose synthase 3 (SUS3). A weak functional allele of WBR7 is beneficial in increasing grain yield and quality. During the domestication of indica rice, a functional G/A variation in the coding region of WBR7 resulted in an E541K amino acid substitution in the GT-4 glycosyltransferase domain, leading to a significant decrease in decomposition activity of WBR7A (allele in cultivar Jin23B) compared with WBR7G (allele in cultivar Beilu130). The NIL(J23B) and knockout line NIL(BL130)KO exhibited lower WBR7 decomposition activity than that of NIL(BL130) and NIL(J23B)COM, resulting in less sucrose decomposition and metabolism in the conducting organs. This caused more sucrose transportation to the endosperm, enhancing the synthesis of storage components in the endosperm and leading to decreased WBR. More sucrose was also transported to the anthers, providing sufficient substrate and energy supply for pollen maturation and germination, ultimately leading to an increase rate of seed setting and increased grain yield. Our findings elucidate a mechanism for enhancing rice yield and quality by modulating sucrose metabolism and allocation, and provides a valuable allele for improved rice quality.

Bridge synergy and simplicial interaction in complex contagions.

Modeling complex contagion in networked systems is an important topic in network science, for which various models have been proposed, including the synergistic contagion model that incorporates coherent interference and the simplicial contagion model that involves high-order interactions. Although both models have demonstrated success in investigating complex contagions, their relationship in modeling complex contagions remains unclear. In this study, we compare the synergy and the simplest form of high-order interaction in the simplicial contagion model, known as the triangular one. We analytically show that the triangular interaction and the synergy can be bridged within complex contagions through the joint degree distribution of the network. Monte Carlo simulations are then conducted to compare simplicial and corresponding synergistic contagions on synthetic and real-world networks, the results of which highlight the consistency of these two different contagion processes and thus validate our analysis. Our study sheds light on the deep relationship between the synergy and high-order interactions and enhances our physical understanding of complex contagions in networked systems.

A Ghd7-centered regulatory network provides a mechanistic approximation to optimal heterosis in an elite rice hybrid.

Heterosis refers to the superior performance of hybrids over their parents, which is a general phenomenon occurring in diverse organisms. Many commercial hybrids produce high yield without delayed flowering, which we refer to as optimal heterosis and is desired in hybrid breeding. Here, we attempted to illustrate the genomic basis of optimal heterosis by reinvestigating the single-locus quantitative trait loci and digenic interactions of two traits, the number of spikelets per panicle (SP) and heading date (HD), using recombinant inbred lines and 'immortalized F2 s' derived from the elite rice (Oryza sativa) hybrid Shanyou 63. Our analysis revealed a regulatory network that may provide an approximation to the genetic constitution of the optimal heterosis observed in this hybrid. In this network, Ghd7 works as the core element, and three other genes, Ghd7.1, Hd1, and Hd3a/RFT1, also have major roles. The effects of positive dominance by Ghd7 and Ghd7.1 and negative dominance by Hd1 and Hd3a/RFT1 in the hybrid background contribute the major part to the high SP without delaying HD; numerous epistatic interactions, most of which involve Ghd7, also play important roles collectively. The results expand our understanding of the genic interaction networks underlying hybrid rice breeding programs, which may be very useful in future crop genetic improvement.

The versatile functions of OsALDH2B1 provide a genic basis for growth-defense trade-offs in rice.

In plants, enhanced defense often compromises growth and development, which is regarded as trade-offs between growth and defense. Here we identified a gene, OsALDH2B1, that functions as a master regulator of the growth-defense trade-off in rice. OsALDH2B1 has its primary function as an aldehyde dehydrogenase and a moonlight function as a transcriptional regulator. Loss of function of OsALDH2B1 greatly enhanced resistance to broad-spectrum pathogens, including fungal blast, bacterial leaf blight, and leaf streak, but caused severe phenotypic changes such as male sterility and reduced plant size, grain size, and number. We showed that its primary function as a mitochondrial aldehyde dehydrogenase conditions male fertility. Its moonlight function of transcriptional regulation, featuring both repressing and activating activities, regulates a diverse range of biological processes involving brassinolide, G protein, jasmonic acid, and salicylic acid signaling pathways. Such regulations cause large impacts on the morphology and immunity of rice plants. The versatile functions of OsALDH2B1 provide an example of the genic basis of growth-defense trade-offs in plants.

Investigation on the influence of heterogeneous synergy in contagion processes on complex networks.

Synergistic contagion in a networked system occurs in various forms in nature and human society. While the influence of network's structural heterogeneity on synergistic contagion has been well studied, the impact of individual-based heterogeneity on synergistic contagion remains unclear. In this work, we introduce individual-based heterogeneity with a power-law form into the synergistic susceptible-infected-susceptible model by assuming the synergistic strength as a function of individuals' degree and investigate this synergistic contagion process on complex networks. By employing the heterogeneous mean-field (HMF) approximation, we analytically show that the heterogeneous synergy significantly changes the critical threshold of synergistic strength σc that is required for the occurrence of discontinuous phase transitions of contagion processes. Comparing to the synergy without individual-based heterogeneity, the value of σc decreases with degree-enhanced synergy and increases with degree-suppressed synergy, which agrees well with Monte Carlo prediction. Next, we compare our heterogeneous synergistic contagion model with the simplicial contagion model [Iacopini et al., Nat. Commun. 10, 2485 (2019)], in which high-order interactions are introduced to describe complex contagion. Similarity of these two models are shown both analytically and numerically, confirming the ability of our model to statistically describe the simplest high-order interaction within HMF approximation.

The rice heavy-metal transporter OsNRAMP1 regulates disease resistance by modulating ROS homoeostasis.

Crop diseases threaten food security and sustainable agriculture. Consumption of crops containing nonessential toxic metals leads to health risks for humans. Therefore, cultivation of disease-resistant and toxic metal-safe crops is a double-gain strategy that can contribute to food security. Here, we show that rice heavy-metal transporter OsNRAMP1 plays an important role in plant immunity by modulating metal ion and reactive oxygen species (ROS) homoeostasis. OsNRAMP1 expression was induced after pathogenic bacteria and fungi infections. The osnramp1 mutants had an increased content of H2 O2 and activity of superoxide dismutase, but decreased activity of catalase, showing enhanced broad-spectrum resistance against bacterial and fungal pathogens. RNA-seq analysis identified a number of differentially expressed genes that were involved in metal ion and ROS homoeostasis. Altered expression of metal ion-dependent ROS-scavenging enzymes genes and lower accumulation of cations such as Mn together induced compromised metal ion-dependent enzyme-catalysing activity and modulated ROS homoeostasis, which together contributed towards disease resistance in osnramp1 mutants. Furthermore, the osnramp1 mutants contained lower levels of toxic heavy metals Cd and Pb and micronutrients Ni and Mn in leaves and grains. Taken together, a proof of concept was achieved that broad-spectrum disease-resistant and toxic heavy-metal-safe rice was engineered by removal of the OsNRAMP1 gene.

qFC6, a major gene for crude fat content and quality in rice.

KEY MESSAGE: qFC6, a major quantitative trait locus for rice crude fat content, was fine mapped to be identical with Wx. FC6 negatively regulates crude fat content and rice quality. Starch, protein and lipids are the three major components in rice endosperm. The lipids content in rice influences both storage and quality. In this study, we identified a quantitative trait locus (QTL), qFC6, for crude fat (free lipids) content through association analysis and linkage analysis. Gene-based association analysis revealed that LOC_Os06g04200, also known as Wx, was the candidate gene for qFC6. Complementation and knockout transgenic lines revealed that Wx negatively regulates crude fat content. Lipid composition and content analysis by gas chromatography and taste evaluation analysis showed that FC6 positively influenced bound lipids content and negatively affected both free lipids content and taste. Besides, higher free lipids content rice varieties exhibit more lustrous appearance after cooking and by adding extra oil during cooking could improve rice luster and taste score, indicating that higher free lipids content may make rice more lustrous and delicious. Together, we cloned a QTL coordinating rice crude fat content and eating quality and assisted in uncovering the genetic basis of rice lipid content and in the improvement of rice eating quality.

Fine mapping of qTGW2b and qGL9, two minor QTL conferring grain size and weight in rice.

Rice grain size is a key determinant of both grain yield and quality. In this study, we conducted QTL mapping on grain size using a recombinant inbred line (RIL) population derived from a cross between japonica variety Beilu130 (BL130) and indica variety Jin23B (J23B). A total of twenty-two QTL related to grain length (GL), grain width (GW), grain length-to-width ratio (LWR), grain thickness (GT), and thousand grain weight (TGW) were detected under two environments, and 14 of them were repeatedly detected. Two minor QTL, qTGW2b and qGL9, were validated and further delimited to regions of 631 kb and 272 kb, respectively. Parental sequence comparison of genes expressed in inflorescence in corresponding candidate regions identified frameshifts in the exons of LOC_Os02g38690 and LOC_Os02g38780, both of which encode protein phosphatase 2C-containing protein, and LOC_Os09g29930, which encodes a BIM2 protein. Scanning electron microscopy (SEM) analysis revealed that the increase of cell size rather than cell number caused the differences in grain size between NILs of qTGW2b and qGL9. Quantitative RT-PCR analysis showed that the expression levels of EXPA4, EXPA5, EXPA6, EXPB3, EXPB4, and EXPB7 were significantly different in both qTGW2b NILs and qGL9 NILs. Our results lay the foundation for the cloning of qTGW2b and qGL9, and provide genetic materials for the improvement of rice yield and quality. Supplementary Information: The online version contains supplementary material available at 10.1007/s11032-022-01328-2.

Patterns of genome-wide allele-specific expression in hybrid rice and the implications on the genetic basis of heterosis.

Utilization of heterosis has greatly increased the productivity of many crops worldwide. Although tremendous progress has been made in characterizing the genetic basis of heterosis using genomic technologies, molecular mechanisms underlying the genetic components are much less understood. Allele-specific expression (ASE), or imbalance between the expression levels of two parental alleles in the hybrid, has been suggested as a mechanism of heterosis. Here, we performed a genome-wide analysis of ASE by comparing the read ratios of the parental alleles in RNA-sequencing data of an elite rice hybrid and its parents using three tissues from plants grown under four conditions. The analysis identified a total of 3,270 genes showing ASE (ASEGs) in various ways, which can be classified into two patterns: consistent ASEGs such that the ASE was biased toward one parental allele in all tissues/conditions, and inconsistent ASEGs such that ASE was found in some but not all tissues/conditions, including direction-shifting ASEGs in which the ASE was biased toward one parental allele in some tissues/conditions while toward the other parental allele in other tissues/conditions. The results suggested that these patterns may have distinct implications in the genetic basis of heterosis: The consistent ASEGs may cause partial to full dominance effects on the traits that they regulate, and direction-shifting ASEGs may cause overdominance. We also showed that ASEGs were significantly enriched in genomic regions that were differentially selected during rice breeding. These ASEGs provide an index of the genes for future pursuit of the genetic and molecular mechanism of heterosis.

Heterogeneity induced splay state of amplitude envelope in globally coupled oscillators.

Splay states of the amplitude envelope are stably observed as a heterogenous node is introduced into the globally coupled identical oscillators with repulsive coupling. With the increment of the frequency mismatches between the heterogenous nodes and the rest identical globally coupled oscillators, the formal stable splay state based on the time series becomes unstable, while a splay state based on the new-born amplitude envelopes of time series is stably observed among the rest identical oscillators. The characteristics of the splay state based on the amplitude envelope are numerically and theoretically presented for different parameters of the coupling strength ϵ and the frequency mismatches Δω for small coupling strength and large frequency mismatches. We expect that all these results could reveal the generality of splay states in coupled nonidentical oscillators and help to understand the rich dynamics of amplitude envelopes in multidisciplinary fields.

GLW7.1, a Strong Functional Allele of Ghd7, Enhances Grain Size in Rice.

Grain size is a key determinant of both grain weight and grain quality. Here, we report the map-based cloning of a novel quantitative trait locus (QTL), GLW7.1 (Grain Length, Width and Weight 7.1), which encodes the CCT motif family protein, GHD7. The QTL is located in a 53 kb deletion fragment in the cultivar Jin23B, compared with the cultivar CR071. Scanning electron microscopy analysis and expression analysis revealed that GLW7.1 promotes the transcription of several cell division and expansion genes, further resulting in a larger cell size and increased cell number, and finally enhancing the grain size as well as grain weight. GLW7.1 could also increase endogenous GA content by up-regulating the expression of GA biosynthesis genes. Yeast two-hybrid assays and split firefly luciferase complementation assays revealed the interactions of GHD7 with seven grain-size-related proteins and the rice DELLA protein SLR1. Haplotype analysis and transcription activation assay revealed the effect of six amino acid substitutions on GHD7 activation activity. Additionally, the NIL with GLW7.1 showed reduced chalkiness and improved cooking and eating quality. These findings provide a new insight into the role of Ghd7 and confirm the great potential of the GLW7.1 allele in simultaneously improving grain yield and quality.

Wide Grain 7 increases grain width by enhancing H3K4me3 enrichment in the OsMADS1 promoter in rice (Oryza sativa L.).

Grain size is a major determinant of grain weight, a key component of grain yield of rice. Here, we identified the grain size gene WIDE GRAIN 7 (WG7) from a T-DNA insertion mutant. The grain size of WG7 knockout mutants and WG7 overexpression lines indicated that WG7 is a positive regulator of grain size. WG7 encodes a cysteine-tryptophan (CW) domain-containing transcriptional activator. EMSAs and ChIP-qPCR assay confirmed that WG7 directly bound to the promoter of OsMADS1, a grain size gene, and thereby significantly activated its expression. Point mutations showed that the cis-element CATTTC motif in the promoter was the binding site of WG7. Compared with the wild-type, deletion mutants of the cis-element motif exhibited lower expression of OsMADS1 and produced narrower grains, implicating the requirement of this motif for WG7 function. ChIP-qPCR assays showed that WG7 enhanced histone H3K4me3 enrichment in the promoter of OsMADS1. WG7 underwent directional selection due to the poor fertility of the non-functional mutant. These findings demonstrated that WG7 upregulated OsMADS1 expression by directly binding to its promoter, enhanced histone H3K4me3 enrichment in the promoter and ultimately increased grain width. This study will enrich the knowledge concerning the regulatory network of grain size formation in rice.

An ethyl methanesulfonate-induced neutral mutant-bridging method efficiently identifies spontaneously mutated genes in rice.

Spontaneous mutants are mainly obtained from tissue culture or natural occurrences in plants. The traditional strategy for identifying spontaneously mutated genes is to continuously backcross these mutants to another variety and develop a near-isogenic F2 population for map-based cloning or bulked segregant analysis. However, this strategy is time-consuming. Here, we have developed a new method to efficiently accelerate the identification process. The chemical mutagen ethyl methanesulfonate was first used to treat the wild type of the spontaneous mutants to induce thousands of neutral mutations. An induced individual without any statistically significant phenotypic changes which was compared with the wild type was chosen as the neutral mutant. The spontaneous mutant was then crossed with the neutral mutant to develop a pseudo-near-isogenic F2 population in which only the induced neutral mutations and the causal mutation were segregated in the genome. This population ensures that the variation of the mutated trait is controlled only by the spontaneously mutated gene. Finally, after sequencing the neutral mutant and the mutant-type DNA pool of the F2 population the spontaneous mutation will be identified quickly by bioinformatics analysis. Using this method, two spontaneously mutated genes were identified successfully. Therefore, the neutral mutant-bridging method efficiently identifies spontaneously mutated genes in rice, and its value in other plants is discussed.

OsHOX1 and OsHOX28 Redundantly Shape Rice Tiller Angle by Reducing HSFA2D Expression and Auxin Content.

Tiller angle largely determines plant architecture, which in turn substantially influences crop production by affecting planting density. A recent study revealed that HEAT STRESS TRANSCRIPTION FACTOR2D (HSFA2D) acts upstream of LAZY1 (LA1) to regulate tiller angle establishment in rice (Oryza sativa). However, the mechanisms underlying transcriptional regulation of HSFA2D remain unknown. In this study, two class II homeodomain-Leu zipper genes, OsHOX1 and OsHOX28, were identified as positive regulators of tiller angle by affecting shoot gravitropism. OsHOX1 and OsHOX28 showed strong transcriptional suppressive activity in rice protoplasts and formed intricate self- and mutual-transcriptional negative feedback loops. Moreover, OsHOX1 and OsHOX28 bound to the pseudopalindromic sequence CAAT(C/G)ATTG within the promoter of HSFA2D, thus suppressing its expression. In contrast to HSFA2D and LA1, OsHOX1 and OsHOX28 attenuated lateral auxin transport, thus repressing the expression of WUSCHEL-RELATED HOMEOBOX 6 (WOX6) and WOX11 in the lower side of the shoot base of plants subjected to gravistimulation. Genetic analysis further confirmed that OsHOX1 and OsHOX28 act upstream of HSFA2D Additionally, both OsHOX1 and OsHOX28 inhibit the expression of multiple OsYUCCA genes and decrease auxin biosynthesis. Taken together, these results demonstrated that OsHOX1 and OsHOX28 regulate the local distribution of auxin, and thus tiller angle establishment, through suppression of the HSFA2D-LA1 pathway and reduction of endogenous auxin content. Our finding increases the knowledge concerning fine tuning of tiller angles to optimize plant architecture in rice.

OsVQ1 links rice immunity and flowering via interaction with a mitogen-activated protein kinase OsMPK6.

KEY MESSAGE: Knocking out OsVQ1 in rice released OsMPK6 for activation and in turn promoted H2O2 accumulation, which repressed the expression of flowering-promoting genes, thus delaying rice flowering but enhancing disease resistance. The valine-glutamine (VQ) protein family, which contains the conserved motif FxxxVQxLTG ("x" represents any amino acid), plays a crucial role in plant growth and immunity along with mitogen-activated protein kinase (MAPK) cascades. However, only a few rice VQ proteins have been functionally characterized, and the roles of the MAPK-VQ module in rice biological processes are not fully understood. Here, we investigated the role of OsVQ1 in rice disease resistance and the control of flowering time. The OsVQ1-knock out (KO) mutants exhibited increased resistance to Xanthomonas oryzae pathovars, accumulated high levels of hydrogen peroxide (H2O2), and showed a late flowering phenotype under natural long-day conditions, while the OsVQ1-overexpressing plants showed phenotypes similar to that of the wild type. Further studies revealed that OsVQ1 physically interacted with and inhibited OsMPK6 activity. In addition, OsVQ1 expression was downregulated by the pathogen-induced OsMPKK10.2-OsMPK6-OsWRKY45 cascade, suggesting a feedback loop between OsVQ1 and OsMPK6. Moreover, the OsVQ1-KO/osmpk6 double-mutant exhibited increased susceptibility to X. oryzae infection and showed an early flowering phenotype, which may partially be attributed to the reduced accumulation of H2O2 and the consequent up-expression of flowering-promoting genes. These results suggested that the OsVQ1-OsMPK6 module was involved in rice immunity and flowering.

Transfer learning of chaotic systems.

Can a neural network trained by the time series of system A be used to predict the evolution of system B? This problem, knowing as transfer learning in a broad sense, is of great importance in machine learning and data mining yet has not been addressed for chaotic systems. Here, we investigate transfer learning of chaotic systems from the perspective of synchronization-based state inference, in which a reservoir computer trained by chaotic system A is used to infer the unmeasured variables of chaotic system B, while A is different from B in either parameter or dynamics. It is found that if systems A and B are different in parameter, the reservoir computer can be well synchronized to system B. However, if systems A and B are different in dynamics, the reservoir computer fails to synchronize with system B in general. Knowledge transfer along a chain of coupled reservoir computers is also studied, and it is found that, although the reservoir computers are trained by different systems, the unmeasured variables of the driving system can be successfully inferred by the remote reservoir computer. Finally, by an experiment of chaotic pendulum, we demonstrate that the knowledge learned from the modeling system can be transferred and used to predict the evolution of the experimental system.

The rice Raf-like MAPKKK OsILA1 confers broad-spectrum resistance to bacterial blight by suppressing the OsMAPKK4-OsMAPK6 cascade.

Mitogen-activated protein kinase kinase kinase (MAPKKK) are the first components of MAPK cascades, which play pivotal roles in signaling during plant development and physiological processes. The genome of rice encodes 75 MAPKKKs, of which 43 are Raf-like MAPKKKs. The functions and action modes of most of the Raf-like MAPKKKs, whether they function as bona fide MAPKKKs and which are their downstream MAPKKs, are largely unknown. Here, we identified the osmapkkk43 mutant, which conferred broad-spectrum resistance to Xanthomonas oryzae pv. oryzae (Xoo), the destructive bacterial pathogen of rice. Oryza sativa (Os)MAPKKK43 encoding a Raf-like MAPKKK was previously known as Increased Leaf Angle 1 (OsILA1). Genetic analysis indicated that OsILA1 functioned as a negative regulator and acted upstream of the OsMAPKK4-OsMAPK6 cascade in rice-Xoo interactions. Unlike classical MAPKKKs, OsILA1 mainly phosphorylated the threonine 34 site at the N-terminal domain of OsMAPKK4, which possibly influenced the stability of OsMAPKK4. The N-terminal domain of OsILA1 is required for its homodimer formation and its full phosphorylation capacity. Taken together, our findings reveal that OsILA1 acts as a negative regulator of the OsMAPKK4-OsMAPK6 cascade and is involved in rice-Xoo interactions.

The origin of Wxla provides new insights into the improvement of grain quality in rice.

Appearance and taste are important factors in rice (Oryza sativa) grain quality. Here, we investigated the taste scores and related eating-quality traits of 533 diverse cultivars to assess the relationships between-and genetic basis of-rice taste and eating-quality. A genome-wide association study highlighted the Wx gene as the major factor underlying variation in taste and eating quality. Notably, a novel waxy (Wx) allele, Wxla , which combined two mutations from Wxb and Wxin , exhibited a unique phenotype. Reduced GBSSI activity conferred Wxla rice with both a transparent appearance and good eating quality. Haplotype analysis revealed that Wxla was derived from intragenic recombination. In fact, the recombination rate at the Wx locus was estimated to be 3.34 kb/cM, which was about 75-fold higher than the genome-wide mean, indicating that intragenic recombination is a major force driving diversity at the Wx locus. Based on our results, we propose a new network for Wx evolution, noting that new Wx alleles could easily be generated by crossing genotypes with different Wx alleles. This study thus provides insights into the evolution of the Wx locus and facilitates molecular breeding for quality in rice.