Labisia pumila (Blume) Fern.-Vill. (Myrsinaceae) is a small understory shrub that is widely distributed in the tropical forests of Malaysia, Indonesia, Thailand, the Philippines, and Myanmar (Sunarno, 2005). Eight varieties are recognized (Sunarno, 2005), of which L. pumila var. pumila, L. pumila var. alata (Scheff.) Mez, and L. pumila var. lanceolata (Scheff.) Mez are commonly found in Malaysia. These varieties are morphologically distinct from one another in terms of their petiole and leaf characteristics. Among the Malay communities, these varieties are collectively known as kacip fatimah, which has long been used as the traditional medicine for the treatment of pre- and postpartum complications, menstrual disorders, dysentery, rheumatism, flatulence, and gonnorrhea (Burkill, 1966; Jaganath, 2000).
To date, the exploitation of L. pumila for commercial demand, particularly in the pharmacological and cosmeceutical applications, is gradually increasing. Attention to conservation should therefore be prioritized to ensure sustainable utilization. Despite the importance of L. pumila, the availability of genetic information for the species is still very limited. Only two genetic variability studies of the species have been reported using dominant markers (Bhore et al, 2009; Farah Fazwa et al, 2013). In this study, we report the development of 10 microsatellite loci in L. pumila var. alata, and we evaluate their polymorphism across var. alata, var. pumila, and var. lanceolata. Microsatellites are preferred markers because of the nature of their codominant inheritance, high abundance, extent of allelic diversity, and the ease of assessing the size variation by PCR with pairs of flanking primers (Weising et al, 2005).
METHODS AND RESULTS
Leaf samples of 25, 20, and four individuals of L. pumila var. alata, var. pumila, and var. lanceolata, respectively, were obtained from Pasoh Forest Reserve (2°58′N, 102°18′E). An additional six individuals of var. lanceolata were obtained from the Ethnobotanical Garden of the Forest Research Institute Malaysia (FRIM). The voucher specimens of these three varieties were deposited in FRIM Herbarium (KEP; barcode numbers 223663-223665). Total genomic DNA was extracted from fresh leaves of L. pumila using a modified cetyltrimethylammonium bromide (CTAB) protocol (Murray and Thompson, 1980) and further purified using the High Pure PCR Template Preparation Kit (Roche Diagnostics GmbH, Penzberg, Germany).
A genomic library enriched for dinucleotide CT and GT was constructed following the approach of Kijas et al. (1994). Approximately 5 µg of genomic DNA was obtained from an individual of L. pumila var. alata from the Ethnobotanical Garden of FRIM. After digestion with NdeII (Promega Corporation, Madison, Wisconsin, USA), the digested fragments were electrophoresed on 2% agarose gels with a 100-bp DNA ladder (New England Biolabs, Ipswich, Massachusetts, USA). Fragments of 300–1000 bp were excised and ligated into Sau3A1 cassettes (TaKaRa Bio, Otsu, Shiga, Japan). After ligation, the nicks were repaired using DNA polymerase I (TaKaRa Bio). The cassette-ligated DNA was enriched for microsatellite repeats via hybridization to 5′-biotinylated (CT)15 and (GT)15 probes and retrieved using magnetic beads coated with streptavidin (Promega Corporation). The selectively recovered hybrids were reamplified using C1 cassette primers, digested with NdeII, cloned into pUC118 Bam H1/BAP vector (TaKaRa Bio), and transformed into QIAGEN EZ Competent Cells (QIAGEN Gmbh, Hilden, Germany). Insert-containing clones were selected by blue/white screening on Luria-Bertani (LB) agar plates containing 100 µg/mL ampicillin, 50 µM isopropyl-β-d-1 -thiogalactopyranoside (IPTG), and 80 µg/mL X-gal. Plasmid DNAs of a total of 608 clones were amplified using the Illustra TempliPhi Amplification Kit (GE Healthcare, Piscataway, New Jersey, USA) and sequenced using BigDye Terminator Sequencing Kit version 3.1 (Applied Biosystems, Foster City, California, USA) on an ABI 3130×1 Genetic Analyzer (Applied Biosystems).
Description of 10 polymorphic and three monomorphic microsatellites screened in Labisia pumila var. alata.
Sequence redundancy was checked using CD-HIT (Li and Godzik, 2006), and microsatellite-containing sequences were then identified from the nonredundant sequences using the MIcroSAtellite Identification Tool (MISA) (Thiel et al., 2003). A total of 56 sequences harbored at least one microsatellite motif. Forty primer pairs of 19–24 nucleotides in length were designed using OLIGO 6.67 (Molecular Biology Insights, Cascade, Colorado, USA), based on the criteria of GC content (40–60%), melting temperature (52–58°C), and the lack of secondary structures. Four unrelated individuals of L. pumila var. alata were used for initial screening of these primer pairs. The PCR amplifications were carried out in 10-µL volumes containing 10 ng of template DNA, 1× GoTaq Flexi Buffer, 1.5 mM MgCl2, 0.3 µM of each primer, 0.2 mM of each dNTP, and 0.5 unit of GoTaq Flexi DNA polymerase (Promega Corporation). The reaction mixture was amplified using a GeneAmp PCR System 9700 (Applied Biosystems) with the following cycling conditions: 4 min at 94°C; 35 cycles of 94°C for 1 min, 50°C for 30 s, and 72°C for 30 s; followed by a final extension of 30 min at 72°C. Twenty-nine primer pairs that showed specific amplification products of expected fragment size were selected for fluorescent labeling at the 5′-end of the forward primers with either 6-FAM or HEX. These primers were further screened using 25 samples of L. pumila var. alata using the same cycling conditions as above. Fragment analysis was conducted using an ABI 3130×1 with GeneScan 400 ROX (Applied Biosystems) as the internal size standard. Allele sizes were assigned using GeneMapper version 4.0 (Applied Biosystems). Characteristics of the microsatellites developed for L. pumila var. alata are shown in Table 1. Thirteen loci yielded consistent and scorable genotypes; of these, 10 were polymorphic and the remaining three were monomorphic. Because the majority of the loci evaluated (Lpu02, Lpu08a, Lpu08b, Lpu13, Lpu16b, Lpu21a, and Lpu21b) exhibited up to four alleles per individual (Fig. 1), our results suggest that the species is likely a tetraploid.
To assess the degree of polymorphism across L. pumila var. alata, var. pumila, and var. lanceolata, the 10 polymorphic microsatellites were further checked for variability in 55 individuals of the three varieties using the same PCR and genotyping protocols as above. Observed heterozygosity (Ho) of each variety was calculated based on the allele phenotypes, following Bever and Felber (1992), with values being one minus the probability that any two alleles drawn at random were identical by descent (AAAA = 0, AAAB = 0.50, AABB = 0.667, AABC = 0.833, and ABCD = 1). Minimum and maximum bounds of Ho were calculated for biallelic phenotypes (AAAB and AABB) (De Walt et al., 2011). Expected heterozygosity (He) was calculated using ATETRA 1.2 (Van Puyvelde et al., 2010), with 10,000 Monte Carlo simulations. A total of 84, 48, and 66 alleles were observed in L. pumila var. alata, var. pumila, and var. lanceolata, respectively (Table 2). He ranged from 0.039 to 0.857, 0.000 to 0.793, and 0.000 to 0.874 in L. pumila var. alata, var. pumila, and var. lanceolata, respectively. Notably, locus Lpu15 was found to be monomorphic in varieties pumila and lanceolata.
Genetic properties of 10 microsatellites of Labisia pumila across varieties alata, pumila, and lanceolata.
This is the first report on the development of microsatellites in L. pumila. The observed levels of polymorphism and genetic diversity suggest that, apart from monomorphic loci, these microsatellites can serve as useful tools for DNA profiling and population genetic studies of the species.
- J. D. Bever , and F. Felber . 1992. The theoretical population genetics of autopolyploidy. In J. Antonovics and D. J. Futuyama [eds.], Oxford surveys in evolutionary biology, 185–217. Oxford University Press, Oxford, United Kingdom. Google Scholar
- S. J. Bhore , A. H. Nurul , and F. H. Shah . 2009. Genetic variability based on randomly amplified polymorphic DNA in Kacip Fatimah (Labisia pumila Benth & Hook f.) collected from Melaka and Negeri Sembilan states of Malaysia. Journal of Forest Science 25: 93–100. Google Scholar
- I. H. Burkill 1966. A dictionary of the economic product of the Malay Peninsula, vol. II (I–Z). Government of Malaysia and Singapore by the Ministry of Agriculture and Cooperative, Kuala Lumpur, Malaysia. Google Scholar
- S. J. De Walt , E. Siemann , and W. E. Rogers . 2011. Geographic distribution of genetic variation among native and introduced populations of Chinese tallow tree, Triadica sebifera (Euphorbiaceae). American Journal of Botany 98: 1128–1138. Google Scholar
- M. A. Farah Fazwa , H. Siti Salwana , H. Maideen , and O. Mohamad . 2013. An assessment of genetic relationship among superior accessions of Labisia pumila analyzed by amplified fragment length polymorphism (AFLP) markers. Open Science Repository Agriculture. doi: 10.7392/Agriculture.70081945 Google Scholar
- N. L. Jaganath 2000. Herbs: The green pharmacy of Malaysia. Vinpress Sdn. Bhd. in collaboration with the Malaysian Agricultural Research and Development (MARDI), Kuala Lumpur, Malaysia. Google Scholar
- J. M. Kijas , J. C. Fowler , C. A. Garbett , and M. R. Thomas . 1994. Enrichment of microsatellites from the citrus genome using biotinylated oligonucleotide sequences bound to streptavidin-coated magnetic particles. BioTechniques 16: 656–662. Google Scholar
- W. Li , and A. Godzik . 2006. CD-HIT: A fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics (Oxford, England) 22: 1658–1659. Google Scholar
- M. Murray , and W. F. Thompson . 1980. Rapid isolation of high molecular weight plant DNA. Nucleic Acids Research 8: 4321–4325. Google Scholar
- B. Sunarno 2005. Revision of the genus Labisia (Myrsinaceae). Blumea 50: 579–597. Google Scholar
- T. Thiel , W. Michalek , R. K. Varshney , and A. Graner . 2003. Exploiting EST databases for the development and characterization of genederived SSR-markers in barley (Hordeum vulgare L.). Theoretical and Applied Genetics 10: 411–422. Google Scholar
- K. Van Puyvelde , A. Van Geert , and L. Triest . 2010. ATETRA, a new software program to analyse tetraploid microsatellite data: Comparison with TETRA and TETRASAT. Molecular Ecology Resources 10: 331–334. Google Scholar
- K. Weising , H. Nybom , K. Wolff , and G. Kahl . 2005. DNA fingerprinting in plants: Principles, methods, and applications, 2nd ed. Taylor & Francis Group, Boca Raton, Florida, USA. Google Scholar
 The authors thank D. Mariam, T. Sharifah, J. Ghazali, M. Yahya, P. Ramli, B. Yasri, C. S. Suryani, and H. M. Nurul for their assistance in the laboratory and in the field. We are also grateful to A. L. Tan for her help in plant identification. This project was supported by the Ministry of Agriculture and Agro-Based Industry, Malaysia, under the ScienceFund (05-03-10-SF1031).