Vijayalakshmi Pujarula
Madhu Pusuluri
Srikanth Bollam
Roma Rani Das
Rambabu Ratnala
Gopikrishna Adapala
Vishnukiran Thuraga
Abhishek Rathore
Rakesh K. Srivastava*
Rajeev Gupta*- International Crops Research Institute for the Semi-Arid Tropics (ICRISAT), Patancheru, India
Nitrogen (N) is one of the primary macronutrients required for crop growth and yield. This nutrient is especially limiting in the dry and low fertility soils where pearl millet [Pennisetum glaucum (L.) R. Br] is typically grown. Globally, pearl millet is the sixth most important cereal grown by subsistence farmers in the arid and semi-arid regions of sub-Saharan Africa and the Indian subcontinent. Most of these agro-ecologies have low N in the root zone soil strata. Therefore, there is an immense need to identify lines that use nitrogen efficiently. A set of 380 diverse pearl millet lines consisting of a global diversity panel (345), parents of mapping populations (20), and standard checks (15) were evaluated in an alpha-lattice design with two replications, 25 blocks, a three-row plot for 11 nitrogen use efficiency (NUE) related traits across three growing seasons (Summer 2017, Rainy 2017, and Summer 2018) in an N-depleted precision field under three different N levels (0%-N0, 50%-N50, 100%-N100 of recommended N, i.e., 100 kg ha–1). Analysis of variance revealed significant genetic variation for NUE-related traits across treatments and seasons. Nitrogen in limited condition (N0) resulted in a 27.6 and 17.6% reduction in grain yield (GY) and dry stover yield (DSY) compared to N50. Higher reduction in GY and DSY traits by 24.6 and 23.6% were observed under N0 compared to N100. Among the assessed traits, GY exhibited significant positive correlations with nitrogen utilization efficiency (NUtE) and nitrogen harvest index (NHI). This indicated the pivotal role of N remobilization to the grain in enhancing yield levels. Top 25 N-insensitive (NIS-top grain yielders) and N-sensitive (NS-poor grain yielders) genotypes were identified under low N conditions. Out of 25 NIS lines, nine genotypes (IP 10820, IP 17720, ICMB 01222-P1, IP 10379, ICMB 89111-P2, IP 8069, ICMB 90111-P2, ICMV IS89305, and ICMV 221) were common with the top 25 lines for N100 level showing the genotype plasticity toward varying N levels. Low N tolerant genotypes identified from the current investigation may help in the identification of genomic regions responsible for NUE and its deployment in pearl millet breeding programs through marker-assisted selection (MAS).
Introduction
Global nitrogen (N) demand, one of the most expensive farm inputs currently stands at about 117 million metric tons with a projected annual increase of ∼1.5% in the future (FAO, 2019). Indian agriculture consumes over 17 million tons of N fertilizer per year. However, plants can utilize only 30–40% of the applied N for food production and the remaining (up to 60%) is lost to the environment by leaching, de-nitrification, and runoff. Surplus nitrogen pollutes freshwater streams and air which is hazardous to the majority of living species (Hickman et al., 2014; Russo et al., 2017).
Besides, excess usage of N fertilizer not only decreases the efficiency of nutrient use but also affects the rate of economic returns per unit of chemical fertilizer applied. The effect of negative environmental and economic impacts could be reduced through better agronomic practices and also by utilizing N efficient lines with improved nitrogen use efficiency (NUE) (Raun et al., 1999). Hence, improving the NUE of crop plants could help in reducing fertilizer input, increased productivity and profitability coupled with a reduced negative impact on the environment.
Globally, pearl millet (Pennisetum glaucum (L.) R. Br) widely grown for food and fodder is considered as the sixth most important cereal crop in terms of area. It is one of the oldest cultivated cereal crops, originated in Africa but later spread to many countries (D’Andrea and Casey, 2002; Manning et al., 2011). This millet crop is vital for food and nutritional security for the world’s poorest people living in different agro-ecological zones. Besides, pearl millet has a rich nutritional profile, wide genetic diversity, high photosynthetic rates being a C4 and also provides a healthy balanced diet, thus contributes to the economic security of poor farmers (Srivastava et al., 2019). In the Indian agriculture scenario, pearl millet is grown by poor and marginal farmers under low fertile and rainfall dependent areas often facing drought during different stages of crop growth. From the past few years, pearl millet demand has been increased for its nutritional characteristics and its adaptability to a wide range of climatic conditions (Tako et al., 2015). Since then farmers started applying a high dose of N fertilizer for maximizing grain yield (GY), but excessive usage of N fertilizer leads to low NUE of the crop along with various environmental hazards. In some parts of India, poor farmers lack knowledge about the ideal dose of fertilizer to harvest the real yield potential. Hence, yield improvement of pearl millet under low nitrogen input is indeed beneficial for economic and environmentally sustainable cultivation.
Nitrogen use efficiency is a complex trait, which is associated with various morphological, physiological, molecular, and biochemical changes in plants throughout the life cycle. For a clear understanding of this complex nature, studies of various physiological traits and their close correlation with one or more economically important traits like GY is foremost critical. Nevertheless, it will help in selecting low N tolerant/high yielding lines at different N conditions (Monostori et al., 2016). In general, plant function is always associated with chlorophyll content, which directly indicates the N status of the leaf (Yang et al., 2014). Leaf chlorophyll content and photosynthetic capacity are appropriate benchmarks for identifying high NUE (HNUE) genotypes under low N conditions during field trials (Vijayalakshmi et al., 2015; Kiran et al., 2016). Leaf N status was usually measured by using a hand-held optical chlorophyll meter to monitor the leaf nitrogen status and chlorophyll content. Moreover, leaf area (LA) is also one of the important physiological traits, plays a critical role in enhancing plant biomass (Gianquinto et al., 2004; Monostori et al., 2016).
Nitrogen use efficiency is defined as the ability to produce GY per unit N available in the soil. NUE mainly depends on the results of two main processes, such as N uptake efficiency (NupE), and nitrogen utilization efficiency (NutE) (Good et al., 2004; Hakeem et al., 2012; Vijayalakshmi et al., 2013). NupE is the ability of the plant to take up N from the soil and NutE is the ability to use N to produce gain yield (Ladha et al., 1998; Hirel et al., 2007). Studies on genotypic variations under low and recommended N for NUE traits at both seedling and maturity stage of plants under controlled and field conditions in various crops resulted in the identification of high NUE genotypes possessing high yield sustainability under low N condition (Muchow, 1998; Inthapanya et al., 2000; Le Gouis et al., 2000; Presterl et al., 2003; Anbessa et al., 2009; Namai et al., 2009; Vijayalakshmi et al., 2015). The high NUE/N-insensitive genotypes (NIS-top grain yielders) are defined as genotypes that give more or equal GY with minimal application of N fertilizer when compared to recommended or standard N fertilizer conditions (Hawkesford, 2017). The understanding of GY and its associated NUE traits performance is still largely lagging in the pearl millet diversity panel. However, few studies had explored genotypic variations for NUE at different N levels. In one such report, 20 diverse pearl millet genotypes and few high-yielding hybrids were screened under field conditions at two N levels (Alagarswamy and Bidinger, 1987). In another study, pearl millet hybrids were grown over 3 years in western and eastern Nebraska at four different N levels to determine the optimum N rate for cultivation and concluded that the NUtE trait was less responsive than the N uptake with increased N levels (Maman et al., 2006).
Nitrogen use efficiency, which is dependent on soil or external N supply, is an output of available N uptake, its efficient utilization and remobilization to grain at end of the season. Several studies revealed that improved NUtE might lead to enhanced NUE under low N conditions, particularly in cereals (Moll et al., 1982; Good et al., 2004). NUtE is defined as the genotype ability to assimilate and remobilize N which ultimately convert into GY. It is an essential physiological parameter that unveiled the positive relationship with GY. Whereas in pearl millet comprehensive studies are required to pinpoint the critical factors underlying NUtE under different N levels. GY is a complex trait controlled by a network of multiple traits and their associations. The GY in pearl millet is a result of many yield components, such as grain number, grain weight, tiller number, and panicle number (Rai et al., 2012; Basava et al., 2019). Hence, uncovering the genetic basis of GY, and other related NUE traits under low and high N conditions is a prerequisite to understand the mechanism and also to identify ideal NUE associated traits for selection. Thus, the current study was aimed to evaluate a diverse set of pearl millet lines for their response to GY and NUE related traits under different N levels and also to determine the suitable NUE traits for selection. This study provides useful information toward uncovering the physiological and genetic basis of GY and its related traits under low N which further facilitates the development of low N stress-resilient pearl millet cultivars.
Materials and Methods
Experimental Details
The experiment was conducted at an N-depleted precision field of the International Crops Research Institute for the Semi-Arid Tropics (ICRISAT), Patancheru, Telangana, India. The farm is geographically situated at an altitude of 545 m above mean sea level on 17.53°′ N latitude and 78.27°′ E longitude. The climate of the location is semi-arid with an average rainfall of 898 mm. The minimum and maximum temperatures ranged from 38 to 42°C which was observed during the experimental seasons. A set of 380 diverse pearl millet lines consisting of a world diversity panel, Pearl Millet Inbred Germplasm Association Panel-(PMiGAP) (Sehgal et al., 2015), mapping population parents and checks were screened during three seasons (summer 2017, rainy 2017, and summer 2018) in the same field with minimal fertility and moisture gradient. Across all the three seasons, field experiments were carried out in a split-plot alpha-lattice design with 2 m and three-row plots with two replications under three N treatments [N0 (without additional N application), N50 (@50 kg ha–), and N 100 (@100 kg ha–)]. Here N source used was urea (46.4% N). Field area with three different N treatments was considered as the main plot and it was divided into six blocks (two replications for each treatment). Each N level block was divided into 25 sub-blocks, and each sub-block consisted of 16 genotypes. Nitrogen depleted plot was developed and maintained at ICRISAT over past many seasons before the inception of the experiment by withholding N application. Before initiating the experiment, important physical and chemical properties of the soil were measured. Soil samples were collected at 15 cm depth from different places and four corners of the plot; mixed composite was used to determine the soil properties (Supplementary Table 1).
Crop Management
Crop management including land preparation, sowing/planting, management of water, nutrient, disease, and pest control were taken care of throughout the crop period. Seeds were sown on the raised furrows by using a four cone planter by maintaining sufficient spacing of 60 cm from row to row and 15 cm from plant to plant. The crop was fertilized by manual broadcasting as per the treatment. N fertilization was done as per the treatments using urea (46.5%) in two equal splits at 20–25 and 40–45 days after sowing. Besides, other nutrients such as single super phosphate @250 kg ha–1 and muriate of potash (60% K2O) @50 kg ha–1 was also applied at the time of sowing. All the pearl millet genotypes were evaluated for a total of 11 traits across three consecutive growing seasons under different N regimes. The measured phenotypic traits were (1) SPAD chlorophyll content (soil plant analytical device), (2) leaf number (LN), (3) LA, (4) GY, (5) dry stover yield (DSY), (6) N percent in dry stover (NPS), (7) N uptake in grain (NUpG), (8) N uptake in stover (NUpS), (9) total nitrogen uptake (TNUp), (10) NUtE, and (11) nitrogen harvest index (NHI). The phenotypic data were averaged across three seasons in each genotype for the first five traits SPAD chlorophyll content, LN, LA, GY, and DSY. For the remaining six NUE related traits data was presented as a pooled mean of two summer seasons, due to lack of quality data in the rainy season.
Physiological Traits
SPAD Chlorophyll Content (Soil Plant Analytical Device)
Soil plant analytical device chlorophyll content was measured by using Minolta Corporation’s Chlorophyll SPAD-502 plus, United States. Plants from all the treatments were marked a day before with different color ribbons. Fully expanded uppermost leaves at 45 days’ stage (from the day of leaf emergence) were selected and an average of three reads was recorded from a total of three plants from the middle row.
Leaf Area and Leaf Number
In order to examine the LA and LN, leaf samples were collected approximately 80 days after emergence for all the treatments and replications. Healthy plant per plot was harvested from the middle row in the field in early morning and on the same day, specific LA and number were measured by passing in all the collected leaves of the plant through LA meter LI-3100C (LI-COR, United States).
Grain and Stover Yield
At the crop harvesting stage, GY was measured after threshing sundried panicles to remove the grains from their central rachis. Harvested plants were selected from the middle row and expressed in grams (g) per 2-meter area. Similarly, DSY was also measured after drying the stover samples of all the middle row plants without panicles and expressed in grams per 2-meter area.
Nitrogen Estimation and Nitrogen Indices
Nitrogen Estimation
Nitrogen content in grain and stover samples was estimated by using the sulfuric acid-selenium digestion method. Grain and stover were made into fine powder by using clone mixture (cyclone sample mill) and 0.250 gm was used for nitrogen estimation in the Charles Renard Analytical Laboratory at ICRISAT, Patancheru. Pre-weighed samples were digested with the sulfuric acid-selenium and then analyzed (Sahrawat et al., 2002) by using an auto-analyzer (Skalar SAN System, AA Breda, Netherlands). Nitrogen concentration is expressed as Nitrogen% (N%).
Different Nitrogen Use Efficiencies
Total N uptake was calculated from the sum of grain uptake and stover N uptake values. Then different N efficiencies were calculated as per the formulas given by Fageria et al. (2010).
Nitrogen utilization efficiency (NUtE)
NUtE (kg grain kg–1 N) = GY/TNUp.
Nitrogen harvest index (NHI)
NHI (%) = (NUpG/TNUp) × 100.
Statistical Analysis
The phenotypic mean values over the replications of each genotype for all the traits were prepared and the best linear unbiased estimates (BLUE) for genotypes were estimated at each and across N levels using Genstat software, 20th edition. The minimum, maximum, and BLUE of each trait at each N level were calculated. ANOVA was performed for the individual and pooled seasons by modeling individual season’s residuals variance using Genstat software, 20th edition (VSN International, 2019). Fisher’s t-test was used to ascertain the significant difference among the genotypes, treatments and interactions. Pearson’s correlations coefficient analysis was performed by using pooled data of two summer seasons for all the 11 traits to identify the relationship between traits at each level of N treatments.
Cluster Analysis
Hierarchical cluster analysis was carried out using the Euclidean distance metric and UPGMA method (un-weighted paired group method and arithmetic averages)1. By using the pooled data of three seasons treatment wise (N0, N50, and N100) cluster analysis was performed for the five traits viz., SPAD chlorophyll content, LN, LA, GY, and DSY.
Results
According to the fertilizer map, the available nitrogen (N) content of the soil sample was divided in to low, medium, and high. Therefore, low N conditions were observed when soil nitrogen content is below 280 kg/ha (critical limits of the soil maintained as per ICRISAT guide limits). The selected experimental field was not used for any cultivation from the past 5–6 years and before starting the trial, soil samples were collected across the field and undergone soil nutrient analysis. Nutrient analysis (Supplementary Table 1) revealed 0.5% organic carbon (OC), indicating a low available N in the soil and the basal N availability was common across all the N regimes before the initiation of the experiment. In addition, the soil micronutrient availability was in the medium range. The basal N availability in the soil is common for all the treatments but coming to the N50 and N100 regimes, an additional N source in the form of urea was applied as per the protocol mentioned in the methodology section. Overall, the N0 regime was maintained strictly without any additional supply of nitrogen.
A total of 380 pearl millet genotypes including the diversity panel, parents of mapping populations, and checks were screened for 11 physio-agronomic and NUE traits at three nitrogen levels (0, 50, and 100% of the recommended N doses) with two replications. Data recorded were averaged across the three seasons for each trait in each N level. Results revealed the native variation across the genotypes toward N response which had given the scope to identify nitrogen use efficient lines under low and high N conditions for the marginal and favorable ecologies, respectively. ANOVA results revealed significant variations among the genotypes, treatments and their interactions (season × treatment (N levels), season × genotype, treatment × genotypes, and season × treatment × genotype). The pooled means over three seasons in each genotype in each N level for five physio agronomic traits were provided in Supplementary Table 2. The pooled mean of over two seasons for six N related traits in each genotypes were provided in Supplementary Table 3.
Physiological Traits
Key physiological traits viz., SPAD chlorophyll content, LN, and LA were recorded for all three seasons. Significant variations across the genotypes under different N regimes were noticed. Under N0, SPAD chlorophyll content was ranged from 28.18 to 49.11 with an average of 37.85. Whereas in N50, it was ranged from 32.16 to 50.1 with an average of 40.59. In N100, the values were ranged from 29.04 to 49.07 with an average of 41.77. In N0, SPAD chlorophyll content values were reduced by 6.8 and 9.4% compared to N50 and N100 conditions, respectively. LN was decreased under N0 conditions by 19.3 to 29.6% compared to N50 and N100 conditions, respectively. The range of LA was lower under N0 compared to N100 conditions. The mean of LA was reduced significantly under the N0 condition compared to the N100 condition. In addition, the analysis of variance (ANOVA) results revealed that the effect of genotype, treatments and their interactions were significant (P < 0.001) for all the three physiological traits studied. The data was averaged across the three seasons for each trait in each N level. Descriptive statistics over the three seasons, and ANOVA results for all the physiological traits were provided in Table 1.
Table 1. Summary statistics and ANOVA results from three seasons pooled data for physiological and agronomic traits under three N levels.
Agronomic Traits
Grain Yield
In N0, the average GY of the tested genotypes was 67.04 g and ranged from 9.64 g (IP13363) to 138.48 g (IP18621). Similarly, significant genotypic variations were observed in other treatments. In N50, GY values varied from 20.72 g (Tift238 D1) to 168.85 g (IP16096) with a mean of 92.61 g. Whereas in N100, GY ranged from 20.48 g (IP13363) to 164.87 g (ICMV-IS89305) with a mean of 88.69 g. Nitrogen in limited condition (N0) resulted in a 27.6 and 24.6% reduction in GY compared to N50, and N100, respectively. The ANOVA results indicated that the effect of genotype, treatments were significant (P < 0.001) and their interactions were also significant between seasons × treatment (P < 0.001), seasons × genotypes (P < 0.001), treatment × genotype (P < 0.001), and season × treatments × genotype.
Based on GY data, top 25 (N-insensitive (NIS-Top grain yielders) and least 25 (N-sensitive (NS-Poor grain yielders) genotypes were identified under N0 conditions (Table 2). Out of 25 NIS lines, nine genotypes (IP10820, IP17720, ICMB01222-P1, IP10379, ICMB89111-P2, IP8069, ICMB90111-P2, ICMV-IS89305, and ICMV221) were common in the top 25 lines at the N100 level which shows the genotype plasticity toward N0 and N100 conditions (Figure 1). Similarly, the top 25 high yielding genotypes were identified in N100. The 25 least grain yielding genotypes [N-sensitive (NS)] in N0 and N100 conditions were presented in Figure 2.
Table 2. Top 25 NIS and NS lines under low nitrogen (N0) with contrasting NUE grain yield derived from the pooled mean over three seasons.
Figure 1. Top 25 high grain yielding genotypes under N0 and N100 conditions. Total 9 genotypes (in red color) viz. IP 10820, IP 17720, ICMB 01222-P1, IP 10379, ICMB 89111-P2, IP 8069, ICMB 90111-P2, ICMV IS89305, and ICMV 221 exhibited high grain yield in both N0 and N100 conditions.
Figure 2. Top 25 poor yielding genotypes for grain yield under N0 and N100 conditions.
Furthermore, agronomic data was provided for the top five lines viz., IP 18621, IP 10820, IP 17720, IPC 804, and ICMB01222-P1, out of selected 25 NIS lines (Supplementary Table 4). Interestingly, days to 50% flowering (DF50%) and plant height (PH) was more in all the tested genotypes at N0, indicating the genotype efficiency in terms of agronomic traits. These genotypes can be used as parents in the pearl millet breeding program to develop N efficient genotypes under low N input conditions.
Dry Stover Yield
In N0, DSY ranged from 58.47 g (LGD-1-B-10) to 806.19 g (IP8863) with a mean of 256.3 g. Whereas in N50, values ranged from 62.42 g (IP13363) to 740.78 g (IP11584) with a mean of 246.16 g. In N100, the average DW of 380 genotypes was 264.41 g, and values varied from 73.92 g (Tift23D2B1-P1-P2) to 678.87 g (IP8786). Application of N significantly increased DSY by 12.6% in N50 and 18.02% in N100 compared to the N0 condition. Significant differences were observed for dry stover weight among the genotypes (P < 0.001), and interactions were also significant between seasons × treatment (P < 0.00), seasons × genotypes (P < 0.001), treatment × genotype (P < 0.001), and season × treatments × genotype.
Interestingly, out of the top 25 high dry stover yielding lines, nine genotypes (IP14398, IP10579, IP15857, IP10394, IP17125, IP13608, IP6098, IP12020, and IP11584) were common in the top 25 lines at N100 condition (Figure 3). At both N levels, the top 25 poor-performing dry stover genotypes were also identified (Figure 4).
Figure 3. Top 25 high dry stover yielding (DSY) genotypes under N0 and N100 conditions. Total 9 genotypes (highlighted in red color) viz., IP 14398, IP 10579, IP 15857, IP 10394, IP 17125, IP 13608, IP 6098, IP 12020, and IP 11584 exhibited high dry stover yield (DSY) in both N0 and N100 conditions.
Figure 4. Least performing genotypes for dry stover yield (DSY) under N0 and N100 conditions.
In addition, three high yielding dual purpose genotypes (IP5560, IP15857, and IP17125) were identified in the top 25 GY and DSY genotypes under N0. Whereas in N100, only one common high yielding genotype (IP3106) was identified. Descriptive statistics over the three seasons, and ANOVA results for the agronomic traits were provided in Table 1.
N Content in Grain and Stover
Nitrogen percent in grain (NPG) varied from 1.18% (IC804) to 2.43% (Tift 186), with an average of 1.61% in N0 condition. Whereas in N50, NPG ranged from 1.19% (9444) to 2.43% (IP4942) with a mean of 1.67%. Under N limiting conditions NPG was reduced by 3.6 and 11.4% as compared with N50 and N100 conditions, respectively. In N50, an average mean of NPG decreased by 8.24% as compared to N100. The response of genotypes in the N100 condition varied significantly from 1.27% (IP10539) to 2.59% (IP4828) with a mean of 1.82%. Significant differences were observed in NPG among the genotypes (P < 0.001), treatments (P < 0.001), and interactions and are also significant between seasons × treatment (P < 0.00), seasons × genotypes (P < 0.001), and season × treatments × genotype except for treatment × genotype (0.677).
Genotypic variations found in NPS in N0 and the minimum and maximum value of NPS ranged from 0.42% (IP8786) to 1.08% (ICMB90111-P5) with a mean of 0.61%. The average NPS was more in N0 compared to N50 and N100 conditions. Significant variations were observed in all the interactions except season × treatment (0.6771) and treatment × genotypes (0.0102).
Nitrogen Uptake for Dry Stover
Genotypic variations were observed for NUpS across the N levels. The mean of all genotypes was 1.65 kgha–1 ranged from 3.34 kgha–1 (ICMB9533-P5) to 8.03 kgha–1 (IP5560) and decreased by 9.8 and 14.45% in N0 compared to N50 and N100, respectively. Whereas, in N50, the values varied from 0.51 kgha–1 (AIMP92901-S1-15-1-B-P3) to 8.03 kgha–1 (IP10543) with a mean of 1.83 kgha–1 and decreased by 5.67% when compared with N100. The minimum and maximum values of NUpS varied from 0.38 kgha–1 (IP155363) to 1.08 kgha–1 (IP11584) with a mean of 6.59 kgha–1 under N100 condition.
Total Nitrogen Uptake
Application of N fertilizer resulted in an increase of TNUp and also revealed significant genotypic variations among the genotypes with in the treatments. Under N0, the values ranged from 1.09 to 10.47 with an average of 3.03 kg ha–1 and decreased by 14.4 and 21% as compared with N50 and N100, respectively. Whereas in N50 and N100, values varied from 1.08 (IP11229) to 9.57 (843B) and 0.58 (IP13363) to 8.97 (IP11584) with a mean of 3.54 and 3.83 kg ha–1, respectively. The results of ANOVA for TNUp indicated significant effects for all the interactions except treatments.
Nitrogen Utilization Efficiency
Large variation was observed for NUtE of genotypes, ranged from 2.40 (IP3110) to 50.60 (IP3557) with an average of 30.81 under N0 and the efficiency was reduced as compared with N50 and N100 by 7.3 and 5.95%, respectively. Whereas in N50, values varied from 10.63 to 50.86 with a mean of 33.22, and a 7.57% reduction was observed as compared with the N100 condition. The application of N fertilizer (N100) resulted in a significant increase in NUtE. The minimum and maximum NUtE were observed in IP31110 (10.06) and IP10456 (56.25), respectively with an average of 32.76. Significant differences were observed among the genotypes (P < 0.001), treatments (P < 0.001), and seasons × treatments (P < 0.001), seasons × genotypes (P < 0.001), treatment × genotypes and interactions were also significant between season × treatments × genotype.
Nitrogen utilization efficiency under N0 was compared with NUtE under N100 to determine the genotypic efficiency and their responsiveness to N. Based on the NUtE, the genotypes were classified into four groups (Figure 5) viz., (1) N efficient non-responsive (NENR), (2) N efficient responsive (NER), (3) N responsive inefficient, and (4) N inefficient non-responsive (Rengel and Graham, 1995; Worku et al., 2007). In this study, NutE data recorded in 374 genotypes across different N levels. The average NUtE of 374 genotypes under N0 (30.81) and N100 condition (32.76) were considered as cut off for the identification of genotypic efficiency and responsiveness for N use. Under low N, above-average genotypes were considered as efficient and below-average genotypes were considered as inefficient. Similarly, in N100, above-average genotypes were considered as responders, and below-average genotypes are considered as non-responders. Overall, efficient genotypes are higher in the utilization of absorbed N over inefficient genotypes.
Figure 5. Relationship between genotypes performance of nitrogen utilization efficiency (NUtE) under low (N0) and recommended N (N100) conditions. The yellow line represents the mean of NUtE at N100 and the purple line represents the mean of NUtE at N0 (pooled data from two summer seasons 2017 and 2018).
Nitrogen Harvest Index (NHI)
The nitrogen harvest index varied from 5.46% (IP3110) to 74.52% (ICMB 95333-P5) with an average of 47.56% in N0 condition and a reduction of about 10.3 and 15.43% compared to N50 and N100 condition, respectively. Whereas in N50, NHI values ranged from 21.10% (IP10488) to 84.02% (IP5923) with a mean of 53.0%. The response of genotypes under the N100 condition varied significantly and ranged from 23.51% (IP12020) to 78.89% (IP6110) with a mean of 56.24%. Further, the NHI reduction under N0 was about 10.3 and 15.43% as compared with N50 and N100, respectively. In N50, a reduction of NHI about 5.46% was observed compared to the N100 condition. Significant differences were observed among the genotypes (P < 0.001), treatments (P < 0.001), and interactions were also significant between seasons × treatment (P < 0.001), seasons × genotypes (P < 0.001), and season × treatments × genotype except for the treatment × genotype (0.677) (Table 3).
Table 3. ANOVA for N-related traits under three nitrogen levels during Summer 2017 and 2018.
Phenotypic Correlation
Correlation coefficient analysis was performed to identify the interrelation among the traits for all the N levels. The complete list of correlation coefficient values among the traits in each N level is provided in Figure 6. GY showed a significant positive correlation with TNUp, NUtE, and NHI. The results revealed that the GY of pearl millet was increased through the selection of higher NUtE and NHI lines. Furthermore, GY negatively correlated with DSY, NPS, NPG, and NUpS at both N0 and N100 conditions. The results suggested that the higher GY was, the lower DSY, NPS, NPG, and NUpS. Further DSY showed a negative significant correlation with all the measured traits except NUpS and TNUp across N levels. Interestingly NUtE showed a positive correlation with GY and NHI across N levels. SPAD chlorophyll content was found to be positively correlated with all traits except LA, DSY, and NUpS across all the treatments. Further LA was negatively correlated with GY, NUtE, and NHI across the treatments.
Figure 6. Pearson correlation for traits underlying NUE in 380 global diverse pearl millet genotypes (PMiGAP) under three different nitrogen (N) regimes carried out in 2017/18 (data pooled from two summer seasons). (A) Correlations among eight traits at N0 condition (N0). (B) Correlations among eight traits at N50 condition. (C) Correlations among eight traits at N100 condition. ∗ indicates < 0.05 significance, ∗∗ indicates < 0.01 significance.
Cluster Analysis
Hierarchical cluster analysis of 380 pearl millet genotypes revealed large genotypic variation in five traits viz., SPAD chlorophyll content, LN LA, GY, and DSY. This analysis is useful for grouping genotypes for NUE. Therefore, clustering 380 genotypes based on five traits identified two to three major clusters. Based on this, genotypes were classified into high, medium, and low-performing genotypes. At both N levels, clusters I and II contained better-performing genotypes and also found the least dissimilarity among the genotypes. Interestingly few common genotypes viz., IP 15536, IP 6098, IP11961, IP8863, IP11378, IP 20679, and IP 12138 were identified in N0 and N100 levels in the first two clusters, implying that genotype plasticity toward N under N0 and N100 conditions (Figures 7, 8). Cluster III contained more number of medium performing genotypes that form a maximum number of sub clusters at both N levels. Finally, the least performing genotypes were present in Cluster IV. Overall, the maximum dissimilarity (from all these matrices) was observed between cluster I and IV.
Figure 7. Dendrogram of 380 pearl millet genotypes for five traits under low nitrogen conditions (N0) by the UPGMA method.
Figure 8. Dendrogram of 380 pearl millet genotypes for five traits under recommended conditions of nitrogen (N100) by the UPGMA method.
Discussion
Nitrogen (N) is an essential macronutrient required for plant growth and is often a limiting factor for crop yield. From the past 50 years of agriculture, augmented food production was attained by the extensive application of N fertilizer in combination with N-responsive cultivars across the globe. Excessive application of N fertilizer is becoming expensive which accounts for the loss of economic profit to the farmers along with negative impacts on the environment (Raun et al., 2002; Hawkesford and Griffiths, 2019). Hence to overcome this problem, a clear understanding of genotype behavior, identification and development of genotypes (without compromising GY) with high NUE under low N conditions is a paramount need for improving NUE.
Therefore, the present study aimed to identify the N-Sensitive (NS) and N-insensitive (NIS) genotypes and also to determine the traits regulating low N tolerance. However, accurate and reliable phenotyping under low N input is challenging and influenced by the genotype (G), environment (E), and G × E interactions (Chen et al., 2014; Rao et al., 2018). Moreover, there are limited studies available in pearl millet toward the identification of nitrogen insensitive genotypes under low N conditions at field level. Based on the current knowledge, this study presents the first report of genetic variations for NUE in the global association panel, Pearl Millet inbred Germplasm Association Panel (PMiGAP) (PMiGAP collected from a large set of diverse germplasm and breeding lines) and from mapping population parents. Notably, very few pearl millet breeding programs are targeting the development of low N tolerance traits which are a must for sustainable agriculture with minimal negative impacts on the environment. Components of NUE were studied in pearl millet hybrids at two N levels found that the efficient genotype (Souna B) had a 32% of higher NUE value than the N inefficient Indian genotype (BJ104) (Alagarswamy and Bidinger, 1987). Another study revealed that NUtE contributes more to genetic variation in NUE (Maman et al., 2006).
In this study, a total of 380 diverse pearl millet genotypes were characterized for SPAD chlorophyll content, LN, LA, GY, DSY, and NUE related traits at three different N levels under field conditions. Notably, diverse responses have been observed among the genotypes across N levels, despite a similar growth conditions and an equal amount of N fertilizer in a given N level. These observed genotypic variations purely reveal the genotype plasticity toward traits. This set of lines used in this study which are collected from different regions across the globe and these lines were used for various abiotic stress studies to identify the tolerant genotypes. Previous studies in rice, maize, wheat, and oilseed rape, etc., have established significant genetic variation for NUE-related traits with large germplasm panels, hybrids, open-pollinated and recombinant inbred line populations (Chen et al., 2014; Li et al., 2015; Vijayalakshmi et al., 2015; Ertiro et al., 2017; He et al., 2017; Rao et al., 2018). Generally, the results of ANOVA explained statistical differences among the treatments and genotypes and their interactions. In the current study, ANOVA results revealed significant variations among the genotypes and treatments for all measured traits except a few traits (Tables 1, 3). These results were indicating that the environment under different N levels was a crucial factor in explaining the genotypic variance of GY and its related traits, which are strongly influenced by the relative contribution of G × E interactions under field evaluation. The results concur with other reported field experiments in rice (Srikanth et al., 2016) and wheat (Sial et al., 2005; Belete et al., 2018).
In the present study, SPAD chlorophyll content significantly increased with nitrogen application, which might be sufficient availability of N in the leaf. Similar results of increased SPAD chlorophyll content by increasing N application was reported by Kitajima and Hogan (2003), Pramanik and Bera (2013). The relation between the slow or controlled release of N fertilizers supporting more N absorption and related physiological mechanisms has been well studied in different crops (Li et al., 2003; Long et al., 2013). LA of all the pearl millet genotypes was significantly increased with N application which might be attributed to translocation of N to leaves, which brings variation in plant architecture and leaf internal structure (Singh et al., 2004; Mhaskar et al., 2005). Reduction in LN under N0 condition was observed, due to deprivation of N, thereby reduced source size and hindering the plant development.
Across the year on average, GY ranged from 9.64 g to 159.71 g (N0), 20.72 g to 236.43 g (N50), and 20.48 g to 242.12 g (N100), indicating the genotypic variability in GY at different N regimes. In general, GY increased was correlated with the N application rate, which might be due to sufficient nitrogen availability. Recently similar kind of results was reported in a diverse set of foxtail millet genotypes (Bandyopadhyay et al., 2020). In wheat, GY and quality traits were evaluated by Šarčević et al. (2014) at low and normal N and found that GY was reduced by 10% under low N condition compared to normal condition. The current results also are in tune with the earlier reports as there was an increased GY with increase N levels. Hawkesford and Griffiths (2019) reported that split application of N is the best option to recover the maximum applied N in the form of harvested grain/increased GY and found a strong significant relationship between GY and physiological traits.
Remarkably, 25 high yielding genotypes were identified under the N0 condition and considered as N-insensitive (NIS). Out of 25, nine genotypes (IP10820, IP17720, ICMB01222-P1, IP10379, ICMB89111-P2, IP8069, ICMB90111-P2, ICMV-IS89305, and ICMV221 = ICMV88904) are common in the top performing 25 lines under N100 condition, which shows the genotype plasticity toward varying N levels in these genotypes. These selected genotypes are good genetic resources to breed for tolerance to low N conditions. Top 25 high yielding genotypes in N100 can be considered to have the best acceptances for cultivation where soils are fertile and when followed the ideal N levels for cultivation. These identified genotypes could be used to improve GY, along with higher NUE. Numerous studies have reported utilization of N efficient lines with enhanced GY in the farmer fields which may help to reduce fertilizer input as well as increase the profitability of farm operations (Würschum, 2012; Vijayalakshmi et al., 2015; He et al., 2017). Likewise, selected NIS lines with high NUE will certainly play a role in reducing environmental pollution and could increase economic profit to farmers.
Pearl millet is an ideal fodder crop for feeding livestock and has unique features like dry season crop, production environments having low nitrogen, high photosynthetic efficiency, and high dry forage capacity (Govintharaj et al., 2018). Several reports depict the dry fodder demand and it would require approximately 568 million tonnes across the globe by the year 2030 (Govintharaj et al., 2018). To tackle this problem, high fodder yielding genotypes need to be identified. Generally, most of the genotypes have lower dry matter production under low N conditions, hence identification of genotypes with high dry matter production under low N conditions is important for sustainable fodder production. In this framework, the top and least performing 25 genotypes for DSY were identified under N0 and N100 conditions (Table 4).
Table 4. Dry stover yield (DSY) for top and least 25 genotypes under N0 condition (derived from pooled data of three seasons).
These contrasting genotypes can be further exploited in the breeding program and molecular studies to develop varieties with high NUE for biomass. Irrespective of N, the common nine genotypes (IP14398, IP10579, IP15857, IP10394, IP17125, IP13608, IP6098, IP12020, and IP11584) were identified and can produce more dry matter production under both N0 and N100 conditions. These genotypes need further physiological and molecular characterization under low nitrogen conditions to understand the molecular basis for low nitrogen tolerance and biomass production. Overall, observed genetic variability in pearl millet genotypes are suited well for improving DSY, can be utilized in the breeding program to enhance livestock productivity. Studies on various crops to determine the relationship between nitrogen levels and biomass traits have shown that N fertilizer application increased the biomass related traits, including DSY (Chan-Navarrete et al., 2014). The present results were tuned with these observations as DSY was more in N100 condition than in N0 condition; the increasing quantity was up to the significant levels. Another study revealed wide genetic variations for biomass traits and stover nitrogen in pearl millet germplasm (Gupta et al., 2015). Few studies revealed that significant genetic variation present in pearl millet for biomass traits in dual-purpose hybrids and open-pollinated varieties (OPVs), populations, and top cross hybrids (Bidinger et al., 2010; Blümmel et al., 2010; Rai et al., 2012). More importantly from this study, dual-purpose genotypes were identified and these can be used as a reference set for the developing high grain and stover yielding lines.
In the present study, a total of 11 traits were recorded for all three seasons. However, due to heavy rain, we are unable to generate quality data for the rainy season. Hence to avoid misinterpretation of results/findings, data for only six NUE traits were pooled from two summer seasons and presented. Moreover, the NPG was slightly increased with the nitrogen application and genetic differences for the traits were also observed. These results are in tune with the report of Belete et al. (2018), who found the higher NPG with increased N levels. NPG in N0 was reduced by 3.6 and 11.4% compared to N50 and N100, respectively. Whereas, in N50 condition, the mean NPG decreased by 8.24% as compared to N100. These results are similar to the finding of Lopez-Bellido et al. (2004), Arduini et al. (2006) who found the genetic variations among wheat genotypes for grain N concentration and increased NPG with increased N levels. In the case of the NPS, the trend was contrasting to that of NPG. Herein, NPS was slightly more in N0 as compared with N50 and N100 but was not up to a significant level. These results concur with the earlier reports of Belete et al. (2018) who found higher NPS with increasing N levels.
Nitrogen utilization efficiency is defined as the genotype ability to assimilate and remobilize N ultimately to produce the GY (Rengel and Graham, 1995; Fageria et al., 2008). Determination of the genetic variations in NUtE is essential for the selection of efficient genotypes and can be used further in breeding programs to develop low N tolerant material. The concept of genotypes grouping is used widely in nutrient use efficiency (Fageria and Baligar, 2003). Based on NUtE efficiency data in N0 versus N100, genotypes are classified into four groups viz., N efficient non-responsive, N efficient responsive (NER), N responsive inefficient, and N inefficient non-responsive under N0 and N100 condition, respectively. In N0 NUtE data, a total of 198 efficient and 176 inefficient genotypes were identified and out of 198 efficient genotypes, 58 genotypes were falling in the first desirable group which is N efficient non-responsive. These genotypes are exhibiting a progressive performance under low N. This may enable breeders to develop efficient genotypes under low input environments for pearl millet breeding activities. The remaining 140 genotypes were falling in the next most desirable group, NER and these genotypes were exhibiting a progressive response to increased N availability. The NER genotypes identified from the present investigation could be the prospective targets for selection toward the genetic improvement of pearl millet for N utilization. Interestingly, in the identified 25 NIS genotypes, 20 genotypes which are falling under N efficient responsive group were showed progressive performance in terms of efficient and responsive use of nitrogen. Genotypes with more NUtE could produce high GY per unit of N consumption (Moll et al., 1982). Various studies on NUtE have already reported that breeding for efficient genotypes under low N could be achievable with high NUtE (Fageria, 2014; He et al., 2017). The third group, responsive inefficient genotypes, can be used in breeding programs. The rest of the genotypes fall into the fourth group and these are less desirable from the NUE point of view. Interestingly, few common genotypes were identified from the non-responsive and inefficient genotypes. Overall, efficient genotypes are higher in the utilization of absorbed nutrients than inefficient genotypes. Vast genotypic variations in NUtE have been reported under field/pot screening in various wheat genotypes and other crops by several researchers (Bouchet et al., 2014; Ma et al., 2015).
The nitrogen harvest index indicates the level of efficiency of plants to use acquired total N for grain formation (He et al., 2017). In the present study, significant genotypic variation was observed for NHI at both N0 and N100 conditions. NHI was positively correlated with NUtE and GY in all three N levels. These findings concur with previous studies in oilseed rape (Ulas et al., 2013; Stahl et al., 2016) and wheat (Monostori et al., 2017).
In the current study, the Pearson correlation coefficient analysis for the GY and NUE related traits indicated that NUtE trait has the largest contribution to the GY followed by NHI and TNUp across N levels, implying that GY improvement in pearl millet could be possible by selecting genotypes that attain a higher NUtE and NHI across N levels. Similarly, a significant correlation between GY and NUtE was reported in maize, wheat, and oilseed rape (He et al., 2017; Belete et al., 2018). Fageria et al. (2010) reported that strong associations between GY and NUE related traits could be a better option for GY improvements under limited N condition. Furthermore, GY was negatively correlated with NPG and NPS across N levels. Previous studies also reported a similar significant negative correlation between grain N concentration and GY (Sinebo et al., 2004). In the present study, a strong positive correlation of DSY with NUpS and TNUp across N levels was observed indicating that more DSY was with higher NUpS and TNUp, where it was negatively correlated with NPS, NUtE, and NHI. The inverse relationship between NPS and DSY was reported in rice (Subudhi et al., 2020). Identified genetic stocks will be further utilized to carry an in-depth investigation to understand the genes and associated pathways related to NUE. Overall, important correlations identified from this study, will certainly help for the future pearl millet NUE breeding programs.
Conclusion
Improving NUE of pearl millet is pivotal for sustainable crop growth and yield especially under low nitrogen soils. Likewise, improving crop productivity using N fertilization is important for achieving climate resilience. Nevertheless, the genetic improvement of pearl millet NUE depends on the nature and extent of variation among the germplasm. As studies on pearl millet NUE are still nascent, this study was aimed to derive morphologic and agronomical traits associated with NUE in a set of 380 diverse lines under different N levels. The first step in this research revealed extensive genetic variations across the diversity panel and mapping population parents under N0 and N100 conditions. Also, large environmental variations (including N inputs) for GY and its related NUE traits were observed. Nitrogen limitation resulted in the reduction of GY and DSY in N0, as compared with N50 and N100. Under low N, the best grain yielding genotypes were identified and considered as nitrogen-insensitive (NIS). Genotypes IP10820, IP17720, ICMB01222-P1, IP10379, ICMB89111-P2, IP8069, ICMB90111-P2, ICMV-IS89305, and ICMV221 (ICMV88904) proved to be the most efficient genotypes in terms of GY at low and high N levels, and indeed shows their inherent genotypic plasticity toward N application. Furthermore, the use of NIS genotypes will help in breeding N-efficient genotypes for arid and marginal agro-ecologies. Overall results suggest that genotypes with more yield and high to moderate NUtE can be used as parents for the breeding of N efficient genotypes. The lines identified from the present study, coupled with multi-omics technologies can help to identify candidate genes for NUE in pearl millet. The available genetic stocks will be useful to carry an in-depth dissection of the genes and pathways related to NUE. These contrasting genetic stocks may help to map the genome, transcriptome, proteome, and metabolome signatures of the traits leading to a better understanding of NUE in pearl millet. These current findings may also help the farmers for optimizing the use of fertilizer inputs for economic and environmentally sustainable food production.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.
Author Contributions
RKS and RG planned and designed this research. VP, MP, SB, RR, and VT performed the field experiments, phenotyping, and analyzed the data. AR, RD, and GA carried out statistical analysis. RKS and RG interpreted the data. VP, MP, RKS, and RG wrote the manuscript. All the authors have made their contribution in editing the manuscript for publication.
Funding
Funding support from the Department of Biotechnology (Grant# BT-IN-UK-VNC-42-RG-2015-16), Government of India, the Biotechnology and Biological Sciences Research Council (GCRF/BBSRC Grant # BB/P027970/1), United Kingdom, and the CGIAR Research Program on Grain Legumes and Dryland Cereals (CRP-GLDC) is gratefully acknowledged.
Conflict of Interest
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.
Acknowledgments
We are grateful to Nagaraj Goud N., Ravinder T., and Shyam U. for data collection at the experimental site. We are also thankful to soil science laboratory technicians for Nitrogen (N) analysis in different samples and tissues.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2021.625915/full#supplementary-material
Supplementary Table 1 | Details of soil properties before establishment of experimental plot.
Supplementary Table 2 | Phenotypic data of 380 pearl millet genotypes for five traits (SPAD, Leaf number, Leaf area, Grain yield, and Dry stover yield) at three N levels. Mean values of over three seasons (summer 2017, rainy 2017, and summer 2018).
Supplementary Table 3 | Phenotypic data of 380 pearl millet genotypes for six traits (1. N percent in grain 2. N percent in stover 3. N uptake in stover 4. Total N uptake 5. Nitrogen utilization efficiency 6. Nitrogen Harvest Index) at three N levels. Mean values of over two seasons (summer 2017 and summer 2018).
Supplementary Table 4 | Agronomy (days to 50% flowering and plant height) data provided for top five NIS genotypes under N0.
Footnotes
References
Alagarswamy, G., and Bidinger, F. R. (1987). “Genotypic variation in biomass production and nitrogen use efficiency in pearl millet [Pennisetum americanum (L.) Leeke],” in Genetic Aspects of Plant Mineral Nutrition, eds W. H. Gabelman and B. C. Loughman (Dordrecht: Springer), 281–286. doi: 10.1007/978-94-009-3581-5_25
Anbessa, Y., Juskiw, P., Good, A., Nyachiro, J., and Helm, J. (2009). Genetic variability in nitrogen use efficiency of spring barley. Crop Sci. 49, 1259–1269. doi: 10.2135/cropsci2008.09.0566
Arduini, I., Masoni, A., Ercoli, L., and Mariotti, M. (2006). Grain yield, and dry matter and nitrogen accumulation and remobilization in durum wheat as affected by variety and seeding rate. Eur. J. Agron. 25, 309–318. doi: 10.1016/j.eja.2006.06.009
Bandyopadhyay, T., Swarbreck, S. M., Jaiswal, V., Gupta, R., Bentley, A. R., Griffiths, H., et al. (2020). Grain number and genotype drive nitrogen-dependent yield response in the C4 model Setaria italica (L.) P. Beauv. bioRxiv [Preprint]. doi: 10.1101/2020.03.23.003004
Basava, R. K., Hash, C. T., Mahendrakar, M. D., Kishor, P. B. K., Satyavathi, C. T., Kumar, S., et al. (2019). Discerning combining ability loci for divergent environments using chromosome segment substitution lines (CSSLs) in pearl millet. PLoS One 14:e0218916. doi: 10.1371/journal.pone.0218916
Belete, F., Dechassa, N., Molla, A., and Tana, T. (2018). Effect of nitrogen fertilizer rates on grain yield and nitrogen uptake and use efficiency of bread wheat (Triticum aestivum L.) varieties on the Vertisols of central highlands of Ethiopia. Agric. Food Secur. 7, 1–12. doi: 10.1186/s40066-018-0231-z
Bidinger, F. R., Blümmel, M., Hash, C. T., and Choudhary, S. (2010). Genetic enhancement for superior food-feed traits in a pearl millet (Pennisetum glaucum (L.) R. Br.) variety by recurrent selection. Anim. Nutr. Feed Tech. 10, 61–68.
Blümmel, M., Khan, A. A., Vadez, V., Hash, C. T., and Rai, K. N. (2010). Variability in stover quality traits in commercial hybrids of pearl millet (Pennisetum glaucum (L.) R. Br.) and grain-stover trait relationships. Anim. Nutr. Feed Tech. 10, 29–38.
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
Chan-Navarrete, R., Kawai, A., Dolstra, O., Van Bueren, E. T. L., and Van der Linden, C. G. (2014). Genetic diversity for nitrogen use efficiency in spinach (Spinacia oleracea L.) cultivars using the Ingestad model on hydroponics. Euphytica 199, 155–166. doi: 10.1007/s10681-014-1186-1
Chen, B., Xu, K., Li, J., Li, F., Qiao, J., Li, H., et al. (2014). Evaluation of yield and agronomic traits and their genetic variation in 488 global collections of Brassica napus L. Genet. Resour. Crop Evol. 61, 979–999. doi: 10.1007/s10722-014-0091-8
D’Andrea, A. C., and Casey, J. (2002). Pearl millet and Kintampo subsistence. Afr. Archaeol. Rev. 19, 147–173. doi: 10.1023/A:1016518919072
Ertiro, B. T., Beyene, Y., Das, B., Mugo, S., Olsen, M., Oikeh, S., et al. (2017). Combining ability and testcross performance of drought-tolerant maize inbred lines under stress and non-stress environments in Kenya. Plant Breed. 136, 197–205. doi: 10.1111/pbr.12464
Fageria, N. K. (2014). Nitrogen harvest index and its association with crop yields. J. Plant Nutr. 37, 795–810. doi: 10.1080/01904167.2014.881855
Fageria, N. K., and Baligar, V. C. (2003). Methodology for evaluation of lowland rice genotypes for nitrogen use efficiency. J. Plant Nutr. 26, 1315–1333. doi: 10.1081/PLN-120020373
Fageria, N. K., Baligar, V. C., and Li, Y. C. (2008). The role of nutrient efficient plants in improving crop yields in the twenty first century. J. Plant Nutr. 31, 1121–1157. doi: 10.1080/01904160802116068
Fageria, N. K., De Morais, O. P., and Dos Santos, A. B. (2010). Nitrogen use efficiency in upland rice genotypes. J. Plant Nutr. 33, 1696–1711. doi: 10.1080/01904167.2010.496892
FAO (2019). FAO World Fertilizer Trends and Outlook to 2020. Available online at: http://www.fao.org/3/a-i6895e.pdf (accessed June 11, 2019).
Gianquinto, G., Goffart, J. P., Olivier, M., Guarda, G., Colauzzi, M., Dalla Costa, L., et al. (2004). The use of hand-held chlorophyll meters as a tool to assess the nitrogen status and to guide nitrogen fertilization of potato crop. Potato Res. 47, 35–80. doi: 10.1007/bf02731970
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
Govintharaj, P., Gupta, S. K., Maheswaran, M., Sumathi, P., and Atkari, D. G. (2018). Correlation and path coefficient analysis of biomass yield and quality traits in forage type hybrid parents of pearl millet. Int. J. Pure Appl. Biosci. 6, 1056–1061. doi: 10.18782/2320-7051.5992
Gupta, S. K., Nepolean, T., Sankar, S. M., Rathore, A., Das, R. R., Rai, K. N., et al. (2015). Patterns of molecular diversity in current and previously developed hybrid parents of pearl millet [Pennisetum glaucum (L.) R. Br.]. Am. J. Plant Sci. 6, 1697–1712. doi: 10.4236/ajps.2015.611169
Hakeem, K. R., Chandna, R., Ahmad, A., Qureshi, M. I., and Iqbal, M. (2012). Proteomic analysis for low and high nitrogen-responsive proteins in the leaves of rice genotypes grown at three nitrogen levels. Appl. Biochem. 168, 834–850. doi: 10.1007/s12010-012-9823-4
Hawkesford, M. J. (2017). Genetic variation in traits for nitrogen use efficiency in wheat. J. Exp. Bot. 68, 2627–2632. doi: 10.1093/jxb/erx079
Hawkesford, M. J., and Griffiths, S. (2019). Exploiting genetic variation in nitrogen use efficiency for cereal crop improvement. Curr. Opin. 49, 35–42. doi: 10.1016/j.pbi.2019.05.003
He, H., Yang, R., Li, Y., Ma, A., Cao, L., Wu, X., et al. (2017). Genotypic variation in nitrogen utilization efficiency of oilseed rape (Brassica napus) under contrasting N supply in pot and field experiments. Front. Plant Sci. 8:1825. doi: 10.3389/fpls.2017.01825
Hickman, J. E., Palm, C. A., Mutuo, P., Melillo, J. M., and Tang, J. (2014). Nitrous oxide (N 2 O) emissions in response to increasing fertilizer addition in maize (Zea mays L.) agriculture in western Kenya. Nutr. Cycling Agroecosyst. 100, 177–187. doi: 10.1007/s10705-014-9636-7
Hirel, B., Le Gouis, J., Ney, B., and Gallais, A. (2007). The challenge of improving nitrogen use efficiency in crop plants: towards a more central role for genetic variability and quantitative genetics within integrated approaches. J. Exp. Bot. 58, 2369–2387. doi: 10.1093/jxb/erm097
Inthapanya, P., Sihavong, P., Sihathep, V., Chanphengsay, M., Fukai, S., and Basnayake, J. (2000). Genotype differences in nutrient uptake and utilisation for grain yield production of rainfed lowland rice under fertilized and non-fertilized conditions. Field Crops Res. 65, 57–68. doi: 10.1016/S0378-4290(99)00070-2
Kiran, T. V., Vijayalakshmi, P., Rao, Y. V., Swamy, K. N., Kondamudi, R., Srikanth, B., et al. (2016). Effects of nitrogen limitation on antioxidant enzymes, chlorophyll content and grain yield of rice genotypes. Asian Res. J. Agri. 1, 1–10. doi: 10.9734/arja/2016/28503
Kitajima, K., and Hogan, K. P. (2003). Increases of chlorophyll a/b ratios during acclimation of tropical woody seedlings to nitrogen limitation and high light. Plant Cell Environ. 26, 857–865. doi: 10.1046/j.1365-3040.2003.01017.x
Ladha, J. K., Kirk, G. J. D., Bennett, J., Peng, S., Reddy, C. K., Reddy, P. M., et al. (1998). Opportunities for increased nitrogen-use efficiency from improved lowland rice germplasm. Field Crops Res. 56, 41–71. doi: 10.1016/S0378-4290(97)00123-8
Le Gouis, J., Béghin, D., Heumez, E., and Pluchard, P. (2000). Genetic differences for nitrogen uptake and nitrogen utilization efficiencies in winter wheat. Eur. J. Agron. 12, 163–173. doi: 10.1016/S1161-0301(00)00045-9
Li, F., Ai, T., Zhou, S., Nie, J., and Liu, F. (2003). Influence of slow-release nitrogen fertilizers on lowland rice yield and nitrogen use efficiency. Chin. J. Soil Sci. 35, 311–315.
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
Long, J. R., Wan, Y. Z., Song, C. F., Jian, S. U. N., and Qin, R. J. (2013). Effects of nitrogen fertilizer level on chlorophyll fluorescence characteristics in flag leaf of super hybrid rice at late growth stage. Rice Sci. 20, 220–228. doi: 10.1016/s1672-6308(13)60138-9
Lopez-Bellido, R. J., Shepherd, C. E., and Barraclough, P. B. (2004). Predicting post-anthesis N requirements of bread wheat with a Minolta SPAD meter. Eur. J. Agron. 20, 313–320. doi: 10.1016/S1161-0301(03)00025-X
Ma, B. L., Biswas, D. K., Herath, A. W., Whalen, J. K., Ruan, S. Q., Caldwell, C., et al. (2015). Growth, yield, and yield components of canola as affected by nitrogen, sulfur, and boron application. J. Plant Nutr. Soil Sci. 178, 658–670. doi: 10.1002/jpln.201400280
Maman, N., Mason, S. C., and Lyon, D. J. (2006). Nitrogen rate influence on pearl millet yield, nitrogen uptake, and nitrogen use efficiency in Nebraska. Commun. Soil Sci. Plant Anal. 37, 127–141. doi: 10.1080/00103620500406112
Manning, K., Pelling, R., Higham, T., Schwenniger, J. L., and Fuller, D. Q. (2011). 4500-Year old domesticated pearl millet (Pennisetum glaucum) from the Tilemsi Valley, Mali: new insights into an alternative cereal domestication pathway. J. Archaeol. Sci. 38, 312–322. doi: 10.1038/nrg2612
Mhaskar, N. V., Dongale, J. H., Dademal, A. A., and Khanvilkar, S. A. (2005). Performance of scented rice varieties under different nitrogen levels in lateritic soil of Konkan. Oryza 42:323.
Moll, R. H., Kamprath, E. J., and Jackson, W. A. (1982). Analysis and interpretation of factors which contribute to efficiency of nitrogen utilization. J. Agron. 74, 562–564. doi: 10.2134/agronj1982.00021962007400030037x
Monostori, I., Árendás, T., Hoffman, B., Galiba, G., Gierczik, K., Szira, F., et al. (2016). Relationship between SPAD value and grain yield can be affected by cultivar, environment and soil nitrogen content in wheat. Euphytica 211, 103–112. doi: 10.1007/s10681-016-1741-z
Monostori, I., Szira, F., Tondelli, A., Arendas, T., Gierczik, K., Cattivelli, L., et al. (2017). Genome-wide association study and genetic diversity analysis on nitrogen use efficiency in a Central European winter wheat (Triticum aestivum L.) collection. PLoS One 12:e0189265. doi: 10.1371/journal.pone.0189265
Muchow, R. C. (1998). Nitrogen utilization efficiency in maize and grain sorghum. Field Crops Res. 56, 209–216. doi: 10.1016/S0378-4290(97)00132-9
Namai, S., Toriyama, K., and Fukuta, Y. (2009). Genetic variations in dry matter production and physiological nitrogen use efficiency in rice (Oryza sativa L.) varieties. Breed. Sci. 59, 269–276. doi: 10.1270/jsbbs.59.269
Pramanik, K., and Bera, A. K. (2013). Effect of seedling age and nitrogen fertilizer on growth, chlorophyll content, yield and economics of hybrid rice (Oryza sativa L.). Int. J. Agron. Plant Prod. 4, 3489–3499.
Presterl, T., Seitz, G., Landbeck, M., Thiemt, E. M., Schmidt, W., and Geiger, H. H. (2003). Improving nitrogen-use efficiency in European maize. Crop Sci. 43, 1259–1265. doi: 10.2135/cropsci2003.1259
Rai, K. N., Blummel, M., Singh, A. K., and Rao, A. S. (2012). Variability and relationships among forage yield and quality traits in pearl millet. Eur. J. Plant Sci. Biotechnol. 6, 118–124.
Rao, I. S., Neeraja, C. N., Srikanth, B., Subrahmanyam, D., Swamy, K. N., Rajesh, K., et al. (2018). Identification of rice landraces with promising yield and the associated genomic regions under low nitrogen. Sci. Rep. 8:9200. doi: 10.1038/s41598-018-27484-0
Raun, W. R., Johnson, G. V., and Westerman, R. L. (1999). Fertilizer nitrogen recovery in long-term continuous winter wheat. Soil Sci. Soc. Am. J. 63, 645–650. doi: 10.2136/sssaj1999.03615995006300030030x
Raun, W. R., Solie, J. B., Johnson, G. V., Stone, M. L., Mullen, R. W., Freeman, K., et al. (2002). Improving nitrogen use efficiency in cereal grain production with optical sensing and variable rate application. J. Agron. 94, 815–820. doi: 10.2134/agronj2002.8150
Rengel, Z., and Graham, R. D. (1995). Wheat genotypes differ in Zn efficiency when grown in chelate-buffered nutrient solution. Plant Soil 176, 307–316. doi: 10.1007/BF00011795
Russo, T. A., Tully, K., Palm, C., and Neill, C. (2017). Leaching losses from Kenyan maize cropland receiving different rates of nitrogen fertilizer. Nutr. Cycling Agroecosyst. 108, 195–209. doi: 10.1007/s10705-017-9852-z
Sahrawat, K. L., Kumar, G. R., and Murthy, K. V. S. (2002). Sulfuric acid–selenium digestion for multi-element analysis in a single plant digest. Commun. Soil Sci. Plan 33, 3757–3765. doi: 10.1081/CSS-120015920
Šarčević, H., Jukić, K., Ikić, I., and Lovrić, A. (2014). Estimation of quantitative genetic parameters for grain yield and quality in winter wheat under high and low nitrogen fertilization. Euphytica 199, 57–67. doi: 10.1007/s10681-014-1154-9
Sehgal, D., Skot, L., Singh, R., Srivastava, R. K., Das, S. P., and Taunk, J. (2015). Exploring potential of pearl millet germplasm association panel for association mapping of drought tolerance traits. PLoS One 10:e0122165. doi: 10.1371/journal.pone.0122165
Sial, M. A., Arain, M. A., Khanzada, S. H. A. M. A. D. A. D., Naqvi, M. H., Dahot, M. U., and Nizamani, N. A. (2005). Yield and quality parameters of wheat genotypes as affected by sowing dates and high temperature stress. Pak. J. Bot. 37, 575–584.
Sinebo, W., Gretzmacher, R., and Edelbauer, A. (2004). Genotypic variation for nitrogen use efficiency in Ethiopian barley. Field Crops Res. 85, 43–60. doi: 10.1016/S0378-4290(03)00135-7
Singh, T., Shivay, Y. S., and Singh, S. (2004). Effect of date of transplanting and nitrogen on productivity and nitrogen use indices in hybrid and non-hybrid aromatic rice. Acta Agron. Hung. 52, 245–252. doi: 10.1556/aagr.52.2004.3.5
Srikanth, B., Rao, I. S., Surekha, K., Subrahmanyam, D., Voleti, S. R., and Neeraja, C. N. (2016). Enhanced expression of OsSPL14 gene and its association with yield components in rice (Oryza sativa) under low nitrogen conditions. Gene 576, 441–450. doi: 10.1016/j.gene.2015.10.062
Srivastava, R., Singh, R. B., Lakshmi, P. V., Srikanth, B., Madhu, P., Satyavathi, C. T., et al. (2019). Genome-wide association studies (GWAS) and genomic selection (GS) in pearl millet: advances and prospects. Front. Genet. 10:1389. doi: 10.3389/fgene.2019.01389
Stahl, A., Friedt, W., Wittkop, B., and Snowdon, R. J. (2016). Complementary diversity for nitrogen uptake and utilisation efficiency reveals broad potential for increased sustainability of oilseed rape production. Plant Soil 400, 245–262. doi: 10.1007/s11104-015-2726-8
Subudhi, H. N., Prasad, K. V. S. V., Ramakrishna, C., Rameswar, P. S., Pathak, H., Ravi, D., et al. (2020). Genetic variation for grain yield, straw yield and straw quality traits in 132 diverse rice varieties released for different ecologies such as upland, lowland, irrigated and salinity prone areas in India. Field Crops Res. 245:107626. doi: 10.1016/j.fcr.2019.107626
Tako, E., Reed, S. M., Budiman, J., Hart, J. J., and Glahn, R. P. (2015). Higher iron pearl millet (Pennisetum glaucum L.) 967 provides more absorbable iron that is limited by increased polyphenolic content. Nutr. J. 14, 1–9. doi: 10.1186/1475-2891-14-11
Ulas, A., Behrens, T., Wiesler, F., and Horst, W. J. (2013). Does genotypic variation in nitrogen remobilisation efficiency contribute to nitrogen efficiency of winter oilseed-rape cultivars (Brassica napus L.). Plant Soil 371, 463–471. doi: 10.1007/s11104-013-1688-y
Vijayalakshmi, P., Kiran, T. V., Rao, Y. V., Srikanth, B., Rao, I. S., Sailaja, B., et al. (2013). Physiological approaches for increasing nitrogen use efficiency in rice. Indian J. Plant Physiol. 18, 208–222. doi: 10.1007/s40502-013-0042-y
Vijayalakshmi, P., Vishnukiran, T., Kumari, B. R., Srikanth, B., Rao, I. S., Swamy, K. N., et al. (2015). Biochemical and physiological characterization for nitrogen use efficiency in aromatic rice genotypes. Field Crops Res. 179, 132–143. doi: 10.1016/j.fcr.2015.04.012
Worku, M., Bänziger, M., Erley, G. S. A. M., Friesen, D., Diallo, A. O., and Horst, W. J. (2007). Nitrogen uptake and utilization in contrasting nitrogen efficient tropical maize hybrids. Crop Sci. 47, 519–528. doi: 10.2135/cropsci2005.05.0070
Würschum, T. (2012). Mapping QTL for agronomic traits in breeding populations. Theor. Appl. Genet. 125, 201–210. doi: 10.1007/s00122-012-1887-6
Keywords: pearl millet, nitrogen use efficiency, genotypic variations, phenotyping, nitrogen-responsive
Citation: Pujarula V, Pusuluri M, Bollam S, Das RR, Ratnala R, Adapala G, Thuraga V, Rathore A, Srivastava RK and Gupta R (2021) Genetic Variation for Nitrogen Use Efficiency Traits in Global Diversity Panel and Parents of Mapping Populations in Pearl Millet. Front. Plant Sci. 12:625915. doi: 10.3389/fpls.2021.625915
Received: 04 November 2020; Accepted: 11 January 2021;
Published: 04 February 2021.
Edited by:
Velu Govindan, International Maize and Wheat Improvement Center, MexicoReviewed by:
Chandirakala Ramakrishnan, Tamil Nadu Agricultural University, IndiaMehanathan Muthamilarasan, University of Hyderabad, India
Copyright © 2021 Pujarula, Pusuluri, Bollam, Das, Ratnala, Adapala, Thuraga, Rathore, Srivastava and Gupta. 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) and the copyright owner(s) 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: Rakesh K. Srivastava, ci5rLnNyaXZhc3RhdmFAY2dpYXIub3Jn; Rajeev Gupta, Zy5yYWplZXZAY2dpYXIub3Jn