Prosopis L. (Fabaceae) contains woody species of arid and semiarid areas of the Americas, Africa, and Asia. There is considerable taxonomic confusion about the American sect. Algarobia DC., due largely to the phenotypic variation within species and the low genetic differentiation and frequent hybridization among species (Hunziker et al., 1986).
Population genetic, conservation, and evolutionary studies using molecular markers have been generated over the past few years in Prosopis (Bessega et al., 2009, 2010). However, the lack of a sufficient number of simple sequence repeat (SSR) markers limits the conclusions that can be inferred. The only available microsatellite markers in Prosopis were six loci developed in P. chilensis (Molina) Stuntz (Mottura et al., 2005). They revealed cross-species affinity when tested in five species of sect. Algarobia and two of sect. Strombocarpa Benth. Although these markers proved to be useful loci for population studies in P. alba Griseb. (Bessega et al., 2009, 2010), at least 10 microsatellite loci are recommended to accurately estimate relatedness coefficients and diverse population parameters depending on parentage estimates (Wang et al., 2010). In the current study, we developed additional SSRs useful for population genetic analysis including effects of fragmentation on loss of diversity, fine-scale population structure, breeding system and pollen dispersal, and detection of recent bottlenecks in P. alba and P. chilensis, the most important native forest resources among Argentinean “algarrobos.” Also, these SSRs will contribute to more accurate relationships estimates, and these estimates may be used to estimate in situ heritability estimates. Taking advantage of the high genetic similarity among species of sect. Algarobia (Bessega et al., 2005) and the observation of cross-amplification in several species of Prosopis (Mottura et al., 2005), we searched for possible microsatellite loci from 454 sequence data from a bulk sample of American mesquites introduced to Australia. We then tested their utility in these species and cross-amplification in five other species of sect. Algarobia.
METHODS AND RESULTS
We extracted DNA from a bulk collection of leaf material from Queensland, Australia, which originated from American species of sect. Algarobia. The sample included P. pallida (Humb. & Bonpl. ex Willd.) Kunth, P. velutina Wooton, P. glandulosa Torr., and putative hybrids involving these species (Appendix 1). Taxonomic identification is difficult due to the high hybridization rate of these invasive populations. Three micrograms of DNA was sent to the Australian Genome Research Facility (www.agrf.org.au/; St. Lucia, Queensland, Australia). The DNA was sheared and one-quarter plate was sequenced on a Roche GS FLX (454 Life Sciences, a Roche Company, Branford, Connecticut, USA) following the manufacturer's protocol. We screened the sequences obtained for perfect di-, tri-, and tetranucleotide repeats using MSATCOMMANDER version 0.8 (Faircloth, 2008). We used the option “Design primer,” in which case the software searches for microsatellite repeats and identifies possible primer annealing sites in one step. Primer3 (Rozen and Skaletsky, 2000) is implemented in MSATCOMMANDER for primer design according to the following criteria: amplification products within the size range of 100–500 bp, optimal melting temperature (range 57–62°C), optimal GC content of 50%, possession of at least 1-bp GC clamp, low levels of self- or pair-complementary, and maximum end stability (ΔG) of 8.0 (Faircloth, 2008). The 454 sequencing from the Australian sample bulk yielded 66,100 reads, which assembled into 2672 contigs of which 795 were longer than 500 bp. We detected 760 sequences with tandem repeats: 695 di-, 31 tri-, and 10 tetranucleotide microsatellites with at least 10, eight, and eight repeats, respectively. From these, 736 sequences contained suitable primer sites.
To test the amplification success of the microsatellites, DNA was extracted from samples belonging to P. alba, P. chilensis, P. ruscifolia Griseb., P. flexuosa DC., P. kuntzei Harms, P. glandulosa, and P. pallida (Appendix 1). Due to limited economic resources, we randomly tested 30 primer pairs (GL1–GL30; Table 1). Untested sequences are planned to be explored for additional primers in the future; in the meantime, the corresponding author can be contacted for information about them. In P. alba and P. chilensis, we genotyped between eight and 18 individual samples; in P. ruscifolia, P. flexuosa, P. glandulosa, P. pallida, and P. kuntzei, we only tested microsatellite amplification (between three and five individuals were screened). To confirm that the amplified products corresponded to the microsatellites found bioinformatically in the Australian sequence data set, one allele from each locus of P. alba was sequenced (Departamento de Ecología, Genética y Evolución, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Argentina). Table 1 presents the characteristics of the 12 loci studied. The PCR amplifications were carried out in a 50-µL reaction volume containing 10–30 ng DNA, 0.6 µM each primer, 0.2 mM dNTPs, 0.3 U Taq DNA polymerase (Invitrogen, Buenos Aires, Argentina), and 1.5 mM MgCl2. A MyCycler Techne thermal cycler (Bio-Rad Laboratories, Hercules, California, USA) was used for amplifications, where the cycling profile was initial denaturation at 94°C for 5 min; followed by 30 cycles at 94°C for 45 s denaturation, primer-specific annealing temperature (58–62°C; see Table 2) for 45 s, and at 72°C for 45 s extension; and a final extension at 72°C for 10 min. PCR products were run in a 6% (w/v) Polyacrylamide gel containing 5 M urea in 1 × TBE buffer (89 mM Tris, 89 mM boric acid, 2 mM EDTA, pH 8) with a 10-bp DNA Ladder (Invitrogen) size marker. Gels were stained with silver nitrate. For all the polymorphic loci that generated clear patterns, we calculated observed, unbiased, and biased heterozygosity and the effective number of alleles using BIOSYS-2 (Swofford and Selander, 1999). The polymorphism information content (PIC) (Botstein et al., 1980) was estimated.
Characteristics of the polymorphic microsatellite markers developed for Prosopis alba and P. chilensis in this study.
From the 30 potential microsatellite markers tested, 18 failed to produce a reproducible polymorphic pattern. Twelve primer pairs revealed clear and reproducible variable loci in P. alba (Table 2). Eight of these were also polymorphic in P. chilensis, GL3 amplified only in these two species, GL6 failed to amplify only in P. kuntzei, and GL9 did not amplify in P. chilensis. Ten primer pairs amplified DNA fragments in all tested species (Table 2). Across the samples, we found between two and seven alleles, with an average effective number of alleles of 2.75 for P. alba and 2.92 for P. chilensis. Average heterozygosity estimations in P. alba were 0.57, 0.60, and 0.54 (heterozygosity [H], unbiased estimated heterozygosity [Hunb], and direct count heterozygosity [HDC], respectively). In P. chilensis, the estimations were similar, yielding average values of 0.61, 0.65, and 0.59, respectively. Individual HDC did not differ significantly from H for any locus. The PIC index describes diversity within populations and characterizes the degree of polymorphism at each locus. Across all loci, the PIC ranged from 0.22 to 0.79 in P. alba and from 0.37 to 0.77 in P. chilensis. PIC values of less than 0.25 indicate low polymorphism, values between 0.25 and 0.5 indicate average polymorphism, and values higher than 0.5 indicate a highly polymorphic locus (Botstein et al., 1980). Eight out of 12 loci in P. alba and five out of eight loci in P. chilensis were highly polymorphic.
We succeeded in developing highly polymorphic SSRs in P alba and P. chilensis. These microsatellite primer pairs are an important addition to the number of genetic markers previously available, increasing from six to 18 loci. They allow accurate estimates of genetic parameters for population structure analysis, including genetic diversity, genetic structure, mating system, pollen dispersal, relatedness among individuals, and heritability of quantitative traits. The cross-amplification observed in P. kuntzei, the most differentiated species so far studied from sect. Algarobia, suggests that the loci described might also be useful for population studies in several other species of the section.
Descriptive statistics for polymorphic microsatellites in Prosopis alba and P. chilensis and evaluation of amplification success in related species.
- C. Bessega , B. O. Saidman , and J. C. Vilardi . 2005. Genetic relations among American species of Prosopis (Leguminosae) based on enzyme markers. Genetics and Molecular Biology 28: 277–286. Google Scholar
- C. Bessega , B. O. Saidman , M. R. Darquier , M. Ewens , L. Sánchez , P. Rozenberg , and J. C. Vilardi . 2009. Consistency between marker and genealogy-based heritability estimates in an experimental stand of Prosopis alba (Leguminosae). American Journal of Botany 96: 458–465. Google Scholar
- C. Bessega , B. O. Saidman , M. R. Darquier , M. Ewens , and J. C. Vilardi . 2010. Accuracy of dominant markers for estimation of relatedness and heritability in an experimental stand of Prosopis alba (Leguminosae). Tree Genetics and Genomes 7: 103–115. Google Scholar
- D. Botstein , R. D. White , M. Skolnick , and R. W. Davis . 1980. Construction of a genetic linkage map in man using restriction fragment length polymorphisms. American Journal of Human Genetics 32: 314–331. Google Scholar
- B. C. Faircloth 2008. MSATCOMMANDER: Detection of microsatellite repeat arrays and automated, locus-specific primer design. Molecular Ecology Resources 8: 92–94. Google Scholar
- J. H. Hunziker , B. O. Saidman , C. A. Naranjo , R. A. Palacios , L. Poggio , and A. D. Burghardt . 1986. Hybridization and genetic variation of Argentine species of Prosopis (Leguminosae, Mimosoideae). Forest Ecology and Management 16: 301–315. Google Scholar
- M. C. Mottura , R. Finkeldey , A. R. Verga , and O. Gailing . 2005. Development and characterization of microsatellite markers for Prosopis chilensis and Prosopis flexuosa and cross-species amplification. Molecular Ecology Notes 5: 487–489. Google Scholar
- M. Nei 1978. Estimation of average heterozygosity and genetic distance from a small number of individuals. Genetics 89: 583–590. Google Scholar
- S. Rozen , and H. J. Skaletsky . 2000. Primer3 on the WWW for general users and for biologist programmers. In S. Misener and S. A. Krawetz [eds.], Methods in molecular biology, vol. 132: Bioinformatics methods and protocols, 365–386. Humana Press, Totowa, New Jersey, USA. Google Scholar
- S. L. Swofford , and R. B. Selander . 1999. BIOSYS-2: A computer program for the analysis of allelic variation in genetics. University of Illinois, Champaign, Illinois, USA. Google Scholar
- J. Wang , P. Brekke , E. Huchard , L. A. Knapp , and G. Cowlishaw . 2010. Estimation of parameters of inbreeding and genetic drift in populations with overlapping generations. Evolution; International Journal of Organic Evolution 64: 1704–1718. Google Scholar
 This research was supported by funding from the Universidad de Buenos Aires (20020100100008) and Consejo Nacional de Investigaciones Científicas y Técnicas (PIP 11220090100147) granted to B.O.S. and J.C.V. R.W. and J.T.M. acknowledge the Taxonomy Research and Information Network (TRIN), which was funded by the Australian Commonwealth Environment Research Facilities program.