Development and characterization of 16 polymorphic microsatellite markers from Taiwan cow-tail fir, Keteleeria davidiana var. formosana (Pinaceae) and cross-species amplification in other Keteleeria taxa
BMC Research Notes volume 7, Article number: 255 (2014)
Keteleeria davidiana var. formosana (Pinaceae), Taiwan cow-tail fir, is an endangered species listed on the IUCN Red List of Threatened Species and only two populations remain, both on the Taiwan Island. Sixteen polymorphic microsatellite loci were developed in an endangered and endemic gymnosperm species, Keteleeria davidiana var. formosana, and were tested in an additional 6 taxa, K. davidiana var. calcarea, K. davidiana var. chienpeii, K. evelyniana, K. fortunei, K. fortunei var. cyclolepis, and K. pubescens, to evaluate the genetic variation available for conservation management and to reconstruct the phylogeographic patterns of this ancient lineage.
Polymorphic primer sets were developed from K. davidiana var. formosana using the modified AFLP and magnetic bead enrichment method. The number of alleles ranged from 3 to 16, with the observed heterozygosity ranging from 0.28 to 1.00. All of the loci were found to be interspecifically amplifiable.
These polymorphic and transferable loci will be potentially useful for future studies that will focus on identifying distinct evolutionary units within species and establishing the phylogeographic patterns and the process of speciation among closely related species.
For conservation management, the evaluation of genetic diversity and the identification of distinct evolutionary units in endangered taxa should allow for the development of more efficient conservation strategies[1–3]. Over the past century, human activities associated with the lumber industry and agriculture have destroyed and fragmented forest habitats at altitudes of 500–2500 m on Taiwan Island. Several gymnosperm species have decreased in population size and are now under threat, including Cycas taitungensis[4–6], Podocarpus nakaii, Amentotaxus formosana, and Keteleeria davidiana var. formosana. Extreme population declines can affect the genetic diversity of a species. The loss of genetic diversity carries serious evolutionary concerns regarding both the ability to adapt to changes on oceanic islands and longer-term interactions with other organisms. Thus, molecular markers are used to evaluate endangered species in greater depth, with the purpose of applying this genetic diversity information to the conservation and restoration of biodiversity.
Keteleeria davidiana var. formosana (Pinaceae), Taiwan cow-tail fir, is an endangered species listed on the IUCN Red List of Threatened Species and derived from K. davidiana, including var. davidiana and var. calcarea, occurred in China by geographical disjunction and leaf scars obscurely or obviously protruding on branchlets. The taxa has only two wild populations, the Dawu Taiwan Keteleeria Forest Reserve (Dawu) and the Pinglin Taiwan Keteleeria Nature Reserve (Pinglin), which are located in Northern and Southern Taiwan, respectively. A previous inventory study has found that less than a hundred individuals remain in the wild, few seedlings are found in the field, and adult trees are vulnerable to lightning damage or death. Additionally, continuous logging over the past century in the low-elevation broadleaf forests has restricted K. davidiana var. formosana to narrow habitats and has decreased its population size. Therefore, the conservation strategies for Keteleeria davidiana var. formosana should include in situ protection and management in conjunction with ex situ approaches such as clonal orchards. Hence, the population genetics of this species should be evaluated before in situ and ex situ conservation management begins. In this study, we developed polymorphic microsatellite loci from Taiwan cow-tail fir to estimate the population structure within the species and to identify distinct evolutionary units within the population. In addition, the transferability of the microsatellite primers developed in this study was tested in other taxa, specifically, three species and three varieties of the Keteleeria genus.
Sampling and DNA extractions
Twenty-five individuals were sampled from each of the two remaining populations of K. davidiana var. formosana from Dawu and Pinglin, Taiwan (Table 1). To test the transferability of the markers, two to three individuals of three species, K. pubescen s, K. fortune, and K. evelyniana, and three varieties, K. fortunei var. cyclolepis, var. calcarea, and K. davidiana var. chienpeii, were sampled from the field. Permission for tissue collection was obtained from the Council of Agriculture, Republic of China, authorities. All taxa of the Keteleeria genus are diploid with 24 chromosomes. The sample size, location, and storage herbarium for the voucher specimens are listed in Table 1. Total genomic DNA was extracted from silica-dried leaf powder using the Plant Genomic DNA Extraction Kit (RBC Bioscience, Taipei, Taiwan).
Isolation of microsatellite DNA loci and identification
To isolate the microsatellites efficiently, we used the modified AFLP and magnetic bead enrichment method[1, 17, 18] to select microsatellite loci. Genomic DNA from one individual of K. davidiana var. formosana was treated with the restriction enzyme Mse I (Promega, Madison, Wisconsin, USA) and separated using 1% Nusieve® 3:1 agarose gel (FMC Bio Products, Rockland, ME, USA) electrophoresis. The cut DNA fragments ranged from 400 to 1000 bps, and the purified fragments were extracted from the agarose gel using the HiYield™ Gel PCR DNA Fragments Extraction Kit (RBC Bioscience). Purified DNA fragments were ligated to a double-stranded Mse I-adaptor (complementary oligo A: 5′-TACTCAGGACTCAT-3′; 5′phosphorylated oligo B: 5′-GACGATGAGTCCTGAG-3′) using the Quick Ligation™ Kit (New England Biolabs, MA, USA) at 25°C for 5 minutes. The ligation product was used as the template DNA for the enrichment of the partial genomic library, and the adaptor-specific primers, named Mse I-N (5′-GATGAGTCCTGAGTAAN-3′), were used to perform 20 cycles of prehybridization PCR amplification. The PCR mixture contained 20 ng template DNA, 10 pmol adapter-specific primer, 2 μL 10 × reaction buffer, 2 mM dNTP mix, 2 mM MgCl2, 0.5 U Taq DNA polymerase (Promega), and sterile water was added to reach a total volume of 20 μL. The PCR protocol was set at 94°C for 5 min, followed by 18 cycles of 94°C for 30 s, 53°C for 1 min, and 72°C for 1 min, using a Labnet MultiGene 96-well Gradient Thermal Cycler (Labnet, Edison, New Jersey, USA). The amplicons were denatured and hybridized to four separated 5′-biotinylated oligonucleotide probes, (AG)15, (AC)15, (TCC)10, and (TTG)10, in 250 μL of hybridization solution at 68°C for 1 h, and 1 mg of Streptavidin MagneSphere Paramagnetic Particles (Promega) was added to the mixture solution to capture the hybridizations at 42°C for 2 h. The four enriched microsatellite DNA fragments were washed with high-salt and low-salt solutions and then used as templates for 25 cycles of PCR amplification using the adaptor-specific primers; the amplification protocol was same as that used for the prehybridization PCR. The PCR products were purified using the HiYield™ Gel PCR DNA Fragments Extraction Kit (RBC Bioscience) and then cloned using the pGEM®-T Easy Vector System (Promega). The plasmid DNAs that were purified from white colonies were digested by the restriction enzyme Eco RI (Promega) and screened using 1% agarose electrophoresis to determine the size of the inserted DNA fragments. In total, 670 white colonies were selected to purify plasmids for sequencing from Keteleeria davidiana var. formosana. The sequencing cocktail contained 200 ng plasmid DNA as the template DNA, 0.5 mM T7 or SP6 primers, 4 μL Big Dye Terminator v3.1 Ready Reaction Premix Reactions (Applied Biosystems, Carlsbad, California, USA), and sequencing buffer in a total volume of 20 μL. Selected plasmid DNAs were sequenced in both directions using an ABI PRISM® 3700 DNA Sequencer (Applied Biosystems). These sequences, including their microsatellite loci, were recognized using Tandem Repeats Finder version 4.07b, and a pair of specific primers for each microsatellite locus was designed using FastPCR software version 6.4.18.
DNA amplification and genotyping
The optimal annealing temperature was evaluated by gradient PCR using a Labnet MultiGene 96-well Gradient Thermal Cycler (Labnet) for a temperature range from 50°C to 65°C. For each Keteleeria taxon, two to three individuals were evaluated for the optimal annealing temperature. The polymorphism of the two remaining populations of K. davidiana var. formosana was evaluated using 25 individuals from each population (Table 1). To test the optimal annealing temperature, a 20 μL reaction mixture containing 20 ng template DNA, 0.2 μM each of the reverse and forward primers, 2 μL 10 × reaction buffer, 2 mM dNTP mix, 2 mM MgCl2, 0.5 U Taq DNA polymerase (Promega), and sterile water was amplified using a Labnet MultiGene 96-well Gradient Thermal Cycler (Labnet). The PCR program was set at 94°C for 5 min, followed by 30 cycles of 94°C for 40 s, a temperature gradient ranging from 50 to 65°C for 60 s, and 72°C for 60 s, and a final extension of 72°C for 10 minutes. Amplicons were checked by 1% agarose electrophoresis to isolate the target DNA bands, which were confirmed by sequencing. To examine the genetic polymorphisms, the PCR cocktail was amplified using the previously described program with the temperature gradient replaced by the optimal annealing temperature (Ta). The amplicons were checked by 1% agarose electrophoresis to isolate the target DNA bands, which were confirmed by sequencing. The amplicons were separated by electrophoresis on a 10% polyacrylamide gel (acrylamide: bisacrylamide 29: 1, 80 V for 14–16 hours) using a 25 or 50 bp DNA Step Ladder (Promega) to determine the allele size. The bands were then imaged under UV light using the Flo Gel FGIS-3 fluorescent gel image system (Top BIO Co., Taipei, Taiwan), and the sizes of the PCR products were identified using Quantity One software version 4.62 (Bio-Rad Laboratories, Hercules, California, USA).
The genetic variation indices, including the number of alleles (Na), the number of effective alleles (Ne), the observed heterozygosity (Ho) and expected heterozygosity (He), Shannon’s information index (H), and fixation index (F IS ) were estimated using GenAlEx version 6.5. The Hardy–Weinberg equilibrium (H WE ) was tested using Arlequin software version 126.96.36.199.
Results and discussion
Development of polymorphic microsatellite markers
A total of 392 microsatellite loci containing repeat motifs were discovered, with maximum and minimum lengths of 992 bps and 119 bps, respectively, and an average sequence length of 477 bps. We designed 102 primer pairs between the up- and down-flanking regions of the motifs based on the primer design parameters computed using FastPCR software version 6.4.18. To test the optimal annealing temperatures, which were obtained using gradient temperature PCRs, template DNA was derived from three individuals of each of the Keteleeria taxa, which are listed in Table 1. We tested the polymorphisms for K. davidiana var. formosana using 25 individuals from each population. Finally, we selected 16 polymorphic loci from the 102 microsatellites based on the detection of unambiguous polymorphic amplicons using fixed annealing temperature PCR and a polyacrylamide gel genotyping protocol. The characteristics of the 16 polymorphic microsatellite loci are listed in Table 2. Of the 16 loci, 11 are complete microsatellite loci, consisting of 5 with a dinucleotide motif, 5 with a trinucleotide motif, and 1 with a tetranucleotide motif. Of the 5 remaining loci, 3 carried a compound motif and 2 carried an interrupted motif. The sequences of 16 loci reported in this paper are available from GenBank (accession numbers: HG518488–HG518503) (Table 2).
Genotyping and population genetics analysis
To examine the extent of genetic polymorphism at each locus, 25 individuals were collected in the field from each of the two remaining populations of K. davidiana var. formosana (Table 1). For the 16 new polymorphic microsatellite loci, the number of alleles per locus (Na) and the number of effective alleles (Ne) ranged from 5 to 14 and 2.59 to 11.25 in the Dawu population, respectively, and from 3 to 16 and 2.14 to 10.59 in the Pinglin population, respectively (Table 3). The observed and expected heterozygosity (Ho and He, respectively) varied from 0.44 to 1.00 (average of 0.68) and 0.61 to 0.91 (average of 0.82) in the Dawu population, respectively, and from 0.28 (average of 0.63) to 0.80 (average of 0.78) in the Pinglin population. The means of Shannon’s information index (H) and the fixation index (F IS ) were 1.93 and 0.175 and 1.73 and 0.188 in the Dawu and Pinglin populations, respectively (Table 3). Significant deviations from Hardy–Weinberg equilibrium (H WE ) were detected at 6 and 3 loci in the two populations, respectively, and these deviations were attributed to the heterozygote deficiency of the endangered species (Table 3). A total of 52 and 29 private alleles were observed in the Dawu and Pinglin populations, respectively, revealing population differentiation between the two remaining populations of the species.
The genetic variability, including the means of the observed and expected heterozygosity (Table 3), of K. davidiana var. formosana was high compared with that of other Pinaceae species such as Pinus koraiensis (0.38 and 0.57) and P. massoniana (0.27 and 0.65)[25, 26], but was similar to that of P. pinaster (0.65 and 0.83). Unfortunately, no data for other Keteleeria taxa are available for the comparison of genetic variability. However, the high observed and expected heterozygosity values implied that the species was historically larger in population size and was more common in low-elevation forests than it is at present.
Test of transferability
To test the transferability of these microsatellite loci, we tested these primers in six other cow-tail fir taxa: K. fortunei var. cyclolepis, K. pubescens, K. davidiana var. calcarea, K. davidiana var. chienpeii, K. fortune, and K. evelyniana (Table 1). Two to three samples of each taxon were used in the evaluation of cross-amplification. All loci were transferable to the six taxa of Keteleeria (Table 4), and the annealing temperatures are listed on Table 2. However, few alleles were recovered in the six other Keteleeria taxa sampled here were caused by a small sampling effect. Nonetheless, the systematic and phylogeographic patterns among the Keteleeria species remain unclear. The transferability of these microsatellite loci among different Keteleeria taxa indicates the usefulness of these molecular genetic markers for interspecific and intraspecific research on such topics as phylogeography, speciation and introgression.
For conservation purposes, 16 new polymorphic microsatellite loci were developed from K. davidiana var. formosana. The genetic variation indices evaluated using these 16 polymorphic microsatellite loci for the two remaining populations of this endangered species are potentially useful for future studies that will focus on identifying distinct evolutionary units within the populations for conservation management. The interspecies transferability of these microsatellite loci may also be useful for future research aiming to reconstruct the phylogeographic patterns and the process of speciation among closely related species.
Availability of supporting data
The microsatellite sequences are available through the National Centre for Biotechnology Information (seehttp://www.ncbi.nlm.nih.gov/). The accession numbers on the repository are the following: GenBank accession number HG518488 through HG518503.
Chiang YC, Shih HC, Chang LW, Li WR, Lin HY, Ju LP: Isolation of 16 polymorphic microsatellite markers from an endangered and endemic species, Podocarpus nakaii (Podocarpaceae). Am J Bot. 2011, 98: e306-e309. 10.3732/ajb.1100229.
Ho CS, Chen CT, Ko YZ, Ge XJ, Hung KH, Hsu TW, Chiang YC: Isolation and characterization of 15 microsatellite loci in four endangered Amentotaxus species (Taxaceae). Am J Bot. 2012, 99: e157-e159. 10.3732/ajb.1100452.
Ju LP, Shih HC, Chiang YC: Microsatellite primers for the endangered beech tree, Fagus hayatae (Fagaceae). Am J Bot. 2012, 99: e453-e456. 10.3732/ajb.1200118.
Huang S, Chiang YC, Schaal BA, Chou CH, Chiang TY: Organelle DNA phylogeography of Cycas taitungensis, a relict species in Taiwan. Mol Ecol. 2001, 10: 2669-2681. 10.1046/j.0962-1083.2001.01395.x.
Huang S, Hsieh HT, Fang K, Chiang YC: Patterns of genetic variation and demography of Cycas taitungensis in Taiwan. Bot Rev. 2004, 70: 86-92. 10.1663/0006-8101(2004)070[0086:POGVAD]2.0.CO;2.
Ju LP, Kuo CC, Chao YS, Cheng YP, Gong X, Chiang YC: Microsatellite primers in the native perennial cycad, Cycas taitungensis (Cycadaceae). Am J Bot. 2011, 98: e84-e86. 10.3732/ajb.1000504.
Fang JY, Chung JD, Chiang YC, Chang CT, Chen CY, Hwang SY: Divergent selection and local adaptation in disjunct populations of an endangered conifer, Keteleeria davidiana var. formosana (Pinaceae). PLoS One. 2013, 8: e70162-10.1371/journal.pone.0070162.
Frankham R, Ballou JD, Briscoe DA: Introduction to Conservation Genetics. 2002, Cambridge, UK:. Cambridge University Press
Schneider D, Wink M, Sporer F, Lounibos P: Cycads: their evolution, toxins, herbivores and insect pollinators. Naturwissenschaften. 2002, 89: 281-294. 10.1007/s00114-002-0330-2.
Frankham R: Genetics and conservation biology. C R Biol. 2003, 326: S22-S29.
Luscombe D, Yang Y: Keteleeria davidiana var. formosana. IUCN 2013. IUCN Red List of Threatened Species. 2013, 1-
Farjon A: A second revision of the genus Keteleeria Carrière (Taxonomic notes on Pinaceae II). Notes Roy Bot Gard Edinburgh. 1989, 46: 81-99.
Li HL, Hsuan K: Pinaceae. Flora of Taiwan, 1, Volume 2. Edited by: Committee of the Flora of Taiwan. 1994, Taipei, Taiwan, 568-569. 2
Chiang YC: Research report for 101 annual science and technology program of forestry bureau. Research Project: Identification the Genetic Units for Conservation and Restoration of Keteleeria davidiana var. formosana. 2012, Taipei, Taiwan: Council of Agriculture
Kanehira R: Pinaceae. Formosan trees: Indigenous to the island (revised). 1936, Taipei, Formosa: Department of Forestry, Government Research Institute, 39-41.
Zane L, Bargelloni L, Patarnello T: Strategies for microsatellite isolation: A review. Mol Ecol. 2002, 11: 1-6. 10.1046/j.0962-1083.2001.01418.x.
Liao PC, Gong X, Shih HC, Chiang YC: Isolation and characterization of eleven polymorphic microsatellite loci from an endemic species, Piper polysyphonum (Piperaceae). Conserv Genet. 2009, 10: 1911-1914. 10.1007/s10592-009-9852-x.
Hsu TW, Shih HC, Kuo CC, Chiang TY, Chiang YC: Characterization of 42 microsatellite markers from poison ivy, Toxicodendron radicans (Anacardiaceae). Int J Mol Sci. 2013, 14: 20414-20426. 10.3390/ijms141020414.
Benson G: Tandem repeats finder: A program to analyze DNA sequences. Nucleic Acids Res. 1999, 27: 573-580. 10.1093/nar/27.2.573.
Kalendar R, Lee D, Schulman AH: Java web tools for PCR, in silico PCR, and oligonucleotide assembly and analysis. Genomics. 2011, 98: 137-144. 10.1016/j.ygeno.2011.04.009.
Chiang YC, Shih HC, Huang MC, Ju LP, Hung KH: Characterization of microsatellite loci from Litsea hypophaea (Lauraceae), a tree endemic to Taiwan. Am J Bot. 2012, 99: e251-e254. 10.3732/ajb.1100551.
Peakall R, Smouse PE: GenAlEx 6.5: genetic analysis in Excel. Population genetic software for teaching and research-an update. Bioinformatics. 2012, 28: 2537-2539. 10.1093/bioinformatics/bts460.
Excoffier L, Lischer HEL: Arlequin suite ver 3.5: A new series of programs to perform population genetics analyses under Linux and Windows. Mol Ecol Resour. 2010, 10: 564-567. 10.1111/j.1755-0998.2010.02847.x.
Yu JH, Chen CM, Tang ZH, Yuan SS, Wang CJ, Zu YG: Isolation and characterization of 13 novel polymorphic microsatellite markers for Pinus koraiensis (Pinaceae). Am J Bot. 2012, 99: e421-e424. 10.3732/ajb.1200145.
Hung KH, Lin CY, Huang CC, Hwang CC, Hsu TW, Kuo YL, Wang WK, Hung CY, Chiang TY: Isolation and characterization of microsatellite loci from Pinus massoniana (Pinaceae). Bot Stud. 2012, 53: 191-196.
Ge XJ, Hsu TW, Hung KH, Lin CJ, Huang CC, Huang CC, Chiang YC, Chiang TY: Inferring mutiple refugia and phylogeographical patterns in Pinus massoniana based on nucleotide sequence variation and DNA fingerprinting. PLoS One. 2012, 7: e43717-10.1371/journal.pone.0043717.
Mariette S, Chagné D, Lézier C, Pastuszka P, Raffin A, Plomion C, Kremer A: Genetic diversity within and among Pinus pinaster populations: comparison between AFLP and microsatellite markers. Heredity. 2001, 86: 469-479. 10.1046/j.1365-2540.2001.00852.x.
We thank Dr. Xun Gong and Dr. Jianwen Dong for their assistance in collecting the related species. We also thank two anonymous reviewers for their helpful comments. This work was supported by grants from the National Science Council, Taiwan (NSC 98-2815-C-020-014-B, 98-2621-B-110-004-MY3, 100-2621-B-110-001-MY3, and NSC 102-2628-B-110 -001 -MY3) and the Forestry Bureau, Council of Agriculture, Republic of China (tfba-1010621) to Y-C Chiang.
The authors declare that they have no competing interests.
T-WH and Y-CC supervised the project. C-SH, H-CS, H-YL, S-TC, L-PJ, C-TC, T-WH and Y-CC collected plant sample in the field. H-CS, M-HC, Y-ZK, Y-SS and Y-CC mined the SSR primers. C-SH, Y-ZK, Y-SS and Y-CC analyzed the data. Y-CC wrote the manuscript. All authors read and approved the final manuscript.
About this article
Cite this article
Ho, CS., Shih, HC., Liu, HY. et al. Development and characterization of 16 polymorphic microsatellite markers from Taiwan cow-tail fir, Keteleeria davidiana var. formosana (Pinaceae) and cross-species amplification in other Keteleeria taxa. BMC Res Notes 7, 255 (2014). https://doi.org/10.1186/1756-0500-7-255
- Simple sequence repeat markers