Phylogeography of the common vampire bat (Desmodus rotundus): Marked population structure, Neotropical Pleistocene vicariance and incongruence between nuclear and mtDNA markers
© Martins et al; licensee BioMed Central Ltd. 2009
Received: 3 June 2009
Accepted: 20 December 2009
Published: 20 December 2009
The common vampire bat Desmodus rotundus is an excellent model organism for studying ecological vicariance in the Neotropics due to its broad geographic range and its preference for forested areas as roosting sites. With the objective of testing for Pleistocene ecological vicariance, we sequenced a mitocondrial DNA (mtDNA) marker and two nuclear markers (RAG2 and DRB) to try to understand how Pleistocene glaciations affected the distribution of intraspecific lineages in this bat.
Five reciprocally monophyletic clades were evident in the mitochondrial gene tree, and in most cases with high bootstrap support: Central America (CA), Amazon and Cerrado (AMC), Pantanal (PAN), Northern Atlantic Forest (NAF) and Southern Atlantic Forest (SAF). The Atlantic forest clades formed a monophyletic clade with high bootstrap support, creating an east/west division for this species in South America. On the one hand, all coalescent and non-coalescent estimates point to a Pleistocene time of divergence between the clades. On the other hand, the nuclear markers showed extensive sharing of haplotypes between distant localities, a result compatible with male-biased gene flow. In order to test if the disparity between the mitochondrial and nuclear markers was due to the difference in mutation rate and effective size, we performed a coalescent simulation to examine the feasibility that, given the time of separation between the observed lineages, even with a gene flow rate close to zero, there would not be reciprocal monophyly for a neutral nuclear marker. We used the observed values of theta and an estimated mutation rate for the nuclear marker gene to perform 1000 iterations of the simulation. The results of this simulation were inconclusive: the number of iterations with and without reciprocal monophyly of one or more clades are similar.
We therefore conclude that the pattern exhibited by the common vampire bat, with marked geographical structure for a mitochondrial marker and no phylogeographic structure for nuclear markers is compatible with a historical scenario of complete isolation of refuge-like populations during the Pleistocene. The results on demographic history on this species is compatible with the Carnaval-Moritz model of Pleistocene vicariance, with demographic expansions in the southern Atlantic forest.
The diversity of life in the Neotropics diversity has fascinated researchers worldwide. Many theories have been created to account for its extraordinary number of species, and these theories have been vigorously debated in the last decades . One of the most cited and controversial is the Refuge Theory. Originally proposed by Haffer , it suggests that the Pleistocene glaciation cycles would create contraction and subsequent expansion of forested areas that in turn would create allopatry between populations of the same forest-dwelling species, leading to intraspecific differentiation and consequently speciation. Paleopalinogical studies demonstrate that the Neotropical forested areas have been very dynamic: the Atlantic forest and the Amazon were connected in the past [3, 4], becoming separated as increasing aridity in the Tertiary triggered the formation of the belt of xeromorphic formations between them . There seemed to be a predominance of arboreal vegetation during most of the Pleistocene with Amazonic and Atlantic forest tree species in areas that today lie in the dry diagonal that separates these two biomes [see [6, 7]]. The extent to which these climatic fluctuations and associated vegetation changes affected the patterns of distribution and diversification of the fauna remains a central question in understanding the evolution of forest-associated taxa.
Species with broad geographical distribution in the Neotropics and a preference for forested habitats are specially suitable for the study of past fragmentation and connectivity between forested areas such as the Amazon and the Atlantic forest. The common vampire bat, Desmodus rotundus Geoffroy, 1810, is a species with a wide geographical distribution: it ranges from southern Mexico to northern Chile in the west, and the entire Brazilian and Uruguayan territory in the east . Throughout its extensive range it occurs from sea level to over 3500 m of altitude and has been captured in habitats as diverse as rainforests and semi-arid vegetation. Due to the nature of previous phylogeographic studies in D. rotundus, it is believed that this species relies on either caves or forested areas for roosting, with the possibility of being captured in open sites while foraging . Even with such a wide distribution, and with the existence of considerable morphological variation, there are no currently recognized subspecies for this taxon .
The common vampire bat feeds preferentially on medium to large sized mammals . It lives in colonies that generally consist of fewer than 100 individuals. These colonies are formed by a dominant male, a number of peripheral males and several groups of unrelated or distantly related females [10–12]. The dominant male expels the juvenile males from its natal colony, and it has been suggested that the females are philopatric [11, 12].
Previous studies that focused on the intraspecific genetic variation of the common vampire bat pointed to different phylogeographic scenarios. Baker et al.  analysed 22 allozyme loci in different colonies of D. rotundus bats in Central America and found low levels of genetic divergence both between individuals from the same and different colonies. Later, in a preliminary study carried in the Atlantic forest of Brazil, Ditchfield  analysed cytochrome b sequences from seven different individuals and found evidence for genetic structure together with sympatric haplotypes with high sequence divergence (4.0%). This value is much higher than the average divergence found for all other phyllostomid bats studied and suggets that, when studied on a larger scale, this bat could show phylogeographic structure for this marker.
Although largely used in phylogeographic studies, the mitochondrial DNA (mtDNA) is a single locus and bears limitations. Genealogies based on one single gene may be problematic because each reconstruction is just one point in the space of all possible genealogies. Depending on the genealogical history of the population studied, this approach may lead to erroneous inferences [15–19]. For example, loci under selection may deviate from the expected patterns of loci solely under demographic pressure and might mimic alternate demographic patterns . Selection acts locally in the genome, while it is expected that demography affect all neutral loci uniformly .
Species that present high differentiation at the mtDNA level at times appear to be one single deme when nuclear loci are included in the analysis. That has been described for Drosophila , reptiles  and bats . In the latter, the authors found population structure when analyzing distinct hibernation colonies of the bat Myotis myotis throughout the Alps using mtDNA, but no differentiation when microsatellites were used. The authors attributed the results to female philopatry and male-biased dispersal. A similar scenario was described for the ghost bat Macroderma gigas .
Given the factors explained above, this work aims at using mitochondrial and nuclear molecular markers to answer the following questions: (I) is there geographic structure for the common vampire bat D. rotundus? (II) Is this structure consistent across all markers? If not, what historical scenarios could be responsible for the incongruence? (III) Does the biogeographical pattern exhibited by the evolutionary lineages of this bat correspond to the expected distribution of lineages under the refuge hypothesis or any other theory created to account for the high diversity in the Neotropics?
Sampling, DNA extraction, amplification and sequencing
The Brazilian samples were collected by the authors under federal license issued by the Instituto Brasileiro do Meio Ambiente e dos Recursos Naturais - IBAMA (Brazilian Institute of Environment and Renewable Natural Resources) with 36 mm mistnets. The animals were sacrificed and all the samples are deposited in the Laboratório de Biologia Evolutiva e Conservação de vertebrados (LABEC) tissue bank. Samples from outside Brazil were donated by the American Museum of Natural History (AMNH), the Museum of Vertebrate Zoology (MVZ) and the Royal Ontario Museum (ROM).
Total DNA was obtained from ethanol-preserved liver or muscle and in two cases dried museum skin following Bruford et al. . The protocol and primers for amplification of the cyt b fragments are described elsewhere . The amplification and sequencing of the RAG2 gene followed the same protocol used for the mitochondrial gene. The primers used for amplification and sequencing of the RAG2 gene are described in Baker et al. . The primers used for amplification and sequencing of the DRB intron 5 are described in Kupfermann et al. . The eletropherograms were analyzed using the Sequence Navigator program and the Se-Al software  was used to align the DNA sequences by eye.
Phasing of nuclear genotypes and recombination tests
Heterozygous nucleotide positions were identified by conspicuous double peaks in the electropherograms of both L and R strands. We used a Bayesian approach implemented in the software PHASE  to identify the haplotypes of heterozygous individuals.
In order to detect possible recombination events in the RAG2 and DRB markers, we used the RDP2 software . This software applies six different methods for detecting recombination (see  and references therein). We used windows of 20, 50 and 100 bp in each of the scans.
Four different methods of phylogenetic inference were used. A maximum parsimony (MP) analysis was implemented, weighting all nucleotide changes and codon positions equally. Heuristic searches were run using random addition of taxa and tree bisection and reconnection algorithm (TBR), as implemented in the PAUP* software . Bootstrap support  was estimated using 1000 replicates with heuristic mode parsimony. The maximum likelihood (ML) analysis used the GTR model of nucleotide substitution with different base frequencies and a gamma shape parameter of 0.3373, chosen using the Modeltest software . The search for the best maximum likelihood tree also used the TBR algorithm. A distance-based analysis was implemented by using the substitution model described above to estimate distances and the neighbour joining algorithm to infer a phylogeny, also using the PAUP* software. Last, we used the Bayesian approach implemented in the Mr. Bayes software . This method used the same substitution model as the ML analyses. The program ran 5 × 106 generations until the two chains had a standard deviation inferior to 0.01. The tree saved with posterior probabilities on the nodes was visualized using the Treeview software .
To complement the tree-based approaches we have also implemented Nested Clade Analysis (NCA) using the TCS software  with a statistical parsimony algorithm [38, 39]. The haplotype network was used as input for the GEODIS software  and the result from this analysis was used on the software's inference key (release date: 11/11/2005).
Population level and coalescent analyses
Analysis of molecular variance (AMOVA ) was used to quantify the extent of population subdivision using the Arlequin software . Other summary statistics, as well as neutrality tests (carried over for testing the hypothesis of recent population expansion as expected under a Pleistocene refuge model) and dN/dS analyses were implemented using the DNAsp software .
To simultaneously estimate several population parameters, we used coalescent-based analyses. The following population parameters were estimated: time of separation between populations (t), θ (4Nμ), and a migration rate (m), when applicable. These analyses were implemented using the software MDIV  with a finite-sites model (HKY ) using the first 5 × 105 cycles as burn in and 5.5 × 106 total cycles. Different runs were carried out assuming either complete isolation between populations (m = 0) or populations exchanging migrants at a rate m (estimated by the software). In order to determine whether the model where m is estimated fits the data significantly better than the model of m = 0 we used the Akaike information criterion  to compare likelihoods, as described in Nielsen & Wakeley . We compared the coalescent-based estimates of t and θ obtained to estimates using a methodology based on the net-sequence divergence measure (calculated using MEGA4 software ; see [48, 49]) as described in Edwards & Beerli . Two different mutation rates estimated for small-bodied mammals were used to calculate divergence times as applied by Hoffman et al.  for phyllostomid bats: 2.6% and 5.0% per million years. Generation time was assumed to be one year. We have also estimated a substitution rate for the RAG2 gene based on the data published by Baker et al.  on the phylogeny of the Phyllostomid family. We used Modeltest to estimate the model of nucleotide substitution and calculated the number of substitution using the origin of the family at 32 My ago as a calibration point as implemented by Teeling et al. . The mutation rate was estimated at 1.94 × 10-9 substitutions per site per year, similar to other mutation rates for nuclear genes in mammals .
Coalescent simulations of nuclear sequence data
According to Hare , the average nuclear marker will not trace Pleistocene events in the same manner as mitochondral markers due to its larger effective size, lower mutation rates and diploid nature. If that statement is correct then the nuclear marker used in this study may be non-informative in phylogeographic analyses, even in the case of null or negligible gene flow between the demes. In order to examine the compatibility of a scenario where (I) there is null or negligible gene flow between two demes and (II) the mitochondrial analysis yields reciprocally monophyletic clades and the nuclear markers shows little to no phylogenetic signal, we decided to conduct a coalescent simulation of the nuclear data. We used the estimates of time of separation between the lineages and of mutation rates to generate DNA data for three populations. The next step is to use this data for generating phylogenetic trees and check for presence/absence of reciprocal monophyly of the populations. If the demographic scenario estimated for the observed data (with the given values of t and θ) generates no monophyly in ≥ 95% of the simulations, we can demonstrate that a pleistocenic scenario of complete isolation generates dissimilar results between mitochondrial and nuclear markers, as suggested by Hare .
We used two softwares that implement coalescent simulations to test this hypothesis. The first one is the SIMCOAL2 software . This program allows for complex demographic models including population split, arbitrary migration rates and population contraction and expansion. The values used for this simulation were the ones estimated for the nuclear loci analyzed in this study: we used a θ value of 3.66, the mutation rate calculated for the RAG2 gene (1.94 × 10-9 substitutions per site per year) and estimated the same number of chromosomes sampled for the actual data (24 for AMC, 50 for SAF and 62 for NAF). We ran 1000 iterations. We also used the MLCOALSIM software . This coalescent simulator does not allow for very complex demographic settings: therefore we used an island model with negligible migration rates (m = 0.01), as null migration rates make the trees infinite and cause the program to crash. This simulation also used 1000 iterations. The outfiles generated by the software were converted to a format that was used as infile for the program TNT  using a LINUX script. The TNT program generated as outfile a matrix showing presence/absence of reciprocal monophyly for each of the simulated demes. On both simulations we used a biogeographic scenario based on the results obtained on the analyses implemented in this work: a split between the Amazon and the Atlantic forest around 900.000 YBP and a split between North and South Atlantic forest circa 500.000 YBP.
For the mitochondrial marker, 118 individuals were sampled from 54 localities. A sequence of 832 bp of the cyt b gene was obtained for each specimen, and 72 haplotypes were described for D. rotundus and two for Diphylla ecaudata (used as outgroup). A total of 233 nucleotide sites were variable, and 152 of these variable sites were parsimony informative.
For the RAG2 gene 88 individuals were sequenced for a fragment 774 bp long, for a total of 176 chromosomes sampled. The individuals sampled belong to all the major localities sampled in the mitochondrial analysis except the Pantanal area. A total of 43 variable sites were identified, fifteen of which were heterozygous. After the analyses with the PHASE software, 45 unique haplotypes were identified. The haplotypic diversity for this dataset is very high (h = 0.896), but the nucleotide diversity was very low (π = 0.00314). The recombination analyses revealed no recombination events for this marker. For the DRB1 gene, only 62 sequences were obtained, each 282 bp long. Again, no samples from the Pantanal region could be amplified and sequenced. Fourteen variable sites were identified, including a two base pair indel at the end of this intron. The analysis with the software PHASE identified 18 unique haplotypes. Once again very low levels of nucleotide variability were observed (π = 0.00458). The table with all individuals, its respective localities, haplotypes and Genbank accession numbers can be found in the Additional file 1.
Phylogenetic inference: mitochondrial marker
Phylogenetic inference: nuclear markers
With the extensive sharing of haplotypes and low nucleotide variability observed for the nuclear markers, all tree-based methods were non-informative. We therefore decided to implement unrooted networks to extract more information from the data. For the RAG2 gene, the analysis with the TCS software generated a very large number of ambiguous connections. In order to carry the subsequent analyses, single-copy non-internal haplotypes with multiple connections were removed from the data matrix. The final data matrix generated the network used as input in the GEODIS software. The network for the DRB intron used all haplotypes. The hierarchical grouping of clades was done manually. The analyses carried with the GEODIS software found no statistical correlation between clades and geography for both markers, meaning that all analyses were inconclusive.
Population-level analyses and coalescent estimates
Pairwise divergence times for the mtDNA clades described.
Time estimate (coalescent)
Time estimate (non-coalescent)
SAF × NAF
SAF × AMC
SAF × PAN
NAF × AMC
NAF × PAN
AMC × PAN
CA × SAF
CA × NAF
CA × AMC
CA × PAN
For the RAG2 gene, the mean nucleotide divergence in the sample was very low (0.6%) and there is an extensive sharing of haplotypes between very distant localities, including the localities belonging to distinct mitochondrial lineages. At the same time, there is an equal number of exclusive haplotypes for each of these regions.
Population genetics analysis results for the nuclear marker RAG2.
The results obtained for dN/dS analyses for the nuclear markers showed approximately ten synonymous substitutions for each non-synonymous, congruent with the expected pattern for neutral markers described in the literature .
The coalescent analyses using the MDIV software was implemented in the RAG2 dataset and yielded a time of divergence of 500,000 years (95% confidence interval 229,000-940,000 years) between the two Atlantic forest clades and 900,000 years (95% confidence interval 524,000-1,400,000 years) between the Atlantic forest clades and the AMC clade, values congruent with the ones obtained for the mitochondrial marker and that were used in the coalescent simulations. We decided to not implement the coalescent analyses on the DRB intron due to a lack of reliable mutation rate estimates.
The simulations generated by the sofware SIMCOAL yield unrealistic results. The outfiles were run in the Arlequin software and in all the outfiles generated the number of variable sites was equal to the number of nucleotide sites. Given these conditions, each chromosome was a distinct haplotype. These results are completely incongruent with the observed data and had to be discarded from the analyses.
Coalescent simulations and clade monophyly.
No reciprocal monophyly
Reciprocal monophyly for SAF
Reciprocal monophyly for NAF
Reciprocal monophyly for AMC
Reciprocal monophyly for the three clades
It comes as no surprise that species with wide geographical distributions consist of two or more evolutionary units when molecular markers are studied. In the case of the common vampire bat, the population subdivision found for the mitochondrial marker does not translate into colonies or hibernation sites, as has been described for other bats (see  and ). The subdivision corresponds to different ecodomains and has correspondence in other vertebrate organisms. Therefore we are dealing with a complex pattern that needs to be carefully discussed.
The common vampire bat can fly 20 km from roost to feeding site in a single night . It inhabits all the biomes that exist in the Neotropics, from seasonally flooded forests to semi-arid environments, and from sea level to 3600 meters of altitude. Therefore there are no identifiable physical barriers for dispersal and gene flow in the distributional range of this bat. In this case, it is likely that ecological barriers serve as a possible explanation to the structure found for this bat. All the divergence times estimated using coalescent and non-coalescent approaches fall within the Pleistocene epoch, suggesting that this bat is indeed susceptible to forest fragmentation.
The east/west separation detected by the mitochondrial marker in D. rotundus is very clear and coincides with the biogeographic provinces described by Koopman  for South America, but unfortunately, when using this marker, the relationship between the western clades remains obscure. On the one hand, the high nucleotide divergence to the outgroups of the analysis (up to 30%) may have influenced the topology of the tree. On the other hand the lack of resolution regarding the phylogenetic relationships between these three clades may reflect concomitant historical events of divergence. In this case, there would be an interesting corridor for gene flow along the eastern slopes of the Andes cordillera, that generated the clade formed by PAN and CA that showed high bootstrap support as shown in Figure 1. The other ancestral clade would be formed along the Guiana and Brazilian shields (that originated the AMC and the Atlantic forest clades).
The high levels of sequence divergence may also be responsible for low bootstrap values showed by the AMC clade; but since this clade was consistently present with the exact same topology in all methods employed in this study, we consider this clade to be reflecting the true relationship between the sampled haplotypes. Since two of the three basal haplotypes are peripheral in the clades geographical distribution (one in Venezuela and one in Marajó Island, at the mouth of the Amazon river), it is possible that fine-scale sampling could reveal a better resolution of this very large geographical area.
The Atlantic forest has become separated of the other area clades during the Pleistocene epoch, a result that is congruent with the appearance of a dry belt separating this forested area from the Amazon. This is the first study to suggest a Pleistocene separation of the Atlantic forest and the Amazon based on the phylogeographic data collected on a vertebrate species using molecular-based estimates of divergence times.
This work describes the Atlantic forest of Brazil as a composite area, with northern and southern components. The latitudinal division of this area has been recognized using parsimonious analysis of endemicity in amphibians [62, 63], reptiles , birds  and harvestmen . More recently, phylogeographic studies that used mtDNA described this structure in organisms as diverse as birds , pit vipers , non-volant small mammals , canids  and the bat species Carollia perspicillata . All the studies that implemented time estimates yielded Pleistocene divergence times for this event. There are numerous paleopalinological and sediment studies that describe that this region has been fragmented with dry open areas related to glaciation-driven events during the Pleistocene [71–74].
In a recent study, Carnaval and Moritz  generated climatic simulation data and cross-referenced their results with phylogeographic and paleopalinological studies. The authors describe a scenario where in its northern portion the Atlantic forest have always supported an evergreen forest even during the driest conditions, while south of the Doce river the climatic conditions would not support a forest formation. The authors suggest that in its current southern distribution, the Atlantic forest was probably fragmented in several small patches in the wettest areas, a scenario that was proposed before by Whitmore and Prance . The results shown here are congruent with this scenario: the SAF clade is the only one with significant evidence of population expansion for two different markers. According to Lessa et al. , refuges would not only create geographic structure that cannot be associated with conspicuous contemporaneous physical barriers on molecular markers, but these same markers should bear the footprints of a population expansion related to the end of the last glacial cycle. The results obtained here reflect the predictions of Pleistocene forest dynamics: the time divergence estimates all fall within this epoch and the estimated historical demography is congruent with refugia.
The existence of a Atlantic forest lineage that is basal to both SAF and NAF comes as a surprise, specially because this lineage is at the southernmost end of the Atlantic forest distribution - where according to the Carnaval-Moritz model, there should be no forests at the last glacial maximum. Two different scenarios can be hypothesized for this data: (I) retention of ancestral polimorphism associated with demographic expansion  from the hypothesized São Paulo refugium [9, 75] or (II) there could be more geographic structure regarding the common vampire bat that the sampling presented here allows discriminating. This particular haplotype could represent another lineage and a possible contact zone between this lineage and the SAF.
Mitochondrial and nuclear incongruence
The incongruence between the results observed for the mitochondrial and nuclear markers could be due to two different scenarios: complete fragmentation and incongruence between markers due to the nature and characteristics of each molecular marker or long-term female philopatry and male-biased dispersal. We will discuss each of these hypotheses in detail.
The first possibility - of complete isolation between populations but no footprints in nuclear markers due to larger effective size and lower mutation rates - was the reason behind the coalescent simulations carried in this study. The results of the simulations have shown that given the mutation rate, time of separation and effective population size calculated for these bats, an average nuclear marker might or might not reflect the true demographic history of this species with similar probabilities. Given that, we believe that the best way to approach the question on whether the structure found is valid or the outcome of long-term female philopatry would be to sequence a larger number of loci - at least 16 according to Moore  - in search of more accurate phylogenetic reconstructions or use other nuclear markers such as the Y chromosome or microsatellites, all of which are beyond the scope of the present analysis. The RAG2 marker was chosen based on the work by Lewis-Oritt et al.  that described relatively high levels of intraspecific divergence for this marker in Moormopidae bats (comparable in some cases to mtDNA intraspecific divergence in bats). The DRB intron was chosen due to the possibility of studying balancing selection on the common vampire bat - as described for many of the genes that comprise the major histocompatibility complex (MHC) in vertebrates. This work shows that these two markers do not seem appropriate for intraspecific phylogeographic studies.
Females of D. rotundus usually remain in their natal colonies after they reach maturity, as observed by Wilkinson , but they may also join a new group not far from their birthplace. In addition, females do occasionally migrate among roosts, so there is evidence of some adult female dispersal. Under these conditions of female philopatry with restricted dispersal, the expected pattern of differentiation would be one of isolation by distance . However, in our analyses we find low sequence divergence among localities at large geographic distances (over 400 km), and large increases in genetic differentiation taking place over relatively short distances. This pattern is particularly clear in the analysis of the Atlantic forest. The NAF and SAF possess similar geographical area and distances between the localities sampled. Both clades show low differentiation (maximum haplotype divergence being under 1% in SAF and under 3% in NAF) over a broad geographical range (see Figure 2), but are highly divergent from one another (mean divergence 6.6%), with haplotypes differing by 7% separated by less than 200 km. The NCA results point towards the same direction in opposition to isolation by distance. We believe that if the sampling gaps were filled, this outcome on the inference key (allopatric fragmentation) would be observed repeatedly among all clades. The results observed for the mtDNA marker show strong evidence of historical fragmentation, even with limited information on the extent of historical and current male-mediated gene flow.
The phylogeographic pattern described for the common vampire bat Desmodus rotundus is characterized by Pleistocene ecological vicariance. The mitochondrial marker showed deep divergence between reciprocally monophyletic clades representing distinct ecodomains within the Neotropics. There was also a clear East/West division within South America, where the coastal Atlantic forest was separated from the Amazon and the Pantanal by the Brazilian dry diagonal of open formations. The times of separation between the lineages are all within the Pleistocene epoch. The phylogenetic pattern is congruent with many other Neotropical clades. In addition, the historical demography, with a population expansion at the southern end of the distribution, is compatible with the Carnaval-Moritz model of historic Atlantic forest dynamics. The coalescent simulations showed that, given the population parameters estimated for this species, a nuclear marker may or may not recover Pleistocene population history with similar probabilities.
The next step in revealing the true nature of the interactions between the mitochondrial clades indentified and the species status within D. rotundus will be finding the exact locations of the contact zones between these clades and studying in detail the different kinds of interactions between the animals in these localities. This study would consider not only a multiloci approach but also field observations on ecology and behavior. We are also on the way to conduct experiments with bats genotyped for the different clades under controlled laboratory conditions, to test for reproductive isolation and the possible ramification of the hybridization between the animals belonging to different clades.
The authors thank the following researchers and/or institutions for collaborating with samples for this study: the Royal Ontario Museum (ROM) and Burton Lim; the American Museum of Natural History (AMNH) and Nancy simmons; Susi Pacheco; Ana Carolina Martins; Caroline Aires; Ives Arnone; Sérgio Althoff; and the agriculture offices of the states of São Paulo, Mato Grosso and Tocantins. We also would like to thank Dr. Fernando Marques for making the simulation analysis using TNT possible by developing the LINUX script; and thank Steven Woolley for assistance with computer-based analyses. Thanks for Adalberto Cesari, Leon Franciatto, Verônica Bueno, Mauro Cardoso Jr. and Paulo Noffs for help with figures. Also thanks for Adalberto Cesari for english revisions. Part of this work was carried out by using the resources of the Computational Biology Service Unit from Cornell University which is partially funded by Microsoft Corporation. We would like to thank two anonymous reviewers for detailed and helpful suggestions. This work became possible through finacial aid from FAPESP grants number 03/01583 -3 and 04/08682-4 and CAPES grant number BEX4687/06-0.
- Patton JL, da Silva MNF: Rivers, refuges ans ridges: the geography of speciation of Amazonian mammals. Endless forms: Species and speciation. Edited by: Berlocher S, Howard D. 1998, Oxford: Oxford University Press, 202-13.Google Scholar
- Haffer J: Speciation in Amazonian forest birds. Science. 1969, 165: 131-7. 10.1126/science.165.3889.131.View ArticlePubMedGoogle Scholar
- Oliveira-Filho AT, Ratter JA: A study of the origin of central Brazilian forests by the analysis of plant species distribution patterns. Edinburgh Journal of Botany. 1995, 52: 141-94.View ArticleGoogle Scholar
- Vivo M: Mammalian evidence of historical ecological change in the Caatinga semiarid vegetation of northeastern Brazil. J Comp Biol. 1997, 2: 65-73.Google Scholar
- Bigarella JJ, Andrade-Lima D, Riehs PJ: Consideracões a respeito das mudanças paleoambientais na distribuição de algumas espécies vegetais e animais no Brasil. An Acad Bras Cienc. 1975, 47: 411-464.Google Scholar
- Ledru M-P: Late quaternary environmental and climatic changes in central Brazil. Quat Res. 1993, 39: 90-8. 10.1006/qres.1993.1011.View ArticleGoogle Scholar
- De Oliveira PE, Franca-Barreto AM, Suguio K: Late Pleistocene/Holocene climatic and vegetational history of the Brazilian caatinga: the fossil dunes of the middle Sao Francisco River. Palaeogeog Palaeoclim Palaeoecol. 1999, 152: 319-37. 10.1016/S0031-0182(99)00061-9.View ArticleGoogle Scholar
- Koopman KF: Systematics and distribution. Natural History of Vampire Bats. Edited by: Greenhall AM, Schmidt U. 1988, Boca Raton: CRC Press, 7-18.Google Scholar
- Martins FM, Meyer D, Ditchfield AD, Morgante JM: Mitochondrial DNA phylogeography reveals marked population structure in the common vampire bat, Desmodus rotundus (Phyllostomidae). J Zool Syst Evol Res. 2007, 45: 372-8. 10.1111/j.1439-0469.2007.00419.x.View ArticleGoogle Scholar
- Wilkinson GS: Reciprocal food sharing in the vampire bat. Nature. 1984, 308: 181-4. 10.1038/308181a0.View ArticleGoogle Scholar
- Wilkinson GS: The social organization of the common vampire bat. II. Mating system, genetic structure, and relatedness. Behav Ecol Sociobiol. 1985, 17: 111-21.Google Scholar
- Wilkinson GS: Social organization and behavior. Natural History of Vampire Bats. Edited by: Greenhall AM, Schmidt U. 1988, Boca Raton: CRC Press, 85-98.Google Scholar
- Baker RJ, Honeycutt RL, Bass RA: Genetics. Natural History of Vampire Bats. Edited by: Greenhall AM, Schmidt U. 1988, Boca Raton: CRC Press, 31-40.Google Scholar
- Ditchfield AD: The comparative phylogeography of Neotropical mammals: patterns of intraespecific mitochondrial DNA variation among bats contrasted to nonvolant small mammals. Mol Ecol. 2000, 9: 1307-18. 10.1046/j.1365-294x.2000.01013.x.View ArticlePubMedGoogle Scholar
- Templeton AR: Out of Africa again and again. Nature. 2002, 416: 45-51. 10.1038/416045a.View ArticlePubMedGoogle Scholar
- Wakeley J: Coalescent Theory: An Introduction. 2008, Greenwood Village: Roberts and PublishersGoogle Scholar
- Pamilo P, Nei M: The relationships between gene trees and species trees. Mol Biol Evol. 1988, 5: 568-83.PubMedGoogle Scholar
- Templeton AR: Using haplotype trees for phylogeographic and species inference in fish populations. Environ Biol Fishes. 2004, 69: 7-20. 10.1023/B:EBFI.0000022870.76912.0c.View ArticleGoogle Scholar
- Hoelzer GA: Inferring phylogenies from mtDNA variation: mitochondrial-gene trees versus nuclear-gene trees revisited. Evolution. 1997, 51: 622-6. 10.2307/2411136.View ArticleGoogle Scholar
- Hare MP: Prospects for nuclear gene phylogeography. TREE. 2001, 16: 700-6.Google Scholar
- Fu Y-X, Li W-H: Coalescing into the 21st century: an overview and prospects of coalescent theory. Theor Pop Biol. 1999, 56: 1-10. 10.1006/tpbi.1999.1421.View ArticleGoogle Scholar
- DeSalle R, Templeton A, Mori I, Pletscher S, Johnston JS: Temporal and spatial heterogeneity of mtDNA polymorphisms in natural populations of Drosophila mercatorum. Genetics. 1987, 116: 215-223.PubMed CentralPubMedGoogle Scholar
- Godinho R, Mendonça B, Crespo EG, Ferrand N: Genealogy of the nuclear b- fibrinogen locus in a highly structured lizard species: comparison with mtDNA and evidence for intragenic recombination in the hybrid zone. Heredity. 2006, 96: 454-63. 10.1038/sj.hdy.6800823.View ArticlePubMedGoogle Scholar
- Castella V, Ruedi M, Excoffier L: Contrasted patters of mitochondrial and nuclear structure among nursery colonies of the bat Myotis myotis. J Evol Biol. 2001, 14: 708-20. 10.1046/j.1420-9101.2001.00331.x.View ArticleGoogle Scholar
- Worthington Wilmer J, Hall L, Barratt E, Moritz C: Genetic structure and male- mediated gene flow in the ghost bat (Macroderma gigas). Evolution. 1999, 53: 1582-91. 10.2307/2640903.View ArticleGoogle Scholar
- Bruford MW, Hanotte O, Brookfield JFY, Burke T: Single-locus and multilocus DNA fingerprinting. Molecular genetic analyses of populations A practical Approach. Edited by: Hoelzel AR. 1992, Oxford: IRL Press, 225-269.Google Scholar
- Baker RJ, Porter CA, Bussche Van Den RA: Systematics of the bats of the family Phyllostomidae based on RAG2 DNA sequences. Occasional Papers, Museum of Texas Tech University. 2000, 202: 1-16.Google Scholar
- Kupfermann H, Satta Y, Takahata N, Tichy H, Klein J: Evolution of MHC-DRB introns: implications for the origin of primates. J Mol Evol. 1999, 48: 663-74. 10.1007/PL00006510.View ArticlePubMedGoogle Scholar
- Rambaut A: Se Al. [http://tree.bio.ed.ac.uk/software/seal/]
- Stephens M, Smith N, Donnelly P: A new statistical method for haplotype reconstruction from population data. Am J Hum Gen. 2001, 68: 978-89. 10.1086/319501.View ArticleGoogle Scholar
- Martin DP, Williamson C, Posada D: RDP2: recombinations detection and analysis from sequence aligments. Bioinformatics. 2005, 21: 260-2. 10.1093/bioinformatics/bth490.View ArticlePubMedGoogle Scholar
- Swofford DL: PAUP*. Philogenetic Analysis Using Parsimony (* and Other Methods), version 4. 1998, Sunderland: Sinauer AssociatesGoogle Scholar
- Felsenstein J: Confidence limits on phylogeny: using the bootstrap. Evolution. 1985, 39: 783-91. 10.2307/2408678.View ArticleGoogle Scholar
- Posada D, Crandall KA: MODELTEST: testing the model of DNA substitution. Bioinformatics. 1998, 19: 817-8. 10.1093/bioinformatics/14.9.817.View ArticleGoogle Scholar
- Huelsenbeck JP, Ronquist F: Bayesian analysis of molecular evolution using MrBayes. Statistical Methods in Molecular Evolution. Edited by: Nielsen R. 2005, New York: Springer, 183-232. full_text.View ArticleGoogle Scholar
- Page RDM: TREEVIEW: An application to display phylogenetic trees on personal computers. Comp App Biosci. 1996, 12: 357-8.Google Scholar
- Clement M, Posada D, Crandall KA: TCS: a computer program to estimate gene genealogies. Mol Ecol. 1993, 9: 1657-9. 10.1046/j.1365-294x.2000.01020.x.View ArticleGoogle Scholar
- Templeton AR, Boerwinkle E, Sing CF: A cladistic analyses of phenotypic associations with haplotypes inferred from restriction endonuclease mapping. I. Basic theory and an analysis of Alcohol Dehydrogenase activity in Drosophila. Genetics. 1987, 117: 343-351.PubMed CentralPubMedGoogle Scholar
- Templeton AR, Sing CF: A cladistic association analysis of phenotypic association with haplotypes inferred from restriction endonuclease mapping. IV. Nested analyses with cladogram uncertainty and recombination. Genetics. 1993, 134: 659-69.PubMed CentralPubMedGoogle Scholar
- Posada D, Crandall KA, Templeton AR: GeoDis: A program for the Cladistic Nested Analysis of the Geographical Distribution of Genetic Haplotypes. Mol Ecol. 2000, 9: 487-8. 10.1046/j.1365-294x.2000.00887.x.View ArticlePubMedGoogle Scholar
- Excoffier L, Smouse PE, Quattro JM: Analysis of molecular variance inferred from metric distances among DNA haplotypes: application to human mitochondrial DNA restriction data. Genetics. 1992, 144: 2001-14.Google Scholar
- Schneider S, Roessli D, Excoffier L: Arlequin ver. 2000: a software for population genetics data analysis. 2000, Switzerland: Genetics and Biometry Laboratory, University of GenevaGoogle Scholar
- Rozas J, Sanchez-DelBarrio JC, Messeguer X, Rozas R: DNAsp, DNA polimorphism analyses by the coalescent and other methods. Bionformatics. 2003, 19: 2496-7. 10.1093/bioinformatics/btg359.View ArticleGoogle Scholar
- Nielsen R, Wakeley J: Distinguishing migration from isolation: a Markov Monte Carlo approach. Genetics. 2001, 158: 885-96.PubMed CentralPubMedGoogle Scholar
- Hasegawa M, Kishino H, Iano T: Dating of the human-ape splitting by a molecular clock of mitochondrial DNA. J Mol Evol. 1985, 21: 160-74. 10.1007/BF02101694.View ArticleGoogle Scholar
- Akaike H: Prediction and entropy. A celebration of statistics. Edited by: Atkinson AC, Fienberg SE. 1985, Berlim: Springer-Verlag, 1-24.View ArticleGoogle Scholar
- Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Mol Biol Evol. 2007, 24: 1596-99. 10.1093/molbev/msm092.View ArticlePubMedGoogle Scholar
- Nei M, Li W-H: Mathematical models for studying genetic variation in terms of restriction endonucleases. PNAS. 1979, 76: 5269-73. 10.1073/pnas.76.10.5269.PubMed CentralView ArticlePubMedGoogle Scholar
- Takahata N, Nei M: Gene genealogy and variance of interpopulational nucleotide sequences. Genetics. 1985, 110: 325-344.PubMed CentralPubMedGoogle Scholar
- Edwards SV, Beerli P: Perspective: gene divergence, population divergence and the variance in coalescence time in phylogeographic studies. Evolution. 2000, 54: 1839-54.PubMedGoogle Scholar
- Hoffmann FG, Owen JG, Baker RJ: mtDNA perspective of chromosomal diversification and hybridization in Peters' tent-making bat (Uroderma bilobatum: Phyllostomidae). Mol Ecol. 2003, 12: 2981-93. 10.1046/j.1365-294X.2003.01959.x.View ArticlePubMedGoogle Scholar
- Teeling EC, Springer MS, Madsen O, Bates P, O'Brien SJ, Murphy WS: A molecular phylogeny for bats illuminates biogeography and the fossil record. Science. 2005, 307: 580-4. 10.1126/science.1105113.View ArticlePubMedGoogle Scholar
- Kumar S, Subramanlan S: Mutation rates in mammalian genomes. PNAS. 2002, 99: 803-8. 10.1073/pnas.022629899.PubMed CentralView ArticlePubMedGoogle Scholar
- Laval G, Excoffier L: SIMCOAL 2.0: a program to simulate genomic diversity over large recombining regions in a subdivided population with a complex history. Bioinformatics. 2004, 20: 2485-7. 10.1093/bioinformatics/bth264.View ArticlePubMedGoogle Scholar
- Ramos-Onsins S, Mitchell-Olds T: Mlcoalsim: Multilocus Coalescent Simulations. Evolutionary Bioinformatics. 2007, 3: 41-4.Google Scholar
- Goloboff P, Farris S, Nixon K: TNT (Tree analysis using New Technology) (BETA). 2000, Tucumán, Published by the authorsGoogle Scholar
- Durand JD, Templeton AR, Guinand B, Imsiridou A, Bouvet Y: Nested clade and phylogeographic analyses of the chub, Leuciscus cephalus (Teleostei, Cyprinidae), in Greece: Implications for Balkan Peninsula biogeography. Mol Phylogenet Evol. 1999, 13: 566-80. 10.1006/mpev.1999.0671.View ArticlePubMedGoogle Scholar
- Li W-H: . Molecular Evolution. 1997, Massachussets: Sinauer Associates, IncGoogle Scholar
- Worthington Wilmer J, Moritz C, Hall L, Toop J: Extreme population structuring in the threatened ghost bat, Macroderma gigas: evidence from mitochondrial DNA. Proc R Soc Lond B. 1994, 257: 193-8. 10.1098/rspb.1994.0115.View ArticleGoogle Scholar
- Greenhall AM, Joermann G, Schmidt U: Desmodus rotundus. Mammalian Species. 1983, Published by the American Society of Mammalogists. New York, 1-6. 10.2307/3503895.Google Scholar
- Koopman KF: Zoogeography. Biology of Bats of the New World Family Phyllostomatidae, Part 1. Edited by: Baker RJ, Jones JK, Carter DC. 1976, Special Publications of The Museum, Texas Tech University, 10: 1-218.Google Scholar
- Müller P: Dispersal centres of terrestrial vertebrates in the Neotropical realm. Hague. Edited by: Junk W. 1973, 2: [Junk W (Series Editor): Biogeographica, vol 2.]Google Scholar
- Lynch JD: The amphibians of the lowland tropical forests. The south american Herpetofauna: Its origin, evolution and dispersal. Edited by: Duellman WE. 1979, Lawrence: University of Kansas Museum of Natural History, 189-216.Google Scholar
- Vanzolini PE: Distributional Patterns of South American lizards. Proceedings of a workshop on Neotropical distribution patterns: 12-16 January 1987; Rio de Janeiro. Edited by: Vanzolini PE, Hayer W. 1988, Academia Brasileira de Ciências, 317-42.Google Scholar
- Bates JM, Hackett SJ, Cracraft J: Area-relationships in the Neotropical lowlands: an hypotheses based on raw distribution of passerine birds. J Biogeography. 1998, 25: 783-93. 10.1046/j.1365-2699.1998.2540783.x.View ArticleGoogle Scholar
- Pinto-da-Rocha R, da Silva MB: Faunistic similarity and historical biogeography of the harvestmen of southern and southeastern Atlantic rainforest of Brazil. J Arach. 2005, 33: 290-9. 10.1636/04-114.1.View ArticleGoogle Scholar
- Cabanne GS, Santos FR, Miyaki CY: Phylogeography of Xiphorhynchus fuscus (Passeriformes, Dendrocolaptidae): vicariance and recent demographic expansion in southern Atlantic forest. Biol J Linn Soc. 2007, 91: 73-84. 10.1111/j.1095-8312.2007.00775.x.View ArticleGoogle Scholar
- Grazziotin FG, Monzel M, Echeverrigarauy S, Bonatto S: Phylogeography of the Bothrops jararaca complex (Serpentes: Viperidae): past fragmentation and island colonization in the Brazilian Atlantic Forest. Mol Ecol. 2006, 15: 3969-82. 10.1111/j.1365-294X.2006.03057.x.View ArticlePubMedGoogle Scholar
- Costa LP: The historical bridge between the Amazon and the Atlantic Forest of Brazil: a study of molecular phylogeography with small mammals. J Biogeography. 2003, 30: 71-86. 10.1046/j.1365-2699.2003.00792.x.View ArticleGoogle Scholar
- Tchaicka L, Eizirik E, Oliveira TG, Candido Jr JF, Freitas TRO: Phylogeography and population history of the crab-eating fox (Cerdocyon thous). Mol Ecol. 2007, 16: 819-38. 10.1111/j.1365-294X.2006.03185.x.View ArticlePubMedGoogle Scholar
- Behling H: Late quaternary vegetational and climatic changes in Brazil. Rev Paleobotany and Palynology. 1999, 99: 143-56. 10.1016/S0034-6667(97)00044-4.View ArticleGoogle Scholar
- Behling H: South and southeast Brazilian grasslands during Late Quaternary times: a synthesis. Palaeogeog Palaeoclim Palaeoecol. 2002, 177: 19-27. 10.1016/S0031-0182(01)00349-2.View ArticleGoogle Scholar
- Behling H, Lichte M: Evidence of dry and cold climatic conditions at glacial times in tropical southeastern Brazil. Quat Res. 1997, 48: 348-58. 10.1006/qres.1997.1932.View ArticleGoogle Scholar
- Lichte M, Behling H: Dry and cold climatic conditions in the formation of the present landscape in Southeastern Brazil: an interdisciplinary approach to a controversial topic. Z Geomorph. 1999, 43: 341-58.Google Scholar
- Carnaval AC, Moritz C: Historical climate modelling predicts patterns of current biodiversity in the Brazilian Atlantic Forest. J Biogeography. 2008, 35: 1187-1201. 10.1111/j.1365-2699.2007.01870.x.View ArticleGoogle Scholar
- Whitmore TC, Prance GT: . Biogeography and quaternary history in tropical America. 1987, Oxford: Oxford University PressGoogle Scholar
- Lessa, EP, Cook JA, Patton JL: Genetic footprints of demographic expansion in North America, but not Amazonia, during the late quaternary. PNAS. 2003, 100: 10331-4. 10.1073/pnas.1730921100.View ArticleGoogle Scholar
- Excoffier L, Ray N: Surfing during population expansion promotes genetic revolution and structurations. TREE. 2008, 23: 347-51.PubMedGoogle Scholar
- Moore WS: Inferring phylogenies from mtDNA variation: mitochondrial-gene versus nuclear-gene trees. Evolution. 1995, 49: 718-26. 10.2307/2410325.View ArticleGoogle Scholar
- Lewis-Oritt N, Porter CA, Baker RJ: Molecular systematics of the family Moormopidae (Chiroptera) based on Cytochrome b and Recombination activator gene 2 sequences. Mol Phyl Evol. 2001, 20: 426-36. 10.1006/mpev.2001.0978.View ArticleGoogle Scholar