Quillaja saponaria Molina is an endemic Chilean tree, known as soapbark, soap bark tree, or quillay (as it is known in Chile). It is an evergreen tree of the family Quillajaceae (Kubitzki, 2007), found from the Coquimbo Region to Arauco Province in the Bío-Bío Region, approximately between 31° and 38° south (García and Ormazabal, 2008). It grows from sea level to 1600 m a.s.l., preferably in dry areas that are poor in nutrients. The family is monotypic with a single genus and two species from warm-temperate South America (Chile, Brazil, and northern Argentina). Quillaja saponaria is considered one of the most important and representative species of the sclerophyllous forest from central Chile, and these communities are part of a biodiversity hotspot called the Chilean Winter Rainfall–Valdivian Forest (Mittermeier et al., 1998; Myers et al., 2000). Moreover, Q. saponaria is a timber species (García and Ormazabal, 2008), its flowers are melliferous (Montenegro et al., 2009), and the bark is rich in saponins and medicinal adjuvants (Kensil et al., 1991; San Martín and Briones, 1999).
Our purpose is to evaluate the effects of anthropogenic fragmentation on the patterns of genetic variation and connectivity in populations of Q. saponaria. For this reason, we isolated and characterized 12 nuclear microsatellite loci that are being successfully applied to describe spatial patterns of genetic structure. These are the first microsatellite markers developed for a Quillaja Molina species.
METHODS AND RESULTS
Microsatellite isolation was performed by the simple sequence repeat (SSR) development company Genetic Marker Services (Brighton, United Kingdom; www.geneticmarkerservices.com). Genomic DNA was extracted from a single Q. saponaria (Qsa) individual collected in the locality of Coya, near the city of Rancagua, O'Higgins Region, Chile (locality 5 in Appendix 1), using a modified cetyltrimethylammonium bromide (CTAB) protocol described by Doyle and Doyle (1990), and used to develop an enriched library to isolate microsatellite (SSR)–containing loci. Enrichment involved incubating adapter-ligated restricted DNA with filter-bonded synthetic repeat motifs: (AG)17, (AC)17, (AAC)10, (CCG)10, (CTG)10, and (AAT)10. Thirty-two motif-positive Escherichia coli JM109 clones were detected and sequenced, of which 23 contained exploitable repeat motifs with sufficient flanking regions to design forward/reverse primer pairs. The online primer design software Primer3 (Rozen and Skaletsky, 1999) was used to develop primer pairs amplifying fragments ranging in size from 100 to 250 bp, to help minimize later multiloading overlap ambiguities during sequencer genotyping. The primers were then tested for successful amplification on one individual from each of seven populations, chosen to represent the whole latitudinal range of the species' distribution (Cuesta el Espino, Cuesta Los Cristales, Santa Marta, Fundo La Rosa, La Gatera, San Fabián de Alico, and Santa Bárbara; Appendix 1), using a touchdown PCR protocol. PCR amplifications were performed in a 25-µL final volume containing 7 pmol of each primer, 1.5 mM of MgCl2, 0.2 mM of each dNTP, 1× PCR buffer, 0.8 µg/µL bovine serum albumin (BSA), 0.5 units Taq polymerase (Promega Corporation, Madison, Wisconsin, USA), and 1.5 µL of DNA diluted 20-fold. For all loci, the same touchdown PCR program was used, consisting of 32 cycles with denaturation at 95°C for 60 s; annealing for 60 s with temperature decreasing 1°C every two cycles from 64°C to 59°C (12 cycles), then 10 cycles at 58°C, and 10 cycles at 57°C; elongation at 72°C for 60 s; and a final extension at 72°C for 5 min. Products were checked for specificity and polymorphism on a cooled high-resolution agarose gel, consisting of 4% MetaPhor (Lonza, Basel, Switzerland) agarose in TAE, run in a cold room at 10°C. Sixteen loci showed clear and specific bands with size variation among the seven individuals assayed (Table 1). Of these, the 12 with the highest degree of apparent polymorphism were chosen for fluorescent labeling (Table 1).
To determine the number of alleles, expected and observed heterozygosity, and to test for deviations from Hardy–Weinberg equilibrium (HWE) and linkage disequilibrium, we tested the 12 fluorescent-labeled primer pairs on 39 individuals from three different populations, situated in the northern and southern parts of the species' range (i.e., Angostura, Astillero, and Cuesta Batuco; localities 1, 8, and 9 in Appendix 1). The PCR amplifications were performed with the Platinum Multiplex PCR Master Mix (Applied Biosystems, Carlsbad, California, USA) using three mixes, each with four fluorescent-labeled primers (Applied Biosystems, Warrington, United Kingdom; Table 1): Mix 1 (Qsa7, Qsa11, Qsa13, Qsa24), Mix 2 (Qsa16, Qsa17, Qsa18, Qsa70), and Mix 3 (Qsa19, Qsa26, Qsa28, Qsa66). The PCR amplifications were performed in a 5-µL final volume containing final concentrations of 1× Platinum Multiplex PCR Master Mix, 2 mM of additional MgCl2, 0.10 µM of primer mix, and 20 ng of template DNA, with the following cycling conditions in all cases: an initial heat activation at 94°C for 5 min; followed by a touchdown PCR consisting of 32 cycles with denaturation at 95°C for 60 s; annealing for 60 s with temperature decreasing 1°C every two cycles from 64°C to 59°C (12 cycles), then 10 cycles at 58°C and 10 cycles at 57°C; elongation at 72°C for 60 s; and a final extension at 72°C for 5 min. The PCR products (2 µL) were mixed with 10 µL of Hi-Di formamide and 0.3 µL of GeneScan 600 LIZ Size Standard before analysis on an automatic ABI-PRISM 3100-Avant sequencer (Applied Biosystems, Carlsbad, California, USA). Fragment analysis and final sizing were performed using Peak Scanner software (version 1.0; Applied Biosystems, Carlsbad, California, USA).
Characteristics of 16 microsatellite loci developed in Quillaja saponaria.
The number of alleles, expected and observed heterozygosity, and the significance of the deviations from HWE and linkage equilibrium were estimated with Arlequin 188.8.131.52 (Excoffier et al., 2005). The average of the total number of alleles per locus was 5.3, with a range between two (Qsa24) and 13 (Qsa19). Significant deviations from expectations under HWE after a Bonferroni correction for multiple comparisons were found in three loci from population Astillero, and in two loci from Cuesta Batuco, probably due to the small sample sizes or to the presence of null alleles, particularly in the case of the heterozygote deficits observed at loci Qsa26 and Qsa66 in population Astillero and at Qsa66 in Cuesta Batuco (Table 2). No instances of significant linkage disequilibrium were found among pairs of loci after Bonferroni correction.
The 12 microsatellite loci characterized here for Q. saponaria are the first developed for this species and genus. These markers will be useful in studies of the population genetic diversity and structure of the species for the purpose of conservation under a scenario of land use and climate change.
Genetic properties of 12 newly developed microsatellites of Quillaja saponaria.a
 The authors thank Chile Abroad Grants program (Doctorate level, Call 2010), Comisión Nacional de Investigación Científica y Tecnológica (CONICYT, Chile), II Concurso Fondo de Investigación del Bosque Nativo (Project no. 041-2011), Corporación Nacional Forestal (CONAF, Chile), and the Postgraduate Program in Biological Sciences of Universidad Nacional Autónoma de México, fulfilling the requirement for the degree of Doctor in Sciences of the Biological Sciences postgraduate program.