Simple Sequence Repeat (SSR) - Based Diversity Analysis of Groundnut (L.) Germplasm Resistant To Bacterial Wilt
Simple Sequence Repeat (SSR) - Based Diversity Analysis of Groundnut (L.) Germplasm Resistant To Bacterial Wilt
Abstract
Groundnut is one of the most important oilseed crops in the world. Bacterial wilt, caused by
Ralstonia solanacearum E. F. Smith, is one of the major biotic constraints to groundnut pro-
duction particularly in South-East Asia and East Africa. Several sources of resistance to bacterial
wilt have been identified through field screening of groundnut germplasm. The aim of the pre-
sent study was to quantify the genetic diversity among selected bacterial wilt-resistant lines, in
comparison with the levels of variation observable within the cultivated A. hypogaea gene
pool. Thirty-two SSR markers were used to assess the degree of molecular polymorphism
between 46 selected genotypes revealing 107 alleles, of which 101 (99.4%) were polymorphic
with gene diversity scores ranging from 0.103 to 0.669, averaging 0.386. Cluster and multidi-
mensional scaling analysis revealed two distinct groups within the germplasm broadly corre-
sponding to the two subspecies (hypogaea and fastigiata) of A. hypogaea. However,
accessions of varieties peruviana and aequatoriana grouped together with the varieties
from subsp. hypogaea, rather than grouping with the other varieties of subsp. fastigiata. Anal-
ysis of molecular variance (AMOVA) revealed that 15% of the total observed variation was
accounted for by disease response groups. This analysis will be useful in the selection of par-
ental genotypes for mapping populations and breeding programmes attempting to broaden
the genetic base of future groundnut cultivars. In particular, this opens up significant opportu-
nities for the development of intraspecific mapping populations that will be highly relevant to
modern groundnut breeding programmes.
production is crushed for edible oil, 32% for confection- for groundnut production should the pathogen overcome
ery consumption and the remaining 15% is used for feed these sources of resistance. In addition, the currently
and seed production. Cultivated groundnut consists of available BW-resistant cultivars are low yielding and
two subspecies, hypogaea and fastigiata, which are have poor tolerance to drought. Identification and utiliz-
further classified into six botanical varieties largely ation of a broad spectrum of genetically diverse sources
based on growth habit, flowering and branching of BW resistance is, therefore, critical for the develop-
patterns, presence of hairs on leaf surfaces and numbers ment of a new generation of broad-based high-yielding
of seeds per pod (Krapovickas and Gregory, 1994) and BW-resistant groundnut cultivars. Limited knowledge
with lesser support by protein and genomic analysis about the genetic diversity of the BW-resistant germplasm
(Smartt and Simmonds, 1995), with a number of recent and deleterious linkage drag have impeded the utilization
genomic studies refuting this classification (He and of a wide spectrum of BW resistance donors.
Prakash, 2001; Raina et al., 2001; Ferguson et al., 2004a; Diversity studies in groundnut have generally revealed
He et al., 2005; Tallury et al., 2005). Botanical varieties extensive phenotypic variation amongst varieties (Upad-
hypogaea (Virginia) and hirsuta (Peruvian) are currently hyaya et al., 2001, 2003) yet limited variation at the mol-
classified under subsp. hypogaea while varieties fasti- ecular level (Halward et al., 1991, 1992; Kochert et al.,
giata (Valencia), peruviana, aequatoriana and vulgaris 1991, Paik-Rao et al., 1992; He and Prakash, 1997;
(Spanish) are assigned within subsp. fastigiata. Subramanian et al., 2000; Moretzsohn et al., 2004). It is
Bacterial wilt (BW), caused by Ralstonia solanacearum hypothesized that this may be due to the selective neu-
E. F. Smith, is a major biotic factor affecting groundnut trality of the molecular markers utilized, while phenoty-
production particularly in South-East Asia and East pic traits have been subjected to intense selection
Africa (Hayward, 1990). It also infects many other crop (He and Prakash, 1997). It has also been suggested
plants including potato (Solanum tuberosum L.), tomato that the lack of molecular polymorphism revealed to
(Lycopersicon esculentum Mill), tobacco (Nicotiana date within the cultivated groundnut gene pool could
spp.), pepper (Capsicum spp.), eggplant (Solanum mel- be due to the inadequacy of the material studied and
ongina L.) and ginger (Zingiber officinale Rosc.). The the range of techniques used (Singh et al., 1998).
bacterial species have been isolated and classified into Microsatellite markers, also known as simple sequence
five races based on host range (Buddenhagen and repeat (SSR) markers, have been reported to detect high
Kelman, 1964; He et al., 1983) and five biovars based levels of polymorphism even amongst closely related cul-
on biochemical characteristics (Hayward, 1964; He et al., tivated germplasm (Gianfranceschi et al., 1998). For this
1983). Race 1 causes wilt in groundnut, in addition to reason, considerable efforts have been recently made to
many other leguminous and solanaceous plants. Biovar develop a large number of SSR markers in groundnut
1 causes wilt in groundnut and predominantly occurs in (Hopkins et al., 1999; He et al., 2003, 2005; Ferguson
America, whereas Biovars 3 and 4 cause wilt of ground- et al., 2004b; Moretzsohn et al., 2004, 2005). It has been
nut in Asia and Africa (Hayward, 1991). demonstrated in previous studies that SSR markers are
Extensive screening of groundnut germplasm, largely more variable within genomes than other marker types
based on field evaluations at disease ‘hot-spots’ in Indo- (e.g. Belaj et al., 2003). Additionally, SSRs have the
nesia and China, has resulted in the identification of advantage of being co-dominant, only requiring very
many BW-resistant lines. Most of these sources of BW small amounts of DNA and hence have been widely
resistance originate from China or Indonesia. Up to a applied in many plant genetics studies, e.g. for evaluating
quarter of germplasm accessions have shown some genetic diversity (Zhebentyayaeva et al., 2003; Fahima
level of resistance to BW (Singh et al., 1997; Pande et al., 1998), genome mapping and gene tagging, e.g.
et al., 1998); resistance has been identified across differ- in rice (Chen et al., 1997), wheat (Röder et al., 1998),
ent botanical types and also in some wild species. barley (e.g. Künzel and Waugh, 2002) and tomato
Despite the broad range of BW-resistant donor geno- (Broun and Tanksley, 1996). The recent development of
types, only a few of the resistant landraces have been groundnut-specific SSRs (Hopkins et al., 1999; He et al.,
successfully used in breeding in China and Indonesia 2003, 2005; Ferguson et al., 2004b; Moretzsohn et al.,
(Liao et al., 1998). The majority of BW-resistant cultivars 2004, 2005) now offers new and exciting opportunities
released in China were derived from just three sources for groundnut genomics.
(Xiekangqing, Taishan Sanlirou or Taishan Zhenzhu) In the present study we have used groundnut-specific
(Liang, 1998; Liao et al., 1998), while a single resistance SSRs to analyse a diverse range of cultivated groundnut
donor source (Schwartz 21) has been the basis of the accessions encompassing all six botanical varieties. The
majority of cultivars released in Indonesia. Thus, there purpose of this study was to investigate the level of mol-
is a progressive narrowing of genetic diversity in BW- ecular polymorphism amongst BW-resistant accessions
resistant breeding programmes that creates a major risk and to compare this with the genetic diversity across
SSR-based diversity analysis of groundnut 29
the cultivated A. hypogaea gene pool. This analysis is (2 18C/cycle) for 1 min and 728C for 1 min 30 s. This
important for the selection of genetically diverse parental was then followed by 20 cycles of 948C for 45 s, 558C
genotypes for mapping populations and BW resistance for 1 min and 728C for 1 min 30 s, followed by a final
breeding programmes aimed at the development of extension step of 10 min (728C).
broad-based cultivars with durable disease resistance.
Electrophoresis and data collection
(pPGPseq-8D9) to 531 bp (A1-193). The GD scores of the termed an amphidiploid (Burow et al., 2001). This has
29 polymorphic SSR loci ranged from 0.103 (pPGPseq- significant implications in the appropriate choice of bio-
13E9) to 0.669 (pPGPseq-3A08) (Table 2). metric analysis. In particular, SSRs may not always
retain their codominant nature as defining allelic relation-
ships becomes difficult when three or more alleles are
Genetic diversity analysis of cultivated groundnut detected in a single individual, however, this does not
germplasm occur in any of the polymorphic loci scored in this
study, with the exception of A1-275. For this reason clus-
Groundnut is a complex polyploid with two distinct gen- ter analysis was undertaken based on two different simi-
omes that largely segregate in isolation and is, therefore, larity measures: Jaccard (for dominant datasets) and
SSR-based diversity analysis of groundnut
Table 2. Details of SSRs used in the groundnut diversity analysis (SSRs were developed by Ferguson et al., 2004a unless indicated otherwise)
Expected
product Total no. of alleles Gene diversity
Primer name Repeat motif Tm (8C) size (bp) Observed product size(s) (bp) observed (GD)
pPGPseq-1B09 GA 64 282 268; 269 2 0.306
pPGPseq-2B10 TAA 58 259 268; 269 2 0.297
pPGPseq-2D12B TAA 60 265 289; 300; 323; 333 4 0.57
pPGPseq-2E06 GA 60 250 269; 289; 306; 323 4 0.504
pPGPseq-2G03 TAA 64 215 254; 269; 281 3 0.459
pPGPseq-2G04 TAA 60 289 269; 289; 300; 333 4 0.47
pPGPseq-3A01 TAA 64 238 257; 269; 277; 289 4 0.304
pPGPseq-3A08 TAA 64 152 173; 178; 191; 197 4 0.669
pPGPseq-3B06 GA 61 244 157 1 0
pPGPseq-3B08 TAA 56 266 289; 300; 314; 323 4 0.555
pPGPseq-3D09 GA,GT 63 292 281 1 0
pPGPseq-4A06 AT 63 126 167; 173; 177 3 0.402
pPGPseq-7G2 TATC 65 225 223; 239; 246; 250; 257; 262 6 0.61
pPGPseq-7H6 CTT 60 300 308; 310 2 0.297
pPGPseq-8D9 CTT 61 132 131; 135; 146 3 0.468
pPGPseq-8E12 TTG,TAA 59 198 204; 207; 210; 214 4 0.485
pPGPseq-10H1A CTT 58 139 193; 200 2 0.296
pPGPseq-12F7 TAA 57 290 305; 310 2 0.375
pPGPseq-13E9 TAA 59 299 323; 333 2 0.103
pPGPseq-13A7 TAA 58 265 289; 291; 293 3 0.44
pPGPseq-14A7 CTT,CTG 60 173 167; 173; 177 3 0.402
pPGPseq-14F4 TAA 60 163 173; 178; 184 3 0.402
pPGPseq-14H6 GT 59 285 269; 271; 284; 297; 306; 310; 312; 319; 333; 348 10 0.618
pPGPseq-15C10 TAA 64 203 212; 220 2 0.427
pPGPseq-16G8 TAA 60 194 214; 217; 227; 229 4 0.32
pPGPseq-18A5A AT,TAA 60 268 300; 328 2 0.468
Lec-1a AT 65– 55 120,125 243; 250; 261; 281; 300 5 0.473
Ah4-26a CT 65– 55 160 173; 178; 184 3 0.616
A1-041b Unknown 65– 55 230,350 269; 281; 293 3 0.351
A1-193b Unknown 65– 55 460 510; 520; 531 3 0.403
A1-275b Unknown 65– 55 190,300 181; 195; 305; 330 4 0
A1-745b Unknown 65– 55 150,250 224; 226; 236 3 0.57
a
Microsatellites developed by Hopkins et al. (1999).
b
Microsatellites developed by Moretzsohn et al. (2004).
31
32 E. S. Mace et al.
six botanical varieties within the A. hypogaea gene pool. A There was some degree of clustering of accessions
clear distinction was observed between the two subspe- from similar geographic origins within subspecies, e.g.
cies, subsp. hypogaea and subsp. fastigiata. Accessions of four landraces and one breeding line from South and
var. hypogaea and hirsuta (subsp. hypogaea) grouped Central America grouped together at a level of approxi-
together in cluster B (Fig. 2), whereas var. vulgaris and fas- mately 50% similarity in cluster B and were clearly differ-
tigiata (subsp. fastigiata) grouped together in cluster entiated from accessions derived from other geographic
A. However, accessions of var. peruviana and var. aequa- regions. Additionally, the BW-resistant accessions Bobai
toriana grouped together with the varieties from subsp. Shiyaodou and Bobai Dahuasheng both originate from
hypogaea in cluster B, rather than grouping with the Bobai County in China and were observed to group
other varieties of subsp. fastigiata in cluster A. These results together at a level of approximately 82% similarity in clus-
support recent conclusions from amplified fragment length ter B. However, AMOVA estimated 86% of the SSR vari-
polymorphism (AFLP) analysis (He and Prakash, 2001) that ation is accounted for within countries. So although
var. aequatoriana and peruviana are much closer to subsp. there may be some influence of geographic isolation on
hypogaea than to subsp. fastigiata. genetic polymorphism, as has been reported previously
Most previous reports of diversity analysis in cultivated (e.g. He and Prakash, 2001), it is likely that breeding
groundnut have readily detected morphological variation selection pressure has had greater impact. However, it
but consistently failed to detect a parallel level of molecu- should be noted that the groundnut accessions studied
lar genetic variability, however recent studies (Ferguson here were selected primarily based on their response to
et al., 2004a; Moretzsohn et al., 2004, 2005) using more BW and thus may not be truly reflective of groundnut
recently developed SSRs have found higher levels of gen- variability in each geographical region.
etic diversity in the six botanical varieties of A. hypogaea Regarding the genetic variation observed within and
than reported previously. Significantly, the SSR screening between the BW-resistant and -susceptible accessions,
reported in this study also detected a substantial level of it was observed that the BW-resistant lines belonging
molecular genetic variation between genotypes of all to subsp. fastigiata and subsp. hypogaea were clearly
botanical varieties. However, AMOVA showed an equal differentiated, and clustered within subspecies. In con-
level of diversity within and between botanical varieties. trast, two BW-susceptible accessions (Zao18 and Zhon-
This may suggest that botanical variety designations gua No. 5) did not fall in the expected cluster (subsp.
(largely based on morphological traits) are not truly fastigiata var. vulgaris), but were grouped in cluster
reflective of gross genetic diversity. B (predominately subsp. hypogaea). This anomaly
A
II
III
II
III
Fig. 2. Dendrogram constructed using Jaccard’s similarity coefficient and UPGMA clustering, for the 46 groundnut geno-
types. Two main clusters (A and B) and subclusters are identified.
34 E. S. Mace et al.
Table 3. AMOVA (Excoffier et al., 1992) for (A) 39 genotypes of two disease response types, bacterial wilt-resistant
and -susceptible, employing 107 SSR alleles; and for (B) 46 genotypes in six botanical varieties employing 107
microsatellite alleles (nested analysis was carried out on populations grouped as above)
Kruskal JB and Wish M (1978) Multidimensional Scaling. cation, and phylogenetic relationships in peanut (Arachis
Newbury Park, CA: Sage. hypogaea) cultivars and wild species. Genome 44: 763– 772.
Künzel G and Waugh R (2002) Integration of microsatellite mar- Röder MS, Korzun V, Wendehake K, Plaschke J, Tixier M-H,
kers into the translocation-based physical RFLP map of Leroy P and Ganal MW (1998) A microsatellite map of
barley chromosome 3H. Theoretical and Applied Genetics wheat. Genetics 149: 2007 –2023.
105: 660–665. Rohlf FJ (2001) Numerical Taxonomy and Multivariate Analysis
Liang XQ (1998) Present status of groundnut bacterial wilt in System. New York: Applied Biostatistics.
south China. In: Pande S, Liao BS, Hong NX, Johansen C Schneider S, Roessli D and Excoffier L (2000) Arlequin: A Soft-
and Gowda CLL (eds) Groundnut Bacterial Wilt in Asia. ware for Population Genetics Data Analysis. version
Proceedings of the Fourth Working Group Meeting. 11–13 2.000. Geneva: Genetics and Biometry Laboratory, Depart-
May 1998. Hanoi, Vietnam: Vietnam Agricultural Science ment of Anthropology, University of Geneva.
Institute, pp. 23– 26. Simpson CE (1991) Pathways for introgression of pest resistance
Liao B, Duan N, Jiang H, Liang X and Gao G (1998) Germplasm into Arachis hypogaea L. Peanut Science 18: 22 –26.
screening and breeding for resistance to bacterial wilt in Simpson CE, Nelson SC, Starr JL, Woodard KE and Smith OD
China. In: Pande S, Liao BS, Hong NX, Johansen C and (1993) Registration of TxAG-6 and TxAG-7 peanut germ-
Gowda CLL (eds) Groundnut Bacterial Wilt in Asia. Pro- plasm lines. Crop Science 33: 1418.
ceedings of the Fourth Working Group Meeting. 11– 13 Singh AK, Mehan VK and Nigam SN (1997) Sources of Resistance
May 1998. Hanoi, Vietnam: Vietnam Agricultural Science to Groundnut Fungal and Bacterial Diseases, an Update
Institute, pp. 75– 81. and Appraisal. ICRISAT Bulletin No. 50. Patancheru,
Mace ES, Buhariwalla HK and Crouch JH (2003) A high India: ICRISAT.
throughput DNA extraction protocol for tropical molecular Singh AK, Smartt J, Simpson CE and Raina SN (1998) Genetic
breeding programs. Plant Molecular Biology Reporter 21: variation vis-à-vis molecular polymorphism in groundnut,
459a – 459h. Arachis hypogaea L. Genetic Resources and Crop Evolution
Mohan M, Nair S, Bhagwat A, Krishna TG, Yano M, Bhatia CR 45: 119– 126.
and Sasaki T (1997) Genome mapping, molecular markers Smartt J and Simmonds NW (1995) Evolution of Crop Plants.
and marker-assisted selection in crop plants. Molecular Harlow: Longman Scientific and Technical.
Breeding 3: 87– 103. Sokal RR and Michener CD (1958) A statistical method for eval-
Moretzsohn MC, Hopkins MS, Mitchell SE, Kresovich S, Valls JFM uating systematic relationships. University of Kansas Scien-
and Ferreira ME (2004) Genetic diversity of peanut (Ara- tific Bulletin 38: 1409 –1438.
chis hypogaea L.) and its wild relatives based on the anal- Subramanian V, Gurtu S, Nageswara Rao RC and Nigam SN
ysis of hypervariable regions of the genome. BMC Plant (2000) Identification of DNA polymorphism in cultivated
Biology 4: 11. groundnut using random amplified polymorphic DNA
Moretzsohn MC, Leoi L, Proite K, Guimãres PM, Leal-Bertioli (RAPD) assay. Genome 43: 656– 660.
SCM, Gimenes MA, Martins WS, Valls JFM, Grattapaglia D Tallury SP, Hilu KW, Milla SR, Friend SA, Alsaghir M, Stalker
and Bertoli DJ (2005) A microsatellite-based, gene-rich HT and Quandt D (2005) Genomic affinities in Arachis
linkage map for the AA genome of Arachis (Fabaceae). section Arachis (Fabaceae): molecular and cytogenetic
Theoretical and Applied Genetics 111: 1060 –1071. evidence. Theoretical and Applied Genetics 111:
Nei M and Li WH (1979) Mathematical model for studying genetic 1229 – 1237.
variation in terms of restriction end nucleases. Proceedings of Upadhyaya HD, Ferguson ME and Bramel PJ (2001) Status of the
the National Academy of Sciences USA 76: 5269– 5273. Arachis germplasm collection at ICRISAT. Peanut Science
Paik-Rao OG, Smith RL and Knauft DA (1992) Restriction frag- 28: 89– 96.
ment length polymorphism evaluation of six peanut Upadhyaya HD, Ortiz R, Bramel PJ and Singh S (2003) Develop-
species within the Arachis section. Theoretical and Applied ment of a groundnut core collection using taxonomical,
Genetics 84: 201– 208. geographical and morphological descriptors. Genetic
Pande S; Liao BS; Hong NX; Johansen C; Gowda CLL, (eds) Resources and Crop Evolution 50: 139–148.
(1998) Groundnut Bacterial Wilt in Asia. Proceedings of Weir BS (1990) Genetic Data Analysis. Methods for Discrete Gen-
the Fourth Working Group Meeting. 11 –13 May 1998. etic Data. Sunderland, MA: Sinauer Associates.
Hanoi, Vietnam: Vietnam Agricultural Science Institute. Zhebentyayaeva TN, Reighard GL, Gorina V and Abbott AG
Raina SN, Kojima T, Ogihara Y, Singh KP and Devarumath RM (2003) Simple sequence repeat (SSR) analysis for assess-
(2001) RAPD and ISSR fingerprints as useful genetic ment of genetic variability in apricot germplasm. Theoreti-
markers for analysis of genetic diversity, varietal identifi- cal and Applied Genetics 106: 435–444.