Rice (Oryza sativa L.) Grain Size, Shape, and Weight-Related QTLs Identified Using GWAS with Multiple GAPIT Models and High-Density SNP Chip DNA Markers

Abstract

This study investigated novel quantitative traits loci (QTLs) associated with the control of grain shape and size as well as grain weight in rice. We employed a joint-strategy multiple GAPIT (Genome Association and Prediction Integrated Tool) models [(Bayesian-information and Linkage-disequilibrium Iteratively Nested Keyway (BLINK)), Fixed and random model Circulating Probability Uniform (FarmCPU), Settlement of MLM Under Progressive Exclusive Relationship (SUPER), and General Linear Model (GLM)]–High-Density SNP Chip DNA Markers (60,461) to conduct a Genome-Wide Association Study (GWAS). GWAS was performed using genotype and grain-related phenotypes of 143 recombinant inbred lines (RILs). Data show that parental lines (Ilpum and Tung Tin Wan Hein 1, TTWH1, Oryza sativa L., ssp. japonica and indica, respectively) exhibited divergent phenotypes for all analyzed grain traits), which was reflected in their derived population. GWAS results revealed the association between seven SNP Chip makers and QTLs for grain length, co-detected by all GAPIT models on chromosomes (Chr) 1–3, 5, 7, and 11, were qGL1-1BFSG (AX-95918134, Chr1: 3,820,526 bp) explains 65.2–72.5% of the phenotypic variance explained (PVE). In addition, qGW1-1BFSG (AX-273945773, Chr1: 5,623,288 bp) for grain width explains 15.5–18.9% of PVE. Furthermore, BLINK or FarmCPU identified three QTLs for grain thickness independently, and explain 74.9% (qGT1Blink, AX-279261704, Chr1: 18,023,142 bp) and 54.9% (qGT2-1Farm, AX-154787777, Chr2: 2,118,477 bp) of the observed PVE. For the grain length-to-width ratio (LWR), the qLWR2BFSG (AX-274833045, Chr2: 10,000,097 bp) explains nearly 15.2–32% of the observed PVE. Likewise, the major QTL for thousand-grain weight (TGW) was detected on Chr6 (qTGW6BFSG, AX-115737727, 28,484,619 bp) and explains 32.8–54% of PVE. The qTGW6BFSG QTL coincides with qGW6-1Blink for grain width and explained 32.8–54% of PVE. Putative candidate genes pooled from major QTLs for each grain trait have interesting annotated functions that require functional studies to elucidate their function in the control of grain size, shape, or weight in rice. Genome selection analysis proposed makers useful for downstream marker-assisted selection based on genetic merit of RILs.