Skip to main content

Genetic mapping of physiological traits associated with terminal stage drought tolerance in rice



Drought during reproductive stage is among the main abiotic stresses responsible for drastic reduction of grain yield in rainfed rice. The genetic mechanism of reproductive stage drought tolerance is very complex. Many physiological and morphological traits are associated with this stress tolerance. Robust molecular markers are required for detection and incorporation of these correlated physiological traits into different superior genetic backgrounds. Identification of gene(s)/QTLs controlling reproductive stage drought tolerance and its deployment in rainfed rice improvement programs are very important.


QTLs linked to physiological traits under reproductive stage drought tolerance were detected by using 190 F7 recombinant inbred lines (RIL) mapping population of CR 143–2-2 and Krishnahamsa. Wide variations were observed in the estimates of ten physiological traits studied under the drought stress. The RIL population was genotyped using the bulk- segregant analysis (BSA) approach. A total of 77 SSR polymorphic markers were obtained from the parental polymorphisms survey of 401 tested primers. QTL analysis using inclusive composite interval mapping detected a total of three QTLs for the physiological traits namely relative chlorophyll content (qRCC1.1), chlorophyll a (qCHLa1.1), and proline content (qPRO3.1) in the studied RIL population. The QTL, qPRO3.1 is found to be a novel one showing LOD value of 13.93 and phenotypic variance (PVE) of 78.19%. The QTL was located within the marker interval of RM22-RM517 on chromosome 3. Another novel QTL, qRCC1.1 was mapped on chromosome 1 at a distance of 142.8 cM and found to control relative chlorophyll content during terminal drought stress. A third novel QTL was detected in the population that controlled chlorophyll a content (qCHLa1.1) under the terminal stress period. The QTL was located on chromosome 1 at a distance of 81.8 cM and showed 64.5% phenotypic variation.


The three novel QTLs, qRCC1.1, qCHLa1.1 and qPRO3.1 controlling relative chlorophyll content, chlorophyll a and proline content, respectively were identified in the mapping population derived from CR 143–2-2 and Krishnahamsa. These 3 QTLs will be useful for enhancement of terminal drought stress tolerance through marker-assisted breeding approach in rice.


Rice is an excellent gift of nature and majority of the global population consume it as staple food. Rice production and its related activities are main source of livelihood and socioeconomic security to the people in South-east Asia. Currently, about 725 million tons of paddy is being produced, worldwide, per year, from 160.8 million hectares of rice area [1]. In recent years, rice cultivation has been challenged by major production constraints mainly because of the adverse effects of climate change in India. Drought is the main yield reducing factor among the major abiotic stresses in rain-fed rice ecologies, worldwide. The frequency and occurrences of this stress is unpredictable. This stress affects rice productions in about 42 mha of rainfed rice including 8 mha upland rice in Asia [2]. In India as well, during the recent years, the rainfall pattern affects rice production to a greater extent. Therefore, high yielding varieties showing tolerance to reproductive stage drought stress need to be developed for the rainfed rice growers. Gene mapping for various traits involved during terminal drought stress tolerance and transfer of the relevant genes into superior backgrounds are needed for the improvement of rainfed rice.

The genetics of drought stress tolerance is very complex in nature and associated with several quantitative traits including various physiological and biochemical traits, involved during vegetative and reproductive stages of rice crop. Various plant traits are associated in plant growth and development during different growth stages of the crop. Majority of these traits are highly affected by drought stress, particularly during terminal stage of the crop. Even though the progress in drought breeding is slow, several traits controlled by genes/QTLs under the stress are available [3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24]. Drought stress during terminal stage is very detrimental to the crop and reduces grain yield severely [25,26,27,28,29,30]. Many physiological traits show strong correlation with drought tolerance in rice [31,32,33,34,35,36,37]. The highly correlated physiological traits to drought tolerance observed during terminal stage need to be mapped for deployment in the drought stress improvement in rice. Few yield QTLs are reported for enhancement of yield under drought stress under terminal stage drought [23, 38,39,40,41,42,43]. However, robust markers linked to various physiological traits involved in the stress tolerance are needed for improvement against the stress in rice.

Several research results on physiological traits responses under drought stress have already been published. Results on chlorophyll content, proline content and leaf area showed positive association with grain yield under terminal drought stress [44]. Also, the higher intensity of drought always correlated with decrease in chlorophyll content and increase in proline content in wheat plants [45]. Relative chlorophyll content (RCC) is an important physiological parameter involved during stress tolerance in rice which measures the greenness of leaves [46, 47]. Nine putative QTLs on seven different chromosomes for the trait cell-membrane stability were reported earlier under vegetative stage drought stress in rice [48]. QTL mapping study in Bala/Azucena reported 24 QTLs controlling various morphological, physiological and root related traits under drought stress, explained 4.6 to 22.3% phenotypic variance [49]. Five consistent QTLs were reported for various morphological and physiological traits linked to drought stress tolerance during terminal stage in rice [50].

Though few QTLs reports are available for drought tolerance, the need of strong molecular markers linked to various physiological traits involved during terminal stage drought tolerance are needed for molecular breeding in rice. Therefore, a RIL population comprising of 190 F7 lines was developed from the cross of a drought tolerant donor (CR 143–2-2) and susceptible genotype (Krishnahamsa) to map the gene(s)/QTL(s) responsible for various physiological traits conferring drought tolerance during terminal stage in rice.


Estimation of physiological traits of the mapping population under terminal stage drought stress

Reproductive stage drought stress affects rice yield drastically as flowering stage is the most critical stage of rice crop. Estimates of the ten physiological traits showed significant variations in the contrasting parents recorded during both the years (Table 1). All the physiological parameters viz., chlorophyll a, chlorophyll b, relative chlorophyll content, chlorophyll a + b, chlorophyll a/b, proline content, cell membrane stability, flag leaf width, biomass and per se yield were comparatively higher in the drought tolerant parent, CR143–2-2 compared to the susceptible parent except the trait, flag leaf length. Therefore, the selection of tolerant and susceptible parents for generation of mapping population may be effective.

Table 1 Mean value of estimates the physiological traits of contrasting parents under normal and reproductive stage drought stress situations during, 2014 and 2015

Leaf area is a major trait for controlling osmosis and photosynthesis in plants. In our study, leaf length (LL) and width (LW) showed significant variations among the RILs (Table 2). LL and LW ranged from 21.5 to 46.24 cm and 0.58 to 1.5 cm, respectively. Coefficient of variations obtained for the LL and LW in the RILs were 9.6 and 9.3, respectively (Table 2). The tolerant parent, CR 143–2-2 had shorter leaf length and wider leaf width of 27.96 cm and 1.29 cm, respectively while Krishnahamsa showed longer leaf length and narrow leaf width (Table 1).

Table 2 Mean statistical parameters of ten physiological traits and grain yield of RILs under reproductive stage drought stress during 2014 and 2015

The relative chlorophyll content of tolerant the parent was higher (36.06) compared to Krishnahamsa (27.47) (Table 1). RCC estimates of RILs showed wide variation under the stress and ranged from 20.14 to 45.95 with a mean of 30.61 (Table 2). The trait showed LSD5% and coefficient of variation to be 8.9 and 15.0, respectively. A wide variation in proline content was also observed in both the parents, CR 143–2-2 (160.62 μm/g) and Krishnahamsa (41.45 μm/g). In addition, proline content in the recombinant inbred lines varied from 6.1 to 228.5 μm/g with a mean of 74.88 μm/g (Table 1). Heritability (broad-sense) was found to be high for the trait showing maximum value of 0.95. The genetic advance and genetic advance over mean exhibited higher values of 119.12 and 149.58, respectively for proline content (Table 2).

Higher estimates of chlorophyll a (CHLa) and chlorophyll b (CHLb) content were estimated from the donor parent, CR 143–2-2 based on fresh weight basis with values of 3.45 and 0.83 mg/g, respectively. The sensitive parent, Krishnahamsa, showed relatively low values of 1.87 and 0.58 mg/g, respectively (Table 1). Additionally, chlorophyll a/b (CHLa/b) and chlorophyll a + b (CHLa+b) were higher in the tolerant parent than in the sensitive one (Table 1). Higher diversity was noticed in the estimates of CHLa, CHLb, CHLa+b and CHLa/b in the RILs (Table 2). Higher PCV, GCV and heritability values were observed for CHLb with value of 55.16, 54.13 and 0.96, respectively (Table 2). Cell membrane stability (CMS) varied widely which ranged from 6.74 to 93.81% with mean value of 63.42% using the RIL population. The donor parent showed higher value of CMS (87.73%) compared to the drought sensitive parent (51.14%). Heritability (broad-sense) and genetic advance for the CMS were found to be 0.95 and 38.95, respectively. High values of phenotypic covariance (PC), environmental covariance (EC) and coefficient of variation for physiological traits were estimated from the mapping population indicated an ideal mapping population used for the physiological traits (Table 2).

Frequency distributions

The distributions of recombinant inbred lines and parents based on the estimates of the ten physiological traits are shown in the Fig. 1. Both the parental lines are in the figures are depicted as P1 (tolerant parent) and P2 (susceptible parent) are placed wide apart from each other based on the studied physiological traits. Kurtosis and skewness values of the ten physiological traits for construction of a normal curve are furnished in Table 2. A positive leptokurtic skewed distribution curve was obtained for seven physiological traits viz., biomass, grain yield, chlorophyll b, chlorophyll a/b, relative chlorophyll content and proline content. However, three traits viz., leaf length, leaf width and chlorophyll a + b showed positive skewness estimates and negative kurtosis values. The leptokurtic distribution observed for cell membrane stability showed negative skewness and positive kurtosis. Negatively skewed platykurtic distribution was observed for the trait, chlorophyll a with both negative skewness and kurtosis value. All the studied traits except chlorophyll a and cell membrane stability showed almost a normal distribution of pattern under the stress condition (Fig.1).

Fig. 1

Frequency histogram and phenotypic distribution curves of ten physiological traits and grain yield generated from the RILs of CR 143–2-2 / Krishnahamsa

Nature of association of physiological traits and grain yield under the terminal drought stress

The correlation coefficient of the 10 physiological traits showed correlations among themselves and with grain yield (Table 3). Out of these correlation values, 26 correlations showed significant values at 0.01 and 8 were at 0.05 probability level. High positive correlations values were found for grain yield and biomass followed by chlorophyll a and chlorophyll a + b at 0.01 level of significance. Chlorophyll b and chlorophyll a/b showed a strong negative correlation (r = − 0.624**) at 0.01 probability level. Relative chlorophyll content, leaf length, leaf width, chlorophyll a and chlorophyll a/b showed significant positive correlation with grain yield under this stress condition (Table 3).

Table 3 Mean correlation coefficients of ten physiological traits and grain yield of RILs along with parents estimated under drought stress pooled over 2014 and 2015

Mapping of physiological traits involved in reproductive stage drought tolerance

In this investigation, a total of four hundred one microsatellite markers were utilized for detection of polymorphic markers between the parents (Table 4). Among the tested primers, 77 were detected to be polymorphic between the both the contrasting parents. Bulk-segregant analysis (BSA) stategy was followed by preparing two extreme inbred lines phenotypes bulks (B1: tolerant bulk and B2: susceptible bulk) and genotyped using the 77 polymorphic primers already obtained. These 77 primers were used in genotyping the RILs for mapping of the traits under the stress (Table 5; Fig. 2). ICIM (inclusive composite interval mapping) analysis revealed the presence three QTLs linked to relative chlorophyll content, proline content and chlorophyll a content under terminal drought stress situation (Table 6; Fig. 3a and b). These three QTLs showed LOD value ≥3.0 and were controlling three different physiological traits and located on two chromosomes. The QTL controlling relative chlorophyll content and chlorophyll a were found to be located on the chromosome 1 (Fig. 3). A QTL, qPRO3.1 controlling the trait, proline content was detected on the chromosome 3 (Fig. 3). High phenotypic variance of 78.19 and LOD value of 13.93 were obtained for proline content in the mapping population (Table 6). The location of qPRO3.1 was mapped on chromosome 3 at 21.2 cM within the marker interval of RM22 and RM517 (Table 6). A clear peak was observed for qPRO3.1 with additive effect of − 61.5. The QTL linked to the trait was detected in both the years’ phenotypic data (2014 and 2015) using ICIM software and showed the same marker interval for the QTL, qPRO3.1 (Table 6).

Table 4 Microsatellite markers obtained through the polymorphic analysis between CR143–2-2 and Krishnahamsa
Table 5 Details of polymorphic SSR markers obtained from bulk segregant analysis detected in the QTL mapping
Fig. 2

Representative electrophoregram obtained in different recombinant inbred lines using SSR markers. The numbers represent different RILs genotyped in the mapping study. Primer names are indicated in the right top corner position in each gel photos. P1: Tolerant parent, P2: Susceptible parent, M: 50 bp DNA ladder

Table 6 QTL identified by inclusive composite interval mapping using the physiological traits estimated during wet season 2014, 2015 and pooled over the years
Fig. 3

a QTL detected on chromosome 1 at LOD above the threshold level (violet colour represents QTL detected the Chlorophyll a and yellow colour represents QTL detected for relative chlorophyll content), b QTL detected on chromosome 3 at LOD above the threshold level (pink colour represents QTL detected for proline content)

A QTL, designated as qCHLa1.1 was found to be linked to the trait chlorophyll ‘a’ and present on chromosome 1 at 81.8 cM in the marker interval of RM495 and RM6703 (Table 6; Fig. 3). A small contribution of − 0.65 additive effect was obtained from this QTL mapping study (Table 6). The linkage was obtained showing phenotypic variance of 64.5% and LOD value of 4.13 for the trait, chlorophyll a (Fig. 3; Table 6). The QTL was detected using using both the years’ phenotypic data (2014 and 2015) by the ICIM software and revealed the same marker interval for the QTL, qCHLa. In addition, qRCC1.1 governing the trait, relative chlorophyll content was also detected from this QTL mapping study as an important QTL showing PVE% of 12.47 and LOD value of 4.76 and detected on the chromosome 1 (Table 6; Fig. 3). A clear peak was detected in map at142.8 cM position and located in the maker interval RM6703-RM3825 (Fig. 3). This linkage region was also observed by using the 2 years’ phenotypic data and revealed the same marker interval for qRCC1.1 (Table 6). An additive effect of − 2.35 was estimated from the analysis which is contributed by the QTL (Table 6).


The target physiological parameters were obtained from 190 inbred and parental lines under the stress showed wide variations among the RILs and between the parents. Few physiological traits showed strong correlations among themselves under the terminal drought stress condition. The frequency distribution curves were continuous for the studied physiological traits under the stress condition (Fig. 1). Hence, the mapping population used for tagging of the genes for the targeted physiological traits will be effective. Existence of genetic variation for relative chlorophyll content in rice genotypes under drought stress condition was also reported by earlier researchers [51]. Physiological traits namely chlorophyll a and chlorophyll b are the major traits which controll photosynthesis in plants. Reduction in photosynthetic rate under drought stress condition is attributed to chlorophyll a and b content in the leaves of rice plant [52].

In this study, large variations for the traits were observed in the RILs. Chlorophyll content of rice plant is an important secondary parameter for selection of suitable genotypes under drought stress condition due to its positive correlation with grain yield [33]. In our study, physiological traits viz., chlorophyll a/b and chlorophyll a showed significant positive correlation with grain yield may be useful for the QTL study. Relative chlorophyll content is an important physiological trait use to measure the greenness that enhances photosynthesis. The analysis of relative chlorophyll content (SPAD reading) showed a positive correlation of RCC under terminal drought stress [53]. In our investigation, a significant positive correlation also observed for relative chlorophyll content with grain yield under the stress (Table 3). Another important physiological trait found to be involved under the stress condition was proline content in leaves. Proline content in the leaves showed significant increase in the tolerant lines under drought stress that enhanced the plant growth and development [54]. Under this investigation, the frequency distribution of RILs showed to be normally distributed for proline content (Fig.1). A positive significant correlation values was observed for proline content and grain yield under the stress condition. Thus, the tolerant lines showed relatively better yield under the stress with increased proline content under the stress. Therefore, increased proline content under the stress may not be antagonistic to the grain yield. In addition, highly significant correlation of proline content with chlorophyll a and chlorophyll a + b was observed. Hence, proline content may be relied as a breeding parameter during selections of drought tolerant genotypes.

Among the physiological traits studied, three parameters viz., relative chlorophyll content, proline content and chlorophyll a content were linked in results mapped by QTL composite interval mapping. It is known that chlorophyll content is important for the photosynthetic process, hence highly correlated with photosynthetic rate. Low photosynthetic rate under drought stress is also associated with decrease in chlorophyll content. Low chlorophyll content in leaf reduces light absorbance in the sensitive genotypes than the tolerant genotypes. The drought sensitive genotypes could be damaged due to the damage in the photosynthetic apparatus, while drought tolerant genotypes were affected much less than tolerant genotypes. Higher lipid fluidity of the thylakoid membrane and lipid/protein ratio are seen in the tolerant genotypes, and thus increases the drought tolerance [55]. QTLs for chlorophyll content qCC-1, qCC-3 and qCC8 were detected on chromosome 1, 3 and 8, respectively within the marker interval of RG541- RG101, G62-G144 and RG598-RG418B under normal field condition in a DH population [56]. In our experiment, the linkage of the trait CHLa was detected at LOD value of 4.18 exhibiting phenotypic variance of 64.5%. The QTL, qCHLa1.1 was located in the SSR marker interval of RM495-RM6703 on chromosome 1. No other QTLs controlling drought tolerance were reported in this chromosomal region under terminal stage drought stress. Therefore, this QTL, qCHLa1.1 is considered as a novel QTL for controlling the trait chlorophyll a content in rice.

Another effective way to measure chlorophyll content via greenness of leaves was through SPAD meter reading. The study on quantitative trait locus for relative chlorophyll content (SPAD reading) was reported earlier by researchers [53, 57, 58]. They reported the QTL in the region of RM3916-RM2431 on the chromosome 4. Another study revealed that marker interval of RM302-RM472 on chromosome 1 was responsible for the physiological parameter for degree of greenness but undernormal field condition [59]. In our study, qRCC1.1 was detected showing a significant LOD value, PVE (%) and additive effect of 4.76, 12.47 and − 2.35, respectively for relative chlorophyll content on Chromosome 1 flanked by marker interval of RM6703-RM3825. However, there was no earlier report of QTLs in this location for relative chlorophyll content (RCC) trait underterminal stage drought stress. Thus, qRCC1.1 on chromosome 1 is a novel QTL controlling the trait, relative chlorophyll content at terminal drought stress situation.

Proline content in the leaves under terminal drought stress exhibited wide variation in the studied recombinant inbred lines (Table 2). High proline content is a good index for selection of tolerant genotypes under terminal drought stress tolerance in rice [60]. Proline content is an adaptive response for the accumulation of proline by plant tissue in stress condition like drought, salt and water stress [61]. Under drought condition, proline accumulation was suggested as a good parameter for drought resistance in plants [62]. The results on linkage analysis detected one QTL for proline content between marker interval of RM22-RM517 at position of 21.2 cM on the chromosome 3. The QTL showed a high LOD value and PVE (%) for the trait (Table 6). However, proline content showed high significant negative correlation with protein content but exhibited a strong positive correlation with catalase activity [60]. Experiment conducted using chromosome substitution lines for drought study revealed that the line with high proline content showed better drought tolerance and reported a segment on chromosome 1 for it [51]. The QTL detected by us was located on chromosome 3. As no earlier reports are available about any locus on the chromosome 3 controlling proline content, the QTL qPRO3.1 identified in our study is a novel QTL that can be useful in drought tolerance improvement in rice.

We observed a significant correlation (r = 0.404**) value between leaf length (LL) and leaf width (LW). These two traits were reported earlier researchers as important traits for drought tolerance which controlled canopy temperature and water use efficiency in rice [63]. But, no QTL was detected by the software for these two traits possibly due to the mapping population and mapping strategy used in this study. Also, no QTL was detected by this mapping technique for cell membrane stability (CMS). But, this trait was reported earlier as an drought breeding selection trait in cereals [48]. In this study, however, a negative correlation of − 0.134 was estimated for CMS with grain yield.


The nature of association of 10 studied physiological traits among themselves and with grain yield under the terminal drought stress were studied. Relative chlorophyll content, chlorophyll a and proline content showed strong asociation with grain yield under the terminal drought stress. Three QTLs viz., qRCC1.1, qCHLa1.1 and qPRO3.1 for relative chlorophyll content, chlorophyll a and proline content, respectively were detected which controlled tolerance under reproductive stage drought stress. These correlated traits will be useful as selection parameters in selecting desirable progenies for enhancement of terminal drought stress tolerance in rice. The QTLs and markers detected will be much useful in molecular breeding programs for enhancement of terminal stage drought tolerance in rice.


Plant materials

A total of 190 RILs along with susceptible (Krishnahamsa) and tolerant (CR 143–2-2) parents were taken as the experimental materials for the mapping study. The investigation was performed under the rain-out shelter, the controlled screening facility of ICAR-National Rice Research Institute (NRRI), Cuttack, Odisha during wet the seasons, 2014 and 2015. CR 143–2-2 is an early duration drought tolerant line developed by Institute for upland ecology used as the donor parent. The susceptible parent, Krishnahamsa (DRR Dhan 20) is an irrigated variety of Andhra Pradesh state of India. Both the contrasting parental lines for drought tolerance were obtained from NRRI, Gene Bank. The developed RIL lines at F7 generation were used for phenotyping of physiological traits and genotyping using microsatellite markers.

Phenotyping of the mapping population for physiological traits under terminal drought stress condition

All the recombinant inbred and parental lines were direct sown in an alpha lattice design using two replications during August month and irrigated up to panicle initiation stage. All the RILs were sown in 6 blocks accommodating 34 entries per block providing spacing of 10 × 15 cm. Both the parents were included in the total entries. Each row contains 25 hills per each recombinant line. Ten hill samples were collected for the evaluation of each RIL line. All the recombinant lines were grouped into three groups according to their flowering duration. Drought stress was applied at the primordium initiation (PI) stage to all the RILs and the parental lines. Fertilizer dose of 40:20:20 N:P:K was used in this the phenotyping experiment. Stress period was maintained throughout the reproductive stage and for excess stress, immediate irrigation was applied to maintain the stress up to -50kPA. Ten physiological traits viz., leaf length, leaf width, biomass, relative chlorophyll content, chlorophyll a, chlorophyll b, chlorophyll a/b, chlorophyll a + b, proline content, cell membrane stability and grain yield were computed under the terminal drought stress situation.

Days to 50% flowering of the parents and recombinant lines were recorded on whole plot basis. All other pre-harvest data viz., leaf length, relative chlorophyll content, and leaf width were recorded at 6–9 growth stage (SES 2014). SPAD-502 m was used to evaluate the greenness due to chlorophyll content under reproductive stage drought stress. Ten hills data were collected for recording of post-harvest data. For estimation of the studied traits namely chlorophyll a, chlorophyll b, chlorophyll a + b, chlorophyll a/b, the leaf samples were collected during mid-day. Chlorophyll content was computed by adopting the standard published method [64]. For cell membrane stability and proline content estimation, the samples were collected during mid-day period at growth stage 7–8 following the established protocols [65, 66], respectively. For estimating the relative chlorophyll content, SPAD meter or chlorophyll meter was used for recording the reading [46, 47].


DNA extraction

Leaves of 20 days old plants were collected aseptically from different pots for extraction of total genomic DNA. Extraction procedure started with the homogenization of leaf samples using liquid nitrogen in micro centrifuge tubes along with pre-warmed (65 °C) CTAB (Cetyltrimethyl ammonium bromide) extraction buffer (2% CTAB, 100 mM Tris pH 8, 1.3 M NaCl, 20 mM Ethylene diamine tetra acetate (EDTA) pH 8) followed by extraction in chloroform isoamyl alcohol, treatment in RNase and precipitating in ethanol as described in the protocol [67]. Final product of extracted DNA was checked by comparing it with λ-DNA on 1% agarose gel for the qualitative and quantitative purpose. Also, DNA isolated was measured for quantification and purity by UV visible spectrophotometer OD at 260 and 280 nm. The DNA samples were diluted uniformly to approximately 30 ng/μl and stored for use.

Polymerase chain reaction (PCR)

The PCR was performed in a thermal cycler (Applied-Biosystems) described in the publication [50] using the simple sequence repeat primers (Tables 4 and 5). The reaction mix for the PCR included 30 ng genomic DNA, 1 X reaction buffer (1.5 mM Tris-HCl, 50 mM KCl, 2 mM MgCl2), 10 mM dNTPs, 1 U Taq polymerase and 5 pmole each of forward and reverse primers. The final reaction mixture volume of 20 μl was performed for polymerase chain reaction (PCR). The thermal cycler settings starts with initial denaturation at 94 °C for 4 min, denaturation at 94 °C for 30 s for 35 cycles, primer annealing at 55 °C for 1 min and extension at 72 °C for 1.30 min; final extension at 72 °C for 10 min. After completion of amplification, PCR products were stored at − 20 °C and the amplified products were analyzed by electrophoresis using 3.5% agarose gel. The DNA fragments were then visualized by using ethidium bromide dye and the banding pattern was documented using gel documentation unit (Syngene GBox).

Bulk-segregant analysis (BSA)

This method is used to tag the presence of major QTLs linked to the trait of interest [2]. According to the phenotypic classification of the recombinant inbred line population, 10 RILs were bulked based on the extreme tolerant and susceptible phenotypes to detect the variation based on the use of polymorphic SSR markers. Using the software ICIM V4.0, the effect of QTLs and their relation with phenotypic and molecular proportion was analyzed.

Statistical analysis

From the physiological trait estimates of 190 RILs and their parents during wet seasons, 2014 and 2015, were used for the analysis of range, mean, skewness and kurtosis to determine the phenotypic distribution, main effect of RILs with the relative traits by employing SPSS v20.0 software [68]. Also the phenotypic correlation analysis and genetic advance among the RILs were performed by INDOSTAT software [69]. Analysis for environmental variance, genotypic coefficient of variance, phenotypic coefficient of variance, heritability, and F values were computed following the previous publications [70,71,72].

Linkage map and QTL analysis

Data on ten physiological traits (leaf length, leaf width, biomass, chlorophyll a, chlorophyll b, chlorophyll a + b, chlorophyll a/b, relative chlorophyll content, proline content and cell membrane stability) and grain yield of 190 recombinant inbred lines and parents were used for construction of linkage map as described in the earlier publication [50]. The linkage map was generated by employing ICIM (inclusive composite interval mapping) v4.0 software [73]. CIM analysis and additive effect were used to calculate the association of phenotypic and molecular proportions for the construction of the map. For map construction of all QTLs, a walking speed of 1.0 cM along the chromosomes, and LOD value of 3.0 was considered as threshold value along with 1000 permutation at P < 0.05. The naming of the QTLs were as per the standard nomenclatural guidelines published [74].

Availability of data and materials

The data generated or analyzed in this study are included in this article.



Bulk-segregant analysis


Cetyltrimethyl ammonium bromide




Ethylene diamine tetra acetate


Inclusive composite interval mapping


Logarithmic of odds




Quantitative trait loci


Relative chlorophyll content


Recombinant inbred line


  1. 1.

    United States department of agriculture (USDA). Foreign Agricultural Service, Office of Global Analysis. 2018.

    Google Scholar 

  2. 2.

    Venuprasad R, Dalid CO, Del Valle M, Zhao D, Espiritu M, Sta Cruz MT, Amante M, Kumar A, Atlin GN. 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–90.

    PubMed  Google Scholar 

  3. 3.

    Kamoshita A, Babu RC, Boopathi NM, Fukai S. Phenotypic and genotypic analysis of drought-resistance traits for development of rice cultivars adapted to rain-fed environments. Field Crops Resour. 2008;109:1–23.

    Google Scholar 

  4. 4.

    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(2):507–18.

    Google Scholar 

  5. 5.

    Li ZK, Xu JL. Breeding for drought and salt tolerant rice (Oryza sativa L.): progress and perspectives. Advances in molecular breeding toward drought and salt tolerant crops. Springer; 2007. p. 531–64.

  6. 6.

    Kumar R, Venuprasad R, Atlin GN. Genetic analysis of rainfed lowland rice drought tolerance under naturally occurring stress in eastern India: heritability and QTL effects. Field Crop Resour. 2007;103:42–52.

    Google Scholar 

  7. 7.

    Vikram P, BPM S, Dixit S, Ahmed HU, Cruz MTS, Singh AK, 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):1–15.

  8. 8.

    Yadaw RB, Dixit S, Raman A, et al. A QTL for high grain yield under lowland drought in the background of popular rice variety Sabitri from Nepal. Field Crop Res. 2013;144:281–7.

    Google Scholar 

  9. 9.

    Kam H, Fofana M, Zongo A, Taiwo S. Grain yield and related traits under drought and irrigated conditions. Int J Curr Adv Res. 2018;07(1):9059–66.

    Google Scholar 

  10. 10.

    Price A, Courtois B. Mapping QTLs associated with drought resistance in rice: progress, problems and prospects. Plant Growth Regul. 1999;29:123–33.

    CAS  Google Scholar 

  11. 11.

    Nguyen TL, Bui CB. Fine mapping for drought tolerance in rice (Oryza sativa L.). Omonrice. 2008;16:9–15.

    Google Scholar 

  12. 12.

    Sarvestani ZT, Pirdashti H, Sanavy SAM, Balouchi H. Study of water stress effects in different growth stages on yield and yield components of differentrice (Oryza sativa L.) cultivar. Pak J Biol Sci. 2008;11(10):1303–9.

    PubMed  Google Scholar 

  13. 13.

    Serraj R, Kumar A, McNally KL, Slamet-Loedin I, Bruskiewich R, Mauleon R, Cairns J, Hijmans RJ. Improvement of droughtresistance in rice. Adv Agronomy. 2009;103:41–99.

    CAS  Google Scholar 

  14. 14.

    Xu Q, Yuan XP, Yu HY, Wang YP, Tang SX, Wei X. Mapping QTLs for drought tolerance at seedling stage in rice using doubled haploid population. Rice Sci. 2011;18(1):23–8.

    Google Scholar 

  15. 15.

    Al-Shugeairy Z, Price AH, Robinson D. Genome wide association mapping for drought recovery trait in rice (Oryza sativa L.). Int J Appl Agric Sci. 2015;1(1):11–8.

    Google Scholar 

  16. 16.

    Fang Y, Xiong L. General mechanisms of drought response and their application in drought resistance improvement in plants. Cell Mol Life Sci. 2015;72(4):673–89.

    CAS  PubMed  Google Scholar 

  17. 17.

    Muthukumar C, Subathra T, Aiswarya J, Gayathri V, Chandra BR. Comparative genome-wide association studies for plant production traits under drought in diverse rice (Oryza sativa L.) lines using SNP and SSR markers. Curr Sci. 2015;109:139–47.

    Google Scholar 

  18. 18.

    Todaka D, Shinozaki K, Yamaguchi-Shinozaki K. Recent advances in the dissection of drought-stress regulatory networks and strategies for development of drought - tolerant transgenic rice plants. Front Plant Sci. 2015;I6:84–104.

    Google Scholar 

  19. 19.

    Sudeshna P, Sekhar Garg H, Mandi V, Sarkar KK, Bhattacharya C. Effect of water stress at tillering stage on different morphological traits of rice (Oryza sativa L.) genotypes. Int J Agric Sci Res. 2017;7(3):471–80.

    Google Scholar 

  20. 20.

    Hoang TH, Dinh LV, Nguyen TT, Ta NK, Gathignol F, Mai CD, Jouannic S, Tran KD, Khuat TH, Do VN, Lebrun M, Courtois B, Gantet P. Genome-wide association study of a panel of Vietnamese rice landraces reveals new QTLs for tolerance to water deficit during the vegetative phase. Rice. 2019;12(4):1–20.

  21. 21.

    Lanceras JC, Pantuwan G, Jongdee B, Toojinda T. Quantitative trait loci associated with drought tolerance at reproductive stage in rice. Plant Physiol. 2004;35:384–99.

    Google Scholar 

  22. 22.

    Yue B, Xue W, Xiong L, Yu X, Luo L, Cui K, Jin D, Xing Y, Zhang Q. Genetic basis of drought resistance at reproductive stage in rice: separation of drought tolerance from drought avoidance. Genetics. 2006;172(2):1213–28.

    PubMed  PubMed Central  Google Scholar 

  23. 23.

    Swamy BPM, Shamsudin NAA, Rahman SNA, Mauleon R, Ratnam W, Cruz MTS, Kumar A. Association Mapping of Yield and Yield-related Traits Under Reproductive Stage Drought Stress in Rice (Oryza sativa L.). Rice. 2017;10:21.

    PubMed  PubMed Central  Google Scholar 

  24. 24.

    Barik SR, Pandit E, Pradhan SK, Singh S, Swain P, Mohapatra T. QTL mapping for relative water content trait at reproductive stage drought stress in rice. Indian J Genet. 2018;78(4):401–8.

    CAS  Google Scholar 

  25. 25.

    Singh CS, Kumar B, Mehandi S, Chandra K. Effect of drought stress in Rice: a review on morphological and physiological characteristics. Trends Biosci. 2012;5(4):261–5.

    Google Scholar 

  26. 26.

    Kumar R, Sarawgi AK, Ramos C, Amarante ST, Ismail AM, Wade WJ. Partioning of dry matter during drought stress in rainfed lowland rice. Field Crop Res. 2006;96:455–65.

    Google Scholar 

  27. 27.

    Davatgar N, Neishabouri MR, Sepaskhah AR, Soltani A. Physiological and morphological responses of rice (Oryza sativa L.) to varying water stress management strategies. Int J Plant Prod. 2009;3(4):19–32.

    Google Scholar 

  28. 28.

    Kumar R, Sreenu K, Singh N, Jain N, Singh NK, Rai V. Effect of drought stress on contrasting cultivars of rice. Int J Tropic Agric. 2015;33(2):1559–64.

    Google Scholar 

  29. 29.

    Torres RO, Henry A. Yield stability of selected rice breeding lines and donors across conditions of mild to moderately severe drought stress. Field Crop Res. 2018;220:37–45.

    Google Scholar 

  30. 30.

    Maisura AC. M, Lubis I, Junaedi a, Ehara H. some physiological character responses of rice under drought conditions in a paddy system. J Int Soc Southeast Asian Agric Sci. 2014;20(1):104–14.

    Google Scholar 

  31. 31.

    Dharmappa PM, Doddaraju P, Malagondanahalli MV, et al. Introgression of Root and Water Use Efficiency Traits Enhances Water Productivity: An Evidence for Physiological Breeding in Rice (Oryza sativa L.). Rice. 2019;12(14)

  32. 32.

    Gomez SM, Satheesh Kumar P, Jeyaprakash P, Suresh R, Biji KR, Boopathi MN, Price AH, Babu RC. Mapping QTLs linked to physio-morphological and plant production traits under drought stress in rice (Oryza sativa L.) in the target environment. Am J Biochem Biotechnol. 2006;2(4):161–9.

    Google Scholar 

  33. 33.

    Lafitte HR, Price AH, Courtois B. Yield response to water deficit in an upland rice mapping population: associations among traits and genetic markers. Theo Appl Genet. 2004;109:1237–46.

    CAS  Google Scholar 

  34. 34.

    Panda RK, Pandit E, Swain A, Mohanty DP, Baig MJ, Kar M, Pradhan SK. Response of physiological and biochemical parameters in deeper rooting rice genotypes under irrigated and water stress conditions. Oryza. 2016;53(4):422–7.

    Google Scholar 

  35. 35.

    Panda RK, Pandit E, Dash SK, Kar M, Pradhan SK. Comparison of morpho-physiological traits and root architecture of tolerant and susceptible rice genotypes both under phosphorus and water stressed and normal condition. Oryza. 2017;54(1):21–8.

    Google Scholar 

  36. 36.

    Pandit E, Panda RK, Sahoo A, Pani DR, Pradhan SK. Genetic relationship and structure analyses of root growth angle for improvement of drought avoidance in early and mid-early maturing rice genotypes. Rice Sci. 2020;27(2):124–32.

    Google Scholar 

  37. 37.

    Srinivasan S, Gomez SM, Kumar SS, Ganesh SK, Biji KR, Senthil A, Babu RC. QTLs linked to leaf epicuticular wax, physio-morphological and plant production traits under drought stress in rice (Oryza sativa L.). Plant Growth Regul. 2008;56(3):245–56.

    CAS  Google Scholar 

  38. 38.

    Swamy BPM, Vikram P, Dixit S, Ahmed HU, Kumar A. Meta-analysis of grain yield QTL identified during agricultural drought in grasses showed consensus. BMC Genomics. 2011;2:319.

    Google Scholar 

  39. 39.

    Ghimire KH. Identification and mapping of QTL (qDTY1.1) with a consistent effect on GY under RS. Field Crop Res. 2012;131:88–96.

    Google Scholar 

  40. 40.

    Swamy BPM, Ahmed HU, Henry A, Mauleon R, Dixit S, Vikram P. Physiological and gene expression analyses reveal multiple QTLs enhance the yield of rice mega-variety IR64 under drought. PLoS One. 2013;8:e62795.

    PubMed  Google Scholar 

  41. 41.

    Dixit S, et al. Multiple major QTL lead to stable yield performance of rice cultivars across varying drought intensities. BMC Genet. 2014;15:16.

    PubMed  PubMed Central  Google Scholar 

  42. 42.

    Kumar A, Dixit S, Ram T, Yadaw RB, Mishra KK, Mandal NP. Breeding high-yielding drought-tolerant rice: genetic variations and conventional and molecular approaches. J Exp Bot. 2014.

  43. 43.

    Noraziyah AAS, Swamy BPM, Wickneswari R, Teressa SC, Anitha R, Kumar A. Marker-assisted pyramiding of drought yield QTLs into a popular Malaysian rice cultivar, MR219. BMC Genet. 2016;17:30.

    Google Scholar 

  44. 44.

    Kumar S, Dwivedi SK, Singh SS, Jha SK, Lekshmy S, Elanchezhian R, Singh ON, Bhatt BP. Identification of drought tolerant rice genotypes by analysing drought tolerance indices and morpho-physiological traits. SABRAO J Breed Genet. 2014;46:217–30.

    Google Scholar 

  45. 45.

    Keyvan S. The effects of drought stress on yield, relative water content, proline, soluble carbohydrates and chlorophyll of bread wheat cultivars. J Anim Plant Sci. 2010;8:1051–60.

    Google Scholar 

  46. 46.

    Kariya K, Matsuzaki A, Machida H. Distribution of chlorophyll content in leaf blade of rice plant. Jpn J Crop Sci. 1982;51:134–5.

    Google Scholar 

  47. 47.

    Inada K. Spectral ratio of reflectance for estimating chlorophyll content of leaf. Jpn J Crop Sci. 1985;54:261–5.

    CAS  Google Scholar 

  48. 48.

    Tripathy JN, Zhang J, Robin S, Nguyen TT, Nguyen HT. QTLs for cell- membrane stability mapped in rice (Oryza sativa L.) under drought stress. Theo Appl Genet. 2000;100:1197–200.

    CAS  Google Scholar 

  49. 49.

    Gomez MS, Kumar SS, Jeyaprakash P, Suresh R, Biji KR, Boopathi AM, Price AH, Babu RC. Mapping QTLs linked to physio-morphological and plant production traits under drought stress in rice (Oryza sativa L.) in the target environment. Amer. J Biochem Biotech. 2006;2:161–9.

    Google Scholar 

  50. 50.

    Barik SR, Pandit E, Pradhan SK, Mohanty SP, Mohapatra T. Genetic mapping of morpho-physiological traits involved during reproductive stage drought tolerance in rice. PLoS One. 2019;14(12):e0214979

    CAS  PubMed  PubMed Central  Google Scholar 

  51. 51.

    Larkunthod P, Nounjan N, Siangliw JL, Oojinda T, Sanitchon J, Jongdee B, Theerakulpisut P. Physiological Responses under drought stress of improved drought tolerant rice lines and their Parents. Not Bot Hort Agrobio. 2018;46:679–87.

    CAS  Article  Google Scholar 

  52. 52.

    Jaleel CA, Manivannan P, Wahid P, Farooq M, Somasundaram R, Panneerselvam R. Drought stress in plants: a review on morphological characteristics and pigments composition. Int J Agric Biol. 2009;11:100–5.

    Google Scholar 

  53. 53.

    Gholizadeh A, Amin MSM, Anuar AR, Aimrun W. Evaluation of SPAD chlorophyll meter in two different rice growth stages and its temporal variability. Euro J Sci Res. 2009;37:591–8.

    Google Scholar 

  54. 54.

    Mishra BK, Chaturvedi GS. Flowering stage drought stress resistance in upland rice in relation to physiological, biochemical traits and yield. Int J Curr Microbiol App Sci. 2018;7(2):781–2.

    Google Scholar 

  55. 55.

    Zaharieva M, Gaulin E, Havaux M, Acevedo E, Monneveux P. Drought and heat responses on the wild wheat relative Segilops geniculata Roth: potential interest for wheat improvement. Crop Sci. 2001;41:1321–9.

    Google Scholar 

  56. 56.

    Teng S, Qian Q, Zeng D, Kunihiro Y, Fujimoto K, Huang D, Zhu L. QTL analysis of leaf photosynthetic rate and related physiological traits in rice (Oryza sativa L.). Euphytica. 2004;135:1–7.

    CAS  Google Scholar 

  57. 57.

    Lafitte HR, Yongsheng G, Yan S, Li ZK. Whole plant responses, key processes, and adaptation to drought stress: the case of rice. J Exp Bot. 2007;58:169–75.

    CAS  PubMed  Google Scholar 

  58. 58.

    Takai T, Kondo M, Yano M, Yamamoto T. A quantitative trait locus for chlorophyll content and its association with leaf photosynthesis in rice. Rice. 2010;3:172–80.

    Article  Google Scholar 

  59. 59.

    Bing Y, Ya XW, Li-Jun L, Yong-Zhong X. QTL analysis for flag leaf characteristics and their relationships with yield and yield traits in rice. Acta Genet Sin. 2006;33:824–32.

    Google Scholar 

  60. 60.

    Roy R, Mazumder PB, Sharma GD. Proline, catalase and root traits as indices of drought resistance in bold grained rice (Oryza sativa L.) genotypes. Afr J Biotech. 2009;8:6521–8.

    CAS  Google Scholar 

  61. 61.

    Izanloo A, Condon AG, Langridge P, Tester M, Schnurbusch T. Different mechanisms of adaptation to cyclic water stress in two south Australian bread wheat cultivars. Jour of Exp Bot. 2008;59:3327–46.

    CAS  Google Scholar 

  62. 62.

    Singh KP, Singh K. Influence of stimulated water stress on free proline accumulation inTriticum aestivum L. Ind J Plant Phys. 1983;26:319–21.

    CAS  Google Scholar 

  63. 63.

    Biswal AK, Kohli A. Cereal flag leaf adaptations for grain yield under drought: Knowledge status and gaps. Mol Breed. 2013;31:749–66.

    Article  Google Scholar 

  64. 64.

    Arnon DI. Copper enzyme in isolated chloroplast polyphenoloxidase in Beta vulgaris L. Plant Physiol. 1949;24:1–15.

    CAS  PubMed  PubMed Central  Google Scholar 

  65. 65.

    Blum A, Ebercon A. Cell membrane stability as a measure of drought andheat tolerance in wheat. Crop Sci. 1981;21(1):43–7.

  66. 66.

    Sadasivam S, Manickam A. Biochemical Methods. New Delhi: New Age International Publishers (P) Ltd; 1996.

    Google Scholar 

  67. 67.

    Murray MG, Thomson WF. Rapid isolation of high molecular weight plant DNA. Nucl Acids Res. 1980;8:4321–5.

    CAS  PubMed  Google Scholar 

  68. 68.

    IBM Corp. IBM SPSS statistics for windows, version 20.0. Armonk: IBM Corp; 2011.

    Google Scholar 

  69. 69.

    Payne RW, Murray DA, Harding SA, Baird DB, Soutar DM. An introduction to Gen-Stat for windows. 14th ed. Hemel Hempstead: VSN International; 2011.

    Google Scholar 

  70. 70.

    Pradhan SK, Pandit E, Pawar S, Baksh SY, Mukherjee AK, Mohanty SP. Development of flash-flood tolerant and durable bacterial blight resistant versions of mega rice variety ‘Swarna’ through marker-assisted backcross breeding. Sci Rep. 2019;9:12810

    PubMed  PubMed Central  Google Scholar 

  71. 71.

    Pradhan SK, Nayak DK, Pandit E, Behera L, Anandan A, Mukherjee AK, Lenka S, Barik DP. Incorporation of bacterial blight resistance genes into lowland rice cultivar through marker-assisted backcross breeding. Phytopathology. 2016;

  72. 72.

    Pradhan SK, Pandit E, Pawar S, Naveenkumar R, Barik S R, Mohanty SP, Nayak DK, Ghritlahre SK, Sanjiba Rao D, Reddy JN, SSC Patnaik. Linkage disequilibrium mapping for grain Fe and Zn enhancing QTLs useful for nutrient dense rice breeding. BMC Plant Biol. 2020;

  73. 73.

    Wang J, Li H, Zhang L, Meng L. Users’ manual of QTL Ici-mapping. The quantitative genetics group, Institute of Crop Science, Chinese Academy of Agricultural Sciences (CAAS), Beijing 100081, China, and Genetic Resources Program, InternationalMaizeand Wheat Improvement Centre (CIMMYT). 2014; Apdo. Postal 6-641, 06600 Mexico, D.F., Mexico.

  74. 74.

    McCouch SR, Cho YG, Yano M, Paul E, Blinstrub M, Morishima H, Kinoshita T. Report on QTL nomenclature. Rice Genet Newsl. 1997;14:11–3.

    Google Scholar 

Download references


The authors are highly grateful to the Director, ICAR-National Rice Research Institute, and Head, Crop Improvement Division of the Institute for providing all the necessary facilities.


This work was internal project of the Institute under program1; project No.1.6 of ICAR- National Rice Research Institute Cuttack, Odisha, India. Institute fund was utilized for completion of the project work. No externally aided fund was received for this study.

Author information




SKP conceived the idea, SKP and SRB wrote, and EP and SKP revised the paper. SRB, EP, DKN and SPM generated phenotypic and genotypic data. SRB and EP performed data analyses. SKP interpreted the data. All authors reviewed and approved the final manuscript.

Authors’ information

Crop Improvement Division, ICAR-National Rice Research Institute, Cuttack-753006, Odisha, India.

Corresponding author

Correspondence to Sharat Kumar Pradhan.

Ethics declarations

Ethics approval and consent to participate

The authors declare that this study complies with the current laws of the countries in which the experiments were performed.

Consent for publication

Not applicable.

Competing interests

We declare that there is no competing interest among the authors.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Barik, S.R., Pandit, E., Mohanty, S.P. et al. Genetic mapping of physiological traits associated with terminal stage drought tolerance in rice. BMC Genet 21, 76 (2020).

Download citation


  • Bulk-segregant analysis
  • Chlorophyll content
  • Proline content
  • Physiological traits
  • QTL mapping
  • Relative chlorophyll content
  • Reproductive stage drought tolerance
  • Rice