Coordinating High Yield and Superior Eating Quality of Rice: A Case Study of Hybrid Varieties Derived from Longke638S and Jing4155S

: The planting area of hybrid rice ( Oryza sativa L.) has continuously decreased in recent years partially because of the low eating quality. Longke638S and Jing4155S are two elite male sterile lines, which are used for developing many hybrid varieties with large planting areas in China. Here, 21 hybrid rice varieties of an incomplete diallel cross population with Longke638S and Jing4155S as female parents were planted under ﬁeld conditions for two consecutive years, aiming to investigate the physiological and transcriptomic characteristics that is required for coordinating high yield and superior eating quality. As a result, grain yield ranged from 7.15 to 9.23 t ha − 1 in 2019 and from 7.49 to 9.23 t ha − 1 in 2020, and a negative relationship was identiﬁed between yield and quality-related traits. Hybrid varieties with Longke638S as the female parent (LLY) had higher yield but lower eating quality than those with Jing4155S as the female parent (JLY). The higher yield of LLY over JLY was due to higher biomass production and N uptake as well as higher yield stability across planting years, while their lower grain quality can be mainly ascribed to lower amylose content and alkali spreading value. Male parents signiﬁcantly affected yield-related traits of LLY and eating quality-related traits of JLY. Candidate genes involved in grain starch biosynthesis were differentially expressed between LLYHZ, JLYHZ, and JLY1212, such as Wx and ALK . Overall, this study examined the relationship between grain yield and eating quality and identiﬁed the physiological and molecular traits that limited grain yield or eating quality of LLY and JLY hybrid varieties.


Introduction
In addition to the first Green Revolution, heterosis has been utilized to significantly improve the yield potential of various crops, which has contributed greatly to food security worldwide [1]. The grain yield of hybrid rice (Oryza sativa L.) is higher than that of inbred rice by approximately 20% in China and by nearly 9% in tropical irrigated lowlands [2,3]. Development of super hybrid rice by combining intersubspecific heterosis and the ideotype approach further increases the yield potential to more than 15 t ha −1 [4]. From 1976 to 2014, hybrid rice accounted for 50-60% of the total rice planting area in China, and increased rice production by 500 billion kg in total based on a 15% higher grain yield of hybrid rice over inbred rice on average [5]. However, despite the advantage of hybrid rice in yield, the planting area of hybrid rice has been continuously decreasing in recent years, which may be partially ascribed to its relatively lower grain quality [5]. Therefore, it is essential to diallel cross population with Longke638S and Jing4155S as female parents were planted at a medium N fertilizer condition for two consecutive years. This study aims to (1) determine the differences in grain yield and eating quality of these elite hybrid varieties across two planting years, and (2) explore the physiological and molecular attributes that can be used for tailoring rice varieties coordinating high grain yield and superior quality.

Site Description
The experiments were conducted in farmers' fields at Maoshi Township, Jianli County, Hubei Province, China (29 • 51 N, 115 • 53 E) during the rice growing season (May to October) in 2019-2020. This area is located in a subtropical monsoon humid climate region. Every year, the daily meteorological data were collected during the rice growing period between transplanting and maturity (based on the test variety with the longest duration) from a weather station (AWS 800, Campbell Scientific, Inc., Logan, UT, USA) located near the experimental site. The seasonal average daily mean temperature, average daily solar radiation, and total precipitation were 27.8 • C, 16.3 MJ m −2 , and 214.4 mm in 2019, respectively, and were 25.7 • C, 11.9 MJ m −2 , and 629.1 mm in 2020, respectively. In general, 2019 was a hotter year than 2020 ( Figure S1). The soil properties before the experiment in 2019 were as follows: pH 7.27, 21.72 g kg −1 organic matter, 1.63 g kg −1 total N, 7.34 mg kg −1 Olsen P, 100.61 mg kg −1 exchangeable K. The soil was sampled at a depth of 20 cm from the surface before the application of basal fertilizer.

Experimental Design and Crop Management
The experiments were conducted in a randomized block design with three replications and a plot size of 20.0 m 2 . Pre-germinated seeds were sown into seedbeds on the 15th and 16th May in 2019 and 2020, respectively. Transplanting was performed at a density of 25 hills m −2 and a hill spacing of 13.3 × 30.0 cm with two seedlings per hill on 16 and 17 June in 2019 and 2020, respectively. In the two years, a total N rate of 100 kg ha −1 was applied as basal, at tillering (10 days after transplanting (DAT)) and panicle initiation (PI) stage at a ratio of 4:2:4. All N fertilizer was applied in the form of urea. The N rate in the present study was around 45% lower than the average N rate of 180 kg N ha −1 in central China to reduce the environmental impact and increase grain quality while maintaining a high grain yield [27]. Phosphorus (40 kg P ha −1 as calcium superphosphate) was applied in all plots one day before transplanting. Potassium (100 kg K ha −1 as potassium chloride) was split equally and applied in all plots at basal and PI. The field was flooded after transplanting, and a floodwater depth of 3-5 cm was maintained until a week before maturity except for water draining at the maximum tillering stage to reduce unproductive tillers. Weeds, pests, and diseases were intensively controlled to avoid yield loss.
The tested varieties were from an incomplete diallel cross population. The female parents were Longke638s and Jing4155S, and the male parents were 12 restorer lines (Table S1). In 2019, 21 hybrid varieties and Fengliangyou4 (FLY4) were planted, while in 2020, 12 restorer lines and Huanghuazhan (HHZ) were added to the experiment. Fengliangyou4 was used as the control variety as it is commonly used as the control in official variety certification in the Middle and Lower Reaches of Yangtze River. HHZ is a widely planted inbred variety in this area. Among these varieties, LLYHZ, LLY1377, LLY1988, LLY1212, JLYHZ, and JLY1212 are designated as super hybrid rice by the Ministry of Agriculture and Rural Affairs (MARA). JLY1212 has superior grain quality (especially taste quality), which meets the standards (NY/T593-2013) of level 2 high-quality rice when planted in this region (China Rice Data Center, https://www.ricedata.cn/, 1 March 2022) [28]. LLYHZ and JLYHZ rank top 5, and JLY1212 ranks top 10 in planting area among all the hybrid varieties of China in recent years.

Plant Sampling and Measurement
Twelve hills were sampled from each plot at the heading (HD) and maturity (MA) stage. The plant samples were then separated into leaves, stems (culm plus sheath), and panicles (when present). The dry biomass at each stage was determined after oven drying at 80 • C to a constant weight. At MA, the panicles were hand threshed after recording the panicle number, and then filled spikelets were separated from unfilled ones by submerging in tap water. Empty spikelets were separated from partially filled spikelets by winnowing. Three subsamples of filled spikelets (30 g for each subsample), three subsamples of empty spikelets (2 g for each subsample), and all partially filled spikelets were considered when recording the number of spikelets. The dry biomass of the rachis, the filled, partially filled, and empty spikelets were measured after oven drying at 80 • C to a constant weight. The dry biomass of the above tissues was added with that of leaves and stems to calculate the total aboveground dry biomass at MA. The number of spikelets per panicle (spikelets m −2 /panicles m −2 ), grain filling percentage (100 × filled spikelets m −2 /spikelets m −2 ), and harvest index (HI) (100 × filled spikelet weight/total dry biomass) were calculated. The grain yield was determined from a 5-m 2 area in the center of each plot at 14% moisture content as determined by a digital moisture tester (DMC-700, Seedburo, Chicago, IL, USA).
The N concentration in each plant part at all growth stages was measured by Elementar vario MAX CNS/CN (Elementar Trading Co., Ltd., Langenselbold, Germany). The N content in each organ was calculated as the product of N concentration and dry biomass. The nitrogen uptake at each growth stage was calculated as the sum of N contents in each organ. The grain protein content was determined by multiplying the N concentration with 5.95 [29]. The nitrogen harvest index (NHI) was calculated as the ratio of N accumulation in grain over the total aboveground N uptake. The nitrogen use efficiency for grain production (NUEg) was the ratio of grain yield over the total aboveground N uptake. The nitrogen use efficiency for biomass production (NUEb) was the ratio of total aboveground biomass over the total aboveground N uptake.

Milling and Appearance Grain Quality
Rice grain quality was assessed after 3 months of storage in a sample drying room, where the temperature was approximately 20 • C and the humidity was approximately 20% for stabilizing the physical and chemical properties of the samples. The milling quality of rice grains was determined according to the national standards of quality evaluation GB/T17891-2017, the People's Republic of China, including the brown rice, milled rice, and head rice rate.
The scanned images of head rice were obtained by the Rice Appearance Quality Scanner (Epson Expressing 1680 Professional, Epson America, Lnc., Los Alamitos, CA, USA). Then, an image analysis software (Image J, the National Institutes of Health, Bethesda, MD, USA) was used to analyze the scanned images based on color to identify the number and area of chalky parts of rice grains and the length and width of head rice. After that, the chalky grain rate and chalkiness were calculated.

Determination of Amylose Content
The amylose content (AC) was measured according to the GB/T 17891-2017.

Determination of Gel Consistency
Gel consistency (GC) was determined according to the GB/T 17891-2017. Briefly, test tubes (13 mm × 100 mm) were filled with 0.01 g of sample rice flour in duplicate, followed by the addition of 0.2 mL 0.025% thymol blue solution. The tubes were then gently shaken and added with 2.0 mL 0.2 mol L −1 potassium hydroxide solution, and then shaken again. Then, the test tube was covered with a glass marble, put in a boiling water bath, and heated for 8 min. The rice glue solution surface in the test tube was controlled at one-half to two-thirds of the height of the test tube during the heating process. The tubes were then taken out, and after the removal of the glass marbles, the tubes were allowed to cool for 5 min and then put in an ice water bath at about 0 • C to cool for 20 min. Subsequently, the tubes were immediately taken out and placed horizontally on an operating table covered with coordinate paper, and let to stand at room temperature (25 ± 2 • C) for 1 h. The length of the rice glue flowing in the test tube (mm) was measured. When the difference between double test results did not exceed 7 mm, the average value was taken as the test result.

Determination of Gelatinization Temperature
The gelatinization temperature was determined indirectly by measuring the alkali spreading value (ASV). Six head rice grains were evenly placed in a covered square box, and 1.7% KOH solution was added. For better dispersion, enough spacing was set between the head rice grains. After capping, the head rice grains were placed at 30 ± 2 • C incubator. After 23 h, the gelatinization state of the rice grains was recorded one by one according to the 7-level classification of alkali spreading value, and the average value was taken as the test result.

RNA Sequencing
Panicles at the flowering stage were sampled and frozen in liquid nitrogen immediately, and then stored at -80 • C in a refrigerator. The total RNA was isolated from the spikelets using a Quick RNA Isolation Kit (Huayueyang, Beijing, China) according to the manufacturer's instructions. The stranded RNA libraries were constructed using the Hieff NGS Ultima Dual-mode mRNA Library Prep Kit for MGI (Yeasen, Shanghai, China) following the manufacturer's recommendations. After library preparation, quality control (QC) and pooling, pair-end sequencing was performed on the MGISEQ-2000RS (MGI, Shenzhen, China) at the Sequencing Platform of the National Key Laboratory of Crop Genetic Improvement in Huazhong Agricultural University (Wuhan, China). The sequence quality of RNA-Seq libraries was evaluated using FastQC, and the adapter sequences and low-quality reads were filtered using Fastp (v0.20.0). Clean reads were aligned against the Oryza sativa Nipponbare reference genome using HISAT2 (v2.2.1) with default parameters [30,31]. Gene expression was quantified using StringTie (v2.1.4) [32]. Differentially expressed genes (DEGs) were identified with the DESeq2 package in R with an adjusted p value < 0.05 and log10 fold-change > 1 [33]. Annotation of DEGs was based on the information from funRiceGenes and MsuNote (MSU v7).

Statistical Analysis
The statistical analysis was performed using analysis of variance (Statistix 9.0, Analytical software, Tallahassee, FL, USA). When significant, the mean values were compared based on the least significant difference (LSD) test at the 0.05 probability level. The Partial Least Squares-Discriminant Analysis (PLS-DA) was performed using the plsr function provided by R pls package [34]. The classification and cross-validation were performed using the corresponding wrapper function offered by the caret package [35]. Figures were generated by Sigmaplot 12.5 (SPSS Inc., Point Richmond, CA, USA).

Grain Yield, Grain Quality and Nitrogen Use Efficiency of Hybrid Varieties
The PLS-DA of traits related with grain yield, grain quality, and NUE indicated that the varieties with Longke638S as the female parent (LLY), those with Jingke4155S as the female parent (JLY), and the control variety (FLY4) were separated along with component 1 ( Figure S2). The detail differences among LLY, JLY, and FLY4 are shown in Figure 1 and Tables 1-3. In terms of average value across the 2 years, LLY exhibited higher (by 5.77%) grain yields than JLY. The higher grain yield of LLY over JLY was mainly due to the higher total dry weight and total nitrogen uptake (Tables 3 and S3; Figure 1). Generally, the test hybrid varieties had higher head rice rate and lower chalkiness, ASV, and AC than the Agronomy 2022, 12, 1628 6 of 16 control FLY4, while LLY had lower ASV and AC but higher GC than JLY ( Figure 1; Table 2). The hybrid varieties with significantly higher yield than FLY4 maintained a high NUEg when N uptake increased. Moreover, higher total N uptake contributed to the higher grain yield of LLY compared with JLY ( Figure 1; Table 3).
The PLS-DA of traits related with grain yield, grain quality, and NUE indicated that the varieties with Longke638S as the female parent (LLY), those with Jingke4155S as the female parent (JLY), and the control variety (FLY4) were separated along with component 1 ( Figure S2). The detail differences among LLY, JLY, and FLY4 are shown in Figure 1 and Tables 1-3. In terms of average value across the 2 years, LLY exhibited higher (by 5.77%) grain yields than JLY. The higher grain yield of LLY over JLY was mainly due to the higher total dry weight and total nitrogen uptake (Tables 3 and S3; Figure 1). Generally, the test hybrid varieties had higher head rice rate and lower chalkiness, ASV, and AC than the control FLY4, while LLY had lower ASV and AC but higher GC than JLY ( Figure 1; Table 2). The hybrid varieties with significantly higher yield than FLY4 maintained a high NUEg when N uptake increased. Moreover, higher total N uptake contributed to the higher grain yield of LLY compared with JLY ( Figure 1; Table 3).  Figure S2.   Figure S2.    For all the varieties, grain yield ranged from 7.15 to 9.23 t ha −1 in 2019, and from 7.49 to 9.23 t ha −1 in 2020 (Table 1). LLYHZ exhibited stably high grain yields across the 2 years with a coefficient of variance (CV) of 0.85%, while the yield of JLY1212 had the highest CV (11.72%). JLYHZ had the highest grain yield among JLY hybrids, which was comparable to that of LLYHZ in 2020 but significantly lower than that of LLYHZ in 2021. Three parameters related to cooking and eating quality exhibited significant variations among these hybrid varieties ( Table 2). LLYHZ had relatively lower ASV and AC but higher protein content. JLYHZ and JLY1212 had the highest ASV and AC among the test hybrids except for FLY4, and JLY1212 also had the lower chalkiness degree and protein content in 2021. Overall, LLYHZ had the highest yield potential and yield stability, JLYHZ had medium yield and quality, while JLY1212 had the best eating quality but lower yield.
N uptake at maturity ranged from 118.6 kg ha −1 to 159.3 kg ha −1 in 2019 and from 118.2 kg ha −1 to 148.5 kg ha −1 in 2020. Most of the LLY hybrid varieties had higher N uptake than the control variety FLY4. Averaged across 2 years, LLYHZ had the highest TN, while LLY8612 had the highest NUEg. Among JLY hybrid varieties, JLYHZ had the highest TN and JLY1206 had the highest NUEg (Table 3).

Correlation between Yield and Quality-Related Traits for Hybrid Varieties
There were significant correlations between yield, NUE, and quality-related parameters ( Figure 2). Grain yield negatively correlated with ASV at p < 0.01 but positively correlated with biomass, HI, TN, NHI, and GC at p < 0.05. For yield components, grain weight positively correlated with chalkiness but negatively correlated with ASV; grain filling percentage positively correlated with NHI, spikelets m −2 and spikelets per panicle positively correlated with NUEg but negatively correlated with protein content and chalkiness degree; and panicle No. positively correlated with ASV but negatively correlated with GC and chalkiness degree. Total N uptake negatively correlated with AC and ASV at p < 0.01, while NUEg positively correlated with AC at p < 0.01. Protein content negatively correlated with AC at p < 0.01. age positively correlated with NHI, spikelets m −2 and spikelets per panicle positively correlated with NUEg but negatively correlated with protein content and chalkiness degree; and panicle No. positively correlated with ASV but negatively correlated with GC and chalkiness degree. Total N uptake negatively correlated with AC and ASV at p < 0.01, while NUEg positively correlated with AC at p < 0.01. Protein content negatively correlated with AC at p < 0.01.

Figure 2.
Correlations between parameters related with grain yield, grain quality, and nitrogen use efficiency. * and ** represent significance at p < 0.05 and p < 0.01, respectively. GF: grain filling percentage; GW: grain weight; GY: grain yield; HI: harvest index; TN: Total N uptake; NUEg: nitrogen use efficiency for grain production; NUEb: nitrogen use efficiency for biomass production; ASV: alkali spreading value; GC: gel consistency; AC: amylose content; PC: protein content.

Grain Yield, Grain Quality, and Nitrogen Use Efficiency of the Male Parents
There were significant differences in parameters related with grain yield, quality, and NUE among the male parents in 2020 (Table 4). There was no significant difference in grain yield between the male parents and HHZ except for R1125, R1377, R3189, and R149. The difference in grain yield of these male parents was mainly due to variation in the total dry weight. Both N uptake and NUEg contributed to the variation in grain yield among the male parents, for example R1206 had the highest N uptake of 139.0 kg ha −1 , while HZ had the highest NUEg of 57.8 kg kg −1 . R1212 had medium N uptake and lower NUEg. Correlations between parameters related with grain yield, grain quality, and nitrogen use efficiency. * and ** represent significance at p < 0.05 and p < 0.01, respectively. GF: grain filling percentage; GW: grain weight; GY: grain yield; HI: harvest index; TN: Total N uptake; NUEg: nitrogen use efficiency for grain production; NUEb: nitrogen use efficiency for biomass production; ASV: alkali spreading value; GC: gel consistency; AC: amylose content; PC: protein content.

Grain Yield, Grain Quality, and Nitrogen Use Efficiency of the Male Parents
There were significant differences in parameters related with grain yield, quality, and NUE among the male parents in 2020 (Table 4). There was no significant difference in grain yield between the male parents and HHZ except for R1125, R1377, R3189, and R149. The difference in grain yield of these male parents was mainly due to variation in the total dry weight. Both N uptake and NUEg contributed to the variation in grain yield among the male parents, for example R1206 had the highest N uptake of 139.0 kg ha −1 , while HZ had the highest NUEg of 57.8 kg kg −1 . R1212 had medium N uptake and lower NUEg.
The test hybrid varieties had significantly lower AC than the control inbred variety HHZ except for R3189. R1199, R1988, R8612, R149, and HZ had lower ASV than the other varieties (Table 4). Correlations between male parents and hybrids in yield and quality related parameters were different depending on the female parent ( Figure 3). Quality-related parameters for LLY hybrid varieties were more related with that of the male parents (Figure 3a-e). However, the male parents had larger effects on yield variation in JLY hybrid varieties (R 2 = 0.40; Figure 3f).

Transcriptomic Analysis of LLYHZ, JLYHZ and JLY1212
Based on the results of yield, NUE, and quality, LLYHZ with high and stable yield JLYHZ with medium yield and eating quality, and JLY1212 with premium eating quality were selected for RNA sequencing (RNA-Seq) to investigate the molecular mechanisms underlying the differences in grain yield and quality (Table S4). The up-regulated DEGs of LLYHZ relative to JLYHZ were mostly related to carbohydrate metabolic process, pho tosynthesis (especially light reaction), and ion transport ( Figure 4A), while those relative to JLY1212 were associated with cell wall organization or biogenesis, carbohydrate meta bolic process, cell death and stress response ( Figure 4B). Most of the up-regulated genes in JLYHZ compared with JLY1212 were involved in response to stress and stimulus ( Figure 4C

Transcriptomic Analysis of LLYHZ, JLYHZ and JLY1212
Based on the results of yield, NUE, and quality, LLYHZ with high and stable yield, JLYHZ with medium yield and eating quality, and JLY1212 with premium eating quality were selected for RNA sequencing (RNA-Seq) to investigate the molecular mechanisms underlying the differences in grain yield and quality (Table S4). The up-regulated DEGs of LLYHZ relative to JLYHZ were mostly related to carbohydrate metabolic process, photosynthesis (especially light reaction), and ion transport ( Figure 4A), while those relative to JLY1212 were associated with cell wall organization or biogenesis, carbohydrate metabolic process, cell death and stress response ( Figure 4B). Most of the up-regulated genes in JLYHZ compared with JLY1212 were involved in response to stress and stimulus ( Figure 4C).  The DEGs with known functions and particularly associations with yield and eating quality were screened out, and the overlapping DEGs between LLYHZ and JLYHZ or JLY1212 are shown in Table S5. The DEGs related to growth and yield were involved in root growth, plant height, nitrogen assimilation, plant type (tiller and leaf angle), and panicle formation (Table S5). The sixteen DEGs related to panicle development can be divided into two categories, namely grain quality and yield components. Most of the DEGs related to grain quality including Wx and ALK were down-regulated, while those related to spikelet number and grain size were mostly up-regulated in LLYHZ (Table 5).  The DEGs with known functions and particularly associations with yield and eating quality were screened out, and the overlapping DEGs between LLYHZ and JLYHZ or JLY1212 are shown in Table S5. The DEGs related to growth and yield were involved in root growth, plant height, nitrogen assimilation, plant type (tiller and leaf angle), and panicle formation (Table S5). The sixteen DEGs related to panicle development can be divided into two categories, namely grain quality and yield components. Most of the DEGs related to grain quality including Wx and ALK were down-regulated, while those related to spikelet number and grain size were mostly up-regulated in LLYHZ (Table 5).

Physiological and Molecular Traits Limiting the Grain Yield of Hybrid Varieties Derived from Jing4155S
In the present study, grain yield and quality were negatively related to each other, and generally, varieties with Longke638S as the female parent had higher yield but lower eating quality compared with those with Jing4155S as the female parent (Tables 1 and 2; Figures 1 and 2). The higher yield resulted from higher dry matter accumulation and higher N uptake (Figures 1 and 2). These results are consistent with the previous finding that genetic improvement in yield potential of rice is associated with concomitant improvement of dry matter accumulation, radiation use efficiency, and N use efficiency [13,36]. Through transcriptomic analysis of three representative super hybrid varieties, namely LLYHZ, JLYHZ, and JLY1212 (Tables 5 and S5), it was found that the larger panicle size of LLYHZ was possibly related to the down-regulated expression of OsHAP3E, which is a negative regulator of panicle and spikelet development [37], and up-regulated expression of positive regulators such as OsFOR1 [38], OsYABBY5 [39], OsEBS [40], and OsJAR1 [41]. Crop yield is largely dependent on N uptake and utilization, which are significantly affected by root architecture and morphology [42]. Here, we found differences in the expression of several genes related to root elongation between LLYHZ and JLYHZ or JLY1212 (Table S5). Moreover, OsNiR, which promotes N uptake and signaling, was also upregulated in LLYHZ (Table S5). Moreover, DEGs related to photosynthesis, especially the light reaction, were up-regulated in LLYHZ compared with JLYHZ ( Figure 4A), which might be related with the higher biomass of LLYHZ.
Yield stability across different planting years is another important factor limiting the yield of hybrid varieties derived from Jing4155S. For example, JLYHZ had the highest yield stability while JLY1212 had the lowest yield stability across years (Table 1). Previous studies have demonstrated that variations in grain yield across different planting years are mainly associated with the variations in temperature in central China [29]. In this study, OsHSBP2, a negative regulator of heat shock response [43], was down-regulated in LLYHZ compared with JLYHZ and JLY1212 (Table S5). Moreover, many other genes related to drought, salinity, and blast resistance were also differentially expressed between LLYHZ and JLYHZ or JLY1212 (Figure 4; Table S5). These molecular characteristics may be utilized to improve the yield stability of high-quality varieties across different environmental and management conditions.

Physiological and Molecular Traits for Improving Grain Quality of Hybrid Varieties Derived from Longke638S
In the present study, compared with FLY4, most of the varieties had higher head rice rate and lower chalkiness, demonstrating genetic improvement in milling and appearance quality ( Table 2). Hybrid varieties derived from Jing4155S had higher ASV and AC, but lower GC than that derived from Longke638S ( Figure 1; Table 2), which resulted from the change in gene expression involved in starch biosynthesis [17,18]. Wx and ALK are the major genes determining AC and ASV, respectively, and they both had minor effects on GC [18,44]. Both of the two genes were down-regulated in LLYHZ compared with JLYHZ and JLY1212 ( Table 5). Mutation of OsPPDKB also resulted in down-regulation of Wx expression, and then reduced the AC [45]. Some other genes involved in starch biosynthesis, including Pho1, OsAGPL2, OsSSIIIa, and OsISA1 were also down-regulated in LLYHZ (Table 5) [17,46]. These traits may serve as potential breeding targets for further improving the eating quality of hybrid varieties derived from Longke638S. For example, Huang et al. [47] found that coordinated regulation of amylose and amylopectin synthesis via manipulation of SSII-2, SSII-3, and Wx expression in endosperm can improve rice eating and cooking quality. Eating quality-related traits of hybrid varieties with Longke638S as the female parent were affected by the male parent significantly ( Figure 3); therefore, selecting a male parent with superior eating quality is essential in hybrid breeding using Longke638S.

Coordination of High Yield and Superior Quality
Balancing grain yield, quality, and stress (abiotic or biotic) tolerance is a great challenge for rice sustainable production in a continuously changing climate scenario. In the present study, LLYHZ, JLYHZ, JLY1212, and several other hybrids are all super hybrid varieties designated by MARA; moreover, JLY1212 had superior eating quality (level 2 in Cooking Rice Variety Quality, NY/T 593-2013). However, it was found that hybrids derived from Longke638S including LLYHZ had high yield potential and stability, but lower grain quality mainly due to the lower AC and ASV. Moreover, ASV was negatively related with yield and yield components such as grain weight, spikelets per panicle, total dry weight and N uptake, and AC negatively correlated with total N uptake and protein content (Figure 2). Wx and ALK are the major genes determining AC and ASV, respectively; however, the effects of natural allelic variation in Wx and ALK on grain yield and the physiological mechanisms deserve further investigation. Although JLY1212 and JLYHZ were super hybrid varieties with high yield potential, both of them had lower grain yield stability. This was mainly due to the lower N uptake and biomass production at a low N fertilizer rate in present study (Tables 3 and S3). Therefore, it is important to improve N uptake capacity of Jing4155S or select restorer lines with high N uptake and biomass production for hybrid breeding with Jing4155S.
The frequently observed negative relationships among grain yield, quality, and stress tolerance are ascribable to the competition in resource allocation and deleterious genetic linkage, since each of the three complex traits is controlled by a number of quantitative genes, some of which exhibit pleiotropic effects [26,48]. In this study, some DEGs between high-yield and high-quality varieties were also found to have pleiotropic effects on several traits related to growth or stress tolerance (Tables 5 and S5). By taking advantage of the allelic variation in qGW7 (a QTL for rice grain quality), Wang et al. [24] realized the improvement of rice yield and quality simultaneously. Xiao et al. [26] also broke the deleterious linkage between yield, taste quality, and blast resistance through the identification of several superior alleles for these traits. Therefore, more studies are required to identify and pyramiding the superior alleles for eating quality-related traits, such as Wx and ALK [19].

Conclusions
Elucidation of the physiological and molecular traits related to high yield and superior quality will facilitate the improvement of the widely planted hybrid rice varieties. In the present study, these traits were investigated using an incomplete diallel cross population with Longke638S and Jing4155S as female parents. The results suggest that for rice hybrid varieties with superior grain quality derived from Jing4155S, sink size (panicle and grain size), biomass production (N uptake and optimized plant type), and yield stability across planting years are essential physiological traits that may serve as the targets to improve their yields; while for high yield hybrid varieties derived from Longke638S, grain quality improvement may be achieved by optimizing the amylose content, ASV, and protein content. These physiological targets can be tailored simultaneously through the pyramiding of superior genes for these traits, which still require more extensive studies with more germplasm resources.
Supplementary Materials: The following supporting information can be downloaded at: https: //www.mdpi.com/article/10.3390/agronomy12071628/s1. Table S1: Information about tested rice varieties used in the experiment. Table S2: Yield components of tested hybrid rice varieties in 2019 and 2020. Table S3: Dry matter accumulation and harvest index of tested hybrid rice varieties in 2019 and 2020. Table S4: Number of differentially expressed genes between LLYHZ, JLYHZ, and JLY1212 in spikelets at flowering stage in 2020. Table S5: The overlapping genes differentially expressed between LLYHZ and JLYHZ and between LLYHZ and JLY1212 in panicle at flowering stage. Figure  S1: Climate conditions during rice growing period in 2019 and 2020. Figure  Author Contributions: S.X., W.Y. and K.W. were responsible for the experimental work, and took joint responsibility with F.W. for writing the manuscript and drawing the figures. Y.F. and T.D. were responsible for determining the grain quality. Y.Y. and F.W. was responsible for the experimental design, analyzing the data, writing the manuscript and drawing the figures with W.Y. All authors have read and agreed to the published version of the manuscript.