In arid and semi-arid regions, the water and nutrients are scarce in the soil. The phosphorus use efficiency between different wheat genotypes varies greatly. Therefore, identification of low phosphorus-tolerant germplasm and mapping of related loci is helpful for genetic improvement of wheat.

Using 282 Shanxi wheat varieties as materials, twelve seedling morphological indicators were investigated under three phosphorus concentrations, including SDW, RDW, DW, SFW, RFW, FW, MRL, TRL, RS, RV, RD, and RN. Principal component analysis, membership function analysis, and cluster analysis were used to comprehensively evaluate the low phosphorus tolerance characteristics of different varieties at the seedling stage. On this basis, the trait evolution trend and biomass allocation at seedling stage were analyzed. At the same time, GWAS was used to identify significant loci related to the low phosphorus-related traits.

The response of different traits to low phosphorus at the seedling stage was different. Lower phosphorus concentrations led to changes in biomass allocation strategy, and shoot growth was less affected by change in phosphorus concentrations than root growth. The decrease in phosphorus concentration inhibited the growth of shoot, and SDW and SFW were significantly reduced. In contrast, low phosphorus promoted root growth, and the indicators of RDW, RFW, MRL, TRL, RV and RN increased significantly. According to the correlation analysis between D-value and morphological indicators, it was found that MRL and RD could be used as selection indicators for low phosphorus tolerance at seedling stage. Based on D-value clustering analysis, 9 low phosphorus tolerant varieties were selected, including Jinmai 46, Jinmai 61, Youmangdahongjing, Hongtumai, Hongheshang, Baikehong, Baixianmai, Huoshaotou, Baishanmai. Analysing trends in trait evolution showed that cultivars were not directly selected for their ability to tolerate low phosphorus. The ability to tolerate low phosphorus decreased first and then increased over time. Before 2010, there was a decreasing trend in the ability of varieties to tolerate low phosphorus, and after 2010, there was an increase in the ability of varieties to tolerate low phosphorus. GWAS stably detected eight loci with R²>10% in three environments, in which 1A_545074550, 2B_489279799, 6A_166899658 and 6A_273060644 were not reported previously.

The MRL and RD can be used as selection indicators for low phosphorus tolerance at seedling stage. A total of nine varieties were selected through comprehensive evaluation of ability in Shanxi wheat to tolerate low phosphorus during seedling stage. Association analysis detected four novel loci associated with low phosphorus tolerance on chromosomes 1A, 2B and 6A, and the results provide germplasm resources and QTL for future low phosphorus tolerance wheat breeding.