- Short Report
- Open Access
Isolation and characterization of 11 novel microsatellite loci in a West African leaf-nosed bat, Hipposideros aff. ruber
BMC Research Notesvolume 7, Article number: 607 (2014)
Noack’s leaf-nosed bat, Hipposideros ruber, is a cryptic species within the Hipposideros caffer species complex. Despite a widespread distribution in Africa and being host to potentially zoonotic viruses, the genetic structure and ecology of H. ruber is poorly known. Here we describe the development of 11 novel polymorphic microsatellite loci to facilitate the investigation of genetic structure.
We selected 20 microsatellite sequences identified from high throughput sequence reads and PCR amplified these for 38 individuals, yielding 11 consistently amplifying and scorable loci. The number of alleles per locus ranged from two to 12, and observed heterozygosities from 0.00 to 0.865. No evidence of linkage disequilibrium was observed, and nine of the markers showed no departure from Hardy-Weinberg equilibrium. We demonstrate successful amplification in two closely related species and two divergent lineages of the H. caffer species complex.
These new markers will provide a valuable tool to investigate genetic structure in the poorly understood Hipposideros caffer species complex.
Noack’s leaf-nosed bat Hipposideros ruber (Noack, 1893) is one of two recognised cryptic species within the Hipposideros caffer (Sundevall, 1846) species complex. These bats are widespread throughout sub-Saharan Africa and among the most abundant mammals on the continent [1, 2]. Mitochondrial evidence has shown the existence of several deeply divergent lineages within the H. caffer complex, which most likely constitute more than the two species . Recently, they have been discovered to host viruses with zoonotic potential , emphasizing the need for knowledge about their ecology in order to gain insight into zoonotic processes and risk factors for public health. Microsatellites provide a powerful tool to investigate the poorly known ecology and life history of these bats, including genetic structure, social arrangements and mating systems. Assessment of nuclear gene flow through microsatellite analysis may thus help to shed light also on the taxonomy of this species complex.
We isolated and characterized 11 microsatellite loci from a single, exclusively West African mitochondrial lineage of Hipposideros ruber, determined by sequencing of the cytochrome b gene [lineage D; 3]. This lineage is henceforth called H. aff. ruber due to its distant evolutionary relationship to H. ruber s. str. from East Africa, and may represent a distinct species . Hipposideros aff. ruber has been previously identified in central Ghana  and seems to be the most abundant of the three main lineages of the H. caffer complex in this region (unpublished data). The markers described herein represent, to our knowledge, the first suite of microsatellites for an African hipposiderid bat.
DNA was extracted from wing tissue from eight individuals sampled from the Brong Ahafo and Volta regions in central Ghana. DNA was extracted using an innuPREP DNA mini kit (Analytik Jena, Jena, Germany). Five micrograms of pooled DNA from eight individuals was sent to AGRF (http://www.agrf.com.au), where high throughput sequencing was performed on a Roche GS FLX 454 sequencing machine as described elsewhere [6–8]. QDD 1.3  was used to screen for di- to hexanucleotide repeat motifs with a minimum of eight repeats. From the 1689 microsatellites identified, a total of 32 primer pairs flanking tetranucleotide repeats with 11–15 repeat motifs were designed using PRIMER3 . Twenty primer pairs for which the annealing temperatures were most similar for each primer were selected for initial amplification trials. Amplification products from these primer pairs were visualised by electrophoresis on an agarose gel, from which 13 pairs with strong, stutter-free amplification bands were selected for optimisation. Forward primers for these 13 loci were directly labelled with a fluorochrome at the 5′ end. Twelve of these loci were successfully amplified by polymerase chain reaction (PCR), with one discarded due to the excessive amplification of non-specific product. PCR conditions for these 12 loci were optimized and genotyping was performed on 38 individuals (16 females, 22 males) sampled in Brong Ahafo and Volta Regions.
PCRs were performed using 10–50 ng of template DNA and reagent concentrations as follows: 200 μM each dNTP, one unit reaction buffer, between 2.0 and 2.5 mM MgCl2, equal concentrations of forward and reverse primer (0.25–1.0 μM) and one unit Taq polymerase (see Table 1). PCR amplification consisted of an initial denaturation at 94°C for 3 min followed by six touchdown cycles of 94°C denaturation for 30 s, annealing for 30 s with temperatures decreased by 2°C per cycle (55–47°C, 60–50°C, or 65–55°C; Table 1), and polymerase extension step at 72°C for 45 s. An additional 35 cycles were conducted, of denaturation (94°C, 30 s), primer annealing (final touchdown temperature, 45 s), and polymerase extension (72°C, 45 s), followed by a final extension (72°C, 10 min). PCR products were electrophoresed using an ABI3130 Genetic Analyzer (Applied Biosystems, Foster City, CA, USA). Allele sizes were determined via manual inspection using the software PEAK SCANNER 1.0 (Applied Biosystems), followed by automated binning performed using TANDEM 1.09 . We reanalyzed 20% of individuals to evaluate data integrity. One locus (Hr3) was discarded due to high rounding error in the TANDEM analysis, indicating poor marker quality. MICRO-CHECKER 2.2.3 was used to assess the probability of scoring errors, allelic dropout and the presence of null alleles . No scoring errors or allelic dropout were detected, although there were potentially null alleles at loci Hr7 and Hr12. Locus Hr13 may also suffer from null alleles, though low allelic variability (one common and one rare allele) did not allow this to be confirmed (Table 1).
The program CERVUS was used to calculate number of alleles, observed (Ho) and expected (He) heterozygosities, and probabilities of identity . All 11 loci were determined to be polymorphic, with a range of 2–12 alleles per locus (Table 1). Tests for pairwise linkage disequilibrium and deviations from Hardy-Weinberg equilibrium with Bonferroni corrections were calculated using FSTAT 2.9.3 . Two loci (Hr7, Hr12) deviated significantly from the Hardy-Weinberg equilibrium with a homozygote excess (Table 1). No linkage disequilibrium was detected between any loci. The probability of identity for the 11 loci was low at 1.6E-10 overall, and 8.7E-10 and 3.1E-9 for the Brong Ahafo and Volta localities, respectively. Probability of sibling identity was 1.4E-4, 2.2E-4 and 2.8E-4 for overall, Brong Ahafo and Volta, respectively.
In order to explore utility in closely related taxa, we tested whether these loci could be amplified across four related taxa in the genus Hipposideros using the PCR conditions specified above (Table 2). All but one locus successfully amplified PCR product across the tested taxa.
These microsatellite loci provide useful resources for the study of population genetic structure of bats in the Hipposideros caffer complex, and likely also related species in this genus. These findings will help to address questions regarding connectivity, social behaviour, and zoonotic disease ecology in African leaf-nosed bats.
All animals were handled in accordance with Ghanaian legislation. Bat capture and sampling were authorized by permit from the Wildlife Division of the Ministry of Lands, Forestry and Mines in Ghana, and approved by the Macquarie University Ethics Committee. Exports were conducted under a state agreement between the Republic of Ghana and the Federal Republic of Germany, and to Australia with permission from the Department of Agriculture, Fisheries and Forestry.
Availability of the supporting data
The microsatellite sequences are available through the National Centre for Biotechnology Information (http://www.ncbi.nlm.nih.gov); GenBank accession numbers KM370156 – KM370168.
Brosset A: Chiroptères d’altitude du mont nimba (guinée). description d’une espèce nouvelle, hipposideros lamottei. Mammalia. 1984, 48: 545-556.
Wright G: Hipposideros caffer (chiroptera: hipposideridae). Mamm Species. 2009, 845: 1-9.
Vallo P, Guillén-Servent A, Benda P, Pires D, Koubek P: Variation of mitochondrial DNA in the hipposideros caffer complex (chiroptera: hipposideridae) and its taxonomic implications. Acta Chiropterologica. 2008, 10: 193-206. 10.3161/150811008X414782.
Pfefferle S, Oppong S, Drexler J, Gloza-Rausch F, Ipsen A, Seebens A, Müller M, Annan A, Vallo P, Adu-Sarkodie Y, Kruppa T, Drosten C: Distant relatives of severe acute respiratory syndrome coronavirus and close relatives of human coronavirus 229E in bats, Ghana. Emerg Infect Dis. 2009, 15: 1377-1384. 10.3201/eid1509.090224.
Vallo P, Benda P, Martínková N, Kaňuch P, Kalko EK, Červený J, Koubek P: Morphologically uniform bats hipposideros aff. ruber (hipposideridae) exhibit high mitochondrial genetic diversity in southeastern senegal. Acta Chiropterologica. 2011, 13: 79-88. 10.3161/150811011X578633.
Margulies M, Egholm M, Altman W, Attiya S, Bader J, Bemben L, Berka J, Braverman M, Chen Y, Chen Z: Genome sequencing in microfabricated high-density picolitre reactors. Nature. 2005, 437: 376-380.
Boomer J, Stow A: Rapid isolation of the first set of polymorphic microsatellite loci from the Australian gummy shark, mustelus antarcticus and their utility across divergent shark taxa. Conserv Gen Res. 2010, 2: 393-395. 10.1007/s12686-010-9274-6.
Gardner MG, Fitch AJ, Bertozzi T, Lowe AJ: Rise of the machines – recommendations for ecologists when using next generation sequencing for microsatellite development. Mol Ecol Res. 2011, 11: 1093-1101. 10.1111/j.1755-0998.2011.03037.x.
Meglécz E, Costedoat C, Dubut V, Gilles A, Malausa T, Pech N, Martin J-F: QDD: a user-friendly program to select microsatellite markers and design primers from large sequencing projects. Bioinformatics. 2010, 26: 403-404. 10.1093/bioinformatics/btp670.
Rozen S, Skaletsky H: Primer3 on the WWW for general users and for biologist programmers. Methods Mol Biol. 2000, 132: 365-386.
Matschiner M, Salzburger W: TANDEM: integrating automated allele binning into genetics and genomics workflows. Bioinformatics. 2009, 25: 1982-1983. 10.1093/bioinformatics/btp303.
van Oosterhout C, Hutchinson W, Wills D, Shipley P: MICRO-CHECKER: software for identifying and correcting genotyping errors in microsatellite data. Mol Ecol Notes. 2004, 4: 535-538. 10.1111/j.1471-8286.2004.00684.x.
Kalinowski ST, Taper ML, Marshall TC: Revising how the computer program CERVUS accommodates genotyping error increases success in paternity assignment. Mol Ecol. 2007, 16: 1099-1106. 10.1111/j.1365-294X.2007.03089.x.
Goudet J: FSTAT (version 1.2): a computer program to calculate F-statistics. J Hered. 1995, 86: 485-486.
This work was supported by the German Research Foundation and funds from Macquarie, Ulm and Bonn Universities. We would like to thank S. Oppong for valuable assistance with permits and logistics. We thank all those involved in fieldwork, in particular E. E. Nkrumah, P. Anti, and E. K. Badu. We thank the communities of Buoyem and Likpe Todome for their hospitality and for providing access to the sites.
The authors declare that they have no competing interests.
AS was responsible for the design of this study, supervision of the work and contributed to the interpretation of results. HB performed field sampling, data analysis and marker validation, and HB and PV drafted the manuscript. MG contributed to analysis of sequences. CD and MT coordinated field sampling and were responsible for the implementation of the study. All authors read and approved the final manuscript.