Coastal plants often exhibit continuous distribution and are widely distributed. Coastal plant populations have been affected not only by Quaternary climatic oscillations but also by modern sea currents (Westberg and Kadereit, 2009), and they provide an opportunity to study various aspects of the phylogeography of flowering plants. However, except for mangrove species (e.g., Dodd et al., 2002), relatively little information is available regarding the patterns, levels of gene flow, and population genetic structure among coastal plants (e.g., Kadereit et al., 2005; Takayama et al., 2008; Westberg and Kadereit, 2009). Canavalia cathartica Thouars occurs in coastal regions from South to Southeast Asia and in East Africa (Sauer, 1964), and its seeds are dispersed by sea currents. The development of highly polymorphic microsatellite markers will provide valuable insight into the population genetic structure and dispersal mechanisms of sea-dispersed plants. In this study, we report the isolation and characterization of 10 microsatellite markers for C. cathartica and their applicability to the closely related sea-dispersed species C. lineata (Thunb.) DC.
METHODS AND RESULTS
Genomic DNA was extracted from leaf samples collected at Takana, Iriomote Island, Japan (see Appendix 1 for voucher information) using a DNeasy Plant Mini Kit (QIAGEN, Hilden, Germany) according to the manufacturer's recommended protocol. An enriched genomic DNA library was constructed as described by Bloor et al. (2001), with slight modifications. The genomic DNA was digested with three restriction enzymes, EcoRVI, AluI, and SspI (Toyobo Co., Osaka, Japan), and the resulting fragments, ranging from 0.4 to 1 kb, were excised from 1% low-melting-point agarose gels and purified. The fragments were then ligated into a blunt-end adapter as described by Bloor et al. (2001). Di- or trinucleotide sequences found within the library were enriched by hybridization to (AG)10, (GT)15, and (CAG)8 biotinylated oligonucleotides bound to streptavidin-coated magnetic beads (Promega Corporation, Madison, Wisconsin, USA). Additionally, to obtain a compound microsatellite primer (Lian et al., 2006), fragments flanked by a microsatellite at one end were amplified from the above libraries using the compound SSR primer (AC)6(AG)5 and oligo A (Bloor et al., 2001). The enriched fragments and amplified fragments were then ligated into pGEM-T Easy Vector (Promega Corporation) and transformed into DH5α-competent cells (Toyobo Co.). The recombinant DNA was sequenced in both directions on an ABI PRISM 310 DNA sequencer (Applied Biosystems, Foster City, California, USA) using M13 primers (F: 5′-GTTGAAAACGACGGCCAGT-3′; R: 5′-GGAAACAGCTATGACCATGA-3′) and a DYEnamic ET Terminator Cycle Sequencing Kit (GE Healthcare, Little Chalfont, Buckinghamshire, United Kingdom). A total of 22 cloned sequences were used to design locus-specific primers with Primer3 Plus (Untergasser et al., 2007). Ten of the 22 primer pairs successfully amplified the target region. To evaluate the loci for polymorphisms, 34 and 43 leaf samples were used for C. cathartica (Ishigaki Island [N = 22] and Iriomote Island [N = 12]) and C. lineata (Nichinan, Miyazaki [N = 27] and Minami, Tokushima [N = 16]), respectively. Voucher specimens representing the sampled populations have been deposited in the herbarium at the College of Education, University of the Ryukyus (URO) (Appendix 1). PCR was performed using 3 ng of template DNA, 0.5 µM each primer, 0.25 mM each dNTP, 1.5 mM MgCl2, 10 mM Tris-HCl (pH 8.0), 50 mM KCl, and 0.25 U of Ex Taq polymerase (TaKaRa Bio Inc., Otsu, Shiga, Japan) in a final volume of 10 µL. The forward primers were labeled with one of three Beckman Coulter WellRED dyes (D2–4; Sigma-Aldrich, St. Louis, Missouri, USA). The cycling conditions were 2 min at 94°C, followed by 30 cycles of 30 s at 94°C, 30 s at the annealing temperature (Ta; Table 1), and 45 s at 72°C, with a final extension of 10 min at 72°C. The reactions were run using a PC-818S program temperature control system (Astec, Fukuoka, Japan). Genotypes were determined using CEQ 8000 Fragment Analysis software (CEQ 8000; Beckman Coulter, Pasadena, California, USA) after electrophoresis with a CEQ 8000 Genetic Analysis System (Beckman Coulter). A GenomeLab DNA Size Standard Kit (400; Beckman Coulter) was used to determine allele size. PCR using the 10 primer sets resulted in polymorphic, single-locus amplification products (Table 1).
Characteristics of the 10 microsatellites developed in Canavalia cathartica.
Results of initial primer screening in Canavalia cathartica and C. lineata.
The results of our initial screening of C. cathartica and C. lineata are summarized in Table 2. Ten loci were polymorphic, with two to 13 and two to 10 alleles for C. cathartica and C. lineata, respectively, except for CANA4 and CANA7 in the Minami population. The observed heterozygosity values ranged from 0.333 to 0.919 and from 0.074 to 0.875 for C. cathartica and C. lineata, respectively. The expected heterozygosity values ranged from 0.375 to 0.870 and 0.071 to 0.877 for C. cathartica and C. lineata, respectively. Analyses of Hardy-Weinberg equilibrium and linkage disequilibrium were performed using GENEPOP version 3.3 (Raymond and Rousset, 1995). No significant heterozygote deficiency (P < 0.05) was detected for any of the loci except for CANA5 in the Iriomote Island population. For each population, no significant linkage disequilibrium was found among the 45 possible pairwise comparisons (P = 0.05 after sequential Bonferroni correction [Rice, 1989]).
The 10 microsatellite markers for C. cathartica revealed high levels of polymorphism, suggesting that they can be used to analyze the population genetic structure and genetic diversity of the species. Moreover, the ability to cross-amplify C. lineata sequences indicates that these markers may be used in the crossamplification of other Canavalia species.
- P. A. Bloor , F. S. Barker , P. C. Watts , H. A. Noyes , and S. J. Kemp . 2001. Microsatellite libraries by enrichment. Website http://www.genomics.liv.ac.uk/animal/RESEARCH/MICROSAT.PDF [accessed 10 December 2012]. Google Scholar
- R. S. Dodd , Z. Afzal-Rafii , N. Kahani , and J. Budrick . 2002. Land barriers and open oceans: Effects on gene diversity and population structure in Avicennia germinans L. (Avicenniaceae). Molecular Ecology 11: 1327–1338. Google Scholar
- J. W. Kadereit , R. Arafeh , G. Somogyi , and E. Westberg . 2005. Terrestrial growth and marine dispersal? Comparative phylogeography of five coastal plant species at a European scale. Taxon 54: 861–876. Google Scholar
- C. Lian , M. A. Wadud , Q. F. Geng , K. Shimatani , and T. Hogetsu . 2006. An improved technique for isolating codominant compound microsatellite markers. Journal of Plant Research 119: 415–417. Google Scholar
- M. Raymond , and F. Rousset . 1995. GENEPOP (version 3.3): Population genetics software for exact tests and ecumenicism. Journal of Heredity 86: 248–249. Google Scholar
- W. Rice 1989. Analyzing tables of statistical tests. Evolution; International Journal of Organic Evolution 43: 223–225. Google Scholar
- J. Sauer 1964. Revision of Canavalia. Brittonia 16: 106–181. Google Scholar
- K. Takayama , Y. Tateishi , Y. J. Murata , and T. Kjita . 2008. Gene flow and population subdivision in a pantropical plant with sea-drifted seeds Hibiscus tiliaceus and its allied species: Evidence from microsatellite analyses. Molecular Ecology 17: 2730–2742. Google Scholar
- A. Untergasser, H. Niiveen, X. Rao, T. Bisseling, R. Geurts, and J. A. M. Leunissen . 2007. Primer3Plus, an enhanced web interface to Primer3. Nucleic Acids Research 35: W71–W74. Google Scholar
- E. Westberg , and J. W. Kadereit . 2009. The influence of sea currents, past disruption of gene flow and species biology on the phylogeographical structure of coastal flowering plants. Journal of Biogeography 36: 1398–1410. Google Scholar