Skip to main content
  • Research article
  • Open access
  • Published:

Genome-wide identification of microsatellite markers from cultivated peanut (Arachis hypogaea L.)

Abstract

Background

Microsatellites, or simple sequence repeats (SSRs), represent important DNA variations that are widely distributed across the entire plant genome and can be used to develop SSR markers, which can then be used to conduct genetic analyses and molecular breeding. Cultivated peanut (A. hypogaea L.), an important oil crop worldwide, is an allotetraploid (AABB, 2n = 4× = 40) plant species. Because of its complex genome, genomic marker development has been very challenging. However, sequencing of cultivated peanut genome allowed us to develop genomic markers and construct a high-density physical map.

Results

A total of 8,329,496 SSRs were identified, including 3,772,653, 4,414,961, and 141,882 SSRs that were distributed in subgenome A, B, and nine scaffolds, respectively. Based on the flanking sequences of the identified SSRs, a total of 973,984 newly developed SSR markers were developed in subgenome A (462,267), B (489,394), and nine scaffolds (22,323), with an average density of 392.45 markers per Mb. In silico PCR evaluation showed that an average of 88.32% of the SSR markers generated only one in silico-specific product in two tetraploid A. hypogaea varieties, Tifrunner and Shitouqi. A total of 39,599 common SSR markers were identified among the two A. hypogaea varieties and two progenitors, A. duranensis and A. ipaensis. Additionally, an amplification effectiveness of 44.15% was observed by real PCR validation. Moreover, a total of 1276 public SSR loci were integrated with the newly developed SSR markers. Finally, a previously known leaf spot quantitative trait locus (QTL), qLLS_T13_A05_7, was determined to be in a 1.448-Mb region on chromosome A05. In this region, a total of 819 newly developed SSR markers were located and 108 candidate genes were detected.

Conclusions

The availability of these newly developed and public SSR markers both provide a large number of molecular markers that could potentially be used to enhance the process of trait genetic analyses and improve molecular breeding strategies for cultivated peanut.

Background

Cultivated peanut or groundnut (Arachis hypogaea L.) is a globally important legume that is widely planted in Asia, Africa, America, and other areas because it is rich in seed oil and protein. Thus, peanut has great significance for fighting malnutrition and ensuring food security. In China, peanut accounts for almost half of the total output of all oil crops and is increasingly important as an oil and protein crop. Therefore, it is critical to improve peanut production and quality to ensure an edible oil supply.

However, peanut production is often constrained by factors such as drought, salinity, and disease [1]. During the past 10 years, with the development of peanut genomics, there have been several successful achievements in peanut trait mapping [2,3,4,5] and molecular breeding [6, 7]. Owing to the limited available genetic markers and low-density genetic maps, those studies could not provide optimal resolution of trait dissection and identify candidate genes. Therefore, the development of a high-density genetic map is particularly urgent for peanut trait mapping and breeding.

Simple sequence repeats (SSRs) are genomic fragments that consist of tandemly repeated units that are present in both coding and non-coding regions of the genome [8, 9]. SSR markers, designed by flanking sequences, are useful for and widely applied in plant genetic analyses and marker-assisted selection breeding. SSRs derived from expressed sequence tags (ESTs), transcriptome sequences, and genomic DNA sequences are referred to as EST-SSRs, transcriptome-SSRs, and g-SSRs, respectively. In the past decade, several hundred EST-SSR markers were developed by investigating ESTs [10, 11], and thousands of transcriptome-SSR markers were identified based on different transcriptome libraries of cultivated peanut [12,13,14]. Meanwhile, two integrated consensus genetic maps with thousands of different types of markers, such as EST-SSRs, transcriptome-SSRs, and g-SSRs, were constructed [15, 16]. Limited by the polymorphisms, these available SSRs were insufficient for constructing high-density genetic maps and enhancing molecular breeding [17].

Cultivated peanut is an allotetraploid (AABB, 2n = 4× = 40) that probably derived from hybridization between two diploids, A. duranensis and A. ipaensis [18, 19]. Recently, breakthroughs have been made in peanut genome sequencing. From 2016 to 2018, the genomes of the two diploid progenitors and one allotetraploid wild species, A.monticola, were successfully sequenced [20,21,22,23]. Based on the genomic sequences of the two diploid progenitors, genome-wide g-SSRs were identified and developed, and a high-density SSR physical map of wild peanut species was constructed [24]. Importantly, in 2019, cultivated peanut genomics research experienced a substantial milestone when genome sequencing was completed for three cultivated peanut : Fuhuasheng [25], Shitouqi [26], and Tifrunner [27]. The high-quality genome assemblies provide the opportunity for developing genome-wide g-SSR markers in cultivated peanut.

Here, we identified genome-wide g-SSRs and developed g-SSR markers from the genome assembly of A. hypogaea cv. Fuhuasheng, a landrace from North China that was sequenced in our previous work [25]. The aims of this study were to: (1) identify genome-wide SSRs and show the distribution of motif length, type, and repeat number between the two subgenomes (A and B); (2) develop g-SSR markers and construct a high-density SSR physical map of cultivated peanut; and (3) evaluate the application of these SSR markers and validate the polymorphisms in different cultivated peanut species. These novel, newly developed g-SSR markers could be helpful for advancing agronomic trait mapping, gene cloning, and molecular breeding of cultivated peanut in the future.

Results and discussion

Whole genome identification of SSRs

In this study, the whole genome sequence of A. hypogaea cv. Fuhuasheng was used to identify SSRs with different repeat motifs, from mono- to hexa-nucleotide. A total of 8,329,496 SSRs were obtained, with a density of ~ 3264.31 SSRs per Mb (Table 1; Additional file 1: Table S1). Penta-nucleotide was the most common type, accounting for more than half of all identified SSRs (57.75%), followed by hexa-nucleotide (26.61%) (Table 1). In addition, more SSRs were identified in subgenome B (4,414,961) than in subgenome A (3,772,653), and 141,882 SSRs were identified on the nine scaffolds of the peanut reference genome assembly (Table 1).

Table 1 Different types of SSRs identified in A. hypogaea L

Analysis of SSR distribution on each chromosome revealed that the largest number of SSRs was present on chromosome A07 (538,928), followed by chromosome B09 (532,326) (Additional file 1: Table S1; Additional file 2: Figure S1). Pearson correlation analysis revealed that chromosome length was significantly positively associated with the number of SSRs of each chromosome (r = 0.996, p < 0.01) (Additional file 1: Table S1). Furthermore, the average density of SSRs was 3283.89 SSRs per Mb, ranging from 3164.75 SSRs per Mb on chromosome A10 to 3843.36 SSRs per Mb on chromosome A08 (Additional file 1: Table S1; Additional file 2: Figure S1). Moreover, the densities of SSRs on chromosomes A08 and B03 were the highest in subgenomes A and B, respectively (Additional file 2: Figure S1).

A total of 501 types of SSR motifs were identified in the peanut genome (Additional file 1: Table S2). The repeat number of all of these motifs ranged from 2 to 335, most of which were concentrated in the top of 50s (Additional file 2: Figure S2a). The proportions of penta- and hexa-nucleotides with two repeats were much greater than those of other types in both subgenome A and B, and ranged from ~ 20 to 50% (Additional file 2: Figure S2b). Of the SSR motif types, the penta-nucleotide type, AAAAT/ATTTT, had the highest occurrence, and accounted for 9.55% of all types, followed by AAATT/AATTT and AAAAG/CTTTT, which accounted for 4.97 and 4.85%, respectively (Additional file 2: Figure S3). For each type of SSR motif, from mono- to hexa-nucleotide, the richest motifs were T (1.48%), TA (1.25%), AAT (0.37%), AAAT (0.34%), AAAAT (1.44%), and ACGCGT (0.64%).

Genome-wide SSR marker development

The flanking sequences of all identified SSRs were used to design suitable forward and reverse primer pairs. A total of 973,984 SSR markers were successfully developed on 20 chromosomes and nine scaffolds in peanut (Additional file 1: Table S1), which accounted for 11.69% of all identified SSRs. Of these newly developed SSR markers, a total of 462,267 and 489,394 SSR markers were located in subgenomes A (Additional file 3: Table S3 and Additional file 4: Table S4) and B (Additional file 5: Table S5 and Additional file 6: Table S6), which accounted for 47.46 and 50.25% of all SSR markers, respectively (Additional file 1: Table S1). In addition, 22,323 SSR markers were mined on the nine scaffolds, and only accounted for 2.29% of all SSR markers (Additional file 7: Table S7).

Based on the start positions of SSR markers, we successfully anchored these markers to the reference genome physical map (Fig. 1). For the physical map, the marker density significantly differed among chromosomes, and ranged from 338.38 per Mb on chromosome B01 to 699.92 per Mb on chromosome A08, with an average density of 392.45 per Mb (Additional file 1: Table S1; Fig. 2). Moreover, the distribution of each chromosome showed that there was lower marker density in the middle of each chromosome except for chromosome B07 (Figs. 1 and 2). In addition, for all SSR markers in subgenome A, subgenome B, and the scaffolds, penta-nucleotide SSRs were the most abundant (491,927) (Additional file 2: Figure S4a) and accounted for more than half of all SSRs (~ 50.51%) (Additional file 2: Figure S4b). Hexa-nucleotide SSRs were the second most abundant (199,408),and represented 20.47% of the SSRs in the two subgenomes and scaffolds, followed by compound and di-nucleotide SSRs. Of the repeat motif types, TA/AT was most abundant, accounting for more than 4.0% of all repeat motifs (Additional file 2: Figure S5). The second highest number of repeat motifs was AAAAT/TTTTA, which represented 3.4% of all repeat motifs. These analyses showed that the AT repeat patterns were the dominant repeat motifs of SSR markers, whereas GC repeat patterns were rare.

Fig. 1
figure 1

Overview of the high-density SSR physical map in peanut (A. hypogaea L.).The bar represents the number of SSR markers within a 1-Mb window

Fig. 2
figure 2

Genome-wide distribution of SSR markers on different chromosomes of the A. hypogaea genome. From the outer edge inward, circles represent the (a) gene position; (b) distribution of all genes; (c) gene density per Mb; (d) distribution of the 1276 integrated public SSR markers; (e) distribution of the 188 SSRs tested by PCR amplification; (f) distribution of all newly developed SSR markers; and (g) SSR marker density per Mb. Blue and red lines represent homologous loci in different and the same chromosomes, respectively

Development of new markers can help improve genetic analysis, gene/QTL mapping, and molecular breeding of crops. Currently, g-SSR markers are common and popular for such analyses, and they have wide applications in molecular genetics and breeding, because they have multiple advantages, such as simplicity, abundance, ubiquity, variation, co-dominance, and multi-allelism [28]. Recently, several studies were devoted to developing different types of SSR markers in peanut, such as EST-SSRs [11, 29], transcriptome-SSRs [12,13,14], and g-SSRs [30], even though the peanut genome had not yet been resolved. With the recent completion of genome sequencing of peanut and two diploid progenitor species, A. duranensis and A. ipaensis, a large number of genome-wide g-SSRs were identified [20,21,22]. Furthermore, tens of thousands of g-SSR markers (51,354 for A. duranensis and 60,893 for A. ipaensis) were also developed from the two progenitor species in 2017 [24]. However, there was limited reports on the development of a large number of genome-wide g-SSR markers from allotetraploid cultivated peanut because of the challenging in its genome sequencing. Fortunately, the genomes of allotetraploid A. hypogaea cv. Fuhuasheng, Shitouqi, and Tifrunner were successfully sequenced [25,26,27]. Here, large-scale genome-wide g-SSR markers were developed from the A. hypogaea cv. Fuhuasheng genome to help enhance genetic and genomic analyses and molecular breeding of peanut.

In silico evaluation of the newly developed SSR markers

To evaluate the amplification specificity of the newly developed SSR markers, the forward and reverse primers of each SSR marker were used for in silico analysis based on the cultivated peanut genome sequences of A. hypogaea cv. Tifrunner and Shitouqi, and genome sequences of its two progenitors, A. duranensis and A. ipaensis (Table 2).

Table 2 In silico PCR products in A. hypogaea and its two progenitors

The number of in silico products indicated that 87.35% (855,422) and 89.29% (874,384) of the SSR markers generated only one in silico-specific product in the twotetraploids A. hypogaea cv. Tifrunner and Shitouqi, respectively. Approximately 4% of markers were mismatched in the two cultivated varieties, and less than 8% of markers generated two products. In addition, less than 0.8% of SSR markers generated more than two products; in particular, only 0.08 and 0.11% of the markers generated ≥6 in silico products. For the two progenitors, a total of 356,778 and 410,764 SSR markers, which accounted for 36.43 and 41.59% of all SSR markers, respectively, generated only one in silico-specific product. However, more than half of total markers were mismatched in the two progenitors. Approximately 3% of all markers generated two in silico products, and less than 0.5% of the markers generated more than two in silico products (Table 2). In total, 1,729,806 and 767,542 markers could generate only one in silico product in cultivated peanut and its two progenitors, respectively. These SSR markers that only generated one in silico product would be potentially useful for molecular breeding in the future.

After clumped, a total of 39,599 SSR markers were shared among the four Arachis species, which indicates that these common markers were very conservative in different Arachis (Fig. 3; Additional file 7: Table S8). Moreover, 9334 and 210,058 markers overlapped between the two progenitors and two cultivated varieties, respectively. In particular, numerous specific SSR markers were also obtained for each genome, which indicates genomic polymorphism among the four different Arachis species.

Fig. 3
figure 3

Venn diagram of SSR markers in different Arachis species

PCR validation of the newly developed SSR markers in different species

To validate amplification of the newly developed SSR markers, a total of 188 SSR markers with motifs ≥5 repeats were arbitrarily and evenly selected for PCR amplification in two cultivars, Fuhuasheng and Yueyouhei4hao (Additional file 7: Table S9). In total, 83 of the SSR markers (44.15%) could amplify clear target products in at least one of the two varieties (Additional file 2: Figure S6). Moreover, 25 of the 83 SSR markers displayed polymorphism between the two varieties. These results indicated that the newly developed SSR markers were available and could be useful for molecular breeding strategies, such as true hybrid F1 offspring selection, in the future.

Integration of publicly available SSR markers

Previously, an integrated consensus genetic map, which contained 5874 markers, was successfully constructed with 20 linkage groups [16]. Based on forward and reverse primer sequences, a total of 5125 markers were mapped to the reference genome, A. hypogaea cv. Fuhuasheng. Finally, 3304 markers with 6838 loci were mined from 20 chromosomes (6699 loci) and nine scaffolds (139 loci) (Additional file 2: Figure S7). In addition, most of the markers (2927), which accounted for 88.6% of all markers, only matched one or two loci. For example, approximately 1832 markers only matched a single locus (Additional file 7: Table S10), and a total of 1095 markers were double loci identification on the reference genome (Additional file 2: Figure S7). The remaining markers had more than three loci, and one marker even matched 258 loci (Additional file 2: Figure S7). Furthermore, most of the multiple loci of these markers were located in the corresponding subgenomic chromosomes, such as A01 vs. B01 and A02 vs. B02 (Fig. 2). This may have been caused by co-evolution of the two peanut subgenomes.

Of these 1832 unique loci that matched markers, a total of 1276 SSR markers were located in the corresponding genetic linkage groups (Additional file 7: Table S10). Based on the physical positions of the 1276 publicly available SSR markers, they were successfully merged with the physical map and anchored to the draft genome sequence of tetraploid cultivated peanut (Fig. 2). This map, which has the highest marker density and the most uniform genetic background, will be of great benefit to molecular breeding, and gene and QTL mining of peanut in the future.

A high-density physical map with uniform genomic positions and coverage is necessary for conducting high-resolution gene/QTL mapping in crops. During the past years, several available QTLs for yield, disease resistance, and other traits were detected in peanut (https://www.peanutbase.org/). However, most of the QTLs were detected across different genetic backgrounds and environments, and only a few QTLs have been used for molecular breeding; therefore, the use of the QTLs in molecular breeding has been limited [31]. Thus, consensus genetic map construction was needed to improve the use of QTLs in molecular breeding. Therefore, two integrated consensus genetic maps with thousands of markers were constructed in 2013 and 2018 [15, 16], and multiple independent consensus QTLs were identified [16]. Here, to obtain a comprehensive consensus map, we integrated the published public SSR markers with the newly developed SSR physical map (Fig. 1). The availability of a high-density physical map would provide an opportunity to generate high-throughput genotyping data for different types of populations, and accelerate mapping and breeding applications of different traits.

Application of SSR markers in peanut

High-density SSR physical maps can be used to screen more markers for QTL fine mapping. For example, a previous study reported that multiple novel QTLs for resistance to leaf spots and tomato spotted wilt virus were identified in peanut based on an improved genetic linkage map with a total of 418 markers, and multiple resistance QTL clusters were detected on linkage A05, which only contained 24 markers [32]. Based on the physical position of the marginal linkage markers (PM179–TC40D04), a total of 6804 new SSR markers were obtained in this region (Additional file 2: Figure S8a). On the A05 linkage, a total of 11 leaf spot QTLs were detected, especially for the QTL qLLS_T13_A05_7, which had the largest LOD value (5.26) and phenotypic variation explanation (15.55%) (Additional file 2: Figure S8b). Moreover, this QTL overlapped with our previous meta-QTL MQTL_LLS_A05.1, sharing the common marker GNB464 (Additional file 2: Figure S8b and c) [16]. Based on the sequences of the flanking markers near this QTL, GNB464 and TC40D04 were successfully anchored in a 1.448-Mb region on chromosome A05 in the A. hypogaea cv. Fuhuasheng reference genome (Additional file 2: Figure S8d). Searching this region in the high-density SSR physical map revealed a total of 819 newly developed markers (Additional file 2: Figure S8d). Gene identification of this region indicated that there were 108 genes, and functional annotation showed that multiple genes were related to disease resistance, such as the TMV resistance protein and NAC domain (Additional file 7: Table S11).

Genetic linkage mapping is an effective strategy for QTL identification in crops. However, a high-density genetic map is a prerequisite for QTL fine mapping of multiple phenotypic traits. In peanut, most QTLs were only preliminarily mapped because of the lack of sufficiently available genetic markers in the last decade [33,34,35,36,37]. Therefore, the development of molecular markers has become increasingly important. Then, high-density integrated consensus maps with thousands of markers were constructed in 2013 and 2018 [15, 16]. However, in cultivated peanut, development of high-density physical markers was rare. Here, we developed millions of SSR markers in cultivated peanut, and this large number of newly developed SSR markers, such as multiple QTLs on chromosome A05 or even on the region of qLLS_T13_A05_7, could provide more possibilities for further fine mapping or even gene cloning.

Conclusions

Development of new genome-wide markers and construction of a physical map with uniformly distributed genome-wide physical markers would aid in the elucidation of complex traits and improvement of molecular breeding. In this study, we identified a total of 8,329,496 genomic SSRs and developed 973,984 newly genomic SSR markers from the cultivated peanut reference genome “Fuhuasheng” using MISA software with default parameters. Moreover, we mined a number of 39,599 common SSR markers from two A. hypogaea cv. Tifrunner and Shitouqi and two progenitors, A. duranensis and A. ipaensis. In addition, we integrated 1276 public SSR loci with the newly developed physical map, and obtained the first high-density genomic physical map for peanut. Finally, we fine mapped a leaf spot quantitative trait locus to a 1.448-Mb region by marker encryption and identified 108 candidate genes. In summary, these newly developed and integrated public SSR markers are an important genomic resource for both accelerating genetic analyses of complex traits and molecular breeding applications in peanut.

Methods

Plant materials and DNA isolation

A cultivated peanut that was de novo sequenced in our previous study [25], A. hypogaea cv. Fuhuasheng, was used as the reference genome for SSR identification. Two peanut varieties, Fuhuasheng and Yueyouhei4hao, were used for amplification and validation of the newly developed SSR markers by polymerase chain reaction (PCR). All of the materials were planted in a field in the summer season (July to November) at the experimental station of Guangdong Academy of Agricultural Sciences, Guangzhou, China. High-quality genomic DNA was extracted from young leaves using a Plant Genomic DNA Extraction Kit (BioTeke Corporation, Beijing, China) according to the manufacturer’s handbook (http://www.bioteke.com/en/?c=show&m=view&id=99).

Identification of SSRs and design for SSR markers

The reference genome sequence of A. hypogaea cv. Fuhuasheng was downloaded from GenBank under accession number SDMP00000000 [25]. Perl scripts from MISA were used to perform SSR identification with the default parameters (http://pgrc.ipk-gatersleben.de/misa/). The identification criteria were as follows: mono-nucleotide repeat motifs with at least 12 repeats, di-nucleotide repeat motifs with five repeats, tri-nucleotide repeat motifs with four repeats, tetra-nucleotide repeat motifs with three repeats, and penta- and hexa-nucleotide repeat motifs with two repeats. Compound SSRs were defined as those with a < 100-nt interval between two repeat motifs [24, 38].

The forward and reverse primers of each SSR were designed based on the flanking sequences of SSR repeat motifs using Primer 3 [39]. Two Perl scripts, p3_in.pl and p3_out.pl, were used for the programmer-to-programmer data interchange between MISA and Primer 3.0 (http://pgrc.ipk-gatersleben.de/misa/primer3.html). The primer design criteria were as follows: primer length was between 18 and 27 bp, melting temperature was 57 to 63 °C, GC content ranged from 30 to 70%, and product size was from 100 to 300 bp.

An R package (http://www.R-project.org), Cmplot, was used to draw the high-density physical map of the newly developed SSR markers, and another R package, Rcircos, was used to visualize various features of each chromosome, such as gene position and distribution, the density of all SSR markers, and the distribution of the tested SSR markers.

In silico evaluation and PCR validation of SSR markers

Four whole genome sequences, including those of two cultivated peanut landraces, A. hypogaea cv. Shitouqi [26] and Tifrunner [27], and their two diploid progenitors, A. duranensis and A. ipaensis [20], were used as templates for in silico evaluation of the designed SSR markers using electronic PCR (e-PCR) with the following parameters: 4-bp mismatch, 1-bp gap, and 0–2000-bp product size [40]. The paired primers should meet the following criteria: (1) only one in silico PCR product generated from any one template; (2) the SSR basic motif was as expected; and (3) the PCR product length was also as expected.

A total of 188 newly developed SSR primer pairs were arbitrarily selected and then synthesized for validation by PCR amplification. The PCR mixture was prepared in a 10-μl volume that contained 1 μl template DNA (~ 100 ng), 0.4 μl of each primer, 3.2 μl ddH2O, and 5 μl 2× Power Taq PCR Master Mix (BioTeke Corporation, Beijing, China). The PCR amplification program was as follows: (1) 95 °C for 5 min; (2) total of 35 cycles, each cycle was as below: 30 s at 95 °C, 30s at 55 °C for annealing, 30 s at 72 °C for extension; (3)10 min at 72 °C, and then storage at 4 °C. PCR products were separated and tested in Fragment Analyzer™ Automated CE System (Advanced Analytical Technologies, Inc., Beijing, China), and PROSize 2.0 was used to analyze and visualize the data (https://www.aati-us.com/support/software/).

Integrating new and previously published SSR markers

In a previous study, we constructed an integrated consensus genetic map with published markers, including EST-SSR, transcriptome-SSR, and g-SSR markers [16]. Based on the forward and reverse primer sequences of each marker, these public SSR markers were mapped to the A. hypogaea cv. Fuhuasheng reference genome for integration with the newly developed SSR markers using e-PCR software with default parameters [40].

Availability of data and materials

The datasets supporting the conclusions of this article are included within the article and its additional files.

Abbreviations

ESTs:

Expressed sequence tags

PCR:

Polymerase chain reaction

QTL:

Quantitative trait locus

SSR:

Simple sequence repeats

References

  1. Krishna G, Singh BK, Kim EK, Morya VK, Ramteke PW. Progress in genetic engineering of peanut (Arachis hypogaea L.)—a review. Plant Biotechnol J. 2015;13:147–62.

    Article  CAS  Google Scholar 

  2. Pandey MK, Wang ML, Qiao L, Feng S, Khera P, Wang H, et al. Identification of QTLs associated with oil content and mapping FAD2 genes and their relative contribution to oil quality in peanut (Arachis hypogaea L.). BMC Genet. 2014;15:133.

    Article  Google Scholar 

  3. Chen Y, Ren X, Zheng Y, Zhou X, Huang L, Yan L, et al. Genetic mapping of yield traits using RIL population derived from Fuchuan Dahuasheng and ICG6375 of peanut (Arachis hypogaea L.). Mol Breed. 2017;37:17.

    Article  Google Scholar 

  4. Shasidhar Y, Vishwakarma MK, Pandey MK, Janila P, Variath MT, Manohar SS, et al. Molecular mapping of oil content and fatty acids using dense genetic maps in groundnut (Arachis hypogaea L.). front. Plant Sci. 2017;8:794.

    Google Scholar 

  5. Mondal S, Badigannavar AM. Mapping of a dominant rust resistance gene revealed two R genes around the major rust_QTL in cultivated peanut (Arachis hypogaea L.). Theor Appl Genet. 2018;131:1671–81.

    Article  CAS  Google Scholar 

  6. Varshney RK, Mohan SM, Gaur PM, Gangarao NV, Pandey MK, Bohra A, et al. Achievements and prospects of genomics-assisted breeding in three legume crops of the semi-arid tropics. Biotechnol Adv. 2013;31:1120–34.

    Article  Google Scholar 

  7. Pandey MK, Roorkiwal M, Singh VK, Ramalingam A, Kudapa H, Thudi M, et al. Emerging genomic tools for legume breeding: current status and future prospects. Front Plant Sci. 2016;7:455.

    PubMed  PubMed Central  Google Scholar 

  8. Gupta PK, Balyan HS, Sharma PC, Ramesh B. Microsatellites in plants: a new class of molecular markers. Curr Sci. 1996;70:45–54.

    CAS  Google Scholar 

  9. Haq SU, Jain R, Sharma M, Kachhwaha S, Kothari SL. Identification and characterization of microsatellites in expressed sequence tags and their cross transferability in different plants. Int J Genomics. 2014; 2014:863948. See: https://www.ncbi.nlm.nih.gov/pubmed/25389527.

  10. Guo B, Chen X, Hong Y, Liang X, Dang P, Brenneman T, et al. Analysis of gene expression profiles in leaf tissues of cultivated peanuts and development of EST-SSR markers and gene discovery. Int J Plant Genomics. 2009;715605. https://doi.org/10.1186/1471-2229-9-35.

  11. Liang X, Chen X, Hong Y, Liu H, Zhou G, Li S, et al. Utility of EST derived SSR in cultivated peanut (Arachis hypogaea L.) and Arachis wild species. BMC Plant Biol. 2009;9:35.

    Article  Google Scholar 

  12. Zhang J, Liang S, Duan J, Wang J, Chen S, Cheng Z, et al. De novo assembly and characterization of the transcriptome during seed development, and generation of genic-SSR markers in Peanut (Arachis hypogaea L.). BMC Genomics. 2012;13:90.

    Article  CAS  Google Scholar 

  13. Huang L, Wu B, Zhao J, Li H, Chen W, Zheng Y, et al. Characterization and transferable utility of microsatellite markers in the wild and cultivated Arachis species. PLoS One. 2016;11:e0156633.

    Article  Google Scholar 

  14. Zhong R, Zhou M, Zhao C, Hou L, Li C, Wang X, et al. SSR marker development from peanut gynophore transcriptome sequencing. Plant Breed. 2016;135:111–7.

    Article  CAS  Google Scholar 

  15. Shirasawa K, Bertioli DJ, Varshney RK, Moretzsohn MC, Leal-Bertioli SC, Thudi M, et al. Integrated consensus map of cultivated peanut and wild relatives reveals structures of the a and B genomes of Arachis and divergence of the legume genomes. DNA Res. 2013;20:173–84.

    Article  CAS  Google Scholar 

  16. Lu Q, Liu H, Hong YB, Li H, Liu HY, Li XY, et al. Consensus map integration and QTL meta analysis narrowed a locus for yield traits to 0.7 cM and refined a region for late leaf spot resistance traits to 0.38 cM on linkage group A05 in peanut (Arachis hypogaea L.). BMC Genomics. 2018;19:887.

    Article  CAS  Google Scholar 

  17. Cho YG, Ishii T, Temnykh S, Chen X, Lipovich L, McCouch SR, et al. Diversity of microsatellites derived from genomic libraries and GenBank sequences in rice (Oryza sativa L.). Theor Appl Genet. 2000;100:713–22.

    Article  CAS  Google Scholar 

  18. Kochert G, Halward T, Branch WD, Simpson CE. RFLP variability in peanut (Arachis hypogaea L.) cultivars and wild species. Theor Appl Genet. 1991;81:565–70.

    Article  CAS  Google Scholar 

  19. Moretzsohn MC, Gouvea EG, Inglis PW, Leal-Bertioli SC, Valls JF, Bertioli DJ. A study of the relationships of cultivated peanut (Arachis hypogaea) and its most closely related wild species using intron sequences and microsatellite markers. Ann Bot. 2013;111:113–26.

    Article  CAS  Google Scholar 

  20. Bertioli DJ, Cannon SB, Froenicke L, Huang G, Farmer AD, Cannon EK, et al. The genome sequences of Arachis duranensis and Arachis ipaensis, the diploid ancestors of cultivated peanut. Nat Genet. 2016;48:438–46.

    Article  CAS  Google Scholar 

  21. Chen X, Li H, Pandey MK, Yang Q, Wang X, Garg V, et al. Draft genome of the peanut A-genome progenitor (Arachis duranensis) provides insights into geocarpy, oil biosynthesis, and allergens. Proc Natl Acad Sci U S A. 2016;113:6785–90.

    Article  CAS  Google Scholar 

  22. Lu Q, Li H, Hong Y, Zhang G, Wen S, Li X, et al. Genome sequencing and analysis of the peanut B-genome progenitor (Arachis ipaensis). Front Plant Sci. 2018;9:1099.

    Article  Google Scholar 

  23. Yin D, Ji C, Ma X, Li H, Zhang W, Li S, et al. Genome of an allotetraploid wild peanut Arachis monticola: a de novo assembly. Gigascience. 2018;7:1–9.

  24. Zhao C, Qiu J, Agarwal G, Wang J, Ren X, Xia H, et al. Genome-wide discovery of microsatellite markers from diploid progenitor species, Arachis duranensis and A. ipaensis, and their application in cultivated peanut (A. hypogaea). Front Plant Sci. 2017;8:1209.

    Article  Google Scholar 

  25. Chen X, Lu Q, Liu H, Zhang J, Hong Y, Lan H, et al. Sequencing of cultivated peanut, Arachis hypogaea, yields insights into genome evolution and oil improvement. Mol Plant. 2019;12:920–34.

    Article  CAS  Google Scholar 

  26. Zhuang W, Chen H, Yang M, Wang J, Pandey MK, Zhang C, et al. The genome of cultivated peanut provides insight into legume karyotypes, polyploid evolution and crop domestication. Nat Genet. 2019;51:865–76.

    Article  CAS  Google Scholar 

  27. Bertioli DJ, Jenkins J, Clevenger J, Dudchenko O, Gao D, Seijo G, et al. The genome sequence of segmental allotetraploid peanut Arachis hypogaea. Nat Genet. 2019;51:877–84.

    Article  CAS  Google Scholar 

  28. Powell W, Machray GC, Provan J. Polymorphism revealed by simple sequence repeats. Trends Plant Sci. 1996;1:215–22.

    Article  Google Scholar 

  29. Song GQ, Li MJ, Xiao H, Wang XJ, Tang RH, Xia H, et al. EST sequencing and SSR marker development from cultivated peanut (Arachis hypogaea L.). ElectrJ Biotechnol. 2010;13:7–8.

    Google Scholar 

  30. Wang H, Penmetsa RV, Yuan M, Gong L, Zhao Y, Guo B, et al. Development and characterization of BAC-end sequence derived SSRs, and their incorporation into a new higher density genetic map for cultivated peanut (Arachis hypogaea L.). BMC Plant Biol. 2012;12:10.

    Article  CAS  Google Scholar 

  31. Zhao T, Wang J, Xu G, Weng B. Advance on peanut breeding of high-yielding and good-quality and the application of biological technology. Chin J Trop Crop. 2011;32:2187–95 (in Chinese with English abstract).

    Google Scholar 

  32. Pandey MK, Wang H, Khera P, Vishwakarma MK, Kale SM, Culbreath AK, et al. Genetic dissection of novel QTLs for resistance to leaf spots and tomato spotted wilt virus in peanut (Arachis hypogaea L.). Front Plant Sci. 2017;8:25.

    PubMed  PubMed Central  Google Scholar 

  33. Hong YB, Liang XQ, Chen XP, Liu HY, Zhou GY, Li SX, et al. Construction of Genetic Linkage Map Based on SSR Markers in Peanut (Arachis hypogaea L.). Agr Sci China. 2008;7:915–21.

    Article  CAS  Google Scholar 

  34. Varshney RK, Bertioli DJ, Moretzsohn MC, Vadez V, Krishnamurthy L, Aruna R, et al. The first SSR-based genetic linkage map for cultivated groundnut (Arachis hypogaea L.). Theor Appl Genet. 2009;118:729–39.

    Article  CAS  Google Scholar 

  35. Khedikar YP, Gowda MV, Sarvamangala C, Patgar KV, Upadhyaya HD, Varshney RK. A QTL study on late leaf spot and rust revealed one major QTL for molecular breeding for rust resistance in groundnut (Arachis hypogaea L.). Theor Appl Genet. 2010;121:971–84.

    Article  CAS  Google Scholar 

  36. 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 Crop Res. 2011;122:49–59.

    Article  Google Scholar 

  37. Sujay V, Gowda MVC, Pandey MK, Bhat RS, Khedikar YP, Nadaf HL, 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 Breeding. 2012;30:773–88.

    Article  CAS  Google Scholar 

  38. Deng P, Wang M, Feng K, Cui L, Tong W, Song W, et al. Genome-wide characterization of microsatellites in Triticeae species: abundance, distribution and evolution. Sci Rep. 2016;6:32224.

    Article  CAS  Google Scholar 

  39. Untergasser A, Cutcutache I, Koressaar T, Ye J, Faircloth BC, Remm M, et al. Primer3—new capabilities and interfaces. Nucleic Acids Res. 2012;40:e115.

    Article  CAS  Google Scholar 

  40. Schuler GD. Electronic PCR: bridging the gap between genome mapping and genome sequencing. Trends Biotechnol. 1998;16:456–9.

    Article  CAS  Google Scholar 

Download references

Acknowledgments

Not applicable.

Funding

This work was supported by the National Natural Science Foundation of China (31771841, 31801401), the National key research and development program (2018YFD0201009), the Science and Technology Program of Guangzhou, China (201904010190), the Open Fund of Guangdong Key Laboratory of Crop Genetic Improvement (2017B030314090), the Natural Science Foundation of Guangdong Province (2017A030311007), the Modern Agro-industry Technology Research System (CARS-13), the teamwork projects funded by Guangdong Natural Science Foundation (No.2017A030312004), “Jinying Star” (2018) and “Youth Associate Research Fellow” (2019) Talent Project of Guangdong Academy of Agricultural Sciences. These funding bodies had no role in the design of the study, extraction, analysis, and interpretation of data, or in writing the manuscript.

Author information

Authors and Affiliations

Authors

Contributions

XPC and XQL designed the experiments. QL, YBH, JNZ and HFL1 analyzed the data. HL and HFL2 developed SSR markers. HYL, XYL and SW performed PCR experiments. QL and YBH and HFJ drafted the manuscript. SXL, GYZ and Varshney RK revised the manuscript, respectively. All authors have read and approved the final manuscript. 1 Haofa Lan; 2 Haifen Li.

Corresponding authors

Correspondence to Xiaoping Chen or Xuanqiang Liang.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1:

Table S1. Number of SSRs and SSR markers identified on different chromosomes. Table S2. Summary of SSR motifs and repeats.

Additional file 2:

Figure S1. Chromosome-wide distribution of SSRs in A. hypogaea cv. Fuhuasheng genome. Figure S2. Number of SSR repeat motifs. Figure S3. Abundance of the top 30 different types of SSR motifs. Figure S4. Number (A) and percentage (B) of different types of SSR markers. Figure S5. Summary of SSR types of the developed SSR markers. Figure S6. Product size of 188 SSR markers tested by PCR amplification. Figure S7. Number of loci in the public SSR markers as determined by e-PCR remapping. Figure S8. Comparison of known QTLs in genetic and physical maps.

Additional file 3:

Table S3. Summary of SSR markers in subgenome A01–A05.

Additional file 4:

Table S4. Summary of SSR markers in subgenome A06–A10.

Additional file 5:

Table S5. Summary of SSR markers in subgenome B01–B05.

Additional file 6:

Table S6. Summary of SSR markers in subgenome B06–B10.

Additional file 7:

Table S7. Summary of SSR markers on nine scaffolds. Table S8. Summary of common SSR markers in four Arachis species. Table S9. Summary of 188 SSR markers tested by PCR amplification. Table S10. Summary of the public SSR markers with only a single locus match to the reference genome. Table S11. Candidate gene identification and annotation for the target QTL, qLLS_T13_A05_7.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Lu, Q., Hong, Y., Li, S. et al. Genome-wide identification of microsatellite markers from cultivated peanut (Arachis hypogaea L.). BMC Genomics 20, 799 (2019). https://doi.org/10.1186/s12864-019-6148-5

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12864-019-6148-5

Keywords