Published November 13, 2021 | Version v1
Dataset Open

Around the world in 10 million years: rapid dispersal of a kleptoparasitoid spider wasp (Pompilidae: Ceropales)

  • 1. Commonwealth Scientific and Industrial Research Organisation
  • 2. University of Göttingen
  • 3. Universidade Vila Velha
  • 4. Utah State University

Description

Aim: Our aim was to estimate the historical biogeography of the kleptoparasitoid genus Ceropales and to determine the processes leading to its current worldwide distribution. We tested hypotheses of dispersal and vicariance scenarios underlying its widespread distribution.

Location: Worldwide.

Methods: Data from two nuclear markers (the D2–D3 regions of the 28S ribosomal RNA and long-wavelength rhodopsin) and one mitochondrial marker (cytochrome c oxidase I) for 52 specimens of Ceropales were used to reconstruct a dated phylogeny based on Bayesian inference. Two calibration points were used from previous analyses including all pompilids under a lognormal relaxed molecular clock to estimate lineage divergence times. We compared the fit of 12 biogeographical models, modifying the base BioGeoBEARS models to include a dispersal adjacency matrix. Base BioGeoBEARS models allow different cladogenetic processes: DEC (subset sympatry, narrow vicariance), DIVALIKE (narrow and wide vicariance), BAYAREALIKE (widespread sympatry), and +J versions of these that allow jump dispersal. Using the model with the best AIC score, we performed Biogeographic Stochastic Mapping (BSM) in order to infer biogeographic processes. We simulated 200 BSM using the DEC+J model and the consensus tree for the BEAST analysis.

Results: The origin of crown-group Ceropales was in the early Miocene, ca. 10.6 Ma (15.7–6.5 95% HPD), and eight dispersal events explain its widespread distribution. A constrained model, where only adjacent areas were allowed for dispersal had the highest likelihood under DEC+J model.

Main Conclusions: The widespread distribution of Ceropales can be explained by eleven jump-dispersal events that took place in a period of ca. 10 million years. Two separate dispersals at different times happened from the Eurasia to the Nearctic. These probably occurred across the Bering land bridge in the late Miocene and Pliocene. Dispersal from North and Mesoamerica to South America took place four independent times from the late Miocene to close to present time. Dispersal to the Ethiopian region from Eurasia occurred in the late Miocene and Pliocene. Dispersal back to Eurasia from the Ethiopian region took place three times independently in the Pliocene to close to present time. Dispersal to the Australian region took place from the late Miocene to the Pleistocene.

Notes

BAYAREALIKE+J_200BSM_stochastic_maps.pdf: The 200 simulations of possible biogeographic processes performed through Biogeographic Stochastic Mapping (BSM) in BioGeoBEARS. 

Cero28SCOIOpsNOintrPOFeb28names.nex: A concatenated supermatrix of all molecular markers used in this study (28S, COI and Opsin). A MrBayes block defines character blocks and models of molecular evolution used. 

CeroLessInformative_2.xml: BEAST .xml file generated using BEAUti. Specifies the data matrix, partitions, models, priors and mcmc chain settings.

Cero-FINAL.tre: Ceropales chronogram produced from the BEAST analysis (see methods). This is also input for BioGeoBEARS script.

CeropalesAreasConstrained_byContinent_Apr2020.R: BioGeoBEARS R script modified from Nicholas J. Matzke, License: GPL-3 http://cran.r-project.org/web/licenses/GPL-3, used for model estimation and generating AIC statistics for each model. This should run using the .tre file provided and the files"multipliers_new_continent.txt" which includes dispersal multipliers used for the analysis and "distributionCeropales_LessAreas_continents.txt" which includes taxa distributions.

Ceropales-BSM-JBI-revision-continents.R: BioGeoBEARS R script used for Biogeographic Stochastic Mapping and summarising events per node.

Ceropales-node-numbers.pdf: Ceropales phylogeny plotted with node numbers (these are referred to in the file "Event-simulation-results.xlsx")

Event-simulation-results.xlsx: Results from BSM summarised per node and statistics for each of them.

Files

BAYAREALIKE+J_200BSM_stochastic_maps.pdf

Files (2.1 MB)

Name Size Download all
md5:33762e7c188606ebd32990c47e05fa68
1.6 MB Preview Download
md5:067e116540f7d3e8beab3cce3fc60788
2.8 kB Download
md5:bc5c1336d0508f1f761fec1656888b22
113.8 kB Download
md5:ae35bf39963fb46a8b097e73d5c69311
162.3 kB Preview Download
md5:6e83d08e336482e7722055db3689d8d5
10.4 kB Download
md5:b0c92fd01c8edfb7d67cc407527fb06b
12.4 kB Preview Download
md5:97357a1e4b485cd44fb5395292406f82
46.8 kB Download
md5:ffb9ad537d9d1453345a6e4ba586ce5e
1.0 kB Preview Download
md5:d74cee5ba97c9d2c814e8de762e4ef5c
186.7 kB Download
md5:0a503d02311f60994a6606e10f1343c2
72 Bytes Preview Download