Next Article in Journal
Mitogenomes Reveal Two Major Influxes of Papuan Ancestry across Wallacea Following the Last Glacial Maximum and Austronesian Contact
Next Article in Special Issue
Identification and Validation of a QTL for Bacterial Leaf Streak Resistance in Rice (Oryza sativa L.) against Thai Xoc Strains
Previous Article in Journal
Open Issues for Protein Function Assignment in Haloferax volcanii and Other Halophilic Archaea
Previous Article in Special Issue
Genome-Wide Analysis of Potassium Channel Genes in Rice: Expression of the OsAKT and OsKAT Genes under Salt Stress
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Molecular Breeding for Improving Productivity of Oryza sativa L. cv. Pusa 44 under Reproductive Stage Drought Stress through Introgression of a Major QTL, qDTY12.1

by
Kyaw Swar Oo
1,
Subbaiyan Gopala Krishnan
1,*,
Kunnummal Kurungara Vinod
1,
Gaurav Dhawan
1,
Priyanka Dwivedi
1,
Pankaj Kumar
1,
Prolay Kumar Bhowmick
1,
Madan Pal
2,
Viswanathan Chinnuswamy
2,
Mariappan Nagarajan
3,
Haritha Bollinedi
1,
Ranjith Kumar Ellur
1 and
Ashok Kumar Singh
1
1
Division of Genetics, ICAR-Indian Agricultural Research Institute, New Delhi 110012, India
2
Division of Plant Physiology, ICAR-Indian Agricultural Research Institute, New Delhi 110012, India
3
Rice Breeding and Genetics Research Centre, ICAR-Indian Agricultural Research Institute, Aduthurai 612101, India
*
Author to whom correspondence should be addressed.
Genes 2021, 12(7), 967; https://doi.org/10.3390/genes12070967
Submission received: 11 February 2021 / Revised: 20 April 2021 / Accepted: 27 April 2021 / Published: 24 June 2021
(This article belongs to the Special Issue Advances in Rice Genetics and Breeding)

Abstract

:
Increasing rice production is quintessential to the task of sustaining global food security, as a majority of the global population is dependent on rice as its staple dietary cereal. Among the various constraints affecting rice production, reproductive stage drought stress (RSDS) is a major challenge, due to its direct impact on grain yield. Several quantitative trait loci (QTLs) conferring RSDS tolerance have been identified in rice, and qDTY12.1 is one of the major QTLs reported. We report the successful introgression of qDTY12.1 into Pusa 44, a drought sensitive mega rice variety of the northwestern Indian plains. Marker-assisted backcross breeding (MABB) was adopted to transfer qDTY12.1 into Pusa 44 in three backcrosses followed by four generations of pedigree selection, leading to development of improved near isogenic lines (NILs). Having a recurrent parent genome (RPG) recovery ranging from 94.7–98.7%, the improved NILs performed 6.5 times better than Pusa 44 under RSDS, coupled with high yield under normal irrigated conditions. The MABB program has been modified so as to defer background selection until BC3F4 to accelerate generational advancements. Deploying phenotypic selection alone in the early backcross generations could help in the successful recovery of RPG. In addition, the grain quality could be recovered in the improved NILs, leading to superior selections. Owing to their improved adaptation to drought, the release of improved NILs for regions prone to intermittent drought can help enhance rice productivity and production.

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.

2. Materials and Methods

2.1. Plant Materials

The popular cultivar, Oryza sativa L. cv. Pusa 44, which is highly sensitive to RSDS, was used as the recipient parent, the seeds of which were obtained from the Division of Genetics, ICAR-IARI, New Delhi. The donor parent for the QTL, qDTY12.1 was IR90019-22-28-2B developed by the International Rice Research Institute from the cross, Vandana/Way Rarem [18]. All the field experiments were implemented at the research farms of ICAR-IARI at New Delhi and at the Rice Breeding and Genetics Research Centre (RBGRC), ICAR-IARI, Aduthurai, Tamil Nadu. During the development of NILs, the plants were grown with standard recommended agronomic management under irrigated ecology.

2.2. Markers for Selection

Marker-assisted selection was carried out in all generations in two steps, foreground and background. For foreground selection, initially, five peak microsatellite (SSR) markers tightly linked to qDTY12.1 namely, RM28099, RM511, RM28130, RM28166, and RM28199 [19] were tested for polymorphism. Among these, RM28130, which provided the best contrast between parents, was further used for selection (Figure 1). Eighty-four polymorphic SSR markers identified as polymorphic from a set of 844 SSR markers across the rice genome were employed for background selection in the BC3F4 generation (Supplementary Table S1).

2.3. Development of Pusa 44 qDTY12.1 Near Isogenic Lines

The hybridization of Pusa 44 with IR90019-22-28-2B was carried out during Rabi season (December–May) of 2014–2015 at Aduthurai. Pusa 44 was used as the female parent throughout the crossing program (Figure 2). The F1 seeds were grown in New Delhi during the next Kharif season (Jun–October) of 2015 and tested for hybridity using the foreground marker, RM28130. True F1s were backcrossed to the recipient parent, Pusa 44, to develop the BC1F1 generation. Foreground selection was carried out with BC1F1 plants and the true BC1F1s with maximum phenotypic similarity to Pusa 44 were chosen for the second round of backcrossing. The same protocol was continued until the development of BC3F1. In the BC3F1 generation, the progenies heterozygous for the target marker, RM28130 were identified and selfed to generate the BC3F2 population. Based on foreground selection in the BC3F2, plants homozygous for the donor allele were selected and subjected to phenotypic selection to identify desirable plants, which were then selfed to generate the BC3F3 generation. Further, phenotypic selection for similarity to the recurrent parent, Pusa 44 was carried out for agronomic, grain, and cooking quality traits, was carried out in all generations beginning from the BC1F1 generation. Selected BC3F3 plants were advanced to BC3F4 families by selfing. Continuing this process, improved Pusa 44 qDTY12.1 NILs in BC3F6 with maximum of recurrent parent genome recovery (RPGR) as well as phenotypic and grain quality attributes similar to Pusa 44 were identified for further evaluation of agronomic superiority over Pusa 44 under RSDS. The RPGR was computed using the genome-wide polymorphic markers in each of the NILs as described earlier [20,21]. Additionally, recombinant selection in the carrier chromosome was also undertaken using the markers flanking qDTY12.1 in the BC3F6 generation.

2.4. Field Screening under Reproductive Stage Drought Stress

Field screening of the improved NILs of Pusa 44 carrying qDTY12.1 was carried out during two consecutive years at ICAR-IARI, New Delhi. During 2018, 7 BC3F4 lines that showed superior yield and grain quality were selected from a set of 18 lines and raised under artificially created RSDS along with the parents, Pusa 44 and IR90019-22-28-2B. The experiment was conducted in randomized complete block design with two replications under field conditions. A parallel set of the same genotypes were grown under the recommended irrigation schedule. Both treatments were raised with standard agronomic management. Both trials were grown with spacing of 20 cm × 15 cm, with each genotype occupying a 2.4 m2 area. In the stressed plots, drought was imposed 25 days after transplanting (50 days after sowing), and one life-saving irrigation was provided 55 days after transplanting by flood irrigation. The unstressed treatment received normal irrigation and was maintained with standing water. The moisture level in the drought stressed plots was monitored by placing tensiometers at regular intervals. At the reproductive stage, the soil moisture tension was monitored till it reached upto −70 kPa followed by lifesaving irrigation at peak stress of −70 kPa. At physiological maturity, the trials were measured for grain yield traits such as single plant yield, and plot yield as well as for agronomic traits such as plant height, number of tillers per hill, and panicle length. The remaining lines in the BC3F4 generation were grown only under irrigation. The selections from the BC3F4 lines after agronomic and grain quality evaluation were sent to RBGRC, ICAR-IARI, Aduthurai, where the lines were grown under irrigated conditions. Further, selections from this trial (Supplementary Table S1) were evaluated at ICAR-IARI, New Delhi during Kharif 2019, subjecting the lines under both imposed drought as well as irrigated conditions. There were 18 NILs selected under this trial, coinciding with the BC3F6 generation.
Additionally, to quantify the stress response of the NILs, stress response related indices [16,22] such as relative stress tolerance (RTOL) [23], geometric mean productivity (GMP), stress-tolerance index (STI) [24], stress susceptibility index (SSI) [25], harmonic mean [26], yield index (YI) [27], and yield stability index [28], were computed for comparative evaluation.

2.5. Grain Quality Analysis

The grain quality of the Pusa 44 qDTY12.1 -NILs was tested under both stressed and unstressed conditions to check if the drought has any influence on the quality parameters. For this, grains were tested for hulling and milling recovery and cooking parameters such as grain dimensions (length and width) before and after cooking, length to width ratio, as well as kernel elongation ratio on cooking.

2.6. Statistical Analyses

All the data recorded from the replicated trials were subjected to analysis of variance (ANOVA) for all the agronomic and grain quality traits. Mean comparison was carried out between the irrigated control and drought stress screening plots.

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 F1 seeds. The cross was named Pusa 3003. A total of 106 F1 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 F1 plants that showed clear heterozygous banding for RM28130 were tagged for backcrossing with the recurrent parent, Pusa 44. The BC1F1 seeds generated from one of the F1 plants were raised in Rabi 2015–2016 at Aduthurai. Foreground selection was carried out among 18 BC1F1′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 BC2F1 seeds. Of the 72 BC2F1 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 BC3F1s were raised during Rabi 2016–2017 at Aduthurai. Foreground selection among the BC3F1 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 BC3F2 populations derived from the best BC3F1 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 BC3F3 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 BC3F3 families. During Rabi 2018–2019, the BC3F5 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 BC3F6 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 BC3F4 (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 BC3F6 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 BC3F6 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 F3 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 BC3F4 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 BC3F4. 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.

5. Conclusions

Through this study, we were able to establish a modified MABB approach with background selection deferred, allowing successful introgression of qDTY12.1 without any undesired effects. The improved grain yield under drought conditions in the NILs corroborated the findings of earlier reports that emphasized compatibility with wider genetic backgrounds belonging to either lowland or upland ecosystems [32]. We presume that the addition of this QTL could leverage the adaptation of Pusa 44 to the changing dynamics of water availability for rice production in northwestern India. That assessment notwithstanding, the consistent yield performance of the improved NILs qualifies them as candidates for pan-India multi-location testing under the all-India coordinated rice improvement program. As the improved Pusa 44 qDTY12.1 -NILs is in an elite background, they can be used also as donors for future RSDS improvements using diverse cultivar backgrounds. The enhanced tolerance of the improved NILs to RSDS may facilitate their adoption elsewhere in the country beyond current niche districts: for example, in eastern India, where intermittent episodes of RSDS are a regular feature. We believe that the development of these NILs could be a step forward in achieving climate resilience in future rice cultivars.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/genes12070967/s1, Table S1: Agronomic characteristics of Pusa 44 qDTY12.1 near isogenic lines during Rabi 2018–2019 at Aduthurai.

Author Contributions

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

Funding

This research received funding from the ICAR funded network project on “National Innovations on Climate Resilient Agriculture (NICRA)”and the 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 Indian Council of Agricultural Research.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data used in this study are presented in the supplementary information.

Acknowledgments

K.S.O. acknowledges the financial support he received from the Indian Council of Cultural Relations, Government of India during his doctoral program at ICAR-IARI, New Delhi. He also acknowledges the Advanced Centre for Agricultural Research and Education, Government of Myanmar for providing his deputation and scholarship during the program.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

ANOVAAnalysis of variance
BCBackcross
bpBase pairs
BSBackground selection
CDCritical difference
CmNNumber of tillers (culms)
CVCoefficient of variation
DFFDays to 50% flowering
ERElongation ratio
FSForeground selection
g.sqm−1Grams per square meter
GMPGeometric mean productivity
GQSelection for grain and cooking qualities
GYGrain yield
HAMHarmonic mean
HtPlant height
HULLHulling recovery
IARIIndian Agricultural Research Institute
ICARIndian Council of Agricultural Research
KERKernel elongation ratio
KLACKernel length after cooking
KLBCKernel length before cooking
kPaKilo Pascal
LBRLength by breadth ratio
MABBMarker assisted backcross breeding
MILLMilling recovery
NILNear isogenic line
PnLPanicle length
PSPhenotypic selection
qDTYQTL for yield under drought
QTLQuantitative trait loci
RBGRCRice breeding and genetics research centre
RMRice microsatellite
RPRecurrent parent
RPGRecurrent parent genome
RPPRecurrent parent phenome
RPGRRecurrent parent genome recovery
RSDSReproductive stage drought stress
RTOLRelative stress tolerance
SStressed
SFSpikelet fertility
SSIStress susceptibility index
SSRSimple sequence repeats
STIStress tolerance index
t.ha−1tonnes per hectare
USUnstressed
YIYield index
YldGrain yield per plot
YSIYield stability index

References

  1. FAO. FAOSTAT. Food and Agriculture Organization of the United Nations: Rome, Italy, 2020. Available online: http://fao.org/faostat/en/ (accessed on 20 January 2021).
  2. Feng, D.; Wang, Y.; Wu, J.; Lu, T.; Zhang, Z. Development and drought tolerance assay of marker-free transgenic rice with OsAPX2 using biolistic particle-mediated co-transformation. Crop J. 2017, 5, 271–281. [Google Scholar] [CrossRef] [Green Version]
  3. Menguer, P.K.; Sperotto, R.A.; Ricachenevsky, F.K. A walk on the wild side: Oryza species as source for rice abiotic stress tolerance. Genet. Mol. Biol. 2017, 40, 238–252. [Google Scholar] [CrossRef] [Green Version]
  4. Moonmoon, S.; Islam, M.T. Effect of drought stress at different growth stages on yield and yield a component of six rice (Oryza sativa L.) genotypes. Fundam. Appl. Agric. 2017, 2, 285–289. [Google Scholar] [CrossRef]
  5. Pinheiro, B.D.S.; Austin, R.B.; Carmo, M.P.D.; Hall, M.A. Carbon isotope discrimination and yield of upland rice as affected by drought at flowering. Pesq. Agropec. Bras. 2000, 35, 1939–1947. [Google Scholar] [CrossRef]
  6. Swapna, S.; Shylaraj, S.K. Screening for osmotic stress responses in rice varieties under drought condition. Rice Sci. 2017, 24, 253–263. [Google Scholar] [CrossRef]
  7. Salehifar, M.; Rabiei, B.; Mohammadian, M.A.; Asghari, J. Physiological and fluorescence reaction of four rice genotypes to exogenous application of IAA and Kinetin under drought stress. Not. Sci. Biol. 2017, 9, 378–385. [Google Scholar] [CrossRef] [Green Version]
  8. 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; Rajpal, V.R., Sehgal, D., Kumar, A., Raina, S.N., Eds.; Springer Nature: Basel, Switzerland, 2019; Volume 2. [Google Scholar] [CrossRef]
  9. Vikram, P.; Swamy, B.P.M.; Dixit, S.; Ahmed, H.U.; Cruz, M.T.S.; Singh, A.K.; Kumar, A. qDTY1.1, a major QTL for rice grain yield under reproductive-stage drought stress with a consistent effect in multiple elite genetic backgrounds. BMC Genet. 2011, 12, 89. [Google Scholar] [CrossRef] [Green Version]
  10. Ghimire, K.H.; Quiatchon, L.A.; Vikram, P.; Swamy, B.P.M.; Dixit, S.; Ahmed, H.; Hernandez, J.E.; Borromeo, T.; Kumar, A. Identification and mapping of a QTL (qDTY1.1) with a consistent effect on grain yield under drought. Field Crop Res. 2012, 131, 88–96. [Google Scholar] [CrossRef]
  11. Venuprasad, R.; Dalid, C.O.; Valle, M.D.; 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]
  12. Bernier, J.; Kumar, A.; Ramaiah, V.; Spaner, D.; Atlin, G. A large-effect QTL for grain yield under reproductive-stage drought stress in upland rice. Crop Sci. 2007, 47, 507–516. [Google Scholar] [CrossRef]
  13. 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]
  14. Kumar, A.; Sandhu, N.; Venkateshwarlu, C.; Priyadarshi, R.; Yadav, S.; Majumder, R.R.; Singh, V.K. Development of introgression lines in high yielding, semi-dwarf genetic backgrounds to enable improvement of modern rice varieties for tolerance to multiple abiotic stresses free from undesirable linkage drag. Sci. Rep. 2020, 10, 13073. [Google Scholar] [CrossRef] [PubMed]
  15. Dwivedi, P.; Ramawat, N.; Dhawan, G.; Krishnan, S.G.; 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. [Google Scholar] [CrossRef]
  16. Oo, K.S.; Dhawan, G.; Kumar, P.; Prakash, V.; Krishnan, S.G.; Bhowmick, P.K.; Pal, M.; Nagarajan, M.; Ellur, R.K.; Bollinedi, H.; et al. Development and evaluation of Pusa 44 backcross derived lines possessing qDTY2.1, qDTY3.1 and their combinations indicate differential yield response under severe reproductive stage drought stress. Indian J. Genet. Plant Breed. 2021. Accepted. [Google Scholar]
  17. Shamsudin, N.A.A.; Swamy, B.P.M.; Ratnam, W.; Cruz, M.T.S.; Raman, A.; Kumar, A. Marker assisted pyramiding of drought yield QTLs into a popular Malaysian rice cultivar, MR219. BMC Genet. 2016, 17, 30. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. 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 a QTL for drought resistance in rice, qtl12.1, over a range of environments in the Philippines and eastern India. Euphytica 2009, 166, 207–217. [Google Scholar] [CrossRef]
  19. 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]
  20. Joseph, M.; Gopalakrishnan, S.; Sharma, R.K.; Singh, V.P.; Singh, A.K.; Singh, N.K.; Mohapatra, T. Combining bacterial blight resistance and Basmati quality characteristics by phenotypic and molecular marker-assisted selection in rice. Mol. Breed. 2004, 13, 377–387. [Google Scholar] [CrossRef]
  21. Gopalakrishnan, S.; Sharma, R.K.; Rajkumar, K.A.; 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]
  22. Ghobadi, M.; Ghobadi, M.E.; Kahrizi, D.; Zebarjadi, A.; Geravandi, M. Evaluation of drought tolerance indices in dryland bread wheat genotypes under post-anthesis drought stress. Int. J. Agric. Biosyst. Eng. 2012, 67, 1257–1261. [Google Scholar] [CrossRef]
  23. Rosielle, A.A.; Hamblin, J. Theoretical aspects of selection for yield in stress and non-stress environments. Crop Sci. 1981, 21, 943–946. [Google Scholar] [CrossRef]
  24. Fernandez, G.C.J. Effective selection criteria for assessing plant stress tolerance. In Proceedings of the International Symposium on Adaptation of Vegetable and Other Food Crops in Temperature and Water Stress, Shanhua, Taiwan, 13–18 August 1992; pp. 257–270. [Google Scholar] [CrossRef]
  25. Fischer, R.A.; Maurer, R. Drought resistance in spring wheat cultivars: I. Grain yield response. Aust. J. Agric. Res. 1978, 29, 897–912. [Google Scholar] [CrossRef]
  26. Schneider, K.A.; Rosales-Serna, R.; Ibarra-Perez, F.; Cazares-Enriquez, B.; Acosta-Gallegos, J.A.; Ramirez-Vallego, P.; Wassimi, N.; Kelly, J.D. Improving common bean performance under drought stress. Crop Sci. 1997, 37, 43–50. [Google Scholar] [CrossRef]
  27. Bouslama, M.; Schapaugh, T. Stress tolerance in soybean. Part I: Evaluation of three screening techniques for heat and drought tolerance. Crop Sci. 1984, 24, 933–937. [Google Scholar] [CrossRef]
  28. Lin, C.S.; Binns, M.R.; Lefkovich, L.P. Stability analysis: Where do we stand? Crop Sci. 1986, 26, 894–900. [Google Scholar] [CrossRef] [Green Version]
  29. 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]
  30. Singh, A.K.; Krishnan, S.G.; Vinod, K.K.; Ellur, R.K.; Bollinedi, H.; Bhowmick, P.K.; Nagarajan, M. Precision breeding with genomic tools: A decade long journey of molecular breeding in rice. Indian J. Genet. Plant Breed. 2019, 79, 181–191. [Google Scholar] [CrossRef]
  31. Dixit, S.; Biswal, A.K.; Min, A.; Henry, A.; Oane, R.H.; Raorane, M.L.; Longkumer, T.; Pabuayon, I.M.; Mutte, S.K.; Vardarajan, A.R.; et al. Action of multiple intra-QTL genes concerted around a co-localized transcription factor underpins a large effect QTL. Sci. Rep. 2015, 5, 15183. [Google Scholar] [CrossRef] [Green Version]
  32. Mishra, K.K.; Vikram, P.; Yadaw, R.B.; Swamy, B.P.M.; 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] [Green Version]
  33. Swamy, B.P.M.; Vikram, P.; Dixit, S.; Ahmed, H.U.; Kumar, A. Meta-analysis of grain yield QTL identified during agricultural drought in grasses showed consensus. BMC Genom. 2011, 12, 319. [Google Scholar] [CrossRef] [Green Version]
  34. Ellur, R.; Khanna, A.; Krishnan, S.G.; 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]
  35. Babu, N.N.; Krishnan, S.G.; 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, 41. [Google Scholar] [CrossRef] [Green Version]
  36. Sagar, V.; Dhawan, G.; Krishnan, S.G.; Vinod, K.K.; Ellur, R.E.; 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]
  37. Yadav, A.K.; Kumar, A.; Grover, N.; Ellur, R.K.; Krishnan, S.G.; Bollinedi, H.; Bhowmick, P.K.; Vinod, K.K.; Nagarajan, M.; Krishnamurthy, S.L.; et al. Marker aided introgression of ‘Saltol’, a major QTL for seedling stage salinity tolerance into an elite Basmati rice variety ‘Pusa Basmati 1509’. Sci. Rep. 2020, 10, 13877. [Google Scholar] [CrossRef] [PubMed]
  38. Grover, N.; Kumar, A.; Yadav, A.K.; Krishnan, S.G.; Ellur, R.E.; Bhowmick, P.K.; Vinod, K.K.; Bollinedi, H.; Nagarajan, M.; Viswanthan, C.; et al. Marker assisted development and characterization of herbicide tolerant near isogenic lines of a mega Basmati rice variety, “Pusa Basmati 1121”. Rice 2020, 13, 68. [Google Scholar] [CrossRef] [PubMed]
  39. Dhawan, G.; Kumar, A.; Dwivedi, P.; Krishnan, S.G.; Pal, M.; Vinod, K.K.; Nagarajan, M.; Bhowmick, P.K.; Bollinedi, H.; Ellur, R.K.; et al. Introgression of qDTY1.1 governing reproductive stage drought tolerance into an elite Basmati rice variety “Pusa Basmati 1” through marker assisted backcross breeding. Agronomy 2021, 11, 202. [Google Scholar] [CrossRef]
  40. Ghosh, A.; Singh, O. Determination of threshold regime of soil moisture tension for scheduling irrigation in tropical aerobic rice for optimum crop and water productivity. Exp. Agric. 2010, 46, 489–499. [Google Scholar] [CrossRef]
  41. Vikram, P.; Swamy, B.P.M.; Dixit, S.; Trinidad, J.; Cruz, M.T.S.; Maturan, P.C.; Amante, M.; Kumar, A. Linkages and interactions analysis of major effect drought grain yield QTLs in rice. PLoS ONE 2016, 11, e0151532. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Representative gel picture showing polymorphism of foreground markers. RM28130 provided the best contrast between Pusa 44 and IR90019-22-28-2B. L: ladder; RP: recurrent parent, Pusa 44; DP: donor parent, IR90019-22-28-2B. RP and DP were replicated 3 times.
Figure 1. Representative gel picture showing polymorphism of foreground markers. RM28130 provided the best contrast between Pusa 44 and IR90019-22-28-2B. L: ladder; RP: recurrent parent, Pusa 44; DP: donor parent, IR90019-22-28-2B. RP and DP were replicated 3 times.
Genes 12 00967 g001
Figure 2. Breeding scheme adopted for the development of Pusa 44 qDTY12.1 -NILs.
Figure 2. Breeding scheme adopted for the development of Pusa 44 qDTY12.1 -NILs.
Genes 12 00967 g002
Figure 3. Graphical genotype of Chromosome 12, the carrier chromosome of the QTL, qDTY12.1 among the Pusa44 qDTY12.1-NILs. 29 NILs showed complete recovery of Chromosome 12, based on the polymorphic flanking markers of the QTL.
Figure 3. Graphical genotype of Chromosome 12, the carrier chromosome of the QTL, qDTY12.1 among the Pusa44 qDTY12.1-NILs. 29 NILs showed complete recovery of Chromosome 12, based on the polymorphic flanking markers of the QTL.
Genes 12 00967 g003
Figure 4. Grain morphology of the selected Pusa 44 qDTY12.1 -NILs, showing remarkable recovery of grain and cooking quality.
Figure 4. Grain morphology of the selected Pusa 44 qDTY12.1 -NILs, showing remarkable recovery of grain and cooking quality.
Genes 12 00967 g004
Table 1. Chromosome wise polymorphism between the parents, Pusa 44 and IR90019-22-28-2B.
Table 1. Chromosome wise polymorphism between the parents, Pusa 44 and IR90019-22-28-2B.
ChromosomeMarkers UsedDiversity (%)
PolymorphicMonomorphicTotal
171221295.43
211566716.42
37293619.44
4576816.17
5684906.67
69627112.68
77344117.07
87941016.93
95313613.89
10467715.63
119475616.07
127 *586510.77
Total847608449.95
* The foreground marker, RM28130, is not included in the list of markers used for background selection across chromosome 12.
Table 2. Progressive selection statistics for Pusa 44 qDTY12.1 near isogenic lines.
Table 2. Progressive selection statistics for Pusa 44 qDTY12.1 near isogenic lines.
GenerationNo. of Plants/LinesSelection MethodRPGR %
RaisedSelected
F110636Hybridity testing50%
BC1F1185FS + PS-
BC2F17231FS + PS-
BC3F1307FS + PS-
BC3F26 (Families)67FS + PS-
BC3F367 (Lines)18 *FS + PS-
BC3F418 (Lines)18 NILsPS + BS + GQ94.7–98.7
BC3F518 NILs6 NILs **PS + BS + GQ94.7–98.7
FS, foreground selection; BS, background selection; PS, phenotypic selection; GQ, selection for grain and cooking qualities; NIL, near isogenic line; * 18 plants belonged to 5 families and 15 lines; ** 6 NILs belonged to 4 families and 6 lines.
Table 3. Combined analysis of variance among the Pusa 44 NILs carrying qDTY12.1 under irrigated (US) and drought stress (S) during Kharif 2018 and Kharif 2019.
Table 3. Combined analysis of variance among the Pusa 44 NILs carrying qDTY12.1 under irrigated (US) and drought stress (S) during Kharif 2018 and Kharif 2019.
TraitsTreatmentGenotypeGenotype x TreatmentCV%CD (5%)
Days to 50% flowering10.51 ns42.92 **15.48 **2.8210.03
Plant height6656.87 **89.60 **29.95 **2.948.49
Effective tillers per hill0.007 ns6.37 *6.94 **12.466.58
Panicle length50.99 **3.91 **1.98 **3.482.93
Plot yield Kharif 20181341072.08 **10927.24 ns21542.61 **4.3344.22
Plot yield Kharif 20196174952.63 **10872.42 **23166.12 **11.42110.02
*, ** Significant at 0.05 and 0.01 probability level, respectively; ns, non-significant; CV, coefficient of variation; CD, the critical difference at 5% level probability, by Tukey’s honestly significant test.
Table 4. Agronomic performance of Pusa 44 qDTY12.1-NILs under stressed and unstressed conditions during Kharif 2018 in New Delhi.
Table 4. Agronomic performance of Pusa 44 qDTY12.1-NILs under stressed and unstressed conditions during Kharif 2018 in New Delhi.
LineDFF (Days)Ht (cm)CmNPnL (cm)Yld (g.sqm−1)
SUSSUSSUSSUSSUSrPLY
Pusa 3003-15-121-30-1491.5 ab91.5 b69.177.6 ab8.9 b9.6 b22.321.9 ab295.2 a716.9 ab2.4
Pusa 3003-15-121-30-2391.0 ab91.5 b69.574.0 ab7.9 b8.1 b22.021.3 b327.4 a707.5 abc2.2
Pusa 3003-15-121-31-195.0 ab95.5 ab66.862.3 c9.0 b9.3 b21.420.4 b288.9 a730.1 a2.5
Pusa 3003-15-121-29-2383.0 b86.5 b77.378.1 a7.9 b8.8 b20.322.4 ab198.7 ab708.8 abc3.6
Pusa 3003-15-121-9-192.0 ab93.0 ab67.573.4 abc7.4 b7.5 b22.820.1 b297.7 a601.9 abc2.0
Pusa 3003-15-121-9-491.5 ab95.5 ab65.471.0 abc6.4 b8.1 b20.620.4 b301.7 a576.0 c1.9
Pusa 3003-15-121-9-2792.0 ab92.0 ab72.376.4 ab8.8 b9.4 b22.122.0 ab249.8 a587.7 bc2.4
Pusa 4496.0 ab106.0 a73.075.9 ab7.3 b7.9 b22.120.9 b75.2 b732.9 a9.8
IR90019-22-28-2-B102.0 a103.5 a65.167.6 bc16.7 a17.0 a22.527.4 a290.6 a437.2 d1.5
CD (5%)12.812.0ns11.76.76.9ns5.7148.3136.4
CV%3.43.15.75.118.517.74.56.414.15.2
DFF, days to 50% flowering; Ht, plant height; CmN, tiller (culm) number; PnL, panicle length; Yld, plot yield; rPLY, the ratio between unstressed and stressed yield (reciprocal of yield stability index); ns, non-significant; CV%, coefficient of variation; CD, critical difference; S, stressed; US, unstressed. Means followed by similar letters are not statistically significant at 5% level, based on Tukey’s honestly significant test.
Table 5. Agronomic characteristics of Pusa 44 qDTY12.1-NILs under stressed and unstressed conditions during Kharif 2019 in New Delhi.
Table 5. Agronomic characteristics of Pusa 44 qDTY12.1-NILs under stressed and unstressed conditions during Kharif 2019 in New Delhi.
PedigreeDFFPHTTLNPNLPLYRPGR (%)
SUSSUSSUSSUSSUSrPLY
Pusa 3003-15-121-17-47-282.0 d92.5 b82.7 a84.4 b–e18.7 a10.821.6 ab20.8 d295.1 a799.1 ab2.797.4
Pusa 3003-15-121-17-47-589.5 bcd93.0 b73.4 abc88.4 a–d13.8 ab13.621.2 ab21.8 bcd176.7 de847.6 ab4.897.4
Pusa 3003-15-121-29-2-193.0 bcd93.0 b61.0 de80.7 def14.2 ab14.321.2 ab23.6 bc141.8 efg589.5 bc4.298.7
Pusa 3003-15-121-29-21-495.0 bc96.5 ab61.5 de80.9 def14.6 ab13.122.3 ab23.7 b121.7 fgh743.2 ab6.198.7
Pusa 3003-15-121-29-23-291.5 bcd94.0 b63.9 cd83.2 cde13.6 ab15.321.4 ab22.0 bcd158.9 ef707.6 abc4.598.7
Pusa 3003-15-121-29-25-289.0 cd93.5 b72.2 bc88.2 a–d12.6 ab11.621.6 ab22.2 bcd178.2 de774.0 ab4.398.7
Pusa 3003-15-121-29-3-188.5 cd93.0 b64.2 cd87.2 b–e13.0 ab11.921.6 ab24.1 b106.3 gh862.2 ab8.198.7
Pusa 3003-15-121-29-32-392.0 bcd91.5 b67.8 bcd96.1 a10.2 b15.319.8 ab23.5 bc107.5 gh752.7 ab7.097.4
Pusa 3003-15-121-31-3290.5 bcd91.0 b65.4 bcd84.8 b–e12.6 ab13.619.6 b22.9 bcd124.5 fgh910.3 a7.398.7
Pusa 3003-15-121-9-37-295.0 bc93.5 b68.6 bcd87.2 b–e11.9 b15.220.9 ab21.9 bcd140.8 efg824.5 ab5.997.4
Pusa 3003-15-121-9-43-293.0 bcd92.5 b65.4 bcd84.9 b–e12.8 ab12.020.8 ab22.2 bcd114.0 fgh776.5 ab6.894.7
Pusa 3003-15-121-9-5-391.5 bcd92.0 b65.5 bcd84.0 b–e13.5 ab14.721.1 ab23.3 bcd95.0 h734.7 abc7.797.4
Pusa 3003-15-121-9-5-598.0 bc93.0 b67.1 bcd86.4 b–e13.8 ab14.422.3 ab21.9 bcd105.7 gh694.6 abc6.697.4
Pusa 3003-15-121-9-6-391.5 bcd92.5 b74.5 ab92.2 ab13.6 ab12.022.0 ab24.2 b207.7 cd843.7 ab4.197.4
Pusa 3003-15-121-9-7-295.0 bc93.0 b64.9 bcd85.8 b–e13.3 ab12.921.5 ab23.9 b178.2 de775.2 ab4.494.8
Pusa 3003-15-121-30-2397.5 bc92.5 b67.4 bcd86.6 b–e14.7 ab11.921.4 ab22.9 bcd235.7 bc771.5 ab3.398.7
Pusa 3003-15-121-31-1100.5 bc94.0 b53.7 e72.9 f14.8 ab16.619.6 b21.1 cd283.6 a592.5 bc2.197.4
Pusa 3003-15-121-9-2798.0 b97.5 ab66.5 bcd80.4 def14.6 ab17.321.6 ab22.7 bcd267.3 ab714.2 abc2.798.7
Pusa 44116.5 a105.0 ab69.3 bcd79.7 ef14.8 ab14.923.4 a23.2 bcd93.3 h743.1 ab8.0-
IR90019-22-28-2-B102.0 b103.5 a65.1 bcd90.8 abc16.7 ab17.022.5 ab27.4 a315.7 a473.5 c1.5-
CV%3.42.13.62.411.713.24.22.76.510.0
CD (5%)12.98.09.88.46.6ns3.62.645.1301.8
DFF, days to 50% flowering; PHT, plant height; TLN, tillers number hill; PNL, panicle length; PLY, plot yield; rPLY, ratio between unstressed and stressed yield (reciprocal of yield stability index); ns, non-significant; CV%, coefficient of variation, CD, critical difference; S, stressed; US, unstressed. Means followed by similar letters are not statistically significant at 5% level, based on Tukey’s honestly significant test.
Table 6. Drought tolerance indices for Pusa 44 qDTY12.1 near isogenic lines based on the performance in Kharif 2019.
Table 6. Drought tolerance indices for Pusa 44 qDTY12.1 near isogenic lines based on the performance in Kharif 2019.
PedigreeDrought Response IndicesAverage Rank
RTOLGMPSSISTIHAMYIYSI
Pusa 3003-15-121-17-47-262.8 d485.3 a0.81 cd0.43 a430.6 a1.73 a0.37 bc2.0
Pusa 3003-15-121-17-47-579.0 ab386.3 abc1.02 ab0.27 b–e291.9 cd1.03 de0.21 cde8.7
Pusa 3003-15-121-29-2-175.9 abc289.0 cde0.98 abc0.15 de228.5 d–g0.83 efg0.24 cde10.9
Pusa 3003-15-121-29-21-483.6 ab300.7 cde1.08 ab0.16 de209.1 efg0.71 fgh0.16 de13.3
Pusa 3003-15-121-29-23-277.6 abc335.3 b–e1.00 abc0.20 b–e259.5 de0.92 ef0.22 cde10.3
Pusa 3003-15-121-29-25-276.9 abc371.3 bcd1.00 abc0.25 b–e289.6 cd1.04 de0.23 cde8.0
Pusa 3003-15-121-29-3-187.7 a302.6 cde1.14 a0.16 cde189.2 efg0.62 gh0.12 e16.7
Pusa 3003-15-121-29-32-385.8 a284.3 de1.12 a0.15 de188.0 efg0.63 gh0.14 e16.1
Pusa 3003-15-121-31-3286.2 a336.3 b–e1.12 a0.20 b–e218.8 d–g0.73 fgh0.13 e13.9
Pusa 3003-15-121-9-37-282.9 ab340.6 b–e1.08 ab0.21 b–e240.4 def0.82 efg0.17 de11.3
Pusa 3003-15-121-9-43-285.2 a297.0 cde1.10 ab0.16 de198.5 efg0.67 fgh0.105 de15.3
Pusa 3003-15-121-9-5-387.1 a264.2 e1.13 a0.12 e168.3 fg0.56 h0.13 e18.3
Pusa 3003-15-121-9-5-584.8 ab271.0 de1.10 ab0.13 e183.5 fg0.61 gh0.15 de16.3
Pusa 3003-15-121-9-6-375.2 abc418.3 ab0.97 abc0.31 abc333.0 bc1.21 cd0.25 cde5.4
Pusa 3003-15-121-9-7-277.0 abc371.6 bcd0.99 abc0.25 b–e289.7 cd1.04 de0.23 cde7.7
Pusa 3003-15-121-30-2369.4 bc426.1 ab0.90 bc0.32 ab360.7 abc1.38 bc0.31 cd4.3
Pusa 3003-15-121-31-150.1 d407.4 ab0.65 d0.30 a–d379.5 ab1.66 a0.50 b3.1
Pusa 3003-15-121-9-2762.6 cd436.8 ab0.81 cd0.34 ab388.9 ab1.56 ab0.37 bc2.7
Pusa 4487.5 a263.2 e1.13 a0.12 e165.7 g0.54 h0.12 e19.1
IR90019-22-28-2-B33.6 e356.4 b–e0.43 e0.23 b–e349.1 bc1.70 a0.66 a4.0
CV%5.17.25.216.76.76.516.1-
CD (5%)15.8101.80.20.1573.30.260.16-
RTOL, relative stress tolerance; GMP, geometric mean productivity; SSI, stress susceptibility index; STI, stress tolerance index; HAM, harmonic mean; YI, yield index; YSI, yield susceptibility index. Rank sum is the sum of individual ranks for different tolerance index for each line. Average rank is the mean of ranks of all the stress response parameters. Correlations between average rank and grain yield under stressed and unstressed conditions were −0.958 and 0.262, respectively. Means followed by similar letters are not statistically significant at 5% level, based on Tukey’s honestly significant test.
Table 7. Grain quality data for Pusa 44 qDTY12.1 NILs under stressed and unstressed treatments.
Table 7. Grain quality data for Pusa 44 qDTY12.1 NILs under stressed and unstressed treatments.
PedigreeKLBCKBBCLBRKLACER
SUSSUSSUSSUSSUS
Pusa 3003-15-121-17-47-25.95 b6.88 a–e2.09 ab2.052.85 ab3.3610.3 b10.5 c–g1.73 abc1.52 c
Pusa 3003-15-121-17-47-55.59 cde7.15 abc1.96 b1.932.84 ab3.709.4 def11.3 a–e1.68 bc1.58 abc
Pusa 3003-15-121-29-2-15.23 fg6.56 de2.01 ab2.062.59 b–e3.219.3 ef10.4 d–g1.75 abc1.58 abc
Pusa 3003-15-121-29-21-45.05 g6.62 cde2.01 ab2.122.52 de3.149.1 ef10.5 c–g1.81 ab1.59 abc
Pusa 3003-15-121-29-23-25.41 def6.64 b–e2.04 ab2.182.64 b–e3.069.3 def10.3 efg1.72 abc1.55 bc
Pusa 3003-15-121-29-25-25.30 efg6.67 b–e1.98 b2.022.67 b–e3.318.9 f10.8 b–g1.67 bc1.62 abc
Pusa 3003-15-121-29-3-15.12 fg6.97 a–e1.97 b2.062.60 b–e3.398.9 f11.5 abc1.73 abc1.64 abc
Pusa 3003-15-121-29-32-35.16 fg7.26 a2.04 ab2.062.52 de3.519.2 def11.7 ab1.79 abc1.61 abc
Pusa 3003-15-121-31-325.22 fg6.91 a–e2.04 ab2.102.56 cde3.299.5 cde12.1 a1.83 a1.75 a
Pusa 3003-15-121-9-37-25.28 efg6.75 a–e2.01 ab2.022.62 b–e3.369.4 def11.0 b–g1.77 abc1.62 abc
Pusa 3003-15-121-9-43-25.34 d–g6.67 b–e1.96 b2.012.72 b–e3.329.5 cde10.2 fg1.77 abc1.53 c
Pusa 3003-15-121-9-5-35.42 def7.06 a–d2.01 ab2.042.70 b–e3.489.5 cde12.3 a1.75 abc1.74 ab
Pusa 3003-15-121-9-5-55.43 def6.48 e1.98 b1.922.73 b–e3.389.7 cd10.1 g1.79 ab1.56 abc
Pusa 3003-15-121-9-6-35.42 def6.81 a–e2.00 ab2.042.71 b–e3.359.8 bcd10.8 b–g1.81 ab1.58 abc
Pusa 3003-15-121-9-7-25.44 def7.17 ab2.00 ab2.102.71 b–e3.429.5 cde11.7 ab1.75 abc1.64 abc
Pusa 3003-15-121-30-235.65 bcd6.89 a–e2.04 ab2.052.76 bcd3.3710.0 bc11.4 a–d1.77 abc1.66 abc
Pusa 3003-15-121-31-15.59 cde6.83 a–e2.12 ab2.132.64 b–e3.219.6 cde11.3 a–f1.71 abc1.65 abc
Pusa 3003-15-121-9-275.81 bc6.67 b–e2.08 ab1.932.79 bc3.479.6 cde10.9 b–g1.65 c1.63 abc
Pusa 445.06 g6.83 a–e2.04 ab2.032.49 e3.379.2 def10.8 b–g1.83 a1.58 abc
IR90019-22-28-2B4.93 h5.44 f1.62 c1.913.04 a3.2611.1 a10.1 g1.70 abc1.54 bc
CV%1.571.981.894.042.364.191.442.371.913.04
CD (5%)0.340.540.15ns0.25ns0.551.050.130.19
KLBC, kernel length before cooking, in mm; KBBC, kernel breadth before cooking, in mm; L/B, length/breadth ratio; KLAC, kernel length after cooking, in mm; KBAC, kernel breadth after cooking, in mm; ER, elongation ratio; HULL, hulling recovery percent; MILL, milling recovery percent; ns, non-significant; CV%, coefficient of variation, CD, critical difference; S, stressed; US, unstressed. Means followed by similar letters are not statistically significant at 5% level, based on Tukey’s honestly significant 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

Oo, K.S.; Krishnan, S.G.; Vinod, K.K.; Dhawan, G.; Dwivedi, P.; Kumar, P.; Bhowmick, P.K.; Pal, M.; Chinnuswamy, V.; Nagarajan, M.; et al. Molecular Breeding for Improving Productivity of Oryza sativa L. cv. Pusa 44 under Reproductive Stage Drought Stress through Introgression of a Major QTL, qDTY12.1. Genes 2021, 12, 967. https://doi.org/10.3390/genes12070967

AMA Style

Oo KS, Krishnan SG, Vinod KK, Dhawan G, Dwivedi P, Kumar P, Bhowmick PK, Pal M, Chinnuswamy V, Nagarajan M, et al. Molecular Breeding for Improving Productivity of Oryza sativa L. cv. Pusa 44 under Reproductive Stage Drought Stress through Introgression of a Major QTL, qDTY12.1. Genes. 2021; 12(7):967. https://doi.org/10.3390/genes12070967

Chicago/Turabian Style

Oo, Kyaw Swar, Subbaiyan Gopala Krishnan, Kunnummal Kurungara Vinod, Gaurav Dhawan, Priyanka Dwivedi, Pankaj Kumar, Prolay Kumar Bhowmick, Madan Pal, Viswanathan Chinnuswamy, Mariappan Nagarajan, and et al. 2021. "Molecular Breeding for Improving Productivity of Oryza sativa L. cv. Pusa 44 under Reproductive Stage Drought Stress through Introgression of a Major QTL, qDTY12.1" Genes 12, no. 7: 967. https://doi.org/10.3390/genes12070967

APA Style

Oo, K. S., Krishnan, S. G., Vinod, K. K., Dhawan, G., Dwivedi, P., Kumar, P., Bhowmick, P. K., Pal, M., Chinnuswamy, V., Nagarajan, M., Bollinedi, H., Ellur, R. K., & Singh, A. K. (2021). Molecular Breeding for Improving Productivity of Oryza sativa L. cv. Pusa 44 under Reproductive Stage Drought Stress through Introgression of a Major QTL, qDTY12.1. Genes, 12(7), 967. https://doi.org/10.3390/genes12070967

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