- Research article
- Open Access
Differentiation and description of aromatic short grain rice landraces of eastern Indian state of Odisha based on qualitative phenotypic descriptors
BMC Ecology volume 16, Article number: 36 (2016)
Speciality rice, in general, and aromatic rice in particular, possess enormous market potential for enhancing farm profits. However, systematic characterization of the diversity present in this natural wealth is a major pre requisite for using it in the breeding programs. This study reports qualitative phenotypic trait based characterization of 126 short grain aromatic rice genotypes, collected from different areas of the state of Odisha, India.
Out of the 24 descriptors employed, highest variability (8 different types) was observed for lemma-palea colour with a genetic diversity index (He) of 0.696. The principal component analysis reveals that the tip colour of lemma, colour of awn and colour of stigma, cumulatively explain 74 % of the total variation. The Population STRUCTURE analysis classified the population into two subpopulations which were subdivided further into four distinct groups. The western and southern districts of Odisha are endowed with maximum diversity in comparison to eastern and northern districts and at district level comparisons, Koraput and Puri districts are rich with a genetic diversity values of 0.324 and 0.303 respectively. With this set of morphological qualitative traits, based on ‘phenoprinting’, a newly proposed bar coding system, unique fingerprints of each genotype can be effectively generated that can help in easy identification of these genotypes.
Though aromatic rices represent a tiny fraction of the total rice germplasm, a small collection of 126 land races did exhibit rich diversity for all the qualitative traits. For lemma-palea colour, eight different types were detected while for tip colour of lemma, six different types were recorded, suggesting the presence of rich variability in short grain aromatic rices that are conserved in this region. The proposed ‘phenoprinting’ can be an effective descriptor with the unique finger prints generated for each genotype and coupled with molecular (DNA) finger printing, we can discriminate and identify each and every aromatic short grain rice genotype. The proposed system not only help in conservation but also can confer IPR protection to these specialty rices.
A large and diverse set of aromatic rice landraces are being maintained in different gene banks of India, which have been assembled by various explorations spread over several decades. Many accessions of these accessions have not yet been characterized . Since landraces are thought to be the intermediate group between the wild ancestors and cultivated varieties , characterization of this wealth is of relevance for isolation of desirable genotypes as donors of useful traits. India has rich genetic diversity of aromatic rice landraces, majority of them having small to medium grains [3–6]. Aromatic rices constitute a small but an important sub-group (Group V) identified by Glaszmann  with a set of 15 isozyme markers. Other than Basmati, the indigenous short grain aromatic rices are grown in localized pockets in almost all states of the country [8, 9]. These rices possessing unique aroma, cooking and eating qualities are consumed locally as a delicacy [10, 11]. The aromatic rice varieties have evolved over thousands of years in nature and in the small farms of local farmers who select different types to suit their local cultivation practices and needs .
Odisha, an ancient state of India, is one of the major producer and consumer of aromatic rices. The state has its own set of aromatic short grain rices, which are being cultivated in almost all districts in different agro-climatic zones [6, 13]. These specialty rices are treated as sacred among people. Many of them like Deulabhog, Durgabhog, Krishnabhog, Prasadbhog etc. are being used in several temple duties . The aromatic short grain rice landraces enriched in cooking and eating quality features hold enormous potential to be utilized as value added products, thereby providing higher economic return to small and marginal farmers . Further, these native short grain aromatic rices hold significance over traditional Basmati rice for their varied and intense aroma with long retention in relatively warmer climate . Characterization of these landraces is a prerequisite for understanding the extent of diversity, identification of valuable traits required for aromatic rice improvement and defining the conservation needs. Jeypore tract of Odisha is known for its biodiversity as secondary centre for origin of rice [16, 17]. Characterization of aromatic rice genotypes from Jeypore and its nearby districts can provide insights into the origin and spread of these rices. Since, rice is a self-pollinated crop and morphological description of rice plant is well established, its phenotypic descriptors could serve as potential marker for varietal identification and differentiation. Phenotypic characterization of the aromatic landraces would be a primary but essential measure leading to cataloguing of these unique high quality and diverse group, thereby strengthening their conservation and utilization strategies.
Despite that Odisha holds significance from the view point of origin and spread of rice, the status of aromatic short grain rice genotypes of this Indian state has not been studied extensively. The present study investigates structure and diversity among the aromatic short grain rice genotypes (ASGs) collected from different districts of Odisha.
Plant material collection and characterization
One hundred and twenty-six short grain aromatic rice lines (grains/panicles) were collected by germplasm experts and plant breeders of Central Rice Research Institute (CRRI), Cuttack from different parts of Odisha, India (Additional file 1: Table S1) during an exploration and collection programme between 2005 and 2009 conducted by the Institute. Further, the genotypes were grouped under 19 geographical districts based on their area of collection (Additional file 2: Figure S1). The genotypes were transplanted in randomised complete block design in different experimental plots of CRRI and the genetic purity of each landrace was maintained by periodic removal of ‘off types’ in four purification cycles. The experiment was conducted under controlled environmental condition with recommended dosage of fertilizers to minimize any environmental error. The purified rice genotypes were grown in 2 m × 2 m plot, from which 10 hills were sampled for characterization. Twenty-four phenotypic characters were measured at different growth stages following the guidelines of International Rice Research Institute for three consecutive years i.e. 2013–2015 .
Data for 24 qualitative phenotypic traits were recorded on 10 random plants per accession for 3 years and their means were calculated for analysis (Table 1). The qualitative phenotypic traits were coded for the presence (1) and absence (0) of each alternate form. Diversity parameters viz., number of alleles/variables (Na), effective number of alleles/variables (Ne), Shannon Index (I) and unbiased Nei’s genetic diversity index (He)  were calculated using POPGENE v 1.32  with 1000 permutations. The neighbour joining dendrogram based on Nei’s unbiased pairwise genetic distances among genotypes was constructed in MEGA 6 . District wise genetic diversity measures were calculated using POPGENE v 1.32 and the neighbour joining tree was constructed with pairwise genetic distance measures between districts by MEGA 6. The molecular variance was analysed by hierarchical AMOVA implemented in the software GeneAlEx6  based on unbiased Fst estimator as described by Weir and Cockerham  to estimate genetic variation among and within districts. Significance of the F-statistics was tested by non-differentiation of probability for 10,000 randomizations. The Bayesian model-based clustering analysis of the genotypes was used for determining the optimal number of genetic clusters found among rice varieties by the software STRUCTURE  using admix approach with 100,000 burn-in (iteration) periods followed by 100,000 Markov Chain Monte Carlo (MCMC) replicates with ten independent runs (K) ranging from 1 to 10. The ΔK based on the change in the log probability of the data between successive K values was estimated using the parameters described by Evanno and colleagues  with Structure Harvester v6.0  and the inferred population clusters were produced by Structure Plot  (http://www.btismysore.in/strplot). Individuals with probability of membership >80 % were assigned to the same group, while those with <80 % probability membership in any single group were assigned to a admixed group . We used a graphic presentation method based on alternate forms for the morphological descriptors, which we describe as ‘Phenoprint’. The Phenoprint was developed using Microsoft Excel 2013, in which the binary data were represented as bar(s); black bar for the presence of the variable and grey bar for the absence to provide uniqueness for each of the genotype. Different morphological characters that contributed most to the observed phenotypic variance were identified by principal component analysis (PCA) in the software SAS 9.2 . The subset regression was used to identify different principal components that contributed to the total variance in the dataset.
Phenotypic trait diversity
Eighteen (75 %) out of 24 phenotypic traits showed variation among the 126 indigenous short grain aromatic rice genotypes, whereas no variation was observed for the traits like intensity of blade colour (LI), leaf anthocyanin colour (LA), pubescence of blade surface (LPB), leaf ligule (LL), colour of ligule (CL) and collar colour (CC) (Table 1). Phenotypic variation for some of the traits are given in Fig. 1. The number of variables per individual trait ranged from 1 (traits having no variation) to a maximum of 8 for lemma-palea colour with a total number of 70 variables with an average of 2.92 per trait. The genetic diversity index (He) for different traits ranged from 0.016 (lowest) to 0.696 (highest) for culm angle and lemma-palea colour, respectively with an average of 0.286. The number of effective variables (Ne) ranged from 1 to 3.292 with an average of 2.958 (Table 2). Further, frequency distribution of the phenotypic traits against the studied genotypes indicated unequal distribution of variables under each trait (Additional file 2: Figure S2). Among the 126 genotypes, presence of auricles (99 %), semi erect culm angle (99 %), split ligule (95 %), absence of awn (91 %), drooping type panicles (91 %), semi erect attitude of panicle branches (86 %), well exserted panicle (76 %), white colour stigma (74 %), light gold colour internodes (71 %), yellowish white colour of awns (64 %), branching of panicles weak (63 %), medium stem thickness (55 %), green basal leaf sheath colour (51 %), lemma having purple tip (45 %) and erect flag leaf attitude of blade (44 %) were predominant. A comparatively equal frequency of distribution for short bold (47 %) and medium bold (51 %) grain type were recorded for the genotypes. In case of lemma-palea colour that showed maximum variation with eight alternative variables, straw colour was most predominant (50 %) followed by purple furrow (14 %), brown (11 %), purple and brown furrow (10 % each), black and red (2 % each) and brown spot (1 %).
The principal component analysis (PCA) showed that 3 of the 24 morphological traits were the most important components for explaining the grouping of genotypes (Fig. 2). Principal component analysis showed that, first, second, third, fourth, and fifth principal components accounted for 46, 20, 7, 5 and 4 % of the total variance, respectively. The first three principal components together explained a cumulative 74 % of the total variance. Further, the subset regression analysis could identify tip colour of lemma, colour of awn and colour of stigma as the three most important traits. Grain type, flag leaf attitude of blade and lemma-palea colour were identified to be secondary traits to explain the maximum available phenotypic variability within the set of studied genotypes. However, ten (tip colour of lemma, colour of awn, colour of stigma, grain type, flag leaf attitude, lemma palea colour, attitude of panicle branches, basal leaf sheath colour, leaf angle and internode colour) out of 24 morphological traits could adequately explain the total variation existing in these short grain aromatic rice genotypes (Table 2).
Genetic relationship and population structure
The pairwise Nei’s genetic distance was calculated to understand the genetic diversity and relatedness among the 126 short grain aromatic rice genotypes, which ranged from 0.426 (lowest) to 0.875 (highest) with an average of 0.347. The neighbour-joining dendrogram constructed based on Nei’s unbiased pairwise genetic distance grouped the rice genotypes into 2 major clusters (Fig. 3a), of which, Cluster I consisted of 115 genotypes while the rest 11 genotypes grouped in Cluster II. Since, the genotypes differed from each other for at least one or more morphological characters, no duplicates with identical features were detected in the entire set of germplasm.
Bayesian analysis of population structure using the model-based approach with the admixture model proposed by Falush and colleagues  provided support for the existence of population structure. The maximum log likelihood (84.76) was observed for K = 2 (Additional file 2: Figure S3). The genotypes were differentiated into two subpopulations i.e. SP 1 with a membership percentage of 54.8 % while SP 2 with 45.2 %. The fixation index (Fst) values of subpopulations were 0.1415 and 0.2794 for SP 1 and 2, respectively, while the pairwise allele frequency divergence value between the subpopulations was 0.0489 (Additional file 1: Table S2), indicating the existence of moderate population structure in these genotypes. Further, 8 and 7 genotypes were identified to be admix type having a membership value of <80 % in the subpopulations SP1 and SP2, respectively. Moreover, at K = 4, the 126 short grain aromatic rice genotypes could be further differentiated into four groups. (Fig. 3c, d). Forty-four genotypes were clustered in Group II followed by 38 in Group III, 23 in Group IV and 21 genotypes in Group I. Relationship among the genotypes obtained from the STRUCTURE based clustering was in line with the NJ dendrogram except a few deviations.
The phenoprints of the aromatic short grain rice genotypes were developed based on 70 variables of the 24 morphological characters (Fig. 3b). The phenoprints uniquely identified each landrace and differentiated them from each other for at least one morphological character. Moreover, it provided a clear depiction of variation for identification of genotypes based on their morphological attributes. For instance, when phenotypic trait of different genotypes were compared, their phenoprints depicted Atmasital-1 having purple lines for basal leaf sheath colour, while light purple in Dangar Basumati, and green in Deulabhog-1. The probability of identical match by chance based on all the polymorphic traits was 2.9 × 10−4, with a differentiation possibility of 104 genotypes. This suggested that the phenoprints can be used as identifying tags for the ASGs (aromatic short grain rices). Different accessions with the same name like Deulabhog, Lektimachi, Acharmati, Atmasital, Ganjamlocal, Kalikati, Basnadhan, Basumati, Laxmibilas, Samleibhog, Basuabhog and Kalajeera collected from different geographical areas differed from each other for more than one qualitative morphological characters. For instance, the two accessions of Samleibhog (collected from Sundargarh area) were different for the characters like panicle attitude of branches, panicle secondary branching, lemma-palea colour and grain type. The three accessions of Deulabhog (collected from Puri and Koraput) showed variation for a maximum of 14 out of 24 morphological traits.
Further, we performed a heuristic search and identified 16 genotypes (12.7 %) i.e. Deulabhog-1 (Puri), Lektimasi (Malkangiri), Jayaphul (Sundargarh), Karpurabasa (Koraput), Lektimachi-1 (Malkangiri), Pirima (Koraput), Kalajeera-1 (Mayurbhanj), Gadakakudinga (Phulbani), Garmatia (Puri), Saragadhuli (Cuttack), Ganjeikali (Dhenkanal), Kalazeera (Dhenkanal), Basanadhan (Koraput), Dhurabahila (Koraput), Sunsuniasunaphul (Deogarh) and Kalakanhu (Bolangir), which represented the maximum available phenotypic variation (more than 98 %) for the 24 phenotypic traits recorded. Higher numbers of genotypes (6) in this group were recorded from Koraput and Malkangiri region. The average Nei’s genetic diversity value of 0.245 detected for this group of genotypes was comparable with genetic diversity value of the complete set (0.286).
Differentiation of aromatic short grain rice genotypes based on geographical districts of collection
The aromatic short grain rice genotypes collected and screened in our study were represented by 19 different districts of Odisha. We grouped the genotypes based on these 19 districts to address the phenotypic variability within a particular district. In our collection, maximum (17) number of genotypes were from Koraput followed by Cuttack (15) and Ganjam (12) whereas, only limited number of genotypes (2 each) were included from Anugul, Jajpur, Kendrapara and Phulbani districts. We could detect a non-uniform but wide distribution of phenotypic characters like basal leaf sheath colour, flag leaf attitude of blade, tip colour of lemma, colour of stigma, lemma-palea colour and grain type in all the districts (Fig. 4). However, a comparatively uniform but narrow distribution of phenotypic characters was recorded for culm angle (exception of Koraput) and leaf auricles (exception of Dhenkanal). In case of panicle secondary branching, all the district except Phulbani contained accession with both the alternate forms. A high level of variation was observed in Koraput and Puri districts, whereas it was significantly low in Anugul, Jajpur, Kendrapara and Phulbani districts. Genetic diversity parameters revealed that effective number of variables was high for Puri district (1.703) as compared to Anugul (1.167) (Table 3). Similarly, the highest value for Nei’s genetic diversity was recorded in the germplasm set of Koraput (0.324) followed by Puri (0.303), Anugul having the lowest (0.083). The values for Shannon’s information index were in line with the genetic diversity values. However, when we calculated the polymorphism percentage of genotypes for the 24 qualitative phenotypic traits, greater part of the variation was observed in Koraput district (67.61 %) followed by Puri (64.79 %), while, the genotypes from Anugul district represented a minimum variation of 8.45 %. The pair-wise comparisons of Nei’s unbiased genetic distance between the 19 districts estimated based on the qualitative morphological traits ranged from 0.370 (between Phulbani and Kendrapara) to 0.013 (between Kalahandi and Cuttack) with an average distance value of 0.097 (Additional file 1: Table S3). The neighbour-joining tree constructed based on Nei’s genetic distance grouped the 19 geographical regions into two major clusters (Fig. 5). The six districts i.e. Sambalpur, Kalahandi, Deogarh, Malkangiri, Bolangir and Sundargarh, which are close to the border of nearby state, Chhattisgarh on the western side of Odisha, were grouped in one cluster (I) with only exception of Nayagarh. The rest of the districts were grouped in cluster II. Genotypes from Puri, Kendrapara, Balasore, Koraput, Mayurbhanj, Dhenkanal, Ganjam, Anugul, Keonjhar, Jajpur and Phulbani showed close relationship among them. Further, analysis of molecular variance (AMOVA) based on geographical districts indicated that on an average, 92 % variation exists within districts and 8 % variation was observed among the district (Table 4).
Assessment of genetic variation in the short grain aromatic rices of a particular region has implications for conservation of these neglected heritage resource and unlocking genes controlling valuable traits for their utilization in rice improvement . Further, the state of Odisha being a major producer of rice in India and having a large assembly of indigenous rice genotypes , there was an imperative need to characterize them systematically.
Overall phenotypic variation
In this study, we could identify eighteen, out of 24 phenotypic traits which showed significant variation in the ASG rice collection. A total of 70 variables were detected by 24 phenotypic traits in our study, which is significantly higher than that in aromatic rice germplasm of Madhya Pradesh and Chhattisgarh of India , in upland rice collection of Japan , in aromatic rice collection from Maharashtra and Karnataka of India , in aromatic short grain rice genotypes of eastern India  and in indigenous rice landraces of Chhattisgarh state of India . Highest genetic diversity value of 0.696 was contributed by lemma-palea colour, where 8 different forms were identified which is higher than the reports of Sarawgi and colleagues , who reported 5 types in 782 rice germplasm collection. Since lemma and palea are sepal equivalent , genetic base of lemma-palea colour differentiation could be utilized as a potential tool to understand the evolutionary aspects of these native genotypes. Lemma-palea pigmentation is the character of wild rice and straw hull is similar to indica and black hull is similar to that in aus groups . Black husk4 (Bh4) gene, encoding an amino acid transporter, fine mapped on chromosome 4, controls the black husk in O. rufipogon and O. nivara. Zhu and colleagues  reported that transition from black husk in wild rice to straw husk in cultivated rice is due to a functional deletion of 22 bp in exon 3 of bh4 variant. The intermediate types, which were observed for this trait, could be useful to establish a lineage between aromatic rice landraces and their wild progenitors to understand their domestication process. Further, red and black husked rices having nutritive and medicinal properties, are more resistant to storage insect and pests as compared to brown husked rices . Ray and colleagues  reported that aromatic short grain rices vary for their grain characters, which was observed in our study also, where 4 different grain types were identified. Similar to our results, high level of variation was reported for grain shape in Brazilian and Pakistan rice collections [33, 43].
In our study, genotypes having same name collected from different districts (both nearby and far apart) showed variation for one or more traits. This observation is interesting. This could be due to the naming practice being followed by traditional communities for the ASGs. Since they are offered to Gods and Goddesses in the temples, the different material in different locations might be given similar/identical names. For example Deula Bhog (Deula meaning temple; Bhog meaning offering) samples collected from Puri and Koraput districts differ for basal leaf sheath colour, flag leaf attitude of blade stem thickness, panicles characteristics, tip colour of lemma, panicle exertion, shape of ligule, colour of stigma, internode colour, leaf angle and grain type. Similarly, Laxmibilas collected from Sambalpur and Deogarh districts differ for basal leaf sheath colour, flag leaf attitude of blade, tip colour of lemma, panicle exertion, colour of stigma, internode colour, panicle secondary branching, leaf angle, lemma-palea colour and grain type. The other reason of this situation arising could be introduction of the same landrace in different locations followed by selection by local growers or farmers leading to fixation of different genotypic constitution at different locations. It has been demonstrated by Kohli and colleagues  using molecular marker and other studies that the landraces do possess composite genetic structure and selection exercised on the base landrace is expected to result in genetically different purelines. However, what remains common in the developed lines is the target trait i.e. short grain and aroma. In our earlier report also we have detected different grades of aroma (mild/moderate/strong) in these genotypes having similar name . Similarly, Bisne and Sarawgi  have found variation for leaf blade colour, lemma-palea colour, apiculus colour and lemma-palea pubescence in Badshah Bhog lines. Genetic variability in Indian scented rice accessions of Hansraj were also reported by Raghunathachari and colleagues  by using RAPD markers.
In 6 traits (leaf intensity of blade colour, leaf anthocyanin colour, leaf pubescence of blade surface, leaf ligule, leaf colour of ligule and collar colour), we could not detect any variation for the set of aromatic rice genotypes used in this study. Of the rest 6 traits, 4 related to colour of blade, leaf, collar and ligule which was green in all the genotypes. Since anthocyanin pigmentation in different parts of rice plant is due to allelic variation and complex organization of gene(s) , this uniformity for colour of the four traits could be due to pleiotropic effect which has remained unchanged during the process of domestication and selection by farmers .
The average genetic diversity value recorded in this study was relatively low (He = 0.286). Similar type of result was also observed in the all India collection of aromatic rice landraces . However, the genetic diversity value was comparatively high than the rice landraces of Nepal . The Shannon diversity index of 0.515 detected in our germplasm set is comparatively high than the rice landraces of Santhal Paraganas of Jharkhand state, India . We could identify six traits i.e. tip colour of lemma, colour of awn, grain type, flag leaf attitude of blade, lemma-palea colour and colour of stigma as the most important traits to explain the variation and these traits in combination could be utilized as key for differentiation of short grain rice genotypes at field level.
Bayesian analysis of STRUCTURE using the model-based approach suggested the presence of two optimum populations i.e. SP1 and SP2 (with a K value of 2) in these natural heritages. Choudhury and colleagues  also inferred 2 subpopulations in the Eastern Himalayan and Northeast Indian indigenous rice with a set of seven microsatellite markers. Similarly, 3 subpopulations were detected in Indian rice core set with molecular markers, which reveals presence of moderate population structure in the rice genotypes of India . A total of 15 genotypes were determined to be admixture type, which might be due to intraspecific crosses between genotypes of SP1 and SP2 . Further, the two identified populations were divided into four different groups similar to the results of Laido and colleagues , where similar type of groups were identified in two STRUCTURE populations in the tetraploid wheat genotypes with SSR and DArT markers.
The aromatic rice collection was found to be a reserve of many useful agronomic traits. Erect flag leaf is one of the important features and related to high yielding ability . In our collection, 44 % of genotypes had erect flag leaf and these could be potential donors of the trait in high yielding aromatic rice breeding programmes. Tall plant type is prone to lodging, unresponsive to fertilizer application and hence has low yielding ability . In our study, 25 genotypes were identified with thick/strong stem, which could be deployed to achieve lodging resistance under higher dose of fertilizers . Similarly, 74 % of the genotypes identified with well exserted panicle type could be utilized for introgession of this trait into high yielding backgrounds. Further, 37 % of the genotypes having strong panicle secondary branching, which is another important trait related to yield, could be utilized in combination with other yield traits to overcome the low yield potential of aromatic rices .
Trait specific characterization of these lines would help in identifying novel genes/alleles to breed improved disease resistant quality rices. This is highly relevant in the context that in the past few decades, significant decline in cultivation and huge erosion of short grain aromatic rice germplasm in Odisha has been a major concern . Artificial selection pressure and promotion of high yielding varieties and hybrids over years has led to reduced genetic base and increased genetic vulnerability in traditional rice germplasm, emphasizing the conservation, characterization and documentation of native aromatic genotypes before extinction [61–63]. The present study that describes variation in ASGs is a relevant step in this direction.
Phenoprint of aromatic short grain rice genotypes
Bar coding tools have proven to be significantly useful for varietal differentiation and their identification, almost all of which are based on DNA markers. After the establishment of a protocol for DNA barcoding of all land plants , this technology has been extensively used by several researchers for barcoding of many edible plants [65–70]. However, potential application of morphology based characterization cannot be ignored, since it is having useful implication to varietal discrimination at field and farm level . And hence, in this study, we have characterized 126 aromatic short grain rice genotypes of Odisha, India based on qualitative phenotypic descriptors and proposed a unique bar coding approach, named ‘phenoprinting’, for varietal description. Once established, phenoprint could be utilized for rice and other crops as well. Given that the probability of identity based on all the polymorphic traits was low, this approach would be of use in varietal identification, farmers’ rights protection and intellectual property rights issues. Since, the phenoprint deals with morphological characters, it could be generated with minimal financial inputs. The phenoprint in our study that could clearly differentiate the genotypes will be a useful guide for visual comparison of additional aromatic short grain genotypes. In fact, the traits having no variation could be easily distinguishable and the traits with maximum variability could also be precisely identified. This approach provided uniqueness to each of the lines tested, providing robustness for trait specific varietal identification. Further, the 16 genotypes identified in our study, containing more than 98 % of the variables for phenotypic traits provided the key, which could be utilized for description of ASGs and also in the rice improvement programme to overcome the genetic bottleneck in years to come.
Phenotypic diversity in different districts
The tribal dominated Koraput is endowed with diverse and valuable flora and fauna . Further, Koraput region has been recently declared as ‘Globally important agricultural heritage systems’ by the Food and Agricultural Organization of the United Nations (FAO) . Our findings, revealing the highest polymorphism percentage of more than 67 % and greater diverse group of aromatic rices in Koraput as compared to other areas, justifies the National and International attention towards Koraput for various biodiversity projects. The genetic diversity parameters were similar between Koraput and Puri but were significantly higher than the other districts. Detection of high genetic diversity in the germplasm set of Puri was interesting. Puri district of Odisha is renowned for the temple of Lord Jagannath where the famous car festival takes place annually. Besides, it also houses several other temples, monasteries and holy ashrams where aromatic rices are used in different temple rituals. Hence, import of aromatic rices from other parts of the state resulting in an assembly of a diverse group of genotypes in this district is possible. Further, 8 % of the total variation observed among districts suggested that each district have their unique group of aromatic short grain rices.
The neighbour joining tree based on pairwise genetic distance estimates, grouped the six districts of western Odisha in one group. Since, western Odisha, Jharkhand and Chhattisgarh are recognized as the centre of origin of aus ecotypes of rice , close relationship among the districts of western Odisha reflects evolution of short grain rice landraces in this region. In fact, critical observations showed that number of aromatic short grain rice genotypes and their genetic diversity is comparatively on higher side in western and southern Odisha than eastern and northern Odisha. Fuller [75, 76] reported an independent domestication of rice—Neolithic in Ganga Valley and the western Odisha with archaeo-botanical evidences. Considering Koraput in western Odisha as the secondary centre of origin for rice, greater diversity in Koraput and adjoining region is expected.
Native aromatic rice landraces that are highly preferred by consumers needs to be characterized which can help in varietal diagnostics purpose and their conservation . Here, we have characterised native short grain aromatic genotypes of Odisha (India), one of the major consumer and producer of these native rices, based on qualitative phenotypic descriptors. For the trait lemma-palea colour, a highest of 8 different forms were detected followed by tip colour of lemma where 6 variables were identified. Tip colour of lemma, colour of awn and colour of stigma were the primary determinants for explaining the existing variation in this group of aromatic rice genotypes. The 126 genotypes were broadly grouped into 2 sub populations and 4 distinct groups by population STRUCTURE analysis, revealing the existence of moderate population structure within them. But, we did not detect any duplicates in our set. Further, western and southern districts of Odisha had maximum diverse aromatic short grain rice genotypes as compared to eastern and northern districts. Moreover, the proposed phenoprinting approach, discriminating the aromatic short grain rices based on qualitative phenotypic descriptors, could provide unique identification and description to the genotypes. Since, genetic erosion over years is a major threat to aromatic rice improvement, genetic differentiation and phenotypic description study may help in preservation of this group of indigenous short grain aromatic rices.
Rana JC, Negi KS, Wani SA, Saxena S, Pradheep K, Pareek SK, Sofi PA. Genetic resources of rice in the Western Himalayan region of India: current status. Genet Resour Crop Evol. 2009;56:963–73.
Pusadee T, Jamjod S, Chiang YC, Rerkasem B, Schaal BA. Genetic structure and isolation by distance in a landrace of Thai rice. Proc Nat Acad Sci USA. 2009;106(33):13880–5.
Singh RK, Gautam PL, Saxena S, Singh S. Scented rice germplasm: conservation, evaluation and utilization. In: Singh RK, Singh US, Khush GS, editors. Aromatic rices. Kalyani: New Delhi; 2000. p. 107–33.
Roy JK, De RN, Ghorai DP, Panda A. Collection and evaluation of genetic resources of rice in India. Phytobreedon. 1983;1:1–9.
Dikshit N, Malik SS, Mohapatra P. Seed protein variability in scented rice. Oryza. 1992;29:65–6.
Malik SS, Dikshit N, Das AB, Lodh SB. Studies on morphological and physio- chemical properties of local scented rice. Oryza. 1994;31:106–10.
Glaszmann JC. Isozymes and classification of Asian rice varieties. Theor Appl Genet. 1987;74:21–30.
Shobha Rani N, Singh RK. Efforts on aromatic rice improvement in India. In: Singh RK, Singh US, editors. A treatise on the scented rices of India. New Delhi: Kalyani Publishers; 2003. p. 23–72.
Shobha Rani N, Krishnaiah K. Current status and future prospects for improvement of aromatic rices in India, In: Specialty rices of the world: Breeding, production and marketing, New York: Science Publishers, Inc.; 2001. p. 49–78.
Bhattacharjee P, Singhal RS, Kulkarni PR. Basmati rice: a review. Int J Food Sci Tech. 2002;37:1–12.
Ahuja SC, Pawar DVS, Ahuja U, Gupta KR. Basmati rice—the scented pearl. Hissar: CCS Haryana Agricultural University; 1995. p. 1–63.
Pachauri V, Singh MK, Singh AK, Singh S, Shakeel NA, Singh VP, Singh NK. Origin and genetic diversity of aromatic rice varieties, molecular breeding and chemical and genetic basis of rice aroma. J Plant Biochem Biotechnol. 2010;19(2):127–43.
Das SR. Rice in Odisha. IRRI Technical Bulletin No. 16. Los Baños: International Rice Research Institute; 2012. p. 1–31.
Gangadharan C. Breeding. In: Jaiswal PL, editor. Rice research in India. New Delhi: Indian Council of Agricultural Research; 1985. p. 73–109.
DelCruz N, Khush GS. Rice grain quality evaluation procedures. In: Singh RK, Singh US, Khush GS, editors. Aromatic rices. New Delhi: Oxford & IBH Publishing Co Pvt Ltd; 2000. p. 16–28.
Ramaiah K, Ghose RLM. Origin and distribution of cultivated plants of South Asian rice. Indian J Genet Plant Breed. 1951;11:7–13.
Arunachalam VA, Chaudhury SS, Sarangi SK, Ray T, Mohanty BP, Nambi VA, Mishra S. Rising on rice: The story of Jeypore. Chennai: MS Swaminathan Research Foundation; 2006. p. 1–39.
IBPGR-IRRI. Descriptors for rice (Oryza sativa L.). IBPGR-IRRI Rice Advisory Committee. Manila: IRRI; 1980.
Nei M. Analysis of gene diversity in subdivided populations. Proc Nat Acad Sci USA. 1973;70:3321–3.
Yeh FC, Yang RC, Boyle TBJ. POPGENE, version 1.32 Microsoft window-based freeware for population genetic analysis. Edmonton: University of Alberta; 1999.
Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013;30(12):2725–9.
Peakall R, Smouse PE. Genalex 6: genetic analysis in Excel. Population genetic software for teaching and research. Mol Ecol Notes. 2006;6:288–95.
Weir BS, Cockerham CC. Estimating F-statistics for the analysis of population structure. Evolution. 1984;38:1358–70. doi:10.2307/2408641.
Pritchard JK, Stephens M, Donnelly P. Inference of population structure using multilocus genotype data. Genetics. 2000;155:945–59.
Evanno G, Regnaut S, Goudet J. Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol Ecol. 2005;14:2611–20.
Earl DA, von Holdt BM. STRUCTURE harvester: a website and program for visualizing STRUCTURE output and implementing the Evanno method. Conserv Genet Resour. 2012;4(2):359–61.
Ramasamy RK, Ramasamy S, Bindroo BB, Naik VG. STRUCTURE plot: a program for drawing elegant STRUCTURE bar plots in user friendly interface. Springer Plus. 2014;3(1):1–3. doi:10.1186/2193-1801-3-431.
Zhang P, Li J, Li X, Liu X, Zhao X, Lu Y. Population structure and genetic diversity in a rice core collection (Oryza sativa L.) investigated with SSR markers. PLoS ONE. 2011;6(12):1–13.
SAS Institute. SAS/STAT Version 9.2. Cary: SAS Institute; 2010.
Falush D, Stephens M, Pritchard JK. Inference of population structure using multilocus genotype data: linked loci and correlated allele frequencies. Genetics. 2003;164:1567–87.
Rabbani MA, Pervaiz ZH, Masood MS. Genetic diversity analysis of traditional and improved cultivars of Pakistani rice (Oryza sativa L.) using RAPD markers. Electron. J Biotechnol. 2008;11:1–10.
Parikh M, Motiramani NK, Rastogi NK, Sharma B. Agro-morphological characterization and assessment of variability in aromatic rice germplasm. Bangladesh J Agric Res. 2012;37(1):1–8.
Nascimento WF, Silva EF, Veasey EA. Agro-morphological characterization of upland rice accessions. Sci Agric. 2011;68(6):652–60.
Mathure S, Shaikh A, Renuka N, Wakte K, Jawali N, Thengane R, Nadaf A. Characterisation of aromatic rice (Oryza sativa L.) germplasm and correlation between their agronomic and quality traits. Euphytica. 2011;179:237–46.
Subudhi HN, Swain D, Samantaray S, Singh ON. Collection and agromorphological characterization of aromatic short grain rice in eastern India. Afr J Agric Res. 2012;7(36):5060–8.
Tandekar K, Koshta N. To Study the agro morphological variation and genetic variability in rice germplasm. Middle East J Sci Res. 2014;20(2):218–24.
Sarawgi AK, Subba Rao LV, Parikh M, Sharma B, Ojha GC. Assessment of variability of Rice (Oryza sativa L.) germplasm using agro-morphological characterization. J Rice Res. 2013;6(1):15–28.
Lombardo F, Yoshida H. Interpreting lemma and palea homologies: a point of view from rice floral mutants. Front Plant Sci. 2015;6:1–6.
Vigueira CC, Li W, Olsen KMJ. The role of Bh4 in parallel evolution of hull colour in domesticated and weedy rice. Evol Biol. 2013;26(8):1738–49.
Zhu BF, Si L, Wang Z, Zhou Y, Zhu J, Shangguan Y, et al. Genetic control of a transition from black to straw-white seed hull in rice domestication. Plant Physiol. 2011;155:1301–11.
Ahuja U, Ahuja SC, Chaudhary N, Thakrar R. Red rices—past, present and future. Asian Agri Hist. 2007;11:291–304.
Ray A, Deb D, Ray R, Chattopadhayay B. Phenotypic characters of rice landraces reveal independent lineages of short-grain aromatic indica rice. AoB Plants. 2013;5:1–9. doi:10.1093/aobpla/plt032.
Siddiqui SU, Kumamaru T, Satoh H. Pakistan rice genetic resources. I. Grain morphological diversity and its distribution. Pakistan J Bot. 2007;39:841–8.
Kohli S, Mohapatra T, Das SR, Singh AK, Tandon V, Sharma RP. Composite genetic structure of rice land races revealed by STMS markers. Curr Sci. 2004;86:850–4.
Roy PS, Jena S, Maharana A, Rao GJN, Patnaik SSC. Molecular characterization of short grain aromatic rice landraces of Odisha for detection of aroma. Oryza. 2014;51(2):116–20.
Bisne R, Sarawgi AK. Agro-morphological and quality characterization of Badshahbhog group from aromatic rice germplasm of Chhattisgarh. Bangladesh J Agril Res. 2008;33(3):479–92.
Raghunathachari P, Khanna VK, Singh US. RAPD analysis of genetic variability in Indian scented rice germplasm (Oryza sativa L.). Curr Sci. 2000;79:994–8.
Reddy VS, Dash S, Reddy AR. Anthocyanin pathway in rice (Oryza sativa L.): identification of a mutant showing dominant inhibition of anthocyanins in leaf and accumulation of proanthocyanidins in pericarp. Theor Appl Genet. 1995;91:301–12.
Ramiah K. Rice breeding and genetics. Science monograph. 19. New Delhi: ICAR; 1953.
Varaprasad GS, Shobha Rani N, Padmavati G, Sesu Madhav M, Bentur JS, Laksmi VJ, et al. Catalogue on aromatic short grain rices of India. DRR Technical Bulletin No. 69. Hyderabad: Directorate of Rice Research Rajendranagar; 2013. p. 1–112.
Bajracharya J, Steele KA, Jarvis DI, Sthapit BR, Witcombe JR. Rice landrace diversity in Nepal: variability of agro-morphological traits and SSR markers in landraces from a high-altitude site. Fields Crop Res. 2006;95:327–35.
Dikshit N, Das AB, Sivaraj N, Kar MK. Phenotypic diversity for agro-morphological traits in 105 landraces of rice (Oryza sativa L.) from Santhal Parganas, Jharkhand, India. Proc Natl Acad Sci India Sect B. 2013;83(3):291–304.
Choudhury B, Khan ML, Dayanandan S. Genetic structure and diversity of indigenous rice (Oryza sativa) varieties in the Eastern Himalayan region of Northeast India. Springer Plus. 2013;2(228):1–10. doi:10.1186/2193-1801-2-228.
Tiwari KK, Singh A, Pattnaik S, Sandhu M, Kaur S, Jain S, Tiwari S, et al. Identification of a diverse mini-core panel of Indian rice germplasm based on genotyping using microsatellite markers. Plant Breed. 2015;134(2):164–71.
Kim HJ, Jeong EG, Ahn SN, Doyle J, Singh N, Greenberg AJ, Won YJ, McCouch SR. Nuclear and chloroplast diversity and phenotypic distribution of rice (Oryza sativa L.) germplasm from the democratic people’s republic of Korea. Rice. 2014;7(7):1–15.
Laido G, Mangini G, Taranto F, Gadaleta A, Blanco A, Cattivelli L, et al. Genetic diversity and population structure of tetraploid wheats (Triticum turgidum L.) estimated by SSR, DArT and pedigree data. PLoS ONE. 2013;8(6):1–17. doi:10.1371/journal.pone.0067280.
Yan W, Hu B, Zhang QJ, Jia L, Jackson A, Pan X, Huang B, Yan Z, Deren C. Short and erect rice (ser) mutant from Khao Dawk Mali 105 improves plant architecture. Plant Breed. 2012;131:282–5. doi:10.1111/j.1439-0523.2011.01943.x.
Ookawa T, Yasuda K, Kato H, Sakai M, Seto M, Sunaga K, Motobayashi T, Tojo S, Hirasawa T. Biomass production and lodging resistance in ‘Leaf Star’, a new long-culm rice forage cultivar. Plant Prod Sci. 2010;13:58–66.
Ogunbayo SA, Si M, Ojo DK, Sanni KA, Akinwale MG, Toulou B, Shittu A, Idehen EO, Popoola AR, Daniel IO, Gregorio GB. Genetic variation and heritability of yield and related traits in promising rice genotypes (Oryza sativa L.). J Plant Breed Crop Sci. 2014;6(11):153–9.
Deb D. Seeds of tradition, seeds of future, folk rice varieties of Eastern India. New Delhi: Research Foundation for Science Technology & Ecology; 2005.
Rabara RC, Ferrer MC, Diaz CL, Newingham MCV, Romero GO. Phenotypic diversity of farmers’ traditional rice varieties in the Philippines. Agronomy. 2014;4:217–41.
Yamasaki M, Tenaillon MI, Bi IV, Schroeder SG, Sanchez-Villeda H, Doebley JF, Gaut BS, McMullen MD. A large-scale screen for artificial selection in maize identifies candidate agronomic loci for domestication and crop improvement. Plant Cell. 2005;17:2859–72.
Samal KC, Rout GR, Das SR. Study of genetic divergence of Indigenous Aromatic Rice (Oryza Sativa L.): potentials and consequences of on-farm management in traditional farming. J Agric Sci. 2014;4(4):176–89.
Chase MW, Cowan RS, Hollingsworth PM, vanden Berg C, Madrinan S, Petersen G, et al. A proposal for a standardised protocol to barcode all land plants. Taxon. 2007;56:295–9.
Pasqualone A, Lotti C, Blanco A. Identification of durum wheat cultivars and monovarietal semolinas by analysis of DNA microsatellites. Eur Food Res Technol. 1999;210:144–7.
Ren X, Zhu X, Warndorff M, Bucheli P, Shu Q. DNA extraction and fingerprinting of commercial rice cereal products. Food Res Int. 2006;39:433–9.
Salem HH, Ali BA, Huang TH, Qin DN, Wang XM, Xie QD. Use of random amplified polymorphic DNA analysis for economically important food crops. J Integr PlANT Biol. 2007;49:1670–80.
Terzi V, Morcia C, Gorrini A, Stanca AM, Shewry PR, Faccioli P. DNA-based methods for identification and quantification of small grain cereal mixtures and fingerprinting of varieties. J Cereal Sci. 2005;41:213–20.
De Mattia F, Bruni I, Galimberti A, Cattaneo F, Casiraghi M, Labra M. A comparative study of different DNA barcoding markers for the identification of some members of Lamiacaea. Food Res Int. 2011;44:693–702.
Nicole S, Erickson DL, Ambrosi D, Bellucci E, Lucchin M, Papa R, Kress WJ, Barcaccia G, Donini P. Biodiversity studies in Phaseolus species by DNA barcoding. Genome. 2011;54:529–45.
Rao LVS, Prasad GS, Chiranjivi M, Chaitanya U, Surendhar R. DUS characterization for farmer varieties of rice. IOSR J Agric Vet Sci. 2013;4(5):35–43.
Panda D, Bisoi SS, Palita SK. Floral diversity conservation through sacred groves in Koraput district, Odisha, India: a case study. Int Res J Environ Sci. 2014;3(9):80–6.
Singh AK. Probable agricultural biodiversity heritage sites in India: XVI. The Koraput region. Asian Agri Hist. 2013;17(2):97–122.
Sharma SD, Tripathy S, Biswal J. Origin of O. sativa and its ecotypes. In: Nanda JS, editor. Rice breeding and genetics: research priorities and challenges. New Delhi: Science Publishers, Enfield and Oxford and IBH; 2000. p. 349–69.
Fuller DQ. Agricultural origins and frontiers in South Asia: a working synthesis. J World Prehist. 2006;20:1–86.
Fuller DQ. Finding plant domestication in the Indian subcontinent. Curr Anthropol. 2011;52:S347–62.
Choudhury PR, Kohli S, Srinivasan K, Mohapatra T, Sharma RP. Identification and classification of aromatic rices based on DNA fingerprinting. Euphytica. 2001;118(3):243–51.
The study was designed by PSR, GJNR and TM. AP and SSCP collected the germplasm. RS collected the data. PSR and GJNR analysed the data. NNJ has helped in statistical analysis of data. PSR, GJNR and TM has drafted the manuscript. All authors read and approved the final manuscript.
We sincerely thank, Dr. Sergey Hegay, Senior Scientist, Biotechnology Institute, Bishkek, Kyrgyz Republic for helping in statistical analysis of data. We also acknowledge the help of Mr. Sunil K Sinha, Sr. Technical Officer, Central Rice Research Institute, Cuttack (Odisha), India, in editing the figures.
The authors declare that they have no competing interests.
Aromatic short grain rice genotypes used in the present study and region of their collection. Table S2. Population statistics of the estimated clusters. *SP 1 and SP 2 are estimated subpopulations. Table S3. Pair-wise Nei’s unbiased genetic distance of 19 geographical districts.
Sampling sites in nineteen districts of Odisha (India) and information on population structure. Figure S2. Frequency distribution of phenotypic traits among 126 short grain aromatic rice landraces. Figure S3. The relationship between ΔK and K showing the highest peak at K = 2.
About this article
Cite this article
Roy, P.S., Samal, R., Rao, G.J.N. et al. Differentiation and description of aromatic short grain rice landraces of eastern Indian state of Odisha based on qualitative phenotypic descriptors. BMC Ecol 16, 36 (2016). https://doi.org/10.1186/s12898-016-0086-8
- Aromatic short grain