Insights into the differential proteome landscape of a newly isolated Paramecium multimicronucleatum in response to cadmium stress.

Employing microbial systems for the bioremediation of contaminated waters represent a potential option, however, limited understanding of the underlying mechanisms hampers the implication of microbial-mediated bioremediation. The omics tools offer a promising approach to explore the molecular basis of the bioremediation process. Here, a mass spectrometry-based quantitative proteome profiling approach was conducted to explore the differential protein levels in cadmium-treated Paramecium multimicronucleatum. The Proteome Discoverer software was used to identify and quantify differentially abundant proteins. The proteome profiling generated 7,416 peptide spectral matches, yielding 2824 total peptides, corresponding to 989 proteins. The analysis revealed that 29 proteins exhibited significant (p ≤ 0.05) differential levels, including a higher abundance of 6 proteins and reduced levels of 23 proteins in Cd2+ treated samples. These differentially abundant proteins were associated with stress response, energy metabolism, protein degradation, cell growth, and hormone processing. Briefly, a comprehensive proteome profile in response to cadmium stress of a newly isolated Paramecium has been established that will be useful in future studies identifying critical proteins involved in the bioremediation of metals in ciliates. SIGNIFICANCE: Ciliates are considered a good biological indicator of chemical pollution and relatively sensitive to heavy metal contamination. A prominent ciliate, Paramecium is a promising candidate for the bioremediation of polluted water. The proteins related to metal resistance in Paramecium species are still largely unknown and need further exploration. In order to identify and reveal the proteins related to metal resistance in Paramecia, we have reported differential protein abundance in Paramecium multimicronucleatum in response to cadmium stress. The proteins found in our study play essential roles during stress response, hormone processing, protein degradation, energy metabolism, and cell growth. It seems likely that Paramecia are not a simple sponge for metals but they could also transform them into less toxic derivatives or by detoxification by protein binding. This data will be helpful in future studies to identify critical proteins along with their detailed mechanisms involved in the bioremediation and detoxification of metal ions in Paramecium species.

Uncoupling programmed DNA cleavage and repair scrambles the Paramecium somatic genome.

In the ciliate Paramecium, precise excision of numerous internal eliminated sequences (IESs) from the somatic genome is essential at each sexual cycle. DNA double-strands breaks (DSBs) introduced by the PiggyMac endonuclease are repaired in a highly concerted manner by the non-homologous end joining (NHEJ) pathway, illustrated by complete inhibition of DNA cleavage when Ku70/80 proteins are missing. We show that expression of a DNA-binding-deficient Ku70 mutant (Ku70-6E) permits DNA cleavage but leads to the accumulation of unrepaired DSBs. We uncoupled DNA cleavage and repair by co-expressing wild-type and mutant Ku70. High-throughput sequencing of the developing macronucleus genome in these conditions identifies the presence of extremities healed by de novo telomere addition and numerous translocations between IES-flanking sequences. Coupling the two steps of IES excision ensures that both extremities are held together throughout the process, suggesting that DSB repair proteins are essential for assembly of a synaptic precleavage complex.

Structural maintenance of chromosomes (SMC) proteins are required for DNA elimination in Paramecium.

Chromosome (SMC) proteins are a large family of ATPases that play important roles in the organization and dynamics of chromatin. They are central regulators of chromosome dynamics and the core component of condensin. DNA elimination during zygotic somatic genome development is a characteristic feature of ciliated protozoa such as Paramecium This process occurs after meiosis, mitosis, karyogamy, and another mitosis, which result in the formation of a new germline and somatic nuclei. The series of nuclear divisions implies an important role of SMC proteins in Paramecium sexual development. The relationship between DNA elimination and SMC has not yet been described. Here, we applied RNA interference, genome sequencing, mRNA sequencing, immunofluorescence, and mass spectrometry to investigate the roles of SMC components in DNA elimination. Our results show that SMC4-2 is required for genome rearrangement, whereas SMC4-1 is not. Functional diversification of SMC4 in Paramecium led to a formation of two paralogues where SMC4-2 acquired a novel, development-specific function and differs from SMC4-1. Moreover, our study suggests a competitive relationship between these two proteins.

Microbial impact to environmental toxicants Ni(II) and Co(II): Joint toxicity and cellular response in Paramecium.

Cobalt (Co) and Nickel (Ni) are increasingly found in our environment. We analysed their combined toxicity and uptake mechanisms in the early food chain by studying bacteria and the bacterivorous ciliate Paramecium as a primary consumer. We exposed both species to these metals to measure the toxicity, uptake and transfer of metals from bacteria to Paramecium. We found that Ni is more toxic than Co, and that toxicity increases for both metals when (i) food bacteria are absent and (ii) both metals are applied in combination. The cellular content in bacteria after exposure shows a concentration dependent bias for either Ni or Co. Comparing single treatment and joint exposure, bacteria show increased levels of both metals when these are both exposed. To imitate the basic level of the food chain, we fed these bacteria to paramecia. The cellular content shows a similar ratio of Nickel and Cobalt as in food bacteria. This is different to the direct application of both metals to paramecia, where Cobalt is enriched over Nickel. This indicates that bacteria can selectively pre-accumulate metals for introduction into the food chain. We also analysed the transcriptomic response of Paramecium to sublethal doses of Nickel and Cobalt to gain insight into their toxicity mechanisms. Gene ontology (GO) analysis indicates common deregulated pathways, such as ammonium transmembrane transport and ubiquitine-associated protein degradation. Many redox-related genes also show deregulation of gene expression, indicating cellular adaptation to increased RONS stress. This suggests that both metals may also target the same cellular pathways and this is consistent with the increased toxicity of both metals when used together. Our data reveal complex ecotoxicological pathways for these metals and highlights the different parameters for their fate in the ecosystem, in the food chain and their ecotoxicological risk after environmental contamination.

Comparative transcriptome and antioxidant biomarker response reveal molecular mechanisms to cope with zinc ion exposure in the unicellular eukaryote Paramecium.

The development of industry has resulted in excessive environmental zinc exposure which has caused various health problems in a wide range of organisms including humans. The mechanisms by which aquatic microorganisms respond to environmental zinc stress are still poorly understood. Paramecium, a well-known ciliated protozoan and a popular cell model in heavy metal stress response studies, was chosen as the test unicellular eukaryotic organism in the present research. In this work, Paramecium cf. multimicronucleatum cells were exposed in different levels of zinc ion (0.1 and 1.0 mg/L) for different periods of exposure (1 and 4 days), and then analyzed population growth, transcriptomic profiles and physiological changes in antioxidant enzymes to explore the toxicity and detoxification mechanisms during the zinc stress response. Results demonstrated that long-term zinc exposure could have restrained population growth in ciliates, however, the response mechanism to zinc exposure in ciliates is likely to show a dosage-dependent and time-dependent manner. The differentially expressed genes (DEGs) were identified the characters by high-throughput sequencing, which remarkably enriched in the phagosome, indicating that the phagosome pathway might mediate the uptake of zinc, while the pathways of ABC transporters and Na+/K+-transporting ATPase contributed to the efflux transport of excessive zinc ions and the maintenance of osmotic balance, respectively. The accumulation of zinc ions triggered a series of adverse effects, including damage to DNA and proteins, disturbance of mitochondrial function, and oxidative stress. In addition, we found that gene expression changed significantly for metal ion binding, energy metabolism, and oxidation-reduction processes. RT-qPCR of ten genes involved in important biological functions further validated the results of the transcriptome analysis. We also continuously monitored changes in activity of four antioxidant enzymes (SOD, CAT, POD and GSH-PX), all of which peaked on day 4 in cells subjected to zinc stress. Collectively, our results indicate that excessive environmental zinc exposure initially causes damage to cellular structure and function and then initiates detoxification mechanisms to maintain homeostasis in P. cf. multimicronucleatum cells.

Developmental mRNA clearance by PIWI-bound endo-siRNAs in Paramecium.

The clearance of untranslated mRNAs by Argonaute proteins is essential for embryonic development in metazoans. However, it is currently unknown whether similar processes exist in unicellular eukaryotes. The ciliate Paramecium tetraurelia harbors a vast array of PIWI-clade Argonautes involved in various small RNA (sRNA) pathways, many of which have not yet been investigated. Here, we investigate the function of a PIWI protein, Ptiwi08, whose expression is limited to a narrow time window during development, concomitant with the start of zygotic transcription. We show that Ptiwi08 acts in an endogenous small interfering RNA (endo-siRNA) pathway involved in the clearance of untranslated mRNAs. These endo-siRNAs are found in clusters that are strictly antisense to their target mRNAs and are a subset of siRNA-producing clusters (SRCs). Furthermore, the endo-siRNAs are 2'-O-methylated by Hen1 and require Dcr1 for their biogenesis. Our findings suggest that sRNA-mediated developmental mRNA clearance extends beyond metazoans and may be a more widespread mechanism than previously anticipated.

Methods for Paramecium tetraurelia ciliary membrane protein identification and function.

In this chapter we provide some tools to study the ciliary proteins that make it possible for Paramecium cells to swim by beating their cilia. These proteins include many ion channels, accessory proteins, peripheral proteins, structural proteins, rootlets of cilia, and enzymes. Some of these proteins are also found in the soma membrane, but their distinct and critical functions are in the cilia. Paramecium has 4000 or more cilia per cell, giving it an advantage for biochemical studies over cells that have one primarily cilium per cell. Nonetheless, a challenge for studies of many ciliary proteins in Paramecium is their low abundance. We discuss here several strategies to overcome this challenge and other challenges such as working with very large channel proteins. We also include for completeness other techniques that are critical to the study of swimming behavior, such as genetic crosses, recording of swimming patterns, electrical recordings, expression of very large channel proteins, RNA Interference, among others.

Chromatin remodeling is required for sRNA-guided DNA elimination in Paramecium.

Small RNAs mediate the silencing of transposable elements and other genomic loci, increasing nucleosome density and preventing undesirable gene expression. The unicellular ciliate Paramecium is a model to study dynamic genome organization in eukaryotic cells, given its unique feature of nuclear dimorphism. Here, the formation of the somatic macronucleus during sexual reproduction requires eliminating thousands of transposon remnants (IESs) and transposable elements scattered throughout the germline micronuclear genome. The elimination process is guided by Piwi-associated small RNAs and leads to precise cleavage at IES boundaries. Here we show that IES recognition and precise excision are facilitated by recruiting ISWI1, a Paramecium homolog of the chromatin remodeler ISWI. ISWI1 knockdown substantially inhibits DNA elimination, quantitatively similar to development-specific sRNA gene knockdowns but with much greater aberrant IES excision at alternative boundaries. We also identify key development-specific sRNA biogenesis and transport proteins, Ptiwi01 and Ptiwi09, as ISWI1 cofactors in our co-immunoprecipitation studies. Nucleosome profiling indicates that increased nucleosome density correlates with the requirement for ISWI1 and other proteins necessary for IES excision. We propose that chromatin remodeling together with small RNAs is essential for efficient and precise DNA elimination in Paramecium.

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.

Endosymbiotic Chlorella variabilis reduces mitochondrial number in the ciliate Paramecium bursaria.

Extant symbioses illustrate endosymbiosis is a driving force for evolution and diversification. In the ciliate Paramecium bursaria, the endosymbiotic alga Chlorella variabilis in perialgal vacuole localizes beneath the host cell cortex by adhesion between the perialgal vacuole membrane and host mitochondria. We investigated whether host mitochondria are also affected by algal endosymbiosis. Transmission electron microscopy of host cells showed fewer mitochondria beneath the algae-bearing host cell cortex than that of alga-free cells. To compare the density and distribution of host mitochondria with or without symbiotic algae, we developed a monoclonal antibody against Paramecium mitochondria. Immunofluorescence microscopy with the monoclonal antibody showed that the mitochondrial density of the algae-bearing P. bursaria was significantly lower than that of the alga-free cells. The total cell protein concentration of alga-free P. bursaria cells was approximately 1.8-fold higher than that of algae-bearing cells, and the protein content of mitochondria was significantly higher in alga-free cells than that in the algae-bearing cells. These results corresponded with those obtained by transmission electron and immunofluorescence microscopies. This paper shows that endosymbiotic algae affect reduced mitochondrial number in the host P. bursaria significantly.

Paramecium Polycomb repressive complex 2 physically interacts with the small RNA-binding PIWI protein to repress transposable elements.

Polycomb repressive complex 2 (PRC2) maintains transcriptionally silent genes in a repressed state via deposition of histone H3K27-trimethyl (me3) marks. PRC2 has also been implicated in silencing transposable elements (TEs), yet how PRC2 is targeted to TEs remains unclear. To address this question, we identified proteins that physically interact with the Paramecium enhancer-of-zeste Ezl1 enzyme, which catalyzes H3K9me3 and H3K27me3 deposition at TEs. We show that the Paramecium PRC2 core complex comprises four subunits, each required in vivo for catalytic activity. We also identify PRC2 cofactors, including the RNA interference (RNAi) effector Ptiwi09, which are necessary to target H3K9me3 and H3K27me3 to TEs. We find that the physical interaction between PRC2 and the RNAi pathway is mediated by a RING finger protein and that small RNA recruitment of PRC2 to TEs is analogous to the small RNA recruitment of H3K9 methylation SU(VAR)3-9 enzymes.

Mechanisms for establishing primary and secondary endosymbiosis in Paramecium.

Primary (eukaryote and procaryote) and secondary (eukaryote and eukaryote) endosymbioses are driving forces in eukaryotic cell evolution. These phenomena are still contributing to acquire new cell structures and functions. To understand mechanisms for establishment of each endosymbiosis, experiments that can induce endosymbiosis synchronously by mixing symbionts isolated from symbiont-bearing host cells and symbiont-free host cells are indispensable. Recent progress on endosymbiosis using Paramecium and their endonuclear symbiotic bacteria Holospora or symbiotic green alga Chlorella has been remarkable, providing excellent opportunities for elucidating host-symbiont interactions. These organisms are now becoming model organisms to know the mechanisms for establishing primary and secondary endosymbioses. Based on experiments of many researchers, we introduce how these endosymbionts escape from the host lysosomal fusion, how they migrate in the host cytoplasm to localize specific locations within the host, how their species specificity and strain specificity of the host cells are controlled, how their life cycles are controlled, how they escape from the host cell to infect more young host cell, how they affect the host viability and gene expression, what kind of substances are needed in these phenomena, and what changes had been induced in the symbiont and the host genomes.

Broad domains of histone marks in the highly compact Paramecium macronuclear genome.

The unicellular ciliate Paramecium contains a large vegetative macronucleus with several unusual characteristics, including an extremely high coding density and high polyploidy. As macronculear chromatin is devoid of heterochromatin, our study characterizes the functional epigenomic organization necessary for gene regulation and proper Pol II activity. Histone marks (H3K4me3, H3K9ac, H3K27me3) reveal no narrow peaks but broad domains along gene bodies, whereas intergenic regions are devoid of nucleosomes. Our data implicate H3K4me3 levels inside ORFs to be the main factor associated with gene expression, and H3K27me3 appears in association with H3K4me3 in plastic genes. Silent and lowly expressed genes show low nucleosome occupancy, suggesting that gene inactivation does not involve increased nucleosome occupancy and chromatin condensation. Because of a high occupancy of Pol II along highly expressed ORFs, transcriptional elongation appears to be quite different from that of other species. This is supported by missing heptameric repeats in the C-terminal domain of Pol II and a divergent elongation system. Our data imply that unoccupied DNA is the default state, whereas gene activation requires nucleosome recruitment together with broad domains of H3K4me3. In summary, gene activation and silencing in Paramecium run counter to the current understanding of chromatin biology.

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.

Evolutionary Plasticity of Mating-Type Determination Mechanisms in Paramecium aurelia Sibling Species.

The Paramecium aurelia complex, a group of morphologically similar but sexually incompatible sibling species, is a unique example of the evolutionary plasticity of mating-type systems. Each species has two mating types, O (Odd) and E (Even). Although O and E types are homologous in all species, three different modes of determination and inheritance have been described: genetic determination by Mendelian alleles, stochastic developmental determination, and maternally inherited developmental determination. Previous work in three species of the latter kind has revealed the key roles of the E-specific transmembrane protein mtA and its highly specific transcription factor mtB: type O clones are produced by maternally inherited genome rearrangements that inactivate either mtA or mtB during development. Here we show, through transcriptome analyses in five additional species representing the three determination systems, that mtA expression specifies type E in all cases. We further show that the Mendelian system depends on functional and nonfunctional mtA alleles, and identify novel developmental rearrangements in mtA and mtB which now explain all cases of maternally inherited mating-type determination. Epistasis between these genes likely evolved from less specific interactions between paralogs in the P. aurelia common ancestor, after a whole-genome duplication, but the mtB gene was subsequently lost in three P. aurelia species which appear to have returned to an ancestral regulation mechanism. These results suggest a model accounting for evolutionary transitions between determination systems, and highlight the diversity of molecular solutions explored among sibling species to maintain an essential mating-type polymorphism in cell populations.

The conformational cycle of a prototypical voltage-gated sodium channel.

Electrical signaling was a dramatic development in evolution, allowing complex single-cell organisms like Paramecium to coordinate movement and early metazoans like worms and jellyfish to send regulatory signals rapidly over increasing distances. But how are electrical signals generated in biology? In fact, voltage-gated sodium channels conduct sodium currents that initiate electrical signals in all kingdoms of life, from bacteria to man. They are responsible for generating the action potential in vertebrate nerve and muscle, neuroendocrine cells, and other cell types1,2. Because of the high level of conservation of their core structure, it is likely that their fundamental mechanisms of action are conserved as well. Here we describe the complete cycle of conformational changes that a bacterial sodium channel undergoes as it transitions from resting to activated/open and inactivated/closed states, based on high-resolution structural studies of a single sodium channel. We further relate this conformational cycle of the ancestral sodium channel to the function of its vertebrate orthologs. The strong conservation of amino acid sequence and three-dimensional structure suggests that this model, at a fundamental level, is relevant for both prokaryotic and eukaryotic sodium channels, as well as voltage-gated calcium and potassium channels.

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.

Comparison of Independent Evolutionary Origins Reveals Both Convergence and Divergence in the Metabolic Mechanisms of Symbiosis.

Through the merger of previously independent lineages, symbiosis promotes the acquisition of new traits and exploitation of inaccessible ecological niches [1, 2], driving evolutionary innovation and important ecosystem functions [3-6]. The transient nature of establishment makes study of symbiotic origins difficult, but experimental comparison of independent origins could reveal the degree of convergence in the underpinning mechanisms [7, 8]. We compared the metabolic mechanisms of two independent origins of Paramecium bursaria-Chlorella photosymbiosis [9-11] using a reciprocal metabolomic pulse-chase method. This showed convergent patterns of nutrient exchange and utilization for host-derived nitrogen in the Chlorella genotypes [12, 13] and symbiont-derived carbon in the P. bursaria genotypes [14, 15]. Consistent with a convergent primary nutrient exchange, partner-switched host-symbiont pairings were functional. Direct competition of hosts containing native or recombined symbionts against isogenic symbiont-free hosts showed that the fitness benefits of symbiosis for hosts increased with irradiance but varied by genotype. Global metabolism varied more between the Chlorella than the P. bursaria genotypes and suggested divergent mechanisms of light management. Specifically, the algal symbiont genotypes either produced photo-protective carotenoid pigments at high irradiance or more chlorophyll, resulting in corresponding differences in photosynthetic efficiency and non-photochemical quenching among host-symbiont pairings. These data suggest that the multiple origins of P. bursaria-Chlorella symbiosis use a convergent nutrient exchange, whereas other photosynthetic traits linked to functioning of photosymbiosis have diverged. Although convergence enables partner switching among diverse strains, phenotypic mismatches resulting from divergence of secondary symbiotic traits could mediate host-symbiont specificity in nature.

Photosynthetic Endosymbionts Benefit from Host's Phagotrophy, Including Predation on Potential Competitors.

In many endosymbioses, hosts have been shown to benefit from symbiosis, but it remains unclear whether intracellular endosymbionts benefit from their association with hosts [1, 2]. This makes it difficult to determine evolutionary mechanisms underlying cooperative behaviors between hosts and intracellular endosymbionts, such as mutual exchange of vital resources. Here, we investigate the fitness effects of symbiosis on the ciliate host Paramecium bursaria and on the algal endosymbiont Chlorella [3, 4], using experimental microcosms that include the free-living alga Chlamydomonas reinhardtii to mimic ecologically realistic conditions. We demonstrate that both host ciliate and the endosymbiotic algae gain fitness benefits from the symbiosis when another alga C. reinhardtii is present in the system. Specifically, the endosymbiotic Chlorella can grow as the host ciliate feeds and grows on C. reinhardtii, whereas the growth of free-living Chlorella is reduced by its competitor, C. reinhardtii. Thus, we propose that the endosymbiotic algae benefit from the host's phagotrophy, which allows the endosymbiont to access particulate nutrient sources and to indirectly prey on the potential competitors competing with its free-living counterparts. Even though the ecological contexts in which each partner receives its benefits differ, both partners would gain net fitness benefits in an ecological timescale. Thus, the cooperative behaviors can evolve through fitness feedback (partner fidelity feedback) between the host and the endosymbiont, without need for special partner control mechanisms. The proposed ecological and evolutionary mechanisms provide a basis for understanding cooperative resource exchanges in endosymbioses, including many photosynthetic endosymbioses widespread in aquatic ecosystems.

Dynamics of a magnetically rotated micro swimmer inspired by paramecium metachronal wave.

In the past few years, a significant body of research has been devoted to designing magnetic micron scale robotic systems for minimally invasive medicine. The motion of different microorganisms is the nature's solution for efficient propulsion of these swimmers. So far, there has been a considerable effort in designing micro swimmers based on the propulsion of bacteria while the motion of numerous other microorganisms has not been a source of inspiration for designing micro swimmers yet. Inspired by propulsion of Paramecium which is a ciliate microorganism, a novel micro swimmer is proposed in this article which is capable of cargo transport. This novel swimmer is composed of multiple equally spaced rigid loxodromic rods spanning the surface of a sphere which can carry a cargo placed inside it. The propulsion of this swimmer is influenced by the geometry of the swimmer (diameter, number of rods, cargo size), therefore, CFD simulations have been performed to investigate it. Finally, the dynamics of this swimmer is investigated analytically which sheds light into the complex dynamics of a swimmer with this geometry.