Historical biogeography of New World emballonurid bats (tribe Diclidurini): taxon pulse diversification

Aim To reconstruct the biogeographical history of New World emballonurid bats (tribe Diclidurini). Although bats are the second most species‐rich order of mammals, they have not contributed substantially to our understanding of the historical biogeography of mammals in the Neotropics because of a poor fossil record. In addition, being the only group of mammals that fly, bats typically have large distributions with relatively few species endemic to restricted areas that are amenable to vicariant biogeographical approaches. Location Central and South America. Methods Phylogenetic analysis for comparing trees (PACT) is a new algorithm that incorporates all spatial information from taxon area cladograms into a general area cladogram. There were nine biogeographical areas identified in Central and South America for New World emballonurid bats. Molecular dating was used to incorporate the temporal aspect of historical biogeography. This method was compared with dispersal–vicariance analysis (DIVA), which assumes vicariance as the default mode of speciation. Results Of the 45 speciation events in a fully resolved phylogeny, eight that were hypothesized by DIVA as vicariance were considered by PACT as two peripheral isolations and six within‐area events. DIVA was less parsimonious because it required six more post‐speciation dispersal events in addition to the 73 hypothesized by PACT. DIVA reconstructed a widely distributed ancestor, suggesting that most dispersal events occurred earlier, whereas the ancestral area for PACT based on character optimization was the Northern Amazon, suggesting that dispersal events were more recent phenomena. Main conclusions The general area cladogram from PACT indicated that within‐area events, and not vicariance, provide the major mode of speciation for New World emballonurid bats. There was no biological evidence supporting or rejecting sympatric speciation in New World emballonurid bats. However, the geological history, combined with fluctuations in temperature and sea level, suggested within‐area speciation in a changing and heterogeneous environment in the Northern Amazon during the Miocene. This scenario is similar to the taxon‐pulse hypothesis of biotic diversification, which posits repeated episodes of range expansions and contractions from a stable core area such as the Guiana Shield within the Northern Amazon.     

PACT in practice: comparative historical biogeographic patterns and species–area relationships of the Greater Antillean and Hawaiian Island terrestrial biotas

Aim To compare the evolutionary and ecological patterns of two extensively studied island biotas with differing geological histories (the Hawaiian Islands and the Greater Antilles). We evaluated the results from PACT (phylogenetic analysis for comparing trees), an innovative approach that has been proposed to reveal general patterns of biotic expansion (between regions) and in situ (within a region) diversification, as well as species–area relationships (SAR) and the taxon pulse dynamic. Location The Hawaiian Islands and Greater Antilles. Methods We used the PACT algorithm to construct general area cladograms and identified biotic expansion and in situ nodes. We analysed the power‐law SAR and relative contribution of biotic expansion and in situ diversification events using power‐law and linear regression analyses. Results Both biotic expansion and in situ nodes were prevalent throughout the PACT general area cladograms (Greater Antilles, 55.9% biotic expansion, 44.1% in situ; Hawaiian Islands, 40.6% biotic expansion, 59.4% in situ). Of the biotic expansion events, both forward and backward events occurred in both regions (Greater Antilles, 85.1% forward, 14.9% backward; Hawaiian Islands, 65% forward, 35% backward). Additionally, there is a power‐law SAR for the Greater Antilles but not for the Hawaiian Islands. However, exclusion of Hawai'i (the youngest, largest Hawaiian Island) produced a power‐law SAR for the Hawaiian Islands. Main conclusions The prevalence of in situ events as well as forward and backward biotic expansion events reveals that both Hawaiian and Greater Antillean biotas have evolved through alternating episodes of biotic expansion and in situ diversification. These patterns are characteristic of the taxon pulse dynamic, for which few data have previously been recorded on islands. Additionally, our analysis revealed that historical influences on the power‐law SARs are pronounced in both assemblages: old, small islands are relatively species rich and young, large islands are relatively species poor. Thus, our PACT results are consistent with hypotheses of geological influence on the evolution of island biotas and also provide greater insight into the role of the taxon pulse dynamic in the formation of island equilibria.     

Time and space in biogeography: response to Parenti & Ebach (2013)

A recent Guest Editorial by Parenti & Ebach (2013, Journal of Biogeography, 40, 813–820) disagrees with the methods or interpretations in two of our recent papers. In addition, the authors open a debate on biogeographical concepts, and present an alternative philosophy for biogeographical research in the context of their recently described biogeographical subregion called ‘Pandora’. We disagree with their approach and conclusions, and comment on several issues related to our differing conceptual approaches for biogeographical research; namely, our use of molecular phylogenetic analyses, including time estimates; and Parenti & Ebach's reliance on taxon/general area cladograms. Finally, we re‐examine their ‘tests’ supporting the existence of ‘Pandora’.     

Cladistic and phylogenetic biogeography: the art and the science of discovery

All methods used in historical biogeographical analysis aim to obtain resolved area cladograms that represent historical relationships among areas in which monophyletic groups of taxa are distributed. When neither widespread nor sympatric taxa are present in the distribution of a monophyletic group, all methods obtain the same resolved area cladogram that conforms to a simple vicariance scenario. In most cases, however, the distribution of monophyletic groups of taxa is not that simple. A priori and a posteriori methods of historical biogeography differ in the way in which they deal with widespread and sympatric taxa. A posteriori methods are empirically superior to a priori methods, as they provide a more parsimonious accounting of the input data, do not eliminate or modify input data, and do not suffer from internal inconsistencies in implementation. When factual errors are corrected, the exemplar presented by M.C. Ebach & C.J. Humphries (Journal of Biogeography, 2002, 29 , 427) purporting to show inconsistencies in implementation by a posteriori methods actually corroborates the opposite. The rationale for preferring a priori methods thus corresponds to ontological rather than to epistemological considerations. We herein identify two different research programmes, cladistic biogeography (associated with a priori methods) and phylogenetic biogeography (associated with a posteriori methods). The aim of cladistic biogeography is to fit all elements of all taxon–area cladograms to a single set of area relationships, maintaining historical singularity of areas. The aim of phylogenetic biogeography is to document, most parsimoniously, the geographical context of speciation events. The recent contribution by M.C. Ebach & C.J. Humphries (Journal of Biogeography, 2002, 29 , 427) makes it clear that cladistic biogeography using a priori methods is an inductivist/verificationist research programme, whereas phylogenetic biogeography is hypothetico‐deductivist/falsificationist. Cladistic biogeography can become hypothetic‐deductive by using a posteriori methods of analysis.     

Escaping the matrix: a new algorithm for phylogenetic comparative studies of co-evolution

An algorithm for generating host cladograms from parasite‐host cladograms derived from parasite phylogenies, Phylogenetic Analysis for Comparing Trees (PACT), is described. PACT satisfies Assumption 0, that all the information in each parasite‐host cladogram must be used in a co‐evolutionary analysis, and that the host relationships depicted in the final host cladogram must be logically consistent with the phylogenetic relationships depicted in every part of every parasite‐host cladogram used to construct the host cladogram. It accounts for cases of speciation by host switching and expansion of host range, and reticulated host relationships, in addition to co‐speciation, sympatric speciation, and extinction in all input parasite‐host cladograms, and does so without a priori weighting schemes and without a posteriori manipulation of the data.     

Gauging the effects of sampling failure in biogeographical analysis

Aim Various methods are employed to recover patterns of area relationships in extinct and extant clades. The fidelity of these patterns can be adversely affected by sampling error in the form of missing data. Here we use simulation studies to evaluate the sensitivity of an analytical biogeographical method, namely tree reconciliation analysis (TRA), to this form of sampling failure. Location Simulation study. Methods To approximate varying degrees of taxonomic sampling failure within phylogenies varying in size and in redundancy of biogeographical signal, we applied sequential pruning protocols to artificial taxon–area cladograms displaying congruent patterns of area relationships. Initial trials assumed equal probability of sampling failure among all areas. Additional trials assigned weighted probabilities to each of the areas in order to explore the effects of uneven geographical sampling. Pruned taxon–area cladograms were then analysed with TRA to determine if the optimal area cladograms recovered match the original biogeographical signal, or if they represent false, ambiguous or uninformative signals. Results The results indicate a period of consistently accurate recovery of the true biogeographical signal, followed by a nonlinear decrease in signal recovery as more taxa are pruned. At high levels of sampling failure, false biogeographical signals are more likely to be recovered than the true signal. However, randomization testing for statistical significance greatly decreases the chance of accepting false signals. The primary inflection of the signal recovery curve, and its steepness and slope depend upon taxon–area cladogram size and area redundancy, as well as on the evenness of sampling. Uneven sampling across geographical areas is found to have serious deleterious effects on TRA, with the accuracy of recovery of biogeographical signal varying by an order of magnitude or more across different sampling regimes. Main conclusions These simulations reiterate the importance of taxon sampling in biogeographical analysis, and attest to the importance of considering geographical, as well as overall, sampling failure when interpreting the robustness of biogeographical signals. In addition to randomization testing for significance, we suggest the use of randomized sequential taxon deletions and the construction of signal decay curves as a means to assess the robustness of biogeographical signals for empirical data sets.     

A historical biogeographical protocol for studying biotic diversification by taxon pulses

Aim To present a historical biogeographical protocol for distinguishing biotic diversification by taxon pulse radiations from biotic diversification by vicariance. Location Mexico and northern Central America. Methods Brooks Parsimony Analysis (BPA), phylogenetic inference, linear correlation analysis. Results The taxon pulse radiation of 33 clades in nine areas of endemism in Mesoamerica is based on nine episodes of biotic expansion from three areas, and six episodes of vicariance, involving four geographical splits. Nineteen per cent of speciation events are due to vicariance, 25% to peripheral isolates speciation and 56% are within‐area events. The species–area curve has a correlation coefficient (r²) of 0.47. Extinction events and species richness are highly correlated (r² = 0.75), but colonization events and species richness are poorly correlated (r² = 0.36), suggesting that colonization is not the main determinant of the species–area relationship. Colonization events are more poorly correlated with size of area (r² = 0.05) than are in situ speciation events (r² = 0.60). Colonization events and in situ events are poorly correlated (r² = 0.02). All areas of endemism have reticulated histories, and have acted as both sources and islands at various times. Main conclusions Taxon pulses can be distinguished from maximum vicariance using this protocol; refining it requires a method for generating area cladograms from complex data and incorporation of direct dating of evolutionary events.     

Parsimony analysis of endemicity (PAE) of Mexican terrestrial mammals at different area units: when size matters

Abstract Aim Parsimony analysis of endemicity (PAE) is a biogeographical method that uses a parsimony algorithm to obtain an area cladogram, based on taxa inhabiting the study areas. We compare its performance at different geographical units (½° and 1° quadrats, ecoregions and biogeographical provinces) to analyse distributional patterns of Mexican terrestrial mammals, in order to assess the importance of the size of area units. Location The area analysed corresponds to Mexico. Methods Parsimony analyses were based on 56,859 collection records, corresponding to 703 genera, species and subspecies. Four data matrices were constructed for: (1) 716 quadrats of ½° latitude × ½° longitude, (2) 230 quadrats of 1° latitude × 1° longitude, (3) forty‐five ecoregions and (4) fourteen biogeographical provinces. Results For the ½° quadrat matrix, we obtained six cladograms of 17,138 steps. For the 1° quadrat matrix, we obtained five cladograms (strict consensus with 9394 steps). For the matrix of ecoregions, we obtained twelve cladograms (strict consensus cladogram with 3009 steps). For the provinces, we obtained a single cladogram with 1603 steps. Main conclusions The best results were obtained with natural areas instead of quadrats. There seems to exist a trend to decrease the absolute number of steps and an increase in the absolute and relative number of synapomorphies as the size of the area units decreases, although this does not necessarily occur for the number of cladograms.     

Do the Oaxacan Highlands represent a natural biotic unit? A cladistic biogeographical test based on vertebrate taxa

Aim  We analysed the distributional patterns of six terrestrial vertebrate taxa from the Oaxacan Highlands (Sierra Mazateca, Nudo de Zempoaltépetl and Sierra de Juárez) through a cladistic biogeographical approach, in order to test their naturalness as a biotic unit.
Location  The Oaxacan Highlands, Mexico.
Methods  The cladistic biogeographical analysis was based on the area cladograms of the Pseudoeurycea bellii species group (Amphibia: Plethodontidae), the genus Chlorospingus (Aves: Thraupidae), the genera Microtus , Reithrodontomys and Habromys , and the Peromyscus aztecus species group (Mammalia: Rodentia). We obtained paralogy-free subtrees, from which the components were coded in a data matrix for parsimony analysis. The data matrix was analysed with N ona through W in C lada .
Results  The parsimony analysis resulted in a single general area cladogram in which areas were fragmented following the sequence Sierra Madre Occidental, Trans-Mexican Volcanic Belt, Chiapas, Sierra Madre Oriental + Sierra Mazateca, Sierra Madre del Sur, Nudo de Zempoaltépetl and Sierra de Juárez.
Main conclusions  The general area cladogram shows that the Oaxacan Highlands do not constitute a natural unit. The Sierra Mazateca is the sister area to the Sierra Madre Oriental, whereas the Nudo de Zempoaltépetl and the Sierra de Juárez are closely related to the Sierra Madre del Sur. The events that might have caused these patterns include cycles of expansion and contraction of mountain pinyon, juniper and oak woodlands during the Pleistocene.     

Dispersal and diversification: macroevolutionary implications of the MacArthur–Wilson model, illustrated by Simulium (Inseliellum) Rubstov (Diptera: Simuliidae)

Aim Provide an empirical test of the ‘radiation zone’ hypothesis of the MacArthur–Wilson theory of island biogeography using the taxon‐pulse hypothesis of Erwin and Brooks Parsimony Analysis (BPA) on Simulium (Inseliellum) Rubstov. Location Micronesia, Cook Islands, Austral Islands, Society Islands, Marquesas Islands, Fiji and New Caledonia. Methods Primary and secondary BPA of the phylogeny of Inseliellum. Results Primary BPA showed that 15% of the taxon area cladogram contained area reticulations. Secondary BPA (invoking the area duplication convention) generated a clear sequence of dispersal for Inseliellum. The sequence follows a Micronesia – Cook Islands – Marquesas Islands – Society Islands dispersal, with a separate dispersal from the Cook Islands to the Austral Islands less than 1 Ma. A radiation in the island of Tahiti (Society Islands) produced numerous dispersals from Tahiti to other islands within the Society Islands system. Islands close to Tahiti (source island) have been colonized from Tahiti more often than islands far from Tahiti, but a higher proportion of those species colonizing distant islands have become distinct species. Main conclusions The dispersal sequence of Inseliellum exhibits both old to young island dispersal and young to old island dispersal. This is due to habitat availability on each island. Inseliellum is a model system in exemplifying the ‘radiation zone’ hypothesis of MacArthur and Wilson. As well, islands close to the source are colonized more often that those far from the source, but colonization of islands far away from the source results in a higher proportion of speciation events than for islands close to the source. The diversification of Inseliellum corresponds to a taxon‐pulse radiation, with a centre of diversification on Tahiti resulting from its large area and abundant freshwater habitats. This study illustrates the utility of BPA in identifying complex scenarios that can be used to test theories about the complementary roles of ecology and phylogeny in historical biogeography.     

The historical biogeography of co-evolution: emerging infectious diseases are evolutionary accidents waiting to happen

Ecological fitting refers to interspecific associations characterized by ecologically specialized, yet phylogenetically conservative, resource utilization. During periods of biotic expansion, parasites and hosts may disperse from their areas of origin. In conjunction with ecological fitting, this sets the stage for host switching without evolving novel host utilization capabilities. This is the evolutionary basis of emerging infectious diseases (EIDs). Phylogenetic analysis for comparing trees (PACT) is a method developed to delineate both general and unique historically reticulated and non‐reticulated relationships among species and geographical areas, or among parasites and their hosts. PACT is based on ‘Assumption 0’, which states that all species and all hosts in each input phylogeny must be analysed without modification, and the final analysis must be logically consistent with all input data. Assumption 0 will be violated whenever a host or area has a reticulated history with respect to its parasites or species. PACT includes a Duplication Rule, by which hosts or areas are listed for each co‐evolutionary or biogeographical event affecting them, which satisfies Assumption 0 even if there are reticulations. PACT maximizes the search for general patterns by using Ockam's Razor – duplicate only enough to satisfy Assumption 0. PACT applied to the host and geographical distributions of members of two groups of parasitic helminths infecting anthropoid primates indicates a long and continuous association with those hosts. Nonetheless, c. 30% of the host associations are due to host switching. Only one of those involves non‐primate hosts, suggesting that most were constrained by resource requirements that are phylogenetically conservative among primates (ecological fitting). In addition, most of the host switches were associated with episodes of biotic expansion, also as predicted by the ecological fitting view of EIDs.     

Input data, analytical methods and biogeography of Elegia (Restionaceae)

Aim The aim of this paper is to determine the optimal methods for delimiting areas of endemism for Elegia L. (Restionaceae), an endemic genus of the Cape Floristic Region. We assess two methods of scoring the data (presence–absence in regular grids, or in irregular eco‐geographical regions) and three methods for locating biogeographical centres or areas of endemism, and evaluate one method for locating biotic elements. Location The Cape Floristic Region (CFR), South Africa. Methods The distribution of all 48 species of Elegia was mapped as presence–absence data on a quarter‐degree grid and on broad habitat units (eco‐geographical areas). Three methods to delimit areas of endemism were applied: parsimony analysis of endemism (PAE), phenetic cluster analysis, and NDM (‘end em ism’). In addition, we used presence–absence clustering (‘Prabclus’) to delimit biotic elements. The performances of these methods in elucidating the geographical patterns in Elegia were compared, for both types of input data, by evaluating their efficacy in maximizing the proportion of endemics and the number of areas of endemism. Results Eco‐geographical areas perform better than quarter‐degree grids. The eco‐geographical areas are potentially more likely to track the distribution of species. The phenetic approach performed best in terms of its ability to delimit areas of endemism in the study area. The species richness and the richness of range‐restricted species are each highest in the south‐western part of the CFR, decreasing to the north and east. The phytogeographical centres identified in the present study are the northern mountains, the southern mountains (inclusive of the Riviersonderend Mountains and the Cape Peninsula), the Langeberg range, the south coast, the Cape flats, and the west coast. Main conclusions This study demonstrates that (1) eco‐geographical areas should be preferred over a grid overlay in the study of biogeographical patterns, (2) phenetic clustering is the most suitable analytical method for finding areas of endemism, and (3) delimiting biotic elements does not contribute to an understanding of the biogeographical pattern in Elegia. The areas of endemism in Elegia are largely similar to those described in other studies, but there are many detailed differences.     

Application of parsimony analysis of endemicity to Mexican gymnosperm distributions: grid-cells, biogeographical provinces and track analysis

Parsimony analysis of endemicity (PAE) was used to analyse the distributional patterns of 124 species of Mexican gymnosperms, using two different sample units: grid-cells and biogeographical provinces. PAE analyses were based on distributional data from herbarium specimens and specialized literature. Two data matrices were constructed for 60 grid-cells of 2° and 14 biogeographical provinces. The analysis of the 2° grid-cell matrix led to 7084 cladograms. The strict consensus cladogram showed several clades equivalent to the results obtained with the biogeographical provinces. Three clades agree with some principal regions of distribution of Mexican pines, previously identified by several authors, located at the northern portion of the Baja California peninsula, the Sierra Madre Occidental, and the Sierra Madre Oriental. These areas represent important centres of species diversity and endemism for Mexican gymnosperms. The analysis of the province matrix led to two most parsimonious cladograms, which only differed in the position of the Sierra Madre Occidental province. The iterative procedure PAE with progressive character elimination was applied to identify generalized tracks, where clades of provinces were considered equivalent to generalized tracks, and each time a cladogram was obtained, species defining its clades were deleted and a new run was undertaken. We found five generalized tracks, mainly located in montane provinces. The distribution patterns of gymnosperms agree with the existence of several Mexican biogeographical provinces, and a different historical biogeography of the Mexican peninsulas from the rest of the country is evident.  © 2007 The Linnean Society of London, Biological Journal of the Linnean Society , 2007, 92 , 405–417.     

Dispersal vs. vicariance in the Mediterranean: historical biogeography of the Palearctic Pachydeminae (Coleoptera, Scarabaeoidea)

Aim The geological evolution of the Mediterranean region is largely the result of the Tertiary collision of the African and Eurasian Plates, but also a mosaic of migrating island arcs, fragmenting tectonic belts, and extending back‐arc basins. Such complex paleogeography has resulted in a ‘reticulate’ biogeographical history, in which Mediterranean biotas repeatedly fragmented and merged as dispersal barriers appeared and disappeared through time. In this study, dispersal‐vicariance analysis (DIVA) is used to assess the relative role played by dispersal and vicariance in shaping distribution patterns in the beetle subfamily Pachydeminae Reitter, 1902 (Scarabaeoidea), an example of east–west Mediterranean disjunction. Location The Mediterranean region, including North Africa, the western Mediterranean, Balkans–Anatolia, Middle East, Caucasus, the Iranian Plateau, and Central Asia. Methods A phylogenetic hypothesis of the Palearctic genera of Pachydeminae in conjunction with distributional data was analysed using DIVA. This method reconstructs the ancestral distribution in a given phylogeny based on the vicariance model, while allowing dispersal and extinction to occur. Unlike other methods, DIVA does not enforce area relationships to conform to a hierarchical ‘area cladogram’, so it can be used to reconstruct ‘reticulate’ biogeographical scenarios. Results Optimal reconstructions, requiring 23 dispersal events, suggest that the ancestor of Pachydeminae was originally present in the south‐east Mediterranean region. Basal splitting within the subfamily was caused by vicariance events related to the late Tertiary collision of the African microplates Apulia and Arabia with Eurasia, and the resultant arise of successive dispersal barriers (e.g. the Red Sea, the Zagros Mountains). Subsequent diversification in Pachydeminae involved multiple speciation events within the Middle East and Iran–Afghanistan regions, which gave rise to the least speciose genera of Pachydeminae (e.g. Otoclinius Brenske, 1896). Finally, the presence of Pachydeminae in the western Mediterranean region seems to be the result of a recent dispersal event. The ancestor of the Iberian genera Ceramida Baraud, 1987 and Elaphocera Gené, 1836 probably dispersed from the Middle East to the Iberian Peninsula across North Africa and the Gibraltar Strait during the ‘Messinian salinity crisis’ at the end of the Miocene. Main conclusions Although the basal diversification of Pachydeminae around the Mediterranean appears to be related to vicariance events linked to the geological formation of the Mediterranean Basin, dispersal has also played a very important role. Nearly 38% of the speciation events in the phylogeny resulted from dispersal to a new area followed by allopatric speciation between lineages. Relationships between western and eastern Mediterranean disjuncts are usually explained by dispersal through Central Europe. The biogeographical history of the Pachydeminae corroborates other biogeographical studies that consider North Africa to be an alternative dispersal route by which Mediterranean taxa could have achieved circum‐Mediterranean distributions.     

Island evolutionary biogeography: analysis of the weevils (Coleoptera: Curculionidae) of the Falkland Islands (Islas Malvinas)

Aim I analysed distributional and phylogenetic information on weevils (Coleoptera: Curculionidae) from the Falklands, and integrated it with molecular, palaeontological and geological information to infer a geobiotic scenario. Location Falkland Islands (Islas Malvinas). Methods The panbiogeographical analysis was based on data on 23 Falkland species and their related taxa from southern South America. For the cladistic biogeographical analysis I analysed six weevil taxa for which phylogenetic hypotheses are available (the generic groups Cylydrorhinus, Strangaliodes and Falklandius, and the genera Antarctobius, Germainiellus and Puranius). Results from this analysis were compared with previous regionalizations. Cenocrons (sets of taxa that share the same biogeographical history) were identified by considering temporal information provided by fossils and molecular clocks. Finally, a geobiotic scenario was proposed by integrating the available information. Results Six generalized tracks were detected: Maule–Valdivian forests, Magellanic forest, Magellanic moorland, Falkland Islands, Magellanic forest–Magellanic moorland, and Magellanic forest–Falkland Islands. A node was identified in the Magellanic forest, based on the overlap of two generalized tracks. A single general area cladogram was obtained, implying the following sequence: (Magellanic moorland (Maule–Valdivian forests (Magellanic forest, Falkland Islands))). The Falklands are classified here as a biogeographical province in the Austral realm, Andean region and Subantarctic subregion. Falkland weevils seem to belong to a single Subantarctic cenocron. The sequence of events deduced implies the following steps: development of the Subantarctic biota in southern South America, arrival of the Falkland crustal block from South Africa in the Early Cretaceous, geodispersal of the Subantarctic cenocron from southern South America to the Falklands during the Early Oligocene, vicariance of the Magellanic moorland, vicariance of the Maule–Valdivian forests, and final vicariance between the Magellanic forest and the Falkland Islands. Main conclusions The biotic components identified support the connection of the Falkland weevils with the Magellanic forest. Falkland weevils belong to a single cenocron, dated to at least the Early Oligocene, when geodispersal from southern South America may have occurred. An older African cenocron may have been replaced completely by the Subantarctic one when the proto‐Falklands made contact with the Patagonian continental shelf. A geobiotic scenario implying vicariance events related to sea‐level variations could explain the distributional patterns analysed herein.     

Around the world in 10 million years: biogeography of the nearly cosmopolitan true toads (Anura: Bufonidae)

Aim The species‐rich family of true toads (Anura: Bufonidae) has been the focus of several earlier studies investigating the biogeography of geographically widespread taxa. Herein, we employ newly developed Bayesian divergence estimate methods to investigate the biogeographical history of this group. Resulting age estimates are used to test several key temporal hypotheses including that the origin of the bufonid clade pre‐dates Gondwanan vicariance (~105 million years ago, Ma). Area cladograms are also invoked to investigate the geographical origin of the family. Location Worldwide, except the Australia–New Guinea plate, Madagascar and the Antarctic. Methods A phylogenetic hypothesis of the relationships among true toads was derived from analysis of 2521 bp of DNA data including fragments from three mitochondrial (12S, tRNA^val, 16S) and two nuclear (RAG‐1, CXCR‐4) genes. Analysis of multiple, unlinked loci with a Bayesian method for estimating divergence times allowed us to address the timing and biogeographical history of Bufonidae. Resulting divergence estimates permitted the investigation of alternative vicariance/dispersal scenarios that have been proposed for true toads. Results Our area cladogram resulting from phylogenetic analysis of DNA data supports a South American origin for Bufonidae. Divergence estimates indicate that the family originated earlier than had been suggested previously (78–99 Ma). The age of the enigmatic Caribbean clade was dated to the late Palaeocene–early Eocene. A return of bufonids to the New World in the Eocene was followed by rapid diversification and secondary expansion into South America by the early Oligocene (Rupelian). Main conclusions The South American origin of Bufonidae in the Upper Cretaceous was followed by relatively rapid expansion and radiation around the globe, ending with a return to the Americas via a Eurasian/North American land bridge in the Eocene. Though the exact route of this dispersal (Beringia or North Atlantic) remains unclear, an argument is made for the less frequently invoked North Atlantic connection. The origin of the enigmatic Caribbean lineage was found to be consistent with colonization following the bolide impact at the K/T boundary. These findings provide the first, firm foundation for understanding true toad divergence times and their truly remarkable and global radiation.     

Spatial congruence between biotic history and species richness of Muscidae (Diptera,Insecta) in the Andean and Neotropical regions

Considering that Earth and life evolve together, the present study aims to verify whether the species richness patterns are spatially congruent to biotic history. Niche conservatism was adopted as a background hypothesis to associate species richness with phylogenetic information. A parallel analysis between this procedure and cladistic biogeography was undertaken. Eleven Muscidae genera that were previously systematically reviewed for phylogenetic hypotheses were chosen for the analysis. The genera were split into ‘basal’ and ‘derived’ species, following terminal taxon root distance within each genus. Richness patterns were contrasted for the most basal and most derived 33% of species, and richness maps were constructed at 220 × 220 km grid size. A difference richness map was drawn by derived minus basal values (=derived?basal). For regions with difference values around zero, a component analysis was performed and compared with relationships established by other studies. Derived and basal species richness showed a very concise richness gradient in the Neotropical region and it was compatible with its known biogeographical history. In the Andean region, richness did not show any pattern. The area cladogram grouped Subantarctic subregion in a polytomy and Central Chile as a paraphyletic group. All hypotheses about area relationship were divergent and no vicariant pattern could be recognized in Andean region. In Neotropical region, Muscidae results corroborated a previous component relationship. The hypothesis that Paleogene climatic changes could drive the biotic component’s split was suggested. In the Andean region, recently ice sheet covering events had driven the species to disperse and/or extinct resulted in absence of pattern seen either in richness analysis or in component analysis. It is believed that species richness is linked to biotic history and this fact may be considered when evaluating hypotheses to explain broad‐scale richness gradients.     

Historical biogeography of some river basins in central Mexico evidenced by their goodeine freshwater fishes: a preliminary hypothesis using secondary Brooks parsimony analysis

Aims Our aim was to uncover and describe patterns of historical biogeography of the main river basins in central Mexico, based on a secondary Brooks parsimony analysis (BPA) of goodeine fishes, and to understand the processes that determine them with respect to the molecular clock of the goodeines and the geological events that have taken place in the region since the Miocene. Location The region covered in this study includes central Mexico, mostly the so‐called Mesa Central of Mexico, an area argued to be a transitional zone comprising several major river drainages from their headwaters at high elevations along the Transmexican Volcanic Belt to the coast of the Gulf of Mexico and the Pacific Ocean. Methods Based on a previous phylogenetic hypothesis regarding the Goodeidae, we built a data matrix using additive binary coding. First, we conducted a primary BPA to provide general explanations of the historical biogeography of Central Mexico. As ambiguity was found, a secondary BPA was conducted, and some areas were duplicated in order to explain the reticulated history of the area. Area cladograms were obtained by running a parsimony analysis. Instances of vicariance and non‐vicariance processes were described with reference to the cladogram obtained from secondary BPA. Results The study area was divided into 18 discrete regions. Primary BPA produced nine equally parsimonious cladograms with 129 steps, and a consistency index (CI) of 0.574. A strict consensus cladogram shows low resolution among some areas, but other area relationships are consistent. For secondary BPA, five of the 18 regions were duplicated (LEA, COT, AYU, CUT, PAN); one was triplicated (BAL); and one was quadruplicated (AME), suggesting that the pattern of distribution of species in these areas reflects multiple independent events. These areas correspond with the regions exhibiting the highest levels of diversification and the most complex geological history, and those for which river piracy events or basin connections have been proposed. The secondary BPA produced a single most parsimonious cladogram with 118 steps, and a CI of 0.858. This cladogram shows that none of the duplicated areas are nested together, reinforcing the idea of a reticulated history of the areas and not a single vicariant event. Main conclusions Although our results are preliminary and we cannot establish this as a general pattern, as the BPA is based on a single‐taxon cladogram, resolution obtained in the secondary BPA provides some insights regarding the historical biogeography of this group of fishes in river basins of central Mexico. Secondary BPA indicates that the historical biogeography of central Mexico, as shown by their goodeine freshwater fishes, is complex and is a result of a series of vicariant and non‐vicariant events such as post‐dispersal speciation and post‐speciation dispersal.     

Biogeographical origins and diversification of the exoneurine allodapine bees of Australia (Hymenoptera,Apidae)

Aim Early diversification of allodapine bees occurred in Africa c. 50 Ma. They are most abundant in sub‐Saharan Africa and Australia, and one of the oldest phylogenetic divergences in the tribe involves a split between an African + Malagasy clade and an Australian clade. The historical biogeographical scenario for this has been highly problematic, entailing an Eocene dispersal from Africa to Australia, followed by an unresolved, and apparently rapid, set of bifurcations leading to the Australian ‘exoneurine’ genera. Here we use an expanded taxon set of Australian species to explore the timing and historical biogeography of the exoneurine radiation. Location Australia, Africa, Madagascar. Methods One nuclear gene (F2 copy of elongation factor 1α) and two mitochondrial genes (cytochrome c oxidase subunit I and cytochrome b) were sequenced for 33 Australian exoneurine species from all five genera found on the continent, as well as for an additional 37 species from all non‐parasitic genera in the remainder of the tribe. We used Bayesian inference analyses to study phylogenetic topology and penalized likelihood analyses to infer key dates of divergence within the tribe. We also used lineage‐through‐time (LTT) analyses and Bayesian analyses to explore the tempo of radiations and biogeographical history of the exoneurines. Results Results from the phylogenetic analyses were congruent with previous studies, indicating a single colonization event c. 34 Ma, too late for Gondwanan vicariance models, and too early for a Laurasian dispersal route. In contrast to earlier studies, we show that this colonization event did not result in an ancient rapid radiation. However, LTT patterns indicated a rapid radiation of the temperate‐adapted genera Exoneura and Brevineura, but not of the xeric‐adapted genus Exoneurella, from 10 to 6 Ma. Main conclusions Our results indicate a trans‐oceanic dispersal event from Africa to Australia, most likely via Antarctica, with an accelerated diversification of temperate‐adapted lineages during the major Late Miocene event referred to as the ‘Hill Gap’. This is the first study to link radiations in Australian bee faunal elements to changing climate, and differs from many other plant and insect phylogenetic studies by showing increased radiation of temperate clades, rather than xeric clades, with increasing aridification of Australia.