Enzymatic and chemical mapping of nucleosome distribution in purified micro- and macronuclei of the ciliated model organism, Tetrahymena thermophila

Group-specific functional patterns of mitochondrion-related organelles shed light on their multiple transitions from mitochondria in ciliated protists

Adaptations of ciliates to hypoxic environments have arisen independently several times. Studies on mitochondrion-related organelle (MRO) metabolisms from distinct anaerobic ciliate groups provide evidence for understanding the transitions from mitochondria to MROs within eukaryotes. To deepen our knowledge about the evolutionary patterns of ciliate anaerobiosis, mass-culture and single-cell transcriptomes of two anaerobic species, Metopus laminarius (class Armophorea) and Plagiopyla cf. narasimhamurtii (class Plagiopylea), were sequenced and their MRO metabolic maps were compared. In addition, we carried out comparisons using publicly available predicted MRO proteomes from other ciliate classes (i.e., Armophorea, Litostomatea, Muranotrichea, Oligohymenophorea, Parablepharismea and Plagiopylea). We found that single-cell transcriptomes were similarly comparable to their mass-culture counterparts in predicting MRO metabolic pathways of ciliates. The patterns of the components of the MRO metabolic pathways might be divergent among anaerobic ciliates, even among closely related species. Notably, our findings indicate the existence of group-specific functional relics of electron transport chains (ETCs). Detailed group-specific ETC functional patterns are as follows: full oxidative phosphorylation in Oligohymenophorea and Muranotrichea; only electron-transfer machinery in Armophorea; either of these functional types in Parablepharismea; and ETC functional absence in Litostomatea and Plagiopylea. These findings suggest that adaptation of ciliates to anaerobic conditions is group-specific and has occurred multiple times. Our results also show the potential and the limitations of detecting ciliate MRO proteins using single-cell transcriptomes and improve the understanding of the multiple transitions from mitochondria to MROs within ciliates.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42995-022-00147-w.

Identification and utilization of a mutated 60S ribosomal subunit coding gene as an effective and cost-efficient selection marker for Tetrahymena genetic manipulation

Since the onset of molecular biology, the ciliate Tetrahymena thermophila has been one of the most convenient single-celled model eukaryotes for genetics, biochemistry, and cell biology. Particularly, thanks to the availability of several different selection markers, it is possible to knock out or knock in genes at multiple genetic loci simultaneously in Tetrahymena, which makes it an excellent model ciliate for tackling complex regulatory mechanisms. Despite these selection markers are efficient for genetic manipulation, the costly drugs used for selection have highlighted the urgent demand for an additional cost-efficient and effective selection marker. Here, we found that a mutated 60S ribosomal subunit component, RPL36A, confers T. thermophila with cycloheximide resistance. On top of that, we developed a cycloheximide cassette and explored suitable transformation and selection conditions. Using the new cassette, we obtained both knockout and knock-in strains successfully at a relatively low cost. This study also provided the first evidence that a cycloheximide resistance gene can be engineered as a selection marker to completely delete a gene from the highly-polyploid somatic nucleus in Tetrahymena.

Case Study of the Response of N6-Methyladenine DNA Modification to Environmental Stressors in the Unicellular Eukaryote Tetrahymena thermophila

Rediscovered as a potential epigenetic mark, N⁶-methyladenine DNA modification (6mA) was recently reported to be sensitive to environmental stressors in several multicellular eukaryotes. As 6mA distribution and function differ significantly in multicellular and unicellular organisms, whether and how 6mA in unicellular eukaryotes responds to environmental stress remains elusive. Here, we characterized the dynamic changes of 6mA under starvation in the unicellular model organism Tetrahymena thermophila. Single-molecule, real-time (SMRT) sequencing reveals that DNA 6mA levels in starved cells are significantly reduced, especially symmetric 6mA, compared to those in vegetatively growing cells. Despite a global 6mA reduction, the fraction of asymmetric 6mA with a high methylation level was increased, which might be the driving force for stronger nucleosome positioning in starved cells. Starvation affects expression of many metabolism-related genes, the expression level change of which is associated with the amount of 6mA change, thereby linking 6mA with global transcription and starvation adaptation. The reduction of symmetric 6mA and the increase of asymmetric 6mA coincide with the downregulation of AMT1 and upregulation of AMT2 and AMT5, which are supposedly the MT-A70 methyltransferases required for symmetric and asymmetric 6mA, respectively. These results demonstrated that a regulated 6mA response to environmental cues is evolutionarily conserved in eukaryotes.

IMPORTANCE Increasing evidence indicated that 6mA could respond to environmental stressors in multicellular eukaryotes. As 6mA distribution and function differ significantly in multicellular and unicellular organisms, whether and how 6mA in unicellular eukaryotes responds to environmental stress remains elusive. In the present work, we characterized the dynamic changes of 6mA under starvation in the unicellular model organism Tetrahymena thermophila. Our results provide insights into how Tetrahymena fine-tunes its 6mA level and composition upon starvation, suggesting that a regulated 6mA response to environmental cues is evolutionarily conserved in eukaryotes.

New contribution to epigenetic studies: Isolation of micronuclei with high purity and DNA integrity in the model ciliated protist, Tetrahymena thermophila

The ciliated protist Tetrahymena thermophila is a well-known model organism with typical nuclear dimorphism containing a somatic macronucleus (MAC) and a germline micronucleus (MIC). The presence in the same cell compartment of two nuclei with distinctly different structural and functional properties provides an ideal model system to explore mechanisms of genome maintenance. Although methods for the isolation of MIC have been available for many years, cross-contamination and DNA degradation remain unresolved. Here, we describe a reliable and quick method to isolate MIC with high purity and DNA integrity in T. thermophila. Different factors are examined to optimize the MIC purification. The MAC contamination ratio in purified MIC is about 0.19% and DNA integrity of purified MIC is maintained. We also establish a more accurate method to detect the contamination rate of nuclei including microscopic observation and PCR detection. This study will facilitate further epigenetic research in Tetrahymena.

An Optimized and Versatile Counter-Flow Centrifugal Elutriation Workflow to Obtain Synchronized Eukaryotic Cells

Cell synchronization is a powerful tool to understand cell cycle events and its regulatory mechanisms. Counter-flow centrifugal elutriation (CCE) is a more generally desirable method to synchronize cells because it does not significantly alter cell behavior and/or cell cycle progression, however, adjusting specific parameters in a cell type/equipment-dependent manner can be challenging. In this paper, we used the unicellular eukaryotic model organism, Tetrahymena thermophila as a testing system for optimizing CCE workflow. Firstly, flow cytometry conditions were identified that reduced nuclei adhesion and improved the assessment of cell cycle stage. We then systematically examined how to achieve the optimal conditions for three critical factors affecting the outcome of CCE, including loading flow rate, collection flow rate and collection volume. Using our optimized workflow, we obtained a large population of highly synchronous G1-phase Tetrahymena as measured by 5-ethynyl-2'-deoxyuridine (EdU) incorporation into nascent DNA strands, bulk DNA content changes by flow cytometry, and cell cycle progression by light microscopy. This detailed protocol can be easily adapted to synchronize other eukaryotic cells.

Further insights into the phylogeny of peniculid ciliates (Ciliophora, Oligohymenophorea) based on multigene data

Peniculids comprise a large order of ciliated protists in Class Oligohymenophorea having many unresolved evolutionary relationships. Herein, we report 27 new sequences, including 18S rRNA, ITS1-5.8S- ITS2 rRNA, 28S rRNA and the mitochondrial cox1 genes of eight peniculids. We conducted phylogenetic analyses based on each these markers and on a four-gene concatenated data set (18S rRNA, ITS1-5.8S- ITS2 rRNA, 28S rRNA, and cox1 gene). The main findings are: 1) subclass Peniculia and family Parameciidae are monophyletic, with genus Frontonia remaining non-monophyletic; 2) Urocentrids have traditionally been regarded as a family, multi-gene analyses support the rank of Urocentrida and consistently recovers this order as sister to Peniculida, and Urocentrida and Peniculida comprise subclass Peniculia in agreement with Lynn's (2008) classification; 3) discrepancies between multiple-gene phylogenies, and conflicts with morphologic data regarding genus Frontonia necessitate expansion and revision of species diagnoses and we propose consideration of Group III of Frontonia (including F. didieri, F. ocularis, F. anatolica, F. pusilla and F. elegans) as incertae sedis in Peniculida; 4) multi-gene analyses of Parameciidae support five previously established subgenera. Paramecium buetschlii is placed in subgenus Chloroparamecium, and P. chlorelligerum into subgenus Viridoparamecium.

A distinct class of eukaryotic MT-A70 methyltransferases maintain symmetric DNA N6-adenine methylation at the ApT dinucleotides as an epigenetic mark associated with transcription

Rediscovered as a potential eukaryotic epigenetic mark, DNA N⁶-adenine methylation (6mA) varies across species in abundance and its relationships with transcription. Here we characterize AMT1—representing a distinct MT-A70 family methyltransferase—in the ciliate Tetrahymena thermophila. AMT1 loss-of-function leads to severe defects in growth and development. Single Molecule, Real-Time (SMRT) sequencing reveals that AMT1 is required for the bulk of 6mA and all symmetric methylation at the ApT dinucleotides. The detection of hemi-methylated ApT sites suggests a semi-conservative mechanism for maintaining symmetric methylation. AMT1 affects expression of many genes; in particular, RAB46, encoding a Rab family GTPase involved in contractile vacuole function, is likely a direct target. The distribution of 6mA resembles H3K4 methylation and H2A.Z, two conserved epigenetic marks associated with RNA polymerase II transcription. Furthermore, strong 6mA and nucleosome positioning in wild-type cells is attenuated in ΔAMT1 cells. Our results support that AMT1-catalyzed 6mA is an integral part of the transcription-associated epigenetic landscape. AMT1 homologues are generally found in protists and basal fungi featuring ApT hyper-methylation associated with transcription, which are missing in animals, plants, and true fungi. This dichotomy of 6mA functions and the underlying molecular mechanisms may have implications in eukaryotic diversification.

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.

Cyclin Cyc2p is required for micronuclear bouquet formation in Tetrahymena thermophila

Meiotic bouquet formation (known as crescent formation in Tetrahymena thermophila) is indispensable for homologous pairing and recombination, but the regulatory mechanism of bouquet formation remains largely unknown. As a conjugation specific cyclin gene, CYC2 knockout mutants failed to form an elongated crescent structure and aborted meiosis progress in T. thermophila. γ-H2A.X staining revealed fewer micronuclear DNA double-strand breaks (DSBs) in cyc2Δ cells than in wild-type cells. Furthermore, cyc2Δ cells still failed to form a crescent structure even though DSBs were induced by exogenous agents, indicating that a lack of DSBs was not completely responsible for failure to enter the crescent stage. Tubulin staining showed that impaired perinuclear microtubule structure may contribute to the blockage in micronuclear elongation. At the same time, expression of microtubule-associated kinesin genes, KIN11 and KIN141, was significantly downregulated in cyc2Δ cells. Moreover, micronuclear specific accumulation of heterochromatin marker trimethylated H3K23 abnormally increased in the cyc2Δ mutants. Together, these results show that cyclin Cyc2p is required for micronuclear bouquet formation via controlling microtubule-directed nuclear elongation in Tetrahymena.

Further analyses of variation of ribosome DNA copy number and polymorphism in ciliates provide insights relevant to studies of both molecular ecology and phylogeny

Sequence-based approaches, such as analyses of ribosome DNA (rDNA) clone libraries and high-throughput amplicon sequencing, have been used extensively to infer evolutionary relationships and elucidate the biodiversity in microbial communities. However, recent studies demonstrate both rDNA copy number variation and intra-individual (intra-genomic) sequence variation in many organisms, which challenges the application of the rDNA-based surveys. In ciliates, an ecologically important clade of microbial eukaryotes, rDNA copy number and sequence variation are rarely studied. In the present study, we estimate the intraindividual small subunit rDNA (SSU rDNA) copy number and sequence variation in a wide range of taxa covering nine classes and 18 orders of the phylum Ciliophora. Our studies reveal that: (i) intra-individual sequence variation of SSU rDNA is ubiquitous in all groups of ciliates detected and the polymorphic level varies among taxa; (ii) there is a most common version of SSU rDNA sequence in each cell that is highly predominant and may represent the germline micronuclear template; (iii) compared with the most common version, other variant sequences differ in only 1-3 nucleotides, likely generated during macronuclear (somatic) amplification; (iv) the intra-cell sequence variation is unlikely to impact phylogenetic analyses; (v) the rDNA copy number in ciliates is highly variable, ranging from 10³ to 10⁶, with the highest record in Stentor roeselii. Overall, these analyses indicate the need for careful consideration of SSU rDNA variation in analyses of the role of ciliates in ecosystems.

Redescription of a Hymenostome Ciliate, Tetrahymena setosa (Protozoa, Ciliophora) Notes on its Molecular Phylogeny

In recent years, Tetrahymena species have been used as model organisms for research in a wide range of fields, highlighting the need for a fuller understanding of the taxonomy of this group. It is in this context that this paper uses living observation and silver staining methods to investigate the morphology and infraciliature of one Tetrahymena species, T. setosa (Schewiakoff 1892 Verh. Naturh. Med. Ver. Heidelb., 4:544) McCoy (1975) Acta Protozool., 14:253; the senior subjective synonym of T. setifera Holz and Corliss (1956) J. Protozool., 3:112; isolated from a freshwater pond in Harbin, north-eastern China. This organism can be distinguished from other described Tetrahymena species mainly by its single caudal cilium, which is about twice the length of the somatic ciliature. While the Harbin isolate appears similar to the population described by Holz and Corliss (1956) J. Protozool., 3:112, an improved diagnosis for T. setosa is given based on the previous descriptions and the Harbin population. In summary, this species can be recognized mainly by the combination of the following characters: body in vivo approximately 40 μm × 25 μm, 21-26 somatic kineties, one to four contractile vacuole pores associated with meridians 6-11 and a single caudal cilium. The small subunit ribosomal (SSU) rRNA gene and the cox1 gene sequences of Harbin population are also characterized in order to corroborate that the isolated species branches in phylogenetic trees as a T. setosa species. The phylogenetic analysis also indicated that sequences of populations of Tetrahymena species should be published with detailed morphological identifications.

New contribution to the species-rich genus Euplotes: Morphology, ontogeny and systematic position of two species (Ciliophora; Euplotia)

The morphology, ontogeny and phylogeny of two Euplotes species, E. estuarinus sp. nov. and a population of E. platystoma Dragesco and Dragesco-Kernéis, 1986, both collected from tropical brackish waters in south China, were investigated based on living morphology, ciliary pattern and molecular data. Euplotes estuarinus sp. nov. is small (about 60 × 40 μm in vivo), has a dargyrome of the double-eurystomus type, and the transverse cirri are arranged in two groups, with two left and three right ones. The original description of the poorly known species, E. platystoma, is brief, and the species was never investigated using live observation and molecular methods Hence, we provided a detailed redescription. Some stages of their morphogenesis were observed which proceed in the same pattern as in their congeners. The new species E. estuarinus sp. nov. clusters with E. curdsi, differing only by 1 bp in their SSU rRNA gene sequences, which is likely due to the recent speciation event and the limited resolution of the SSU rRNA gene at species level in this group as the two species are clearly morphologically distinct.

N6-adenine DNA methylation is associated with the linker DNA of H2A.Z-containing well-positioned nucleosomes in Pol II-transcribed genes in Tetrahymena

DNA N6-methyladenine (6mA) is newly rediscovered as a potential epigenetic mark across a more diverse range of eukaryotes than previously realized. As a unicellular model organism, Tetrahymena thermophila is among the first eukaryotes reported to contain 6mA modification. However, lack of comprehensive information about 6mA distribution hinders further investigations into its function and regulatory mechanism. In this study, we provide the first genome-wide, base pair-resolution map of 6mA in Tetrahymena by applying single-molecule real-time (SMRT) sequencing. We provide evidence that 6mA occurs mostly in the AT motif of the linker DNA regions. More strikingly, these linker DNA regions with 6mA are usually flanked by well-positioned nucleosomes and/or H2A.Z-containing nucleosomes. We also find that 6mA is exclusively associated with RNA polymerase II (Pol II)-transcribed genes, but is not an unambiguous mark for active transcription. These results support that 6mA is an integral part of the chromatin landscape shaped by adenosine triphosphate (ATP)-dependent chromatin remodeling and transcription.

Parallel mapping with site-directed hydroxyl radicals and micrococcal nuclease reveals structural features of positioned nucleosomes in vivo

Micrococcal nuclease (MNase) has been widely used for analyses of nucleosome locations in many organisms. However, due to its sequence preference, the interpretations of the positions and occupancies of nucleosomes using MNase have remained controversial. Next-generation sequencing (NGS) has also been utilized for analyses of MNase-digests, but some technical biases are commonly present in the NGS experiments. Here, we established a gel-based method to map nucleosome positions in Saccharomyces cerevisiae, using isolated nuclei as the substrate for the histone H4 S47C-site-directed chemical cleavage in parallel with MNase digestion. The parallel mapping allowed us to compare the chemically and enzymatically cleaved sites by indirect end-labeling and primer extension mapping, and thus we could determine the nucleosome positions and the sizes of the nucleosome-free regions (or nucleosome-depleted regions) more accurately, as compared to nucleosome mapping by MNase alone. The analysis also revealed that the structural features of the nucleosomes flanked by the nucleosome-free region were different from those within regularly arrayed nucleosomes, showing that the structures and dynamics of individual nucleosomes strongly depend on their locations. Moreover, we demonstrated that the parallel mapping results were generally consistent with the previous genome-wide chemical mapping and MNase-Seq results. Thus, the gel-based parallel mapping will be useful for the analysis of a specific locus under various conditions.

A Contribution to the Morphology and Phylogeny of Chlamydodon, with Three New Species from China (Ciliophora, Cyrtophoria)

Three new cyrtophorian ciliates isolated from coastal areas of China were described based on morphological and genetic data. The Chlamydodon mnemosyne-like species Chlamydodon similis sp. n. differs from its congeners mainly by its number of somatic kineties. Chlamydodon oligochaetus sp. n. is distinguished from its congeners mainly by having fewer somatic kineties, and/or an elongated body shape. Chlamydodon crassidens sp. n. is characterized mainly by an inverted triangular body shape, a posteriorly interrupted cross-striated band (5-6 μm wide), and a large cytostome. Moreover, we provided small-subunit (SSU) rDNA sequences of C. similis sp. n. and C. oligochaetus sp. n. Maximum likelihood (ML) and Bayesian inference (BI) consistently placed C. similis sp. n. as a sister to C. paramnemosyne, but showed different branching position of C. oligochaetus sp. n., which may be due to a low taxon sampling in the Chlamydodontidae and/or an insufficient resolution of the marker gene at species level.