- Short report
- Open Access
The characterisation of microsatellite markers reveals tetraploidy in the Greater Water Parsnip, Sium latifolium (Apiaceae)
BMC Research Notesvolume 10, Article number: 204 (2017)
The Greater Water Parsnip, Sium latifolium (Apiaceae), is a marginal aquatic perennial currently endangered in England and consequently the focus of a number of conservation translocation projects. Microsatellite markers were developed for S. latifolium to facilitate comparison of genetic diversity and composition between natural and introduced populations.
We selected 65 S. latifolium microsatellite (MiSeq) sequences and designed primer pairs for these. Primer sets were tested in 32 individuals. We found 15 polymorphic loci that amplified consistently. For the selected 15 loci, the number of alleles per locus ranged from 8 to 17. For all loci, S. latifolium individuals displayed up to four alleles indicating polyploidy in this species.
These are the first microsatellite loci developed for S. latifolium and each individual displayed 1–4 alleles per locus, suggesting polyploidy in this species. These markers provide a valuable resource in evaluating the population genetic composition of this endangered species and thus will be useful for guiding conservation and future translocations of the species.
Plant translocation is a common occurrence, with an estimated 600 species of plants having been relocated as population introduction, re-introduction or augmentation [1, 2]. Whilst a tactic for large scale habitat restoration is through the planting of multiple species, translocation is also an important conservation strategy for specific plants at risk . Guidance on plant translocations recommends consideration of genetic composition  however projects infrequently utilise genetic techniques in planning and evaluating reintroductions (; although see [6, 7] as examples).
One species that has been widely translocated in the UK is Sium latifolium L., the Greater Water Parsnip. S. latifolium is a herbaceous, marginal aquatic perennial in the plant family Apiaceae, tribe Oenantheae; one of nine species within the genus, it is found across Europe and Asia . With large, conspicuous, umbel inflorescences and growing to 2 m tall , S. latifolium was once a noticeable dominant in wetland areas of England, where it grows in habitats of fen, pond margins and grazing marsh ditches . However, the population of S. latifolium has much declined over the past 40 years, due to habitat loss and change in wetland management . It is now classified as ‘endangered’ on the vascular plant red list for England . As a response to the marked decline in populations, conservation projects involving translocations of S. latifolium have occurred independently in at least seven counties of England, re-introducing the species in regions where it has been lost or declined, however the success of these translocations has been mixed.
The goal of this study was to generate a suite of microsatellite markers specifically developed for S. latifolium in order to evaluate and compare the genetic composition of populations, both old and new, with the view to guide practitioners in the best approaches for further translocations of this species. With many independent reintroductions it can also be used as a case study for exploring broader questions relating to genetic management of plant translocations.
Samples of S. latifolium were collected in May 2012 and August 2013 (Table 1), permission for sampling was obtained from the landowner of each site. Three leaflets per plant were preserved in silica gel and stored at room temperature. Prior to extraction, 10–20 mg of leaf tissue was frozen overnight at −80 °C before being homogenised at 1000 Hz for 3 min using a GenoGrinder 2000 (Spex CertiPrep, Metuchen, NJ USA). Genomic DNA was isolated employing a cetyltrimethyl ammonium bromide (CTAB) protocol , with the addition of 1% polyvinyl pyrrolidone (PVP) to the isolation buffer to remove polyphenols . Once washed and air-dried, DNA was re-suspended in 100 µl low TE (10 mM Tris–HCl, 0.1 mM EDTA, pH 8.4) and subsequently diluted to 100 ng/µl with low TE.
The microsatellite library was prepared from one individual sampled at Wickhampton Marshes, Norfolk, UK (52°35′N 1°35′E; sample identification code = I50). The library was enriched for microsatellites, using magnetic beads in the hybridisation [15, 16]. An Illumina paired-end library was created using 1 µg of the repeat-enriched genomic DNA. The SureSelect Library Prep Kit, ILM (Agilent Technologies Inc. Santa Clara, California) protocol was followed and 2 × 250 bp paired-end sequencing conducted using a MiSeq Benchtop Sequencer (Illumina Inc. San Diego, California).
Sequences with at least ten repeats were selected for primer design; primer sets were designed to amplify the microsatellite regions using PRIMER3 v 0.4.0 . Specifications for primer selection were set at a primer length of 16–36 base pairs (optimum 20 bp), an optimal primer melting temperature of 60 °C, (min–max of 59–61 °C), a maximum of 0.5 °C between primers, presence of a 3′ GC clamp, a maximum poly-X of three and the default settings for all other parameters. Sixty-five primer sets were designed. The 5′ end of each forward primer was fluorescently-labelled with HEX or 6-FAM.
Microsatellites were amplified in 2-µl PCRs, including 1 µl (100 ng) genomic DNA (air dried), 2 µl primer mix (forward and reverse primer at 0.2 µM) and 1 µl Qiagen Multiplex PCR Master Mix including HotStar Taq DNA polymerase (Qiagen Inc.). Covered with a thin layer of mineral oil, products were amplified under the following profile: incubate at 95 °C for 15 min, followed by 35 cycles of 94 °C for 30 s, selected primer temperature (51, 53 or 58 °C, see Table 2) for 1 min 30 s and 72 °C for 1 min 30 s, and finally incubated at 72 °C for 10 min. The optimum annealing temperature for each primer set was initially selected by testing a temperature gradient on two samples (Table 1), this varied the annealing temperature for each well across 12 rows from 50 to 70 °C. PCR products were diluted with double-deionized H2O (1:160). They were visualised on an ABI 3730 48-well capillary DNA Analyser (Applied Biosystems Inc. California, USA) and sized with a ROX-labelled size standard. Allele sizes were scored using GENEMAPPER v3.7 software (Applied Biosystems Inc. California, USA).
All primer sets were initially tested in six unrelated individuals (Table 1), each from a different geographic population in the UK. Markers failing to amplify or appearing monomorphic at this stage were discarded. The remaining primer sets were then tested in a further 24 individuals from the same population as the individual sequenced to isolate the microsatellites (Wickhampton Marshes, England; I50; Table 1) to fully evaluate their characteristics and usefulness. Overall, of the 65 primer pairs tested, 15 (23%) loci were polymorphic and easily scoreable (Table 2). The remainder were monomorphic (18%), not useable due to stutter and scoring difficulty (31%) or had poor/no amplification (28%).
To estimate genotyping error, extraction and scoring for a proportion of individuals was repeated to compare the data. The mean scoring error was found to be 0.02% (calculated as per ). All of the 15 markers tested displayed more than 2 alleles in multiple individuals and all individuals tested displayed more than 2 alleles in several markers, suggesting S. latifolium is polyploid (for data, see Additional file 1). A maximum of 4 alleles were observed per individual indicating tetraploidy in this species (see Additional file 2). Characteristics of each microsatellite locus were calculated for S. latifolium samples using the R package polysat [19, 20]. The number of alleles per locus ranged from 8 to 17 and the mean average number of alleles was 12 (Table 3). Observed heterozygosity per locus ranged from 0.88 to 1.00, with a mean average of 0.99 (Table 3). Due to polyploidy and unknown inheritance patterns, deviation from Hardy–Weinberg equilibrium could not be calculated nor could the frequency of null alleles be estimated .
Initial measures of genetic diversity were calculated for the genotyped population (Wickhampton Marshes) using the programme GenoDive . In this population, the mean average number of alleles per locus was 9.13 and observed heterozygosity was 0.976. Genetic distances between individuals within the library population were calculated (Bruvo distance, R package polysat [20, 23]) and visualised by ordination (R package Vegan ). The microsatellite markers revealed variation in the genetic distance between individuals within a single population and identified clusters of individuals with similar genotypes (Fig. 1).
We have successfully developed the first set of microsatellite markers for S. latifolium. The 15 loci amplified reliably and have been shown to be sufficiently variable for distinguishing between individuals (Fig. 1). These will be helpful in providing a genetic context for planning and managing further reintroductions of S. latifolium. Additionally, using S. latifolium as an example species, these microsatellite loci will also be helpful in interpreting the effects of genetic diversity and source population composition on plant reintroductions.
We also found each S. latifolium individual genotyped displayed 1–4 alleles. We conclude that this is evidence of tetraploidy, a trait not previously reported in this species. Polyploidy occurs occasionally through the Apiaceae family, in just over 10% of species . In other species of Sium intraspecific variation in ploidy levels has been recorded, with local polyploid cytotypes found within a diploid species . A chromosome count of 12 or 20 has been reported in S. latifolium . As these previous cytological studies used specimens from continental Europe, the chromosomal characteristics of UK S. latifolium is unknown. Differences in records suggests that there may be variation within the species and all reported counts are a multiple of 4, indicating that tetraploidy is possible. Additional cytological analyses would also consider historical polyploidy or aneuploidy as causes of the multiple alleles observed. Further work on S. latifolium is needed to determine the nature of the ploidy (i.e. the inheritance type) and the patterns of ploidy throughout the species’ geographic range.
cetyltrimethyl ammonium bromide
ethylenediamine tetraacetic acid
polymerase chain reaction
Dalrymple SE, Stewart GB, Pullin AS. Are reintroductions an effective way of mitigating against plant extinctions? CEE review 07-008 (SR32). Collaboration for environmental evidence. 2011. http://www.environmentalevidence.org/SR32.html.
Godefroid S, Vanderborght T. Plant reintroductions: the need for a global database. Biodivers Conserv. 2011;20:3683–8.
Maunder M. Plant reintroduction: an overview. Biodivers Conserv. 1992;1:51–61.
IUCN/SSC. Guidelines for reintroductions and other conservation translocations. Gland: IUCN Species Survival Commission; 2013.
Weeks AR, Sgro CM, Young AG, Frankham R, Mitchell NJ, Miller KA, Byrne M, Coates DJ, Eldridge MDB, Sunnucks P, Breed MF, James EA, Hoffmann AA. Assessing the benefits and risks of translocations in changing environments: a genetic perspective. Evol Appl. 2011;4(6):709–25.
Gonzalez-Perez MA, Lledo MD, Lexer C, Fay M, Marrero M, Banares-Baudet A, Carque E, Sosa PA. Genetic diversity and differentiation in natural and reintroduced populations of Bencomia exstipulata and comparisons with B. caudata (Rosaceae) in the Canary Islands: an analysis using microsatellites. Bot J Linn Soc. 2009;160:429–41.
Lloyd MW, Burnett RK, Engelhardt KAM, Neel MC. Does genetic diversity of restored sites differ from natural sites? A comparison of Vallisneria americana (Hydrocharitaceae) populations within the Chesapeake Bay. Conserv Genet. 2012;13:753–65.
Spalik K, Downie SR. The evolutionary history of Sium sensu lato (Apiaceae): dispersal, vicariance, and domestication as inferred from ITS rDNA phylogeny. Am J Bot. 2006;93:747–61.
Stace C. New flora of the British Isles. 3rd ed. Cambridge: Cambridge University Press; 2010.
Preston C, Pearman D, Dines T. New atlas of the British and Irish flora. Oxford: Oxford University Press; 2002.
Stewart ADP, Preston C. Scarce plants in Britain. Peterborough: Joint Nature Conservation Committee; 1994.
Stroh P, Leach S, August T, Walker K, Pearman D, Rumsey F, Harrower C, Fay M, Martin J, Pankhurst T, Preston C, Taylor I. A vascular plant red list for England. Bristol: Botanical Society of Britain and Ireland; 2014.
Doyle JJ, Doyle JL. A rapid DNA isolation procedure for small quantities of fresh leaf tissue. Phytochem Bull. 1987;19:11–5.
Maliyakal EJ. An efficient method for isolation of RNA and DNA from plants containing polyphenolics. Nucleic Acids Res. 1992;20:2381.
Armour JAL, Neumann R, Gobert S, Jeffreys AJ. Isolation of human simple repeat loci by hybridization selection. Hum Mol Genet. 1994;3:599–605.
Glenn TC, Schable NA. Isolating microsatellite DNA loci. Method Enzymol. 2005;395:202–22.
Rozen S, Skaletsky HJ. Primer 3 on the WWW for general users and for biologist programers. In: Misener S, Krawetz S, editors. Methods in molecular biology: bioinformatics methods and protocols. Totawa: Humana Press Inc.; 1999. p. 1999.
Hoffman JI, Amos W. Microsatellite genotyping errors: detection approaches, common sources and consequences for paternal exclusion. Mol Ecol. 2005;14:599–612.
R Development Core Team R. R. A language and environment for statistical computing. R foundation for statistical. Computing. 2011;1:409.
Clark LV, Jasieniuk M. Polysat: an R package for polyploid microsatellite analysis. Mol Ecol Resour. 2011;11(3):562–6.
Dufresne F, Stift M, Vergilino R, Mable BK. Recent progress and challenges in population genetics of polyploid organisms: an overview of current state-of-the-art molecular and statistical tools. Mol Ecol. 2014;23:40–69.
Meirmans PG, Van Tienderen PH. GENOTYPE and GENODIVE: two programs for the analysis of genetic diversity of asexual organisms. Mol Ecol Notes. 2004;4:792–4.
Bruvo R, Michiels NK, D’Souza TG, Schulenburg H. A simple method for calculation of microsatellite genotype distances irrespective of ploidy level. Mol Ecol. 2004;13:2101–6.
Oksanen J, Blanchet FG, Kindt R, Legendre P, O’Hara RB, Simpson GL, Solymos P, Stevens MHH, Wagner H. Vegan: community ecology package, R package version 2.0–4. 2012. http://CRAN.R-project.org/package=vegan. Accessed 7 Apr 2016.
Wood TE, Takebayashi N, Barker MS, Mayrose I, Greenspoon PB, Rieseberg LH. The frequency of polyploid speciation in vascular plants. P Natl Acad Sci USA. 2009;106(33):13875–9.
Kumar S, Jeelani SM, Rani S, Kumari S, Gupta RC. Cytological evaluation of Apiaceae Lindl from Western Himalayas. Cytol Genet. 2014;48(4):250–7.
Rice A, Glick L, Abadi S, Einhorn M, Kopelman N, Salman-Minkov A, Mayzel J, Chay O, Mayrose I. The Chromosome Counts Database (CCDB)—a community resource of plant chormosome numbers. New Phytol. 2015;206:19–26.
ND performed DNA extraction, designed primers, optimised PCR reactions, selected and validated markers, conducted data analysis, and drafted the manuscript. GH constructed the microsatellite-enriched genomic library, designed primers and helped in optimising PCR reactions and interpreting results. DD participated in primer design, marker selection and validation, discussion of results, and revised all drafts of the manuscript. All authors read and approved the final manuscript.
Microsatellite isolation and genotyping was performed at the NERC Biomolecular Analysis Facility at the University of Sheffield, UK. The authors thank J. Dawe and D. Grafham of Sheffield Diagnostics Genetics Service for performing the MiSeq sequencing at The Children’s Hospital Sheffield supported by the Sheffield Children’s NHS Trust, UK.
The authors declare that they have no competing interests.
Availability of data and materials
The microsatellite sequences are available through the European Molecular Biology Laboratory (EMBL/EBI) European Nucleotide Archive (see http://www.ebi.ac.uk/ena); ENA Accession Numbers LN849725 to LN849739. The data generated and analysed during the study (sample genotypes) are included in the Additional files of this report.
This research was supported by the Natural Environment Research Council (NERC), UK (Award NBAF769). ND was supported by NERC Grant NE/1018336/1 and the Somerset Wildlife Trust, UK.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.