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.

Clarification of the Taxonomic Position of Paramecium caudatum Micronucleus Symbionts.

Bacteria of genus Holospora (order Holosporales, class Alphaproteobacteria) are obligate intranuclear symbionts of ciliates Paramecium spp. with strict host species and nuclear (macronucleus or micronucleus) specificity. However, three species under study Holospora undulata, Holospora elegans and 'Holospora recta' occupy the same ecological niche-micronucleus of Paramecium caudatum and demonstrate some differences in morphology of infectious form. The genetic diversity of holosporas by rrs and rpoB sequence analysis was determined. Phylogenetic and phylogenomic analysis of Holospora spp., as well as some phenotypic features indicate that there is no distinctive difference supporting studied micronuclear endosymbionts as distinct species. Therefore, Holospora elegans and 'Holospora recta' should be considered subspecies of Holospora undulata (ex Haffkine 1890) Gromov and Ossipov 1981, which was described first. Thus, we confirmed the evolutionary aspects of the development of symbiotic relationships: holosporas have a strict specificity to the host species and the type of nucleus.

Parasitism and host dispersal plasticity in an aquatic model system.

Dispersal is a central determinant of spatial dynamics in communities and ecosystems, and various ecological factors can shape the evolution of constitutive and plastic dispersal behaviours. One important driver of dispersal plasticity is the biotic environment. Parasites, for example, influence the internal condition of infected hosts and define external patch quality. Thus, state-dependent dispersal may be determined by infection status and context-dependent dispersal by the abundance of infected hosts in the population. A prerequisite for such dispersal plasticity to evolve is a genetic basis on which natural selection can act. Using interconnected microcosms, we investigated dispersal in experimental populations of the freshwater protist Paramecium caudatum in response to the bacterial parasite Holospora undulata. For a collection of 20 natural host strains, we found substantial variation in constitutive dispersal and to a lesser degree in dispersal plasticity. First, infection tended to increase or decrease dispersal relative to uninfected controls, depending on strain identity, indicative of state-dependent dispersal plasticity. Infection additionally decreased host swimming speed compared to the uninfected counterparts. Second, for certain strains, there was a weak negative association between dispersal and infection prevalence, such that uninfected hosts dispersed less when infection was more frequent in the population, indicating context-dependent dispersal plasticity. Future experiments may test whether the observed differences in dispersal plasticity are sufficiently strong to be picked up by natural selection. The evolution of dispersal plasticity as a strategy to mitigate parasite effects spatially may have important implications for epidemiological dynamics.

An evolutionary trade-off between parasite virulence and dispersal at experimental invasion fronts.

Exploitative parasites are predicted to evolve in highly connected populations or in expanding epidemics. However, many parasites rely on host dispersal to reach new populations, potentially causing conflict between local transmission and global spread. We performed experimental range expansions in interconnected microcosms of the protozoan Paramecium caudatum, allowing natural dispersal of hosts infected with the bacterial parasite Holospora undulata. Parasites from range front treatments facilitated host dispersal and were less virulent, but also invested less in horizontal transmission than parasites from range cores. These differences were consistent with parameter estimates derived from an epidemiological model fitted on population-level time-series data. Our results illustrate how dispersal selection can have profound consequences for the evolution of parasite life history and virulence. Decrypting the eco-evolutionary processes that shape parasite 'dispersal syndromes' may be important for the management of spreading epidemics in changing environments, biological invasions or in other spatial non-equilibrium settings.

DNA-Based Sensor Particles Enable Measuring Light Intensity in Single Cells.

"Lab on a particle" architecture is employed in designing a light nanosensor. Light-sensitive protecting groups are installed on DNA, which is encapsulated in silica particles, qualifying as a self-sufficient light sensor. The nanosensors allow measuring light intensity and duration in very small volumes, such as single cells, and store the irradiation information until readout.

Long-term selection experiment produces breakdown of horizontal transmissibility in parasite with mixed transmission mode.

Evolutionary transitions from parasitism toward beneficial or mutualistic associations may encompass a change from horizontal transmission to (strict) vertical transmission. Parasites with both vertical and horizontal transmission are amendable to study factors driving such transitions. In a long-term experiment, microcosm populations of the protozoan Paramecium caudatum and its bacterial parasite Holospora undulata were exposed to three growth treatments, manipulating vertical transmission opportunities over ca. 800 host generations. In inoculation tests, horizontal transmission propagules produced by parasites from a "high-growth" treatment, with elevated host division rates increasing levels of parasite vertical transmission, showed a near-complete loss of infectivity. A similar reduction was observed for parasites from a treatment alternating between high growth and low growth (i.e., low levels of population turn-over). Parasites from a low-growth treatment had the highest infectivity on all host genotypes tested. Our results complement previous findings of reduced investment in horizontal transmission and increased vertical transmissibility of high-growth parasites. We explain the loss of horizontal transmissibility by epidemiological feedbacks and resistance evolution, reducing the frequency of susceptible hosts in the population and thereby decreasing the selective advantage of horizontal transmission. This illustrates how environmental conditions may push parasites with a mixed transmission mode toward becoming vertically transmitted nonvirulent symbionts.

Identification of two nickel ion-induced genes, NCI16 and PcGST1, in Paramecium caudatum.

Here, we describe the isolation of two nickel-induced genes in Paramecium caudatum, NCI16 and PcGST1, by subtractive hybridization. NCI16 encoded a predicted four-transmembrane domain protein (∼16 kDa) of unknown function, and PcGST1 encoded glutathione S-transferase (GST; ∼25 kDa) with GST and glutathione peroxidase (GPx) activities. Exposing cells to cobalt chloride also caused the moderate upregulation of NCI16 and PcGST1 mRNAs. Both nickel sulfate and cobalt chloride dose dependently induced NCI16 and PcGST1 mRNAs, but with different profiles. Nickel treatment caused a continuous increase in PcGST1 and NCI16 mRNA levels for up to 3 and 6 days, respectively, and a notable increase in H2O2 concentrations in P. caudatum. NCI16 expression was significantly enhanced by incubating cells with H2O2, implying that NCI16 induction in the presence of nickel ions is caused by reactive oxygen species (ROS). On the other hand, PcGST1 was highly induced by the antioxidant tert-butylhydroquinone (tBHQ) but not by H2O2, suggesting that different mechanisms mediate the induction of NCI16 and PcGST1. We introduced a luciferase reporter vector with an ∼0.42-kb putative PcGST1 promoter into cells and then exposed the transformants to nickel sulfate. This resulted in significant luciferase upregulation, indicating that the putative PcGST1 promoter contains a nickel-responsive element. Our nickel-inducible system also may be applicable to the efficient expression of proteins that are toxic to host cells or require temporal control.

Variability in secondary structure of 18S ribosomal RNA as topological marker for identification of Paramecium species.

Besides cytological and molecular applications, Paramecium is being used in water quality assessment and for determination of saprobic levels. An unambiguous identification of these unicellular eukaryotes is not only essential, but its ecological diversity must also be explored in the local environment. 18SrRNA genes of all the strains of Paramecium species isolated from waste water were amplified, cloned and sequenced. Phylogenetic comparison of the nucleotide sequences of these strains with 23 closely related Paramecium species from GenBank Database enabled identification of Paramecium multimicronucleatum and Paramecium jenningsi. Some isolates did not show significant close association with other Paramecium species, and because of their unique position in the phylogenetic tree, they were considered new to the field. In the present report, these isolates are being designated as Paramecium caudatum pakistanicus. In this article, secondary structure of 18SrRNA has also been analyzed as an additional and perhaps more reliable topological marker for species discrimination and for determining possible phylogenetic relationship between the ciliate species. On the basis of comparison of secondary structure of 18SrRNA of various isolated Paramacium strains, and among Paramecium caudatum pakistanicus, Tetrahymena thermophila, Drosophila melanogaster, and Homo sapiens, it can be deduced that variable regions are more helpful in differentiating the species at interspecific level rather than at intraspecific level. It was concluded that V3 was the least variable region in all the organisms, V2 and V7 were the longest expansion segments of D. melanogaster and there was continuous mutational bias towards G.C base pairing in H. sapiens.

Insights into three whole-genome duplications gleaned from the Paramecium caudatum genome sequence.

Paramecium has long been a model eukaryote. The sequence of the Paramecium tetraurelia genome reveals a history of three successive whole-genome duplications (WGDs), and the sequences of P. biaurelia and P. sexaurelia suggest that these WGDs are shared by all members of the aurelia species complex. Here, we present the genome sequence of P. caudatum, a species closely related to the P. aurelia species group. P. caudatum shares only the most ancient of the three WGDs with the aurelia complex. We found that P. caudatum maintains twice as many paralogs from this early event as the P. aurelia species, suggesting that post-WGD gene retention is influenced by subsequent WGDs and supporting the importance of selection for dosage in gene retention. The availability of P. caudatum as an outgroup allows an expanded analysis of the aurelia intermediate and recent WGD events. Both the Guanine+Cytosine (GC) content and the expression level of preduplication genes are significant predictors of duplicate retention. We find widespread asymmetrical evolution among aurelia paralogs, which is likely caused by gradual pseudogenization rather than by neofunctionalization. Finally, cases of divergent resolution of intermediate WGD duplicates between aurelia species implicate this process acts as an ongoing reinforcement mechanism of reproductive isolation long after a WGD event.

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.

[Analysis of G.F. Gause experimental time-series by means of continuous-time models].

Four population dynamics models, namely Verhulst, Gompertz, Rosenzweig, and Svirezhev ones, have been used to approximate two well-known time-series of Paramecia aurelia and P. caudatum population size (Gause, 1934). The parameters are estimated for each of the models by the least-square method (with global fitting) in two different ways: with and without an additional upper bound for a parameter value. In the latter (traditional) case, when the deviations of theoretical (model) trajectories from experimental time-series have been tested for normality (Kolmogorov-Smirnov test, Shapiro-Wilk test) with zero average, and for the presence/absence of serial correlations (Durbin-Watson criteria), the best results are obtained for the Gompertz and Verhulst models. In the former, more realistic, case (when we impose an additional constraint that the parameter meaning the carrying capacity of the environment has to be greater than any element in the sample), the best results are observed for the Gompertz model. Under this constraint, the canonical technique for deviation analysis can be applied in a restricted version only.

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.

Genetic influence on disease spread following arrival of infected carriers.

Epidemiology in host meta-populations depends on parasite ability to disperse between, establish and persist in distinct sub-populations of hosts. We studied the genetic factors determining the short-term establishment, and long-term maintenance, of pathogens introduced by infected hosts (i.e. carriers) into recipient populations. We used experimental populations of the freshwater ciliate Paramecium caudatum and its bacterial parasite Holospora undulata. Parasite short-term spread (approximately one horizontal transmission cycle) was affected mainly by carrier genotype, and its interactions with parasite and recipient genotypes. By contrast, parasite longer term spread (2-3 horizontal transmission cycles) was mostly determined by parasite isolate. Importantly, measures of parasite short-term success (reproductive number, R) were not good predictors for longer term prevalence, probably because of the specific interactions between host and parasite genotypes. Analogous to variation in vectorial capacity and super-spreader occurrence, two crucial components of epidemiology, we show that carrier genotype can also affect disease spread within meta-populations.

Reverse evolution: selection against costly resistance in disease-free microcosm populations of Paramecium caudatum.

Evolutionary costs of parasite resistance arise if genes conferring resistance reduce fitness in the absence of parasites. Thus, parasite-mediated selection may lead to increased resistance and a correlated decrease in fitness, whereas relaxed parasite-mediated selection may lead to reverse evolution of increased fitness and a correlated decrease in resistance. We tested this idea in experimental populations of the protozoan Paramecium caudatum and the parasitic bacterium Holospora undulata. After eight years, resistance to infection and asexual reproduction were compared among paramecia from (1) "infected" populations, (2) uninfected "naive" populations, and (3) previously infected, parasite-free "recovered" populations. Paramecia from "infected" populations were more resistant (+12%), but had lower reproduction (-15%) than "naive" paramecia, indicating an evolutionary trade-off between resistance and fitness. Recovered populations showed similar reproduction to naive populations; however, resistance of recently (<3 years) recovered populations was similar to paramecia from infected populations, whereas longer (>3 years) recovered populations were as susceptible as naive populations. This suggests a weak, convex trade-off between resistance and fitness, allowing recovery of fitness, without complete loss of resistance, favoring the maintenance of a generalist strategy of intermediate fitness and resistance. Our results indicate that (co)evolution with parasites can leave a genetic signature in disease-free populations.

New details indicated by different stainings during conjugation of ciliated protozoa Paramecium.

During conjugation of Paramecium caudatum, nuclear events occur in a scheduled program. Morphological studies on nuclear behavior during conjugation of P. caudatum have been performed since the end of the 19th century. Here we report on new details concerning the conjugation of P. caudatum through the staining of conjugating cells with protargol, carbol fuchsin solution, Hoechst 33342 and immunofluorescence labeling with monoclonal antibody of anti-α tubulin. 1) The crescent nucleus is a characteristic of the meiotic prophase of P. caudatum, has an unstained area. We stained this area with protargol, which was separated from the chromatin area and was not detected by the other stainings. 2) In regards to the four meiotic products, it has long been considered that only one product enters the paroral cone region (PC) and survives after meiosis. However, our protargol and immunofluorescence labeling results indicated that PC entrance of the meiotic product happened before the completion of meiosis instead of after. 3) In our previous study, protargol staining indicated the presence of a swollen structure around the central part of the "U" and "V" shaped spindles connecting the two types of prospective pronuclei. However, immunofluorescence labeling with anti-α tubulin antibodies gave a different image from protargol. All these observations form the basis for further studies of their molecular mechanisms.

Gene transport and expression by arginine-rich cell-penetrating peptides in Paramecium.

Owing to the cell membrane barriers, most macromolecules and hydrophilic molecules could not freely enter into living cells. However, cell-penetrating peptides (CPPs) have been discovered that can translocate themselves and associate cargoes into the cytoplasm. In this study, we demonstrate that three arginine-rich CPPs (SR9, HR9 and PR9) can form stable complexes with plasmid DNA at the optimized nitrogen/phosphate ratio of 3 and deliver plasmid DNA into Paramecium caudatum in a noncovalent manner. Accordingly, the transported plasmid encoding the green fluorescent protein (GFP) gene could be expressed in cells functionally assayed at both the protein and DNA levels. The efficiency of gene delivery varied among these CPPs in the order of HR9>PR9>SR9. In addition, these CPPs and CPP/DNA complexes were not cytotoxic in Paramecium detected by the 3-(4,5-dimethylthiazol-2-yl)-2,5-diohenyltetrazolium bromide (MTT) assay. Thus, these results suggest that the functionality of arginine-rich CPPs offers an efficient and safe tool for transgenesis in eukaryotic protozoans.

Localization of stationary pronuclei during conjugation of Paramecium as indicated by immunofluorescence staining.

After the third prezygotic division during conjugation of Paramecium caudatum, migratory and stationary pronuclei are produced. The migratory pronuclei remain in the paroral region tightly against the conjugating boundaries; while the stationary pronuclei are located beside the migratory pronuclei. To date, however, it is not clear what causes this close side-by-side localization between migratory and stationary pronuclei. In the current study, immunofluorescence staining with monoclonal antibody of anti-α tubulin indicated that ''U'' or ''V'' shaped spindles connected the migratory and stationary pronuclei during the third prezygotic division. This observation accounts for the close localization between these two types of pronuclei.

The mitochondrial genome sequence of the ciliate Paramecium caudatum reveals a shift in nucleotide composition and codon usage within the genus Paramecium.

BACKGROUND: Despite the fact that the organization of the ciliate mitochondrial genome is exceptional, only few ciliate mitochondrial genomes have been sequenced until today. All ciliate mitochondrial genomes are linear. They are 40 kb to 47 kb long and contain some 50 tightly packed genes without introns. Earlier studies documented that the mitochondrial guanine + cytosine contents are very different between Paramecium tetraurelia and all studied Tetrahymena species. This raises the question of whether the high mitochondrial G+C content observed in P. tetraurelia is a characteristic property of Paramecium mtDNA, or whether it is an exception of the ciliate mitochondrial genomes known so far. To test this question, we determined the mitochondrial genome sequence of Paramecium caudatum and compared the gene content and sequence properties to the closely related P. tetraurelia. RESULTS: The guanine + cytosine content of the P. caudatum mitochondrial genome was significantly lower than that of P. tetraurelia (22.4% vs. 41.2%). This difference in the mitochondrial nucleotide composition was accompanied by significantly different codon usage patterns in both species, i.e. within P. caudatum clearly A/T ending codons dominated, whereas for P. tetraurelia the synonymous codons were more balanced with a higher number of G/C ending codons. Further analyses indicated that the nucleotide composition of most members of the genus Paramecium resembles that of P. caudatum and that the shift observed in P. tetraurelia is restricted to the P. aurelia species complex. CONCLUSIONS: Surprisingly, the codon usage bias in the P. caudatum mitochondrial genome, exemplified by the effective number of codons, is more similar to the distantly related T. pyriformis and other single-celled eukaryotes such as Chlamydomonas, than to the closely related P. tetraurelia. These differences in base composition and codon usage bias were, however, not reflected in the amino acid composition. Most probably, the observed picture is best explained by a hitherto unknown (neutral or adaptive) mechanism that increased the guanine + cytosine content in P. tetraurelia mtDNA on the one hand, and strong purifying selection on the ancestral amino acid composition on the other hand. These contradicting forces are counterbalanced by a considerably altered codon usage pattern.

Temporal variation in temperature determines disease spread and maintenance in Paramecium microcosm populations.

The environment is rarely constant and organisms are exposed to temporal and spatial variations that impact their life histories and inter-species interactions. It is important to understand how such variations affect epidemiological dynamics in host-parasite systems. We explored effects of temporal variation in temperature on experimental microcosm populations of the ciliate Paramecium caudatum and its bacterial parasite Holospora undulata. Infected and uninfected populations of two P. caudatum genotypes were created and four constant temperature treatments (26°C, 28°C, 30°C and 32°C) compared with four variable treatments with the same mean temperatures. Variable temperature treatments were achieved by alternating populations between permissive (23°C) and restrictive (35°C) conditions daily over 30 days. Variable conditions and high temperatures caused greater declines in Paramecium populations, greater fluctuations in population size and higher incidence of extinction. The additional effect of parasite infection was additive and enhanced the negative effects of the variable environment and higher temperatures by up to 50 per cent. The variable environment and high temperatures also caused a decrease in parasite prevalence (up to 40%) and an increase in extinction (absence of detection) (up to 30%). The host genotypes responded similarly to the different environmental stresses and their effect on parasite traits were generally in the same direction. This work provides, to our knowledge, the first experimental demonstration that epidemiological dynamics are influenced by environmental variation. We also emphasize the need to consider environmental variance, as well as means, when trying to understand, or predict population dynamics or range.

Dynamic behaviour of stationary pronuclei during their positioning in Paramecium caudatum.

During conjugation of Paramecium caudatum, there are two well-known stages when nuclear migration occurs. What happens to the nuclei is closely related to their localisations in cells. The first of these stages is the entrance of one meiotic product into the paroral region. This nucleus survives, while the remaining three outside this area degenerate. The second stage is the antero-posterior localisation of eight synkaryon division products. Four posterior nuclei are differentiated into macronuclear anlagen, whereas four anterior nuclei remain as the presumptive micronuclei. In this experiment, the process of the third prezygotic division of P. caudatum was studied with the help of protargol staining. Here, a third nuclear migration was discovered. By two spindle turnings and two spindle elongations, stationary pronuclei were positioned near migratory pronuclei. This positioning of stationary pronuclei could shorten the distance for transferred migratory pronuclei to recognise and reach the stationary pronuclei. This fosters the synkaryon formation of P. caudatum.