- Research article
- Open Access
Evolution of nectarivory in phyllostomid bats (Phyllostomidae Gray, 1825, Chiroptera: Mammalia)
BMC Evolutionary Biologyvolume 10, Article number: 165 (2010)
Bats of the family Phyllostomidae show a unique diversity in feeding specializations. This taxon includes species that are highly specialized on insects, blood, small vertebrates, fruits or nectar, and pollen. Feeding specialization is accompanied by morphological, physiological and behavioural adaptations. Several attempts were made to resolve the phylogenetic relationships within this family in order to reconstruct the evolutionary transitions accompanied by nutritional specialization. Nevertheless, the evolution of nectarivory remained equivocal.
Phylogenetic reconstructions, based on a concatenated nuclear-and mitochondrial data set, revealed a paraphyletic relationship of nectarivorous phyllostomid bats. Our phylogenetic reconstructions indicate that the nectarivorous genera Lonchophylla and Lionycteris are closer related to mainly frugivorous phyllostomids of the subfamilies Rhinophyllinae, Stenodermatinae, Carolliinae, and the insectivorous Glyphonycterinae rather than to nectarivorous bats of the Glossophaginae. This suggests an independent origin of morphological adaptations to a nectarivorous lifestyle within Lonchophyllinae and Glossophaginae. Molecular clock analysis revealed a relatively short time frame of about ten million years for the divergence of subfamilies.
Our study provides strong support for diphyly of nectarivorous phyllostomids. This is remarkable, since their morphological adaptations to nutrition, like elongated rostrums and tongues, reduced teeth and the ability to use hovering flight while ingestion, closely resemble each other. However, more precise examinations of their tongues (e.g. type and structure of papillae and muscular innervation) revealed levels of difference in line with an independent evolution of nectarivory in these bats.
The diversity of feeding specialization of phyllostomid bats are unique among all mammals [1–7]. They range from insect-to diverse vegetable-feeding strategies, as well as omnivory, carnivory, and even blood-feeding [8–16]. This ecological diversification is accompanied by morphological, behavioural and physiological adaptations [4, 9, 17–32]. A striking example is specialization for nectarivory, with several species feeding primarily on nectar. These bats have the ability to hover in front of a plant, while drinking nectar with their elongated and extensile tongues adorned with brush-like papillae and grooves for ingestion of nectar [3, 26, 29, 30, 33–37]. They digest and metabolize nectar and pollen quickly [32, 38–44]. Phyllostomid bats represent the second largest chiropteran family after the vesper bats (Vespertilionidae Gray, 1821), with more than 150 species in at least 49 genera. Their distribution ranges from southern Arizona and the West Indies to northern Argentina .
Although phylogenetic analyses of morphology, chromosomes, and molecules have helped to illuminate relationships among many genera and subfamilies of phyllostomid bats, relationships among nectarivorous genera are still unclear. Many phylogenies based on morphological characters suggest a monophyletic origin for all specialized nectarivorous phyllostomids [5, 46, 47]. We provide a well-supported phylogenetic estimate of phyllostomid bats based on a large molecular data set, comprising 10396 bp from a total of twelve nuclear-and mitochondrial genes, and try to clarify phylogenetic relationships among nectarivorous taxa by testing whether they share a close common ancestry. Furthermore, we used a molecular clock approach to evaluate the evolutionary time frame of diversification in phyllostomid bats.
Phylogeny of the Phyllostomidae
Figure 1 shows our maximum-likelihood reconstruction (b) compared to the latest comprehensive analysis of phyllostomid phylogeny (a) after Baker et al. (2003) . Baker and colleagues used sequences of 12S and 16S rRNA, tRNA Valin and the rag2 gene for their inference. Our reconstruction shows high congruence, even though it is completely based on independent genes (see methods section: Alignment 1). Although no members of the subfamilies Lonchorhinae, Glyphonycterinae and Rhinophyllinae were included (because of incomplete data for these taxa), major branching patterns were consistently reconstructed. Our reconstruction received good bootstrap support and is in line with Baker et al.'s phylogeny of phyllostomid bats. Therefore, we combined our data with the data from Baker et al. (2003) .
A separate analysis of all mitochondrial and nuclear loci (Alignment 2&3) resulted in high congruent phylogenies (Figure 2). Among the frugivorous species relationships changed between the independent inferences. A sister-group relationship between Carolliinae and Glyphonycterinae could not be inferred from the mitochondrial data set. In this reconstruction glyphonycterids were found basal to all frugivores. But this relationship obtained low support (BS 50) compared to the reconstruction based on nuclear loci, where Carolliinae is sister taxon to them (BS 73).
Maximum likelihood (ML) analysis based on our supermatrix (see methods section: Alignment 5) revealed a well-resolved phylogeny for the Phyllostomidae (Figure 3), with most nodes receiving high bootstrap support (BS > 90). Monophyly of all subfamilies recognized by Baker et al. (2003)  was verified, and relevant nodes were highly supported by different measurements (Table 1).
Three basal lineages, comprising the taxa Macrotus (1), Micronycteris (2), and the vampire bats Desmodus and Diaemus (3), were confirmed (Figure 3). A bifurcation in more or less omnivorous bats (Phyllostominae) and predominantly vegetarian species followed. Within the frugivores a sister-group relationship between Rhinophyllinae and the Stenodermatinae was well-supported (BS 99). However, support for a sister-group relationship of Carolliinae and Glyphonycterinae was weak (BS 48). The highly specialized nectarivorous taxa Lonchophylla and Lionycteris do not align closely with other nectarivorous phyllostomids (Glossophaginae). Instead, they shared a common ancestor with the frugivores, as previously proposed by Baker et al. (2003)  and others [33, 49, 50].
Divergence time estimation and model decision
The analysis under the lognormal relaxed clock model (UCLN) produced the smallest confidence intervals compared to the exponential-(UCED) or strict clock model (CLOC). Estimates of mean likelihood, substitution rate, and node age were most accurately inferred under the UCLN model (Table 2). The assumption of the relaxed clock, that branches differ in their substitution rates, was confirmed. A coefficient of variation of 0.405 indicated moderate rate variation . Figure 4 shows the dated Bayesian tree inferred with BEAST under the UCLN model. The common ancestor of all phyllostomids was dated to the Middle Eocene (42 MYA), with a confidence interval between 49- and 37 MYA. Basal lineages within the phyllostomids arose shortly thereafter in the Late Eocene or Early Oligocene (35-32 MYA). The prominent amount of the remaining lineages emerged in a time frame of about ten million years at the transition from Oligocene to Miocene (29-20 MYA), with 21 out of 33 lineages already present in the Early Miocene (20 MYA).
Reconstruction of ancestral states
Figure 3 shows the reconstruction of ancestral states by the maximum-likelihood approach under the Markov k-state model. Only relevant nodes, which will be used in the discussion chapter, are shown. Reconstructed feeding specialization of the common ancestor of all phyllostomids and of the common ancestor of important clades were mapped on the tree (Figure 3). Unambiguous character states were assigned to nodes with a probability of more than 90% for one reconstructed state. The feeding specialization of the common ancestor of all omnivorous and predominantly vegetarian phyllostomid species could not be resolved, as the reconstruction was ambiguous for this node (marked with ?). We obtained probability values of 47% for a nectarivorous-, 39% for an omnivorous-, and 12% for an insectivorous state at this node.
Phylogeny of the Phyllostomidae
Our molecular phylogenetic reconstructions based on more than 10 kb DNA sequences obtained high bootstrap support for almost all nodes and challenges several phylogenetic relationships derived from morphological data sets. Our results partly disagree with recent classifications of phyllostomid bats [45, 49] including: (1) placement of insectivorous genera Macrotus, Micronycteris, Glyphonycteris and Trinycteris within the Phyllostomidae; (2) relationship of the fruit-eating genus Rhinophylla to other frugivores; and (3) relationships among nectarivorous phyllostomids.
The molecular data suggest that the genera Macrotus and Micronycteris do not belong to the subfamily Phyllostominae as proposed by Koopman (1994) , McKenna and Bell (1997) , Wetterer et al. (2000) , and Jones et al. (2002) . Instead, they form two divergent basal lineages within phyllostomid bats (Figure 3). Our data are in line with the findings of Baker et al. (2003) . The authors proposed a classification of two different subfamilies Macrotinae and Micronycterinae. Three studies placed the genera Glyphonycteris and Trinycteris within the subfamily Phyllostominae [5, 45, 49]. In contrast, our data revealed a close relationship of Glyphonycteris and Trinycteris with frugivorous species of the subfamily Carolliinae. Despite low support for this sister-group relationship (BS 45), our data support a closer relationship of Glyphonycteris and Trinycteris to fruit-eating species (BS 99) than to omnivorous phyllostomids of the subfamily Phyllostominae.
The genus Rhinophylla does not belong to the subfamily Carolliinae, as proposed by McKenna and Bell (1997) , Wetterer et al. (2000) , and Jones et al. (2002) . Our data support a sister-group relationship between Rhinophylla and the subfamily Stenodermatinae, as proposed by Baker et al. (2003) .
Many authors excluded the genera Phyllonycteris, Erophylla and Brachyphylla, all endemic to the West Indies, from other nectarivorous phyllostomids and placed them mostly into the subfamilies Phyllonycterinae and Brachyphyllinae [3, 5, 19, 33, 45, 49, 50, 52–59]. In contrast, our data show that these nutritionally more generalized bats belong to the Glossophaginae (BS 100). The three genera are closely related to more specialized nectarivorous bats of the genera Glossophaga, Leptonycteris and Monophyllus. This is in line with an earlier molecular phylogeny of Baker et al. (2003) . The phylogenetic position of highly specialized nectarivorous bats of the genera Lonchophylla and Lionycteris is controversial. Several studies, primarily based on morphological analyses, placed them within the Glossophaginae [5, 46, 47, 52]. Instead, our molecular data suggest that they are closer related to Rhinophyllinae, Stenodermatinae, Carolliinae, and Glyphonycterinae than to the Glossophaginae. This finding is in line with previous studies of Koopman (1994) , Jones et al. (2002)  and Baker et al. (2003) . The distinctness of Lonchophyllinae is also supported by fixed differences in the tongue morphology (see below) between representatives of the Lonchophyllinae and Glossophaginae .
In summary, our study supports the classification of phyllostomid bats after Baker et al. (2003) . Their division into more subfamilies, compared to Koopman (1994)  and Simmons (2005) , seems justifiable, because this better reflects the remarkable ecological diversity of this family.
The vast majority of bats feed on insects . This includes the family Mormoopidae, which represents the sister group of the Phyllostomidae. In addition, the diet of the most basal subfamilies Macrotinae and Micronycterinae consists mainly of insects (Figure 3). These findings indicate, that the common ancestor of phyllostomid bats was an insect-feeder. This supposition is also supported by the maximum-likelihood reconstruction of the ancestral state (Figure 3).
Members of the Phyllostominae have a mixed diet. The reconstruction of the ancestral state for this group revealed that their physiological pre-adaptations to omnivory could have evolved only once, and involved metabolic changes from insectivorous to an omnivorous diet. However, too little is known about the diet of these bats. A high spacial and seasonal plasticity is observed . A few members of the Phyllostominae are carnivorous and feed on small vertebrates . For example, Trachops cirrhosus is specialized on tungara frogs [62, 63]. Such a unique specialization likely evolved in a formerly insectivorous/omnivorous species. It was shown, for the seasonally carnivorous Greater Noctule bat (Nyctalus lasiopterus, Vespertilionidae), that only minor changes are needed to switch from insectivory to carnivory (inclusion of small vertebrates in the diet) [64, 65]. The transition from large-bodied insects to small vertebrates as prey does not need any major adaptations and occurred several times independently in different bats and is correlated with an increase in body size .
A large number of phyllostomid species have a vegetarian diet. They form a monophyletic clade (BS 99), comprising the subfamilies Rhinophyllinae, Stenodermatinae, Carolliinae, Glossophaginae, Lonchophyllinae, and surprisingly the Glyphonycterinae. The last subfamily includes several strict insectivorous species; thus, a shift from a vegetarian diet back to insectivory seems to be the most plausible scenario. Alternatively, the Glyphonycterinae retained the ancestral insectivorous lifestyle. This assumption would require that the frugivorous subfamilies Rhinophyllinae, Stenodermatinae and Carolliinae have evolved their feeding specialization independently from each other. However, the relevant node is weakly supported in our phylogenetic reconstruction. It is also possible (see Figure 2 based on mitochondrial data) that the Glyphonycterinae represents a basal lineage to all frugivores and therefore possess the plesiomorphic state for this group. The common ancestry of all frugivore species was also postulated by previous studies [5, 46, 48, 50, 59]. However, there is a high dietary plasticity in this group. For example, Carollia is a known switch hitter between fruits and insects depending on the season (summarized in ).
The diphyly of the nectarivorous Lonchophyllinae and Glossophaginae is surprising, since they resemble each other in many morphological, behavioural, ecological, and physiological traits (e.g. skull elongation, reduction of dentition, hovering flight, forest foraging behaviour and ability to metabolize pollen). Accordingly, these similarities have evolved independently by natural selection during the adaptation to a nectar-feeding lifestyle. This hypothesis is supported by some obvious differences in these adaptive traits : The lonchophyllines have a deep longitudinal groove in their tongue, lined dorsal and ventral with hairlike papillae. This groove is missing in the glossophagines and hairlike papillae are distributed anterodorsal, forming a brush tip. Furthermore, the lonchophyllines lost most types of papillae found on the tongues of other phyllostmids, including the glossophagines. Also, the internal tongue structure is very different. The lonchophyllines have complex, omnidirectional bundles of muscles within the tongue, while glossophagines have predominantly horizontal skeletal muscle bundles. The complex orientated muscles in the lonchophyllines are supposed to control the shape of the groove during nectar feeding . Drinking behaviour varies widely between both subfamilies (Marco Tschapka, pers. comm., ). Other characters show similar apomorphic states in lonchophyllines and some glossophagines (e.g. posterior shift of sternohyoid origin, xiphoid origin of sternohyoid, elongated hyoglossus and loss of connection to hyoid bone, double insertion of geniohyoid, posterior shift of genioglossus insertion ), however, there are no consistent patterns. The endemic West Indian genera, Brachyphylla, Erophylla, and Phyllonycteris, show many plesiomorphic characters. It seems that functional constraints on the muscular innervation of the tongue curtain the evolutionary signals of these characters. Hence, it is possible that lonchophyllines and glossophagines may have evolved these adaptations for nectar-feeding independently (but see also [67, 68]).
The large number of species within the clade of frugivorous and nectarivorous bats (Figure 2) suggests, that a shift to a vegetarian diet accelerated the diversification rates in this group. The majority of phyllostomid bats, 117 out of 158 listed species , i.e. 74%, belong to this clade. Possibly the presence of numerous vacant ecological niches in tropical and subtropical regions of America (see also ) resulted in allopatric speciation.
Time frame of evolution
Our analysis revealed a time frame of ten million years (29-20 MYA) from Oligocene to Early Miocene, in which all prominent lineages evolved (Figure 4). Most of the species diversification occurred subsequent to the Oligocene epoch (since 23 MYA). During the Miocene substantial changes of the landscape occurred in Tropical America due to massive plate tectonics [70–75]. Global climate cooled and resulted in an increase in aridity [76–78]. Frequent isolation events could have resulted in allopatric populations and thus promoted speciation . Interestingly, the radiation of extant hummingbirds (Trochilidae), another alimentary competitor, shows a similar pattern of diversification in the Middle Miocene . Geologic upheavals as well as the ability of ecologically generalized species to invade new regions were considered as major forces promoting hummingbird radiations in newly arisen montane regions. In order to test whether these factors has also promoted speciation in bats, and to infer other underlying evolutionary mechanisms, a much denser taxon sampling is required.
Our analysis of more than 10.000 base pairs of concatenated DNA sequences reveals a strongly supported phyllostomid phylogeny, thus allowing for clear predictions about the evolution of feeding specialization of these bats. Several morphological and even molecular studies were unable to resolve the specific branches with sufficient support, either due to the convergent nature of the analyzed characters or insufficient amount of sequence data. Our multi-gene approach, combined with a relaxed clock analysis, detected and dated major splitting events within this family. This study gives support for the classification of phyllostomid bats after Baker et al. (2003) . All prominent lineages with diverse feeding strategies evolved within a relatively short time frame of about ten million years from Oligocene to Early Miocene. Geological and climate changes as well as the shift to a vegetarian diet may have promoted the radiation into diverse lineages. In this context, the diphyly of the nectarivorous Lonchophyllinae and Glossophaginae is remarkable. Despite many similarities between both groups, it seems plausible, that they evolved their adaptations to nectarivory independently from each other. This would represent an example of convergent evolution within bats that led to very similar features, which play a major role in food acquisition.
Thirty-seven phyllostomid species of 29 genera were analyzed. Our sampling comprises members of all extant subfamilies [45, 49], except bats of the subfamily Lonchorhinae. We used species and subfamily assignments according to Baker et al. (2003) . One representative each of the families Mormoopidae, Furipteridae, Noctilionidae, Molossidae and Vespertilionidae were used as outgroup taxa. Two closely related specimens were used for the family Rhinolophidae, because we were not able to analyze all loci entirely for one taxon. GenBank accession numbers are given in addtitional file 1. Tissue samples were provided by cooperation partners (see acknowledgements). The name of the body which gave approval and corresponding reference numbers could be obtained from them.
Extraction of total genomic DNA was done by Chloroform-Isoamyl-Phenol precipitation. A 1.3 kb fragment of the exon 28 of the von Willebrand factor gene (vwf) was amplified with the primers vWF-A and vWF-B , or with vWF-A and vWF-B2  within a Nested PCR. Primer vWF-B2 anneals 139 bp upstream from vWF-B. An approximately 1.4 kb fragment of the recombination activating gene 2 (rag2) was amplified with the primers RAG2-F1 and RAG2-R2 , or with RAG2-F1B and RAG2-R2 . The PCR Mastermix (25 μl final reaction volume) included 2 μl of total genomic DNA extract, 1.25 μl of each primer (10 μM), 1 μl of MgCl2 (25 mM), 1 μl of a dNTP-Mix (10 mM) and 1 unit of Peglab Taq polymerase. Nested PCR was performed using 2 μl from a 1:40 delution of the first PCR reaction. The fragments were amplified following a Two-Step protocol. Thermocycling consisted of a 3 min initial denaturation at 95°C, followed by 5 cycles of 30s at 95°C, 50s at 65°C (for the vwf), or 30s at 60°C (for the rag2), and 90s at 72°C. 35 cycles with 50s annealing at 62°C (for the vwf) and 30s at 57°C (for the rag2) were performed, followed by a final extension of 6 min at 72°C. A fragment of exon 11 of the breast cancer susceptibility gene (brca1) was amplified with the primers BRCA1-F126  and a newly designed (ER 515: 5'- AAGTGTTGGAAGCAGGGAAGCTCTTC-3'). The PCR-Mastermix (50 μl final reaction volume) included 2 μl of total DNA extract, 2.5 μl of each primer (10 μM) and 25 μl Phusion Mastermix. Thermocycling consisted of a 30s initial denaturation at 98°C, followed by 5 cycles of 10s of 98°C, 25s at 66°C, and 90s at 72°C. 30 cycles with 25s annealing at 63°C were performed, followed by a final extension of 6 min at 72°C. Two non-coding nuclear loci were also analyzed: 3'-UTR region of the phospholipase C beta 4 gene (plcb4)  and short intron of the phosphoenolpyruvate carboxykinase gene (pepck) ;
We amplified a mitochondrial fragment of the NADH dehydrogenase subunit 1 gene (nd1) and the tRNA Leucin, using the primers ER 65 and ER 66 . Published sequences of five additional mitochondrial loci (COI, Cytb, 12S rRNA, 16S rRNA and tRNA Valin) were incorporated. For all analyses, the ribosomal RNAs and the tRNA Valin were combined (12StRNA16S). Accession numbers are given in Additional file 1. It also includes an overview of all incorporated taxa, loci and sequences and the percentage of missing data per species, as well the geographic origin of our samples. The overall amount of missing data is about 30%.
Alignments and model selection
All alignments were done with Sequencher v4.7  and Bioedit v7.0.9 [87, 88] and checked manually by eye. We performed bootstrap analyses of each individual loci to check for compatibility of their individual phylogenetic signal. Because none of the strongly supported clades based on individual loci were mutually incompatible, we concatenated all loci except the ribosomal RNAs, tRNA Valin and the rag2 gene. These loci were already used by Baker and colleagues to infer a molecular phylogeny of phyllostomid bats . We avoided in a first step the inclusion of them to get an independent data set [Alignment 1]. In a second step we concatenated all mitochondrial loci (this time with the inclusion of the ribosomal RNAs and the tRNA Valin) [Alignment 2] and also all nuclear loci (with rag2) [Alignment 3]. We concatenated all loci into one supermatrix for the final analyses. The supermatrix contained three nuclear protein-coding genes (rag2, vwf and brca1), two non-coding nuclear markers (pepck, plcb4), three mitochondrial protein-coding genes (co1, cytb and nd1), two tRNAs (Valin, Leucin) and two mitochondrial rRNAs (12S, 16S). For the Bayesian analyses, we excluded all 3rd codon positions in the mitochondrial protein-coding genes because they showed a high degree of homoplasy (homoplasy index, HI = 0.75 - parsimony analysis of the 3rd codon positions in PAUP 4.0 beta ). Such high homoplastic characters give a misleading phylogenetic signal and lead especially to an underestimation of real branch lengths. Therefore, we excluded them from the analyses. This resulted in a final length of 10396bp, including 2761 parsimony informative characters [Alignment 4 - Additional file 2]. For the maximum-likelihood analyses, we used a second alignment, in which the mitochondrial protein-coding sequences were translated in amino acids and combined with the remaining DNA sequences [Alignment 5 - Additional file 3]. The best fitting evolutionary model for the protein data was inferred with Prottest v1.4 . The MTMAM model, designed for the evolution of mitochondrial proteins of mammals , showed the highest fit. We ran jModelTest  for the remaining DNA sequences separate for the alignments 1-5. Except for alignment 3, GTR+Γ  was proposed to be the best fitting evolutionary model according to Akaike- (AIC) and Bayesian (BIC) information criterion [94, 95]. The slightly simpler Symmetrical Model SYM+Γ  was proposed for alignment 3 by jModelTest. However, we also used the GTR+Γ model for this data set for general compatibility among the inferences. Genes could have a different sequence evolution. Therefore, we generated five partitioning schemes  for alignment 5 to decide, which is the best adjustment for our analysis: (1) no partitioning; (2) mitochondrial- and nuclear loci separately; (3) three partitions; (4) eight partitions; and (5) 14 partitions with partitioning into codon positions for all nuclear genes. According to AIC and BIC, scheme 5 was preferred.
Equal weighted maximum-parsimony (MP) analyses were performed with PAUP 4.0 beta  with a heuristic search using the TBR (tree-bisection-reconnection) algorithm for branch swapping. Bootstrap inferences were conducted separately for each loci with 500 pseudoreplicates.
Maximum-likelihood (ML) inferences were performed with RAxML v7.0.4 [98–100]. ML searches were conducted with the rapid hill-climbing algorithm [101, 102] under GTR+Γ with four rate categories as model of evolution. Multiple independent runs were started to get an impression of the robustness of the phylogenetic reconstruction. Support values were obtained through a full non-parametric bootstrap- or rapid bootstrap inference (stated for each analysis).
Reconstruction of ancestral states
Ancestral character states were reconstructed in Mesquite v2.71 . Observed character states (insectivore, sanguinivore, omnivore, frugivore and nectarivore) of the main diet were mapped on the original maximum-likelihood tree (Figure 3). We used the "Trace Character History" analysis with a symmetric, one-parameter Markov k-state model [104, 105], which computes likelihoods for categorical characters, and reconstructs ancestral states by the maximum marginal probability (MLE) criterion.
Bayesian inferences were performed with BEAST v1.4.8 . The searches were conducted under Hasegawa-Kishino-Yano HKY+Γ  with four rate categories as model of evolution. We chose a simpler model of sequence evolution for the Bayesian analyses as proposed by jModelTest because there was a trade-off between computational power and model complexity. It was not possible to get a consistent phylogenetic reconstruction between different runs under the GTR+Γ model in reasonable time. Bayes factor analysis between these runs under the GTR+Γ model resulted always in values far above 20 and sampling efficiency was drastically reduced compared to the HKY+Γ model. A high Bayes factor is a sign for incompatibility and poor convergence among the trees gathered in independent runs.
Calibration of the molecular clock
We incorporated three different calibration points including: (1) divergence between Vespertilionidae and Molossidae set at 37 million years ago (MYA) in the Middle Eocene ; (2) age of the Mormoopidae oldest fossils from Whitneyan (30-32 MYA) land deposits in Florida ; and (3) age of the oldest crown group fossils of the phyllostomids in the Laventan about 11.8 to 13.8 MYA  and age of the oldest stem group fossils in the Whitneyan within the Early Oligocene . We used the proposed age of the fossils and lognormal distributions to model minimum age constraints for the specific nodes (1,2). Maxmimum age constraints were set to the Cretaceous-Tertiary boundary at 65 MYA (1,2). Additional, a maximum age constraint for the phyllostomids (3) was set with an exponential distribution to 34 MYA with an arbitrarily lower limit of 11.5 MYA.
Model- and prior decision
We performed several Bayesian inferences under one strict (CLOC) and two relaxed (UCLN, UCED) molecular clock models [111, 112]. Always 10 million steps were performed. We examined the joint influence of the calibrations on the divergence time estimates by running a strict clock model with fixed topology, but with no sequence data. Further, we examined the influence of each individual calibration by running several inferences under an uncorrelated lognormal relaxed clock model (UCLN) with all possible combinations of the three calibrations. A precise examination and comparison of the results were performed in Tracer v1.4 . An overview of important parameters for model comparison is given in Table 2. Important parameters, such as mean likelihood value, substitution rate, and node age, were calculated for every inference and compared with each other. Confidence intervals measured as 95% highest posterior density interval (HPD) were also computed. The clock model that produced the smallest confidence intervals altogether was considered most appropriate for the data .
Estimation of divergence times
We conducted three independent runs for the final divergence time estimates under the UCLN model with 20 million inferences and a sample frequency of 1000 steps. We used always the same parsimony tree as starting point. We compared the results and calculated pairwise Bayes factors for the difference in their marginal likelihoods. The first 4 million steps were cut off as burnin for each comparison. Low Bayes factors are a sign for high convergence of the values and compatibility of the inferences, while high Bayes factors indicate incompatibility. Individual runs were combined with LogCombiner, TreeAnnotator and analyzed with Tracer v1.4 and FigTree v1.1.2 . TreeAnnotator and LogCombiner are provided as part of the BEAST package.
Gardner AL: Chromosomal variation in Vampyressa and a review of chromosomal evolution in the Phyllostomidae (Chiroptera). Syst Zool. 1977, 26: 300-318. 10.2307/2412677.
Findley JS: Bats: a community perspective. 1993, Cambridge, Cambridge Univ Press
Gimenez EDA, Ferrarezzi H, Taddei VA: Lingual morphology and cladistic analysis of the New World nectar-feeding bats (Chiroptera: Phyllostomidae). J Comp Biol. 1996, 1: 41-64.
Freeman PW: Macroevolution in microchiroptera: recoupling morphology and ecology with phylogeny. Evol Ecol Res. 2000, 2: 317-335.
Wetterer AL, Rockman MV, Simmons NB: Phylogeny of phyllostomid bats (Mammalia: Chiroptera): Data from diverse morphological systems, sex chromosomes, and restriction sites. Bull Amer Mus Nat Hist. 2000, 248: 1-200. 10.1206/0003-0090(2000)248<0001:POPBMC>2.0.CO;2.
Giannini NP, Kalko KV: Trophic structure in a large assemblage of phyllostomid bats in Panama. Oikos. 2004, 105: 209-220. 10.1111/j.0030-1299.2004.12690.x.
Rex K, Kelm DH, Wiesner K, Matt FG, Kunz TH, Voigt CC: Species richness and structure of three Neotropical bat assemblages. Biol J Linn Soc. 2008, 94: 617-629. 10.1111/j.1095-8312.2008.01014.x.
Humphrey SR, Bonaccorso FJ, Zinn TL: Guild structure of surface-gleaning bats in Panama. Ecology. 1983, 64: 284-294. 10.2307/1937076.
Dobat K, Peikert-Holle T: Blüten und Fledermäuse: Bestäubung durch Fledermäuse und Flughunde (Chiropterophilie). 1985, Frankfurt am Main, Germany: Waldemar Kramer
Fleming TH: The short-tailed fruit bat: A study in plant-animal interactions. 1988, Chicago Univ Press, 1-380. [Wildlife Behavior and Ecology series].
Freeman PW: Frugivorous and animalivorous bats (Microchiroptera): Dental and cranial adaptations. Biol J Linn Soc. 1988, 33: 249-272. 10.1111/j.1095-8312.1988.tb00811.x.
Medellin RA: Prey of Chrotopterus auritus, with notes on feeding behavior. J Mammal. 1988, 69: 841-844. 10.2307/1381644.
Handley COJ, Wilson DE, Gardner AL: Demography and natural history of the common fruit bat, Artibeus jamaicensis, on Barro Colorado Island, Panama. Smithson Contrib Zool. 1991, 511: 1-173.
von Helversen O: Adaptations of flowers to the pollination by glossophagine bats. Animal-plant interactions in tropical environments. Edited by: Barthlott W, Naumann CM, Schuchmann SLK, Schuchmann KL. 1993, Bonn, Germany, Museum Koenig, 41-59.
Freeman PW: Nectarivorous feeding mechanisms in bats. Biol J Linn Soc. 1995, 56: 439-463. 10.1111/j.1095-8312.1995.tb01104.x.
Proctor M, Yeo P, Lack A: The natural history of pollination. 1996, Portland, Oregon, Timber Press
Vogel S: Chiropterophilie in der neotropischen Flora. Neue Mitteilungen III. Flora, B. 1969, 158: 289-323.
Phillips CJ: The dentition of glossophagine bats: development, morphological characteristics, variation, pathology, and evolution. Miscellaneous Publications of the Museum of Natural History. 1971, University of Kansas, 54: 1-138.
Smith JD: Chiropteran evolution. Biology of the Bats of New World Phyllostomidae. Part I. Edited by: Baker RJ, Jones JK, Carter DC. 1976, Spec Publ Mus Tex Tech Univ, 49-69.
Baker RJ, Bass RA: Evolutionary relationship of the Brachyphyllinae to the glossophagine genera Glossophaga and Monophyllus. J Mammal. 1979, 60: 364-372. 10.2307/1379808.
Hopkins HC: Floral biology and pollination ecology of the Neotropical species of Parkia. Ecology. 1984, 72: 1-23. 10.2307/2260003.
Eguiarte L, Búrquez A: Reproductive ecology of Manfreda brachystachya, an iteroparous species of Agavaceae. Southwest Nat. 1987, 32: 169-178. 10.2307/3671560.
Gribel R, Hay JD: Pollination ecology of Caryocar brasiliense (Caryocaraceae) in Central Brazil cerrado vegetation. J Trop Ecol. 1993, 9: 199-211. 10.1017/S0266467400007173.
von Helversen O: Blumenfledermäuse und Fledermausblumen - Wechselbeziehungen zwischen Blüte und Bestäuber und energetische Grenzbedingungen. Rundgespräche der Kommission für Ökologie. Band 10, Tropenforschung. 1995, Bayr Akad Wiss, 217-229.
Luckow M, Hopkins HCF: A cladistic analysis of Parkia (Leguminosae: Mimosoideae). Am J Bot. 1995, 82: 1300-1320. 10.2307/2446253.
Winter Y, Voigt CC, von Helversen O: Gas exchange during hovering flight in nectar-feeding bat Glossophaga soricina. J Exp Biol. 1998, 201: 237-244.
Dumont ER: The effect of food hardness on feeding behaviour in frugivorous bats (Phyllostomidae): an experimental study. J Zool. 1999, 248: 219-229. 10.1111/j.1469-7998.1999.tb01198.x.
Carstens BC, Lundrigan BL, Myers P: A Phylogeny of the Neotropical nectar-feeding bats (Chiroptera: Phyllostomidae) based on morphological and molecular data. J Mamm Evol. 2002, 9: 23-39. 10.1023/A:1021331711108.
von Helversen O, Winter Y: Glossophagine bats and their flowers: costs and benefits for plants and pollinators. Bat Ecology. Edited by: Kunz TH, Fenton MB. Chicago and London, Chicago Univ Press, 2003: 346-397. 2
Winter Y, von Helversen O: Operational tongue length in phyllostomid nectar-feeding bats. J Mammal. 2003, 84: 886-896. 10.1644/BWG-032.
Kelm DH, von Helversen O: How to budget metabolic energy - torpor in a small Neotropical mammal. J Comp Physiol B. 2007, 177 (6): 667-677. 10.1007/s00360-007-0164-5.
Voigt CC, Speakman JR: Nectar-feeding bats fuel their high metabolism directly with exogenous carbohydrates. Funct Ecol. 2007, 21: 913-921. 10.1111/j.1365-2435.2007.01321.x.
Griffiths TA: Systematics of the New World nectar-feeding bats (Mammalia: Phyllostomidae) based on the morphology of the hyoid and lingual regions. Am Mus Novit. 1982, 27 (42): 1-45.
Heithaus ER, Stashko E, Anderson PK: Cumulative effects of plant- animal interactions on seed production by Bauhinia ungulata, a Neotropical legume. Ecology. 1982, 63: 1294-1302. 10.2307/1938857.
Winter Y: Energetic cost of hovering flight in a nectar-feeding bat measured with fast-response respirometry. J Comp Physiol B. 1998, 168: 434-444. 10.1007/s003600050163.
Norberg UM, Kunz TH, Steffensen JF, Winter Y, von Helversen O: The cost of hovering and forward flight in a nectar-feeding bat, Glossophaga soricina, estimated from aerodynamic theory. J Exp Biol. 1993, 182: 207-227.
Hedenström A, Johansson LC, Wolf M, Busse R, Winter Y, Spedding GR: Bat flight generates complex aerodynamic tracks. Science. 2007, 316 (5826): 894-897. 10.1126/science.1142281.
Howell D: Bats and pollen: physiological aspects of the syndrome of chiropterophily. Comp Biochem Physiol. 1974, 48: 263-276. 10.1016/0300-9629(74)90707-5.
Stanley RG, Linskens HG: Pollen biology, biochemistry, management. 1974, Berlin, Heidelberg, New York, Springer-Verlag
Law BS: Physiological factors affecting pollen use by Queensland blossom bats, Syconycteris australis. Funct Ecol. 1992, 65: 634-648.
Fleming TH: The use of stable isotopes to study the diets of plant-visiting bats. Bats: Ecology, Behavior, and Evolution. Edited by: Racey PA, McDonnell U, Swift S. 1995, Oxford, Oxford Univ Press, 99-110.
Herrera LGM, del Río CM: Pollen digestion by New World Bats: Effects of processing time and feeding habits. Ecology. 1998, 79: 2828-2838. 10.1890/0012-9658(1998)079[2828:PDBNWB]2.0.CO;2.
Schondube JE, Herrera LGM, del Río CM: Diet and the evolution of digestion and renal function in phyllostomid bats. Zoology. 2001, 104: 59-73. 10.1078/0944-2006-00007.
Mirón LLM, Herrera LGM, Ramírez PN, Hobson KA: Effect of diet quality on carbon and nitrogen turnover and isotopic discrimination in blood of a New World nectarivorous bat. J Exp Biol. 2006, 209: 541-548. 10.1242/jeb.02016.
Simmons NB: Order Chiroptera. Mammalian Species of the World: A Taxonomic and Geographic Reference. Edited by: Wilson DE, Reeder DM, Baltimore. 2005, Maryland: Johns Hopkins Univ Press, 312-529. 3
Baker RJ, Hood CS, Honeycutt RL: Phylogenetic relationship and classification of the higher categories of the New World bat family Phyllostomidae. Syst Zool. 1989, 38: 228-238. 10.2307/2992284.
McKenna MC, Bell SK: Classification of mammals above the species level. 1997, New York, NY, Columbia Univ Press
Baker RJ, Hoofer SR, Porter CA, Van Den Bussche RA: Diversification among New World leaf-nosed bats: An evolutionary hypothesis and classification inferred from digenomic congruence of DNA sequence. Occ Pap Mus Tex Tech Univ. 2003, 230: 1-32.
Koopman KF: Chiroptera: Systematics. Handbook of Zoology, Vol 8, Part 60. 1994, Berlin, Germany: Walter de Gruyter
Jones KE, Purvis A, MacLarnon A, Bininda-Emonds ORP, Simmons NB: A phylogenetic supertree of the bats (Mammalia: Chiroptera). Biol Rev. 2002, 77: 223-259. 10.1017/S1464793101005899.
Drummond AJ, Ho SYW, Rawlence N, Rambaut A: A rough guide to BEAST 1.4. 2007
Miller GSJ: The families and genera of bats. Bull US Nat Mus. 1907, 57: 1-282.
Simpson GG: The principles of classification and a classification of the mammals. Bull Am Mus Nat Hist. 1945, 85: 1-350.
de la Torre L: The evolution, variation, and systematics of the Neotropical bats of the genus Sturnira. PhD thesis. 1961, Univ Illinois, Urbana
Silva-Taboda GS, Pine RH: Morphological and behavioral evidence for the relationship between the bats genus Brachyphylla and the Phylonycterinae. Biotropica. 1969, 1: 10-19. 10.2307/2989745.
Koopman KF, Jones JKJ: A classification of bats. About Bats. Edited by: Slaughter BH, Walton WD, Dallas, TX. 1970, South Meth Univ Press, 22-28.
Corbet GB, Hill JE: A world list of mammalian species. 1980, London: British Museum, 1
Baker RJ, Honeycutt RL, Arnold ML, Sarich VM, Genoways AA: Electrophoretic and immunological studies of the relationship of the Brachyphyllinae and the Glossophaginae. J Mammal. 1981, 62: 665-672. 10.2307/1380587.
Baker RJ, Porter CA, Patton JC, Van Den Bussche RA: Systematics of the family Phyllostomidae based on RAG2 DNA sequences. Occ Pap Mus Tex Tech Univ. 2000, 202: 1-16.
Rex K, Czaczkes BI, Michener R, Kunz TH, Voigt CC: Specialization and omnivory in diverse mammalian assemblages. Ecoscience. 2010, 17 (1): 37-46. 10.2980/17-1-3294.
Norberg UM, Fenton MB: Carnivorous bats?. Biol J Linn Soc. 1988, 33: 383-394. 10.1111/j.1095-8312.1988.tb00451.x.
Barclay RMR, Fenton MB, Tuttle MD, Ryan MJ: Echolocation calls produced by Trachops cirrhosus while hunting frogs. Can J Zool. 1981, 59: 750-753. 10.1139/z81-107.
Tuttle MD, Ryan MJ: Bat predation and the evolution of frog vocalizations in the Neotropics. Science. 1981, 214: 677-678. 10.1126/science.214.4521.677.
Ibáñez C, Juste J, García-Mudarra JL, Agirre-Mendi PT: Bat predation on nocturnally migrating birds. PNAS. 2001, 98: 9700-9702. 10.1073/pnas.171140598.
Popa-Lisseanu AG, Delgado-Huertas A, Forero MG, Rodríguez A, Arlettaz R, Ibáñez C: Bats' conquest of a formidable foraging niche: The myriads of nocturnally migrating songbirds. PLoS ONE. 2007, 2 (2): e205-10.1371/journal.pone.0000205.
Cloutier D, Thomas DW: Carollia perspicillata. Mammal species. 1992, 417: 1-9.
Haiduk MW, Baker RJ: Scientific method, opinion, phylogenetic reconstruction and nectar- feeding bats: A response to Griffiths and Warner. Syst Zool. 1984, 33: 343-350. 10.2307/2413082.
Smith JD, Hood CS: Genealogy of the New World nectar-feeding bats revisited: a reply to Griffiths. Syst Zool. 1984, 33: 435-460. 10.2307/2413096.
Jones KE, Bininda-Emonds ORP, Gittleman JL: Bats, clocks, and rocks: Diversification patterns in Chiroptera. Evolution. 2005, 59 (10): 2243-2255.
Allmon WD, Bottjer DJ: Evolutionary Paleoecology: The Ecological Context of Macroevolutionary Change. 2001, New York, Columbia Univ Press
Rogers RD, Karason H, Van Der Hilst RD: Epeirogenic uplift above a detached slab in northern Central America. Geology. 2002, 30 (11): 1031-1034. 10.1130/0091-7613(2002)030<1031:EUAADS>2.0.CO;2.
MacFadden BJ, Higgins P: Ancient ecology of 15-million-year-old browsing mammals within C3 plant communities from Panama. Oecologia. 2004, 140: 169-182. 10.1007/s00442-004-1571-x.
MacMillan I, Gans PB, Alvarado G: Middle Miocene to present plate tectonic history of the southern Central American Volcanic Arc. Tectonophysics. 2004, 392 (1-4): 325-348. 10.1016/j.tecto.2004.04.014.
Kirby MX, Jones DS, MacFadden BJ: Lower Miocene stratigraphy along the Panama Canal and its bearing on the Central American peninsula. PloS ONE. 2008, 3 (7): e2791-10.1371/journal.pone.0002791.
Meschede M, Barckhausen U, Engels M, Weinrebe W: The trace of the Pacific-Cocos-Nazca triple junction in the Central Pacific and the formation of an overlapping spreading centre. Terra Nova. 2008, 20 (3): 246-251. 10.1111/j.1365-3121.2008.00812.x.
Shevenell AE, Kennett JP, Lea DW: Middle Miocene southern ocean cooling and Antarctic cryosphere expansion. Science. 2004, 17: 1766-1770. 10.1126/science.1100061.
Singh VP, Harmon RS: The Río Chagres, Panama: A multidisciplinary profile of a Tropical watershed. Water Science and Technology Library. 2005, Netherlands: Springer, 45-62.
Kürschner WM, Kvacek Z, Dilcher DL: The impact of Miocene atmospheric carbon dioxide fluctuations on climate and the evolution of terrestrial ecosystems. PNAS. 2008, 105 (2): 449-453. 10.1073/pnas.0708588105.
Coyne JA, Orr HA: Speciation. 2004, Sunderland, Massachusetts USA, Sinauer Associates, Inc
Bleiweiss R: Tempo and mode of hummingbird evolution. Biol J Linn Soc. 2008, 65: 63-76. 10.1111/j.1095-8312.1998.tb00351.x.
Porter CA, Goodman M, Stanhope MJ: Evidence on mammalian phylogeny from sequences of exon 28 of the von Willebrand factor gene. Mol Phylogenet Evol. 1996, 5: 89-101. 10.1006/mpev.1996.0008.
Teeling CE, Scally M, Kao DJ, Romagnoli ML, Springer MS, Stanhope MJ: Molecular evidence regarding the origin of echolocation and flight in bats. Nature. 2000, 403 (13): 188-192.
Murphy WJ, Eizirik E, Johnson WE, Zhang JP, Ryder OA, OBrien SJ: Molecular phylogenetics and the origins of placental mammals. Nature. 2001, 409: 614-618. 10.1038/35054550.
Piaggio AJ, Perkins SL: Molecular phylogeny of North American long-eared bats (Vespertilionidae: Corynorhinus); inter- and intraspecific relationships inferred from mitochondrial and nuclear DNA sequences. Mol Phylogenet Evol. 2005, 37 (3): 762-775. 10.1016/j.ympev.2005.03.029.
Mayer F, von Helversen O: Cryptic diversity in European bats. Proc R Soc Lond. 2001, 268B: 1825-1832. 10.1098/rspb.2001.1744.
SEQUENCHER 4.7 DNA sequence assembly software: Gene Codes Corporation. [http://www.genecodes.com/]
Thompson JD, Higgins DG, Gibson TJ: CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, positions-specific gap penalties and weight matrix choice. Nucl Acids Res. 1994, 22: 4673-4680. 10.1093/nar/22.22.4673.
Hall TA: BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucl Acids Symp Ser. 1999, 41: 95-98.
PAUP*: Phylogenetic Analysis Using Parsimony (and other methods) 4.0 Beta. [http://www.sinauer.com/detail.php?id=8060]
Abascal F, Zardoya R, Posada D: ProtTest: Selection of best-fit models of protein evolution. Bioinformatics. 2005, 21 (9): 2104-2105. 10.1093/bioinformatics/bti263.
Yang Z: Synonymous and nonsynonymous rate variation in nuclear genes of mammals. Mol Evol. 1998, 46: 409-418. 10.1007/PL00006320.
Posada D: jModelTest: Phylogenetic model averaging. Mol Biol Evol. 2008, 25 (7): 1253-1256. 10.1093/molbev/msn083.
Rodriguez F, Oliver JF, Marin A, Medina JR: The general stochastic model of nucleotide substitutions. J Theor Biol. 1990, 142: 485-501. 10.1016/S0022-5193(05)80104-3.
Akaike H: A new look at the statistical model identification. IEEE Trans Aut Cont. 1974, 19 (6): 716-723. 10.1109/TAC.1974.1100705.
Schwarz GE: Estimating the dimension of a model. Ann Stat. 1978, 6 (2): 461-464. 10.1214/aos/1176344136.
Zharkikh A: Estimation of evolutionary distances between nucleotide sequences. J Mol Evol. 1994, 39: 315-329. 10.1007/BF00160155.
McGuire JA, Witt CC, Altshuler DL, Remsen JV: Phylogenetic systematics and biogeography of hummingbirds: Bayesian and maximum likelihood analyses of partitioned data and selection of an appropriate partitioning strategy. Syst Biol. 2007, 56: 837-856. 10.1080/10635150701656360.
Stamatakis A, Ludwig T, Meier H: A program for sequential, parallel distributed inference of large phylogenetic trees. Concurrency and Computation: Pract Exper. 2005, 17 (14): 1705-1723. 10.1002/cpe.954.
Stamatakis A: RAxML-VI-HPC: Maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics. 2006, 22 (21): 2688-2690. 10.1093/bioinformatics/btl446.
Stamatakis A, Hoover P, Rougemont J: A rapid bootstrap algorithm for the RAxML web servers. Syst Biol. 2008, 57: 758-771. 10.1080/10635150802429642.
Stamatakis A, Ludwig T, Meier H: RAxML-III: a fast program for maximum likelihood-based inference of large phylogenetic trees. Bioinformatics. 2005, 21 (4): 456-463. 10.1093/bioinformatics/bti191.
Stamatakis A, Ott M, Ludwig T: RAxML-OMP: An efficient program for phylogenetic inference on SMPs. Proceedings of 8th International Conference on Parallel Computing Technologies (PaCT2005). 2005, 3606: 288-302.
Maddison WP, Maddison DR: Mesquite: a modular system for evolutionary analysis. 2009
Cunningham CW, Omland KE, Oakley TH: Reconstructing ancestral character states: A critical reappraisal. Tree. 1998, 13 (9): 361-366.
Lewis PO: A likelihood approach to estimating phylogeny from discrete morphological character data. Syst Biol. 2001, 50: 913-925. 10.1080/106351501753462876.
Drummond AJ, Rambaut A: BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol Biol. 2007, 7: 214-10.1186/1471-2148-7-214.
Hasegawa M, Kishino H, Yano T: Dating of the human-ape splitting by a molecular clock of mitochondrial DNA. Mol Evol. 1985, 22: 160-174. 10.1007/BF02101694.
Morgan GS: New bats in the Neotropical families Emballonuridae and Mormoopidae from the Oligocene and Miocene of Florida, and the biochronology of Florida Whitneyan, Arikareean, and Hemingfordian faunas. J Vert Paleont. 2002, 22 (Suppl 3): 90A-
Czaplewski NJ, Takai M, Naeher TM, Shigehara N, Setoguchi T: Additional bats from the Middle Miocene La Venta fauna of Colombia. Rev Acad Colombiana de Cienc Exactas, Fis Nat. 2003, 27: 263-282.
Czaplewski NJ, Morgan GS, McLeod SA: Evolution of Tertiary mammals of North America, Vol 2, small mammals, Xenarthrans, and marine mammals. 2008, Cambridge, Cambridge Univ Press
Drummond AJ, Nicholls GK, Rodrigo AG, Solomon W: Estimating mutation parameters, population history and genealogy simultaneously from temporally spaced sequence data. Genetics. 2002, 161: 1307-1320.
Drummond AJ, Ho SYW, Phillips MJ, Rambaut A: Relaxed phylogenetics and dating with confidence. PLoS Biology. 2006, 4 (5): e88-10.1371/journal.pbio.0040088.
Molecular evolution, phylogenetics and epidemiology. [http://tree.bio.ed.ac.uk/software/tracer/]
Molecular evolution, phylogenetics and epidemiology. [http://tree.bio.ed.ac.uk/software/figtree/]
Ferrarezzi H, Gimenez EDA: Systematic patterns and the evolution of feeding habits in Chiroptera (Archonta: Mammalia). J Comp Biol. 1996, 1: 75-94.
Dumont ER: Feeding mechanisms in bats: Variation within the constraints of flight. Integr Comp Biol. 2007, 47: 137-146. 10.1093/icb/icm007.
The Sorenson lab Boston University. [http://people.bu.edu/msoren/TreeRot.html]
Pons J, Barraclough TG, Gomez-Zurita J, Cardoso A, Duran DP, Hazell S, Kamoun S, Sumlin WD, Vogler AP: Sequence-based species delimitation for the DNA taxonomy of undescribed insects. Syst Biol. 2006, 55: 595-609. 10.1080/10635150600852011.
Fontaneto DE, Herniou C, Boschetti M, Caprioli G, Melone C, Ricci C, Barraclough TG: Independently evolving species in asexual bdelloid rotifers. PLoS Biology. 2007, 5: e87-10.1371/journal.pbio.0050087.
Monaghan MT, Wild R, Elliot M, Fujisawa T, Balke M, Inward DJG, Lees DC, Ranaivosolo R, Eggleton P, Barraclough TG, Vogler AP: Accelerated species inventory on Madagascar using coalescent-based models of species delineation. Syst Biol. 2009, 58 (3): 298-311. 10.1093/sysbio/syp027.
We thank Marco Tschapka, Martina Nagy, Mirjam Knörnschild, and Ralph Simon for fruitful discussions and for providing tissue samples. Wolfram Schulze, Jana Ustinova, Dagmar Dachlauer, Andrea Ross, Claudius Kerth, Corinna Koch von Helversen, Sebastian Ebert, Heiko Stuckas, Kai Drilling, and Dirk Berger supported the study at different stages. We would like to thank Christian Voigt, Simon Ghanem, Markus Zweier, Inge Müller, Corine Schoebel, Ulrich Marckmann, Volker Runkel, Dina Dechmann, Kirsten Jung, Stefan Prost, Sebastian Lippold, Michael Knapp, and Maria Helbig for their willingness to discuss diverse aspects and giving advice. We thank the CBSU Web Computing Resources (BIOHPC) for their free web-based service for phylogenetic analysis with BEAST. Fruitful discussions occurred with Alexandros Stamatakis, Ignacio Gonzalo Bravo, and Olaf Bininda-Emonds during the Carolinensiel summer school in 2008. Finally we thank two anonymous reviewers for their excellent comments and help to improve this manuscript. The study was financed by the Luise-Prell-Stiftung and the Schmauser-Stiftung at the University of Erlangen- Nürnberg.
OvH was the initiator of this study. FM supervised the whole project, gave many ideas, helped to evaluate the results and to draw up the manuscript. TD made all the lab work, performed the phylogenetic analyses and wrote the manuscript.