1. Introduction
Rice is grown in more than 100 countries across the world on a combined total of 162.06 million hectares [
1]. Prominently cultivated in Asia, which accounts for about 90% of total world cultivation, rice is grown in varying ecologies ranging from upland, hill, lowland, and deep water. Grown either under irrigated or rainfed conditions, rice suffers yield loss up to 65.0% under rainfed upland and rainfed lowland ecosystems due to the intermittent occurrence seasonal drought stress. Besides, significant yield reductions occur under prolonged severe stress.
Even though rice suffers drought stress in all developmental stages, stress at the reproductive stage is particularly significant in causing yield loss. Reproductive stage drought stress (RSDS) results in shorter panicles with incomplete panicle exsertion and increased spikelet sterility due to pollen and stigma desiccation. Stress during grain filling reduces grain weight, and increases chalkiness due to poor starch granule packing. Therefore, RSDS remains a major challenge for rice production in major rice growing zones [
2]. At the seedling stage, however, drought stress leads to poor seedling establishment, growth retardation, very slow growth rate, and poor tiller initiation along with reduced plant biomass development [
3]. Drought stress during the vegetative growth stage, causes plant height reduction with less effective number of tillers per hill, reduced leaf area, and general stunting [
4].
Rice exhibits significant heritable variation in drought response, with genotypes behaving differently in different rice ecosystems. This diversity in response is wide between the rainfed upland and irrigated lowland systems. Drought response is often related to the severity of drought, with mild stress commonly well-tolerated by most of the genotypes that are either tolerant or moderately tolerant. Severe stress although affects the yield significantly in all genotypes, drought sensitive genotypes suffers the most than the tolerant ones. Several traits have been shortlisted by the breeders as key indicators of selection for drought tolerance. Pinheiro et al. [
5] identified uppermost internode elongation, flag leaf traits such as green color, carbon dioxide assimilation, chlorophyll content, and photosynthesis rate as traits significantly correlated with grain yield, panicle length, filled grain per panicle, spikelet fertility, and 1000 grain weight, making these traits indicators for selection under drought stress. Other important traits associated with drought adaptation are a deeper root system, profuse tillering, thicker wax layer, and wider leaf area. Additionally, physiological traits such as higher relative water content and high membrane stability index are used as indicators for cellular level drought tolerance [
6,
7].
Attempts to breed for RSDS tolerance in rice through conventional means have resulted in limited success due to poor selection efficiency for this trait, confounded by low heritability and genotype x environment interaction. However, significant progress could be achieved through molecular marker-based interventions, wherein several quantitative trait loci (QTLs) governing grain yield under drought (qDTY) have been mapped [
8]. Among these, prominent QTLs such as
qDTY1.1 [
9,
10],
qDTY3.1 [
11] and
qDTY12.1 [
12] have been used for marker-assisted improvement of several elite rice cultivars that are drought sensitive. Currently, improved cultivars having qDTYs either singly or in combination have already been released in several countries [
13,
14].
In India, Pusa 44 is a medium- to long-duration and highly productive mega variety cultivated across the northwestern plain zone, particularly in Punjab and Haryana. Pusa 44 was developed and released by the Indian Council of Agricultural Research (ICAR)-Indian Agricultural Research Institute (IARI), New Delhi solely for the irrigated lowland ecology [
15]. Catering to the livelihood of thousands of rice farmers, Pusa 44 is grown extensively in northwestern India, and is widely popular due to its high yield potential (7.0–8.0 t.ha
−1). A non-aromatic rice variety with translucent long slender grains, Pusa 44 has semi-dwarf stature, strong culm with erect plant type, suitability to mechanical harvesting, medium duration (145–150 days), higher number of productive tillers, and superior cooking quality. However, owing to its high susceptibility to drought stress, especially for RSDS, the continuation of its cultivation has become a major topic of debate in northwestern India, where the receding groundwater level poses a serious challenge to rice production. This has led to the initiation of research to augment Pusa 44 with RSDS tolerance by targeted introgression of qDTYs, either as mono-QTL near isogenic lines (NILS) or QTL stacked NILs. There are recent reports of successful introgression of qDTYs such as
qDTY2.1 and
qDTY3.1 into Pusa 44 [
15,
16]. Introgressing
qDTY12.1, a robust QTL, into Pusa 44 could also be another step toward redesigning this popular cultivar to combat RSDS. This would enable the development of improved lines with a minimum yield potential of one to two tonnes per hectare (t.ha
−1) under severe RSDS, while ensuring the higher yield potential under irrigated lowland conditions [
17]. Cultivars with potential for dual or multiple adaptation are required in the future in areas where climate is fast changing. Therefore, climate resilience will be the key attribute contributing to sustained crop production and ensuring food security under changing climate. Moreover, accelerated development of cultivars using the marker-assisted breeding strategy would be a boon over classical approaches, which normally take 10–12 years for cultivar development.
In the present study, we report the successful introgression of qDTY12.1 into Pusa 44 through marker-assisted backcross breeding, leading to the development of improved NILs showing improved adaptation to RSDS under two contrasting environments in India.
3. Results
3.1. Marker Polymorphism between the Parents
To identify markers for foreground selection, five SSR markers linked to the QTL
qDTY12.1, namely, RM28099, RM511, RM28130, RM28166, and RM28199, were screened for parental polymorphism; of which, RM28130 showed better resolution between parents, making it the ideal choice for further foreground selections. All other markers were found to be monomorphic between parents. RM28130 produced an amplicon size of 176 bp in the donor parent, while the amplicon size in the recurrent parent, Pusa 44, was 180 bp. For background selection, parental polymorphism was carried out using 844 genome-wide SSR markers. The background genomes of the parents were found to be diverse for 84 markers, indicating a diversity of 10.0% between them (
Table 1). Chromosome-wise genome-wide polymorphism indicated that chromosome 3 showed the maximum diversity of 19.4%, and chromosome 1 was the least polymorphic (5.4%). Chromosome 12, the carrier chromosome of
qDTY12.1, had a polymorphism of 10.8%. Among the 12 chromosomes, as many as six showed a diversity above 10%.
3.2. Development of NILs through Marker Assisted Breeding
The hybridization of parents Pusa 44 and IR90019-22-28-2B was attempted during
Rabi 2014–2015 season at the offseason nursery at RBGRC, ICAR-IARI, Aduthurai to generate 200 F
1 seeds. The cross was named Pusa 3003. A total of 106 F
1 plants were raised at ICAR-IARI, New Delhi during
Kharif 2015 and were subjected to hybridity testing using the foreground marker, RM 28130, at the active tillering stage (
Table 2). Thirty-six true F
1 plants that showed clear heterozygous banding for RM28130 were tagged for backcrossing with the recurrent parent, Pusa 44. The BC
1F
1 seeds generated from one of the F
1 plants were raised in
Rabi 2015–2016 at Aduthurai. Foreground selection was carried out among 18 BC
1F
1′s to identify five plants heterozygous for the
qDTY12.1 linked marker RM28130, which were subjected to phenotypic selection. The plants with maximum phenotypic similarity to Pusa 44 were used for generating the BC
2F
1 seeds. Of the 72 BC
2F
1 plants raised during the
Kharif 2016 at New Delhi, 31 were selected based on the foreground and phenotypic selection for further backcrossing with Pusa 44. Thirty BC
3F
1s were raised during
Rabi 2016–2017 at Aduthurai. Foreground selection among the BC
3F
1 plants could identify seven plants heterozygous for the QTL linked SSR marker RM28130, and higher phenotypic resemblance with Pusa 44, after rigorous phenotypic selection. The selection parameters included visual assessment of plant statures, such as plant height, tillering habit, grain morphology, panicle architecture, and grain yield. All the selected plants were selfed during the flowering stage and harvested individually. In
Kharif 2017, six BC
3F
2 populations derived from the best BC
3F
1 plants were raised at New Delhi. Among these, 67 plants were selected from across the six families by foreground selection, wherein only plants with homozygous donor (IR90019-22-28-2B) alleles (176 bp) were selected for generation advancement. The 67 BC
3F
3 families were grown at Aduthurai during the
Rabi 2017–2018, from which progenies similar to Pusa 44 were shortlisted for further evaluation of agronomic performance, resulting in the selection of 18 plants. All 18 selected lines carried the donor allele for
qDTY12.1 in the homozygous state. Subsequently, the selected lines were raised during
Kharif 2018 at New Delhi and subjected to phenotypic and background selections. All 84 polymorphic markers were used for background selection in 18 NILs, which originated from five BC
3F
3 families. During
Rabi 2018–2019, the BC
3F
5 generation was raised at Aduthurai, where further phenotypic selection for agro-morphological as well as grain and cooking quality was carried out to identify six potential NILs for testing. All the six NILs were subsequently evaluated during
Kharif 2019 in New Delhi. The RPGR of the NILs ranged between 94.7 and 98.7%. By the BC
3F
6 generation, the recovery on chromosome 12 was 100% in as many as 29 NILs out of 41 tested, as determined by the flanking markers (
Figure 3).
3.3. Agronomic Evaluation for Differential Response under Stressed and Unstressed Treatments
Evaluation of the NILs for agronomic traits under irrigated and imposed drought conditions for two consecutive years revealed significant variation in components such as genotype and genotype x treatment for most of the traits (
Table 3), except for days to 50% flowering and effective tiller number. During 2018, the genotypic variation for grain yield was not apparent, but between treatments the variation was remarkable. During 2019, however, all three components showed clear variation for grain yield.
3.4. Drought Response of Promising Pusa 44 qDTY12.1-NILs in BC3F4
Compared to the unstressed control, the recurrent parent, Pusa 44, showed significant delay in flowering, with reduced plant height, tiller number, and grain yield under drought stress (
Table 4). The donor parent, IR90019-22-28-2-B, expressed a striking improvement in grain yield compared to other traits such as plant height, panicle length, and tiller number. The yield reduction in Pusa 44 was 9.5 times lower than when irrigated, whereas in IR90019-22-28-2-B it was only about 1.5 times lower than when irrigated. On the other hand, NILs showed a yield reduction level almost similar to that of the donor parent under drought stress conditions, while yielding similarly to Pusa 44 under irrigated conditions, except that five NILs had significantly lower yields under irrigated control. The ratio between unstressed yield to stressed yield among the NILs ranged from 1.9 (Pusa 3003-15-121-9-4) to 3.6 (Pusa 3003-15-121-29-23).
The highest grain yield under drought among the NILs was recorded in Pusa 3003-15-121-30-23 (327.39 g.sqm−1), and all the other NILs were at par with the highest yielding NIL under drought stress. Under irrigated condition, all the NILs, except for Pusa 3003-15-121-9-4, and Pusa 3003-15-121-9-27 were at par with Pusa 44 for grain yield. For other traits, there was no significant variation among the NILs both under stressed and unstressed conditions, except for reduced plant height in Pusa 3003-15-121-31-1 in the unstressed plots.
3.5. Agronomic Performance of BC3F5 Generation NILs under Irrigation
The NILs selected in BC
3F
4 (including those tested under drought as well as those tested only under irrigated conditions in New Delhi) were evaluated at RBGRC, Aduthurai during
Rabi 2018–2019. The data on agronomic performance (
Supplementary Table S1) indicated that all NILs had statistically similar performance for agronomic traits such as days to 50% flowering, plant height, effective tillers number per hill, filled grains per panicle, total number and fertility of spikelets, and grain yield. However, in some of the NILs, the 1000 grain weight was higher than that of Pusa 44. Average grain yield of Pusa 44 was 402.7 g.sqm
−1, while NILs on average produced a grain yield of 384.0 g.sqm
−1.
3.6. Response of BC3F6 Selections under Stressed and Unstressed Conditions
The agronomic performance of the Pusa 44 NILs in the BC
3F
6 generation, both under stressed and unstressed conditions, is presented in
Table 5.
In general, under unstressed treatment, the NILs demonstrated performance at par with Pusa 44, by showing statistically indifferent measurements for most of the agronomic traits. However, agronomic performance of the NILs significantly deviated from that of Pusa 44, particularly for the grain yield. Nevertheless, within each treatment, traits such as days to 50% flowering, plant height, tiller number, panicle length, etc. did not show much variation. Under drought, the days to 50% flowering, however, showed a mixed pattern as one of the parents (Pusa 44) had a delay in flowering under stress, while the other parent had a tendency to flower early under stress. There were nine NILs that showed earliness in flowering under drought, while another nine showed delayed flowering behavior. Among the NILs that showed earliness in flowering, Pusa 3003-15-121-17-47-2 flowered 9.5 days earlier than its irrigated system behavior. This was followed by the NIL, Pusa 3003-15-121-17-47-5, which exhibited an earliness of 4.5 days compared to that under unstressed control. However, in the recurrent parent, the delay in flowering under stress was 11.5 days, whereas the donor parent had closely similar flowering times under both conditions, with a subtle tendency to flower early with a 1.5 days advantage. Grain yields under drought among the NILs varied between 95.0 g.sqm−1 and 295.1 g.sqm−1, but under unstressed conditions, ranged from 589.5 g.sqm−1 to 910.3 g.sqm−1. However, the NILs that produced higher grain yields under irrigation as well as drought were Pusa 3003-15-121-9-6-3, Pusa 3003-15-121-30-23, Pusa 3003-15-121-17-47-2, and Pusa 3003-15-121-9-27 followed by Pusa 3003-15-121-31-1. Among the genotypes tested, the recurrent parent, Pusa 44, was the one that was most affected by drought, whereas the donor parent, IR90019-22-28-2-B, was the least affected. Comparing the ratio between grain yields under unstressed and stressed plots, the largest ratio was shown by the NIL, Pusa 3003-15-121-29-3-1 (8.1), followed by Pusa 44, which had a reduced yield under stress that was eight times lower than its irrigated system yield.
3.7. Indicators of Drought Tolerance among the Pusa 44 qDTY12.1 -NILs
A set of drought response indicators were calculated for assessing the drought tolerance of the NILs in relation to their parents (
Table 6). The cumulative tolerance of each NIL was computed using the average ranks of all the indices. The highest average of 19.1 was recorded for the recurrent parent, Pusa 44 while the donor, IR90019-22-28-2-B recorded a rank average of 4.0. Among the NILs, there were three, Pusa 3003-15-121-17-47-2, Pusa 3003-15-121-9-27, and Pusa 3003-15-121-31-1, that had better average ranks than the donor, while an additional six NILs had average ranks lower than 10.5, the median average rank. Comparing the absolute values for grain yield and average rank, it was found that the yield under stress conditions was highly correlated with a coefficient of −0.958, whereas the yield under irrigated conditions correlated with a coefficient of 0.262. Furthermore, the grain yield of the better ranking lines under drought conditions was found to range between 158.9 and 295.1 g.sqm
−1, whereas their unstressed yield ranged from 592.45 to 847.6 g.sqm
−1.
3.8. Grain Quality of Pusa 44 qDTY12.1 Near Isogenic Lines
The grain quality of the Pusa 44
qDTY12.1 -NILs developed in this study, in BC
3F
6 generation is presented in
Table 7. The parents, Pusa 44 and IR90019-22-28-2B, had uncooked grain lengths of 6.83 and 6.53 mm, respectively, which were statistically at par. Similarly, all the NILs possesses uncooked grain lengths statistically similar with the values for both parents. Additionally, their grain width and the length-width ratios were also found to be non-significantly different from one other. A parallel pattern was observed for grain traits such as length and width when cooked. Nevertheless, only marginal variations between the parents and NILs were observed for grain processing properties such as hulling percentage and milling percentage.
4. Discussion
As hypothesized in this study, future rice cultivars should possess multiple stress tolerance combined with higher yield and desirable grain quality. Such resilient varieties will gain popularity among farmers, as they will contribute to better harvest yields in the event of a climatic anomaly. One of the foremost targets for abiotic tolerance is endurance under RSDS. There are several yield contributing QTLs identified for RSDS tolerance in rice, such as
qDTY12.1, which was established to work under varying genetic backgrounds [
29]. Selection of a consistent QTL is very important to the ultimate success of the backcross breeding program, because introgression of unreliable QTLs/genes may result in considerable losses in time and money. Equally important is the selection of parental lines, because the improved lines are to be used for quick replacement of existing sensitive cultivars, and should enable ready adoption by farmers. In addition, care should be exercised in the choice of the donor, because the ultimately improved line should not carry any undesirable trait from the donor, only the trait desired through introgression. Among undesirable traits, grain quality is a complex of multiple factors that can be affected when the donor is inferior in grain quality. In practical molecular breeding, stringent phenotypic selection has been recognized as a key criterion for imparting all the desirable qualities such as grain quality [
21,
30].
Pusa 44 was selected as the recurrent parent because it has several desirable features such as amenability to mechanical harvest, easily processable grains, and excellent cooking characteristics for consumption. The donor, IR90019-22-28-2-B, which carries the
qDTY12.1 allele, was a derivative of Vandana, an early duration cultivar popular in eastern India. The donor line carried
qDTY12.1 from Way Rarem, an Indonesian cultivar known for its susceptibility to drought [
29,
31]. Mapped earlier in an F
3 population of Vandana/Way Rarem [
12],
qDTY12.1 has been identified as a major QTL contributing as much as 51% of the phenotypic variation for yield under RSDS. This QTL could improve the drought tolerance of an already tolerant cultivar such as Vandana, and is effective under varying genetic backgrounds adapted to different ecosystems such as upland and lowland [
31]. When present,
qDTY12.1 is known to increase tillering, increase in biomass, reduce crop duration, and improve the harvest index [
32]. Furthermore, its prevalence among 85% of the random drought cultivars is remarkable, and evidence indicates that other cereal genomes carry homologous sequences to this QTL [
33]. Due to its versatility,
qDTY12.1 has been a choice QTL for introgression into multiple backgrounds across the globe for improved RSDS tolerance [
14] and is known to integrate high grain yields under mild to severe stress conditions.
We employed MABB for introgression of
qDTY12.1 into Pusa 44. MABB has been established as a successful approach to the introgression of target genes/QTLs, particularly in rice, with high precision, very good success rates, and a short developmental period [
34]. This can be construed from the fact that, within the decade from 2008 to 2018, eight varieties were developed and released by ICAR-IARI using the MABB approach [
30]. Additionally, several other improved lines targeting multiple features are also under development. Among these, Pusa 44 has been a prominent choice among the non-Basmati rice varieties for improvement due to the advantages described earlier. Recent developments in the improvement of Pusa 44 by the introgression of QTLs imparting RSDS tolerance such as
qDTY2.1 and
qDTY3.1 were reported by Dwivedi et al. [
15] and Oo et al. [
16].
Throughout the developmental process for Pusa 44
qDTY12.1 -NILs, we performed foreground selection in all generations but resorted to phenotype selection up until the BC
3F
4 generation. This has been a conspicuous deviation from our earlier approaches wherein background selection was augmented with phenotypic selection to ensure better RPG as well as recurrent parent phenome (RPP) recovery [
34,
35,
36]. Ultimately, deferred background selection provided a recovery of 94.7–98.7%, closely similar to the level of recovery that was previously reported [
37,
38]. This provided us with two advantages: (i) considerable time could be saved on screening the background markers, which helped us to shuttle the breeding materials between New Delhi and Aduthurai in time, and (ii) considerable resources could be saved as the background recovery was tested on a small panel of selected lines at BC
3F
4. Although the approach resulted in gaining the desired RPGR, the uncertainty of gain in RPG remained until the final background selection. This is feasible only when the breeder has excellent expertise in phenotypic selection. Another noticeable feature has been the complete recovery of the target chromosome 12, among several of the NILs.
The NILs generated in this study performed well agronomically under imposed drought stress as well as unstressed conditions. Evaluated for two years during the
Kharif seasons of 2018 and 2019, the NILs showed significant yield reduction under drought stress when compared to yields under unstressed, irrigated conditions. A similar reduction in yield under imposed drought conditions was reported among the qDTY introgressed NILs of mega cultivars such as Pusa Basmati 1 [
39] and Pusa 44 [
15,
16]. However, the level of yield reduction was significantly lower among the NILs, with 6.5 times higher yields than those with the recurrent parent, Pusa 44. On average, the yield under drought conditions among the NILs was 2.3 times lower than that under unstressed conditions. For comparison, in Pusa 44 the yield due to RSDS was 9.7 times lower than the irrigated yield. This indicated that the
qDTY12.1 introgression could manifest consistent and discernible drought tolerance among the improved NILs. Considering the intensity of drought imposed, up to a significant soil moisture tension of −70 kPa, we concluded that the level of tolerance achieved was significant. Ghosh and Singh [
40] reported that in aerobic rice, soil moisture tension of −60 kPa could result in a yield loss of 42.8% and concluded that −40 kPa could be used as the threshold for scheduling irrigation.
The remarkable gain in yield response in the Pusa 44 derived
qDTY12.1 -NILs could not be attributed to component traits. Comparing agronomic performance across years and between stressed and unstressed treatments, the only trait that showed striking variation was plant height. Plant height was, however, found to decrease under drought in all the genotypes, including the recurrent as well as donor parents. However, a distinguishable pattern in the reduction in height could not be drawn between the sensitive parent, Pusa 44, and the NILs. This could imply that the effect of
qDTY12.1 in imparting drought endurance could vary between the genetic backgrounds. Similar differences in the agronomic advantage of
qDTY12.1 have already been reported [
41]. When compared to agronomic changes that are conditioned by different qDTY QTLs in rice, we could see that plant height reduction was a common response under drought conditions, as reported in previous works, wherein
qDTY2.1 and
qDTY3.1 were introgressed into Pusa 44 [
15,
16] and Pusa Basmati 1 [
39]. In addition, a delay in flowering and a reduction in the number of tillers were also reported in the improved NILs. Interestingly, in the present case, the Pusa 44
qDTY12.1 NILs did not show any significant delay in flowering or reduction in tiller numbers. These observations provide possible support to the earlier reports that
qDTY12.1 improved tillering and flowering response in rice [
29]. Furthermore, a remarkable recovery of grain quality in terms of cooking properties, grain dimensions, and processing properties has been realized among the Pusa 44
qDTY12.1 lines (
Figure 4).
The present investigation identified some of the outstanding NILs, with potential for improved yields under both drought and irrigated environments. NILs such as Pusa 3003-15-121-17-47-2, Pusa 3003-15-121-9-27, Pusa 3003-15-121-31-1, Pusa 3003-15-121-30-23, and Pusa 3003-15-121-9-6-3 have better drought response indices when compared to the recurrent parent, Pusa 44. Drought response indices used in this study placed significantly more weight on drought response than yield under normal management. This was evident from the strong association between the average rank and yield under drought conditions (−0.958). The correlation to the yield under irrigated conditions was 0.262. This provided us with the opportunity to select the best drought-tolerant NIL based on the lowest average rank. In breeding for drought tolerance, selection of the line with the best drought tolerance should be given preference over yield potential in irrigated environments to ensure reasonable grain yield in the event of a severe RSDS episode.
Grain quality remains an integral goal of crop improvement, especially for a crop such as rice, irrespective of the target trait of introgression as rice is consumed with minimal processing. Desirable grain quality in the improved NILs was accomplished through rigorous selection for quality components during the breeding process. It is noteworthy that the processing properties of the improved lines did not deteriorate under stressed treatment. However, a slight reduction in grain dimensions was noticed among the NILs, similarly to the reduction noticed in Pusa 44. A similar reduction in grain dimensions was noticed by Dhawan et al. [
39] among the Pusa Basmati 1
qDTY1.1 NILs. Desired grain quality in the selected NILs would facilitate their ready adoption by farmers who are already familiar with the grain quality of Pusa 44, augmented by the potential of yielding well under RSDS.