Next Article in Journal
Bisphenol A Modulates Autophagy and Exacerbates Chronic Kidney Damage in Mice
Next Article in Special Issue
Role of Reactive Oxygen Species and Hormones in Plant Responses to Temperature Changes
Previous Article in Journal
The OsERF115/AP2EREBP110 Transcription Factor Is Involved in the Multiple Stress Tolerance to Heat and Drought in Rice Plants
Previous Article in Special Issue
Elucidating the Response of Crop Plants towards Individual, Combined and Sequentially Occurring Abiotic Stresses
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genetic Dissection of Seedling Root System Architectural Traits in a Diverse Panel of Hexaploid Wheat through Multi-Locus Genome-Wide Association Mapping for Improving Drought Tolerance

1
Division of Genetics, Indian Council of Agricultural Research (ICAR)—Indian Agricultural Research Institute, New Delhi 110012, India
2
Division of Germplasm Evaluation, ICAR-National Bureau of Plant Genetic Resources, New Delhi 110012, India
3
Division of Genomic Resources, ICAR-National Bureau of Plant Genetic Resources, New Delhi 110012, India
4
ICAR-National Institute of Plant Biotechnology, New Delhi 110012, India
5
Division of Agricultural Economics, ICAR-Indian Agricultural Research Institute, New Delhi 110012, India
6
Department of Agronomy, Kansas State University, Manhattan, KS 66506, USA
*
Author to whom correspondence should be addressed.
These authors contributed equally to the work.
Int. J. Mol. Sci. 2021, 22(13), 7188; https://doi.org/10.3390/ijms22137188
Submission received: 8 June 2021 / Revised: 27 June 2021 / Accepted: 28 June 2021 / Published: 2 July 2021
(This article belongs to the Special Issue Mechanisms of Drought, Temperature and Salinity Tolerance in Plants)

Abstract

:
Cultivars with efficient root systems play a major role in enhancing resource use efficiency, particularly water absorption, and thus in drought tolerance. In this study, a diverse wheat association panel of 136 wheat accessions including mini core subset was genotyped using Axiom 35k Breeders’ Array to identify genomic regions associated with seedling stage root architecture and shoot traits using multi-locus genome-wide association studies (ML-GWAS). The association panel revealed a wide variation of 1.5- to 50-fold and were grouped into six clusters based on 15 traits. Six different ML-GWAS models revealed 456 significant quantitative trait nucleotides (QTNs) for various traits with phenotypic variance in the range of 0.12–38.60%. Of these, 87 QTNs were repeatedly detected by two or more models and were considered reliable genomic regions for the respective traits. Among these QTNs, eleven were associated with average diameter and nine each for second order lateral root number (SOLRN), root volume (RV) and root length density (RLD). A total of eleven genomic regions were pleiotropic and each controlled two or three traits. Some important candidate genes such as Formin homology 1, Ubiquitin-like domain superfamily and ATP-dependent 6-phosphofructokinase were identified from the associated genomic regions. The genomic regions/genes identified in this study could potentially be targeted for improving root traits and drought tolerance in wheat.

1. Introduction

Globally, wheat (Triticum aestivum L.) is one of the most important cultivated food grain crops with an area and production of 218.54 million hectares and 771.71 million tons, respectively [1]. Wheat ranks second after rice (Oryza sativa L.) as a major staple food crop in India with an acreage of 30.55 million hectares, and total production and productivity of 107.18 million tons and 3508 kg/ha, respectively [2]. Wheat is adapted to various climatic conditions including areas vulnerable to heat, salt and drought stress [3,4]. Drought is one of the major challenges that limits crop growth and yield in different areas of the world [5,6]. Drought stress can cause significant yield reduction in wheat and its impact varies with intensity, timing and duration of stress relative to crop growth stages. Yield reduction to the tune of 21.0%, 25.8% and 32.0% were reported in wheat under mild (55–70% relative soil water content), moderate (35–55%, relative soil water content) and severe (<35% relative soil water content) drought stress, respectively [7].
Root traits are important not only for water absorption from a drought stress perspective, but also for soil nutrient uptake and environmental stress tolerance [8,9,10]. Morphology of the root system includes basic features such as root length, root length density, root diameter, surface area and volume which influence the root structure’s spatial arrangement and are significantly correlated with water uptake and nutrient absorption [11,12,13]. Drought tolerance is directly associated with root diameter, as thicker roots with large xylem vessels are more effective towards extraction of water and nutrients from deep soil layers under rainfed condition [12]. Fine roots and the number of root tips, the major components of root systems, often increase the water and nutrient uptake efficiency by increasing the root surface area and volume [8,10]. Further, deep rooting, an essential root feature allowing access to water from deeper soil profiles, improves crop productivity [14]. Therefore, one approach to improve drought tolerance is to identify and introgress deeper rooting alleles in shallow rooted and drought prone cultivars [15,16,17]. Thus, the study of root system architecture (RSA) features can help define proxy traits to enhance tolerance to various soil types and status of moisture and nutrient stress conditions [18]. RSA characteristics are governed by polygenes with additive cumulative gene action and are influenced by the environment [12]. Due to the difficulty of high-throughput phenotyping of RSA characteristics in field conditions, their optimization was ignored [10,12,19]. However, seedling level root architecture is associated with adult plant root architecture and is being studied to translate the knowledge for later growth stages [20,21].
The knowledge on extent of genetic variability in wheat germplasm collections could pave the way for successful conservation and utilization in crop improvement. In recent years, cost effective advanced high throughput genotyping approaches, such as single nucleotide polymorphism (SNP) arrays and sequencing based genotyping methods such as genotyping by sequencing (GBS), restricted site associated DNA sequencing (RAD-seq) and whole genome re-sequencing, have greatly accelerated genomic level characterization of genetic diversity of crop germplasm resources in many species. Availability of these techniques has made it feasible to conduct whole genome scan studies, such as association mapping (AM) and genomic selection in a variety of crops including wheat, maize (Zea mays L.), and rice [22,23]. Association mapping is an alternative approach to mapping quantitative trait loci (QTLs) controlling complex traits which is based on the principle of linkage disequilibrium between marker and trait. In contrast to the conventional QTL approach, the AM approach enables screening of a large number of alleles present in natural populations of species, thus providing greater resolution mapping of a target trait. In the past couple of years, AM has been widely applied to map traits in many crop species including wheat, rice, barley (Hordeum vulgare L.), maize and soybean (Glycine max L. Merril) [21,24,25,26,27] with the availability of high-throughput genotyping and phenotyping platforms [22,23].
Exploring genetic variability of root traits in a diverse germplasm collection could assist in development of varieties with desired root features for drought tolerance or target environments. To date, only a few reports are available on genetic variability and QTL mapping of root traits in wheat [28,29,30]. The wheat core and mini core collections constituted from the large wheat accessions conserved in the national and international gene banks could be a valuable genetic resource for mapping complex traits using a genome wide association study (GWAS). The GWAS could facilitate uncovering genomic regions/genes controlling various root architecture traits. In our earlier study, a core set of wheat genotypes was identified from the entire wheat collection in the national gene bank of India based on 35 agro-morphological traits [31]. This core set was subsequently reduced to a mini core collection that captured huge variability for various traits and could be used for mapping root associated traits. The AM analysis on this mini core collection could facilitate identification of genes/genomic regions associated with RSA traits that could be exploited in a breeding program to develop drought tolerant cultivars. Keeping this in mind, the objectives of this research were: (1) to study genetic variation of seedling root system architecture in a diverse set of wheat genotypes including a subset of wheat mini core collection; and (2) to identify genomic regions/candidate regions linked with these traits using association mapping.

2. Results

2.1. Phenotypic Analyses

There were large differences for various root traits (Table 1) in the wheat association panel owing to their diverse geographical origin (Supplementary Table S1). The abbreviations used here are listed and expanded in Table 1. A significant difference (p ≤ 0.001) was observed among the genotypes for each phenotypic trait except for lateral root density, LRD (Supplementary Table S2). Mean values of the traits (mean ± SE) were as follows: total root size, TRS (242.94 ± 5.35); first order lateral root number, FOLRN (95.22 ± 3.58); second order lateral root number, SOLRN (165.57 ± 8.97); lateral root size, LRS (0.39 ± 0.01); lateral root density, LRD (2.81 ± 0.08); root length density, RLD (0.49 ± 0.01); average diameter, AD (0.38 ± 0.01); root volume, RV (0.16 ± 0.01); seminal root number, SRN (5.06 ± 0.04); shoot dry weight, SDW (35.58 ± 0.88); root dry weight, RDW (10.76 ± 0.32); shoot length, SL (30.61 ± 0.41); root length, RL (24.81 ± 0.54); root shoot dry weight ratio, RSDWR (0.31 ± 0.01) and specific root length, SRL (2.52 ± 0.07) (Table 1). We observed a high range of coefficient variation (%CV) for all the phenotypic traits. The highest CV was observed for SOLRN (64.10%) followed by FOLRN (44.50%) and LRS (38.39%). The lowest CV was observed for AD (8.46%) followed by SRN (10.29%). The estimates of broad-sense heritability were low to moderately high, ranging from 15.16% for LRD, to 78.95% for AD (Supplementary Table S2).
The frequency distribution pattern of traits is presented in Supplementary Figure S1. Most of the traits followed normal distribution, except SRN and SL which were negatively skewed, while SOLRN and SRL were positively skewed. Selected wheat genotypes with high and low mean performance of root and shoot morphological traits are shown in Supplementary Table S3. The wheat genotype IC28755 exhibited maximum root length of 37.50 cm, whereas the genotype EC6903, a collection from the USA, showed a maximum shoot length of 40.76 cm. The genotype IC128151 (C306) had highest value of 60.41 mg for SDW and the genotype EC576889 had maximum value for RDW at 23.06 mg. The trait root-shoot dry weight ratio indicated the amount of root dry mass to shoot dry mass. The genotype IC532019 had the highest ratio of 0.59 contributing higher root dry mass than shoot at seedling stage. Donors for multiple traits were also identified. Genotype IC82425A exhibited the highest mean values for SOLRN, LRS and RV. Further the genotype IC128151 (C306) was among the top ten genotypes for seven traits, namely TRS, FOLRN, LRD, RLD, RV, SDW and RDW. The genotype IC542076 ranked in the top ten genotypes for seven traits: TRS, SOLRN, RLD, RV, SRN, SDW and RDW, and was therefore considered promising for tolerating drought by extracting water stored in the deeper soil layers. The genotypes IC128151, IC82425A and IC542076 were superior for several root and shoot traits in the drought experiment carried out at the seedling stage under pot screening during the winter season 2019–2020(data not presented). In addition, genotypes with low mean performance were also identified so they could be utilized to generate bi-parental mapping populations or multiparent advanced generation intercross (MAGIC) populations. An exotic collection from the USA, EC339632, was observed to have lowest mean performance for TRS, FOLRN, SOLRN, RLD, RV, SRN, RDW and RL, whereas EC187159 showed lower mean values for TRS, SOLRN, LRS, RLD, SRN, SDW, RDW and SL.

2.2. Multivariate Analyses of Phenotypic Traits

2.2.1. Association Analysis

Association analysis revealed the influence of a correlated response from different traits (Figure 1). A highly significant positive association was observed between TRS with FOLRN (r = 0.50), SOLRN (r = 0.65), LRS (r = 0.42), RLD (r = 0.99), RV (r = 0.77), SRN (r = 0.42), SDW (r = 0.52), RDW (r = 0.70) and RL (r = 0.71). TRS showed a negative association with AD (r = −0.37), SRL (−0.27) and LRD (r = −0.01). The traits FOLRN, SOLRN, LRS and RV had a highly significant positive association with all traits except AD and SRL. Further, LRD showed a highly significant positive correlation with FOLRN (r = 0.75), SOLRN (r = 0.36) and LRS (r = 0.36) and a positive significant correlation with RV (r = 0.21), SRN (r = 0.23), SDW (r = 0.20) and RDW (r = 0.21). The RLD had a positive and highly significant association with all traits, while there were negative associations with AD (r = −0.37), LRD (r = −0.01) and SRL (r = −0.27). In contrast, AD was significantly correlated with all the traits in a negative direction except SL. Similarly, SRL was also negatively correlated with all traits.

2.2.2. Principal Component Analysis

Principal component analysis (PCA) was performed for 15 RSA and shoot traits to determine their contribution to overall variability. The first four principal components (PCs) in the PCA analysis with eigenvalues >1 contributed a cumulative variance of 79.24% among the genotypes (Supplementary Table S4). The scree plot showing the eigenvalues on the y-axis and the number of components on the x-axis also depicted the number of relevant PCs to explain maximum variance (Supplementary Figure S2A). Based on the contribution towards the principal components, RV, RDW, SOLRN, RLD, TRS, RL, LRS, FOLRN, SRN and RSDWR were the major contributors toward PC1 and were most significantly associated with the genotypic variations. Likewise, SDW and SL were major contributors to PC2, LRD and FOLRN to PC3 and SRL to PC4 (Supplementary Table S4). Length of coordinates or distance from the origin revealed the amount of variation contributed by each trait (Supplementary Figure S2B) which also ascertained about their major contributions. The PCA also indicated trait association, as the correlated traits have fewer angles between them or are mostly parallel to each other. In this case, the positive quadrant had correlated traits SDW and SRN. The results indicated that AD and SRL were opposite to other traits in the second quadrant, hence their negative association with them.

2.2.3. Clustering of Genotypes

The diverse association panel of the wheat mini core subset (136 accessions, four had a poor SNP call and were therefore excluded from the analysis) was divided into six clusters (Supplementary Figure S3) when subjected to hierarchical analysis based on the 15 root and shoot traits, which showed that the wheat genotypes were sufficiently diverse. Among different classes, the largest cluster (II) contained 38 member genotypes and the smallest cluster (III) had 2 members (EC339632, a USDA collection and EC578153 and a synthetic line). Clusters six and three consisted of genotypes with high and low performance respectively for many traits such as TRS, FOLRN and SOLRN (Supplementary Table S5). A heatmap generated by two-way clustering also depicted that cluster six had high performing genotypes, whereas cluster 1 had low performing genotypes for many of the traits (Supplementary Figure S3).

2.3. Marker Distribution, Linkage Distribution and Population Structure

The SNP genotyping of 136 wheat genotypes yielded a total of 16,616 polymorphic SNPs with exact chromosomal position and minor allele frequency (MAF) ≥ 10. The SNP markers were proportionally distributed across the 21 chromosomes and their number ranged from 287 for Chr4D to 1179 for Chr2B, and the average SNP markers per chromosome was 791. The polymorphic SNPs covered 14066.28 Mb of wheat genome. The sub-genome A included 5355 SNPs, the sub-genome B had 6343 SNPs and 4918 SNPs were located on the sub-genome D. On sub-genome A, maximum SNPs were found on 2A (952) followed by 1A (911). Further, on sub-genome B, maximum SNPs were found on Chr2B (1178), followed by Chr5B (1125). Chr2D of Sub-genome D had maximum SNPs (1037) followed by Chr1D (895) (Table 2). The SNP marker density per chromosome per Mb ranged from 0.56 for 4D to 1.59 for 2D, with a mean of 1.182. These polymorphic 16,616 SNP markers were used to estimate genome-wide linkage disequilibrium (LD) and the pairwise pattern of LD was calculated using r2 of allele pairs between two loci. Among these, a total of 19,528 SNP pairs were in perfect LD (r2 = 1), with B sub-genome having a maximum (7287) number followed by D (7100) and A (5141) sub-genomes. Overall, in the current association panel, low mean LD was observed across the 21 chromosomes with values ranging from 0.087 (4D) to 0.298 (1D). The average LD (r2) per chromosome ranged from 0.087 for 4D to 0.298 for 1D. The measures of LD decay for all 21 chromosomes are summarized in Supplementary Figure S4. The LD decay in the AM panel was determined at r2 = 0.157 (background LD). LD decayed the fastest in sub-genome A followed by sub-genome D and B. In sub-genome A, LD decayed at 3 Mb compared to 3.4 Mb in sub-genome B and 4.2 Mb in sub-genome D (Table 2 and Supplementary Figure S4).
Population analysis of the association panel was performed using both STRUCTURE and PCA which indicated optimal number of population (K) as two (Figure 2). Of the total accessions, 40 accessions belonged to subpopulation 1 and 96 accessions belonged to subpopulation 2. Further, in subpopulation 1 (SP1), 30.50% of the accessions were admixtures and the remaining 62.50% of the genotypes were pure (>80% similarity) which included 25% EC (exotic) collection and 75% IC (indigenous) collection. The IC belonged to multiple states in India (Rajasthan, Gujarat, Haryana, West Bengal, Uttarakhand, Uttar Pradesh, Himachal Pradesh, Madhya Pradesh and Maharashtra). The EC belonged to Israel, China, the USA and Australia. In subpopulation 2 (SP2), 42.70% of the genotypes were admixtures and 57.29% of the genotypes were pure which included 31.25% EC and 68.75% IC collection. The EC mainly belonged to Mexico, Russia, Finland and the USA. Based on PCA analysis, genotypes contributing more variation to the PC1 were parallel to the PC1 axis and genotypes contributing more variation to the PC2 were parallel to the PC2 axis (Figure 2). A total of 88.2% cumulative variance were explained by two PCs. It also indicated that genotypes used in this study were divided into two structured subpopulations.

2.4. Genome-Wide Association Studies

Marker-trait association analyses using six different ML-GWAS models (mrMLM, FASTmrMLM, FASTmrEMMA, pLARmEB, ISIS EM-BLASSO and pKWmEB) revealed 456 significant quantitative trait nucleotides (QTNs) for 15 seedling stage root architecture and shoot traits at a logarithm of the odd (LOD) score ≥3 with phenotypic variance in the range of 0.12–38.60%. Among these, 18, 59, 93, 92 and 194 significant QTNs were detected using FASTmrEMMA, mrMLM, FASTmrMLM, ISIS EM-BLASSO and pLARmEB, respectively (Table 3 and Figure 3). In this case, pKWmEB model could not detect any significant QTN. Maximum number of QTNs were recorded for AD (68) and RLD (68), followed by LRS (46), SOLRN (42) and RV (37). The extent of phenotypic variation of the QTNs varied from R2 = 1.8% for QTN-SOLRN at Chr 6A- AX-95244609 to R2 = 38.6% for QTN-SOLRN at Chr 3B- AX-95150666. The highest number of QTNs was identified in chromosome 7B (44) followed by 6A (34), 3B (33), 7A (31), 1D (31), 2B (29), 1B (29) and 6B (27). Further, among the total 456 QTNs, 87 QTNs and 27 QTNs were repeatedly detected by two, three or more ML-GWAS models, respectively. These were considered reliable QTNs for the respective traits. The reliable QTN that explained phenotypic variation ≥ 10% for the trait in at least one of the GWAS models was considered a major QTN. Among these 87 reliable QTNs, 11 QTNs were found to be associated with AD and 9 each for SOLRN, RV and RL, respectively. Further, a total of 12 QTNs for traits including SOLRN: AX-95077960, AX-94448890, AX-95150666 and AX-95244609; LRS: AX-95123855; AD: AX-95244609; RV: AX-94439232 and AX-94664277; SRN: AX-94551988; SDW: AX-94470023, RL: AX-94502864 and SRL: AX-94457792 were detected in four or five models (Table 3). Manhattan plots depicting identification of a reliable QTN for SDW (Q.SDW-4AL) and RL (Q.RL-6BL) by five ML-GWAS methods are presented in Supplementary Figures S5 and S6. Further, graphical representation of the genetic position of identified reliable QTNs on a physical map was given in Figure 4. It was revealed that the top five genomic regions that contained a maximum number of significant QTNs were found on chromosomes Chr1B, Chr3B, Chr4A, Chr6A and Chr7B.

2.4.1. Colocalization of Root Architecture Loci

Different root traits are mostly correlated and complex biological mechanisms might be involved in coordination for their expression. The pleiotropic action of genetic loci on different traits or their tight linkage results in a correlation between traits. A total of 11 QTNs were found to be pleiotropic markers (Table 4). A locus on 1D (AX-94448890) was associated with root morphology (SOLRN and RL) and another locus on 2A (AX-94952472) was associated with TRS and RLD. Similarly, other loci were also identified as pleiotropic on 2B (AX-95227366) for RV and RL, 3B (AX-94516395) for SOLRN and RLD, 4A (AX-95105488) for TRS, RLD and RV, 6A (AX-95244609) for SOLRN, LRS, AD and 6B (X-94502864) for RLD, RV and RL. We also identified two pleiotropic loci on the short arm of 7B (AX-95123855) associated with root morphology traits (FOLRN, SOLRN, LRS) and QTN-7B-AX-95119337 associated with TRS and RLD whereas its long arm had a locus (AX-94528392) associated with RL and LRS. Other pleiotropic QTNs were 1DS- AX-94448890 (SOLRN and RL) and 7AL- AX-94664277 (RV and SDW) (Figure 4 and Table 4). These identified loci influencing several traits could be potential markers for marker-assisted breeding after validation.

2.4.2. Allelic Effects of Identified Genomic Regions on Respective Phenotypes

Seven major loci detected in four or more ML-GWAS methods were analyzed to determine the range of phenotypic variations for all traits (Figure 5). Association panel genotypes were divided into two classes according to allele types. It was observed that all seven QTNs demonstrated a significant effect on respective traits (p ≤ 0.01). Among these seven QTNs, three showed a significant effect on SOLRN (Q.SOLRN-1DS: AX-94448890, Q.SOLRN-1BL: AX-95077960, Q.SOLRN-6AL: AX-95244609) and one each on LRS (Q.LRS-7BS: AX-95123855), RV (Q.RV-5BL: AX-94439232), SDW (Q.SDW-4AL: AX-94470023) and AD (Q.AD-6AL: AX-95244609). The marker, AX-95244609 at Chr6AL was pleiotropic for SOLRN, LRS and AD. These significant associations suggest their plausible role in determining root architecture traits.

2.4.3. Genes Linked to Quantitative Trait Nucleotides

Out of the total SNPs associated with various root and shoot traits, 52 SNP markers were found to be located within genes (Supplementary Table S6). Some important proteins encoded by these genes included glycoside hydrolase superfamily protein, carbonic anhydrase superfamily protein, auxin response factor, nitrogen regulatory protein, transmembrane protein, enhancer of polycomb-like protein and glutamate receptor. Moreover, we found that some of the loci corresponding to large effect QTNs and detected by multiple methods harbored genes for the growth and development of wheat root and shoot organs. For instance, QTN locus Chr2B: AX-94457792 for SRL is comprised of genes TraesCS2B02G406800 and is predicted to encode Formin-like protein which regulates root-hair elongation in rice [32]. This locus also controls the actin cytoskeleton in root hair in wheat [33]. The other largely effected QTN associated with SL (Chr7BS: AX-94876335) was linked to a transcriptional responses of maize seedling root to phosphorus starvation [34]. Most QTNs associated with SRN were at loci containing genes predicted to encode Histone-lysine N-methyltransferase ATXR2, ATPase, F1/V1/A1 complex, alpha/beta subunit and the N-terminal domain superfamily. Other harboring genes encode nucleotide-diphospho-sugar transferases, ribosome-inactivating protein superfamily and papain-like cysteine peptidase superfamily. The large effect QTNs Chr7AS: AX-95249973 and Chr7BS: AX-95123855 associated with FOLRN were involved in plant growth, development and adaptation (Supplementary Table S6).

3. Discussion

Root traits play a major role in resource uptake, particularly with respect to water access and uptake in the context of drought tolerance. In this way, root traits help maintain crop yield under limited water environments. To realize maximum yield in wheat, a well-developed root system is warranted [28,35]. The discovery of a novel source of germplasm and new alleles to improve RSA and introgression of new traits into adapted but otherwise susceptible phenotypes are a desirable approach toward breeding for drought tolerance. Different strategies were adopted for early screening of the wheat RSA, assuming that genotypes with a diverse nature of root architecture at the seedling stage would also respond in a similar way at the adult stage when water and/or nutrients become scarce for grain yield [36]. We studied the root and shoot behaviors at the seedling stage under controlled conditions in perlite vermiculite mixtures in pots. This method of screening was considered a reliable method to examine the root system variations compared to field conditions, where there are several confounding impacts and extracting and measuring roots is difficult. Genetic variability in root systems during seed germination and seedling development is a key trait for seedling establishment and early vigor under water stress and resource poor environments. Improved root traits during this key stage will impact crop emergence and seedling establishment resulting in a good plant population which is one of the key components of yield under both optimal conditions and stress conditions.

3.1. Phenotypic Variability

The significant variation observed for various RSA features and shoot traits could be attributed to diversity among the genotypes of the association panel due to their diverse genetic background and wide geographical distributions. The association panel genotypes were collected from various parts of India and around the world. Further, a broad range of variation was detected for individual RSA traits. The highest and lowest morphometric values for various RSA traits (Table 1) reflected the level of diversity within the gene pool and could be promising material to improve root traits through breeding. High coefficient of variation (%) for FOLRN, SOLRN, LRS, LRD, RV, RDW and SRL (>30%) revealed wide variability of these traits in the AM panel. The broad sense heritability (H2) of the trait is a critical parameter determining their utility in breeding. We observed high heritability for traits such as SOLRN, SDW, RDW, SL, RL, AD, RV and moderate heritability for traits such as, LRS, RLD, SRN, RSDWR, SRL, TRS (Supplementary Table S2). Li et al. [37] reported moderate to high H2 for root and shoot traits, ranging from 56.0% to 94.6%.
For acquisition of water and nutrients from the deep zone of soil, RV with a greater number of root tips is essentially required [8,10]. Further, it would be desirable to have genotypes with longer RV, RDW, and other associated parameters so they could extract water and maintain growth under moisture stress environments. The genotype IC406521 from Uttarakhand (India) showed the highest mean value of RV (0.35 cm3), and lowest AD (0.31 mm), thus its root distribution occupied a larger area of soil. Further, the genotype IC539574 ranked one with the value of 240.66 for FOLRN whereas the genotype IC82425A recorded the highest mean value of 490.33 for SOLRN. Similarly, top ranking genotypes for RLD, LRS, LRD and AD were IC542076 (0.75 cm−2), IC82425A (0.86), IC29008 (5.89 cm−1) and EC187159 (0.49 mm), respectively. These genotypes could serve as donors of these traits for crop improvement to combat drought stress. The genotypes with the highest AD are considered best for drought stress tolerance due to large xylem vessels which can extract more water and nutrients from deep soil layers [12,38]. We recorded variation for SRN in our diverse panel, ranging from 2 to 6 which was also observed by Djanaguiraman et al. [14]. However, Cane et al. [36] reported 4 to 6 seminal roots in elite wheat varieties which might be due to less variability. We identified donors for multiple traits. For example, the genotype EC426644, an Australian cultivar (Tincurrin), showed high performance for traits, namely TRS, RLD, FOLRN and RSDWR. Bustos et al. [39] also reported Tincurrin cultivar with high root system size in their study.

3.2. Multivariate Analyses

Correlations between root and shoot traits shows the balance between the organs of roots and shoots and resource partitioning between the above ground and below ground plant parts [40,41]. In this case, FOLRN and SOLRN exhibited strong correlations with all the traits except AD and SRL. FOLRN directly or indirectly contributes to the variability in root morphological traits which increase water use efficiency at critical stages of plant growth through more absorption of water and nutrients from soil sub layers. The inverse relationship between AD with other traits observed in the present study was in line with the previous studies on flax (Linum usitatissimum L.) and Arabidopsis thaliana [42,43]. This is because fine roots are the main root feature, and the number of root tips increases the root surface area and volume, hence enhancing the absorption efficiency of water and nutrients [8,10]. A high correlation between RV and RDW is obvious. Moreover, a high correlation between shoot and root dry weights observed might be due to supply of nutrients from root to the shoot parts, as is evident in the case of rice [44]. One of the easily scorable key traits, SRN, showed a high correlation with a majority of the root traits, suggesting its utility to provide a broad idea about the root system. The high correlation between root and shoot traits is due to reliance of crops on the root system for water and nutrient uptake.
Based on PCA, RV, RDW, SOLRN, RLD, TRS, RL, LRS, FOLRN, SRN and RSDWR were major contributors to PC1 and most significantly associated with genotypic variations. Likewise, SDW and SL are major contributors to PC2. Hence these traits need to be focused on root studies. RV indicates the proliferation of roots which is essential for plants to uptake more water and nutrients from the soil, whereas linear root elongation may be quite useful in extracting water from deeper soil layers in the event of drought. However, RDW varies among genotypes with varying RL. Therefore, it would be desirable to have genotypes with high RL as well as RDW for a more efficient root system. Higher root density is found to be useful for the plants to uptake more nutrients and therefore the plants with a deeper root system would be more advantageous for drought tolerant genotypes [45,46]. Clustering of genotypes in six groups was also evident in a two-dimensional PCA biplot based on root and shoot morphological traits. Genotypes from two contrasting clusters 1 and 6 could be used in crossing program for genetics studies, mapping population development and trait introgression.

3.3. Genome Wide Association and Candidate Genes Identification

Identification of a novel source of genes and genomic regions associated with RSA traits at the seedling stage is expected to boost development of high yielding drought tolerant varieties. In this context, historical wheat germplasm collections maintained in gene banks could prove valuable genetic resource for searching genes/genomic regions using association genetics approach. In this study, we used a diverse panel of 136 wheat genotypes including landraces, varieties, local collection, synthetic germplasm and elite breeding line for mapping RSA and shoot traits. For conducting GWAS, ML-GWAS models were used as these are considered superior to SL-GWAS models for mapping complex traits due to the fact that in ML-GWAS models, all-marker effects are simultaneously estimated. Moreover, unlike SL-GWAS models, these do not require testing of identified associations using stringent multiple testing corrections that generally result in rejection of significant associations [47]. Using six different ML-GWAS models, 456 QTNs were detected for 15 RSA features and shoot traits, distributed over 21 chromosomes. The phenotypic variation estimated (PVE) (r2) for the traits ranged from 0.12% to 38.67%, indicating that RSA traits are controlled by multiple loci with small to moderate effects. This also revealed the complex genetic control of these traits at an early stage of crop growth. Among the six ML-GWAS models used in our analysis, the pLARmEB model was the most powerful, which revealed the maximum number of associations (194) whereas the FASTmrEMMA was the least powerful as it detected the lowest number of MTAs (18). These findings were consistent with the observation of the Safdar et al. [48] who used ML-GWAS models to dissect agronomic traits in bread wheat.
Among the root traits, TRS was found to be associated with drought tolerance in wheat due to spreading of roots in the soil and its effect on the resource uptake [49]. The SNP marker AX-94952472, AX-95105488 and AX-95119337 on chr2AS, 4AL and 7BS were associated with TRS as well as RLD with phenotypic variance in the range of 3.29–11.23%. Root length is an important parameter that determines the ability of a plant to capture water from deeper soil layer [50]. Eight QTNs for RL were detected with more than two models, out of which Q.RL-6BL was consistently detected in five models. Further, this marker was pleiotropic and was associated with RV and RLD which were positively associated traits. Two reliable QTNs for RL (Q.RL-6BL and Q.RL-7AL) were not documented in earlier studies. There, they are considered novel QTNs, whereas other QTNs were reported to be associated with wheat seedling root development under abiotic stresses and phosphorus starvation [51,52]. For traits, FOLRN and SOLRN, two and nine QTNs were identified, respectively. Lateral roots play an important role for foraging water from a shallow depth. Wheat cultivation in the dry area requires more lateral roots at deeper layers for absorption of nutrients and water from deep soil layers. Our analysis revealed a pleiotropic marker, AX-95123855 on 7BS associated with RSA traits FOLRN, SOLRN and LRS. This marker was located within a gene (TraesCS7B02G086500) that encodes Calcineurin which was reported to control inward ion flux in the root, a process essential for plant survival and growth [53].
The SRL is a widely used morphological parameter for root trait which characterizes economic aspects of root system, demonstrating the role of root mass for nutrient acquisition. Genotypes with large SRL have thin roots that increase the surface area per unit root volume enhancing water and nutrient uptake efficiencies [54]. Furthermore, genotypes with small root diameter and large RLD at depth are better adapted to drought conditions [55]. Our study identified 2 QTNs for SRL, one each on 5AL and 2BL. Of these, QTN on 2BL (Q.SRL-2BL) was associated with a gene encoding for Formin-like protein, a family of actin organizing protein involved in root hair development [33]. Further, RLD is considered a primary driver of drought avoidance and enables complete uptake of the soil moisture. Breeding for plants with increased RLD in medium and deep layers and less RLD in shallow soil layers have been proposed as an efficient growth strategy where deep water could be available to crops during late maturity [45,56,57,58]. Nine QTNs were identified for RLD in our study. Root diameter is another critical parameter determining nutrient uptake and transport. In this study, we identified as many as eleven consistent QTNs for AD and all of them had minor effects, suggesting complex genetic regulation of this trait. Among these, two QTNs Q.AD-1BL (AX-94620468) and Q.AD-6AL (AX-95244609) were detected in three or more models and could represent important genomic regions controlling this trait. Moreover, we found that marker AX-95244609 for QTN Q.AD-6AL was associated with other traits such as SOLRN and LRS, and thus could be a potential target for improvement of more than one RSA traits.
Another important trait is SRN, which emerges first from coleorhiza of seed embryo. The gravitropic response of wheat seminal roots was heritable and suggested to be under control of a single dominant gene [59]. Soriano et al. [60] reported meta-QTNs for SRN in wheat on 2A, 2B, 3B, 4A, 5A, 7A and 7B. In this study, a total of six novel QTNs for SRN were identified with the QTNs on 2D (Q.SRN-2DL: AX-94551988) having consistency in four methods. We found that SNP marker AX-94551988 associated with Q.SRN-2DL was located within gene TraesCS2D02G478700 encoding for histone-lysine N-methyltransferase. In a previous study on Arabidopsis, histone-lysine N-methyltransferase gene was reported to inhibit lateral root development [61]. Based on this, it could be suggested that identified Histone-lysine N-methyltransferase candidate gene might participate in regulating SRN through modification of histone proteins.
Moreover, root biomass accumulation at the seedling stage is beneficial for drought escape as it could enable efficient uptake of water and nutrient from soil. In previous studies, RDW QTLs have been identified under both normal and drought stress conditions [60]. Among the four QTNs for RDW identified in this study, the one on 7AS (Q.RDW-7AS) explained maximum percentage of phenotypic variance (18.35–18.60%) and coincided with a meta-QTL, root_MQTL_68 [60], suggesting this is an important genomic region associated with root trait and its role can be further validated. Among the studied traits, RSDWR is considered the most important as it indicates ability of a genotype to absorb water and nutrients from the soil. Generally, plants have higher root shoot ratio under nutrient deficit soils, suggesting more allocation of resources for the development of root so the plant can extract nutrients from deep and wider zone [62]. Studies have identified genomic regions for the RSDWR trait on various chromosomes of wheat under both control as well as nutrient deficit conditions [60]. Our analysis revealed three QTNs for RSDWR, one each on 2A, 2B and 3B respectively, however each of them have small effect, indicating limited potential for their utilization in breeding programs.
The underground root traits significantly affect development of above ground traits including shoot traits such as, SDW and SL, and thus grain yield. We could identify six and five reliable QTNs for SDW and SL, respectively. The strongest QTN for SDW was identified on 4A, Q.SDW-4AL which explained phenotypic variance in the range of 8.09–19.05%, and was found reliable, as detected by using all five models. The SNP marker associated with this QTN was located within a gene TraesCS4B02G095100 that encodes for an F-box-like domain superfamily. F-box family proteins play a diverse role in plant growth and development and could be critical for shoot growth as their expression has also been detected in the shoot tissues of plant species [49].
Applications of root phenotype-genotype association through GWAS has enabled the identification of important QTNs for root traits that impact shoot traits including yield [63,64]. Markers significantly associated with these traits can be used in marker assisted backcross breeding for varietal improvement. Pleiotropic loci with consistent effects should be amenable to MAS for many traits together. Further, our study provides understanding into phenotype–genotype associations for early root and shoot traits of diverse wheat genotypes by identifying QTNs and proposing plausible candidate genes for future investigations.

4. Materials and Methods

4.1. Experimental Material and Design

Phenotyping of 140 wheat germplasm of diverse panel including mini core subset was carried out at the Indian Council of Agricultural Research—National Institute of Plant Biotechnology (NIPB), New Delhi (India) as pot screening in growth chambers under controlled environmental conditions during Kharif season (June–September), 2019. Detailed information about the origin and pedigree of the material is given in Supplementary Table S1. Root morphology and root system architecture (RSA) along with shoot traits were studied by destructive methods at early seedling stage (Fifteen-day-old seedlings). Figure 6 describes the steps and methodology used in this experiment.
In this experiment, five germinated wheat seeds were sown initially in pots (4-inch diameter) containing perlite and vermiculite (1:2 ratio v/v). Germinated seeds that were two days old with primary roots of 1 cm in length were placed inside the pot with the root radical facing down at equal depth. The Murashige–skoog (MS) liquid media were supplemented with 8 mM NO3 as nutrient media at uniform intervals during the entire growth period. The calcium concentration in all media was kept constant by adjusting the amount of CaNO3. Murashige–skoog (MS) N+ nutrient solution was prepared as per the given composition and stored in 4 °C. For one liter of working solution preparation, 100 mL of MS media, 4 mL of CaNO3 and the remaining 896 mL of distilled water constituted 1 L of solution, which were added in equal proportion to all the genotypes grown in controlled conditions at equal intervals. In the growth chamber, 150–200 μmol photon/m2/s light intensity, 10/14 dark/light hours, 70% relative humidity (RH) and at 22 ± 1 °C conditions were maintained as described by Sinha et al. [65].

4.2. Root and Shoot Traits Measured for Phenotyping

After completing 15 days of growth in the controlled growth conditions, the seedlings were harvested from the pots by removing the perlite + vermiculite mixture and plants were separated carefully without damaging the roots and shoots. The plants were immediately placed on a tray to wash out the adhering perlite + vermiculite particles by running tap water with the utmost care so that to avoid breakage of any roots. The entire seedling was carefully spread on blotting paper and maximum root length and shoot length were measured with the help of a length measuring scale. The root was then extracted intact by cutting at the collar region using a sharp blade and the roots were placed in a tray containing distilled water. The roots were then individually scanned in an Epson Perfection V 700 Photo® flatbed scanner at a resolution of 400 dpi modified for this purpose (Regent Instruments Inc., Quebec, QC, Canada) as per the manufacturer’s guidelines. The root images from the scanner were analyzed with customized software WinRHIZO™ (Regent Instruments Inc., Quebec, QC, Canada) [66]. Various root trait data were recorded by software which were later transformed in major RSA traits measured manually using published protocols [32,67,68,69] (Table 1).

4.3. Statistical Analyses of Phenotypic Data

Mean data across the replications were used as the input data for statistical analyses and GWAS analysis. Descriptive statistics and frequency distribution were analyzed to check range of variability among the traits. Pearson’s correlation coefficient, cluster analysis and PCA was performed using SAS software version 9.3 (JMP) program (SAS Institute, Cary, NC, USA). Heritability (H2) was estimated from the analysis of variance.

4.4. DNA Extraction and SNP Genotyping

Genomic DNA was extracted from 15-day-old wheat seedlings using CTAB methods [70]. DNA quality was checked using NanodropTM 2000 (Thermo Fisher Scientific, Wilmington, DE, USA). Samples with good quality DNA were genotyped using Axiom® Wheat Breeders’ Array (Thermo Fisher Scientific, Wilmington, DE, USA) according to the procedure described by Affymetrix (Axiom® 2.0 Assay for 384 samples P/N 703,154 Rev. 2). SNP markers with >10% missing data and <10% MAF were excluded. Identified SNPs were localized on a wheat genome assembly International Wheat Genome Sequencing Consortium (IWGSC) RefSeq version 1.0 using BLASTn program with default parameters.

4.5. Population Structure and LD

Population structure was estimated with 525 unlinked SNP markers nearly uniformly distributed across the wheat genome using the Bayesian model-based approach implemented in STRUCTURE program version 2.2 (Pritchard Lab, Stanford University, Stanford, CA, USA) [71]. A burn-in of 20,000 iterations followed by 50,000 Monte Carlo Markov Chain (MCMC) was run to estimate the number of subpopulation (k) in a putative range of k = 1 to 10. The subpopulation number was estimated using an ad hoc statistic delta k based on the rate of change in log probability of data between successive values [72]. The squared allele frequency correlation (r2) between SNP markers was used to estimate LD across sub-genome A, B, and D using TASSEL v5.0 (Buckler’s Lab, Ithaca, NY, USA) [73]. LD decay across the three sub-genomes and whole genome level was estimated as the physical distance between SNPs where average r2 reduced to half of the maximum LD value.

4.6. Genome Wide Association Analysis

The use of multi-locus methods that capture small effect loci in complex polygenic traits such as in plant roots and shoots have recently become a popular and feasible approach. To benefit the algorithmic merits of different models and support results of one by another, it is also advantageous to apply multiple methods [49]. GWAS analysis was performed using six multi locus GWAS methods within mrMLM [74], FASTmrMLM [75], FASTmrEMMA [76], pKWmEB [77], pLARmEB [78] and ISIS EMBLASSO [79], which were included in the R package mrMLM v3.1 [80]. All parameters were set at default values in this GWAS. The critical thresholds of significant association for the six methods were set as LOD score 3.00 or >3.00. The most significant QTNs, detected in at least two methods, were considered as reliable QTNs. The associated SNPs and their putative underlying genes were illustrated on the wheat chromosomes using Map Chart 2.3, (https://www.wur.nl/en/show/Mapchart.htm accessed on 8 June 2021) [81]. The favorability of alleles at QTNs detected by at least four of the multi-locus models was illustrated using a box plot based on mean phenotypic value of genotypes with each allele.

4.7. Identification of Potential Candidate Genes

SNPs (probe sequences) that were significantly associated with root architecture traits were searched against the Triticum aestivum genome assembly IWGSC-refseq version1.0 in online web resource Ensemble plants (https://plants.ensembl.org/Triticum_aestivum/Tools/Blast, accessed on 14 November 2020) using BLASTn with default parameters to identify potential candidate genes. BLAST2GO tool was used to get annotation of expressed transcripts [82].

5. Conclusions

This study revealed wide variability for RSA and shoot traits at the seedling stage in the studied association panel. Ten top and bottom performing lines for 15 traits were identified for use in genetics, development of mapping population and introgression studies. Based on PCA, several traits such as RV, RDW, SOLRN, RLD, TRS, RL, LRS, FOLRN, SRN and RSDWR, were the most influential traits for phenotypic variations. The GWAS analysis enabled genetic detection of RSA traits and revealed 11 pleotropic loci associated with correlated traits. Further, putative candidate genes were identified from the associated genomic region that could be validated using a functional genomics approach. Development of wheat cultivars possessing superior root traits will play an important role in enhancing drought tolerance under water stress conditions. The large variability found in Indian germplasm for RSA traits and the novel genomic regions regulating them makes this germplasm a valuable source for improving root architecture, which plays a significant role in absorption and uptake of water and nutrients and increases crop productivity.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/ijms22137188/s1.

Author Contributions

J.K., A.K.S. and S.K.S. conceived the idea and provided critical inputs to the concept and the experiment. M.C.Y. and A.K.S. generated the genotyping data. T.D., J.K., S.K.S., S.S., R.B. generated phenotyping data. G.K.J., A.K.P., J.K., analyzed the data and interpret the result, J.K., A.K.S., T.D., S.K.J. wrote the manuscript. S.K. and P.V.V.P. edited the final version of manuscript. All authors contributed to the final reading and approved the submitted version. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded jointly by Indian Council of Agricultural Research (ICAR)-National Innovations in Climate Resilient Agriculture (NICRA)—National Bureau of Plant Genetic Resources (NBPGR) project (Grant number: 1006607) and Post Graduate School, ICAR-Indian Agricultural Research Institute.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data available upon request from corresponding author.

Acknowledgments

The authors are thankful to the director, ICAR-National Institute of Plant Biotechnology, New Delhi and ICAR-NBPGR, New Delhi for extending infrastructure facility. J.K. appreciates the opportunity to visit the Department of Agronomy at Kansas State University through support from ICAR. Contribution number 21-336-J from Kansas Agricultural Experiment Station.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. FAOSTAT. 2017. Available online: www.fao.org/faostat (accessed on 27 September 2020).
  2. ICAR-IIWBR. Director’s Report of AICRP on Wheat and Barley 2019-2020; Singh, G.P., Ed.; ICAR-Indian Institute of Wheat and Barley Research: Karnal, Haryana, India, 2020; p. 76. [Google Scholar]
  3. Tardif, G.; Kane, N.A.; Adam, H.; Labrie, L.; Major, G.; Gulick, P.; Sarhan, F.; Laliberte, J.F. Interaction network of proteins associated with abiotic stress response and development in wheat. Plant Mol. Biol. 2007, 63, 703–718. [Google Scholar] [CrossRef] [PubMed]
  4. Monneveux, P.; Jing, R.; Misra, S. Phenotyping for drought adaptation in wheat using physiological traits. Front. Physiol. 2012, 3, 429. [Google Scholar] [CrossRef] [Green Version]
  5. Fleury, D.; Jefferies, S.; Kuchel, H.; Langridge, P. Genetic and genomic tools to improve drought tolerance in wheat. J. Exp. Bot. 2010, 61, 3211–3222. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Zhao, T.; Dai, A. The Magnitude and causes of global drought changes in the twenty-first century under a low–moderate emissions scenario. J. Clim. 2015, 28, 4490–4512. [Google Scholar] [CrossRef]
  7. Zhang, J.; Zhang, S.; Cheng, M.; Jiang, H.; Zhang, X.; Peng, C.; Lu, X.; Zhang, M.; Jin, J. Effect of drought on agronomic traits of rice and wheat: A meta-analysis. Int. J. Environ. Res. 2018, 15, 839. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Osmont, K.S.; Sibout, R.; Hardtke, C.S. Hidden branches: Developments in root system architecture. Annu. Rev. Plant Biol. 2007, 58, 93–113. [Google Scholar] [CrossRef]
  9. Smith, S.; De Smet, I. Root system architecture: Insights from Arabidopsis and cereal crops. Phil. Trans. R. Soc. B 2012, 367, 1441–1452. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Bishopp, A.; Lynch, J.P. The hidden half of crop yields. Nat. Plants 2015, 1, 15117. [Google Scholar] [CrossRef]
  11. 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]
  12. Maccaferri, M.; El-Feki, W.; Nazemi, G.; Salvi, S.; Cane, M.A.; Colalongo, M.C.; Tuberosa, R. Prioritizing quantitative trait loci for root system architecture in tetraploid wheat. J. Exp. Bot. 2016, 67, 1161–1178. [Google Scholar] [CrossRef]
  13. Djanaguiraman, M.; Prasad, P.V.V.; Kumari, J.; Rengel, Z. Root length and root lipid composition contribute to drought tolerance of winter and spring wheat. Plant Soil 2019, 439, 57–73. [Google Scholar] [CrossRef] [Green Version]
  14. Djanaguiraman, M.; Prasad, P.V.V.; Kumari, J.; Sehgal, S.K.; Friebe, B.; Djalovic, I.; Chen, Y.; Siddique, K.H.M.; Gill, B.S. Alien chromosome segment from Aegilops speltoides and Dasypyrum villosum increases drought tolerance in wheat via profuse and deep root system. BMC Plant Biol. 2019, 19, 242. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Steele, K.A.; Virk, D.S.; Kumar, R.; Prasad, S.C.; Witcombe, J.R. Field evaluation of upland rice lines selected for QTLs controlling root traits. Field Crops Res. 2007, 101, 180–186. [Google Scholar] [CrossRef]
  16. Ehdaie, B.; Merhaut, D.J.; Ahmadian, S.; Hoops, A.C.; Khuong, T.; Layne, A.P.; Waines, J.G. Root system size influences water-nutrient uptake and nitrate leaching potential in wheat. J. Agron. Crop Sci. 2010, 196, 455–466. [Google Scholar] [CrossRef]
  17. Uga, Y.; Sugimoto, K.; Ogawa, S.; Rane, J.; Ishitani, M.; Hara, N.; Inoue, H. Control of root system architecture by DEEPER ROOTING 1 increases rice yield under drought conditions. Nat. Genet. 2013, 45, 1097. [Google Scholar] [CrossRef]
  18. Yu, G.R.; Zhuang, J.; Nakayama, K.; Jin, Y. Root water uptake and profile soil water as affected by vertical root distribution. Plant Ecol. 2007, 189, 15–30. [Google Scholar] [CrossRef]
  19. Mace, E.S.; Singh, V.; Van Oosterom, E.J.; Hammer, G.L.; Hunt, C.H.; Jordan, D.R. QTL for nodal root angle in sorghum (Sorghum bicolor L. Moench) co-locate with QTL for traits associated with drought adaptation. Theor. Appl. Genet. 2012, 124, 97–109. [Google Scholar] [CrossRef] [Green Version]
  20. Sertse, D.; You, F.M.; Ravichandran, S.; Cloutier, S. The complex genetic architecture of early root and shoot traits in flax revealed by genome-wide association analyses. Front. Plant Sci. 2019, 10, 1483. [Google Scholar] [CrossRef] [Green Version]
  21. Jia, Z.; Liu, Y.; Gruber, B.D.; Neumann, K.; Kilian, B.; Graner, A.; von Wirén, N. Genetic dissection of root system architectural traits in spring barley. Front. Plant Sci. 2019, 10, 400. [Google Scholar] [CrossRef]
  22. Collard, B.C.; Mackill, D.J. Marker-assisted selection: An approach for precision plant breeding in the twenty-first century. Philos. Trans. R. Soc. B Biol. Sci. 2008, 363, 557–572. [Google Scholar] [CrossRef] [Green Version]
  23. Zhu, X.; Tang, H.; Risch, N. Admixture mapping and the role of population structure for localizing disease genes. In Genetic Dissection of Complex Traits, 2nd ed.; Rao, D.C., Gu, C.C., Eds.; Academic Press: San Diego, CA, USA, 2008; pp. 547–569. [Google Scholar]
  24. Hu, J.; Guo, C.; Wang, B.; Ye, J.; Liu, M.; Wu, Z.; Liu, K. Genetic properties of a nested association mapping population constructed with semi-winter and spring oilseed rapes. Front. Plant Sci. 2018, 9, 1740. [Google Scholar] [CrossRef] [Green Version]
  25. Bollinedi, H.; Yadav, A.K.; Vinod, K.K.; Krishnan, S.G.; Bhowmick, P.K.; Nagarajan, M.; Singh, A.K. Genome-wide association study reveals novel marker-trait associations (MTAs) governing the localization of Fe and Zn in the rice grain. Front. Plant Sci. 2020, 11, 213. [Google Scholar]
  26. Chaurasia, S.; Singh, A.K.; Songachan, L.S.; Sharma, A.D.; Bhardwaj, R.; Singh, K. Multi-locus genome-wide association studies reveal novel genomic regions associated with vegetative stage salt tolerance in bread wheat (Triticum aestivum L.). Genomics 2020, 112, 4608–4621. [Google Scholar] [CrossRef] [PubMed]
  27. Kumar, M.; Rani, K.; Ajay, B.C.; Patel, M.S.; Mungra, K.D.; Patel, M.P. Multivariate diversity analysis for grain micronutrients concentration, yield and agro-morphological traits in pearl millet (Pennisetum glaucum (L) R. Br.). Int. J. Curr. Microbiol. Appl. Sci. 2020, 9, 2209–2226. [Google Scholar] [CrossRef]
  28. Narayanan, S.; Mohan, A.; Gill, K.S.; Prasad, P.V.V. Variability of root traits in spring wheat germplasm. PLoS ONE 2014, 9, e100317. [Google Scholar] [CrossRef] [PubMed]
  29. Rosello, 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]
  30. Chen, Y.; Palta, J.; Prasad, P.V.V.; Siddique, K.H. Phenotypic variability in bread wheat root systems at the early vegetative stage. BMC Plant Biol. 2020, 20, 185. [Google Scholar] [CrossRef] [PubMed]
  31. Phogat, B.; Kumar, S.; Kumari, J.; Kumar, N.; Pandey, A.C.; Singh, T.; Tyagi, R.; Jacob, S.R.; Singh, A.K.; Srinivasan, K.; et al. Characterization of wheat germplasm conserved in the Indian National Genebank and establishment of a composite core collection. Crop Sci. 2020, 61, 604–620. [Google Scholar] [CrossRef]
  32. Huang, J.; Kim, C.M.; Xuan, Y.H.; Liu, J.; Kim, T.H.; Kim, B.K.; Han, C.D. Formin homology 1 (OsFH1) regulates root-hair elongation in rice (Oryza sativa). Planta 2013, 237, 1227–1239. [Google Scholar] [CrossRef]
  33. Pei, W.; Du, F.; Zhang, Y.; He, T.; Ren, H. Control of the actin cytoskeleton in root hair development. Plant Sci. 2012, 187, 10–18. [Google Scholar] [CrossRef]
  34. Lin, H.J.; Gao, J.; Zhang, Z.M.; Shen, Y.O.; Lan, H.; Liu, L.; Gao, S.B. Transcriptional responses of maize seedling root to phosphorus starvation. Mol. Biol. Rep. 2013, 40, 5359–5379. [Google Scholar] [CrossRef] [PubMed]
  35. Muñoz-Romero, V.; Benítez-Vega, J.; López-Bellido, L.; López-Bellido, R.J. Monitoring wheat root development in a rainfed vertisol: Tillage effect. Eur. J. Agron. 2010, 33, 182–187. [Google Scholar] [CrossRef]
  36. 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] [Green Version]
  37. Li, L.; Peng, Z.; Mao, X.; Wang, J.; Chang, X.; Reynolds, M.; Jing, R. Genome-wide association study reveals genomic regions controlling root and shoot traits at late growth stages in wheat. Ann. Bot. 2019, 124, 993–1006. [Google Scholar] [CrossRef] [PubMed]
  38. Rowse, H.R.; Goodman, D. Axial resistance to water movement in broad bean (Vicia faba) roots. J. Exp. Bot. 1981, 32, 591–598. [Google Scholar] [CrossRef]
  39. Figueroa-Bustos, V.; Palta, J.A.; Chen, Y.; Siddique, K.H. Characterization of root and shoot traits in wheat cultivars with putative differences in root system size. Agronomy 2018, 8, 109. [Google Scholar] [CrossRef] [Green Version]
  40. Davidson, R.L. Effect of root/leaf temperature differentials on root/shoot ratios in some pasture grasses and clover. Ann. Bot. 1969, 33, 561–569. [Google Scholar] [CrossRef]
  41. Garnier, E. Resource capture, biomass allocation and growth in herbaceous plants. Trends Ecol. Evol. 1991, 6, 126–131. [Google Scholar] [CrossRef]
  42. Qian, H.; Lu, H.; Ding, H.; Lavoie, M.; Li, Y.; Liu, W.; Fu, Z. Analyzing Arabidopsis thaliana root proteome provides insights into the molecular bases of enantioselective imazethapyr toxicity. Sci. Rep. 2015, 5, 11975. [Google Scholar] [CrossRef] [Green Version]
  43. Soto-Cerda, B.J.; Cloutier, S.; Gajardo, H.A.; Aravena, G.; Quian, R. Identifying drought-resilient flax genotypes and related-candidate genes based on stress indices, root traits and selective sweep. Euphytica 2019, 215, 41. [Google Scholar] [CrossRef]
  44. Zhao, Z.; Ge, T.; Gunina, A.; Li, Y.; Zhu, Z.; Peng, P.; Kuzyakov, Y. Carbon and nitrogen availability in paddy soil affect rice photosynthate allocation, microbial community composition, and priming: Combining continuous 13 C labeling with PLFA analysis. Plant Soil 2019, 445, 137–152. [Google Scholar] [CrossRef]
  45. Manschadi, A.M.; Christopher, J.; Devoil, P.; Hammer, G.L. The role of root architectural traits in adaptation of wheat to water-limited environments. Funct. Plant Biol. 2006, 33, 823–837. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Christopher, J.T.; Manschadi, A.M.; Hammer, G.L.; Borrell, A.K. Developmental and physiological traits associated with high yield and stay-green phenotype in wheat. Aust. J. Agric. Res. 2008, 59, 354–364. [Google Scholar] [CrossRef]
  47. Zhang, Y.M.; Jia, Z.; Dunwell, J.M. The applications of new multi-locus GWAS methodologies in the genetic dissection of complex traits. Front. Plant Sci. 2019, 10, 100. [Google Scholar] [CrossRef] [Green Version]
  48. Safdar, L.B.; Andleeb, T.; Latif, S.; Umer, M.J.; Tang, M.; Li, X.; Quraishi, U.M. Genome-wide association study and QTL meta-analysis identified novel genomic loci controlling potassium use efficiency and agronomic traits in bread wheat. Front. Plant Sci. 2020, 11, 70. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Roy, R.; Purty, R.S.; Agrawal, V.; Gupta, S.C. Transformation of tomato cultivar ‘Pusa Ruby’ with bspA gene from Populus tremula for drought tolerance. Plant Cell Tissue Organ Cult. 2006, 84, 56–68. [Google Scholar] [CrossRef]
  50. Manske, G.G.B.; Vlek, P.L.G. Root architecture-wheat as a model plant. In Plant Roots: The Hidden Half; Waisel, Y., Eshel, A., Kafkafi, U., Eds.; Marcel Dekker: New York, NY, USA, 2002; pp. 249–259. [Google Scholar]
  51. Wei, K.; Li, Y. Functional genomics of the protein kinase superfamily from wheat. Mol. Breed. 2019, 39, 141. [Google Scholar] [CrossRef]
  52. Li, Y.; Wei, K. Comparative functional genomics analysis of cytochrome P450 gene superfamily in wheat and maize. BMC Plant Biol. 2020, 20, 2887. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Trewavas, A. How plants learn. Proc. Natl. Acad. Sci. USA 1999, 96, 4216–4218. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Saradadevi, R.; Bramley, H.; Palta, J.A.; Siddique, K.H. Stomatal behaviour under terminal drought affects post-anthesis water use in wheat. Funct. Plant Biol. 2017, 44, 279–289. [Google Scholar] [CrossRef]
  55. Schoppach, R.; Wauthelet, D.; Jeanguenin, L.; Sadok, W. Conservative water use under high evaporative demand associated with smaller root metaxylem and limited trans-membrane water transport in wheat. Funct. Plant Biol. 2014, 41, 257–269. [Google Scholar] [CrossRef]
  56. Kirkegaard, J.A.; Lilley, J.M.; Howe, G.N.; Graham, J.N. Impact of subsoil water use on wheat yield. Aust. J. Agric. Res. 2007, 58, 303–315. [Google Scholar] [CrossRef]
  57. Wasson, A.P.; Richards, R.A.; Chatrath, R.; Misra, S.C.; Prasad, S.S.; Rebetzke, G.J.; Watt, M. Traits and selection strategies to improve root systems and water uptake in water-limited wheat crops. J. Exp. Bot. 2012, 63, 3485–3498. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Lynch, J.P. Steep, cheap and deep: An ideotype to optimize water and N acquisition by maize root systems. Ann. Bot. 2013, 112, 347–357. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. Oyanagi, A.; Sato, A.; Wada, M. Varietal differences in geotropic response of primary seminal root in Japanese wheat. Jpn. J. Crop. Sci. 1991, 60, 312–319. [Google Scholar] [CrossRef] [Green Version]
  60. 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] [Green Version]
  61. Gu, C.; Begley, T.J.; Dedon, P.C. tRNA modifications regulate translation during cellular stress. FEBS Lett. 2014, 588, 4287–4296. [Google Scholar] [CrossRef] [Green Version]
  62. Bai, C.; Liang, Y.; 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]
  63. Hochholdinger, F.; Yu, P.; Marcon, C. Genetic control of root system development in maize. Trends Plant Sci. 2018, 23, 79–88. [Google Scholar] [CrossRef]
  64. Reinert, S.; Kortz, A.; Leon, J.; Naz, A.A. Genome-wide association mapping in the global diversity set reveals new QTL controlling root system and related shoot variation in barley. Front. Plant Sci. 2016, 7, 1061. [Google Scholar] [CrossRef] [Green Version]
  65. Sinha, S.K.; Rani, M.; Kumar, A.; Kumar, S.; Venkatesh, K.; Mandal, P.K. Natural variation in root system architecture in diverse wheat genotypes grown under different nitrate conditions and root growth media. Theor. Exp. Plant Physiol. 2018, 30, 223–234. [Google Scholar] [CrossRef]
  66. Arsenault, J.L.; Poulcur, S.; Messier, C.; Guay, R. WinRHlZO™, a root-measuring system with a unique overlap correction method. Hort. Sci. 1995, 30, 906. [Google Scholar] [CrossRef] [Green Version]
  67. Sinha, S.K.; Rani, M.; Bansal, N.; Venkatesh, K.; Mandal, P.K. Nitrate starvation induced changes in root system architecture, carbon: Nitrogen metabolism, and miRNA expression in nitrogen-responsive wheat genotypes. Appl. Biochem. Biotechnol. 2015, 177, 1299–1312. [Google Scholar] [CrossRef] [PubMed]
  68. Sinha, S.K.; Kumar, A.; Tyagi, A.; Venkatesh, K.; Paul, D.; Singh, N.K.; Mandal, P.K. Root architecture traits variation and nitrate-influx responses in diverse wheat genotypes under different external nitrogen concentrations. Plant Physiol. Biochem. 2020, 148, 246–259. [Google Scholar] [CrossRef] [PubMed]
  69. Nagar, C.K.; Gayatri, A.B.; Mandal, P.K. Nitrogen stress induced changes in root system architecture (RSA) in diverse wheat (T. aestivum L.) genotypes at seedling stage. Wheat Barley Res. 2018, 10, 93–101. [Google Scholar] [CrossRef] [Green Version]
  70. Doyle, J.J.; Doyle, J.L. Isolation of plant DNA from fresh tissue. Focus 1990, 12, 39–40. [Google Scholar]
  71. Pritchard, J.K.; Stephens, M.; Rosenberg, N.A.; Donnelly, P. Association mapping in structured populations. Am. J. Hum. Genet. 2000, 67, 170–181. [Google Scholar] [CrossRef] [Green Version]
  72. Evanno, G.; Regnaut, S.; Goudet, J. Detecting the number of clusters of individuals using the software STRUCTURE: A simulation study. Mol. Ecol. 2005, 14, 2611–2620. [Google Scholar] [CrossRef] [Green Version]
  73. Bradbury, P.J.; Zhang, Z.; Kroon, D.E.; Casstevens, T.M.; Ramdoss, Y.; Buckler, E.S. TASSEL: Software for association mapping of complex traits in diverse samples. Bioinformatics 2007, 23, 2633–2635. [Google Scholar] [CrossRef]
  74. Wang, S.B.; Feng, J.Y.; Ren, W.L.; Huang, B.; Zhou, L.; Wen, Y.J.; Zhang, Y.M. Improving power and accuracy of genome-wide association studies via a multi-locus mixed linear model methodology. Sci. Rep. 2016, 6, 19444. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  75. Tamba, C.L.; Zhang, Y.M. A fast mrMLM algorithm for multi-locus genome-wide association studies. bioRxiv 2018, 341784. [Google Scholar] [CrossRef]
  76. Wen, Y.J.; Zhang, H.; Ni, Y.L.; Huang, B.; Zhang, J.; Feng, J.Y.; Wu, R. Methodological implementation of mixed linear models in multi-locus genome-wide association studies. Brief. Bioinform. 2018, 19, 700–712. [Google Scholar] [CrossRef] [Green Version]
  77. Ren, W.L.; Wen, Y.J.; Dunwell, J.M.; Zhang, Y.M. pKWmEB: Integration of Kruskal–Wallis test with empirical Bayes under polygenic background control for multi-locus genome-wide association study. Heredity 2018, 120, 208–218. [Google Scholar] [CrossRef]
  78. Zhang, J.; Feng, J.Y.; Ni, Y.L.; Wen, Y.J.; Niu, Y.; Tamba, C.L.; Zhang, Y.M. pLARmEB: Integration of least angle regression with empirical Bayes for multilocus genome-wide association studies. Heredity 2017, 118, 517–524. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  79. Tamba, C.L.; Ni, Y.L.; Zhang, Y.M. Iterative sure independence screening EM-Bayesian LASSO algorithm for multi-locus genome-wide association studies. PLoS Comput. Biol. 2017, 13, e1005357. [Google Scholar] [CrossRef] [PubMed]
  80. Wen, Y.J.; Zhang, H.; Ni, Y.L.; Huang, B.; Zhang, J.; Feng, J.Y.; Wang, S.B.; Dunwell, J.M.; Zhang, Y.M.; Wu, R. mrMLM: Multi-Locus Random-SNP-Effect Mixed Linear Model Tools for Genome-Wide Association Study. Available online: https://cran.r-project.org/package=mrMLM (accessed on 2 October 2020).
  81. Voorrips, R.E. MapChart: Software for the graphical presentation of linkage maps and QTLs. J. Hered. 2002, 93, 77–78. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  82. Conesa, A.; Götz, S.; García-Gómez, J.M.; Terol, J.; Talón, M.; Robles, M. Blast2GO: A universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics 2005, 21, 3674–3676. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Correlations between root system architecture traits and shoot traits. TRS, total root size (cm); FOLRN, first order lateral root number; SOLRN, second order lateral root number; LRS, lateral root size; LRD, lateral root density (cm−1); RLD, root length density cm−2; AD, average diameter (mm); RV, root volume (cm3); SRN, seminal root number; SDW, shoot dry weight (mg); RDW, root dry weight (mg); SL, shoot length (cm); RL, root length (cm); RSDWR, root shoot dry weight ratio; SRL, specific root length (cm mg−1). Red boundary for each pairwise correlation displays density ellipses covering 95% points between each variable. The density ellipses are a graphical indicator of the correlation between two variables. It collapses diagonally as the correlation between the two variables approaches either 1 or −1. It is more circular if the two variables are less correlated.
Figure 1. Correlations between root system architecture traits and shoot traits. TRS, total root size (cm); FOLRN, first order lateral root number; SOLRN, second order lateral root number; LRS, lateral root size; LRD, lateral root density (cm−1); RLD, root length density cm−2; AD, average diameter (mm); RV, root volume (cm3); SRN, seminal root number; SDW, shoot dry weight (mg); RDW, root dry weight (mg); SL, shoot length (cm); RL, root length (cm); RSDWR, root shoot dry weight ratio; SRL, specific root length (cm mg−1). Red boundary for each pairwise correlation displays density ellipses covering 95% points between each variable. The density ellipses are a graphical indicator of the correlation between two variables. It collapses diagonally as the correlation between the two variables approaches either 1 or −1. It is more circular if the two variables are less correlated.
Ijms 22 07188 g001
Figure 2. Population stratification of wheat association mapping panel. (A) Delta K, rate of change from 2 to 4. (B) Biplot of principal component analysis (PCA), PCA 1 and PCA 2. (C) Bar plot of the AM panel showing the structuring of two subpopulations (SP1 and SP2) in different colors, viz. SP1 (red, 40 genotypes), SP2 (green, 96 genotypes).
Figure 2. Population stratification of wheat association mapping panel. (A) Delta K, rate of change from 2 to 4. (B) Biplot of principal component analysis (PCA), PCA 1 and PCA 2. (C) Bar plot of the AM panel showing the structuring of two subpopulations (SP1 and SP2) in different colors, viz. SP1 (red, 40 genotypes), SP2 (green, 96 genotypes).
Ijms 22 07188 g002
Figure 3. Distribution of the QTNs (quantitative trait nucleotides) identified using five different multi-locus GWAS (genome wide association study) models. (A) Number of significant QTNs detected for 15 traits across five multi-locus GWAS methods. (B) Number of significant QTNs detected using each of five multi-locus GWAS methods. (C) Number of QTNs per chromosome.
Figure 3. Distribution of the QTNs (quantitative trait nucleotides) identified using five different multi-locus GWAS (genome wide association study) models. (A) Number of significant QTNs detected for 15 traits across five multi-locus GWAS methods. (B) Number of significant QTNs detected using each of five multi-locus GWAS methods. (C) Number of QTNs per chromosome.
Ijms 22 07188 g003
Figure 4. Chromosome location of 87 reliable QTNs (quantitative trait nucleotides) for 15 root and shoot traits. The color bars represent QTNs identified for different traits, red = SOLRN (second order lateral root number); dark green = TRS; (total root size) light green = RV (root volume); orange = AD (average diameter); pink = RDW (root dry weight); light blue = RLD (root length density); blue = SL (shoot length). (For interpretation of color references in this figure legend, refer to the web version of this article).
Figure 4. Chromosome location of 87 reliable QTNs (quantitative trait nucleotides) for 15 root and shoot traits. The color bars represent QTNs identified for different traits, red = SOLRN (second order lateral root number); dark green = TRS; (total root size) light green = RV (root volume); orange = AD (average diameter); pink = RDW (root dry weight); light blue = RLD (root length density); blue = SL (shoot length). (For interpretation of color references in this figure legend, refer to the web version of this article).
Ijms 22 07188 g004aIjms 22 07188 g004b
Figure 5. Boxplot for seven reliable QTNs with significant effects (p < 0.01) on corresponding root and shoot traits. For each QTNs, the germplasm lines were divided into two groups according to superior and inferior allele type. The x-axis represents the two alleles for each QTNs, while the y-axis corresponds to phenotypic value. Subfigures, (AC) represent allelic differences for SOLRN and superior alleles for QTNs Q.SOLRN-1DS, Q.SOLRN-1BL and Q.SOLRN-6AL were T, A and T respectively. Subfigures (DG) represent allelic differences for LRS, RV, SDW and AD for QTNs Q.LRS-7BS, Q.RV-5BL, Q.SDW-4AL and Q.AD-6AL with superor alleles as T, C, T and C respectively.
Figure 5. Boxplot for seven reliable QTNs with significant effects (p < 0.01) on corresponding root and shoot traits. For each QTNs, the germplasm lines were divided into two groups according to superior and inferior allele type. The x-axis represents the two alleles for each QTNs, while the y-axis corresponds to phenotypic value. Subfigures, (AC) represent allelic differences for SOLRN and superior alleles for QTNs Q.SOLRN-1DS, Q.SOLRN-1BL and Q.SOLRN-6AL were T, A and T respectively. Subfigures (DG) represent allelic differences for LRS, RV, SDW and AD for QTNs Q.LRS-7BS, Q.RV-5BL, Q.SDW-4AL and Q.AD-6AL with superor alleles as T, C, T and C respectively.
Ijms 22 07188 g005
Figure 6. Flow diagram of phenotyping protocol. (A) Sterilized wheat seeds for germination; (B) germinated seeds after two days; (C) transferring healthy seedlings to pot containing perlite and vermiculite; (D) transplanted seedlings; (E) different stages of seedlings growth under controlled condition; (F) manual measurement of root and shoot length; (G) root images captured by using flatbed scanner; (H) scanned images of roots. The experimental design adopted for screening was a completely randomized design (CRD) with three replications including two checks (C306 and HD2967). For each replication five biological replicates were used in each pot. The experiment was done in batches for 20 genotypes at each time. Before germination, 20 seeds of each accession were washed carefully and thoroughly with double distilled water and surface sterilized with 0.5% Sodium hypochlorite solution for 30 s. The sterilized seeds were then washed with double distilled water three to four times to remove any trace of adhering chemicals. The seeds were placed well spread in a thoroughly moist germination paper/filter paper taken in a petri dish and allowed to germinate under a growth chamber at 22 ± 1 °C room temperature in the dark.
Figure 6. Flow diagram of phenotyping protocol. (A) Sterilized wheat seeds for germination; (B) germinated seeds after two days; (C) transferring healthy seedlings to pot containing perlite and vermiculite; (D) transplanted seedlings; (E) different stages of seedlings growth under controlled condition; (F) manual measurement of root and shoot length; (G) root images captured by using flatbed scanner; (H) scanned images of roots. The experimental design adopted for screening was a completely randomized design (CRD) with three replications including two checks (C306 and HD2967). For each replication five biological replicates were used in each pot. The experiment was done in batches for 20 genotypes at each time. Before germination, 20 seeds of each accession were washed carefully and thoroughly with double distilled water and surface sterilized with 0.5% Sodium hypochlorite solution for 30 s. The sterilized seeds were then washed with double distilled water three to four times to remove any trace of adhering chemicals. The seeds were placed well spread in a thoroughly moist germination paper/filter paper taken in a petri dish and allowed to germinate under a growth chamber at 22 ± 1 °C room temperature in the dark.
Ijms 22 07188 g006
Table 1. Root system architecture and shoot traits along with their basic statistics of the wheat accessions in the diverse association panel.
Table 1. Root system architecture and shoot traits along with their basic statistics of the wheat accessions in the diverse association panel.
S. No.Trait Name AbbreviationTrait DescriptionUnitRangeMean ± SECV (%)
1Total root size TRSSum of the path length of seminal including primary roots and lateral roots (LRs)cm65.10–378.63242.94 ± 5.3526.06
2First order lateral root numberFOLRNNumber of first order LRs (emerging from primary and seminal roots)-14.33–240.6795.22 ± 3.5844.50
3Second order lateral root numberSOLRNNumber of second order LRs (emerging from first-order LRs)-6.67–490.33165.57 ± 8.9764.10
4Lateral root sizeLRS First order length + second order length divided by total root size (TRS)-0.02–0.870.39 ± 0.0138.39
5Lateral root densityLRDFirst-order LR no. divided by PRPcm−10.81–5.892.81 ± 0.0835.67
6Root length densityRLDTRS divided by volume of pot (500 cm3)cm−20.11–0.760.49 ± 0.0126.30
7Average diameterADProjected area divided by TRSmm 0.31–0.490.38 ± 0.018.46
8Root volumeRVπ × (Half of the avg. diameter/2) 2 × TRScm30.02–0.350.16 ± 0.0137.62
9Seminal root numberSRNPrimary root bursts through the coleorhiza -2.67–6.005.06 ± 0.0410.29
10Shoot dry weightSDWShoot dry weight mg5.02–60.4135.58 ± 0.8829.18
11Root dry weight RDWRoot dry weight mg1.07–23.0610.76 ± 0.3235.14
12Shoot lengthSLShoot length cm13.57–40.7730.61 ± 0.4115.97
13Root lengthRLRoot length cm9.67–37.5024.81 ± 0.5425.83
14Root shoot dry weight ratioRSDWRRatio of root to shoot dry weight-0.16–0.590.31 ± 0.0126.07
15Specific root lengthSRLTotal root length per unit root dry masscm mg−11.32–9.102.52 ± 0.0734.02
Table 2. Distribution of 16,616 SNPs and identified linkage disequilibrium (LD) on 21 chromosomes in 136 wheat genotypes.
Table 2. Distribution of 16,616 SNPs and identified linkage disequilibrium (LD) on 21 chromosomes in 136 wheat genotypes.
ChromosomeSize (Mb)No. of SNPAverage Number of SNPs per (Mb)Average LD (r2)No. of Marker Pairs in Perfect LD (r2 = 1)
1A594.19111.530.2062134
1B689.859871.430.2683511
1D495.458951.800.2983923
2A780.89521.220.241999
2B801.2611781.470.159819
2D651.8510371.590.2261813
3A750.846630.880.114274
3B830.8310031.200.149368
3D615.555790.940.098121
4A744.595730.760.122310
4B673.624470.660.111197
4D509.862870.560.08773
5A709.777371.030.122368
5B713.1511251.570.1891769
5D566.087691.350.147713
6A618.086571.060.157655
6B720.998461.170.147354
6D473.596221.310.112255
7A736.718621.170.122401
7B750.627571.000.126269
7D638.697291.140.099202
Table 3. Common quantitative trait nucleotides (QTNs) of root and shoot traits identified using different multi-locus genome-wide association studies (GWAS) methods.
Table 3. Common quantitative trait nucleotides (QTNs) of root and shoot traits identified using different multi-locus genome-wide association studies (GWAS) methods.
S. NoQTNTraitMarkerAlleleCHRPhysical Position (bp)LOD ScoreR2 (%)Method
1Q.TRS-2ASTRSAX-94952472G/C2AS81817943.29–3.759.77–17.721,5
2Q.TRS-4ALTRSAX-95105488A/C4AL6319033543.38–5.018.82–10.852,4
3Q.TRS-7BSTRSAX-95119337G/A7BS1388827916.81–10.7317.10–22.222,4
4Q.FOLRN-7ASFOLRNAX-95249973G/A7AS549979933.10–4.925.22–6.831,2
5Q.FOLRN-7BSFOLRNAX-95123855T/C7BS996350313.59–4.416.12–8.624,5
6Q.SOLRN-1ALSOLRNAX-95102105A/G1AL4745732303.24–4.581.68–3.651,2,4
7Q.SOLRN-1BLSOLRNAX-95077960G/A1BL5723014733.55–8.972–4.151,2,4,5
8Q.SOLRN-1DSSOLRNAX-94448890T/C1DS107416984.44–6.922.47–4.121,2,4,5
9Q.SOLRN-3BLSOLRNAX-94516395C/T3BL7386992254.34–5.582.61–3.991,2,4
10Q.SOLRN-3BLSOLRNAX-95150666G/C3BL66401944110.31–14.8427.09–38.601,2,4,5
11Q.SOLRN-3DLSOLRNAX-94450588A/G3DL3673098938.17–10.8421.14–26.164,5
12Q.SOLRN-6ALSOLRNAX-95244609T/C6AL5990351594.61–6.781.88–3.291,2,4,5
13Q.SOLRN-7BLSOLRNAX-94564853A/G7BL6739626834.30–5.989.92–18.501,2,4
14Q.SOLRN-7BSSOLRNAX-95123855T/C7BS996350314.49–4.691.60–2.174,5
15Q.LRS-2DSLRSAX-94471646A/G2DS139297123.31–7.301–2.402,4
16Q.LRS-3ASLRSAX-94961347C/T3AS1846352017.15–15.487.47–12.411,4,5
17Q.LRS-3BLLRSAX-94480990G/T3BL4172613676.63–15.771–6.152,4
18Q.LRD-5BSLRSAX-94706358A/C5BS691342153.41–13.253.47–8.201,4,5
19Q.LRS-6ALLRSAX-95244609T/C6AL5990351593.27–4.210.12–0.24,5
20Q.LRD-7BLLRSAX-94528392G/A7BL6753144954.92–9.1816.43–30.671,4,5
21Q.LRD-7BSLRSAX-95123855T/C7BS996350313.35–17.244.17–8.351,2,4,5
22Q.LRD-4ALLRDAX-94841365G/A4AL7257514994.90–5.5612.04–13.404,5
23Q.LRD-5DLLRDAX-94571501T/C5DL4726357823.46–5.695.74–11.603,4,5
24Q.LRD-7BLLRDAX-94763902T/C7BL5521208413.58–6.436.83–10.022,4,5
25Q.RLD-1BLRLDAX-94728516T/C1BL4912414766.01–28.190.45–4.154,5
26Q.RLD-2ASRLDAX-94952472G/C2AS81817943.67–9.520.98–4.234,5
27Q.RLD-3BSRLDAX-94516395C/T3BL7386992253.36–3.461.66–1.934,5
28Q.RLD-3DLRLDAX-94544285C/G3DL5944788105.88–21.543.94–6.332,4
29Q.RLD-4ALRLDAX-95105488A/C4AL6319033543.26–8.065.97–15.441,4
30Q.RLD-5ALRLDAX-95659861G/A5AL5657542303.49–6.681.89–2.572,5
31Q.RLD-5DSRLDAX-95160709G/T5DS104719643.49–13.730.16–4.384,5
32Q.RLD-6BLRLDAX-94502864G/A6BL7103245429.73–14.541.14–8.224,5
33Q.RLD-7BSRLDAX-95119337G/A7BS1388827916.71–11.2312.20–17.081,4,5
34Q.AD-1ALADAX-94932678A/G1AL5579540563.67–6.822.42–5.924,5
35Q.AD-1BLADAX-94620468A/G1BL5863026985.16–8.863.31–7.792,4,5
36Q.AD-1BSADAX-94668789G/A1BS965155303.67–5.970.17–1.532,5
37Q.AD-2DLADAX-94647816A/C2DL6419801264.58–22.753.53–7.082,4
38Q.AD-3ASADAX-94996868C/T3AS10305066.25–16.451.67–8.542,4
39Q.AD-3BSADAX-94886515T/C3BS1294379303.90–9.271.55–2.722,5
40Q.AD-3DSADAX-94726849C/T3DS813910816.70–7.870.21–2.182,4
41Q.AD-4ALADAX-94740863T/C4AL7258026753.10–5.140.12–3.544,5
42Q.AD-6ALADAX-95244609T/C6AL5990351593.38–10.623.03–5.882,3,4,5
43Q.AD-6BSADAX-94691127G/A6BS2884195504.97–5.800.18–1.284,5
44Q.AD-6DLADAX-94853162G/A6DL4702945933.70–4.510.39–1.654,5
45Q.RV-2BSRVAX-95227366G2BS460892203.25–1152.6311.40–13.561,2,4
46Q.RV-4ALRVAX-95105488A/C4AL6319033543.97–1166.092.06–8.081,2,3
47Q.RV-4BSRVAX-94648074C/T4BS1409297243.45–1122.582.72–4.811,2,4
48Q.RV-5BLRVAX-94439232C/A5BL2876241523.80–1152.241.47–5.311,2,4,5
49Q.RV-6ALRVAX-95162623C/T6AL6118514053.76–7.560.18–2.614,5
50Q.RV-6ASRVAX-94687596A/G6AS9590415.10–7.870.97–2.204,5
51Q.RV-6ASRVAX-94893701A/G6AS248574775.21–9.450.70–1.234,5
52Q.RV-6BLRVAX-94502864G/A6BL7103245424.54–1136.810.52–5.932,4,5
53Q.RV-7ALRVAX-94664277C/T7AL6899207464.61–1155.165.53–11.201,2,4,5
54Q.SRN-2DLSRNAX-94551988A/G2DL5799651384.27–7.946.83–17.741,3,4,5
55Q.SRN-3BSSRNAX-94486290T/C3BL4778858355.68–6.685.82–6.511,2
56Q.SRN-3DLSRNAX-94382595C/T3DL3929566613.19–4.400.31–8.322,4
57Q.SRN-4DLSRNAX-95157799A/G4DL4612180175.66–7.1917.38–20.911,2
58Q.SRN-6ALSRNAX-94637211A/G6AL6115766063.08–4.8711.77–14.951,2
59Q.SRN-6ALSRNAX-94790960C/A6AL6001310556.19–9.1323.20–28.081,2
60Q.SDW-3DSSDWAX-94439998T/C3DS850609154.29–4.298.15–9.741,2
61Q.SDW-4ALSDWAX-94470023T/C4AL5812164904.49–6.618.09–19.051,2,3,4,5
62Q.SDW-4BSSDWAX-95630040A/G4BS128936143.25–8.405.43–11.573,4,5
63Q.SDW-5BLSDWAX-94484139T/C5BL2895670144.72–4.875.64–5.924,5
64Q.SDW-7ALSDWAX-94664277C/T7AL6899207463.36–4.034.39–5.304,5
65Q.SDW-7BLSDWAX-95165787C/A7BL5910534493.11–3.136.54–7.064,5
66Q.RDW-1BLRDWAX-94925225T/C1BL3729962414.49–4.974.87–7.594,5
67Q.RDW-5ALRDWAX-94501549C/T5AL6723395964.29–4.914.71–7.452,3,4
68Q.RDW-6ASRDWAX-94809955C/T6AS6316984.67–5.394.37–5.672,4
69Q.RDW-7ASRDWAX-94386260G/A7AS2690939213.65–4.3518.35–18.601,2
70Q.SL-1ALSLAX-94929421T/C1AL5503310925.63–5.665.16–5.391,2
71Q.SL-2ALSLAX-95076063T/C2AL5558203553.52–4.913.68–8.241,3,5
72Q.SL-4BSSLAX-95012217C/T4BS979225354.71–6.365.06–5.741,5
73Q.SL-5BLSLAX-95092434C/G5BL5791535794.16–5.205.57–8.611,2
74Q.SL-7BSSLAX-94876335G/A7BS1444849194.04–5.487.23–8.271,2
75Q.RL-1DSRLAX-94448890T/C1DS107416984.02–5.523.67–7.891,2,5
76Q.RL-2BSRLAX-95227366G2BS460892203.68–5.778.59–15.761,4
77Q.RL-6BLRLAX-94502864G/A6BL7103245424.63–6.064.68–9.411,2,3,4,5
78Q.RL-6BSRLAX-94539094A/C6BS83875283.52–6.772.83–17.151,2,5
79Q.RL-6BSRLAX-94900754A/G6BS1577928433.27–5.445.26–9.161,3
80Q.RL-7ALRLAX-95241843G/T7AL6104980246.31–6.505.55–5.971,2
81Q.RL-7BLRLAX-94528392G/A7BL6753144953.21–8.5215.42–30.831,5
82Q.RL-7DLRLAX-94861078A/G7DL6142760513.36–5.013.12–4.243,4
83Q.RSDWR-2ALRSDWRAX-94430108C/T2AL5093577865.19–14.680.55–1.382,4
84Q.RSDWR-2BLRSDWRAX-94856412C/T2BL7548643634.85–103.960.58–1.954,5
85Q.RSDWR-3BSRSDWRAX-94691217A/G3BS80299678.14–83.910.02–3.222,4
86Q.SRL-5ALSRLAX-94748697A/G5AL4065348463.37–6.014.81–8.102,4,5
87Q.SRL-2BLSRLAX-94457792T/C2BL5760834718.12–19.4064.33–31.921,2,4,5
TRS, total root size (cm); FOLRN, first order lateral root number; SOLRN, second order lateral root number; LRS, lateral root size; LRD, lateral root density (cm−1); RLD, root length density cm−2; AD, average diameter (mm); RV, root volume (cm3); SRN, seminal root number; SDW, shoot dry weight (mg); RDW, root dry weight (mg); SL, shoot length (cm); RL, root length (cm); RSDWR, root shoot dry weight ratio; SRL, specific root length (cm mg−1). ML-GWAS methods (MrMLM-1, FASTMrMLM-2, FASTMrEMMA-3, pLARmEB-4, ISIS EM-BLASSO-5).
Table 4. Pleiotropic quantitative trait nucleotides (QTNs).
Table 4. Pleiotropic quantitative trait nucleotides (QTNs).
S. No.MarkerTraitsChromosomePosition (Mb)
1AX-94448890SOLRN, RL1DS10.7417
2AX-94502864RLD, RV, RL6BL710.3245
3AX-94516395SOLRN, RLD3BL738.6992
4AX-94528392LRS, RL7BL675.3145
5AX-94664277RV, SDW7AL689.9207
6AX-94952472TRS, RLD2AS8.181794
7AX-95105488TRS, RLD, RV4AL631.9034
8AX-95119337TRS, RLD7BS138.8828
9AX-95123855FOLRN, SOLRN, LRS7BS99.63503
10AX-95227366RV, RL2BS46.08922
11AX-95244609SOLRN, LRS, AD6AL599.0352
TRS, total root size (cm); FOLRN, first order lateral root number; SOLRN, second order lateral root number; LRS, lateral root size; RLD, root length density cm−2; AD, average diameter (mm); RV, root volume (cm3); SDW, shoot dry weight (mg); RL, root length (cm).
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Danakumara, T.; Kumari, J.; Singh, A.K.; Sinha, S.K.; Pradhan, A.K.; Sharma, S.; Jha, S.K.; Bansal, R.; Kumar, S.; Jha, G.K.; et al. Genetic Dissection of Seedling Root System Architectural Traits in a Diverse Panel of Hexaploid Wheat through Multi-Locus Genome-Wide Association Mapping for Improving Drought Tolerance. Int. J. Mol. Sci. 2021, 22, 7188. https://doi.org/10.3390/ijms22137188

AMA Style

Danakumara T, Kumari J, Singh AK, Sinha SK, Pradhan AK, Sharma S, Jha SK, Bansal R, Kumar S, Jha GK, et al. Genetic Dissection of Seedling Root System Architectural Traits in a Diverse Panel of Hexaploid Wheat through Multi-Locus Genome-Wide Association Mapping for Improving Drought Tolerance. International Journal of Molecular Sciences. 2021; 22(13):7188. https://doi.org/10.3390/ijms22137188

Chicago/Turabian Style

Danakumara, Thippeswamy, Jyoti Kumari, Amit Kumar Singh, Subodh Kumar Sinha, Anjan Kumar Pradhan, Shivani Sharma, Shailendra Kumar Jha, Ruchi Bansal, Sundeep Kumar, Girish Kumar Jha, and et al. 2021. "Genetic Dissection of Seedling Root System Architectural Traits in a Diverse Panel of Hexaploid Wheat through Multi-Locus Genome-Wide Association Mapping for Improving Drought Tolerance" International Journal of Molecular Sciences 22, no. 13: 7188. https://doi.org/10.3390/ijms22137188

APA Style

Danakumara, T., Kumari, J., Singh, A. K., Sinha, S. K., Pradhan, A. K., Sharma, S., Jha, S. K., Bansal, R., Kumar, S., Jha, G. K., Yadav, M. C., & Prasad, P. V. V. (2021). Genetic Dissection of Seedling Root System Architectural Traits in a Diverse Panel of Hexaploid Wheat through Multi-Locus Genome-Wide Association Mapping for Improving Drought Tolerance. International Journal of Molecular Sciences, 22(13), 7188. https://doi.org/10.3390/ijms22137188

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