A Set of Novel Microsatellite Markers Developed for a Distylous Species Luculia gratissima (Rubiaceae)

Luculia gratissima (Wall.) Sweet (Rubiaceae) is a perennial shrub distributed in the southeast margin of the Tibetan Plateau in southwest China and adjacent region of Nepal and Myanmar. The plant is a distylous species with reciprocally placed stigmas and anthers in each floral morph. By using the Fast Isolation by Amplified Fragment Length Polymorphism (AFLP) of Sequences Containing (FIASCO) repeats protocol, 19 primer sets were identified in two wild populations. Of these primers, 10 displayed polymorphisms and nine were monomorphic. The number of alleles per locus ranged from two to five, values for observed and expected heterozygosities ranged from 0.000 to 1.000 and from 0.289 to 0.760, with averages of 0.303 and 0.555, respectively. These microsatellite loci will facilitate further studies on breeding system, gene flow patterns, and population structure of L. gratissima and its allied species.


Introduction
Luculia gratissima (Wall.) Sweet (Rubiaceae) is an evergreen perennial shrub distributed in the southeast margin of the Tibetan Plateau in southwest China and adjacent region of Nepal and Myanmar at altitudes between 800 and 2400 m [1,2]. The species flowers from August to December, and produces several compact pink to white tubular flowers on each tight inflorescence that are sweetly fragrant. It is very interesting that the species is distylous with complementary positioning of stigmas and anthers in the two floral forms [3]. Previous studies suggest that the floral syndrome of distylous species would facilitate disassortative mating during pollinator visitation [4,5]. However, the potential functions of distyly are not well understood, and DNA microsatellites would provide a precise molecular marker to trace the characteristic of pollination patterns by parentage analysis. In this study, we have developed and characterized 19 microsatellite markers for L. gratissima using the Fast Isolation by Amplified Fragment Length Polymorphism (AFLP) of sequences containing repeats (FIASCO) [6], which will be used for further studies of breeding system, gene flow patterns, and population structure.

Results and Discussion
A total of 275 positive clones were captured, among these 102 clones (37%) were found to contain simple sequence repeats (SSR). Finally, 33 sequences contained SSR loci were selected for primer design. Nineteen microsatellite loci successfully amplified in L. gratissima for 33 microsatellite loci and 10 of them were polymorphic amplification, the remaining nine microsatellite loci were monomorphic as a result of capillary gel (Table 1). LG17 The number of alleles ranged from two to five in 24 individuals of the species sampled from the two natural populations. Values for H O and H E ranged from 0.000 to 1.000 and from 0.289 to 0.760, with averages of 0.303 (SD = 0.285) and 0.555 (SD = 0.133), respectively ( Table 2). For all 10 microsatellite loci, except LG2 and LG25, the genotypic frequencies showed significant deviation from Hardy-Weinberg equilibrium (HWE) (P < 0.01) indicating the possibility of null alleles and the disassortative mating of this distylous species. These microsatellite markers developed in our study will be a useful tool for further studies of population genetics, and will be used to assign parentage to seeds which will help us understand the characteristic of pollination patterns for this distylous plant.

Experimental Section
Genomic DNA samples of L. gratissima were extracted from silica-gel-dried leaves of three different individuals using a modified hexadecyltrimethylammonium bromide (CTAB) method [7]. The extracted DNA was dissolved in 30 µL TE buffer. The fast isolation by AFLP of sequences containing repeats (FIASCO) [6] was performed in this study. Total genomic DNA (approximate 250-500 ng) was completely digested with 2.5 U of Mse I restriction enzyme (New England Biolabs, Beverly, MA, USA), and then ligated to an Mse I AFLP adaptor (5′-TAC TCA GGA CTC AT-3′/5′-GAC GAT GAG TCC TGA G-3′) using T4 DNA ligase (Fermentas, Burlington, ON, Canada). The digested-ligated fragments were diluted in a ratio of 1:10, and 5 µL of them were used amplification reaction with adaptor-specific primers (5′-GAT GAG TCC TGA GTA AN-3′/5′-TTA CTC AGG ACT CAT CN-3′). The amplified DNA fragments (200-800 bp) were enriched by magnetic bead selection with a 5-biotinylated (AG) 15 , (AAG) 10 and (AC) 15 probe, respectively [6]. The Recovered DNA fragments were reamplified with Mse I-N primers. The purified PCR products using EZNA Gel Extraction Kit (Omega Bio-Tek, Guangzhou, China), were ligated into pBS-T II vector (Tiangen, Beijing, China), and then transformed into E. coli strain DH5α competent cells (TaKaRa, Dalian, Liaoning, China). The positive clones were picked out and tested using vector primers T3/T7 and primer (AC) 10 /(AG) 10 /(AAG) 7 respectively to select appropriate fragments which contained SSR. In other words, a set of tested PCR included three reactions was performed using T3 and T7, T3 and (AC) 10 , (AC) 10 and T7 as primers, respectively. The second set of tested PCRs was done using T3 and T7, T3 and (AG) 10 , (AG) 10 and T7 as primers, respectively. The last set of tested PCRs was done using T3 and T7, T3 and (AAG) 7 , (AAG) 7 and T7 as primers, respectively. All these PCR reactions had the same conditions: 95 °C for 3 min followed by 30 cycles at 94 °C for 45 s, 52 °C for 1 min, 72 °C for 1 min, and a final extension step at 72 °C for 7 min. The positive clones were captured for sequencing with an ABI PRISM 3730XL DNA sequencer (Applied Biosystems, Foster City, CA, USA). Sequences contained simple sequence repeat and enough flanking regions were selected for primer design using Primer Premier 5.0 program [8].
The data was analyzed by GENEPOP 4.0 [9], which included test of observed heterozygosity (H O ), expected heterozygosity (H E ), and departure from Hardy-Weinberg equilibrium (HWE) for the 10 polymorphic microsatellite loci.

Conclusions
In summary, 19 microsatellite markers have been specifically developed for L. gratissima in this study. The high discriminatory power of 10 polymorphic loci suggests that they should be suitable for survey of population structure and parentage analysis in this distylous species. These developed and characterized SSR markers for L. gratissima would also be useful for exploring genetic diversity and genetic structure of other species in Luculia.