QTL Mapping and Favorable Allele Mining of Nitrogen Deficiency Tolerance Using an Interconnected Breeding Population in Rice

Nitrogen is one of the most important nutrients for rice growth and development. Breeding of nitrogen deficiency tolerance (NDT) variety is considered to be the most economic measure to solve the constrain of low nitrogen stress on grain yield in rice. An interconnected breeding (IB) population of 497 lines developed using Huanghuazhan (HHZ) as the recurrent parent and eight elite lines as the donor parents were tested for five traits including grain yield, biomass, harvest index, thousand grain weight, and spikelet fertility under two nitrogen treatments in three growing seasons. Association analysis using 7,388 bins generated by sequencing identified a total of 14, 14, and 12 QTLs for the five traits under low nitrogen (LN), normal nitrogen (NN), and LN/NN conditions, respectively, across three seasons. Favorable alleles were dissected for the 40 QTLs at the 10 NDT regions, and OM1723 was considered as the most important parent with the highest frequency of favorable alleles contributing to NDT-related traits. Six superior lines all showed significantly higher GY in LN environments and similar GY under NN environments except for H10. Substitution mapping using near-isogenic introgression lines delimited the qTGW2-1, which was identified on chromosome 2 under LN, NN, and LN/NN conditions into two QTLs, which were located in the two regions of about 200 and 350 kb with different favorable alleles. The bins 16, 1301, 1465, 1486, 3464, and 6249 harbored the QTLs for NDT detected in this study, and the QTLs/genes previously identified for NDT or nitrogen use efficiency (NUE) could be used for enhancing NDT and NUE by marker-assisted selection (MAS).


INTRODUCTION
Rice (Oryza sativa L.) is one of the most important staple crops in the world and also the main calorie source for more than 65% of the population in China. Nitrogen is one of the most important nutrients in the growth stage of crops (Yoshida, 1981). In China, about 70% of the rice paddy fields have low or moderate nitrogen supply, which limits yield potential of rice production (Feng et al., 2018). Therefore, development of rice varieties with good nitrogen deficiency tolerance (NDT) is considered as a key method of sustainable agriculture for food security.
In recent years, many QTLs/genes for NDT, measured as the relative trait values under low nitrogen (LN) stress to normal nitrogen (NN) conditions have been reported by QTL mapping and gene expression methods. QTLs were detected for the traits or the relative trait values of seedling height, shoot dry weight, chlorophyll content (Tong et al., 2011), maximum root length, root dry weight, plant dry weight (Lv et al., 2010;Zhao et al., 2014) at seedling stage, spikelet fertility percentage (Shan et al., 2005), grain yield (Tong et al., 2006), panicle number and total panicle weight , grain number per panicle, spikelet fertility percentage and 100-grain weight (Tong et al., 2011), and nitrogen efficiency (Tang et al., 2011) at maturity stage using populations such as backcross, chromosome segment substitution line, and recombinant inbred lines under LN stress or different nitrogen level conditions. A total of 14 QTLs for NDT-related traits (relative shoot and root biomass, and relative plant height) were identified by Lian et al. (2005) in an RIL population derived from Zhenshan97 × Minghui63. Fifteen QTLs were identified for the four NDT traits including relative grain yield, relative biomass, relative grain nitrogen, and relative biomass nitrogen (Wei et al., 2012), and some of which are close to genes controlling nitrogen cycle.
Most of NDT QTLs reported in rice were conducted using populations derived from two parents and sparse linkage maps constructed using restriction fragment length polymorphism (RFLP) or simple sequence repeat (SSR). It is very hard to obtain precise information about the QTLs using a small number of markers, which were coarsely located. Fine mapping using a large secondary population and new markers are needed (Wang et al., 2014). With the development of sequencing technology, highdensity single-nucleotide polymorphism (SNP) markers can be easily and quickly generated, which has been widely applied to genome-wide association studies (GWAS) and QTL mapping in rice and many other crops (Yan et al., 2010;Wang et al., 2014;Chen et al., 2020).
The past decade has seen the rise of multiparental populations as a study design offering great advantages for genetic studies in plants. Multiparent mapping populations such as nested association mapping (NAM) population (Yu et al., 2008) and multiparent advanced generation inter-cross (MAGIC) population (Cavanagh et al., 2008) have been developed for many crops. A NAM population is usually generated by crossing multiple genotypes with a single genotype. Using multiparental population for QTL mapping provides an opportunity to test pleiotropy, genetic background effect, and the genetic overlap between different complex traits (Qu et al., 2020). Buckler et al. (2009) developed a NAM population that was consisted of 25 families with 200 lines per family in maize. The NAM population has been used for studying the complex traits in maize, such as pathogen resistances (Kump et al., 2011), morphological traits (Tian et al., 2015), and kernel composition (Cook et al., 2012). Jordan et al. (2011) developed a NAM population containing 56 families for detecting QTLs in sorghum. However, the application of NAM population for QTL mapping in rice is rarely reported.
A backcross (BC)-based breeding strategy has been adopted by our team to improve multiple abiotic stress tolerance for many years (Ali et al., 2006;He et al., 2010;Meng et al., 2013;Wang et al., 2013;Feng et al., 2018). A few outstanding varieties were used as recurrent parents, and many varieties and landraces were used as donor parents in combination with a selection of tolerance to multiple stresses and grain yield. As a result, many small-to-medium size populations sharing a common parent were developed (Ali et al., 2017). Mapping QTL for a range of complex traits using these populations has been carried out Cui et al., 2015;Feng et al., 2018). However, the interconnectedness between the populations has not been well exploited. A set of breeding populations linked together by a common parent can be regarded as a single interconnected breeding (IB) population for mapping purpose, which is similar to a NAM population that consists of subpopulations sharing a common parent. It is expected that IB population could significantly increase the mapping resolution and power by exploring multiple populations simultaneously. An IB population consisting of highly selected introgression lines derived using HHZ as recurrent parent and eight elite lines as donors has been successfully used in identifying QTL for cold tolerance at the booting stage (Zhu et al., 2015). The fact that a QTL for cold tolerance was fine mapped to a 192-kb region encouraged us to explore this IB population for other traits.
The current study was to identify QTL and favorable alleles for NDT at the reproductive stage. The objectives of this study were (1) to identify the QTLs affecting NDT under different nitrogen conditions; and (2) to gain a better understanding of the genetic relationships between NDT and GY at the QTL level. The results will provide a good example of QTL mapping using breeding population and useful information for rice breeding of NDT by marker-assisted selection (MAS).

Plant Materials and Field Experiments
The HHZ IB population used in the study was derived from recurrent parent Huanghuazhan (HHZ) crossed with eight donor parents IR50, IR64, Teqing, PSBRc28, PSBRc66, CDR22, OM1723, and Phalguna. Eight F 1 s were backcrossed once with HHZ to harvest the BC 1 F 1 . Through multiple generations selfing and screenings for yield and abiotic stresses (drought, salinity, and submergence), a total of 496 BC 1 F 5 lines were produced (Ali et al., 2017). A total of 496 lines of the IB population plus the HHZ were field tested under different nitrogen conditions in the early and late seasons of 2013 (2013E and 2013L) and the early season of 2014 (2014E) in Shenzhen (22 • 52 N 113 • 46 E), Guangdong province, China. Two nitrogen conditions were applied, they were (1) low nitrogen condition (LN), no chemical nitrogen fertilizer was applied to the LN paddy since 2012; (2) normal nitrogen condition (NN), 150 kg N ha −1 was applied with two splits. About 70% nitrogen fertilizers were used as basal, and 30% nitrogen fertilizers were applied at 15 days after transplanting under the normal nitrogen condition. Phosphorus (40 kg ha −1 ) and potassium (40 kg ha −1 ) were applied as basal under two nitrogen conditions. All fertilizers were broadcast by hand. The LN paddy soil had the following properties: pH 6.18, organic matter of 4.32 g kg −1 , total N of 380 mg kg −1 , available P of 59.2 mg kg −1 , and available K of 165 mg kg −1 . The key properties of the NN paddy soil were pH 5.97, organic matter of 10.9 g kg −1 , total N of 920 mg kg −1 , available P of 88.8 mg kg −1 , and available K of 155 mg kg −1 .
The experiments were conducted using a randomized block design with two replications. The seeds were sowed on March 5 in 2013E and 2014E, and on July 13 in 2013L. The seedlings were transplanted into three rows with eight plants per row at a planting density of 20 cm × 16.5 cm on April 5 in 2013E and 2014E, and on August 5 in 2013L. The field management was according to the local rice production, and the pests, diseases, birds, weeds, and rats were intensively controlled during the whole growth duration.

Sampling and Trait Measurement
At maturity stage, six uniform plants from the middle of each plot were harvested. The plants were cut at ground level, and grains were separated from other parts (leaves and sheath, and stems and panicle branches). All samples were oven-dried at 80 • C until a constant weight was achieved. Grain yield per plant (GY), biomass per plant (BM), thousand grain weight (TGW), and spikelet fertility (SF) were measured. Harvest index (HI) was calculated by GY/BM. Also, specific parameters were calculated using the following equations:

SNP Genotyping and Bin Markers
To get the consensus sequence that can be used to compare the differential bases between the HHZ and the donor parent, all reads from each parent with 30× re-sequencing were first aligned to the Nipponbare reference genome. Based on the differential SNP results, the pseudo-molecules of both parents were generated in pairs by using Perl script for eight mapping populations, respectively. A total of 400 k high-quality SNPs was developed based on whole-genome re-sequencing of 2× in IB population. To improve analytical accuracy and the resolution of re-sequencing data, we used the sliding window approach to determine the recombination breakpoint based on the SNP ratios (Huang et al., 2009). Based on the recombination breakpoint, the 496 lines were aligned to all the chromosomes and compared for the minimal of 50-kb intervals. Adjacent 50-kb intervals with the same genotype across the entire IB population were recognized as a single recombination bin. A total of 7,388 bins covered in all of 12 chromosomes were obtained for association analysis.

Statistical and Association Analyses
We used the R project for statistical analysis. The basic statistics of phenotypic values were calculated using the R package named EnvStats. Correlation coefficients were obtained from the R function called corr. Two-way ANOVA was conducted by the R function called ANOVA. Based on the adjusted means and the trait ratios of LN to NN of ILs, association analyses were conducted by R-package MAGICqtl 1 . This package implemented random model-based methods proposed by Wei and Xu (2016) for QTL detection using a multi-parental population. The initial threshold to declare a significant association was set at p = 1.0 × 10 −4 . The threshold value of the main effect QTL was LOD > 3.2. The gene effects of multiple alleles in each QTL were estimated by the difference between the mean phenotypic value of one allele and mean phenotypic value of all the alleles.

Phenotypic Performance
Grain yield and BM were much higher in the NN than in the LN, while HI and SF were similar in the two nitrogen treatments ( Table 1). TGW was slightly higher in the NN. The average GY, BM, HI, SF, and TGW of the IB population were similar to those of HHZ. The estimates of skewness ranged from −2.873 to 1.783, and the estimates of kurtosis ranged from 0.059 to 16.65. The heritability in the NN was similar to or higher than that in the LN in all three seasons, ranging from 0.647 of BM under LN in 2013E to 0.929 of HI under NN in 2014E.
Two-way ANOVA (genotype and season) indicated that season was highly significant for GY, BM, HI, TGW, and SF in the two nitrogen treatments (Supplementary Table 1). Genotype was highly significant for GY, HI, TGW, and SF in the two nitrogen treatments. Genotype-by-season interaction was highly significant for HI and SF in the two nitrogen treatments. For the five traits, the main source of the variation was the season. The effect of genotype was significant for all five traits in the two conditions except for BM in the NN treatment.
In the LN, GY was highly, significantly, and positively correlated with other traits in three seasons except TGW in 2014E (Supplementary Table 2). In the NN, GY was highly significantly correlated with BM, HI, and SF, and not significantly correlated with TGW in the three seasons.

QTL Mapping
A total of 40 QTLs were identified for the five NDT-related traits across the six environments (nitrogen-season combinations) (Table 2 and Figure 1).

QTLs Detected for Traits in Low Nitrogen (LN) Condition
One QTL for GY, qGY1 was detected on chromosome 1 in 2013L and explained 15.0% of the phenotypic variation ( Table 2). For BM, three QTLs (qBM1, qBM2-1, and qBM2-2) were identified on chromosomes 1, 2, and 2 in 2013L and 2013E, and accounted for phenotypic variations of 9.7%, 3.3%, and 8.2%, respectively. One QTL named qHI8-2 was identified for HI, which explained a phenotypic variation of 7.9%. Three QTLs were identified for TGW. Specifically, qTGW2-1 was simultaneously detected on chromosome 2 in all three conditions with phenotypic variations from 8.3 to 13.4%. qTGW2-2 was detected on chromosome 2 in 2013L and 2014E with phenotypic variations of 14.2% and 10.7%, while qTGW10 was detected on chromosome 10 in 2013L and 2014E with phenotypic variations of 1.9-3.5%. Two QTLs (qSF2, qSF8) were mapped to chromosomes 2 and 8 for SF in 2013L and 2013E, which accounted for phenotypic variations of 6.1% and 18.4%, respectively.

QTLs Detected for Traits in Normal Nitrogen (NN) Condition
One QTL named qBM1 was detected on chromosome 1 for BM in 2014E (Table 2), which explained 8.3% of the phenotypic variation. For HI, three QTLs (qHI1, qHI6-1, and qHI8-1) were identified. The QTL on chromosome 1, qHI1, was identified in 2013E and explained 11.6% of the phenotypic variation. The QTL, qHI6-1, on chromosome 6 was simultaneously identified in 2013E and 2014E and explained 20.1% and 19.0% of the phenotypic variation. The QTL qHI8-1 on chromosome 8 was detected in 2013E and explained 13.5% of the phenotypic variation. Four QTLs were identified for TGW. qTGW2-2 was simultaneously identified on chromosome 2 in all three conditions and explained 6.5-9.1% of phenotypic variation. qTGW2-1 and qTGW10 were mapped to chromosomes 2 and 10 in 2014E and explained 6.4% and 2.0% of the phenotypic variation, respectively. qTGW7 was detected on chromosome 7 in 2013E with a phenotypic variation of 5.0%. For SF, two QTLs (qSF2 and qSF5) were mapped on chromosomes 2 and 5 in 2013L with phenotypic variations of 8.5% and 12.7%, respectively, while another one (qSF8-1) was mapped on chromosome 8 in 2013E with phenotypic variations of 13.9%.

Favorable Allele Mining of Important QTLs
Ten regions were detected for NDT-related traits in multiple seasons (Table 3). Among them, five were reported near the cloned genes related to nitrogen use efficiency (NUE) and/or functional ammonium transporter. The region bin16 on chromosome 1 was detected and associated with multiple traits, such as GY, RGY, BM, HI, and RHI, and the best favorable alleles all came from the donor parent OM1723 except the HI in 2013E NN environment. The region bin1301 was identified in relation to the traits BM, RBM, and SF, which were located very close to the cloned gene LOC_Os02g38230 with functional annotation of partner protein for high-affinity nitrate transport . The best favorable alleles for BM and RBM came from the donor Teqing, while the best favorable alleles for SF came from OM1723 and Phalguna in 2013L LN and 2013L NN environments, respectively. Even though the best favorable alleles came from different donors, they are all carrying higher positive effects than the recurrent parent HHZ. The bin1465 and bin1486 were very near the cloned genes LOC_Os02g47280 or GRF4 (Hu et al., 2015;Li et al., 2018a)

Selecting Super Nitrogen Deficiency Tolerant Lines
Based on the favorable allele dissection, six superior lines with four favorable introgressed alleles were selected (       NN. b Data of GY (grain yield per plant) and TGW (1000-grain weight) were calculated based on the three seasons. * and ** represent significant levels at P < 0.05 and 0.01, respectively. c The environment where the bins were detected; the parents in the brackets are those whose favorable alleles came from.

Fine Mapping of qTGW2-1
The qTGW2-1, which was simultaneously detected in 2014E in three conditions (LN, NN, and LN/NN) and also in 2013E and 2013L in LN conditions, most likely associated with NDT and the favorable alleles mainly came from OM1723 and PSBRc28 in different conditions. To fine map the stable qTGW2-1 in the region of bin1456-bin1466 on chromosome 2 by means of substitution mapping, ILs with overlapped introgressed segments and covered the genomic region of qTGW2-1 were selected from the IB population. Seven selected ILs had the recurrent genome percentage higher than 94.0%. The graphic genotypes of these seven ILs are given in Figure 2. Three, three, and one lines had the alleles from the donor parents PSBRc66, IR50, and Phalguna, respectively. The lines of H129, H342, H358, and H359 had significantly higher TGW than the recipient parent HHZ, while H180 and H181 had significantly lower TGW than the recipient parent. The TGW of H165 was very similar to that of HHZ. Therefore, the qTGW2-1 was separated into two QTLs for TGW, one was located in about the 200-kb region flanked by bin1456 and bin1460 with the positive alleles for increased TGW from the donor parent PSBRc66 and Phalguna, and the other one was located in about the 350-kb region delimited by bin1459 and 1466 with the negative allele for decreased TGW from the donor parent IR50. The GRF4 for NUE (Hu et al., 2015;Li et al., 2018a) was around 633-kb apart from the QTL flanked by bin1456 and bin1460 (Figure 2), suggesting the two QTLs separated from qTGW2-1 were new loci for NDT.

Comparisons of the NDT QTLs Identified in This Study With Previously Reported QTLs and Cloned Genes
Nitrogen-deficiency tolerance (NDT) traits were measured by the ratios of the trait values under LN to those under NN and could be used as a criterion for selecting genotypes for the tolerance of low nitrogen. Many agronomic traits such as effective tiller number, biomass, and grain yield have been studied under low-N stress in rice, and some QTLs associated with N translocation were found to locate in the same chromosomal regions (Shan et al., 2005;Tong et al., 2006;Senapathy et al., 2008). Among 10 QTL regions identified for the NDT-related traits in this study (Table 3), the region bin16 on chromosome 1 harboring QTLs for BM, GY, and HI was associated with the QTLs affecting shoot dry weight and plant dry weight at seedling stage under low nitrogen condition identified by Lian et al. (2005), and the QTLs for relative shoot weight, relative plant height, and relative biomass yield at seedling stage (Feng et al., 2010). QTLs affecting BM and SF in the bin1301 identified in this study were mapped together with the QTLs for dry shoot weight under low nitrogen condition (Lian et al., 2005), grain nitrogen use efficiency and grain yield (Wei et al., 2011), soluble protein content detected in the backcross inbred lines of Nipponbare/Kasalath (Ishimaru et al., 2001), and the cloned gene LOC_Os02g38230 for partner protein for high-affinity nitrate transport . The region bin1465 harboring QTLs for BM and for TGW on chromosome 2 identified in the present study were mapped together with the QTL affecting root dry weight under low nitrogen condition (Lian et al., 2005) and near the LOC_Os02g47280 for nitrogen use efficiency and growthregulating factor (Hu et al., 2015;Li et al., 2018a), and the gy2b affecting GY under low nitrogen condition (Cho et al., 2007). The bin1486 harboring QTLs for TGW and RTGW on chromosome 2 in this study were mapped together to the QTL for the ratio of plant height under low nitrogen stress to normal conditions (Feng et al., 2010) and the LOC_Os02g53130, a nitrate reductase gene for nitrogen use efficiency (Gao et al., 2019). The bin3464 harboring QTLs for RSF and SF on chromosome 5 in this study was mapped near the qHGW-5b for grain weight under low nitrogen input (Tong et al., 2011) and the LOC_Os05g39240 for functional ammonium transporter, which was constitutively expressed in roots and shoots (Gaur et al., 2012). The bin4899 harboring QTLs for HI and SF on chromosome 8 in this study was associated with the QTL for plant dry weight under low nitrogen condition (Lian et al., 2005). The bin6249 on chromosome 10 affecting TGW detected under LN, NN, and LN/NN conditions in this study harbored the QTL (qFGPP-10b) for filled grains per panicle under low nitrogen input (Tong et al., 2011), the qNR10 for nitrogen response (Wei et al., 2012), and the LOC_Os10g40600, a nitrate-transporter gene for nitrogen use efficiency Zhang et al., 2019). Allelisms of the above QTLs for NDT-related traits identified in this study with the previously reported QTLs or cloned genes for NDT or NUE require further verification via fine mapping and QTL cloning. Nitrogen deficiency tolerance is defined as the plants' ability to maintain normal growth and good yield when the soil N content is low. Thus, NDT could have something to do with NUE to some extent. Based on comparisons of QTLs affecting NDT detected in this study with those previously reported, five regions underlying NDT at bins 1301, 1465, 1486, 3464, and 6249 on chromosomes 2, 2, 2, 5, and 10, respectively, were found to map together with or near the five cloned NUE-related genes. Using 127 recombinant inbred lines (RILs) derived from the cross of Zhenshan97/Minghui63, Wei et al. (2012) also identified four genomic regions containing QTLs for NDT and NUE traits with the same additive effect directions. It was suggested that partial genetic overlaps at least exist between NDT and NUE, which is beneficial for rice breeding for simultaneously improving NDT and NUE.

Mining Favorable Alleles for NDT
From the breeding point of view, the major use of identified QTL or cloned gene is to mine favorable alleles for efficient molecular breeding. So far, over 2,000 genes controlling important agronomic traits have been cloned in rice (Li et al., 2018b). Unfortunately, the phenotypic effects of different alleles in rice germplasm accessions at these cloned rice genes and their values in rice breeding remain largely unknown. This is one of the most important reasons why so many cloned genes have never been exploited in breeding through MAS so far. Using eight BC1 populations derived from a widely adaptable recipient and eight donors plus three rounds of phenotypic selection, 496 trait-specific introgression lines (ILs) in Huanghuazhan (HHZ) background have been developed (Ali et al., 2017), forming a material platform for discovery of QTLs underlying the target and non-target traits (Zhu et al., 2015;Feng et al., 2018). The IL population has characteristics of similar developmental stage and plant height in elite rice background, thus very suitable for evaluation of abiotic stress tolerances such as drought, salt, and cold whose tolerance performances largely depend on developmental stage. In this study, favorable alleles at QTLs for different NDT-related traits in the 10 NDT-related regions were dissected ( Table 3). Among the total 40 loci affecting different NDT-related traits at the 10 bins under different conditions (

Implications for Breeding
Historically, many breeding programs took yield potential as a primary target, particularly in China in the past few decades. In China, pursuing high yield of super rice cultivar has been accompanied with excess uses of chemical fertilizers and pesticides in rice production, resulting in serious environmental pollutions and largely reduced rate of fertilizer utilization by the crops (Zhang, 2007). Nevertheless, modern semidwarf rice cultivars have rarely achieved their yield potentials in farmers' fields because of many abiotic and biotic stresses. The main reason for this huge yield gap is primarily due to the fact that most super rice cultivars developed under the high input conditions do not perform well in more than 70% of the moderate-and low-yielding fields in China, much of which suffer more frequently inadequate fertilizers and other abiotic stresses (Pandey, 1997;Feng et al., 2018). Consequently, developing crops that are less dependent on the heavy application of N fertilizers is essential for the sustainability of agriculture. NDT traits have been considered as indirect selection criteria for the improvement of NUE (Lian et al., 2005;Namai et al., 2009;Feng et al., 2010;Wei et al., 2012). As a matter of experience, the lines with superior yield and yield-related traits under both LN and NN conditions, or superior under LN and normal or average under NN, could be preferentially selected because such lines can use the nitrogen efficiently to produce BM and GY under limited nitrogen supply (Lafitte and Edmeades, 1994).
Nitrogen deficiency tolerance and NUE are complex quantitative traits controlled by multiple genes, which makes it difficult to improve these complex traits using conventional breeding approaches. Recently, high-throughput SNP genotyping based on re-sequencing and gene chips promises to greatly accelerate QTL mapping and pyramiding on the whole genome (Thomson, 2014;Feng et al., 2018). Ten important NDT regions were identified in the breeding population with elite variety background in this study. Some of them were co-located with the previously reported NUE-related QTLs or genes, which are useful for rice breeding for high NDT or both for NDT and NUE by introgressing or pyramiding of favorable alleles at those important loci by MAS. Specifically, the Bin16 affecting GY, RGY, BM, and RHI under different nitrogen conditions, and other five bins 1301, 1465, 1486, 3464, and 6249 shared with the previously reported QTLs, or genes underlying NDT-and NUE-related traits could be used in MAS-based or QTL-designed breeding for enhancing NDT and/or NUE in rice. As indicated in Table 4, two promising NDT lines, H10 carrying favorable alleles from OM1723 at bins 16 and 1465, and H11 with favorable alleles from OM1723 at bins 16, 1465, and 1486, can be used as donor parents to improve NDT of an elite variety by marker-assisted selection against the flanking markers linked to respective bins. Of course, important NDT loci identified herein could be also pyramided with other previously identified QTLs or genes controlling NDT and NUE as most of the NDT and NUE QTLs identified previously under different environmental conditions were genetically independent (Lian et al., 2005;Senapathy et al., 2008).

CONCLUSION
A total of 14, 14, and 12 QTLs for the five traits was identified under LN, NN, and LN/NN conditions, respectively, in an interconnected breeding population across three seasons. Among them, 10 NDT-bin regions were identified, and the favorable alleles contributing to NDT-related traits were primarily from OM1723, secondarily from Phalguna. Six superior lines showed significantly higher GY in LN environments and similar GY under NN environments due to introgressing the favorable alleles at the related NDT-QTLs. The bins 16,1301,1465,1486,3464, and 6249 shared the QTLs or genes for NDT identified in this study and NUE previously reported could be used for improvement of NDT and NUE by MAS.

DATA AVAILABILITY STATEMENT
All datasets presented in this study are included in the article/Supplementary Material.