Next Article in Journal
Organic Matter Properties of Spent Button Mushroom Substrate in the Context of Soil Organic Matter Reproduction
Next Article in Special Issue
Brassinosteroid Signaling Pathways Interplaying with Diverse Signaling Cues for Crop Enhancement
Previous Article in Journal
Unmanned Aerial Vehicles in Agriculture: A Survey
Previous Article in Special Issue
Combined Linkage Mapping and Genome-Wide Association Study Identified QTLs Associated with Grain Shape and Weight in Rice (Oryza sativa L.)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Introgression of qDTY1.1 Governing Reproductive Stage Drought Tolerance into an Elite Basmati Rice Variety “Pusa Basmati 1” through Marker Assisted Backcross Breeding

by
Gaurav Dhawan
1,2,
Aruna Kumar
2,
Priyanka Dwivedi
1,
Subbaiyan Gopala Krishnan
1,*,
Madan Pal
3,
Kunnummal Kurungara Vinod
1,
Mariappan Nagarajan
4,
Prolay Kumar Bhowmick
1,
Haritha Bollinedi
1,
Ranjith Kumar Ellur
1,
Kuram Tirumala Ravikiran
1,†,
Pankaj Kumar
1 and
Ashok Kumar Singh
1
1
Division of Genetics, ICAR-Indian Agricultural Research Institute, New Delhi 110012, India
2
Amity Institute of Biotechnology, Amity University, Noida 201303, India
3
Division of Plant Physiology, ICAR-Indian Agricultural Research Institute, New Delhi 110012, India
4
Rice Breeding and Genetics Research Centre, ICAR-Indian Agricultural Research Institute, Aduthurai 612101, India
*
Author to whom correspondence should be addressed.
Present address: Regional Research station, ICAR-Central Soil Salinity Research Institute, Lucknow 226002, India.
Agronomy 2021, 11(2), 202; https://doi.org/10.3390/agronomy11020202
Submission received: 18 December 2020 / Revised: 14 January 2021 / Accepted: 15 January 2021 / Published: 21 January 2021
(This article belongs to the Special Issue Rice Genetics: Trends and Challenges for the Future Crops Production)

Abstract

:
Drought stress at the reproductive stage in rice is one of the most important cause for yield reduction, affecting both productivity and quality. All Basmati rice varieties, including the popular cultivar “Pusa Basmati 1 (PB1)” is highly sensitive to reproductive stage drought stress (RSDS). We report for the first time, improvement of a Basmati cultivar for RSDS tolerance, with the introgression of a major quantitative trait locus (QTL), “qDTY1.1” into PB1. The QTL was sourced from an aus variety, Nagina 22 (N22). A microsatellite (simple sequence repeat (SSR)) marker “RM 431” located at telomeric end (38.89 mb) of chromosome 1, and located within a 1.04 mb QTL region was employed for foreground selection for qDTY1.1 in the marker assisted backcross breeding process. A set of 113 SSR markers polymorphic between N22 and PB1 were utilized for background selection to ensure higher recurrent parent genome recovery. After three backcrosses followed by five generations of selfing, eighteen near isogenic lines (NILs) were developed, through combinatory selection for agro-morphological, grain and cooking quality traits. The NILs were evaluated for three consecutive Kharif seasons, 2017, 2018 and 2019 under well-watered and drought stress conditions. RSDS tolerance and yield stability indicated that P1882-12-111-3, P1882-12-111-5, P1882-12-111-6, P1882-12-111-7, P1882-12-111-12, P1882-12-111-15 and P1882-12-111-17 were best in terms of overall agronomic and grain quality under RSDS. Additionally, NILs exhibited high yield potential under normal condition as well. The RSDS tolerant Basmati NILs with high resilience to water stress, is a valuable resource for sustaining Basmati rice production under water limiting production environments.

1. Introduction

India is gifted with a vast rice varietal diversity spread across diverse ecosystems. The region spanning from the Himalayan foothills traversing through north-western Indo-Gangetic plains is particularly bestowed with rices of incredible quality, popularly known as Basmati [1]. Over the time, Basmati has become pride possession of Indian subcontinent, serving exquisite cuisine and thereby an export commodity for trade in the world market. Basmati export from India alone, has earned foreign exchange worth of INR 310,255.6 million (US$ 4330.6 million) during 2019–2020 from an annual milled rice production of more than eight million tons [2]. Besides this, almost an equal amount of as much as exported, is consumed in the domestic market. Since Basmati rice is confined to a specific geographic boundary earmarked as Basmati geographical indication (GI) area, their natural genetic diversity remains low rendering them sensitive to all major biotic and abiotic stress factors. With the ensuing threat of global climate change, rice is poised to suffer the onslaught of several stresses, particularly abiotic stresses. Among these, drought would be a major factor as global warming is recognized as one of the important causes of climate change [3]. Global warming results from atmospheric increase in greenhouse gases propelling vagaries of temperature and rainfall [4]. Furthermore, climate change envisages exacerbation of drought in the coming years [5] and the strength and frequency of drought are expected to become worse [6]. Rice requires about 3000 to 5000 L of water to yield one kilogram of rice [7], and erratic rainfall can drive rice production ambivalent in rainfed regions. Estimates indicate that almost 34 million ha of rice in rainfed lowland areas in Asia alone recurrently suffers from water stress, in addition to an area of eight million ha under upland ecology [8].
Rice is grown throughout the year in India depending on the geographical region, in different seasons such as Kharif and Rabi. Kharif season coincides sowing in late summer, and is the principal rice growing season. Rabi which coincides with the winter sowing, is however limited to peninsular India, where winter is mild. Rabi, however, forms the second most important rice growing season with plantings restricted to southern and eastern India. Depending on geography, the weather parameters of Kharif and Rabi season vary, including intermittent spells of stress such as drought. Onset of drought can happen in all stages of rice crop. However, sensitivity to the drought stress depends on its duration and intensity. Although rice is affected badly in all phenological stages of growth, the stress occurring during the reproductive stage is particularly perilous, leading to significant reduction in grain yield [9,10]. The impact on rice plant is multifaceted, ranging from the reduction in dry matter accumulation to poor apportioning of metabolites from the source (stem and leaves) to sink (grain) [11]. This results in reduced number of filled grains per panicle, reduced grain weight and, ultimately, the grain yield. Varietal diversity in rice embodies various mechanisms for drought tolerance particularly at the reproductive stage. Tolerance is established through manifestations such as reduction in the tiller number, leaf area, shortening and thickening of leaves, leaf rolling and promotion leaf senescence [12,13]. These mechanisms offer opportunities for varietal improvement targeting drought prone areas. Additionally, breeding of futuristic cultivars resilient to multiple stresses requires enhanced drought tolerance, especially at the reproductive stage. Drought stress at reproductive stage leads to delay in flowering which slows down the panicle development, impairs flowering and grain development leading to crop loss, that can be total at severe situations. Stress occurring at anthesis stage leads to sterile panicle, because the panicle that is still inside the leaf sheath, fails to develop. Failure of reproduction ultimately translates into significant yield decline. The reproductive stage drought stress (RSDS) tolerance in rice is considered as a complex and multifarious trait, governed by several genes, major and minor. These may include a large number of quantitative trait loci (QTLs) that may comprise of structural and functional genes as well as regulatory elements such as transcription factors, controlling a multitude of morpho-physiological and biochemical responses [14].
Exploiting molecular markers of various kind, several studies have reported QTLs governing drought tolerance in rice, primarily using grain yield under drought stress (DTY) as the surrogate trait. Among these, qDTY1.1 is a major QTL for grain yield under drought stress (explaining about 12.6% to 16.9% of phenotypic variance), mapped on chromosome 1 independently from two tolerant cultivars such as Nagina 22 (N22) [15] and Dhagaddeshi [16]. qDTY1.1 was the first QTL reported to have consistent effect in different genetic backgrounds and was shown to affect many drought related traits like root length, relative water content, root thickness, osmotic adjustment and biomass [16,17]. Previous studies have demonstrated that qDTY1.1 was flanked by markers RM431 on telomeric end and RM11943 on the centromeric end. Additionally, Bernier et al. [17] identified a different major QTL on chromosome 12, qDTY12.1, from the cross between Vandana and Way Rarem explaining approximately 51% of phenotypic variation. There are also other major QTLs reported which were demonstrated effective either under upland or lowland situations [14,18].
Marker assisted introgression of major-effect QTLs could be a proficient and rapid approach for breeding rice varieties tolerance to drought stress [17]. Consequent attempts to introgress/pyramid these QTLs by marker assisted selection, primarily into mega-varieties has found significant advancements towards breeding climate-adaptive cultivars [19,20], such as Sabitri (qDTY3.2 and qDTY12.1) [21], IR 64 (qDTY2.2 and qDTY4.1), Vandana (qDTY12.1) [22] and Pusa 44 [23]. In the past twelve years, there has been about sixty-six RSDS tolerant varieties released around the world, which involve several of these QTLs [14]. Although important, there has not been a previous attempt to improve Basmati cultivars for RSDS tolerance. One of the major reasons for this hiatus was the lack of donors from the Basmati group. Use of donors from the non-Basmati backgrounds for Basmati improvement, present a major challenge of loss of grain quality of the Basmati parent, while introgression [24]. However, marker assisted backcross breeding coupled with phenotypic selection, has been demonstrated to effectively address this problem [25,26,27].
Developed by the Indian Council of Agricultural Research (ICAR)-Indian Agricultural Research Institute (ICAR-IARI), New Delhi and released for commercial cultivation during year 1989, Pusa Basmati 1 (PB1) is the world’s first semi-dwarf (105 cm) Basmati variety with high yielding potential. This variety, that showcased a tremendous level of improvement over the traditional Basmati cultivars that were low yielding (~2.3 tons/ha), photosensitive, tall (>150 cm), lodging with weaker stem and of long duration (>150 days), became popular among the farmers in no time. Besides the high yield (5.0 tons/ha), PB1 possesses insensitivity to photoperiodism, resistance to lodging, superior grain quality, semi-dwarf stature and a duration of 135–140 days. The grain quality of PB1 was adorned with strong aroma with explicit cooking qualities such as high kernel elongation ratio of 1.8, an average milled grain length of 7.4 mm and cooked kernel length of 13.7 mm [19]. However, as that of traditional cultivars, PB1 too is highly sensitive to several stress factors such as pests and diseases as well as drought. This study forms the maiden attempt to improve any Basmati cultivar towards climate resilience, by introgressing RSDS tolerance into PB1. The augmented objectives were to recover the Basmati grain quality traits as well as the recovery of all the agronomic traits including yield in the near isogenic lines (NILs). Further, evaluation of the improved NILs across multiple locations and environments to ascertain the stability in yield and drought tolerance to deploy as commercial Basmati rice cultivar.

2. Materials and Methods

2.1. Plant Material and Experimental Sites

PB1, the elite Basmati cultivar was selected as the recurrent parent. An aus cultivar, N22, was used as donor parent for the QTL, qDTY1.1. N22 is a pure line selection from a landrace, Rajbhog. Well-known as a universal donor for high temperature tolerance, drought tolerance and grain dormancy [28], N22 is a tall (~120 cm), short duration (90–95 days) cultivar, with short bold non-aromatic grains. The entire study was conducted between two locations, New Delhi and Aduthurai. Experiments at New Delhi was conducted in the research fields of the Genetics Division, ICAR-IARI, situated at 28°35′ N latitude, 77°12′ E longitude with an altitude of 228.16 m above mean sea level. The other experimental site at Aduthurai was the research farm of IARI-Rice Genetics and Breeding Research Centre situated at 11°00′ N latitude, 79°28′ E longitude and 20 m above MSL. The whole generation advancement between New Delhi and Aduthurai was carried out through shuttle breeding approach. Final drought stress evaluation of NILs was done at New Delhi for three subsequent crop seasons.

2.2. Development of NILs through Marker Assisted Backcrossing

Hybridization was taken up between PB1 and N22 at New Delhi, using PB1 as recipient (female) and N22 as donor (male) in Kharif 2013. Hybrid seeds (F1s) were harvested and propagated in the late Rabi season of 2013–2014 at Aduthurai, Tamil Nadu. F1s were backcrossed to PB1 as recurrent parent (RP), to generate BC1F1. During subsequent seasons, two successive backcrosses were done to generate BC2F1 and BC3F1 during 2014 and 2014–2015. In the BC3 generation, from Kharif 2015, four successive selfings were carried out to generate BC3F5 NILs by 2016–2017 late Rabi season at Aduthurai. All through during different generations, selections were carried out to recover maximum recurrent parent genome recovery (RPGR) among the progenies. The selected NILs, beginning from the BC3F4 generation, were subjected to three consecutive drought screenings at New Delhi during Kharif seasons 2017 to 2019.
For molecular marker analyses, the genomic DNA was isolated from fresh leaf samples using CTAB method with minor modifications [29]. Polymerase chain reaction was performed with the volume of 10 μL containing 20–30 ng template DNA, 5 pmol of each primer, 0.05 mM dNTPs (MBI, Fermentas, Lithuania, USA), 10× PCR buffer (10 mM Tris, pH 8.4, 50 mM KCl and 1.8 mM MgCl2) and 0.5 U of Taq DNA polymerase (Bangalore Genei Pvt. Ltd., Bangalore, India), running a amplification profile consisting of one cycle of initial denaturation at 94 °C for 5 min; followed by 35 cycles containing denaturation with at 94 °C for 30 s, annealing at 55 °C for 30 s, extension at 72 °C for 1 min; and a final extension at 72 °C for 7 min. The amplified PCR products were resolved by electrophoresis having 3.5% agarose gel and fluoro-stained with ethidium bromide. The amplicon resolution was pictographed using a Gel Doc XR+® (BioRad Laboratories, Hercules, CA, USA) gel documentation system.
For the marker analyses, marker sequences were downloaded from the marker database at Gramene (http://www.gramene.org) and the primers were custom synthesized from Sigma-Aldrich (Bengaluru, India). The qDTY1.1 linked simple sequence repeat (SSR) marker RM431 (F: tcctgcgaactgaagagttg; R: agagcaaaaccctggttcac) was polymorphic between the parents and was thus used for foreground selection [15]. For background analysis, the genome wide polymorphism between PB1 and N22 was determined using a genome wide survey employing 651 SSR markers, which resulted in identification of 113 markers that could distinguish both the genotypes. These polymorphic SSR markers spanned across the rice genome and were used for background selection in various backcross generations (Supplementary Table S1), using the reductive screening approach [30]. The recovery of PB1 type alleles among the backcross progenies during each generation was used for computing the RPGR [31]. For the graphical comparison of the genomes of NILs and parents, graphical genotypes were drawn using GGT v.2.0 [32].

2.3. Field Management under Stressed and Unstressed Conditions

In order to assess the efficiency of qDTY1.1 introgression into PB1, the selected NILs along with their parents and checks were evaluated for three consecutive Kharif seasons of 2017 to 2019 period at ICAR-IARI, New Delhi, under an irrigated ecology with drought stress imposed at reproductive stage. The drought tolerant check was IR86918-B-B-305, a backcross inbred line carrying qDTY1.1 in the background of IR64 developed from the cross IR64/N22. The sensitive check used in this experiment was IR64. Two treatments were maintained, stressed and unstressed. Experiments were laid out using a randomized complete block design (RCBD) with two replications for each stress treatment. The plot size was 6.5 m2 with a spacing of 20 by 15 cm. Approximately 25 g of seeds were nursery sown to raise the population. In the plots designated for stress treatment, tensiometers were installed at every six NILs to monitor and characterize the soil moisture status. Initially, for both the treatments, seedlings were raised in a wet-bed nursery and after 21 days transplanted into a flooded field with 5 cm standing water. To ensure uniform establishment of the transplanted seedlings, all the plants were maintained under irrigation for 30 days post transplanting. At the 31st day after transplanting, water from the stress treatment plots were drained to initiate the stress. The stressed plots were left un-irrigated until the soil moisture tension reached −70 kPa at 30 cm depth. Severe leaf rolling and leaf drying were observed at this soil moisture level. At this severe stress, a flash life-saving irrigation was provided, and the excess water was drained out approximately after 24 h. This cycle was constantly repeated until harvest [15,33]. The unstressed plots were maintained with normal irrigation and the plants were maintained in the standing water. Altogether, three irrigations were given into stressed plots, while the unstressed plots were irrigated six times. During the crop duration, from June to November, a total of 839.8, 913.4 and 608.1 mm of rainfall was received in the years 2017, 2018 and 2019, respectively.

2.4. Phenotypic Data Collection

From each of the NILs under both treatments, phenotypic data were recorded from five randomly tagged plants on days to 50% flowering (DF)—the number of days from sowing to flowering in 50% plants/tillers was recorded, plant height at maturity (PH)—was measured from ground level to the tip of main panicle, panicle length (PL)—average length of the primary panicles was taken from peduncle base to the tip, number of reproductive tillers (NT)—average number of tillers containing grain filled panicles per plant, spikelet fertility % (SF)—the proportion of filled grains to the total number grains per panicle and grain yield (GY)—the harvested grain from the plant were dried to optimum moisture level of approximately 14% after which they were weighed. The NILs were also characterized for kernel length and breadth before and after cooking. From these measurements, length/breadth ratio (LBR) and cooking quality characteristics, such as kernel length elongation ratio (KER) were computed. Further, alkali spreading value (AS) and aroma were determined using standard protocols [34].

2.5. Statistical Analyses

Initially, independent analyses of variance (ANOVA) was carried for each season to identify significant responses among the NILs under stressed and unstressed situations and to compare them with the checks. A subsequent combined ANOVA was carried out using linear mixed model approach with genotypes as fixed factor and seasons and stress as random factors. The trait predictions obtained from the model was saved as best linear unbiased predictors (BLUPs). The data were analyzed using STAR package version 2.0.1. (http://bbi.irri.org/products). To identify the genotypes with stable performance across the environments, seasons under drought condition, an additive main effective and multiplicative interaction (AMMI) model was constructed using BLUPs and AMMI stability value (AST) [35] and yield stability index (YSI) [36] were generated, as follows:
AMMI   stability   value ,   AST = SS IPC 1 SS IPC 2 × IPC 1 2 + ( IPC 2 ) 2
Yield   stability   index ,   YSI   =   R AST +   R Y
where, SSIPC1 and SSIPC2 are the sum of squares of interaction principal component axes (IPCA) 1 and 2, respectively, and IPC1 and IPC2 are the respective IPCA scores. Similarly, RAST and RY are the respective genotypes ranks based on AST and yield (Y).
To compare the effect of drought on different traits, a forward stepwise regression analysis was carried out for delineating the traits which contributed significantly to yield. Based on the BLUPs for grain yield under stressed (S) and unstressed (NS) treatments, different indices namely, drought yield index [37], stress tolerance index [38], stress susceptibility index [39], and percent reduction in yield were worked out as follows:
Drought   yield   index ,   DYI   =   ( Y i ) NS / ( Y i ) S ( Y G ) NS / ( Y G ) S
Stress   tolerance   index ,   STI   =   ( Y i ) NS × ( Y i ) S ( Y NS ) 2
Stress   susceptibility   index ,   SSI   =   ( Y i ) NS ( Y i ) S Y NS Y S × Y NS ( Y i ) NS
Percent   reduction   in   yield ,   % R   =   ( Y i ) NS ( Y i ) S ( Y i ) NS ×   100
where, Yi represents the mean yield of a genotype i on the untransformed scale, YG refers to geometric mean across genotypes, the suffices, NS and S represent unstressed and stressed conditions, respectively and Y refers to the arithmetic mean yield across genotypes.

3. Results

3.1. Parental Polymorphism

The genome wide polymorphism survey using 651 SSR markers, revealed a diversity of 17.4% between the parents PB1 and N22, identifying 113 markers to be polymorphic between them (Table 1). Out of 101 markers tested on chromosome 1, the carrier chromosome of qDTY1.1, 14 markers were observed polymorphic including the linked marker, RM431. The target chromosome diversity was 13.9%. Of the remaining chromosomes, highest diversity was found for chromosome 8 (40.9%) while the chromosome 6 (9.2%) indicated low diversity.

3.2. Introgression of qDTY1.1 into PB 1

The breeding scheme for development of PB1 NILs carrying qDTY1.1 is given in Figure 1. From a total of 23 seeds initially collected of the cross PB1/ N22, five plants were found to be pure hybrids showing heterozygosity for the foreground marker, RM431 linked to qDTY1.1. These F1s were backcrossed to the recurrent parent, PB1 to obtain 39 BC1F1 seeds. Out of these, 15 seedlings were confirmed heterozygous for qDTY1.1 by foreground selection. These 15 plants were further analyzed with 113 polymorphic background SSRs (Table 2). Based on the recovery of PB1 alleles, the RPGR among the BC1F1 plants was estimated to range between 78.0% and 86.7%. The plant with highest recovery (86.7%) had 83 background markers in homozygous state for recurrent parent (RP) allele, while 30 markers were heterozygous. Further, these plants were also tested for their agronomic and grain quality similarities with PB1. The plant showing phenotypic resemblance and high RPGR (86.7%) was used to backcross with PB1 to generate 37 BC2F1 seeds. However, only one BC2F1 plant was found to be heterozygous for RM431. Background analysis on this plant using 30 markers heterozygous from the previous backcross generation, showed the RPGR of 94.7%. Subsequent backcrossing of this plant to PB1 yielded 25 seeds. On raising the BC3F1 generation from these seeds, five plants were found to be heterozygous for qDTY1.1 linked marker, RM431. Further the one F1 out of five with maximum RPGR and excellent cooking quality was selected for further generation advancement.
Background selection with the remaining unrecovered heterozygous markers among the BC2F1 plant, resulted in an RPGR ranging from 95.6% to 97.4%. These BC3F1 plants were further subjected to phenotype matching with PB1, and the plant with maximum genome recovery (97.4%) as well as phenotype similarity was used for further step of selfing. Selfing the selected BC3F1, yielded 229 BC3F2 seeds. Through foreground followed by background selections of these BC3F2 plants, 55 plants homozygous for qDTY1.1 with an RPGR of 97.8–98.7% were identified. All these 55 plants were advanced to BC3F3 families which underwent a rigorous phenotype selection for agro-morphological and grain quality traits to select 24 families for further advancement to BC3F4 generation. The BC3F4 families were grown under both irrigated and stressed conditions. Subsequent selection among the 24 BC3F4 families was carried out for agronomic performance, drought response and grain quality. This resulted in 18 families of near isogenic lines, having close similarity to PB1 and comparable grain and cooking qualities. Phenotype matching of these 18 NILs with PB1, both in terms of agronomic features and grain quality, identified them closer to PB1. All the 18 NILs had target marker pattern of the donor (250 bp) parent for qDTY1.1 (Figure 2). The RPGR at BC3F5 generation, ranged between 98.2% and 99.1%. All along the selfing generation, selection was carried out with major emphasis on recovery of grain quality, agro-morphological features, and yield of PB1, among the NILs. The graphical genotype of targeted segment for qDTY1.1 from donor and recovery of background genome on chromosome 1 in 18 NILs are depicted in Figure 3.

3.3. Recovery of RP Alleles in Carrier Chromosome

The polymorphism survey between PB1 and N22 on chromosome 1, the carrier chromosome of qDTY1.1, was performed using 102 markers, of which 14 were found polymorphic including the foreground marker, RM431. There were two flanking markers for qDTY1.1 in the marker array, RM431 and RM11943. Among these, only RM431, the telomeric end marker was found polymorphic between the parents. RM11943 was found monomorphic. Further, among the downstream markers towards the telomeric end, only one marker RM6840 was polymorphic. However, RM6840 was at the distal end, 4.27 Mb away from RM431. On the centromeric end, however, the next upstream polymorphic marker was RM3825. RM3825 was used for recombinant selection. These markers, except for RM431, were included in the whole genome background recovery analysis. All the 14 polymorphic markers in chromosome 1 showed complete recovery of RP alleles among the 18 NILs by BC3F5 generation (Supplementary Figure S1).

3.4. Per se Performance of the NILs under Reproductive Stage Drought Stress and Unstress Treatments

Agronomic evaluation of NILs along with parents and checks over the three Kharif seasons under field-imposed drought stress as well as under normal conditions, indicated significant variation for several traits (Table 3). These evaluations belonged to BC3F4, BC3F5 and BC3F6 generations. ANOVA over individual years revealed significant variation for genotypes, treatment and genotype × treatment components, for most of the agronomic characters studied. However, only two traits, SF and GY showed consistently significant variation for genotype × treatment component in all the seasons. For the remaining traits, during 2018, genotype effects for NT and genotype × treatment effects for DF, NT and PL were found non-significant, as that of the genotype × treatment effect of NT in 2017. Combined ANOVA, indicated a uniform pattern across the traits, particularly for GY and SF. NT had non-significant effects for various sources of variation except genotype component. Traits that showed significant variation for all the sources except for one component, included DF and PH where in year × treatment interaction was non-significant, while PL had non-significant year × genotype effect. Except for PL and NT, year × treatment × genotype interaction was significant for all the other traits. Lease significant difference (LSD) values generated from the individual ANOVA in different years were utilized for mean comparison of NILs with checks.
Considering the individual traits (Supplementary Table S2), DF exhibited uniform pattern across three experimental years, but a general delay was observed in flowering under reproductive stage drought stress. The delay was conspicuous in IR64 and was particularly apparent during 2019 period. Among the NILs, P1882-12-111-9 and P1882-12-111-10 showed significant delay in DF under stress during year 2017 than rest of the NILs and PB1. In the same year, there were five NILs that showed significant late flowering than PB1. During 2018, both under stressed as well unstressed conditions, DF of all the NILs were at par with PB1. Similarly, during 2019 season too, DF of all the NILs were similar to PB1, except for P1882-12-111-1 under stress, that showed delayed flowering. On an average, the flowering delay under stress was between one to five days among the NILs. Overall, it was observed that drought stress increased the vegetative period in most of the NILs and checks across all the three years of evaluation (Table 4).
Significant reduction in PH was observed under drought stress in all the years. When stressed, among the NILs, PH ranged from 69 cm to 79 cm during 2017, 83 cm to 92 cm during 2018 and 74 cm to 86 cm during 2019. However, in the first two seasons, the PH among NILs were on par with PB1, but during 2019, six NILs showed significantly low PH. Whereas, under normal conditions all the NILs showed PH on par with PB1 during all the seasons, except for P1882-12-111-1 and P1882-12-111-9 during 2017. However, the average performance across the seasons revealed no significant difference among the NILs under both the treatments. No significant difference was noticed for PL and NT during 2017 and 2018 seasons under both unstressed and stressed situations. However, during 2019, differences were noticed in few of the NILs, for NT alone. In case of PL too, average values within each of the treatments showed no difference among the NILs. Reduction in SF percentage of all the NILs and checks was found under stress during all three years. Two NILs out of18, P1882-12-111-14 and P1882-12-111-15, had significantly better SF than PB1 under stress of 2017, which turned insignificant during 2018. Whereas during 2019, significantly better SF was found among thirteen NILs than PB1 under drought stress, which was on par for all the NILs and PB1 during unstressed treatment.
Most of the NILs out-yielded PB1 under stress in all the experimental years 2017, 2018 and 2019 under drought stress clearly reflecting the effect of the qDTY1.1 introgression (Supplementary Table S3). The GY of NILs under stress during 2017 ranged from 463 kg/ha (P1882-12-111-10) to 1015 kg/ha (P1882-12-111-20) in comparison to PB1 (309 kg/ha). In the same season, under unstressed treatment, the yield ranged from 5541 kg/ha (P1882-12-111-9) to 6355 kg/ha (P1882-12-111-5) in comparison to PB1 yield of 5581 kg/ha. When subjected to drought stress during 2018 season, the yield of NILs ranged from 325 kg/ha (P1882-12-111-1) to 1654 kg/ha (P1882-12-111-17) in comparison to PB1 (573 kg/ha), however, under unstressed field the yield was in the range of 4733 kg/ha (P1882-12-111-9) to 5955 kg/ha (P1882-12-111-3) as against PB1 (5291 kg/ha). During the next season, 2019, the yield of NILs was in the range of 727 kg/ha (P1882-12-111-1) to 2516 kg/ha (P1882-12-111-17) under stress as against the yield of 692 kg/ha for PB1. However, under the unstressed condition, the yield of all the NILs was at par with PB1. The GY of donor parent N22 was found 889 kg/ha 1796 kg/ha and 2632 kg/ha under stress during year 2017, 2018 and 2019, respectively.

3.5. Drought Tolerance Level of NILs Judged through Percent Reduction of Yield and Stress Indices

In order to assess the true tolerance of NILs, percent reduction in yield (%R) and stress indices were calculated. The average values are provided in Table 5. Further, the data for individual lines across different years is provided in Supplementary Table S3. The unstressed yield levels of the 18 NILs were non-significantly different from that of the recurrent parent, PB1 and the donor parent N22 in several cases. However, under drought, the yield level was significantly superior than PB1 in seventeen NILs, and seven NILs showed yield similar to N22. The %R was found to be maximum in PB1 in all the experimental years. The minimum and maximum %R among NILs was found to be 82% in P1882-12-111-20 and 93% in P1882-12-111-8 during year 2017. The donor parent, N22 maintained the least %R of 64% and minimum and maximum %R among NILs was for P1882-12-111-17 (72%) and P1882-12-111-1 (93%) in year, 2018. The %R was the lowest in 2019, among NILs as compared to previous years 2017 and 2018; it ranged from minimum 59% of P1882-12-111-17 to maximum 86% of P1882-12-111-1. The drought yield index (DYI) ranged from 0.98 (P1882-12-111-20) to 3.14 (PB1) in year 2017; DYI was between 0.49 (N22) to 2.79 (P1882-12-111-1) in year 2018 and in year 2019 DYI range was from 0.35 (N22) to 1.37 (PB1). Lowest DYI (0.7) was found for P1882-12-111-17 followed by six other NILs (0.9). The stress tolerance index (STI) was in the range 0.05 (PB1) to 0.19 (P1882-12-111-17) during 2017; during 2018 STI was ranged from 0.05 (P1882-12-111-1) to 0.31 (P1882-12-111-17) and during 2019 its range was from 0.12 (PB1) to 0.47 (P1882-12-111-17). P1882-12-111-17 had the highest STI in all the experimental seasons. SSI values are in range of 1.04 (P1882-12-111-17) to 1.18 (PB1) during 2017; while its range in 2018 was 0.90 (P1882-12-111-17) to 1.17 (P1882-12-111-1) and during 2019 it was 0.75 (P1882-12-111-17) to 1.09 (PB1). The SSI values have shown same pattern during 2018 and 2019, and on average P1882-12-111-17 showed the lowest SSI value of 0.90 followed by P1882-12-111-6 (0.94) and P1882-12-111-7 (0.95).

3.6. Stability Analysis of NILs under Stress Situation for Yield and Spikelet Fertility

To identify the NILs which could maintain superior performance under drought stress, an additive main effects and multiplicative interactions stability analysis across the season under drought condition, particularly for plot yield and spikelet fertility and stability indices, AMMI stability value and yield stability index was generated based on the AMMI model formulated (Table 5). For plot yield, the AST ranged from 2.80 (P1882-12-111-11) to 69.10 (IR86918-B-B-305) and the YSI range was from 12 (P1882-12-111-3, P1882-12-111-11, P1882-12-111-12, P1882-12-111-15) to 42 (P1882-12-111-1). Considering AST and YSI values with above average plot yields, P1882-12-111-11, P1882-12-111-8, P1882-12-111-15 and P1882-12-111-12 are the most stable among all NILs. For spikelet fertility, the AST range was from 0.61 (P1882-12-111-12) to 8.61 (P1882-12-111-4) and the YSI range was from 7 (P1882-12-111-15) to 40 (P1882-12-111-12). Considering AST, YSI and above average spikelet fertility values, NILs P1882-12-111-12, P1882-12-111-20 and P1882-12-111-13 were found most stable. Keeping in perspective the above two traits, it is apparent that the P1882-12-111-12 could maintain superior plot yield coupled with spikelet fertility across seasons under drought situation.

3.7. Stepwise Forward Regression Analysis to Identify the Significantly Affected Yield Contributing Trait

To find out the trait which has contributed significantly to yield and hence conditioned primarily by the QTL transferred, a stepwise forward regression was done with BLUP values for individual years and across years by taking the grain yield as the dependent variable (Table 6). SF stood out as the major determinant trait contributing to grain yield consistently across seasons, particularly in 2017, where it alone contributed for 51% of the variation in grain yield under drought conditions. It was also clear that the contribution of SF was more apparent under drought stress than unstressed conditions, from its significant positive coefficients across the seasons. Remaining yield contributing traits did not show strong and stable influence on yield.

3.8. Quality Assessment of NILs under Stressed and Unstressed Situations

The recovery grain and cooking quality traits is a major criterion in improvement of Basmati rice. Analysis of the grain quality recovered among the NILs under both stressed as well as unstressed conditions during 2019 season indicated that all the improved lines were as aromatic as PB1, the recurrent parent having a panel score of 2.0 for aroma (Table 7). Visually, the milled grains of the NILs appeared similar to that of PB1. However, some of the NILs showed non-significant variation for kernel length before cooking (KLBC) with respect to PB 1. By and large, similar pattern was observed for all the raw grain quality traits. After cooking, no significant variation was observed for kernel length after cooking (KLAC) than PB1 (Figure 4), among several of the NILs under unstressed conditions. However, under stress there were six NILs that showed significantly superior elongation than PB1, such as P1882-12-111-3, P1882-12-111-5, P1882-12-111-6, P1882-12-111-7, P1882-12-111-8 and P1882-12-111-11. Interestingly, under drought, the elongation ratio was found significantly good in all NILs except in seven NILs, that were on par with PB1. Furthermore, the alkali spreading value (ASV) of all the NILs were similar to PB1 (score of 7) as against the score of 5 recorded among the other checks.

4. Discussion

Climate change is predicted to affect rice cultivation worldwide adversely, through various impacts, such as drought, excess rainfall, temperature fluctuations and also predisposing it to several biotic stresses. Any adverse effect on rice production would threaten world food security, because more than half of the world population is rice dependent. This is particularly important to rice consuming countries like India, where about 90% of the total rice production is internally consumed [40]. Among these, biotic threats due to diseases and pests, as well as abiotic stresses such as drought, submergence and salinity are particularly relevant under climate change scenario. Predominantly in the rainfed environment of South and Southeast Asia, drought and submergence are more frequently encountered during the crop growing season [41]. Most of the prominent and popular rice varieties are vulnerable to these abiotic stresses [20]. Since climatic vagaries occur spontaneously and intermittently, growing tolerant rice varieties for wide range of stresses is the only economically viable option to manage abiotic stresses. This option is particularly important in the case of drought, because it remains as the most frequently occurring stress that bears the potential to fail a rice crop. Judging drought tolerance of genotypes and its transfer into elite backgrounds is quite laborious due to complex nature of tolerance. It is this complexity that renders the improvement for drought resilience through conventional breeding tardy. However, efforts can be remarkably hastened by leveraging molecular markers that are linked to various traits associated with drought tolerance. Moreover, use of markers can aid a relatively cost-effective and environment neutral selection, while improving the accuracy and reducing the turnover time. Molecular marker assisted introgression of traits provides various advantages such as easy recovery of an otherwise difficult phenotypic trait, easy selection, increased accuracy and shortening the breeding time. Therefore, marker assisted backcross breeding (MABB) was proven to be an efficient strategy for incorporation of desired trait associated genes/QTLs in numerous prominent rice varieties [21,42,43]. MABB along with stringent phenotypic selection was successfully employed for Basmati rice improvement for different biotic stresses [26,30,44,45,46,47,48]. MABB is also proven to a superior method for QTLs transfer into desired cultivar/variety for different abiotic stresses as well, such as drought. Out of the several QTLs reported for reproductive stage drought tolerance in rice, qDTY1.1 [15,16], qDTY2.1 [33], qDTY3.1 [33,49] and qDTY12.1 [50,51], showed consistent grain yield under drought across different genetic backgrounds and has been used in breeding applications. Recently, Dwivedi et al. [23] reported successful pyramiding of two QTLs, qDTY2.1 and qDTY3.1 into the megavariety, Pusa 44 significantly enhancing yield under drought. Despite the proven advantage, however, only few success stories in improving rice varieties for drought tolerance using MABB has been reported rice.
In the present study, MABB was utilized to introgress a RSDS tolerance QTL qDTY1.1 from N22 into PB1, one of the popular Basmati rice varieties of India. Popular for its high yield, better-quality grain, and excellent cooking quality and pleasing aroma, PB1 is the first ever semi-dwarf Basmati rice variety developed. Bred through conventional convergent breeding procedures, PB1 was derived from a cross between Pusa 150/Karnal local. Pusa 150 has Basmati 370 in its lineage, while Karnal local was a landrace with better aroma and cooking quality. Karnal local was a selection from Haryana Basmati collection 19 (HBC19), that was later released as Taraori Basmati [52]. Despite its excellent yield gain over the conventional Basmati cultivars, PB1 was semi-dwarf with excellent plant architecture, photo insensitivity, high yielding and with unparalleled grain quality [53]. However, like its congeners, PB1 is also sensitive to drought as well as to many biotic stresses. Released during 1989, PB1 remains popular even today among the Basmati farmers of India, despite above limitations. Furthermore, till date there is no report on improvement any Basmati rice variety for drought tolerance. Therefore, the present was carried out with the objective to incorporate drought tolerance in PB1. The donor parent, N22 is a tall upland rice variety belonging to aus group. N22 is a pure line selection from a landrace, Rajbhog [54] that possesses deeper root system, shorter duration and tolerant to heat and drought [55]. This variety has short bold non-aromatic grains with very less elongation on cooking.
The genome wide average polymorphism between the parents, PB1 and N22 was 17.35%, which was more than that of the carrier chromosome of qDTY1.1. The molecular analyses on the background recovery has revealed maximum recovery of chromosome 1, thereby eliminating the chances for linkage drag. The fact that maximum background markers (101 markers) were surveyed on chromosome 1, substantiate our claim for the elimination of linkage drag. Furthermore, the agronomic recovery of the NILs in general suggests that no undesirable trait has been incorporated by chance from other donor fragments, that might have escaped marker-based background selection. In this study, qDTY1.1 was found providing increased level of drought tolerance in the NILs vis-à-vis PB1. It is also known that QTLs from a particular genetic background usually show minor effects or may remain completely silent in diverse genetic backgrounds [56]. In the case of drought, interaction between QTL × genetic background has been a major bottleneck limiting the use of QTLs for MABB in rice [50,57,58]. Nevertheless, qDTY1.1 derived from N22 has been demonstrated effective in multiple non-Basmati backgrounds such as MTU1010, IR64 and Swarna, which showed its role as a major effective QTL governing grain yield under reproductive stage water stress [15]. This justifies the use of qDTY1.1 as a most desirable candidate for MABB to improve prominent high-yielding varieties with augmented drought tolerance.
Notwithstanding, drought has a significant role in reducing the yield under drought relative to irrigated cultivation. Under drought, yield reduction was found majorly influenced by spikelet fertility. In comparative assessment, there was no significant yield reduction noticed under unstressed condition among the NILs as a consequence of introgressing the qDTY1.1. In an earlier report on introgression of another QTL, qDTY3.1, a negative effect on grain yield was encountered by Venuprasad et al. [33] under unstressed treatment. The NILs evolved under this program, showed significant higher gain for grain yield over PB1, with P1882-12-111-2, P1882-12-111-3, P1882-12-111-5, P1882-12-111-6, P1882-12-111-7, P1882-12-111-10, P1882-12-111-12, P1882-12-111-15 and P1882-12-111-17 having clear advantage under both stressed and unstressed treatments. This also indicated that negative effect of QTL × background interaction was not observed under PB1 background. P1882-12-111-17 had the maximum grain yield among the NILs followed by P1882-12-111-3, P1882-12-111-5, P1882-12-111-6, P1882-12-111-7 and P1882-12-111-12 under both stressed, and unstressed conditions during all the years. This signified the effectiveness of qDTY1.1 in enhancing yield under drought stress. However, the effect of the introgressed QTL among the PB1 NILs was more conspicuous in the 2019 than in the previous years. In the recent report from World Meteorological Organization [59], it has been reported that climate change is accelerating by the last five-year period between 2015–2019. Although 2016, has been identified as the most erratic year, 2019 also showed significantly high incidence of drought than 2017 and 2018. Therefore, 2019 season data were more explicit to showcase the drought tolerance response of the NILs. This was evident from performance of NILs vis-à-vis PB1, which showed significant deviation for traits such as SF, GY, DF and PH particularly during 2019 and 2017 seasons. In all the seasons, however, the donor (N22) and the positive control (IR86918-B-B-305) remained the best performers under drought possibly indicating the necessity to transfer additional QTLs from multiple donors to further push the tolerance threshold of PB1. Similar observations were made by Vikram et al. [15] that qDTY1.1 had a considerable effect on enhancing the GY under both stressed and unstressed conditions. Similar effect of qDTY1.1 on improving the yield was earlier reported, conditioned through several adventitious drought-linked traits such as root characteristics, relative water content, biomass, and osmotic adjustment [17,60,61]. Meta QTLs have also been documented for maximum root length and GY within the qDTY1.1 region [62,63]. In the PB1 background, qDTY1.1 seems to enhance the yield under drought indirectly by modulating spikelet fertility as inferred from stepwise regression analysis. This implies on the importance of maintaining higher spikelet filling in achieving higher yields. Therefore, we hypothesize that grain filling could be the putative trait, qDTY1.1 is associated with, that require further validation.
Stability of yield performance under stress situations across the seasons is also important for a NIL to be considered for varietal evaluation prior to release as a cultivar across different rice growing areas of the country. The AMMI stability analysis carried out towards assessing the stable yielding potential, revealed that the NILs showed differential stability pattern. Two concomitant stability indices AST and YSI were used for assessing the yield stability, while AST alone was used for SF. Considering the stability indices and mean performance, one of the NILs, P1882-12-111-12 was adjudged stable for both GY and SF. This line can be a potential candidate for varietal evaluation pipeline in the future.
Among the agronomic traits, it was found that there was a delay in DF found across the NILs under stressed conditions. Delay in the flowering time under drought was also reported earlier in rice [16,64,65], which is due to delay in flower development and slow rate of panicle elongation under stress. The delay was however conspicuous in IR64, the short duration high yielding check variety used in the study. In short duration varieties with little drought adaptation, flower development delay occurs more significantly than longer duration cultivars. Similarly, there was a general decrease in PH among all the lines under stress, the degree of decrease was slightly reduced among the NILs. However, there was no increase in height observed among the NILs, although they have all possessed qDTY1.1, attributable to the tight linkage of qDTY1.1 with SD1 gene, the gene responsible for tallness in N22 [66]. In the present study, we surmise that linkage between the qDTY1.1 and sd1 allele might have broken resulting in semi-dwarf NILs with height similar to that of PB1.
Tolerance indices are better yardsticks for judging the tolerance of genotypes than the yield per se under stress condition. Three popular stress indices were utilized in the present study to identify best performing NILs. DYI is based on the mixed model which accounts for the genotype × stress level interaction across different environments. For drought prone areas, the combination of DYI with deviations in genotype performance under irrigated conditions may facilitate breeders to select genotypes with no yield reduction under favorable environment in comparison to currently cultivated varieties [37]. STI sorts out the entries which perform well under nonstress and fairly well under stress condition, high value of STI for a genotype infers higher tolerance level to the drought stress [39]. Fisher and Maurer [39] proposed SSI that estimates the yield reduction due to unfavorable environment in comparison to favorable environment. Lower the SSI value lesser the yield difference between stress and non-stress condition and hence greater is the tolerance of the genotype [67]. Thus, DYI, STI and SSI favored P1882-12-111-17 as best genotype for drought situation due to its highest grain yield under both stressed and unstressed conditions followed by P1882-12-111-3, P1882-12-111-6, P1882-12-111-7, P1882-12-111-12, P1882-12-111-15 and P1882-12-111-20. Further, P1882-12-111-17 out yielded majority of the remaining NILs. Among various traits measured in the present study, SF stood out as the most significant trait influencing grain yield under drought stress. Successful reproduction and efficient grain filling are the important determinants of ultimate grain yield under stress which is reflected as higher spikelet fertility. Thus, the NILs showed superior yields due to the maintenance of spikelet fertility under drought stress.
PB1 being a Basmati rice variety having specialty grain characteristics, transfer of drought tolerance from a non-aromatic variety with totally different grain architecture was a great challenge. Since the QTL transfer was into a Basmati background from a non-Basmati donor, a critical comprehensive quality check was done under both non-stress and stress situations in 2019. The NILs generated, possessed all the grain and cooking quality attributes of PB1 along with drought tolerance (Figure 4).
The recovery of quality in the NILs was achieved by stringent phenotypic selection integrated to the marker-based selection for the target trait, RSDS. NILs under well-watered treatment has grain quality similar to that of PB1, while an inconspicuous quality variation was noticed under drought situation. In spite of the successful recovery of PB1 genome component governing grain quality among the NILs, there was a marginal reduction in head rice recovery noticed under severe stress, although statistically insignificant, together with a slight increase in chalkiness. It is well known that drought at grain filling interferes with starch packaging in the grains, leading to chalkiness. Besides, the grain quality under stressed treatment was not significantly different from that under unstressed, indicating that drought had little effect of grain quality in the NILs. Grain quality recovery using MABB with augmented phenotypic selection, particularly in Basmati cultivars, was earlier reported from several studies [27,30,68].

5. Conclusions

In the present investigation, we have developed NILs of the popular Basmati rice variety, PB1 carrying a major QTL, qDTY1.1 for RSDS tolerance. This paper also forms the first ever report of successful transfer of a drought tolerance QTL into a Basmati rice cultivar. The improved drought tolerant NILs of PB1, were also combined with high yield and grain quality. These NILs could be a good alternative for the Basmati growing regions with limited soil moisture regimes. Additionally, these may help farmers to reduce number of irrigations without foregoing yield and grain quality potential of Basmati rice. These NILs are assessed for three consecutive seasons to assess their stability and drought response to identify potential candidates to be deployed for varietal testing and cultivar release. Moreover, these NILs can also serve as improved donor lines for imparting drought tolerance in future Basmati breeding programs.

Supplementary Materials

Supplementary materials can be found at https://0-www-mdpi-com.brum.beds.ac.uk/2073-4395/11/2/202/s1. Figure S1. Graphical genotypes of Pusa Basmati 1 near isogenic lines showing recovery of carrier chromosome of qDTY1.1 (Chromosome 1); Table S1: List of markers used for foreground and background selection in the present study; Table S2: Agro-morphological parameters indicating the relative performance PB1 NILs and parents under stressed (S) and unstressed (NS) conditions during 2017, 2018 and 2019; Table S3: Drought tolerance indices and grain yield of PB1 NILs along with checks and parents under stressed and unstressed conditions during 2017, 2018 and 2019.

Author Contributions

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

Funding

This research was funded by 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 Materials.

Acknowledgments

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

Conflicts of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Abbreviations

%RPercent reduction in yield
AMMIAdditive main effects and multiplicative interaction
ANOVAAnalysis of variance
ASAlkali spreading value
ASTAMMI stability value
BLUPBest linear unbiased predictor
DFDays to 50% flowering
DTYYield under drought
DYIDrought yield index
GYGrain yield
ICARIndian Council of Agricultural Research
KERKernel elongation ratio
KLACKernel length after cooking
KLBCKernel length before cooking
LBRLength by breadth ratio
MABBMarker assisted backcross breeding
N22Nagina 22
NILNear isogenic line
NSUnstressed
NTNumber of tillers
PB1Pusa Basmati 1
PHPlant height
PLPanicle length
QTLQuantitative trait loci
RPRecurrent parent
RPGRRecurrent parent genome recovery
RSDSReproductive stage drought stress
SStressed
SFSpikelet fertility
SSIStress susceptibility index
SSRSimple sequence repeats
STIStress tolerance index
YSIYield stability index

References

  1. Khanna, A.; Sharma, V.; Ellur, R.K.; Shikari, A.B.; Gopala Krishnan, S.; Singh, U.D.; Prakash, G.; Sharma, T.R.; Rathour, R.; Variar, M.; et al. Marker assisted pyramiding of major blast resistance genes Pi9 and Pita in the genetic background of an elite Basmati rice variety, Pusa Basmati 1. Indian J. Genet. 2015, 75, 417–425. [Google Scholar] [CrossRef] [Green Version]
  2. APEDA. India Export of Agro Food Products—Basmati Rice. 2020. Available online: https://agriexchange.apeda.gov.in/indexp/Product_description_32head.aspx?gcode=0601 (accessed on 11 December 2020).
  3. Trenberth, K.E.; Dail, A.; van der Schrier, G.; Jones, P.D.; Barichivich, J.; Briffa, K.R.; Sheffield, J. Global warming and changes in drought. Nat. Clim. Chang. 2013, 4. [Google Scholar] [CrossRef]
  4. Masson-Delmotte, V.; Zhai, P.; Pörtner, H.-O.; Roberts, D.; Skea, J.; Shukla, P.R.; Pirani, A.; Moufouma-Okia, W.; Péan, C.; Pidcock, R.; et al. Global Warming of 1.5 °C, An IPCC Special Report; IPCC Secretariat: Geneva, Switzerland, 2019; p. 616. [Google Scholar]
  5. Wassmann, R.; Jagadish, S.V.K.; Sumfleth, K.; Pathak, H.; Howell, G.; Ismail, A.; Serraj, R.; Redoña, E.; Singh, R.K.; Heuer, S. Regional vulnerability of climate change impacts on Asian rice production and scope for adaptation. Adv. Agron. 2009, 102, 91–133. [Google Scholar] [CrossRef]
  6. Bates, B.C.; Kundzewicz, Z.W.; Wu, S.; Palutikof, J.P. (Eds.) Climate Change and Water, Technical Paper of the Intergovernmental Panel on Climate Change; IPCC Secretariat: Geneva, Switzerland, 2008; p. 210. [Google Scholar]
  7. Pathak, H.; Tewari, A.N.; Sankhyan, S.; Dubey, D.S.; Mina, U.; Singh, V.K.; Jain, N.; Bhatia, A. Direct-seeded rice: Potential, performance and problems—A review. Curr. Adv. Agric. Sci. 2011, 3, 77–88. [Google Scholar]
  8. Dar, M.H.; Waza, S.A.; Shukla, S.; Zaidi, N.W.; Nayak, S.; Hossain, M.; Kumar, A.; Ismail, A.B.; Singh, U.S. Drought tolerant rice for ensuring food security in eastern India. Sustainability 2020, 12, 2214. [Google Scholar] [CrossRef] [Green Version]
  9. O’Toole, J.C. Adaptation of rice to drought prone environments. In Drought Resistance in Crops with Emphasis on Rice; IRRI: Los Baños, Philippines, 1982; pp. 195–213. [Google Scholar]
  10. Venuprasad, R.; Lafitte, H.R.; Atlin, G.N. Response to direct selection for grain yield under drought stress in rice. Crop Sci. 2007, 47, 285–293. [Google Scholar] [CrossRef]
  11. Kumar, R.; Sarawgi, A.K.; Ramos, C.; Amarante, S.T.; Ismail, A.M.; Wade, L.J. Partitioning of dry matter during drought stress in rainfed lowland rice. Field Crops Res. 2006, 96, 455–465. [Google Scholar] [CrossRef]
  12. Hsiao, T.C.; O’Toole, J.C.; Yambao, E.B.; Turner, N.C. Influence of osmotic adjustment on leaf rolling and tissue death in rice (Oryza sativa L.). Plant Physiol. 1984, 75, 338–341. [Google Scholar] [CrossRef] [Green Version]
  13. Turner, N.C.; O’Toole, J.C.; Cruz, R.T.; Namuco, O.S.; Ahmad, S. Response of seven diverse rice cultivars to water deficits. I. Stress development, canopy temperature, leaf rolling and growth. Field Crops Res. 1986, 13, 257–271. [Google Scholar] [CrossRef]
  14. Vinod, K.K.; Gopala Krishnan, S.; 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 II. [Google Scholar] [CrossRef]
  15. Vikram, P.; Swamy, B.P.M.; Dixit, S.; Ahmed, H.U.; Sta Cruz, M.T.; 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]
  16. Ghimire, K.H.; Quiatchon, L.A.; Vikram, P.; Swamy, B.P.M.; Dixit, S.; Ahmed, H.; Hernandez, J.E.; Borromeo, T.H.; Kumar, A. Identification and mapping of a QTL (qDTY1.1) with a consistent effect on grain yield under drought. Field Crops Res. 2012, 131, 88–96. [Google Scholar] [CrossRef]
  17. Bernier, J.; Kumar, A.; Venuprasad, R.; Spaner, D.; Atlin, G.N. A large-effect QTL for grain yield under reproductive-stage drought stress in upland rice. Crop Sci. 2007, 47, 507–516. [Google Scholar] [CrossRef]
  18. Kumar, A.; Sandhu, N.; Dixit, S.; Yadav, S.; Swamy, B.P.M.; Shamsudin, N.A.A. Marker-assisted selection strategy to pyramid two or more QTLs for quantitative trait-grain yield under drought. Rice 2018, 11, 35. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Sandhu, N.; Dixit, S.; Swamy, B.P.M.; Raman, A.; Kumar, S.; Singh, S.P.; Yadaw, R.B.; Singh, O.N.; Reddy, J.N.; Anandan, A.; et al. Marker assisted breeding to develop multiple stress tolerant varieties for flood and drought prone areas. Rice 2019, 12, 8. [Google Scholar] [CrossRef] [Green Version]
  20. Muthu, V.; Abbai, R.; Nallathambi, J.; Rahman, H.; Ramasamy, S.; Kambale, R. Pyramiding QTLs controlling tolerance against drought, salinity, and submergence in rice through marker assisted breeding. PLoS ONE 2020, 15. [Google Scholar] [CrossRef] [Green Version]
  21. Dixit, S.; Singh, A.; Sandhu, N.; Bhandari, A.; Vikram, P.; Kumar, A. Combining drought and submergence tolerance in rice: Marker-assisted breeding and QTL combination effects. Mol. Breed. 2017, 37, 143. [Google Scholar] [CrossRef] [Green Version]
  22. 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] [Green Version]
  23. Dwivedi, P.; Ramawat, N.; Dhawan, G.; Gopala Krishnan, S.; Vinod, K.K.; Singh, M.P.; Nagarajan, M.; Bhowmick, P.K.; Mandal, N.P.; Perraju, P.; et al. Drought Tolerant near Isogenic Lines (NILs) of Pusa 44 Developed through Marker Assisted Introgression of qDTY2.1 and qDTY3.1 Enhances Yield under Reproductive Stage Drought Stress. Agriculture 2021, 11, 64. [Google Scholar] [CrossRef]
  24. Singh, V.K.; Singh, A.; Singh, S.P.; Ellur, R.K.; Singh, D.; Gopala Krishnan, S.; Bhowmick, P.K.; Nagarajan, M.; Vinod, K.K.; Singh, U.D.; et al. Marker assisted simultaneous but stepwise backcross breeding for pyramiding blast resistance genes Pi2 and Pi54 into an elite Basmati rice restorer line PRR78. Plant Breed. 2013, 132, 486–495. [Google Scholar] [CrossRef]
  25. Singh, A.K.; Gopala Krishnan, S.; Singh, V.P.; Mohapatra, T.; Prabhu, K.V.; Singh, N.K.; Sharma, T.R.; Nagarajan, M.; Vinod, K.K.; Singh, D.; et al. Marker assisted selection: A paradigm shift in Basmati breeding. Indian J. Genet. 2011, 71, 120–128. [Google Scholar]
  26. Khanna, A.; Sharma, V.; Ellur, R.K.; Shikari, A.B.; Gopala Krishnan, S.; Singh, U.D.; Prakash, G.; Sharma, T.R.; Rathour, R.; Variar, M.; et al. Development and evaluation of near-isogenic lines for major blast resistance gene(s) in Basmati rice. Theor. Appl. Genet. 2015, 128, 1243–1259. [Google Scholar] [CrossRef] [PubMed]
  27. Sagar, V.; Dhawan, G.; Gopala Krishnan, S.; Vinod, K.K.; Ellur, R.K.; Mondal, K.K.; Rathour, R.; Prakash, G.; Nagarajan, M.; Bhowmick, P.K.; et al. Marker assisted introgression of genes governing resistance to bacterial blight and blast diseases into an elite Basmati rice Variety, ‘Pusa Basmati 1509’. Euphytica 2020, 216, 16. [Google Scholar] [CrossRef]
  28. Patra, B.C.; Ray, S.; Ngangkham, U.; Mohapatra, T. Rice. In Genetic and Genomic Resources for Grain Cereals Improvement; Academic Press: Cambridge, MA, USA, 2016; pp. 1–80. [Google Scholar] [CrossRef]
  29. Murray, H.G.; Thompson, W.F. Rapid isolation of high molecular weight DNA. Nucleic Acids Res. 1980, 8, 4321–4325. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Ellur, R.K.; Khanna, A.; Yadav, A.; Pathania, S.; Rajashekara, H.; Singh, V.K.; Gopala Krishnan, S.; Bhowmick, P.K.; Nagarajan, M.; Vinod, K.K.; et al. Improvement of Basmati rice varieties for resistance to blast and bacterial blight diseases using marker assisted backcross breeding. Plant Sci. 2016, 242, 330–341. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Sundaram, R.M.; Vishnupriya, M.R.; Biradar, S.K.; Laha, G.S.; Reddy, A.G.; Rani, N.S.; Sarma, N.P.; Sonti, R.V. Marker assisted introgression of bacterial blight resistance in Samba Mahsuri, an elite indica rice variety. Euphytica 2008, 160, 411–422. [Google Scholar] [CrossRef]
  32. Van Berloo, R. GGT 2.0: Versatile software for visualization and analysis of genetic data. J. Hered. 2008, 99, 232–236. [Google Scholar] [CrossRef] [Green Version]
  33. Venuprasad, R.; Dalid, C.O.; Del Valle, M.; Zhao, D.; Espiritu, M.; Sta Cruz, M.T.; 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]
  34. Sood, B.C.; Siddiq, E.A. A rapid technique for scent determinations in rice. Indian J. Genet. 1978, 38, 2268–2271. [Google Scholar]
  35. Purchase, J.L.; Hatting, H.; Van Deventer, C.S. Genotype × environment interaction of winter wheat (Triticum aestivum L.) In South Africa: II. Stability analysis of yield performance. S. Afr. J. Plant Soil 2000, 17, 101–107. [Google Scholar] [CrossRef]
  36. Sabaghnia, N.; Dehghani, H.; Sabaghpour, S.H. Graphic analysis of genotype by environment interaction for lentil yield in Iran. Agron. J. 2008, 100, 3. [Google Scholar] [CrossRef]
  37. Raman, A.; Verulkar, S.B.; Mandal, N.P.; Variar, M.; Shukla, V.D.; Dwivedi, J.L.; Singh, B.N.; Singh, O.N.; Swain, P.; Mall, A.K.; et al. Drought yield index to select high yielding rice lines under different drought stress severities. Rice 2012, 5, 31. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Fernandez, G.C.J. Effective selection criteria for assessing plant stress tolerance. In Adaptation of Food Crop Temperature and Water Stress, Proceeding of 4th International Symposium, Shanhua and Taipei, Taiwan, 13–18 August 1992; Asian Vegetable and Research and Development Center: Taipei, Taiwan, 1992; pp. 257–270. [Google Scholar]
  39. Fischer, R.A.; Maurer, R. Drought resistance in spring wheat cultivars. I. Grain yield responses. Aust. J. Agric. Res. 1978, 29, 897–912. [Google Scholar] [CrossRef]
  40. Khush, G.S. What it will take to Feed 5.0 billion rice consumers in 2030. Plant Mol. Biol. 2005, 59, 1–6. [Google Scholar] [CrossRef] [PubMed]
  41. Dilley, M.; Chen, R.S.; Deichmann, U.; Lerner-Lam, A.L.; Arnold, M. Multi hazard exposure analysis. In Natural Disaster Hotspots: A Global Risk Analysis; The World Bank: Washington, DC, USA; Columbia University: New York, NY, USA, 2005; pp. 47–52. [Google Scholar] [CrossRef]
  42. Shamsudin, N.A.A.; Swamy, B.M.; Ratnam, W.; Sta Cruz, M.T.; Sandhu, N.; Raman, A.K.; Kumar, A. Pyramiding of drought yield QTLs into a high-quality Malaysian rice cultivar MRQ74 improves yield under reproductive stage drought. Rice 2016, 9, 21. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Septiningsih, E.M.; Pamplona, A.M.; Sanchez, D.L.; Neeraja, C.N.; Vergara, G.V.; Heuer, S.; Ismail, A.M.; Mackill, D.J. Development of submergence-tolerant rice cultivars: The Sub1 locus and beyond. Ann. Bot. 2009, 103, 151–160. [Google Scholar] [CrossRef] [Green Version]
  44. Joseph, M.; Gopala Krishnan, S.; Sharma, R.K.; Singh, A.K.; Singh, V.P.; 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]
  45. Gopala Krishnan, 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]
  46. Singh, A.; Singh, V.K.; Singh, S.P.; Ellur, R.K.; Singh, D.; Bhowmick, P.K.; Gopala Krishnan, S.; Nagarajan, M.; Vinod, K.K.; Mohapatra, T.; et al. Molecular breeding for the development of multiple disease resistance in Basmati rice. AoB Plants 2012, 2012. [Google Scholar] [CrossRef]
  47. Ellur, R.K.; Khanna, A.; Gopala Krishnan, S.; Bhowmick, P.K.; Vinod, K.K.; Nagarajan, M.; Mondal, K.K.; Singh, N.K.; Singh, K.; Prabhu, K.V.; et al. Marker-aided incorporation of Xa38, a novel bacterial blight resistance gene, in PB1121 and comparison of its resistance spectrum with xa13 + Xa21. Sci. Rep. 2016, 6, 29188. [Google Scholar] [CrossRef]
  48. Sagar, V.; Gopala Krishnan, S.; Dwivedi, P.; Mondal, K.K.; Prakash, G.; Nagarajan, M.; Singh, A.K. Development of Basmati rice genotypes with resistance to both bacterial blight and blast diseases using marker assisted restricted backcross breeding. Indian J. Genet. 2018, 78, 36–47. [Google Scholar] [CrossRef]
  49. Dixit, S.; Singh, A.; Sta Cruz, M.T.; Maturan, P.T.; Amante, M.; Kumar, A. Multiple major QTL lead to stable yield performance of rice cultivars across varying drought intensities. BMC Genet. 2014, 15, 16. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Bernier, J.; Kumar, A.; Spaner, D.; Verulkar, S.; Mandal, N.P.; Sinha, P.K.; Peeraju, P.; Dongre, P.R.; Mahto, R.N.; Atlin, G.N. Characterization of the effect of a QTL for drought tolerance in rice, qtl12.1, over a range of environments in the Philippines and eastern India. Euphytica 2009, 166, 207–217. [Google Scholar] [CrossRef]
  51. Mishra, K.K.; Vikram, P.; Yadav, R.B.; Swamy, B.P.M.; Dixit, S.; Sta Cruz, M.T.; Maturan, P.; Marker, S.; Kumar, A. qDTY 12.1: A locus with a consistent effect on grain yield under drought in rice. BMC Genet. 2013, 14, 12. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Singh, V.P.; Pratik, S.; Gopala Krishnan, S.; Singh, A.K. Role of Indian Agricultural Research Institute in genetic improvement of rice varieties in India. In Genetic Improvement of Rice Varieties in India; Sharma, S.D., Prasada, R.U., Eds.; Today and Tomorrow’s Printers and Publications: New Delhi, India, 2004; pp. 141–187. [Google Scholar]
  53. Singh, V.P.; Singh, A.K.; Mohapatra, T.; Gopala Krishnan, S.; Ellur, R.K. Pusa Basmati 1121—A rice variety with exceptional kernel elongation and volume expansion after cooking. Rice 2018, 11, 19. [Google Scholar] [CrossRef] [Green Version]
  54. Vikram, P.; Kadam, S.; Singh, B.P.; Lee, Y.J.; Pal, J.K.; Singh, S.; Singh, O.N.; Swamy, B.P.M.; Thiyagarajan, K.; Singh, S.; et al. Genetic diversity analysis reveals importance of green revolution gene (sd1 locus) for drought tolerance in rice. Agric. Res. 2016, 5. [Google Scholar] [CrossRef]
  55. Jagadish, S.V.K.; Craufurd, P.Q.; Wheeler, T.R. High temperature stress and spikelet fertility in rice (Oryza sativa L.). J. Exp. Bot. 2007, 58, 1627–1635. [Google Scholar] [CrossRef] [Green Version]
  56. Collins, N.C.; Tardieu, F.; Toberosa, R. Quantitative trait loci and crop performance under abiotic stress: Where do we stand? Plant Physiol. 2008, 147, 469–486. [Google Scholar] [CrossRef] [Green Version]
  57. Price, A.H.; Steele, K.A.; Moore, B.J.; Jones, R.G.W. Upland rice grown in soil filled chambers and exposed to contrasting water- deficit regimes: II. Mapping quantitative trait loci for root morphology and distribution. Field Crops Res. 2002, 76, 25–43. [Google Scholar] [CrossRef]
  58. Courtois, B.; Shen, L.; Petalcorin Carandang, S.; Mauleon, R.; Li, Z.K. Locating QTLs controlling constitutive root traits in the rice population IAC 165 × Co39. Euphytica 2003, 134, 335–345. [Google Scholar] [CrossRef]
  59. WMO. The Global Climate in 2015–2019; World Meteorological Organization: Geneva, Switzerland, 2019; p. 21. [Google Scholar]
  60. Gomez, S.M.; Boopathi, N.M.; Kumar, S.S.; Ramasubramanian, T.; Chengsong, Z.; Jeyaprakash, P.; Senthil, A.; Babu, R.C. Molecular mapping and location of QTLs for drought-resistance traits in indica rice (Oryza sativa L.) lines adapted to target environments. Acta Physiol. Plant. 2010, 32, 355–364. [Google Scholar] [CrossRef]
  61. Kanbar, A.; Shashidhar, H.E.; Hittalmani, S. Mapping QTL associated with root and related traits in DH population of rice. Indian J. Genet. 2003, 62, 287–290. [Google Scholar]
  62. Courtois, B.; Ahmadi, N.; Khowaja, F.; Price, A.H.; Rami, J.F.; Frouin, J. Rice root genetic architecture: Meta-analysis from a drought QTL database. Rice 2009, 2, 115–128. [Google Scholar] [CrossRef] [Green Version]
  63. 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]
  64. Lafitte, H.R.; Ismail, A.; Bennet, J. Abiotic stress tolerance in rice for Asia: Progress and the future. In New Directions for a Diverse Planet, Proceedings of the 4th International Crop Science Congress, Brisbane, QLD, Australia, 26 September–1 October 2004; Fischer, T., Turner, N., Angus, J., McIntyre, L., Robertson, M., Borrell, A., Lloyd, D., Eds.; The Regional Institute Limited: Gosford, NSW, Australia, 2004; p. 17. [Google Scholar]
  65. Pantuwan, G.; Fukai, S.; Cooper, M.; Rajatasereekul, S.; O’Toole, J.C. Yield response of rice (Oryza sativa L.) genotypes to drought under rainfed lowlands: 2. Selection of drought resistant genotypes. Field Crops Res. 2002, 73, 169–180. [Google Scholar] [CrossRef]
  66. Vikram, P.; Swamy, B.P.M.; Dixit, S.; Singh, R.; Singh, B.P.; Miro, B.; Kohli, A.; Henry, A.; Singh, N.K.; Kumar, A. Drought susceptibility of modern rice varieties: An effect of linkage of drought tolerance with undesirable traits. Sci. Rep. 2015, 5, 14799. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  67. Chauhan, J.S.; Tyagi, M.K.; Kumar, A.; Nashaat, N.I.; Singh, M.; Singh, N.B.; Jakhar, M.L.; Welham, S.J. Drought effects on yield and its components in Indian mustard (Brassica juncea L.). Plant Breed. 2007, 126, 399–402. [Google Scholar] [CrossRef]
  68. Singh, V.K.; Singh, B.D.; Kumar, A.; Maurya, S.; Gopala Krishnan, S.; Vinod, K.K.; Pal, M.; Ellur, R.K.; Bhowmick, P.K.; Singh, A.K. Marker-assisted introgression of Saltol QTL enhances seedling stage salt tolerance in the rice variety ‘Pusa Basmati 1’. Int. J. Genom. 2018, 2018, 8319879. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Breeding scheme used in the marker-assisted backcross program for the transfer of qDTY1.1 quantitative trait locus (QTL) in the background of the elite rice variety, Pusa Basmati 1. Recurrent parent genome recovery (RPGR) denotes recurrent parent genome recovery.
Figure 1. Breeding scheme used in the marker-assisted backcross program for the transfer of qDTY1.1 quantitative trait locus (QTL) in the background of the elite rice variety, Pusa Basmati 1. Recurrent parent genome recovery (RPGR) denotes recurrent parent genome recovery.
Agronomy 11 00202 g001
Figure 2. Foreground selection of PB1 near isogenic lines (NILs) for qDTY1.1 using the marker, RM 431.
Figure 2. Foreground selection of PB1 near isogenic lines (NILs) for qDTY1.1 using the marker, RM 431.
Agronomy 11 00202 g002
Figure 3. Graphical representation of the genotypes of 18 qDTY1.1 introgressed NILs of PB1. All the NILs had maximum recovery on the carrier chromosome 1. CHR: Chromosome.
Figure 3. Graphical representation of the genotypes of 18 qDTY1.1 introgressed NILs of PB1. All the NILs had maximum recovery on the carrier chromosome 1. CHR: Chromosome.
Agronomy 11 00202 g003
Figure 4. Grain and cooking quality of the representative NILs of Pusa Basmati 1 carrying qDTY 1.1.
Figure 4. Grain and cooking quality of the representative NILs of Pusa Basmati 1 carrying qDTY 1.1.
Agronomy 11 00202 g004
Table 1. Genome wide marker polymorphism between the recurrent parent, Pusa Basmati 1 (PB1) and donor parent, Nagina 22 (N22).
Table 1. Genome wide marker polymorphism between the recurrent parent, Pusa Basmati 1 (PB1) and donor parent, Nagina 22 (N22).
ChromosomeNo. of MarkersMarker Diversity * (%)
Total SurveyedPolymorphic
1101 1413.86
2741216.21
3641218.75
4611016.39
571912.67
68789.19
7551323.63
822940.90
918422.22
1013538.46
1134926.47
1249816.32
Total65111317.35
* Marker diversity = Number of polymorphic markers × 100/ Total number of markers. Foreground marker, RM431 not counted among the chromosome 1 markers.
Table 2. Progressive selection statistics for the development of PB1 near isogenic lines introgressed with the reproductive stage drought stress (RSDS) tolerance QTL, qDTY1.1.
Table 2. Progressive selection statistics for the development of PB1 near isogenic lines introgressed with the reproductive stage drought stress (RSDS) tolerance QTL, qDTY1.1.
GenerationNo. of Seedlings/Families TestedRPGR (%)
TotalQTL Positive SelectedRangeMeanAverage Gain
F123 (s)5150.050.00.0
BC1F139 (s)15178.0–86.773.423.4
BC2F137 (s)1194.794.721.3
BC3F125 (s)5195.6–97.496.51.8
BC3F2229 (s)955597.8–98.798.21.7
BC3F355 (f)5524*--
BC3F424 (f)2418*--
BC3F518 (f)181898.2–99.199.00.8
* Not tested; RPGR, Recurrent parent genome recovery; s, seeds; f, families.
Table 3. Analysis of variance (ANOVA) for agronomic performance of NILs in individual years and combined under two stress treatments, stressed and unstressed.
Table 3. Analysis of variance (ANOVA) for agronomic performance of NILs in individual years and combined under two stress treatments, stressed and unstressed.
TraitsSeasonsC.V (%)Component Variance
BlockGenotype (G)Stress (S)G × SErrorYear (Y)Y × SY × GY × S × G
DF20171.12.6148.5 **102.6 **27.9 **1.1
20181.91.9175.2 **23.0 *5.03.7
20191.12.9278.2 **192.0 **2.6 *1.4
Combined1.31.2276.1 *548.5 **15.6 **1.61443.4 **20.726.7 **9.9 **
PH20174.89.0220.8 **1.7E+04 **26.318.1
20182.82.4458.1 **1.1E+04 **22.5 **7.2
20192.08.2325.0 **1.2E+04 **14.8 *3.4
Combined2.86.43.9E+04 **930.2 **20.0 **7.23850.5 **204.236.9 **21.7 **
NT201714.717.04.8 *46.8 **2.22.6
201813.225.82.5255.7 **2.73.5
201910.614.64.9 *207.1 **4.1 *2.1
Combined11.99.7184.85.3 **2.12.3243.8162.33.53.4
PL20172.60.019.8 **58.5 **2.0 **0.5
20184.30.110.9 **536.1 **1.41.0
20193.41.714.7 **253.6 **0.60.8
Combined3.30.7728.0 **42.2 **1.8 **0.734.0 **60.1 **1.51.1
SF20175.0126.389.4 **2.1E+04 **37.0 **11.9
20186.823.986.9 **1.1E+04 **58.3 **23.8
20193.79.556.4 **4866.1 **21.5 **5.3
Combined7.21.7E+056.7E+08 **2.5E+05 **1.1E+05 **2.5E+043.0E+08 **1.7E+08 **3.4E+05 **1.6E+05 **
GY20177.43.3E+052.3E+05 **6.1E+08 **1.8E+05 **5.8E+04
20185.72.2E+046.9E+05 **4.1E+08 **2.5E+05 **3.2E+04
20193.13.1E+044.4E+05 **3.9E+08 **3.8E+05 **1.9E+04
Combined4.81.2E+041.3E+08 **1.5E+05 **1.3E+05 **3731.63.9E+08 **1.3E+08 **1.5E+05 **1.3E+05 **
DF, days to 50% flowering; PH, plant height in cm; NT, number of tillers; PL, panicle length in cm; SF, spikelet fertility in %; GY, grain yield in kg/ha. *, ** indicates significance at 5% and 1% level, respectively.
Table 4. Average agronomic performance of NILs, parents and checks under drought stressed (S) and unstressed (NS) conditions over three seasons.
Table 4. Average agronomic performance of NILs, parents and checks under drought stressed (S) and unstressed (NS) conditions over three seasons.
NILsDFPHNTPLSF
SNSSNSSNSSNSSNS
P1882-12-111-1108.5 a103.7 ab78.2 cd108.0 c12.0 a14.6 a24.1 b28.9 a50.1 c80.1 bc
P1882-12-111-2103.5 bc100.8 cd80.9 c104.9 cd12.4 a14.0 a26.6 a28.5 a62.3 ab83.7 abc
P1882-12-111-3103.7 bc100.2 d79.9 c103.2 cd11.8 a13.7 a25.7 ab28.2 a60.1 ab82.8 abc
P1882-12-111-4104.7 bc104.5 a80.6 c106.6 c11.5 a14.9 a25.3 ab28.2 a64.1 ab82.6 abc
P1882-12-111-5103.3 bc102.5 a–d80.7 c105.0 c11.3 a14.0 a25.5 ab29.0 a65.1 a85.0 ab
P1882-12-111-6102.2 c101.2 b–d80.0 c105.9 c11.6 a14.7 a25.5 ab28.5 a57.7 b86.7 ab
P1882-12-111-7102.7 bc101.2 b–d81.8 c105.6 c13.5 a15.7 a26.4 a29.1 a57.8 b84.8 abc
P1882-12-111-8105.0 bc102.7 a–d79.5 c105.7 c12.3 a14.5 a25.6 ab28.5 a59.9 ab79.9 bc
P1882-12-111-9103.7 bc103.0 a–d80.2 c106.9 c11.6 a13.8 a25.4 ab28.7 a59.5 ab83.2 abc
P1882-12-111-10105.5 b103.3 a–c78.7 cd103.3 cd10.9 a14.1 a24.9 ab28.4 a63.3 ab83.1 abc
P1882-12-111-11104.5 bc103.2 a–c79.0 c104.9 cd10.6 a15.5 a25.3 ab29.5 a57.2 b81.4 abc
P1882-12-111-12103.3 bc102.0 a–d82.2 c106.4 c10.6 a13.4 a25.9 ab29.4 a62.5 ab80.7 bc
P1882-12-111-13103.7 bc101.5 b–d81.3 c104.4 cd10.6 a13.4 a26.3 a29.0 a57.0 b80.1 bc
P1882-12-111-14103.0 bc100.7 cd77.4 cd100.1 d11.1 a13.9 a25.5 ab28.8 a60.3 ab82.2 abc
P1882-12-111-15102.3 c101.3 b–d81.5 c104.2 cd11.2 a13.1 a25.1 ab28.6 a65.5 a80.5 bc
P1882-12-111-16103.8 bc103.7 ab80.8 c106.1 c10.9 a12.8 a25.0 ab29.4 a57.3 b82.0 abc
P1882-12-111-17103.8 bc102.7 a-d81.4 c107.1 c12.1 a15.2 a25.7 ab29.6 a63.6 ab84.1 abc
P1882-12-111-20104.7 bc104.0 ab82.0 c106.2 c10.8 a13.8 a26.0 a29.0 a61.7 ab81.7 abc
PB1102.3 c100.7 cd80.4 c103.9 cd10.2 a14.3 a25.2 ab29.1 a50.9 c77.6 c
N2284.2 e82.3 f103.2 b126.5 b13.4 a14.3 a18.2 d23.0 b62.4 ab88.4 a
IR86918-B-B-30579.5 f78.8 g115.0 a136.2 a12.2 a14.1 a21.2 c24.3 b62.0 ab86.4 ab
IR6498.7 d87.7 e73.9 d89.9 e11.1 a14.3 a21.9 c23.9 b60.7 ab82.5 abc
DF, days to 50% flowering; PH, plant height in cm; NT, number of tillers; PL, panicle length in cm; SF, spikelet fertility in %; S, stressed; NS, unstressed. Means followed by same letters are statistically at par based on the least significant difference test at 95% confidence level.
Table 5. Drought tolerance and stability indices for grain yield among the NILs and checks averaged across seasons, with stability for spikelet fertility under stress.
Table 5. Drought tolerance and stability indices for grain yield among the NILs and checks averaged across seasons, with stability for spikelet fertility under stress.
GenotypesGrain Yield (kg ha−1)IndicesSpikelet Fertility (%)
StressedUnstressed%RDYISTI SSIASTYSIMeanAST
P1882-12-111-1572.3 h5545.0 abc89.81.90.11.145.842.052.66.1
P1882-12-111-21090.6 d–g5757.1 abc80.91.10.21.09.719.062.36.2
P1882-12-111-31295.8 a–f5943.7 abc78.20.90.21.09.112.060.23.0
P1882-12-111-4901.6 d–g5542.2 abc83.41.30.21.019.427.064.18.6
P1882-12-111-51236.6 b–g5839.8 abc78.31.00.21.011.017.065.12.3
P1882-12-111-61408.8 a–d5841.0 abc75.70.90.30.929.218.057.71.8
P1882-12-111-71393.3 a–e5871.8 abc76.41.00.31.034.923.057.85.4
P1882-12-111-8879.6 efg5656.5 abc84.01.40.11.03.119.059.93.3
P1882-12-111-9806.9 fg5324.4 bc84.81.30.11.116.728.059.54.4
P1882-12-111-10725.9 g5765.5 abc87.31.60.11.117.530.063.35.9
P1882-12-111-111122.3 d–g5678.1 abc80.11.00.21.02.812.057.22.4
P1882-12-111-121254.6 a–f5722.2 abc78.00.90.21.06.512.062.50.6
P1882-12-111-13866.2 fg5582.3 abc84.41.30.21.113.126.057.00.9
P1882-12-111-141050.0 d–g5547.9 abc81.00.90.21.033.430.060.47.5
P1882-12-111-151238.8 b–g5879.9 abc78.90.90.21.05.112.065.53.8
P1882-12-111-16955.8 d–g5938.4 abc83.81.10.21.034.732.057.45.3
P1882-12-111-171715.8 ab6122.9 a72.00.70.30.927.815.063.55.9
P1882-12-111-201105.2 d–g5622.7 abc80.40.90.21.029.827.061.80.8
PB1531.0 h5475.6 abc90.32.00.11.127.834.049.52.3
N221772.8 a5228.6 c65.80.60.30.838.020.062.44.6
IR86918-B-B-3051706.3 abc5748.6 abc71.10.80.30.969.125.062.01.1
IR641188.9 c–g6109.2 ab80.71.40.21.038.326.059.95.1
%R, percentage reduction in yield; DYI, drought yield index; STI, stress tolerance index; SSI, stress susceptibility index; AST, additive main effective and multiplicative interaction (AMMI) stability value; YSI, yield stability index. Means followed by same letters are statistically at par based on the least significant difference test at 95% confidence level.
Table 6. Stability model for PB1 NILs along with checks under stressed and unstressed treatments during Kharif seasons of 2017, 2018 and 2019.
Table 6. Stability model for PB1 NILs along with checks under stressed and unstressed treatments during Kharif seasons of 2017, 2018 and 2019.
TreatmentSeasonModelR2F-ValueRMSE
Unstressed20171874.43 + 35.49DF − 15.32PH + 25.66SF0.7114.8210.5
201814079.10 − 43.43DF − 52.14SF0.355.3430.2
20193621.95 − 13.90PH − 44.44TN + 51.84PL + 34.79SF0.524.7161.2
Pooled7453.59 − 11.62PH − 36.23TN0.176.4383.4
Stressed2017−920.62 + 29.31SF0.5121.4146.9
20183619.85 − 46.53DF − 17.27PH + 229.72TN + 16.22SF0.514.6303.5
2019−998.34 − 62.25DF + 139.83TN + 113.11PL + 71.43SF0.8220.6267.3
Pooled−4850.56 + 11.73PH + 165.36TN + 50.29SF0.6336.8348.3
F-value, variance ratio; RMSE, root mean square error; DF, days to 50% flowering; PH, plant height in cm; TN, number of tillers; PL, panicle length in cm; SF, spikelet fertility in %; R2, proportion of variation explained by the model.
Table 7. Grain and cooking quality traits of the NILs and the recurrent parent, PB1 along with checks during Kharif 2019 under stressed (S) and unstressed (NS) conditions.
Table 7. Grain and cooking quality traits of the NILs and the recurrent parent, PB1 along with checks during Kharif 2019 under stressed (S) and unstressed (NS) conditions.
EntryKLBCKBBCLBRKLACKERASV
NSSNSSNSSNSSNSSNS/S
P1882-12-111-17.80 a6.93 b–e1.73 cd1.60 de4.50 abc4.33 ef13.47 def12.53 bc1.73 f1.81 cd7
P1882-12-111-27.53 bcd7.2 ab1.67 d1.40 g4.52 ab5.14 ab14.20 a–d12.27 b–e1.88 cde1.70 fgh7
P1882-12-111-37.53 bcd6.6 fg1.73 d1.60 de4.35 a–d4.13 fg14.27 a–d12.67 b1.89 cd1.92 a7
P1882-12-111-47.40 c–f6.47 g1.80 c1.67 d4.11 d3.88 g14.00 a–e11.93 cde1.89 cd1.85 bc7
P1882-12-111-57.67 ab6.93 b–e1.67 d1.40 g4.60 a4.95 abc14.67 ab12.87 ab1.91 a–d1.86 abc7
P1882-12-111-67.53 bcd7.00 bcd1.67 d1.53 ef4.52 ab4.57 de14.80 a12.67 b1.96 ab1.81 cd7
P1882-12-111-77.47 b–e7.07 abc1.67 d1.40 g4.48 abc5.05 ab14.80 a13.53 a1.98 a1.92 a7
P1882-12-111-87.27 e–h7.13 ab1.67 d1.47 fg4.36 a–d4.86 bcd14.13 a–e12.93 ab1.94 abc1.81 cd7
P1882-12-111-97.33 d–g7.13 ab1.67 d1.40 g4.40 abc5.1 ab13.93 b–e11.93 cde1.9 bcd1.67 ghi7
P1882-12-111-107.40 c–f7.07 abc1.67 d1.47 fg4.44 abc4.82 bcd13.80 cde11.6 e1.86 de1.64 hi7
P1882-12-111-117.47 b–e7.27 a1.67 d1.40 g4.48 abc5.19 a14.00 a–e13.53 a1.88 cde1.86 abc7
P1882-12-111-127.47 b–e7.2 ab1.67 d1.47 fg4.48 abc4.91 a–d13.80 cde11.73 de1.85 de1.63 i7
P1882-12-111-137.60 abc7.07 abc1.67 d1.47 fg4.56 a4.82 bcd13.87 b–e12.53 bc1.82 e1.77 de7
P1882-12-111-147.13 gh7.00 bcd1.67 d1.60 de4.28 bcd4.38 ef13.60 c–f12.4 bcd1.91 a–d1.77 de7
P1882-12-111-157.20 fgh6.93b–e1.67 d1.67 d4.32 a–d4.16 fg13.33 ef11.87 de1.85 de1.71 efg7
P1882-12-111-167.07 h7.07 abc1.67 d1.53 ef4.24 cd4.61 cde12.87 f12.33 b–e1.82 e1.75 def7
P1882-12-111-177.20 fgh6.53 fg1.67 d1.53 ef4.32 a–d4.26 ef13.60 c–f12.47 bcd1.89 cd1.91 ab7
P1882-12-111-207.27 e–h6.8 def1.60 d1.53 ef4.54 ab4.43 ef14.33 abc11.6 e1.97 a1.71 efg7
PB17.47 b–e7.07 abc1.67 d1.47 fg4.48 abc4.82 bcd14.13 a–e11.93 cde1.89 cd1.69 f–i7
N225.53 j4.93 h2.33 a2.33 a2.37 g2.11 j8.00 h6.67 h1.45 h1.35 k5
IR86918-B-B-3055.67 j5.2 h2.00 b1.93 bc2.83 f2.69 i9.00 g8.13 g1.59 g1.56 j5
IR646.53 i6.53 fg2.00 b2.00 b3.27 e3.27 h9.67 g9.00 f1.48 h1.38 k5
S.Em0.120.130.040.050.130.170.40.360.030.03
KLBC, kernel length before cooking in mm; KBBC, kernel breadth before cooking in mm; LBR, length by breadth ratio; KLAC, kernel length after cooking in mm; KBAC, kernel breadth after cooking in mm; KER, kernel elongation ratio; ASV, alkali spreading value; SEm, standard error of mean; S, stressed; NS, unstressed. Means followed by same letters are statistically at par based on the least significant difference test at 95% confidence level.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Dhawan, G.; Kumar, A.; Dwivedi, P.; Gopala Krishnan, S.; 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. https://0-doi-org.brum.beds.ac.uk/10.3390/agronomy11020202

AMA Style

Dhawan G, Kumar A, Dwivedi P, Gopala Krishnan S, Pal M, Vinod KK, Nagarajan M, Bhowmick PK, Bollinedi H, Ellur RK, 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(2):202. https://0-doi-org.brum.beds.ac.uk/10.3390/agronomy11020202

Chicago/Turabian Style

Dhawan, Gaurav, Aruna Kumar, Priyanka Dwivedi, Subbaiyan Gopala Krishnan, Madan Pal, Kunnummal Kurungara Vinod, Mariappan Nagarajan, Prolay Kumar Bhowmick, Haritha Bollinedi, Ranjith Kumar Ellur, and et al. 2021. "Introgression of qDTY1.1 Governing Reproductive Stage Drought Tolerance into an Elite Basmati Rice Variety “Pusa Basmati 1” through Marker Assisted Backcross Breeding" Agronomy 11, no. 2: 202. https://0-doi-org.brum.beds.ac.uk/10.3390/agronomy11020202

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