Water scarcity and drought have seriously threatened traditional rice cultivation practices in several parts of the world, including India. Aerobic rice that uses significantly less water than traditional flooded systems has emerged as a promising water-saving technology. The identification of QTL conferring improved aerobic adaptation may facilitate the development of high-yielding aerobic rice varieties. In this study, experiments were conducted for mapping QTL for yield, root-related traits, and agronomic traits under aerobic conditions using HKR47 × MAS26 and MASARB25 × Pusa Basmati 1460 F2:3 mapping populations.
A total of 35 QTL associated with 14 traits were mapped on chromosomes 1, 2, 5, 6, 8, 9, and 11 in MASARB25 x Pusa Basmati 1460 and 14 QTL associated with 9 traits were mapped on chromosomes 1, 2, 8, 9, 10, 11, and 12 in HKR47 × MAS26. Two QTL (qGY8.1 with an R2 value of 34.0% and qGY2.1 with an R2 value of 22.8%) and one QTL (qGY2.2 with an R2 value of 43.2%) were identified for grain yield under aerobic conditions in the mapping populations MASARB25 × Pusa Basmati 1460 and HKR47 × MAS26, respectively.
A number of breeding lines with higher yield per plant, root length, dry biomass, length-breadth ratio, and with Pusa Basmati 1460-specific alleles in a homozygous or heterozygous condition at the BAD2 locus were identified that will serve as novel material for the selection of stable aerobic Basmati rice breeding lines.
Our results identified positive correlation between some of the root traits and yield under aerobic conditions, indicating the role of root traits for improving yield under aerobic situations possibly through improved water and nutrient uptake. Co-localization of QTL for yield, root traits, and yield-related agronomic traits indicates that the identified QTL may be immediately exploited in marker-assisted-breeding to develop novel high-yielding aerobic rice varieties.
Keywords:Oryza sativa; QTL; Rice; Aerobic; Dry direct seeded
Water scarcity is one of the most pressing issues facing agriculture today. In many countries, water for agriculture consumes about 70% of the total freshwater use. To meet the needs of a growing population, more food needs to be produced with less water . Rice (Oryza sativa L.) is the primary source of food for more than half of the world’s population. Rice is cultivated in highly diverse situations that range from flooded wetland to rainfed dryland . Irrigated rice, which accounts for 55% of the world rice area, provides 75% of global rice production and consumes about 90% of the freshwater resources used for agriculture in Asia . Water deficit is a key constraint that affects rice production in different countries. The affinity of the rice crop with water is universally known. Conventional rice production ecosystems (puddled transplanted) require an average of 2,500 liters of water to produce 1 kg of rough rice , which is 2-3 times more than what other cereals require [5,6]. It is also reported that about 50% of the diverted fresh water in Asia is used to irrigate rice fields . Seasonal water inputs for puddled transplanted rice vary from 660 to 5,280 mm depending on growing season, climatic conditions, soil type, and hydrological conditions . Although evapo-transpirational losses of water in rice are similar to those in wheat , the higher water requirements for rice cultivation are mainly due to the water required for puddling, and high seepage and percolation losses associated with continuous flooding .
Rice can be established by transplanting seedlings in puddled fields or by direct seeding in dry or puddled fields . Although the intensive water and labor requirements in transplanting of rice in puddled fields are well known, technologies such as dry and wet direct seeding and alternate wetting and drying (AWD) could be an option to produce rice in both irrigated and rainfed rice ecosystems. Aerobic rice is one such extensive water-saving technology for rice. Aerobic rice refers to a cultivation system in which rice is dry direct seeded in well-tilled leveled fields with uniform slope under unpuddled conditions. The crop is cultivated under aerobic conditions with no standing water throughout the season. In this system, rice can be established using different systems such as broadcasting, drilling, or dibbling in well-prepared fields and direct seeding with zero tillage using a mechanical seed drill or raised beds . Traditionally, this method has been practiced in rainfed upland and rainfed shallow lowland areas of Asia . However, the high proportion of water savings associated with this method compared with conventional rice-growing practices has made this method increasingly popular in irrigated areas where the problem of water shortage is also expanding . Under favorable irrigated conditions, rice is drilled in well-prepared fields, and cultivated with efficient weed control and uniform irrigation throughout the crop season. Aerobic rice systems use less water than conventional flooded rice  through the use of rice varieties capable of responding well to reduced water inputs in non-puddled and non-saturated soils [11,12]. In the context of current and predicted water scarcity, increasing irrigation is generally not a viable option for alleviating drought problems in rainfed rice-growing systems . It is therefore critical that genetic management strategies be undertaken for cultivating rice with less water and focusing on maximum extraction of available soil moisture and its efficient use in crop establishment and growth to maximize biomass and yield. With rice being grown extensively in the upland ecosystem, extensive genetic variation for aerobic adaptation exists in rice germplasm. However, the current challenge is to decipher the complexities of aerobic adaptation in rice and exploit all available genetic resources to produce rice varieties combining aerobic adaptation with high yield potential, good grain quality, and resistance to biotic stresses. This involves the development of high-throughput, high-precision phenotyping systems to allow genes for yield components under aerobic conditions to be efficiently mapped and their effects assessed for a range of traits, and then move the most promising genes into widely grown rice mega-varieties, while scaling up gene detection and delivery for use in marker-aided breeding.
One important aspect for understanding aerobic adaptation is the response of root growth and root development under aerobic conditions, including efficient uptake of water and nutrients. Roots are important for maintaining crop yields, which is vital when plants are grown in soils containing insufficient supplies of water or nutrients , and roots are one of the primary sites for stress signal perception that initiates a cascade of gene expression responses to water deficit . Previous studies showed that plant growth largely depends on the severity of the stress; mild water deficit leads to growth inhibition of leaves and stems, whereas roots may continue to elongate . Furthermore, root architecture is a key trait for dissecting the genotypic differences in rice responses to water deficit .
An understanding of the underlying physiological and molecular mechanisms is necessary to improve the adaptation of rice varieties to aerobic conditions. For water-deficit conditions, Price et al.  reported that a long and thick root system, the ratio of root weight to shoot weight, and root penetration ability of upland rice contribute greatly to drought tolerance. Several of the QTL identified for root length are consistent across mapping populations  and common genomic regions across populations and even across species have been identified for root thickness, root penetration, and stomatal behavior . Progress has been made in detecting large-effect QTL conferring drought tolerance in lowland and irrigated rice [17,20]. Several QTL for grain yield under drought stress have been reported for both upland and lowland rice [21,22]. However, very few reports identify the genomic regions responsible for increased aerobic adaptation of rice. In our study, two aerobic × non-aerobic populations were studied for root traits, grain yield, and yield-related agronomic traits to identify QTL conferring a potential yield advantage under aerobic conditions; identify QTL for yield-related agronomic traits under aerobic conditions; identify root traits having a positive correlation with grain yield under aerobic conditions; identify QTL for root-related traits with a positive correlation with grain yield under aerobic conditions; and identify QTL for yield, yield-related agronomic traits, and root traits that co-localize with the aim to introgress identified QTL in HKR47 and Pusa Basmati 1460 rice varieties using marker-assisted breeding to develop rice varieties better adapted to aerobic conditions.
Two mapping populations (HKR47 × MAS26 and MASARB25 × Pusa Basmati 1460) derived from crosses of HKR47 and Pusa Basmati 1460, lowland high-yielding rice cultivars unadapted to cultivation in aerobic conditions, with upland and aerobic adapted genotypes MAS26 (developed at the University of Agricultural Sciences, Bangalore) and MASARB25 (developed at IRRI, Manila, Philippines), respectively, were used in our study.
Development and management of mapping populations
The F2:3 populations were developed from crosses involving MAS26 and MASARB25 as male parents and HKR47 and Pusa Basmati 1460 as female parents. True F1 seeds from the cross were grown to get F2 seeds. Seeds harvested from HKR47 × MAS26 and MASARB25 × Pusa Basmati 1460 F1 plants were germinated in petri plates and then a total of 235 and 250 germinated seeds from the two populations, respectively, were raised to maturity in pots (25 cm diameter × 25 cm height) with one plant/pot in a net house in 2010. The pots were irrigated with 1 liter of water for the first 15 days, and then with 1 liter after every third day up to panicle emergence. After every fifth day, the pots were irrigated with full-strength nutrient solution for the first 21 days, and with half-strength nutrient solution thereafter.
Seeds obtained from each F2 plant were evaluated in fields at the Rice Research Station, Kaul, Karnal, India, during 2011 and 2012 under dry direct-seeded aerobic cultivation practices that involved dry seeding at approximately 2-cm depth in dry-ploughed and harrowed aerobic plots with row spacing of 30 cm (2011) and 25 cm (2012), resulting in a seed rate of approximately 305 seeds m-2 in 2011 and 365 seeds m-2 in 2012. The seeds obtained from each selected F2 plant (20% best and 20% worst, 94 plants in the case of HKR47 × MAS26 and 100 plants in the case of MASARB25 × Pusa Basmati 1460) were grown in 2-m single-row plots with two replications using a seeding density of 2 g per linear meter of row to record observations for the root study, grain yield, agronomic traits, and genotyping. For the yield trial, the plots were randomized by using Crop Stat version 7.2. Aerobic fields were irrigated for about 1 week with a 2-3-cm water layer to facilitate crop establishment; thereafter, the fields were re-irrigated once at a 10-day interval. Nitrogen was applied at 60 kg ha-1 basal after sowing, 60 kg ha-1 25 days after sowing, and 60 kg ha-1 55 days after sowing. In addition, 30 kg P ha-1, 40 kg K ha-1, and 5 kg Zn ha-1 were applied as basal and 30 kg P ha-1 at 25 days after sowing. The plots were kept weed-free by manual weeding. At physiological maturity or 80-85% of the panicles turned into golden yellow and the panicles at the base were already at hard dough stage, data were recorded on agronomic traits, plant height in cm (PH), effective number of tillers per plant (TN), panicle length in cm (PL), number of panicles per plant (P/P), number of grains per panicle (S/P), 1,000-grain weight (1,000 GW), length/breadth ratio (L/B), and grain yield in Kg ha-1 (GY). Grain yield per plot was recorded after harvesting, threshing, and drying to moisture content adjusted to 14% and converted to kg ha-1. The data on root morphological traits—root length (RL, cm), root thickness (RT), root number (RN), root volume (RV), fresh and dry root weight (FRW and DRW in g), and fresh and dry shoot weight (FSW and DSW in g) from six plants from each line at maturity were recorded and analyzed. For the measurement of root traits, plants were removed from the soil. Soil sections that contained roots were removed by digging a hole and put into a sieve and gently washed with a hose until only the roots remained. The roots and shoots were then separated by cutting from the top soil line. For fresh root weight, the roots were blotted gently with a soft paper towel to remove any free surface moisture. Then, the roots and shoots were weighed immediately (plants have a high composition of water, so waiting to weigh them may lead to some drying and therefore produce inaccurate data). Root volume was measured by the actual volume displacement analysis, which measures the volume of water displaced when roots are submerged in a vessel of water. Volume displacement analysis has the advantage of providing a fast measure of new root production. In addition, its non-destructive nature permits repeated measures over time. The thickness of the root crown was measured using a vernier caliper. Root numbers were counted manually. Root length was measured using a centimeter scale. For dry root weight, the roots and shoots were dried in an oven set to low heat (100°F) for 3 days, and then cooled in a dry environment. Once cooled, they were weighed on a scale.
Statistical analysis of the data on individual characters was performed using SAS v9.1 (SAS Institute, Inc., 2002-2003). Line means were estimated using the REML option of the SAS MIXED procedure taking lines as fixed and replicates and blocks within replicates as random. The phenotypic r(P) correlation coefficients for all possible pairs of characters were calculated from the variance and covariance was estimated:
σ x y (P) = Phenotypic covariance between characters x and y
σ2 x (P) = Phenotypic variance of character x
σ2 y (P) = Phenotypic variance of character y
The phenotypic correlation coefficients were tested against a standardized tabulated significant value of r with (n-2) degree of freedom as per the procedure of Fisher and Yates .
Genotyping of mapping populations
All DNA marker work was conducted in the Department of Molecular Biology and Biotechnology, CCSHAU, Hisar, India. Six fresh leaves from each line were collected in bulk and genomic DNA was isolated from 1-month-old plants using the CTAB method . DNA quantity was estimated by ethidium bromide staining on 1% agarose gels using a standard containing 100 ng/μl genomic DNA. PCR amplification, denaturing polyacrylamide gel electrophoresis, and silver staining were essentially carried out as described earlier by Jain et al. .
A total of 300 mapped SSR and aroma gene-specific BAD2 markers (for which basmati serves as a parent) were screened for polymorphism between the parents. The markers were obtained based on published rice genome maps (IRGSP 2005) and their physical position (Mb) on the indica genome (http://www.gramene.org webcite) was used as a reference. A total of 91 and 112 markers showed polymorphism and were run on 94 and 100 plants for HKR47 x MAS26 and MASARB25 × Pusa Basmati 1460 populations, respectively. Genetic similarities between the cultivars were measured by the similarity coefficient based on the proportion of shared electromorphs using the 'Simqual’ subprogram of the NTSYS-PC (Version 2.02 Exeter Software, Setauket, NY, USA) software package . The resultant distance matrix data were used for two-dimensional scaling of rice genotypes by principal component analysis (PCA). Linkage analysis was performed using the Map Manager/QTX computer program  using the Kosambi function and linkage evaluation of P = 0.001. The ripple command was used to verify the marker order. QTL analysis was performed using Windows QTL Cartographer version 2.5 . For interval mapping (IM) analysis, an LOD threshold score of 2.5 was selected. The proportion of the total phenotypic variation explained by each QTL was calculated as R2 value (R2 = ratio of the sum of squares explained by the QTL to the total sum of squares). To more accurately determine QTL positions, composite interval mapping (CIM) was performed with default parameters (permutation time 300, significance level 0.05, model 6; standard model, method 3; forward and backward method, walk speed 2 cM, etc.).
Phenotypic analysis of parents and populations
MAS26 and MASARB25, the aerobic genotypes, outyielded the lowland parents, HKR47 and Pusa Basmati 1460, in both years under dry direct-seeded aerobic conditions. Root length, root volume, root thickness, root number, and root biomass of the tolerant parent were found to be higher than those of the susceptible parent. Mean values of RL, RV, RT, RN, FRW, DRW, FSW, and DSW for the two populations are presented in Table 1. Mean and range values of PH, TN, PL, P/P, S/P, 1,000 GW, and GY on two populations under dry direct-seeded aerobic conditions in 2011 and 2012 are presented in Table 2.
Table 1. Mean and range for various root traits in two populations under pot-house conditions
Table 2. Mean and range for various agronomic traits in two populations under direct-seeded aerobic field conditions
In both mapping populations, grain yield was found to be positively correlated with root length, root number, root volume, and root biomass, suggesting that these traits might be the key factor for improving grain yield under aerobic conditions (Table 3).
Table 3. Phenotypic correlation coefficients between agronomic and root traits in MASARB25 × Pusa Basmati and HKR47 × MAS26-derived populations
Phenotypic correlation coefficient analysis of MASARB25 × Pusa Basmati 1460 showed that GY is positively and significantly correlated with the effective number of tillers (0.442, p = 0.01), panicle length (0.411, p = 0.01), panicles plant-1 (0.427, p = 0.01), seeds panicle-1 (0.377, p = 0.01), and length/breadth ratio (0.689, p = 0.01). In this population, root number is positively and significantly correlated with panicles plant-1 (0.365, p = 0.01), panicle length (0.324, p = 0.01), seeds panicle-1 (0.367, p = 0.01), GY (0.285, p = 0.05), dry root weight (0.575, p = 0.01), fresh root weight (0.630, p = 0.01), and root length (0.602, p = 0.01). Root thickness has been found to be positively and significantly correlated with fresh root weight (0.544, p = 0.01), root length (0. 0.352, p = 0.01), root number (0.770, p = 0.01), dry root weight (575, p = 0.01), panicles plant-1 (0.335, p = 0.01), and GY (0.278, p = 0.05). Root volume is positively and significantly correlated with panicle length (0.358, p = 0.01), dry root weight (0.575, p = 0.01), fresh root weight (0.621, p = 0.01), root thickness (0.587, p = 0.01), root number (0.740, p = 0.01), and root length (0.700, p = 0.01). Dry root weight is positively and significantly correlated with GY (0.330, p = 0.01) and seeds panicle-1 with panicle length (0.604, p = 0.01). Root length was found to be negatively correlated with plant height (Table 3).
Phenotypic correlation coefficient analysis of HKR47 × MAS26 showed that GY is positively and significantly correlated with panicle length (0.512, p = 0.01), panicles plant-1 (0.294, p = 0.05), seeds panicle-1 (0.476, p = 0.01), and length/breadth ratio (0.725, p = 0.01). In this population, seeds panicle-1 is positively and significantly correlated with panicle length (0.694, p = 0.01) and panicles plant-1 (0.529, p = 0.01). Root number is positively and significantly correlated with GY (0.386, p = 0.05) and fresh root weight (0.482, p = 0.01). Root thickness has been found to be positively and significantly correlated with fresh root weight (0.626, p = 0.01), root number (0.746, p = 0.01), and GY (0.408, p = 0.01). Root volume is positively and significantly correlated with fresh root weight (0.412, p = 0.01), root thickness (0.728, p = 0.01), and GY (0.455, p = 0.01). Root length is positively and significantly correlated with GY (0.317, p = 0.05) and seeds panicle-1 (0.355, p = 0.01). Fresh root weight is positively and significantly correlated with panicles plant-1 (0.412, p = 0.01), GY (0.346, p = 0.01), and effective number of tillers (0.539, p = 0.01) (Table 3).
Genotypic analysis of mapping populations
Polymorphism of HKR47 with MAS26 and of MASARB25 with Pusa Basmati 1460 was 31.7% and 37.3%, respectively. A total of 226 and 184 alleles were identified in the HKR47 × MAS26 and MASARB25 Pusa Basmati 1460 populations by using 300 SSR markers, respectively. MASARB25 × Pusa Basmati 1460-derived mapping populations were also evaluated for aroma using aroma gene-specific BAD2 markers. A total of 13 lines were found to be homozygous and 24 lines were found to be heterozygous for aroma.
Similarity coefficient data based on the proportion of shared alleles using 91 and 112 SSR markers were used to calculate the coefficient values among the selected HKR47 × MAS26 and MASARB25 × Pusa Basmati 1460 lines, respectively, through UPGMA tree cluster analysis. Allelic diversity was used to produce a dendrogram (cluster tree analysis, NTSYS-PC) to demonstrate the genetic relationship among selected plants and the parental rice varieties. All 94 and 100 plants from the two populations, HKR47 × MAS26 and MASARB25 × Pusa Basmati 1460, respectively, clustered in two major groups at the similarity coefficients of 0.520 and 0.540, respectively. Genetic relationship was also assessed by PCA analysis (NTSYS-PC). Two-dimensional PCA scaling showed that the parental genotypes were quite distinct, whereas the mapping population genotypes were interspersed between the two parental lines (Figures 1 and 2).
Figure 1. Two-dimensional PCA scaling of MASARB25 × Pusa Basmati 1460 based on SSR diversity data.
Figure 2. Two-dimensional PCA scaling of HKR47 × MAS26 based on SSR diversity data.
Several QTL associated with root traits and agronomically important traits in the two mapping populations were detected in the MASARB25 x Pusa Basmati 1460 (Table 4) and HKR47 × MAS26 (Table 5) populations. A total of 35 QTL associated with 14 traits were mapped on chromosomes 1, 2, 5, 6, 8, 9, and 11 in MASARB25 × Pusa Basmati 1460-derived population (Table 4, Additional file 1: Figure S1). A total of 14 QTL associated with nine traits were mapped on chromosomes 1, 2, 8, 9, 10, 11, and 12 in HKR47 × MAS26-derived population (Table 5, Additional file 1: Figure S2).
Table 4. QTL associated with agronomic and aerobic root traits of MASARB25 × Pusa Basmati 1460 population
Table 5. QTL associated with agronomic and aerobic root traits of HKR47 × MAS26 population
Additional file 1: Figure S1. QTL associated with agronomic and aerobic root traits of MASARB25 × Pusa Basmati 1460 population. Figure S2. QTL associated with agronomic and aerobic root traits of HKR47 × MAS26 population.
Format: DOC Size: 204KB Download file
This file can be viewed with: Microsoft Word Viewer
QTL qRL8.1 with a peak at RM544 on chromosome 8 with an R2 value of 25.1% and qRL8.2 with a peak at RM547 on chromosome 8 with an R2 value of 38.8% and QTL qRL8.2 with a peak at RM547 on chromosome 8 with an R2 value of 19.8% and qRL9.1 with a peak at RM524 on chromosome 9 with an R2 value of 31.6% were identified associated with root length in the mapping populations MASARB25 × Pusa Basmati 1460 and HKR47 × MAS26, respectively. qRL8.2 was identified in both populations.
In addition, QTL qGY8.1 with RM339 on chromosome 8 being the peak marker with an R2 value of 34.0% and qGY2.1 on chromosome 2 with RM7562 being the peak marker with an R2 value of 22.8% (Figure 3) and QTL qGY2.2 on chromosome 2 with RM526 being the peak marker with an R2 value of 43.2% (Figure 4) were identified for grain yield under aerobic conditions in the mapping populations MASARB25 × Pusa Basmati 1460 and HKR47 × MAS26, respectively.
Figure 3. QTL likelihood curves of LOD scores for grain yield showing significant regions within qGY2.1 and qGY8.1 under aerobic conditions in MASARB25 × Pusa Basmati 1460.
Figure 4. QTL likelihood curves of LOD scores for grain yield showing significant regions within qGY2.2 under aerobic conditions in HKR47 × MAS26.
For the MASARB25 × Pusa Basmati 1460 population on chromosome 2 in a region of 0-10.2 cM, seven QTL were reported to be significantly associated with seven different root traits and agronomically important traits (RV, FRW, FSW, DSW, PH, PL, and GY). These coexisting chromosomal regions governing root and yield traits provide a unique opportunity for breeders to introgress such regions together as a unit into lowland varieties through MAS/MAB and to develop cultivars possessing higher yield and increased adaptation to aerobic conditions.
In HKR47 × MAS26-derived mapping population, RM256 on chromosome 8, RM28048 on chromosome 12 and RM258 on chromosome 10, and RM3331 on chromosome 8 showed significant association with plant height, panicle length, and number of tillers per plant, respectively (Table 5). RM339 on chromosome 8 and RM204 on chromosome 6 showed significant association with grain yield and dry shoot weight in MASARB25 × Pusa Basmati 1460-derived mapping population. On chromosome 2 in a region of 0-10.2 cM, seven QTL were reported to be significantly associated with seven different root traits and agronomically important traits in the case of MASARB25 × Pusa Basmati 1460 population (Table 4).
In this study, three QTL for root thickness (chromosomes 1, 5, and 6), two for root volume (chromosomes 2 and 6), three for fresh root weight (chromosomes 2 and 6), three for dry root weight (chromosomes 2 and 8), three for fresh shoot weight (chromosomes 2 and 6), two for dry shoot weight (chromosomes 2 and 6), six for plant height (chromosomes 2, 6, 8, 9, and 11), seven for panicle length (chromosomes 2, 8, 9, 10, 11, and 12), five for tiller number (chromosomes 1, 2, and 8), four for number of panicles per plant (chromosomes 2 and 8), three for number of seeds per panicle (chromosomes 8 and 11), one for length/breadth ratio (chromosome 11), and three for grain yield (chromosomes 2 and 8) were identified (Tables 4 and 5). qGY2.1 and qGY8.1, individually and combined, and qGY2.2 exhibited a grain yield improvement of 24%, 28%, 37%, and 36%, respectively, in 2011 and 18%, 25%, 39%, and 24%, respectively, in 2012 under aerobic conditions (Table 6).
Table 6. Yield improvement of lines possessing QTL (QTL+) for grain yield under aerobic conditions over lines not possessing QTL (QTL-) for the two populations
QTL mapping of drought resistance traits in the background of locally adapted indica/Basmati rice lines is reported recently from China and India . However, the identification of QTL linked to yield under aerobic conditions in target populations of an environment is critical. Marker-assisted selection using consistent-effect QTL is an efficient approach for developing appropriate aerobic rice varieties. Trait selection is another important concern in molecular breeding for aerobic traits. As in the case of drought, QTL for a number of physiological and morphological traits have been identified but they have limited implications in breeding drought-tolerant rice varieties. Large-effect QTL for grain yield under drought have been recently identified [21,22] and their successful introgression has established a yield advantage under drought . For developing aerobic rice with high yield potential, three QTL with a large effect on grain yield under aerobic conditions have been identified in our study.
The study also identified large-effect QTL for different root traits that may increase plant ability to uptake nutrients under aerobic conditions. Reduced nutrient uptake, especially of nitrogen and phosphorus under aerobic conditions vis-à-vis flooded conditions, has been the most important factor for lower yield in aerobic systems than in flooded systems of rice cultivation. The identification of such root-related traits and QTL associated with these root traits that increase nutrient uptake under aerobic conditions can help develop aerobic rice varieties with high yield potential.
The root system is considered as one important component for solving the problem of water scarcity. The improvement of upland rice through a deeper root system is considered to be a promising way to increase water uptake, reduce lower canopy temperature, increase stomatal conductance, and ultimately increase grain yield under water-stress conditions. A deeper and thickened root system having large xylem vessels has been shown to allow upland rice varieties to extract more water from the soil, resulting in a higher yield potential under water-scarce conditions . Such root traits are highly likely to have a positive effect on increasing the adaptation of rice genotypes to aerobic conditions. Growth of the rice root, in terms of total dry matter, maximum root depth, root length density, root number, and root volume, increases until flowering stage and then decreases sharply to maturity . Kawata and Soejima  indicated that roots produced after flowering may play an important role during the grain-filling period that ultimately leads to an increase in yield. Root traits are critical for increasing yield under soil-related stresses.
Correlation analysis was carried out to identify how root morphological characters influence the grain yield and yield morphological traits under aerobic conditions. Our study identified a significant and positive correlation between some of the root traits (root number, root volume, root thickness and dry root weight) and yield under aerobic conditions, indicating the role of root traits for improving yield under aerobic situations possibly through improved water and nutrient uptake. This shows that a well developed root system will help the plant in maintaining high plant water status which ultimately leads to increase in yield potential under water deficit conditions.
Correlated characters are of interest for three main reasons: in connection with the genetic causes of correlation through the pleiotropic action of genes, to know how selection for one character will cause a simultaneous change in other characters, and to determine the relationship between character and fitness. In classical quantitative genetics, trait correlations are attributed to the effect of pleiotropy or very close linkage of genes  but unwanted traits might also be selected during MAS because of the co-localization of QTL. The magnitude and direction of influence of these loci on the different phenotypes will markedly affect the utility of such loci in selection for simultaneous improvement of these traits . In our study, on chromosome 2 in a region of 0-10.2 cM, seven QTL were reported to be significantly associated with seven different root traits and agronomically important traits (RV, FRW, FSW, DSW, PH, PL, and GY) in the case of MASARB25 × Pusa Basmati 1460 population (Figure 5). Two QTL, qGY2.1 in the 10-cM region between RM7562 and RM279 on chromosome 2 in MASARB25 × Pusa Basmati 1460 population and qGY2.2 in the 118.5-cM region between RM475 and RM526 on chromosome 2 in HKR47 × MAS26, for grain yield under aerobic conditions were identified. qGY2.1 and qGY8.1 in MASARB25 × Pusa Basmati 1460 showed a stable high over two different years and combined over two years individually, and an increased stable effect over two different years and combined over two years (Table 6). Similarly, in HKR47 × MAS26, qGY2.2 showed a stable high effect over individual years and combined over two years. These QTL were found to be adjacent to the earlier reported QTL qDTY2.2 and qDTY2.3[21,36] for grain yield under drought stress. Although the effect of the two identified QTL needs to be validated in different genetic backgrounds, under the present scenario, successful introgression of the identified QTL following marker-assisted backcrossing can be used to improve the popular varieties Pusa Basmati 1460 and HKR47.
Figure 5. QTL likelihood curves of LOD score of a segment of 10.2 cM on chromosome 2 in MASARB25 × Pusa Basmati 1460 population.
Similarly, QTL for root length within the qGY2.1 region  and for root thickness adjacent to qGY2.2 have been reported. Coexisting chromosomal regions/loci governing different traits provide a unique opportunity for breeders to introgress such regions together as a unit into high-yielding lowland varieties through MAS/MAB and to develop cultivars possessing increased adaptation to aerobic conditions.
Two QTL, qRL8.2 in both mapping populations and qRL8.1 in MASARB25 × Pusa Basmati 1460, for root length have been reported in our investigation. Qu et al.  reported QTL for fresh root weight and root number coinciding with qRL8.2 linked to marker RM331. Studies also compiled numerous earlier reports of QTL for root length , root thickness, and root number  adjacent to or coinciding with qRL9.1. In HKR47 × MAS26 mapping population, qRT1.1 was reported for root thickness with an LOD of 4.4 and R2 value of 38.4%. These QTL may confer a grain yield advantage under direct-seeded conditions and this is supported by the earlier reported large-effect consistent QTL (qDTY1.1) for grain yield under drought  in the region of qRT1.1.
Under dry direct-seeded conditions, the QTL peak (qGY8.1) for grain yield was seen in MASARB25 × Pusa Basmati 1460 at RM339 with an R2 value of 34%. At the same position and linked to the same marker (RM339), qDTY8.1 in Basmati334/Swarna was reported by Vikram et al.  for grain yield under drought. Hanamaratti,  reported RM339 on chromosome 8 associated with relative yield and drought susceptibility index in IR64 × Binam-derived NILs under drought stress. Adjacent to the QTL qGY8.1 within a region of 14.4 cM, other QTL (qTN8.2, qNPP8.2, and qRL8.2) have been reported (Figure 6). In a segment of 21.7 cM on chromosome 8 (101.5-123.2 cM), QTL for DRW, TN, NPP, PL, NSP, PH, and GY were also reported in our study. The effect of these regions on a number of traits that is likely to impart improved aerobic adaptation. GY strongly suggests the presence of more than one gene within these QTL affecting a wide range of traits. These genes conferring a GY advantage under aerobic adaptation may have undergone strong natural selection to stay together in the course of evolution. These sub-QTL with a discernible phenotypic effect on GY may affect the same/different physiological traits in response to different severities of stress, leading to a GY advantage.
Figure 6. QTL likelihood curves of LOD score of coexisting qtl on chromosome 8 in MASARB25 × Pusa Basmati 1460 population.
QTL qPL12.1 (RM28048) identified in our study in HKR47 × MAS26 mapping population showed an effect on a yield-attributing trait (i.e., panicle length) under dry direct-seeded conditions. Bernier et al. [21,45] reported QTL qtl12.1 linked to the same marker. Bernier et al.  and Dixit et al.  showed a consistent large effect on grain yield under drought over a wide range of environments in the Philippines and eastern India. This confirms that this QTL is effective not only at IRRI, where it was initially detected, but also in at least a part of the target environment of eastern and northern India, where improved drought resistance in upland rice cultivars is an important breeding objective. This is also supported by the study of Bernier et al.  that mentioned that this locus is associated with a low QTL × environment interaction under severely stressed conditions, which is one of the two major requirements for the use of a QTL in MAS. The QTL for biomass, harvest index, days to flower, plant height, flowering delay, drought response index, and panicle number under stress were mapped in the same region. qtl12.1 also influences water uptake under upland stress and increases harvest index .
In our study, we reported a 25.1-cM segment between RM589 and RM314 on chromosome 6 affecting different root (RV, RT, and FRW) and shoot (FSW and DSW) traits under aerobic conditions in the case of MASARB25 × Pusa Basmati 1460 mapping population (Figure 7). In this region, a large-effect QTL (qDTY6.1) associated with grain yield in aerobic environments was identified in a total of 20 hydrological environments over a period of five seasons and in five populations in three genetic backgrounds using bulk-segregant analysis . Co-localization of qDTY6.1 with a region identified to govern several root and shoot traits is an implication that this is an important region for improving yield under aerobic conditions.
Figure 7. QTL likelihood curves of LOD score of coexisting qtl on chromosome 6 in MASARB25 × Pusa Basmati 1460 population.
Our study identified a number of lines with higher yield under aerobic conditions, higher root length and dry biomass, GY, and length/breadth ratio, and with Pusa Basmati 1460-specific alleles in homozygous or heterozygous condition at the BAD2 locus (Table 7). These lines will serve as novel materials for the selection of stable aerobic Basmati rice varieties.
Table 7. Values of mean yield and mean root traits of selected promising lines under aerobic conditions
Our study identified a total of 35 QTL associated with 14 traits on chromosomes 1, 2, 5, 6, 8, 9, and 11 in MASARB25 × Pusa Basmati 1460 population and 14 QTL associated with nine traits on chromosomes 1, 2, 8, 9, 10, 11, and 12 in HKR47 × MAS26-derived population. These identified QTL included three large-effect stable QTL for increased yield under aerobic conditions and QTL for several root-related traits likely to increase water and nutrient uptake under aerobic conditions. A mechanism associated with higher yield of promising lines under dry direct-seeded conditions, indicating relocation of resources during grain filling, is suggested by coexisting QTL for root and yield-attributing traits. Results from our study suggest that reported QTL are complex loci where multiple genes may be working independently or in coordination with each other, leading to an increase in GY under drought. Coexisting chromosomal regions/loci governing different traits for aerobic adaptation provide a unique opportunity for breeders to introgress such regions together as a unit into high-yielding drought-susceptible varieties through MAS/MAB and to develop cultivars possessing increased tolerance of varying stress severities.
cm: Centimeter; CTAB: Cetyltrimethyl ammonium bromide; DNA: Deoxyribonucleic acid; DPE: Direction of phenotypic effect; DRW: Dry root weight; DSW: Dry shoot weight; FRW: Fresh root weight; FSW: Fresh shoot weight; g: Gram; 1,000 GW: 1,000-grain weight; GY: Grain yield; IRRI: International Rice Research Institute; Kg ha-1: Kilogram per hectare; L/B: Length/breadth ratio; LOD: Logarithm of odds; LR: Likelihood ratio; MAB: Marker-assisted breeding; MAS: Marker-assisted selection; ml: Milliliter; mm: Millimeter; NPK: Nitrogen, phosphorus, and potassium; P/P: Number of panicles plant-1; PL: Panicle length; PH: Plant height; PAGE: Polyacrylamide gel electrophorosis; PCR: Polymerase chain reaction; RIL: Recombinant inbred line; R2: Percent phenotypic variance; RL: Root length; RN: Root number, RV, Root volume; RT: Root thickness; S/P: Seeds panicle-1; TN: Effective number of tillers plant-1.
In the past five years have you received reimbursements, fees, funding, or salary from an organization that may in any way gain or lose financially from the publication of this manuscript, either now or in the future? Is such an organization financing this manuscript (including the article-processing charge)? If so, please specify.–No
Do you hold any stocks or shares in an organization that may in any way gain or lose financially from the publication of this manuscript, either now or in the future? If so, please specify.–No
Do you hold or are you currently applying for any patents relating to the content of the manuscript? Have you received reimbursements, fees, funding, or salary from an organization that holds or has applied for patents relating to the content of the manuscript? If so, please specify.–No
Do you have any other financial competing interests? If so, please specify.–No
Non-financial competing interests
Are there any non-financial competing interests (political, personal, religious, ideological, academic, intellectual, commercial or any other) to declare in relation to this manuscript? If so, please specify.–No
NS was involved in the conception of the experiment, analysis, interpretation of the data, and drafting the article and final approval of the version to be published; SJ was involved with conception and design of the experiment, revising manuscript content critically and final approval of the version to be published; AK was involved with the analysis, interpretation of the data, in the critical revision of the manuscript and final approval of the version to be published; BSM was involved in the conception of the experiment, in the critical revision of the manuscript and final approval of the version to be published; RKJ was involved in the design of the experiment, the critical revision of the manuscript and final approval of the version to be published.
The authors thank the Department of Science and Technology, New Delhi, India (DST No. SR/SO/PS/0041/2009), and Monsanto Beachell Borlaug International Scholarship Program (MBBISP) for providing financial support for this study.
What is the challenge program on water and food.
Barker R, Dawe D, Tuong TP, Bhuiyan SI, Guerra LC: The outlook for water resources in the year 2020: challenges for research on water management in rice production. In “Assessment and orientation towards the 21st Century”, 7-9 September 1998. Cairo, Egypt: FAO; 1998:96-109. [Proceedings of 19th Session of the International Rice Commission]
Plant Prod Sci 2005, 8:231-241. Publisher Full Text
Tuong TP, Bouman BAM: Rice production in water scarce environments. In “Water productivity in agriculture: limits and opportunities for improvement”. Edited by Kijne JW, Barker R, Molden D. Wallingford, UK: CABI Publishing; 2003:53-67.
Agric Water Manage 2007, 92:81-89. Publisher Full Text
Field Crops Res 2006, 97:43-52. Publisher Full Text
Field Crops Res 2006, 96:252-259. Publisher Full Text
Rabello AR, Guimaraes CM, Rangel PHN, Silva FR, Seixas D, Souza E, Brasileiro ACM, Spehar CR, Ferreira ME, Mehta M: Identification of drought responsive genes in roots of upland rice (Oryza sativa L).
Field Crops Res 2011, 120:205-214. Publisher Full Text
Theor Appl Genet 1997, 95:143-152. Publisher Full Text
Rice 2009, 2:115-128. Publisher Full Text
Zhang J, Zheng HG, Aarti A, Pantuwan G, Nguyen TT, Tripathy JN, Sarial AK, Robin S, Babu RC, Nguyen BD, Sarkarung S, Blum A, Nguyen NT: Locating genomic regions associated with components of drought resistance in rice: comparative mapping within and across species.
Theor Appl Genet 2001, 103:19-29. Publisher Full Text
Plant Prod Sci 2011, 14:1-14. Publisher Full Text
Crop Sci 2007, 47:507-516. Publisher Full Text
Euphytica 2009, 167:261-269. Publisher Full Text
Euphytica 2006, 152:259-273. Publisher Full Text
Gomez SM, Bhoopati NM, Kumar SS, Ramasubramaniam T, Chengsong Z, Jayaprakash P, Senthil A, Babu RC: Molecular mapping and location of QTL for drought resistance traits in Indica rice (Oryza sativa L.) lines adapted to target environments.
Acta Physiol Plant 2010, 32:355-364. Publisher Full Text
Crop Sci 1997, 37:1426-1434. Publisher Full Text
Proc Crop Sci Soc Jpn 1974, 43:354-374. Publisher Full Text
Hittalmani S, Huang N, Courtious B, Venuprasad R, Shashidar HE, Bagali GG, Li ZK, Zhuang JY, Zheng KL, Liu GF, Wang GC, Singh VP, Sidhu JS, Srivantaneeyakul S, McLaren G, Khush GS: Identification of QTL for growth and grain yield related traits in rice across nine locations in Asia.
Dixit S, Swamy BPM, Vikram P, Ahmed HU, Sta Cruz MT, Amante M, Atri D, Leung H, Kumar A: Fine mapping of QTL for rice grain yield under drought reveals sub-QTL conferring a response to variable drought severities.
MacMillan K, Emrich K, Piepho HP, Mullins CE, Price AH: Assessing the importance of genotype x environmental interaction for root traits in rice using a mapping population and conventional QTL analysis.
Theor Appl Genet 1995, 90:969-981. PubMed Abstract
Mol Breed 2010, 26:533-538. Publisher Full Text
Ghimire KH, Quiatchon LA, Vikram P, Swamy BPM, Dixit S, Ahmed HU, Hernandez JE, Borromeo TH, Kumar A: Identification and mapping of a QTL (qDTY1.1) with a consistent effect on grain yield under drought.
Vikram P, Mallikarjuna Swami BP, Dixit S, Helaluddin A, Sta Cruz MT, Singh AK, Guoyou Y, Kumar A: Bulk segregant analysis: an effective approach for mapping consistent-effect drought grain yield QTL in rice.
Hanamaratti NG: Identification of QTL for physiological and productivity traits under drought stress and stability analysis in upland rice (Oryza sativa L.). Dharwad: University of Agricultural Sciences; 2007.
Bernier J, Kumar A, Spaner D, Verulkar S, Mandal NP, Sinha PK, Peeraju P, Dongre PR, Mahto RN, Atlin GN: Characterization of the effect of rice drought tolerance qtl12.1 over a range of environments in Philippines and Eastern India.
Euphytica 2009, 166:207-217. Publisher Full Text