Vis/NIR hyperspectral imaging distinguishes sub-population, production environment, and physicochemical grain properties in rice AbstractRice grain quality is a multifaceted quantitative trait that impacts crop value and is influenced by multiple genetic and environmental factors. Chemical, physical, and visual analyses are the standard methods for measuring grain quality. In this study, we evaluated high-throughput hyperspectral imaging for quantification of rice grain quality and classification of grain samples by genetic sub-population and production environment. Whole grain rice samples from the USDA mini-core collection grown in multiple locations were evaluated using hyperspectral imaging and compared with results from standard phenotyping. Loci associated with hyperspectral values were mapped in the mini-core with 3.2 million SNPs in a genome-wide association study (GWAS). Our results show that visible and near infra-red (Vis/NIR) spectroscopy can classify rice according to sub-population and production environment based on differences in physicochemical grain properties. The 702鈥?00鈥塶m range of the NIR spectrum was associated with the chalky grain trait. GWAS revealed that grain chalk and hyperspectral variation share genomic regions containing several plausible candidate genes for grain chalkiness. Hyperspectral quantification of grain chalk was validated using a segregating bi-parental mapping population. These results indicate that Vis/NIR can be used for non-destructive high throughput phenotyping of grain chalk and potentially other grain quality properties. IntroductionRice grain quality influences crop value and is important to growers, millers, and processors as well as consumers1. Grain quality in rice is determined by multiple factors including starch composition, cooking quality, and grain size, shape, and translucency (chalky appearance)2. High grain chalk causes grain breakage during milling and loss of crop value impacting domestic and export markets3. Molecular markers are sought as tools for marker-assisted selection (MAS) in rice breeding for traits like grain quality that are complex, difficult to phenotype and are influenced by the production environment4.Rice has a well-defined population structure with two distinct sub-species that are further divided into sub-populations5. The indica (IND) and aus (AUS) sub-populations belong to the sub-species Indica (AUS-IND) and the tropical japonica (TRJ), temperate japonica (TEJ), and aromatic (ARO) sub-populations belong to the Japonica sub-species (TEJ-TRJ). Even though there is wide variation within each of these categories, rice sub-species and sub-populations have phenotypic and genetic differences that influence their adaptation to different environments and are associated with physicochemical traits6. Having knowledge of sub-species and sub-populations in rice is important for breeding programs targeting different production environments and grain market (i.e. short, medium, and long grain markets) classes.Genome-wide association mapping studies (GWAS) have been used in rice to map a wide range of traits7,8. Several rice diversity panels exist that are genotyped at a high density and are suitable for GWAS such as the 3000 rice genomes9, the High Density Rice Array (HDRA)10, and the USDA rice mini-core collection8,11,12. The USDA rice mini-core germplasm collection of 217 accessions is selected to be phenotypically and genotypically representative of the USDA worldwide rice collection11,13, and includes the five sub-populations of O. sativa (AUS, IND, TRJ, TEJ, and ARO). It has been re-sequenced to an average depth of 1.5脳8 and has a filtered genomic dataset of 3.2 million single nucleotide polymorphic (SNP) markers8,12. Therefore, it is an excellent genetic resource to identify chromosomal regions associated with various phenotypic traits. However, analysis of diversity panels can have confounding factors that may mask or produce false associations with the phenotype of interest. Therefore, quantitative trait loci (QTL) mapping using bi-parental recombinant inbred line populations is a means to validate findings from GWAS studies.One of the bottlenecks in mapping of genes for grain quality traits is the intensive labor, time, and expense required to phenotype the diversity of physicochemical traits impacting rice quality. Near-infrared (NIR) spectroscopy has been widely used to determine protein, amylose, oil, and moisture contents in rice whole grain and flour14,15,16 as well as other quality traits (i.e. sensory and starch pasting properties)17 and in other crops18,19. However, Vis/NIR spectroscopy is a rapid analytical tool that assesses samples by utilizing visible and near-infrared regions of the spectrum and has been demonstrated to be useful in detecting biotic and abiotic stress factors in plant products20,21.The aims of this study were to (1) determine if Vis/NIR hyperspectral imaging of whole grain rice can accurately classify samples according to sub-population or production environment, (2) to determine if variation in Vis/NIR wavelengths correlates with grain quality traits, including chalkiness, and (3) to identify specific genomic regions that are associated with the variation measured through Vis/NIR that can be used to identify corresponding candidate genes.Materials and MethodsGrain productionGrain samples used for hyperspectral imaging came from the USDA rice mini-core collection (221 accessions) including 38 AUS, 86 IND, 33 TEJ, 40 TRJ, 6 ARO, and 18 admixed accessions. The samples were grown in three environments, the USDA-ARS Rice Research Unit/Texas A M Agrilife Research Center located in Beaumont, Texas in 2008 (TX08) and in 2009 (AR09) and 2010 (AR10) at the USDA-ARS Dale Bumpers National Rice Research Center located in Stuttgart, Arkansas. The cultural management practices for the TX08 study (League clay soil, fine, smectitic, hypothermic Oxyaquic Dystrudert) have been described by Pinson et al.22. Li et al.23 described the cultural practices used for the AR09 study which were the same for AR10 (unpublished) (Dewitt silt loam soil, fine, smectitic, thermic, Typic Albaqualf). Field plots were drill seeded to a depth of approximately 2鈥塩m in hill plots in TX08 and in row plots in AR09 and AR10. At both locations, plots were irrigated after seeding to enhance uniform germination. When seedlings reached a height of approximately 9鈥塩m, a permanent flood was applied for the remainder of the season. At the TX08 location, the total fertilizer applied was 73鈥塳g鈥塰a鈭? N as urea and 33.6鈥塳g鈥塰a鈭? P; whereas the fertilizers applied in AR09 were 55鈥塳g鈥塰a鈭? of N, 34鈥塒鈥塳g ha鈭?, 67鈥塊鈥塳g ha鈭? and 11鈥塟n聽kg鈥塰a鈭? ; and in AR10 were 55鈥塳g鈥塰a鈭? of N, 56鈥塒鈥塳g ha鈭?, and 67鈥塊鈥塳g ha鈭? according to soil analyses and local recommendations. As plots reached maturity they were harvested by hand, threshed and dried to approximately 12% moisture prior to storage at 4鈥壜癈 and 50% humidity as rough rice. Prior to imaging, grain samples stored at 4鈥壜癈 were transferred to a desiccator for 2鈥? days to prevent accumulation of聽moisture on the grain as it equilibrated to room temperature (10鈥壜癈) prior to imaging. Arkansas weather data in Stuttgart was provided from the United States Department of Agriculture (USDA) Agricultural Research Service Dale Bumpers National Rice Research Center (https://www.ars.usda.gov/southeast-area/stuttgart-ar/dale-bumpers-national-rice-research-center/docs/weather-data-archives/), and Texas weather information in Beaumont was from the National Weather Service (NWS) Cooperative Observer Program (COOP) (http://www.ncdc.noaa.gov/data-access/land-based-station-data/land-based-datasets/cooperative-observer-network-coop).In addition, samples came from a field study using the bi-parental long grain mapping population KBNT lpa1鈥? x Zhe733 recombinant inbred lines (KZ-RILs) that was conducted at the USDA-ARS Dale Bumpers National Rice Research Center/University of Arkansas Rice Research and Extension Center, Stuttgart, Arkansas24. The population of 187 KZ-RILs was planted on two planting dates, about 30 days apart, during 2013 and 2014 using a randomized complete block design with three replications. Details of the study are described by Edwards et al. (2017) and were similar to the AR09 and AR10 studies in terms of crop management, harvest methods, and grain analysis. For the current study, only samples from the second planting date in 2013 and the first planting date in 2014, were used. These two environments were determined to be the most diverse in terms of climatic differences and one field replication was used from each year. All rough rice samples were dehulled (Satake Rice Machine, Satake Engineering Co., Ltd, Tokyo, Japan) to produce at least 100 grains of brown rice for hyperspectral imaging. The seeds used for hyperspectral imaging came from the same seed source as those used for percent chalk but were different sub-samples.Hyperspectral imaging system and image acquisitionThe reflectance values of brown rice samples were obtained using a visible and near infra-red (Vis/NIR) hyperspectral imaging system that was developed at the Environmental Microbiological and Food Safety Laboratory, Agricultural Research Service, USDA in Beltsville, MD, USA25. It consists of an electron-multiplying charge-coupled-device (EMCCD) camera (Luca, Andor Technology Inc. CT, USA), a spectrograph, six 100-watt halogen bulbs at a distance of 50鈥塩m to the rice samples with 15掳 angle, and a translational stage25,26. The system acquires spectral wavelengths in the range of 400鈥塶m to 1004鈥塶m. Grain samples were dried under vacuum desiccator for 2鈥? days prior to imaging to avoid any inconsistent moisture content across the samples. Ten samples (about 100 grains per sample) were taken per image.聽The exposure time for each sample was set to 5 milliseconds. The final size of the hyperspectral image was 502 脳 600 pixels. The original hyperspectral image was calibrated for white and dark balance using the following equation:$${rm{I}}=,frac{{{rm{I}}}_{0}-{rm{D}}}{{rm{W}}-{rm{D}}}$$ where I is the relative reflectance of the hyperspectral image (scaled from no reflectance at a value of 0 to 100% reflectance at a value of 1), I0 is the original image data, W is the white-reflectance image data, and D is the dark-current image data. The region of interest (ROI) for the hyperspectral images was manually extracted at 750鈥塶m, as this wavelength had the highest intensity with intact brown rice kernel. The binary image without background noise from non-grain areas of the image was used for obtaining the spectral information corresponding to each pixel in the hyperspectral data. Eight mathematical pre-processing methods; smoothing, mean normalization, maximum normalization, range normalization, multiplicative scatter correction, standard normal variate (SNV), Savitzky-Golay first derivative and Savitzky-Golay second derivative; were applied to the raw data before developing the optimal model27,28,29,30 Image correction, segmentation, and spectral-data extraction were customized using the MATLAB software31.Grain quality traitsData on the physicochemical traits for the USDA rice mini-core collection used in this study were previously described in Huggins et al. 2019. This included apparent amylose content (AAC), ASV (an indicator of gelatinization temperature), grain length, width and thickness (mm) (N=197). Brown (unpolished) rice was used to determine percent chalk, and only non-pigmented bran (i.e. white, light brown, and brown) accessions of the mini-core (N= 137) were used because purple and red pericarp can mask the chalk phenotype. For the mini-core, an image analysis system as described in Edwards et al. (2017) was used to determine grain length and grain width on brown rice, whereas hand calipers were used to determine grain thickness on 20 kernels. Milled rice was used to determine AAC and ASV values. Grain quality trait data for the bi-parental mapping population KBNT lpa1鈥? x Zhe733 have been previously described by Edwards et al. (2017) and include percent chalk measured using brown rice.Scanning electron microscopy analysis of grain chalkScanning electron microscopy (SEM) was used to validate the percent chalk as determined by the image analysis phenotyping method (described above) and determine if chalk differences in rice grains are associated with spectral differences. Brown rice samples from KZ-RILs (N=4) with diverse chalk phenotypes were analyzed as these did not differ dramatically for grain dimension or bran color like the mini-core accessions. Whole grains were cut longitudinally with a razor blade and placed on 15鈥壝椻€?0鈥塵m copper plates using ultra smooth, round (12鈥塵m diameter) carbon adhesive tabs (Electron Microscopy Sciences, Inc., Hatfield, PA, USA). The samples were transferred to the Quorum PP2000 cryo-prep-chamber (Quorum Technologies, East Sussex, UK) attached to an S-4700 field emission scanning electron microscope (Hitachi High Technologies America, Inc., Dallas, TX, USA). The specimens were coated with a 10鈥塶m layer of platinum using a magnetron sputter head equipped with a platinum target in the Quorum PP2000. After coating, the specimens were transferred to the SEM for observation. An accelerating voltage of 5鈥塳V with a working distance of 10鈥塵m was used to view specimens. Images were captured using a 4pi Analysis System (Durham, NC, USA).Genome-wide association analysisWhole-genome SNP data for the USDA mini-core diversity panel was obtained from resequencing by Wang et al., (2016). The raw reads generated by Wang et al., (2016) for 203 mini-core accessions were downloaded from the sequence read archive and called against the Michigan State University version 7 (MSU7) rice pseudomolecules32 using the Genome Analysis Toolkit (GATK)33 as described in Huggins et al., (2019). A set of approximately 3.2 million SNPs were generated after filtering loci with a missing rate over 20% and minor allele frequency (MAF) less than 0.05. Principal components (PC) and a kinship matrix were generated using the 3.2 million SNPs in TASSEL version 534. The first three PCs and the kinship matrix were incorporated into a mixed linear model (MLM) to account for relatedness and population stratification. The GWAS analysis of the first principal component for spectral values in the range 702鈥?22鈥塶m spectral range was performed with an MLM that included the options 鈥渆ach marker鈥?and 鈥渘o-compression鈥?to determine trait marker associations in TASSEL version 534,35. Only the non-pigmented bran subset (N=132) was used for GWAS. The GWAS significance threshold for p-value was determined through false discovery rate (FDR) analysis using the R package 鈥榪value鈥?sup>36,37. Based on FDR analysis, the p-value threshold for genome-wide significance was set at 10鈭? and only SNPs that met or exceeded this value were considered as significant. The p-values of the GWAS was visualized as Manhattan and Q-Q plots using the R package 鈥榪qman鈥?sup>38 (Fig.聽6; Supplementary Fig.聽S4).Figure 1Regression beta coefficient and variable importance plot (VIP) of PLS-DA models that discriminate TRJ (cyan), TEJ (blue), and IND (orange) using centered and scaled Vis/NIR spectra with 24 selected predictive wavelengths indicated (A). Population classification of TRJ, TEJ, and IND using PLS-DA models and percent variance explained (PVE) for each population (B). Canonical plot for linear discriminant analysis of TRJ, TEJ, and IND using the 24 predictive wavelengths selected from the PLS-DA model (C).Full size imageCandidate genomic regionsThe GWAS output was used to determine chromosomal regions that are significantly associated with grain quality traits. Chromosomal regions (loci) containing putative QTLs were defined as spanning from 100鈥塳b on either side of a significant SNP and extended if additional significant SNPs were contained in the 100鈥塳b and calculated with a Perl script as described in Huggins et al. (2019). The chromosome, start and stop positions, the most significant SNP (peak SNP) and the p-value of the region were output to a text file for each hyperspectral trait. Text files of identified chromosomal segments for hyperspectral traits and grain quality traits (previously published in Huggins et al., 2019) were used as input for a Perl script to identify overlapping chromosomal regions. The script compared identified significant segments for hyperspectral traits with grain quality traits from three environments (TX08, AR09, and AR10). The identified overlapping chromosomal segments between hyperspectral traits and grain quality traits were processed with another Perl script to detect candidate genes. The above output text files, as well as the MSU7 gene annotation32 of the Oryza sativa genome was added as input for a Perl script to detect genes within a 150鈥塳b distance to either side of the peak SNP in each region. The chromosome number, peak SNP position, gene name and distance from peak SNP were output to a text file. Additionally, Ricebase (https://ricebase.org)39, Oryzabase (https://shigen.nig.ac.jp/rice/oryzabase/), and SNPSeek40 resources were used to inspect significant segment regions and identify candidate genes.Statistical analysisIn this study, Vis/NIR images were analyzed to detect differences in (i) population structure, (ii) environmental differences, and (iii) quality traits. Adjusted means across years and locations of the original spectra were used in Partial Least Squares Discriminant Analysis (PLS-DA) to predict the quality variable Y matrix (group) using the process variable X matrix (Vis/NIR imaging)41. For the mini-core population structure聽study (i), the Vis/NIR hyperspectral image data of five sub-populations (P5), i.e. 鈥楢RO鈥? 鈥楢US鈥? 鈥業ND鈥? 鈥楾EJ鈥? and 鈥楾RJ鈥? and 2 sub-species (P2), 鈥楢US-IND (INDICA)鈥?and 鈥楾EJ-TRJ (JAPONICA)鈥?were collected and the dependent values of each category were coded as binary variables corresponding to membership (coded as 1) or non-membership (coded as zero) in the group. The ARO and AUS groups were later excluded from statistical analysis due to small sample size (4 for ARO, and 23 for AUS). Selected highly informative variables from the PLS-DA were used in a linear discriminant analysis to develop and validate a model for predicting sub-population. These variables were selected based on Variable Importance Plot (VIP) scores and coefficient values. The Vis/NIR image data for the environmental comparisons (TX08, AR09, and AR10) were also coded as binary variables corresponding to membership (coded as 1) or non-membership (coded as zero) in the group and analyzed. For both the population structure and environmental difference PLS-DA, a holdback set of 10% of the randomized samples, was used to validate the PLS-DA model developed from the training set. The training set consisted of all samples not included in the 10% holdback (validation) set. To identify specific spectral regions associated with grain traits such as AAC, ASV, bran color, grain length, width and thickness (mm) (N=197), and % grain chalk (N=137), regression analysis was performed using the imaging and grain trait data from AR09 mini-core samples. The chalk phenotype was further investigated because it showed the best correlation with imaging data. The selected range of wavelengths based on the PLS-DA model was verified by two-way clustering analysis performed by the MeV program42. The results were presented as a heatmap with clustering distances (Pearson鈥檚 correlation). Prior to the clustering analysis, the data for spectra and percent chalk were normalized as 0 to 1.Results and DiscussionVis/NIR phenotypic differences by sub-populationTo explore whether phenotypic characteristics evaluated through a聽Vis/NIR hyperspectral imaging system can discern population structure, the mean percent reflectance across the spectra among the five sub-populations (abbreviated as P5) and that of the 2 sub-species (abbreviated as P2) were calculated (Supplementary Fig.聽S1A,B). Regression analyses of the pre-processed values were found to be similar to results using raw values since image acquisition was performed in a controlled manner, i.e. having a consistent air temperature (10鈥壜癈) and completing the imaging process within a 2-week period for all of the mini-core sets and for the bi-parental population. Therefore, these steps were omitted in further data processing. In the P5 comparisons, the reflectance of AUS was lower than the other four sub-populations in the visible (VIS) (400鈥?00鈥塶m) and near-infrared (NIR) (700鈥?000鈥塶m) regions, however, the difference was less pronounced between 700鈥?20鈥塶m in the NIR region. The lower reflectance of AUS is due to a high frequency of red bran accessions in this sub-population.PLS-DA was used to examine population structure differences within the IND, TEJ, and TRJ sub-populations using VIS (400鈥?00鈥塶m) and NIR (鈥?gt;鈥?00鈥塶m) imaging data (Fig.聽1). Because the red and purple grains may reduce reflectance across the spectra, accessions with these bran colors were excluded from the sub-population differentiation analysis. The sample sizes of the ARO group and the AUS group (after excluding red bran) were insufficient, therefore these groups were not tested for hyperspectral differentiation. Regression beta coefficient plots of the PLS-DA model show peaks in both positive and negative directions for classifying IND, TEJ, and TRJ groups, where peak size relates to the relationship of the wavelength predictor variable to the population classification response variable41, and the variable importance plot (VIP) indicates the influence of each wavelength on the predicted (sub-population) outcome (Fig.聽1A). The PLS-DA actual vs predicted group classifications are shown in Fig.聽1B with percent variance explained (PVE) of 29%, 62%, and 73% for TRJ, TEJ, and IND respectively. Based on coefficients and the VIP, 24 wavelengths were selected for differentiating sub-populations, as indicated by vertical lines in Fig.聽1A. Of the selected wavelengths, 13 were in the visible region (400鈥?00鈥塶m) and 11 were in the NIR region (700鈥?000鈥塶m). Within the visible spectrum, four selected wavelength peaks were found in the violet range (380鈥?50), two in the blue range (450鈥?95鈥塶m), two in the green range (495鈥?70鈥塶m), zero in the yellow range (570鈥?90鈥塶m), two in the orange range (590鈥?20鈥塶m), and five in the red color range (620鈥?50鈥塶m). Selected wavelengths in portions of the NIR spectrum are in a range that is聽indicative of N-H stretching and C-H stretching43. All 24 selected wavelengths contributed to differentiating IND from TEJ and TRJ, and wavelengths at 558.6, 601.7, 659.2, 711.9, 846.0, 898.7, 917.9, and 956.2鈥塶m contributed to differentiating TRJ from TEJ. The greater number of wavelengths differentiating IND likely is due to the greater genetic distance of IND to the TEJ and TRJ sub-populations which both belong to the Japonica subspecies. These observations were expected and are consistent with genetic differentiation between rice sub-populations and sub-species based on phenotypes, DNA markers, and sequencing9,44. The large oscillations in the beta coefficients between similar wavelengths, e.g. the cluster of 419.7, 434.1, and 448.5鈥塶m and the cluster of 759.8, 769.4, and 783.8鈥塶m, were unexpected and the reason for the pattern is unclear.Figure 2Regression beta coefficient and variable importance plot (VIP) of PLS-DA models that discriminate TX08 (red), AR09 (blue), and AR10 (green) using centered and scaled Vis/NIR spectra with the selected predictive wavelengths indicated (A). Classification of TX08, AR09, and AR10 using PLS-DA models and percent variance explained (PVE) for each environment with regression line (red) (B). Canonical plot for linear discriminant analysis of TX08 vs. AR09 and AR10 using the 9 predictive wavelengths selected from the PLS-DA model (C), AR09 vs. TX08 and AR10 using the 8 predictive wavelengths selected from the PLS-DA model (D), and AR10 vs. TX08 and AR09 using the 3 predictive wavelengths selected from the PLS-DA model (E).Full size imageThe 24 predictive wavelengths from the PLS-DA were used in linear discriminant analysis to develop a model for predicting sub-population from Vis/NIR data. A canonical plot from the discriminant analysis shows distinct clustering by sub-population (Fig.聽1C). The discriminant analysis model achieved 70.4% prediction accuracy in the training set and 73.3% accuracy in a holdback validation set consisting of 10 randomly selected accessions from each of the three modeled sub-populations at (Supplementary Table聽S1). In both the holdback set, IND had the greatest prediction accuracy at 90%, followed by TRJ at 70% and TEJ at 60% correctly classified. These results indicate that genomic differences due to population structure can be detected by a combination of visible and NIR hyperspectral imaging phenotypes. The mini-core includes accessions with white, light brown, brown, red, and purple colored bran. With purple and red bran accessions excluded, we speculate that the predictive wavelengths are differentiating variation in the range of white, light brown or brown colored bran. Sub-population prediction accuracy may be increased with a larger number of accessions or greater replication per accession. Testing for hyperspectral differentiation of all five rice sub-populations would require greater representation of non-red accessions from the AUS and ARO groups.Vis/NIR phenotypic difference by environmentTo determine whether the growing environment resulted in significant differences in grain physicochemical profiles, PLS-DA was used to examine differences between the TX08, AR09, and AR10 environments using VIS (400鈥?00鈥塶m) and NIR (鈥?gt;鈥?00鈥塶m) imaging data (Fig.聽2A). Because production environment does not change bran color, all grain samples were used to evaluate any environmental differences. The wavelengths differentiating among TX08, AR09 or AR10 groups were selected based on coefficient values of contrasting peaks among groups having a VIP threshold of 0.8. Wavelengths at 405.4, 414.9, 424.5, 438.9, 458.1, 755, 831.7, 870, and聽937鈥塶m contributed to differentiating TX08 from AR09 and AR10, wavelengths at 745.4, 769.4, 783.8, 879.6, 913, 917.9, and 956.2鈥塶m contributed to differentiating AR09 from AR10 and TX08, and wavelengths at 750.2, 850.8, and 893.9鈥塶m contributed to differentiating AR10 from AR09 and TX08. The PLS-DA actual vs predicted group classifications with PVE of 98%, 96%, and 98% for TX08, AR09, and AR10, respectively (Fig.聽2B). Of the selected wavelengths, five were in the visible region (400鈥?00鈥塶m), and 15 were in the NIR region (700鈥?000鈥塶m). Within the visible spectrum all selected wavelength peaks were found in the violet range (380鈥?50), and discriminated TX08 from AR09 and AR10 groups. Wavelengths in the visible spectrum are known to be associated with traits such as response to biotic and abiotic stress20,45.Figure 3Beta coefficients of the PLS-DA model representing % chalk difference (A), and a heatmap of hyperspectral relative reflectance showing absolute intensity values across populations and production environments in the Vis/NIR spectra with individuals on the y-axis. (B). The x-axes for both A and B figures are wavelengths ranging from 400 to 1000鈥塶m.Full size imageThe predictive wavelengths from the PLS-DA were used in linear discriminant analysis to develop a model for predicting environmental response from Vis/NIR data. A canonical plot from the discriminant analysis shows distinct clustering of TX08 from AR09 and AR10 groups using the 9 selected predictive wavelengths (Fig.聽2C), AR09 from TX08 and AR10 groups using the 8 selected predictive wavelengths (Fig.聽2D), and AR10 from TX08 and AR09 groups using the 3 selected predictive wavelengths (Fig.聽2E). The discriminant analysis model was used to determine % prediction accuracy in the training set and % accuracy in a holdback set consisting of 10 randomly selected accessions from each of the three modeled environments at (Supplementary Table聽S2, A to C). In both the training and holdback sets, TX08 vs. AR09 and AR10 comparison had the greatest prediction accuracy at 73.7%, followed by AR10 vs. TX08 and AR09 comparison at 53.9% and AR09 vs. TX08 and AR10 comparison at 32% correctly classified. These results indicate that different environmental responses can be detected by a combination of visible and NIR hyperspectral imaging phenotypes.Soil type and weather are potential major factors that can explain observed environmental differences. In Arkansas, the soil is a Dewitt silt loam soil, whereas in Texas it is a League clay soil, and this could have produced a difference between Arkansas and Texas environments. Coupled with the weather during the rice growing season, May to November (5鈥? months long because the mini-core collection contains both early and late maturing accessions) in Texas is generally warmer and extends the growing season longer than in Arkansas. The weather data for TX in 2008 and Arkansas in 2009 and 2010 indicated that the average air temperatures were higher in TX08 and AR10 compared to AR09 (Analysis of Variance, Prob F, p=0.0002) and the accumulative solar radiation was higher in AR10 than in TX08 and AR09 (Analysis of Variance, Prob F, p鈥?lt;鈥?.0001) during the growing season (Supplementary Fig.聽S3, and Table聽S3). The accumulative rainfall during the growing season was higher in AR09 and TX08 compared to AR10 (Analysis of Variance, Prob F, p=0.012) (Supplementary Fig.聽S3, and Supplementary Table聽S3).Hyperspectral analysis and genotyping of grain chalkPrevious studies identified specific NIR spectral regions detecting starch, protein, and fat content, and weight of grains14,15,16, demonstrating hyperspectral imaging systems as a potential means of detecting grain quality traits. We questioned whether variation in various rice grain quality traits in the mini-core can be detected using the hyperspectral imaging system. The chalk phenotype became the primary focus among grain traits we tested, because: 1) in rice it is a very important grain trait that affects crop value 1,2,聽2) it showed the most significant correlation with hyperspectral data, and 3) a bi-parental mapping population segregating for percent chalk that was previously phenotyped24 allowed us to verify specific spectral regions identified from the聽USDA mini-core collection GWAS.To identify Vis/NIR spectral regions that are associated with grain quality traits, we first performed K-means clustering analysis of hyperspectral imaging data from a single environment (AR09 as the representative environment). The Vis/NIR spectral regions were divided into several groups based on clustering of wavelengths. This resulted in the identification of five distinctive wavelength groups (Groups 1 to 5), which included spectral regions of 395鈥?24, 429鈥?87, 592鈥?02, 702鈥?22, and 927鈥?004鈥塶m, respectively (data not shown).PCA was performed using the hyperspectral image data for Groups 1 to 5, and the resulting first principal component (PC1) (accounting for 95% of the variation) values were regressed with actual grain quality trait values. Significant correlations (p鈥?lt;鈥?.05) were found with percent chalk, percent amylose, ASV, and kernel bran color (white, light brown, and brown classes only), width, length, and thickness traits, with R2 values ranging from 0.09 to 0.89. The highest R2 was found between the聽chalk phenotype and the PC1 for hyperspectral data in the 702鈥?22鈥塶m region (Group 4) (R2 = 0.93, p鈥?lt;鈥?.0001) (Table聽1). To further evaluate the 702鈥?22聽nm region that is closely associated with percent chalk, a PLS-DA model was developed, and it was determined that wavelengths in the range of 690鈥?20鈥塶m were highly correlated with percent chalk (Fig.聽3A). This region was also consistent across production environments, i.e. TX08, AR09, and AR10 (Fig.聽3B) and encompassed the NIR regions that distinguished the three environments (Fig.聽2A).Table 1 ANOVA result of five Vis/NIR spectral groups displaying phenotypic variation for grain traits, % chalk, % amylose, ASV, bran color, width, length, and thickness among white grain samples from the AR09 environment. R2 values are presented between PC1 components per wavelength group and grain quality trait.Full size tableFigure 4Hierarchical clustering heatmap of the hyperspectral images for AR09 samples and their % chalk in Vis/NIR spectra (A) and the selected region (692鈥?94鈥塶m) associated with % chalk clustered by Euclidean distance of individuals on the x-axis and wavelengths on the y-axis with subpopulation assignments of individual lines indicated on the far right聽(B). r.u. stands for relative unit.Full size imageOne-way hierarchical clustering analysis using AR09 samples further verified spectral regions which showed a similar pattern of percent chalk in the mini-core (Fig.聽4A). Three distinct clusters at 448鈥?49聽nm, 1004鈥塶m, 692鈥?94鈥塶m, and 904鈥?99鈥塶m were identified. Among those clusters, the 692鈥?94鈥塶m region was most strongly associated with percent chalk phenotype (Fig.聽4A). Region 692鈥?94鈥塶m was further studied to determine more specifically which wavelengths were associated with percent chalk when categorized as high, intermediate, or low chalk using the Euclidean distance聽metric. Hierarchical cluster results showed that region 692鈥?74鈥塶m better separated low chalk accessions from intermediate and high chalk accessions, while the 779鈥?94鈥塶m region distinctly separated low from intermediate, and intermediate from high chalk accessions (Fig.聽4B). There was no obvious relationship between sub-species and percent chalk categories (Fig.聽4B).Figure 5PLS image of a subset of KZ-RILs with diverse levels of grain chalk (A), linear regression of predicted values of a PLS-DA model using the 702鈥?22鈥塶m range as a function of actual % chalk observed in brown rice of the KZ-RIL population (B), and scanning electron microscope (SEM) images of low (C), intermediate (D), and high chalk KZ-RILs (E). a.u. in (A) stands for arbitrary unit. Black and yellow dots in (B) display calibrated and predicted values, respectively. Grey dotted line/red solid line, grey solid line/blue solid line, and black lines in (B) represent 95% prediction and confidence bands and regression of calibrated/predicted values, respectively.Full size imageDiversity panels can be an excellent genetic resource to identify chromosomal regions associated with various phenotypic traits. However, diverse genetic backgrounds can have confounding factors that may mask or produce false associations with the phenotype of interest. Therefore, we used the bi-parental recombinant inbred line population, KBNT-lpa1 x ZHE733 (KZ-RIL), that is segregating for percent chalk to verify that the 702鈥?22鈥塶m region (from Fig.聽4A) can differentiate chalk phenotypes from a common genetic background. Hyperspectral images of the grain samples of the聽KZ-RIL population grown in two different years were captured (Fig.聽5A, selected extremes), and the 702鈥?22鈥塶m region was used to develop a PLS-DA model. Based on this model, the calibrated correlation coefficient was 0.91 and the predicted correlation coefficient was 0.94 (Fig.聽5B). Among the selected KZ-RILs having divergent percent chalk (Fig.聽5A), SEM images showed that the RIL with high chalk (51.1%) had an overall disorganized and irregular packing of starch granules (Fig.聽5E). Conversely, the RIL with low chalk (5.6%) had a regular and organized cellular structure (Fig.聽5C), and the RIL with intermediate chalk (21.7%) had a cellular structure that was less organized than the RIL with high chalk and more irregular than the RIL with low chalk (Fig.聽5D). Notably, the high chalk RIL had more air聽spaces between starch granules than the low chalk RIL as has been reported previously46.Figure 6Manhattan plots for GWAS of the first principal component for the hyperspectral region 707鈥?22鈥塶m for environments TX08, AR09 and AR10. Only non-pigmented bran accessions are presented in this analysis. The x-axis displays chromosome pseudomolecule coordinates and the y-axis displays the 鈥搇og10(p) value for each SNP across chromosomes. The dark horizontal line represents the genome-wide significance threshold. The labeled SNPs are discussed as potential candidates in the text. GW5 鈥?grain weight 5 gene; LTPL164 鈥?Protease inhibitor/seed storage gene.Full size imagePC1 calculated from the 702鈥?22鈥塶m wavelengths of the hyperspectral image data collected from the mini-core accessions was used for GWAS analysis. Because bran color can mask chalk and affect wavelength absorption, reflectance and transmittance (Figure聽S3), the analysis of only non-pigmented bran accessions was performed to reduce the possibility of confounding effects and false associations. A total of 44 chromosomal candidate regions were identified as associated with spectral PC1 values across two or more environments (Fig.聽6; Supplementary File聽S1). Seventeen segments associated with the first principal component of the hyperspectral values (hyperspectral PC1) were identified for TX08, 49 for AR09, and 50 for AR10. Forty-eight segments were common between AR09 and AR10. Five hyperspectral PC1-associated segments were common between TX08, AR09, and AR10 including two of the segments located on chromosome 4, one on chromosome 5, one on chromosome 7, and the other on chromosome 9 (Supplementary File聽S1). The 68 non-overlapping chromosomal segments from all three environments associated with hyperspectral PC1 were examined for overlap with percent chalk, AAC, ASV, and grain length and width segments detected by GWAS in the mini-core by Huggins et al., (2019) to identify shared chromosomal segments. Grain chalk-associated segments overlapped with other hyperspectral grain trait-associated segments from two or more environments for 21 different loci. The hyperspectral-chalk associated segments also overlapped with 13 AAC-associated segments, one ASV segment, and two grain dimension segments previously reported by Huggins et al., (2019). Some of the significant hyperspectral chalk segments contained characterized known genes or were proximal to them. On chromosome 3, a segment was proximal to the grain size/plant height TIFY11b gene and another co-located with a phosphatidylinositol-4-phosphate 5-kinase gene (PIPK) for AR09 and AR10 (Fig.聽6; Supplementary File聽S1). A segment identified on chromosome 5 for AR09 and AR10 co-located with the characterized chalk5 gene, an inorganic H鈥?鈥塸yrophosphatase. Another segment on chromosome 5, detected in TX08, AR09, and AR10, was proximal to the characterized grain dimension gene, Grain Weight 5 (DQ991205). A segment on chromosome 7, also detected in the three environments, co-located with a seed storage/protease inhibitor gene (LTPL164) (Fig.聽6; Supplementary File聽S1).The shared significant chromosomal segments detected between spectral regions associated with chalk and other grain quality traits were used to identify associated candidate genes. Gene annotations, biological and molecular functions were used to propose candidate rice genes that fell within 150鈥塳b on either side of peak SNP in significant segments. We identified 12 potential candidate genes related to starch and sucrose that met the criteria and overlapped with grain chalk. Several of these candidate genes have vital roles in starch or sucrose biosynthesis. The candidate LOC_Os01g54560, is a trehalose phosphate synthase gene on chromosome 1, ~117鈥塳b upstream of SNP-31,494,903 (the number included in the SNP name represents the physical position in base pairs of the SNP on the MSU7 pseudomolecule assembly) (Table聽2). The synthesis of trehalose can occur via multiple pathways, but the best known involves trehalose-6-phosphatase, which is a known regulator of plant sucrose47,48. Trehalose synthesis and accumulation in plant tissues induce sucrose synthase activity49, thus affecting sucrose and starch biosynthetic activity. Notably, trehalose accumulation in varying plant tissues can be triggered by abiotic stresses such as oxidation, heat, and drought50, where they mitigate the stress effects. The candidate LOC_Os03g07480, is a sucrose transporter (SUT1) gene, ~19鈥塳b upstream of SNP-3,823,250 on chromosome 3 (Table聽2). Sucrose transporters act as distribution centers for photo-assimilates and move sucrose into the phloem of most plants by active transport51,52,53. However, in rice, the SUT1 not only functions in phloem loading but may be involved in sucrose transportation to the sink (grain)53,54,55,56. The SUT1 gene is usually expressed in panicles, leaf sheaths, and leaf blades after heading, and promotes the mobilization of starch to sink tissue57,58,59, thus playing a key role in grain filling and quality. Candidate LOC_Os03g40270, ~49鈥塳b upstream of SNP-22,433,619, is an 伪-1,4-glucan protein synthase (伪-1,4-glucanotransferase) gene located on chromosome 3 (Table聽2). The 伪-1,4-glucan protein synthase gene belongs to a class of enzymes called disproportionating enzymes (DPE)60,61,62. Previous studies have linked this enzyme with starch granule biosynthesis, more specifically, amylopectin structure and a role in starch degradation62,63,64. When this disproportionating enzyme was suppressed, it resulted in reduced amylopectin chains and increased amylose, leading to loosely packed starch granules in the endosperm65. Hence, 伪-1,4-glucan protein synthase contributes to the complex processes regulating grain filling. The candidate LOC_Os10g08022赂 is characterized as a fructose-bisphosphate aldose (FBA) gene, ~31鈥塳b downstream of SNP-4,312,440 on chromosome 10 (Table聽2). FBAs catalyze D-glyceraldehyde-3-phosphate (GAP) and dihydroxyacetone phosphate (DHAP) from D-fructose-1,6-bisphosphate (FBP) and is vital to glycolysis and gluconeogenesis66,67,68. FBAs also play a key role in regulating both starch and sucrose biosynthesis in plants69,70,71. Two forms of FBA exist, cytosolic FBA which is involved in sucrose synthesis71,72, and chloroplastic FBA which is vital in starch biosynthesis69. Moreover, inhibition of cytosolic FBA led to decreased sucrose biosynthesis but increased starch biosynthesis72,73. Notably, FBAs not only regulate starch and sucrose biosynthesis but play pivotal roles in various biological, metabolic and physiological pathways that require sucrose, including response to abiotic stresses71,74,75,76.Table 2 Summary of common candidate loci between hyperspectral results and chalk trait in the mini-core diversity panel.Full size tableThe results of this study demonstrate that hyperspectral imaging offers a high-throughput, efficient method of assessing rice grain traits that otherwise require labor and time-consuming assays. With this method, large amounts of data are captured, and multiple traits can be characterized simultaneously. Hyperspectral imaging technique provides a quantitative assessment of grain quality that is repeatable and not dependent on subjective ratings聽(e.g., ASV). NIR hyperspectral imagery consists of numerous bands with small spectrum gaps (every 4鈥塶m in our study) and can assess grain traits such as fat, starch, protein, moisture, color, and many other physicochemical compounds at once. The involvement of multiple bands in the prediction of chalk suggests that there may be multiple mechanisms leading to chalky grains that each have different impacts on the hyperspectral profile. The spectral regions associated with chalk could be further refined for use in a multispectral apparatus for high-throughput quantification of chalk in rice grains. Such an apparatus may improve the efficiency of selection for grain quality in rice breeding programs.GWAS was used to confirm known genes and to identify novel candidate genes affecting grain quality traits using hyperspectral imaging. The PLS-DA models of hyperspectral data identify spectral ranges that distinguish genetic and production environment differences, and this information may help to resolve the genetics of complex traits such as rice grain quality. References1.Fitzgerald, M. A. Resurreccion, A. P. Maintaining the yield of edible rice in a warming world. Funct. Plant Biol. 36, 1037鈥?045 (2009).Article聽Google Scholar聽 2.Lisle, A. J., Martin, M. Fitzgerald, M. A. Chalky and translucent rice grains differ in starch composition and structure and cooking properties. Cereal Chem. J. 77, 627鈥?32 (2000).CAS聽 Article聽Google Scholar聽 3.Bennet, D. Quality paramount to importers of U.S. rice. Delta FarmPress blog (2013). Available at: https://www.deltafarmpress.com/blog/quality-paramount-importers-us-rice.4.Zhao, X. et al. Climate Change: Implications for the yield of edible rice. PLoS One 8, e66218 (2013).ADS聽 CAS聽 PubMed聽 PubMed Central聽 Article聽Google Scholar聽 5.Sweeney, M. McCouch, S. The complex history of the domestication of rice. Ann. Bot. 100, 951鈥?57 (2007).PubMed聽 PubMed Central聽 Article聽Google Scholar聽 6.Ali, L. M. et al. A rice diversity panel evaluated for genetic and agro-morphological diversity between subpopulations and its geographic distribution. Crop Sci. 51, 2021鈥?035 (2011).Article聽Google Scholar聽 7.Zhao, K. et al. Genome-wide association mapping reveals a rich genetic architecture of complex traits in Oryza sativa. Nat. Commun. 2, 467 (2011).8.Wang, H. et al. The power of inbreeding: NGS-based GWAS of rice reveals convergent evolution during rice domestication. Mol. Plant 9, 975鈥?85 (2016).CAS聽 PubMed聽 Article聽Google Scholar聽 9.Wang, W. et al. Genomic variation in 3,010 diverse accessions of Asian cultivated rice. Nature 557, 43 (2018).ADS聽 CAS聽 PubMed聽 PubMed Central聽 Article聽Google Scholar聽 10.McCouch, S. R. et al. Open access resources for genome-wide association mapping in rice. Nat. Commun. 7, 10532 (2016).ADS聽 CAS聽 PubMed聽 PubMed Central聽 Article聽Google Scholar聽 11.Agrama, H. A. et al. Genetic assessment of a mini-core subset developed from the USDA rice genebank. Crop Sci. 49, 1336鈥?346 (2009).Article聽Google Scholar聽 12.Huggins, T. D. et al. Association analysis of three diverse rice (Oryza sativa L.) germplasm collections for loci regulating grain quality traits. Plant Genome 12, 1 (2019).13.Li, X. et al. Genotypic and phenotypic characterization of genetic differentiation and diversity in the USDA rice mini-core collection. Genetica 138, 1221鈥?230 (2010).CAS聽 PubMed聽 Article聽Google Scholar聽 14.Sohn, M., Barton, F. E., McClung, A. M. Champagne, E. T. Near-infrared spectroscopy for determination of protein and amylose in rice flour through use of derivatives. Cereal Chem. 81, 341鈥?44 (2004).CAS聽 Article聽Google Scholar聽 15.Wang, H. L. et al. Quantitative analysis of fat content in rice by near-infrared spectroscopy technique. Cereal Chem. 83, 402鈥?06 (2006).ADS聽 CAS聽 Article聽Google Scholar聽 16.Wu, J. G. Shi, C. H. Prediction of grain weight, brown rice weight and amylose content in single rice grains using near-infrared reflectance spectroscopy. F. Crop. Res. 87, 13鈥?1 (2004).Article聽Google Scholar聽 17.Bao, J. S., Cai, Y. Z. Corke, H. Prediction of rice starch quality parameters by near-infrared reflectance spectroscopy. J. Food Sci. 66, 936鈥?39 (2001).CAS聽 Article聽Google Scholar聽 18.Osborne, B. G. Near-Infrared Spectroscopy in Food Analysis. Encyclopedia of Analytical Chemistry (2006). https://doi.org/10.1002/9780470027318.a101819.Armstrong, P. R., Maghirang, E. B., Xie, F. Dowell, F. E. Comparison of dispersive and Fourier-transform NIR instruments for measuring grain and flour attributes. Appl. Eng. Agric. 22, 453鈥?57 (2006).Article聽Google Scholar聽 20.Baek, I. et al. Selection of optimal hyperspectral wavebands for detection of discolored, diseased rice seeds. Appl. Sci. 9, 1027 (2019).Article聽Google Scholar聽 21.Gao, J., Li, X., Zhu, F. He, Y. Application of hyperspectral imaging technology to discriminate different geographical origins of Jatropha curcas L. seeds. Comput. Electron. Agric. 99, 186鈥?93 (2013).22.Pinson, S. R. M. et al. Worldwide genetic diversity for mineral element concentrations in rice grain. Crop Sci. 55, 294鈥?11 (2015).CAS聽 Article聽Google Scholar聽 23.Li, X. et al. Unraveling the complex trait of harvest index with association mapping in rice (Oryza sativa L.). PLoS One 7, e29350 (2012).ADS聽 CAS聽 PubMed聽 PubMed Central聽 Article聽Google Scholar聽 24.Edwards, J. D., Jackson, A. K. McClung, A. M. Genetic architecture of grain chalk in rice and interactions with a low phytic acid locus. F. Crop. Res. 205, 116鈥?23 (2017).Article聽Google Scholar聽 25.Lee, H. et al. Detection of Cracks on Tomatoes Using a Hyperspectral Near-Infrared Reflectance Imaging System. Sensors 14, 18837鈥?8850 (2014).PubMed聽 Article聽Google Scholar聽 26.Denvir, D. J. Conroy, E. Electron-multiplying CCD: the new ICCD. in In Low-Light-Level and Real-Time Imaging Systems, Components, and Applications 4796, 164鈥?74 (International Society for Optics and Photonics, 2003).27.Barnes, R. J., Dhanoa, M. S. Lister, S. J. Standard normal variate transformation and de-trending of near-infrared diffuse reflectance spectra. Appl. Spectrosc. 43, 772鈥?77 (1989).ADS聽 CAS聽 Article聽Google Scholar聽 28.Rinnan, 脜., van den Berg, F. Engelsen, S. B. Review of the most common pre-processing techniques for near-infrared spectra. TrAC Trends Anal. Chem. 28, 1201鈥?222 (2009).29.Lohumi, S., Lee, S., Lee, H. Cho, B. K. A review of vibrational spectroscopic techniques for the detection of food authenticity and adulteration. Trends in Food Science and Technology 46, 85鈥?8 (2015).CAS聽 Article聽Google Scholar聽 30.Lee, H. et al. Prediction of crude protein and oil content of soybeans using Raman spectroscopy. Sensors Actuators B Chem. 185, 694鈥?00 (2013).CAS聽 Article聽Google Scholar聽 31.MATLAB 8.0 and Statistics Toolbox 8.1. The MathWorks Inc., Natick, Massachusetts, United States.32.Ouyang, S. et al. The TIGR rice genome annotation resource: improvements and new features. Nucleic Acids Res. 35, D883鈥揇887 (2006).PubMed聽 PubMed Central聽 Article聽Google Scholar聽 33.McKenna, A. et al. The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 20, 1297鈥?303 (2010).CAS聽 PubMed聽 PubMed Central聽 Article聽Google Scholar聽 34.Bradbury, P. J. et al. TASSEL: software for association mapping of complex traits in diverse samples. Bioinformatics 23, 2633鈥?635 (2007).CAS聽 PubMed聽 PubMed Central聽 Article聽Google Scholar聽 35.Zhang, Z. et al. Mixed linear model approach adapted for genome-wide association studies. Nat. Genet. 42, 355 (2010).CAS聽 PubMed聽 PubMed Central聽 Article聽Google Scholar聽 36.R Core Team. R: A language and environment for statistical computing. R Found. Stat. Comput. Vienna, Austria (2018).37.Storey, J. D. Tibshirani, R. Statistical significance for genomewide studies. Proc. Natl. Acad. Sci. 100, 9440鈥?445 (2003).ADS聽 MathSciNet聽 CAS聽 PubMed聽 MATH聽 Article聽Google Scholar聽 38.Turner, S. D. qqman: an R package for visualizing GWAS results using QQ and manhattan plots. BioRxiv 005165 (2014). https://doi.org/10.21105/joss.0073139.Edwards, J. D., Baldo, A. M. Mueller, L. A. Ricebase: a breeding and genetics platform for rice, integrating individual molecular markers, pedigrees and whole-genome-based data. Database 2016, (2016).40.Mansueto, L. et al. Rice SNP-seek database update: new SNPs, indels, and queries. Nucleic Acids Res. 45, D1075鈥揇1081 (2016).PubMed聽 PubMed Central聽 Article聽 CAS聽Google Scholar聽 41.Haaland, D. M. Thomas, E. V. Partial least-squares methods for spectral analyses. 1. Relation to other quantitative calibration methods and the extraction of qualitative information. Anal. Chem. 60, 1193鈥?202 (1988).CAS聽 Article聽Google Scholar聽 42.Saeed, A. I. et al. TM4: a free, open-source system for microarray data management and analysis. Biotechniques 34, 374鈥?78 (2003).CAS聽 PubMed聽 PubMed Central聽 Article聽Google Scholar聽 43.Metrohm. Monograph. NIR Spectroscopy: A guide to near-infrared spectroscopic analysis of industrial manufacturing processes. Metrohm AG, CH-9101 Herisau, Switzerland, 46 pp (2014).44.Garris, A. J., Tai, T. H., Coburn, J., Kresovich, S. McCouch, S. Genetic structure and diversity in Oryza sativa L. Genetics 169, 1631鈥?638 (2005).45.Lowe, A., Harrison, N. French, A. P. Hyperspectral image analysis techniques for the detection and classification of the early onset of plant disease and stress. Plant Methods 13, 80 (2017).PubMed聽 PubMed Central聽 Article聽Google Scholar聽 46.Tashiro, T. Wardlaw, I. F. The effect of high temperature on kernel dimensions and the type and occurrence of kernel damage in rice. Aust. J. Agric. Res. 42, 485鈥?96 (1991).Article聽Google Scholar聽 47.Wingler, A., Fritzius, T., Wiemken, A., Boller, T. Aeschbacher, R. A. Trehalose induces the ADP-glucose pyrophosphorylase gene, ApL3, and starch synthesis in Arabidopsis. Plant Physiol. 124, 105鈥?14 (2000).CAS聽 PubMed聽 PubMed Central聽 Article聽Google Scholar聽 48.Lunn, J. E. et al. Sugar-induced increases in trehalose 6-phosphate are correlated with redox activation of ADPglucose pyrophosphorylase and higher rates of starch synthesis in Arabidopsis thaliana. Biochem. J. 397, 139鈥?48 (2006).49.Muller, J., Boller, T. Wiemken, A. Trehalose affects sucrose synthase and invertase activities in soybean (Glycine max [L.] Merr.) roots. J. Plant Physiol. 153, 255鈥?57 (1998).50.Benaroudj, N., Lee, D. H. Goldberg, A. L. Trehalose accumulation during cellular stress protects cells and cellular proteins from damage by oxygen radicals. J. Biol. Chem. 276, 24261鈥?4267 (2001).CAS聽 PubMed聽 Article聽Google Scholar聽 51.Hackel, A. et al. Sucrose transporter LeSUT1 and LeSUT2 inhibition affects tomato fruit development in different ways. Plant J. 45, 180鈥?92 (2006).CAS聽 PubMed聽 Article聽Google Scholar聽 52.Schmitt, B., Stadler, R. Sauer, N. Immunolocalization of solanaceous SUT1 proteins in companion cells and xylem parenchyma: new perspectives for phloem loading and transport. Plant Physiol. 148, 187鈥?99 (2008).CAS聽 PubMed聽 PubMed Central聽 Article聽Google Scholar聽 53.Braun, D. M. Slewinski, T. L. Genetic control of carbon partitioning in grasses: roles of sucrose transporters and tie-dyed loci in phloem loading. Plant Physiol. 149, 71鈥?1 (2009).CAS聽 PubMed聽 PubMed Central聽 Article聽Google Scholar聽 54.Scofield, G. N., Hirose, T., Aoki, N. Furbank, R. T. Involvement of the sucrose transporter, OsSUT1, in the long-distance pathway for assimilate transport in rice. J. Exp. Bot. 58, 3155鈥?169 (2007).CAS聽 PubMed聽 Article聽Google Scholar聽 55.K眉hn, C. Grof, C. P. L. Sucrose transporters of higher plants. Current Opinion in Plant Biology 13, 287鈥?97 (2010).Article聽 CAS聽Google Scholar聽 56.Ayre, B. G. Membrane-transport systems for sucrose in relation to whole-plant carbon partitioning. Molecular Plant 4, 377鈥?94 (2011).CAS聽 PubMed聽 Article聽Google Scholar聽 57.Hirose, T., Imaizumi, N., Scofield, G. N., Furbank, R. T. Ohsugi, R. cDNA cloning and tissue specific expression of a gene for sucrose transporter from rice (Oryza sativa L.). Plant Cell Physiol. 38, 1389鈥?396 (1997).58.Aoki, N. et al. Three sucrose transporter genes are expressed in the developing grain of hexaploid wheat. Plant Mol. Biol. 50, 453鈥?62 (2002).CAS聽 PubMed聽 Article聽Google Scholar聽 59.Aoki, N., Hirose, T., Scofield, G. N., Whitfeld, P. R. Furbank, R. T. The sucrose transporter gene family in rice. Plant Cell Physiol. 44, 223鈥?32 (2003).CAS聽 PubMed聽 Article聽Google Scholar聽 60.Ball, S. et al. From glycogen to amylopectin: a model for the biogenesis of the plant starch granule. Cell 86, 349鈥?52 (1996).CAS聽 PubMed聽 Article聽Google Scholar聽 61.Takaha, T., Yanase, M., Takata, H., Okada, S. Smith, S. M. Potato D-enzyme catalyzes the cyclization of amylose to produce cycloamylose, a novel cyclic glucan. J. Biol. Chem. 271, 2902鈥?908 (1996).CAS聽 PubMed聽 Article聽Google Scholar聽 62.Colleoni, C. et al. Genetic and biochemical evidence for the involvement of 伪-1,4 glucanotransferases in amylopectin synthesis. Plant Physiol. 120, 993鈥?004 (1999).CAS聽 PubMed聽 PubMed Central聽 Article聽Google Scholar聽 63.Colleoni, C. et al. Biochemical characterization of the Chlamydomonas reinhardtii 伪-1,4 glucanotransferase supports a direct function in amylopectin biosynthesis. Plant Physiol. 120, 1005鈥?014 (1999).64.Critchley, J. H., Zeeman, S. C., Takaha, T., Smith, A. M. Smith, S. M. A critical role for disproportionating enzyme in starch breakdown is revealed by a knock-out mutation in Arabidopsis. Plant J. 26, 89鈥?00 (2001).65.Dong, X. et al. Plastidial disproportionating enzyme participates in starch synthesis in rice endosperm by transferring maltooligosyl groups from amylose and amylopectin to amylopectin. Plant Physiol. 169, 2496鈥?512 (2015).CAS聽 PubMed聽 PubMed Central聽Google Scholar聽 66.Rutter, W. J. Evolution of aldolase. Fed. Proc. Am. Soc. Exp. Biol. 23, 1248鈥?257 (1964).CAS聽Google Scholar聽 67.Berg, I. A. et al. Autotrophic carbon fixation in archaea. Nat. Rev. Microbiol. 8, 447 (2010).CAS聽 PubMed聽 Article聽Google Scholar聽 68.Lv, G.-Y. et al. Molecular characterization, gene evolution, and expression analysis of the fructose-1, 6-bisphosphate aldolase (FBA) gene family in wheat (Triticum aestivum L.). Front. Plant Sci. 8, 1030 (2017).69.Sonnewald, U., Lerchl, J., Zrenner, R. Frommer, W. Manipulation of sink-source relations in transgenic plants. Plant Cell Environ. 17, 649鈥?58 (1994).CAS聽 Article聽Google Scholar聽 70.Anderson, L. E., Bryant, J. A. Carol, A. A. Both chloroplastic and cytosolic phosphoglycerate kinase isozymes are present in the pea leaf nucleus. Protoplasma 223, 103鈥?10 (2004).CAS聽 PubMed聽Google Scholar聽 71.Fan, W., Zhang, Z. Zhang, Y. Cloning and molecular characterization of fructose-1,6-bisphosphate aldolase gene regulated by high-salinity and drought in Sesuvium portulacastrum. Plant Cell Rep. 28, 975鈥?84 (2009).72.Zrenner, R., Krause, K.-P., Apel, P. Sonnewald, U. Reduction of the cytosolic fructose-1,6-bisphosphatase in transgenic potato plants limits photosynthetic sucrose biosynthesis with no impact on plant growth and tuber yield. Plant J. 9, 671鈥?81 (1996).CAS聽 PubMed聽 Article聽Google Scholar聽 73.Strand, A. et al. Decreased expression of two key enzymes in the sucrose biosynthesis pathway, cytosolic fructose-1,6-bisphosphatase and sucrose phosphate synthase, has remarkably different consequences for photosynthetic carbon metabolism in transgenic Arabidopsis thaliana. Plant J. 23, 759鈥?70 (2000).74.Haake, V., Zrenner, R., Sonnewald, U. Stitt, M. A moderate decrease of plastid aldolase activity inhibits photosynthesis, alters the levels of sugars and starch, and inhibits growth of potato plants. Plant J. 14, 147鈥?57 (1998).CAS聽 PubMed聽 Article聽Google Scholar聽 75.Michelis, R. Gepstein, S. Identification and characterization of a heat-induced isoform of aldolase in oat chloroplast. Plant Mol. Biol. 44, 487鈥?98 (2000).CAS聽 PubMed聽 Article聽Google Scholar聽 76.Henkes, S., Sonnewald, U., Badur, R., Flachmann, R. Stitt, M. A small decrease of plastid transketolase activity in antisense tobacco transformants has dramatic effects on photosynthesis and phenylpropanoid metabolism. Plant Cell 13, 535鈥?51 (2001).CAS聽 PubMed聽 PubMed Central聽 Article聽Google Scholar聽 Download referencesAcknowledgementsMention of a trademark or proprietary product does not constitute a guarantee or warranty of the product by the U.S. Department of Agriculture and does not imply its approval to the exclusion of other products that also can be suitable. USDA is an equal opportunity provider and employer. All experiments complied with the current laws of the United States, the country in which they were performed. The authors thank Lorie Bernhardt for providing the mini-core accessions seed for AR10 through the Genetic Stock Oryza (GSOR) (www.ars.usda.gov/GSOR) repository. We also would like to thank Dr. Stephen R. Delwiche and Dr. Brook T. Moyers for many insightful comments and suggestions.Author informationAuthor notesThese authors contributed equally: Jinyoung Y. Barnaby and Trevis D. Huggins.AffiliationsDale Bumpers National Rice Research Center, United States Department of Agriculture - Agricultural Research Service, Stuttgart, AR, 72160, USAJinyoung Y. Barnaby,聽Trevis D. Huggins,聽Anna M. McClung,聽Shannon R. M. Pinson聽 聽Jeremy D. EdwardsEnvironmental Microbial and Food Safety Laboratory, United States Department of Agriculture - Agricultural Research Service, Beltsville, MD, 20705, USAHoonsoo Lee,聽Mirae Oh聽 聽Moon S. KimDepartment of Biosystems Engineering, Chungbuk National University, Cheongju, 28644, Republic of KoreaHoonsoo LeeElectron Confocal Microscopy Unit, United States Department of Agriculture - Agricultural Research Service, Beltsville, MD, 20705, USAGary R. BauchanGrassland and Forages Division, National Institute of Animal Science, Rural Development Administration, Cheonan, 31000, Republic of KoreaMirae OhTexas A M AgriLife Research Center, Texas A M University System, Beaumont, TX, 77713, USALee TarpleyNational Institute of Horticultural and Herbal Sciences, Rural Development Administration, Haman, 52054, Republic of KoreaKangjin LeeAuthorsJinyoung Y. BarnabyView author publicationsYou can also search for this author in PubMed聽Google ScholarTrevis D. HugginsView author publicationsYou can also search for this author in PubMed聽Google ScholarHoonsoo LeeView author publicationsYou can also search for this author in PubMed聽Google ScholarAnna M. McClungView author publicationsYou can also search for this author in PubMed聽Google ScholarShannon R. M. PinsonView author publicationsYou can also search for this author in PubMed聽Google ScholarMirae OhView author publicationsYou can also search for this author in PubMed聽Google ScholarGary R. BauchanView author publicationsYou can also search for this author in PubMed聽Google ScholarLee TarpleyView author publicationsYou can also search for this author in PubMed聽Google ScholarKangjin LeeView author publicationsYou can also search for this author in PubMed聽Google ScholarMoon S. KimView author publicationsYou can also search for this author in PubMed聽Google ScholarJeremy D. EdwardsView author publicationsYou can also search for this author in PubMed聽Google ScholarContributionsJ.B. and A.M. conceived and designed the experiments; J.B. and M.O. performed the imaging experiments, and G.B. conducted the SEM analysis; A.M., S.P., and L.T. conducted the field experiments for grain samples and grain analysis; J.B., H.L., K.L., and M.K. analyzed the imaging data; T.H. and J.E. analyzed the GWAS data; J.B., T.H., S.P., A.M., and J.E. wrote the manuscript; J.B., T.H., H.L., A.M., M.O., G.B., S.P., L.T., K.L., M.K., and J.E. approved manuscript for publication.Corresponding authorCorrespondence to Jeremy D. Edwards.Ethics declarations Competing interests 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. Additional informationPublisher鈥檚 note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.Supplementary information CommentsBy submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate. Sign up for the Nature Briefing newsletter 鈥?what matters in science, free to your inbox daily.