# Genetic Diversity of Upland Rice Germplasm in Malaysia Based on Quantitative Traits

Authors: M. Sohrabi; M. Y. Rafii; M. M. Hanafi; A. Siti Nor Akmar; M. A. Latif
Journal: The Scientific World Journal
 (2012)

---

## Abstract

Genetic diversity is prerequisite for any crop improvement program as it helps in the development of superior recombinants. Fifty Malaysian upland rice accessions were evaluated for 12 growth traits, yield and yield components. All of the traits were significant and highly significant among the accessions. The higher magnitudes of genotypic and phenotypic coefficients of variation were recorded for flag leaf length-to-width ratio, spikelet fertility, and days to flowering. High heritability along with high genetic advance was registered for yield of plant, days to flowering, and flag leaf length-to-width ratio suggesting preponderance of additive gene action in the gene expression of these characters. Plant height showed highly significant positive correlation with most of the traits. According to UPGMA cluster analysis all accessions were clustered into six groups. Twelve morphological traits provided around 77% of total variation among the accessions.


---

## Body

## 1. Introduction

Around 3 billion people of the world use rice as a critical or basic food that provides 50 to 80% of their daily calories. Rice is cultivated on more than 150 million hectares, and annual world production is around 600 million tons [1–3]. Upland rice comprises eleven percent of global rice production and is cultivated on around 14 million hectares. Upland rice has a small role in total production but is major food in some tropical countries [4]. Bangladesh, Indonesia, and Philippines are the areas that plant the most upland rice, but the yield is so low (about 1 t/ha on average) and highly variable [5, 6].In Malaysia, two types of rice are cultivated: wetland rice in Peninsular Malaysia (503,184 ha) and upland rice in Sabah and Sarawak (165,888 ha). The average yield of wetland rice is around 3.3 t/ha; in good conditions, however, it can increase to around 10 t/ha. In contrast, the average yield of upland rice ranges from 0.46 to 1.1 t/ha. In 2005, the total national rice production was roughly 2.24 million metric tons. In Malaysia, upland rice is usually cultivated for home consumption by rural people living in Sabah and Sarawak [7].Genetic diversity is the basis of plant breeding, so understanding and assessing it is important for crop management, crop improvement by selection, use of crop germplasm, detection of genome structure, and transfer of desirable traits to other plants [8, 9]. Rice is one of the best plants for the study of genome structure and genetic diversity because it is diploid and has a small genome size of 430 Mb [10], a significant level of genetic polymorphism [11, 12], and a large amount of well-conserved genetically diverse material.The breeders are interested to evaluate genetic diversity based on morphological traits because they are inexpensive, rapid, and simple to score. The study of these traits needs neither sophisticated methods nor complicated equipments, and also these traits can be inherited without either specific biochemical or molecular techniques. Until now scientific classification of plant was based on morphological traits [13, 14]. The rice plant (Oryza sativa) shows great morphological variation, especially in vegetative traits such as plant height and leaf length. Therefore, the present study was undertaken to assess the genetic diversity of upland rice genotypes in Malaysia.


## 2. Materials and Methods

### 2.1. Plant Material and Experimental Design

Fifty accessions of upland rice were selected from MARDI (24 from Peninsular Malaysia and 26 from Sabah). The accessions were cultivated in experimental field of Universiti Putra Malaysia. Sprouted seeds were sown in the pots (Table1). Randomized complete block design (RCBD) with three replications was used with 50 pots for each replication.Table 1
Information on locations, seasons of seed collection and local name of the upland rice accessions.
SLAccessionsLocationSeasonLocal name16040Peninsular MalaysiaMain season (MS)Bedor26041Peninsular MalaysiaOff season (OS)Berjer36043Peninsular MalaysiaMain season (MS)Buih46044Peninsular MalaysiaMain season (MS)Gemalah56045Peninsular MalaysiaOff season (OS)Kura66048Peninsular MalaysiaMain season (MS)Piya76050Peninsular MalaysiaOff season (OS)Ulat86059Peninsular MalaysiaOff season (OS)Rengan bembang96067Peninsular MalaysiaMain season (MS)Lumut/Kuku balam106068Peninsular MalaysiaMain season (MS)Padi Kuku balam116070Peninsular MalaysiaMain season (MS)Selayang126071Peninsular MalaysiaMain season (MS)Lalang137531SabahMain season (MS)Kungkuling A147534SabahMain season (MS)Bukit157535SabahMain season (MS)Pagalan167537SabahMain season (MS)Sibuku177538SabahMain season (MS)Lapaung187539SabahMain season (MS)Sanding197540SabahMain season (MS)Putus tunang207541SabahMain season (MS)Ruabon217543SabahMain season (MS)Semiali227544SabahMain season (MS)Tadaong237545SabahMain season (MS)Tayakon kecil247546SabahMain season (MS)Teun257597SabahMain season (MS)Batangan267596SabahMain season (MS)Kaca277595SabahMain season (MS)Turayo287594SabahMain season (MS)Tarakan297590SabahMain season (MS)Dinabor307589SabahMain season (MS)Rangayat317588SabahMain season (MS)Turakin327585SabahMain season (MS)Peturu337576SabahMain season (MS)Pagalan347571SabahMain season (MS)Turayan357574SabahMain season (MS)Dedawar367575SabahMain season (MS)Lelangsat373824Peninsular MalaysiaMain season (MS)Huma kuning lenggong383825Peninsular MalaysiaMain season (MS)Huma wangi lenggong393826Peninsular MalaysiaOff season (OS)Jarom mas403828Peninsular MalaysiaMain season (MS)Kunyit413830Peninsular MalaysiaMain season (MS)Langsat423831Peninsular MalaysiaOff season (OS)Lenggong433832Peninsular MalaysiaMain season (MS)Puteh perak443834Peninsular MalaysiaOff season (OS)Rambut453833Peninsular MalaysiaOff season (OS)Putih463837Peninsular MalaysiaMain season (MS)Tangkai langsat473838Peninsular MalaysiaMain season (MS)Wangi puteh483835Peninsular MalaysiaOff season (OS)Rengan wangi497508SabahMain season (MS)Beliong507509SabahMain season (MS)Bedumpok


### 2.2. Data Collection

Twelve quantitative traits were recorded for all accessions at each replication: plant height (cm), days to flowering (day), days to maturity (day), flag-leaf-length-to-width ratio (cm), number of tillers per hill (no.), number of grains per panicle (no.), one thousand grains weight (g), yield of plant per pot (g), number of panicles per hill (no.), panicle length (cm), spikelet per panicle (no.), and spikelet fertility (%) (Table2).Table 2
List of quantitative traits of upland rice.
TraitsMethod of evaluationPlant height (PH, cm)The average of height from the base to the tip of last leaf (flag leaf)Days to flowering (DF, days)The number of days from seeding to flowering dayDays to maturity (DM, days)The number of days from seeding to maturing dayFlag leaf length to width ratio (FLR, cm)Dividing the flag leaf length to widthNumber of tillers per hill (NT, no.)Counting of the tillers per hillNumber of grains per panicle (NG, no.)Counting the number of grains on per panicleOne thousand grain weight (1000 GW, g)200 grains were weighted then 1000 weight grains were calculated from these weightsYield of plant per pot (YP, g)Weighting total grains per potNumber of panicles per hill (NP, no.)Counting the panicles per hillPanicle length (PL, cm)From base of the lowest spikelet to the top of latest spikelet on panicleSpikelet per panicle (SP, no.)Counting the spikelet per panicleSpikelet fertility (SF, %)Dividing ripped spikelet to all spikelet


### 2.3. Statistical Analysis

The analysis of variance (ANOVA) revealed the main interaction effects. Least significant difference (LSD) was calculated using Statistical analysis system software (SAS version 9.1) (Table3). Genetic parameters were estimated to identify genetic variation among accessions and to determine genetic and environmental effects on various characters. These genetic parameters were calculated by the formula given by Burton [15], Burton and De Vane [16], and Johnson et al. [17]. These parameters include the following.(a)
Genotypic variance:(1)σg2=MSG-MSEr,
where MSG is the mean square of accessions, MSE is mean square of error, and r is number of replications.(b)
Phenotypic variance:(2)σp2=σg2+σe2,
where σg2 is the genotypic variance and σe2 is the mean squares of error.(c)
Phenotypic coefficient of variance (PCV):(3)PCV(%)=σp2X̅×100,
where σp2 is the phenotypic variance and X̅ is the mean of trait.(d)
Genotypic coefficient of variance (GCV):(4)GCV(%)=σg2X̅×100,
where σg2 is the genotypic variance and X̅ is the mean of character.(e)
Heritability (broad sense):(5)hB2=σg2σp2,
where σg2 is the genotypic variance and σp2 is the phenotypic variance.(f)
Expected genetic advance (GA):(6)GA(%)=K×σp2×hB2×100.
GA is a percent of the mean assuming selection of the superior 5% of accession:
(7)GA(%)=K×σp2X̅×hB2×100,
where K is a constant,σp2/X̅ is the phenotypic standard deviation,hB2 is the heritability, andX̅ is the mean of traits.Table 3
ANOVA showing source of variation, degrees of freedom, means square, and error mean square.
Source of variationdfMSEMSBlocks (r)r-1MSBσe2+gσr2Accessions (g)g-1MSGσe2+rσg2Groups (t)[t-1]MSTσe2+rσg/t2+rgσt2Accessions/groups[t(g-1)]MSG/Tσe2+rσg/t2Error(r-1)(g-1)MSEσe2r: blocks, g: accessions, t: groups, e: error, df: degree of freedom, MS: mean squares, EMS: expected mean squares.The correlation coefficient was analyzed to evaluate the relationships among the different variables in the experiment using SAS software (version 9.1). Data were also analyzed based on Jaccardʼs similarity coefficient by NTSYS-pc software (version 2.1). UPGMA algorithm and SAHN clustering were applied for calculating genetic relationships. The PCA of fifty accessions was calculated by EIGEN and PROJ modules of NTSYS-pc and Minitab software (version 15).


## 2.1. Plant Material and Experimental Design

Fifty accessions of upland rice were selected from MARDI (24 from Peninsular Malaysia and 26 from Sabah). The accessions were cultivated in experimental field of Universiti Putra Malaysia. Sprouted seeds were sown in the pots (Table1). Randomized complete block design (RCBD) with three replications was used with 50 pots for each replication.Table 1
Information on locations, seasons of seed collection and local name of the upland rice accessions.
SLAccessionsLocationSeasonLocal name16040Peninsular MalaysiaMain season (MS)Bedor26041Peninsular MalaysiaOff season (OS)Berjer36043Peninsular MalaysiaMain season (MS)Buih46044Peninsular MalaysiaMain season (MS)Gemalah56045Peninsular MalaysiaOff season (OS)Kura66048Peninsular MalaysiaMain season (MS)Piya76050Peninsular MalaysiaOff season (OS)Ulat86059Peninsular MalaysiaOff season (OS)Rengan bembang96067Peninsular MalaysiaMain season (MS)Lumut/Kuku balam106068Peninsular MalaysiaMain season (MS)Padi Kuku balam116070Peninsular MalaysiaMain season (MS)Selayang126071Peninsular MalaysiaMain season (MS)Lalang137531SabahMain season (MS)Kungkuling A147534SabahMain season (MS)Bukit157535SabahMain season (MS)Pagalan167537SabahMain season (MS)Sibuku177538SabahMain season (MS)Lapaung187539SabahMain season (MS)Sanding197540SabahMain season (MS)Putus tunang207541SabahMain season (MS)Ruabon217543SabahMain season (MS)Semiali227544SabahMain season (MS)Tadaong237545SabahMain season (MS)Tayakon kecil247546SabahMain season (MS)Teun257597SabahMain season (MS)Batangan267596SabahMain season (MS)Kaca277595SabahMain season (MS)Turayo287594SabahMain season (MS)Tarakan297590SabahMain season (MS)Dinabor307589SabahMain season (MS)Rangayat317588SabahMain season (MS)Turakin327585SabahMain season (MS)Peturu337576SabahMain season (MS)Pagalan347571SabahMain season (MS)Turayan357574SabahMain season (MS)Dedawar367575SabahMain season (MS)Lelangsat373824Peninsular MalaysiaMain season (MS)Huma kuning lenggong383825Peninsular MalaysiaMain season (MS)Huma wangi lenggong393826Peninsular MalaysiaOff season (OS)Jarom mas403828Peninsular MalaysiaMain season (MS)Kunyit413830Peninsular MalaysiaMain season (MS)Langsat423831Peninsular MalaysiaOff season (OS)Lenggong433832Peninsular MalaysiaMain season (MS)Puteh perak443834Peninsular MalaysiaOff season (OS)Rambut453833Peninsular MalaysiaOff season (OS)Putih463837Peninsular MalaysiaMain season (MS)Tangkai langsat473838Peninsular MalaysiaMain season (MS)Wangi puteh483835Peninsular MalaysiaOff season (OS)Rengan wangi497508SabahMain season (MS)Beliong507509SabahMain season (MS)Bedumpok


## 2.2. Data Collection

Twelve quantitative traits were recorded for all accessions at each replication: plant height (cm), days to flowering (day), days to maturity (day), flag-leaf-length-to-width ratio (cm), number of tillers per hill (no.), number of grains per panicle (no.), one thousand grains weight (g), yield of plant per pot (g), number of panicles per hill (no.), panicle length (cm), spikelet per panicle (no.), and spikelet fertility (%) (Table2).Table 2
List of quantitative traits of upland rice.
TraitsMethod of evaluationPlant height (PH, cm)The average of height from the base to the tip of last leaf (flag leaf)Days to flowering (DF, days)The number of days from seeding to flowering dayDays to maturity (DM, days)The number of days from seeding to maturing dayFlag leaf length to width ratio (FLR, cm)Dividing the flag leaf length to widthNumber of tillers per hill (NT, no.)Counting of the tillers per hillNumber of grains per panicle (NG, no.)Counting the number of grains on per panicleOne thousand grain weight (1000 GW, g)200 grains were weighted then 1000 weight grains were calculated from these weightsYield of plant per pot (YP, g)Weighting total grains per potNumber of panicles per hill (NP, no.)Counting the panicles per hillPanicle length (PL, cm)From base of the lowest spikelet to the top of latest spikelet on panicleSpikelet per panicle (SP, no.)Counting the spikelet per panicleSpikelet fertility (SF, %)Dividing ripped spikelet to all spikelet


## 2.3. Statistical Analysis

The analysis of variance (ANOVA) revealed the main interaction effects. Least significant difference (LSD) was calculated using Statistical analysis system software (SAS version 9.1) (Table3). Genetic parameters were estimated to identify genetic variation among accessions and to determine genetic and environmental effects on various characters. These genetic parameters were calculated by the formula given by Burton [15], Burton and De Vane [16], and Johnson et al. [17]. These parameters include the following.(a)
Genotypic variance:(1)σg2=MSG-MSEr,
where MSG is the mean square of accessions, MSE is mean square of error, and r is number of replications.(b)
Phenotypic variance:(2)σp2=σg2+σe2,
where σg2 is the genotypic variance and σe2 is the mean squares of error.(c)
Phenotypic coefficient of variance (PCV):(3)PCV(%)=σp2X̅×100,
where σp2 is the phenotypic variance and X̅ is the mean of trait.(d)
Genotypic coefficient of variance (GCV):(4)GCV(%)=σg2X̅×100,
where σg2 is the genotypic variance and X̅ is the mean of character.(e)
Heritability (broad sense):(5)hB2=σg2σp2,
where σg2 is the genotypic variance and σp2 is the phenotypic variance.(f)
Expected genetic advance (GA):(6)GA(%)=K×σp2×hB2×100.
GA is a percent of the mean assuming selection of the superior 5% of accession:
(7)GA(%)=K×σp2X̅×hB2×100,
where K is a constant,σp2/X̅ is the phenotypic standard deviation,hB2 is the heritability, andX̅ is the mean of traits.Table 3
ANOVA showing source of variation, degrees of freedom, means square, and error mean square.
Source of variationdfMSEMSBlocks (r)r-1MSBσe2+gσr2Accessions (g)g-1MSGσe2+rσg2Groups (t)[t-1]MSTσe2+rσg/t2+rgσt2Accessions/groups[t(g-1)]MSG/Tσe2+rσg/t2Error(r-1)(g-1)MSEσe2r: blocks, g: accessions, t: groups, e: error, df: degree of freedom, MS: mean squares, EMS: expected mean squares.The correlation coefficient was analyzed to evaluate the relationships among the different variables in the experiment using SAS software (version 9.1). Data were also analyzed based on Jaccardʼs similarity coefficient by NTSYS-pc software (version 2.1). UPGMA algorithm and SAHN clustering were applied for calculating genetic relationships. The PCA of fifty accessions was calculated by EIGEN and PROJ modules of NTSYS-pc and Minitab software (version 15).


## 3. Result

### 3.1. Variation and Genetic Parameters among Accessions

Eight traits including plant height, days to flowering, flag leaf length-to-width ratio, 1000-GW, yield per pot, panicle length, spikelet per panicle, and spikelet fertility showed highly significant (P≤0.01) variation and the rest of them such as days to maturity, number of tillers per hill, number of grains per panicle, and number of panicles per hill were significant (P≤0.05) among all accessions (Table 4).Table 4
(a)
Mean squares of analysis of variance for 5 traits among 50 accessions of upland rice.
Source of variationdfPHDFDMFLRNTBlocks20.01ns6.61ns41.33ns0.01ns0.0027nsAccessions491295.77**3236.17**5410.90*617.37**2.12*Groups[1]4977.62**15272.12**37117.26**214.76**0.35**Groups/Accessions[48]85.94**276.13**537.95**62.83**1.10**Error9853.43172.58397.5618.220.35*Significant at 0.05. **Highly significant at 0.01. PH: plant height, DF: days to flowering, DM: days to maturing, FLR: flag leaf length-to-width ratio, and NT: number of tillers per hill.(b)
Mean squares of analysis of variance for 7 traits among 50 accessions of upland rice.
Source of variationdfNG1000 GWYPNPPLSPSFBlocks232.03ns0.007ns0.13ns0.090ns0.03ns145.68ns8.39nsAccessions491775.43*155.001**373.65**1.92*129.23**5644.79**4051.39**Groups[1]3998**160.58**2434.58**0.31*41.10**1607.47**351.94**Groups/Accessions[48]454.84**154.88**330.71**1.046**13.18**1269.94**318.59**Error98167.460.1150.1990.3463.78400.5383.21*Significant at 0.05. **Significant at 0.01. NG: number of grains per panicle, 1000 GW: one thousand grain weight, YP: yield per pot, NP: number of panicles per hill, PL: panicle length, SP: spikelet per panicle, and SF: spikelet fertility.In this study, most of the growth traits showed higher PCV compared to yield and yield component traits. However, lower PCV belonged to plant height (15.85%) while flag leaf length-to-width ratio (69.63%) was recorded as higher value. Spikelet fertility (47.31%), days to flowering (40.94%), and days to maturity (40.77%) were recorded as higher values of PCV and number of grains per panicles (21.27%), number of panicle (24.54%), and panicle length (24.63%) showed lower values. The higher GCV was recorded at flag leaf length-to-width ratio (66.66%) and the lower was found at plant height (14.92%). GCV value was low in yield and yield components compared to growth characters. Board sense heritability ranged from 60.26 to 99.84%. The highest and the lowest amount of heritability was recorded at yield of plant and number of panicles, respectively. The estimates of heritability were high for 1000 GW (99.76%), spikelet fertility (94.08%), panicle length (91.69%), flag leaf length-to-width ratio (91.63%), plant height (88.57%), days to flowering (85.54%), spikelet per panicle (81.35%), and days to maturity (80.28%) whereas other characters showed relatively low heritability. GA ranged from 28.93% for plant height to 131.45% for flag leaf length-to-width ratio. The average of GA value in growth traits was higher than the average of GA value in yield and yield components (Table5).Table 5
Genetic variance of 12 morphological characteristics.
TraitsMEANMSGMSEσg2σp2PCV (%)GCV (%)hB2 (%)GA (%)PH136.361295.7753.43414.11467.5415.8614.9288.5728.93DF84.393236.17172.591021.21193.7840.9437.8785.5472.15DM111.555410.90397.571671.112068.6840.7736.6580.7867.85FLR21.20617.3718.23199.71217.9469.6466.6691.63131.45NT3.862.120.350.590.9425.1519.9262.7232.49NG124.651775.43167.47535.98703.4521.2818.5776.1933.401000 GW24.06155.000.1251.6251.7429.8929.8699.7761.44YP43.45373.650.20124.48124.6825.7025.6899.8452.85NP3.811.930.350.520.8724.5519.0660.2630.47PL27.41129.233.7941.8145.6024.6423.5991.6946.54SP159.145644.79400.541748.092148.6229.1326.2781.3548.82SF79.254051.3983.211322.731405.9447.3145.8994.0891.70PH: plant height, DF: days to flowering, DM: days to maturing, FLR: flag leaf length to width ratio, NT: number of tiller per hill   NG: number of grains per panicle, 1000 GW: one thousand grain weight, YP: yield per pot, NP: number of panicles per hill, PL: panicle length, SP: spikelet per Panicle, and SF: spikelet fertility, MSG: mean square of accessions, MSE: mean square of error, PCV: phenotypic coefficient of variation, GCV: genotypic coefficient of variationhB2: board sense heritability, GA: genetic advance, σg2: genotypic variance, and σp2: phenotypic variance.


### 3.2. Association between Traits

Pearson’s correlation coefficient was computed between 12 quantitative traits among 50 accessions of upland rice (Table6). Positive correlation was found between most of traits. Plant height was highly significant and positively correlated with most of traits such as days to flowering, days to maturity, flag leaf width-to-length ratio, number of grains per panicle, yield of plant, panicle length, and spikelet fertility. Yield of plant had highly significant (P<0.01) and positively correlated with plant height (r=0.38), days to maturity (r=0.36), and number of panicles (r=0.48) at 1% probability level and also significant (P≤0.05) and positively correlated with days to flowering (r=0.31) and 1000-grain weight (r=0.34).Table 6
Pearson’s correlation coefficient among 12 quantitative traits of upland rice.
PHDFDMFLRNTNG1000 GWYPNPPLSPSFPH1.00DF0.77**1.00DM0.76**0.90**1.00FLR0.48**0.64**0.48**1.00NT0.120.220.220.271.00NG0.77**0.93**0.99**0.74**0.281.001000 GW0.280.190.150.02−0.020.131.00YP0.38**0.31*0.36**0.060.48**0.34*0.34*1.00NP0.110.200.220.250.99**0.28−0.010.48**1.00PL0.46**0.52**0.45**0.53**0.140.51**0.250.200.151.00SP−0.08−0.13−0.12−0.04−0.23−0.12−0.28*−0.08−0.24−0.101.00SF0.39**0.37**0.29*0.260.180.31*0.230.260.180.26−0.67**1.00*Significantly at 0.05. **Significantly at 0.01. PH: plant height, DF: days to flowering, DM: days to maturing, FLR: flag Leaf length to width ratio, NT: number of tillers per hill, NG: number of grains per panicle, 1000 GW: one thousand grain weight, YP: yield per pot, NP: number of panicles per hill, PL: panicle length, SP: spikelet per panicle, and SF: spikelet fertility.


### 3.3. Cluster Analysis

Fifty accessions of upland rice were clustered into six groups by 12 quantitative traits. As evident from Figure1 and Table 7 cluster III was the biggest (27 accessions) and cluster VI was the smallest (only one member) group. Cluster I, II, IV, and V consisted of 6, 10, 2, and 4 members, respectively. The first group had the highest average in comparison with the other five groups considering five traits (Table 8) such as plant height (147.9 cm), days to flowering (112.8 days), days to maturity (144 days), flag leaf length-to-width ratio (31.1 cm), and panicle length (30.07 cm). Group VI included the highest average for four traits such as number of tillers (4.7), 1000 GW (33 g), yield of plant (55.1 g), and spikelet fertility (95.8%). On the other hand, accessions having this group (VI) showed the lowest average values in the characters such as plant height, days to maturity, flag leaf length-to width ratio, number of panicles, panicle length, and spikelet per panicle.Table 7
Groups of upland rice accessions according to cluster analysis.
GroupsAccessionsGroup I06040, 06041, 06048, 06068, 06070, 06067Group II06043, 06050, 06044, 06059, 06071, 06045, 03835, 03824, 03834, 03838Group III07537, 07571, 07544, 07539, 07538, 07546, 07589, 07597, 07576,07574, 07595, 03825, 03830, 07531, 07535, 07585, 07534, 03831,07509, 07540, 07543, 07590, 07594, 07588, 07575,03832, 03833Group IV03826, 03837Group V07541, 07596, 03828, 07545Group VI07508Table 8
Mean value of 12 quantitative traits for six groups by cluster analysis on 50 upland rice accessions.
GroupPHDFDMFLRNTNG1000 GWYPNPPLSPSFI147.97112.87144.0031.124.50134.3625.9952.424.4030.06148.3090.70II142.8188.73121.8519.263.53141.9726.3248.403.4927.79164.5586.84III132.0477.65102.4019.593.71116.6623.7139.103.6826.97153.2876.68IV142.1096.64129.5824.964.00117.3223.3247.923.9230.07187.4362.62V131.1072.0794.7421.134.40128.7816.0341.904.2925.44199.1664.63VI128.0676.8992.2817.314.75107.0333.0355.134.6722.20111.7095.86PH: plant height, DF: days to flowering, DM: days to maturing, FLR: flag leaf length-to-width ratio, NT: number of tiller per hill NG: number of grains per panicle, 1000 GW: one thousand grain weight, YP: yield per pot, NP: number of panicles per hill, PL: panicle length, SP: spikelet per panicle, and SF: spikelet fertility.Figure 1
The dendrogram of 50 accessions of upland rice based on 12 quantitative traits.


### 3.4. Principal Component Analysis (PCA)

PCA approximately confirmed the cluster analysis for distant accession, 07508, and it was clustered alone in cluster VI (Figure2). But on the other hand, some accessions were close together in PCA such as accessions 06040, 06041, 06070, 06067, 06043, 06050, and 06059, whereas they were clustered into 2 groups (group I and II) in cluster analysis. According to PCA, the first four principal components accounted for around 76.7% of total variation of all morphological traits. The analysis of eigenvectors indicated the information of morphological traits for percentage of variation to the first four principal components, which were 36.4, 17.9, 12.8, and 9.6%, respectively (Table 9).Table 9
Eigenvectors and eigenvalues of the first four principal components.
VariableEigenvectorsPC 1PC 2PC 3PC 4Eigenvalue4.372.141.531.14Variation (%)36.417.912.89.6Cumulative (%)36.454.567.176.7PH0.388−0.227−0.004−0.081DF0.416−0.1790.0760.167DM0.394−0.1400.0930.115FLR0.317−0.1220.2060.399NT0.2270.5410.269−0.014NG0.180−0.3530.114−0.4731000 GW0.1520.007−0.454−0.288YP0.2660.2280.042−0.591NP0.2240.5490.259−0.017PL0.306−0.1350.0080.197SP−0.141−0.3030.607−0.303SF0.2810.038−0.463−0.071PH: plant height, DF: days to flowering, DM: days to maturing, FLR: flag leaf length to width ratio, NT: number of tillers per hill NG: number of grains per panicle, 1000 GW: one thousand grain weight, YP: yield per pot, NP: number of panicles per hill, PL: panicle length, SP: spikelet per panicle, and SF: spikelet fertility, and PC: principal components.Figure 2
Three-dimensional graph of 50 upland rice accessions based on 12 quantitative traits.


## 3.1. Variation and Genetic Parameters among Accessions

Eight traits including plant height, days to flowering, flag leaf length-to-width ratio, 1000-GW, yield per pot, panicle length, spikelet per panicle, and spikelet fertility showed highly significant (P≤0.01) variation and the rest of them such as days to maturity, number of tillers per hill, number of grains per panicle, and number of panicles per hill were significant (P≤0.05) among all accessions (Table 4).Table 4
(a)
Mean squares of analysis of variance for 5 traits among 50 accessions of upland rice.
Source of variationdfPHDFDMFLRNTBlocks20.01ns6.61ns41.33ns0.01ns0.0027nsAccessions491295.77**3236.17**5410.90*617.37**2.12*Groups[1]4977.62**15272.12**37117.26**214.76**0.35**Groups/Accessions[48]85.94**276.13**537.95**62.83**1.10**Error9853.43172.58397.5618.220.35*Significant at 0.05. **Highly significant at 0.01. PH: plant height, DF: days to flowering, DM: days to maturing, FLR: flag leaf length-to-width ratio, and NT: number of tillers per hill.(b)
Mean squares of analysis of variance for 7 traits among 50 accessions of upland rice.
Source of variationdfNG1000 GWYPNPPLSPSFBlocks232.03ns0.007ns0.13ns0.090ns0.03ns145.68ns8.39nsAccessions491775.43*155.001**373.65**1.92*129.23**5644.79**4051.39**Groups[1]3998**160.58**2434.58**0.31*41.10**1607.47**351.94**Groups/Accessions[48]454.84**154.88**330.71**1.046**13.18**1269.94**318.59**Error98167.460.1150.1990.3463.78400.5383.21*Significant at 0.05. **Significant at 0.01. NG: number of grains per panicle, 1000 GW: one thousand grain weight, YP: yield per pot, NP: number of panicles per hill, PL: panicle length, SP: spikelet per panicle, and SF: spikelet fertility.In this study, most of the growth traits showed higher PCV compared to yield and yield component traits. However, lower PCV belonged to plant height (15.85%) while flag leaf length-to-width ratio (69.63%) was recorded as higher value. Spikelet fertility (47.31%), days to flowering (40.94%), and days to maturity (40.77%) were recorded as higher values of PCV and number of grains per panicles (21.27%), number of panicle (24.54%), and panicle length (24.63%) showed lower values. The higher GCV was recorded at flag leaf length-to-width ratio (66.66%) and the lower was found at plant height (14.92%). GCV value was low in yield and yield components compared to growth characters. Board sense heritability ranged from 60.26 to 99.84%. The highest and the lowest amount of heritability was recorded at yield of plant and number of panicles, respectively. The estimates of heritability were high for 1000 GW (99.76%), spikelet fertility (94.08%), panicle length (91.69%), flag leaf length-to-width ratio (91.63%), plant height (88.57%), days to flowering (85.54%), spikelet per panicle (81.35%), and days to maturity (80.28%) whereas other characters showed relatively low heritability. GA ranged from 28.93% for plant height to 131.45% for flag leaf length-to-width ratio. The average of GA value in growth traits was higher than the average of GA value in yield and yield components (Table5).Table 5
Genetic variance of 12 morphological characteristics.
TraitsMEANMSGMSEσg2σp2PCV (%)GCV (%)hB2 (%)GA (%)PH136.361295.7753.43414.11467.5415.8614.9288.5728.93DF84.393236.17172.591021.21193.7840.9437.8785.5472.15DM111.555410.90397.571671.112068.6840.7736.6580.7867.85FLR21.20617.3718.23199.71217.9469.6466.6691.63131.45NT3.862.120.350.590.9425.1519.9262.7232.49NG124.651775.43167.47535.98703.4521.2818.5776.1933.401000 GW24.06155.000.1251.6251.7429.8929.8699.7761.44YP43.45373.650.20124.48124.6825.7025.6899.8452.85NP3.811.930.350.520.8724.5519.0660.2630.47PL27.41129.233.7941.8145.6024.6423.5991.6946.54SP159.145644.79400.541748.092148.6229.1326.2781.3548.82SF79.254051.3983.211322.731405.9447.3145.8994.0891.70PH: plant height, DF: days to flowering, DM: days to maturing, FLR: flag leaf length to width ratio, NT: number of tiller per hill   NG: number of grains per panicle, 1000 GW: one thousand grain weight, YP: yield per pot, NP: number of panicles per hill, PL: panicle length, SP: spikelet per Panicle, and SF: spikelet fertility, MSG: mean square of accessions, MSE: mean square of error, PCV: phenotypic coefficient of variation, GCV: genotypic coefficient of variationhB2: board sense heritability, GA: genetic advance, σg2: genotypic variance, and σp2: phenotypic variance.


## 3.2. Association between Traits

Pearson’s correlation coefficient was computed between 12 quantitative traits among 50 accessions of upland rice (Table6). Positive correlation was found between most of traits. Plant height was highly significant and positively correlated with most of traits such as days to flowering, days to maturity, flag leaf width-to-length ratio, number of grains per panicle, yield of plant, panicle length, and spikelet fertility. Yield of plant had highly significant (P<0.01) and positively correlated with plant height (r=0.38), days to maturity (r=0.36), and number of panicles (r=0.48) at 1% probability level and also significant (P≤0.05) and positively correlated with days to flowering (r=0.31) and 1000-grain weight (r=0.34).Table 6
Pearson’s correlation coefficient among 12 quantitative traits of upland rice.
PHDFDMFLRNTNG1000 GWYPNPPLSPSFPH1.00DF0.77**1.00DM0.76**0.90**1.00FLR0.48**0.64**0.48**1.00NT0.120.220.220.271.00NG0.77**0.93**0.99**0.74**0.281.001000 GW0.280.190.150.02−0.020.131.00YP0.38**0.31*0.36**0.060.48**0.34*0.34*1.00NP0.110.200.220.250.99**0.28−0.010.48**1.00PL0.46**0.52**0.45**0.53**0.140.51**0.250.200.151.00SP−0.08−0.13−0.12−0.04−0.23−0.12−0.28*−0.08−0.24−0.101.00SF0.39**0.37**0.29*0.260.180.31*0.230.260.180.26−0.67**1.00*Significantly at 0.05. **Significantly at 0.01. PH: plant height, DF: days to flowering, DM: days to maturing, FLR: flag Leaf length to width ratio, NT: number of tillers per hill, NG: number of grains per panicle, 1000 GW: one thousand grain weight, YP: yield per pot, NP: number of panicles per hill, PL: panicle length, SP: spikelet per panicle, and SF: spikelet fertility.


## 3.3. Cluster Analysis

Fifty accessions of upland rice were clustered into six groups by 12 quantitative traits. As evident from Figure1 and Table 7 cluster III was the biggest (27 accessions) and cluster VI was the smallest (only one member) group. Cluster I, II, IV, and V consisted of 6, 10, 2, and 4 members, respectively. The first group had the highest average in comparison with the other five groups considering five traits (Table 8) such as plant height (147.9 cm), days to flowering (112.8 days), days to maturity (144 days), flag leaf length-to-width ratio (31.1 cm), and panicle length (30.07 cm). Group VI included the highest average for four traits such as number of tillers (4.7), 1000 GW (33 g), yield of plant (55.1 g), and spikelet fertility (95.8%). On the other hand, accessions having this group (VI) showed the lowest average values in the characters such as plant height, days to maturity, flag leaf length-to width ratio, number of panicles, panicle length, and spikelet per panicle.Table 7
Groups of upland rice accessions according to cluster analysis.
GroupsAccessionsGroup I06040, 06041, 06048, 06068, 06070, 06067Group II06043, 06050, 06044, 06059, 06071, 06045, 03835, 03824, 03834, 03838Group III07537, 07571, 07544, 07539, 07538, 07546, 07589, 07597, 07576,07574, 07595, 03825, 03830, 07531, 07535, 07585, 07534, 03831,07509, 07540, 07543, 07590, 07594, 07588, 07575,03832, 03833Group IV03826, 03837Group V07541, 07596, 03828, 07545Group VI07508Table 8
Mean value of 12 quantitative traits for six groups by cluster analysis on 50 upland rice accessions.
GroupPHDFDMFLRNTNG1000 GWYPNPPLSPSFI147.97112.87144.0031.124.50134.3625.9952.424.4030.06148.3090.70II142.8188.73121.8519.263.53141.9726.3248.403.4927.79164.5586.84III132.0477.65102.4019.593.71116.6623.7139.103.6826.97153.2876.68IV142.1096.64129.5824.964.00117.3223.3247.923.9230.07187.4362.62V131.1072.0794.7421.134.40128.7816.0341.904.2925.44199.1664.63VI128.0676.8992.2817.314.75107.0333.0355.134.6722.20111.7095.86PH: plant height, DF: days to flowering, DM: days to maturing, FLR: flag leaf length-to-width ratio, NT: number of tiller per hill NG: number of grains per panicle, 1000 GW: one thousand grain weight, YP: yield per pot, NP: number of panicles per hill, PL: panicle length, SP: spikelet per panicle, and SF: spikelet fertility.Figure 1
The dendrogram of 50 accessions of upland rice based on 12 quantitative traits.


## 3.4. Principal Component Analysis (PCA)

PCA approximately confirmed the cluster analysis for distant accession, 07508, and it was clustered alone in cluster VI (Figure2). But on the other hand, some accessions were close together in PCA such as accessions 06040, 06041, 06070, 06067, 06043, 06050, and 06059, whereas they were clustered into 2 groups (group I and II) in cluster analysis. According to PCA, the first four principal components accounted for around 76.7% of total variation of all morphological traits. The analysis of eigenvectors indicated the information of morphological traits for percentage of variation to the first four principal components, which were 36.4, 17.9, 12.8, and 9.6%, respectively (Table 9).Table 9
Eigenvectors and eigenvalues of the first four principal components.
VariableEigenvectorsPC 1PC 2PC 3PC 4Eigenvalue4.372.141.531.14Variation (%)36.417.912.89.6Cumulative (%)36.454.567.176.7PH0.388−0.227−0.004−0.081DF0.416−0.1790.0760.167DM0.394−0.1400.0930.115FLR0.317−0.1220.2060.399NT0.2270.5410.269−0.014NG0.180−0.3530.114−0.4731000 GW0.1520.007−0.454−0.288YP0.2660.2280.042−0.591NP0.2240.5490.259−0.017PL0.306−0.1350.0080.197SP−0.141−0.3030.607−0.303SF0.2810.038−0.463−0.071PH: plant height, DF: days to flowering, DM: days to maturing, FLR: flag leaf length to width ratio, NT: number of tillers per hill NG: number of grains per panicle, 1000 GW: one thousand grain weight, YP: yield per pot, NP: number of panicles per hill, PL: panicle length, SP: spikelet per panicle, and SF: spikelet fertility, and PC: principal components.Figure 2
Three-dimensional graph of 50 upland rice accessions based on 12 quantitative traits.


## 4. Discussion

All traits showed highly significant (P<0.01) and significant (P<0.05) variations among 50 accessions, which originated in Peninsular Malaysia and Sabah. Pandey et al. [18] recorded highly significant difference among 40 genotypes of rice with 12 quantitative traits. Wang et al. [19] observed 95% differentiation among 5 populations of rice by 20 morphological traits. Caldo et al. [20] measured highly significant difference (P<0.01) in 41 morphological characteristics between 81 ancestors of rice and also CV ranged from 2.0% for grain length, grain width, and 1000 grains weight to 22.1% for culm number. Chandra et al. [21] and Abarshahr et al. [22] measured highly significant variation at 0.01 revealed by 14 and 19 quantitative traits among 57 accessions of upland rice and 30 genotypes of rice, respectively.Correlation between traits is so important because it helps the breeder to select important characters from the studied traits. Most of the traits such as yield and yield component traits are influenced by interaction of genotype and environment, and, therefore, selection based on correlation coefficient makes it easy for plant breeders [23]. As mentioned, in this assay yield of plant had positive correlation with 12 morphological traits. Lasalita-Zapico et al. [24] evaluated correlation coefficient of 10 quantitative traits for 32 upland rice varieties. In this distinguish significant positive correlation among the majority of the morphological traits was recorded except flag leaf angle that had negative correlation with most of characters such as panicle length, leaf length, leaf width, ligule length, leaf area, and culm length. Zafar et al. [25] recorded positive correlation coefficient of panicle length (yield component) with tillers of plant and 100 grains weight and also significantly positive correlation with grain length (0.278).The computing of heritability and genetic advance useful for selection on phenotypic expression [17]. Therefore, high amount of heritability and genetic advance can be the base of selection according to morphological traits. In present study, flag leaf length-to-width ratio, spikelet fertility, yield of plant, and days to flowering indicated both high heritability and genetic advance. Thus, selection based on these traits would bring about improvement in the genotypes. In previous studies, Sedeek et al. [26] reported both high heritability and high genetic advance for days to heading, flag leaf area, number of filled grains per panicle, and grain yield per plant. The heritability ranged between 86% and 99.4%, and for genetic advance was ranged from 17.81% for number of panicles per plant to 46.16% for grain yield per plant among 24 of rice varieties. Laxuman et al. [27] recorded high heritability (more than 60%) and high genetic advance (more than 20%) for chlorophyll meter reading, number of productive tillers per plant, panicle weight and number of grains per panicle, and 1000 grain weight. Pandey et al. [18] recorded high broad sense heritability among 40 rice varieties for plant height (99.8%), biological yield (99.6%), harvest index (99%), test weight (98.8 g), number of panicles per hills (98.5%), number of spikelets per panicle (98.3%), and grain yield (98.11 g).In our study, fifty accessions of the upland rice were clustered into six groups based on 12 quantitative traits. Ahmadikhah et al. [23] clustered 58 rice varieties into four groups based on 18 morphological traits, and genetic distance was around 0.75. Group A composed of only one member and groups B, C, and D contained 14, 20, and 23 members, respectively. Veasey et al. [28] computed clustering  for 23 populations of rice by 20 morphological characteristics. So the varieties were clustered into 10 groups of the last group was the biggest group with seven members and groups 1, 2, 7, and 8 were the smallest groups including only one variety. So, genotypes having distant clusters could be hybridized to get the higher heterotic responses. The similar studies were reported by several authors [29–31].Principal component analysis indicated diversity among 50 accessions of upland rice by a few eigenvectors. In the present study, the first four principal components indicated 76.4% of total variation for which PC1 showed 36.4% of the variation PC2, PC3, and PC4 explained 17.9%, 12.8%, and 9.6% of total variation, respectively. Lasalita-Zapico et al. [24] computed approximately 82.7% of total variation among 32 upland rice varieties, 66.9% variation for PC1 and 15.87% for PC2. Caldo et al. [20] recorded the first 10 principal components accounting for 67% of total variation. This suggested a strong correlation among characters being examined. Rajiv et al. [32] reported the first two principal components accounting for 82.1% of total variation in control and 68.6% in the stress-induced genotypes.


## 5. Conclusion

Fifty accessions of upland rice were clustered into six main groups. To achieve a wide spectrum of variation among the segregates, genotypes having distant cluster, group I (accessions 6040, 6041, 6048, 6068, 6070, and 6067) could be hybridized with group V (accessions 7541, 7596, 3828, and 7545) and group VI (accession 7508). Principal component analysis indicated 76.4% of the total variation. PCA and cluster analysis complemented each other with some slight inconsistencies in terms of cluster composition. Heritability is one of the most important factors in statistical analysis. Separation and selection of varieties based on high heritability of traits make it easy for breeders. Most researchers agree that high heritability alone is not enough; both high heritability and high genetic advance are needed. In this experiment, flag leaf length-to-width ratio, spikelet fertility, yield of plant, and days to flowering had high heritability and high genetic advance. Most traits such as plant height, yield of plant, panicle length, number of panicles, and days to flowering had positive correlations among each other, which suggested that utilization of these traits could improve the genotype by selection of desirable varieties.



---

*Source: 416291-2012-04-30.xml*