Physiological Response of Advanced Lines of Rice under Water Deficit in Campeche, Mexico ()
1. Introduction
Cultivation of rice (Oryza sativa L.) in Mexico occupies the fourth place after corn, bean and wheat in terms of surface, production and consumption, and is grown in an approximate area of 42,310 ha with a national average grain yield of 5.8 t·ha−1 [1]. Of this surface, 75% is cultivated under conditions of irrigation or permanent flooding during the cycle, with an average yield of 6.4 t·ha−1 and the remaining 25% is cultivated under rainfed conditions with average yields of 3.8 t·ha−1 [1]. However, since 2001, the national production has decreased as a result of the disarticulation of the rice chain, which is why there is use for importing it, and 85% of the national consumption of thin grain rice was imported in 2015 [2] [3] [4]. Therefore, the National Institute of Forest, Agricultural and Livestock Research (Instituto Nacional de Investigaciones Forestales Agrícolas y Pecuarias, INIFAP) in Mexico and the Latin American Fund for Irrigated Rice (Fondo Latinoamericano de Arroz de Riego, FLAR) in Colombia work together for the stockpiling, selection and biometric evaluation of materials of rice with thing long grain under the environmental conditions of the dry and humid tropics in Mexico [4].
Sowing rainfed rice is carried out in the humid tropics of the southeastern region of Mexico, which stretches from the Papaloapan Basin in part of the states of Veracruz and Oaxaca, the region of La Chontalpa (Tabasco) to the areas of Palizada and Yohaltun in the state of Campeche, zone considered as vulnerable to climate change due to the frequent losses of rice harvests caused by the “intraestival” summer drought and/or torrential rains in the autumn that cause flooding and lodging of rice crops [5]. Drought is considered as a deficit of sufficient water availability to cause a decrease in the yield of crops [6]. On the other hand, the magnitude of the decrease in the grain yield depends on the duration of the drought and severity of stress [7], in addition to the growth stage of the crop [8]. Thus, the sensitivity of rice (Oryza sativa L.) to drought increases, when the drought takes place in the flowering period [9]. As consequence, the water content is considerably reduced in the reproductive stage in plant cells, showing symptoms of withering and loss of turgidity, which has an effect on a decrease of the gas exchange and photosynthesis of the plant and, finally, on the grain yield [10].
The studies in advanced lines of rice from FLAR are directed toward parameters of stability of yields [11]. However, the low yields are related to climate change, and studies of the materials with different abiotic stresses are required for this, in this case, water stress to guarantee availability of materials with resistance or tolerance to drought [12]. Therefore, the objective of this study was to research the response of eight advanced lines and a control variety of rice (Oryza sativa L.), on grain yield and its components, leaf area, total transpiration and transpiration efficiency under irrigation (I) and drought (D) in a greenhouse, with the aim of selecting high yield materials with tolerance to drought and making recommendations to farmers of the region.
2. Materials and Methods
2.1. Soil and Climate in the Study Area
The experiment was established in the 2015 autumn-winter cycle in a greenhouse of Colegio de Postgraduados, Campus Campeche, five meters high with a white mesh roof, however, to avoid the incidence of rain, a polyethylene cover was placed at a height of three meters. To allow the plants to be exposed to room temperature, air was allowed into the greenhouse, lifting the mesh at the ends. Colegio de Postgraduados, Campus Campeche, located in the Haltunchén-Edzná road, km 17.5, Sihochac, municipality of Champotón, Campeche, Mexico (17˚49'N, 91˚08' and 27 masl). The locality of Sihochac has warm subhumid type climate with summer rains [Awo (w)], with a mean precipitation of 1500 mm and mean annual temperature of 26˚C, with maximum temperatures of 42˚C and minimum of 10˚C [13]. Soil analysis was carried out in the Soil Laboratory of Universidad Autónoma Chapingo (Table 1).
2.2. Genetic Material in Evaluation
Eight advanced lines (F6) of rice and a control variety (El Silverio) were included. The lines presented characteristics of high yield and resistance to the vector of the white leaf virus (Tagosodes oryzicolus) and rice “burn” (Magnaporthe oryzae, anamorphic Pyricularia oryzae), higher quality of grain, defined by its amylose content, appearance of the processed grain, high recovery of full grains and tolerance to the delay in harvest [14]. The lines came from the nurseries of the Latin American Fund for Irrigated Rice (FLAR, Palmira, Colombia). The genetic material was preselected by grain yield in humid, subhumid and dry tropical environments in the spring-summer cycle (S-S) of 2012 by researchers from INIFAP in different localities (Palizada, Campeche; Zacatepec, Morelos; Cárdenas, Tabasco; and Tecomán, Colima, Mexico). The control variety was chosen for the following characteristics: The Silverio is a medium grain variety with 10% of “chalkiness” and highly versatile, it can be grown under rainfed conditions or rainfed with auxiliary irrigation in the south-east humid tropic (Veracruz, Oaxaca, Tabasco, Campeche and Chiapas) and sub-humid northeast (Tamaulipas), also under irrigation in the Pacific dry tropic (Nayarit, Jalisco, Colima and Michoacán). Its average yield potential is around 6 t·ha−1 under rainfed conditions, of 7 t·ha−1 in rainfed with auxiliary irrigation and 8 t·ha−1 under irrigation. Due to its growth cycle characteristics and intermediate height (130 to 135 days and 91 cm), is well suited to rainfed conditions [4] (Table 2).
2.3. Study Factors and Experimental Design
The experiment consisted in two treatments of soil moisture; irrigation (I) that consisted in maintaining the soil moisture close to field capacity (FC) and drought (D) with the suspension of irrigation at 65 days after sowing (das) and until physiological maturity. The day of onset of drought was chosen based on
Table 1. Physical and chemical characteristics of the soil used in the experiment to study the physiological response of advanced lines of rice to water stress. Autumn-winter cycle 2015. Campeche, Mexico.
The texture was determined, Bouyoucos hydrometer method; AD: apparent density, test tube method (t·m−3); FC: field capacity, pressure cooker method equivalent to 0.3 atm PPM: percentage of permanent withering, pressure membrane method equivalent to 15 atm; pH, electronic pH gauge; EC: electric conductivity, bridge method of electric conductivity; OM: organic matter, Walkley-Black method; N: total nitrogen, vapor dragging; P: phosphorus, Olsen method; K: potassium, flame emission spectrophotometry method.
Table 2. Genetic material of rice preselected by yield, used in the experiment to study the physiological response of advanced lines of rice to water stress. Autumn-winter 2015 cycle. Campeche, Mexico.
The lines come from the nurseries of the Latin American Fund for Irrigated Rice (FLAR), Palmira, Colombia.
the fact that the reproductive stage of the crop begins at 65 days. Sowing was done on February 12, 2015, with a rice seed of each variety per experimental unit, with a fertilization dose of 180-40-40; a third of the nitrogen and all the phosphorus and potassium were applied during sowing; a second third of the nitrogen at the beginning of the stage of clustering, 25 days after sowing (das), and the rest of the nitrogen during filling, 78 days after sowing (das), with urea as a nitrogen source, triple calcium superphosphate and potassium chloride as sources of phosphorus and potassium, respectively [14].
A complete random block design with three repetitions was used in irrigation and three in drought; the experimental unit consisted in an individual plant per PVC tube of 1 m high and 10 cm of diameter. A cylindrical polyethylene bag of the same dimensions was placed in each PVC tube and it was filled with previously sterilized soil; the soil was taken to FC and the initial weight at FC was defined (IWFC) in each experimental unit. In the irrigation treatment (I) the tubes were weighed every third day, to calculate the amount of water lost from evapotranspiration (ET) and add the water required, to take the soil from each tube to IWFC and maintain the level of moisture close to FC from sowing to physiological maturity of plants; the same procedure was done in the drought treatment (D) as in irrigation; water was added to the tubes, to take the soil to FC and the soil water content was maintained close to FC until 65 days after sowing (das); from this date the application of water was suspended and the weight was recorded only at 114 das, date when recovery irrigation was applied. The amount of ET (g of water) in each tube was calculated as the difference between initial weights of the tube to FC (kg) − weight of the tube in each date of weighing (kg). The transpiration per plant (T) was calculated as the difference between ET − direct evaporation from the soil (ES) estimated in control tubes without plant in each weighing; the total accumulated transpiration was calculated when adding the T obtained in each weighing from sowing to physiological maturity. The plants were kept under a polyethylene cover at a height of 3.0 m, to avoid the incidence of rain on the experiment and allow the plants to be exposed to the environmental temperature.
2.4. Study Variables and Statistical Analysis
The following variables were evaluated: grain yield (GY, g) which was obtained when determining the weight of the grain after dekernelling of all the panicles of the plants harvested; the final aerial biomass (FABM, g) was defined when obtaining the total dry weight of the aerial part of the plant when reaching physiological maturity; the harvest index (HI) was calculated as the quotient between grain yield (GY, g) and final aerial biomass (g):
(HI = GY/BM) (2.4)
the number of panicles per plant (NPP) was determined when counting the number of panicles in each experimental unit; the number of grains per panicle (NGP−1) was determined when counting the number of grains in each panicle per plant and averaging the number of grains between the total panicles per plant; the weight of 1000 grains (WTG, g) was obtained when weighing 1000 grains taken randomly from each experimental unit of the sample used to determine the grain yield; the days until flowering or anthesis (A) was recorded when 50% of the panicles of each experimental unit showed spikelets with exposed anthers and in dehiscence; the days until physiological maturity (PM) were found when 80% of the plants in each experimental unit had lost the green color of their aerial organs and presented a characteristic hay color of each genotype; the plant height (PH, cm) was determined by measuring the height from the soil surface to the terminal spikelet in each experimental unit in the stage of physiological maturity; for the green leaf area (GLA, cm2), the area (cm2) of each leaf from the main stalk of each experimental unit was determined by multiplying the length (cm) × the width (cm) × 0.75, every five days from 18 to 109 days after sowing (das); the total transpiration per plant (TT, kg) was calculated when adding the transpiration determined between consecutive measurements during the biological cycle of each plant; and the transpiration efficiency to produce total biomass of the plant (ETBMT, g of BMT/kg of H2O transpired) was calculated by dividing the total biomass between the amount of water transpired by the plant at the end of the biological cycle.
(ETBMT = g of BMT/kg of H2O transpired) (2.5)
ANOVA was carried out with the data for all the variables with the SAS software, Version 9.1, for Windows [15] individually for each treatment of soil moisture and in a combined way, to determine the significance of genotype x environment interaction. The means comparison was carried out with DMS (p ≤ 0.05).
3. Results
3.1. Soil Moisture Content
The soil moisture content in the irrigation treatment was kept close to field capacity (FC), during the entire experiment; in the drought treatment the soil water content was kept close to FC until beginning the drought treatment at 65 days after sowing (das). In drought the soil water content reached values under the permanent wilting Point (PWP), as the cultivation cycle happened (Figure 1).
3.2. Grain Yield and Its Components
In irrigation, the statistical analysis of the data detected highly significant differences (p ≤ 0.01) between genotypes for grain yield (GY), final aerial biomass (FABM), harvest index (HI), number of panicles per plant (NPP), number of grains per panicle (NGP), weight of 1000 grains (WTG), days to anthesis (DA) and plant height (PH) (Table 3). The grain yield and its components in irrigation presented high variability between genotypes. The greatest variation was observed in the grain yield, final aerial biomass, number of grains per panicle and plant height. Genotypes 2 and 7 presented the highest grain yield in addition to final aerial biomass, harvest index, number of panicles per plant, weight of one thousand grains, and plant height (Table 3).
In drought the statistical analysis of the data detected highly significant differences (p ≤ 0.01) between genotypes for grain yield (GY), final aerial biomass (FABM), harvest index (HI), number of panicles per plant (NPP), number of grains per panicle (NGP), weight of 1000 grains (WTG), days to anthesis (DA), days to physiological maturity (DPM) and plant height (Table 4). Genotypes 2 and 7 presented the highest grain yield accompanied by the harvest index, number of panicles per plant, and plant height. Genotype 2 presented a higher final aerial biomass, number of grains per panicle and weight of one thousand grains, compared to the rest of the genotypes (Table 4).
The combined statistical analysis of irrigation and drought detected highly significant differences (p ≤ 0.01) between levels of humidity, between genotypes and interaction levels of humidity-genotypes in all the variables assessed (Table 5). Drought reduced grain yield, final aerial biomass, harvest index, number of panicles per plant, number of grains per panicle, weight of one thousand grains, and plant height (Table 5). The reduction resulting from the water stress effect was higher in the final aerial biomass (52%), grain yield (45%), number of panicles per plant (36%), weight of one thousand grains (31%), plant height (28%), number of grains per panicle (27%), and harvest index (15%) (Table 5).
Figure 1. Water content of the soil in irrigation and drought from 8 to 120 days after sowing (das), under greenhouse conditions. Field Capacity (FC), Permanent Wilting Point (PWP) Autumn-winter cycle 2015. Campeche, Mexico.
Table 3. Grain yield and its components, days to anthesis and plant height of rice genotypes under irrigation. Autumn-winter cycle 2015. Campeche, Mexico.
Means with different letter in each column are statistically different (DMS, P ≤ 0.05). GY: grain yield per plant; FABM: final aerial biomass; HI: harvest index; NPP: number of panicles per plant; NGP: number of grains per panicle; WTG: weight of one thousand grains; DA: days until flowering and and PH: plant height. *The genotypes come from nurseries of the Latin American Fund for Irrigated Rice (Fondo Latinoamericano para Arroz de Riego, FLAR), Palmira, Colombia.
Table 4. Grain yield and its components, days to flowering, days to physiological maturity and plant height of rice genotypes in drought. Autumn-winter cycle 2015. Campeche, Mexico.
Means with different letter in each column are statistically different (DMS, P ≤ 0.05). GY: grain yield per plant; FABM: final aerial biomass; HI: harvest index; NPP: number of panicles per plant; NGP: number of grains per panicle; WTG: weight of one thousand grains; DA: days until flowering; DPM: days until physiological maturity and PH: plant height. *The genotypes come from nurseries of the Latin American Fund for Irrigated Rice (Fondo Latinoamericano para Arroz de Riego, FLAR), Palmira, Colombia.
Table 5. Grain yield and its components: days until flowering, days until physiological maturity and plant height, average of rice genotypes in irrigation and drought. Autumn-winter cycle 2015. Sihochac, Champotón, Campeche-Mexico.
Means with different letter in each column are statistically different (DMS, P ≤ 0.05). GY: grain yield (g); FABM: final aerial biomass (g); HI: harvest index (%); NP P-1: number of panicles per plant-1; NGP-1: number of grains per panicle-1; WTG: weight of one thousand grains; DA: days until anthesis; DPM: days until physiological maturity; and PH: plant height.
3.3. Leaf Area
In irrigation, the statistical analysis detected highly significant differences (p ≤ 0.01) between genotypes for the leaf area and statistical differences (p ≤ 0.05) at 25 days after sowing (das) (Figure 2). Figure 2 shows the trend of the genotypes in irrigation on the increase of leaf area throughout development to later decrease because of the withering of lower leaves. In average of all the measurements carried out, genotypes 1, 2 and 7 were detected with the largest leaf area; on the contrary, genotypes 3, 6, 8 and 9 presented the smallest leaf area.
Figure 2. Green leaf area of rice genotypes, in irrigation from 18 to 109 days after sowing (das). Autumn-winter cycle 2015. Campeche, Mexico. **(P ≤ 0.01), *(P ≤ 0.05), ns = Not significant. *The genotypes come from the nurseries of the Latin American Fund for Irrigated Rice (Fondo Latinoamericano para Arroz de Riego, FLAR), Palmira, Colombia.
In drought the statistical analysis detected that highly significant differences (p ≤ 0.01) were present between genotypes for all the dates of measurement except at 25 das, when no significant differences were found, and at 46 das when significant differences between genotypes were detected (p ≤ 0.05) (Figure 3). Figure 3 shows the trend of the genotypes in the leaf area; it is observed that starting at 65 das, the leaf area decreases due to the effect of water stress. The average of all the measurements carried out was detected on genotype 1 with the highest values in leaf area, and genotypes 8 and 9 with the lowest values of leaf area.
In the combined statistical analysis of irrigation and drought, highly significant differences were detected between moisture levels of the soil at 32 and 53 das, and from 67 to 109 das. Highly significant differences (p ≤ 0.01) were detected between genotypes in all the measurement dates, except at 25 das when significant differences (p ≤ 0.05) were found. The leaf area under irrigation was higher in irrigation than in drought, because starting at 65 das the water deficit began at the time of suspending water in the plants in drought treatment (Figure 4).
3.4. Total Transpiration and Transpiration Efficiency of the Plant
In irrigation, the statistical analysis of the data detected highly significant differences (p ≤ 0.01) between genotypes for transpiration efficiency and total transpiration, average of all the measurements performed. In drought the statistical analysis of the data detected significant differences (p ≤ 0.05) between genotypes for transpiration efficiency and no differences were found for total transpiration per plant. The transpiration efficiency in irrigation was lower (5.3 g·kg−1 H2O p−1), compared to drought (9.8 g·kg−1 H2O p−1), genotypes 2, 6 and 7 presented the highest TE compared to the rest of the genotypes and in drought it was genotypes 2 and 7 that outperformed the rest of the genotypes (Table 6). For total transpiration, it was higher in irrigation (8.1 kg) than in drought (3.8 kg); in irrigation genotypes 7 and 8 showed the highest values; no significant differences were found in drought (Table 6).
In the combined analysis of irrigation and drought, highly significant differences (p ≤ 0.01) were detected between environments of humidity and between genotypes for transpiration efficiency and total transpiration per plant; highly significant differences were found in the interaction of humidity environments-genotypes for total transpiration per plant and no significant differences were found for transpiration efficiency (Table 7).
Figure 3. Green leaf area of rice genotypes, in drought from 18 to 109 days after sowing (das). Autumn-winter cycle 2015. Campeche, Mexico. **(P ≤ 0.01), *(P ≤ 0.05), ns = Not significant. *The genotypes come from the nurseries of the Latin American Fund for Irrigated Rice (Fondo Latinoamericano para Arroz de Riego, FLAR), Palmira, Colombia.
Figure 4. Average green leaf area of nine rice genotypes, in irrigation and drought from 18 to 109 days after sowing (das). Autumn-winter cycle 2015. Campeche, Mexico. **(P ≤ 0.01), *(P ≤ 0.05), ns = Not significant. *The genotypes come from the nurseries of the Latin American Fund for Irrigated Rice (Fondo Latinoamericano para Arroz de Riego, FLAR), Palmira, Colombia.
Table 6. Transpiration efficiency to produce total biomass (ETBMT) and total transpiration (TT) of rice genotypes in irrigation and drought. Autumn-winter cycle 2015. Campeche, Mexico.
Means with different letter in each column are statistically different (DMS, P ≤ 0.05). *The genotypes come from the nurseries of the Latin American Fund for Irrigated Rice (Fondo Latinoamericano para Arroz de Riego, FLAR), Palmira, Colombia.
Table 7. Transpiration efficiency to produce total biomass (ETBMT) and total transpiration (TT) of the average plant in rice genotypes for irrigation and drought. Autumn-winter 2015 cycle. Campeche, Mexico.
Means with different letter in each column are statistically different (DMS, P ≤ 0.05). *The genotypes come from the nurseries of the Latin American Fund for Irrigated Rice (Fondo Latinoamericano para Arroz de Riego, FLAR), Palmira, Colombia.
4. Discussion
When comparing the result from nine rice genotypes under conditions of irrigation and drought, it was determined that the genotypes presented high genetic variability in grain yield, biomass, harvest index, number of panicles per plant, number of grains per panicle, and weight of one thousand grains, as well as in days until anthesis, plant height, leaf area, transpiration efficiency and total transpiration, this clearly influenced by the environment, with higher grain yield and its components in irrigation. These results agree with what has been reported by various authors [16] [17] inferring that water stress affected negatively the grain yield in rice cultivation. Of the genotypes studied, lines 2 and 7 stood out both under conditions of irrigation and in drought; it should be mentioned that the rest of the advanced lines from FLAR were superior to the control variety and showed adaptation capacity to the zone. It is important to mention that a deficit in water availability causes a decrease in the yield from crops, which was observed in the treatment under drought; when the drought takes place during the flowering period, the water content of the plant cells is reduced, showing symptoms of withering and loss of turgidity that have an impact on the decrease of the plant’s gas exchange and photosynthesis, and therefore on the grain yield [9] [10] [18] [19]. The variable height is an aspect to consider in the selection of rice materials, in this experiment it was higher in irrigation than in drought and fluctuated between 67 to 87 cm, the height is associated with lodging problems which is why a taller height represents higher risk in the plant’s fall, and what is sought in a material of long and thin grain is a compact height to avoid the loss of the yield from lodging [3] [12]. The least affected variables were days to anthesis (9%) and days to physiological maturity (12%), with a difference of eight and five days until flowering and days until physiological maturity, respectively, between irrigation and drought, with them being higher in drought; it was mentioned that drought before flowering delays the anthesis and this delay reduces the percentage of fertile spikelets and full grains [20]; as a result there is a lower photosynthetic rate, relative water content and stomata conductance, and these characteristics show a close relationship with the efficiency in water use, which is reflected in a reduction in grain yield [21]. The importance of leaf development for the crops depends on the interception of solar radiation, key element for the photosynthetic activity that is of vital importance in biomass production with a considerable contribution to the yield [22]. The leaf area in irrigation was higher than in drought, because starting at 65 days after sowing, the water deficit began when watering was suspended in the plants; this is because the plants generally limit the area and number of leaves in response to the water stress to reduce the consumption of water and avoid losses in the yield [23]. The transpiration efficiency was higher under drought (9.8 g·kg−1) compared to irrigation (5.3 g·kg−1), and this agrees with what was obtained in other studies where an increase in transpiration efficiency is observed under conditions of water deficit [24] [25]. On the contrary, the total transpiration per plant was 8.1 kg of water, higher than the treatment under drought that showed 3.8 kg of water; when the soil moisture begins to decrease, the most effective response of the plant is to reduce transpiration. The first signs of water decrease in the soil are leaf rolling and closing of stomata, and they are the basic mechanism to reduce the impact of drought and are induced first in the vegetative phase and then in the reproductive phase [26] [27].
5. Conclusion
The advanced lines of long grain rice generally presented a higher capacity for adaptation to the producing zone of Campeche outperforming the control variety; the genotypes identified as one and two presented higher grain yield and its components, as well as plant height, leaf area, transpiration efficiency, both in irrigation and in drought, which makes them strong candidates to be used as progenitors in breeding programs and as commercial varieties in irrigation and drought production systems.