Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Genetic Diversity and Population Structure in Aromatic and Quality Rice (Oryza sativa L.) Landraces from North-Eastern India

  • Somnath Roy ,

    Contributed equally to this work with: Somnath Roy, Amrita Banerjee

    sroypbr@gmail.com

    Affiliation ICAR-National Bureau of Plant Genetic Resources, Regional Station, Umiam, Meghalaya, India

  • Amrita Banerjee ,

    Contributed equally to this work with: Somnath Roy, Amrita Banerjee

    Affiliation ICAR Research Complex for North Eastern Hill Region, Umiam, Meghalaya, India

  • Bandapkuper Mawkhlieng,

    Affiliation ICAR-National Bureau of Plant Genetic Resources, Regional Station, Umiam, Meghalaya, India

  • A. K. Misra,

    Affiliation ICAR-National Bureau of Plant Genetic Resources, Regional Station, Umiam, Meghalaya, India

  • A. Pattanayak,

    Current address: ICAR-Vivekananda Parvatiya Krishi Anusandhan Sansthan, Almora, Uttarakhand, India

    Affiliation ICAR Research Complex for North Eastern Hill Region, Umiam, Meghalaya, India

  • G. D. Harish,

    Affiliation ICAR-National Bureau of Plant Genetic Resources, Regional Station, Umiam, Meghalaya, India

  • S. K. Singh,

    Affiliation ICAR-National Bureau of Plant Genetic Resources, Regional Station, Umiam, Meghalaya, India

  • S. V. Ngachan,

    Affiliation ICAR Research Complex for North Eastern Hill Region, Umiam, Meghalaya, India

  • K. C. Bansal

    Affiliation ICAR-National Bureau of Plant Genetic Resources, Pusa Campus, New Delhi, India

Correction

20 Oct 2015: Roy S, Banerjee A, Mawkhlieng B, Misra AK, Pattanayak A, et al. (2015) Correction: Genetic Diversity and Population Structure in Aromatic and Quality Rice (Oryza sativa L.) Landraces from North-Eastern India. PLOS ONE 10(10): e0141405. https://doi.org/10.1371/journal.pone.0141405 View correction

Abstract

The North-eastern (NE) India, comprising of Arunachal Pradesh, Assam, Manipur, Meghalaya, Mizoram, Nagaland, Sikkim and Tripura, possess diverse array of locally adapted non-Basmati aromatic germplasm. The germplasm collections from this region could serve as valuable resources in breeding for abiotic stress tolerance, grain yield and cooking/eating quality. To utilize such collections, however, breeders need information about the extent and distribution of genetic diversity present within collections. In this study, we report the result of population genetic analysis of 107 aromatic and quality rice accessions collected from different parts of NE India, as well as classified these accessions in the context of a set of structured global rice cultivars. A total of 322 alleles were amplified by 40 simple sequence repeat (SSR) markers with an average of 8.03 alleles per locus. Average gene diversity was 0.67. Population structure analysis revealed that NE Indian aromatic rice can be subdivided into three genetically distinct population clusters: P1, joha rice accessions from Assam, tai rices from Mizoram and those from Sikkim; P2, chakhao rice germplasm from Manipur; and P3, aromatic rice accessions from Nagaland. Pair-wise FST between three groups varied from 0.223 (P1 vs P2) to 0.453 (P2 vs P3). With reference to the global classification of rice cultivars, two major groups (Indica and Japonica) were identified in NE Indian germplasm. The aromatic accessions from Assam, Manipur and Sikkim were assigned to the Indica group, while the accessions from Nagaland exhibited close association with Japonica. The tai accessions of Mizoram along with few chakhao accessions collected from the hill districts of Manipur were identified as admixed. The results highlight the importance of regional genetic studies for understanding diversification of aromatic rice in India. The data also suggest that there is scope for exploiting the genetic diversity of aromatic and quality rice germplasm of NE India for rice improvement.

Introduction

Rice (Oryza sativa L.) is rich in genetic diversity at both inter- and intra-specific levels. The genetic structure of rice is well characterized (see reviews in [1,2]). In addition to the two major subspecies, Indica and Japonica, five genetically defined groups, indica, aus, aromatic, temperate japonica and tropical japonica have been identified with genetic markers [3,4,5,6,7]. Aromatic/ quality rice is a special class of rice with high market value due to its superior grain qualities and pleasant aroma. Aromatic accessions have been identified within Group V (i.e., ‘Basmati’ and ‘Sadri’ varieties), indica (i.e., ‘Jasmine’ varieties) and tropical japonica [8,9]. The centre of diversity of aromatic rice in India is the foothills of Himalayas, covering Uttar Pradesh, Bihar and the Tarai region of Nepal [10]. From this centre of diversity, aromatic rice germplasm dispersed to the other states of India and neighbouring countries and adapted to the local environments. The genetic structure of Indian aromatic and quality rice has been studied in the past [11,12,13]. All these studies included very few aromatic rice accessions from the North-eastern (NE) states of India and thus did not provide the actual picture of genetic diversity in aromatic and quality rice germplasm prevailing in this region.

The NE India is a large geographical area comprising of eight states namely Arunachal Pradesh, Assam, Manipur, Meghalaya, Mizoram, Nagaland, Sikkim and Tripura. This region is surrounded by Bhutan and China in the North, Myanmar in the East and Bangladesh in the Southern side. The NE India is characterized by high rainfall, humidity, varied topography, heavy natural selection pressures and environmental stresses. Being a biodiversity hotspot [14], the region is rich both in floristic and crop diversities. Rice is the principal crop of this region possessing around 72% of the total cultivated area. Rice is cultivated in upland, lowland and deep water conditions. It has been estimated that at least 10,000 indigenous cultivars are prevailing in this region [15]. The farmers of this region still grow their heirloom cultivars which not only suit to their taste but also provide crop security. The NE India is also the home to many locally adapted aromatic and quality rice land races. Despite their low-yield potential, these cultivars are grown for their high market and social values. There are prominent cultivar groups within the aromatic rice gene pool of NE India. For example, joha, chakhao and tai cultivars grown in the states of Assam, Manipur and Mizoram, respectively. The traditional aromatic rice cultivars of Assam are locally known as joha. It is a special class under sali or winter rice. The joha rice cultivars are popular for its unique aroma, good cooking qualities and excellent palatability [16]. In the state of Manipur, the aromatic rice cultivars are locally known as chakhao (chak = rice and ahoba or hao = delicious). These landraces are grown both under upland and lowland ecologies. The population genetic analysis of chakhao rice accessions of Manipur have been conducted recently [17]. The tai rice cultivars of Mizoram are aromatic in nature and are mostly cultivated under wetland conditions. The aromatic rice germplasm collected from the state of Nagaland has been phenotypically characterized. This rice germplasm has considerable variation in grain size, shape, awn and glume characteristics [18]. A considerable number of aromatic rice cultivars, namely Krishna Bhog, Brimphul, Kalanunia and local Basmati are grown in the states of Sikkim and Tripura. The aromatic and quality rice cultivars are mainly used for preparing special dishes in festive occasions and marriage ceremonies.

Knowledge regarding the extent of genetic variation and genetic relationships between genotypes are vital for designing effecting breeding and conservation strategies. Population genetic analysis of rice landraces collected from a region also helps in understanding the complex interaction between rice diversity and human cultivation practices and culture, as the cultivar structure is shaped by the interplay between adaptation to the local environment and artificial selection imposed by the farmers. In contrast, the global studies though provided an excellent overview of the genetic structure of cultivated rice, very often include only a few cultivars from a region and thus unable to explain the rice diversity at a local scale. A large collection of rice germplasm including the aromatic and quality rice landraces from NE India is conserved at National Genebank, National Bureau of Plant Genetic Resources (NBPGR), New Delhi. These valuable aromatic rice germplasm has not been characterized at molecular level.

In the present study, a diverse set of 107 aromatic and quality rice accessions collected from various parts of the NE India were analysed using 40 genome-wide SSR markers. The objectives were (i) to estimate the level of genetic diversity and population structure and (ii) to address the genetic classification of the NE Indian rice accession in the context of structured global rice cultivars.

Materials and Methods

Plant materials

A total of 107 aromatic/ quality rice accessions (Assam: 54, Nagaland: 21, Manipur 16, Mizoram: 9, Sikkim: 5 and Arunachal Pradesh: 2) were sampled from the National Genebank, NBPGR based on characterization data and passport information. The information on the names, accession number, collection sites, kernel length, kernel shape, aroma level and inferred group as determined with STRUCTURE [19] of each accession is presented in S1 Table. The collection sites of these rice accessions are shown in Fig 1. In addition, a set of 67 structured global rice accessions (representing indica, aus, aromatic, temperate japonica, tropical japonica and admixed cultivars originating from India and its neighbouring countries) were selected from the study conducted by Garris et al. [3] (S2 Table). The genotyping data for the 30 SSR markers (common with the present study) were taken from Garris et al. [3] and used for a combined analysis.

thumbnail
Fig 1. Geographical location of North-eastern region of India.

The collection sites of rice landraces are indicated by circles. Smaller and bigger dots denote the collection sites with and without actual geographical coordinates, respectively.

https://doi.org/10.1371/journal.pone.0129607.g001

DNA extraction and SSR genotyping

Total genomic DNA was extracted from five day-old seedlings germinated from five well developed seeds of each accession using Qiagen DNeasy plant mini kit (Qiagen, CA, USA) following the protocol provided by manufacturer. Forty SSR markers distributed across the rice genome were designed from the ‘Gramene’ marker database (http://www.gramene.org/markers/microsat/). Detailed information about these markers is given in Table 1.

thumbnail
Table 1. Summary statistics of the 40 SSR markers used in this study.

https://doi.org/10.1371/journal.pone.0129607.t001

Polymerase chain reaction (PCR) was performed in a total volume of 25 μl containing 50 ng of template DNA, 0.2 mM of dNTPs, 2.5 mM MgCl2, 2.5 μl 10X PCR buffer, 2.5 pmol of each primer, 1.0 unit of Taq DNA polymerase and double-distilled water to make the final volume up to 25 μl. Amplification was carried out using a thermocycler (Mastercycler, Eppendorf, Hamburg, Germany) and PCR conditions were: 5 min at 94°C; 35 cycles of 30 s at 94°C 30 s at annealing temperature and 30 s at 72°C; followed by 7 min at 72°C. The amplified products were visualized by ethydium bromide stained 2.5% agarose gels in BioRad GelDoc XR system.

The molecular sizes (in nucleotide) of the amplified alleles were determined based on migration relative to the DNA ladder using the PyElph-1.4 gel image analysis software [20]. The band sizes of each marker were checked from the Gramene database (http://www.gramene.org/marker/microsat/). In case of non-amplification, we repeated the PCR to exclude technical failure and in case of failure in both PCRs a null allele was recorded. Genotype data of all accessions are given in S3 Table.

Data analysis

The average number of alleles per locus, major allele frequency, gene diversity, heterozygosity and polymorphism information content (PIC) were calculated using PowerMarker V3.25 [21].

Population structure was examined using the Bayesian model-based approach implemented in STRUCTURE V2.3.4 [19,22]. The number of clusters (K) evaluated here ranged from 1 to 10. The analysis was performed using five replicate runs per K value, a burn-in period length of 5000, a run length of 50000, 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 (K = 3 was optimum for this analysis) based on both the LnP(D) and Evanno’s ΔK [23]. Subsequently, 10 simulations at K = 3 were then performed with a burn-in period of 10000 and a run length of 100000. At K = 3, the membership coefficient from the run with the lowest likelihood value (-9442.1) was used to assign each accession to the K = 1 to 3 populations based on the highest membership coefficient. The STRUCTURE analysis of the combined dataset (NE Indian and global accessions) was performed following same procedures and run conditions. The graphical display of the structure results was generated using DISTRUCT programme [24].

Pair-wise values of F-statistics (FST) were estimated using Arlequin V3.5 [25]. Genetic distance was calculated using the C. S. Chord distances [26], and the phylogenetic tree was constructed based on neighbour-joining (NJ) method implemented in PowerMarker. The tree was visualized using Dendroscope V3 [27]. To summarize the major patterns of variation in multi-locus data set, principal coordinate analysis (PCoA) was performed using GenAlEx V6.5 [28]. The molecular variance of sub-populations and accessions within the sub-populations were calculated using an Analysis of Molecular Variance (AMOVA) approach in Arlequin.

Results

Microsatellite diversity

In this study, 40 SSR markers distributed over all 12 chromosomes of rice were genotyped in 107 accessions of aromatic and quality rice landraces of NE India. The summary statistics of the 40 SSR markers are given in Table 1. All the markers showed polymorphism and a total of 322 alleles were identified. The average number of alleles per locus (AN) was 8.03, ranging from two (RM338 and RM507) to 20 (RM207). The gene diversity or expected heterozygosity (He) ranged from 0.27 to 0.88, with an average of 0.67. The observed heterozygosity (Ho) was 0.01. The polymorphism information content (PIC) for the SSR loci ranged from 0.26 (RM178) to 0.87 (RM259).

Population structure and genetic relationships

Bayesian analysis of population structure using the model-based approach in STRUCTURE provided support for the existence of population structure in aromatic and quality rice accessions of NE India. The LnP(D) and Evanno’s ΔK method supported the presence of three genetically distinct clusters (i.e., K = 3; S1 Fig), here denoted as P1, P2 and P3, respectively. The inferred population structure is given in Fig 2a. P1 contained 69 accessions, among which 54 were from Assam (joha accessions), nine were from Mizoram (tai accessions), five were from Sikkim, and one (cv. Amo) from Arunachal Pradesh. This group was dominated by accessions with short to medium kernel length, bold to medium kernel shape and high aroma. P2 had 22 accessions, most from Nagaland except one (cv. Khaw) from Arunachal Pradesh. They have long to extra-long kernel length, predominantly bold kernel shape and medium aroma. P3 contained 16 accessions, all from Manipur (chakhao accessions). This group was dominated by medium to long kernel length, medium kernel shape and medium to high aroma. The three sub-populations (P1–P3) had an Fst (population-specific Fst values calculated in STRUCTURE) of 0.208, 0.67 and 0.183, respectively with an average value of 0.354, indicating moderate population structure. Overall, 12 (11.2%) admixed accessions were found at K = 3.

thumbnail
Fig 2. Population structure of 107 aromatic and quality rice accessions of North-eastern India.

(a) Model-based clustering using STRUCTURE analysis. (b) NJ tree based on C. S. Chord genetic distance. (c) Principal coordinate analysis.

https://doi.org/10.1371/journal.pone.0129607.g002

The neighbour-joining (NJ) tree and principal coordinate analysis (PCoA) further confirmed the STRUCTURE results. The tree model-based groups (P1–P3) clearly separated in NJ tree (Fig 2b). Note that the accessions of P1 further divided into four sub-clusters. The largest genetic distance (0.66) was noted between P2 and P3, while, the genetic distances between P1 vs P2 and P1 vs P3 were similar (0.53 and 0.55, respectively). In the PCoA (Fig 2c), the first two eigenvectors clearly separated the accessions into three groups, similar to those observed in STRUCTURE analysis and NJ tree.

The analysis of molecular variance (AMOVA) among populations indicated that 30.5% of the variation was due to differences among groups with the remaining 68.2% due to difference within groups (Table 2). Around 1% of the total genetic variance was explained by differences at the individual level. Pair-wise estimates of FST (S4 Table) indicated that the highest degree of differentiation was between P2 and P3 (0.453), and the lowest was between P1 and P2 (0.223).

thumbnail
Table 2. Analysis of molecular variance of aromatic and quality rice landraces of North-eastern India.

https://doi.org/10.1371/journal.pone.0129607.t002

Genetic diversity

The amount and organization of genetic diversity differed among the model-based populations (Table 3). The highest mean major allele frequency was present in P2 (0.819) and the lowest was found in P3 (0.507). The AN was the highest in accessions of P1 (6.13) and the least in accessions of P2 (2.45). Genetic diversity varied slightly among P1 and P3 (He = 0.601 and 0.61, respectively; PIC = 0.56 in each population). The genetic diversity in P2 was comparatively low (He = 0.253; PIC = 0.227).

thumbnail
Table 3. Summary statistics of microsatellite diversity in 107 aromatic and quality rice landraces of North-eastern India.

https://doi.org/10.1371/journal.pone.0129607.t003

Combined analysis of NE Indian and global rice accessions

To establish genetic grouping of the aromatic and quality rice landraces from NE India, a combined dataset including genotyping data from Garris et al. [3] was subjected to STRUCTURE analysis (174 genotypes, 30 SSR markers). Since a total of 169 SSR markers were originally used by Garris et al. [3] to infer the population structure, and only 30 of them were included in the combined dataset, the 67 global rice accessions were independently reanalysed to assess the impact of the reduction in marker number. No significant differences were observed in the structure outputs produced by the 30 SSR data matrices (K = 1–6) and the accessions were divided into same groupings described by Garris et al. [3] (data not shown). For all subsequent analyses, the global accessions were labelled following the grouping proposed by Garris et al. [3].

It was observed that LnP(D) value increased with K from 1 to 10 (S2a Fig), with highest log likelihood score at K = 2. The procedure of Evanno et al. [23] also detected the maximal ΔK at K = 2 followed by K = 3. The ΔK value decreased with increased K, and no peak of ΔK was observed at K > 4 (S2b Fig). Since there was no clear indication of which K value (except K = 2) providing the best fit for the data, the population sub-groups were examined for biological relevance from K = 2–5 (Fig 3a). At K = 2, two major groups corresponding to two rice subspecies (Indica and Japonica) were apparent within the population. Out of 174 accessions, 100 (57%) were Indica, 42 (24%) were Japonica and 32 (18%) were admixed. Whereas, within NE Indian accessions 68 (63%) were Indica, 22 (21%) were Japonica and 17 (15%) were admixed (S1 Table). All the aromatic rice accessions from Nagaland were identified as japonica. The aromatic rice accessions from Mizoram were found to be admixed. Some of chakhao accessions (Kathaibuw, Napnag hangmei, Chakhao phou, Makrei and Manui maa) also identified as admixed. Note that these chakhao accessions were collected from the hill districts of Manipur. The clustering pattern obtained for K = 3 revealed an association between: (i) group P3 (chakhao) and indica accessions; and (ii) group P2 and japonica accessions. Admixed individuals seem to be prevalent in group P1 (Fig 3a). The grouping of indica,aus, aromatic and japonica accessions becomes apparent at K = 5. At this level, group P1 showed a closer affiliation to the aus and aromatic accessions. Most of the chakhao accessions of P3 showed clear association with indica accessions, while few identified as indica/japonica admixture. The aromatic accessions of Nagaland differentiated from the japonica cultivars at K = 4 and 5 (Fig 3a).

thumbnail
Fig 3. Population structure of the combined samples of North-eastern India and global rice cultivars.

(a) Model-based population assignment using STRUCTURE analysis. (b) NJ tree based on C. S. Chord genetic distance. (c) Principal coordinate analysis.

https://doi.org/10.1371/journal.pone.0129607.g003

The grouping of the NE Indian aromatic rices in relation to global rice cultivars was also analysed in NJ tree and PCoA. In the NJ Tree the rice cultivars of model-based groups P1 showed a closer affiliation to the aus and indica group (Fig 3b). The accessions of P3 were clustered close to aromatic group. The accessions of P2 showed closer affiliation to the temperate japonica group and admixed accessions (Fig 3b). The results of PCoA supported the STRUCTURE results obtained at K = 2 (Fig 3c). Moreover, the cluster analysis of 107 NE Indian aromatic rice accessions using 30 SSR data matrices showed similar grouping to that of previous analysis using 40 SSRs (data not shown).

Discussion

Genetic diversity

Genetic diversity and population structure of NE Indian rice germplasm has been analyzed in past using molecular markers [29,30,31,32]. In a landmark study, Choudhury et al. [33] characterized 6,984 NE Indian rice accessions from National Genebank using 36 SNP markers. Despite the availability of tremendous diversity, very little is known about the amount and distribution of genetic diversity and population structure within aromatic and quality rice germplasm from NE India. Knowledge of the genetic structure within this rice germplasm will be of great value for utilization and development of improved cultivars. The past studies on the genetic diversity of Indian aromatic and quality rice cultivars included only a few (3–5) genotypes from NE India [11,13]. Thus these studies were unable to explain the actual diversity prevailing in this region. Recently, we have reported the population structure of aromatic (chakhao) rice landraces collected from Manipur [17]. Therefore, the present study is the first effort to characterize a major portion of aromatic and quality rice diversity available in NE India.

Analysis of genetic diversity parameters indicated a high level of variability in the NE Indian aromatic and quality rice collection. The mean number of alleles per locus (8.03) was similar to the estimates reported for 52 Indian aromatic rice cultivars (7.8) [13] and 83 accessions of rice cultivars from Eastern and NE Indian states (7.9) [32]. However, the average number of alleles per locus for this study was lower than the estimated recorded in the diverse sets of global rice accessions [34,35,36]. In this study the average PIC value (0.62) was similar to the values obtained for Indian aromatic and quality rice cultivars (0.6) [13] and for a subset of 26 aromatic rice cultivars from West Bengal, India (0.65) [32]. The gene diversity index in the NE Indian rice landraces (0.67) was at par with the values reported for indigenous rice varieties of NE India (0.78) [31] and 234 global rice accessions from 34 countries (0.70) [3]. Although allelic diversity and gene diversity index could be used as indicators of genetic variation, such estimates are relative and largely depend on the number of polymorphic loci and relatedness of genotypes included [37]. Despite their marginal location, aromatic and quality rice landraces from NE India exhibited high levels of genetic diversity. The diverse nature of these landraces may be a reflection of the prevalent diverse agro-climatic, ethnic and eco-geographical features of the region.

Genetic structure

The Indian aromatic and quality rice germplasm is comprised of small-, medium-, and long-grain types with mild to strong aroma [38]. Based on conventional taxonomy, Indian aromatic rices have been classified as indicas [10]. Subsequent studies based on isozyme and SSR markers have shown that most of the aromatic rice cultivars of the Indian sub-continent, including Basmati types, are identifiable as a genetically distinct cluster [3,8,11,12]. Jain et al. [13] reported two major groups among a set of 52 Indian aromatic and quality rice accessions based on SSR marker analysis. They classified short-grained joha rice cultivars of Assam as aromatic indicas and proposed that these cultivars might have evolved via hybridization between indica and Basmati rice varieties. In the present study, population structure analysis revealed three subpopulations (P1–P3) out of which the majority of accession included in P1. This grouping is well in agreement with the genetic distance based clustering and PCoA. Analysis of rice accessions, both aromatic and non-aromatic, from NE Indian states revealed 2–4 clusters. Based on morphological analysis, Vairavan et al. [39] reported three clusters within a set of 400 accessions from NE India. Similar grouping was obtained by analysing 6,984 NE Indian rice accessions using SNP markers [33]. Choudhury et al. [31] reported two clusters within 24 indigenous and improved rice varieties of NE India, while Das et al. [32] found four groups among a set of 26 rice cultivars. The rice population dynamics is inextricably linked to human activities. In addition to the biological determinants, farmers’ management practices and preferences greatly influence the distribution of genetic diversity and the levels of gene flow. In the present study, it is apparent that aromatic landraces grown in Assam (joha), Mizoram (tai) and Sikkim are mostly short to medium in kernel length, bold in kernel shape and strongly aromatic. These landraces are both genetically and morphologically distinct from those of Manipur (chakhao) and Nagaland. Moreover, the aromatic landraces of Manipur are distinct from that of Nagaland. Combined analysis of NE Indian and global rice accessions using multiple approaches including model-based grouping, NJ tree and PCoA confirmed the assignment of the aromatic accessions from Assam, Manipur and Sikkim to the Indica group, while the accessions from Nagaland exhibited a close association with Japonica. The japonica cultivars have overall lower genetic diversity than indica [8,40,41]. In the present study, the genetic diversity within Nagaland accessions was also lower than the values recorded in other groups. The tai accessions of Mizoram along with few chakhao accessions collected from the hill districts of Manipur have been identified as admixed. Although STRUCTURE analysis clearly identified the uppermost level of structure for the populations studied here, the difficulties encountered in determining the optimal value of K may obscure the detection of further sub-structuring with real biological significance. From the overall results, it is proposed that the joha, tai and aromatic accessions of Sikkim are intermediate accessions between aus and aromatic. The chakhao accessions are aromatic indicas. The aromatic landraces from Nagaland can be classified as temperate japonicas based on their position in the NJ tree and PCoA and their lower pair-wise FST values with temperate japonicas in comparison to tropical japonicas.

In conclusion, SSR marker based molecular characterization of NE Indian aromatic and quality rice landraces revealed that large variation exists among the accessions. The population structure analysis identified three major groups within these landraces which largely correlate with the state-wise grouping as well as farmers’ classification of the landraces. Future population genetics-based studies including extensive collections of rice genetic resources from all the states of NE India will help in exploiting this rice gene pool more effectively in rice improvement programmes. In both domestic and international markets Basmati rices command premium prices, typically two to three times higher than traditional aromatic and quality rice cultivars due to consumer quality preference. Some of the traditional aromatic rice of NE India such as chakhao cultivars of the state of Manipur has export potentiality in South-east Asian countries where consumers prefer sticky rice. In contrast to Basmati rices, the traditional aromatic landraces can be used for a variety of cuisine starting from table rice to deserts, from offerings in various religious ceremonies and social gatherings to diet for the convalescent [32]. These landraces are also used for making popped rice. Greater appreciation of the genetic diversity preserved in the aromatic rice gene pool of NE India will facilitate proper maintenance, conservation and utilization of this valuable resource.

Supporting Information

S1 Fig. Evaluation of STRUCTURE outputs for aromatic and quality rice landraces from North-eastern India.

(a) Mean LnP(D) over five runs for each K value. (b) Rate of change in the log probability of data between successive K values (ΔK).

https://doi.org/10.1371/journal.pone.0129607.s001

(TIF)

S2 Fig. Evaluation of STRUCTURE outputs for the combined dataset including aromatic and quality rice landraces from North-eastern India and global rice cultivars.

(a) Mean LnP(D) over five runs for each K value. (b) Rate of change in the log probability of data between successive K values (ΔK).

https://doi.org/10.1371/journal.pone.0129607.s002

(TIF)

S1 Table. Information on the aromatic and quality rice accessions of North-eastern India used in this study.

https://doi.org/10.1371/journal.pone.0129607.s003

(DOC)

S2 Table. Details of global rice cultivars taken from Garris et al., 2005.

https://doi.org/10.1371/journal.pone.0129607.s004

(DOC)

S3 Table. Genotyping data of 107 aromatic and quality rice landraces of North-eastern India.

https://doi.org/10.1371/journal.pone.0129607.s005

(XLS)

Acknowledgments

We thank Dr DK Hore, Former Officer-In-Charge, ICAR-National Bureau of Plant Genetic Resources (ICAR-NBPGR), Regional Station, Umiam, Meghalaya for critical reading and constructive suggestions on the manuscript. The help rendered by Mr E Khraw, Senior Research Fellow of ICAR Research Complex for NEH Region, Umiam, Meghalaya during gel electrophoreses is also acknowledged.

Author Contributions

Conceived and designed the experiments: SR. Performed the experiments: AB BM SR SKS GDH. Analyzed the data: SR AB AP. Contributed reagents/materials/analysis tools: KCB SVN AP AKM. Wrote the paper: SR AB AP.

References

  1. 1. Sang T, Ge S (2007) Genetics and phylogenetics of rice domestication. Curr Opin Genet Dev 17: 533–538. pmid:17988855
  2. 2. Sweeney M, McCouch S (2007) The complex history of the domestication of rice. Ann Bot-London 100: 951–957. pmid:17617555
  3. 3. Garris AJ, Tai TH, Coburn J, Kresovich S, McCouch S (2005) Genetic structure and diversity in Oryza sativa L. Genetics 169: 1631–1638. pmid:15654106
  4. 4. Caicedo AL, Williamson SH, Hernandez RD, Boyko A, Fledel-Alon A, York TL, et al. (2007) Genome-wide patterns of nucleotide polymorphism in domesticated rice. PLoS Genet 3: 1745–1756. pmid:17907810
  5. 5. Zhao KY, Wright M, Kimball J, Eizenga G, McClung A, Kovach M, et al. (2010) Genomic diversity and introgression in O. sativa reveal the impact of domestication and breeding on the rice genome. PLoS One 5: e10780. pmid:20520727
  6. 6. Huang XH, Zhao Y, Wei XH, Li CY, Wang A, Zhao Q, et al. (2012) Genome-wide association study of flowering time and grain yield traits in a worldwide collection of rice germplasm. Nat Genet 44: 32–53. pmid:22623593
  7. 7. Wang C-H, Zheng X-M, Xu Q, Yuan X-P, Huang L, Zhou H-F, et al. (2014) Genetic diversity and classification of Oryza sativa with emphasis on Chinese rice germplasm. Heredity 112: 489–496. pmid:24326293
  8. 8. Glaszmann JC (1987) Isozymes and classification of Asian rice varieties. Theor Appl Genet 74: 21–30. pmid:24241451
  9. 9. Kovach MJ, Calingacion MN, Fitzgerald MA, McCouch SR (2009) The origin and evolution of fragrance in rice (Oryza sativa L.) Proc Natl Acad Sci USA 106: 14444–14449. pmid:19706531
  10. 10. Khush GS (2000) Taxonomy and origin of rice. In: Singh RK, Singh US, Khush GS (eds) Aromatic rices. Oxford & IBH, New Delhi, pp 5–13.
  11. 11. Aggarwal RK, Shenoy VV, Ramadevi J, Rajkumar R, Singh L (2002) Molecular characterization of some Indian Basmati and other elite rice genotypes using fluorescent-AFLP. Theor Appl Genet 105: 680–690. pmid:12582481
  12. 12. Nagaraju J, Kathirvel M, Kumar RR, Siddiq EA, Hasnain SE (2002) Genetic analysis of traditional and evolved Basmati and non-Basmati rice varieties by using fluorescence-based ISSR-PCR and SSR markers. Proc Natl Acad Sci USA 99: 5836–5841. pmid:11959900
  13. 13. Jain S, Jain KR, McCouch SR (2004) Genetic analysis of Indian aromatic and quality rice (Oryza sativa L.) germplasm using panels of fluorescently-labeled microsatellite markers. Theor Appl Genet 109: 965–977. pmid:15309297
  14. 14. Myers N, Mittermeier RA, Mittermeier CG, da Fonseca GAB, Kent J (2000) Biodiversity hotspots for conservation priorities. Nature 403: 853–858. pmid:10706275
  15. 15. Hore DK (2005) Rice diversity collection, conservation and management in North-eastern India. Genet Resour Crop Evol 52: 1129–1140.
  16. 16. Das A, Kesari TV, Rangan L (2010) Aromatic joha rice of Assam-A review. Agric Rev 31: 1–10.
  17. 17. Roy S, Banerjee A, Pattanayak A, Roy SS, Rathi RS, Misra AK, et al. (2014a). Chakhao (delicious) rice landraces (Oryza sativa L.) of North-east India: collection, conservation and characterization of genetic diversity. Plant Genet Resour: Charact Util 12: 264–272.
  18. 18. Roy S, Rathi RS, Misra AK, Bhatt BP Bhandari DC (2014b) Phenotypic characterization of indigenous rice (Oryza sativa L.) germplasm collected from the state of Nagaland. Plant Genet Resour: Charact Util 12: 58–66.
  19. 19. Pritchard JK, Stephens M, Donnelly P (2000) Inference of population structure using multilocus genotype data. Genetics 155: 945–959. pmid:10835412
  20. 20. Pavel AB, Vasile CI (2012) PyElph—a software tool for gel images analysis and phylogenetics. BMC Bioinformatics 13: 9. pmid:22244131
  21. 21. Liu KJ, Muse SV (2005) PowerMarker: an integrated analysis environment for genetic marker analysis. Bioinformatics 21: 2128–2129. pmid:15705655
  22. 22. Falush D, Stephens M, Pritchard JK (2003) Inference of population structure using multilocus genotype data: linked loci and correlated allele frequencies. Genetics 164: 1567–1587. pmid:12930761
  23. 23. Evanno G, Regnaut S, Goudet J (2005) Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol Ecol 14: 2611–2620. pmid:15969739
  24. 24. Rosenberg NA (2004) DISTRUCT: a program for the graphical display of population structure Molecular Ecology Notes 4: 137–138.
  25. 25. Excoffier L, Lischer HEL (2010) Arlequin suite ver 3.5: A new series of programs to perform population genetics analyses under Linux and Windows. Mol Ecol Resour 10: 564–567. pmid:21565059
  26. 26. Cavalli-Sforza LL, Edwards AWF (1967) Phylogentic analysis. Models and estimation procedures. Am J Hum Genet 19: 233–257. pmid:6026583
  27. 27. Huson DH, Scornavacca C (2012) Dendroscope 3: An interactive tool for rooted phylogenetic trees and networks. Syst Biol 61: 1061–1067. pmid:22780991
  28. 28. Peakall R, Smouse PE (2012) GenAlEx 6.5: genetic analysis in Excel. Population genetic software for teaching and research—an update. Bioinformatics 28: 2537–2539. pmid:22820204
  29. 29. Sarma RN, Bahar B (2005) Genetic variation of bora rice (glutinous rice) of Assam as revealed by RAPDs. Plant Genet Resour Newsl FAO, IPGRI 144: 34–38.
  30. 30. Bhuyan N, Borah BK, Sarma RN (2007) Genetic diversity analysis in traditional lowland rice (Oryza sativa L.) of Assam using RAPD and ISSR markers. Curr Sci 93: 967–972.
  31. 31. Choudhury B, Khan ML, Dayanandan S (2013) Genetic structure and diversity of Indigenous Rice varieties (Oryza sativa) in Eastern Himalayan region of Northeast India. Springer Plus 2: 228–237. pmid:23741655
  32. 32. Das B, Sengupta S, Parida SK, Roy B, Ghosh M, Prasad M, et al. (2013) Genetic diversity and population structure of rice landraces from Eastern and North Eastern states of India. BMC Genet 14: 71. pmid:23945062
  33. 33. Choudhury DR, Singh N, Singh AK, Kumar S, Srinivasan K, Tyagi RK, et al. (2014) Analysis of genetic diversity and population structure of rice germplasm from North-eastern region of India and development of a core germplasm set. PLOS One 9: e113094. pmid:25412256
  34. 34. Xu YB, Beachell H, McCouch SR (2004) A marker-based approach to broadening the genetic base of rice in the USA. Crop Sci 44: 1947–1959.
  35. 35. Thomson MJ, Septiningsih EM, Suwardjo F, Santoso TJ, Silitonga TS,McCouch SR (2007) Genetic diversity analysis of traditional and improved Indonesian rice (Oryza sativa L.) germplasm using microsatellite markers. Theor Appl Genet 14: 559–568.
  36. 36. Borba TCDO, Brondani RPV, Rangei PHN, Brondani C (2009) Microsatellite marker-mediated analysis of the EMBRAPA rice core collection genetic diversity. Genetics 137: 293–304.
  37. 37. Ren X, Jiang H, Yan Z, Chen Y, Zhou X, Huang L, et al. (2014) Genetic Diversity and Population Structure of the Major Peanut (Arachis hypogaea L.) Cultivars Grown in China by SSR Markers. PLOS One 9: e88091. pmid:24520347
  38. 38. Singh RK, Singh US, Khush GS, Rohilla R, Singh JP, Singh G, et al. (2000) Small and medium grained aromatic rices of India. In: Singh RK, Singh US, Khush GS (eds) Aromatic rices. Oxford & IBH, New Delhi, pp 155–177.
  39. 39. Vairavan S, Siddiq EA, Arunachalam V, Swaminathan MS (1973) A study on the nature of genetic divergence in rice from Assam and Northeast Himalayas. Theor Appl Genet 43: 213–221. pmid:24425072
  40. 40. Zhang QF, Maroof MAS, Lu TU, Shen BZ (1992) Genetic diversity and differentiation of Indica and Japonica rice detected by RFLP analysis. Theor Appl Genet 83: 495–499. pmid:24202597
  41. 41. Ni J, Colowit PM, Mackill DJ (2002) Evaluation of genetic diversity in rice subspecies using microsatellite markers. Crop Sci 42: 601–607.