Next Article in Journal
Therapeutic Interventions to Mitigate Mitochondrial Dysfunction and Oxidative Stress–Induced Damage in Patients with Bipolar Disorder
Next Article in Special Issue
Linking Multi-Omics to Wheat Resistance Types to Fusarium Head Blight to Reveal the Underlying Mechanisms
Previous Article in Journal
Autophagy-Related Gene PlATG6a Is Involved in Mycelial Growth, Asexual Reproduction and Tolerance to Salt and Oxidative Stresses in Peronophythora litchii
Previous Article in Special Issue
Genome-Wide Identification and Characterization of RNA/DNA Differences Associated with Drought Response in Wheat
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genome-Wide Association Study on Root System Architecture and Identification of Candidate Genes in Wheat (Triticum aestivum L.)

1
College of Life Science, Henan Normal University, Xinxiang 453007, China
2
State Key Laboratory of Crop Stress Biology for Arid Areas, Northwest A&F University, Yangling, Xianyang 712100, China
3
Key Laboratory of Crop Gene Resources and Germplasm Enhancement, Ministry of Agriculture, Institute of Crop Sciences, Chinese Academy of Agricultural Sciences, Beijing 100081, China
*
Authors to whom correspondence should be addressed.
Int. J. Mol. Sci. 2022, 23(3), 1843; https://doi.org/10.3390/ijms23031843
Submission received: 8 January 2022 / Revised: 2 February 2022 / Accepted: 4 February 2022 / Published: 6 February 2022
(This article belongs to the Special Issue Wheat Genetics and Genomics)

Abstract

:
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 −log10 (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.

1. Introduction

Bread wheat (Triticum aestivum L.) is one of the most important staple crops worldwide. The cultivated area of wheat is approximately 217 million hectares, annually producing approximately 765 million tonnes [1]. However, the demand for wheat continues to increase owing to a rapid increase in population. According to the speculation by the Food and Agriculture Organisation of the United Nations, approximately 840 million tonnes of wheat will be required annually by 2050 [2]. However, wheat production is facing serious challenges with the increasing global climate change in recent years [3,4]. Wheat breeding for high yield and stress tolerance will remain the main strategy for a long time. Roots, the underground parts of plants, absorb nearly all of the nutrients required for plant development and are primarily exposed to osmotic stress [5]. Optimal RSA improves the wheat yield and stress tolerance [6,7]. Therefore, it is important to study the genetic mechanisms of crop RSA traits to optimise the root system for further breeding.
RSA traits include root length, root surface area, root diameter, root tip number, and root volume [8]. RSA traits play important roles in determining crop development [9,10], yield formation [11,12], and stress tolerance [13]. For marker-assisted selection and key gene identification, the QTLs for RSA traits were explored in different crops [14,15]. Zhang et al. [16] analysed the seminal root length in rice using a recombinant inbred line (RIL) population of 150 lines under different water supply conditions, and identified two QTLs that explained 10% and 11% of phenotypic variation. Li et al. [17] identified 46 environment-specific QTLs for the crown root traits in maize, among which 12 QTLs were detected in multiple environments. Wu et al. [18] identified two major QTLs of swollen roots in turnips, providing a new idea regarding the molecular mechanism underlying storage root formation. In addition, Ren et al. [19] collated the QTLs for RSA and other traits, and found overlaps between QTLs for the root traits and those for nutrient uptake and productivity in many crop species, suggesting that RSA traits are closely related to crop development. These researches provided abundant information for further marker-assisted selection. Steele et al. [20] constructed introgressed lines using the previously reported root QTLs, which significantly improved the rice yield. Therefore, studying RSA traits is of great significance.
Bread wheat is consumed by 30–40% of humans [1]. Many researchers have paid attention to wheat root tissues. Manschadi et al. [10] analysed drought adaptation using wheat seedlings, and revealed that the root tissues are more closely related to drought tolerance than the aboveground tissues. Li et al. [21] found that deep roots, optimal root length density and xylem diameter, and increased root surface area contribute to drought tolerance. To elucidate the mechanisms of root development, many studies on QTL mapping of RSA traits have been performed in wheat. Liu et al. [22] analysed RSA traits using a diverse panel of 165 wheat cultivars, and identified 28 and 4 QTLs in hydroponic and field environments, respectively. Alemu et al. [23] analysed RSA traits using 192 Ethiopian durum wheat accessions, and identified 4 major (–log10p ≥ 4) and 34 nominal (–log10p ≥ 3.5) QTLs. Fan et al. [24] analysed the RSA traits in the seedlings of 188 RILs subjected to different nitrogen concentrations, and detected 53 QTLs responsible for RSA traits. In total, 26 published reports performed QTL exploration for RSA traits [5,19,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45], which provided ample information on RSA traits in wheat.
The wheat genome is enormous and complex. Research on the quantitative traits in wheat lags behind the research on other crops [28]. Especially for the RSA traits, each trait is controlled by multiple genes. Although researchers identified many QTLs for RSA traits, a few genes were functionally verified to regulate root development using the forward genetic approaches in wheat [46,47]. The mechanism for the RSA traits remains unclear, and more loci need to be explored in wheat. For the completion of genome sequencing, many types of SNP chips were developed, and the 660K SNP chip was considered to be the optimal chip for GWAS in wheat [48]. In this study, a diversity panel of 388 wheat accessions, based on pre-diversity assessment, was used to analyse RSA traits at the seedling stage under two environments. The 660K SNP chip was used for GWAS. We expect to identify vital QTLs and candidate genes to improve RSA.

2. Results

2.1. Population Structure and Linkage Disequilibrium (LD) Decay Analysis

The population structure was determined using STRUCTURE V2.3.4 according to the ∆K method of Bayesian clustering. We determined the slope breakpoint at K = 8. Thus, the panel used in this study was divided into eight subpopulations (Sp1-Sp8, Figure 1a, Table S1). Sp1 comprised 11 exotic varieties, 6 landraces, and 49 local cultivars. Sp2 consisted of 7 exotic varieties, 4 landraces, and 29 local cultivars. Sp3 comprised 5 exotic varieties, 1 landrace, and 58 local cultivars. Sp4 contained 13 exotic varieties, 15 landraces, and 19 local cultivars. Sp5 was composed of 10 exotic varieties, 3 landraces, and 28 local cultivars. Sp7 consisted of 1 exotic variety, 1 landrace, and 28 local cultivars. Sp8 contained 19 exotic varieties, 7 landraces, and 40 local cultivars (Table 1, Figure 1b).
We estimated the degree of LD using 161,801, 184,658, and 58,479 SNP markers in the A, B, and D sub-genomes, respectively. LD was analysed based on 717,701,068 pair-wise comparisons of 404,938 SNPs. Pair-wise LD was estimated using the squared allele-frequency correlation (r2). A plot of the LD estimates (r2) as a function of physical distance (Mb) indicated clear LD decay with increasing physical distance (Figure 1c). A comparison of LD between the sub-genomes and chromosomes showed variances in the LD decay. Overall, the average LD decay distance in the whole genome was approximately 3.6 Mb. LD decayed faster in the D genome (1.6 Mb) than in the A (3 Mb) and B (5.6 Mb) genomes (Figure 1c). The faster LD decay in the D genome is compatible with the evolutionary history of wheat.

2.2. Phenotypic Variation of the RSA Traits

A total of 4656 (388 wheat accessions × 6 replications × 2 environments) seedling roots were collected and observed for the nine RSA traits namely, total root length (TRL), total root surface area (TRA), average root diameter (ARD), total root volume (TRV), number of root tips (NRT), root length (RL), root surface area (RA), root volume (RV), and root tips (RT) (Table S1). The frequency distribution of the investigated traits was normal, as indicated by the kurtosis (bk) and skewness (bs) (Figure 2). This means that the traits were quantitative and suitable for GWAS. Phenotypic variation analysis showed that the nine RSA traits exhibited greater variation under OPC than under IPC, indicated by the coefficient of variation (Table 2). In addition, the nine RSA traits were significantly higher under OPC than under IPC (p ≤ 0.001). In other words, the root development was faster under OPC than that under IPC. The broad-sense heritability (H2) of the nine RSA traits ranged from 66.78% to 77.83%, indicating that the variations observed in these RSA traits were mainly regulated by genotype.

2.3. Correlations between RSA Traits

The relationships of the nine RSA traits between the two environments were analysed. The results showed that each RSA trait showed a significant correlation (p ≤ 0.001) between the two environments. This means that these traits were mainly regulated by genotype. The relationships between the nine traits were calculated (Figure 3). The results showed that significant positive correlations were observed between eight RSA traits, except for ARD. ARD showed a significant negative correlation with NRT, RL, RA, RV, and RT in the two environments. The correlation coefficients between RL, RA, and RV were higher than 0.94 in the two environments, indicating that these three traits were strongly correlated.

2.4. GWAS of RSA Traits

After filtering low-quality SNPs (minor allele frequency < 0.05 and missing data > 0.1), a total of 411,605 polymorphic SNP markers with effective chromosome information were available for GWAS using a univariate mixed linear model of GEMMA (Table S2). Of these, 404,938 SNPs were mapped to 21 wheat chromosomes, with 161,801, 184,658, and 58,479 SNPs in the A, B, and D sub-genomes, respectively. The marker density varied among different chromosomes (Figure S1, Table 3). We discovered that the minimum marker density was 8.1 markers per Mb on chromosome 4D, whereas the maximum was 55.09 markers per Mb on chromosome 3B. Compared with the A and B sub-genomes, the D sub-genome had fewer SNP markers and effective markers, especially chromosome 4D, which contained 1702 SNPs. The markers of known chromosome positions were used to analyse the genetic diversity. We found that there was little difference in the genetic diversity among the three sub-genomes. The mean values of the expected heterozygosity (He) and polymorphism information content (PIC) were 0.66–0.74 and 0.27–0.30 among the three sub-genomes, respectively.
We performed a GWAS using the best linear unbiased predictions (BLUPs) of the nine RSA traits. Significant SNPs for each trait were selected with a threshold of –log10 (p-value) ≥ 3.5 (Table S3, Figure S2). After analysing the data, 133 QTLs evenly distributed on 21 chromosomes were identified for the nine RSA traits, among which 38 were associated with multiple RSA traits. Eight QTLs were overlapped with previously identified QTLs (Table S3).
We focused on significant SNPs with a threshold of –log10 (p-value) ≥ 4. For convenience, we selected tagged SNPs for each QTL exhibiting the strongest association with RSA traits (Table 4), yielding 36 QTLs on 17 chromosomes (Figure 4), of which 13 were associated with more than one RSA trait. Among these QTLs, five QTLs, mapped on chromosomes 5B, 6B, 6D, 7A, and 7B, were co-localised with those identified in previous studies (Table 4), indicating that our results were reliable.

2.5. QTL Clusters for Root System Architecture Traits

QTL clusters refer to the QTLs that were associated with more than one RSA trait. We identified 13 QTL clusters located on chromosomes 2A, 2B, 3A, 3B, 3D, 4A, 5B, 5D, 6B, and 7B. Chromosomes 2A, 3B, and 7B contained two QTL clusters, and chromosomes 2B, 3A, 3D, 4A, 5B, 5D, and 6B harboured a single QTL cluster. One QTL cluster on chromosome 5B (571.23 Mb) was associated with five RSA traits: RL, RA, RV, RT, and NRT. One QTL cluster on chromosome 3D (549.74 Mb) was associated with four RSA traits: RL, RA, RV, and TRL. Two QTL clusters on chromosome 2A (61.28 Mb and 740.01 Mb) and one QTL cluster on chromosome 2B (694.26 Mb) were associated with TRA and TRV. Two QTL clusters on chromosome 7B were associated with TRL and TRA. One QTL cluster on chromosome 3A, two on chromosome 3B, one on chromosome 4A, one on chromosome 5D, and one on chromosome 6B were associated with RT and NRT. Therefore, these QTL clusters showed a pleiotropic effect.

2.6. Candidate Genes Associated with RSA

A total of 649 genes were identified in the 13 QTL clusters (Table S4). The genes encoding the bHLH transcription factor, WUSCHEL-RELATED HOMEOBOX 4 (WOX4), ROOT HAIR DEFECTIVE 6 (RHD6), GLABRA 3 (GL3), and ROOTHAIRLESS-LIKE 1 (LRL1), which are closely related to RSA traits, were found in the QTL clusters. We analysed the significant SNPs in thirteen QTL clusters, and identified nine genes with non-synonymous polymorphisms across the wheat panel (Table S5). The phenotypic variation related to each haplotype of the above nine genes was characterised, and the phenotypic differences between the haplotypes of two genes, TraesCS2A01G516200 with the SNP AX-110564036 (GG/AA) and TraesCS7B01G036900 with the SNP AX-112287343 (CC/TT), were significant (p ≤ 0.05) in more than one RSA trait.
TraesCS2A01G516200 (AX-110564036), mapped to chromosome 2A between 739.8 Mb and 740.8 Mb, encoded a short-chain dehydrogenase/reductase (Figure 5a). AX-110564036 SNP caused an amino acid change from isoleucine (Ile) to valine (Val) at 167 bp in this gene (Figure 5b). The 388 accessions were divided into two haplotypes (GG or AA) according to their genotype at this SNP. We determined that 258 and 93 accessions harboured the GG and AA haplotypes, respectively, whereas the remaining 37 accessions were heterozygous or lacked genotype information. Under the two environments, the accessions carrying the AA allele showed significantly higher values of the RSA traits TRL, TRA, and NRT than those carrying the GG allele (p ≤ 0.05), indicating that the AA haplotype had a positive effect on RSA (Figure 5c). RT-qPCR analysis showed that this gene was preferentially expressed in the root tissues at the jointing and heading stages (Figure 5d). These results showed that TraesCS2A01G516200 plays an important role in the RSA traits of wheat.
TraesCS7B01G036900 (AX-112287343), mapped to chromosome 7B between 35 Mb and 36 Mb, encoded a development and cell death domain protein (Figure 6a). AX-112287343 SNP caused an amino acid change from valine (Val) to alanine (Ala) at 396 bp in this gene (Figure 6b). The 388 accessions were divided into two haplotypes (CC or TT) according to their genotype at this SNP. We determined that 348 and 36 accessions harboured the CC and TT haplotypes, respectively, while the remaining four accessions were heterozygous or lacked genotype information. Under the two environments, the accessions with the CC haplotype exhibited markedly higher RSA trait values for TRL, TRA, TRV, RL, RA, and RV than those with the TT haplotype (p ≤ 0.05), indicating that the CC haplotype had a positive effect on RSA (Figure 6c). The RT-qPCR analysis indicated that this gene was preferentially expressed in the root tissues at the jointing and heading stages (Figure 6d). These results showed that TraesCS7B01G036900 plays an important role in the RSA traits of wheat.

3. Discussion

3.1. Soil Culture and 660K SNp Chip Are Ideal Methods for GWAS of RSA Traits in Wheat

Bread wheat is a staple crop worldwide, whose RSA traits strongly affect crop development and yield [49]. The spatial distribution of roots affects the absorption of water and nutrients and plant development. Zhu et al. [50] established that the length and number of seminal roots are particularly vital for the uptake of immobile nutrients. Devaiah et al. [51] found that the differences between RSA traits lead to changes in nutrient absorption. In addition, a deeper root system has higher water absorption capacity owing to its better access to soil water [52]. More importantly, White et al. [53] indicated that shoot and root architectures are inherited independently, and could be improved to increase the absorption and utilisation of water and nutrients. Therefore, it is important to analyse the RSA traits and explore the QTLs for these characteristics for further screening of key genes.
Some studies have identified several QTLs for RSA traits in wheat. To conveniently obtain whole roots, plants were usually cultivated in Hoagland’s nutrient solution, sand, or agarose gel [45,54]. Hydroponics is the most common approach used to obtain the whole roots in wheat [22,31,32,42]. Although the QTLs for RSA traits that were identified through hydroponic cultivation coincided with those for yield and nutrient absorption in wheat [43,44,55], we were not sure whether the phenotype of RSA traits in hydroponics was similar to that in soil cultivation. Therefore, soil cultivation is reliable to investigate RSA traits, although it is time-consuming and difficult. In this study, 388 wheat accessions were cultivated in soil under OPC and IPC, and the whole roots were sampled and cleaned up at the three-leaf stage to obtain the RSA traits for further analysis. The nine RSA traits showed high phenotypic variation at the seeding stage, indicating that the wheat panel used in this study exhibited high genetic diversity. The nine RSA traits under the two environments were significantly correlated, and the correlation between these RSA traits was significant, indicating that some QTLs may contribute to more than one RSA trait. The mean values of these RSA traits under OPC were significantly higher than those under IPC. We found that the soil under the two environments was comparable, and the mean temperature under IPC was higher than that under OPC. However, the root development was faster under OPC than that under IPC. It may be due to the climatic conditions from 15 October 2019 to 9 November 2019 that were optimal for wheat seedlings in Xinxiang (35.29° N, 113.83° E), Henan, China. The broad-sense heritability of the nine RSA traits ranged from 66.78–77.83%, indicating that these traits were suitable for QTL exploration.
Because of the vastness of genomic data, it is expensive to explore QTLs through genome re-sequencing in wheat. A series of high-density SNP chips, including 820K [56], 50K [57], 55K [58], 90K [33], and 660K chips, were designed and utilised for marker-assisted breeding in bread wheat. The 660K SNP chip contains 229,266 SNPs in the gene or promoter regions of 66,834 genes, representing 63.52% of the annotated genes in wheat. It is an optimal chip for marker-assisted breeding and has been widely used for QTL exploration for quality traits, agronomic traits, and disease resistance in bread wheat [48]. In the present study, the wheat panel was genotyped for GWAS using the 660K SNP chip. It is reliable to screen for significant QTLs associated with RSA traits.

3.2. QTL Analysis for RSA Traits

A substantial number of studies have performed QTL mapping for RSA traits in different plants, of which 26 studies focused on RSA traits in wheat. They identified 410 QTLs for RSA traits (Table 5). In this study, 36 QTLs with a threshold of −log10 (p-value) ≥ 4 were identified for nine RSA traits. Among these, five QTLs overlapped with previously reported QTLs, and 31 new QTLs were identified. This indicates that these genomic regions for the RSA traits are reliable and provide new information for further studies.
This study identified thirteen QTL clusters associated with more than one RSA trait that were distributed on ten chromosomes, indicating that there were co-localisations among these RSA traits. From the relationship among the nine RSA traits, we found that eight RSA traits showed significant positive correlations, and ARD showed a significant negative correlation with NRT, RL RA, RV, and RT in the two environments. From these results, we found that root traits are closely associated with each other, and some traits could be controlled by the same QTL. This might be because of the pleiotropic effects of the same gene or some RSA traits that might have a common genetic basis.

3.3. QTL Analysis Reveals the Key Genes for RSA

Key genes regulate RSA and root development. Therefore, the screening of key genes is important for RSA improvement. Zhang et al. [59] found that WOX4, encoded by AT1G46480, plays a central role in cambium development in Arabidopsis roots. Here, we found TraesCS2A01G514000, a homologue of AT1G46480, in the QTL on chromosome 2A, which was associated with TRA and TRV. Menand et al. [60] found that RHD6, encoded by AT1G66470, was a positive regulator for root hair development. TraesCS3B01G007700 was found to be homologous with AT1G66470 and mapped to the QTL of chromosome 3B that was associated with NRT and RT. Bernhardt et al. [61] found that a GL3 protein, encoded by AT5G41315 in Arabidopsis, partially regulates root epidermal cell fates. Here, we found that TraesCS2D01G575600, a homologue of AT5G41315, was located in the QTL of chromosome 2D that was associated with ARD, which might influence root epidermal cell fate. Karas et al. [62] found that LRL1, encoded by AT2G24260, is a partial regulator of root hair development in Arabidopsis. Here, we found that TraesCS7D01G158300, a homologue of AT2G24260, was located in the QTL of chromosome 7D that was associated with ARD. Therefore, these genes were considered candidates for the respective RSA traits.
In addition, we manually screened the significant SNPs for haplotype analysis, and found that nine genes showed non-synonymous polymorphisms across the wheat panel. Two genes, TraesCS2A01G516200 (AX-110564036) and TraesCS7B01G036900 (AX-112287343), have attracted our attention. The different combinations of haplotypes for the two genes showed significant differences for RSA traits. Importantly, the two genes were highly expressed in the root tissues. Functional annotation analysis revealed that TraesCS2A01G516200 encoded a short-chain dehydrogenase/reductase and TraesCS7B01G036900 encoded a development and cell death domain protein. These results indicate that these two genes may be vital candidate genes for root architecture.

4. Materials and Methods

4.1. Plant Materials and Experimental Treatment

Based on pre-diversity assessment of more than 5000 wheat accessions using KASP markers or 660 K SNP chip, a diversity panel of 388 wheat accessions was collected, including 38 local landraces, 66 exotic cultivars, and 284 cross-derived cultivars in China. The name and origin of each wheat accession are listed in Table S1. Wheat seeds of uniform size were washed thrice with distilled water and germinated on wet filter papers in Petri dishes for germination at 20 ℃. After 2 days, nine germinating seeds of each accession were evenly planted in one pot, and each pot was thinned to six seedlings after 10 days of sowing. Each pot (height 15 cm; diameter 20 cm) was filled with 4.5 kg of soil (nutritional soil: sieved tillage soil, 1:1) to easily obtain the whole root system.
For OPC, the pots were randomly buried in a field on 15 October 2019 at Xinxiang (35.29° N, 113.83° E), Henan, China. To facilitate management, 10 blocks were set in the field, each containing approximately 40 pots. The plants were watered to 70–80% of the field capacity. The seedlings developed to the three-leaf stage after 25 days of sowing with a mean temperature of 14.85 °C and 0 mm of rainfall. For IPC, the pots were randomly placed in a phytotron. The seedlings were grown under 80% relative humidity, a 16 h light (20 °C) / 8 h darkness (18 °C) photoperiod, and 8000 Lux light intensity. The seedlings developed to the three-leaf stage after 30 days of sowing. The roots were cleaned with water. For each wheat accession under the two environments, six biological replications were sampled for further RSA trait measurement.

4.2. RSA Trait Measurement

The root samples were placed in plexiglass trays (200 mm × 250 mm) containing 4–6 mm of distilled water. The roots were manipulated to minimise tangling and overlap. Nine RSA traits (Table 6), TRL, TRA, ARD, TRV, NRT, RL, RA, RV, and RT, were measured using a recording scanner (Epson 1680, Suwa, Japan). The images were analysed using Win RHIZO Pro Vision 5.0a (Regent Instruments Inc., Quebec, QC, Canada), which was operationally semi-automated.

4.3. Phenotypic Data Analysis for RSA Traits

SPSS (Version 17.0) and R packages (R Version 4.1.1) were used for data analysis. Genotype, environments, and genotype-by-environment interactions were considered as random effects in a linear mixed model (LMM_4 K). The lme4 package in the R 4.1.1 program was used to estimate the BLUPs. The H2 of the RSA traits was calculated using the mean values of each experiment among six replications using the formula:
H 2 = σ G 2 σ G 2 + σ G E 2 n + σ e 2 n r    
where σ G 2 is the genetic variance, σ G E 2 is the genotype × environment, σ e 2 is the error variance, n represents the number of environments, and r represents the number of replications. The mean values were used to analyse the differences using a t-test. The CV was calculated for all the RSA traits [63].

4.4. Population Structure

Principal component analysis of the population was performed using the GCTA software [64]. For population structure analysis, a Bayesian model-based (Markov Chain Monte Carlo) clustering approach was used in the STRUCTURE software v2.3.4 with unlinked markers (r2 = 0) [65]. The K-value, representing the number of subgroups, was set from 2 to 20, and five independent runs for each K were performed with the burn-in length set to 20,000 and iterations set to 10,000. The most likely number of subpopulations was determined using the ∆K method based on the rate of change in LnP (D) between successive K-values [66].

4.5. LD Analysis

Genome-wide LD analysis for the A, B, and D genomes was performed using the PLINK software. The squared allele-frequency correlations (r2) were used to estimate the LD by applying pair-wise comparisons among the filtered SNP markers. The values for genomes A, B, and D were plotted against the genetic distance to determine LD decay. The parameters for calculating r2 were set to --r2 --ld-window-kb 30000 --ld-window 1000 --ld-window-r2 0. A locally weighted polynomial regression (LOESS) curve was drawn to fit the data using second-order local weighted scatter plot smoothing in the R program. The critical distance spanning the QTL was defined based on the intersection point of the fitted LOESS curve with the LD (r2 = 0.1).

4.6. Genotyping Data Using 660K SNp Chip

The wheat panel was genotyped using the Affymetrix Wheat 660K SNP Chip [48], which was provided by the Beijing CapitalBio Technology Company (http://www.capitalbiotech.com, accessed on 8 January 2022). SNP genotyping calling and allele clustering were performed using the polyploid version of the Affymetrix Genotyping Console software. To ensure genotyping accuracy, quality control, and filtering of raw data with missing data > 10%, a minor allele frequency < 5%, or Hardy-Weinberg Equilibrium > 0.01 were performed. The remaining high-quality SNP markers were used for GWAS. The sequence of each polymorphic SNP was identified using Basic Local Alignment Search Tool (BLAST)n 2.2.29+, and the physical positions of all SNPs were determined by blasting them with the wheat reference genome IWGSC RefSeq v1.0. The parameters of the blast result were set to: ratio (AL/QL) ≥ 0.95, identity ≥ 0.95, and E-value < 1 × e−10. Two genetic diversity parameters, PIC and He, were calculated for each SNP marker and each chromosome using the formulas described by Botsteinet et al. [67] and Nei [68], respectively.

4.7. GWAS

GWAS was conducted using a univariate linear mixed model in the GEMMA software [69]. The suggestive threshold for p-value (p = 1/Ne) was calculated using a modified Bonferroni correction (genetic type 1 Error calculator, version 0.2), where Ne represents the effective number of SNP markers [70]. A threshold was established at −log10p ≥ 4 to determine the marker–trait association significance. Moreover, the p-value of −log10p ≥ 3.5 was used to establish a second, less restrictive threshold. Significant markers from GWAS were visualised in Manhattan plots, and important p-value distributions were visualised using quantile-quantile (Q-Q) plots drawn with the qqman package in R 4.1.1. The phenotypic variance (R2) explained by significant SNPs was evaluated using the GCTA software [64].

4.8. Candidate Gene Identification, RNA Extraction, and RT-qPCR Verification

Candidate genes were selected using functional annotation and haplotype analysis. For haplotype analysis, the significant SNPs in QTL clusters were used to identify non-synonymous polymorphisms. The phenotypic differences between the haplotypes of the nine RSA traits were analysed to identify the candidate genes.
The root, stem, and leaf tissues were sampled from wheat plants of cv. Chinese Spring grown in the field at the jointing and heading stages, respectively. The total RNA was extracted using the TRIzol reagent (TaKaRa, Ohtsu, Japan) according to the manufacturer’s instructions. One microgram of the total RNA was used for first-strand cDNA synthesis using the HiScript III First-Strand cDNA Synthesis Kit (Vazyme, Nanjing, China). qPCR was performed on an ABI 7500 real-time PCR system (Applied Biosystems, Forest City, CA, USA) using SYBR Premix Ex Taq Kit (TaKaRa, Japan). Three biological replications were used for each experiment. TaActin3 was used as an internal control [71], and the relative expression value was calculated using the 2−∆∆Ct method. All primers used for RT-qPCR were listed in Table S6.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms23031843/s1.

Author Contributions

Conceptualization, J.M., T.L. and L.J.; investigation, J.M., D.Z. (Dongyang Zhao), X.T., M.Y., D.Z. (Daijing Zhang), M.X., Y.D. and H.R.; writing—original draft preparation, J.M., Q.Z., J.W. and D.H.; writing—review and editing, J.M., T.L. and L.J. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Key Research and Development Program of China (2017YFD0301101 and 2016YFD0300203-3).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ray, D.K.; Mueller, N.D.; West, P.C.; Foley, J.A. Yield Trends Are Insufficient to Double Global Crop Production by 2050. PLoS ONE 2013, 8, e66428. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Sharma, I.; Tyagi, B.S.; Singh, G.; Venkatesh, K.; Gupta, O.P. Enhancing wheat production-A global perspective. Indian J. Agric. Sci. 2015, 85, 3–13. [Google Scholar]
  3. Wheeler, T.; Braun, J.V. Climate Change Impacts on Global Food Security. Science 2013, 341, 508–513. [Google Scholar] [CrossRef] [PubMed]
  4. Zampieri, M.; Ceglar, A.; Dentener, F.; Toreti, A. Wheat yield loss attributable to heat waves, drought and water excess at the global, national and subnational scales. Environ. Res. Lett. 2017, 12, 064008. [Google Scholar] [CrossRef]
  5. Maccaferri, M.; El-Feki, W.; Nazemi, G.; Salvi, S.; Canè, M.A.; Colalongo, M.C.; Stefanelli, S.; Tuberosa, R. Prioritizing quantitative trait loci for root system architecture in tetraploid wheat. J. Exp. Bot. 2016, 67, 1161–1178. [Google Scholar] [CrossRef]
  6. Atkinson, J.A.; Pound, M.P.; Bennett, M.J.; Wells, D.M. Uncovering the hidden half of plants using new advances in root phenotyping. Curr. Opin. Biotechnol. 2019, 55, 1–8. [Google Scholar] [CrossRef]
  7. Li, X.K.; Guo, Z.L.; Lv, Y.; Cen, X.; Ding, X.P.; Wu, H.; Li, X.H.; Huang, J.P.; Xiong, L.Z. Genetic control of the root system in rice under normal and drought stress conditions by genome-wide association study. PLOS Genet. 2017, 13, e1006889. [Google Scholar] [CrossRef] [Green Version]
  8. Chen, X.X.; Ding, Q.S.; Błaszkiewicz, Z.; Sun, J.A.; Sun, Q.; He, R.Y.; Li, Y.N. Phenotyping for the dynamics of field wheat root system architecture. Sci. Rep. 2017, 7, 37649. [Google Scholar] [CrossRef] [Green Version]
  9. Huang, C.Y.; Kuchel, H.; Edwards, J.; Hall, S.; Parent, B.; Eckermann, P.; Herdina; Hartley, D.M.; Langridge, P.; McKay, A.C. A DNA-based method for studying root responses to drought in field-grown wheat genotypes. Sci. Rep. 2013, 3, 3194. [Google Scholar] [CrossRef] [Green Version]
  10. Manschadi, A.M.; Hammer, G.L.; Christopher, J.T.; Devoil, P. Genotypic variation in seedling root architectural traits and implications for drought adaptation in wheat (Triticum aestivum L.). Plant Soil 2008, 303, 115–129. [Google Scholar] [CrossRef]
  11. Voss-Fels, K.P.; Qian, L.; Parra-Londono, S.; Uptmoor, R.; Frisch, M.; Keeble-Gagnère, G.; Apples, R.; Snowdon, R.J. Linkage drag constrains the roots of modern wheat. Plant Cell Environ. 2017, 40, 717–725. [Google Scholar] [CrossRef] [PubMed]
  12. Man, J.; Shi, Y.; Yu, Z.; Zhang, Y. Dry Matter Production, Photosynthesis of Flag Leaves and Water Use in Winter Wheat Are Affected by Supplemental Irrigation in the Huang-Huai-Hai Plain of China. PLoS ONE 2015, 10, e0137274. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Wasaya, A.; Zhang, X.Y.; Fang, Q.; Yan, Z.Z. Root Phenotyping for Drought Tolerance: A Review. Agronomy 2018, 8, 241. [Google Scholar] [CrossRef] [Green Version]
  14. Christopher, J.; Christopher, M.; Jennings, R.; Jones, S.; Fletcher, S.; Borrell, A.; Manschadi, B.A.; Jordan, D.; Mace, E.; Hammer, G. QTL for root angle and number in a population developed from bread wheats (Triticum aestivum) with contrasting adaptation to water-limited environments. Theor. Appl. Genet. 2013, 126, 1563–1574. [Google Scholar] [CrossRef] [PubMed]
  15. Bai, C.H.; Liang, Y.L.; Hawkesford, M.J. Identification of QTLs associated with seedling root traits and their correlation with plant height in wheat. J. Exp. Bot. 2013, 64, 1745–1753. [Google Scholar] [CrossRef] [Green Version]
  16. Zhang, W.P.; Shen, X.Y.; Wu, P.; Hu, B.; Liao, C.Y. QTLs and epistasis for seminal root length under a different water supply in rice (Oryza sativa L.). Theor. Appl. Genet. 2001, 103, 118–123. [Google Scholar] [CrossRef]
  17. Li, P.C.; Fan, Y.Y.; Yin, S.Y.; Wang, Y.Y.; Wang, H.M.; Xu, Y.; Yang, Z.F.; Xu, C.W. Multi-environment QTL mapping of crown root traits in a maize RIL population. Crop J. 2020, 8, 645–654. [Google Scholar] [CrossRef]
  18. Wu, Y.D.; Zhang, S.F.; Zhang, H.; Li, F.; Li, G.L.; Fan, C.C.; Sun, R.F.; Zhang, S.J. QTL mapping and candidate gene identification of swollen root formation in Turnip. Int. J. Mol. Sci. 2021, 22, 653. [Google Scholar] [CrossRef]
  19. Ren, Y.; Qian, Y.; Xu, Y.; Zou, C.; Liu, D.; Zhao, X.; Zhang, A.; Tong, Y. Characterization of QTLs for Root Traits of Wheat Grown under Different Nitrogen and Phosphorus Supply Levels. Front. Plant Sci. 2017, 8, 2096. [Google Scholar] [CrossRef] [Green Version]
  20. Steele, K.A.; Price, A.H.; Witcombe, J.R.; Shrestha, R.; Singh, B.N.; Gibbons, J.M.; Virk, D.S. QTLs associated with root traits increase yield in upland rice when transferred through marker-assisted selection. Theor. Appl. Genet. 2012, 126, 101–108. [Google Scholar] [CrossRef]
  21. Li, C.N.; Li, L.; Reynolds, M.P.; Wang, J.Y.; Chang, X.P.; Mao, X.G.; Jing, R.L. Recognizing the hidden half in wheat: Root system attributes associated with drought tolerance. J. Exp. Bot. 2021, 72, 5117–5133. [Google Scholar] [CrossRef] [PubMed]
  22. Liu, P.; Jin, Y.R.; Liu, J.D.; Liu, C.Y.; Yao, H.P.; Luo, F.Y.; Guo, Z.H.; Xia, X.C.; He, Z.H. Genome-wide association mapping of root system architecture traits in common wheat (Triticum aestivum L.). Euphytica 2019, 215, 121. [Google Scholar] [CrossRef]
  23. Alemu, A.; Feyissa, T.; Maccaferri, M.; Sciara, G.; Tuberosa, R.; Ammar, K.; Badebo, A.; Acevedo, M.; Letta, L.; Abeyo, B. Genome-wide association analysis unveils novel QTLs for seminal root system architecture traits in Ethiopian durum wheat. BMC Genom. 2021, 22, 20. [Google Scholar] [CrossRef] [PubMed]
  24. Fan, X.L.; Zhang, W.; Zhang, N.; Chen, M.; Zheng, S.S.; Zhao, C.H.; Han, J.; Liu, J.J.; Zhang, X.L.; Song, L.Q.; et al. Identification of QTL regions for seedling root traits and their effect on nitrogen use efficiency in wheat (Triticum aestivum L.). Theor. Appl. Genet. 2018, 131, 2677–2698. [Google Scholar] [CrossRef]
  25. Salarpour, M.; Pakniyat, H.; Abdolshahi, R.; Heidari, B.; Razi, H.; Afzali, R. Mapping QTL for agronomic and root traits in the Kukri/RAC875 wheat (Triticum aestivum L.) population under drought stress conditions. Euphytica 2020, 216, 105. [Google Scholar] [CrossRef]
  26. Huang, F.; Chen, Z.Y.; Du, D.J.; Guan, P.F.; Chai, L.L.; Guo, W.L.; Hu, Z.R.; Xin, M.M.; Peng, H.R.; Yao, Y.Y.; et al. Genome-wide linkage mapping of QTL for root hair length in a Chinese common wheat population. Crop J. 2020, 8, 1049–1056. [Google Scholar] [CrossRef]
  27. Yang, M.J.; Wang, C.R.; Hassan, M.A.; Li, F.J.; Xia, X.C.; Shi, S.B.; Xiao, Y.G.; He, Z.H. QTL mapping of root traits in wheat under different phosphorus levels using hydroponic culture. BMC Genom. 2021, 22, 174. [Google Scholar] [CrossRef]
  28. Xu, F.D.; Chen, S.L.; Yang, X.W.; Zhou, S.M.; Wang, J.S.; Zhang, Z.L.; Huang, Y.; Song, M.; Zhang, J.; Zhan, K.H.; et al. Genome-Wide Association Study on Root Traits Under Different Growing Environments in Wheat (Triticum aestivum L.). Front. Genet. 2021, 12, 646712. [Google Scholar] [CrossRef]
  29. Li, Z.X.; Ni, Z.F.; Peng, H.R.; Liu, Z.Y.; Nie, X.L.; Xu, S.B.; Liu, G.; Sun, Q.X. Molecular mapping of QTLs for root response to phosphorus deficiency at seedling stage in wheat (Triticum aestivum L.). Prog. Nat. Sci. 2007, 17, 1177–1184. [Google Scholar]
  30. Yang, M.J.; Wang, C.R.; Hassan, M.A.; Wu, Y.Y.; Xia, X.C.; Shi, S.B.; Xiao, Y.G.; He, Z.H. QTL mapping of seedling biomass and root traits under different nitrogen conditions in bread wheat (Triticum aestivum L.). J. Integr. Agric. 2021, 20, 1180–1192. [Google Scholar] [CrossRef]
  31. Li, T.; Ma, J.; Zou, Y.Y.; Chen, G.D.; Ding, P.Y.; Zhang, H.; Lan, X.J. Quantitative trait loci for seeding root traits and the relationships between root and agronomic traits in common wheat. Genome 2019, 63, 27–36. [Google Scholar] [CrossRef]
  32. Zheng, X.W.; Wen, X.J.; Qiao, L.; Zhao, J.J.; Zhang, X.J.; Li, X.; Zheng, J. A novel QTL QTrl.saw-2D.2 associated with the total root length identified by linkage and association analyses in wheat (Triticum aestivum L.). Planta 2019, 250, 129–143. [Google Scholar] [CrossRef]
  33. Kabir, M.R.; Liu, G.; Guan, P.; Wang, F.; Khan, A.A.; Ni, Z.; Sun, Q. Mapping QTLs associated with root traits using two different populations in wheat (Triticum aestivum L.). Euphytica 2015, 206, 175–190. [Google Scholar] [CrossRef]
  34. Soriano, J.M.; Alvaro, F. Discovering consensus genomic regions in wheat for root-related traits by QTL meta-analysis. Sci. Rep. 2019, 9, 10537. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Roselló, M.; Royo, C.; Sanchez-Garcia, M.; Soriano, J.M. Genetic Dissection of the Seminal Root System Architecture in Mediterranean Durum Wheat Landraces by Genome-Wide Association Study. Agronomy 2019, 9, 364. [Google Scholar] [CrossRef] [Green Version]
  36. Iannucci, A.; Marone, D.; Russo, M.A.; De Vita, P.; Miullo, V.; Ferragonio, P.; Mastrangelo, A.M. Mapping QTL for Root and Shoot Morphological Traits in a Durum Wheat × T. dicoccum Segregating Population at Seedling Stage. Int. J. Genom. 2017, 2017, 1–17. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Ibrahim, S.E.; Schubert, A.; Pillen, K.; Léon, J. QTL analysis of drought tolerance for seedling root morphological traits in an advanced backcross population of spring wheat. Int. J. Agric. Sci. 2012, 2, 619–629. [Google Scholar]
  38. Xie, Q.; Fernando, K.M.C.; Mayes, S.; Sparkes, D.L. Identifying seedling root architectural traits associated with yield and yield components in wheat. Ann. Bot. 2017, 119, 1115–1129. [Google Scholar] [CrossRef] [PubMed]
  39. Ma, J.; Luo, W.; Zhang, H.; Zhou, X.H.; Qin, N.N.; Wei, Y.M.; Liu, Y.X.; Jiang, Q.T.; Lan, X.J. Identification of quantitative trait loci for seedling root traits from Tibetan semi-wild wheat (Triticum aestivum subsp. tibetanum). Genome 2017, 60, 1068–1075. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Atkinson, J.A.; Wingen, L.U.; Griffiths, M.; Pound, M.P.; Gaju, O.; Foulkes, M.J.; Wells, D.M. Phenotyping pipeline reveals major seedling root growth QTL in hexaploid wheat. J. Exp. Bot. 2015, 66, 2283–2292. [Google Scholar] [CrossRef] [Green Version]
  41. Merchuk-Ovnat, L.; Fahima, T.; Ephrath, J.E.; Krugman, T.; Saranga, Y. Ancestral QTL Alleles from Wild Emmer Wheat Enhance Root Development under Drought in Modern Wheat. Front. Plant Sci. 2017, 8, 1–12. [Google Scholar] [CrossRef] [PubMed]
  42. Halder, T.; Liu, H.; Chen, Y.L.; Yan, G.J.; Siddique, K.H. Identification of Candidate Genes for Root Traits Using Genotype-Phenotype Association Analysis of Near-Isogenic Lines in Hexaploid Wheat (Triticum aestivum L.). Int. J. Mol. Sci. 2021, 22, 3579. [Google Scholar] [CrossRef]
  43. Canè, M.A.; Maccaferri, M.; Nazemi, G.; Salvi, S.; Francia, R.; Colalongo, C.; Tuberosa, R. Association mapping for root architectural traits in durum wheat seedlings as related to agronomic performance. Mol. Breed. 2014, 34, 1629–1645. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Horn, R.; Wingen, L.U.; Snape, J.W.; Dolan, L. Mapping of quantitative trait loci for root hair length in wheat identifies loci that co-locate with loci for yield components. J. Exp. Bot. 2016, 67, 4535–4543. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Liu, X.L.; Li, R.Z.; Chang, X.P.; Jing, R.L. Mapping QTLs for seedling root traits in a doubled haploid wheat population under different water regimes. Euphytica 2012, 189, 51–66. [Google Scholar] [CrossRef]
  46. Kai, V.F.; Robinson, H.; Mudge, S.R.; Richard, C.; Hickey, L.T. Vernalization1 modulates root system architecture in wheat and barley. Mol. Plant 2018, 11, 226. [Google Scholar] [CrossRef] [Green Version]
  47. Wang, J.; Li, L.; Li, C.; Yang, X.; Xue, Y.; Zhu, Z.; Mao, X.; Jing, R. A transposon in the vacuolar sorting receptor gene TaVSR1-B promoter region is associated with wheat root depth at booting stage. Plant Biotechnol. J. 2021, 19, 1456–1467. [Google Scholar] [CrossRef]
  48. Sun, C.W.; Dong, Z.D.; Zhao, L.; Ren, Y.; Zhang, N.; Chen, F. The Wheat 660K SNP array demonstrates great potential for marker-assisted selection in polyploid wheat. Plant Biotechnol. J. 2020, 18, 1354–1360. [Google Scholar] [CrossRef]
  49. Fang, C.; Ma, Y.M.; Wu, S.W.; Liu, Z.; Wang, Z.; Yang, R.; Hu, G.H.; Zhou, Z.K.; Yu, H.; Tian, Z.X. Genome-wide association studies dissect the genetic networks underlying agronomical traits in soybean. Genome Biol. 2017, 18, 1–4. [Google Scholar] [CrossRef]
  50. Zhu, J.; Mickelson, S.M.; Kaeppler, S.M.; Lynch, J.P. Detection of quantitative trait loci for seminal root traits in maize (Zea mays L.) seedlings grown under differential phosphorus levels. Theor. Appl. Genet. 2006, 113, 1–10. [Google Scholar] [CrossRef]
  51. Devaiah, B.N.; Nagarajan, V.K.; Raghothama, K.G. Phosphate Homeostasis and Root Development in Arabidopsis Are Synchronized by the Zinc Finger Transcription Factor ZAT6. Plant Physiol. 2007, 145, 147–159. [Google Scholar] [CrossRef] [Green Version]
  52. Tambussi, E.A.; Bort, J.; Araus, J.L. Water use efficiency in C3 cereals under Mediterranean conditions: A review of physiological aspects. Ann. Appl. Biol. 2007, 150, 307–321. [Google Scholar] [CrossRef]
  53. White, P.J. Root traits benefitting crop production in environments with limited water and nutrient availability. Ann. Bot. 2019, 124, 883–890. [Google Scholar] [CrossRef]
  54. Ayalew, H.; Desalegn, T.; Hui, L.; Yan, G. Performance of Ethiopian bread wheat (Tritium aestivum L.) genotypes under contrasting water regimes: Potential sources of variability for drought resistance breeding. Aust. J. Crop Sci. 2016, 10, 370–376. [Google Scholar] [CrossRef]
  55. An, D.G.; Su, J.Y.; Liu, Q.Y.; Li, B.; Jing, R.L.; Li, B.; Li, Z.Z. Mapping QTLs for nitrogen uptake in relation to the early growth of wheat (Triticum aestivum L.). Plant Soil 2006, 284, 73–84. [Google Scholar] [CrossRef]
  56. Winfield, M.O.; Allen, A.M.; Burridge, A.J.; Barker, G.L.A.; Benbow, H.R.; Wilkinson, P.A.; Edwards, K.J. High-density SNP genotyping array for hexaploid wheat and its secondary and tertiary gene pool. Plant Biotechnol. J. 2015, 14, 1195–1206. [Google Scholar] [CrossRef]
  57. Rasheed, A.; Hao, Y.F.; Xia, X.C.; Khan, A.; Xu, Y.B.; Varshney, R.K.; He, Z.H. Crop Breeding Chips and Genotyping Platforms: Progress, Challenges, and Perspectives. Mol. Plant 2017, 10, 1047–1064. [Google Scholar] [CrossRef] [Green Version]
  58. Ye, X.L.; Li, J.; Cheng, Y.K.; Yao, F.J.; Long, L.; Wang, Y.Q.; Chen, G.Y. Genome-wide association study reveals new loci for yield-related traits in Sichuan wheat germplasm under stripe rust stress. BMC Genom. 2019, 20, 640. [Google Scholar] [CrossRef] [Green Version]
  59. Zhang, J.; Eswaran, G.; Alonso-Serra, J.; Kucukoglu, M.; Xiang, J.; Yang, W.; Helariutta, Y. Transcriptional regulatory framework for vascular cambium development in Arabidopsis roots. Nat. Plants 2019, 5, 1033–1042. [Google Scholar] [CrossRef]
  60. Menand, B.; Yi, K.; Jouannic, S.; Hoffmann, L.; Ryan, E.; Linstead, P.; Dolan, L. An Ancient Mechanism Controls the Development of Cells with a Rooting Function in Land Plants. Science 2007, 316, 1477–1480. [Google Scholar] [CrossRef] [Green Version]
  61. Bernhardt, C. The bHLH genes GLABRA3 (GL3) and ENHANCER OF GLABRA3 (EGL3) specify epidermal cell fate in the Arabidopsis root. Development 2003, 130, 6431–6439. [Google Scholar] [CrossRef] [Green Version]
  62. Karas, B.; Amyot, L.; Johansen, C.; Sato, S.; Tabata, S.; Kawaguchi, M.; Szczyglowski, K. Conservation of Lotus and Arabidopsis Basic Helix-Loop-Helix Proteins Reveals New Players in Root Hair Development. Plant Physiol. 2009, 151, 1175–1185. [Google Scholar] [CrossRef] [Green Version]
  63. Nyquist, W.E.; Baker, R.J. Estimation of heritability and prediction of selection response in plant populations. Crit. Rev. Plant Sci. 1991, 10, 235–322. [Google Scholar] [CrossRef]
  64. Yang, J.; Lee, S.H.; Goddard, M.E.; Visscher, P.M. GCTA: A Tool for Genome-wide Complex Trait Analysis. Am. J. Hum. Genet. 2011, 88, 76–82. [Google Scholar] [CrossRef] [Green Version]
  65. Pritchard, J.K.; Stephens, M.; Donnelly, P. Inference of Population Structure Using Multilocus Genotype Data. Genetics 2000, 155, 945–959. [Google Scholar] [CrossRef]
  66. Earl, D.A.; VonHoldt, B.M. STRUCTURE HARVESTER: A website and program for visualizing STRUCTURE output and implementing the Evanno method. Conserv. Genet. Resour. 2011, 4, 359–361. [Google Scholar] [CrossRef]
  67. Botstein, D.; White, R.L.; Skolnick, M.; Davis, R.W. Construction of a genetic linkage map in man using restriction fragment length polymorphisms. Am. J. Hum. Genet. 1980, 32, 314. [Google Scholar] [CrossRef]
  68. Nei, M. Estimation of average heterozygosity and genetic distance from a small number of individuals. Genetics 1978, 89, 583–590. [Google Scholar] [CrossRef]
  69. Zhou, X.; Stephens, M. Genome-wide efficient mixed-model analysis for association studies. Nat. Genet. 2012, 44, 821–824. [Google Scholar] [CrossRef] [Green Version]
  70. Li, M.X.; Yeung, J.M.Y.; Cherny, S.S.; Sham, P.C. Evaluating the effective numbers of independent tests and significant p-value thresholds in commercial genotyping arrays and public imputation reference datasets. Hum. Genet. 2011, 131, 747–756. [Google Scholar] [CrossRef] [Green Version]
  71. Wang, X.; Ren, Y.; Li, J.; Wang, Z.; Xin, Z.; Lin, T. Knock-down the expression of TaH2B-7D using virus-induced gene silencing reduces wheat drought tolerance. Biol. Res. 2019, 52, 14. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Population structure and LD analysis of the wheat panel. (a) Subpopulations inferred by K-means structure analysis. (b) Principal component analysis of all wheat accessions. (c) Genome-wide average LD decay over physical distance. Pair-wise single-nucleotide polyLD (r2) values based on the physical positions from the IWGSC RefSeq v.1.0 reference genome are plotted as a function of mapping distance (Mb) between markers. The different-coloured curves represent LD decay fits for sub-genomes A (purple), B (blue), D (green), and the whole genome (red). The thick horizontal black line represents the population-specific critical r2 value (0.1) above which LD may be due to linkage.
Figure 1. Population structure and LD analysis of the wheat panel. (a) Subpopulations inferred by K-means structure analysis. (b) Principal component analysis of all wheat accessions. (c) Genome-wide average LD decay over physical distance. Pair-wise single-nucleotide polyLD (r2) values based on the physical positions from the IWGSC RefSeq v.1.0 reference genome are plotted as a function of mapping distance (Mb) between markers. The different-coloured curves represent LD decay fits for sub-genomes A (purple), B (blue), D (green), and the whole genome (red). The thick horizontal black line represents the population-specific critical r2 value (0.1) above which LD may be due to linkage.
Ijms 23 01843 g001
Figure 2. Distribution frequency of the RSA traits of 388 wheat accessions under OPC. The black curve represents the density curve.
Figure 2. Distribution frequency of the RSA traits of 388 wheat accessions under OPC. The black curve represents the density curve.
Ijms 23 01843 g002
Figure 3. Correlation analysis of the nine RSA traits in 388 wheat accessions under OPC (a) and IPC (b). Intensities of the blue and red colours indicate the degree of positive and negative correlations, respectively. Circle size indicates low to high significance.
Figure 3. Correlation analysis of the nine RSA traits in 388 wheat accessions under OPC (a) and IPC (b). Intensities of the blue and red colours indicate the degree of positive and negative correlations, respectively. Circle size indicates low to high significance.
Ijms 23 01843 g003
Figure 4. Location of the QTLs for the RSA traits (−log10 [p-value] ≥ 4).
Figure 4. Location of the QTLs for the RSA traits (−log10 [p-value] ≥ 4).
Ijms 23 01843 g004aIjms 23 01843 g004b
Figure 5. GWAS results of the nine RSA traits and identification of candidate gene for the SNP AX-110564036 on chromosome 2A. (a) Manhattan plot for the candidate region on chromosome 2A. The green dot represents the significant SNP AX-110564036. The LD block analysis of the SNPs in this region is shown below. The degree of linkage is represented by the coefficient of r2. (b) Gene structure of TraesCS2A01G516200 and its non-synonymous SNP. Yellow rectangles and black lines represent exons and introns, respectively. The arrow indicates the position of amino acid variations in TraesCS2A01G516200. (c) Box plots for RSA traits based on the allele at SNP AX-110564036. Differences between the haplotypes were statistically analysed using the Student’s t-test. The upper and lower edges of the box represent the 75th and 25th quantiles, respectively, and the whiskers show the 90th and 10th quantiles; the horizontal solid line represents the average; * Significant at p ≤ 0.05, *** Significant at p ≤ 0.001. (d) Expression profile of TraesCS2A01G516200 in different tissues, as determined by RT-qPCR. The data are the means ± Sd of n = 3. ‘ns’ means ‘not significant’.
Figure 5. GWAS results of the nine RSA traits and identification of candidate gene for the SNP AX-110564036 on chromosome 2A. (a) Manhattan plot for the candidate region on chromosome 2A. The green dot represents the significant SNP AX-110564036. The LD block analysis of the SNPs in this region is shown below. The degree of linkage is represented by the coefficient of r2. (b) Gene structure of TraesCS2A01G516200 and its non-synonymous SNP. Yellow rectangles and black lines represent exons and introns, respectively. The arrow indicates the position of amino acid variations in TraesCS2A01G516200. (c) Box plots for RSA traits based on the allele at SNP AX-110564036. Differences between the haplotypes were statistically analysed using the Student’s t-test. The upper and lower edges of the box represent the 75th and 25th quantiles, respectively, and the whiskers show the 90th and 10th quantiles; the horizontal solid line represents the average; * Significant at p ≤ 0.05, *** Significant at p ≤ 0.001. (d) Expression profile of TraesCS2A01G516200 in different tissues, as determined by RT-qPCR. The data are the means ± Sd of n = 3. ‘ns’ means ‘not significant’.
Ijms 23 01843 g005
Figure 6. GWAS results of the nine RSA traits and identification of the candidate gene for the SNP AX-112287343 on chromosome 7B. (a) Manhattan plot for the candidate region on chromosome 7B. The green dot represents the significant SNP AX-112287343. The LD block analysis of the SNPs in this region is shown below. The degree of linkage is represented by the coefficient of r2. (b) Gene structure of TraesCS7B01G036900 and its non-synonymous SNP. Yellow rectangles and black lines represent exons and introns, respectively. The arrow indicates the position of amino acid variations in TraesCS7B01G036900. (c) Box plots for RSA traits based on the allele at SNP AX-112287343. Differences between the haplotypes were statistically analysed using the Student’s t-test. The upper and lower edges of the box represent the 75th and 25th quantiles, respectively, and the whiskers show the 90th and 10th quantiles; the horizontal solid line represents the average; * Significant at p ≤ 0.05, ** Significant at p ≤ 0.01, *** Significant at p ≤ 0.001. (d) Expression profile of TraesCS7B01G036900 in different tissues, as determined by RT-qPCR. The data are the means ± Sd of n = 3. ‘ns’ means ‘not significant’.
Figure 6. GWAS results of the nine RSA traits and identification of the candidate gene for the SNP AX-112287343 on chromosome 7B. (a) Manhattan plot for the candidate region on chromosome 7B. The green dot represents the significant SNP AX-112287343. The LD block analysis of the SNPs in this region is shown below. The degree of linkage is represented by the coefficient of r2. (b) Gene structure of TraesCS7B01G036900 and its non-synonymous SNP. Yellow rectangles and black lines represent exons and introns, respectively. The arrow indicates the position of amino acid variations in TraesCS7B01G036900. (c) Box plots for RSA traits based on the allele at SNP AX-112287343. Differences between the haplotypes were statistically analysed using the Student’s t-test. The upper and lower edges of the box represent the 75th and 25th quantiles, respectively, and the whiskers show the 90th and 10th quantiles; the horizontal solid line represents the average; * Significant at p ≤ 0.05, ** Significant at p ≤ 0.01, *** Significant at p ≤ 0.001. (d) Expression profile of TraesCS7B01G036900 in different tissues, as determined by RT-qPCR. The data are the means ± Sd of n = 3. ‘ns’ means ‘not significant’.
Ijms 23 01843 g006
Table 1. Components of eight subpopulations.
Table 1. Components of eight subpopulations.
Sub-PopulationNumber of AccessionExotic VarietyLocal LandraceLocal Cultivar
Sp16611649
Sp2417529
Sp3645158
Sp447131519
Sp54110328
Sp6330033
Sp7301128
Sp86619740
Table 2. Phenotypic analysis of nine RSA traits for this wheat panel under two environments.
Table 2. Phenotypic analysis of nine RSA traits for this wheat panel under two environments.
TraitEMeanSDMinMaxbkbsCV (%)H2 (%)p-Value<0.001
TRLOPC236.0694.9664.25645.141.841.1440.22%77.83%7.96175 × 10−51***
IPC147.9150.0353.91336.760.650.8333.82%
TRAOPC31.6812.989.5086.521.250.9140.97%76.86%2.91719 × 10−81***
IPC16.454.937.1136.000.550.8029.94%
ARDOPC0.430.070.270.58−1.06−0.3216.91%66.78%9.8363 × 10−58***
IPC0.360.030.280.45−0.290.378.52%
TRVOPC0.350.170.071.091.510.9547.94%76.86%4.92408 × 10−89***
IPC0.150.040.070.446.181.6029.81%
NRTOPC360.80186.97113.001703.678.552.2051.82%72.38%1.24864 × 10−28***
IPC233.30110.0878.67699.832.141.3647.18%
RLOPC113.9163.2918.27429.164.951.8355.56%75.71%3.14966 × 10−23***
IPC75.9935.8219.24225.261.060.9847.14%
RAOPC6.903.591.1524.263.781.6252.08%76.15%6.77309 × 10−28***
IPC4.482.111.0111.820.440.8947.13%
RVOPC0.040.020.00670.12143.041.4950.12%76.06%7.24198 × 10−31***
IPC0.020.010.00530.06260.250.8748.38%
RTOPC311.12187.6445.331641.838.302.2060.31%71.15%4.80007 × 10−17***
IPC216.96106.3167.33678.002.221.3649.00%
E, Environments; bk, kurtosis; bs, skewness; CV, coefficient of variation; H2, heritability; *** Significant at p ≤ 0.001.
Table 3. Summary of the genetic diversity in three sub-genomes and chromosomes of this wheat panel and evaluation of the effective number of independent SNPs, including the suggested p-value thresholds.
Table 3. Summary of the genetic diversity in three sub-genomes and chromosomes of this wheat panel and evaluation of the effective number of independent SNPs, including the suggested p-value thresholds.
ChromosomeNo. of MarkersEffective NumberEffective RatioSuggested p-Value% MarkersLength (Mb)Marker DensityHePIC
1A287935643.190.201.77 × 10−47.00594.0248.470.650.26
2A286766088.170.211.64 × 10−46.97780.7636.730.710.28
3A193664607.430.242.17 × 10−44.70750.7325.800.680.28
4A178564241.960.242.36 × 10−44.34744.5423.980.660.27
5A224824935.880.222.03 × 10−45.46709.7631.680.780.31
6A164383851.640.232.60 × 10−43.99617.9726.600.730.29
7A281906297.360.221.59 × 10−46.85736.6938.270.680.27
1B205844750.10.232.11 × 10−45.00689.3829.860.750.30
2B287156701.210.231.49 × 10−46.98801.2535.840.730.29
3B457667920.520.171.26 × 10−411.12830.7055.090.780.31
4B131303016.410.233.32 × 10−43.19673.4719.500.720.29
5B338746596.790.191.52 × 10−48.23713.0247.510.790.31
6B255495554.710.221.80 × 10−46.21720.9535.440.690.28
7B170404491.90.262.23 × 10−44.14750.6122.700.720.29
1D105973058.850.293.27 × 10−42.57495.4421.390.680.28
2D104303731.570.362.68 × 10−42.53651.8116.000.690.28
3D72912668.70.373.75 × 10−41.77615.4811.850.670.27
4D41281702.710.415.87 × 10−41.00509.858.100.660.27
5D84373285.780.393.04 × 10−42.05566.0414.910.650.26
6D74862940.590.393.40 × 10−41.82473.5615.810.640.26
7D101103669.40.362.73 × 10−42.46638.6515.830.640.26
A genome161801 39.314934.4732.790.700.28
B genome184658 44.865719.3832.290.740.30
D genome58479 14.213950.8314.800.660.27
Total40493895754.87 14,064.6828.790.700.28
Average 2.51 × 10−4
Table 4. QTLs and QTL clusters with RSA traits and overlapping with QTLs from previous studies.
Table 4. QTLs and QTL clusters with RSA traits and overlapping with QTLs from previous studies.
TraitsSNPCHRMbp−log10 (p)R2 (%)Reference
ARDAX-947981622B245.898.53 × 10−54.071.23
ARDAX-947578652B797.546.71 × 10−54.17–4.382.24–2.76
ARDAX-1108878442D380.987.38 × 10−54.132.57
ARDAX-1091253752D640.049.95 × 10−54.00–5.381.82–3.41
TRVAX-1088924353A21.287.86 × 10−54.102.58
ARDAX-1114909343B49.299.67 × 10−54.013.59
TRVAX-948979064A102.399.21 × 10−54.042.28
TRVAX-1107126454B606.857.29 × 10−54.14–4.744.77–5.11
TRVAX-950737365A632.605.59 × 10−54.252.64
TRLAX-1115370955B639.689.81 × 10−54.01–4.304.16–4.95
ARDAX-950931235D385.392.69 × 10−54.572.31
ARDAX-949995826A104.366.99 × 10−54.164.23
ARDAX-948273036A115.188.43 × 10−54.072.84
ARDAX-945028646B710.327.22 × 10−65.143.17[27]
ARDAX-948268246D60.145.84 × 10−54.232.36
ARDAX-944680396D86.377.64 × 10−54.123.18
ARDAX-1108377686D464.846.32 × 10−54.202.29[28]
ARDAX-948498697A19.017.99 × 10−54.103.46
RTAX-1088529937A40.267.97 × 10−54.103.76
NRTAX-949709017A692.058.53 × 10−54.074.03[27]
RTAX-951819927B228.815.95 × 10−54.233.84
ARDAX-1095815437B403.543.66 × 10−54.443.03[29]
ARDAX-1103604357D104.099.83 × 10−54.012.24
TRA/TRVAX-1116305482A61.288.82 × 10−54.05–4.292.68–3.20
TRA/TRVAX-1105640362A740.018.78 × 10−54.06–4.113.41–3.69
TRA/TRVAX-1095372442B694.268.70 × 10−54.06–4.292.52–2.64
NRT/RTAX-951936483A0.562.39 × 10−54.62–4.923.97–4.29
NRT/RTAX-947110433B0.991.17 × 10−54.93–5.144.28–4.42
NRT/RTAX-1093371923B2.283.94 × 10−54.40–4.803.77–4.16
TRL/RL/RA/RVAX-1092375103D549.749.90 × 10−54.00–4.163.68–4.58
NRT/RTAX-1104448574A584.719.92 × 10−54.00–4.374.64–4.78
NRT/RT/RL/RA/RVAX-861794465B571.239.98 × 10−54.00–5.733.99–6.77[28]
NRT/RTAX-1105226575D464.078.90 × 10−54.05–4.663.97–5.60
NRT/RTAX-1110608336B127.235.66 × 10−54.25–4.424.00–4.15
TRL/TRAAX-951558297B36.494.72 × 10−54.33–4.584.51–4.80
TRL/TRAAX-947594167B594.418.96 × 10−54.05–4.234.00–5.28
Table 5. Published studies on RSA traits and their chromosome localisations.
Table 5. Published studies on RSA traits and their chromosome localisations.
ChrTraitsReference
1ARGA, TRL, NRT, TRV, RHL, MRL, TRA, ARD[5,22,23,26,29,30,32,34,35,44,45]
1BTRL, NRT, RHL, TRA, ARD, TRV, MRL[5,22,23,26,27,29,32,34,35,36,45]
1DNRT, TRA, TRL, ARD, NRT, TRA[28,30,32,37]
2ATRL, TRV, ARD, TRA, NRT, RGA, MRL[5,23,30,32,33,34,35,36,37,38,44]
2BRGA, TRL, TRV, TRA, ARD, NRT, MRL[5,19,22,23,24,29,30,31,32,33,34,35,37,38,44,45]
2DTRL, RHL, TRV, TRA, ARD, NRT, MRL[22,24,26,31,32,33,37,45]
3ARGA, TRL, NRT, RHL, MRL, TRV, TRA, ARD[5,22,23,24,26,30,32,33,34,35,36,43,45]
3BRGA, TRL, NRT, MRL, ARD, TRV, TRA[5,22,23,24,30,33,34,35,38,39,45]
3DARD, TRL, NRT[32,38,39,45]
4ARGA, TRL, NRT, TRV, TRA, ARD[22,23,25,28,30,32,33,34,35,37,38,42,45]
4BTRL, NRT, ARD, TRA, RGA, TRV[5,22,23,28,30,31,32,34,35,36,37,39,42,43,45]
4DTRL, TRA, NRT, MRL, TRV[31,33,34,38]
5ATRL, ARD, TRV, NRT, TRA, MRL[5,19,22,23,25,30,33,34,35,38,45]
5BRGA, TRL, NRT, TRA, MRL, TRV, TRA, ARD[22,23,28,29,30,31,32,33,34,35,45,56,58]
5DTRV, TRL, TRA, ARD, MRL[22,28,29,32,34,37,45]
6ARGA, NRT, MRL, TRL, ARD, TRA[5,22,23,29,31,34,35,36,44]
6BRGA, TRL, NRT, MRL, RHL, TRA, ARD, TRV[5,19,22,23,24,26,29,34,35,36,39,43]
6DNRT, TRL, TRA, ARD, NRT, MRL, RHL[24,26,28,29,33,34,40,46]
7ARGA, TRL, NRT, TRA, MRL, ARD, TRV[5,22,23,24,29,30,31,32,33,34,35,37,39,40,41,42,45]
7BRGA, NRT, TRL, TRA, MRL, ARD, RHL, TRV[5,19,22,23,24,26,29,30,31,33,34,35,36,39,45]
7DNRT, TRA, TRL, TRV, MRL[22,29,30,32,34,37,38,40,45]
Table 6. Abbreviations for RSA traits and their units.
Table 6. Abbreviations for RSA traits and their units.
AcronymRSA TraitsUnits
TRLTotal root lengthCentimetres (cm)
TRATotal root surface areaSquare centimetres (cm2)
ARDAverage root diameterMillimetres (mm)
TRVTotal root volumeCubic centimetres (cm3)
NRTNumber of root tipsNumber (no.)
RLRoot length (root diameter ≤ 0.3 mm)Centimetres (cm)
RARoot surface area (root diameter ≤ 0.3 mm)Square centimetres (cm2)
RVRoot volume (root diameter ≤ 0.3 mm)Cubic centimetres (cm3)
RTRoot tips (root diameter ≤ 0.3 mm)Number (no.)
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ma, J.; Zhao, D.; Tang, X.; Yuan, M.; Zhang, D.; Xu, M.; Duan, Y.; Ren, H.; Zeng, Q.; Wu, J.; et al. Genome-Wide Association Study on Root System Architecture and Identification of Candidate Genes in Wheat (Triticum aestivum L.). Int. J. Mol. Sci. 2022, 23, 1843. https://doi.org/10.3390/ijms23031843

AMA Style

Ma J, Zhao D, Tang X, Yuan M, Zhang D, Xu M, Duan Y, Ren H, Zeng Q, Wu J, et al. Genome-Wide Association Study on Root System Architecture and Identification of Candidate Genes in Wheat (Triticum aestivum L.). International Journal of Molecular Sciences. 2022; 23(3):1843. https://doi.org/10.3390/ijms23031843

Chicago/Turabian Style

Ma, Jianhui, Dongyang Zhao, Xiaoxiao Tang, Meng Yuan, Daijing Zhang, Mengyuan Xu, Yingze Duan, Haiyue Ren, Qingdong Zeng, Jianhui Wu, and et al. 2022. "Genome-Wide Association Study on Root System Architecture and Identification of Candidate Genes in Wheat (Triticum aestivum L.)" International Journal of Molecular Sciences 23, no. 3: 1843. https://doi.org/10.3390/ijms23031843

APA Style

Ma, J., Zhao, D., Tang, X., Yuan, M., Zhang, D., Xu, M., Duan, Y., Ren, H., Zeng, Q., Wu, J., Han, D., Li, T., & Jiang, L. (2022). Genome-Wide Association Study on Root System Architecture and Identification of Candidate Genes in Wheat (Triticum aestivum L.). International Journal of Molecular Sciences, 23(3), 1843. https://doi.org/10.3390/ijms23031843

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop