Skip to main content

Molecular species identification boosts bat diversity

Abstract

The lack of obvious morphological differences between species impedes the identification of species in many groups of organisms. Meanwhile, DNA-based approaches are increasingly used to survey biological diversity. In this study we show that sequencing the mitochondrial protein-coding gene NADH dehydrogenase, subunit 1 (nd1) from 534 bats of the Western Palaearctic region corroborates the promise of DNA barcodes in two major respects. First, species described with classical taxonomic tools can be genetically identified with only a few exceptions. Second, substantial sequence divergence suggests an unexpected high number of undiscovered species.

Findings

A total of at least 1.5 million species of organisms are described, but an estimated number of 10 to 100 million species may live on our planet [1, 2]. This huge gap gave rise to search for new technologies to describe the biological diversity. Meanwhile, genetic screening approaches involving a single or a small number of genes – often referred to as DNA barcoding – were started in diverse taxonomic groups for two major aims: (i) to establish a standardized technique to identify species, and (ii) to detect new species in an efficient manner to bring us closer to the true number of species [36]. So far, DNA barcodes were primarily used for genetic species identification of taxonomically poorly studied taxa and geographic regions like the tropics [711]. Here we tested the applicability of DNA sequencing for species delineation and species identification in one of the best-known taxonomic groups in an intensively sampled geographic region.

A recently compiled list of all mammal species of the world specifies a total of 44 vespertilionid bats species for the Western Palaearctic region [12]. One additional species was recently described (Pipistrellus hanaki [13]) and another was given species rank (Plecotus christii [14]). In this study we sequenced 900 base pairs of the mitochondrial nd1 gene [15] of 534 bats from 41 of the 46 described species [see Additional file 1]. This revealed species-specific DNA sequences that differed at least by 4% (uncorrected p distance) from sequences from other species in 816 of 820 pairwise species comparisons. Intraspecific variation was much lower than interspecific variation, which allowed a reliable genetic identification of most species (figure 1). Only four pairs of species showed genetic distances below 2.5% and species did not split in two monophyletic groups in all four pairs of species. Individuals of different species within these four species pairs occasionally shared the same sequence type. This lack of substantial sequence divergence between two species can have different reasons. Recent speciation or introgression of mitochondrial haplotypes are likely explanations for similar mitochondrial DNA sequences within the three species pairs Eptesicus serotinus/E. nilssonii [15], Myotis myotis/M. oxygnathus [16, 17] and possibly also for Pipistrellus kuhlii/P. deserti. In contrast, the observed genetic similarity of only 0.3% sequence divergence and the similarity in morphology and echolocation calls between Hypsugo ariel and H. bodenheimeri is due to the fact that both taxa are likely to represent one single species as morphological differences are restricted to an additional cusp at the second upper incisor in H. bodenheimeri [18]. Therefore, we treated both taxa as a single species. In total, nd1 sequences of 40 investigated species allowed the unambiguous genetic identification of 34 species, which equals to 85%.

Figure 1
figure 1

Phylogeny of western Palaearctic vespertilionid bats. Neighbor-joining (NJ) tree based on Kimura 2-parameter distances between mitochondrial nd1 sequences from 216 vespertilionid bats of the Western Palaearctic region. Thick branches lead to the different species while thin branches indicate intraspecific distances. Species names are encircled for species that were recently discovered with the help of DNA sequencing. Names of taxa for which a status as separate species is first proposed herein are highlighted by a black rectangle. Tadarida teniotis of the family Molossidae was used as an outgroup. Bootstrap support values of at least 70 are given above branches (5000 bootstrap replicates).

Substantial nd1 sequence divergence allowed also the distinction of several morphologically characterized taxa whose species status was questioned until recently [12, 19]. Sequences of Barbastella leucomelas from the type locality on the Sinai Peninsula differed on average by 13.3% from sequences obtained from European B. barbastellus. In addition, morphologically characterized subspecies within four species (Eptesicus bottae, E. serotinus, Plecotus auritus and P. kolombatovici [see Additional file 2]) showed mean genetic distances between 5.8 and 15.2%. These examples illustrate that morphologically similar taxa can be genetically well differentiated and that the combination of morphological and mitochondrial DNA sequence datasets is a powerful approach to resolve taxonomic uncertainties. The congruence between slight morphological differences and substantial sequence divergence of at least 5% indicates distinct species. Nevertheless, only the sympatric occurrence of both taxa can finally proof the existence of true biological species [20].

Highly divergent mitochondrial lineages were found in four species. Most astounding were two bats in central Europe (southern Austria and northern Italy) that had been identified as Myotis nattereri according to external morphology. They differ in their sequences from other central European M. nattereri by on average 9.7%. Among 108 individuals of Myotis mystacinus from all over Europe we found two individuals in Bulgaria that carried another mtDNA haplotype that differed on average by 10.2% from typical M. mystacinus. Bats originally assigned to Hypsugo savii split even into three highly divergent lineages [see Additional file 2]. Finally, in Pipistrellus kuhlii a western (Europe, Morocco and Libya) and eastern (Israel, Syria and Iran) lineage can be distinguished, which differ on average by 5.2%.

Highly divergent mitochondrial lineages provide strong evidence for cryptic species diversity in western Palaearctic bats. Sequence differences of at least 5% clearly fall into the range of interspecific and even intergeneric differences in bats and other mammals [21]. They are concordant with genetic distances between recently discovered cryptic species pairs in bats [15, 2225]. A comparison of inter-specific genetic variation in the mitochondrial cytochrome b gene (cob) between sister taxa in four mammalian orders including bats led to the conclusion that sequence divergences of >5% are a powerful first indicator for unrecognized cryptic species [21]. The two mitochondrial protein coding genes nd1 and cob are known to evolve at a similar rate [16]. We thus applied a somehow arbitrary species delineation threshold of 5% nd1 sequence divergence, which suggests nine new bat species in the Western Palaearctic. This value drops to still six yet new species, if a more conservative threshold of 9% is applied. Both thresholds indicate at least three new species for Europe, one of them even occurring in central Europe. The validity of these newly proposed bat species will now be tested by analysing other characters like morphological, behavioural or ecological traits. However, the mitochondrial sequences constitute a reliable sorting criterion to search for species-specific phenotypic traits.

Although sequencing of the mitochondrial nd1 gene revealed evidence for substantial cryptic bat diversity in Europe, it is likely that further species exist that are morphologically and genetically difficult to recognize. First, mitochondrial lineages with less than 5% DNA sequence divergence as found within the species Plecotus macrobullaris, Pipistrellus pipistrellus and Otonycteris hemprichii might also represent distinct species [see Additional file 2]. Second, it is well known that DNA barcodes fail to distinguish recently diverged species [5]. A possible example are specimens of Myotis mystacinus bulgaricus from southeastern Europe (Myotis aurascens sensu Benda and Tsytsulina [26]) that differ from the nominate central European M. mystacinus in the chromosomal locations of nucleolus organizer regions [27] and in some morphological characters [26]. However, they show similar nd1 sequences that do not split into two separate lineages. Therefore, sequencing parts of the mitochondrial genome will always provide only the lower limit of true species diversity.

Prior to the emergence of DNA sequencing in taxonomy 37 morphological defined vespertilionid bat species were acknowledged for the Western Palaearctic realm. Now a total of 51 species can be distinguished according to mitochondrial nd1 sequences (figure 2). However, neither the morphological nor the genetic approach alone detected all species. It is the combination of both that results in a preliminary total of 54 species in the Western Palaearctic (figure 2). This clearly shows the need for an integrative taxonomy approach in which phenotypic and genetic data are combined [28]. Only external morphological characters will enable rapid species identification in the field, which is an essential requirement for behavioural, ecological and evolutionary studies or conservation programs (e.g. European Bat Agreement, EU Habitats Directive).

Figure 2
figure 2

Cryptic species diversity in western Palaearctic vespertilionid bats. Number of vespertilionid bat species in the Western Palaearctic according to a morphological, a mtDNA sequencing or an integrative approach. Patterns distinguish genetically not investigated species (grey), morphologically differentiated species with distinct DNA sequences (black), species without a distinct DNA sequences (hatched) and morphologically cryptic species (white).

The raise in species number from 37 to 54 corresponds to an increase by almost 50%. New species were discovered in all major genera. Since Western Palaearctic bat species are among the best-studied taxonomic groups in terms of number of studied individuals, geographic coverage and applied techniques, we anticipate an enormous rise in bat species diversity in less studied areas like for example the tropics.

Methods

Sampling

Morphological measurements and a wing tissue sample with a diameter of 5 mm were taken from each individual before it was released. Tissue samples were immediately stored in 80% ethanol and kept at room temperature. The number of specimens and geographic coverage varied among taxa and ranged from one to 112 individuals [see Additional file 1]. To avoid a bias due to unequal sampling intensity, we kept only five sequences per major mitochondrial lineage. Sequences from different localities were given priority and sequences from one locality were chosen at random. This resulted in a total of 217 sequences that were included in genetic distance calculations.

Genetic analysis

We sequenced 900 bp of the mitochondrial gene nd1 (NADH dehydrogenase subunit 1) because a large data set was already available and the nd1 gene in vespertilionid bats evolves at a similar evolutionary rate as the mitochondrial protein coding gene cytb [16]. Unfortunately, no data from bats were available for the mitochondrial gene for cytochrome oxidase, subunit 1 (cox1), which is commonly used in DNA barcoding. The protocols of DNA isolation, amplification and sequencing were previously published [15]. Sequence alignments were unequivocal. We found neither insertions, deletions nor stop codons that would have been expected for nuclear pseudogenes. Sequence differences of all divergent lineages were confirmed by sequencing 550 base pairs of a second mitochondrial gene, the gene for 16S rRNA, which was proposed as a standard DNA barcoding marker for vertebrates [29] [see Additional file 3]. Genetic distances (uncorrected p and Kimura 2-parameter with a gamma distribution shape parameter of 0.5) were calculated and illustrated in a neighbor-joining tree by using the computer program PAUP* 4.0b10 [30]. In the text we always refer to uncorrected p distances. All DNA sequences are available from GenBank: Accession numbers of previous studies [GenBank:AF065106, GenBank:AF401362, GenBank:AF401363–AF401366, GenBank:AF401369, GenBank:AF401372, GenBank:AF401373, GenBank:AF401375, GenBank:AF401376, GenBank:AF401378, GenBank:AF401382, GenBank:AF401390, GenBank:AF401393, GenBank:AF401395–AF401398, GenBank:AF401400, GenBank:AF401407, GenBank:AF401409–AF401411, GenBank:AF401414–AF401419, GenBank:AF401421, GenBank:AF401427, GenBank:AF401431–AF401436, GenBank:AF401439–AF401442, GenBank:AF401444, GenBank:AF401446, GenBank:AF401448, GenBank:AF401449, GenBank:AF401451, GenBank:AF401453–AF401457, GenBank:AF401461–AF401466, GenBank:AF401468–AF401477, GenBank:AF516269, GenBank:AF516270, GenBank:AF516272–AF516275, GenBank:AF516277, GenBank:AY027834–AY027837, GenBank:AY027840, GenBank:AY027846, GenBank:AY027848, GenBank:AY027851, GenBank:AY027853, GenBank:AY027856, GenBank:AY027858, GenBank:AY027859, GenBank:AY033950, GenBank:AY033955, GenBank:AY033984–AY033987, GenBank:AY699856, GenBank:AY699858, GenBank:AY699860] and of this study [GenBank:DQ914967–DQ915087].

References

  1. Wilson EO: The encyclopedia of life. Trends Ecol Evol. 2003, 18: 77-80. 10.1016/S0169-5347(02)00040-X.

    Article  Google Scholar 

  2. Heywood VH, Watson RT: Global Biodiversity Assessment. 1995, Cambridge, U.K. , Cambridge University Press

    Google Scholar 

  3. Blaxter ML: The promise of a DNA taxonomy. Philos Trans R Soc Lond B Biol Sci. 2004, 359: 669-679. 10.1098/rstb.2003.1447.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  4. Tautz D, Arctander P, Minelli A, Thomas RH, Vogler AP: A plea for DNA taxonomy. Trends Ecol Evol. 2003, 18: 70-74. 10.1016/S0169-5347(02)00041-1.

    Article  Google Scholar 

  5. Moritz C, Cicero C: DNA barcoding: promise and pitfalls. PLoS Biol. 2004, 2 (10): e354-10.1371/journal.pbio.0020354.

    Article  PubMed Central  PubMed  Google Scholar 

  6. Hebert PD, Gregory TR: The promise of DNA barcoding for taxonomy. Syst Biol. 2005, 54 (5): 852-859. 10.1080/10635150500354886.

    Article  PubMed  Google Scholar 

  7. Monaghan MT, Balke M, Gregory TR, Vogler AP: DNA-based species delineation in tropical beetles using mitochondrial and nuclear markers. Philos Trans R Soc Lond B Biol Sci. 2005, 360 (1462): 1925-1933. 10.1098/rstb.2005.1724.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  8. Hebert PD, Penton EH, Burns JM, Janzen DH, Hallwachs W: Ten species in one: DNA barcoding reveals cryptic species in the neotropical skipper butterfly Astraptes fulgerator. Proc Natl Acad Sci U S A. 2004, 101 (41): 14812-14817. 10.1073/pnas.0406166101.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  9. Smith MA, Fisher BL, Hebert PD: DNA barcoding for effective biodiversity assessment of a hyperdiverse arthropod group: the ants of Madagascar. Philos Trans R Soc Lond B Biol Sci. 2005, 360 (1462): 1825-1834. 10.1098/rstb.2005.1714.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  10. Smith MA, Woodley NE, Janzen DH, Hallwachs W, Hebert PD: DNA barcodes reveal cryptic host-specificity within the presumed polyphagous members of a genus of parasitoid flies (Diptera: Tachinidae). Proc Natl Acad Sci U S A. 2006, 103 (10): 3657-3662. 10.1073/pnas.0511318103.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  11. Blaxter M, Elsworth B, Daub J: DNA taxonomy of a neglected animal phylum: an unexpected diversity of tardigrades. Proc Roy Soc London, Biol Sci. 2004, 271 Suppl 4: S189-92.

    Article  Google Scholar 

  12. Simmons NB: Order Chiroptera. Mammal species of the World: a taxonomic and geographic reference. Edited by: Wilson DE, Reeder DM. 2005, Baltimore, MD , Johns Hopkins University Press, 312-529.

    Google Scholar 

  13. Benda P, Hulva P, Gaisler J: Systematic status of African populations of Pipistrellus pipistrellus complex (Chiroptera: Vespertilionidae), with a description of a new species from Cyrenaica, Libya. Acta Chiropterologica. 2004, 6 (2): 193-217.

    Article  Google Scholar 

  14. Benda P, Kiefer A, Hanak V, Veith M: Systematic status of African populations of long-eared bats, genus Plecotus (Mammalia: Chrioptera). Folia Zoologica. 2004, 53: 1-47.

    Google Scholar 

  15. Mayer F, von Helversen O: Cryptic diversity in European bats. Proc Roy Soc London, Biol Sci. 2001, 268 (1478): 1825-1832. 10.1098/rspb.2001.1744.

    Article  CAS  Google Scholar 

  16. Ruedi M, Mayer F: Molecular systematics of bats of the genus Myotis (Vespertilionidae) suggests deterministic ecomorphological convergences. Mol Phylogenet Evol. 2001, 21 (3): 436-448. 10.1006/mpev.2001.1017.

    Article  CAS  PubMed  Google Scholar 

  17. Berthier P, Excoffier L, Ruedi M: Recurrent replacement of mtDNA and cryptic hybridization between two sibling bat species Myotis myotis and Myotis blythii. Proc Roy Soc London, Biol Sci. 2006, 273 (1605): 3101-3109. 10.1098/rspb.2006.3680.

    Article  CAS  Google Scholar 

  18. Harrison DL: A new species of pipistrelle bat (Chiroptera Pipistrellus) from South Israel. Durban Mus Novit. 1960, 5 (19): 261-267.

    Google Scholar 

  19. Horácek I, Hanak V, Gaisler J: Bats of the Palaearctic region: A taxonomic and biogeographic review. Proceedings of the VIIIth European bat research symposium. Edited by: Woloszyn BW. 2000, Krakow, Poland , Institute of Systematics and Evolution of Animals, Vol. I Approaches to biogeography and ecology of bats: 11-157.

    Google Scholar 

  20. Mayr E: Animal species and evolution. 1963, Cambridge, Mass. , Harvard University Press

    Chapter  Google Scholar 

  21. Baker RJ, Bradley RD: Speciation in mammals and the genetic species concept. Journal of Mammalogy. 2006, 87: 643-662. 10.1644/06-MAMM-F-038R2.1.

    Article  PubMed Central  PubMed  Google Scholar 

  22. Barratt EM, Deaville R, Burland TM, Bruford MW, Jones G, Racey PA, Wayne RK: DNA answers the call of pipistrelle bat species. Nature. 1997, 387: 138-139. 10.1038/387138b0.

    Article  CAS  PubMed  Google Scholar 

  23. Castella V, Ruedi M, Excoffier L, Ibanez C, Arlettaz R, Hausser J: Is the Gibraltar strait a barrier to gene flow for the bat Myotis myotis (Chiroptera: Vespertilionidae)?. Mol Ecol. 2000, 9 (11): 1761-1772. 10.1046/j.1365-294x.2000.01069.x.

    Article  CAS  PubMed  Google Scholar 

  24. von Helversen O, Heller KG, Mayer F, Nemeth A, Volleth M, Gombkötö P: Cryptic mammalian species: a new species of whiskered Bat (Myotis alcathoe n. sp.) in Europe. Naturwissenschaften. 2001, 88 (5): 217-223. 10.1007/s001140100225.

    Article  CAS  PubMed  Google Scholar 

  25. Kiefer A, Mayer F, Kosuch J, von Helversen O, Veith M: Conflicting molecular phylogenies of European long-eared bats (Plecotus) can be explained by cryptic diversity. Mol Phylogenet Evol. 2002, 25 (3): 557-566. 10.1016/S1055-7903(02)00299-3.

    Article  CAS  PubMed  Google Scholar 

  26. Benda P, Tsytsulina KA: Taxonomic revision of Myotis mystacinus group (Mammalia: Chiroptera) in the western Palearctic. Acta Soc Zool Bohem. 2000, 64: 331-398.

    Google Scholar 

  27. Volleth M: Differences in the location of nucleolus organizer regions in European vespertilionid bats. Cytogenetics and Cell Genetics. 1987, 44: 186-197.

    Article  CAS  PubMed  Google Scholar 

  28. Will KW, Mishler BD, Wheeler QD: The perils of DNA barcoding and the need for integrative taxonomy. Syst Biol. 2005, 54 (5): 844-851. 10.1080/10635150500354878.

    Article  PubMed  Google Scholar 

  29. Vences M, Thomas M, van der Meijden A, Chiari Y, Vieites DR: Comparative performance of the 16S rRNA gene in DNA barcoding of amphibians. Front Zool. 2005, 2 (1): 5-10.1186/1742-9994-2-5.

    Article  PubMed Central  PubMed  Google Scholar 

  30. Swofford DL: PAUP*: Phylogenetic analysis using parsimony (*and other methods). 2002, Sunderland, MA , Sinauer, 4th

    Google Scholar 

Download references

Acknowledgements

We are especially grateful to the numerous people who provided tissue samples. Andrea Ross assisted in DNA sequencing. Otto von Helversen, Martina Nagy, Stefan Schuster, Björn Siemers, Michael Veith and two anonymous referees participated in early discussions or made important comments on the manuscript.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Frieder Mayer.

Additional information

Competing interests

The author(s) declare that they have no competing interests.

Electronic supplementary material

12983_2006_48_MOESM1_ESM.rtf

Additional file 1: Investigated species. A list of species, sample size and geographic origin of analysed bats from the Western Palaearctic realm. (RTF 135 KB)

12983_2006_48_MOESM2_ESM.rtf

Additional file 2: Taxonomic comments. Taxonomic details are provided for newly proposed species and divergent mitochondrial lineages. (RTF 85 KB)

12983_2006_48_MOESM3_ESM.pdf

Additional file 3: DNA sequence divergence. Mean pair-wise sequence divergences (uncorrected p distances) between Western Palaearctic vespertilionid bat species measured by sequencing the two mitochondrial genes nd1 and 16S rDNA. Sequence divergences between newly proposed cryptic species pairs are indicated. The two lowest genetic distances originate from the mean values of the species comparisons Eptesicus serotinus/E. nilssonii and Myotis myotis/M. oxygnathus. (PDF 381 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Mayer, F., Dietz, C. & Kiefer, A. Molecular species identification boosts bat diversity. Front Zool 4, 4 (2007). https://doi.org/10.1186/1742-9994-4-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1742-9994-4-4

Keywords