Massive colonization of protein-coding exons by selfish genetic elements in Paramecium germline genomes.

Ciliates are unicellular eukaryotes with both a germline genome and a somatic genome in the same cytoplasm. The somatic macronucleus (MAC), responsible for gene expression, is not sexually transmitted but develops from a copy of the germline micronucleus (MIC) at each sexual generation. In the MIC genome of Paramecium tetraurelia, genes are interrupted by tens of thousands of unique intervening sequences called internal eliminated sequences (IESs), which have to be precisely excised during the development of the new MAC to restore functional genes. To understand the evolutionary origin of this peculiar genomic architecture, we sequenced the MIC genomes of 9 Paramecium species (from approximately 100 Mb in Paramecium aurelia species to >1.5 Gb in Paramecium caudatum). We detected several waves of IES gains, both in ancestral and in more recent lineages. While the vast majority of IESs are single copy in present-day genomes, we identified several families of mobile IESs, including nonautonomous elements acquired via horizontal transfer, which generated tens to thousands of new copies. These observations provide the first direct evidence that transposable elements can account for the massive proliferation of IESs in Paramecium. The comparison of IESs of different evolutionary ages indicates that, over time, IESs shorten and diverge rapidly in sequence while they acquire features that allow them to be more efficiently excised. We nevertheless identified rare cases of IESs that are under strong purifying selection across the aurelia clade. The cases examined contain or overlap cellular genes that are inactivated by excision during development, suggesting conserved regulatory mechanisms. Similar to the evolution of introns in eukaryotes, the evolution of Paramecium IESs highlights the major role played by selfish genetic elements in shaping the complexity of genome architecture and gene expression.

"Candidatus Intestinibacterium parameciiphilum"-member of the "Candidatus Paracaedibacteraceae" family (Alphaproteobacteria, Holosporales) inhabiting the ciliated protist Paramecium.

Protists frequently host diverse bacterial symbionts, in particular those affiliated with the order Holosporales (Alphaproteobacteria). All characterised members of this bacterial lineage have been retrieved in obligate association with a wide range of eukaryotes, especially multiple protist lineages (e.g. amoebozoans, ciliates, cercozoans, euglenids, and nucleariids), as well as some metazoans (especially arthropods and related ecdysozoans). While the genus Paramecium and other ciliates have been deeply investigated for the presence of symbionts, known members of the family "Candidatus Paracaedibacteraceae" (Holosporales) are currently underrepresented in such hosts. Herein, we report the description of "Candidatus Intestinibacterium parameciiphilum" within the family "Candidatus Paracaedibacteraceae", inhabiting the cytoplasm of Paramecium biaurelia. This novel bacterium is almost twice as big as its relative "Candidatus Intestinibacterium nucleariae" from the opisthokont Nuclearia and does not present a surrounding halo. Based on phylogenetic analyses of 16S rRNA gene sequences, we identified six further potential species-level lineages within the genus. Based on the provenance of the respective samples, we investigated the environmental distribution of the representatives of "Candidatus Intestinibacterium" species. Obtained results are consistent with an obligate endosymbiotic lifestyle, with protists, in particular freshwater ones, as hosts. Thus, available data suggest that association with freshwater protists could be the ancestral condition for the members of the "Candidatus Intestinibacterium" genus.

Population Genomics of Paramecium Species.

Population-genomic analyses are essential to understanding factors shaping genomic variation and lineage-specific sequence constraints. The dearth of such analyses for unicellular eukaryotes prompted us to assess genomic variation in Paramecium, one of the most well-studied ciliate genera. The Paramecium aurelia complex consists of ∼15 morphologically indistinguishable species that diverged subsequent to two rounds of whole-genome duplications (WGDs, as long as 320 MYA) and possess extremely streamlined genomes. We examine patterns of both nuclear and mitochondrial polymorphism, by sequencing whole genomes of 10-13 worldwide isolates of each of three species belonging to the P. aurelia complex: P. tetraurelia, P. biaurelia, P. sexaurelia, as well as two outgroup species that do not share the WGDs: P. caudatum and P. multimicronucleatum. An apparent absence of global geographic population structure suggests continuous or recent dispersal of Paramecium over long distances. Intergenic regions are highly constrained relative to coding sequences, especially in P. caudatum and P. multimicronucleatum that have shorter intergenic distances. Sequence diversity and divergence are reduced up to ∼100-150 bp both upstream and downstream of genes, suggesting strong constraints imposed by the presence of densely packed regulatory modules. In addition, comparison of sequence variation at non-synonymous and synonymous sites suggests similar recent selective pressures on paralogs within and orthologs across the deeply diverging species. This study presents the first genome-wide population-genomic analysis in ciliates and provides a valuable resource for future studies in evolutionary and functional genetics in Paramecium.

Evolutionary rescue and local adaptation under different rates of temperature increase: a combined analysis of changes in phenotype expression and genotype frequency in Paramecium microcosms.

Evolutionary rescue (ER) occurs when populations, which have declined due to rapid environmental change, recover through genetic adaptation. The success of this process and the evolutionary trajectory of the population strongly depend on the rate of environmental change. Here we investigated how different rates of temperature increase (from 23 to 32 °C) affect population persistence and evolutionary change in experimental microcosms of the protozoan Paramecium caudatum. Consistent with theory on ER, we found that those populations experiencing the slowest rate of temperature increase were the least likely to become extinct and tended to be the best adapted to the new temperature environment. All high-temperature populations were more tolerant to severe heat stress (35, 37 °C), indicating a common mechanism of heat protection. High-temperature populations also had superior growth rates at optimum temperatures, leading to the absence of a pattern of local adaptation to control (23 °C) and high-temperature (32 °C) environments. However, high-temperature populations had reduced growth at low temperatures (5-9 °C), causing a shift in the temperature niche. In part, the observed evolutionary change can be explained by selection from standing variation. Using mitochondrial markers, we found complete divergence between control and high-temperature populations in the frequencies of six initial founder genotypes. Our results confirm basic predictions of ER and illustrate how adaptation to an extreme local environment can produce positive as well as negative correlated responses to selection over the entire range of the ecological niche.

Disentangling the Taxonomy of Rickettsiales and Description of Two Novel Symbionts ("Candidatus Bealeia paramacronuclearis" and "Candidatus Fokinia cryptica") Sharing the Cytoplasm of the Ciliate Protist Paramecium biaurelia.

In the past 10 years, the number of endosymbionts described within the bacterial order Rickettsiales has constantly grown. Since 2006, 18 novel Rickettsiales genera inhabiting protists, such as ciliates and amoebae, have been described. In this work, we characterize two novel bacterial endosymbionts from Paramecium collected near Bloomington, IN. Both endosymbiotic species inhabit the cytoplasm of the same host. The Gram-negative bacterium "Candidatus Bealeia paramacronuclearis" occurs in clumps and is frequently associated with the host macronucleus. With its electron-dense cytoplasm and a distinct halo surrounding the cell, it is easily distinguishable from the second smaller symbiont, "Candidatus Fokinia cryptica," whose cytoplasm is electron lucid, lacks a halo, and is always surrounded by a symbiontophorous vacuole. For molecular characterization, the small-subunit rRNA genes were sequenced and used for taxonomic assignment as well as the design of species-specific oligonucleotide probes. Phylogenetic analyses revealed that "Candidatus Bealeia paramacronuclearis" clusters with the so-called "basal" Rickettsiales, and "Candidatus Fokinia cryptica" belongs to "Candidatus Midichloriaceae." We obtained tree topologies showing a separation of Rickettsiales into at least two groups: one represented by the families Rickettsiaceae, Anaplasmataceae, and "Candidatus Midichloriaceae" (RAM clade), and the other represented by "basal Rickettsiales," including "Candidatus Bealeia paramacronuclearis." Therefore, and in accordance with recent publications, we propose to limit the order Rickettsiales to the RAM clade and to raise "basal Rickettsiales" to an independent order, Holosporales ord. nov., inside Alphaproteobacteria, which presently includes four family-level clades. Additionally, we define the family "Candidatus Hepatincolaceae" and redefine the family Holosporaceae IMPORTANCE: In this paper, we provide the characterization of two novel bacterial symbionts inhabiting the same Paramecium host (Ciliophora, Alveolata). Both symbionts belong to "traditional" Rickettsiales, one representing a new species of the genus "Candidatus Fokinia" ("Candidatus Midichloriaceae"), and the other representing a new genus of a "basal" Rickettsiales According to newly characterized sequences and to a critical revision of recent literature, we propose a taxonomic reorganization of "traditional" Rickettsiales that we split into two orders: Rickettsiales sensu stricto and Holosporales ord. nov. This work represents a critical revision, including new records of a group of symbionts frequently occurring in protists and whose biodiversity is still largely underestimated.

From genes to toxins: Profiling Prymnesium parvum during a riverine harmful algal bloom.

Blooms of Prymnesium parvum, a unicellular alga globally distributed in marine and brackish environments, frequently result in massive fish kills due to the production of toxins called prymnesins by this haptophyte. In August 2022, a harmful algal bloom (HAB) of this species occurred in the lower Oder River (Poland and Germany), which caused mass mortalities of fish and other organisms. This HAB was linked to low discharge of the Oder and mining activities that caused a significant increase in salinity. In this context, we report on the molecular detection and screening of this haptophyte and its toxins in environmental samples and clonal cultures derived thereof. Both conventional PCR and droplet digital PCR assays reliably detected P. parvum in environmental samples. eDNA metabarcoding using the V4 region of the 18S rRNA gene revealed a single Prymnesium sequence variant, but failed to identify it to species level. Four clonal cultures established from environmental samples were unambiguously identified as P. parvum by molecular phylogenetics (near full-length 18S rRNA gene) and light microscopy. Phylogenetic analysis (ITS1-5.8S-ITS2 marker region) placed the cultured phylotype within a clade containing other P. parvum strains known to produce B-type prymnesins. Toxin-screening of the cultures using liquid chromatography-electrospray ionization - time of flight mass spectrometry identified B-type prymnesins, which were also detected in extracts of filter residues from water samples of the Oder collected during the HAB. Overall, our investigation provides a detailed characterization of P. parvum, including their prymnesins, during this HAB in the Oder River, contributing valuable insights into this ecological disaster. In addition, the droplet digital PCR assay established here will be useful for future monitoring of low levels of P. parvum on the Oder River or any other salt-impacted and brackish water bodies.

Convergent evolution of heat-inducibility during subfunctionalization of the Hsp70 gene family.

BACKGROUND: Heat-shock proteins of the 70 kDa family (Hsp70s) are essential chaperones required for key cellular functions. In eukaryotes, four subfamilies can be distinguished according to their function and localisation in different cellular compartments: cytosol, endoplasmic reticulum, mitochondria and chloroplasts. Generally, multiple cytosol-type Hsp70s can be found in metazoans that show either constitutive expression and/or stress-inducibility, arguing for the evolution of different tasks and functions. Information about the hsp70 copy number and diversity in microbial eukaryotes is, however, scarce, and detailed knowledge about the differential gene expression in most protists is lacking. Therefore, we have characterised the Hsp70 gene family of Paramecium caudatum to gain insight into the evolution and differential heat stress response of the distinct family members in protists and to investigate the diversification of eukaryotic hsp70s focusing on the evolution of heat-inducibility. RESULTS: Eleven putative hsp70 genes could be detected in P. caudatum comprising homologs of three major Hsp70-subfamilies. Phylogenetic analyses revealed five evolutionarily distinct Hsp70-groups, each with a closer relationship to orthologous sequences of Paramecium tetraurelia than to another P. caudatum Hsp70-group. These highly diverse, paralogous groups resulted from duplications preceding Paramecium speciation, underwent divergent evolution and were subject to purifying selection. Heat-shock treatments were performed to test for differential expression patterns among the five Hsp70-groups as well as for a functional conservation within Paramecium. These treatments induced exceptionally high mRNA up-regulations in one cytosolic group with a low basal expression, indicative for the major heat inducible hsp70s. All other groups showed comparatively high basal expression levels and moderate heat-inducibility, signifying constitutively expressed genes. Comparative EST analyses for P. tetraurelia hsp70s unveiled a corresponding expression pattern, which supports a functionally conserved evolution of the Hsp70 gene family in Paramecium. CONCLUSIONS: Our analyses suggest an independent evolution of the heat-inducible cytosol-type hsp70s in Paramecium and in its close relative Tetrahymena, as well as within higher eukaryotes. This result indicates convergent evolution during hsp70 subfunctionalization and implies that heat-inducibility evolved several times during the course of eukaryotic evolution.

Predicting evolution in experimental range expansions of an aquatic model system.

Predicting range expansion dynamics is an important goal of both fundamental and applied research in conservation and global change biology. However, this is challenging if ecological and evolutionary processes occur on the same time scale. Using the freshwater ciliate Paramecium caudatum, we combined experimental evolution and mathematical modeling to assess the predictability of evolutionary change during range expansions. In the experiment, we followed ecological dynamics and trait evolution in independently replicated microcosm populations in range core and front treatments, where episodes of natural dispersal alternated with periods of population growth. These eco-evolutionary conditions were recreated in a predictive mathematical model, parametrized with dispersal and growth data of the 20 founder strains in the experiment. We found that short-term evolution was driven by selection for increased dispersal in the front treatment and general selection for higher growth rates in all treatments. There was a good quantitative match between predicted and observed trait changes. Phenotypic divergence was further mirrored by genetic divergence between range core and front treatments. In each treatment, we found the repeated fixation of the same cytochrome c oxidase I (COI) marker genotype, carried by strains that also were the most likely winners in our model. Long-term evolution in the experimental range front lines resulted in the emergence of a dispersal syndrome, namely a competition-colonization trade-off. Altogether, both model and experiment highlight the potential importance of dispersal evolution as a driver of range expansions. Thus, evolution at range fronts may follow predictable trajectories, at least for simple scenarios, and predicting these dynamics may be possible from knowledge of few key parameters.

Among-Strain Variation in Resistance of Paramecium caudatum to the Endonuclear Parasite Holospora undulata: Geographic and Lineage-Specific Patterns.

Resistance is a key determinant in interactions between hosts and their parasites. Understanding the amount and distribution of variation in this trait between strains can provide insights into (co)evolutionary processes and their potential to shape patterns of diversity in natural populations. Using controlled inoculation in experimental mass cultures, we investigated the quantitative variation in resistance to the bacterial parasite Holospora undulata across a worldwide collection of strains of its ciliate host Paramecium caudatum. We combined the observed variation with available information on the phylogeny and biogeography of the strains. We found substantial variation in resistance among strains, with upper-bound values of broad-sense heritability >0.5 (intraclass correlation coefficients). Strain estimates of resistance were repeatable between laboratories and ranged from total resistance to near-complete susceptibility. Early (1 week post inoculation) measurements provided higher estimates of resistance heritability than did later measurements (2-3 weeks), possibly due to diverging epidemiological dynamics in replicate cultures of the same strains. Genetic distance (based on a neutral marker) was positively correlated with the difference in resistance phenotype between strains (r = 0.45), essentially reflecting differences between highly divergent clades (haplogroups) within the host species. Haplogroup A strains, mostly European, were less resistant to the parasite (49% infection prevalence) than non-European haplogroup B strains (28%). At a smaller geographical scale (within Europe), strains that are geographically closer to the parasite origin (Southern Germany) were more susceptible to infection than those from further away. These patterns are consistent with a picture of local parasite adaptation. Our study demonstrates ample natural variation in resistance on which selection can act and hints at symbiont adaptation producing signatures in geographic and lineage-specific patterns of resistance in this model system.

"Candidatus Trichorickettsia mobilis", a Rickettsiales bacterium, can be transiently transferred from the unicellular eukaryote Paramecium to the planarian Dugesia japonica.

Most of the microorganisms responsible for vector-borne diseases (VBD) have hematophagous arthropods as vector/reservoir. Recently, many new species of microorganisms phylogenetically related to agents of VBD were found in a variety of aquatic eukaryotic hosts; in particular, numerous new bacterial species related to the genus Rickettsia (Alphaproteobacteria, Rickettsiales) were discovered in protist ciliates and other unicellular eukaryotes. Although their pathogenicity for humans and terrestrial animals is not known, several indirect indications exist that these bacteria might act as etiological agents of possible VBD of aquatic organisms, with protists as vectors. In the present study, a novel strain of the Rickettsia-Like Organism (RLO) endosymbiont "Candidatus (Ca.) Trichorickettsia mobilis" was identified in the macronucleus of the ciliate Paramecium multimicronucleatum. We performed transfection experiments of this RLO to planarians (Dugesia japonica) per os. Indeed, the latter is a widely used model system for studying bacteria pathogenic to humans and other Metazoa. In transfection experiments, homogenized paramecia were added to food of antibiotic-treated planarians. Treated and non-treated (i.e. control) planarians were investigated at day 1, 3, and 7 after feeding for endosymbiont presence by means of PCR and ultrastructural analyses. Obtained results were fully concordant and suggest that this RLO endosymbiont can be transiently transferred from ciliates to metazoans, being detected up to day 7 in treated planarians' enterocytes. Our findings might offer insights into the potential role of ciliates or other protists as putative vectors for diseases caused by Rickettsiales or other RLOs and occurring in fish farms or in the wild.

The Genome Sequence of "Candidatus Fokinia solitaria": Insights on Reductive Evolution in Rickettsiales.

"Candidatus Fokinia solitaria" is an obligate intracellular endosymbiont of a unicellular eukaryote, a ciliate of the genus Paramecium. Here, we present the genome sequence of this bacterium and subsequent analysis. Phylogenomic analysis confirmed the previously reported positioning of the symbiont within the "Candidatus Midichloriaceae" family (order Rickettsiales), as well as its high sequence divergence from other members of the family, indicative of fast sequence evolution. Consistently with this high evolutionary rate, a comparative genomic analysis revealed that the genome of this symbiont is the smallest of the Rickettsiales to date. The reduced genome does not present flagellar genes, nor the pathway for the biosynthesis of lipopolysaccharides (present in all the other so far sequenced members of the family "Candidatus Midichloriaceae") or genes for the Krebs cycle (present, although not always complete, in Rickettsiales). These results indicate an evolutionary trend toward a stronger dependence on the host, in comparison with other members of the family. Two alternative scenarios are compatible with our results; "Candidatus Fokinia solitaria" could be either a recently evolved, vertically transmitted mutualist, or a parasite with a high host-specificity.

"Candidatus Fokinia solitaria", a Novel "Stand-Alone" Symbiotic Lineage of Midichloriaceae (Rickettsiales).

Recently, the family Midichloriaceae has been described within the bacterial order Rickettsiales. It includes a variety of bacterial endosymbionts detected in different metazoan host species belonging to Placozoa, Cnidaria, Arthropoda and Vertebrata. Representatives of Midichloriaceae are also considered possible etiological agents of certain animal diseases. Midichloriaceae have been found also in protists like ciliates and amoebae. The present work describes a new bacterial endosymbiont, "Candidatus Fokinia solitaria", retrieved from three different strains of a novel Paramecium species isolated from a wastewater treatment plant in Rio de Janeiro (Brazil). Symbionts were characterized through the full-cycle rRNA approach: SSU rRNA gene sequencing and fluorescence in situ hybridization (FISH) with three species-specific oligonucleotide probes. In electron micrographs, the tiny rod-shaped endosymbionts (1.2 x 0.25-0.35 µm in size) were not surrounded by a symbiontophorous vacuole and were located in the peripheral host cytoplasm, stratified in the host cortex in between the trichocysts or just below them. Frequently, they occurred inside autolysosomes. Phylogenetic analyses of Midichloriaceae apparently show different evolutionary pathways within the family. Some genera, such as "Ca. Midichloria" and "Ca. Lariskella", have been retrieved frequently and independently in different hosts and environmental surveys. On the contrary, others, such as Lyticum, "Ca. Anadelfobacter", "Ca. Defluviella" and the presently described "Ca. Fokinia solitaria", have been found only occasionally and associated to specific host species. These last are the only representatives in their own branches thus far. Present data do not allow to infer whether these genera, which we named "stand-alone lineages", are an indication of poorly sampled organisms, thus underrepresented in GenBank, or represent fast evolving, highly adapted evolutionary lineages.

Rediscovering the genus Lyticum, multiflagellated symbionts of the order Rickettsiales.

Among the bacterial symbionts harbored by the model organism Paramecium, many still lack a recent investigation that includes a molecular characterization. The genus Lyticum consists of two species of large-sized bacteria displaying numerous flagella, despite their inability to move inside their hosts' cytoplasm. We present a multidisciplinary redescription of both species, using the deposited type strains as well as newly collected material. On the basis of 16S rRNA gene sequences, we assigned Lyticum to the order Rickettsiales, that is intensely studied because of its pathogenic representatives and its position as the extant group most closely related to the mitochondrial ancestor. We provide conclusive proofs that at least some Rickettsiales possess actual flagella, a feature that has been recently predicted from genomic data but never confirmed. We give support to the hypothesis that the mitochondrial ancestor could have been flagellated, and provide the basis for further studies on these ciliate endosymbionts.

Coping with temperature at the warm edge--patterns of thermal adaptation in the microbial eukaryote Paramecium caudatum.

BACKGROUND: Ectothermic organisms are thought to be severely affected by global warming since their physiological performance is directly dependent on temperature. Latitudinal and temporal variations in mean temperatures force ectotherms to adapt to these complex environmental conditions. Studies investigating current patterns of thermal adaptation among populations of different latitudes allow a prediction of the potential impact of prospective increases in environmental temperatures on their fitness. METHODOLOGY/PRINCIPAL FINDINGS: In this study, temperature reaction norms were ascertained among 18 genetically defined, natural clones of the microbial eukaryote Paramecium caudatum. These different clones have been isolated from 12 freshwater habitats along a latitudinal transect in Europe and from 3 tropical habitats (Indonesia). The sensitivity to increasing temperatures was estimated through the analysis of clone specific thermal tolerances and by relating those to current and predicted temperature data of their natural habitats. All investigated European clones seem to be thermal generalists with a broad thermal tolerance and similar optimum temperatures. The weak or missing co-variation of thermal tolerance with latitude does not imply local adaptation to thermal gradients; it rather suggests adaptive phenotypic plasticity among the whole European subpopulation. The tested Indonesian clones appear to be locally adapted to the less variable, tropical temperature regime and show higher tolerance limits, but lower tolerance breadths. CONCLUSIONS/SIGNIFICANCE: Due to the lack of local temperature adaptation within the European subpopulation, P. caudatum genotypes at the most southern edge of their geographic range seem to suffer from the predicted increase in magnitude and frequency of summer heat waves caused by climate change.

Thermal performance curves of Paramecium caudatum: a model selection approach.

The ongoing climate change has motivated numerous studies investigating the temperature response of various organisms, especially that of ectotherms. To correctly describe the thermal performance of these organisms, functions are needed which sufficiently fit to the complete optimum curve. Surprisingly, model-comparisons for the temperature-dependence of population growth rates of an important ectothermic group, the protozoa, are still missing. In this study, temperature reaction norms of natural isolates of the freshwater protist Paramecium caudatum were investigated, considering nearly the entire temperature range. These reaction norms were used to estimate thermal performance curves by applying a set of commonly used model functions. An information theory approach was used to compare models and to identify the best ones for describing these data. Our results indicate that the models which can describe negative growth at the high- and low-temperature branch of an optimum curve are preferable. This is a prerequisite for accurately calculating the critical upper and lower thermal limits. While we detected a temperature optimum of around 29 °C for all investigated clonal strains, the critical thermal limits were considerably different between individual clones. Here, the tropical clone showed the narrowest thermal tolerance, with a shift of its critical thermal limits to higher temperatures.

Intraspecific genetic variation in Paramecium revealed by mitochondrial cytochrome C oxidase I sequences.

Studies of intraspecific genetic diversity of ciliates, such as population genetics and biogeography, are particularly hampered by the lack of suitable DNA markers. For example, sequences of the non-coding ribosomal internal transcribed spacer (ITS) regions are often too conserved for intraspecific analyses. We have therefore identified primers for the mitochondrial cytochrome c oxidase I (COI) gene and applied them for intraspecific investigations in Paramecium caudatum and Paramecium multimicronucleatum. Furthermore, we obtained sequences of the ITS regions from the same strains and carried out comparative sequence analyses of both data sets. The mitochondrial sequences revealed substantially higher variation in both Paramecium species, with intraspecific divergences up to 7% in P. caudatum and 9.5% in P. multimicronucleatum. Moreover, an initial survey of the population structure discovered different mitochondrial haplotypes of P. caudatum in one pond, thereby demonstrating the potential of this genetic marker for population genetic analyses. Our primers successfully amplified the COI gene of other Paramecium. This is the first report of intraspecific variation in free-living protozoans based on mitochondrial sequence data. Our results show that the high variation in mitochondrial DNA makes it a suitable marker for intraspecific and population genetic studies.