Identification of Candidate Genes for Root Traits Using Genotype–Phenotype Association Analysis of Near-Isogenic Lines in Hexaploid Wheat (Triticum aestivum L.)

Global wheat (Triticum aestivum L.) production is constrained by different biotic and abiotic stresses, which are increasing with climate change. An improved root system is essential for adaptability and sustainable wheat production. In this study, 10 pairs of near-isogenic lines (NILs)—targeting four genomic regions (GRs) on chromosome arms 4BS, 4BL, 4AS, and 7AL of hexaploid wheat—were used to phenotype root traits in a semi-hydroponic system. Seven of the 10 NIL pairs significantly differed between their isolines for 11 root traits. The NIL pairs targeting qDSI.4B.1 GR varied the most, followed by the NIL pair targeting qDT.4A.1 and QHtscc.ksu-7A GRs. For pairs 5–7 targeting qDT.4A.1 GR, pair 6 significantly differed in the most root traits. Of the 4 NIL pairs targeting qDSI.4B.1 GR, pairs 2 and 4 significantly differed in 3 and 4 root traits, respectively. Pairs 9 and 10 targeting QHtscc.ksu-7A GR significantly differed in 1 and 4 root traits, respectively. Using the wheat 90K Illumina iSelect array, we identified 15 putative candidate genes associated with different root traits in the contrasting isolines, in which two UDP-glycosyltransferase (UGT)-encoding genes, TraesCS4A02G185300 and TraesCS4A02G442700, and a leucine-rich repeat receptor-like protein kinase (LRR-RLK)-encoding gene, TraesCS4A02G330900, also showed important functions for root trait control in other crops. This study characterized, for the first time, that these GRs control root traits in wheat, and identified candidate genes, although the candidate genes will need further confirmation and validation for marker-assisted wheat breeding.


Introduction
Wheat (Triticum aestivum L.) is an important cereal crop, which is cultivated worldwide and significantly for international food markets. The ever-changing global climate has increased abiotic and biotic stresses, constraining wheat production. Considerable attention has been paid to wheat root systems in recent years in an attempt to adapt to climate change and achieve sustainable yields. Roots are directly associated with soil, water, and nutrient uptake, anchoring plants in the substrate, interacting with plant microbes, and storing resources; therefore, root architecture greatly contributes to crop yield [1]. Wheat genotypes with modified root mass had higher water uptake than those with standard root mass, resulting in a 200 kg ha −1 yield advantage in drought conditions [2]. Under nutrientdeficit conditions in the soil, fast-and early-proliferating roots with long and dense root hairs enhanced macro-and micronutrient uptake in wheat; vigorous root growth helps in early nutrient acquisition [3]. Therefore, the root system is essential for improving stress tolerance and yield potential. Identifying genetic variation and genes responsible for root traits in wheat is of great interest.
Single-nucleotide polymorphism (SNP) is a powerful tool for studying diversity within the individual species at the genomic level. Due to its abundance, ubiquity, and amenability to high-throughput automation [4], SNP is now widely used in quantitative trait loci (QTL) mapping, marker-assisted selection (MAS), genome-wide association studies (GWAS), and genomic selection (GS) [5,6]. The wheat 90K Infinium iSelect SNP array has been used in many wheat studies; several grain-related genes and potential markers have been identified using this approach [7,8]. However, few studies have focused on identifying genes associated with root traits in wheat. Rufo et al. (2020) found that higher genetic diversity in root traits is mostly associated with the A (46%) and B (48%) genomes compared with the D (6%) genome in wheat. They identified 31 candidate genes for root system architecture in landraces and modern wheat varieties using a GWAS study [9]. Two root trait architecture candidate genes, AX-110092021 (encoding protein kinase) and AX-95681859 (encoding E3 ubiquitin-protein ligase), were identified on chromosomes 5BL and 1BL in Chinese wheat cultivars [10]. Protein kinase genes were expressed stably on wheat roots at different stages, and played important roles in root anisotropic growth [11]. VERNALIZATION1 (VRN1), a MADS-box protein-encoding gene on chromosome 5B, is associated with root length variations at different wheat growth stages [12]. An earlier study reported 63 marker-trait associations with different root traits, including rooting depth (RD), root mass (RM), and root diameter (Rdia); the height-reducing gene Rht encoding a DELLA protein played an important role in increasing RD but decreasing RM and Rdia [13]. A malate transporter gene, ALMT1 (encoding a novel aluminum (Al)-activated malate transport protein) on chromosome 4DL, is expressed in root apices and improves Al stress tolerance [14]. The recently available reference genome has opened an opportunity to identify trait-controlling genes in wheat, including their sequence and diversity in molecular breeding [15]. Due to difficulties in studying root systems, wheat root phenotyping studies are limited to early growth stage under controlled environments [16]. Potential genes for many root traits from pure lines are yet to be discovered for root improvement.
Near-isogenic lines (NILs) are pairs of genetically identical lines, except for target genetic loci, which are important genetic resources for identifying candidate genes and their functions associated with phenotypic traits. The NILs used in the present study were developed targeting different genomic regions (GRs) on wheat chromosome arms 4BS [17], 4AS [18], 4BL [19], and 7AL [20]. The NILs were developed using the heterogeneous inbred family (HIF) method [19], and therefore, the genetic backgrounds of each pair are different even if they target the same GR, which allowed comparisons among the NIL pairs to narrow down the candidates of the trait-linked markers and genes. In this study, we investigated root traits in 10 pairs of wheat NILs using the wheat 90K Infinium iSelect SNP array as a genotyping method. The study aimed to (1) characterize phenotypic variabilities between isolines for different root traits, (2) identify candidate genes associated with root traits contrasting in NIL pairs, and (3) investigate functional relationships between the identified candidate genes and previously reported genes on the respective GRs.
For Pairs 9 and 10, NIL9a had significantly higher RDCL1 than NIL9b, and NIL10a had significantly higher RL (Supplementary Figure S1), root diameter (Rdia), and RDCL2 than NIL10b. However, NIL10b had significantly more nodal root number per plant (NNR) than NIL10a (Table 2).    Significant correlation was found between different root and shoot traits of the NIL pairs. According to the correlation coefficient values, strong (r ≥ 0.7) and weak (r ≤ 0.3) correlations were denoted [24]. Correlation coefficient values between 0.3 and 0.7 were denoted as moderate correlation. Total root length showed strong positive correlation with RM, RSA2, RDCL1, and NLP; moderate positive correlation with RD, RSR, NNR, and SH; and moderate but negative correlation with Rdia. Root mass showed strong positive correlation with RD, RDCL1, and NNR, and moderate positive correlation with RSR, RSA2, SH, and NLP. Strong positive correlation was found between RD and NNR, whereas moderate but negative correlation was found between RD and Rdia. Other traits, such as RSA1 and RDCL2, RSA2 and RDCL1, RSA2 and NLP, and RDCL1 and NLP, also showed strong positive correlations (Supplementary Figure S2).

SNPs and Candidate Genes
Of all the contrasting NIL pairs, a total of 67 SNPs on the targeted chromosomes were identified as candidate tightly linked markers for the studied root traits; for each investigated GR, they showed common polymorphism between the isolines among all the pairs targeting the same GR (Supplementary Table S2); specifically, pairs 1-4, 5-7, and 9-10 showed polymorphism for 11, 50, and 6 SNPs, respectively (Supplementary Table S2). In pairs 1-4, 5-7, and 9-10, 5, 11 and 2 SNPs showed polymorphism within the respective GRs. Using those 67 polymorphic SNP markers, from the nine NIL pairs, a total of 181 candidate genes were identified (Supplementary Table S2). However, only those genes with functions that have been previously reported to be responsible for root traits in wheat and other crops (especially cereal crops) were considered as important candidates (Table 3). Fifteen candidate genes tightly linked with the SNPs were identified (Table 3). On GR qDSI.4B.1, TraesCS4B01G114500 and TraesCS4B01G114800, associated with the marker Excalibur_c869_2091 ( Figure 2), were the genes annotated to encode the Phox-associated domain and bHLH-MYC and R2R3-MYB transcription factors N-terminal proteins, respectively (Table 3). Table 3. Candidate genes identified in the targeted QTL regions are strongly suggested to be responsible for the tested root and shoot traits in the nine pairs of wheat near-isogenic lines (NILs) used in this study.  For pairs 5-7, the TraesCS4A02G163700 gene encoded the inositol-1-monophosphatase family protein. An LRR-RLK-encoding gene, TraesCS4A02G330900, was located very close (7.78 kb away) to the marker BobWhite_c27287_232 (Table 3; Figure 3). The gene TraesCS4A02G331000 was also linked (17.30 kb away) to BobWhite_c27287_232. Two genes, TraesCS4A02G215800 and TraesCS4A02G215700, were 202.63 and 214.39 kb away from BS00098868_51, respectively (Table 3; Figure 3); these genes encode the F-box family protein. The UDP-glucoronosyl and/or UDP-glucosyl transferase-encoding gene TraesCS4A02G185300 was associated with BS00094406_51 (962.14 kb away), and TraesCS4A02G442700 was associated with CAP12_c3789_390 (11.58 kb away) and RAC875_c51781_771 (10.29 kb away) ( Table 3). TraesCS4A02G323500 and TraesCS4A02G323400 were closely associated with Tdurum_contig12899_342 (5.86 and 10.11 kb away, respectively), and TraesCS4A02G000300 was closely associated with Jagger_c9935_144 (90.30 kb away), encoding serine-threonine/tyrosine-protein kinase (Table 3; Figure 3).

NIL
For pairs 9-10, two genes linked to SNP CAP8_c3496_118, namely, TraesCS7A02G431600 and TraesCS7A02G431500, had root-trait-related functions as F-box protein and serinethreonine/tyrosine-protein kinase, respectively (Table 3; Figure 4).  The genes were identified based on the 90K wheat microarray assessment of near-isogenic lines with contrasting root traits. SNP markers are located in different positions of the genes. Structural information of the genes and SNPs was from a wheat genome database (https://urgi.versailles.inra.fr/jbrowseiwgsc/ gmod_jbrowse/, accessed on 2 March 2021). The physical positions of the linking markers of the QTL were found from the original mapping study using parental cultivar; therefore, in the reference genome (Chinese Spring), the QTL extends in both the long and short arm of chromosome 7A, which is unusual. Chromosomal fragment translocation might be the possible reason for the difference between the reference genome and the parental cultivars used in the original mapping study [25]. The measuring bar indicates the DNA length in kilo base pairs (kb).

NIL Pairs Show Different Phenotypic Performances
Among the 10 NIL pairs, 7 pairs showed significant differences between the isolines for the studied root traits. The HIF method for NIL development allowed the isolines of a NIL pair to be homozygous apart from the target locus [20]. Therefore, the isolines showed contrasting performances, suggesting that they contain the gene(s) controlling the corresponding root traits. NIL pairs that were developed from the same parents showed different phenotypic performances among the pairs. Due to the use of the HIF method in NIL development, each NIL pair was from a single seed descent (SSD) of F2 individuals, and for that reason [17], the NIL pairs had completely different genetic backgrounds from each other. Therefore, despite having the same parents, individual pairs showed different root phenotypic performances. By comparison, among all the pairs targeting the same GR, we narrowed down in each GR the candidate SNPs that could be linked to the root traits.

Targeted Genomic Regions of NILs Overlap Previously Reported QTL for Root Traits and Other Yield-Related Traits
The NILs used in this study to investigate root traits were initially developed to target different stress tolerance QTL ( Table 1). Some of the QTL overlapped previously reported QTL for root or yield-related traits. Under drought stress, rooting depth QTL (qMRL.4B.1), root mass QTL (qTRB.4B.1), and grain yield QTL (qMRL.4B.2)-identified through linkage mapping of the hexaploid wheat population C306 × WL711-overlapped with qDSI.4B.1 in the current study [21]. Under well-watered and low-nitrogen (N) conditions, two root mass QTL (physical position, 75.74-222.28 Mb) were reported in the Xiaoyan 54 × Jing 411 population [26], which are very close to qDSI.4B.1 (59.61-75.74 Mb). From the meta-QTL analysis in durum wheat, Root_MQTL_52 (35.7-74.7 Mb) controlling rooting depth, root volume, root mass, root angle, root number, and root/shoot ratio overlapped the 4BS QTL [27]. Root_MQTL_51 for total root length overlapped the same QTL [27]. A QTL (QRv.sau-4B) controlling root volume was reported on chromosome 4BS (105.15 Mb) in bread wheat, only 29.4 Mb away from our studied qDSI.4B.1 [28]. Besides the root trait QTL, a few other yield-related QTL, such as QTL for thousand-kernel weight, were also found on chromosome 4BS as reported earlier [29,30].
A meta-QTL, Root_MQTL_86 for root mass, overlapped QHtscc.ksu-7A targeted in the current study, and another meta-QTL, Root_MQTL_85, was 10 Mb away from it [27]. Other grain-related QTL in cereals overlapped marker Xbarc121, a flanking marker of QHtscc.ksu-7A. Several QTL controlling thousand-grain weight and grain number in bread wheat for grain weight per plant in durum wheat, panicle number, seed setting rate, and grain weight per plant in rice also overlapped Xbarc121 and therefore QHtscc.ksu-7A [35,36].
The overlapping of our targeted QTL with other reported QTL for root, stress tolerance, and yield-related traits indicates that the genes underlying these GRs might have multiple functions in controlling root, yield, and stress tolerance traits in wheat. Therefore, characterization of these GRs might lead to the identification of important candidate genes for the traits.

Putative Candidate Genes Controlling Wheat Root Traits
The identified candidate genes reportedly have functional similarities for controlling root traits in different crops, including wheat ( Table 4). The most outstanding candidates were the UDP-glycosyltransferase (UGT)-encoding genes TraesCS4A02G185300 and TraesCS4A02G442700 and the LRR-RLK-encoding gene TraesCS4A02G330900. UGT has the potential to improve crop root growth and development through the glycosylation of different hormones, including auxin, cytokinin. and abscisic acid [37]. Lignification of a crop cell is linked to monolignol glycosylation and co-expressed by the UGT gene, which had more lateral root formation and higher primary root elongation than the wild type despite the application of a strong root growth inhibitor (cytokinin) [38]. Overexpression of the UGT-encoding gene OsIAGT1 reduced root length due to the upregulation of auxin synthesis genes and decreased endogenous indole acetic acid in rice seedlings [39]. Like UGT, LRR-RLK plays a crucial role in root growth and development, but there are no reports on the LRR-RLK gene controlling root traits in wheat. However, the LRR-RLK genes TaLRRKs are expressed highly in wheat roots [40]. Homozygous T-DNA insertion lines of Arabidopsis developed for LRR-RLK expression demonstrated root growth control through cytokinin mediation [41]; in other cases, a mutant with two T-DNA insertions (BAK1/SERK3 and IRK genes) had inhibited root growth and rlk902 had reduced root length [41]. Transgenic Arabidopsis developed from overexpression of a somatic embryogenesis receptor-like kinase (SERK5-Ler), an LRR-RLK, had improved rooting depth [42]. Mutant rice with loss of function of an LRR-RLK-encoding gene, defective in outer cell layer specification 1 (DOCS1), had modified root gravitropism that resulted in a more open root cone angle (angle between most right and left of two external roots) and more radial root system, but reduced relative root elongation rate and shorter root hairs than the wild type [41]. Therefore, candidate genes encoding UGT and LRR-RLK could improve root traits in wheat, especially total RL and SRL. In addition, as the development of the 4AS NILs targeted drought tolerance QTL (Table 1), the genomic region might play a role in drought tolerance through improved root traits.
TraesCS4B01G114500 identified in qDSI.4B.1 encodes the Phox-associated domain. There are few reports on the role of this protein in controlling root traits in wheat. Lee et al. (2008) found that the homolog of the protein phosphatidylinositol 3-phosphate (PtdIns(3)P) regulated root hair growth in Arabidopsis [43]. The protein also controlled root growth and development and lateral root length in Arabidopsis [44].
TraesCS4B01G114800 identified in pairs 1-4 encodes bHLH-MYC and R2R3-MYB transcription factors N-terminal, an important protein for crop growth and development and abiotic stress tolerance in different crops, including wheat [45]. The regulatory genes encoding the proteins are poorly investigated in wheat [46]. MYC2 is essential for controlling root phenotype, and jasmonate-associated MYC2-LIKE (JAM) transcription factors inhibited root growth in Arabidopsis [47]. A cluster of MYC2, MYC3, and MYC4 played a role in root growth inhibition of Arabidopsis [48]. Wheat genes that encode MYC (TaMyc-B1) improved drought and salinity stress tolerance [49]. In wheat, apart from TaMyc1 and TaMyc-B1, MYC genes have not been explored [49]. Thus, there is an excellent scope to study the mechanism of TraesCS4B01G114800 in root trait control and drought stress tolerance in the future.
MYB transcription factor (TF) encoding genes in wheat, such as TaPHR3-A1 and TaSIM, improved P uptake, salt stress tolerance, grain number, and notably root length in Arabidopsis [35,50]. A mutant of soybean overexpressing wheat R2R3-MYB TF (GmMYB84) had greater primary rooting depth at optimum reactive oxygen species (ROS) level and drought tolerance [51]. Under low phosphate conditions, overexpression of the rice gene OsMYB2P-1 in Arabidopsis had greater lateral root density and longer primary and adventitious roots than the wild type [52]. In contrast, MYB TF AftMYB93 negatively regulated lateral root development and played a role in root hair growth in Arabidopsis [53]. The sesame MYB gene SiMYB75, a root-specific gene, promoted root growth and abiotic stress (drought and salt) tolerance in Arabidopsis and sesame [54]. In this study, it was shown that TraesCS4B01G114800 may play a role in RSR, RSA1, RSA2, and SRL, where the NIL isolines significantly differed in the traits. Table 4. Proteins and genes identified as responsible for root traits in crops and Arabidopsis.

TraesCS7A02G101000
Root growth, development, root length in rice, Arabidopsis; abiotic stress tolerance [68,69] Three of the candidate genes (TraesCS4A02G331000, TraesCS4A02G215800, and TraesCS4A02G215700) identified in pairs 5-7 and one gene (TraesCS7A02G431600) identified in pairs 9-10 encode F-box protein. Among the subunit proteins of Skp1-Cullin-F-box (SCF), F-box protein plays an important role in abiotic stress tolerance; in particular, its C-terminal structural domains interact with other proteins to identify substrates [56]. An F-box gene, TaFBA1, was identified in wheat organs, including roots [57]. In other few studies, overexpression of the gene showed drought tolerance, oxidative stress tolerance, heat tolerance, and longer root lengths relative to wild-genotype transgenic tobacco and Arabidopsis [56][57][58][59]. In rice, overexpression of MAIF1 improved abiotic stress tolerance and helped to control root growth through multiple signal pathways, including negative signaling with abscisic acid (ABA), auxin, and ABA interactions, and the sucrose signaling pathway [62]. In Arabidopsis, the F-box protein At1g08710 was found in roots, with longer roots under drought stress, and thus showed drought tolerance [60], whereas the F-box protein SKP2B (S-phase kinase-associated protein 2B) reduced lateral root formation by suppressing cell division in meristematic and founder cells [70]. In soybean, GmFBX176 increased root length to tolerate salt stress [61]. A meta-QTL analysis found that TraesCS7A01G481200 on chromosome 7A regulated root mass by encoding F-box protein [27].
TraesCS7A02G431500 identified in pairs 9-10 encode protein kinase. Alpha4 (α4) subunit of casein kinase 2, CK2 α4, a protein kinase regulated primary root elongation, lateral root formation, and root development in Arabidopsis [63,66]. Overexpression of stress/ABA-activated protein kinase 10 (SAPK10) in transgenic rice produced longer root hairs than the wild type [71]. In soybean, a root-specific serine/threonine kinase WNK (with no lysine K) homolog, GmWNK1, played an important role in root system architecture; overexpression of 35S-GmWNK1 in transgenic soybean produced shallower and fewer lateral roots than the wild type [67]. Furthermore, overexpression of GmWNK1 improved salt stress tolerance in Arabidopsis, possibly due to improved root traits [72]. In rice, phosphorus-starvation tolerance 1 (PSTOL1), a protein kinase gene, played an important role in improving root traits, including early growth of primary roots and increased rooting depth, root mass, and surface area; consequently, the gene improved N and P uptake [64]. In wheat, besides serine/threonine protein kinase genes, other protein kinase genes, such as TaSnRK2.4, TaCIPK29, and TaMPK4, affected root traits (including root length, root volume, and root surface area), grain yield, and abiotic stress tolerance [73][74][75]. Interestingly, an ortholog of an abiotic stress (drought and salt)-tolerant gene (OsMAPK5) in rice was found in wheat root tissue, indicating the importance of the protein-coding gene in abiotic stress tolerance [65]. Despite many protein kinases reports on different root traits in wheat, none has associated the protein with root mass. However, our study suggests that protein kinase-encoding genes contribute to genetic variations in root mass. Therefore, there is great scope for using the genes to improve root mass and other traits.
The seven NIL pairs significantly differed between the isolines for 11 important root traits. We identified 15 important genes on chromosomes 4BS, 4AS, and 7AL that encode proteins involved in controlling root traits. Of these, several have functions frequently reported as important proteins functioning in the controlling root traits of crops other than wheat, such as TraesCS4B01G114500 encoding the Phox-associated domain, TraesCS4A02G185300 and TraesCS4A02G442700 encoding UDP-glycosyltransferase (UGT), and TraesCS4A02G330900 encoding leucine-rich repeat receptor-like protein kinases (LRR-RLK). These candidate genes are important for further functional validation studies and marker developments to improve root traits in wheat.

Phenotyping for Root Traits
A randomized block design was used in a semi-hydroponic phenotyping system [76] for the 10 NIL pairs with three replications. The experiment was conducted in a temperaturecontrolled (10-24 • C) glasshouse at the University of Western Australia, Perth, from mid-June to late August 2019. The semi-hydroponic system included a 240 L (751 × 850 mm top, 1080 mm depth) mobile plastic bin, accommodating 16 growth units made from acrylic panels (260 mm wide × 480 mm long × 5 mm thick) wrapped in black calico cloth. The design, assembly, and operation of the semi-hydroponic system are described in Chen et al. (2011 [16,76].  (20) [16,76]. Wheat seeds were sterilized with 1% sodium hypochlorite (NaOCl) and germinated in multiple-welled plastic trays containing washed river sand (<2 mm). Seedlings with 4−5 cm long roots were washed in deionized (DI) water and transplanted into the bins. The plants in the growth units were watered with a controlled irrigation system described in   [16]. Based on the evapotranspiration rate calculated during the trial experiment, 7 L of nutrient solution was added to each bin weekly. The position of each bin was rotated weekly to minimize environmental effects.

Trait Measurements
All plants were harvested 42 days after transplanting at tiller onset (Zadoks 2.4) [77] when the deeper root systems reached the surface of the nutrient solution. Leaf number per plant (NLP) was counted the day before harvest. At harvest, each growth panel was removed from the bin and laid out on a flat bench. After removing the wrapped cloth, maximal shoot height (SH) and rooting depth (RD; maximum root depth of an individual plant measured from the crown) were measured with a ruler, and the number of nodal roots (NNR) counted. A portable photographing system, including a fluorescence light and Nikon D5100 camera, was set up above the bench [78] to photograph the root systems. After photographing, the roots were separated from the shoots at the crown, and the number of nodal roots counted. Shoots were cut from the crown and dried in an air-forced oven at 65 • C for 72 h to obtain shoot dry weight. Each root system was cut into 20 cm sections for optical scanning at 400 dpi using a desktop scanner (Epson Perfection V800/850). Root morphological traits (total root length (RL), root diameter (Rdia), root surface areas (RSAs), and root diameter length (RDCL) were determined for each root section after analyzing the root images in WinRHIZO Pro software (v2009, Regent Instruments Inc., Montreal, QC, Canada) [16]. The total root length (RL) was the sum of all root length sections from the same plant. Root diameter length and surface areas in two root diameter classes (0-0.25 mm and >0.25 mm) were measured. All root sections were combined for the same plant and dried in an air-forced oven at 65 • C for 72 h to obtain root dry mass (RM). Root/shoot ratio (RSR) was the ratio of RM and shoot dry mass. Specific root length (SRL) was the total root length per unit of biomass.

Statistics
A two-sample Mann-Whitney U-test was conducted using R package (R version 4.0.3) to detect significant differences between the isolines. A correlation analysis was done using RStudio Version 1.2.5033.

90K SNP Genotyping and Candidate Gene Identification
The wheat 90K Illumina iSelect array containing 81,587 SNPs was used to genotype genomic DNA samples of the NILs showing contrasting phenotypes [6]. SNP genotype calling and clustering were done with GenomeStudio 2.0 software (San Diego, California, United States) (Illumina, https://www.illumina.com, accessed on 2 March 2021). SNPs with no polymorphism, or those with more than 20% missing values, minor allele frequency (MAF) <0.05, or heterozygous calls >0.25 were deleted. The polymorphic SNP loci between each NIL pair were then analyzed to identify candidate genes associated with root traits.
The sequence of each polymorphic SNP was identified from the Triticeae Toolbox (T3) (https://triticeaetoolbox.org/wheat/, accessed on 2 March 2021) database and blasted with the wheat reference genome RefSeq v1.0 to locate the physical position [79]. To identify candidate genes, JBrowse (http://www.wheatgenome.org/Tools-and-Resources/ Sequences, accessed on 2 March 2021) databasewas used. Genes closely associated with the screening markers of the targeted QTL region and those genes located up to two mega base pairs (Mb) away from the SNPs [80] were scrutinized to identify candidate genes. The functions of the candidate high-confidence genes were searched on the International Wheat Genome Sequencing Consortium (IWGSC) RefSeq v1.0 website (https://wheaturgi.versailles.inra.fr/Seq-Repository/Annotation, accessed on 2 March 2021) [15]. Genes previously related to the growth and development of root traits in plants were considered potential candidate genes.
Supplementary Materials: The following are available online at https://www.mdpi.com/article/10 .3390/ijms22073579/s1: Supplementary Table S1. All replicated data, mean values, and standard deviations of different root and shoot traits of wheat near-isogenic lines used in this study. Supplementary Table S2. Genes within and near the markers flanking QTL qDSI.4B.1, QTL qDT.4A.1, and QHtscc.ksu-7A on chromosomes 4BS, 4AS, and 7AL of near-isogenic lines (NILs) of hexaploid wheat. Figure S1. Differences in root traits in the pairs of near-isogenic lines (NILs) grown in a semi-hydroponic system under controlled conditions: (a) total root length (m), (b) root mass (mg); "*" indicates significant differences between isolines (NIL isoline a and NIL isoline b) at p < 0.05. The bars indicate the standard errors. Figure S2. Correlogram across phenotypic root traits. Correlated traits are corresponded to significant correlation test (Pearson's correlation at 5% level of significance); blue and light salmon colors indicate positive and negative correlation, respectively. The color intensity indicates the strength of the correlation between two traits: darker boxes are corresponded to stronger correlation (i.e., closer to -1 or 1), whereas light-colored boxes indicate weaker correlation. RL = total root length (cm), RM = root dry mass (g), RD = rooting depth (cm), RSR = root/shoot ratio, RSA1 = root surface area (>0. 25