Timelines of rice breeding activities 2006–2012.
Until recently, there was limited research on breeding upland rice varieties. Moreover, there is an increasing expansion of rice production from traditional irrigated production areas to rain‐fed environments in the East African region, where drought problem is a serious challenge. To date, several initiatives aimed at increasing rice production have been made. Of the initiatives, promotion of upland rice production has been the most important in Uganda, but yield penalty due to drought continued to be a major drawback. This article traces progress in the upland rice breeding that started with improvement of late maturing varieties that had nonpreferred cooking qualities. Initially, introduced lines were evaluated and released. These varieties are the ‘New Rice for Africa’ (NERICA) that had been generated from interspecific crosses involving Oryza glaberrima and Oryza sativa. Several studies to understand the mode of gene action and modified pedigree breeding approaches for drought tolerance were conducted and used to develop new rice varieties. Up to 11 improved upland rice varieties were released and deployed in the country from 2002 to 2011 as a result of this initiative.
- drought tolerance
- gene action
- modified pedigree breeding
Rice is an important food crop that is consumed mostly outside its major production areas in Uganda, with over 90% of production marketed to urban areas and major institutions within the country. This aspect makes rice to have a long value chain engaging several players. Rice is cultivated under rain‐fed upland conditions, partly rain‐fed lowland conditions and irrigated conditions in Uganda, taking advantage of diverse ecosystems in Uganda . Since the introduction of rice in 1904, Uganda had production under different agro‐ecological conditions covering rain‐fed upland conditions, partly rain‐fed lowland conditions and irrigated conditions. Various production challenges are faced in these production areas.
New technologies are valuable for use in developing the new rice varieties globally including Africa. However, new tools can be most helpful if the existing varieties and candidate lines are properly characterized and documented. The purpose of this paper is to trace upland rice breeding efforts in Uganda and present for alignment, learning, and application in the new technology in their rice breeding programs with focus on breeding for drought tolerance and other stresses. This is critical considering that there is limited research on development of upland rice varieties suitable for production under mild drought conditions in the East African region and other similar Agro ecologies. There is also increasing expansion of rice production from traditional irrigated production areas to rain‐fed environments, where drought problem is an inherent challenge. Indeed, drought emerged as a critical rice production constraint in East Africa [1, 2], particularly in Uganda , as promotion of upland rice was growing in the country.
Many upland rice varieties, earlier introduced in the country, were late maturing and did not have preferred cooking qualities. Later, more introduced lines were evaluated and released. These varieties had been generated through interspecific crossing involving
|Selection of parental lines|
|Making crosses and generating F2 (genetic studies)|
|Evaluation of F2 and harvesting F3|
|Evaluating F3 under rain‐fed lowland and upland conditions (660)|
|Evaluation F4 upland conditions|
|Evaluation F5 in Advanced Trials|
|Evaluation F6 in NPT (12 lines)|
|Evaluation F6 in NPT (10 lines)|
|Testing for grain traits|
|Organoleptic testing of the new varieties|
2.1. Screening of introductions for drought tolerance
A total of 191 rice introductions from major rice breeding centers were evaluated. Of the 191 materials, 77 were
This experiment was conducted at National Crop Resources Institute (NaCRRI), at Namulonge in central Uganda, at 00°32’ N latitude and 32°53’ E longitudes with altitude of 1150 m above sea level during dry season. The soils of the place are clay loam. The period December to March is characteristically the long dry season, but mean long term annual rainfall is 1270 mm.
In order to assess drought stress during reproductive growth stage, drought stress was imposed by terminating irrigation, when about 50% of the population had reached a point where interauricular distance between the flag leaf and penultimate leaf was zero . It is the period when it is about 10 days before anthesis. It is the time when the penultimate leaves were fully expanded. Rainfall during the trial period was recorded. Irrigation was 14 days later, when 30% of the available water had been lost from the soil at 20‐cm depth. The available soil moisture was taken using the ECHO soil moisture tester (Decagon Devices, Inc., Pullman, Washington, USA). All the grains from each panicle were hand threshed and dried. The filled and unfilled grains were then separated using floatation methods.
2.2. Mode of gene action for drought tolerance
This study investigated the nature of inheritance of drought tolerance in crosses between interspecific and intraspecific rice genotypes using secondary traits. Two separate experiments were conducted, using
The aim of the first experiment was to investigate the inheritance of drought tolerance at reproductive growth stage. Eighteen crosses were generated from two sets of 3 × 3 parents using the North Carolina mating design II (NCD II). All the 18 F2 and the 12 parents were evaluated in a 2 × 15 alpha lattice design with two replicates under a rain‐out shelter and nonstress conditions in the field.
Thirty genotypes comprising 18 F2 progenies from sets, A and B, along with the 12 parents were used in this experiment (Table 2). The 30 entries were established in a rain‐out shelter at National Crops Resources Research Institute (NaCRRI), Namulonge. The rain‐out shelter was constructed using translucent sheets for the roof and wire mesh on the sides of the structure to prevent rain water and to allow free air circulation, respectively. In the rain‐out shelter, standard troughs that are 1m wide, 8m long, and 1.5m deep were filled with soil for fallow field from Namulonge. Four troughs were made and filled with the soil, referred to as strips. The seeds were planted in a 2 × 15 alpha lattice design. Two strips represented a replicate. The 12 parental genotypes were planted in three rows planted across the 1 m width strips, while the F2 populations were planted in six rows. The plant to plant spacing was 15 cm making plant population to be 36 for the parental lines and 72 for the F2 lines.
|Experiment||Crossing set||Genotype no||Breeding line||Type||Parent type|
|1||SET A||18||CT 16334(2)‐CA‐2‐M||Interspecific||Male|
|SET B||2||CT 16346‐CA‐20‐M||Interspecific||Male|
|9||CT 16350‐ CA‐5‐M||Interspecific||Male|
|2||SET C||18||CT 16334 (2)‐CA‐2‐M||Interspecific||Female|
A second set of the 30 entries were planted in the field under optimal conditions. These conditions involved irrigating the field at 20 mm per week, during the period when there was no rain. In both trials, a 2 × 15 alpha lattice design planted in two replicates was used. Two seeds from each generation were drilled at a depth of 3 cm at spacing of 20 × 20 cm in each plot. In order to reduce border effects, 20 cm was left between plots. The 12 parental genotypes were planted in 5‐row and 3‐column plots, while the F2 populations were planted in 5‐row and 6‐column plots. Overall, there were 15 plants per replicate of the parents and 30 plants per replicate of each F2 genotype, thus the total number of plants were 30 and 60 for the parents and F2, respectively. The plants were thinned to one plant per hill. Standard cultural practices including hand planting and hand weeding were followed. The crops were fertilized with 25 kg N ha-1 at 20–25 days after transplanting (DAT) and the same rate at 40–45 DAT to enhance plant vigor.
Drought stress was imposed by terminating irrigation, when about 50% of the populations had attained an interauricular distance between the flag leaf and penultimate leaf of zero, that is the period about 10 days before anthesis [4, 5]. This method of identifying the stage of imposing drought was applied both in the field and in the rain‐out shelter. In general, this is the time when the penultimate leaves were fully expanded. Rainfall during the trial period was recorded.
In the field experiment, irrigation was applied using sprinkler irrigation. The field was irrigated, every three days before imposing drought stress. On the day the irrigation was terminated, the field was irrigated to field capacity in the evening between 5:00 and 6:00 pm, which was resumed 14 days after its termination using sprinkler irrigation. The duration of drought stress was determined by testing the level of soil moisture daily, using the ECHO soil moisture tester (Decagon Devices, Inc Pullman, Washington USA). On the day, when 30% of the available water had been lost from the soil at 20‐cm depth, irrigation was resumed. In the rain‐out shelter, water was applied using hand irrigation cans but water was calculated for each strip at 140 L per week, which is equivalent to 20 mm per week.
The number of filled grains was counted per panicle at grain maturity period. Two panicles from each plant were randomly collected and record of number of filled grains was determined using floatation method described by these authors .
In the second experiment, the magnitude and direction of gene action for drought tolerances at reproductive stage was determined in five populations P1, P2, F1, F2, and F3 generated from a drought tolerant × susceptible cross using generation mean analysis (GMA). They are in set C (Table 2). The materials were planted in the dry season and drought was imposed by terminating at the stage of panicle initiation.
In this experiment, all the five populations generated from crossing; parents P1 and P2, their F1, F2 and F3 genotypes were planted following a randomized complete block design (RCBD) with two replicates. Two seeds from each generation were drilled at a depth of 3 cm at spacing of 20 × 20 cm in each experimental unit (plot) in the field at NaCRRI. The generations P1, P2, and F1 were planted in 5‐row and 3‐column plots, while F2 and F3 were planted in 5‐row and 6‐column plots. Overall, there were 15 plants per replicate of the parents, 30 plants per replicate of each F2 genotype, thus the total numbers of plants were 30 and 60 for the parents, F2, respectively. The cultural practice in experiment 1 was followed, and drought stress was imposed following procedures in experiment 1.
2.2.1. Data analysis
Data was analyzed in three parts, namely analyses of variance, residual maximum likelihood (REML), regression, and generation means. The analysis of variance was performed for different traits associated with drought tolerance in the two sets of populations, A and B, pooling for both stress and nonstress environments. Using REML, the separate sets were analyzed for each trait. The analyses of the variance components of genotypes were further partitioned into variations, due to parents and crosses.
General analyses of variance were performed for filled grains, grains per panicle, leaf area, plant height, tiller number, and panicle number of all hybrids including checks. Genetic analyses for the six parameters of experimental hybrids were then performed in GenStat  as a fixed effects model across two locations  as follows:
Generation mean analysis of the genotypes CT 16334 (2)‐CA‐2‐M crossed with WAB 450‐1‐BL1‐136‐HB was used to determine additive, dominant, and epistatic effects following the model . The various generations did not have equal variances; therefore, weighted inverse of the variances was used in subsequent analysis according to these authors in Ref. . Regression analysis procedures were used to find the best fit model. It is a graphical method used to compare the additive model with additive‐dominance models. Any effect that was not significant at 5% level was excluded from the model. The parameters were fitted using weighted mean squares as described by Ref. .
A scaling test was conducted using linear combinations of various means according to Refs. [9–12] to detect the presence of nonallelic interactions that are known to bias estimates of additive and dominance components in the populations when present. However, in this case, where F3 populations are used instead of backcross populations, the additive effects estimate is for both additive effects and additive × additive interaction effects. Similarly, the dominance effect combined both dominance effects and dominance × dominance interaction effects as a single estimate . This is not a major drawback considering that most breeding work exploits additive effects and dominance effects. Standard errors of generation means were computed by performing nested analysis of variance following methods used in Ref. .
In order to verify the number of genes involved in the transmission of traits associated with drought tolerance, Castle‐Wrights formulae described in Ref.  was used.
2.3. Evaluation of segregating lines
2.3.1. Preliminary evaluation 1
Preliminary yield trials were conducted on station with objective of varietal screening, evaluation, and seed increase. These were F3 selections from the previous experiment. Overall 660 genotypes were selected from the F3 generation based on field performance. All the seed from each of the 660 hills of F3 genotypes were divided into three sets. One set was remnant; a second set was planted in the rain‐fed low‐land environment, while the third set was planted under rain‐fed upland conditions, all on station within Namulonge. The planting was in November, 2010. All the seed from each hill was planted to 5‐m long rows and evaluated.
The second set was planted under rain‐fed low‐land with ample moisture throughout the growth period of the lines. The evaluation focused on maturity period, tillering capacity, presence of foliar diseases, and physical grain characteristics. Lines that had longer maturity period than the variety NERICA‐4, number of reproductive tillers less than 5 per hill, presence of foliar diseases, and grain discoloration were eliminated. Besides, infection by common pathogen namely rice blast, bacterial leaf blight, grain discoloration, and sheath rot were used to eliminate lines. The team that evaluated the materials comprised of scientists, farmers, and rice field workers.
The third set comprising all the 660 lines was planted under rain‐fed upland conditions. Selection was made as previously stated for rain‐fed low‐land conditions. Unlike selection under rain‐fed lowland conditions where a minimum of seven productive tillers was considered acceptable, in this production environment, five productive tillers was considered the minimum.
2.3.2. Preliminary evaluation 2
2.4. Evaluation of promising lines
2.4.1. On‐farm evaluation
On‐farm trials were conducted through participatory and multilocational testing of selected upland varieties. Selection of sites for participatory and multilocational testing considered the following: (i) key representative ecological zones, (ii) participation of stakeholders, and (iii) availability of resources to effectively conduct the exercise. Seed companies were invited and a proposed method of allowing most seed stakeholders to participate in variety evaluation was adhered to. Twenty lines that were tested in 2011B in two locations namely Namulonge and Kibaale. Subsequently, 12 lines were tested in 2012A and finally 8–10 in 2012 B.
In order to identify suitable upland rice varieties, the best rice lines from preliminary trials were submitted to advanced yield trial (AYT). These genotypes were WAB 95 B‐B‐40‐HB (the best performing line among lines received through STRASA and two best lines selected from Upland Regional performance trial (ART3‐11L1P1‐B‐B‐2 and ART8‐L15P14‐1‐2‐1) as well as three genotypes that performed well among 600 new lines developed at Namulonge. The 2011 season II was suitable for selecting high yielding diseases resistant. For instance, WAB95‐B‐B‐40‐HB and WAB788‐16‐3‐2‐1‐HB earlier selected had considerable symptoms of BLS and narrow leaf spot.
2.5. Variety release and status
In the year 2013, six best performing rice varieties were presented for release to the Variety Release Committee in Uganda. Among the traits and Characteristics that was provided as evidence of superiority to the existing rice varieties were, higher yield, preferred grain and cooking qualities, maturity, tolerance to stresses especially drought.
3.1. Screening introductions for drought tolerance
A list of only 30 genotypes, including top 20 and bottom 10 least performing genotypes, in terms of filled grains are presented in Table 3. Among the top 20 genotypes, three namely NERICA 7, CO 39, and VANDANA were reference materials for high drought tolerance at reproductive growth stage. There were nine out of the 20 lines from the CT breeding lines and five from the NERICA generations.
|No.||Genotypes||Filled grains (%)|
|177||WBK 35 (F3)||84.9|
|Bottom 10 genotypes|
|169||IR 57514‐PMI 5‐B‐1‐2||49.4|
|104||ITA 123 (FKR 28)||47.9|
3.2. The mode of gene action for drought tolerance
3.2.1. Gene action
Generalized linear analysis for different traits pooled across sets and sites are presented in Table 4. Results showed that both GCA and SCA effects within sets for filled grains, grains per panicle, leaf area tiller number, and number of panicles per plant were significant (P = 0.001), while only the GCA effects within sets for tiller number were significant (P = 0.001) but not the SCA effects.
|Source of variation||Mean square value|
|d.f||Spikelet fertility||Grains per panicle||Leaf area||Plant height||Tiller number||Panicle number|
|Env × Set||1||14.9***||15.4***||35.4***||1.9||6.3**||13.2***|
|Env × Set/GCAf||4||14.8***||8.4***||14.7***||1.4||4.5***||0.3|
|Env × Set/GCAm||4||5.6***||7.9***||11.9***||3.7**||2.5**||0.1|
|Env × Set/SCA||8||2.8**||7.2***||8.3***||1.8||1.6||0.5|
The male and female mean squares were all significant (P = 0.05) for the filled grains under drought stress (DS) and nondrought stress (NDS) conditions for the A set population (Table 5). There was significant (P < 0.05) mean square for male × female interaction for the filled grains under NDS for set A. In the case of the B crossing set, the male × female interaction mean squares were significant under DS and NDS conditions. In addition, the mean square of male and female were significant under DS but not under NDS. The male, female, and male × female interaction mean squares were all highly significant (P = 0.001) for the total number of grains per panicle under NDS conditions for the A set and significant (P = 0.05) under DS conditions. In the case of B crossing set, male, female, and the male × female interaction, mean squares were significant under NDS conditions, but not the case under DS conditions. The set A had highly significant (P < 0.001) male, female, and male × female mean squares for leaf area under NDS conditions. Mean squares for male and male × female interactions were significant for leaf area under DS conditions, but not the case of female mean square. The results of the B crossing set revealed that male, female, and the male × female interaction mean squares were all highly significant (P < 0.001) under the NDS and DS conditions.
|Mean square values for sets A and B|
|Set A||Set B|
|Source||d.f||Drought stress||Nondrought stress||Drought stress||Nondrought stress|
|Male × female||4||1.56||6.84*||3.29*||0.80**|
|Male × female||4||2.64*||9.40***||0.62||5.83***|
|Male × female||4||4.76***||9.74**||17.74***||5.34***|
|Male × female||4||8.14***||8.14***||15.70***||14.70***|
|Male × female||4||2.17***||13.97||15.70***||1.67|
The male × female interaction mean squares were all highly significant (P = 0.001) for the plant height under DS, and NDS conditions for the A and B populations. There was significant mean square for female effects under NDS, but not under DS conditions for the A and B populations. In the case of B crossing set, male, female, and the male × female interaction mean squares were all highly significant under both NDS and DS conditions. Mean squares for male, female, and male × female interaction mean squares were not significant for the tiller number under NDS conditions for the A populations. There was, however, significant mean square for male and the male × female interactions but not female mean squares under DS conditions. On the other hand, the B crossing set had male, female, and the male × female interaction mean squares, all highly significant, under both DS and only the male under NDS conditions.
3.2.2. Relative contribution of GCA and SCA
General combining ability (GCA) for female and male (GCAf and GCAm) in A populations under DS and NDS conditions are presented in Figure 1. The total GCA for both male and female parents (GCAt) under DS were more than 55% for all five traits except leaf area that had 42%. All the SCA values were less than 50%. The SCA effects of tiller number, under DS, was 14% and the male × female interaction was not significant (Table 5). Similarly, the SCA effects for filled grains were not important under DS conditions because there were lack of significance (Table 5). This finding, therefore implies that the additive effects was more important than nonadditive effects for filled grains, total number of grains per panicle, plant height, and tiller number.
Under NDS, however, filled grains, leaf area, and tiller number had GCAt more than 55%. The total number of grains per panicle and the plant height had nearly equal GCAt, when compared with SCA. This finding implies that the additive effects are more important than nonadditive effects for filled grains, leaf area, and tiller number, while additive and nonadditive effects had nearly equal effects for total number of grains per panicle. However, the lack of significance in the male × female interactions for tiller numbers (Table 5) makes the importance of SCA not valid.
Results of the GCA effects for female and male (GCAf and GCAm) for B populations, under DS and NDS conditions, are shown in Figure 2. The GCA total (GCAt) under drought was more than 55% for filled grains, total number of grains per panicle, and tiller number under DS. The GCAt and SCA for plant height were nearly equal, while a very high SCA value of 68% was found for leaf area. All the SCA values were less than 50%; moreover, the SCA effects for tiller number were not significant (Table 5) and that of total number of grains per panicle were also not significant (Table 5). This finding implies that the additive effects were more important than nonadditive effects for filled grains, total number of grains per panicle, and tiller number, while additive and nonadditive effects had nearly equal effects for plant height. Under NDS conditions, however, filled grains, total number of grains per panicle, leaf area, and tiller number had GCAt more than 55%. The plant height had nearly equal GCAt when compared with the SCA. This finding implies that the additive effects were more important than nonadditive effects for filled grains, total number of grains per panicle, leaf area and tiller number, while additive and nonadditive effects had nearly equal effects for plant height.
3.2.3. General combining ability effects
Table 6 showed the GCA effects for filled grains for interspecific and intraspecific rice. The GCA values for filled grains were the only one presented, because other secondary traits had weak correlation with the filled grains, which is a trait associated with drought tolerance. Positive GCA effect is desirable in breeding for improved drought tolerance. Strong negative values of GCA effects of parents show contribution of GCA towards low filled grains, while high positive values show high filled grains. Since both GCA effects and SCA effects were significant for filled grains, the individual values for both GCA and SCA effects are presented (Tables 6 and 7). Parents CT 16350‐ CA‐5‐M, IRAT 325, and NERICA 9 had positive and significant scores of filled grains under NDS conditions. In the DS environment, CK 73 was highly positive and significant at P = 0.001, while CT 16344‐CA‐9‐M, NERICA 9, and CT 16346‐CA‐20‐M had positive and significant filled grain scores at P = 0.01.
|Nondrought stress||Drought stress|
|CT 16350‐ CA‐5‐M||5.26***||-1.16|
|Nondrought stress||Drought stress|
|NERICA 8 × CT 16334(2)‐CA‐2‐M||-3.46||-12.86***|
|NERICA 13 × CT 16334(2)‐CA‐2‐M||1.47||4.48*|
|IRAT 325 × CT 16334(2)‐CA‐2‐M||0.73||7.21**|
|NERICA 8 × WAB 365‐B‐1H1‐HB||-0.59||8.53***|
|NERICA 13 × WAB 365‐B‐1H1‐HB||-2.01||-4.93*|
|IRAT 325 × WAB 365‐B‐1H1‐HB||1.90||-5.71*|
|NERICA 8 × NERICA 9||3.60||0.92|
|NERICA 13 × NERICA 9||-1.12||3.91|
|IRAT 325 × NERICA 9||-0.46||-1.52|
|Bonanca × CT 16346‐CA‐20‐M||6.15**||-3.19|
|WITA 1 × CT 16346‐CA‐20‐M||-6.39**||-0.88|
|CK 73 × CT 16346‐CA‐20‐M||-2.32||1.86|
|Bonanca × CT 16350‐ CA‐5‐M||1.30||-1.74|
|WITA 2 × CT 16350‐ CA‐5‐M||7.12**||-1.63|
|CK 73 × CT 16350‐ CA‐5‐M||-3.17||2.26|
|Bonanca × CT 16344‐CA‐9‐M||-10.27***||6.02**|
|WITA 2 × CT 16344‐CA‐9‐M||4.24*||-3.67|
|CK 73 × CT 16344‐CA‐9‐M||3.30***||1.12|
3.2.4. Specific combining ability effects
Superior crosses were observed, with positive SCA effects (Table 7). Under nondrought stress conditions, crosses WITA 1 × CT 16350‐ CA‐5‐M, Bonanca × CT 16346‐CA‐20‐M, and CK 73 × CT 16344‐CA‐9‐M had significant filled grain score of 0.01%, 0.01%, and 0.001%, respectively. The cross WITA 1 × CT 16344‐CA‐9‐M had significant filled grain score at 0.05 level of significance. In the drought stress conditions, the cross NERICA 8 × WAB 365‐B‐1H1‐HB were highly significant at P = 0.001, and Bonanca × CT 16344‐CA‐9‐M and IRAT 325 × CT 16334(2)‐CA‐2‐M were positive and significant at 0.01.
3.2.5. Summary of analysis of generation of means
The mean, variance, and mean variance of filled grains for P1, P2, F1, F2, and F3 are shown in Table 8. The F2 populations had the highest variance followed by F3 and F1. Scaling tests for dominance × dominance and additive × additive interactions were nonsignificant for both levels. Dominance main effects were not significant, but additive main effects were significant at P = 0.01. When the mean scores were fitted to an additive model, it fitted with r2 = 0.77 (Figure 3). Mean filled grains score was best linear unbiased estimator (BLUE) of the traits
|Descriptive summary of generations|
|dominance × dominance||-15||3.559||146||-1.184NS|
|additive × additive||3||3.361||176||-0.893NS|
|Gene effects||Expectation estimates||SE||t = (component/SE)||d.f|
The estimate of the number of genes that control filled grains trait based on Castle‐Wrights method was 0.9 ≈ 1 gene. Estimate of the degree of dominance in the F1 and F2 generation based on the  method was -3 ≈ 0 and 0.9 ≈ 1 level of dominance, respectively.
The narrow sense heritability in the generations from the cross between CT 16334 (2)‐CA‐2‐M and WAB 450‐1‐BL1‐136‐HB using regression of F1 on mid‐parents and F2 to F1 based on single seed decent are shown in Figures 4 and 5 respectively. In the F1 to midparent regression, heritability of 60% was realized, but when F2 was regressed onto F1 means, the heritability estimate was 74%.
3.3. Evaluation of segregating lines
3.3.1. Preliminary evaluation of 660 F3
Results of evaluation of two sets of new 660 at NaCRRI are presented in this section. The first set grown under optimum moisture throughout the growth period is presented in Table 9. The selection pressure was 11.4% (75 out of 660 rows selected) for rain‐fed lowland conditions and 9.85% (65 out of 660 rows selected) for rain‐fed upland conditions. Candidate varieties CAIAPO/CT 16324‐CA‐9‐M, WAB 450‐1‐BL1‐136‐B/WAB 450‐B‐136‐HB, CT 16317‐CA‐4‐M/WAB 365‐B‐1H1‐HB, IRAT 325/WAB 450‐B‐136‐HB, and CT 16342‐CA‐25‐M/CK 73 are among the lines. Overall, 84 genotypes were selected for further evaluation.
|Selection under rain‐fed upland conditions||Selection under rain‐fed lowland conditions|
|One row crosses selected at F4||One row crosses selected at F4|
|1||Bonanca × WAB 881‐10‐37‐18‐3‐P1‐HB||1||WAB 56‐104 × CT 16324‐CA‐9‐M|
|2||IRAT 325 × WAB 450‐B‐136‐HB||2||CT 16350‐ CA‐5‐M × WITA 2|
|3||CT 16355‐CA‐15‐M × IRAT 112||3||CT 16355‐CA‐15‐M × IRAT 112|
|4||WAB 365‐B‐1H1‐HB × WAB 450‐1‐BL1‐136‐HB||4||CT 16317‐CA‐4‐M × IRAT 104|
|5||WAB 450‐B‐136‐HB × IRAT 325||5||WAB 450‐B‐136‐HB × IRAT 325|
|6||CT 16344‐CA‐9‐M × WITA 2||6||WAB 365‐B‐1H1‐HB × IRAT 325|
|7||WAB 365‐B‐1H1‐HB × IRAT 325||7||CT 16313‐CA‐4‐M × Caiapo|
|8||WBK 35 (F3) × WAB 450‐1‐BL1‐136‐HB||8||WAB 56‐104 × CT 16313‐CA‐4‐M|
|9||Bonanca × CT 16346‐CA‐20‐M||9||CK 73 × CT 16346‐CA‐20‐M||9|
|10||CT 16342‐CA‐25‐M × IRAT 257|
|11||WAB 450‐B‐136‐HB × WAB 365‐B‐1H1‐HB||1||IRAT 13 × CT 16342‐CA‐25‐M|
|12||CT 16334(2)‐CA‐2‐M × IRAT 325||2||CT 16346‐CA‐20‐M × Bonanca|
|13||CT 16344‐CA‐9‐M × CK 73||3||CT 16342‐CA‐25‐M × CK 73|
|14||CT 16344‐CA‐9‐M × Bonanca||4||WAB 365‐B‐1H1‐HB × IRAT 325|
|15||Bonanca × WAB 450‐I‐B‐P‐38‐HB||15||5||IRAT 112 × WAB 365‐B‐1H1‐HB||10|
|1||Caiapo × CT 16324‐CA‐9‐M||1||Bonanca × WAB 881‐10‐37‐18‐3‐P1‐HB|
|2||CT 16313‐CA‐4‐M × WAB 56‐104||2||WAB 365‐B‐1H1‐HB × IRAT 325||6|
|3||IRAT 325 × WAB 365‐B‐1H1‐HB|
|4||CT 16324‐CA‐9‐M × WAB 56‐104||1||Caiapo × CT 16324‐CA‐9‐M|
|5||WAB 365‐B‐1H1‐HB × IRAT 325||2||WAB 450‐1‐BL1‐136‐HB × WAB 450‐B‐136‐HB|
|6||WAB 450‐B‐136‐HB × IRAT 112||3||CT 16324‐CA‐9‐M × WAB 56‐104||12|
|7||CT 16334(2)‐CA‐2‐M × IRAT 325||14|
|1||IRAT 325 × WAB 450‐B‐136‐HB||5|
|1||WAB 450‐1‐BL1‐136‐HB × WAB 450‐B‐136‐HB|
|2||CT 16317‐CA‐4‐M × WAB 365‐B‐1H1‐HB||1||WAB 450‐B‐136‐HB × WAB 365‐B‐1H1‐HB||6|
|3||IRAT 112 × WAB 365‐B‐1H1‐HB|
|4||CT 16334(2)‐CA‐2‐M × WAB 450‐1‐BL1‐136‐HB||1||WAB 365‐B‐1H1‐HB × WAB 450‐1‐BL1‐136‐HB|
|5||WAB 56‐104 × CT 16313‐CA‐4‐M||2||Bonanca × CT 16346‐CA‐20‐M||14|
|6||CK 73 × CT 16350‐CA‐5‐M|
|7||IRAT 257 × CT 16355‐CA‐15‐M||21||1||CT 16317‐CA‐4‐M × WAB 365‐B‐1H1‐HB||13|
|1||IRAT 112 × WAB 450‐B‐136‐HB|
|2||IRAT 13 × CT 16342‐CA‐25‐M|
|3||CT 16342‐CA‐25‐M × IRAT 13||15|
3.3.2. Evaluation of 84 F4‐F5 lines
Results of evaluation of 84 rain‐fed segregating lines showed that 20 genotypes were resistant to RYMV, blast, and BLB in all the five locations, namely Namulonge, Kigumba, Kibaale, Lira, and Doho (Table 10). Results of yield in Lira are presented in Table 11. The best six genotypes in yield in descending order are P27‐H14 (11,950 kg/ha), P29‐H4 (9750 kg/ha), P36‐H17 (9313 kg/ha), P5‐H1 (9111 kg/ha), P36‐H9 (8688 kg/ha), and P36‐H4 (8417 kg/ha). When yield, pest and disease resistance, plant height, and panicle length was considered a total of nine lines were nominated for National Performance evaluation.
|1||P 22 H13 WAB 450‐1‐BL1‐136‐HB × WAB 450‐B‐136‐HB||v||0||0|
|2||P 36 H1 WAB 365‐B‐1H1‐HB × WAB 450‐1‐BL1‐136‐HB||0||0||0|
|3||16‐16 CT 16344‐CA‐9‐M × Bonanc||0||0||0|
|4||13‐13 CT 16344‐CA‐9‐M × CK 73||0||0||0|
|6||P 25 H1 CT 16346‐CA‐20‐M × Bonanca||0||0||0|
|7||P 8 H2 Caiapo × CT 16324‐CA‐9‐M||0||0||0|
|11||P 24 H8 IRAT 13 × CT 16342‐CA‐25‐M||0||0||0|
|12||P 1 H14 Bonanca × WAB 881‐10‐37‐18‐3‐P1‐HB||0||0||0|
|13||P 4 H6 CT 16350‐CA‐5‐M × WITA 2||0||0||0|
|14||P 29 H1 CT 16342‐CA‐25‐M × CK 73||0||0||0|
|15||P 23 H1 CT 16346‐CA‐20‐M × WITA 2||0||0||0|
|16||P 45 H15 WAB 365‐B‐1H1‐HB × IRAT 325||0||0||0|
|15||P 24 H9 IRAT 13 × CT 16342‐CA‐25‐M||0||0||0|
|18||P 27 H10 CT 16317‐CA‐4‐M × WAB 365‐B‐1H1‐HB||0||0||0|
|19||P 5 H2 IRAT 325 × WAB 450‐B‐136‐HB||0||0||0|
|20||P 29 H4 CT 16342‐CA‐25‐M × CK 73||0||0||0|
|Rank||Acc no.||Yield||Rank||Acc no.||Yield||Rank||Acc no.||Yield||Rank||Acc no.||Yield|
3.4. Evaluation of promising lines
Results of evaluation of nine selected lines along with two earlier selected lines and a local check is presented in Table 12. Six lines were selected and presented for release to the National Variety Release Committee of Uganda
|Genotype||Arua||NaCRRI||Masindi||Soroti||Kibaale||Kanungu||Mean yield||Mean rank||Yield under optimal condition|
Seven genotypes namely 1. P5H2 (IRAT 325/WAB 450‐B‐136‐HB‐F6), 2. P29H4 (CT 16342‐CA‐25‐M/CK 73‐F6), 3. P8H2 (Caiapo/CT 16324‐CA‐9‐M‐F6), 4. ART3‐11L1P1‐B‐B‐2, ([WAB56‐104/(WAB56‐104/CG14)]/Moroberekan), 5. P27H1 (CT 16317‐CA‐4‐M/WAB 365‐B‐1H1‐HB‐F6), 6. WAB 95 B‐B‐40‐HB (ITA257/(IDSA6/ROK16)), and 7. P24H9 (IRAT 13/CT 16342‐CA‐25‐M‐F6), higher yields that NERICA‐4 the local check. Under optimum conditions, six genotypes had higher yield than NERICA‐4.
3.5. Varietal release and status of release
Breeding background, characteristics, and selected agronomic information on six varieties were presented to the variety release committee. These varieties were released based on important characteristics detailed in Table 13. The names proposed and accepted by the variety released committee were NamChe‐1. NamChe‐2, NamChe‐3, NamChe‐4, NamChe‐5, and NamChe‐6.
|Variety name||NamChe 1||NamChe 2||NamChe 3||NamChe 4||NamChe 5||NamChe 6|
|Year of release||2013||2013||2013||2013||2013||2013|
|NamChe 1||NamChe 2||NamChe 3||NamChe 4||NamChe 5||NamChe 6|
|ITA257/(IDSA6/ROK16)||Caiapo/CT 16324-CA-9-M-F6||CT 16342‐CA‐25‐M/CK 73||[WAB56‐104/(WAB56‐104/CG14)]/Moroberekan||CT 16317‐CA‐4‐M/WAB 365‐B‐1H1‐HB||IRAT 325/WAB 450‐B‐136‐HB–|
|AfricaRice, Senegal||NaCRRI, Uganda||NaCRRI, Uganda||AfricaRice, Ibadan||NaCRRI, Uganda||NaCRRI, Uganda|
4.1. Screening introductions for drought tolerance
Results that nine out of 20 best lines were from the CT breeding work imply there was adequate variability in the selected set for selection of suitable genotypes. There were five out of 20 genotypes also indentified as drought tolerant. Also, result that three reference genotypes namely NERICA 7, CO 39, and VANDANA were suitable for identified as drought tolerant in this study implies that the method of screening was acceptable.
4.2. Genetic analysis for filled grains and other agronomic traits
The analysis of F2 crosses revealed the various components of gene action controlling various drought tolerance traits in rice. Both male GCA and female GCA effects were significant for filled grains under both DS and NDS, for the A populations and the B populations under DS. The finding that the additive effects were more important than nonadditive effects for total number of grains filled, grains in set A under NDS and B under NDS and DS, implies that additive effects control the traits in different populations and the nonadditive effects varied with populations under study. This result is contrary to the finding by Mohapatra and Mohanty  that filled grains was predominantly controlled by nonadditive gene effects under drought stress. However, in the study by Mohapatra and Mohanty , the populations were generated by crossing
Findings of this study that nonadditive effects for total number of grains per panicle was important in both set A and set B under NDS conditions, implies that breeding methods that involve late selection could improve drought tolerance under NDS using number of grains per panicle trait. The use of yield components including grains per panicle has been demonstrated to be effective in improving yield under drought stress by selecting under NDS conditions . The differences between the responses under DS and NDS conditions for total number of grains per panicle could be due to fewer loci within the set B that could segregate for the trait than in A. Set B comprised of lines with more susceptibility to drought stress than those in set A.
Additive effects for number of grains per panicle were important in all the population in set A and set B, under DS and NDS. This implies that breeding methods that involve selection in the early generations especially, single seed decent, pedigree selection, and modified bulk methods could improve drought tolerance through selection of number of grains per panicle. In another study involving
Findings of this study that nonadditive effects for leaf area were more important than additive effects in both set A and set B under DS conditions, suggests that late selection could improve drought tolerance. In addition, the findings that additive effects were more important than nonadditive effects for the populations in sets A and B under NDS implied that selection methods that involve early selection could be employed under NDS. In the populations in sets A and B, interspecific rice genotypes generated from
Results of this study that additive and nonadditive effects for plant height were nearly equal with contribution for total GCA, varying between 45 and 55% for both set A and set B under DS and NDS conditions, implied that that breeding methods that involve both early and late selection could be employed in the improvement of drought tolerance using this trait. Modified bulk method of selection method could be appropriate. In another study involving
The current study found that additive effects were the more important in the transmission of drought tolerance using tiller number as evidenced by the lack of significance for male × female interaction effects for tiller number. This finding is contrary to the work reported by other scientists that nonadditive effects were more important under drought stress conditions [22, 23]. In another study, however, expression of tiller number, under both NDS and DS situations, was found to involve nonallelic gene interactions .
Overall, in situations where nonadditive effects are more important, selection should be delayed until later generations. In these types of populations, repeated crossing in the segregating generations may be useful to pool all the desirable genes in one genotypes according to Ref. . The modified bulk method is another useful method of improvement. However, when additive affects are more important, then a modified pedigree method that involves bulking germplasm before evaluation is appropriate. However, when both additive and nonadditive effects are important, two options can be taken depending on the objective of the breeding and the relative importance of the additive or nonadditive effects. In case, if the objective is to develop hybrid rice, as it is planned in Uganda, then pure line selection should be employed. In this approach, additive effects will be extracted because rice is autogamous . In a situation, where both additive and nonadditive gene action are to be exploited, a modified bulk breeding method would hasten the rate of genetic improvement. Similar exploitation of both additive and nonadditive gene action has been conducted in the improvement of cold tolerance  and sodicity tolerance in rice .
4.3. Combining abilities filled grains under drought stress
Generally, there was no clear distinction in combining ability between
Specific combining ability analysis revealed that crosses WITA 1 × CT 16350‐CA‐5‐M, Bonanca × CT 16346‐CA‐20‐M, and CK 73 × CT 16344‐CA‐9‐M were best under NDS condition. The cross CK 73 × CT 16344‐CA‐9‐M had both parents as good combiners indicating additive × additive type of gene action. It is expected that these crosses could provide transgressive segregants that could be selected using pedigree methods . The others crosses had mixed combiners, therefore additive and nonadditive gene action could be the major contributors. In such crosses, bulk breeding methods could exploit both gene actions.
4.4. Generation means for filled grains under drought stress
There were significant differences among generations for filled grains indicating the presence of sufficient genetic variability. Variability for various traits of rice has been reported [29–32]. The scaling test showed that additive genetic effects but not dominance and epistatic genetic effects were important in the inheritance of filled grains. Fitting means of filled grains on the additive model showed that additive effects accounted for 77% of the genetic variation. In addition, the finding that dominance level was 0 in the F1 population showed that there were no dominance effects.
The generation means analysis confirmed that additive effects were significant in the transmission of filled grains in the populations generated. This study had no inconsistencies in detecting that additive effects were the most important genetic factor in the population under study. In addition, results where narrow sense heritability was high indicated that a high proportion of genetic components of variance can be fixed in segregating generations. Since the selection was conducted under drought stress, it is appropriate that selection for improved drought stress is conducted as early as F2 in the study location. According to Ref. , it is appropriate that selection for improved drought stress is conducted using heritability estimates for target traits. There is limited information on the inheritance of filled grains trait under drought stress. However, various reports indicated that additive effects were the main components that controlled the transmission of this trait under high temperature [33, 34]. A single gene pair was estimated to control filled grains under drought stress. A single gene was found to be responsible for the transmission of filled grains under high temperatures [33, 34].
4.5. Evaluation of segregating lines
Results of evaluation of two sets of new 660 genotypes showed that CT lines namely CAIAPO, CT 16324-CA-9- CT 16317-CA-4-M and CT 16342-CA-25-M had the highest number of parents that could improve the landraces. These lines were developed for drought tolerance through CIAT Colombia Breeding program.
Results of evaluation of 84 rain-fed genotypes that P27-H14 P29-H4 P36-H17, P5-H1 P36-H9 and P36-H4 were the preferred genotypes based on resistance to diseases and yield concurs with other reports (3, 23) that rice varieties with tropical Japonica have higher resistance to RYMV and other diseases.
4.6. Evaluation of promising lines
Results of evaluation of nine selected lines along with two earlier selected lines and a local check is presented in Table 12. Although seven varieties were more had higher yields than NERICA‐4, only six were presented for release when information on milling and cooking qualities were considered. These genotypes were:
4.7. Varietal release and status of release
Breeding background, characteristics, and selected agronomic information on six varieties were presented to the variety release committee. These varieties were released based on important characteristics summarized in Table 13. Information from genetic studies during F2 generation guided selection of promising lines from F2 through F6. Subsequently, promising varieties were nominated for National Performance Trials and eventually released. Four new varieties were released namely, NM7‐8‐2‐B‐P‐11‐6 generated from CAIAPO/CT 16324‐CA‐9‐M cross, NM7‐29‐4‐B‐P‐80‐8 (CT 16342‐CA‐25‐M/CK 73), NM7‐5‐2‐ B‐P‐79‐7 (IRAT 325/WAB 450‐B‐136‐HB), NM7‐27‐1‐ B‐P‐77‐6 (CT 16317‐CA‐4‐M/WAB 365‐B‐1H1‐HB), and NM7‐5‐2‐ B‐P‐79‐7 (IRAT 325/WAB 450‐B‐136‐HB). These varieties were assigned release names, where WAB95‐B‐B‐40‐HB was named NamChe‐1 at the release in Uganda and ARICA 5 by the AfricaRice Breeding Task Force. ARICA acronym means advanced rice for Africa, implying that the harmonized names are to be used by all parties involved. Another variety bred by AfricaRice is NamChe‐1 (ARICA‐5) with designation ART3‐11L1P1‐B‐B‐2. Of the six varieties released, four were bred from Uganda with support from Alliance for Green Revolution in Africa (AGRA) and the other two were developed by AfricaRice through the AfricaWide Rice Breeding Task Force with support from Stress-tolerant rice for poor farmers in Africa and South Asia. These were NamChe‐2 (NM7‐8‐2‐B‐P‐11‐6), NamChe 3 (NM7‐29‐4‐B‐P‐80‐8), NamChe 5 (NM7‐27‐1‐ B‐P‐77‐6), and NamChe 6 (NM7‐5‐2‐ B‐P‐79‐7). The acronym NamChe means
This research found that there was adequate variability in the rice population studied for secondary traits for drought tolerance namely, leaf roll and filled grains. However, the filled grains were found to be more informative and therefore recommended for further studies. Of the three rice groups
The genetic studies for drought provided information on the gene action for drought tolerance at reproductive stage of crosses between interspecific and
Follow up of their performance in countries in the region shows that NamChe‐3 (NM7‐29‐4‐B‐P‐80‐8) and NamChe‐2 (NM7‐8‐2‐B‐P‐11‐6) could be mega variety and a major source of disease resistance. In 2015, over 20,000 ha was under production based on figures of direct seed sale by different producers. This is a success story demonstrating the benefit of collaboration and rigorous breeding in the development of locally adapted rice varieties.
This work was supported by the Government of Uganda, AGRA‐PASS program, and AfricaRice Breeding Task Force program of AfricaRice and University of Kwa‐Zulu Natal, in South Africa.
Lamo J, Imanywoha J, Bigirwa G, Walusimbi M, Kyetere D, Kikafunda J, and Kalule T. 2010. First NERICA rice released in Uganda tops farmer's rankings. International Rice Research Notes, 35:2.
ECARRN. 2005. Five Year Priority setting for the East, Central and Southern Africa Rice Research Network( EACRRN) 140 p. Mikocheni Agricultural Research Institute, Mwenge Coca Cola Road P O Box 6226, Dae salaam, Tanzania.
Lamo J, Asea G, Otim M, Stella Adur, Serrumaga J, Onaga G, Tsuboi, Baboucarr M, Mande S, Moussa S, Cesear M, Okanya S, Ekebu J, and Ochen S. 2014. Release of 6 upland rice varieties in Uganda in 2014. 41 pages. NARO Annual Report 2014.
Ji X, Ende WE, Schroeven L, Clerens S, Geuten K, Cheng S, and Bennett J. 2007. The rice genome encodes two vacoular invertases with fructan exohydrolase activity but lacks the related fructan biosyntheis genes of the Pooidaee. The New Phytologist, 173:50–163.
Itoh JI, Nonomura KI, Ikeda K, Yamaki S, Inukai Y, Yamagishi H, Kitano H, and Nagato Y. 2005. Rice plant development: from zygote to Spikelet Plant Cell Physiol 2005;46:23–47. doi: 10.1093/pcp/pci501
IRRI. 2002. Field evaluation hand book of the International Rice Research Institute. IRRI Publications, Manila, Philippines.
Payne RW, Murray DA, Harding SA, Baird DB, and Soutar DM. 2007. GenStat for Windows, 10th Ed. Introduction, VSN International, Hemel Hempstead.
Hallauer AR and Miranda F. 1988. Quantitative genetics in maize breeding, Iowa State University Press, Ames, Iowa.
Hayman BI. 1958. The separation of epistatic from additive and dominance variation in generation means. Heredity, 12:371–390.
Nigam SN, Upadhyaya S, Chandra S, Nageswara Rao RC, Wright GC, and Reddy AGS. 2001. Gene effects for specific leaf area and harvest index in three crosses of groundnuts, Arachis hypogea. Annuals of Applied Biology. 139(3):301–306. DOI: 10.1111/j.1744-7348.2001.tb00143.x.
Rowe KE and Alexander WL. 1980. Computations for estimating the genetic parameters in joint‐scaling test. Crop Science, 20:109–110.
Jinks JL and Jones RM. 1958. Estimation of components of heterosis. Genetics, 43:223–234.
Nunir M, Chowdhry MA, and Ahsan M. 2007. Generation mean studies in bread wheat under drought condition. International Journal of Agriculture and Biology, 9:282–286.
Sharma JR. 1995. Statistical and biometrical techniques in plant breeding, New Age International Limited, New Delhi, 432 pp.
Peter FC and Frey KJ. 1966. Genotypic correlations, dominance and heritability of quantitative characters in oats. Crop Science, 6:259–262.
Mohapatra K and Mohanty HK. 1985. Inheritance of some quantitative characters including heterosis in rice by combining ability analysis, pp. 579–591, In.Swaminathan MS, ed. Rice Genetics. Proceedings of International Rice Genetics Symposium, Island Publishing House, Manila, Philippines.
Fujii M, Andoh C, and Ishihara S. 2005. Drought resistance of NERICA (New Rice for Africa) compared with Oryza sativaL. and millet evaluated by stomatal conductance and soil content. Proceeding of the fourth International Crop Science Conference.
Atlin G. 2003. Improving drought tolerance by selecting for yield, pp. 14–22, In.Fischer KS, Lafitte R, Fukai S, Atlin Q, and Hardy B, eds. Breeding rice for drought‐prone environments. IRRI Publication, Los Bańos, Philippines.
Ahmad L, Zakri AH, Jalani BS, and Omar D. 1985. Detection of additive and nonadditive variation in rice pp. 555–564, In.Swaminathan MS, ed. Rice Genetics. Proceedings of International Rice Genetics Symposium, Island Publishing House, Manila, Philippines.
Narayanan KK and Sree Rangasamy SR. 1990. Genetic analysis for salt tolerance in rice, pp. 167–174, In.Lampe KJ, ed. Rice Genetics. Proceedings of Second International Rice Genetics Symposium. Island Publishing House Manila, Philippines, International Rice Research Institute (IRRI), Manila, Philippines.
Jones MP, Mande S, Aluko K (1997). Diversity and potential of Oryza glaberrimaSteud in upland rice breeding. Breed Sci. 47:395–398.
Sarker U, Biswas PS, Prasad B, and Mian MAK. 2002. Heterosis and genetic analysis in rice hybrids. Pakistan Journal of Biological Sciences, 5:1–5.
Efisue A, Tongoona P, Derera J, Langyintuo A, Laing M, and Ubi B. 2008. Farmers perceptions on rice varieties in Sikasso region of Mali and their implications to rice breeding. Journal of Agronomy and Crop Science, 4:212–218.
Manonmani S and Fazlullah Kahn AK. 2003. Studies on combining ability and heterosis in rice. Madras Agricultural Journal, 90:228–231.
Chandraratna MF. 1964. Genetics and breeding of rice p. 389 Tropical Science Series, Longman, London, U.K.
Acharya S. 1987. Genetic parameters and their implications in breeding cold tolerance varieties of rice ( Oryza sativaL.). Crop Improvement, 14:100–103.
Geetha S, Shanthi P, Jebaraj S, and Mohammed SEN. 2006. Gene action for sodicity tolerance in rice. Indian Journal of Crop Science, 1:201–202.
Chakraborty S, Hazarika M, and Hazarika G. 1994. Combining ability analysis in rice. Oryza31:281–283.
Singh BN and Mackill DJ. 1990. Genetics of leaf rolling under drought stress. Rice Genetics, Proceedings of the Second International Rice Genetics Symposium II:159–166.
Garrity DP and O’Toole JC. 1994. Screening rice for drought resistance at the reproductive phase. Field Crops Research, 39:99–110.
Fukai S and M Cooper. 1995. Development of drought‐resistant cultivars using physio‐morphological traits in rice. Field Crops Research, 40:67–86.
Lafitte R, Blum A, and Courtois G. 2003. Secondary traits to help identify drought‐tolerant genotypes, pp. 37–48, In. Fischer KS, Lafitte R, Fukai S, Atlin Q, and Hardy B, eds. Breeding rice for drought‐prone environments, IRRI, Los Bańos, Philippines.
Sun ZX, Xiong ZM, Min SK, and M SH. 1989. Identification of the temperature‐sensitive male sterile rice. Chinese Journal of Rice Science, 3:49–55.
Ali J, Siddiq EA, Zaman FU, Abraham MJ, and Ahmed IM. 1995. Identification and characterization of temperature sensitive genic male sterile sources in rice ( Oryza sativaL.). Indian Journal of Genetics, 55:243–259.