Ciliate SSU-rDNA reference alignments and trees for phylogenetic placements of metabarcoding data

File S1: Description of the relationships within the main lineages in the reference trees Ľubomír Rajter, Micah Dunthorn

KAYORELICTEA encompasses around 170 morphospecies belonging to 6 families (Cryptopharyngidae, Geleiidae, Kentrophoridae, , Trachelocercidae, and Wilbertomorphidae) (Lynn 2008). Those families represent the main molecular clades in recent SSU-rRNA karyorelictean phylogenetic studies (e.g., Campello-Nunes et al., 2015; Xu et al., 2015, 2013a, 2013b; Yan et al., 2016). The only exception is Cryptopharyngidae with an unknown evolutionary position, as molecular data for this family are missing. This class was monophyletic in the unconstrained tree with high statistical support. We also recognized all five main karyorelictean clades/families with the consistent topology in all our trees.

HETEROTRICHEA is a diverse class consisting of ten monophyletic families (Shazib et al. 2014). Each of this family represents well-supported distinct clade in heterotrichean SSU- rRNA phylogeny (Schmidt et al. 2007, Fernandes et al. 2016, Chen et al. 2020). In our trees were all the heterotrichean families also revealed monophyletic. Specifically, the family Peritromidae branched off first in all our trees. Climacostomidae was a sister to Spirostomidae in all trees except the constrained masked tree, where Climacostomidae was a sister to Condylostomatidae. The next branched off a lineage that united four smaller heterotrichean families (Fabreidae, Folliculinidae, Gruberiidae, and Maristentoridae) with two crown species-rich families Blepharismidae and Stentoridae. This topology was also revealed by Shazib et al. (2014) based on multiple genes phylogenetic analyses.

SPIROTRICHEA is one of the most diverse Ciliophora class containing three main clades – Euplotia, Hypotrichia, Oligotrichea. Relationships between these three clades had rather paraphyletic character in all our trees. Specifically, the clade Euplotia contained Hypotrichia+Oligotrichea and Hypotrichia contained monophyletic Oligotrichea. This topology is also common in other spirotrichean phylogenetic studies (Kim et al. 2005, Lyu et

1

Ciliate SSU-rDNA reference alignments and trees for phylogenetic placements of metabarcoding data al. 2018, Song et al. 2015). Inconsistent is only the clade Euplotia which is mostly monophyletic (Gao et al 2016), but becomes paraphyletic if sequences from Diophrys genera are included (Gao et al. 2014, this study). Nevertheless, the nodes defining the relationships between these main clades had only low statistical support in the unconstrained tree. Only well supported deep-nodes were for the Oligotrichean’s monophyly (96 bootstrap value) and monophyly of the whole Spirotrichea class (88 bootstrap value).

Except for these three core spirotrichean clades, we also included in Spirotrichea two other groups – Licnophoridae and Kitrichidae – following the Adl et al. (2019) system. The affinity of these groups to the class Spirotrichea is still discussed. In all our constrained trees, Licnophoridae and Kitrichidae were monophyletic and branched off first within sprirotricheans, although they were forced to be within the class. In the unconstrained tree, Kitrichidae was in the vicinity of the class Spirotrichea, but Licnophoridae branched outside of Spirotrichea as the earliest branch of the whole SAL lineage.

ARMOPHOREA (incl. Cariacotrichea, Muranotrichea and Parablepharismea). This class encompassed four main lineages – Caenomorphidae, Odontostomatida, Metopida, and Clevelandellida – based on the Adl et al. (2019) system. As in our trees are also in the close vicinity of Armophorea clades Cariacotrichea, Muranotrichea, and Parablepharismea, we discuss those here as well. Cariacotrichea always grouped with Parablepharismea and Muranotrichea, creating a sister clade to Armophorea in all our trees. This grouping was also inferred in Rotterová et al. (2020), where the authors described these two new classes: Muranotrichea and Parablepharismea.

Within the class Armophorea, Caenomorphidae branched off first in constrained unmasked alignment or they represented a sister group to the Parablepharismidae+Muranotrichidae+Cariacotrichea clade in our masked and V4 unmasked alignments. In the unconstrained tree, Caenomorphids were sister to Litostomatea but only with low statistical support. But this is not unusual as caenomorphids branch outside of the Armophorea in all recent molecular studies (Campello-Nunes et al. 2020, Fernandes et al. 2018 and reference therein, Rotterová et al. 2020), so their definitive phylogenetic home has not been resolved yet. The phylogenetic position of Odontostomatida is stable in all our trees, always as sister to Armophorea within APM-clade, same as in other molecular studies (e.g.,

2

Ciliate SSU-rDNA reference alignments and trees for phylogenetic placements of metabarcoding data

Campello-Nunes et al. 2020). It’s worth mentioning that Fernandes et al. (2018) recently proposed odontostomatids as separate class, but we follow the Adl et al. (2018) system where they belong to Armophorea. Metopida contained all its families (Apometopidae, Metopidae, and Tropidoatractidae), but it is paraphyletic as it also encompassed monophyletic Clevelandellida. This result corresponds with other recent molecular studies (e.g., Bourland et al. 2020).

LITOSTOMATEA contains two free-living subclasses (Rhynchostomatia and Haptoria) and an endosymbiotic subclass (Trichostomatia). But many litostomatean taxa also have an uncertain phylogenetic position, and intraclass phylogeny is not entirely resolved yet (Vďačný et al. 2014). The phylogenetic position of enigmatic Mesodiniidae is discussed in the incertae sedis part below.

In all the trees, branched off the first two long-branched litostomatean taxa: Helicoprorodon and then Chaenea. Lacrymariida and Rhynchostomatia clustered together in all the trees except the V4 unmasked alignment in which Lacrymariida clustered with Homalozoon, and Rhynchostomatia represented a separate clade. But the bootstrap support for Lacrymariida+Rhynchostomatia grouping is low (39) in the unconstrained tree. Didiniida and Pleurostomatida grouped together in all trees except the masked alignment tree in which are both as distinct clades, but the statistical support for sister relationship of these two clades is low (33) in the unconstrained tree. Finally, two haptorian species-rich orders Spathidiida and Haptorida are mixed and contained endosymbiotic subclass Trichostomatia as in all the other phylogenetic studies (e.g., Vďačný et al. 2014).

Within Trichostomatia branched off first two entodiniomorphids groups (Spirodiniidae and Buetscheliida) and then Vestibulifera that contained Macropodiniida. The crown trichostomatian groups consisted of the rest of the entodiniomorphids (i.e., Isotrichidae, Blepharocorythina, Polydiniellidae, Spirodiniida, Troglodytellidae, Cycloposthiidae, and Ophyoscolecidae).

COLPODEA. All its four main lineages - Platyophryida, Bursariomorphida, Cyrtolophosidida, and Colpodida - were recognized in all constrain trees. The only difference

3

Ciliate SSU-rDNA reference alignments and trees for phylogenetic placements of metabarcoding data was in the topology. Platyophryida branched off first, then Bursariomorphida as a sister to Cyrtolophosidida+Colpodida in both unmasked alignments. On the contrary, Bursariomorphida branched off first, then Platyophryida as a sister to Cyrtolophosidida+Colpodida in masked and V4 unmasked alignments. Though both topology scenarios are plausible (Vďačný and Foissner 2019), platyophryids-branching-first topology has slightly more molecular support (Rajter et al. 2020).

OLIGOHYMENOPHOREA. Oligohymenophorea is a diverse class that encompasses six main lineages: Apostomatia, Astomatia, Hymenostomatia, Peniculia, Peritrichia, and Scuticociliatia. Relationships among these main lineages varied based on the masking strategy used. We can also see in the unconstrained trees very low bootstrap support for nodes defining these relationships. The monophyly of Peniculia had full-statistical support but excluding Urocentrum turbo that cluster outside of this clade in all our trees. Peritrichia is bifurcated to its two primary groups – Mobilida and Sessilida – in all our trees. The statistical supports for the monophyly of both groups are high (96, 100) although the support for the monophyly of Peritrichia itself is low (65). This could by the reason that in the masked alignment are Mobilida and Sessilida not grouped together. Astomatia is monophyletic with moderate statistical support (90) although the Haptophrya planarium and Durchoniella brasili may clustered outside (see the trees from masked and V4 unmasked alignments). Apostomatia and Hymenostomatia were always monophyletic with full-statistical support in unconstrained tree. Finally, Scuticociliatia was non-monophyletic and got only low statistical support in the unconstrained tree. But their three main lineages – Loxocephalida, Pleuronematida and Philasterida – were recognized in all our trees and had moderate to full support (91, 95, and 100) in the unconstrained tree.

NASSOPHOREA. Nassophorea represents a smaller group compare to other ciliate classes and contains two main lineages – Nassulida and Microthoracida. We forced these two lineages to cluster together as they always represent separate clades in the SSU-rRNA phylogenies. In our masked and V4 unmasked alignments was Nassulida paraphyletic as contained monophyletic Microthoracida. In unmasked alignment were both lineages monophyletic. In the unconstrained tree, Microthoracida had full statistical support but

4

Ciliate SSU-rDNA reference alignments and trees for phylogenetic placements of metabarcoding data branched outside of Nassophorea as sister to Phyllopharyngea, same as in other phylogenetic studies (e.g., Pan et al. 2019). PHYLLOPHARYNGEA Phyllopharyngea is a species-rich class but can be split into four main lineages: Synhymeniida, Rhynchodia, Suctoria, and Cyrtophoria. All four lineages were monophyletic with identical topologies in all our trees. Synhymeniida branched off first and had a full statistical support in the unconstrained tree. Rhynchodia, Suctoria, and Cyrtophoria grouped together with the full statistical support. Rhynchodia, which is represented only by the one sequence (Hypocoma acinetarum), was a sister to Suctoria. And finally, the monophyly of Cyrtophoria was moderately suppoted (88) in the unconstrained tree. This lineage is further divided into two subgroups: paraphyletic Chlamidodontida that contained monophyletic Dysteriida. All of these evolutionary relationships correspond with those from other phyllopharyngean phylogenetic studies (e.g., Wang et al. 2017).

PLAGIOPYLEA. Plagiopylea represents a small ciliophoran class and contains four main lineages/families: Epalxellidae, Trimyemidae, Sonderiidae, Plagiopylidae. All four were recognized with high to full-statistical support. The topology was consistent in all our trees: Epalxellidae branched off first. Sonderiidae+Plagiopylidae grouped together as a sister clade to Trimyemidae.

PROSTOMATEA. Prostomatea is represented here by its eight families (Balanionidae, Colepidae, Cryptocaryonidae, Holophryidae, Metacystidae, Placidae, Plagiocampidae, and Urotrichidae), which were recognized in all trees. Topologies slightly differ based on the alignment used and also deep-nodes within this class had only low statistical support in the unconstrained tree, so relationships among those families are difficult to reconstructed. Only exceptions are Holophryidae+Metacystidae that always cluster together and had a high bootstrap value (97), and Cryptocaryonidae+Plagiocampidae that also always cluster together and had a full statistical support. These results correspond with other prostomatea studies (e.g., Zhang et al. 2014).

5

Ciliate SSU-rDNA reference alignments and trees for phylogenetic placements of metabarcoding data

References

Adl SM, Bass D, Lane CE, Lukeš J, Schoch CL, Smirnov A, Agatha S, Berney C, Brown MW, Burki F, Cárdenas P, Čepička I, Chistyakova L, Campo J del, Dunthorn M, Edvardsen B, Eglit Y, Guillou L, Hampl V, Heiss AA, Hoppenrath M, James TY, Karnkowska A, Karpov S, Kim E, Kolisko M, Kudryavtsev A, Lahr DJG, Lara E, Gall LL, Lynn DH, Mann DG, Massana R, Mitchell EAD, Morrow C, Park JS, Pawlowski JW, Powell MJ, Richter DJ, Rueckert S, Shadwick L, Shimano S, Spiegel FW, Torruella G, Youssef N, Zlatogursky V, Zhang Q (2019) Revisions to the Classification, Nomenclature, and Diversity of Eukaryotes. Journal of Eukaryotic Microbiology 66: 4–119. https://doi.org/10.1111/jeu.12691 Bourland W, Rotterová J, Čepička I (2020) Description of Three New Genera of Metopidae (Metopida, Ciliophora): Pileometopus gen. nov., Castula gen. nov., and Longitaenia gen. nov., with Notes on the Phylogeny and Cryptic Diversity of Metopid . Protist 171: 125740. https://doi.org/10.1016/j.protis.2020.125740 Campello-Nunes PH, Fernandes N, Schlegel M, Silva-Neto ID (2015) Description and phylogenetic position of Corlissina maricaensis gen. nov., sp. nov. (, Geleiidae), a novel interstitial ciliate from Brazil, with redefinition of the family Geleiidae. International Journal of Systematic and Evolutionary Microbiology 65: 4297–4308. https://doi.org/10.1099/ijsem.0.000579 Campello-Nunes PH, Fernandes NM, Szokoli F, Fokin SI, Serra V, Modeo L, Petroni G, Soares CAG, Paiva T da S, Silva-Neto ID da (2020) Parablepharisma (Ciliophora) is not a Heterotrich: A Phylogenetic and Morphological Study with the Proposal of New Taxa. Protist 171: 125716. https://doi.org/10.1016/j.protis.2020.125716 Chen X, Li J, Xu K (2020) Multigene-based phylogeny analyses of the controversial family Condylostomatidae (Ciliophora, Heterotrichea). Zoologica Scripta 49: 250–264. https://doi.org/10.1111/zsc.12383 Fernandes NM, Paiva T da S, da Silva-Neto ID, Schlegel M, Schrago CG (2016) Expanded phylogenetic analyses of the class Heterotrichea (Ciliophora, ) using five molecular markers and morphological data. Molecular Phylogenetics and Evolution 95: 229–246. https://doi.org/10.1016/j.ympev.2015.10.030 Fernandes NM, Vizzoni VF, Borges B do N, Soares CAG, da Silva-Neto ID, Paiva T da S (2018) Molecular phylogeny and comparative morphology indicate that odontostomatids (Alveolata, Ciliophora) form a distinct class-level taxon related to Armophorea. Molecular Phylogenetics and Evolution 126: 382–389. https://doi.org/10.1016/j.ympev.2018.04.026 Gao F, Katz LA (2014) Phylogenomic analyses support the bifurcation of ciliates into two major clades that differ in properties of nuclear division. Molecular Phylogenetics and Evolution 70: 240–243. https://doi.org/10.1016/j.ympev.2013.10.001 Gao F, Warren A, Zhang Q, Gong J, Miao M, Sun P, Xu D, Huang J, Yi Z, Song W (2016) The All- Data-Based Evolutionary Hypothesis of Ciliated Protists with a Revised Classification of the Phylum Ciliophora (Eukaryota, Alveolata). Scientific Reports 6: 24874. https://doi.org/10.1038/srep24874 Kim JS, Jeong HJ, Strüeder-Kypke MC, Lynn DH, Kim S, Kim JH, Lee SH (2005) Parastrombidinopsis shimi n. gen., n. sp. (Ciliophora: Choreotrichia) from the coastal waters of Korea: morphology and small subunit ribosomal DNA sequence. The Journal of Eukaryotic Microbiology 52: 514–522. https://doi.org/10.1111/j.1550-7408.2005.00062.x 6

Ciliate SSU-rDNA reference alignments and trees for phylogenetic placements of metabarcoding data

Lynn D (2008) The Ciliated Protozoa: Characterization, Classification, and Guide to the Literature. 3rd ed. Springer Netherlands. https://doi.org/10.1007/978-1-4020-8239-9

Lyu Z, Wang J, Huang J, Warren A, Shao C (2018) Multigene-based phylogeny of Urostylida (Ciliophora, Hypotrichia), with establishment of a novel family. Zoologica Scripta 47: 243– 254. https://doi.org/10.1111/zsc.12267 Pan B, Chen X, Hou L, Zhang Q, Qu Z, Warren A, Miao M (2019) Comparative Genomics Analysis of Ciliates Provides Insights on the Evolutionary History Within “Nassophorea–Synhymenia– Phyllopharyngea” Assemblage. Frontiers in Microbiology 10. https://doi.org/10.3389/fmicb.2019.02819 Rajter Ľ, Ewers I, Graupner N, Vďačný P, Dunthorn M (2021) Colpodean ciliate phylogeny and reference alignments for phylogenetic placements. European Journal of Protistology 77: 125747. https://doi.org/10.1016/j.ejop.2020.125747 Rotterová J, Salomaki E, Pánek T, Bourland W, Žihala D, Táborský P, Edgcomb VP, Beinart RA, Kolísko M, Čepička I (2020) Genomics of New Ciliate Lineages Provides Insight into the Evolution of Obligate Anaerobiosis. Current Biology 30: 2037-2050.e6. https://doi.org/10.1016/j.cub.2020.03.064 Schmidt SL, Foissner W, Schlegel M, Bernhard D (2007) Molecular phylogeny of the Heterotrichea (Ciliophora, Postciliodesmatophora) based on small subunit rRNA gene sequences. The Journal of Eukaryotic Microbiology 54: 358–363. https://doi.org/10.1111/j.1550- 7408.2007.00269.x Shazib SUA, Vďačný P, Kim JH, Jang SW, Shin MK (2014) Phylogenetic relationships of the ciliate class Heterotrichea (Protista, Ciliophora, Postciliodesmatophora) inferred from multiple molecular markers and multifaceted analysis strategy. Molecular Phylogenetics and Evolution 78: 118–135. https://doi.org/10.1016/j.ympev.2014.05.012 Song W, Zhao X, Liu W, Hu X, Al-Farraj SA, Al-Rasheid KAS, Song W, Warren A (2015) Biodiversity of oligotrich ciliates in the South China Sea: description of three new Strombidium species (Protozoa, Ciliophora, Oligotrichia) with phylogenetic analyses. Systematics and Biodiversity 13: 608–623. https://doi.org/10.1080/14772000.2015.1081992 Vďačný P, Foissner W (2019) Re-analysis of the 18S rRNA gene phylogeny of the ciliate class Colpodea. European Journal of Protistology 67: 89–105. https://doi.org/10.1016/j.ejop.2018.11.003 Vďačný P, Breiner H-W, Yashchenko V, Dunthorn M, Stoeck T, Foissner W (2014) The chaos prevails: molecular phylogeny of the Haptoria (Ciliophora, Litostomatea). Protist 165: 93–111. https://doi.org/10.1016/j.protis.2013.11.001 Wang P, Wang Y, Wang C, Zhang T, Al-Farraj SA, Gao F (2017) Further consideration on the phylogeny of the Ciliophora: Analyses using both mitochondrial and nuclear data with focus on the extremely confused class Phyllopharyngea. Molecular Phylogenetics and Evolution 112: 96–106. https://doi.org/10.1016/j.ympev.2017.04.018 Xu Y, Li J, Song W, Warren A (2013a) Phylogeny and establishment of a new ciliate family, Wilbertomorphidae fam. nov. (Ciliophora, Karyorelictea), a highly specialized taxon represented by Wilbertomorpha colpoda gen. nov., spec. nov. The Journal of Eukaryotic Microbiology 60: 480–489. https://doi.org/10.1111/jeu.12055

7

Ciliate SSU-rDNA reference alignments and trees for phylogenetic placements of metabarcoding data

Xu Y, Gao S, Hu X, Al-Rasheid KAS, Song W (2013b) Phylogeny and systematic revision of the karyorelictid genus (Ciliophora, Karyorelictea) with descriptions of two new species. European Journal of Protistology 49: 438–452. https://doi.org/10.1016/j.ejop.2012.12.001 Xu Y, Pan H, Miao M, Hu X, Al-Farraj SA, Al-Rasheid KAS, Song W (2015) Morphology and phylogeny of two species of (Ciliophora, Karyorelictea), with description of a new subspecies, Loxodes striatus orientalis subsp. n. The Journal of Eukaryotic Microbiology 62: 206–216. https://doi.org/10.1111/jeu.12162 Yan Y, Xu Y, Al‐Farraj SA, Al‐Rasheid KAS, Song W (2016) Morphology and phylogeny of three trachelocercids (Protozoa, Ciliophora, Karyorelictea), with description of two new species and insight into the evolution of the family Trachelocercidae. Zoological Journal of the Linnean Society 177: 306–319. https://doi.org/10.1111/zoj.12364 Zhang Q, Yi Z, Fan X, Warren A, Gong J, Song W (2014) Further insights into the phylogeny of two ciliate classes Nassophorea and Prostomatea (Protista, Ciliophora). Molecular Phylogenetics and Evolution 70: 162–170. https://doi.org/10.1016/j.ympev.2013.09.015

8