ParameciumDB 2019: integrating genomic data across the genus for functional and evolutionary biology

Application of RNA interference and protein localization to investigate housekeeping and developmentally regulated genes in the emerging model protozoan Paramecium caudatum

Unicellular eukaryotes represent tremendous evolutionary diversity. However, the molecular mechanisms underlying this diversity remain largely unexplored, partly due to a limitation of genetic tools to only a few model species. Paramecium caudatum is a well-known unicellular eukaryote with an unexpectedly large germline genome, of which only two percent is retained in the somatic genome following sexual processes, revealing extensive DNA elimination. However, further progress in understanding the molecular mechanisms governing this process is hampered by a lack of suitable genetic tools. Here, we report the successful application of gene knockdown and protein localization methods to interrogate the function of both housekeeping and developmentally regulated genes in P. caudatum. Using these methods, we achieved the expected phenotypes upon RNAi by feeding, and determined the localization of these proteins by microinjection of fusion constructs containing fluorescent protein or antibody tags. Lastly, we used these methods to reveal that P. caudatum PiggyMac, a domesticated piggyBac transposase, is essential for sexual development, and is likely to be an active transposase directly involved in DNA cleavage. The application of these methods lays the groundwork for future studies of gene function in P. caudatum and can be used to answer important biological questions in the future.

Programmed chromosome fragmentation in ciliated protozoa: multiple means to chromosome ends

SUMMARYCiliated protozoa undergo large-scale developmental rearrangement of their somatic genomes when forming a new transcriptionally active macronucleus during conjugation. This process includes the fragmentation of chromosomes derived from the germline, coupled with the efficient healing of the broken ends by de novo telomere addition. Here, we review what is known of developmental chromosome fragmentation in ciliates that have been well-studied at the molecular level (Tetrahymena, Paramecium, Euplotes, Stylonychia, and Oxytricha). These organisms differ substantially in the fidelity and precision of their fragmentation systems, as well as in the presence or absence of well-defined sequence elements that direct excision, suggesting that chromosome fragmentation systems have evolved multiple times and/or have been significantly altered during ciliate evolution. We propose a two-stage model for the evolution of the current ciliate systems, with both stages involving repetitive or transposable elements in the genome. The ancestral form of chromosome fragmentation is proposed to have been derived from the ciliate small RNA/chromatin modification process that removes transposons and other repetitive elements from the macronuclear genome during development. The evolution of this ancestral system is suggested to have potentiated its replacement in some ciliate lineages by subsequent fragmentation systems derived from mobile genetic elements.

Potential role of N6-adenine DNA methylation in alternative splicing and endosymbiosis in Paramecium bursaria

N⁶-adenine DNA methylation (6mA), a rediscovered epigenetic mark in eukaryotic organisms, diversifies in abundance, distribution, and function across species, necessitating its study in more taxa. Paramecium bursaria is a typical model organism with endosymbiotic algae of the species Chlorella variabilis. This consortium therefore serves as a valuable system to investigate the functional role of 6mA in endosymbiosis, as well as the evolutionary importance of 6mA among eukaryotes. In this study, we report the first genome-wide, base pair-resolution map of 6mA in P. bursaria and identify its methyltransferase PbAMT1. Functionally, 6mA exhibits a bimodal distribution at the 5' end of RNA polymerase II-transcribed genes and possibly participates in transcription by facilitating alternative splicing. Evolutionarily, 6mA co-evolves with gene age and likely serves as a reverse mark of endosymbiosis-related genes. Our results offer new insights for the functional diversification of 6mA in eukaryotes as an important epigenetic mark.

iGDP: An integrated genome decontamination pipeline for wild ciliated microeukaryotes

Ciliates are a large group of ubiquitous and highly diverse single-celled eukaryotes that play an essential role in the functioning of microbial food webs. However, their genomic diversity is far from clear due to the need to develop cultivation methods for most species, so most research is based on wild organisms that almost invariably contain contaminants. Here we establish an integrated Genome Decontamination Pipeline (iGDP) that combines homology search, telomere reads-assisted and clustering approaches to filter contaminated ciliate genome assemblies from wild specimens. We benchmarked the performance of iGDP using genomic data from a contaminated ciliate culture and the results showed that iGDP could recall 91.9% of the target sequences with 96.9% precision. We also used a synthetic dataset to offer guidelines for the application of iGDP in the removal of various groups of contaminants. Compared with several popular metagenome binning tools, iGDP could show better performance. To further validate the effectiveness of iGDP on real-world data, we applied it to decontaminate genome assemblies of three wild ciliate specimens and obtained their genomes with high quality comparable to that of previously well-studied model ciliate genomes. It is anticipated that the newly generated genomes and the established iGDP method will be valuable community resources for detailed studies on ciliate biodiversity, phylogeny, ecology and evolution. The pipeline (https://github.com/GWang2022/iGDP) can be implemented automatically to reduce manual filtering and classification and may be further developed to apply to other microeukaryotes.

Paramecium bursaria as a Potential Tool for Evaluation of Microplastics Toxicity

Microplastics (MPs) are normally defined as small plastic wastes with a size of 1 μm to 5 mm in diameter. This tiny plastic debris is abundant in aquatic systems and poses a great threat to aquatic biota. To date, toxicological assessment of MPs is predominantly dependent on metazoan animals, although their applications are sometimes limited due to the high cost, narrow ecological niche, or ethical considerations. In this regard, unicellular eukaryotes (i.e., protozoa) that are ubiquitously present in nature represent a promising alternative for evaluating the toxicity of MPs. In this study, we selected Paramecium bursaria (P. bursaria) as a representative of protozoa and further investigated behavioral and molecular changes in MPs-exposed P. bursaria. Our results showed that following MPs uptake, P. bursaria exhibited various changes, including anomalies in swimming patterns, reduction in moving speed, impairment of avoidance behavior, elevation of oxidative stress, and potential disturbance of endosymbiosis. These elicited changes in P. bursaria in response to MPs exposure were pronounced and measurable. Overall, this study demonstrated that P. bursaria could serve as a promising alternative for the toxicological assessment of MPs and may be further applied to evaluate the toxicity of other environmental contaminants.

Calmodulin in Paramecium: Focus on Genomic Data

Calcium (Ca²⁺) is a universal second messenger that plays a key role in cellular signaling. However, Ca²⁺ signals are transduced with the help of Ca²⁺-binding proteins, which serve as sensors, transducers, and elicitors. Among the collection of these Ca²⁺-binding proteins, calmodulin (CaM) emerged as the prototypical model in eukaryotic cells. This is a small protein that binds four Ca²⁺ ions and whose functions are multiple, controlling many essential aspects of cell physiology. CaM is universally distributed in eukaryotes, from multicellular organisms, such as human and land plants, to unicellular microorganisms, such as yeasts and ciliates. Here, we review most of the information gathered on CaM in Paramecium, a group of ciliates. We condense the information here by mentioning that mature Paramecium CaM is a 148 amino acid-long protein codified by a single gene, as in other eukaryotic microorganisms. In these ciliates, the protein is notoriously localized and regulates cilia function and can stimulate the activity of some enzymes. When Paramecium CaM is mutated, cells show flawed locomotion and/or exocytosis. We further widen this and additional information in the text, focusing on genomic data.

The evolutionary conserved proteins CEP90, FOPNL, and OFD1 recruit centriolar distal appendage proteins to initiate their assembly

In metazoa, cilia assembly is a cellular process that starts with centriole to basal body maturation, migration to the cell surface, and docking to the plasma membrane. Basal body docking involves the interaction of both the distal end of the basal body and the transition fibers/distal appendages, with the plasma membrane. Mutations in numerous genes involved in basal body docking and transition zone assembly are associated with the most severe ciliopathies, highlighting the importance of these events in cilium biogenesis. In this context, the ciliate Paramecium has been widely used as a model system to study basal body and cilia assembly. However, despite the evolutionary conservation of cilia assembly events across phyla, whether the same molecular players are functionally conserved, is not fully known. Here, we demonstrated that CEP90, FOPNL, and OFD1 are evolutionary conserved proteins crucial for ciliogenesis. Using ultrastructure expansion microscopy, we unveiled that these proteins localize at the distal end of both centrioles/basal bodies in Paramecium and mammalian cells. Moreover, we found that these proteins are recruited early during centriole duplication on the external surface of the procentriole. Functional analysis performed both in Paramecium and mammalian cells demonstrate the requirement of these proteins for distal appendage assembly and basal body docking. Finally, we show that mammalian centrioles require another component, Moonraker (MNR), to recruit OFD1, FOPNL, and CEP90, which will then recruit the distal appendage proteins CEP83, CEP89, and CEP164. Altogether, we propose that this OFD1, FOPNL, and CEP90 functional module is required to determine in mammalian cells the future position of distal appendage proteins.

CEP90, FOPNL and OFD1 form an evolutionary conserved module which promotes the assembly of centriolar distal appendages. This study uses ultrastructure expansion microscopy to reveal the recruitment of this module on early-born procentrioles to in turn recruit centriolar distal appendage proteins, proposing that this dictates the future location of distal appendages.

GC content, but not nucleosome positioning, directly contributes to intron splicing efficiency in Paramecium

Eukaryotic genes are interrupted by introns that must be accurately spliced from mRNA precursors. With an average length of 25 nt, the more than 90,000 introns of Paramecium tetraurelia stand among the shortest introns reported in eukaryotes. The mechanisms specifying the correct recognition of these tiny introns remain poorly understood. Splicing can occur cotranscriptionally, and it has been proposed that chromatin structure might influence splice site recognition. To investigate the roles of nucleosome positioning in intron recognition, we determined the nucleosome occupancy along the P. tetraurelia genome. We show that P. tetraurelia displays a regular nucleosome array with a nucleosome repeat length of ∼151 bp, among the smallest periodicities reported. Our analysis has revealed that introns are frequently associated with inter-nucleosomal DNA, pointing to an evolutionary constraint favoring introns at the AT-rich nucleosome edge sequences. Using accurate splicing efficiency data from cells depleted for nonsense-mediated decay effectors, we show that introns located at the edge of nucleosomes display higher splicing efficiency than those at the center. However, multiple regression analysis indicates that the low GC content of introns, rather than nucleosome positioning, is associated with high splicing efficiency. Our data reveal a complex link between GC content, nucleosome positioning, and intron evolution in Paramecium.

Developmental timing of programmed DNA elimination in Paramecium tetraurelia recapitulates germline transposon evolutionary dynamics

With its nuclear dualism, the ciliate Paramecium constitutes a unique model to study how host genomes cope with transposable elements (TEs). P. tetraurelia harbors two germline micronuclei (MICs) and a polyploid somatic macronucleus (MAC) that develops from one MIC at each sexual cycle. Throughout evolution, the MIC genome has been continuously colonized by TEs and related sequences that are removed from the somatic genome during MAC development. Whereas TE elimination is generally imprecise, excision of approximately 45,000 TE-derived internal eliminated sequences (IESs) is precise, allowing for functional gene assembly. Programmed DNA elimination is concomitant with genome amplification. It is guided by noncoding RNAs and repressive chromatin marks. A subset of IESs is excised independently of this epigenetic control, raising the question of how IESs are targeted for elimination. To gain insight into the determinants of IES excision, we established the developmental timing of DNA elimination genome-wide by combining fluorescence-assisted nuclear sorting with high-throughput sequencing. Essentially all IESs are excised within only one endoreplication round (32C to 64C), whereas TEs are eliminated at a later stage. We show that DNA elimination proceeds independently of replication. We defined four IES classes according to excision timing. The earliest excised IESs tend to be independent of epigenetic factors, display strong sequence signals at their ends, and originate from the most ancient integration events. We conclude that old IESs have been optimized during evolution for early and accurate excision by acquiring stronger sequence determinants and escaping epigenetic control.

The Cilioprotist Cytoskeleton , a Model for Understanding How Cell Architecture and Pattern Are Specified: Recent Discoveries from Ciliates and Comparable Model Systems

The cytoskeletons of eukaryotic, cilioprotist microorganisms are complex, highly patterned, and diverse, reflecting the varied and elaborate swimming, feeding, reproductive, and sensory behaviors of the multitude of cilioprotist species that inhabit the aquatic environment. In the past 10-20 years, many new discoveries and technologies have helped to advance our understanding of how cytoskeletal organelles are assembled in many different eukaryotic model systems, in relation to the construction and modification of overall cellular architecture and function. Microtubule organizing centers, particularly basal bodies and centrioles, have continued to reveal their central roles in architectural engineering of the eukaryotic cell, including in the cilioprotists. This review calls attention to (1) published resources that illuminate what is known of the cilioprotist cytoskeleton; (2) recent studies on cilioprotists and other model organisms that raise specific questions regarding whether basal body- and centriole-associated nucleic acids, both DNA and RNA, should continue to be considered when seeking to employ cilioprotists as model systems for cytoskeletal research; and (3) new, mainly imaging, technologies that have already proven useful for, but also promise to enhance, future cytoskeletal research on cilioprotists.

Ciliary Ca2+ pumps regulate intraciliary Ca2+ from the action potential and may co-localize with ciliary voltage-gated Ca2+ channels

Calcium ions (Ca2+) entering cilia through the ciliary voltage-gated calcium channels (CaV) during the action potential causes reversal of the ciliary power stroke and backward swimming in Paramecium tetraurelia. How calcium is returned to the resting level is not yet clear. Our focus is on calcium pumps as a possible mechanism. There are 23 P. tetraurelia genes for calcium pumps that are members of the family of plasma membrane Ca2+ ATPases (PMCAs). They have domains homologous to those found in mammalian PMCAs. Of the 13 pump proteins previously identified in cilia, ptPMCA2a and ptPMCA2b are most abundant in the cilia. We used RNAi to examine which PMCA might be involved in regulating intraciliary Ca2+ after the action potential. RNAi for only ptPMCA2a and ptPMCA2b causes cells to significantly prolong their backward swimming, which indicates that Ca2+ extrusion in the cilia is impaired when these PMCAs are depleted. We used immunoprecipitations (IP) to find that ptPMCA2a and ptPMCA2b are co-immunoprecipitated with the CaV channel α1 subunits that are found only in the cilia. We used iodixanol (OptiPrep) density gradients to show that ptPMCA2a and ptPMCA2b and CaV1c are found in the same density fractions. These results suggest that ptPMCA2a and ptPMCA2b are located in the proximity of ciliary CaV channels.

Characterization of the RNA-interference pathway as a tool for reverse genetic analysis in the nascent phototrophic endosymbiosis, Paramecium bursaria

Endosymbiosis was fundamental for the evolution of eukaryotic complexity. Endosymbiotic interactions can be dissected through forward- and reverse-genetic experiments, such as RNA-interference (RNAi). However, distinguishing small (s)RNA pathways in a eukaryote-eukaryote endosymbiotic interaction is challenging. Here, we investigate the repertoire of RNAi pathway protein-encoding genes in the model nascent endosymbiotic system, Paramecium bursaria-Chlorella spp. Using comparative genomics and transcriptomics supported by phylogenetics, we identify essential proteome components of the small interfering (si)RNA, scan (scn)RNA and internal eliminated sequence (ies)RNA pathways. Our analyses reveal that copies of these components have been retained throughout successive whole genome duplication (WGD) events in the Paramecium clade. We validate feeding-induced siRNA-based RNAi in P. bursaria via knock-down of the splicing factor, u2af1, which we show to be crucial to host growth. Finally, using simultaneous knock-down 'paradox' controls to rescue the effect of u2af1 knock-down, we demonstrate that feeding-induced RNAi in P. bursaria is dependent upon a core pathway of host-encoded Dcr1, Piwi and Pds1 components. Our experiments confirm the presence of a functional, host-derived RNAi pathway in P. bursaria that generates 23-nt siRNA, validating the use of the P. bursaria-Chlorella spp. system to investigate the genetic basis of a nascent endosymbiosis.

Morphology, ultrastructure, genomics, and phylogeny of Euplotes vanleeuwenhoeki sp. nov. and its ultra-reduced endosymbiont "Candidatus Pinguicoccus supinus" sp. nov

Taxonomy is the science of defining and naming groups of biological organisms based on shared characteristics and, more recently, on evolutionary relationships. With the birth of novel genomics/bioinformatics techniques and the increasing interest in microbiome studies, a further advance of taxonomic discipline appears not only possible but highly desirable. The present work proposes a new approach to modern taxonomy, consisting in the inclusion of novel descriptors in the organism characterization: (1) the presence of associated microorganisms (e.g.: symbionts, microbiome), (2) the mitochondrial genome of the host, (3) the symbiont genome. This approach aims to provide a deeper comprehension of the evolutionary/ecological dimensions of organisms since their very first description. Particularly interesting, are those complexes formed by the host plus associated microorganisms, that in the present study we refer to as "holobionts". We illustrate this approach through the description of the ciliate Euplotes vanleeuwenhoeki sp. nov. and its bacterial endosymbiont "Candidatus Pinguicoccus supinus" gen. nov., sp. nov. The endosymbiont possesses an extremely reduced genome (~ 163 kbp); intriguingly, this suggests a high integration between host and symbiont.

The Paramecium histone chaperone Spt16-1 is required for Pgm endonuclease function in programmed genome rearrangements

In Paramecium tetraurelia, a large proportion of the germline genome is reproducibly removed from the somatic genome after sexual events via a process involving small (s)RNA-directed heterochromatin formation and DNA excision and repair. How germline limited DNA sequences are specifically recognized in the context of chromatin remains elusive. Here, we use a reverse genetics approach to identify factors involved in programmed genome rearrangements. We have identified a P. tetraurelia homolog of the highly conserved histone chaperone Spt16 subunit of the FACT complex, Spt16-1, and show its expression is developmentally regulated. A functional GFP-Spt16-1 fusion protein localized exclusively in the nuclei where genome rearrangements take place. Gene silencing of Spt16-1 showed it is required for the elimination of all germline-limited sequences, for the survival of sexual progeny, and for the accumulation of internal eliminated sequence (ies)RNAs, an sRNA population produced when elimination occurs. Normal accumulation of 25 nt scanRNAs and deposition of silent histone marks H3K9me3 and H3K27me3 indicated that Spt16-1 does not regulate the scanRNA-directed heterochromatin pathway involved in the early steps of DNA elimination. We further show that Spt16-1 is required for the correct nuclear localization of the PiggyMac (Pgm) endonuclease, which generates the DNA double-strand breaks required for DNA elimination. Thus, Spt16-1 is essential for Pgm function during programmed genome rearrangements. We propose a model in which Spt16-1 mediates interactions between the excision machinery and chromatin, facilitating endonuclease access to DNA cleavage sites during genome rearrangements.

The genome is generally similar in all the cells of an organism. However, in the ciliate Paramecium tetraurelia, massive and reproducible programmed DNA elimination leads to a highly streamlined somatic genome. In eukaryotes, DNA is packaged into nucleosomes, which ensure genome integrity but act as a barrier to enzymes acting on DNA. How the endonuclease PiggyMac gains access to the genome to initiate DNA elimination remains elusive. Here, we identified four P. tetraurelia genes encoding homologs of the conserved histone chaperone Spt16, which can modulate access to DNA by promoting nucleosome assembly and disassembly. We demonstrated that the most divergent gene, SPT16-1, has a highly specialized expression pattern, similar to that of PiggyMac, and a specific role in programmed DNA elimination. We show that the Spt16-1 protein, like PiggyMac, is exclusively localized in the differentiating somatic nucleus, and is also required for the dramatic elimination of germline-limited sequences. We further show that Spt16-1 directs the correct nuclear localization of the PiggyMac endonuclease. Thus, Spt16-1 is essential for PiggyMac function during programmed DNA elimination. We propose that Spt16-1 mediates the interaction between PiggyMac and chromatin or DNA, facilitating endonuclease access to DNA cleavage sites.

The Obligate Symbiont "Candidatus Megaira polyxenophila" Has Variable Effects on the Growth of Different Host Species

"Candidatus Megaira polyxenophila" is a recently described member of Rickettsiaceae which comprises exclusively obligate intracellular bacteria. Interestingly, these bacteria can be found in a huge diversity of eukaryotic hosts (protist, green algae, metazoa) living in marine, brackish or freshwater habitats. Screening of amplicon datasets revealed a high frequency of these bacteria especially in freshwater environments, most likely associated to eukaryotic hosts. The relationship of "Ca. Megaira polyxenophila" with their hosts and their impact on host fitness have not been studied so far. Even less is known regarding the responses of these intracellular bacteria to potential stressors. In this study, we used two phylogenetically close species of the freshwater ciliate Paramecium, Paramecium primaurelia and Paramecium pentaurelia (Ciliophora, Oligohymenophorea) naturally infected by "Ca. Megaira polyxenophila". In order to analyze the effect of the symbiont on the fitness of these two species, we compared the growth performance of both infected and aposymbiotic paramecia at different salinity levels in the range of freshwater and oligohaline brackish water i.e., at 0, 2, and 4.5 ppt. For the elimination of "Ca. Megaira polyxenophila" we established an antibiotic treatment to obtain symbiont-free lines and confirmed its success by fluorescence in situ hybridization (FISH). The population and infection dynamics during the growth experiment were observed by cell density counts and FISH. Paramecia fitness was compared applying generalized additive mixed models. Surprisingly, both infected Paramecium species showed higher densities under all salinity concentrations. The tested salinity concentrations did not significantly affect the growth of any of the two species directly, but we observed the loss of the endosymbiont after prolonged exposure to higher salinity levels. This experimental data might explain the higher frequency of "Ca. M. polyxenophila" in freshwater habitats as observed from amplicon data.

Paramecium Diversity and a New Member of the Paramecium aurelia Species Complex Described from Mexico

Paramecium (Ciliophora) is an ideal model organism to study the biogeography of protists. However, many regions of the world, such as Central America, are still neglected in understanding Paramecium diversity. We combined morphological and molecular approaches to identify paramecia isolated from more than 130 samples collected from different waterbodies in several states of Mexico. We found representatives of six Paramecium morphospecies, including the rare species Paramecium jenningsi, and Paramecium putrinum, which is the first report of this species in tropical regions. We also retrieved five species of the Paramecium aurelia complex, and describe one new member of the complex, Paramecium quindecaurelia n. sp., which appears to be a sister species of Paramecium biaurelia. We discuss criteria currently applied for differentiating between sibling species in Paramecium. Additionally, we detected diverse bacterial symbionts in some of the collected ciliates.

Functional diversification of Paramecium Ku80 paralogs safeguards genome integrity during precise programmed DNA elimination

Gene duplication and diversification drive the emergence of novel functions during evolution. Because of whole genome duplications, ciliates from the Paramecium aurelia group constitute a remarkable system to study the evolutionary fate of duplicated genes. Paramecium species harbor two types of nuclei: a germline micronucleus (MIC) and a somatic macronucleus (MAC) that forms from the MIC at each sexual cycle. During MAC development, ~45,000 germline Internal Eliminated Sequences (IES) are excised precisely from the genome through a 'cut-and-close' mechanism. Here, we have studied the P. tetraurelia paralogs of KU80, which encode a key DNA double-strand break repair factor involved in non-homologous end joining. The three KU80 genes have different transcription patterns, KU80a and KU80b being constitutively expressed, while KU80c is specifically induced during MAC development. Immunofluorescence microscopy and high-throughput DNA sequencing revealed that Ku80c stably anchors the PiggyMac (Pgm) endonuclease in the developing MAC and is essential for IES excision genome-wide, providing a molecular explanation for the previously reported Ku-dependent licensing of DNA cleavage at IES ends. Expressing Ku80a under KU80c transcription signals failed to complement a depletion of endogenous Ku80c, indicating that the two paralogous proteins have distinct properties. Domain-swap experiments identified the α/β domain of Ku80c as the major determinant for its specialized function, while its C-terminal part is required for excision of only a small subset of IESs located in IES-dense regions. We conclude that Ku80c has acquired the ability to license Pgm-dependent DNA cleavage, securing precise DNA elimination during programmed rearrangements. The present study thus provides novel evidence for functional diversification of genes issued from a whole-genome duplication.

Epidemiology of Nucleus-Dwelling Holospora: Infection, Transmission, Adaptation, and Interaction with Paramecium

The chapter describes the exceptional symbiotic associations formed between the ciliate Paramecium and Holospora, highly infectious bacteria residing in the host nuclei. Holospora and Holospora-like bacteria (Alphaproteobacteria) are characterized by their ability for vertical and horizontal transmission in host populations, a complex biphasic life cycle, and pronounced preference for host species and colonized cell compartment. These bacteria are obligate intracellular parasites; thus, their metabolic repertoire is dramatically reduced. Nevertheless, they perform complex interactions with the host ciliate. We review ongoing efforts to unravel the molecular adaptations of these bacteria to their unusual lifestyle and the host's employment in the symbiosis. Furthermore, we summarize current knowledge on the genetic and genomic background of Paramecium-Holospora symbiosis and provide insights into the ecological and evolutionary consequences of this interaction. The diversity and occurrence of symbioses between ciliates and Holospora-like bacteria in nature is discussed in connection with transmission modes of symbionts, host specificity and compatibility of the partners. We aim to summarize 50 years of research devoted to these symbiotic systems and conclude trying to predict some perspectives for further studies.