The root tissues play important roles in water and nutrient acquisition, environmental adaptation, and plant development. In this study, a diversity panel of 388 wheat accessions was collected to investigate nine root system architecture (RSA) traits at the three-leaf stage under two growing environments: outdoor pot culture (OPC) and indoor pot culture (IPC). Phenotypic analysis revealed that root development was faster under OPC than that under IPC and a significant correlation was observed between the nine RSA traits. The 660K single-nucleotide polymorphism (SNP) chip was used for a genome-wide association study (GWAS). Significant SNPs with a threshold of −log₁₀ (p-value) ≥ 4 were considered. Thus, 36 quantitative trait loci (QTLs), including 13 QTL clusters that were associated with more than one trait, were detected, and 31 QTLs were first identified. The QTL clusters on chromosomes 3D and 5B were associated with four and five RSA traits, respectively. Two candidate genes, TraesCS2A01G516200 and TraesCS7B01G036900, were found to be associated with more than one RSA trait using haplotype analysis, and preferentially expressed in the root tissues. These favourable alleles for RSA traits identified in this study may be useful to optimise the root system in wheat.

The NAC genes, a large plant-specific family of transcription factors, regulate a wide range of pathways involved in development and response to biotic and abiotic stress. In this study, the NAC transcription factors were identified in 27 green plants, and the results showed that NAC transcription factors in plants undergo an appearance stage from water to land and a number expansion stage from gymnosperm to angiosperm. Investigating the evolutionary process of the NAC transcription factors from diploid species to hexaploid wheat revealed that tandem replications during the polyploidization process is an important event for increasing the number of NAC transcription factors in wheat. Then, the molecular characteristics, phylogenetic relationships, and expression patterns of 462 NAC transcription factors of hexaploid wheat (TaNACs) were analyzed. The protein structure results showed that TaNAC was relatively conservative at the N-terminal that contains five subdomains. All these TaNACs were divided into Group I and Group II by phylogenetic analysis, and the TaNACs in Group I should undergo strong artificial selection based on single nucleotide polymorphism (SNP) analysis. Through genome synteny and phylogenetic analysis, these TaNACs were classified into 88 groups and 9 clusters. The biased expression results of these TaNACs showed that there are 24 groups and 67 groups of neofunctionalization genes under biotic and abiotic stress, respectively, and 16 groups and 59 groups of subfunctionalization genes. This shows that neofunctionalization plays an important role in coping with different stresses. Our study provides new insights into the evolution of NAC transcription factors in hexaploid wheat.