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.

Changes in Organic Acids, Phenolic Compounds, and Antioxidant Activities of Lemon Juice Fermented by Issatchenkia terricola.

High content of citric acid in lemon juice leads to poor sensory experience. The study aimed to investigate the dynamics changes in organic acids, phenolic compounds, and antioxidant activities of lemon juice fermented with Issatchenkia terricola WJL-G4. The sensory evaluation of fermented lemon juice was conducted as well. Issatchenkia terricola WJL-G4 exhibited a potent capability of reducing the contents of citric acid (from 51.46 ± 0.11 g/L to 8.09 ± 0.05 g/L within 60 h fermentation) and increasing total phenolic level, flavonoid contents, and antioxidant activities compared to those of unfermented lemon juice. A total of 20 bioactive substances, including 10 phenolic acids and 10 flavonoid compounds, were detected both in fermented and unfermented lemon juice. The lemon juice fermented for 48 h had better sensory characteristics. Our findings demonstrated that lemon juice fermented with Issatchenkia terricola exhibited reduced citric acid contents, increased levels of health-promoting phenolic compounds, and enhanced antioxidant activities.

Hidden genomic evolution in a morphospecies-The landscape of rapidly evolving genes in Tetrahymena.

A morphospecies is defined as a taxonomic species based wholly on morphology, but often morphospecies consist of clusters of cryptic species that can be identified genetically or molecularly. The nature of the evolutionary novelty that accompanies speciation in a morphospecies is an intriguing question. Morphospecies are particularly common among ciliates, a group of unicellular eukaryotes that separates 2 kinds of nuclei-the silenced germline nucleus (micronucleus [MIC]) and the actively expressed somatic nucleus (macronucleus [MAC])-within a common cytoplasm. Because of their very similar morphologies, members of the Tetrahymena genus are considered a morphospecies. We explored the hidden genomic evolution within this genus by performing a comprehensive comparative analysis of the somatic genomes of 10 species and the germline genomes of 2 species of Tetrahymena. These species show high genetic divergence; phylogenomic analysis suggests that the genus originated about 300 million years ago (Mya). Seven universal protein domains are preferentially included among the species-specific (i.e., the youngest) Tetrahymena genes. In particular, leucine-rich repeat (LRR) genes make the largest contribution to the high level of genome divergence of the 10 species. LRR genes can be sorted into 3 different age groups. Parallel evolutionary trajectories have independently occurred among LRR genes in the different Tetrahymena species. Thousands of young LRR genes contain tandem arrays of exactly 90-bp exons. The introns separating these exons show a unique, extreme phase 2 bias, suggesting a clonal origin and successive expansions of 90-bp-exon LRR genes. Identifying LRR gene age groups allowed us to document a Tetrahymena intron length cycle. The youngest 90-bp exon LRR genes in T. thermophila are concentrated in pericentromeric and subtelomeric regions of the 5 micronuclear chromosomes, suggesting that these regions act as genome innovation centers. Copies of a Tetrahymena Long interspersed element (LINE)-like retrotransposon are very frequently found physically adjacent to 90-bp exon/intron repeat units of the youngest LRR genes. We propose that Tetrahymena species have used a massive exon-shuffling mechanism, involving unequal crossing over possibly in concert with retrotransposition, to create the unique 90-bp exon array LRR genes.

Tetrahymena Comparative Genomics Database (TCGD): a community resource for Tetrahymena.

Ciliates are a large and diverse group of unicellular organisms characterized by having the following two distinct type of nuclei within a single cell: micronucleus (MIC) and macronucleus (MAC). Although the genomes of several ciliates in different groups have been sequenced, comparative genomics data for multiple species within a ciliate genus are not yet available. Here we collected the genome information and comparative genomics analysis results for 10 species in the Tetrahymena genus, including the previously sequenced model organism Tetrahymena thermophila and 9 newly sequenced species, and constructed a genus-level comparative analysis platform, the Tetrahymena Comparative Genomics Database (TCGD). Genome sequences, transcriptomic data, gene models, functional annotation, ortholog groups and synteny maps were built into this database and a user-friendly interface was developed for searching, visualizing and analyzing these data. In summary, the TCGD (http://ciliate.ihb.ac.cn) will be an important and useful resource for the ciliate research community.

Metal-Organic Framework Template Synthesis of NiCo2S4@C Encapsulated in Hollow Nitrogen-Doped Carbon Cubes with Enhanced Electrochemical Performance for Lithium Storage.

Owing to its richer redox reaction and remarkable electrical conductivity, bimetallic nickel cobalt sulfide (NiCo2S4) is considered as an advanced electrode material for energy-storage applications. Herein, nanosized NiCo2S4@C encapsulated in a hollow nitrogen-doped carbon cube (NiCo2S4@D-NC) has been fabricated using a core@shell Ni3[Co(CN)6]2@polydopamine (PDA) nanocube as the precursor. In this composite, the NiCo2S4 nanoparticles coated with conformal carbon layers are homogeneously embedded in a 3D high-conduction carbon shell from PDA. Both the inner and the outer carbon coatings are helpful in increasing the electrical conductivity of the electrode materials and prohibit the polysulfide intermediates from dissolving in the electrolyte. When researched as electrode materials for lithium storage, owing to the unique structure with double layers of nitrogen-doped carbon coating, the as-obtained NiCo2S4@D-NC electrode maintains an excellent specific capacity of 480 mAh g-1 at 100 mA g-1 after 100 cycles. Even after 500 cycles at 500 mA g-1, a reversible capacity of 427 mAh g-1 can be achieved, suggesting an excellent rate capability and an ultralong cycling life. This remarkable lithium storage property indicates its potential application for future lithium-ion batteries.

Epidermal growth factor-induced stimulation of proliferation and gene expression changes in the hypotrichous ciliate, Stylonychia lemnae.

Epidermal growth factor (EGF) induces proliferation of epidermal and epithelial tissues in mammals. However, the effect of EGF on the single-celled eukaryotes is not well characterized, especially in the protists. Ciliates, an important group of protists, are well characterized as both pollution indicators and model organisms for research. Stylonychia lemnae, is one of the most common free-living ciliates, widely distributed in ponds, rivers and marshes. Here, we report the role of EGF on cell proliferation stimulation in S. lemnae. The growth curve of S. lemnae was established, and the stimulation effect of EGF on the proliferation of S. lemnae was investigated. Based on the results, potential EGF receptors were identified in S. lemnae according to the conserved domains and gene expression. Differential gene expression revealed that EGF-induced genes in other organisms (e.g. antioxidant) also up-regulated in S. lemnae cells at propagation stages. In addition, our results showed that EGF could up-regulate the signal transduction-related processes in the decline stage of S. lemnae cells, indicating its potential function in apoptosis inhibition. In summary, this study reports findings of the first investigation of EGF effects in hypotrich ciliates, and establishes an additional system for the study of the molecular mechanisms of EGF actions in eukaryotic cell division and proliferation.

Tracing the structural evolution of eukaryotic ATP binding cassette transporter superfamily.

The ATP binding cassette (ABC) transporters superfamily is one of the largest classes of membrane proteins. The core of the ABC transporter protein is composed of transmembrane domains (TMDs) and nucleotide binding domains (NBD). Eukaryotes ABC transporters are classified into seven main families (ABCA to ABCG) based on sequence similarity and domain organizations. With different domain number and domain organizations, eukaryote ABC transporters show diverse structures: the single structure (NBD or TMD), the ABC2 structure (NBD-NBD), the half structure (TMD-NBD or NBD-TMD) and the full structure (TMD-NBD-TMD-NBD or NBD-TMD-NBD-TMD). However, studies on how various ABC transporter gene structures evolved is still absent. Therefore, in this study, we comprehensively investigated the structural evolution of eukaryotic ABC transporters. The seven eukaryote ABC transporter families (A to G) fell into three groups: A&G group, B,C&D group and E&F group. There were at least four times the number of NBD and TMD fusion events in the origin of the half structure transporter. Two fusion modes were found in the full and ABC2 structure origination. Based on these findings, we present a putative structural evolutionary path of eukaryote ABC transporters that will increase our understanding on their origin, divergence and function.

Genome of the facultative scuticociliatosis pathogen Pseudocohnilembus persalinus provides insight into its virulence through horizontal gene transfer.

Certain ciliates of the subclass Scuticociliatia (scuticociliates) are facultative parasites of fishes in which they cause a suite of diseases collectively termed scuticociliatosis. Hitherto, comparatively little was known about genetics and genomics of scuticociliates or the mechanism of scuticociliatosis. In this study, a laboratory culture of the facultatively pathogenic scuticociliate Pseudocohnilembus persalinus was established and its genome sequenced, giving the first genome of a marine ciliate. Genome-wide horizontal gene transfer (HGT) analysis showed P. persalinus has acquired many unique prokaryote-derived genes that potentially contribute to the virulence of this organism, including cell adhesion, hemolysis and heme utilization genes. These findings give new insights into our understanding of the pathology of scuticociliates.

ATP-binding cassette transporter enhances tolerance to DDT in Tetrahymena.

The reuse of dichlorodiphenyltrichloroethane (DDT) as an indoor residual spray was permitted by the World Health Organization in 2007, and approximately 14 countries still use DDT to control disease vectors. The extensive exposure of insects to DDT has resulted in the emergence of DDT resistance, especially in mosquitoes, and the mechanism for this resistance in mosquitoes has been widely reported. Spraying can also introduce DDT directly into surface water, and DDT can subsequently accumulate in microorganisms, but the mechanism for the resistance to DDT degradation in microorganisms is unclear. Using whole-genome microarray analysis, we detected an abcb15 gene that was up-regulated in a specific manner by DDT treatment in T. thermophile. The deduced ABCB15 peptide sequence had two transmembrane domains (TMDs) and two nucleotide-binding domains (NBDs) to form the structure TMD-NBD-TMD-NBD, and each NBD contained three conserved motifs: Walker-A, C-loop, and Walker-B, which indicated the T. thermophila abcb15 was a typical ABC transporter gene. The expression of ABCB15 fused with a C-terminal green fluorescent protein was found to be on the periphery of the cell, suggesting that ABCB15 was a membrane pump protein. In addition, cells with abcb15 partially knocked down (abcb15-KD) grew slower than wild-type cells in the presence of 256 mg L(-1) DDT, indicating the tolerance of abcb15-KD strain to DDT exposure was decreased. Thus, we suggest that in Tetrahymena, the membrane pump protein encoded by ABCT gene abcb15 can enhance the tolerance to DDT and protect cells from this exogenous toxin by efficiently pumping it to the extracellular space.

Cd-metallothioneins in three additional tetrahymena species: intragenic repeat patterns and induction by metal ions.

Ciliate metallothioneins (MTs) possess many unique features compared to the "classic" MTs in other organisms, but they have only been studied in a small number of species. In this study, we investigated cDNAs encoding subfamily 7a metallothioneins (CdMTs) in three Tetrahymena species (T. hegewischi, T. malaccensis, and T. mobilis). Four CdMT genes (ThegMT1, ThegMT2, TmalMT1, and TmobMT1) were cloned and characterized. They share high sequence similarity to previously identified subfamily 7a MT members. Tetrahymena CdMTs exhibit a remarkably regular intragenic repeat homology. The CdMT sequences were divided into two main types of modules, which had been previously described, and which we name "A" and "B". ThegMT2 was identified as the first MT isoform solely composed of module "B". A phylogenetic analysis of individual modules of every characterized Tetrahymena CdMT rigorously documents the conclusion that modules are important units of CdMT evolution, which have undergone frequent and rapid gain/loss and shuffling. The transcriptional activity of the four newly identified genes was measured under different heavy metal exposure (Cd, Cu, Zn, Pb) using real-time quantitative PCR. The results showed that these genes were differentially induced after short (1 h) or long (24 h) metal exposure. The evolutionary diversity of Tetrahymena CdMTs is further discussed with regard to their induction by metal ions.

A P450 gene associated with robust resistance to DDT in ciliated protozoan, Tetrahymena thermophila by efficient degradation.

Analysis of metabolic mechanisms of dichlorodiphenyltrichloroethane (DDT) accumulation and degradation in microorganisms, which could be used to reduce its hazard to higher organisms at the higher in the food chain, have not been investigated. Robust resistance to DDT (grows well in 256 mg/L DDT) and a surprising ability to degrade DDT (more than 70% DDT within 4h) were found in the ciliated protozoan Tetrahymena thermophila. A P450 gene (CYP5013C2) was found to respond specifically to DDT treatment. In the presence of 256 mg/L DDT, cells with overexpressing CYP5013C2 (p450-OE) grew faster and degraded DDT more efficiently than wild-type (WT) cells, while cells with CYP5013C2 partially knocked down (p450-KD) grew slower and exhibited reduced ability to degrade DDT compared to WT cells. Both dichlorodiphenyldichloroethylene (DDE) and dichlorodiphenyldichloroethane (DDD) were detected in cells after exposure to DDT, and the concentration of DDD in the p450-OE strain gradually decreased from 0.5 to 4h. Thus, we argue that this P450 gene (CYP5013C2), by efficiently degrading DDT to DDD, is associated with robust resistance to DDT in Tetrahymena, and that a strain overexpressing this gene has the potential to serve as bioreactor that degrades environmental DDT.

Identification and characterization of the arsenite methyltransferase from a protozoan, Tetrahymena pyriformis.

Arsenic (As) methylation in aquatic microbes plays a major role in the biogeochemistry of As. Protozoa, especially the free-living freshwater species, are important players in aquatic ecological health. In this study, an arsenite (As(III)) methyltransferase, TpyArsM, was identified and characterized in a free-living protozoan, Tetrahymena pyriformis. In order to confirm its function, TpyarsM gene was knocked-out in Tetrahymena and was also heterologously expressed in hypersensitive E. coli; these events resulted in expected decreases in As tolerance and methylation ability, respectively. In-vitro tests revealed that purified TpyArsM protein methylated inorganic As to mono- and di- methylarsenate, and also had the novel property of producing trimethylarsenite (TMA(III)) and dimethylarsine (Me2AsH) gases. This new methyltransferase gene, identified in a species near the base of the food web, has enriched our knowledge of As methyltransferases and has great potential for bioremediation of As-contaminated environments.

Impaired replication elongation in Tetrahymena mutants deficient in histone H3 Lys 27 monomethylation.

Replication of nuclear DNA occurs in the context of chromatin and is influenced by histone modifications. In the ciliate Tetrahymena thermophila, we identified TXR1, encoding a histone methyltransferase. TXR1 deletion resulted in severe DNA replication stress, manifested by the accumulation of ssDNA, production of aberrant replication intermediates, and activation of robust DNA damage responses. Paired-end Illumina sequencing of ssDNA revealed intergenic regions, including replication origins, as hot spots for replication stress in ΔTXR1 cells. ΔTXR1 cells showed a deficiency in histone H3 Lys 27 monomethylation (H3K27me1), while ΔEZL2 cells, deleting a Drosophila E(z) homolog, were deficient in H3K27 di- and trimethylation, with no detectable replication stress. A point mutation in histone H3 at Lys 27 (H3 K27Q) mirrored the phenotype of ΔTXR1, corroborating H3K27me1 as a key player in DNA replication. Additionally, we demonstrated interactions between TXR1 and proliferating cell nuclear antigen (PCNA). These findings support a conserved pathway through which H3K27me1 facilitates replication elongation.

Selecting one of several mating types through gene segment joining and deletion in Tetrahymena thermophila.

The unicellular eukaryote Tetrahymena thermophila has seven mating types. Cells can mate only when they recognize cells of a different mating type as non-self. As a ciliate, Tetrahymena separates its germline and soma into two nuclei. During growth the somatic nucleus is responsible for all gene transcription while the germline nucleus remains silent. During mating, a new somatic nucleus is differentiated from a germline nucleus and mating type is decided by a stochastic process. We report here that the somatic mating type locus contains a pair of genes arranged head-to-head. Each gene encodes a mating type-specific segment and a transmembrane domain that is shared by all mating types. Somatic gene knockouts showed both genes are required for efficient non-self recognition and successful mating, as assessed by pair formation and progeny production. The germline mating type locus consists of a tandem array of incomplete gene pairs representing each potential mating type. During mating, a complete new gene pair is assembled at the somatic mating type locus; the incomplete genes of one gene pair are completed by joining to gene segments at each end of germline array. All other germline gene pairs are deleted in the process. These programmed DNA rearrangements make this a fascinating system of mating type determination.

Tetrahymena functional genomics database (TetraFGD): an integrated resource for Tetrahymena functional genomics.

The ciliated protozoan Tetrahymena thermophila is a useful unicellular model organism for studies of eukaryotic cellular and molecular biology. Researches on T. thermophila have contributed to a series of remarkable basic biological principles. After the macronuclear genome was sequenced, substantial progress has been made in functional genomics research on T. thermophila, including genome-wide microarray analysis of the T. thermophila life cycle, a T. thermophila gene network analysis based on the microarray data and transcriptome analysis by deep RNA sequencing. To meet the growing demands for the Tetrahymena research community, we integrated these data to provide a public access database: Tetrahymena functional genomics database (TetraFGD). TetraFGD contains three major resources, including the RNA-Seq transcriptome, microarray and gene networks. The RNA-Seq data define gene structures and transcriptome, with special emphasis on exon-intron boundaries; the microarray data describe gene expression of 20 time points during three major stages of the T. thermophila life cycle; the gene network data identify potential gene-gene interactions of 15 049 genes. The TetraFGD provides user-friendly search functions that assist researchers in accessing gene models, transcripts, gene expression data and gene-gene relationships. In conclusion, the TetraFGD is an important functional genomic resource for researchers who focus on the Tetrahymena or other ciliates. Database URL: http://tfgd.ihb.ac.cn/

Transcriptome analysis of the model protozoan, Tetrahymena thermophila, using Deep RNA sequencing.

BACKGROUND: The ciliated protozoan Tetrahymena thermophila is a well-studied single-celled eukaryote model organism for cellular and molecular biology. However, the lack of extensive T. thermophila cDNA libraries or a large expressed sequence tag (EST) database limited the quality of the original genome annotation. METHODOLOGY/PRINCIPAL FINDINGS: This RNA-seq study describes the first deep sequencing analysis of the T. thermophila transcriptome during the three major stages of the life cycle: growth, starvation and conjugation. Uniquely mapped reads covered more than 96% of the 24,725 predicted gene models in the somatic genome. More than 1,000 new transcribed regions were identified. The great dynamic range of RNA-seq allowed detection of a nearly six order-of-magnitude range of measurable gene expression orchestrated by this cell. RNA-seq also allowed the first prediction of transcript untranslated regions (UTRs) and an updated (larger) size estimate of the T. thermophila transcriptome: 57 Mb, or about 55% of the somatic genome. Our study identified nearly 1,500 alternative splicing (AS) events distributed over 5.2% of T. thermophila genes. This percentage represents a two order-of-magnitude increase over previous EST-based estimates in Tetrahymena. Evidence of stage-specific regulation of alternative splicing was also obtained. Finally, our study allowed us to completely confirm about 26.8% of the genes originally predicted by the gene finder, to correct coding sequence boundaries and intron-exon junctions for about a third, and to reassign microarray probes and correct earlier microarray data. CONCLUSIONS/SIGNIFICANCE: RNA-seq data significantly improve the genome annotation and provide a fully comprehensive view of the global transcriptome of T. thermophila. To our knowledge, 5.2% of T. thermophila genes with AS is the highest percentage of genes showing AS reported in a unicellular eukaryote. Tetrahymena thus becomes an excellent unicellular model eukaryote in which to investigate mechanisms of alternative splicing.

[Alternative splicing of an ATP-binding cassette transporter ABCC10 in Tetrahymena thermophila].

ATP-binding cassette transporters (ABCT) could generate multiple transcripts through alternative splicing (AS) in mammalian. Some AS introduced PTC (premature terminal codon)-containing isoforms of ABCT couple with NMD (nonsense-mediated mRNA decay) to regulate relevant functions. However, there are no similar reports in lower organisms. This paper focuses on the unicellular protozoa Tetrahymena thermophila, based on the RNA-seq data of Tetrahymena thermophila, identified two alternative splicing variants of gene ABCC10 (SV1 and SV2). The SV2 contained an intron retention event at the fifth intron, and this 49 bp intron resulted in shift-frame and introduced PTC. Then, a knock-down Tetrahymena strain of gene UPF1 which is a key factor of NMD was constructed, and the expression levels of SV2 were performed using a real-time quantitative PCR. The results showed the expression levels of SV2 were up-regulated significantly in knock-down strain, indicating that SV2 was targeted by NMD, which is consistent to the mechanism which the AS introduced PTC-containing isoforms of ABCC proteins can be targeted by NMD in mammalian. Thus, we infer that this mechanism is highly evolutionary conserved in eukaryotes and was already functional in the last eukaryotic common ancestor.

Gene network landscape of the ciliate Tetrahymena thermophila.

BACKGROUND: Genome-wide expression data of gene microarrays can be used to infer gene networks. At a cellular level, a gene network provides a picture of the modules in which genes are densely connected, and of the hub genes, which are highly connected with other genes. A gene network is useful to identify the genes involved in the same pathway, in a protein complex or that are co-regulated. In this study, we used different methods to find gene networks in the ciliate Tetrahymena thermophila, and describe some important properties of this network, such as modules and hubs. METHODOLOGY/PRINCIPAL FINDINGS: Using 67 single channel microarrays, we constructed the Tetrahymena gene network (TGN) using three methods: the Pearson correlation coefficient (PCC), the Spearman correlation coefficient (SCC) and the context likelihood of relatedness (CLR) algorithm. The accuracy and coverage of the three networks were evaluated using four conserved protein complexes in yeast. The CLR network with a Z-score threshold 3.49 was determined to be the most robust. The TGN was partitioned, and 55 modules were found. In addition, analysis of the arbitrarily determined 1200 hubs showed that these hubs could be sorted into six groups according to their expression profiles. We also investigated human disease orthologs in Tetrahymena that are missing in yeast and provide evidence indicating that some of these are involved in the same process in Tetrahymena as in human. CONCLUSIONS/SIGNIFICANCE: This study constructed a Tetrahymena gene network, provided new insights to the properties of this biological network, and presents an important resource to study Tetrahymena genes at the pathway level.

[Expression analysis of 5 hsp70 genes in Tetrahymena thermophila].

Thirteen hsp70 genes with complete conserved domains were identified in Tetrahymena thermophila, and expression of five similar and non-intron hsp70 genes were analyzed. Under heat shock conditions of 37, 39 and 41oC, hsp70-2 mRNA had the highest relative expression level, suggesting it is closely related to heat shock. The basal level of constitutive T. thermophila hsp70-4 gene was high during 20 physiological/developmental stages of growth, starvation and conjugation, and it changed little upon exposure to heat shock: evidence that hsp70-4 is an hsc70 gene. The hsp70-4 cDNA is 2 208 bp long, and contains an open reading frame of 1 959 bp encoding 635 amino acids. Microarray data indicated that T. thermophila hsp70-3 gene probably participated in early starvation (0-12 h) stress and late conjugation (6-10 h) events, such as new macronuclear and micronuclear anlagen formation and old macronuclear elimination. However, hsp70-5 gene possibly participates in late starvation (12-15 h) stress and early conjugation (0-6 h) events such as micronuclear meiosis, micronuclear exchange and pronuclear fusion. Blast2GO indicated that they participated in dissimilar biological processes, suggesting hsp70-3 and hsp70-5 perform different functions.