Are interphylum spiralian relationships resolvable?

  1. Ana Serra Silva
  2. Maximilian J Telford  Is a corresponding author
  1. Centre for Life’s Origins and Evolution, Department of Genetics, Evolution and Environment, University College London, United Kingdom
7 figures, 12 tables and 2 additional files

Figures

There are 15 unrooted trees linking the 5 major spiralian phyla and each has been supported in at least one publication (Supplementary file 1).

Tree numbering throughout this study refers to the trees in this figure. For all 105 rooted trees, the rooting schemes follow that depicted in unrooted T1.

Branches linking spiralian phyla are shorter than both the short, disputed deuterostome branch (Kapli et al., 2021; Serra Silva et al., 2025) and the short branch leading to Lophotrochozoa (Telford et al., 2015).

(A and B) are data from SERR. (C and D) are data from MARL. Compared to the site-homogeneous model (A, C), internal branches (orange) become shorter under the site-heterogeneous model while the long branches leading to Platyhelminthes (blue) and Mollusca (dark green) become even longer (B, D). Circles on the y-axis indicate the median branch lengths of Deuterostomia (D), Ecdysozoa (E), Lophotrochozoa (L), and Protostomia (P). Triangles correspond to the median spiralian internal branches (Int), values given in Appendix 2—table 1. Scales are identical on either side of the y-axis breaks.

Discordance-aware summary tree analyses show extensive gene-tree conflict within Spiralia.

For SERR (A), the Astral analysis recovers all spiralian phyla as monophyletic but no internal Spiralia branch with quartet scores (QS) over 0.50. The exception is the branch separating Platyhelminthes + outgroups from the remaining Spiralia phyla; however, Platyhelminthes are recovered within paraphyletic Gnathifera. For MARL (B), the Astral analysis yields very similar results, except that Annelida is not monophyletic. Branch-lengths are in coalescent units and monophyletic phyla have been collapsed, with the QS (Q1/Q2/Q3) for the branch subtending each phylum listed underneath the phylum names. Detailed QS can be found in the Astral.zip folder in the Dryad data repository (https://doi.org/10.5061/dryad.280gb5n3j).

Relative median resampling estimated log-likelihood (RELL) pseudo-bootstraps (A–D) and randomised taxon-jackknives (E–H) log-likelihood ranks for all rooted trees, see Appendix 2—table 3 and Appendix 2—table 5 for all ranks.

Highest- and lowest-ranked topologies labelled for all analyses, along with every seventh topology from highest- to lowest-ranked. Topologies ranked by their difference in log-likelihood from the mean log-likelihood of each RELL or jackknife distribution. The RELL analyses of both datasets under the site-homogeneous LG + F + G4 model (A, B) occupy the widest ranges of log-likelihoods, which is consistent with the statistically significant differences found by the Kruskal–Wallis (KW) test. The RELL analyses under the site-heterogeneous EDM + F + G4 model (C, D) occupy a much narrower range of log-likelihoods. In practical terms, the wider the log-likelihood range, the more dissimilar scores are and the more we would expect to recover significantly different scores between topologies. The jackknife analyses show a similar trend, wider range of median log-likelihood scores under LG + F + G4 (E, F) than under EDM + F + G4 (G, H); however, the log-likelihood ranges of the jackknives under LG + F + G4 are very similar to that occupied by RELL under EDM + F + G4, consistent with the lack of statistically significant differences between topologies in the jackknife analyses.

With simulated data inadequate models incorrectly support a platyhelminth root.

Using data simulated on an unresolved Spiralia tree, we ranked rooted spiralian topologies by difference to mean log-likelihood score. Highest- and lowest-ranked topologies are labelled for all analyses, along with every seventh-topology from highest- to lowest-ranked. For datasets simulated under EDM (complex) but scored under LG (simple) the 15 topologies with Spiralia rooted on Platyhelminthes are highest ranked (A, B). The plots for the ‘correctly’ specified models (C, D) show the expected approximately equal log-likelihood scores for all 105 rooted topologies. See Appendix 2—table 7 for complete topology rankings.

Figure 6 with 1 supplement
Most fixed-taxon resampling estimated log-likelihood (RELL) bootstrap and randomised taxon-jackknife replicates, under both site-homogeneous (LG) and site-heterogeneous (EDM) models, recover unrooted T4 as the best-scoring topology.

Heights of bars correspond to the number of replicates where each topology had the highest log-likelihood score.

Figure 6—figure supplement 1
Relative median resampling estimated log-likelihood (RELL) pseudo-bootstraps (A–D) and randomised taxon-jackknives (E–H) log-likelihood ranks for all unrooted trees, see Appendix 2—table 8 and Appendix 2—table 10 for all ranks.

Topologies ranked by their difference in log-likelihood to the mean log-likelihood of each RELL or jackknife distribution. The RELL analyses of both datasets under the site-homogeneous LG + F + G4 model (A, B) occupy the widest ranges of log-likelihoods, which is consistent with the statistically significant differences found by the Kruskal-Wallis test. The RELL analyses under the site-heterogeneous EDM + F + G4 model (C, D) occupy a much narrower range of log-likelihoods. The jackknife analyses show a similar trend, wider range of median log-likelihood scores under LG + F + G4 (E, F) than under EDM + F + G4 (G, H); however, the log-likelihood ranges of the jackknives under LG + F + G4 are very similar to that occupied by RELL under EDM + F + G4, consistent with the lack of statistically significant differences between topologies in the jackknife analyses.

With simulated data, misspecified models recover unrooted T4 within the set of highest-scoring topologies.

Using data simulated on an unresolved Spiralia tree we ranked unrooted spiralian topologies by difference to mean log-likelihood score. For datasets simulated under EDM (complex) but scored under LG (simple), three topologies (uT4, uT5, and uT14) are within the highest-ranked topologies for both empirical datasets (A, B). The plots for the ‘correctly’ specified models (C, D) show the expected approximately equal log-likelihood scores for all 15 unrooted topologies.

Tables

Key resources table
Reagent type (species) or resourceDesignationSource or referenceIdentifiersAdditional information
OtherAlignments, phylogenetic analyses, simulations, and ancillary scriptsDRYAD; https://doi.org/10.5061/dryad.280gb5n3jPlease see Materials and methods for data processing details
Software and algorithmsIQ-TREE2 (v.2.3.5)Minh et al., 2020; http://www.iqtree.org/
Software and algorithmsEDClusterSchrempf et al., 2020; Schrempf, 2020; https://github.com/dschrempf/EDCluster
Software and algorithmsASTRAL-III (v.5.1.3)Zhang et al., 2017a; Zhang et al., 2017b; https://github.com/smirarab/ASTRAL
Software and algorithmsMAFFT (v.7.455)Katoh and Standley, 2013; https://mafft.cbrc.jp/alignment/software/
Software and algorithmstrimAl (v.1.4.rev15)Capella-Gutiérrez et al., 2009; https://vicfero.github.io/trimal/
Software and algorithmsAmino acid diversity testsGiacomelli et al., 2022; Giacomelli et al., 2025; Giacomelli, 2026; https://github.com/mgiacom/tardigrades_catpmsf
Software and algorithmsAdditional data processing and analysis scriptsDRYAD; https://doi.org/10.5061/dryad.280gb5n3j
Appendix 2—table 1
Median branch lengths for the 105 rooted topologies.

Phylum-level branch lengths shown in Figure 2.

BranchSERR + LGSERR + EDMMARL + LGMARL + EDM
Deuterostomia0.031390.028120.023050.01782
Ecdysozoa0.047620.062570.038350.03941
Lophotrochozoa0.024400.026070.023820.02613
Protostomia0.064200.093250.048520.06641
Internal0.016160.013930.013160.009462
Appendix 2—table 2
Fixed-taxon RELL pseudo-bootstraps rankings for 105 rooted topologies.
TopologySERR + LGSERR + EDMMARL + LGMARL + EDM
140342612
25328
33835157
473774627
571573821
692817777
7101916180
856668024
93531486
1027224011
1141425662
1219215155
1330246657
14221340
1572706823
165036395
173323349
1841610
1943254219
2077547071
2164458267
223938374
232218221
241712202
2521163218
2618131215
2745263549
2864529
2957625333
3046475237
3134282514
322319143
332014813
3455332951
35105330
3697966959
3790867161
3875603331
3978683152
4048401020
4149392126
4276122
43100929572
44897810179
4570534943
46138928
4774515953
4893696568
4967464132
50821029999
518110085103
5247886787
5329494544
5416325883
5526437288
56171966
571021047386
581031056094
5966954469
6052755070
6136551634
6237522842
6339435
649810396101
659610186104
6660906490
678201165
6854896384
6951635546
7059769091
7165647441
7211111725
7363617945
7458445417
7584597538
76998210598
7710593102102
7888839474
7914172758
80857910082
81918010393
8294719776
8386657856
84104879896
8587988995
8615272478
87839792100
88959988105
8980948792
9069728154
91798410497
9225293016
9331415747
9432563639
9542744348
9624301836
9762674785
981215763
9953506250
10076738489
10161588364
10268859375
10328377673
10444489181
1059102360
Appendix 2—table 3
Fixed-taxon RELL pseudo-bootstraps best-ranked topology frequencies for 105 rooted topologies.

Omitted topologies were never recovered as the best-ranked tree.

TopologySERR + LGSERR + EDMMARL + LGMARL + EDM
28.411.9627.620.43
30000.59
148.0915.7400
1811.2351.781.720.34
220000.04
2300057.2
2400.28036.52
2600.10.010.05
280.413.84.720
290000.01
310000.67
320003.53
3300.040.360.5
350.266.1513.150
380000.01
390000
400000.01
420.345.2650.330.05
460.011.4900.02
5663.532.2300
610000.02
637.650.542.030
670.03000
7200.0700
920000.01
98000.060
1050.050.500
Appendix 2—table 4
Randomised taxon-jackknife rankings for 105 rooted topologies.
TopologySERR + LGSERR + EDMMARL + LGMARL + EDM
139183619
2715324
32214610
463643013
554132817
680547583
7100886572
858458511
93312526
10202354
1148415053
121574550
1334167769
14491970
1577587112
165211487
17244265
1823114
193554320
2091617976
2167279177
225533383
23386222
24191111
252783123
262310916
2756374760
289241254
2972596326
3061406440
314322209
323131158
3336651018
3475874467
3513801358
3693816946
3779637863
3864552525
3992922742
403056722
4142463228
42842437
4396769671
44907010290
4560385844
4611261652
4768476757
48101937278
4969535730
50821009985
51831048694
5229695456
5321213932
5412174966
5528307682
561201784
57951016061
58991055575
5946772338
6037504255
611639533
6226342443
63332250
6497989287
651021028895
6644714668
67660865
6845685164
6950355639
7049439089
7170757034
7210572141
7366796145
7453255321
7571236236
769862104105
7710589100102
7887849874
7914784191
8081869393
818474103104
8288669579
8386728049
84104959792
8589978796
86179434100
87859982101
8810310384103
8973837381
9076448148
91785210599
9240193315
9347286647
9441492929
9562854031
9632731427
9774965988
9818901886
9957676835
10094918973
10159518359
10265829462
10325297480
104514810198
1055363797
Appendix 2—table 5
Randomised taxon-jackknife best-ranked topology frequencies for 105 rooted topologies.

Omitted topologies were never recovered as the best-ranked tree.

TopologySERR + LGSERR + EDMMARL + LGMARL + EDM
215000
80005
90000
1001005
120500
145500
170500
18305250
2305030
24040055
260550
310500
320505
3300150
4000100
420050
5410000
5630000
6100100
63100200
6700100
920500
1030500
Appendix 2—table 6
Simulations on rooted Spiralia star-tree rankings for 105 rooted topologies.

For columns SERR + EDM and MARL +EDM, proportions are based on 10 alignments.

TopologySERR + LG EDM simSERR + EDMMARL + LG EDM simMARL + EDM
12492471
24247
321411834
457443744
554893646
6691037047
781915661
8958610593
9649010178
1049469380
1151618687
1223195856
1326537386
1465148
1594998199
1668577792
17479555102
181031911
19489748100
2073778898
2162739596
22103817421
2378876828
245584924
255322509
263732132
27394233101
281147616
29102296523
3090586620
3150712925
3246703235
3334501748
34387923105
35848219
3687655742
3767266227
3841832750
3943763055
4019681637
4122402282
423711
4396629875
4479459960
4542783979
46960713
4744984590
48619463103
4952806176
50913410464
51856710372
5231208467
5330357877
5417105369
5520396765
56241312
5780367151
5874756459
5929164354
6028244441
6116371945
6218132681
631132
6499288795
6598848597
6632125483
675111018
6833824794
6935438388
7036309285
71601003584
721310583
7356212543
7466145153
7582234649
76100699129
77105748270
7863884274
791463126
8059593436
8176258015
8284647639
8392529058
84104859468
8572963889
8615101115
8770562840
88101515963
8989155238
9086387933
9193338930
9275726014
9383177210
9477324017
9597494126
9640182052
97459331104
981266522
99659210073
1008810210291
101581049666
10271559757
10325276931
10427547562
10576154
Appendix 2—table 7
Simulations on rooted Spiralia star-tree best-ranked topology frequencies for 105 rooted topologies.

For columns SERR + EDM and MARL + EDM, proportions are based on 10 alignments. Omitted topologies were never recovered as the best-ranked tree.

TopologySERR + LG EDM simSERR + EDMMARL + LG EDM simMARL + EDM
1000020
1802000
2302000
2800150
330010
3500840
492000
5698000
6002000
7200020
8500020
8600020
8902000
9102000
9800020
Appendix 2—table 8
Fixed-taxon RELL pseudo-bootstraps best-ranked topology frequencies for 15 unrooted topologies.
TopologySERR + LGSERR + EDMMARL + LGMARL + EDM
10000
219203
30000
481986096
500351
60050
70000
80000
90000
100000
110000
120000
130000
140000
150000
Appendix 2—table 9
Compositional diversity model adequacy tests on the fixed-taxon data matrices.

Data simulated under the EDM + G4 model are closest in amino acid compositional diversity to the empirical data. Bold values correspond to the simulated data that most closely resemble the empirical data.

ModelAmino acid diversityZ-score
Empirical dataSimulated data
MARL
LG + G45.15925.4859 ± 0.00585955.7581
EDM + G45.15925.3082 ± 0.00577125.8182
SERR
LG + G46.29506.5769 ± 0.00728638.6836
EDM + G46.29506.1735 ± 0.006021–20.1752
Appendix 2—table 10
Randomised taxon-jackknife best-ranked topology frequencies for 15 unrooted topologies.
TopologySERR + LGSERR + EDMMARL + LGMARL + EDM
12040
2231102
30002
438614392
506261
600120
70000
830001
951151
100000
110000
120000
130000
1422001
150100
Appendix 2—table 11
Simulations on 5-taxon Spiralia star-tree best-ranked topology frequencies for 15 unrooted topologies.

For columns SERR + EDM and MARL + EDM, proportions are based on 10 alignments.

TopologySERR + LG EDM simSERR + EDMMARL +LG EDM simMARL + EDM
10010
2732000
30000
40000
5040080
6009920
70000
822000
90000
100000
110000
1202000
130000
140000
1552000

Additional files

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Ana Serra Silva
  2. Maximilian J Telford
(2026)
Are interphylum spiralian relationships resolvable?
eLife 15:RP110607.
https://doi.org/10.7554/eLife.110607.3