Rice is the main staple food and source of income for most people in the world, particularly in Asia and Africa. In 2014, the world’s rice production was estimated at 750.9 million tons, of which 678.8 million tons and 28.3 million tons were harvested in Asia and Africa, respectively (Food and Agriculture Organization 2014).
Rice is generally cultivated in irrigated conditions for ease of weed control and better establishment. Declining water availability combined with erratic climatic conditions (El Niño) is becoming a major concern that threatens the sustainability of irrigated rice ecosystems (Rahman et al. Reference Rahman, Juraimi, Suria, Man and Anwar2012). Thus, the aerobic rice system and water-saving techniques are being developed and promoted to mitigate the problem of water scarcity (Belder et al. Reference Belder, Bouman, Spiertz, Peng and Castañeda2005; Bouman Reference Bouman2001). However, weed problems present a major constraint that contributes to the greatest yield loss in aerobic rice systems (Chauhan et al. Reference Chauhan, Opeña and Ali2015; Jabran and Chauhan Reference Jabran and Chauhan2015; West Africa Rice Development Association 1996). Aerobic soil conditions are highly conducive to weed seed germination and growth, which result in greater weed pressure and higher yield losses than in flood-irrigated rice (Balasubramanian and Hill Reference Balasubramanian and Hill2002; Rao et al. Reference Rao, Johnson, Sivaprasad, Ladha and Mortimer2007). Failure to control weeds in aerobic rice often results in very low or zero yield (Johnson Reference Johnson2009). In sub-Saharan Africa, weeds lead to an estimated rice yield loss of 2.2 million tons per year (Rodenburg and Johnson Reference Rodenburg and Johnson2009), whereas in the tropics, yield losses are estimated to be at 35% (Oerke and Dehne Reference Oerke and Dehne2004).
The most common control measure against weeds is hand weeding, which is laborious and expensive. Hand-weeding operations (two to three times per cropping season) comprise 15% of the total farming operation cost (VanDevender et al. Reference VanDevender, Costello and Smith1997). Nowadays, the use of herbicides to control weeds is considered to be the most practical, effective, and economical approach (Rahman et al. Reference Rahman, Juraimi, Suria, Man and Anwar2012). However, the intensive use of herbicides has led to negative results, such as herbicide resistance, environmental contamination, and weed shifts toward tolerant ecotypes (Carey et al. Reference Carey, Hoagland and Talbert1995; Fisher et al. Reference Fisher, Granados and Trojillo1993; Lemerle et al. Reference Lemerle, Gill, Murphy, Walker, Cousens, Mokhtari, Peltzer, Coleman and Luchett2001; Valverde et al. Reference Valverde, Riches and Caseley2000).
Breeding for weed-competitive rice cultivars (with high yield ability) can be an important strategy for reducing hand weeding and herbicide inputs (Zhao et al. Reference Zhao, Atlin, Bastiaans and Spiertz2006a). However, limited studies have been carried out to evaluate weed competitiveness in various crops, including corn (Zea mays L.) (Silva et al. Reference Silva, Silva, Silva, Oliveira and Ferreira2010), sorghum [Sorghum bicolor (L.) Moench] (Wu et al. Reference Wu, Walker, Osten and Robinson2010), canola (Brassica napus L.) (Lemerle et al. Reference Lemerle, Lockley, Luckett and Wu2010), soybean [Glycine max (L.) Merr.] (Jannink et al. Reference Jannink, Orf, Jordan and Shaw2000; Place et al. Reference Place, Reberg-Horton, Dickey and Carter2011), wheat (Triticum aestivum L.) (Huel and Hucl Reference Huel and Hucl1996; Lanning et al. Reference Lanning, Talbert, Martin, Blake and Bruckner1997; Lemerle et al. Reference Lemerle, Verbe, Cousens and Coombes1996, Reference Lemerle, Gill, Murphy, Walker, Cousens, Mokhtari, Peltzer, Coleman and Luchett2001), and rice (Moukoumbi et al. Reference Moukoumbi, Sie, Vodouhe, Bonou, Toulou and Ahanchede2011; Namuco et al. Reference Namuco, Cairns and Johnson2009; Zhao et al. Reference Zhao, Atlin, Bastiaans and Spiertz2006a, Reference Zhao, Atlin, Bastiaans and Spiertz2006b).
Weed-competitive ability is a complex trait that cannot be attributed to a single characteristic but is a result of the interaction among several desirable traits (Caton et al. Reference Caton, Mortimer, Foin, Hill, Gibson and Fischer2000; Chauhan et al. Reference Chauhan, Opeña and Ali2015; Kruepl et al. Reference Kruepl, Hoad, Davies, Bertholdsson and Paolini2007), thus making it difficult for plant breeders to breed for weed-competitive crop cultivars. Consequently, it is imperative to determine the complexity of the mechanisms and agro-morphological traits that confer weed competitiveness to fast-forward breeding efforts for weed-competitive rice cultivars. This can be done by identifying suitable germplasm for weed competitiveness and standardizing the screening protocol for phenotyping traits related to weed competitiveness.
This study was conducted to identify suitable donors and recurrent parents for weed competitiveness, identify and standardize the traits related to weed competitiveness, and identify promising introgression lines (ILs) with superior and competitive traits.
Materials and Methods
Place and Date of Study
Laboratory, greenhouse (NG-04), and field experiments (lowland/irrigated: 14.23°N, 121.43°E; upland/dry seeded: 14.22°N, 121.43°E and 23 m elevation) were conducted at the International Rice Research Institute (IRRI) Los Baños, Laguna, Philippines, during the 2012 dry season (DS), 2013 wet season (WS), and 2014 DS. The soil type on the experimental farm is Maahas clay loam (isohyperthermic mixed Typic Tropudalf).
Selection of Plant Materials
Four early-generation, backcross populations (BC1F2) derived from one common recipient parental varietal background of Weed Tolerant Rice 1 (WTR1) and four different donor parents (Y134, Zhong 143, Khazar, and Cheng Hui-448) were tested to identify a suitable donor parent and recipient parent for weed competitiveness. The germplasm lines were selected based on grain yield performance and selection criteria for weed competitiveness (Table 1). The cultivars used as a control in the investigation of weed-competitive ability traits are presented in Table 2.
a Conditions: L, low-input; RF, rainfed conditions; Y, irrigated.
a Abbreviations: GSR, Green Super Rice; IRRI, International Rice Research Institute; NCT, National Rice Cooperative Tests.
Selection Scheme
In the 2011 WS (before this study), four early-generation backcross populations (BC1F2) derived from one common recipient background (WTR1) and four different donor parents (Y134, Zhong 143, Khazar, and Cheng Hui-448) were initially tested to identify suitable donor and recipient parents for weed competitiveness. These BC1F2 populations were designated by a Green Super Rice (GSR) name as follows: G-4 (WTR1/Cheng Hui-448//WTR1), G-6 (WTR1/Y134//WTR1), G-7 (WTR1/Zhong 413//WTR1), and G-8 (WTR1/Khazar//WTR1). These BC1F2 populations were evaluated and selected under irrigated (Y), rainfed (RF), and low-input (L) conditions to produce the BC1F3 generation. Promising ILs were selected from each population by single plant selection. The G-4 population has 46 promising ILs: 16 (Y), 14 (L), and 16 (RF); G-6 has 48 ILs: 16 (Y), 15 (L), and 17 (RF); G-7 has 43 ILs: 14 (Y), 15 (L), and 14 (RF); and G-8 has 47 ILs: 16 (Y), 15 (L), and 16 (RF). In total, 184 BC1F3 ILs were used for the 2012 DS selection process under three conditions: upland weed-free (UPWF), upland weedy (UPW), and lowland weedy (LLW). In the 2013 WS, the BC1F4 populations were reduced to G-4 (138 ILs) and G-6 (144 ILs) (Figure 1). Ultimately, the G-6 population was advanced to the 2014 DS and was used to standardize the evaluation criteria for weed-competitive traits in rice.
Selection intensity (SI) is the number of single plants selected from a given population size of segregating plants planted in a condition and is expressed as a percent. G-6 BC1F2 population size in the 2011 WS was 576 to create an SI of 2.8% (Y), 2.6% (L), and 3.0% (RF) (Figure 2). In the 2012 DS, SI was 3.6% from a population size of 1,344 across all conditions, whereas in the 2013 WS, SI was 2.2% in UPW (89 ILs) and 2.7% in LLW (107 ILs) (Figure 2). In the 2014 DS, 196 ILs from the G-6 BC1F5 population were planted in UPW, UPWF, and LLW and two batches were grown under greenhouse conditions to determine their weed-competitive ability and identify promising ILs to be developed as competitive rice varieties.
Experimental Design and Treatments
BC1F3 to BC1F5 populations from four different donor parents derived in one common WTR1 recipient background (subplots) were tested for weed competitiveness under main plots: UPWF, UPW, and LLW conditions. The experimental design is an augmented randomized complete block design (RCB) with two replications (due to the limited number of seeds). A greenhouse experiment, seed germination test, and seedling vigor test (replicated twice in an RCBD) were also conducted in the 2014 DS to strengthen selection for weed-competitive traits.
Field Experiments
The G populations with a seeding rate of 90 viable seeds were manually drilled in four rows (7 seeds row−1 in 20-cm row spacing by 15-cm plant-to-plant spacing). The UPWF conditions were treated with the PRE herbicide oxadiazon (Ronstar EC25, 0.75 kg ai ha−1, Bayer CropScience, Canlubang Industrial Estate, Laguna, Philippines) using a backpack sprayer (TeeJet® 200015, TeeJet Technologies, Wheaton, IL) mounted with a flat-fan nozzle calibrated to deliver 160 L ha−1 at 138 kPa. Herbicide applications were made just after seeding. Plots were irrigated using a sprinkler immediately after herbicide application and were maintained weed-free during the critical period of competition. The LLW and UPW conditions were completely hand weeded only once at 21 d after sowing (DAS). Weeds were allowed to grow thereafter. To augment and homogenize the natural weed population in the LLW trial, barnyardgrass [Echinochloa crus-galli (L.) Beauv.] seeds were broadcast at a rate of 15 kg ha‒1, whereas in the upland weedy trial, itchgrass [Rottboellia cochinchinensis (Lour.) W. D. Clayton] seeds were broadcast at a rate of 10 kg ha‒1 in the 2012 DS field experiment. Weed seeds were broadcast at 28 d after rice was sown.
In each field, a compound N–P–K fertilizer (14:14:14) was applied at 200 kg ha−1, and two additional applications of urea were topdressed each at 60 kg ha−1 at 28 and 56 DAS, respectively. The UPWF and UPW fields were irrigated regularly with an overhead sprinkler and were maintained under nonsaturated aerobic conditions to facilitate the growth of weeds.
Rice agro-morphological traits, yield, and weed-related component traits were evaluated for weed competitiveness under field conditions. For plant height and tiller number, five random plants were measured for each IL. For yield traits, three random plants from each IL were selected and measured for panicle number, flag-leaf length, and flag-leaf width. Also, the three selected plants were threshed and separated into filled and unfilled grains, which were then counted and weighed.
The predominant weed species under weedy conditions (UPW and LLW) were determined by computing the summed dominance ratio (SDR) of each species based on weed density, relative density, biomass, relative biomass, and frequency per square meter. Three samples (1-m2 quadrats) were taken in UPW and two in LLW conditions.
The most dominant weed species in UPW conditions based on SDR values were purple nutsedge (Cyperus rotundus L.) (53.06%), goosegrass [Eleusine indica (L.) Gaertn.] (44.63%), guineagrass [Urochloa maxima (Jacq.) R. Webster] (37.97%), eclipta [Eclipta prostrata (L.) L.] (37.47%), and knotgrass (Paspalum distichum L.) (26.97%). The most dominant weed species in LLW conditions based on SDR were fimbry (Fimbristylis littoralis Gaudich.) (83.89%), junglerice [Echinochloa colona (L.) Link] (43.14%), Chinese sprangletop [Leptochloa chinensis (L.) Nees] (34.96%), eclipta (34.22%), and guineagrass (34.01%). The externally seeded weed species barnyardgrass had a 17.8% SDR value, but no data were collected for itchgrass, because random sampling (sample size of 3) was done across the experimental plot.
Greenhouse Experiments
Two greenhouse experiments were conducted to determine the traits related to weed competitiveness in the G-6 population. In the first experiment, two replications of 5 seeds of G-6 ILs (BC1F5 population), parents WTR1 and Y134, and checks PSB Rc82, NSIC Rc222, NSIC Rc192, GSR IR1-8-S6-S3-Y2, IR74371-70-1-1, and Apo were direct seeded in plastic pots filled with soil, thinned to 1 seedling pot−1, and maintained until grain maturity. One hundred seeds of junglerice were sown in each pot simultaneously with rice seeds to homogenize weed emergence. The initial germination rate of each IL, parents, and checks was measured at 7 DAS. Plant height of rice seedlings was recorded at 7, 14, 21, and 28 DAS. At maturity, data pertaining to maximum plant height, tiller number, panicle fresh weight, flag-leaf length, flag-leaf width, leaf area, number of panicles, number of filled and unfilled grains, spikelet fertility, and grain weight of each IL, parents, and checks were collected.
In the second experiment, two replications of 5 seeds of G-6 ILs (BC1F5 population), parents WTR1 and Y134, and checks were direct seeded and maintained in plastic pots filled with field soil for 28 DAS to investigate traits related to seedling vigor. One hundred seeds of junglerice were sown in each pot simultaneously with rice seeds to homogenize weed emergence. The initial number of germinated seeds and percent germination of each IL, parents, and checks were gathered. Rice seedling height and number of leaves were measured at 7, 14, 21, and 28 DAS. At 28 DAS, the number of tillers and leaf chlorophyll content (Soil-Plant Analyses Development [SPAD] meter readings) were recorded. After each observation, the rice seedlings were pulled out carefully for measurement of shoot and root length. The lengths of shoots and roots were measured from the collar region down to the tip of the longest root and shoot of each seedling. The leaves and roots of each seedling were separated for leaf fresh and dry weights, root fresh and dry weights, and total fresh and dry weights. A vigor index was computed by multiplying percent germination by seedling total dry weight (Diwan Reference Diwan2006). Weed density, fresh weight, and dry weight data were also collected to correlate with seedling vigor performance as related to weed competitiveness. Weed biomass was clipped at the soil surface from each pot, weighed, and oven-dried at 70 C for 5 d to obtain dry weight data.
Seed Germination Tests
Seed dormancy was broken by subjecting the seeds to 50 C temperature for 4 d (Jennings and de Jesus Reference Jennings and de Jesus1964). Two replications of 25 seeds were germinated in a 9-cm-diameter petri dish lined with a filter paper and placed in a germination chamber for 14 d. Seeds showing 2-mm radicle length were considered germinated. Seeds that germinated after 48 h were counted and recorded as first-count germination. Seeds that germinated after 7 d were counted and recorded as second-count germination. The rate of germination was calculated as the ratio of the final count of germination to the total number of initial seeds and expressed as a percentage. On the day 14, leaf and root lengths were recorded from five normal germinated seeds randomly selected from each replication. Leaf length was measured from the collar region to the tip of the longest leaf and expressed in centimeters. Root length was measured from the collar region down to the tip of the longest root. Average fresh and dry weights were derived from the total fresh and dry weights of germinated seeds (in milligrams) divided by the number of germinated seeds. Vigor index was calculated by multiplying the rate of germination by total dry weight.
Statistical Analysis
Statistical analysis was performed using R v. 1.5 (R Core Team 2016) and PBTools 1.4 (IRRI 2017). Descriptive statistics, correlation analysis, and ANOVA were performed on all data collected (Table 3). Pearson’s correlation test was carried out to assess possible variable linear relationships across weedy conditions. Analysis of variance was performed to verify whether the weed-competitive ability traits of at least one of the ILs varies with weed regimes. The ILs were analyzed for the significance of genotypic effect per weedy condition. Pairwise mean comparison using Fisher’s least significant different test was used to analyze the yield differences among the top-performing ILs, parents, and check cultivars. Regression analysis was carried out to determine which traits were significantly affecting the grain yield of the newly developed ILs tested under weedy conditions. Using R v. 1.5, the regression model was selected using all possible regression models with leaps and the coefficient of variation (R2) as the criterion.
Results and Discussion
Novel Weed-Competitive Rice-Breeding Strategy
Breeding for weed-competitive rice is becoming an important area of interest to plant breeders and weed scientists because of the increasing weed problems in aerobic rice ecosystems, adverse effects of extensive herbicide use, increasing reports of herbicide resistance, herbicide contamination of the environment, lack of labor for hand weeding, and inaccessibility to herbicides in marginal areas. To overcome this, we need to breed weed-competitive cultivars that could suppress or tolerate weeds without compromising grain yield. Weed-competitive cultivars will complement other weed control methods by suppressing weeds while maintaining acceptable yields (Worthington and Reberg-Horton Reference Worthington and Reberg-Horton2013). Successful breeding of weed-competitive rice cultivars will significantly reduce rice production costs in aerobic rice systems through a reduction in herbicide use and hand weeding. This strategy could further augment water-saving techniques to reduce production costs and improve yields in upland rice production.
Breeding programs for cultivars that possess weed-competitive ability will be advanced by satisfactory evaluation and germplasm screening for weed competitiveness (Kruepl et al. Reference Kruepl, Hoad, Davies, Bertholdsson and Paolini2007). To achieve success in germplasm screening, evaluation criteria should be measured in a nondestructive and rapid process (Caton et al. Reference Caton, Cope and Mortimer2003). Moreover, the key traits that confer weed competitiveness need to be correctly identified and elucidated (Gill and Coleman Reference Gill and Coleman1999).
Selection for weed-competitive ability can be implemented directly (under weedy conditions) or indirectly (under weed-free conditions) for secondary traits related to weed competitiveness (Zhao et al. Reference Zhao, Atlin, Bastiaans and Spiertz2006a). Traits under weedy and weed-free conditions can be considered as correlated traits, expressed by a single genotype in two diverse environments (Zhao Reference Zhao2006). All traits gathered under weed-free conditions were closely correlated with the same traits gathered under weedy conditions (Zhao et al. Reference Zhao, Atlin, Bastiaans and Spiertz2006b). Moreover, traits that are considered useful in indirect selection for weed competitiveness should be highly correlated with yield in weedy conditions and practical for use in large breeding populations to attain adequate SI (Atlin et al. Reference Atlin, Cooper and Bjornstad2001). Indirect selection for weed competitiveness under weed-free conditions permits germplasm selection to be implemented at the beginning of the breeding program, while direct selection under weedy conditions can be employed only in the later phase, when a satisfactory number of seeds are available (Wall Reference Wall1983). Based on these criteria, the candidate populations in this experiment were selected through direct (under UPW and LLW) and indirect (under UPWF) selection procedures with satisfactory SI. The weed-competitive traits in rice were further confirmed and standardized in a controlled greenhouse experiment for seedling vigor tests and analyzed for correlation with weed density, weed fresh weight, and weed dry weight.
Severe SI is extremely useful for identifying and selecting elite ILs only in each stage and in different conditions as compared with the performance of both the parents and check varieties. In our breeding program on weed competitiveness, we ensured a stronger SI of less than 3.6% under all conditions and stages. SI of 3.6% is the percentage of the number of single plants selected (48) from a given population size of segregating plants (1,344) planted in each weedy condition. Each elite selection for weed competitiveness was based on seedling vigor, rapid growth, effective tiller number, and single-plant grain yield. The parent Y134 (donor) and WTR1 (recipient) were identified as the most suitable parents for weed competitiveness. A single backcross population (G-6) from a cross between Y134 and WTR1 was identified as the best population, as it expressed both weed competitiveness and high yield potential under weedy conditions in lowland and upland environments.
Among the selections, G-6-L2-WL-3, G-6-RF6-WL-3,G-6-L15-WU-1, G-6-Y16-WL-2, and G-6-L6-WU-3 were the top ILs in LLW conditions, whereas G-6-Y7-WL-3, G-6-Y6-WU-3, G-6-Y3-WL-3, and G-6-Y8-WU-1 were the highest-yielding ILs under upland weedy conditions.
Selection Pressure by Environment
Initial screening under Y, RF, and L conditions allowed the selection of 48 promising G-6 ILs that gave relatively higher BC1F2 bulk population grain weights under both rainfed and irrigated environments. This correlates to their desirable response to indirect selection for weed competitiveness while having desirable yield potential. G-7 ILs gave higher BC1F2 bulk population grain weights under irrigated low-input environments, while G-8 ILs had higher BC1F2 bulk population grain weights under rainfed conditions only. However, the ILs of G-7 and G-8 did not perform well when tested under UPWF, UPW, and LLW conditions in the 2012 DS. Thus, these two populations were dropped in the 2013 WS. G-4 ILs gave higher BC1F2 bulk population grain weights under irrigated conditions only, signifying that these ILs have no desirable response to indirect selection for weed-competitive ability. When tested under UPWF, UPW, and LLW conditions in the 2012 DS, both G-4 and G-6 performed well; thus, selected ILs from these populations were advanced for further screening in the 2013 WS. Finally, G-6 ILs were selected as the best population demonstrating weed-competitive ability under LLW and UPW conditions and expressing higher grain yield potential than both the parents and inbred checks.
Performance of G-6 ILs in UPWF, UPW, and LLW Conditions in the 2013 WS
Figure 3 shows the grain yield of the top 30 G-6 ILs, WTR1 (recipient parent), Y134 (donor parent), and checks PSB Rc82 and Apo in the 2013 WS. In UPW conditions, the ILs had the highest grain yield, but there was no significant difference in the grain yields of WTR1, Y134, and the checks, which all performed significantly below the ILs. The ILs from the G-6 populations were better weed competitors than both the parents and the checks. Thorough study of the different traits of these identified ILs will help elucidate the key traits for weed competitiveness in rice.
In UPWF conditions, the checks surpassed the grain yields of both the parents and ILs. This is expected, since the checks were adapted to weed-free environments. However, with the presence of weeds in upland conditions, the yield of the checks and both WTR1 and Y134 decreased significantly, with greater yield reductions in the checks compared with WTR1 or Y134. The checks were sensitive to weeds; therefore, their yields declined drastically under weedy conditions.
In LLW conditions, the grain yield of the ILs was greater than that of both the parents and checks, although significantly lower than the grain yields obtained under upland conditions. The checks still outperformed both parents (WTR1 and Y134), while the recipient parent (WTR1) had the lowest grain yield.
There was no significant difference in the grain yields of ILs subjected to UPWF and UPW conditions in the 2013 WS (Table 4). The ILs had strong weed-competitive ability and were able to produce comparable yields even with the presence of weeds. This observation was further evaluated in field and greenhouse experiments in the 2014 DS.
a Weedy conditions with the same letters are not significant at the 5% level.
Performance of G-6 ILs in UPWF, UPW, and LLW Conditions in the 2014 DS
Figure 4 shows the grain yield of the top 30 G-6 ILs, WTR1 (recipient parent), Y134 (donor parent), and checks PSB Rc82 and Apo in the 2014 DS. In UPWF conditions, the ILs had the highest grain yield, whereas Y134 produced only 0.2 g grain yield plant−1 due to bacterial blight disease. The grain yields of the recurrent parent WTR1 and the checks showed nonsignificant differences.
In UPW conditions, the grain yields of the ILs were significantly lower in the presence of weeds. As expected, the selected ILs had the highest grain yield, whereas Y134 had the lowest yield. The grain yields of WTR1 and the checks did not differ and were significantly lower than those of the ILs. In LLW conditions, the grain yield of the ILs declined significantly. Recurrent parent WTR 1 had the highest grain yield, whereas Y134 had the lowest. However, the grain yield of Y134 did not differ from that of the check cultivars. There was a significant difference in the grain yields of ILs subjected to UPWF, UPW, and LLW conditions (Table 4). The highest grain yields were observed in UPWF conditions. However, even with the presence of weeds under UPW conditions, the ILs produced high yields, suggesting that the selected ILs have a strong weed-competitive ability. Consequently, even with less hand weeding or herbicide application, the yields of the ILs will just be the same or higher. In the long run, this is efficient in terms of lower labor cost, less herbicide use, and more profit. However, the ILs were not performing well in LLW conditions. The selected ILs seemed to be more adapted to upland conditions, and the selection for weed-competitive cultivars is stronger for UPW conditions. The first few cycles of selection under LLW conditions and using parental lines that perform well under LLW conditions may be important. These UPW and LLW conditions are two extreme conditions, and selection needs to be done separately in a more systematic manner involving crosses with ideally suitable parental lines.
The different agronomic and competitive traits of these ILs from the 2014 DS field experiment were further investigated under controlled environments including a greenhouse experiment, a seedling vigor test, and a germination test.
Key Traits for Weed-Competitive Ability
Correlation Analysis for Weed-Competitive Ability Traits under Field Conditions in the 2013 WS
High yield capacity despite the presence of weeds is a strong indication of weed-competitive ability in cultivars. In the 2013 WS, grain yield was regarded as the key factor in considering selection for weed competitiveness. Correlation analysis showed that tiller number (P<0.001; r=0.49) and panicle number (P<0.001; r=0.64) are significantly and positively correlated with grain yield (Table 5). This indicates that ILs with more tillers and panicles tend to have higher yield. In contrast, plant height (r=−0.02), leaf chlorophyll content (r=−0.06), or panicle length (r=−0.13) did not show a significant correlation (P>0.05) with grain yield. Plant height is a key trait for weed competitiveness (Lanning et al. Reference Lanning, Talbert, Martin, Blake and Bruckner1997; Lemerle et al. Reference Lemerle, Verbe, Cousens and Coombes1996). However, the result for the 2013 WS showed that taller plants will not necessarily give higher grain yield. Fischer et al. (Reference Fischer, Ramirez and Lozano1997, Reference Fischer, Ramirez, Gibson and Pinheiro2001) reported that there is no clear association between plant height and weed-suppressive ability in irrigated and upland rice. Similarly, having higher leaf chlorophyll content (SPAD reading) and longer panicles is not associated with grain yield. Moukoumbi et al. (Reference Moukoumbi, Sie, Vodouhe, Bonou, Toulou and Ahanchede2011) had observed weed-competitive cultivars showing a SPAD unit of less than 30, whereas in this study, results of the 2013 WS experiment showed a SPAD mean unit of 37.12. Negative correlation exists between panicle number and panicle length, but tiller number is positively correlated with panicle number (P<0.001; r=0.55) and negatively correlated with panicle length (P<0.05; r=−0.26). Conversely, leaf chlorophyll content (SPAD reading) is positively correlated with panicle length (P<0.05; r=0.23).
a Significance codes: * 0.05≥P≥0.01; ** 0.01≥P≥0.001; *** P≤0.001.
Correlation Analysis for Weed-Competitive Ability Traits under Field Conditions in the 2014 DS
Traits reported for weed-suppressive ability and weed tolerance are plant height, tiller number per plant, panicle number, flag-leaf length, flag-leaf width, leaf area, number of filled grains, number of unfilled grains, and grain yield per plant (Huel and Hucl Reference Huel and Hucl1996; Lanning et al. Reference Lanning, Talbert, Martin, Blake and Bruckner1997; Lemerle et al. Reference Lemerle, Verbe, Cousens and Coombes1996). However, growth duration and flowering date were not associated with weed suppression (Zhao et al. Reference Zhao, Atlin, Bastiaans and Spiertz2006a); thus, these traits were not considered in the 2014 DS. Ability to produce yield and to suppress weeds is moderately heritable and closely associated under weed competition (Zhao et al. Reference Zhao, Atlin, Bastiaans and Spiertz2006a). Results of correlation analysis in this study showed that grain yield was positively and significantly correlated with all the traits and that all the traits were significantly and positively correlated with one another (P<0.001) (Table 6). This also proves the claims of Caton et al. (Reference Caton, Mortimer, Foin, Hill, Gibson and Fischer2000) and Kruepl et al. (Reference Kruepl, Hoad, Davies, Bertholdsson and Paolini2007) that the complexity of traits for weed-competitive ability is controlled by the interaction of several traits rather than a single distinguishable trait.
a Significance codes: * 0.05≥P≥0.01; ** 0.01≥P≥0.001; *** P≤0.001.
The highest correlation (P<0.001) with grain yield was recorded for filled grains (r=0.99), followed by panicle number (r=0.84) and percent spikelet fertility (r=0.83) (Table 6), indicating a very strong correlation of these traits with yield, which is typically expected. Plant height and tiller number had the highest correlation with flag-leaf length (r=0.83 and 0.81, respectively), followed by both flag-leaf width and leaf area (r=0.82 and 0.78, respectively). Leaf area can be related to having a wider and longer flag leaf, which is essential for weed-competitive ability (Huel and Hucl Reference Huel and Hucl1996). Panicle number had the highest correlation with unfilled grains (r=0.88) and filled grains (r=0.86), whereas percent spikelet fertility had the highest correlation with flag-leaf length and flag-leaf width (r=0.91). Percent spikelet fertility is a key trait for yield; thus, correlation analysis suggests that flag-leaf length and flag-leaf width are strong traits for yield in weed-competitive rice. Flag-leaf length and orientation could affect crop competitiveness because of their effects on light interception by the crop canopy (Huel and Hucl Reference Huel and Hucl1996).
Correlation Analysis for Weed-Competitive Ability Traits under Greenhouse Conditions
Results of correlation analysis under greenhouse conditions showed that grain weight was positively and significantly correlated with all the traits and that all the traits were significantly and positively correlated with one another (P<0.001), except for plant height at 7 DAS (P>0.05; r=0.09) (unpublished data). The highest correlation with grain yield was recorded for maturity and yield-related traits: filled grains (r=0.98) and panicle fresh weight (r=0.77). The trait with the lowest correlation with grain weight was plant height at 7 DAS (r=0.16). This implies that plant height at a very early stage (7 DAS) cannot be associated with the yield potential of ILs. Plant height at 7 DAS also has low correlation with panicle fresh weight and number of filled grains (r=0.16), maximum plant height, flag-leaf length, and panicle number (r=0.18).
Correlation Analysis for Seedling Vigor Traits
The traits associated with rapid seedling biomass accumulation are strongly associated with weed suppression and yielding ability under weed competition (Zhao et al. Reference Zhao, Atlin, Bastiaans and Spiertz2006b). Results of correlation analysis in Table 7 showed that vigor index was positively and significantly correlated with all the traits and that all the traits were significantly and positively correlated with one another, except for root dry weight and percent germination (P>0.05; r=0.07). All the rice traits were negatively correlated with weed density, weed fresh weight, and weed dry weight, except for root length, root fresh weight, and root dry weight. Leaf number at 28 DAS, tiller number, leaf dry weight, and total dry weight of ILs did not have any significant correlation with weed parameters. Likewise, leaf number at 7 DAS and leaf chlorophyll content (SPAD reading) were not associated with weed density. Few studies have been conducted on the contribution of root growth and activity to crop–weed competition. In rice, however, Fofana and Rauber (Reference Fofana and Rauber2000) found that competitive varieties tended to have desirable root growth. In this study, root length and root dry weight had no significant correlation with weed fresh weight (r=−0.01) and weed dry weight (r=0.1), respectively.
a Abbreviations: %G, percent germination; 7PH, plant height at 7 d after sowing (DAS); 14PH, plant height at 14 DAS; 21PH, plant height at 21 DAS; 28PH, plant height at 28 DAS; 7LN, no. of leaves at 7 DAS; 14LN, no. of leaves at 14 DAS; 21LN, no. of leaves at 21 DAS; 28LN, no. of leaves at 21 DAS; LCC, leaf chlorophyll content; LFW, leaf fresh weight.; LDW, leaf dry weight; RFW, root fresh weight; RDW, root dry weight; RL, root length; TDW, total dry weight; TFW, total fresh weight; TN, tiller number; VI, vigor index; WD, weed density; WDW, weed dry weight; WFW, weed fresh weight.
b Significance codes: * 0.05≥P≥0.01; ** 0.01≥P≥0.001; *** P≤0.001.
Correlation Analysis for Seed Germination Traits
Results of correlation analysis in Table 8 showed that first-count germination (number of germinated seeds after 48 h) is positively and significantly correlated with second-count germination (number of germinated seeds after 7 d), percent germination, total fresh weight, total dry weight, and vigor index. However, first-count germination has a significant negative correlation with shoot and root length, average fresh weight, and average dry weight. Second-count germination and percent germination have a significant positive correlation with total fresh weight, total dry weight, and vigor index, but a nonsignificant negative correlation with shoot and root length.
a Abbreviations: %G, percent germination; 1st count, first-count germination; 2nd count, second-count germination; AFW, average fresh weight; ADW, average dry weight; Root lt, root length; Shoot lt, shoot length; TFW, total fresh weight; TDW, total dry weight.
b Significance codes: * 0.05≥P≥0.01; ** 0.01≥P≥0.001; *** P≤0.001.
c Vigor index was measured by computing percent germination×seedling dry weight.
Shoot and root length are significantly and positively correlated. These two traits have a nonsignificant positive correlation with total dry weight and nonsignificant negative correlation with vigor index. Finally, average dry weight is not correlated with root length and total fresh weight, whereas average fresh weight is significantly and negatively correlated with total dry weight. Diwan (Reference Diwan2006) reported that first-count germination, rate of germination, and seedling dry weight were significantly correlated with vigor index. In this study, the same result was observed. Vigor index is strongly correlated with first-count germination, second-count germination, percent germination, total fresh weight, and total dry weight. However, vigor index has a strong negative correlation with average fresh weight and average dry weight, and a nonsignificant negative correlation with shoot length and root length.
Regression Analysis for Weed-Competitive Ability Traits under Field Conditions in the 2013 WS
Regression analysis was carried out to determine which traits were significantly affecting the grain yield of the newly developed introgression lines tested under weedy conditions in the 2013 WS. Using R version 1.5, the regression model was selected using all possible regression with leaps and the coefficient of variation (r2) as the criterion. The final regression model was YLD~Max.PH+Max.TN+SPD+PN+PL, which was significant at the 1% level (P<0.001) with an r2 value of 0.4083 (Table 9).
This suggests that 40.83% of the variation in the yield of ILs can be attributed to the regression model. The variables with the highest r2 were tiller number and panicle number, indicating the importance of these variables in rice grain yield and weed-competitive ability of ILs in weedy environments. The fitted regression model showed that only number of panicles was significant at the 1% level using the standard t-test, indicating the importance of this trait to yield (Table 10). The regression coefficients for leaf chlorophyll content (SPAD reading) and panicle length were both negative. This implies that even very low values of these two traits can significantly increase the grain yield of the ILs. On the other hand, plant height, tiller number, and panicle number have positive regression coefficients, which means that these variables directly affect grain yield and are important traits under weed competition. These results coincide with the reports of Lemerle et al. (Reference Lemerle, Verbe, Cousens and Coombes1996) and Lanning et al. (Reference Lanning, Talbert, Martin, Blake and Bruckner1997) that the characteristics that appear to confer competitiveness include crop height, tiller number, and percentage light interception.
Regression Analysis for Weed-Competitive Ability Traits under Field Conditions in the 2014 DS
Regression analysis was performed to discover and identify the traits that were significantly affecting the grain yield of the ILs grown in the 2014 DS under weedy conditions. The final regression model was YLD~Max.PH+PN+FLL+FLW+LA+Filled+Unfilled+SPF, which was significant at the 1% level (P<0.001) with an r2 value of 0.9866 (Table 11). This signifies that 98.66% of the variation in the yield of the ILs can be attributed to the regression model. The variables with the highest r2 were panicle number, leaf area, filled grains, unfilled grains, and percent spikelet fertility, which might be an indication of the importance of these traits in the rice grain yield and weed-competitive ability of ILs in weedy environments. The fitted regression model shows that only the number of filled grains was significant at the 1% level using the standard t-test, indicating the significant contribution of this trait to grain yield (Table 12). The regression coefficients for the variables flag-leaf length, flag-leaf width, number of unfilled grains, and percent spikelet fertility are all negative, but only the number of unfilled grains was significant, which indicates that very low values of this variable can significantly affect the yield of the ILs. On the other hand, plant height, panicle number, leaf area, and number of filled grains are all positive, but only the number of filled grains is significant. This signifies that this variable is directly affecting grain yield and weed-competitive ability of the ILs.
Phenotypic Variations in Traits for Weed-Competitive Ability under Field Conditions
To successfully improve a trait through breeding, variation must exist in the trait among the available germplasm and the trait must be heritable. A significant portion of the observable phenotypic variation expressed among genotypes must be attributed to genotypic differences (Fehr et al. Reference Fehr, Fehr and Jessen1987). There is high variation in all the phenotypic traits observed, which implies high genotypic variation for the investigated traits among the ILs. The traits that showed the highest coefficient of variation (CV) values were number of unfilled grains (79.6%), grain yield per plant (79.3%), number of filled grains (78.2%), panicle number (65.2%), and percent spikelet fertility (56.3%). The traits with CV values lower than 50% were plant height (34.9%), tiller number (43.2%), flag-leaf length (46.3%), flag-leaf width (48.7%), and leaf area (49.8%). The mean grain yield per plant in the 2014 DS was 4.03 g (SD=3.19), with an average of four tillers and three panicles. The yield performance of each IL varied with weed regime. To verify this observation, the ILs were analyzed for the significance of genotypic effect per weedy condition, specifically to determine whether each genotype has a differential yield response because of varying weed conditions. The fixed-effects model was used for the ANOVA with the Satterthwaite denominator. The summary of the ANOVA (Table 13) shows a significant genotypic effect on UPW (F=2.47, P=0.0095), UPWF (F=2.74, P=0.0000), and LLW (F=1.78, P=0.0000) treatments at the 1% level of significance. The summary of ANOVA (Table 14) shows significant genotypic effects on plant height, flag-leaf length, flag-leaf width, and leaf area for all weed regimes. Yield-related traits such as number of tillers, number of panicles, number of filled and unfilled grains, and grain yield were significantly affected by genotype only for UPWF and UPW conditions. In LLW conditions, no significant genotypic effect was found on yield-related traits: number of tillers, number of panicles, number of filled and unfilled grains, and grain yield. This suggests that selection for the weed-competitive genotypes was stronger for upland conditions than for lowland conditions.
a MODEL: YLD ~ Max.PH + Max.TN + SPD + PN + PL.Abbreviations: Max.PH, maximum plant height (cm); Max.TN, maximum number of tillers; PL, panicle length (cm); PN, number of panicles; SPD, Soil-Plant Analysis (SPAD) reading.
a Significance codes: * 0.05 ≥P≥0.01; ** 0.01≥P≥0.001; ***P≤0.001.
Phenotypic Variations in Traits for Weed-Competitive Ability under Greenhouse Conditions
High phenotypic variations were observed under greenhouse conditions, indicating high genotypic variation for the investigated traits among the ILs. The highest CV values were observed for panicle fresh weight (77.1%), number of unfilled grains (72.8%), grain yield per plant (68.9%), number of filled grains (67.7%), and percent germination (52.2%). Traits with CV values lower than 50% were plant height at 21 DAS (21.7%), maximum plant height (23.7%), plant height at 14 DAS (26.3%), flag-leaf length (29.4%), and flag-leaf width (33.0%). The summary of ANOVA (unpublished data) shows significant genotypic effects on all traits measured in the greenhouse experiment at the 1% level of significance. The genotypes (ILs) differed significantly in terms of germination, seedling vigor, canopy traits, and maturity and yield-related traits.
Phenotypic Variations in Seedling Vigor Traits
High variation in all the phenotypic traits was observed in the seedling vigor experiment. This implies high genotypic variation for the investigated traits between the ILs. Traits that showed the highest CV values were vigor index (50.1%), initial germination and percent germination (45.4%), root fresh weight (45.3%), root dry weight (40.5%), and leaf fresh weight (39.7%). Traits with CV values lower than 30% were plant height at 28 DAS (15.0%), number of tillers (22.3%), SPAD reading (23.3%), and number of leaves at 28 DAS (25.8%).
The summary of ANOVA (Table 15) shows significant genotypic effects on all traits measured (P≤0.0001) at the 1% level of significance, except for plant height at 7 DAS (P=0.174). Seedling vigor traits with highly significant variations include plant height at 28 DAS, leaf number at 28 DAS, number of tillers, root fresh weight, and root dry weight. Plant height at 21 DAS was significant only at the 5% level.
a MODEL: YLD ~ Max.PH + PN + FLL + FLW + LA + Filled + Unfilled + SPF.Abbreviations: Filled, number of filled grains; FLL, flag-leaf length (cm); FLW, flag-leaf width (cm); LA, leaf area (cm2); Max.PH, maximum plant height (cm); PN, number of panicles; SPF, percent spikelet fertility; Unfilled, number of unfilled grains.
a Significance codes: * 0.05 ≥P ≥0.01; ** 0.01≥P≥0.001; *** P≤0.001.
a Significance codes: * 0.05≥P≥0.01; ** 0.01≥P≥0.001; *** P≤0.001.
a Significance codes: * 0.05≥P≥0.01; ** 0.01≥P≥0.001; *** P<0.001.
b Abbreviations: LLW, lowland weedy; UPW, upland weedy; UPWF, upland weed-free.
a Significance codes: * 0.05≥P≥0.01; ** 0.01≥P≥0.001; *** P≤0.001.
b Abbreviation: DAS, days after sowing.
Phenotypic Variation in Seed Germination Traits
High phenotypic variation for early germination was observed in the seed germination test, indicating high genotypic variation within the ILs. Traits that showed the highest CV values were relative fresh weight (48.3%), first-count germination (44.1%), relative dry weight (33.3%), and vigor index (31.5%). Traits with CV values lower than 25% were total dry weight (17.8%), shoot length (18.6%), root length (24.5%), total fresh weight (17.8%), and second-count germination (21.9%). The summary of ANOVA shows significant genotypic effects on all traits measured (P≤0.0001) at the 1% level, except for root length (P=0.0834) (Table 16). Shoot length was significant only at the 5% level. This implies that the ILs had comparable shoot and root lengths without any distinct differences. Seed germination and vigor traits that showed highly significant variations include first-count germination, second-count germination, percent germination, total fresh weight, total dry weight, and vigor index. Diwan (Reference Diwan2006) reported similar results in a doubled-haploid population of IR64/‘Azucena’.
a Significance codes: * 0.05≥P≥0.01; ** 0.01≥P≥0.001; *** P≤0.001.
Promising ILs with Superior Traits and Weed-Competitive Ability
Figure 5 shows that, under UPWF and UPW conditions, the top five ILs had higher grain yields than the recurrent parent (WTR1), donor parent (Y134), and checks (PSB Rc82 and Apo). This suggests that these ILs have strong weed-competitive ability and are able to produce comparable yields despite the presence of competing weeds. In LLW conditions, the ILs did not perform very well, having very low grain yield. This implies that selection for weed-competitive rice genotypes is stronger in upland conditions than in lowland conditions. Recipient parent WTR1 had the highest grain yield under LLW conditions, while the donor parent Y134 had the lowest grain yield across all weedy conditions. However, the grain yield of the checks was higher than that of both parents under UPWF and UPW conditions.
Table 17 shows the grain yield of the top five ILs in UPWF, UPW, LLW, and greenhouse conditions and their corresponding yield difference with WTR1 (recipient parent), Y134 (donor parent), and checks PSB Rc82 and Apo. The ILs had the highest grain yield across weed regimes. Check cultivar Apo had higher grain yield than the parents (WTR1 and Y134). G-6-Y9-WU-2 (21.13 g plant−1), G-6-Y7-WL3 (12.43 g plant−1), and G-6-L2-WL-3 (15.28 g plant−1) were the highest-yielding ILs under UPWF, UPW, and LLW conditions, respectively (Table 17).
a Italics signify mean or average value. Asterisks denote significant at the 5% level.
In UPWF conditions, only G-6-Y9-WU-2 had a significant yield difference with WTR1, whereas the yield differences between the top five ILs and Y134 were significant at the 0.5% level. G-6-Y9-WU-2, G-6-Y5-WU-3, and G-6-L6-WU-1 had significant yield differences with PSB Rc82. The yield differences between the top five ILs and check cultivar Apo were not significant.
In UPW conditions, G-6-Y7-WL3, G-6-Y7-WL2, G-6-Y6-WU-3, and G-6-Y3-WL3 had significant yield differences with WTR1, Y134, and PSB Rc82. The yield differences between the top five ILs and check cultivar Apo were not significant.
In LLW conditions, G-6-L2-WL-3, G-6-R6-WL3, G-6-L15-WU1, and G-6-L15-WU-1 had significant yield differences with all the parents and checks. In greenhouse conditions, only Y134 had a significant yield difference with the top four ILs: G-6-L9-WU-2, G-6-RF13-WU-1, G-6-L5-WU-2, and G-6-L6-WU-1.
The maximum plant height of ILs, parents, and checks did not differ across weed regimes. However, G-6-Y9-WU-2, G-6-Y7-WL-3, G-6-RF7-WL-3, and G-6-RF13-WU1 were the tallest in UPWF, UPW, LLW, and greenhouse conditions, respectively.
Leaf area is an important trait for weed competitiveness, as it affects competition for light. Flag-leaf length and width are two important traits for determining leaf area. G-6-L6-WU-1, G-6-Y6-WU-3, G-6-L2-WL-3, and G-6-L5-WU-2 had the largest leaf area under UPWF, UPW, LLW, and greenhouse conditions, respectively.
Compared with the parents and checks, the ILs had a higher number of filled grains and percent spikelet fertility and a lower number of unfilled grains. These traits are important components in attaining high grain yield. This suggests that, despite the presence of weed competition, these top five ILs were the most promising for producing higher yields and better morphological and agronomic traits in each condition.
Promising ILs for Replicated Yield Trials
Across the three conditions (UPWF, UPW, and LLW), 21 lines were selected as the most promising ILs based on yield ability and percent spikelet fertility (Table 18). Among the ILs, G-6-Y5-WU-3, G-6-L6-WU-1, and G-6-Y10-WU-3 were the three highest-yielding lines in UPWF conditions. G-6-L2-WL-3, G-6-RF6-WL-3, G-6-L15-WU-1, G-6-Y16-WL-2, and G-6-L6-WU-3 were the top ILs in LLW conditions, whereas G-6-Y7-WL-3, G-6-Y6-WU-3, G-6-Y3-WL-3, and G-6-Y8-WU-1 were the highest yielding in UPW conditions. In greenhouse conditions, G-6-L9-WU-2, G-6-RF13-WU-1, G-6-L5-WU-2, G-6-L6-WU-1, and G-6-RF12-WU-1 were the most promising.
a Selected based on grain yield and percent spikelet fertility.
In summary, correlation analysis on G-6 ILs grown under weedy conditions showed that grain yield was positively and significantly correlated with all the traits and that all the traits were significantly and positively correlated with one another. This confirms previous reports that the complexity of traits for weed-competitive ability is controlled by the interaction of several traits rather than a single desirable trait. Regression analysis indicated that panicle number, leaf area, number of filled grains, number of unfilled grains, and percent spikelet fertility are the important traits for grain yield of G-6 ILs under weed competition. Selection for weed-competitive genotypes was stronger in upland conditions than in lowland conditions. Significant genotypic effects were observed for plant height, flag-leaf length, flag-leaf width, and leaf area for all weed regimes. However, yield-related traits such as number of tillers, number of panicles, number of filled and unfilled grains, and grain yield were significantly affected by genotype only for UPWF and UPW conditions.
In weedy conditions, the mean yield of the G-6 ILs was significantly higher than the mean yield of WTR1, Y134, and checks PSB Rc82 and Apo. The yield performance of ILs was directly affected by traits for weed-competitive ability; thus, G-6 ILs are better weed competitors and can produce comparable yields despite the presence of weeds.
The height of plants at a very early stage (7 DAS) cannot be clearly associated with the yield ability of ILs and other important agronomic traits. A seedling vigor test showed that all traits are significantly and positively correlated with one another, except for root dry weight and percent germination. All the traits are negatively correlated with weed density, weed fresh weight, and weed dry weight, except for root length, root fresh weight, and root dry weight.
A germination test of G-6 ILs indicated that vigor index is significantly and strongly associated with first-count germination, second-count germination, percent germination, total fresh weight, and total dry weight. However, vigor index had a significant and strong negative correlation with relative fresh and dry weights and a nonsignificant correlation with shoot and root lengths.
After three cycles of selection process across weedy conditions, 21 most-promising introgression fixed lines showing superior traits and weed-competitive ability were identified. The successful development of weed-competitive rice cultivars will provide an effective new strategy to manage weeds with a single hand weeding or herbicide application. This will result in lower labor costs, less herbicide use, prevention of possible development of herbicide-resistant weeds, and increased profit for rice farmers. Furthermore, the use of weed-competitive rice varieties will significantly improve yields in aerobic rice systems and can become an important strategy that complements water-saving techniques for promoting the adoption of upland rice production.
Acknowledgments
We would like to acknowledge the Green Super Rice team of IRRI, Graduate School of the University of the Philippines–Los Baños, and the Chinese Academy of Agricultural Sciences for their collaborative efforts in conducting this study and the Bill & Melinda Gates Foundation for funding this research under ID OPP1130530.