Genetic diversity and structure in hill rice (Oryza sativa L.) landraces from the North-Eastern Himalayas of India
© The Author(s). 2016
Received: 9 February 2016
Accepted: 7 July 2016
Published: 13 July 2016
Hill rices (Oryza sativa L.) are direct seeded rices grown on hill slopes of different gradients. These landraces have evolved under rainfed and harsh environmental conditions and may possess genes governing adaptation traits such as tolerance to cold and moisture stress. In this study, 64 hill rice landraces were collected from the state of Arunachal Pradesh of North-Eastern region of India, and assessed by agro-morphological variability and microsatellite markers polymorphism. Our aim was to use phenotypic and genetic diversity data to understand the basis of farmers’ classification of hill rice landraces into two groups: umte and tening. Another goal was to understand the genetic differentiation of hill rices into Indica or japonica subspecies.
According to farmers’ classification, hill rices were categorized into two groups: umte (large-grained, late maturing) and tening (small-grained, early maturing). We did not find significant difference in days to 50 % flowering between the groups. Principal component analysis revealed that two groups can be distinguished on the basis of kernel length-to-width ration (KLW), kernel length (KL), grain length (GrL), grain length-to-width ration (GrLW) and plant height (Ht). Stepwise canonical discriminant analysis identified KL and Ht as the main discriminatory characters between the cultivar groups. Genetic diversity analysis with 35 SSR markers revealed considerable genetic diversity in the hill rice germplasm (gene diversity: 0.66; polymorphism information content: 0.62). Pair-wise allelic difference between umte and tening groups was not statistically significant. The model-based population structure analysis showed that the hill rices were clustered into two broad groups corresponding to Indica and Japonica. The geographic distribution and cultivars grouping of hill rices were not congruent in genetic clusters. Both distance- and model-based approaches indicated that the hill rices were predominantly japonica or admixture among the groups within the subspecies. These findings were further supported by combined analysis hill rices with 150 reference rice accessions representing major genetic groups of rice.
This study collected a valuable set of hill rice germplasm for rice breeding and for evolutionary studies. It also generated a new set of information on genetic and phenotypic diversity of hill rice landraces in North-Eastern region of India. The collected hill rices were mostly japonica or admixture among the subpopulations of Indica or Japonica. The findings are useful for utilization and conservation of hill rice germplasm.
Rice is the primary food source for more than half of the world’s population . Over 90 % of the world’s rice is produced and consumed in the Asia-Pacific region. Asian cultivated rice (Oryza sativa L.) is one of the important crops in the world. Being a model organism with fully sequenced genome, it also affords unique opportunities to use genomic approaches to study its domestication, adaptive selection, and the history of crop improvement [2, 3]. The genetic structure of rice on a global scale is well characterized (see reviews in [4–6]). In addition to the two widely accepted major subspecies, Indica and Japonica, three to seven genetically distinct groups have been identified within global rice (Oryza sativa) germplasm sets in different studies [7–14]. Numerous studies on genetic structure of rice cultivars at a local scale (within a country) have also been conducted [15–22]. Such local-scale studies not only provide a detail view of rice genetic diversity within a country, but also allow for a better understanding of complex interaction between rice genetic diversity and human cultivation practices , and for formulating in situ conservation strategies . The cultivated rice genetic diversity includes indigenous landraces and commercially bred improved varieties. As a consequence of both natural and artificial selection under diverse habitats, high levels of morphological, physiological and genetic diversity exist within O. sativa, resulting in >120,000 distinct rice varieties [1, 23]. An assessment of this diversity is important for making decision about the conservation and effective utilization in breeding programmes.
The North-Eastern (NE) region of India, due to its unique eco-geographical features and high ethnic diversity, possesses at least 10,000 indigenous rice landraces which cultivated under upland, lowland and deep water conditions . The impact of Green Revolution on the varietal landscape of rice in this region was insignificant, particularly in the hilly states, and the farmers still grow their heirloom varieties which not only suit their taste but also provide crop security. The current study estimates the phenotypic and genetic diversity present in a set of hill rice landraces collected from temperate to sub-temperate mountain region of the state of Arunachal Pradesh. The state of Arunachal Pradesh is situated on the eastern most corner of India sharing international borders with Bhutan, China and Myanmar. It has a geographical area of 83,743 sq. km. The state has a population density of 17 persons per sq. km, and around 68.8 % of the population is tribal . The tribal communities of Arunachal Pradesh are of Tibeto-Burman linguistic origin, and there are as many as 21 tribes and 50 sub-tribes. The topography of the state is mostly hilly terrain and the agro-climate varies from tropical to alpine. The farmers follow traditional shifting or swidden cultivation (locally known as jhum) and sedentary agriculture. Shifting or ‘slash-and-burn’ cultivation is the earliest form of agriculture and is still practiced in vast areas of the state. About 76 % of the total cropped area in the state is under jhum cultivation.
Rice is the principal crop of Arunachal Pradesh and it is grown in 44.7 % of the total cropped land and accounts for 70.3 % of the total food grain production . Rice productivity in the state is very low (2065 kg/ha). Here rice is grown up to an altitude of 2000 m as a rainfed crop in hill slopes of newly cleared forests (jhum fields), plateaus, terraces and river beds. A wide range of rice landraces were found throughout the state. Sometimes it depends on the tribal communities who inhabit the area. It was estimated that more than 70 % of the total rice varieties grown in the state are traditional cultivars or landraces having low to very low productivity . Despite having low yield potential, the hill rice landraces grown in the mountains under jhum cultivation system possess many adaptive traits such as cold tolerance important for crop improvement programmes. Cold-tolerance in hill rice landraces was also observed in high altitude (1300–2200 m) rice areas in Nepal . Similarly, in the high altitude rice areas of Arunachal Pradesh cold injury due to cool air temperature, particularly during flowering stage, is a major constraint to improving rice productivity. In general, hill rice landraces are grown on the cleared mountain sides for one year and then left fallow for 4–5 years before the farmers return to the same location (Additional file 1). The farmers usually grow few maize and millets along with rice in the same field. The hill rice landraces have been evolved by the interplay between adaptation to the harsh environment and selection imposed by the farmers who determine which varieties will be grown under a particular agro-ecological condition. Therefore, it will be interesting to study the impact of farmers’ practices and local culture on the conservation, exchange and genetic structure of traditional landrace varieties, by studying rice landraces collected from a unique ecological niche like mountain sides. There is no previous study done on genetic diversity of the hill (jhum) rice cultivars of Arunachal Pradesh.
In 2012, an exploration trip to three districts viz, East Kameng, Papum Pare and Kurung Kumey of the state of Arunachal Pradesh was conducted for the collection of hill rice landraces. These rices are direct seeded, grown on upland fields, usually on sloping hill sides (jhum lands), and are grown as rainfed during June to September. During the collection trip, we interviewed local farmers to gather information on significance of naming a landrace, its cultural value and popularity, and production practices of local farmers. The indigenous farmers of the surveyed region rely on gathering forest products, hunting and cultivation of rice and maize for their subsistence. In the present study, we phenotyped these rice landraces for 16 agro-morphological and grain characteristics, and genotyped them with 35 SSR markers. The specific objectives of this study were to: (i) collect information on the meanings of the variety names, special uses or properties, and important agronomic characteristics; (ii) study the agro-morphological variability; and (iii) understand the genetic relatedness and structure of these landraces.
Collection of hill rice
Out of these 67 hill rice accessions, altogether 64 were selected for agro-morphological and genetic characterization in this study. Information about the accessions such as name, accession number, place of collection/seed sources, farmers’ grouping of landraces and subpopulation ancestry based on STRUCTURE is given in Additional file 3. The seeds of the hill rice accessions used in this study are conserved in the National Genebank of ICAR-National Bureau of Plant Genetic Resources (NBPGR), New Delhi, India, and are publicly available for research purposes upon request with a material transfer agreement. The rice accessions were grown in 2013 and 2014 for purification. In each year, seeds from a single panicle were kept for sowing in the next year.
Three check varieties such as ‘Bhalum-1’, ‘Bhalum-3’ and ‘Sukhadhan-1’ were used during phenotypic characterization trial. ‘Bhalum-1’ and ‘Bhalum-3’ varieties are suited for cultivation under rainfed upland conditions up to an altitude of 1000 m. ‘Sukhadhan-1’ is a drought tolerant upland rice variety of Nepal. Seed of these accessions were obtained from ICAR Research Complex for North-Eastern Hill Region, Umiam, India.
A total of 15 control varieties were selected for SSR genotyping along with the hill rices (Additional file 3). These checks included five indica (‘IR8’, ‘IR36’, ‘Jaya’, ‘Taichung Native 1’ and ‘TKM9’), two aus (‘FR 13A’ and ‘Kalamkati’), three Basmati (‘Basmati370’, ‘Pakistani Basmati’ and ‘Pusa Basmati 1’), three short-grain traditional aromatic (‘Badshabhog’, ‘Kataribhog’ and ‘Mohanbhog’) and a single accession each of tropical japonica (‘Azucena’) and japonica admixture (‘Dular’). Pure seeds of these accessions were obtained from ICAR-National Rice Research Institute, Cuttack, India. These varieties were used to understand genetic grouping of the hill rice accessions, and also to act as a molecular weight reference in DNA gel electrophoresis. We did not include these varieties in phenotypic evaluation as the plant growth was very poor for most of the accessions.
To provide a reference for the rice subpopulation classification in the current study, we also used the SSR fingerprinting data of 150 rice accessions (here we denoted these as reference accessions) previously reported in the genetic diversity study by Garris et al. . The band size data were obtained from ‘Gramene genetic diversity database’ (http://archive.gramene.org/db/diversity/diversity_view). We sampled 150 accessions from the originally reported 234 accessions used in Garris et al.  based on the recent report on open access resources for genome-wide association mapping in rice .
Range of variation in agro-morphological characteristics of hill rice landraces and the loading of the measured traits on five principal components
Mean ± SD
115.78 ± 13.05
120.84 ± 11.79
81.11 ± 23.52
69.86 ± 17.90
53.76 ± 11.04
59.23 ± 10.87
1.86 ± 0.21
1.92 ± 0.25
1.65 ± 0.53
1.56 ± 0.38
5.19 ± 0.97
5.39 ± 0.97
23.17 ± 2.87
22.16 ± 2.65
143.07 ± 52.53
159.89 ± 57.13
6.94 ± 0.70
6.19 ± 0.65
3.39 ± 0.37
3.46 ± 0.37
2.07 ± 0.31
1.81 ± 0.34
2.31 ± 0.50
2.19 ± 0.30
23.03 ± 4.32
22.85 ± 3.09
5.44 ± 0.46
4.58 ± 0.42
2.9 ± 0.35
3.02 ± 0.24
1.91 ± 0.34
1.53 ± 0.25
Genotyping with SSR markers
Thirty five SSR markers were designed from the ‘Gramene’ marker database (http://www.gramene.org/markers/microsat/). Information about these markers is given in Additional file 4. Extraction of total genomic DNA, polymerase chain reaction (PCR) with SSR markers and molecular size determination of amplified band were performed as described in Roy et al. . During gel electrophoresis at least one control accession of ‘IR8’, ‘IR36’, ‘Taichung Native 1’, ‘FR 13A’, ‘Kalamkati’ or ‘Dular’ was used as a molecular weight reference in each gel because a reference allele size in these control accessions is available . The original SSR genotype data of all rice accessions used in this study are provided in Additional file 5.
Agro-morphological data analysis
Summary statistics such as mean, standard deviation (SD), minimum and maximum values, and coefficient of variation (CV) were determined. The mean data were standardized by deriving Z scores for further analyses. ANOVA was carried out to test the significance of variation between hill rice cultivar groups for different traits. Principal component analysis (PCA) was performed on the correlation matrix of the data to understand the most important variables contributing to the total phenotypic variation among the accessions. PCA represents the simplest and most commonly used multivariate method to visualize the grouping of accessions based on component loadings. Ward’s hierarchical clustering was used to understand the relationships among rice accessions based on phenotypic data. All these analyses were performed using IBM-SPSS Statistics version 20.0 .
Population structure analysis
Population structure of 64 hill and 15 control rice accessions was examined using the Bayesian model-based approach implemented in STRUCTURE V2.3.4 . The number of clusters (K) evaluated here ranged from 1 to 8. The analysis was performed using five replicate runs per K value, a burn-in period length of 5000, a run length of 50,000, and a model allowing for admixture and correlated allele frequency. ‘Structure harvester’ programme (http://taylor0.biology.ucla.edu) was used to determine the final K value(s) based on both the LnP(D) and Evanno’s ΔK . Subsequently, ten simulations at K = 2–4 were then performed with a burn-in period of 10,000 and a run length of 100,000. The membership coefficient from the run with the lowest likelihood value was used to assign each accession to the K = 1 to 4 subpopulations based on the estimated membership coefficients.
Summary statistics and genetic diversity analysis
The average number of alleles per locus (AN), major allele frequency (MAF), gene diversity (H e ), heterozygosity (H o ) and polymorphism information content (PIC) were calculated using PowerMarker V3.25 . Average allelic richness (R s ) and Wright’s fixation index (F st ) values were calculated using FSTAT V22.214.171.124 . The molecular variance of subpopulations and accessions within the subpopulations were calculated using an Analysis of Molecular Variance (AMOVA) approach in GenAlEx V6.5 . Separate analyses were conducted by classifying the rice accessions into: districts, farmers’ classified groups and STRUCTURE subpopulations. The statistical test of the differences in allelic count among different groups with variable numbers of accessions was conducted in R  using FPTestR, the R version of FPTest method reported in Fu et al. . The FPTestR package (communicated for publication) was kindly provided by Dr Yong-Bi Fu, Plant Genetic Resources of Canada, Agriculture and Agri-Food Canada. The allele frequency data from PowerMarker was used to export the data in binary format (1/0) for analysis with NTSYS-pc V2.2 . A neighbour-joining (NJ) cluster diagram was constructed with the NJOIN sub-programme using the genetic dissimilarity matrix calculated in SIMGEND sub-programme with Nei72 coefficient. To summarize the patterns of variation in multi-locus dataset, principal coordinate analysis (PCoA) was performed in GenAlex software using the genetic distance matrix among the accessions.
Combined analysis of hill and reference accessions
Cluster analysis was conducted on hill and reference accessions to understand the grouping of hill rices in respect to different genetically defined groups (indica, aus, aromatic, tropical japonica and temperate japonica). The reference set also included some ‘admixture’ accessions. For clustering, a pair-wise genetic distance matrix was calculated in PowerMarker following C. S. Chord distance method . The NJ method was used for phylogenetic reconstruction. The unrooted NJ tree was visualized using Dendroscope V3 . Pair-wise F st and allelic differences (FPTest) between the groups were calculated as stated earlier.
Variability in agro-morphological traits
Based on farmers’ classification of the hill rices, the collected rice germplasm was categorized into two cultivar types: (i) umte (um = rice; te = large), large grained-late maturing rice, and (ii) tening, small grained-early maturing types (Additional file 3). The umte cultivars are high yielding and good to taste. In this study, we have assessed the agro-morphological and genetic diversity between these two groups.
Descriptive statistics for 16 agro-morphological characteristics in umte and tening groups of hill rice are given in Table 1. The variation observed between two groups was significant for the traits evaluated such as plant height (Ht), leaf length (LL), grain length (GrL), grain length-to-width ratio (GrLW), kernel length (KL) and kernel length-to-width ratio (KLW). There was no significant difference in days to 50 % flowering (DtF) between gropus. The range of variation in DtF in both the groups was similar, and interestingly tening group showed higher average value of DtF (120.8 d) than umte (115.8 d). However, the difference was prominent for grain and kernel characteristics. The umte accessions had longer grains with an average value of 6.9 mm than that of tening (6.2 mm). Similarly, the value of KL in umte (5.4 mm) was higher than that recorded in tening (4.6 mm). Both the highest values for GrLW and KLW were recorded in umte. On average, the plants of umte accessions (81.1 cm) were taller than the tenings (69.9 cm). The tenning accessions produced larger leaves (59.7 cm).
Ward’s hierarchical cluster analysis revealed a similar grouping of accessions as found in PCA (Fig. 2b). Overall, 67 accessions were divided into two groups. All check varieties formed a distinct cluster (C1) along with ten hill accessions. The C2 was further divided into two sub-clusters. There was no clear grouping of umte and tening accessions, except a sub-group (C2B) of 15 tenning accessions which had highest values for KW, GrLW, DtF, LL and leaf width (LW).
The summary statistics of 35 SSR markers are given in Additional file 4. Overall, 297 alleles were detected at the 35 SSR markers, ranging from 2 alleles (RM338) up to 21 alleles (RM259), with an average of 8.49 alleles across the loci. The gene diversity or expected heterozygosity (H e ) ranged from 0.41 (RM55) to 0.94 (RM259) and the average observed heterozygosity (H o ) was 0.051. Allelic richness (R s ) varied from 2.0 (RM338) to 14.6 (RM259) with an average value of 6.04 across the study. The polymorphism information content (PIC) values ranged from 0.37 (RM338 and RM507) up to 0.93 (RM259) with an average of 0.65.
Genetic structure and subdivision of hill rice accessions
Genetic diversity statistics of hill rice landraces at different sub-population levels
Farmers’ classified groups
sub-populations (K = 4)
All hill accessions
The number of alleles in a sample (allelic count or allelic richness) is an important measure of diversity. Pair-wise tests of allelic differences between different groups of accessions are given in Fig. 3b. Allelic counts of Papum Pare and Kurung Kumey was significantly different with that of the control. The allelic differences between the districts were also significant, except the case between East Kameng and Kurung Kumey. No significant allelic differences were noted between umte and tening, or between controls and either of the two groups. Among the four STRUCTURE subpopulations, allelic count of StrGr3 was significantly different with that of StrGr1 and StrGr2 (Fig. 3b). Pair-wise allelic difference between StrGr3 and StrGr4 was also significant.
Pair-wise estimates of F st values among and between subpopulations are given in Additional file 8. The level of differentiation between umte and tening was very low (0.008). Overall, the rice accessions from three districts had a low level of differentiation with F st values ranging from 0.043 to 0.126. Among the STRUCTURE groups, the greatest differentiation was observed between StrGr1 and 2 (0.290), and the lowest was between StrGr3 and 4 (0.086).
Genetic grouping in relation to reference rice accessions
Our characterization of hill rice germplasm of North-Eastern region of India has revealed that a diverse array of rice landraces was grown by the farmers which traditionally classified as umte and tenning based on grain morphology and crop duration. These two groups of hill rice germplasm varied significantly for the traits such as plant height, leaf length, grain/kernel length and grain/kernel length-to-width ratio, but non-significant genetic differentiation between the two groups were observed. Population genetic analysis showed that these hill rices were mostly japonica or admixture among the subpopulations of Indica or Japonica. These results are useful for understanding the agro-morphological variability, genetic diversity and structure of hill rice landraces, and significant for utilization of hill rices in crop improvement programmes. Our study also provides baseline information for planning in-situ conservation of the hill rice landraces, which are adapted to local environments in harmony with cultural preferences of the farmers.
Hill rice landrace
The traditional names of rice cultivars are generally found to be a valuable concept in describing the landraces worldwide. In the present study, the names of rice cultivars varied among the districts. Few rice accessions ‘Taba’, ‘Kilung’ and ‘Sarpung’ were cultivated in all three districts. The farmers in Arunachal Pradesh most often name the rice landraces after their place of origin. For instance, ‘Nepali dhan’ is originated in Nepal. Likewise, ‘Manipuri dhan’ and ‘Naga dhan’ are from the state of Manipur and Nagaland, respectively. Rice landraces often named after a community/tribe. For example, ‘Taba dugu’ and ‘Adi’ named after ‘Taba’ and ‘Adi’ communities of Arunachal Pradesh, respectively. During collection trip the farmers expressed that hill rice cultivars generally taste better than the low land cultivars, and thus used for making local delicacies like pitha and sweets. Therefore, knowledge about the local names of rice landraces is useful for understanding their origin as well as grain quality and other characteristics. The traditional rice landraces possess wide genetic variability valuable for rice breeding [42, 43]. The broad genetic base of traditional rice cultivars is suitable for subsistence farming practiced by the farming communities. Rice germplasm collected from geographically isolated and ecologically important sites need to be characterized both at phenotypic and genetic levels to understand the interaction and genetic controls for traits influencing adaptation of rice crop to harsh environments, and to identify sources for important adaptive traits . In Nepal, the rice accessions sampled from mid- and high-hill sites were found to possess tolerance to cold stress. However, interestingly rice diversity in the high-hill regions was lower than those of mid-hill and low-hill sites . It was proposed that rice diversity in the mid-hills is linked with the range in altitude in these regions that results in great environmental heterogeneity and diverse agro-ecosystems, and great diversity in the socio-economic structure of the farming communities .
The tribal farmers of Arunachal Pradesh usually classify hill rice cultivars on the basis of grain morphology and/or quality characteristics, in addition to other agronomic traits such as crop duration. In the present study, we have assessed the variability in agronomic and grain/ kernel morphological traits to understand the morphological factors differentiating the two groups (umte and tening). Our analysis apparently showed that two groups of hill rice are mainly differentiated by the traits such as grain/kernel length and length-to-width ratio. We did not find significant variation in days to flowering between umte and tening. However, this is a key trait considered by the farmers in categorizing the hill rices. The non-significant variation in days to flowering in the present study may be due to change in crop growing conditions including soil and other environmental factors. We also performed a step-wise canonical discriminant analysis (CDA) taking umte and tening as priori groups to identify the most important discriminating traits contributing to the grouping of hill rices. Kernel length and plant height were found to be the best predictors for classifying the cultivar groups (results not shown). In both PCA and cluster analysis, we found that the accessions with similar names formed close groups, indicating the consistency in naming the cultivars by traditional farmers. However, in a different study with 130 indigenous rice accessions from a neighbouring state (Nagaland) in NE India, we found inconsistency in naming the cultivars . Similar observation was also reported in upland rices of Nepal . The grouping of hill rice accessions in both PCA and Ward’s cluster did not follow geographical origin for a genotype. Similar trend was also recorded in previous studies [47, 49]. The average number of panicles per plant in hill rice accessions of Arunachal Pradesh (5.30) is lower than those reported in rice landraces from Nagaland state (7.9)  and rice cultivars adapted to high-altitude environments in Nepal (8.0) . Extremely strong selection pressure for adaptation to the harsh high-altitude environments and direct seeding method of planting could have limited the number of productive tillers in hill rice germplasm.
Genetic diversity and structure
In the present study, we have used 35 SSR markers for genotyping rice landraces. A majority of these markers were chosen from the panel of 50 standard SSR markers recommended for rice diversity analysis by the Generation Challenge Programme by CGIAR. Among these SSRs, 33 markers also represent a subset of 169 markers originally used by Garris et al. . Grouping 150 reference rice accessions from Garris et al.  using these 33 SSRs, was nearly the same. This suggests that the 35 SSR markers used in the present study can sufficiently resolve the genetic structure in rice.
Generally, the rice varieties grown under upland culture in the hill areas of Southeast Asia as a component of shifting cultivation are japonicas . Earlier studies also indicated that a majority of upland rice landraces grown in upland fields are japonica . However, upland indica cultivars were also recorded in different studies [20, 50]. An earlier study targeting indica-japonica classification of Asian rice ecotypes and Japanese rice cultivars proposed that lowland rices can be clearly classified into indica or japonica, while, upland cultivars cannot . In the present study, we attempted to classify the hill rice accessions by grouping them in relation to the control varieties. The STRUCTURE analysis revealed that considering the optimal grouping of the genotypes into two clusters, 79 rice accessions can be divided in to two broad groups: Indica and Japonica. Among the hill rices, only 18 accessions grouped with indica and aus control varieties, suggesting that the majority of hill rices are Japonica. Based on the STRUCTURE grouping at K = 4, 50 % of hill rice accessions were identified as japonica. The number of admixed individuals at K = 4 was considerable higher than that at K = 2, which suggests that the majority of hill rice cultivars are derived from admixture among the groups within either Indica or Japonica rather than between two subspecies. Similar findings were also reported in the Chinese rice collection . Genetic grouping of rice accessions in both NJ cluster and PCoA also supported the STRUCTURE results. The genetic grouping of rice accessions in this study did not support the farmer’s classification of hill rices. This discrepancy may be due to the limited number of SSR loci used in diversity analysis. These 35 SSR markers might reflect only a small part of the genotype-phenotype association of otherwise complex traits like grain and kernel length and shape, which are actually attributed to multiple loci with small effect [52, 53]. Similar observations were also reported in upland rice accessions of Nepal when analysed with 36 SSR markers , and in a set of 91 rice accessions from Eastern and NE India using 23 SSR markers .
In the present study, the genetic grouping of hill rices based on SSR data was also inconsistent with the geographical distribution (district-wise) and naming of rice cultivars. We cannot explain why geographical distribution did not match with genetic relationship among accessions. Seed exchange between the villages might be the reason. We noted that the farmers in Arunachal Pradesh most often gift seeds of good rice varieties during wedding ceremonies. It has been observed that the landraces with the same name collected from different villages/district are not always genetically identical. Similar observations were also reported in other studies [18, 48]. Therefore, during the germplasm collection, cultivars with the same names need to be collected as each accession is inherently valuable.
The current hill rice accessions revealed overall genetic diversity (H e = 0.66; PIC = 0.62) similar to the 107 aromatic rice cultivars of NE India (H e = 0.67; PIC = 0.62) . Slightly lower genetic diversity was reported in a subset of 26 rice accessions collected from NE India (PIC = 0.57) . The genetic diversity in the current hill rice accessions was considerably higher than that found in high-altitude rice cultivars of Nepal (PIC = 0.17) . Although allelic diversity indices could be used as indicators of genetic variation in germplasm collection, such estimates are relative and largely depend on the number of polymorphic loci and relatedness of genotypes included . The diverse nature of Arunachal Pradesh hill rices may be a reflection of the prevalent diverse agro-climatic, ethno-cultural and eco-geographical features of the state. Interestingly, the genetic diversity between the umte and tening groups of hill rice was at par, and pair-wise allelic difference was non-significant. This is perhaps due to similar genetic makeup of the accessions of two morphological groups. The genetic diversity estimates among four STRUCTURE subpopulations were higher for group 3 and 4 (H e = 0.69 and 0.59, respectively), which can be categorized as japonica groups. However, the gene diversity in nine accessions of groups 2 was slightly lower (H e = 0.53). This group was characterized by indica accessions, and higher genetic diversity of these accessions supports the fact that indica accessions are more diverse than japonica accessions [15, 56].
The present research, although with a limited number of SSR markers, indicated that hill rices were mostly japonica or admixture among the subpopulations of Indica or Japonica. Both FPTest and F st results indicated that the hill accessions are genetically diverse from temperate japonica and indica. The rice germplasm prevailing in NE India has great importance in studying the evolution and population genetics of cultivated rice. This region is considered to be the domestication origin of Indica rices . While, the adjoining Southern China is the proposed origin of Japonica rice [1, 57]. Therefore, the bordering region of China and India can be considered as the differentiating zone of the two subspecies. Many studies in the past have reported the occurrence of diverse rice ecotypes, belonging to various genetic groups of rice, in NE India and Bangladesh . Prevalence of Japonica rices has also been recorded in NE Indian rice germplasm [30, 59]. The hill rice germplasm reported in the present study may be a valuable resource for studying Indica-Japonica evolution.
The findings presented here also have some implications for the utilization and conservation of hill rice germplasm. First, considerable diversity, both in agro-morphological and genetic terms, has been maintained by the traditional farmers. This is most likely the result of maintaining different sets of landraces in different regions in response to local heterogeneities in climate and soil type, as well as indigenous culture, as demonstrated by Fu et al.  in Thai cassava farming. However, due to increased susceptibility to pests and diseases, some landraces have become obsolete. Second, our findings also are encouraging for on-farm conservation programme in the study location, along with the ex situ conservation efforts. The genetic variation among the cultivars from different districts was not significant, which indicates that on farm conservation programme can be initiated in any of the three districts of Arunachal Pradesh. Third, significant variations for many agro-morphological traits can be exploited in the rice breeding programmes. The information on genetic distance among the landraces could be used to guide the selection of valuable germplasm in rice breeding.
AMOVA, analysis of molecular variance; AN, average number of alleles per locus; ANOVA, analysis of variance; CDA, canonical discriminant analysis; CGIAR, Consultative Group on International Agricultural Research; CV, coefficient of variation; ICAR, Indian Council of Agricultural research; MAF, major allele frequency; NJ, neighbour joining; PCA, principal component analysis; PCoA, principal coordinate analysis; PCR, polymerase chain reaction; PIC, polymorphism information content; SD, standard deviation; SSR, Simple sequence repeats
The work was supported by Indian Council of Agricultural Research under the project IXX04376. The exploration trip for collecting hill rice landraces from Arunachal Pradesh, India was undertaken as per approved exploration programme of ICAR-National Bureau of Plant Genetic Resources, New Delhi, India. The authors wish to thank Dr SV Ngachan, Director, ICAR Research Complex for NEH Region for providing laboratory facilities to carry out the SSR genotyping. They also wish to thank Dr Trilochan Mohapatra, former Director, ICAR-National Rice Research Institute, Cuttack, Odisha for providing the control varieties of rice. The authors are thankful to Dr Mo-Hua Yang, Associate Professor, Central South University of Forestry and Technology, China for performing the FPTest analysis in R.
Availability of data and materials
The data sets supporting the results of this article are included within the article and its additional files.
SR, BCM conducted the exploration trip to Arunachal Pradesh for hill rice collection. SR conceived the experiment. SR, AB designed the experiment. BM, SR, AB, RMY conducted all the experiments pertaining to DNA extraction, PCR and SSR genotyping. SR, RMY conducted field experiments and phenotyping of rice accessions. SR conducted the data analysis and drafted the manuscript. KCB, AB, AKM revised the manuscript. All authors have read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Khush GS. Origin, dispersal, cultivation and variation of rice. Plant Mol Biol. 1997;35:25–34.View ArticlePubMedGoogle Scholar
- Wing RA, Ammiraju JSS, Luo MZX, Kim H, Yu YS, Kudrna D, et al. The Oryza Map Alignment Project: the golden path to unlocking the genetic potential of wild rice species. Plant Mol Biol. 2005;59:53–62.View ArticlePubMedGoogle Scholar
- Zhang QF, Li JY, Xue YB, Han B, Deng XW. Rice 2020: a call for an international coordinated effort in rice functional genomics. Mol Plant. 2008;1:715–9.View ArticlePubMedGoogle Scholar
- Oka HI. Origin of Cultivated Rice. Tokyo: Japan Scientific Societies Press; 1988.Google Scholar
- Sang T, Ge S. Genetics and phylogenetics of rice domestication. Curr Opin Genet Dev. 2007;17:533–8.View ArticlePubMedGoogle Scholar
- Sweeney M, McCouch S. The complex history of the domestication of rice. Ann Bot. 2007;100:951–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Glaszmann JC. Isozymes and classification of Asian rice varieties. Theor Appl Genet. 1987;74:21–30.View ArticlePubMedGoogle Scholar
- Ni J, Colowit PM, Mackill DJ. Evaluation of genetic diversity in rice subspecies using microsatellite markers. Crop Sci. 2002;42:601–7.View ArticleGoogle Scholar
- Garris AJ, Tai TH, Coburn J, Kresovich S, McCouch S. Genetic structure and diversity in Oryza sativa L. Genetics. 2005;169:1631–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Caicedo AL, Williamson SH, Hernandez RD, Boyko A, Fledel-Alon A, York TL, et al. Genome-wide patterns of nucleotide polymorphism in domesticated rice. PLoS Genet. 2007;3:1745–56.View ArticlePubMedGoogle Scholar
- Ebana K, Yonemaru J, Fukuoka S, Iwata H, Kanamori H, Namiki N, et al. Genetic structure revealed by a whole-genome single-nucleotide polymorphism survey of diverse accessions of cultivated Asian rice (Oryza sativa L.). Breeding Sci. 2010;60:390–7.View ArticleGoogle Scholar
- Zhao KY, Wright M, Kimball J, Eizenga G, McClung A, Kovach M, et al. Genomic diversity and introgression in O. sativa reveal the impact of domestication and breeding on the rice genome. PLoS ONE. 2010;5:e10780.View ArticlePubMedPubMed CentralGoogle Scholar
- Huang XH, Zhao Y, Wei XH, Li CY, Wang A, Zhao Q, et al. Genome-wide association study of flowering time and grain yield traits in a worldwide collection of rice germplasm. Nat Genet. 2012;44:32–53.View ArticleGoogle Scholar
- Wang C-H, Zheng X-M, Xu Q, Yuan X-P, Huang L, Zhou H-F, et al. Genetic diversity and classification of Oryza sativa with emphasis on Chinese rice germplasm. Heredity. 2014;112:489–96.View ArticlePubMedGoogle Scholar
- Gao LZ, ChiHong Z, LiPing C, JiZeng J, ZongEn Q, Dong YS. Microsatellite diversity within Oryza sativa with emphasis on indica-japonica divergence. Genet Res. 2005;85:1–14.View ArticlePubMedGoogle Scholar
- Prashanth SR, Parani M, Mohanty BP, Talame V, Tuberosa R, Parida A. Genetic diversity in cultivars and landraces of Oryza sativa subsp. indica as revealed by AFLP markers. Genome. 2002;45:451–9.View ArticlePubMedGoogle Scholar
- Barry MB, Pham JL, Noyer JL, Courtois B, Billot C, Ahmadi N. Implications for in situ genetic resource conservation from the ecogeographical distribution of rice genetic diversity in Maritime Guinea. Plant Genet Resour. 2007;5:45–54.View ArticleGoogle Scholar
- Thomson MJ, Polato NR, Prasetiyono J, Trijatmiko KR, Silitonga TS, McCouch SR. Genetic diversity of isolated populations of Indonesian landraces of rice (Oryza sativa L.) collected in East Kalimantan on the island of Borneo. Rice. 2009;2:80–92.View ArticleGoogle Scholar
- de Oliveira Borba TC, Brondani RP, Rangel PH, Brondani C. Microsatellite marker mediated analysis of the EMBRAPA Rice Core Collection genetic diversity. Genetica. 2009;137:293–304.View ArticlePubMedGoogle Scholar
- Zhang H, Sun J, Wang M, Liao D, Zeng Y, Shen S, et al. Genetic structure and phylogeography of rice landraces in Yunnan, China, revealed by SSR. Genome. 2007;50:72–83.View ArticlePubMedGoogle Scholar
- Zhang DL, Zhang HL, Wang MX, Sun JL, Qi YW, Wang F, et al. Genetic structure and differentiation of Oryza sativa L. in China revealed by microsatellites. Theor Appl Genet. 2009;119:1105–17.View ArticlePubMedGoogle Scholar
- Li X, Yan W, Agrama H, Hu B, Jia L, Jia M, et al. Genotypic and phenotypic characterization of genetic differentiation and diversity in the USDA rice mini-core collection. Genetica. 2010;138:1221–30.View ArticlePubMedGoogle Scholar
- Vaughan DA, Lu BR, Tomooka N. The evolving story of rice evolution. Plant Sci. 2008;174:394–408.View ArticleGoogle Scholar
- Hore DK. Rice diversity collection, conservation and management in North-eastern India. Genet Resour Crop Evol. 2005;52:1129–40.View ArticleGoogle Scholar
- Government of India. Basic statistics of North Eastern region 2015. Government of India, North Eastern Council Secretariat, Nongrim Hills, Shillong, India; 2015.p.411Google Scholar
- Directorate of Rice Development (2015) State-wise rice productivity analysis, State: Arunachal Pradesh. Directorate of rice development, Department of Agriculture and Co-operation, Patna, Bihar; http://drdpat.bih.nic.in. Accessed 30 December, 2015.
- Sthapit BR. Chilling injury of rice crop in Nepal: A review. J Inst Agric Anim Sci. 1992;13:1–32.Google Scholar
- McCouch SR, Wright MH, Tung C-W, Maron LG, McNally KL, Fitzgerald M, et al. Open access resources for genome-wide association mapping in rice. Nat Commun. 2016;7:10532.View ArticlePubMedPubMed CentralGoogle Scholar
- IBPGR-IRRI. Descriptors for rice (Oryza sativa L.) IBPGR-IRRI Rice Advisory Committee. Manila: IRRI; 1980.Google Scholar
- Roy S, Banerjee A, Mawkhlieng B, Misra AK, Pattanayak A, Harish GD, et al. Genetic diversity and population structure in aromatic and quality rice (Oryza sativa L.) landraces from North-Eastern India. PLoS ONE. 2015;10:e0129607.View ArticlePubMedPubMed CentralGoogle Scholar
- IBM Corp. IBM SPSS Statistics for Windows, version 20.0. Armonk: IBM Corp; 2011.Google Scholar
- Pritchard JK, Stephens M, Donnelly P. Inference of population structure using multilocus genotype data. Genetics. 2000;155:945–59.PubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Liu KJ, Muse SV. PowerMarker: an integrated analysis environment for genetic marker analysis. Bioinformatics. 2005;21:2128–9.View ArticlePubMedGoogle Scholar
- Goudet J. FSTAT version 1.2: a computer program to calculate F statistics. J Hered. 1995;86:485–6.Google Scholar
- Peakall R, Smouse PE. GenAlEx 6.5: genetic analysis in Excel. Population genetic software for teaching and research-an update. Bioinformatics. 2012;28:2537–9.View ArticlePubMedPubMed CentralGoogle Scholar
- R Core Team. R: A language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2016.Google Scholar
- Fu YB, Peterson GW, Scoles G, Rossnagel B, Schoen DJ, Richards KW. Allelic diversity changes in 96 Canadian oat cultivars released from 1886 to 2001. Crop Sci. 2003;43:1989–95.View ArticleGoogle Scholar
- Rohlf F. NTSYS-pc: Numerical Taxonomy and Multivariate Analysis System, Version 2.00. Setauket: Exeter Software; 1997.Google Scholar
- Cavalli-Sforza LL, Edwards AWF. Phylogentic analysis. Models and estimation procedures. Am J Hum Genet. 1967;19:233–57.PubMedPubMed CentralGoogle Scholar
- Huson DH, Scornavacca C. Dendroscope 3 - An interactive viewer for rooted phylogenetic trees and networks. Syst Biol. 2012;61(6):1061–7.View ArticlePubMedGoogle Scholar
- Fukuoka S, Suu T, Ebana K, Trinh L, Nagamine T, Okuno K. Diversity in phenotypic profiles in landrace populations of Vietnamese rice: a case study of agronomic characters for conserving crop genetic diversity on farm. Genet Resour Crop Evol. 2006;53:753–61.View ArticleGoogle Scholar
- Rana RB, Garforth C, Sthapit B, Jarvis D. Influence of socio-economic and cultural factors in rice varietal diversity management on-farm in Nepal. Agric Human Values. 2007;24:461–72.View ArticleGoogle Scholar
- Sharma RC, Chaudhary NK, Ojha B, Yadav L, Pandey MP, Shrestha SM. Variation in rice landraces adapted to the lowlands and hills in Nepal. Plant Genet Resour: Charact Util. 2007;5:120–7.View ArticleGoogle Scholar
- Bajracharya J, Rana RB, Gauchan D, Sthapit BR, Jarvis DI, Witcombe JR. Rice landrace diversity in Nepal. Socio-economic and ecological factors determining rice landrace diversity in three agro-ecozones of Nepal based on farm surveys. Genet Resour Crop Evol. 2010;57:1013–22.View ArticleGoogle Scholar
- Sthapit BR, Joshi KD, Rana RB, Upadhayay MP, Eyzaguirre P, Jarvis D. Enhancing biodiversity and production through participatory plant breeding: setting breeding goals. In: An exchange of experiences from South and South East Asia. Proceedings of the International Symposium on Participatory Plant Genetic Resource Enhancement. 2000.Google Scholar
- Roy S, Rathi RS, Misra AK, Bhatt BP, Bhandari DC. Phenotypic characterization of indigenous rice (Oryza sativa L.) germplasm collected from the state of Nagaland. Plant Genet Resour: Charact Util. 2014;12:58–66.View ArticleGoogle Scholar
- 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. Field Crop Res. 2006;95(2–3):327–35.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- Pessoa-Filho M, Beló A, Alcochete A, Rangel P, Ferreira M. A set of multiplex panels of microsatellite markers for rapid molecular characterization of rice accessions. BMC Plant Biol. 2007;7:23.View ArticlePubMedPubMed CentralGoogle Scholar
- Ueno K, Sato T, Takahashi N. The indica-japonica classification of Asian rice ecotypes and Japanese lowland and upland rice (Oryza sativa L.). Euphytica. 1990;46:161–4.View ArticleGoogle Scholar
- Ge XJ, Xing YZ, Xu CG, He YQ. QTL analysis of cooked rice grain elongation, volume expansion, and water absorption using a recombinant inbred population. Plant Breed. 2005;124:121–6.View ArticleGoogle Scholar
- Bai X, Luo L, Yan W, Rao MK, Zhan W, Xing Y. Genetic dissection of rice grain shape using a recombinant inbred line population derived from two contrasting parents and fine mapping a pleiotropic quantitative trait locus qGL7. BMC Genet. 2010;11:16.View ArticlePubMedPubMed CentralGoogle Scholar
- Das B, Sengupta S, Parida SK, Roy B, Ghosh M, Prasad M, et al. Genetic diversity and population structure of rice landraces from Eastern and North Eastern states of India. BMC Genet. 2013;14:71.View ArticlePubMedPubMed CentralGoogle Scholar
- Ren X, Jiang H, Yan Z, Chen Y, Zhou X, Huang L, et al. Genetic diversity and population structure of the major peanut (Arachis hypogaea L.) cultivars grown in China by SSR markers. PLoS ONE. 2014;9:e88091.View ArticlePubMedPubMed CentralGoogle Scholar
- Ali ML, McClung AM, Jia MH, Kimball JA, McCouch SR, Eizenga GC. A rice diversity panel evaluated for genetic and agro-morphological diversity between subpopulations and its geographic distribution. Crop Sci. 2011;51:2021–35.View ArticleGoogle Scholar
- Huang X, Kurata N, Wei X, Wang Z-X, Wang A, Zhao Q, et al. A map of rice genome variation reveals the origin of cultivated rice. Nature. 2012;490:497–501.View ArticlePubMedGoogle Scholar
- Travis AJ, Norton GJ, Datta S, Sarma R, Dasgupta T, Savio FL, et al. Assessing the genetic diversity of rice originating from Bangladesh, Assam and West Bengal. Rice. 2015;8:35.View ArticlePubMedPubMed CentralGoogle Scholar
- Choudhury B, Khan ML, Dayanandan S. Genetic structure and diversity of indigenous rice varieties (Oryza sativa) in Eastern Himalayan region of Northeast India. Springer Plus. 2013;2:228–37.View ArticlePubMedPubMed CentralGoogle Scholar
- Fu YB, Wangsomnuk PP, Ruttawa B. Thai elite cassava genetic diversity was fortuitously conserved through farming with different sets of varieties. Conserv Genet. 2014;15:1463–78.View ArticleGoogle Scholar