Absolute quantification of chromosome copy numbers in the polyploid macronucleus of Tetrahymena thermophila at the single-cell level

Importin-9 and a TPR domain protein MpH drive periodic patterning of ciliary arrays in Tetrahymena

We explored how the number of structures is determined in an intracellular organelle series. In Tetrahymena, the oral apparatus contains three diagonal ciliary rows: M1, M2, and M3. During development, the M rows emerge by sequential segmentation of a group of basal bodies, starting with the longest and most anterior M1 and ending with the shortest and most posterior M3. The mpD-1 and mpH-1 alleles increase and decrease the number of M rows, respectively. We identify MpH as a TPR protein and MpD as an importin-9. Both proteins localize to the M rows and form concentration gradients. MpH is a row elongation factor whose loss shortens all M rows and often prevents the formation of M3. MpD limits row initiation after the emergence of M2. MpD could be a part of a negative feedback loop that limits row initiation when M1 assembly is properly advanced. We conclude that the forming oral apparatus has properties of a semi-autonomous intracellular developmental field.

Unveiling an ancient whole-genome duplication event in Stentor, the model unicellular eukaryotes

Whole-genome duplication (WGD) events are widespread across eukaryotes and have played a significant role in moulding the genetic architectures of diverse organisms. In the present study, the newly sequenced genome of a giant ciliated protist, Stentor roeselii, provides an opportunity for the analysis of the collinearity and retention of reciprocal best-hit genes between two Stentor species. As a main result, we have unveiled a previously undetected ancient WGD event shaping the genome of its congener, Stentor coeruleus, a model protist used in cytological and evolutionary studies. Genomes of two congeners, S. coeruleus and S. roeselii, are compared and analyzed, revealing that: (i) the former exhibits a much higher retention rate of colinear-gene pairs (28%) than does S. roeselii, and in S. coeruleus, 75% of genes that have a RBH hit in S. roeselii, have paralogs with high amino-acid identity, consistent with a WGD event in the lineage leading to S. coeruleus; (ii) the S. roeselii genome possesses extremely short intergenic regions, implying that the lengths of intergenic regions are under strong selection; (iii) the unique characteristics of introns may have been shaped in the common ancestor of heterotrichs; (iv) gene families that play a role in activities of multiple protein kinases and voltage-gated ion channels expanded rapidly in the ancestor of both taxa, possibly relating to the remarkable regenerative ability in Stentor. This study offers new insights into the evolutionary dynamics of ciliate genomes, with implications for understanding of the processes underlying the evolution of genomic complexity.

Cip1, a CDK regulator, determines heterothallic mating or homothallic selfing in a protist

Mating type (sex) plays a crucial role in regulating sexual reproduction in most extant eukaryotes. One of the functions of mating types is ensuring self-incompatibility to some extent, thereby promoting genetic diversity. However, heterothallic mating is not always the best mating strategy. For example, in low-density populations or specific environments, such as parasitic ones, species may need to increase the ratio of potential mating partners. Consequently, many species allow homothallic selfing (i.e., self-fertility or intraclonal mating). Throughout the extensive evolutionary history of species, changes in environmental conditions have influenced mating strategies back and forth. However, the mechanisms through which mating-type recognition regulates sexual reproduction and the dynamics of mating strategy throughout evolution remain poorly understood. In this study, we show that the Cip1 protein is responsible for coupling sexual reproduction initiation to mating-type recognition in the protozoal eukaryote Tetrahymena thermophila. Deletion of the Cip1 protein leads to the loss of the selfing-avoidance function of mating-type recognition, resulting in selfing without mating-type recognition. Further experiments revealed that Cip1 is a regulatory subunit of the Cdk19-Cyc9 complex, which controls the initiation of sexual reproduction. These results reveal a mechanism that regulates the choice between mating and selfing. This mechanism also contributes to the debate about the ancestral state of sexual reproduction.

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.

Global and local functions of the Fused kinase ortholog CdaH in intracellular patterning in Tetrahymena

Ciliates assemble numerous microtubular structures into complex cortical patterns. During ciliate division, the pattern is duplicated by intracellular segmentation that produces a tandem of daughter cells. In Tetrahymena thermophila, the induction and positioning of the division boundary involves two mutually antagonistic factors: posterior CdaA (cyclin E) and anterior CdaI (Hippo kinase). Here, we characterized the related cdaH-1 allele, which confers a pleiotropic patterning phenotype including an absence of the division boundary and an anterior-posterior mispositioning of the new oral apparatus. CdaH is a Fused or Stk36 kinase ortholog that localizes to multiple sites that correlate with the effects of its loss, including the division boundary and the new oral apparatus. CdaH acts downstream of CdaA to induce the division boundary and drives asymmetric cytokinesis at the tip of the posterior daughter. CdaH both maintains the anterior-posterior position of the new oral apparatus and interacts with CdaI to pattern ciliary rows within the oral apparatus. Thus, CdaH acts at multiple scales, from induction and positioning of structures on the cell-wide polarity axis to local organelle-level patterning.

A seven-sex species recognizes self and non-self mating-type via a novel protein complex

Although most species have two sexes, multisexual (or multi-mating type) species are also widespread. However, it is unclear how mating-type recognition is achieved at the molecular level in multisexual species. The unicellular ciliate Tetrahymena thermophila has seven mating types, which are determined by the MTA and MTB proteins. In this study, we found that both proteins are essential for cells to send or receive complete mating-type information, and transmission of the mating-type signal requires both proteins to be expressed in the same cell. We found that MTA and MTB form a mating-type recognition complex that localizes to the plasma membrane, but not to the cilia. Stimulation experiments showed that the mating-type-specific regions of MTA and MTB mediate both self- and non-self-recognition, indicating that T. thermophila uses a dual approach to achieve mating-type recognition. Our results suggest that MTA and MTB form an elaborate multifunctional protein complex that can identify cells of both self and non-self mating types in order to inhibit or activate mating, respectively.

Genome-wide identification and in-silico analysis of papain-family cysteine protease encoding genes in Tetrahymena thermophila

Tetrahymena thermophila is a promising host for recombinant protein production, but its utilization in biotechnology is mostly limited due to the presence of intracellular and extracellular papain-family cysteine proteases (PFCPs). In this study, we employed bioinformatics approaches to investigate the T. thermophila PFCP genes and their encoded proteases (TtPFCPs), the most prominent protease family in the genome. Results from the multiple sequence alignment, protein modeling, and conserved motif analyses revealed that all TtPFCPs showed considerably high homology with mammalian cysteine cathepsins and contained conserved amino acid motifs. The total of 121 TtPFCP-encoding genes, 14 of which were classified as non-peptidase homologs, were found. Remaining 107 true TtPFCPs were divided into four distinct subgroups depending on their homology with mammalian lysosomal cathepsins: cathepsin L-like (TtCATLs), cathepsin B-like (TtCATBs), cathepsin C-like (TtCATCs), and cathepsin X-like (TtCATXs) PFCPs. The majority of true TtPFCPs (96 out of the total) were in TtCATL-like peptidase subgroup. Both phylogenetic and chromosomal localization analyses of TtPFCPs supported the hypothesis that TtPFCPs likely evolved through tandem gene duplication events and predominantly accumulated on micronuclear chromosome 5. Additionally, more than half of the identified TtPFCP genes are expressed in considerably low quantities compared to the rest of the TtPFCP genes, which are expressed at a higher level. However, their expression patterns fluctuate based on the stage of the life cycle. In conclusion, this study provides the first comprehensive in-silico analysis of TtPFCP genes and encoded proteases. The results would help designing an effective strategy for protease knockout mutant cell lines to discover biological function and to improve the recombinant protein production in T. thermophila.

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.

Recombinant production of hormonally active human insulin from pre-proinsulin by Tetrahymena thermophila

Alternative cell factories, such as the unicellular ciliate eukaryotic Tetrahymena thermophila, may be required for the production of protein therapeutics that are challenging to produce in conventional expression systems. T. thermophila (Tt) can secrete proteins with the post-translational modifications necessary for their function in humans. In this study, we tested if T. thermophila could process the human pre-proinsulin to produce hormonally active human insulin (hINS) with correct modifications. Flask and bioreactor culture of T. thermophila were used to produce the recombinant Tt-hINS either with or without an affinity tag from a codon-adapted pre-proinsulin sequence. Our results indicate that T. thermophila can produce a 6 kDa Tt-hINS monomer with the appropriate disulfide bonds after removal of the human insulin signal sequence or endogenous phospholipase A signal sequence, and the C-peptide of the human insulin. Additionally, Tt-hINS can form 12 kDa dimeric, 24 kDa tetrameric, and 36 kDa hexameric complexes. Tt-hINS-sfGFP fusion protein was localized to the vesicles within the cytoplasm and was secreted extracellularly. Assessing the affinity-purified Tt-hINS activity using the in vivo T. thermophila extracellular glucose drop assay, we observed that Tt-hINS induced a significant reduction (approximately 21 %) in extracellular glucose levels, indicative of its functional insulin activity. Our results demonstrate that T. thermophila is a promising candidate for the pharmaceutical and biotechnology industries as a host organism for the production of human protein drugs.

Comparative Genomic and Transcriptomic Profiling Revealed the Molecular Basis of Starch Promoting the Growth and Proliferation of Balantioides coli

Carbohydrates are the main source of nutrition for B. coli, supplying energy for cell growth and development. The research aimed at investigating the mechanism of starch on the growth and replication of B. coli. Single-cell separation was used to isolate single trophozoites of B. coli under a stereomicroscope, transcriptomic profiling was conducted based on the SMART-seq2 single-cell RNA-seq method. Comparative genomic analysis was performed on B. coli and eight other ciliates to obtain specific and expanded gene families of B. coli. GO and KEGG enrichment analysis were used to analyze the key genes of B. coli under the action of starch in the present study. The results of single-cell RNA-seq depicts starch affected the growth and replication of B. coli in two ways: (1) the cell cycle was positively promoted by the activation of the cAMP/PKA signaling pathway via glycolysis; (2) the cell autophagy was suppressed through the PI3K/AKT/mTOR pathway. Genes involved in endocytosis, carbohydrate utilization, and the cAMP/PKA signaling pathway were highly enriched in both specific and expanded gene families of B. coli. Starch can be ingested and hydrolyzed into glucose, in turn affecting various biological processes of B. coli. The molecular mechanism of the effect of starch on the growth and proliferation of B. coli by promoting cell cycle and inhibiting the autophagy of trophozoites has been elucidated in our study.

iGDP: An integrated genome decontamination pipeline for wild ciliated microeukaryotes

Ciliates are a large group of ubiquitous and highly diverse single-celled eukaryotes that play an essential role in the functioning of microbial food webs. However, their genomic diversity is far from clear due to the need to develop cultivation methods for most species, so most research is based on wild organisms that almost invariably contain contaminants. Here we establish an integrated Genome Decontamination Pipeline (iGDP) that combines homology search, telomere reads-assisted and clustering approaches to filter contaminated ciliate genome assemblies from wild specimens. We benchmarked the performance of iGDP using genomic data from a contaminated ciliate culture and the results showed that iGDP could recall 91.9% of the target sequences with 96.9% precision. We also used a synthetic dataset to offer guidelines for the application of iGDP in the removal of various groups of contaminants. Compared with several popular metagenome binning tools, iGDP could show better performance. To further validate the effectiveness of iGDP on real-world data, we applied it to decontaminate genome assemblies of three wild ciliate specimens and obtained their genomes with high quality comparable to that of previously well-studied model ciliate genomes. It is anticipated that the newly generated genomes and the established iGDP method will be valuable community resources for detailed studies on ciliate biodiversity, phylogeny, ecology and evolution. The pipeline (https://github.com/GWang2022/iGDP) can be implemented automatically to reduce manual filtering and classification and may be further developed to apply to other microeukaryotes.

Developmentally Programmed Switches in DNA Replication: Gene Amplification and Genome-Wide Endoreplication in Tetrahymena

Locus-specific gene amplification and genome-wide endoreplication generate the elevated copy number of ribosomal DNA (rDNA, 9000 C) and non-rDNA (90 C) chromosomes in the developing macronucleus of Tetrahymena thermophila. Subsequently, all macronuclear chromosomes replicate once per cell cycle during vegetative growth. Here, we describe an unanticipated, programmed switch in the regulation of replication initiation in the rDNA minichromosome. Early in development, the 21 kb rDNA minichromosome is preferentially amplified from 2 C to ~800 C from well-defined origins, concurrent with genome-wide endoreplication (2 C to 8-16 C) in starved mating Tetrahymena (endoreplication (ER) Phase 1). Upon refeeding, rDNA and non-rDNA chromosomes achieve their final copy number through resumption of just the endoreplication program (ER Phase 2). Unconventional rDNA replication intermediates are generated primarily during ER phase 2, consistent with delocalized replication initiation and possible formation of persistent RNA-DNA hybrids. Origin usage and replication fork elongation are affected in non-rDNA chromosomes as well. Despite the developmentally programmed 10-fold reduction in the ubiquitous eukaryotic initiator, the Origin Recognition Complex (ORC), active initiation sites are more closely spaced in ER phases 1 and 2 compared to vegetative growing cells. We propose that initiation site selection is relaxed in endoreplicating macronuclear chromosomes and may be less dependent on ORC.