BACKGROUND
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  <cit> .
since landraces are thought to be the intermediate group between the wild ancestors and cultivated varieties  <cit> , 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 .
aromatic rices constitute a small but an important sub-group  identified by glaszmann  <cit>  with a set of  <dig> isozyme markers.
other than basmati, the indigenous short grain aromatic rices are grown in localized pockets in almost all states of the country  <cit> .
these rices possessing unique aroma, cooking and eating qualities are consumed locally as a delicacy  <cit> .
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  <cit> .
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  <cit> .
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  <cit> .
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  <cit> .
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  <cit> .
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  <cit> .
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  collected from different districts of odisha.
methods
plant material collection and characterization
one hundred and twenty-six short grain aromatic rice lines  were collected by germplasm experts and plant breeders of central rice research institute , cuttack from different parts of odisha, india  during an exploration and collection programme between  <dig> and  <dig> conducted by the institute.
further, the genotypes were grouped under  <dig> geographical districts based on their area of collection .
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  <dig> 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– <dig>  <cit> .
data analysis
data for  <dig> qualitative phenotypic traits were recorded on  <dig> random plants per accession for 3 years and their means were calculated for analysis .
the qualitative phenotypic traits were coded for the presence  and absence  of each alternate form.
diversity parameters viz., number of alleles/variables , effective number of alleles/variables , shannon index  and unbiased nei’s genetic diversity index   <cit>  were calculated using popgene v  <dig>   <cit>  with  <dig> permutations.
the neighbour joining dendrogram based on nei’s unbiased pairwise genetic distances among genotypes was constructed in mega  <dig>  <cit> .
district wise genetic diversity measures were calculated using popgene v  <dig>  and the neighbour joining tree was constructed with pairwise genetic distance measures between districts by mega  <dig>  the molecular variance was analysed by hierarchical amova implemented in the software genealex <dig>  <cit>  based on unbiased fst estimator as described by weir and cockerham  <cit>  to estimate genetic variation among and within districts.
significance of the f-statistics was tested by non-differentiation of probability for  <dig>  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  <cit>  using admix approach with  <dig>  burn-in  periods followed by  <dig>  markov chain monte carlo  replicates with ten independent runs  ranging from  <dig> to  <dig>  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  <cit>  with structure harvester v <dig>   <cit>  and the inferred population clusters were produced by structure plot  <cit>  .
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  <cit> .
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  <dig>  in which the binary data were represented as bar; 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  in the software sas  <dig>   <cit> .
the subset regression was used to identify different principal components that contributed to the total variance in the dataset.table  <dig> phenotypic traits used in the present study for diversity analysis of  <dig> aromatic short grain rice genotypes
basal leaf sheath colour
flag leaf attitude of blade
of stigma
culm angle
of branches
stem thickness
panicle: awns
panicle: colour of awns
branching
panicle: curvature of main axis
spikelet: tip colour of lemma
panicle exertion
leaf intensity of blade colour
leaf: anthocyanin colour
leaf: pubescence of blade surface
leaf: ligule
leaf: colour of ligule
variables only recorded in the present set of germplasm are given.
numbers in parenthesis in the last column represents order of variable for individual trait in the phenoprint of fig.
3b
RESULTS
phenotypic trait diversity
eighteen  out of  <dig> phenotypic traits showed variation among the  <dig> indigenous short grain aromatic rice genotypes, whereas no variation was observed for the traits like intensity of blade colour , leaf anthocyanin colour , pubescence of blade surface , leaf ligule , colour of ligule  and collar colour  .
phenotypic variation for some of the traits are given in fig.
<dig>  the number of variables per individual trait ranged from  <dig>  to a maximum of  <dig> for lemma-palea colour with a total number of  <dig> variables with an average of  <dig>  per trait.
the genetic diversity index  for different traits ranged from  <dig>   to  <dig>   for culm angle and lemma-palea colour, respectively with an average of  <dig> .
the number of effective variables  ranged from  <dig> to  <dig>  with an average of  <dig>  .
further, frequency distribution of the phenotypic traits against the studied genotypes indicated unequal distribution of variables under each trait .
among the  <dig> genotypes, presence of auricles , semi erect culm angle , split ligule , absence of awn , drooping type panicles , semi erect attitude of panicle branches , well exserted panicle , white colour stigma , light gold colour internodes , yellowish white colour of awns , branching of panicles weak , medium stem thickness , green basal leaf sheath colour , lemma having purple tip  and erect flag leaf attitude of blade  were predominant.
a comparatively equal frequency of distribution for short bold  and medium bold  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  followed by purple furrow , brown , purple and brown furrow , black and red  and brown spot .fig.
<dig> variations in aromatic short grain rice genotypes for different phenotypic characters.
a shape of ligule ; b leaf auricle ; c internode colour ; d basal leaf sheath colour ; e tip colour of lemma ; f grain type ; g lemma-palea colour
table  <dig> diversity parameters of phenotypic traits
na
ne
he
i
na number of variables, ne effective number of variables, he nei’s genetic diversity, i shannon diversity index
the principal component analysis  showed that  <dig> of the  <dig> morphological traits were the most important components for explaining the grouping of genotypes .
principal component analysis showed that, first, second, third, fourth, and fifth principal components accounted for  <dig>   <dig>   <dig>   <dig> 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  out of  <dig> morphological traits could adequately explain the total variation existing in these short grain aromatic rice genotypes .fig.
<dig> scree plot of principal component analysis  showing the number of morphological traits and their importance for grouping of aromatic short grain rice landraces
genetic relationship and population structure
the pairwise nei’s genetic distance was calculated to understand the genetic diversity and relatedness among the  <dig> short grain aromatic rice genotypes, which ranged from  <dig>   to  <dig>   with an average of  <dig> .
the neighbour-joining dendrogram constructed based on nei’s unbiased pairwise genetic distance grouped the rice genotypes into  <dig> major clusters , of which, cluster i consisted of  <dig> genotypes while the rest  <dig> 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.fig.
<dig> genetic diversity, population structure and phenoprint representation of  <dig> indigenous aromatic short grain rice genotypes based on  <dig> variables for  <dig> phenotypic traits.
a neighbour-joining dendrogram based on nei’s genetic distance showing genetic relationship among indigenous aromatic short grain rice genotypes; b phenoprint illustration of indigenous aromatic short grain rice genotypes ; c population structure for indigenous aromatic short grain rice genotypes estimated by the structure programme at k = 2; d population structure at k = 4; e names of the corresponding genotypes
bayesian analysis of population structure using the model-based approach with the admixture model proposed by falush and colleagues  <cit>  provided support for the existence of population structure.
the maximum log likelihood  was observed for k =  <dig> .
the genotypes were differentiated into two subpopulations i.e.
sp  <dig> with a membership percentage of  <dig>  % while sp  <dig> with  <dig>  %.
the fixation index  values of subpopulations were  <dig>  and  <dig>  for sp  <dig> and  <dig>  respectively, while the pairwise allele frequency divergence value between the subpopulations was  <dig>  , indicating the existence of moderate population structure in these genotypes.
further,  <dig> and  <dig> genotypes were identified to be admix type having a membership value of <80 % in the subpopulations sp <dig> and sp <dig>  respectively.
moreover, at k =  <dig>  the  <dig> short grain aromatic rice genotypes could be further differentiated into four groups.
.
forty-four genotypes were clustered in group ii followed by  <dig> in group iii,  <dig> in group iv and  <dig> 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.
phenoprint development
the phenoprints of the aromatic short grain rice genotypes were developed based on  <dig> variables of the  <dig> morphological characters .
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- <dig> having purple lines for basal leaf sheath colour, while light purple in dangar basumati, and green in deulabhog- <dig>  the probability of identical match by chance based on all the polymorphic traits was  <dig>  × 10− <dig>  with a differentiation possibility of  <dig> genotypes.
this suggested that the phenoprints can be used as identifying tags for the asgs .
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  were different for the characters like panicle attitude of branches, panicle secondary branching, lemma-palea colour and grain type.
the three accessions of deulabhog  showed variation for a maximum of  <dig> out of  <dig> morphological traits.
further, we performed a heuristic search and identified  <dig> genotypes  i.e.
deulabhog- <dig> , lektimasi , jayaphul , karpurabasa , lektimachi- <dig> , pirima , kalajeera- <dig> , gadakakudinga , garmatia , saragadhuli , ganjeikali , kalazeera , basanadhan , dhurabahila , sunsuniasunaphul  and kalakanhu , which represented the maximum available phenotypic variation  for the  <dig> phenotypic traits recorded.
higher numbers of genotypes  in this group were recorded from koraput and malkangiri region.
the average nei’s genetic diversity value of  <dig>  detected for this group of genotypes was comparable with genetic diversity value of the complete set .
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  <dig> different districts of odisha.
we grouped the genotypes based on these  <dig> districts to address the phenotypic variability within a particular district.
in our collection, maximum  number of genotypes were from koraput followed by cuttack  and ganjam  whereas, only limited number of genotypes  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 .
however, a comparatively uniform but narrow distribution of phenotypic characters was recorded for culm angle  and leaf auricles .
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  as compared to anugul  .
similarly, the highest value for nei’s genetic diversity was recorded in the germplasm set of koraput  followed by puri , anugul having the lowest .
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  <dig> qualitative phenotypic traits, greater part of the variation was observed in koraput district  followed by puri , while, the genotypes from anugul district represented a minimum variation of  <dig>  %.
the pair-wise comparisons of nei’s unbiased genetic distance between the  <dig> districts estimated based on the qualitative morphological traits ranged from  <dig>   to  <dig>   with an average distance value of  <dig>  .
the neighbour-joining tree constructed based on nei’s genetic distance grouped the  <dig> geographical regions into two major clusters .
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  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  based on geographical districts indicated that on an average, 92 % variation exists within districts and 8 % variation was observed among the district .fig.
<dig> distribution of qualitative phenotypic traits in  <dig> geographical districts of collection
table  <dig> diversity parameters of  <dig> aromatic short grain rice genotypes based on  <dig> geographical districts of collection
ne
he
i
% p
ne average effective number of variables, he average nei’s genetic diversity, i average shannon’s information index, % p percent polymorphism
fig.
<dig> unrooted neighbour-joining tree showing genetic relationship among  <dig> geographical districts based on phenotypic of aromatic short grain rice landraces
table  <dig> analysis of molecular variance for  <dig> short grain aromatic rice genotypes from  <dig> geographical districts
f
st
df degree of freedom, ss sum of squares, ms mean squares, ev estimated variance
discussion
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  <cit> .
further, the state of odisha being a major producer of rice in india and having a large assembly of indigenous rice genotypes  <cit> , there was an imperative need to characterize them systematically.
overall phenotypic variation
in this study, we could identify eighteen, out of  <dig> phenotypic traits which showed significant variation in the asg rice collection.
a total of  <dig> variables were detected by  <dig> phenotypic traits in our study, which is significantly higher than that in aromatic rice germplasm of madhya pradesh and chhattisgarh of india  <cit> , in upland rice collection of japan  <cit> , in aromatic rice collection from maharashtra and karnataka of india  <cit> , in aromatic short grain rice genotypes of eastern india  <cit>  and in indigenous rice landraces of chhattisgarh state of india  <cit> .
highest genetic diversity value of  <dig>  was contributed by lemma-palea colour, where  <dig> different forms were identified which is higher than the reports of sarawgi and colleagues  <cit> , who reported  <dig> types in  <dig> rice germplasm collection.
since lemma and palea are sepal equivalent  <cit> , 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  <cit> .
black husk <dig>  gene, encoding an amino acid transporter, fine mapped on chromosome  <dig>  controls the black husk in o. rufipogon and o. nivara.
zhu and colleagues  <cit>  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  <dig> of bh <dig> 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  <cit> .
ray and colleagues  <cit>  reported that aromatic short grain rices vary for their grain characters, which was observed in our study also, where  <dig> different grain types were identified.
similar to our results, high level of variation was reported for grain shape in brazilian and pakistan rice collections  <cit> .
in our study, genotypes having same name collected from different districts  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  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  <cit>  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  in these genotypes having similar name  <cit> .
similarly, bisne and sarawgi  <cit>  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  <cit>  by using rapd markers.
in  <dig> traits , we could not detect any variation for the set of aromatic rice genotypes used in this study.
of the rest  <dig> traits,  <dig> 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  <cit> , 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  <cit> .
the average genetic diversity value recorded in this study was relatively low .
similar type of result was also observed in the all india collection of aromatic rice landraces  <cit> .
however, the genetic diversity value was comparatively high than the rice landraces of nepal  <cit> .
the shannon diversity index of  <dig>  detected in our germplasm set is comparatively high than the rice landraces of santhal paraganas of jharkhand state, india  <cit> .
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.
sp <dig> and sp <dig>  in these natural heritages.
choudhury and colleagues  <cit>  also inferred  <dig> subpopulations in the eastern himalayan and northeast indian indigenous rice with a set of seven microsatellite markers.
similarly,  <dig> subpopulations were detected in indian rice core set with molecular markers, which reveals presence of moderate population structure in the rice genotypes of india  <cit> .
a total of  <dig> genotypes were determined to be admixture type, which might be due to intraspecific crosses between genotypes of sp <dig> and sp <dig>  <cit> .
further, the two identified populations were divided into four different groups similar to the results of laido and colleagues  <cit> , 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  <cit> .
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  <cit> .
in our study,  <dig> genotypes were identified with thick/strong stem, which could be deployed to achieve lodging resistance under higher dose of fertilizers  <cit> .
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  <cit> .
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  <cit> .
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 .
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  <cit> , this technology has been extensively used by several researchers for barcoding of many edible plants .
however, potential application of morphology based characterization cannot be ignored, since it is having useful implication to varietal discrimination at field and farm level  <cit> .
and hence, in this study, we have characterized  <dig> 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  <dig> 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  <cit> .
further, koraput region has been recently declared as ‘globally important agricultural heritage systems’ by the food and agricultural organization of the united nations   <cit> .
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  <cit> , 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  <cit>  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.
CONCLUSIONS
native aromatic rice landraces that are highly preferred by consumers needs to be characterized which can help in varietal diagnostics purpose and their conservation  <cit> .
here, we have characterised native short grain aromatic genotypes of odisha , 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  <dig> different forms were detected followed by tip colour of lemma where  <dig> 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  <dig> genotypes were broadly grouped into  <dig> sub populations and  <dig> 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.
additional files
<dig> /s12898-016-0086- <dig> aromatic short grain rice genotypes used in the present study and region of their collection.
table s <dig>  population statistics of the estimated clusters.
*sp  <dig> and sp  <dig> are estimated subpopulations.
table s <dig>  pair-wise nei’s unbiased genetic distance of  <dig> geographical districts.
<dig> /s12898-016-0086- <dig> sampling sites in nineteen districts of odisha  and information on population structure.
figure s <dig>  frequency distribution of phenotypic traits among  <dig> short grain aromatic rice landraces.
figure s <dig>  the relationship between Δk and k showing the highest peak at k =  <dig>
authors’ contributions
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.

 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. aromatic rices constitute a small but an important sub-group identified by glaszmann with a set of isozyme markers. other than basmati, the indigenous short grain aromatic rices are grown in localized pockets in almost all states of the country. these rices possessing unique aroma, cooking and eating qualities are consumed locally as a delicacy. 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. 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. 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 collected from different districts of odisha. methods plant material collection and characterization one hundred and twenty-six short grain aromatic rice lines were collected by germplasm experts and plant breeders of central rice research institute, cuttack from different parts of odisha, india during an exploration and collection programme between and conducted by the institute. further, the genotypes were grouped under geographical districts based on their area of collection. 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 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–. data analysis data for qualitative phenotypic traits were recorded on random plants per accession for 3 years and their means were calculated for analysis. the qualitative phenotypic traits were coded for the presence and absence of each alternate form. diversity parameters viz., number of alleles/variables, effective number of alleles/variables, shannon index and unbiased nei’s genetic diversity index were calculated using popgene v with permutations. the neighbour joining dendrogram based on nei’s unbiased pairwise genetic distances among genotypes was constructed in mega. district wise genetic diversity measures were calculated using popgene v and the neighbour joining tree was constructed with pairwise genetic distance measures between districts by mega the molecular variance was analysed by hierarchical amova implemented in the software genealex 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 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 burn-in periods followed by markov chain monte carlo replicates with ten independent runs ranging from to 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 v and the inferred population clusters were produced by structure plot. 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 in which the binary data were represented as bar; 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 in the software sas. the subset regression was used to identify different principal components that contributed to the total variance in the dataset.table  phenotypic traits used in the present study for diversity analysis of aromatic short grain rice genotypes basal leaf sheath colour flag leaf attitude of blade of stigma culm angle of branches stem thickness panicle: awns panicle: colour of awns branching panicle: curvature of main axis spikelet: tip colour of lemma panicle exertion leaf intensity of blade colour leaf: anthocyanin colour leaf: pubescence of blade surface leaf: ligule leaf: colour of ligule variables only recorded in the present set of germplasm are given. numbers in parenthesis in the last column represents order of variable for individual trait in the phenoprint of fig. 3b  phenotypic trait diversity eighteen out of phenotypic traits showed variation among the indigenous short grain aromatic rice genotypes, whereas no variation was observed for the traits like intensity of blade colour, leaf anthocyanin colour, pubescence of blade surface, leaf ligule, colour of ligule and collar colour. phenotypic variation for some of the traits are given in fig.  the number of variables per individual trait ranged from to a maximum of for lemma-palea colour with a total number of variables with an average of per trait. the genetic diversity index for different traits ranged from to for culm angle and lemma-palea colour, respectively with an average of. the number of effective variables ranged from to with an average of. further, frequency distribution of the phenotypic traits against the studied genotypes indicated unequal distribution of variables under each trait. among the genotypes, presence of auricles, semi erect culm angle, split ligule, absence of awn, drooping type panicles, semi erect attitude of panicle branches, well exserted panicle, white colour stigma, light gold colour internodes, yellowish white colour of awns, branching of panicles weak, medium stem thickness, green basal leaf sheath colour, lemma having purple tip and erect flag leaf attitude of blade were predominant. a comparatively equal frequency of distribution for short bold and medium bold 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 followed by purple furrow, brown, purple and brown furrow, black and red and brown spot.fig.  variations in aromatic short grain rice genotypes for different phenotypic characters. a shape of ligule; b leaf auricle; c internode colour; d basal leaf sheath colour; e tip colour of lemma; f grain type; g lemma-palea colour table  diversity parameters of phenotypic traits na ne he i na number of variables, ne effective number of variables, he nei’s genetic diversity, i shannon diversity index the principal component analysis showed that of the morphological traits were the most important components for explaining the grouping of genotypes. principal component analysis showed that, first, second, third, fourth, and fifth principal components accounted for 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 out of morphological traits could adequately explain the total variation existing in these short grain aromatic rice genotypes.fig.  scree plot of principal component analysis showing the number of morphological traits and their importance for grouping of aromatic short grain rice landraces genetic relationship and population structure the pairwise nei’s genetic distance was calculated to understand the genetic diversity and relatedness among the short grain aromatic rice genotypes, which ranged from to with an average of. the neighbour-joining dendrogram constructed based on nei’s unbiased pairwise genetic distance grouped the rice genotypes into major clusters, of which, cluster i consisted of genotypes while the rest 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.fig.  genetic diversity, population structure and phenoprint representation of indigenous aromatic short grain rice genotypes based on variables for phenotypic traits. a neighbour-joining dendrogram based on nei’s genetic distance showing genetic relationship among indigenous aromatic short grain rice genotypes; b phenoprint illustration of indigenous aromatic short grain rice genotypes; c population structure for indigenous aromatic short grain rice genotypes estimated by the structure programme at k = 2; d population structure at k = 4; e names of the corresponding genotypes 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 was observed for k = . the genotypes were differentiated into two subpopulations i.e. sp with a membership percentage of  % while sp with  %. the fixation index values of subpopulations were and for sp and respectively, while the pairwise allele frequency divergence value between the subpopulations was, indicating the existence of moderate population structure in these genotypes. further, and genotypes were identified to be admix type having a membership value of <80 % in the subpopulations sp and sp respectively. moreover, at k =  the short grain aromatic rice genotypes could be further differentiated into four groups.. forty-four genotypes were clustered in group ii followed by in group iii, in group iv and 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. phenoprint development the phenoprints of the aromatic short grain rice genotypes were developed based on variables of the morphological characters. 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- having purple lines for basal leaf sheath colour, while light purple in dangar basumati, and green in deulabhog- the probability of identical match by chance based on all the polymorphic traits was  × 10− with a differentiation possibility of genotypes. this suggested that the phenoprints can be used as identifying tags for the asgs. 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 were different for the characters like panicle attitude of branches, panicle secondary branching, lemma-palea colour and grain type. the three accessions of deulabhog showed variation for a maximum of out of morphological traits. further, we performed a heuristic search and identified genotypes i.e. deulabhog-, lektimasi, jayaphul, karpurabasa, lektimachi-, pirima, kalajeera-, gadakakudinga, garmatia, saragadhuli, ganjeikali, kalazeera, basanadhan, dhurabahila, sunsuniasunaphul and kalakanhu, which represented the maximum available phenotypic variation for the phenotypic traits recorded. higher numbers of genotypes in this group were recorded from koraput and malkangiri region. the average nei’s genetic diversity value of detected for this group of genotypes was comparable with genetic diversity value of the complete set. 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 different districts of odisha. we grouped the genotypes based on these districts to address the phenotypic variability within a particular district. in our collection, maximum number of genotypes were from koraput followed by cuttack and ganjam whereas, only limited number of genotypes 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. however, a comparatively uniform but narrow distribution of phenotypic characters was recorded for culm angle and leaf auricles. 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 as compared to anugul. similarly, the highest value for nei’s genetic diversity was recorded in the germplasm set of koraput followed by puri, anugul having the lowest. 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 qualitative phenotypic traits, greater part of the variation was observed in koraput district followed by puri, while, the genotypes from anugul district represented a minimum variation of  %. the pair-wise comparisons of nei’s unbiased genetic distance between the districts estimated based on the qualitative morphological traits ranged from to with an average distance value of. the neighbour-joining tree constructed based on nei’s genetic distance grouped the geographical regions into two major clusters. 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 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 based on geographical districts indicated that on an average, 92 % variation exists within districts and 8 % variation was observed among the district.fig.  distribution of qualitative phenotypic traits in geographical districts of collection table  diversity parameters of aromatic short grain rice genotypes based on geographical districts of collection ne he i % p ne average effective number of variables, he average nei’s genetic diversity, i average shannon’s information index, % p percent polymorphism fig.  unrooted neighbour-joining tree showing genetic relationship among geographical districts based on phenotypic of aromatic short grain rice landraces table  analysis of molecular variance for short grain aromatic rice genotypes from geographical districts f st df degree of freedom, ss sum of squares, ms mean squares, ev estimated variance discussion 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 phenotypic traits which showed significant variation in the asg rice collection. a total of variables were detected by 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 was contributed by lemma-palea colour, where different forms were identified which is higher than the reports of sarawgi and colleagues, who reported types in 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 husk gene, encoding an amino acid transporter, fine mapped on chromosome 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 of bh 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 different grain types were identified. similar to our results, high level of variation was reported for grain shape in brazilian and pakistan rice collections. in our study, genotypes having same name collected from different districts 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 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 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 traits, we could not detect any variation for the set of aromatic rice genotypes used in this study. of the rest traits, 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, 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. 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 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. sp and sp in these natural heritages. choudhury and colleagues also inferred subpopulations in the eastern himalayan and northeast indian indigenous rice with a set of seven microsatellite markers. similarly, 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 genotypes were determined to be admixture type, which might be due to intraspecific crosses between genotypes of sp and sp. 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, 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. 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. 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 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 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. 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 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, 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 different forms were detected followed by tip colour of lemma where 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 genotypes were broadly grouped into sub populations and 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. additional files /s12898-016-0086- aromatic short grain rice genotypes used in the present study and region of their collection. table s population statistics of the estimated clusters. *sp and sp are estimated subpopulations. table s pair-wise nei’s unbiased genetic distance of geographical districts. /s12898-016-0086- sampling sites in nineteen districts of odisha and information on population structure. figure s frequency distribution of phenotypic traits among short grain aromatic rice landraces. figure s the relationship between Δk and k showing the highest peak at k = authors’ contributions 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. 
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. this study reports qualitative phenotypic trait based characterization of short grain aromatic rice genotypes, collected from different areas of the state of odisha, india. out of the descriptors employed, highest variability was observed for lemma-palea colour with a genetic diversity index of. 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. speciality rice, in general, and aromatic rice in particular, possess enormous market potential for enhancing farm profits. though aromatic rices represent a tiny fraction of the total rice germplasm, a small collection of land races did exhibit rich diversity for all the qualitative traits. 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 and respectively. the proposed system not only help in conservation but also can confer ipr protection to these specialty rices.
 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 short grain aromatic rice genotypes, collected from different areas of the state of odisha, india.  out of the descriptors employed, highest variability was observed for lemma-palea colour with a genetic diversity index of. 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 and 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 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 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. electronic supplementary material the online version of this article contains supplementary material, which is available to authorized users. keywords ricelandracesaromatic short grainphenotypictraitcharacterizationdiversityissue-copyright-statement© the author 2016 