While molecular analyses have provided insight into the phylogeny of ciliates, the few studies assessing intraspecific variation have largely relied on just a single locus (e.g. nuclear small subunit rDNA (nSSU-rDNA) or mitochondrial cytochrome oxidase I). In this study, we characterize the diversity of several nuclear protein-coding genes plus both nSSU-rDNA and mitochondrial small subunit rDNA (mtSSU-rDNA) of five isolates of the ciliate morphospecies Chilodonella uncinata. Although these isolates have nearly identical nSSU-rDNA sequences, they differ by up to 8.0% in mtSSU-rDNA. Comparative analyses of all loci, including β-tubulin paralogs, indicate a lack of recombination between strains, demonstrating that the morphospecies C. uncinata consists of multiple cryptic species. Further, there is considerable variation in substitution rates among loci as some protein-coding domains are nearly identical between isolates while others differ by up to 13.2% at the amino acid level. Combining insights on macronuclear variation among isolates, the focus of this study, with published data from the micronucleus of two of these isolates indicates that C. uncinata lineages are able to maintain both highly divergent and highly conserved genes within a rapidly evolving germline genome.

The multiple species concepts currently in use by the scientific community (e.g. Morphological, Biological, Phylogenetic) are united in that they all aim to capture the process of divergence between populations. For example, the Biological Species Concept (BSC) defines a species as a natural group of organisms that is reproductively isolated from other such groups. Here we synthesize nearly a century of research on the ciliate genus Paramecium that highlights the shortcomings of our prevailing notions on the nature of species. In this lineage, there is discordance between morphology, mating behavior, and genetics, features assumed to be correlated, at least after sufficient time has passed, under all species concepts. Intriguingly, epigenetic phenomena are well documented in ciliates where they influence features such as germline/soma differentiation and mating type determination. Consequently, we hypothesize that divergence within ciliate populations is due to a dynamic interaction between genetic and epigenetic factors. The growing list of examples of epigenetic phenomena that potentially impact speciation (i.e. by influencing the dynamics of sex chromosomes, fate of hybrids, zygotic drive and genomic conflicts) suggests that interactions between genetics and epigenetics may also drive divergence in other eukaryotic lineages.