Identifiability of species network topologies from genomic sequences using the logDet distance

Inference of network-like evolutionary relationships between species from genomic data must address the interwoven signals from both gene flow and incomplete lineage sorting. The heavy computational demands of standard approaches to this problem severely limit the size of datasets that may be analyzed, in both the number of species and the number of genetic loci. Here we provide a theoretical pointer to more efficient methods, by showing that logDet distances computed from genomic-scale sequences retain sufficient information to recover network relationships in the level-1 ultrametric case. This result is obtained under the Network Multispecies Coalescent model combined with a mixture of General Time-Reversible sequence evolution models across individual gene trees. It applies to both unlinked site data, such as for SNPs, and to sequence data in which many contiguous sites may have evolved on a common tree, such as concatenated gene sequences. Thus under standard stochastic models statistically justifiable inference of network relationships from sequences can be accomplished without consideration of individual genes or gene trees.

Automated summary

Inference of network-like evolutionary relationships between species from genomic data must consider both gene flow and incomplete lineage sorting. Standard methods have high computational demands and limit the size of analyzed datasets. This study shows that logDet distances computed from genomic-scale sequences can efficiently recover network relationships in the level-1 ultrametric case. It applies to both unlinked site data and sequence data.