Carex helodes Link (sect. Spirostachyae (Drejer) L. H. Bailey, Cyperaceae) is a diploid, wind-pollinated, perennial herb with a minimum generation time of two years. The species is endemic to the western Mediterranean, being locally distributed in southern Portugal and southwestern Spain, and rare in northern Morocco (Escudero et al., 2008a). This sedge occurs in temporarily inundated acidic soils in open cork oak woodlands. Despite its well-characterized morphology, C. helodes has been misidentified as C. laevigata Sm. by some authors (see Luceño et al., 2009). Recent cytotaxonomic and nuclearand plastid-based phylogenetic studies have revealed the monophyly of C. helodes populations and its taxonomic independence within sect. Spirostachyae (Escudero et al., 2008a, 2008b; Escudero and Luceño, 2009). Carex helodes is an endangered species in Spain because the extent of its severely fragmented occurrence is less than 100 km2, with continued loss of area and habitat quality (Moreno, 2008; Bañares et al., 2010).
Our aim is to develop molecular markers for further studies of gene flow among and within populations. Nuclear microsatellites have been proven to be highly variable and very suitable to the study of recent gene flow between populations (Ouborg et al., 1999). To accomplish our task, we isolated and characterized 109 nuclear microsatellites.
METHODS AND RESULTS
We extracted genomic DNA using a DNeasy Plant Mini Kit (QIAGEN, Valencia, California, USA). We used ∼5 µg from one C. helodes individual collected in Madroñalejo (Aznalcóllar, Seville, Spain; see Appendix 1 for GPS coordinates and voucher specimens) to construct a shotgun genomic library that was sequenced on l/4th of a plate using 454 GS FLX Titanium chemistry (Roche Applied Science, Indianapolis, Indiana, USA) at the University of Arizona Genetic Core (Tucson, Arizona, USA). We generated 108.3 Mb of quality-filtered data, distributed over 221,198 unique reads with an average length of 490 bp after quality filtering (quality score [Q] ≥ 20 using a 10-bp sliding window). We searched for all possible microsatellite loci containing at least six perfect repeats for hexa-, penta-, tetra-, and trinucleotides or 12 perfect repeats for dinucleotides and designed primers using the software QDD version 3 (Meglécz et al., 2014). We used the unique reads as input to detect microsatellite sequences. The reads were used to build contigs using QDD version 3 (default options were used: sequence set limit of 80 bp, 95% minimum identity between two sequences to make a consensus, and 66% as the proportion of sequences that must have the same base on the aligned site to accept it as a consensus). For primer design, the default options were also used but the minimum size of the PCR product was set to 100 bp and the maximum to 450 bp.
We found a total of 3985 microsatellite loci, including 51 hexa-, 58 penta-, 78 tetra-, 406 tri-, and 3392 dinucleotide loci. We selected 27 hexa-, 26 penta-, 23 tetra-, 206 tri-, and 152 dinucleotide loci that met our criteria (at least 12 repeats for dinucleotides and six for the rest) and tested a total of 132 loci. Specifically, we tested the eight hexa-, eight penta-, eight tetra-, 18 tri-, and 90 dinucleotide loci with the highest numbers of repeats (Table 1).
For primer testing, DNA was isolated from silica gel–dried leaves using a modified cetyltrimethylammonium bromide (CTAB) extraction method (Milligan, 1998) that included tissue grinding in a Mixer Mill MM301 (Retsch GmbH, Haan, Germany) and resuspension in TLE buffer (10 mM Tris-HCl [pH 8.0], 0.1 mM EDTA). We sampled a total of 64 individuals from five different populations (Appendix 1). Initially, 19 individuals from the five populations were sampled and genotyped for 109 loci. Nine individuals from the Madroñalejo population, four individuals from the Barraçao-Caldas de Monchique population, one individual from the Cortelha population, four individuals from population 1 (Pop. 1) at Ksar el Kebir, and one individual from population 2 (Pop. 2) at Ksar el Kebir (Appendix 1). Finally, 38 additional individuals from the Madroñalejo population (Spain) and seven additional individuals from the Cortelha population (Portugal, Table 2) were genotyped for the 34 most polymorphic loci. The 34 most polymorphic loci were estimated based on the initial screening using 19 individuals from five populations.
Characteristics of 109 microsatellite markers isolated from populations of Carex helodes.
PCR amplifications were performed in a 20-µL final volume containing 1× buffer (67 mM Tris-HCl [pH 8.8], 16 mM (NH4)2SO4, 10 mM KCl, 0.01% stabilizer), 2.5 mM MgCl2, 0.01% bovine serum albumin (BSA; Roche Diagnostics, Mannheim, Germany), 0.25 mM dNTP, 0.40 µM dye-labeled M13 primer, 0.40 µM PIG-tailed reverse primer, 0.04 µM M13-tailed forward primer (see M13 and PIG-tail sequences in Table 1), 0.5 units Taq DNA polymerase (Bioline, London, United Kingdom), and approximately 50–70 ng genomic DNA. Reactions were undertaken in a touchdown PCR protocol in a Bio-Rad DNA Engine Peltier Thermal Cycler (Bio-Rad Laboratories, Hercules, California, USA), with an initial 2 min of denaturation at 94°C; 17 cycles at 92°C for 30 s, annealing at 60–44°C for 30 s (1°C decrease in each cycle), and extension at 72°C for 30 s; 25 cycles at 92°C for 30 s, 44°C for 30 s, and 72°C for 30 s; and a final extension of 5 min at 72°C. PCR products were labeled using FAM, VIC, NED, or PET dyes (Applied Biosystems, Foster City, California, USA) on an additional 19-bp M13 primer (5′-CACGACGTTGTAAAACGAC-3′) according to the methods of Boutin-Ganache et al. (2001). Moreover, a palindromic sequence tail (5′-GTGTCTT-3′) was added to the 5′ end of the reverse primer to improve adenylation and facilitate genotyping. Amplified fragments were analyzed on an ABI 3130x1 Genetic Analyzer (Applied Biosystems) and sized using GeneMapper 4.0 and GeneScan 500 LIZ Size Standard (Applied Biosystems). No multiplexing was attempted at the PCR stage.
From a total of 132 loci tested, 18 were monomorphic, 14 showed complex or nonspecific amplification, and nine failed to amplify. The remaining 91 loci were polymorphic in Portugal and Morocco (Table 1). Sequences of 109 loci (18 monomorphic and 91 polymorphic) from a shotgun genomic library were submitted to the National Center for Biotechnology Information (NCBI) Sequence Read Archive (SRA; accession SRP062192). We observed a total of 275 alleles for our initial sampling of C. helodes (nine individuals from Madroñalejo, five from Portugal, and five from Morocco), averaging overall 3.23 alleles per locus (SD = 1.15). Gametic disequilibrium (GENEPOP 4.1.4; Rousset, 2008) and the presence of null alleles according to the Oosterhout method (MICRO-CHECKER 2.2.3; van Oosterhout et al., 2004) were checked using Bonferroni-corrected P values (P < 0.05 / 34 = 0.0015) to assess the significance of the results obtained. No significant gametic disequilibrium was detected for any pair of loci (P > 0.01), and four loci (Cahe587, CaheWE1, CaheGVB, and CaheQWV) showed signs of the presence of null alleles suggested by an excess of homozygotes. Some loci showed one to three alleles with a size difference of one to two base pairs relative to contiguous alleles (see Table 1). Nevertheless, no variation outside of the microsatellite region was shown by these loci, and no scoring errors were detected by MICRO-CHECKER (van Oosterhout et al., 2004). These polymorphisms could result from nonstepwise mutations in the repeat array.
The 47 individuals from the Madroñalejo population (Spain) showed an identical genotype (nine individuals were genotyped using 109 loci, and 38 additional individuals were genotyped using the 34 most variable loci). This seems to indicate that this population has probably suffered a recent founder or bottleneck event.
We used Arlequin 22.214.171.124 (Excoffier and Lischer, 2010) to calculate number of alleles, observed (Ho) and expected heterozygosities (He), and to test for deviations from Hardy–Weinberg equilibrium (HWE; Table 2) for the 34 most polymorphic loci in the eight individuals from the Cortelha population (Portugal). The eight individuals from the Cortelha population had a mean Ho of 0.357 (SD = 0.292) and a mean He of 0.384 (SD = 0.255). Six loci were monomorphic in this population, and only one (CaheWE1, P < 0.05; Table 2) of the 28 remaining loci deviated from HWE (Table 2).
This study provides 109 simple sequence repeat markers to quantify the degree of genetic diversity in the endangered C. helodes. These markers may also help disentangle the phylogeographic history of C. helodes as well as gene flow among populations. This next-generation sequencing approach may also be useful for linkage mapping studies of experimental crosses between populations in a nonmodel organism. However, the low observed levels of variation within populations are not enough to genotype reproductive individuals and their seeds within populations and study gene interchange patterns within populations.
Results of screening 34 loci in the Cortelha population of Carex helodes.a
 The authors are grateful for logistical support provided by the Molecular Ecology Laboratory, Estación Biológica de Doñana (LEM-EBD-CSIC) for primer testing. This study was supported by a postdoctoral 2011 Marie Curie International Outgoing Fellowship for Career Development (IOF; HoloChromEvol 301119) to M.E. and a Junta de Andalucía grant (RNM-5731) to P.J., as well as a Severo Ochoa Excellence Award (SEV-2012-0262).