Journal of Plant Pathology Research

 Editor-in-chief

  Dr. Petersen Dan
  University of Cincinnati,   USA

Research Article | Volume 5 | Issue 1 | DOI: 10.36959/394/631 Open Access

Variability and Correlation Study on Lowland Green Super Rice in Western Ethiopia

Mulugeta Bitew, Desta Abebe, Desalegn Wondiferaw and Abebaw Desie

  • Mulugeta Bitew 1*
  • Desta Abebe 1
  • Desalegn Wondiferaw 1
  • Abebaw Desie 1
  • Department of Plant Breeding, Debere Markos Agricultural Research Center Crop Research Process, Pawe Agricultural Research Centre, Ethiopia

Bitew M, Abebe D, Wondiferaw D, et al. (2023) Variability and Correlation Study on Lowland Green Super Rice in Western Ethiopia. J Plant Pathol Res 5(1):99-103

Accepted: March 11, 2023 | Published Online: March 13, 2023

Variability and Correlation Study on Lowland Green Super Rice in Western Ethiopia

Abstract


The success of plant breeding activities entirely depends on the existence of genetic variability with respect to desired traits and selection skills of plant breeders. A field experiment was conducted to evaluate the extent of variability, heritability genetic advance, and character association among yield and yield-related traits in lowland GSR. The experiment was conducted using 20 low-land green supper rice genotypes in the Northwestern part of Ethiopia at Pawe agricultural research center for three consecutive years from 2016-2018 main cropping seasons. A randomized complete block design with three replications was used. The results of the analysis of variance for yield and yield-related characters revealed that highly significant mean square due to genotypes observed for all traits studied, Higher phenotypic co-efficient of variation, compared to genotypic co-efficient of variation was recorded for all the traits studied indicating the variation is not only genetic but also influence of growing environment. The highest PCV was recorded for the trait plant height (32.6) followed by panicle length (23.3). Broad sense heritability analysis showed that days to heading (56.3), plant height (49.8), Days to maturity (42.4%), and panicle length (42.3%) were high. High heritability with the high genetic advance in percent of mean was recorded by plant height (49.8 & 33.54) followed by Panicle length (42.3 & 36.55) which confirms their additive gene action and direct selection of these traits depending on phenotypic expression would be effective for yield improvement. The positive and highly significant phenotypic correlation was observed between days to heading and days to maturity (r = 9385***) plant height was positively correlated with thousand-grain weight (0.50436*) and Panicle length (0.48312*). Grain yield showed negative correlation between Days to maturity (-0.044662*), but weak positive correlation was observed for days to heading (r = 0.33472), plant height (r = 0.18729) and panicle length (r = 0.08771). Therefore, the traits recorded high heritability and higher genetic advance; and the traits showed positive correlation and could be considered as suitable traits for selection and improvement to achieve high-yielding green super rice varieties in rain-fed environments.

Keywords


GSR, Grain yield, Heritability, Variability, Yield related traits

Introduction


In Ethiopia, rice production has started three decades ago and the country has proved to have reasonable potential to grow different rice types in rain-fed lowland, upland, and irrigated ecosystems. The estimated potential of Ethiopia for upland rice production is 30 million hectares [1]. It is among the important cereal crops grown in different parts of Ethiopia as a food crop next to teff, maize, wheat and sorghum. A wide production of rice in the country is believed to make a great contribution to food security. Considering the importance and potential of the crop, it has been recognized by the Government as "the new millennium crop of Ethiopia" to attain food security [2]. However, the productivity of rice in Ethiopia is 3.04 t ha-1 which was much lower than the world's average 4.67 t ha-1 [3] due to severe disease and insect pest occurrence, inspired yield increment and over application of fertilizer, water stress and terminal drought occurrence and soil fertility depilation due to continuous cropping are the major challenges that obstruct rice production in the country. As a food security crop in the country, rice production must meet increasing requirements, in both quantity and quality, and should be in harmony with the environment thus ensuring a proper level of sustainability. This requires a gradual reduction in application of pesticides, fertilizers, and water while still achieving continuous yield increase and quality improvement [4]. Recently, green super rice (GSR), defined as rice varieties that can produce high and stable yields under fewer inputs (water, nutrients, and pesticides) and adverse conditions [5].

Global rice production requirements in the future need to come from increased productivity per se of improved rice varieties and require major shifts in breeding goals and objectives [6]. To achieve sustainable rice production in environmental and economic terms, GSR varieties have numerous properties, high efficiency of fertilizer use, resistance to drought and other abiotic stresses, good grain quality, and high yield potential. Development of high yielding varieties requires the existence of genetic variability and genetic variability of agronomic traits is the key component of breeding programs for broadening the gene pool [7,8]. Therefore, success on plant breeding activities entirely depends on the existence of genetic variability with respect to desired traits and selection skill of plant breeder [9]. Different genetic variability parameters, namely, Genotypic Coefficient of Variability (GCV), Phenotypic Coefficient of Variability (PCV), heritability and genetic advance for yield attributing traits is a major concern for any plant breeder for crop improvement programs [10]. Similarly, information on correlation coefficients between grain yield and its component characters is essential since grain yield in rice is a complex character and is highly influenced by several component characters [11]. While selecting the suitable plant type, correlation studies would provide reliable information in nature, extent, and the direction of the selection, especially when the breeder needs to combine high yield potentials with desirable agronomic traits and grain quality characters. It is important that breeders understand the amount of variation, correlation and inheritance of important agronomic traits. Therefore, the current study was conducted to evaluate the extent of variability, heritability genetic advance and character association among yield and yield related traits on lowland GSR.

Materials and Methods


Experimental sites

The experiment was conducted at Pawe Agricultural Research center, which is located in Metekel zone of the Benishangul Gumuze Regional State in the lowlands of the Northwestern part of the country. Its geographical location is between 11°15'and 11°23" North latitude and 36°30' East longitudes at an altitude of 1120 m.a.s. The mean annual temperature ranges 32 ℃ and 16 ℃. Pawe is characterized as a lowland area with very high rainfall (unique agro ecology) and has unimodal rainfall pattern with extended rainy season from May to October, having mean annual rain fall of 1659 mm per year. Its average humidity is 54.7%, with Wind speed 53.6 km/day with sunshine rate of 7.73 hours.

Experimental materials, design and procedures

The experiment was conducted using 20 low land green supper rice genotypes in the Northwestern part of Ethiopia at Pawe agricultural research center for three consecutive years from 2016-2018 main cropping seasons. Randomized complete block design with three replications was used. Each experimental plot had a total area of 6 m2 (1.2 m × 5 m) and six rows at 0.20 m interval. There was a 0.5 m distance between two consecutive plots within a replication. Seeds were sown in rows with manual drilling at a rate of 60 kg % ha. Fertilizer was applied at a rate of 100 kg DAP and 100 kg Urea per ha. Recommended amount of DAP was applied during planting while urea was applied in three splits at planting, tillering and at panicle initiation stages.

Data collection and analysis

Data was collected on number of fertile tillers per plant, plant height, panicle length, number of filled grains per panicle and number of unfilled grains per panicle were collected on plant base from five randomly selected plants, and days to heading and days to maturity on plot base; harvest index, thousand-grain weight, biological yield and grain yield were collected on plot base from four harvestable rows. Analysis of variance (ANOVA) was carried out on the data to assess the genotypic effects and their interaction using general linear model (GLM) procedure for randomized complete blocks design (RCBD) using SAS (2004) version 9.0. Mean comparisons among treatment means were conducted by the least significant difference (LSD) test at 5% levels of significance.

Result and Discussion


Analysis of variance

The results on analysis of variance for yield and yield related characters are presented in Table 1. Highly significant mean square due to genotypes observed for traits like days to heading, days to 85% maturity, plant height and panicle length at (p ≤ 0.01) level and significant difference (p ≤ 0.05) among genotypes were observed for grain yield and thousand grain weight, indicating the existence of sufficient variability among genotypes for yield and yield related traits studied in the present findings, therefore; there is higher opportunities for effective selection of green super rice for the study area. Significant genetic variation among rice genotypes for yield and yield attributing traits were reported by Bandhi, et al. [12]. Berhanu, et al. [13] reported that Significant genetic variation among genotypes for various characters suggested that the genotypes were genetically diverse and could be a good opportunity for breeders to select genotypes for trait of interest. The combined analysis of variance across years showed highly significant difference (p ≤ 0.01) for days to heading, days to maturity, plant height, and panicle length; and significant difference (p ≤ 0.05) was observed for grain yield and thousand seed weight among the tested varieties. The effect of year was also highly significant on all measured traits (p ≤ 0.01) except plant height, panicle length and thousand seed weight. The interaction of variety and year revealed a significant difference (p ≤ 0.05) for plant height, grain yield and thousand seed weight, whereas days to heading, days to maturity panicle length number of filled grains per panicle, didn't show significant difference.

Estimation of variability parameters

The estimates of genetic parameters of quantitative trait in Green super rice genotypes is presented in Table 1. The result revealed a wide range of variation among 20 low land green super rice for 6 quantitative traits studied. The phenotypic variance of all traits was higher than genotypic variance. Similarly, higher phenotypic co-efficient of variation, compared to genotypic co-efficient of variation were recorded for all the traits studied in the present investigation except numbers of effective tillers per plant, indicating the variation is not only genetic but also influence of growing environment. The GCV and PCV values were categorized as low (< 10%), moderate (10 to 20%) and high (> 20) as indicated by Deshmukh, et al. [14]. The difference between PCV and GCV values were as low as (1.3) for days to maturity and as higher panicle length (17.2). The lower difference may be the lower sensitivity of the trait days to maturity (3.8 & 2.5). And days to heading (5.8 & 4.4) to the environment which indicates that these traits are less affected by the environment for expression of the trait. Sravan, et al. [15] also reported the lowest difference between PCV and GCV value for days to heading and days to maturity; Therefore, selection based on phenotypic performance for days to heading and days to maturity will be effective to bring significant genetic improvement. The highest PCV was recorded for the trait plant height (32.6) followed by panicle length (23.3). This result was in agreement with the study of Anis, et al. [16] for plant height and panicle length. Moderate PVC was recorded for the trait grain yield and thousand seed weight. In contrast the lowest PVC were recorded for the trait days to maturity (3.8) and days to heading (5.8). All traits recorded the lowest GCV. This may be due to the narrow genetic base of the trait. The lower GCV and PVC were recorded for days to heading and days to maturity. The greater difference between GCV and PCV was detected in panicle length, plant height, grain yield and thousand grain weight. The result indicates that the influence of environmental factors affects the expression of these traits than genetic factors and can be improved by growing within a proper environment. Generally, high coefficient of variability shows scope of selection in favors of traits of interest and low coefficient of variability indicates the need for creation of variability and selection Adhikari, et al. [9].

Estimates of heritability and genetic advance as percent of mean

Knowledge of heritability assists plant breeders in making selection more effective and simpler and together with heritability, genetic advance offers more attainable gain at specific selection intensity Iqbal, et al. [17]. In the present study, genotypic variances were higher in magnitude than corresponding environmental variances for all the studied traits except grain yield, resulting in high heritability for most of the traits.

Among the studied traits, estimates of broad sense heritability ranged from 20.6 for grain yield to 56.3% for days to heading. Heritability is grouped as low (< 0.2), moderate (0.2-0.4) and high (> 0.4). Broad sense heritability analysis showed that days to heading (56.3), plant height (49.8), Days to maturity (42.4%) and panicle length (42.3%) was high. These findings were in agreement with Adhikari, et al. [9] for days to heading, plant height and days to maturity. Gyawali, et al. [18] reported the same findings for days to flowering and maturity. Thousand seed weight and grain yield showed moderate heritability estimate. Similar results were reported by Gyawali, et al. [18] for grain yield this may be due to the influence of environment on the yield as it is polygenic character. High heritability indicates higher response to selection for specific traits. Very low heritability indicates a greater role of environment on the expression of the trait.

Among the studied traits, estimates of broad sense heritability ranged from 20.6 for grain yield to 56.3% for days to heading. Heritability is grouped as low (< 0.2), moderate (0.2-0.4) and high (> 0.4). Broad sense heritability analysis showed that days to heading (56.3), plant height (49.8), Days to maturity (42.4%) and panicle length (42.3%) was high. These findings were in agreement with Adhikari, et al. [9] for days to heading, plant height and days to maturity. Gyawali, et al. [18] reported the same findings for days to flowering and maturity. Thousand seed weight and grain yield showed moderate heritability estimate. Similar results were reported by Gyawali, et al. [18] for grain yield this may be due to the influence of environment on the yield as it is polygenic character. High heritability indicates higher response to selection for specific traits. Very low heritability indicates a greater role of environment on the expression of the trait.

Correlation coefficient analysis

Correlation coefficient is a measure of the degree association and relationship between two variables. It is important in plant breeding as it can be used for indirect selection. The study of correlation between different characters may help the plant breeder to know how the improvement of one character will bring simultaneous changes in other characters. Positive highly significant phenotypic correlation was observed between days to heading and days to maturity (r = 9385***) indicating selection based on early heading types can be used to select early maturing genotypes, highly significant phenotypic correlation between heading date and maturity date was also reported by Nidhi Tiwari, et al. [19] and Adhikari, et al. [19], but days to heading was negatively associated with thousand grain weight (-0.45979*) and plant height (r = -0.30419. Days to maturity was also negatively correlated with thousand grain weight (-0.52622*) panicle length (-0.49403*), grain yield (-0.044662*) and plant height (-0.33102*).

The weight of a thousand grain (0.50436*) and Panicle length (0.48312*) were positively correlated with plant height. Iqbal, et al. (2018) [17] and Ratna, et al. [20] also reported significant positive correlation between plant height and panicle length. Rice grain yield showed a negative correlation between Days to maturity (-0.044662*), but weak positive correlation was observed for days to heading (r = 0.33472), plant height (r = 0.18729) and panicle length (r = 0.08771). Similar results were found in experiments conducted by Kishore, et al. [10] for panicle length.

Conclusion and Recommendation


Development of high yielding varieties requires the existence of genetic variability and genetic variability of agronomic traits. The success of plant breeding activities entirely depends on the existence of genetic variability with respect to desired traits and selection skill of plant breeder. The results on analysis of variance for yield and yield related characters reveled that highly significant mean square due to genotypes observed for all traits studied, indicating the existence of sufficient variability among genotypes; therefore, there is higher opportunities for effective selection of green super rice for the study area. Higher phenotypic co-efficient of variation, compared to genotypic co-efficient of variation were recorded for all the traits studied indicating the variation is not only genetic but also influence of growing environment. The highest PCV was recorded for the trait plant height followed by panicle length. The greater difference between GCV and PCV was detected in panicle length, plant height, grain yield and thousand grain weight. The result indicated that the influence of environmental factors affects the expression of these traits than genetic factors and can be improved by growing within a proper environment. Broad sense heritability analysis showed that days to heading, plant height, days to maturity and panicle length was high. Thousand seed weight and grain yield showed moderate heritability estimate due to the influence of environment on the yield as it is polygenic character. High heritability with high genetic advance in percent of mean was recorded by plant height followed by Panicle length which confirms their additive gene action and direct selection of these trait depending on phenotypic expression would be effective for yield improvement. Positive and highly significant phenotypic correlation was observed between days to heading and days to maturity Thousand grain weight and Panicle length were positively correlated with plant height. Grain yield showed negative correlation between days to maturity, but weak positive correlation was observed for days to heading, plant height and panicle length Therefore, the traits recorded high heritability and higher genetic advance; and the traits showed positive correlation could be considered as suitable trait for selection and improvement to achieve high yielding green super rice varieties under rain fed environment.

Acknowledgments


The authors are grateful to Fogera Rice research and training Center for supplying green super rice genotypes varieties and financial assistance respectively.

Conflict of Interest


The Author has declared that no conflict of interest exists.

Ethics Approval


Not applicable.

References


  1. Alemu D (2015) Rice in Ethiopia: Progress in production increase and success factors. In: 6th CARD General Meeting 1-22.
  2. MoARD-Ministry of Agriculture and Rural Development (2010) National Rice Research and Development Strategy of Ethiopia. Addis Ababa, Ethiopia 48.
  3. FAO-Food and Agriculture Organization of the United Nations (2018) FAOSTAT Database, Rome, Italy.
  4. Zhang Q (2007) Strategies for developing green super rice. Proc Natl Acad Sci USA 104: 16402-16409.
  5. Li Z, Ali J (2017) Breeding green super rice (GSR) Varieties for sustainable rice cultivation. In: Burleigh Dodds Science Publishing, Cambridge.
  6. Li ZK, Zhang F (2013) Rice breeding in the post-genomics era: From concept to practice. Curr Opin Plant Biol 16: 261-269.
  7. Aditya JP, Bhartiya A (2013) Genetic variability, correlation and path analysis for quantitative characters in rainfed upland rice of Uttarakhand Hills. Journal of Rice Research 6: 24-34.
  8. North S (2013) Estimation of genetic variability and correlation for grain yield components in rice (Oryza sativa L.). Glob J Plant Ecophysiol 3: 1-6.
  9. Adhikari BN, Joshi BP, Shrestha J, et al. (2018) Genetic variability, heritability, genetic advance and correlation among yield and yield components of rice (Oryza sativa L.). Journal of Agriculture and Natural Resources 1: 149-160.
  10. Kishore SN, Srinivas T, Nagabhushanam U, et al. (2015) Genetic variability, correlation and path analysis for yield and yield components in promising rice (Oryza sativa L.) Genotypes. SAARC J Agri 13: 99-108.
  11. Hossain S, Haque M, Rahman J (2015) Genetic variability, correlation and path coefficient analysis of morphological traits in some extinct local Aman rice (Oryza sativa L.). Rice Res 3: 158.
  12. Bandi HRK, Satyanarayana PV, Babu DR, et al. (2018) Genetic variability estimates for yield and yield components traits and quality traits in rice (Oryza sativa L.). Int J Curr Microbiol App Sci 7: 551-559.
  13. Meles B, Mohammed W, Tsehaye Y (2017) Genetic variability, correlation and path analysis of yield and grain quality traits in bread wheat (Tritium aestivum L.) genotypes at Axum, Northern Ethiopia. J Plant Breed Crop Sci 9: 175-185.
  14. Iqbal T, Majeed A, Khattak SI, et al. (2018) Genetic variability and interrelationship studies in green super rice. Res Plant Biol 8: 37-41.
  15. Deshmukh SN, Basu MS, Reddi PS (1986) Genetic variability, character association and path coefficient of quantitative traits in Virginia bunch varieties of ground nut. Indian J Agric Sci 56: 816-821.
  16. Sravan T, Rangare N, Suresh B, et al. (2012) Genetic variability and character association in rainfed upland rice (Oryza sativa L.). J Rice Res 5: 24-29.
  17. Anis G, El-sabagh A, Ghareb A, et al. (2016) Evaluation of promising lines in rice (Oryza sativa L.) to agronomic and genetic performance under Egyptian condition. Int J Agron Agric Res 8: 52-57.
  18. Gyawali S, Poudel A, Poudel S (2018) Genetic variability and association analysis in different rice genotypes in mid hill of western Nepal. Acta Sci Agric 2: 69-76.
  19. Tiwari DN, Tripathi SR, Tripathi MP, et al. (2019) Genetic variability and correlation coefficients of major traits in early maturing rice under rainfed lowland environments of Nepal. Adv Agric 2019: 1-9.
  20. Ratna M, Begum S, Husna A, et al. (2015) Correlation and path coefficients analyses in basmati rice. Bangladesh J Agric Res 40: 153-161.

Abstract


The success of plant breeding activities entirely depends on the existence of genetic variability with respect to desired traits and selection skills of plant breeders. A field experiment was conducted to evaluate the extent of variability, heritability genetic advance, and character association among yield and yield-related traits in lowland GSR. The experiment was conducted using 20 low-land green supper rice genotypes in the Northwestern part of Ethiopia at Pawe agricultural research center for three consecutive years from 2016-2018 main cropping seasons. A randomized complete block design with three replications was used. The results of the analysis of variance for yield and yield-related characters revealed that highly significant mean square due to genotypes observed for all traits studied, Higher phenotypic co-efficient of variation, compared to genotypic co-efficient of variation was recorded for all the traits studied indicating the variation is not only genetic but also influence of growing environment. The highest PCV was recorded for the trait plant height (32.6) followed by panicle length (23.3). Broad sense heritability analysis showed that days to heading (56.3), plant height (49.8), Days to maturity (42.4%), and panicle length (42.3%) were high. High heritability with the high genetic advance in percent of mean was recorded by plant height (49.8 & 33.54) followed by Panicle length (42.3 & 36.55) which confirms their additive gene action and direct selection of these traits depending on phenotypic expression would be effective for yield improvement. The positive and highly significant phenotypic correlation was observed between days to heading and days to maturity (r = 9385***) plant height was positively correlated with thousand-grain weight (0.50436*) and Panicle length (0.48312*). Grain yield showed negative correlation between Days to maturity (-0.044662*), but weak positive correlation was observed for days to heading (r = 0.33472), plant height (r = 0.18729) and panicle length (r = 0.08771). Therefore, the traits recorded high heritability and higher genetic advance; and the traits showed positive correlation and could be considered as suitable traits for selection and improvement to achieve high-yielding green super rice varieties in rain-fed environments.

References

  1. Alemu D (2015) Rice in Ethiopia: Progress in production increase and success factors. In: 6th CARD General Meeting 1-22.
  2. MoARD-Ministry of Agriculture and Rural Development (2010) National Rice Research and Development Strategy of Ethiopia. Addis Ababa, Ethiopia 48.
  3. FAO-Food and Agriculture Organization of the United Nations (2018) FAOSTAT Database, Rome, Italy.
  4. Zhang Q (2007) Strategies for developing green super rice. Proc Natl Acad Sci USA 104: 16402-16409.
  5. Li Z, Ali J (2017) Breeding green super rice (GSR) Varieties for sustainable rice cultivation. In: Burleigh Dodds Science Publishing, Cambridge.
  6. Li ZK, Zhang F (2013) Rice breeding in the post-genomics era: From concept to practice. Curr Opin Plant Biol 16: 261-269.
  7. Aditya JP, Bhartiya A (2013) Genetic variability, correlation and path analysis for quantitative characters in rainfed upland rice of Uttarakhand Hills. Journal of Rice Research 6: 24-34.
  8. North S (2013) Estimation of genetic variability and correlation for grain yield components in rice (Oryza sativa L.). Glob J Plant Ecophysiol 3: 1-6.
  9. Adhikari BN, Joshi BP, Shrestha J, et al. (2018) Genetic variability, heritability, genetic advance and correlation among yield and yield components of rice (Oryza sativa L.). Journal of Agriculture and Natural Resources 1: 149-160.
  10. Kishore SN, Srinivas T, Nagabhushanam U, et al. (2015) Genetic variability, correlation and path analysis for yield and yield components in promising rice (Oryza sativa L.) Genotypes. SAARC J Agri 13: 99-108.
  11. Hossain S, Haque M, Rahman J (2015) Genetic variability, correlation and path coefficient analysis of morphological traits in some extinct local Aman rice (Oryza sativa L.). Rice Res 3: 158.
  12. Bandi HRK, Satyanarayana PV, Babu DR, et al. (2018) Genetic variability estimates for yield and yield components traits and quality traits in rice (Oryza sativa L.). Int J Curr Microbiol App Sci 7: 551-559.
  13. Meles B, Mohammed W, Tsehaye Y (2017) Genetic variability, correlation and path analysis of yield and grain quality traits in bread wheat (Tritium aestivum L.) genotypes at Axum, Northern Ethiopia. J Plant Breed Crop Sci 9: 175-185.
  14. Iqbal T, Majeed A, Khattak SI, et al. (2018) Genetic variability and interrelationship studies in green super rice. Res Plant Biol 8: 37-41.
  15. Deshmukh SN, Basu MS, Reddi PS (1986) Genetic variability, character association and path coefficient of quantitative traits in Virginia bunch varieties of ground nut. Indian J Agric Sci 56: 816-821.
  16. Sravan T, Rangare N, Suresh B, et al. (2012) Genetic variability and character association in rainfed upland rice (Oryza sativa L.). J Rice Res 5: 24-29.
  17. Anis G, El-sabagh A, Ghareb A, et al. (2016) Evaluation of promising lines in rice (Oryza sativa L.) to agronomic and genetic performance under Egyptian condition. Int J Agron Agric Res 8: 52-57.
  18. Gyawali S, Poudel A, Poudel S (2018) Genetic variability and association analysis in different rice genotypes in mid hill of western Nepal. Acta Sci Agric 2: 69-76.
  19. Tiwari DN, Tripathi SR, Tripathi MP, et al. (2019) Genetic variability and correlation coefficients of major traits in early maturing rice under rainfed lowland environments of Nepal. Adv Agric 2019: 1-9.
  20. Ratna M, Begum S, Husna A, et al. (2015) Correlation and path coefficients analyses in basmati rice. Bangladesh J Agric Res 40: 153-161.