Skip to main content

Towards a synthesis of the Caribbean biogeography of terrestrial arthropods



The immense geologic and ecological complexity of the Caribbean has created a natural laboratory for interpreting when and how organisms disperse through time and space. However, competing hypotheses compounded with this complexity have resulted in a lack of unifying principles of biogeography for the region. Though new data concerning the timing of geologic events and dispersal events are emerging, powerful new analytical tools now allow for explicit hypothesis testing. Arthropods, with varying dispersal ability and high levels of endemism in the Caribbean, are an important, albeit understudied, biogeographic model system. Herein, we include a comprehensive analysis of every publicly available genetic dataset (at the time of writing) of terrestrial Caribbean arthropod groups using a statistically robust pipeline to explicitly test the current extent of biogeographic hypotheses for the region.


Our findings indicate several important biogeographic generalizations for the region: the South American continent is the predominant origin of Caribbean arthropod fauna; GAARlandia played a role for some taxa in aiding dispersal from South America to the Greater Antilles; founder event dispersal explains the majority of dispersal events by terrestrial arthropods, and distance between landmasses is important for dispersal; most dispersal events occurred via island hopping; there is evidence of ‘reverse’ dispersal from islands to the mainland; dispersal across the present-day Isthmus of Panama generally occurred prior to 3 mya; the Greater Antilles harbor more lineage diversity than the Lesser Antilles, and the larger Greater Antilles typically have greater lineage diversity than the smaller islands; basal Caribbean taxa are primarily distributed in the Greater Antilles, the basal-most being from Cuba, and derived taxa are mostly distributed in the Lesser Antilles; Jamaican taxa are usually endemic and monophyletic.


Given the diversity and deep history of terrestrial arthropods, incongruence of biogeographic patterns is expected, but focusing on both similarities and differences among divergent taxa with disparate life histories emphasizes the importance of particular qualities responsible for resulting diversification patterns. Furthermore, this study provides an analytical toolkit that can be used to guide researchers interested in answering questions pertaining to Caribbean biogeography using explicit hypothesis testing.


The Caribbean region (Fig. 1) holds a long and rich history of entomological research and discovery. After nearly two centuries of entomological studies [1], we are still collecting and describing many new species of terrestrial arthropods (e.g., [2], and references therein, [3,4,5,6,7]). The > 700 islets and islands of the Caribbean (~ 240,000 km2), their dramatic elevational gradients (− 39 m to + 3098 m), and their proximity to two continents (North and South America) have resulted in a hyperdiversity of arthropods that can be both a boon and a pitfall for research. Unlike island systems such as Hawaii and French Polynesia where the terrestrial arthropod fauna is well-known [8], the Caribbean fauna is more diverse and less well-known, hampering areas of research such as biogeography.

Fig. 1
figure 1

Map of the study region. The Caribbean Basin. Colors correspond to those used in subsequent figures of ancestral range estimations from BioGeoBEARS analyses. For information about geographic areas and geologic dates, see Additional file 1: Table S1. Base map created in ArcGIS v9.2 and areas shaded in Adobe Illustrator Creative Cloud. Red = South America (SA); Orange = Central America (CA); Yellow-Green = North America (NA); Purple = Florida (FL); Bright Pink = Bahamas and Turks and Caicos Islands (BA, TCI); Green = Greater Antilles (GA); Yellow = Jamaica (JA); Light Blue = Northern Lesser Antilles (NLA); Dark Blue = Southern Lesser Antilles (SLA); Blue-Green = Barbados (BAR)

The Caribbean Basin is a geologically complex region, bounded on the north, west, and south by continental landmasses. The Antillean Archipelago within the basin comprises islands with broad ages, size ranges, and elevational gradients. To the east is a volcanic arc comprising older volcanics as well as active volcanoes, and to the north are the Greater Antilles and the Bahamas Bank. While land may have been continuously available in the Greater Antilles for nearly 40 million years (my), islands such as those on the Bahama Bank are just over 100 thousand years old (kyo). Additionally, some islands, such as the Greater Antilles, may have been connected to one another and to the continental landmasses by land bridges (Wallacean islands), whereas the Lesser Antilles have never been connected to one another or the continents (de novo islands). For reviews of Caribbean geology, see [9,10,11,12,13] and references therein.

The primary biogeographic questions concerning Caribbean biota are how it arrived to the islands and where it is derived from. The ages of the present-day island biota have also been debated. For example, Hedges et al. [14] suggested that the bolide impact at the K-T boundary would have wiped out all organisms present on the islands, whereas Crother and Guyer [15] presented evidence to the contrary. However, geologic data [11, 12, 16] indicate that continuous land in the Greater Antilles was unavailable until after ~ 40 mya. More recent evidence about the region’s geology has shifted the discussion from whether organisms became established on islands via vicariance or dispersal to when organisms arrived and what happened after [17].

Competing hypotheses

In the last 25 years, new geologic and biologic evidence has allowed researchers to propound two somewhat controversial hypotheses (Fig. 2). The first being the GAARlandia (Greater Antilles-Aves Ridge) hypothesis of Iturralde-Vinent and MacPhee [11]. According to geologic and paleobiogeographic data, a landspan may have connected northern South America with the Greater Antilles 35–32 mya. This hypothesis has been invoked to explain the distribution of mammals [11], including bats [22, 23], as well as toads [24], spiders [25, 26], butterflies [27], fish [28], and plants [29, 30], although only two of these studies [25, 27] explicitly tested the GAARlandia hypothesis using statistical methods. Contrarily, some researchers (e.g., [31]) remain unconvinced that the Aves Ridge was subaerial based on the available geologic data and suggest the need for more sea-floor drilling to elucidate the extent of the land bridge or island chain.

Fig. 2
figure 2

Historical Caribbean landmasses and major competing hypotheses for dispersal. Schematic representations of the Caribbean land areas hypothesized to be available at certain time periods (see Additional file 1: Table S1 for geologic references) and the major competing hypotheses examined herein (see Table 1). Arrows indicate dispersal, dates next to the images depicting the Isthmus of Panama relate to the hypothesized Central American Seaway closure(s)

A second equally controversial hypothesis that could greatly impact the origins of Caribbean fauna is the timing of the Central American Seaway (CAS) closure. The prevailing paradigm from both biologic and geologic data has been that the Isthmus of Panama (IoP) was completely formed (the CAS was completely closed) by the Late Pliocene, approximately 3 mya, but that the entire process occurred gradually over tens of millions of years. O’Dea et al. [18, 32] and Leigh et al. [33] and references therein review the geologic and biologic evidence for CAS closure at various times. In general, collision between the southern tip of Central America and the South American landmass began 24–23 mya ([19] and references therein), and it is thought that most of the volcanic arc was submerged; however, there is evidence that indicates significant land may have been above sea level during the Mid Miocene (~ 15–13 mya) in the form of an archipelago or a peninsula. Coates and Stallard [20] suggest a closure at 15 mya, and Osborne et al. [21] suggest the CAS closed between 12 and 7 mya. Leigh et al. [33] suggest that by 10 mya the land bridge was nearly complete, allowing some biotic interchange, and by 7 mya, most deep water populations were separated. After O’Dea et al. [18] exhaustively reviewed all available data, they concluded final closure nearer to the Late Pliocene. Bacon et al. [34, 35] used biologic data to infer geologic hypotheses, concluding an older age for IoP formation/CAS closure.

Terrestrial arthropods in biogeographic modeling

The first study of Caribbean biogeography focused on vertebrates [36], though soon after, studies including terrestrial arthropods followed suit [37]. Arthropods are excellent study organisms for biogeography because they have relatively large populations and resilient life stages [38, 39]. Their resilience additionally means that there are often older lineages which can help clarify the effects of older geologic events [40]. Also, because similar patterns are often found between vertebrate and invertebrate taxa (e.g., [25, 41, 42]), such congruence only serves to reinforce the significance of the results (e.g., [43, 44]).

Great advances in our understanding of Caribbean arthropods have been made in the last 30 years. Liebherr [2] edited the first synthesis of Caribbean insect biogeography. Soon after, the first Caribbean biogeography studies using molecular data from insects began to emerge [45,46,47,48]. Over the past 10 years, numerous papers featuring terrestrial arthropods and focusing on historical biogeography of the Caribbean have been published (e.g., [25,26,27, 49,50,51,52]) as well as several phylogeographic studies (e.g., [41, 53, 54]). However, as noted by Rosen [55] and emphasized by Slater [56], “a stupendous multidisciplinary effort to resolve decisive patterns for the [Caribbean] region” is needed.

In this study we have examined all publicly available datasets of terrestrial arthropod taxa that met our criteria for selection: a taxon with a distribution across a significant part of the Caribbean Basin (> 75% land areas) with some endemic species or populations; mostly complete sampling of known lineages and distributions (> 75% across the distribution); genetic differentiation sufficient for phylogenetic analysis. Our analyses consist of ten taxa: Platythyrea ants, Heterotermes and Nasutitermes termites, butterflies from the genera Calisto and Papilio, Drosophila flies, centruroidine scorpions, and Micrathena, Spintharus, and Selenops spiders. These taxa provide a good representation of sedentary and more vagile taxa with various life histories. Our strategy was to maximize the sampling effort, including geographic, taxonomic, and phylogenetic coverage in order to test all of the prevailing hypotheses regarding the timing and route of dispersal into the Caribbean Archipelago. We time-calibrated phylogenies to explicitly test up to 252 biogeographic models per taxon in order to elucidate the overarching patterns of Caribbean biogeography for terrestrial arthropods.

All models are wrong, but some are useful

We would like to make a call to biologists to interact more with geologists when studying biogeography so that the maximum information may be gleaned from the geological evidence and unlikely or impossible conclusions aren’t invoked to explain the results. Using biology to infer specific geologic events should only be done with extreme caution (e.g. – [43, 44, 57]) as the work of researchers outside of a specific discipline can be easily misinterpreted and subsequently mis-cited repeatedly in the literature. Most Caribbean biogeographic research has tried to invoke support for or against the GAARlandia land bridge (35–32 mya) using diversification rates coupled with the timing of diversification events. While the authors acknowledge GAARlandia is a controversial hypothesis, in this study the models assume that GAARlandia existed as a geologic entity, and we have evaluated the extent to which it would have played a role in dispersal into the Caribbean from South America.

There is also ongoing debate [58] about whether there are fundamental problems with biogeographic models including the ‘jump dispersal’ parameter. However, for island taxa, jump dispersal is known to be very important [59]. Portions of the Caribbean archipelago were never connected or proximal to other islands (Darwinian islands); therefore, the only route to these islands is via jump dispersal. Following this logic, jump dispersal is expected to be common in island taxa, especially those that fly or balloon, and any potential bias introduced is insignificant compared to the bias introduced by removing the ‘jump’ parameter. To fully examine this effect, we have tested 54–252 models per dataset. Similarly, we have incorporated our current understanding of fossils, molecular clocks, and geology to obtain dated phylogenies for hypothesis testing considering the current ongoing debate [60].

A note about the term ‘colonization’

We would like to add a note about our omission of the term ‘colonization’ in this paper and make a plea to the biogeography community that this terminology be used with caution. The term ‘colonization’ has been used outside of the biology literature to refer to the arrival of European culture outside of Europe. The process of ‘colonization’ of the Americas (and many other parts of the world) by Europeans resulted in the death and enslavement of millions of people. This is particularly true in the Caribbean region. In an acknowledgement of the negative connotations associated with this term, and because none of the taxa we examined are invasive, we have chosen to omit it and instead discuss the process by which populations arrive to, and become established in, a new area by simply using the term ‘dispersal’.


For each arthropod taxon (Platythyrea ants, Heterotermes termites, Nasutitermes termites, Calisto butterflies, Papilio butterflies, Drosophila fruit flies, centruroidine scorpions, Micrathena spiders, Spintharus spiders, and Selenops spiders), we followed a pipeline of analyses, and the detailed results of all analyses can be found in Additional file 1. Briefly, for each, a phylogeny was constructed in MrBayes, and if the resulting topology was congruent with a published phylogeny for that taxon, then it was used as the constraint for dating in BEAST 2 and subsequent analyses in BioGeoBEARS (Papilio, centruroidine scorpions, Heterotermes). If the resulting topology was incongruent with the published phylogeny, both a topology constrained to the previously published topology and a topology constrained to our Bayesian consensus tree were used for all downstream dating and analyses (Platythyrea, Nasutitermes). In some cases, there was incongruence among published phylogenies (Drosophila), or irreconcilable differences between the published phylogeny and our own topology (Micrathena, Spintharus). For the former, we analyzed multiple phylogenies with various constraints based on previous studies, and for the latter, we used our Bayesian consensus tree for downstream analysis. Finally, some taxa had no published topology for which we could compare our results due to the addition of data, so the Bayesian consensus tree was used for all downstream dating and analysis (Selenops, Calisto).

Dating analyses for each taxon established minimum ages for the group. For those taxa that were younger than the proposed date of GAARlandia, the set of models which allow for the putative landspan were excluded from testing (Platythyrea, Nasutitermes, Papilio, Drosophila, Selenops – younger dated tree, see Table 2 and Additional file 1). Those taxa that were old enough to show any potential effect of GAARlandia were tested for the full suite of 252 models (Heterotermes, Centruroidinae, Micrathena, Spintharus, Selenops – older dated tree, see Table 2 and Additional file 1). Finally, Calisto is endemic to the Greater Antilles, so could only be analyzed under 18 models.

A total of 2282 biogeographic models were tested across all 10 taxa in BioGeoBEARS, and the results of biogeographic model testing for all taxa are summarized in Table 2. Resulting phylograms are shown in Figs. 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13 and 14, and phylograms representing probabilities as pie charts can be found in Additional file 1: Figures S1–S5 and S7–S13 along with tables listing included taxa and respective GenBank reference numbers for samples (Additional file 1: Tables S5, S7–S15). All input and output files for all analyses can be found on figshare. Results from individual taxa are detailed in Additional file 1. In a few instances, the BioGeoBEARS analyses for certain models simply would not run, even when we changed starting parameters, probably because the likelihoods were too low for that particular model (Additional file 1: Table S3 and Table S4).

Fig. 3
figure 3

Phylogeny and ancestral range estimation for Platythyrea ants. BioGeoBEARS phylogram corresponding to the DIVALIKE C1a model (Table 1). The tree shows the single most probable geographic range at each node pre- and post-split. Colors correspond to Fig. 1. SA = South America; CA = Central America; NA = North America; GA = Greater Antilles; SLA = Southern Lesser Antilles; NLA = Northern Lesser Antilles; BA = Bahamas; FL = Florida. Ant photo by Judy Gallagher

Fig. 4
figure 4

Phylogeny and ancestral range estimation for Heterotermes termites. BioGeoBEARS phylogram corresponding to the DIVALIKE A1a model (Table 1). The tree shows the single most probable geographical range at each node pre- and post-split. Colors correspond to Fig. 1. SA = South America; CA = Central America; NA = North America; GA = Greater Antilles; JA = Jamaica; SLA = Southern Lesser Antilles; NLA = Northern Lesser Antilles; TCI = Turks and Caicos Islands; FL = Florida. Heterotermes photo provided by Rudolf Scheffrahn

Fig. 5
figure 5

Phylogeny and ancestral range estimation for Nasutitermes termites. BioGeoBEARS phylogram corresponding to the DEC D1b model (Table 1). The tree shows the single most probable geographical range at each node pre- and post-split. Colors correspond to Fig. 1. SA = South America; CA = Central America; NA = North America; GA = Greater Antilles; JA = Jamaica; SLA = Southern Lesser Antilles; NLA = Northern Lesser Antilles; BA = Bahamas. Nasutitermes photo by Bernard Dupont

Fig. 6
figure 6

Phylogeny and ancestral range estimation for Calisto butterflies (Euptychia as the outgroup). BioGeoBEARS phylogram corresponding to the DEC + J A1b model (Table 1). The tree shows the single most probable geographical range at each node pre- and post-split. Colors correspond to Additional file 1: Figure S6. SA = South America; CU = Cuba; JA = Jamaica; HI = Hispaniola; PR = Puerto Rico; BA = Bahamas. Photo by S. Crews

Fig. 7
figure 7

Phylogeny and ancestral range estimation for Calisto butterflies (Euptychia is not the sister taxon). BioGeoBEARS phylogram corresponding to the DEC + J A2a model (B2a, D2a, E2a, and G2a being equiprobable as the best models based on the AICc weights) (Table 1). The tree shows the single most probable geographical range at each node pre- and post-split. Colors correspond to Additional file 1: Figure S6. SA = South America; NA = North America; CU = Cuba; JA = Jamaica; HI = Hispaniola; PR = Puerto Rico; BA = Bahamas. Photo by S. Crews

Fig. 8
figure 8

Phylogeny and ancestral range estimation for Papilio butterflies. BioGeoBEARS phylogram corresponding to the BAYAREALIKE + J C1a model (D1a is equally likely) (Table 1). The tree shows the single most probable geographical range at each node pre- and post-split. Colors correspond to Fig. 1. SA = South America; CA = Central America; NA = North America; GA = Greater Antilles; JA = Jamaica; SLA = Southern Lesser Antilles; NLA = Northern Lesser Antilles; FL = Florida; BA = Bahamas. Papilio photo by Eduardo Lopez

Fig. 9
figure 9

Phylogeny and ancestral range estimation for Drosophila flies. BioGeoBEARS phylogram corresponding to the BAYAREALIKE + J D1a model (Table 1) for the tree with the relationship of D. belladunni and D. acutilabella constrained (see Additional file 1). Asterisks indicate differences between the constrained and unconstrained trees in the ancestral range estimation. The tree shows the single most probable geographical range at each node pre- and post-split. Colors correspond to Fig. 1. SA = South America; CA = Central America; NA = North America; GA = Greater Antilles; JA = Jamaica; SLA = Southern Lesser Antilles; NLA = Northern Lesser Antilles; FL = Florida; BAR = Barbados. Drosophila photo by Mark Yokoyama

Fig. 10
figure 10

Phylogeny and ancestral range estimation from BioGeoBEARS for Drosophila flies without constraining the relationship of D. belladunni and D. acutilabella. The best model based on AICc weights is BAYAREALIKE + J D1a (Table 1). The tree shows the single most probable geographical range at each node pre- and post-split. Colors correspond to Fig. 1. SA = South America; CA = Central America; NA = North America; GA = Greater Antilles; JA = Jamaica; SLA = Southern Lesser Antilles; NLA = Northern Lesser Antilles; FL = Florida; BAR = Barbados. Drosophila photo by Mark Yokoyama

Fig. 11
figure 11

Phylogeny and ancestral range estimation for centruroidine scorpions. BioGeoBEARS phylogram corresponding to the BAYAREALIKE + J B1a and B2a models (Table 1). The tree shows the single most probable geographical range at each node pre- and post-split. Colors correspond to Fig. 1. SA = South America; CA = Central America; NA = North America; GA = Greater Antilles; JA = Jamaica; SLA = Southern Lesser Antilles; NLA = Northern Lesser Antilles; FL = Florida; BA = Bahamas. Photo by S. Crews

Fig. 12
figure 12

Phylogeny and ancestral range estimation for Micrathena spiders. BioGeoBEARS phylogram corresponding to the DIVALIKE + J B2a model (Table 1). The tree shows the single most probable geographical range at each node pre- and post-split. Colors correspond to Fig. 1. SA = South America; CA = Central America; NA = North America; GA = Greater Antilles; JA = Jamaica; SLA = Southern Lesser Antilles; NLA = Northern Lesser Antilles. Photo by S. Crews


Our analyses indicate some shared patterns and processes of dispersal to and subsequent diversification on islands in the Caribbean biodiversity hotspot, including dispersal across the IoP/CAS (Table 2). Below we discuss some of the major themes surrounding the evolution and biogeography in the Caribbean. These data provide a baseline for indicating the types of questions that can be confidently answered given the available data and highlight where knowledge is lacking and can be improved upon in future research. We base our discussion on the phylogeny that shows the single most probable ancestral range for each taxon; for trees that show the probabilities for each possible range, see Additional file 1.

Origins of Caribbean taxa

South America is the most common place of origin for the Caribbean fauna examined, followed by Central America and North America (Table 2). This is congruent with the generalized track analysis and cladistic biogeography of Morrone [61] and Echeverry and Morrone [62] for many arthropod (and other) taxa. There does not seem to be a shared evolutionary history, ecology, or life history among groups that share a point of origin, and this does not appear to coincide with the timing of dispersal or age of the group in question. Some taxa appear to have dispersed to the islands from more than one area of the mainland, including a relatively young group of good dispersers (Drosophila) [63], and a relatively old group of poor dispersers (Micrathena) [26] (Figs. 9, 10, 12). These results are congruent with other literature on Caribbean arthropods; for example, beetles are also thought to have dispersed to the islands from both Central America and South America [40, 64], and ants [65] and auchenorrhynchous homoptera [66] are thought to have dispersed to the Caribbean from both Central America and North America.

Number of dispersal events to Caribbean islands

We find that all but the poorest dispersers have arrived in the Caribbean multiple times (Table 2). The primary predictor of the number of dispersal events into the Caribbean appears to be dispersal ability, including mode of locomotion, as well as adaptability and ecological plasticity, which are important for establishment after dispersal. However, extinction events could also produce patterns of multiple dispersal events. It should be noted that the plot of most-probable areas is not the same as event counts [67]. Simple counts give some idea of event number but do not consider events occurring along a branch; for better estimates, biogeographic stochastic mapping is required.

Among flying insects, there is a general pattern of multiple dispersal events. Drosophila fruit flies, Papilio butterflies, and Nasutitermes termites each dispersed at least 4 times (Figs. 5, 8, 9 and 10), and Heterotermes termites and Calisto butterflies dispersed twice (Figs. 4, 6, 7). Drosophila is thought to be a good disperser [63], and the groups examined here are thought to have dispersed relatively recently (~ 4 my) to the Caribbean, perhaps indicating that dispersal ability has played a role in the number of dispersal events to the islands. Papilio’s dispersal abilities have been noted by Lewis et al. [52] who also found several separate dispersal events rather than a single dispersal followed by diversification (Fig. 8). In contrast, Calisto butterfly species have only dispersed one or two times (depending on the outgroup used) to the Greater Antilles from South America (Figs. 6, 7). Among the termites, Heterotermes dispersed at least 3 times from South America and North America (Fig. 4). However, according to our results, Nasutitermes is relatively new in the Caribbean (~ 15 my) but appears to be a good disperser (Fig. 5). In at least one species (N. corniger), Scheffrahn et al. [68] noted that kings and queens are quite good flyers and can adapt to many different types of habitats, from urban to more natural, and are good at taking over and defending territories.

Selenops spiders have dispersed to Caribbean islands only twice, from South America to the Lesser Antilles for one event, and to the Greater Antilles for the second event, although the point of origin remains unclear (Fig. 14). Within Centruroides scorpions, there are also two dispersal events, one from South America and one from North America (Fig. 11). Both Selenops spiders and Centruroides scorpions likely have a similar ecology as nocturnal somewhat stationary sit-and-wait predatory arachnids. Spintharus spiders have seemingly only dispersed to the Caribbean once, and although this species is thought to be a poor disperser [51], it does appear capable of dispersal to proximal islands (Fig. 13).

Fig. 13
figure 13

Phylogeny and ancestral range estimation for Spintharus spiders. BioGeoBEARS phylogram corresponding to the DIVALIKE + J A1b model (Table 1). The tree shows the single most probable geographical range at each node pre- and post-split. Colors correspond to Fig. 1. SA = South America; CA = Central America; NA = North America; GA = Greater Antilles; JA = Jamaica; SLA = Southern Lesser Antilles; NLA = Northern Lesser Antilles. Photo by Judy Gallagher

Somewhat puzzling are Micrathena spiders, considered to be a poor dispersersing taxon [26], they have dispersed to the Caribbean multiple times (6 in our analyses) (Fig. 12). One possible explanation is, as our analyses indicate, that they are a rather old group which would allow for more time in which dispersal events could occur. Missing taxa from the mainland could also produce a pattern of multiple dispersal events. Nevertheless, other studies with more taxa sampled also indicated multiple dispersal events [69], leading us to conclude that this taxon is in fact a rather good disperser like many other orb-weaving spiders.

Szalanski et al. [48] suggested that the present distribution of Heterotermes in the Caribbean is likely due to overwater dispersal, however he emphasized that anthropogenic introductions couldn’t be ruled out. In our analyses, it does appear that most species lineages are at least a few million years old (Fig. 4), but that some of the intraspecific distributions may be anthropogenically caused. Human-mediated dispersal, even intraspecific dispersal within the native range, could have a significant impact on biogeographic results. In some cases, this may be a very evident pattern which emerges from the data, however, there may be examples such as this where it is cryptic and confounds the conclusions that can be drawn about the taxon.

Islands as sources of diversity

Islands are not usually thought of as sources for diversity, but rather sinks [70, 71]. Recently, however, more examples of dispersal from an island source in the Caribbean has become more well-known, particularly in Anolis lizards, including from a large island to smaller islands or even from islands to the mainland [72, 73]. Dispersal from the islands to the mainland is also seen in weevils [74] and orioles [75]. One reason we may not be as aware of this phenomenon is because mainland sampling from Mexico, Central America, and South America is often poor and therefore unable to provide the phylogenetic resolution required to detect such events; even excluding a single species could change the results.

There are many examples from our analyses of dispersal from Cuba to the Bahamas, perhaps due to these land areas being particularly close to one another during times of low sea level. Although we only used a single specimen of Nasutitermes ripperti from the Bahamas in our analyses, this species is also found in Cuba, so perhaps dispersed from Cuba to the Bahamas. In Calisto butterflies, dispersal from Cuba to the Bahamas has occurred at least twice (Fig. 13). In Selenops spiders and Centruroides scorpions, taxa are shared between Cuba and the Bahamas and presumably dispersed from Cuba to the Bahamas (Figs. 11, 13, 14). Similarly, the large island of Hispaniola is a source of fauna for the small Turks and Caicos Islands in one Heterotermes termite lineage, Selenops spiders, and Centruroides scorpions (Figs. 4, 11, 13, 14).

Fig. 14
figure 14

Phylogeny and ancestral range estimation for Selenops spiders. BioGeoBEARS phylogram corresponding to the DIVALIKE + J B1a model (Table 1). The tree shows the single most probable geographical range at each node pre- and post-split. Colors correspond to Fig. 1. SA = South America; NA = North America; CA = Central America; GA = Greater Antilles; JA = Jamaica; SLA = Southern Lesser Antilles; NLA = Northern Lesser Antilles; FL = Florida; TCI = Turks and Caicos Islands; BA = Bahamas. Photo by S. Crews

Contrarily, there are instances of smaller islands as sources of taxa for larger islands, or even the mainland. Drosophila is unique in that there seems to have been a small species radiation in the Southern Lesser Antilles (Figs. 9, 10), and from this, two dispersal events occurred to the Greater Antilles (D. ornatipennis and D. dunni). Additionally, it appears that one lineage (D. acutilabella) emerged from a dispersal event to the Greater Antilles from Jamaica. Our analysis of Heterotermes indicates that one species (H. tenuis) from Guadeloupe island is basal to the South American lineage, with the caveats that sampling across the range of the species is poor, and these data are from a single gene (Fig. 4). In Platythyrea ants, there have been at least 2 dispersal events from the Greater Antilles to North America (Fig. 3).

Notably, the two butterfly species we examined show contrasting patterns of dispersal. Papilio are good dispersers and don’t seem to be affected by allopatric speciation or within-island vicariance, so whereas there are many species on the islands, ranges are broad and large species radiations do not appear to have occurred (Fig. 8). Calisto butterflies, however, are small and poor fliers. These butterflies show no evidence of dispersal from islands to the mainland and very few examples of dispersal from one island to the next; however, they are very diverse within islands likely due to being habitat specialists as well as poor dispersers (Figs. 6, 7) [27]. For Neotropical butterflies, Condamine et al. [76] emphasized the importance of islands as sources of diversity for the mainland. In Papilio butterflies there have been multiple dispersal events from the Greater Antilles to North America or to South America. This result differs from those of Lewis et al. [52] who found several “reverse” dispersal events, or events where the taxon is returning to the ancestral range of origin (dispersal from islands to the mainland); they found the entire clade of Heraclides to have a Caribbean origin whereas our analyses reveal a South American origin. Nevertheless, both analyses found that islands can be a source of diversity for other islands or even the mainland. Lewis et al. [52] attributed this in part to Papilio being big, strong fliers that favor edge habitats.

Founder event dispersal

The ability of BioGeoBEARS to incorporate “jump” dispersal into models has been advantageous in biogeographic studies (but see [58]). The program allows for models that emphasize the importance of founder event dispersal in island taxa [59]. In our analyses, models incorporating jump dispersal were favored in 8/12 analyses and were more common in taxa with older origins in the Caribbean (pre-35 my) (Table 2). However, in selenopid spiders, it was favored in the analysis with the younger tree, but not the older tree (see Additional file 1). The only analyses where a model incorporating jump dispersal was not favored were of ants and termites—organisms that only fly for a part of their adult life. Our results indicate that jump dispersal plays a role for both good dispersers, like Papilio, and poor dispersers, like spiders and scorpions. Founder event dispersal was always the favored model of dispersal for spiders and animals that fly during their entire adult life stage (Calisto, Papilio, and Drosophila). This finding corroborates results from Matos-Maraví et al. [27] for Calisto butterflies and Zhang et al. [74] for weevils.

Distance between areas of dispersal

To understand the influence of distance between neighboring landmasses (islands or continents) on dispersal probability we used manual dispersal multiplier (mdm) matrices in BioGeoBEARS analyses. Despite incorporating jump dispersal into our models, a model in which dispersal was very difficult (mdm values set to 0.01) was never favored (Table 2). We did find that in most taxa, distance-dependent models were favored (Platythyrea ants, Heterotermes termites, Papilio butterflies, centruroidine scorpions, Micrathena and Selenops spiders). For Calisto butterflies, distance did not affect dispersal when Euptychia was used as the sister taxon but did affect dispersal when other taxa were used. In Drosophila, two models were favored, and in one, distance affected dispersal, but in the other it did not.

Geologically, the Lesser Antilles islands can be divided into volcanic islands and two limestone volcanic arcs with the Southern Lesser Antilles older than the Northern Lesser Antilles (Additional file 1: Table S1). There were many episodes of gradual emergence, even within single islands, so a strict northern vs. southern boundary is elusive. These islands are younger and smaller than the Greater Antilles and have never been connected to the mainland, lowering overall dispersal probability. Previous biogeographic studies have recovered a clear pattern of a genetic split between organisms, congruent with Northern/Southern Lesser Antilles geology. This faunal-island split isn’t always between sister taxa, such as in Anolis lizards, Selenops spiders, and Centruroides scorpions where it appears there were two dispersal events: one from the Greater Antilles south to the Northern Lesser Antilles and one from South America north to the Southern Lesser Antilles, but without any discreet faunal break between islands. In this and a previous analysis [41], selenopid spider distribution agrees with the geologic data in that dispersal from South America to the older islands occurred sometime after 15 mya, whereas taxa didn’t disperse to the Northern Lesser Antilles until after 8 mya and are instead derived from the Greater Antillean fauna. The precise location of the split is also highly variable; for example, in Anolis lizards [77], it occurs between Martinique and Dominica, however it occurs further south between St. Vincent and St. Lucia [46] in lygaeid bugs [56], carabid beetles [40], Eleutherodactylus frogs [78], and bananaquits [79]. The location of the split recovered in Selenops spiders occurs between Dominica and Les Saintes, Guadeloupe. In the termite Heterotermes tenuis, we recovered some population structure between northern (Guadeloupe) and more southern islands (St. Lucia, St. Vincent, Grenada, Trinidad) [48]. For Dryas iulia butterflies [46] there have been multiple splits inferred, between both Martinique/Dominica and St. Vincent/St. Lucia. Multiple splits were also recovered in this study in Spintharus spiders: one between St. Kitts and Nevis, and another between Dominica and more southern islands.

Dispersal across the Isthmus of Panama/Central American Seaway

It should be noted that we are not aiming to disprove or provide evidence for a geologic hypothesis using biologic data, nor are we suggesting that our results indicate that the Central American Seaway had completely receded > 3 mya. The taxa we have examined do not require a land connection to be able to disperse and diversify across aquatic barriers. Most of our results indicate that the favored biogeographic model includes a single dispersal event between North and South America that occurred prior to 3 mya (Table 2). This result was also found for Hercules beetles [80], where it was suggested that the beetles could fly over water for short distances or that larvae could raft in downed wood. Taxa in which older dispersals (23 mya) were favored are centruroidine scorpions (Fig. 11), Micrathena spiders (Fig. 12), and Selenops spiders (Fig. 14). Selenops spiders and centruroidine scorpions are known from Chiapas and Dominican amber, indicating they have been present in southern North America and the Caribbean for at least 20 my.

The only taxon with a favored model of earliest dispersal at 15 my is Platythyrea ants (Table 2; Fig. 3). Seal et al. [53] mention that the ants probably disperse by walking over land, so perhaps they dispersed via a land connection or floated on flotsam. For some taxa, the models couldn’t be differentiated to determine a favored one. This was true in the analysis of Calisto butterflies, where 5 models were equiprobable (Table 2). In Nasutitermes and Papilio, both 15 and 8 my were favored (Table 2, Figs. 5, 8).

Several taxa were recovered with a favored model of 8 my or younger: Drosophila fruit flies, Spintharus spiders, and Heterotermes termites (Figs. 4, 9, 10, 13). Other studies have determined dispersal across the IoP/CAS to have occurred prior to 3 my, including those of butterflies [76] and Loxosceles spiders [81], although in the latter, the authors suggest dispersal occurred via a land bridge, which according to geologic data is not possible [18]. Bagley and Johnson [17] suggest that earlier dispersals were bidirectional and could have occurred via island hopping. Having better Central American sampling where possible, especially in taxa where we know there are species missing from the analyses, would certainly improve the level of uncertainty.

The role of GAARlandia in shaping Caribbean arthropod biodiversity

Several of the datasets we analyzed were younger than 32 mya, precluding GAARlandia as a transport path from South America to the Greater Antilles for these taxa. GAARlandia also did not play a role in the dispersal of Heterotermes termites, Selenops spiders, or Spintharus spiders (Figs. 4, 13, 14). In many cases our findings contrast those of previous research, although these studies lacked thorough sampling for parts of the distribution [25, 51] or hypothesis testing was not conducted [51]. For example, Dziki et al. [51] suggested that GAARlandia played a role in the dispersal of Spintharus from South America to the Greater Antilles. However, this was based on the coincidental timing of the diversification of a Caribbean clade rather than hypothesis testing. Our results, based on the relative AICc weights of the BioGeoBEARS analyses, favored a dispersal model without GAARlandia and also recovered a later date of dispersal into the Caribbean (~ 20 my). The results we obtained for Selenops spiders contradict those of Crews and Gillespie [25] in which similar hypothesis testing methods using Lagrange [82, 83] (which only uses the DEC model) found that models including GAARlandia were supported. Ancestral range estimation (Fig. 3) suggests that Selenops was in the Greater Antilles shortly after 40 my, before GAARlandia is hypothesized to have existed, and the model favored includes the jump dispersal parameter (Fig. 14, Table 2). This difference could be due to better dating methods now available, the use of different geologic dates/time slices, differences in the programs Lagrange and BioGeoBEARS, including the incorporation of jump dispersal, or the addition of more taxa/data, indicating how missing taxa affect biogeographic conclusions. For Micrathena spiders, the favored model indicated that GAARlandia played a role in dispersal to the Greater Antilles (Fig. 12). This is slightly confounding given the ancestral range estimation. Although there was dispersal to the Greater Antilles around 35–32 mya, these taxa came from North American ancestors. These findings contrast the results of McHugh et al. [26], and more sampling is required to explain these results. In centruroidine scorpions, two models had equal relative probabilities, one that included GAARlandia and another that did not (Fig. 11, Table 2). Esposito and Prendini [84], using Lagrange, previously found that the best fitting model was a distance-dependent dispersal model that included GAARlandia. According to the ancestral range estimation (Fig. 11), the timing of Heteroctenus scorpions entering the Caribbean is nearly simultaneous with the proposed timing of GAARlandia; however, the ancestral range estimation indicates a North American origin (although with some uncertainty), contradicting the results of Esposito and Prendini [84], who recovered a South American + Greater Antilles ancestral range using RASP [85, 86]. In Calisto butterflies, the results depend on whether Euptychia is the sister taxon or if other outgroups are used, emphasizing the importance of outgroup selection. For the former, distance-independent dispersal without GAARlandia was chosen as the favored model (Fig. 6). For the latter, distance-dependent dispersal with GAARlandia was the favored model (Fig. 7).

Species-area relationships and endemism

MacArthur and Wilson [70] suggested that island area may be responsible for lineage diversity or species numbers on islands. We know that the way in which island area is measured, habitat diversity, island age, and distance from the mainland are also important determiners [8, 87,88,89,90]. Generally, for all taxa examined, the Greater Antilles harbor more species than the Lesser Antilles, and the larger Greater Antilles islands, like Cuba and Hispaniola, typically have more species than the smaller islands of Puerto Rico and Jamaica. This expected pattern of Cuba > Hispaniola > Jamaica > Puerto Rico generally holds for the bug family Lyagaeidae [56], pseudoscorpions [91], scaritine beetles [92], Centruroides scorpions (Fig. 11), Micrathena spiders (Fig. 12), and Papilio butterflies (Fig. 8). However, in many taxa there are more species in Hispaniola than Cuba, such as in Platynus carabid beetles and Amphiacusta crickets [40, 49], the latter with 18 species in the Lesser Antilles. In Calisto butterflies, Hispaniola has nearly double the species of Cuba, and there are 2 species in the Bahamas but only 1 each in Puerto Rico and Jamaica. Grimaldi [63] suggested that there may be less endemism in fruit flies because they are good dispersers whose ranges usually cover large areas. This appears to be true in the cardini and dunni groups examined here, which have only a few species in the Greater Antilles and few single island endemics in the Lesser Antilles (Figs. 9, 10). For selenopid spiders, there are more endemics in Hispaniola than in Cuba (Fig. 12), although both islands harbor the same total number of species. Reasons for discrepancies from the typical pattern are worth a second look, and in the case of Cuba vs. Hispaniola could relate to differences in the toplogical complexity.

Lineage age and lineage radiations

Older lineages found on the Greater Antilles are generally ancestral to those of the Lesser Antilles; and Cuba, Hispaniola, or Puerto Rico lineages diverge earlier than Jamaican lineages. Species occurring in Cuba are the most basal for Selenops and the Caribbean endemic clade of Centruroides scorpions (Figs. 11, 14). Hispaniolan lineages are basal in Platythyrea ants and Calisto butterflies when Euptychia is used as a sister group; alternatively, Puerto Rican lineages are basal, which is also found in weevils [74] (Figs. 3, 6). Jamaican lineages are basal in Argiope spiders [93], however, it is possible that Argiope butchko is not a valid species but another lineage of Argiope argentata, in which case Cuban lineages would be basal. If we accept the proposal of Agnarsson et al. [93], then A. argentata has managed to disperse throughout the Caribbean to the exception of Cuba. Although competition could perhaps prevent A. argentata from becoming established in Cuba, multiple species of Argiope can be found in sympatry all over the world (e.g., [94]). Argiope butchko was described based solely on a few Cox1 differences and appears to show no morphological differences from A. argentata).

Lineage radiations have occurred in Cuba and Hispaniola in the butterfly species Calisto and Heteroctenus scorpion species, and all are single island endemics (Figs. 6, 7, 11). In species from the Caribbean clade of Centruroides scorpions, most exhibit patterns of within-island diversification or are single island endemics (Fig. 11). Considering that Calisto and scorpions are poor dispersers, the pattern of a single dispersal event followed by within-island radiation is expected. Conversely, Papilio, a good disperser, has not diversified on islands. Selenops spiders have undergone multiple lineage radiations in Hispaniola and Jamaica, and likely in Cuba, but without improved sampling it is unclear how many.

Island monophyly

Some groups exhibit patterns of dispersing to an island only once, whereas others have dispersed to a single island multiple times. Sampling is obviously important when examining island monophyly because missing species could cause false negatives as well as false positives. There does not appear to be island monophyly for the majority of taxa, including Papilio butterflies (Fig. 8), Drosophila flies (Figs. 9 and 10), Micrathena (Fig. 12) and Selenops spiders (Fig. 14), Platythyrea ants (Fig. 3), termites (Figs. 4 and 5), or weevils [74]. In the analysis of Calisto with Euptychia as the sister taxon, the butterflies have dispersed to Cuba at least 3 times from Hispaniola, and to Hispaniola at least twice, but to Jamaica and Puerto Rico only once (Fig. 6). However, in the analysis without Euptychia as the sister taxon, the butterflies dispersed to Cuba and Hispaniola only once each (Fig. 7). In Spintharus spiders there appears to be 1 Jamaican lineage, 1 Hispaniolan lineage, and 1 Puerto Rican lineage, but 2 Cuban lineages (Fig. 13). In Argiope argentata, lineages on Hispaniola and Jamaica are monophyletic [93]. Among scorpions, taxa are almost always monophyletic on islands, with the exception of 2 Cuban lineages of Centruroides.

Although Jamaica is considered one of the Greater Antilles, it has a rather different geologic history than Cuba, Hispaniola, and Puerto Rico. The island began emerging 40–30 mya and is thought to have had continuous available land since 15 mya, whereas the other Greater Antilles may have had continuous available land for much longer (~ 40 my). It is plausible that the Blue Mountains Block of eastern Jamaica has been above sea level for much longer and/or the part of western Jamaica that was attached to the Nicaraguan Rise of Central America was never submerged. Jamaica appears to have had a unique biogeographic history compared to the other Greater Antilles [95]. Of taxa included in this study, there are Jamaican endemic lineages in Calisto (1 species), Papilio (2 species), fruit flies (2 species), Centruroides scorpions (clade of 2 endemic species), Micrathena (1 species), Argiope (clade), Spintharus (clade), and Selenopid spiders (clade of 4 endemic species).


Terrestrial arthropods represent an immense wealth of biodiversity, dispersal abilities, life history traits, and ecologies. Their species assemblages in the Caribbean region underpin the importance of the Caribbean as both a major gateway for faunal exchange between the continents of the Western Hemisphere as well as a hotspot of endemic biodiversity and diversification. Our results indicate several biogeographic generalizations for terrestrial arthropods: South America is the predominant origin of the Caribbean arthropod fauna, and GAARlandia may have played an important role in aiding dispersal from South America to the Greater Antilles for some taxa but not universally. Founder (“jump” dispersal) events explain the majority of dispersal into the Caribbean, and most dispersal events occurred from the mainland to the islands and subsequent island hopping; this holds true for both good and poor dispersers. Distance is an important predictor of dispersal between landmasses, and there is evidence of islands as sources of diversity for the mainland in some taxa. Dispersal across what is now the Isthmus of Panama generally occurred prior to closure of the Central American Seaway 3 mya. The Greater Antilles harbor more species than the Lesser Antilles, and the larger Greater Antilles islands, Cuba and Hispaniola, typically have more species than the smaller islands of Puerto Rico and Jamaica. Among the Caribbean endemic taxa, basal taxa are most often distributed in the Greater Antilles, with the basal-most being Cuba, and derived taxa are most often distributed in the Lesser Antilles. Jamaican taxa are usually younger, endemic, and monophyletic for the island, whereas this is not the case for the other Greater Antillean islands.

The data presented herein provide a good baseline for understanding the biogeography of terrestrial arthropods in the Caribbean through explicit hypothesis testing. These methods allow for reproducible science and an improvement of our understanding of geologically complex regions. As genomic data continue to become more commonplace, genetic data from older museum specimens can be obtained, filling in gaps and helping resolve phylogenies and improve dating inferences. We hope that this paper will serve as a guide to help researchers determine what sorts of questions they can ask and answer with the data available to them.


Selection of datasets

We selected datasets using the following criteria: 1) Availability of a phylogenetic hypothesis. Because we are not attempting to represent ourselves as taxonomic experts for all terrestrial arthropods, we required a published phylogeny, corresponding datasets, and information that could be used to build time calibrated trees. Sometimes the datasets and calibration information were available in the publications’ supplementary materials, but sometimes only parts of datasets were available, or datasets did not correspond to the published trees, in which case we began the analyses from scratch. 2) Mostly complete sampling of species and distributions. As discussed in the introduction, the number of species of certain taxa and their distributions are often unknown for terrestrial arthropods. However, in some taxa, we know the true diversity, or at least that there are many unsampled lineages. We chose to omit these datasets from our analyses at this time and hope that they will be subjected to explicit hypothesis testing when more data become available. Datasets we were unable to use because of too little Caribbean or mainland sampling or too many missing taxa were Newportia centipedes [96], onychophora [97], Loxosceles spiders [81], and Amphiacusta crickets [49]. 3) Genetic differentiation sufficient for phylogenetic analysis. In the case of Argiope argentata spiders [93] and Dryas iulia and Heliconius charithonia butterflies [46], the taxa lacked enough genetic differentiation to get at the questions we are asking.

Using these criteria, we were left with 10 datasets, some better than others for addressing the hypotheses we wished to test. Our analyses consist of Platythyrea ants, Heterotermes and Nasutitermes termites, butterflies from the genera Calisto and Papilio, Drosophila flies, centruroidine scorpions, and Micrathena, Spintharus, and Selenops spiders. These taxa provide a good representation of sedentary and more vagile taxa with various life histories. When possible, molecular data were supplemented with new data from GenBank (Additional file 1: Tables S5–S15), including both markers and species not available in previous publications. In some cases, the addition of taxa for better calibration precluded the use of certain markers – our strategy was to maximize the sampling effort, including geographic, taxonomic, and phylogenetic coverage. New molecular sequence data were collected only in the case of Selenops spiders.

Phylogenetic and dating analysis

For each dataset (Platythyrea ants, Heterotermes termites, Nasutitermes termites, Calisto butterflies, Papilio butterflies, Drosophila fruit flies, centruroidine scorpions, Micrathena spiders, Spintharus spiders, and Selenops spiders) a phylogeny was estimated in MrBayes [98] (details of each analysis can be found in Additional file 1). Missing data were treated as “?”, and the BIC was used to select the best partitioning scheme and models in PartitionFinder2 [99] and jModelTest2 [100, 101]. Additionally, all data blocks were separated by gene (for non-coding genes) and by codon (for coding genes) for the input files. Further information about each dataset (e.g., genes used, partitioning schemes, models of molecular evolution, generations) can be found Additional file 1. Analyses were run with the default setting of 2 analyses with 4 chains, sampling every 1000 generations. Clades with support values ≥0.95 are considered well-supported. Tracer v. 1.6 [102] was used to examine ESS values to ensure they were > 200 and to examine stationarity plots. All trees were viewed in FigTree v1.4.3 [103]. Each phylogeny was compared to the published phylogeny, then used as a topological constraint for the dating analysis in BEAST 2 [104]. If our results were inconsistent with those previously published or certain relationships were unresolvable, multiple possibilities were tested.

BEAST 2 input files were constructed using BEAUti. Constrained MrBayes molecular phylogenies were dated in BEAST 2 using either fossil calibrations or a molecular clock rate. For information on the dating parameters, please refer to the detailed methods for each taxon in Additional file 1. Burn-in for Tracer 1.6 analyses was 10%, and for trees in TreeAnnotator, it was set to 25% unless otherwise noted. All MrBayes and BEAST 2 analyses were run using XSEDE [105] on the CIPRES Science Gateway [106]. The .xml files and output trees used for subsequent BioGeoBEARS analyses are available on figshare.

Biogeographic model testing

We conducted biogeographic hypothesis testing and ancestral range reconstruction using the BEAST 2 time-calibrated phylogenies in BioGeoBEARS [107]. The models used (Table 1) to test the biogeographic hypotheses generally follow those of Crews and Gillespie [25] and consider GAARlandia, the various proposed closure dates of the IoP/CAS, and various multipliers of “dispersal ability” in shaping biogeographic histories. The most basic model considered the probability of dispersal without GAARlandia vs. a model that considers GAARlandia. Manual dispersal multiplier matrices were employed to further model the effect of distance and dispersal ability [42] as follows:

  1. a)

    Distance dependent dispersal. The probability of reaching a land area decreased as distance between land areas increased. If there was a land connection available during a particular time period, the manual dispersal multiplier (mdm) value was increased (Note: mdm values were only set to 1 in very closely adjacent areas).

  2. b)

    Distance independent dispersal. Distance did not play a role in dispersal and establishment probability (presence absence of landmasses only). If there was a hypothesized or known land connection, mdm values were set to 1 (land) and with no connection, set to 0.5. During the proposed timing of GAARlandia, the probability of dispersal from South America to the Greater Antilles increased from 0.5 to 1.

  3. c)

    Vicariance only model. A very low probability of dispersal. Land connections remained set to 1, but overwater dispersal was lowered to 0.01. The mdm values from South America to the Greater Antilles during the proposed timing of GAARlandia were increased from 0.01 to 1. When there was no land available for a designated landmass, probabilities were set to 0.0000001 [67].

Table 1 Biogeographic models analyzed in BioGeoBEARS. Each dispersal or vicariance scenario was tested using the six models available in BioGeoBEARS (DEC, DEC+J, DIVALIKE, DIVALIKE+J, BAYAREALIKE, BAYAREALIKE+J)

Finally, we considered various geologic times of taxa moving across the IoP. (See Table 1 and Fig. 2 for details and citations.) All data were analyzed using all models available in the BioGeoBEARS package (DEC, DEC + J, DIVALIKE, DIVALIKE+J, BAYAREALIKE, BAYAREALIKE+J) [107,108,109]. This resulted in a maximum of 252 models per dataset. The best model was evaluated using the AICc weights [59].

Model limitations

The Caribbean islands have moved and increased in area and connectedness through time. We acknowledge that these processes are gradual, and geologic dates may not be exact. We’ve therefore elected to examine dispersal vs. vicariance using large generalized areas for both geologically-grouped islands and continents (e.g., “Northern Lesser Antilles” vs. each individual island, “Greater Antilles” vs. each individual island) in spite of potential criticisms of oversimplification (e.g., [110] of [25]). We’ve opted to employ a strategy that avoids over-parameterization while utilizing the data available. We encourage other authors to use the generalized biogeoraphic models we have provided as input files on figshare as a baseline – using different dispersal probabilities, different time periods, different areas, etc. – to test their own hypotheses.

For the time periods and geographic areas used for each analysis, see Tables 1 and 2. We relied on authors and experts for taxa distribution information. While it is possible that taxa may have broader or even narrower distributions, this potentially would not have a large effect on the results because we used broader land areas as our geographic areas for analyses. Geologic dates based on the current published literature for the first appearance of geographic areas guided our manual dispersal multiplier values and are provided in Table 1 and Additional file 1: Table S1.

Table 2 Major biogeographic conclusions and summary of analyses for terrestrial arthropod taxa in the Caribbean biodiversity hotpot. (SA = South America, NA = North America, CA = Central America, GA = Greater Antilles, SLA = Southern Lesser Antilles, NLA = Northern Lesser Antilles, BA = Bahamas and/or Turks and Caicos Islands, FL = Florida, CU = Cuba, HI = Hispaniola, JA=Jamaica, PR = Puerto Rico, BAR=Barbados)

Availability of data and materials

All data generated or analyzed during this study are included in the additional files, are available on GenBank, or are available as a downloadable package at



Bahamas and/or Turks and Caicos Islands




Central America


Central American Seaway






Greater Antilles


Greater Antilles Aves Ridge




Isthmus of Panama




Million years




North America


Northern Lesser Antilles


Puerto Rico


South America


Southern Lesser Antilles


  1. 1.

    Poey F. Memorias sobre la historia natural de la isla de Cuba, acompañadas de sumarios latinos y extractos en frances, por Felipe Poey Habana. Cuba: Impr De Barcina; 1861. p. 1851–8.

    Google Scholar 

  2. 2.

    Liebherr JK. Zoogeography of Caribbean insects. 1st ed. Ithaca: Cornell University Press; 1988.

  3. 3.

    Alayón G. La familia Selenopidae (Arachnida: Araneae) en Cuba. Solenodon. 2005;5:10–52.

    Google Scholar 

  4. 4.

    Scheffrahn RH, Křeček J, Chase JA, Maharajh B, Mangold JR. Taxonomy, biogeography, and notes on termites (Isoptera: Kalotermitidae, Rhinotermitidae, Termitidae) of the Bahamas and Turks and Caicos Islands. Ann Entomol Soc Amer. 2006;99:464–86.

    Article  Google Scholar 

  5. 5.

    Crews SC. A revision of the spider genus Selenops Latreille, 1819 (Arachnida, Araneae, Selenopidae) in North America, Central America and the Caribbean. ZooKeys. 2011:1–82.

    Article  Google Scholar 

  6. 6.

    Bloom T, Binford G, Alayón G, Esposito L, Peterson I, Nishida A, et al. Discovery of two new species of eyeless spiders within a single Hispaniola cave. J Arachnol. 2014:148–54.

    Article  Google Scholar 

  7. 7.

    Crews SC, Esposito LA, Cala-Riquelme F. The arachnids of Hellshire Hills, Jamaica. Caribbean Nat. 2015;28:1–14.

    Google Scholar 

  8. 8.

    Gillespie RG, Roderick GK. Arthropods on islands: colonization, speciation, and conservation. Annu Rev Entomol. 2002:595–632.

    CAS  PubMed  Article  Google Scholar 

  9. 9.

    Donnelly TW. In: Liebherr JK, editor. Geologic constraints on Caribbean biogeography. Ithaca: Comstock Publishing Associates, Cornell University Press; 1988. p. 15–37.

    Google Scholar 

  10. 10.

    Dengo G, Case JE. The geology of North America, vol. H, the Caribbean region (a decade of north American geology) Boulder. Colorado: Geologic Society of America; 1990.

    Google Scholar 

  11. 11.

    Iturralde-Vinent MA, MacPhee RDE. Paleogeography of the Caribbean region: implications for Cenozoic biogeography. Bull Amer Mus Nat Hist. 1999:238;1–95.

  12. 12.

    Iturralde-Vinent MA. Meso-Cenozoic Caribbean paleogeography: implications for the historical biogeography of the region. Int Geol Rev. 2006:791–827.

    Article  Google Scholar 

  13. 13.

    James KH. In situ origin of the Caribbean: discussion of data. Geol Soc, London, Special Publications. 2009:77–125.

    Article  Google Scholar 

  14. 14.

    Hedges SB, Hass CA, Maxson LR. Caribbean biogeography: molecular evidence for dispersal in West Indian terrestrial vertebrates. P Natl Acad Sci USA. 1992:1909–13.

    CAS  Article  Google Scholar 

  15. 15.

    Crother BI, Guyer C. Caribbean historical biogeography: was the dispersal vicariance debate eliminated by an extraterrestrial bolide? Herpetologica. 1996:52;440–65.

  16. 16.

    MacPhee RDE, Iturralde-Vinent MA. The interpretation of Caribbean paleogeography: reply to Hedges. Proc Int Symposium on Insular Vert Evol: the paleontological approach; 2005.

    Google Scholar 

  17. 17.

    Bagley JC, Johnson JB. Phylogeography and biogeography of the lower central American Neotropics: diversification between two continents and between two seas. Biol Rev. 2014:767–90.

    PubMed  Article  Google Scholar 

  18. 18.

    O'Dea A, Lessios HA, Coates AG, Eytan RI, Restrepo-Moreno SA, Cione AL, et al. Formation of the Isthmus of Panama. Sci Adv. 2016;2:e1600883.

    Article  PubMed  PubMed Central  Google Scholar 

  19. 19.

    Montes C, Cardona A, McFadden R, Morón SE, Silva CA, Restrepo-Moreno S, et al. Evidence for middle Eocene and younger land emergence in Central Panama: implications for isthmus closure. Geol Soc Amer Bull. 2012;124:780–99.

    CAS  Article  Google Scholar 

  20. 20.

    Coates AG, Stallard RF. How old is the isthmus of Panama? B Mar Sci. 2013;89:801–13.

    Article  Google Scholar 

  21. 21.

    Osborne AH, Newkirk DR, Groeneveld J, Martin EE, Tiedemann R, Frank M. The seawater neodymium and lead isotope record of the final stages of central American seaway closure. Paleoceanography. 2014;29:715–29.

    Article  Google Scholar 

  22. 22.

    Tavares VC, Simmons NB. Phylogenetic relationships and biogeography of Short-faced stenodermatine bats: preliminary results and hypotheses. Bat Res News. 2000;41:143.

    Google Scholar 

  23. 23.

    Dávalos LM. Phylogeny and biogeography of Caribbean mammals. Biol J Linn Soc. 2004:373–94.

    Article  Google Scholar 

  24. 24.

    Alonso R, Crawford AJ, Bermingham E. Molecular phylogeny of an endemic radiation of Cuban toads (Bufonidae: Peltophryne) based on mitochondrial and nuclear genes. J Biogeogr. 2011:434–51.

    Article  Google Scholar 

  25. 25.

    Crews SC, Gillespie RG. Molecular systematics of Selenops spiders (Araneae: Selenopidae) from north and Central America: implications for Caribbean biogeography. Biol J Linn Soc. 2010:288–322.

    Article  Google Scholar 

  26. 26.

    McHugh A, Yablonsky C, Binford G, Agnarsson I. Molecular phylogenetics of Caribbean Micrathena (Araneae: Araneidae) suggests multiple colonisation events and single island endemism. Invertebr Syst. 2014:337–49.

    Article  Google Scholar 

  27. 27.

    Matos-Maraví P, Águila RN, Peña C, Miller JY, Sourakov A, Wahlberg N. Causes of endemic radiation in the Caribbean: evidence from the historical biogeography and diversification of the butterfly genus Calisto (Nymphalidae: Satyrinae: Satyrini). BMC Evol Biol. 2014:199.

  28. 28.

    Říčan O, Piálek L, Zardoya R, Doadrio I, Zrzavý J. Biogeography of the Mesoamerican Cichlidae (Teleostei: Heroini): colonization through the GAARlandia land bridge and early diversification. J Biogeogr. 2013:579–93.

    Article  Google Scholar 

  29. 29.

    Fritsch PW. Multiple geographic origins of Antillean Styrax. Syst Bot. 2003:28;421–30.

  30. 30.

    van Ee B, Berry PE, Riina R, Gutiérrez Amaro JE. Molecular phylogenetics and biogeography of the Caribbean-centered Croton subgenus Moacroton (Euphorbiaceae s.s.). Bot Rev. 2008;74:132–65.

    Article  Google Scholar 

  31. 31.

    Ali JR. Colonizing the Caribbean: is the GAARlandia land-bridge hypothesis gaining a foothold? J Biogeogr. 2011;39:431–3.

    Article  Google Scholar 

  32. 32.

    O'Dea A, Hoyos N, Rodríguez F, Degarcia B, De Gracia C. History of upwelling in the tropical eastern Pacific and the paleogeography of the isthmus of Panama. Palaeogeogr Palaeocl. 2012;348:59–66.

    Article  Google Scholar 

  33. 33.

    Leigh EG, O’Dea A, Vermeij GJ. Historical biogeography of the isthmus of Panama. Biol Rev. 2014;89:148–72.

    Article  PubMed  Google Scholar 

  34. 34.

    Bacon CD, Mora A, Wagner WL, Jaramillo CA. Testing geologic models of evolution of the isthmus of Panama in a phylogenetic framework. Bot J Linn Soc. 2013:287–300.

    Article  Google Scholar 

  35. 35.

    Bacon CD, Silvestro D, Jaramillo C, Smith BT, Chakrabarty P, Antonelli A. Biologic evidence supports an early and complex emergence of the isthmus of Panama. P Natl Acad Sci USA. 2015;112:6110–5.

    CAS  Article  Google Scholar 

  36. 36.

    Barbour T. A contribution to the zoogeography of the West Indies, with special reference to amphibians and reptiles. Mem Mus Comp Zool Harvard College. 1914:44;205–359.

  37. 37.

    Darlington PJ. The origin of the fauna of the greater Antilles, with discussion of dispersal of animals over water and through the air. Q Rev Biol. 1938;13:274–300.

    Article  Google Scholar 

  38. 38.

    Ricklefs RE, Lovette IJ. The roles of island area per se and habitat diversity in the species-area relationships of four lesser Antillean faunal groups. J Anim Ecol. 1999;68:1142–60.

    Article  Google Scholar 

  39. 39.

    Penney D, Wheater CP, Selden PA. Resistance of spiders to cretaceous–tertiary extinction events. Evolution. 2003;57:2599–607.

    Article  PubMed  Google Scholar 

  40. 40.

    Liebherr JK. In: Liebherr JK, editor. Biogeographic patterns of West Indian Platynus carabid beetles (Coleoptera). Ithaca: Comstock Publishing Associates, Cornell University Press; 1988. p. 121–52.

    Google Scholar 

  41. 41.

    Crews SC, Puente-Rolón AR, Rutstein E, Gillespie RG. A comparison of populations of island and adjacent mainland species of Caribbean Selenops (Araneae: Selenopidae) spiders. Mol Phylogenet Evol. 2010;54:970–83.

    Article  PubMed  Google Scholar 

  42. 42.

    Van Dam MH, Matzke NJ. Evaluating the influence of connectivity and distance on biogeographic patterns in the south-western deserts of North America. J Biogeogr. 2016;43:1514–32.

    Article  Google Scholar 

  43. 43.

    Crews SC, Hedin MC. Studies of morphological and molecular phylogenetic divergence in spiders (Araneae: Homalonychus) from the American southwest, including divergence along the Baja California peninsula. Mol Phylogenet Evol. 2006:470–87.

    CAS  PubMed  Article  Google Scholar 

  44. 44.

    Leaché AD, Crews SC, Hickerson MJ. Two waves of diversification in mammals and reptiles of Baja California revealed by hierarchical Bayesian analysis. Biol Lett. 2007:646–50.

    PubMed  PubMed Central  Article  Google Scholar 

  45. 45.

    Hollocher H, Williamson M. Island hopping in Drosophila: patterns and processes [and discussion]. Phil Trans R Soc London B: Biol Sci. 1996;351:735–43.

    CAS  Article  Google Scholar 

  46. 46.

    Davies N, Bermingham E. The historical biogeography of two Caribbean butterflies (Lepidoptera: Heliconiidae) as inferred from genetic variation at multiple loci. Evolution. 2002;56:573–89.

    CAS  Article  PubMed  Google Scholar 

  47. 47.

    Wilder JA, Hollocher H. Recent radiation of endemic Caribbean Drosophila of the dunni subgroup inferred from multilocus DNA sequence variation. Evolution. 2003;57:2566–79.

    CAS  Article  PubMed  Google Scholar 

  48. 48.

    Szalanski AL, Scheffrahn RH, Austin JW, Křeček J, Su N-Y. Molecular phylogeny and biogeography of Heterotermes (Isoptera: Rhinotermitidae) in the West Indies. Ann Entomol Soc Amer. 2004;97:556–66.[0556:MPABOH]2.0.CO;2.

    CAS  Article  Google Scholar 

  49. 49.

    Oneal E, Otte D, Knowles LL. Testing for biogeographic mechanisms promoting divergence in Caribbean crickets (genus Amphiacusta). J Biogeogr. 2010;37:530–40.

    Article  Google Scholar 

  50. 50.

    Sourakov A, Zakharov EV. “Darwin’s butterflies”? DNA barcoding and the radiation of the endemic Caribbean butterfly genus Calisto (Lepidoptera, Nymphalidae, Satyrinae). Comp Cytogenet. 2011;5:191–210.

    Article  PubMed  PubMed Central  Google Scholar 

  51. 51.

    Dziki A, Binford GJ, Coddington JA, Agnarsson I. Spintharus flavidus in the Caribbean—a 30 million year biogeographical history and radiation of a ‘widespread species’. PeerJ. 2015;3:e1422.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  52. 52.

    Lewis DS, Sperling FA, Nakahara S, Cotton AM, Kawahara AY, Condamine FL. Role of Caribbean Islands in the diversification and biogeography of Neotropical Heraclides swallowtails. Cladistics. 2015;31:291–314.

    Article  Google Scholar 

  53. 53.

    Seal JN, Kellner K, Trindl A, Heinze J. Phylogeography of the parthenogenic ant Platythyrea punctata: highly successful colonization of the West Indies by a poor disperser. J Biogeogr. 2011;38:868–82.

    Article  Google Scholar 

  54. 54.

    Esposito LA, Bloom T, Caicedo-Quiroga L, Alicea-Serrano AM, Sánchz-Ruíz JA, May-Collado LJ, et al. Islands within islands: diversification of tailless whip spiders (Amblypygi, Phrynus) in Caribbean caves. Mol Phylogenet Evol. 2015;93:107–17.

    Article  PubMed  Google Scholar 

  55. 55.

    Rosen DE. Geologic hierarchies and biogeographic congruence in the Caribbean. Ann Mo Bot Gard. 1985:636–59.

    Article  Google Scholar 

  56. 56.

    Slater J. In: Liebherr JK, editor. Zoogeography of West Indian Lygaeidae (Hemiptera). Ithaca: Comstock Publishing Associates, Cornell University Press; 1988. p. 38–60.

    Google Scholar 

  57. 57.

    Riddle BR, Hafner DJ, Alexander LF, Jaeger JR. Cryptic vicariance in the historical assembly of a Baja California Peninsular Desert biota. P Natl Acad Sci USA. 2000:14438–43.

    CAS  Article  Google Scholar 

  58. 58.

    Ree R, 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(4):741–9.

    Article  Google Scholar 

  59. 59.

    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.

    Article  PubMed  Google Scholar 

  60. 60.

    Lieberman BS. Phylogenetic biogeography with and without the fossil record: gauging the effects of extinction and paleontological incompleteness. Palaeogeogr Palaeoclimatol Palaeoecol. 2002;178:39–52.

    Article  Google Scholar 

  61. 61.

    Morrone JJ. Biogeographic areas and transition zones of Latin America and the Caribbean islands based on panbiogeographic and cladistic analyses of the entomofauna. Annu Rev Entomol. 2006;51:467–94.

    CAS  Article  PubMed  Google Scholar 

  62. 62.

    Echeverry A, Morrone JJ. Generalized tracks, area cladograms and tectonics in the Caribbean. J Biogeogr. 2013;40:1619–37.

    Article  Google Scholar 

  63. 63.

    Grimaldi DA. In: Liebherr JK, editor. Relicts in the Drosophilidae (Diptera). Ithaca: Comstock Publishing Associates, Cornell University Press; 1988. p. 183–213.

    Google Scholar 

  64. 64.

    Browne DJ, Peck SB, Ivie MA. The longhorn beetles (Coleoptera Cerambycidae) of the Bahama Islands with an analysis of species-area relationships, distribution patterns, origin of the fauna and an annotated species list. Trop Zool. 1993;6:27–53.

    Article  Google Scholar 

  65. 65.

    Wilson EO. In: Liebherr JK, editor. The biogeography of the West Indian ants (Hymenoptera: Formicidae). Ithaca: Comstock Publishing Associates, Cornell University Press; 1988. p. 214–30.

    Google Scholar 

  66. 66.

    Ramos JA. In: Liebherr JK, editor. Zoogeography of the auchenorrhynchous Homoptera of the Greater Antilles. Ithaca: Comstock Publishing Associates, Cornell University Press; 1988. p. 61–70.

    Google Scholar 

  67. 67.

    Matzke NJ.; Accessed 31 Aug 2017.

  68. 68.

    Scheffrahn RH, Krecek J, Szalanski AL, Austin JW, Roisin Y. Synonymy of two arboreal termites (Isoptera: Termitidae: Nasutitermitinae): Nasutitermes corniger from the Neotropics and N. polygynus from New Guinea. Fla Entomol. 2005:28–33.;2.

    Article  Google Scholar 

  69. 69.

    Magalhães ILF, Santos AJ. Phylogenetic analysis of Micrathena and Chaetacis spiders (Araneae: Araneidae) reveals multiple origins of extreme sexual size dimorphism and long abdominal spines. Zool J Linn Soc-Lond. 2012:14–53.

  70. 70.

    MacArthur RH, Wilson EO. An equilibrium theory of insular zoogeography. Evolution. 1963:373–87.

    Article  Google Scholar 

  71. 71.

    MacArthur RH, Wilson EO. The theory of island biogeography. 1st ed. Princeton: Princeton University Press; 1967.

  72. 72.

    Nicholson KE, Glor RE, Kolbe JJ, Larson A, Hedges SB, Losos JB. Mainland colonization by island lizards. J Biogeogr. 2005;32:929–38.

    Article  Google Scholar 

  73. 73.

    Glor RE, Losos JB, Larson A. Out of Cuba: overwater dispersal and speciation among lizards in the Anolis carolinensis subgroup. Mol Ecol. 2005:2419–32.

    CAS  PubMed  Article  Google Scholar 

  74. 74.

    Zhang G, Basharat U, Matzke N, Franz NM. Model selection in statistical historical biogeography of Neotropical insects—the Exophthalmus genus complex (Curculionidae: Entiminae). Mol Phylogenet Evol. 2016:226–39.

    PubMed  Article  Google Scholar 

  75. 75.

    Sturge RJ, Jacobsen F, Rosensteel BB, Neale RJ, Omland KE. Colonization of South America from Caribbean Islands confirmed by molecular phylogeny with increased taxon sampling. Condor. 2009;111:575–9.

    Article  Google Scholar 

  76. 76.

    Condamine FL, Silva-Brandao KL, Kergoat GJ, Sperling FA. Biogeographic and diversification patterns of Neotropical Troidini butterflies (Papilionidae) support a museum model of diversity dynamics for Amazonia. BMC Evol Biol. 2012;12:82.

    Article  PubMed  PubMed Central  Google Scholar 

  77. 77.

    Losos JB, Thorpe RS.Evolutionary Diversification of Caribbean Anolis Lizards. In: Adaptive Speciation, eds. Dieckmann U, Doebeli M, Metz JAJ & Tautz D. International Institute for Applied Systems Analysis. Cambridge: Cambridge University Press.; 2004:322–344.

  78. 78.

    Kaiser H, Hardy JD Jr, Green DM. Taxonomic status of Caribbean and south American frogs currently ascribed to Eleutherodactylus urichi (Anura: Leptodactylidae). Copeia. 1994.

    Article  Google Scholar 

  79. 79.

    Seutin G, Klein NK, Ricklefs RE, Bermingham E. Historical biogeography of the bananaquit (Coereba flaveola) in the Caribbean region: a mitochondrial DNA assessment. Evolution. 1994:1041–61.

    CAS  PubMed  Google Scholar 

  80. 80.

    Huang J-P. The great American biotic interchange and diversification history in Dynastes beetles (Scarabaeidae; Dynastinae). Zool J Linnean Soc. 2016:88–96.

    Article  Google Scholar 

  81. 81.

    Binford GJ, Callahan MS, Bodner MR, Rynerson M, Núñez PB, Ellison CE, et al. Phylogenetic relationships of Loxosceles and Sicarius spiders are consistent with Western Gondwanan vicariance. Mol Phylogenet Evol. 2008:538–53.

    PubMed  Article  Google Scholar 

  82. 82.

    Ree RH, Moore BR, Webb CO, Donoghue MJ. A likelihood framework for inferring the evolution of geographic range on phylogenetic trees. Evolution. 2005:2299–311.

    PubMed  Article  Google Scholar 

  83. 83.

    Ree RH, Smith SA. Maximum likelihood inference of geographic range evolution by dispersal, local extinction, and cladogenesis. Syst Biol. 2008:4–14.

    PubMed  Article  Google Scholar 

  84. 84.

    Esposito LA, Prendini L. Island ancestors and New World biogeography: a Case study from the scorpions (Buthidae: Rhopalurusinae). Sci Rep. 2019:1–11.

  85. 85.

    Yu Y, Harris AJ, He XJ. RASP (reconstruct ancestral state in phylogenies): a tool for historical biogeography. Mol Phylogenet Evol. 2010:848–50.

    PubMed  Article  Google Scholar 

  86. 86.

    Yu Y, Harris AJ, Blair C, He XJ. S-DIVA (statistical dispersal-Vicariance analysis): a tool for inferring biogeographic histories. Mol Phylogenet Evol. 2015:46–9.

    PubMed  Article  Google Scholar 

  87. 87.

    Losos JB. Phylogenetic perspectives on community ecology. Ecology. 1996:1344–54.

    Article  Google Scholar 

  88. 88.

    Ricklefs RE, Bermingham E. The concept of the taxon cycle in biogeography. Global Ecol and Biogeogr. 2002:353–61.

    Article  Google Scholar 

  89. 89.

    Ricklefs R, Bermingham E. The West Indies as a laboratory of biogeography and evolution. Rev Zool Bot Afr. 2008:2393–413.

    Article  Google Scholar 

  90. 90.

    Triantis KA, Nogués-Bravo D, Hortal J, Borges PA, Andersen H, Fernández-Palacios J. Measurements of area and the (island) species–area relationship: new directions for an old pattern. Oikos. 2008:1555–9.

    Article  Google Scholar 

  91. 91.

    Cosgrove JG, Agnarsson I, Harvey MS, Binford GJ. Pseudoscorpion diversity and distribution in the West Indies: sequence data confirm single island endemism for some clades, but not others. J Arachnol. 2016:257–71.

    Article  Google Scholar 

  92. 92.

    Cenzi De Ré F, Loreto EL, Robe LJ. Gene and species trees reveal mitochondrial and nuclear discordance in the Drosophila cardini group (Diptera: Drosophilidae). Invertebr Biol. 2010; 353–367. doi:

    Article  Google Scholar 

  93. 93.

    Agnarsson I, LeQuier SM, Kuntner M, Cheng R-C, Coddington JA, Binford G. Phylogeography of a good Caribbean disperser: Argiope argentata (Araneae, Araneidae) and a new ‘cryptic’ species from Cuba. ZooKeys. 2016:25–44.

    Article  Google Scholar 

  94. 94.

    Levi HW. The spider genera Gea and Argiope in America (Araneae: Araneidae). Bull Mus Comp Zool. 1968:136;319–52.

  95. 95.

    Buskirk RE. Zoogeographic patterns and tectonic history of Jamaica and the northern Caribbean. J Biogeogr. 1985.

    Article  Google Scholar 

  96. 96.

    Edgecombe GD, Vahtera V, Giribet G, Kaunisto P. Species limits and phylogeography of Newportia (Scolopendromorpha) and implications for widespread morphospecies. ZooKeys. 2015:65–77.

    Article  Google Scholar 

  97. 97.

    Murienne J, Daniels SR, Buckley TR, Mayer G, Giribet G. A living fossil tale of Pangaean biogeography. P Roy Soc Lond B Bio. 2014;281:20132648.

    Article  Google Scholar 

  98. 98.

    Ronquist F, Huelsenbeck JP. MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003;19:1572–4.

    CAS  Article  Google Scholar 

  99. 99.

    Lanfear R, Frandsen PB, Wright AM, Senfeld T, Calcott B. PartitionFinder 2: new methods for selecting partitioned models of evolution for molecular and morphological phylogenetic analyses. Mol Biol Evol. 2016;34:772–3.

    CAS  Article  Google Scholar 

  100. 100.

    Guindon S, Gascuel O. A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003;52:696–794.

    Article  PubMed  Google Scholar 

  101. 101.

    Darriba D, Taboada GL, Doallo R, Posada D. jModelTest 2: more models, new heuristics and parallel computing. Nat Methods. 2012;9:772.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  102. 102.

    Rambaut A, Suchard MA, Xie D, Drummond AJ. Tracer v1.6. 2014;

  103. 103.

    Rambaut A. FigTree. Accessed 2017.

  104. 104.

    Boukaert R, Heled J, Kühnert D, Vaughan T, Wu C-H, Suchard MA, et al. BEAST 2: a software platform for Bayesian evolutionary analysis. PLoS Comput Biol. 2014;10:e1003537.

    CAS  Article  Google Scholar 

  105. 105.

    Towns J, Cockerill T, Dahan M, Foster I, Gaither K, Grimshaw A, et al. XSEDE: accelerating scientific discovery. Comput Sci Eng. 2014;16:62–74.

    CAS  Article  Google Scholar 

  106. 106.

    Miller MA, Pfeiffer W, Schwartz T. Creating the CIPRES Science Gateway for inference of large phylogenetic trees, Proc of the Gateway Computing Env Workshop (GCE); 2010. p. 1–8.

    Google Scholar 

  107. 107.

    Matzke NJ. BioGeoBEARS: BioGeography with Bayesian (and Likelihood) Evolutionary Analysis in R Scripts. R package, version 0.2.1, published July 27, 2013 at: 2013.

    Google Scholar 

  108. 108.

    Matzke NJ. rexpokit: R wrappers for EXPOKIT. R package, version 0.24.2, published July 15, 2013 at: 2013.

    Google Scholar 

  109. 109.

    Matzke NJ, Sidje RB. cladoRcpp: C++ implementations of phylogenetic calculations. R package, version 0.14.2, published July 15, 2013 at: 2013.

    Google Scholar 

  110. 110.

    Santamaria CA, Mateos M, Hurtado LA. Diversification at the narrow sea-land interface in the Caribbean: phylogeography of endemic supralittoral Ligia isopods. Front Ecol and Evol. 2014;2:42.

    Article  Google Scholar 

  111. 111.

    Nichols SW. In: Liebherr JK, editor. Kaleidoscopic biogeography of West Indian Scaratinae (Coleoptera: Carabidae). Ithaca: Comstock Publishing Associates, Cornell University Press; 1988. p. 71–121.

    Google Scholar 

  112. 112.

    Matzke NJ. Probabilistic historical biogeography: new models for founder-event speciation, imperfect detection, and fossils allow improved accuracy and model-testing. Front Biogeogr. 2013:8;242–8.

  113. 113.

    Maddison WP, Maddison DR. Mesquite: a modular system for evolutionary analysis. Version 3.2. 2017.

    Google Scholar 

  114. 114.

    Zhang J, Kapli P, Pavlidis P, Stamatakis A. A general species delimitation method with applications to phylogenetic placements. Bioinformatics. 2013:2869–76.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  115. 115.

    Schmidt C. Molecular phylogenetics of ponerine ants (Hymenoptera: Formicidae: Ponerinae). Zootaxa. 2013:201–50.

    PubMed  Article  Google Scholar 

  116. 116.

    Brady SG, Larkin L, Danfort BN. In: Hedges SB, Kumar S, editors. Bees, ants, and stinging wasps (Aculeata). Oxford: Oxford University Press; 2009.

  117. 117.

    Katoh K, Misawa K, Kuma K, Miyata T. MAFFT: a novel method for rapid multiple sequence alignment based on fast Fourier transform. Nucleic Acids Res. 2002:3059–66.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  118. 118.

    Ware JL, Grimaldi DA, Engel MS. The effects of fossil placement and calibration on divergence times and rates: an example from the termites (Insecta: Isoptera). Arthropod Struct Dev. 2010:204–19.

    Article  Google Scholar 

  119. 119.

    Vaidya G, Lohman DJ, Meier R. SequenceMatrix: concatenation software for the fast assembly of multi-gene datasets with character set and codon information. Cladistics. 2011:171–80.

    Article  Google Scholar 

  120. 120.

    Roy V, Constantino R, Chassany V, Giusti-Miller S, Diouf M, Mora P, et al. Species delimitation and phylogeny in the genus Nasutitermes (Termitidae: Nasutitermitinae) in French Guiana. Mol Ecol. 2014:902–20.

    CAS  PubMed  Article  Google Scholar 

  121. 121.

    Sukumaran J, Holder MT. DendroPy: a Python library for phylogenetic computing. Bioinformatics. 2010:1569–71.

    CAS  PubMed  Article  Google Scholar 

  122. 122.

    Sukumaran J, Holder MT. SumTrees: Phylogenetic Tree Summarization. 2015;

    Google Scholar 

  123. 123.

    Peña C, Nylin S, Wahlberg N. The radiation of Satyrini butterflies (Nymphalidae: Satyrinae): a challenge for phylogenetic methods. Zool J Linnean Soc. 2011:64–87.

    Article  Google Scholar 

  124. 124.

    Lotts K, Naberhaus T. Butterflies and moths of North America. 2017.

    Google Scholar 

  125. 125.

    Obbard DJ, Maclennan J, Kim KW, Rambaut A, O'Grady PM, Jiggins FM. Estimating divergence dates and substitution rates in the Drosophila phylogeny. Mol Biol Evol. 2012:3459–73.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  126. 126.

    Brisson JA, Wilder J, Hollocher H. Phylogenetic analysis of the cardini group of Drosophila with respect to changes in pigmentation. Evolution. 2006:1228–41.

    CAS  PubMed  Article  Google Scholar 

  127. 127.

    Viera J. Molecular evolution group: Drosophila distributions. 2016.

    Google Scholar 

  128. 128.

    Drosophila Species Stock Center Geography. 2017. Accessed 2017.

  129. 129.

    Ashburner M, Carson HL, Thompson JN. The Genetics and Biology of Drosophila, vol. Volume 3a. London: Academic Press, London; 1981.

  130. 130.

    Dimitrov D, Lopardo L, Girbet G, Arnedo MA, Álvarez-Padilla F, Hormiga G. Tangled in a sparse spider web: single origin of orb weavers and their spinning work unravelled by denser taxonomic sampling. P Roy Soc Lond B Bio. 2011.

    PubMed  Article  Google Scholar 

  131. 131.

    Wheeler WC, Coddington JA, Crowley LM, Dimitrov D, Goloboff P, Griswold C, et al. The spider tree of life: phylogeny of Araneae based on target-gene analyses from an extensive taxon sampling. Cladistics. 2016:1–43.

    Article  Google Scholar 

  132. 132.

    Kallal RJ, Hormiga G. Systematics, phylogeny and biogeography of the Australasian leaf-curling orb-weaving spiders (Araneae: Araneidae: Zygiellinae), with a comparative analysis of retreat evolution. Zool J Linn Soc-Lond In press.

  133. 133.

    Gregorič M, Agnarsson I, Blackledge TA, Kuntner M. Phylogenetic position and composition of Zygiellinae and Caerostris, with new insight into orb-web evolution and gigantism. Zool J Linn Soc-Lond. 2015:225–43.

    Article  Google Scholar 

  134. 134.

    Stucky BJ. SeqTrace: a graphical tool for rapidly processing DNA sequencing chromatograms. J Biomol Tech. 2012;23:90.

    Article  PubMed  PubMed Central  Google Scholar 

  135. 135.

    Stamatakis A. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics. 2014:1312–3.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  136. 136.

    Li W, Godzik A. Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics. 2006:1658–9.

    CAS  PubMed  Article  Google Scholar 

  137. 137.

    Huang Y, Niu B, Gao Y, Fu L, Li W. CD-HIT suite: a web server for clustering and comparing biologic sequences. Bioinformatics. 2010:680–2.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  138. 138.

    Fu L, Niu B, Zhu Z, Wu S, Li W. CD-HIT: accelerated for clustering the next-generation sequencing data. Bioinformatics. 2012:3150–2.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  139. 139.

    Polotow D, Carmichael A, Griswold CE. Total evidence analysis of the phylogenetic relationships of Lycosoidea spiders (Araneae, Entelegynae). Invertebr Syst. 2015:124–63.

  140. 140.

    Ramírez MJ. The morphology and phylogeny of dionychan spiders (Araneae: Araneomorphae). Bull Amer Mus Nat Hist. 2014:1–374.

    Article  Google Scholar 

  141. 141.

    Bidegaray-Batista L, Arnedo MA. Gone with the plate: the opening of the Western Mediterranean basin drove the diversification of ground-dweller spiders. BMC Evol Biol. 2011;317.

    PubMed  PubMed Central  Article  Google Scholar 

  142. 142.

    Selden PA, Wang Y. Fossil spiders (Araneae) from the Eocene Green River formation of Colorado. Arthropoda Selecta. 2014:23;207–19.

  143. 143.

    Fouke BW. Quaternary geology and depositional history of Provodenciales, Turks and Caicos Islands, British West Indies. M.Sc. Thesis; 1984.

    Google Scholar 

  144. 144.

    Cuezzo C, Carrijo TF, Cancello EM. Transfer of two species from Nasutitermes Dudley to Cortaritermes Mathews (Isoptera: Termitidae: Nasutitermitinae). Austral Entomol. 2015:172–9.

    Article  Google Scholar 

  145. 145.

    Futch DG. Hybridization tests within the cardini species group of the genus Drosophila. Univ Texas Publ. 1962;6205:539–54.

    Google Scholar 

  146. 146.

    Heed WB. Genetic characteristics of island populations. Univ Texas Publ. 1962;6205:502–18.

    Google Scholar 

  147. 147.

    Garrison NL, Rodriguez J, Agnarsson I, Coddington JA, Griswold CE, Hamilton CA, et al. Spider phylogenomics: untangling the spider tree of life. PeerJ. 2016;4:e1719.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  148. 148.

    Degnan JH, Salter LA. Gene tree distributions under the coalescent process. Evolution. 2005:24–37.

    PubMed  Article  Google Scholar 

  149. 149.

    Kubatko LS, Degnan JH. Inconsistency of phylogenetic estimates from concatenated data under coalescence. Syst Biol. 2007:17–24.

    CAS  PubMed  Article  Google Scholar 

  150. 150.

    Iturralde-Vinent MA. Aspectos geológicos de la biogeografía de Cuba. Ciencias de la Tierra y el Espacio. 1982:5;85–100.

  151. 151.

    Bouysse P, Westercamp D. Subduction of Atlantic aseismic ridges and late Cenozoic evolution of the Lesser Antilles island arc. Tectonophysics. 1990:349357–55380.

    Article  Google Scholar 

  152. 152.

    Maury RC, Westbrook GK, Baker PE, Bouysse Ph, Westercamp D. Geology of the Lesser Antilles. In: Dengo G and Case JE, editors. Boulder. Geologic Society of America. 1990.

    Google Scholar 

  153. 153.

    Stephan J-F, Mercier-de-Lepinay B, Calais E, Tardy M, Beck C, Carfantan J-C, et al. Paleogeodynamic maps of the Caribbean: 14 steps from Lias to present. Bull Soc Géol Fr. 1990;6:6–919.

    Google Scholar 

  154. 154.

    Mann P, Schubert C, Burke K. In: Dengo G, Case JE, editors. Review of Caribbean neotectonics. Boulder: Geologic Society of America; 1990.

    Google Scholar 

  155. 155.

    Randazzo AF, Jones DS. The geology of Florida; 1997.

    Google Scholar 

  156. 156.

    Hine AC. Geologic history of Florida: major events that formed the sunshine state; 2013.

    Google Scholar 

  157. 157.

    Speed RC, Larue DK. Barbados: Architecture and implications for accretion. J Geophys Res-Sol Ea. 1982;87(B5):3633–43.

    Article  Google Scholar 

  158. 158.

    Fouke BW. Quaternary geology and depositional history of Providenciales, Turks and Caicos Islands, British West Indies. M.Sc. Thesis; 1984.

    Google Scholar 

  159. 159.

    Wanless HR, Dravis JJ. Carbonate environments and sequences of Caicos platform. 1989; 28th International Geologic Congress, Fieldtrip Guidebook T-374.

    Google Scholar 

  160. 160.

    Morgan WA, Harris PM. Developing models and analogs for isolated carbonate platforms - Holocene and Pleistocene carbonates of Caicos Platform, British West Indies. 2008; SEPM Core Workshop 22.

    Book  Google Scholar 

  161. 161.

    Smith AM, Fisher BL, Hebert PD. DNA barcoding for effective biodiversity assessment of a hyperdiverse arthropod group: the ants of Madagascar. Phil Trans R Soc London B: Biol Sci. 2005;360:1825–34.

    CAS  Article  Google Scholar 

  162. 162.

    Moreau CS, Bell CD, Vila R, Archibald SB, Pierce NE. Phylogeny of the ants: diversification in the age of angiosperms. Science. 2006:101–4.

    CAS  PubMed  Article  Google Scholar 

  163. 163.

    Spagna JC, Vakis AI, Schmidt CA, Patek SN, Zhang X, Tsutsui N, et al. Phylogeny, scaling, and the generation of extreme forces in trap-jaw ants. J Exp Biol. 2008:2358–68.

    PubMed  Article  Google Scholar 

  164. 164.

    Marino A, Sabatini AG, Carpana E, Lodesani M, Mantovani B. ND2 and CO1 mitochondrial genes in Apis mellifera L.: A Molecular Approach to Mediterranean Population Monitoring. Proc 6th Euro Bee Conf. 2002:4;157.

  165. 165.

    Field J, Ohl M, Kennedy M. A molecular phylogeny for digger wasps in the tribe Ammophilini (Hymenoptera, Apoidea, Sphecidae). Syst Entomol. 2011;36:732–40.

    Article  Google Scholar 

  166. 166.

    Heraty J, Ronquist F, Carpenter JM, Hawks D, Schulmeister S, Dowling AP, et al. Evolution of the hymenopteran megaradiation. Mol Phylogenet Evol. 2011:73–88.

    PubMed  Article  Google Scholar 

  167. 167.

    Eaton TD, Jones SC, Jenkins TM. Species diversity of Puerto Rican Heterotermes (Dictyoptera: Rhinotermitidae) revealed by phylogenetic analyses of two mitochondrial genes. J Insect Sci. 2016;111.

  168. 168.

    Saldarriaga JF, Gile GH, James ER, Horak A, Scheffrahn RH, Keeling PJ. Morphology and molecular phylogeny of Pseudotrichonympha hertwigi and Pseudotrichonympha Paulistana (Trichonymphea, Parabasalia) from neotropical rhinotermitids. J Eukaryot Microbiol. 2011:487–96.

    PubMed  Article  Google Scholar 

  169. 169.

    Lo N, Kitade O, Miura T, Constantino R, Matsumoto T. Molecular phylogeny of the Rhinotermitidae. Insect Soc. 2004:365–71.

    Article  Google Scholar 

  170. 170.

    Casalla R, Scheffrahn R, Korb J. Cryptotermes colombianus a new drywood termite and distribution record of Cryptotermes in Colombia. ZooKeys. 2016:39–52.

    Article  Google Scholar 

  171. 171.

    Buijang NS, Harrison NA, Su NY. A phylogenetic study of endo-beta-1, 4-glucanase in higher termites. Insect Soc. 2014:29–40.

    Article  Google Scholar 

  172. 172.

    Dejernaes M, Klass KD, Picker MD, Damgaard J. Phylogeny of cockroaches (Insecta, Dictyoptera, Blattodea), with placement of aberrant taxa and exploration of out-group sampling. Syst Entomol. 2012:65–83.

    Article  Google Scholar 

  173. 173.

    Legendre F, Whiting MF, Bordereau C, Cancello EM, Evans TA, Grandcolas P. The phylogeny of termites (Dictyoptera: Isoptera) based on mitochondrial and nuclear markers: implications for the evolution of the worker and pseudergate castes, and foraging behaviors. Mol Phylogenet Evol. 2008:615–27.

    CAS  PubMed  Article  Google Scholar 

  174. 174.

    Reid NM, Addison SL, West MA, Lloyd-Jones G. The bacterial microbiota of Stolotermes ruficeps (Stolotermitidae), a phylogenetically basal termite endemic to New Zealand. FEMS Microbiol Ecol. 2014:678–88.

    CAS  PubMed  Article  Google Scholar 

  175. 175.

    Kambhampati S. A phylogeny of cockroaches and related insects based on DNA sequence of mitochondrial ribosomal RNA genes. Proc Natl Acad Sci. 1995:2017–20.

    CAS  Article  Google Scholar 

  176. 176.

    Liu C. Molecular phylogeny of an enigmatic Heterochaeta Wood-Mason, 1889 (Mantodea: Toxoderidae: Oxyothespinae), with description of one new species from China. unpubl. data 2015; GenBank.

  177. 177.

    Atkinson L, Adams ES, Crozier RH. Microsatellite markers for the polygamous termite Nasutitermes corniger (Isoptera: Termitidae). Mol Ecol Resour. 2007:299–301.

    Article  CAS  Google Scholar 

  178. 178.

    Miura T, Roisin Y, Matsumoto T. Molecular phylogeny and biogeography of the nasute termite genus Nasutitermes (Isoptera: Termitidae) in the Pacific tropics. Mol Phylogenet Evol. 2000:1–10.

    CAS  PubMed  Article  Google Scholar 

  179. 179.

    Akoth P. Genetic characterization and distribution of termites in Taita Taveta County, Kenya. Ph.D. dissertation, College of Pure Science, Jomo Kenyatta University of Agriculture and Technology. Available from: Accessed 2017.

  180. 180.

    Svenson GJ, Whiting MF. Reconstructing the origins of praying mantises (Dictyoptera, Mantodea): the roles of Gondwanan vicariance and morphological convergence. Cladistics. 2009;25:468–514.

    Article  Google Scholar 

  181. 181.

    Pellens R, D'Haese CA, Bellés X, Piulachs MD, Legendre F, Wheeler WC, et al. The evolutionary transition from subsocial to eusocial behaviour in Dictyoptera: phylogenetic evidence for modification of the “shift-in-dependent-care” hypothesis with a new subsocial cockroach. Mol Phylogenet Evol. 2007;43:616–26.

    CAS  Article  PubMed  Google Scholar 

  182. 182.

    Inward D, Beccaloni G, Eggleton P. Death of an order: a comprehensive molecular phylogenetic study confirms that termites are eusocial cockroaches. Biol Lett. 2007:331–5.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  183. 183.

    Fukuda H, Koshikawa S, Kusaka Y, Toda M, Katoh T. Molecular phylogeny of the immigrans-tripuncata lineage (Diptera: Drosophilidae). unpubl. data 2014; GenBank.

  184. 184.

    Hollocher H, Hatcher JL, Dyreson EG. Evolution of abdominal pigmentation differences across species in the Drosophila dunni subgroup. Evolution. 2000:2046–56.

    CAS  PubMed  Article  Google Scholar 

  185. 185.

    Colón-Parilla WV, Pérez-Chiesa I. Partial characterization and evolution of Adh-Adhr in Drosophila dunni. Biochem Genet. 2007:225–38.

    PubMed  Article  CAS  Google Scholar 

  186. 186.

    Robe LJ, Loreto EL, Valente VL. Radiation of the, Drosophila subgenus (Drosophilidae, Diptera) in the Neotropics. J Zool Syst Evol Res. 2010:310–21.

    Article  Google Scholar 

  187. 187.

    O'Grady P, DeSalle R. Out of Hawaii: the origin and biogeography of the genus Scaptomyza (Diptera: Drosophilidae). Biol Lett. 2008:195–9.

    PubMed  PubMed Central  Article  Google Scholar 

  188. 188.

    Hatadani LM, McInerney JO, Fonseca de Medeiros H, ACM J, de Azeredo-Espin AM, Klaczko LB. Molecular phylogeny of the Drosophila tripunctata and closely related species groups (Diptera: Drosophilidae). Mol Phylogenet Evol. 2009:595–600.

    CAS  PubMed  Article  Google Scholar 

  189. 189.

    Chaumot A, Da Lage JL, Maestro O, Martin D, Iwema T, Brunet F, et al. Molecular adaptation and resilience of the insect’s nuclear receptor USP. BMC Evol Biol. 2012;12:199.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  190. 190.

    Dyer KA, Jaenike J. Evolutionarily stable infection by a male-killing endosymbiont in Drosophila innubila. Genetics. 2004:1443–55.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  191. 191.

    Remsen J, O'Grady P. Phylogeny of Drosophilinae (Diptera: Drosophilidae), with comments on combined analysis and character support. Mol Phylogenet Evol. 2002:249–64.

    PubMed  Article  Google Scholar 

  192. 192.

    Katoh TK, Komoda M, Watada M. A survey of the drosophilid fauna and genetic diversity of two dominant species of Drosophila in Ehime prefecture, Japan. Biogeogr. 2016;18:35–46.

    Google Scholar 

  193. 193.

    Albalat R, Gonzàlez-Duarte R. Adh and Adh-dup sequences of Drosophila lebanonensis and D. immigrans: interspecies comparisons. Gene. 1993;126:171–8.

    CAS  PubMed  Article  Google Scholar 

  194. 194.

    Wolff JN, Camus MF, Clancy DJ, Dowling DK. Complete mitochondrial genome sequences of thirteen globally sourced strains of fruit fly (Drosophila melanogaster) form a powerful model for mitochondrial research. Mitochondr DNA Part A. 2016:4672–2674.

    Article  CAS  Google Scholar 

  195. 195.

    Kellerman V, van Heerwaarden B, Sgrò CM, Hoffmann AA. Fundamental evolutionary limits in ecological traits drive Drosophila species distributions. Science. 2009:1244–6.

    CAS  PubMed  Article  Google Scholar 

  196. 196.

    Li F-B. The complete mitochondrial genome of Drosophila melanogaster unpubl. data 2014; GenBank.

  197. 197.

    Yeh S-D, von Grotthuss M, Gandasetiawan KA, Jayasekera S, Xia XQ, Chan C, et al. Functional divergence of the miRNA transcriptome at the onset of Drosophila metamorphosis. Mol Biol Evol. 2014:2557–72.

    CAS  PubMed  Article  Google Scholar 

  198. 198.

    Gramates LS, Marygold SJ, dos Santos G, Urbano J-M, Antonazzo G, Matthews BB, et al. FlyBase at 25: looking to the future. Nucleic Acids Res. 2017:D663–71.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  199. 199.

    Yotoko KS, Medeiros HF, Solferini VN, Klaczko LB. A molecular study of the systematics of the Drosophila tripunctata group and the tripunctata radiation. Mol Phylogenet Evol. 2003:614–9.

    CAS  PubMed  Article  Google Scholar 

  200. 200.

    Jaenike J, Dyer KA, Cornish C, Minhas MS. Asymmetrical reinforcement and Wolbachia infection in Drosophila. PLoS Biol. 2006:e325.

    PubMed  PubMed Central  Article  Google Scholar 

  201. 201.

    Satta Y, Takahata N. Evolution of Drosophila mitochondrial DNA and the history of the melanogaster subgroup. Proc Natl Acad Sci. 1990:9558–62.

    CAS  Article  Google Scholar 

  202. 202.

    Meiklejohn CD, Holmbeck MA, Siddiq MA, Abt DN, Rand DM, Montooth KL. An incompatibility between a mitochondrial tRNA and its nuclear-encoded tRNA synthetase compromises development and fitness in Drosophila. PLoS Genet. 2013;9:e1003238.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  203. 203.

    Junqueira ACM, Azeredo-Espin AML, Paulo DF, Marinho MAT, Tomsho LP, Drautz-Moses DI, et al. Large-scale mitogenomics enables insights into Schizophora (Diptera) radiation and population diversity. Sci Rep-UK. 2016;6:21762.

    CAS  Article  Google Scholar 

  204. 204.

    Clary DO, Wahleithner JA, Wolstenholme DR. Transfer RNA genes in Drosophila mitochondrial DNA: related 5’flanking sequences and comparisons to mammalian mitochondrial tRNA genes. Nucleic Acids Res. 1983;11:2411–25.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  205. 205.

    Llopart A, Herrig D, Brud E, Stecklein Z. Sequential adaptive introgression of the mitochondrial genome in Drosophila yakuba and Drosophila santomea. Mol Ecol. 2014:1124–36.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  206. 206.

    Cheng R-C, Kuntner M. Phylogeny suggests nondirectional and isometric evolution of sexual size dimorphism in argiopine spiders. Evolution. 2014:2861–72.

    PubMed  Article  Google Scholar 

  207. 207.

    Agnarsson I, Blackledge TA. Can a spider web be too sticky? Tensile mechanics constrains the evolution of capture spiral stickiness in orb-weaving spiders. J Zool. 2009:134–40.

    Article  Google Scholar 

  208. 208.

    Álvarez-Padilla F, Dimitrov D, Giribet G, Hormiga G. Phylogenetic relationships of the spider family Tetragnathidae (Araneae, Araneoidea) based on morphological and DNA sequence data. Cladistics. 2009:109–46.

    Article  Google Scholar 

  209. 209.

    Blackledge TA, Scharff N, Coddington JA, Szuts T, Wenzel JW, Hayashi CY, Agnarsson I. Reconstructing web evolution and spider diversification in the molecular era. Proc Natl Acad Sci. 2009:106;5229–34.

  210. 210.

    Coddington JA, Agnarsson I, Cheng R-C, Čandek K, Driskell A, Frick H, et al. DNA barcode data accurately identify higher taxa. PeerJ. 2016.

    PubMed  PubMed Central  Article  Google Scholar 

  211. 211.

    Arnedo MA, Coddington J, Agnarsson I, Gillespie RG. From a comb to a tree: phylogenetic relationships of the comb-footed spiders (Araneae, Theridiidae) inferred from nuclear and mitochondrial genes. Mol Phylogenet Evol. 2004:225–45.

    CAS  PubMed  Article  Google Scholar 

  212. 212.

    Simon C, Frati F, Beckenbach A, Crespi B, Liu H, Flook P. Evolution, weighting, and phylogenetic utility of mitochondrial gene sequences and a compilation of conserved polymerase chain reaction primers. Ann Entomol Soc Amer. 1994:651–701.

    CAS  Article  Google Scholar 

  213. 213.

    Colgan DJ, McLauchlan A, Wilson GDF, Livingston SP, Edgecombe GD, Macaranas J, et al. Histone H3 and U2 snRNA DNA sequences and arthropod molecular evolution. Aust J Zool. 1998:419–37.

    Article  Google Scholar 

Download references


We would like to acknowledge Peter Roopnarine for helpful suggestions that improved the manuscript. We thank Matt Van Dam for answering an endless stream of questions about BioGeoBEARS analyses, and Delano Lewis for providing us with BEAST input files. We are grateful to Bob Kallal for his advice regarding araneid outgroups and BEAST calibration points. We acknowledge Corrie Moreau, Barbara Thorne, Rudolf Scheffrahn, and Michele Esposito for answering questions about specimen localities and distributions. We thank Boni Cruz and Anna Sellas for help with lab work and acknowledge the CarBio project and National Museum of Natural History for providing selenopid extractions. We are grateful to Paul Mann, Bruce Fouke, and René Maury for answering lots of questions about geology. Finally, we are grateful to Judy Gallagher, Mark Yokoyama, Eduardo Lopez, Bernard Dupont and Rudolf Scheffrahn for providing images and allowing us to alter them for some of our figures.



Author information




SCC and LAE contributed to the conception and design of the study; SCC acquired the data and conducted analyses; SCC and LAE contributed to interpreting the data and drafting and revising the manuscript. All authors have read and approved the manuscript.

Corresponding author

Correspondence to Lauren A. Esposito.

Ethics declarations

Ethics approval and consent to participate


Consent for publication


Competing interests

The authors declare that they have no competing interests.

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: Table S1.

Geologic dates used in BioGeoBEARS analyses. Table S2. Partitions and models used for each dataset. Table S3. Summary of untested models for datasets from older taxa. Table S4. Summary of untested models for datasets from younger taxa. Figure S1. Ancestral range estimation for Platythyrea ants. Table S5. Sequences for Platythyrea analyses. Table S6. Bayes Factors from Platythyrea analysis. Figure S2. Phylogeny and ancestral range estimation for Heterotermes termites. Table S7. Sequences for Heterotermes analyses. Figure S3. Phylogeny and ancestral range estimation for Nasutitermes termites. Table S8. Sequences for Nasutitermes analyses. Figure S4. Phylogeny and ancestral range estimation for Calisto butterflies (Euptychia as the outgroup). Figure S5. Ancestral range estimation for Calisto butterflies (Euptychia is not the sister taxon). Figure S6. Map of the geographical areas used in the Calisto BioGeoBEARS analyses. Table S9. Calisto butterfly sample information. Figure S7. Ancestral range estimation for Papilio butterflies. Table S10. Papilio butterfly sample information Figure S8. Ancestral range estimation for Drosophila flies. Table S11. Sequences for Drosophila analyses. Figure S9. Ancestral range estimation for centruroidine scorpions. Table S12. Sequences for centruroidine scorpion analyses. Figure S10. Phylogeny and ancestral range estimation for Micrathena spiders. Table S13. Sequences for Micrathena analyses. Figure S11. Ancestral range estimation for Spintharus spiders. Table S14. Sequences for Spintharus analyses. Figure S12. Ancestral range estimation for Selenops spiders using more recent dated phylogenies. Figure S13. Ancestral range estimation for Selenops spiders using more recent dated phylogenies. Figure S14. Selenopidae tree from the Bayesian analysis using no redundant haplotypes. Figure S15. Selenopidae tree from the RAxML analysis using no redundant haplotypes. Table S15. Sequences for Selenops analyses. Table S16. Primers and PCR protocols used for amplification of DNA from Selenops species [111,112,113,114,115,116,117,118,119,120,121,121,122,123,124,125,126,127,128,129,130,131,132,133,134,135,136,137,138,139,140,141,142,143,144,145,146,147,148,149,150,151,152,153,154,155,156,157,158,159,160,161,162,163,164,165,167,168,169,170,171,172,173,174,175,176,177,178,179,180,181,182,183,184,185,186,187,188,189,190,191,192,193,194,195,196,197,198,199,200,201,202,203,204,205,206,207,208,209,210,211,212,213,213].

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Crews, S.C., Esposito, L.A. Towards a synthesis of the Caribbean biogeography of terrestrial arthropods. BMC Evol Biol 20, 12 (2020).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Biogeography
  • Insects
  • Arachnids
  • BioGeoBEARS
  • Isthmus of Panama
  • Central American seaway
  • Islands
  • Probabilistic inference