Next Article in Journal
Recommended vs. Practice: Smallholder Fertilizer Decisions in Central Myanmar
Next Article in Special Issue
Genome-Wide Association Study Using a Multiparent Advanced Generation Intercross (MAGIC) Population Identified QTLs and Candidate Genes to Predict Shoot and Grain Zinc Contents in Rice
Previous Article in Journal
Research on Measurement Method of Leaf Length and Width Based on Point Cloud
Previous Article in Special Issue
Epistasis and Quantitative Resistance to Pyricularia oryzae Revealed by GWAS in Advanced Rice Breeding Populations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Drought Tolerant near Isogenic Lines (NILs) of Pusa 44 Developed through Marker Assisted Introgression of qDTY2.1 and qDTY3.1 Enhances Yield under Reproductive Stage Drought Stress

1
Division of Genetics, ICAR-Indian Agricultural Research Institute, New Delhi 110 012, India
2
Amity Institute of Organic Agriculture (AIOA), Amity University, Noida 201303, Uttar Pradesh, India
3
Division of Plant Physiology, ICAR-Indian Agricultural Research Institute, New Delhi 110 012, India
4
Rice Breeding and Genetics Research Centre, ICAR- Indian Agricultural Research Institute, Aduthurai 612 101, Tamil Nadu, India
5
Central Rainfed Upland Rice Research Station, ICAR-National Rice Research Institute, Post Box 48, Hazaribagh 825 301, Jharkhand, India
6
College of Agriculture, Jawaharlal Nehru Krishi Vishwa Vidyalaya, Rewa 486 004, Madhya Pradesh, India
*
Author to whom correspondence should be addressed.
Submission received: 15 December 2020 / Revised: 8 January 2021 / Accepted: 9 January 2021 / Published: 14 January 2021
(This article belongs to the Special Issue Rice Breeding and Genetics)

Abstract

:
Reproductive stage drought stress (RSDS) is detrimental for rice, which affects its productivity as well as grain quality. In the present study, we introgressed two major quantitative trait loci (QTLs), namely, qDTY2.1 and qDTY3.1, governing RSDS tolerance in a popular high yielding non-aromatic rice cultivar, Pusa 44, through marker-assisted backcross breeding (MABB). Pusa 44 is highly sensitive to RSDS, which restricts its cultivation across drought-prone environments. Foreground selection was carried out using markers, RM520 for qDTY3.1 and RM 521 for qDTY2.1. Background selection was achieved with 97 polymorphic SSR markers in tandem with phenotypic selection to achieve faster recurrent parent genome (RPG) recovery. Three successive backcrosses followed by three selfings aided RPG recoveries of 98.6% to 99.4% among 31 near isogenic lines (NILs). Fourteen NILs were found to be significantly superior in yield and grain quality under RSDS with higher drought tolerance efficiency (DTE) than Pusa 44. Among these, the evaluation of two promising NILs in the multilocational trial during Kharif 2019 showed that they were significantly superior to Pusa 44 under reproductive stage drought stress, while performing on par with Pusa 44 under normal irrigated conditions. These di-QTL pyramided drought-tolerant NILs are in the final stages of testing the All India Coordinated Rice Improvement Project varietal trials for cultivar release. Alternately, the elite drought-tolerant Pusa 44 NILs will serve as an invaluable source of drought tolerance in rice improvement.

1. Introduction

Rice (Oryza sativa L.) is an important food crop of Asian countries, including India. Globally, rice is widely consumed as a staple food by more than half of the human population, contributing over 20 percent of the total calorie intake. To meet the growing food demand, rice production has to be boosted by 0.6–0.9% annually [1]. Among the rice-growing countries of the world, India holds the first position in acreage (44 mha), but is second in rice production. India witnessed a 226% increase in rice production with a corresponding 28.3% increase in area, and a 152% increase in productivity during the green revolution, mainly due to the substantial increase in rice productivity in India through the cultivation of modern semi-dwarf high-yielding rice cultivars. However, the average annual productivity is still low compared to the world average. One of the major factors for lower productivity is the uncertain water availability in major rice-growing regions, followed by biotic stresses.
During the past sixty years, over 1200 rice varieties have been released for commercial cultivation; however, only a few of them have gained large-scale popularity among the farmers and other stakeholders associated with rice. One such variety is Pusa 44, developed and released in 1994 by the ICAR-Indian Agricultural Research Institute (ICAR-IARI), New Delhi. Even though Pusa 44 was released for commercial cultivation in Karnataka and Kerala, this variety became immensely popular elsewhere, growing into a megavariety in northern India [2], particularly due to its higher productivity and suitability for mechanical harvesting. Pusa 44 is a semi-dwarf indica rice variety, possessing sturdy culm, with long slender grains and high head rice recovery. With average productivity of 8–10 tonnes/hectare, Pusa 44 became popular in Punjab and quickly spread across to neighboring states of Haryana and Uttar Pradesh. Pusa 44 contributes a significant share in supply towards the public distribution system (PDS) from Punjab, and also contributes to the non-Basmati rice export from Northern India [3].
Water scarcity for agriculture purposes is a growing concern in the contemporary world. The increasing demand for available water is exacerbating water shortage in irrigated ecosystems, predisposing drought-like situations while bolstering the intensity of the prevailing drought in rainfed ecosystems [4]. Rice is vulnerable to drought stress, and drought has been predicted to occur more frequently due to climate change. Drought stress is detrimental to rice production and yield stability. Therefore, there is a need to develop and disseminate ecosystem-specific technologies that can reduce the excessive use of groundwater for rice production without compromising the productivity. To enable this, breeding rice varieties with resilience to drought is an important step forward in addressing this challenge. One of the major bottlenecks in analyzing the drought-related response in rice genotypes is the inherent complexity of the experimental system. This is because the drought itself is manifested in several ways, and is often complicated with the presence of high temperatures [5]. Additionally, the quick drought progression also affects the chance of crop survivability, as the plants get relatively less time to recover against the rapid damage induced due to drought stress. Screening for drought tolerance under artificially managed stress levels is comparatively hard compared to managing stress during the dry season, as well as under rainfed upland situations [6,7,8,9]. Therefore, screening for drought tolerance is carried out under managed conditions where facilities for controlling irrigations are available. Additionally, rainout shelters are also used to prevent precipitation interference in the experimentation [2]. In recent times, northern India has been facing issues of receding water levels bringing in episodes of intermittent drought when rainfall patterns fluctuate. Since Pusa 44 is highly susceptible to drought, especially at the reproductive stage, its yield and grain quality are compromised when exposed to reproductive stage drought stress (RSDS). This has raised concerns on the continuity of Pusa 44 cultivation in its niche areas of adoption. On account of being a lowland cultivar suitable for the irrigated rice ecosystem, it is obvious that the productivity of Pusa 44 can be significantly reduced due to water stress.
Genetically, drought tolerance in rice is a complex response. Several large-effect and small-effect QTLs and meta-QTLs have been reported for grain yield under reproductive-stage drought stress [10,11]. A large-effect QTL, qDTY2.2, on chromosome 2 was identified in the MTU 1010/Kali Aus population under upland [12] and lowland [13] conditions. Mishra et al. [14] reported another QTL, qDTY12.1, on chromosome 12 that was significantly associated with reproductive-stage drought stress. qDTY12.1 explains 42% of the genetic variation with an additive effect of 172 kg/ha for grain yield under drought with no significant penalty under unstressed conditions [8]. Another major-effect QTL, qDTY1.1 was identified on chromosome 1 between the marker intervals RM431 and RM12091 in the Swarna/IR 64 population [15]. Venuprasad et al. [16] identified two major QTLs, qDTY2.1 and qDTY3.1, on chromosome 2 and 3, respectively, from the population of Apo/2*Swarna, that explained 13–16% and 31% of the genetic variances, respectively. The subsequent use of these QTLs, particularly of qDTY3.1, has been demonstrated to impart improved drought tolerance among several genetic backgrounds, such as improved White Ponni [17], MR219 [18], Samba Mahsuri [19], Sabitri [20], IR64 and Vandana [21], and several others [22].
Marker-assisted selection (MAS) approach to transferring these QTLs into drought-sensitive genetic backgrounds of high-yielding varieties has been widely adopted in rice [22]. Jongdee et al. [23] advocated the use of molecular breeding to accelerate product development, especially for characters which are complex and have low heritability, or if the breeding process is time-consuming and expensive. The marker-assisted incorporation of major- and minor-effect QTLs for grain yield under drought stress has been recognized as time- and cost-effective, as well as a fast-track approach for breeding drought-tolerant rice varieties [24]. Dixit et al. [25] developed drought-tolerant near isogenic lines (NILs) of Savitri (CR1009) with the combination of two large-effect QTLs, qDTY3.2 and qDTY12.1. Marker-assisted backcross breeding (MABB) has been successful in transferring the genes/QTLs governing tolerance to biotic/abiotic stresses in rice [26,27,28,29,30,31,32,33,34,35].
Considering the popularity of Pusa 44 and its high sensitivity to RSDS, the present study aimed at the introgression of two QTLs, qDTY2.1 and qDTY3.1, governing tolerance to RSDS into Pusa 44 through MABB (Supplementary Figure S1). The improved drought tolerant NILs with tolerance to RSDS were tested for two years on the station and for another two years in multilocation trials, the results of which are discussed below.

2. Materials and Methods

2.1. Plant Material

Pusa 44 used as the recurrent parent and a drought-tolerant NIL of Swarna, namely, IR81896-B-B-142 was used as a donor for the two QTLs, qDTY2.1 and qDTY3.1, for tolerance to RSDS. Pusa 44 was crossed with the donor and the F1 progenies were tested for hybridity. The true F1 was backcrossed to the recurrent parent to generate BC1 generation. Two more backcrosses were made in subsequent generations using selected progenies based on marker-assisted foreground selection for the two QTLs, qDTY2.1 and qDTY3.1, followed by background and phenotypic selection to generate BC3 generation. Further selections on the segregating populations of BC3 were used to develop NILs of Pusa 44 possessing these two QTLs. All the early generations were developed by growing under irrigated conditions, with recommended agronomic maintenance. However, phenotypic selection was carried out starting with BC1F1, for various agro-morphological traits in comparison with Pusa 44. All the field experiments were conducted at the research farm of Division of Genetics, ICAR-Indian Agricultural Research Institute (ICAR-IARI), New Delhi (28°38′ N; 77°10′ E; 223 m AMSL), and the offseason crops were taken up at the Rice Breeding and Genetics Research Centre (RBGRC)- IARI, Aduthurai (11°00′ N; 79°28′ E; 19 m AMSL).

2.2. Marker-Assisted Selection

For marker analysis, DNA was extracted from leaf samples collected from field-grown plants as per the standard procedure [36]. PCR amplification was conducted as per the protocol standardized in our laboratory [33]. Foreground selection was carried out with QTL-linked markers, RM520 for qDTY3.1 and RM 521 for qDTY2.1 (Table 1). The parental lines were screened for the genome-wide polymorphism using 460 microsatellite (SSR) markers distributed across the rice genome. Among these polymorphic markers (Supplementary Figure S2), SSR markers on the target chromosomes 2 and 3 flanking the QTL linked markers and were used for effecting recombinant selection. After the generation of F1, the hybridity of the F1s was confirmed with the foreground markers by confirming the heterozygosity for the QTL-linked markers. One of the true F1 plants was backcrossed to Pusa 44 to generate BC1F1 plants. The BC1F1 plants were initially subjected to foreground selection using RM520, and the heterozygous plants thus identified were subsequently genotyped using RM521 to identify plants heterozygous for both the markers. Such double heterozygotic BC1F1s identified from foreground selection were screened for the background recovery using the genome-wide polymorphic SSR markers. The recovery of Pusa 44-specific alleles in these polymorphic markers was used for computing the recovery of the recurrent parent genome (RPG). Additionally, the selected BC1F1s were also screened for morphological similarity to Pusa 44, for identifying plants with higher RPG and recurrent parent phenome (RPP). One BC1F1 plant with desirable RPG and RPP recovery was used to generate BC2F1, by backcrossing it with Pusa 44. The selection strategy followed earlier in BC1F1 was adopted until the BC3F1 was generated. From BC3F1 onwards, the plants were advanced to BC3F4 to identify genotypes which are homozygous for both the target QTLs. Schematic selection (Figure 1) for agronomical, morphological, grain and cooking quality traits, coupled with the genotypic selection at every selection stage, was done to ensure maximum recovery of RPP. In each step, background SSR markers that showed non-recovery/heterozygosity for recurrent parent alleles in the selected progenies were only used for background selections in the succeeding backcross generations. From the BC3F3 generation onwards, the selected families were evaluated in field trials under stressed as well as unstressed conditions. Finally, at the BC3F5 stage, the selected NILs were reconfirmed for the presence of two target QTLs, qDTY2.1 and qDTY3.1. The agronomic evaluation was done for morphological and yield component traits as well as grain quality traits. RPG recovery in the NILs was depicted using Graphical GenoTypes (GGT) Version 2.0 software [37].

2.3. Evaluation of the NILs Under Stressed and Unstressed Conditions

A preliminary field evaluation of BC3F4 NILs was carried out along with both the parents, Pusa 44 and IR81896-B-B-142, together with two other checks, namely, IR81896-B-B-195 and IR87728-59-B-B. During Kharif 2016, agro-morphological evaluation was carried out using an augmented randomized block design (ARBD) with six blocks. Two treatments, stressed and unstressed, were maintained, with stressed plots having irrigation withheld during the entire heading period, beginning from 30 days after transplanting. The unstressed treatment had the normal recommended irrigation schedule. Barring irrigation, both the treatments had similar agronomic management. The NILs were grown in plots of 10.3 m2 following a spacing of 20 cm × 15 cm. The data were recorded on five uniform-looking healthy plants from each of the NILs for the characteristics, namely, days to 50% flowering (DtF), plant height (PtH), number of productive tillers (NpT), panicle length (PnL), grain yield (GrY), number of filled grains per panicle (FdG), spikelet fertility (SpF) and thousand grain weight (GrW). Plants were harvested at physiological maturity. The mean of five plants constitutes yield per plant for the genotype and plot yield was also measured for all the NILs.
Based on the agronomic and quality performance of the NILs under both stress and non-stress conditions, the set of 31 best-performing NILs in BC3F5 generation was evaluated in RBD along with parents and two checks in two replications during Kharif 2017 with two treatments, namely, moisture stressed and normal irrigated (unstressed). The plot size was 10.3 m2 with a spacing of 20 × 15 cm. Recommended agronomic management was followed in the trial, except for the irrigation schedule. Under unstressed conditions, irrigation was scheduled to maintain a 10 cm water level in the fields for the entire duration from transplanting to harvest. In the moisture-stressed plots, approximately 5 cm of standing water was maintained till panicle initiation, which coincided with 30 days after transplanting. After booting, irrigation was withheld totally until maturity. In the drought-stress plots, the moisture stress imposed was monitored using tensiometers and the irrigation was withheld till the stress level reached −70 KPa, to ensure stringent evaluation of the drought-tolerance in these NILs. Data recording was done for all the traits that were recorded earlier in the ARBD trial, from five random uniform-looking plants from each plot at maturity.
Both the trials during Kharif 2016 and Kharif 2017 were carried out at the research farm of Division of Genetics, ICAR-Indian Agricultural Research Institute (ICAR-IARI), New Delhi. In New Delhi, the average temperature during Kharif 2017 ranged between 38 °C in June to 15 °C in November, with a total rainfall of 856.0 mm.

2.4. Drought Tolerance Indices

Based on the performance of NILs under both stressed and unstressed conditions, the following drought indices were worked out for a comprehensive evaluation of the drought tolerance.
  • Drought susceptibility index, DSI   =   Y ns Y s Y ns [38]
  • Relative decrease yield, RDY   = 100 - ( Y i ) s ( Y i ) ns ×   100 [39]
  • Drought tolerance efficiency, DTE   ( % )   = Y s   Y ns   ×   100 [40]
where Ys and Yns represent the yields of all genotypes evaluated under drought stressed and unstressed conditions, respectively. (Yi)s denotes the yield of the ith genotype under drought stress and (Yi)ns the yield of the ith genotype under unstressed conditions.

2.5. Grain and Cooking Quality

The grain samples from the NILs were analyzed for the grain dimension and quality parameters as described by Ellur et al. [34]. The traits observed were hulling percentage (HgP), milling percentage (MgP), kernel length before cooking (KLBC), kernel width before cooking (KWBC) and length-width ratio (LWB), kernel length after cooking (KLAC), kernel width after cooking (KWAC), elongation ratio (ER) and alkali spreading value (ASV). All the measurements of physical parameters of the milled rice were done using a photo analyzer.

2.6. Statistical Analyses

The agro-morphological and grain quality data from the improved drought tolerant NILs of Pusa 44 were subjected to statistical analyses of variance and means [41] and analysis of variance (ANOVA). The mean comparison between the genotypes was carried out using the appropriate statistical assumption of normal distribution. The comparison was done using Tukey’s honestly significant test at the 5% confidence level. The conformability of genotype performance across both the seasons was computed using the Pearson correlation coefficient. All the data analyses were carried out using Analysis Toolpack plugin in Microsoft Excel, and STAR statistical analysis package [42]. The hierarchical agglomerative clustering of the multi-location data was performed under the R statistical environment using the “Stats” package [43].

3. Results

3.1. Introgression of the QTLs, qDTY2.1 and qDTY3.1

Screening of the genome-wide polymorphism between the recurrent parent, Pusa 44, and the donor parent, IR81896-B-B-142, using 460 SSR markers revealed 97 markers to be polymorphic, indicating 21.1% of genetic variability between the parents (Supplementary Figure S1). These polymorphic markers were further used for background selection. Among the polymorphic markers, two markers from chromosome 3, RM15791 (28.56 Mb) and RM16033 (32.56 Mb), were identified to be closely flanking qDTY3.1. Similarly, RM5791 (10.74 Mb) and RM324 (11.4 Mb) were identified flanking qDTY2.1 on chromosome 2. These markers were used for recombinant selection to help minimize linkage drag (Table 1).
The initial hybridization of Pusa 44 with IR81896-B-B-142 produced 11 seedlings, out of which 9 were found to be heterozygous for the QTL-linked markers during foreground selection (Table 2). The cross was designated as Pusa 1823. Two F1 plants were selected out of the nine true F1s and were backcrossed to the recurrent parent, Pusa 44, to generate forty-eight BC1F1 plants. Foreground selection using the RM520 identified 18 plants heterozygous for this marker. These plants on further genotyping using RM521 could identify eight plants heterozygous for both the QTLs. All the 8 BC1F1 plants were subjected to background selection with 97 polymorphic SSR markers, which indicated RPG recovery ranging between 73.08 and 82.87%. Further, the phenotypic selection of BC1F1 with Pusa 44 helped in identifying two BC1F1 plants with the RPG recovery of 82.0 and 82.7, respectively. These plants were also found to show the highest recurrent parent phenome (RPP) recovery. They were further backcrossed to Pusa 44 to generate 32 BC2F1 plants. Foreground selection using the marker RM520 among the BC2F1 plants revealed that nine plants were heterozygous for qDTY3.1. Additional foreground selection among these nine plants for the second QTL, qDTY2.1, identified three plants heterozygous for both the QTLs. Based on agro-morphological similarity to Pusa 44 and RPG recovery (98.5%), backcrossing was done using two BC2F1 plants. Among the 18 BC3F1 plants produced, only 1 BC3F1 plant was found to be heterozygous for both the QTLs, qDTY2.1 and qDTY3.1. The single BC3F1 plant on background selection revealed that it had an RPG recovery of 98.6%. Being a single progeny, this line was not subjected to morphological selection. The selfing of this BC3F1 plant produced 400 BC3F2 plants, of which 54 plants were identified as homozygous for both the markers linked to qDTY2.1 and qDTY3.1, and were selected from the BC3F2 population. From the selected plants, a large population of 9488 BC3F3 was raised, out of which a total of 242 single plant selections showing agronomic resemblance to Pusa 44 were initially made. Following the grain and cooking quality evaluations, the selections were subsequently reduced to 108 lines. In the BC3F4 generation, progenies from these selected lines were evaluated under both stress and non-stress conditions, reducing the selection to 31 BC3F5 families. Finally, the 31 NILs were further evaluated for drought stress response under replicated trials for two years.
Recombinant selection for the target QTLs was conducted to eliminate the undesirable effect of linkage drag, if any, that may adversely affect the agronomic recovery of the NILs. The recombinant selection conducted at BC3F4 generation revealed that the markers RM15791 and RM16033 on chromosome 3 have recovered Pusa 44 alleles at these loci, indicating a introgression of the QTL, qDTY3.1 free from linkage drag. Similarly, for qDTY2.1, the recovery of Pusa 44 alleles at the closest flanking markers, RM5791 and RM324, was complete (Figure 1).

3.2. Performance of NILs under Stressed vis-à-vis Unstressed Conditions in BC3F4

The preliminary screening for agronomic performance under stressed and unstressed conditions, with BC3F4 generation NILs, revealed significant variation for yield and related agronomic traits (Figure 2). The analysis of variance (ANOVA) for the ARBD is given in Table 3. The variation among the NILs was higher for most of the traits than that under the unstressed conditions. The ratio of variances indicated greater variation for traits such as hulling and milling recovery. Importantly, grain yield and plant height showed higher variance under unstressed situations than under stress. Similarly, the increased variance for grain quality traits was also noticed under unstressed conditions. The variation in the panicle weights of the NILs, under both the treatments, was relatively similar.
The average performances of NILs under stress and unstressed conditions (Table 4) indicated that NILs regained almost all agronomic parameters of the recurrent parent, Pusa 44 (Supplementary Table S1). The average days to fifty percent flowering among the NILs without stress was 110.5 days, with a range of 106 to 114 days, as against 108.5 days in Pusa 44. Under stress, however, there was a marginal increase of about 2–3 days for flowering in all the genotypes, including the recurrent parent. In the case of plant height, NILs varied between 92.7 cm and 109.9 cm under unstressed situations, as against the height of 97.8 cm for Pusa 44. Under drought stress, there was a decrease in height averaging about 4.0 cm. However, the decrease in tillering was remarkable under drought, to the tune of 15.4% over the tiller number under unstressed conditions. Similarly, other traits that showed significant reductions under stressed conditions were grain yield (32.9%) and weight of 1000 grains (16.7%). Grain yield per plant under unstressed treatment ranged between 15.2 g and 32.9 g among the NILs, with an average of 24.6 g as against 25.1 g in Pusa 44. Under RSDS, the grain yield in Pusa 44 was 40.6% of that under unstressed conditions, while among the NILs, 67.7% of the normal yield was realized under RSDS on an average. The yield under stress in NILs ranged between 7.7 g and 26.6 g, with a mean of 16.5 g. However, Pusa 44 yielded 10.2 g of grains per plant under RSDS.
Another yield component that was found largely affected under RSDS was grain weight. However, comparative grain weight reduction among the NILs and the recurrent parent was found to be almost similar, with a reduction of around 4.0 g per 1000 grains. Interestingly, grain dimension traits including raw and cooked kernel traits were found to show an inconspicuous difference between stressed and unstressed treatments, whereas milling recovery had an apparent reduction under drought. Milling percentage was found reduced by around 5% under RSDS in both NILs and Pusa 44, which however was considerably less in the donor parent, IR81896-B-B-142 (Supplementary Table S2).
The relative influence of RSDS on BC3F4 NILs as determined by three drought indices, namely, drought tolerance efficiency (DTE), relative decrease in yield (RDY) and drought susceptibility index (DSI), reflected that the 31 lines selected from 108 NILs were better performing over the recurrent parent, Pusa 44. The DTE of the NILs ranged between 62.74% (P1823–12-15) and 98.15% (P1823-12-83), which was significantly improved over Pusa 44, having a DTE of 40.60%. Similarly, RDY was in the range from 1.85% to 37.26% among the NILs, and DSI was ranging between 0.02% and 0.37%. These parameters again emphasized that the selection of 31 BC3F4 NILs for better adaptability and performance under RSDS was effective, whereas Pusa 44 indicated sensitive values for these parameters. The RDY of Pusa 44 was 59.4, while it showed DSI of 0.59. The NILs carrying the QTLs qDTY2.1 + qDTY3.1 therefore indicated that they can tolerate RSDS significantly as compared to the recipient parent (Supplementary Table S3), indicating that these lines can be further tested for their RSDS endurance. However, as far as the grain dimensions and cooking quality characteristics are concerned, the selected NILs were on par with the recipient parent Pusa 44 except for five NILs which possessed significantly lower kernel widths.

3.3. Yield under Drought and No Drought Conditions in Large-Scale Screening

The large-scale evaluation carried out during the 2017 Kharif season explicated the agronomic response of 31 NILs under both RSDS and unstressed treatments. Under unstressed conditions, the agronomic performances of the 31 NILs did not show any significant differences among themselves, or between the recurrent parent, Pusa 44, for traits such as plant height, number of productive tillers or panicle length. There was no significant variation for grain quality characteristics such as grain weight and kernel breadth before cooking as well. Similarly, under RSDS, the grain quality parameters showed non-significant deviation from Pusa 44. However, agronomic performance of NILs and Pusa 44 was on par for most of the traits under unstressed conditions, which under stress situation showed up a significant variation. This enabled the comparison of NILs for yield-related traits, unravelling the drought response patterns of different yield component traits. Moreover, the pattern of yield response of NILs both under stressed and unstressed conditions showed a good agreement between BC3F4 and BC3F5 generations (Table 5).
This table clearly indicates that the percent yield reduction under stress conditions was highest in the RP and was lower in the DP, as well as other checks. Moreover, the reduction in yield among the NILs (78.9%) was lower as compared to the RP (86.7%). The relative reduction under drought compared to unstressed among NILs was found to be 79.2% for yield per plant, and 88% for grain yield, which is very high, indicating the severity of the drought during Kharif 2017. The correlation coefficient for both the years of trial, Kharif 2016 and Kharif 2017, was found to be highly significant for all the characteristics, such as yield per plant, grain yield, days to 50% flowering and spikelet fertility, which reveals that the NILs were showing consistent performance during both the years of evaluation.

3.4. Performance under Stress Conditions

Under drought stress conditions, significant phenotypic variations were found for all morphological traits, particularly in some of the NILs (Supplementary Table S4). Besides this, no significant differences could be observed for grain and cooking quality traits for almost all the lines, denoting the recovery of quality parameters in the selected NILs. For days to 50% flowering, eleven NILs were found to be significantly earlier than Pusa 44, with earliness values of 1.5 to 2.5 days. In the case of plant height, most of the NILs had similar heights to Pusa 44, except a few that showed slightly taller plant stature. Most of the NILs possessed longer panicles than Pusa 44, whereas, only a few NILs showed panicle length similar to or lower than the RP. In general, the NILs were found to perform better than the recipient parent, Pusa 44, for other phenotypic characteristics, such as productive tiller number, plant biomass, spikelet fertility and panicle weight. A similar situation could be noticed with all the quality parameters. For grain yield per plant, all the NILs recorded significantly higher values than Pusa 44, except 4 NILs, and 17 NILs significantly out-yielded the recipient parent for plot yield (Table 6). Among the outperforming NILs, P1823-12-82 and P1823-12-122 were the top yielders under RSDS, which had shown 157% and 111% more yield than Pusa 44 (Figure 3). Besides this, there were 15 other NILs that showed more than a 50% increased yield over Pusa 44 under RSDS. The NILs that performed well under drought also showed a significantly high number of productive tillers and high spikelet fertility, indicating their improved adaptability. The grain quality parameters, however, did not show significant deviation from those of Pusa 44 under drought conditions (Supplementary Table S5). However, there was a significant reduction in the milling properties of grains when compared to unstressed conditions, specifically for traits such as hulling and milling percentages and head rice recovery.

3.5. Performance under Unstressed Conditions

Under unstressed conditions, plants were maintained under sufficient moisture with recommended irrigation schedules and did not experience any stress at any crop stage. The agronomic performance of the majority of the NILs showed no statistically significant deviation from Pusa 44 for all the traits (Supplementary Table S6). The on-par performance as compared to Pusa 44 indicated that the NILs have acquired all the agronomic characteristics of the recurrent parent. Concerning grain yield per plant, the plot yield of the NILs showed an advantage over Pusa 44 under unstressed conditions that was also, however, insignificant. However, there were two NILs that had lower plot yields than Pusa 44. A similar pattern was also found with single plant yield, as expected. The grain quality under unstressed conditions showed a significant increase in hulling and milling percentage and head rice recovery (Supplementary Table S7). Other grain-related parameters, such as the number of filled grains per panicle, grain weight and spikelet fertility, had significant deviations among some of the NILs, as compared to Pusa 44. In general, the plant height of the NILs was greater under unstressed conditions when compared to that under stressed situations. Similarly, there was a reduction in days to heading under drought. Twelve NILs showed significant deviation for days to 50% flowering, of which five showed a reduction in the days. In the case of filled grains per panicle, 23 out of 31 NILs released lesser grain number; however, spikelet fertility was better among the 17 NILs.

3.6. Comparative Performance in Terms of Drought Indices

A clear pattern of the drought endurance features of the NILs was revealed on comparative evaluation using the drought-related indices (Table 6). Among the 31 NILs that were selected during the previous season, 28 NILs showed better DTE values when subjected to a large-scale replicated trial in Kharif 2017. The DTE of Pusa 44 was only 13.3%, while the tolerant donor (IR81896-B-B-142) and the tolerant check (IR81896-B-B-195) had DTE values of 34.4% and 43.9%, respectively. Among the NILs, the DTE ranged between 10.1% and 30.4%, with an average of 20.7%. The highest DTE among the NILs was recorded in P1823-12-82 followed by P1823-12-10. Further, there were 15 NILs that showed above-average DTE values. The other drought impact measure, the relative decrease in yield (RDY), showed a decrease of 86.7% in Pusa 44, while IR 64 had a reduction of 78.4%. The tolerance checks had RDY values ranging between 56.1 and 65.7%. Among the NILs, however, the RDY values ranged between 69.7% and 89.9%, with an average RDY of 79.3%. When compared to the previous season, the reduction in yield was higher in Kharif 2017, which was due to the higher severity of drought in the later season. The third measure, DSI, that ranges between zero (0) and one (1), indicated a continuous measure of drought response from sensitivity to tolerance. However, the DSI showed the same pattern as RDY. Relatively, the NILs had an average advantage of more than 100% over Pusa 44. Although there were few NILs that showed drought responses similar to or lower than that of Pusa 44, the proportion of such NILs was relatively very low among the NILs. Based on the agronomic performance and grain quality parameters, as well as considering the drought response indices, two NILs, Pusa 1823-12-62 and Pusa 1823-12-82, were nominated under drought screening trials under the All India Coordinated Rice Improvement Project (AICRIP) for national testing. These NILs showed DTE values of 21.1% and 35.4%, respectively, with plot yields of 5279.0 kg ha−1 and 5454.8 kg ha−1 under unstressed conditions. Additionally, these NILs also showed better or similar grain quality compared to Pusa 44, with an added advantage of acceptable head rice recovery.

3.7. Multilocation Evaluation for Varietal Identification

Two of the top performing NILs that were evaluated in the national testing of drought tolerant rice cultivars performed better than the recurrent parent, Pusa 44, at two locations when grown in the two rainfed locations which are prone to drought stress for two years (Figure 4). During Kharif 2019, the average performance of the NILs at both Hazaribagh and Rewa was significantly superior to the recurrent parent, Pusa 44. However, under a rainout shelter, screening during 2018 at Hazaribagh indicated a slightly different pattern of RSDS response among the NILs. This notwithstanding, during both the seasons and at both the locations, the NILs yields were significantly better than Pusa 44. Among the locations, the performance of the NILs was better than those of both the parents, Pusa 44 and IR81896-B-B-142, except for Hazaribagh in 2018. Although the NILs produced significantly higher yields as compared to the recurrent Pusa 44 at Hazaribagh during 2018, there was a marginally but insignificantly lower yield than the donor parent, IR81896-B-B-142.

4. Discussion

In the wake of climate change, drought has undoubtedly become the most severe production constraint of all food crops, especially rice. In recent years, there have been frequent episodes of water scarcity in rice production worldwide, that have resulted in adverse effects. The consistent occurrence of drought challenges the sustenance of rice cultivation, particularly of grain yield. Drought is often dubious and occurs either with a sudden shift in climatic pattern, especially in rainfall pattern, or with a prolonged interval without rains. Coupled with high temperatures, the adverse effects of drought become multifold, often leading to crop failures. Without adequate resources, which is a major constraint in the economically weaker and agriculture-dependent nations, managing drought in physical terms is an almost impractical proposition. Harnessing external sources of water supply, which requires huge infrastructural improvements, becomes seldom practicable under erratic weather occurrences. Therefore, a more sustainable solution is to alter the cropping pattern to coincide with the probable drought occurrence, or to grow varieties that have more drought endurance features.
Drought tolerance in crops can be of two types—escape or endurance. Escape is often managed with the use of short duration varieties of crops that can complete their life cycle before the onset of drought. However, a lack of adequate short duration varieties is a major challenge in annual crops such as rice, and is not an option for perennial species. A prolonged period without precipitation may also affect the escape mechanism, leading to failure. Therefore, endurance can be considered as a better alternate for drought management in a wider perspective. Drought endurance can assure some grain production, where nothing is expected, by growing a drought-tolerant cultivar. Fortunately, there are several mechanisms of drought endurance in crops that are genetically determined. This offers us the opportunity of transferring the corresponding genes into a sensitive variety. With the characterization of these genes, or the genomic locations they reside in, it is now possible to integrate them into a new varietal background using DNA-based detectable markers. Marker-assisted breeding, particularly using the backcross method, has been extensively used for introgressing genes governing resistance/tolerance biotic as well as abiotic stresses in rice [17,44].
In rice, drought can affect crop growth in various stages, ultimately resulting in different outcomes. In the seedling stage, drought hinders population establishment, while that at the vegetative stage affects crop proliferation and tillering. Although drought at early growth stages may indirectly affect the crop yield, drought stress at the reproductive stage is more crucial because it directly affects grain development and ultimately the yield. Therefore, to breed rice varieties with drought endurance, RSDS tolerance is the most important trait to be harnessed. There are many QTLs governing tolerance to RSDS that have been mapped in rice [15,16,17,24], and many of them have been transferred into the well-known varieties [2]. Among these, qDTY3.1 and qDTY2.1 are two large-effect QTLs governing tolerance under severe lowland drought stress [16]. Identified on chromosome 3 from Apo, qDTY3.1 was found to explain 31% of genetic variance for grain yield, while qDTY2.1 on chromosome 2 had a significant effect on grain yield (13–16%). In this study, we have targeted these two QTLs for introgression into the popular rice variety, Pusa 44, using a NIL developed in the genetic background of Swarna with both qDTY2.1 and qDTY3.1 (IR81896-B-B-142) as the donor parents [16]. However, QTLs, such as qDTY3.1, have been reported as linked to undesirable traits, such as yield decline under non-stress [16,45,46]. This necessitates the breaking of such linkages for the effective transfer of the trait by marker-assisted backcross breeding (MABB) program [25]. The use of combined selection for the marker as well as the phenotype on a large breeding population can ensure the breaking of such linkages.

4.1. The severity of Drought Stress

During the two years of field evaluation of the Pusa 44 NILs, the severity of the drought was much higher during Kharif 2017 than during the previous season of Kharif 2016. During 2017, the RDY was found to be more than 69.7% for all the NILs, while during Kharif 2016 the maximum RDY was found to be 35.8%, indicating an almost twofold reduction during 2017. During 2016, however, the NILs showed better yield, suggesting a lesser drought impact during this year. Although less, the relative mean yield reduction was apparent in 2016 also, under the stress conditions; however, it was more prominent during 2017, indicating the more severe nature of drought during this season. Indicating this, during this experiment, the soil developed deep cracks due to insufficient soil moisture. In another study, Yambao et al. [47] reported a reduction of up to 70% in yield upon imposing drought for 15 days at the panicle initiation stage, and 88 and 52% reductions during stress imposition in flowering and grain filling stage, respectively. Kumar et al. [48] recommended that a screening method that could lessen the mean yield of the lines under study by at least 65% under severe stress as compared to irrigated non-stress is necessary to identify true drought-tolerant lines.

4.2. Agronomic Performance of Pusa 44 NILs under Stress and Non-Stress Conditions

Days to fifty percent flowering is a characteristic which is extremely sensitive to drought. The shift in flowering duration under stress conditions in the present experiment indicates that water stress at the initiation of the reproductive stage affected flowering time. Similar results were reported by earlier studies [15,23,45,49,50,51,52,53,54,55,56]. Bernier et al. [9] and Venuprasad et al. [16] reported that qDTY12.1 and qDTY3.1 influenced both grain yield and days to flowering under stress conditions, evincing that genes present within the QTLs are likely to be allied with early flowering. In the present experiment, stress was imposed 30 days after transplanting at the initiation of reproductive stage to ensure that the lines with the earliest flowering cannot escape drought, and only genuine drought-tolerant lines are selected. Similarly, there was also a significant suppression of plant height under severe stress, particularly during the 2017 season. The growth suppression further indicated that the plants did begin to experience stress much before entering into the reproductive phase. This could help us in divulging the relative endurance of the plants under prolonged drought, beginning from floral primordial initiation. We therefore presume that the drought tolerance exhibited by the NILs in the present study would be more pragmatic and would not be accounted for due to drought escape, as escape can occur solely under a transient drought rather than a prolonged one.

4.3. Relative Drought Tolerance of the Improved NILs

In the present study, all the selected lines showed a significant yield advantage over the recurrent parent under severe drought conditions. They also possessed grain and cooking quality on par with Pusa 44 during both years of evaluation. The performance of the NILs under unstressed situations proved that they recovered the maximum agronomic and grain quality attributes of Pusa 44 during the selection process. As successfully reported earlier [50,51,52,53], a stringent phenotypic selection along with marker-based background selection ensured an RPG recovery ranging between 98.6 and 99.4% after three successive backcrosses. This was further apparent from the recombinant selection data obtained later in the selection cycle, during the BC3F4 stage, that indicated complete recovery of carrier chromosomes.
Another trait that showed significant improvement among the Pusa 44 NILs was spikelet fertility. There were eight NILs that showed better, higher spikelet fertility than Pusa 44. A higher spikelet fertility in combination with a greater number of spikelets per panicle can ultimately result in better yield under RSDS, as observed in some of the NILs in this study. Although the relative drought response of the NILs was tested using three independently derived parameters, in practice, only two parameters were employed. This is because two of the indices, DSI and RDY, were the same, having been expressed by similar formulae. The DTE of the NILs showed significantly higher values as compared to Pusa 44, indicating the advantage they acquired by integrating both the QTLs, qDT2.1 and qDTY3.1. Further, we could achieve an effective transfer of two QTLs together in three years of successive backcrossing along with phenotypic selection for agromorphological as well as grain quality traits. Additionally, selection among a large BC3F3 population could aid us in eliminating the undesirable effects associated with QTL, such as qDTY3.1, which was reported to be associated with yield decline under normal conditions. Among the NILs tested, except for eight NILs, all the remaining lines showed improved yield compared to Pusa 44, most of which was statistically insignificant. Compared by drought indices, 14 NILs out of 31 showed significantly superior performance as compared to the Pusa 44 parent, among which P1823-12-1, P1823-12-14, P1823-12-30, P1823-12-55, P1823-12-62, P1823-12-77, P1823-12-82, P1823-12-122, P1823-12-133, P1823-12-137 and P1823-12-140 showed the least RDI and DSI during both years of evaluation. Furthermore, P1823-12-14, P1823-12-77, P1823-12-82, P1823-12-122, P1823-12-137 and P1823-12-140 also had significantly higher spikelet fertility under stress conditions, and significantly lower days to flowering than the recipient parent. The NIL P1823-12-82 produced the highest yield per plant, as well as plot yield, in the replicated trial under the stress conditions. Additionally, P1823-12-14, P1823-12-55 and P1823-12-62 also produced comparatively good yields under both drought-stressed as well as unstressed conditions, along with equivalent grain and cooking as compared to Pusa 44.
We could demonstrate that QTL introgression using MABB is the most efficient way to transfer RSDS tolerance into mega/popular varieties, such as Pusa 44. Since they have acquired most of the recurrent parent features together with the target trait, these lines can be directly tested under the target environments to which the recurrent parent was originally released or widely cultivated. Moreover, the drought-tolerant NILs developed in the present study could be effectively used as potential donors for future breeding, and can be evaluated in multi-location trials to be released for commercial cultivation.

5. Conclusions

The marker-assisted transfer of two major QTLs, qDTY2.1 and qDTY3.1, for grain yield under reproductive stage drought stress, into the variety Pusa 44 was accomplished, which will help in managing yield losses under drought stress, which is a major stress affecting rice production, especially under the changing climate and water limitations imposed by this change. These NILs are potential candidates for release as a variety after the required testing for replacing Pusa 44, which is highly susceptible to drought and still occupies a substantial area under cultivation. Additionally, these NILs can be used as donor lines for drought tolerance and as a component of a multiline in breeding programs. The deployment of drought-tolerant Pusa 44 NILs could reduce the water requirement and could practically cope with the present scenario of drastic climatic changes without compromising the yield and grain quality of the variety. The improved lines will also serve as an invaluable source of drought tolerance in rice improvement.

Supplementary Materials

Supplementary Materials can be found at https://0-www-mdpi-com.brum.beds.ac.uk/2077-0472/11/1/64/s1. Supplementary Figure S1: Marker-assisted backcross breeding scheme adopted for introgression of qDTY3.1 and qDTY2.1 in Pusa 44; Supplementary Figure S2: Graphical representation of all polymorphic markers used for background recovery; Supplementary Table S1: Yield performance of Pusa 44 qDTY NILs under drought-stressed and unstressed conditions in the BC3F4 generation (Kharif 2016); Supplementary Table S2: Grain dimensions and quality characteristics of Pusa 44 qDTY NILs under BC3F4 generation; Supplementary Table S3: Yield differential under stressed and unstressed conditions and drought indices of Pusa 44 qDTY NILs; Supplementary Table S4: Performance of selected NILs under drought stress conditions in Kharif 2017 season (BC3F5 generation); Supplementary Table S5: Quality data of Pusa 44 qDTY NILs under drought stress during Kharif 2017; Supplementary Table S6: Performance of Pusa 44 qDTY NILs under unstressed conditions in Kharif 2017; Supplementary Table S7: Quality parameters of selected NILs under unstressed conditions in Kharif 2017.

Author Contributions

Conceptualization, S.G.K.; methodology, S.G.K., A.K.S., M.P.S.; software, K.K.V., P.K.B., R.K.E., S.G.K.; validation, P.D., G.D.; formal analysis, P.D., G.D., K.K.V., H.B.; investigation, P.D., S.G.K., N.P.M., P.P.; resources, S.G.K., A.K.S., M.P.S., M.N.; data curation, P.D., K.K.V.,S.G.K.; writing—original draft preparation, P.D., K.K.V., S.G.K.; writing—review and editing, K.K.V., S.G.K.; visualization, K.K.V., P.D., S.G.K.; supervision, S.G.K., A.K.S., N.R.; project administration, S.G.K., M.P.S., A.K.S.; funding acquisition, S.G.K., M.P.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the ICAR-funded network project on National Innovations on Climate Resilient Agriculture (NICRA) and National Agricultural Higher Education Project (NAHEP)—Centre for Advanced Agricultural Science and Technology (CAAST) project on “Genomics Assisted Breeding for Crop Improvement”, of the World Bank and the Indian Council of Agricultural Research.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All the data used are provided in the supplementary files.

Acknowledgments

The study is part of the Ph.D. thesis research of the first author (P.D.) and the first author wishes to acknowledge ICAR-Indian Agricultural Research Institute, New Delhi and Indian Council of Agriculture Research (ICAR) for providing all the necessary facilities.

Conflicts 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. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Abbreviations

ANOVAAnalysis of variance
ARBDAugmented randomized block design
ASVAlkali spreading value
BmPBiomass per plant
DPDonor parent
DSIDrought susceptibility index
DTEDrought tolerance efficiency
DtFDays to 50% flowering
EgRElongation ratio
FdGNumber of whole grains per panicle
GGTGraphical Genotypes
GrWWeight of 1000 grains
GrYGrain yield
HgPHulling percent
IARIIndian Agricultural Research Institute
ICARIndian Council of Agricultural Research
KLACKernel length after cooking
KLBCKernel length before cooking
KWACKernel width after cooking
KWBCKernel width before cooking
LWRLength/width ratio
MABBMarker-assisted backcross breeding
MASMarker-assisted selection
MgPMilling percentage
NILNear isogenic lines
NpTNumber of panicle-bearing tillers
PDSPublic distribution system
PnLLength of panicle
PnWPanicle weight
PtHPlant height
QTLQuantitative trait loci
RBDRandomized block design
RBGRCRice Breeding and Genetics Research Centre
RDYRelative decrease in yield
RPRecurrent parent
RPGRecurrent parent genome
RPPRecurrent parent phenome
RSDSReproductive stage drought stress
SStressed
SpFSpikelet fertility
SSRSimple sequence repeat
USUnstressed
YdPYield per plant

References

  1. Carriger, S.; Vallee, D. More crop per drop. Rice Today 2007, 6, 10–13. [Google Scholar]
  2. Singh, R.; Singh, Y.; Xalaxo, S.; Verulkar, S.; Yadav, N.; Singh, S.; Singh, N.; Prasad, K.S.N.; Kondayya, K.; Rao, P.V.R.; et al. From QTL to variety-harnessing the benefits of QTLs for drought, flood and salt tolerance in mega rice varieties of India through a multi-institutional network. Plant Sci. 2016, 242, 278–287. [Google Scholar] [CrossRef] [PubMed]
  3. Joshi, K.; Joshi, P.K.; Khan, M.T.; Kishore, A. ‘Sticky Rice’: Variety inertia and groundwater crisis in a technologically progressive state of India. In IFPRI Discuss; Paper 01766; International Food Policy Research Institute: Washington, DC, USA, 2018; Available online: http://ebrary.ifpri.org/cdm/singleitem/collection/p15738coll2/id/132924 (accessed on 13 January 2021).
  4. Nairizi, S. Drought and Water Scarcity Drought and Water Scarcity; International Commission on Irrigation and Drainage: New Delhi, India, 2017; Available online: https://www.icid.org/drought_pub2017.pdf (accessed on 13 January 2021).
  5. Lamaoui, M.; Jemo, M.; Datla, M.; Bekaoui, F. Heat and Drought Stresses in Crops and Approaches for Their Mitigation. Front. Chem. 2018, 6. [Google Scholar] [CrossRef] [PubMed]
  6. Courtois, B.; McLaren, G.; Sinha, P.K.; Prasad, K.; Yadav, R.; Shen, L. Mapping QTL associated with drought avoidance in upland rice. Mol. Breed. 2000, 6, 55–66. [Google Scholar] [CrossRef]
  7. Lafitte, H.R.; Li, Z.K.; Vijayakumar, C.H.M.; Gao, Y.M.; Shi, Y.; Xu, J.L. Improvements of rice drought tolerance through backcross breeding: Evaluation of donors and selection in drought nurseries. Field Crops Res. 2006, 97, 77–86. [Google Scholar] [CrossRef]
  8. Venuprasad, R.; Lafitte, H.R.; Atlin, G.N. Response to direct selection for grain yield under drought stress in rice. Crop Sci. 2009, 47, 285–293. [Google Scholar] [CrossRef]
  9. Bernier, J.; Kumar, A.; Venuprasad, R.; Spaner, D.; Verulkar, S.; Mandal, N.P.; Sinha, P.K.; Peeraju, P.; Dongre, P.R.; Mahto, R.N.; et al. Characterization of the effect of rice drought tolerance QTL12.1 over a range of environments in the Philippines and eastern India. Euphytica 2008, 166, 207–217. [Google Scholar] [CrossRef]
  10. Swamy, M.; Ahmed, H.U.; Henry, A.; Mauleon, R.; Dixit, S.; Vikram, P.; Tilatto, R.; Verulkar, S.B.; Perraju, P.; Mandal, N.P.; et al. Genetic, Physiological, and Gene Expression Analyses Reveal That Multiple QTL Enhance Yield of Rice Mega-Variety IR64 under Drought. PLoS ONE 2013, 8. [Google Scholar] [CrossRef]
  11. Sandhu, N.; Kumar, A. Bridging the Rice Yield Gaps under Drought: QTLs, Genes, and Their Use in Breeding Programs. Agronomy 2017, 7, 27. [Google Scholar] [CrossRef]
  12. Sandhu, N.; Singh, A.; Dixit, S.; Cruz, M.T.S.; Maturan, P.C.; Jain, R.K.; Kumar, A. Identification and mapping of stable QTL with main and epistasis effect on rice grain yield under upland drought stress. BMC Genet. 2014, 15, 63. [Google Scholar] [CrossRef] [Green Version]
  13. Palanog, A.; Swamy, B.P.M.; Shamsudin, N.A.A.; Dixit, S.; Hernandez, J.E.; Boromeoa, T.H.; Cruza, P.S.; Kumar, A. Grain yield QTLs with consistent effect under reproductive-stage drought stress in rice. Field Crops Res. 2014, 161, 46–54. [Google Scholar] [CrossRef]
  14. Mishra, K.K.; Vikram, P.; Yadaw, R.B.; Mallikarjuna Swamy, B.G.; Dixit, S.; Cruz, M.T.S.; Maturan, P.; Marker, S.; Kumar, A. qDTY12.1: A locus with a consistent effect on grain yield under drought in rice. BMC Genet. 2013, 14, 12. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Ghimire, K.H.; Lenie, A.Q.; Vikram, P.; Swamy, B.P.M.; Dixit, S.; Ahmed, H. Identification and mapping of a QTL (qDTY1.1) with a consistent effect on grain yield under drought. Field Crops Res. 2012, 131, 88–96. [Google Scholar] [CrossRef]
  16. Venuprasad, R.; Dalid, C.O.; Del Valle, M.; Zhao, D.; Espiritu, M.; Cruz, M.T.S.; Amante, M.; Kumar, A.; Atlin, G.N. Identification and characterization of large-effect quantitative trait loci for grain yield under lowland drought stress in rice using bulk-segregant analysis. Theor. Appl. Genet. 2009, 120, 177–190. [Google Scholar] [CrossRef] [PubMed]
  17. Muthu, V.; Abbai, R.; Nallathambi, J.; Rahman, H.; Ramasamy, S.; Kambale, R.; Thulasinathan, T.; Ayyenar, B.; Muthurajan, R. Pyramiding QTLs controlling tolerance against drought, salinity, and submergence in rice through marker assisted breeding. PLoS ONE 2020, 15, e0227421. [Google Scholar] [CrossRef] [Green Version]
  18. Shamsudin, N.A.A.; Mallikarjuna Swamy, B.P.; Ratnam, W.; Cruz, M.T.S.; Sandhu, N.; Raman, A.K.; Kumar, A. Marker assisted pyramiding of drought yield QTLs into a popular Malaysian rice cultivar, MR219. BMC Genet. 2016, 9, 21. [Google Scholar] [CrossRef] [Green Version]
  19. Sandhu, N.; Dixit, S.; Swamy, B.P.M.; Vikram, P.; Venkateswaralu, C.; Catolos, M.; Kumar, A. Positive interactions of major-effect QTLs with genetic background that enhances rice yield under drought. Sci. Rep. 2018, 8, 1626. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Yadaw, R.B.; Dixit, S.; Raman, A.; Mishra, K.K.; Vikram, P.; Swamy, B.P.M.; Cruz, M.T.S.; Maturan, P.T.; Pandey, M.; Kumar, A. A QTL for high grain yield under lowland drought in the background of popular rice variety Sabitri from Nepal. Field Crops Res. 2013, 144, 281–287. [Google Scholar] [CrossRef]
  21. Kumar, A.; Dixit, S.; Ram, T.; Yadaw, R.B.; Mishra, K.K.; Mandal, N.P. Breeding high-yielding drought-tolerant rice: Genetic variations and conventional and molecular approaches. J. Exp. Bot. 2014, 65, 6265–6278. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Price, A.H.; Courtois, B. Mapping QTLs associated with drought resistance in rice progress, problems and prospects. Plant Growth Regul. 1999, 29, 123–133. [Google Scholar] [CrossRef]
  23. Jongdee, B.; Fukai, S.; Cooper, M. Leaf water potential and osmotic adjustment as physiological traits to improve drought tolerance in rice. Field Crops Res. 2002, 76, 153–163. [Google Scholar] [CrossRef]
  24. Bernier, J.; Kumar, A.; Venuprasad, R.; Spaner, D.; Atlin, G.N. A large-effect QTL for grain yield under reproductive-stage drought stress in upland rice. Crop Sci. 2007, 47, 507–516. [Google Scholar] [CrossRef]
  25. Dixit, S.; Yadaw, R.B.; Mishra, K.K.; Kumar, A. Marker-assisted breeding to develop the drought-tolerant version of Sabitri, a popular variety from Nepal. Euphytica 2017, 213, 184. [Google Scholar] [CrossRef] [PubMed]
  26. Singh, S.; Sidhu, J.S.; Huang, N.; Vikal, Y.; Li, Z.K.; Brar, D.S.; Dhaliwal, H.S.; Khush, G.S. Pyramiding three bacterial blight resistance genes (xa5, xa13 and Xa21) using marker-assisted selection into indica rice cultivar PR-106. Theor. Appl. Genet. 2001, 102, 1011–1015. [Google Scholar] [CrossRef]
  27. Gopalakrishnan, S.; Sharma, R.K.; Anand, R.K.; Joseph, M.; Singh, V.P.; Singh, A.K.; Bhat, K.V.; Singh, N.K.; Mohapatra, T. Integrating marker assisted background analysis with foreground selection for identification of superior bacterial blight resistant recombinants in Basmati rice. Plant Breed. 2008, 127, 131–139. [Google Scholar] [CrossRef]
  28. Basavaraj, S.H.; Singh, V.K.; Singh, A.; Singh, D.; Nagarajan, M.; Mohapatra, T. Marker aided improvement of Pusa6B, the maintainer parent of hybrid Pusa RH10, for resistance to bacterial blight. Indian J. Genet. 2009, 69, 10–16. [Google Scholar]
  29. Basavaraj, S.H.; Singh, V.K.; Singh, A.; Singh, A.; Singh, A.; Yadav, S. Marker-assisted improvement of bacterial blight resistance in parental lines of Pusa RH10, a superfine grain aromatic rice hybrid. Mol. Breed. 2010, 26, 293–305. [Google Scholar] [CrossRef]
  30. Singh, A.K.; Gopalakrishnan, S.; Singh, V.P.; Prabhu, K.V.; Mohapatra, T.; Singh, N.K. Marker assisted selection: A paradigm shift in Basmati breeding. Indian J. Genet. 2011, 71, 120–128. [Google Scholar]
  31. Singh, A.; Singh, V.K.; Singh, S.P.; Pandian, R.T.P.; Ellur, R.K.; Singh, D. Molecular breeding for the development of multiple disease resistance in Basmati rice. AoB Plants 2012. [Google Scholar] [CrossRef]
  32. Singh, V.K.; Singh, A.; Singh, S.P.; Ellur, R.K.; Singh, D.; Gopalakrishnan, S.G. Marker assisted simultaneous but stepwise backcross breeding for pyramiding blast resistance genes Piz5 and Pi54 into an elite Basmati rice restorer line “PRR78”. Plant Breed. 2013, 132, 486–495. [Google Scholar] [CrossRef]
  33. Khanna, A.; Sharma, V.; Ellur, R.K.; Shikari, A.B.; Gopalakrishnan, S.; Singh, U.D. Development and evaluation of near isogenic lines for major blast resistance gene(s) in Basmati rice. Theor. Appl. Genet. 2015, 128, 1243–1259. [Google Scholar] [CrossRef] [PubMed]
  34. Ellur, R.K.; Khanna, A.; Gopalakrishnan, S.; Bhowmick, P.K.; Vinod, K.K.; Nagarajan, M.; Mondal, K.K.; Singh, N.K.; Singh, K.; Prabhu, K.V.; et al. Marker-aided Incorporation of Xa38, a Novel Bacterial Blight Resistance Gene, in PB1121 and Comparison of its Resistance Spectrum with xa13 + Xa21. Sci. Rep. 2016, 6, 29188. [Google Scholar] [CrossRef] [PubMed]
  35. Babu, N.N.; Gopalakrishnan, S.; Vinod, K.K.; Krishnamurthy, S.L.; Singh, V.K.; Singh, M.P.; Singh, R.; Ellur, R.K.; Rai, V.; Bollinedi, H.; et al. Marker Aided Incorporation of Saltol, a Major QTL Associated with Seedling Stage Salt Tolerance, into Oryza sativa ‘Pusa Basmati 1121’. Front. Plant Sci. 2017, 8. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Murray, H.G.; Thompson, W.F. Rapid isolation of high molecular weight DNA. Nucleic Acids Res. 1980, 8, 4321–4325. [Google Scholar] [CrossRef] [Green Version]
  37. VanBerloo, R. GGT: Software for display of graphical genotypes. J. Hered. 1999, 90, 328–329. [Google Scholar] [CrossRef] [Green Version]
  38. Fischer, R.A.; Maurer, R. Drought resistance in spring wheat cultivars. 1. Grain yield responses. Aust. J. Agric. Res. 1978, 29, 897–912. [Google Scholar] [CrossRef]
  39. İlker, E.; Tatar, Ö.; Aykut-Tonk, F.; Tosun, M.; Turk, J. Determination of Tolerance Level of Some Wheat Genotypes to Post-Anthesis Drought. Turk. J. Field Crops 2011, 6, 59–63. [Google Scholar]
  40. Fischer, K.S.; Wood, G. Breeding and selection for drought tolerance in tropical maize. In Proceedings of the Symposium on Principles and Methods in Crop Improvement for Drought Resistance with Emphasis on Rice; International Rice Research Institute: Los Baños, Philippines, 1981. [Google Scholar]
  41. Gomez, K.A.; Gomez, A.A. Statistical Procedures for Agricultural Research; John Wiley and Sons: New York, NY, USA, 1984. [Google Scholar]
  42. IRRI. Statistical Tools for Agricultural Research, STAR, version 2.0.1; Biometrics and Breeding Informatics; International Rice Research Institute: Los Baños, Philippines, 2014. [Google Scholar]
  43. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2020; Available online: www.R-project.org (accessed on 13 January 2021).
  44. Septiningsih, E.M.; Pamplona, A.M.; Sanchez, D.L.; Neeraja, C.N.; Vergara, G.V.; Heuer, S.; Ismail, A.M.; Mackill, D.J. Development of submergence-tolerant rice cultivars: The Sub1 locus and beyond. Ann. Bot. 2009, 103, 151–160. [Google Scholar] [CrossRef] [Green Version]
  45. Dixit, S.; Swamy, B.P.M.; Vikram, P.; Ahmed, H.U.; Cruz, S.T.A.; Amante, M.; Atri, D.; Leiung, H.; Kumar, A. Fine mapping of QTLs for rice grain yield under drought reveals sub-QTLs conferring a response to variable drought intensities. Theor. Appl. Genet. 2012, 125, 155–169. [Google Scholar] [CrossRef]
  46. Vikram, P.; Swamy, B.P.M.; Dixit, S.; Trinidad, J.; Cruz, S.T.A.; Maturan, P.C.; Amante, M.; Kumar, A. Linkages and interaction analysis of major effect drought grain yield QTLs in rice. PLoS ONE 2016, 11, e0151532. [Google Scholar] [CrossRef] [Green Version]
  47. Yambao, E.B.; Ingram, K.T. Drought stress index for rice. Philipp. J. Crop Sci. 1988, 13, 105–111. [Google Scholar]
  48. Kumar, A.; Bernier, J.; Verulkar, S.; Lafitte, H.R.; Atlin, G.N. Breeding for drought tolerance: Direct selection for yield, response to selection and use of drought-tolerant donors in upland and lowland-adapted populations. Field Crops Res. 2008, 7, 221–231. [Google Scholar] [CrossRef]
  49. Vikram, P.; Swamy, B.P.M.; Dixit, S.; Ahmad, H.; Cruz, M.T.S.; Singh, A.K. qDTY1.1, a major QTL for rice grain yield under drought with a consistent effect in multiple elite genetic backgrounds. BMC Genet. 2011, 12, 1–30. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Ellur, R.K.; Khanna, A.; Yadav, A.; Pathania, S.; Rajashekara, H.; Singh, V.K. Improvement of Basmati rice varieties for resistance to blast and bacterial blight diseases using marker assisted backcross breeding. Plant Sci. 2016, 242, 330–341. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Sagar, V.; Gopala Krishnan, S.; Dwivedi, P.; Mondal, K.K.; Prakash, G.; Nagarajan, M.; Singh, A.K. Development of Basmati rice genotypes with resistance to both bacterial blight and blast diseases using marker assisted restricted backcross breeding. Indian J. Genet 2018, 78, 36–47. [Google Scholar] [CrossRef]
  52. Sagar, V.; Dhawan, G.; Krishnan, S.G.; Vinod, K.K.; Ellur, R.K.; Mondal, K.K.; Rathour, R.; Prakash, G.; Nagarajan, M.; Bhowmick, P.K.; et al. Marker assisted introgression of genes governing resistance to bacterial blight and blast diseases into an elite Basmati rice Variety, ‘Pusa Basmati 1509’. Euphytica 2020, 216, 16. [Google Scholar] [CrossRef]
  53. Vinod, K.K.; Krishnan, S.G.; Thribhuvan, R.; Singh, A.K. Genetics of Drought Tolerance, Candidate Genes and Their Utilization in Rice Improvement. In Genomics Assisted Breeding of Crops for Abiotic Stress Tolerance, Vol. II, Sustainable Development and Biodiversity 21; Rajpal, V.R., Ed.; Springer Nature Switzerland AG: Cham, Switzerland, 2019. [Google Scholar] [CrossRef]
  54. Ouk, M.; Basnayake, J.; Tsubo, M.; Fukai, S.; Fischer, K.S.; Cooper, M. Use of drought response index for identification of drought tolerant genotypes in rain fed lowland rice. Field Crops Res. 2006, 99, 48–58. [Google Scholar] [CrossRef]
  55. Zhao, X.; Qin, Y.; Sohn, J.K. Identification of main effects, epistatic effects and their environmental interactions of QTLs for yield traits in rice. Genes Genom. 2010, 32, 37–45. [Google Scholar] [CrossRef]
  56. Blum, A. Plant Water Relations, Plant Stress and Plant Production. In Plant Breeding for Water-Limited Environments; Blum, A., Ed.; Springer: New York, NY, USA, 2011; pp. 11–52. [Google Scholar] [CrossRef]
Figure 1. Graphical representation of recurrent parent genome recovery and quantitative trait locus (QTL) introgression among the near-isogenic lines of Pusa 44 in BC3F4 generation. Chromosome segments in red represent Pusa 44 type, while blue segments represent that from the donor, IR81896-B-B-142. Heterozygous segments are shown in a dull grey color. The complete recovery of recurrent parent genome underpins the importance of augmenting phenotypic selection with the marker-assisted background selection.
Figure 1. Graphical representation of recurrent parent genome recovery and quantitative trait locus (QTL) introgression among the near-isogenic lines of Pusa 44 in BC3F4 generation. Chromosome segments in red represent Pusa 44 type, while blue segments represent that from the donor, IR81896-B-B-142. Heterozygous segments are shown in a dull grey color. The complete recovery of recurrent parent genome underpins the importance of augmenting phenotypic selection with the marker-assisted background selection.
Agriculture 11 00064 g001
Figure 2. Field view of the selected Pusa 44 qDTY NILs in BC3F4 generation under irrigated conditions. NILs, in general, showed agronomic superiority over the recurrent parent, Pusa 44 under drought compared to irrigated conditions. The grain quality of NILs was simialr to Pusa 44, in most of the NILs.
Figure 2. Field view of the selected Pusa 44 qDTY NILs in BC3F4 generation under irrigated conditions. NILs, in general, showed agronomic superiority over the recurrent parent, Pusa 44 under drought compared to irrigated conditions. The grain quality of NILs was simialr to Pusa 44, in most of the NILs.
Agriculture 11 00064 g002
Figure 3. Comparative field view of Pusa 44 near isogenic line P1823-12-82 carrying two QTLs for reproductive stage drought stress tolerance under drought-stressed conditions. Note that P1823-12-82 had better grain filling and spikelet fertility than Pusa 44 under drought, and yielded 157% more than Pusa 44.
Figure 3. Comparative field view of Pusa 44 near isogenic line P1823-12-82 carrying two QTLs for reproductive stage drought stress tolerance under drought-stressed conditions. Note that P1823-12-82 had better grain filling and spikelet fertility than Pusa 44 under drought, and yielded 157% more than Pusa 44.
Agriculture 11 00064 g003
Figure 4. Heatmap of relations between the selected NILs for yield under drought stress, at two test locations, Hazaribag and Rewa, under varietal testing.
Figure 4. Heatmap of relations between the selected NILs for yield under drought stress, at two test locations, Hazaribag and Rewa, under varietal testing.
Agriculture 11 00064 g004
Table 1. qDTY3.1 and qDTY2.1 linked markers for foreground selection (FS) and recombinant selection (RS).
Table 1. qDTY3.1 and qDTY2.1 linked markers for foreground selection (FS) and recombinant selection (RS).
QTLMarkerChromosomePhysical LocationTypeReference
qDTY3.1RM520330.71 MbFS[16]
RM15791 28.56 MbRS
RM16033 32.56 MbRS
qDTY2.1RM521210.8 MbFS[16]
RM5791 10.74 MbRS
RM324 11.4 MbRS
FS, foreground selection; RS, recombinant selection.
Table 2. Progressive selection statistics of backcross generations and the recovery of the recurrent parent genome during the development of Pusa 44 near isogenic lines.
Table 2. Progressive selection statistics of backcross generations and the recovery of the recurrent parent genome during the development of Pusa 44 near isogenic lines.
GenerationNo. of Plants RaisedQTL Positive Progenies Plants SelectedRPG Recovery (%)
F11192*
BC1F1488273.1–82.9
BC2F1323286.1–95.9
BC3F1181198.6
BC3F2400545498.6–99.0
BC3F39488 (54 families)24210898.8–99.4
BC3F4108 NILs108 NILs31 NILs>99.6
* Not estimated; RPG, recurrent parent genome; number of plants/ lines positive for both qDTY3.1 and qDTY2.1.
Table 3. The analysis of variance (ANOVA) for agronomic and grain quality traits showing significant components of variance under stressed and unstressed treatments in BC3F4 generation.
Table 3. The analysis of variance (ANOVA) for agronomic and grain quality traits showing significant components of variance under stressed and unstressed treatments in BC3F4 generation.
TraitVariance under Stress (Vs)Variance under Unstress (Vus)(Vs/Vus)NIL
NILsChecksNIL × CheckNILsChecksNIL × Check
DtF2.40 *906.89 *41.08 *2.20 *844.15 *91.01 *1.09
PtH6.90 ns1315.62 *1487.98 *8.64 *1861.47 *2818.28 *0.80
NpT2.56 *22.72 *1.61 *1.92 *11.60 *0.96 ns1.33
PnL1.07 ns16.57 *41.66 *0.97 ns25.06 *3.39 *1.10
BmP110.34 *611.94 *1998.34 *69.36 *181.41 *3.34 ns1.59
PnW29.92 ns156.80 *65.03 ns30.68 *158.37 *811.58 *0.98
FdG668.29 *1610.25 *28,600.08 *535.97 *9693.43 *1771.72 *1.25
SpF55.60 *248.72 *241.68 *37.45 *43.02 *1059.19 *1.48
GrY9.79 *11.00 *511.57 *11.22 *63.72 *364.31 *0.87
GrW1.50 *45.37 *27.56 *2.66 *14.54 *31.54 *0.56
HgP84.27 *12.87 *75.44 *0.73 *7.92 *0.09 ns115.44
MgP60.58 *42.97 *242.92 *1.31 *6.58 *11.06 *46.24
KLBC0.04 *9.01 *0.00 ns0.03 *9.36 *0.85 *1.33
KWBC0.01 *0.15 *0.02 *0.15 *0.17 *0.16 *0.07
LWR0.03 *1.87 *0.02 *0.06 *1.78 *0.01 ns0.50
KLAC0.12 *21.37 *0.19 *0.94 *23.48 *0.09 ns0.13
KWAC0.04 *0.54 *1.19 *0.13 *0.84 *1.50 *0.31
ER0.01 *0.02 *0.00 ns0.03 *0.01 *0.09 *0.33
DtF, days to 50% flowering; PtH, plant height in cm; NpT, number of panicle-bearing tillers; PnL, length of panicle in cm; BmP, biomass per plant in g; PnW, panicle weight in g; FdG, number of whole grains per panicle; SpF, spikelet fertility in %; GrY, grain yield per plant in g; GrW, weight of 1000 grains in g; HgP, hulling percent; MgP, milling percentage; KLBC, kernel length before cooking in mm; KWBC, kernel width before cooking in mm; LWR, length/width ratio; KLAC, kernel length after cooking in mm; KWAC, kernel width after cooking in mm; ER, elongation ratio; * Significant at 5% level; ns, non-significant.
Table 4. Comparison of average agronomic and grain quality of the NILs vis-à-vis parents in the BC3F4 generation under stressed and unstressed treatments.
Table 4. Comparison of average agronomic and grain quality of the NILs vis-à-vis parents in the BC3F4 generation under stressed and unstressed treatments.
TraitEnvRP (P1)DP (P2)NILsRUsCVLSD
MeanRange
DtFUS108.5120.3110.5106.2–114.2−2.440.61.7
S111.6122.6113.2108.7–115.2-0.51.3
PtHUS97.8136.6102.692.7–109.94.291.23.2
S93.3126.798.288.4–102.4-2.56.1
NpTUS13.711.411.78.2–14.615.387.32.1
S9.08.09.97.2–13.4-4.61.1
PnLUS28.228.326.524.3–29.1−0.382.61.7
S26.327.026.623.7–30.1-3.22.1
GrYUS25.119.924.615.2–32.932.936.84.2
S10.213.416.57.7–26.6-4.01.7
FdGUS202.6168.5145.475.6–194.5−7.697.236.0
S128.6123.2156.688.9–215.3 4.724.0
SpFUS89.688.181.764.3–92.72.943.77.6
S67.973.579.354.1–93.8-2.65.2
GrWUS22.420.522.217.4–25.216.673.41.9
S18.216.518.515.0–21.2-4.32.0
HgPUS79.778.378.575.9–81.42.550.30.6
S76.779.976.571.4–81.6-0.40.8
MgPUS72.871.772.267.1–75.16.230.50.9
S67.072.867.762.4–72.4-1.83.0
KLBCUS6.36.95.95.5–6.35.081.10.2
S6.26.75.65.1–6.3-2.10.3
LWRUS3.43.53.02.8–3.4−3.331.90.1
S3.53.73.12.7–3.4-1.60.1
KLACUS10.711.29.59.0–10.63.162.10.5
S10.110.69.28.3–9.8-1.10.3
ERUS1.71.61.61.4–1.8−6.252.20.1
S1.61.61.71.5–1.9-2.10.1
RP, recurrent parent (Pusa 44); DP, donor parent (IR81896-B-B-142); Env, environments (stressed and unstressed); DtF, days to 50% flowering; PtH, plant height in cm; NpT, number of panicle-bearing tillers; PnL, length of panicle in cm; BmP, biomass per plant in g; PnW, panicle weight in g; FdG, number of whole grains per panicle; SpF, spikelet fertility in %; GrY, grain yield per plant in g; GrW, weight of 1000 grains in g; MgP, milling percentage; HgP, hulling percent; KLBC, kernel length before cooking in mm; LWR, length/width ratio; KLAC, kernel length after cooking in mm; EgR, elongation ratio; RUs, relative reduction under drought over unstressed in % (RUs = (trait under stress—trait under unstress) × 100/trait under unstress); CV, coefficient of variation; SE, standard error of difference; LSD, least significant difference at 5% level of significance.
Table 5. Statistics of yield performance of Pusa 44 NILs along with parents and checks during Kharif 2017 under stressed (S) and unstressed (US) conditions.
Table 5. Statistics of yield performance of Pusa 44 NILs along with parents and checks during Kharif 2017 under stressed (S) and unstressed (US) conditions.
TraitEnvRP (P1)DP (P2)CH1CH2NILsRangeCV(RUs)NILPr
YdPUS6186.14581.25427.44014.25747.0 *4235.0–6543.08.688.00.45 **
S466.7766.71050.0416.7690.8 **410.0–1200.018.9-0.34 **
GrYUS21.414.414.015.718.3 *15.1–23.412.379.20.54 **
S2.95.06.23.43.8 **1.7–5.923.7-0.46 **
DtFUS113.5120.0119.0100.5113.9 **111.0–116.00.5−4.00.85 **
S119.5123.5121.5108.0118.4 **117.0–120.50.5-0.82 **
SpFUS84.079.381.678.879.5 **66.7–88.12.229.90.81 **
S48.761.259.543.055.7 **43.8–79.311.3-0.81 **
YdRS86.765.756.178.478.9----
RP, recurrent parent (Pusa 44); DP, donor parent (IR81896-B-B-142); CH1, check 1 (IR81896-B-B-195); CH2, check 2 (IR64); Env, environments (stressed and unstressed); YdP, grain yield in kg/ha; GrY, grain yield per plant in g; DtF, days to 50% flowering; SpF, spikelet fertility in %; YdR, percentage of yield reduction under stress; (RUs)NIL, relative reduction under drought over unstressed among NILs expressed in % ((RUs)NIL = (trait under stress—trait under unstress) × 100/trait under unstress); Pr, Pearson correlation between 2016 and 2017 data; CV, coefficient of variation; *,** Significant at 5% and 1% P levels, respectively.
Table 6. Drought indices and yield performance of the Pusa 44 qDTY NILs during Kharif 2017.
Table 6. Drought indices and yield performance of the Pusa 44 qDTY NILs during Kharif 2017.
ENTRIESGrY YdP DTERDYDSI
USSUSS
P1823-12-116.84 d–i4.80 a–f6176.00 a–d816.66 b–e28.5071.500.71
P1823-12-1018.81 b–h5.65 abc6330.65 ab766.66 c–f30.0469.960.70
P1823-12-1219.18 a–g2.20 ij6235.99 abc700.00 d–h11.4788.530.89
P1823-12-1418.40 b–i4.85 a–f6539.12 a733.33 c–g26.3673.640.74
P1823-12-1917.39 c–i4.35 a–g5851.04 a–e616.66 e–i25.0174.990.75
P1823-12-3018.56 b–i4.00 c–i6064.90 a–d666.66 d–i21.5578.450.78
P1823-12-3323.45 a3.00 f–j6543.06 a750.00 c–f12.7987.210.87
P1823-12-3517.73 c–i3.50 d–j6090.25 a–d733.33 c–g19.7480.260.80
P1823-12-3816.86 d–i2.45 hij5324.64 c–f716.66 d–h14.5385.470.85
P1823-12-4215.11 ghi3.20 e–j5442.00 b–f716.66 d–h21.1878.820.79
P1823-12-4519.20 a–g3.15 e–j5850.59 a–e750.00 c–f16.4183.590.84
P1823-12-5019.25 a–g4.25 b–h6335.81 ab816.66 b–e22.0877.920.78
P1823-12-5320.94 a–d3.75 d–i5959.70 a–e733.33 c–g17.9182.090.82
P1823-12-5518.05 b–i4.00 c–i6153.59 a–d666.67 d–i22.1677.840.78
P1823-12-6215.42 ghi3.25 e–j5279.02 c–f816.66 b–e21.0878.920.79
P1823-12-6320.68 a–d3.95 c–i5979.95 a–e616.66 e–i19.1080.900.81
P1823-12-6420.48 a–e4.00 c–i5594.44 a–e650.00 d–i19.5380.470.80
P1823-12-6620.31 a–e3.15 e–j5971.05 a–e700.00 d–h15.5184.490.84
P1823-12-6822.30 ab3.90 c–i6272.35 abc600.00 e–i17.4982.510.83
P1823-12-6916.76 d–i1.70 j6202.13 a–d433.33 i10.1489.860.90
P1823-12-7220.61 a–e4.00 c–i6059.02 a–d483.33 ghi19.4180.590.81
P1823-12-7620.05 a–f3.17 e–j6222.65 abc533.33 f–i15.8184.190.84
P1823-12-7718.23 b–i5.35 a–d5492.99 b–f883.33 bcd29.3570.650.71
P1823-12-8017.23 c–i2.40 hij5917.69 a–e410.00 i13.9386.070.86
P1823-12-8120.02 a–f3.15 e–j5605.51 a–e433.33 i15.7384.270.84
P1823-12-8219.44 a–g5.90 ab5454.75 b–f1200.00 a30.3569.650.70
P1823-12-8316.77 d–i3.15 e–j5208.29 d–g600.00 e–i18.7881.220.81
P1823-12-12215.24 ghi4.55 a–g4235.52 gh983.33 abc29.8670.140.70
P1823-12-13315.24 ghi4.25 b–h5573.36 a–f666.67 d–i27.8972.110.72
P1823-12-13716.06 e–i3.40 e–j5028.36 efg733.33 c–g21.1778.830.79
P1823-12-14019.14 a–h4.85 a–e5943.33 a–e550.00 f–i25.3474.660.75
PUSA 4421.43 abc2.85 g–j6186.13 a–d466.67 hi13.3086.700.87
IR81896-B-B-14214.41 hi4.95 a–e4581.15 fgh766.66 c–f34.3565.650.66
IR81896-B-B-19514.00 i6.15 a5427.44 b–f1050.00 ab43.9356.070.56
IR6415.71 f–i3.40 ef4014.19 h416.67 i21.6478.360.78
GrY, grain yield per plant in grams; YdP, grain yield in kg per hectare; DTE, drought tolerance efficiency in percent; RDY, relative decrease in yield in percent; DSI, drought susceptibility index; US, unstressed; S, stressed. Means followed by similar letters are statistically non-significant at 5% level of confidence as detected by Tukey’s honestly significant difference (HSD) test.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Dwivedi, P.; Ramawat, N.; Dhawan, G.; Gopala Krishnan, S.; Vinod, K.K.; Singh, M.P.; Nagarajan, M.; Bhowmick, P.K.; Mandal, N.P.; Perraju, P.; et al. Drought Tolerant near Isogenic Lines (NILs) of Pusa 44 Developed through Marker Assisted Introgression of qDTY2.1 and qDTY3.1 Enhances Yield under Reproductive Stage Drought Stress. Agriculture 2021, 11, 64. https://0-doi-org.brum.beds.ac.uk/10.3390/agriculture11010064

AMA Style

Dwivedi P, Ramawat N, Dhawan G, Gopala Krishnan S, Vinod KK, Singh MP, Nagarajan M, Bhowmick PK, Mandal NP, Perraju P, et al. Drought Tolerant near Isogenic Lines (NILs) of Pusa 44 Developed through Marker Assisted Introgression of qDTY2.1 and qDTY3.1 Enhances Yield under Reproductive Stage Drought Stress. Agriculture. 2021; 11(1):64. https://0-doi-org.brum.beds.ac.uk/10.3390/agriculture11010064

Chicago/Turabian Style

Dwivedi, Priyanka, Naleeni Ramawat, Gaurav Dhawan, Subbaiyan Gopala Krishnan, Kunnummal Kurungara Vinod, Madan Pal Singh, Mariappan Nagarajan, Prolay Kumar Bhowmick, Nimai Prasad Mandal, Puvvada Perraju, and et al. 2021. "Drought Tolerant near Isogenic Lines (NILs) of Pusa 44 Developed through Marker Assisted Introgression of qDTY2.1 and qDTY3.1 Enhances Yield under Reproductive Stage Drought Stress" Agriculture 11, no. 1: 64. https://0-doi-org.brum.beds.ac.uk/10.3390/agriculture11010064

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop