Data from: Quartet-based species tree methods enable fast and consistent tree of blobs reconstruction under the network multispecies coalescent
Data files
Jun 25, 2026 version files 9.31 GB
-
README.md
2.86 KB
-
simulation-estimated-gene-trees.tar.gz
2.13 GB
-
simulation-estimated-tobs-benchmarking.tar.gz
1.21 MB
-
simulation-true-gene-trees.tar.gz
2.68 GB
-
simulation-true-msas.tar.gz
4.50 GB
-
simulation-true-networks.tar.gz
2.07 MB
-
simulation-true-tobs.tar.gz
321.99 KB
Abstract
Hybridization between species is an important force in evolution, commonly modeled by the network multispecies coalescent. Reconstructing evolutionary histories under this model is computationally challenging, even for level-1 networks where hybridization events are isolated. Divide-and-conquer is a promising path forward, but current methods with statistical guarantees rely on an estimated tree of blobs (TOB) for the network, which compresses the non-tree-like parts into single vertices. TOB reconstruction is itself challenging, with the only available method TINNiK having time complexity O(n5 + n4k) for k genes and n species. Here, we present a new framework for scalable TOB reconstruction with statistical guarantees. Our approach operates by (1) seeking a refinement of the TOB and then (2) contracting edges in it. For step (1), we show that any optimal solution to Weighted Quartet Consensus is a TOB refinement almost surely, as the number of genes goes to infinity, motivating the use of methods, such as ASTRAL or TREE-QMC. For step (2), we show that applying the same hypothesis tests as TINNiK to just O(n) four-taxon subsets around each edge is sufficient for statistical consistency when the underlying network is level-1. Leveraging TREE-QMC for the first step gives our method time complexity O(n3k) and its name: TOB-QMC. On simulated data, TOB-QMC typically matches or exceeds TINNiK in accuracy while being more scalable. TOB-QMC also enables fast exploration of non-tree-like evolution, as demonstrated through re-analysis of three phylogenomic data sets. Lastly, our study clarifies the theoretical utility of quartet-based species tree methods in the context of hybridization, which is critical given the recent result that ASTRAL can be misleading.
This repository stores the simulate data sets generated as part of the study presenting TOB-QMC.
The simulated data are defined by the following variables:
LEVindicating the level of the true species network (1 or 2)TAXindicating the number of species (taxa) in the true network (50, 100, or 200)REPindicating the replicate identification number of the networkILSindicating the branch length scaling factor (2.0, 1.0, 0.5, 0.25) used to change the ILS level
See the paper for detailed information.
To access data, download and decompress each of the following files with command: tar -xzvf <file>.tar.gz.
simulation-true-networks.tar.gzcontains the species networks- Files are labeled
LEVEL_{LEV}/TAXA_{TAX}/network_{REP}/network_{REP}/network_{REP}_ils_{ILS}.newick - Each file is in newick format
- Files are labeled
simulation-true-tobs.tar.gzcontains the tree of blobs computed from the species network- Files are labeled
LEVEL_{LEV}/TAXA_{TAX}/network_{REP}/network_{REP}/network_{REP}_ils_1.0.blob_tree - Each file is in newick format (note that the true tree of blobs does not depend on the ILS level)
- Files are labeled
simulation-true-gene-trees.tar.gzcontains the gene trees simulated from the species network- Files are labeled
LEVEL_{LEV}/TAXA_{TAX}/network_{REP}/network_{REP}_ils_{ILS}.true_gene_trees - Each file contains 1000 gene trees, one newick string per line
- Files are labeled
simulation-true-msas.tar.gzcontains the multiple sequence alignments (MSAs) simulated from each of the 1000 gene trees- Files are labeled
LEVEL_{LEV}/TAXA_{TAX}/network_{REP}/network_{REP}_ils_{ILS}_split_MSA.tar.gz - Each file can be decompressed with the same command as above, producing 1000 files (fast format), each labeled
<GEN>.fawhereGENindicates the MSA was simulated from the gene tree on line numberGEN
- Files are labeled
simulation-estimated-gene-trees.tar.gzcontains the gene trees estimated from the MSAs (Bewick format)- Files are labeled
LEVEL_{LEV}/TAXA_{TAX}/network_{REP}/network_<REP>_ils_<ILS>.estimated_gene_trees - Each file contains up to 1000 gene trees, one newick string per line (note that gene tree estimation with IQ-TREE sometimes failed, as noted in the supplement)
- Note that the line numbers between true and estimated gene trees are not necessarily consistent
- Files are labeled
simulation-estimated-tobs-benchmarking.tar.gzcontains the tree of blobs reconstructed from estimated gene trees as part of method comparison benchmarking- Files are labeled
LEVEL_{LEV}/TAXA_{TAX}/network_{REP}/network_<REP>_ils_<ILS>.estimated_gene_trees. {MET}_ALPHA1e-7_BETA0.95.tree METindicates the method used to estimate the tree of blobs (TINNIK, TOB-QMC-fast, TOB-QMC-default). All methods were benchmarked with the same alpha value (1e-7) and the same beta value (0.95)- Each file is in newick format
- Files are labeled
