Abstract

The associations among yield-related traits and the pattern of influence on rice grain yield were investigated. This evaluation is important to determine the direct and indirect effects of various traits on yield to determine selection criteria for higher grain yield. Fifteen rice genotypes were evaluated under tropical condition at five locations in two planting seasons. The experiment was laid out in a randomized complete block design with three replications across the locations. Data were collected on vegetative and yield components traits. The pooled data based on the analysis of variance revealed that there were significant differences () among the fifteen genotypes for all the characters studied except for panicle length and 100-grain weight. Highly significant and positive correlations at phenotypic level were observed in grain weight per hill (0.796), filled grains per panicle (0.702), panicles per hill (0.632), and tillers per hill (0.712) with yield per hectare, while moderate positive correlations were observed in flag leaf length to width ratio (0.348), days to flowering (0.412), and days to maturity (0.544). By contrast, unfilled grains per panicle (-0.225) and plant height (-0.342) had a negative significant association with yield per hectare. Filled grains per panicle (0.491) exhibited the maximum positive direct effect on yield followed by grain weight per hill (0.449), while unfilled grain per panicle (-0.144) had a negative direct effect. The maximum indirect effect on yield per hectare was recorded by the tillers per hill through the panicles per hill. Therefore, tillers per hill, filled grains per panicle, and grain weight per hill could be used as selection criteria for improving grain yield in rice.

1. Introduction

Rice (Oryza sativa L.) is an important staple food that constitutes a dominant portion of a world standard diet. In spite of its position among the highly rated cereals, the geometric growth rate of the global population has called for improving the current yield for this extremely important cereal. Thus, several methods have been attempted by scientists to combat this perennial problem. Some researchers have attempted nutritional methods such as physiological methods, breeding, and pest and disease control [13]. Among these methods, breeding for high yield traits has been established as the most sustainable method because the traits are heritable. However, yield is a complex trait which is controlled by many factors such as polygene, environment, and genetic variability [4]. Selection for increased grain yield should not be based on yield only because of its complexity and relationship with other yield components. Therefore, other yield-related traits should be taken into consideration.

Path coefficient analysis is a reliable statistical technique used for organising and presenting the cause and effect relationship between the predictor character and response character based on experimental results. It quantifies the interrelationships of different yield components, whether it has a direct influence on the yield or takes another pathway for ultimate effects so that the contribution of each character to yield could be estimated [5]. Path analysis was developed by Wright [6], with a major advantage of portioning the correlation coefficient into direct and indirect effects components. According to Deway and Lu [7], the first component of path analysis is the direct effect of a predictor character upon its response character, while the second component is the indirect effect of a predictor character upon its response character through another predictor character(s).

Path analysis is used in agriculture by plant breeders for identification of characters that can be used as selection criteria for improving grain yield [8]. Except for Bagheri et al. [8], most of the path analyses in rice focusing on grain yield and yield traits considered only a few yield traits. In addition, previous research in path analyses treated yield components traits as the first-order component which resulted in the presence of multicollinearity in path coefficients with values greater than one. The path analysis of grain yield and yield component traits that includes at least second-order components variables in their path diagram is limited in rice. Therefore, the objectives of this study were to acquire and deduce information on the nature of interrelationships between yield and yield component traits among fifteen rice varieties that have been prearranged into first- and second-order predictor character in a path diagram.

Yield is a character that resulted from the association and expression of different yield-related components [4]. Hence, knowledge of the degree of this association through correlation studies can identify traits that could be used as indirect selection criteria for yield or as secondary traits, improving the efficiency of the selection process. If the cause and effect relationship is well understood, then presenting a whole system of the variable in diagram form known as path diagram will be easier [4]. Since correlation coefficient measures only the relationship between two characters and does not actually reveal the relative importance of each trait, this study was conducted to determine the nature of the relationship between grain yield and yield components.

Yield component traits such as plant height, tillers per hill, panicles per hills, panicle length, numbers of filled and unfilled grains per panicle, 100-grain weight, and total grain weight per panicle are important and also essential fundamental task before making any successful breeding program. Path analysis can be used to calculate the quantitative impact on yield of direct or indirect effects caused by one or the other components [9]. Scientists in chilli [4], wheat [10], rapeseed [11], chickpea [5], and cowpea [12] commonly use path coefficient analysis to explain clearly the relations among yield components. Grain yield increment using yield component breeding could be met if the yield component characters are highly heritable and have a positive correlation with the total grain yield. The present investigation was conducted with the following objectives: (i) to find the genetic variability among different plant traits and (ii) to study the interrelationships and influential patterns of some important yield components on rice grain yield by adopting path coefficient analysis tools. The results will then be used as selection criteria for improvement of grain yield in rice under tropical conditions.

2. Materials and Methods

2.1. Plant Husbandry

Fifteen genotypes comprised of six advanced mutant lines, six mutant varieties from abroad (Vietnam and Bangladesh), and three commercial varieties were used in this study. The commercial varieties of MR219, MR220, and MR253 were developed by the Malaysian Agricultural Research and Development Institute (MARDI) and officially released in 2001. These varieties were the first set of varieties to be developed by means of a direct seeding planting system. The emphasis was on the panicle characteristics, short life cycle (105–111 days), fairly long but strong culms, and being tolerant to blast and bacterial leaf blight. VN121, VN124, and VN001 were developed by Vietnam Atomic Energy Institute (VINATOM), and the main characteristics of these varieties include good ratoon potential and semidwarf stature. Three mutant varieties were obtained from Bangladesh Institute of Nuclear Agriculture (BINA), namely, Binadhan4, Binadhan7, and Iratom. These varieties possess high-tillering capacity but are prone to lodging. Six advanced mutant lines (ML4, ML6, ML9, ML10, ML21, and ML24) were promoted from preliminary studies of ion beam irradiation [13] with an extremely high yielding potential. These mutant lines along with mutant varieties from Vietnam and Bangladesh were selected as a representative of combinations of low and high levels of five important yield determinants. These traits include plant height, days to maturity, tillers per hill, number of panicles, and panicle length. The separation between the low- and high-level genotypes was based on our previous research [14].

The field trials were conducted in five locations in two different cropping seasons in peninsular Malaysia, namely, Kota Sarang Semut, Seberang Perai, Tanjung Karang, Sekinchan, and Serdang. The locations covered a wide range of environmental conditions differing in management practices (farmers’ fields versus research field), water system (supplementary irrigation versus rainfed), and temperatures (moderate and warm climate conditions). The locations were chosen to represent major rice producing areas in Malaysia. In each environment, the experiment was laid out in a randomized complete block design with three replications. Plot size was 342 m2, with subplot size of 4 m2 unit for each genotype in each replication. Optimum date for transplanting at each location was followed according to the farmer’s schedule.

Fertilizer application was applied following the recommendation by Malaysian Agricultural Research and Development Institute (MARDI) as nitrogen was applied in form of urea in 30, 55, and 75 days after transplanting at 80, 12, and 20 kg per hectare. NPK fertilizer was also applied in triplicate starting on days 15, 55, and 75 after transplanting at the rate of 140, 107, and 50 kg ha−1. Phosphorus (in form of triple superphosphate) was applied at 15 days at the rate of 57kg per ha−1 and potassium (applied in form of Murate of potash) at 42kg ha−1.

2.2. Data Collection

Five hills were randomly selected for each genotype in each replication to record observations for plant height, flag leaf length to width ratio, days to 50% flowering, days to maturity, tillers per hill, panicles per hill, panicle length, filled grains per panicle, unfilled grains per panicle, 100-grain weight, and grain weight per hill. Yield in ton per hectare (ton ha−1) was estimated from the weight of threshed grains from all panicles in 1.5 × 1.5 m2, excluding border rows.

2.3. Coefficient Analysis

The genotypic and phenotypic correlation coefficient estimates were carried out using SAS, version 9.4 (SAS Institute Inc., Cary, NC, USA). Association of the various characters with yield per hectare was worked out at genotypic and phenotypic levels as described by Kashiani and Saleh [15]. The phenotypic correlations were further partitioned into components of direct and indirect effects using path coefficient analysis according to the method given by Wright [6]. The path coefficients were calculated following Usman et al. [4], where sets of simultaneous equations were arranged in matrix notation that reveals the relationships between correlations and path coefficients as shown in the equation below. In these equations, represents the phenotypic correlations values between variables, while values are the direct effects of one variable upon another and ijij values are the indirect effects. Each observation is defined according to their serial number.1 = Plant height2 = Flag leaf length to width ratio3 = Days to flowering4 = Days to maturity5 = Tillers per hill6 = Panicles per hill7 = Panicle lengths8 = filled grains per panicle9 = Unfilled grains per panicle10 = 100-grain weight11 = Grain weight per hill12 = Yield in t ha−1.

Effects of Vegetative and Yield Component Variables on Yield per Hectarer112=P112+r12P212+r13P312+r14P412+r15P512+r16P612+r17P712+r18P812+r19P912+r110P1012+r111P1112r212=r21P112+P212+r23P312+r24P412+r25P512+r26P612+r27P712+r28P812+r29P912+r210P1012+r211P1112r312=r31P112+r32P212+P312+r34P412+r35P512+r36P612+r37P712+r38P812+r39P912+r310P1012+r311P1112r412=r41P112+r42P212+r43P312+P412+r45P512+r46P612+r47P712+r48P812+r49P912+r410P1012+r411P1112r512=r51P112+r52P212+r53P312+r54P412+P512+r56P612+r57P712+r58P812+r59P912+r510P1012+r511P1112r612=r61P112+r62P212+r63P312+r64P412+r65P512+P612+r67P712+r68P812+r69P912+r610P1012+r611P1112r712=r71P112+r72P212+r73P312+r74P412+r75P512+r76P612+P712+r78P812+r79P912+r710P1012+r711P1112r812=r81P112+r82P212+r83P312+r84P412+r85P512+r86P612+r87P712+P812+r89P912+r810P1012+r811P1112r912=r91P112+r92P212+r93P312+r94P412+r95P512+r96P612+r97P712+r98P812+P912+r910P1012+r911P1112r1012=r101P112+r102P212+r103P312+r104P412+r105P512+r106P612+r107P712+r108P812+r109P912+P1012+r1011P1112r1112=r111P112+r112P212+r113P312+r114P412+r115P512+r116P612+r117P712+r118P812+r119P912+r1110P1012+P1112.

The studied traits were further subdivided into two-stage relations: first-order components and second-order components. The first-order components include plant height, flag leaf length to width ratio, days to flowering, days to maturity, and tillers per hill. The second-order components are panicles per hill, panicle lengths, filled grains per panicle, unfilled grains per panicle, 100-grain weight, and grain weight per hill. The cause and effect relationships between the two components were worked out using additional simultaneous equations arranged in matrix notation as indicated in the equations below.

Effects of First-Order Components on the Panicles per Hill, Panicle Lengths, Grains per Panicle, Unfilled Grains per Panicle, 100-Grain Weight, and Weight per HillPanicles per Hill   r12 r13 r14 r15P56 r21   r23 r24 r25P56 r31 r32 P36 + r34 r35P56 r41 r42 r43   r45P56r56=r51 r52 r53 r54 P56.Panicle Lengths   r12 r13 r14 r15P57 r21   r23 r24 r25P57 r31 r32   r34 r35P57 r41 r42 r43   r45P57 r51 r52 r53 r54 P57.Grains per Panicle   r12 r13 r14 r15P58 r21   r23 r24 r25P58 r31 r32   r34 r35P58 r41 r42 r43   r45P58 r51 r52 r53 r54 P58.Unfilled Grains per Panicle   r12 r13 r14 r15P59 r21   r23 r24 r25P59 r31 r32   r34 r35P59 r41 r42 r43   r45P59 r51 r52 r53 r54 P59.Unfilled Grains per Panicle   r12 r13 r14 r15P59 r21   r23 r24 r25P59 r31 r32   r34 r35P59 r41 r42 r43   r45P59 r51 r52 r53 r54 P59.100-Grain Weight   r12 r13 r14 r15P510 r21   r23 r24 r25P510 r31 r32   r34 r35P510 r41 r42 r43   r45P510 r51 r52 r53 r54 P510.100-Grain Weight   r12 r13 r14 r15P510 r21   r23 r24 r25P510 r31 r32   r34 r35P510 r41 r42 r43   r45P510 r51 r52 r53 r54 P510.Grain Weight per Hill   r12 r13 r14 r15P511 r21   r23 r24 r25P511 r31 r32   r34 r35P511 r41 r42 r43   r45P511 r51 r52 r53 r54 P51.

Effects of Second-Order Components on Yield per Hectare   r67 r68 r69 r610 r611P1112 r76   r78 r79 r710 r711P1112 r86 r87   r89 r810 r811P1112 r96 r97 r98   r910 r911P1112 r106 r107 r108 r109   r1011P1112 r116 r117 r118 r119 r1110P1112.

3. Results and Discussion

3.1. Analysis of Variance

The data presented in Table 1 represent the pooled analysis of variance for the vegetative, yield, and yield component traits for all the genotypes across the entire locations. There was a significant difference among the genotypes, environment, and the genotype by environment interaction as indicated in Table 1. The observed significant differences indicate the presence of considerable amount of genetic variation that exists between the genotypes evaluated. The genetic variability observed in any breeding materials shows the brighter chances of producing desirable traits of plant and perhaps can be used in heterosis breeding [13]. The differences showed by the genotypes could be because of their genetic background and their origination from different source. In this direction, several reports have been published on phenotypic variation among rice genotypes. Pandey et al. [16] reported highly significant differences among 40 rice accessions with the use of 12 quantitative characters. Similarly, using 20 morphological characters, Rao [17] discovered 95% differences among five rice populations.

3.2. Phenotypic and Genotypic Correlation

In this study, the phenotypic and genotypic correlation coefficient of vegetative, yield, and yield component are separated for a clear understanding as shown in Table 2. The proc corr in SAS program provides the r-values and the test of its significance as present in Table 2. The result showed that yield per hectare had a positive correlation with other traits except for unfilled grain per panicle and plant height. Grain yield trait does not exist in isolation but rather as a result of an association with other traits that form a complex relationship that ultimately affects the yield. This association may be either positive or negative. The r-value for Karl Pearson’s correlation coefficient helps in identification of an association between two distinct traits, although it does not measure the magnitude of association but it does give the idea of the relationship. For the correlation coefficient interpretations, Ratner [18] gives a standard accepted guideline. The r-value of 0, +1, and -1 indicates no linear relationship, a perfect positive linear relationship, and negative linear relationship, respectively. The values that range from 0 to 0.3, 0.3 to 0.7, and 0.7 to 1 indicate a low, moderate, and strong positive linear relationships, respectively, while the values that range from 0 to -0.3, -0.3 to -0.7, and -0.7 to -1 indicate a low, moderate, and strong negative linear relationships, respectively.

The correlation coefficient of phenotypic characters ranges from 0.026 to 0.818 while the genotypic character correlation coefficient ranges from 0.04 to 0.919. This indicates that there is higher magnitude at the genotypic level in most cases as compared with the corresponding phenotypic level. The yield per hectare shows a strong, positive, and highly significant correlation with grain weight per hill, tillers per hill, filled grains per panicle, and panicles per hill, respectively. Selection based on these four components is effective due to their equal contribution towards grain yield increment.

3.3. Direct and Indirect Effects of Vegetative Traits on Yield per Hectare

Grain yield is considered as the artefact of all the contributory traits (plant height, flag leaf length to width ratio, days to flowering, days to maturity, tillers per hill, panicles per hill, panicle length, total numbers of grains per panicle, unfilled grains per panicle, and weight per hill), and the correlation coefficient of these contributory factors with final yield are partitioned into direct and indirect effects as presented in Figure 1 and Table 3. These allow the separation of direct and indirect effect through other traits by allotting the correlations for better interpretation of cause and effect relationship [6].

The present study revealed a significant interrelationship among various vegetative and yield components. These traits define the limitation to yield per hectare and the component characters through the direct and indirect effects as a result of interrelationships between them. Therefore, the use of path coefficient analysis investigates the direct and indirect relationships among the component characters through the partitioning of correlation coefficients [6]. The path analysis revealed that the maximum direct effect on yield per hectare was exerted by grain weight per hill followed by tillers per hill and filled grains per panicle; however, the maximum indirect effect on yield per hectare was recorded by the tillers per hill through the panicles per hill (Figure 1). These characters can be used to develop an optimally reliable selection index for realizing improvements in grain yield in rice.

3.4. Two-Stage Relations

The vegetative traits were grouped as the first-order component which includes the following: plant height, flag leaf length to width ratio, days to flowering, days to maturity, and tillers per hill. The second-order component was regarded as yield component traits which were regarded as the principal yield determining factors in rice, and these comprise panicles per hill, panicle length, filled grains per panicle, unfilled grains per panicle, 100-grain weight, and grain weight per hill. The interrelation between these two components is presented in Tables 4 and 5.

3.5. Effect of First-Order Component Relation on Second-Order Component

The path of the influence of first-order component on panicles per hill (Table 4) revealed that only plant height had a negative correlation while other parameters showed a positive relationship. The days to maturity had the highest value of direct correlation with panicles per hill followed by tillers per hill. The path analysis of first-order component with panicle length revealed a negative direct effect on plant height and day to maturity while other parameters such as flag leaf length to width ratio, days to flowering, and tillers per hill show a positive direct effect (Table 4). The interrelationships of plant height, days to flowering, and day to maturity with the grains per panicle showed a negative direct effect while tillers per hill and flag leaf length to width ratio showed a positive direct effect (Table 4). It was reported by Rahman et al. [19] that the flag leaf area was directly related to the yield components. Also, flag leaves play a significant role in enhancing rice yield because these leaves remain the only source of assimilating during the grain filling stage [19]. The larger the flag leaf area is, the more the solar interception and photosynthate productions are, provided that all other factors of production are not limiting. Moreover, tillers per hill had positive and the highest direct effect (0.928) on filled grains per panicle. This is because the tillers per hill play a significant role in determining the yield of the rice grain since it is directly related to panicle number that will be produced per unit ground area. Fewer tillers result in fewer panicles; excess tillers cause high tiller abortions, small panicles, poor grain filling, and reduction in grain yield [14]. The path analysis relationship of first-order component with the unfilled grains per panicle showed that only flag leaf length to width ratio had a positive direct effect while other parameters show a negative direct effect (Table 4). Also, the influence of first-order component on 100-grain weight revealed that plant height and days to maturity showed a negative direct effect while flag leaf length to width ratio, days to flowering, and tillers per hill had a positive and direct effect. The relationship of plant height, flag leaf length to width ratio, days to flowering, days to maturity, and tillers per hill relation with grain weight per hill show positive and direct effect except for plant height which had a negative and direct effect (Table 4). As indirectly pointed out earlier, rice yield is indirectly related to its height. This is due to sink competition for the limited photosynthates produced by limited sources. So what will be used for yield increase will be unnecessarily used for somatic cell enlargement that results in dense vegetative growth and enhanced height.

3.6. Second-Order Component on Yield

The path analysis of second-order component on yield was presented in Table 5. It is revealed that filled grains per panicle (0.491) exhibited the maximum positive direct effect on yield. This was supported by Seyoum et al. [20] who have stated that cultivars with higher grains per panicle showed higher grain yield in rice. Although number of grains per panicle depends on tillers per hill and panicle per hill, its suitability as selection criteria in crop improvement program is mainly dependent on the percentage of filled grains. Grain weight per hill (0.449) also had a positive direct effect on yield followed panicle per hill (0.339). Number of panicles per plant is an important trait contributing to total grain yield as it has a positive relationship to yield. Generally, cultivars with higher panicle number are desirable towards crop improvement. Panicle length (0.056) and 100-grain weight (0.228) also had a positive direct effect but weak relationship while unfilled grain per panicle (-0.144) had a negative direct effect. It could be recalled that grain weight per hill showed the highest genotypic correlation (0.871) with the yield. This strong correlation was as a result of a high positive direct effect on yield. Many previous research works conducted on rice showed similar result [2022]. On the other hand, the highest positive and direct effect has been reported on panicle number [23], days to maturity [24], and spikelet fertility [25]. It could be concluded that the grains per panicle and weight per hill should be given priority during selection in rice improvement program because of their influence on grain yield.

4. Conclusion

The critical analysis of path coefficient analysis and partitioning of correlation reveals that tillers per hill and grain weight per hill possessed positive direct effect and positive association with yield per hectare. Selection for the improvement of grain yield can be efficient if it is based on tillers per hill and grain weight per hill because of their contribution directly towards grain yield. Breeders in these areas should, therefore, develop early maturing genotypes focusing on tillers per hill and grain weight per hill for improving the grain yield per plant in both rainfed and irrigated area.

Data Availability

All data and information used to support the findings of this study are included within the article’s tables and figure.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The authors would like to acknowledge the financial support of IAEA’s Coordinated Research Project (CRP) for Food and Feed (CRP Code D2.30.30) for enhancing sustainable rice production in Malaysia.