Genetic Diversity Analysis of Colored and White Rice Genotypes Using Microsatellite (SSR) and Insertion-Deletion (INDEL) Markers
Genetic Diversity Analysis of Colored and White Rice Genotypes Using Microsatellite (SSR) and Insertion-Deletion (INDEL) Markers
Genetic Diversity Analysis of Colored and White Rice Genotypes Using Microsatellite (SSR) and Insertion-Deletion (INDEL) Markers
26 (6): 497-507
doi: 10.9755/ejfa.v26i6.15722
http://www.ejfa.info/
Regular Article
PLANT MOLECULAR BIOLOGY
Department of Agricultural Botany, Anand Agricultural University, Anand 388110, Gujarat, India
Abstract
Genetic diversity analysis of 19 coloured and white rice genotypes were conducted using 14 Simple Sequence
Repeat (SSR) and 21 Insertion Deletion (INDEL) markers. Among them, polymorphic results were obtained for
9 SSR and 12 INDEL markers. The PIC values ranged from 0.36(RM484) to 0.78 (RM 167) and 0.50 (R9M20)
and 0.81 (R9M10) for SSR and INDEL markers respectively. In the case of both SSR and INDEL markers
highest genetic diversity was observed between Krishna Kamod (White pericarp) and IRST 1 (Red pericarp)
whereas least genetic diversity was observed between Lal Kada (Red pericarp) and Krishna Kamod (White
pericarp). However, it was also found that brown, black and red pericarp share more similarity among
themselves. It was also observed that INDEL markers reveals greater diversity among the genotypes as
compared to SSR markers which was indent from the low average similarity index observed in the former. On
the whole INDEL markers were found to be more efficient than SSR markers for diversity analysis.
497
Siddhi Patel et al.
RAPD, AFLP are extensively been used to quantify Materials and Methods
the inter-specific, intra-specific and inter-generic The study was conducted at Plant Biotechnology
variability in different plant groups and crop Laboratory, Department of Agricultural Botany, B.A
varieties (Pawar et al., 2013). Kibria et al. (2009) College of Agriculture, Anand Agricultural
has screened several rice varieties for studying the University, Anand during 2011. The seeds of 19 rice
genetic diversity by using SSR and RAPD markers. genotypes comprising of coloured rice (red, black,
SSR markers can estimate genetic diversity brown and white) used in the present study was
between cultivars e.g. between parents of genepool obtained from the Main Rice Research Station
or between plants extracted from a population or (MRRS), Anand Agricultural University, Nawagam
between populations. Microsatellites are more (Table 1). Few healthy seeds were sown and allowed
powerful for the identification of within cultivar to grow for three weeks in the pots. For the proper
variation (Lapitan et al., 2007). Zhou et al. (2003) growth and emergence of healthy seedlings, these pots
investigated the genetic diversity and genetic were watered regularly and proper light and
structure of natural populations of O. rufipogon in temperature conditions were maintained. Total
China using SSR markers and information was genomic DNA extraction from leaves of three weeks
found to be significant. Steele et al. (2008) reported old seedlings were carried out by Cetyl trimethyl
that a genome wide rice polymorphism database ammonium bromide (CTAB) method (Ahmadikhah
developed by Shen et al. (2004) contains more than et al., 2006). 300mg tissue samples were
4,00,000 insertion/ deletion polymorphisms homogenized in liquid nitrogen prior mixing with 800
(InDels). They used 50 insertion/ deletion markers µl of extraction buffer and were incubated for 1 hour
(InDels) for distinguishing between indica and at 65º C in water bath. Later on Chloroform: isoamyl
japonica variety. These are co- dominant markers alcohol (24:1) was added in the tubes and centrifuged
that give two possible alleles according to presence at 4ºC for 15 minutes at 12,000rpm. The supernatant
or absence of insertion sequence, situated between was collected and washed 1-2 times with Chloroform:
the primers. isoamyl alcohol (24:1) and kept for precipitation with
The present research is envisaged with an absolute alcohol for overnight. The samples were
objective to access genetic diversity among 19 centrifuged at 5,000 rpm for 15 min to get DNA in
colored rice genotypes (selected according to the pellet form. The pellets were further washed with
availability of genotypes) linked SSR markers as 70% alcohol and kept for drying and the quality was
these markers are highly polymorphic and easy to confirmed through Nanodrop N.D. 1000 software (ver
detect. The INDELs selected were identified from 3.7.1). For SSR and INDEL analysis DNA
12 chromosomal locations. amplification was performed in 25µl reaction volume
containing 50ng genomic DNA, 10X polymerase
buffer, 25mM dNTPs, 0.5 µl of each primer
(10pmol), 1 unit of Taq DNA polymerase using
Table 1. List of genotypes used in the present study. Applied Biosystem Thermal Cyclers. The cycling
conditions were: 1 cycle of 94°C for 7min followed
Sr. No Genotype Pericarp color
by 40 cycles of 45 sec each of 94°C, 48°C, 72°C and
1 RRT IRST 1 RED
finally 1 cycle of 72°C for 5min. Total of 10 µl
2 RRT IRST 2 RED
3 RRT IRST4 RED aliquots of the amplification products loaded in 3%
4 RRT IRST7 RED (w/v) agarose gel for electrophoresis in 1X TBE
5 RRT IRST 10 RED buffer and stained with Ethidium bromide and
6 RRT IRST 13 RED documented using SYNGENE GENESNAP G-BOX
7 RRT IRST 14 RED gel documentation system. These photographs were
8 RRT IRST 16 RED used to score the DNA bands for analysis. The gels
9 RRT IRST 19 RED were scored for computer analysis on the basis of the
10 RRT IRST 41 RED presence and absence of the amplified products. If a
11 RRT IRST 44 RED product was present in a genotype, it was designated
12 RRT IRST 47 RED
as ‘1’ and if absent, it was designated as ‘0’. A total of
13 RRT IRST 50 RED
14 TSP RED 14 SSR and 21 INDEL markers used for
15 A004 BLACK characterizing nineteen rice genotypes. Both SSR and
16 077 BLACK INDEL data were analyzed using NTSYS-PC
17 BC1 BROWN (Numerical Taxonomical and Multivariant Analysis
18 LAL KADA RED System computer package). The data were used to
19 KRISHNA KAMOD WHITE generate Jaccard’s similarity coefficients for SSR and
498
Emir. J. Food Agric. 2014. 26 (6): 497-507
http://www.ejfa.info/
INDEL bands. The Jaccard’s coefficients between polymorphism among the varieties was evaluated by
each pair of accessions were used to construct a calculating the PIC values for each of the nine
dendrogram using the unweighted pair group method microsatellite markers. Among the nine
with arithmetic averages (UPGMA). The PIC, microsatellite markers used lowest PIC of 0.36 was
hetrozygosity, and allelic diversity measures marker recorded in RM 484 while highest PIC of 0.78 was
informativeness and allelic and genotypic frequency obtained with RM 167. The Cluster analysis divided
were calculated. These measures are calculated using the genotypes into two major groups A and B. In this
PROC ALLELE procedure of SAS/GENETICS 9.3 group TSP and IRST50 showed highest genetic
version. similarity. First cluster (A) was further sub-divided
into two minor clusters i.e., A1 and A2. Second
Results and Discussion
major cluster B included two minor clusters B1 and
The present investigation on assessment of
B2 (Figure 1). The second major cluster indicated
pigmented colored rice varieties was carried out in
that black and brown rice are genetically more
the Biotechnology Laboratory of Department of
closely related to red rice. Highest genetic diversity
Agricultural Botany, B. A. College of Agriculture,
was shown between Krishna Kamod and IRST1
Anand Agriculture University, Anand, Gujarat,
whereas least genetic diversity was observed
India. Total of 14 SSR and 21 INDEL markers
between Lal Kada and Krishna Kamod, IRST1 and
were used for the study out of which 9 SSR
IRST2. In the present study multiple loci were
markers and 12 INDEL markers gave polymorphic
detected by RM 166, RM484, RM154, RM174,
results (Table 2).
RM167 and RM102. Occurrence of rare alleles
Microsatellite (SSR) markers are the PCR based
might have resulted from unequal crossing over,
markers that have been developed in many plant
translocation, other types of mutation or residual
species; they have an advantage of being multi
heterozygosity (Lapitan et al., 2007). Similar work
allelic, highly polymorphic and codominant. In the
has been carried out by Rajguru et al. (2005) where
present study total of 14 primers were used to
genetic diversity between white and red rice varieties
generate fingerprint of 19 genotypes of Oryza sativa
was studied through SSR markers. They also found
among which nine succeeded to produce
substantial difference between colored rice with
polymorphic or monomorphic alleles when applied
microsatellite markers.
with rice cultivars generating 129 bands. The level of
499
Siddhi Patel et al.
Figure 1. Dendrogram of genetic relationship among rice genotypes based on SSR markers.
A total of 21 INDEL markers were screened third minor cluster (B3) consisted of 077, BC1, Lal
out of which only 12 INDEL markers showed Kada, Krishna Kamod. Second major cluster
expected results. Highest numbers of polymorphic depicted black and brown rice, which were
bands were observed in R5M13, R8M46 and genetically closely related to red rice (Figure 2).
R9M10. The highest PIC of 0.81 was recorded in Similar work has been carried out by Yu et al.
R9M10 while lowest of 0.50 was obtained in (2002) where SSR and INDEL markers were used
R9M20. The level of polymorphism among the for diversity between several rice genotypes.
varieties was evaluated by calculating the number Xingxing et al. (2007) validated the INDEL primer
and PIC values for each of the 9 microsatellite pairs based on the comparative genomic study on
markers evaluated. The dendrogram divides the DNA sequences between several rice species. The
varieties into two major groups i.e. A and B. Line result revealed that out of 45 INDEL primers used
between the genotypes IRST2 and IRST4 indicates 41 could accurately identify indica and japonica
that they have highest genetic similarity. First rice varieties with reliability over 80%.
cluster (A) was further sub-divided into two minor All genotypes clustered at a similarity index of
clusters, first minor cluster (A1) comprised of 0.05 for SSR and 0.03 for INDEL. The similarity
IRST1, IRST2, IRST4, IRST7 and IRST10. Second index values for SSR ranged from 0.56 to 0.95
A2 comprised of IRST13 and IRST 14. Second (Table 3). Highest similarity (0.95) observed
major cluster (B) included three minor clusters, first between IRST41 and IRST44, IRST50 and TSP
minor cluster (B1) comprised IRST16, IRST19, whereas lowest similarity (0.56) obtained between
IRST41 and IRST44. Second minor cluster (B2) IRST47 and LAL KADA, IRST4 and LAL KADA.
included IRST47, IRST50, TSP, and A004. While Similarity index for INDEL ranged from 0.00 to
500
Emir. J. Food Agric. 2014. 26 (6): 497-507
http://www.ejfa.info/
0.88 (Table 4). Highest similarity (0.88) found alleles that a marker has and the frequency of each
between IRST2 and IRST4.Average similarity of the alleles in the subset of germplasm tested.
index observed was 0.74 in SSR whereas 0.17 in Since, a marker with fewer alleles has less power to
INDEL. Thus, it can be inferred that less diversity distinguish several samples, and alleles present at
was detected using SSR markers as it is evident low frequency have less power to distinguish, a
from its similarity value. higher PIC is assigned to a marker with many
Arithmetic mean heterozygosity Allelic alleles and with alleles present at roughly equal
diversity, polymorphism information content (PIC) proportions in the population (Jiang et al., 2010).
and marker indices (MI) are measures for detecting Among 9 SSR markers and 12 INDEL markers
polymorphism and comparing efficiency of highest PIC value was obtained for INDEL marker
markers. PIC is a quantification of the number of R9M10 0.81.
Figure 2. Dendrogram of genetic relationship among rice genotypes based on INDEL markers.
501
Emir. J. Food Agric. 2014. 26 (6): 497-507
doi: 10.9755/ejfa.v26i6.15722
http://www.ejfa.info/
502
Emir. J. Food Agric. 2014. 26 (6): 497-507
http://www.ejfa.info/
503
Emir. J. Food Agric. 2014. 26 (6): 497-507
doi: 10.9755/ejfa.v26i6.15722
http://www.ejfa.info/
A more specific comparison of overall heterozygosity across all loci was 0.61 and 0.60 for
efficiency of two marker systems was provided by SSR and INDEL marker respectively. All the
marker index (MI). Comparison of the marker measures used to compare the efficiency of the two
index values revealed high level of polymorphism markers indicated that there was higher level of
in INDEL (0.23) than SSR (0.19) (Tables 5, 6). polymorphism detected by INDEL than SSR
Mean allelic diversity for 9 SSR markers were markers, thus it can be concluded that INDEL
found to be 0.67 and 0.64 for 12 INDEL markers markers was more efficient than SSR in the present
which indicated that the diversity detected by study.
INDEL was more than SSR. The arithmetic mean
Heterozygosity
Individuals
Polymorph
Frequency
Frequency
Genotype
Diversity
Number
Number
Content
Alleles
Allelic
Allele
Count
Count
Locus
Info
of
of
504
Emir. J. Food Agric. 2014. 26 (6): 497-507
http://www.ejfa.info/
Heterozygosity
Polymorph
Frequency
Frequency
Frequency
Genotype
Genotype
Diversity
Content
Allelic
Count
Count
Count
Locus
Info
R3M37 1 1 A 2 1.0000 0.0000 0.0000 0.0000 A/A 1 1.0000
A 2 0.1429 A/A 1 0.1429
B 4 0.2857 B/B 2 0.2857
R5M13 7 4 0.6414 0.0000 0.6939
C 2 0.1429 C/C 1 0.1429
D 6 0.4286 D/D 3 0.4286
A 4 0.2500 A/A 2 0.2500
R6M30 8 3 B 4 0.2500 0.5547 0.0000 0.6250 B/B 2 0.2500
C 8 0.5000 C/C 4 0.5000
A 20 0.5882 A/A 10 0.5882
B 2 0.0588 B/B 1 0.0588
R7M20 17 4 0.5290 0.0000 0.5813
C 8 0.2353 C/C 4 0.2353
D 4 0.1176 D/D 2 0.1176
A 2 0.2500 A/A 1 0.2500
B 2 0.2500 B/B 1 0.2500
R8M33 4 4 0.7031 0.0000 0.7500
C 2 0.2500 C/C 1 0.2500
D 2 0.2500 D/D 1 0.2500
A 14 0.3684 A/A 7 0.3684
B 10 0.2632 B/B 5 0.2632
R8M46 19 5 C 4 0.1053 0.7085 0.0000 0.7479 C/C 2 0.1053
D 4 0.1053 D/D 2 0.1053
E 6 0.1579 E/E 3 0.1579
A 4 0.2857 A/A 2 0.2857
R9M10 7 3 B 2 0.1429 0.5015 0.0000 0.5714 B/B 1 0.1429
C 8 0.5714 C/C 4 0.5714
R9M20 14 A 2 0.0714 A/A 1 0.0714
B 10 0.3571 B/B 5 0.3571
4 0.5407 0.0000 0.6122
C 14 0.5000 C/C 7 0.5000
D 2 0.0714 D/D 1 0.0714
A 2 0.1000 A/A 1 0.1000
B 2 0.1000 B/B 1 0.1000
C 2 0.1000 C/C 1 0.1000
R10M10 10 7 D 6 0.3000 0.7978 0.0000 0.8200 D/D 3 0.3000
E 2 0.1000 E/E 1 0.1000
F 2 0.1000 F/F 1 0.1000
G 4 0.2000 G/G 2 0.2000
A 2 0.0714 A/A 1 0.0714
C 2 0.0714 C/C 1 0.0714
D 2 0.0714 D/D 1 0.0714
E 8 0.2857 E/E 4 0.2857
R10M30 14 8 0.8058 0.0000 0.8265
F 6 0.2143 F/F 3 0.2143
G 2 0.0714 G/G 1 0.0714
H 2 0.0714 H/H 1 0.0714
I 4 0.1429 I/I 2 0.1429
5 A 10 0.4167 A/A 5 0.4167
B 4 0.1667 B/B 2 0.1667
R11M23 12 C 4 0.1667 0.6990 0.0000 0.7361 C/C 2 0.1667
D 4 0.1667 D/D 2 0.1667
E 2 0.0833 E/E 1 0.0833
A 6 0.1667 A/A 3 0.1667
B 16 0.4444 B/B 8 0.4444
R12M43 18 5 C 4 0.1111 0.6859 0.0000 0.7222 C/C 2 0.1111
D 4 0.1111 D/D 2 0.1111
E 6 0.1667 E/E 3 0.1667
505
Siddhi Patel et al.
506
Emir. J. Food Agric. 2014. 26 (6): 497-507
http://www.ejfa.info/
Zeng, Y. X., Z. H. Wen, L.Y. Ma, Z. J. Ji, X. M. Li Zhou, H., Z. Xie and S. Ge. 2003. Microsatellite
and C. D. Yang. 2013. Development of 1047 analysis of genetic diversity and population
insertion deletion markers for rice genetic genetic structure of a wild rice (Oryza
studies and breeding. Genet. Mol. Res. rufipogon Griff.) in China. Theor. Appl.
12(4):5226-5235. Genet. 107:332-339.
507