- Short Report
- Open Access
Microsatellite marker development for the rubber tree (Hevea brasiliensis): characterization and cross-amplification in wild Hevea species
BMC Research Notes volume 5, Article number: 329 (2012)
The rubber tree (Hevea brasiliensis) is native to the Amazon region and it is the major source of natural rubber in the world. Rubber tree breeding is time-consuming and expensive. However, molecular markers such as microsatellites can reduce the time required for these programs. This study reports new genomic microsatellite markers developed and characterized in H. brasiliensis and the evaluation of their transferability to other Hevea species.
We constructed di- and trinucleotide-enriched libraries. From these two libraries, 153 primer pairs were designed and initially evaluated using 9 genotypes of H. brasiliensis. A total of 119 primer pairs had a good amplification product, 90 of which were polymorphic. We chose 46 of the polymorphic markers and characterized them in 36 genotypes of H. brasiliensis. The expected and observed heterozygosities ranged from 0.1387 to 0.8629 and 0.0909 to 0.9167, respectively. The polymorphism information content (PIC) values ranged from 0.097 to 0.8339, and the mean number of alleles was 6.4 (2–17). These 46 microsatellites were also tested in 6 other Hevea species. The percentage of transferability ranged from 82% to 87%. Locus duplication was found in H. brasiliensis and also in 5 of other species in which transferability was tested.
This study reports new microsatellite markers for H. brasiliensis that can be used for genetic linkage mapping, quantitative trait loci identification and marker- assisted selection. The high percentage of transferability may be useful in the evaluations of genetic variability and to monitor introgression of genetic variability from different Hevea species into breeding programs.
Hevea brasiliensis (Willd. ex Adr. de Juss.) Muell. -Arg., native to the Amazon rainforest, is a diploid (2n = 36, x = 9), perennial, monoecious and cross-pollinated tree species. It belongs to the genus Hevea and the botanical family Euphorbiaceae.
The genus Hevea comprises 11 inter-crossable species [1, 2] (H. benthamiana Muell.- Arg., H. brasiliensis H. carmagoana Pires, H. camporum Ducke, H. guianensis Aubl, H. microphylla Ule, H. nitida Mart. ex-Muel.-Arg., H. pauciflora (Spruce ex-Benth.) Muell.-Arg., H. rigidifolia (Spruce ex-Benth.) Muell.-Arg., H. spruceana (Benth.) Muell.-Arg. and H. paludosa Ule), which have evolved in the Amazon rainforest over 100,000 years .
Of all the species in the genus Hevea H. brasiliensis is the most economically important, because it is the major source of natural rubber worldwide. Natural rubber is important mainly in the tire industry but also in many other sectors because it is flexible, resistant, impermeable to liquids and abrasion resistant . These singular properties make natural rubber both complementary and competitive to synthetic rubber and furthermore superior to it in varied applications. As is the case for many other perennial trees, rubber tree breeding is time-consuming and expensive. An average duration of 20–25 years of field experiments in large areas is generally required to obtain a new cultivar with reasonably low risks .
Molecular markers, such as amplified fragment length polymorphisms (AFLPs), restriction fragment length polymorphisms (RFLPs), random amplified polymorphic DNA (RAPD), simple sequence repeats (SSRs or microsatellites) and single nucleotide polymorphisms (SNPs), have successively become increasingly important in plant breeding. These markers are efficient tools for the assessment of genetic diversity, the identification of quantitative trait loci (QTLs) and/or gene mapping, variety protection and marker-assisted selection (MAS) .
Of these markers, microsatellite markers are considered the most suitable for genetic studies, because they combine co-dominance and high polymorphism with abundance, locus specificity and uniform dispersion in plant genomes. Moreover, microsatellite markers can discriminate closely related individuals . In addition, microsatellite analysis is inexpensive, high reproducible and highly transferable across related species.
Following biochemical markers such as isozymes , molecular markers have been developed and used since the middle of 1990’s for diversity studies [9, 10], genetic mapping  and the identification of genetic loci implicated in the expression of agronomic traits in H. brasiliensis[11–14]. However, most of the markers used have been isozyme, RAPD, RFLP or AFLP markers.
To date, few studies using genomic microsatellites [10, 15, 16] or microsatellites from expressed sequence tags (EST-SSRs) have been published . Accordingly, we present the development of genomic microsatellites using dinucleotide- and trinucleotide- enriched libraries, the characterization of these microsatellite markers in multiple accessions of H. brasiliensis and test their transferability in six other Hevea species.
Microsatellite-enriched library analysis
Di- and trinucleotide-enriched libraries were constructed, and the clones from each library were sequenced (576 and 288 clones, respectively).
A total of 291 (50.5%) clones from the dinucleotide-enriched library contained microsatellite sequences. Since dinucleotide probes were used for the library enrichment, these motifs were the most abundant comprising 324 (90.2%) SSRs followed by 20 tetranucleotide (5.6%), 12 trinucleotide (3.3%) and 3 pentanucleotide (0.8%) (Table 1).
In the trinucleotide-enriched library, 62 (21.5%) clones contained microsatellites sequences. Trinucleotides were the most frequent motif with 45 (72.5%) SSRs followed by 13 dinucleotide (20.9%), 9 tetranucleotide (14.5%) and 3 pentanucleotide (4.8%) SSRs (Table 1).
For some sequences, both the dinucleotide- and trinucleotide-enriched libraries had more than one microsatellite, which explains the greater number of SSRs found compared to the number of sequences analyzed.
A total of 130 and 32 primer pairs were designed based on the di- and trinucleotide-enriched libraries, respectively. To remove possible redundancies with published SSRs, each sequence containing SSR was compared against GenBank database using BLASTN. Nine sequences, only from the dinucleotide library, were identical to previously published sequences, which already had primer pairs. These sequences were removed from this study.
Dinucleotide motifs have been reported as the most abundant type of microsatellite in plant genomes . Recent studies relating to expressed sequence tags revealed that trinucleotide motifs are the most abundant motifs in ESTs in many plants, such as sugarcane , barley , grapes , rice , wheat  and citrus , whereas in other plants dinucleotide motifs are the most abundant in ESTs, such as kiwi , coffee  and apricot and peach  In rubber tree, Feng and co-workers analyzed all of the ESTs in NCBI database in search for microsatellites. They found that dinucleotide motifs were the most abundant and were three times more abundant than trinucleotides . The low efficiency of the trinucleotide enrichment in our study could result from the low frequency of trinucleotide motifs in the rubber tree genome.
Polymorphism analysis and cross-species transferability
In total, 153 specific primer pairs were designed. A total of 119 primer pairs produced good amplification products (Additional file 1: Table S1) and 90 of these products were polymorphic among a set of 9 rubber clones indicated in Table 2.
The observed and expected heterozygosities ranged from 0.1111 to 1 and 0.1111 to 0.9150, respectively, and the PIC values ranged between 0.0994 and 0.8496. The mean number of alleles was 4.46 (2–9 alleles).
Of the 90 polymorphic markers we chose 46 SSRs for characterization among a set of 36 genotypes of H. brasiliensis (Table 2). The PIC values of these markers ranged from 0.097 to 0.8339, and the observed and expected heterozygosities ranged from 0.0909 to 0.9167 and 0.1387 to 0.8629, respectively. The mean number of alleles was 6.4 (2–17 alleles) (Table 3).
Six other species from the genus Hevea (H. guianensis, H. rigidifolia, H. nitida, H. pauciflora, H. benthamiana and H. camargoana) being two different genotypes of H. pauciflora, were used to evaluate the transferability of the markers (Table 2). All loci were tested under the same PCR conditions used for H. brasiliensis.
Of the 46 loci tested, 40 (87%) were amplified for H. guianensis and H. pauciflora- (112CNSG), 39 (85%) were amplified for H. camargoana, H. nitida and H. pauciflora-(116CNSG) and 38 (82%) were amplified for H. benthamiana (Table 4).
For most of the SSR loci, which cross-species amplification were tested, the number of alleles found in the H. brasiliensis genotypes in conjunction with the tested species (from 2 to 21, with a mean number of 9.5) was larger when the same loci were analyzed in H. brasiliensis alone (Table 3) which means that other species revealed some novel alleles.
Three primers pairs, HB54, HB62 and HB69, did not produce amplification products for the six other species tested, whereas 27 loci were amplified for all species (Table 4).
Saha and co-workers first observed SSR cross-amplification in H. benthamiana and H. spruceana. Together with the cross-fertility potential, this high SSR transferability supports the consideration of the Hevea genus as a species complex with moderate differentiation among the species. These aspects appear to be favorable for genetic introgressions using other Hevea species to the rubber breeding population, which is mainly based on H. brasiliensis.
Members of the Euphorbiaceae family have a basic number of chromosomes between 6 and 11; thus, any species with more than this number of chromosomes could be amphidiploid in origin [28, 29]. Similar to cassava, plants in the genus Hevea have 36 chromosomes and behave as diploid. In these species, it has been assumed that n = 18 and the basic number of chromosomes is x = 9 [30, 31].
In different Hevea species, chromosomes mainly formed bivalents and tetravalents are rarely produced as a result of pairing between non-homologous chromosomes during prophase I and metaphase I of meiosis [31, 32].
Moreover, cytogenetic studies revealed two distinct loci on two different chromosomes, bearing the same18S-5.8 S-25 S rDNA sequence which may have arisen by the hybridization of two unknown diploid species (n = 9), thus suggesting a possible allotetraploid origin, however no potencial diploid ancestor has been described to date .
Locus duplication, as revealed by molecular markers, had been reported in H. brasiliensis 11] and H. guianensis H. rigidifolia and H. pauciflora.
In H. brasiliensis, we observed locus duplication for the markers HB45 and HB109 for clones RRIM 701 and Fx4098, respectively. Cases of locus duplication were observed in all Hevea species tested, with the exception of H. rigidifolia. We also observed locus duplication of the HB36, HB68, HB100 and HB105 markers in H. nitida; HB68 and HB105 markers in H. benthamiana; HB101 and HB105 in H. guianensis. The marker HB105 exhibited one duplicated locus for H. pauciflora (112 CNSG), H. pauciflora (116 CNSG) and H. camargoana.
Although we have not been mapped the loci we cited as duplicated in the wild species, we observed more than 2 alleles (excluding the stutter bands) in the same base pair range compared with the other genotypes used for characterization (Figure 1). The presence of allele duplication, as determined by polyacrylamide gel electrophoresis suggests that the loci are duplicated for the species involved.
Although locus duplication had been described in H. brasiliensis, H. guianinsis, H. rigidifolia and H. pauciflora, this is the first report for H. benthamiana, H. camargoana and H. nitida.
Herein, we report the development of new SSR markers for H. brasiliensis, representing a powerful resource for genetic diversity studies and genetic breeding techniques, such as molecular genetic mapping, QTLs identification and MAS.
Due to observed high percentage of cross-amplification and the absence of reproductive barriers between the species within this genus, these markers may be important tools to monitor the genetic variability from other Hevea species into the current breeding programs.In addition, these SSR markers can be helpful for the identification of important agronomic characteristics in different Hevea species.
Plant material and DNA extraction
We used 36 genotypes of H. brasiliensis to characterize the SSRs (Table 2). These genotypes were kindly provided by the Agronomic Institute of Campinas (IAC) and Michelin Plantation (Brasil). We also used 6 other species of the genus Hevea (H. guianensis H. rigidifolia H. nitida H. pauciflora H. benthamiana and H. camargoana) that were kindly provided by the Brazilian Agricultural Research Corporation (EMBRAPA - Amazônia Ocidental) in Manaus, AM to test the transferability of the SSRs (Table 2). The genomic DNA samples were extracted from lyophilized leaf tissues using a modified CTAB method  and their quality and quantity were assessed using 1% agarose gel electrophoresis.
Construction of microsatellite-enriched libraries and sequence analysis
The microsatellite-enriched libraries for H. brasiliensis were constructed using the RRIM 600 clone according to the methodology described by Billote and co-workers . The DNA samples were digested with AFAI and enriched using streptavidin magnetic-coated beads (Streptavidin MagneSphere Paramagnetic Particles, Promega, Madison, WI) and (CT)8 and (GT)8 biotinylated microsatellite probes for the dinucleotide-enriched library or (ATC)8 and (CCT)8 for the trinucleotide-enriched library.
Selected DNA fragments were amplified by PCR and then cloned into the pGEM-T vector (Promega, Madison, WI). Competent XL1-blue Escherichia coli cells were transformed with the recombinant plasmids and cultivated on agar medium containing ampicillin and 100 μg/ml of X-galactosidase. The clones containing the insert were sequenced using the Big Dye Terminator v3.1 Cycle Sequencing Kit (Applied Biosystems, Foster City, CA) and an automated ABI 377 sequencer (Applied Biosystems, Foster City, CA).
All the sequences obtained were aligned, edited and eliminated if redundant using SeqMan (DNASTAR, Madison, WI). The sequences were also evaluated using Microsat software (A. M. Risterucci, CIRAD, personal communication), which removes the adapters and verifies the presence of restriction sites.
The microsatellites were identified using a specific research tool, SSRIT – (Simple Sequence Repeat Identification Tool) . As a criterion for the SSR selection, the sequences that showed at least five dinucleotide repeats; four trinucleotide repeats; and three tetra-, penta- and hexanucleotide repeats were selected. Primers complementary to the sequences flanking the microsatellites were designed using Primer Select Program (DNAStar, Madison, WI) and Primer 3 . To eliminate possible redundancies with the published SSRs, each sequence containing microsatellite was compared against the GenBank database using BLASTN. The redundant sequences which already had primer pairs published were eliminated.
The microsatellite fragments were PCR amplified in a 15 μl reaction containing 25 ng template DNA, 0.5 μM each primer, 100 μM each dNTP, 1.5 mM MgCl2, 20 mM Tris–HCl, 50 mM KCl and 0.5 U Taq DNA Polymerase. The PCR amplifications were performed with the following four programs: (1) initial denaturation at 94°C for 4 min, 30 amplification cycles (1 min at 94°C, 45 s at the specific annealing temperature and 1 min at 72°C), and a final extension at 72°C for 10 min; (2) TD-1, initial denaturation for 3 min at 94°C, 10 amplification cycles with a 0.5°C decrease in annealing temperature per cycle starting at 55°C or 60°C (94°C for 1 min, 60°C or 55°C for 30 s, and 72°C for 1 min 15 s); followed by 20 cycles with annealing at 50°C (94°C for 1 min, 50°C for 30 s and 72°C for 30 s) and a final elongation step at 72°C for 10 min; (3) TD-2, previously described by Le Guen and co-workers ; and (4) TD-3, initial denaturation at 94°C for 2 min, 2× 10 cycles with a 1°C decrease in annealing temperature per cycle from 65°C to 55°C (94°C for 1 min, 65°C for 1 min and 72°C for 1 min), followed by 18 cycles at 55°C (94°C for 1 min, 55°C for 1 min and 72°C for 2 min), and a final elongation step at 72°C for 5 min.
The amplification products were resolved by electrophoresis through 3% agarose gels prior to vertical electrophoresis using 6% denaturing polyacrylamide gels and were subsequently silver stained . The product sizes were determined by comparison with a 10 bp DNA ladder (Invitrogen, Carlsbad, CA).
The allelic polymorphic information content of each SSR was calculated using the formula, where n is the number of alleles of the marker among the set of genotypes used for characterizing the SSR polymorphism, and pi and pj are the frequencies of alleles i and j. The observed and expected heterozygosities were calculated using the TFPGA program .
Pires JM, Secco RS, Gomes JI: Taxonomia e fitogeografia das seringueiras Hevea spp. 2002, Belém: Embrapa Amazônia Oriental, 103-
Gonçalves P, Cardoso M, Ortolani AA: Origem, variabilidade e domesticação da Hevea; uma revisão. Pesquisa agropecuária brasileira. 1990, 25 (2): 135-156.
Clément-Demange A, Priyadarshan PM, Thuy Hoa TT, Venkatachalam P: Hevea Rubber Breeding and Genetics. Plant Breeding Reviews. Edited by: Janick J. 2007, New Jersey: John Wiley & Sons, Inc, ch 4-Volume 29
Cornish K: Similarities and differences in rubber biochemistry among plant species. Phytochem. 2001, 57: 1123-1134. 10.1016/S0031-9422(01)00097-8.
Webster CC, Baulkwill WJ: Rubber Longman edn. 1989, Longman Singapore Publishers Ltd
Feng SP, Li WG, Huang HS, Wang JY, Wu YT: Development, characterization and cross-species/genera transferability of EST-SSR markers for rubber tree (Hevea brasiliensis). Mol Breeding. 2009, 23: 85-97. 10.1007/s11032-008-9216-0.
Kalia RK, Rai MK, Kalia S, Singh R, Dhawan AK: Microsatellite markers: an overview of the recent progress in plants. Euphytica. 2011, 177: 309-334. 10.1007/s10681-010-0286-9.
Chevallier MH: Genetic variability of Hevea brasiliensis germplasm using isozyme markers. J Nat Rubber Res. 1988, 3: 42-53.
Besse P, Seguin M, Lebrun P, Chevallier MH, Nicolas D, Lanaud C: Genetic diversity among wild and cultivated populations of Hevea brasiliensis assessed by nuclear RFLP analysis. Theor Appl Genet. 1994, 88: 199-207.
Le Guen V, Gay C, Xiong TC, Souza LM, Rodier-Goud M, Seguin M: Development and characterization of 296 new polymorphic microsatellite markers for rubber tree (Hevea brasiliensis). Plant Breeding. 2011, 130: 294-296. 10.1111/j.1439-0523.2010.01774.x.
Lespinasse D, Rodier-Goud M, Grivet L, Leconte A, Legnate H, Seguin M: A saturated genetic linkage map of rubber tree (Hevea spp.) based on RFLP, AFLP, microsatellite, and isozyme markers. Theor Appl Genet. 2000, 100: 127-138. 10.1007/s001220050018.
Le Guen V, Garcia D, Mattos CRR, Doaré F, Lespinasse D, Seguin M: Bypassing of a polygenic Microcyclus ulei resistance in rubber tree, analyzed by QTL detection. New Phytologist. 2007, 173: 335-345. 10.1111/j.1469-8137.2006.01911.x.
Le Guen V, Lespinasse D, Oliver G, Rodier-Goud M, Pinard F, Seguin M: Molecular mapping of genes conferring field resistance to South American Leaf Blight (Microcyclus ulei) in rubber tree. Theor Appl Genet. 2003, 108: 160-167. 10.1007/s00122-003-1407-9.
Lespinasse D, Grivet L, Troispoux V, Rodier-Goud M, Pinard F, Seguin M: Identification of QTLs involved in the resistance to South American leaf blight (Microcyclus ulei) in the rubber tree. Theor Appl Genet. 2000, 100: 975-984. 10.1007/s001220051379.
Souza LM, Mantello CC, Santos MO, Gonçalves PO, Souza AP: Microsatellites from rubber tree (Hevea brasiliensis) for genetic diversity analysis and cross amplification in six Hevea wild species. Cons Genet Res. 2009, 1: 75-79. 10.1007/s12686-009-9018-7.
Yu F, Wang BH, Feng SP, Wang JY, Li WG, Wu YT: Development, characterization, and cross-species/genera transferability of SSR markers for rubber tree (Hevea brasiliensis). Plant Cell Rep. 2010, 30: 335-344.
Morgante M, Hanafey M, Powell W: Microsatellites are preferentially associated with nonrepetitive DNA in plant genomes. Nat Genet. 2002, 30: 194-200. 10.1038/ng822.
Marconi TG, Costa EA, Miranda HRCAN, Mancini MC, Cardoso-Silva CB, Oliveira KM, Pinto LR, Mollinari M, Garcia AAF, Souza AP: Functional markers for gene mapping and gnetic diversity studies in sugarcane. BMC Research Notes. 2011, 4: 264-10.1186/1756-0500-4-264.
Scott KD, Eggler P, Seaton G, Rossetto M, Ablett EM, Lee LS, Henry RJ: Analysis of SSRs derived from grape ESTs. Theor Appl Genet. 2000, 100: 723-726. 10.1007/s001220051344.
Thiel T, Michalek W, Varshney RK, Graner A: Exploiting EST databases for the development and characterization of genederived SSR-markers in barley (Hordeum vulgare L.). Theor Appl Genet. 2003, 106: 411-422.
Han Z, Wang C, Song X, Guo W, Gou J, Li C, Chen X, Zhang T: Characteristics, development and mapping of Gossypium hirsutum derived EST-SSRs in allotetraploid cotton. Theor Appl Genet. 2006, 112: 430-439. 10.1007/s00122-005-0142-9.
La Rota M, Kantety RV, Yu JK, Sorrells ME: Nonrandom distribution and frequencies of genomic and EST-derived microsatellite markers in rice, wheat, and barley. BMC Genomics. 2005, 6: 2-10.1186/1471-2164-6-2.
Chen C, Zhou P, Choi YA, Huang S, Gmitter FG: Mining and characterizing microsatellites from citrus ESTs. Theor Appl Genet. 2006, 112: 1248-1257. 10.1007/s00122-006-0226-1.
Fraser LG, Harvey CF, Crowhurst RN, De Silva HN: EST-derived microsatellites from Actinidia species and their potential for mapping. Theor Appl Genet. 2004, 108: 1010-1016. 10.1007/s00122-003-1517-4.
Aggarwal RK, Hendre PS, Varshney RK, Bhat PR, Krishnakumar V, Singh L: Identification, characterization and utilization of EST-derived genic microsatellite markers for genome analyses of coffee and related species. Theor Appl Genet. 2007, 114: 359-372. 10.1007/s00122-006-0440-x.
Jung S, Albert A, Jesudurai C, Tomkins J, Main D: Frequency, type, distribution and annotation of simple sequence repeats in Rosaceae ESTs. Funct Integr Genomics. 2005, 5: 136-143. 10.1007/s10142-005-0139-0.
Saha T, Roy CB, Nazeer MA: Microsatellite variability and its use in the characterization of cultivated clones of Hevea brasiliensis. Plant Breeding. 2005, 124: 86-92. 10.1111/j.1439-0523.2004.01053.x.
Perry AB: Chromosome Number and Phylogenetic Relationships in the Euphorbiaceae. Am J Bot. 1943, 30: 527-543. 10.2307/2437291.
Leitch AR, Lim KY, Leitch IJ, O’Neill M, Chye M, Low F: Molecular cytogenetic studies in rubber, Hevea brasiliensis Muell. Arg. (Euphorbiaceae). Genome. 1998, 41: 464-467.
Majumder SK: Chromosome studies of some species of Hevea. J. Rubber Res. Inst. Malays. 1964, 18: 269-273.
Ong SN: Chromosome morphology at pachytene stage in Hevea brasiliensis: a preliminary report. Proceedings of the International Rubber Conference. Edited by: Sripathi B. 1975, Kuala Lumpur: Rubber Research Institue of Malaysia, 3-12.
Bouharmont J: Recherches taxonomiques et caryologiques chez quelques espèces du genre Hevea. Série Scientifique. 1960, 64-Volume 85
Doyle JJ, Doyle JL: A rapid DNA isolation procedure for small quantities of fresh leaf tissue. Phytochem Bull. 1987, 19: 11-15.
Billotte N, Lagoda P, Risterucci AM, Baurens FC: Microsatellite-enriched libraries: applied methodology for the development of SSR markers in tropical crops. Fruits. 1999, 54: 277-288.
Temnykh S, Declerck G, Lukashova A, Lipovich L, Cartinhour S, McCouch S: 2001 . Computational and experimental analysis of microsatellites in rice (Oryza sativaL.): Frequency, length variation, transposon associations, and genetic marker potential. Genome Res. 2001, 11: 1441-1452. 10.1101/gr.184001.
Untergasser A, Nijveen H, Rao X, Bisseling T, Geurts R, Leunissen JAM: Primer3Plus, an enhanced web interface to Primer3. Nucleic Acids Res. 2007, 35: W 71-W74.
Creste S, Tulmann-Neto A, Figueira A: Detection of single repeat polymorphisms in denaturing polyacrilamide sequencing gels by silver staining. Plant Mol Biol Rep. 2001, 19: 229-306.
Miller MP: Tools for population genetic analysis (TFPGA) 1.3: A Windows program for the analysis of allozyme and molecular population genetic data. 1997, Computer Software distributed by the author
The authors thank the Brazilian Agricultural Research Corporation (EMBRAPA - Amazônia Ocidental) and the Agronomic Institute for donating the analyzed Hevea germplasm. This work was supported by grants from Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP, proc. 2007/50562-4). CCM and FIS received undergraduate scholarships (2008/55974 and 2009/51757-9) and, LMS and CCM received graduate scholarships (2007/50392-1 and 2011/50188-0) from FAPESP. APS and PSG are recipient of a research fellowship from Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq).
The authors declare that they have no competing interests.
CCM and FIS developed the microsatellite-enriched libraries, performed the computational microsatellite identification, designed the flanking primers and performed the microsatellite marker validation. CCM performed the statistical analysis and drafted the manuscript. LMS participated in the design and implementation of the study and the microsatellite identification. APS and PSG conceived the study and participated in its design and coordination. APS helped to draft the manuscript. All of the authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1 : Table S1. Characterization of the developed SSR markers. The table presents the 119 SSR markers developed for H. brasiliensis, including the primers sequence, annealing temperature, number of alleles, expected size in base pair, allelic range, observed and expected heterozigosity and polymorphism information content. The nine accession indicated in Table 1 were used for the SSR characterization. (XLS 71 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.