The genus Ruta L. (Rutaceae) includes nine species of perennial shrubs with three species endemic to the Canary Islands: R. oreojasme Webb. & Berthel. in Gran Canaria, R. pinnata L. in La Palma and Tenerife, and R. microcarpa Svent. in La Gomera. The three Canarian endemics have a monophyletic origin, R. montana Mill. being the sister species of the Canarian group (Salvo et al., 2010). All Canarian species are characterized by small population size and a highly fragmented distribution. Both R. oreojasme and R. microcarpa are listed in the Red List of Spanish vascular flora (Moreno, 2008) as endangered and critically endangered, respectively.
Gaining a better understanding of the genetic structure of the three Canarian endemics is fundamental for preserving these species. Within a conservation genetics project, microsatellite markers were newly developed for R. oreojasme to assess the genetic structure of populations of this species. The applicability of the markers has also been tested for the other two Canarian endemics and for the closely related R. montana.
METHODS AND RESULTS
Voucher specimens were deposited at the herbarium of the University of Zurich (Z) and at the herbarium of the Botanical Garden of Gran Canaria (LPA). Total genomic DNA was isolated from dried leaves using the DNeasy Plant Mini Kit (QIAGEN, Hombrechtikon, Switzerland) according to the manufacturer's protocol. Because the plants generated very viscous cell lysate, minor modifications were applied to the protocol to optimize results. Modifications included an increase in volume of buffer AP1 (from 400 to 600 µL), buffer AP2 (from 130 to 200 µL), and RNase A (from 4 to 6 µL), as well as a longer incubation time with buffer AP1 (at 15 min) for cell lysis.
DNA isolated from one specimen of R. oreojasme (voucher no. Z-85501) was used by Genetic Marker Services (Brighton, United Kingdom; http://www.geneticmarkerservices.com) to develop an enriched library and to design and test microsatellite primer pairs. 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. Eighteen positive library colonies were selected for sequencing, from which 11 microsatellites were designed and tested. The primer sets were designed using the primer design software Primer3 (Rozen and Skaletsky, 2000), and were developed to amplify products ranging from 100 to 300 bp to help minimize later overlap ambiguities during multiloading genotyping projects. Each primer pair was tested for specificity, polymorphism, and cross-species amplification on high-resolution agarose gels before continuing to fluorescent-labeled genotyping. Screening of microsatellites was carried out using eight individuals of R. oreojasme from three different populations. PCRs were performed in 25 µL and contained approximately 50 ng of DNA, 5 pmol of each unlabeled primer, 1.5 mM of MgCl2, 0.2 mM of each dNTP, 1×PCR buffer, and 0.5 U of SupraTherm DNA Polymerase (GeneCraft, Cologne, Germany). All PCRs were carried out under the following conditions: 60 s denaturation at 95°C, followed by 25 cycles of 95°C for 60 s, annealing at 55°C for 60 s and 72°C for 60 s, and a final extension at 72°C for 5 min. Ten microsatellite markers were easily amplified. The markers were then thoroughly screened by amplification with fluorescent-labeled primers and separation of PCR products by capillary electrophoresis. The screening was carried out using 30 individuals of R. oreojasme from three different populations: Arteara (27°50′36.7″N, 15°34′02″W; voucher no. Z-85501; 10 individuals), Montaña La Gorra (27°49′18.2″N, 15°36′07.6″W; voucher no. LPA-28995; 10 individuals), and Barranco de Las Palmas (27°50′47.3″N, 15°30′39.1″W; voucher no. LPA-29007; 10 individuals). Amplifications were performed following the two-step method described by Schuelke (2000). PCR was carried out in a final volume of 25 µL using the following components: approximately 20 ng of genomic DNA, 2.5 µL of 10× reaction buffer, 0.5 µL of each dNTP (10 mM), 1 µL of MgCl2 (50 mM), 0.2 µL of the forward primer with M13(-21 ) tail at the 5′ end (10 µM), 0.5 µL of the reverse primer (10 µM), 0.5 µL of the fluorescently labeled M13(-21) primer (FAM, NED, VIC, PET; 10 µM), and 0.1 µL of Taq DNA polymerase (5 U/µL; Bioline GmbH, Luckenwalde, Germany). PCRs were performed in a T1 Thermocycler (Biometra GmbH, Goettingen, Germany) under the following conditions: initial denaturation at 94°C for 3 min followed by 30 cycles of 94°C for 30 s, at the annealing temperature for 45 s (see Table 1), and 72°C for 1 min. The incorporation of the fluorescently labeled M13(-21) primer occurred in the following eight cycles of 94°C for 30 s, 53°C for 45 s, and 72°C for 1 min, followed by a final extension step of 72°C for 5 min. Up to four PCR products of different primer sets were pooled for each individual and separated by capillary electrophoresis on an ABI 3130xl Genetic Analyzer (Applied Biosystems, Foster City, California, USA). Alleles were sized against the internal size standard GeneScan LIZ500 (Applied Biosystems) and scored using GeneMapper software version 4.0 (Applied Biosystems). The observed allele size of all individuals genotyped was reduced by 18 bp to account for the M13(-21) universal sequence tag (Table 1). All 10 microsatellites showed a clear amplification pattern after capillary electrophoresis. Even though R. oreojasme is a tetraploid species (Stace et al., 1993), only up to two alleles were detected for individuals in all populations and for all loci, suggesting that a diploidization had occurred in their genome. For this reason, the species was treated as diploid for statistical analyses. The program GENEPOP 4.0 (Rousset, 2008) was used to calculate the observed and expected heterozygosity at each locus in the three populations (Table 2). The total number of observed alleles per locus per population ranged from 1 (locus RO78 in populations Montaña La Gorra and Barranco de Las Palmas) to 11 (locus RO62 in populations Montaña La Gorra and Barranco de Las Palmas, locus RO66 in population Montaña La Gorra; Table 2). Observed heterozygosity per locus ranged from 0.00 (locus RO78 in populations Montaña La Gorra and Barranco de Las Palmas) to 0.91 (locus RO78 in population Montaña La Gorra; Table 2). Gene diversity, inferred from Nei's heterozygosity (He), ranged from 0.00 (locus RO78 in populations Montaña La Gorra and Barranco de Las Palmas) to 0.92 (locus RO78 in population Montaña La Gorra; Table 2). Preliminary analyses suggest that these loci reveal high levels of genetic diversity in the three studied populations of R. oreojasme.
Characteristics of 10 microsatellite markers developed in Ruta oreojasme.
Results of primer screening in three populations of Ruta oreojasme. a
Cross-species amplification was tested in three related taxa using 15 individuals per species (R. microcarpa, voucher no. Z-85504; R. pinnata, voucher no. Z-85503; R. montana, voucher no. Renaux349-Z). All markers except for RO78 amplified R. microcarpa and R. pinnata samples; eight markers showed polymorphism in R. microcarpa, and seven markers showed polymorphism in R. pinnata (no data on polymorphism are available for locus RO66 in R. pinnata). Seven markers amplified R. montana samples, no data on polymorphism are available for this species. Results on cross-amplification tests are shown in Table 3.
The present investigation developed and tested 10 new unique polymorphic dinucleotide and trinucleotide nuclear DNA microsatellite loci. These markers represent useful tools for studying genetic diversity and population structure of R. oreojasme, R. microcarpa, R. pinnata, and R. montana. Information on the distribution of genetic variation in these species will contribute to their management and conservation as well as to the study of their history and evolution.
The cross-amplification of 10 microsatellite markers developed for Ruta oreojasme in three closely related species of the same genus.
- J. C. Moreno [ed.]. 2008. Lista Roja 2008 de la flora vascular española. Dirección General de Medio Natural y Política Forestal (Ministerio de Medio Ambiente, y Medio Rural y Marino, y Sociedad Española de Biología de la Conservación de Plantas), Madrid, Spain. Google Scholar
- F. Rousset 2008. GENEPOP'007: A complete re-implementation of the GENEPOP software for Windows and Linux. Molecular Ecology Resources 8: 103–106. Google Scholar
- S. Rozen , and H. Skaletsky . 2000. Primer3 on the WWW for general users and for biologist programmers In S. Krawetz and S. Misener [eds.], Methods in molecular biology, vol. 132: Bioinformatics methods and protocols, 365–386. Humana Press, Totowa, New Jersey, USA. Google Scholar
- G. Salvo , S. Y. W. Ho , G. Rosembaum , R. Ree , and E. Conti . 2010. Tracing the temporal and spatial origins of island endemics in the Mediterranean region: A case study from the Citrus family (Ruta L., Rutaceae). Systematic Biology 59: 705–722. Google Scholar
- M. Schuelke 2000. An economic method for the fluorescent labeling of PCR fragments. Nature Biotechnology 18: 233–234. Google Scholar
- H. M. Stace , J. A. Armstrong , and S. H. James . 1993. Cytoevolutionary pattern in Rutaceae. Plant Systematics and Evolution 187: 1–28. Google Scholar