Reproduced from Crop Science. Published by Crop Science Society of America. All copyrights reserved.
RESEARCH
Genetic Relatedness of Portuguese Rice Accessions from Diverse Origins as Assessed by Microsatellite Markers P. Jayamani, S. Negrão, M. Martins, B. Maçãs, and M. M. Oliveira*
ABSTRACT Simple sequence repeat (SSR) markers detect a significantly high degree of polymorphism in rice (Oryza sativa L.) and are particularly suitable for evaluating genetic diversity among closely related cultivars. A total of 176 rice accessions originating from 19 countries in the Portuguese working germplasm collection and two standard rice varieties (IR36-indica and Nipponbarejaponica) were analyzed for DNA profile using 24 SSR loci covering two loci per chromosome. A total of 184 alleles were detected. The number of alleles per locus ranged from 3 to 16, with an average of 7.7, and the PIC value ranged from 0.179 to 0.894 with an average of 0.667. All the loci were polymorphic among the accessions and clearly distinguished the indica and japonica subspecies. At 20% similarity, cluster analysis of the 178 accessions revealed three major groups, japonica, basmati, and indica (Groups I, II, and III, respectively). The japonica group contained 87% of the accessions and showed a wide range of similarity values (0.21–0.92), revealing a high degree of diversity among the accessions. Many of the accessions included in this study are morphologically similar and lack pedigree information. Hence, identification of genetic distances among the accessions should improve their use in breeding programs. As a result of this study, genetically diverse parents can be identified, increasing the usefulness of germplasm collections by broadening the genetic base of rice varieties.
P. Jayamani, S. Negrão, M. Martins, and M.M. Oliveira, ITQB/IBET, Quinta do Marquês, 2784-505 Oeiras, Portugal; P. Jayamani, CPBG, Tamil Nadu Agricultural Univ., Coimbatore 641003, India. M. Martins, Instituto Gulbenkian de Ciência, Apartado 14, P-2781-901, Oeiras, Portugal; B. Maçãs, ENMP Apartado 6, 7350-951 Elvas, Portugal; M.M. Oliveira, Univ. Lisboa, Fac. Ciências, Dep. Biologia Vegetal, 1749-016 Campo Grande, Lisboa Portugal. Received 11 Apr. 2006. *Corresponding author (
[email protected]). Abbreviations: COTArroz, Centro Operativo e Tecnológico do Arroz; EAN, Estação Agronómica Nacional; IRRI, International Rice Research Institute; NPT, new plant type; PCA, principal component analysis; PCR, polymerase chain reaction; PIC, polymorphism information content; RAPD, randomly amplified polymorphic DNA; SSR, simple sequence repeat; UPGMA, unweighted pair group mean average method.
R
ice (Oryza sativa L.) has one of the largest ex situ germplasm collections in the world (Jackson and Juggan, 1993; Ni et al., 2002), with 108 000 accessions, mostly landraces or breeding materials of O. sativa plus a small number of O. glaberrima, 22 wild species of Oryza, and a few other related genera being housed at the International Rice Research Institute (IRRI) in the Philippines (Khush, 1997; Vaughan et al., 2003). To meet the continuously expanding needs of rice varietal improvement, the assembling, evaluation, and preservation of germplasm collections is essential. In addition to phenotypic characterization, DNA markers are being used to enhance the characterization and use of germplasm resources. For rice McCouch et al. (2002) reported that 2240 microsatellite markers or simple sequence repeats (SSR) were available, and this number has now doubled (http://www.gramene.org; verified 19 Feb. 2007).
Published in Crop Sci 47:879–886 (2007). doi: 10.2135/cropsci2006.04.0236 © Crop Science Society of America 677 S. Segoe Rd., Madison, WI 53711 USA All rights reserved. No part of this periodical may be reproduced or transmitted in any form or by any means, electronic or mechanical, including photocopying, recording, or any information storage and retrieval system, without permission in writing from the publisher. Permission for printing and for reprinting the material contained herein has been obtained by the publisher.
CROP SCIENCE, VOL. 47, MARCH–APRIL 2007
879
Reproduced from Crop Science. Published by Crop Science Society of America. All copyrights reserved.
SSR markers have been effectively used to study genetic diversity among closely related rice cultivars (Yu et al., 2003; Saini et al., 2004; Spada et al., 2004; Zhu et al., 2004; Garris et al., 2005). Cultivated rice can be divided in two different subspecies, japonica and indica, that can be distinguished based on agromorphological traits. Indica rice has longer and narrower grains that tend to remain separate after cooking, while japonica has shorter and rounded grains that tend to stick together after cooking. Rice was domesticated in Asia and brought to Europe (Mediterranean region) by returning members of Alexander the Great’s expedition to India (324 bc). It was only in the 15th or 16th century, however, that rice became an established crop in the region (Khush, 1997). In Portugal, rice production was introduced by the Moors in the 10th century. Today the Portuguese consume the most rice in Europe, with 17.8 kg capita–1 yr–1, followed by the Spanish (7.1 kg capita–1 yr–1) (Vasconcelos et al., 2002). The Portuguese rice germplasm has mostly japonica varieties, which presently occupies 80% of the rice-growing area in Portugal, with the remaining area having indica-like varieties (indica grain type with japonica genetic background). In Portugal, the rice breeding program was resurrected about 3 yr ago after being discontinued in the mid-1980s, with the aim of improving traditional japonica varieties. Traditional Portuguese landraces are very tall, have low productivity, and are sensitive to diseases, especially blast (Magnaporthe oryzae B. Couch), but the landraces have good grain quality and are well adapted to Portuguese preferences. The objective of the rice breeding work is to use traditional Portuguese japonica rice varieties as the recipient and IRRI varieties (indica varieties) as donors to introduce blast resistance, increase yield through introgression of the semidwarf gene (sd1), and further improve grain quality traits. Presently, all commercial varieties are either from Italy or from overseas. The renewed interest in rice led to the formation of COTArroz (Centro Operativo e Tecnológico do Arroz) in Salvaterra de Magos, Portugal, and included a range of people from producers to those in the rice industry. The working rice germplasm collection includes 176 rice accessions that are maintained in the Department of Breeding and Genetic Resources of Estação Agronómica Nacional (EAN), Oeiras, Portugal. The objective of our study was to assess the genetic diversity of rice accessions in the working germplasm collection using microsatellite markers so that diverse parents could be selected to broaden the genetic base of Portuguese rice varieties.
MATERIALS AND METHODS Plant Materials The Portuguese germplasm collection has 1200 accessions, with several introductions made in the beginning of the 20th century. The accessions in collection were introduced from Asia and 880
other European countries, as well as from USA and the IRRI. The materials include traditional varieties of Portugal, breeding lines, cultivars, etc., but most of them lack pedigree information. A total of 176 accessions from this set, corresponding to the Portuguese working germplasm collection and originated from 19 different countries, were used in this study. The material includes indica and japonica subspecies, basmati and new plant type (NPT) accessions that are accessions recently developed at IRRI, with low tillering and high-density panicle (Table 1). Two rice varieties, IR36 and Nipponbare, recently received from the IRRI in the Philippines were used as standard varieties representing typical indica and japonica subspecies, respectively.
DNA Extraction and PCR Analysis Thirty-day-old whole seedlings of 178 accessions were collected from the field (COTArroz, Salvaterra de Magos). The seeds were sown continuously in two rows with 3-m length. Leaf samples were prepared in duplicate from mixtures of leaves taken from 10 seedlings per accession. For DNA extraction, after collection, the leaf samples were maintained in cold for up to 1 h until freezing in liquid N and preservation at –80ºC. The DNA was isolated from all the accessions using the CTAB method (Doyle and Doyle, 1987). DNA from all samples was quantified by spectrophotometry. The quality of random samples was checked by agarose gel electrophoresis. A total number of 24 SSR primer pairs (2 chromosome−1) were selected on the basis of a previous study we conducted involving seven japonica and two indica accessions and testing 165 SSRs (S. Negrão, unpublished data, 2004). The selected primers were originally mapped by Panaud et al. (1996), Chen et al. (1997) and Temnykh et al. (2000, 2001) (Table 2). The primer sequences and PCR conditions can be found in the rice database (http://www.gramene.org). Polymerase chain reaction amplification (Thermocycler Biometra UNO II; Biometra, Germany) using the SSR primers (Illumina, Inc., San Diego, CA) was conducted in a 25-μL reaction mixture containing 40 ng of template DNA, 1× PCR buffer (Invitrogen, Carlsbad, CA), 2 mM of magnesium chloride (MgCl 2), 400μM of dNTPs, 0.3 μM of each primer, and 1.5 units of Taq DNA polymerase (Invitrogen, Carlsbad, CA). The PCR conditions were set at 94ºC for 5 min, followed by 35 cycles of 94ºC for 1 min, 55ºC for 1 min, 72ºC for 2 min, and 7 min at 72ºC for the final extension. The same PCR conditions were adopted for all primers except for RM186 where the annealing temperature was set at 61ºC as recommended in the Gramene Database. The products were detected using 6 or 8% polyacrylamide gels in 1× Tris-Borate EDTA buffer according to the size of the bands. The gel system used was the Double Wide Mini-Vertical Gel unit (16 by 33 cm), from CBS Scientific (Solana Beach, CA). After electrophoresis the gels were stained with ethidium bromide solution and were photographed under ultraviolet light using Gel-Doc 1000 (Biorad, Hercules, CA). The molecular size of the amplification products was estimated using a 25-bp ladder (Invitrogen, Carlsbad, CA). The experiments were repeated to confirm cases of null alleles and to verify reproducibility.
Data Analysis Differences in molecular weight for products amplified from SSR markers were manually measured using image analysis
WWW.CROPS.ORG
CROP SCIENCE, VOL. 47, MARCH–APRIL 2007
Table 1. List of accessions used for simple sequence repeat analysis. Standard varieties are indicated in italic type. Accessions in the table are ordered according to the dendrogram (Fig. 1). Origin†/ Sl. no.
Reproduced from Crop Science. Published by Crop Science Society of America. All copyrights reserved.
CU 1
Origin†/ Sl. no.
Name V.226 C.I.700
PT 151 ‡
Origin†/ Sl. no.
Name Lido‡
IT 24 ‡
Name Sésia‡
Origin†/ Sl. no.
Name
PT GI-12 153 Safari‡ ‡
ES 20
Precoz Nº6 Italiano
PT 152
Valtego
PT 77
Harém
RU 9
M-196
PT 30
Banata 8-Estrela-A
PT 167
404 (Elvas) ‡
PT 128
Borracho‡
RU 130
D-291
PT 6
Secano‡
PT 169
Riva‡
PT 89
Super‡
RU 133
M-2233
RU 134
M-1062
‡
‡
ES 58
Pataca
ES 121
Girona
ES 56
Nano de Alcantara‡
EG 122
GZ 4120-205-2
BG 70
Plovdiv
EG GI-5 123 GZ4120-298-4-1
PT 75
Jamba Sel.1 ‡
PT 26
Almarelo‡
‡
PT 93
Diamante
PT 126
Agata‡
RU 138
M-2313
PT 127
Gritna Sel.1
US 108
Nato CI 8998
US 107
Vista CI 9628-2
PT 145
Mars
ES 21
Pataca-Nano
RU 27
Hodzaha Akhmat FR 114
Stirsia K
US 104
Blue Belle CI 9544
PT 31
Balzzaretti‡
PT 28
Arbelo Sel.1
PT 115
Topázio‡
US 106
Dawn CI 9534
IT 32
Baraggia‡
US 36
Colusa180
PT 116
Granada‡
US 105
Belle Patna CI 9433
PT 50
Greggio
PT 157
Estrela A‡
PT 117
Ribão‡
PT 154
VB1–26 VB7
‡
‡
JP 49
Fusisaka 2
PT 57
Oldenico
PT 109
Emiliano
PT 156
PT 62
Prits‡
IT 60
Precoce Rossi‡
MA 124
446
PT GI-13 155 VB36-49-3-3
PT 64
IBO 400‡
PT 81
Nora‡
PT 120
Ribelo‡
PT 25
Praganudo Rajado‡
PT 68
Redi‡
PT GI-6 82
Maioral‡
RU 136
M-2147
RU 132
M-554 M-1941
RO 69
Krasnoarmeisk 313
RO GI-1 71 L. Bannat 2951
PT 4
‡
Amarelo
PT 131
D-454 Sel 1
RU 135
PT 54
Lusitano‡
PT 147
Smeraldo‡
RU GI-14 137 M-7733
ES 7
Agostano de Espanha‡ MO 59
Precoce 6‡
PT 144
Venérea–S1‡
BR GI-15 48
PT 33
Belém‡
HU 174
Média
PT 148
Venérea-S2‡
PH GI-16 164 IR68552-100-1-2-2§
FR 34
Cesariot‡
HU GI-7 175 Rotundus
FR 139
Pygmalion
PK 40
Basmati 6129¶
FR 140
Deltaribe V
PK 41
Basmati 6131¶
PT 8
Ribatejo Sel.1
PT 72
PT 18
Precoz Nº6‡
PT 85
Onça‡
FR 141
Deltaribe A
PK 42
Basmati 6113¶
RU 10
Krasnodarsky 424
PT 95
Liso‡
FR 142
Deltaribe Y
PK 44
Basmati 5853¶
PT 143
‡
PK 45
Basmati 5836¶
PK 46
Basmati C.621¶
ES 19
Nano-Sollana
PT 76
‡
‡
Jardim
PT 97
JP GI-2 178 Nipponbare‡ PT 15 ES 39 PT GI-3 98
PT 84
Shiu-ei‡ ‡
EAS Nº2
Norinkelle ‡
RU 129
D-212
PT 149
Silla (Novo) ‡
PK 43
Basmati 5875¶
PT 73
Silla‡
IT 150
Rosa Marchetti‡
PK 47
Basmati C.622¶
PH GII 102
IR 56#
Estrela
‡
PH 99
IR 50#
Gazela‡
PT 3
Real‡
PH 100
IR 52#
Restauração‡
BG 51
Iskra
PH 103
IR 58#
PH 166
IR36 (EAN)#
PH 177
IR36 (IRRI)#
Nahda‡
PT 79
MA 13
Triomphe‡
PT 88
PT GI-4 87 Pirisca
PT 91
Patna
Livorno
MA 11
‡
MA 14 ‡
‡
PT GI-9 53
Cachaça
Judas
‡
Strella
PT 96
PT 16
Tório
PT 146
S. Domenico
‡
‡
Raffaelo
PT 80 ‡
Sesilla
‡
Dourado (DO-88)
‡
Tarado
GR 37
‡
Corbetta
PT GI-8 92
Rocca
‡
FR 38
Delta
‡
‡
PT 74
Astro‡
JP 5
Rikuto-Norin 6
GR 35
C.I. 825
PH 101
IR 54#
PT 125
Santamaro Sel. 1
PT 22
Mantova‡
IT 12
RB‡
KR 110
Suweon 281
PT 158
Muga‡
IT 23
Roma‡
PT 17
Carol‡
KR 113
Milyang 40
Ferónio
PT 63
Pancas
PT GI-10 86 Ourique‡
PH 111
IRI 339#
PT 163
EAN Nº3‡
PT 29
Baldo‡
PT 2
Preto de Pragana Algarvio‡ PH 112
IRI 340#
PT 172
Ispaniki A‡
PT 94
Rubi‡
IT 55
Maratelli‡
PH 159
IR72#
PT 162
PT 168 ES 170
‡
‡
Rio
Niva
‡
‡
PT 61
Santamaro
PT 65
‡
‡
IBO 380/33
‡
PT 52
Laboratório 7
PH 161
IR68144-2B-2-2-3-3#
PT 67
126/34
PH 160
IR58614-B-B-2-2#
PH GIII 165
IR74963-262-5 -1-3-3§
ES 171
Sénia‡
PT 78
Faísca‡
IT GI-11 173 Panda
HU 176
Sel. of Agostano
PT 90
Sagres‡
RU 118
PT GI-9 66
295/33
RU GI-12 119 Usbros 59
PT GI-5 83 Mestre‡ †
Usbros 269
Country codes as per the ISO 3166 (BG, Bulgaria; BR, Brazil; CU, Cuba; EG, Egypt; ES, Spain; FR, France; GR, Greece; HU, Hungary; IT, Italy; JP, Japan; KR, Korea; MA, Macedonia; MO, Mozambique; PH, Philippines; PK, Pakistan; PT, Portugal; RO, Romania; RU, Russia; US, USA).
‡
japonica.
§
New plant type variety (NPT).
¶
basmati.
#
indica.
CROP SCIENCE, VOL. 47, MARCH–APRIL 2007
WWW.CROPS.ORG
881
Table 2. Data on number of alleles, allele size range, number of rare alleles, number of null alleles, number of multiple alleles, highest frequency alleles and polymorphism information content (PIC) found among 178 rice accessions for 24 microsatellites. Marker
Chromosome no.
Repeat type
Annealing No. of temp. alleles
Reproduced from Crop Science. Published by Crop Science Society of America. All copyrights reserved.
ºC
Size range
No. of rare alleles
No. of No. of accessions Alleles with high PIC frequency null with multiple bands value alleles 2 bands 3 bands Size Frequency bp
%
RM84
1
(TCT)10
55
4
109–127
bp 2
2
0
0
109
90.3
0.179
RM297
1
(GA)13
55
5
148–188
0
1
130
0
180
50.8
0.665 0.735
RM341
2
(CTT)20
55
8
130–172
2
1
0
31
134
40.7
RM208
2
(CT)17
55
3
170–180
0
1
0
5
170
51.4
0.564
RM186
3
(CGG)5
61
4
116–128
1
0
0
0
120
54.5
0.545
RM227
3
(CT)10
55
5
100–124
1
4
0
0
106
55.2
0.587
RM335
4
(CTT)25
55
11
104–174
3
0
0
92
116
35.4
0.801
RM255
4
(AGG)5(AG)2-(GA)16
55
5
133–157
1
0
178
0
151
44.4
0.630
RM574
5
(GA)11
55
7
151–173
2
0
0
4
159
47.2
0.644
RM249
5
(AG)5-A2(AG)14
55
16
117–171
4
0
0
61
127
21.3
0.894
RM585
6
(TC)45
55
12
179–245
5
1
166
0
241
37.3
0.761
RM539
6
(TAT)21
55
7
266–300
2
2
178
0
272
39.8
0.753
RM481
7
(CAA)12
55
15
141–207
6
0
178
0
167
59.0
0.636
RM234
7
(CT)25
55
5
130–156
1
0
0
20
134
41.0
0.685
RM72
8
(TAT)5C(ATT)15
55
16
146–200
6
1
178
0
178
27.1
0.852
RM223
8
(CT)25
55
6
151–169
1
1
0
0
151
47.5
0.671
RM434
9
(TC)12
55
6
146–162
1
1
0
0
152
57.1
0.618
RM242
9
(CT)26
55
3
197–255
0
0
0
15
201
54.5
0.555
RM222
10
(CT)18
55
6
195–217
2
0
0
0
199
68.0
0.484
RM228
10
(CA) 6 (GA) 36
55
6
106–154
0
0
0
11
110
37.6
0.703
RM286
11
(GA)16
55
10
98–138
1
1
0
0
130
23.7
0.851
RM536
11
(CT)16
55
6
237–259
0
0
0
0
253
42.1
0.650
RM20A
12
(ATT)59
55
13
216–290
4
0
178
0
220
27.5
0.842
RM101
12
(CT)37
55
5
324–372
1
3
155
0
348
40.0
0.712
46
19 45.5
0.667
Total Mean
184 7.7
1.9
software (Mocha 1.2, Jandal Scientific, San Rafael, CA). The molecular weights for the standard varieties IR36 and Nipponbare were confi rmed with the sizes in the database (http:// www.gramene.org). The molecular size was also confi rmed by manual scoring. The polymorphism information content (PIC) for each microsatellite marker was calculated according to Nei (1973). A matrix was formed based on the allele size for each locus and accession. Two alleles per locus for SSR marker were considered. Genetic similarity was calculated with the bandsharing coefficient (Lynch 1990), using the formula: SD = 2Nab/ Na + N b, where Na and N b are the number of bands in accessions A and B; Nab, the number of bands shared by the two accessions; and SD, standard deviation. The dendrogram was constructed after cluster analysis of the similarity coefficient by means of the unweighted pair group mean average (UPGMA) method. A principal component analysis (PCA) was performed to highlight the resolving power of the ordination. These calculations were performed with a statistical software package, NTSYS-pc version 2.1 (Rohlf, 2000). One hundred resamples were used to estimate bootstrap values, using the software PowerMarker Version 3.25 (Liu and Muse, 2005; www.powermarker.net), and a consensus tree with the bootstrap values was generated in PHYLIP Version 3.5 software (Felsenstein, 1993).
882
0.8
RESULTS AND DISCUSSION SSR Analysis In the present investigation, 24 SSR markers were used to assess the genetic diversity of 178 accessions originating from 19 different countries (Table 2). This study was the first step for the characterization of the molecular diversity of a large set of rice accessions present in the Portuguese national germplasm collection. All the 24 SSR loci were polymorphic and produced a total of 184 alleles among the 178 accessions analyzed, with product sizes ranging from 98 to 372 bp. The PIC value ranged from a minimum of 0.179 (RM84) to a maximum of 0.894 (RM249), with an average of 0.667. The average number of alleles per locus was 7.7, with a range of 3 (RM208 and RM242) to 16 (RM249 and RM72). This is significantly higher than the average number of alleles reported by Zhu et al. (2004) (4.37) and Spada et al. (2004) (7.2). However, these authors used fewer accessions (85 and 96 accessions, respectively) and SSR loci (19 and 12 loci, respectively) than we used in the present work. Yu et al. (2003) studied 193 rice accessions from 26 countries with 101 SSR
WWW.CROPS.ORG
CROP SCIENCE, VOL. 47, MARCH–APRIL 2007
Reproduced from Crop Science. Published by Crop Science Society of America. All copyrights reserved.
primer pairs and detected an average allele number of 6.3 per locus, which is also lower than the value reported here. Luce et al. (2001) analyzed 419 rice accessions from gene banks in five European countries, including Portugal, with 16 SSR loci (different from the ones we selected) and reported an average of 9.1 alleles per locus. This value is higher than we observed (7.7) and probably a result of the larger number of accessions used by these authors. Still, the high level of SSR polymorphism that we detected in the Portuguese collection could be attributed primarily to the diverse origins of the japonica accessions. An allele observed in less than 5% of the 178 accessions was considered to be rare. A total of 46 rare alleles were observed. All but five loci (RM297, RM208, RM242, RM228, and RM536) exhibited one or more rare allele. In general, markers detecting a higher number of alleles per locus also detected more rare alleles. RM72 and RM481 produced the maximum of six rare alleles, with 16 and 15 alleles per locus, respectively. An accession was assigned a null allele for a locus whenever an amplification product was not observed for the particular marker. Of the 24 SSR loci used in this study, 12 had null alleles in one to four of the 178 accessions. Multiple alleles (2 or 3 alleles locus−1) were detected at one or more loci per accession even in standard varieties such as IR36 and Nipponbare. Accessions with 2 alleles locus−1 were identified when two different bands had the same intensity. Whenever the two bands had different intensities, the stronger band was considered as varietal norm. Eight loci produced two bands with different intensities (Table 2). Out of these, RM255, RM539, RM481, and RM72 had two bands with different intensities in all the accessions analyzed. RM20A produced an independent segregation of two allelic bands with equal intensity in all the accessions analyzed, which were considered as two different alleles. Multiple alleles with 3 bands locus−1 were produced by eight different SSR markers, ranging from 4 (RM574) to 92 (RM335) accessions, with the most intense band considered as the variety norm. Genetic heterogeneity is a common feature of rice accessions, despite the high inbred nature of the species (Jain et al., 2004). Most of the varieties were derived from pedigree breeding programs. Pedigree selection involves homozygous breeding; this should fix the genetic material revealing no more than 1 allele locus−1. Nonetheless, in this study, we found that nearly all accessions had two bands in one or more loci and the two bands did not segregate, indicating the accessions were not mixed. An independent segregation of two loci (two alleles) produced by RM20A may be due to the presence of duplicated regions within the rice genome. Panaud et al. (1996) observed an independent segregation of two loci (RM20A and RM20B), which were mapped on different chromosomes in a subset of doubledhaploid lines. Three alleles locus−1 (multiple alleles) were CROP SCIENCE, VOL. 47, MARCH–APRIL 2007
detected in 8 of the 24 loci. The number of multiple alleles varied among the accessions and SSR loci. Several genetic hypotheses can explain the observed polymorphism (Jain et al., 2004). The major hypotheses are residual heterozygosity, accidental seed mixtures, and mixtures of landraces. Mutation and outcrossing also may contribute to heterozygosity in the genetic accessions. In this study, intravarietal polymorphism could not be detected because DNA samples were extracted from bulk leaf samples. The average percentage of the high frequency alleles was 45.5. It ranged from 21.3 (RM249) to 90.3% (RM84). In general, a locus having a low frequency of the most common allele produced a higher number of alleles per locus. RM249 and RM72 produced the maximum number of alleles and they had the minimum frequency of a common allele, with 21.3 and 27.1%, respectively. All the frequent alleles were from japonica accessions. The most common allele at any given locus among the japonica accessions differed from the alleles in the indica and basmati accessions.
Genetic Relatedness of Accessions We used UPGMA cluster analysis based on genetic similarity values for SSR alleles from all the rice accessions to construct a dendrogram and used 100 resamplings to estimate the bootstrap values (Fig. 1). The correlation between the similarity index and cophenetic value was estimated at r = 0.867, indicating a high level of reliability. The cluster analysis showed a significant genetic variation among the rice accessions studied, with a similarity coefficient varying between 0.09 and 1.00 (Fig. 1). The dendrogram revealed two distinct groups (indica plus basmati and japonica) at the similarity coefficient of 0.09. The first group clustered with 87% of accessions from japonica subspecies; a second group included 13% of the accessions corresponding to the indica subspecies and basmati accessions (Table 1). This is in agreement with results reported by Luce et al. (2001) who described the structure of rice from the European collection as almost exclusively bipolar, with 85 and 13% of the accessions belonging to the japonica and indica subspecies, respectively. The remaining 2% aggregated in an intermediate group between indica and japonica groups. In our study, the SSRs clearly showed the well-documented indica–japonica differentiation that was reflected by subspecies-specific alleles at all the loci studied. However, the results were analyzed using 20% similarity value as the threshold for clustering where three major groups can be observed (Fig. 1). The selection of this threshold is supported by the bootstrap analysis, since the consensus tree obtained by PHYLIP, with the bootstrap values (data not shown) generated unambiguously the same three main groups that are shown in the dendrogram (Fig. 1). Group I contained mainly japonica accessions. Groups II and III represented the basmati and indica accessions, respectively.
WWW.CROPS.ORG
883
Reproduced from Crop Science. Published by Crop Science Society of America. All copyrights reserved. Figure 1. Dendrogram built from the simple sequence repeat data of 178 rice accessions of the working germplasm collection. Analyses were made by the unweighted pair group mean average method and were based on Lynch’s coefficient. Bootstrap values (out of 100) are indicated at the main branch points. Subgroups in Group I (GI-1–GI-16) although not supported by bootstrap analyses were indicated for easier correlation with the list of accessions in Table 1. Names of the rice accessions refer to those in Table 1. Most accessions in Group I are japonicas, in GII are basmati, and in GIII are indicas.
884
WWW.CROPS.ORG
Group I ( japonica) consisted of 155 accessions and had all the japonica accessions, including Nipponbare, a temperate japonica included as standard variety in the molecular analysis. The similarity coefficient among the accessions in Group I ranged from 0.21 (IR68552-100-12-2) to 0.92 (Ferónio and EAN Nº3). Dourado (DO-88) from Brazil and the new plant type accession developed in the Philippines (IR68552-100-1-2-2) were the most genetically diverse accessions in Group I, with each one clustering independently from all other accessions. The presence of this high genetic diversity in Group I may be because of the multiple introductions of genotypes from 17 different countries into the Portuguese germplasm. Genetic diversity partitioning within and among different geographic regions was observed by Yu et al. (2003). These authors found a clear distinction in the japonica group between East and South Asia, although with a reduced number of plants tested (15 and 8 accessions, respectively). In our work, ecogeographical specific differentiation was not resolved, and accessions from different origins were mixed in the japonica group. Japonica accessions originating from Portugal were also distributed throughout Group I. Although several subgroups were formed, an intermediate subgroup (indica–japonica cross derivatives) could not be identified. However, the hybrid derivatives of indica and Portuguese japonica accessions, VB7, VB1-26, and VB36-49-3-3, grouped at a low similarity level (35%) in their subgroup. The low similarity level in intermediate accessions compared with japonica accessions was also reported by Luce et al. (2001). Furthermore, Dourado from Brazil and the intermediate new plant type accession (IR68552-100-12-2) developed at IRRI showed the lowest similarity coefficients (0.28 and 0.21, respectively) and were placed between the major japonica and indica groups. This also suggests the presence of intermediate accessions in the last few subgroups of the japonica group. The higher genetic distance observed for the new plant type was expected considering its highly complex genealogy (Virk et al., 2004). Group II comprised eight basmati accessions originating from Pakistan (Table 1) and one solitary indica accession (IR56) developed in the Philippines. The similarity coefficient of this group ranged from 0.24 CROP SCIENCE, VOL. 47, MARCH–APRIL 2007
Reproduced from Crop Science. Published by Crop Science Society of America. All copyrights reserved.
previously clustered with the basmati group (Fig. 1), it had to 0.92, with Basmati 6129 and Basmati 6131 having the a low similarity value (0.24). In fact, IR56 was expected highest value. This result supports the concept that basto group with indica accessions, since it does not have a mati rice had a long independent evolution and is genetibasmati genetic background in its pedigree. The reason for cally distinct from other groups within O. sativa, namely the separation may be due to mislabeling or mixtures with indica and japonica. Based on isozyme analysis, Glaszmann basmati seeds. The remaining 14 indica accessions origi(1987) reported that the basmati genotypes are genetically nated from the Philippines and Korea formed a separate distinct from the other groups. The same was observed cluster identified as Group III in the dendrogram, except by Jain et al. (2004) and Saini et al. (2004) while studying for IR52, which was separated. In fact, IR52 was found to Indian basmati varieties using SSR markers. The acceshave a few different alleles from those of the other indica sion IR56 had the lowest similarity coefficient (0.24) comaccessions studied here. pared with other accessions in this group and was situated The correct determination of genetic diversity of between the basmati and indica groups (Fig. 2). varieties in the germplasm is important for rice breedGroup III consisted of 14 indica accessions clustering at ing programs, allowing selection of the desired rice accesthe similarity level of 0.39. In this group, 12 indica accessions for crossing. This will maximize the probability of sions were developed at IRRI, including IR74963-262-5transgressive segregation and increase the probability that 1-3-3 (NPT) and two entries of IR36 variety. One entry unrelated accessions contribute with positive alleles at difof IR 36 was available in the collection and the other one ferent loci. Usually, it is difficult to accurately identify was recently obtained from IRRI and used as standard varieties in the indica and japonica subspecies by their morvariety. These two accessions of IR 36 showed a 100% phological characteristics; however the SSR markers used similarity. The other two accessions in this group were in this study detected a high level of polymorphism and Milyang 40 and Suweon 281, originating from Korea. were successful in distinguishing the indica and japonica Varieties like Milyang 23 or Suweon 213, considered a accessions. Three major groups were identified based on Tongil-type, seem to have strong similarities to indicas. the cluster analysis: (i) the japonica group (GI), which was From a molecular study using SSRs and RAPDs, Kwon et the largest with 87% of the accessions and exhibited great al. (1999) considered the rice accessions of the Tongil-type diversity due to the diverse origins of japonica accessions, (obtained from a crossing scheme indica/japonica/indica) as (ii) the basmati accessions (GII), and (iii) the indica accesindicas. Our results with Milyang 40 and Suweon 281 are sions (GIII). Many of the accessions included in the study consistent with the data of Kwon et al. (1999). are morphologically similar and lack the pedigree inforThe groupings identified by PCA analysis (Fig. 2) were mation. Thus, the identification of the genetic distance comparable to those identified by the UPGMA tree clusamong the accessions will be important to maximize their ter analysis (Fig. 1). The first and second principal cooruse in breeding programs. Furthermore, the assessment of dinates, namely, Dim-1 and Dim-2, account for 3.23 and genetic diversity of rice accessions present in the working 2.67% respectively, of the total variation in the SSR data. germplasm collection will help the breeders to formulate Three major clusters were formed by japonica, basmati, and indica genetic accessions. Unlike the UPGMA analysis, where all accessions were assigned to a group, Ribatejo Sel.1, Rikuto-Norin 6, Riva, Sénia, Ribelo, IR56, and IR52 were out of the major clusters and appeared to be distinct from other accessions in the PCA, possibly indicating genetic differentiation. The first cluster included 150 of the 155 japonica accessions of Group I (Fig. 1). A second cluster observed in the PCA had eight basmati accessions from Pakistan (Fig. 2) as in the dendrogram (Group II). The accession IR56 was separated and placed between basmati and indica groups in the dendrogram, which is consis- Figure 2. Two-dimensional scaling resulting from principal component analysis (NTSYS-pc) of tent with the result obtained in the 178 accessions using genetic diversity data for 184 alleles at 24 simple sequence response PCA. Although this accession has loci. Circles indicate the three major clusters formed. CROP SCIENCE, VOL. 47, MARCH–APRIL 2007
WWW.CROPS.ORG
885
crosses by choosing accessions with different genetic backgrounds and will assist in the development of gene-mapping populations with greater marker polymorphism.
Reproduced from Crop Science. Published by Crop Science Society of America. All copyrights reserved.
Acknowledgments We gratefully acknowledge financial support from the Fundação para a Ciência e a Tecnologia, Lisboa, Portugal, through postdoctoral and Ph.D. research fellowships to P. Jayamani (SFRH/ BPD/14542/2003) and S. Negrão (SFRH/BD/10613/2002), respectively. We thank Isabel Marques (IGC) and Carla Ribeiro (EFN) for their help in the statistical analysis.We also thank scientists from the Departamento de Recursos Genéticos e Melhoramento, Estação Agronómica Nacional, Portugal, and COTArroz for providing plant materials and making available field conditions. IRRI is also acknowledged for providing seed material.We thank Dr. Dave Mackill (IRRI) for critical suggestions to improve this manuscript and Dr. Margarida O. Krause for revising the English. Finally, the reviewers of our paper are gratefully acknowledged for their comments.
References Chen, X., S. Temnykh, Y. Xu, G. Cho, and S.R. McCouch. 1997. Development of a micosatellite framework map providing genome wide coverage in rice. Theor. Appl. Genet. 95:553–567. Doyle, J.J., and J.L. Doyle. 1987. A rapid DNA isolation procedure for small quantities of fresh leaf tissue. Phytochem. Bull. 19:1–15. Felsenstein, J. 1993. PHYLIP (Phylogenetic Inference Package) version 3.5c. Dep. of Genetics, Univ. of Washington, Seattle. Garris, A.J., T.H. Tai, J. Coburn, S. Kresovich, and S.R. McCouch. 2005. Genetic structure and diversity in Oryza sativa L. Genetics 169:1631–1638. Glaszmann, J.C. 1987. Isozymes and classification of Asian rice varieties. Theor. Appl. Genet. 74:21–30. Jackson, M.T., and R. Juggan. 1993. Sharing the diversity of rice to feed the world. Diversity 9:22–25. Jain, S., R.K. Jain, and S.R. McCouch. 2004. Genetic analysis of Indian aromatic and quality rice (Oryza sativa L.) germplasm using panels of fluorescently labelled microsatellite markers. Theor. Appl. Genet. 109:965–977. Khush, G.S. 1997. Origin, dispersal, cultivation, and variation of rice. Plant Mol. Biol. 35:25–34. Kwon, S.J., S.N. Ahn, H.C. Hong, H.C. Choi, and H.P. Moon. 1999. Development of PCR markers for indica and japonica differentiation in Oryza sativa L. Rice Genet. Newsl. 16:21–23. Liu, K., and S.V. Muse. 2005. PowerMarker: An integrated analysis environment for genetic marker analysis. Bioinformatics 21:2128–2129. Luce, C., J.L. Noyer, D. Tharreau, N. Ahmadi, and H. Feyt. 2001. The use of microsatellite markers to examine the diversity of the genetic resources of rice (Oryza sativa) adapted to European conditions. Acta Hortic. 546:221–235.
886
Lynch, M. 1990. The similarity index and DNA fi ngerprinting. Mol. Biol. Evol. 7:478–484. McCouch, S.R., L. Teytelman, Y. Xu, K. Lbos, K. Clare, M. Walton, B. Fu, R. Maghirang, Z. Li, Y. Xing, Q. Zhang, I. Kono, M. Yano, R. Fjellstrom, G. DeClerck, D. Schneider, S. Cartinhour, D. Ware, and L. Stein. 2002. Development and mapping of 2240 new SSR markers for rice (Oryza sativa L.). DNA Res. 9:199–207. Nei, M. 1973. The analysis of gene diversity in subdivided populations. Proc. Natl. Acad. Sci. USA 70:3321–3323. Ni, J., P.M. Colowit, and D.J. Mackill. 2002. Evaluation of genetic diversity in rice sub species using microsatellite markers. Crop Sci. 42:601–607. Panaud, O., X. Chen, and S.R. McCouch. 1996. Development of microsatellite markers and characterisation of simple sequence length polymorphism (SSLP) in rice (Oryza sativa L.). Mol. Gen. Genet. 252:597–607. Rohlf, F.J. 2000. NTSYS-pc: Numerical taxonomy and multivariate analysis system, version 2.1. Exeter Software, Setauket, NY. Saini, N., N. Jain, S. Jain, and R.K. Jain. 2004. Assessment of genetic diversity within and among Basmati and non Basmati rice varieties using AFLP, ISSR, and SSR markers. Euphytica 140:133–146. Spada, A., R. Mantegazza, M. Biloni, E. Capoali, and F. Sala. 2004. Italian rice varieties: Historical data, molecular markers, and pedigrees to reveal their genetic relationships. Plant Breed. 123:105–111. Temnykh, S., G. De Clerck, A. Lukashova, L. Lipovich, S. Cartinhour, and S.R. McCouch. 2001. Computational and experimental analysis of microsatellites in rice (Oryza sativa L.): Frequency, length variation, transposon associations and genetic marker potential. Genome Res. 11:1441–1452. Temnykh, S., W.D. Park, N. Ayres, S. Cartinhour, N. Hauck, L. Lipovich, Y.G. Cho, T. Ishii, and S.R. McCouch. 2000. Mapping and genome organisation of microsatellite sequences in rice. Theor. Appl. Genet. 100:697–712. Vasconcelos, M., A. Farinha, S. Negrão, M. Romera, C. Ricardo, M. Pagés, and S.K. Datta. 2002. Rice research produced in Portugal: Present and perspectives. p. 14–15. Proc. Conf. Current European Res. on Rice, Torino, Italy. 6–8 June 2002. Fondazione per le Biotechnologie, Torino, Italy. Vaughan, D.A., H. Morishima, and K. Kadowaki. 2003. Diversity in the Oryza genus. Curr. Opin. Plant Biol. 6:139–146. Virk, P.S., G.S. Kush, and S. Peng. 2004. Breeding to enhance yield potential of rice at IRRI: The ideotype approach. Int. Rice Res. Notes 29:5–9. Yu, S.B., W.J. Xu, C.H. Vijayakumar, J. Ali, B.Y. Fu, J.L. Xu, Y.Z. Jiang, R. Marghirang, J. Domingo, C. Aquino, S. Virmani, and Z.K. Li. 2003. Molecular diversity and multilocus organization of the parental lines used in the International Rice Molecular Breeding Program. Theor. Appl. Genet. 108:131–140. Zhu, M., Y. Wang, Y. Zhu, and B.R. Lu. 2004. Estimating genetic diversity of rice land races from Yunnan by SSR assay and its implication for conservation. Acta Bot. Sin. 46:1458–1467.
WWW.CROPS.ORG
CROP SCIENCE, VOL. 47, MARCH–APRIL 2007