Three species of Selliera Cav. (Goodeniaceae), a small genus of rhizomatous perennial herbs, are currently recognized. The most common, S. radicans Cav., was originally described from Australia and also occurs in Chile and New Zealand. Within New Zealand, the species is relatively common along much of the coast and occurs less frequently in inland freshwater habitats. Selliera radicans was described by Allan (1961) as “polymorphic, with a considerable range of leaf form and size.” Ensuing investigation (Ogden, 1974) into the polymorphic nature of S. radicans distinguished an estuarine form and a dune form based on differences in rhizome formation, growth form, and leaf shape, the last of which is the most conspicuous. The estuarine form is typical of S. radicans, having elongated spatulate leaves, whereas the dune form has shorter rotund leaves. Ogden (1974) performed common garden studies and determined that the leaf form difference is strongly genetically determined and, on this basis, suggested that the two forms be regarded as distinct ecotypes. Based on further taxonomic investigations, Heenan (1997) later raised the dune ecotype to a distinct species, S. rotundifolia Heenan; however, obvious hybrid swarms between S. radicans and S. rotundifolia have been observed at sites of sympatry.
The third species in the genus, S. microphylla Colenso, was described from two regions of New Zealand in 1890 (Colenso, 1890). This species is morphologically similar to but distinguished from S. radicans primarily by a smaller form; however, this difference appears to be a plastic developmental response to the environment, as it disappears when field-collected individuals are grown in a common greenhouse (Symonds and Pilkington, pers. obs.). Based on a single individual, S. microphylla has a distinct chromosome number (2n = 56; Murray and de Lange, 2013) relative to S. radicans and S. rotundifolia (both 2n = 16; Dawson, 2000).
Given various degrees of morphological overlap, developmental plasticity, sympatry, and hybridization, genetic distinction among these species and, therefore, taxonomic status warrant further investigation. Here, using 454 pyrosequencing, microsatellite markers were developed for Selliera species for use in assessing genetic structure within and investigating hybridization among Selliera species in future work.
METHODS AND RESULTS
DNA from S. radicans collected from a population at Moana Roa beach (Appendix 1) was chosen for 454 sequencing. Genomic DNA was extracted from silica gel–dried leaf tissue using a modified cetyltrimethylammonium bromide (CTAB) method with an initial sucrose-Tris-EDTA (STE) wash (Shepherd and McLay, 2011) and an additional RNase step. The resulting DNA was dissolved in 100 µL of TE buffer. The sample had a concentration of 84.8 ng/µL and a 260/280 absorbance reading of 1.99 as measured on a NanoDrop 2000 (Thermo Fisher Scientific, Waltham, Massachusetts, USA). The DNA was run on a 1% agarose gel to assess DNA quality and ensure that RNA had been removed successfully. Approximately 5 µg of this DNA was used to construct a shotgun genomic DNA library that was sequenced in a full run on a 454 GS FLX system (454 Life Sciences, a Roche Company, Branford, Connecticut, USA) by New Zealand Genomics Ltd.
The 454 sequencing run generated more than 23 Mb of quality data, with 57,561 sequences averaging 407 bp in length. The sequence data were assembled into contigs in Geneious (version 5.6.7; Kearse et al., 2012) to increase the efficiency of microsatellite detection and to prevent locus duplication. The assembly yielded 8101 contigs with an average sequence length of 672 bp. MSATCOMMANDER version 0.8.2 (Faircloth, 2008) was used to search the Selliera contigs for di-, tri-, and tetranucleotide repeat motifs with a minimum of seven uninterrupted repeats and with the requirement to design primers at least 50 bp from the repeat region using Primer3 (Rozen and Skaletsky, 1999). Criteria for primer pair design included PCR product size between 150–350 bp with no long repeats (>4 bp) in the region surrounding the microsatellite (e.g., mononucleotide repeats) and primers optimally with 60% GC content and a GC clamp at the 3′ end.
Table 1.
Characteristics of 15 microsatellite marker primer pairs developed from Selliera radicans.
![t01_01.gif](ContentImages/Journals/apps/4/6/apps.1600012/graphic/WebImages/t01_01.gif)
From the contigs, MSATCOMMANDER (Faircloth, 2008) detected 227 repeat motifs; of these, 196 were dinucleotide repeats (86%), of which there was a high frequency of AT repeats (51%). 30 were trinucleotide repeats (13%), and one was a tetranucleotide repeat (0.4%). Given our criteria. Primer3 (Rozen and Skaletsky, 1999) successfully designed primers for 107 of the 227 repeat regions detected; 90 were designed for dinucleotide repeats. 17 were designed for trinucleotide repeats, and no primer pairs could be designed for the one tetranucleotide repeat.
Table 2.
Results of primer screening in populations of Selliera radicans.a
![t02_01.gif](ContentImages/Journals/apps/4/6/apps.1600012/graphic/WebImages/t02_01.gif)
From the 107 primer pairs, 43 were selected for initial testing based on a refinement of criteria, including a maximum number of uninterrupted repeats (12), primer melting temperatures, and overall maximum repeat length. Selected primer pairs were manufactured by IDT (Coralville, Iowa, USA) and screened initially on 15 individuals representing multiple populations of S. radicans, S. rotundifolia, and S. microphylla. PCR amplification was performed in a volume of 10 µL with 1× buffer BD (Solis BioDyne, Tartu, Estonia), 50 µM of each dNTP, 2.5 µM MgCl2, 0.5 units of FIREPol DNA polymerase (Solis BioDyne), 20 nM of forward primer. 450 nM of reverse primer, and 450 nM M13 tail primer labeled with FAM (see Schuelke, 2000 for M13-tailed PCR). Amplification by PCR was attained by: 95°C for 3 min; 35 cycles of 95°C for 30 s. 53°C for 30 s, and 72°C for 1 min; and a final 72°C hold step for 20 min. PCR products were separated on an ABI 3730 DNA Analyzer (Applied Biosystems, Waltham, Massachusetts, USA) at the Massey Genome Service (Palmerston North, New Zealand) and analyzed in GeneMapper (version 4.0; Applied Biosystems) using CASS size standard (Symonds and Lloyd, 2004). Fifteen loci amplified consistently and generated easily interpretable results (Table 1), 10 of which were polymorphic. The 10 polymorphic loci were further tested across six populations (Table 2) using the methods above, except that for each sample, three markers were pooled together, each labeled with either FAM, NED, or VIC dyes. Twenty individuals were sampled from each population from New Zealand representing all three species, and 10 individuals were sampled from one population of S. radicans from Australia (see Appendix 1).
Table 3.
Results of primer screening in populations of Selliera microphylla and S. rotundifolia.a
![t03_01.gif](ContentImages/Journals/apps/4/6/apps.1600012/graphic/WebImages/t03_01.gif)
The observed and expected heterozygosity and the number of alleles per locus were calculated using GenAlEx 6.5 (Peakall and Smouse, 2012). All 10 polymorphic markers amplified consistently and yielded two to seven alleles per locus in S. radicans. Across S. radicans populations and loci, observed heterozygosity ranged from 0 to 1.00, averaging 0.23 (SE = 0.07), and expected heterozygosity ranged from 0 to 0.66, averaging 0.47 (SE = 0.08) (Table 2). The single population of S. rotundifolia had observed heterozygosities ranging by locus from 0 to 0.80, averaging 0.23 (SE = 0.09), and expected heterozygosities of 0 to 0.61, averaging 0.26 (SE = 0.08) (Table 3). The S. microphylla population had observed heterozygosities ranging by locus from 0 to 0.45, averaging 0.10 (SE = 0.05), and expected heterozygosities of 0 to 0.59, averaging 0.17 (SE = 0.06) (Table 3).
CONCLUSIONS
Here we describe the development of 15 new markers from microsatellite loci isolated from S. radicans. Ten of the new markers are polymorphic within S. radicans and also amplify from the New Zealand congeners S. rotundifolia and S. microphylla. These 10 markers will be used in future studies of population structure and hybridization in the genus.
ACKNOWLEDGMENTS
The authors thank the Institute of Fundamental Sciences at Massey University for funding.