RNAi-dependent Polycomb repression controls transposable elements in Tetrahymena

Conserved functions of chromatin regulators in basal Archaeplastida

Chromatin is a dynamic network that regulates genome organization and gene expression. Different types of chromatin regulators are highly conserved among Archaeplastida, including unicellular algae, while some chromatin genes are only present in land plant genomes. Here, we review recent advances in understanding the function of conserved chromatin factors in basal land plants and algae. We focus on the role of Polycomb-group genes which mediate H3K27me3-based silencing and play a role in balancing gene dosage and regulating haploid-to-diploid transitions by tissue-specific repression of the transcription factors KNOX and BELL in many representatives of the green lineage. Moreover, H3K27me3 predominantly occupies repetitive elements which can lead to their silencing in a unicellular alga and basal land plants, while it covers mostly protein-coding genes in higher land plants. In addition, we discuss the role of nuclear matrix constituent proteins as putative functional lamin analogs that are highly conserved among land plants and might have an ancestral function in stress response regulation. In summary, our review highlights the importance of studying chromatin regulation in a wide range of organisms in the Archaeplastida.

Chromatin profiling identifies putative dual roles for H3K27me3 in regulating transposons and cell type-specific genes in choanoflagellates

Gene expression is tightly controlled during animal development to allow the formation of specialized cell types. Our understanding of how animals evolved this exquisite regulatory control remains elusive, but evidence suggests that changes in chromatin-based mechanisms may have contributed. To investigate this possibility, here we examine chromatin-based gene regulatory features in the closest relatives of animals, choanoflagellates. Using Salpingoeca rosetta as a model system, we examined chromatin accessibility and histone modifications at the genome scale and compared these features to gene expression. We first observed that accessible regions of chromatin are primarily associated with gene promoters and found no evidence of distal gene regulatory elements resembling the enhancers that animals deploy to regulate developmental gene expression. Remarkably, a histone modification deposited by polycomb repressive complex 2, histone H3 lysine 27 trimethylation (H3K27me3), appeared to function similarly in S. rosetta to its role in animals, because this modification decorated genes with cell type-specific expression. Additionally, H3K27me3 marked transposons, retaining what appears to be an ancestral role in regulating these elements. We further uncovered a putative new bivalent chromatin state at cell type-specific genes that consists of H3K27me3 and histone H3 lysine 4 mono-methylation (H3K4me1). Together, our discoveries support the scenario that gene-associated histone modification states that underpin development emerged before the evolution of animal multicellularity.

The genome-wide meiotic recombination landscape in ciliates and its implications for crossover regulation and genome evolution

Meiotic recombination is essential for sexual reproduction and its regulation has been extensively studied in many taxa. However, genome-wide recombination landscape has not been reported in ciliates and it remains unknown how it is affected by the unique features of ciliates: the synaptonemal complex (SC)-independent meiosis and the nuclear dimorphism. Here, we show the recombination landscape in the model ciliate Tetrahymena thermophila by analyzing single-nucleotide polymorphism datasets from 38 hybrid progeny. We detect 1021 crossover (CO) events (35.8 per meiosis), corresponding to an overall CO rate of 9.9 cM/Mb. However, gene conversion by non-crossover is rare (1.03 per meiosis) and not biased towards G or C alleles. Consistent with the reported roles of SC in CO interference, we find no obvious sign of CO interference. CO tends to occur within germ-soma common genomic regions and many of the 44 identified CO hotspots localize at the centromeric or subtelomeric regions. Gene ontology analyses show that CO hotspots are strongly associated with genes responding to environmental changes. We discuss these results with respect to how nuclear dimorphism has potentially driven the formation of the observed recombination landscape to facilitate environmental adaptation and the sharing of machinery among meiotic and somatic recombination.

Small-RNA-guided histone modifications and somatic genome elimination in ciliates

Transposable elements and other repeats are repressed by small-RNA-guided histone modifications in fungi, plants and animals. The specificity of silencing is achieved through base-pairing of small RNAs corresponding to the these genomic loci to nascent noncoding RNAs, which allows the recruitment of histone methyltransferases that methylate histone H3 on lysine 9. Self-reinforcing feedback loops enhance small RNA production and ensure robust and heritable repression. In the unicellular ciliate Paramecium tetraurelia, small-RNA-guided histone modifications lead to the elimination of transposable elements and their remnants, a definitive form of repression. In this organism, germline and somatic functions are separated within two types of nuclei with different genomes. At each sexual cycle, development of the somatic genome is accompanied by the reproducible removal of approximately a third of the germline genome. Instead of recruiting a H3K9 methyltransferase, small RNAs corresponding to eliminated sequences tether Polycomb Repressive Complex 2, which in ciliates has the unique property of catalyzing both lysine 9 and lysine 27 trimethylation of histone H3. These histone modifications that are crucial for the elimination of transposable elements are thought to guide the endonuclease complex, which triggers double-strand breaks at these specific genomic loci. The comparison between ciliates and other eukaryotes underscores the importance of investigating small-RNAs-directed chromatin silencing in a diverse range of organisms. This article is categorized under: Regulatory RNAs/RNAi/Riboswitches > RNAi: Mechanisms of Action.

Impact of Germline Depletion of Bonus on Chromatin State in Drosophila Ovaries

Gene expression is controlled via complex regulatory mechanisms involving transcription factors, chromatin modifications, and chromatin regulatory factors. Histone modifications, such as H3K27me3, H3K9ac, and H3K27ac, play an important role in controlling chromatin accessibility and transcriptional output. In vertebrates, the Transcriptional Intermediary Factor 1 (TIF1) family of proteins play essential roles in transcription, cell differentiation, DNA repair, and mitosis. Our study focused on Bonus, the sole member of the TIF1 family in Drosophila, to investigate its role in organizing epigenetic modifications. Our findings demonstrated that depleting Bonus in ovaries leads to a mild reduction in the H3K27me3 level over transposon regions and alters the distribution of active H3K9ac marks on specific protein-coding genes. Additionally, through mass spectrometry analysis, we identified novel interacting partners of Bonus in ovaries, such as PolQ, providing a comprehensive understanding of the associated molecular pathways. Furthermore, our research revealed Bonus's interactions with the Polycomb Repressive Complex 2 and its co-purification with select histone acetyltransferases, shedding light on the underlying mechanisms behind these changes in chromatin modifications.

A practical reference for studying meiosis in the model ciliate Tetrahymena thermophila

Meiosis is a critical cell division program that produces haploid gametes for sexual reproduction. Abnormalities in meiosis are often causes of infertility and birth defects (e.g., Down syndrome). Most organisms use a highly specialized zipper-like protein complex, the synaptonemal complex (SC), to guide and stabilize pairing of homologous chromosomes in meiosis. Although the SC is critical for meiosis in many eukaryotes, there are organisms that perform meiosis without a functional SC. However, such SC-less meiosis is poorly characterized. To understand the features of SC-less meiosis and its adaptive significance, the ciliated protozoan Tetrahymena was selected as a model. Meiosis research in Tetrahymena has revealed intriguing aspects of the regulatory programs utilized in its SC-less meiosis, yet additional efforts are needed for obtaining an in-depth comprehension of mechanisms that are associated with the absence of SC. Here, aiming at promoting a wider application of Tetrahymena for meiosis research, we introduce basic concepts and core techniques for studying meiosis in Tetrahymena and then suggest future directions for expanding the current Tetrahymena meiosis research toolbox. These methodologies could be adopted for dissecting meiosis in poorly characterized ciliates that might reveal novel features. Such data will hopefully provide insights into the function of the SC and the evolution of meiosis from a unique perspective.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42995-022-00149-8.

A small RNA-guided PRC2 complex eliminates DNA as an extreme form of transposon silencing

In animal germlines, transposons are silenced at the transcriptional or post-transcriptional level to prevent deleterious expression. Ciliates employ a more direct approach by physically eliminating transposons from their soma, utilizing piRNAs to recognize transposons and imprecisely excise them. Ancient, mutated transposons often do not require piRNAs and are precisely eliminated. Here, we characterize the Polycomb Repressive Complex 2 (PRC2) in Paramecium and demonstrate its involvement in the removal of transposons and transposon-derived DNA. Our results reveal a striking difference between the elimination of new and ancient transposons at the chromatin level and show that the complex may be guided by Piwi-bound small RNAs (sRNAs). We propose that imprecise elimination in ciliates originates from an ancient transposon silencing mechanism, much like in plants and metazoans, through sRNAs, repressive methylation marks, and heterochromatin formation. However, it is taken a step further by eliminating DNA as an extreme form of transposon silencing.

Comparative Genome Analysis Reveals Cis-Regulatory Elements on Gene-Sized Chromosomes of Ciliated Protists

Gene-sized chromosomes are a distinct feature of the macronuclear genome in ciliated protists known as spirotrichs. These nanochromosomes are often only several kilobase pairs long and contain a coding region for a single gene. However, the ways in which transcription is regulated on nanochromosomes is still largely unknown. Here, we generated macronuclear genome assemblies for two species of Pseudokeronopsis ciliates to better understand transcription regulation on gene-sized chromosomes. We searched within the short subtelomeric regions for potential cis-regulatory elements and identified distinct AT-rich sequences conserved in both species, at both the 5’ and 3’ end of each gene. We further acquired transcriptomic data for these species, which showed the 5’ cis-regulatory element is associated with active gene expression. Gene family evolution analysis suggests nanochromosomes in spirotrichs may originated approximately 900 million years ago. Together our comparative genomic analyses reveal novel insights into the biological roles of cis-regulatory elements on gene-sized chromosomes.

Bacterial N4-methylcytosine as an epigenetic mark in eukaryotic DNA

DNA modifications are used to regulate gene expression and defend against invading genetic elements. In eukaryotes, modifications predominantly involve C5-methylcytosine (5mC) and occasionally N6-methyladenine (6mA), while bacteria frequently use N4-methylcytosine (4mC) in addition to 5mC and 6mA. Here we report that 4mC can serve as an epigenetic mark in eukaryotes. Bdelloid rotifers, tiny freshwater invertebrates with transposon-poor genomes rich in foreign genes, lack canonical eukaryotic C5-methyltransferases for 5mC addition, but encode an amino-methyltransferase, N4CMT, captured from bacteria >60 Mya. N4CMT deposits 4mC at active transposons and certain tandem repeats, and fusion to a chromodomain shapes its “histone-read-DNA-write” architecture recognizing silent chromatin marks. Furthermore, amplification of SETDB1 H3K9me3 histone methyltransferases yields variants preferentially binding 4mC-DNA, suggesting “DNA-read-histone-write” partnership to maintain chromatin-based silencing. Our results show how non-native DNA methyl groups can reshape epigenetic systems to silence transposons and demonstrate the potential of horizontal gene transfer to drive regulatory innovation in eukaryotes.

Eukaryotic DNA can be methylated as 5-methylcytosine and N6-methyladenine, but whether other forms of DNA methylation occur has been controversial. Here the authors show that a bacterial DNA methyltransferase was acquired >60 Mya in bdelloid rotifers that catalyzes N4-methylcytosine addition and is involved in suppression of transposon proliferation.

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.

Polycomb Repressive Complex 2 in Eukaryotes-An Evolutionary Perspective

Polycomb repressive complex 2 (PRC2) represents a group of evolutionarily conserved multi-subunit complexes that repress gene transcription by introducing trimethylation of lysine 27 on histone 3 (H3K27me3). PRC2 activity is of key importance for cell identity specification and developmental phase transitions in animals and plants. The composition, biochemistry, and developmental function of PRC2 in animal and flowering plant model species are relatively well described. Recent evidence demonstrates the presence of PRC2 complexes in various eukaryotic supergroups, suggesting conservation of the complex and its function. Here, we provide an overview of the current understanding of PRC2-mediated repression in different representatives of eukaryotic supergroups with a focus on the green lineage. By comparison of PRC2 in different eukaryotes, we highlight the possible common and diverged features suggesting evolutionary implications and outline emerging questions and directions for future research of polycomb repression and its evolution.

Silencing and anti-silencing mechanisms that shape the epigenome in plants

Epigenome information mediates genome function and maintenance by regulating gene expression and chromatin organization. Because the epigenome pattern can change in response to internal and external environments, it may underlie an adaptive genome response that modulates phenotypes during development and in changing environments. Here I summarize recent progress in our understanding of how epigenome patterns are shaped and modulated by concerted actions of silencing and anti-silencing factors mainly in Arabidopsis thaliana. I discuss the dynamic nature of epigenome regulation, which is realized by cooperation and counteraction among silencing and anti-silencing factors, and how the dynamic epigenome mediates robust and plastic responses of plants to fluctuating environments.

Combinations of maternal-specific repressive epigenetic marks in the endosperm control seed dormancy

Polycomb Repressive Complex 2 (PRC2)-mediated trimethylation of histone H3 on lysine 27 (H3K27me3) and methylation of histone 3 on lysine 9 (H3K9me) are two repressive epigenetic modifications that are typically localized in distinct regions of the genome. For reasons unknown, however, they co-occur in some organisms and special tissue types. In this study, we show that maternal alleles marked by H3K27me3 in the Arabidopsis endosperm were targeted by the H3K27me3 demethylase REF6 and became activated during germination. In contrast, maternal alleles marked by H3K27me3, H3K9me2, and CHG methylation (CHGm) are likely to be protected from REF6 targeting and remained silenced. Our study unveils that combinations of different repressive epigenetic modifications time a key adaptive trait by modulating access of REF6.

Role of Polycomb in the control of transposable elements

It is generally considered that Polycomb Repressive Complex (PRC)2 deposits the histone mark H3K27me3 on silent protein-coding genes, while transposable elements are repressed by DNA and/or H3K9 methylation. Yet, there is increasing evidence that PRC2 also targets and even silences transposable elements in representatives of several distantly related eukaryotic lineages. In plants and animals, H3K27me3 is present on transposable elements in mutants and specific cell types devoid of DNA methylation. In this Opinion, we summarize the experimental evidence for this phenomenon across the eukaryotic kingdom, and discuss its functional and evolutionary significance. We hypothesize that an ancestral role of Polycomb group (PcG) proteins was to silence transposable elements.

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.

Disruption of a ∼23-24 nucleotide small RNA pathway elevates DNA damage responses in Tetrahymena thermophila

Endogenous RNA interference (RNAi) pathways regulate a wide range of cellular processes in diverse eukaryotes, yet in the ciliated eukaryote, Tetrahymena thermophila, the cellular purpose of RNAi pathways that generate ∼23–24 nucleotide (nt) small (s)RNAs has remained unknown. Here, we investigated the phenotypic and gene expression impacts on vegetatively growing cells when genes involved in ∼23–24 nt sRNA biogenesis are disrupted. We observed slower proliferation and increased expression of genes involved in DNA metabolism and chromosome organization and maintenance in sRNA biogenesis mutants RSP1Δ, RDN2Δ, and RDF2Δ. In addition, RSP1Δ and RDN2Δ cells frequently exhibited enlarged chromatin extrusion bodies, which are nonnuclear, DNA-containing structures that may be akin to mammalian micronuclei. Expression of homologous recombination factor Rad51 was specifically elevated in RSP1Δ and RDN2Δ strains, with Rad51 and double-stranded DNA break marker γ-H2A.X localized to discrete macronuclear foci. In addition, an increase in Rad51 and γ-H2A.X foci was also found in knockouts of TWI8, a macronucleus-localized PIWI protein. Together, our findings suggest that an evolutionarily conserved role for RNAi pathways in maintaining genome integrity may be extended even to the early branching eukaryotic lineage that gave rise to Tetrahymena thermophila.

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.

Morphology, morphogenesis and phylogeny of Anteglaucoma orientalis n. sp. (Ciliophora, Oligohymenophorea), a freshwater hymenostomatid from northeastern China

A new hymenostomatid ciliate, Anteglaucoma orientalis n. sp., isolated from a freshwater pond in Harbin, northeastern China, was investigated using live observation and silver staining methods. Anteglaucoma orientalis is characterized as follows: size in vivo about 50-60 × 30-35 μm; oval body shape; buccal area occupies about 25% of body length; 28-36 somatic kineties; membranelle 1 having six or seven basal body rows, membranelle 2 five to seven rows, and membranelle 3 three rows; single macronucleus with one micronucleus attached. Morphogenesis of the genus Anteglaucoma is revealed for the first time. The main events during binary fission are as follows: morphogenesis begins with proliferation of kinetosomes in the middle part of postoral kinety 1, and kinetosomes of this primordial field multiply and organize to finally form the paroral membrane and membranelles 1-3 of the opisthe; the parental apparatus in the proter does not take part in the stomatogenetic process. Phylogenetic analyses based on SSU rRNA gene sequences show that Anteglaucoma orientalis n. sp. clusters with the type species, A. harbinensis Pan et al., 2017, with full support.

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.

Characterization and Comparative Analyses of Mitochondrial Genomes in Single-Celled Eukaryotes to Shed Light on the Diversity and Evolution of Linear Molecular Architecture

Determination and comparisons of complete mitochondrial genomes (mitogenomes) are important to understand the origin and evolution of mitochondria. Mitogenomes of unicellular protists are particularly informative in this regard because they are gene-rich and display high structural diversity. Ciliates are a highly diverse assemblage of protists and their mitogenomes (linear structure with high A+T content in general) were amongst the first from protists to be characterized and have provided important insights into mitogenome evolution. Here, we report novel mitogenome sequences from three representatives (Strombidium sp., Strombidium cf. sulcatum, and Halteria grandinella) in two dominant ciliate lineages. Comparative and phylogenetic analyses of newly sequenced and previously published ciliate mitogenomes were performed and revealed a number of important insights. We found that the mitogenomes of these three species are linear molecules capped with telomeric repeats that differ greatly among known species. The genomes studied here are highly syntenic, but larger in size and more gene-rich than those of other groups. They also all share an AT-rich tandem repeat region which may serve as the replication origin and modulate initiation of bidirectional transcription. More generally we identified a split version of ccmf, a cytochrome c maturation-related gene that might be a derived character uniting taxa in the subclasses Hypotrichia and Euplotia. Finally, our mitogenome comparisons and phylogenetic analyses support to reclassify Halteria grandinella from the subclass Oligotrichia to the subclass Hypotrichia. These results add to the growing literature on the unique features of ciliate mitogenomes, shedding light on the diversity and evolution of their linear molecular architecture.

Large-scale phylogenomic analysis provides new insights into the phylogeny of the class Oligohymenophorea (Protista, Ciliophora) with establishment of a new subclass Urocentria nov. subcl

The class Oligohymenophorea is one of the most diverse assemblage of ciliated protists, which are particularly important in fundamental biological studies including understanding the evolutionary relationships among the lineages. Phylogenetic relationships within the class remain largely elusive, especially within the subclass Peniculia, which contains the long-standing problematic taxa Urocentrum and Paranassula. In the present study, we sequenced the genomes and/or transcriptomes of six non-culturable oligohymenophoreans using single-cell sequencing techniques. Phylogenomic analysis was performed based on expanded taxon sampling of 85 taxa, including 157 nuclear genes encoding 36,953 amino acids. The results indicate that: (1) urocentrids form an independent branch that is sister to the clade formed by Scuticociliatia and Hymenostomatia, which, together with the morphological data, supports the establishment of a new subclass, Urocentria n. subcl., within Oligohymenophorea; (2) phylogenomic analysis and ortholog comparison reveal a close relationship between Paranassula and peniculines, providing corroborative evidence for removing Paranassula from Nassulida and elevating it as an order, Paranassulida, within the subclass Peniculia; (3) based on the phylogenomic analyses and morphological data, we hypothesize that Peritrichia is the earliest diverging clade within Oligohymenophorea while Scuticociliatia and Hymenostomatia share the most common ancestor, followed successively by Urocentria and Peniculia. In addition, stop codon analyses indicate that oligohymenophoreans widely use UGA as the stop codon, while UAR are reassigned to glutamate (peritrichs) or glutamine (others), supporting the evolutionary hypothesis.

A Polycomb repressive complex is required for RNAi-mediated heterochromatin formation and dynamic distribution of nuclear bodies

Polycomb group (PcG) proteins are widely utilized for transcriptional repression in eukaryotes. Here, we characterize, in the protist Tetrahymena thermophila, the EZL1 (E(z)-like 1) complex, with components conserved in metazoan Polycomb Repressive Complexes 1 and 2 (PRC1 and PRC2). The EZL1 complex is required for histone H3 K27 and K9 methylation, heterochromatin formation, transposable element control, and programmed genome rearrangement. The EZL1 complex interacts with EMA1, a helicase required for RNA interference (RNAi). This interaction is implicated in co-transcriptional recruitment of the EZL1 complex. Binding of H3K27 and H3K9 methylation by PDD1-another PcG protein interacting with the EZL1 complex-reinforces its chromatin association. The EZL1 complex is an integral part of Polycomb bodies, which exhibit dynamic distribution in Tetrahymena development: Their dispersion is driven by chromatin association, while their coalescence by PDD1, likely via phase separation. Our results provide a molecular mechanism connecting RNAi and Polycomb repression, which coordinately regulate nuclear bodies and reorganize the genome.

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.

Programmed genome rearrangements in ciliates

Ciliates are a highly divergent group of unicellular eukaryotes with separate somatic and germline genomes found in distinct dimorphic nuclei. This characteristic feature is tightly linked to extremely laborious developmentally regulated genome rearrangements in the development of a new somatic genome/nuclei following sex. The transformation from germline to soma genome involves massive DNA elimination mediated by non-coding RNAs, chromosome fragmentation, as well as DNA amplification. In this review, we discuss the similarities and differences in the genome reorganization processes of the model ciliates Paramecium and Tetrahymena (class Oligohymenophorea), and the distantly related Euplotes, Stylonychia, and Oxytricha (class Spirotrichea).

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 domesticated transposase ALP2 mediates formation of a novel Polycomb protein complex by direct interaction with MSI1, a core subunit of Polycomb Repressive Complex 2 (PRC2)

A large fraction of plant genomes is composed of transposable elements (TE), which provide a potential source of novel genes through “domestication”–the process whereby the proteins encoded by TE diverge in sequence, lose their ability to catalyse transposition and instead acquire novel functions for their hosts. In Arabidopsis, ANTAGONIST OF LIKE HETEROCHROMATIN PROTEIN 1 (ALP1) arose by domestication of the nuclease component of Harbinger class TE and acquired a new function as a component of POLYCOMB REPRESSIVE COMPLEX 2 (PRC2), a histone H3K27me3 methyltransferase involved in regulation of host genes and in some cases TE. It was not clear how ALP1 associated with PRC2, nor what the functional consequence was. Here, we identify ALP2 genetically as a suppressor of Polycomb-group (PcG) mutant phenotypes and show that it arose from the second, DNA binding component of Harbinger transposases. Molecular analysis of PcG compromised backgrounds reveals that ALP genes oppose silencing and H3K27me3 deposition at key PcG target genes. Proteomic analysis reveals that ALP1 and ALP2 are components of a variant PRC2 complex that contains the four core components but lacks plant-specific accessory components such as the H3K27me3 reader LIKE HETEROCHROMATION PROTEIN 1 (LHP1). We show that the N-terminus of ALP2 interacts directly with ALP1, whereas the C-terminus of ALP2 interacts with MULTICOPY SUPPRESSOR OF IRA1 (MSI1), a core component of PRC2. Proteomic analysis reveals that in alp2 mutant backgrounds ALP1 protein no longer associates with PRC2, consistent with a role for ALP2 in recruitment of ALP1. We suggest that the propensity of Harbinger TE to insert in gene-rich regions of the genome, together with the modular two component nature of their transposases, has predisposed them for domestication and incorporation into chromatin modifying complexes.

A large part of the genomes of plants and animals consists of transposable elements (TE), which are usually considered as selfish or parasitic as they encode proteins (transposases) which promote TE proliferation but not functions useful for their hosts. As a result, hosts have evolved ways of reducing TE proliferation, usually by modifying the DNA or chromatin of TE so that their transposases are no longer produced. Once the TE are inactivated they can no longer proliferate and over time they accumulate mutations and can evolve new functions, often beneficial for their hosts. This process is known as domestication and is increasingly recognised as a potent source of evolutionary novelty. For example, the CRISPR/Cas system that has provided the basis for a revolution in genetic engineering (“genome editing”) has evolved via domestication of transposons in bacteria. We have identified the ALP proteins, two domesticated transposases which function as components of an enzyme complex (PRC2) involved in modifying chromatin and regulating host gene activity in plants. Here we show how ALPs contact PRC2 and direct formation of a novel complex that lacks several of the usual components. The ALPs and related proteins will provide valuable tools for manipulating plant chromatin.

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.

Comparative Studies on the Polymorphism and Copy Number Variation of mtSSU rDNA in Ciliates (Protista, Ciliophora): Implications for Phylogenetic, Environmental, and Ecological Research

While nuclear small subunit ribosomal DNA (nSSU rDNA) is the most commonly-used gene marker in studying phylogeny, ecology, abundance, and biodiversity of microbial eukaryotes, mitochondrial small subunit ribosomal DNA (mtSSU rDNA) provides an alternative. Recently, both copy number variation and sequence variation of nSSU rDNA have been demonstrated for diverse organisms, which can contribute to misinterpretation of microbiome data. Given this, we explore patterns for mtSSU rDNA among 13 selected ciliates (representing five classes), a major component of microbial eukaryotes, estimating copy number and sequence variation and comparing to that of nSSU rDNA. Our study reveals: (1) mtSSU rDNA copy number variation is substantially lower than that for nSSU rDNA; (2) mtSSU rDNA copy number ranges from 1.0 × 104 to 8.1 × 105; (3) a most common sequence of mtSSU rDNA is also found in each cell; (4) the sequence variation of mtSSU rDNA are mainly indels in poly A/T regions, and only half of species have sequence variation, which is fewer than that for nSSU rDNA; and (5) the polymorphisms between haplotypes of mtSSU rDNA would not influence the phylogenetic topology. Together, these data provide more insights into mtSSU rDNA as a powerful marker especially for microbial ecology studies.

Chromatin Organization in Early Land Plants Reveals an Ancestral Association between H3K27me3, Transposons, and Constitutive Heterochromatin

Genome packaging by nucleosomes is a hallmark of eukaryotes. Histones and the pathways that deposit, remove, and read histone modifications are deeply conserved. Yet, we lack information regarding chromatin landscapes in extant representatives of ancestors of the main groups of eukaryotes, and our knowledge of the evolution of chromatin-related processes is limited. We used the bryophyte Marchantia polymorpha, which diverged from vascular plants circa 400 mya, to obtain a whole chromosome genome assembly and explore the chromatin landscape and three-dimensional genome organization in an early diverging land plant lineage. Based on genomic profiles of ten chromatin marks, we conclude that the relationship between active marks and gene expression is conserved across land plants. In contrast, we observed distinctive features of transposons and other repetitive sequences in Marchantia compared with flowering plants. Silenced transposons and repeats did not accumulate around centromeres. Although a large fraction of constitutive heterochromatin was marked by H3K9 methylation as in flowering plants, a significant proportion of transposons were marked by H3K27me3, which is otherwise dedicated to the transcriptional repression of protein-coding genes in flowering plants. Chromatin compartmentalization analyses of Hi-C data revealed that repressed B compartments were densely decorated with H3K27me3 but not H3K9 or DNA methylation as reported in flowering plants. We conclude that, in early plants, H3K27me3 played an essential role in heterochromatin function, suggesting an ancestral role of this mark in transposon silencing.

Morphogenetic characters of the model ciliate Euplotes vannus (Ciliophora, Spirotrichea): Notes on cortical pattern formation during conjugational and postconjugational reorganization

Ciliated protists represent a morphologically and genetically distinct group of single-celled eukaryotes which can reproduce asexually and sexually. Morphogenesis occurs in both asexual and sexual modes of reproduction which is of interest for researchers investigating cell differentiation, regeneration, systematics and evolution. However, studies of morphogenesis have concentrated almost entirely on the asexual mode. Here we use protargol staining to investigate the morphogenetic processes during sexual reproduction in the model species Euplotes vannus (Müller). The major events include: (1) two rounds of morphogenesis occur during sexual reproduction, i.e., conjugational and postconjugational reorganization; (2) in both processes the oral primordium is generated de novo in a pouch beneath the cortex; (3) the frontoventral-transverse cirri anlagen are formed de novo and fragment in a 3:3:3:3:2 pattern; (4) the leftmost cirrus and the paroral membrane do not change during conjugational morphogenesis, but reorganize de novo during postconjugational morphogenesis; (5) marginal cirral anlagen are formed de novo in both morphogenetic processes; (6) two or three caudal cirri are formed at the ends of the rightmost two or three old dorsal kineties; (7) the dorsal kineties are retained entirely. These results can serve as reference to investigate the morphogenetic events in the different stages of sexual reproduction.

Conjugation in Euplotes raikovi (Protista, Ciliophora): New Insights into Nuclear Events and Macronuclear Development from Micronucleate and Amicronucleate Cells

Ciliates form a distinct group of single-celled eukaryotes that host two types of nuclei (micro and macronucleus) in the same cytoplasm and have a special sexual process known as conjugation, which involves mitosis, meiosis, fertilization, nuclear differentiation, and development. Due to their high species diversity, ciliates have evolved different patterns of nuclear events during conjugation. In the present study, we investigate these events in detail in the marine species Euplotes raikovi. Our results indicate that: (i) conjugation lasts for about 50 h, the longest stage being the development of the new macronucleus (ca. 36 h); (ii) there are three prezygotic micronuclear divisions (mitosis and meiosis I and II) and two postzygotic synkaryon divisions; and (iii) a fragment of the parental macronucleus fuses with the new developing macronucleus. In addition, we describe for the first time conjugation in amicronucleate E. raikovi cells. When two amicronucleate cells mate, they separate after about 4 h without evident nuclear changes; when one amicronucleate cell mates with a micronucleate cell, the micronucleus undergoes regular prezygotic divisions to form migratory and stationary pronuclei, but the two pronuclei fuse in the same cell. In the amicronucleate cell, the parental macronucleus breaks into fragments, which are then recovered to form a new functional macronucleus. These results add new information on the process of conjugation in both micronucleate and amicronucleate Euplotes cells.

Comparative Transcriptome Analyses during the Vegetative Cell Cycle in the Mono-Cellular Organism Pseudokeronopsis erythrina (Alveolata, Ciliophora)

Studies focusing on molecular mechanisms of cell cycles have been lagging in unicellular eukaryotes compared to other groups. Ciliates, a group of unicellular eukaryotes, have complex cell division cycles characterized by multiple events. During their vegetative cell cycle, ciliates undergo macronuclear amitosis, micronuclear mitosis, stomatogenesis and somatic cortex morphogenesis, and cytokinesis. Herein, we used the hypotrich ciliate Pseudokeronopsis erythrina, whose morphogenesis has been well studied, to examine molecular mechanisms of ciliate vegetative cell cycles. Single-cell transcriptomes of the growth (G) and cell division (D) stages were compared. The results showed that (i) More than 2051 significantly differentially expressed genes (DEGs) were detected, among which 1545 were up-regulated, while 256 were down-regulated at the D stage. Of these, 11 randomly picked DEGs were validated by reverse transcription quantitative polymerase chain reaction (RT-qPCR); (ii) Enriched DEGs during the D stage of the vegetative cell cycle of P. erythrina were involved in development, cortex modifications, and several organelle-related biological processes, showing correspondence of molecular evidence to morphogenetic changes for the first time; (iii) Several individual components of molecular mechanisms of ciliate vegetative division, the sexual cell cycle and cellular regeneration overlap; and (iv) The P. erythrina cell cycle and division have the same essential components as other eukaryotes, including cyclin-dependent kinases (CDKs), cyclins, and genes closely related to cell proliferation, indicating the conserved nature of this biological process. Further studies are needed focusing on detailed inventory and gene interactions that regulate specific ciliated cell-phase events.

Further analyses on the phylogeny of the subclass Scuticociliatia (Protozoa, Ciliophora) based on both nuclear and mitochondrial data

So far, the phylogenetic studies on ciliated protists have mainly based on single locus, the nuclear ribosomal DNA (rDNA). In order to avoid the limitations of single gene/genome trees and to add more data to systematic analyses, information from mitochondrial DNA sequence has been increasingly used in different lineages of ciliates. The systematic relationships in the subclass Scuticociliatia are extremely confused and largely unresolved based on nuclear genes. In the present study, we have characterized 72 new sequences, including 40 mitochondrial cytochrome oxidase c subunit I (COI) sequences, 29 mitochondrial small subunit ribosomal DNA (mtSSU-rDNA) sequences and three nuclear small subunit ribosomal DNA (nSSU-rDNA) sequences from 47 isolates of 44 morphospecies. Phylogenetic analyses based on single gene as well as concatenated data were performed and revealed: (1) compared to mtSSU-rDNA, COI gene reveals more consistent relationships with those of nSSU-rDNA; (2) the secondary structures of mtSSU-rRNA V4 region are predicted and compared in scuticociliates, which can contribute to discrimination of closely related species; (3) neither nuclear nor mitochondrial data support the monophyly of the order Loxocephalida, which may represent some divergent and intermediate lineages between the subclass Scuticociliatia and Hymenostomatia; (4) the assignments of thigmotrichids to the order Pleuronematida and the confused taxon Sulcigera comosa to the genus Histiobalantium are confirmed by mitochondrial genes; (5) both nuclear and mitochondrial data reveal that the species in the family Peniculistomatidae always group in the genus Pleuronema, suggesting that peniculistomatids are more likely evolved from Pleuronema-like ancestors; (6) mitochondrial genes support the monophyly of the order Philasterida, but the relationships among families of the order Philasterida remain controversial due to the discrepancies between their morphological and molecular data.

Diversification of Transposable Elements in Arthropods and Its Impact on Genome Evolution

Transposable elements (TEs) are ubiquitous in arthropods. However, analyses of large-scale and long-term coevolution between TEs and host genomes remain scarce in arthropods. Here, we choose 14 representative Arthropoda species from eight orders spanning more than 500 million years of evolution. By developing an unbiased TE annotation pipeline, we obtained 87 to 2266 TE reference sequences in a species, which is a considerable improvement compared to the reference TEs previously annotated in Repbase. We find that TE loads are diversified among species and were previously underestimated. The highly species- and time-specific expansions and contractions, and intraspecific sequence diversification are the leading driver of long terminal repeat (LTR) dynamics in Lepidoptera. Terminal inverted repeats (TIRs) proliferated substantially in five species with large genomes. A phylogenetic comparison reveals that the loads of multiple TE subfamilies are positively correlated with genome sizes. We also identified a few horizontally transferred TE candidates across nine species. In addition, we set up the Arthropod Transposable Elements database (ArTEdb) to provide TE references and annotations. Collectively, our results provide high-quality TE references and uncover that TE loads and expansion histories vary greatly among arthropods, which implies that TEs are an important driving force shaping the evolution of genomes through gain and loss.