Estimation of divergence times in litostomatean ciliates (Ciliophora: Intramacronucleata), using Bayesian relaxed clock and 18S rRNA gene.

The class Litostomatea comprises a diverse assemblage of free-living and endosymbiotic ciliates. To understand diversification dynamic of litostomateans, divergence times of their main groups were estimated with the Bayesian molecular dating, a technique allowing relaxation of molecular clock and incorporation of flexible calibration points. The class Litostomatea very likely emerged during the Cryogenian around 680 Mya. The origin of the subclass Rhynchostomatia is dated to about 415 Mya, while that of the subclass Haptoria to about 654 Mya. The order Pleurostomatida, emerging about 556 Mya, was recognized as the oldest group within the subclass Haptoria. The order Spathidiida appeared in the Paleozoic about 442 Mya. The three remaining haptorian orders evolved in the Paleozoic/Mesozoic periods: Didiniida about 419 Mya, Lacrymariida about 269 Mya, and Haptorida about 194 Mya. The subclass Trichostomatia originated from a spathidiid ancestor in the Mesozoic about 260 Mya. A further goal of this study was to investigate the impact of various settings on posterior divergence time estimates. The root placement and tree topology as well as the priors of the rate-drift model, birth-death process and nucleotide substitution rate, had no significant effect on calculation of posterior divergence time estimates. However, removal of calibration points could significantly change time estimates at some nodes.