Genetic basis for the establishment of endosymbiosis in Paramecium

Complex evolutionary patterns within the tubulin gene family of ciliates, unicellular eukaryotes with diverse microtubular structures

BACKGROUND

Tubulins are major components of the eukaryotic cytoskeletons that are crucial in many cellular processes. Ciliated protists comprise one of the oldest eukaryotic lineages possessing cilia over their cell surface and assembling many diverse microtubular structures. As such, ciliates are excellent model organisms to clarify the origin and evolution of tubulins in the early stages of eukaryote evolution. Nonetheless, the evolutionary history of the tubulin subfamilies within and among ciliate classes is unclear.

RESULTS

We analyzed the evolutionary pattern of ciliate tubulin gene family based on genomes/transcriptomes of 60 species covering 10 ciliate classes. Results showed: (1) Six tubulin subfamilies (α_Tub, β_Tub, γ_Tub, δ_Tub, ε_Tub, and ζ_Tub) originated from the last eukaryotic common ancestor (LECA) were observed within ciliates. Among them, α_Tub, β_Tub, and γ_Tub were present in all ciliate species, while δ_Tub, ε_Tub, and ζ_Tub might be independently lost in some species. (2) The evolutionary history of the tubulin subfamilies varied. Evolutionary history of ciliate γ_Tub, δ_Tub, ε_Tub, and ζ_Tub showed a certain degree of consistency with the phylogeny of species after the divergence of ciliate classes, while the evolutionary history of ciliate α_Tub and β_Tub varied among different classes. (3) Ciliate α- and β-tubulin isoforms could be classified into an "ancestral group" present in LECA and a "divergent group" containing only ciliate sequences. Alveolata-specific expansion events probably occurred within the "ancestral group" of α_Tub and β_Tub. The "divergent group" might be important for ciliate morphological differentiation and wide environmental adaptability. (4) Expansion events of the tubulin gene family appeared to be consistent with whole genome duplication (WGD) events in some degree. More Paramecium-specific tubulin expansions were detected than Tetrahymena-specific ones. Compared to other Paramecium species, the Paramecium aurelia complex underwent a more recent WGD which might have experienced more tubulin expansion events.

CONCLUSIONS

Evolutionary history among different tubulin gene subfamilies seemed to vary within ciliated protists. And the complex evolutionary patterns of tubulins among different ciliate classes might drive functional diversification. Our investigation provided meaningful information for understanding the evolution of tubulin gene family in the early stages of eukaryote evolution.

Exploring the landscape of symbiotic diversity and distribution in unicellular ciliated protists

BACKGROUND

The eukaryotic-bacterial symbiotic system plays an important role in various physiological, developmental, and evolutionary processes. However, our current understanding is largely limited to multicellular eukaryotes without adequate consideration of diverse unicellular protists, including ciliates.

RESULTS

To investigate the bacterial profiles associated with unicellular organisms, we collected 246 ciliate samples spanning the entire Ciliophora phylum and conducted single-cell based metagenome sequencing. This effort has yielded the most extensive collection of bacteria linked to unicellular protists to date. From this dataset, we identified 883 bacterial species capable of cohabiting with ciliates, unveiling the genomes of 116 novel bacterial cohabitants along with 7 novel archaeal cohabitants. Highlighting the intimate relationship between ciliates and their cohabitants, our study unveiled that over 90% of ciliates coexist with bacteria, with individual hosts fostering symbiotic relationships with multiple bacteria concurrently, resulting in the observation of seven distinct symbiotic patterns among bacteria. Our exploration of symbiotic mechanisms revealed the impact of host digestion on the intracellular diversity of cohabitants. Additionally, we identified the presence of eukaryotic-like proteins in bacteria as a potential contributing factor to their resistance against host digestion, thereby expanding their potential host range.

CONCLUSIONS

As the first large-scale analysis of prokaryotic associations with ciliate protists, this study provides a valuable resource for future research on eukaryotic-bacterial symbioses. Video Abstract.

Cellular interactions and evolutionary origins of endosymbiotic relationships with ciliates

As unicellular predators, ciliates engage in close associations with diverse microbes, laying the foundation for the establishment of endosymbiosis. Originally heterotrophic, ciliates demonstrate the ability to acquire phototrophy by phagocytizing unicellular algae or by sequestering algal plastids. This adaptation enables them to gain photosynthate and develop resistance to unfavorable environmental conditions. The integration of acquired phototrophy with intrinsic phagotrophy results in a trophic mode known as mixotrophy. Additionally, ciliates can harbor thousands of bacteria in various intracellular regions, including the cytoplasm and nucleus, exhibiting species specificity. Under prolonged and specific selective pressure within hosts, bacterial endosymbionts evolve unique lifestyles and undergo particular reductions in metabolic activities. Investigating the research advancements in various endosymbiotic cases within ciliates will contribute to elucidate patterns in cellular interaction and unravel the evolutionary origins of complex traits.

Elucidating Multiple Electron-Transfer Pathways for Metavanadate Bioreduction by Actinomycetic Streptomyces microflavus

While microbial reduction has gained widespread recognition for efficiently remediating environments polluted by toxic metavanadate [V(V)], the pool of identified V(V)-reducing strains remains rather limited, with the vast majority belonging to bacteria and fungi. This study is among the first to confirm the V(V) reduction capability of Streptomyces microflavus, a representative member of ubiquitous actinomycetes in environment. A V(V) removal efficiency of 91.0 ± 4.35% was achieved during 12 days of operation, with a maximum specific growth rate of 0.073 d-1. V(V) was bioreduced to insoluble V(IV) precipitates. V(V) reduction took place both intracellularly and extracellularly. Electron transfer was enhanced during V(V) bioreduction with increased electron transporters. The electron-transfer pathways were revealed through transcriptomic, proteomic, and metabolomic analyses. Electrons might flow either through the respiratory chain to reduce intracellular V(V) or to cytochrome c on the outer membrane for extracellular V(V) reduction. Soluble riboflavin and quinone also possibly mediated extracellular V(V) reduction. Glutathione might deliver electrons for intracellular V(V) reduction. Bioaugmentation of the aquifer sediment with S. microflavus accelerated V(V) reduction. The strain could successfully colonize the sediment and foster positive correlations with indigenous microorganisms. This study offers new microbial resources for V(V) bioremediation and improve the understanding of the involved molecular mechanisms.

Paramecium Genetics, Genomics, and Evolution

The ciliate genus Paramecium served as one of the first model systems in microbial eukaryotic genetics, contributing much to the early understanding of phenomena as diverse as genome rearrangement, cryptic speciation, cytoplasmic inheritance, and endosymbiosis, as well as more recently to the evolution of mating types, introns, and roles of small RNAs in DNA processing. Substantial progress has recently been made in the area of comparative and population genomics. Paramecium species combine some of the lowest known mutation rates with some of the largest known effective populations, along with likely very high recombination rates, thereby harboring a population-genetic environment that promotes an exceptionally efficient capacity for selection. As a consequence, the genomes are extraordinarily streamlined, with very small intergenic regions combined with small numbers of tiny introns. The subject of the bulk of Paramecium research, the ancient Paramecium aurelia species complex, is descended from two whole-genome duplication events that retain high degrees of synteny, thereby providing an exceptional platform for studying the fates of duplicate genes. Despite having a common ancestor dating to several hundred million years ago, the known descendant species are morphologically indistinguishable, raising significant questions about the common view that gene duplications lead to the origins of evolutionary novelties.

Taming the perils of photosynthesis by eukaryotes: constraints on endosymbiotic evolution in aquatic ecosystems

An ancestral eukaryote acquired photosynthesis by genetically integrating a cyanobacterial endosymbiont as the chloroplast. The chloroplast was then further integrated into many other eukaryotic lineages through secondary endosymbiotic events of unicellular eukaryotic algae. While photosynthesis enables autotrophy, it also generates reactive oxygen species that can cause oxidative stress. To mitigate the stress, photosynthetic eukaryotes employ various mechanisms, including regulating chloroplast light absorption and repairing or removing damaged chloroplasts by sensing light and photosynthetic status. Recent studies have shown that, besides algae and plants with innate chloroplasts, several lineages of numerous unicellular eukaryotes engage in acquired phototrophy by hosting algal endosymbionts or by transiently utilizing chloroplasts sequestrated from algal prey in aquatic ecosystems. In addition, it has become evident that unicellular organisms engaged in acquired phototrophy, as well as those that feed on algae, have also developed mechanisms to cope with photosynthetic oxidative stress. These mechanisms are limited but similar to those employed by algae and plants. Thus, there appear to be constraints on the evolution of those mechanisms, which likely began by incorporating photosynthetic cells before the establishment of chloroplasts by extending preexisting mechanisms to cope with oxidative stress originating from mitochondrial respiration and acquiring new mechanisms.

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

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

Synthetic symbiosis between a cyanobacterium and a ciliate toward novel chloroplast-like endosymbiosis

Chloroplasts are thought to have co-evolved through endosymbiosis, after a cyanobacterial-like prokaryote was engulfed by a eukaryotic cell; however, it is impossible to observe the process toward chloroplasts. In this study, we constructed an experimental symbiosis model to observe the initial stage in the process from independent organisms to a chloroplast-like organelle. Our system of synthetic symbiosis is capable of long-term coculture of two model organisms: a cyanobacterium (Synechocystis sp. PCC6803) as a symbiont and a ciliate (Tetrahymena thermophila) as a host with endocytic ability. The experimental system was clearly defined, because we used a synthetic medium and the cultures were shaken to avoid spatial complexity. We determined the experimental conditions for sustainable coculture, by analyzing population dynamics using a mathematical model. We experimentally demonstrated that the coculture was sustainable for at least 100 generations, through serial transfers. Moreover, we found that cells isolated after the serial transfer improved the probability of coexistence of both species without extinction in re-coculture. The constructed system will be useful for understanding the initial stage of primary endosymbiosis from cyanobacteria to chloroplasts, i.e., the origin of algae and plants.

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.

An extremely streamlined macronuclear genome in the free-living protozoan Fabrea salina

Ciliated protists are among the oldest unicellular organisms with a heterotrophic lifestyle and share a common ancestor with Plantae. Unlike any other eukaryotes, there are two distinct nuclei in ciliates with separate germline and somatic cell functions. Here, we assembled a near-complete macronuclear genome of Fabrea salina, which belongs to one of the oldest clades of ciliates. Its extremely minimized genome (18.35 Mb) is the smallest among all free-living heterotrophic eukaryotes and exhibits typical streamlined genomic features, including high gene density, tiny introns, and shrinkage of gene paralogs. Gene families involved in hypersaline stress resistance, DNA replication proteins, and mitochondrial biogenesis are expanded, and the accumulation of phosphatidic acid may play an important role in resistance to high osmotic pressure. We further investigated the morphological and transcriptomic changes in the macronucleus during sexual reproduction and highlighted the potential contribution of macronuclear residuals to this process. We believe that the minimized genome generated in this study provides novel insights into the genome streamlining theory and will be an ideal model to study the evolution of eukaryotic heterotrophs.

Toxic Effects of Copper Nanoparticles on Paramecium bursaria–Chlorella Symbiotic System

Although many reports have demonstrated that nanoparticles can have a negative effect on aquatic organisms, the toxic effects on symbiotic organisms remain poorly understood. The present study conducts ultrastructure, enzyme activity, and transcriptomics to assess the toxic effects to the Paramecium bursaria–Chlorella symbiotic system from exposure to copper nanoparticles (CuNPs) for 24 h. We found that in both the host and symbiotic algae, CuNP exposure induced high reactive oxygen species level, which leads to oxidative damage and energy metabolism disorder. Moreover, transmission electron micrographs (TEMs) showed that the symbiotic algae in the cytoplasm of P. bursaria were enveloped in the digestive vacuole and digested, and the level of acid phosphatase activity increased significantly within 24 h, which indicated that the stability of the symbiotic system was affected after CuNP exposure. We speculated that the increased energy demand in the host and symbiotic algae resulted from oxidative stress, precipitating the decrease of the photosynthetic products provided to the host, the digestion of the symbiont, and the destruction of the stable symbiotic relationship. The study provides the first insight into the mechanisms of nanoparticles’ toxicity to the symbiotic relationship in the ecosystem, which may help to understand the environmental effects and toxicological mechanisms of nanoparticles.

Natural Shifts in Endosymbionts' Occurrence and Relative Frequency in Their Ciliate Host Population

The role of bacterial endosymbionts harbored by heterotrophic Paramecium species is complex. Obligate intracellular bacteria supposedly always inflict costs as the host is the only possible provider of resources. However, several experimental studies have shown that paramecia carrying bacterial endosymbionts can benefit from their infection. Here, we address the question which endosymbionts occur in natural paramecia populations isolated from a small lake over a period of 5 years and which factors might explain observed shifts and persistence in the symbionts occurrence. One hundred and nineteen monoclonal strains were investigated and approximately two-third harbored intracellular bacteria. The majority of infected paramecia carried the obligate endosymbiotic "Candidatus Megaira polyxenophila", followed by Caedimonas varicaedens, and Holospora undulata. The latter was only detected in a single strain. While "Ca. M. polyxenophila" was observed in seven out of 13 samplings, C. varicaedens presence was limited to a single sampling occasion. After the appearance of C. varicaedens, "Ca. M. polyxenophila" prevalence dramatically dropped with some delay but recovered to original levels at the end of our study. Potential mechanisms explaining these observations include differences in infectivity, host range, and impact on host fitness as well as host competitive capacities. Growth experiments revealed fitness advantages for infected paramecia harboring "Ca. M. polyxenophila" as well as C. varicaedens. Furthermore, we showed that cells carrying C. varicaedens gain a competitive advantage from the symbiosis-derived killer trait. Other characteristics like infectivity and overlapping host range were taken into consideration, but the observed temporal persistence of "Ca. M. polyxenophila" is most likely explained by the positive effect this symbiont provides to its host.

The ciliate Paramecium bursaria allows budding of symbiotic Chlorella variabilis cells singly from the digestive vacuole membrane into the cytoplasm during algal reinfection

The ciliate Paramecium bursaria harbors several hundred symbiotic Chlorella spp. cells in the cytoplasm. Algal re-endosymbiosis can be artificially induced using alga-removed P. bursaria. During algal re-endosymbiosis, algae ingested into the host digestive vacuoles (DVs) avoid digestion by the host lysosomal enzymes and then escape into the cytoplasm by budding off of the DV membrane. The budded alga-enclosing DV membrane then differentiates into the symbiosome or perialgal vacuole (PV) membrane and is localized beneath the host cell cortex. In this study, we determined whether the PV membrane has the ability to recognize the symbiotic alga singly by eliminating other small microspheres in the same DV. To clarify the accuracy of the budding process, we mixed fluorescent-labeled microspheres of diameter 0.20 µm with isolated symbiotic algae during algal re-endosymbiosis. No fluorescence was observed from the PV membrane, as expected, and the budding DVs that enclosed both undigested and digested algae. Additionally, the algal re-endosymbiosis rate was significantly reduced in the presence of microspheres. These observations showed that the host P. bursaria allowed budding of the algae singly from the membranes of DVs without microspheres and this process required close contact between the DV membrane and the algal cell wall.

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

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

Freshwater sponge hosts and their green algae symbionts: a tractable model to understand intracellular symbiosis

In many freshwater habitats, green algae form intracellular symbioses with a variety of heterotrophic host taxa including several species of freshwater sponge. These sponges perform important ecological roles in their habitats, and the poriferan:green algae partnerships offers unique opportunities to study the evolutionary origins and ecological persistence of endosymbioses. We examined the association between Ephydatia muelleri and its chlorophyte partner to identify features of host cellular and genetic responses to the presence of intracellular algal partners. Chlorella-like green algal symbionts were isolated from field-collected adult E. muelleri tissue harboring algae. The sponge-derived algae were successfully cultured and subsequently used to reinfect aposymbiotic E. muelleri tissue. We used confocal microscopy to follow the fate of the sponge-derived algae after inoculating algae-free E. muelleri grown from gemmules to show temporal patterns of symbiont location within host tissue. We also infected aposymbiotic E. muelleri with sponge-derived algae, and performed RNASeq to study differential expression patterns in the host relative to symbiotic states. We compare and contrast our findings with work in other systems (e.g., endosymbiotic Hydra) to explore possible conserved evolutionary pathways that may lead to stable mutualistic endosymbioses. Our work demonstrates that freshwater sponges offer many tractable qualities to study features of intracellular occupancy and thus meet criteria desired for a model system.

Taxonomy, Ontogenesis and Evolutionary Relationships of the Algae-Bearing Ciliate Bourlandella viridis () comb. nov., With Establishment of a New Genus and New Family (Protista, Ciliophora, Hypotrichia)

Hypotrichs are the most complex and highly differentiated ciliate lineages and have great potential for evolutionary novelties. Problems in hypotrich systematics are mainly due to discordance between the morphological and genetic data (mainly small subunit rRNA gene). Species with morphologies that are characteristic of two or more higher rank taxa are probably a major contributing factor to these conflicts. The present study describes a Chinese population of a poorly known organism with numerous endosymbiotic zoochlorellae, the morphology of which corresponds well with the type population of Limnoholosticha viridis (Kahl, 1932) Li et al., 2017. Newly obtained information shows this species has a zigzag midventral cirral pattern that is diagnostic of urostylids, whereas the dorsal ciliature shares features (presence of dorsomarginal kinety and dorsal kinety 3 fragmentation) that are typical of oxytrichids. Molecular phylogenetic analyses reveal a close relationship with oxytrichids. An integrative approach combining morphological, morphogenetic, phylogenetic and ecological analyses indicates that L. viridis represents a new genus and new family which might be an intermediate form between uorstylids and dorsomarginalians. Thus, Bourlandellidae fam. nov. and Bourlandella gen. nov. are here established. Lastly, we speculate that phenotypic convergence and mixtrophy might confer on the new combination, Bourlandella viridis (Kahl, 1932) comb. nov., the ability to adapt to a wide range of environmental conditions.

Autolysis of Chlorella variabilis in Starving Paramecium bursaria Help the Host Cell Survive Against Starvation Stress

The endosymbiosis between Paramecium bursaria and Chlorella spp. is mutualistic. Symbiotic algae localize beneath the host Paramecium cell cortex compete for their attachment sites with preexisting organelle trichocysts. To examine the relationship between P. bursaria trichocysts and their symbiotic algae, algae-bearing or alga-free P. bursaria were starved for several days and the changes in the number of Chlorella sp. and presence or absence of trichocysts were evaluated. We conducted an indirect immunofluorescence microscopy with an anti-trichocyst monoclonal antibody against P. bursaria cells. Indirect immunofluorescence microscopy demonstrated that under starvation and darkness conditions, the immunofluorescence of trichocysts in alga-free P. bursaria decreased much faster than that in the normal algae-bearing P. bursaria. In the latter case, our observations proposed the possibility that the nutrition obtained from symbiotic algal digestion may promote trichocysts synthesis. This algal digestion mechanism may permit host P. bursaria cells to survive for a longer time under starvation condition. To the best of our knowledge, this may be a new benefit that host P. bursaria gain from harboring symbiotic algae.

Genome plasticity in Paramecium bursaria revealed by population genomics

Background

Ciliates are an ancient and diverse eukaryotic group found in various environments. A unique feature of ciliates is their nuclear dimorphism, by which two types of nuclei, the diploid germline micronucleus (MIC) and polyploidy somatic macronucleus (MAC), are present in the same cytoplasm and serve different functions. During each sexual cycle, ciliates develop a new macronucleus in which newly fused genomes are extensively rearranged to generate functional minichromosomes. Interestingly, each ciliate species seems to have its way of processing genomes, providing a diversity of resources for studying genome plasticity and its regulation. Here, we sequenced and analyzed the macronuclear genome of different strains of Paramecium bursaria, a highly divergent species of the genus Paramecium which can stably establish endosymbioses with green algae.

Results

We assembled a high-quality macronuclear genome of P. bursaria and further refined genome annotation by comparing population genomic data. We identified several species-specific expansions in protein families and gene lineages that are potentially associated with endosymbiosis. Moreover, we observed an intensive chromosome breakage pattern that occurred during or shortly after sexual reproduction and contributed to highly variable gene dosage throughout the genome. However, patterns of copy number variation were highly correlated among genetically divergent strains, suggesting that copy number is adjusted by some regulatory mechanisms or natural selection. Further analysis showed that genes with low copy number variation among populations tended to function in basic cellular pathways, whereas highly variable genes were enriched in environmental response pathways.

Conclusions

We report programmed DNA rearrangements in the P. bursaria macronuclear genome that allow cells to adjust gene copy number globally according to individual gene functions. Our results suggest that large-scale gene copy number variation may represent an ancient mechanism for cells to adapt to different environments.

Supplementary information

The online version contains supplementary material available at 10.1186/s12915-020-00912-2.

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

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

Algal Diversity in Paramecium bursaria: Species Identification, Detection of Choricystis parasitica, and Assessment of the Interaction Specificity

The ‘green’ ciliate Paramecium bursaria lives in mutualistic symbiosis with green algae belonging to the species Chlorella variabilis or Micractinium conductrix. We analysed the diversity of algal endosymbionts and their P. bursaria hosts in nine strains from geographically diverse origins. Therefore, their phylogenies using different molecular markers were inferred. The green paramecia belong to different syngens of P. bursaria. The intracellular algae were assigned to Chl. variabilis, M. conductrix or, surprisingly, Choricystis parasitica. This usually free-living alga co-occurs with M. conductrix in the host’s cytoplasm. Addressing the potential status of Chor. parasitica as second additional endosymbiont, we determined if it is capable of symbiosis establishment and replication within a host cell. Symbiont-free P. bursaria were generated by cycloheximid treatment. Those aposymbiotic P. bursaria were used for experimental infections to investigate the symbiosis specificity not only between P. bursaria and Chor. parasitica but including also Chl. variabilis and M. conductrix. For each algae we observed the uptake and incorporation in individual perialgal vacuoles. These host-symbiont associations are stable since more than five months. Thus, Chor. parasitica and P. bursaria can form an intimate and long-term interaction. This study provides new insights into the diversity of P. bursaria algal symbionts.

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

ParameciumDB (https://paramecium.i2bc.paris-saclay.fr) is a community model organism database for the genome and genetics of the ciliate Paramecium. ParameciumDB development relies on the GMOD (www.gmod.org) toolkit. The ParameciumDB web site has been publicly available since 2006 when the P. tetraurelia somatic genome sequence was released, revealing that a series of whole genome duplications punctuated the evolutionary history of the species. The genome is linked to available genetic data and stocks. ParameciumDB has undergone major changes in its content and website since the last update published in 2011. Genomes from multiple Paramecium species, especially from the P. aurelia complex, are now included in ParameciumDB. A new modern web interface accompanies this transition to a database for the whole Paramecium genus. Gene pages have been enriched with orthology relationships, among the Paramecium species and with a panel of model organisms across the eukaryotic tree. This update also presents expert curation of Paramecium mitochondrial genomes.

The completed macronuclear genome of a model ciliate Tetrahymena thermophila and its application in genome scrambling and copy number analyses

The ciliate Tetrahymena thermophila has been a powerful model system for molecular and cellular biology. However, some investigations have been limited due to the incomplete closure and sequencing of the macronuclear genome assembly, which for many years has been stalled at 1,158 scaffolds, with large sections of unknown sequences (available in Tetrahymena Genome Database, TGD, http://ciliate.org/ ). Here we completed the first chromosome-level Tetrahymena macronuclear genome assembly, with approximately 300× long Single Molecule, Real-Time reads of the wild-type SB210 cells-the reference strain for the initial macronuclear genome sequencing project. All 181 chromosomes were capped with two telomeres and gaps were entirely closed. The completed genome shows significant improvements over the current assembly (TGD 2014) in both chromosome structure and sequence integrity. The majority of previously identified gene models shown in TGD were retained, with the addition of 36 new genes and 883 genes with modified gene models. The new genome and annotation were incorporated into TGD. This new genome allows for pursuit in some underexplored areas that were far more challenging previously; two of them, genome scrambling and chromosomal copy number, were investigated in this study. We expect that the completed macronuclear genome will facilitate many studies in Tetrahymena biology, as well as multiple lines of research in other eukaryotes.

Genome analyses of the new model protist Euplotes vannus focusing on genome rearrangement and resistance to environmental stressors

As a model organism for studies of cell and environmental biology, the free-living and cosmopolitan ciliate Euplotes vannus shows intriguing features like dual genome architecture (i.e., separate germline and somatic nuclei in each cell/organism), "gene-sized" chromosomes, stop codon reassignment, programmed ribosomal frameshifting (PRF) and strong resistance to environmental stressors. However, the molecular mechanisms that account for these remarkable traits remain largely unknown. Here we report a combined analysis of de novo assembled high-quality macronuclear (MAC; i.e., somatic) and partial micronuclear (MIC; i.e., germline) genome sequences for E. vannus, and transcriptome profiling data under varying conditions. The results demonstrate that: (a) the MAC genome contains more than 25,000 complete "gene-sized" nanochromosomes (~85 Mb haploid genome size) with the N50 ~2.7 kb; (b) although there is a high frequency of frameshifting at stop codons UAA and UAG, we did not observe impaired transcript abundance as a result of PRF in this species as has been reported for other euplotids; (c) the sequence motif 5'-TA-3' is conserved at nearly all internally-eliminated sequence (IES) boundaries in the MIC genome, and chromosome breakage sites (CBSs) are duplicated and retained in the MAC genome; (d) by profiling the weighted correlation network of genes in the MAC under different environmental stressors, including nutrient scarcity, extreme temperature, salinity and the presence of ammonia, we identified gene clusters that respond to these external physical or chemical stimulations, and (e) we observed a dramatic increase in HSP70 gene transcription under salinity and chemical stresses but surprisingly, not under temperature changes; we link this temperature-resistance to the evolved loss of temperature stress-sensitive elements in regulatory regions. Together with the genome resources generated in this study, which are available online at Euplotes vannus Genome Database (http://evan.ciliate.org), these data provide molecular evidence for understanding the unique biology of highly adaptable microorganisms.