- Research article
- Open access
- Published:
The evolutionary history and ancestral biogeographic range estimation of old-world Rhinolophidae and Hipposideridae (Chiroptera)
BMC Ecology and Evolution volume 22, Article number: 112 (2022)
Abstract
Background
Family Rhinolophidae (horseshoe bats), Hipposideridae (leaf-nosed bats) and Rhinonycteridae (trident bats) are exclusively distributed in the Old-World, and their biogeography reflects the complex historic geological events throughout the Cenozoic. Here we investigated the origin of these families and unravel the conflicting family origin theories using a high resolution tree covering taxa from each zoogeographic realm from Africa to Australia. Ancestral range estimations were performed using a probabilistic approach implemented in BioGeoBEARS with subset analysis per biogeographic range [Old-World as whole, Australia–Oriental–Oceania (AOO) and Afrotropical–Madagascar–Palearctic (AMP)].
Result
Our result supports an Oriental origin for Rhinolophidae, whereas Hipposideridae originated from the Oriental and African regions in concordance with fossil evidence of both families. The fossil evidence indicates that Hipposideridae has diversified across Eurasia and the Afro-Arabian region since the Middle Eocene. Meanwhile, Rhinonycteridae (the sister family of Hipposideridae) appears to have originated from the Africa region splitting from the common ancestor with Hipposideridae in Africa. Indomalaya is the center of origin of Rhinolophidae AOO lineages, and Indomalayan + Philippines appears to be center of origin of Hipposideridae AOO lineage indicating allopatric speciation and may have involved jump-dispersal (founder-event) speciation within AOO lineage. Wallacea and the Philippines may have been used as stepping stones for dispersal towards Oceania and Australia from the Oriental region. Multiple colonization events via different routes may have occurred in the Philippines (i.e., Palawan and Wallacea) since the Late Miocene. The colonization of Rhinolophidae towards Africa from Asia coincided with the estimated time of Tethys Ocean closure around the Oligocene to Miocene (around 27 Ma), allowing species to disperse via the Arabian Peninsula. Additionally, the number of potential cryptic species in Rhinolophidae in Southeast Asia may have increased since Plio-Pleistocene and late Miocene.
Conclusion
Overall, we conclude an Oriental origin for Rhinolophidae, and Oriental + African for Hipposideridae. The result demonstrates that complex historical events, in addition to species specific ecomorphology and specialization of ecological niches may shape current distributions.
Background
Bats constitute the second most diverse group of mammals after rodents, with over 1400 species recognized to date [1]. The rapid diversification of this order may reflect their capacity for powered flight and echolocation, which has allowed them to colonize a wide range of ecological niches [2, 3], and is responsible for bats frequently represent the only native mammals on oceanic islands. In the Old World, insectivorous bats are split into those within the Yangochiroptera (of which most families are not limited to the Old World), and the Old-World endemic Yinpterochiroptera, which includes Rhinolophoidea as well as the frugivorous Pteropodidae, (with Old-World insectivorous bat communities frequently dominated by Rhinolophoidea). Here, we focus exclusively on Old-World insectivorous bats, the superfamily Rhinolophoidea and to date, it encompasses six families of insectivorous bats (Crasseonycteridae, Megadermatidae, Rhinopomatidae, Rhinolophidae, Hipposideridae and Rhinonycteridae [2, 4], but in this study we only focus on the three latter families due to their high diversity (Crasseonycteridae only has one described species, and Megadermatidae only six species). Rhinolophidae, Hipposideridae and Rhinonycteridae are monophyletic and are closely related, with Hipposideridae as sister family of Rhinolophidae and Rhinonycteridae recently resolved as separate family from Hipposideridae.
Rhinolophidae (horseshoe bats) consist of a single genus Rhinolophus Lacépède, 1799. They are insectivorous bats distributed throughout the Old World, primarily in tropical regions, from Africa through Eurasia, Oceania and Australia [5]. Around 106 species have been described to date [1], although this number is likely to be an underestimate given that many species are cryptic [6]. The common name of Rhinolophidae is derived from specialized horseshoe-shaped noseleaf which is used to emit acoustic calls emission [7,8,9,10,11,12,13]. The distinctive noseleaf morphological components, include the sella, lancet, furrows of lancet, internarial cup and ears shape of each rhinolophids species plays important role in determining call structure, and may provide useful cues in the identification of cryptic species [6].
Hipposideridae is the sister family of Rhinolophidae. Commonly known as leaf-nosed bats, they are distributed in the same range as Rhinolophidae across the Old World [14]. They consist of seven genera: Hipposideros, Anthops, Asellia, Aselliscus, Coelops, Doryrhina, Macronycteris, including 90 described species [1]. Recently Cloeotis, Paratriaenops, Rhinonycteris, and Triaenops were elevated to a separate family, Rhinonycteridae (the Trident bats) [15], and is sister family to Hipposideridae [14]. Hipposideros is the most diverse genus within Hipposideridae, which includes almost 80% of the total species in the family, and many species are cryptic [13, 16, 17], and thus true diversity may be considerably higher [1].
The systematics of these groups have not been well resolved. Although recent years have seen significant progress towards resolving the systematics of bats using an abundance of phylogenetic datasets [4], these analyses generally had low systematic coverage and have only included a limited number of species within each family, and with limited spatial coverage [14, 18]. As a result, the understanding of systematics and evolution in these predominantly tropical groups are limited by the number of genes in the study, and systematic biases from limited taxon sampling and low geographical coverage [14, 18,19,20,21,22,23,24,25]. As the consequence of contradictory dating estimations, the interfamilial time divergence estimation within Rhinolophidae and Hipposideridae still has not been reliably estimated, thus assessing the evolutionary history and the biogeographic origin is still challenging (i.e., [14, 18, 22,23,24,25].
Accurate inferrence of the biogeographic history of a taxon is heavily dependant on estimated ages and relationship between taxa, but often lacks rigour when it is presented as a narrative in addition to phylogenetic studies without explicit analysis [26,27,28]. The geographic origin of ancestral Rhinolophidae remains unresolved due to the limited fossil record, lack of representative taxonomic sampling across their distribution [14], and conflicting phylogenies [20].
The biogeographic history of bats is particularly interest due to their capability of true flight and ability to disperse over a wide geographic area. Several alternative hypotheses of the Rhinolophidae family origins have been proposed based on different data types and methods, for instance European origin [29], Asian origin [18, 23, 30], African origin [14, 22], and Middle-Eastern origin [20]. Flight ability relates to wing shape for each species and thus has direct implications for the probability of these routes (as reviewed in Norberg and Rayner 1987) [2, 31,32,33]. Bats with pointed wings and a high aspect ratio are expected to have energetically inexpensive flight and be able to travel long distances but have low manouverability. Wings with high wing loading enables fast flight to reduce the time invested to achieve long distance flight, thus suitable for migration and dispersal over long-distances [34]. Conversely, broad-short wings with rounded tips (low wing loading, low aspect ratio) which are present through most rhinolophids species may limit dispersal of species across large water bodies or other open areas [35] and current distribution of taxa may have follow vicariance events in the past. Furthermore, different methods have different historical biogeographic event assumptions, which has a significant influence on the inferences. Thus, to distinguish between these competing hypothesis, testing and comparison between statistical approaches is needed to optimize model selection and conduct meaningful and well-supported biogeographic analyses [36].
In this study, we use probabilistic modeling of geographic range evolution, which allows users to statistically choose the number of biogeographic models based on Maximum Likelihood and Bayesian methods implemented in the R package “BioGeoBEARS” [36,37,38,39]. However, to date, no study has performed this analysis in Rhinolophidae and Hipposideridae. This approach unifies multiple models and provides a flexible framework for comparing alternative models in a statistical context. Here we use representative species distributed in each biogeographic realm (Fig. 1), ranging from the Afrotropical, Palearctic (Europe and Mediterranean), Oriental, to Oceanic and Australian realms. To the best of our knowledge, this study is the first to infer the ancestral ranges of Rhinolophidae and Hipposideridae using explicit biogeographic analysis. In addition, here we attempt to infer the evolutionary history and biogeographic ranges of potential cryptic species of Asian Rhinolophidae (revised in [6]), in particular Southeast Asian region. We aim to understand historical biogeography of a subset of the superfamily Rhinolophoidea (Rhinolophidae, Hipposideridae and Rhinonycteridae) of Old World and the evolutionary history of potential cryptic species of Rhinolophidae.
Results
Historical biogeography of Rhinolophoidea superfamily
The time-tree pruned from Álvarez-Carretero et al. (2021) is provided in Additional file 1: Fig. S1, further, the time divergence estimations mentioned in this result were cited based on Álvarez-Carretero et al. 2021 study. In general, our result showed Likelihood Ratio Test (LRT) using chi-squared one-tailed test conferred higher log likelihood (natural log) for + j models (Table 1). However, the likelihood ratio test (LRT) shows there was no significant difference between DEC* + J and its nested models, DEC* (Lnl = − 302.6891; Lnl = − 302.6879, respectively, p = 1), shown the two models has the same likelihood on the data. While the DEC* + J model for Rhinolophoidea has some AIC support, its estimate of j is almost 0, therefore there is no support of founder events being important in this clade. To ensure the same likelihood in DEC* and DEC* + J was not due to the maximum likelihood optimization problems and computational issues, we ran the analysis three times and similar results were acquired. Here we focus on reporting the result under the best-fit model given the data, less-fitted models are provided in Additional file 1: Figs. S2–S9. The best fit model given the data is DEC* and followed with DEC* + J (Table 1). In DEC* model, extinction e value is higher than dispersal d (d < e; d = 0.0078; e = 0.0685) indicating the range contraction rate is higher than range expansion event in this superfamily (Fig. 2).
Our result suggests that the common ancestor of Rhinolophoidea superfamily as a whole [Rhinolophidae, (Hipposideridae, Rhinonycteridae)] were from the Oriental region (Fig. 2), diverged around middle Eocene at 43.25 Ma (95% highest posterior density/HPD = 39.74–46.76 Ma) with the percentage of relative probability was slightly higher from Oriental region than the model which showed an origin in the Afrotropical region (Fig. 2). Similarly, more than 50% relative probability showed that the ancestor of Rhinolophidae family also from Oriental region, and then later expanded to colonized Afrotropical region. The common ancestor of Hipposideridae and Rhinonycteridae were suggested to be from Oriental and Afrotropical regions (with almost the same relative probability for each region but slightly higher from Oriental), where Rhinonycteridae split earlier around middle Eocene at 40.55 Ma (95% HPD = 36.03–44.14 Ma) from the rest of Hipposideridae of Afrotropical lineages (Fig. 2).
Historical biogeography of family Rhinolophidae (Horseshoe bats).
The Rhinolophidae family split into two major lineages in the early Oligocene at 27.68 Ma (95% HPD = 23.96–31.54 Ma): Rhinolophus Afrotropical–Palearctic lineages and Oriental lineages (with Australian and Oceanic species). R. trifoliatus, R. luctus, R. formosae and R. hipposideros are sisters with respect to the above two lineages, and split earlier in the early Oligocene at 31.24 Ma (95% HPD = 27.36–35.16 Ma). The best fit model given the data was DEC* model followed with DEC* + J (Table 1). Similarly with Rhinolophoidea superfamily, range expansion is lower than extinction (d < e; d = 0.0156, e = 0.0444) suggesting that range constriction is higher compared to dispersal rate in historical colonization Rhinolophidae.
The result from DEC* suggests that ancestor of Rhinolophidae species the Oriental–Oceania–Australia lineage originated from Indomalayan region and species diversification began around Early Miocene around 22.77 Ma (95% HPD = 19.25–26.51 Ma) (Fig. 3). The ancestral range expanded from the Indomalayan region to the north to Sino-Japanese and west to India (within Asia continent) then later colonized towards the Philippine archipelago around end of Middle Miocene at 12.64 Ma (95% HPD = 9.37–16.19 Ma) (Fig. 3). The result shows ancestral Japanese rhinolophids originated from the Indomalayan region, the species thus diversified around the Middle Pliocene at 3.43 Ma (95% = 1.75–5.67 Ma) with high endemicity of extant species such as R. cornutus, R. pumilus, and R. perditus (Fig. 3a). The rhinolophid diversification in the Philippine archipelagoes is estimated to have occurred around the Late Miocene–Early Pleistocene at 6.09 Ma (95% HPD = 3.31–4.53 Ma). The Philippine rhinolophids appears to be the result from multiple colonization events originating from the Indomalayan and Wallacean regions. Rhinolophids colonization in Wallacea region, also occurred in multiple colonization event from the Indomalayan region through greater Sunda islands and Philippine archipelago, around Late Miocene at 7.12 Ma (95% HPD = 4.46–10.2 Ma). Australian and Oceanian rhinolophids appears to colonized from the Indomalayan region through two pathways, from the Indomalayan region–Sundaland–Wallacea and Indomalayan–Philippine–Wallacea which occurred in Late Miocene to Early Pliocene (Fig. 3a).
Moreover, the Afrotropical–Palearctic Rhinolophidae diversified around Early Miocene at 21.38 Ma (95% HPD = 17.74–25.36 Ma), the result exhibit BAYAREALIKE has 72% being the best fit model given the data, and ancestral range contraction rate is higher than dispersal or range expansion (d < e; d = 0.0137; e = 0.042755) (Table 1). BAYAREALIKE assumes the ancestral range process is predominantly explained by sympatry (subset and widespread) and no vicariance events occured. It is important to note that the relative probabilities in the nodes and the corner showed less than 50% probabilities for each states. Therefore, the ancestral species may be widely distributed in most of Afro-Palearctic region including Europe-Mediterranean, Sahara-Arabian, Sudanian-Somalia-Ethiopia, Zambezian-Kalahari-Southern Africa region (ABDE, Fig. 3b), or the combination of these regions. The ancestral ranges of Rhinolophidae predominantly occupied Sudanian + Somalia + Ethiopia region and Zambezian + Kalahari + Southern Africa around Middle Miocene to Pliocene (15–3.62 Ma). Re-colonization to Sahara-Arabian and Europe-Mediterranean might have occurred around Late Miocene at 7.63 Ma (95% HPD = 3.9–12. 14 Ma) (Fig. 3b).
Historical biogeography of family Hipposideridae (leaf-nosed bats) and Rhinonycteridae (the Trident bats)
The ancestors of the Hipposideridae family appears to have a wide distribution from the Oriental to the African region (see above results). The subset analysis results showed that DEC + J has a 47% probability of being the best fit model in the hipposiderids Oriental–Oceania–Australian lineages, alternatively, the other two models had strong support: DEC* (39% probability) and DEC* + J (13% probability). All trees generated in these models generally recovered the same ancestral range for each node, therefore indicated the similar scenarios. The single-most probable states in the output of the three models indicated Indomalayan + Philippines origins, however the relative probabilities given in early ancestral ranges were poorly supported (Fig. 4a). The result showed the dispersal rate is higher than range contraction (d > e; d = 0.00611, e = 0, j = 0.0288), which is in contrast with Rhinolophidae family, where dispersal rate is smaller than extinction rate. The hipposiderids ancestral range then expanded to India, Greater Sunda Islands, Wallacea and Oceania in the early Oligocene (30.52 Ma (95% HPD = 25.26–34.92 Ma). Colonization towards the Greater Sunda islands was directly from the Indomalayan region which occurred since the early Miocene. Oceanian hipposiderids appears to have colonized the region since early Oligocene and may use the Philippines as a stepping stone. The hipposiderids in Wallacea and Australia regions may have been colonized through greater Sunda Islands, as this coincided with historical presence of land bridges since Oligocene to Pliocene (Fig. 4a). Paraphyly within Hipposideridae is shown in African species H. jonesi (A = Sudanian + Somalia + Ethiopia) which was recovered within the rest of Oriental–Oceania–Australian taxa, diverged around early Oligocene at 28.12 Ma (95% HPD = 24.45–31.82 Ma), which may indicated earlier Asian–African colonization.
Furthermore, we found that DEC* has 79% probability of being the best fit model for Hipposideridae and Rhinonycteridae Afrotropical lineages, followed by DEC* + J (20% probability) with range contraction is higher than range expansion (d < e; d = 0.54, e = 2.43) (Table 1). The result suggested that the ancestral of Hipposideridae Afrotropical lineage were from the Sahara-Arabian (A) and Sudan-Somalian-Ethiopian (B) regions (Fig. 4b). The diversification of Hipposideridae Afrotropical lineages began around middle Eocene with splitting of Asellia tridens with the rest of Hipposideros genera (38.18 Ma (95% HPD = 36.93–41.79 Ma), originated from Sahara-Arabian region (however with little percentage of relative probability, Fig. 4b). Multiple colonization occurred to Guinean-Congolian region via range expansion from the Sahara-Arabian region through Sudanian + Somalia + Ethiopia and from Southern part (Zambezian + Kalahari + Southern Africa). Madagascar Hipposideridae ancestral originated from African continent (Fig. 4b). Similarly, the ancestral of Rhinonycteridae may have originated from Africa continent then colonized Madagascar which began at early Miocene at 21.88 Ma (95% HPD = 16.54–28.48 Ma). Our result suggests over-land dispersal event might occurred to expand Rhinonycteridae ancestral range from the Africa continents to Madagascar, however, it is necessary to consider that each of ancestral state in the result of this lineage were having almost similar relative probability percentage for each state and we discussed below.
The evolutionary history and historical biogeography range of potential cryptic species of Rhinolophidae Asia lineages
The result showed the estimation of time divergences between Rhinolophidae and Hipposideridae began in late Eocene, with posterior age 40.26 Ma (95% HPD = 37–43.5 Ma), slightly younger but in general fall within the dating range estimation in previous studies as 42 Ma in [14], 45.47 Ma in [4] and 43 Ma (95% HPD = 39–46 Ma) in [43]. Our analysis showed the diversification of Rhinolophidae in the region began around late Oligocene till early Miocene. R. JLEsp (undescribed species [69]; the acronym from [70] denoting Judith L. Eger), R. sedulus, R. trifoliatus and R. luctus represent the oldest lineages, separated from other Rhinolophus species around 24.52 Ma (95% HPD = 16.5–27 Ma), in concordance with Álvarez-Carretero et al. (2021) positioned R. luctus, R. trifoliatus and R. hipposideros (not included in this analysis) as the sister of Rhinolophidae species but their estimation is higher when Rhinolophidae Afrotropical lineages is included as of 31.24 Ma (95% HPD = 27.36–35.16 Ma)). Rapid speciation and species diversification in Indochina region happened approximately since late Miocene (23 Ma) (Fig. 5), coincided with date estimated from with Álvarez-Carretero et al. (2021) [22.77 Ma (95% HPD = 19.25–26.51 Ma)]. The majority of potential cryptic species or incipient species within rhinolophids species complexes diverged in mid-Pliocene within the last 2 Ma indicates rapid radiation in Plio-Pleistocene, with the exception of the R. pearsonii complex that diverged during mid-Miocene 7.38 Ma (HPD = 5.2–9 Ma) indicated a support to split R. pearsonii into multiple different species. The time divergence construction using secondary calibration at multiple nodes in BEAST supported the evolutionary distinction between the clades within species complexes and the incipient species might further be considered as distinct species.
The result of historical geographic ranges estimation for cryptic Rhinolophidae suggests DEC* + J has 82% probability of being the best fit model given the data and DEC* (Dispersal–Extinction–Cladogenetic) has 18% probability of best fit model, and other models shown < 1% of probabilities as fit model given the dataset. DEC* + J showed the estimate extinction rate (range contraction) events was higher compared to dispersal (range expansion) and jump-dispersal events (d < e; d = 0.0205, e = 0.298, j = 0.0228).
Cryptic Asian Rhinolophidae likely originated from the Indomalayan region, here in particular Indochina (Fig. 5), in concordance with our result in ancestral biogeographic of Rhinolophidae Oriental–Oceanian–Australian lineages, then expanded to current geographic ranges subsequently. The differences in the results we present here with other section (Rhinolophidae Oriental–Oceanian–Australian lineage) was the biogeographic areas and systematic coverages we used in the analysis (see “Materials and methods”).
The origin of R. luctus, R. sedulus and R. trifoliatus ancestors were in the Indomalayan region then occupied Borneo island via vicariance events since middle Miocene (13 Ma). Similarly, colonization to Borneo Islands from Indochina also present in R. philippinensis and R. borneensis around Pliocene (2–5 Ma). Our result suggests a rapid diversification began in the late Miocene around 7 Ma within R. pearsonii as the ancestral group where the ancestral ranges subset-sympatry into two largely ranges; some remained in Indochina and then subset-sympatry in northern areas (Himalaya region and China).
Our analysis indicates a high probability of recent diversification through subset-sympatry in R. affinis and R. pusillus group around Isthmus of Kra, Kangar-Pattani line to the tip of Malay Peninsula from the rest of northern population. Genetic drift followed by isolation and divergence may explain the emergence of separate geographical ranges of R. affinis in Sumatra Island after separation from mainland which shows a closer relationship to Isthmus of Kra, Kangar Pattani and Malay Peninsula lineages dated around the Pleistocene (2 Ma) (Fig. 5).
Discussion
The ancestral range estimation for superfamily Rhinolophoidaea [Rhinolophidae, (Hipposideridae, Rhinonycteridae)] suggests that the ancestors of family were had an Oriental origin, then diversified to other biogeographic regions. African region were colonized through Sahara-Arabian regions, ancestral range expanded to Madagascar. The Oriental region experienced complex colonizations related to geological events in the region, including colonization toward Philippines, Indonesia (Greater Sunda Islands, Wallacea) and Oceania. Australian rhinolophoids may have colonized through the Philippine, Greater Sunda Islands, Wallacea and Oceania. Our results are congruent with the previous hypothesis that Rhinolophidae bats may have originated from the old world tropics in Asia [18, 23, 30, 73] then the range expanded to East Asia region including India, Africa and islands in the vicinity, but, in contrast with [14, 22] which hypothesized that the ancestral of Rhinolophus species were of African origin based on LAGRANGE biogeographic analysis. Rhinolophidae of Oriental–Oceanian–Australian lineages originated from the Indomalayan region, and the ancestral range of Afrotropical lineages were widespread in the Afro-Palearctic. Allozyme variability also suggests colonization from Eurasia toward North Africa and that subsequent diversification took place in Africa [74], and morphological studies suggest that plesiomorph Oriental rhinolophids were basal and Afro-Palearctic species were more derived [30]. The widespread ancestral and current distribution of taxa in best-fit model BAYAREALIKE indicated the ancestral ranges were similar to those occupied by their descendants, and suggests over-land range expansion in the past around Late Oligocene coincided with diversification of Afrotropical lineages at 27.68–24.38 Ma. The result showed that dispersal rate was constantly lower than extinction rate in Rhinolophidae (d < e). Similarly with Hipposideridae–Rhinonycteridae Afro-Palearctic lineages (d < e), but in contrast with Hipposideridae Oriental lineages where dispersal or range expansion is higher (d > e). The estimated rate of dispersal is low suggesting that species have mostly retained the same geographic ranges as their ancestors [36, 38].
The common ancestor of Hipposideridae originated from Oriental or Afrotropical regions (similar posterior probability, Fig. 2), suggesting possible widespread ancestral ranges of this family. The ancestor of Hipposideridae Oriental–Oceania–Australian lineages originated from the Indomalayan region and the Philippines, which suggests multiple colonizations may occured during range expansion. The common ancestors of Hipposideridae Afro-Palearctic lineages originated from Sahara-Arabian and Sudanian–Somalia–Ethiopia suggested the early colonization from northern part of Africa and Arabian Peninsula (Fig. 4b). Even though the inferences of ancestral ranges of Hipposideridae lineages are not strongly supported, but the results are in agreement with previous studies based on fossils which suggest that major the dispersal axis of Hipposideridae was from North Africa toward South Europe during the Middle Eocene [75]. The result also suggests that Rhinonycteridae in Madagascar split from their common ancestor (Hipposideridae Afrotropical lineages) around Middle Eocene (~ 40.55 Ma) and ancestral range originated from the Sudanian–Somalia–Ethiopia regions.
Our result supports the importance of statistical model testing as each model has strong assumptions which impact on inferences, as highlighted in [36]. Although, recent results by [71] revealed potential statistical problem for the comparison between DEC and DEC + J models and to alleviate this issue they proposed the use of geographic state speciation and extinction models (ClaSSE) to study range evolution. However, [72] argues that DEC and DEC + J can be compared using likelihood and AIC, and that the comparison is equivalent to ClaSSE comparisons. Despite these controversies, our results showed that model with + J were not always selected as the best fit models and the simulation results in [72] indicated better fit models will result in more accurate inference of parameters and ancestral ranges. However, it also important to consider that biogeographic range estimation can be biased by model choice.
Even thought our result indicates an Oriental (Indomalayan) origin for the Rhinolophidae, there is lack of Paleogene fossil found in the Indomalayan region, though variable taphonomy means that the majority of the fossil record is missing across taxa. The only known possible record of Eocene bats in this region is Megachiroptera from Krabi Mine in Thailand. However, fossil evidence of bats is relatively rare due to delicate skeletons and are therefore rarely preserved, thus, leaving only teeth and postcranial fragments for identifcation [76]. Fossils of ancestral Rhinolophidae include Protorhinolophus shanghuangensis from Shanghuang fissure, Jiangsu (northern part of Asia) [73], highlighting ancient distrobution in the Oriental region. The oldest bat fossils are from the early Eocene, and are known from North America, Europe, Africa and Australia [76, 77]. However, geographic origin of bats ancestor still being debated with two hypothesis of Laurasia or Gondwana origin [78]. The initial explosive radiation of bats occurred in Eocene, the extinct families such as “Eochiroptera” sensu Van Valen (1977) found in most of continents except Antartica [25, 77]. Modern radiation of extant taxa appears to have begun at least by Middle Eocene or earlier, a period characterized by a significant global rise temperature after K-Pg (Cretaceous–Paleogene) mass extinction event [19, 76, 79], which coincided with 43.25 Ma divergence between Hipposideridae and Rhinolophidae.
The oldest Rhinolophus fossils was Protorhinolophus shanghuangensis fossil found in Shanghuang fissures, China aged Middle Eocene [73], making it older than the Vaylatsia prisca (Rhinolophus priscus) fossil found in Europe (early Late Eocene to Early Oligocene). The oldest known fossil of Rhinolophus priscus is dated to the Late Eocene to Oligocene in the Quercy of France (Europe), but Protorhinolophus showed more primitive dental patterns, which indicates the genus is older than Rhinolophus [73]. The oldest fossil known in Africa, R. mellali of Bani Mellal in Morocco, North Africa were dated in Late Miocene to Early Pliocene, which is probably close related with R. ferrumequinum [30]. The estimated split between Rhinolophidae and its sister lineages, Hipposideridae in the middle Eocene were supported with Hipposideros fossils in North Africa Pseudorhinolophus africanum indicated the major dispersal axis of the family from North Africa to South Europe [75] (Figs. 6, 7). In contrast with Rhinolophidae, the Hipposideridae have diversified in Europe, Africa and Arabia since the middle Eocene [73], this agrees with our ancestral range estimation result, and the widespread ancestral range of Hipposideridae (Figs. 2, 6).
Another published study hypothesized that European origin of Rhinolophidae [80] but had weaker support for species in basal clades. Bogdanowicz and Owen (1992) also hypothesized that rhinolophids dispersed from Asia to Australia through Indonesia and New Guinea based on morphological data, which indicates the ancestral of Rhinolophus were from Asia, which is consistent with our result, and the higher rhinolophid diversity in the Asian region [6]. The colonization of Rhinolophus is also limited by over-water barriers, therefore limits the dispersal from Eurasia by the presence of Tethys oceans in Early Miocene (Fig. 7). In the late Oligocene and Early Miocene, the Tethys was connected the Proto-Mediterranean Sea to the Indian Ocean [81] (Fig. 7), thus may have inhibited the ancestral range expansion from Eurasia continent to Africa continent. The Afro-Arabian and Eurasian plates collided during the Middle Miocene, caused the presence of land-bridges Gomphoterium landbridge due to Tethys Ocean shrinkage around Late Miocene [81,82,83,84] (Fig. 7), this event may allow species to disperse from Eurasia to the Afro-Arabian continent. In general, our result suggested the ancestral Rhinolophidae originated from the Oriental region then subsequently the range expanded towards Africa through Europe, and the ancestral Hipposideridae were widespread from Afro-Arabian region to Eurasia, consistent with Pseudorhinolophus africanum fossils in Middle Eocene in North Africa [75].
Ancestral biogeography range of horseshoe bats (Rhinolophidae)
Four subset analysis (Rhinolophoidea, Rhinolophidae Oriental–Oceanian–Australia lineages and Rhinolophidae-Afro-Palearctic lineage, and cryptic Rhinolophidae) suggests estimation of dispersal rate was lower than range constriction (d < e). The high e in best fit model demonstrated all range-changes effectively occurred anagenetically along the branches [38]. Thus, we assume the distribution of Rhinolophidae is driven primarily by dispersal (i.e., over-land range expansion, land-bridge colonization, and stepping-stone events) mixture with vicariance events. Additionally, jump-dispersal or founder-events was not well supported for Rhinolophidae ancestral ranges especially when we used bigger biogeographic ranges. However, jump-dispersal events were chosen as the best fit model for explaining oceanic-dispersal pattern when we separate each islands as different biogeographic units (i.e., Greater Sunda Island: Sumatra, Borneo and Java). Range expansion by dispersal, extinction, sympatry (subset and narrow) and vicariance events are useful to understand historical biogeographic of Rhinolophidae in Oriental–Oceania–Australia lineages, however, ancestral range of Rhinolophidae Afro-Palearctic is a better fit with dispersal, extinction, sympatry (narrow and widespread) events (without vicariance). The understanding of how organisms came to be distributed as they are also related to historical events involving complex geological history, such as glacial-inter glaciation, continental drift, biotic turnover and long-distance colonization [28]. The main event in Cenozoic era that shapes the continent including collision between the Indian plate and the Eurasian plate, created mountain ranges in Himalaya regions and acted as physical barrier for species dispersal between Indochina and India in late Eocene (around 50 Ma till present). The continuation of India attachment to Eurasia hugely affects the movement of other plates and influenced the shape of continents, archipelago and affects climatic condition in the region (Fig. 7) [46, 50, 57, 58, 85, 86].
The rapid diversification during Miocene coincided with the Mid-Miocene Climatic Optimum (around 15–18 Ma) which may provide a favourable climate in support for the evolution [87] of the Rhinolophus lineages. The high diversification in the Miocene not only occurred in bats, but also in the diversification of modern bird genera in Southeast Asia [61]. The basal position of two lineages in Rhinolophidae are the three species belonging to the trifoliatus group (R. trifoliatus, R. luctus, R. formosae) and R. hipposideros. The three trifoliatus species are Oriental species distributed from Indian sub-continent, Southeast Asia and Eastern Asia [5], and R. hipposideros is distributed throughout the Europe from Ireland in the northwest to Pakistan in the east, and south into northern regions of Africa and Saudi Arabia [88]. Phylogeographic studies suggests early colonization event of R. hipposideros and R. ferrumequinum from the east (west Asian refugium) and both of species used multiple glacial refugia across Mediterranean during the ice age [88, 89]. Last Glacial Maximum in late Pleistocene have impact on current distribution of R. ferrumequinum, with secondary contact was identified between Central/East China and East China/Japan [89].
Our results also show colonization of Rhinolophidae species toward Japanese archipelago may have occurred since the Middle Pliocene, which coincided with fossils from caves and fissure deposits of Middle Pleistocene, Late Pleistocene and Early Holocene in Honshu and Kyushu Island (see [90]). The colonization from continental population toward Japanese archipelago may occurred via a continuous Korean Peninsula-Japanese land bridge due to lower sea levels than at present [89], similar to the colonisation route for Asian black bear (Ursus thibetanus) [91].
Our result showed the colonization towards the Philippines appears to be multiple colonization events via different routes originating from the Indomalayan region (which may be via Palawan) and Wallacea region in late Miocene (Fig. 7). Similarly with Wallacean and Oceanian region, the ancestral range and colonization events in these region are complex. The thousands of islands in these regions have different continental, oceanic and volcanic origin, with many of them undergone rapid tectonic movement since Cenozoic [92, 93]. The archipelagoes (including Philippine and Indonesia) with the combination of tectonic movement, climatological oscillations and Pleistocene sea level fluctuations causing the changes of island size, connectivity and boundaries [94,95,96]. These important geological events may have contributed to current species distributions [97, 98], and as the effect of presence or absence of the Pleistocene land-bridge connection [95]. In addition, current distribution of species in the archipelago may be a direct result of when bats reaching the older islands longer ago and younger island more recently which leading to allopatric and vicariance speciation [99].
Two hypothesis of species biogeography in the Philippines including Pleistocene aggregate island complexes (PAICs) suggesting land exposure in the Pleistocene due to glaciation allowed species to expand their range inter five major islands [100]. Palawan Ark Hypothesis suggested species “rafted” with the North Palawan block since the separation from mainland Asia by Early Oligocene 30 Ma [50], in contact with Borneo around 15 Ma [96] and then move northward toward present position [101, 102]. The multi-route colonization of Rhinolophidae species toward the Philippines is similar to other taxa, for instance Begonia [103], Cynopterus and Macroglossus [104, 105]. This coincides with the theory of biotic colonization of Philippines, postulated as submerged land bridges, with many taxa known to have colonized Philippine through northern Philippines (from mainland Asia and Palawan) and through south route (via Sulu archipelago and Sulawesi) [106]. Previous studies also suggested that colonization of the Philippines may have taken place from the Sunda Shelf (Sumatra, Java and Borneo) and Wallacea [95].
Wallacea, including Sulawesi Island, and many small islands surrounding (i.e., Outer Banda Arc (Sumba, Timor, Babar, Yamdena, Kai, Seram), Inner Banda Arc (Bali, Lombok, Sumbawa, Flores, Alor and Wetar), Halmahera etc.) result from complex tectonic plate movement from Australian and Asian plate. Some of the Island were never connected in the past (e.g. Inner and Outer Banda Arc), and Inner Banda Arc exposed the land above sea-level and therefore permitted the colonization in Lesser Sunda Islands around 3 Ma [107]. The current arrangement of the islands provides a series of stepping stones facilitating movement of terrestrial mammals, which may include volant mammals, to colonize the Australian region [47].
Furthermore, glacial and sea-level fluctuation repeatedly formed land bridges in Pleistocene [108] and the landmass between Asia, Sumatra, Borneo and Java islands (Sundaland) [50] allowed colonialization and range expansion of Rhinolophidae from Indochina towards Sumatra and Borneo, which is in line with our results here. The diversification of Rhinolophus species in Borneo Island may coincide with separation of Sundaic and Indochinese rainforests [53]. The rainforest refugia in some parts of Borneo and Sumatra may allowed populations of forest-species to diverge and adapted with local climatic conditions and environment [53, 109, 110]. The glaciation and interglacial events during Plio-Pleistocene caused dramatic changes in climate, forest cover and the connection between land areas allowed species to colonized different geographic ranges. Recent diversification within species groups in Pleistocene may coincide with climatic fluctuations which affects the vegetation transition in the region, and the possible savanna corridor in part of the region [109, 111], thus indirectly influencing diversification of forest-dwelling mammals [112] and insects [113]. The divergences of R. affinis and R. pusillus lineages in Malay Peninsula during Plio-Pleistocene may coincide with major events in the peninsula. The possible flooding of Isthmus of Kra during Pliocene, and persistent climatic transition in this region [110], the adaptation to specific climatic conditions and long term ecological differences, combined with the peninsula effect may cause the faunal transition between Indochinese taxa and Sundaic taxa, with major transitions at the Isthmus of Kra, and the Kangar Pattani line [56, 110, 112, 114] (Fig. 7). Furthermore, the sea-level rises during Pliocene isolated Sumatra, Borneo and Java Islands in Indonesia created physical barriers between Indochina and species from the islands [50, 61]. Borneo is the largest landmass of former Sundaland and has served as stable land for at least 20 million years and was less affected by sea-levels changes compared with Sumatra (which has come together as a stable from 5 to 10 Ma) and Java (2–5 Ma) [53, 61].
High rhinolophid diversity in Asia compared to the other biogeographic regions may be expected as many species are restricted to islands or group of islands, for instance four endemic species in Japan (i.e., R. cornutus, R. pumilus, R. perditus, R.imaizumii) (Ohdachi et al. 2015) and at least four endemic to the Philippines (i.e., R. inops, R. rufus, R. subrufus, R. virgo) and various species endemic in Indonesia islands (i.e., R. nereis, R. madurensis, R, keyensis, R. montanus, R. euryotis, R. celebensis, R. canuti [115], which represents physical barriers for bat dispersal such as water and mountain ranges. In contrast, Palearctic regions were influenced with repeated glaciations and in Africa, that may cause high species turn-over and high rates of gene flow which decelerating speciation compared to complex biogeographic areas in Asia [74]. The dispersal of Rhinolophidae species to Africa may be via forest corridors that appeared during middle Eocene because of the warm climate [116], and the dispersal of ancient species throughout the southern Palearctic and Mediterranean, through Sahara Arabian and eventually into Africa [18]. The basal lineages within African radiation are R. landeri and R. alcyone, which occur in rainforests and may indicate early colonization around the late Oligocene (~ 20 Ma) through forest affinities as predicted in [20]. The closure of the Tethys Ocean us estimated around the Oligocene (around 27) Ma to Miocene [81, 117, 118]. A land bridge in Arabian Peninsula which connected Asia and Africa, and may have facilitated dispersal of many animals lineages in the Miocene such as lizards [119], frogs [120], chameleons [121], and butterflies [122, 123]. As the consequences of shrinkage of Tethys sea, desert and arid conditions expanded across North Africa in the Late Miocene (around 7 Ma), marking the origin of Sahara Desert and the Middle East Desert and the Arabian Peninsula [82]. Arid adapted species colonized Africa and now inhabit most of savanna region in Sudanian, Somalia, Ethiopia and deciduous woodland in Southern Africa [21, 30].
Unlike Hipposideridae, there are currently no records of Rhinolophidae species in Madagascar. Madagascar + India + Africa are ancient fragments of Gondwana and has been separated from Gondwana since 120 Ma and, and Madagascar separated from India by 90 Ma [124]. Madagascar started to break away from Africa around 165 Ma [125] and to become isolated in Cretaceous [126], and the invasion from Africa continent toward the island may not occurred due to broad watergaps as physical barriers for dispersal as Rhinolophidae are weak fliers [30]. Additionally, higher diversity richness and endemism in Madagascar appears as a result of dispersal from Africa and followed with diversification [127, 128], and typically reflects more recent events around Plio-Pleistocene [129].
The d and e parameters may show different rates depending on the bats family, for instance in Pteropodidae [130]. We assume the historical biogeography in bats varied between families which may relates to flight performance between species and differences in dispersal ability. For instance bat colonization to Madagascar includes most bat families (Pteropodidae, Emballonuridae, Hipposideridae, Vespertillionidae, Nycteridae, Molossidae and Myzopodidae), with a notable of absence Rhinolophidae species (Racey et al. 2009). Wing-loading, wing-span, aspect ratio and wing shape are the main aerodynamic variables in determining the flight performance and flight efficiency of species. Rhinolophids bats generally possess broad and short, low wing loading and aspect ratio adapted for good maneuverability in foraging as slow aerial hawkers, perch hunters and gleaners (Norberg and Rayner 1987; Amador et al. 2020). However, flight would be expensive for a long distance flight and therefore limiting the ability for dispersal. Contrary to Rhinolophidae, most Pteropodidae possess high wing loading and large wing spans, with low aspect ratios that support for long distance travel but are less maneuverable [3], as a consequence many Pteropodids are restricted to various oceanic islands where they are often the only native mammal species.
True-flight abilities together with echolocation have been long considered as remarkable evolutionary features which have driven the success of bat species and enabled them to occupy wide range of ecological niches, but the morphological features have evolved into current form from a common ancestor [2]. Other factors such as dispersal filters (isolation, geologic history), environmental filters (present and past climate, and environmental heterogeneity) have also influenced the diversity of vascular plants [131] and may also have contributed in shaping the current biogeographic ranges of Rhinolophidae.
The evolutionary history of species divergence in potential cryptic Rhinolophidae and the comparisons with previous study
Indomalaya is a region with particularly high species diversity of Rhinolophidae and Hipposideridae. However, many species are cryptic, suggesting the number of species currently underestimated [70]. Around 40% of rhinolophids species are potentially cryptic based on integrative taxonomic approaches, with estimated around 44 potential cryptic species from total of 10 Rhinolophus species sensu lato [6]. However, a comparison of time divergence estimation across the genes is challenging because of the limited number of species used in previous studies, and the paucity of the bat fossil record. We acknowledge the disconcordance of time tree in this study with other studies, which was also due to differences in systematic taxonomic coverage and genes being used (as described in the “Introduction” section). Although we found some discordance between tree topologies we generated in this study and from previous studies, we attained similar results across subset analysis related to ancestral range estimation in Rhinolophidae. In terms of species coverage from South China and the Southeast Asian region, we covered most described species from each geographical region. Our results suggest current species sensu lato started to diverge in the late Oligocene and Miocene, meanwhile the potential cryptic species in the region diverged in the Plio-Pleistocene epoch. This analysis indicates that geological events during the epoch contributed in shaping current cryptic diversity patterns seen today. Future directions should aim to include all Rhinolophoid species distributed in the Old-World tropics, as we currently cover around 108 (52%) species of total species identified to date (210 species; Rhinolophidae = 110, Hipposideridae = 90, Rhinonycteridae = 9 species) (Simmons and Cirranello, 2021).
The systematic and ancestral biogeography range of Hipposideridae and Rhinonycteridae
Hipposideridae and Rhinonycteridae are sister families within Rhinolophoidea superfamily. Amador et al. (2018), suggested the phylogenetic tree for superfamily as (Rhinonycteridae, (Hipposideridae, Rhinolophidae)) but Foley et al. (2015) resolved the relationship as (Rhinolophidae, (Hipposideridae, Rhinonycteridae)), similar with species tree in this study [43]. Foley et al. (2015) suggested the African species of H. abae, H. caffer and H. jonesi were within Asian Hipposideros, and Coelops frithii, Aselliscus stoliczkanus were sister to Hipposideros African species (H. commersoni and H. vittatus). Nevertheless, this arrangement has low branch support, and low taxonomic coverage which only includes total of 13 species of Hipposideridae Asian lineages and Afrotropical lineages. In the species tree use in this study, [43] resolved the relationship and H. jonesi (Afrotropical species) species fall within Oriental–Australian lineages, similar with tree topology of Amador et al. (2018), paraphyly with Aselliscus sp., Coelops frithii and Anthops ornatus in the same lineage, which is in agreement with [132]. H. abae and H. caffer as Afrotropical species are within Afrotropical lineages (in contrast with Foley et al. 2015, but in concordance with [60]) with Asellia tridens in the basal of Afrotropical and Oriental–Australian tree.
For Rhinonycteridae, here we include seven species out of nine species within the family, belonging to four genera, Triaenops, Paratriaenops, Cloeotis, distributed in Madagascar and Afrotropical, and Rhinonycteris endemic to Australia with tree topology as [(Cloeotis, Paratriaenops), (Rhinonycteris, Triaenops)]. This arrangement also in contrasts with Foley et al. (2015) with tree topology as [(Paratriaenops, ((Triaenops, Cloeotis), Rhinonycteris)] while using six species, and Amador et al. (2018) [Paratriaenops, (Triaenops, Cloeotis)] with absence of genus Rhinonycteris. Certainly, missing taxa and gene choice compromise these results to at least some degree, thus the tree generated in Álvarez-Carretero et al. 2021 using supermatrix genome-scale of 182 genes may assume as the latest update tree with smaller uncertainties which facilitates precise testing of historical biogeographic analysis.
The historical biogeography of family Hipposideridae is also related to complex historical geology of the Old-World regions (see discussion above) (Fig. 7). The ancestral ranges of Hipposideridae suggests an Oriental and African origin, the Oriental–Oceania–Australia lineages supports the jump-dispersal for species colonization. This explains how new lineages colonized between regions (such as island), and allowed inter-continental disjunction or oceanic-dispersal pattern (dispersal without range expansion) (Matzke 2014).
In contrast with Rhinolophidae, their dispersal rate is higher than range constriction in Hipposideridae. This indicates that ancestral Hipposideridae were able to disperse across water-bodies or used land-bridges and stepping stones. The colonization of Hipposideridae towards Greater Sunda Island may have occurred when the land was connected with Peninsula since early Miocene, followed with multiple fluctuation of sea level until the Sunda Island formed separately at 5 Ma [50, 57, 58, 86], similarly with Rhinolophidae. The Wallacea region might have played an important role as a stepping stone to colonize Oceania and Australia regions, but migration, population exchanges and secondary contact may have occurred in the past between Oceanian-Wallacean region. In concordance with Rhinolophidae ancestral ranges, the range expansion of Hipposideridae to Africa may also occurs via Arabian Peninsula and Sahara-Arabian around Oligocene periods (27–30 Ma). Even though the Gomphoterium landbridge and Tethys Ocean closed around the Late Miocene, the marine barriers did not totally prevent mammalia exchanges between Eurasia and Afro-Arabian, for example for proboscideans [133]. Ancestral Hipposideridae were widespread in Eurasia and Africa since the end of Early Eocene with the evidence of fossils records in Africa and Arabia [75].
The colonization of Hipposideridae towards the Philippines, Wallacea, Oceania and Australia coincided with acceleration of the orogeny of Philippine archipelago, Wallacea with the orogeny of Sulawesi and the main stages of the New Guinean orogeny in Oligocene [50, 86, 122]. Our result suggests Oceanian colonization since the early Miocene (~ 20 Ma), similarly with passerine birds by [134]. However, central range of present day of Papua New-Guinea likely did not begin to appear as land until the early-middle Miocene (14–16 Ma) [135] and the present form is predicted since 4–5 Ma. Therefore, founder-events may have involved in island-hopping across the final fragments of a proto-Papuan archipelago in Hipposideridae [136].
Rhinonycteridae split earlier from the common ancestor Hipposideridae Afrotropical lineage around Middle Miocene 40.55 Ma (95% HPD = 36.98–44.14 Ma), but species diversification occurred since 21.88 Ma (95% HPD = 16.54–28.48 Ma) with most extant species distributed in Madagascar, except Rhinonycteris aurantia which is distributed in Australia. Our result suggests the origin of Rhinonycteridae in Madagascar was Sudanian + Somalia + Ethiopia with Paratriaenops and Triaenops was recently emerged in Plio-Pleistocene epoch (1–3 Ma). This coincided with previous studies, which stated that most the present-day organism in Madagascar are of predominantly African origin [127]. Dispersal over-water may explained the colonization of ancient Rhinonycteridae in Madagascar, followed with diversification and in-situ radiation, which coincided with Plio-Pleistocene climate cycles [129]. The ancestral Cloeotis percivali may have high dispersal ability allowed the ancestor to exchange from Madagascar to Southern Africa. This associated with our result that suggests higher extinction rate in these lineages, indicating that the extinction of species in the ancestral range was followed by the colonization of descendants to the new area. The endemic species in Australia belongs to Rhinonycteridae, Rhinonycteris aurantia is the only species that currently distributed outside Africa continent [15]. The colonization event of this species towards Australia is challenging to explain, our result shows an African origin but it is almost impossible to explain range extension from Africa to Australia with diversification around Miocene. Some possible hypothesis maybe long distance dispersal over water barrier from Africa to Australia, however the dispersal mechanisms are unknown considering the species are weak fliers. Some of hypotheses suggest that bats colonized Australia by storm-blown to continental shore, for example red flying foxes (Pteropus scapulatus) [137], though landbridges between North Australia and Papua existed for extended periods of time. Ancestral Rhinonycteris may have entered Australia before the Miocene, with fossil evidence of Rhinonycteris tedfordi dating around Miocene from Riversleigh, Northwestern Queensland [138]. Alternative hypothesis suggests waif dispersal and stepping stones through Asia, Sundaland, Wallacea and Oceania [137], therefore, how this species colonized the region still debatable, but low diversity means that species in intermediate regions may have become extinct.
Finally, as we stated in the results above, the inferences of early Hipposideridae ancestral ranges are weakly supported for each state. This issue may can be alleviated if the Hipposideridae are analyzed as a whole, however some extant species were widespread (i.e., H. ater, H. cervinus, H. galeritus, H. diadema) in 6–7 zoogeographic regions (India, Indomalayan, Greater Sunda Islands, Philippines, Wallacea, Oceania and some in Australia, though further analysis may reveal undescribed species in these groups), causing large computational demands to analyze whole family. Reducing number of zoogeographic regions would facilitate processing, but would give the similar results as ancestral range reconstructions presented in Fig. 2. In addition, further systematic research needed for these widespread species that may contain undetected cryptic diversity [139]. For instance, there are multiple potential species listed as H. ater within Southeast Asia and Australia [132]. Resolving the systematic and taxonomic issues, which if resolved would strengthen the inferences of ancestral range estimation within this family.
Conclusion
The results presented here provide an overview of biogeographic history of extant Rhinolophoidea superfamily (Rhinolophidae, Hipposideridae and Rhinonycteridae). The result suggests that the Oriental origin, and the Indomalayan is the ancestral origin of Rhinolophidae and Hipposideridae, but Rhinonycteridae ancestors originated in Africa. The result showed a lower dispersal rate (range expansion) than extinction rate (range constriction) in Rhinolophidae, but in contrasts with Hipposideridae Oriental–Oceania–Australia lineage. Jump dispersal events explain the ancestral range in Hipposideridae Oriental lineages. Current distribution of extant taxa appears to be a result of a combination of dispersal, extinction, sympatry and vicariance events, followed with complex geological history of the Old-World regions. Further study combining complete species coverage of Rhinolophidae, Hipposideridae and Rhinonycteridae may increase the resolution of the historical biogeography of these taxa.
Methods
Sampling, systematic coverage and biogeographic sampling
In total we included 104 species of Rhinolophidae (n = 47; 43% from 110 total described species), Hipposideridae (n = 50; 56% of 90 species) and Rhinonycteridae (n = 7; 78% of 9 species) in this study. Taxonomic sampling coverage was based on number of species used in the tree published in Álvarez-Carretero et al. 2021. The total species coverage per biogeographic area was; Australia–Oriental–Oceania (thereafter, AOO) = 31; Afrotropical–Madagascar–Palearctic (thereafter AMP) = 16 for Rhinolophidae and AOO = 38; AMP = 12 for Hipposideridae. The term of AMP is used for all subset analysis for consistency, but it is important to note there is no record of Rhinolophidae species in Madagascar to date (see details in “Discussion”). In addition, total of 26 Rhinolophidae sensu lato mainly distributed in Southeast Asia and their potential cryptic species (44 potential cryptic [6]) was used to assess their evolutionary history and ancestral biogeographic range.
Biogeographic mapping
To assess species routes based on fossil data and assess the correspondence between fossil data (where present) and ancient coastlines we downloaded recent files showing coastlines for periods between the present (0 Ma) and 45 Ma, including 37 Ma, 30 Ma, 22 Ma, 13 Ma based on [40]. Fossil data from Rhinolophids was then downloaded from the Fossilworks (http://www.fossilworks.org/cgi-bin/bridge.pl?a=taxonInfo&taxon_no=40644) and Global Biodiversity Information Facility (GBIF), using a criteria of fossil species and the search term Rhinolophidae (https://doi.org/10.15468/dl.uy9gx5). This was repeated for Hipposideridae (https://doi.org/10.15468/dl.2826pj). Data was cleaned to give putative species name at the highest available level, and age of fossil in millions of years, dates were also checked from associated publications when not available in the spreadsheets. This data was rounded up to the nearest million years (due to uncertainty and ease of analysis and display), and where not available for any given fossil the dates for the same species in other parts of the region was used where present for mapped ranges of species.
Biogeographic analysis of old-World Rhinolophidae and Hipposideridae
We performed ancestral geographic range analyses using probabilistic modelling in R package “BioGeoBEARS version 1.1.2” [39]. We statistically compared the likelihood-based model of geographic range evolution of DEC model (Dispersal–Extinction–Cladogenesis) of LAGRANGE [41], a likelihood implementation of the processes assumed by parsimony of DIVA [42] (therefore named DIVALIKE), a likelihood version of the range evolution model of BAYAREA [37] (therefore named BAYAREALIKE) and a modification of DEC model by prohibiting the transition into null-range (DEC*) [38]. Each model is fully parameterized in BioGeoBEARS supermodels with different assumptions about anagenetic and cladogenetic change processes. A free parameter of “j” (jump dispersal or founder-event speciation) was added in DEC + J, DEC* + J, DIVALIKE + J, BAYAREALIKE + J models and nested with another two free parameters (d and e) within DEC + J, DEC* is nested within DEC* + J, DIVALIKE is nested within DIVALIKE + J and BAYAREALIKE is nested within BAYAREALIKE + J [39].
To infer the ancestral biogeographic ranges of Rhinolophidae and Hipposideridae, we used a dated-tree generated from [43] which used a supermatrix of 33.2 × 106 alignment length from 72 genomes (15,268 genes). The date estimation in this published study produce smaller uncertainties in dates estimates, thus are recommended for macro-evolutionary studies [43], including ancestral biogeography (Additional file 1: Fig. S1). There are some limitations in statistical analyses of ancestral range evolution. The 14 areas coded here would produce a huge transition matrix (214 = 16,384 possible area combination; 16,384 × 16,384 transition matrix (268,435,456 total) that has to exponentiated across each tree branch for maximum likelihood optimization. This large matrix is currently not feasible with matrix-handling algorithms. Therefore, here we analyzed the ancestral ranges in 4 subset analyses, including: (1) global Rhinolophoidea superfamily (Rhinolophidae, Hipposideridae and Rhynonycteridae), we then subdivided the analysis into sub-regions for each family into: (2) Rhinolophidae Asia–Australia lineages and Rhinolophidae Afrotropical–Palearctic lineages, (3) Hipposideridae Asia–Australia lineages and Hipposideridae–Rhinonycteridae Afrotropical lineages, (4) Potential cryptic species of Rhinolophidae in Asia lineages (potential cryptic species were delineated based on integrative taxonomic approaches as reviewed in [6]). This analysis strategy was to optimize the computational efficiency by considering the possible numbers of states for each sub-analysis. The input files for analyses required a dated-tree and geography file. We pruned the species tree published from Álvarez-Carretero et al. 2021 to Yinpterochiroptera lineages and trimmed the branches to Hipposideridae, Rhinonycteridae, Rhinolophidae using the function drop.tip in R package ‘phytools’ v1.0-1 [44], and removed the outgroups. We merged the subspecies in the tree tip into species level (OTUs) with taxonomy following Simmons and Cirranello 2021, using R package BioGeoBEARS v1.1.2 [39] in the function prune_specimens_to_species. The geography files are PHYLIP-formatted files, and were also used for C++ LAGRANGE editor file, which was generated in R using function save_tipranges_to_lagrangePHYLIP. To determine the biogeographic region, we used the current distribution of extant taxa from the data provided in previous literature and combined this with data from the Global Biodiversity Information Facility GBIF (https://www.gbif.org/) (https://doi.org/10.15468/39omei), as IUCN ranges of bats may be unrepresentative and risk being particularly inaccurate [45]. We mapped the current occurences record of each species and codes area based on the distribution in each zoogeographic zones. The division of zoogeographic zones were defined based on geologic history of the zones and areas of endemism which was defined based on previous studies [26, 46,47,48,49,50,51,52,53,54,55,56] (Fig. 1, modified from Holt et al. 2013).
Ancestral biogeography ranges estimation of Rhinolophoidea superfamily
For biogeographic regions of Rhinolophoidea superfamily, we divided the current known species distribution into nine zoogeographic areas using an updated version of Wallace’s zoogeographic regions of the world constructed on the distribution and phylogenetic relationship on 21,037 species of amphibians, birds and mammals [55]. The regions are shown in Figs. 1 and 5: Afrotropical + Madagascar (A); Oceanian (B); Philippine (C); Sino-Japanese (D); Oriental (including India + Indo-Malaya (Southeast Asia + Greater Sunda Islands (Sumatra, Java, Borneo)) (E); Palearctic (Europe, Mediterranean) (F); Sahara-Arabian (G); Australia (H); and Wallacea (Sulawesi, Lesser Sunda Islands and area in vicinity) (I). We assessed the biogeography of Rhinolophoidea across the tree, with 104 species at the tips (outgroup removed). We set the maximum area range to five (according to number of maximum areas occupied by extant taxa) resulting in 382 and 381 possible ranges with and without null-range, respectively.
Ancestral biogeography ranges estimation of Rhinolophidae family
Secondly, we pruned the Rhinolophus species from the tree of Asia lineages, consisting of 27 species sensu lato. The pruned tree did not include three species of R. hipposideros, R. formosae, R. luctus which are diverged early in the Oriental and Afrotropical lineages in the initial trees (Fig. 5). We attempted to prune Rhinolophidae as a whole family including all lineages. We included 47 species (OTUs) distributed in 13 biogeographic areas and 5 areas maximum occupied with extant taxa resulting in 2380 of possible geographic ranges. This caused computational issues, thus analysis was run separately for each lineage. We divided the analysis into Rhinolophus of AOO and the AMP subdivision. For the AOO clade, we divided the biogeographic regions into a smaller divisions (as shown in Fig. 2a) consisting of: Sino-Japanese (A); India (B); Indomalayan (C); Greater Sunda (Sumatra + Java + Borneo) (D); Philippine (E); Wallacea (Sulawesi + the Lesser Sunda Islands + Halmahera (from Wallace line to Lydekker line); Oceania (G); and Australia (H) [46, 47, 53, 55, 57, 58]. The total maximum area range was set to five, based on the number of areas occupied by extant species, (i.e. R. affinis with current distribution in India, Indomalaya, Greater Sunda, Philippines, and Wallacea), resulting in a total of 219 areas occupied.
In addition, we pruned Rhinolophidae in the AMP region consisting of 16 species, and number of maximum areas occupied by one species were set to 4, resulting 57 possible geographic ranges. The geographic areas included: Sahara-Arabian (A); Sudanian + Somalia + Ethiopia (B); Guinean + Congolian (C); Zambezian + Kalahari + Southearn Africa (D); Europe and Mediterranean (E); and Sino-Japanese (F) following biogeographic division of Africa, Palearctic and Sino-Japanese (as shown in Fig. 2b) [55, 59].
Ancestral biogeography ranges estimation of Hipposideridae and Rhinonycteridae family
We first ran Hipposideridae and Rhinonycteridae analysis across the entire old-world biogeographic region from 13 total biogeographic regions, of 57 species (OTUs) and seven areas that can be occupied by a single species. This calculated the total number of possible geographic ranges of 5812 states and created a large matrix (33,779,344 total) that could not be handled with current hardware and matrix-handling algorithms. Thus, we analyzed the AOO lineages separately from AMP lineages. For Hipposideridae AOO lineages, total of 39 species distributed in nine biogeographic regions were included in the analysis and a species can occupied a maximum of seven biogeographic areas, resulting in 502 possible geographic ranges. The biogeographic areas are shown in Fig. 3a including: Sudanian + Somalia + Ethiopia (A); Sino-Japanese (B); India (C); Indomalayan (D); Greater Sunda islands (Sumatra + Java + Borneo) (E); Philippine (F), Wallacean (G); Oceanian (H); and Australia (I). The biogeographic region of Sudanian + Somalia + Ethiopia (A) was being included in the Asian–Australian clades as H. jonesi appears to be sister to the Asian clades, and is likely to reflect an earlier African–Asian colonization event [60].
We pruned the AMP branches of Hipposideridae and the sister family Rhinonycteridae [14, 15]. The total number of species for this subset analysis is 18 distributed in six biogeographic regions, and five total area occupied by single extant species (numbers of possible ranges = 63). The biogeographic areas are shown in Fig. 3b including: Sahara-Arabian (A); Sudanian + Somalia + Ethiopia (B); Guinean + Congolian (C); Zambezian + Kalahari + Southern Africa (D), Madagascar (E); and Australia (F).
Ancestral biogeography ranges estimation of Cryptic Rhinolophidae in Asia
Additionally, we use the maximum clade credibility (MCC) tree from BEAST analysis (see below) to estimate the ancestral range of potential cryptic rhinolophids species in Asia (reviewed and delineated in [6]). In this subset analyses, we seek further detail of how the ancestors of recently diverged lineages expanded and contracted. The biogeographic units were split into eight “areas of endemism” of Southeast Asia and the West Pacific [47, 54, 61] including India (I), Himalaya (H); Southeast Asia (S); China (C); Isthmus of Kra region to Kangar-Pattani (K); and Kangar-Pattani line to the tip of Malay Peninsula (P); Borneo (B); and Sumatra (S) (as shown in Fig. 4). Maximum range size in analysis were set to eight (to match the number of regions) based on the assumption that OTUs (species) could occur in all areas. These areas represent the biogeographic zones used in former biogeographic analyses [47, 54, 61]. Area assignments for each OTU were based on current species distribution known from sequenced individuals which were cross-referenced with each biogeographic region they were recorded in ArcMap 10.3 to determine the range each OTU occurred in. The OTUs as tip nodes represent the species and cryptic species assigned from phylogenetic tree based on Maximum Likelihood and Bayesian Inferences.
For all above analyses, we use R package GenSA [62] (Generalized Simulated Annealing) to optimize the maximum likelihood calculation for all models in BioGeoBears. All BioGeoBEARS supermodels run under non-time-stratified analysis. We use statistical model comparison to compare the best-fit model given the data, a likelihood ratio test (LRT) for nested models and Akaike Information Criterion corrected (AICc) and weighted (AICw) were used for the non-nested models to observe the best model among all biogeographic scenarios.
Evolutionary history of potential cryptic Rhinolophidae species in Asia
To assess the evolutionary relationship and to estimate the time divergences between potential cryptic species, we ran BEAST v2.6.3 [63]. Potentially cryptic species in the region were delineated using integrative taxonomic approaches by combined phenotypic, acoustic and genetic data (as detailed in Chornelia et al. 2022) [6]. The tree was constructed based on 26 Rhinolophidae species senso lato distributed in Southeast Asia and India, using mtDNA COI 680 bp. The sequences were acquired from this study and mined from GenBank databases and from our previous work in the region (GenBank accession number: OK483366; OK483495–OK483509;OK562850–OK563036; OK563109; and OK563727) (see Chornelia et al. 2022 for all details of samples). Inferring evolutionary history using a single genes can be problematic for a variety of reasons, nonetheless, this is the available data in this region with higher systematic coverage [6]. Therefore, a care is needed in inferring the result and broad comparison with previous studies is crucial to provide an accurate understanding of species shifting ranges. Alignments were conducted in MAFFT using G-INS-i strategies, the output (in FASTA file) then converted into NEXUS file using Mesquite v3.6 [64], prior setting in.xml file was generated in BEAUti v2.6.3 [63] and SSM (Standard Substitution Models) package was loaded in BEAUti prior the analysis,. Relaxed clock log normal [65] was used to allow the clock rate to vary across the tree branches. We used a Birth–Death prior in the tree model which provided better accuracy and provides more precise result in all speciation scenarios based on previous studies [66]. The secondary calibration dates were taken from published papers due to the scarcity of bat fossils. Here we; (1) estimated the time divergence between Rhinolophidae and Hipposideridae was between 39 and 45 Ma, based on previous estimates including 41 Ma (95% highest posterior density (HPD) = 37–47 Ma) [22]; 39 Ma (95% HPD = 37–43 Ma) [67]; 42 Ma (95% HPD = 39–45 Ma) [14]. (2) To estimate time divergence between closely related species, we use Foley et al. (2015) for node calibration between R. shameli and R. creaghi (4 Ma (95% HPD = 3–5 Ma); R. trifoliatus-R. luctus (3 Ma (95% HPD = 2–4 Ma) and in addition, we calibrated node between H. armiger and A. stoliczkanus (31 Ma (26.5–31.5 Ma)). Thus, we set the priors to calibrate the nodes in BEAUti with parameterization as follows: (1) For Rhinolophidae and Hipposideridae, we selected H. armiger and R. sedulus, with prior distribution means (M) = 42 and standard deviations (S) = 0.045 which specifies that the distribution of priors is centred at 42 Ma and 95% probability range covering at 39–45 Ma; (2) R. shameli and R. creaghi prior, M = 4, S = 0.15 (median = 3.92 Ma. 95% probability range = 3.09–5.06 Ma); (3) R. trifoliatus and R. luctus prior, M = 3, S = 0.23 (median = 2.92 Ma, 95% probability range 2–4.27 Ma); (4) H. armiger and A. stoliczkanus prior, M = 30, S = 0.6 (median = 29.9 Ma, 95% probability range 27.1–33.1 Ma). All calibration priors are set to log-normal, each configuration of 2.5% quantile and 97.5% quantile set in ‘mean in real space’ and none of the prior settings enforce to be monophyletic as each “species” consist of several clades (as many of these are species complexes). MCMC algorithms were set to run for 50,000,000 cycles, trees stored every 1000 cycles. The trace files generated from BEAST were analyzed in TRACER1.7 [68]. 10% of initial tree were discarded in TreeAnnotator v2.6.3 [63]. The final tree was visualized in FigTree and finalized the graphics in Adobe Illustrator. The maximum clade credibility tree from the output was used for BioGeoBEARS analysis.
Availability of data and materials
All data used in this study are provided in this document.
Abbreviations
- AOO:
-
Australia–Oriental–Oceanian
- AMP:
-
Afrotropical–Madagascar–Palearctic
- BioGeoBEARS:
-
Biogeography with Bayesian (and likelihood) evolutionary analysis in R script
- Ma:
-
Million of years ago
- GBIF:
-
Global Biodiversity Information Facility
- DEC:
-
Dispersal–extinction cladogenetic
- d :
-
Dispersal
- e :
-
Extinction
- j :
-
Jump-dispersal or founder-event
- BEAST:
-
Bayesian Evolutionary Analysis by Sampling Trees
- GenSA:
-
Generalized simulated annealing
- LRT:
-
Likelihood Ratio Test
- AICc, AICc_wt:
-
Akaike Informatio Criterion corrected and weighted
- HPD:
-
Highest posterior density
References
Simmons NB, Cirranello AL. Bat Species of the World: a taxonomic and geographic database. https://batnames.org/. 2021. https://batnames.org/. Accessed 1 Jun 2021.
Anderson SC, Ruxton GD. The evolution of flight in bats: a novel hypothesis. Mamm Rev. 2020;50:426–39.
Amador LI, Almeida FC, Giannini NP. Evolution of traditional aerodynamic variables in bats (Mammalia: Chiroptera) within a comprehensive phylogenetic framework. J Mamm Evol. 2020;27:549–61.
Amador LI, Moyers Arévalo RL, Almeida FC, Catalano SA, Giannini NP. Bat systematics in the light of unconstrained analyses of a comprehensive molecular supermatrix. J Mamm Evol. 2018;25:37–70.
Csorba G, Ujhelyi P, Thomas N. Horseshoe bats of the world (Chiroptera: Rhinolophidae). Shropshire: Alana Books; 2003.
Chornelia A, Lu J, Hughes AC. How to accurately delineate morphologically conserved taxa and diagnose their phenotypic disparities: species delimitation in cryptic Rhinolophidae (Chiroptera). Front Ecol Evol. 2022;10:1–18.
Vanderelst D, Reijniers J, Firzlaff U, Peremans H. Dominant glint based prey localization in horseshoe bats: a possible strategy for noise rejection. PLoS Comput Biol. 2011;7(12):e1002268.
Vanderelst D, Reijniers J, Steckel J, Peremans H. Information generated by the moving pinnae of rhinolophus rouxi: tuning of the morphology at different harmonics. PLoS ONE. 2011;6(6):e20627.
Reijniers J, Vanderelst D, Peremans H. Morphology-induced information transfer in bat sonar. Phys Rev Lett. 2010;105:1–4.
Zhang Z, Truong SN, Müller R. Acoustic effects accurately predict an extreme case of biological morphology. Phys Rev Lett. 2009;103:1–4.
Zhuang Q, Müller R. Noseleaf furrows in a horseshoe bat act as resonance cavities shaping the biosonar beam. Phys Rev Lett. 2006;97(21):218701.
Vanderelst D, Lee YF, Geipel I, Kalko EKV, Kuo YM, Peremans H. The noseleaf of Rhinolophus formosae focuses the frequency modulated (FM) component of the calls. Front Physiol. 2013;4:1–9.
Kingston T, Lara MC, Jones G, Akbar Z, Kunz TH, Schneider CJ. Acoustic divergence in two cryptic Hipposideros species: a role for social selection ? Proc R Soc Lond. 2001;268:1381–6.
Foley NM, Thong VD, Soisook P, Goodman SM, Armstrong KN, Jacobs DS, et al. How and why overcome the impediments to resolution: lessons from rhinolophid and hipposiderid bats. Mol Biol Evol. 2015;32:313–33.
Armstrong KN, Goodman SM, Benda P, Hand SJ. A common name for the bat family Rhinonycteridae-the Trident Bats. Zootaxa. 2016;4179:115–7.
Thabah A, Rossiter SJ, Kingston T, Zhang S, Parsons S, Mya KMYA, et al. Genetic divergence and echolocation call frequency in cryptic species of Hipposideros larvatus s l (Chiroptera : Hipposideridae) from the Indo-Malayan region. Biol J Linn Soc. 2006;88(1):119–30.
Vallo P, Guillén-Servent A, Benda P, Pires DB, Koubek P. Variation of mitochondrial DNA in the Hipposideros caffer complex (Chiroptera: Hipposideridae) and its taxonomic implications. Acta Chiropterol. 2008;10:193–206.
Stoffberg S, Jacobs DS, Mackie IJ, Matthee CA. Molecular phylogenetics and historical biogeography of Rhinolophus bats. Mol Phylogenet Evol. 2010;54:1–9.
Teeling EC, Springer MS, Madsen O, Bates P, O’Brien SJ, Murphy WJ. A molecular phylogeny for bats illuminates biogeography and the fossil record. Science (80-). 2005;307:580–4.
Dool SE, Puechmaille SJ, Foley NM, Allegrini B, Bastian A, Mutumi GL, et al. Nuclear introns outperform mitochondrial DNA in inter-specific phylogenetic reconstruction: lessons from horseshoe bats (Rhinolophidae: Chiroptera). Mol Phylogenet Evol. 2016;97:196–212.
Demos TC, Webala PW, Goodman SM, Kerbis Peterhans JC, Bartonjo M, Patterson BD. Molecular phylogenetics of the African horseshoe bats (Chiroptera: Rhinolophidae): expanded geographic and taxonomic sampling of the Afrotropics. BMC Evol Biol. 2019;19:1–14.
Eick GN, Jacobs DS, Matthee CA. A nuclear DNA phylogenetic perspective on the evolution of echolocation and historical biogeography of extant bats (Chiroptera). Mol Biol Evol. 2005;22:1869–86.
Teeling EC. A molecular phylogeny for bats illuminates biogeography and the fossil record a molecular phylogeny for bats illuminates biogeography and the fossil record. Science (80-). 2005;307:580–4.
Simmons NB. Evolution. An Eocene big bang for bats. Science. 2005;307:527–8.
Jones K, Bininda-Emonds O, Gittleman J. Bats, clocks, and rocks: diversification patterns in chiroptera. Evolution (N Y). 2005;59:2243.
Sanmartín I. Historical biogeography: evolution in time and space. Evol Educ Outreach. 2012;5:555–68.
Crisci JV, Katinas L. Darwin, historical biogeography, and the importance of overcoming binary opposites. J Biogeogr. 2009;36:1027–32.
Crisp MD, Trewick SA, Cook LG. Hypothesis testing in biogeography. Trends Ecol Evol. 2011;26:66–72.
Guillen-Servent A, Francis CM, Ricklefs RE. Phylogeography and bioegeography of the horseshoe bats (Chiroptera: Rhinolpohidae). In: Csorba G, Ujhelyi P, Thomas N, editors. Horseshoe bats of the World (Chiroptera: Rhinolophidae). Shropshire: Alana Books; 2003. p. xii–xxii.
Bogdanowicz W, Owen RD. Phylogenetic analyses of the bat family Rhinolophidae. J Zool Syst Evol Res. 1992;30:142–60.
Norberg UM, Rayner JMV. Ecological morphology and flight in bats (Mammalia; Chiroptera): wing adaptations, flight performance, foraging strategy and echolocation. Philos Trans R Soc Lond B. 1987;316:335–427.
Luo B, Santana SE, Pang Y, Wang M, Xiao Y, Feng J. Wing morphology predicts geographic range size in vespertilionid bats. Sci Rep. 2019;9:1–6.
Moyers Arévalo RL, Amador LI, Almeida FC, Giannini NP. Evolution of body mass in bats: insights from a large supermatrix phylogeny. J Mamm Evol. 2020;27:123–38.
Norberg UM, Rayner JMV. Ecological morphology and flight in bats (mammalia: chiroptera): wing adaptations, flight performance, foraging strategy and echolocation. Philos Trans R Soc Lond B. 1987;316:335–427.
Burland TM, Wilmer JW, Burland TM, Wilmer JW. Seeing in the dark: molecular approaches to the study of bat populations. Biol Rev. 2007;76:389–409.
Matzke NJ. Model selection in historical biogeography reveals that founder-event speciation is a crucial process in island clades. Syst Biol. 2014;63:951–70.
Landis MJ, Matzke NJ, Moore BR, Huelsenbeck JP. Bayesian analysis of biogeography when the number of areas is large. Syst Biol. 2013;62:789–804.
Massana KA, Beaulieu JM, Matzke NJ, O’Meara BC. Non-null effects of the null range in biogeographic models: exploring parameter estimation in the DEC model. bioRxiv. 2015. https://doi.org/10.1101/026914.
Matzke NJ. Probabilistic historical biogeography: new models for founder-event speciation, imperfect detection, and fossils allow improved accuracy and model-testing. Front Biogeogr. 2013. https://doi.org/10.21425/F55419694.
Heine C, Yeo LG, Müller RD. Evaluating global paleoshoreline models for the Cretaceous and Cenozoic. Aust J Earth Sci. 2015;62:275–87.
Ree RH, Smith SA. Maximum likelihood inference of geographic range evolution by dispersal, local extinction, and cladogenesis. Syst Biol. 2008;57:4–14.
Ronquist F. Dispersal-vicariance analysis: a new approach to the quantification of historical biogeography. Syst Biol. 1997;46:195–203.
Álvarez-Carretero S, Tamuri AU, Battini M, Nascimento FF, Carlisle E, Asher RJ, et al. A species-level timeline of mammal evolution integrating phylogenomic data. Nature. 2021. https://doi.org/10.1038/s41586-021-04341-1.
Revell LJ. phytools: an R package for phylogenetic comparative biology (and other things). Methods Ecol Evol. 2012;3:217–23.
Hughes AC, Orr MC, Ma K, Costello MJ, Waller J, Provoost P, et al. Sampling biases shape our view of the natural world. Ecography (Cop). 2021;44:1259–69.
Holloway JD, Hall R. SE Asian geology and biogeography: an introduction. Leiden: Backhuys Publishers; 1996.
Turner H, Hovenkamp P, Van Welzen PC. Biogeography of Southeast Asia and the West Pacific. J Biogeogr. 2001;28:217–30.
Huggett R. Biogeographical regions. Chichester: eLS; 2011. p. 1–9.
Parnell J. The biogeography of the Isthmus of Kra region: a review. Nord J Bot. 2013;31:001–15.
Hall R. Reconstructing Cenozoic SE Asia. Geol Soc Spec Publ. 1996;106:153–84.
Sieh K, Natawidjaja D. Neotectonics of the Sumatran fault, Indonesia. J Geophys Res. 2000;105:28295.
Stauffer PH. Malaya and Southeast Asia in the pattern of continental drift. Bull Geol Soc Malaysia. 1974;7:89–138.
De Bruyn M, Stelbrink B, Morley RJ, Hall R, Carvalho GR, Cannon CH, et al. Borneo and Indochina are major evolutionary hotspots for Southeast Asian biodiversity. Syst Biol. 2014;63:879–901.
Ung V, Zaragueta-Bagils R, Williams DM. Comparative biogeography of Southeast Asia and the West Pacific region. Biol J Linn Soc. 2016;117:372–85.
Holt BG, Lessard J-P, Borregaard MK, Fritz SA, Araujo MB, Dimitrov D, et al. An update of Walace’s zoogeographic regions of the world. Science (80-). 2013;339:74–9.
Woodruff DS, Turner LM. The Indochinese-Sundaic zoogeographic transition: a description and analysis of terrestrial mammal species distributions. J Biogeogr. 2009;36:803–21.
Lohman DJ, de Bruyn M, Page T, von Rintelen K, Hall R, Ng PKL, et al. Biogeography of the Indo-Australian Archipelago. Annu Rev Ecol Evol Syst. 2011;42:205–26.
Hall R. Southeast Asia’s changing palaeogeography. Blumea J Plant Taxon Plant Geogr. 2009;54:148–61.
Linder HP, de Klerk HM, Born J, Burgess ND, Fjeldså J, Rahbek C. The partitioning of Africa: statistically defined biogeographical regions in sub-Saharan Africa. J Biogeogr. 2012;39:1189–205.
Patterson BD, Webala PW, Lavery TH, Agwanda BR, Goodman SM, Peterhans JCK, et al. Evolutionary relationships and population genetics of the afrotropical leaf-nosed bats (Chiroptera, Hipposideridae). Zookeys. 2020;2020:117–61.
Sheldon FH, Lim HC, Moyle RG. Return to the Malay Archipelago: the biogeography of Sundaic rainforest birds. J Ornithol. 2015;156:90–113.
Xiang Y, Gubian S, Suomela B, Hoeng J. Generalized simulated annealing for efficient global optimization: the GenSA package for R. R J. 2013;5:2013.
Bouckaert R, Vaughan TG, Barido-Sottani J, Duchêne S, Fourment M, Gavryushkina A, et al. BEAST 2.5: an advanced software platform for Bayesian evolutionary analysis. PLoS Comput Biol. 2019;15:1–28.
Maddison W, Maddison DR. Mesquite: a modular system for evolutionary analysis. Version 3.70. 2021. http://www.mesquiteproject.org.
Drummond AJ, Ho SYW, Phillips MJ, Rambaut A. Relaxed phylogenetics and dating with confidence. PLoS Biol. 2006;4:699–710.
Ritchie AM, Lo N, Ho SYW. The impact of the tree prior on molecular dating of data sets containing a mixture of inter- and intraspecies sampling. Syst Biol. 2017;66:413–25.
Teeling EC, Springer MS, Madsen O. A molecular phylogeny for bats illuminates biogeography and the fossil record. Science (80-). 2005;307:580–4.
Rambaut A, Drummond AJ, Xie D, Baele G, Suchard MA. Posterior summarization in Bayesian phylogenetics using Tracer 1.7. Syst Biol. 2018;67:901–4.
Soisook P, Struebig MJ, Noerfahmy S, Bernard H, Maryanto I, Chen SF, et al. Description of a new species of the Rhinolophus trifoliatus-group (Chiroptera: Rhinolophidae) from Southeast Asia. Acta Chiropterol. 2015;17:21–36.
Francis CM, Borisenko AV, Ivanova NV, Eger JL, Lim BK, Guillén-Servent A, et al. The role of DNA barcodes in understanding and conservation of mammal diversity in Southeast Asia. PLoS ONE. 2010;5:1–12.
Ree RH, Sanmartín I. Conceptual and statistical problems with the DEC+J model of founder-event speciation and its comparison with DEC via model selection. J Biogeogr. 2018;45:741–9.
Matzke NJ. Statistical comparison of DEC and DEC+J is identical to comparison of two ClaSSE submodels, and is therefore valid. OSF Prepr. 2021; April 26.
Ravel A, Marivaux L, Qi T, Wang YQ, Beard KC. New chiropterans from the middle Eocene of Shanghuang (Jiangsu Province, Coastal China): new insight into the dawn horseshoe bats (Rhinolophidae) in Asia. Zool Scr. 2014;43:1–23.
Maree S, Grant WS. Origins of HORSESHOE BATS ( Rhinolophus, Rhinolophidae ) in Southern Africa: evidence from allozyme variability. J Mamm Evol. 1997;4:195–215.
Ravel A, Adaci M, Bensalah M, Charruault A-L, Essid EM, Ammar HK, et al. Origine et radiation initiale des chauves-souris modernes: nouvelles découvertes dans l’Éocène d’Afrique du Nord. Geodiversitas. 2016;38:355–434.
Gunnell GF, Simmons NB. Fossil evidence and the origin of bats. J Mamm Evol. 2005;12:209–46.
Simmons NB, Seymour KL, Habersetzer J, Gunnell GF. Primitive early Eocene bat from Wyoming and the evolution of flight and echolocation. Nature. 2008;451:818–21.
Springer MS, Meredith RW, Janecka JE, Murphy WJ. The historical biogeography of mammalia. Philos Trans R Soc B. 2011;366:2478–502.
Ravel A, Marivaux L, Tabuce R, Adaci M, Mahboubi M, Mebrouk F, et al. The oldest African bat from the early Eocene of El Kohol (Algeria). Naturwissenschaften. 2011;98:397–405.
Servent AG, Francis CM, Ricklefs RE. Phylogeny and Biogeography of the Horseshoe Bats. In: Horseshoe bats of the World (Chiroptera: Rhinolophidae). 2003. p. xii–xxiv.
Rögl F. Palaeogeographic considerations for Mediterranean and Paratethys Seaways (Oligocene to Miocene). Ann des Naturhistorischen Museums Wien. 1998;99:279–310.
Zhang Z, Ramstein G, Schuster M, Li C, Contoux C, Yan Q. Aridification of the Sahara desert caused by Tethys Sea shrinkage during the late Miocene. Nature. 2014;513:401–4.
Ricou LE, Mercier De Lepinay B, Marcoux J. Evolution of the tethyan seaways and implications for the oceanic circulation around the Eocene–oligocene boundary. Dev Palaeontol Stratigr. 1986;9:387–94.
Rögl F. Mediterranean and paratethys. Facts and hypotheses of an Oligocene to Miocene paleogegoraphy (short overview). Cambridge: Cambridge University Press; 1999. p. 339–49.
Fleming TH, Racey P, Woodruff DS, Turner LM, Wilting A, Sollmann R, Meijaard E, et al. The Indochinese-Sundaic zoogeographic transition: a description and analysis of terrestrial mammal species distributions. J Biogeogr. 2009;42:803–21.
Hall R. Late Jurassic-Cenozoic reconstructions of the Indonesian region and the Indian Ocean. Tectonophysics. 2012;570–571:1–41.
Zhang L, Sun K, Feng J. Complete mitochondrial genome of the big-eared horseshoe bat Rhinolophus macrotis (Chiroptera, Rhinolophidae) Mitochondrial DNA Part A DNA Mapping. Seq Anal. 2016;27:4078–9.
Dool SE, Puechmaille SJ, Dietz C, Juste J, Ibáñez C, Hulva P, et al. Phylogeography and postglacial recolonization of Europe by Rhinolophus hipposideros: evidence from multiple genetic markers. Mol Ecol. 2013;22:4055–70.
Flanders J, Wei L, Rossiter SJ, Zhang S. Identifying the effects of the Pleistocene on the greater horseshoe bat, Rhinolophus ferrumequinum, in East Asia using ecological niche modelling and phylogenetic analyses. J Biogeogr. 2011;38:439–52.
Ohdachi S, Ishibashi Y, Iwasa M, Fukui D. The Wild Mammals of Japan. 2nd ed. Kyoto: Shoukadoh Book Sellers and The Mammals Society of Japan; 2015. p. 2015.
Ohnishi N, Uno R, Ishibashi Y, Tamate HB, Oi T. The influence of climatic oscillations during the Quaternary Era on the genetic structure of Asian black bears in Japan. Heredity (Edinb). 2009;102:579–89.
Hall M, Scotia N, Zealand N. Reconstructing recent divergence: evaluating nonequilibrium population structure in New Zealand chinook salmon. Mol Ecol. 2002;11:739–54.
Mitchell AHG, Hernandez F, dela Cruz AP. Cenozoic evolution of the Philippine archipelago. J Southeast Asian Earth Sci. 1986;1:3–22.
Voris HK. Maps of Pleistocene sea levels in SoutheastAsia: shorelines, river systems and time durations. J Biogeogr. 2000;27:1153–67.
Heaney LR. Biogeography of mammals in SE Asia: estimates of rates of colonization, extinction and speciation. Biol J Linn Soc. 1986;28:127–65.
Hall R. Cenozoic geological and plate tectonic evolution of SE Asia and the SW Pacific: computer-based reconstructions, model and animations. J Asian Earth Sci. 2002;20:353–431.
Woodruff DS. Biogeography and conservation in Southeast Asia: how 2.7 million years of repeated environmental fluctuations affect today’s patterns and the future of the remaining refugial-phase biodiversity. Biodivers Conserv. 2010;19:919–41.
Esselstyn JA, Brown RM. The role of repeated sea-level fluctuations in the generation of shrew (Soricidae: Crocidura) diversity in the Philippine Archipelago. Mol Phylogenet Evol. 2009;53:171–81.
Heaney LR, Robert T. New perspectives on the long-term biogeographic dynamics and conservation of Philippine fruit bats. In: Fleming TH, Racey PA, editors. Island bats: evolution, ecology and conservation. Chichago: The University of Chicago Press; 2009.
Heaney LR. Zoogeographic evidence for Middle and Late Pleistocene land bridges to the Philippine Islands. Mod Quarternary Res South East Asia. 1985;9:127–43.
Blackburn DC, Bickford DP, Diesmos AC, Iskandar DT, Brown RM. An ancient origin for the enigmatic flat-headed frogs (Bombinatoridae: Barbourula) from the Islands of Southeast Asia. PLoS ONE. 2010;5:e12090.
Siler CD, Oaks JR, Welton LJ, Linkem CW, Swab JC, Diesmos AC, et al. Did geckos ride the Palawan raft to the Philippines? J Biogeogr. 2012;39:1217–34.
Hughes M, Rubite RR, Blanc P, Chung K-F, Peng C-I. The Miocene to Pleistocene colonization of the Philippine archipelago by Begonia sect. Baryandra (Begoniaceae). Am J Bot. 2015;102:695–706.
Roberts TE, Museum F, Shore SL. Multiple levels of allopatric divergence in the endemic Philippine fruit bat Haplonycteris fischeri (Pteropodidae). Biol J Linn Soc. 2006;88:329–49.
Schmitt L, Kitchener D, How R. A genetic perspective of mammalian variation and evolution in the Indonesian archipelago: biogeographic correlates in the fruit bats genus Cynopterus. Evolution (N Y). 1995;49:399–412.
Jones AW, Kennedy RS. Evolution in a tropical archipelago: comparative phylogeography of Philippine fauna and flora reveals complex patterns of colonization and diversification. Biol J Linn Soc. 2008;95:620–39.
Schmitt L, Hisheh S, Suyanto A, Maharadatunkamsi, Newbound C, Kitchener D, et al. Crossing the line: the impact of contemporary and historical sea barriers on the population structure of bats in Southern Wallacea. In: Fleming TH, Racey PA, editors., et al., Island bats: evolution, ecology and conservation. Fleming. Chichago: The University of Chicago Press; 2009.
Mao XG, Zhu GJ, Zhang S, Rossiter SJ. Pleistocene climatic cycling drives intra-specific diversification in the intermediate horseshoe bat (Rhinolophus affinis) in Southern China. Mol Ecol. 2010;19:2754–69.
Wurster CM, Bird MI, Bull ID, Creed F, Bryant C, Dungait JAJ, et al. Forest contraction in north equatorial Southeast Asia during the last glacial period. Proc Natl Acad Sci USA. 2010;107:15508–11.
Hughes AC, Satasook C, Bates PJJ, Bumrungsri S, Jones G. Explaining the causes of the zoogeographic transition around the Isthmus of Kra: using bats as a case study. J Biogeogr. 2011;38:2362–72.
Hope G, Kershaw AP, van der Kaars S, Xiangjun S, Liew PM, Heusser LE, et al. History of vegetation and habitat change in the Austral-Asian region. Quat Int. 2004;118–119:103–26.
Patou ML, Wilting A, Gaubert P, Esselstyn JA, Cruaud C, Jennings AP, et al. Evolutionary history of the Paradoxurus palm civets—a new model for Asian biogeography. J Biogeogr. 2010;37:2077–97.
Morgan K, O’Loughlin SM, Chen B, Linton YM, Thongwat D, Somboon P, et al. Comparative phylogeography reveals a shared impact of pleistocene environmental change in shaping genetic diversity within nine Anopheles mosquito species across the Indo-Burma biodiversity hotspot. Mol Ecol. 2011;20:4533–49.
Meijaard E. Solving mammalian riddles along the Indochinese-Sundaic zoogeographic transition: new insights from mammalian biogeography. J Biogeogr. 2009;36:801–2.
Wiles GJ, Brooke AP, Wiles GJ, Brooke AP. Conservation threats to bats in the tropical Pacific Islands and insular Southeast Asia. In: Fleming TH, Racey PA, editors. Island bats: evolution, ecology and conservation. Chicago: THe University of Chicago Press; 2009. p. 2009.
Janis CM. Tertiary mammal evolution in the context of changing climates, vegetation, and tectonic events. Annu Rev Ecol Syst. 1993;24:467–500.
Bosworth W, Huchon P, McClay K. The Red Sea and Gulf of Aden Basins. J Afr Earth Sci. 2005;43:334–78.
Pirouz M, Avouac JP, Hassanzadeh J, Kirschvink JL, Bahroudi A. Early Neogene foreland of the Zagros, implications for the initial closure of the Neo-Tethys and kinematics of crustal shortening. Earth Planet Sci Lett. 2017;477:168–82.
Tamar K, Metallinou M, Wilms T, Schmitz A, Crochet PA, Geniez P, et al. Evolutionary history of spiny-tailed lizards (Agamidae: Uromastyx) from the Saharo-Arabian region. Zool Scr. 2018;47:159–73.
Yuan Z-Y, Zhang B-L, Raxworthy CJ, Weisrrock DW, Hime PM, Jin J-Q, et al. Natatanuran frogs used the Indian plate to step-stone disperse and radiate across the Indian Ocean. Natl Sci Rev. 2018;0:1–5.
Georgalis GL, Villa A, Delfino M. First description of a fossil chamaeleonid from greece and its relevance for the european biogeographic history of the group. Sci Nat. 2016. https://doi.org/10.1007/s00114-016-1336-5.
Toussaint EFA, Chiba H, Yago M, Dexter KM, Warren AD, Storer C, et al. Afrotropics on the wing: phylogenomics and historical biogeography of awl and policeman skippers. Syst Entomol. 2021;46:172–85.
Sahoo RK, Warren AD, Collins SC, Kodandaramaiah U. Hostplant change and paleoclimatic events explain diversification shifts in skipper butterflies (Family: Hesperiidae). BMC Evol Biol. 2017;17:1–9.
Pyron RA. Biogeographic analysis reveals ancient continental vicariance and recent oceanic dispersal in amphibians. Syst Biol. 2014;63:779–97.
Rabinowitz PD, Coffin MF, Falvey D. The separation of Madagascar and Africa. Science (80-). 1983;220:67–9.
Ali JR, Krause DW. Late Cretaceous bioconnections between Indo-Madagascar and Antarctica: Refutation of the Gunnerus Ridge causeway hypothesis. J Biogeogr. 2011;38:1855–72.
Yoder AD, Nowak MD. Has vicariance or dispersal been the predominant biogeographic force in Madagascar? Only time will tell. Annu Rev Ecol Evol Syst. 2006;37:405–31.
Crottini A, Madsen O, Poux C, Strauß A, Vieites DR, Vences M. Vertebrate time-tree elucidates the biogeographic pattern of a major biotic change around the K–T boundary in Madagascar. Proc Natl Acad Sci USA. 2012;109:5358–63.
Wood HM, Gillespie RG, Griswold CE, Wainwright PC. Why is Madagascar special? the extraordinarily slow evolution of pelican spiders (Aranaeae, Archaeidae). Evolution (N Y). 2015;69:462–81.
Tsang SM, Wiantoro S, Veluz MJ, Sugita N, Nguyen YL, Simmons NB, et al. Dispersal out of Wallacea spurs diversification of Pteropus flying foxes, the world’s largest bats (Mammalia: Chiroptera). J Biogeogr. 2020;47:527–37.
Weigelt P, Daniel Kissling W, Kisel Y, Fritz SA, Karger DN, Kessler M, et al. Global patterns and drivers of phylogenetic structure in island floras. Sci Rep. 2015;5:1–13.
Lavery TH, Leung LKP, Seddon JM. Molecular phylogeny of hipposiderid bats (Chiroptera: Hipposideridae) from Solomon Islands and Cape York Peninsula, Australia. Zool Scr. 2014;43:429–42.
Sen S. Dispersal of African mammals in Eurasia during the Cenozoic: ways and whys. Geobios. 2013;46:159–72.
Jønsson KA, Fabre PH, Ricklefs RE, Fjeldså J. Major global radiation of corvoid birds originated in the proto-Papuan archipelago. Proc Natl Acad Sci USA. 2011;108:2328–33.
van Ufford AQ, Cloos M. Cenozoic tectonics of New Guinea. Am Assoc Pet Geol Bull. 2005;89:119–40.
Schweizer M, Wright TF, Peñalba JV, Schirtzinger EE, Joseph L. Molecular phylogenetics suggests a New Guinean origin and frequent episodes of founder-event speciation in the nectarivorous lories and lorikeets (Aves: Psittaciformes). Mol Phylogenet Evol. 2015;90:34–48.
Hand SJ. On the winds of fortune. Aust Nat Hist. 1989;23:131–8.
Hand SJ. New Miocene leaf-nosed bats (Microchiroptera: HIpposideridae) from Riversleigh, Northwestern Queensland. In: Memoirs of the Queensland Museum. 2nd edition. Brisbane: Queensland Museum; 1997. p. 335–49.
Esselstyn JA, Evans BJ, Sedlock JL, Ali F, Khan A, Heaney LR. Single-locus species delimitation: a test of the mixed Yule—coalescent model, with an empirical application to Philippine round-leaf bats. Proc R Soc B. 2012. https://doi.org/10.1098/rspb.2012.0705.
Barbolini N, Woutersen A, Dupont-Nivet G, Silvestro D, Tardif D, Coster PMC, et al. Cenozoic evolution of the steppe-desert biome in Central Asia. Sci Adv. 2020;6:1–17.
De Bruyn M, Rüber L, Nylinder S, Stelbrink B, Lovejoy NR, Lavoué S, et al. Paleo-drainage basin connectivity predicts evolutionary relationships across three southeast asian biodiversity hotspots. Syst Biol. 2013;62:398–410.
Acknowledgements
We thank many researchers who made this study possible through their detailed studies of bats, mammals and biogeography and uploading their sequences to GenBank, and publicly shared the published data used in this study. We thank for the input of all reviewers on this manuscript. We especially thanks to Xing Yiaowu for suggestions, Jianmei Lu, Ade Prasetyo Agung, Thilina Nimalrathna, Supparat Tongkok and Arisa Kaewmano for assistance.
Funding
This project was supported by the Chinese National Natural Science Foundation (Grant No. U1602265, Mapping Karst Biodiversity in Yunnan), the Strategic Priority Research Program of the Chinese Academy of Sciences (Grant No. XDA20050202), and the Chinese Academy of Sciences Southeast Asia Biodiversity Research Center fund (Grant No. Y4ZK111B01).
Author information
Authors and Affiliations
Contributions
AC and ACH conceived the study, edited and revised the manuscript. AC compiled the data, conducted the analyses and drafted the initial draft of the manuscript, ACH compiled and processed the fossil data. Both of authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
This article does not involved any studies with human participants, and no animals were collected in this study. Field studies were conducted in accordance of local legislation and the study permit and ethics was reviewed and approved by Xishuangbanna Tropical Botanical Garden Ethics Committee.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests and the research conducted in the absence of any commercial or financial relationship that could be construed as a potential conflict of interest.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Additional file 1.
Supplementary figures.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Chornelia, A., Hughes, A.C. The evolutionary history and ancestral biogeographic range estimation of old-world Rhinolophidae and Hipposideridae (Chiroptera). BMC Ecol Evo 22, 112 (2022). https://doi.org/10.1186/s12862-022-02066-x
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12862-022-02066-x