Cultivated peanut or groundnut (Arachis hypogaea L.) is an important oilseed crop with an allotetraploid genome (AABB, 2n = 4x = 40). Both the low level of genetic variation within the cultivated gene pool and its polyploid nature limit the utilization of molecular markers to explore genome structure and facilitate genetic improvement. Nevertheless, a wealth of genetic diversity exists in diploid Arachis species (2n = 2x = 20), which represent a valuable gene pool for cultivated peanut improvement. Interspecific populations have been used widely for genetic mapping in diploid species of Arachis. However, an intraspecific mapping strategy was essential to detect chromosomal rearrangements among species that could be obscured by mapping in interspecific populations. To develop intraspecific reference linkage maps and gain insights into karyotypic evolution within the genus, we comparatively mapped the A- and B-genome diploid species using intraspecific F2 populations. Exploring genome organization among diploid peanut species by comparative mapping will enhance our understanding of the cultivated tetraploid peanut genome. Moreover, new sources of molecular markers that are highly transferable between species and developed from expressed genes will be required to construct saturated genetic maps for peanut.
A total of 2,138 EST-SSR (expressed sequence tag-simple sequence repeat) markers were developed by mining a tetraploid peanut EST assembly including 101,132 unigenes (37,916 contigs and 63,216 singletons) derived from 70,771 long-read (Sanger) and 270,957 short-read (454) sequences. A set of 97 SSR markers were also developed by mining 9,517 genomic survey sequences of Arachis. An SSR-based intraspecific linkage map was constructed using an F2 population derived from a cross between K 9484 (PI 298639) and GKBSPSc 30081 (PI 468327) in the B-genome species A. batizocoi. A high degree of macrosynteny was observed when comparing the homoeologous linkage groups between A (A. duranensis) and B (A. batizocoi) genomes. Comparison of the A- and B-genome genetic linkage maps also showed a total of five inversions and one major reciprocal translocation between two pairs of chromosomes under our current mapping resolution.
Our findings will contribute to understanding tetraploid peanut genome origin and evolution and eventually promote its genetic improvement. The newly developed EST-SSR markers will enrich current molecular marker resources in peanut.
Peanut or groundnut (Arachis hypogaea L.) is both an important oilseed crop and a direct protein source for human nutrition and is the only domesticated species in the genus Arachis. It is an allotetraploid (2n = 4x = 40, AABB) with one pair of distinctively small chromosomes
, and was believed to have originated from a single hybridization event at least 3,500 years ago
[2–4]. This genus also contains additional tetraploid, diploid (2n = 2x = 20) and aneuploid (2n = 2x = 18) species. Arachis hypogaea is a member of section Arachis, which contains species with A, B, or D genomes. The Arachis genome is characterized by the presence of a small pair of chromosomes with a lower level of euchromatin condensation
, whereas the B genome is distinguished by the presence of a pair of chromosomes with a secondary constriction proximal to the centromere
, and has recently been divided into subgroups
. Arachis duranensis and A. ipaënsis are most likely the ancestral A- and B-genome species of cultivated peanut, respectively
[2, 5, 8–10]. The D genome is represented by the single species A. glandulifera[11, 12].
Due to both the low level of genetic variation within the cultivated gene pool and the polyploid nature of cultivated peanut, relatively fewer numbers of genetic linkage maps have been constructed as compared to many other crop species
. Considering its relatively large genome size (2800Mb/1C), there is still great need to improve marker density and genome coverage
[14, 15]. For example, mapping of quantitative trait loci (QTL) for late leaf spot, rust
 and seed quality
 traits has only been possible at low marker density with 225 and 45 SSR loci, respectively.
A wealth of genetic diversity exists in diploid Arachis species
 with the potential to introgress novel alleles into cultivated peanut
[4, 19]. The diploid wild species are not commercially important in the food industry, but they provide a valuable gene pool for sources of resistance to many economically important pathogens and pests
[20–23]. Several linkage maps have been constructed in both A- and B-genome diploid species
[24–28], with the first A-genome linkage map reported being composed of 11 linkage groups with 117 RFLP loci and with a genome coverage of 1,063 cM
. Subsequent A-genome linkage maps consisted of 167 RAPD and 39 RFLPs, 170 SSR, or 369 assorted markers
[25, 27, 28]. An interspecific B-genome linkage map has been constructed with 149 SSR loci covering 1,294 cM
. Synteny between A and B genomes was compared using diploid as well as synthetic amphidiploid linkage maps
[24, 29, 30].
The above-referenced diploid linkage maps were generated from interspecific crosses. Thus, the detection of chromosomal rearrangements within species could be obscured, and comparative A- and B-genome linkage maps could be complicated due to chromosomal rearrangements associated with speciation
[31–35]. Previous comparisons were based on limited numbers of orthologous loci and markers that were primarily derived from genomic sequences. Therefore, map coverage and resolution need to be improved using more function-related and highly transferable markers such as EST-SSRs that facilitate comparative and evolutionary genomics studies. To enrich the currently available SSR marker resources in peanut and gain clearer insights into karyotypic evolution within the genus, a new set of EST-SSR markers was developed and mapped in an intraspecific B-genome mapping population. Comparative mapping to intraspecific A-genome linkage maps revealed a high degree of macrosynteny.
Results and discussion
A total of 101,132 unigenes representing ca. 37 Mb of the A. hypogaea genome (Additional file
1) were mined for SSRs. We found 7,413 perfectly repeated di-, tri-, tetra-, penta-, and hexa-nucleotide motifs (7.3% of the unigenes contained SSRs). The SSR frequency in the above EST resources is comparable with previous reports in cultivated peanut
[36, 37], and wild Arachis species
. The overall SSR density was 3,190 bp per Mb and corresponded to approximately 1 per 5 kb of the genic region, which is similar to a previous report of 1/5.5 kb in cultivated peanut
. higher than Arabidopsis, and barley
 but slightly lower than that reported for rice
 and pepper
, yet within the range of most other plant species (~5%)
. The average SSR length was about 16 bp with almost 88% of SSRs shorter than 22 bp. Among repeat motifs, dinucleotides were predominant (53.3%), which was inconsistent with the study of Koilkonda et al.
, who found trinucleotide repeats to be the most abundant (66.8%). Discrepancies observed in various studies could be explained by the degree of representation of dinucleotide rich UTRs in the genic sequences used or by EST database mining software and SSR search criteria
. Dinucleotide repeat motifs were predominantly distributed in the UTRs while trinucleotide repeat motifs were more frequent in exons. Considering the mode of slippage-mediated mutations, it is unlikely that a large proportion of the dinucleotides would be present in the coding regions; mutations in trinucleotides or their multiples would only be tolerated if they do not disturb the open reading frame.
In our study, the most common dinucleotide repeat motif class was (AG)n (61.0%), while the least common repeat class was (CG)n (1.2%). Similarly, the most abundant trinucleotide motif class was (AAG)n (35.2%), and the least frequent was (CGA)n (0.9%). Motif classes (CG)n and (CGA)n are also relatively infrequent in other plants as well as animals
. The motif types (AG)n and (AAG)n have been reported as the most common di- and tri-nucleotide repeats identified in other plant EST databases
[43, 45–47], including peanut (A. hypogaea)
[15, 30, 37, 40, 48, 49].
We designed 2,138 EST-SSR primer pairs from the identified SSRs (Additional file
2), with 94.6% of them targeted to amplify perfect repeats and over 70% of the primers targeted at trinucleotide motif types. By screening a panel containing eight genotypes, (Additional file
3), 15.3% of the 2,138 primer pairs didn't amplify any interpretable fragments, and 82.7% of them were fully transferable between tetraploid and diploid species. As expected, we observed a relatively higher transferability for EST-SSRs than for genomic SSRs; this is most likely due to greater sequence conservation within expressed regions among related species compared to non-coding regions
[25, 36, 43, 50, 51]. The frequency of polymorphism among the four tetraploid genotypes was 11.2% and less than 10.0% between paired tetraploids. The polymorphism between two A. duranensis accessions was 41.9% while between two A. batizocoi accessions was 21.3% (Additional file
4). The polymorphisms for our intraspecific diploid mapping population parents were comparable to the previous reports of interspecific diploid mapping populations for EST-SSR markers
[24, 25, 30]. We also observed that the polymorphism was higher in the A-genome species than in the B-genome species
Our findings also support the general theory that the degree of polymorphism of the SSR marker increases with the total length of the repeat
[25, 53, 54]. A positive correlation was observed between repeat length and polymorphism rate, but the trend seemed more obvious in tetraploid than in diploid genotypes (Additional file
5), which was supported by the previous observations of Moretzsohn et al.
. For example, SSRs with a repeat length >26 bp showed up to 30% polymorphism between the four tetraploid genotypes, while less than 15% polymorphism was observed for SSRs with repeat length <20 bp. However, when all eight genotypes were combined, no trends were observed between the polymorphism and repeat length. For tetraploid genotypes, AG/CT repeats were more polymorphic than GT/CA, while in diploid genotypes, this effect can only be observed between interspecific genotypes in spite of the fact that the polymorphisms among four diploid genotypes were overall high (around 70%) (Table
1). Similar observations that AG/CT repeats were more polymorphic than GT/CA were previously reported in peanut
[25, 36, 48]. The polymorphisms for dinucleotide repeat motif types were generally higher than for trinucleotide repeat motif types. However, no consistent pattern emerged for ranking of polymorphism rate by motif type.
Effect of SSR repeat motif types on frequency of polymorphism among tetraploid and diploid genotypes
No. of markers
Between Tif-runner and GTC20 (%)
Between NC94022 and SunOleic (%)
Among four tetraploid (%)
Between 30081 and 9484 (%)
Between Grif 15036 and PI 475887 (%)
Among four diploid (%)
Among eight genotypes (%)
In addition to ESTs, a total of 9,517 genome survey sequences (GSSs) representing ca. 5.5 Mb of the Arachis genome were mined for SSRs. Overall, 1,168 perfectly repeated di-, tri-, or tetra-nucleotide motifs were identified from 960 unique sequences. By SSCP screening, 97 SSR primer pairs can produce reliable amplification across a panel of 12 genotypes representing different species (Additional file
Collectively, 2,138 newly developed EST-SSR primer pairs (Additional file
2), 97 genomic SSR markers developed from genome survey sequences and 612 genomic SSR primer pairs in the public domain (Additional file
7) were screened for polymorphisms between the parents of an intraspecific A. batizocoi (BB, 2n = 2x = 20) F2 mapping population. Although A. ipaënsis is the more likely B-genome donor than A. batizocoi to tetraploid peanut species A. hypogaea[52, 55, 56], A. batizocoi retains a high level of similarity to the B subgenome of cultivated peanut
. For example, the F1 plant derived from crossing A. hypogaea by a synthetic amphidiploid (A. batizocoi × A. duranensis) produced bivalents, and a few F4 plants from this cross were even able to produce two-seeded pods
. Furthermore, a diversity study indicated that among all the B-genome species, A. batizocoi showed the second closest relationship to A. hypogaea, after A. ipaënsis
. Although the hypothetical B-genome donor A. ipaënsis was not used for linkage mapping in this report because only a single accession is available in the U.S. germplasm collection, our A. batizocoi intraspecific map should still provide a very close representation of the B-genome donor of tetraploid peanut.
The screening of A. batizocoi accessions 9484 and 30081 produced 455 polymorphic EST-SSR and 171 polymorphic genomic SSR markers. After excluding makers with numerous and/or faint bands and abnormal segregation ratios (markers that showed extreme segregation ratios were assumed to be caused by loci with indistinguishable bands), a total of 481 markers were used for linkage map construction. Of these, 449 loci (including 347 loci from the newly developed EST-SSR markers, 14 loci from genomic SSR markers developed from GSS sequences, and 88 loci from the genomic SSR markers already reported
[18, 25, 48, 49, 59–62]) were mapped into 16 linkage groups (LGs), 14 of which aligned with the 10 chromosome pairs of diploid peanut numbered according to colinearity with A. duranensis (A- genome) linkage groups from Nagy et al.
. The remaining two small linkage groups had no common markers with A- genome linkage groups, thus their chromosomal locations are unknown. Their lengths were 1.1 cM and 9.2 cM, respectively. One LG consisted of markers GM2227 and GM1611, and another LG was composed of markers GM1241 and GM748.
Overall, the linkage map covered 1,278.6 cM, with marker densities ranging from 1.1 cM/locus in LG11 to 9.2 cM/locus in LG12, giving an average density of 2.9 cM/locus for the entire map (Table
1). The LGs ranged from 1.1 to 210 cM in length, and had two (4/9B.1, 4/9B.2, and the two linkage groups described above) to 80 (4/9B) marker loci. Gaps larger than 30.0 cM were observed only at the end of 1B (31.9 cM) and 8B (33.6 cM) (Figure
1). The number of linkage groups observed in this study is larger than the number of haploid chromosomes in the diploid species (n = 10), which may be due to insufficient markers for the chromosome coverage. Furthermore, since the B-genome linkage groups were named by comparing common markers to A-genome linkage groups, the chromosome location of some small linkage groups could not be identified if they lacked anchor markers. The map length was comparable to previously published diploid peanut genetic maps while the density was the highest among all B-genome linkage maps constructed to date
[24–26]. This is by far the most saturated map constructed in a B-genome diploid peanut, and also represents the first intraspecific map of a B-genome species.
Number of loci, map length and density of each linkage group in the 9484x30081 map
No. of marker/LG
Number of distorted loci (α = 0.05)
Significant segregation distortion (α ≤ 0.05) was observed for 18 (4.0%) marker loci, distributed in seven linkage groups representing six B-genome chromosomes (Additional file
8), which was much lower than previously reported segregation distortion (21.5%) in an interspecific B-genome linkage map
. Similarly, high levels of skewed segregation in interspecific crosses compared to intraspecific crosses have been reported in both barley and cowpea
[64, 65]. In the A. batizocoi map, a total of 11 markers skewed towards parental line 9484 and two towards 30081. Linkage group 1B contained the largest number of distorted marker loci (10 marker loci comprising 55.6% of the total distorted markers). All the marker loci on this linkage group skewed towards 9484 and clustered within a 37.1 cM genomic region, which is 37.8% of the total mapping distance for 1B and 26.3% of the marker loci for this chromosome. Linkage groups 3B, 5B, and 8B had no markers showing segregation distortion (Additional file
Synteny between A-genome and B-genome linkage maps
The B-genome linkage groups were aligned to A-genome linkage groups by 167 putative orthologous marker loci, 10 to 25 in each chromosome (Table
3). Compared to previous reports, which identified 51 shared markers between the A- and B-genome maps
, or 53 SSR markers that mapped on both A and B subgenomes in the tetraploid map
, our comparisons are based on a higher density, higher information content map, and should more precisely position chromosomal rearrangement events within defined genomic regions. Syntenic segments were supported by multiple markers in the same linear order with some incongruities reflecting simple translocations and inversions. A high degree of macrosynteny was observed when comparing the nine major linkage groups identified in the B genome to the 10 major linkage groups (1A-10A) in the A genome (Figure
1). Linkage groups 5B and 5B.1 can both be aligned with linkage group 5A, with putative orthologous markers showing colinearity. Similarly, LGs 8B and 8B.1 can both be aligned with linkage group 8A. Since the subgroups couldn’t be merged into one linkage group with a minimum LOD threshold of 3.0 without a >50 cM gap between the adjacent subgroup loci, the two subgroups were displayed separately.
Number of common markers between the corresponding linkage groups from A- and B-genome genetic maps
No. of marker loci
No. of marker loci
*4/9B including 4/9B, 4/9B.1, 4/9B.2, and 4/9B.3.
**5B including 5B and 5B.1.
***8B including 8B and 8B.1.
Completely syntenic chromosomes
For LGs 2A and 2B, all 13 putative orthologous markers were colinear, covering a map distance of 138.3 cM (100%) on B2 and 103.7 cM (100%) on A2. Therefore, 2B and 2A showed a high degree of macrosynteny as previously reported
, although we did not detect any split correspondence relationship with A2 and B10 as in the previous report
. Ten putative orthologous markers between 10A and 10B indicated well defined macrosynteny. The common markers spanned 81.1 cM on 10B and 75.0 cM on 10A, accounting for 75.7% and 68.0% of the total linkage group length, respectively.
Chromosomes with inversions
The 12 putative orthologous markers between 1A and 1B were clustered into two chromosome segments. One segment was colinear, containing 10 putative orthologous markers with 28.2 cM map coverage on 1B and 34.5 cM on 1A, accounting for 28.7% and 37.7% of the total lengths, respectively. The other chromosome fragment had a reversed colinear order, with two putative orthologous markers spanning 58.5 cM on 1B and 40.7 cM on 1A, accounting for 59.6% and 44.5% of the total lengths, respectively. The reversed region defined by two putative orthologous markers on both 1A and 1B indicated a likely inversion between these two chromosomes. This observation was similar to previous reports
For LGs 3B and 3A, 22 common markers were clustered into two chromosomal segments. The first of these was extensively colinear between 3A and 3B, with 12 putative orthologous markers spanning 88.7 cM (75.3%) on 3B and 74.9 cM (51.4%) on 3A. The other fragment had a reversed order with 10 putative orthologous markers spanning 20.9 cM (17.7%) on 3B and 43.2 cM (29.7%) on 3A. The chromosome segment inversion between 3A and 3B has not been previously reported.
For LGs 5A and 5B, there were 25 putative orthologous markers in total, spanning the entire 5B, 4.1 cM (13.6%) of 5B.1, and 108.0 cM (92.8%) of 5A, with a generally colinear order, except for a chromosome segment with nine putative orthologous marker loci spanning 31.4 cM (27.4%) in 5B, and 49.0 cM (42.1%) in the reversed order for 5A. Since the inverted chromosome segments accounted for nearly 30% of the total linkage groups’ lengths and 36.0% (9 out of 25) of the total putative orthologous markers in both A and B genomes, there could be an inverted chromosome segment between 5A and 5B, which also has not been detected from previous reports.
Between chromosomes 6A and 6B, there were 15 putative orthologous markers that spanned 95.6 cM in 6B and 69.1 cM in 6A, accounting for 58.2% and 65.6% of the total lengths on the linkage groups, respectively. Seven putative orthologous markers were mapped in a colinear manner at the top part of both linkage groups. Segments containing eight putative orthologous markers at the bottom part were involved in an inversion event. The inverted chromosome segment on 6B was 40.5 cM long and accounted for 24.7% of the entire 6B. The corresponding segment on 6A was 23.2 cM, accounting for 22.0% of the entire 6A. This inversion was reported previously
, but they also reported the split correspondence relationships between B6 and A10, which was not revealed in our study.
Complex chromosome rearrangements
In addition to the observed simple inversion events, we also found more complex chromosome rearrangements. This intraspecific A. batizocoi genetic map has nine major linkage groups instead of the expected 10. When the A- and B-genome linkage groups were aligned, one of the major B-genome linkage groups was found to correspond to both 4A and 9A and was therefore named as 4/9B. Another smaller linkage group also had putative orthologous markers with 4A and 9A therefore it was named 4/9B.3. The linkage group 4/9B contains 80 markers spanning 209.9 cM, while 4/9B.3 contains 7 markers spanning 47.1 cM. In addition, there were two smaller LGs, 4/9B.1 and 4/9B.2, containing two markers each with genetic distances of 3.4 cM and 8.6 cM, respectively. They were also designated as fragments of 4/9B according to their putative orthologous markers with 4A and/or 9A.
There were several reasons that we did not further separate 4/9B into 4B and 9B. Firstly, these linkage groups remained inseparable even when increasing the LOD threshold to 20. Secondly, the markers were located densely and evenly, with no obvious gaps between two chromosome fragments. The average marker density was 2.7 cM/locus on 4/9B and 7.9 cM on 4/9B.3, which was comparable with the rest of the genome. Thirdly, the effects on map order due to distorted segregation were minimal. There was only one marker that was distorted on each of 4/9B and 4/9B.3. Double crossover events were evaluated alongside 4/9B and 4/9B.3 but no unusual segregating markers or marker blocks were observed. Fourthly, when aligning 4/9B, 4/9B.1, 4/9B.2 and 4/9B.3 with 4A and 9A separately, by using only putative orthologous markers, 4A with 4/9B, 4/9B.1, 4/9B.2 and 4/9B.3 showed complete colinearity, while 9A and 4/9B, 4/9B.1, 4/9B.2 and 4/9B.3 showed a possible inversion. Synteny between 4A and 4B had been reported by Moretzsohn et al.
 and Fonceka et al.
. An inversion between 9A and 9B was also consistent with a previous report
. Lastly, all putative orthologous markers between 4A and 9A were interspersed alongside 4/9B, 4/9B.1, 4/9B.2 and 4/9B.3. From a previous cytological study of the intraspecific variability of A. batizocoi using five accessions
, hybrids between accessions 30081 and 9484 had reduced pollen stainability (88.6%). Moreover, cytogenetic analyses of F1s from the same cross showed a low frequency of quadrivalents
[11, 66], indicating a reciprocal translocation that would cause pairing between two non-homologous chromosomes during meiosis. The similarity between the two chromosomes involved in reciprocal translocation could explain the integrated linkage group 4/9B in our study.
A quadrivalent relationship was observed when aligning 7B and 8B with 7A and 8A. On linkage group 7B, there were two colinearly located putative orthologous markers with 7A, covering 3.3 cM and accounting for 9.0% of the total length on 7B. The bottom segment had seven colinearly located putative orthologous markers with 8A, spanning 21.4 cM, about 58.5% of the entire 7B. On 8B, a colinear fragment with 7A was located from 2.2 cM (from the top) to 91.2 cM (from the top), with nine putative orthologous markers spanning a genetic distance of 89.0 cM (48.1% of 8B), whereas the other 18 putative orthologous markers with 8A were colinearly located from 67.2 cM to 184.0 cM (also from the top of the linkage group), spanning a genetic distance of 116.8 cM (63.1% of 8B). Moreover, three colinear putative orthologous markers with 8A also were found on 8B.1. This may indicate a reciprocal translocation between 7B and 8B. Similar translocations on corresponding linkage groups were found both at the diploid level
 and tetraploid level (from the crosses of a tetraploid variety with a tetraploid AABB amphidiploid)
. However, the previous reports did not identify the correspondence between LGs a07 (corresponding to LG 7A in our study) and b08 (corresponding to LG 8B in our study), which might due to their relatively shorter chromosome coverage on LG 8B (29.8 cM in the synthetic amphidiploid map and 86.4 cM in the diploid map versus 192.2 cM in this report)
[24, 30]. Therefore, the rearrangement is most likely a reciprocal translocation, and may have contributed to the divergence of A and B genomes, perhaps as an ancient event that occurred before peanut polyploidization and remained subsequent to tetraploidization of cultivated peanut
. Although the chromosomes of A. hypogaea have differentiated botanical varieties and individual lines that can be separated based on location of the secondary constriction and symmetry of chromosomes
, introgression of disease and insect resistance traits from wild species into the cultivated peanut has been successful
Comparative mapping summary
Comparison of the genetic linkage maps of A and B genomes indicated that the chromosomal differences between these two species could be explained by a total of five inversions and one reciprocal translocation under our current mapping resolution. Because of the possible minor differences in ordering of tightly linked markers, confident identification of small inversions is more difficult than translocations. To detect the confidence of inversion, we compared our predicted inversions with previously published maps. Moretzsohn et al.
, used different A- and B-genome diploid species and showed four inversions and one translocation by comparing diploid A- and B-genome maps from interspecific mapping populations. Fonceka et al.
 used a synthetic allotetraploid population to compare linkage maps of the A and B subgenomes, which revealed at least three inversions, while Burow et al.
 revealed four inversions. Our study identified additional translocation events by using intraspecific diploid mapping populations, and also detected a reciprocal translocation within the B-genome species A. batizocoi.
Chromosomal rearrangements are common within and among A- and B-genome diploid species. For example, in an analysis to determine the intraspecific variability within the B-genome species A. batizocoi, quadrivalents, hexavalents and octavalents were observed during meiosis in F1 hybrids of different accessions, indicating one to three reciprocal translocations that differentiate these accessions. Thus, karyotypic evolution via translocations was considered to be an important mechanism for species differentiation
. In our study, the linkage map from A. duranensis was used as the reference map from which the chromosome rearrangements between A and B genomes were inferred
. However, chromosomal rearrangements within the A-genome could exist when considering the higher genetic diversity among the various accessions of the A-genome species, A. duranensis. Previous research found a low frequency (0.01-0.26/PMC) of multivalents in 12 of 27 hybrids from crosses of A. duranensis accessions
, and quadrivalents were identified in all these 12 hybrids, which likely represents chromosome translocations within this species. The asymmetrical chromosomes found in different accessions, furthermore, indicated the presence of translocations. In addition, varied fertility of F1s (from less than 4.7% to greater than 95%) and the diverse morphological traits also indicated wide genetic diversity in A. duranensis[12, 69]. Univalents, laggards, and multivalents can all be detected in intraspecific A-genome and interspecific A- by B-genome F1 hybrids, indicating the prevalence of chromosomal rearrangements in peanut diploid species. We herein verify by genetic mapping that a chromosomal translocation has occurred within the B-genome species A. batizocoi. Our findings of one reciprocal translocation between chromosomes 7 and 8 in A-genome species A. duranensis and B-genome species A. batizocoi is consistent with previous reports
[24, 30]. The slight discrepancy of the number of inversions between A- and B-genome chromosomes could be due to the genetic variation of the different accessions/species used in the mapping population development.
The present study developed and characterized an extensive set of EST-SSR and genomic SSR markers. Comparative mapping of our intraspecific A- and B-genome populations showed a high degree of macrosynteny between A- and B-genome diploid species of peanut. Consistent with previous cytological studies, it was evident that chromosomal rearrangements were common within and between both A- and B-genome diploid species. Karyotypic evolution via translocations could be an important mechanism for differentiation of the species. Our findings will facilitate an understanding of tetraploid peanut genome origin and eventually promote its genetic improvement.
An F2 population consisting of 94 plants was developed by selfing four F1 plants from the intraspecific cross of two A. batizocoi lines PI 298639 (accession no. K 9484) and PI 468327 (accession no. GKBSPSc 30081). The A-genome F2 mapping population also consisted of 94 plants from the intraspecific cross of A. duranensis PI 475887 and Grif 15036
Plant materials for EST-SSR characterization
Eight genotypes including four tetraploid and four diploid accessions were used to screen all 2,138 markers. The four tetraploid genotypes in the screening panel included Tifrunner (A. hypogaea subsp. hypogaea var. hypogaea), a runner-type peanut cultivar; GT-C20 (A. hypogaea subsp. fastigiata var. vulgaris), a Spanish exotic accession with reduced aflatoxin contamination; NC94022 (A. hypogaea subsp. hypogaea), an exotic accession with higher resistance to TSWV (tomato spotted-wilt virus) derived from var. hirsuta; and SunOleic 97R (A. hypogaea subsp. hypogaea var. hypogaea), a high oleic runner-type cultivar. For the diploid genotypes, PI 475887 and Grif 15036 are A-genome germplasm accessions of A. duranensis, while accessions 9484 and 30081 are B-genome germplasm accessions in A. batizocoi. They are parental lines of the respective A- and B-genome mapping populations.
SSR discovery, marker development, and length polymorphism screening
Mining the peanut EST database for SSRs
The sequence database used for SSR marker development harbors a total of 70,771 long-read (Sanger) ESTs and 270,957 short-read (454) ESTs assembled into 101,132 unigenes (Accession: PRJNA49471; Additional file
1). Unigenes in the transcript assembly were screened for perfect repeat motifs using SSR-IT (
http://www.gramene.org/db/markers/ssrtool) and for imperfect motifs using FastPCR (
http://primerdigital.com/fastpcr.html). The repeat count (n) threshold for each motif type was set for n ≥ 5. Information on repeat motif, repeat number, and SSR start and end positions within the respective ESTs were extracted from the SSR-IT output. The grouping of SSR motifs into respective repeat classes was performed following the method of Jurka and Pethiyagoda
. Flanking forward and reverse primers were designed using Primer3 (
. The parameters were set as follows: primer length from 19 to 23 with 21 nucleotides as optimum; amplification size of 100 to 400 base pairs; annealing temperatures from 59°C to 63°C with a maximum difference of 3°C; and GC contents from 25% to 45%.
Mining genome survey sequences (GSSs) for SSRs
Methylation filtered (ME) and unfiltered (U) genome libraries were constructed from A. duranensis, A. batizocoi, and A. hypogaea by Orion Genomics (Saint Louis, Missouri)
. A total of 9,517 unique genome survey sequences (GSS) were used for mining SSRs.
Similarly, SSR-IT was used to screen for perfect repeat motifs (
http://www.gramene.org/db/markers/ssrtool), and Primer3 (
 was used for primer design. Primers were screened for overall amplification quality against 12 genotypes, including four diploid and eight tetraploid accessions (Additional file
6), and with SSCP gels by silver staining according to protocols described previously
SSR marker genotyping
Genomic DNA was isolated from young leaves by a modified cetyltrimethylammonium bromide (CTAB) method
. SSR markers were genotyped on an ABI3730XL Capillary DNA Sequencer (Applied Biosystems, Foster City, CA) using forward primers labeled with FAM, HEX, or TAMRA fluorophores. PCR was performed in a 12 μL reaction mixture containing 1.0× PCR buffer, 2.5 mM Mg++, 0.2 mM each of dNTPs, 5.0 pmol of each primer, 0.5 unit of Taq polymerase, and 10 ng of genomic DNA. Touchdown PCR was used to reduce spurious amplification
. The SSR markers were screened for amplification and length polymorphisms using GeneMapper 3.0 software (Applied Biosystems, Foster City, CA).
Genetic mapping, macrosynteny analysis, and cMap database construction
A total of 481 polymorphic markers were used to screen 94 F2 progenies for map construction. Segregation distortion at each marker locus was tested against the expected segregation ratios (1:2:1 for codominant markers and 3:1 for dominant markers) using a chi-square goodness of fit test. Genetic maps were constructed using Mapmaker 3.0, with error detection on
[76, 77]. The initial linkage groups were first determined using the “group” application with a minimum likelihood odds (LOD) threshold of 15 and a maximum recombination fraction (θ) of 0.35. After aligning the draft map with the A-genome map
, the LOD score was relaxed to 5 with θ =0.35 for a second analysis to merge the linkage groups that could align with the homoeologous A-genome linkage groups. Map distances (cM) were calculated using the Kosambi mapping function
. The “try”, “compare”, and “ripple” commands were used to confirm the marker order. Mapchart 2.2 was used for the graphic visualization of the linkage groups
The A. batizocoi (B-genome) linkage groups were numbered based on colinearity to A. duranensis (A-genome) linkage groups except with the suffix “B”. The colinear subgroups in A. batizocoi were named by identical numbers with numerical suffixes, while an un-separated linkage group was identified by using linkage group numbers from the fused groups.
This research was supported by funding from the USDA National Institute of Food and Agriculture National Research Initiative Competitive Grants Program (#2006-35604-17242) awarded to SJK and POA and by the National Peanut Board, the Peanut Foundation, the Georgia Seed Development Commission, and Georgia Research Alliance endowment funding awarded to SJK.
Institute of Plant Breeding, Genetics, and Genomics, The University of Georgia
Department of Horticulture, The University of Georgia
Department of Crop Science, North Carolina State University
Smartt J, Gregory WC: Interspecific cross-compatibility between the cultivated peanut and Arachis hypogaea L. and its behavior in backcrosses.Oleάgineux 1967, 22:455–459.
Kochert G, Stalker H, Gimenes M, Galgaro M, Lopes C, Moore K: RFLP and cytogenetic evidence on the origin and evolution of allotetraploid domesticated peanut, Arachis hypogaea (Leguminosae).Am J Bot 1996,83(10):1282–1291.View Article
Singh AH, Simpson CE: Biosystematics and genetic resources, Chap. 4. In The Groundnut Crop: A Scientific Basis for Improvement. Edited by: Smartt J. London: Chapman & Hall; 1994.
Bertioli DJ, Seijo G, Freitas FO, Valls JFM, Leal-Bertioli SCM, Moretzsohn MC: An overview of peanut and its wild relatives.Plant Genetic Resources 2011,9(1):134–149.View Article
Seijo J, Lavia G, Fernandez A, Krapovickas A, Ducasse D, Moscone E: Physical mapping of the 5S and 18S–25S rRNA genes by FISH as evidence that Arachis duranensis and A. ipaensis are the wild diploid progenitors of A. hypogaea (Leguminosae).Am J Bot 2004, 91:1294–1303.PubMedView Article
Smartt J, Gregory WC, Gregory MP: The genomes of Arachis hypogaea 1. Cytogenetics studies of putative genome donors.Euphytica 1978, 27:665–675.View Article
Robledo G, Seijo G: Species relationships among the wild B genome of Arachis species (section Arachis) based on FISH mapping of rDNA loci and heterochromatin detection: a new proposal for genome arrangement.Theor Appl Genet 2010, 121:1033–1046.PubMedView Article
Seijo G, Lavia GI, Fernandez A, Krapovickas A, Ducasse DA, Bertioli DJ, Moscone EA: Genomic relationships between the cultivated peanut (Arachis hypogaea, Leguminosae) and its close relatives revealed by double GISH.Am J Bot 2007,94(12):1963–1971.PubMedView Article
Halward TM, Stalker HT, Larue EA, Kochert G: Genetic variation detectable with molecular markers among unadapted germplasm resources of cultivated peanut and related wild species.Genome 1991,34(6):1013–1020.View Article
Ramos M, Fleming G, Chu Y, Akiyama Y, Gallo M, Ozias-Akins P: Chromosomal and phylogenetic context for conglutin genes in Arachis based on genomic sequence.Mol Genet Genomics 2006,275(6):578–592.PubMedView Article
Stalker HT, Dhesi JS, Parry DC, Hahn JH: Cytological and interfertility relationships of Arachis section Arachis.Am J Bot 1991,78(2):238–246.View Article
Stalker HT: A new species in section Arachis of peanuts with a D genome.Am J Bot 1991,78(5):630–637.View Article
Pandey MK, Monyo E, Ozias-Akins P, Liang X, Guimarães P, Nigam SN, Upadhyaya HD, Janila P, Zhang X, Guo B, et al.: Advances in Arachis genomics for peanut improvement.Biotechnol Adv 2012,30(3):639–651.PubMedView Article
Temsch EM, Greilhuber J: Genome size variation in Arachis hypogaea and A. monticola re-evaluated.Genome 2000,43(3):449–451.PubMed
Koilkonda P, Sato S, Tabata S, Kenta S, Hirakawa H, Sakai H, Sasamoto S, Watanabe A, Wada T, Kishida Y, et al.: Large-scale development of expressed sequence tag-derived simple sequence repeat markers and diversity analysis in Arachis spp.Mol Breed 2012,30(1):125–138.PubMedView Article
Sujay V, Gowda MVC, Pandey MK, Bhat RS, Khedikar YP, Nadaf HL, Gautami B, Sarvamangala C, Lingaraju S, Radhakrishan T, et al.: Quantitative trait locus analysis and construction of consensus genetic map for foliar disease resistance based on two recombinant inbred line populations in cultivated groundnut (Arachis hypogaea L.).Mol Breed 2012,30(2):773–788.PubMedView Article
Sarvamangala C, Gowda MVC, Varshney RK: Identification of quantitative trait loci for protein content, oil content and oil quality for groundnut (Arachis hypogaea L.).Field Crops Res 2011,122(1):49–59.View Article
Moretzsohn M, Hopkins M, Mitchell S, Kresovich S, Valls J, Ferreira M: Genetic diversity of peanut (Arachis hypogaea L.) and its wild relatives based on the analysis of hypervariable regions of the genome.BMC Plant Biol 2004, 4:11–11.View Article
Simpson CE: Use of wild Arachis species/introgression of genes into A. hypogaea L.Peanut Sci 2001,28(2):114–116.View Article
Abdou Y, Gregory W, Cooper W: Sources and nature of resistance to Cercospora arachidicola Hori and Cercosporodium personatum (Beck & Curtis) Deighton in Arachis species.Peanut Sci 1974, 1:6–11.View Article
Herbert TT, Stalker HT: Resistance to peanut stunt virus in cultivated and wild Arachis species.Peanut Sci 1981, 8:45–47.View Article
Pande S, Narayana Rao J: Resistance of wild Arachis species to late leaf spot and rust in greenhouse trials.Plant Dis 2001, 85:851–855.View Article
Subrahmanyam P, Moss J, Rao V: Resistance to peanut rust in wild Arachis species.Plant Dis 1983, 67:209–212.View Article
Moretzsohn M, Barbosa A, Alves-Freitas D, Teixeira C, Leal-Bertioli S, Guimaraes P, Pereira R, Lopes C, Cavallari M, Valls J, et al.: A linkage map for the B-genome of Arachis (Fabaceae) and its synteny to the A-genome.BMC Plant Biol 2009,9(1):40.PubMedView Article
Moretzsohn M, Leoi L, Proite K, Guimaraes P, Leal-Bertioli S, Gimenes M, Martins W, Valls J, Grattapaglia D, Bertioli D: A microsatellite-based, gene-rich linkage map for the AA genome of Arachis (Fabaceae).Theor Appl Genet 2005, 111:1060–1071.PubMedView Article
Halward T, Stalker HT, Kochert G: Development of an RFLP linkage map in diploid peanut species.Theor Appl Genet 1993,87(3):379–384.View Article
Garcia GM, Stalker HT, Schroeder E, Lyerly JH, Kochert G: A RAPD-based linkage map of peanut based on a backcross population between the two diploid species Arachis stenosperma and A. cardenasii.Peanut Sci 2005,32(1):1–8.View Article
Leal-Bertioli S, Jose AC, Alves-Freitas D, Moretzsohn M, Guimaraes P, Nielen S, Vidigal B, Pereira R, Pike J, Favero A, et al.: Identification of candidate genome regions controlling disease resistance in Arachis.BMC Plant Biol 2009,9(1):112.PubMedView Article
Burow M, Simpson C, Starr J, Paterson A: Transmission genetics of chromatin from a synthetic amphidiploid to cultivated peanut (Arachis hypogaea L.): Broadening the gene pool of a monophyletic polyploid species.Genetics 2001, 2:823–837.
Fonceka D, Hodo-Abalo T, Rivallan R, Faye I, Sall M, Ndoye O, Favero A, Bertioli D, Glaszmann JC, Courtois B, et al.: Genetic mapping of wild introgressions into cultivated peanut: a way toward enlarging the genetic basis of a recent allotetraploid.BMC Plant Biol 2009,9(1):103.PubMedView Article
Lai Z, Nakazato T, Salmaso M, Burke JM, Tang S, Knapp SJ, Rieseberg LH: Extensive chromosomal repatterning and the evolution of sterility barriers in hybrid sunflower species.Genetics 2005,171(1):291–303.PubMedView Article
Strasburg JL, Scotti-Saintagne C, Scotti I, Lai Z, Rieseberg LH: Genomic patterns of adaptive divergence between chromosomally differentiated sunflower species.Mol Biol Evol 2009,26(6):1341–1355.PubMedView Article
Rieseberg LH: Chromosomal rearrangements and speciation.Trends Ecol Evol 2001,16(7):351–358.PubMedView Article
Noor MAF, Grams KL, Bertucci LA, Reiland J: Chromosomal inversions and the reproductive isolation of species.Proc Nat Acad Sci USA 2001,98(21):12084–12088.PubMedView Article
Liang X, Chen X, Hong Y, Liu H, Zhou G, Li S, Guo B: Utility of EST-derived SSR in cultivated peanut (Arachis hypogaea L.) and Arachis wild species.BMC Plant Biol 2009,9(1):35.PubMedView Article
Song GQ, Li MJ, Xiao H, Wang XJ, Tang RH, Xia H, Zhao CZ, Bi YP: EST sequencing and SSR marker development from cultivated peanut (Arachis hypogaea L.).Electron J Biotechnol 2010,13(3):1–9.View Article
Proite K, Leal-Bertioli S, Bertioli D, Moretzsohn M, da Silva F, Martins N, Guimaraes P: ESTs from a wild Arachis species for gene discovery and marker development.BMC Plant Biol 2007,7(1):7.PubMedView Article
Cardle L, Ramsay L, Milbourne D, Macaulay M, Marshall D, Waugh R: Computational and experimental characterization of physically clustered simple sequence repeats in plants.Genetics 2000,156(2):847–854.PubMed
Thiel T, Michalek W, Varshney R, Graner A: Exploiting EST databases for the development and characterization of gene-derived SSR-markers in barley (Hordeum vulgare L.).Theor Appl Genet 2003,106(3):411–422.PubMed
Varshney R, Thiel T, Stein N, Langridge P, Graner A: In silico analysis on frequency and distribution of microsatellites in ESTs of some cereal species.Cell Mol Biol Lett 2002,7(2A):537–546.PubMed
Portis E, Nagy I, Sasvári Z, Stágel A, Barchi L, Lanteri S: The design of Capsicum spp. SSR assays via analysis of in silico DNA sequence, and their potential utility for genetic mapping.Plant Sci 2007,172(3):640–648.View Article
Varshney RK, Graner A, Sorrells ME: Genic microsatellite markers in plants: Features and applications.Trends Biotechnol 2005,23(1):48–55.PubMedView Article
Tóth G, Gáspári Z, Jurka J: Microsatellites in different eukaryotic genomes: survey and analysis.Genome Research 2000,10(7):967–981.PubMedView Article
Kantety RV, La Rota M, Matthews DE, Sorrells ME: Data mining for simple sequence repeats in expressed sequence tags from barley, maize, rice, sorghum and wheat.Plant Mol Biol 2002,48(5):501–510.PubMedView Article
Morgante M, Hanafey M, Powell W: Microsatellites are preferentially associated with nonrepetitive DNA in plant genomes.Nat Genet 2002,30(2):194–200.PubMedView Article
Temnykh S, DeClerck G, Lukashova A, Lipovich L, Cartinhour S, McCouch S: Computational and experimental analysis of microsatellites in rice (Oryza sativa L.): Frequency, length variation, transposon associations, and genetic marker potential.Genome Res 2001,11(8):1441–1452.PubMedView Article
Ferguson M, Burow M, Schulze S, Bramel P, Paterson A, Kresovich S, Mitchell S: Microsatellite identification and characterization in peanut (A. hypogaea L.).Theor Appl Genet 2004,108(6):1064–1070.PubMedView Article
Hopkins M, Casa A, Wang T, Mitchell S, Dean R, Kochert G, Kresovich S: Discovery and characterization of polymorphic simple sequence repeats (SSRs) in peanut.Crop Sci 1999, 39:1243–1247.View Article
Luo M, Dang P, Guo BZ, He G, Holbrook CC, Bausher MG, Lee RD: Generation of expressed sequence tags (ESTs) for gene discovery and marker development in cultivated peanut.Crop Sci 2005,45(1):346–353.View Article
Cuc L, Mace E, Crouch J, Quang V, Long T, Varshney R: Isolation and characterization of novel microsatellite markers and their application for diversity assessment in cultivated groundnut (Arachis hypogaea).BMC Plant Biol 2008,8(1):55–55.PubMedView Article
Milla S, Isleib T, Stalker H: Taxonomic relationships amongArachissect. Arachis species as revealed by AFLP markers.Genome 2005,48(1):1–11.PubMedView Article
Gupta PK, Varshney RK: The development and use of microsatellite markers for genetic analysis and plant breeding with emphasis on bread wheat.Euphytica 2000,113(3):163–185.View Article
Weber JL: Informativeness of human (dC-dA)n · (dG-dT)n polymorphisms.Genomics 1990,7(4):524–530.PubMedView Article
Koppolu R, Upadhyaya H, Dwivedi S, Hoisington D, Varshney R: Genetic relationships among seven sections of genus Arachis studied by using SSR markers.BMC Plant Biol 2010,10(1):15.PubMedView Article
Krapovickas A, Gregory W: Taxonomia del genero Arachis (Leguminosae).Bonplandia 1994, 8:1–186.
Simpson CE: Pathways for introgression of pest resistance into Arachis hypogaea L.Peanut Sci 1991,18(1):22–26.View Article
Singh AK: Utilization of wild relatives in the genetic-improvement of Arachis hypogaea L. 8. Synthetic amphidiploids and their importance in interspecific breeding.Theor Appl Genet 1986,72(4):433–439.View Article
Krishna G, Zhang J, Burow M, Pittman R, Delikostadinov S, Lu Y, Puppala N: Genetic diversity analysis in valencia peanut (Arachis hypogaea L.) using microsatellite markers.Cell Mol Biol Lett 2004, 9:685–697.PubMed
Palmieri D, Bechara M, Curi R, Gimenes M, Lopes C: Novel polymorphic microsatellite markers in section Caulorrhizae (Arachis, Fabaceae).Mol Ecol Notes 2005, 5:77–79.View Article
Palmieri D, Hoshino A, Bravo J, Lopes C, Gimenes M: Isolation and characterization of microsatellite loci from the forage species Arachis pintoi (genus Arachis).Mol Ecol Notes 2002, 2:551–553.View Article
He G, Meng R, Newman M, Gao G, Pittman R, Prakash C: Microsatellites as DNA markers in cultivated peanut (Arachis hypogaea L.).BMC Plant Biol 2003, 3:3–3.PubMedView Article
Nagy ED, Guo Y, Tang S, Bowers JE, Okashah RA, Taylor CA, Zhang D, Khanal S, Heesacker AF, Khalilian N, et al.: A high-density genetic map of Arachis duranensis, a diploid ancestor of cultivated peanut.BMC Genomics 2012,13(1):469.PubMedView Article
Becker J, Vos P, Kuiper M, Salamini F, Heun M: Combined mapping of AFLP and RFLP markers in barley.Mol Genet Genomics 1995,249(1):65–73.View Article
Menéndez CM, Hall AE, Gepts P: A genetic linkage map of cowpea (Vigna unguiculata) developed from a cross between two inbred, domesticated lines.Theor Appl Gen 1997,95(8):1210–1217.View Article
Stalker HT, Dhesi JS, Parry DC: An analysis of the B genome species Arachis batizocoi (Fabaceae).Plant Syst Evol 1991, 174:159–169.View Article
Stalker HT, Dalmacio RD: Karyotype analysis and relationships among varieties of Arachis hypogaea L.Cytologia 1986, 58:617–629.View Article
Stalker HT, Dhesi JS, Kochert D: Genetic diversity within the species Arachis duranensis Krapov. & W.C. Gregory, a possible progenitor of cultivated peanut.Genome 1995,38(6):1201–1212.PubMedView Article
Jurka J, Pethiyagoda C: Simple repetitive DNA sequences from primates: Compilation and analysis.J Mol Evol 1995,40(2):120–126.PubMedView Article
Rozen S, Skaletsky H: Primer3 on the WWW for general users and for biologist programmers.Meth Mol Biol 2000, 132:365–386.
Whitelaw CA, Barbazuk WB, Pertea G, Chan AP, Cheung F, Lee Y, Zheng L, van Heeringen S, Karamycheva S, Bennetzen JL, et al.: Enrichment of gene-coding sequences in maize by genome filtration.Science 2003,302(5653):2118–2120.PubMedView Article
Orita M, Iwahana H, Kanazawa H, Hayashi K, Sekiya T: Detection of polymorphisms of human DNA by gel-electrophoresis as single-strand conformation polymorphisms.Proc Nat Acad Sci USA 1989,86(8):2766–2770.PubMedView Article
Murray MG, Thmposon WF: Rapid isolation of high molecular weight plant DNA.Nucleic Acids Res 1980,8(19):4321–4325.PubMedView Article
Don R, Cox P, Wainwright B, Baker K, Mattick J: ‘Touchdown’ PCR to circumvent spurious priming during gene amplification.Nucleic Acids Res 1991, 19:4008.PubMedView Article
Lander ES, Green P, Abrahamson J, Barlow A, Daly MJ, Lincoln SE, Newburg L: MAPMAKER: An interactive computer package for constructing primary genetic linkage maps of experimental and natural populations.Genomics 1987, 1:174–181.PubMedView Article
Lincoln SE, Daly MJ, Lander ES: Constructing genetic maps with MAPMAKER/EXP 3.0 Manual. Cambridge: Whitehead Institute; 1992.
Kosambi DD: The estimation of map distance from recombination values.Ann Eugen 1944, 12:172–175.
Voorrips RE: MapChart: Software for the graphical presentation of linkage maps and QTLs.J Hered 2002, 93:77–78.PubMedView Article
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (
http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.