Phosphorus (P) deficiency is an important challenge the world faces while having to increase crop yields. It is therefore necessary to select maize (Zea may L.) genotypes with high phosphorus use efficiency (PUE). Here, we extensively analyzed the biomass, grain yield, and PUE-related traits of 359 maize inbred lines grown under both low-P and normal-P conditions. A significant decrease in grain yield per plant and biomass, an increase in PUE under low-P condition, as well as significant correlations between the two treatments were observed. In a genome-wide association study, 49, 53, and 48 candidate genes were identified for eleven traits under low-P, normal-P conditions, and in low-P tolerance index (phenotype under low-P divided by phenotype under normal-P condition) datasets, respectively. Several gene ontology pathways were enriched for the genes identified under low-P condition. In addition, seven key genes related to phosphate transporter or stress response were molecularly characterized. Further analyses uncovered the favorable haplotype for several core genes, which is less prevalent in modern lines but often enriched in a specific subpopulation. Collectively, our research provides progress in the genetic dissection and molecular characterization of PUE in maize.

Genotype-by-environment interaction (G-by-E) is a common but potentially problematic phenomenon in plant breeding. In this study, we investigated the genotypic performance and two measures of plasticity on a phenotypic and genetic level by assessing 234 maize doubled haploid lines from six populations for 15 traits in seven macro-environments with a focus on varying soil phosphorus levels. It was found intergenic regions contributed the most to the variation of phenotypic linear plasticity. For 15 traits, 124 and 31 quantitative trait loci (QTL) were identified for genotypic performance and phenotypic plasticity, respectively. Further, some genes associated with phosphorus use efficiency, such as Zm00001eb117170, Zm00001eb258520, and Zm00001eb265410, encode small ubiquitin-like modifier E3 ligase were identified. By significantly testing the main effect and G-by-E effect, 38 main QTL and 17 interaction QTL were identified, respectively, in which MQTL38 contained the gene Zm00001eb374120, and its effect was related to phosphorus concentration in the soil, the lower the concentration, the greater the effect. Differences in the size and sign of the QTL effect in multiple environments could account for G-by-E. At last, the superiority of G-by-E in genomic selection was observed. In summary, our findings will provide theoretical guidance for breeding P-efficient and broadly adaptable varieties.