ABSTRACT: BACKGROUND:Grain appearance quality is a main determinant of market value in rice and one of the highly important traits requiring improvement in breeding programs. The genetic basis of grain shape and endosperm chalkiness have been given significant attention because of their importance in affecting grain quality. Meanwhile, the introduction of NGS (Next Generation Sequencing) has a significant part to play in the area of genomics, and offers the possibility for high-resolution genetic map construction, population genetics analysis and systematic expression profile study. RESULTS:A RIL population derived from an inter-subspecific cross between indica rice PYZX and japonica rice P02428 was generated, based on the significant variations for the grain morphology and cytological structure between these two parents. Using the Genotyping-By-Sequencing (GBS) approach, 2711 recombination bin markers with an average physical length of 137.68 kb were obtained, and a high-density genetic map was constructed. Global genetic mapping of QTLs affecting grain shape and chalkiness traits was performed across four environments and the newly identified stable loci were obtained. Twelve important QTL clusters were detected, four of which were coincident with the genomic regions of cloned genes or fine mapped QTL reported. Eight novel QTL clusters (including six for grain shape, one for chalkiness, and one for both grain shape and chalkiness) were firstly obtained and highlighted the value and reliability of the QTL analysis. The important QTL cluster on chromosome 5 affects multiple traits including circularity (CS), grain width (GW), area size of grain (AS), percentage of grains with chalkiness (PGWC) and degree of endosperm chalkiness (DEC), indicating some potentially pleiotropic effects. The transcriptome analysis demonstrated an available gene expression profile responsible for the development of chalkiness, and several DEGs (differentially expressed genes) were co-located nearby the three chalkiness-related QTL regions on chromosomes 5, 7, and 8. Candidate genes were extrapolated, which were suitable for functional validation and breeding utilization. CONCLUSION:QTLs affecting grain shape (grain width, grain length, length-width ratio, circularity, area size of grain, and perimeter length of grain) and chalkiness traits (percentage of grains with chalkiness and degree of endosperm chalkiness) were mapped with the high-density GBS-SNP based markers. The important differentially expressed genes (DEGs) were co-located in the QTL cluster regions on chromosomes 5, 7 and 8 affecting PGWC and DEC parameters. Our research provides a crucial insight into the genetic architecture of rice grain shape and chalkiness, and acquired potential candidate loci for molecular cloning and grain quality improvement.