Next Article in Journal
Utilization of Palm Oil Midrib Biochar as Soil Amendment with a Newly Isolated Bacillus sp. SM11 for Growth Enhancement and Nitrate Reduction in Romaine Lettuce (Lactuca sativa L. var. longifolia)
Previous Article in Journal
Recombinase Polymerase Amplification Assay for Rapid Field Diagnosis of Stewart’s Wilt of Corn Pathogen Pantoea stewartii subsp. stewartii
Previous Article in Special Issue
Design and Experiment of Double-Row Seed-Metering Device for Buckwheat Large Ridges
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Specialty Rice (Oryza sativa L.) with High and Stable Grain Yield under Rainfed Lowland Conditions

by
Ace Mugssy L. Agustin
1,2,
Jeremias L. Ordonio
1,
Marie Bie S. Natividad
1,
Nonawin B. Lucob-Agustin
3,
Roel R. Suralta
3,
Hiroshi Ehara
4,
Shiro Mitsuya
2 and
Mana Kano-Nakata
4,*
1
Department of Crop Science, College of Agriculture, Central Luzon State University, Science City of Muñoz 3120, Philippines
2
Graduate School of Bioagricultural Sciences, Nagoya University, Nagoya 464-8601, Japan
3
Crop Biotechnology Center, Philippine Rice Research Institute, Science City of Muñoz 3119, Philippines
4
International Center for Research and Education in Agriculture, Nagoya University, Nagoya 464-8601, Japan
*
Author to whom correspondence should be addressed.
Agriculture 2023, 13(10), 1985; https://doi.org/10.3390/agriculture13101985
Submission received: 18 September 2023 / Revised: 8 October 2023 / Accepted: 9 October 2023 / Published: 12 October 2023
(This article belongs to the Special Issue Improvement of the Technology of Cereal Production)

Abstract

:
This study aimed to identify superior genotypes of specialty rice (SR) with comparable or higher grain yield than the drought-tolerant check variety under rainfed and controlled-drought conditions. A total of 17 SR varieties (six aromatic, six pigmented, five glutinous) and a drought-tolerant check variety with ordinary grain quality were evaluated under rainfed lowland and controlled-drought conditions from 2019 to 2021 at Central Luzon State University in the Philippines. Among the SR varieties, aromatic NSIC Rc344, pigmented Black rice, and glutinous NSIC Rc15 had comparable or higher grain yield than the drought-tolerant check variety under both rainfed and controlled-drought conditions. These selected genotypes were classified as the highest yielding, with a more stable yield than the drought-tolerant check variety across the hydrological conditions based on the BLUPs productivity and stability test and drought tolerance indices. The selected SR varieties had a greater panicle number (NSIC Rc344), more grains per panicle (NSIC Rc15), and a higher 1000-grain weight and harvest index (Black rice). In comparison to a higher yield but with a higher market price due to the superior grain quality of the identified SR than the drought-tolerant check variety, the net income in rainfed lowland conditions significantly increased by 69–108%. These results suggest that planting good-performing SR in rainfed lowlands can increase profitability in this ecosystem due to the higher market price compared to ordinary drought-tolerant varieties.

1. Introduction

Rainfed lowland ecosystems occupy around 36% (52 M ha) of the total global rice area, with more than 50% (~27 M ha) vulnerable to drought stress [1], a major limiting factor in this ecosystem. Drought stress is further aggravated by changing climates, specifically the changing rainfall patterns and faster depletion of soil moisture. It directly affects rice productivity by inhibiting growth and development, causing a significant reduction in grain yield [2,3,4]. Drought stress during the reproductive phase is the most detrimental to yield formation [5,6]. This is primarily attributed to its adverse effects on both the spikelet quantity [4] and fertility [5]. The decrease in the spikelet number per panicle arises from reduced panicle branching [4], while the decrease in spikelet fertility is due to reduced panicle exsertion and anther dehiscence [6]. Dry matter production is also adversely affected, mainly due to a reduced photosynthetic rate [7] from the limited supply of CO2 [8], which is the consequence of conserving water loss through stomatal closure [9]. Drought may also affect the implementation of production management activities, such as delayed crop establishment and fertilizer application, as well as the effective control of weeds through water management [10]. Low rice productivity leads to low income for the farmers and, consequently, low input use in the following cropping season [11], a cycle that traps farmers into poverty. Furthermore, the low productivity itself also discourages farmers from investing more inputs resulting in lower productivity, even in years when rainfall is favorable for a higher yield [12]. Thus, farmers in this ecosystem are generally suffering from poverty and malnutrition as compared to those in irrigated ecosystems [11,12]. The absolute poverty rate is more severe in rainfed areas of sub-Saharan Africa, but more farmers are also affected in rainfed areas of Asia [11]. Additionally, even productivity and economic progress in irrigated rice areas are at high risk due to the increasing frequency of water shortages in the world. The decreasing water reserved in dams affects rice areas at the tail end of the irrigation system. Consequently, some farmers may leave rice farming [13] and convert their rice fields for another enterprise, thereby reducing rice cultivated areas and production volumes. Therefore, considering the extent of the negative impact of drought, increasing rice productivity and/or profitability in water-scarce areas is a serious concern.
Similar to other crop research and development programs, increasing the yield has been the major goal of countless explorations into rainfed lowland ecosystems. The development of a high-yielding and drought-tolerant rice variety is one of the keys to achieving this goal. However, combining high-yielding and drought-tolerant traits remains a major challenge due to the complexity of both traits [14,15]. The plant selection through secondary traits, involving either shoots or roots, cannot effectively support the expected progress in breeding for drought tolerance. Other plant selection, such as the direct selection for grain yield instead of using secondary traits has been successful [16,17], if the release of several rainfed rice varieties in many countries is the gauge [18]. However, despite these advancements, the reported genetic gain from the drought-tolerant rice breeding program by the International Rice Research Institute (IRRI) for 17 years (2003–2019) was only 0.13–0.55%, which is equivalent to 2.29–9.52 kg ha−1 yr−1, and below the 1.5% target to meet future global demand [19]. The insufficiency of rice supply in this resource-poor ecosystem may lead to widespread malnutrition, hunger, and potential social unrest.
Considering a low-yielding environment like a rainfed area, planting of a rice variety that commands a higher market price could potentially improve farming income. Specialty rice, such as aromatic, pigmented, and glutinous rice, represents unique types of rice with a high market price due to their exceptional grain quality, such as their aroma, kernel color, and chemical composition [20,21,22]. Aromatic rice is known for its nutty aroma and taste caused by the main aroma compound 2-acetyl-1-pyroline or 2AP [23], with the Basmati type from India and the Jasmine type from Thailand being prime examples. Pigmented rice, on the other hand, has the highest nutritional value because it contains pigments (anthocyanins), which are potent antioxidants [22,24]. There is increasing interest in this type of rice because of its significant health and nutritional benefits. Pigmented rice is categorized as black and red rice. Furthermore, glutinous rice has sticky characteristics because of its high amylopectin content, thus making it suitable as a raw material for rice-based products and rice delicacies, such as butchi in China, mochi in Japan, and kakanin in the Philippines [21]. The superior features and quality of specialty rice make them distinguishable from ordinary white rice, thus commanding a better price on the market. In Thailand, planting specialty rice in rainfed environments, specifically, the aromatic type, generated a higher income even when compared to that of favorable, irrigated conditions planted with ordinary rice [25]. The better market price for aromatic rice can compensate for the low yield, which made Thailand cultivate its large rainfed areas using the Jasmine type of aromatic rice. This aromatic rice is also popular on the international market [26], implying that a rainfed ecosystem can be utilized to produce high-quality rice. In the Philippines, some farmers also prefer to plant rice varieties with higher market prices not only in low-yielding environments (rainfed) but also during low-yielding seasons, such as the wet season in irrigated areas. For instance, some farmers prefer popular varieties (NSIC Rc160 and Rc216) with intermediate amylose content (AC) because of their higher market price rather than high-yielding varieties with high AC (NSIC Rc222) [27]. In the case of specialty rice, however, local production under a rainfed lowland system has not yet been reported on. In recent years, the development and release of modern specialty rice varieties was targeted in irrigated lowland (aromatic and pigmented types) and upland conditions (glutinous type), but not in rainfed lowland conditions [28]. Furthermore, the grain quality of the majority of the released rainfed rice (inbred) varieties in the Philippines was inferior to that of specialty rice [28]. In the present study, therefore, we investigated specialty rice varieties for yield, stability, and drought tolerance under rainfed lowland cultivation conditions, and identified superior varieties better suited to these conditions.

2. Materials and Methods

2.1. Plant Materials

A total of 17 SR genotypes pre-selected based on the reported yield and popularity [28,29], together with the drought-tolerant (DT) check variety, NSIC Rc192, were used in this study (Table 1). The SR genotypes were composed of 6 aromatic (NSIC Rc34, Rc218, Rc342, Rc344, Basmati 370, and CLRice-1), 6 pigmented (Black rice, Red rice, Calatrava, Pinilisa, CLRice-2 and 3), and 5 glutinous (NSIC Rc15, 17, 19, 21, and 31) genotypes. Those with NSIC codes are released varieties intended for irrigated lowland (aromatic) and upland cultivation (glutinous), while the pigmented rice except for those with CLRice codes are traditional varieties. Rice genotypes with CLRice codes are advanced breeding lines developed by the Central Luzon State University (CLSU), Muñoz, Nueva Ecija (15°44′ N, 120° 56′ E). Additionally, Black rice and Red rice were sourced from Sta. Rosa, Nueva Ecija (15°25′ N, 120°56′ E), while Calatrava was collected from Calatrava, Negros Occidental (10°36′ N, 123°29′ E). Both types of rice were cultivated in irrigated lowland conditions. Pinilisa, on the other hand, was collected from Jones, Isabela (16°33′ N, 121°42′ E), and cultivated in upland conditions.

2.2. Field Evaluation

To analyze the responses of the genotypes under drought-prone rainfed lowland conditions, the above 18 genotypes were evaluated in rainfed and controlled drought conditions for two years, 2019 wet season (WS)—2021 dry season (DS) (June 2019–April 2021), in the experimental area of the Department of Crop Science, the College of Agriculture, CLSU, Muñoz, Nueva Ecija, the Philippines (15°44′ N, 120°56′ E, 80 masl). In the WS experiments, normal and late planting in rainfed conditions were conducted. Late planting was included in the rainfed conditions to expose the plants to terminal drought stress. To ensure successful crop establishment, irrigation was supplemented during land preparation until the imposition of the drought treatment at 14 days after transplanting (DAT). Thereafter, the source of irrigation up to maturity was solely from rainfall.
In the 2019 WS, the normal schedule of planting involved the rice being sown on June 25, while the late planting was one month later (July 25). In the 2020 WS, normal planting involved the rice being sown on June 26, while the late planting was 15 days later (July 11) or 15 days earlier than the late planting for the previous year. The changes in the late planting schedule were conducted due to the early cessation of rainfall in the previous year. In the DS experiments, the genotypes were further evaluated under intermittent-mild and severe-drought conditions in the 2020 DS, while only under mild-drought conditions during the 2021 DS. Likewise, irrigated conditions served as the control. Similar to the WS experiments, the fields were irrigated to facilitate land preparation and for faster crop recovery from transplanting shock. Thereafter, the surface water from the field was drained at 14 DAT. Re-irrigation was conducted when the susceptible check reached a leaf rolling (LR) score of 3 (deep V) under intermittent-mild drought and an LR score of 7 (O shape) for severe drought in the morning. The seeds were sown on 26 December 2019, for the 2020 DS, while it was 22 December 2020, for the 2021 DS controlled-drought screening. Soil moisture content (SMC, %) and water table depth (WTD, cm) were recorded weekly, using the gravimetric method and using an installed piezometer with a 1 m depth, respectively. Daily rainfall was recorded from the weather station at the Philippine Atmospheric, Geophysical and Astronomical Services Administration (PAGASA), which was located 350 m northwest from the experimental area.

2.3. Experimental Design and Management

The seeds were sown in a raised seedbed, sprinkled with carbonized rice hull, then covered with sacks to protect the seeds from rain drops, birds, and rodents. The sack covers were removed at 5–7 days after sowing (DAS) or when seedlings had a shoot length of approximately 5 cm. Complete fertilizer (14-14-14) at a rate of 35.71 g m−2 was applied at 14 DAS. The seedlings were transplanted at 21 DAS at a distance of 20 cm × 20 cm within a 5.0 m × 4.0 m plot. The susceptible check was planted at both ends and at the middle of each block within a 5.0 m × 0.6 m plot. The entries were arranged in an RCBD, with three replications for each water treatment. A fertilizer rate of 90-30-30 kg NPK ha−1 was applied in both seasons. So, each plot received 428.57 g 14-14-14 at 7 DAT, and 130.43 g 46-0-0 both at 21–28 DAT and at 35–42 DAT.

2.4. Analysis of Drought Tolerance Indices

The grain yield data at 14% MC were used to compute different drought tolerance indices such as geometric mean productivity index (GMP), stress tolerance index (STI) [30], mean productivity index (MP) [31], harmonic mean index (HM) [32], tolerance index (TOL) [31], and stress susceptibility index (SSI) [33]. The formula used to compute these indices are shown below:
GMP   = Ypi   ×   Ysi
MP   = Ypi + Ysi 2
HM   = 2   Ysi   ×   Ypi Ysi + Ypi
STI   = Ysi   ×   Ypi Yp 2
TOL   = Ypi Ysi
SSI   = 1 Ysi Ypi SI
where SI = 1 − (YsYp), Ypi is the yield under non-stress for the genotype, Ysi is the yield under stress for the genotype, Yp is the mean grain yield under non-stress, and Ys is the mean grain yield under stress.

2.5. Analysis of Milling and Grain Quality

The grains were immediately sundried between 1 to 5 d after harvest, depending on the grain moisture content and weather conditions. Upon reaching 14% MC, the grains were milled using 3 HP portable electric rice milling machines (SATO, Sendai, Japan). Milling recovery (MR, %) was computed using the formula:
Milling   recovery = weight   of   milled   rice weight   of   grain   samples × 100  
The milled rice with at least 75% of whole grains, considered as the head rice (HR), was manually separated. The HR recovery was computed using the formula:
Head   rice   recovery = weight   of   head   rice weight   of   grain   samples × 100  

2.6. Statistical Analysis

To compare the performance between the water treatments within season, the data were subjected to a combined ANOVA. The analyses were performed using Statistical Tools for Agricultural Research (STAR) software version 2.0.1 developed by IRRI. Treatment means were compared using the least significant difference (LSD) at the 5% level of significance.
To evaluate the adaptability, productivity, and stability across the water treatments, the which-won-where GGE biplot was used [34], using the ‘metan’ package [35], and the weighted average of absolute scores (WAASB) for the best linear unbiased predictions (BLUPs) [36] were performed using the packages in RStudio, R v. 4.3.0 [37].

3. Results and Discussion

3.1. SMC, WTD, and Agro-Climatic Condition

The soil moisture content (SMC), water table depth (WTD), and daily rainfall in the rainfed and controlled-drought conditions from 14 DAT up to maturity and other agro-climatic parameters are presented in Figure 1. The rainfed field condition during the 2019 WS received less than half of the cumulative rainfall (1126 mm), but it was evenly distributed throughout the growing period compared to 2020 WS (2823 mm) (Figure 1A,C,E). In the 2019 WS normal schedule of planting, most of the genotypes were exposed to 23–46% SMC or an average of 39% SMC from the vegetative to early ripening phase. A low SMC of 14% and a WTD of 59 cm below the soil surface were recorded only at 126 DAS. In the late schedule of planting, on the other hand, the evaluation of terminal drought stress was affected by the early onset of the dry season, thus a low amount of rainfall (4 mm d−1) was received starting from early October, resulting in a low SMC of 14%. Concurrently, increasing solar radiation and temperature and a decreasing RH were also evident (Supplemental Table S1). The low rainfall and SMC affected the test entries from the reproductive to the ripening phase, resulting in the failure in the seed setting of most of the test entries, including the drought-tolerant (DT) check variety. Thus, only the analysis on the normal schedule of planting is available for the 2019 WS results. In the 2020 WS, late planting was adjusted 15 d earlier than the previous year. Both the late and normal planting in the rainfed conditions resulted in reductions in the grain yields relative to their irrigated counterparts. This reduction could be due to the uneven distribution of rainfall despite the higher precipitation (above 2800 mm) compared to the previous year (Figure 1). Around 72% of the total rainfall was received during only three periods (37–39 DAT, 65–67 DAT, and 113 DAT in normal planting or 26–28 DAT, 54–56 DAT, and 102 DAT in late planting). Consequently, normal and late planting had an average of 32% and 29% SMC, respectively, which were lower than the observed 35% in the 2019 WS. Moreover, between normal and late planting in the 2020 WS, the latter had a lower average SMC (26%) from flowering to ripening than the former (32% SMC). On the other hand, despite the unexpected occurrence of rainfall during the 2020 DS, a low SMC and WTD (24% and 48 cm for mild and 20% and 56 cm for severe) were recorded in the controlled-drought conditions. In the 2021 DS, the agro-climatic condition became drier with almost no rainfall occurrence and a lower RH. The recorded average SMC in mild-drought conditions was 23%, with the WTD at 55 cm below the soil surface. These SMC values are close to the reported SMC values in rice drought-tolerance screening during the WS and the DS in the Philippines [3,38].

3.2. Grain Yield

The genotype, water treatment, and genotype-by-water treatment interactions on the yield response were significant over the seasons and years (except for the interaction in the 2019 WS) (Supplemental Table S2). The consistently higher phenotypic coefficient of variation (PCV) than the genotypic coefficient of variation (GCV) in all water treatments (Supplemental Table S3) suggests that the environment influenced the observed variation in the grain yield. The disparity between the PCV and GCV became more pronounced under water-limited conditions, with at least a 50% yield reduction in all controlled-drought conditions and rainfed, except for the 2019 WS. This indicates increasing environmental influence on the grain yield. Consequently, the heritability values were reduced to 36–54%, but still classified as moderate heritability [39]. A moderate heritability of grain yield under the drought stress condition was the basis for the proposed direct selection for grain yield [16], which has proven effective in developing several rainfed rice varieties in different countries [18].
The present study also focused on grain yield and employed a proven method in selecting superior genotypes. The grain yield ranged from 0.86 (Pinilisa) to 6.90 t ha1 (NSIC Rc192) under irrigated conditions, and 1.45 (Calatrava) to 6.78 t ha1 (NSIC Rc342) under the rainfed conditions in the 2019–2020 WS (Table 2). In the DS, the grain yield ranged from 3.52 (Pinilisa) to 9.50 t ha1 (NSIC Rc15) under irrigated conditions, whereas the grain yield was 0.87 (Pinilisa) to 4.21 t ha1 (NSIC Rc192) and 0.00 (Pinilisa) to 2.26 t ha1 (NSIC Rc342) under mild and severe drought conditions, respectively (Table 2). The DT check variety, NSIC Rc192, obtained a grain yield of 1.78–6.33 t ha1 under rainfed conditions and 0.98–4.21 t ha1 under controlled-drought conditions with an average of 4.32 and 2.18 t ha1, respectively. The average grain yield in the rainfed conditions was relatively higher than those observed in other studies, 3.63 t ha1 in other regions in the Philippines [40,41] and 3.50 t ha1 in Nepal [42,43]. The grain yields for the DT check variety were reduced by 8–67% and 47–88% under the rainfed and controlled-drought conditions, respectively, relative to the irrigated control. Despite these reductions, it had a significantly higher grain yield than most of the SR. However, there were some genotypes from each SR group with a comparable or higher yield than the check variety under rainfed and controlled-drought conditions. Under rainfed conditions in the 2019 WS, two aromatic (NSIC Rc344 and Rc342), two pigmented (Black rice and CLRice-2), and one glutinous (NSIC Rc15) variety obtained a comparable grain yield to the DT check variety, NSIC Rc192 (6.33 t ha1). NSIC Rc344 (6.39 t ha1) had the highest grain yield among the aromatic rice, which was comparable to NSIC Rc342 (5.64 t ha1). Among the pigmented rice, Black rice produced the highest grain yield (5.98 t ha1) and was also comparable to CLRice-2 (5.91 t ha1). In the glutinous group, only NSIC Rc15 obtained the highest grain yield (6.85 t ha1). In contrast to the 2019 WS, the following year (the 2020 WS), specifically during the normal schedule of planting in the rainfed conditions, only one SR, the aromatic NSIC Rc342, obtained a significantly higher grain yield (6.78 t ha1) than the DT check variety (4.84 t ha1), whereas the rest of the aromatic rice together with the three pigmented (Black rice, CLRice-2, and Red rice) varieties and all the glutinous rice obtained a comparable grain yield (3.88–5.68 t ha1). During the late planting, which was exposed to terminal drought stress with an average yield reduction of 50% in all entries, three aromatic (NSIC Rc344 with 3.66 t ha1, Basmati 370 with 3.18 t ha1, and CLRice-1 with 2.60 t ha1) and one glutinous rice (NSIC Rc15 with 3.30 t ha1) produced a higher grain yield than the DT check variety, while the rest of the SR had a comparable yield. The DT check variety obtained a grain yield of only 1.78 t ha1, despite its shorter growth duration (111 DAS) compared to the majority of the SR (115–125 DAS). Overall, in both years, two aromatic (NSIC Rc344 and Rc342), two pigmented (Black rice and CLRice-2), and one glutinous (NSIC Rc15) rice had a consistently comparable if not higher grain yield than the DT check variety under rainfed conditions. Although conducted in different locations, the recorded grain yields of 1.80–6.78 t ha1 from the two high-yielding aromatic rice (NSIC Rc344 and Rc342) were almost similar to the reported grain yields of 1.40–6.29 t ha1 from the progenies of the aromatic and drought-tolerant rice genotypes [44,45]. Moreover, the maximum grain yields that were observed in the present study were close to the reported maximum yield under rainfed lowland conditions (4.34–5.40 t ha1) [3,38,46,47], even for hybrid rice varieties (5.8 t ha1) [48]. This reconfirms the persistently low yield in rainfed lowland conditions, underscoring the significance of the high market price for SR to improve farmers’ income.
Controlled-drought conditions during the dry season were conducted to validate the performance of the SR in rainfed field conditions. Among the identified genotypes under the rainfed conditions, only the aromatic NSIC Rc344 rice (3.34 t ha1) and pigmented Black rice (3.19 t ha1) obtained a comparable grain yield to the DT check variety (4.21 t ha1) in the mild drought during the 2020 DS (Table 2). These SR have numerically highest grain yields in their respective groups, whereas the NSIC Rc15 rice (2.04 t ha1) had the highest grain yield in the glutinous group. When drought stress became severe, the aromatic NSIC Rc342 rice (2.26 t ha1) and Black rice (2.11 t ha1) had higher grain yields than the DT check variety (0.98 t ha1) by 115–131%, while the rest of the SR had a comparable yield (0–1.81 t ha1) to the DT check variety (0.98 t ha1). In the following year, the 2021 DS, almost all the SR (0.87–2.51 t ha1) produced a comparable grain yield to the DT check variety (1.34 t ha1), except a higher grain yield was recorded for the aromatic NSIC Rc344 rice (2.60 t ha1). Collectively, the aromatic NSIC Rc344 rice and Black rice had at least a comparable grain yield to the DT check variety under rainfed and controlled-drought conditions, while the glutinous NSIC Rc15 rice together with the aromatic NSIC Rc342 rice and pigmented CLRice-2 rice also had at least a comparable grain yield to the DT check variety, except under mild drought stress in the 2020 DS.

3.3. Evaluation of Yield Performance Using GGE Biplot Analyses

Since the SR were evaluated for different water treatments and years, a discriminativeness and representativeness analysis was conducted using GGE biplot analysis. This is to determine the discrimination capability and representativeness of the water treatments [34], particularly within and between the rainfed and controlled-drought conditions. In these important water treatments, the largest angle vector was an acute angle with a closer angle or distance within the rainfed conditions and between the normal planting schedule in the rainfed and controlled-drought conditions (Figure 2A). An acute angle between the two environments indicates a positive correlation or the two environments provide the same information about the genotypes [34]; thus, the results for the rainfed conditions were validated or supported by the controlled-drought conditions. Interestingly, the two most important water treatments, namely normal planting in rainfed conditions for both years, were the closest to the average environment axis and were, thus, the best representative of all the water treatments. More importantly, the 2019 normal planting in the rainfed conditions was one of the water treatments with the longest vector or strongest discriminating power, suggesting its effectiveness in differentiating specialty rice in terms of the yield levels.
The grain yield of the SR were further analyzed using GGE biplot analyses, specifically the which-won-where was used to select the superior SR suitable for rainfed or drought stress conditions. The which-won-where biplot analysis was performed to simultaneously evaluate the grouping or similarity between the different water treatments and the corresponding adaptability of the SR genotypes. The analysis showed that the different water treatments could be grouped into three distinct mega environments (Figure 2B). The two normal plantings in the rainfed conditions during the 2019 WS and the 2020 WS, together with the irrigated conditions in the 2020 DS and mild drought stress in the 2021 DS, belonged to one mega environment, indicating the similarity of these water treatments in ranking the yield performance of the SR [34]. Three out of the five SR identified above, namely the aromatic NSIC Rc342 rice and pigmented Black rice and CLRice-2, were suitable in this group of water treatments, in which the first genotype had the highest grain yield. The other rainfed conditions, namely the 2020 WS late planting, were included in another group, together with the irrigated conditions in the 2020 WS and the 2021 DS. The suitable genotypes included here were the glutinous NSIC Rc15 rice and Rc21, in which the former was the winning genotype and was among the five SR identified earlier as high yielding. On the other hand, the two levels of drought stress (mild and severe) in the 2020 DS, together with the irrigated conditions in the 2019 WS, were clustered into one mega environment. The adapted genotypes included here were the DT check variety and three aromatic rice, namely NSIC Rc344, Rc218, and Rc34, whereby the DT check variety was the best-performing genotype.
Furthermore, biplot analysis using the weighted average of absolute scores (WAASB) for the best linear unbiased predictions (BLUPs) was also conducted to evaluate productivity and stability over a wide range of hydrological conditions [36]. This method is known for identifying genotypes with both high and stable yield across diverse growing conditions [49,50,51]. The biplot is composed of four quadrants or groupings. Quadrant I groups genotypes that have both low productivity and stability, quadrant II groups genotypes with high productivity but low stability, quadrant III groups genotypes with low productivity but high stability, and quadrant IV groups genotypes with both high productivity and stability, which is the most desirable. Genotypes clustered in quadrant IV include the pigmented Black rice and CLRice-2, the aromatic NSIC Rc344, Rc218, Rc34, and Basmati 370 rice, and the glutinous NSIC Rc15 and Rc21 rice, which have a higher productivity and stability relative to the other genotypes (Figure 2C). Within the SR group in this quadrant, Black rice had a higher and more stable grain yield than CLRice-2; NSIC Rc34 had the most stable grain yield among the aromatic group, but NSIC Rc344 had the highest productivity; and glutinous NSIC Rc21 rice was more stable than NSIC Rc15, although the latter had a higher grain yield. Interestingly, the DT check variety, as well as the aromatic NSIC Rc342 rice, were found in quadrat II indicating a high but unstable grain yield across the water treatments. In summary, the biplot analyses confirmed the superior performance of the two aromatic (NSIC Rc342 and Rc344), two pigmented (Black rice and CLRice-2), and one glutinous (NSIC Rc15) rice across the water treatments.

3.4. Drought Tolerance Indices

To further select the best genotype from each group of SR, especially the identification of drought-tolerant ones, the evaluation of drought tolerance indices was employed using the grain yield data. The indices used included GMP, MP, HMI, and STI to identify the genotype/s with a high grain yield under the irrigated conditions and a reasonably high yield in both the rainfed and drought stress conditions [52]. These indices were found to be effective in identifying genotypes with a high yield in both favorable and stress conditions in numerous studies [52,53,54,55,56] and were also found to be strongly associated with yield under favorable, as well as rainfed and drought stress conditions in this study (r ≥ 0.96 **). The results showed that four (aromatic NSIC Rc344 and Rc342, pigmented Black rice, and glutinous NSIC Rc15) out of the five identified SR and the DT check variety were the top five genotypes in the GMP, MP, HMI, and STI under the rainfed conditions, with the four SRs achieving better rankings than the DT check variety (Table 3). This suggests that these four SRs have a higher grain yield in rainfed and irrigated conditions during the wet season than the rest of the SRs evaluated and the DT check variety. Other drought tolerance indices, TOL and SSI, were used to identify the genotypes with a stable grain yield across the hydrological conditions [31,33]. Among the tested genotypes, only the glutinous NSIC Rc15 rice and four other genotypes, namely the aromatic CLRice-1, pigmented CLRice-3, glutinous NSIC Rc21, and Rc17 rice, were in the top five for TOL and SSI, suggesting that the glutinous NSIC Rc15 rice had a high and stable grain yield across the water treatments during the wet season. Furthermore, the other identified SRs (aromatic NSIC Rc344 and Rc342, and pigmented Black rice) had better TOL and SSI rankings than the DT check variety, suggesting that they had a more stable yield than the DT check variety. On the other hand, the top five genotypes for the GMP, MP, HMI, and STI under the rainfed conditions were almost the same as under the controlled-drought conditions, except for the glutinous NSIC Rc15 rice, which was displaced by the aromatic NSIC Rc34 rice. Among the glutinous rice, however, the NSIC Rc15 rice still had the most desirable drought tolerance indices under the controlled-drought conditions. Moreover, the DT check variety, the aromatic NSIC Rc344 and Rc342 rice, and the Black rice were in the top five for the SSI under the controlled-drought conditions, with the SR showing better rankings or a more stable grain yield than the DT check variety. Between the aromatic NSIC Rc344 and Rc342 rice, the former had better TOL and SSI rankings than the latter under the rainfed and controlled-drought conditions indicating a more stable grain yield, which supports the results of the productivity and stability test involving the WAASB biplot. Taken together, the drought tolerance indices support the results of the biplot analysis, specifically the WAASB test in which the aromatic NSIC Rc344 rice, pigmented Black rice, and glutinous NSIC Rc15 rice have a higher and more stable grain yield than the DT check variety. Aside from the high yield of the identified SR, stability is advantageous from the farmer’s perspective [57,58], especially in high-risk production areas in the face of a changing climate [59], like rainfed rice areas [60].

3.5. Unique Yield Characteristics of Selected Specialty Rice

The superior yield performance of the identified SR was further examined by analyzing the yield-contributing traits across the water treatments using the ranking of genotypes in biplot analysis. The analysis showed that among the SR, the aromatic NSIC Rc344 rice had the highest panicle number, the glutinous NSIC Rc15 rice had the highest number of filled grains per panicle, while the Black rice had the heaviest 1000-grain weight (except when compared to NSIC Rc34) and the highest harvest index across the water treatments (Figure 3). This result supports the high grain yield of the selected genotypes across the water treatments. The yield-contributing traits of the genotypes differ from each other. This likely genetic complexity behind the regulation of yield trait expression is a necessary consideration for effective selection and crop improvement efforts for SR. Furthermore, such yield traits have been reported to be related with other desirable traits for drought resistance. For instance, the yield attributed by the high number of filled grains per panicle was closely associated with the production of more nodal and lateral roots [61], which are important for water uptake under drought or fluctuating soil moisture conditions [62,63]. Such trait associations or complementarities may result in a better understanding of the improved performance of SR. Therefore, further studies are also required to uncover the mechanisms of drought resistance involving these genetic materials, possibly by avoidance, tolerance, and other known mechanisms, studies of which are currently ongoing.

3.6. Head Rice Recovery

To determine the profitability of planting SR with a higher market price than the ordinary-grain quality DT check variety under rainfed conditions, the head rice recovery (HRR) was evaluated, since the price of SR is more established in the milled form on the local and international markets [21]. Only the head rice of those genotypes planted in normal planting schedules in both years are presented because of the insufficient amount of grain samples for each variety from the late planting schedule. The ANOVA showed that the year and genotype x year interaction was not significant on the HRR (Table 4), indicating that this trait was stable. The HRR was influenced only by the genotype. The HRR of the released varieties was higher than those of improved lines and traditional varieties. This is because the HRR is one of the important criteria in the approval and release of new rice varieties, of which the recommended value is at least 48% [64]. Most of the HRR of the newly released rice varieties were higher than the recommended value. The aromatic NSIC Rc344 rice and the glutinous NSIC Rc15 rice are released varieties with a HRR of 53% and 56%, respectively. On the other hand, Black rice, a traditional variety, had only a 41% HRR, which was statistically comparable to those HRR of other pigmented rice, such as Red rice (48%) and Calatrava (47%). However, it is important to note that this value for Black rice is underestimated because all the test entries were milled using only one type of equipment designed for polished rice. Nevertheless, the milled rice from pigmented entries still had significant level of antioxidant activity (2.05–5.70% w/w), with Black rice containing one of the highest values (Supplemental Table S5).

3.7. Profitability Analysis

The HRR value of the selected SR, as well as the DT check variety, was multiplied together with its grain yield (ha1) to estimate the gross income (Figure 4A,B). The lowest-yielding genotype from each group was also included to provide the range of income from planting SR under rainfed conditions. The market prices for SR in milled form were based on the reported prevailing price (kg1), namely USD 1.46 (PHP 72.00) for aromatic rice, USD 1.56 (PHP 77.00) for pigmented rice, and USD 1.20 (PHP 59.00) for glutinous rice [21], while USD 0.91 (PHP 45.00) for rice with ordinary grain quality was used [65,66]. The net income was estimated as the difference between the gross income and the total production cost from crop establishment to milling (Supplemental Tables S6–S9). The net income (ha1) was as follows: USD 760.42 to 2748.33 for aromatic rice, USD 208.17–2229.59 for pigmented rice, and USD 1336.69–2410.70 for glutinous rice, while it was USD 1320.42 for the DT check variety (Figure 4A). Accordingly, these gave a comparative net income of −42–108% for aromatic rice, −83–69% for pigmented rice, and 1–83% for glutinous rice (Figure 4B). Despite the low HRR of Black rice, it still provided a >69% net income than the DT check variety, attributed to the higher yield and higher market price (USD 1.56 kg1 vs. USD 0.91 kg1). Also, given that the other two selected genotypes had a comparable/higher yield and HRR when compared to the DT check variety, the higher market prices of USD 1.46 kg−1 for NSIC Rc344 rice and USD 1.20 kg−1 for NSIC Rc15 rice increased the net income by 108% and 83%, respectively. SR also has a higher market price than ordinary rice in other countries, such as Thailand, India, Indonesia, Laos, Cambodia, and the international market [20,21,22]; thus the use of SR may improve the profitability of global rainfed lowland rice environments. The improvement in net income means more investment in farming that may help address other problems to improve the productivity of low yielding rainfed lowland environments in the long term. Therefore, contributing to breaking the cycle of low yield, low income, and low input use, which may help alleviate poverty in rainfed agriculture.

4. Conclusions

In a low-yielding environment, such as in rainfed lowland environments, the planting of specialty rice with superior grain quality and nutritional values could contribute to the feasibility of rice farming amidst the changing climate. In this study, we showed the potential of specialty rice cultivation in rainfed lowland conditions, which has not yet been reported on by any other studies. We identified the superior variety in each group of specialty rice, with a comparable or higher grain yield than the drought-tolerant check variety (NSIC Rc192) under rainfed or water-limited conditions, and we found three specialty rice genotypes, namely aromatic NSIC Rc344 rice, pigmented Black rice, and glutinous NSIC Rc15 rice, under rainfed field conditions. This selection was further supported by GGE biplot analyses and drought tolerance indices, whereby the selected specialty rice were classified as the highest yielding among the specialty rice with a more stable grain yield than the drought-tolerant check variety across a wide range of hydrological conditions. The high and stable grain yield of the selected specialty rice was contributed by the high number of panicles for the NSIC Rc344 rice, the high 1000-grain weight and harvest index for the Black rice, and the high number of filled grains per panicle for the NSIC Rc15 rice. Given the comparable to higher grain yield of the identified specialty rice, but with a higher market price than the drought-tolerant check variety, the planting of the selected specialty rice under rainfed conditions increased the net income by 69–108% as compared to the planting of the drought-tolerant check variety with ordinary grain quality.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/agriculture13101985/s1. Table S1: Agroclimatic condition of the area. Table S2: Combined ANOVA for grain yield of all entries under different water treatments in wet and dry seasons from 2019 to 2021. Table S3. Estimation of some genetic parameters for grain yield in different water treatments. Table S4: Hydrological parameters, mean grain yield and grain yield reduction of specialty rice and drought-tolerant check under rainfed and controlled drought conditions. Table S5: Antioxidant activity (%w/w) of pigmented rice under rainfed condition. Table S6: Cost and return per hectare of drought tolerant check NSIC Rc192 under rainfed condition of 2019 WS. Table S7: Cost and return per hectare of aromatic NSIC Rc344 under rainfed condition of 2019 WS. Table S8: Cost and return per hectare of drought tolerant check NSIC Rc192 under rainfed condition of 2020 WS. Table S9: Cost and return per hectare of aromatic NSIC Rc344 under rainfed condition of 2020 WS.

Author Contributions

Conceptualization, A.M.L.A. and J.L.O.; methodology, A.M.L.A., J.L.O. and M.B.S.N.; formal analysis, A.M.L.A. and M.K.-N.; writing—original draft preparation, A.M.L.A.; writing—review and editing, N.B.L.-A., R.R.S., M.K.-N., H.E. and S.M.; funding acquisition, A.M.L.A., J.L.O. and M.K.-N. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Bureau of Agricultural Research at the Department of Agriculture of the Philippines (DA-BAR), grant number RICE-2019-14. The APC was funded by Nagoya University, Asian Satellite Campus Institute (the Philippines).

Institutional Review Board Statement

Not applicable.

Data Availability Statement

Data are accessible from the main author on reasonable request.

Acknowledgments

The authors extend their appreciation to the DA-BAR of the Republic of the Philippines for providing funds for the study. We are also grateful to Via Ann C. Marcelo for her help in the GGE biplot analysis.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. GRiSP (Global Rice Science Partnership). Rice Almanac, 4th ed.; International Rice Research Institute: Los Baños, CA, USA, 2013. [Google Scholar]
  2. Lafitte, H.R.; Yongsheng, G.; Yan, S.; Li, Z.-K. Whole plant responses, key processes, and adaptation to drought stress: The case of rice. J. Exp. Bot. 2007, 58, 169–175. [Google Scholar] [CrossRef] [PubMed]
  3. Torres, R.O.; McNally, K.L.; Vera Cruz, C.; Serraj, R.; Henry, A. Screening of rice Gene bank germplasm for yield and selection of new drought tolerance donors. Field Crops Res. 2013, 147, 12–22. [Google Scholar] [CrossRef]
  4. Quiñones, C.; Mattes, N.; Faronilo, J.; Yadav, S.; Jagadish, K.S.V. Drought stress reduces grain yield by altering floral meristem development and sink size under dry-seeded rice cultivation. Crop Sci. 2017, 57, 2098–2108. [Google Scholar] [CrossRef]
  5. Zhang, Y.; Zhao, X.; Liu, F.; Zhu, L.; Yu, H. Effect of Different Water Stress on Growth Index and Yield of Semi-Late Rice. Environ. Sci. Proc. 2023, 25, 84. [Google Scholar] [CrossRef]
  6. Praba, M.L.; Cairns, J.E.; Babu, R.C.; Lafitte, H.R. Identification of Physiological Traits Underlying Cultivar Differences in Drought Tolerance in Rice and Wheat. J. Agron. Crop Sci. 2009, 195, 30–46. [Google Scholar] [CrossRef]
  7. Salgotra, R.K.; Chauhan, B.S. Ecophysiological Responses of Rice (Oryza sativa L.) to Drought and High Temperature. Agronomy 2023, 13, 1877. [Google Scholar] [CrossRef]
  8. Lawlor, D.W. Limitation to Photosynthesis in Water-stressed Leaves: Stomata vs. Metabolism and the Role of ATP. Ann. Bot. 2002, 89, 871–885. [Google Scholar] [CrossRef]
  9. Wang, X.; Huang, J.; Peng, S.; Xiong, D. Leaf rolling precedes stomatal closure in rice (Oryza sativa) under drought conditions. J. Exp. Bot. 2023, erad316. [Google Scholar] [CrossRef]
  10. Mackill, D.J.; Coffman, W.R.; Garrity, D.P. Rainfed Lowland Rice Improvement; International Rice Research Institute: Manila, Philippines, 1996; 242p. [Google Scholar]
  11. Department of Agriculture-Bureau of Agricultural Research. Philippine Rainfed Agriculture Research; Development and Extension Program (PhiRARDEP) Framework and Action Agenda: Quezon City, Philippines, 2011. [Google Scholar]
  12. Pandey, S.; Bhandari, H. Drought: An overview. In Economic Costs of Drought and Rice Farmers’ Coping Mechanisms: A Cross-Country Comparative Analysis; Pandey, S., Bhandari, H., Hardy, B., Eds.; International Rice Research Institute: Los Baños, Philippines, 2007; 203p. [Google Scholar]
  13. Manalo, J.A., IV; van de Fliert, E.; Fielding, K. Rice farmers adapting to drought in the Philippines. Int. J. Agric. Sustain. 2020, 18, 594–605. [Google Scholar] [CrossRef]
  14. Li, G.; Tang, J.; Zheng, J.; Chu, C. Exploration of rice yield potential: Decoding agronomic and physiological traits. Crop J. 2021, 9, 577–589. [Google Scholar] [CrossRef]
  15. Panda, D.; Mishra, S.S.; Behera, P.K. Drought tolerance in rice: Focus on recent mechanisms and approaches. Rice Sci. 2021, 28, 119–132. [Google Scholar] [CrossRef]
  16. Venuprasad, R.; Lafitte, H.R.; Atlin, G.N. Response to direct selection for grain yield under drought stress in rice. Crop Sci. 2007, 47, 285–293. [Google Scholar] [CrossRef]
  17. Kumar, A.; Bernier, J.; Verulkar, S.; Lafitte, H.R.; Atlin, G.N. Breeding for drought tolerance: Direct selection for yield, response to selection and use of drought-tolerant donors in upland and lowland-adapted populations. Field Crops Res. 2008, 107, 221–231. [Google Scholar] [CrossRef]
  18. Kumar, A.; Dixit, S.; Ram, T.; Yadaw, R.B.; Mishra, K.K.; Mandal, N.P. Breeding high-yielding drought-tolerant rice: Genetic variations and conventional and molecular approaches. J. Exp. Bot. 2014, 65, 6265–6278. [Google Scholar] [CrossRef] [PubMed]
  19. Khanna, A.; Anumalla, M.; Catolos, M.; Bartholomé, J.; Fritsche-Neto, R.; Platten, J.D.; Pisano, D.J.; Gulles, A.; Sta Cruz, M.T.; Ramos, J.; et al. Genetic Trends Estimation in IRRIs Rice Drought Breeding Program and Identification of High Yielding Drought-Tolerant Lines. Rice 2022, 15, 14. [Google Scholar] [CrossRef] [PubMed]
  20. Juliano, B.O.; Villareal, C.P. Grain Quality Evaluation of World Rices; International Rice Research Institution: Los Baños, Philippines, 1993. [Google Scholar]
  21. Beltran, J.C.; Bordey, F.H.; Moya, P.F.; Relado-Sevilla, R.Z.; Romero, M.V. Philippine Specialty Rice: Understanding Production, Culture, Quality, and Market; Philippine Rice Research Institute and Manila (Philippines); International Rice Research Institute: Science City of Muñoz, Philippines, 2020; 380p. [Google Scholar]
  22. Yamuangmorn, S.; Prom-u-Thai, C. The potential of high-anthocyanin purple rice as a functional ingredient in human health. Antioxidants 2021, 10, 833. [Google Scholar] [CrossRef] [PubMed]
  23. Buttery, R.G.; Ling, L.C.; Juliano, B.O. 2-acetyl-1-pyrroline: An important aroma component of cooked rice. Chem. Ind. 1982, 12, 958–959. [Google Scholar]
  24. Rathna Priya, T.; Eliazer Nelson, A.R.L.; Ravichandran, K. Nutritional and functional properties of coloured rice varieties of South India: A review. J. Ethn. Food 2019, 6, 11. [Google Scholar] [CrossRef]
  25. Suebpongsang, P.; Ekasingh, B.; Cramb, R. Commercialisation of rice farming in Northeast Thailand. In White Gold: The Commercialisation of Rice Farming in the Lower Mekong Basin; Cramb, R., Ed.; Palgrave Macmillan, Gateway East Singapore: Singapore, 2019. [Google Scholar] [CrossRef]
  26. Grandstaff, T.B.; Grandstaff, S.; Limpinuntana, V.; Suphanchaimat, N. Rainfed revolution in Northeast Thailand. Southeast Asian Stud. 2008, 46, 289–376. [Google Scholar]
  27. Padolina, T.F.; (Retired Plant Breeder, Philippine Rice Research Institute Affiliation, Science City of Muñoz, Nueva Ecija, Philippines). Personal communication, 2018.
  28. PhilRice. Pinoy Rice Knowledge Bank. 2019. Available online: https://www.pinoyrice.com/rice-varieties/ (accessed on 18 March 2019).
  29. Arocena, E.C.; (Plant Breeder, Philippine Rice Research Institute Affiliation, Science City of Muñoz, Nueva Ecija, Philippines). Personal communication, 2019.
  30. Fernandez, G.C.J. Effective selection criteria for assessing plant stress tolerance. In Proceeding of the International Symposium on Adaptation of Vegetables and Other Food Crops in Temperature and Water Stress, Shanhua, Taiwan, 13–18 August 1992. [Google Scholar] [CrossRef]
  31. Rosielle, A.A.; Hamblin, J. Theoretical aspects of selection for yield in stress and non-stress environments. Crop Sci. 1981, 21, 943–946. [Google Scholar] [CrossRef]
  32. Farshadfar, E.; Mohammadi, R.; Farshadfar, M.; Dabiri, S. Relationships and repeatability of drought tolerance indices in wheat-rye disomic addition lines. Aust. J. Crop Sci. 2013, 7, 130–198. [Google Scholar]
  33. Fischer, R.; Maurer, R. Drought Resistance in Spring Wheat Cultivars. I. Grain Yield Responses. Aust. J. Agric. Res. 1978, 29, 897–912. [Google Scholar] [CrossRef]
  34. Yan, W.; Tinker, N.A. Biplot analysis of multi-environment trial data: Principles and applications. Can. J. Plant Sci. 2006, 86, 623–645. [Google Scholar] [CrossRef]
  35. Olivoto, T.; Lúcio, A.D. Metan: An R package for multi-environment trial analysis. Methods Ecol. Evol. 2020, 11, 783–789. [Google Scholar] [CrossRef]
  36. Olivoto, T.; Lúcio, A.D.; da Silva, J.A.; Sari, B.G.; Diel, M.I. Mean Performance and Stability in Multi-Environment Trials II: Selection Based on Multiple Traits. J. Agron. 2019, 111, 2961–2969. [Google Scholar] [CrossRef]
  37. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2023; Available online: https://www.R-project.org/ (accessed on 16 March 2023).
  38. Torres, R.O.; Henry, A. Yield stability of selected rice breeding lines and donors across conditions of mild to moderately severe drought stress. Field Crops Res. 2018, 220, 37–45. [Google Scholar] [CrossRef] [PubMed]
  39. Robinson, H.F.; Comstock, R.E.; and Harvey, P.H. Breeding procedure designed to make maximum use of both general and specific combining ability. Agron. J. 1949, 41, 360–367. [Google Scholar] [CrossRef]
  40. Banayo, N.P.M.; Rahon, R.E.; Santa Cruz, P.; Kato, Y. Fertilizer responsiveness of high-yielding drought-tolerant rice in rainfed lowlands. Plant Prod. Sci. 2020, 24, 279–289. [Google Scholar] [CrossRef]
  41. Yamane, K.; Garcia, R.; Imayoshi, K.; Mabesa-Telosa, R.C.; Banayo, N.P.M.C.; Vergara, G.; Yamauchi, A.; Cruz, P.S.; Kato, Y. Seed vigour contributes to yield improvement in dry direct-seeded rainfed lowland rice. Ann. Appl. Biol. 2018, 172, 100–110. [Google Scholar] [CrossRef]
  42. Adhikari, M.; Adhikari, N.R.; Sharma, S.; Gairhe, J.; Bhandari, R.R.; Paudel, S. Evaluation of Drought Tolerant Rice Cultivars Using Drought Tolerant Indices under Water Stress and Irrigated Condition. Am. J. Clim. Change 2019, 8, 228–236. [Google Scholar] [CrossRef]
  43. Gauchan, D.; Sapkota, B.; Gautam, S.; Thapa Magar, D.B.; Sharma, B.; Amatya, S.; Sapkota, S.; Nepali, M.B.; Singh, S.; Singh, U.S. Development and Dissemination of Stress-tolerant Rice Varieties in Nepal; IRRI/NARC STRASA Project; Socioeconomics and Agricultural Research Policy Division (SARPOD), Nepal Agricultural Research Council: Khumaltar, Nepal, 2014. [Google Scholar]
  44. Dhawan, G.; Kumar, A.; Dwivedi, P.; Gopala Krishnan, S.; Pal, M.; Vinod, K.K.; Nagarajan, M.; Bhowmick, P.K.; Bollinedi, H.; Ellur, R.K.; et al. 2021. Introgression of qDTY1.1 Governing Reproductive Stage Drought Tolerance into an Elite Basmati Rice Variety “Pusa Basmati 1” through Marker Assisted Backcross Breeding. Agronomy 2021, 11, 202. [Google Scholar] [CrossRef]
  45. Waheed, R.; Ignacio, J.C.; Arbelaez, J.D.; Juanillas, V.M.; Asif, M.; Henry, A.; Kretzschmar, T.; Arif, M. Drought Response QTLs in a Super Basmati × Azucena Population by High-density GBS-based SNP Linkage Mapping. Plant Breed. 2021, 140, 758–774. [Google Scholar] [CrossRef]
  46. Henry, A.; Gowda, V.R.P.; Torres, R.O.; McNally, K.L.; Serraj, R. Variation in root system architecture and drought response in rice (Oryza sativa): Phenotyping of the Oryza SNP panel in rainfed lowland fields. Field Crops Res. 2011, 120, 205–214. [Google Scholar] [CrossRef]
  47. Xangsayasane, P.; Jongdee, B.; Pantuwan, G.; Fukai, S.; Mitchell, J.H.; Inthapanya, P.; Jothiyangkoon, D. Genotypic performance under intermittent and terminal drought screening in rainfed lowland rice. Field Crops Res. 2014, 156, 281–292. [Google Scholar] [CrossRef]
  48. Villa, J.E.; Henry, A.; Xie, F.; Serraj, R. Hybrid rice performance in environments of increasing drought severity. Field Crops Res. 2012, 125, 14–24. [Google Scholar] [CrossRef]
  49. Ahakpaz, F.; Abdi, H.; Neyestani, E.; Hesami, A.; Mohammadi, B.; Mahmoudi, K.N.; Abedi-Asl, G.; Noshabadi, M.R.J.; Ahakpaz, F.; Alipour, H. Genotype by environment interaction analysis for grain yield of barley genotypes under dryland conditions and the role of monthly rainfall, Agric. Water Manag. 2021, 245, 106665. [Google Scholar] [CrossRef]
  50. Singamsetti, A.; Shahi, J.P.; Zaidi, P.H.; Seetharam, K.; Vinayan, M.T.; Kumar, M.; Singla, S.; Shikha, K.; Madankar, K. Genotype× environment interaction and selection of maize (Zea mays L.) hybrids across moisture regimes. Field Crops Res. 2021, 270, 108224. [Google Scholar] [CrossRef]
  51. Pour-Aboughadareh, A.; Barati, A.; Koohkan, S.A.; Jabari, M.; Marzoghian, A.; Gholipoor, A.; Shahbazi-Homonloo, K.; Zali, H.; Poodineh, O.; Kheirgo, M. Dissection of genotype-by-environment interaction and yield stability analysis in barley using AMMI model and stability statistics. Bull. Natl. Res. Cent. 2022, 46, 19. [Google Scholar] [CrossRef]
  52. Raman, A.; Verulkar, S.; Mandal, N.; Variar, M.; Shukla, V.; Dwivedi, J.; Singh, B.; Singh, O.; Swain, P.; Mall, A.; et al. Drought yield index to select high yielding rice lines under different drought stress severities. Rice 2012, 5, 31. [Google Scholar] [CrossRef]
  53. Bennani, S.; Nsarellah, N.; Jlibene, M.; Tadesse, W.; Birouk, A.; Ouabbou, H. Efficiency of drought tolerance indices under different stress severities for bread wheat selection. AJCS 2017, 11, 395–405. [Google Scholar] [CrossRef]
  54. Mau, Y.S.; Ndiwa, A.S.S.; Oematan, S.S.; Markus, J.E.R. Drought tolerance indices for selection of drought tolerant, high yielding upland rice genotypes. AJCS 2019, 13, 170–178. [Google Scholar] [CrossRef]
  55. Hussain, T.; Hussain, N.; Ahmed, M.; Nualsri, C.; Duangpan, S. Responses of lowland rice genotypes under terminal water stress and identification of drought tolerance to stabilize rice productivity in Southern Thailand. Plants 2021, 10, 2565. [Google Scholar] [CrossRef] [PubMed]
  56. Kandel, B.P.; Joshi, L.P.; Sharma, S.; Adhikari, P.; Koirala, B.; Shrestha, K. Drought tolerance screening of rice genotypes in mid-hills of Nepal using various drought indices, Acta Agriculturae Scandinavica. Sect. B-Soil Plant Sci. 2022, 72, 744–750. [Google Scholar] [CrossRef]
  57. Asrat, S.; Yesuf, M.; Carlsson, F.; Wale, E.; Farmers’ Preferences for Crop Variety Traits: Lessons for on-Farm Conservation and Technology Adoption. Environment for Development Initiative. 2009. Available online: http://www.jstor.org/stable/resrep14912 (accessed on 25 August 2023).
  58. Maligalig, R.; Demont, M.; Umberger, W.J.; Peralta, A. Understanding Filipino Rice Farmer Preference Heterogeneity for Varietal Trait Improvements: A Latent Class Analysis. J. Agric. Econ. 2021, 72, 134–157. [Google Scholar] [CrossRef] [PubMed]
  59. Powell, N.; Ji, X.; Ravash, R.; Edlington, J.; Dolferus, R. Yield stability for cereals in a changing climate. Funct. Plant Biol. FPB 2012, 39, 539–552. [Google Scholar] [CrossRef] [PubMed]
  60. Stuecker, M.F.; Tigchelaar, M.; Kantar, M.B. Climate variability impacts on rice production in the Philippines. PLoS ONE 2018, 13, e0201426. [Google Scholar] [CrossRef] [PubMed]
  61. Niones, J.M.; Suralta, R.R.; Inukai, Y.; Yamauchi, A. Field evaluation on functional roles of root plastic responses on dry matter production and grain yield of rice under cycles of transient soil moisture stresses using chromosome segment substitution lines. Plant Soil 2012, 359, 107–120. [Google Scholar] [CrossRef]
  62. Suralta, R.R.; Lucob, N.B.; Perez, L.M. Shoot and root development in rice (Oryza sativa L.) genotypes during progressive drying in soils with varying moisture regimes. PJCS 2012, 37, 1–12. [Google Scholar]
  63. Kano-Nakata, M.; Gowda, V.R.P.; Henry, A.; Serraj, R.; Inukai, Y.; Fujita, D.; Kobayashi, N.; Suralta, R.R.; Yamauchi, A. Functional roles of the plasticity of root system development in biomass production and water uptake under rainfed lowland conditions. Field Crops Res. 2013, 144, 288–296. [Google Scholar] [CrossRef]
  64. RTWG (Rice Technical Working Group). National Cooperative Testing Manual for Rice: Guidelines and Policies; RTWG, National Seed Industry Council, Department of Agriculture: Quezon City, Philippines, 1997; 113p. [Google Scholar]
  65. PSA. Price Situationer of Selected Agricultural Commodities First Phase of June 2021. Available online: https://psa.gov.ph/statistics/price-situationer/selected-agri-commodities/node/164718 (accessed on 30 June 2021).
  66. Exchange Rates. 2021. Available online: https://exchangerate.org.uk (accessed on 23 June 2021).
Figure 1. Weekly soil moisture content (SMC) and water table depth (WTD), and daily, maximum and minimum temperatures, RH, solar radiation (SR), and rainfall in the rainfed and controlled-drought field conditions from transplanting to harvesting. Left side, (A,C,E), rainfed conditions of the 2019 WS, 2020 WS normal planting, and 2020 WS late planting, respectively; right side, (B,D,F) controlled-drought conditions of the 2020 DS mild drought, 2020 DS severe drought, and 2021 DS mild drought.
Figure 1. Weekly soil moisture content (SMC) and water table depth (WTD), and daily, maximum and minimum temperatures, RH, solar radiation (SR), and rainfall in the rainfed and controlled-drought field conditions from transplanting to harvesting. Left side, (A,C,E), rainfed conditions of the 2019 WS, 2020 WS normal planting, and 2020 WS late planting, respectively; right side, (B,D,F) controlled-drought conditions of the 2020 DS mild drought, 2020 DS severe drought, and 2021 DS mild drought.
Agriculture 13 01985 g001
Figure 2. Biplot view for the discriminativeness vs. representativeness (A), which-won-where (B), and grain yield x weighted average of the absolute scores (C), across the water treatments. Green diamonds represent the treatment or environment in the year–season–water treatment combination. In the treatment labels, 19, 20, and 21 indicate the years, namely 2019, 2020, and 2021, respectively; D and W denote the dry and wet season, respectively; the last letters represent the water treatments: I—irrigated, R—rainfed, RL—rainfed-late, M—mild drought, and S—severe drought. Blue circles or diamonds represent genotypes. For the released varieties with NSIC codes, the suffix numbers were retained, e.g., 192 for the drought-tolerant check variety NSIC Rc192. For the other genotypes, BR—Black rice, RR—Red rice, Cal—Calatrava, Pin—Pinilisa, BAS—Basmati, CL1—CLRice-1, CL2—CLRice-2, and CL3—CLRice-3. In discriminativeness vs. representativeness, lines connecting water treatments to the biplot origin are called environment vectors. In which-won-where, dashed line means sectors of convex hull, polygon blue line is the convex hull, while green ellipse shape represents the so-called megaenvironment.
Figure 2. Biplot view for the discriminativeness vs. representativeness (A), which-won-where (B), and grain yield x weighted average of the absolute scores (C), across the water treatments. Green diamonds represent the treatment or environment in the year–season–water treatment combination. In the treatment labels, 19, 20, and 21 indicate the years, namely 2019, 2020, and 2021, respectively; D and W denote the dry and wet season, respectively; the last letters represent the water treatments: I—irrigated, R—rainfed, RL—rainfed-late, M—mild drought, and S—severe drought. Blue circles or diamonds represent genotypes. For the released varieties with NSIC codes, the suffix numbers were retained, e.g., 192 for the drought-tolerant check variety NSIC Rc192. For the other genotypes, BR—Black rice, RR—Red rice, Cal—Calatrava, Pin—Pinilisa, BAS—Basmati, CL1—CLRice-1, CL2—CLRice-2, and CL3—CLRice-3. In discriminativeness vs. representativeness, lines connecting water treatments to the biplot origin are called environment vectors. In which-won-where, dashed line means sectors of convex hull, polygon blue line is the convex hull, while green ellipse shape represents the so-called megaenvironment.
Agriculture 13 01985 g002
Figure 3. Ranking of genotypes based on mean yield-related traits across water treatments. Number of panicles per hill (A), number of filled grains per panicle (B), 1000-grain weight (C), and harvest index (D). Green diamonds represent the treatment or environment in the year-season-water treatment combination. In the treatment labels, 19, 20, and 21 indicate the years 2019, 2020, and 2021, respectively; D and W denote the dry and wet season, respectively; the last letters represent the water treatments: I—irrigated, R—rainfed, RL—rainfed-late, M—mild drought, and S—severe drought. Blue circles or diamonds represent genotypes. For released varieties with NSIC codes, the suffix numbers were retained, e.g., 192 for the drought-tolerant check variety NSIC Rc192. For the other genotypes, BR—Black rice, RR—Red rice, Cal—Calatrava, Pin—Pinilisa, BAS—Basmati, CL1—CLRice-1, CL2—CLRice-2, and CL3—CLRice-3. The blue arrow suroounded by concentric circles is the ideal genotype.
Figure 3. Ranking of genotypes based on mean yield-related traits across water treatments. Number of panicles per hill (A), number of filled grains per panicle (B), 1000-grain weight (C), and harvest index (D). Green diamonds represent the treatment or environment in the year-season-water treatment combination. In the treatment labels, 19, 20, and 21 indicate the years 2019, 2020, and 2021, respectively; D and W denote the dry and wet season, respectively; the last letters represent the water treatments: I—irrigated, R—rainfed, RL—rainfed-late, M—mild drought, and S—severe drought. Blue circles or diamonds represent genotypes. For released varieties with NSIC codes, the suffix numbers were retained, e.g., 192 for the drought-tolerant check variety NSIC Rc192. For the other genotypes, BR—Black rice, RR—Red rice, Cal—Calatrava, Pin—Pinilisa, BAS—Basmati, CL1—CLRice-1, CL2—CLRice-2, and CL3—CLRice-3. The blue arrow suroounded by concentric circles is the ideal genotype.
Agriculture 13 01985 g003
Figure 4. Average net income (ha1) for specialty rice with the highest and lowest grain yield (A), and their comparative net income relative to the drought-tolerant (DT) check variety (NSIC Rc192) (B), under rainfed conditions in the 2019 and 2020 wet seasons (WS). The market prices for milled rice were set to USD 0.91 (PHP 45.00 kg1) for the drought-tolerant check variety [65], USD 1.46 (PHP 72.00 kg1) for aromatic rice, USD 1.56 (PHP 77.00 kg1) for pigmented rice, and USD 1.20 (PHP 59.00 kg1) for glutinous rice [21].
Figure 4. Average net income (ha1) for specialty rice with the highest and lowest grain yield (A), and their comparative net income relative to the drought-tolerant (DT) check variety (NSIC Rc192) (B), under rainfed conditions in the 2019 and 2020 wet seasons (WS). The market prices for milled rice were set to USD 0.91 (PHP 45.00 kg1) for the drought-tolerant check variety [65], USD 1.46 (PHP 72.00 kg1) for aromatic rice, USD 1.56 (PHP 77.00 kg1) for pigmented rice, and USD 1.20 (PHP 59.00 kg1) for glutinous rice [21].
Agriculture 13 01985 g004
Table 1. List of specialty rice genotypes that were used in the study.
Table 1. List of specialty rice genotypes that were used in the study.
NameCountry of
Origin
Year
Acquired
Collection
Number
Grain Yield
(t ha−1)
Remarks
NSIC Rc192 (DT check)Philippines2016CSRc-1013.7Modern variety
Aromatic
Basmati 370India2017CSRc-1233.1Traditional variety
CLRice-1Philippines2017CSRc-1153.5Improved line
NSIC Rc34Philippines2016CSRc-1024.8Modern variety
NSIC Rc218Philippines2016CSRc-1033.8Modern variety
NSIC Rc342Philippines2016CSRc-1044.3Modern variety
NSIC Rc344Philippines2016CSRc-1055.3Modern variety
Pigmented
Black ricePhilippines2016CSRc-1123.8Traditional variety
CalatravaPhilippines2016CSRc-1143.4Traditional variety
CLRice-2Philippines2017CSRc-1213.6Improved line
CLRice-3Philippines2017CSRc-1223.5Improved line
PinilisaPhilippines2016CSRc-1113.0Traditional variety
Red ricePhilippines2016CSRc-1133.2Traditional variety
Glutinous
NSIC Rc15Philippines2016CSRc-1065.4Modern variety
NSIC Rc17Philippines2016CSRc-1074.2Modern variety
NSIC Rc19Philippines2016CSRc-1084.5Modern variety
NSIC Rc21Philippines2016CSRc-1094.5Modern variety
NSIC Rc31Philippines2016CSRc-1104.7Modern variety
Table 2. Grain yield (t ha−1) of specialty rice using different water treatments.
Table 2. Grain yield (t ha−1) of specialty rice using different water treatments.
GenotypeWater TreatmentWet SeasonWater TreatmentDry Season
2019202020202021
Aromatic
NSIC Rc192 (DT check)Irrigated6.90 a5.32 cdIrrigated7.89 abc7.10 b
Rainfed6.33 ab4.84 cdefMild drought4.21 e1.34 e
Rainfed-laten.a.1.78 kSevere drought0.98 jn.a.
NSIC Rc344Irrigated6.28 ab5.22 cdeIrrigated7.50 bcd6.93 b
Rainfed6.39 ab5.02 cdefMild drought3.34 ef2.60 d
Rainfed-laten.a.3.66 ghiSevere drought1.80 ghijn.a.
NSIC Rc342Irrigated6.09 b6.30 abIrrigated8.93 a8.84 a
Rainfed5.64 bc6.78 aMild drought2.77 fg1.87 de
Rainfed-laten.a.2.54 jkSevere drought2.26 fghin.a.
NSIC Rc218Irrigated4.77 de4.72 cdefIrrigated8.28 ab6.97 b
Rainfed5.25 cd4.78 cdefMild drought2.69 fgh2.51 de
Rainfed-laten.a.2.09 kSevere drought1.61 ghijn.a.
NSIC Rc34Irrigated4.56 def5.63 bcIrrigated8.35 ab7.07 b
Rainfed5.16 cde4.12 fghMild drought3.02 ef1.79 de
Rainfed-laten.a.2.74 ijkSevere drought1.81 ghijn.a.
Basmati 370Irrigated3.93 fg5.27 cdIrrigated7.01 cd7.13 b
Rainfed4.38 ef4.28 efgMild drought2.78 fg2.17 de
Rainfed-laten.a.3.18 hijSevere drought1.25 ijn.a.
CLRice-1Irrigated3.46 g3.25 hijIrrigated6.66 d5.34 c
Rainfed3.80 fg4.47 defgMild drought1.53 hij1.58 de
Rainfed-laten.a.3.18 hijSevere drought1.72 ghijn.a.
Pigmented
NSIC Rc192 (DT check)Irrigated6.90 a5.32 bcIrrigated7.89 a7.10 b
Rainfed6.33 a4.84 cMild drought4.21 cd1.34 ef
Rainfed-laten.a.1.78 fgSevere drought0.98 hijn.a.
Black riceIrrigated4.83 bc7.13 aIrrigated7.79 a6.77 b
Rainfed5.98 ab5.68 bcMild drought3.19 de1.97 e
Rainfed-laten.a.2.78 defSevere drought2.11 fgn.a.
CLRice-2Irrigated4.68 cd6.07 abIrrigated8.01 a7.93 a
Rainfed5.91 ab5.58 bcMild drought2.45 ef1.28 ef
Rainfed-laten.a.2.14 defgSevere drought2.01 fghn.a.
Red riceIrrigated4.50 cde5.12 bcIrrigated6.18 b5.77 c
Rainfed4.68 cd5.04 bcMild drought2.21 ef 1.65 e
Rainfed-laten.a.2.57 defgSevere drought1.13 ghin.a.
CalatravaIrrigated3.38 ef4.52 cIrrigated5.24 bc4.63 d
Rainfed3.53 def3.23 dMild drought1.99 fgh1.39 ef
Rainfed-laten.a.1.45 gSevere drought1.63 fghin.a.
CLRice-3Irrigated3.23 f2.61 defgIrrigated5.15 bc5.14 cd
Rainfed3.17 f3.22 dMild drought2.30 ef1.28 ef
Rainfed-laten.a.1.49 gSevere drought1.54 fghin.a.
PinilisaIrrigated0.86 g4.62 cIrrigated3.52 d7.09 b
Rainfed1.52 g3.12 deMild drought0.89 ij0.87 f
Rainfed-laten.a.1.93 efgSevere drought0.00 jn.a.
Glutinous
NSIC Rc192 (DT check)Irrigated6.90 a5.32 aIrrigated7.89 a7.10 b
Rainfed6.33 ab4.84 abMild drought4.21 c1.34 cd
Rainfed-laten.a.1.78 eSevere drought0.98 en.a.
NSIC Rc15Irrigated4.56 cd5.54 aIrrigated7.67 a9.50 a
Rainfed6.85 a5.42 aMild drought2.05 d2.47 c
Rainfed-laten.a.3.30 cdSevere drought1.76 den.a.
NSIC Rc21Irrigated4.32 cd5.04 abIrrigated7.77 a7.35 b
Rainfed5.22 bc4.91 abMild drought1.62 de1.66 cd
Rainfed-laten.a.2.95 deSevere drought0.85 en.a.
NSIC Rc19Irrigated4.23 cd5.43 aIrrigated6.33 b7.50 b
Rainfed5.09 c4.61 abMild drought1.65 de1.06 d
Rainfed-laten.a.2.95 deSevere drought1.29 den.a.
NSIC Rc31Irrigated3.50 d4.88 abIrrigated7.08 ab6.89 b
Rainfed4.83 c4.68 abMild drought2.04 d1.93 cd
Rainfed-laten.a.2.69 deSevere drought1.02 en.a.
NSIC Rc17Irrigated3.49 d4.41 abcIrrigated6.12 b6.66 b
Rainfed4.65 cd3.88 bcdMild drought2.03 d1.42 cd
Rainfed-laten.a.2.69 deSevere drought1.75 den.a.
Values followed by the same letter in a column within the specialty rice group are not significantly different (p < 0.05). Abbreviations: n.a.: data not available due to extremely low seed setting, DT: drought tolerant.
Table 3. Drought tolerance indices of specialty rice varieties and check under rainfed and controlled drought conditions.
Table 3. Drought tolerance indices of specialty rice varieties and check under rainfed and controlled drought conditions.
VarietyRainfedControlled Drought
GMPMPHMISTITOLSSIGMPMPHMISTITOLSSI
NSIC Rc192 (DT check)5.13 (4)5.21 (4)5.06 (5)1.35 (4)1.79 (18)2.10 (17)3.65 (4)4.14 (4)3.21 (4)0.22 (3)3.93 (19)0.13 (1)
Aromatic
NSIC Rc3445.38 (2)5.39 (3)5.36 (2)1.48 (2)0.73 (13)0.91 (7)3.85 (1)4.17 (3)3.56 (1)0.25 (1)3.17 (12)0.13 (1)
NSIC Rc3425.56 (1)5.59 (1)5.52 (1)1.58 (1)1.21 (17)1.40 (14)3.78 (3)4.25 (1)3.36 (3)0.24 (2)3.89 (18)0.13 (1)
NSIC Rc2184.38 (11)4.39 (12)4.37 (12)0.98 (10)0.70 (11)1.06 (10)3.28 (7)3.51 (8)3.07 (7)0.18 (5)2.48 (6)0.15 (3)
NSIC Rc344.52 (7)4.55 (8)4.49 (9)1.04 (7)1.09 (15)1.53 (16)3.35 (5)3.65 (5)3.08 (5)0.19 (4)2.89 (8)0.15 (3)
Basmati 3704.26 (12)4.27 (13)4.25 (13)0.93 (11)0.66 (9)1.14 (13)3.09 (10)3.34 (9)2.85 (9)0.16 (7)2.53 (7)0.16 (4)
CLRice-13.49 (14)3.49 (15)3.48 (15)0.62 (13)−0.27 (1)−0.57 (1)2.33 (16)2.48 (16)2.18 (15)0.09 (12)1.74 (2)0.19 (7)
Pigmented
Black rice5.36 (3)5.40 (2)5.33 (3)1.47 (3)1.17 (16)1.40 (14)3.80 (2)4.20 (2)3.45 (2)0.24 (2)3.56 (17)0.13 (1)
CLRice -24.94 (6)4.96 (7)4.92 (6)1.25 (6)0.83 (14)1.11 (12)3.21 (9)3.65 (6)2.82 (10)0.17 (6)3.46 (15)0.14 (2)
Red rice4.44 (9)4.45 (10)4.42 (10)1.01 (8)0.72 (12)1.07 (11)2.83 (11)3.24 (10)2.47 (11)0.14 (8)3.15 (11)0.16 (4)
Calatrava3.29 (15)3.34 (16)3.23 (16)0.55 (14)1.21 (17)2.21 (18)2.57 (14)2.81 (15)2.35 (14)0.11 (11)2.28 (5)0.18 (6)
CLRice -32.77 (16)2.77 (17)2.76 (17)0.39 (15)0.29 (4)0.72 (5)2.23 (17)2.31 (17)2.16 (16)0.08 (13)1.21 (1)0.19 (7)
Pinilisa2.45 (17)2.46 (18)2.43 (18)0.31 (16)0.55 (7)1.45 (15)1.27 (18)1.66 (18)0.96 (19)0.03 (14)2.15 (3)0.25 (10)
Glutinous
NSIC Rc155.12 (5)5.12 (5)5.12 (4)1.34 (5)−0.14 (2)−0.20 (2)3.25 (8)3.57 (7)2.96 (8)0.18 (5)2.96 (9)0.15 (3)
NSIC Rc214.52 (7)4.52 (9)4.51 (7)1.04 (7)0.32 (5)0.49 (4)2.54 (15)3.03 (13)2.13 (17)0.11 (11)3.30 (14)0.16 (4)
NSIC Rc194.51 (8)4.52 (9)4.50 (8)1.04 (7)0.61 (8)0.91 (8)2.54 (15)3.08 (12)2.09 (18)0.11 (11)3.49 (16)0.16 (4)
NSIC Rc314.40 (10)4.41 (11)4.39 (11)0.99 (9)0.69 (10)1.04 (9)2.81 (12)3.21 (11)2.46 (12)0.13 (9)3.09 (10)0.16 (4)
NSIC Rc173.84 (13)3.84 (14)3.84 (14)0.76 (12)0.21 (3)0.38 (3)2.62 (13)2.84 (14)2.41 (13)0.12 (10)2.22 (4)0.17 (5)
Abbreviations: GMP: geometric mean productivity, MP: mean productivity index, HMI: harmonic mean index, STI: stress tolerance index, TOL: tolerance index, SSI: stress susceptibility index, DT: drought tolerant. The numbers in parentheses indicate the rank of the genotype in each index.
Table 4. Head rice recovery of specialty rice and drought tolerance check under rainfed conditions in the 2019 and 2020 wet seasons.
Table 4. Head rice recovery of specialty rice and drought tolerance check under rainfed conditions in the 2019 and 2020 wet seasons.
VarietyHead Rice Recovery (%)
19WS20WSMean
NSIC Rc192
(Drought tolerant check)
51.94 abc54.38 ab53.16
Aromatic
NSIC Rc34452.99 abc53.97 ab53.48
NSIC Rc34249.60 abcd53.02 ab51.31
NSIC Rc21855.72 a49.00 abc52.36
NSIC Rc3454.25 ab52.69 abc53.47
Basmati 37040.17 de47.32 bcd43.74
CLRice-134.31 ef38.46 de36.39
Pigmented
Black rice40.18 de42.66 cd41.42
CLRice-229.62 f32.64 ef31.13
Red rice48.26 abcd47.39 bcd47.83
Calatrava46.08 bcd47.22 bcd46.65
CLRice-331.86 ef33.62 ef32.74
Pinilisa43.96 cd31.30 f37.63
Glutinous
NSIC Rc1553.12 abc58.47 a55.79
NSIC Rc2153.90 ab55.75 ab54.82
NSIC Rc1950.45 abc52.34 abc51.40
NSIC Rc3146.21 abcd52.74 abc49.47
NSIC Rc1746.19 abcd58.54 a52.37
Year (Y)ns
Genotype (G)**
Y × Gns
Values followed by the same letter in a column are not significantly different (p < 0.05). ** and ns represent highly significant and not significant, respectively.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Agustin, A.M.L.; Ordonio, J.L.; Natividad, M.B.S.; Lucob-Agustin, N.B.; Suralta, R.R.; Ehara, H.; Mitsuya, S.; Kano-Nakata, M. Specialty Rice (Oryza sativa L.) with High and Stable Grain Yield under Rainfed Lowland Conditions. Agriculture 2023, 13, 1985. https://doi.org/10.3390/agriculture13101985

AMA Style

Agustin AML, Ordonio JL, Natividad MBS, Lucob-Agustin NB, Suralta RR, Ehara H, Mitsuya S, Kano-Nakata M. Specialty Rice (Oryza sativa L.) with High and Stable Grain Yield under Rainfed Lowland Conditions. Agriculture. 2023; 13(10):1985. https://doi.org/10.3390/agriculture13101985

Chicago/Turabian Style

Agustin, Ace Mugssy L., Jeremias L. Ordonio, Marie Bie S. Natividad, Nonawin B. Lucob-Agustin, Roel R. Suralta, Hiroshi Ehara, Shiro Mitsuya, and Mana Kano-Nakata. 2023. "Specialty Rice (Oryza sativa L.) with High and Stable Grain Yield under Rainfed Lowland Conditions" Agriculture 13, no. 10: 1985. https://doi.org/10.3390/agriculture13101985

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop