ORIGINAL RESEARCH article
Genetic Dissection of Root Morphological Traits Related to Nitrogen Use Efficiency in Brassica napus L. under Two Contrasting Nitrogen Conditions
- Key Laboratory of Biology and Genetic Improvement of Oil Crops, Ministry of Agriculture, Oil Crops Research Institute of the Chinese Academy of Agricultural Sciences, Wuhan, China
As the major determinant for nutrient uptake, root system architecture (RSA) has a massive impact on nitrogen use efficiency (NUE). However, little is known the molecular control of RSA as related to NUE in rapeseed. Here, a rapeseed recombinant inbred line population (BnaZNRIL) was used to investigate root morphology (RM, an important component for RSA) and NUE-related traits under high-nitrogen (HN) and low-nitrogen (LN) conditions by hydroponics. Data analysis suggested that RM-related traits, particularly root size had significantly phenotypic correlations with plant dry biomass and N uptake irrespective of N levels, but no or little correlation with N utilization efficiency (NUtE), providing the potential to identify QTLs with pleiotropy or specificity for RM- and NUE-related traits. A total of 129 QTLs (including 23 stable QTLs, which were repeatedly detected at least two environments or different N levels) were identified and 83 of them were integrated into 22 pleiotropic QTL clusters. Five RM-NUE, ten RM-specific and three NUE-specific QTL clusters with same directions of additive-effect implied two NUE-improving approaches (RM-based and N utilization-based directly) and provided valuable genomic regions for NUE improvement in rapeseed. Importantly, all of four major QTLs and most of stable QTLs (20 out of 23) detected here were related to RM traits under HN and/or LN levels, suggested that regulating RM to improve NUE would be more feasible than regulating N efficiency directly. These results provided the promising genomic regions for marker-assisted selection on RM-based NUE improvement in rapeseed.
As the key element of protein, nucleic acid etc., nitrogen (N) is one of the essential mineral nutrients for crop growth and development. Application of sufficient synthetic N fertilizer at the appropriate time can significantly improve crop yield and quality, but only 30–50% of the applied N fertilizer is taken up from soil by crops (Smil, 1999; Cassman et al., 2002). The overuse of N fertilizer globally causes an increased cost as well as serious environmental problems, such as soil acidification, groundwater and air pollution, and so on (Zhu and Chen, 2002; Good et al., 2004; Galloway et al., 2008; Wuebbles, 2009; Guo et al., 2010). Improving the N use efficiency (NUE) and reducing residual N in environment is therefore the longstanding, worldwide challenge for economical and sustainable agricultural production.
In general, NUE referring to plant productivity can be divided into two components: N uptake efficiency and N utilization efficiency (Good et al., 2004; Chardon et al., 2010). Root system is the key place for nutrient and water acquisition and its ability on soil exploring is the major determinant of N uptake efficiency (Li X. et al., 2016). Root system architecture (RSA), often defined as the spatial configuration and distribution of root system in the growth medium, determines the soil exploration in time and space (Lynch, 1995). Considerable researches have shown that RSA was closely related with N uptake, for example, plants with a steeper and deeper root can absorbs N more efficiently in deep soil layers (Lynch, 2013; Trachsel et al., 2013; Zhan and Lynch, 2015). As an important part of RSA, root morphology (RM, including root length, root surface area, root number etc.) plays an important role in crop N acquisition under field or artificial growth medium condition (Liu et al., 2008; Li et al., 2015; Mu et al., 2015). Correspondingly, RM is in turn significantly influenced by N availability. Therefore, heterogeneity of N supply permits plants to optimize N acquisition from growth medium by modulating their root morphology (Linkohr et al., 2002; Walch-Liu et al., 2006; Yu et al., 2014; Li P. et al., 2016).
Both NUE and root system are complex traits and cannot be measured directly; especially the hidden feature of root makes it difficult to be investigated. Thus, assessment the features of NUE and root system under field conditions for conventional genetic improvement is slow, imprecise and expensive (Hochholdinger and Tuberosa, 2009). In recent decades, various glasshouse- or laboratory-based new technologies (such as hydroponic culture, agar plate culture and paper culture system) have been developed and applied for seedlings root system studying (Yang et al., 2010; Shi et al., 2013; Li P. et al., 2016). Furthermore, quantitative trait locus (QTL) analysis based on genome-wide molecular linkage mapping has become a powerful approach for dissecting the genetic basis of root system as related to nutrient-deficiency tolerance, including several RM-related traits such as lateral root density, total/mean lateral root length, crown/seminal root number, etc. (Zhang et al., 2014, 2016; Li et al., 2015; Gu et al., 2016). With the easiness of phenotype identification, those glasshouse- or laboratory-based new technologies were successfully applied for QTL mapping of root traits and/or NUE traits (Coque et al., 2008; Li et al., 2015). For example, many studies on QTL analysis for RM and/or NUE under low-nitrogen (LN) and high-nitrogen (HN) conditions have been reported in rice (Huang et al., 2004; Lian et al., 2005), maize (Coque et al., 2008; Liu et al., 2008; Li et al., 2015; Li P. et al., 2016; Pestsova et al., 2016), wheat (An et al., 2006), barley (Hoffmann et al., 2012), and cotton (Shang et al., 2016). Hydroponic culture system was used for most of the above studies, while agar plate culture system (Huang et al., 2004) and paper culture system (Pestsova et al., 2016) were applied for the rest two studies. Partial of these studies uncovered a significantly genetic relationship between RM and NUE traits and provided the most promising genomic regions for marker-assisted selection of RM to improve NUE. Importantly, Li et al. (2015) identified 53 advanced backcross-derived lines (ABLs) containing RSA-NUE QTL clusters (root traits included seminal root length, crown/seminal root number, etc.) via marker-assisted selection and the grain yield (GY)/NUE of these ABLs showed apparent mean increases of 13.8% under HN and of 15.9% under LN conditions, compared to recurrent background Wu312, providing a successful study case for the manipulation of RM to improve NUE via marker-assisted selection QTLs in maize.
Rapeseed (Brassica napus L., AACC, 2n = 38) is one of the most important oil crops in the world, whereas N fertilization is a limiting factor in rapeseed productivity (Rathke et al., 2005). Understanding the genetic control of RM and improving N uptake ability in low N environment is believed to be an effective way to improve NUE and maintain sustainable production in rapeseed. Although, several root traits QTLs as related to phosphorus uptake efficiency (Yang et al., 2010; Shi et al., 2013; Zhang et al., 2016) and QTLs of NUE traits at agronomic and physiological levels (Bouchet et al., 2014, 2016) have recently been identified, few QTLs of root morphological traits responsible for NUE have been reported in rapeseed. These hinder the progress in nitrogen-efficient genetic improvement through marker-assisted selection of effective RM in rapeseed. In the present study, a rapeseed recombinant inbred line (RIL) population (BnaZNRIL) was used to investigate RM- and NUE-related traits of seedling plants in hydroponics under two contrasting N levels. The main objectives of this study were as follows: (i) to investigate the genetic relationship between RM and NUE under HN and LN conditions; (ii) to identify QTL clusters with pleiotropy or specificity for RM- and NUE-related traits through QTL mapping and QTL meta-analysis; (iii) to provide the promising genomic regions for marker-assisted selection on NUE improvement in rapeseed. Further, our study proposed two distinct NUE-improving approaches (RM-based and N utilization-based directly) in rapeseed, and indicated that RM-based NUE improvement would be more effective than NUE improvement directly via marker-assisted selection in rapeseed.
Materials and Methods
The BnaZNRIL population, consisting of 184 F7 lines, was derived by single-seed descent from a cross between two sequenced rapeseed cultivars, Zhongshuang11 (de novo sequencing) and No. 73290 (resequencing) (Yang et al., 2016). As 9 lines were removed from this population for the lack of seeds, 175 lines plus the two parents were used in this study.
Plump and uniform rapeseed seeds were sowed on medical gauze that was fixed to a blue plastic basin (60 × 40 × 15 cm, length × width × height) filled with quarter-strength modified Hoagland's solution (Hoagland and Arnon, 1950). The modified Hoagland's solution (the concentration of N was 15 mM) consisted of: 5 mM Ca(NO3)2·4H2O, 5 mM KNO3, 2 mM MgSO4·7H2O, 1 mM KH2PO4, 0.05 mM EDTA-Fe, 46 μM H3BO3, 9.14 μM MnCl2·4H2O, 0.77 μM ZnSO4·7H2O, 0.37 μM NaMoO4·2H2O, and 0.32 μM CuSO4·5H2O. Six days after sowing as described by Dun et al. (2016), uniform seedlings were selected and transplanted into smaller blue plastic basins (34 × 26 × 12 cm, length × width × height) containing quarter-strength nutrient solution (two N treatments, HN and LN) under the natural condition with a removable rain-shelter. Each basin contained 24 seedlings of 4 lines (six seedlings for each line). As a total, 90 basins were used in each independent experiment. Nutrient solution was renewed once a week. For HN treatment, the quarter-strength and half-strength nutrient solution was used at the first 2 weeks respectively, and full-strength nutrient solution was used until harvest. For LN treatment, the concentration of N was adjusted to 0.5 mM by reducing KNO3 and replacing Ca(NO3)2 by CaCl2, while K+ was complemented by adding K2SO4 (refer to Stahl, 2015). The pH value was adjusted to 5.8 ± 0.2 with NaOH or HCl.
Three independent hydroponic culture experiments with a completely random design were carried out at Oil Crops Research Institute of the Chinese Academy of Agricultural Sciences, Wuhan, PR China. The plants were harvested with five fully expanding leaves under HN condition (~37 days after sowing). Accordingly, the plants under LN condition showed typical N-deficiency symptoms. The first experiment (E1) was carried out during the period from October 18 to November 23 in 2015 (37 days), the second (E2) from March 7 to April 14 in 2016 (39 days), and the third (E3) from October 2 to November 4 in 2016 (34 days). During the three experiments, the average temperatures and the average humidities in Wuhan were 19/11°C, 19/10°C, 21/15°C (day/night), and 79, 78, 83%, respectively. Besides, the total sunshine times in Wuhan were 158, 187, and 207 h, respectively, among the three experiments and the average light intensities were about 600–800 μmol•m−2•s−1.
At harvest, four uniform plants of each line were removed from the basin. After they had been sampled, the total roots were separated from the shoot base and primary root length (PRL) was investigated manually using a ruler, while shoots were over dried at 80°C until a constant weight to evaluate shoot dry weight (SDW). The intact root system were immersed and dispersed in a transparent plastic tray with water for scanning with a scanner (EPSON V700, Japan), and total root length (TRL), total root surface area (TSA), total root volume (TRV), total root number (TRN) were analyzed using WinRHIZO software (Pro, 2012b, Canada). Finally, roots were over dried at 80°C to evaluate root dry weight (RDW). Besides, total dry weight (TDW) and root-shoot ratio in dry weight (RSRD) were calculated.
The dried root and shoot samples were separately ground into powder, about 0.1 g were weighted and dissolved in H2SO4-H2O2, and then diluted with pure water to 1.25 L. Subsequently, N concentration (mg/L) were analyzed using Smartchem 200 automatic analyzer (Westco Scientific Instruments, Westco). Each sample was measured with three repetitions. RNC or SNC was N content (mg) in per unit weight (g) of root or shoot sample respectively and was calculated as follows: RNC or SNC = (sample N concentration × 1.25) / sample weight. The N uptake of root (RNU) and shoot (SNU) were calculated by multiplying weight by N concentration, respectively. Total N uptake (TNU) or N utilization efficiency (NUtE) was obtained by adding RNU and SNU together or by dividing TDW by TNU. All the traits investigated in this study were summarized in Table 1.
A total of 15 phenotypic traits, including 7 RM- and 8 NUE-related traits (Table 1) represented by the means of four plants for each genotype, were used for phenotypic analysis and QTL analysis. The broad-sense heritability (h2) of each measured trait was calculated as h2 = /( + /n + /nr), where , , and are the variance of genotype, genotype × environment, and error, respectively, and n and r are the number of independent experiments and replications, respectively. The estimation of , , and were obtained with the software SAS 9.2 (SAS Institute Inc., NC, USA) using the GLM procedure. Correlation analysis and principal component analysis (PCA) were calculated with software SAS 9.2 (SAS Institute Inc., NC, USA) using the PROC CORR and PROC PRINCOMP procedure, respectively.
The linkage mapping of QTLs was performed by composite interval mapping program (Zeng, 1994) using the Windows QTL Cartographer version 2.5 software (http://statgen.ncsu.edu/qtlcart/WQTLCart.htm). The corresponding genetic map consisted of 2264 unique loci/bins, which covered a total length of 2,107 cM distributed on 19 linkage groups (Yang et al., 2016). In this study, a walk speed of 1 cM, 5 control markers, a window size of 10 cM and forward regression method were used. The LOD threshold was determined by permutation analysis with 1,000 repetitions (Churchill and Doerge, 1994). The LOD threshold was set at 3.2–6.2 (p = 0.05) to identify significant QTLs. To avoid missing QTLs with a relatively small effect, a lower LOD threshold was set at 1.7–2.2 (p = 0.50) to detect suggestive QTLs. Both significant QTLs and overlapping suggestive QTLs were admitted (Long et al., 2007) and named as “identified QTL” (Shi et al., 2009).
In this study, “stable QTL (sQTL)” was defined, which was repeatedly detected with more than half of 2-LOD confidence interval overlapping at least two environments or different N levels and has the same additive-effect direction. In addition, sQTLs detected under both HN and LN conditions were called as “constitutive sQTL,” while sQTLs that were detected in at least two environments under either HN or LN conditions, respectively were named as “HN-specific sQTL” or “LN-specific sQTLs,” respectively (Li et al., 2015). These sQTLs were divided into two types: QTLs detected at least once with phenotypic variation explained (R2) ≥ 20% or at least twice with R2 ≥ 10% were named as “major sQTL”, and the remainder were named as “minor sQTL” (Price, 2006; Maccaferri et al., 2008; Shi et al., 2009). And QTL meta-analysis was performed using BioMercator 4.2 software (Arcade et al., 2004) to estimate the coincidences of several QTLs for different traits (at least two different traits), which were integrated into the “QTL cluster.” For a QTL cluster, the coincidence of QTLs for two or more traits with same additive-effect directions was considered to be positive, while the coincidence of QTLs with opposite additive-effect directions was considered negative (Coque et al., 2008).
Each identified QTL or QTL cluster was named as “q” + “the name of the trait abbreviation” or “qc” + “the linkage group,” respectively. Arabic numerals were added, if more than one QTL cluster was located on the same linkage group.
Phenotypic Variation for RM- and NUE-Related Traits in the BnaZNRIL Mapping Population of “Zhongshuang11” × “No. 73290” under Two Contrasting N Levels
Three independent hydroponic experiments for the parental lines and the RIL population were performed to evaluate 7 RM- and 8 NUE-related traits (Table 1) under both HN and LN growth conditions. Although the two parents (Zhongshuang11 and No.73290) showed no significant difference for most of RM- and NUE-related traits under HN level (except in E2, Zhongshuang11 showed significant higher values than No. 73290 in RDW, TRL, TSA, TRV, RNU, SNU, and TNU), they displayed significant differences for PRL, TRL, TSA, TRV, TRN, RDW, and NUtE in all three environments or a single environment under LN level (Table S1).
Continuous phenotypic distribution values with obvious kurtosis among these lines suggesting a quantitative inheritance pattern suitable for QTL identification. Minimum, maximum, mean values and coefficient of variations (CVs) for all investigated traits from each experiment were listed in Table 2. Most traits had considerable phenotypic variation within the BnaZNRIL population as suggested by the CVs ranging from 7.6 to 33.8% (Table 2). Increased PRL (1.6–89.9%), RSRD (54.5–172.7%), TRL (5.2–14.8%), and NUtE (11.5–200%) were response to LN stress, compared with those under HN levels (Table 2), indicating that LN stress stimulated root growth for more N nutrition uptake and improved N utilization, though shoot growth was obviously inhibited. The heritability (h2) of RM-related traits was relatively high, ranging from 0.38 to 0.67 under HN and from 0.48 to 0.60 under LN condition (Table 2). However, the heritability (h2) of NUE-related traits was moderate except SNC and NUtE (low under HN condition, 0.07 and 0.06, respectively), ranging from 0.23 to 0.48 under HN and from 0.26 to 0.59 under LN condition (Table 2). This indicated that RM-related traits are less affected than NUE-related traits on changeable environment under both HN and LN conditions. Although significant genotype × environment was observed, the relatively high heritability (h2) in the majority of the investigated traits indicated the genetic stability of these traits in a genotype among the three repetitions.
Table 2. Statistics for all investigated traits of BnZNRIL lines across three independent experiments (E1-E3).
Further, for each trait, part of RIL lines had values outside the range between parental lines, suggesting the presence of alleles with positive effects for root development and N utilization in both parents. For individual RM- or NUE-related traits, several of the transgressive lines always displayed extreme values lower or higher than parental lines over the three experiments (for example, B031 had 7–32% lower TRL and RDW than No. 73290, while B152 had 11–41% higher TRL and RDW than ZS11 under LN conditions), providing useful lines for RM or NUE genetic improvement and corresponding molecular mechanism dissection.
Genetic Relevance between RM- and NUE-Related Traits in the BnaZNRIL Mapping Population
The phenotypic correlations among all examined traits were calculated and the principal component analyses (PCA) for all investigated traits were performed to reveal the genetic correlation among RM- and NUE-related traits. For RM-related traits, the majority of traits had strong and significantly correlations (r = 0.71–0.93 under HN levels and r = 0.61–0.95 under LN conditions, P < 0.0001) with each other irrespective of N levels with the exception of PRL and RSRD (Table 3), indicating developmental relevance among these RM-related traits. PRL and RSRD showed smaller or no significant correlation with other RM-related traits under HN condition (r = 0.01–0.44 for PRL, r = 0.02–0.13 for RSRD with others). However, PRL and RSRD had significantly correlations with other RM-related traits under LN condition (r = 0.17–0.54 for PRL, P < 0.05; r = 0.25–0.38 for RSRD, P < 0.05; except r between RSRD and TRN) (Table 3), revealing the change of root system response to N stress, including PRL increasing and faster root development. For NUE-related traits, the dry biomass traits (DB, including SDW and TDW) had high values of correlation coefficients with N uptake traits (NU), SNU (r = 0.82–0.95), RNU (r = 0.61–0.71), and TNU (r = 0.77–0.95) irrespective of N levels, as expected (Table 3), suggesting that plant growth depend heavily on N uptake. On other hand, DB had no significant correlation with NUtE under HN level, while positive correlation with NUtE (r = 0.28) under LN level, indicating that plants rely on the higher N uptake and NUtE to maintain larger biomass under LN level (Table 3).
Table 3. Phenotypic correlation coefficients (r) among all investigated traits in the BnaZNRIL population grown under HN and LN levels.
Pearson's correlations were also calculated comparing the RM- and NUE-related traits. It was obviously that RM-related traits (except for PRL and RSRD) was significant correlated with NU and DB irrespective of N levels (r = 0.51–0.86 under HN condition, r = 0.43–0.82 under LN condition) (Table 3), indicating the phenotypic relationship between root system, plant growth and N use efficiency. Interestingly, although SNC and RNC showed no significant correlation with RM-related traits under HN condition, they are significantly correlated with several RM-related traits (TRL, TSA, and TRV) under LN condition.
Principle component analysis (PCA) also showed that RM-related traits except PRL and RSRD were all closely related and combined into two groups distributed in two different N levels, respectively (Figure 1). In NUE-related traits, only DB and NU were both assigned into the two RM-related groups. Thus, all these results showed that RM-related traits (except for PRL and RSRD) were more likely to be genetically associated with DB and NU, but not with NUtE, indicating that RM-related traits may play an extremely important role in plant N uptake rather than in plant N utilization.
Figure 1. PCA of the RIL population for RM- and NUE-related traits evaluated in hydroponics under the conditions of HN (triangle and circle respectively) and LN (cross and horizontal line respectively) levels. Thirty traits were projected onto the first and second principal components. Two obvious groups indicated by red circles: group1 or group2 included all of the RM- and NUE-related traits (except PRL, RSRD, RNC, SNC, and NUtE) under LN or HN conditions, respectively. The traits had the greatest positive influence over the determination of PC1 (>0.2) were all located in the two groups. Group1 had positive influence over the determination of PC2 while group2 had negative influence. Traits in each group were all closely related and combined into two clusters distributed in two different N levels, respectively.
Detection of QTLs Associated With RM- and NUE-Related Traits
Using the WinQTL cartographer software, a total of 129 identified QTLs (104 significant QTLs and 25 overlapping suggestive QTLs) associated with the investigated traits was detected on three environments, respectively. Among these, 91 QTLs were assumed to influence RM-related traits and 38 QTLs were assumed to influence NUE-related traits (Figure 2, Tables S2, S3). These QTLs were located across 12 of the 19 chromosomes in B. napus (including A01-A07, A09, C04, C05, C08, and C09). Of these, 45 (~35%) were located on the A07 linkage. More QTLs for NUE-related traits was detected under LN condition than that under HN level (7 under HN condition, 31 under LN condition), while similar number of QTLs for RM-related traits was detected under two N levels (43 under HN condition, 48 under LN condition; Table S3). Consequently, 56 out of 129 identified QTLs (~43% of the total number QTLs) were integrated into 23 sQTLs, which were repeatedly detected at least two environments or different N levels (Figure 2 and Table 4). Among these, 20 sQTLs were associated with RM-related traits and 3 sQTLs were associated with NUE-related traits. Notably, 16 out of 23 sQTLs were considered as “constitutive sQTL,” because they were detected under both HN and LN conditions. These results suggested that QTLs for RM-related traits are less affected than NUE-related traits on changeable environment under both HN and LN conditions.
Figure 2. Genomic regions of 129 QTLs detected for RM-related traits and NUE-related traits across three independent hydroponics experiments (E1-E3) under HN and LN levels. For each QTL symbol, red sign (+) or blue sign (–), represents positive additive effect from ZS11 or No. 73290 alleles respectively and the horizontal ordinate indicate the peak position. The rectangles represent the sQTL: red for LN-specific sQTL, blue for HN-specific sQTL and green for constitutive sQTL.
A total of 70 significant QTLs and 21 overlapping suggestive QTLs, explained 4.3–17.1% of the phenotypic variance, were detected with LOD scores ranged from 2.2 to 8.0 for RM-related traits. 23 QTLs were detected for PRL, 57 for root size (RS: including TRL, TSA, TRV, TRN, and RDW) and 9 for RSRD (Tables S2, S3). These QTLs were integrated into 20 sQTLs, including 14 constitutive sQTLs, 3 HN-specific sQTLs and 3 LN-specific sQTLs. Six sQTLs were found for PRL that explained 4.9–13.4% phenotypic variation, including 2 HN-specific sQTLs detected on chromosome A09 (Table 4). 13 sQTLs were found from RS (3 for TRL, 2 for TSA, 2 for TRV, 3 for TRN and 3 for RDW), which included 4 major sQTLs (all sQTLs for TRN and one for TSA), 3 LN-specific sQTLs (for TSA, TRN and RDW, respectively) and one HN-specific sQTL for TRL. Besides, only one sQTL located on chromosome A02 was found for RSRD. Notably, most of these sQTLs (~69%) were densely distributed on 57.4–111.2 cM of chromosome A07, suggesting the presence of major genes regulating RM in this region. In particular, the LN-specific sQTLs for RDW (sqRDW.A09) was co-localized with the LN-specific sQTL for TSA (sqTSA.A09) on chromosome A09 (68.9–104.9 cM), though the two sQTLs had minor genetic effects explaining 4.4–7.9% of the phenotypic variations.
Relatively fewer identified QTLs were detected for NUE-related traits, including 34 significant QTLs and 4 overlapping suggestive QTLs, explained 5.5–12.6% of the phenotypic variance, with LOD scores ranged from 2.3 to 5.8. These QTLs included 6 QTLs for SDW, 2 for TDW, 10 for N concentration (NC, including RNC and SNC), 17 for N uptake (NU, including RNU, SNU, and TNU) and 3 for NUtE (Tables S2, S3). Of these QTLs, only three sQTLs were found, including one LN-specific sQTLs for SDW and two constitutive sQTLs (one each for SDW and SNU; Table 4). This LN-specific sQTL for SDW (sqSDW.A01) was located on 86.6–93 cM of chromosome A01 and explained 5.5–7.4% of phenotypic variation. The two constitutive sQTLs were co-localized on chromosome A04 (77.3–83 cM) and explained 6.2–9.1% of phenotypic variation.
Detection of Pleiotropic or Specific QTL Clusters for RM- and NUE-Related Traits Revealed Two NUE-Improving Approaches in Rapeseed
Most pairs of RM-and NUE-related traits showed significantly genetic correlation in Table 3, providing the potential to identify QTLs with pleiotropy or specificity for these traits. QTL meta-analysis was performed to find pleiotropic “QTL cluster.” As a result, a total of 83 identified QTLs (~64% of the total number) were integrated into 22 distinct pleiotropic QTL clusters (Figure 3, Table 5). The 22 distinct QTL clusters were distributed on eight chromosomes (A01, A02, A04-A07, A09, and C04) with a range of 2–13 identified QTLs. About 68% of these were distributed on chromosome A01 (four), A02 (three), and A07 (eight). Besides, QTL clusters detected for both RM-related traits and NUE-related traits were defined as “RM-NUE QTL cluster,” while QTL clusters that were exclusively detected for RM-related traits or NUE-related traits were named as “RM-specific QTL cluster” or “NUE-specific QTL cluster,” respectively. Thus, all of the 22 distinct QTL clusters were divided into 10 (~45%) RM-specific QTL clusters, 6 (~27%) NUE-specific QTL clusters, and 6 (~27%) RM-NUE QTL clusters.
Figure 3. Identified QTL clusters for all investigated traits as revealed by meta-QTL analysis. The genetic distance where the QTL cluster located is represented by a vertical ruler and the chromosome number is marked on top. Each QTL cluster is indicated by the opening brace across all chromosomes and its name is give on the left of the chromosome. Red opening brace is marked if QTLs for QTLs for RM-related traits and NUE-related traits are located in a same cluster else black. The QTLs contained in a QTL cluster are given on the right of the chromosome and QTLs detected at HN levels are indicated in black, and at LN levels in red. The number between parentheses indicates how many experiments this QTL was detected in. The sign, (+) or (–), represents positive additive effects from ZS11 and No. 73290 alleles, respectively.
Table 5. Information of the detected QTL clusters for all investigated traits in different environments.
Pleiotropic QTL Clusters for RM-Related Traits
Among the ten RM-specific QTL clusters, half of them were located on chromosome A07 and the remainders were on A01 (two), A02 (one), A05 (one) and A09 (one) (Figure 3, Table 5). On chromosome A07, the five QTL clusters were densely distributed on 62.46–74.22 cM, comprising seven sQTLs (including three major sQTLs, one LN-specific sQTLs, one HN-specific sQTL as well as five constitutive sQTLs). Notably, the largest RM-specific QTL clusters, qcA07-5 included up to 13 identified RM-related QTLs and was co-localized with five stable QTLs (two major sQTLs, one LN-specific sQTL as well as four constitutive sQTLs), indicating that qcA07-5 have a higher possibility of contributing specifically to root development under both HN and LN levels. Besides, qcA09-2 was detected only under LN level and co-localized with two LN-specific sQTL (sqTSA.A09 and sqRDW.A09), suggesting that qcA09-2 may be closely related to root growth under LN stress.
Pleiotropic QTL Clusters for NUE-Related Traits
For the six NUE-specific QTL clusters, five of them was detected only under LN level, indicating the presence of LN stress stimulated genetic mechanism for NUE (Figure 3, Table 5). Among them, the largest QTL cluster, qcA04, was co-localized with two stable QTLs (constitutive sQTLs for SDW and SNU, respectively), showing that qcA04 may specifically contribute to shoot development and N uptake. Besides, qcA01-3 was co-located with the LN-specific sQTL, sqSDWA01.
Pleiotropic QTL Clusters for Both RM- and NUE-Related Traits
Of the six RM-NUE QTL cluster, five clusters was detected under both HN and LN condition (Figure 3, Table 5). On the other hand, four RM-NUE QTL clusters were overlapped between RM and NU traits; two were between RM and DB and others were between RM and NC, but no were between RM and NUtE. These results indicated that the co-localization of RM-related traits and NUE-related traits was existed under both HN and LN conditions, while RM-related traits were closely related with plant N uptake instead of N utilization efficiency. Five RM-NUE QTL clusters (qcA02-1, qcA07-6, qcA07-7, qcA07-8, and qcA09-1) had same directions of additive-effects. Among them, three QTL clusters had the favorable alleles coming from No. 73290 (qcA07-6, qcA07-7, and qcA07-8), while two QTL clusters had the favorable alleles derived from ZS11 (qcA02-3 and qcA09-1). The three largest RM-NUE QTL clusters (qcA07-6, qcA07-7, and qcA07-8) were densely distributed on 90.12–107.84 cM of chromosome A07 and co-localized with two constitutive stable QTLs, sqTRV.A07-2 and sqRDW.A07-2, implicating the presence of major genes regulating root growth and plant N uptake in this region. In addition, qcA02-1 and qcA09-1 were co-localized with two constitutive stable QTLs, sqRSRD.A02 and sqPRL.A09-3, respectively.
Taken together, three types of QTL clusters (RM-NUE QTL clusters and RM or NUE-specific QTL clusters) implied two NUE improvement approaches were existed in rapeseed: RM-based and N utilization-based direct approaches. The credible genetic relationship between RM and NUE in rapeseed laid a solid theoretical foundation for RM-based approach to NUE genetic improvement. More importantly, more QTLs for RM-related traits (91) were identified than that for NUE-related traits (38) in this study, and all of the four major sQTLs and most of stable QTLs (20 out of 23) were detected to be related to RM traits under HN and/or LN levels, indicated that the manipulation of RM for improving NUE would be more feasible and reliable than regulating nitrogen efficiency directly in rapeseed.
Genetic Relationship between RM and NUE Traits Implied Two NUE Improvement Approaches in Rapeseed: RM-Based and N Utilization-Based Direct Approaches
Root system is the key place for N acquisition (Lynch, 1995). Illuminating the genetic relationship between RM and NUE traits could provide the basis of RM-based approach for NUE genetic improvement. Our study uncovered the tight connection between RM and NUE traits from two aspects. On the one hand, phenotypic correlation analysis and principal component analysis showed high genetic correlation between RM and NUE traits (Figure 1, Table 3). Generally, RM-related traits (except for PRL and RSRD) had high positive correlations with DB and NU irrespective of N levels, but no correlation or little correlation with NUtE under HN or LN level respectively. PCA also showed that RM-related traits (except for PRL and RSRD) and NUE-related traits (except for NC and NUtE) were assigned into the same group under HN and LN conditions, respectively (Figure 1, group 1 and group 2). These results indicated that RM, especially “root size,” play a critical role in plant N uptake, rather than in N utilization, which was consistent with the previous report on maize (Li et al., 2015). On the other hand, the co-localization of QTLs for RM- and NUE-related traits in this study further supports the tight connection between RM and NUE traits. Among the six RM-NUE QTL clusters, four overlapped between RM and NU, two between RM and DB and others between RM and NC, but no between RM and NUtE, and five of them were positive in which all the QTLs had same additive-effect directions (Figure 3, Table 5). These proved that N acquisition, rather than N utilization, is more likely correlated with root morphology.
In conclusion, the significantly phenotypic correlation and the co-localization of QTLs fully demonstrated the presence of credible genetic relationship between RM and NUE in rapeseed. However, the process of plant N utilization is complex, mainly including N uptake by roots and translocation, assimilation and remobilization inside the plant (Xu et al., 2012). RM and root uptake activity (including root vigor, root N transport, etc.) are the major determinants of N acquisition (Glass, 2003; Garnett et al., 2009; Xu et al., 2012), while the N translocation, assimilation and remobilization refer to N metabolism process (Xu et al., 2012). So, given the presence of six NUE-specific QTL clusters in this study (Figure 3, Table 5) and the credible genetic relationship between RM and NUE, our results provided evidences for two NUE improvement approaches in rapeseed: RM-based and N utilization-based direct approaches.
Important Pleiotropic QTL Clusters Used for NUE-Improving in Rapeseed
Both RM and NUE are complex traits, susceptible to growth environment, and difficult in directly measuring that restrict genetic improvement of NUE and RM in crops (Hochholdinger and Tuberosa, 2009). Evaluating these traits using the genomic approach and marker-assisted selection are worldwide challenges for NUE improvement (Coque et al., 2008; Liu et al., 2008; Li et al., 2015). In this study, 129 identified QTLs (including 23 stable QTLs and 4 major QTLs) were detected and 83 QTLs of them were integrated into 22 pleiotropic QTL clusters (Figures 2, 3, Tables 4, 5, Tables S2, S3). Five RM-NUE (qcA02-1, qcA07-6 to qcA07-8 and qcA09-1), ten RM-specific (qcA01-1, qcA01-2, qcA02-3, qcA07-1 to qcA07-5, and qcA09-2) and three NUE-specific (qcA01-3, qcA04, and qcA06) QTL clusters with same directions of additive-effect provided the potential to identify pleiotropic genes or different closely-linked genes controlling root development and N uptake. More importantly, the favorable alleles among these QTL clusters can be used directly for genetic improvement of NUE via a marker-assisted selection approach in rapeseed.
The three largest RM-NUE QTL clusters (qcA07-6, qcA07-7, and qcA07-8) were closely distributed on 90.12–107.84 cM of chromosome A07 and all had the favorable alleles coming from No. 73290. These QTL clusters co-localized with two constitutive stable QTLs, sqTRV.A07-2, and sqRDW.A07-2, suggesting that this region may play a key role in root growth and plant N uptake. Besides, three single-nucleotide polymorphisms (SNPs), significantly associated with RSRD, RDW and/or PRL, which were detected under low phosphorus (P) condition (Wang et al., 2017) were located in the similar genomic region with the three QTL clusters, respectively. In addition, Wang et al. (2017) identified another significantly associated SNP for root length, which was most likely co-localized with qcA02-1. Another QTL cluster, qcA09-1, located on 28.96–38.45 cM, was associated with PRL and RNU/TNU, whereas a QTL for PRL have previously been identified nearby (Zhang et al., 2016). Further, a significantly associated SNP for shoot concentration of Na discovered by Bus et al. (2014) was also localized in this genomic region. Consequently, we proposed that these positive RM-NUE QTL clusters were likely to be the common regions for regulating rapeseed RM traits across different genetic backgrounds, and further involved in plant N, P, and Na uptake.
Five RM-specific QTL clusters on chromosome A07 (from qcA07-1 to qcA07-5) were distributed on 62.46–74.22 cM. Seven sQTLs (including three major sQTLs) for RM-related traits were co-localized with these QTL clusters, suggesting the presence of major genes regulating RM in this region. Besides, under low P condition, Wang et al. (2017) identified six significantly associated SNPs with several root traits and Zhang et al. (2016) detected two QTLs for PRL and SDW that were located in the similar genomic region. Interestingly, all the QTLs distributed on 62.46–74.22 cM of chromosome A07 had the favorable alleles coming from the parental line ZS11, while on 90.12–107.84 cM from the parental line No. 73290, indicating that there are two distinct but important regulatory mechanisms for root development on chromosome A07 under this mapping population. On chromosome A01, qcA01-1 and qcA01-2 were distributed on 43.05–50.41 cM and the major stable QTL sqTRN.A01 was co-localized with these QTL clusters, suggesting the presence of major genes regulating RM in this region. The QTL uq.A1_1 for TRL detected by Zhang et al. (2016) using another genetic population was located in the similar genomic region. Thus, given the credible genetic relationship between RM and NUE, these important QTL clusters would have a considerable breeding value for RM-based NUE improvement in rapeseed, though no co-localization with NUE detected in this study.
RM-Based NUE Improvement Would be More Effective than N Utilization-Based Direct NUE Improvement Via Marker-Assisted Selection in Rapeseed
Our study suggested two different approaches for NUE improvement via marker-assisted selection in rapeseed: RM-based and N utilization-based direct approaches. Previous reports showed that RM-related traits investigated at an early development stage had significant links with NUE, yield and yield components (Li et al., 2015; Xie et al., 2017). These suggested that breeding varieties with large “root size” at seedling stage may be a promising approach to optimize N uptake efficiency in crops. RM-based NUE improvement provided a novel and exciting opportunity for the manipulation of RM via marker-assisted selection QTLs for improving NUE in rapeseed, though direct NUE improvement had made some achievements as the major approach for NUE breeding previously (Xu et al., 2012; Han et al., 2015). The credible genetic relationship between RM and NUE in this study laid a theoretical foundation for RM-based approach on NUE genetic improvement in rapeseed. Furthermore, in this study, ~2.4 times more QTLs were identified for RM-related traits (91) than that for NUE-related traits (38) and all of the four major QTLs and most of stable QTLs (20 out of 23) were detected to be related to RM traits under HN and/or LN levels, suggested that QTLs for RM-related traits are less affected than NUE-related traits on changeable environment under both HN and LN conditions. Thus, regulating RM to improve NUE via marker-assisted selection would be more feasible and reliable than regulating nitrogen efficiency directly in rapeseed. On the other hand, many co-localized QTLs associated with RM and NUE traits have been reported in crops (An et al., 2006; Coque et al., 2008; Liu et al., 2008; Li et al., 2015) and manipulation of RM via marker-assisted selection to improve nutrient efficiency has been successfully used in rice (Wissuwa, 2005; Pariasca-Tanaka et al., 2014; all for phosphorus use efficiency, PUE) and maize (Li et al., 2015, for NUE; Gu et al., 2016, for PUE). Moreover, Li et al. (2015) identified 53 advanced backcross-derived lines (ABLs) containing RM-NUE QTL clusters via marker-assisted selection of RM and the GY/NUE of these ABLs showed apparent mean increases of 13.8% under HN and of 15.9% under LN conditions, compared to recurrent background Wu312. All this indicated that manipulation of RM via marker-assisted selection to improve NUE would be more effective than marker-assisted selection for NUE directly in rapeseed.
Though nutrient absorption and root development of plants are different in hydroponic culture and field experiment (Chapman et al., 2012; Watt et al., 2013; Petrarulo et al., 2015), some of the experimental evidence available has shown their correlations (Lynch and Brown, 2001; Cui et al., 2008). Moreover, QTLs for root traits detected in hydroponics have been successfully used in nutrient efficiency improvement in field (Li et al., 2015, for NUE; Gu et al., 2016, for PUE). We deduced that the QTLs identified in this study could be used for marker-assisted selection in field, and corresponding examination will be conducted in the field in progress.
Development of N efficient cultivars is beneficial for environment-friendly and sustainable rapeseed production. As shown in this study, RM had a significant positive phenotypic correlation with NUE and many RM-related QTLs were co-localized with QTLs for NUE, providing the credible genetic evidence for the significant associations between RM and NUE. In addition, 129 identified QTLs (including 23 stable QTLs and 4 major sQTLs) were detected and 83 QTLs of them were integrated into 22 pleiotropic QTL clusters. Five RM-NUE, ten RM-specific and three NUE-specific QTL clusters with same directions of additive-effect implied two distinct NUE-improving approaches (RM-based and N utilization-based directly) via marker-assisted selection in rapeseed. Importantly, all of four major QTLs and most of stable QTLs (20 out of 23) were detected to be related to RM traits under HN and/or LN levels, suggested that regulating RM to improve NUE would be more effective than regulating N efficiency directly. These pleiotropic QTL clusters also provided valuable genomic regions for RM-based or direct NUE-improvement. Whether the co-localization of QTLs was originated from pleiotropic genes or different closely-linked genes, fine mapping, RNA-seq and bioinformatics analysis should be conducted to confirm corresponding gene or genes in the future, which will lay a solid foundation for uncovering the genetic and molecular mechanism on NUE improvement in rapeseed.
JW performed the data analysis and wrote the main manuscript text. XD designed and managed the experiments. JS contributed to the construction of molecular genetic map. XW and GL prepared the plant materials. HW reviewed the manuscript. All of the authors have read and approved the final manuscript.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
This work was supported by the National Key Basic Research Program of China (2015CB150200), the National Natural Science Foundation of China (31501820), the National Key Research and Development Program of China (2016YFD0101300), Rapeseed Industry Technology System (CARS-13), and the Agricultural Science and Technology Innovation Project (CAAS-ASTIP-2013-OCRI).
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fpls.2017.01709/full#supplementary-material
DB, dry biomass traits; HN, high-nitrogen; LN, low-nitrogen; NC, N concentration; NU, N uptake traits; NUE, nitrogen use efficiency; NUtE, N utilization efficiency; PRL, primary root length; SNP, single-nucleotide polymorphism; sQTL, stable QTL; RDW, root dry weight; RIL, recombinant inbred line; RM, root morphology; RNC, root N concentration; RNU, root N uptake; RSRD, root-shoot ratio in dry weight; SDW, shoot dry weight; SNC, shoot N concentration; SNU, shoot N uptake; TDW, total dry weight; TNU, total N uptake; TRL, total root length; TRN, total root number; TRV, total root volume; TSA, total root surface area.
An, D., Su, J., Liu, Q., Zhu, Y., Tong, Y., Li, J., et al. (2006). Mapping QTLs for nitrogen uptake in relation to the early growth of wheat (Triticum aestivum L.). Plant Soil 284, 73–84. doi: 10.1007/s11104-006-0030-3
Arcade, A., Labourdette, A., Falque, M., Mangin, B., Chardon, F., Charcosset, A., et al. (2004). BioMercator: integrating genetic maps and QTL towards discovery of candidate genes. Bioinformatics 20, 2324–2326. doi: 10.1093/bioinformatics/bth230
Bouchet, A. S., Laperche, A., Bissuel-Belaygue, C., Baron, C., Morice, J., Rousseau-Gueutin, M., et al. (2016). Genetic basis of nitrogen use efficiency and yield stability across environments in winter rapeseed. BMC Genet. 17:131. doi: 10.1186/s12863-016-0432-z
Bouchet, A. S., Nesi, N., Bissuel, C., Bregeon, M., Lariepe, A., Navier, H., et al. (2014). Genetic control of yield and yield components in winter oilseed rape (Brassica napus L.) grown under nitrogen limitation. Euphytica 199, 183–205. doi: 10.1007/s10681-014-1130-4
Bus, A., Korber, N., Parkin, I. A., Samans, B., Snowdon, R. J., Li, J., et al. (2014). Species- and genome-wide dissection of the shoot ionome in Brassica napus and its relationship to seedling development. Front. Plant Sci. 5:485. doi: 10.3389/fpls.2014.00485
Chardon, F., Barthelemy, J., Daniel-Vedele, F., and Masclaux-Daubresse, C. (2010). Natural variation of nitrate uptake and nitrogen use efficiency in Arabidopsis thaliana cultivated with limiting and ample nitrogen supply. J. Exp. Bot. 61, 2293–2302. doi: 10.1093/jxb/erq059
Coque, M., Martin, A., Veyrieras, J. B., Hirel, B., and Gallais, A. (2008). Genetic variation for N-remobilization and postsilking N-uptake in a set of maize recombinant inbred lines. 3. QTL detection and coincidences. Theor. Appl. Genet. 117, 729–747. doi: 10.1007/s00122-008-0815-2
Cui, K., Huang, J., Xing, Y., Yu, S., Xu, C., and Peng, S. (2008). Mapping QTLs for seedling characteristics under different water supply conditions in rice (Oryza sativa). Physiol. Plant 132, 53–68. doi: 10.1111/j.1399-3054.2007.00991.x
Dun, X., Tao, Z., Wang, J., Wang, X., Liu, G., and Wang, H. (2016). Comparative transcriptome analysis of primary roots of Brassica napus seedlings with extremely different primary root lengths using RNA sequencing. Front. Plant Sci. 7:1238. doi: 10.3389/fpls.2016.01238
Galloway, J. N., Townsend, A. R., Erisman, J. W., Bekunda, M., Cai, Z., Freney, J. R., et al. (2008). Transformation of the nitrogen cycle: recent trends, questions, and potential solutions. Science 320, 889–892. doi: 10.1126/science.1136674
Good, A. G., Shrawat, A. K., and Muench, D. G. (2004). Can less yield more? is reducing nutrient input into the environment compatible with maintaining crop production? Trends Plant Sci. 9, 597–605. doi: 10.1016/j.tplants.2004.10.008
Gu, R., Chen, F., Long, L., Cai, H., Liu, Z., Yang, J., et al. (2016). Enhancing phosphorus uptake efficiency through QTL-based selection for root system architecture in maize. J. Genet. Genomics 43, 663–672. doi: 10.1016/j.jgg.2016.11.002
Han, M., Okamoto, M., Beatty, P. H., Rothstein, S. J., and Good, A. G. (2015). The genetics of nitrogen use efficiency in crop plants. Annu. Rev. Genet. 49, 269–289. doi: 10.1146/annurev-genet-112414-055037
Hoffmann, A., Maurer, A., and Pillen, K. (2012). Detection of nitrogen deficiency QTL in juvenile wild barley introgression lines growing in a hydroponic system. BMC Genet. 13:88. doi: 10.1186/1471-2156-13-88
Huang, G., Yi, K. K., Wu, Y. R., Zhu, L., Mao, C. Z., and Wu, P. (2004). QTLs for nitrate induced elongation and initiation of lateral roots in rice (Oryza sativa L.). Plant Soil 263, 229–237. doi: 10.1023/B:PLSO.0000047736.66712.c6
Li, P., Chen, F., Cai, H., Liu, J., Pan, Q., Liu, Z., et al. (2015). A genetic relationship between nitrogen use efficiency and seedling root traits in maize as revealed by QTL analysis. J. Exp. Bot. 66, 3175–3188. doi: 10.1093/jxb/erv127
Li, P., Zhuang, Z., Cai, H., Cheng, S., Soomro, A. A., Liu, Z., et al. (2016). Use of genotype-environment interactions to elucidate the pattern of maize root plasticity to nitrogen deficiency. J. Integr. Plant Biol. 58, 242–253. doi: 10.1111/jipb.12384
Lian, X., Xing, Y., Yan, H., Xu, C., Li, X., and Zhang, Q. (2005). QTLs for low nitrogen tolerance at seedling stage identified using a recombinant inbred line population derived from an elite rice hybrid. Theor. Appl. Genet. 112, 85–96. doi: 10.1007/s00122-005-0108-y
Linkohr, B. I., Williamson, L. C., Fitter, A. H., and Leyser, H. M. O. (2002). Nitrate and phosphate availability and distribution have different effects on root system architecture of Arabidopsis. Plant J. 29, 751–760. doi: 10.1046/j.1365-313X.2002.01251.x
Liu, J., Li, J., Chen, F., Zhang, F., Ren, T., Zhuang, Z., et al. (2008). Mapping QTLs for root traits under different nitrate levels at the seedling stage in maize (Zea mays L.). Plant Soil 305, 253–265. doi: 10.1007/s11104-008-9562-z
Long, Y., Shi, J., Qiu, D., Li, R., Zhang, C., Wang, J., et al. (2007). Flowering time quantitative trait loci analysis of oilseed brassica in multiple environments and genomewide alignment with Arabidopsis. Genetics 177, 2433–2444. doi: 10.1534/genetics.107.080705
Maccaferri, M., Sanguineti, M. C., Corneti, S., Ortega, J. L., Salem, M. B., Bort, J., et al. (2008). Quantitative trait loci for grain yield and adaptation of durum wheat (Triticum durum Desf.) across a wide range of water availability. Genetics 178, 489–511. doi: 10.1534/genetics.107.077297
Mu, X., Chen, F., Wu, Q., Chen, Q., Wang, J., Yuan, L., et al. (2015). Genetic improvement of root growth increases maize yield via enhanced post-silking nitrogen uptake. Eur. J. Agron. 63, 55–61. doi: 10.1016/j.eja.2014.11.009
Pariasca-Tanaka, J., Chin, J. H., Drame, K. N., Dalid, C., Heuer, S., and Wissuwa, M. (2014). A novel allele of the P-starvation tolerance gene OsPSTOL1 from African rice (Oryza glaberrima Steud) and its distribution in the genus Oryza. Theor. Appl. Genet. 127, 1387–1398. doi: 10.1007/s00122-014-2306-y
Pestsova, E., Lichtblau, D., Wever, C., Presterl, T., Bolduan, T., Ouzunova, M., et al. (2016). QTL mapping of seedling root traits associated with nitrogen and water use efficiency in maize. Euphytica 209, 585–602. doi: 10.1007/s10681-015-1625-7
Petrarulo, M., Marone, D., Ferragonio, P., Cattivelli, L., Rubiales, D., De Vita, P., et al. (2015). Genetic analysis of root morphological traits in wheat. Mol. Genet. Genomics 290, 785–806. doi: 10.1007/s00438-014-0957-7
Rathke, G. W., Christen, O., and Diepenbrock, W. (2005). Effects of nitrogen source and rate on productivity and quality of winter oilseed rape (Brassica napus L.) grown in different crop rotations. Field Crops Res. 94, 103–113. doi: 10.1016/j.fcr.2004.11.010
Shang, L., Cai, S., Ma, L., Wang, Y., Abduweli, A., Wang, M., et al. (2016). Seedling root QTLs analysis on dynamic development and upon nitrogen deficiency stress in upland cotton. Euphytica 207, 645–663. doi: 10.1007/s10681-015-1564-3
Shi, J., Li, R., Qiu, D., Jiang, C., Long, Y., Morgan, C., et al. (2009). Unraveling the complex trait of crop yield with quantitative trait loci mapping in Brassica napus. Genetics 182, 851–861. doi: 10.1534/genetics.109.101642
Shi, L., Shi, T., Broadley, M. R., White, P. J., Long, Y., Meng, J., et al. (2013). High-throughput root phenotyping screens identify genetic loci associated with root architectural traits in Brassica napus under contrasting phosphate availabilities. Ann. Bot. 112, 381–389. doi: 10.1093/aob/mcs245
Stahl, A. (2015). Detailed Phenotypic Assessment of Genetic Variation for Nitrogen Uptake and Utilisation Efficiency Traits in a Diverse Panel of Brassica napus L. Ph.D. Thesis, Justus Liebig University.
Wang, X., Chen, Y., Thomas, C. L., Ding, G., Xu, P., Shi, D., et al. (2017). Genetic variants associated with the root system architecture of oilseed rape (Brassica napus L.) under contrasting phosphate supply. DNA Res. 24, 407–417. doi: 10.1093/dnares/dsx013
Watt, M., Moosavi, S., Cunningham, S. C., Kirkegaard, J. A., Rebetzke, G. J., and Richards, R. A. (2013). A rapid, controlled-environment seedling root screen for wheat correlates well with rooting depths at vegetative, but not reproductive, stages at two field sites. Ann. Bot. 112, 447–455. doi: 10.1093/aob/mct122
Wissuwa, M. (2005). Combining a modelling with a genetic approach in establishing associations between genetic and physiological effects in relation to phosphorus uptake. Plant Soil 269, 57–68. doi: 10.1007/s11104-004-2026-1
Xie, Q., Fernando, K. M., Mayes, S., and Sparkes, D. L. (2017). Identifying seedling root architectural traits associated with yield and yield components in wheat. Ann. Bot. 119, 1115–1129. doi: 10.1093/aob/mcx001
Yang, M., Ding, G., Shi, L., Feng, J., Xu, F., and Meng, J. (2010). Quantitative trait loci for root morphology in response to low phosphorus stress in Brassica napus. Theor. Appl. Genet. 121, 181–193. doi: 10.1007/s00122-010-1301-1
Yang, Y., Shi, J., Wang, X., Liu, G., and Wang, H. (2016). Genetic architecture and mechanism of seed number per pod in rapeseed: elucidated through linkage and near-isogenic line analysis. Sci. Rep. 6:24124. doi: 10.1038/srep24124
Yu, P., White, P. J., Hochholdinger, F., and Li, C. (2014). Phenotypic plasticity of the maize root system in response to heterogeneous nitrogen availability. Planta 240, 667–678. doi: 10.1007/s00425-014-2150-y
Zhang, D., Hua, Y., Wang, X., Zhao, H., Shi, L., and Xu, F. (2014). A high-density genetic map identifies a novel major QTL for boron efficiency in oilseed rape (Brassica napus L.). PLoS ONE 9:e112089. doi: 10.1371/journal.pone.0112089
Zhang, Y., Thomas, C. L., Xiang, J., Long, Y., Wang, X., Zou, J., et al. (2016). QTL meta-analysis of root traits in Brassica napus under contrasting phosphorus supply in two growth systems. Sci. Rep. 6:33113. doi: 10.1038/srep33113
Keywords: root morphology, nitrogen use efficiency, genetic relationship, QTL clusters, rapeseed
Citation: Wang J, Dun X, Shi J, Wang X, Liu G and Wang H (2017) Genetic Dissection of Root Morphological Traits Related to Nitrogen Use Efficiency in Brassica napus L. under Two Contrasting Nitrogen Conditions. Front. Plant Sci. 8:1709. doi: 10.3389/fpls.2017.01709
Received: 03 July 2017; Accepted: 19 September 2017;
Published: 29 September 2017.
Edited by:Maoteng Li, Huazhong University of Science and Technology, China
Reviewed by:Sheng Chen, University of Western Australia, Australia
Fangsen Xu, Huazhong Agricultural University, China
Copyright © 2017 Wang, Dun, Shi, Wang, Liu and Wang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Hanzhong Wang, email@example.com
†These authors have contributed equally to this work.