Genome plasticity in Paramecium bursaria revealed by population genomics

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.

Characterization of Crystals in Ciliate Paramecium bursaria Harboring Endosymbiotic Chlorella variabilis

Protists, including ciliates retain crystals in their cytoplasm. However, their functions and properties remain unclear. To comparatively analyze the crystals of Paramecium bursaria, a ciliate, associated with and without the endosymbiotic Chlorella variabilis, we investigated the isolated crystals using a light microscope and analyzed their length and solubility. A negligible number of crystals was found in P. bursaria cells harboring symbiotic algae. The average crystal length in alga-free and algae-reduced cells was about 6.8 μm and 14.4 μm, respectively. The crystals of alga-free cells were spherical, whereas those of algae-reduced cells were angular in shape. The crystals of alga-free cells immediately dissolved in acids and bases, but not in water or organic solvents, and were stable at - 20 °C for more than 3 weeks. This study, for the first time, reveals that the characteristics of crystals present in the cytoplasm of P. bursaria vary greatly depending on the amount of symbiotic algae.

Mutualism on the edge: Understanding the Paramecium-Chlorella symbiosis

Exploring the mechanisms that underpin symbiosis requires an understanding of how these complex interactions are maintained in diverse model systems. The ciliate protist, Paramecium bursaria, offers a valuable insight into how emergent endosymbiotic interactions have evolved.

Comparative Analyses of the Symbiotic Associations of the Host Paramecium bursaria with Free-Living and Native Symbiotic Species of Chlorella

Each symbiotic Chlorella variabilis associated with the ciliate Paramecium bursaria is enclosed in a symbiosome called the perialgal vacuole. Various potential symbionts, such as bacteria, yeasts, other algae, and free-living Chlorella spp., can infect P. bursaria. However, the detailed infection process of each of them in algae-free P. bursaria is unknown. Here, we aimed to elucidate the difference of the infection process between the free-living C. sorokiniana strain NIES-2169 and native symbiotic C. variabilis strain 1N. We investigated the fate of ingested algae using algae-free P. bursaria exposed separately to three types of algal inocula: NIES-2169 only, 1N only, or a mixture of NIES-2169 and 1N. We found that (1) only one algal species, preferably the native one, was retained in host cells, indicating a type of host compatibility and (2) the algal localization style beneath the host cell cortex varied between different Chlorella spp. showing various levels of host compatibilities, which was prospectively attributable to the difference in the formation of the perialgal vacuole membrane.

The macronuclear genomic landscape within Tetrahymena thermophila

The extent of intraspecific genomic variation is key to understanding species evolutionary history, including recent adaptive shifts. Intraspecific genomic variation remains poorly explored in eukaryotic micro-organisms, especially in the nuclear dimorphic ciliates, despite their fundamental role as laboratory model systems and their ecological importance in many ecosystems. We sequenced the macronuclear genome of 22 laboratory strains of the oligohymenophoran Tetrahymena thermophila, a model species in both cellular biology and evolutionary ecology. We explored polymorphisms at the junctions of programmed eliminated sequences, and reveal their utility to barcode very closely related cells. As for other species of the genus Tetrahymena, we confirm micronuclear centromeres as gene diversification centres in T. thermophila, but also reveal a two-speed evolution in these regions. In the rest of the genome, we highlight recent diversification of genes coding for extracellular proteins and cell adhesion. We discuss all these findings in relation to this ciliate's ecology and cellular characteristics.

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

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

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.

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.

Origins of genome-editing excisases as illuminated by the somatic genome of the ciliate Blepharisma

Massive DNA excision occurs regularly in ciliates, ubiquitous microbial eukaryotes with somatic and germline nuclei in the same cell. Tens of thousands of internally eliminated sequences (IESs) scattered throughout the ciliate germline genome are deleted during the development of the streamlined somatic genome. The genus Blepharisma represents one of the two high-level ciliate clades (subphylum Postciliodesmatophora) and, unusually, has dual pathways of somatic nuclear and genome development. This makes it ideal for investigating the functioning and evolution of these processes. Here we report the somatic genome assembly of Blepharisma stoltei strain ATCC 30299 (41 Mbp), arranged as numerous telomere-capped minichromosomal isoforms. This genome encodes eight PiggyBac transposase homologs no longer harbored by transposons. All appear subject to purifying selection, but just one, the putative IES excisase, has a complete catalytic triad. We hypothesize that PiggyBac homologs were ancestral excisases that enabled the evolution of extensive natural genome editing.

Role of host ciliate Paramecium bursaria mitochondria and trichocysts for symbiotic Chlorella variabilis attachment beneath the host cell cortex

Symbiotic Chlorella variabilis is encased in the perialgal vacuole (PV) membrane of ciliate Paramecium bursaria. The PV membrane is stably anchored below the host cell cortex by adhesion to host mitochondria. Host trichocysts, which are defensive organelles against predators, are present in the mitochondria and PV membrane vicinity. The mechanism by which PV attaches beneath the host cell cortex remains unknown. When P. bursaria is centrifuged at high speed, the symbiotic algae are displaced from the host cell cortex and concentrate at the posterior end. When centrifugation is stopped, the dislocated algae reattach beneath the host cell cortex with fast cytoplasmic streaming. The densities of mitochondria and trichocysts before and after centrifugation were compared using indirect immunofluorescence microscopy with monoclonal antibodies. Almost all trichocysts were shed by high-speed centrifugation, but dislocated algae could reattach even in the absence of trichocysts. In contrast, host mitochondria were unaffected in localization and number, and the dislocated algae also reattached. These findings suggest trichocysts are unnecessary for algal relocalization and that mitochondria are colocalized with the algae. However, many mitochondria were also present in the cell's anterior region without symbiotic algae. Therefore, not all areas with mitochondria contained algae, but there was an algal localization bias within the host cell.

Paramecium bursaria as a Potential Tool for Evaluation of Microplastics Toxicity

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

Species-Specific Duplication of Surface Antigen Genes in Paramecium

Paramecium is a free-living ciliate that undergoes antigenic variation and still the functions of these variable surface antigen coats in this non-pathogenic ciliate remain elusive. Only a few surface antigen genes have been described, mainly in the two model species P. tetraurelia strain 51 and P. primaurelia strain 156. Given the lack of suitable sequence data to allow for phylogenetics and deeper sequence comparisons, we screened the genomes of six different Paramecium species for serotype genes and isolated 548 candidates. Our approach identified the subfamilies of the isogenes of individual serotypes that were mostly represented by intrachromosomal gene duplicates. These showed different duplication levels, and chromosome synteny suggested rather young duplication events after the emergence of the P. aurelia species complex, indicating a rapid evolution of surface antigen genes. We were able to identify the different subfamilies of the surface antigen genes with internal tandem repeats, which showed consensus motifs across species. The individual isogene families showed additional consensus motifs, indicating that the selection pressure holds individual amino acids constant in these repeats. This may be a hint of the receptor function of these antigens rather than a presentation of random epitopes, generating the variability of these surface molecules.

Macronuclear development in ciliates, with a focus on nuclear architecture

Ciliates are defined by the presence of dimorphic nuclei as they have both a somatic macronucleus and germline micronucleus within each individual cell. The size and structure of both germline micronuclei and somatic macronuclei varies tremendously among ciliates. Except just after conjugation (i.e. the nuclear exchange in sexual cycle), the germline micronucleus is transcriptionally-inactive and contains canonical chromosomes that will be inherited between generations. In contrast, the transcriptionally-active macronucleus contains chromosomes that vary in size in different classes of ciliates, with some lineages having extensively-fragmented gene-sized somatic chromosomes while others contain longer multigene chromosomes. Here, we describe the variation in somatic macronuclear architecture in lineages sampled across the ciliate tree of life, specifically focusing on lineages with extensively fragmented chromosomes (e.g. the classes Phyllopharyngea and Spirotrichea). Further, we synthesize information from the literature on the development of ciliate macronuclei, focusing on changes in nuclear architecture throughout life cycles. These data highlight the tremendous diversity among ciliate nuclear cycles, extend our understanding of patterns of genome evolution, and provide insight into different germline and somatic nuclear features (e.g. nuclear structure and development) among eukaryotes.

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.

Phenotypic plasticity through disposable genetic adaptation in ciliates

Ciliates have an extraordinary genetic system in which each cell harbors two distinct kinds of nucleus, a transcriptionally active somatic nucleus and a quiescent germline nucleus. The latter undergoes classical, heritable genetic adaptation, while adaptation of the somatic nucleus is only short-term and thus disposable. The ecological and evolutionary relevance of this nuclear dimorphism have never been well formalized, which is surprising given the long history of using ciliates such as Tetrahymena and Paramecium as model organisms. We present a novel, alternative explanation for ciliate nuclear dimorphism which, we argue, should be considered an instrument of phenotypic plasticity by somatic selection on the level of the ciliate clone, as if it were a diffuse multicellular organism. This viewpoint helps to put some enigmatic aspects of ciliate biology into perspective and presents the diversity of ciliates as a large natural experiment that we can exploit to study phenotypic plasticity and organismality.