Abstract: 【Objective】Hardness, a structural feature of seed physical dormancy, is an important trait in soybean domestication. Although hardness is beneficial for seeds to survive in unfavorable environments, it will seriously reduce the emergence rate of soybean in the field, and detrimental to yield and processing quality. Analyzing the QTL and candidate genes using bulked segregant analysis sequencing (BSA-Seq), can provide a theoretical reference for understanding the molecular mechanism of hard seededness in soybean. 【Method】The hard seed mutant M_zp661 was obtained from the seeds of Zhongpin 661 induced by ethyl methane sulfonate (EMS), and was crossed with cultivated soybean Zhonghuang 13 (male parent) to construct recombinant inbred line (RIL) population. The progeny lines were investigated for seed hardness, water absorption capacity and anatomical structure of seed coats. Two types of extreme lines in the RIL population, with hard seeds or with imbibed seeds, were selected to construct DNA mixed pools respectively, and then BSA-Seq technology was used to detect genotype differences in extreme-mixed pools and parents. Euclidean distance (ED), delta SNP-index, and delta InDel-index methods were applied to associate hard seed genetic loci of soybean. Combining with bioinformatics analysis, transcriptome data of different soybean tissues and gene annotation information, candidate genes within significant association regions were predicted. 【Result】In the progenies of M_zp661, all areas of imbibitive seeds had the penetration ability, and the seed volume increased continuously with the soaking time. However, no changes were observed for hard seeds over 36 hours. With the prolonged of soaking time, the seed coat of hard seeds began to shrink locally and gradually spread to other parts, and finally cotyledons recovered their imbibition ability. The hard seed not only has smooth and compact seed coat, but also has regular network structure of cuticle and thicker palisade layer, while numbers of stomata and loose structures, tiny cracks and thinner palisade layer were existed in the imbibed seeds. These results suggest that the seed hardness of M_zp661 may be caused by the impermeability of the seed coat. ED, delta SNP-index and delta InDel-index association analysis methods not only identified the reported seed physical dormancy locus qHS1, but also simultaneously detected the candidate region Chr06: 45897227-47746047, which contains a total of 189 genes. Further, transcriptome data and gene annotation predicted that Glyma.06G275300, which is specifically and highly expressed in seeds, might be the candidate gene for this associated region to regulate soybean seed hardness. 【Conclusion】Seed hardness of soybean mutant M_zp661 was caused by the impermeability of the seed coat, and Glyma.06G275300 was predicted as a candidate gene affecting the structure of seed coat using BSA-Seq.

Key words: soybean, seed hardness, seed coat structure, BSA-Seq, candidate gene