Hassaan Ahmad1‡
Syed Adeel Zafar1†‡
Muhammad Kashif Naeem1
Sajid Shokat2
Safeena Inam1
Malik Attique ur Rehman1
Shahzad Amir Naveed3
Jianlong Xu3Zhikang Li3Ghulam Muhammad Ali1
Muhammad Ramzan Khan1*- 1National Institute for Genomics and Advanced Biotechnology (NIGAB), National Agricultural Research Centre, Islamabad, Pakistan
- 2Nuclear Institute for Agriculture and Biology, Faisalabad, Pakistan
- 3Institute of Crop Science, Chinese Academy of Agricultural Sciences, Beijing, China
Optimum soil water availability is vital for maximum yield production in rice which is challenged by increasing spells of drought. The reproductive stage drought is among the main limiting factors leading to the drastic reduction in grain yield. The objective of this study was to investigate the molecular and morphophysiological responses of pre-anthesis stage drought stress in green super rice. The study assessed the performance of 26 rice lines under irrigated and drought conditions. Irrigated treatment was allowed to grow normally, while drought stress was imposed for 30 days at the pre-anthesis stage. Three important physiological traits including pollen fertility percentage (PFP), cell membrane stability (CMS), and normalized difference vegetative index (NDVI) were recorded at anthesis stage during the last week of drought stress. Agronomic traits of economic importance including grain yield were recorded at maturity stage. The analysis of variance demonstrated significant variation among the genotypes for most of the studied traits. Correlation and principal component analyses demonstrated highly significant associations of particular agronomic traits with grain yield, and genetic diversity among genotypes, respectively. Our study demonstrated a higher drought tolerance potential of GSR lines compared with local cultivars, mainly by higher pollen viability, plant biomass, CMS, and harvest index under drought. In addition, the molecular basis of drought tolerance in GSR lines was related to upregulation of certain drought-responsive genes including OsSADRI, OsDSM1, OsDT11, but not the DREB genes. Our study identified novel drought-responsive genes (LOC_Os11g36190, LOC_Os12g04500, LOC_Os12g26290, and LOC_Os02g11960) that could be further characterized using reverse genetics to be utilized in molecular breeding for drought tolerance.
Introduction
Rice (Oryza sativa L.) is one of the primary staple food crops for nearly 50% of the world population (Zafar et al., 2018). The countries located in East Asia, South Asia, and Southeast Asia are dominant in production and consumption of rice across the globe. Historically, more than 90% of world rice production is contributed from these countries (Barker et al., 1999). Its production is needed to increase by 0.6%–0.9% per year until 2050 to feed the further 2 billion people (Desa, 2015). However, different abiotic and biotic stresses are major limiting factors for obtaining higher yield in rice (Zafar et al., 2017; Oliva et al., 2019; Ahmed et al., 2021). Being a water-loving plant, rice is highly sensitive to drought stress, which significantly affects its grain yield (Shuxing, 2014). Drought is becoming a serious yield constraint for various major crops due to global water scarcity (Wattoo et al., 2018; Hussain et al., 2019; Shokat et al., 2020a). A recent study using the yield and metrological data from 1980 to 2015 reported the yield decline up to 21% in wheat (Triticum aestivum L.) and 40% in maize (Zea mays L.) due to drought on a global scale (Daryanto et al., 2016). In rice, mild-drought stress reduced grain yield by 31%–64%, while severe stress reduced 65%–85% yield compared with normal conditions (Kumar et al., 2008). It affects the yield by altering different agronomic and physiological traits including plant height, number of tillers, leaf area, leaf rolling, transpiration rate, accumulation of osmoprotectants, root system, and stomatal closure (Nakashima et al., 2007; Islam et al., 2009; Tong et al., 2009). Anthesis stage drought stress can interrupt flowering, floret initiation (Bajji et al., 2002), pollen fertility (Zhou et al., 2011), and grain filling, resulting in poor paddy yield. Rice growth is affected by drought at different stages including booting (Shao et al., 2014), flowering (Liu et al., 2006), and grain filling stage (Zhang et al., 2018). However, drought stress at anthesis stage restricts the availability of photosynthates by disturbing the sink capacity (Do et al., 2010) and reduces the grain yield, plant biomass, and ultimately the harvest index (Blum, 2018). It also impairs anther dehiscence, pollen viability, and pollen germination in rice resulting in spikelet sterility and more sterile grains in the panicles (Prasad et al., 2017). Drought induced spikelet sterility is considered as one of the major causes of yield reduction.
To address the challenge, natural variation in rice germplasm for drought tolerance could be exploited to identify the drought-tolerant genotypes, the associated traits, and underlying genes (Panda et al., 2021). In addition, induced variation via hybridization and mutagenesis could serve as an important genetic resource for target breeding (Zafar et al., 2020c). For the purpose, scientists have started to put efforts to breed green super rice (GSR), an elite rice type that could withstand multiple stresses with high nutrient-use efficiency (Wing et al., 2018; Jewel et al., 2019). The idea was given by a famous rice geneticist Qifa Zhang in 2007 (Zhang, 2007), which was later implemented by a team of international scientists from China and the International Rice Research Institute (IRRI), Philippines (Yu et al., 2020). The present study was conducted to evaluate 22 selected GSR lines along with four local rice cultivars for drought tolerance in Pakistan, and identify agronomic and physiological traits associated with drought tolerance in GSR. In addition, the contrasting drought tolerant and sensitive lines were assessed for gene expression profile to identify underlying genes related to drought tolerance in GSR. This study identified high-yielding drought-tolerant GSR lines and provided us knowledge about drought tolerance-related traits, and novel drought-related genes.
Materials and methods
Experimental site
The field experiment was conducted at the National Institute for Genomics and Advanced Biotechnology, NARC, Pakistan (33.684°N and 73.048°E) during rice growing period (May–October, 2020). To minimize the water infiltration from control to drought plot, a 6- to 8-feet path was made between both plots, and furthermore, plastic film was applied under the soil surface with a depth of 60 cm.
Experimental design
The 22 diverse GSR lines were selected based on diverse phenotypic characteristics from the 552 GSR genotypes (Supplementary Table S1). Twenty-two GSR lines along with four checks were evaluated using split plot randomized complete block design with two treatments (well-watered and drought) each having three replications. Seeds were sown in nursery trays and 30-day-old seedlings were transplanted in the field. Each plot consisted of five rows of 10 plants with 30-cm row/row and plant/plant distance (Yugandhar et al., 2017). Both plots were irrigated normally (8–10 cm) until anthesis stage. Fertilizer, weedicide, and insecticide application was done according to recommended dosage. Crop cultivation was carried out according to normal cultural practices.
Drought was imposed for 30 days by withholding the applied water at the beginning of anthesis stage. Physiological traits were recorded during the last week of stress. After 30 days, the field was rewatered. At physiological maturity, five representative plants were selected for the measurement of agronomic traits from the three middle rows of each replication to avoid confounding border effects (Chaturvedi et al., 2017).
Physiological measurements
Cell membrane stability
Leaf samples were collected at the last week of drought stress to examine the cell membrane stability by recording the electrolyte conductivity using and electrical conductivity meter (HI 9811-5 Portable EC meter HANNA® Instruments, USA). Flag leaves from three plants per replicate (of each genotype) were collected from both control and drought stress fields in 20-ml glass vials. Further measurement was recorded as proposed by Tripathy et al. (2000). CMS was formulated as the reciprocal of cell membrane injury by using the following formula (Blum and Ebercon, 1981):
CMS% = {[1 − (T1 / T2)] / [1−(C1 / C2)]} × 100
where, T and C refer to stressed and controlled, respectively. C1 (initial control), T1 (initial stress), and after autoclave, C2 (final control), T2 (final stress) were the assumed conductance.
Normalized difference vegetation index
Normalized difference vegetation index (NDVI) is a spectral reflectance-based measure of the density of green vegetation on a land area. NDVI measurements were taken using GreenSeeker™ Handheld Optical Sensor Unit (NTech Industries, Inc., USA), keeping the sensor at 0.5–1 m above the central rows of all the genotypes individually in three replications of both control and stress field plots (Govaerts and Verhulst, 2010).
Pollen fertility test
About five to eight mature spikelets from five panicles (one from each plant) were collected in the morning before anthesis. Spikelets were fixed in FAA solution (formaldehyde:ethanol:acetic acid with a ratio of 1:18:1, respectively) until staining. Anthers were crushed with forceps on glass slide to release pollens, which were immersed in 1% potassium iodide (I2-KI) solution followed by observation under a light microscope (NIKON DIGITAL SIGHT DS-Fi2). Pollens that stained black and circular were considered fertile, while those stained red-orange and of irregular shape were considered sterile (Zafar et al., 2020b). Pollen fertility percentage (PFP) was calculated using the following formula:
Measurement of agronomic traits
Agronomic traits, including plant height per plant (PH), tillers per plant (TPP), grain yield per plant (GY), straw yield per plant (SY), total biomass per plant (TBM), 1,000-grain weight (TGW), and grain length (GL) were recorded manually. Harvest index (HI) was calculated as the ratio of GY to TBM. Drought susceptibility index (DSI) was calculated as [(1−Y / YP) / D] as described earlier (Khanna-Chopra and Viswanathan, 1999; Zafar et al., 2020a). Here, Y is the grain yield under stress conditions, and YP is the grain yield under normal conditions, while D represents the stress intensity, which was calculated as D = (1 − X / XP), where X and XP are means of Y and YP, respectively. Measurements for these traits were carried out on five randomly selected plants of each genotype from each replication by following the method (IRRI, 2002).
RNA extraction and cDNA synthesis
Total RNA was extracted from the panicles of selected drought-tolerant and -sensitive genotypes from both well-watered (WW) and drought-stressed plants. Panicles were harvested from plants and immediately kept in liquid nitrogen followed by storage at −80°C to avoid the denaturation of RNA. The PureLink RNA Mini kit (Thermo Fisher Scientific) was used to extract the total RNA, in accordance with the manufacturer’s protocol. The quality of isolated RNA was observed on 1.5% RNase-free agarose gel and quantified using the BioSpec-nano spectrophotometer. One microgram of total RNA was used to reverse transcribe into cDNA using RevertAid Reverse Transcriptase kit (Thermo-Fisher Scientific) following the manufacturer’s instructions.
Differentially expressed gene selection and quantitative real-time PCR
Ten differentially expressed genes (DEGs) under drought stress were selected from a comparative transcriptome study in rice (Huang et al., 2014). To our knowledge, these genes have not been studied before specifically for drought response. In addition, we studied the expressions of three known drought tolerance-related genes: Oryza sativa Salt-, ABA- and Drought-Induced RING Finger Protein 1 (OsSADR1) (Park et al., 2018), Drought-Hypersensitive Mutant1 (DSM1) (Ning et al., 2009), Drought tolerance 11 (OsDT11) (Li et al., 2017), OsDREB1E, and OsDREB2B (Chen et al., 2008). Selected genes are listed in Supplementary Table S2. Coding sequences (CDS) of the selected DEGs were retrieved from the Rice Genome Annotation project (http://rice.plantbiology.msu.edu/cgi-bin/gbrowse/rice/). A gene-specific pair of primers was designed using AmplifX version 1.7.0 software, and primer sequences are listed in Supplementary Table S3.
Quantitative real-time PCR (qRT-PCR) was carried out to determine the relative expression levels of 13 selected genes on StepOne™ Real-Time PCR System (Thermo Fisher Scientific) using Maxima SYBR Green. The delta cT method was used to calculate the relative expression level of each gene, and rice Actin1 gene was used to normalize the expression (Fang et al., 2021; Zafar et al., 2021).
Statistical analysis
Morphophysiological traits data were analyzed by analysis of variance using SPSS software according to split plot randomized complete block design. Principal component analysis (PCA) was done through the XL-STAT software (ver. 2018) to categorize various physiological and morphological traits (Mohammadi and Prasanna, 2003). Pearson’s correlation matrix analysis was done using the “cor” package in R studio. The p-values for the coefficient of correlation (r) were obtained by applying Student’s t-test with the “cor.test” function in R-studio. In the correlation matrix plot, only significant relationships were labeled with stars. Expression pattern significance was calculated using t-test.
Results
Analysis of variance showed significant variation among green super rice accessions under drought stress
Analysis of variance (ANOVA) was performed to see the significant differences of variation among the genotypes and water treatments for physiological and yield-related traits. ANOVA showed significant variation (p < 0.01) among the tested genotypes for PH, GY, HI, TGW, GL, and NDVI (Table 1), while nonsignificant differences were observed for TPP, SY, and PFP. There was no significant effect on the studied traits among the replications, which strengthen the reliability of this experiment. Drought significantly affected PH, GY, SY, HI, GL, PFP, and NDVI, while traits, such as TPP, TBM, and TGW, were not affected by drought. The genotype × environment interaction was also significant for PH, TPP, SY, HI, TGW, PFP, and NDVI (Table 1) where pronounced reduction was recorded under drought conditions. Since PH, GY, HI, GL, and NDVI displayed significant differences for genotypes as well as drought treatment, these traits could be key selection markers for drought tolerance screening in rice.
TABLE 1. Mean square values from the analysis of variance for the effect of genotype, environment, and their interaction on agronomic and physiological traits.
Significance levels are indicated: *p < 0.05; **p < 0.01. df, degrees of freedom; PH, plant height; TTP, tillers per plant; GY, grain yield; SY, straw yield; TBM, total biomass; HI, harvest index; TGW, thousand-grain weight; GL, grain length; PFP, panicle fertility percentage; NDVI, normalized difference vegetation index.
Principal component analysis revealed genetic variation among green super rice accessions under drought
Separated PCA analyses were performed to develop a trait–genotype (T–G) biplot and to detect genetic variation among the studied genotypes for various morphophysiological traits under well-watered and drought-stressed conditions. Under WW environment, a biplot was drawn between PC1 and PC2 explaining 31.3% and 25.6%, of total variation, respectively (Figure 1A). Our results indicate that SY, GY, HI, and TBM were in the opposite direction of NDVI, TGW, and GL indicating their opposite relationship with each other. In addition, the GSR lines were mostly clustered near the origin and show less genetic variability, while checks Kashmir Basmati, Kissan Basmati, and IR-64 were widely distributed apart from the origin and showed remarkable genetic variability (Figure 1A). In case of drought treatment, the PC1 alone accounted for 51.10% of the total variability, while PC2 shared 15.20% (Figure 1B). Results of this experiment show GY were clustered closer to PFP, CMS, TBM, and HI, while it was in opposite direction of NDVI and DSI. In contrast with the WW treatment, many GSR lines, namely, NGSR-3, NGSR-15, NGSR-18, NGSR-13, NGSR-21, and NIAB-IR-9, fall near the apex of the biplot and show remarkable genetic variation under drought stress (Figure 1B). The check varieties Kissan Basmati, NIAB-IR-9, and Kashmir Basmati also showed considerable genetic variability and reputation of these accessions for further selection in breeding programs.
FIGURE 1. PCA showing biplot for genotypes and studied traits (A) under normal condition and (B) under drought stress condition.
Mean performance of green super rice accessions for studied traits
Drought stress showed a remarkable reduction in grain yield and yield-related traits in all studied genotypes except NGSR-15 and NGSR-18 (Figures 2 and 3). Among the 22 GSR lines, the minimum PH (86.6 cm) under drought condition was attained by NGSR-8 and the maximum (103.7 cm) was recorded by NGSR-14, whereas among the four checks, the maximum PH was recorded for Kashmir Basmati (113.6 cm), and the minimum was depicted by NIAB-IR-9 (52.3 cm). All GSR lines demonstrated higher PH than the drought-sensitive check NIAB-IR-9 (Figure 2A). These results suggested that the GSR lines were comparatively less affected by drought stress and maintained the normal plant growth.
FIGURE 2. Effect of drought stress on (A) plant height/plant, (B) tillers/plant, (C) grain yield/plant, and (D) straw yield/plant. Values are means ± SD.
FIGURE 3. Effect of drought stress on (A) harvest index/plant, (B) total biomass/plant, (C) 1,000-grain weight, and (D) grain length. Values are means ± SD.
Overall, TPP were not significantly affected under drought stress except for Kashmir Basmati, while a few of the GSR lines showed increased TPP under drought (Figure 2B). GY is the most important agronomic trait of economic importance, and drought stress affected the GY in most genotypes except NGSR-15 and NGSR-18 (Figure 2C). Under drought stress, maximum GY was reported by NGSR-15 (90.7 g) and the minimum by NGSR-3 (48.6 g). Among the checks, IR-64 showed the maximum GY−1 (53.1 g), whereas the minimum was shown by the drought-sensitive NIAB-IR-9 (16.8 g). All the GSR lines (except NGSR-3) demonstrated higher grain yield than the check varieties; even the drought-sensitive NGSR-3 accounted for higher grain yield than the sensitive check NIAB-IR-9 (Figure 2C). These results suggested that the grain yield of GSR lines was less affected by drought stress compared with local checks.
SY was generally increased in most GSR accessions along with two check varieties IR-64 and Kissan Basmati under drought stress (Figure 2D). The maximum increase in SY was observed in NGSR-15, NGSR-6, NGSR18, and NGSR-12. However, NGSR-3, Kashmir Basmati, and NIAB-IR-9 showed a decrease in SY under drought. It is noteworthy that SY was only decreased in the most drought-sensitive GSR line and checks, thus, it is considered an important trait for drought escape at the flowering stage. This is because plants tend to continue their vegetative stage bypassing the flowering stage until they got favorable conditions.
Generally, drought stress negatively impacted the HI in all genotypes, but the non-GSR lines showed a higher decrease compared with the GSR lines with the highest decrease observed in our drought-sensitive check NIAB-IR-9 (Figure 3A). These results suggest that GSR lines have the potential to maintain the HI under drought stress conditions (Figure 2A).
TBM was not significantly affected under drought stress except for Kashmir Basmati and NIAB-IR-9 (Figure 3B). Similarly, TGW was not significantly affected under drought stress in the tested genotypes (Figure 3C).
Drought stress had a significant effect on GL; however, differences for genotypes were nonsignificant (Figure 3D). Three genotypes, including NGSR-3, NGSR-1, and NGSR-15, showed longer GLs. The genotypes NGSR-3 and NGSR-21 showed the highest (11.1 mm) and the lowest (7.3 mm) GL, respectively. Among the experimental checks, the maximum (10.9 mm) and the minimum (8.03 mm) GLs were depicted by Kissan Basmati and Kashmir Basmati, respectively. Overall, GSR lines maintained the grain length under drought stress compared with sensitive checks (Figure 3D), except NGSR-21, which showed a reduced grain length (7.3 mm).
Drought susceptibility index
DSI is an important indicator of drought tolerance, and a lower value indicates better tolerance. Overall, the GSR lines showed lower DSI compared with local varieties, and the genotypes NGSR-18, NGSR-15, and NGSR-16 exhibited the lowest DSIs among 22 GSR lines, showing their potential toward drought tolerance (Figure 4). In contrast, the highest DSI was recorded for NIAB-IR-9 followed by Kashmir Basmati, indicating the least drought tolerance among the tested genotypes. Overall, these results demonstrated that GSR lines generally performed better than the checks.
FIGURE 4. Frequency distribution of drought susceptibility index for grain yield showing the degree of susceptibility to drought stress. Genotypes below the line are declared as most drought-tolerant genotypes.
Pollen fertility percentage
Pollen fertility is an important indicator of drought tolerance as it directly affects the seed setting and ultimately the grain yield. The microscopic analyses of potassium iodide (I2-KI)-stained anthers revealed significant differences in PFPs between tolerant and sensitive genotypes. Overall, the GSR lines maintained higher PFPs under drought stress compared with non-GSR checks (Figure 5). While most of the GSR lines showed completely fertile pollens under drought, a higher sterility up to 57.1% was recorded in NGSR-3 (Figures 5 and 6A). In contrast, check varieties except the Kissan Basmati showed lower PFP compared with GSR under drought where NIAB-IR-9 showed only a 3.6% PFP (Figure 6A). These findings suggest that PFP could be a good indicator for drought tolerance in rice.
FIGURE 5. Examination of pollen fertility of the 22 GSR lines and 4 checks with I2-KI solution staining of the mature pollen grains. The sterile pollen grains failed to be stained or stained weakly, indicating that they did not contain starch or contained irregularly distributed starch, whereas the viable pollen grains were stained deep brown. Scale bars are 100 µm.
FIGURE 6. Effect of drought stress on (A) pollen fertility percentage, (B) normalized difference vegetative index (NDVI), and (C) cell membrane stability (CMS). Values (except CMS) are means ± SD.
Normalized difference vegetative index
NDVI has recently emerged as an indicator of plant health. We observed a considerable decrease in the value of NDVI under drought stress for GSR and check varieties (Figure 6B). The genotypes NGSR-13, IR-64, and Kissan Basmati showed the highest NDVIs (>0.7) under drought. Among GSR lines, the maximum NDVI (0.7) was reported in NGSR-13, whereas the minimum (0.59) was reported in NGSR-1. Similarly, among the check varities, the highest NDVI was depicted by IR-64 (0.72), whereas it was minimum by Kashmir Basmati (0.6) (Figure 6B). Since drought often causes leaf yellowing in plants, the reduced NDVI values under drought could be associated with yellow leaves.
Cell membrane stability
CMS indicates the stress tolerance ability of plant cells. Again, GSR lines showed higher CMS% than non-GSR, where NGSR-15 showed the highest (81.1%) CMS followed by NGSR-18 (Figure 6C), while the lowest was measured in NIAB-IR-9 (29.5%).
Correlation of grain yield with other agronomic traits
Understanding the correlation of grain yield with other agronomic and physiological traits is of prime importance as it helps to identify certain prebreeding traits that could be best indicators of grain yield. Under WW environment, GY has shown a significant positive correlation with SY (r = 0.59**), TBM (r = 0.90**), and HI (r = 0.60**) (Figure 7). In addition, a significant negative correlation was found for TPP with HI (r = −0.57**), which suggests the importance of optimum number of tillers for better HI and GY.
FIGURE 7. The scatter matrix below the histogram and correlation coefficient value with p-value above the histogram calculated from the means of all the studied traits under well-watered (WW) environmental condition. The p-values of all correlations were 0.05* and 0.01**.
Under drought stress, GY has shown a significant positive correlation with HI (r = 0.89***), CMS (r = 0.88***), TBM (r = 0.85***), PFP (r = 0.80***), and SY (r = 0.79***), while a significant negative correlation was found with DSI (r = −0.72***) and GL (r = −0.43*) (Figure 8). In addition to GY, PFP has shown a significant positive correlation with HI (r = 0.82**), CMS (r = 0.77**), SY (r = 0.71**), and TBM (r = 0.68**), while a significant negative correlation of PFP was observed with DSI (r = −0.69**). Notably, DSI had significant negative correlations with GY (r = −0.72**), HI (r = −0.72**), CMS (r = −0.70**), SY (r = −0.69**), PFP (r = −0.69**), TBM (r = −0.63**), and PH (r = −0.50**), which suggest the importance of DSI being an important indicator of drought susceptibility in rice (Figure 7). These findings revealed important agronomic and physiological traits to be considered as reliable selection criteria for screening rice germplasm against drought tolerance.
FIGURE 8. The scatter matrix below the histogram and correlation coefficient value with p-value above the histogram calculated from the means of all the studied traits under drought stress condition. The p-values of all correlations were 0.05* and 0.01**.
Expression analysis of drought-related genes
To see the role of drought responsive genes in drought tolerance of GSR, we analyzed the expression pattern in selected drought-tolerant and drought-sensitive genotypes using quantitative real-time PCR (Figure 9). The genotype NGSR-15 was selected as drought tolerant, and NGSR-3 and NIAB-IR-9 were chosen as drought-susceptible genotypes (Figure 10). Ten differentially expressed genes (DEGs) were selected for qRT-PCR analysis from the comparative transcriptome dataset between drought-sensitive (HHZ) and -tolerant (H471) genotypes (Huang et al., 2014). These genes have not been studied previously for their role in drought tolerance, except the transcriptome analysis. In addition, we analyzed the expression of five previously characterized genes for drought tolerance in rice (OsSADRI, OsDSM1, OsDT11, OsDREB1E, and OsDREB2B).
FIGURE 9. (qRT)-PCR analysis of drought-responsive genes in NGSR-3, NGSR-15, and NIAB-IR-9 revealed the relative expression in terms of fold change (log2FC). Young panicle tissues (∼1.5 cm) of three selected genotypes were employed in this analysis. Rice actin gene (OsACT1) was the internal control gene. Values of three biological replicates (n = 3) were expressed as mean ± SD.
FIGURE 10. Phenotypic comparison of plants, grain length, and shape (husked and de-husked), pollen viability, and panicle fertility of NGSR-15, NGSR-3, and NIAB-IR-9 under well-watered (WW) and drought stress. Abbreviations: F, fertile spikelets; S, sterile spikelets. White arrow heads denote fertile spikelets and yellow arrow heads denote sterile spikelets. Spikelets with open tips represent sterile spikelets with no seed set. Scale bars are approximately 15 cm (A,G,M), 5 mm (B,C,H,I,N,O), 5 cm (F,L,R), and 50 µm (D,E,J,K,P,Q).
Our results indicated that four genes, namely, LOC_Os11g36190 (receptor kinase), LOC_Os12g04500 (response regulator receiver domain-containing protein), LOC_Os12g26290 (alpha-DOX2), and LOC_Os02g11960 (ABC transporter, ATP-binding protein), were upregulated in drought-tolerant genotypes (NGSR-15) and downregulated in drought-sensitive genotypes (NGSR-3 and NIAB-IR-9) (Figure 9). This suggest that these genes may positively regulate drought tolerance in rice. Three genes, namely, LOC_Os05g23880 (lipoxygenase), LOC_Os05g08480 (cytokinin-O-glucosyltransferase 1), and LOC_Os01g28030 (peroxidase precursor), were downregulated in NGSR-15, while they were upregulated in drought-sensitive genotypes (NGSR-3 and NIAB-IR-9), suggesting a negative regulation for drought tolerance (Figure 9). Two genes, LOC_Os07g43560 (protein kinase-like) and LOC_Os07g48450 (no apical meristem protein), were downregulated under drought stress in all genotypes and, thus, may not be related to drought tolerance (Figure 9). The gene LOC_Os12g43450 (thaumatin family domain-containing protein) was upregulated in NIAB-IR-9, while it was downregulated in both NGSR-3 and NGSR-15, suggesting that this gene might be related to nonGSR rice. These results revealed a differential expression pattern of genes among drought-tolerant and -sensitive genotypes and, thus, could be employed for molecular identification of drought-tolerant rice genotypes at large scale. Notably, we observed an increased expression of previously known drought tolerance-related genes (OsSADRI, OsDSM1, and OsDT11) in NGSR-15, while the opposite was observed for NGSR-3 and NIAB-IR-9, which clearly indicated the role of these genes in drought tolerance. In addition, OsDREB1E has shown a sharp increase in expression under drought stress in all the three genotypes, while OsDREB2B has shown a sharp decrease in expression under drought stress in NGSR-15 and NIAB-IR-9 (Figures 9N, O). This suggest that DREB genes are probably not directly involved in drought tolerance in these tested genotypes, and their expression is modulated by some unknown genetic factors.
Discussion
Rice being a prime diet of 50% of the global population and the staple food of many countries is an important grain crop. However, growing rice requires high delta of water where limited water conditions affect its growth and grain yield. Water stress at anthesis directly affects seed setting and grain filling and, ultimately, the grain yield (Shokat et al., 2020a; Shokat et al., 2020b). Green super rice (GSR) was developed by combining the best global germplasm and has the potential to maintain the optimum grain yield under different stress conditions (Jewel et al., 2019). Moreover, this germplasm has never been evaluated for pre-anthesis drought stress. In the current experiment, 22 GSR genotypes and four local lines of Pakistan were used to understand the mechanism of yield reduction at pre-anthesis stages of drought stress. This germplasm was characterized for different agrophysiological traits, and then the most diverse genotypes were further evaluated by novel drought-responsive genes. Yield-related traits are important indications of final grain yield (Zafar et al., 2020a; Waqas et al., 2021). Studies reported that plant genotypes that maintained higher plant biomass under drought stress conditions often maintain higher grain number or weight and ultimately the grain yield (Shokat et al., 2021). In the current study, we identified NGSR-15 as a drought-tolerant line as it maintained higher CMS, PFP, TBM, and GY. In contrast, NGSR-3 and NIAB-IR-9 were ranked as drought-sensitive lines since they showed significant reductions in GY probably due to reduced CMS, PFP, and TBM. A phenotypic presentation of the performance of these genotypes under drought stress is shown in Figure 10. Higher biomass is usually linked with higher photosynthetic rate of the genotypes (Morinaka et al., 2006). Our results indicate that biomass partitioning toward grain filling was limited due to flowering stage drought stress, and there could be a possibility that GSR-15 has maintained a higher grain yield due to better seed setting under moisture stress conditions. To explain the possible mechanism of higher and lower grain yield for the genotype GSR-15 and NIAB-IR-9, respectively, we associated yield data with the few parameters of physiology to understand the physiological basis of yield reduction at flowering stage drought stress.
DSI indicates the extent of susceptibility by drought stress in terms of economically important traits particularly the grain yield. In this study, genotypes NGSR-15 and NGSR-18 showed the lowest susceptibility with values of −0.04 and −0.2, respectively, whereas NIAB-IR-9 (check) showed the highest DSI value of 3.6 (Figure 4). Under drought stress, permeability of membranes and leakage of ions occur from the weak or unstable membranes (Bajji et al., 2002). Likewise, seed setting is dependent on the viability of pollen, while limited water availability at the flowering stage can cause pollen abortion in sensitive genotypes (Mehri et al., 2020). In contrast, plant genotypes that show better cell membrane stability (CMS) or pollen fertility could perform better under flowering stage drought stress. In the current experiment, better cell membrane stability and pollen fertile percentage (PFP) was exhibited by the genotype NGSR-15, while the lowest values were recorded for NGSR-3 and NIAB-IR-9 indicating the physiological basis of drought tolerance and drought susceptibility of these genotypes respectively. A correlation and PCA was drawn to test the significance of these parameters in relation to yield and yield-related traits, and we found a strong significant and positive correlation of CMS and PFP with grain yield (Figure 1B). In contrast, association of grain yield was significant but negatively associated with DSI (Figure 1B) indicating that these traits could be selected as prebreeding traits for flowering stage drought stress in rice. To understand the molecular basis of drought tolerance, these three genotypes were further tested through gene expression.
Stress conditions change the expressions of the stress-induced regulatory or effector genes, which are usually involved in the regulation of normal processes of the plants (Ouyang et al., 2010). We investigated different categories of DEGs, controlling drought tolerance and sensitivity by up-/downregulation of DEGs. Furthermore, this analysis relied on two GSR genotypes and one locally developed genotype, NIAB-IR-9, in order to provide an accurate estimate of expression by comparing GSRs with traditional cultivars under flowering stage drought stress. In agreement with published literature, our expression results suggest the involvement of DEGs in drought tolerance or sensitivity (Chen et al., 2009; Narsai et al., 2013; Zhang et al., 2015b). Leucine-rich RLKs, play a key role in the regulation of plant growth under various abiotic stresses, and gene LOC_Os11g36190 (a leucine-rich receptor-like kinase) is predicted to be upregulated for bacterial leaf blight in rice (Zhang et al., 2015a; Ahsan et al., 2019). LOC_Os12g04500 and LOC_Os12g26290 are also reported as the core of the jasmonic acid (JA) signaling pathway, and in the current experiment, expression of these two genes was increased significantly under prolonged drought period. Moreover, JA signaling genes are also reported to be involved under critical phases of drought stress (Du et al., 2007). We found that four genes, i.e., LOC_Os11g36190, LOC_Os12g04500, LOC_Os12g26290, and LOC_Os02g119600, were upregulated in drought-tolerant genotypes (NGSR-15) and downregulated in drought-sensitive genotypes (NGSR-3 and NIAB-IR-9) indicating their positive relationship with drought tolerance. Likewise, an increased expression of previously known drought tolerance-related genes (OsSADRI, OsDSM1, and OsDT11) (Ning et al., 2009; Li et al., 2017; Park et al., 2018) was observed in NGSR-15, while an opposite trend was observed for NGSR-3 and NIAB-IR-9. This change in expression in the tested genes could be due to a sequence variation in their promoter region or mutation in a major upstream regulator, which is currently unknown to us. Apart from gene expression, these genotypes also showed a contrast for PFP, CMS, and DSI along with clear differences in grain yield suggesting their role in terminal-stage drought tolerance. Drought-responsive element-binding proteins (DREBs) are known to play important roles in abiotic stresses especially drought (Chen et al., 2008). Interestingly, in our study, expression of the DREB gene was either increased in both tolerant and sensitive genotypes, or decreased under drought stress (Figures 9M, O). This suggests that DREB genes are probably not directly involved in drought tolerance in these tested genotypes, and their expression is modulated by some unknown genetic factors.
Conclusion
Through this study, we identified molecular and physiological basis of higher grain yield at the flowering stage drought stress and the role of novel drought-responsive genes in drought tolerance. Importantly, various morphophysiological traits (PFP,CMS, DSI, and HI) had strong association with drought-responsive genes, and ultimately, the grain yield indicating these parameters could be used as prebreeding traits for drought tolerance. Our results also indicate that genotype NGSR-15 was the most drought tolerant, while NGSR-3 and NIAB-IR-9 were the most sensitive genotypes. These genotypes can further be used to improve rice yield under drought stress; however, in-depth mechanism is required to confirm our findings.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.
Author Contributions
SZ, MN, and MK designed the study. HA, SZ, and MN planted the fields and collected the data. SZ, MN, and SS analyzed the data. HA, SZ, and MN wrote the manuscript. SS, SI, GA, and MK provided the technical assistance. SN, JX, and ZL developed the GSR material. MK supervised the whole study and provided the technical and financial support.
Funding
This work was supported by the Ministry of Food Security and Research of Pakistan via the PSDP project: Green Super Rice in Pakistan (contract: 743).
Conflict 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.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
The authors would like to thank the PSDP and Ministry of Food Security and Research, Government of Pakistan for the financial support through the project "Green Super Rice in Pakistan: a component of productivity enhancement of rice".
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2022.832542/full#supplementary-material
References
Ahmed, S., Rashid, M. A. R., Zafar, S. A., Azhar, M. T., Waqas, M., and Uzair, M. (2021). Genome-wide Investigation and Expression Analysis of APETALA-2 Transcription Factor Subfamily Reveals its Evolution, Expansion and Regulatory Role in Abiotic Stress Responses in Indica Rice (Oryza Sativa L. Ssp. Indica). Genomics 113 (1 Pt 2), 1029–1043. doi:10.1016/j.ygeno.2020.10.037
Ahsan, A., Monir, M., Meng, X., Rahaman, M., Chen, H., and Chen, M. (2019). Identification of Epistasis Loci Underlying rice Flowering Time by Controlling Population Stratification and Polygenic Effect. DNA Res. 26 (2), 119–130. doi:10.1093/dnares/dsy043
Bajji, M., Kinet, J.-M., and Lutts, S. (2002). The Use of the Electrolyte Leakage Method for Assessing Cell Membrane Stability as a Water Stress Tolerance Test in Durum Wheat. Plant Growth Regul. 36 (1), 61–70. doi:10.1023/A:1014732714549
Barker, R., Dawe, D., Tuong, T., Bhuiyan, S., and Guerra, L. (1999). The Outlook for Water Resources in the Year 2020: Challenges for Research on Water Management in rice Production. Southeast Asia 1, 1–5.
Blum, A., and Ebercon, A. (1981). Cell Membrane Stability as a Measure of Drought and Heat Tolerance in Wheat 1. Crop Sci. 21 (1), 43–47. doi:10.2135/cropsci1981.0011183x002100010013x
Chaturvedi, A. K., Bahuguna, R. N., Shah, D., Pal, M., and Jagadish, S. K. (2017). High Temperature Stress during Flowering and Grain Filling Offsets Beneficial Impact of Elevated CO 2 on Assimilate Partitioning and Sink-Strength in rice. Scientific Rep. 7 (1), 1–13. doi:10.1038/s41598-017-07464-6
Chen, J. Q., Meng, X. P., Zhang, Y., Xia, M., and Wang, X. P. (2008). Over-expression of OsDREB Genes lead to Enhanced Drought Tolerance in rice. Biotechnol. Lett. 30 (12), 2191–2198. doi:10.1007/s10529-008-9811-5
Chen, S., Li, X.-Q., Zhao, A., Wang, L., Li, X., Shi, Q., et al. (2009). Genes and Pathways Induced in Early Response to Defoliation in rice Seedlings. Curr. Issues Mol. Biol. 11 (2), 81. doi:10.21775/cimb.011.081
Daryanto, S., Wang, L., and Jacinthe, P.-A. (2016). Global Synthesis of Drought Effects on Maize and Wheat Production. PLOS ONE 11 (5), e0156362. doi:10.1371/journal.pone.0156362
Desa, U. (2015). United Nations Department of Economic and Social Affairs, Population Division, World Population Prospects: The 2015 Revision, Key Findings and advance Tables. New York: Online Edition UN DESA.
Do, P. T., Prudent, M., Sulpice, R., Causse, M., and Fernie, A. R. (2010). The Influence of Fruit Load on the Tomato Pericarp Metabolome in a Solanum Chmielewskii Introgression Line Population. Plant Physiol. 154 (3), 1128–1142. doi:10.1104/pp.110.163030
Du, L., Jiao, F., Chu, J., Jin, G., Chen, M., and Wu, P. (2007). The Two-Component Signal System in rice (Oryza Sativa L.): a Genome-wide Study of Cytokinin Signal Perception and Transduction. Genomics 89 (6), 697–707. doi:10.1016/j.ygeno.2007.02.001
Fang, J., Guo, T., Xie, Z., Chun, Y., Zhao, J., Peng, L., et al. (2021). The URL1-ROC5-TPL2 Transcriptional Repressor Complex Represses the ACL1 Gene to Modulate Leaf Rolling in rice. Plant Physiol. 185 (4), 1722–1744. doi:10.1093/plphys/kiaa121
Govaerts, B., and Verhulst, N. (2010). The Normalized Difference Vegetation index (NDVI) Greenseeker (TM) Handheld Sensor: Toward the Integrated Evaluation of Crop Management Part A: Concepts and Case Studies. Mexico: CIMMYT.
Huang, L., Zhang, F., Wang, W., Zhou, Y., Fu, B., and Li, Z. (2014). Comparative Transcriptome Sequencing of Tolerant rice Introgression Line and its Parents in Response to Drought Stress. BMC genomics 15 (1), 1–16. doi:10.1186/1471-2164-15-1026
Hussain, H. A., Men, S., Hussain, S., Chen, Y., Ali, S., Zhang, S., et al. (2019). Interactive Effects of Drought and Heat Stresses on Morpho-Physiological Attributes, Yield, Nutrient Uptake and Oxidative Status in maize Hybrids. Scientific Rep. 9 (1), 3890. doi:10.1038/s41598-019-40362-7
Irri, I. (2002). Standard Evaluation System for rice. Philippine: International Rice Research Institute, 1–45.
Islam, M. A., Du, H., Ning, J., Ye, H., and Xiong, L. (2009). Characterization of Glossy1-Homologous Genes in rice Involved in Leaf Wax Accumulation and Drought Resistance. Plant Mol. Biol. 70 (4), 443–456. doi:10.1007/s11103-009-9483-0
Jewel, Z. A., Ali, J., Pang, Y., Mahender, A., Acero, B., Hernandez, J., et al. (2019). Developing green Super rice Varieties with High Nutrient Use Efficiency by Phenotypic Selection under Varied Nutrient Conditions. Crop J. 7 (3), 368–377. doi:10.1016/j.cj.2019.01.002
Khanna-Chopra, R., and Viswanathan, C. (1999). Evaluation of Heat Stress Tolerance in Irrigated Environment of T. aestivum and Related Species. I. Stability in Yield and Yield Components. Euphytica 106 (2), 169–180. doi:10.1023/a:1003531722420
Kumar, A., Bernier, J., Verulkar, S., Lafitte, H., and Atlin, G. (2008). Breeding for Drought Tolerance: Direct Selection for Yield, Response to Selection and Use of Drought-Tolerant Donors in upland and lowland-adapted Populations. Field Crops Res. 107 (3), 221–231. doi:10.1016/j.fcr.2008.02.007
Li, X., Han, H., Chen, M., Yang, W., Liu, L., Li, N., et al. (2017). Overexpression of OsDT11, Which Encodes a Novel Cysteine-Rich Peptide, Enhances Drought Tolerance and Increases ABA Concentration in rice. Plant Mol. Biol. 93 (1), 21–34. doi:10.1007/s11103-016-0544-x
Liu, J., Liao, D., Oane, R., Estenor, L., Yang, X., Li, Z., et al. (2006). Genetic Variation in the Sensitivity of Anther Dehiscence to Drought Stress in rice. Field Crops Res. 97 (1), 87–100. doi:10.1016/j.fcr.2005.08.019
Mehri, N., Fotovat, R., Mirzaei, M., Fard, E. M., Parsamatin, P., Hasan, M. T., et al. (2020). Proteomic Analysis of Wheat Contrasting Genotypes Reveals the Interplay between Primary Metabolic and Regulatory Pathways in Anthers under Drought Stress. J. Proteomics 226, 103895. doi:10.1016/j.jprot.2020.103895
Mohammadi, S. A., and Prasanna, B. (2003). Analysis of Genetic Diversity in Crop Plants—Salient Statistical Tools and Considerations. Crop Sci. 43 (4), 1235–1248. doi:10.2135/cropsci2003.1235
Morinaka, Y., Sakamoto, T., Inukai, Y., Agetsuma, M., Kitano, H., Ashikari, M., et al. (2006). Morphological Alteration Caused by Brassinosteroid Insensitivity Increases the Biomass and Grain Production of Rice. Plant Physiol. 141 (3), 924–931. doi:10.1104/pp.106.077081
Nakashima, K., Tran, L. S. P., Van Nguyen, D., Fujita, M., Maruyama, K., Todaka, D., et al. (2007). Functional Analysis of a NAC‐type Transcription Factor OsNAC6 Involved in Abiotic and Biotic Stress‐responsive Gene Expression in rice. Plant J. 51 (4), 617–630. doi:10.1111/j.1365-313x.2007.03168.x
Narsai, R., Wang, C., Chen, J., Wu, J., Shou, H., and Whelan, J. (2013). Antagonistic, Overlapping and Distinct Responses to Biotic Stress in rice (Oryza Sativa) and Interactions with Abiotic Stress. BMC genomics 14 (1), 1–21. doi:10.1186/1471-2164-14-93
Ning, J., Li, X., Hicks, L. M., and Xiong, L. (2009). A Raf-like MAPKKK Gene DSM1 Mediates Drought Resistance through Reactive Oxygen Species Scavenging in Rice. Plant Physiol. 152 (2), 876–890. doi:10.1104/pp.109.149856
Oliva, R., Ji, C., Atienza-Grande, G., Huguet-Tapia, J. C., Perez-Quintero, A., Li, T., et al. (2019). Broad-spectrum Resistance to Bacterial Blight in rice Using Genome Editing. Nat. Biotechnol. 37 (11), 1344–1350. doi:10.1038/s41587-019-0267-z
Ouyang, S. Q., Liu, Y. F., Liu, P., Lei, G., He, S. J., Ma, B., et al. (2010). Receptor‐like Kinase OsSIK1 Improves Drought and Salt Stress Tolerance in rice (Oryza Sativa) Plants. Plant J. 62 (2), 316–329. doi:10.1111/j.1365-313x.2010.04146.x
Panda, D., Mishra, S. S., and Behera, P. K. (2021). Drought Tolerance in Rice: Focus on Recent Mechanisms and Approaches. Rice Sci. 28 (2), 119–132. doi:10.1016/j.rsci.2021.01.002
Park, Y. C., Chapagain, S., and Jang, C. S. (2018). A Negative Regulator in Response to Salinity in Rice: Oryza Sativa Salt-, ABA- and Drought-Induced RING Finger Protein 1 (OsSADR1). Plant Cel Physiol 59 (3), 575–589. doi:10.1093/pcp/pcy009
Prasad, P. V., Bheemanahalli, R., and Jagadish, S. K. (2017). Field Crops and the Fear of Heat Stress—Opportunities, Challenges and Future Directions. Field Crops Res. 200, 114–121. doi:10.1016/j.fcr.2016.09.024
Shao, G.-C., Deng, S., Liu, N., Yu, S.-E., Wang, M.-H., and She, D.-L. (2014). Effects of Controlled Irrigation and Drainage on Growth, Grain Yield and Water Use in Paddy rice. Eur. J. Agron. 53, 1–9. doi:10.1016/j.eja.2013.10.005
Shokat, S., Großkinsky, D. K., and Liu, F. (2021). Impact of Elevated CO2 on Two Contrasting Wheat Genotypes Exposed to Intermediate Drought Stress at Anthesis. J. Agron. Crop Sci. 207 (1), 20–33. doi:10.1111/jac.12442
Shokat, S., Großkinsky, D. K., Roitsch, T., and Liu, F. (2020a). Activities of Leaf and Spike Carbohydrate-Metabolic and Antioxidant Enzymes Are Linked with Yield Performance in Three spring Wheat Genotypes Grown under Well-Watered and Drought Conditions. BMC Plant Biol. 20 (1), 400. doi:10.1186/s12870-020-02581-3
Shokat, S., Sehgal, D., Vikram, P., Liu, F., and Singh, S. (2020b). Molecular Markers Associated with Agro-Physiological Traits under Terminal Drought Conditions in Bread Wheat. Int. J. Mol. Sci. 21 (9), 3156. doi:10.3390/ijms21093156
Shuxing, L. (2014). The Responding of Rice after Water Stress in Young Panicle Formation Stage. ACTA Agriculturae Boreali-Sinica.
Tong, H., Jin, Y., Liu, W., Li, F., Fang, J., Yin, Y., et al. (2009). Dwarf and Low‐Tillering, A New Member of the GRAS Family, Plays Positive Roles in Brassinosteroid Signaling in rice. Plant J. 58 (5), 803–816. doi:10.1111/j.1365-313x.2009.03825.x
Tripathy, J. N., Zhang, J., Robin, S., Nguyen, T. T., and Nguyen, H. T. (2000). QTLs for Cell-Membrane Stability Mapped in rice (Oryza Sativa L.) under Drought Stress. Theor. Appl. Genet. 100 (8), 1197–1202. doi:10.1007/s001220051424
Waqas, M. A., Wang, X., Zafar, S. A., Noor, M. A., Hussain, H. A., Azher Nawaz, M., et al. (2021). Thermal Stresses in Maize: Effects and Management Strategies. Plants (Basel) 10 (2). doi:10.3390/plants10020293
Wattoo, F. M., Rana, R. M., Fiaz, S., Zafar, S. A., Noor, M. A., Hassan, H. M., et al. (2018). Identification of Drought Tolerant maize Genotypes and Seedling Based Morpho-Physiological Selection Indices for Crop Improvement. Sains Malays 47, 295–302.
Wing, R. A., Purugganan, M. D., and Zhang, Q. (2018). The rice Genome Revolution: from an Ancient Grain to Green Super Rice. Nat. Rev. Genet. 19 (8), 505–517. doi:10.1038/s41576-018-0024-z
Yu, S., Ali, J., Zhang, C., Li, Z., and Zhang, Q. (2020). Genomic Breeding of Green Super Rice Varieties and Their Deployment in Asia and Africa. Theor. Appl. Genet. 133 (5), 1427–1442. doi:10.1007/s00122-019-03516-9
Yugandhar, P., Basava, R. K., Desiraju, S., Voleti, S. R., Sharma, R., and Neelamraju, S. (2017). Identifying Markers Associated with Yield Traits in Nagina22 rice Mutants Grown in Low Phosphorus Field or in Alternate Wet/dry Conditions. Aust. J. Crop Sci. 11 (5), 548–556. doi:10.21475/ajcs.17.11.05.p372
Zafar, S. A., Hameed, A., Ashraf, M., Khan, A. S., Qamar, Z. U., Li, X., et al. (2020a). Agronomic, Physiological and Molecular Characterisation of rice Mutants Revealed the Key Role of Reactive Oxygen Species and Catalase in High-Temperature Stress Tolerance. Funct. Plant Biol. 47 (5), 440–453. doi:10.1071/fp19246
Zafar, S. A., Hameed, A., Khan, A. S., and Ashraf, M. (2017). Heat Shock Induced Morpho-Physiological Response in Indica rice (Oryza Sativa L.) at Early Seedling Stage. Pak. J. Bot. 49 (2), 453–463.
Zafar, S. A., Hameed, A., Nawaz, M. A., Ma, W., Noor, M. A., Hussain, M., et al. (2018). Mechanisms and Molecular Approaches for Heat Tolerance in rice (Oryza Sativa L.) under Climate Change Scenario. J. Integr. Agric. 17 (4), 726–738. doi:10.1016/S2095-3119(17)61718-0
Zafar, S. A., Patil, S. B., Uzair, M., Fang, J., Zhao, J., Guo, T., et al. (2020b). Degenerated Panicle and Partial Sterility 1 (Dps1) Encodes a Cystathionine β-synthase Domain Containing Protein Required for Anther Cuticle and Panicle Development in rice. New Phytol. 225 (1), 356–375. doi:10.1111/nph.16133
Zafar, S. A., Uzair, M., Khan, M. R., Patil, S. B., Fang, J., Zhao, J., et al. (2021). DPS1 Regulates Cuticle Development and Leaf Senescence in rice. Food and Energy Security 10 (1), e273. doi:10.1002/fes3.273
Zafar, S. A., Zaidi, S. S., Gaba, Y., Singla-Pareek, S. L., Dhankher, O. P., Li, X., et al. (2020c). Engineering Abiotic Stress Tolerance via CRISPR/Cas-Mediated Genome Editing. J. Exp. Bot. 71 (2), 470–479. doi:10.1093/jxb/erz476
Zhang, F., Huang, L.-Y., Ali, J., Cruz, C. V., Zhuo, D.-L., Du, Z.-L., et al. (2015a). Comparative Transcriptome Profiling of a rice Line Carrying Xa39 and its Parents Triggered by Xanthomonas Oryzae Pv. Oryzae Provides Novel Insights into the Broad-Spectrum Hypersensitive Response. BMC genomics 16 (1), 1–14. doi:10.1186/s12864-015-1329-3
Zhang, F., Huang, L., Zhang, F., Hu, D., Wu, W., Wang, W., et al. (2015b). Interacting Transcriptomes Revealing Molecular Mechanisms Underlying Xa39 Mediated Broad Spectrum Resistance of rice to Bacterial Blight. The Plant Genome 8 (3), 0094. plantgenome2014.2012. doi:10.3835/plantgenome2014.12.0094
Zhang, J., Zhang, S., Cheng, M., Jiang, H., Zhang, X., Peng, C., et al. (2018). Effect of Drought on Agronomic Traits of Rice and Wheat: A Meta-Analysis. Int. J. Environ. Res. Public Health 15 (5), 839. doi:10.3390/ijerph15050839
Zhang, Q. (2007). Strategies for Developing Green Super Rice. Proc. Natl. Acad. Sci. U S A. 104 (42), 16402–16409. doi:10.1073/pnas.0708013104
Keywords: drought, anthesis, pollen fertility, grain yield, correlation, drought-responsive genes
Citation: Ahmad H, Zafar SA, Naeem MK, Shokat S, Inam S, Rehman MAu, Naveed SA, Xu J, Li Z, Ali GM and Khan MR (2022) Impact of Pre-Anthesis Drought Stress on Physiology, Yield-Related Traits, and Drought-Responsive Genes in Green Super Rice. Front. Genet. 13:832542. doi: 10.3389/fgene.2022.832542
Received: 10 December 2021; Accepted: 27 January 2022;
Published: 24 March 2022.
Edited by:
Alvina Gul, National University of Sciences and Technology (NUST), PakistanReviewed by:
Muhammad Uzair, the University of Tennessee, Knoxville, United StatesSuyash Patil, International Rice Research Institute, IRRI, India
Copyright © 2022 Ahmad, Zafar, Naeem, Shokat, Inam, Rehman, Naveed, Xu, Li, Ali and Khan. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Muhammad Ramzan Khan, drmrkhan_nigab@yahoo.com
†Present address: Syed Adeel Zafar, Department of Botany and Plant Sciences, University of California, Riverside, Riverside, CA, United States
‡These authors share first authorship