Development of Microsatellite Markers for Croomia Japonica and Cross
Total Page:16
File Type:pdf, Size:1020Kb
Scientia Horticulturae 161 (2013) 228–232 Contents lists available at ScienceDirect Scientia Horticulturae journal homepage: www.elsevier.com/locate/scihorti Development of microsatellite markers for Croomia japonica and cross-amplification in its congener a a,1 b c d a,∗ Ming Fang , Chen-Xi Fu , Cheng-Xin Fu , You-Lin Zhu , Akiyo Naiki , En-Xiang Li a Key Laboratory of Plant Resources, College of Life Sciences and Food Engineering, Nanchang University, Nanchang 330031, China b Key Laboratory of Conservation Biology for Endangered Wildlife of the Ministry of Education, College of Life Sciences, Zhejiang University, Hangzhou 310058, China c Key Laboratory of Molecular Biology and Gene Engineering, College of Life Sciences and Food Engineering, Nanchang University, Nanchang 330031, China d Laboratory of Botany, Graduate School of Education (Science), Okayama University, Okayama 700-8530, Japan a r t i c l e i n f o a b s t r a c t Article history: Croomia is a small monocotyledonous genus with eastern Asian–eastern North American floristic disjunc- Received 4 January 2013 tion. The Croomia species are endangered and in urgent need of conservation. In this study, we report Received in revised form 11 July 2013 the development and characterization of 11 polymorphic compound microsatellite markers from Croo- Accepted 12 July 2013 mia japonica. Moreover, transferability and polymorphism of these primers was tested across Croomia heterosepala and Croomia pauciflora. All of them were transferable and polymorphic in these species. The Keywords: number of alleles per locus ranged from 2 to 14 (mean: 7.7), 2 to 20 (mean: 7.3) and 2 to 16 (mean: 8.2), Croomia while the observed (and expected) heterozygosities ranged from 0.053 to 1.000 (0.053–0.918), 0.000 Genetic diversity Transferability to 1.000 (0.656–0.945) and 0.190 to 0.905 (0.251–0.937) in populations of C. japonica, C. heterosepala and C. pauciflora, respectively. Most loci in the Croomia populations deviated significantly from the HWE expectations. Significant heterozygosity deficiency was found but no bottleneck was detected in Croomia. These polymorphic markers will be useful tool to study the genetic diversity and the population genetic structure, evolution of Croomia species and for establishment of effective conservation strategies. © 2013 Elsevier B.V. All rights reserved. 1. Introduction systematic relationship among these Croomia species needed to be studied. In our previous study, inter-simple sequence repeat Croomia Torr. is a genus of a small monocotyledonous family (ISSR) markers and chloroplast DNA (cpDNA) haplotypes were Stemonaceae (APGIII, 2009). The genus Croomia comprises three used to study the genetic diversity and phylogeography of two species with biogeographically interesting distribution patterns: Asian species (Li et al., 2008). As dominant markers may cause Croomia pauciflora (Nut.) Tor. is distributed in Southeast of North biases in the estimates of genetic diversity and genetic differentia- America and the other two species, Croomia japonica Miq. and Croo- tion (Nybom, 2004) and cpDNA sequences usually cannot provide mia heterosepala (Bak.) Oku., are in East Asia (Rogers, 1982; Ji and enough polymorphic loci for analyzing, other polymorphic codom- Duyfjes, 2000). The two Asian species have adjacently distribution inant markers are required to study the genetic structure and in South Japan, and C. japonica also distribution on adjacent Asiatic evolution in Croomia. Owing to high level of polymorphism and mainland in East China. Nowadays, all the three species of Croomia codominant inheritance, microsatellites (SSR) are very suitable to are surviving in small range size and with small number of popula- assess population genetic diversity and gene flow (Liu et al., 2009). tions, so they are treated as “endangered” or “threatened” (Patrick In order to obtain sufficient working SSR primers pairs, a variety et al., 1995; Estill and Cruzan, 2001; Wang and Xie, 2004). With of methods for SSR isolation have been developed (Squirrell et al., elegant figure, Croomia species are capable to be cultivated as pre- 2003). Compound SSR primers have proved to be valuable tools cious ornamentals. Based on slight difference, two new species of genetic studies (Hayden et al., 2004), because they can be used Croomia were published (Kadota and Saito, 2010). Therefore, the for different markers anchored by the same type of compound repeat sequence (Lian et al., 2006). An approach for developing compound microsatellite markers, with substantial time and cost savings, was introduced (Lian et al., 2006). This method is increas- ∗ Corresponding author. Present address: College of Life Sciences and Food Engi- ingly used in studying population genetics (Inoue et al., 2012). Here, neering, Nanchang University, Nanchang 330031, China. we characterize 11 new compound microsatellite markers for Croo- E-mail addresses: [email protected], [email protected] (E.-X. Li). 1 mia species. Contributed equally to this work. 0304-4238/$ – see front matter © 2013 Elsevier B.V. All rights reserved. http://dx.doi.org/10.1016/j.scienta.2013.07.014 M. Fang et al. / Scientia Horticulturae 161 (2013) 228–232 229 Table 1 Location and sampling size of Croomia populations in this study. Population code Location Latitude, longitude, altitude (m) Sample size C. japonica ◦ ◦ JFL Fuliang County, Jingdezhen, China 29 36 11 N, 117 39 04 E, 650 20 ◦ ◦ JTO Tosashimizu, Japan 32 52 04 N, 132 58 59 E, 30 19 C. heterosepala ◦ ◦ HZJ Numakubo Town, Fujinomiya, Japan 35 11 62 N, 138 35 10 E, 200 24 ◦ ◦ HNB Nanbu Town, Fujinomiya, Japan 35 11 01 N, 138 27 58 E, 600 24 C. pauciflora ◦ ◦ PLO Lowndes County, Alabama, American 32 03 47 N, 86 44 11 W, 83 20 ◦ ◦ PHE Henry County, AL, USA 31 37 35 N, 85 12 56 W, 105 21 2. Materials and methods Gorley, 2001). The primer pairs of specific primer (IP1) and com- pound SSR primer were used as a compound SSR marker (Table 2). 2.1. Plant materials To obtain accurate data, the compound SSR primer (AC)6(AG)5, (TC)6(AC)5 or (GT)6(TC)5 was labeled with a fluorescent dye (6- One hundred and twenty eight individuals from 6 Croomia popu- FAM, HEX or TAMARA). Thirty individuals from six populations lations (Table 1), including two C. japonica populations (JFL and (five individuals per population) were used in the initial screen of JTO), two C. heterosepala populations (HZJ and HNB) and two C. the 126 primers. Polymerase chain reactions were performed in pauciflora populations (PLO and PHE), were used to evaluate lev- 15 l of reaction mixture containing 40 ng of genomic DNA, 1U of els of polymorphism to examine the effectiveness of the designed Taq polymerase (Takara, Dalian, Liaoning, China), 1.5 l of 10 × PCR primer pairs. Silica-dried samples of leaf material were placed in buffer, 0.8 l of dNTPs (2.5 mM each), and 0.5 l of each IP1 (10 M) ◦ − zip-lock plastic bags and stored at 76 C. and a compound SSR primer (AC)6(AG)5, (TC)6(AC)5 or (GT)6(TC)5. The PCR amplification conditions were as follows: initial dena- ◦ ◦ turation at 95 C for 5 min; 35 cycles of 30 s at 94 C, 30 s at the 2.2. Development of microsatellite markers optimized annealing temperature (Table 2), 45 s of elongation at ◦ ◦ 72 C, ending with a 10-min extension at 72 C. Fragment analysis Total genomic DNA was extracted from silica-gel dried leaf was performed on the MegaBACE1000 autosequencer (GE Health- material using the CTAB (Cetyltrimethylammonium Bromide) care Biosciences, Pittsburgh, PA, USA), and the data were scored method (Doyle, 1991). The compound microsatellite marker and compiled using Genetic Profiler version 2.2 (GE Healthcare technique based on a dual-suppression-PCR method (Lian Biosciences). et al., 2006) was used for developing microsatellite markers in this study. First, an adaptor-ligated DNA library was con- 2.3. Data analysis structed according to the protocol of Lian et al. (2001). Then genomic DNA (ca.1000 ng) from C. japonica digested with HaeIII The number of alleles (Na), observed (Ho) and expected restriction enzyme (Takara Biotechnology Co., Dalian, Liao- (He) heterozygosities, linkage disequilibrium (LD), and deviations ning, China). The restricted fragments were then ligated with a from Hardy–Weinberg equilibrium (HWE) were analyzed using specific blunt unequal-length adaptor (consisting of a 48-mer: 5 - GTAATACGACTCACTATAGGGCACGCGTGGTCGACGGCCCGGGCTG- GENEPOP version 4.0.7 (Rousset, 2008). CERVUS version 3.0.3 (Kalinowski et al., 2007) was employed to calculate the value of GT-3 and an 8-mer the 3 -end capped by an amino residue: polymorphic information content (PIC). An UPGMA (unweighted 5 -ACCAGCCC-NH2-3 ) by use of a DNA Ligation Kit (Takara pair group method with arithmetic mean) tree was constructed Biotechnology Co.). Subsequently, the fragments were ampli- with Poptree2 (Takezaki et al., 2010) software. Bootstrap resam- fied from the HaeIII DNA library using compound SSR primers pling (n = 1000) was performed to test dendrogram robustness. The (AC)6(AG)5, (TC)6(AC)5 or (GT)6(TC)5 and an adaptor primer AP2 Cornuet and Luikart (1996) program BOTTLENECK v1.2.02 was used (5 -CTATAGGGCACGCGTGGT-3 ). Polymerase chain reactions were to detect the bottleneck hypothesis. performed in a 50 l reaction mixture containing 1 l of the ligated products, 1.5U of Taq polymerase (TaKaRa Biotechnology Co.), 5 l of 10 × PCR buffer with MgCl2, 5 l of dNTPs (2.5 mM 3. Results each), 0.5 l of compound SSR primer (10 M) and AP2 primer (10 M) for each. PCR amplification conditions were as follows: 1 A total of 126 primer pairs were designed from microsatellite ◦ ◦ ◦ cycle of 9 min at 94 C, 30 s at 62 C and 1 min at 72 C; 5 cycles of sequences isolated from the microsatellite-enriched libraries.