Expanded phylogenetic analyses of the class Heterotrichea (Ciliophora, Postciliodesmatophora) using five molecular markers and morphological data.

Most studies of the molecular evolution of Heterotrichea have been based solely on the 18S-rDNA gene, which were inconsistent with morphological classification. Because of the limitations of single locus phylogenies and the recurring problem of lack of resolution of deeper nodes found in previous studies, we present hypotheses of the evolution of internal groups of the class Heterotrichea based on multi-loci analyses (18S-rDNA, 28S-rDNA, ITS1-5.8S-ITS2 region, COI and alpha-tubulin) and morphological data. Phylogenetic trees from protein coding gene data are presented for Heterotrichea for the first time. Phylogenetic analyses included Bayesian inference, maximum likelihood, maximum parsimony methods, and optimal trees were statistically compared to alternative topologies from the literature. Additionally, the Bayesian concordance approach (BCA algorithm) was used to assess the concordance factor between topologies obtained from isolated analyses. Because different loci may evolve at different rates, resulting in different gene topologies, we also estimated a species tree for Heterotrichea using the STAR coalescence-based method. The results show that: (1) single gene trees are inconsistent regarding the position of some heterotrichean families; (2) the concatenation of all data in a total-evidence tree improved the resolution of deep nodes among the heterotrichean families and genera; (3) the coalescent-based species tree is consistent with phylogenies based on the 18S-rDNA gene and shows Spirostomidae as the stem group of Heterotrichea; (4) however, the total-evidence tree suggests that the large Heterotrichea cluster is divided into nine lineages in which Peritromidae diverges at the base of the Heterotrichea tree.