Translator Disclaimer
24 July 2017 Development of EST-Derived Microsatellite Markers in the Aquatic Macrophyte Ranunculus bungei (Ranunculaceae)
Author Affiliations +

Ranunculus bungei Steud. (section Batrachium DC., Ranunculaceae) is a perennial submerged plant that can proliferate vegetatively through rhizomes or sexually via selfing or outcrossing (Cook, 1966). Ranunculus bungei is widely distributed in heterogeneous environments within the temperate and alpine regions of China and is significant for studies of the adaptation to aquatic habitats in angiosperms (Chen et al., 2015). In recent years, investigations have been carried out to examine genetic variation and population structure in R. bungei with intersimple sequence repeat (ISSR) markers or chloroplast noncoding spacers (Wang et al., 2010; Chen et al., 2014), but these markers are less powerful in studies on reproductive system, hybridization, patterns of gene flow, and fine-scale population structure. The development of suitable markers can provide a better understanding of the evolutionary progress and the underlying ecological factors of R. bungei and its related species.

Microsatellite or simple sequence repeat (SSR) markers are molecular markers with many desirable genetic attributes (e.g., codominant inheritance and hypervariability), which have been used to reveal genetic patterns in a wide variety of species (Kalia et al., 2011). For clonal plants, estimation of genetic variation is often biased with markers of low discriminatory ability (Arnaud-Haond et al., 2005); therefore, genetic studies on aquatic macrophytes, which are characterized by limited sexual proliferation (Barrett et al., 1993), should be assessed using appropriate polymorphic markers. Although a large number of SSR loci for Ranunculus L. species have been identified (e.g., Noel et al., 2005; Matter et al., 2012), we found that cross-species amplification was rarely successful in R. bungei based on preliminary experiments. Therefore, we developed 22 novel EST-SSR markers from R. bungei for use in investigations of population and landscape genetics of this widely distributed submerged species.


A mixture of tissues from roots, stems, and leaves was used for the transcriptome sequencing of R. bungei, conducted by Chen et al. (2015) using the Illumina HiSeq 2000 sequencer (Illumina, San Diego, California, USA). A total of 5,312,841 clean reads of R. bungei deposited in the National Center for Biotechnology Information (NCBI) Sequence Read Archive (accession no. SRR1822529; Chen et al., 2015) were used for de novo transcriptome assembly using Trinity version 2.1.0 with default parameters (Grabherr et al., 2011). The longest sequence was chosen for transcripts with several isoforms, as identified with a perl script (available at The MIcroSAtellite identification tool (MISA) Perl script (Thiel et al., 2003) was then used to screen for microsatellite motifs from total unigenes, and the minimum number of each type of repeat was set to six. MISA recovered a total of 9903 SSR motifs for R. bungei, and 50 unigenes were randomly chosen for the EST-SSR development.

PCR primer design for the targeted unigenes was performed with Primer Premier 5.0 (PREMIER Biosoft International, Palo Alto, California, USA). An initial evaluation of the primers was facilitated in 20 individuals randomly selected from four Chinese populations of R. bungei (Appendix 1), and 22 primer pairs showing unique products ranging from 100–500 bp were individually labeled with the fluorescent dyes 6-FAM or HEX (Table 1). Characterization of the EST-SSR loci was estimated in the four populations of R. bungei, with 22, 20, 15, and 16 individuals, respectively (Appendix 1). Genomic DNA was extracted from the freeze-dried leaves of R. bungei individuals using the DNAsecure Plant Kit (Tiangen Biotech, Beijing, China). PCR amplifications were performed in 20-µL reaction mixtures containing 1.5 µL of genomic DNA (∼30 ng/µL), 0.5 µL of each primer (10 µM), and 10 µL 2× Master PCR Mix (Tiangen Biotech). Microsatellites were amplified under the following PCR conditions: a 5-min initial denaturation step at 95°C; followed by 30–35 cycles of 30 s at 95°C, 30 s at 50–58°C (Table 1), and 1 min at 72°C; and a final extension at 72°C for 7 min. PCR products differed in fluorescent label or length (>80 bp) and were multiplexed and analyzed on the ABI 3730XL sequencer (Applied Biosystems, Foster City, California, USA) with GeneScan 500 LIZ Size Standard (Applied Biosystems). Microsatellite genotyping was performed using GeneMarker version 1.5 software (SoftGenetics, State College, Pennsylvania, USA). The number of alleles, observed and expected heterozygosities, and deviations from Hardy–Weinberg equilibrium (HWE) at each locus were estimated using GenAlEx 6.5 (Peakall and Smouse, 2012). Linkage disequilibrium of locus pairs was tested using Arlequin version (Excoffier et al., 2005).

Table 1.

Characteristics of 22 EST-SSR markers developed in Ranunculus bungei.


Cross-species amplification was conducted in two other species of Ranunculus section Batrachium (R. aquatilis L. var. eradicatus Laest. [n = 11] and R. trichophyllus Chaix ex Vill. [n = 14]), as well as six riparian and aquatic species of Ranunculaceae (R. cheirophyllus Hayata [n = 5], R. natans C. A. Mey. [n = 8], Halerpestes tricuspis (Maxim.) Hand.-Mazz. [n = 12], H. ruthenica (Jacq.) Ovcz. [n = 6], Caltha palustris L. [n = 3], and C. natans Pall. [n = 5]) (Appendix 1).

The characteristics of 22 EST-SSR loci are presented in Table 1. Fourteen loci were polymorphic, two of which were fixed for different alleles in multiple populations (Table 2). The loci BatrB6, BatrB9, and BatrB12 showed the highest number of alleles (five), and eight loci were monomorphic among all individuals (Table 2). The expected and observed heterozygosity ranged from 0.0 to 0.5 and 0.0 to 1.0 per locus, and all polymorphic loci showed significant deviation from HWE (Table 2). Significant linkage disequilibrium (P < 0.05) was observed among seven locus pairs in four R. bungei populations (BatrB5 and BatrB6, BatrB6 and BatrB10, BatrB6 and BatrB12, BatrB9 and BatrB15, BatrB10 and BatrB13, BatrB10 and BatrB15, and BatrB13 and BatrB15). The deviation from HWE and significant linkage disequilibrium could be explained by the small population size, inbreeding, and clonal reproduction in R. bungei. All of the loci were amplified successfully in two Ranunculus species from section Batrachium, and four loci could be amplified successfully in all eight species (Table 3). The allele size ranges were similar among Ranunculus section Batrachium species, but greatly differed among the taxa from different genera in loci BatrB1, BatrB9, BatrB11–13, BatrB17, and BatrB21.

Table 2.

Results of initial primer screening in four populations of Ranunculus bungei.a


Table 3.

Cross-amplification of 22 EST-SSR markers developed in Ranunculus bungei across eight other species of Ranunculaceae. The number of alleles in populations of each species is presented for the loci that could be successfully amplified.a



Fourteen polymorphic and eight monomorphic microsatellite loci were developed in R. bungei. The polymorphism observed for the SSRs in R. bungei is moderate when compared with other aquatic plants (Nies and Reusch, 2004; Wu et al., 2013). Cross-species amplification also indicates that these markers may be widely used in related Ranunculaceae species. We conclude that the EST-SSRs described here will facilitate ecological and evolutionary studies of R. bungei as well as related species.


The authors thank Zuyu Chen and Chun Gong for help in de novo transcriptome assembly, and Xinwei Xu and Elizabeth Schultz for plant materials and assistance with language revision. This study was supported by the State Key Laboratory of Freshwater Ecology and Biotechnology (grant no. 2016FB04).



Arnaud-Haond, S., F. Alberto, S. Teixeira, G. Procaccini, E. A. Serrão, and C. M. Duarte. 2005. Assessing genetic diversity in clonal organisms: Low diversity or low resolution? Combining power and cost efficiency in selecting markers. Journal of Heredity 96: 434–440. Google Scholar


Barrett, S. C. H., C. G. Eckert, and R. C. Husband. 1993. Evolutionary processes in aquatic plant populations. Aquatic Botany 44: 105–145. Google Scholar


Chen, J.-M., Z.-Y. Du, Y-Y. Yuan, and Q.-F. Wang. 2014. Phylogeography of an alpine aquatic herb Ranunculus bungei (Ranunculaceae) on the Qinghai–Tibet Plateau. Journal of Systematics and Evolution 52: 313–325. Google Scholar


Chen, L.-Y., S.-Y. Zhao, Q.-F. Wang, and M. L. Moody. 2015. Transcriptome sequencing of three Ranunculus species (Ranunculaceae) reveals candidate genes in adaptation from terrestrial to aquatic habitats. Scientific Reports 5: 10098. Google Scholar


Cook, C. D. K. 1966. A monographic study of Ranunculus subgenus Ranunculus (DC.) A. Gray. Mitteilungen der Botanischen Staatssammlung München 6: 47–237. Google Scholar


Excoffier, L., G. Laval, and S. Schneider. 2005. Arlequin (version 3.0): An integrated software package for population genetics data analysis. Evolutionary Bioinformatics Online 1: 47–50. Google Scholar


Grabherr, M. G., R. J. Haas, M. Yassour, J. Z. Levin, D. A. Thompson, I. Amit, X. Adiconis, et al. 2011. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nature Biotechnology 29: 644–652. Google Scholar


Kalia, R., M. Rai, S. Kalia, R. Singh, and A. Dhawan. 2011. Microsatellite markers: An overview of the recent progress in plants. Euphytica 177: 309–334. Google Scholar


Matter, P., A. R. Pluess, J. Ghazoul, and C. J. Kettle. 2012. Eight microsatellite markers for the bulbous buttercup Ranunculus bulbosus (Ranunculaceae). American Journal of Botany 99: e399–e401. Google Scholar


Nies, G., and T. B. H. Reusch. 2004. Nine polymorphic microsatellite loci for the fennel Pondweed Potamogeton pectinatus L. Molecular Ecology Notes 4: 563–565. Google Scholar


Noel, F., M.-C. Boisselier-Dubayle, J. Lambourdiere, N. Machon, J. Moret, and S. Samadi. 2005. Characterization of seven polymorphic microsatellites for the study of two Ranunculaceae: Ranunculus nodiflorus L., a rare endangered species and Ranunculus flammula L., a common closely related species. Molecular Ecology Notes 5: 827–829. Google Scholar


Peakall, R., and P. E. Smouse. 2012. GenAlEx 6.5: Genetic analysis in Excel. Population genetic software for teaching and research—an update. Bioinformatics (Oxford, England) 28: 2537–2539. Google Scholar


Thiel, T., W. Michalek, R. Varshney, and A. Graner. 2003. Exploiting EST databases for the development and characterization of gene-derived SSR-markers in barley (Hordeum vulgare L.). Theoretical and Applied Genetics 106: 411–422. Google Scholar


Wang, Y., J. Chen, C. Xu, X. Liu, Q. Wang, and T. J. Motley. 2010. Population genetic structure of an aquatic herb Ranunculus bungei (Ranuculaceae) [sic] in the Hengduan Mountains of China. Aquatic Botany 92: 221–225. Google Scholar


Wu, Z., D. Yu, and X. Xu. 2013. Development of microsatellite markers in the hexaploid aquatic macrophyte, Myriophyllum spicatum (Haloragaceae). Applications in Plant Sciences 1: 1200230. Google Scholar


Appendix 1.

Geographic information of Ranunculus, Halerpestes, and Caltha populations in this study.a

Zhigang Wu, Jinwei Wu, Yalin Wang, and Hongwei Hou "Development of EST-Derived Microsatellite Markers in the Aquatic Macrophyte Ranunculus bungei (Ranunculaceae)," Applications in Plant Sciences 5(7), (24 July 2017).
Received: 17 March 2017; Accepted: 1 April 2017; Published: 24 July 2017

Back to Top