Data and code from: Diversification dynamics in the global radiation of gobies
Data files
Jul 03, 2026 version files 33.42 MB
-
ASTRAL.zip
6.40 MB
-
BAMM.zip
381.22 KB
-
BEAST.zip
1.11 MB
-
CRABS.zip
2.33 MB
-
Gobiiformes_script_R.R
28.94 KB
-
IQtree.zip
7.88 MB
-
README.md
7.67 KB
-
SSE.zip
11.60 KB
-
TACT.zip
15.05 MB
-
TESS.zip
234.50 KB
Abstract
We examine the dynamics and drivers of diversification in gobies, an exceptionally species-rich group with a cosmopolitan distribution in aquatic habitats that exhibits extensive ecological variation coupled with low phenotypic diversity. We establish their evolutionary timescale using a phylogenomic hypothesis of relationships calibrated with multiple fossils and identify a pronounced increase in diversification rate in the clade containing Gobiidae and Oxudercidae. This shift is correlated with both phenotypic and ecological traits; elevated diversification is found among marine lineages and those on reefs as well as in those possessing a fused ventral pelvic disc, a key innovation that facilitates stability on benthic substrates. Our analyses additionally reveal diversification fluctuations at the Eocene-Oligocene transition and in the mid-Miocene, periods of extinction and recovery in global marine habitats. Achieving even near-complete taxon sampling is unlikely in hyperdiverse groups such as gobies, but we show that evolutionary dynamics may still be reliably inferred with less than complete taxon sampling by validating our diversification shift patterns using completely sampled, stochastically resolved trees and confirming the robustness of our rate estimates across congruence classes.
Dataset DOI: 10.5061/dryad.c59zw3rq8
Description of the data and file structure
This dataset provides the input files, trees, matrices, graphics of results, and R code for the paper "Diversification dynamics in the global radiation of gobies" by C. E. Thacker, W. T. McCraney, and M. E. Alfaro.
README.md is this file, describing the Dryad file structure and contents
Gobiiformes_script_R.R is the R code used to run the TESS CoMET, CRABS, BiSSE, MuSSE, and HISSE analyses
ASTRAL.zip folder:
- GobiiformesUCE270genetrees.tre is the concatenated set of trees for individual UCE loci
- GobiiformesUCE270Astraltree.tre is the ASTRAL tree inferred from the set of UCE locus trees
- GobiiformesUCE270Astraltree.nex is the ASTRAL tree in nexus format
- UCElocusalignments folder contains the alignments of individual UCE loci and trees for each inferred by IQtree. For each locus, the file uce-##.iqtree is the results summary, uce-##.treefile is the maximum likelihood tree with bootstrap branch supports, and uce-##.contree is the 50% majority rule consensus tree of the bootstrap replicates.
BAMM.zip folder:
- Gobiiformes75p270calibrated.nwk is the fossil-calibrated BEAST tree for Gobiiformes in newick format
- Gobiiformes75p229calibrated.nwk is the fossil-calibrated BEAST tree for Gobioidei in newick format
- bammGobiiformes75p229.txt is the control file for BAMM
- sample_probs.txt is the file of lineage-specific sampling fractions
- Gobiiformes75p229BAMMbestshift.pdf shows the most favored shift configuration inferred by BAMM
- Gobiiformes75p229BAMMshiftset.pdf shows the set of shift configurations inferred by BAMM
BEAST.zip folder:
These files are for calibrating the Gobiiformes phylogeny using three nonoverlapping subsets of 10 UCE loci:
-
Gobiiformes270UCE10.xml, Gobiiformes270UCE10II.xml, and Gobiiformes270UCE10III.xml are the .xml input files for the BEAST calibration analyses.
-
Gobiiformes270UCE10combined.con is the fossil-calibrated BEAST tree combined from the three replicate runs. It is the same as Gobiiformes75p270calibrated.nwk used in the BAMM and SSE analyses.
CRABS.zip folder:
BAMM subfolder contains output plots from the CRABS analysis of BAMM output:
- Gobioidei229crabsbammmodels.pdf: Diversification rates in the reference model
- Gobioidei229crabsbammextinct.pdf: Diversification rate functions over a range of constant extinction rates
- Gobioidei229crabsbammextinctexpincrease.pdf: Diversification rate functions over a range of exponentially increasing extinction rates
- Gobioidei229crabsbammextinctincrease.pdf: Diversification rate functions over a range of linearly increasing extinction rates
- Gobioidei229crabsbammextinctoppositeadd.pdf: Diversification rate functions over a range of values for extinction rates that vary in opposition to the speciation rates.
TESS subfolder contains output plots from the CRABS analysis of TESS output:
- Gobioidei229crabstessmodels.pdf: Diversification rates in the reference model
- Gobioidei229crabstessextinct.pdf: Diversification rate functions over a range of constant extinction rates
- Gobioidei229crabstessextinctexpincrease.pdf: Diversification rate functions over a range of exponentially increasing extinction rates
- Gobioidei229crabstessextinctincrease.pdf: Diversification rate functions over a range of linearly increasing extinction rates
- Gobioidei229crabstessextinctoppositeadd.pdf: Diversification rate functions over a range of values for extinction rates that vary in opposition to the speciation rates.
IQtree.zip folder:
- Gobiiformes-clean-75p-270.phy is the 75% complete aligned UCE matrix in phylip format
- Gobiiformes-clean-75p-270charsets.nex is the UCE matrix partition file in nexus format
Partitioned output folder contains the IQtree output for the partitioned analysis, including three files:
- Gobiiformes75-270part.merge.iqtree is the results summary
- Gobiiformes75-270part.merge.treefile is the maximum likelihood tree with bootstrap branch supports
- Gobiiformes75-270part.merge.contree is the 50% majority rule consensus tree of the bootstrap replicates
Unpartitioned output folder contains the IQtree output for the unpartitioned analysis, including three files:
- Gobiiformes75-270.merge.iqtree is the results summary
- Gobiiformes75-270.merge.treefile is the maximum likelihood tree with bootstrap branch supports
- Gobiiformes75-270.merge.contree is the 50% majority rule consensus tree of the bootstrap replicates
SSE.zip folder:
Gobiiformes75p229calibrated.nwk is the fossil-calibrated BEAST tree for Gobioidei in newick format
Gobiiformes229SSEtraits.csv is the input file for traits scored in the SSE analyses. The columns are:
- Taxon: the currently valid taxonomic name
- Disc: whether the ventral pelvic disc is absent (state 0) or present (state 1)
- Cryptobenthic: whether the species is not cryptobenthic (state 0) or cryptobenthic (state 1)
- Reef: whether the species is not reef associated (state 0) or reef associated (state 1)
- MAREUFW: salinity tolerance of the species, being exclusively marine (state 1), eurhaline (state 2), or exclusively freshwater (state 3).
- EUMARZERO: salinity tolerance recoded as a binary character, with marine and euryhaline (state 0) contrasted to exclusively freshwater (state 1).
- MARZERO: salinity tolerance recoded as a binary character, with exclusively marine (state 0) contrasted to euryhaline or freshwater (state 1).
TACT.zip folder:
Gobioideispecieslist.csv is the list of all Gobioidei species in a taxonomy, used to construct Gobioidei.taxonomy.tre. Taxonomic assignments are given in five columns representing Order, Family, Lineage, Genus, and the Genus species combination.
Gobioidei.taxonomy.tre is the partially resolved phylogeny that reflects the taxonomic groups in Gobioideispecieslist.csv
Gobiiformes229TACT.nwk is the incompletely sampled UCE tree with corrected taxon names, used as a backbone phylogeny for TACT.
trees subfolder contains 100 completely sampled TACT trees for Gobioidei in newick and nexus format, used to validate the BAMM and TESS analyses on the UCE tree.
TESS.zip folder:
Gobiiformes75p229calibrated.nwk is the fossil-calibrated BEAST tree for Gobioidei in newick format
Gobiiformes75p229TESSspeciationextinction.pdf shows the graphs of results from the TESS CoMET analysis, including inferred speciation rates and shift times, extinction rates and shift times, and mass extinction Bayes factors and times for Gobioidei.
Gobiiformes75p229TESSnetdiversification.pdf shows the graph of results from the TESS CoMET analysis of inferred net diversification rates through time for Gobioidei.
Gobiiformes75p229TESSpriorsspeciationextinction.pdf shows the graphs of results from the TESS CoMET analysis using TESS-inferred priors, including inferred speciation rates and shift times, extinction rates and shift times, and mass extinction Bayes factors and times for Gobioidei.
Gobiiformes75p229TESSpriorsnetdiversification.pdf shows the graph of results from the TESS CoMET analysis using TESS-inferred priors of inferred net diversification rates through time for Gobioidei.
Gobiiformes229TESSnoextinctrates.pdf shows the graphs of results from the TESS CoMET analysis without mass extinctions, including inferred speciation rates and shift times and extinction rates and shift times.
Code/software
Gobiiformes_script_R.R is the R code used to run the TESS CoMET, CRABS, BiSSE, MuSSE, and HISSE analyses
