Translator Disclaimer
4 December 2015 Development of EST-SSR Markers in Barringtonia racemosa (Lecythidaceae) and Cross-Amplification in Related Species
Hongxian Xie, Yang Yuan, Xiaoting Fang, Ying Liu, Chao Yang, Jianhua Jin, Fengxiao Tan, Yelin Huang
Author Affiliations +

Barringtonia racemosa (L.) Spreng. (Lecythidaceae), the most widespread species of Barringtonia, is characteristic of coastal communities in the Paleotropics (Tomlinson, 1986; Prance, 2012). The flowers open at night and are pollinated by night-flying animals. The bark and crushed fruit contain saponins and are used as a fish poison (Payens, 1967). This tree has a distribution from the Pacific to East Africa (Prance, 2012); it is variously recorded as a mangrove associate and is abundant along tidal rivers and in areas subject to tide and salinity. It also occurs inland on the edge of peat swamp forests and on hillsides to altitudes of 200 m (Tomlinson, 1986). The fruit are buoyant because of the spongy fibrous pericarp, and the wide distribution is on account of its dispersal by water (Prance, 2012). However, the actual level of interpopulation migration has never been studied. Recent advances in molecular techniques, such as microsatellite markers, have created new opportunities for population genetic research that can be used to estimate the genetic variation and clarify the level of migration among populations of B. racemosa across its range. Transcriptome sequencing is an efficient method for acquiring expressed sequence tag–simple sequence repeat (EST-SSR) markers (Bouck and Vision, 2007) based on its low cost and high reliability. Here, we describe the development of 15 nuclear microsatellite markers for ongoing population genetic research in B. racemosa. Additionally, we tested the transferability of these markers for studying two related species, B. asiatica (L.) Kurz and B. acutangula (L.) Gaertn.; these three species comprise all of the wide-ranging, water-dispersed species in Barringtonia (Prance, 2012).

METHODS AND RESULTS

In this study, a three-month-old seedling of B. racemosa (from Dongzhai harbor, Hainan, China) was immediately frozen in liquid nitrogen and stored at −80°C for RNA extraction. The total RNA of B. racemosa was extracted using an improved cetyltrimethylammonium bromide (CTAB) method (Fu et al., 2005). The RNA quality and quantity were measured using an Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, California, USA). The mRNA was isolated using Oligotex-dT30 (TaKaRa Biotechnology Co., Dalian, China) after extraction and then fragmented ultrasonically. The cDNA libraries were prepared for sequencing following the Illumina protocol, and paired-end short read sequencing was done using the Illumina Genome Analyzer II system (Illumina, San Diego, California, USA; sequencing performed by BGI, Shenzhen, Guangdong Province, China) (Huang et al., 2012; Qiang et al., 2015). A total of 13,510,000 90-nucleotide paired-end reads were obtained and assembled using Trinity (release 20140413) with the default parameters (Grabherr et al., 2011). CAP3 (Huang and Madan, 1999) was further used to cluster similar contigs and obtain a high-quality reference genome with nonredundant unigenes (minimum identity = 99%). The analysis yielded 63,932 contigs with an average length of 589 bp, an N50 length of 808, and an average depth of coverage of 18.14×. The software MISA (Thiel et al., 2003) was used to detect microsatellites from all unigenes obtained in the study. We screened for SSR motifs containing two to six nucleotides with the minimum number of repeats as follows: six for dinucleotide and five for trinucleotide, tetranucleotide, pentanucleotide, and hexanucleotide. Altogether, 6524 SSR motifs were found, and 30 of them were selected at random to design primers using Primer3 (Rozen and Skaletsky, 1999), with the optimum conditions set at a length of 20 bp (18–22 bp), a temperature of 60.0°C (56–62°C), and a product size range of 120–250 bp.

Table 1.

Characteristics of 15 microsatellite loci isolated from Barringtonia racemosa.

t01_01.gif

Forty-three individuals of B. racemosa from three natural populations (Appendix 1) were used to evaluate the polymorphisms of the target microsatellite loci. In addition, individuals from one population of B. asiatica and one population of B. acutangula were also sampled to detect the efficiency of these markers in cross-species amplification (Appendix 1). Genomic DNA was extracted from silica gel–dried leaves using the DNA Extraction Kit (Magen, Guangzhou, China) following the manufacturer's protocol. PCR amplifications were performed in a final volume of 30 µL, containing 60 ng of genomic DNA, 1× PCR buffer (10 mM Tris-HCl [pH 8.4] and 1.5 mM MgCl2; TransGen Biotech Co., Beijing, China), 0.2 mM dNTPs (TransGen Biotech Co.), 0.5 µM of each primer (Life Technologies, Shanghai, China), and 1 unit EasyTaq DNA polymerase (TransGen Biotech Co.). The PCR reactions were carried out under standard conditions for all primers in a Bio-Rad PTC-200 thermocycler (Bio-Rad Laboratories, Hercules, California, USA) under the following conditions: initial denaturation was at 94°C for 4 min, followed by 35 cycles at 94°C for 40 s, then annealing for 45 s at the optimal temperature for each primer pair (from 56–60°C, see Table 1), and 72°C for 1 min, with a final extension of 10 min at 72°C. To test the utility of the primers, PCR products were detected on 1% agarose gels. Finally, a total of 15 of the 30 primer pairs were successfully amplified. The sequences of microsatellite loci were deposited into GenBank (accession no. KR822404–KR822418; Table 1). The SSR genotyping was performed using a Fragment Analyzer Automated CE System (Advanced Analytical Technologies [AATI], Ames, Iowa, USA) using the Quant-iT PicoGreen dsDNA Reagent Kit, 35–500 bp (Invitrogen, Carlsbad, California, USA). Raw data were analyzed and allele sizes and number of alleles per locus were called using PROSize version 2.0 software (AATI). GenAlEx version 6.5 software (Peakall and Smouse, 2012) was used to calculate the average number of alleles per locus (A), the observed heterozygosity (Ho), and the expected heterozygosity (He) of each microsatellite locus. Hardy–Weinberg equilibrium (HWE) and linkage disequilibrium were tested using GENEPOP 4.3 (Rousset, 2008). MICRO-CHECKER (van Oosterhout et al., 2004) was employed for testing scoring errors and null alleles.

In total, 15 out of 30 primer pairs amplified the expected products and displayed clear polymorphisms. In B. racemosa, the number of alleles per locus ranged from two to seven, the observed heterozygosity ranged from 0 to 0.933, and the expected heterozygosity ranged from 0 to 0.772 (Table 2). All 15 primer sets also successfully amplified SSR loci in B. asiatica and B. acutangula.

In B. racemosa, five out of the 15 polymorphic microsatellite loci showed significant deviations from HWE in population BRSH and BRKL, respectively (Table 2). Only one locus showed significant deviation from HWE in population BRDR (Table 2). Deviations from HWE may have been related to the presence of null alleles, although we found no null homozygotes. There was also no significant linkage equilibrium (P < 0.05) between locus pairs (Table 2). In B. asiatica, two out of the 15 polymorphic microsatellite loci showed significant deviations from HWE (Table 2), and no linkage disequilibrium was detected between any pair of loci. In B. acutangula, no loci deviated from HWE, and no linkage disequilibrium was detected between any pair of loci.

CONCLUSIONS

The 15 microsatellite markers developed here have been proven to be useful in the evaluation of the genetic diversity of B. racemosa. Cross-amplification of these microsatellite loci in B. asiatica and B. acutangula suggests that they will also be useful in studies of other species within Barringtonia. We are currently using these markers for inferring genetic diversity and spatial genetic structure in these species. In particular, we are using these markers to assess the impact of seed dispersal by water and to estimate seed dispersal distances for these species in coastal tropical landscapes of the Indo-West Pacific.

Table 2.

Results of initial primer screening in populations of Barringtonia racemosa, B. asiatica, and B. acutangula.a

t02_01.gif

LITERATURE CITED

1.

A. Bouck , and T. Vision . 2007. The molecular ecologist's guide to expressed sequence tags. Molecular Ecology 16: 907–924. Google Scholar

2.

X. H. Fu , Y. L. Huang , S. L. Deng , R. C. Zhou , G. L. Yang , X. W. Ni , W. J. Li , and S. H. Shi . 2005. Construction of a SSH library of Aegiceras corniculatum under salt stress and expression analysis of four transcripts. Plant Science 169: 147–154. Google Scholar

3.

M. G. Grabherr , B. 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

4.

X. Q. Huang , and A. Madan . 1999. CAP3: A DNA sequence assembly program. Genome Research 9: 868–877. Google Scholar

5.

Y. L. Huang , X. T. Fang , L. Lu , Y. B. Yan , S. F. Chen , C. C. Zhu , L. Hu , et al. 2012. Transcriptome analysis of an invasive weed Mikania micrantha. Biologia Plantarum 56: 111–116. Google Scholar

6.

J. P. D. W. Payens 1967. A monograph of the genus Barringtonia (Lecythidaceae). Blumea 15: 157–263. Google Scholar

7.

R. Peakall , 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

8.

G. T. Prance 2012. A revision of Barringtonia (Lecythidaceae). Allertonia 12: 1–16. Google Scholar

9.

Y. M. Qiang , H. X. Xie , S. T. Qiao , Y. Yuan , Y. Liu , X. G. Shi , M. Su , et al. 2015. Development of microsatellite markers for Carallia brachiata (Rhizophoraceae). Applications in Plant Sciences 3: 1400125. Google Scholar

10.

F. Rousset 2008. GENEPOP'007: A complete reimplementation of the GENEPOP software for Windows and Linux. Molecular Ecology Resources 8: 103–106. Google Scholar

11.

S. Rozen , and H. Skaletsky . 1999. Primer3 on the WWW for general users and for biologist programmers. In S. Misener and S. A. Krawetz [eds.], Methods in molecular biology, vol. 132: Bioinformatics methods and protocols, 365–386. Humana Press, Totowa, New Jersey, USA. Google Scholar

12.

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

13.

P. B. Tomlinson 1986. The botany of mangroves. Cambridge University Press, Cambridge, United Kingdom. Google Scholar

14.

C. van Oosterhout , W. F. Hutchinson , D. P. M. Wills , and P. Shipley . 2004. MICRO-CHECKER: Software for identifying and correcting genotyping errors in microsatellite data. Molecular Ecology Notes 4: 535–538. Google Scholar

Appendices

Appendix 1.

Voucher and location information for the species and populations used in this study. All voucher specimens are deposited at the herbarium of Sun Yat-sen University (SYSU), Guangzhou, China.

tA01_01.gif

Notes

[1] The authors thank S. Shi, S. He, and C. Zhong for their assistance in collecting plant materials. This study is supported by grants from the National Natural Science Foundation of China (41276107, 31200466), the Natural Science Foundation of Guangdong Province (2015A030313136), and Chang Hungta Science Foundation of Sun Yat-sen University.

Hongxian Xie, Yang Yuan, Xiaoting Fang, Ying Liu, Chao Yang, Jianhua Jin, Fengxiao Tan, and Yelin Huang "Development of EST-SSR Markers in Barringtonia racemosa (Lecythidaceae) and Cross-Amplification in Related Species," Applications in Plant Sciences 3(12), (4 December 2015). https://doi.org/10.3732/apps.1500080
Received: 2 July 2015; Accepted: 1 July 2015; Published: 4 December 2015
JOURNAL ARTICLE
PAGES


SHARE
ARTICLE IMPACT
Back to Top