Utility of EST-derived SSR in cultivated peanut (Arachis hypogaea L.) and Arachiswild species

Type and frequency of peanut EST-SSRs

Primer design and validation

EST-SSR polymorphism

In the present study, 251 valid EST-SSR primer pairs were used for assessment of the polymorphism among cultivated and wild Arachis species. Within cultivated peanuts, 26 (10.3%) EST-SSRs exhibited polymorphism (Table 3). A total of 55 alleles were detected and the average number of alleles per SSR marker was 2.1 with a range of 2–4 alleles based on the dominant scoring of the SSR bands characterized by the presence or absence of a particular band (Additional File 1). The PIC values ranged from 0.09 to 0.69 with an average value of 0.33. The greatest variation of SSR alleles was found for EM-78, which interacted with 4 alleles in 22 cultivated peanuts genotypes and the PIC value was 0.69.

The polymorphism of 251 cultivated peanut-derived EST-SSR in 16 accessions of wild species was evaluated. The results showed that 221 of 251 EST-SSR loci (88%) were polymorphic (Table 3), with a total of 867 alleles (Additional File 1). Allelic diversity was estimated for those polymorphic EST-SSR markers. The number of alleles detected among 16 wild species ranged from 2 to 9, with an average of 3.9 alleles per locus (Additional File 1). A maximum of 9 alleles were observed for primer EM-71. The PIC values varied between 0.594 and 0.820 with an average value of 0.721.

Sequence comparison of SSR bands

For further understanding of the EST-SSR polymorphism at the nucleotide level, the amplified products of primer EM-31 from two genotypes of cultivated peanuts and three accessions of wild species were cloned and sequenced (Figure 1, Figure 2). All the sequenced alleles from both cultivars and wild species were highly identical to the original locus (EST sequence) from which the EST-SSR marker EM-31 was mined. Sequence alignment showed that all the primer-binding regions are highly conserved. Allelic diversity could be attributed mainly to differences in repeat type and length in the microsatellite regions, although some variations such as repeat number or insertions of additional motifs were observed in the microsatellite regions. In addition, a few single base substitutions were observed in the microsatellite flanking regions. Out of them, one occurred in A. cardenasii, one in A. duranensis, and two in A. pintoi.

Frequency and distribution of EST-SSRs

The frequency of SSRs in SSR-ESTs more accurately reflects the density of SSRs in the transcribed region of the genome. However, random sequencing within cDNA libraries usually resulted in a high proportion of redundant ESTs. In this study, to reduce the dataset size and avoid overestimation of the EST-SSR frequency, SSR search were performed following redundancy elimination. A total of 11,432 potential unique EST sequences (about 6.4 Mb) were used for SSR search and 6.8% (780) of ESTs contained specified SSR motifs, generating 881 unique SSRs. This is a relatively higher abundance of SSRs for peanut ESTs, compared to the previous reports for maize (1.4%), barley (3.4%), wheat (3.2%), soyghum (3.6%), rice (4.7%) [18], Medicago truncatula (3.0%) [23] and wild Arachis species [34]. The different abundance of SSRs was known to be dependent on the SSR search criteria, the size of the dataset, the database-mining tools and different species [22]. In this work, the frequency of occurrence for EST-derived SSRs was one EST-SSR in every 7.3 kb. In previous reports, an EST-SSR occurs every 13.8 kb in Arabidopsis thaliana, 3.4 kb in rice, 8.1 kb in maize, 7.4 kb in soybean, 11.1 kb in tomato, 20.0 kb in cotton and 14.0 kb in poplar [37]. The variations of frequencies among different studies were mainly due to the criteria used to identify SSR in the database mining.

In earlier reports, tri-nucleotide repeats were generally the most common motif found in both monocots [22] and dicots [23]. During the process of mining EST-SSRs in the various plant species, tri-nucleotide was also observed to be most frequent [26], regardless of the EST-SSR search criteria. Until now, only one report described that di-nucleotide repeats were most abundant followed by tri- or mono-nucleotide repeats in dicots [38]. In the present investigation, tri-nucleotide repeat was found to be abundant followed by di-nucleotide. In term of single SSR motif, the di-nucleotide motif (AG/TC) _n was highest frequent [18, 39]. Among the di-nucleotide motifs, the two most dominant motif types were AG and AT, representing an average frequency of 24.7% and 6.4%, respectively. This was in agreement with recent studies in cultivated peanut (Arachis hypogaea L.) [35] and wild Arachis species [34]. In this work, the AAG with 17.4% of frequency following di-nucleotide motif AG was the most abundant in the ten tri-nucleotide motifs. In other plant species, the most frequent tri-nucleotide repeat motifs were (AAC/TTG) _n in wheat, (AGG/TCC) _n in rice, (CCG/GGC) _n in maize, (AAG/TTC) _n in soybean, and (CCG/GGC) _n in barley and sorghum [18, 19, 39, 40]. The previous studies of Arabidopsis [37] and soybean [24] also suggested that the tri-nucleotide AAG motif may be common motif in dicots. In contrast, the abundance of the tri-nucleotide CCG repeat motif was favored overwhelmingly in cereal species [18, 19, 32] and also considered as a specific feature of monocot genome, which may be due to increasing the G + C content [26].

Validation and polymorphism of EST-SSR markers

In this study, a total of 290 designed primer pairs were used for validation of the EST-SSR markers. Of these, 251 (86.5%) yielded amplicons in both cultivated peanut and wild species. This result was similar to previous studies in which a success rate of 60–90% amplification has been reported [21, 25, 40–42]. In those studies, they also reported a similar success rate of amplification for both genomic SSRs and EST-SSRs. However, EST-SSRs were reported to be less polymorphic than genomic SSRs in crop plants due to greater DNA sequence conservation in transcribed regions [17, 28, 43–46]. Previous studies highlighted the fact that EST-SSR markers have higher transferability and better applicability than genomic SSR markers [17, 47–49]. In addition to high transferability, EST-SSRs were good candidates for the development of conserved orthologous markers for genetic analysis and breeding of different species [22]. Pervious reports showed that the transferability of EST-SSRs from one species to another ranged from 40–89% [21, 23, 24, 27, 29, 40, 41, 50, 51]. Our results indicated that 100% of EST-SSR amplifiable primers for cultivated peanut can produce amplicons in Arachis wild species.

In the present investigation, the mean percentage of polymorphic loci of EST-SSR markers was 9.96% in cultivated peanuts. This value was lower than those of genomic SSR found in earlier studies [7, 12, 47], but higher than the percentage of polymorphic loci in cultivated peanut observed using RAPD (6.6%) [5] and AFLP (6.7%) [4]. No major difference was observed in terms of allele numbers and PIC values for the EST-SSR markers among the cultivated genotypes, while significant difference was observed among wild species. Therefore, the low level of EST-SSR polymorphism detected in cultivated peanuts may be compensated by their higher potential for cross-species transferability to wild species. In the present study, 100% transferability of EST-SSR with 86.6% polymorphism from cultivated peanut to wild Arachis species was observed. The value is higher than that of genomic SSR cross-transferability [10]. The high level of transferability indicated that these markers would be highly effective for molecular study of Arachis species. Since current molecular markers display a low level of genetic polymorphism in cultivated peanuts [2–4, 6, 52, 53], it is difficult to construct a high-density genetic linkage map for cultivated peanut which could be used in breeding programs. However, a genetic map constructed using wild species together with transferable molecular markers derived from cultivated peanuts would contribute to understanding the introgression of genes from wild species to cultivated peanuts [10, 13]. Therefore, the development of a set of transferable EST-SSR markers from cultivated peanuts will be a great benefit to construct a high-density genetic map of wild species. The map would allow the identification of markers, especially transferable EST-SSR markers, associated with resistance or other agronomic traits in wild species, and in turn, help to discover corresponding markers or genes in cultivated peanuts.

Additionally, a comparison of sequences of cross-species amplicons generated by primer EM-31 further confirmed the conservation and transferability of the developed EST-SSR loci. In general, the amplified regions were found to be similar to the original peanut EST sequences from which the SSRs were developed and their comparisons across species (Figure 2) correlated the observed 'cross-species alleles' precisely with the expected length variations. Furthermore, in addition to the variation of the number of SSR repeat, the allele sequences also indicated that a few additional point mutation in the SSR motifs flanking regions. Similar variation has been reported in earlier studies [39, 47, 54, 55]. This phenomenon is supposed to be the innate evolving nature of the genome, and thus can be indicative of the evolutionary relationships of the tested taxa [47].

Plant materials and DNA extraction

In the present study, twenty-two accessions of cultivated peanut(A.hypogaea L.) corresponding to two subspecies (hypogaea and fastigiata) and sixteen accessions of wild species from seven sections of the genus Arachis were used (Additional File 2). The leaf samples of each accession were collected from Peanut Germplasm Bank located in Crops Research Institute, Guangdong Academy of Agriculture, Guangzhou, China. the genomic DNA was extracted as described by Sharma [56].

Data mining for SSR marker

A total of 24,238 EST sequences including 20,160 developed by Guo et al (2008)[57] and 4078 retrieved from National Center of Biotechnology Information (NCBI) were used in this study. These ESTs were assembled using the TGICL program [36]. A Perl script known as MIcroSAtellite (MISA http://pgrc.ipk-gatersleben.de/misa/.) was used to mine microsatellites. In this work, SSRs were considered for primer design that fitted the following criteria: a minimum length of 14 bp, excluding polyA and polyT repeat, at least 7 repeat units in case of di-nucleotide and at least 5 repeat units for tri-, tetra-, penta- and hexa-nucleotide SSRs. Therefore, the paired numbers representing SSR motif length and the minimum repeat number in the MISA configuration file (misa.ini) were modified to 2–7, 3–5, 4–5, 5-5 and 6-5 (mono-type excluded).

Primer design and PCR amplification

Using Primer Premier 5 program (Whitehead Institute for Biomedical Research, Cambridge, Mass), primers were designed based on the following core criteria: (1) melting temperature (Tm) between 52°C and 63°C with 60°C as optimum; (2) product size ranging from 100 bp to 350 bp; (3) primer length ranging from 18 bp to 24 bp with amplification rate larger than 80%; (4) GC% content between 40% and 60%. The parameters were modified when unsuitable primer pairs were retrieved by the program. PCR analysis was performed in a total volume of 20 μl with the following cycling profile: 1 cycle of 5 min at 94°C, an annealing temperature of 55°C for 35 cycles (1 min at 94°C, 30 s at 55°C, 45 s at 72°C) and an additional cycle of 10 min at 72°C. Each of the primer pairs was screened twice to confirm the repeatability of the observed bands in each genotype. PCR products were separated on 6% polyacrylamide denaturing gels. The gels were silver stained for SSR bands detection.

Sequencing of PCR bands

The SSR alleles amplified in two cultivars and three wild species for EM-31 primer were individually cloned and sequenced. PCR amplification products were separated by 6% polyacrylamide gel and target allele bands were excised and dipped in 10 μl of nuclease free water for 30 min. Another round of PCR was made following the same protocol with recycled DNA as template. The second-round PCR products were separated in a 2% agarose gel and the target band was purified using TIANGEN Gel Extracting Kit (TIANGEN Inc. Beijing China). The purified PCR fragment from agarose gel was cloned using the Takara TA cloning kit pMD-18 (Takara, Dalian, China). The ligation product was transformed into competent Escherichia coli cells. The positive clones identified by PCR were sequenced by Invitrogen Company. The final edited sequences belonging to different genotypes were compared with the original SSR containing EST sequence using Omiga program [58], and the exported multiple sequence alignment was modified by Genedoc http://www.nrbsc.org/gfx/genedoc/index.html.

Data scoring and statistical analysis

The allelic and genotypic frequencies were calculated for the samples analyzed. The genetic diversity of the samples as a whole was estimated based on the number of alleles per locus (total number of alleles/number of loci), the percentage of polymorphic loci (number of polymorphic loci/total number of loci analyzed) and polymorphism information content (PIC). The polymorphism was determined according to the presence or absence of the SSR locus. The value of PIC was calculated using the formula

where P _i is the frequency of an individual genotype generated by a given EST-SSR primer pair and summation extends over n alleles.