Figures
Abstract
Background
Rice (Oryza sativa L.) is a highly drought sensitive crop, and most semi dwarf rice varieties suffer severe yield losses from reproductive stage drought stress. The genetic complexity of drought tolerance has deterred the identification of agronomically relevant quantitative trait loci (QTL) that can be deployed to improve rice yield under drought in rice. Convergent evidence from physiological characterization, genetic mapping, and multi-location field evaluation was used to address this challenge.
Methodology/Principal Findings
Two pairs of backcross inbred lines (BILs) from a cross between drought-tolerant donor Aday Sel and high-yielding but drought-susceptible rice variety IR64 were produced. From six BC4F3 mapping populations produced by crossing the +QTL BILs with the −QTL BILs and IR64, four major-effect QTL - one each on chromosomes 2, 4, 9, and 10 - were identified. Meta-analysis of transcriptome data from the +QTL/−QTL BILs identified differentially expressed genes (DEGs) significantly associated with QTL on chromosomes 2, 4, 9, and 10. Physiological characterization of BILs showed increased water uptake ability under drought. The enrichment of DEGs associated with root traits points to differential regulation of root development and function as contributing to drought tolerance in these BILs. BC4F3-derived lines with the QTL conferred yield advantages of 528 to 1875 kg ha−1 over IR64 under reproductive-stage drought stress in the targeted ecosystems of South Asia.
Conclusions/Significance
Given the importance of rice in daily food consumption and the popularity of IR64, the BC4F3 lines with multiple QTL could provide higher livelihood security to farmers in drought-prone environments. Candidate genes were shortlisted for further characterization to confirm their role in drought tolerance. Differential yield advantages of different combinations of the four QTL reported here indicate that future research should include optimizing QTL combinations in different genetic backgrounds to maximize yield advantage under drought.
Citation: Swamy B. P. M, Ahmed HU, Henry A, Mauleon R, Dixit S, Vikram P, et al. (2013) Genetic, Physiological, and Gene Expression Analyses Reveal That Multiple QTL Enhance Yield of Rice Mega-Variety IR64 under Drought. PLoS ONE 8(5): e62795. https://doi.org/10.1371/journal.pone.0062795
Editor: Turgay Unver, Cankiri Karatekin University, Turkey
Received: October 3, 2012; Accepted: March 25, 2013; Published: May 8, 2013
Copyright: © 2013 Swamy B. et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: The authors thank the Bill & Melinda Gates Foundation (BMGF), USA, and Generation Challenge Program (GCP), Mexico, for providing financial support to this study. The funders had no role in study design, data collection, and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Among cereals, rice (Oryza sativa L.) is the most drought-sensitive crop. Even a mild drought stress during the reproductive stage results in severe yield losses [1]–[3]. Most of the semi-dwarf high-yielding varieties developed during the green revolution era were meant for irrigated ecosystems and are highly susceptible to drought [4]. Since high-yielding drought-tolerant cultivars are not available, farmers in drought-prone areas cultivate either high-yielding cultivars with good grain quality that are drought susceptible or low-yielding traditional cultivars that are drought tolerant but have poor grain quality and also less input-use efficiency [5]–[7].
An understanding of the sources of genetic variation and physiological mechanisms involved facilitates the development of an appropriate strategy to breed drought-tolerant cultivars [8], [9]. Deep root growth, which may increase water uptake during progressive soil drying, is suggested to be a likely mechanism to confer increased yield under drought. However, there is little direct evidence in the literature of deep root growth conferring a yield advantage under drought [10]. A drought-yield effect of QTLs for deep roots and improved soil penetration [11]–[14] is yet to be confirmed.
Recent studies have identified QTL for yield under drought in rice [15]–[18]. Some of these QTL were derived from traditional donors and carry linkages for undesirable traits along with an effect on grain yield under drought [18]. The advanced backcross QTL (AB-QTL) approach involves two or more backcrosses to the improved recurrent parent to simultaneously identify and introgress QTL in the recurrent parent and to reduce undesirable linkages [19], [20]. AB-QTL analysis on lines with similar agro-morphological characters also provides the opportunity to impose uniform drought stress on all lines and to control differences due to phenology, leading to the detection of more reliable QTL. However, the genetic mapping of complex traits from parents with similar genetic backgrounds is difficult due to low polymorphism.
Expression profiling of contrasting parents under drought stress helps to identify differentially expressed genes and their regions in the genome [21]. The regions enriched with differentially expressed genes can be further genotyped with polymorphic molecular markers to detect the loci for complex traits. The differential expression patterns of drought-responsive genes in different plant tissues at different growth stages could provide an opportunity to characterize the traits associated with yield advantage under drought and to understand the physiological and molecular mechanisms that confer increased drought tolerance.
In this study, major QTL for grain yield under drought were narrowly delimited by expression polymorphism, and then identified in multiple mapping populations by genotyping and phenotyping under managed drought stress. We report physiological differences in backcross inbred lines (BILs) that were genetically similar but showed contrasting responses in yield under drought. The study identified lines with different combinations of QTL in the IR64 background that showed enhanced grain yield under drought in multi-location evaluations in the target environment, thereby confirming the value of these QTL for sustainable yield under drought stress.
Results
Four QTL for Grain Yield under Drought Identified
To define the QTL regions responsible for improved grain yield under drought in BILs derived from and IR64×Aday Sel cross [22] (Table S1), we used Affymetrix Rice Chip analysis to identify genome polymorphism. This approach was chosen after attempts to characterize the QTL regions with SSR markers did not reveal sufficient polymorphism between the parents. Four polymorphic regions were found at 6.8–7.3, 6.7–7.2, 14.6–16.5, and 18.6–19.3 Mb on chromosomes 2, 4, 9, and 10, respectively. In total, 5, 3, 8, and 5 polymorphic SSR markers in the regions detected by the chip-based analysis on chromosomes 2, 4, 9, and 10, respectively, were run on the whole population to detect QTL for grain yield (GY) (Table 1) and related traits (days to 50% flowering, DTF; plant height, PH; Table S2). The mapping results validated the regions predicted by the Affymetrix Rice Chip analysis. Four QTL were identified for GY under drought stress; qDTY2.2, qDTY4.1, qDTY9.1, and qDTY10.1, by both Interval and Composite Interval mapping methods (Table 1 and Fig. S1). The phenotypic variance explained by the QTL ranged from 6 to 19.2% (Table 1). None of these four qDTY QTL was associated with GY under non-stress conditions. QTL analyses in connected populations also consistently detected the major effect QTL qDTY2.2 and qDTY9.1 (Table S3 and Fig. S2).
Physiological Characterization of BILs Contrasting for Yield under Drought
Two pairs of genetically similar BILs (+QTL and –QTL lines from an IR64 × Aday Sel cross; pair 1: +QTL IR77298-5-6-B-18 and −QTL IR77298-5-6-B-11, and pair 2: +QTL IR77298-14-1-2-B-10 and –QTL IR77298-14-1-2-B-13) that were contrasting for yield under drought were characterized in the field to study the physiological mechanisms associated with increased yield under drought. +QTL lines showed cooler canopy temperature and greater stomatal conductance than −QTL lines and IR64 under the most severe drought stress, but not under mild drought stress or non-stress conditions (Fig. 1, Fig. S3). BILs within each pair did not differ significantly in shoot mass or NDVI (Normalized Difference Vegetation Index; Fig. 2), but shoot mass was greater towards the end of the 2010DS and NDVI was consistently greater in 2010DS-ROS in one pair of BILs (+QTL IR77298-14-1-2-B-10 and –QTL IR77298-14-1-2-B-13) compared to the other pair. Root growth at depth in terms of root length density was not greater in +QTL lines in any experiment (Table 2).
Significant differences among lines are indicated by *(p<0.05) and **(p<0.01). D. Infrared thermal image showing contrasting canopy temperatures of +QTL and –QTL lines, taken 98 days after sowing (DAS) in 2010DS-ROS.
Arrows indicate dates that the experiment was re-watered. No significant differences were observed between lines within a BIL pair on any sampling date.
Meta-analysis of Differentially Expressed Genes Relative to the Entire Genome and within the Drought Yield QTL Regions
Transcriptome data of the same parental +QTL and –QTL BILs under two water stress treatments (0.5 FTSW, 0.2 FTSW) from a previous study by Moumeni et al. [21] were re-analyzed to determine their association with QTL detected in this study. For the root transcriptome, highly contrasting counts of DEGs were detected from comparison of +QTL and –QTL lines between the two BIL pairs used, with 570 DEGs in pair 1 and 2,127 DEGs in pair 2, an almost fourfold difference in DEG counts between the two BIL pairs. For the leaf transcriptome, a similar number of DEGs was detected from comparison of +QTL and –QTL lines between the two BIL pairs (748 DEGs in pair 1; 779 in pair 2). For the panicle transcriptome, both +QTL and –QTL BIL pairs had similar but low DEG counts (240 DEGs in pair 1; 201 in pair 2).
For the first meta-analysis, the number of DEGs was counted within 1 MB genome blocks, with sliding window blocks of 500 kb, for the entire genome, to determine whether DEGs were aggregating (having a significantly higher number of DEGs than anywhere else in the genome at p<0.01) in blocks of the genome, adapting the genomic method of Bruce et al. [23] for gene expression data. For pair 1, aggregation analysis of the root, leaf, and panicle transcriptomes all pointed to the five genome regions in chromosomes 5, 9, 10, and 12. For pair 2, transcriptomes from leaf and panicle tissues showed five overlapping regions of DEG aggregates in chromosomes 2, 8, and 11, whereas, for the root transcriptome, only the chromosome 8 DEG aggregation region overlapped with the regions from leaf and panicle tissues, and a unique aggregation region was found in chromosome 5. A total of 5 and 9 distinct DEG aggregation regions were determined for pairs 1 and 2, respectively, and no common regions exist between the two BIL pairs (Table S4).
The second meta-analysis used QTL region information from this study. Using all the DEGs from the significance (re)analysis, we found that there was significant association of 96 DEGs in QTL on chromosomes 2, 4, 9, and 10 (p<0.05, Table 3), with 91 of these DEGs located within the DEG aggregation regions found in chromosomes 2, 9, and 10. Further analysis for enrichment of biological themes/categories in this subset of DEGs, however, showed no explicitly drought-responsive categories being identified using the GO-SLIM and Mapman gene classification (Table 4). However, interesting significant associations of these DEGs with previously identified QTL (p<0.05, using the Gramene QTL categories) were found, such as QTL for root length, deep root dry weight, root weight, penetrated root thickness, and root penetration index, as well as QTL associated with yield traits (spikelet density, spikelet fertility, seed weight, and other grain/seed-related QTL; Table S5).
Improved Lines with QTL Introgressions Identified
From the six mapping populations, 4 lines with four QTL combinations, 15 lines with three QTL combinations, 29 lines with two QTL combinations, and 19 lines with different single QTL were identified. These lines were further selected for yield under drought and under well-watered conditions, high phenotypic and genetic similarity to IR64, and grain quality traits similar to those of IR64 (Table 5, Table S6). Before advancing to evaluation in the target environment, the IR64 QTL lines were screened under managed drought stress in large plots at IRRI, in which all lines showed a yield advantage of 194 to 1920 kg ha−1 over IR64. Subsequently, the three most promising lines, IR87707-445-B-B, IR87707-446-B-B, and IR87707-182-B-B, were evaluated in target drought-prone ecosystems in Bangladesh, India, and Nepal. The three lines showed yield advantages of 528 to 1875 kg ha−1 over IR64 under drought, and produced either similar or higher yields than IR64 under well-watered conditions (Table 6).
Discussion
Crop genetic improvement for drought is challenging because of its complex genetic nature and poor understanding of the physiological and molecular mechanisms associated with drought tolerance [8], [9]. We applied multiple approaches including genetic mapping, physiological characterization, and expression analyses to identify major-effect drought grain yield QTL, and successfully deployed them to improve grain yield under drought in the background of rice mega-variety IR64.
QTL, when mapped back to the physical genome, often translate to tens of megabases in size, which is not precise for identification of gene(s) underlying the QTL effect. We used a meta-analysis approach that combines results from the analysis of transcriptome data (aggregation of DEGs within QTL, and gene set enrichment tests) and QTL position in the physical genome to overcome low resolution in genetic mapping and enable us to identify a smaller set of candidate genes within the QTL region.
The consistent effect of four drought-yield QTL in the background of popular variety IR64 indicates their suitability for marker-assisted breeding (MAB) to improve the drought tolerance of IR64. The region on chromosome 2 in which qDTY2.2 was detected has been reported previously to have an effect on drought-related traits other than GY, including drought tolerance index, canopy temperature, osmotic adjustment, and leaf water content [24]–[27]. The clustering of a number of physiological traits in the same region suggests its importance for GY under severe drought stress. Differential expression of qDTY4.1 was observed under different severities of stress, indicating its usefulness over a wide range of stress severities. Lanceras et al. [28] also reported that different QTL were detected on different chromosomes depending on the severity of stress, and that three different grain yield QTL on chromosome 4 were identified in the CT9993/IR62266 doubled-haploid (DH) population under non-stress, mild, and severe drought stress conditions. QTL for GY and for grains per panicle were also detected on chromosome 4 under drought stress conditions in the CT9993/IR62266 DH population at the distal end of chromosome 4 [29] and for spikelet fertility and for grain weight in the Bala/Azucena DH population [30]. qDTY9.1 showed high phenotypic variance in both years, indicating that this QTL was equally effective under both severe and moderate drought stress. This region on chromosome 9 was reported to have a QTL for spikelet fertility (27), GY, and plant height under non-stress conditions [31], and for grain weight under non-stress conditions [32]–[34]. In joint analyses of connected populations, the major effect QTL qDTY2.2 and qDTY9.1 were also detected, indicating their reliability and worthiness for use in MAS. qDTY10.1 with high additive effect and a high phenotypic variance of 17% was identified under severe drought stress in the population P3 in 2010DS. A QTL for GY was also reported in a nearby region on chromosome 10 under severe drought stress [25], [28], [29] and under non-stress conditions [35]. A QTL for grain weight was also found near qDTY10.1 [36] under stress conditions. In all of these studies, including the present study, the yield-increasing allele was contributed by the drought-tolerant parent.
The physiological characterization of the BILs showed that +QTL lines maintained higher transpiration rates under drought, as evidenced by cooler canopy temperature and higher stomatal conductance than −QTL lines as soil water availability decreased (Fig. 1). We expected that this transpiration advantage in +QTL lines would be conferred by greater root growth at depth, but this was not observed. The similar shoot and root growth among +QTL and −QTL lines within BIL pairs suggests that the yield advantages under drought are not due to architectural or allometric differences in plant growth. Therefore, these results point to differences in water uptake among QTL lines; the QTL effect in the BILs appears to be associated with root anatomy/development or regulation of root function rather than deep root growth.
Gene expression in leaf, panicle, and root tissues under different drought stresses also pointed to differences in root development between the +QTL and −QTL BILs. Aggregation analysis of transcriptome data already available from genome-wide comparisons [21] allowed us to identify a subset of DEGs that aggregated in 14 blocks in the genome, which we then overlayed with QTL information independently generated in the present study, as well as from other published studies. The co-localization of DEGs in this study with published QTL suggests that some of the DEGs could be related to root structure and yield parameters. The convergent analysis of QTL, expression, and previously published data has enabled us to nominate a relatively small number of candidate genes for further investigation (Table S7). A caveat is that this analysis relied on a japonica reference genome, which may not contain the actual genes that are related to the QTL found. The availability of an IR64 reference genome and genome sequences for the +QTL BILs will provide more accurate positions of QTL and DEGs and is likely to improve the power of QTL enrichment analysis.
A literature survey was carried out to know the specific/broad biological functions of the differentially expressed genes/gene families between the +QTL and –QTL lines. In terms of gene function, most of the DEGs were associated with growth and development processes, stress responses and signaling, and hormonal regulation (Table S8). The growth- and development-related genes/gene families were associated with meristematic tissue activity in the roots, inflorescence, and leaf; genes for cell wall synthesis; root elongation; formation of lateral roots; pollen fertility; embryogenesis; seed development; and seed germination. The genes for hormonal regulation were associated with auxin and cytokinin homeostasis, the two most important hormones involved in growth and development processes. Many stress-signaling and stress-responsive genes were differentially expressed between the +QTL and −QTL lines, and an interesting abundance of four transcription factor genes all on chromosome 9 was observed (Os09 g24660: zinc finger motif, C2HC5-type family protein, putative, expressed; Os09 g24800: myb-related protein Myb4, putative; Os09 g25420: zinc finger, C2H2-type family protein, expressed; and Os09 g26180: transcription initiation factor TFIID subunit 10, putative, expressed). Although these expressed genes show some consistencies with the physiological responses observed at the whole-plant level in the field, further research is necessary to dissect the water uptake attributes associated with low canopy temperature under drought in these lines, as well as to identify the key genes within the QTL regions associated with drought response.
The higher yield of IR64+QTL NILs compared with that of IR64 under drought stress in two seasons indicates the consistent effect of the QTL. Among the lines with individual QTL, qDTY4.1 lines performed better than lines with other QTL. Three lines with qDTY2.2, qDTY4.1, and qDTY10.1, and two lines with qDTY2.2 and qDTY4.1 QTL combinations performed better than lines with other QTL combinations, indicating a differential synergistic relationship between QTL combinations [37]. A better understanding of the complementary effects among the four QTL could enable pyramiding of QTL in appropriate combinations to maximize the yield advantage under drought. Such optimization of QTL combinations may be required for introgression into different genetic backgrounds. The best-performing IR64 QTL lines at IRRI were also the highest performing lines in Bangladesh, India, and Nepal, indicating stable performance of the IR64 QTL lines across environments and that results from managed dry-season drought screening at IRRI are relevant to the target drought-prone environments.
In this study, IR64 QTL lines were developed that confer a yield advantage of 500 to 1800 kg ha−1 under different severities of drought over IR64, a popular rice mega-variety grown on millions of hectares in Asia. These lines can be immediately disseminated for cultivation by farmers in drought-prone environments. The identification of multiple QTL that collectively enhance yield under drought in rice confirms the effectiveness of direct selection for increased yield under drought. Although canopy temperature and transpiration results indicated increased water uptake under drought, we did not see differences in root length density in the BILs, suggesting a role of the QTL in root function rather than root architecture. Hypotheses on the functional roles of these QTL were provided by convergent analysis of mapping and gene expression patterns, with candidate genes significantly enriched in the mapped QTL regions being nominated for further characterization. Differential yield advantages of different combinations of the four QTL reported here indicate that future research should include optimizing QTL combinations in different genetic backgrounds to maximize yield advantage under drought.
Materials and Methods
Plant Materials and Mapping Populations for Physiological and Genetic Study
Two BILs, IR77298-14-1-2 (BC3F3∶4) and IR77298-5-6-B-18 (BC3F2∶3), were derived from a cross between drought-tolerant traditional donor Aday Sel and drought-susceptible recurrent parent IR64. These backcross lines had been fixed for most of the agro-morphological traits such as days to 50% flowering (DTF), plant height (PH), grain yield (GY), and biomass (BIO), but showed differential performance under drought [22]. A set of around 60 single panicles from the two lines with differential performance under drought was selected and evaluated under drought stress in lowland conditions. For physiological characterization, we used two pairs of +QTL and −QTL BILs (IR77298-5-6-B-18 (+QTL) and IR77298-5-6-B-11 (−QTL) and IR77298-14-1-2-B-10 (+QTL), and IR77298-14-1-2-B-13 (−QTL) that showed contrasting performance for grain yield under drought. We presumed that the drought tolerant high yielding BILs contained QTL and drought susceptible low yielding BILs did not, and designated them as +QTL and –QTL lines respectively. Subsequently, +QTL lines (IR77298-5-6-B-18, IR77298-14-1-2-B-10) as well as IR77298-14-1-2, a drought tolerant line from the same population, were crossed with –QTL lines and with IR64 in different combinations to develop six mapping populations (Table S1) to identify, introgress, and pyramid QTL for grain yield under drought in the IR64 background.
Phenotypic Characterization
All field-based phenotyping experiments were conducted at the International Rice Research Institute in lowland transplanted conditions (IRRI, Los Baños, Philippines, 14°30′N, 121°15′E). Drought stress (stress) and well-watered control (non-stress) experiments were managed as described by Venuprasad et al. [17].
To characterize the physiological response of +QTL and −QTL BILs, three lowland field experiments were conducted (2009 wet season (2009WS), 2010 dry season (2010DS), and 2010DS in a rainout shelter (2010DS-ROS). The soil was classified as an Isohyperthermic Typic Hapludalf in 2009WS and 2010DS, and Aquandic Epiaquall at the site used in 2010DS-ROS, ranging in bulk density at 30 cm from 0.86 to 1.3 g cm−3. Four replicates per genotype were planted in randomized compete block designs in plots of 3-m length with 4 rows per plot. Soil water potential was monitored in all seasons with three tensiometers installed in each experimental field (Soilmoisture Equipment Corp, CA, USA) and it reached −70 kPa at a depth of 20 cm from 60–75 days after sowing in 2010DS, indicating that the most severe drought treatment was applied in that season (Fig. S2). Canopy temperature was measured at mid-day from a 3-m ladder by infrared thermography (NEC Avio Infrared Technologies Co. Ltd., Tokyo, Japan). Leaf stomatal conductance was measured at mid-day with a porometer (AP4, Delta-T Devices, Cambridge, UK). Shoot growth dynamics were monitored by destructive harvests in 2009WS and 2010DS, and with NDVI measured around mid-day in 2010DS-ROS (Greenseeker Hand-held Sensor, NTech Industries, CA, USA). Root samples were taken at 94 DAS in 2009WS, at 73 DAS in 2010DS, and at 94 DAS in 2010DS-ROS to a depth of 60 cm using a 4-cm-diameter core sampler according to Henry et al. [38].
Six BC4F3 mapping populations were screened for grain yield under stress and non-stress conditions during the dry seasons (Table S1). An alpha lattice design with two replications and two 4-m rows per plot was used in all experiments. Based on cumulative rainfall, rainfall distribution, and yield reductions compared to the non-stress treatment, a range of drought severities was achieved. Severe droughts occurred in 2008DS, 2010DS, 2011DS, and 2012DS. In all experiments, data on DTF, PH, and GY per plot (normalized to 12% moisture content) were recorded as described in the standard evaluation system of IRRI [39].
Defining QTL Regions using Affymetrix Rice Genome Array DNA Analysis
Based on screening with 600 SSR markers, low polymorphism was detected between the +QTL and –QTL BILs due to their similar genetic background. We used genomic DNA hybridization on an Affymetrix GeneChip® Rice Genome Array to determine regions with significant differences between BILs of one pair (IR77298-5-6-B-18 (+QTL) and IR77298-5-6-B-11 (−QTL)). Probeset intensity and differentially hybridized probesets (at p<0.05) between the +QTL and –QTL BILs were determined using the R/Bioconductor packages Affy and R/MAANOVA. To identify genome blocks that were polymorphic between the lines of interest, aggregation of differentially hybridized probesets along the genome was determined at 1000 kb, with sliding windows of 500 kb as described by Bruce et al. [23]. Out of 17 candidate regions, 4 regions that showed significant differences between the +QTL and –QTL BILs in the GeneChip results were selected for further QTL analysis. A polymorphism survey between parents of the six mapping populations was carried out on the four regions. A total of 15, 29, 24, and 33 markers in the four differentially expressed regions on chromosomes 2, 4, 9, and 10, respectively, were surveyed for polymorphism and we detected 3, 6, 4 and 5 polymorphic markers in those regions. Six populations, P1, P2, P3, P4, P5 and P6, with 487, 478, 455, 286, 485 and 488 individuals, respectively, were genotyped using the polymorphic markers and phenotyped for yield and yield-related traits under drought stress and non-stress conditions.
Statistical Analysis
Statistical analysis was carried out using SAS v.9.1.3 [40] for the phenotypic data of mapping populations, and means of trait values were estimated using the REML algorithm of PROC MIXED. Statistical analyses for the physiology experiments were performed in R v. 2.8.0 [41] using ANOVA, and Tukey’s HSD was used for mean comparison.
QTL Identification, and Differential Expressed Gene Analyses
QTL analysis was carried out following Interval mapping and Composite Interval Mapping (CIM) using Q gene software [42] and results were validated by QTL Cartographer [43]. A significance level of 0.01 was used to detect putative QTL. The minimal LOD value required to declare a QTL was obtained from 1000 permutation tests. QTL analyses were confirmed by performing CIM and iQTLM analysis in connected populations using MC-QTL software [44].
To identify differentially expressed genes within the four QTL regions, a reanalysis of the Agilent 44k oligoarray data was conducted using R/MAANOVA (significance p<0.001) to compare +QTL and –QTL BILs rather than drought treatments (FTSW 0.5 and 0.2; moderate and severe drought stress) as reported by Moumeni et al. [21]. The DEGs determined from the reanalysis were tested for enrichment in the four QTL regions (gene models from the Michigan State University Rice Genome Annotation Project, release 6.1: http://rice.plantbiology.msu.edu/) using a Fisher exact test (significance p<0.05, method from McNally et al. [45]. The same enrichment test was then made for all DEGs against the GO-SLIM, Mapman (http://mapman.gabipd.org/), and Gramene QTL (http://gramene.org/qtl/) annotation of the MSU 6.1 genome. All microarray data and information on the transcriptomes of the three tissues are available at NCBI GEO accessions GSE30463 (root expression data), GSE30449 (shoot expression data) and GSE30462 (panicle expression data) (http://www.ncbi.nlm.nih.gov/geo/).
Identification of IR64 Introgression Lines
From among 2806 lines in six populations, a set of 84 lines with high yield under both stress and non-stress conditions, and with phenotypic similarity to IR64, was genotyped with foreground markers for QTL qDTY2.2, qDTY4.1, qDTY9.1, and qDTY10.1 to identify lines with combinations of one, two, three, and four QTL. Background genotyping of the 84 lines was carried out with 580 SSR markers randomly distributed throughout the genome to identify introgression lines with high genetic similarity to IR64. In 2011DS, 84 lines were evaluated for yield under stress and non-stress conditions. These 84 lines were analyzed for amylose content (AC), gelatinization temperature (GT), milling potential (MP), and chalkiness (CS) following the protocol as described in the standard evaluation system of IRRI (39) with the aim to identify lines with quality traits similar to those of IR64. The three best identified IR64 QTL lines (Table 6) were evaluated under non-stress and drought stress conditions at one site each in Bangladesh (Rajshahi) and Nepal (Nepalganj) and at four sites in India (Hyderabad, Rewa, Raipur, and Hazaribagh).
Supporting Information
Figure S1.
QTLs identified for grain yield under drought stress.
https://doi.org/10.1371/journal.pone.0062795.s001
(DOCX)
Figure S2.
Major effect QTL DTY2.2 and DTY9.1 identified in multiple populations.
https://doi.org/10.1371/journal.pone.0062795.s002
(DOCX)
Figure S3.
Soil water potential measured by tensiometers in each field study season.
https://doi.org/10.1371/journal.pone.0062795.s003
(DOCX)
Table S1.
Mean for yield and related traits under drought stress and non-stress conditions during 2007DS to 2010DS.
https://doi.org/10.1371/journal.pone.0062795.s004
(DOCX)
Table S2.
QTLs for yield related traits under drought stress and non-stress in IR64×Aday Sel derived populations.
https://doi.org/10.1371/journal.pone.0062795.s005
(DOCX)
Table S3.
QTLs for grain yield under drought identified by multi-population joint analyses.
https://doi.org/10.1371/journal.pone.0062795.s006
(DOCX)
Table S4.
Regions of aggregation of differentially expressed genes from the re-analyzed transcriptome dataset.
https://doi.org/10.1371/journal.pone.0062795.s007
(DOCX)
Table S5.
Enriched Gramene QTL categories of the subset candidate genes in IR77298-14-1-2-B-10 (+QTL) and IR77298-14-1-2-B-13 (−QTL).
https://doi.org/10.1371/journal.pone.0062795.s008
(DOCX)
Table S6.
Performance of 84 IR64-NILs under non –stress (NS) and drought stress (S) conditions.
https://doi.org/10.1371/journal.pone.0062795.s009
(DOCX)
Table S7.
List of candidate genes associated with yield-under-drought QTLs.
https://doi.org/10.1371/journal.pone.0062795.s010
(DOCX)
Table S8.
Candidate genes/gene families and their biological functions under stress conditions.
https://doi.org/10.1371/journal.pone.0062795.s011
(DOCX)
Acknowledgments
The authors acknowledge the development of initial BIL pairs by R. Venuprasad and G. Atlin, and technical support from M. Del Valle, L. Quiatchon, J. Guevarra, M.T. Sta Cruz, and P. Maturan for the QTL study and R. Torres, L. Holongbayan, N. Turingan, L. Satioquia, for the physiological study.
Author Contributions
Conceived and designed the experiments: AK MSBP AH RM. Performed the experiments: MSBP AH HUA RM RT SBV PP NPM MV RS RC ONS JLD SPD KKM RBY TLA BK KS. Analyzed the data: BPMS HUA AH RM SD PV. Contributed reagents/materials/analysis tools: AK KS AM SK HL. Wrote the paper: MSBP AH HL RM AK.
References
- 1.
O’ Toole JC (1982) Adaptation of rice to drought prone environments. In: Drought resistance in crops with emphasis on rice. Los Baños (Philippines): International Rice Research Institute. 195–213.
- 2. Verulkar SB, Mandal NP, Dwivedi JL, Singh BN, Sinha PK, et al. (2010) Breeding resilient and productive genotypes adapted to drought prone rainfed ecosystems of India. Field Crop Res 117: 197–208.
- 3. Venuprasad R, Bool ME, Dalid CO, Bernier J, Kumar A, et al. (2009a) Genetic loci responding to two cycles of divergent selection for grain yield under drought stress in a rice breeding population. Euphytica 167: 261–269.
- 4. Lafitte HR, Yongsheng G, Yan S, Li ZK (2007) Whole-plant responses, key processes, and adaptation to drought stress: the case of rice. J Exp Bot 58: 169–175.
- 5.
Mackill D, Coffman W, Garrity D (1996) Rainfed Lowland Rice Improvement. International Rice Research Institute. Manila, Philippines.
- 6. Courtois B, McLaren G, Sinha PK, Prasad K, Yadava R, et al. (2000) Mapping QTLs associated with drought avoidance in upland rice. Mol Breed6: 55–66.
- 7. Mandal NP, Sinha PK, Variar M, Shukla VD, Perraju P, et al. (2010) Implications of genotype×input interactions in breeding superior genotypes for favorable and unfavorable rainfed upland environments. Field Crop Res 118: 135–144.
- 8.
Blum A (1988) Plant breeding for stress environments. Boca Raton, Florida (USA): CRC Press. 222 pp.
- 9. Sinclair TR (2011) Challenges in breeding for yield increase for drought. Trends Plant Sci 16: 289–293.
- 10. Kumar A, Bernier J, Verulkar S, Lafitte HR, Atlin GN (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 Crop Res 107: 221–231.
- 11. Shen L, Courtois B, McNally KI, Robin S, Li Z (2001) Evaluation of near-isogenic lines of rice introgressed with QTLs for root depth through marker-aided selection. Theor Appl Genet 103: 75–83.
- 12. Price AH, Steele KA, Moore BJ, Barraclough PB, Clark LJ (2000) A combined RFLP and AFLP linkage map of upland rice (Oryza sativa L.) used to identify QTLs for root-penetration ability. Theor Appl Genet 100: 49–56.
- 13. Steele KA, Price AH, Shashidhar HE, Witcombe JR (2006) Marker assisted selection to introgress rice QTLs controlling root traits into an Indian upland rice variety. TheorAppl Genet 112: 208–221.
- 14. Steele KA, Virk DS, Kumar R, Prasad SC, Witcombe JR (2007) Field evaluation of upland rice lines selected for QTLs controlling root traits. Field Crop Res 101: 180–186.
- 15. Bernier J, Kumar A, Venuprasad R, Spaner D, Atlin G (2007) A large-effect QTL for grain yield under reproductive-stage drought stress in upland rice. Crop Sci 47: 505–516.
- 16. Kumar R, Venuprasad R, Atlin GN (2007) Genetic analysis of rainfed lowland rice drought tolerance under naturally-occurring stress in eastern India: heritability and QTL effects. Field Crop Res 103: 42–52.
- 17. Venuprasad R, Dalid CO, Delvalle M, Zhao D, Espiritu M, et al. (2009b) Identification and characterization of large-effect quantitative trait loci (QTL) for grain yield under lowland drought stress and aerobic conditions in rice using bulk-segregant analysis (BSA). Theor Appl Genet120: 177–190.
- 18. Vikram P, Swamy BPM, Dixit S, Ahmed HU, Sta Cruz MT, et al. (2011) qDTY1.1, a major QTL for rice grain yield under drought with consistent effect in multiple elite genetic backgrounds. BMC Genet 12: 89.
- 19. Xiao J, Li J, Yuan L, Tanksley SD (1996) Identification of QTLs affecting traits of agronomic importance in a recombinant inbred population derived from a sub specific cross. TheorApplGenet 92: 230–244.
- 20. Tanksley SD, Nelson JC (1996) Advanced backcross QTL analysis: a method for the simultaneous discovery and transfer of valuable QTLs from unadapted germplasm into elite breeding lines. Theor Appl Genet 92: 191–203.
- 21. Moumeni A, Satoh KK, Kondoh HH, Asano TT, Hosaka AA, et al. (2011) Comparative analysis of root transcriptome profiles of two pairs of drought-tolerant and susceptible rice near-isogenic lines under different drought stress. BMC Plant Biol 11: 174.
- 22. Venuprasad R, Impa SM, Veeresh Gowda RP, Atlin GN, Serraj R (2011) Rice near- isogenic lines (NILs) contrasting for grain yield under lowland drought stress. Field Crop Res 123: 38–46.
- 23. Bruce MA, Hess A, Bai J, Mauleon R, Diaz MG, et al. (2009) Detection of genomic deletions in rice using oligonucleotide microarrays. BMC Genomics 10: 129.
- 24. Nguyen HT, Babu RC, Blum A (1997) Breeding for drought resistance in rice: physiology and molecular genetics considerations. Crop Sci 37: 1426–1434.
- 25. Zhang J, Zheng HG, Aarti AA, Pantuwan G, Nguyen TT, et al. (2001) Locating genomic regions associated with components of drought resistance in rice: comparative mapping within and across species. TheorAppl Genet 103: 19–29.
- 26. Robin S, Pathan MS, Courtois B, Lafitte R, Carandang S, et al. (2003) Mapping osmotic adjustment in an advanced back-cross inbred population of rice. Theor Appl Genet 107: 1288–1296.
- 27. Yue B, Xiong L, Xue W, Xing Y, Luo L, et al. (2005) Genetic analysis for drought resistance of rice at reproductive stage in field with different types of soil. Theor ApplGenet 111: 1127–1136.
- 28. Lanceras JC, Pantuwan G, Jongdee B, Toojinda T (2004) Quantitative trait loci associated with drought tolerance at reproductive stage in rice. Plant Physiol 135: 384–399.
- 29. Babu RC, Nguyen BD, Chamarerk V, Shanmugasundaram P, Chezhian P, et al. (2003) Genetic analysis of drought resistance in rice by molecular markers: Association between secondary traits and field performance. Crop Sci 43: 1457–1469.
- 30. Lafitte HR, Price AH, Courtois B (2004) Yield response to water deficit in an upland rice mapping population: association among traits and genetic markers. Theor ApplGenet 109: 1237–1246.
- 31. Srividya A, Reddy VLN, Hariprasad AS, Jayaprada M, Sridhar S, et al. (2010) Identification and mapping of landrace derived QTL associated with yield and its components in rice under different nitrogen levels and environments. Int J Plant Breeding Genet 4: 210–227.
- 32. Brondani C, Rangel PHN, Brondani RPV, Ferreira ME (2002) QTL mapping and introgression of yield-related traits from Oryza glumaepatula to cultivated rice (Oryza sativa) using microsatellite markers. Theor Appl Genet 104: 1192–1203.
- 33. Hittalmani S, Huang N, Courtois B, Venuprasad R, Shashidhar HE, et al. (2003) Identification of QTL for growth and grain yield-related traits in rice across nine locations of Asia. TheorAppl Genet 107: 679–690.
- 34. Thomson MJ, Tai TH, McClung AM, Lai XH, Hinga ME, et al. (2003) Mapping quantitative trait loci for yield, yield components and morphological traits in an advanced backcrossing population between Oryza rufipogon and the Oryza sativa cultivar Jefferson. Theor Appl Genet 107: 479–493.
- 35. Zou GH, Mei HW, Liu HY, Liu GL, Hu SP, et al. (2005) Grain yield response to moisture regimes in a rice population: association among traits and genetic markers. Theor Appl Genet 112: 106–113.
- 36. Lin MH, Lin CW, Chen JC, Lin YC, Cheng SY, et al. (2007) Tagging rice drought-related QTL with SSR DNA markers. Crop, Environ Bioinformatics 4: 65–76.
- 37. Dixit S, Swamy BPM, Vikram P, Ahmed HU, StaCruz MT, et al. (2011) Fine mapping four large effect QTLs for rice grain yield under drought. Theor Appl Genet 125: 155–169.
- 38. Henry A, Veeresh Gowda RP, Torres RO, McNally KL, Serraj R (2011) Genetic variation in root architecture and drought response in Oryza sativa: rainfed lowland field studies of the Oryza SNP panel. Field Crop Res 120: 205–214.
- 39.
IRRI (International Rice Research Institute) (1996) SES (Standard Evaluation System for Rice). International Network for Genetic Evaluation of Rice. IRRI, Los Baños, Philippines.
- 40.
SAS Institute, Inc. (2004) SAS online doc 9.1.3. SAS Institute Inc., Cary, NC.
- 41.
R Development Core Team (2008) A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria.
- 42. Nelson JC (1997) Software for marker-based genomic analysis and breeding. Mol Breeding 3: 239–245.
- 43.
Wang S, Basten CJ, Zeng ZB (2010) Windows QTL Cartographer 2.5. Department of Statistics. Available: http://statgen.ncsu.edu/qtlcart/winqtl cart.htm.
- 44. Jourjon MF, Jasson S, Marcel J, Ngom B, Mangin B (2005) MCQTL: multi-allelic QTL mapping in multi-cross design. Bioinformatics 21: 128–130.
- 45. McNally KL, Childs KL, Bohnert R, Davidson RM, Zhao K, et al. (2009) Genome wide SNP variation reveals relationships among landraces and modern varieties of rice. Proc Natl Aca Sci USA 106: 12273–12278.