Phylogenomic analyses support the bifurcation of ciliates into two major clades that differ in properties of nuclear division

Abstract

Ciliates are a diverse assemblage of eukaryotes that have been the source of many discoveries including self-splicing RNAs, telomeres and trans-splicing. While analyses of ciliate morphology have given rise to robust hypotheses on relatively shallow level relationships, the deeper evolutionary history of ciliates is largely unknown. This is in part because studies to date have focused on only a single locus, small subunit ribosomal DNA (SSU-rDNA). In the present study, we use a taxon-rich strategy based on multiple loci from GenBank and recently completed transcriptomes to assess deep phylogenetic relationships among ciliates. Our phylogenomic data set includes up to 537 taxa, all of which have been sampled for SSU-rDNA and a subset of which have LSU-rDNA and up to 7 protein-coding sequences. Analyses of these data support the bifurcation of ciliates as suggested by SSU-rDNA, with one major clade defined by having somatic macronuclei that divide with intranuclear microtubules (Intramacronucleata) and the other clade containing lineages that either divide their macronuclei with microtubules external to the macronucleus or are unable to divide their macronuclei (Postciliodesmatophora). These multigene phylogenies provide a robust framework for interpreting the evolution of innovations across the ciliate tree of life.

Keywords: Ciliophora; Intramacronucleata; Macronucleus; Phylogenomic analysis; Postciliodesmatophora.

Figure 1

Maximum likelihood (ML) tree reconstructed using 537 ciliates and all 9 genes (all: 9; SSU-rDNA, LSU-rDNA plus 7 protein genes) yields a well-resolved bifurcation of ciliates phylogenies, Intramacronucleata (in maroon) and Postciliodesmatophora (in blue). Numbers at the nodes represent ML support values. The scale bar corresponds to 0.05 expected substitutions per site.

Figure 2

Maximum likelihood (ML) tree reconstructed using 68 ciliates that had at least 2 proteins of the 9 genes (2P:7; SSU-rDNA, LSU-rDNA plus 7 protein genes; see Table 1). Numbers at the main nodes represent ML support values. The scale bar corresponds to 0.05 expected substitutions per site.